Uneven selection pressure accelerating divergence of Populus and Salix

Abstract

Populus (poplars) and Salix (willows) are sister genera in the Salicaceae family that arise from a common tetraploid ancestor. The karyotypes of these two lineages are distinguished by two major interchromosomal and some minor intrachromosomal rearrangements, but which one is evolutionarily more primitive remains debatable. In this study, we compare the selection pressure acting on the paralogous genes resulting from salicoid duplication (PGRS) within and between the genomes of the two lineages. Purifying selection was determined to act more strongly on the PGRS in willow than on those in poplar, which would cause a faster loss of paralogous duplicates in willow. Therefore, Salix species are supposed to evolve faster than Populus species, which is consistent with the observation that the former are taxonomically and morphologically more diverse than the latter. In these two lineages, different autosomes were found to have been evolving into sex chromosomes. Examining the ω ratio and the PGRS in the sex determination regions in willow and poplar revealed higher convergent selection pressure and a faster loss of PGRS in the sex determination regions of both lineages. At the chromosome level, the sex chromosome in poplar is characterized by the lowest gene density among all chromosome members, while this feature is not observed on the sex chromosome in willow, suggesting that Populus species may inherit the more incipient sex chromosome from their progenitor. Taken together, Salix is supposed to be the nascent lineage arising from the additional round of genome reorganization that distinguishes the karyotypes of the two sister genera. In this study, assessment of ω ratios also detected a list of paralogous genes under unusual selection pressure, which could have special consequences for the adaptive evolution of Salicaceae species. In conclusion, the results of this study provide unique information for better understanding the genetic mechanism accelerating the divergence of these two closely related lineages.

References

[1] Otto, S. P. & Whitton, J. Polyploid incidence and evolution. Annu. Rev. Genet. 34, 401-437 (2000).

[2] Wood, T. E. et al. The frequency of polyploid speciation in vascular plants. P Natl. Acad. Sci. USA 106, 13875-13879 (2009).

[3] Cui, L. et al. Widespread genome duplications throughout the history of flowering plants. Genome Res. 16, 738-749 (2006).

[4] Jaillon, O. et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 449, 463-467 (2007).

[5] Jiao, Y. et al. Ancestral polyploidy in seed plants and angiosperms. Nature 473, 97-100 (2011).

[6] He, X. & Zhang, J. Rapid subfunctionalization accompanied by prolonged and substantial neofunctionalization in duplicate gene evolution. Genetics 169, 1157-1164 (2005).

[7] Force, A. & Lynch, M. Preservation of duplicate genes by complementary, degenerative mutations. Genetics 151, 1531-1545 (1999).

[8] Prince, V. E. & Pickett, F. B. Splitting pairs: the diverging fates of duplicated genes. Nat. Rev. Genet. 3, 827-837 (2002).

[9] Adams, K. L., Cronn, R., Percifield, R. & Wendel, J. F. Genes duplicated by polyploidy show unequal contributions to the transcriptome and organ-specific reciprocal silencing. P Natl. Acad. Sci. USA 100, 4649-4654 (2003).

[10] Ohno, S. (eds) Evolution by Gene Duplication (Berlin Heidelberg, Springer-Verlag Press, 1970).

[11] deWet, J. M. Origins of polyploids. Basic Life Sci. 13, 3-15 (1979).

[12] Liu, B. & Wendel, J. F. Epigenetic phenomena and the evolution of plant allopolyploids. Mol. Phylogenet. Evol. 29, 365-379 (2003).

[13] Dickmann, D. I. & Kuzovkina, J. in Poplars and willows in the world, “ Poplars and willows of the world, with emphasis on silviculturally important species” Working paper IPC/9-2 (eds Isebrands, J.G. & Richardson, J.) Ch. 2 (FAO, Rome, Italy, 2008).

[14] Tuskan, G. A. et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray) . Science 313, 1596-1604 (2006).

[15] Dai, X. et al. The willow genome and divergent evolution from poplar after the common genome duplication. Cell Res. 24, 1274-1277 (2014).

[16] Hou, J. et al. Major chromosomal rearrangements distinguish willow and poplar after the ancestral “Salicoid” genome duplication. Genome Biol. Evol. 8, 1868-1875 (2016).

[17] Blackburn, K. B. & Harrison, J. W. H. A preliminary account of the chromosomes and chromosome behaviour in the Salicaceae. Ann. Bot.-Lond. 38, 361-378 (1924).

[18] Town, C. D. et al. Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy . Plant Cell 18, 1348-1359 (2006).

[19] Velasco, R. et al. The genome of the domesticated apple (Malus×domestica Borkh.) . Nat. Genet. 42, 833-839 (2010).

[20] Cheng, F. et al. Deciphering the diploid ancestral genome of the mesohexaploid Brassica rapa . Plant Cell 25, 1541-1554 (2013).

[21] Pucholt, P., Rönnberg-Wästljung, A. C. & Berlin, S. Single locus sex determination and female heterogamety in the basket willow (Salix viminalis L.) . Heredity 114, 575-583 (2015).

[22] Yin, T. et al. Genome structure and emerging evidence of an incipient sex chromosome in Populus . Genome Res. 18, 422-430 (2008).

[23] Hou, J. et al. Different autosomes evolved into sex chromosomes in the sister genera of Salix and Populus . https://doi.org/10.1038/srep09076 (2015).

[24] Bradshaw, H. D. & Stettler, R. F. Molecular genetics of growth and development in Populus. I. Triploidy in hybrid poplars . Theor. Appl. Genet. 86, 301-307 (1993).

[25] Thibault, J. Nuclear DNA amount in pure species and hybrid willows (Salix): a flow cytometric investigation . Can. J. Bot. 76, 157-165 (1998).

[26] Ma, T. et al. Genomic insights into salt adaptation in a desert poplar. Nat. Commun. 4, 2797 (2013).

[27] Ma, J. et al. Genome sequence and genetic transformation of a widely distributed and cultivated poplar. Plant Biotech. J. https://doi.org/10.1111/pbi.12989 (2018).

[28] Zhou, R. et al. Characterization of a large sex determination region in Salix purpurea L. (Salicaceae) . Mol. Genet. Genom. https://doi.org/10.1007/s00438-018-1473-y (2018).

[29] Pál, C., Papp, B. & Lercher, M. J. An integrated view of protein evolution. Nat. Rev. Genet. 7, 337-348 (2006).

[30] Clark, A. G. et al. Inferring nonneutral evolution from human-chimp-mouse orthologous gene trios. Science 302, 1960-1963 (2003).

[31] Kondrashov, F. A., Rogozin, I. B., Wolf, Y. I. & Koonin, E. V. Selection in the evolution of gene duplications. Genome Boil. 3, 1-9 (2002).

[32] Rodgers-Melnick, E. et al. Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus . Genome Res. 22, 95-105 (2012).

[33] Gaudet, M. et al. Genetic linkage maps of Populus nigra L. including AFLPs, SSRs, SNPs, and sex trait . Tree Genet Genome 4, 25-36 (2008).

[34] Hanley, S. J., Mallott, M. D. & Karp, A. Alignment of a Salix linkage map to the Populus genomic sequence reveals macrosynteny between willow and poplar genomes . Tree Genet Genomes 3, 35-48 (2006).

[35] Berlin, S., Lagercrantz, U., von Arnold, S., Öst, T. & Rönnberg-Wästljung, A. C. High-density linkage mapping and evolution of paralogs and orthologs in Salix and Populus . BMC Genom. 11, 129 (2010).

[36] Lynch, M. & Conery, J. S. The evolutionary fate and consequences of duplicate genes. Science 290, 1151-1155 (2000).

[37] Eckenwalder, J. E. Systematics and evolution of Populus . in Biology of Populus and its Implications for Management and Conservation. (eds Stettler, R. F., Bradshaw, H. D., Jr., Heilman, P. E. & Hinckley, T. M.) ON. p. 7-32 (Ottawa, NRC Research Press, 1996).

[38] Argus, G. W. Infrageneric classification of Salix (Salicaceae) in the new world . Syst. Bot. Monogr. 52, 1-121 (1997).

[39] Skvortsov, A. K. Willows of Russia and adjacent countries. Taxonomical and Geographical Review (eds Zinovjev, A. G., Argus, G. W., Tahvanainen, J. & Roininen, H.). (Joensuu, University of Joensuu Press, 1999).

[40] Dorn, R. D. A synopsis of American Salix . Can. J. Bot. 54, 2769-2789 (1976).

[41] Markussen, T., Pakull, B. & Fladung, M. Positioning of sex-correlated markers for Populus in a AFLP-and SSR-Marker based genetic map of Populus tremula×tremuloides . Silvae Genet. 56, 180-184 (2007).

[42] Paolucci, I. et al. Genetic linkage maps of Populus alba L. and comparative mapping analysis of sex determination across Populus species . Tree Genet Genomes 6, 863-875 (2010).

[43] Pakull, B. et al. Genetic mapping of linkage group XIX and identification of sex-linked SSR markers in a Populus tremula×Populus tremuloides cross . Can. J. For. Res. 41, 245-253 (2011).

[44] Tuskan, G. A. et al. The obscure events contributing to the evolution of an incipient sex chromosome in Populus: a retrospective working hypothesis . Tree Genet Genomes 8, 559-571 (2012).

[45] Xin, H. et al. Chromosome painting and comparative physical mapping of the sex chromosomes in Populus tomentosa and Populus deltoides . Chromosoma 127, 313-321 (2018).

[46] Telgmann-Rauber, A., Jamsari, A., Kinney, M. S., Pires, J. C. & Jung, C. Genetic and physical maps around the sex-determining M-locus of the dioecious plant asparagus. Mol. Genet. Genom. 278, 221-234 (2007).

[47] Freeling, M. Bias in plant gene content following different sorts of duplication: tandem, whole-genome, segmental, or by transposition. Annu Rev. Plant Boil. 60, 433-453 (2009).

[48] Wu, P., Wang, W., Li, Y. & Hou, X. Divergent evolutionary patterns of the MAPK cascade genes in Brassica rapa and plant phylogenetics . Hortic Res. https://doi.org/10.1038/hortres.2017.79 (2017).

[49] Cheng, F., Wu, J. & Wang, X. Genome triplication drove the diversification of Brassica plants . Hortic Res. https://doi.org/10.1038/hortres.2014.24 (2014).

[50] Liu, S. et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes . Nat Commun. https://doi.org/10.1038/ncomms4930 (2014).

[51] Kimura, M. (eds) The Neutral Theory of Molecular Evolution (London, Cambridge Univ. Press, 1983).

[52] Altschul, S. F., Gish, W., Miller, W., Meyers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403-410 (1990).

[53] Tang, H. et al. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 18, 1944-1954 (2008).

[54] Wang, Y. et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, e49 (2012).

[55] Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792-1797 (2004).

[56] Suyama, M., Torrents, D. & Bork, P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 34, W609-W612 (2006).

[57] Wang, D., Zhang, Y., Zhang, Z., Zhu, J. & Yu, J. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genom. Proteom. Bioinf. 8, 77-80 (2010).

[58] Tallarida, R. J. & Murray, R. B. (eds) Mann-Whitney Test. in Manual of Pharmacologic Calculations (New York, Springer-Verlag Press, 1987).