Historical Biogeography of Melanthiaceae: A Case of Out-of-North America Through the Bering Land Bridge

First Record of Comparative Plastid Genome Analysis and Phylogenetic Relationships among Corylopsis Siebold & Zucc. (Hamamelidaceae)

Corylopsis Siebold & Zucc. (Hamamelidaceae) is widely used as a horticultural plant and comprises approximately 25 species in East Asia. Molecular research is essential to distinguish Corylopsis species, which are morphologically similar. Molecular research has been conducted using a small number of genes but not in Corylopsis. Plastid genomes of Corylopsis species (Corylopsis gotoana, Corylopsis pauciflora, and Corylopsis sinensis) were sequenced using next-generation sequencing techniques. Repeats and nucleotide diversity that could be used as DNA markers were also investigated. A phylogenetic investigation was carried out using 79 protein-coding genes to infer the evolutionary relationships within the genus Corylopsis. By including new plastomes, the overall plastid genome structure of Corylopsis was similar. Simple sequence repeats of 73-106 SSRs were identified in the protein-coding genes of the plastid genomes, and 33-40 long repeat sequences were identified in the plastomes. The Pi value of the rpl33_rps18 region, an intergenic spacer, was the highest. Phylogenetic analysis demonstrated that Corylopsis is a monophyletic group and Loropetalum is closely related to Corylopsis. C. pauciflora, C. gotoana, and C. spicata formed a clade distributed in Japan, whereas C. sinensis, C. glandulifera, and C. velutina formed a clade that was distributed in China.

East Asian-North American disjunctions and phylogenetic relationships within subtribe Nepetinae (Lamiaceae)

Biogeographic disjunctions, including intercontinental disjunctions, are frequent across plant lineages and have been of considerable interest to biologists for centuries. Their study has been reinvigorated by molecular dating and associated comparative methods. One of the "classic" disjunction patterns is that between Eastern Asia and North America. It has been speculated that this pattern is the result of vicariance following the sundering of a widespread Acrto-Teritary flora. Subtribe Nepetinae in the mint family (Lamiaceae) is noteworthy because it contains three genera with this disjunction pattern: Agastache, Dracocephalum, and Meehania. These disjunctions are ostensibly the result of three separate events, allowing for concurrent testing of the tempo, origin, and type of each biogeographic event. Using four plastid and four nuclear markers, we estimated divergence times and analyzed the historical biogeography of Nepetinae, including comprehensive sampling of all major clades for the first time. We recover a well-supported and largely congruent phylogeny of Nepetinae between genomic compartments, although several cases of cyto-nuclear discordance are evident. We demonstrate that the three disjunctions are pseudo-congruent, with unidirectional movement from East Asia at slightly staggered times during the late Miocene and early Pliocene. With the possible exception of Meehania, we find that vicariance is likely the underlying driver of these disjunctions. The biogeographic history of Meehania in North America may be best explained by long-distance dispersal, but a more complete picture awaits deeper sampling of the nuclear genome and more advanced biogeographical models.

Holocentromeres can consist of merely a few megabase-sized satellite arrays

The centromere is the chromosome region where microtubules attach during cell division. In contrast to monocentric chromosomes with one centromere, holocentric species usually distribute hundreds of centromere units along the entire chromatid. We assembled the chromosome-scale reference genome and analyzed the holocentromere and (epi)genome organization of the lilioid Chionographis japonica. Remarkably, each of its holocentric chromatids consists of only 7 to 11 evenly spaced megabase-sized centromere-specific histone H3-positive units. These units contain satellite arrays of 23 and 28 bp-long monomers capable of forming palindromic structures. Like monocentric species, C. japonica forms clustered centromeres in chromocenters at interphase. In addition, the large-scale eu- and heterochromatin arrangement differs between C. japonica and other known holocentric species. Finally, using polymer simulations, we model the formation of prometaphase line-like holocentromeres from interphase centromere clusters. Our findings broaden the knowledge about centromere diversity, showing that holocentricity is not restricted to species with numerous and small centromere units.

Out of the Qinghai-Tibetan plateau: Origin, evolution and historical biogeography of Morchella (both Elata and Esculenta clades)

Introduction

Morchella has become a research hotspot because of its wide distribution, delicious taste, and phenotypic plasticity. The Qinghai-Tibet Plateau subkingdoms (QTPs) are known as the cradle of Ice age biodiversity. However, the diversity of Morchella in the QTPs has been poorly investigated, especially in phylogenetic diversity, origin, and biogeography.

Methods

The genealogical concordance phylogenetic species recognition (GCPSR, based on Bayesian evolutionary analysis using sequences from the internal transcribed spacer (ITS), nuclear large subunit rDNA (nrLSU), translation elongation factor 1-α (EF1-α), and the largest and second largest subunits of RNA polymerase II (RPB1 and RPB2)), differentiation time estimation, and ancestral region reconstruction were used to infer Morchella's phylogenetic relationships and historical biogeography in the QTPs.

Results

Firstly, a total of 18 Morchella phylogenetic species are recognized in the QTPs, including 10 Elata clades and 8 Esculenta clades of 216 individuals Secondly, the divergences of the 18 phylogenetic species were 50.24-4.20 Mya (Eocene-Pliocene), which was closely related to the geological activities in the QTPs. Furthermore, the ancestor of Morchella probably originated in the Northern regions (Qilian Shan, Elata cade) and southwestern regions (Shangri-La, Esculenta clade) of QTPs and might have migrated from North America (Rufobrunnea clade) via Beringian Land Bridge (BLB) and Long-Distance Dispersal (LDD) expansions during the Late Cretaceous. Moreover, as the cradle of species origin and diversity, the fungi species in the QTPs have spread out and diffused to Eurasia and South Africa starting in the Paleogene Period.

Conclusion

This is the first report that Esculenta and Elata clade of Morchella originated from the QTPs because of orogenic, and rapid differentiation of fungi is strongly linked to geological uplift movement and refuge in marginal areas of the QTPs. Our findings contribute to increasing the diversity of Morchella and offer more evidence for the origin theory of the QTPs.

Molecular Phylogeny and Historical Biogeography of Goodyera R. Br. (Orchidaceae): A Case of the Vicariance Between East Asia and North America

Understanding of intercontinental distribution in the Northern Hemisphere has attracted a lot of attention from botanists. However, although Orchidaceae is the largest group of angiosperms, biogeographical studies on the disjunctive pattern have not been sufficient for this family. Goodyera R. Br. (tribe Cranichideae, subfamily Orchidoideae, family Orchidaceae) is widely distributed in temperate and tropical regions. Although the phylogenetic relationship of Goodyera inferred from both morphological and molecular data has been conducted, the sampled taxa were mainly distributed in Asia regions that resulted in non-monophyly of this genus. In this study, the complete plastid genomes of Goodyera, generated by next-generation sequencing (NGS) technique and sampled in East Asia and North America, were used to reconstruct phylogeny and explore the historical biogeography. A total of 18 Goodyera species including seven newly sequenced species were analyzed. Based on 79 protein-coding genes, the phylogenetic analysis revealed that Goodyera could be subdivided into four subclades with high support values. The polyphyletic relationships among Goodyera taxa were confirmed, and the unclear position of G. foliosa was also resolved. The datasets that are composed of the 14 coding sequences (CDS) (matK, atpF, ndhK, accD, cemA, clpP, rpoA, rpl22, ndhF, ccsA, ndhD, ndhI, ndhA, and ycf 1) showed the same topology derived from 79 protein-coding genes. Molecular dating analyses revealed the origin of Goodyera in the mid-Miocene (15.75 Mya). Nearctic clade of Goodyera was diverged at 10.88 Mya from their most recent common ancestor (MRCA). The biogeographical reconstruction suggests that subtropical or tropical Asia is the origin of Goodyera and it has subsequently spread to temperate Asia during the Miocene. In addition, Nearctic clade is derived from East Asian species through Bering Land Bridge (BLB) during the Miocene. The speciation of Goodyera is most likely to have occurred during Miocene, and climatic and geological changes are thought to have had a part in this diversification. Our findings propose both origin and vicariance events of Goodyera for the first time and add an example for the biogeographical history of the Northern Hemisphere.

Evolution of floral characters and biogeography of Heloniadeae (Melanthiaceae): an example of breeding system shifts with inflorescence change

Heloniadeae (Melanthiaceae) presents an East Asia–North America disjunct distribution. Different molecular and morphological data nevertheless support the tribe as a monophyletic group. However, their phylogenetic relationships and biogeographic history, together with the character evolution, are not clear. Therefore, we constructed a Bayesian phylogenetic tree for Heloniadeae using cpDNA and inferred the historical biogeography and floral character evolution. The results revealed that Heloniadeae was distributed in high-latitudes of East Asia and North America, originating since 22.2 mya. The East Asia clade migrated into southwest China, and subsequently colonized the Korean Peninsula, Taiwan, the Ryukyus, and spread northward to Japan and southern Sakhalin. The evolution of the inflorescence and number of flowers were phylogenetically conserved, associated with the historical biogeography of Heloniadeae. The inflorescences transferred from raceme to sub-umbel, and the number of flowers decreased during the dispersal process, which may be accompanied by changes in the breeding system. Besides, the anthesis period was more affected by the habitat environment than phylogenetic constraints. The flowering temperature of Heloniadeae was below 20 °C in most species, except H. kawanoi. Such a low temperature might not be conductive to pollinator activities, but it could be compensated by sustaining seed production with long-lasting flowers.

Genome Size Doubling Arises From the Differential Repetitive DNA Dynamics in the Genus Heloniopsis (Melanthiaceae)

Plant genomes are highly diverse in size and repetitive DNA composition. In the absence of polyploidy, the dynamics of repetitive elements, which make up the bulk of the genome in many species, are the main drivers underpinning changes in genome size and the overall evolution of the genomic landscape. The advent of high-throughput sequencing technologies has enabled investigation of genome evolutionary dynamics beyond model plants to provide exciting new insights in species across the biodiversity of life. Here we analyze the evolution of repetitive DNA in two closely related species of Heloniopsis (Melanthiaceae), which despite having the same chromosome number differ nearly twofold in genome size [i.e., H. umbellata (1C = 4,680 Mb), and H. koreana (1C = 2,480 Mb)]. Low-coverage genome skimming and the RepeatExplorer2 pipeline were used to identify the main repeat families responsible for the significant differences in genome sizes. Patterns of repeat evolution were found to correlate with genome size with the main classes of transposable elements identified being twice as abundant in the larger genome of H. umbellata compared with H. koreana. In addition, among the satellite DNA families recovered, a single shared satellite (HeloSAT) was shown to have contributed significantly to the genome expansion of H. umbellata. Evolutionary changes in repetitive DNA composition and genome size indicate that the differences in genome size between these species have been underpinned by the activity of several distinct repeat lineages.

Insights into phylogenetic relationships and genome evolution of subfamily Commelinoideae (Commelinaceae Mirb.) inferred from complete chloroplast genomes

Background

Commelinaceae (Commelinales) comprise 41 genera and are widely distributed in both the Old and New Worlds, except in Europe. The relationships among genera in this family have been suggested in several morphological and molecular studies. However, it is difficult to explain their relationships due to high morphological variations and low support values. Currently, many researchers have been using complete chloroplast genome data for inferring the evolution of land plants. In this study, we completed 15 new plastid genome sequences of subfamily Commelinoideae using the Mi-seq platform. We utilized genome data to reveal the structural variations and reconstruct the problematic positions of genera for the first time.

Results

All examined species of Commelinoideae have three pseudogenes (accD, rpoA, and ycf15), and the former two might be a synapomorphy within Commelinales. Only four species in tribe Commelineae presented IR expansion, which affected duplication of the rpl22 gene. We identified inversions that range from approximately 3 to 15 kb in four taxa (Amischotolype, Belosynapsis, Murdannia, and Streptolirion). The phylogenetic analysis using 77 chloroplast protein-coding genes with maximum parsimony, maximum likelihood, and Bayesian inference suggests that Palisota is most closely related to tribe Commelineae, supported by high support values. This result differs significantly from the current classification of Commelinaceae. Also, we resolved the unclear position of Streptoliriinae and the monophyly of Dichorisandrinae. Among the ten CDS (ndhH, rpoC2, ndhA, rps3, ndhG, ndhD, ccsA, ndhF, matK, and ycf1), which have high nucleotide diversity values (Pi > 0.045) and over 500 bp length, four CDS (ndhH, rpoC2, matK, and ycf1) show that they are congruent with the topology derived from 77 chloroplast protein-coding genes.

Conclusions

In this study, we provide detailed information on the 15 complete plastid genomes of Commelinoideae taxa. We identified characteristic pseudogenes and nucleotide diversity, which can be used to infer the family evolutionary history. Also, further research is needed to revise the position of Palisota in the current classification of Commelinaceae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07541-1.

Complete chloroplast genomes shed light on phylogenetic relationships, divergence time, and biogeography of Allioideae (Amaryllidaceae)

Allioideae includes economically important bulb crops such as garlic, onion, leeks, and some ornamental plants in Amaryllidaceae. Here, we reported the complete chloroplast genome (cpDNA) sequences of 17 species of Allioideae, five of Amaryllidoideae, and one of Agapanthoideae. These cpDNA sequences represent 80 protein-coding, 30 tRNA, and four rRNA genes, and range from 151,808 to 159,998 bp in length. Loss and pseudogenization of multiple genes (i.e., rps2, infA, and rpl22) appear to have occurred multiple times during the evolution of Alloideae. Additionally, eight mutation hotspots, including rps15-ycf1, rps16-trnQ-UUG, petG-trnW-CCA, psbA upstream, rpl32-trnL-UAG, ycf1, rpl22, matK, and ndhF, were identified in the studied Allium species. Additionally, we present the first phylogenomic analysis among the four tribes of Allioideae based on 74 cpDNA coding regions of 21 species of Allioideae, five species of Amaryllidoideae, one species of Agapanthoideae, and five species representing selected members of Asparagales. Our molecular phylogenomic results strongly support the monophyly of Allioideae, which is sister to Amaryllioideae. Within Allioideae, Tulbaghieae was sister to Gilliesieae-Leucocoryneae whereas Allieae was sister to the clade of Tulbaghieae- Gilliesieae-Leucocoryneae. Molecular dating analyses revealed the crown age of Allioideae in the Eocene (40.1 mya) followed by differentiation of Allieae in the early Miocene (21.3 mya). The split of Gilliesieae from Leucocoryneae was estimated at 16.5 mya. Biogeographic reconstruction suggests an African origin for Allioideae and subsequent spread to Eurasia during the middle Eocene. Cool and arid conditions during the late Eocene led to isolation between African and Eurasian species. African Allioideae may have diverged to South American taxa in the late Oligocene. Rather than vicariance, long-distance dispersal is the most likely explanation for intercontinental distribution of African and South American Allioideae species.

Chromosomal Evolution and Apomixis in the Cruciferous Tribe Boechereae

The mustard family (Brassicaceae) comprises several dozen monophyletic clades usually ranked as tribes. The tribe Boechereae plays a prominent role in plant research due to the incidence of apomixis and its close relationship to Arabidopsis. This tribe, largely confined to western North America, harbors nine genera and c. 130 species, with >90% of species belonging to the genus Boechera. Hundreds of apomictic diploid and triploid Boechera hybrids have spurred interest in this genus, but the remaining Boechereae genomes remain virtually unstudied. Here we report on comparative genome structure of six genera (Borodinia, Cusickiella, Phoenicaulis, Polyctenium, Nevada, and Sandbergia) and three Boechera species as revealed by comparative chromosome painting (CCP). All analyzed taxa shared the same seven-chromosome genome structure. Comparisons with the sister Halimolobeae tribe (n = 8) showed that the ancestral Boechereae genome (n = 7) was derived from an older n = 8 genome by descending dysploidy followed by the divergence of extant Boechereae taxa. As tribal divergence post-dated the origin of four tribe-specific chromosomes, it is proposed that these chromosomal rearrangements were a key evolutionary innovation underlaying the origin and diversification of the Boechereae in North America. Although most Boechereae genera exhibit genomic conservatism, intra-tribal cladogenesis has occasionally been accompanied by chromosomal rearrangements (particularly inversions). Recently, apomixis was reported in the Boechereae genera Borodinia and Phoenicaulis. Here, we report sexual reproduction in diploid Nevada, diploid Sandbergia, and tetraploid Cusickiella and aposporous apomixis in tetraploids of Polyctenium and Sandbergia. In sum, apomixis is now known to occur in five of the nine Boechereae genera.

Implications of plastome evolution in the true lilies (monocot order Liliales)

The families of the monocot order Liliales exhibit highly contrasting characteristic of photosynthetic and mycoheterotrophic life histories. Although previous phylogenetic and morphological studies of Liliales have been conducted, they have not examined molecular evolution associated with this contrasting phenomenon. Here, we conduct the first comparative plastome study of all ten families of Liliales using 29 newly sequenced plastid genomes analyzed together with previously published data. We also present a phylogenetic analysis for Liliales of 78 plastid genes combined with 22 genes from all three genomes (nuclear 18S rDNA and phyC; 17 plastid genes; and mitochondrial matR, atpA, and cob). Within the newly generated phylogenetic tree of Liliales, we evaluate the ancestral state changes of selected morphological traits in the order. There are no significant differences in plastid genome features among species that show divergent characteristics correlated with family circumscriptions. However, the results clearly differentiate between photosynthetic and mycoheterotrophic taxa of Liliales in terms of genome structure, and gene content and order. The newly sequenced plastid genomes and combined three-genome data revealed Smilacaceae as sister to Liliaceae instead of Philesiaceae and Ripogonaceae. Additionally, we propose a revised familial classification system of Liliales that consists of nine families, considering Ripogonaceae a synonym of Philesiaceae. The ancestral state reconstruction indicated synapomorphies for each family of Liliales, except Liliaceae, Melanthiaceae and Colchicaceae. A taxonomic key for all nine families of Liliales is also provided.

Plastome phylogenomics, biogeography, and clade diversification of Paris (Melanthiaceae)

BACKGROUND

Paris (Melanthiaceae) is an economically important but taxonomically difficult genus, which is unique in angiosperms because some species have extremely large nuclear genomes. Phylogenetic relationships within Paris have long been controversial. Based on complete plastomes and nuclear ribosomal DNA (nrDNA) sequences, this study aims to reconstruct a robust phylogenetic tree and explore historical biogeography and clade diversification in the genus.

RESULTS

All 29 species currently recognized in Paris were sampled. Whole plastomes and nrDNA sequences were generated by the genome skimming approach. Phylogenetic relationships were reconstructed using the maximum likelihood and Bayesian inference methods. Based on the phylogenetic framework and molecular dating, biogeographic scenarios and historical diversification of Paris were explored. Significant conflicts between plastid and nuclear datasets were identified, and the plastome tree is highly congruent with past interpretations of the morphology. Ancestral area reconstruction indicated that Paris may have originated in northeastern Asia and northern China, and has experienced multiple dispersal and vicariance events during its diversification. The rate of clade diversification has sharply accelerated since the Miocene/Pliocene boundary.

CONCLUSIONS

Our results provide important insights for clarifying some of the long-standing taxonomic debates in Paris. Cytonuclear discordance may have been caused by ancient and recent hybridizations in the genus. The climatic and geological changes since the late Miocene, such as the intensification of Asian monsoon and the rapid uplift of Qinghai-Tibet Plateau, as well as the climatic fluctuations during the Pleistocene, played essential roles in driving range expansion and radiative diversification in Paris. Our findings challenge the theoretical prediction that large genome sizes may limit speciation.

A fossil-calibrated phylogeny reveals the biogeographic history of the Cladrastis clade, an amphi-Pacific early-branching group in papilionoid legumes

The early-branching Cladrastis clade of papilionoid legumes (Leguminosae, Papilionoideae) has an intriguing amphi-Pacific disjunct distribution in eastern Asia and temperate-tropical Americas. Here we used nuclear and three plastid regions to reconstruct the phylogenetic relationships and divergence times in the Cladrastis clade, as well as the evolution of morphological characters that might have been key in its biogeographic history. The ancestral character state estimation revealed that the most recent common ancestor of the Cladrastis clade was deciduous trees possessing compressed, winged fruits. The Cladrastis clade was inferred to have originated in the mid-latitude thermophilic forests of North America in the early Eocene, followed by the split between ancestors of wing-fruited Platyosprion and the non-wing-fruited group, and later the divergence of Cladrastis s.s. from the non-wing-fruited group in middle Eocene. Platyosprion and Cladrastis s.s. display an "out-of-North-America" biogeographic pattern and might have migrated to Asia via the Bering land bridge (BLB) or the North Atlantic land bridges (NALB) during middle to late Eocene. Our results, coupled with the relatively well documented fossil record for the clade, suggest that Platyosprion experienced an extinction event in North America caused by climatic cooling around the Eocene-Oligocene transition, which drove a major vegetation shift in western North America, in turn serving as a barrier for the vicariance of Pickeringia and Styphnolobium. The evolution of shrubby habit and sclerophyllous leaves in the former might be adaption to the chaparral vegetation in southwestern North America; the latter gained the trait of moniliform, succulent fruit. Styphnolobium further dispersed southward to tropical North America in the Oligocene, and eastward to Asia through BLB during middle Miocene. Subsequent sundering of BLB facilitated the vicariance of St. affine and St. japonicum.