Small RNAs as big players in plant abiotic stress responses and nutrient deprivation

In Situ Localization of Small RNAs in Plants

Small RNAs have vital roles in numerous aspects of plant biology. Deciphering their precise contributions requires knowledge of a small RNA's spatiotemporal pattern of accumulation. The in situ hybridization protocol described here takes advantage of locked nucleic acid (LNA) oligonucleotide probes to visualize small RNA expression at the cellular level with high sensitivity and specificity. The procedure is optimized for paraffin-embedded plant tissue sections, is applicable to a wide range of plants and tissues, and can be completed within 2-6 days.

Identification of miRNAs and their target genes in Larix olgensis and verified of differential expression miRNAs

BACKGROUND

MiRNAs (microRNA) are 18-24 nt endogenous noncoding RNAs that regulate gene expression at the post-transcriptional level, including tissue-specific, developmental timing and evolutionary conservation gene expression.

RESULTS

This study used high-throughput sequencing technology for the first time in Larix olgensis, predicted 78 miRNAs, including 12,229,003 reads sRNA, screened differentially expressed miRNAs. Predicting target genes was helpful for understanding the miRNA regulation function and obtained 333 corresponding target genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation were analysed, mostly including nucleic acid binding, plant hormone signal transduction, pantothenate and CoA biosynthesis, and cellulose synthase. This study will lay the foundation for clarifying the complex miRNA-mediated regulatory network for growth and development. In view of this, spatio-temporal expression of miR396, miR950, miR164, miR166 and miR160 were analysed in Larix olgensis during the growth stages of not lignified, beginning of lignification, and completely lignified in different tissues (root, stem, and leaf) by quantitative real-time PCR (qRT-PCR). There were differences in the expression of miRNAs in roots, stems and leaves in the same growth period. At 60 days, miR160, miR166 and miR396-2 exhibited the highest expression in leaves. At 120 days, most miRNAs in roots and stems decreased significantly. At 180 days, miRNAs were abundantly expressed in roots and stems. Meanwhile, analysis of the expression of miRNAs in leaves revealed that miR396-2 was reduced as time went on, whereas other miRNAs increased initially and then decreased. On the other hand, in the stems, miR166-1 was increase, whereas other miRNAs, especially miR160, miR164, miR396 and miR950-1, first decreased and then increased. Similarly, in the roots, miR950-2 first decreased and then increased, whereas other miRNAs exhibited a trend of continuous increase.

CONCLUSIONS

The present investigation included rapid isolation and identification of miRNAs in Larix olgensis through construction of a sRNA library using Solexa and predicted 78 novel miRNAs, which showed differential expression levels in different tissues and stages. These results provided a theoretical basis for further revealing the genetic regulation mechanism of miRNA in the growth and development of conifers and the verification of function in target genes.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Conserved miR396b-GRF Regulation Is Involved in Abiotic Stress Responses in Pitaya (Hylocereus polyrhizus)

MicroRNA396 (miR396) is a conserved microRNA family that targets growth-regulating factors (GRFs), which play significant roles in plant growth and stress responses. Available evidence justifies the idea that miR396-targeted GRFs have important functions in many plant species; however, no genome-wide analysis of the pitaya (Hylocereus polyrhizus) miR396 gene has yet been reported. Further, its biological functions remain elusive. To uncover the regulatory roles of miR396 and its targets, the hairpin sequence of pitaya miR396b and the open reading frame (ORF) of its target, HpGRF6, were isolated from pitaya. Phylogenetic analysis showed that the precursor miR396b (MIR396b) gene of plants might be clustered into three major groups, and, generally, a more recent evolutionary relationship in the intra-family has been demonstrated. The sequence analysis indicated that the binding site of hpo-miR396b in HpGRF6 is located at the conserved motif which codes the conserved "RSRKPVE" amino acid in the Trp-Arg-Cys (WRC) region. In addition, degradome sequencing analysis confirmed that four GRFs (GRF1, c56908.graph_c0; GRF4, c52862.graph_c0; GRF6, c39378.graph_c0 and GRF9, c54658.graph_c0) are hpo-miR396b targets that are regulated by specific cleavage at the binding site between the 10th and 11th nucleotides from the 5' terminus of hpo-miR396b. Furthermore, quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that hpo-miR396b is down-regulated when confronted with drought stress (15% polyethylene glycol, PEG), and its expression fluctuates under other abiotic stresses, i.e., low temperature (4 ± 1 °C), high temperature (42 ± 1 °C), NaCl (100 mM), and abscisic acid (ABA; 0.38 mM). Conversely, the expression of HpGRF6 showed the opposite trend to exposure to these abiotic stresses. Taken together, hpo-miR396b plays a regulatory role in the control of HpGRF6, which might influence the abiotic stress response of pitaya. This is the first documentation of this role in pitaya and improves the understanding of the molecular mechanisms underlying the tolerance to drought stress in this fruit.

Identification and Characterization of Salt-Responsive MicroRNAs in Vicia faba by High-Throughput Sequencing

Salt stress has detrimental effects on plant growth and development. MicroRNAs (miRNAs) are a class of noncoding RNAs that are involved in post-transcriptional gene expression regulation. In this study, small RNA sequencing was employed to identify the salt stress-responsive miRNAs of the salt-sensitive Hassawi-3 and the salt-tolerant ILB4347 genotypes of faba bean, growing under salt stress. A total of 527 miRNAs in Hassawi-3 plants, and 693 miRNAs in ILB4347 plants, were found to be differentially expressed. Additionally, 284 upregulated and 243 downregulated miRNAs in Hassawi-3, and 298 upregulated and 395 downregulated miRNAs in ILB4347 plants growing in control and stress conditions were recorded. Target prediction and annotation revealed that these miRNAs regulate specific salt-responsive genes, which primarily included genes encoding transcription factors and laccases, superoxide dismutase, plantacyanin, and F-box proteins. The salt-responsive miRNAs and their targets were functionally enriched by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, which showed that the miRNAs were involved in salt stress-related biological pathways, including the ABC transporter pathway, MAPK signaling pathway, plant hormone signal transduction, and the phosphatidylinositol signaling system, among others, suggesting that the miRNAs play an important role in the salt stress tolerance of the ILB4347 genotype. These results offer a novel understanding of the regulatory role of miRNAs in the salt response of the salt-tolerant ILB4347 and the salt-sensitive Hassawi-3 faba bean genotypes.

Quantitative expression of microRNAs in Brassica oleracea infected with Xanthomonas campestris pv. campestris

Brassica oleracea var. capitata (cabbage) is an economically important crop affected by black rot disease caused by Xanthomonas campestris pv. campestris (Xcc). MicroRNAs (miRNAs) play an important role in plant defense modulation and therefore the analysis of these molecules can help better understand plant-pathogen interactions. In this study, we report the differential expression of four miRNAs that seem to participate in the plant response to Xcc infection. Northern Blot and RT-qPCR techniques were used to measure miRNA expression in resistant (União) and susceptible (Kenzan) cultivars. From 6 miRNAs analyzed, 4 were detected and differentially expressed, showing a down- and upregulated expression profile in susceptible and resistant cultivars, respectively. These results suggest that miR156, miR167, miR169, and miR390 could play a role in B. oleracea resistance enhancement against Xcc and could be explored as potential resistance markers in B. oleracea-Xcc interaction.

Transcriptome-wide identification of miRNA targets and a TAS3-homologous gene in Populus by degradome sequencing

BACKGROUND

Degradome sequencing has been applied to identify miRNA-directed mRNA cleavage and understand the biological function of miRNAs and their target genes in plants defense to stress. miRNAs involved in the response to cold stress have been identified in Populus, however, there are few reports about the validated targets of miRNAs in Populus under cold stress.

OBJECTIVES

The primary objective of this investigation was to globally identify and validate the targets of the miRNAs and regulatory components in Populus under cold stress.

METHODS

Populus plantlets grown in vitro were treated with cold (4 °C for 8 h) and total RNA was extracted using Trizol reagent. Approximately 200 µg total RNA was used for the construction of the degradome library, and degradome sequencing was conducted on an Illumina HiSeq 2000. The sequences were mapped to Populus genome using SOAP 2.0 and then were collected for degradome analysis. Additionally, trans-acting siRNA sequences from transacting siRNA gene 3 sequences and mature miRNAs cleaved from precursor miRNAs of Populus were analyzed. 5' RNA ligase-mediated-RACE (5'-RACE) were further conducted.

RESULTS

80 genes were experimentally determined to be the target of 51 unique miRNAs, including three down-regulated miRNAs (pto-miR156k, pto-miR169i-m, and pto-miR394a-5p/b-5p) and two up-regulated miRNAs (pto-miR167a-d and pto-miR167f/g). The specificity and diversity of cleavage sites of miRNA targets were validated through 5'-RACE experiment and the results were similar with that of degradome sequencing, further supporting the empirical cleavage of miRNAs on targets in vivo in Populus. Interestingly, the TAS-homologous gene pto-TAS3 (EF146176.1) was identified and 11 potential ta-siRNAs [D1(+)-D11(+)] and their possible biogenesis sites within the pto-TAS3 transcript sequence were predicted in Populus. In addition, the biosynthesis of miRNA from precursor miRNA (pre-miRNA) was also validated through the detection of a total of 17 pre-miRNAs.

CONCLUSION

Our investigation expands the application of degradome sequencing for evaluating miRNA regulatory elements and evidence of the miRNA synthesis process, and provides empirical evidence of bona fide cleavage of target genes by miRNAs in Populus, which might be used for the research of miRNA-mediated regulation mechanism and molecular improvement of resistance to cold stress.

Profiling the Abiotic Stress Responsive microRNA Landscape of Arabidopsis thaliana

It is well established among interdisciplinary researchers that there is an urgent need to address the negative impacts that accompany climate change. One such negative impact is the increased prevalence of unfavorable environmental conditions that significantly contribute to reduced agricultural yield. Plant microRNAs (miRNAs) are key gene expression regulators that control development, defense against invading pathogens and adaptation to abiotic stress. Arabidopsis thaliana (Arabidopsis) can be readily molecularly manipulated, therefore offering an excellent experimental system to alter the profile of abiotic stress responsive miRNA/target gene expression modules to determine whether such modification enables Arabidopsis to express an altered abiotic stress response phenotype. Towards this goal, high throughput sequencing was used to profile the miRNA landscape of Arabidopsis whole seedlings exposed to heat, drought and salt stress, and identified 121, 123 and 118 miRNAs with a greater than 2-fold altered abundance, respectively. Quantitative reverse transcriptase polymerase chain reaction (RT-qPCR) was next employed to experimentally validate miRNA abundance fold changes, and to document reciprocal expression trends for the target genes of miRNAs determined abiotic stress responsive. RT-qPCR also demonstrated that each miRNA/target gene expression module determined to be abiotic stress responsive in Arabidopsis whole seedlings was reflective of altered miRNA/target gene abundance in Arabidopsis root and shoot tissues post salt stress exposure. Taken together, the data presented here offers an excellent starting platform to identify the miRNA/target gene expression modules for future molecular manipulation to generate plant lines that display an altered response phenotype to abiotic stress.

Genome-wide comparative analysis in Solanaceous species reveals evolution of microRNAs targeting defense genes in Capsicum spp

MicroRNAs (miRNAs) play roles in various biological processes in plants including growth, development, and disease resistance. Previous studies revealed that some plant miRNAs produce secondary small interfering RNAs (siRNAs) such as phased, secondary siRNAs (phasiRNAs), and they regulate a cascade of gene expression. We performed a genome-wide comparative analysis of miRNAs in Solanaceous species (pepper, tomato, and potato), from an evolutionary perspective. Microsynteny of miRNAs was analysed based on the genomic loci and their flanking genes and most of the well-conserved miRNA genes maintained microsynteny in Solanaceae. We identified target genes of the miRNAs via degradome analysis and found that several miRNAs target many genes encoding nucleotide-binding leucine-rich repeat (NLR) or receptor-like proteins (RLPs), which are known to be major players in defense responses. In addition, disease-resistance-associated miRNAs trigger phasiRNA production in pepper, indicating amplification of the regulation of disease-resistance gene families. Among these, miR-n033a-3p, whose target NLRs have been duplicated in pepper, targets more NLRs belonging to specific subgroup in pepper than those in potato. miRNAs targeting resistance genes might have evolved to regulate numerous targets in Solanaceae, following expansion of target resistance genes. This study provides an insight into evolutionary relationship between miRNAs and their target defense genes in plants.

Inferring the regulatory network of the miRNA-mediated response to biotic and abiotic stress in melon

BACKGROUND

MiRNAs have emerged as key regulators of stress response in plants, suggesting their potential as candidates for knock-in/out to improve stress tolerance in agricultural crops. Although diverse assays have been performed, systematic and detailed studies of miRNA expression and function during exposure to multiple environments in crops are limited.

RESULTS

Here, we present such pioneering analysis in melon plants in response to seven biotic and abiotic stress conditions. Deep-sequencing and computational approaches have identified twenty-four known miRNAs whose expression was significantly altered under at least one stress condition, observing that down-regulation was preponderant. Additionally, miRNA function was characterized by high scale degradome assays and quantitative RNA measurements over the intended target mRNAs, providing mechanistic insight. Clustering analysis provided evidence that eight miRNAs showed a broad response range under the stress conditions analyzed, whereas another eight miRNAs displayed a narrow response range. Transcription factors were predominantly targeted by stress-responsive miRNAs in melon. Furthermore, our results show that the miRNAs that are down-regulated upon stress predominantly have as targets genes that are known to participate in the stress response by the plant, whereas the miRNAs that are up-regulated control genes linked to development.

CONCLUSION

Altogether, this high-resolution analysis of miRNA-target interactions, combining experimental and computational work, Illustrates the close interplay between miRNAs and the response to diverse environmental conditions, in melon.

Nutrient Deficiency Tolerance in Citrus Is Dependent on Genotype or Ploidy Level

Plants require essential minerals for their growth and development that are mainly acquired from soil by their roots. Nutrient deficiency is an environmental stress that can seriously affect fruit production and quality. In citrus crops, rootstock/scion combinations are frequently employed to enhance tolerance to various abiotic stresses. These tolerances can be improved in doubled diploid genotypes. The aim of this work was to compare the impact of nutrient deficiency on the physiological and biochemical response of diploid (2x) and doubled diploid (4x) citrus seedlings: Volkamer lemon, Trifoliate orange × Cleopatra mandarin hybrid, Carrizo citrange, Citrumelo 4475. Flhorag1 (Poncirus trifoliata + and willow leaf mandarin), an allotetraploid somatic hybrid, was also included in this study. Our results showed that depending on the genotype, macronutrient and micronutrient deficiency affected certain physiological traits and oxidative metabolism differently. Tetraploid genotypes, mainly Flhorag1 and Citrumelo 4475, appeared resistant compared to the other genotypes as indicated by the lesser decrease in photosynthetic parameters (P net, F v/F m, and G s) and the lower accumulation of oxidative markers (MDA and H2O2) in roots and leaves, especially after long-term nutrient deficiency. Their higher tolerance to nutrient deficiency could be explained by better activation of their antioxidant system. For the other genotypes, tetraploidization did not induce greater tolerance to nutrient deficiency.

Genetic Engineering of the Biosynthesis of Glycine Betaine Modulates Phosphate Homeostasis by Regulating Phosphate Acquisition in Tomato

Glycine betaine (GB), as a putative compatible substance, protects plants against the damaging effects of abiotic stresses. Phosphorus deficiency is one type of abiotic stress that is detrimental to plant growth. Maintenance of phosphate (Pi) homeostasis is crucial. This study demonstrates GB-regulated phosphate homeostasis in the tomato (Solanum lycopersicum cv. 'Moneymaker') transformed with the choline oxidase gene codA from Arthrobacter globiformis. The codA-transgenic lines displayed more resistance to low-phosphate stress. The data revealed that the wild-type plants were stunted and consistently retained less Pi than transgenic lines, especially when grown under low-phosphate conditions. This difference in Pi retention was attributable to the enhanced Pi uptake ability in the transgenic lines. The transgenic plants translocated more Pi into the plant cell due to the enhanced enzymatic activity of plasma membrane H+-ATPase and increased Pi/H+ co-transport, which improved Pi uptake. The differential expression of 'PHO regulon' genes further maintained intracellular Pi homeostasis. Furthermore, GB maintained a higher photosynthesis rate, thus increasing the production and translocation of sucrose via phloem loading to enhance plant response to low-phosphate stress. We conclude that GB mediates Pi uptake and translocation by regulating physiological and biochemical processes that promote adaptation to environmental changes in Pi availability. These processes eventually lead to better growth and development of the codA-transgenic lines. This finding will help to further elucidate the signaling mechanism of how GB perceives and transmits low-phosphate signals to alleviate Pi nutritional stress.

MicroR408 regulates defense response upon wounding in sweet potato

MiRNAs play diverse roles in plant development and defense responses by binding to their mRNA targets based on sequence complementarity. Here, we investigated a wound-related miR408 and its target genes in sweet potato (Ipomoea batatas) by small RNA deep sequencing and transcriptome analysis. The expression patterns of miR408 and the miR408 precursor were significantly repressed by wounding and jasmonate (JA). In contrast, expression of the putative target genes IbKCS (3-ketoacyl-CoA synthase 4), IbPCL (plantacyanin), and IbGAUT (galacturonosyltransferase 7-like) of miR408 was increased following wounding, whereas only IbKCS was increased after JA treatment. Target cleavage site mapping and Agrobacterium-mediated transient assay demonstrated that IbKCS, IbPCL, and IbGAUT were the targets of miR408. The expression of miR408 target genes was repressed in transgenic sweet potatoes overexpressing miR408. These data indicated a relationship between miR408 and its target genes. Notably, miR408-overexpressing plants showed a semi-dwarf phenotype and attenuated resistance to insect feeding, while transgenic plants overexpressing IbKCS exhibited more insect resistance than plants overexpressing only the empty vector. Collectively, sweet potato reduces the abundance of miR408 upon wounding to elevate the expression of IbKCS, IbPCL, and IbGAUT. The expression of IbKCS enhances the defense system against herbivore wounding.

The Altered Expression of microRNA408 Influences the Arabidopsis Response to Iron Deficiency

MicroRNAs contribute to the adaptation of plants to varying environmental conditions by affecting systemic mineral nutrient homeostasis. Copper and iron deficiencies antagonistically control the expression of Arabidopsis thaliana microRNA408 (miR408), which post-transcriptionally regulates laccase-like multicopper oxidase family members LAC3, LAC12, and LAC13. In this work, we used miR408 T-DNA insertion mutants (408-KO1 and 408-KO2) and a previously characterized transgenic line overexpressing miR408 (35S:408-14) to explore how miR408 influences copper- and iron-dependent metabolism. We observed that the altered expression of miR408 diminished plant performance and the activation of the iron-regulated genes under iron-deficient conditions. Consistently with the low expression of the miR408-target laccases, we showed that the vascular bundle lignification of the 35S:408-14 plants diminished. The decrease in the phenoloxidase and ferroxidase activities exhibited by wild-type plants under iron deficiency did not occur in the 408-KO1 plants, probably due to the higher expression of laccases. Finally, we observed that the hydrogen peroxide levels under iron starvation were altered in both the 408-KO1 and 35S:408-14 lines. Taken together, these results suggest that Arabidopsis plants with modified miR408 levels undergo multiple deregulations under iron-deficient conditions.

T7 exo-mediated FRET-breaking combined with DSN–RNAse–TdT for the detection of microRNA with ultrahigh signal-amplification

A DSN–RNAse–TdT–T7 exo probing system allows the detection of miRNA 21 with very high sensitivity (LOD = 2.57 fM) and selectivity—the result of (i) avoiding the false-positive signal from miRNA reacting with TdT polymerase and (ii) signal amplification occurring through a FRET-breaking mechanism involving T7 exo.

Chloroplast-to-Nucleus Signaling Regulates MicroRNA Biogenesis in Arabidopsis

As integral regulators in plant development and stress response, microRNAs (miRNAs) themselves need to be tightly regulated. Here, we show that tocopherols (vitamin E), lipid-soluble antioxidants synthesized from tyrosine in chloroplasts, positively regulate the biogenesis of miRNAs. Tocopherols are required for the accumulation of 3'-phosphoadenosine 5'-phosphate (PAP), a retrograde inhibitor of the nuclear exoribonucleases (XRN), which may protect primary miRNAs from being degraded and promote mature miRNA production. Such regulation is involved in heat-induced accumulation of miR398 and plant acquisition of heat tolerance. Our study reveals a chloroplast-to-nucleus signaling mechanism that favors miRNA biogenesis under heat and possibly other environmental perturbations.

Novel microRNAs are associated with population divergence in transcriptional response to thermal stress in an intertidal copepod

The role of gene expression in adaptation to differing thermal environments has been assayed extensively. Yet, in most natural systems, analyses of gene expression reveal only one level of the complexity of regulatory machineries. MicroRNAs (miRNAs) are small noncoding RNAs which are key components of many gene regulatory networks, and they play important roles in a variety of cellular pathways by modulating post-transcriptional quantities of mRNA available for protein synthesis. The characterization of miRNA loci and their regulatory dynamics in nonmodel systems are still largely understudied. In this study, we examine the role of miRNAs in response to high thermal stress in the intertidal copepod Tigriopus californicus. Allopatric populations of this species show varying levels of local adaptation with respect to thermal regimes, and previous studies showed divergence in gene expression between populations from very different thermal environments. We examined the transcriptional response to temperature stress in two populations separated by only 8 km by utilizing RNA-seq to quantify both mRNA and miRNA levels. Using the currently available genome sequence, we first describe the repertoire of miRNAs in T. californicus and assess the degree to which transcriptional response to temperature stress is governed by miRNA activity. The two populations showed large differences in the number of genes involved, the magnitude of change in commonly used genes and in the number of miRNAs involved in transcriptional modulation during stress. Our results suggest that an increased level of regulatory network complexity may underlie improved survivorship under thermal stress in one of the populations.

Regulation through MicroRNAs in Response to Low-Energy N+ Ion Irradiation in Oryza sativa

MicroRNAs (miRNAs) are a non-coding regulatory RNAs that play significant roles in plant growth and development, especially in the stress response. Low-energy ion radiation, a type of environmental stress, can cause multiple biological effects. To understand the roles of miRNAs in response to low-energy N⁺ ion radiation in Oryza sativa, high-throughput sequencing of small RNAs was carried out to detect the expression of miRNAs in the shoots of the rice after 2 × 10¹⁷ N+/cm² irradiation. The differentially expressed 28 known miRNAs were identified, 17 of these identified miRNAs were validated by real-time quantitative fluorescent PCR (q-PCR), including 9 up-regulated miRNAs (miR1320-3p, miR1320-5p, miR156b-3p, miR156c-5p, miR156c-3p/g-3p, miR1561-5p, miR398b and miR6250) and 8 down-regulated miRNAs (miR156a/e/i, miR156k, miR160f-5p, miR166j-5p, miR1846e and miR399d). In addition, 45 novel radiationresponsive miRNAs were predicted, and 8 of them were verified by q-PCR. The target genes of radiation-responsive miRNAs were predicted and gene function enrichment analysis was performed with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of 9 targets of 4 known miRNA families (miR156, miR399, miR1320 and miR398) and 2 targets of 2 novel miRNAs were quantified by q-PCR, and a strong negative regulation relation between miRNAs and their targets were observed. Those targets including SQUAMOSA promoterbinding-like protein (SPL) genes, copper/zinc superoxide dismutase (Cu/Zn-SOD), copper chaperone for SOD (CCS1) and electron transporter/ heat-shock protein binding protein (HSP), which are involved in growth and defense against various stresses, especially associated with reactive oxygen species (ROS) scavenging. This work provides important information for understanding the ROS generation and elimination mechanisms closely related to miRNAs in rice seedlings after low-energy N⁺ radiation exposure.

It's Hard to Avoid Avoidance: Uncoupling the Evolutionary Connection between Plant Growth, Productivity and Stress "Tolerance"

In the last 100 years, agricultural developments have favoured selection for highly productive crops, a fact that has been commonly associated with loss of key traits for environmental stress tolerance. We argue here that this is not exactly the case. We reason that high yield under near optimal environments came along with hypersensitization of plant stress perception and consequently early activation of stress avoidance mechanisms, such as slow growth, which were originally needed for survival over long evolutionary time periods. Therefore, mechanisms employed by plants to cope with a stressful environment during evolution were overwhelmingly geared to avoid detrimental effects so as to ensure survival and that plant stress "tolerance" is fundamentally and evolutionarily based on "avoidance" of injury and death which may be referred to as evolutionary avoidance (EVOL-Avoidance). As a consequence, slow growth results from being exposed to stress because genes and genetic programs to adjust growth rates to external circumstances have evolved as a survival but not productivity strategy that has allowed extant plants to avoid extinction. To improve productivity under moderate stressful conditions, the evolution-oriented plant stress response circuits must be changed from a survival mode to a continued productivity mode or to avoid the evolutionary avoidance response, as it were. This may be referred to as Agricultural (AGRI-Avoidance). Clearly, highly productive crops have kept the slow, reduced growth response to stress that they evolved to ensure survival. Breeding programs and genetic engineering have not succeeded to genetically remove these responses because they are polygenic and redundantly programmed. From the beginning of modern plant breeding, we have not fully appreciated that our crop plants react overly-cautiously to stress conditions. They over-reduce growth to be able to survive stresses for a period of time much longer than a cropping season. If we are able to remove this polygenic redundant survival safety net we may improve yield in moderately stressful environments, yet we will face the requirement to replace it with either an emergency slow or no growth (dormancy) response to extreme stress or use resource management to rescue crops under extreme stress (or both).

MicroRNAs and their targets in cucumber shoot apices in response to temperature and photoperiod

BACKGROUND

The cucumber is one of the most important vegetables worldwide and is used as a research model for study of phloem transport, sex determination and temperature-photoperiod physiology. The shoot apex is the most important plant tissue in which the cell fate and organ meristems have been determined. In this study, a series of whole-genome small RNA, degradome and transcriptome analyses were performed on cucumber shoot apical tissues treated with high vs. low temperature and long vs. short photoperiod.

RESULTS

A total of 164 known miRNAs derived from 68 families and 203 novel miRNAs from 182 families were identified. Their 4611 targets were predicted using psRobot and TargetFinder, amongst which 349 were validated by degradome sequencing. Fourteen targets of six miRNAs were differentially expressed between the treatments. A total of eight known and 16 novel miRNAs were affected by temperature and photoperiod. Functional annotations revealed that "Plant hormone signal transduction" pathway was significantly over-represented in the miRNA targets. The miR156/157/SBP-Boxes and novel-mir153/ethylene-responsive transcription factor/senescence-related protein/aminotransferase/acyl-CoA thioesterase are the two most credible miRNA/targets combinations modulating the plant's responsive processes to the temperature-photoperiod changes. Moreover, the newly evolved, cucumber-specific novel miRNA (novel-mir153) was found to target 2087 mRNAs by prediction and has 232 targets proven by degradome analysis, accounting for 45.26-58.88% of the total miRNA targets in this plant. This is the largest sum of genes targeted by a single miRNA to the best of our knowledge.

CONCLUSIONS

These results contribute to a better understanding of the miRNAs mediating plant adaptation to combinations of temperature and photoperiod and sheds light on the recent evolution of new miRNAs in cucumber.

An integrated analysis of mRNA and sRNA transcriptional profiles in tomato root: Insights on tomato wilt disease

Tomato wilt disease caused by Fusarium oxysporum f. sp. lycopersici (FOL) is a worldwide destructive disease of tomato. As exploring gene expression and function approaches constitute an initial point for investigating pathogen-host interaction, we performed RNA-seq and sRNA-seq analysis to investigate the transcriptome of tomato root under FOL infection. Differentially expressed (DE) protein-coding gene and miRNA gene profiles upon inoculation with FOL were presented at twenty-four hours post-inoculation in four treatments. A total of more than 182.6 million and 132.2 million high quality clean reads were obtained by RNA-seq and sRNA-seq, respectively. A large overlap was found in DE mRNAs between susceptible cultivar Moneymaker and resistant cultivar Motelle. Gene Ontology terms were mainly classified into catalytic activity, metabolic process and binding. Combining with qRT-PCR and Northern blot, we validated the transcriptional profile of five genes and five miRNAs conferred to FOL infection. Our work allowed comprehensive understanding of different transcriptional reaction of genes/miRNAs between the susceptible and resistant cultivars tomato to the FOL challenge, which could offer us with a future direction to generate models of mediated resistance responses.

Comparative profiling of microRNAs and their effects on abiotic stress in wild-type and dark green leaf color mutant plants of Anthurium andraeanum 'Sonate'

MicroRNAs (miRNAs) are a class of non-coding small RNAs that play important roles in the regulation of gene expression. Although plant miRNAs have been extensively studied in model systems, less is known in other plants with limited genome sequence data, including Anthurium andraeanum. To identify miRNAs and their target genes in A. andraeanum and study their responses to abiotic stresses, we conducted deep-sequencing of two small RNA (sRNA) libraries prepared from young leaves of wild-type (WT) and dark green (dg) leaf color mutant plants of A. andraeanum 'Sonate'. A total of 53 novel miRNAs were identified, 32 of which have been annotated to 18 miRNA families. 10 putative miRNAs were found to be differentially expressed in WT and dg, among which two miRNAs were significantly up-regulated and eight down-regulated in dg relative to WT. One differentially expressed miRNA, Aa-miR408, was dramatically up-regulated in dg. qRT-PCR analysis and heterologous expression of Aa-miR408 in Arabidopsis under different stress treatments suggest that Aa-miR408 is involved in abiotic stress responses in A. andraeanum. Our results provide a foundation for further dissecting the roles of miRNAs and their targets in regulating abiotic stress tolerance in A. andraeanum.

Ascophyllum nodosum extract mitigates salinity stress in Arabidopsis thaliana by modulating the expression of miRNA involved in stress tolerance and nutrient acquisition

Ascophyllum nodosum extract (ANE) contains bioactive compounds that improve the growth of Arabidopsis in experimentally-induced saline conditions; however, the molecular mechanisms through which ANE elicits tolerance to salinity remain largely unexplored. Micro RNAs (miRNAs) are key regulators of gene expression, playing crucial roles in plant growth, development, and stress tolerance. Next generation sequencing of miRNAs from leaves of control Arabidopsis and from plants subjected to three treatments (ANE, NaCl and ANE+NaCl) was used to identify ANE-responsive miRNA in the absence and presence of saline conditions. Differential gene expression analysis revealed that ANE had a strong effect on miRNAs expression in both conditions. In the presence of salinity, ANE tended to reduce the up-regulation or the down-regulation trend induced caused by NaCl in miRNAs such as ath-miR396a-5p, ath-miR399, ath-miR2111b and ath-miR827. To further uncover the effects of ANE, the expression of several target genes of a number of ANE-responsive miRNAs was analyzed by qPCR. NaCl, but not ANE, down-regulated miR396a-5p, which negatively regulated the expression of AtGRF7 leading to a higher expression of AtDREB2a and AtRD29 in the presence of ANE+NaCl, as compared to ANE alone. ANE+NaCl initially reduced and then enhanced the expression of ath-miR169g-5p, while the expression of the target genes AtNFYA1 and ATNFYA2, known to be involved in the salinity tolerance mechanism, was increased as compared to ANE or to NaCl treatments. ANE and ANE+NaCl modified the expression of ath-miR399, ath-miR827, ath-miR2111b, and their target genes AtUBC24, AtWAK2, AtSYG1 and At3g27150, suggesting a role of ANE in phosphate homeostasis. In vivo and in vitro experiments confirmed the improved growth of Arabidopsis in presence of ANE, in saline conditions and in phosphate-deprived medium, further substantiating the influence of ANE on a variety of essential physiological processes in Arabidopsis including salinity tolerance and phosphate uptake.

Elucidating the Small Regulatory RNA Repertoire of the Sea Anemone Anemonia viridis Based on Whole Genome and Small RNA Sequencing

Cnidarians harbor a variety of small regulatory RNAs that include microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs), but detailed information is limited. Here, we report the identification and expression of novel miRNAs and putative piRNAs, as well as their genomic loci, in the symbiotic sea anemone Anemonia viridis. We generated a draft assembly of the A. viridis genome with putative size of 313 Mb that appeared to be composed of about 36% repeats, including known transposable elements. We detected approximately equal fractions of DNA transposons and retrotransposons. Deep sequencing of small RNA libraries constructed from A. viridis adults sampled at a natural CO₂ gradient off Vulcano Island, Italy, identified 70 distinct miRNAs. Eight were homologous to previously reported miRNAs in cnidarians, whereas 62 appeared novel. Nine miRNAs were recognized as differentially expressed along the natural seawater pH gradient. We found a highly abundant and diverse population of piRNAs, with a substantial fraction showing ping–pong signatures. We identified nearly 22% putative piRNAs potentially targeting transposable elements within the A. viridis genome. The A. viridis genome appeared similar in size to that of other hexacorals with a very high divergence of transposable elements resembling that of the sea anemone genus Exaiptasia. The genome encodes and expresses a high number of small regulatory RNAs, which include novel miRNAs and piRNAs. Differentially expressed small RNAs along the seawater pH gradient indicated regulatory gene responses to environmental stressors.

Identification and functional characterization of intermediate-size non-coding RNAs in maize

BACKGROUND

The majority of eukaryote genomes can be actively transcribed into non-coding RNAs (ncRNAs), which are functionally important in development and evolution. In the study of maize, an important crop for both humans and animals, aside from microRNAs and long non-coding RNAs, few studies have been conducted on intermediate-size ncRNAs.

RESULTS

We constructed a homogenized cDNA library of 50-500 nt RNAs in the maize inbred line Chang 7-2. Sequencing revealed 169 ncRNAs, which contained 58 known and 111 novel ncRNAs (including 70 snoRNAs, 27 snRNAs, 13 unclassified ncRNAs and one tRNA). Forty of the novel ncRNAs were specific to the Panicoideae, and 24% of them are located on sense-strand of the 5' or 3' terminus of protein coding genes on chromosome. Target site analysis found that 22 snoRNAs can guide to 38 2'-O-methylation and pseudouridylation modification sites of ribosomal RNAs and small nuclear RNAs. Expression analysis showed that 43 ncRNAs exhibited significantly altered expression in different tissues or developmental stages of maize seedlings, eight ncRNAs had tissue-specific expression and five ncRNAs were strictly accumulated in the early stage of leaf development. Further analysis showed that 3 of the 5 stage-specific ncRNAs (Zm-3, Zm-18, and Zm-73) can be highly induced under drought and salt stress, while one snoRNA Zm-8 can be repressed under PEG-simulated drought condition.

CONCLUSIONS

We provided a genome-wide identification and functional analysis of ncRNAs with a size range of 50-500 nt in maize. 111 novel ncRNAs were cloned and 40 ncRNAs were determined to be specific to Panicoideae. 43 ncRNAs changed significantly during maize development, three ncRNAs can be strongly induced under drought and salt stress, suggesting their roles in maize stress response. This work set a foundation for further study of intermediate-size ncRNAs in maize.

Differentially Expressed MicroRNAs Link Cellular Physiology to Phenotypic Changes in Rice Under Stress Conditions

Plant microRNAs (miRNAs) and their target genes have important functional roles in nutrition deficiency and stress response. However, the underlying mechanisms relating relative expression of miRNAs and target mRNAs to morphological adjustments are not well defined. By combining miRNA expression profiles, corresponding target genes and transcription factors that bind to computationally identified over-represented cis-regulatory elements (CREs) common in miRNAs and target gene promoters, we implement a strategy that identifies a set of differentially expressed regulatory interactions which, in turn, relate underlying cellular mechanisms to some of the phenotypic changes observed. Integration of experimentally reported individual interactions with identified regulatory interactions explains how (i) during mineral deficiency osa-miR167 inhibits shoot growth but activates adventitious root growth by influencing free auxin content; (ii) during sulfur deficiency osa-miR394 is involved in adventitious root growth inhibition, sulfur and iron homeostasis, and auxin-mediated regulation of sulfur homeostasis; (iii) osa-miR399 contributes to cross-talk between cytokinin and phosphorus deficiency signaling; and (iv) a feed-forward loop involving the osa-miR166, trihelix and HD-ZIP III transcription factors may regulate leaf senescence during drought. This strategy not only identifies various regulatory interactions connecting phenotypic changes with cellular or molecular events triggered by stress, but also provides a framework to deepen our understanding of stress cellular physiology.

sRNAome and transcriptome analysis provide insight into chilling response of cowpea pods

Cowpea is an important horticultural crop in tropical and subtropical areas of Asia, Africa, and Latin America, as well as parts of southern Europe and Central and South America. Chilling injury is a common physiological hazard of cowpea in cold chain logistics which reduce the cowpea pods nutritional quality and product value. However, the molecular mechanism involved in chilling injury remains unclear in cowpea pods. RNA-Seq and sRNA-Seq technologies were employed to decipher the miRNAs and mRNAs expression profiles and their regulatory networks in cowpea pods involved in chilling stress. Differentially expressed miRNAs and mRNA profiles were obtained based on cluster analysis, miRNAs and target genes were found to show coherent relationships in the regulatory networks of chilling injury. Furthermore, we found that numerous miRNAs and nat-siRNAs' targets were predicted to be key enzymes involved in the redox reactions such as POD, CAT, AO and LOX, energy metabolism such as ATPase, FAD and NAD related enzymes and different transcription factors such as WRKY, bHLH, MYB, ERF and NAC which play important roles in chilling injury.

Nitrogen Supply and Leaf Age Affect the Expression of TaGS1 or TaGS2 Driven by a Constitutive Promoter in Transgenic Tobacco

Glutamine synthetase (GS) plays a key role in nitrogen metabolism. Here, two types of tobacco transformants, overexpressing Triticum aestivum GS1 (TaGS1) or GS2 (TaGS2), were analysed. Four independent transformed lines, GS1-TR1, GS1-TR2, GS2-TR1 and GS2-TR2, were used for the nitrogen treatment. Under nitrogen-sufficient conditions, the leaves of GS2-TR showed high accumulation of the TaGS2 transcript, while those of GS1-TR showed a low TaGS1 transcript levels. However, compared with nitrogen-sufficient conditions, the TaGS1 transcript level increased in the leaves under nitrogen starvation, but the TaGS2 transcript level decreased. In addition, the TaGS1 and TaGS2 transcript levels were highest in the middle leaves under nitrogen-sufficient and starvation conditions. These results show that nitrogen supply and leaf age regulate TaGS expression, even when they are driven by a super-promoter. Additionally, in regard to nitrogen metabolism level, the lower leaves of the GS1-TR exhibited lower NH₄⁺ and higher amino acid contents, while the upper leaves exhibited higher amino acid, soluble protein and chlorophyll contents. The leaves of the GS2-TR exhibited lower NH₄⁺ but higher amino acid, soluble protein and chlorophyll contents. Given the role that GS isoforms play in nitrogen metabolism, these data suggest that TaGS1 overexpression may improve nitrogen transport, and that TaGS2 overexpression may improve nitrogen assimilation under nitrogen stress.

Identification of exosome-like nanoparticle-derived microRNAs from 11 edible fruits and vegetables

Edible plant-derived exosome-like nanoparticles (EPDELNs) are novel naturally occurring plant ultrastructures that are structurally similar to exosomes. Many EPDELNs have anti-inflammatory properties. MicroRNAs (miRNAs) play a critical role in mediating physiological and pathological processes in animals and plants. Although miRNAs can be selectively encapsulated in extracellular vesicles, little is known about their expression and function in EPDELNs. In this study, we isolated nanovesicles from 11 edible fruits and vegetables and subjected the corresponding EPDELN small RNA libraries to Illumina sequencing. We identified a total of 418 miRNAs-32 to 127 per species-from the 11 EPDELN samples. Target prediction and functional analyses revealed that highly expressed miRNAs were closely associated with the inflammatory response and cancer-related pathways. The 418 miRNAs could be divided into three classes according to their EPDELN distributions: 26 "frequent" miRNAs (FMs), 39 "moderately present" miRNAs (MPMs), and 353 "rare" miRNAs (RMs). FMs were represented by fewer miRNA species than RMs but had a significantly higher cumulative expression level. Taken together, our in vitro results indicate that miRNAs in EPDELNs have the potential to regulate human mRNA.

In silico identification and computational characterization of endogenous small interfering RNAs from diverse grapevine tissues and stages

Small interfering RNAs (siRNAs) are effectors of regulatory pathways underlying plant development, metabolism, and stress- and nutrient-signaling regulatory networks. The endogenous siRNAs are generally not conserved between plants; consequently, it is necessary and important to identify and characterize siRNAs from various plants. To address the nature and functions of siRNAs, and understand the biological roles of the huge siRNA population in grapevine (Vitis vinifera L.). The high-throughput sequencing technology was used to identify a large set of putative endogenous siRNAs from six grapevine tissues/organs. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to classify the target genes of siRNA. In total, 520,519 candidate siRNAs were identified and their expression profiles exhibited typical temporal characters during grapevine development. In addition, we identified two grapevine trans-acting siRNA (TAS) gene homologs (VvTAS3 and VvTAS4) and the derived trans-acting siRNAs (tasiRNAs) that could target grapevine auxin response factor (ARF) and myeloblastosis (MYB) genes. Furthermore, the GO and KEGG analysis of target genes showed that most of them covered a broad range of functional categories, especially involving in disease-resistance process. The large-scale and completely genome-wide level identification and characterization of grapevine endogenous siRNAs from the diverse tissues by high throughput technology revealed the nature and functions of siRNAs in grapevine.

Perspectives on Structural, Physiological, Cellular, and Molecular Responses to Desiccation in Resurrection Plants

Resurrection plants possess a unique ability to counteract desiccation stress. Desiccation tolerance (DT) is a very complex multigenic and multifactorial process comprising a combination of physiological, morphological, cellular, genomic, transcriptomic, proteomic, and metabolic processes. Modification in the sugar composition of the hemicellulosic fraction of the cell wall is detected during dehydration. An important change is a decrease of glucose in the hemicellulosic fraction during dehydration that can reflect a modification of the xyloglucan structure. The expansins might also be involved in cell wall flexibility during drying and disrupt hydrogen bonds between polymers during rehydration of the cell wall. Cleavages by xyloglucan-modifying enzymes release the tightly bound xyloglucan-cellulose network, thus increasing cell wall flexibility required for cell wall folding upon desiccation. Changes in hydroxyproline-rich glycoproteins (HRGPs) such as arabinogalactan proteins (AGPs) are also observed during desiccation and rehydration processes. It has also been observed that significant alterations in the process of photosynthesis and photosystem (PS) II activity along with changes in the antioxidant enzyme system also increased the cell wall and membrane fluidity resulting in DT. Similarly, recent data show a major role of ABA, LEA proteins, and small regulatory RNA in regulating DT responses. Current progress in "-omic" technologies has enabled quantitative monitoring of the plethora of biological molecules in a high throughput routine, making it possible to compare their levels between desiccation-sensitive and DT species. In this review, we present a comprehensive overview of structural, physiological, cellular, molecular, and global responses involved in desiccation tolerance.

Integration of proteomic and transcriptomic profiles reveals multiple levels of genetic regulation of salt tolerance in cotton

BACKGROUND

Salinity is a major abiotic stress that limits upland cotton growth and reduces fibre production worldwide. To reveal genetic regulation via transcript and protein levels after salt stress, we comprehensively analysed the global changes in mRNA, miRNA, and protein profiles in response to salt stress in two contrasting salt-tolerant cotton genotypes.

RESULTS

In the current study, proteomic and mRNA-seq data were combined to reveal that some genes are differentially expressed at both the proteomic and mRNA levels. However, we observed no significant change in mRNA corresponding to most of the strongly differentially abundant proteins. This finding may have resulted from global changes in alternative splicing events and miRNA levels under salt stress conditions. Evidence was provided indicating that several salt stress-responsive proteins can alter miRNAs and modulate alternative splicing events in upland cotton. The results of the stringent screening of the mRNA-seq and proteomic data between the salt-tolerant and salt-sensitive genotypes identified 63 and 85 candidate genes/proteins related to salt tolerance after 4 and 24 h of salt stress, respectively, between the tolerant and sensitive genotype. Finally, we predicted an interaction network comprising 158 genes/proteins and then discovered that two main clusters in the network were composed of ATP synthase (CotAD_74681) and cytochrome oxidase (CotAD_46197) in mitochondria. The results revealed that mitochondria, as important organelles involved in energy metabolism, play an essential role in the synthesis of resistance proteins during the process of salt exposure.

CONCLUSION

We provided a plausible schematic for the systematic salt tolerance model; this schematic reveals multiple levels of gene regulation in response to salt stress in cotton and provides a list of salt tolerance-related genes/proteins. The information here will facilitate candidate gene discovery and molecular marker development for salt tolerance breeding in cotton.

Identification of microRNAs and long non-coding RNAs involved in fatty acid biosynthesis in tree peony seeds

MicroRNAs (miRNAs) and long noncoding RNAs (lncRNAs) act as important molecular regulators in a wide range of biological processes during plant development and seed formation, including oil production. Tree peony seeds contain >90% unsaturated fatty acids (UFAs) and high proportions of α-linolenic acid (ALA, > 40%). To dissect the non-coding RNAs (ncRNAs) pathway involved in fatty acids synthesis in tree peony seeds, we construct six small RNA libraries and six transcriptome libraries from developing seeds of two cultivars (J and S) containing different content of fatty acid compositions. After deep sequencing the RNA libraries, the ncRNA expression profiles of tree peony seeds in two cultivars were systematically and comparatively analyzed. A total of 318 known and 153 new miRNAs and 22,430 lncRNAs were identified, among which 106 conserved and 9 novel miRNAs and 2785 lncRNAs were differentially expressed between the two cultivars. In addition, potential target genes of the microRNA and lncRNAs were also predicted and annotated. Among them, 9 miRNAs and 39 lncRNAs were predicted to target lipid related genes. Results showed that all of miR414, miR156b, miR2673b, miR7826, novel-m0027-5p, TR24651|c0_g1, TR24544|c0_g15, and TR27305|c0_g1 were up-regulated and expressed at a higher level in high-ALA cultivar J when compared to low-ALA cultivar S, suggesting that these ncRNAs and target genes are possibly involved in different fatty acid synthesis and lipid metabolism through post-transcriptional regulation. These results provide a better understanding of the roles of ncRNAs during fatty acid biosynthesis and metabolism in tree peony seeds.

Identification of miRNAs linked with the drought response of tef [Eragrostis tef (Zucc.) Trotter]

Tef [Eragrostis tef (Zucc.) Trotter], a staple food crop in the Horn of Africa and particularly in Ethiopia, has several beneficial agronomical and nutritional properties, including waterlogging and drought tolerance. In this study, we performed microRNA profiling of tef using the Illumina HiSeq 2500 platform, analyzing both shoots and roots of two tef genotypes, one drought-tolerant (Tsedey) and one drought-susceptible (Alba). We obtained more than 10 million filtered reads for each of the 24 sequenced small cDNA libraries. Reads mapping to known miRNAs were more abundant in the root than shoot tissues. Thirteen and 35 miRNAs were significantly modulated in response to drought, in Alba and Tsedey roots, respectively. One miRNA was upregulated under drought conditions in both genotypes. In shoots, nine miRNAs were modulated in common between the two genotypes and all showed similar trends of expression. One-hundred and forty-seven new miRNA mature sequences were identified in silico, 22 of these were detected in all relevant samples and seven were differentially regulated when comparing drought with normal watering. Putative targets of the miRNA regulated under drought in root and shoot tissues were predicted. Among the targets were transcription factors such as CCAAT-HAP2, MADS and NAC. Verification with qRT-PCR revealed that five of six potential targets showed a pattern of expression that was consistent with the correspondent miRNA amount measured by RNA-Seq. In general, candidate miRNAs involved in the post-transcriptional regulation of the tef response to drought could be included in next-generation breeding programs.

Recent Advances in the Transcriptional Regulation of Secondary Cell Wall Biosynthesis in the Woody Plants

Plant cell walls provide structural support for growth and serve as a barrier for pathogen attack. Plant cell walls are also a source of renewable biomass for conversion to biofuels and bioproducts. Understanding plant cell wall biosynthesis and its regulation is of critical importance for the genetic modification of plant feedstocks for cost-effective biofuels and bioproducts conversion and production. Great progress has been made in identifying enzymes involved in plant cell wall biosynthesis, and in Arabidopsis it is generally recognized that the regulation of genes encoding these enzymes is under a transcriptional regulatory network with coherent feedforward and feedback loops. However, less is known about the transcriptional regulation of plant secondary cell wall (SCW) biosynthesis in woody species despite of its high relevance to biofuels and bioproducts conversion and production. In this article, we synthesize recent progress on the transcriptional regulation of SCW biosynthesis in Arabidopsis and contrast to what is known in woody species. Furthermore, we evaluate progress in related emerging regulatory machineries targeting transcription factors in this complex regulatory network of SCW biosynthesis.

Evolutionary Origin, Gradual Accumulation and Functional Divergence of Heat Shock Factor Gene Family with Plant Evolution

Plants, as sessile organisms, evolved a complex and functionally diverse heat shock factor (HSF) gene family to cope with various environmental stresses. However, the limited evolution studies of the HSF gene family have hindered our understanding of environmental adaptations in plants. In this study, a comprehensive evolution analysis on the HSF gene family was performed in 51 representative plant species. Our results demonstrated that the HSFB group which lacks a typical AHA activation domain, was the most ancient, and is under stronger purifying selection pressure in the subsequent evolutionary processes. While, dramatic gene expansion and functional divergence occurred at evolution timescales corresponding to plant land inhabit, which contribute to the emergence and diversification of the HSFA and HSFC groups in land plants. During the plant evolution, the ancestral functions of HSFs were maintained by strong purifying pressure that acted on the DNA binding domain, while the variable oligomerization domain and motif organization of HSFs underwent functional divergence and generated novel subfamilies. At the same time, variations were further accumulated with plant evolution, and this resulted in remarkable functional diversification among higher plant lineages, including distinct HSF numbers and selection pressures of several HSF subfamilies between monocots and eudicots, highlighting the fundamental differences in different plant lineages in response to environmental stresses. Taken together, our study provides novel insights into the evolutionary origin, pattern and selection pressure of plant HSFs and delineates critical clues that aid our understanding of the adaptation processes of plants to terrestrial environments.

Chinese Wheat Mosaic Virus-Induced Gene Silencing in Monocots and Dicots at Low Temperature

Virus-induced gene silencing (VIGS) is an important tool for functional genomics studies in plants. With this method, it is possible to target most endogenous genes and downregulate the messenger RNA (mRNA) in a sequence-specific manner. Chinese wheat mosaic virus (CWMV) has a bipartite, single-strand positive RNA genome, and can infect both wheat and Nicotiana benthamiana, and the optimal temperature for systemic infection in plants is 17°C. To assess the potential of the virus as a vector for gene silencing at low temperature, a fragment of the N. benthamiana or wheat phytoene desaturase (PDS) gene was expressed from a modified CWMV RNA2 clone and the resulting photo bleaching in infected plants was used as a reporter for silencing. Downregulation of PDS mRNA was also measured by quantitative reverse-transcriptase polymerase chain reaction (RT-qPCR). In experiments using fragments of PDS ranging from 500 to 1500 nucleotides, insert length influenced the stability and the efficiency of VIGS. The CWMV induced silencing system was also used to suppress miR165/166 and miR3134a through expression of miRNA target mimics. The relative expression levels of mature miR165/166 and miR3134a decreased whereas the transcript levels of their target genes increased. Interestingly, we also found the CWMV-induced silencing system was more efficient compare with the vector based on Barley stripe mosaic virus (BSMV) or Foxtail mosaic virus (FoMV) in wheat or the vector based on TRV in N. benthamiana at 17°C. In summary, the CWMV vector is effective in silencing endogenous genes and miRNAs at 17°C, thereby providing a powerful tool for gene function analysis in both N. benthamiana and wheat at low temperature.

Regulatory Gene Networks in Drought Stress Responses and Resistance in Plants

Plant responses to drought stress have been analyzed extensively to reveal complex regulatory gene networks, including the detection of water deficit signals, as well as the physiological, cellular, and molecular responses. Plants recognize water deficit conditions at their roots and transmit this signal to their shoots to synthesize abscisic acid (ABA) in their leaves. ABA is a key phytohormone that regulates physiological and molecular responses to drought stress, such as stomatal closure, gene expression, and the accumulation of osmoprotectants and stress proteins. ABA transporters function as the first step for propagating synthesized ABA. To prevent water loss, ABA influx in guard cells is detected by several protein kinases, such as SnRK2s and MAPKs that regulate stomatal closure. ABA mediates a wide variety of gene expression machineries with stress-responsive transcription factors, including DREBs and AREBs, to acquire drought stress resistance in whole tissues. In this chapter, we summarize recent advances in drought stress signaling, focusing on gene networks in cellular and intercellular stress responses and drought resistance.

Genome-Wide Identification and Characterization of Long Non-Coding RNA in Wheat Roots in Response to Ca2+ Channel Blocker

It remains unclear whether plant lncRNAs are responsive to Ca²⁺-channel blocking. When using the Ca²⁺-channel blocker, LaCl₃, to treat germinated wheat seeds for 24 h, we found that both root length and mitosis were inhibited in the LaCl₃-treated groups. The effect of the Ca²⁺-channel blocker was verified in three ways: a [Ca²⁺]_cyt decrease detected using Fluo-3/AM staining, a decrease in the Ca content measured using inductively coupled plasma mass spectrometry, and an inhibition of Ca²⁺ influx detected using Non-invasive Micro-test Technology. Genome-wide high throughput RNA-seq and bioinformatical methods were used to identify lncRNAs, and found 177 differentially expressed lncRNAs that might be in responsive to Ca²⁺-channel blocking. Among these, 108 were up-regulated and 69 were down-regulated. The validity of identified lncRNAs data from RNA-seq was verified using qPCR. GO and KEGG analysis indicated that a number of lncRNAs might be involved in diverse biological processes upon Ca²⁺-channel blocking. Further GO analysis showed that 23 lncRNAs might play roles as transcription factor (TF); Moreover, eight lncRNAs might participate in cell cycle regulation, and their relative expressions were detected using qPCR. This study also provides diverse data on wheat lncRNAs that can deepen our understanding of the function and regulatory mechanism of Ca²⁺-channel blocking in plants.

FRG3, a Target of slmiR482e-3p, Provides Resistance against the Fungal Pathogen Fusarium oxysporum in Tomato

The vast majority of plant disease resistance (R) genes encode nucleotide binding site-leucine-rich repeat (NBS-LRR) proteins, which specifically determine the plant immune response and have been demonstrated to be targets of several microRNA (miRNA) families. The fungus Fusarium oxysporum f. sp. lycopersici (FOL) causes vascular wilt disease in tomato worldwide. Here, we explored a possible role for FGR3 in tomato defense against FOL. FRG3 is a predicted NBS-LRR like gene that is targeted by slmiR482e-3p, a member of slmiR482 miRNA family. Northern blot data demonstrated that all seven members of the slmiR482 family were regulated in diverse ways after infection by FOL. The ability of FRG3 to be regulated by slmiR482e-3p was confirmed at the transcript level by co-expression studies in Nicotiana benthamiana. A virus-induced gene silencing (VIGS) approach revealed that FRG3 confers resistance to the Motelle tomato cultivar. Taken together, our study has identified a novel R gene, FRG3, which is targeted by slmiR482e-3p at the transcript level, and is necessary for resistance to tomato wilt disease in planta.

Cold-responsive miRNAs and their target genes in the wild eggplant species Solanum aculeatissimum

Background

Low temperature is an important abiotic stress in plant growth and development, especially for thermophilic plants. Eggplants are thermophilic vegetables, although the molecular mechanism of their response to cold stress remains to be elucidated. MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play an essential role during plant development and stress responses. Although the role of many plant miRNAs in facilitating chilling tolerance has been verified, little is known about the mechanisms of eggplant chilling tolerance.

Results

Here, we used high-throughput sequencing to extract the miRNA and target genes expression profiles of Solanum aculeatissimum (S. aculeatissimum) under low temperature stress at different time periods(0 h, 2 h, 6 h, 12 h, 24 h). Differentially regulated miRNAs and their target genes were analyzed by comparing the small RNA (sRNA) and miRBase 20.0 databases using BLAST or BOWTIE, respectively. Fifty-six down-regulated miRNAs and 28 up-regulated miRNAs corresponding to 220 up-regulated mRNAs and 94 down-regulated mRNAs, respectively, were identified in S. aculeatissimum. Nine significant differentially expressed miRNAs and twelve mRNAs were identified by quantitative Real-time PCR and association analysis, and analyzed for their GO function enrichment and KEGG pathway association.

Conclusions

In summary, numerous conserved and novel miRNAs involved in the chilling response were identified using high-throughput sequencing, which provides a theoretical basis for the further study of low temperature stress-related miRNAs and the regulation of cold-tolerance mechanisms of eggplant at the miRNA level.

Electronic supplementary material

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

The Importance of ncRNAs as Epigenetic Mechanisms in Phenotypic Variation and Organic Evolution

Neo-Darwinian explanations of organic evolution have settled on mutation as the principal factor in producing evolutionary novelty. Mechanistic characterizations have been also biased by the classic dogma of molecular biology, where only proteins regulate gene expression. This together with the rearrangement of genetic information, in terms of genes and chromosomes, was considered the cornerstone of evolution at the level of natural populations. This predominant view excluded both alternative explanations and phenomenologies that did not fit its paradigm. With the discovery of non-coding RNAs (ncRNAs) and their role in the control of genetic expression, new mechanisms arose providing heuristic power to complementary explanations to evolutionary processes overwhelmed by mainstream genocentric views. Viruses, epimutation, paramutation, splicing, and RNA editing have been revealed as paramount functions in genetic variations, phenotypic plasticity, and diversity. This article discusses how current epigenetic advances on ncRNAs have changed the vision of the mechanisms that generate variation, how organism-environment interaction can no longer be underestimated as a driver of organic evolution, and how it is now part of the transgenerational inheritance and evolution of species.