Identification and characterization of cold-responsive microRNAs in tea plant (Camellia sinensis) and their targets using high-throughput sequencing and degradome analysis

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.

Non-coding RNAs (ncRNAs) in plant: Master regulators for adapting to extreme temperature conditions

Unusual daily temperature fluctuations caused by climate change and climate variability adversely impact agricultural crop production. Since plants are immobile and constantly receive external environmental signals, such as extreme high (heat) and low (cold) temperatures, they have developed complex molecular regulatory mechanisms to cope with stressful situations to sustain their natural growth and development. Among these mechanisms, non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), small-interfering RNAs (siRNAs), and long-non-coding RNAs (lncRNAs), play a significant role in enhancing heat and cold stress tolerance. This review explores the pivotal findings related to miRNAs, siRNAs, and lncRNAs, elucidating how they functionally regulate plant adaptation to extreme temperatures. In addition, this review addresses the challenges associated with uncovering these non-coding RNAs and understanding their roles in orchestrating heat and cold tolerance in plants.

Phased secondary small interfering RNAs in Camellia sinensis var. assamica

Phased secondary small interfering RNAs (phasiRNAs) in plants play important roles in regulating genome stability, plant development and stress adaption. Camellia sinensis var. assamica has immense economic, medicinal and cultural significance. However, there are still no studies of phasiRNAs and their putative functions in this valuable plant. We identified 476 and 43 PHAS loci which generated 4290 twenty one nucleotide (nt) and 264 twenty four nt phasiRNAs, respectively. Moreover, the analysis of degradome revealed more than 35000 potential targets for these phasiRNAs. We identified several conserved 21 nt phasiRNA generation pathways in tea plant, including miR390 → TAS3, miR482/miR2118 → NB-LRR, miR393 → F-box, miR828 → MYB/TAS4, and miR7122 → PPR in this study. Furthermore, we found that some transposase and plant mobile domain genes could generate phasiRNAs. Our results show that phasiRNAs target genes in the same family in cis- or trans-manners, and different members of the same gene family may generate the same phasiRNAs. The phasiRNAs, generated by transposase and plant mobile domain genes, and their targets, suggest that phasiRNAs may be involved in the inhibition of transposable elements in tea plant. To summarize, these results provide a comprehensive view of phasiRNAs in Camellia sinensis var. assamica.

Functional role of microRNA in the regulation of biotic and abiotic stress in agronomic plants

The increasing demand for food is the result of an increasing population. It is crucial to enhance crop yield for sustainable production. Recently, microRNAs (miRNAs) have gained importance because of their involvement in crop productivity by regulating gene transcription in numerous biological processes, such as growth, development and abiotic and biotic stresses. miRNAs are small, non-coding RNA involved in numerous other biological functions in a plant that range from genomic integrity, metabolism, growth, and development to environmental stress response, which collectively influence the agronomic traits of the crop species. Additionally, miRNA families associated with various agronomic properties are conserved across diverse plant species. The miRNA adaptive responses enhance the plants to survive environmental stresses, such as drought, salinity, cold, and heat conditions, as well as biotic stresses, such as pathogens and insect pests. Thus, understanding the detailed mechanism of the potential response of miRNAs during stress response is necessary to promote the agronomic traits of crops. In this review, we updated the details of the functional aspects of miRNAs as potential regulators of various stress-related responses in agronomic plants.

Phyto-microbiome to mitigate abiotic stress in crop plants

Plant-associated microbes include taxonomically diverse communities of bacteria, archaebacteria, fungi, and viruses, which establish integral ecological relationships with the host plant and constitute the phyto-microbiome. The phyto-microbiome not only contributes in normal growth and development of plants but also plays a vital role in the maintenance of plant homeostasis during abiotic stress conditions. Owing to its immense metabolic potential, the phyto-microbiome provides the host plant with the capability to mitigate the abiotic stress through various mechanisms like production of antioxidants, plant growth hormones, bioactive compounds, detoxification of harmful chemicals and toxins, sequestration of reactive oxygen species and other free radicals. A deeper understanding of the structure and functions of the phyto-microbiome and the complex mechanisms of phyto-microbiome mediated abiotic stress mitigation would enable its utilization for abiotic stress alleviation of crop plants and development of stress-resistant crops. This review aims at exploring the potential of phyto-microbiome to alleviate drought, heat, salinity and heavy metal stress in crop plants and finding sustainable solutions to enhance the agricultural productivity. The mechanistic insights into the role of phytomicrobiome in imparting abiotic stress tolerance to plants have been summarized, that would be helpful in the development of novel bioinoculants. The high-throughput modern approaches involving candidate gene identification and target gene modification such as genomics, metagenomics, transcriptomics, metabolomics, and phyto-microbiome based genetic engineering have been discussed in wake of the ever-increasing demand of climate resilient crop plants.

Genome-wide identification of nitrate-responsive microRNAs by small RNA sequencing in the rice restorer cultivar Nanhui 511

Rice productivity relies heavily on nitrogen fertilization, and improving nitrogen use efficiency (NUE) is important for hybrid rice breeding. Reducing nitrogen inputs is the key to achieving sustainable rice production and reducing environmental problems. Here, we analyzed the genome-wide transcriptomic changes in microRNAs (miRNAs) in the indica rice restorer cultivar Nanhui 511 (NH511) under high (HN) and low nitrogen (LN) conditions. The results showed that NH511 is sensitive to nitrogen supplies and HN conditions promoted the growth its lateral roots at the seedling stage. Furthermore, we identified 483 known miRNAs and 128 novel miRNAs by small RNA sequencing in response to nitrogen in NH511. We also detected 100 differentially expressed genes (DEGs), including 75 upregulated and 25 downregulated DEGs, under HN conditions. Among these DEGs, 43 miRNAs that exhibited a 2-fold change in their expression were identified in response to HN conditions, including 28 upregulated and 15 downregulated genes. Additionally, some differentially expressed miRNAs were further validated by qPCR analysis, which showed that miR443, miR1861b, and miR166k-3p were upregulated, whereas miR395v and miR444b.1 were downregulated under HN conditions. Moreover, the degradomes of possible target genes for miR166k-3p and miR444b.1 and expression variations were analyzed by qPCR at different time points under HN conditions. Our findings revealed comprehensive expression profiles of miRNAs responsive to HN treatments in an indica rice restorer cultivar, which advances our understanding of the regulation of nitrogen signaling mediated by miRNAs and provides novel data for high-NUE hybrid rice cultivation.

Whole-Transcriptome Sequencing Reveals a ceRNA Regulatory Network Associated with the Process of Periodic Albinism under Low Temperature in Baiye No. 1 (Camellia sinensis)

The young shoots of the tea plant Baiye No. 1 display an albino phenotype in the early spring under low environmental temperatures, and the leaves re-green like those of common tea cultivars during the warm season. Periodic albinism is precisely regulated by a complex gene network that leads to metabolic differences and enhances the nutritional value of tea leaves. Here, we identified messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs) to construct competing endogenous RNA (ceRNA) regulatory networks. We performed whole-transcriptome sequencing of 12 samples from four periods (Bud, leaves not expanded; Alb, albino leaves; Med, re-greening leaves; and Gre, green leaves) and identified a total of 6325 differentially expressed mRNAs (DEmRNAs), 667 differentially expressed miRNAs (DEmiRNAs), 1702 differentially expressed lncRNAs (DElncRNAs), and 122 differentially expressed circRNAs (DEcircRNAs). Furthermore, we constructed ceRNA networks on the basis of co-differential expression analyses which comprised 112, 35, 38, and 15 DEmRNAs, DEmiRNAs, DElncRNAs, and DEcircRNAs, respectively. Based on the regulatory networks, we identified important genes and their interactions with lncRNAs, circRNAs, and miRNAs during periodic albinism, including the ceRNA regulatory network centered on miR5021x, the GAMYB-miR159-lncRNA regulatory network, and the NAC035-miR319x-circRNA regulatory network. These regulatory networks might be involved in the response to cold stress, photosynthesis, chlorophyll synthesis, amino acid synthesis, and flavonoid accumulation. Our findings provide novel insights into ceRNA regulatory mechanisms involved in Baiye No. 1 during periodic albinism and will aid future studies of the molecular mechanisms underlying albinism mutants.

Recent progress and perspectives on physiological and molecular mechanisms underlying cold tolerance of tea plants

Tea is one of the most consumed and widely planted beverage plant worldwide, which contains many important economic, healthy, and cultural values. Low temperature inflicts serious damage to tea yields and quality. To cope with cold stress, tea plants have evolved a cascade of physiological and molecular mechanisms to rescue the metabolic disorders in plant cells caused by the cold stress; this includes physiological, biochemical changes and molecular regulation of genes and associated pathways. Understanding the physiological and molecular mechanisms underlying how tea plants perceive and respond to cold stress is of great significance to breed new varieties with improved quality and stress resistance. In this review, we summarized the putative cold signal sensors and molecular regulation of the CBF cascade pathway in cold acclimation. We also broadly reviewed the functions and potential regulation networks of 128 cold-responsive gene families of tea plants reported in the literature, including those particularly regulated by light, phytohormone, and glycometabolism. We discussed exogenous treatments, including ABA, MeJA, melatonin, GABA, spermidine and airborne nerolidol that have been reported as effective ways to improve cold resistance in tea plants. We also present perspectives and possible challenges for functional genomic studies on cold tolerance of tea plants in the future.

Effect of shading on physiological attributes and comparative transcriptome analysis of Camellia sinensis cultivar reveals tolerance mechanisms to low temperatures

Tea is a vital beverage crop all over the world, including in China. Low temperatures restrict its growth, development, and terrestrial distribution, and cold event variability worsens cold damage. However, the physiological and molecular mechanisms of Camellia sinensis under shade in winter remain unclear. In our study, tea leaves were utilized for physiological attributes and transcriptome analysis in November and December in three shading groups and no-shade control plants. When compared to the no-shade control plants, the shading group protected tea leaves from cold damage, increased photochemical efficiency (Fv/Fm) and soil plant analysis development (SPAD), and sustained chlorophyll a, chlorophyll b, chlorophyll, and carotenoid contents by physiological mean. Then, transcriptome analysis revealed 20,807 differentially expressed genes (DEGs) and transcription factors (TFs) in November and December. A comparative study of transcriptome resulted in 3,523 DEGs and many TFs under SD0% vs. SD30%, SD0% vs. SD60%, and SD0% vs. SD75% of shading in November and December. Statistically, 114 DEGs were downregulated and 72 were upregulated under SD0% vs. SD30%. SD0% vs. SD60% resulted in 154 DEGs, with 60 downregulated and 94 upregulated. Similarly, there were 505 DEGs of which 244 were downregulated and 263 were upregulated under SD0% vs. SD75% of shading throughout November. However, 279 DEGs were downregulated and 105 were upregulated under SD0% vs. SD30%. SD0% vs. SD60% resulted in 296 DEGs, with 172 downregulated and 124 upregulated. Finally, 2,173 DEGs were regulated in December, with 1,428 downregulated and 745 upregulated under SD0% vs. SD75%. These indicate that the number of downregulated DEGs in December was higher than the number of upregulated DEGs in November during low temperatures. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of differentially expressed genes were highly regulated in the photosynthesis, plant hormone signal transduction, and mitogen-activated protein kinase (MAPK) signaling pathways. However, qRT-PCR and RNA-seq relative expression of photosynthetic (DEGs) Lhcb2 in both November and December, plant hormone (DEGs) BRI1 and JAZ in November and IAA and ERF1 in December, and key DEGs of MAPK signal transduction FLS2, CHIB, and MPK4 in November and RBOH, MKK4_5, and MEKK1 in December in three shading groups and no-shade control plants responded to tea cold tolerance. The enhanced expression of light-harvesting photosystem I gene Lhca5, light-harvesting photosystem II gene Lhcb2, and mitogen-activated protein kinases MEKK1 and MPK4/6 enhance the cold-tolerance mechanism of C. sinensis. These comprehensive transcriptomic findings are significant for furthering our understanding of the genes and underlying regulatory mechanisms of shade-mediated low-temperature stress tolerance in horticultural crops.

Evolutionary Landscape of Tea Circular RNAs and Its Contribution to Chilling Tolerance of Tea Plant

Chilling stress threatens the yield and distribution pattern of global crops, including the tea plant (Camellia sinensis), one of the most important cash crops around the world. Circular RNA (circRNA) plays roles in regulating plant growth and biotic/abiotic stress responses. Understanding the evolutionary characteristics of circRNA and its feedbacks to chilling stress in the tea plant will help to elucidate the vital roles of circRNAs. In the current report, we systematically identified 2702 high-confidence circRNAs under chilling stress in the tea plant, and interestingly found that the generation of tea plant circRNAs was associated with the length of their flanking introns. Repetitive sequences annotation and DNA methylation analysis revealed that the longer flanking introns of circRNAs present more repetitive sequences and higher methylation levels, which suggested that repeat-elements-mediated DNA methylation might promote the circRNAs biogenesis in the tea plant. We further detected 250 differentially expressed circRNAs under chilling stress, which were functionally enriched in GO terms related to cold/stress responses. Constructing a circRNA-miRNA-mRNA interaction network discovered 139 differentially expressed circRNAs harboring potential miRNA binding sites, which further identified 14 circRNAs that might contribute to tea plant chilling responses. We further characterized a key circRNA, CSS-circFAB1, which was significantly induced under chilling stress. FISH and silencing experiments revealed that CSS-circFAB1 was potentially involved in chilling tolerance of the tea plant. Our study emphasizes the importance of circRNA and its preliminary role against low-temperature stress, providing new insights for tea plant cold tolerance breeding.

MicroRNAs Mediated Plant Responses to Salt Stress

One of the most damaging issues to cultivatable land is soil salinity. While salt stress influences plant growth and yields at low to moderate levels, severe salt stress is harmful to plant growth. Mineral shortages and toxicities frequently exacerbate the problem of salinity. The growth of many plants is quantitatively reduced by various levels of salt stress depending on the stage of development and duration of stress. Plants have developed various mechanisms to withstand salt stress. One of the key strategies is the utilization of microRNAs (miRNAs) that can influence gene regulation at the post-transcriptional stage under different environmental conditions, including salinity. Here, we have reviewed the miRNA-mediated adaptations of various plant species to salt stress and other abiotic variables. Moreover, salt responsive (SR)-miRNAs, their targets, and corresponding pathways have also been discussed. The review article concludes by suggesting that the utilization of miRNAs may be a vital strategy to generate salt tolerant crops ensuring food security in the future.

Application of RNAi technology: a novel approach to navigate abiotic stresses

BACKGROUND

Due to the rising population globally, and the demand for food, it is critical to significantly increase crop production by 2050. However, climate change estimates show that droughts and heatwaves will become more prevalent in many parts of the world, posing a severe danger to food output.

METHODS

Selective breeding based on genetic diversity is falling short of meeting the expanding need for food and feed. However, the advent of modern plant genetic engineering, genome editing, and synthetic biology provides precise techniques for producing crops capable of sustaining yield under stress situations.

RESULTS

As a result, crop varieties with built-in genetic tolerance to environmental challenges are desperately needed. In the recent years, small RNA (sRNA) data has progressed to become one of the most effective approaches for the improvement of crops. So many sRNAs (18-30nt) have been found with the use of hi-tech bioinformatics and sequencing techniques which are involved in the regulation of sequence specific gene noncoding RNAs (short ncRNAs) i.e., microRNA (miRNA) and small interfering RNA (siRNA). Such research outcomes may advance our understanding of the genetic basis of adaptability of plants to various environmental challenges and the genetic variation of plant's tolerance to a number of abiotic stresses.

CONCLUSION

The review article highlights current trends and advances in sRNAs' critical role in responses of plants to drought, heat, cold, and salinity, and also the potential technology that identifies the abiotic stress-regulated sRNAs, and techniques for analyzing and validating the target genes.

Development and verification of SSR markers from drought stress-responsive miRNAs in Dongxiang wild rice (Oryza rufipogon Griff.)

Rice production worldwide has continued to decline due to various environmental stresses, with drought stress being a prominent factor, as rice is a semi-aquatic plant. Thus, development of drought stress-resistant rice varieties is of great importance for rice production. In our previous study, we found that microRNAs (miRNAs) play a crucial role in the response to drought stress in Dongxiang wild rice (DXWR) (Oryza rufipogon Griff.). Developing drought stress-responsive miRNA-based single sequence repeat (SSR) markers for DXWR will be of great value for the efficient identification and utilization of miRNA genes to breed drought stress-resistant rice varieties. In this study, ninety-nine novel SSR markers were developed based on the drought stress-responsive miRNAs of DXWR. These markers were distributed in all 12 rice chromosomes, and most were in chromosomes 2 and 6, with di- and tri-nucleotides being the most abundant repeat motifs. Twelve out of fourteen synthesized markers displayed high levels of genetic diversity in the genomes of three populations of DXWR and 40 modern rice varieties worldwide. The number of alleles per locus ranged from 2 to 7, with an average of 4.67; the genetic diversity index ranged from 0.21 to 0.76, with an average of 0.58; and the polymorphism information content value ranged from 0.18 to 0.72, with an average of 0.53. These novel molecular markers developed from the drought stress-responsive miRNAs of DXWR could be additional tools for mapping elite miRNA genes and breeding drought stress-resistant rice varieties.

HSP70 Gene Family in Brassica rapa: Genome-Wide Identification, Characterization, and Expression Patterns in Response to Heat and Cold Stress

Heat shock proteins protect plants from abiotic stress, such as salt, drought, heat, and cold stress. HSP70 is one of the major members of the heat shock protein family. To explore the mechanism of HSP70 in Brassica rapa, we identified 28 putative HSP70 gene family members using state-of-the-art bioinformatics-based tools and methods. Based on chromosomal mapping, HSP70 genes were the most differentially distributed on chromosome A03 and the least distributed on chromosome A05. Ka/Ks analysis revealed that B. rapa evolution was subjected to intense purifying selection of the HSP70 gene family. RNA-sequencing data and expression profiling showed that heat and cold stress induced HSP70 genes. The qRT-PCR results verified that the HSP70 genes in Chinese cabbage (Brassica rapa ssp. pekinensis) are stress-inducible under both cold and heat stress. The upregulated expression pattern of these genes indicated the potential of HSP70 to mitigate environmental stress. These findings further explain the molecular mechanism underlying the responses of HSP70 to heat and cold stress.

MTAGCN: predicting miRNA-target associations in Camellia sinensis var. assamica through graph convolution neural network

Background

MircoRNAs (miRNAs) play a central role in diverse biological processes of Camellia sinensis var.assamica (CSA) through their associations with target mRNAs, including CSA growth, development and stress response. However, although the experiment methods of CSA miRNA-target identifications are costly and time-consuming, few computational methods have been developed to tackle the CSA miRNA-target association prediction problem.

Results

In this paper, we constructed a heterogeneous network for CSA miRNA and targets by integrating rich biological information, including a miRNA similarity network, a target similarity network, and a miRNA-target association network. We then proposed a deep learning framework of graph convolution networks with layer attention mechanism, named MTAGCN. In particular, MTAGCN uses the attention mechanism to combine embeddings of multiple graph convolution layers, employing the integrated embedding to score the unobserved CSA miRNA-target associations.

Discussion

Comprehensive experiment results on two tasks (balanced task and unbalanced task) demonstrated that our proposed model achieved better performance than the classic machine learning and existing graph convolution network-based methods. The analysis of these results could offer valuable information for understanding complex CSA miRNA-target association mechanisms and would make a contribution to precision plant breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04819-3.

miR319a targeting of CsTCP10 plays an important role in defense against gray blight disease in tea plant (Camellia sinensis)

Gray blight disease occurs widely in major tea-producing areas and harms the leaves of tea trees, which affects the quality and yield of processed tea. According to an analysis of previous sequencing data, miR319a may be important in the resistance of tea plants to gray blight disease. In this study, based on 5'RLM-RACE, qRT-PCR, sODN, CIN and transient transformation experiments in tobacco, CsTCP10 and CsTCP4 were found to be cleaved by miR319a. qRT-PCR and northern blotting also revealed that the expression pattern of CsTCP10 in tea leaves was opposite to that of miR319a, while that of CsTCP4 displayed no similar change. Furthermore, a large amount of reactive oxygen species was found to accumulate in tea leaves in the antisense oligodeoxynucleotide experiment, while the expression of CsTCP10 was inhibited. These results suggest that CsTCP10 is a positive regulator of the resistance of tea plants to gray blight disease. Compared with the wild-type, the expression of AtTCP10 in transgenic Arabidopsis plants was downregulated. After infection with the pathogen, the transgenic plants were more severely damaged. Our results suggest that miR319a facilitates Pestalotiopsis infection by suppressing the expression of CsTCP10 in tea plants.

Comprehensive Genome-Wide Identification and Expression Profiling of Eceriferum (CER) Gene Family in Passion Fruit (Passiflora edulis) Under Fusarium kyushuense and Drought Stress Conditions

Plant surfaces are covered with cuticle wax and are the first barrier between a plant and environmental stresses. Eceriferum (CER) is an important gene family involved in wax biosynthesis and stress resistance. In this study, for the first time, 34 CER genes were identified in the passion fruit (Passiflora edulis) genome, and PeCER proteins varied in physicochemical properties. A phylogenetic tree was constructed and divided into seven clades to identify the evolutionary relationship with other plant species. Gene structure analyses revealed that conserved motifs ranged from 1 to 24, and that exons ranged from 1 to 29. The cis-element analysis provides insight into possible roles of PeCER genes in plant growth, development and stress responses. The syntenic analysis revealed that segmental (six gene pairs) and tandem (six gene pairs) gene duplication played an important role in the expansion of PeCER genes and underwent a strong purifying selection. In addition, 12 putative ped-miRNAs were identified to be targeting 16 PeCER genes, and PeCER6 was the most targeted by four miRNAs including ped-miR157a-5p, ped-miR164b-5p, ped-miR319b, and ped-miR319l. Potential transcription factors (TFs) such as ERF, AP2, MYB, and bZIP were predicted and visualized in a TF regulatory network interacting with PeCER genes. GO and KEGG annotation analysis revealed that PeCER genes were highly related to fatty acid, cutin, and wax biosynthesis, plant-pathogen interactions, and stress response pathways. The hypothesis that most PeCER proteins were predicted to localize to the plasma membrane was validated by transient expression assays of PeCER32 protein in onion epidermal cells. qRT-PCR expression results showed that most of the PeCER genes including PeCER1, PeCER11, PeCER15, PeCER17, and PeCER32 were upregulated under drought and Fusarium kyushuense stress conditions compared to controls. These findings provide a foundation for further studies on functions of PeCER genes to further facilitate the genetic modification of passion fruit wax biosynthesis and stress resistance.

Small RNA and Degradome Deep Sequencing Reveals the Roles of microRNAs in Peanut (Arachis hypogaea L.) Cold Response

Cold stress is a major environmental factor that affects plant growth and development, as well as fruit postharvest life and quality. MicroRNAs (miRNAs) are a class of non-coding small RNAs that play crucial roles in various abiotic stresses. Peanuts (Arachis hypogaea L.), one of the most important grain legumes and source of edible oils and proteins, are cultivated in the semi-arid tropical and subtropical regions of the world. To date, there has been no report on the role of miRNAs in the response to cold stress in cultivated peanuts. In this study, we profiled cold-responsive miRNAs in peanuts using deep sequencing in cold-sensitive (WQL20) alongside a tolerant variety (WQL30). A total of 407 known miRNAs and 143 novel peanut-specific miRNAs were identified. The expression of selected known and novel miRNAs was validated by northern blotting and six known cold-responsive miRNAs were revealed. Degradome sequencing identified six cold-responsive miRNAs that regulate 12 target genes. The correlative expression patterns of several miRNAs and their target genes were further validated using qRT-PCR. Our data showed that miR160-ARF, miR482-WDRL, miR2118-DR, miR396-GRF, miR162-DCL, miR1511-SRF, and miR1511-SPIRAL1 modules may mediate cold stress responses. Transient expression analysis in Nicotiana benthamiana found that miR160, miR482, and miR2118 may play positive roles, and miR396, miR162, and miR1511 play negative roles in the regulation of peanut cold tolerance. Our results provide a foundation for understanding miRNA-dependent cold stress response in peanuts. The characterized correlations between miRNAs and their response to cold stress could serve as markers in breeding programs or tools for improving cold tolerance of peanuts.

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.

Genome-Wide Identification and Expression Profiling of KCS Gene Family in Passion Fruit (Passiflora edulis) Under Fusarium kyushuense and Drought Stress Conditions

Plant and fruit surfaces are covered with cuticle wax and provide a protective barrier against biotic and abiotic stresses. Cuticle wax consists of very-long-chain fatty acids (VLCFAs) and their derivatives. β-Ketoacyl-CoA synthase (KCS) is a key enzyme in the synthesis of VLCFAs and provides a precursor for the synthesis of cuticle wax, but the KCS gene family was yet to be reported in the passion fruit (Passiflora edulis). In this study, thirty-two KCS genes were identified in the passion fruit genome and phylogenetically grouped as KCS1-like, FAE1-like, FDH-like, and CER6-like. Furthermore, thirty-one PeKCS genes were positioned on seven chromosomes, while one PeKCS was localized to the unassembled genomic scaffold. The cis-element analysis provides insight into the possible role of PeKCS genes in phytohormones and stress responses. Syntenic analysis revealed that gene duplication played a crucial role in the expansion of the PeKCS gene family and underwent a strong purifying selection. All PeKCS proteins shared similar 3D structures, and a protein-protein interaction network was predicted with known Arabidopsis proteins. There were twenty putative ped-miRNAs which were also predicted that belong to nine families targeting thirteen PeKCS genes. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation results were highly associated with fatty acid synthase and elongase activity, lipid metabolism, stress responses, and plant-pathogen interaction. The highly enriched transcription factors (TFs) including ERF, MYB, Dof, C2H2, TCP, LBD, NAC, and bHLH were predicted in PeKCS genes. qRT-PCR expression analysis revealed that most PeKCS genes were highly upregulated in leaves including PeKCS2, PeKCS4, PeKCS8, PeKCS13, and PeKCS9 but not in stem and roots tissues under drought stress conditions compared with controls. Notably, most PeKCS genes were upregulated at 9th dpi under Fusarium kyushuense biotic stress condition compared to controls. This study provides a basis for further understanding the functions of KCS genes, improving wax and VLCFA biosynthesis, and improvement of passion fruit resistance.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Comparative analysis identifies micro-RNA associated with nutrient homeostasis, development and stress response in Arabidopsis thaliana upon high Zn and metal hyperaccumulator Arabidopsis halleri

miRNAs have been found to be key players in mineral homeostasis, both in the control of nutrient balance and in the response to toxic trace elements. However, the effect of Zn excess on miRNAs has not been elucidated; moreover, no data are present regarding miRNAs in hyperaccumulator species, where metal homeostasis is tightly regulated. Therefore, expression levels of mature miRNAs were measured by RNA-Seq in Zn-sensitive Arabidopsis thaliana grown in control conditions and upon high Zn, in soil and in Zn-hyperaccumulator Arabidopsis halleri grown in control conditions. Differential expression of notable miRNAs and their targets was confirmed by real-time RT-PCR. The comparison in A. thaliana revealed a small subset modulated upon Zn treatment that is associated with stress response and nutrient homeostasis. On the other hand, a more consistent group of miRNAs was differentially expressed in A. halleri compared with A. thaliana, reflecting inherent differences in nutritional requirements and response to stresses and plant growth and development. Overall, these results confirm the involvement of miRNAs in Zn homeostasis and support the hypothesis of distinct regulatory pathways in hyperaccumulator species.

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.

MicroRNA Omics Analysis of Camellia sinesis Pollen Tubes in Response to Low-Temperature and Nitric Oxide

Nitric oxide (NO) as a momentous signal molecule participates in plant reproductive development and responds to various abiotic stresses. Here, the inhibitory effects of the NO-dominated signal network on the pollen tube growth of Camellia sinensis under low temperature (LT) were studied by microRNA (miRNA) omics analysis. The results showed that 77 and 71 differentially expressed miRNAs (DEMs) were induced by LT and NO treatment, respectively. Gene ontology (GO) analysis showed that DEM target genes related to microtubules and actin were enriched uniquely under LT treatment, while DEM target genes related to redox process were enriched uniquely under NO treatment. In addition, the target genes of miRNA co-regulated by LT and NO are only located on the cell membrane and cell wall, and most of them are enriched in metal ion binding and/or transport and cell wall organization. Furthermore, DEM and its target genes related to metal ion binding/transport, redox process, actin, cell wall organization and carbohydrate metabolism were identified and quantified by functional analysis and qRT-PCR. In conclusion, miRNA omics analysis provides a complex signal network regulated by NO-mediated miRNA, which changes cell structure and component distribution by adjusting Ca²⁺ gradient, thus affecting the polar growth of the C. sinensis pollen tube tip under LT.

Divergent DNA methylation contributes to duplicated gene evolution and chilling response in tea plants

The tea plant (Camellia sinensis) is a thermophilic cash crop and contains a highly duplicated and repeat-rich genome. It is still unclear how DNA methylation regulates the evolution of duplicated genes and chilling stress in tea plants. We therefore generated a single-base-resolution DNA methylation map of tea plants under chilling stress. We found that, compared with other plants, the tea plant genome is highly methylated in all three sequence contexts, including CG, CHG and CHH (where H = A, T, or C), which is further proven to be correlated with its repeat content and genome size. We show that DNA methylation in the gene body negatively regulates the gene expression of tea plants, whereas non-CG methylation in the flanking region enables a positive regulation of gene expression. We demonstrate that transposable element-mediated methylation dynamics significantly drives the expression divergence of duplicated genes in tea plants. The DNA methylation and expression divergence of duplicated genes in the tea plant increases with evolutionary age and selective pressure. Moreover, we detect thousands of differentially methylated genes, some of which are functionally associated with chilling stress. We also experimentally reveal that DNA methyltransferase genes of tea plants are significantly downregulated, whereas demethylase genes are upregulated at the initial stage of chilling stress, which is in line with the significant loss of DNA methylation of three well-known cold-responsive genes at their promoter and gene body regions. Overall, our findings underscore the importance of DNA methylation regulation and offer new insights into duplicated gene evolution and chilling tolerance in tea plants.

Multi-omics analysis to visualize the dynamic roles of defense genes in the response of tea plants to gray blight

Gray blight (GB) is one of the most destructive diseases of tea plants, causing considerable damage and productivity losses; however, the dynamic roles of defense genes during pathogen infection remain largely unclear. To explore the numerous molecular interactions associated with GB stress in tea plants, we employed transcriptome, sRNAome and degradome sequencing from 1 to 13 days post-inoculation (dpi) at 3-day intervals. The transcriptomics results showed that differentially expressed genes (DEGs) related to flavonoid synthesis, such as chalcone synthase (CHS) and phenylalanine ammonia-lyase (PAL), were particularly induced at 4 dpi. Consistent with this, the contents of catechins (especially gallocatechin), which are the dominant flavonoids in tea plants, also increased in the leaves of tea plants infected with GB. Combined analysis of the sRNAome and degradome revealed that microRNAs could mediate tea plant immunity by regulating DEG expression at the post-transcriptional level. Co-expression network analysis demonstrated that miR530b-ethylene responsive factor 96 (ERF96) and miRn211-thaumatin-like protein (TLP) play crucial roles in the response to GB. Accordingly, gene-specific antisense oligonucleotide assays suggested that suppressing ERF96 decreased the levels of reactive oxygen species (ROS), whereas suppressing TLP increased the levels of ROS. Furthermore, ERF96 was induced, but TLP was suppressed, in susceptible tea cultivars. Our results collectively demonstrate that ERF96 is a negative regulator and TLP is a positive regulator in the response of tea plants to GB. Taken together, our comprehensive integrated analysis reveals a dynamic regulatory network linked to GB stress in tea plants and provides candidate genes for improvement of tea plants.

CsTCPs regulate shoot tip development and catechin biosynthesis in tea plant (Camellia sinensis)

The growth of leaves and biosynthesis of characteristic secondary metabolites are critically important for tea production and quality control. However, little is known about the coordinated regulation of leaf development and catechin biosynthesis in tea plants. Here, we reported that TCP TFs are involved in both catechin biosynthesis and leaf development. An integrated analysis of catechin profiling and CsTCP expression in different tissues of plants under various environmental conditions at different developmental stages indicated significant correlations between the transcript levels of CIN-type TCPs and catechin production. CIN-type CsTCP3 and CsTCP4 and PCF-type CsTCP14 interacted with the MYB-bHLH-WD40 repeat (MBW) complex by forming a CsTCP3-CsTT8 heterodimer and modulating the transactivation activity of the promoters of anthocyanin synthase (CsANS1) and anthocyanidin reductase (CsANR1). Four types of microRNA/target modules, miR319b/CsTCP3-4, miR164b/CsCUC, miR396/CsGRF-GIF, and miR165b/HD-ZIPIII ones, were also identified and characterized for their functions in the regulation of the development of tea plant shoot tips and leaf shape. The results of these modules were reflected by their different expression patterns in developing buds and leaves that had distinctly different morphologies in three different tea plant varieties. Their roles in the regulation of catechin biosynthesis were also further verified by manipulation of microRNA319b (miR319b), which targets the transcripts of CsTCP3 and CsTCP4. Thus, CsTCPs represent at least one of these important groups of TFs that can integrate tea plant leaf development together with secondary metabolite biosynthesis. Our study provides new insight into shoot tip development and catechin production in tea plants and lays a foundation for further mechanistic understanding of the regulation of tea plant leaf development and secondary metabolism.

Genome-wide identification, evolutionary relationship and expression analysis of AGO, DCL and RDR family genes in tea

Three gene families in plants viz. Argonaute (AGOs), Dicer-like (DCLs) and RNA dependent RNA polymerase (RDRs) constitute the core components of small RNA mediated gene silencing machinery. The present study endeavours to identify members of these gene families in tea and to investigate their expression patterns in different tissues and various stress regimes. Using genome-wide analysis, we have identified 18 AGOs, 5 DCLs and 9 RDRs in tea, and analyzed their phylogenetic relationship with orthologs of Arabidopsis thaliana. Gene expression analysis revealed constitutive expression of CsAGO1 in all the studied tissues and stress conditions, whereas CsAGO10c showed most variable expression among all the genes. CsAGO10c gene was found to be upregulated in tissues undergoing high meristematic activity such as buds and roots, as well as in Exobasidium vexans infected samples. CsRDR2 and two paralogs of CsAGO4, which are known to participate in biogenesis of hc-siRNAs, showed similarities in their expression levels in most of the tea plant tissues. This report provides first ever insight into the important gene families involved in biogenesis of small RNAs in tea. The comprehensive knowledge of these small RNA biogenesis purveyors can be utilized for tea crop improvement aimed at stress tolerance and quality enhancement.

Catalase (CAT) Gene Family in Rapeseed (Brassica napus L.): Genome-Wide Analysis, Identification, and Expression Pattern in Response to Multiple Hormones and Abiotic Stress Conditions

Catalase (CAT) is an antioxidant enzyme expressed by the CAT gene family and exists in almost all aerobic organisms. Environmental stresses induce the generation of reactive oxygen species (ROS) that eventually hinder plant growth and development. The CAT enzyme translates the hydrogen peroxide (H₂O₂) to water (H₂O) and reduce the ROS levels to shelter the cells’ death. So far, the CAT gene family has not been reported in rapeseed (Brassica napus L.). Therefore, a genome-wide comprehensive analysis was conducted to classify the CAT genes in the rapeseed genome. The current study identified 14 BnCAT genes in the rapeseed genome. Based on phylogenetic and synteny analysis, the BnCATs belong to four groups (Groups I–IV). A gene structure and conserved motif analysis showed that Group I, Group II, and Group IV possess almost the same intron/exon pattern, and an equal number of motifs, while Group III contains diverse structures and contain 15 motifs. By analyzing the cis-elements in the promoters, we identified five hormone-correlated responsive elements and four stress-related responsive elements. Further, six putative bna-miRNAs were also identified, targeting three genes (BnCAT4, BnCAT6, and BnCAT8). Gene ontology (GO) enrichment analysis showed that the BnCAT genes were largely related to cellular organelles, ROS response, stimulus response, stress response, and antioxidant enzymes. Almost 10 BnCAT genes showed higher expression levels in different tissues, i.e., root, leaf, stem, and silique. The expression analysis showed that BnCAT1–BnCAT3 and BnCAT11–BnCAT13 were significantly upregulated by cold, salinity, abscisic acid (ABA), and gibberellic acid (GA) treatment, but not by drought and methyl jasmonate (MeJA). Notably, most of the genes were upregulated by waterlogging stress, except BnCAT6, BnCAT9, and BnCAT10. Our results opened new windows for future investigations and provided insights into the CAT family genes in rapeseed.

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.

A comprehensive analysis of cotton VQ gene superfamily reveals their potential and extensive roles in regulating cotton abiotic stress

Background

Valine-glutamine (VQ) motif-containing proteins play important roles in plant growth, development and abiotic stress response. For many plant species, the VQ genes have been identified and their functions have been described. However, little is known about the origin, evolution, and functions (and underlying mechanisms) of the VQ family genes in cotton.

Results

In this study, we comprehensively analyzed the characteristics of 268 VQ genes from four Gossypium genomes and found that the VQ proteins evolved into 10 clades, and each clade had a similar structural and conservative motif. The expansion of the VQ gene was mainly through segmental duplication, followed by dispersal. Expression analysis revealed that many GhVQs might play important roles in response to salt and drought stress, and GhVQ18 and GhVQ84 were highly expressed under PEG and salt stress. Further analysis showed that GhVQs were co-expressed with GhWRKY transcription factors (TFs), and microRNAs (miRNAs) could hybridize to their cis-regulatory elements.

Conclusions

The results in this study broaden our understanding of the VQ gene family in plants, and the analysis of the structure, conserved elements, and expression patterns of the VQs provide a solid foundation for exploring their specific functions in cotton responding to abiotic stresses. Our study provides significant insight into the potential functions of VQ genes in cotton.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s12864-020-07171-z.

Integrated Transcriptome, microRNA, and Phytochemical Analyses Reveal Roles of Phytohormone Signal Transduction and ABC Transporters in Flavor Formation of Oolong Tea (Camellia sinensis) during Solar Withering

The unique aroma and flavor of oolong tea develop during the withering stage of postharvest processing. We explored the roles of miRNA-related regulatory networks during tea withering and their effects on oolong tea quality. We conducted transcriptome and miRNA analyses to identify differentially expressed (DE) miRNAs and target genes among fresh leaves, indoor-withered leaves, and solar-withered leaves. We identified 32 DE-miRNAs and 41 target genes involved in phytohormone signal transduction and ABC transporters. Further analyses indicated that these two pathways regulated the accumulation of flavor-related metabolites during tea withering. Flavonoid accumulation was correlated with the miR167d_1-ARF-GH3, miR845-ABCC1-3/ABCC2, miR166d-5p_1-ABCC1-2, and miR319c_3-PIF-ARF modules. Terpenoid content was correlated with the miR171b-3p_2-DELLA-MYC2 and miR166d-5p_1-ABCG2-MYC2 modules. These modules inhibited flavonoid biosynthesis and enhanced terpenoid biosynthesis in solar-withered leaves. Low auxin and gibberellic acid contents and circRNA-related regulatory networks also regulated the accumulation of flavor compounds in solar-withered leaves. Our analyses reveal how solar withering produces high-quality oolong tea.

Strontium stress disrupts miRNA biogenesis by reducing HYL1 protein levels in Arabidopsis

Strontium (Sr) is an emerging environmental pollutant that has become a major global concern after the nuclear accident at the Fukushima Daiichi Nuclear Power Plant in 2011. Although many studies have demonstrated the harmful effects of Sr on plant growth and development at the physiological level, knowledge regarding how plants sense and respond to Sr stress at the molecular level is limited. Recent studies have suggested that microRNAs (miRNAs) function as key regulators of plant growth and development as well as in the responses of plants to environmental stresses, including salinity, drought, cold, nutrient starvation, and heavy metals. In this study, we examined the global expression profile of miRNAs under Sr stress using small RNA sequencing analysis in Arabidopsis to better understand the molecular basis of plant responses to Sr stress. To identify specific Sr-responsive miRNAs, we performed comparative miRNA expression profiling analysis using control, CaCl₂-, and SrCl₂-treated seedlings. Compared to the control treatment, the expressions of most miRNAs were considerably decreased in the Sr-treated seedlings. However, under Sr stress, the expressions of primary miRNAs (pri-miRNAs) and their target genes were significantly increased; the protein levels of HYPONASTIC LEAVES 1 (HYL1), one of the core components of the microprocessor complex, were strongly reduced despite the increased HYL1 mRNA expression. In addition, hyl1-2 mutant plants were shown to be more sensitive to Sr stress than wild-type plants. Collectively, our results strongly suggested that Sr stress may be associated with the disruption of miRNA biogenesis by reducing the protein level of HYL1, which is required to maintain proper growth and development for plants. Our findings further indicated that some miRNAs may play important roles in plant responses to Sr stress.

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.

Late Embryogenesis Abundant Protein-Client Protein Interactions

The intrinsically disordered proteins belonging to the LATE EMBRYOGENESIS ABUNDANT protein (LEAP) family have been ascribed a protective function over an array of intracellular components. We focus on how LEAPs may protect a stress-susceptible proteome. These examples include instances of LEAPs providing a shield molecule function, possibly by instigating liquid-liquid phase separations. Some LEAPs bind directly to their client proteins, exerting a holdase-type chaperonin function. Finally, instances of LEAP-client protein interactions have been documented, where the LEAP modulates (interferes with) the function of the client protein, acting as a surreptitious rheostat of cellular homeostasis. From the examples identified to date, it is apparent that client protein modulation also serves to mitigate stress. While some LEAPs can physically bind and protect client proteins, some apparently bind to assist the degradation of the client proteins with which they associate. Documented instances of LEAP-client protein binding, even in the absence of stress, brings to the fore the necessity of identifying how the LEAPs are degraded post-stress to render them innocuous, a first step in understanding how the cell regulates their abundance.

Identification of small RNAs during cold acclimation in Arabidopsis thaliana

BACKGROUND

Cold stress causes dynamic changes in gene expression that are partially caused by small non-coding RNAs since they regulate protein coding transcripts and act in epigenetic gene silencing pathways. Thus, a detailed analysis of transcriptional changes of small RNAs (sRNAs) belonging to all known sRNA classes such as microRNAs (miRNA) and small interfering RNA (siRNAs) in response to cold contributes to an understanding of cold-related transcriptome changes.

RESULT

We subjected A. thaliana plants to cold acclimation conditions (4 °C) and analyzed the sRNA transcriptomes after 3 h, 6 h and 2 d. We found 93 cold responsive differentially expressed miRNAs and only 14 of these were previously shown to be cold responsive. We performed miRNA target prediction for all differentially expressed miRNAs and a GO analysis revealed the overrepresentation of miRNA-targeted transcripts that code for proteins acting in transcriptional regulation. We also identified a large number of differentially expressed cis- and trans-nat-siRNAs, as well as sRNAs that are derived from long non-coding RNAs. By combining the results of sRNA and mRNA profiling with miRNA target predictions and publicly available information on transcription factors, we reconstructed a cold-specific, miRNA and transcription factor dependent gene regulatory network. We verified the validity of links in the network by testing its ability to predict target gene expression under cold acclimation.

CONCLUSION

In A. thaliana, miRNAs and sRNAs derived from cis- and trans-NAT gene pairs and sRNAs derived from lncRNAs play an important role in regulating gene expression in cold acclimation conditions. This study provides a fundamental database to deepen our knowledge and understanding of regulatory networks in cold acclimation.

Genome-wide characterization of tea plant (Camellia sinensis) Hsf transcription factor family and role of CsHsfA2 in heat tolerance

BACKGROUND

Heat stress factors (Hsfs) play vital roles in signal transduction pathways operating in responses to environmental stresses. However, Hsf gene family has not been thoroughly explored in tea plant (Camellia sinensis L.).

RESULTS

In this study, we identified 25 CsHsf genes in C. sinensis that were separated by phylogenetic analysis into three sub-families (i.e., A, B, and C). Gene structures, conserved domains and motifs analyses indicated that the CsHsf members in each class were relatively conserved. Various cis-acting elements involved in plant growth regulation, hormone responses, stress responses, and light responses were located in the promoter regions of CsHsfs. Furthermore, degradome sequencing analysis revealed that 7 CsHsfs could be targeted by 9 miRNAs. The expression pattern of each CsHsf gene was significantly different in eight tissues. Many CsHsfs were differentially regulated by drought, salt, and heat stresses, as well as exogenous abscisic acid (ABA) and Ca²⁺. In addition, CsHsfA2 was located in the nucleus. Heterologous expression of CsHsfA2 improved thermotolerance in transgenic yeast, suggesting its potential role in the regulation of heat stress response.

CONCLUSIONS

A comprehensive genome-wide analysis of Hsf in C. sinensis present the global identification and functional prediction of CsHsfs. Most of them were implicated in a complex gene regulatory network controlling various abiotic stress responses and signal transduction pathways in tea plants. Additionally, heterologous expression of CsHsfA2 increased thermotolerance of transgenic yeast. These findings provide new insights into the functional divergence of CsHsfs and a basis for further research on CsHsfs functions.

Integrative transcriptome, proteome, and microRNA analysis reveals the effects of nitrogen sufficiency and deficiency conditions on theanine metabolism in the tea plant (Camellia sinensis)

Nitrogen (N) is associated with amino acid metabolism in higher plants. Theanine is an important amino acid in tea plants. To explore the relationship between theanine metabolism and N conditions, we examined the differentially expressed genes (DEGs), proteins (DEPs), and microRNAs (DEMs) involved in theanine metabolism in tea plant shoots and roots under N sufficiency and deficiency conditions. Transcriptome, proteome, and microRNA analyses were performed on tea plant shoots and roots under N sufficiency and deficiency conditions. The contents of theanine, expression levels of genes involved in theanine metabolism, contents of proteinogenic amino acids, and activity of enzymes were analyzed. The DEP-DEG correlation pairs and negative DEM-DEG interactions related to theanine metabolism were identified based on correlation analyses. The expression profiles of DEGs and negative DEM-DEG pairs related to theanine biosynthesis were consistent with the sequencing results. Our results suggest that the molecular and physiological mechanism of theanine accumulation is significantly affected by N sufficiency and deficiency conditions. The DEGs, DEPs, and DEMs and the activity of the enzymes involved in theanine biosynthesis might play vital roles in theanine accumulation under N sufficiency and deficiency conditions in the shoots and roots of tea plants.

Elucidating the regulatory roles of microRNAs in maize (Zea mays L.) leaf growth response to chilling stress

MAIN CONCLUSION: miRNAs control leaf size of maize crop during chilling stress tolerance by regulating developmentally important transcriptional factors and sustaining redox homeostasis of cells. Chilling temperature (0-15 °C) is a major constraint for the cultivation of maize (Zea mays) which inhibits the early growth of maize leading to reduction in leaf size. Growth and development take place in meristem, elongation, and mature zones that are linearly located along the leaf base to tip. To prevent shortening of leaf caused by chilling, this study aims to elucidate the regulatory roles of microRNA (miRNA) genes in the controlling process switching between growth and developmental stages. In this respect, hybrid maize ADA313 seedlings were treated to the chilling temperature which caused 26% and 29% reduction in the final leaf length and a decline in cell production of the fourth leaf. The flow cytometry data integrated with the expression analysis of cell cycle genes indicated that the reason for the decline was a failure proceeding from G2/M rather than G1/S. Through an miRNome analysis of 321 known maize miRNAs, 24, 6, and 20 miRNAs were assigned to putative meristem, elongation, and mature zones, respectively according to their chilling response. To gain deeper insight into decreased cell production, in silico, target prediction analysis was performed for meristem specific miRNAs. Among the miRNAs, miR160, miR319, miR395, miR396, miR408, miR528, and miR1432 were selected for confirming the potential of negative regulation with their predicted targets by qRT-PCR. These findings indicated evidence for improvement of growth and yield under chilling stress of the maize.

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.

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.

Tissue-specific transposon-associated small RNAs in the gymnosperm tree, Norway spruce

BACKGROUND

Small RNAs (sRNAs) are regulatory molecules impacting on gene expression and transposon activity. MicroRNAs (miRNAs) are responsible for tissue-specific and environmentally-induced gene repression. Short interfering RNAs (siRNA) are constitutively involved in transposon silencing across different type of tissues. The male gametophyte in angiosperms has a unique set of sRNAs compared to vegetative tissues, including phased siRNAs from intergenic or genic regions, or epigenetically activated siRNAs. This is contrasted by a lack of knowledge about the sRNA profile of the male gametophyte of gymnosperms.

RESULTS

Here, we isolated mature pollen from male cones of Norway spruce and investigated its sRNA profiles. While 21-nt sRNAs is the major size class of sRNAs in needles, in pollen 21-nt and 24-nt sRNAs are the most abundant size classes. Although the 24-nt sRNAs were exclusively derived from TEs in pollen, both 21-nt and 24-nt sRNAs were associated with TEs. We also investigated sRNAs from somatic embryonic callus, which has been reported to contain 24-nt sRNAs. Our data show that the 24-nt sRNA profiles are tissue-specific and differ between pollen and cell culture.

CONCLUSION

Our data reveal that gymnosperm pollen, like angiosperm pollen, has a unique sRNA profile, differing from vegetative leaf tissue. Thus, our results reveal that angiosperm and gymnosperm pollen produce new size classes not present in vegetative tissues; while in angiosperm pollen 21-nt sRNAs are generated, in the gymnosperm Norway spruce 24-nt sRNAs are generated. The tissue-specific production of distinct TE-derived sRNAs in angiosperms and gymnosperms provides insights into the diversification process of sRNAs in TE silencing pathways between the two groups of seed plants.

Identification of microRNAs and relative target genes in Moringa oleifera leaf and callus

MicroRNAs, a class of small, non-coding RNAs, play important roles in plant growth, development and stress response by negatively regulating gene expression. Moringa oleifera Lam. plant has many medical and nutritional uses; however, little attention has been dedicated to its potential for the bio production of active compounds. In this study, 431 conserved and 392 novel microRNA families were identified and 9 novel small RNA libraries constructed from leaf, and cold stress treated callus, using high-throughput sequencing technology. Based on the M. oleifera genome, the microRNA repertoire of the seed was re-evaluated. qRT-PCR analysis confirmed the expression pattern of 11 conserved microRNAs in all groups. MicroRNA159 was found to be the most abundant conserved microRNA in leaf and callus, while microRNA393 was most abundantly expressed in the seed. The majority of predicted microRNA target genes were transcriptional factors involved in plant reproduction, growth/development and abiotic/biotic stress response. In conclusion, this is the first comprehensive analysis of microRNAs in M. oleifera leaf and callus which represents an important addition to the existing M. oleifera seed microRNA database and allows for possible exploitation of plant microRNAs induced with abiotic stress, as a tool for bio-enrichment with pharmacologically important phytochemicals.

Genome-wide investigation of superoxide dismutase (SOD) gene family and their regulatory miRNAs reveal the involvement in abiotic stress and hormone response in tea plant (Camellia sinensis)

Superoxide dismutases (SODs), as a family of metalloenzymes related to the removal of reactive oxygen species (ROS), have not previously been investigated at genome-wide level in tea plant. In this study, 10 CsSOD genes were identified in tea plant genome, including 7 Cu/Zn-SODs (CSDs), 2 Fe-SODs (FSDs) and one Mn-SOD (MSD), and phylogenetically classified in three subgroups, respectively. Physico-chemical characteristic, conserved motifs and potential protein interaction analyses about CsSOD proteins were carried out. Exon-intron structures and codon usage bias about CsSOD genes were also examined. Exon-intron structures analysis revealed that different CsSOD genes contained various number of introns. On the basis of the prediction of regulatory miRNAs of CsSODs, a modification 5’ RNA ligase-mediated (RLM)-RACE was performed and validated that csn-miR398a-3p-1 directly cleaves CsCSD4. By prediction of cis-acting elements, the expression patterns of 10 CsSOD genes and their regulatory miRNAs were detected under cold, drought, exogenous methyl jasmonate (MeJA) and gibberellin (GA₃) treatments. The results showed that most of CsSODs except for CsFSD2 were induced under cold stress and CsCSDs may play primary roles under drought stress; exogenous GA₃ and MeJA could also stimulated/inhibited distinct CsSODs at different stages. In addition, we found that csn-miR398a-3p-1 negatively regulated the expression of CsCSD4 may be a crucial regulatory mechanism under cold stress. This study provides a certain basis for the studies about stress resistance in tea plants, even provide insight into comprehending the classification, evolution, diverse functions and influencing factors of expression patterns for CsSOD genes.

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.

Identification of Regulatory Networks of MicroRNAs and Their Targets in Response to Colletotrichum gloeosporioides in Tea Plant (Camellia sinensis L.)

Anthracnose disease is caused by Colletotrichum gloeosporioides, and is common in leaves of the tea plant (Camellia sinensis). MicroRNAs (miRNAs) have been known as key modulators of gene expression in response to environmental stresses, disease resistance, defense responses, and plant immunity. However, the role of miRNAs in responses to C. gloeosporioides remains unexplored in tea plant. Therefore, in the present study, six miRNA sequencing data sets and two degradome data sets were generated from C. gloeosporioides-inoculated and control tea leaves. A total of 485 conserved and 761 novel miRNAs were identified. Of those, 239 known and 369 novel miRNAs exhibited significantly differential expression under C. gloeosporioides stress. One thousand one hundred thirty-four and 596 mRNAs were identified as targets of 389 conserved and 299 novel miRNAs by degradome analysis, respectively. Based on degradome analysis, most of the predicted targets are negatively correlated with their corresponding conserved and novel miRNAs. The expression levels of 12 miRNAs and their targets were validated by quantitative real-time PCR. A negative correlation between expression profiles of five miRNAs (PC-5p-80764_22, csn-miR160c, csn-miR828a, csn-miR164a, and csn-miR169e) and their targets (WRKY, ARF, MYB75, NAC, and NFY transcription factor) was observed. The predicted targets of five interesting miRNAs were further validated through 5'RLM-RACE. Furthermore, Gene Ontology and metabolism pathway analysis revealed that most of the target genes were involved in the regulation of auxin pathway, ROS scavenging pathway, salicylic acid mediated pathway, receptor kinases, and transcription factors for plant growth and development as well as stress responses in tea plant against C. gloeosporioides stress. This study enriches the resources of stress-responsive miRNAs and their targets in C. sinensis and thus provides novel insights into the miRNA-mediated regulatory mechanisms, which could contribute to the enhanced susceptibility of C. gloeosporioides in tea plant.

Transcriptomic and metabolomic profiling of Camellia sinensis L. cv. 'Suchazao' exposed to temperature stresses reveals modification in protein synthesis and photosynthetic and anthocyanin biosynthetic pathways

To determine the mechanisms in tea plants responding to temperature stresses (heat and cold), we examined the global transcriptomic and metabolomic profiles of the tea plant cultivar 'Suchazao' under moderately low temperature stress (ML), severely low temperature stress (SL), moderately high temperature stress (MH) and severely high temperature stress (SH) using RNA-seq and high performance liquid chromatography tandem mass spectrometry/mass spectrometry (HPLC-MS/MS), respectively. The identified differentially expressed genes indicated that the synthesis of stress-resistance protein might be redirected to cope with the temperature stresses. We found that heat shock protein genes Hsp90 and Hsp70 played more critical roles in tea plants in adapting to thermal stress than cold, while late embryogenesis abundant protein genes (LEA) played a greater role under cold than heat stress, more types of zinc finger genes were induced under cold stress as well. In addition, energy metabolisms were inhibited by SH, SL and ML. Furthermore, the mechanisms of anthocyanin synthesis were different under the cold and heat stresses. Indeed, the CsUGT75C1 gene, encoding UDP-glucose:anthocyanin 5-O-glucosyl transferase, was up-regulated in the SL-treated leaves but down-regulated in SH. Metabolomics analysis also showed that anthocyanin monomer levels increased under SL. These results indicate that the tea plants share certain foundational mechanisms to adjust to both cold and heat stresses. They also developed some specific mechanisms for surviving the cold or heat stresses. Our study provides effective information about the different mechanisms tea plants employ in surviving cold and heat stresses, as well as the different mechanisms of anthocyanin synthesis, which could speed up the genetic breeding of heat- and cold-tolerant tea varieties.

Comparative transcriptomic analysis reveals gene expression associated with cold adaptation in the tea plant Camellia sinensis

BACKGROUND

Low temperature restricts the planting range of all crops, but cold acclimation induces adaption to cold stress in many plants. Camellia sinensis, a perennial evergreen tree that is the source of tea, is mainly grown in warm areas. Camellia sinensis var. sinensis (CSS) has greater cold tolerance than Camellia sinensis var. assamica (CSA). To gain deep insight into the molecular mechanisms underlying cold adaptation, we investigated the physiological responses and transcriptome profiles by RNA-Seq in two tea varieties, cold resistant SCZ (classified as CSS) and cold susceptible YH9 (classified as CSA), during cold acclimation.

RESULTS

Under freezing stress, lower relative electrical conductivity and higher chlorophyll fluorescence (Fv/Fm) values were detected in SCZ than in YH9 when subjected to freezing acclimation. During cold treatment, 6072 and 7749 DEGs were observed for SCZ and YH9, respectively. A total of 978 DEGs were common for both SCZ and YH9 during the entire cold acclimation process. DEGs were enriched in pathways of photosynthesis, hormone signal transduction, and transcriptional regulation of plant-pathogen interactions. Further analyses indicated that decreased expression of Lhca2 and higher expression of SnRK2.8 are correlated with cold tolerance in SCZ.

CONCLUSIONS

Compared with CSA, CSS was significantly more resistant to freezing after cold acclimation, and this increased resistance was associated with an earlier expression of cold-induced genes. Because the greater transcriptional differentiation during cold acclimation in SCZ may contribute to its greater cold tolerance, our studies identify specific genes involved in photoinhibition, ABA signal conduction, and plant immunity that should be studied for understanding the processes involved in cold tolerance. Marker-assisted breeding focused on the allelic variation at these loci provides an avenue for the possible generation of CSA cultivars that have CSS-level cold tolerance.