Piece and parcel of gymnosperm organellar genomes.

MAIN CONCLUSION: Significant past, present, and potential future research into the organellar (plastid and mitochondrial) genomes of gymnosperms that can provide insight into the unknown origin and evolution of plants is highlighted. Gymnosperms are vascular seed plants that predominated the ancient world before their sister clade, angiosperms, took over during the Late Cretaceous. The divergence of gymnosperms and angiosperms took place around 300 Mya, with the latter evolving into the diverse group of flowering plants that dominate the plant kingdom today. Although gymnosperms have reportedly made some evolutionary innovations, the literature on their genome advances, particularly their organellar (plastid and mitochondrial) genomes, is relatively scattered and fragmented. While organellar genomes can shed light on plant origin and evolution, they are frequently overlooked, due in part to their limited contribution to gene expression and lack of evolutionary dynamics when compared to nuclear genomes. A better understanding of gymnosperm organellar genomes is critical because they reveal genetic changes that have contributed to their unique adaptations and ecological success, potentially aiding in plant survival, enhancement, and biodiversity conservation in the face of climate change. This review reveals significant information and gaps in the existing knowledge base of organellar genomes in gymnosperms, as well as the challenges and research needed to unravel their complexity.

Complete organelle genomes of Korean fir, Abies koreana and phylogenomics of the gymnosperm genus Abies using nuclear and cytoplasmic DNA sequence data.

Abies koreana E.H.Wilson is an endangered evergreen coniferous tree that is native to high altitudes in South Korea and susceptible to the effects of climate change. Hybridization and reticulate evolution have been reported in the genus; therefore, multigene datasets from nuclear and cytoplasmic genomes are needed to better understand its evolutionary history. Using the Illumina NovaSeq 6000 and Oxford Nanopore Technologies (ONT) PromethION platforms, we generated complete mitochondrial (1,174,803 bp) and plastid (121,341 bp) genomes from A. koreana. The mitochondrial genome is highly dynamic, transitioning from cis- to trans-splicing and breaking conserved gene clusters. In the plastome, the ONT reads revealed two structural conformations of A. koreana. The short inverted repeats (1186 bp) of the A. koreana plastome are associated with different structural types. Transcriptomic sequencing revealed 1356 sites of C-to-U RNA editing in the 41 mitochondrial genes. Using A. koreana as a reference, we additionally produced nuclear and organelle genomic sequences from eight Abies species and generated multiple datasets for maximum likelihood and network analyses. Three sections (Balsamea, Momi, and Pseudopicea) were well grouped in the nuclear phylogeny, but the phylogenomic relationships showed conflicting signals in the mitochondrial and plastid genomes, indicating a complicated evolutionary history that may have included introgressive hybridization. The obtained data illustrate that phylogenomic analyses based on sequences from differently inherited organelle genomes have resulted in conflicting trees. Organelle capture, organelle genome recombination, and incomplete lineage sorting in an ancestral heteroplasmic individual can contribute to phylogenomic discordance. We provide strong support for the relationships within Abies and new insights into the phylogenomic complexity of this genus.

Disparity of cycad leaves dispels the living fossil metaphor.

The living fossil metaphor is tightly linked with the cycads. This group of gymnosperms is supposed to be characterised by long-term morphological stasis, particularly after their peak of diversity and disparity in the Jurassic. However, no formal test of this hypothesis exists. Here, we use a recent phylogenetic framework and an improved character matrix to reconstruct the Disparity Through Time for cycad leaves using a Principal Coordinate Analysis and employing Pre-Ordination Ancestral State Reconstruction to test the impact of sampling on the results. Our analysis shows that the cycad leaf morsphospace expanded up to the present, with numerous shifts in its general positioning, independently of sampling biases. Moreover, they also show that Zamiaceae expanded rapidly in the Early Cretaceous and continued to expand up to the present, while now-extinct clades experienced a slow contraction from their peak in the Triassic. We also show that rates of evolution were constantly high up to the Early Cretaceous, and then experienced a slight decrease in the Paleogene, followed by a Neogene acceleration. These results show a much more dynamic history for cycads, and suggest that the 'living fossil' metaphor is actually a hindrance to our understanding of their macroevolution.

Hybridization mediated range expansion and climate change resilience in two keystone tree species of boreal forests.

Current global climate change is expected to affect biodiversity negatively at all scales leading to mass biodiversity loss. Many studies have shown that the distribution of allele frequencies across a species' range is often influenced by specific genetic loci associated with local environmental variables. This association reflects local adaptation and allele changes at those loci could thereby contribute to the evolutionary response to climate change. However, predicting how species will adapt to climate change from this type of data alone remains challenging. In the present study, we combined exome capture sequences and environmental niche reconstruction, to test multiple methods for assessing local adaptation and climate resilience in two widely distributed conifers, Norway spruce and Siberian spruce. Both species are keystone species of the boreal forest and share a vast hybrid zone. We show that local adaptation in conifers can be detected through allele frequency variation, population-level ecological preferences, and historical niche movement. Moreover, we integrated genetic and ecological information into genetic offset predictive models to show that hybridization plays a central role in expanding the niche breadth of the two conifer species and may help both species to cope better with future changing climates. This joint genetic and ecological analysis also identified spruce populations that are at risk under current climate change.

Early-acting inbreeding depression can evolve as an inbreeding avoidance mechanism.

Despite the potential for mechanical, developmental and/or chemical mechanisms to prevent self-fertilization, incidental self-fertilization is inevitable in many predominantly outcrossing species. In such cases, inbreeding can compromise individual fitness. Unquestionably, much of this inbreeding depression is maladaptive. However, we show that when reproductive compensation allows for the replacement of inviable embryos lost early in development, selection can favour deleterious recessive variants that induce 'self-sacrificial' death of inbred embryos. Our theoretical results provide numerous testable predictions which could challenge the assumption that inbreeding depression is always maladaptive. Our work is applicable any species that cannot fully avoid inbreeding, exhibits substantial inbreeding depression, and has the potential to compensate embryos lost early in development. In addition to its general applicability, our theory suggests that self-sacrificial variants might be responsible for the remarkably low realized selfing rates of gymnosperms with high primary selfing rates, as gymnosperms exhibit strong inbreeding depression, have effective reproductive compensation mechanisms, and cannot evolve chemical self-incompatibility.

Nitrogen addition affected the root competition in Cunninghamia lanceolata-Phoebe chekiangensis mixed plantation.

Little is known about below-ground competition in mixed-species plantations under increasing nitrogen (N) deposition. This study aims to determine the effects of N addition on root competition in coniferous and broad-leaved species mixed plantations. A pot experiment was conducted using the coniferous species Cunninghamia lanceolata and the broad-leaved species Phoebe chekiangensis planted in mixed plantations with different competition intensities under N addition (0 or 45 kg N ha-1 yr-1). Biomass allocation, root morphology, root growth level, and competitive ability were determined after five months of treatment. Our findings indicated that root interactions in mixed plantations did not influence biomass allocation in either C. lanceolata or P. chekiangensis but promoted growth in C. lanceolata when no N was added. However, N addition decreased biomass accumulation in both species in the mixed plantation and had a negative effect on the root growth of C. lanceolata due to intensified competition. Addition of N increased the relative importance of root predatory competition in P. chekiangensis, and increased the allelopathic competitive advantage in C. lanceolata. This suggests that N addition causes a shift in the root competitive strategy from tolerance to competition. Overall, these findings highlight the significant impact that the addition of N can have on plant interactions in mixed plantations. Our results provide implications for the mechanisms of root competition in response to increasing atmospheric N deposition in mixed plantations.

Epigenetic stress memory in gymnosperms.

Gymnosperms are long-lived, cone-bearing seed plants that include some of the most ancient extant plant species. These relict land plants have evolved to survive in habitats marked by chronic or episodic stress. Their ability to thrive in these environments is partly due to their phenotypic flexibility, and epigenetic regulation likely plays a crucial part in this plasticity. We review the current knowledge on abiotic and biotic stress memory in gymnosperms and the possible epigenetic mechanisms underlying long-term phenotypic adaptations. We also discuss recent technological improvements and new experimental possibilities that likely will advance our understanding of epigenetic regulation in these ancient and hard-to-study plants.

Conserved autophagy and diverse cell wall composition: unifying features of vascular tissues in evolutionarily distinct plants.

BACKGROUND AND AIMS: The formation of multifunctional vascular tissues represents a significant advancement in plant evolution. Differentiation of conductive cells is specific, involving two main pathways, namely protoplast clearance and cell wall modification. In xylogenesis, autophagy is a crucial process for complete protoplast elimination in tracheary elements, whose cell wall also undergoes strong changes. Knowledge pertaining to living sieve elements, which lose most of their protoplast during phloemogenesis, remains limited. We hypothesized that autophagy plays a crucial role, not only in complete cytoplasmic clearance in xylem but also in partial degradation in phloem. Cell wall elaborations of mature sieve elements are not so extensive. These analyses performed on evolutionarily diverse model species potentially make it possible to understand phloemogenesis to an equal extent to xylogenesis. METHODS: We investigated the distribution of ATG8 protein, which is an autophagy marker, and cell wall components in the roots of ferns, gymnosperms and angiosperms (monocots, dicot herbaceous plants and trees). Furthermore, we conducted a bioinformatic analysis of complete data on ATG8 isoforms for Ceratopteris richardii. KEY RESULTS: The presence of ATG8 protein was confirmed in both tracheary elements and sieve elements; however, the composition of cell wall components varied considerably among vascular tissues in the selected plants. Arabinogalactan proteins and ß-1,4-galactan were detected in the roots of all studied species, suggesting their potential importance in phloem formation or function. In contrast, no evolutionary pattern was observed for xyloglucan, arabinan or homogalacturonan. CONCLUSIONS: Our findings indicate that the involvement of autophagy in plants is universal during the development of tracheary elements that are dead at maturity and sieve elements that remain alive. Given the conserved nature of autophagy and its function in protoplast degradation for uninterrupted flow, autophagy might have played a vital role in the development of increasingly complex biological organizations, including the formation of vascular tissues. However, different cell wall compositions of xylem and phloem in different species might indicate diverse functionality and potential for substance transport, which is crucial in plant evolution.

Molecular regulatory mechanisms of staminate strobilus development and dehiscence in Torreya grandis.

Gymnosperms are mostly dioecious, and their staminate strobili undergo a longer developmental period than those of angiosperms. However, the underlying molecular mechanisms remain unclear. This study aimed to identify key genes and pathways involved in staminate strobilus development and dehiscence in Torreya grandis. Through weighted gene co-expression network analysis (WGCNA), we identified fast elongation-related genes enriched in carbon metabolism and auxin signal transduction, whereas dehiscence-related genes were abundant in alpha-linolenic acid metabolism and the phenylpropanoid pathway. Based on WGCNA, we also identified PHYTOCHROME-INTERACTING FACTOR4 (TgPIF4) as a potential regulator for fast elongation of staminate strobilus and 2 WRKY proteins (TgWRKY3 and TgWRKY31) as potential regulators for staminate strobilus dehiscence. Multiple protein-DNA interaction analyses showed that TgPIF4 directly activates the expression of TRANSPORT INHIBITOR RESPONSE2 (TgTIR2) and NADP-MALIC ENZYME (TgNADP-ME). Overexpression of TgPIF4 significantly promoted staminate strobilus elongation by elevating auxin signal transduction and pyruvate content. TgWRKY3 and TgWRKY31 bind to the promoters of the lignin biosynthesis gene PHENYLALANINE AMMONIA-LYASE (TgPAL) and jasmonic acid metabolism gene JASMONATE O-METHYLTRANSFERASE (TgJMT), respectively, and directly activate their transcription. Overexpression of TgWRKY3 and TgWRKY31 in the staminate strobilus led to early dehiscence, accompanied by increased lignin and methyl jasmonate levels, respectively. Collectively, our findings offer a perspective for understanding the growth of staminate strobili in gymnosperms.

Middle Jurassic insect mines on gymnosperms provide missing links to early mining evolution.

We investigated the mining mode of insect feeding, involving larval consumption of a plant's internal tissues, from the Middle Jurassic (165 million years ago) Daohugou locality of Northeastern China. Documentation of mining from the Jurassic Period is virtually unknown, and results from this time interval would address mining evolution during the temporal gap of mine-seed plant diversifications from the previous Late Triassic to the subsequent Early Cretaceous. Plant fossils were examined with standard microscopic procedures for herbivory and used the standard functional feeding group-damage-type system of categorizing damage. All fossil mines were photographed and databased. We examined 2014 plant specimens, of which 27 occurrences on 14 specimens resulted in eight, new, mine damage types (DTs) present on six genera of bennettitalean, ginkgoalean, and pinalean gymnosperms. Three conclusions emerge from this study. First, these mid-Mesozoic mines are morphologically conservative and track plant host anatomical structure rather than plant phylogeny. Second, likely insect fabricators of these mines were three basal lineages of polyphagan beetles, four basal lineages of monotrysian moths, and a basal lineage tenthredinoid sawflies. Third, the nutrition hypothesis, indicating that miners had greater access to nutritious, inner tissues of new plant lineages, best explains mine evolution during the mid-Mesozoic.

Mongolitria: A new Early Cretaceous three-valved seed from Northeast Asia.

PREMISE: Fossil seeds recovered from the Early Cretaceous of Mongolia and Inner Mongolia, China, are described and assigned to Mongolitria gen. nov., a new genus of gymnosperm seed. METHODS: Abundant lignitized seeds along with compression specimens isolated from the matrix were studied using a combination of scanning electron microscopy, anatomical sectioning, light microscopy, synchrotron radiation X-ray microtomography, and cuticle preparations. A single permineralized seed was examined by light microscopy of cellulose acetate peels and X-ray microtomography. RESULTS: Two species are recognized, Mongolitria friisae sp. nov. and Mongolitria exesum sp. nov. Both seeds are orthotropous with a short apical micropyle and a small, basal, circular attachment scar. The thick sclerenchymatous integument has a consistently three-parted organization and about 20 conspicuous longitudinal ribs on the surface. Mongolitria exesum differs from M. friisae primarily in its much larger size and thicker seed coat, which also preserves clear evidence of insect damage. CONCLUSIONS: Mongolitria is similar to other fossil seeds that have been assigned to Cycadales, but displays a unique combination of characters not found in any living or extinct cycadaceous plant, leaving its higher-level systematic affinities uncertain. Germination apparently involved splitting of the integument into three valves. Mongolitria was prominent among the plant parts accumulating in peat swamps in eastern Asia during the Early Cretaceous.

Nitrogen isotopes reveal independent origins of N2-fixing symbiosis in extant cycad lineages.

Cycads are ancient seed plants (gymnosperms) that emerged by the early Permian. Although they were common understory flora and food for dinosaurs in the Mesozoic, their abundance declined markedly in the Cenozoic. Extant cycads persist in restricted populations in tropical and subtropical habitats and, with their conserved morphology, are often called 'living fossils.' All surviving taxa receive nitrogen from symbiotic N2-fixing cyanobacteria living in modified roots, suggesting an ancestral origin of this symbiosis. However, such an ancient acquisition is discordant with the abundance of cycads in Mesozoic fossil assemblages, as modern N2-fixing symbioses typically occur only in nutrient-poor habitats where advantageous for survival. Here, we use foliar nitrogen isotope ratios-a proxy for N2 fixation in modern plants-to probe the antiquity of the cycad-cyanobacterial symbiosis. We find that fossilized cycad leaves from two Cenozoic representatives of extant genera have nitrogen isotopic compositions consistent with microbial N2 fixation. In contrast, all extinct cycad genera have nitrogen isotope ratios that are indistinguishable from co-existing non-cycad plants and generally inconsistent with microbial N2 fixation, pointing to nitrogen assimilation from soils and not through symbiosis. This pattern indicates that, rather than being ancestral within cycads, N2-fixing symbiosis arose independently in the lineages leading to living cycads during or after the Jurassic. The preferential survival of these lineages may therefore reflect the effects of competition with angiosperms and Cenozoic climatic change.

Quantifying the complexity of plant reproductive structures reveals a history of morphological and functional integration.

Vascular plant reproductive structures have undoubtedly become more complex through time, evolving highly differentiated parts that interact in specialized ways. But quantifying these patterns at broad scales is challenging because lineages produce disparate reproductive structures that are often difficult to compare and homologize. We develop a novel approach for analysing interactions within reproductive structures using networks, treating component parts as nodes and a suite of physical and functional interactions among parts as edges. We apply this approach to the plant fossil record, showing that interactions have generally increased through time and that the concentration of these interactions has shifted towards differentiated surrounding organs, resulting in more compact, functionally integrated structures. These processes are widespread across plant lineages, but their extent and timing vary with reproductive biology; in particular, seed-producing structures show them more strongly than spore or pollen-producing structures. Our results demonstrate that major reproductive innovations like the origin of seeds and angiospermy were associated with increased integration through greater interactions among parts. But they also reveal that for certain groups, particularly Mesozoic gymnosperms, millions of years elapsed between the origin of reproductive innovations and increased interactions among parts within their reproductive structures.

Museomics unveil systematics, diversity and evolution of Australian cycad-pollinating weevils.

Weevils have been shown to play significant roles in the obligate pollination of Australian cycads. In this study, we apply museomics to produce a first molecular phylogeny estimate of the Australian cycad weevils, allowing an assessment of their monophyly, placement and relationships. Divergence dating suggests that the Australian cycad weevils originated from the Late Oligocene to the Middle Miocene and that the main radiation of the cycad-pollinating groups occurred from the Middle to the Late Miocene, which is congruent with the diversification of the Australian cycads, thus refuting any notion of an ancient ciophilous system in Australia. Taxonomic studies reveal the existence of 19 Australian cycad weevil species and that their associations with their hosts are mostly non-species-specific. Co-speciation analysis shows no extensive co-speciation events having occurred in the ciophilous system of Australian cycads. The distribution pattern suggests that geographical factors, rather than diversifying coevolution, constitute the overriding process shaping the Australian cycad weevil diversity. The synchronous radiation of cycads and weevil pollinators is suggested to be a result of the post-Oligocene diversification common in Australian organisms.

Large-Scale Identification and Characterization Analysis of VQ Family Genes in Plants, Especially Gymnosperms.

VQ motif-containing (VQ) proteins are a class of transcription regulatory cofactors widely present in plants, playing crucial roles in growth and development, stress response, and defense. Although there have been some reports on the member identification and functional research of VQ genes in some plants, there is still a lack of large-scale identification and clear graphical presentation of their basic characterization information to help us to better understand this family. Especially in gymnosperms, the VQ family genes and their evolutionary relationships have not yet been reported. In this study, we systematically identified 2469 VQ genes from 56 plant species, including bryophytes, gymnosperms, and angiosperms, and analyzed their molecular and evolutionary features. We found that amino acids are only highly conserved in the VQ domain, while other positions are relatively variable; most VQ genes encode relatively small proteins and do not have introns. The GC content in Poaceae plants is the highest (up to 70%); these VQ proteins can be divided into nine subgroups. In particular, we analyzed the molecular characteristics, chromosome distribution, duplication events, and expression levels of VQ genes in three gymnosperms: Ginkgo biloba, Taxus chinensis, and Pinus tabuliformis. In gymnosperms, VQ genes are classified into 11 groups, with highly similar motifs in each group; most VQ proteins have less than 300 amino acids and are predicted to be located in nucleus. Tandem duplication is an important driving force for the expansion of the VQ gene family, and the evolutionary processes of most VQ genes and duplication events are relatively independent; some candidate VQ genes are preliminarily screened, and they are likely to be involved in plant growth and stress and defense responses. These results provide detailed information and powerful references for further understanding and utilizing the VQ family genes in various plants.

Proline concentrations in seedlings of woody plants change with drought stress duration and are mediated by seed characteristics: a meta-analysis.

Proline accumulation represents one of mechanisms used by plants to prevent the adverse consequences of water stress. The effects of increased proline levels in response to drought differ among species. Trees are exposed to the long-term effects of climate change. The reproductive success of species in a specific environment depends on the functional trait of tree seeds. We conducted a meta-analysis to evaluate the effects of drought stress on the proline concentrations in seedling leaf tissues of woody plant species and their relationships to drought duration, seed mass, seed category and coniferous/deciduous classification. Drought duration exhibited a nonlinear effect on proline accumulations. The drought effect on proline accumulations is greater for deciduous than for coniferous species and is higher for orthodox seed species than for recalcitrant. The seedlings of large-seeded species showed greater effect sizes than those of small-seeded species. Our results suggest that there is an optimum level at which proline accumulations under the influence of drought are the highest. A link between seed functional traits, as well as the coniferous/deciduous classification, and proline concentrations in tree seedlings during water stress were determined for the first time. Proline may help to identify high-quality seeds of trees used for reforestation.

Does a gender of Welwitschia mirabilis plants influence their photosynthetic activity?

Welwitschia mirabilis Hook.f. (Welwitschiaceae, Gnetales) is a gymnosperm plant unique in its habit with an isolated taxonomic position. This species is dioecious, but no studies of its photosynthetic activity were conducted with examination of differences among male and female plants. To fill this gap, the day and night photosynthetic activity of male and female specimens of Welwitschia mirabilis cultivated in the botanical garden was studied in controlled conditions. Photosynthetic activity was studied using net photosynthetic rate (PN), stomatal conductance (gs) and intercellular CO2 concentration (Ci) parameters. Additionally, a normalized difference vegetation index (NDVI) was used to assess the condition among male and female plants in full sunlight. The studied Welwitschia plants revealed variability in photosynthetic activity both during the day and the night. The photosynthetic activity was low in the morning hours and higher in the afternoon. There is a difference in the photosynthetic activity during the night between sexes, being higher in female specimens. Stomatal density was evaluated separately for adaxial and abaxial leaf surfaces. Statistically significant differences in the stomatal density on abaxial and adaxial leaf surfaces were observed in both sexes, especially distinctive in female specimens. NDVI has revealed that there were weak differences between male and female plants.

Insights from multi-omics integration into seed germination of Taxus chinensis var mairei.

The transition from deep dormancy to seed germination is essential for the life cycle of plants, but how this process occurs in the gymnosperm Chinese yew (Taxus chinensis var mairei), the natural source of the anticancer drug paclitaxel, remains unclear. Herein, we analyse the transcriptome, proteome, spatial metabolome, and spatial lipidome of the Chinese yew and present the multi-omics profiles of dormant and germinating seeds. Our results show that abscisic acid and gibberellic acid 12 homoeostasis is closely associated with gene transcription and protein translation, and the balance between these phytohormones thereby determines if seeds remain dormant or germinate. We find that an energy supply of carbohydrates from glycolysis and the TCA cycle feed into the pentose phosphate pathway during seed germination, and energy supplied from lipids are mainly derived from the lipolysis of triacylglycerols. Using mass spectrometry imaging, we demonstrate that the spatial distribution of plant hormones and phospholipids has a remarkable influence on embryo development. We also provide an atlas of the spatial distribution of paclitaxel C in Chinese yew seeds for the first time. The data from this study enable exploration of the germination mechanism of Chinese yew seeds across several omics levels.

Origin and Evolution of Plant Long Terminal Repeat Retrotransposons with Additional Ribonuclease H.

Retroviruses originated from long terminal repeat retrotransposons (LTR-RTs) through several structural adaptations. One such modification was the arrangement of an additional ribonuclease H (aRH) domain next to native RH, followed by degradation and subfunctionalization of the latter. We previously showed that this retrovirus-like structure independently evolved in Tat LTR-RTs in flowering plants, proposing its origin from sequential rearrangements of ancestral Tat structures identified in lycophytes and conifers. However, most nonflowering plant genome assemblies were not available at that time, therefore masking the history of aRH acquisition by Tat and challenging our hypothesis. Here, we revisited Tat's evolution scenario upon the aRH acquisition by covering most of the extant plant phyla. We show that Tat evolved and obtained aRH in an ancestor of land plants. Importantly, we found the retrovirus-like structure in clubmosses, hornworts, ferns, and gymnosperms, suggesting its ancient origin, broad propagation, and yet-to-be-understood benefit for the LTR-RTs' adaptation.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing.

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.