Genome-wide identification of conserved and novel microRNAs in one bud and two tender leaves of tea plant (Camellia sinensis) by small RNA sequencing, microarray-based hybridization and genome survey scaffold sequences

Identification and Characterisation of the CircRNAs Involved in the Regulation of Leaf Colour in Quercus mongolica

Circular RNAs (circRNAs) are important regulatory molecules involved in various biological processes. However, the potential function of circRNAs in the turning red process of Quercus mongolica leaves is unclear. This study used RNA-seq data to identify 6228 circRNAs in leaf samples from four different developmental stages and showed that 88 circRNAs were differentially expressed. A correlation analysis was performed between anthocyanins and the circRNAs. A total of 16 circRNAs that may be involved in regulating the colour of Mongolian oak leaves were identified. CircRNAs may affect the colour of Q. mongolica leaves by regulating auxin, cytokinin, gibberellin, ethylene, and abscisic acid. This study revealed the potential role of circRNAs in the colour change of Q. mongolica leaves.

Identification and Validation of the miR156 Family Involved in Drought Responses and Tolerance in Tea Plants (Camellia sinensis (L.) O. Kuntze)

The microRNA156 (miR156) family, one of the first miRNA families discovered in plants, plays various important roles in plant growth and resistance to various abiotic stresses. Previously, miR156s were shown to respond to drought stress, but miR156s in tea plants (Camellia sinensis (L.) O. Kuntze) have not been comprehensively identified and analyzed. Herein, we identify 47 mature sequences and 28 precursor sequences in tea plants. Our evolutionary analysis and multiple sequence alignment revealed that csn-miR156s were highly conserved during evolution and that the rates of the csn-miR156 members' evolution were different. The precursor sequences formed typical and stable stem-loop structures. The prediction of cis-acting elements in the CsMIR156s promoter region showed that the CsMIR156s had diverse cis-acting elements; of these, 12 CsMIR156s were found to be drought-responsive elements. The results of reverse transcription quantitative PCR (RT-qPCR) testing showed that csn-miR156 family members respond to drought and demonstrate different expression patterns under the conditions of drought stress. This suggests that csn-miR156 family members may be significantly involved in the response of tea plants to drought stress. Csn-miR156f-2-5p knockdown significantly reduced the Fv/Fm value and chlorophyll content and led to the accumulation of more-reactive oxygen species and proline compared with the control. The results of target gene prediction showed that csn-miR156f-2-5p targeted SQUAMOSA promoter binding protein-like (SPL) genes. Further analyses showed that CsSPL14 was targeted by csn-miR156f-2-5p, as confirmed through RT-qPCR, 5' RLM-RACE, and antisense oligonucleotide validation. Our results demonstrate that csn-miR156f-2-5p and CsSPL14 are involved in drought response and represent a new strategy for increasing drought tolerance via the breeding of tea plants.

Comparative genome-wide characterization of salt responsive micro RNA and their targets through integrated small RNA and de novo transcriptome profiling in sugarcane and its wild relative Erianthus arundinaceus

Soil salinity and saline irrigation water are major constraints in sugarcane affecting the production of cane and sugar yield. To understand the salinity induced responses and to identify novel genomic resources, integrated de novo transcriptome and small RNA sequencing in sugarcane wild relative, Erianthus arundinaceus salt tolerant accession IND 99-907 and salt-sensitive sugarcane genotype Co 97010 were performed. A total of 362 known miRNAs belonging to 62 families and 353 miRNAs belonging to 63 families were abundant in IND 99-907 and Co 97010 respectively. The miRNA families such as miR156, miR160, miR166, miR167, miR169, miR171, miR395, miR399, miR437 and miR5568 were the most abundant with more than ten members in both genotypes. The differential expression analysis of miRNA reveals that 221 known miRNAs belonging to 48 families and 130 known miRNAs belonging to 42 families were differentially expressed in IND 99-907 and Co 97010 respectively. A total of 12,693 and 7982 miRNA targets against the monoploid mosaic genome and a total of 15,031 and 12,152 miRNA targets against the de novo transcriptome were identified for differentially expressed known miRNAs of IND 99-907 and Co 97010 respectively. The gene ontology (GO) enrichment analysis of the miRNA targets revealed that 24, 12 and 14 enriched GO terms (FDR < 0.05) for biological process, molecular function and cellular component respectively. These miRNAs have many targets that associated in regulation of biotic and abiotic stresses. Thus, the genomic resources generated through this study are useful for sugarcane crop improvement through biotechnological and advanced breeding approaches.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03867-7.

CircRNAs responsive to winter dormancy and spring flushing conditions of tea leaf buds

Circular RNAs (circRNAs) are important regulators of diverse biological processes of plants. However, the evolution and potential functions of circRNAs during winter dormancy and spring bud flushing of tea plant is largely unknown. Using RNA-seq data, a total of 1184 circRNAs were identified in the winter dormant and spring bud flushing leaf samples of tea plants in two different cultivars exhibiting different duration of winter dormancy. A total of 156 circRNAs are found to be differentially expressed and the weighted gene co-expression network (WGCNA) analysis revealed that 22 and 20 differentially expressed-circRNAs (DE-circRNAs) positively correlated with the flushing and dormant leaf traits, respectively, in both the tea cultivars used. Some transcription factors (TFs) viz. MYB, WRKY, ERF, bHLH and several genes related to secondary metabolite biosynthetic pathways are found to co-express with circRNAs. DE-circRNAs also predicted to interact with miRNAs and can regulate phytohormone biosynthesis and various signalling pathways in tea plant. This study uncovers the potential roles of circRNAs to determine winter dormancy and spring bud flushing conditions in tea plants.

An Integrated Analysis of microRNAs and the Transcriptome Reveals the Molecular Mechanisms Underlying the Regulation of Leaf Development in Xinyang Maojian Green Tea (Camellia sinensis)

Xinyang Maojian (XYMJ) tea is one of the world's most popular green teas; the development of new sprouts directly affects the yield and quality of tea products, especially for XYMJ, which has hairy tips. Here, we used transcriptome and small RNA sequencing to identify mRNAs and miRNAs, respectively, involved in regulating leaf development in different plant tissues (bud, leaf, and stem). We identified a total of 381 conserved miRNAs. Given that no genomic data for XYMJ green tea are available, we compared the sequencing data for XYMJ green tea with genomic data from a closely related species (Tieguanyin) and the Camellia sinensis var. sinensis database; we identified a total of 506 and 485 novel miRNAs, respectively. We also identified 11 sequence-identical novel miRNAs in the tissues of XYMJ tea plants. Correlation analyses revealed 97 miRNA-mRNA pairs involved in leaf growth and development; the csn-miR319-2/csnTCP2 and miR159-csnMYB modules were found to be involved in leaf development in XYMJ green tea. Quantitative real-time PCR was used to validate the expression levels of the miRNAs and mRNAs. The miRNAs and target genes identified in this study might shed new light on the molecular mechanisms underlying the regulation of leaf development in tea plants.

Ribosome footprint profiling enables elucidating the systemic regulation of fatty acid accumulation in Acer truncatum

BACKGROUND

The accumulation of fatty acids in plants covers a wide range of functions in plant physiology and thereby affects adaptations and characteristics of species. As the famous woody oilseed crop, Acer truncatum accumulates unsaturated fatty acids and could serve as the model to understand the regulation and trait formation in oil-accumulation crops. Here, we performed Ribosome footprint profiling combing with a multi-omics strategy towards vital time points during seed development, and finally constructed systematic profiling from transcription to proteomes. Additionally, we characterized the small open reading frames (ORFs) and revealed that the translational efficiencies of focused genes were highly influenced by their sequence features.

RESULTS

The comprehensive multi-omics analysis of lipid metabolism was conducted in A. truncatum. We applied the Ribo-seq and RNA-seq techniques, and the analyses of transcriptional and translational profiles of seeds collected at 85 and 115 DAF were compared. Key members of biosynthesis-related structural genes (LACS, FAD2, FAD3, and KCS) were characterized fully. More meaningfully, the regulators (MYB, ABI, bZIP, and Dof) were identified and revealed to affect lipid biosynthesis via post-translational regulations. The translational features results showed that translation efficiency tended to be lower for the genes with a translated uORF than for the genes with a non-translated uORF. They provide new insights into the global mechanisms underlying the developmental regulation of lipid metabolism.

CONCLUSIONS

We performed Ribosome footprint profiling combing with a multi-omics strategy in A. truncatum seed development, which provides an example of the use of Ribosome footprint profiling in deciphering the complex regulation network and will be useful for elucidating the metabolism of A. truncatum seed oil and the regulatory mechanisms.

Genome-Wide Identification of 2-Oxoglutarate and Fe (II)-Dependent Dioxygenase (2ODD-C) Family Genes and Expression Profiles under Different Abiotic Stresses in Camellia sinensis (L.)

The 2-oxoglutarate and Fe (II)-dependent dioxygenase (2ODD-C) family of 2-oxoglutarate-dependent dioxygenases potentially participates in the biosynthesis of various metabolites under various abiotic stresses. However, there is scarce information on the expression profiles and roles of 2ODD-C genes in Camellia sinensis. We identified 153 Cs2ODD-C genes from C. sinensis, and they were distributed unevenly on 15 chromosomes. According to the phylogenetic tree topology, these genes were divided into 21 groups distinguished by conserved motifs and an intron/exon structure. Gene-duplication analyses revealed that 75 Cs2ODD-C genes were expanded and retained after WGD/segmental and tandem duplications. The expression profiles of Cs2ODD-C genes were explored under methyl jasmonate (MeJA), polyethylene glycol (PEG), and salt (NaCl) stress treatments. The expression analysis showed that 14, 13, and 49 Cs2ODD-C genes displayed the same expression pattern under MeJA and PEG treatments, MeJA and NaCl treatments, and PEG and NaCl treatments, respectively. A further analysis showed that two genes, Cs2ODD-C36 and Cs2ODD-C21, were significantly upregulated and downregulated after MeJA, PEG, and NaCl treatments, indicating that these two genes played positive and negative roles in enhancing the multi-stress tolerance. These results provide candidate genes for the use of genetic engineering technology to modify plants by enhancing multi-stress tolerance to promote phytoremediation efficiency.

Small RNAs and Karma methylation in Elaeis guineensis mother palms are linked to high clonal mantling

The mantled phenotype is an abnormal somaclonal variant arising from the oil palm cloning process and severe phenotypes lead to oil yield losses. Hypomethylation of the Karma retrotransposon within the B-type MADS-box EgDEF1 gene has been associated with this phenotype. While abnormal Karma-EgDEF1 hypomethylation was detected in mantled clones, we examined the methylation state of Karma in ortets that gave rise to high mantling rates in their clones. Small RNAs (sRNAs) were proposed to play a role in Karma hypomethylation as part of the RNA-directed DNA methylation process, hence differential expression analysis of sRNAs between the ortet groups was conducted. While no sRNA was differentially expressed at the Karma-EgDEF1 region, three sRNA clusters were differentially regulated in high-mantling ortets. The first two down-regulated clusters were possibly derived from long non-coding RNAs while the third up-regulated cluster was derived from the intron of a DnaJ chaperone gene. Several predicted mRNA targets for the first two sRNA clusters conversely displayed increased expression in high-mantling relative to low-mantling ortets. These predicted mRNA targets may be associated with defense or pathogenesis response. In addition, several differentially methylated regions (DMRs) were identified in Karma and its surrounding regions, mainly comprising subtle CHH hypomethylation in high-mantling ortets. Four of the 12 DMRs were located in a region corresponding to hypomethylated areas at the 3'end of Karma previously reported in mantled clones. Further investigations on these sRNAs and DMRs may indicate the predisposition of certain ortets towards mantled somaclonal variation.

Sequence, Secondary Structure, and Phylogenetic Conservation of MicroRNAs in Arabidopsis thaliana

MicroRNAs are small non-coding RNA molecules that are produced in a cell endogenously. They are made up of 18 to 26 nucleotides in strength. Due to their evolutionary conserved nature, most of the miRNAs provide a logical basis for the prediction of novel miRNAs and their clusters in plants such as sunflowers related to the Asteraceae family. In addition, they participate in different biological processes of plants, including cell signaling and metabolism, development, growth, and tolerance to (biotic and abiotic) stresses. In this study profiling, conservation and characterization of novel miRNA possessing conserved nature in various plants and their targets annotation in sunflower (Asteraceae) were obtained by using various computational tools and software. As a result, we looked at 152 microRNAs in Arabidopsis thaliana that had already been predicted. Drought tolerance stress is mediated by these 152 non-coding RNAs. Following that, we used local alignment to predict novel microRNAs that were specific to Helianthus annuus. We used BLAST to do a local alignment, and we chose sequences with an identity of 80% to 100%. MIR156a, MIR164a, MIR165a, MIR170, MIR172a, MIR172b, MIR319a, MIR393a, MIR394a, MIR399a, MIR156h, and MIR414 are the new anticipated miRNAs. We used MFold to predict the secondary structure of new microRNAs. We used conservation analysis and phylogenetic analysis against a variety of organisms, including Gossypium hirsutum, H. annuus, A. thaliana, Triticum aestivum, Saccharum officinarum, Zea mays, Brassica napus, Solanum tuberosum, Solanum lycopersicum, and Oryza sativa, to determine the evolutionary history of these novel non-coding RNAs. Clustal W was used to analyze the evolutionary history of discovered miRNAs.

Revealing the mechanisms of the bioactive ingredients accumulation in Polygonatum cyrtonema by multiomics analyses

Polygonatum cyrtonema is a medicinal and edible herb rich in polysaccharides, steroidal saponins, and flavonoids that has been widely used as a food, vegetable, and medicine over the years. Although previous studies have preliminarily explored the metabolic and transcriptional regulatory mechanisms of the main secondary metabolites in P. cyrtonema, the complex mechanism of microRNA (miRNA)-mediated posttranscriptional regulation remains unclear. Metabolome analysis showed that iso-ophiopogonanone B, (25S)-pratioside D1, disporopsin, and isodiosgenin-Glc-Glc, which are associated with intermediates in the flavonoids and saponins pathways, were significantly upregulated in the stem and leaf compared with the rhizome, and most saccharides, including arabinose, cellobiose, maltotetraose, and panose, showed the opposite trend, suggesting that they may contribute to the formation and accumulation of the main active ingredients in P. cyrtonema. We found that 4-hydroxymandelonitrile have a relatively good inhibitory effect on α-glucosidase, indicating that it may play a role in hypoglycemic functions. Transcriptome and weighted gene coexpression network analysis (WGCNA) were combined to reveal several candidate genes involved in the accumulation of polysaccharides, saponins, and flavonoids, including PcSQLE, PcCYP71A1, PcSUS, PcFK, and PcMYB102. Integrated analyses of miRNAs and messengerRNAs (mRNAs) showed that novel_miR14, novel_miR49, novel_miR75, and aof_miR164 were negatively correlated with alpha-linolenic acid metabolism and the mitogen activated protein kinase (MAPK) signaling pathway, including PcAOS, PcSPLA2, PcFRK1, and PcDELLA, indicating that these miRNAs may coordinately regulate the biosynthesis of other secondary metabolites in P. cyrtonema. These findings will facilitate in-depth research on the functions of these miRNAs and mRNAs related to the main active substances for pathological and biological regulation, which will be beneficial to provide theoretical guidance for the molecular breeding of P. cyrtonema.

Comparative transcriptome analysis reveals key pathways and genes involved in trichome development in tea plant (Camellia sinensis)

Trichomes, which develop from epidermal cells, are considered one of the important characteristics of the tea plant [Camellia sinensis (L.) O. Kuntze]. Many nutritional and metabolomic studies have indicated the important contributions of trichomes to tea products quality. However, understanding the regulation of trichome formation at the molecular level remains elusive in tea plants. Herein, we present a genome-wide comparative transcriptome analysis between the hairless Chuyeqi (CYQ) with fewer trichomes and the hairy Budiaomao (BDM) with more trichomes tea plant genotypes, toward the identification of biological processes and functional gene activities that occur during trichome development. In the present study, trichomes in both cultivars CYQ and BDM were unicellular, unbranched, straight, and soft-structured. The density of trichomes was the highest in the bud and tender leaf periods. Further, using the high-throughput sequencing method, we identified 48,856 unigenes, of which 31,574 were differentially expressed. In an analysis of 208 differentially expressed genes (DEGs) encoding transcription factors (TFs), five may involve in trichome development. In addition, on the basis of the Gene Ontology (GO) annotation and the weighted gene co-expression network analysis (WGCNA) results, we screened several DEGs that may contribute to trichome growth, including 66 DEGs related to plant resistance genes (PRGs), 172 DEGs related to cell wall biosynthesis pathway, 29 DEGs related to cell cycle pathway, and 45 DEGs related to cytoskeleton biosynthesis. Collectively, this study provided high-quality RNA-seq information to improve our understanding of the molecular regulatory mechanism of trichome development and lay a foundation for additional trichome studies in tea plants.

An insight into microRNA biogenesis and its regulatory role in plant secondary metabolism

The present review highlights the regulatory roles of microRNAs in plant secondary metabolism and focuses on different bioengineering strategies to modulate secondary metabolite content in plants. MicroRNAs (miRNAs) are the class of small endogenous, essential, non-coding RNAs that riboregulate the gene expression involved in various biological processes in most eukaryotes. MiRNAs has emerged as important regulators in plants that function by silencing target genes through cleavage or translational inhibition. These miRNAs plays an important role in a wide range of plant biological and metabolic processes, including plant development and various environmental response controls. Several important plant secondary metabolites like alkaloids, terpenoids, and phenolics are well studied for their function in plant defense against different types of pests and herbivores. Due to the presence of a wide range of biological and pharmaceutical properties of plant secondary metabolites, it is important to study the regulation of their biosynthetic pathways. The contribution of miRNAs in regulating plant secondary metabolism is not well explored. Recent advancements in molecular techniques have improved our knowledge in understanding the molecular function of genes, proteins, enzymes, and small RNAs involved in different steps of secondary metabolic pathways. In the present review, we have discussed the recent progress made on miRNA biogenesis, its regulation, and highlighted the current research developed in the field of identification, analysis, and characterizations of various miRNAs that regulate plant secondary metabolism. We have also discussed how different bioengineering strategies such as artificial miRNA (amiRNA), endogenous target mimicry, and CRISPR/Cas9 could be utilized to enhance the secondary metabolite production in plants.

Global DNA Methylation and mRNA-miRNA Variations Activated by Heat Shock Boost Early Microspore Embryogenesis in Cabbage (Brassica oleracea)

Microspore culture, a type of haploid breeding, is extensively used in the cultivation of cruciferous crops such as cabbage. Heat shock (HS) treatment is essential to improve the embryo rate during the culture process; however, its molecular role in boosting early microspore embryogenesis (ME) remains unknown. Here we combined DNA methylation levels, miRNAs, and transcriptome profiles in isolated microspores of cabbage ‘01-88’ under HS (32 °C for 24 h) and normal temperature (25 °C for 24 h) to investigate the regulatory roles of DNA methylation and miRNA in early ME. Global methylation levels were significantly different in the two pre-treatments, and 508 differentially methylated regions (DMRs) were identified; 59.92% of DMRs were correlated with transcripts, and 39.43% of miRNA locus were associated with methylation levels. Significantly, the association analysis revealed that 31 differentially expressed genes (DEGs) were targeted by methylation and miRNA and were mainly involved in the reactive oxygen species (ROS) response and abscisic acid (ABA) signaling, indicating that HS induced DNA methylation, and miRNA might affect ME by influencing ROS and ABA. This study revealed that DNA methylation and miRNA interfered with ME by modulating key genes and pathways, which could broaden our understanding of the molecular regulation of ME induced by HS pre-treatment.

Genome-Wide Investigation of the MiR166 Family Provides New Insights into Its Involvement in the Drought Stress Responses of Tea Plants (Camellia sinensis (L.) O. Kuntze)

MicroRNA166 (miR166) is a highly conserved plant miRNA that plays a crucial role in plant growth and the resistance to various abiotic stresses. However, the miR166s in tea (Camellia sinensis (L.) O. Kuntze) have not been comprehensively identified and analyzed. This study identified 30 mature miR166s and twelve pre-miR166s in tea plants. An evolutionary analysis revealed that csn-miR166s originating from the 3′ arm of their precursors were more conserved than the csn-miR166s derived from the 5′ arm of their precursors. The twelve pre-miR166s in tea were divided into two groups, with csn-MIR166 Scaffold364-2 separated from the other precursors. The Mfold-based predictions indicated that the twelve csn-MIR166s formed typical and stable structures comprising a stem-loop hairpin, with minimum free energy ranging from −110.90 to −71.80 kcal/mol. An analysis of the CsMIR166 promoters detected diverse cis-acting elements, including those related to light responses, biosynthesis and metabolism, abiotic stress defenses, and hormone responses. There was no one-to-one relationship between the csn-miR166s and their targets, but most csn-miR166s targeted HD-Zip III genes. Physiological characterization of tea plants under drought stress showed that leaf water content proportionally decreased with the aggravation of drought stress. In contrast, tea leaves’ malondialdehyde (MDA) content proportionally increased. Moreover, the cleavage site of the ATHB-15-like transcript was identified according to a modified 5′ RNA ligase-mediated rapid amplification of cDNA ends. The RT-qPCR data indicated that the transcription of nine csn-miR166s was negatively correlated with their target gene.

Identification of Novel miRNAs and Their Target Genes in the Response to Abscisic Acid in Arabidopsis

miRNAs are involved in various biological processes, including adaptive responses to abiotic stress. To understand the role of miRNAs in the response to ABA, ABA-responsive miRNAs were identified by small RNA sequencing in wild-type Arabidopsis, as well as in abi1td, mkkk17, and mkkk18 mutants. We identified 10 novel miRNAs in WT after ABA treatment, while in abi1td, mkkk17, and mkkk18 mutants, three, seven, and nine known miRNAs, respectively, were differentially expressed after ABA treatment. One novel miRNA (miRn-8) was differentially expressed in the mkkk17 mutant. Potential target genes of the miRNA panel were identified using psRNATarget. Sequencing results were validated by quantitative RT-PCR of several known and novel miRNAs in all genotypes. Of the predicted targets of novel miRNAs, seven target genes of six novel miRNAs were further validated by 5' RLM-RACE. Gene ontology analyses showed the potential target genes of ABA-responsive known and novel miRNAs to be involved in diverse cellular processes in plants, including development and stomatal movement. These outcomes suggest that a number of the identified miRNAs have crucial roles in plant responses to environmental stress, as well as in plant development, and might have common regulatory roles in the core ABA signaling pathway.

Analysis of miRNA-mediated regulation of flowering induction in Lilium × formolongi

BACKGROUND

MicroRNAs play pivotal roles in plant vegetative phase change and flowering induction via integrating into multiple flowering pathways. Lilium × formolongi is an important ornamental lily cultivar that can flower within one year after sowing. However, it remains unresolved how miRNA-mediated regulation networks contribute to the L. × formolongi characteristics of a short vegetative growth period and rapid flowering.

RESULTS

In this study, the small RNA libraries and one degradome library were constructed for L. × formolongi during vegetative growth and flowering initiation, and 366 conserved miRNAs and 32 novel miRNAs were identified. Additionally, 84 miRNAs were significantly differentially expressed during development. A total of 396 targets of 185 miRNAs were identified and validated through degradome sequencing. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that functions of the targets were top enriched in the cold and cadmium ion responses, pentose phosphate pathway and carbon fixation in photosynthetic organisms. Furthermore, among 23 differentially expressed miRNA-target pairs, the miR156s-LfSPL2, miR172a-LfAP2 and miR164a-LfNAC pairs as well as miR159a-LfSPL2 were found to be relevant to flowering based on the correlation analysis of expression profiles in the miRNA libraries, degradome and transcriptome. A coexpression regulatory network focused on differentially expressed pairs was also constructed by WGCNA, and 14 miRNAs were considered putative key miRNAs during vegetative development and flowering induction. miR156a/ d/ e showed particularly strong relationships with other miRNAs in the coexpression network.

CONCLUSIONS

This study provides cues for the further exploration of the regulatory mechanisms of short vegetative development and flowering in L. × formolongi.

Understanding the role of miRNAs for improvement of tea quality and stress tolerance

MicroRNAs (miRNAs) are an emerging class of small non-coding RNAs that exhibit important role in regulation of gene expression, mostly through the mechanism of cleavage and/or inhibition of translation of target mRNAs during or after transcription. Although much has been unravelled about the role of miRNAs in diverse biological processes like maintenance of functional integrity of genes and genome, growth and development, metabolism, and adaptive responses towards biotic and abiotic stresses in plants, not much is known on their specific roles in majority of cash crops - an area of investigation with potentially significant and gainful economic implications. Tea (Camellia sinensis) is globally the second most consumed beverage after water and its cultivation has major agro-economic and social ramifications. In recent years, global tea production has been greatly challenged by many biotic and abiotic stress factors and a deeper understanding of molecular processes regulating stress adaptation in this largely under investigated crop stands to significantly facilitate potential crop improvement strategies towards durable stress tolerance. This review endeavours to highlight recent advances in our understanding of the role of miRNAs in regulating stress tolerance traits in tea plant with additional focus on their role in determining tea quality attributes.

Characterizing microRNAs and their targets in different organs of Camellia sinensis var. assamica

To comprehensively annotate miRNAs and their targets in tea plant, Camellia sinensis, we sequenced small and messenger RNAs of 9 samples of Camellia sinensis var. assamica (YK-10), a diploid elite cultivar widely grown in southwest China. In order to identify targets of miRNAs, we sequenced two degradome sequencing profiles from leaves and roots of YK-10, respectively. By analyzing the small RNA-Seq profiles, we newly identified 137 conserved miRNAs and 23 species specific miRNAs in the genome of YK-10, which significantly improved the annotation of miRNAs in tea plant. Approximately 2000 differently expressed genes were identified when comparing RNA-Seq profiles of any two of the three organs selected in the study. Totally, more than 5000 targets of conserved miRNAs were identified in the two degradome profiles. Furthermore, our results suggest that a few miRNAs play roles in the biosynthesis pathways of theanine, caffeine and flavonoid. These results enhance our understanding of small RNA guided gene regulations in different organs of tea plant.

Transcriptomic and Translatomic Analyses Reveal Insights into the Developmental Regulation of Secondary Metabolism in the Young Shoots of Tea Plants (Camellia sinensis L.)

Accumulation of secondary metabolites in the young shoots of tea plants is developmentally modulated, especially flavonoids. Here, we investigate the developmental regulation mechanism of secondary metabolism in the developing leaves of tea plants using an integrated multiomic approach. For the pair of Leaf2/Bud, the correlation coefficient of the fold change of mRNA and RPFs abundances involved in flavonoid biosynthesis was 0.9359, being higher than that of RPFs and protein (R² = 0.6941). These correlations were higher than the corresponding correlation coefficients for secondary metabolisms and genome-wide scale. Metabolomic analysis demonstrates that the developmental modulations of the structural genes for flavonoid biosynthesis-related pathways align with the concentration changes of catechin and flavonol glycoside groups. Relatively high translational efficiency (TE > 2) was observed in the four flavonoid structural genes (chalcone isomerase, dihydroflavonol 4-reductase, anthocyanidin synthase, and flavonol synthase). In addition, we originally provided the information on identified small open reading frames (small ORFs) and main ORFs in tea leaves and elaborated that the presence of upstream ORFs may have a repressive effect on the translation of downstream ORFs. Our data suggest that transcriptional regulation coordinates with translational regulation and may contribute to the elevation of translational efficiencies for the structural genes involved in the flavonoid biosynthesis pathways during tea leaf development.

Systematic Identification and Analysis of Light-Responsive Circular RNA and Co-expression Networks in Lettuce (Lactuca sativa)

Circular RNA (circRNA) is a covalently-closed single-stranded RNA molecule that plays an important role in transcriptional regulation of gene expression in a variety of species. Light intensity is a pivotal environmental factor affecting plant growth and development. However, little is known regarding photoresponsive plant circRNAs. Here, we aimed to investigate the expression and function of circRNAs in lettuce leaves in response to different light intensity treatments. We performed RNA sequencing (RNA-Seq) on leaves of lettuce (Lactuca sativa) to determine circRNA expression profiles and reverse-transcription polymerase chain reaction (PCR) to validate the candidate circRNA molecules. We then combined bioinformatics approach to explore the function of the parental genes of circRNA, including network, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes pathway analysis. We identified 1650 circRNAs in lettuce, of which 1508 (86.40%) were derived from exons. Using real-time PCR, we characterized 10 validated differentially expressed circRNAs and their parental genes, all of which showed expression patterns consistent with RNA-Seq data. Interestingly, the expression of circRNA was, in some cases, inversely correlated with the expression of the parental gene. Furthermore, analysis of the circRNA-microRNA-mRNA network suggests that circRNAs may be involved in plant hormone signaling and chlorophyll metabolism during photoreactivity. These findings provide an essential reference basis for studying circRNAs' biological mechanisms in light-treated plants.

Utilization of microRNAs and their regulatory functions for improving biotic stress tolerance in tea plant [Camellia sinensis (L.) O. Kuntze]

MicroRNAs play a central role in responses to biotic stressors through their interactions with their target mRNAs. Tea plant (Camellia sinensis L.), an important beverage crop, is vulnerable to tea geometrid and anthracnose disease that causes considerable crop loss and tea production worldwide. Sustainable production of tea in the current scenario to biotic factors is major challenges. To overcome the problem of biotic stresses, high-throughput sequencing (HTS) with bioinformatics analyses has been used as an effective approach for the identification of stress-responsive miRNAs and their regulatory functions in tea plant. These stress-responsive miRNAs can be utilized for miRNA-mediated gene silencing to enhance stress tolerance in tea plant. Therefore, this review summarizes the current understanding of miRNAs regulatory functions in tea plant responding to Ectropis oblique and Colletotrichum gloeosporioides attacks for future miRNA research. Also, it highlights the utilization of miRNA-mediated gene silencing strategies for developing biotic stress-tolerant tea plant.

Mining of miRNAs from EST data in Dendrobium nobile

Dendrobium nobile is an orchid species highly popular for its therapeutic properties and is often used as a medicinal herb. Documenting miRNA-target associations in D. nobile is an important step to facilitate functional genomics studies in this species. Therefore, it is of interest to identify miRNA sequences from EST data available in public databases using known techniques and tools. We report 14 potential miRNAs from three ESTs of D. nobile. They belong to 3 miRNA families (miR390, miR528 and miR414) linking to transcription factor regulation, signal transduction, DNA and protein binding, and various cellular processes covering 34 different metabolic networks in KEGG. These results help in the understanding of miRNA-mRNAs functional networks in Dendrobium nobile.

Identification of Grafting-Responsive MicroRNAs Associated with Growth Regulation in Pecan [Carya illinoinensis (Wangenh.) K. Koch]

Pecan [Carya illinoinensis (Wangenh.) K. Koch] is an economically important nut tree and grafting is often used for clonal propagation of cultivars. However, there is a lack of research on the effects of rootstocks on scions, which are meaningful targets for directed breeding of pecan grafts. MicroRNAs (miRNAs) play an important role in many biological processes, but the mechanism underlying the involvement of miRNAs in grafting-conferred physiological changes is unclear. To identify the grafting-responsive miRNAs that may be involved in the regulation of growth in grafted pecan, six small RNA libraries were constructed from the phloem of two groups of grafts with significantly different growth performance on short and tall rootstocks. A total of 441 conserved miRNAs belonging to 42 miRNA families and 603 novel miRNAs were identified. Among the identified miRNAs, 24 (seven conserved and 17 novel) were significantly differentially expressed by the different grafts, implying that they might be responsive to grafting and potentially involved in the regulation of graft growth. Ninety-five target genes were predicted for the differentially expressed miRNAs; gene annotation was available for 33 of these. Analysis of their targets suggested that the miRNAs may regulate auxin transport, cell activity, and inorganic phosphate (Pi) acquisition, and thereby, mediate pecan graft growth. Use of the recently-published pecan genome enabled identification of a substantial population of miRNAs, which are now available for further research. We also identified the grafting-responsive miRNAs and their potential roles in pecan graft growth, providing a basis for research on long-distance regulation in grafted pecan.

Transgenerational Effects of Water-Deficit and Heat Stress on Germination and Seedling Vigour-New Insights from Durum Wheat microRNAs

Water deficiency and heat stress can severely limit crop production and quality. Stress imposed on the parents during reproduction could have transgenerational effects on their progeny. Seeds with different origins can vary significantly in their germination and early growth. Here, we investigated how water-deficit and heat stress on parental durum wheat plants affected seedling establishment of the subsequent generation. One stress-tolerant and one stress-sensitive Australian durum genotype were used. Seeds were collected from parents with or without exposure to stress during reproduction. Generally, stress on the previous generation negatively affected seed germination and seedling vigour, but to a lesser extent in the tolerant variety. Small RNA sequencing utilising the new durum genome assembly revealed significant differences in microRNA (miRNA) expression in the two genotypes. A bioinformatics approach was used to identify multiple miRNA targets which have critical molecular functions in stress adaptation and plant development and could therefore contribute to the phenotypic differences observed. Our data provide the first confirmation of the transgenerational effects of reproductive-stage stress on germination and seedling establishment in durum wheat. New insights gained on the epigenetic level indicate that durum miRNAs could be key factors in optimising seed vigour for breeding superior germplasm and/or varieties.

Integrated Analysis of Long Non-coding RNAs (lncRNAs) and mRNAs Reveals the Regulatory Role of lncRNAs Associated With Salt Resistance in Camellia sinensis

Tea plant (Camellia sinensis), an important economic crop, is seriously affected by various abiotic stresses, including salt stress, which severely diminishes its widespread planting. However, little is known about the roles of long non-coding RNAs (lncRNAs) in transcriptional regulation under salt stress. In this study, high-throughput sequencing of tea shoots under salt-stress and control conditions was performed. Through sequencing analysis, 16,452 unique lncRNAs were identified, including 172 differentially expressed lncRNAs (DE-lncRNAs). The results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of their cis- and trans-target genes showed that these DE-lncRNAs play important roles in many pathways such as the galactinol synthase (GOLS), calcium signaling pathway, and interact with transcription factors (TFs) under salt stress. The data from the gene-specific antisense oligodeoxynucleotide-mediated reduction in the lncRNA MSTRG.139242.1 and its predicted interacting gene, TEA027212.1 (Ca²⁺-ATPase 13), in tea leaves revealed that MSTRG.139242.1 may function in the response of tea plants to high salinity. In addition, 12 lncRNAs were predicted to be target mimics of 17 known mature miRNAs, such as miR156, that are related to the salt-stress response in C. sinensis. Our results provide new insights into lncRNAs as ubiquitous regulators in response to salt stress in tea plants.

Transcriptome-wide miRNA identification of Bacopa monnieri: a cross-kingdom approach

Bacopa monnieri known as 'Brahmi' is a well-known medicinal plant belonging to Scrophulariaceae family for its nootropic properties. To the best of our knowledge, no characterization data is available on the potential role of micro RNAs (miRNAs) from this plant till date. We present here the first report of computational characterizations of miRNAs from B. monnieri. Owing to the high conservation of miRNAs in nature, new and potential miRNAs can be identified in plants using in silico techniques. Using the plant miRNA sequences present in the miRBase repository, a total of 12 miRNAs were identified from B. monnieri which pertained to 11 miRNA families from the shoot and root transcriptome data. Furthermore, gene ontology analysis of the identified 68 human target genes exhibited significance in various biological processes. These human target genes were associated with signaling pathways like NF-kB and MAPK with TRAF2, CBX1, IL1B, ITGA4 and ITGB1BP1 as the top five hub nodes. This cross-kingdom study provides initial insights about the potential of miRNA-mediated cross-kingdom regulation and unravels the essential target genes of human with implications in numerous human diseases including cancer.

Tea plant genomics: achievements, challenges and perspectives

Tea is among the world's most widely consumed non-alcoholic beverages and possesses enormous economic, health, and cultural values. It is produced from the cured leaves of tea plants, which are important evergreen crops globally cultivated in over 50 countries. Along with recent innovations and advances in biotechnologies, great progress in tea plant genomics and genetics has been achieved, which has facilitated our understanding of the molecular mechanisms of tea quality and the evolution of the tea plant genome. In this review, we briefly summarize the achievements of the past two decades, which primarily include diverse genome and transcriptome sequencing projects, gene discovery and regulation studies, investigation of the epigenetics and noncoding RNAs, origin and domestication, phylogenetics and germplasm utilization of tea plant as well as newly developed tools/platforms. We also present perspectives and possible challenges for future functional genomic studies that will contribute to the acceleration of breeding programs in tea plants.

MicroRNAs and their targeted genes associated with phase changes of stem explants during tissue culture of tea plant

Elucidation of the molecular mechanism related to the dedifferentiation and redifferentiation during tissue culture will be useful for optimizing regeneration system of tea plant. In this study, an integrated sRNAome and transcriptome analyses were carried out during phase changes of the stem explant culture. Among 198 miRNAs and 8001 predicted target genes, 178 differentially expressed miRNAs and 4264 potential targets were screened out from explants, primary calli, as well as regenerated roots and shoots. According to KEGG analysis of the potential targets, pathway of "aminoacyl-tRNA biosynthesis", "proteasome" and "glutathione metabolism" was of great significance during the dedifferentiation, and pathway of "porphyrin and chlorophyll metabolism", "mRNA surveillance pathway", "nucleotide excision repair" was indispensable for redifferentiation of the calli. Expression pattern of 12 miRNAs, including csn-micR390e, csn-miR156b-5p, csn-miR157d-5p, csn-miR156, csn-miR166a-3p, csn-miR166e, csn-miR167d, csn-miR393c-3p, csn-miR394, csn-miR396a-3p, csn-miR396 and csn-miR396e-3p, was validated by qRT-PCR among 57 differentially expressed phase-specific miRNAs. Validation also confirmed that regulatory module of csn-miR167d/ERF3, csn-miR156/SPB1, csn-miR166a-3p/ATHB15, csn-miR396/AIP15A, csn-miR157d-5p/GST and csn-miR393c-3p/ATG18b might play important roles in regulating the phase changes during tissue culture of stem explants.

Integrated analysis of miRNAs and their targets reveals that miR319c/TCP2 regulates apical bud burst in tea plant (Camellia sinensis)

MAIN CONCLUSION

The roles of microRNA-mediated epigenetic regulation were highlighted in the bud dormancy-activity cycle, implying that certain differentially expressed miRNAs play crucial roles in apical bud burst, such as csn-miR319c/TCP2. microRNAs (miRNAs) are a class of small non-coding RNAs that regulate gene expression by targeting mRNA transcripts for cleavage or directing translational inhibition. To investigate whether miRNAs regulate bud dormancy-activation transition in tea plant, which largely affects the yield and price of tea products and adaptability of tea trees, we constructed small RNA libraries from three different periods of bud dormancy-burst transition. Through sequencing analysis, 262 conserved and 83 novel miRNAs were identified, including 118 differentially expressed miRNAs. Quantitative RT-PCR results for randomly selected miRNAs exhibited that our comprehensive analysis is highly reliable and accurate. The content of caffeine increased continuously from the endodormancy bud to flushing bud, and differentially expressed miRNAs coupling with their targets associated with bud burst were identified. Remarkably, csn-miR319c was downregulated significantly from the quiescent bud to burst bud, while its target gene CsnTCP2 (TEOSINTE BRANCHED/CYCLOIDEA/PROLIFERATING CELL FACTOR 2) displayed opposite expression patterns. Co-transformation experiment in tobacco demonstrated that csn-miR319c can significantly suppress the functions of CsnTCP2. This study on miRNAs and the recognition of target genes could provide new insights into the molecular mechanism of the bud dormancy-activation transition in tea plant.

Long noncoding RNAs and miRNAs regulating terpene and tartaric acid biosynthesis in rose-scented geranium

This study aimed to explore the noncoding RNAs, which have emerged as key regulatory molecules in biological processes, in rose-scented geranium. We analyzed RNA-seq data revealing 26 784 long noncoding RNAs (lncRNAs) and 871 miRNAs in rose-scented geranium. A total of 466 lncRNAs were annotated using different plant lncRNA public databases. Furthermore, 372 lncRNAs and 99 miRNAs were detected that target terpene and tartarate biosynthetic pathways. An interactome, comprising of lncRNAs, miRNAs, and mRNAs, was constructed that represents a noncoding RNA regulatory network of the target mRNAs. Real-time quantitative PCR expression validation was done for selected lncRNAs involved in the regulation of terpene and tartaric acid pathways. This study provides the first insights into the regulatory functioning of noncoding RNAs in rose-scented geranium.

Uncovering anthocyanin biosynthesis related microRNAs and their target genes by small RNA and degradome sequencing in tuberous roots of sweetpotato

BACKGROUND

Compared with white-fleshed sweetpotato (WFSP), purple-fleshed sweetpotato (PFSP) is a desirable resource for functional food development because of the abundant anthocyanin accumulation in its tuberous roots. Some studies have shown that the expression regulation mediated by miRNA plays an important role in anthocyanin biosynthesis in plants. However, few miRNAs and their corresponding functions related to anthocyanin biosynthesis in tuberous roots of sweetpotato have been known.

RESULTS

In this study, small RNA (sRNA) and degradome libraries from the tuberous roots of WFSP (Xushu-18) and PFSP (Xuzishu-3) were constructed, respectively. Totally, 191 known and 33 novel miRNAs were identified by sRNA sequencing, and 180 target genes cleaved by 115 known ib-miRNAs and 5 novel ib-miRNAs were identified by degradome sequencing. Of these, 121 miRNAs were differently expressed between Xushu-18 and Xuzishu-3. Integrated analysis of sRNA, degradome sequencing, GO, KEGG and qRT-PCR revealed that 26 differentially expressed miRNAs and 36 corresponding targets were potentially involved in the anthocyanin biosynthesis. Of which, an inverse correlation between the expression of ib-miR156 and its target ibSPL in WFSP and PFSP was revealed by both qRT-PCR and sRNA sequencing. Subsequently, ib-miR156 was over-expressed in Arabidopsis. Interestingly, the ib-miR156 over-expressing plants showed suppressed abundance of SPL and a purplish phenotype. Concomitantly, upregulated expression of four anthocyanin pathway genes was detected in transgenic Arabidopsis plants. Finally, a putative ib-miRNA-target model involved in anthocyanin biosynthesis in sweetpotato was proposed.

CONCLUSIONS

The results represented a comprehensive expression profiling of miRNAs related to anthocyanin accumulation in sweetpotato and provided important clues for understanding the regulatory network of anthocyanin biosynthesis mediated by miRNA in tuberous crops.

Integrated analysis of mRNA-seq and miRNA-seq in calyx abscission zone of Korla fragrant pear involved in calyx persistence

BACKGROUND

The objective of this study was to characterize molecular mechanism of calyx persistence in Korla fragrant pear by transcriptome and small RNA sequencing. Abscission zone tissues of flowers at three stages (the first, fifth and ninth days of the late bloom stage), with 50 mg/L GA₃ (calyx persistence treatment, C_1, C_5, C_9) or 500 mg/L PP₃₃₃ (calyx abscission treatment, T_1, T_5, T_9), were collected and simultaneously conducted transcriptome and small RNA sequencing.

RESULTS

Through association analysis of transcriptome and small RNA sequencing, mRNA-miRNA network was conducted. Compared calyx persistence groups with calyx abscission groups during the same stage, 145, 56 and 150 mRNA-miRNA pairs were obtained in C_1 vs T_1, C_5 vs T_5 and C_9 vs T_9, respectively; When C_1 compared with C_5 and C_9, 90 and 506 mRNA-miRNA pairs were screened respectively, and 255 mRNA-miRNA pairs were obtained from the comparison between C_5 and C_9; When T_1 compared with the T_5 and T_9, respectively, 206 and 796 mRNA-miRNA pairs were obtained, and 383 mRNA-miRNA pairs were obtained from the comparison between T_5 and T_9. These mRNAs in miRNA-mRNA pairs were significantly enriched into the terpenoid backbone biosynthesis, photosynthesis - antenna proteins, porphyrin and chlorophyll metabolism, carotenoid biosynthesis, zeatin biosynthesis and plant hormone signal transduction. In addition, we obtained some key genes from miRNA-mRNA pairs that may be associated with calyx abscission, including protein phosphatase 2C (psi-miR394a-HAB1), receptor-like protein kinase (psi-miR396a-5p-HERK1), cellulose synthase-like protein D3 (psi-miR827-CSLD3), beta-galactosidase (psi-miR858b-β-galactosidase), SPL-psi-miR156j/157d, abscisic acid 8'-hydroxylase 1 (psi-miR396a-5p-CYP707A1) and auxin response factor (psi-miR160a-3p-ARF6, psi-miR167d-ARF18, psi-miR167a-5p-ARF25), etc. CONCLUSION: By integrated analysis mRNA and miRNA, our study gives a better understanding of the important genes and regulation pathway related to calyx abscission in Korla fragrant pear. We have also established the network of miRNA-mRNA pairs to learn about precise regulation of miRNA on calyx abscission.

Circular RNA architecture and differentiation during leaf bud to young leaf development in tea (Camellia sinensis)

Circular RNA (circRNA) discovery, expression patterns and experimental validation in developing tea leaves indicates its correlation with circRNA-parental genes and potential roles in ceRNA interaction network. Circular RNAs (circRNAs) have recently emerged as a novel class of abundant endogenous stable RNAs produced by circularization with regulatory potential. However, identification of circRNAs in plants, especially in non-model plants with large genomes, is challenging. In this study, we undertook a systematic identification of circRNAs from different stage tissues of tea plant (Camellia sinensis) leaf development using rRNA-depleted circular RNA-seq. By combining two state-of-the-art detecting tools, we characterized 3174 circRNAs, of which 342 were shared by each approach, and thus considered high-confidence circRNAs. A few predicted circRNAs were randomly chosen, and 20 out of 24 were experimental confirmed by PCR and Sanger sequencing. Similar in other plants, tissue-specific expression was also observed for many C. sinensis circRNAs. In addition, we found that circRNA abundances were positively correlated with the mRNA transcript abundances of their parental genes. qRT-PCR validated the differential expression patterns of circRNAs between leaf bud and young leaf, which also indicated the low expression abundance of circRNAs compared to the standard mRNAs from the parental genes. We predicted the circRNA-microRNA interaction networks, and 54 of the differentially expressed circRNAs were found to have potential tea plant miRNA binding sites. The gene sets encoding circRNAs were significantly enriched in chloroplasts related GO terms and photosynthesis/metabolites biosynthesis related KEGG pathways, suggesting the candidate roles of circRNAs in photosynthetic machinery and metabolites biosynthesis during leaf development.

Revealing of MicroRNA Involved Regulatory Gene Networks on Terpenoid Biosynthesis in Camellia sinensis in Different Growing Time Points

Tea, made from leaves of Camellia sinensis, has long been consumed worldwide for its unique taste and aroma. Terpenoids play important roles not only in tea beverage aroma formation, but also in the productivity and quality of tea plantation due to their significant contribution to light harvesting pigments and phytohormones. To date, however, the regulation of terpenoid synthase genes remains unclear. Herein, the analyses of metabolomics, sRNAs, degradome, and transcriptomics were performed and integrated for identifying key regulatory miRNA-target circuits on terpenoid biosynthesis in leaf tissues over five different months in which the amount of terpenoids in tea leaves varies greatly. Four classes of miRNA-TF pairs that might play a central role in the regulation of terpenoid biosynthesis were also uncovered. Ultimately, a hypothetical model was proposed that mature miRNAs maintained by light regulator at both the transcriptional and posttranscriptional levels negatively regulate the targets to control terpenoid biosynthesis.

Complementary iTRAQ Proteomic and Transcriptomic Analyses of Leaves in Tea Plant ( Camellia sinensis L.) with Different Maturity and Regulatory Network of Flavonoid Biosynthesis

The quality of tea is highly related with the maturity of the fresh tea leaves at harvest. The present study investigated the proteomic and transcriptomic profiles of tea leaves with different maturity, using iTRAQ and RNA-seq technologies. A total of 4455 proteins and 27 930 unigenes were identified, with functional enrichment analyses of GO categorization and KEGG annotation. The compositions of flavonoids (catechins and flavonols) in tea leaves were determined. The total content of flavonoids decreased with leaf maturity, in accordance with the protein regulation patterns of shikimate, phenylpropanoid, and flavonoid pathways. The abundance of ANR had a positive correlation with epi-catechin content, while LAR abundance was positively related with catechin content ( P < 0.05). The biosynthetic network of flavonoid biosynthesis was discussed in combination with photosynthesis, primary metabolism, and transcription factors. Bud had the lowest activities of photosynthesis and carbon fixation but the highest flavonoid biosynthesis ability in opposite to mature leaf. SUS-INV switch might be an important joint for carbon flow shifting into the follow-up biochemical syntheses. This work provided a comprehensive overview on the functional protein profile changes of tea leaves at different growing stages and also proposed a research direction regarding the correlations between primary metabolism and flavonoid biosynthesis.