The Small Nuclear Genomes of Selaginella Are Associated with a Low Rate of Genome Size Evolution

The Possible Earliest Allopolyploidization in Tracheophytes Revealed by Phylotranscriptomics and Morphology of Selaginellaceae

Selaginellaceae, originated in the Carboniferous and survived the Permian-Triassic mass extinction, is the largest family of lycophyte, which is sister to other tracheophytes. It stands out from tracheophytes by exhibiting extraordinary habitat diversity and lacking polyploidization. The organelle genome-based phylogenies confirmed the monophyly of Selaginella, with six or seven subgenera grouped into two superclades, but the phylogenetic positions of the enigmatic Selaginella sanguinolenta clade remained problematic. Here, we conducted a phylogenomic study on Selaginellaceae utilizing large-scale nuclear gene data from RNA-seq to elucidate the phylogeny and explore the causes of the phylogenetic incongruence of the S. sanguinolenta clade. Our phylogenetic analyses resolved three different positions of the S. sanguinolenta clade, which were supported by the sorted three nuclear gene sets, respectively. The results from the gene flow test, species network inference, and plastome-based phylogeny congruently suggested a probable hybrid origin of the S. sanguinolenta clade involving each common ancestor of the two superclades in Selaginellaceae. The hybrid hypothesis is corroborated by the evidence from rhizophore morphology and spore micromorphology. The chromosome observation and Ks distributions further suggested hybridization accompanied by polyploidization. Divergence time estimation based on independent datasets from nuclear gene sets and plastid genome data congruently inferred that allopolyploidization occurred in the Early Triassic. To our best knowledge, the allopolyploidization in the Mesozoic reported here represents the earliest record of tracheophytes. Our study revealed a unique triad of phylogenetic positions for a hybrid-originated group with comprehensive evidence and proposed a hypothesis for retaining both parental alleles through gene conversion.

Hybridization and reticulate evolution in Diphasiastrum (flat-branched clubmosses, Lycopodiaceae) - New data from the island of Taiwan and Vietnam

In the species groups related to Diphasiastrum multispicatum and D. veitchii, hybridization was investigated in samples from northern and southern Vietnam and the island of Taiwan, including available herbarium specimens from southeast Asia. The accessions were analyzed using flow cytometry (living material only), Sanger sequencing and multiplexed inter-simple sequence repeat genotyping by sequencing. We detected two cases of ancient hybridization involving different combinations of parental species; both led via subsequent duplication to tetraploid taxa. A cross D. multispicatum × D. veitchii from Malaysia represents D. wightianum, a tetraploid taxon according to reported DNA content measurements of dried material (genome formulas MM, VV and MMVV, respectively). The second case involves D. veitchii and an unknown diploid parent (genome formula XX). Three hybridogenous taxa (genome formulas VVX, VVXX, VVVX) were discernable by a combination of flow cytometry and molecular data. Taxon I (VVX, three clones found on Taiwan island) is apparently triploid. Taxon II represents another genetically diverse and sexual tetraploid species (VVXX) and can be assigned to D. yueshanense, described from Taiwan island but occurring as well in mainland China and Vietnam. Taxon III is as well most likely tetraploid (VVVX) and represented by at least one, more likely two, clones from Taiwan island. Taxa I and III are presumably asexual and new to science. Two independently inherited nuclear markers recombine only within, not between these hybrids, pointing towards reproductive isolation. We present an evolutionary scheme which explains the origin of the hybrids and the evolution of new and fully sexual species by hybridization and subsequent allopolyploidization in flat-branched clubmosses.

Nuclear DNA Amounts in Chinese Bryophytes Estimated by Flow Cytometry: Variation Patterns and Biological Significances

There exists an obvious gap in our knowledge of the nuclear DNA amount of bryophytes, not only in terms of the low number of species represented, but also in systematic and geographic representation. In order to increase our knowledge of nuclear DNA amounts and variation patterns in bryophytes, and their potential phylogenetic significances and influences on phenotypes, we used flow cytometry to determine the DNA 1C values of 209 bryophyte accessions, which belong to 145 mosses and 18 liverworts collected from China, by using Physcomitrella patens as a standard. We quantified the differences in DNA 1C values among different orders and families and constructed a phylogenetic tree of 112 mosses with four gene sequences (nad5, rbcL, trnL-F, and 18S-ITS1-5.8S-ITS2-26S). DNA 1C values were mapped onto the phylogenetic tree to test a potential phylogenetic signal. We also evaluated the correlations of the DNA 1C value with the sizes of individuals, leaves, cells, and spores by using a phylogenetically controlled analysis. New estimates of nuclear DNA amounts were reported for 145 species. The DNA 1C values of 209 bryophyte accessions ranged from 0.422 pg to 0.860 pg, with an average value of 0.561 pg, and a 2.04-fold variation covered the extremes of all the accessions. Although the values are not significantly different (p = 0.355) between mosses (0.528 pg) and liverworts (0.542 pg), there are variations to varying extents between some families and orders. The DNA 1C value size exerts a positive effect on the sizes of plants, leaves, and cells, but a negative effect on spore size. A weak phylogenetic signal is detected across most moss species. Phylogenetic signals are comparatively strong for some lineages. Our findings show that bryophytes have very small and highly constrained nuclear DNA amounts. There are nucleotype effects of nuclear DNA amounts for bryophytes at the individual, organ, and cell levels. We speculate that smaller nuclear DNA amounts are advantageous for bryophytes in dry environments. Significant differences in the DNA 1C values among some moss families and orders, as well as phylogenetic signals for some lineages, imply that nuclear DNA amount evolution in mosses seems to be unidirectional.

Genome Evolution in Plants: Complex Thalloid Liverworts (Marchantiopsida)

Why do some genomes stay small and simple, while others become huge, and why are some genomes more stable? In contrast to angiosperms and gymnosperms, liverworts are characterized by small genomes with low variation in size and conserved chromosome numbers. We quantified genome evolution among five Marchantiophyta (liverworts), measuring gene characteristics, transposable element (TE) landscape, collinearity, and sex chromosome evolution that might explain the small size and limited variability of liverwort genomes. No genome duplications were identified among examined liverworts and levels of duplicated genes are low. Among the liverwort species, Lunularia cruciata stands out with a genome size almost twice that of the other liverwort species investigated here, and most of this increased size is due to bursts of Ty3/Gypsy retrotransposons. Intrachromosomal rearrangements between examined liverworts are abundant but occur at a slower rate compared with angiosperms. Most genes on L. cruciata scaffolds have their orthologs on homologous Marchantia polymorpha chromosomes, indicating a low degree of rearrangements between chromosomes. Still, translocation of a fragment of the female U chromosome to an autosome was predicted from our data, which might explain the uniquely small U chromosome in L. cruciata. Low levels of gene duplication, TE activity, and chromosomal rearrangements might contribute to the apparent slow rate of morphological evolution in liverworts.

Genome size evolution of the extant lycophytes and ferns

Ferns and lycophytes have remarkably large genomes. However, little is known about how their genome size evolved in fern lineages. To explore the origins and evolution of chromosome numbers and genome size in ferns, we used flow cytometry to measure the genomes of 240 species (255 samples) of extant ferns and lycophytes comprising 27 families and 72 genera, of which 228 species (242 samples) represent new reports. We analyzed correlations among genome size, spore size, chromosomal features, phylogeny, and habitat type preference within a phylogenetic framework. We also applied ANOVA and multinomial logistic regression analysis to preference of habitat type and genome size. Using the phylogeny, we conducted ancestral character reconstruction for habitat types and tested whether genome size changes simultaneously with shifts in habitat preference. We found that 2C values had weak phylogenetic signal, whereas the base number of chromosomes (x) had a strong phylogenetic signal. Furthermore, our analyses revealed a positive correlation between genome size and chromosome traits, indicating that the base number of chromosomes (x), chromosome size, and polyploidization may be primary contributors to genome expansion in ferns and lycophytes. Genome sizes in different habitat types varied significantly and were significantly correlated with habitat types; specifically, multinomial logistic regression indicated that species with larger 2C values were more likely to be epiphytes. Terrestrial habitat is inferred to be ancestral for both extant ferns and lycophytes, whereas transitions to other habitat types occurred as the major clades emerged. Shifts in habitat types appear be followed by periods of genomic stability. Based on these results, we inferred that habitat type changes and multiple whole-genome duplications have contributed to the formation of large genomes of ferns and their allies during their evolutionary history.

Ammonium metabolism in Selaginella bryopteris in response to dehydration-rehydration and characterisation of desiccation tolerant, thermostable, cytosolic glutamine synthetase from plant

Water deficit (WD) has adverse effects on plant growth, and acclimation requires responses allowing primary metabolism to continue. Resurrection plants can serve as model system to gain insight into metabolic regulation during WD. We herein report the response of a resurrection lycophyte, Selaginella bryopteris, to dehydration-rehydration cycle with emphasis on ammonium metabolism. Dehydration of S. bryopteris fronds resulted in decrease of total protein and increase of free ammonium levels and the effect was reversed on rehydration. The proline content increased twice after 24 h of dehydration, which again recovered to background levels comparable to that at full turgor state. The specific activity of glutamine synthetase (GS) didn't change significantly till 6 h and then declined by 21% after 24 h of dehydration, whereas specific activities of glutamate synthase (GOGAT) and aminating glutamate dehydrogenase (GDH) were enhanced significantly during dehydration. The deaminating activity of GDH also increased during dehydration albeit at a slower rate. Immunoblot analysis indicated overexpression of GS and GDH polypeptides during dehydration and their levels declined on rehydration. The results suggested significant role of GDH along with GS/GOGAT in production of nitrogen-rich amino acids for desiccation tolerance. Unlike higher plants S. bryopteris expressed GS only in cytosol. The enzyme had pH and temperature optima of 5.5 and 60°C, respectively, and it retained 96% activity on preincubation at 60°C for 30 min indicating thermostability. Hence, like higher plants the cytosolic GS from S. bryopteris has a conserved role in stress tolerance.

Comparative transcriptome analysis suggests convergent evolution of desiccation tolerance in Selaginella species

BACKGROUND

Desiccation tolerant Selaginella species evolved to survive extreme environmental conditions. Studies to determine the mechanisms involved in the acquisition of desiccation tolerance (DT) have focused on only a few Selaginella species. Due to the large diversity in morphology and the wide range of responses to desiccation within the genus, the understanding of the molecular basis of DT in Selaginella species is still limited.

RESULTS

Here we present a reference transcriptome for the desiccation tolerant species S. sellowii and the desiccation sensitive species S. denticulata. The analysis also included transcriptome data for the well-studied S. lepidophylla (desiccation tolerant), in order to identify DT mechanisms that are independent of morphological adaptations. We used a comparative approach to discriminate between DT responses and the common water loss response in Selaginella species. Predicted proteomes show strong homology, but most of the desiccation responsive genes differ between species. Despite such differences, functional analysis revealed that tolerant species with different morphologies employ similar mechanisms to survive desiccation. Significant functions involved in DT and shared by both tolerant species included induction of antioxidant systems, amino acid and secondary metabolism, whereas species-specific responses included cell wall modification and carbohydrate metabolism.

CONCLUSIONS

Reference transcriptomes generated in this work represent a valuable resource to study Selaginella biology and plant evolution in relation to DT. Our results provide evidence of convergent evolution of S. sellowii and S. lepidophylla due to the different gene sets that underwent selection to acquire DT.

Paleo-polyploidization in Lycophytes

Lycophytes and seed plants constitute the typical vascular plants. Lycophytes have been thought to have no paleo-polyploidization although the event is known to be critical for the fast expansion of seed plants. Here, genomic analyses including the homologous gene dot plot analysis detected multiple paleo-polyploidization events, with one occurring approximately 13-15 million years ago (MYA) and another about 125-142 MYA, during the evolution of the genome of Selaginella moellendorffii, a model lycophyte. In addition, comparative analysis of reconstructed ancestral genomes of lycophytes and angiosperms suggested that lycophytes were affected by more paleo-polyploidization events than seed plants. Results from the present genomic analyses indicate that paleo-polyploidization has contributed to the successful establishment of both lineages-lycophytes and seed plants-of vascular plants.

Sequencing and Analyzing the Transcriptomes of a Thousand Species Across the Tree of Life for Green Plants

The 1,000 Plants (1KP) initiative was the first large-scale effort to collect next-generation sequencing (NGS) data across a phylogenetically representative sampling of species for a major clade of life, in this case theViridiplantae, or green plants. As an international multidisciplinary consortium, we focused on plant evolution and its practical implications. Among the major outcomes were the inference of a reference species tree for green plants by phylotranscriptomic analysis of low-copy genes, a survey of paleopolyploidy (whole-genome duplications) across the Viridiplantae, the inferred evolutionary histories for many gene families and biological processes, the discovery of novel light-sensitive proteins for optogenetic studies in mammalian neuroscience, and elucidation of the genetic network for a complex trait (C₄ photosynthesis). Altogether, 1KP demonstrated how value can be extracted from a phylodiverse sequencing data set, providing a template for future projects that aim to generate even more data, including complete de novo genomes, across the tree of life.

Genome size and endopolyploidy evolution across the moss phylogeny

BACKGROUND AND AIMS

Compared with other plant lineages, bryophytes have very small genomes with little variation across species, and high levels of endopolyploid nuclei. This study is the first analysis of moss genome evolution over a broad taxonomic sampling using phylogenetic comparative methods. We aim to determine whether genome size evolution is unidirectional as well as examine whether genome size and endopolyploidy are correlated in mosses.

METHODS

Genome size and endoreduplication index (EI) estimates were newly generated using flow cytometry from moss samples collected in Canada. Phylogenetic relationships between moss species were reconstructed using GenBank sequence data and maximum likelihood methods. Additional 1C-values were compiled from the literature and genome size and EI were mapped onto the phylogeny to reconstruct ancestral character states, test for phylogenetic signal and perform phylogenetic independent contrasts.

KEY RESULTS

Genome size and EI were obtained for over 50 moss taxa. New genome size estimates are reported for 33 moss species and new EIs are reported for 20 species. In combination with data from the literature, genome sizes were mapped onto a phylogeny for 173 moss species with this analysis, indicating that genome size evolution in mosses does not appear to be unidirectional. Significant phylogenetic signal was detected for genome size when evaluated across the phylogeny, whereas phylogenetic signal was not detected for EI. Genome size and EI were not found to be significantly correlated when using phylogenetically corrected values.

CONCLUSIONS

Significant phylogenetic signal indicates closely related mosses have more similar genome sizes and EI values. This study supports that DNA content in mosses is defined by small genomes that are highly endopolyploid, suggesting strong selective pressure to maintain these features. Further research is needed to understand the functional significance of DNA content evolution in mosses.

GingerRoot: A Novel DNA Transposon Encoding Integrase-Related Transposase in Plants and Animals

Transposable elements represent the largest components of many eukaryotic genomes and different genomes harbor different combinations of elements. Here, we discovered a novel DNA transposon in the genome of the clubmoss Selaginella lepidophylla. Further searching for related sequences to the conserved DDE region uncovered the presence of this superfamily of elements in fish, coral, sea anemone, and other animal species. However, this element appears restricted to Bryophytes and Lycophytes in plants. This transposon, named GingerRoot, is associated with a 6 bp (base pair) target site duplication, and 100-150 bp terminal inverted repeats. Analysis of transposase sequences identified the DDE motif, a catalytic domain, which shows similarity to the integrase of Gypsy-like long terminal repeat retrotransposons, the most abundant component in plant genomes. A total of 77 intact and several hundred truncated copies of GingerRoot elements were identified in S. lepidophylla. Like Gypsy retrotransposons, GingerRoots show a lack of insertion preference near genes, which contrasts to the compact genome size of about 100 Mb. Nevertheless, a considerable portion of GingerRoot elements was found to carry gene fragments, suggesting the capacity of duplicating gene sequences is unlikely attributed to the proximity to genes. Elements carrying gene fragments appear to be less methylated, more diverged, and more distal to genes than those without gene fragments, indicating they are preferentially retained in gene-poor regions. This study has identified a broadly dispersed, novel DNA transposon, and the first plant DNA transposon with an integrase-related transposase, suggesting the possibility of de novo formation of Gypsy-like elements in plants.

Giberellic Acid-Stimulated Transcript Proteins Evolved through Successive Conjugation of Novel Motifs and Their Subfunctionalization

Gibberellic Acid Stimulated Transcript (GAST)-like genes encode small polypeptides, some of which have been implicated in diverse biological processes regulating plant growth and development. However, the occurrence of GASTs among plants, their protein structures, and the mechanisms by which they evolved remain elusive. Here, using a customized workflow, we report genes encoding GAST proteins, identify novel motifs and evolutionary patterns contributing to the subfunctionalization of GAST domains, and explore functional conservation across diverse plant groups. We show that GAST-like sequences evolved initially in the vascular plant Selaginella moellendorffii, after the divergence from bryophytes, and later emerged in gymnosperms and angiosperms. GASTs in angiosperms are characterized by four conserved novel motifs; however, relatively fewer conserved motifs exist in pteridophytes and gymnosperms. Phylogenetic analysis revealed that the GAST-Cysteine Rich1 motif evolved early in the S. moellendorffii GAST, which further acquired subfunctionalization through successive conjugation of other motifs and remained conserved across plants, as supported by their collinearity. Functional characterization of two orthologs from the dicot Arabidopsis (Arabidopsis thaliana; Gibberellic Acid-Stimulated Arabidopsis 10) and the monocot rice (Oryza sativa; Gibberellic Acid Stimulated Transcript-Related 9) suggests hormonal regulation, novel roles in seed germination, and functional conservation among diverse plant groups. Computational modeling predicts that these GAST genes are regulated by several factors, including the phytohormones gibberellin and abscisic acid, through conserved cis-motifs present in their promoters, and that they might act as signaling molecules in a complex feedback loop. Thus, our study identifies GASTs and their encoded proteins, uncovers their structure, novel motifs, and evolutionary pattern among plants, and suggests their functional conservation.

Phylogenetic trends and environmental correlates of nuclear genome size variation in Helianthus sunflowers

Flowering plants serve as a powerful model for studying the evolution of nuclear genome size (GS) given the tremendous GS variation that exists both within and across angiosperm lineages. Helianthus sunflowers consist of c. 50 species native to North America that occupy diverse habitats and vary in ploidy level. In the current study, we generated a comprehensive GS database for 49 Helianthus species using flow cytometric approaches. We examined variability across the genus and present a comparative phylogenetic analysis of GS evolution in diploid Helianthus species. Results demonstrated that different clades of diploid Helianthus species showed evolutionary patterns of GS contraction, expansion and relative stasis, with annual diploid species evolving smaller GS with the highest rate of evolution. Phylogenetic comparative analyses of diploids revealed significant negative associations of GS with temperature seasonality and cell production rate, indicating that the evolution of larger GS in Helianthus diploids may be more permissible in habitats with longer growing seasons where selection for more rapid growth may be relaxed. The Helianthus GS database presented here and corresponding analyses of environmental and phenotypic correlates will facilitate ongoing and future research on the ultimate drivers of GS evolution in this well-studied North American plant genus.

Desiccation Tolerance Evolved through Gene Duplication and Network Rewiring in Lindernia

Although several resurrection plant genomes have been sequenced, the lack of suitable dehydration-sensitive outgroups has limited genomic insights into the origin of desiccation tolerance. Here, we utilized a comparative system of closely related desiccation-tolerant (Lindernia brevidens) and -sensitive (Lindernia subracemosa) species to identify gene- and pathway-level changes associated with the evolution of desiccation tolerance. The two high-quality Lindernia genomes we assembled are largely collinear, and over 90% of genes are conserved. L. brevidens and L. subracemosa have evidence of an ancient, shared whole-genome duplication event, and retained genes have neofunctionalized, with desiccation-specific expression in L. brevidens Tandem gene duplicates also are enriched in desiccation-associated functions, including a dramatic expansion of early light-induced proteins from 4 to 26 copies in L. brevidens A comparative differential gene coexpression analysis between L. brevidens and L. subracemosa supports extensive network rewiring across early dehydration, desiccation, and rehydration time courses. Many LATE EMBRYOGENESIS ABUNDANT genes show significantly higher expression in L. brevidens compared with their orthologs in L. subracemosa Coexpression modules uniquely upregulated during desiccation in L. brevidens are enriched with seed-specific and abscisic acid-associated cis-regulatory elements. These modules contain a wide array of seed-associated genes that have no expression in the desiccation-sensitive L. subracemosa Together, these findings suggest that desiccation tolerance evolved through a combination of gene duplications and network-level rewiring of existing seed desiccation pathways.

Identification and Analysis of OVATE Family Members from Genome of the Early Land Plants Provide Insights into Evolutionary History of OFP Family and Function

Mosses, liverworts, hornworts and lycophytes represent transition stages between the aquatic to terrestrial/land plants. Several morphological and adaptive novelties driven by genomic components including emergence and expansion of new or existing gene families have played a critical role during and after the transition, and contributed towards successful colonization of terrestrial ecosystems. It is crucial to decipher the evolutionary transitions and natural selection on the gene structure and function to understand the emergence of phenotypic and adaptive diversity. Plants at the "transition zone", between aquatic and terrestrial ecosystem, are also the most vulnerable because of climate change and may contain clues for successful mitigation of the challenges of climate change. Identification and comparative analyses of such genetic elements and gene families are few in mosses, liverworts, hornworts and lycophytes. Ovate family proteins (OFPs) are plant-specific transcriptional repressors and are acknowledged for their roles in important growth and developmental processes in land plants, and information about the functional aspects of OFPs in early land plants is fragmentary. As a first step towards addressing this gap, a comprehensive in silico analysis was carried out utilizing publicly available genome sequences of Marchantia polymorpha (Mp), Physcomitrella patens (Pp), Selaginella moellendorffii (Sm) and Sphagnum fallax (Sf). Our analysis led to the identification of 4 MpOFPs, 19 PpOFPs, 6 SmOFPs and 3 SfOFPs. Cross-genera analysis revealed a drastic change in the structure and physiochemical properties in OFPs suggesting functional diversification and genomic plasticity during the evolutionary course. Knowledge gained from this comparative analysis will form the framework towards deciphering and dissection of their developmental and adaptive role/s in early land plants and could provide insights into evolutionary strategies adapted by land plants.

Genome Analysis of the Ancient Tracheophyte Selaginella tamariscina Reveals Evolutionary Features Relevant to the Acquisition of Desiccation Tolerance

Resurrection plants, which are the "gifts" of natural evolution, are ideal models for studying the genetic basis of plant desiccation tolerance. Here, we report a high-quality genome assembly of 301 Mb for the diploid spike moss Selaginella tamariscina, a primitive vascular resurrection plant. We predicated 27 761 protein-coding genes from the assembled S. tamariscina genome, 11.38% (2363) of which showed significant expression changes in response to desiccation. Approximately 60.58% of the S. tamariscina genome was annotated as repetitive DNA, which is an almost 2-fold increase of that in the genome of desiccation-sensitive Selaginella moellendorffii. Genomic and transcriptomic analyses highlight the unique evolution and complex regulations of the desiccation response in S. tamariscina, including species-specific expansion of the oleosin and pentatricopeptide repeat gene families, unique genes and pathways for reactive oxygen species generation and scavenging, and enhanced abscisic acid (ABA) biosynthesis and potentially distinct regulation of ABA signaling and response. Comparative analysis of chloroplast genomes of several Selaginella species revealed a unique structural rearrangement and the complete loss of chloroplast NAD(P)H dehydrogenase (NDH) genes in S. tamariscina, suggesting a link between the absence of the NDH complex and desiccation tolerance. Taken together, our comparative genomic and transcriptomic analyses reveal common and species-specific desiccation tolerance strategies in S. tamariscina, providing significant insights into the desiccation tolerance mechanism and the evolution of resurrection plants.

Genome Size Diversity and Its Impact on the Evolution of Land Plants

Genome size is a biodiversity trait that shows staggering diversity across eukaryotes, varying over 64,000-fold. Of all major taxonomic groups, land plants stand out due to their staggering genome size diversity, ranging ca. 2400-fold. As our understanding of the implications and significance of this remarkable genome size diversity in land plants grows, it is becoming increasingly evident that this trait plays not only an important role in shaping the evolution of plant genomes, but also in influencing plant community assemblages at the ecosystem level. Recent advances and improvements in novel sequencing technologies, as well as analytical tools, make it possible to gain critical insights into the genomic and epigenetic mechanisms underpinning genome size changes. In this review we provide an overview of our current understanding of genome size diversity across the different land plant groups, its implications on the biology of the genome and what future directions need to be addressed to fill key knowledge gaps.

Extreme haplotype variation in the desiccation-tolerant clubmoss Selaginella lepidophylla

Plant genome size varies by four orders of magnitude, and most of this variation stems from dynamic changes in repetitive DNA content. Here we report the small 109 Mb genome of Selaginella lepidophylla, a clubmoss with extreme desiccation tolerance. Single-molecule sequencing enables accurate haplotype assembly of a single heterozygous S. lepidophylla plant, revealing extensive structural variation. We observe numerous haplotype-specific deletions consisting of largely repetitive and heavily methylated sequences, with enrichment in young Gypsy LTR retrotransposons. Such elements are active but rapidly deleted, suggesting "bloat and purge" to maintain a small genome size. Unlike all other land plant lineages, Selaginella has no evidence of a whole-genome duplication event in its evolutionary history, but instead shows unique tandem gene duplication patterns reflecting adaptation to extreme drying. Gene expression changes during desiccation in S. lepidophylla mirror patterns observed across angiosperm resurrection plants.