Patterns of genomic divergence in sympatric and allopatric speciation of three Mihoutao (Actinidia) species

Abstract

Isolation by geographic distance is involved in the formation of potential genomic islands and the divergence of genomes, which often result in speciation. The mechanisms of sympatric and allopatric speciation associated with geographic distance remain a topic of interest to evolutionary biologists. Here, we examined genomic divergence in three Actinidia species from large-scale sympatric and allopatric regions. Genome sequence data revealed that hexaploid Actinidia deliciosa originated from Actinidia chinensis and supported the speciation-with-gene-flow model in sympatric regions. The common ancestor of Actinidia setosa and A. deliciosa migrated from the mainland to the Taiwan Island ∼2.91 Mya and formed A. setosa ∼0.92 Mya, and the speciation of A. setosa is consistent with the divergence-after-speciation model with selective sweeps. Geographic isolation resulted in population contraction and accelerated the process of lineage sorting and speciation due to natural selection. Genomic islands contained genes associated with organ development, local adaptation, and stress resistance, indicating selective sweeps on a specific set of traits. Our results highlight the patterns of genomic divergence in sympatric and allopatric speciation, with the mediation of geographic isolation in the formation of genomic islands during Actinidia speciation.

References

[1] Veltsos P, Ritchie MG . Sympatric Speciation. In: eLS . Chichester: John Wiley & Sons, Ltd, 2011.

[2] Coyne JA . Sympatric speciation. Curr Biol. 2007; 17: 787-8.

[3] Fahrig L . Ecological responses to habitat fragmentation per se. Annu Rev Ecol Evol Syst. 2017; 48: 1-23.

[4] Sexton JP, Hangartner SB, Hoffmann AA . Genetic isolation by environment or distance: which pattern of gene flow is most common? Evolution . 2014; 68: 1-15.

[5] Rangel-Villalobos H, Martínez-Sevilla VM, Martínez-Cortés G et al. Importance of the geographic barriers to promote gene drift and avoid pre- and post-Columbian gene flow in Mexican native groups: evidence from forensic STR loci. Am J Phys Anthropol. 2016; 160: 298-316.

[6] Wolf JBW, Ellegren H . Making sense of genomic islands of differentiation in light of speciation. Nat Rev Genet. 2016; 18: 87-100.

[7] Van Doren BM, Campagna L, Helm B et al. Correlated patterns of genetic diversity and differentiation across an avian family. Mol Ecol. 2017; 26: 3982-97.

[8] Han F, Lamichhaney S, Grant BR et al. Gene flow, ancient polymorphism, and ecological adaptation shape the genomic landscape of divergence among Darwin’s finches. Genome Res. 2017; 27: 1004-15.

[9] Wang B et al. Ancient polymorphisms contribute to genome-wide variation by long-term balancing selection and divergent sorting in Boechera stricta . Genome Bio. 2019; 20: 126.

[10] Stankowski S, Chase MA, Fuiten AM et al. Widespread selection and gene flow shape the genomic landscape during a radiation of monkeyflowers. PLoS Biol. 2019; 17: e3000391.

[11] Ma T, Wang K, Hu Q et al. Ancient polymorphisms and divergence hitchhiking contribute to genomic islands of divergence within a poplar species complex. Proc Natl Acad Sci U S A. 2018; 115: E236-43.

[12] Cruickshank TE, Hahn MW . Reanalysis suggests that genomic islands of speciation are due to reduced diversity, not reduced gene flow. Mol Ecol. 2014; 23: 3133-57.

[13] Hazzi NA, Moreno JS, Ortiz-Movliav C, Palacio RD . Biogeographic regions and events of isolation and diversification of the endemic biota of the tropical Andes. Proc Natl Acad Sci U S A. 2018; 115: 7985-90.

[14] Craw D, Upton P, Burridge CP et al. Rapid biological speciation driven by tectonic evolution in New Zealand. Nat Geo. 2015; 9: 140-4.

[15] Cornejo-Romero A, Medina-Sánchez J, Hernández-Hernández T et al. Quaternary origin and genetic divergence of the endemic cactus Mammillaria pectinifera in a changing landscape in the Tehuacán Valley, Mexico. Gen Mol Res. 2014; 13: 73.

[16] Wang Q, Zhang ML, Yin LK . Genetic diversity and population differentiation of Capparis spinosa (Capparaceae) in northwestern China. Biochem Syst Ecol. 2016; 66: 1-7.

[17] Yang F-S, Qin A-L, Li Y-F, Wang X-Q . Great genetic differentiation among populations of Meconopsis integrifolia and its implication for plant speciation in the Qinghai-Tibetan plateau. PLoS One. 2012; 7: e37196.

[18] Crowhurst RN, Gleave AP, MacRae EA et al. Analysis of expressed sequence tags from Actinidia: applications of a cross species EST database for gene discovery in the areas of flavor, health, color and ripening. BMC Genomics. 2008; 9: 351.

[19] Liang C, Francisco AR . Emendation of the latin name of Actinidia chinensis Pl. Var Hispida C F Liang Guihaia. 1984; 4: 181-2.

[20] Liang C . Classification of Actinidia. J Syst Evol. 1975; 13: 32-5.

[21] Liang C, Francisco AR . Revision of infraspecific taxa of Actinidia chinensis planch. Guihaia. 1985; 5: 71-2.

[22] Huang H. Actinidia: Classification, Resources, Domestication and Cultivation. Beijing: Science China Press, 2014.

[23] He Z-C, Zhong Y, Liiu H-T, Tang X-H . Quantitative taxonomic analyses of Actinidia (Actinidiaceae) in China based on micromorphological characters of foliar trichomes. J Syst Evol. 2000; 38: 121-36.

[24] Chat J, Jáuregui B, Petit RJ, Nadot S . Reticulate evolution in kiwifruit (Actinidia, Actinidiaceae) identified by comparing their maternal and paternal phylogenies . Am J Bot. 2004; 91: 736.

[25] FOC. Flora of China . Beijing: Science China Press, 2007. http://www.efloras.org.

[26] Phivnil K, Beppu K, Takamura T et al. Flow cytometric assessment of ploidy in native resources of Actinidia in Japan. J Am Pomol Soc. 2005; 59: 44-9.

[27] Huang S, Ding J, Deng D et al. Draft genome of the kiwifruit Actinidia chinensis. Nat Commun. 2013; 4: 2640.

[28] Li H, Durbin R . Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009; 25: 1754-60.

[29] Li H, Handsaker B, Wysoker A et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009; 25: 2078-9.

[30] Liu Y, Li D, Zhang Q et al. Rapid radiations of both kiwifruit hybrid lineages and their parents shed light on a two-layer mode of species diversification. New Phytol. 2017; 215: 877-90.

[31] Wang K, Li M, Hakonarson H . ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010; 38: e164.

[32] Vilella AJ, Severin J, Ureta-Vidal A et al. EnsemblCompara GeneTrees: complete, duplication-aware phylogenetic trees in vertebrates. Genome Res. 2009; 19: 327-35.

[33] Tang H, Peng J, Wang P, Risch NJ . Estimation of individual admixture: analytical and study design considerations. Genet Epidemiol. 2005; 28: 289-301.

[34] Yang JA, Lee SH, Goddard ME, Visscher PM . GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011; 88: 76-82.

[35] Barrett JC, Fry B, Maller J, Daly MJ . Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005; 21: 263-5.

[36] Li H, Durbin R . Inference of human population history from individual whole-genome sequences. Nature. 2011; 475: 493-6.

[37] Liang Z, Duan S, Sheng J et al. Whole-genome resequencing of 472 Vitis accessions for grapevine diversity and demographic history analyses. Nat Commun. 2019; 10: 1190.

[38] An L, Zhou Z, Yan A, Gan Y . Progress on trichome development regulated by phytohormone signaling. Plant Signal Behav. 2014; 6: 1959-62.

[39] Bouckaert R, Heled J, Kühnert D et al. BEAST 2: a software platform for bayesian evolutionary analysis. PLoS Comput Biol. 2014; 10: e1003537.

[40] Helfrich P, Rieb E, Abrami G et al. TreeAnnotator: Versatile Visual Annotation of Hierarchical Text Relations. In: Proceedings of the Eleventh International Conference on Language Resources and Evaluation (LREC 2018), Miyazaki, Japan. European Language Resources Association (ELRA), 2018.

[41] Pickrell JK, Pritchard JK . Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012; 8: e1002967.

[42] Danecek P, Auton A, Abecasis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[43] Chan AH, Jenkins PA, Song YS . Genome-wide fine-scale recombination rate variation in Drosophila melanogaster . PLoS Genet. 2012; 8: e1003090.

[44] Ashburner M, Ball CA, Blake JA et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000; 25: 25-9.

[45] Zhou X, Stephens M . Genome-wide efficient mixed-model analysis for association studies. Nat Genet. 2012; 44: 821-4.

[46] Cui ZX . Actinidia in China . Shangdong, China: Shandong Science and Technology Press, 1993.

[47] Crowhurst RN, Gardner RC . A genome-specific repeat sequence from kiwifruit (Actinidia deliciosa var. deliciosa) . Theor Appl Genet. 1991; 81: 71-8.

[48] Ferguson AR, Huang H . Genetic resources of kiwifruit: domestication and breeding. Hortic Rev. 2007; 33: 1-121.

[49] Pilkington SM, Crowhurst R, Hilario E et al. A manually annotated Actinidia chinensis var. chinensis (kiwifruit) genome highlights the challenges associated with draft genomes and gene prediction in plants . BMC Genomics. 2018; 19: 257.

[50] Mallet J . Hybrid speciation. Nature. 2007; 446: 279.

[51] Liu Y, Li D, Yan L, Huang H . The microgeographical patterns of morphological and molecular variation of a mixed ploidy population in the species complex Actinidia chinensis . PLoS One. 2015; 10: e0117596.

[52] Huang B-H, Lin Y-C, Huang C-W et al. Differential genetic responses to the stress revealed the mutation-order adaptive divergence between two sympatric ginger species. BMC Genomics. 2018; 19: 692.

[53] Ellegren H, Galtier N . Determinants of genetic diversity. Nat Rev Genet. 2016; 17: 422-33.

[54] Kaplan N, Hudson R, Langley C . The hitchhiking effect revisited. Genetics. 1989; 123: 887-99.

[55] Aris-Brosou S, Excoffier L . The impact of population expansion and mutation rate heterogeneity on DNA sequence polymorphism. Mol Biol Evol. 1996; 13: 494-504.

[56] Choi JY, Purugganan M, Stacy EA . Divergent selection and primary gene flow shape incipient speciation of a riparian tree on Hawaii island. Mol Biol Evol. 2020; 37: 695-710.

[57] Chiang Y-C, Hung K-H, Moore S-J et al. Paraphyly of organelle DNAs in Cycas sect. Asiorientales due to ancient ancestral polymorphisms . BMC Evol Biol. 2009; 9: 161.

[58] Tavares H, Whibley A, Field DL et al. Selection and gene flow shape genomic islands that control floral guides. Proc Natl Acad Sci U S A. 2018; 115: 11006-11.

[59] Bijlsma R, Loeschcke V . Genetic erosion impedes adaptive responses to stressful environments. Evol Appl. 2012; 5: 117-29.

[60] Lienert J. Habitat fragmentation effects on fitness of plant populations-a review. J Nat Conserv. 2004; 12: 53-72.

[61] Browne L, Karubian J . Habitat loss and fragmentation reduce effective gene flow by disrupting seed dispersal in a neotropical palm. Mol Ecol. 2018; 27: 3055-69.

[62] Gao Y, Ai B, Kong H et al. Geographical pattern of isolation and diversification in karst habitat islands: a case study in the Primulina eburnea complex. J Biogeogr. 2015; 42: 2131-44.

[63] Wiberg R, Scobie A, Ahara S et al. The genetic consequences of long term habitat fragmentation on a self-incompatible clonal plant, Linnaea borealis L. Biol Conserv. 2016; 201: 405-13.

[64] Mes THM, Brederode J, Hart HT . Origin of the woody Macaronesian sempervivoideae and the phylogenetic position of the east African species of Aeonium . Plant Biol. 1996; 109: 477-91.

[65] Ravinet M, Yoshida K, Shigenobu S et al. The genomic landscape at a late stage of stickleback speciation: high genomic divergence interspersed by small localized regions of introgression. PLoS Genet. 2018; 14: e1007358.

[66] Rosche C, Schrieber K, Lachmuth S et al. Sex ratio rather than population size affects genetic diversity in Antennaria dioica . Plant Biol. 2018; 20: 789-96.

[67] Hermansen T, Minchinton T, Ayre D . Habitat fragmentation leads to reduced pollinator visitation, fruit production and recruitment in urban mangrove forests. Oecologia. 2017; 185: 221-31.

[68] Broeck A, Tobias C, Gunter K et al. Dispersal constraints for the conservation of the grassland herb Thymus pulegioides L. in a highly fragmented agricultural landscape. Conserv Genet. 2015; 16: 765-76.

[69] Llorens TM, Byrne M, Yates CJ et al. Evaluating the influence of different aspects of habitat fragmentation on mating patterns and pollen dispersal in the bird-pollinated Banksia sphaerocarpa var. caesia. Mol Ecol. 2012; 21: 314-28.

[70] Yu W, Wu B, Wang X et al. Scale-dependent effects of habitat fragmentation on the genetic diversity of Actinidia chinensis populations in China. Hortic Res. 2020; 7: 172.

[71] Kuussaari M, Bommarco R, Heikkinen RK et al. Extinction debt: a challenge for biodiversity conservation. Trends Ecol Evol. 2009; 24: 564-71.

[72] Pires HR, Monfared MM, Shemyakina EA, Fletcher JC . ULTRAPETALA trxG genes interact with KANADI transcription factor genes to regulate Arabidopsis gynoecium patterning. Plant Cell. 2014; 26: 4345-61.

[73] Ding ZJ, Yan JY, Li GX et al. WRKY41 controls Arabidopsis seed dormancy via direct regulation of ABI3 transcript levels not downstream of ABA. Plant J. 2014; 79: 810-23.

[74] Zhou Z, An L, Sun L, Gan Y . ZFP5 encodes a functionally equivalent GIS protein to control trichome initiation. Plant Signal Behav. 2012; 7: 28-30.

[75] Bischoff V, Nita S, Neumetzler L et al. TRICHOME BIREFRINGENCE and its homolog AT5G01360 encode plant-specific DUF231 proteins required for cellulose biosynthesis in Arabidopsis. Plant Physiol. 2010; 153: 590-602.

[76] Zhao M, Morohashi K, Hatlestad G et al. The TTG1-bHLH-MYB complex controls trichome cell fate and patterning through direct targeting of regulatory loci. Development. 2008; 135: 1991-9.

[77] Huchelmann A, Boutry M, Hachez C . Plant glandular trichomes: natural cell factories of high biotechnological interest. Plant Physiol. 2017; 175: 6-22.

[78] Yang C, Gao Y, Gao S et al. Transcriptome profile analysis of cell proliferation molecular processes during multicellular trichome formation induced by tomato Wo^v gene in tobacco. BMC Genomics. 2015; 16: 868.

[79] Vadde BVL, Challa KR, Nath U . The TCP4 transcription factor regulates trichome cell differentiation by directly activating GLABROUS INFLORESCENCE STEMS in Arabidopsis thaliana . Plant J. 2018; 93: 259-69.

[80] Garcia-Molina A, Xing S, Huijser P . A conserved KIN17 curved DNA-binding domain protein assembles with SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE7 to adapt Arabidopsis growth and development to limiting copper availability . Plant Physiol. 2014; 164: 828-40.

[81] Bae MS, Cho EJ, Choi EY, Park OK . Analysis of the Arabidopsis nuclear proteome and its response to cold stress. Plant J. 2003; 36: 652-63.

[82] Bayer EM, Bottrill AR, Walshaw J et al. Arabidopsis cell wall proteome defined using multidimensional protein identification technology . Proteomics. 2006; 6: 301-11.