De Novo Transcriptome Assembly of Two Microsorum Fern Species Identifies Enzymes Required for Two Upstream Pathways of Phytoecdysteroids.

Microsorum species produce a high amount of phytoecdysteroids (PEs), which are widely used in traditional medicine in the Pacific islands. The PEs in two different Microsorum species, M. punctatum (MP) and M. scolopendria (MS), were examined using high-performance liquid chromatography (HPLC). In particular, MS produces a high amount of 20-hydroxyecdysone, which is the main active compound in PEs. To identify genes for PE biosynthesis, we generated reference transcriptomes from sterile frond tissues using the NovaSeq 6000 system. De novo transcriptome assembly after deleting contaminants resulted in 57,252 and 54,618 clean transcripts for MP and MS, respectively. The clean Microsorum transcripts for each species were annotated according to gene ontology terms, UniProt pathways, and the clusters of the orthologous group protein database using the MEGAN6 and Sma3s programs. In total, 1852 and 1980 transcription factors were identified for MP and MS, respectively. We obtained transcripts encoding for 38 and 32 enzymes for MP and MS, respectively, potentially involved in mevalonate and sterol biosynthetic pathways, which produce precursors for PE biosynthesis. Phylogenetic analyses revealed many redundant and unique enzymes between the two species. Overall, this study provides two Microsorum reference transcriptomes that might be useful for further studies regarding PE biosynthesis in Microsorum species.

A flavin-dependent monooxygenase catalyzes the initial step in cyanogenic glycoside synthesis in ferns.

Cyanogenic glycosides form part of a binary plant defense system that, upon catabolism, detonates a toxic hydrogen cyanide bomb. In seed plants, the initial step of cyanogenic glycoside biosynthesis-the conversion of an amino acid to the corresponding aldoxime-is catalyzed by a cytochrome P450 from the CYP79 family. An evolutionary conundrum arises, as no CYP79s have been identified in ferns, despite cyanogenic glycoside occurrence in several fern species. Here, we report that a flavin-dependent monooxygenase (fern oxime synthase; FOS1), catalyzes the first step of cyanogenic glycoside biosynthesis in two fern species (Phlebodium aureum and Pteridium aquilinum), demonstrating convergent evolution of biosynthesis across the plant kingdom. The FOS1 sequence from the two species is near identical (98%), despite diversifying 140 MYA. Recombinant FOS1 was isolated as a catalytic active dimer, and in planta, catalyzes formation of an N-hydroxylated primary amino acid; a class of metabolite not previously observed in plants.

Isolation and functional characterization of four microbial type terpene synthases from ferns.

Typical plant terpene synthases (TPSs) are responsible for the production of terpenes, a major class of plant secondary metabolites. However, various nonseed plants also harbor genes encoding microbial terpene synthase-like (MTPSL) enzymes. Here, a scan of 31 ferns transcriptomes revealed 40 sequences putatively encoding MTPSLs. Two groups of sequences were recognized based on the key conserved motifs. Four representative genes were isolated from each of the four species Adiantum capillus-veneris, Cyclosorus parasiticus, Drynaria bonii and Microlepia platyphylla. Following their heterologous expression in E. coli, the recombinant proteins were tested for monoterpene synthase and sesquiterpene synthase activity. These enzymatic products were typical monoterpenes and sesquiterpenes that have been previous shown to be generated by classical plant TPSs when provided with GPP and FPP as substrates. Subcellular localization experiments in the leaf epidermis of Nicotiana benthamiana and onion (Allium cepa) inner epidermal cells indicated that AcMTPSL1 and DbMTPSL were deposited in both the cytoplasm and nucleus, whereas CpMTPSL1 and MpMTPSL were localized in the cytoplasm, chloroplasts and nucleus. AcMTPSL1 was up-regulated in plants exposed to methyl jasmonate treatment, suggesting a role for this gene in host defense. This study provides more information about the catalytic function of MTPSLs in nonseed plants and for the first time, the subcellular localization of MTPSLs was experimentally characterized.

Bioinformatics Analysis of Plant Cell Wall Evolution.

In the past hundreds of millions of years, from green algae to land plants, cell walls have developed into a highly complex structure that is essential for plant growth and survival. Plant cell wall diversity and evolution can be directly investigated by chemically profiling polysaccharides and lignins in the cell walls of diverse plants and algae. With the increasingly low cost and high throughput of DNA sequencing technologies, cell wall evolution can also be studied by bioinformatics analysis of the occurrence of cell wall synthesis-related enzymes in the genomes and transcriptomes of different species. This chapter presents a bioinformatics workflow running on a Linux platform to process genomic data for such gene occurrence analysis. As a case study, cellulose synthase (CesA) and CesA-like (Csl) protein families are mined for in two newly sequenced organisms: the charophyte green alga Klebsormidium flaccidum (renamed as Klebsormidium nitens) and the fern Lygodium japonicum.

Identification and evolutionary analysis of chalcone isomerase-fold proteins in ferns.

The distribution of type I and II chalcone isomerases (CHIs) in plants is highly family specific. We have previously reported that ancient land plants, such as the liverworts and Selaginella moellendorffii, harbor type II CHIs. To better understand the function and evolution of CHI-fold proteins, transcriptomic data obtained from 52 pteridophyte species were subjected to sequence alignment and phylogenetic analysis. The residues determining type I/II CHI identity in the pteridophyte CHIs were identical to those of type I CHIs. The enzymatic characterization of a sample of 24 CHIs, representing all the key pteridophyte lineages, demonstrated that 19 of them were type I enzymes and that five exhibited some type II activity due to an amino acid mutation. Two pteridophyte chalcone synthases (CHSs) were also characterized, and a type IV CHI (CHIL) was demonstrated to interact physically with CHSs and CHI, and to increase CHS activity by decreasing derailment products, thus enhancing flavonoid production. These findings suggest that the emergence of type I CHIs may have coincided with the divergence of the pteridophytes. This study deepens our understanding of the molecular mechanism of CHIL as an enhancer in the flavonoid biosynthesis pathway.

Insecticidal fern protein Tma12 is possibly a lytic polysaccharide monooxygenase.

MAIN CONCLUSION: Amino acid sequence and crystal structure analyses of Tma12, an insecticidal protein isolated from the fern Tectaria macrodonta, identify it as a carbohydrate-binding protein belonging to the AA10 family of lytic polysaccharide monooxygenases, and provide the first evidence of AA10 proteins in plants. Tma12, isolated from the fern Tectaria macrodonta, is a next-generation insecticidal protein. Transgenic cotton expressing Tma12 exhibits resistance against whitefly and viral diseases. Beside its insecticidal property, the structure and function of Tma12 are unknown. This limits understanding of the insecticidal mechanism of the protein and targeted improvement in its efficacy. Here we report the amino acid sequence analysis and the crystal structure of Tma12, suggesting that it is possibly a lytic polysaccharide monooxygenase (LPMO) of the AA10 family. Amino acid sequence of Tma12 shows 45% identity with a cellulolytic LPMO of Streptomyces coelicolor. The crystal structure of Tma12, obtained at 2.2 Å resolution, possesses all the major structural characteristics of AA10 LPMOs. A H2O2-based enzymatic assay also supports this finding. It is the first report of the occurrence of LPMO-like protein in a plant. The two facts that Tma12 possesses insecticidal activity and shows structural similarity with LPMOs collectively advocate exploration of microbial LPMOs for insecticidal potential.

Anti-melanization effects and inhibitory kinetics of tyrosinase of bird's nest fern (Asplenium australasicum) frond extracts on melanoma and human skin.

Some bioactive properties of p-coumaric acid and fucose-rich polysaccharide in skin health have been studied, including melanogenesis inhibition of the phenolic acid and growth inhibitory effects of the polysaccharide on melanoma. The dermatological benefits of bird's nest fern extracts (BNFE), containing both substantial fucose-rich polysaccharide and p-coumaric acid, like promoting collagen production and growth of fibroblast cell and further improving the elasticity and dryness of human skins have been demonstrated in our previous study. Besides, the anti-melanization effects of various BNFE on B16-F10 melanoma and human skin were first studied here. The promising extracts revealed that the main phenolic acid, p-coumaric acid, in BNFE resulted in suppression against tyrosinase activity from melanogenesis. The inhibitory kinetics on the diphenolase activity indicated that AE40 was a noncompetitive inhibitor of mushroom tyrosinase. On the other hand, the fucose-rich mucilage of BNFE showed pronouncedly suppressing effect on B16-F10 melanoma viability. Clinical trial was performed by recruiting 46 female volunteers and the results indicated that the lotions with 1% of BNFE was non-irritant and reduced effectively the pigmentation on human skin after 7-14 days of continuous application. It was suggested that the fucose-rich mucilage and p-coumaric acid in BNFE may have potential for nutricosmetics and phytotherapy applications as a natural hypopigmenting agent.

Separate and combined effects of glyphosate and copper on growth and antioxidative enzymes in Salvinia natans (L.) All.

The coexistence of glyphosate and copper is widely found in bodies of water and terrestrial ecosystems due to widespread application of herbicides and heavy metal. However, their joint ecotoxicological risks in aquatic environments remain unknown. The experiment investigated the individual and combined effects of glyphosate and copper on the growth and physiological response in Salvinia natans (L.) All. The results showed that their joint toxicity is related to concentration. Antagonistic effects were induced when plants were exposed to low concentrations of glyphosate and copper (≤1 + 0.2 mg l-1). Synergistic effects were elicited at higher doses (≥5 + 1 mg l-1). In addition, increased hydrogen peroxide levels indicated the occurrence of oxidative stress at individual or combined exposures. To cope with oxidative stress, S. natans can activate the antioxidant defense systems, including increased superoxide dismutase and changes in peroxidase, ascorbate peroxidase and catalase. High concentrations of combined pollution exceed the oxidative defense capabilities of plants, and therefore, malondialdehyde content increased significantly. Our results indicated that the ecotoxicity of glyphosate or copper may be exacerbated in aquatic environments and caused obvious damage to S. natans.

Evolutionary History of the Glycoside Hydrolase 3 (GH3) Family Based on the Sequenced Genomes of 48 Plants and Identification of Jasmonic Acid-Related GH3 Proteins in Solanum tuberosum.

Glycoside Hydrolase 3 (GH3) is a phytohormone-responsive family of proteins found in many plant species. These proteins contribute to the biological activity of indolacetic acid (IAA), jasmonic acid (JA), and salicylic acid (SA). They also affect plant growth and developmental processes as well as some types of stress. In this study, GH3 genes were identified in 48 plant species, including algae, mosses, ferns, gymnosperms, and angiosperms. No GH3 representative protein was found in algae, but we identified 4 genes in mosses, 19 in ferns, 7 in gymnosperms, and several in angiosperms. The results showed that GH3 proteins are mainly present in seed plants. Phylogenetic analysis of all GH3 proteins showed three separate clades. Group I was related to JA adenylation, group II was related to IAA adenylation, and group III was separated from group II, but its function was not clear. The structure of the GH3 proteins indicated highly conserved sequences in the plant kingdom. The analysis of JA adenylation in relation to gene expression of GH3 in potato (Solanum tuberosum) showed that StGH3.12 greatly responded to methyl jasmonate (MeJA) treatment. The expression levels of StGH3.1, StGH3.11, and StGH3.12 were higher in the potato flowers, and StGH3.11 expression was also higher in the stolon. Our research revealed the evolution of the GH3 family, which is useful for studying the precise function of GH3 proteins related to JA adenylation in S. tuberosum when the plants are developing and under biotic stress.

Production of Hydroxynitrile Lyase from Davallia tyermannii (DtHNL) in Komagataella phaffii and Its Immobilization as a CLEA to Generate a Robust Biocatalyst.

Hydroxynitrile lyase from the white rabbit's foot fern Davallia tyermannii (DtHNL) catalyzes the enantioselective synthesis of α-cyanohydrins, which are key building blocks for pharmaceutical and agrochemical industries. An efficient and competitive process necessitates the availability and robustness of the biocatalyst. Herein, the recombinant production of DtHNL1 in Komagataella phaffii, yielding approximately 900 000 U L-1 , is described. DtHNL1 constitutes approximately 80 % of the total protein content. The crude enzyme was immobilized. Crosslinked enzyme aggregates (CLEAs) resulted in significant enhancement of the biocatalyst stability under acidic conditions (activity retained after 168 h at pH 2.4). The DtHNL1-CLEA was employed for (R)-mandelonitrile synthesis (99 % conversion, 98 % enantiomeric excess) in a biphasic system, and evaluated for the synthesis of (R)-hydroxypivaldehyde cyanohydrin under reaction conditions that immediately inactivated non-immobilized DtHNL1. The results show the DtHNL1-CLEA to be a stable biocatalyst for the synthesis of enantiomerically pure cyanohydrins under acidic conditions.

The ALDH21 gene found in lower plants and some vascular plants codes for a NADP+ -dependent succinic semialdehyde dehydrogenase.

Lower plant species including some green algae, non-vascular plants (bryophytes) as well as the oldest vascular plants (lycopods) and ferns (monilophytes) possess a unique aldehyde dehydrogenase (ALDH) gene named ALDH21, which is upregulated during dehydration. However, the gene is absent in flowering plants. Here, we show that ALDH21 from the moss Physcomitrella patens codes for a tetrameric NADP+ -dependent succinic semialdehyde dehydrogenase (SSALDH), which converts succinic semialdehyde, an intermediate of the γ-aminobutyric acid (GABA) shunt pathway, into succinate in the cytosol. NAD+ is a very poor coenzyme for ALDH21 unlike for mitochondrial SSALDHs (ALDH5), which are the closest related ALDH members. Structural comparison between the apoform and the coenzyme complex reveal that NADP+ binding induces a conformational change of the loop carrying Arg-228, which seals the NADP+ in the coenzyme cavity via its 2'-phosphate and α-phosphate groups. The crystal structure with the bound product succinate shows that its carboxylate group establishes salt bridges with both Arg-121 and Arg-457, and a hydrogen bond with Tyr-296. While both arginine residues are pre-formed for substrate/product binding, Tyr-296 moves by more than 1 Å. Both R121A and R457A variants are almost inactive, demonstrating a key role of each arginine in catalysis. Our study implies that bryophytes but presumably also some green algae, lycopods and ferns, which carry both ALDH21 and ALDH5 genes, can oxidize SSAL to succinate in both cytosol and mitochondria, indicating a more diverse GABA shunt pathway compared with higher plants carrying only the mitochondrial ALDH5.

Enzyme discovery beyond homology: a unique hydroxynitrile lyase in the Bet v1 superfamily.

Homology and similarity based approaches are most widely used for the identification of new enzymes for biocatalysis. However, they are not suitable to find truly novel scaffolds with a desired function and this averts options and diversity. Hydroxynitrile lyases (HNLs) are an example of non-homologous isofunctional enzymes for the synthesis of chiral cyanohydrins. Due to their convergent evolution, finding new representatives is challenging. Here we show the discovery of unique HNL enzymes from the fern Davallia tyermannii by coalescence of transcriptomics, proteomics and enzymatic screening. It is the first protein with a Bet v1-like protein fold exhibiting HNL activity, and has a new catalytic center, as shown by protein crystallography. Biochemical properties of D. tyermannii HNLs open perspectives for the development of a complementary class of biocatalysts for the stereoselective synthesis of cyanohydrins. This work shows that systematic integration of -omics data facilitates discovery of enzymes with unpredictable sequences and helps to extend our knowledge about enzyme diversity.

Migrated Hopene Synthases from Colysis pothifolia and Identification of a Migration Switch Controlling the Number of 1,2-Hydride and Methyl Shifts.

Ferns are known to produce triterpenes, derived from squalene, that are synthesized by squalene cyclases (SCs). Among these, Colysis species produce onoceroids, the bis-cyclic skeleton of which can be cyclized from both termini of squalene. To gain insight into the molecular basis of triterpene structural diversity, cDNA cloning of SCs from C. elliptica and C. pothifolia was performed; this leads to the isolation of five SC cDNAs. Functional analysis of these clones revealed their enzymatic products to be hop-22(29)-ene, α-polypodatetraene, and hop-17(21)-ene. One of these clones (CPF) is a transcribed pseudogene with a 22-nucleotide deletion causing a nonsense mutation. To predict the inherent function of CPF, we constructed an insertion mutant of CPF that successfully converted inert CPF to the active SC, the product of which is fern-9(11)-ene. Subsequent mutations identified active-site residues that control the number of 1,2-hydride and methyl shifts.

Biochemical and growth performance of the aquatic macrophyte Azolla filiculoides to sub-chronic exposure to cylindrospermopsin.

Physiological and biochemical effects of cylindrospermopsin (CYN), a cyanobacterial toxin that inhibits protein synthesis and released during a harmful cyanobacterial bloom, has been overlooked in plants. Therefore, at the present research, the toxic effects (physiological and biochemical) of a crude extract containing CYN were assessed in the aquatic fern Azolla filiculoides exposed to three concentrations (0.05, 0.5 and 5 µg CYN mL(-1)). At 5 µg CYN mL(-1), fern growth rate has showed a drastic decrease (0.001 g g(-1) day(-1)) corresponding to a 99.8% inhibition, but at the concentrations of 0.05 and 0.5 µg CYN mL(-1) the growth rate was similar to the control plants. Growth rate also indicated a IC50 of 2.9 µg CYN mL(-1). Those data point to the presence of other compounds in the crude extract may stimulate the fern growth and/or the fern is tolerant to CYN. Chlorophyll (a and b), carotenoids and protein content as well as the activities of glutathione reductase (GR) and glutathione-S-transferase (GST) has increased at 5 µg CYN mL(-1) which may indicate that photosynthesis and protein synthesis are not affected by CYN and the probable activation of defense and detoxifying mechanisms to overcome the effects induced by the presence of CYN. Low uptake of cylindrospermopsin (1.314 µg CYN g(-1) FW) and low bioconcentration factor (0.401) point towards to a safe use of A. filiculoides as biofertilizer and as food source, but also indicate that the fern is not suitable for CYN phytoremediation.

Feasibility Study of Phragmites karka and Christella dentata Grown in West Bengal as Arsenic Accumulator.

A survey was undertaken, in arsenic (As) contaminated area of the Nadia district, West Bengal, India, to find native As accumulator plants. As was determined both in soil and plant parts. The results showed that the mean translocation factor of Pteris vittata L, Phragmites karka (Cav.) Trin. Ex. Steud and Christella dentata Forssk were higher than 1. It thus appeared that these plants can be efficient accumulators of As. Phytoremediation ability of C. dentata and P. karka was evaluated and compared with known As-hyperaccumulators -P. vittata and Adiantum capillus veneris L. Plants were grown in the As spiked soil (25, 50, 75 and 100 mg kg(-1)). As accumulation was found to be highest in P. vittata, 117.18 mg kg(-1) in leaf at 100 mg kg(-1) As treatment, followed by A. capillus veneris, P. karka and C. dentata being 74, 83.87 and 40.36 mg kg(-1), respectively. Lipid peroxidation increased after As exposure in all plants. However, the antioxidant enzyme activity and molecules concentration also increased which helped the plants to overcome As-induced oxidative stress. The study indicates that P. karka and C. dentata could be considered as As-accumulators and find application for As-phytoextraction in field conditions.

Structure and function of CrACA1, the major PM-type Ca2+-ATPase, expressed at the peak of the gravity-directed trans-cell calcium current in spores of the fern Ceratopteris richardii.

Spores of the fern Ceratopteris richardii have proven to be a valuable single-cell system for studying gravity responses. The earliest cellular change directed by gravity in these cells is a trans-cell calcium current, which peaks near 10 h after the spores are induced to germinate. This current is needed for gravity-directed axis alignment, and its peak is coincident with the time period when gravity polarises the direction of subsequent nuclear migration and rhizoid growth. Transcriptomic analysis of genes expressed at the 10-h time point revealed several that encode proteins likely to be key components that either drive the current or regulate it. Notable among these is a plasma membrane (PM)-type Ca(2+) ATPase, CrACA1, whose activity pumping Ca(2+) out of cells is regulated by gravity. This report provides an initial characterisation of the structure and expression of this protein, and demonstrates its heterologous function complementing the K616 mutant of yeast, which is deficient in PM-type Ca(2+) pump activity. Gravity-induced changes in the trans-cell Ca(2+) current occur within seconds, a result consistent with the hypothesis that the force of gravity can rapidly alter the post-translational state of the channels and pumps that drive this current across spore cells. This report identifies a transporter likely to be a key driver of the current, CrACA1, and characterises the role of this protein in early germination and gravity-driven polarity fixation through analysis of expression levels, functional complementation and pharmacological treatments. These data, along with newly available transcriptomic data obtained at the 10-h time point, indicate that CrACA1 is present, functional and likely a major contributing component of the trans-cell Ca(2+) efflux. CrACA1 is not necessary for polar axis alignment, but pharmacological perturbations of it disrupt rhizoid development. These data support and help refine the post-translational modification model for gravity responses.

Cyclization of all-E- and 2Z-geranylfarnesols by a bacterial triterpene synthase: insight into sesterterpene biosynthesis in Aleuritopteris ferns.

Aleuritopteris ferns produce triterpenes and sesterterpenes with tricyclic cheilanthane and tetracyclic 18-episcalarane skeletons. The structural and mechanistic similarities between both classes of fern terpene suggest that their biosynthetic enzymes may be closely related. We investigate here whether a triterpene synthase is capable of recognizing geranylfarnesols as a substrate, and is able to convert them to cyclic sesterterpenes. We found that a bacterial triterpene synthase converted all-E-geranylfarnesol (1b) into three scalarane sesterterpenes with 18αH stereochemistry (5, 7 and 8), as well as mono- and tricyclic sesterterpenes (6 and 9). In addition, 2Z-geranylfarnesol (4) was converted into an 18-episcalarane derivative (10), whose skeleton can be found in sesterterpenes isolated from Aleuritopteris ferns. These results provide insight into sesterterpene biosynthesis in Aleuritopteris ferns.

Gene flow among populations of two rare co-occurring fern species differing in ploidy level.

Differences in ploidy levels among different fern species have a vast influence on their mating system, their colonization ability and on the gene flow among populations. Differences in the colonization abilities of species with different ploidy levels are well known: tetraploids, in contrast to diploids, are able to undergo intra-gametophytic selfing. Because fertilization is a post-dispersal process in ferns, selfing results in better colonization abilities in tetraploids because of single spore colonization. Considerably less is known about the gene flow among populations of different ploidy levels. The present study examines two rare fern species that differ in ploidy. While it has already been confirmed that tetraploid species are better at colonizing, the present study focuses on the gene flow among existing populations. We analyzed the genetic structure of a set of populations in a 10×10 km study region using isoenzymes. Genetic variation in tetraploid species is distributed mainly among populations; the genetic distance between populations is correlated with the geographical distance, and larger populations host more genetic diversity than smaller populations. In the diploid species, most variability is partitioned within populations; the genetic distance is not related to geographic distance, and the genetic diversity of populations is not related to the population size. This suggests that in tetraploid species, which undergo selfing, gene flow is limited. In contrast, in the diploid species, which experience outcrossing, gene flow is extensive and the whole system behaves as one large population. Our results suggest that in ferns, the ability to colonize new habitats and the gene flow among existing populations are affected by the mating system.

Dynamics of polyploid formation and establishment in the allotetraploid rock fern Asplenium majoricum.

BACKGROUND AND AIMS: Successful establishment of newly formed polyploid species depends on several interlinked genetic and ecological factors. These include genetic diversity within and among individuals, chromosome behaviour and fertility, novel phenotypes resulting from novel genomic make-up and expression, intercytotypic and interspecific competition, and adaptation to distinct habitats. The allotetraploid rock fern Asplenium majoricum is known from one small population in Valencia, Spain, and several larger populations on the Balearic island of Majorca. In Valencia, it occurs sympatrically with its diploid parents, A. fontanum subsp. fontanum and A. petrarchae subsp. bivalens, and their diploid hybrid A. × protomajoricum. This highly unusual situation allowed the study of polyploid genetic diversity and its relationship to the formation and establishment of nascent polyploid lineages. METHODS: Genetic variation for isozyme and chloroplast DNA markers was determined for A. majoricum and A. × protomajoricum sampled thoroughly from known sites in Majorca and Valencia. Results were compared with variation determined previously for the diploid parent taxa. KEY RESULTS: A highly dynamic system with recurring diploid hybrid and allotetraploid formation was discovered. High diversity in the small Valencian A. majoricum population indicates multiple de novo origins from diverse parental genotypes, but most of these lineages become extinct without becoming established. The populations on Majorca most probably represent colonization(s) from Valencia rather than an in situ origin. Low genetic diversity suggests that this colonization may have occurred only once. CONCLUSIONS: There is a striking contrast in success of establishment of the Majorcan and Valencian populations of A. majoricum. Chance founding of populations in a habitat where neither A. fontanum subsp. fontanum nor A. petrarchae subsp. bivalens occurs appears to have been a key factor enabling the establishment of A. majoricum on Majorca. Successful establishment of this polyploid is probably dependent on geographic isolation from diploid progenitor competition.

Evolution of the rpoB-psbZ region in fern plastid genomes: notable structural rearrangements and highly variable intergenic spacers.

BACKGROUND: The rpoB-psbZ (BZ) region of some fern plastid genomes (plastomes) has been noted to go through considerable genomic changes. Unraveling its evolutionary dynamics across all fern lineages will lead to clarify the fundamental process shaping fern plastome structure and organization. RESULTS: A total of 24 fern BZ sequences were investigated with taxon sampling covering all the extant fern orders. We found that: (i) a tree fern Plagiogyria japonica contained a novel gene order that can be generated from either the ancestral Angiopteris type or the derived Adiantum type via a single inversion; (ii) the trnY-trnE intergenic spacer (IGS) of the filmy fern Vandenboschia radicans was expanded 3-fold due to the tandem 27-bp repeats which showed strong sequence similarity with the anticodon domain of trnY; (iii) the trnY-trnE IGSs of two horsetail ferns Equisetum ramosissimum and E. arvense underwent an unprecedented 5-kb long expansion, more than a quarter of which was consisted of a single type of direct repeats also relevant to the trnY anticodon domain; and (iv) ycf66 has independently lost at least four times in ferns. CONCLUSIONS: Our results provided fresh insights into the evolutionary process of fern BZ regions. The intermediate BZ gene order was not detected, supporting that the Adiantum type was generated by two inversions occurring in pairs. The occurrence of Vandenboschia 27-bp repeats represents the first evidence of partial tRNA gene duplication in fern plastomes. Repeats potentially forming a stem-loop structure play major roles in the expansion of the trnY-trnE IGS.