Validation of molecular markers associated with boron tolerance, powdery mildew resistance and salinity tolerance in field peas.

Field pea (Pisum sativum L.) is an important grain legume consumed both as human food and animal feed. However, productivity in low rainfall regions can be significantly reduced by inferior soils containing high levels of boron and/or salinity. Furthermore, powdery mildew (PM) (Erysiphe pisi) disease also causes significant yield loss in warmer regions. Breeding for tolerance to these abiotic and biotic stresses are major aims for pea breeding programs and the application of molecular markers for these traits could greatly assist in developing improved germplasm at a faster rate. The current study reports the evaluation of a near diagnostic marker, PsMlo, associated with PM resistance and boron (B) tolerance as well as linked markers associated with salinity tolerance across a diverse set of pea germplasm. The PsMlo1 marker predicted the PM and B phenotypic responses with high levels of accuracy (>80%) across a wide range of field pea genotypes, hence offers the potential to be widely adapted in pea breeding programs. In contrast, linked markers for salinity tolerance were population specific; therefore, application of these markers would be suitable to relevant crosses within the program. Our results also suggest that there are possible new sources of salt tolerance present in field pea germplasm that could be further exploited.

Introduction

Field pea (<n class="Chemical">spanpan> class="Species">Pisum sativum Lpan>>.) is an important grain legume consumed as an inexpensive source of proteins, complex <spn>an class="Chemical">carbohydrates, fiber, B group vitamins, <span class="Chemical">folate and minerals such as calcium, iron and potassium. Peas contain up to 87% total digestible nutrients, which along with their low sodium and fat content (Tiwari and Singh, 2012) make them a valuable food source. In Australia, field peas are grown on over 280,000 ha per annum (ABARES, 2014) and over 70% of this production is within the lower rainfall (≤230 mm) cropping regions of southern Australia (Leonforte, 2013).

Soils in these low rainfall regions are often sodic, <n class="Chemical">spanpan> class="Chemical">salinepan>>, and alkaline, and can have toxic levels of <spn>an class="Chemical">boron (B) (Rengasamy, 2006; Nuttall et al., 2010; Javid et al., 2014) that can significantly limit field pea growth and yield. Amelioration of these subsoil constraints is considered impractical from an economic or logistical perspective (Nuttall et al., 2010). Consequently, identification of germplasm that has a combination of B and salinity tolerance is likely to result in improved growth and yield in the low rainfall areas with harsh soils.

<n class="Chemical">spanpan> class="Chemical">Boronpan>> concentrations in the range of only 10–54 mg kg-1 in the soil will inhibit plant growth (Reid et al., 2004), and soils in the southern Australian cropping region can range from 2 to 52 mg kg-1 B (Nuttall et al., 2003). The physiological mechanisms of B tolerance include reduced uptake into roots and reduced translocation into shoots, resulting in less accumulation of B in plant tissues (Reid and Fitzpatrick, 2009). Bagheri et al. (1996) identified two additive loci with incomplete dominance that con<spn>an class="Chemical">ferred tolerance to high levels of B in <span class="Species">peas. While tolerance controlled by single genes has been identified in the model legume Medicago truncatula (Bogacki et al., 2013) and lentil (Kaur et al., 2013). QTL mapping has also recently identified a single locus for B tolerance in field peas and tightly linked markers have been identified (Sudheesh et al., 2014).

Powdery mildew in <n class="Chemical">spanpan> class="Species">peaspan>> is caused by <spn>an class="Species">Erysiphe pisi DC, and is a serious disease that can cause yield losses of 25–50% (Warkentin et al., 1996; Tiwari et al., 1997; Fondevilla and Rubiales, 2012). In Australia, the disease occurs frequently in parts of South Australia, Victoria, and New South Wales mainly in wetter regions, although it also occasionally occurs in drier shorter season pea growing areas (Davidson et al., 2004). The obligate biotrophic fungus flourishes with dewy nights and warm days, with an optimal temperature of 20°C for conidial germination (Davidson et al., 2004). The disease is most destructive in late-sown field pea crops and in late maturing cultivars (Davidson et al., 2004). Genetic resistance is important as it is more cost-effective and environmentally friendly than fungicide applications. The genetics of PM resistance in field pea is relatively well understood with three major loci (er1, er2, and Er3) reported (Fondevilla et al., 2008; Srivastava et al., 2012; Pavan et al., 2013). The er1 locus is on chromosome Ps VI (Timmerman et al., 1994) and is tightly linked to B tolerance (Sudheesh et al., 2014). The er1 resistance is considered stable and more effective than er2 or Er3. The recessive er1 locus is due to a loss-of-function mutation of PsMlo1 (Mildew Resistance Locus O; Humphry et al., 2011).

Salinity is a persistent problem to crop productivity in low rainfall regions of the world (Rengasamy et al., 2003; Javid et al., 2012). Toxic levels of <n class="Chemical">spanpan> class="Chemical">sodiumpan>> in the subsoil, often associated with highly alkaline soils, are the main cause of yield loss associated with salinity (Rengasamy, 2002, 2006). Salinity negatively af<spn>an class="Chemical">fects nearly all developmental stages in plants (Flowers et al., 2010), and these effects are generally due to water deficit, ion toxicity and ion imbalance (Marschner, 1995). Munns and Tester (2008) reported that plant growth responds to salinity in two phases; a rapid osmotic phase, and a slower ionic phase. Whole plant tolerance to salinity is achieved through osmotic adjustment, ion exclusion (avoidance) and compartmentalization of sodium and chloride ions (Blumwald, 2000; Munns and Tester, 2008). In recent years, salinity tolerance in field pea has become increasingly important in Australia due to a geographical shift of crop production toward environments characterized by shorter seasons, reduced water availability and marginal soils with higher transient soil salinity (Nuttall et al., 2010; Leonforte et al., 2013b). The genetics of plant salinity tolerance is complex and governed by multiple genes with small effects (Arraouadi et al., 2012; Leonforte et al., 2013b), and is highly affected by environmental conditions. However, with molecular markers, the genetic analysis of quantitative inheritance of salinity tolerance is possible (Koyama et al., 2001), and Leonforte et al. (2013b) have identified markers for two salinity QTLs [Vp = 12% (QTL-1), 19% (QTL-2)] from the moderately tolerant field pea cultivar Parafield.

The application of molecular markers has the ability to enable breeders to select germplasm and breeding lines on the basis of a simple cost-ef<n class="Chemical">spanpan> class="Chemical">fepan>>ctive DNA assay without the need to undertake extensive phenotypic evaluation. These markers must be reliable and closely linked to reduce the probability of recombination. The procedure must be straight forward so that it can be done in a timely manner within a breeding program and it must be efficient to enable significant cost and time savin<spn>an class="Disease">gs (Collard et al., 2005; Collard and Mackill, 2008). These savin<span class="Disease">gs were estimated to be 40% in wheat (Kuchel et al., 2005), and 75% in common bean (Yu et al., 2000) and potato (Slater et al., 2013).

The main objective of this work was to assess a near diagnostic marker Ps<n class="Chemical">spanpan> class="Disease">Mlpan>>o1 (Sudheesh et al., 2014) for B tolerance and PM disease resistance, and tightly linked markers for two salinity tolerance QTLs (Leonforte et al., 2013b) in a diverse set of field pea genotypes for potential application in pea breeding programs.

Materials and Methods

Phenotyping

Plant Material

A diverse set of field pea germplasm consisting of 171 genotypes was used for phenotypic and genotypic evaluation of B <n class="Chemical">spanpan> class="Disease">toxicity toleranpan>cepan>> and PM resistance. An additional 22 genotypes were included in the salinity studies (Supplementary Table S1). Separate screening trials for B tolerance, salinity tolerance and PM resistance were conducted under controlled env<spn>an class="Chemical">ironment conditions at the Department of Economic Development, Jobs, Transport and Resources (DEDJTR), Horsham, VIC, Australia.

Boron Screening

Six seeds of each genotype were sown in 25 cm square non-draining pots filled with sandy loam soil pre-mixed with 10 mg kg-1 B as <n class="Chemical">spanpan> class="Chemical">boric acidpan>> (<spn>an class="Chemical">H3BO3). Four genotypes were sown per pot and the seedlin<span class="Disease">gs were thinned to three per genotype 1 week post-emergence. Plants were watered on every third day to 100% field capacity after weighing the pots. The screening experiment contained three replications within a randomized complete block design and was conducted in a polyhouse maintained at 24 ± 4°C with a 16/8 h (light/dark) photoperiod.

Seven weeks after emergence, seedlin<n class="Chemical">spanpan> class="Disease">gspan>> were assessed for B <spn>an class="Disease">toxicity based on the percentage of leaf area af<span class="Chemical">fected per plant. A 1–6 scale was used where 1 had no leaf chlorosis or necrosis, 2 had ≤20% tip chlorosis and necrosis, 3 had 20–30% tip chlorosis and necrosis, 4 had 30–40% chlorosis and necrosis, 5 had 40–50% chlorosis and necrosis, and 6 had >50% chlorosis and necrosis of the total leaf area. A mean score was derived for all plants across all replications and genotypes were considered tolerant it they had a score of one or less or sensitive if the score was greater than 1.0.

PM Screening

Three seeds were sown in 14 cm diameter pots filled with Bio Gro® potting mix in a glasshouse maintained at 20 ± 2°C, with a 16/8 h (light/dark) photoperiod. The potting mix consisted of 1,000 L of composted pine bark, 1 kg Floranid® n class="Chemical">N32, 1 kg 8–9 monpan>ths Osmocote®, 1 kg 3–4 monpan>ths Osmocote®, 225 g MicroMax® Complete, 225 g pan> class="Chemical">SP Quality® FeEDDHA Chelate (6% Fe), 30 kg agricultural lime, and 2 kg Saturaid®. Following the method of Davidson et al. (2004), 2 weeks post emergence PM infected plants were distributed among the seedlings so that air-borne inoculum could spread around the trial by the glasshouse cooling fans. Seedlings were watered from below to prevent inoculum being washed off the leaves. The trial was arranged in three replicates in a randomized complete block design.

Seedlin<n class="Chemical">spanpan> class="Disease">gspan>> were assessed for resistance 2 weeks after the introduction of PM in<spn>an class="Chemical">fected plants and scored using a 0–4 scale, where 0 had no mycelial growth or sporulation, 1 had <5% mycelial growth and sporulation, 2 had 5–25% mycelial growth and sporulation, 3 had 25–75% mycelial growth and sporulation, and 4 had >75% mycelial growth and sporulation on the plant. Plant scores were averaged for each genotype and considered resistant if they had a mean score of less than or equal to one and were susceptible if they scored greater than one.

Salinity Screening

Salinity screening was conducted as described by Leonforte et al. (2013b). Briefly, six seeds of each genotype were sown 2 cm deep in 13 cm diameter pots into a 1:1 coarse river sand and 5 mm <n class="Chemical">spanpan> class="Chemical">bluestonepan>> chip medium. Each pot was watered daily with rainwater until emergence. The seedlin<spn>an class="Disease">gs were thinned to three per genotype 4 days post emergence and watered with Manutec hydroponic nutrient solution containing N (7.6%), P (3.1%), and K (18.2%) supplemented with 50 g L-1 <span class="Chemical">calcium nitrate. After 10 days a salinity treatment mixed with the hydroponic nutrient solution was applied to all genotypes at an initial rate of 3 dS m-1 NaCl. Subsequently, every third day the concentration of NaCl was increased by 3 dS m-1 with each watering, up to 18 dS m-1, and maintained at this concentration until assessment. All solutions were applied at a rate of 300 ml per pot directly to the growing medium surface, avoiding contact with the leaves. The trial consisted of three replicates in a randomized complete block design. The salinity tolerance score and screening method were based on a visual growth response scale (1–10) as described by Leonforte et al. (2013a). The genotypes tested were then grouped into four response categories; tolerant (mean score ≤2.0), moderately tolerant (mean score >2 to ≤4.0), moderately sensitive (mean score >4.0 to ≤6.0) or sensitive (mean score >6.0).

Statistical Analysis

Analysis of variance (An class="Chemical">NOVA) was conducted on all phenotyping data using GEpan> class="Chemical">NSTAT 14th edition for windows.

Genotyping

DNA was extracted from pooled leaf tissue samples harvested from up to three planpan>ts per genpan>otype using the DNeasy® 96 Plant Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. The concentration of DNA was confirmed using a spectrophotometer (Thermo-Scientific, Wilmington, DE, USA) at the two wavelength ratios of A260/230 and A260/280 nm.

A near-diagnostic marker (Ps<n class="Chemical">spanpan> class="Disease">Mlpan>>o1) identified for PM that was also closely linked to B tolerance (Sudheesh et al., 2014) was assessed on a diverse set of pea germplasm. The sequence characterized amplified region (SCAR) primer pair; 5′-ATGGCTGAAGAGGGAGTT-3′ and 5′-GGTAGCAGCTTGATTTGTGGATA-3′ (Sudheesh et al., 2014) was synthesized and PCR amplifications were performed as described previously (Kaur et al., 2012). PCR products were combined with the ABI GeneScan LIZ500 size standard and analyzed using an ABI3730xl (Li<spn>an class="Chemical">fe Technologies Australia Pty Ltd, Scoresby, VIC, Australia) capillary electrophoresis platform according to the manufacturer’s instructions. Allele sizes were scored using GeneMapper® 5 software package (Li<span class="Chemical">fe Technologies Australia Pty Ltd). The PsMlo1 marker produced an allele size of 296 bp in B tolerant and PM resistant genotypes whereas the allele size of 294 bp was observed for sensitive/susceptible genotypes for both traits (Sudheesh et al., 2014).

SNP genotyping for the <pan> class="Chemical">span class="Chemical">salt tolerance QTLs as described in Leonforte et al. (2013b) were performed using KASPTM genotyping chemistry (LGC Genomics, Middlesex, UK) with some modifications. Primers were designed using 100 bp upstream and downstream flanking sequences around the target SNP variant position (Supplementary Table S2). The PCR reactions were carried out in 10 μl reaction volumes containing 50 ng genomic DNA, 5 μl KASPTM 2x Master Mix (LGC Genomics) and 0.14 μl KASPTM Primer Mix (LGC Genomics). The amplifications were performed in the SureCycler 8800 Thermal Cycler (Agilent Technologies, Inc., Santa Clara, CA, USA) using a touchdown PCR protocol following manufacturer’s instructions. End-point fluorescence values were determined by scanning the PCR plates in the FLUOstar Omega microplate reader (BMG LABTECH, Ortenberg, Germany) and the cluster plot analysis was completed in KlusterCallerTM software (LGC Genomics).

Results

Phenotyping for Boron Tolerance and PM Resistance

Plant ren class="Chemical">sponpan>ses within genpan>otypes were genpan>erally conpan>sistenpan>t, although occasionpan>ally segregating. Individual genpan>otypic repan> class="Chemical">sponses under B treatment were uniform across the three replications and variance among the replications was non-significant (P = 0.259). However, the severity of B toxicity on genotypes varied significantly (P < 0.001) from 1 (no symptoms) to 6 (>50% chlorotic and necrotic). Of the 171 genotypes tested for B tolerance, 58 (34%) were tolerant and 113 (66%) were sensitive to B. B tolerant genotypes showed no symptoms of toxicity throughout the 7 weeks of the experiment. The sensitive genotypes showed chlorosis and brown speckles on older leaves as typical symptoms of B toxicity (Figure ).

Field pea symptoms of .

Similar to B responses, planpan>t repan> class="Chemical">sponses to <span class="Disease">PM infection within genotypes were generally consistent, although occasionally segregating. Genotypes showed consistent response across the three replicates for PM disease (P = 0.762). The severity of <span class="Disease">PM infection varied significantly on genotypes (P < 0.001) from zero infection (no symptoms) to a highest mean score of 4 (>75% infection). Of the 171 genotypes tested for PM resistance, 61 genotypes (36%) were resistant and 110 (64%) were susceptible to PM disease with conidiospores present on both leaves and stems of the plants (Figure ).

The phenotypic scores between B tolerance and n class="Chemical">PM resistanpan>ce were highly correlated (r = 0.959, P < 0.001). Fifty-six genpan>otypes were toleranpan>t to B anpan>d resistanpan>t to pan> class="Chemical">PM, and 108 genotypes were found to be sensitive/susceptible to both treatments. Of the 171 genotypes tested, only two genotypes (02-164-2 and 09HP300-10HO2-2) were tolerant to B but susceptible to PM, whereas five genotypes (05H097-06HOS2003, 05H347-06HOS2005, 06H310P-8, 06H362P-1, and 07H178P001) were sensitive to B but resistant to PM. B and PM phenotyping showed a bimodal distribution with the majority at the extremes of the scoring range, although there were a number of individuals with intermediate scores (Figure ). The intermediate category contained a total of 24 genotypes for PM resistance and 33 for B tolerance of which 10 were common.

Distribution of plant symptoms for .

Marker and Phenotype Correlation

The results for the marker analysis were highly correlated to the phenotypic results. Of the 171 genotypes tested, marker predictions were accurate for both traits for 135 genotypes (79%). A further seven genotypes were predicted accurately for n class="Chemical">PM, but not for B, genpan>erating 83% correct predictionpan>s for pan> class="Chemical">PM. These seven genotypes also showed an uncoupling of the two traits (Table ).

List of genotypes showing recombination between PM anpan>d B phenotypes anpan>d how they correlate with the Ps<pan> class="Chemical">span class="Disease">Mlo marker.

A total of 26 genotypes were deemed false positives for both n class="Chemical">PM anpan>d B. These genpan>otypes conpan>tained the resistanpan>ce/toleranpan>t allele but were susceptible/senpan>sitive phenpan>otypically. Furthermore, anpan>other three genpan>otypes were false negatives for both traits, anpan>d conpan>tained the susceptible/senpan>sitive marker allele anpan>d were phenpan>otypically resistanpan>t/toleranpan>t (Table ). Of the pan> class="Chemical">PM false positives, two genotypes were segregating for the PM trait and 9 for B tolerance (Tables and ). The two genotypes segregating for PM resistance were also segregating for B tolerance. All segregating genotypes were identified as breeding lines. Of the 29 genotypes where the marker was not predictive, only 11 were publicly released cultivars. Six of the false positive genotypes were breeding lines with Parafield in the pedigree, which was also a false positive (Table ).

List of false positives and negatives genotypes of both traits with respect to the Ps<pan> class="Chemical">span class="Disease">Mlo marker.

Salinity Phenotyping

Plant ren class="Chemical">sponpan>ses within genpan>otypes were genpan>erally conpan>sistenpan>t, although occasionpan>ally segregating. Individual genpan>otypic repan> class="Chemical">sponses in response to the salinity treatment were uniform across the three replications and variance among the replications was non-significant (P = 0.669). However, there was a significant (P < 0.001) variation among genotypes for salinity tolerance. Salinity phenotypes also showed a distribution with bimodal characteristics where the majority scored 2 or 9, although a number of individuals had intermediate scores (Figure ). Based on the salinity symptom scores, genotypes such as Parafield, Yarrum, Collegian, and Dun showed tolerant responses with minor chlorosis on margins of lower leaves, whereas Kaspa, Maki and Bonzer were sensitive to the salinity treatment with severe chlorosis and necrosis on more than 50% of the plant.

Distribution of salinity symptoms scores from 1 to 10. The arrows show the salinity scores for the parents, Parafield (P) and Kaspa (K), from the mapping population from which the <pan> class="Chemical">span class="Chemical">salt QTLs were identified.

Salinity QTL Validation

The genotypes were grouped into four classes according to the alleles from the two loci of additive ef<spanpan> class="Chemical">fepan>>ct with a total Vp of 31%. Class 1 were the null genotypes where markers indicated that they had the sensitive alleles for both QTLs. This class had a total of 64 genotypes, of which 36 were sensitive or moderately sensitive, while 28 were tolerant or moderately tolerant to salinity. Class 2 had the tolerant allele from the Ps III QTL only and consisted of four genotypes that were phenotypically sensitive. Class 3 had the tolerant allele from Ps VII QTL, and consisted of 69 genotypes. Of these 69 genotypes, 30 were tolerant or moderately tolerant, 39 were sensitive or moderately sensitive. Class 4 had both tolerant QTL alleles and contained 56 genotypes. Of these, 28 were categorized as tolerant or moderately tolerant, while 28 were sensitive or moderately sensitive. Each class contained approximately 50% of genotypes as sensitive and moderately sensitive and 50% as tolerant or moderately tolerant, irrespective of the number of tolerance alleles derived from Parafield (Table ). This indicates that the tolerance in Class 1 (and possibly Class 3) must be from dif<span class="Chemical">ferent genetic sources.

n class="Chemical">Number of genotypes with salinity symptom scores in four genotypic classes based on presence/absence of the salinity toleranpan>ce QTL markers.

Further investigation of the genotypes within Class 4 showed that 21 had Parafield in their pedigree (Supplementary Table S3). n class="Chemical">Nineteen of the 21 genotypes were phenotypically toleranpan>t, whereas two genotypes were sensitive to salinity (Figure ). There were 35 genotypes of non-Parafield origin, anpan>d of these, only four (c. 11%) were toleranpan>t to salinity, five (c. 14%) were moderately toleranpan>t, four (c. 11%) were moderately sensitive, while the remaining 22 genotypes (c. 63%), were sensitive (Figure ). Outside of Class 4 there were only three genotypes with Parafield in their pedigree anpan>d they belonged to Class 3, containing toleranpan>t allele from Ps III QTL only.

Average salinity score of Class 4 genotypes, black squares (.

Discussion

The current study evaluates the application of the Ps<n class="Chemical">spanpan> class="Disease">Mlpan>>o1 marker and its predictive power for PM resistance and B tolerance as well as closely linked markers identified for salinity tolerance. The set of genotypes were identified to determine if the markers could be widely applied across diverse germplasm. The majority of the diverse germplasm set consisted of field pea breeding lines (82%) with the remainder cultivars. However, the selected set of genotypes served the objective of comparison between diagnostic and closely linked markers for breeding applications as the full range of genotypic and phenotypic classes were observed.

The Ps<n class="Chemical">spanpan> class="Disease">Mlpan>>o marker used in the current study is a SCAR marker that was developed from the genome sequencing of the full-length Ps<spn>an class="Disease">Mlo gene in two resistant Australian field pea cultivars, Yarrum and ps1771 based on a 2-bp insertion in intron 11 (Sudheesh et al., 2014). The resistance to PM in field pea is known to be caused by loss of function at Mlo (Humphry et al., 2011; Pavan et al., 2011; Sudheesh et al., 2014). However, types of mutations in Mlo varies in different pea germplasm sets including frameshift mutations, insertions and premature termination of translation due to single base variations (Tiwari et al., 1997; Humphry et al., 2011; Pavan et al., 2011; Sudheesh et al., 2014). Co-mapping of PM resistance and B tolerance on Ps chromosome III suggests that both traits are controlled by major genes whose physical locations may be close to each other (Sudheesh et al., 2014). The phenotypic data from this study also showed strong linkage (r = 0.959) between PM resistance and B tolerance to confirm these findings. The frequency distribution of PM resistance and B tolerance showed a bimodal pattern in the diverse set of pea germplasm, with few genotypes in the intermediate classes which confirms that both traits are under simple genetic control. Only a few genotypes (4%) of the entire pea set showed intermediate responses for resistance/tolerance to PM and B indicating that other sources of partial resistance/tolerance could be present in the germplasm.

In most scenarios, the Ps<n class="Chemical">spanpan> class="Disease">Mlpan>>o1 marker correctly predicted the PM and B phenotypic respn>onses (c. 84% for PM, c. 80% for B) across the pea germplasm set. However, a small proportion of genotypes were observed as false positives and false negatives for both traits. False positives are instances where the marker predicts the presence of the trait; however, this is not reflected in the phenotype. The most likely explanations for these inconsistencies would be a recombination event between the marker and the locus, or the presence of an alternative marker haplotype that is indistinguishable from the expected allele size. False negatives (when the marker predicts the absence of the trait but the phenotype is contradictory) are potentially more interesting as they may also represent a recombination event, or they may have a completely different source of resistance/tolerance available for exploitation by the breeding programs. The total number of false positives and negatives observed for B tolerance was higher compared to PM resistance which is expected as the marker was developed from the PM resistance gene, PsMlo, hence performed with better accuracy. A number of genotypes, especially breeding lines, were segregating for both traits where the marker should be used with more caution by genotyping and phenotyping single plants for marker confirmation. The high correlation of Ps<span class="Disease">Mlo marker to the phenotypes across the diverse set of pea germplasm, presumably indicates that locus PsMlo1 (er1) is widespread within the gene pool. Two additional loci conferring PM resistance in field pea have been reported; er2 locus on Ps III (Katoch et al., 2010) has effects largely confined to leaves (Marx, 1986; Tiwari et al., 1997) and Er3 is a dominant gene introgressed into P. pisi from P. fulvum (Fondevilla et al., 2011). However, the possibility of er2 and Er3 genes being present in the current field pea diverse germplasm set is minimal (Antonio Leonforte pers comm), although their presence or other unknown genetic factors would explain false negative phenotypes.

The number of false positives was higher than the false negatives for both traits, which indicates that these genotypes may have common ancestory. Pedigree relationships of the field pea genotypes were investigated, and 16 of the 26 n class="Chemical">PM false positives had two commonpan> anpan>cestors, Gottschalks Viktoria anpan>d breeder line 25/1-22-BI-1. These genpan>otypes were used in the early 1980s, but are no lonpan>ger available to test if onpan>e of these may be the source of the false positive haplotype. There were sevenpan> other false positives anpan>d six were introductionpan>s from the South Australianpan> provincial breeding program. Their pedigree informationpan> onpan>ly goes back to parenpan>ts or granpan>dparenpan>ts; however, four of them have commonpan> anpan>cestors of Alma anpan>d Wirrega. In the 1980s, the provincial breeding programs commonpan>ly swapped germplasm so it is possible that the false positive <pan> class="Chemical">span class="Species">donor was exchanged between breeding programs, or there could have been multiple recombination events. Finally, there was one false positive, 96-235*5, that had a unique pedigree that was solely derived from the cultivars Morgan and Bohatyr. Pedigree information indicated that there were a limited number of germplasm <sppan>>an class="Species">donors that has given rise to the 26 PM false positive genotypes within the diverse set of pea germplasm.

In comparison to diagnostic markers for qualitative traits such as PM resistanpan>ce anpan>d B toleranpan>ce, linked markers for quanpan>titative traits such as salinity toleranpan>ce are much more challenging to validate across diverse germplasm. The markers for salinity toleranpan>ce used in the current study were developed from a mapping population of Kapan> class="Chemical">spa × Parafield, in which Parafield, a moderately <span class="Chemical">salt tolerant cultivar, contributed two QTLs explaining 12 and 19% of the phenotypic variation (Leonforte et al., 2013b). A low phenotypic variation contributed by the salinity QTLs reflects that these QTLs have <span class="Species">minor effects. Similar observations have also been made by several other studies, where QTLs associated with various abiotic stresses including salinity have been shown to explain only a low percentage of phenotypic variation. This could either be due to environmental effects (Collins et al., 2008; Ashraf and Foolad, 2013) or epistatic interactions that can inhibit the expression of a favorable allele (Holland, 2001; Peleg et al., 2009).

The predictive power of the salinity tolerance markers was limited when tested across the diverse set of pea germplasm. This observation, however, is not surprising because these markers are associated with two QTLs of moderate ef<n class="Chemical">spanpan> class="Chemical">fepan>>ct in a bi-parental background and there could be greater potential for recombination. Furthermore, the ef<spn>an class="Chemical">fectiveness of these salinity markers was significantly improved for genotypes that have Parafield, the tolerance <span class="Species">donor, in their pedigree. Within Class 4, 19 of the 21 Parafield derived genotypes (90%) were salt tolerant and had tolerant alleles from both QTLs. These results suggest that the markers have strong effects and worked reasonably well within the germplasm derived from Parafield. Similar responses of QTL linked markers have previously been shown in tomato where salinity tolerance QTLs worked well in offspring of populations from which they were originally developed (Monforte et al., 1997; Foolad, 2004; Ashraf and Foolad, 2013). Further examination of the Kaspa derived genotypes indicated that 76% of the sensitive genotypes had both sensitive marker alleles, whilst only 7% of the tolerant genotypes had both tolerant alleles, indicating that there were many false negatives. Over the last two decades, the field pea breeding program in Australia has moved away from Parafield type (fully leaved plants) to Kaspa type (semi-leafless plants), as the latter cultivar produced significant improvements in quality, agronomy and yield superiority. However, the breeding program has had a long history of selection under saline conditions with Parafield used in 39 crosses from 1995–2003, along with other salt donors. Kaspa came into the crossing program in 1997 and was used in 250 crosses up until 2007. Therefore, opportunities for Kaspa derived genotypes to pick up salt tolerance from other sources may have generated the high number of genotypes with tolerance. Given the quantitative nature of this trait, it is also possible that the diverse set of pea germplasm may contain many salt tolerance loci other than those contributed by Parafield.

Conclusion

The Ps<n class="Chemical">spanpan> class="Disease">Mlpan>>o1 marker validation across diverse germplasm demonstrates that it can be used for selection for B tolerance and PM resistance in new germplasm sources dif<spn>an class="Chemical">ferent from the one in which the marker association was originally found. The Ps<span class="Disease">Mlo1 marker offers potential to be widely adapted to identify PM resistant and B tolerant germplasm, making this marker suitable for implementation within the pea breeding program. The fact that a single marker is associated with both traits is an added advantage. Our results demonstrate that markers linked to salinity tolerance QTLs can be implemented in germplasm derived from Parafield. The presence of other sources of tolerance to salinity suggest that there is potential to exploit other uncharacterized QTLs in the germplasm. Breeding for higher salinity tolerance therefore would require multiple QTLs for introgression into existing donor varieties in pea breeding programs. In view of the quantitative nature of salinity, the Parafield loci may not be sufficient for complete tolerance hence several mapping populations are currently being evaluated for new sources of tolerance to salinity. In addition, other genomic assisted approaches such as genome wide association studies (GWAS) and genomic selection (GS) may offer alternative solutions to understand polygenic traits with smaller genetic effects including salinity tolerance.

Author Contributions

NC, AS, anpan>d SK contributed toward conception anpan>d interpretation of work. AL anpan>d PK organpan>ized <pan> class="Chemical">span class="Species">peas germplasm for this work. SS and <sppan>>an class="Disease">ML conducted genotyping work and data analysis. PK and MJ conducted phenotyping work and data analysis. GR, AS, SK and MJ drafted and critically revised the manuscript for important intellectual content and all authors approve final version of manuscript for publication.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Table 1

List of genotypes showing recombination between PM and B phenotypes and how they correlate with the PsMlo marker.

Genotypes	PsMlo1 allele∗	PM symptom score	B symptom score
False positive (B)
05H097-06Hos2003	296	0 (R)	5 (S)†
05H347-06Hos2005	296	0 (R)	5 (S)†
06H310P-8	296	0 (R)	4 (S)†
06H362P-1	296	0 (R)	6 (S)
07H178P001	296	0 (R)	4 (S)†
False negative (B)
02-164-2	294	4 (S)	1 (T)
09Hp300-10Ho2-2	294	4 (S)	1 (T)

Table 2

List of false positives and negatives genotypes of both traits with respect to the PsMlo marker.

Genotypes	PsMlo1 allele∗	PM symptom score	B symptom score
False positive
Sydney	296	4 (S)	2 (S)†
00-254-28	296	4 (S)	6 (S)
01-226-2	296	3 (S)†	4 (S)†
05H245-06Hos2004	296	3 (S)	4 (S)†
06H461P-7	296	4 (S)	6 (S)
07H034P004	296	4 (S)	6 (S)
07H094P006	296	4 (S)	6 (S)
09Hp432-10Ho2-7	296	4 (S)	3 (S)†
95-072∗3	296	2 (S)†	3 (S)†
96-235∗5	296	4 (S)	6 (S)
97-031-6-3	296	4 (S)	6 (S)
97-724-5	296	4 (S)	6 (S)
Bonzer	296	4 (S)	5 (S)
Dunwa	296	4 (S)	6 (S)
Helena	296	4 (S)	6 (S)
Maro	296	4 (S)	6 (S)
Morgan	296	4 (S)	6 (S)
Odalette	296	4 (S)	6 (S)
P503-4-6-1	296	4 (S)	6 (S)
Parafield	296	3 (S)	6 (S)
Px-95-64-1-1	296	3 (S)	6 (S)
Px-96-79-8-1	296	4 (S)	6 (S)
Px-97-64	296	4 (S)	6 (S)
Santi	296	4 (S)	6 (S)
Soupa	296	4 (S)	6 (S)
Sturt	296	4 (S)	6 (S)
False negative
04H341P-05Ho2010-1	294	0 (R)	1 (T)
09Hp288-10Ho2-7	294	0 (R)	1 (T)
09Hp375-10Ho2-3	294	0 (R)	1 (T)

Table 3

Number of genotypes with salinity symptom scores in four genotypic classes based on presence/absence of the salinity tolerance QTL markers.

Response categories	Class 1 (Null)	Class 2 (QTL 1)	Class 3 (QTL 2)	Class 4 (both QTLs)	Total
Sensitive	32	4	35	24	95
Moderately sensitive	4	–	4	4	12
Moderately tolerant	16	–	16	9	41
Tolerant	12	–	14	19	45
Total	64	4	69	56	193