Identifying the factors driving divergence is a key research topic of speciation genomics (Wolf and Ellegren, 2017). Although there has been progress in identification of individual speciation genes and genomic islands of divergence in speciation (Nosil and Feder, 2013; Renaut et al., 2013), it remains an open question as to whether there are general genomic patterns and processes of molecular evolution accompanying the divergence of species. Natural selection and drift are recognized as the major processes promoting divergence. However, the relative roles of these processes, the proportion of the genome affected by natural selection, and the strength of selection in generating genomic divergence have remained poorly understood (Kimura, 1979; Noor and Bennett, 2009; Barrett and Hoekstra, 2011; Nosil and Feder, 2013). A few recent studies have uncovered repeated patterns of genomic divergence driven by selection in species formation (Gagnaire et al., 2013; Arnegard et al., 2014; Renaut et al., 2014; Soria-Carrasco et al., 2014) and suggested that divergence may often result from changes at a relatively small subset of genes (e.g., Conte et al., 2012).

Geographic speciation is a common driving force in generating biodiversity (Coyne and Allen Orr, 1998). Given sufficient time and/or selection pressures, genomic divergence between allopatric populations is expected (Coyne and Allen Orr, 1998; Nosil, 2008). However, how genes diverge over time at a genomic scale after major geographic isolation that completely interrupts gene flow has remained an unanswered question. Meta-analyses of closely related species that span a range of time scales of geographic isolation will be, therefore, particularly beneficial for gaining insights into this question (Nosil et al., 2009). Such studies might allow us to reconstruct how genomic divergence unfolds as speciation proceeds through time.

Putatively orthologous (single-copy) genes (POGs) identified in RNA-seq data sets of closely related species offer a valuable genome-wide window on molecular divergence of a unique subset of the genome following speciation (Feder et al., 2012; De Smet et al., 2013). Investigating how these genes diverge during speciation and subsequent evolution contributes an important perspective to speciation genomics. Such studies may illuminate the functional divergence of the allopatric genomes and the underlying ecological forces that shape them. The forests of eastern Asia (EA) and eastern North America (ENA) share 65 genera of seed plants with closely related species occurring in the two areas (Li, 1952; Boufford and Spongberg, 1983; Wu, 1983; Wen, 1999; Wen et al., 2010; Wen et al., 2016). This well-known floristic disjunction is a major phytogeographic pattern of the Northern Hemisphere and has been known since the time of Linnaeus (Wen, 1999). The origins of the disjunction and impact on plant evolution have been investigated over the past 30 years (Xiang et al., 1998; Wen, 1999; Xiang et al., 2000; Wen, 2001; Xiang and Soltis, 2001; Donoghue and Smith, 2004; Harris et al., 2013; Manos and Meireles, 2015; Qian et al., 2017). Many of these genera are small with one to a few species. Previous analyses using mainly plastid and internal transcribed spacer (ITS) DNA sequence markers revealed that the divergence times of disjunct sister species in EA and ENA varied greatly from the Miocene to the Pleistocene or even more recently (15 mya to <2.0 mya), although some disjunct clades in the two areas diverged earlier (Xiang et al., 2000; Wen, 2001; Milne and Abbott, 2002; Donoghue and Smith, 2004; Harris et al., 2013). These genera therefore represent different stages of allopatric divergence and are ideal for studying genomic divergence associated with geographic speciation.

In this paper, we compare the divergence patterns of genes identified as POGs in leaf transcriptomes of 20 allopatric species/variety pairs (or taxon pairs), all but four of which represent the EA-ENA floristic disjunction (Table 1). The 16 EA-ENA disjunct species or subspecies/variety pairs represent diverse clades of flowering genera with two to several species: Campsis (two spp.) – Lamiales, Convallaria (three vars. of 1 sp.) – Asparagales, Cornus (two subclades, Cornus-1 with eight spp., Cornus-2 with two spp.) – Cornales, Cotinus (two spp.) – Sapindales, Croomia (three spp.) – Pandanales, Gelsemium (three spp.) – Gentianales, Hamamelis (six spp.) – Saxifagales, Liquidambar (4–15 spp.) – Saxifragales, Liriodendron (two spp.) – Magnoliales, Meehania (seven spp.) – Lamiales, Menispermum (two spp.) – Ranunculales, Nelumbo (two spp.) – Proteales, Penthorum (two spp.)– Saxifragales, Phryma (two spp. or two vars. of 1 sp.) – Lamiales, Sassafras (three spp.) – Magnoliales, and Saururus (two spp.) – Piperales (Appendix 1—figure 1). These EA-ENA disjunct genera are among the classic examples of the EA-ENA floristic disjunction (Li, 1952; Wu, 1983; Wen, 1999) and have one to a few species isolated in the two areas, and in a few cases, a species also occurs in western North America (e.g., Cornus-1) and/or southwestern Asia and southeastern Europe (e.g., Convallaria, Liquidambar, Nelumbo). The species pairs we compared in these genera (Table 1) span a range of divergence times (Xiang et al., 2000; Wen, 2001; Milne and Abbott, 2002; Wen et al., 2010).

As noted, four of the 20 pairs represent other allopatric pairs (i.e., Acorus, two spp. – Acorales, Calycanthus, three spp. – Laurales, Dysosma, seven spp. – Ranunculales). These additional pairs were included in the study to increase our sampling of geographically separated taxon pairs for comparison and to cover a wide window of times, from recent divergence of speciation to further divergence of the descendant species. The species pairs analyzed follow: one EA-Western North American (WNA) pair (Calycanthus chinensis – C. occidentalis), one Asia/S. Europe-ENA pair (Cotinus coggygria – C. obovatus), one EA-EA/NA pair (Acorus gramineus – A. calamus), and one Asian pair (Dysosma versipellis – D. pleiantha; Ranunculales) (Table 1). Species of Acorus and Dysosma partially overlap in their geographic ranges, but have diverged in habitat and do not occur together. Acorus calamus naturally occurs throughout China and adjacent countries, and in North America in swamps, pond sides, and standing water below 2800 m. Acorus gramineus occurs in eastern, southern, and western China, and adjacent countries of EA in dense forests, moist rocky stream banks, and meadows below 2600 m. Dysosma versipellis and D. pleiantha both occur in eastern, central, and southern China, but D. versipellis is found in deciduous forests while D. pleiantha occurs in evergreen forests.

Additional information on the 20 taxon pairs and their distributions can be obtained from the Flora of China (http://www.efloras.org) and the Flora of North America (http://floranorthamerica.org/). The 20 taxon pairs we compared span a range of divergence times (20-23, 26-28), mostly representing sister species or varieties; however, in a few cases, the selected species represent sister clades of one vs. two species (Calycanthus, Gelsemium, Sassafras), one vs. a few species (Hamamelis), or two or three vs. a few species (Cornus-1, Liquidambar).

We examined the pattern of variation of synonymous substitutions (Ks; the number of synonymous substitutions per synonymous site), nonsynonymous substitutions (Ka; the number of nonsynonymous substitutions per nonsynonymous site), and the ratio (Ka/Ks) for each POG within each taxon pair using a custom pipeline (Appendix 1—figure 2) and compared the patterns among taxon pairs representing these diverse angiosperm lineages. The ratio of nonsynonymous to synonymous nucleotide substitutions (Ka/Ks) is probably the most common measure of selection pressure and for identifying positive selection (Ka/Ks > 1) (Yang, 2003) in molecular evolutionary studies. The method is known to be a conservative approach to call positive selection on a gene as it uses an average Ka/Ks ratio across all nucleotide sites. We also examined the relative proportions of genes under different levels of selection forces and identify and annotate genes under strong positive selection based on the Ka/Ks values (>2). Finally, we tested if the level of molecular divergence at synonymous sites (as measured by Ks at the peak frequency) reflects the relative divergence time of the species pair. Results from such analyses will illuminate the genomic pattern of genetic divergence of these genes under geographic isolation, and may also shed light on a small window of genomic divergence that unfolds as speciation proceeds through time. Furthermore, this study will provide a ‘genome-wide’ view of evolutionary divergence of species exhibiting the EA-ENA phytogeographic pattern, in comparison to previous understanding based on a few gene regions. The POGs examined in this study represented 52–68% of the total orthologous gene families in the transcriptome data (Supplementary file 1).

We stress that the POGs in leaf transcriptomes analyzed here are putative orthologs between the two species of each biogeographic pair found in the transcriptomes by OrthoMCL followed by customized scripts (see Methods below). It is possible, therefore, that some of these genes may have more than one copy in the genome (as a small gene family), but only one of them was expressed to a detectable level in the leaf. In such cases, it is possible that the two single-copy transcripts from the taxon pair may represent paralogous copies of the low-copy gene family, making the gene comparison non-orthologous. However, we think the likelihood of non-orthology of these genes is probably very low because in each case the two taxa are either sister species or intraspecific taxa, or represent the products of recently diverged sister clades where conservation of expression of a gene copy may be expected. We therefore do not think that the low probability of non-homology will affect our conclusions derived from thousands of putatively orthologous comparisons. However, we do want to make it clear that the discussion and conclusions about ‘genome-wide’ pattern of evolutionary divergence and processes in this paper are restricted to these genes of the genome. Nonetheless these POGs represent a significant portion of the gene families expressed in the transcriptome of each species (52%–68%; Supplementary file 1).

Comparison of the rates of nonsynonymous and synonymous substitutions provides a useful tool for understanding the mechanisms of DNA sequence evolution. A ratio of nonsynonymous (Ka) to synonymous (Ks) nucleotide substitutions greater than one is a common method for identifying positive selection in molecular evolutionary studies (Yang, 2003). We are aware that uncertainty in Ka/Ks ratios is high when Ks is zero or close to zero. This uncertainty may lead to inflation of Ka/Ks ratios or false identification of genes under strong positive selection due to a high, infinite, or undefined Ka/Ks value. To minimize this problem, we removed genes with Ks = 0 and genes with Ka/Ks = 99 for all analyses and only focused on genes with Ka/Ks > 2 as candidates of strong positive selection for gene annotation analysis. Our detailed examination of the Ka and Ks data (see Discussion below) provided evidence supporting the candidate genes under positive selection detected in the study are unlikely artifacts of small Ks. Further, the Ka/Ks > 2 ratios in these genes are unlikely a result of saturation of synonymous substitutions. Saturation of synonymous substitutions is often expected at deeper phylogenetic levels. We further used alternative method (Wagner, 2007) to check for signals of positive selection in these genes with Ka/Ks > 2. This method detects variation clusters of aggregated nucleotide substitutions that are too closely spaced to be observed by chance alone (violating the prediction from neutral evolution). Both the Ka/Ks and variation cluster methods may underestimate the number of genes under positive selection because some genes with Ka/Ks < 1 may be under positive selection at only a few sites and these sites may not be necessarily clustered in the gene. However, our goal is not to determine exactly how many POGs are under positive selection, and this caveat does not affect our comparisons of average selection pressures among genes or comparisons of general patterns of gene divergence across species pairs. The genes identified as under positive selection by both methods are likely the true positive selection genes and will provide top candidates for future studies for their functions in driving allopatric divergence in these taxa.

The number of POGs in the leaf transcriptomes of the 20 biogeographic taxon pairs varies from 8596 to 15161 (Table 1). Among these POGs, 2241 are shared by over 50% of the taxa, 333 are shared by over 80% of the taxa, 79 are shared by over 90% of the taxa, and seven are shared by all of the taxa. Most of these genes exhibit low divergence, with synonymous substitution (Ks) values less than 0.1 and nonsynonymous substitution values (Ka) less than 0.04. Less than 1% of the genes exhibit greater divergence with Ks values greater than 0.15 (Figure 1) and Ka values greater than 0.05 (Figure 2). The Ks values at peak frequency when plotted in intervals of 0.01 range among the taxon pairs from <0.01 to<0.06 (Figure 1). In Convallaria, Liriodendron, Cotinus, Menispermum, and Campsis, the Ks values at peak frequency are less than 0.01, while in Cornus-2, Dysosma, Nelumbo, Penthorum, and Sassafras, they are between 0.01 and 0.02. In Calycanthus, Liquidambar, Meehania, Menispermum, and Phryma, they are between 0.02 and 0.03, and in Cornus-1, Croomia, Hamamelis, and Saururus, they are between 0.03 and 0.04. The largest values are observed in Acorus and Gelsemium, between 0.05 and 0.06 (Supplementary file 2). A portion of the POGs, 53–237 among the genera, were aligned to putative single-copy or low-copy genes shared by seed plants (Li et al., 2017), and 75–365 were mapped to strict single-copy genes shared by 20 angiosperm genomes (De Smet et al., 2013) (for method of the mapping analysis, see Methods).

Figure 1

Download asset Open asset

Figure 2

Download asset Open asset

The Ka/Ks ratios of the POGs exhibit high similarity in patterns of variation among the 20 taxon pairs. In all pairs, the ratios at peak frequencies are highly similar and more than 90% of the genes are overall under purifying selection (referred to as purifying selection in short) based on their average Ka/Ks ratios < 0.9 across all sites within the genes. The exception was Nelumbo, which has 87.6% of the genes with Ka/Ks ratios < 0.9. The greatest proportion of genes under purifying selection was observed in Gelsemium, which had 94.9% of the genes with Ka/Ks ratios < 0.9. A very small proportion, from 1.6% (in Acorus) to 3.7% (in Nelumbo), is evolving nearly neutrally with Ka/Ks ratios in the range of 0.9–1.1, and another small proportion, from 3.1% (in Gelsemium) to 8.8% (in Nelumbo), is under putative positive selection based on values of Ka/Ks ratios greater than 1.1. A very small proportion, from 0.6% (in Gelsemium) to 2.0% (in Nelumbo), is under putative strong positive selection with Ka/Ks ratios ranging from 2 to 7 and rarely to 17 (Figure 3). The proportion of genes under truly positive selection may be smaller the numbers reported here due to uncertainty in estimation of Ka/Ks ratios when Ks is small.

Figure 3

Download asset Open asset

We arbitrarily divided the genes under purifying selection into three categories, Ka/Ks < 0.1 as genes under strong purifying selection, 0.1 < Ka/Ks < 0.5 as genes under moderate purifying selection, and 0.5 < Ka/Ks < 0.9 as genes under relaxed purifying selection. Among the taxon pairs, we observed 41% (in Convallaria) to 68% (in Gelsemium) of the POGs fall into the category of moderate purifying selection, 11% (in Gelsemium) to 36% (in Convallaria) fall into the category of strong purifying selection, while 9.7% (in Acorus) to 19.5% (in Nelumbo) fall into the category of relaxed purifying selection. In all taxon pairs, the relative abundance of POGs under the different categories of selection pressures follows the same order, from high to low, as moderate purifying selection, strong purifying selection, relaxed purifying selection, weak/moderate positive selection, putatively neutral, and strong positive selection (Figure 4; Supplementary file 3).

Figure 4

Download asset Open asset

To investigate if variation in Ks values at peak frequency reflects variation in biogeographic isolation times of taxon pairs, phylogenetic analyses were conducted using the seven POGs shared by all taxa and the 79 POGs shared by >90% of the taxa, respectively. The results recovered the relationships of the 20 lineages congruent with the APG IV summary phylogeny and classification (The Angiosperm Phylogeny Group, 2016) (Figure 5; Figure 5—figure supplement 1). Divergence times of the biogeographic pairs estimated using total substitutions of these two sets of gene sequence data and BEAST (Drummond and Rambaut, 2007) were similar, with the results from 79 genes without the fossil constraint from Cornus slightly younger (Figure 5; Figure 5—figure supplement 1). The estimated divergence times suggested that isolation of the species pairs occurred in an interval between the late Miocene and Pleistocene (~10.67 –~1.69 mya; Figure 5). Species divergence time was positively correlated with the Ks value at peak frequency across taxon pairs (Figure 5; r = 0.8034; p<0.001; for values used for the correlation analysis, see Supplementary file 4). On the other hand, the Ks values at peak frequency of the taxon pairs are: (1) positively correlated with the abundance of genes under moderate purifying selection (0.1 < Ka/Ks <0.5; r = 0.9324; p<0.001), (2) negatively correlated with the abundance of genes under putative positive selection (Ka/Ks >=1.0, r = - 0.5344; p=0.015; 1.1 <= Ka/Ks <2, r = - 0.62; p=0.004), and (3) negatively correlated with the abundance of genes under putative strong purifying selection (Ka/Ks <0.1; r = - 0.7940; p<0.001) (Figure 6; Supplementary file 3, Supplementary file 10). The divergence times of taxon pairs are also positively correlated with the abundance of genes under moderate purifying selection (0.1 < Ka/Ks <0.5; r = 0.73; p<0.001) and negatively correlated with the abundance of genes under strong purifying selection (Ka/Ks <0.1; r = - 0.64; p<0.009; Figure 5; Supplementary file 4). The number of genes under strong positive selection (Ka/Ks >2) varies among the taxon pairs from 72 (in Campsis) to 184 (in Nelumbo) (Table 2). In each pair, approximately half of these genes have GO annotations (Table 2) that, in combination among all pairs, mapped to a total of 223 GO terms, varying from 21 (in Cornus-2) to 46 (in Convallaria) in each taxon pair (Supplementary file 5). Many of these GO terms are unique to a single lineage (100/223 = 44.8%).

Figure 5 with 1 supplement see all

Download asset Open asset

Figure 6

Download asset Open asset

A majority (59) of these lineage-specific GO terms belongs to the Biological Process (BP) categorization, while 23 belong to the Molecular Function (MF) ontology, and 16 belong to the Cell Components (CC) ontology (Supplementary file 5). Fourteen of the BP GO terms were annotated to responses to various stimuli (biotic or abiotic, external or internal stimuli), defense response, immune process, or signal transduction and were distributed among 11 taxon pairs, with each term was present in 1–3 pairs. In all, 124 of the 223 GO terms (55.2%) are shared by two or more of the 20 taxon pairs. Among these, 12 were shared by 10–15 taxon pairs. These 12 genes were annotated to the nucleus, protein complex, hydrolase activity, zinc ion binding, RNA binding, ATP binding, structural constituent of ribosome, oxidoreductase activity, cellular protein modification process, translation, single-organism biosynthetic process, and oxidation-reduction process (Supplementary file 5). A single GO term is shared by all 20 taxon pairs and is annotated to integral component of membrane (ICM) in the CC category (Supplementary file 5). The specific genes and the relative abundance of genes within each of the CC, BP, and MP categories vary among the taxon pairs (Figure 7; Figure 7—figure supplement 1; Figure 7—figure supplement 2). However, within the CC category, all pairs have the greatest proportion annotated to ICM (Figure 7), although they vary in the exact number (Table 3). The analysis of the Ka/Ks > 2 genes using the variation cluster method (Wagner, 2007) showed most of them (>57% in all taxon pairs, >70% in 15 pairs) were also under positive selection (indicated by significant p value for test of clustered nucleotide substitutions) (Supplementary file 6).

Figure 7 with 2 supplements see all

Download asset Open asset

Sequences of ortologous gene families and pairs of POGs sequences (prior to EST Scan) used for calculation of Ka and Ks have been submitted to Dryad (https://doi.org/10.5061/dryad.f1f0q44). Raw transcriptome data have been submitted to NCBI SRA database with BioProject number PRJNA508825 and BioSample number from SAMN10534244 to SAMN10534283 (Supplementary File 11).

1. 1. Dong YB
   2. Chen SC
   3. Cheng SF
   4. Zhou WB
   5. Ma Q
   6. Chen ZD
   7. Fu CX
   8. Liu X
   9. Zhao YP
   10. Soltis PS
   11. Soltis DE
   12. Xiang QY

   (2019) NCBI BioProject

   ID PRJNA508825. Natural selection and repeated patterns of molecular evolution following allopatric divergence.

   https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA508825
2. 1. Dong Y
   2. Chen S
   3. Cheng S
   4. Zhou W
   5. Ma Q
   6. Chen Z
   7. Fu C
   8. Liu X
   9. Zhao Y
   10. Soltis PS
   11. Wong GK
   12. Soltis DE
   13. Xiang Q

   (2019) Dryad Digital Repository

   Data from: Natural selection and repeated genome-wide patterns of molecular evolution following allopatric divergence.

   https://doi.org/10.5061/dryad.f1f0q44

1. 1. Altenhoff AM
   2. Dessimoz C

   (2009) Phylogenetic and functional assessment of orthologs inference projects and methods

   PLOS Computational Biology 5:e1000262.

   https://doi.org/10.1371/journal.pcbi.1000262

   • PubMed
   • Google Scholar
2. 1. Arnegard ME
   2. McGee MD
   3. Matthews B
   4. Marchinko KB
   5. Conte GL
   6. Kabir S
   7. Bedford N
   8. Bergek S
   9. Chan YF
   10. Jones FC
   11. Kingsley DM
   12. Peichel CL
   13. Schluter D

   (2014) Genetics of ecological divergence during speciation

   Nature 511:307–311.

   https://doi.org/10.1038/nature13301

   • PubMed
   • Google Scholar
3. 1. Atkinson BA
   2. Stockey RA
   3. Rothwell GW

   (2016) Cretaceous origin of dogwoods: an anatomically preserved Cornus (Cornaceae) fruit from the campanian of Vancouver island

   PeerJ 4:e2808.

   https://doi.org/10.7717/peerj.2808

   • PubMed
   • Google Scholar
4. 1. Barrett RD
   2. Hoekstra HE

   (2011) Molecular spandrels: tests of adaptation at the genetic level

   Nature Reviews Genetics 12:767–780.

   https://doi.org/10.1038/nrg3015

   • PubMed
   • Google Scholar
5. 1. Boufford DE
   2. Spongberg SA

   (1983) Eastern Asian-Eastern north american phytogeographical Relationships-A history from the time of linnaeus to the twentieth century

   Annals of the Missouri Botanical Garden 70:423–439.

   https://doi.org/10.2307/2992081

   • Google Scholar
6. 1. Chen F
   2. Mackey AJ
   3. Vermunt JK
   4. Roos DS

   (2007) Assessing performance of orthology detection strategies applied to eukaryotic genomes

   PLOS ONE 2:e383.

   https://doi.org/10.1371/journal.pone.0000383

   • PubMed
   • Google Scholar
7. 1. Conesa A
   2. Götz S
   3. García-Gómez JM
   4. Terol J
   5. Talón M
   6. Robles M

   (2005) Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research

   Bioinformatics 21:3674–3676.

   https://doi.org/10.1093/bioinformatics/bti610

   • PubMed
   • Google Scholar
8. 1. Conte GL
   2. Arnegard ME
   3. Peichel CL
   4. Schluter D

   (2012) The probability of genetic parallelism and convergence in natural populations

   Proceedings of the Royal Society B: Biological Sciences 279:5039–5047.

   https://doi.org/10.1098/rspb.2012.2146

   • Google Scholar
9. 1. Coyne JA
   2. Allen Orr H

   (1998) The evolutionary genetics of speciation

   Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences 353:287–305.

   https://doi.org/10.1098/rstb.1998.0210

   • Google Scholar
10. 1. De Smet R
    2. Adams KL
    3. Vandepoele K
    4. Van Montagu MC
    5. Maere S
    6. Van de Peer Y

    (2013) Convergent gene loss following gene and genome duplications creates single-copy families in flowering plants

    PNAS 110:2898–2903.

    https://doi.org/10.1073/pnas.1300127110

    • PubMed
    • Google Scholar
11. 1. De Vos JM
    2. Joppa LN
    3. Gittleman JL
    4. Stephens PR
    5. Pimm SL

    (2015) Estimating the normal background rate of species extinction

    Conservation Biology 29:452–462.

    https://doi.org/10.1111/cobi.12380

    • PubMed
    • Google Scholar
12. Software

    1. Dong Y

    (2018) GDMETS, version 1

    GitHub.

    https://github.com/ybdong919/GDMETS/releases
13. 1. Donoghue MJ
    2. Smith SA

    (2004) Patterns in the assembly of temperate forests around the northern hemisphere

    Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences 359:1633–1644.

    https://doi.org/10.1098/rstb.2004.1538

    • Google Scholar
14. 1. Drummond AJ
    2. Ho SY
    3. Phillips MJ
    4. Rambaut A

    (2006) Relaxed phylogenetics and dating with confidence

    PLOS Biology 4:e88.

    https://doi.org/10.1371/journal.pbio.0040088

    • PubMed
    • Google Scholar
15. 1. Drummond AJ
    2. Rambaut A

    (2007) BEAST: bayesian evolutionary analysis by sampling trees

    BMC Evolutionary Biology 7:214.

    https://doi.org/10.1186/1471-2148-7-214

    • PubMed
    • Google Scholar
16. 1. E. Ricklefs R
    2. Qian H
    3. S. White P

    (2004) The region effect on mesoscale plant species richness between eastern Asia and eastern north america

    Ecography 27:129–136.

    https://doi.org/10.1111/j.0906-7590.2004.03789.x

    • Google Scholar
17. 1. Elyashiv E
    2. Bullaughey K
    3. Sattath S
    4. Rinott Y
    5. Przeworski M
    6. Sella G

    (2010) Shifts in the intensity of purifying selection: an analysis of genome-wide polymorphism data from two closely related yeast species

    Genome Research 20:1558–1573.

    https://doi.org/10.1101/gr.108993.110

    • PubMed
    • Google Scholar
18. 1. Feder JL
    2. Egan SP
    3. Nosil P

    (2012) The genomics of speciation-with-gene-flow

    Trends in Genetics 28:342–350.

    https://doi.org/10.1016/j.tig.2012.03.009

    • PubMed
    • Google Scholar
19. 1. Gagnaire P-A
    2. Pavey SA
    3. Normandeau E
    4. Bernatchez L

    (2013) The genetic architecture of reproductive isolation during speciation-with-gene-flow in lake whitefish species pairs assessed by rad sequencing

    Evolution 67:2483–2497.

    https://doi.org/10.1111/evo.12075

    • Google Scholar
20. 1. Grabherr MG
    2. Haas BJ
    3. Yassour M
    4. Levin JZ
    5. Thompson DA
    6. Amit I
    7. Adiconis X
    8. Fan L
    9. Raychowdhury R
    10. Zeng Q
    11. Chen Z
    12. Mauceli E
    13. Hacohen N
    14. Gnirke A
    15. Rhind N
    16. di Palma F
    17. Birren BW
    18. Nusbaum C
    19. Lindblad-Toh K
    20. Friedman N
    21. Regev A

    (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome

    Nature Biotechnology 29:644–652.

    https://doi.org/10.1038/nbt.1883

    • PubMed
    • Google Scholar
21. 1. Graham A

    (2018) The role of land bridges, ancient environments, and migrations in the assembly of the north american flora

    Journal of Systematics and Evolution 56:405–429.

    https://doi.org/10.1111/jse.12302

    • Google Scholar
22. 1. Harris AJ
    2. Wen J
    3. Xiang Q-YJ

    (2013) Inferring the biogeographic origins of inter-continental disjunct endemics using a Bayes-DIVA approach

    Journal of Systematics and Evolution 51:117–133.

    https://doi.org/10.1111/jse.12007

    • Google Scholar
23. 1. Ian Milne R

    (2006) Northern hemisphere plant disjunctions: a window on tertiary land bridges and climate change?

    Annals of Botany 98:465–472.

    https://doi.org/10.1093/aob/mcl148

    • PubMed
    • Google Scholar
24. 1. Jordon-Thaden IE
    2. Chanderbali AS
    3. Gitzendanner MA
    4. Soltis DE

    (2015) Modified CTAB and TRIzol protocols improve RNA extraction from chemically complex embryophyta

    Applications in Plant Sciences 3:1400105.

    https://doi.org/10.3732/apps.1400105

    • Google Scholar
25. 1. Katoh K
    2. Standley DM

    (2013) MAFFT multiple sequence alignment software version 7: improvements in performance and usability

    Molecular Biology and Evolution 30:772–780.

    https://doi.org/10.1093/molbev/mst010

    • PubMed
    • Google Scholar
26. 1. Kimura M

    (1979) The neutral theory of molecular evolution

    Scientific American 241:98–126.

    https://doi.org/10.1038/scientificamerican1179-98

    • PubMed
    • Google Scholar
27. 1. Kosiol C
    2. Vinar T
    3. da Fonseca RR
    4. Hubisz MJ
    5. Bustamante CD
    6. Nielsen R
    7. Siepel A

    (2008) Patterns of positive selection in six mammalian genomes

    PLOS Genetics 4:e1000144.

    https://doi.org/10.1371/journal.pgen.1000144

    • PubMed
    • Google Scholar
28. 1. Li H-L

    (1952) Floristic relationships between eastern Asia and eastern north america

    Transactions of the American Philosophical Society 42:371–429.

    https://doi.org/10.2307/1005654

    • Google Scholar
29. 1. Li L
    2. Stoeckert CJ
    3. Roos DS

    (2003) OrthoMCL: identification of ortholog groups for eukaryotic genomes

    Genome Research 13:2178–2189.

    https://doi.org/10.1101/gr.1224503

    • PubMed
    • Google Scholar
30. 1. Li Z
    2. De La Torre AR
    3. Sterck L
    4. Cánovas FM
    5. Avila C
    6. Merino I
    7. Cabezas JA
    8. Cervera MT
    9. Ingvarsson PK
    10. Van de Peer Y

    (2017) Single-Copy genes as molecular markers for phylogenomic studies in seed plants

    Genome Biology and Evolution 9:1130–1147.

    https://doi.org/10.1093/gbe/evx070

    • PubMed
    • Google Scholar
31. 1. Lottaz C
    2. Iseli C
    3. Jongeneel CV
    4. Bucher P

    (2003) Modeling sequencing errors by combining hidden markov models

    Bioinformatics 19 Suppl 2:ii103–ii112.

    https://doi.org/10.1093/bioinformatics/btg1067

    • PubMed
    • Google Scholar
32. 1. Maddison WP
    2. Maddison DR

    (2017)

    Mesquite: a modular system for evolutionary analysis

    Mesquite: a modular system for evolutionary analysis, 3.31.

    • Google Scholar
33. 1. Magallón S
    2. Gómez-Acevedo S
    3. Sánchez-Reyes LL
    4. Hernández-Hernández T

    (2015) A metacalibrated time-tree documents the early rise of flowering plant phylogenetic diversity

    New Phytologist 207:437–453.

    https://doi.org/10.1111/nph.13264

    • PubMed
    • Google Scholar
34. 1. Manos PS
    2. Meireles JE

    (2015) Biogeographic analysis of the woody plants of the southern appalachians: implications for the origins of a regional flora

    American Journal of Botany 102:780–804.

    https://doi.org/10.3732/ajb.1400530

    • PubMed
    • Google Scholar
35. 1. Marincovich LJ
    2. Brouwers EM
    3. Hopkins DM
    4. McKenna MC

    (1990)

    The Geology of North America

    Late Mesozoic and Cenozoic paleogeographic and paleoclimatic history of the Arctic Ocean Basin, based on shallow marine faunas and terrestrial vertebrates. Pages 403–426 in A, The Geology of North America, Boulder, The Arctic Ocean region. Geological Society of America, 10.1130/DNAG-GNA-L.403.

    • Google Scholar
36. 1. Marincovich L
    2. Gladenkov AY

    (1999) Evidence for an early opening of the bering strait

    Nature 397:149–151.

    https://doi.org/10.1038/16446

    • Google Scholar
37. 1. Milne RI
    2. Abbott RJ

    (2002) The origin and evolution of tertiary relict floras

    Advances in Botanical Research 38:281–314.

    https://doi.org/10.1016/S0065-2296(02)38033-9

    • Google Scholar
38. 1. Noor MA
    2. Bennett SM

    (2009) Islands of speciation or mirages in the desert? examining the role of restricted recombination in maintaining species

    Heredity 103:439–444.

    https://doi.org/10.1038/hdy.2009.151

    • PubMed
    • Google Scholar
39. 1. Nosil P

    (2008) Speciation with gene flow could be common

    Molecular Ecology 17:2103–2106.

    https://doi.org/10.1111/j.1365-294X.2008.03715.x

    • PubMed
    • Google Scholar
40. 1. Nosil P
    2. Harmon LJ
    3. Seehausen O

    (2009) Ecological explanations for (incomplete) speciation

    Trends in Ecology & Evolution 24:145–156.

    https://doi.org/10.1016/j.tree.2008.10.011

    • PubMed
    • Google Scholar
41. 1. Nosil P
    2. Feder JL

    (2013) Genome evolution and speciation: toward quantitative descriptions of pattern and process

    Evolution 67:2461–2467.

    https://doi.org/10.1111/evo.12191

    • Google Scholar
42. 1. Pereira RJ
    2. Barreto FS
    3. Pierce NT
    4. Carneiro M
    5. Burton RS

    (2016) Transcriptome-wide patterns of divergence during allopatric evolution

    Molecular Ecology 25:1478–1493.

    https://doi.org/10.1111/mec.13579

    • PubMed
    • Google Scholar
43. 1. Posada D

    (2008) jModelTest: phylogenetic model averaging

    Molecular Biology and Evolution 25:1253–1256.

    https://doi.org/10.1093/molbev/msn083

    • PubMed
    • Google Scholar
44. 1. Qian H
    2. Jin Y
    3. Ricklefs RE

    (2017) Phylogenetic diversity anomaly in angiosperms between eastern Asia and eastern north america

    PNAS 114:11452–11457.

    https://doi.org/10.1073/pnas.1703985114

    • PubMed
    • Google Scholar
45. Software

    1. Rambaut A

    (2010)

    FigTree, version 1. 3.1

    Institute of Evolutionary Biology, University of Edinburgh, Edinburgh.
46. Software

    1. Rambaut A
    2. Drummond AJ

    (2009)

    Tracer: MCMC Trace Analysis Tool

    Oxford University Press, Oxford.
47. 1. Renaut S
    2. Grassa CJ
    3. Yeaman S
    4. Moyers BT
    5. Lai Z
    6. Kane NC
    7. Bowers JE
    8. Burke JM
    9. Rieseberg LH

    (2013) Genomic islands of divergence are not affected by geography of speciation in sunflowers

    Nature Communications 4:1827.

    https://doi.org/10.1038/ncomms2833

    • PubMed
    • Google Scholar
48. 1. Renaut S
    2. Owens GL
    3. Rieseberg LH

    (2014) Shared selective pressure and local genomic landscape lead to repeatable patterns of genomic divergence in sunflowers

    Molecular Ecology 23:311–324.

    https://doi.org/10.1111/mec.12600

    • PubMed
    • Google Scholar
49. 1. Smith TF
    2. Waterman MS

    (1981) Identification of common molecular subsequences

    Journal of Molecular Biology 147:195–197.

    https://doi.org/10.1016/0022-2836(81)90087-5

    • PubMed
    • Google Scholar
50. 1. Soria-Carrasco V
    2. Gompert Z
    3. Comeault AA
    4. Farkas TE
    5. Parchman TL
    6. Johnston JS
    7. Buerkle CA
    8. Feder JL
    9. Bast J
    10. Schwander T
    11. Egan SP
    12. Crespi BJ
    13. Nosil P

    (2014) Stick insect genomes reveal natural selection's role in parallel speciation

    Science 344:738–742.

    https://doi.org/10.1126/science.1252136

    • PubMed
    • Google Scholar
51. 1. Stadler T

    (2010) Sampling-through-time in birth-death trees

    Journal of Theoretical Biology 267:396–404.

    https://doi.org/10.1016/j.jtbi.2010.09.010

    • PubMed
    • Google Scholar
52. 1. The Angiosperm Phylogeny Group

    (2016) An update of the angiosperm phylogeny group classification for the orders and families of flowering plants: apg IV

    Botanical Journal of the Linnean Society 181:1–20.

    https://doi.org/10.1111/boj.12385

    • Google Scholar
53. 1. Tiffney BH

    (2000)

    Geographic and climatic influences on the Cretaceous and Tertiary history of Euramerican floristic similarity

    Acta Universitatis Carolinae - Geologica  44:5–16.

    • Google Scholar
54. 1. Wagner A

    (2007) Rapid detection of positive selection in genes and genomes through variation clusters

    Genetics 176:2451–2463.

    https://doi.org/10.1534/genetics.107.074732

    • PubMed
    • Google Scholar
55. 1. Wen J

    (1999) Evolution of eastern asian and eastern north american disjunct distributions in flowering plants

    Annual Review of Ecology and Systematics 30:421–455.

    https://doi.org/10.1146/annurev.ecolsys.30.1.421

    • Google Scholar
56. 1. Wen J

    (2001) Evolution of eastern Asian–Eastern North American Biogeographic Disjunctions: A Few Additional Issues

    International Journal of Plant Sciences 162:S117–S122.

    https://doi.org/10.1086/322940

    • Google Scholar
57. Book

    1. Wen JS
    2. Ickert-Bond Z-L
    3. Li R

    (2010)

    Timing and modes of evolution of eastern Asian-North American biogeographic disjunctions in seed plants

    In: Long M, Gu H, Zhou Z, editors. Darwin’s Heritage Today-Proceeding of the Darwin 200 Beijing International Conference. Beijing: Higher Education Press. pp. 252–269.

    • Google Scholar
58. 1. Wen J
    2. Nie Z-L
    3. Ickert-Bond SM

    (2016) Intercontinental disjunctions between eastern Asia and western north america in vascular plants highlight the biogeographic importance of the bering land bridge from late cretaceous to Neogene

    Journal of Systematics and Evolution 54:469–490.

    https://doi.org/10.1111/jse.12222

    • Google Scholar
59. 1. Wolf JB
    2. Ellegren H

    (2017) Making sense of genomic islands of differentiation in light of speciation

    Nature Reviews Genetics 18:87–100.

    https://doi.org/10.1038/nrg.2016.133

    • PubMed
    • Google Scholar
60. 1. Wu ZY

    (1983) On the significance of pacific intercontinental discontinuity

    Annals of the Missouri Botanical Garden 70:577–590.

    https://doi.org/10.2307/2398977

    • Google Scholar
61. 1. Wyatt R
    2. Broyles SB
    3. Hamrick JL
    4. Stoneburner A

    (1993) Systematic Relationships Within Gelsemium (Loganiaceae): Evidence from Isozymes and Cladistics

    Systematic Botany 18:345–355.

    https://doi.org/10.2307/2419408

    • Google Scholar
62. 1. Xiang Q-Y
    2. Soltis DE
    3. Soltis PS

    (1998) The eastern asian and eastern and western north american floristic disjunction: congruent phylogenetic patterns in seven diverse genera

    Molecular Phylogenetics and Evolution 10:178–190.

    https://doi.org/10.1006/mpev.1998.0524

    • PubMed
    • Google Scholar
63. 1. Xiang Q-Y
    2. Soltis DE
    3. Soltis PS
    4. Manchester SR
    5. Crawford DJ

    (2000) Timing the eastern Asian-eastern north american floristic disjunction: molecular clock corroborates paleontological estimates

    Molecular Phylogenetics and Evolution 15:462–472.

    https://doi.org/10.1006/mpev.2000.0766

    • PubMed
    • Google Scholar
64. 1. Xiang Q-Y
    2. Thomas DT
    3. Xiang QP

    (2011) Resolving and dating the phylogeny of cornales – Effects of taxon sampling, data partitions, and fossil calibrations

    Molecular Phylogenetics and Evolution 59:123–138.

    https://doi.org/10.1016/j.ympev.2011.01.016

    • Google Scholar
65. 1. Xiang Q‐Y (J)
    2. Soltis DE

    (2001) Dispersal‐vicariance analyses of intercontinental disjuncts: historical biogeographical implications for angiosperms in the northern hemisphere

    International Journal of Plant Sciences 162:S29–S39.

    https://doi.org/10.1086/323332

    • Google Scholar
66. 1. Xiang Q-Y
    2. Thomas DT

    (2008)

    Tracking character evolution and biogeographic history through time in Cornaceae - Does choice of methods matter?

    Journal of Systematics and Evolution 46:349–374.

    • Google Scholar
67. 1. Yang Z

    (1997) PAML: a program package for phylogenetic analysis by maximum likelihood

    Bioinformatics 13:555–556.

    https://doi.org/10.1093/bioinformatics/13.5.555

    • Google Scholar
68. Book

    1. Yang Z

    (2003)

    Adaptive molecular evolution

    In: Balding D, Bishop M, Cannings C, editors. Handbook of Statistical Genetics (2nd Edition). New York: Wiley. pp. 229–254.

    • Google Scholar
69. 1. Yang Z

    (2007) PAML 4: phylogenetic analysis by maximum likelihood

    Molecular Biology and Evolution 24:1586–1591.

    https://doi.org/10.1093/molbev/msm088

    • PubMed
    • Google Scholar
70. 1. Yang Z
    2. Nielsen R

    (2000) Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models

    Molecular Biology and Evolution 17:32–43.

    https://doi.org/10.1093/oxfordjournals.molbev.a026236

    • PubMed
    • Google Scholar

Author details

1. Yibo Dong

   1. Department of Plant and Microbial Biology, North Carolina State University, Raleigh, United States
   2. Plant Biology Division, Noble Research Institute, Ardmore, United States

   Contribution

   Software, Formal analysis, Investigation, Visualization, Methodology, Writing—original draft, Writing—review and editing

   Contributed equally with

   Shichao Chen

   Competing interests

   No competing interests declared

2. Shichao Chen

   1. Florida Museum of Natural History, University of Florida, Gainesville, United States
   2. Department of Biology, University of Florida, Gainesville, United States
   3. School of Life Sciences and Technology, Tongji University, Shanghai, China

   Contribution

   Formal analysis, Investigation, Methodology, Writing—original draft, Writing—review and editing

   Contributed equally with

   Yibo Dong

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-7160-4245

3. Shifeng Cheng

   Beijing Genomics Institute, Shenzhen, China

   Contribution

   Data curation, Software, Formal analysis, Investigation, Methodology, Writing—review and editing

   Competing interests

   No competing interests declared

4. Wenbin Zhou

   Department of Plant and Microbial Biology, North Carolina State University, Raleigh, United States

   Contribution

   Software, Formal analysis, Investigation, Visualization, Writing—review and editing

   Competing interests

   No competing interests declared

5. Qing Ma

   Department of Plant and Microbial Biology, North Carolina State University, Raleigh, United States

   Contribution

   Investigation, Writing—review and editing, Collection of experimental materials

   Competing interests

   No competing interests declared

6. Zhiduan Chen

   State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing, China

   Contribution

   Resources, Investigation, Writing—review and editing

   Competing interests

   No competing interests declared

7. Cheng-Xin Fu

   Laboratory of Systematic & Evolutionary Botany and Biodiversity, College of Life Sciences, Zhejiang University, Hangzhou, China

   Contribution

   Resources, Investigation, Project administration, Writing—review and editing

   Competing interests

   No competing interests declared

8. Xin Liu

   Beijing Genomics Institute, Shenzhen, China

   Contribution

   Conceptualization, Resources, Funding acquisition, Investigation, Writing—review and editing

   For correspondence

   liuxin@genomics.cn

   Competing interests

   No competing interests declared

9. Yun-peng Zhao

   Laboratory of Systematic & Evolutionary Botany and Biodiversity, College of Life Sciences, Zhejiang University, Hangzhou, China

   Contribution

   Conceptualization, Resources, Funding acquisition, Investigation, Project administration, Writing—review and editing

   For correspondence

   ypzhao913@gmail.com

   Competing interests

   No competing interests declared

10. Pamela S Soltis

    Florida Museum of Natural History, University of Florida, Gainesville, United States

    Contribution

    Conceptualization, Resources, Supervision, Funding acquisition, Investigation, Methodology, Project administration, Writing—review and editing

    For correspondence

    psoltis@flmnh.ufl.edu

    Competing interests

    No competing interests declared

11. Gane Ka-Shu Wong

    1. Beijing Genomics Institute, Shenzhen, China
    2. Department of Biological Sciences, University of Alberta, Edmonton, Canada
    3. Department of Medicine, University of Alberta, Edmonton, Canada

    Contribution

    Conceptualization, Funding acquisition, Investigation, Writing—review and editing

    For correspondence

    gane@ualberta.ca

    Competing interests

    No competing interests declared

12. Douglas E Soltis

    1. Florida Museum of Natural History, University of Florida, Gainesville, United States
    2. Department of Biology, University of Florida, Gainesville, United States

    Contribution

    Conceptualization, Resources, Data curation, Supervision, Funding acquisition, Investigation, Methodology, Project administration, Writing—review and editing

    For correspondence

    dsoltis@ufl.edu

    Competing interests

    No competing interests declared

13. Jenny Xiang

    Department of Plant and Microbial Biology, North Carolina State University, Raleigh, United States

    Contribution

    Conceptualization, Resources, Data curation, Supervision, Funding acquisition, Investigation, Methodology, Writing—original draft, Project administration, Writing—review and editing

    For correspondence

    jenny_xiang@ncsu.edu

    Competing interests

    No competing interests declared

    "This ORCID iD identifies the author of this article:" 0000-0002-9016-0678

• 4,281

  views

• 513

  downloads

• 20

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.