Building a cluster of NLR genes conferring resistance to pests and pathogens: the story of the Vat gene cluster in cucurbits

Abstract

Most molecularly characterized plant resistance genes (R genes) belong to the nucleotide-binding-site-leucine-rich-repeat (NLR) receptor family and are prone to duplication and transposition with high sequence diversity. In this family, the Vat gene in melon is one of the few R genes known for conferring resistance to insect, i.e., Aphis gossypii, but it has been misassembled and/or mispredicted in the whole genomes of Cucurbits. We examined 14 genomic regions (about 400 kb) derived from long-read assemblies spanning Vat-related genes in Cucumis melo, Cucumis sativus, Citrullus lanatus, Benincasa hispida, Cucurbita argyrosperma, and Momordica charantia. We built the phylogeny of those genes. Investigating the paleohistory of the Vat gene cluster, we revealed a step by step process beginning from a common ancestry in cucurbits older than 50 my. We highlighted Vat exclusively in the Cucumis genera, which diverged about 20 my ago. We then focused on melon, evaluating a minimum duplication rate of Vat in 80 wild and cultivated melon lines using generalist primers; our results suggested that duplication started before melon domestication. The phylogeny of 44 Vat-CDS obtained from 21 melon lines revealed gain and loss of leucine-rich-repeat domains along diversification. Altogether, we revealed the high putative recognition scale offered in melon based on a combination of SNPs, number of leucine-rich-repeat domains within each homolog and number of homologs within each cluster that might jointly confer resistance to a large pest and pathogen spectrum. Based on our findings, we propose possible avenues for breeding programs.

References

[1] Zhang, R. Z., Murat, F., Pont, C., Langin, T. & Salse, J. Paleo-evolutionary plasticity of plant disease resistance genes. BMC Genomics 15, https://doi.org/10.1186/1471-2164-15-187 (2014).

[2] Kapos, P., Devendrakumara, K. T. & Li, X. Plant NLRs: from discovery to application. Plant Sci. 279, 3-18 (2019).

[3] Dangl, J. L., Horvath, D. M. & Staskawicz, B. J. Pivoting the plant immune system from dissection to deployment. Science 341, 746-751 (2013).

[4] Dogimont, C., Chovelon, V., Pauquet, J., Boualem, A. & Bendahmane, A. The Vat locus encodes for a CC-NBS-LRR protein that confers resistance to Aphis gossypii infestation and A. gossypii-mediated virus resistance. Plant J. 80, 993-1004 (2014).

[5] Boissot, N., Thomas, S., Chovelon, V. & Lecoq, H. NBS-LRR-mediated resistance triggered by aphids: viruses do not adapt; aphids adapt via different mechanisms. BMC Plant Biol. 16, 25 (2016).

[6] Boissot, N., Schoeny, A. & Vanlerberghe-Masutti, F. Vat, an amazing gene conferring resistance to aphids and viruses they carry: from molecular structure to field effects. Front. Plant Sci. 7, 1420 (2016).

[7] Maekawa, T., Kufer, T. A. & Schulze-Lefert, P. NLR functions in plant and animal immune systems: so far and yet so close. Nat. Immunol. 12, 818-826 (2011).

[8] Dogimont, C. et al. Molecular diversity at the Vat/Pm-W resistance locus in melon. in Cucurbitaceae 2008, Proceedings of the IXth EUCARPIA meeting on genetics and breeding of Cucurbitaceae (ed. Pitrat, M.) 219-227 (Avignon, 2008).

[9] Gonzalez, V., Aventin, N., Centeno, E. & Puigdomenech, P. Interspecific and intraspecific gene variability in a 1-Mb region containing the highest density of NBS-LRR genes found in the melon genome. BMC Genomics 15, 1131 (2014).

[10] Garcia-Mas, J. et al. The genome of melon (Cucumis melo L.). Proc. Natl Acad. Sci. USA 109, 11872-11877 (2012).

[11] Wu, S. et al. The bottle gourd genome provides insights into Cucurbitaceae evolution and facilitates mapping of a Papaya ring spot virus resistance locus. Plant J. 92, 963-975 (2017).

[12] Sun, H. H. et al. Karyotype stability and unbiased fractionation in the paleo-allotetraploid cucurbita genomes. Mol. Plant 10, 1293-1306 (2017).

[13] Montero-Pau, J. et al. De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus. Plant Biotechnol. J. 16, 1161-1171 (2018).

[14] Toda, N. et al. NLGenomeSweeper: a tool for genome-wide NBS-LRR resistance gene identification. Genes 11, 333 (2020).

[15] Thimmappa, R., Geisler, K., Louveau, T., O’Maille, P. & Osbourn, A. Triterpene biosynthesis in plants. Annu. Rev. Plant Biol. 65, 225-257 (2014).

[16] Quesneville, H. Twenty years of transposable element analysis in the Arabidopsis thaliana genome. Mobile DNA 11, 28 (2020).

[17] Castanera, R., Ruggieri, V., Pujol, M., Garcia-Mas, J. & Casacuberta, J. M. An improved melon reference genome with single-molecule sequencing uncovers a recent burst of transposable elements with potential impact on genes. Front. Plant Sci. 10, https://doi.org/10.3389/fpls.2019.01815 (2020).

[18] Sebastian, P., Schaefer, H., Telford, I. R. H. & Renner, S. S. Cucumber (Cucumis sativus) and melon (C. melo) have numerous wild relatives in Asia and Australia, and the sister species of melon is from Australia. Proc. Natl Acad. Sci. USA 107, 14269-14273 (2010).

[19] Hedges, S. B., Marin, J., Suleski, M., Paymer, M. & Kumar, S. Tree of life reveals clock-like speciation and diversification. Mol. Biol. Evol. 32, 835-845 (2015).

[20] Thomas, S., Boissot, N. & Vanlerberghe-Masutti, F. What do spring migrants reveal about sex and host selection in the melon aphid? BMC Evol. Biol. 12, https://doi.org/10.1186/1471-2148-12-47 (2012).

[21] Dixon, M. S., Hatzixanthis, K., Jones, D. A., Harrison, K. & Jones, J. D. G. The tomato Cf-5 disease resistance gene and six homologs show pronounced allelic variation in leucine-rich repeat copy number. Plant Cell 10, 1915-1925 (1998).

[22] Anderson, P. A. et al. Inactivation of the flax rust resistance gene M associated with loss of a repeated unit within the leucine-rich repeat coding region. Plant Cell 9, 641-651 (1997).

[23] Seifi, A. et al. Linked, if not the same, Mi-1 homologues confer resistance to tomato powdery mildew and root-knot nematodes. Mol. Plant-Microbe Interact. 24, 441-450 (2011).

[24] Seah, S., Telleen, A. C. & Williamson, V. M. Introgressed and endogenous Mi-1 gene clusters in tomato differ by complex rearrangements in flanking sequences and show sequence exchange and diversifying selection among homologues. Theor. Appl. Genet. 114, 1289-1302 (2007).

[25] Diaz, A. et al. A consensus linkage map for molecular markers and Quantitative Trait Loci associated with economically important traits in melon (Cucumis melo L.). BMC Plant Biol. 11, (2011).

[26] Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945-959 (2000).

[27] Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611-2620 (2005).

[28] Zhang, H. et al. A high-quality melon genome assembly provides insights into genetic basis of fruit trait improvement. Iscience 22, 16 (2019).

[29] Yano, R. et al. Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression. Commun. Biol. 3, 432 (2020).

[30] Yang, J. et al. The chromosome-scale genome of melon dissects genetic architecture of important agronomic traits. iScience 23, https://doi.org/10.1016/j.isci.2020.101422 (2020).

[31] Lagnel, J., Tsigenopoulos, C. S. & Iliopoulos, I. NOBLAST and JAMBLAST: New options for BLAST and a Java application manager for BLAST results. Bioinformatics 25, 824-826 (2009).

[32] Li, H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094-3100 (2018).

[33] Li, Q. et al. A chromosome-scale genome assembly of cucumber (Cucumis sativus L.). Gigascience 8, https://doi.org/10.1093/gigascience/giz072 (2019).

[34] Guo, S. G. et al. Resequencing of 414 cultivated and wild watermelon accessions identifies selection for fruit quality traits. Nat. Genet. 51, 1616 (2019).

[35] Xie, D. S. et al. The wax gourd genomes offer insights into the genetic diversity and ancestral cucurbit karyotype. Nat. Commun. 10, https://doi.org/10.1038/s41467-019-13185-3 (2019).

[36] Barrera-Redondo, J. et al. The genome of Cucurbita argyrosperma (silver-seed gourd) reveals faster rates of protein-coding gene and long noncoding RNA turnover and neofunctionalization within cucurbita. Mol. Plant 12, 506-520 (2019).

[37] Zhang, T. et al. Long-read sequencing and de novo assembly of the Luffa cylindrica (L.) Roem. genome. Mol. Ecol. Resour. 20, https://doi.org/10.1111/1755-0998.13129 (2020).

[38] Matsumura, H. et al. Long-read bitter gourd (Momordica charantia) genome and the genomic architecture of domestication. bioRxvi https://doi.org/10.1101/822411 (2019).

[39] Darling, A. E., Mau, B. & Perna, N. T. progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS ONE 5, https://doi.org/10.1371/journal.pone.0011147 (2010).

[40] Robinson, J. T. et al. Integrative genomics viewer. Nat. Biotechnol. 29, 24-26 (2011).

[41] Harris, R. S. Improved Pairwise Alignment of Genomic DNA (The Pennsylvania State University, 2007).

[42] Gouy, M., Guindon, S. & Gascuel, O. SeaView Version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol. Biol. Evolution 27, 221-224 (2010).

[43] Kumar, S., Stecher, G., Li, M., Knyaz, C. & Tamura, K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evolution 35, 1547-1549 (2018).

[44] Murrell, B. et al. FUBAR: a fast, unconstrained bayesian approximation for inferring selection. Mol. Biol. Evolution 30, 1196-1205 (2013).

[45] Campbell, M. S., Holt, C., Moore, B. & Yandell, M. Genome annotation and curation using MAKER and MAKER-P. Curr. Protoc. Bioinformatics 48, 4.11.1-39 (2014).