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

Transcriptomic analyses in the gametophytes of the apomictic fern Dryopteris affinis

MAIN CONCLUSION

A novel genomic map of the apogamous gametophyte of the fern Dryopteris affinis unlocks oldest hindrance with this complex plant group, to gain insight into evo-devo approaches. The gametophyte of the fern Dryopteris affinis ssp. affinis represents a good model to explore the molecular basis of vegetative and reproductive development, as well as stress responses. Specifically, this fern reproduces asexually by apogamy, a peculiar case of apomixis whereby a sporophyte forms directly from a gametophytic cell without fertilization. Using RNA-sequencing approach, we have previously annotated more than 6000 transcripts. Here, we selected 100 of the inferred proteins homolog to those of Arabidopsis thaliana, which were particularly interesting for a detailed study of their potential functions, protein-protein interactions, and distance trees. As expected, a plethora of proteins associated with gametogenesis and embryogenesis in angiosperms, such as FERONIA (FER) and CHROMATING REMODELING 11 (CHR11) were identified, and more than a dozen candidates potentially involved in apomixis, such as ARGONAUTE family (AGO4, AGO9, and AGO 10), BABY BOOM (BBM), FASCIATED STEM4 (FAS4), FERTILIZATION-INDEPENDENT ENDOSPERM (FIE), and MATERNAL EFFECT EMBRYO ARREST29 (MEE29). In addition, proteins involved in the response to biotic and abiotic stresses were widely represented, as shown by the enrichment of heat-shock proteins. Using the String platform, the interactome revealed that most of the protein-protein interactions were predicted based on experimental, database, and text mining datasets, with MULTICOPY SUPPRESSOR OF IRA4 (MSI4) showing the highest number of interactions: 16. Lastly, some proteins were studied through distance trees by comparing alignments with respect to more distantly or closely related plant groups. This analysis identified DCL4 as the most distant protein to the predicted common ancestor. New genomic information in relation to gametophyte development, including apomictic reproduction, could expand our current vision of evo-devo approaches.

Proteogenomics-based functional genome research: approaches, applications, and perspectives in plants

Proteogenomics (PG) integrates the proteome with the genome and transcriptome to refine gene models and annotation. Coupled with single-cell (SC) assays, PG effectively distinguishes heterogeneity among cell groups. Affiliating spatial information to PG reveals the high-resolution circuitry within SC atlases. Additionally, PG can investigate dynamic changes in protein-coding genes in plants across growth and development as well as stress and external stimulation, significantly contributing to the functional genome. Here we summarize existing PG research in plants and introduce the technical features of various methods. Combining PG with other omics, such as metabolomics and peptidomics, can offer even deeper insights into gene functions. We argue that the application of PG will represent an important font of foundational knowledge for plants.

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.

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.

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.

A drought-driven model for the evolution of obligate apomixis in ferns: evidence from pellaeids (Pteridaceae)

PREMISE

Xeric environments impose major constraints on the fern life cycle, yet many lineages overcome these limitations by evolving apomixis. Here, we synthesize studies of apomixis in ferns and present an evidence-based model for the evolution and establishment of this reproductive strategy, focusing on genetic and environmental factors associated with its two defining traits: the production of "unreduced" spores (n = 2n) and the initiation of sporophytes from gametophyte tissue (i.e., diplospory and apogamy, respectively).

METHODS

We evaluated existing literature in light of the hypothesis that abiotic characteristics of desert environments (e.g., extreme diurnal temperature fluctuations, high light intensity, and water limitation) drive the evolution of obligate apomixis. Pellaeid ferns (Cheilanthoideae: Pteridaceae) were examined in detail, as an illustrative example. We reconstructed a plastid (rbcL, trnG-trnR, atpA) phylogeny for the clade and mapped reproductive mode (sexual versus apomictic) and ploidy across the resulting tree.

RESULTS

Our six-stage model for the evolution of obligate apomixis in ferns emphasizes the role played by drought and associated abiotic conditions in the establishment of this reproductive approach. Furthermore, our updated phylogeny of pellaeid ferns reveals repeated origins of obligate apomixis and shows an increase in the frequency of apomixis, and rarity of sexual reproduction, among taxa inhabiting increasingly dry North American deserts.

CONCLUSIONS

Our findings reinforce aspects of other evolutionary, physiological, developmental, and omics-based studies, indicating a strong association between abiotic factors and the establishment of obligate apomixis in ferns. Water limitation, in particular, appears critical to establishment of this reproductive mode.

Differential Proteomic Analysis to Identify Proteins Associated with Apomeiosis in Boehmeria tricuspis (Hance) Makino Using an iTRAQ-Based Strategy

Numerous candidate genes related to apomixis have been identified through transcriptomics; however, the molecular mechanism underlying apomixis remains unclear. Elucidation of the underlying mechanisms is essential to expand its application in crop breeding. Therefore, here, we employed the isobaric tags for a relative and absolute quantification labeling technology to investigate the protein expression in Boehmeria tricuspis generated through different reproductive modes at the functional megaspore stage. We identified 40 differential abundance proteins associated with apomeiosis, most of which were involved in "response to stress". Functional analysis suggested that lower levels of reactive oxygen species (ROS) play a role in inducing the development of apomeiosis. Proteins related to ROS regulation, cell wall modifications, and stability under heat stress play a crucial role in the development of diplosporic apomeiosis. Our results give evidence to the insight that stress can induce a switch from apomixis to sexuality by ROS content, and an increased composition of stress tolerance as well as secondary metabolites can buffer ROS effects. Precise coordination of these proteins involved in inter-related regulatory control mechanisms may act together in the transition from the sexual to apomixis development.

Controlling Apomixis: Shared Features and Distinct Characteristics of Gene Regulation

In higher plants, sexual and asexual reproduction through seeds (apomixis) have evolved as alternative strategies. As apomixis leads to the formation of clonal offspring, its great potential for agricultural applications has long been recognized. However, the genetic basis and the molecular control underlying apomixis and its evolutionary origin are to date not fully understood. Both in sexual and apomictic plants, reproduction is tightly controlled by versatile mechanisms regulating gene expression, translation, and protein abundance and activity. Increasing evidence suggests that interrelated pathways including epigenetic regulation, cell-cycle control, hormonal pathways, and signal transduction processes are relevant for apomixis. Additional molecular mechanisms are being identified that involve the activity of DNA- and RNA-binding proteins, such as RNA helicases which are increasingly recognized as important regulators of reproduction. Together with other factors including non-coding RNAs, their association with ribosomes is likely to be relevant for the formation and specification of the apomictic reproductive lineage. Subsequent seed formation appears to involve an interplay of transcriptional activation and repression of developmental programs by epigenetic regulatory mechanisms. In this review, insights into the genetic basis and molecular control of apomixis are presented, also taking into account potential relations to environmental stress, and considering aspects of evolution.

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.

Proteome Analysis of the Gametophytes of a Western Himalayan Fern Diplazium maximum Reveals Their Adaptive Responses to Changes in Their Micro-Environment

Ferns have survived changing habitats and environmental extremes of different eras, wherein, the exploratory haploid gametophytes are believed to have played a major role. Therefore, the proteome of in vitro grown gametophytes of a temperate Himalayan fern, Diplazium maximum in response to 0 (G0), 1 (G1), and 3% (G3) sucrose was studied. A total of 110 differentially abundant protein spots (DAPs) were obtained. Among these, only 67 could be functionally categorized as unique proteins involved in various metabolic processes. Calcium dependent proteins, receptor like kinases, G proteins, proteins related to hormonal signaling and their interaction with other pathways, and regulatory proteins were recorded indicating the involvement of five different signaling pathways. DAPs involved in the activation of genes and transcription factors of signaling and transduction pathways, transport and ion channels, cell-wall and structural proteins, defense, chaperons, energy metabolism, protein synthesis, modification, and turnover were identified. The gametophytes responded to changes in their micro-environment. There was also significant increase in prothallus biomass and conversion of two-dimensional prothalli into three-dimensional prothallus clumps at 3% sucrose. The three-D clumps had higher photosynthetic surface area and also closer proximity for sexual reproduction and sporophyte formation. Highest accumulation of proline, enhanced scavenging of reactive oxygen species (ROS) and DAPs of mostly, abiotic stress tolerance, secondary metabolite synthesis, and detoxification at 3% sucrose indicated an adaptive response of gametophytes. Protein Protein Interaction network and Principal Component analyses, and qRT-PCR validation of genes encoding 12 proteins of various metabolic processes indicated differential adjustment of gametophytes to different levels of sucrose in the culture medium. Therefore, a hypothetical mechanism was proposed to show that even slight changes in the micro-environment of D. maximum gametophytes triggered multiple mechanisms of adaptation. Many DAPs identified in the study have potential use in crop improvement and metabolic engineering programs, phytoremediation and environmental protection.

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.

A Plant Biologist's Toolbox to Study Translation

Across a broad range of species and biological questions, more and more studies are incorporating translation data to better assess how gene regulation occurs at the level of protein synthesis. The inclusion of translation data improves upon, and has been shown to be more accurate than, transcriptional studies alone. However, there are many different techniques available to measure translation and it can be difficult, especially for young or aspiring scientists, to determine which methods are best applied in specific situations. We have assembled this review in order to enhance the understanding and promote the utilization of translational methods in plant biology. We cover a broad range of methods to measure changes in global translation (e.g., radiolabeling, polysome profiling, or puromycylation), translation of single genes (e.g., fluorescent reporter constructs, toeprinting, or ribosome density mapping), sequencing-based methods to uncover the entire translatome (e.g., Ribo-seq or translating ribosome affinity purification), and mass spectrometry-based methods to identify changes in the proteome (e.g., stable isotope labeling by amino acids in cell culture or bioorthogonal noncanonical amino acid tagging). The benefits and limitations of each method are discussed with a particular note of how applications from other model systems might be extended for use in plants. In order to make this burgeoning field more accessible to students and newer scientists, our review includes an extensive glossary to define key terms.

Proteomics in non-human primates: utilizing RNA-Seq data to improve protein identification by mass spectrometry in vervet monkeys

Background

Shotgun proteomics utilizes a database search strategy to compare detected mass spectra to a library of theoretical spectra derived from reference genome information. As such, the robustness of proteomics results is contingent upon the completeness and accuracy of the gene annotation in the reference genome. For animal models of disease where genomic annotation is incomplete, such as non-human primates, proteogenomic methods can improve the detection of proteins by incorporating transcriptional data from RNA-Seq to improve proteomics search databases used for peptide spectral matching. Customized search databases derived from RNA-Seq data are capable of identifying unannotated genetic and splice variants while simultaneously reducing the number of comparisons to only those transcripts actively expressed in the tissue.

Results

We collected RNA-Seq and proteomic data from 10 vervet monkey liver samples and used the RNA-Seq data to curate sample-specific search databases which were analyzed in the program Morpheus. We compared these results against those from a search database generated from the reference vervet genome. A total of 284 previously unannotated splice junctions were predicted by the RNA-Seq data, 92 of which were confirmed by peptide spectral matches. More than half (53/92) of these unannotated splice variants had orthologs in other non-human primates, suggesting that failure to match these peptides in the reference analyses likely arose from incomplete gene model information. The sample-specific databases also identified 101 unique peptides containing single amino acid substitutions which were missed by the reference database. Because the sample-specific searches were restricted to actively expressed transcripts, the search databases were smaller, more computationally efficient, and identified more peptides at the empirically derived 1 % false discovery rate.

Conclusion

Proteogenomic approaches are ideally suited to facilitate the discovery and annotation of proteins in less widely studies animal models such as non-human primates. We expect that these approaches will help to improve existing genome annotations of non-human primate species such as vervet.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4279-0) contains supplementary material, which is available to authorized users.