Analysis of auxin responses in the fern Ceratopteris richardii identifies the developmental phase as a major determinant for response properties.

The auxin signaling molecule regulates a range of plant growth and developmental processes. The core transcriptional machinery responsible for auxin-mediated responses is conserved across all land plants. Genetic, physiological and molecular exploration in bryophyte and angiosperm model species have shown both qualitative and quantitative differences in auxin responses. Given the highly divergent ontogeny of the dominant gametophyte (bryophytes) and sporophyte (angiosperms) generations, however, it is unclear whether such differences derive from distinct phylogeny or ontogeny. Here, we address this question by comparing a range of physiological, developmental and molecular responses to auxin in both generations of the model fern Ceratopteris richardii. We find that auxin response in Ceratopteris gametophytes closely resembles that of a thalloid bryophyte, whereas the sporophyte mimics auxin response in flowering plants. This resemblance manifests both at the phenotypic and transcriptional levels. Furthermore, we show that disrupting auxin transport can lead to ectopic sporophyte induction on the gametophyte, suggesting a role for auxin in the alternation of generations. Our study thus identifies developmental phase, rather than phylogeny, as a major determinant of auxin response properties in land plants.

A conserved GRAS-domain transcriptional regulator links meristem indeterminacy to sex determination in Ceratopteris gametophytes.

Most land plants alternate between generations of sexual gametophytes and asexual sporophytes. Unlike seed plants, fern gametophytes are free living and grow independently of their sporophytes. In homosporous ferns such as Ceratopteris, gametophytes derived from genetically identical spores exhibit sexual dimorphism, developing as either males or hermaphrodites. Males lack meristems and promote cell differentiation into sperm-producing antheridia. In contrast, hermaphrodites initiate multicellular meristems that stay undifferentiated, sustain cell division and prothallus expansion, and drive the formation of egg-producing archegonia. Once initiating the meristem, hermaphrodites secrete the pheromone antheridiogen, which triggers neighboring slower-growing gametophytes to develop as males, while the hermaphrodites themselves remain insensitive to antheridiogen. This strategy promotes outcrossing and prevents all individuals in the colony from becoming males. This study reveals that an evolutionarily conserved GRAS-domain transcriptional regulator (CrHAM), directly repressed by Ceratopteris microRNA171 (CrmiR171), promotes meristem development in Ceratopteris gametophytes and determines the male-to-hermaphrodite ratio in the colony. CrHAM preferentially accumulates within the meristems of hermaphrodites but is excluded from differentiated antheridia. CrHAM sustains meristem proliferation and cell division through conserved hormone pathways. In the meantime, CrHAM inhibits the antheridiogen-induced conversion of hermaphrodites to males by suppressing the male program expression and preventing meristem cells from differentiating into sperm-producing antheridia. This finding establishes a connection between meristem indeterminacy and sex determination in ferns, suggesting both conserved and diversified roles of meristem regulators in land plants.

Control of leaf development in the water fern Ceratopteris richardii by the auxin efflux transporter CrPINMa in the CRISPR/Cas9 analysis.

BACKGROUND: PIN-FORMED genes (PINs) are crucial in plant development as they determine the directionality of auxin flow. They are present in almost all land plants and even in green algae. However, their role in fern development has not yet been determined. This study aims to investigate the function of CrPINMa in the quasi-model water fern Ceratopteris richardii. RESULTS: CrPINMa possessed a long central hydrophilic loop and characteristic motifs within it, which indicated that it belonged to the canonical rather than the non-canonical PINs. CrPINMa was positioned in the lineage leading to Arabidopsis PIN6 but not that to its PIN1, and it had undergone numerous gene duplications. CRISPR/Cas9 genome editing had been performed in ferns for the first time, producing diverse mutations including local frameshifts for CrPINMa. Plants possessing disrupted CrPINMa exhibited retarded leaf emergence and reduced leaf size though they could survive and reproduce at the same time. CrPINMa transcripts were distributed in the shoot apical meristem, leaf primordia and their vasculature. Finally, CrPINMa proteins were localized to the plasma membrane rather than other cell parts. CONCLUSIONS: CRISPR/Cas9 genome editing is feasible in ferns, and that PINs can play a role in fern leaf development.

The size diversity of the Pteridaceae family chloroplast genome is caused by overlong intergenic spacers.

BACKGROUND: While the size of chloroplast genomes (cpDNAs) is often influenced by the expansion and contraction of inverted repeat regions and the enrichment of repeats, it is the intergenic spacers (IGSs) that appear to play a pivotal role in determining the size of Pteridaceae cpDNAs. This provides an opportunity to delve into the evolution of chloroplast genomic structures of the Pteridaceae family. This study added five Pteridaceae species, comparing them with 36 published counterparts. RESULTS: Poor alignment in the non-coding regions of the Pteridaceae family was observed, and this was attributed to the widespread presence of overlong IGSs in Pteridaceae cpDNAs. These overlong IGSs were identified as a major factor influencing variations in cpDNA size. In comparison to non-expanded IGSs, overlong IGSs exhibited significantly higher GC content and were rich in repetitive sequences. Species divergence time estimations suggest that these overlong IGSs may have already existed during the early radiation of the Pteridaceae family. CONCLUSIONS: This study reveals new insights into the genetic variation, evolutionary history, and dynamic changes in the cpDNA structure of the Pteridaceae family, providing a fundamental resource for further exploring its evolutionary research.

Adaptive evolution and co-evolution of chloroplast genomes in Pteridaceae species occupying different habitats: overlapping residues are always highly mutated.

BACKGROUND: The evolution of protein residues depends on the mutation rates of their encoding nucleotides, but it may also be affected by co-evolution with other residues. Chloroplasts function as environmental sensors, transforming fluctuating environmental signals into different physiological responses. We reasoned that habitat diversity may affect their rate and mode of evolution, which might be evidenced in the chloroplast genome. The Pteridaceae family of ferns occupy an unusually broad range of ecological niches, which provides an ideal system for analysis. RESULTS: We conducted adaptive evolution and intra-molecular co-evolution analyses of Pteridaceae chloroplast DNAs (cpDNAs). The results indicate that the residues undergoing adaptive evolution and co-evolution were mostly independent, with only a few residues being simultaneously involved in both processes, and these overlapping residues tend to exhibit high mutations. Additionally, our data showed that Pteridaceae chloroplast genes are under purifying selection. Regardless of whether we grouped species by lineage (which corresponded with ecological niches), we determined that positively selected residues mainly target photosynthetic genes. CONCLUSIONS: Our work provides evidence for the adaptive evolution of Pteridaceae cpDNAs, especially photosynthetic genes, to different habitats and sheds light on the adaptive evolution and co-evolution of proteins.

N-terminal region is required for functions of the HAM family member.

Shoot meristems contain stem cells, and they sustain growth and development of the above-ground tissues in land plants. The HAIRY MERISTEM (HAM) family genes, encoding GRAS-domain transcriptional regulators, play essential roles in the control of shoot meristem development and stem cell homeostasis in several flowering plants. Similar to other GRAS proteins, the C-terminal regions of HAM family proteins across land plants are conserved, containing signature motifs that define the GRAS domain. In contrast, the N-terminal regions of HAM family proteins display substantial divergence in sequence and length. Whether the variable and divergent N-termini are required for the conserved functions of HAM proteins is unknown. Our recent work showed that CrHAM - the HAM homolog in the fern Ceratopteris richardii was able to replace the role of type-II HAM genes in Arabidopsis, maintaining established shoot apical meristems and promoting the initiation of new stem cell niches. Here, we provide additional information and show that CrHAM contains a much longer N-terminal region compared to Arabidopsis HAM proteins, which is conserved among different fern HAM homologs. The deletion of this region largely compromises the ability of CrHAM to replace the function of Arabidopsis HAM proteins in shoot meristems. These new data together with previous results suggest that, although lacking the sequence conservation among HAM homologs from different plant lineages, the N-termini are important for the conserved functions of HAM family proteins.

A De Novo Transcriptome Assembly of Ceratopteris richardii Provides Insights into the Evolutionary Dynamics of Complex Gene Families in Land Plants.

As the closest extant sister group to seed plants, ferns are an important reference point to study the origin and evolution of plant genes and traits. One bottleneck to the use of ferns in phylogenetic and genetic studies is the fact that genome-level sequence information of this group is limited, due to the extreme genome sizes of most ferns. Ceratopteris richardii (hereafter Ceratopteris) has been widely used as a model system for ferns. In this study, we generated a transcriptome of Ceratopteris, through the de novo assembly of the RNA-seq data from 17 sequencing libraries that are derived from two sexual types of gametophytes and five different sporophyte tissues. The Ceratopteris transcriptome, together with 38 genomes and transcriptomes from other species across the Viridiplantae, were used to uncover the evolutionary dynamics of orthogroups (predicted gene families using OrthoFinder) within the euphyllophytes and identify proteins associated with the major shifts in plant morphology and physiology that occurred in the last common ancestors of euphyllophytes, ferns, and seed plants. Furthermore, this resource was used to identify and classify the GRAS domain transcriptional regulators of many developmental processes in plants. Through the phylogenetic analysis within each of the 15 GRAS orthogroups, we uncovered which GRAS family members are conserved or have diversified in ferns and seed plants. Taken together, the transcriptome database and analyses reported here provide an important platform for exploring the evolution of gene families in land plants and for studying gene function in seed-free vascular plants.

Cryptic diversity in the model fern genus Ceratopteris (Pteridaceae).

Cryptic species are present throughout the tree of life. They are especially prevalent in ferns, because of processes such hybridization, polyploidy, and reticulate evolution. In addition, the simple morphology of ferns limits phenotypic variation and makes it difficult to detect cryptic species. The model fern genus Ceratopteris has long been suspected to harbor cryptic diversity, in particular within the highly polymorphic C. thalictroides. Yet no studies have included samples from throughout its pan-tropical range or utilized genomic sequencing, making it difficult to assess the full extent of cryptic variation within this genus. Here, we present the first multilocus phylogeny of the genus using reduced representation genomic sequencing (RADseq) and examine population structure, phylogenetic relationships, and ploidy level variation. We recover similar species relationships found in previous studies, find support for the cryptic species C. gaudichaudii as genetically distinct, and identify novel genomic variation within two of the mostly broadly distributed species in the genus, C. thalictroides and C. cornuta. Finally, we detail the utility of our approach for working on cryptic, reticulate groups of ferns. Specifically, it does not require a reference genome, of which there are very few available for ferns. RADseq is a cost-effective way to work with study groups lacking genomic resources, and to obtain the thousands of nuclear markers needed to untangle species complexes.

Molecular Evolution of Auxin-Mediated Root Initiation in Plants.

The root originated independently in euphyllophytes (ferns and seed plants) and lycophytes; however, the molecular evolutionary route of root initiation remains elusive. By analyses of the fern Ceratopteris richardii and seed plants, here we show that the molecular pathway involving auxin, intermediate-clade WUSCHEL-RELATED HOMEOBOX (IC-WOX) genes, and WUSCHEL-clade WOX (WC-WOX) genes could be conserved in root initiation. We propose that the "auxin>IC-WOX>WC-WOX" module in root initiation might have arisen in the common ancestor of euphyllophytes during the second origin of roots, and that this module has further developed during the evolution of different root types in ferns and seed plants.

An ontogenetic framework for functional studies in the model fern Ceratopteris richardii.

As the sister group to seed plants, ferns are a phylogenetically informative lineage. Functional studies in representatives of the fern lineage are helping bridge the knowledge gap in developmental mechanisms between angiosperms and non-vascular plants. The fern life cycle has the advantage of combining a sizable free-living haploid gametophyte, more amenable for developmental studies than the reduced seed plant gametophyte, with an indeterminate and complex diploid sporophyte. Ceratopteris richardii has long been proposed as a model fern and has recently become tractable due to stable transgenesis and increasing genomic resources, allowing researchers to test explicit questions about gene function in a fern for the first time. As with any model system, a detailed understanding of wild-type morphology and a staged ontogeny are indispensable for the characterization of mutant phenotypes resulting from genetic manipulations. Therefore, the goal of this study is to provide a unified reference ontogeny for this emerging model fern as a tool for comparative evolutionary and developmental studies. It complements earlier research by filling gaps in major stages of development of the haploid gametophyte and diploid sporophyte generations, and provides additional descriptions of the shoot apical meristem and early leaf development. This resource is meant to facilitate not only studies of candidate genes within C. richardii, but also broader ontogenetic comparisons to other model plants.

A fern WUSCHEL-RELATED HOMEOBOX gene functions in both gametophyte and sporophyte generations.

BACKGROUND: Post-embryonic growth of land plants originates from meristems. Genetic networks in meristems maintain the stem cells and direct acquisition of cell fates. WUSCHEL-RELATED HOMEOBOX (WOX) transcription factors involved in meristem networks have only been functionally characterized in two evolutionarily distant taxa, mosses and seed plants. This report characterizes a WOX gene in a fern, which is located phylogenetically between the two taxa. RESULTS: CrWOXB transcripts were detected in proliferating tissues, including gametophyte and sporophyte meristems of Ceratopteris richardii. In addition, CrWOXB is expressed in archegonia but not the antheridia of gametophytes. Suppression of CrWOXB expression in wild-type RN3 plants by RNAi produced abnormal morphologies of gametophytes and sporophytes. The gametophytes of RNAi lines produced fewer cells, and fewer female gametes compared to wild-type. In the sporophyte generation, RNAi lines produced fewer leaves, pinnae, roots and lateral roots compared to wild-type sporophytes. CONCLUSIONS: Our results suggest that CrWOXB functions to promote cell divisions and organ development in the gametophyte and sporophyte generations, respectively. CrWOXB is the first intermediate-clade WOX gene shown to function in both generations in land plants.

Proteomic analysis of gametophytic sex expression in the fern Ceratopteris thalictroides.

Ceratopteris thalictroides, a model fern, has two kinds of gametophytes with different sex expression: male and hermaphrodite. Hermaphroditic gametophytes have one or several archegonia beneath the growing point and a few antheridia at the base or margin. Male gametophytes show a spoon-like shape with much longer than the width and produce many antheridia at the margin and surface. The results of chlorophyll fluorescence detection showed that the photochemical efficiency of hermaphrodites was higher than that of males. By using two-dimensional electrophoresis and mass spectrometry, the differentially abundant proteins in hermaphroditic and male gametophytes were identified. A total of 1136 ± 55 protein spots were detected in Coomassie-stained gels of proteins from hermaphroditic gametophytes, and 1130 ± 65 spots were detected in gels of proteins from male gametophytes. After annotation, 33 spots representing differentially abundant proteins were identified. Among these, proteins involved in photosynthesis and chaperone proteins were over-represented in hermaphrodites, whereas several proteins involved in metabolism were increased in male gametophytes in order to maintain their development under relatively nutritionally deficient conditions. Furthermore, the differentially abundant cytoskeletal proteins detected in this study, such as centrin and actin, may be involved in the formation of sexual organs and are directly related to sex expression. These differentially abundant proteins are important for maintaining the development of gametophytes of different sexes in C. thalictroides.

Low-copy nuclear sequence data confirm complex patterns of farina evolution in notholaenid ferns (Pteridaceae).

Notholaenids are an unusual group of ferns that have adapted to, and diversified within, the deserts of Mexico and the southwestern United States. With approximately 40 species, this group is noted for being desiccation-tolerant and having "farina"-powdery exudates of lipophilic flavonoid aglycones-that occur on both the gametophytic and sporophytic phases of their life cycle. The most recent circumscription of notholaenids based on plastid markers surprisingly suggests that several morphological characters, including the expression of farina, are homoplasious. In a striking case of convergence, Notholaena standleyi appears to be distantly related to core Notholaena, with several taxa not before associated with Notholaena nested between them. Such conflicts can be due to morphological homoplasy resulting from adaptive convergence or, alternatively, the plastid phylogeny itself might be misleading, diverging from the true species tree due to incomplete lineage sorting, hybridization, or other factors. In this study, we present a species phylogeny for notholaenid ferns, using four low-copy nuclear loci and concatenated data from three plastid loci. A total of 61 individuals (49 notholaenids and 12 outgroup taxa) were sampled, including 31 out of 37 recognized notholaenid species. The homeologous/allelic nuclear sequences were retrieved using PacBio sequencing and the PURC bioinformatics pipeline. Each dataset was first analyzed individually using maximum likelihood and Bayesian inference, and the species phylogeny was inferred using *BEAST. Although we observed several incongruences between the nuclear and plastid phylogenies, our principal results are broadly congruent with previous inferences based on plastid data. By mapping the presence of farina and their biochemical constitutions on our consensus phylogenetic tree, we confirmed that the characters are indeed homoplastic and have complex evolutionary histories. Hybridization among recognized species of the notholaenid clade appears to be relatively rare compared to that observed in other well-studied fern genera.

Meristem development and activity in gametophytes of the model fern, Ceratopteris richardii.

Ceratopteris richardii is a model fern species widely used to analyze various developmental processes and their regulation in gametophytes. The form of mature C. richardii gametophytes depends on the activity of the marginal meristem, but knowledge on meristem formation and structure is limited. Therefore, we analyzed cellular events accompanying the development of gametophytes using cell lineage and proliferation analyses to explain the establishment and functioning of the marginal meristem. We show that: i) gametophytes are devoid of the apical initial cell or the apical cell-based meristem in the early developmental stages; ii) the cells that are predestined to form the marginal meristem divide according to a stable pattern; iii) only one transient initial cell is present in the marginal meristem, and the selection of a new functioning initial cell is related to a stable sequence of its divisions. Our results contribute to a better understanding of the developmental events underlying gametophyte growth and marginal meristem functioning in Ceratopteris. The principles, which were established in this study and enabled the identification of functioning initial cells, can be applied to analyze genetic and/or physiological mechanism(s) governing meristem maintenance in vascular plants, both in developmental and evolutionary contexts.

Mobile Elements Shape Plastome Evolution in Ferns.

Plastid genomes display remarkable organizational stability over evolutionary time. From green algae to angiosperms, most plastid genomes are largely collinear, with only a few cases of inversion, gene loss, or, in extremely rare cases, gene addition. These plastome insertions are mostly clade-specific and are typically of nuclear or mitochondrial origin. Here, we expand on these findings and present the first family-level survey of plastome evolution in ferns, revealing a novel suite of dynamic mobile elements. Comparative plastome analyses of the Pteridaceae expose several mobile open reading frames that vary in sequence length, insertion site, and configuration among sampled taxa. Even between close relatives, the presence and location of these elements is widely variable when viewed in a phylogenetic context. We characterize these elements and refer to them collectively as Mobile Open Reading Frames in Fern Organelles (MORFFO). We further note that the presence of MORFFO is not restricted to Pteridaceae, but is found across ferns and other plant clades. MORFFO elements are regularly associated with inversions, intergenic expansions, and changes to the inverted repeats. They likewise appear to be present in mitochondrial and nuclear genomes of ferns, indicating that they can move between genomic compartments with relative ease. The origins and functions of these mobile elements are unknown, but MORFFO appears to be a major driver of structural genome evolution in the plastomes of ferns, and possibly other groups of plants.

Phylogeny and systematics of the brake fern genus Pteris (Pteridaceae) based on molecular (plastid and nuclear) and morphological evidence.

The brake fern genus Pteris belongs to Pteridaceae subfamily Pteridoideae. It is one of the largest fern genera and has been estimated to contain 200-250 species distributed on all continents except Antarctica. Previous studies were either based on plastid data only or based on both plastid and nuclear data but the sampling was small. In addition, an infrageneric classification of Pteris based on morphological and molecular evidence has not been available yet. In the present study, based on molecular data of eight plastid markers and one nuclear marker (gapCp) of 256 accessions representing ca. 178 species of Pteris, we reconstruct a global phylogeny of Pteris. The 15 major clades identified earlier are recovered here and we further identified a new major clade. Our nuclear phylogeny recovered 11 of these 16 major clades, seven of which are strongly supported. The inclusion of Schizostege in Pteris is confirmed for the first time. Based on the newly reconstructed phylogeny and evidence from morphology, distribution and/or ecology, we classify Pteris into three subgenera: P. subg. Pteris, P. subg. Campteria, and P. subg. Platyzoma. The former two are further divided into three and 12 sections, respectively.

Phylogeny of the fern subfamily Pteridoideae (Pteridaceae; Pteridophyta), with the description of a new genus: Gastoniella.

As the second most genera-rich fern family, Pteridaceae contain more than 1000 species contributing to ca. 10% of extant leptosporangiate fern diversity. The subfamily Pteridoideae is one of the five subfamilies often recognized. The circumscription of Pteridoideae has not been clear. A large number of species have not yet been included in any molecular analyses before. In this study, DNA sequences of six plastid loci of 154 accessions representing ca. 87 species in 14 genera of Pteridaceae subfam. Pteridoideae and four accessions representing two species in subfam. Parkerioideae and one species of subfam. Adiantoideae as outgroups were used to infer a phylogeny using maximum likelihood and maximum parsimony. Our analyses show that (1) Pteridoideae is monophyletic and the newly defined subfamily is composed of 14 genera including a newly described genus; (2) Pteridoideae is resolved into four strongly supported monophyletic clades: the Pteris clade, the Actiniopteris+Onychium clade, the JAPSTT clade, and the GAPCC clade, these being supported by not only molecular data but also morphological features and distribution information; (3) Onychium is confirmed as monophyletic and accessions of Onychium are resolved into two strongly supported clades, the O. cryptogrammoides clade and the O. siliculosum clade; and (4) Accessions of the traditionally defined Anogramma are resolved as paraphyletic in relation to Cerosora, Cosentinica, and Pityrogramma. Three species traditionally treated in Anogramma are in fact more closely related to Cerosora and Pityrogramma than they are to Anogramma. Gastoniella Li Bing Zhang & Liang Zhang, gen. nov. is described to accommodate these species and three new combinations are provided. Three currently known species of Gastoniella are distributed in the Ascension Island in South Atlantic Ocean, central Mexico, and tropical America, respectively. The new genus is distinct from Anogramma s.s. in having ultimate segments linear not obviously broadening toward the upper portion.

Not only in the temperate zone: independent gametophytes of two vittarioid ferns (Pteridaceae, Polypodiales) in East Asian subtropics.

Independent gametophyte ferns are unique among vascular plants because they are sporophyteless and reproduce asexually to maintain their populations in the gametophyte generation. Such ferns had been primarily discovered in temperate zone, and usually hypothesized with (sub)tropical origins and subsequent extinction of sporophyte due to climate change during glaciations. Presumably, independent fern gametophytes are unlikely to be distributed in tropics and subtropics because of relatively stable climates which are less affected by glaciations. Nonetheless, the current study presents cases of two independent gametophyte fern species in subtropic East Asia. In this study, we applied plastid DNA sequences (trnL-L-F and matK + ndhF + chlL datasets) and comprehensive sampling (~80%) of congeneric species for molecular identification and divergence time estimation of these independent fern gametophytes. The two independent gametophyte ferns were found belonging to genus Haplopteris (vittarioids, Pteridaceae) and no genetic identical sporophyte species in East Asia. For one species, divergence times between its populations imply recent oversea dispersal(s) by spores occurred during Pleistocene. By examining their ex situ and in situ fertility, prezygotic sterility was found in these two Haplopteris, in which gametangia were not or very seldom observed, and this prezygotic sterility might attribute to their lacks of functional sporophytes. Our field observation and survey on their habitats suggest microhabitat conditions might attribute to this prezygotic sterility. These findings point to consideration of whether recent climate change during the Pleistocene glaciation resulted in ecophysiological maladaptation of non-temperate independent gametophyte ferns. In addition, we provided a new definition to classify fern gametophyte independences at the population level. We expect that continued investigations into tropical and subtropical fern gametophyte floras will further illustrate the biogeographic significance of non-temperate fern gametophyte independence.

Plasma membrane-anchored chloroplasts are necessary for the gravisensing system of Ceratopteris richardii prothalli.

The prothalli of the fern Ceratopteris richardii exhibit negative gravitropism when grown in darkness. However, no sedimentable organelles or substances have been detected in the prothallial cells, suggesting that a non-sedimentable gravisensor exists. We investigated whether chloroplasts are involved in the gravisensing system of C. richardii prothalli. We used a clumped-chloroplast mutant, clumped chloroplast 1 (cp1), in which the chloroplasts are detached from the plasma membrane and clustered around the nucleus likely because of a partial deletion in the KINESIN-LIKE PROTEIN FOR ACTIN-BASED CHLOROPLAST MOVEMENT 1 gene. The cp1 mutation resulted in prothalli that had a significantly diminished gravitropic response, while the phototropic response occurred normally. These results suggest that plasma membrane-anchored chloroplasts in prothallial cells function as one of the gravisensors in C. richardii prothalli.

Convergent Evolution of Fern-Specific Mitochondrial Group II Intron atp1i361g2 and Its Ancient Source Paralogue rps3i249g2 and Independent Losses of Intron and RNA Editing among Pteridaceae.

Mitochondrial intron patterns are highly divergent between the major land plant clades. An intron in the atp1 gene, atp1i361g2, is an example for a group II intron specific to monilophytes (ferns). Here, we report that atp1i361g2 is lost independently at least 4 times in the fern family Pteridaceae. Such plant organelle intron losses have previously been found to be accompanied by loss of RNA editing sites in the flanking exon regions as a consequence of genomic recombination of mature cDNA. Instead, we now observe that RNA editing events in both directions of pyrimidine exchange (C-to-U and U-to-C) are retained in atp1 exons after loss of the intron in Pteris argyraea/biaurita and in Actiniopteris and Onychium We find that atp1i361g2 has significant similarity with intron rps3i249g2 present in lycophytes and gymnosperms, which we now also find highly conserved in ferns. We conclude that atp1i361g2 may have originated from the more ancestral rps3i249g2 paralogue by a reverse splicing copy event early in the evolution of monilophytes. Secondary structure elements of the two introns, most characteristically their domains III, show strikingly convergent evolution in the monilophytes. Moreover, the intron paralogue rps3i249g2 reveals relaxed evolution in taxa where the atp1i361g2 paralogue is lost. Our findings may reflect convergent evolution of the two related mitochondrial introns exerted by co-evolution with an intron-binding protein simultaneously acting on the two paralogues.