Dynamic evolution of MADS-box genes in extant ferns via large-scale phylogenomic analysis

Introduction

Several studies of MADS-box transcription factors in flowering plants have been conducted, and these studies have indicated that they have conserved functions in floral organ development; MIKC-type MADS-box genes has been proved to be expanded in ferns, however, few systematic studies of these transcription factors have been conducted in non-seed plants. Although ferns and seed plants are sister groups, they exhibit substantial morphological differences.

Methods

Here, we clarified the evolution of MADS-box genes across 71 extant fern species using available transcriptome, genome, and gene expression data.

Results

We obtained a total of 2,512 MADS-box sequences, ranging from 9 to 89 per species. The most recent common ancestor (MRCA) of ferns contained approximately three type I genes and at least 5-6 type II MADS-box genes. The domains, motifs, expression of type I and type II proteins, and the structure of the both type genes were conserved in ferns as to other land plants. Within type II genes, MIKC*-type proteins are involved in gametophyte development in ferns; MIKCC-type proteins have broader expression patterns in ferns than in seed plants, and these protein sequences are likely conserved in extant seed plants and ferns because of their diverse roles in diploid sporophyte development. More than 90% of MADS-box genes are type II genes, and MIKCC genes, especially CRM1 and CRM6-like genes, have undergone a large expansion in leptosporangiate ferns; the diverse expression patterns of these genes might be related to the fuctional diversification and increased complexity of the plant body plan. Tandem duplication of CRM1 and CRM6-like genes has contributed to the expansion of MIKCC genes.

Conclusion or Discussion

This study provides new insights into the diversity, evolution, and functions of MADS-box genes in extant ferns.

Proteome and Interactome Linked to Metabolism, Genetic Information Processing, and Abiotic Stress in Gametophytes of Two Woodferns

Ferns and lycophytes have received scant molecular attention in comparison to angiosperms. The advent of high-throughput technologies allowed an advance towards a greater knowledge of their elusive genomes. In this work, proteomic analyses of heart-shaped gametophytes of two ferns were performed: the apomictic Dryopteris affinis ssp. affinis and its sexual relative Dryopteris oreades. In total, a set of 218 proteins shared by these two gametophytes were analyzed using the STRING database, and their proteome associated with metabolism, genetic information processing, and responses to abiotic stress is discussed. Specifically, we report proteins involved in the metabolism of carbohydrates, lipids, and nucleotides, the biosynthesis of amino acids and secondary compounds, energy, oxide-reduction, transcription, translation, protein folding, sorting and degradation, and responses to abiotic stresses. The interactome of this set of proteins represents a total network composed of 218 nodes and 1792 interactions, obtained mostly from databases and text mining. The interactions among the identified proteins of the ferns D. affinis and D. oreades, together with the description of their biological functions, might contribute to a better understanding of the function and development of ferns as well as fill knowledge gaps in plant evolution.

Comparison between the Gametophyte and the Sporophyte Transcriptomes of the Endangered Fern Vandenboschia speciosa

Genomic resources are essential to understanding the evolution and functional biology of organisms. Nevertheless, generating genomic resources from endangered species may be challenging due to the scarcity of available specimens and sampling difficulties. In this study, we compare the transcriptomes of the sporophyte and the gametophyte of the endangered fern Vandenboschia speciosa. After Illumina sequencing and de novo transcriptome assembly of the gametophyte, annotation proved the existence of cross-species contamination in the gametophyte sample. Thus, we developed an in silico decontamination step for the gametophyte sequences. Once the quality check of the decontaminated reads passed, we produced a de novo assembly with the decontaminated gametophyte reads (with 43,139 contigs) and another combining the sporophyte and in silico decontaminated gametophyte reads (with 42,918 contigs). A comparison of the enriched GO terms from the top 1000 most expressed transcripts from both tissues showed that the gametophyte GO term set was enriched in sequences involved in development, response to stress, and plastid organization, while the sporophyte GO term set had a larger representation of more general metabolic functions. This study complements the available genomic resources on the life cycle of the endangered fern Vandenboschia speciosa.

The Shared Proteome of the Apomictic Fern Dryopteris affinis ssp. affinis and Its Sexual Relative Dryopteris oreades

Ferns are a diverse evolutionary lineage, sister to the seed plants, which is of great ecological importance and has a high biotechnological potential. Fern gametophytes represent one of the simplest autotrophic, multicellular plant forms and show several experimental advantages, including a simple and space-efficient in vitro culture system. However, the molecular basis of fern growth and development has hardly been studied. Here, we report on a proteomic study that identified 417 proteins shared by gametophytes of the apogamous fern Dryopteris affinis ssp. affinis and its sexual relative Dryopteris oreades. Most proteins are predicted to localize to the cytoplasm, the chloroplast, or the nucleus, and are linked to enzymatic, binding, and structural activities. A subset of 145 proteins are involved in growth, reproduction, phytohormone signaling and biosynthesis, and gene expression, including homologs of SHEPHERD (SHD), HEAT SHOCK PROTEIN 90-5 (CR88), TRP4, BOBBER 1 (BOB1), FLAVONE 3'-O-METHYLTRANSFERASE 1 (OMT1), ZEAXANTHIN EPOXIDASE (ABA1), GLUTAMATE DESCARBOXYLASE 1 (GAD), and dsRNA-BINDING DOMAIN-LIKE SUPERFAMILY PROTEIN (HLY1). Nearly 25% of the annotated proteins are associated with responses to biotic and abiotic stimuli. As for biotic stress, the proteins PROTEIN SGT1 HOMOLOG B (SGT1B), SUPPRESSOR OF SA INSENSITIVE2 (SSI2), PHOSPHOLIPASE D ALPHA 1 (PLDALPHA1), SERINE/THREONINE-PROTEIN KINASE SRK2E (OST1), ACYL CARRIER PROTEIN 4 (ACP4), and NONHOST RESISTANCE TO P. S. PHASEOLICOLA1 (GLPK) are worth mentioning. Regarding abiotic stimuli, we found proteins associated with oxidative stress: SUPEROXIDE DISMUTASE[CU-ZN] 1 (CSD1), and GLUTATHIONE S-TRANSFERASE U19 (GSTU19), light intensity SERINE HYDROXYMETHYLTRANSFERASE 1 (SHM1) and UBIQUITIN-CONJUGATING ENZYME E2 35 (UBC35), salt and heavy metal stress included MITOCHONDRIAL PHOSPHATE CARRIER PROTEIN 3 (PHT3;1), as well as drought and thermotolerance: LEA7, DEAD-BOX ATP-DEPENDENT RNA HELICASE 38 (LOS4), and abundant heat-shock proteins and other chaperones. In addition, we identified interactomes using the STRING platform, revealing protein-protein associations obtained from co-expression, co-occurrence, text mining, homology, databases, and experimental datasets. By focusing on ferns, this proteomic study increases our knowledge on plant development and evolution, and may inspire future applications in crop species.

Genome Assembly and Analysis of the Flavonoid and Phenylpropanoid Biosynthetic Pathways in Fingerroot Ginger (Boesenbergia rotunda)

Boesenbergia rotunda (Zingiberaceae), is a high-value culinary and ethno-medicinal plant of Southeast Asia. The rhizomes of this herb have a high flavanone and chalcone content. Here we report the genome analysis of B. rotunda together with a complete genome sequence as a hybrid assembly. B. rotunda has an estimated genome size of 2.4 Gb which is assembled as 27,491 contigs with an N50 size of 12.386 Mb. The highly heterozygous genome encodes 71,072 protein-coding genes and has a 72% repeat content, with class I TEs occupying ~67% of the assembled genome. Fluorescence in situ hybridization of the 18 chromosome pairs at the metaphase showed six sites of 45S rDNA and two sites of 5S rDNA. An SSR analysis identified 238,441 gSSRs and 4604 EST-SSRs with 49 SSR markers common among related species. Genome-wide methylation percentages ranged from 73% CpG, 36% CHG and 34% CHH in the leaf to 53% CpG, 18% CHG and 25% CHH in the embryogenic callus. Panduratin A biosynthetic unigenes were most highly expressed in the watery callus. B rotunda has a relatively large genome with a high heterozygosity and TE content. This assembly and data (PRJNA71294) comprise a source for further research on the functional genomics of B. rotunda, the evolution of the ginger plant family and the potential genetic selection or improvement of gingers.

Involvement of Auxin Biosynthesis and Transport in the Antheridium and Prothalli Formation in Lygodium japonicum

The spores of Lygodium japonicum, cultured in the dark, form a filamentous structure called protonema. Earlier studies have shown that gibberellin (GA) induces protonema elongation, along with antheridium formation, on the protonema. In this study, we have performed detailed morphological analyses to investigate the roles of multiple phytohormones in antheridium formation, protonema elongation, and prothallus formation in L. japonicum. GA4 methyl ester is a potent GA that stimulates both protonema elongation and antheridium formation. We found that these effects were inhibited by simultaneous application of abscisic acid (ABA). On the other hand, IAA (indole-3-acetic acid) promoted protonema elongation but reduced antheridium formation, while these effects were partially recovered by transferring to an IAA-free medium. An auxin biosynthesis inhibitor, PPBo (4-phenoxyphenylboronic acid), and a transport inhibitor, TIBA (2,3,5-triiodobenzoic acid), both inhibited protonema elongation and antheridium formation. L. japonicum prothalli are induced from germinating spores under continuous white light. Such development was negatively affected by PPBo, which induced smaller-sized prothalli, and TIBA, which induced aberrantly shaped prothalli. The evidence suggests that the crosstalk between these plant hormones might regulate protonema elongation and antheridium formation in L. japonicum. Furthermore, the possible involvement of auxin in the prothalli development of L. japonicum is suggested.

Sexual and Apogamous Species of Woodferns Show Different Protein and Phytohormone Profiles

The gametophyte of ferns reproduces either by sexual or asexual means. In the latter, apogamy represents a peculiar case of apomixis, in which an embryo is formed from somatic cells. A proteomic and physiological approach was applied to the apogamous fern Dryopteris affinis ssp. affinis and its sexual relative D. oreades. The proteomic analysis compared apogamous vs. female gametophytes, whereas the phytohormone study included, in addition to females, three apogamous stages (filamentous, spatulate, and cordate). The proteomic profiles revealed a total of 879 proteins and, after annotation, different regulation was found in 206 proteins of D. affinis and 166 of its sexual counterpart. The proteins upregulated in D. affinis are mostly associated to protein metabolism (including folding, transport, and proteolysis), ribosome biogenesis, gene expression and translation, while in the sexual counterpart, they account largely for starch and sucrose metabolism, generation of energy and photosynthesis. Likewise, ultra-performance liquid chromatography-tandem spectrometry (UHPLC-MS/MS) was used to assess the levels of indol-3-acetic acid (IAA); the cytokinins: 6-benzylaminopurine (BA), trans-Zeatine (Z), trans-Zeatin riboside (ZR), dyhidrozeatine (DHZ), dyhidrozeatin riboside (DHZR), isopentenyl adenine (iP), isopentenyl adenosine (iPR), abscisic acid (ABA), the gibberellins GA₃ and GA₄, salicylic acid (SA), and the brassinosteroids: brassinolide (BL) and castasterone (CS). IAA, the cytokinins Z, ZR, iPR, the gibberellin GA₄, the brassinosteoids castasterone, and ABA accumulated more in the sexual gametophyte than in the apogamous one. When comparing the three apogamous stages, BA and SA peaked in filamentous, GA₃ and BL in spatulate and DHRZ in cordate gametophytes. The results point to the existence of large metabolic differences between apogamous and sexual gametophytes, and invite to consider the fern gametophyte as a good experimental system to deepen our understanding of plant reproduction.

De Novo Sporophyte Transcriptome Assembly and Functional Annotation in the Endangered Fern Species Vandenboschia speciosa (Willd.) G. Kunkel

We sequenced the sporophyte transcriptome of Killarney fern (Vandenboschia speciosa (Willd.) G. Kunkel). In addition to being a rare endangered Macaronesian-European endemism, this species has a huge genome (10.52 Gb) as well as particular biological features and extreme ecological requirements. These characteristics, together with the systematic position of ferns among vascular plants, make it of high interest for evolutionary, conservation and functional genomics studies. The transcriptome was constructed de novo and contained 36,430 transcripts, of which 17,706 had valid BLAST hits. A total of 19,539 transcripts showed at least one of the 7362 GO terms assigned to the transcriptome, whereas 6547 transcripts showed at least one of the 1359 KEGG assigned terms. A prospective analysis of functional annotation results provided relevant insights on genes involved in important functions such as growth and development as well as physiological adaptations. In this context, a catalogue of genes involved in the genetic control of plant development, during the vegetative to reproductive transition, in stress response as well as genes coding for transcription factors is given. Altogether, this study provides a first step towards understanding the gene expression of a significant fern species and the in silico functional and comparative analyses reported here provide important data and insights for further comparative evolutionary studies in ferns and land plants in general.

Conservation and diversification of HAIRY MERISTEM gene family in land plants

The shoot apical meristems (SAMs) of land plants are crucial for plant growth and organ formation. In several angiosperms, the HAIRY MERISTEM (HAM) genes function as key regulators that control meristem development and stem cell homeostasis. To date, the origin and evolutionary history of the HAM family in land plants remains unclear. Potentially shared and divergent functions of HAM family members from angiosperms and non-angiosperms are also not known. In constructing a comprehensive phylogeny of the HAM family, we show that HAM proteins are widely present in land plants and that HAM proteins originated prior to the divergence of bryophytes. The HAM family was duplicated in a common ancestor of angiosperms, leading to two distinct groups: type I and type II. Type-II HAM members are widely present in angiosperms, whereas type-I HAM members were independently lost in different orders of monocots. Furthermore, HAM members from angiosperms and non-angiosperms (including bryophytes, lycophytes, ferns and gymnosperms) are able to replace the role of the type-II HAM genes in Arabidopsis, maintaining established SAMs and promoting the initiation of new stem cell niches. Our results uncover the conserved functions of HAM family members and reveal the conserved regulatory mechanisms underlying HAM expression patterning in meristems, providing insight into the evolution of key stem cell regulators in land plants.

Origin and diversification of ECERIFERUM1 (CER1) and ECERIFERUM3 (CER3) genes in land plants and phylogenetic evidence that the ancestral CER1/3 gene resulted from the fusion of pre-existing domains

ECERIFERUM1 (CER1) and ECERIFERUM3 (CER3) are key genes in synthesis of alkanes, a major component of cuticular waxes in land plants. The genes share extensive similarity, including the N-terminal (ERG3/FAH) and C-terminal (WAX2) domains. This study, traces the origin, evolutionary history, phylogenetic relationships and variation in copy number of the two genes within and beyond the Viridiplantae (green plants). Protein homologs of both CER1 and CER3 were identified across most Embryophyta (land plants), a single homolog (CER1/3) in charophytes and prasinophytes, and none in the other green, red or brown algae. Ancestral state reconstructions in 100 sequenced Archaeplastida using presence/absence of CER1/3 family genes revealed that the CER1/3 gene probably originated in the common ancestor of Viridiplantae. Phylogenetic analysis of CER1 and CER3 protein sequences from 146 plant species strongly suggests that the two genes originated by duplication of CER1/3 in the ancestral embryophyte. The evolution of CER1 and CER3 genes involved differential divergence of the two domains. Outside Embryophyta, CER1/3 similar sequences identified in diatoms and a cryptophyte, were the closest relatives of the CER1/3 family proteins. Proteins harbouring WAX2-wxAR (WAX2 associated region) similar regions were identified in proteins of bacteria, Archaea, cryptophytes, dinoflagellates and Stramenopiles. The independent existence of both ERG3/FAH and WAX2-wxAR domains in diverse lineages strongly points to the origin of CER1/3 gene in green plants by the fusion of pre-existing domains.

Leaf dorsoventrality candidate gene CpARF4 has conserved expression pattern but divergent tasiR-ARF regulation in the water fern Ceratopteris pteridoides

PREMISE

Leaves are traditionally classified into microphylls and megaphylls, and recently have been regarded as independently originating in lycophytes, ferns, and seed plants. The developmental genetics of leaf dorsoventrality, a synapomorphy in vascular plants, has been extensively studied in flowering plants. AUXIN RESPONSE FACTOR4 (ARF4) genes are key to leaf abaxial identity in flowering plants, but whether they exist in ferns is still an open question.

METHODS

ARF4 genes from Ceratopteris pteridoides, Cyrtomium guizhouense, and Parathelypteris nipponica were mined from transcriptomes and investigated in terms of evolutionary phylogeny and sequence motifs, with a focus on the tasiR-ARF binding site. In situ hybridization was used to localize expression of CpARF4 in Ceratopteris pteridoides. 5'RNA ligase-mediated-RACE was employed to verify whether CpARF4 transcripts were sliced by tasiR-ARF.

RESULTS

ARF4 genes exist in ferns, and this lineage originates from a gene duplication in the common ancestor of ferns and seed plants. ARF4 genes are of a single copy in the ferns studied here, and they contain divergent and, at most, one tasiR-ARF binding site. CpARF4 is expressed in the abaxial but not the adaxial domain of leaf primordia at various developmental stages. Transcript slicing guided by tasiR-ARF is active in C. pteridoides, but CpARF4 probably has not been affected by it.

CONCLUSIONS

Fern ARF4 genes differ in copy number and tasiR-ARF regulation relative to flowering plants, though they can be similarly expressed in the abaxial domain of leaves, revealing a key role for ARF4 genes in the evolution of leaf dorsoventrality of vascular plants.

Phylogenomics indicates the "living fossil" Isoetes diversified in the Cenozoic

The fossil record provides an invaluable insight into the temporal origins of extant lineages of organisms. However, establishing the relationships between fossils and extant lineages can be difficult in groups with low rates of morphological change over time. Molecular dating can potentially circumvent this issue by allowing distant fossils to act as calibration points, but rate variation across large evolutionary scales can bias such analyses. In this study, we apply multiple dating methods to genome-wide datasets to infer the origin of extant species of Isoetes, a group of mostly aquatic and semi-aquatic isoetalean lycopsids, which closely resemble fossil forms dating back to the Triassic. Rate variation observed in chloroplast genomes hampers accurate dating, but genome-wide nuclear markers place the origin of extant diversity within this group in the mid-Paleogene, 45-60 million years ago. Our genomic analyses coupled with a careful evaluation of the fossil record indicate that despite resembling forms from the Triassic, extant Isoetes species do not represent the remnants of an ancient and widespread group, but instead have spread around the globe in the relatively recent past.

Differential gene expression profiling of one- and two-dimensional apogamous gametophytes of the fern Dryopteris affinis ssp. affinis

Apomixis was originally defined as the replacement of sexual reproduction by an asexual process that does not involve fertilization but, in angiosperms, it is often used in the more restricted sense of asexual reproduction through seeds. In ferns, apomixis combines the production of unreduced spores (diplospory) and the formation of sporophytes from somatic cells of the prothallium (apogamy). The genes that control the onset of apogamy in ferns are largely unknown. In this study, we describe the gametophyte transcriptome of the apogamous fern Dryopteris affinis ssp. affinis using an RNA-Seq approach to compare the gene expression profiles of one- and two-dimensional gametophytes, the latter containing apogamic centers. After collapsing highly similar de novo transcripts, we obtained 166,191 unigenes, of which 30% could be annotated using public databases. Multiple quality metrics indicate a good quality of the de novo transcriptome with a low level of fragmentation. Our data show a total of 10,679 genes (6% of all genes) to be differentially expressed between gametophytes of filamentous (one-dimensional) and prothallial (two-dimensional) architecture. 6,110 genes were up-regulated in two-dimensional relative to one-dimensional gametophytes, some of which are implicated in the regulation of meristem growth, auxin signaling, reproduction, and sucrose metabolism. 4,570 genes were down-regulated in two-dimensional versus one-dimensional gametophytes, which are enriched in stimulus and defense genes, as well as genes involved in epigenetic gene regulation and ubiquitin degradation. Our results provide insights into free-living gametophyte development, focusing on the filamentous-to-prothallus growth transition, and provide a useful resource for further investigations of asexual reproduction.

The evolutionary trajectory of root stem cells

Root stem cells are crucial for the establishment of roots and are therefore a major evolutionary innovation that enabled land plants to spread on land. Despite their importance, not too much is known about the origin and the molecular players installing and maintaining them. Although still fragmentary, the recent availability of new data for early land plants can be used to identify and analyze the conservation of key regulators of root meristems. In this review, we evaluate the possible conservation of important root stem cell regulators to suggest pathways that might have been important at the origin of roots.

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.

Comparative transcriptome analysis of two reproductive modes in Adiantum reniforme var. sinense targeted to explore possible mechanism of apogamy

Background

Apogamy is a unique asexual reproduction in the ferns, in which somatic cells of gametophytes go through dedifferentiation and then differentiate into haploid sporophytes bypassing fertilization. Restricted to the lack of genomic information, molecular mechanisms of apogamy have remained unclear. Comparative transcriptome analysis was conducted at six stages between sexual reproduction and apogamy in the fern Adiantum reniforme var. sinense, in an effort to identify genes and pathways that might initiate the asexual reproduction.

Results

Approximately 928 million high-quality clean reads were assembled into 264,791 unigenes with an average length of 615 bp. A total of 147,865 (55.84%) unigenes were successfully annotated. Differential genes expression analysis indicated that transcriptional regulation was more active in the early stage of apogamy compared to sexual reproduction. Further comparative analysis of the enriched pathways between the early stages of the two reproductive modes demonstrated that starch and sucrose metabolism pathway responsible for cell wall was only significantly enriched in asexual embryonic cell initiation. Furthermore, regulation of plant hormone related genes was more vigorous in apogamy initiation.

Conclusion

These findings would be useful for revealing the initiation of apogamy and further understanding of the mechanisms related to asexual reproduction.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0762-8) contains supplementary material, which is available to authorized users.

Overlapping Patterns of Gene Expression Between Gametophyte and Sporophyte Phases in the Fern Polypodium amorphum (Polypodiales)

Ferns are unique among land plants in having sporophyte and gametophyte phases that are both free living and fully independent. Here, we examine patterns of sporophytic and gametophytic gene expression in the fern Polypodium amorphum, a member of the homosporous polypod lineage that comprises 80% of extant fern diversity, to assess how expression of a common genome is partitioned between two morphologically, ecologically, and nutritionally independent phases. Using RNA-sequencing, we generated transcriptome profiles for three replicates of paired samples of sporophyte leaf tissue and whole gametophytes to identify genes with significant differences in expression between the two phases. We found a nearly 90% overlap in the identity and expression levels of the genes expressed in both sporophytes and gametophytes, with less than 3% of genes uniquely expressed in either phase. We compare our results to those from similar studies to establish how phase-specific gene expression varies among major land plant lineages. Notably, despite having greater similarity in the identity of gene families shared between P. amorphum and angiosperms, P. amorphum has phase-specific gene expression profiles that are more like bryophytes and lycophytes than seed plants. Our findings suggest that shared patterns of phase-specific gene expression among seed-free plants likely reflect having relatively large, photosynthetic gametophytes (compared to the gametophytes of seed plants that are highly reduced). Phylogenetic analyses were used to further investigate the evolution of phase-specific expression for the phototropin, terpene synthase, and MADS-box gene families.

The Physcomitrella patens gene atlas project: large-scale RNA-seq based expression data

High-throughput RNA sequencing (RNA-seq) has recently become the method of choice to define and analyze transcriptomes. For the model moss Physcomitrella patens, although this method has been used to help analyze specific perturbations, no overall reference dataset has yet been established. In the framework of the Gene Atlas project, the Joint Genome Institute selected P. patens as a flagship genome, opening the way to generate the first comprehensive transcriptome dataset for this moss. The first round of sequencing described here is composed of 99 independent libraries spanning 34 different developmental stages and conditions. Upon dataset quality control and processing through read mapping, 28 509 of the 34 361 v3.3 gene models (83%) were detected to be expressed across the samples. Differentially expressed genes (DEGs) were calculated across the dataset to permit perturbation comparisons between conditions. The analysis of the three most distinct and abundant P. patens growth stages - protonema, gametophore and sporophyte - allowed us to define both general transcriptional patterns and stage-specific transcripts. As an example of variation of physico-chemical growth conditions, we detail here the impact of ammonium supplementation under standard growth conditions on the protonemal transcriptome. Finally, the cooperative nature of this project allowed us to analyze inter-laboratory variation, as 13 different laboratories around the world provided samples. We compare differences in the replication of experiments in a single laboratory and between different laboratories.

Sequencing of small RNAs of the fern Pleopeltis minima (Polypodiaceae) offers insight into the evolution of the microrna repertoire in land plants

MicroRNAs (miRNAs) are short, single stranded RNA molecules that regulate the stability and translation of messenger RNAs in diverse eukaryotic groups. Several miRNA genes are of ancient origin and have been maintained in the genomes of animal and plant taxa for hundreds of millions of years, playing key roles in development and physiology. In the last decade, genome and small RNA (sRNA) sequencing of several plant species have helped unveil the evolutionary history of land plants. Among these, the fern group (monilophytes) occupies a key phylogenetic position, as it represents the closest extant cousin taxon of seed plants, i.e. gymno- and angiosperms. However, in spite of their evolutionary, economic and ecological importance, no fern genome has been sequenced yet and few genomic resources are available for this group. Here, we sequenced the small RNA fraction of an epiphytic South American fern, Pleopeltis minima (Polypodiaceae), and compared it to plant miRNA databases, allowing for the identification of miRNA families that are shared by all land plants, shared by all vascular plants (tracheophytes) or shared by euphyllophytes (ferns and seed plants) only. Using the recently described transcriptome of another fern, Lygodium japonicum, we also estimated the degree of conservation of fern miRNA targets in relation to other plant groups. Our results pinpoint the origin of several miRNA families in the land plant evolutionary tree with more precision and are a resource for future genomic and functional studies of fern miRNAs.

What are the evolutionary origins of stomatal responses to abscisic acid in land plants?

The evolution of active stomatal closure in response to leaf water deficit, mediated by the hormone abscisic acid (ABA), has been the subject of recent debate. Two different models for the timing of the evolution of this response recur in the literature. A single-step model for stomatal control suggests that stomata evolved active, ABA-mediated control of stomatal aperture, when these structures first appeared, prior to the divergence of bryophyte and vascular plant lineages. In contrast, a gradualistic model for stomatal control proposes that the most basal vascular plant stomata responded passively to changes in leaf water status. This model suggests that active ABA-driven mechanisms for stomatal responses to water status instead evolved after the divergence of seed plants, culminating in the complex, ABA-mediated responses observed in modern angiosperms. Here we review the findings that form the basis for these two models, including recent work that provides critical molecular insights into resolving this intriguing debate, and find strong evidence to support a gradualistic model for stomatal evolution.

Proteogenomic Analysis Greatly Expands the Identification of Proteins Related to Reproduction in the Apogamous Fern Dryopteris affinis ssp. affinis

Performing proteomic studies on non-model organisms with little or no genomic information is still difficult. However, many specific processes and biochemical pathways occur only in species that are poorly characterized at the genomic level. For example, many plants can reproduce both sexually and asexually, the first one allowing the generation of new genotypes and the latter their fixation. Thus, both modes of reproduction are of great agronomic value. However, the molecular basis of asexual reproduction is not well understood in any plant. In ferns, it combines the production of unreduced spores (diplospory) and the formation of sporophytes from somatic cells (apogamy). To set the basis to study these processes, we performed transcriptomics by next-generation sequencing (NGS) and shotgun proteomics by tandem mass spectrometry in the apogamous fern D. affinis ssp. affinis. For protein identification we used the public viridiplantae database (VPDB) to identify orthologous proteins from other plant species and new transcriptomics data to generate a "species-specific transcriptome database" (SSTDB). In total 1,397 protein clusters with 5,865 unique peptide sequences were identified (13 decoy proteins out of 1,410, protFDR 0.93% on protein cluster level). We show that using the SSTDB for protein identification increases the number of identified peptides almost four times compared to using only the publically available VPDB. We identified homologs of proteins involved in reproduction of higher plants, including proteins with a potential role in apogamy. With the increasing availability of genomic data from non-model species, similar proteogenomics approaches will improve the sensitivity in protein identification for species only distantly related to models.

Orchidstra 2.0-A Transcriptomics Resource for the Orchid Family

Orchidaceae, the orchid family, encompasses more than 25,000 species and five subfamilies. Due to their beautiful and exotic flowers, distinct biological and ecological features, orchids have aroused wide interest among both researchers and the general public. We constructed the Orchidstra database, a resource for orchid transcriptome assembly and gene annotations. The Orchistra database has been under active development since 2013. To accommodate the increasing amount of orchid transcriptome data and house more comprehensive information, Orchidstra 2.0 has been built with a new database system to store the annotations of 510,947 protein-coding genes and 161,826 noncoding transcripts, covering 18 orchid species belonging to 12 genera in five subfamilies of Orchidaceae. We have improved the N50 size of protein-coding genes, provided new functional annotations (including protein-coding gene annotations, protein domain/family information, pathways analysis, Gene Ontology term assignments, orthologous genes across orchid species, cross-links to the database of model species, and miRNA information), and improved the user interface with better website performance. We also provide new database functionalities for database searching and sequence retrieval. Moreover, the Orchidstra 2.0 database incorporates detailed RNA-Seq gene expression data from various tissues and developmental stages in different orchid species. The database will be useful for gene prediction and gene family studies, and for exploring gene expression in orchid species. The Orchidstra 2.0 database is freely accessible at http://orchidstra2.abrc.sinica.edu.tw.

Transcriptome analyses provide insights into the phylogeny and adaptive evolution of the mangrove fern genus Acrostichum

The mangrove fern genus Acrostichum grows in the extremely unstable marine intertidal zone under harsh conditions, such as high salt concentrations, tidal rhythms and long-term climate changes. To explore the phylogenetic relationships and molecular mechanisms underlying adaptations in this genus, we sequenced the transcriptomes of two species of Acrostichum, A. aureum and A. speciosum, as well as a species in the sister genus, Ceratopteris thalictroides. We obtained 47,517, 36,420 and 60,823 unigenes for the three ferns, of which 24.39-45.63% were annotated using public databases. The estimated divergence time revealed that Acrostichum adapted to the coastal region during the late Cretaceous, whereas the two mangrove ferns from the Indo West-Pacific (IWP) area diverged more recently. Two methods (the modified branch-site model and the K_h method) were used to identify several positively selected genes, which may contribute to differential adaptation of the two Acrostichum species to different light and salt conditions. Our study provides abundant transcriptome data and new insights into the evolution and adaptations of mangrove ferns in the inhospitable intertidal zone.

Transcriptome sequencing reveals genome-wide variation in molecular evolutionary rate among ferns

BACKGROUND

Transcriptomics in non-model plant systems has recently reached a point where the examination of nuclear genome-wide patterns in understudied groups is an achievable reality. This progress is especially notable in evolutionary studies of ferns, for which molecular resources to date have been derived primarily from the plastid genome. Here, we utilize transcriptome data in the first genome-wide comparative study of molecular evolutionary rate in ferns. We focus on the ecologically diverse family Pteridaceae, which comprises about 10 % of fern diversity and includes the enigmatic vittarioid ferns-an epiphytic, tropical lineage known for dramatically reduced morphologies and radically elongated phylogenetic branch lengths. Using expressed sequence data for 2091 loci, we perform pairwise comparisons of molecular evolutionary rate among 12 species spanning the three largest clades in the family and ask whether previously documented heterogeneity in plastid substitution rates is reflected in their nuclear genomes. We then inquire whether variation in evolutionary rate is being shaped by genes belonging to specific functional categories and test for differential patterns of selection.

RESULTS

We find significant, genome-wide differences in evolutionary rate for vittarioid ferns relative to all other lineages within the Pteridaceae, but we recover few significant correlations between faster/slower vittarioid loci and known functional gene categories. We demonstrate that the faster rates characteristic of the vittarioid ferns are likely not driven by positive selection, nor are they unique to any particular type of nucleotide substitution.

CONCLUSIONS

Our results reinforce recently reviewed mechanisms hypothesized to shape molecular evolutionary rates in vittarioid ferns and provide novel insight into substitution rate variation both within and among fern nuclear genomes.

Assembly, Assessment, and Availability of De novo Generated Eukaryotic Transcriptomes

De novo assembly of a complete transcriptome without the need for a guiding reference genome is attractive, particularly where the cost and complexity of generating a eukaryote genome is prohibitive. The transcriptome should not however be seen as just a quick and cheap alternative to building a complete genome. Transcriptomics allows the understanding and comparison of spatial and temporal samples within an organism, and allows surveying of multiple individuals or closely related species. De novo assembly in theory allows the building of a complete transcriptome without any prior knowledge of the genome. It also allows the discovery of alternate splice forms of coding RNAs and also non-coding RNAs, which are often missed by proteomic approaches, or are incompletely annotated in genome studies. The limitations of the method are that the generation of a truly complete assembly is unlikely, and so we require some methods for the assessment of the quality and appropriateness of a generated transcriptome. Whilst no single consensus pipeline or tool is agreed as optimal, various algorithms, and easy to use software do exist making transcriptome generation a more common approach. With this expansion of data, questions still exist relating to how do we make these datasets fully discoverable, comparable and most useful to understand complex biological systems?

Pleurochrysome: A Web Database of Pleurochrysis Transcripts and Orthologs Among Heterogeneous Algae

Pleurochrysis is a coccolithophorid genus, which belongs to the Coccolithales in the Haptophyta. The genus has been used extensively for biological research, together with Emiliania in the Isochrysidales, to understand distinctive features between the two coccolithophorid-including orders. However, molecular biological research on Pleurochrysis such as elucidation of the molecular mechanism behind coccolith formation has not made great progress at least in part because of lack of comprehensive gene information. To provide such information to the research community, we built an open web database, the Pleurochrysome (http://bioinf.mind.meiji.ac.jp/phapt/), which currently stores 9,023 unique gene sequences (designated as UNIGENEs) assembled from expressed sequence tag sequences of P. haptonemofera as core information. The UNIGENEs were annotated with gene sequences sharing significant homology, conserved domains, Gene Ontology, KEGG Orthology, predicted subcellular localization, open reading frames and orthologous relationship with genes of 10 other algal species, a cyanobacterium and the yeast Saccharomyces cerevisiae. This sequence and annotation information can be easily accessed via several search functions. Besides fundamental functions such as BLAST and keyword searches, this database also offers search functions to explore orthologous genes in the 12 organisms and to seek novel genes. The Pleurochrysome will promote molecular biological and phylogenetic research on coccolithophorids and other haptophytes by helping scientists mine data from the primary transcriptome of P. haptonemofera.

Evolution of the Cp-Actin-based Motility System of Chloroplasts in Green Plants

During the course of green plant evolution, numerous light responses have arisen that optimize their growth under fluctuating light conditions. The blue light receptor phototropin mediates several photomovement responses at the tissue, cellular and organelle levels. Chloroplast photorelocation movement is one such photomovement response, and is found not only in most green plants, but also in some red algae and photosynthetic stramenopiles. In general, chloroplasts move toward weak light to maximally capture photosynthetically active radiation (the chloroplast accumulation response), and they move away from strong light to avoid photodamage (the avoidance response). In land plants, chloroplast movement is dependent on specialized actin filaments, chloroplast-actin filaments (cp-actin filaments). Through molecular genetic analysis using Arabidopsis thaliana, many molecular factors that regulate chloroplast photorelocation were identified. In this Perspective, we discuss the evolutionary history of the molecular mechanism for chloroplast photorelocation movement in green plants in view of cp-actin filaments.

Telomere dynamics in the lower plant Physcomitrella patens

A comparative approach in biology is needed to assess the universality of rules governing this discipline. In plant telomere research, most of the key principles were established based on studies in only single model plant, Arabidopsis thaliana. These principles include the absence of telomere shortening during plant development and the corresponding activity of telomerase in dividing (meristem) plant cells. Here we examine these principles in Physcomitrella patens as a representative of lower plants. To follow telomerase expression, we first characterize the gene coding for the telomerase reverse transcriptase subunit PpTERT in P. patens, for which only incomplete prediction has been available so far. In protonema cultures of P. patens, growing by filament apical cell division, the proportion of apical (dividing) cells was quantified and telomere length, telomerase expression and activity were determined. Our results show telomere stability and demonstrate proportionality of telomerase activity and expression with the number of apical cells. In addition, we analyze telomere maintenance in mre11, rad50, nbs1, ku70 and lig4 mutants of P. patens and compare the impact of these mutations in double-strand-break (DSB) repair pathways with earlier observations in corresponding A. thaliana mutants. Telomere phenotypes are absent and DSB repair kinetics is not affected in P. patens mutants for DSB factors involved in non-homologous end joining (NHEJ). This is compliant with the overall dominance of homologous recombination over NHEJ pathways in the moss, contrary to the inverse situation in flowering plants.