Genetic Regulatory Pathways Research Articles

Optimizing nutrient transporters to enhance disease resistance in rice.

Fertilizers and plant diseases contribute positively and negatively to crop production, respectively. Macro- and micronutrients provided by the soil and fertilizers are transported by various plant nutrient transporters from the soil to the roots and shoots, facilitating growth and development. However, the homeostasis of different nutrients has different effects on plant disease. This review is aimed at providing insights into the interconnected regulation between nutrient homeostasis and immune responses, and it highlights strategies to enhance disease resistance by optimal manipulation of nutrient transporters in rice. First, we highlight the essential roles of six macronutrients (nitrogen, phosphorus, potassium, sulfur, calcium, magnesium) and eight micronutrients (iron, manganese, zinc, copper, boron, molybdenum, silicon, nickel), and summarize the diverse effects of each on rice diseases. We then systematically review the molecular mechanisms of immune responses modulated by nutrient transporters and the genetic regulatory pathways that control the specific nutrient-mediated immune signaling that is regulated by the pathogens and the host plant. Finally, we discuss putative strategies for breeding disease-resistant rice by genetic engineering of nutrient transporters.

Salivary gland bioengineering - yesterday, today, tomorrow!

Salivary glands are specialized structures developed as an extensively compact, arborized design through classical embryogenesis, accompanied by a cascade of events channelized by numerous growth factors and genetic regulatory pathways. Salivary secretions maintain oral homeostasis and, when diminished in certain conditions, present as xerostomia or salivary hypofunction, adversely impacting the patient's quality of life. The current available treatments primarily aim at tackling the immediate symptoms providing temporary relief to the patient. Despite scientific efforts to develop permanent and effective solutions to restore salivation, a significant permanent treatment is yet to be established. Tissue engineering has proven as a promising remedial tool in several diseases, as well as in xerostomia, and aims to restore partial loss of organ function. Recapitulating the physiological cellular microenvironment to in vitro culture conditions is constantly evolving. Replicating the dynamic multicellular interactions, genetic pathways, and cytomorphogenic forces, as displayed during salivary gland development have experienced considerable barriers. Through this review, we endeavour to provide an outlook on the evolution of in vitro salivary gland research, highlighting the key bioengineering advances and the challenges faced with the current therapeutic strategies for salivary hypofunction, with an insight into our team's scientific contributions.

Whole-genome sequencing unravels novel genetic determinants and regulatory pathways associated with triamcinolone acetonide-induced ocular hypertension.

Glucocorticosteroids commonly used to treat certain ocular inflammatory conditions cause an unwarranted elevation in intraocular pressure (IOP) leading to steroid-induced ocular hypertension (OHT). This study aims to identify novel genetic variants in the Indian population associated with steroid responsiveness, specifically to that of intravitreal Triamcinolone acetonide (TA) injections, which leads to OHT in 27% of the TA-treated Indian subjects. Genetic determinants and pathways regulating TA-OHT progression were investigated by applying whole-genome sequencing (WGS) on DNA extracted from 53 blood samples that included TA responders and non-responders. Sequencing analysis yielded 45 intronic and 49 exonic variants to be associated with TA-OHT, which are known to play a vital role in eye, heart, brain, and bone deformities. Of these, the most significant genetic variant associated with TA-OHT was further considered for molecular dynamics (MD) simulation studies. Variants in the CRPPA, PLOD1, ARHGAP1, TIMELESS and TNFSF4 genes were found to be directly implicating TA-OHT. Furthermore, these genes were enriched in pathways associated with cardiomyopathy, focal adhesion, extracellular matrix, and actin cytoskeleton reorganization. MD simulation studies revealed that the top significant variant (rs141625803) in the CRPPA gene possesses a high pathogenic and structurally destabilizing effect. Thus, novel genetic variants that could be significantly associated with the TA-OHT progression were identified in this study. Validation of these targets in a larger cohort of patients along with their functional analysis would inform on the disease, thereby adding to the existing knowledge on the pathophysiology of TA-OHT.

ACSL3 is a potential prognostic biomarker for immune infiltration in clear cell renal cell carcinoma.

Long-chain acyl-coenzyme A synthases (ACSLs) catalyze the activation of fatty acid and are often dysregulated in malignancies. The purpose of this research was to figure out the ACSL3's prognostic value and mechanism in clear cell renal cell carcinoma (ccRCC). The expression of ACSL3 in ccRCC was investigated in this work using data from the GEO, TCGA, Oncomine and HPA databases. The expression differences of ACSL3 in the cell lines were further detected by qPCR and Western blot. GEPIA, MethSurv, cBioPortal, and the TIMER were used to perform survival and correlation analysis on ACSL3. GO and KEGG analyses were carried out in R using clusterProfiler and GOplot. Protein-protein interactions (PPI) are constructed from Strings website, and the results were visualized in Cytoscape software. The expression level of ACSL3 was significantly reduced in ccRCC tissues, and its mRNA and protein expression were also significantly lower in both renal cancer cell lines. ACSL3 is significantly related to clinical stage, OS, DFS, DNA methylation, and immune-cell infiltration. Our findings demonstrated that data mining was capable of eliciting information on ACSL3 levels and its role in genetic regulatory pathways in ccRCC.

Effects of Combined Abiotic Stresses Related to Climate Change on Root Growth in Crops

Climate change is a major threat to crop productivity that negatively affects food security worldwide. Increase in global temperatures are usually accompanied by drought, flooding and changes in soil nutrients composition that dramatically reduced crop yields. Against the backdrop of climate change, human population increase and subsequent rise in food demand, finding new solutions for crop adaptation to environmental stresses is essential. The effects of single abiotic stress on crops have been widely studied, but in the field abiotic stresses tend to occur in combination rather than individually. Physiological, metabolic and molecular responses of crops to combined abiotic stresses seem to be significantly different to individual stresses. Although in recent years an increasing number of studies have addressed the effects of abiotic stress combinations, the information related to the root system response is still scarce. Roots are the underground organs that directly contact with the soil and sense many of these abiotic stresses. Understanding the effects of abiotic stress combinations in the root system would help to find new breeding tools to develop more resilient crops. This review will summarize the current knowledge regarding the effects of combined abiotic stress in the root system in crops. First, we will provide a general overview of root responses to particular abiotic stresses. Then, we will describe how these root responses are integrated when crops are challenged to the combination of different abiotic stress. We will focus on the main changes on root system architecture (RSA) and physiology influencing crop productivity and yield and convey the latest information on the key molecular, hormonal and genetic regulatory pathways underlying root responses to these combinatorial stresses. Finally, we will discuss possible directions for future research and the main challenges needed to be tackled to translate this knowledge into useful tools to enhance crop tolerance.

Exploring epigenetic reprogramming during central nervous system infection.

Epigenetics involves the study of various modes of adaptable transcriptional regulation, contributing to cell identity, characteristics, and function. During central nervous system (CNS) infection, epigenetic mechanisms can exert pronounced control over the maturation and antimicrobial properties of nearly every immune cell type. Epigenetics is a relatively new field, with the first mention of these marks proposed only a half-century ago and a substantial body of immunological epigenetic research emerging only in the last few decades. Here, we review the best-characterized epigenetic marks and their functions as well as illustrate how various immune cell populations responding to CNS infection utilize these marks to organize their activation state and inflammatory processes. We also discuss the metabolic and clinical implications of epigenetic marks and the rapidly expanding set of tools available to researchers that are enabling elucidation of increasingly detailed genetic regulatory pathways. These considerations paint an intricate picture of inflammatory regulation, where epigenetic marks influence genetic, signaling, and environmental elements to orchestrate a tailored immunological response to the threat at hand, cementing epigenetics as an important player in immunity.

Intergenerational changes in hippocampal transcription in an animal model of maternal depression.

Chronic stress during early life, such as exposure to social conflict or deficits in parental care, can have persistent adverse behavioural effects. Offspring in a rodent model of maternal depression and early life stress have increased susceptibility to maternal depression themselves, suggesting a pathway by which maternal stress could be intergenerationally inherited. The overall aim of this study was to explore the genetic regulatory pathways underlying how maternal social stress and reduced care mediates stress-related behavioural changes in offspring across generations. This study investigated a social stress-based rat model of postpartum depression and the intergenerational inheritance of depressed maternal care where F0 (dams exposed to male intruder stress during lactation) and F1 offspring are directly exposed to social stress. RNASeq was used to investigate genome-wide transcriptome changes in the hippocampus of F1 and F2 generations. Transcriptome analyses revealed differential expression of 69 genes in the F1 generation and 14 in the F2 between controls versus social stress differences. Many of these genes were receptors and calcium-binding proteins in the F1 and involved in cellular oxidant detoxification in F2. The present data identify and characterize changes in the neural expression of key genes involved in the regulation of depression maintained between the generations, suggesting a potential neural pathway for the intergenerational transmission of depressed maternal care and maternal anxiety in the CSS model. Further work is needed to understand to what extent these results are due to molecular germline inheritance and/or the social propagation of deficits in maternal care.

Genome-wide association studies reveal the coordinated regulatory networks underlying photosynthesis and wood formation in Populus.

Photosynthesis and wood formation underlie the ability of trees to provide renewable resources and perform ecological functions; however, the genetic basis and regulatory pathways coordinating these two linked processes remain unclear. Here, we used a systems genetics strategy, integrating genome-wide association studies, transcriptomic analyses, and transgenic experiments, to investigate the genetic architecture of photosynthesis and wood properties among 435 unrelated individuals of Populus tomentosa, and unravel the coordinated regulatory networks resulting in two trait categories. We detected 222 significant single-nucleotide polymorphisms, annotated to 177 candidate genes, for 10 traits of photosynthesis and wood properties. Epistasis uncovered 74 epistatic interactions for phenotypes. Strikingly, we deciphered the coordinated regulation patterns of pleiotropic genes underlying phenotypic variations for two trait categories. Furthermore, expression quantitative trait nucleotide mapping and coexpression analysis were integrated to unravel the potential transcriptional regulatory networks of candidate genes coordinating photosynthesis and wood properties. Finally, heterologous expression of two pleiotropic genes, PtoMYB62 and PtoMYB80, in Arabidopsis thaliana demonstrated that they control regulatory networks balancing photosynthesis and stem secondary cell wall components, respectively. Our study provides insights into the regulatory mechanisms coordinating photosynthesis and wood formation in poplar, and should facilitate genetic breeding in trees via molecular design.

Congenital Malformations in Sea Turtles: Puzzling Interplay between Genes and Environment.

Simple SummaryCongenital malformations can lead to embryonic mortality in many species, and sea turtles are no exception. Genetic and/or environmental alterations occur during early development in the embryo, and may produce aberrant phenotypes, many of which are incompatible with life. Causes of malformations are multifactorial; genetic factors may include mutations, chromosomal aberrations, and inbreeding effects, whereas non-genetic factors may include nutrition, hyperthermia, low moisture, radiation, and contamination. It is possible to monitor and control some of these factors (such as temperature and humidity) in nesting beaches, and toxic compounds in feeding areas, which can be transferred to the embryo through their lipophilic properties. In this review, we describe possible causes of different types of malformations observed in sea turtle embryos, as well as some actions that may help reduce embryonic mortality.The completion of embryonic development depends, in part, on the interplay between genetic factors and environmental conditions, and any alteration during development may affect embryonic genetic and epigenetic regulatory pathways leading to congenital malformations, which are mostly incompatible with life. Oviparous reptiles, such as sea turtles, that produce numerous eggs in a clutch that is buried on the beach provide an opportunity to study embryonic mortality associated with malformations that occur at different times during development, or that prevent the hatchling from emerging from the nest. In sea turtles, the presence of congenital malformations frequently leads to mortality. A few years ago, a detailed study was performed on external congenital malformations in three species of sea turtles from the Mexican Pacific and Caribbean coasts, the hawksbill turtle, Eretmochelys imbricata (n = 23,559 eggs), the green turtle, Chelonia mydas (n = 17,690 eggs), and the olive ridley, Lepidochelys olivacea (n = 20,257 eggs), finding 63 types of congenital malformations, of which 38 were new reports. Of the three species, the olive ridley showed a higher incidence of severe anomalies in the craniofacial region (49%), indicating alterations of early developmental pathways; however, several malformations were also observed in the body, including defects in the carapace (45%) and limbs (33%), as well as pigmentation disorders (20%), indicating that deviations occurred during the middle and later stages of development. Although intrinsic factors (i.e., genetic mutations or epigenetic modifications) are difficult to monitor in the field, some environmental factors (such as the incubation temperature, humidity, and probably the status of feeding areas) are, to some extent, less difficult to monitor and/or control. In this review, we describe the aetiology of different malformations observed in sea turtle embryos, and provide some actions that can reduce embryonic mortality.

Covid-19 pathogenesis in prostatic cancer and TMPRSS2-ERG regulatory genetic pathway.

Members of Coronaviridae family have been the source of respiratory illnesses. The outbreak of SARS-CoV-2 that produced a severe lung disease in afflicted patients in China and other countries was the reason for the incredible attention paid toward this viral infection. It is known that SARS-CoV-2 is dependent on TMPRSS2 activity for entrance and subsequent infection of the host cells and TMPRSS2 is a host cell molecule that is important for the spread of viruses such as coronaviruses.Different factors can increase the risk of prostate cancer, including older age, a family history of the disease. Androgen receptor (AR) initiates a transcriptional cascade which plays a serious role in both normal and malignant prostate tissues. TMPRSS2 protein is highly expressed in prostate secretory epithelial cells, and its expression is dependent on androgen signals. One of the molecular signs of prostate cancer is TMPRSS2-ERG gene fusion. In TMPRSS2-ERG-positive prostate cancers different patterns of changed gene expression can be detected. The possible molecular relation between fusion positive prostate cancer patients and the increased risk of lethal respiratory viral infections especially SARS-CoV-2 can candidate TMPRSS2 as an attractive drug target. The studies show that some molecules such as nicotinamide, PARP1, ETS and IL-1R can be studied deeper in order to control SARS-CoV-2 infection especially in prostate cancer patients.This review attempts to investigate the possible relation between the gene expression pattern that is produced through TMPRSS2-ERG fusion positive prostate cancer and the possible influence of these fluctuations on the pathogenesis and development of viral infections such as SARS-CoV-2.

Quantifying the distribution of protein oligomerization degree reflects cellular information capacity

The generation of information, energy and biomass in living cells involves integrated processes that optimally evolve into complex and robust cellular networks. Protein homo-oligomerization, which is correlated with cooperativity in biology, is one means of scaling the complexity of protein networks. It can play critical roles in determining the sensitivity of genetic regulatory circuits and metabolic pathways. Therefore, understanding the roles of oligomerization may lead to new approaches of probing biological functions. Here, we analyzed the frequency of protein oligomerization degree in the cell proteome of nine different organisms, and then, we asked whether there are design trade-offs between protein oligomerization, information precision and energy costs of protein synthesis. Our results indicate that there is an upper limit for the degree of protein oligomerization, possibly because of the trade-off between cellular resource limitations and the information precision involved in biochemical reaction networks. These findings can explain the principles of cellular architecture design and provide a quantitative tool to scale synthetic biological systems.

Up-regulation of autophagy by low concentration of salicylic acid delays methyl jasmonate-induced leaf senescence

Crosstalk between salicylic acid (SA) and jasmonic acid (JA) signaling plays an important role in regulation of plant senescence. Our previous work found that SA could delay methyl jasmonate (MeJA)-induced leaf senescence in a concentration-dependent manner. Here, the effect of low concentration of SA (LCSA) application on MeJA-induced leaf senescence was further assessed. High-throughput sequencing (RNA-Seq) results showed that LCSA did not have dominant effects on the genetic regulatory pathways of basal metabolism like nitrogen metabolism, photosynthesis and glycolysis. The ClusterONE was applied to identify discrete gene modules based on protein–protein interaction (PPI) network. Interestingly, an autophagy-related (ATG) module was identified in the differentially expressed genes (DEGs) that exclusively induced by MeJA together with LCSA. RT-qPCR confirmed that the expression of most of the determined ATG genes were upregulated by LCSA. Remarkably, in contrast to wild type (Col-0), LCSA cannot alleviate the leaf yellowing phenotype in autophagy defective mutants (atg5-1 and atg7-2) upon MeJA treatment. Confocal results showed that LCSA increased the number of autophagic bodies accumulated in the vacuole during MeJA-induced leaf senescence. Collectively, our work revealed up-regulation of autophagy by LCSA as a key regulator to alleviate MeJA-induced leaf senescence.

Zróżnicowane znaczenie witaminy E, jej występowanie oraz funkcje regulacyjne w różnych tkankach roślinnych

Fruits and vegetables contain vitamins, dietary fibre and phytochemicals, which promote human health. Specifically vitamins are major bioactive compounds. In plant raw materials there are water-soluble (vitamins B and C) and lipid-soluble (vitamins A, E and K) ones, having strong antioxidant potential and limiting several diseases, including cancer and cardio-vascular diseases. Furthermore, vitamins play a key role in plant protection against abiotic and biotic stress factors. Vitamin E is a group of eight fat-soluble compounds known as tocochromanols, which are synthesized exclusively by photosynthetic organisms. Tocochromanols are categorised into two groups – tocopherols and tocotrienols. Each group contains four forms identified by prefixes α-, β-, γ- and δ-, yielding a total of eight forms. Vitamin E is an essential, lipid-soluble antioxidant in the human diet, protecting from cancer, diabetes, cardiovascular and neurological diseases. This review summarizes the recent developments in the understanding the variety of roles of vitamin E as well as its major genetic regulatory pathways in plants. In addition, a comprehensive mapping of vitamin E occurrence is presented in different plant tissues, organelles and horticultural crop species, as well as throughout different stages of fruit development and after postharvest.

Genetic regulatory pathways of split-hand/foot malformation.

Split-hand/foot malformation (SHFM) is caused by mutations in TP63, DLX5, DLX6, FGF8, FGFR1, WNT10B, and BHLHA9. The clinical features of SHFM caused by mutations of these genes are not distinguishable. This implies that in normal situations these SHFM-associated genes share an underlying regulatory pathway that is involved in the development of the central parts of the hands and feet. The mutations in SHFM-related genes lead to dysregulation of Fgf8 in the central portion of the apical ectodermal ridge (AER) and subsequently lead to misexpression of a number of downstream target genes, failure of stratification of the AER, and thus SHFM. Syndactyly of the remaining digits is most likely the effects of dysregulation of Fgf-Bmp-Msx signaling on apoptotic cell death. Loss of digit identity in SHFM is hypothesized to be the effects of misexpression of HOX genes, abnormal SHH gradient, or the loss of balance between GLI3A and GLI3R. Disruption of canonical and non-canonical Wnt signaling is involved in the pathogenesis of SHFM. Whatever the causative genes of SHFM are, the mutations seem to lead to dysregulation of Fgf8 in AER cells of the central parts of the hands and feet and disruption of Wnt-Bmp-Fgf signaling pathways in AER.

Highlights from Recent Cancer Literature

Highlights from Recent Cancer Literature

Gene Silencing In Vitro and In Vivo Using Intronic MicroRNAs.

MicroRNAs (miRNAs), small single-stranded regulatory RNAs capable of interfering with intracellular messenger RNAs (mRNAs) that contain either complete or partial complementarity, are useful for the design of new therapies against cancer polymorphism and viral mutation. Numerous miRNAs have been reported to induce RNA interference (RNAi), a post-transcriptional gene-silencing mechanism. Recent evidence also indicates that they are involved in the transcriptional regulation of genome activities. They were first discovered in Caenorhabditis elegans as native RNA fragments that modulate a wide range of genetic regulatory pathways during embryonic development, and are now recognized as small gene silencers transcribed from the noncoding regions of a genome. In humans, nearly 97% of the genome is noncoding DNA, which varies from one individual to another, and changes in these sequences are frequently noted to manifest in clinical and circumstantial malfunction; for example, type 2 myotonic dystrophy and fragile X syndrome were found to be associated with miRNAs derived from introns. Intronic miRNA is a new class of miRNAs derived from the processing of non-protein-coding regions of gene transcripts. The intronic miRNAs differ uniquely from previously described intergenic miRNAs in the requirement of RNA polymerase (Pol)-II and spliceosomal components for its biogenesis. Several kinds of intronic miRNAs have been identified in C. elegans, mouse, and human cells; however, their functions and applications have not been reported. Here, we show for the first time that intron-derived miRNA is not only able to induce RNAi in mammalian cells but also in fish, chicken embryos, and adult mice cells, demonstrating the evolutionary preservation of this gene regulation system in vivo. These miRNA-mediated animal models provide artificial means to reproduce the mechanisms of miRNA-induced disease in vivo and will shed further light on miRNA-related therapies.

Combined analysis of mRNA and miRNA identifies dehydration and salinity responsive key molecular players in citrus roots

Citrus is one of the most economically important fruit crops around world. Drought and salinity stresses adversely affected its productivity and fruit quality. However, the genetic regulatory networks and signaling pathways involved in drought and salinity remain to be elucidated. With RNA-seq and sRNA-seq, an integrative analysis of miRNA and mRNA expression profiling and their regulatory networks were conducted using citrus roots subjected to dehydration and salt treatment. Differentially expressed (DE) mRNA and miRNA profiles were obtained according to fold change analysis and the relationships between miRNAs and target mRNAs were found to be coherent and incoherent in the regulatory networks. GO enrichment analysis revealed that some crucial biological processes related to signal transduction (e.g. ‘MAPK cascade’), hormone-mediated signaling pathways (e.g. abscisic acid- activated signaling pathway’), reactive oxygen species (ROS) metabolic process (e.g. ‘hydrogen peroxide catabolic process’) and transcription factors (e.g., ‘MYB, ZFP and bZIP’) were involved in dehydration and/or salt treatment. The molecular players in response to dehydration and salt treatment were partially overlapping. Quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) analysis further confirmed the results from RNA-seq and sRNA-seq analysis. This study provides new insights into the molecular mechanisms how citrus roots respond to dehydration and salt treatment.

A New Mechanism for Mendelian Dominance in Regulatory Genetic Pathways: Competitive Binding by Transcription Factors.

We report a new mechanism for allelic dominance in regulatory genetic interactions that we call binding dominance. We investigated a biophysical model of gene regulation, where the fractional occupancy of a transcription factor (TF) on the cis-regulated promoter site it binds to is determined by binding energy (-ΔG) and TF dosage. Transcription and gene expression proceed when the TF is bound to the promoter. In diploids, individuals may be heterozygous at the cis-site, at the TF's coding region, or at the TF's own promoter, which determines allele-specific dosage. We find that when the TF's coding region is heterozygous, TF alleles compete for occupancy at the cis-sites and the tighter-binding TF is dominant in proportion to the difference in binding strength. When the TF's own promoter is heterozygous, the TF produced at the higher dosage is also dominant. Cis-site heterozygotes have additive expression and therefore codominant phenotypes. Binding dominance propagates to affect the expression of downstream loci and it is sensitive in both magnitude and direction to genetic background, but its detectability often attenuates. While binding dominance is inevitable at the molecular level, it is difficult to detect in the phenotype under some biophysical conditions, more so when TF dosage is high and allele-specific binding affinities are similar. A body of empirical research on the biophysics of TF binding demonstrates the plausibility of this mechanism of dominance, but studies of gene expression under competitive binding in heterozygotes in a diversity of genetic backgrounds are needed.

Expression Profiling of Circulating MicroRNAs in Canine Myxomatous Mitral Valve Disease.

MicroRNAs (miRNAs) are small non-coding RNAs that have shown promise as noninvasive biomarkers in cardiac disease. This study was undertaken to investigate the miRNA expression profile in dogs with myxomatous mitral valve disease (MMVD). 277 miRNAs were quantified using RT-qPCR from six normal dogs (American College of Veterinary Internal Medicine Stage A), six dogs with MMVD mild to moderate cardiac enlargement (ACVIM Stage B1/B2) and six dogs with MMVD and congestive heart failure (ACVIM Stage C/D). Eleven miRNAs were differentially expressed (False Discovery Rate < 0.05). Dogs in Stage B1/B2 or C/D had four upregulated miRNAs, including three cfa-let-7/cfa-miR-98 family members, while seven others were downregulated, compared to Stage A. Expression of six of the 11 miRNAs also were significantly different between dogs in Stage C/D and those in Stage B1/B2. The expression changes were greater as disease severity increased. These miRNAs may be candidates for novel biomarkers and may provide insights into genetic regulatory pathways in canine MMVD.

Integrated RNA-Seq and sRNA-Seq Analysis Identifies Chilling and Freezing Responsive Key Molecular Players and Pathways in Tea Plant (Camellia sinensis).

Tea [Camellia sinensis (L) O. Kuntze, Theaceae] is one of the most popular non-alcoholic beverages worldwide. Cold stress is one of the most severe abiotic stresses that limit tea plants’ growth, survival and geographical distribution. However, the genetic regulatory network and signaling pathways involved in cold stress responses in tea plants remain unearthed. Using RNA-Seq, DGE and sRNA-Seq technologies, we performed an integrative analysis of miRNA and mRNA expression profiling and their regulatory network of tea plants under chilling (4℃) and freezing (-5℃) stress. Differentially expressed (DE) miRNA and mRNA profiles were obtained based on fold change analysis, miRNAs and target mRNAs were found to show both coherent and incoherent relationships in the regulatory network. Furthermore, we compared several key pathways (e.g., ‘Photosynthesis’), GO terms (e.g., ‘response to karrikin’) and transcriptional factors (TFs, e.g., DREB1b/CBF1) which were identified as involved in the early chilling and/or freezing response of tea plants. Intriguingly, we found that karrikins, a new group of plant growth regulators, and β-primeverosidase (BPR), a key enzyme functionally relevant with the formation of tea aroma might play an important role in both early chilling and freezing response of tea plants. Quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) analysis further confirmed the results from RNA-Seq and sRNA-Seq analysis. This is the first study to simultaneously profile the expression patterns of both miRNAs and mRNAs on a genome-wide scale to elucidate the molecular mechanisms of early responses of tea plants to cold stress. In addition to gaining a deeper insight into the cold resistant characteristics of tea plants, we provide a good case study to analyse mRNA/miRNA expression and profiling of non-model plant species using next-generation sequencing technology.