The first homosporous lycophyte genome revealed the association between the recent dynamic accumulation of LTR-RTs and genome size variation

The contrasting genome size between homosporous and heterosporous plants is fascinating. Different from the heterosporous seed plants and mainly homosporous ferns, the lycophytes are either heterosporous (Isoetales and Selaginellales) or homosporous (Lycopodiales). Many lycophytes are the resource plants of Huperzine A (HupA) which is invaluable for treating Alzheimer's disease. For the seed-free vascular plants, several high-quality genomes of heterosporous Selaginella, homosporous ferns (maidenhair fern, monkey spider tree fern), and heterosporous ferns (Azolla) have been published and provided important insights into the origin and evolution of early land plants. However, the homosporous lycophyte genome has not been decoded. Here, we assembled the first homosporous lycophyte genome and conducted comparative genomic analyses by applying a reformed pipeline for filtering out non-plant sequences. The obtained genome size of Lycopodium clavatum is 2.30 Gb, distinguished in more than 85% repetitive elements of which 62% is long terminal repeat (LTR). This study disclosed a high birth rate and a low death rate of the LTR-RTs in homosporous lycophytes, but the opposite occurs in heterosporous lycophytes. we propose that the recent activity of LTR-RT is responsible for the immense genome size variation between homosporous and heterosporous lycophytes. By combing Ks analysis with a phylogenetic approach, we discovered two whole genome duplications (WGD). Morover, we identified all the five recognized key enzymes for the HupA biosynthetic pathway in the L. clavatum genome, but found this pathway incomplete in other major lineages of land plants. Overall, this study is of great importance for the medicinal utilization of lycophytes and the decoded genome data will be a key cornerstone to elucidate the evolution and biology of early vascular land plants.

Mitogenome-based phylogenomics provides insights into the positions of the enigmatic sinensis group and the sanguinolenta group in Selaginellaceae (Lycophyte)

Spikemoss (Selaginellaceae) is one of the basal lineages of vascular plants. This family has a single genus Selaginella which consists of about 750 extant species. The phylogeny of Selaginellaceae has been extensively studied mainly based on plastid DNA and a few nuclear sequences. However, the placement of the enigmatic sinensis group is a long-term controversy because of the long branch in the plastid DNA phylogeny. The sanguinolenta group is also a phylogenetically problematic clade owing to two alternative positions resulted from different datasets. Here, we newly sequenced 34 mitochondrial genomes (mitogenomes) of individuals representing all seven subgenera and major clades in Selaginellaceae. We assembled the draft mitogenomes and annotated the genes and performed phylogenetic analyses based on the shared 17 mitochondrial genes. Our major results include: (1) all the assembled mitogenomes have complicated structures, unparalleled high GC content and a small gene content set, and the positive correlations among GC content, substitution rates and the number of RNA editing sites hold; (2) the sinensis group was well supported as a member of subg. Stachygynandrum; (3) the sanguinolenta group was strongly resolved as sister to all other Selaginella species except for subg. Selaginella. This study demonstrates the potential of mitogenome data in providing novel insights into phylogenetically recalcitrant problems.

The evolution of extremely diverged plastomes in Selaginellaceae (lycophyte) is driven by repeat patterns and the underlying DNA maintenance machinery

Two factors are proposed to account for the unusual features of organellar genomes: the disruptions of organelle-targeted DNA replication, repair, and recombination (DNA-RRR) systems in the nuclear genome and repetitive elements in organellar genomes. Little is known about how these factors affect organellar genome evolution. The deep-branching vascular plant family Selaginellaceae is known to have a deficient DNA-RRR system and convergently evolved organellar genomes. However, we found that the plastid genome (plastome) of Selaginella sinensis has extremely accelerated substitution rates, a low GC content, pervasive repeat elements, a dynamic network structure, and it lacks direct or inverted repeats. Unexpectedly, its organelle DNA-RRR system is short of a plastid-targeted Recombinase A1 (RecA1) and a mitochondrion-targeted RecA3, in line with other explored Selaginella species. The plastome contains a large collection of short- and medium-sized repeats. Given the absence of RecA1 surveillance, we propose that these repeats trigger illegitimate recombination, accelerated mutation rates, and structural instability. The correlations between repeat quantity and architectural complexity in the Selaginella plastomes support these conclusions. We, therefore, hypothesize that the interplay of the deficient DNA-RRR system and the high repeat content has led to the extraordinary divergence of the S. sinensis plastome. Our study not only sheds new light on the mechanism of plastome divergence by emphasizing the power of cytonuclear integration, but it also reconciles the longstanding contradiction on the effects of DNA-RRR system disruption on genome structure evolution.

Fire-prone Rhamnaceae with South African affinities in Cretaceous Myanmar amber

The rapid Cretaceous diversification of flowering plants remains Darwin's 'abominable mystery' despite numerous fossil flowers discovered in recent years. Wildfires were frequent in the Cretaceous and many such early flower fossils are represented by charcoalified fragments, lacking complete delicate structures and surface textures, making their similarity to living forms difficult to discern. Furthermore, scarcity of information about the ecology of early angiosperms makes it difficult to test hypotheses about the drivers of their diversification, including the role of fire in shaping flowering plant evolution. We report the discovery of two exquisitely preserved fossil flower species, one identical to the inflorescences of the extant crown-eudicot genus Phylica and the other recovered as a sister group to Phylica, both preserved as inclusions together with burned plant remains in Cretaceous amber from northern Myanmar (~99 million years ago). These specialized flower species, named Phylica piloburmensis sp. nov. and Eophylica priscastellata gen. et sp. nov., exhibit traits identical to those of modern taxa in fire-prone ecosystems such as the fynbos of South Africa, and provide evidence of fire adaptation in angiosperms.

Ellipinema and ×Ellipisorus? Just Lepisorus (Polypodiaceae)!

The establishment of a segregate lepisoroid fern genus Ellipinema was mainly to accommodate the isolated position of Lepisorus jakonensis (Polypodiaceae) recovered in plastid gene tree. Using newly obtained nuclear data, we recovered that Ellipinema and allied genera, such as Lepidomicrosorium, Lemmaphyllum, Neolepisorus, Paragramma, Tricholepidium and Weatherbya are deeply nested within Lepisorus. The nuclear phylogeny showing incongruent phylogenetic placement in comparison with plastid results perhaps indicated ancient hybridization events. The diagnostic morphology characterizing Ellipinema - elliptic scale-like paraphyses, which is shared by all the taxa of sect. Lepisorus and sect. Hymenophyton - falls within the range of continuous variation in the type species Ellipinema jakonense (=Lepisorus jakonensis). Our study, which integrated molecular and morphological data, demonstrates that the segregation of Ellipinema and ×Ellipisorus (= ×Lepinema Li Bing Zhang & Liang Zhang, nom. illeg.) from Lepisorus should be rejected.

Distinctive evolutionary pattern of organelle genomes linked to the nuclear genome in Selaginellaceae

Plastids and mitochondria are endosymbiotic organelles that store genetic information. The genomes of these organelles generally exhibit contrasting patterns regarding genome architecture and genetic content. However, they have similar genetic features in Selaginellaceae, and little is known about what causes parallel evolution. Here, we document the multipartite plastid genomes (plastomes) and the highly divergent mitochondrial genomes (mitogenomes) from spikemoss obtained by combining short- and long-reads. The 188-kb multipartite plastome has three ribosomal operon copies in the master genomic conformation, creating the alternative subgenomic conformation composed of 110- and 78-kb subgenomes. The long-read data indicated that the two different genomic conformations were present in almost equal proportions in the plastomes of Selaginella nipponica. The mitogenome of S. nipponica was assembled into 27 contigs with a total size of 110 kb. All contigs contained directly arranged repeats at both ends, which introduced multiple conformations. Our results showed that plastomes and mitogenomes share high tRNA losses, GC-biased nucleotides, elevated substitution rates and complicated organization. The exploration of nuclear-encoded organelle DNA replication, recombination and repair proteins indicated that, several single-targeted proteins, particularly plastid-targeted recombinase A1, have been lost in Selaginellaceae; conversely, the dual-targeted proteins remain intact. According to the reported function of recombinase A1, we propose that the plastomes of spikemoss often fail to pair homologous sequences during recombination, and the dual-targeted proteins play a key role in the convergent genetic features of plastomes and mitogenomes. Our results provide a distinctive evolutionary pattern of the organelle genomes in Selaginellaceae and evidence of their convergent evolution.

Directed Repeats Co-occur with Few Short-Dispersed Repeats in Plastid Genome of a Spikemoss, Selaginella vardei (Selaginellaceae, Lycopodiopsida)

BACKGROUND

It is hypothesized that the highly conserved inverted repeats (IR) structure of land plant plastid genomes (plastomes) is beneficial for stabilizing plastome organization, whereas the mechanism of the occurrence and stability maintenance of the recently reported direct repeats (DR) structure is yet awaiting further exploration. Here we describe the DR structure of the Selaginella vardei (Selaginellaceae) plastome, to elucidate the mechanism of DR occurrence and stability maintenance.

RESULTS

The plastome of S. vardei is 121,254 bp in length and encodes 76 genes, of which 62 encode proteins, 10 encode tRNAs, and four encode rRNAs. Unexpectedly, the two identical rRNA gene regions (13,893 bp) are arranged in a direct orientation (DR), rather than inverted. Comparing to the IR organization in Isoetes flaccida (Isoetaceae, Lycopodiopsida) plastome, a ca. 50-kb trnN-trnF inversion that spans one DR copy was found in the plastome of S. vardei, which might cause the orientation change. In addition, we find extremely rare short dispersed repeats (SDRs) in the plastomes of S. vardei and its closely related species S. indica.

CONCLUSIONS

We suggest that the ca. 50-kb inversion resulted in the DR structure, and the reduction in SDRs plays a key role in maintaining the stability of plastomes with DR structure by avoiding potential secondary recombination. We further confirmed the presence of homologous recombination between DR regions, which are able to generate subgenomes and form diverse multimers. Our study deepens the understanding of Selaginella plastomes and provides new insights into the diverse plastome structures in land plants.

The Unique Evolutionary Trajectory and Dynamic Conformations of DR and IR/DR-Coexisting Plastomes of the Early Vascular Plant Selaginellaceae (Lycophyte)

Both direct repeats (DR) and inverted repeats (IR) are documented in the published plastomes of Selaginella species indicating the unusual and diverse plastome structure in the family Selaginellaceae. In this study, we newly sequenced complete plastomes of seven species from five main lineages of Selaginellaceae and also resequenced three species (Selaginella tamariscina, Selaginella uncinata, and Selaginella moellendorffii) to explore the evolutionary trajectory of Selaginellaceae plastomes. Our results showed that the plastomes of Selaginellaceae vary remarkably in size, gene contents, gene order, and GC contents. Notably, both DR and IR structures existed in the plastomes of Selaginellaceae with DR structure being an ancestral state. The occurrence of DR structure was at ∼257 Ma and remained in most subgenera of Selaginellaceae, whereas IR structure only reoccurred in Selaginella sect. Lepidophyllae (∼143 Ma) and Selaginella subg. Heterostachys (∼19 Ma). The presence of a pair of large repeats psbK-trnQ, together with DR/IR region in Selaginella bisulcata, Selaginella pennata, S. uncinata, and Selaginella hainanensis, could frequently mediate diverse homologous recombination and create approximately equal stoichiometric isomers (IR/DR-coexisting) and subgenomes. High proportion of repeats is presumably responsible for the dynamic IR/DR-coexisting plastomes, which possess a lower synonymous substitution rate (dS) compared with DR-possessing and IR-possessing plastomes. We propose that the occurrence of DR structure, together with few repeats, is possibly selected to keep the stability of plastomes and the IR/DR-coexisting plastomes also reached an equilibrium in plastome organization through highly efficient homologous recombination to maintain stability.

Plastid Phylogenomics Resolve Deep Relationships among Eupolypod II Ferns with Rapid Radiation and Rate Heterogeneity

The eupolypods II ferns represent a classic case of evolutionary radiation and, simultaneously, exhibit high substitution rate heterogeneity. These factors have been proposed to contribute to the contentious resolutions among clades within this fern group in multilocus phylogenetic studies. We investigated the deep phylogenetic relationships of eupolypod II ferns by sampling all major families and using 40 plastid genomes, or plastomes, of which 33 were newly sequenced with next-generation sequencing technology. We performed model-based analyses to evaluate the diversity of molecular evolutionary rates for these ferns. Our plastome data, with more than 26,000 informative characters, yielded good resolution for deep relationships within eupolypods II and unambiguously clarified the position of Rhachidosoraceae and the monophyly of Athyriaceae. Results of rate heterogeneity analysis revealed approximately 33 significant rate shifts in eupolypod II ferns, with the most heterogeneous rates (both accelerations and decelerations) occurring in two phylogenetically difficult lineages, that is, the Rhachidosoraceae–Aspleniaceae and Athyriaceae clades. These observations support the hypothesis that rate heterogeneity has previously constrained the deep phylogenetic resolution in eupolypods II. According to the plastome data, we propose that 14 chloroplast markers are particularly phylogenetically informative for eupolypods II both at the familial and generic levels. Our study demonstrates the power of a character-rich plastome data set and high-throughput sequencing for resolving the recalcitrant lineages, which have undergone rapid evolutionary radiation and dramatic changes in substitution rates.

Plastid Phylogenomics Resolve Deep Relationships among Eupolypod II Ferns with Rapid Radiation and Rate Heterogeneity

The eupolypods II ferns represent a classic case of evolutionary radiation and, simultaneously, exhibit high substitution rate heterogeneity. These factors have been proposed to contribute to the contentious resolutions among clades within this fern group in multilocus phylogenetic studies. We investigated the deep phylogenetic relationships of eupolypod II ferns by sampling all major families and using 40 plastid genomes, or plastomes, of which 33 were newly sequenced with next-generation sequencing technology. We performed model-based analyses to evaluate the diversity of molecular evolutionary rates for these ferns. Our plastome data, with more than 26,000 informative characters, yielded good resolution for deep relationships within eupolypods II and unambiguously clarified the position of Rhachidosoraceae and the monophyly of Athyriaceae. Results of rate heterogeneity analysis revealed approximately 33 significant rate shifts in eupolypod II ferns, with the most heterogeneous rates (both accelerations and decelerations) occurring in two phylogenetically difficult lineages, that is, the Rhachidosoraceae-Aspleniaceae and Athyriaceae clades. These observations support the hypothesis that rate heterogeneity has previously constrained the deep phylogenetic resolution in eupolypods II. According to the plastome data, we propose that 14 chloroplast markers are particularly phylogenetically informative for eupolypods II both at the familial and generic levels. Our study demonstrates the power of a character-rich plastome data set and high-throughput sequencing for resolving the recalcitrant lineages, which have undergone rapid evolutionary radiation and dramatic changes in substitution rates.

Elevation Shift in Abies Mill. (Pinaceae) of Subtropical and Temperate China and Vietnam-Corroborative Evidence from Cytoplasmic DNA and Ecological Niche Modeling

The "elevational shift" scenario has been proposed as a model to explain the response of cold-adapted organisms to Quaternary climatic fluctuations in Europe and North America. However, the elevational shift model has not been well-explored in eastern Asia, which is more topographically complex than the other Northern Hemisphere biogeographic regions. Here, we evaluated the role of elevational shift in the closely related firs, or Abies Mill., of subtropical and temperate China. These firs are typical alpine trees with sensitivity to climate change. We tested the elevational shift hypothesis in firs of China using phylogeographic methods and ecological niche models. Our phylogeographic analyses comprised mitochondrial and chloroplast polymorphisms surveyed across 479 individuals from 43 populations representing 11 species. M1 of the 11 mitotypes and C1 of the 25 chlorotypes were inferred as the ancestral haplotype, and they had the widest distribution. The results of our phylogeographic survey revealed multiple centers of genetic diversity in distinct geographic regions and no latitudinal trend. Moreover, our results showed range expansions for seven taxa during the last glacial (64.9-18.2 or 32.5-9.1 kya), and this was consistent with the Quaternary fossil record of Abies in China. Taken together, our findings support a historical biogeographic pattern in firs of glacial expansions, probably through corridors at lower elevation, and interglacial fragmentations, through isolation at higher elevation peaks. Therefore, Abies in China probably undergoes elevational shift in response to climate change. Facing the forecasting global warming, the risk of several critically endangered firs was further enhanced as these species would have little escape space in situ to higher altitudes. According to our ENMs, we proposed an ex situ conservation strategy in the southern Hengduan Mountains region of south western China.

Isolation and characterization of microsatellite markers in the Lepisorus clathratus complex (Polypodiaceae)

PREMISE OF THE STUDY

Microsatellites were designed and characterized in the Sino-Himalayan fern Lepisorus clathratus complex (Polypodiaceae) to further study the phylogeography and reproductive ecology of this species.

METHODS AND RESULTS

From a genomic library obtained by next-generation sequencing, 10 polymorphic and six monomorphic microsatellite loci were developed. In one population of L. clathratus from Taibaishan in central China, the number of alleles observed for these microsatellites ranged from seven to 29, and observed and expected heterozygosity ranged from 0.463 to 0.919 and from 0.797 to 0.947, respectively. Cross-amplification in other taxa within this complex was successful, but cross-amplification was poor for other congeneric species.

CONCLUSIONS

This set of newly developed microsatellite markers will be useful for assessing genetic diversity, population structure, and mating system, and to infer polyploid origin in the L. clathratus complex.

The rise of the Himalaya enforced the diversification of SE Asian ferns by altering the monsoon regimes

BACKGROUND

The rise of high mountain chains is widely seen as one of the factors driving rapid diversification of land plants and the formation of biodiversity hotspots. Supporting evidence was reported for the impact of the rapid rise of the Andean mountains but this hypothesis has so far been less explored for the impact of the "roof of the world". The formation of the Himalaya, and especially the rise of the Qinghai-Tibetan Plateau in the recent 20 million years, altered the monsoon regimes that dominate the current climates of South East Asia. Here, we infer the hypothesis that the rise of Himalaya had a strong impact on the plant diversity in the biodiversity hotspot of the Southwest Chinese Mountains.

RESULTS

Our analyses of the diversification pattern of the derived fern genus Lepisorus recovered evidence for changes in plant diversity that correlated with the strengthening of South East Asian monsoon. Southwest China or Southwest China and Japan was recovered as the putative area of origin of Lepisorus and enhancing monsoon regime were found to shape the early diversification of the genus as well as subsequent radiations during the late Miocene and Pliocene.

CONCLUSIONS

We report new evidence for a coincidence of plant diversification and changes of the climate caused by the uplift of the Himalaya. These results are discussed in the context of the impact of incomplete taxon sampling, uncertainty of divergence time estimates, and limitations of current methods used to assess diversification rates.

Phylogeography of the Sino-Himalayan fern Lepisorus clathratus on "the roof of the world"

BACKGROUND

The Qinghai-Tibetan Plateau (QTP) and its southern and southeastern mountain ranges, Himalaya-Hengduan Mountains (HHM), are one of the most extensive habitats for alpine plants in the world. How ferns occurring in QTP and HHM changed their distribution ranges in response to Quaternary climatic oscillations remains almost unknown.

METHODOLOGY AND RESULTS

We employed sequences of two chloroplast DNA regions, rps4-trnS and trnL-trnF, to reconstruct phylogeography of the Sino-Himalayan fern Lepisorus clathratus, occurring mainly in the QTP and HHM. Individuals of this species have either dehiscent or indehiscent sporangia with the latter evolved from the plesiomorphic dehiscent forms. Based on a range-wide sampling, we detected 27 cpDNA haplotypes that were divided into five groups by network analyses. Populations in the Hengduan Mountains possess the highest genetic diversity, while a single haplogroup is detected across the north-central region. A distinct phylogeographical subdivision was detected between the Hengduan Mountains and north-central region by AMOVA analysis. The haplogroup distribution pattern, coalescence and AMOVA analysis suggest that a long term survival area (refugia) of the species was located in the Hengduan Mountains during glaciations, with probable range expansions into north-central regions during interglacial periods. Populations with indehiscent sporangium can carry private haplotypes and are inclined to maintain genetic homogeneity. One group with indehiscent sporangia most likely survived in situ on the QTP during glaciations.

CONCLUSIONS/SIGNIFICANCE

This study for the first time sheds light on the response of alpine ferns in the QTP and HHM to the Quaternary climatic oscillations.