First genetic linkage map of Lathyrus cicera based on RNA sequencing-derived markers: Key tool for genetic mapping of disease resistance

Abstract

The Lathyrus cicera transcriptome was analysed in response to rust (Uromyces pisi) infection to develop novel molecular breeding tools with potential for genetic mapping of resistance in this robust orphan legume species. One RNA-seq library each was generated from control and rust-inoculated leaves from two L. cicera genotypes with contrasting quantitative resistance, de novo assembled into contigs and sequence polymorphisms were identified. In toto, 19,224 SNPs differentiate the susceptible from the partially resistant genotype’s transcriptome. In addition, we developed and tested 341 expressed E-SSR markers from the contigs, of which 60.7% varied between the two L. cicera genotypes. A first L. cicera linkage map was created using part of the developed markers in a RIL population from the cross of the two genotypes. This map contains 307 markers, covered 724.2 cM and is organised in 7 major and 2 minor linkage groups, with an average mapping interval of 2.4 cM. The genic markers also enabled us to compare their position in L. cicera map with the physical position of the same markers mapped on Medicago truncatula genome, highlighting a high macrosyntenic conservation between both species. This study provides a large new set of genic polymorphic molecular markers with potential for mapping rust resistances. It represents the first step towards genomics-assisted precision breeding in L. cicera.

References

[1] Kenicer, G. J., Kajita, T., Pennington, R. T. & Murata, J. Systematics and biogeography of Lathyrus (Leguminosae) based on internal transcribed spacer and cpDNA sequence data. Am. J. Bot. 92, 1199-1209 (2005).

[2] Schaefer, H. et al. Systematics, biogeography, and character evolution of the legume tribe Fabeae with special focus on the middle-Atlantic island lineages. BMC Evol. Biol. 12, 250 (2012).

[3] Hanbury, C. D., Siddique, K. H. M., Galwey, N. W. & Cocks, P. S. Genotype-environment interaction for seed yield and ODAP concentration of Lathyrus sativus L. and L. cicera L. in Mediterranean-type environments. Euphytica 110, 45-60 (1999).

[4] Peña-Chocarro, L. & Peña, L. Z. History and traditional cultivation of Lathyrus sativus L. and Lathyrus cicera L. in the Iberian peninsula. Veg. Hist. Archaeobot. 8, 49-52 (1999).

[5] Silvestre, S., de Sousa Araújo, S., Vaz Patto, M. C. & Marques da Silva, J. Performance index: an expeditious tool to screen for improved drought resistance in the Lathyrus genus. J. Integr. Plant Biol. 56, 610-621 (2014).

[6] Vaz Patto, M. C., Fernández-Aparicio, M., Moral, A. & Rubiales, D. Pre and posthaustorial resistance to rusts in Lathyrus cicera L.. Euphytica 165, 27-34 (2009).

[7] Vaz Patto, M. C., Fernández-Aparicio, M., Moral, A. & Rubiales, D. Resistance reaction to powdery mildew (Erysiphe pisi) in a germplasm collection of Lathyrus cicera from Iberian origin. Genet. Resour. Crop Evol. 54, 1517-1521 (2007).

[8] Martín-Sanz, A., Pérez de la Vega, M. & Caminero, C. Resistance to Pseudomonas syringae in a collection of pea germplasm under field and controlled conditions. Plant Pathol. 61, 375-387 (2012).

[9] Fernández-Aparicio, M., Flores, F. & Rubiales, D. Field response of Lathyrus cicera germplasm to crenate broomrape (Orobanche crenata). Field Crops Res. 113, 321-327 (2009).

[10] Vaz Patto, M. C. & Rubiales, D. Lathyrus diversity: Available resources with relevance to crop improvement - L. sativus and L. cicera as case studies. Ann. Bot. 113, 895-908 (2014).

[11] Almeida, N. F. et al. Transferability of molecular markers from major legumes to Lathyrus spp. for their application in mapping and diversity studies. Mol. Biol. Rep. 41, 269-283 (2014).

[12] Barilli, E., Sillero, J. C., Serrano, A. & Rubiales, D. Differential response of pea (Pisum sativum) to rusts incited by Uromyces viciae-fabae and U. pisi. Crop Prot. 28, 980-986 (2009).

[13] Vaz Patto, M. C. & Rubiales, D. Identification and characterization of partial resistance to rust in a germplasm collection of Lathyrus sativus L.. Plant Breed. 128, 495-500 (2009).

[14] Barilli, E., Moral, A., Sillero, J. C. & Rubiales, D. Clarification on rust species potentially infecting pea (Pisum sativum L.) crop and host range of Uromyces pisi (Pers.) Wint.. Crop Prot. 37, 65-70 (2012).

[15] Rubiales, D., Sillero, J. C. & Emeran, A. A. Response of vetches (Vicia spp.) to specialized forms of Uromyces vicia-fabae and to Uromyces pisi. Crop Prot. 46, 38-43 (2013).

[16] Hopf, M. in The Origin and domestication of cultivated plants, Vol. 16 (ed Barigozzi, C.) Pages 35-60 (Elsevier Science, 1986).

[17] Almeida, N. F. et al. Allelic diversity in the transcriptomes of contrasting rust-infected genotypes of Lathyrus sativus, a lasting resource for smart breeding. BMC Plant Biol. 14, 376 (2014).

[18] Thimm, O. et al. MAPMAN: A user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914-939 (2004).

[19] Mayer, C. Phobos 3.3.11. http://www.rub.de/ecoevo/cm/cm_phobos.htm (2010).

[20] Xu, Y., Zhu, L., Xiao, J., Huang, N. & McCouch, S. R. Chromosomal regions associated with segregation distortion of molecular markers in F2, backcross, doubled haploid, and recombinant inbred populations in rice (Oryza sativa L.). Mol. Gen. Genet. 253, 535-545 (1997).

[21] Fernández-Aparicio, M. & Rubiales, D. Characterisation of resistance to crenate broomrape (Orobanche crenata Forsk.) in Lathyrus cicera L.. Euphytica 173, 77-84 (2010).

[22] Varshney, R. K., Graner, A. & Sorrells, M. E. Genic microsatellite markers in plants: features and applications. Trends Biotechnol. 23, 48-55 (2005).

[23] Zhang, X. & Borevitz, J. O. Global analysis of allele-specific expression in Arabidopsis thaliana. Genetics 182, 943-954 (2009).

[24] Vaz Patto, M. C. & Rubiales, D. Unveiling common responses of Medicago truncatula to appropriate and inappropriate rust species. Front. Plant Sci. 5, 618 (2014).

[25] Jørgensen, I. H. Discovery, characterization and exploitation of Mlo powdery mildew resistance in barley. Euphytica 63, 141-152 (1992).

[26] Chen, Z. et al. Expression analysis of the AtMLO gene family encoding plant-specific seven-transmembrane domain proteins. Plant Mol. Biol. 60, 583-597 (2006).

[27] Humphry, M., Reinstädler, A., Ivanov, S., Bisseling, T. & Panstruga, R. Durable broad-spectrum powdery mildew resistance in pea er1 plants is conferred by natural loss-of-function mutations in PsMLO1. Mol. Plant Pathol. 12, 866-878 (2011).

[28] Hacquard, S. et al. The poplar-poplar rust interaction: insights from genomics and transcriptomics. J. Pathog. 2011, 1-11 (2011).

[29] Cheng, R., Kleinhofs, A. & Ukai, Y. Method for mapping a partial lethal-factor locus on a molecular-marker linkage map of a backcross and doubled-haploid population. Theor. Appl. Genet. 97, 293-298 (1998).

[30] Törjék, O. et al. Segregation distortion in Arabidopsis C24/Col-0 and Col-0/C24 recombinant inbred line populations is due to reduced fertility caused by epistatic interaction of two loci. Theor. Appl. Genet. 113, 1551-1561 (2006).

[31] Aubert, G. et al. Functional mapping in pea, as an aid to the candidate gene selection and for investigating synteny with the model legume Medicago truncatula. Theor. Appl. Genet. 112, 1024-1041 (2006).

[32] Choi, H. K. et al. A sequence-based genetic map of Medicago truncatula and comparison of marker colinearity with M. sativa. Genetics 166, 1463-1502 (2004).

[33] Phan, H. T. T. et al. Extensive macrosynteny between Medicago truncatula and Lens culinaris ssp. culinaris. Theor. Appl. Genet. 114, 549-558 (2007).

[34] Weeden, N. F., Muehlbauer, F. J. & Ladizinski, G. Extensive conservation of linkage relationships between pea and lentil genetic maps. J. Hered. 83, 123-129 (1992).

[35] Kaló, P. et al. Comparative mapping between Medicago sativa and Pisum sativum. Mol. Genet. Genom. 272, 235-46 (2004).

[36] Choi, H.-K. et al. Estimating genome conservation between crop and model legume species. Proc. Natl Acad. Sci. 101, 15289-15294 (2004).

[37] Young, N. D. et al. Sequencing the genespaces of Medicago truncatula and Lotus japonicus. Plant Physiol. 137, 1174-1181 (2005).

[38] Torres, A. M., Weeden, N. F. & Martín, A. Linkage among isozyme, RFLP and RAPD markers in Vicia faba. Theor. Appl. Genet. 85, 937-945 (1993).

[39] Huang, X. & Madan, A. CAP3: A DNA sequence assembly program. Genome Res. 9, 868-877 (1999).

[40] Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460-2461 (2010).

[41] Waterhouse, R. M. et al. BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 35, 543-548 (2017).

[42] Marioni, J. C., Mason, C. E., Mane, S. M., Stephens, M. & Gilad, Y. RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 18, 1509-1517 (2008).

[43] Wang, L., Feng, Z., Wang, X., Wang, X. & Zhang, X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 26, 136-138 (2009).

[44] Lohse, M. et al. Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data. Plant, Cell Environ. 37, 1250-1258 (2014).

[45] Untergasser, A. et al. Primer3-new capabilities and interfaces. Nucleic Acids Res. 40, e115 (2012).

[46] Andersen, C. L., Jensen, J. L. & Ørntoft, T. F. Normalization of real-time quantitative reverse transcription-{PCR} data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets TL - 64. Cancer Res 64 VN-r, 5245-5250 (2004).

[47] Roth, A. et al. JointSNVMix: a probabilistic model for accurate detection of somatic mutations in normal/tumour paired next-generation sequencing data. Bioinformatics 28, 907-913 (2012).

[48] Kwok, S., Chang, S. Y., Sninsky, J. J. & Wang, A. A guide to the design and use of mismatched and degenerate primers. Genome Res. 3, S39-47 (1994).

[49] Liu, J. et al. An improved allele-specific PCR primer design method for SNP marker analysis and its application. Plant Methods 8, 34 (2012).

[50] Drummond, A. J. et al. Geneious v5.4 http://www.geneious.com (2011).

[51] Kosambi, D. D. The geometric method in mathematical statistics. Am. Math. Mon. 51, 382 (1944).

[52] Tang, H. et al. An improved genome release (version Mt4.0) for the model legume Medicago truncatula. BMC Genom. 15, 312 (2014).

[53] Bayer, M. et al. Comparative visualization of genetic and physical maps with strudel. Bioinformatics 27, 1307-1308 (2011).

[54] Van Ooijen, J. W. JoinMap 4, Software for the calculation of genetic linkage maps in experimental populations. Kyazma B. V., Wageningen, Netherlands (2006).