Genetic analysis of tolerance to the root lesion nematode Pratylenchus neglectus in the legume Medicago littoralis
BMC Plant Biology volume 14, Article number: 100 (2014)
The nematode Pratylenchus neglectus has a wide host range and is able to feed on the root systems of cereals, oilseeds, grain and pasture legumes. Under the Mediterranean low rainfall environments of Australia, annual Medicago pasture legumes are used in rotation with cereals to fix atmospheric nitrogen and improve soil parameters. Considerable efforts are being made in breeding programs to improve resistance and tolerance to Pratylenchus neglectus in the major crops wheat and barley, which makes it vital to develop appropriate selection tools in medics.
A strong source of tolerance to root damage by the root lesion nematode (RLN) Pratylenchus neglectus had previously been identified in line RH-1 (strand medic, M. littoralis). Using RH-1, we have developed a single seed descent (SSD) population of 138 lines by crossing it to the intolerant cultivar Herald. After inoculation, RLN-associated root damage clearly segregated in the population. Genetic analysis was performed by constructing a genetic map using simple sequence repeat (SSR) and gene-based SNP markers. A highly significant quantitative trait locus (QTL), QPnTolMl.1, was identified explaining 49% of the phenotypic variation in the SSD population. All SSRs and gene-based markers in the QTL region were derived from chromosome 1 of the sequenced genome of the closely related species M. truncatula. Gene-based markers were validated in advanced breeding lines derived from the RH-1 parent and also a second RLN tolerance source, RH-2 (M. truncatula ssp. tricycla). Comparative analysis to sequenced legume genomes showed that the physical QTL interval exists as a synteny block in Lotus japonicus, common bean, soybean and chickpea. Furthermore, using the sequenced genome information of M. truncatula, the QTL interval contains 55 genes out of which five are discussed as potential candidate genes responsible for the mapped tolerance.
The closely linked set of SNP-based PCR markers is directly applicable to select for two different sources of RLN tolerance in breeding programs. Moreover, genome sequence information has allowed proposing candidate genes for further functional analysis and nominates QPnTolMl.1 as a target locus for RLN tolerance in economically important grain legumes, e.g. chickpea.
Pasture legumes are used in crop rotations with cereals and oilseeds and bring the benefit of fixing atmospheric nitrogen, improving soil organic matter, soil structure and functioning as a break crop for diseases . Neutral to alkaline soils are widespread in the southern cropping regions of Australia and annual medics (Medicago spp.) are widely grown on these soils with Medicago littoralis and M. truncatula (barrel medic) as the predominant annual medics [2, 3]. These conditions appear favourable to the root lesion nematode (RLN) Pratylenchus neglectus that is widely distributed in South Australian cropping soils  and affects all major legumes and cereal crops . Pratylenchus spp. are migratory endoparasitic nematodes that feed and migrate within root cortical tissue causing necrosis and reduced lateral branching of roots upon infection . Damaging the root system, RLN has the potential to reduce water use efficiency, particularly in low rainfall environments  and can exacerbate infections by other soil pathogens. Due to its wide host range, resistance and tolerance to RLN has become a breeding priority in wheat and legume breeding as yield losses can be severe. In wheat, losses between 10-30% are attributed to P. neglectus, costing the wheat industry in South and Western Australia over $190 M each year . Although being moderately resistant to RLN [4, 5], annual medics are intolerant, which means that they are able to inhibit multiplication of nematode numbers but suffer considerable production losses (8-20%) in fields with high nematode numbers . In 2003–2004, the analysis of 389 soil samples from South Australia suggested that over 90% of the soil samples were infested with P. neglectus and 27% of the samples contained more than 20 nematodes per gram of soil . Experience based on field assays suggests that each nematode per gram soil translates into one percent of loss in yield in annual medics and severe damage to wheat, especially under low rainfall conditions where nematode-affected root systems are unable to use the sparsely available water [6, 9].
The control of RLN through use of nematicides is not a viable option as agrochemicals increase production costs and due to their non-specific side effects represent risks to human health and the wider environment. Chemicals such as dibromochloropropane (DBCP), a common soil fumigant for nematode control, was widely used until the mid 1980s when it was banned after being linked to causing sterility among male workers. Soil steaming is another method to control nematodes but is not practical in broad acre agriculture. Thus, the use of new cultivars that carry genetic resistance and tolerance are the most favoured option to reduce the impact on crop yield due to RLN.
A previous germplasm screen of 225 strand medic lines led to the identification of the P. neglectus tolerant strand medic, RH-1 . In order to understand the genetics of the strong nematode tolerance observed in RH-1, we have established a single seed descent (SSD) population of 138 lines segregating for the tolerance trait. Closely linked molecular markers were identified by a dual approach, (1) genetic map construction and (2) QTL analysis and fine mapping using gene-based SNPs. Here, we report on the results of both analyses and the alignment of the genetic information to the fully sequenced M. truncatula genome encompassing a physical interval of about fifty genes.
Phenotypic analysis of tolerance to Pratylenchus neglectus
Plants that were not inoculated with nematodes had negligible root damage (Score < 0.06) indicating that root damage symptoms observed in inoculated treatments was caused by the nematode. All 138 SSD lines were scored for root damage caused by P. neglectus with a broad range of damage observed using a rating scale of 1 (tolerant) to 10 (intolerant) as shown in Figure 1. Highly significant differences (P < 0.01) in root damage occurred between the SSD lines. The tolerant (RH-1) and intolerant (Herald) parents responded as expected.
Tolerant parent RH-1 (score 1.73 ± 0.04) along with three SSD lines had root disease scores less than 2 (Figure 2A). Most commonly, the SSD lines (n = 71) had root damage scores between 2 and 4. The average damage score of the SSD population was 4.24 ± 0.07. The most intolerant SSD line had a root damage score of 8.2 ± 0.09 and was statistically similar to the intolerant parent Herald with score 7.81 ± 0.15 (Figure 2A). The visual RDS was highly correlated to the shoot dry weight (R2 = 0.8862) (Figure 2B). The highly significant phenotypic variation between SSD lines, combined with the disparate performance of the parents made the population well suited for use in the subsequent genetic studies. The observed phenotypic distribution in the population follows a bimodal profile describing a quantitative trait controlled by a major and some minor QTL effects.
By definition, tolerance reactions enable a plant to reduce or prevent damage by nematodes while resistance hinders nematodes to multiply , which would be noticeable as reduced numbers of nematodes in healthy roots with a low RDS value. To prove that the range of observed root damage in the parents and population lines was due to different levels of tolerance rather than resistance, P. neglectus DNA was quantified in the roots of previously scored lines. The quantification of nematode DNA in roots of inoculated plants showed very similar values for both parents, Herald and RH-1 (about 900/plant) (Figure 3). Nematode numbers in the roots of tested SSD lines all exceeded 1000, regardless of RDS score (Figure 3), supporting that the different root damage levels were caused by a tolerance mechanism of M. littoralis to P. neglectus.
Genetic map construction and QTL identification
The here developed genetic map of RH-1 x Herald contained 130 markers with 112 M. truncatula SSRs and 18 InDel and CAPS markers. The InDel and CAPS markers were derived from polymorphisms identified when sequencing PCR products from RH-1 and Herald using primers that we had designed on published M. truncatula BAC sequences at the tolerance locus. A highly significant QTL with a LOD score of 19.9 was identified on a linkage group that contained 30 SSRs and gene-derived markers that, apart from SSR (AC145024), were all based on sequence information of M. truncatula chromosome 1 (Figure 4). The QTL was called QPnTolMl.1, for tolerance to Pratylenchus neglectus in M. littoralis and explained 49% of the phenotypic variation in the RH-1 x Herald SSD population and the tolerance allele was donated by the tolerant parent RH-1. Information on the genetic linkage and physical position of markers that are closely linked to the QTL is summarised in Table 1. The most closely linked marker was the gene-derived CAPS marker STkin1BsrI. The genetic distance between the most closely QTL-flanking gene-derived CAPS markers STkin1BsrI and SatNlaIII was 0.6 cM in the M. littoralis population RH-1 x Herald corresponding to a physical interval of ca. 36 kb. The order of all gene-derived markers on this linkage group reflected the gene order in the sequenced genome of M. truncatula (Mt4.0 at http://www.jcvi.org/cgi-bin/medicago/browse.cgi?page=assembly_stats). The physical distance between the genes Medtr1g071480 and Medtr1g072420 from which STkinBsrI and XylnHincII were derived, respectively, spanned 412 kb harbouring 55 genes according to the most recent release (Mt4.0, August 2013).
Marker validation in advanced breeding lines
The four QTL-flanking markers STkin1BsrI, SatNlaIII, XylnHincII and STkin2HinfI (Figure 4) were validated in 9 advanced breeding lines that resulted from phenotypic selection from F4 to F6 for P. neglectus tolerance and other agronomical important traits. The lines carry two different P.n. tolerance sources, from RH-1 (M. littoralis) or RH-2 (M. truncatula ssp. tricycla).
Parent RH-2 showed the same tolerance allele as RH-1 for the four tested markers as did all available tolerant breeding lines, 5 M. littoralis lines derived from RH-1 x Herald and 4 M. truncatula lines originating from RH-2. The matching genotypes between the tolerant M. littoralis and tolerant M. truncatula ssp. tricycla breeding lines suggest that the same tolerance locus to P. neglectus, QPnTolMl.1, is present in both related species.
Comparative analysis of QTL QPnTolMl.1 in Medicagoand grain legumes
The genomic region of QPnTolMl.1 in M. littoralis had all gene-based markers derived from chromosome 1 of M. truncatula. When this region was aligned to the genomes of other legumes using the Legume Information System (LIS, at http://www.comparative-legumes.org), continuous synteny blocks were identified in Phaseolus vulgaris (common bean) chromosomes 01 and 07, Lotus japonicus chromosome 5, Cicer arietinum (chickpea) chromosome 4 and Glycine max (soybean) chromosomes 2, 3, 10 and 19. Transcripts of the interval-flanking genes (Medtr1g071480 and Medtr1g072420) and the discussed five candidate genes within the physical interval of the tolerance QTL in Medicago were also present in the grain legumes chickpea  and soybean, and to a lesser degree in common bean (Additional file 1: Table S1) (http://www.comparative-legumes.org).
The wide distribution and host range of Pratylenchus neglectus that includes wheat, canola, chickpea and pasture legumes has made tolerance and resistance target traits in crop breeding programs [4, 8, 13]. To date, genetic analysis has lead to the identification of five resistance QTL in wheat [14, 15] and also in barley  but no resistance or tolerance QTL to this pathogen has been described in legumes. To our knowledge, this is the first report on a genetic locus linked to P. neglectus tolerance in a legume species. We refer to the observed host response as a tolerance reaction since plants with low root damage scores do not restrict or prevent nematode multiplication, which would be the case for resistance .
The distribution of the tolerance phenotypes (root damage scores) in the SSD population showed neither positive nor negative transgressive segregation with both parents as the extreme tolerant or intolerant lines. A single SSD line had a mean RDS of 8.2, which exceeded that of the intolerant parent Herald but the difference was not significant. The genetic analysis of the population revealed the very strong QTL, QPnTolMl.1, explaining 49% of the genetic variation. Thanks to the M. truncatula genome sequencing efforts (; Mt4.0 at http://www.jcvi.org/cgi-bin/medicago/browse.cgi?page=assembly_stats) and user-friendly presentation and access to the genome information, a link between genetic and physical map was straight forward.
At the tolerance QTL, the order of the M. truncatula SSRs, originally derived from BAC clones, as well as the gene-based markers developed from primers designed on M. truncatula gene sequences, matched the marker order in M. littoralis. The conserved synteny between these two species of the same genus is expected as it has previously been suggested that M. littoralis is a subspecies of M. truncatula though now accepted as separate species . M. littoralis and M. truncatula can be hybridised  and fertile hybrids have been generated to introgress desirable traits from one to the other species and cultivars have been released [21, 22]. A sequence comparison of the genomic tolerance region between M. truncatula and other legumes suggested that macrosynteny is conserved between medics, Lotus japonicus and the grain legumes chickpea, common bean and soybean. In soybean, comparative analysis revealed synteny blocks to four chromosomes due to the two genome duplication events 44 and 15 million years ago [23, 24].
Assuming that each gene in the QTL interval flanked by STkinBsrI and XylnHincII is also present in strand medic, the annotated gene functions (Mt4.0 at http://www.jcvi.org/cgi-bin/medicago/browse.cgi?page=assembly_stats) can be used to screen for candidates that might be responsible for the observed tolerance phenotype. To date, no resistance gene to P. neglectus has been cloned in any plant species but a few resistance genes to sedentary nematode species have been isolated. Examples of cloned nematode resistance genes are the Mi-1.2 gene from tomato conferring resistance to three species of root knot nematodes Meloidogyne incognita, M. arenaria and M. javanica; the Gpa2 gene in potato mediating resistance to the cyst nematode Globodera pallida and the Hero gene from tomato conferring resistance to two potato cyst nematodes Globodera rostochiensis and G. pallida. The structure of the mentioned nematode resistance genes is comparable to those described for viral, fungal or bacterial resistance genes, encoding nucleotide-binding sites (NBS) and leucine-rich-repeat (LRR) domains . Notable is that the reported resistance loci contain clusters of homologous genes that can confer resistance to diseases other than nematodes, e.g. the homologous gene Rx1 at the Gpa2 locus mediates resistance to potato virus X . Clusters of 3, 4 and 14 homologous genes within a 52 kb, 115 kb and a 118 kb region were observed for the Mi, Gpa2 and Hero locus, respectively [25–27]. In this study, the physical QTL interval with its predicted 55 genes showed no cluster of genes resembling the aforementioned resistance genes or other known resistance gene classes . Given that the observed phenotype is tolerance and not resistance, different types of genes can be expected. The absence of typical R genes might be due to the different mode of action between resistance and tolerance. Gene products that lead to tolerance might have a more direct activity on the invading nematode, e.g. by inhibiting nematode penetration or movement through the root. Tolerance genes do not have to mobilize a defence mechanism that affects reproduction of the nematode in contrast to resistance, as recently observed in the resistance response of wheat to Pratylenchus thornei. Tolerance is simply an ability to reduce damage to the plant likely by genes encoding for proteins that can activate or are directly involved in plant defence to pathogens. Whereas no typical R genes were present within the tolerance QTL, candidate genes and gene clusters that could play a role in the reduction of damage caused by invading P. neglectus were present in the interval. Based on reported characterizations of homologues, the relatively small number of five genes remained to be discussed as candidate tolerance genes. Two genes involved in defence-related signalling were the patatin-like phospholipase gene Medtr1g072190 and the Rho GTPase gene (Medtr1g072280). Patatin-like phospholipases have been shown to play an active role in establishing successful defence responses against bacterial attacks. For example, silencing the PLP1 gene in Capsicum enhanced susceptibility while overexpression of CaPLP1 in Arabidopsis improved the plant’s resistance to bacteria . Rho GTPases have been linked to plant defence responses, susceptibility and cytoskeleton organization. Examples of pathosystems where a role for Rho GTPases has been discussed are the rice and the blast fungus Magnaporthe grisea and barley against the powdery mildew fungus Blumeria graminis f.sp. hordei. The three other candidate genes for the improvement of nematode tolerance are the two QTL flanking genes sat-1 and xyloglucanase-specific endoglucanase inhibitor, as well as lipid transfer protein genes. The sulfate/bicarbonate/oxalate exchanger and transporter sat-1 could be involved in delivering a higher level of the poisonous oxalate to the site of infection or maintain a higher level in the tolerant roots constitutively. An effective defence against chewing insects due to high levels of calcium oxalate crystals has been described in Medicago truncatula where mutant plants deficient in oxalate crystal formation showed clear preference and subsequent severe damage by caterpillar larvae . Both, xyloglucanase-specific endoglucanase inhibitor and lipid transfer protein genes appeared in small clusters at the tolerance QTL interval. LTPs are secreted proteins and have been proposed to support the transport of hydrophobic monomers like cutin into the cell wall to form protective layers against pathogens [34, 35]. A slightly larger gene cluster of three members was found for the gene encoding a xyloglucanase-specific endoglucanase inhibitor, referred to as XEGIP. This protein family and its inhibiting activity on fungal endoglucanases that target xyloglucans in plants had only fairly recently been described [36–38] and a large cluster has been detected in the potato genome .
The most closely QTL-linked marker in the SSD population was derived from the gene encoding a Serine/threonine protein kinase Nek2 but Nek2 and related kinases have only ever been reported in the context of cell cycle regulation  excluding it as a primary candidate for tolerance to nematodes.
The in silico screening for candidate genes conferring tolerance to P. neglectus using an advanced genome sequence draft with excellent embedded information has provided a basis for future functional analysis. Candidate re-sequencing in M. littoralis, gene expression analysis, transgenesis and a detailed analysis of the infection mechanism, similar to the one described for resistance to Pratylenchus thornei in wheat  are to follow.
At this stage, practical outputs of this research are four gene-based molecular markers that are available to track the two tolerance sources, RH-1 and RH-2, in Medicago breeding programs. The two closely linked markers flanking the P. neglectus tolerance QTL can be added to the list of breeder markers previously reported for sulfonylurea tolerance  and boron tolerance  and help to ensure that the rotation crop Medicago is a competitive option in environmentally challenged low rainfall zones. Moreover, with P. neglectus causing severe losses in grain legumes like chickpea  the identified QTL in M. littoralis represents a candidate locus for current genetics research in chickpea.
In the present study, a genetic mapping approach led to the identification of a major QTL for tolerance to the root lesion nematode Pratylenchus neglectus in Medicago littoralis. SSR and gene-based CAPS markers at the QTL interval mapped physically to chromosome 1 in M. truncatula and comparative genome analysis suggested that the region also exists as synteny blocks in other sequenced legume genomes. Thus, the identified QTL is a candidate locus for RLN tolerance in pasture and grain legume species. Supportive of this idea is that the linked gene-based markers derived from M. littoralis can also be used to select RLN tolerance derived from M. truncatula ssp. tricycla. Moreover, five genes that are predicted in the corresponding M. truncatula genome interval are likely candidates involved in reducing the root damage caused by nematodes.
A cross between the two strand medic (M. littoralis) genotypes RH-1 (tolerant) and Herald (intolerant) was made in 2004 to develop a single seed descent (SSD) population with a final number of 138 lines that were used for genetic analysis. The population was specifically developed by SARDI’s pasture breeding group to investigate the observed strong tolerance phenotype of line RH-1 (see below) to Pratylenchus neglectus. The SSD lines were at stage F5:6 when we analysed tolerance phenotypically and extracted DNA for genetic analysis.
In addition, 9 advanced breeding lines with quantified P. neglectus tolerance levels were used for genotyping with four PCR markers spanning the tolerance QTL. Advanced breeding lines are defined as largely homozygous breeding lines (F6) that combine different desirable traits and are potential future cultivars. Here, 9 advanced breeding lines were derived either from additional RLN tolerant lines of the mapping population, RH-1 x Herald, or from crosses between RH-2 and different M. truncatula cultivars. Herald was bred to generate a strand medic cultivar with resistance to the pests spotted alfalfa aphid and bluegreen aphid . RH-1 is a wild accession collected from Cyprus in 1987 and classified as strand medic (M. littoralis) whilst the tolerant line RH-2, also a wild accession that was collected from Cyprus in 1983 has been classified as a barrel medic subspecies, M. truncatula ssp. tricycla. The advanced breeding lines were derived from crosses between the RH-2 and P. neglectus intolerant progeny from crosses between the barrel medics Caliph, Jemalong and the strand medic Angel. The progeny were repeatedly phenotyped, tolerant lines were kept and the best five lines based on repeated agronomical performance in growth room and field were included in the linkage analysis between markers and P. neglectus tolerance.
Plant growth conditions
Seeds from lines for nematode testing were surface sterilised in 80% ethanol for 15 s, followed by 3% hypochlorite solution for 3 min and then thoroughly rinsed with sterile water. Seeds were stored at 4°C on moist filter paper in petri plates for 2 days and then placed in a 25°C incubator for 16 hrs to germinate.
Germinated seeds were planted into a mixture of approximately 200 g of coarse sand and vermiculite (50:50 vol:vol) contained in square plastic pots (50 × 50 mm diameter × 120 mm depth) and moistened with 40 mL of nutrient solution  containing a small amount of N (6 ug N/mL as CaNO3). Plants were grown for 44 days in a controlled environment room maintained at 20/15°C day/night temperature and 14/10 hr light/dark regime. Plants were watered with reverse osmosis water as required and additional N supplied (8 mL per pot of 19 mM CaNO3.4H20) at 18, 25, 32 and 38 days after sowing.
Pratylenchus neglectus was obtained from the SARDI population #99 (originally sourced from a soil collected from Cambrai, South Australia in 1991) maintained on carrot callus as previously described for P. vulnus. Nematodes were collected by placing carrot callus in funnels in a misting chamber under intermittent aqueous mist of 10 s every 10 min for 96 h  at room temperature (22°C). Nematodes extracted were counted and diluted with water to the inoculum concentration of 3000 nematodes per 1 mL. Nematodes were applied to each seedling in two inoculations, at 4 and 10 days after sowing. At each inoculation, 3000 nematodes were applied as two 500 μL aliquots, each dispensed into a 25-mm-deep hole in the soil mix on either side of the seedling, formed by gently pushing a 1000 μL Gilson disposable pipette tip into the soil. Nil nematode treatments received two 500 μL aliquots of water only. After inoculation, holes on either side of each seedling were filled using a small amount of the surface soil in each pot.
Phenotypic assessment of P. neglectustolerance
The tolerance response of the 138 SSD lines to P. neglectus was assessed in a growth chamber experiment. Eight repeats (each pot containing one seedling) of the 138 SSD lines and 32 repeats of each parental line (RH-1 and Herald) were inoculated with nematodes. An additional 8–16 repeats of each parental line were not inoculated with nematodes (negative controls). After the second nematode inoculation (10 days after sowing) all pots were arranged in a randomised block design. Repeats from each treatment were distributed and randomised within eight blocks in the growth chamber. Thirty-two days after the second inoculation, plants were removed from the pots and roots washed free of sand. Individual root systems were scored visually on a scale from 0 to 10 with low scores indicating low root damage (or high nematode tolerance) and vice versa (Figure 1). Differences between population individuals were assessed using one-way analysis of variance (ANOVA).
The nine validation lines resulted from the same phenotypic screen as used for the SSD population applied on progeny that was derived from crosses between the two tolerance sources, RH-1 and RH-2. Both sources were used to introgress tolerance into agronomically important varieties. Plants at F4, F5 and F6 were repeatedly assessed for tolerance and only tolerant plants (RDS ≤ 4) kept for variety development.
All roots of the inoculated mapping population lines were dried at 60°C. Nematode tolerance was tested by quantifying P. neglectus DNA in the roots of the parental and selected SSD lines that had been scored with different levels of tolerance. Total root DNA was extracted by the SARDI Root Testing Service [48–50] and the amount of P. neglectus DNA was quantified using a real-time TaqMan PCR system with primers specific to the internal transcribed spacer (ITS) region of P. neglectus (unpublished data). Nematode DNA was quantified and used with a standard curve that was established with known amounts of DNA per nematodes. The quantified levels of nematodes in the roots were expressed as number of nematodes per plant (Figure 3).
Genetic map construction
For map construction, microsatellite markers were obtained from primer sequence information at the Medicago HapMap project website (http://www.medicagohapmap.org/) and the list of reported markers [51, 52]. Out of the 240 screened SSRs, 110 showed a polymorphism between the parental lines RH-1 and Herald and were subsequently used for genotyping of the entire SSD population. These SSRs were used to construct a genetic linkage map using MapManager QTXb20  with the Kosambi mapping function  and a linkage criterion of P = 10-4. The marker order was finalized using RECORD . After QTL identification, 18 additional markers were developed on the basis of BAC sequences and later on using genes at the critical QTL region.
SSR markers were assayed in a 9.5 μL reaction mixture containing 0.2 mM dNTP, 1x PCR buffer, 1.5 mM MgCl2, 10 μM each of forward and reverse primer, and 2.5 U Taq DNA polymerase (Qiagen). A touchdown profile was used for PCR cycling comprising an initial denaturation step of 94°C for 2 min, followed by a total of 37 cycles of 94°C for 30 s, an annealing step for 30 s and 72°C for 30 s. Initial annealing temperature was 59°C and was reduced by 0.5°C for each of the next 8 cycles. The remaining 29 cycles had an annealing temperature of 55°C, the program ended with a 5 min extension at 72°C. The SSR-marker-amplified products were separated on 8% non-denaturing polyacrylamide gels (Sigma Aldrich, Australia) at constant 200 V for 180 min, and visualized with ethidium bromide staining.
Development and use of DNA markers
Additional markers were developed after QTL analysis based on the genome sequence of M. truncatula (Mt3.5). The sequence of co-located BACs and annotated genes therein were used to design primers for further SNP identification and development of CAPS (cleaved amplified polymorphism) markers. An additional 18 markers, 3 based on InDels and 15 being CAPS, were added to the genetic map. PCR reactions were run at 94°C for 2 min, followed by 94°C for 30 s, 59°C for 30–60 s and 72°C for 60 s. Amplicons were either directly separated on a 2% agarose gel (InDels) or first digested by a SNP specific restriction enzyme (as indicated in the marker name in Table 1) according to the manufacturer’s protocol. Marker names, primer sequences and amplicon sizes of the mapped loci close to QPnTolMl.1 (Figure 4) are provided in Table 1.
To compare the QTL region of QPnTolMl.1 among different legumes genomes, the positions of predicted genes in M. truncatula were aligned with corresponding regions in common bean (Phaseolus vulgaris), soybean (Glycine max), and chickpea (Cicer arietinum) (Additional file 1: Table S1). Databases used for comparisons were Medicago truncatula genome sequence release version 4.0 (jcvi.org/medicago/jbrowse), Legume Information System, LIS (medtr.comparative-legumes.org/gb2/gbrowse/Mt3.5.1) and Gramene release 39 (gramene.org), the latter two databases using M. truncatula genome sequence release version 3.5.
Root lesion nematode
Quantitative trait locus
Single Seed Descent
- SSR marker:
Simple sequence repeat marker.
Angus JF: Nitrogen supply and demand in Australian agriculture. Aust J Exp Agr. 2001, 41: 277-288. 10.1071/EA00141.
Cocks PS, Mathison MJ, Crawford EJ: From Wild Plants to Pasture Cultivars: Annual Medics and Subterranean Clover in Southern Australia. Advances in Legume Science. Edited by: Summerfield RJ, Bunting AH. 1980, London: Ministry of Agriculture and Fisheries, 569-596. 1
Nichols PGH, Revell CK, Humphries AW, Howie JH, Hall EJ, Sandral GA, Ghamkhar K, Harris CA: Temperate pasture legumes in Australia - their history, current use, and future prospects. Crop Pasture Sci. 2012, 63: 691-725. 10.1071/CP12194.
Ballard RA, Hutton RE, Taylor SP, McKay AC, Howie JH: Field resistance of annual pasture legumes to the root lesion nematode, Pratylenchus neglectus. Australas Plant Path. 2006, 35: 303-308. 10.1071/AP06019.
Taylor SP, Hollaway GJ, Hunt CH: Effect of field crops on population densities of Pratylenchus neglectus and P. thornei in southeastern Australia; Part 1: P. neglectus. J Nematol. 2000, 32: 591-599.
Vanstone VA, Rathjen AJ, Ware AH, Wheeler RD: Relationship between root lesion nematodes (Pratylenchus neglectus and P. thornei) and the performance of wheat varieties. Aust J Exp Agr. 1998, 38: 181-188. 10.1071/EA97109.
Richards RA: Genetic opportunities to improve cereal root systems for dryland agriculture. Plant Prod Sci. 2008, 11: 12-16. 10.1626/pps.11.12.
Vanstone VA, Hollaway GJ, Stirling GR: Managing nematode pests in the southern and western regions of the Australian cereal industry: continuing progress in a challenging environment. Australas Plant Path. 2008, 37: 220-234. 10.1071/AP08020.
Taylor SP, Vanstone VA, Ware AH, McKay AC, Szot D, Russ MH: Measuring yield loss in cereals caused by root lesion nematodes (Pratylenchus neglectus and P. thornei) with and without nematicide. Aust J Agr Res. 1999, 50: 617-622. 10.1071/A98103.
Hutton RE, Ballard RA: Identification of Strand Medic (Medicago Littoralis) With Improved Tolerance of Pratylenchus Neglectus. Proceedings of the 3rd Australasian Soilborne Disease Symposium: 8–11 February 2004; Adelaide. Edited by: Ophel-Keller K, Hall B. 2004, South Australian Research and Development Institute, 212-213.
Trudgill DL: Resistance to and tolerance of plant parasitic nematodes in plants. Annu Rev Phytopathol. 1991, 29: 167-192. 10.1146/annurev.py.29.090191.001123.
Varshney RK, Song C, Saxena RK, Azam S, Yu S, Sharpe AG, Cannon S, Baek J, Rosen BD, Tar'an B, Millan T, Zhang X, Ramsay LD, Iwata A, Wang Y, Nelson W, Farmer AD, Gaur PM, Soderlund C, Penmetsa RV, Xu C, Bharti AK, He W, Winter P, Zhao S, Hane JK, Carrasquilla-Garcia N, Condie JA, Upadhyaya HD, Luo MC, et al: Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement. Nat Biotechnol. 2013, 31: 240-248. 10.1038/nbt.2491.
Thompson JP, Owen KJ, Stirling GR, Bell MJ: Root-lesion nematodes (Pratylenchus thornei and P. neglectus): a review of recent progress in managing a significant pest of grain crops in northern Australia. Australas Plant Path. 2008, 37: 235-242. 10.1071/AP08021.
Williams KJ, Taylor SP, Bogacki P, Pallotta M, Bariana HS, Wallwork H: Mapping of the root lesion nematode (Pratylenchus neglectus) resistance gene Rlnn1 in wheat. Theor Appl Genet. 2002, 104: 874-879. 10.1007/s00122-001-0839-3.
Zwart RS, Thompson JP, Godwin ID: Identification of quantitative trait loci for resistance to two species of root-lesion nematode (Pratylenchus thornei and P. neglectus) in wheat. Aust J Agric Res. 2005, 56: 345-352. 10.1071/AR04223.
Sharma S, Sharma S, Kopisch-Obuch FJ, Keil T, Laubach E, Stein N, Graner A, Jung C: QTL analysis of root-lesion nematode resistance in barley: 1. Pratylenchus neglectus. Theor Appl Genet. 2011, 122: 1321-1330. 10.1007/s00122-011-1533-8.
Young ND, Debelle F, Oldroyd GED, Geurts R, Cannon SB, Udvardi MK, Benedito VA, Mayer KFX, Gouzy J, Schoof H, Van de Peer Y, Proost S, Cook DR, Meyers BC, Spannagl M, Cheung F, De Mita S, Krishnakumar V, Gundlach H, Shiguo Zhou S, Mudge J, Bharti AK, Murray JD, Naoumkina MA, Rosen B, Silverstein KAT, Tang H, Rombauts S, Zhao PX, Zhou P, et al: The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature. 2011, 480: 520-524.
Small E, Brookes B: A numerical taxonomic analysis of the Medicago littoralis Medicago truncatula complex. Can J Bot. 1990, 68: 1667-1674. 10.1139/b90-214.
Small E: Alfalfa and Relatives: Evolution and Classification of Medicago. Ottawa: NRC Research Press 2011.
Simon JP, Millington AJ: Relationship in annual spceies of Medicago. III. The complex M. littoralis Rhode - M. truncatula Gaertn. Aust J Bot. 1967, 15: 35-73. 10.1071/BT9670035.
Lake AWH: Medicago truncatula Gaertn. (barrel medic) cv. Mogul. Aust J Exp Agr. 1993, 33: 823-824. 10.1071/EA9930823.
Peck DM, Howie JH: Development of an early season barrel medic (Medicago truncatula Gaertn.) with tolerance to sulfonylurea herbicide residues. Crop Pasture Sci. 2012, 63: 866-874. 10.1071/CP12130.
Choi HK, Mun JH, Kim DJ, Zhu H, Baek JM, Mudge J, Roe B, Ellis N, Doyle J, Kiss GB, Young ND, Cook DR: Estimating genome conservation between crop and model legume species. Proc Natl Acad Sci U S A. 2004, 101: 15289-15294. 10.1073/pnas.0402251101.
Zhu H, Choi H-K, Cook DR, Shoemaker RC: Bridging model and crop legumes through comparative genomics. Plant Physiol. 2005, 137: 1189-1196. 10.1104/pp.104.058891.
Milligan SB, Bodeau J, Yaghoobi J, Kaloshian I, Zabel P, Williamson VM: The root knot nematode resistance gene Mi from tomato is a member of the leucine zipper, nucleotide binding, leucine-rich repeat family of plant genes. Plant Cell. 1998, 10: 1307-1319. 10.1105/tpc.10.8.1307.
Van der Vossen EAG, Van der Voort J, Kanyuka K, Bendahmane A, Sandbrink H, Baulcombe DC, Bakker J, Stiekema WJ, Klein-Lankhorst RM: Homologues of a single resistance-gene cluster in potato confer resistance to distinct pathogens: A virus and a nematode. Plant J. 2000, 23: 567-576. 10.1046/j.1365-313x.2000.00814.x.
Ernst K, Kumar A, Kriseleit D, Kloos DU, Phillips MS, Ganal MW: The broad-spectrum potato cyst nematode resistance gene (Hero) from tomato is the only member of a large gene family of NBS-LRR genes with an unusual amino acid repeat in the LRR region. Plant J. 2002, 31: 127-136. 10.1046/j.1365-313X.2002.01341.x.
Boyes DC, McDowell JM, Dangl JL: Plant pathology: many roads lead to resistance. Curr Biol. 1996, 6: 634-637. 10.1016/S0960-9822(09)00435-7.
Linsell KJ, Riley IT, Davies KA, Oldach KH: Characterization of resistance to Pratylenchus thornei (Nematoda) in wheat (Triticum aestivum); attraction, penetration, motility and reproduction. Phytopathology. 2013, doi:10.1094/PHYTO-12-12-0345-R
Kim DS, Jeun Y, Hwang BK: The pepper patatin-like phospholipase CaPLP1 functions in plant cell death and defense signaling. Plant Mol Biol. 2013, doi 10.1007/s11103-013-0137-x
Chen L, Shiotani K, Togashi T, Miki D, Aoyama M, Wong HL, Kawasaki T, Shimamoto K: Analysis of the Rac/Rop small GTPase family in rice: expression, subcellular localization and role in disease resistance. Plant Cell Physiol. 2010, 51: 585-595. 10.1093/pcp/pcq024.
Hoefle C, Huesmann C, Schultheiss H, Börnke F, Hensel G, Kumlehn J, Hückelhoven R: A barley ROP GTPase activating protein associates with microtubules and regulates entry of the barley powdery mildew fungus into leaf epidermal cells. Plant Cell. 2011, 23: 2422-2439. 10.1105/tpc.110.082131.
Korth KL, Doege SJ, Park S-H, Goggin FL, Wang Q, Gomez SK, Liu G, Jia L, Nakata PA: Medicago truncatula mutants demonstrate the role of plant calcium oxalate crystals as an effective defense against chewing insects. Plant Physiol. 2006, 141: 188-195. 10.1104/pp.106.076737.
Blein J-P, Coutos-Thévenot P, Marion D, Ponchet M: From elicitins to lipid-transfer proteins: a new insight in cell signalling involved in plant defence mechanisms. Trends Plant Sci. 2002, 7: 293-296. 10.1016/S1360-1385(02)02284-7.
Kader JC: Lipid-transfer proteins in plants. Annu Rev Plant Physiol Plant Pathol. 1996, 47: 627-654. 10.1146/annurev.arplant.47.1.627.
Qin Q, Bergmann C, Rose J, Saladie M, Kolli V, Albersheim P, Darvill A, York W: Characterization of a tomato protein that inhibits a xyloglucan-specific endoglucanase. Plant J. 2003, 34: 327-338. 10.1046/j.1365-313X.2003.01726.x.
Costanzo S, Ospina-Giraldo MD, Deahl KL, Baker CJ, Jones RW: Gene duplication event in family 12 glycosylhydrolase from Phytophthora spp. Fungal Genet Biol. 2006, 43: 707-714. 10.1016/j.fgb.2006.04.006.
Jones RW, Ospina-Giraldo M, Deahl K: Gene silencing indicates a role for potato endoglucanase inhibitor protein in germplasm resistance to late blight. Am J Potato Res. 2006, 83: 41-46. 10.1007/BF02869608.
Jones RW: Multiple copies of genes encoding XEGIPs are harbored in an 85-kB region of the potato genome. Plant Mol Biol Rep. 2012, 30: 1040-1046. 10.1007/s11105-012-0416-8.
Fry AM, O’Regan L, Sabir SR, Bayliss R: Cell cycle regulation by the NEK family of protein kinases. J Cell Sci. 2012, 125: 4423-4433. 10.1242/jcs.111195.
Oldach KH, Peck DM, Cheong J, Williams KJ, Nair RM: Identification of a chemically induced point mutation mediating herbicide tolerance in annual medics (Medicago spp.). Ann Bot. 2008, 101: 997-1005. 10.1093/aob/mcn028.
Bogacki P, Peck DM, Nair RM, Howie J, Oldach KH: Genetic analysis of tolerance to Boron toxicity in the legume Medicago truncatula. BMC Plant Biol. 2013, 13: 54-10.1186/1471-2229-13-54.
Thompson JP, Reen RA, Clewett TG, Sheedy JG, Kelly AM, Gogel BJ, Knights EJ: Hybridisation of Australian chickpea cultivars with wild Cicer spp. increases resistance to root-lesion nematodes (Pratylenchus thornei and P. neglectus). Australas Plant Pathol. 2011, 40: 601-611. 10.1007/s13313-011-0089-z.
Lake AWH, Howie JH, Drewry RE, Hill JR, Robinson SS, Schutz PR, Hammer A, Heinrich NB: Register of Australian Herbage Plant Cultivars B. Legumes9. Annual Medics (b) Medicago littoralis Rhode (strand medic) cv. Herald. Aust J Exp Agr. 1997, 37: 609-610.
McKnight T: Efficiency of isolates of Rhizobium in the cowpea group, with proposed additions to this group. Qld J Agric Sci. 1949, 6: 61-76.
Moody EH, Lownsbery BF, Ahmed JM: Culture of the root-lesion nematode Pratylenchus vulnus on carrot disks. J Nematol. 1973, 19: 125-134. 10.1163/187529273X00268.
Southey JF: Principles of Sampling for Nematodes. Laboratory Methods for Work With Plant and Soil Nematodes. Edited by: Southey JF. 1986, London: Ministry of Agriculture, Fisheries and Food, 1-4.
Haling RE, Simpson RJ, McKay AC, Hartley D, Lambers H, Ophel-Keller K, Wiebkin S, Herdina , Riley IT, Richardson AE: Direct measurement of roots in soil for single and mixed species using a quantitative DNA-based method. Plant Soil. 2011, 348: 123-137. 10.1007/s11104-011-0846-3.
Ophel-Keller K, McKay A, Hartley D, Herdina , Curran J: Development of a routine DNA-based testing service for soilborne diseases in Australia. Australas Plant Pathol. 2008, 37: 243-253. 10.1071/AP08029.
Riley IT, Wiebkin S, Hartley D, McKay AC: Quantification of roots and seeds in soil with real-time PCR. Plant Soil. 2010, 331: 151-163. 10.1007/s11104-009-0241-5.
Eujayl I, Sledge MK, Wang L, May GD, Chekhovskiy K, Zwonitzer JC, Mian MAR: Medicago truncatula EST-SSRs reveal cross-species genetic markers for Medicago spp. Theor Appl Genet. 2004, 108: 414-422. 10.1007/s00122-003-1450-6.
Mun J-H, Kim D-J, Choi H-K, Gish J, Debelle F, Mudge J, Denny R, Endre G, Saurat O, Dudez A-M, Kiss GB, Roe B, Young ND, Cook DR: Distribution of microsatellites in the genome of Medicago truncatula: a resource of genetic markers that integrate genetic and physical maps. Genetics. 2006, 172: 2541-2555.
Manly KF, Olson JM: Overview of QTL mapping software and introduction to map manager QT. Mamm Genome. 1999, 10: 327-334. 10.1007/s003359900997.
Kosambi DD: The estimation of map distances from recombination values. Ann Eugenic. 1944, 12: 172-175.
Van Os H, Stam P, Visser RGF, Van Eck HJ: RECORD: A novel method for ordering loci on a genetic linkage map. Theor Appl Genet. 2005, 112: 30-40. 10.1007/s00122-005-0097-x.
We acknowledge the contributions by R.E. Hutton and B. Morgan through earlier nematode tolerance assessments. This work was supported by the Rural Industries Research and Development Corporation of Australia (RIRDC project PRJ-005062 to KHO) and the South Australian Government through SARDI.
The authors declare that they have no competing interests.
KHO conceptualised the idea, carried out all genomics and bioinformatics part of marker development and drafted the manuscript. RB was involved in conceptualising the idea, carried out all the phenotypic analysis and helped progressing the SSD population. DMP derived the breeding lines and developed the SSD population from the initial cross by RMN who was also involved in the conceptualisation of the project. MS carried out all of the SSR and gene-based mapping supported by breeding line genotyping by JH with PB having carried out the polymorphism checks of the SSR markers. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Comparative analysis of physical position of QTL QPnTolMl.1 (flanking genes) and discussed candidates (bold font) in M. truncatula and grain legumes according to LIS. (DOC 78 KB)
About this article
Cite this article
Oldach, K.H., Peck, D.M., Nair, R.M. et al. Genetic analysis of tolerance to the root lesion nematode Pratylenchus neglectus in the legume Medicago littoralis. BMC Plant Biol 14, 100 (2014). https://doi.org/10.1186/1471-2229-14-100
- Root disease
- Gene-based markers
- Genetic map
- Comparative analysis
- Candidate genes