Genome-wide discovery and validation of Eucalyptus small RNAs reveals variable patterns of conservation and diversity across species of Myrtaceae

Genome-wide characterization and analysis of Golden 2-Like transcription factors related to leaf chlorophyll synthesis in diploid and triploid Eucalyptus urophylla

Golden 2-Like (GLK) transcription factors play a crucial role in chloroplast development and chlorophyll synthesis in many plant taxa. To date, no systematic analysis of GLK transcription factors in tree species has been conducted. In this study, 40 EgrGLK genes in the Eucalyptus grandis genome were identified and divided into seven groups based on the gene structure and motif composition. The EgrGLK genes were mapped to 11 chromosomes and the distribution of genes on chromosome was uneven. Phylogenetic analysis of GLK proteins between E. grandis and other species provided information for the high evolutionary conservation of GLK genes among different species. Prediction of cis-regulatory elements indicated that the EgrGLK genes were involved in development, light response, and hormone response. Based on the finding that the content of chlorophyll in mature leaves was the highest, and leaf chlorophyll content of triploid Eucalyptus urophylla was higher than that of the diploid control, EgrGLK expression pattern in leaves of triploid and diploid E. urophylla was examined by means of transcriptome analysis. Differential expression of EgrGLK genes in leaves of E. urophylla of different ploidies was consistent with the trend in chlorophyll content. To further explore the relationship between EgrGLK expression and chlorophyll synthesis, co-expression networks were generated, which indicated that EgrGLK genes may have a positive regulatory relationship with chlorophyll synthesis. In addition, three EgrGLK genes that may play an important role in chlorophyll synthesis were identified in the co-expression networks. And the prediction of miRNAs targeting EgrGLK genes showed that miRNAs might play an important role in the regulation of EgrGLK gene expression. This research provides valuable information for further functional characterization of GLK genes in Eucalyptus.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Genome-wide identification of microRNAs involved in the somatic embryogenesis of Eucalyptus

MicroRNAs (miRNAs) are small noncoding RNAs (18-24 nt) and function in many biological processes in plants. Although Eucalyptus trees are widely planted across the world, our understanding of the miRNA regulation in the somatic embryogenesis (SE) of Eucalyptus is still poor. Here we reported, for the first time, the miRNA profiles of differentiated and dedifferentiated tissues of two Eucalyptus species and identified miRNAs involved in SE of Eucalyptus. Stem and tissue culture-induced callus were obtained from the subculture seedlings of E. camaldulensis and E. grandis x urophylla and were used as differentiated and dedifferentiated samples, respectively. Small RNA sequencing generated 304.2 million clean reads for the Eucalyptus samples (n = 3) and identified 888 miRNA precursors (197 known and 691 novel) for Eucalyptus. These miRNAs were mainly distributed in chromosomes Chr03, Chr05, and Chr08 and can produce 46 miRNA clusters. Then, we identified 327 and 343 differentially expressed miRNAs (DEmiRs) in the dedifferentiation process of E. camaldulensis and E. grandis x urophylla, respectively. DEmiRs shared by the two Eucalyptus species might be involved in the development of embryonic callus, such as MIR156, MIR159, MIR160, MIR164, MIR166, MIR169, MIR171, MIR399, and MIR482. Notably, we identified 81 upregulated and 67 downregulated miRNAs specific to E. camaldulensis, which might be associated with the high embryogenic potential. Target prediction and functional analysis showed that they might be involved in longevity regulating and plant hormone signal transduction pathways. Further, using the gene expression profiles, we observed the negative regulation of miRNA-target pairs, such as MIR160~ARF18, MIR396~GRF6, MIR166~ATHB15/HD-ZIP, and MIR156/MIR157~SPL1. Interestingly, transcription factors such as WRKY, MYB, GAMYB, TCP4, and PIL1 were found to be regulated by the DEmiRs. The genes encoding PIL1 and RPS21C, regulated by upregulated miRNAs (e.g., egd-N-miR63-5p, egd-N-miR63-5p, and MIR169,) were downregulated exclusively in the dedifferentiation of E. camaldulensis. This is the first time to study the miRNA regulation in the dedifferentiation process of Eucalyptus and it will provide a valuable resource for future studies. More importantly, it will improve our understanding of miRNA regulation during the somatic embryogenesis of Eucalyptus and benefit the Eucalyptus breeding program.

De novo transcriptome analysis of white teak (Gmelina arborea Roxb) wood reveals critical genes involved in xylem development and secondary metabolism

BACKGROUND

Gmelina arborea Roxb is a fast-growing tree species of commercial importance for tropical countries due to multiple industrial uses of its wood. Wood is primarily composed of thick secondary cell walls of xylem cells which imparts the strength to the wood. Identification of the genes involved in the secondary cell wall biosynthesis as well as their cognate regulators is crucial to understand how the production of wood occurs and serves as a starting point for developing breeding strategies to produce varieties with improved wood quality, better paper pulping or new potential uses such as biofuel production. In order to gain knowledge on the molecular mechanisms and gene regulation related with wood development in white teak, a de novo sequencing and transcriptome assembly approach was used employing secondary cell wall synthesizing cells from young white teak trees.

RESULTS

For generation of transcriptome, RNA-seq reads were assembled into 110,992 transcripts and 49,364 genes were functionally annotated using plant databases; 5071 GO terms and 25,460 SSR markers were identified within xylem transcripts and 10,256 unigenes were assigned to KEGG database in 130 pathways. Among transcription factor families, C2H2, C3H, bLHLH and MYB were the most represented in xylem. Differential gene expression analysis using leaves as a reference was carried out and a total of 20,954 differentially expressed genes were identified including monolignol biosynthetic pathway genes. The differential expression of selected genes (4CL, COMT, CCoAOMT, CCR and NST1) was validated using qPCR.

CONCLUSIONS

We report the very first de novo transcriptome of xylem-related genes in this tropical timber species of commercial importance and constitutes a valuable extension of the publicly available transcriptomic resource aimed at fostering both basic and breeding studies.

Long Non-Coding RNA and Its Regulatory Network Response to Cold Stress in Eucalyptus urophylla S.T.Blake

Long non-coding RNA (lncRNA) plays an important regulatory role in plant growth and development, but its systematic identification and analysis in Eucalyptus has not yet been reported. Cold stress has a huge impact on the survival and yield of Eucalyptus seedlings, but the regulatory mechanism of lncRNA in Eucalyptus in response to cold stress is still unclear. In this study, the transcriptomes of young leaves of Eucalyptus urophylla S.T.Blake under low-temperature treatment and restoration were analyzed by RNA-seq. A total of 11,394 lncRNAs and 46,276 mRNAs were identified, of which 300 were differentially expressed lncRNAs (DE_lncRNAs) and 5606 were differentially expressed target genes of lncRNAs under cold stress, with the total number of target genes of DE_lncRNAs being 1681. A total of 677 differentially expressed transcription factors (TFs) were also identified, mainly including ERF, MYB and the NAC transcription factor family. Gene ontology (GO) analysis of the differentially expressed genes (DEGs) and target genes of DE_lncRNAs was mostly related to the response to cold stress and external stimuli. Furthermore, lncRNA–miRNA–mRNA regulatory networks were constructed, and 22 DE_lncRNAs were predicted to be targets or targeting mimics of 20 miRNAs. A qRT-PCR was used to verify the relative expression of genes in the regulatory EuGBF3-EUC_00002677-MSTRG.7690 network, and it matched the transcriptome data, indicating that it may play an important role in the response to cold stress in E. urophylla. This study provides a new insight into lncRNA and its regulatory network under abiotic stress, especially cold stress in E. urophylla.

Identification of multiple RNAs using feature fusion

Detection of novel transcripts with deep sequencing has increased the demand for computational algorithms as their identification and validation using in vivo techniques is time-consuming, costly and unreliable. Most of these discovered transcripts belong to non-coding RNAs, a large group known for their diverse functional roles but lacks the common taxonomy. Thus, upon the identification of the absence of coding potential in them, it is crucial to recognize their prime functional category. To address this heterogeneity issue, we divide the ncRNAs into three classes and present RNA classifier (RNAC) that categorizes the RNAs into coding, housekeeping, small non-coding and long non-coding classes. RNAC utilizes the alignment-based genomic descriptors to extract statistical, local binary patterns and histogram features and fuse them to construct the classification models with extreme gradient boosting. The experiments are performed on four species, and the performance is assessed on multiclass and conventional binary classification (coding versus no-coding) problems. The proposed approach achieved >93% accuracy on both classification problems and also outperformed other well-known existing methods in coding potential prediction. This validates the usefulness of feature fusion for improved performance on both types of classification problems. Hence, RNAC is a valuable tool for the accurate identification of multiple RNAs .

Biology and Function of miR159 in Plants

MicroR159 (miR159) is ancient, being present in the majority of land plants where it targets a class of regulatory genes called GAMYB or GAMYB-like via highly conserved miR159-binding sites. These GAMYB genes encode R2R3 MYB domain transcription factors that transduce the gibberellin (GA) signal in the seed aleurone and the anther tapetum. Here, GAMYB plays a conserved role in promoting the programmed cell death of these tissues, where miR159 function appears weak. By contrast, GAMYB is not involved in GA-signaling in vegetative tissues, but rather its expression is deleterious, leading to the inhibition of growth and development. Here, the major function of miR159 is to mediate strong silencing of GAMYB to enable normal growth. Highlighting this requirement of strong silencing are conserved RNA secondary structures associated with the miR159-binding site in GAMYB mRNA that promotes miR159-mediated repression. Although the miR159-GAMYB pathway in vegetative tissues has been implicated in a number of different functions, presently no conserved role for this pathway has emerged. We will review the current knowledge of the different proposed functions of miR159, and how this ancient pathway has been used as a model to help form our understanding of miRNA biology in plants.

Genome-wide analysis in Hevea brasiliensis laticifers revealed species-specific post-transcriptional regulations of several redox-related genes

MicroRNA-mediated post-transcriptional regulation has been reported on ROS production and scavenging systems. Although microRNAs first appeared highly conserved among plant species, several aspects of biogenesis, function and evolution of microRNAs were shown to differ. High throughput transcriptome and degradome analyses enable to identify small RNAs and their mRNA targets. A non-photosynthetic tissue particularly prone to redox reactions, laticifers from Hevea brasiliensis, revealed species-specific post-transcriptional regulations. This paper sets out to identify the 407 genes of the thirty main redox-related gene families harboured by the Hevea genome. There are 161 redox-related genes expressed in latex. Thirteen of these redox-related genes were targeted by 11 microRNAs. To our knowledge, this is the first report on a mutation in the miR398 binding site of the cytosolic CuZnSOD. A working model was proposed for transcriptional and post-transcriptional regulation with respect to the predicted subcellular localization of deduced proteins.

TaNTF2, a contributor for wheat resistance to the stripe rust pathogen

Nuclear Transport Factor 2 (NTF2) functions as a critical regulator in balancing the GTP-and GDP-bound forms of Ran, a class of evolutionarily conserved small GTP-binding protein. During the incompatible interaction between wheat-Puccinia striiformis f. sp. tritici (Pst), a cDNA fragment encoding a putative wheat NTF2 gene was found to be significantly induced, suggesting a potential role in wheat resistance to Pst. In this work, the full length of TaNTF2 was obtained, with three copies located on 7A, 7B and 7D chromosomes, respectively. QRT-PCR further verified the up-regulated expression of TaNTF2 in response to avirulent Pst. In addition, TaNTF2 was also induced by exogenous hormone applications, especially JA treatment. Transient expression of TaNTF2 in tobacco cells confirmed its subcellular localization in the cytoplasm, perinuclear area and nucleus. And virus induced gene silencing (VIGS) was used to identify the function of TaNTF2 during an incompatible wheat-Pst interaction. When TaNTF2 was knocked down, resistance of wheat to avirulentPst was decreased, with a bigger necrotic spots, and higher numbers of hyphal branches and haustorial mother cells. Our results demonstrated that TaNTF2 was a contributor for wheat resistance to the stripe rust pathogen, which will help to comprehensively understand the NTF2/Ran modulating mechanism in wheat-Pst interaction.