MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes

Integrated mRNA and miRNA transcriptome analysis reveals a regulatory network for tuber expansion in Chinese yam (Dioscorea opposita)

BACKGROUND

Yam tuber is a storage organ, derived from the modified stem. Tuber expansion is a complex process, and depends on the expressions of genes that can be influenced by environmental and endogenous factors. However, little is known about the regulatory mechanism of tuber expansion. In order to identify the genes and miRNAs involved in tuber expansion, we examined the mRNAs and small RNAs in Dioscorea opposita (Chinese yam) cv. Guihuai 16 tuber during its initiation and expansion stages.

RESULTS

A total of 14,238 differentially expressed genes in yam tuber at its expansion stage were identified by using RNA sequencing technology. Among them, 5723 genes were up-regulated, and 8515 genes were down-regulated. Functional analysis revealed the coordination of tuber plant involved in processes of cell events, metabolism, biosynthesis, and signal transduction pathways at transcriptional level, suggesting that these differentially expressed genes are somehow involved in response to tuber expansion, including CDPK, CaM, CDL, SAUR, DELLA, SuSy, and expansin. In addition, 541 transcription factor genes showed differential expression during the expansion stage at transcriptional level. MADS, bHLH, and GRAS were involved in cell differentiation, division, and expansion, which may relate to tuber expansion. Noteworthy, data analysis revealed that 22 known tuber miRNAs belong to 10 miRNA families, and 50 novel miRNAs were identified. The integrated analysis of miRNA-mRNA showed that 4 known miRNAs and 11 genes formed 14 miRNA-target mRNA pairs were co-expressed in expansion stage. miRNA160, miRNA396, miRNA535 and miRNA5021 may be involved in complex network to regulate cell division and differentiation in yam during its expansion stage.

CONCLUSION

The mRNA and miRNA datasets presented here identified a subset of candidate genes and miRNAs that are putatively associated with tuber expansion in yam, a hypothetical model of genetic regulatory network associated with tuber expansion in yam was put forward, which may provide a foundation for molecular regulatory mechanism researching on tuber expansion in Dioscorea species.

Integrated mRNA and miRNA transcriptome analysis reveals a regulatory network for tuber expansion in Chinese yam (Dioscorea opposita)

Background

Yam tuber is a storage organ, derived from the modified stem. Tuber expansion is a complex process, and depends on the expressions of genes that can be influenced by environmental and endogenous factors. However, little is known about the regulatory mechanism of tuber expansion. In order to identify the genes and miRNAs involved in tuber expansion, we examined the mRNAs and small RNAs in Dioscorea opposita (Chinese yam) cv. Guihuai 16 tuber during its initiation and expansion stages.

Results

A total of 14,238 differentially expressed genes in yam tuber at its expansion stage were identified by using RNA sequencing technology. Among them, 5723 genes were up-regulated, and 8515 genes were down-regulated. Functional analysis revealed the coordination of tuber plant involved in processes of cell events, metabolism, biosynthesis, and signal transduction pathways at transcriptional level, suggesting that these differentially expressed genes are somehow involved in response to tuber expansion, including CDPK, CaM, CDL, SAUR, DELLA, SuSy, and expansin. In addition, 541 transcription factor genes showed differential expression during the expansion stage at transcriptional level. MADS, bHLH, and GRAS were involved in cell differentiation, division, and expansion, which may relate to tuber expansion. Noteworthy, data analysis revealed that 22 known tuber miRNAs belong to 10 miRNA families, and 50 novel miRNAs were identified. The integrated analysis of miRNA-mRNA showed that 4 known miRNAs and 11 genes formed 14 miRNA-target mRNA pairs were co-expressed in expansion stage. miRNA160, miRNA396, miRNA535 and miRNA5021 may be involved in complex network to regulate cell division and differentiation in yam during its expansion stage.

Conclusion

The mRNA and miRNA datasets presented here identified a subset of candidate genes and miRNAs that are putatively associated with tuber expansion in yam, a hypothetical model of genetic regulatory network associated with tuber expansion in yam was put forward, which may provide a foundation for molecular regulatory mechanism researching on tuber expansion in Dioscorea species.

The peu-miR160a-PeARF17.1/PeARF17.2 module participates in the adventitious root development of poplar

Deep roots give rise to flourishing leaves, and the two complement each other. However, the genetic mechanisms underlying adventitious rooting for forest trees have remained elusive. In this study, we verified that peu-miR160a targets six poplar genes AUXIN RESPONSE FACTORS (ARFs), PeARF10.1, PeARF16.1, PeARF16.2, PeARF16.3, PeARF17.1 and PeARF17.2, using 5'RLM-RACE. Interaction experiments with peu-miR160a and PeARFs in poplar protoplasts further confirmed that peu-miR160a targets and negatively regulates the six PeARFs. Peu-miR160a and its target genes exhibited obvious temporal expression in different stages of adventitious root development, and they could also be induced by IAA and abscisic acid. Peu-miR160a-overexpressing lines exhibited a significant shortening of adventitious root length, an increase in the number of lateral roots, severe dwarfing and shortened internodes. In addition, the overexpression of PeARF17.1 or mPeARF17.2 (peu-miR160a-resistant version of PeARF17.2) significantly increased the number of adventitious roots. Furthermore, PeARF17.1-overexpressing lines had multiple branches with no visible trunk, although the adventitious root length of the PeARF17.1-overexpressing lines was significantly increased. Our findings reveal that the peu-miR160a - PeARF17.1/PeARF17.2 module is an important regulator involved in the development of the adventitious roots of poplar.

MiRNA-target pairs regulate adventitious rooting in Populus: a functional role for miR167a and its target Auxin response factor 8

The ability of a plant to form roots from its non-root tissues is ecologically advantageous during rapid adaptation to a changing environment. Although this biological phenomenon has been widely utilized for cuttings in many economically important agronomic and tree species, its genetic and developmental mechanisms have been poorly understood. In this study, we conducted an association analysis of small RNAs, the degradome and the transcriptome of adventitious rooting in poplar softwood cuttings, which revealed that 373 miRNA-target pairs were detected. Of these, 72 significantly differentially expressed targets were screened as likely to modulate adventitious root (AR) development, in conjunction with plant hormone signal transduction. Poplar miR167a and its targets PeARF6s and PeARF8s were subjected to functional verification of their ability to mediate plant growth and hormone signal transduction. Overexpression of miR167a inhibited target transcripts and improved lateral root (LR) development in poplar, while overexpressing PeARF8.1mut increased AR numbers and slightly inhibited LR development. Taken together, these results suggest that miR167a-PeARF8.1 modules play crucial roles in regulating AR and LR development in poplar and improve the adaptation of poplar to more complex environments.

Analysis of the MIR160 gene family and the role of MIR160a_A05 in regulating fiber length in cotton

The MIR160 family in Gossypium hirsutum and G. barbadense was characterized, and miR160a_A05 was found to increase cotton-fiber length by downregulating its target gene (ARF17) and several GH3 genes. Cotton fiber is the most important raw material for the textile industry. MicroRNAs are involved in regulating cotton-fiber development, but a role in fiber elongation has not been demonstrated. In this study, miR160a was found to be differentially expressed in elongating fibers between two interspecific (between Gossypium hirsutum and G. barbadense) backcross inbred lines (BILs) with different fiber lengths. The gene MIR160 colocalized with a previously mapped fiber-length quantitative trait locus. Its target gene ARF17 was differentially expressed between the two BILs during fiber elongation, but in the inverse fashion. Bioinformatics was used to analyze the MIR160 family in both G. hirsutum and G. barbadense. Moreover, qRT-PCR analysis identified MIR160a as the functional MIR160 gene encoding the miR160a precursor during fiber elongation. Using virus-induced gene silencing and overexpression, overexpressed MIR160a_A05 resulted in significantly longer fibers compared with wild type, whereas suppression of miR160 resulted in significantly shorter fibers. Expression levels of the target gene auxin-response factor 17 (ARF17) and related genes GH3 in the two BILs and/or the virus-infected plants demonstrated similar changes in response to modulation of miR160a level. Finally, overexpression or suppression of miR160 increased or decreased, respectively, the cellular level of indole-3-acetic acid, which is involved in fiber elongation. These results describe a specific regulatory mechanism for fiber elongation in cotton that can be utilized for future crop improvement.

MicroRNA Dynamics and Functions During Arabidopsis Embryogenesis

MicroRNAs (miRNAs) are short noncoding RNAs that mediate the repression of target transcripts in plants and animals. Although miRNAs are required throughout plant development, relatively little is known regarding their embryonic functions. To systematically characterize embryonic miRNAs in Arabidopsis (Arabidopsis thaliana), we developed or applied high-throughput sequencing-based methods to profile hundreds of miRNAs and associated targets throughout embryogenesis. We discovered dozens of miRNAs that dynamically cleave and repress target transcripts, including 30 that encode transcription factors. Transcriptome analyses indicated that these miRNA:target interactions have profound effects on embryonic gene expression programs. Moreover, we demonstrated that the miRNA-mediated repression of six transcription factors are individually required for proper division patterns of various embryonic cell lineages. These data indicate that the miRNA-directed repression of multiple transcription factors is critically important for the establishment of the plant body plan, and they provide a foundation to further investigate how miRNAs contribute to these initial cellular differentiation events.

A transcriptomic profile of topping responsive non-coding RNAs in tobacco roots (Nicotiana tabacum)

BACKGROUND

Non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long ncRNAs (lncRNAs) and circular RNAs (circRNAs), accomplish remarkable variety of biological functions. However, the composition of ncRNAs and their interactions with coding RNAs in modulating and controlling of cellular process in plants is largely unknown. Using a diverse group of high-throughput sequencing strategies, the mRNA, miRNA, lncRNA and circRNA compositions of tobacco (Nicotiana tabacum) roots determined and their alteration and potential biological functions in response to topping treatment analyzed.

RESULTS

A total of 688 miRNAs, 7423 non-redundant lncRNAs and 12,414 circRNAs were identified, among which, some selected differentially expressed RNAs were verified by quantitative real-time PCR. Using the differentially expressed RNAs, a co-expression network was established that included all four types of RNAs. The number of circRNAs identified were higher than that of miRNAs and lncRNAs, but only two circRNAs were present in the co-expression network. LncRNAs appear to be the most active ncRNAs based on their numbers presented in the co-expression network, but none of them seems to be an eTM (endogenous Target Mimicry) of miRNAs. Integrated with analyses of sequence interaction, several mRNA-circRNA-miRNA interaction networks with a potential role in the regulation of nicotine biosynthesis were uncovered, including a QS-circQS-miR6024 interaction network. In this network miR6024 was significantly down-regulated, while the expression levels of its two targets, circQS and its host gene QS, were sharply increased following the topping treatment.

CONCLUSIONS

These results illustrated the transcriptomic profiles of tobacco roots, the organ responsible for nicotine biosynthesis. mRNAs always play the most important roles, while ncRNAs are also expressed extensively for topping treatment response, especially circRNAs are the most activated in the ncRNA pool. These studies also provided insights on the coordinated regulation module of coding and non-coding RNAs in a single plant biological sample. The findings reported here indicate that ncRNAs appear to form interaction complex for the regulation of stress response forming regulation networks with transcripts involved in nicotine biosynthesis in tobacco.

Transcriptome Analysis Implicates Involvement of Long Noncoding RNAs in Cytoplasmic Male Sterility and Fertility Restoration in Cotton

The cytoplasmic male sterility (CMS)/restorer-of-fertility system is an important tool to exploit heterosis during commercially hybrid seed production. The importance of long noncoding RNAs (lncRNAs) in plant development is recognized, but few analyses of lncRNAs during anther development of three-line hybrid cotton (CMS-D2 line A, maintainer line B, restorer-of-fertility line R) have been reported. Here, we performed transcriptome sequencing during anther development in three-line hybrid cotton. A total of 80,695 lncRNAs were identified, in which 43,347 and 44,739 lncRNAs were differentially expressed in A–B and A–R comparisons, respectively. These lncRNAs represent functional candidates involved in CMS and fertility restoration. GO analysis indicated that cellular hormone metabolic processes and oxidation–reduction reaction processes might be involved in CMS, and cellular component morphogenesis and small molecular biosynthetic processes might participate in fertility restoration. Additionally, 63 lncRNAs were identified as putative precursors of 35 miRNAs, and quantitative reverse transcription polymerase chain reaction (qRT-PCR) showed a similar expression pattern to RNA-seq data. Furthermore, construction of lncRNA regulatory networks indicated that several miRNA–lncRNA–mRNA networks might be involved in CMS and fertility restoration. Our findings provide systematic identification of lncRNAs during anther development and lays a solid foundation for the regulatory mechanisms and utilization in hybrid cotton breeding.

The F-box protein HAWAIIAN SKIRT is required for mimicry target-induced microRNA degradation in Arabidopsis

Mimicry target-directed microRNA degradation is widespread and highly conserved among eukaryotes. However, little is known about its mechanism of action. In this letter, by using STTM160 (target mimic of miR160) as a reporter, we show that dysfunction of HAWAIIAN SKIRT (HWS) suppresses the pleiotropic phenotype of STTM160. Small RNA sequencing and Northern blot analyses suggested that HWS only affects a subset of microRNAs. Intriguingly, we identified a stable coexistence of miR160/miR399 and their mimicry targets within the AGO1 complex when HWS is compromised, pointing to a possible role of HWS in the clearance of RNA-induced silencing complexes associated with mimicry target.

Enhancement of ABA Sensitivity Through Conditional Expression of the ARF10 Gene in Brassica juncea Reveals Fertile Plants with Tolerance Against Alternaria brassicicola

Concomitant increase of auxin-responsive factors ARF16 and ARF17, along with enhanced expression of ARF10 in resistant Sinapis alba compared with that in susceptible Brassica juncea upon challenge with Alternaria brassicicola, revealed that abscisic acid (ABA)-auxin crosstalk is a critical factor for resistance response. Here, we induced the ABA response through conditional expression of ARF10 in B. juncea using the A. brassicicola-inducible GH3.3 promoter. Induced ABA sensitivity caused by conditional expression of ARF10 in transgenic B. juncea resulted in tolerance against A. brassicicola and led to enhanced expression of several ABA-responsive genes without affecting the auxin biosynthetic gene expression. Compared with ABI3 and ABI4, ABI5 showed maximum upregulation in the most tolerant transgenic lines upon pathogen challenge. Moreover, elevated expression of ARF10 by different means revealed a direct correlation between ARF10 expression and the induction of ABI5 protein in B. juncea. Through in vitro DNA-protein experiments and chromosome immunoprecipitation using the ARF10 antibody, we demonstrated that ARF10 interacts with the auxin-responsive elements of the ABI5 promoter. This suggests that ARF10 may function as a modulator of ABI5 to induce ABA sensitivity and mediate the resistance response against A. brassicicola.

Identification of Fusarium oxysporum f. sp. cubense tropical race 4 (Foc TR4) responsive miRNAs in banana root

The fungus, Fusarium oxysporum f. sp. cubense (Foc), is the causal agent of Fusarium wilt disease, which is the most serious disease affecting the whole banana industry. Although extensive studies have characterized many Foc-responsive genes in banana, the molecular mechanisms on microRNA level underlying both banana defense and Foc pathogenesis are not yet fully understood. In this study, we aimed to reveal the role of miRNA during banana-Foc TR4 interactions. Illumina sequencing was used to reveal the changes in small RNAome profiles in roots of Foc TR4-inoculated ‘Tianbaojiao’ banana (Musa acuminata cv. Tianbaojiao) in the early stages (i.e. 5 h, 10 h and 25 h post Foc TR4 inoculation, respectively). The expression of some differentially expressed (DE) miRNAs and their predicted target genes was studied by using quantitative real time PCR (qRT-PCR). Totally, 254 known miRNAs from 31 miRNA families and 28 novel miRNAs were identified. Differential expression analysis identified 84, 77 and 74 DE miRNAs at the three respective Foc TR4 infection time points compared with control healthy banana (CK). GO and KEGG analysis revealed that most of the predicted target genes of DE miRNAs (DET) were implicated in peroxisome, fatty acid metabolism, auxin-activated signaling pathway, sulfur metabolism, lignin metabolism and so on, and many known stress responsive genes were identified to be DETs. Moreover, expected inverse correlations were confirmed between some miRNA and their corresponding target genes by using qRT-PCR analysis. Our study revealed that miRNA play important regulatory roles during the banana-Foc TR4 interaction by regulating peroxidase, fatty acid metabolism, auxin signaling, sulfur metabolism, lignin metabolism related genes and many known stress responsive genes.

Design and Synthesis of Novel 4-Hydroxyl-3-(2-phenoxyacetyl)-pyran-2-one Derivatives for Use as Herbicides and Evaluation of Their Mode of Action

In order to develop a novel herbicide containing the β-triketone motif, a series of 4-hydroxyl-3-(2-phenoxyacetyl)-pyran-2-one derivatives were designed and synthesized. The bioassay results showed that compound II15 had good pre-emergent herbicidal activity even at a dosage of 187.5 g ha^-1. Moreover, compound II15 showed a broader spectrum of weed control when compared with a commercial herbicide 2,4-dichlorophenoxyacetic acid (2,4-D), and displayed good crop safety to Triticum aestivum L. and Zea mays Linn. when applied at 375 g ha^-1 under pre-emergence conditions, which indicated its great potential as a herbicide. More importantly, studying the molecular mode of action of compound II15 revealed that the novel triketone structure is a proherbicide of its corresponding phenoxyacetic acid auxin herbicide, which has a herbicidal mechanism similar to that of 2,4-D. The present work indicates that the 4-hydroxyl-3-(2-phenoxyacetyl)-pyran-2-one motif may be a potential lead structure for further development of novel auxin-type herbicides.

Small RNA profiling from meiotic and post-meiotic anthers reveals prospective miRNA-target modules for engineering male fertility in sorghum

Understanding male gametophyte development is essential to augment hybrid production in sorghum. Although small RNAs are known to critically influence anther/pollen development, their roles in sorghum reproduction have not been deciphered yet. Here, we report small RNA profiling and high-confidence annotation of microRNAs (miRNAs) from meiotic and post-meiotic anthers in sorghum. We identified 262 miRNAs (82 known and 180 novel), out of which 58 (35 known and 23 novel) exhibited differential expression between two stages. Out of 35 differentially expressed known miRNAs, 13 are known to regulate anther/pollen development in other plant species. We also demonstrated conserved spatiotemporal patterns of 21- and 24-nt phasiRNAs and their respective triggers, miR2118 and miR2275, in sorghum anthers as evidenced in other monocots. miRNA target identification yielded 5622 modules, of which 46 modules comprising 16 known and 8 novel miRNA families with 38 target genes are prospective candidates for engineering male fertility in grasses.

Gene Expression Profiles and microRNA Regulation Networks in Tiller Primordia, Stem Tips, and Young Spikes of Wheat Guomai 301

Tillering and spike differentiation are two key events for wheat (Triticum aestivum L.). A study on the transcriptomes and microRNA group profiles of wheat at the two key developmental stages will bring insight into the molecular regulation mechanisms. Guomai 301 is a representative excellent new high yield wheat cultivar in the Henan province in China. The transcriptomes and microRNA (miRNA) groups of tiller primordia (TPs), stem tips (STs), and young spikes (YSs) in Guomai 301 were compared to each other. A total of 1741 tillering specifically expressed and 281 early spikes differentiating specifically expressed differentially expressed genes (DEGs) were identified. Six major expression profile clusters of tissue-specific DEGs for the three tissues were classified by gene co-expression analysis using K-means cluster. The ribosome (ko03010), photosynthesis-antenna proteins (ko00196), and plant hormone signal transduction (ko04075) were the main metabolic pathways in TPs, STs, and YSs, respectively. Similarly, 67 TP specifically expressed and 19 YS specifically expressed differentially expressed miRNAs were identified, 65 of them were novel. The roles of 3 well known miRNAs, tae-miR156, tae-miR164, and tae-miR167a, in post-transcriptional regulation were similar to that of other researches. There were 651 significant negative miRNA-mRNA interaction pairs in TPs and YSs, involving 63 differentially expressed miRNAs (fold change > 4) and 416 differentially expressed mRNAs. Among them 12 key known miRNAs and 16 novel miRNAs were further analyzed, and miRNA-mRNA regulatory networks during tillering and early spike differentiating were established.

Arabidopsis ANAC092 regulates auxin-mediated root development by binding to the ARF8 and PIN4 promoters

Auxin is an important plant hormone that is essential for growth and development due to its effects on organogenesis, morphogenesis, tropisms, and apical dominance. The functional diversity of auxin highlights the importance of its biosynthesis, transport, and associated responses. In this study, we show that a NAC transcription factor, ANAC092 (also named AtNAC2 and ORESARA1), known to positively regulate leaf senescence and contribute to abiotic stress responses, also affects primary root development. Plants overexpressing ANAC092 had altered root meristem lengths and shorter primary roots compared with the wild-type control. Additionally, expression of the proANAC092::GUS was strongly induced by indole-3-acetic acid. Quantitative real-time RT-PCR (qRT-PCR) analysis revealed that the YUCCA2, PIN, and ARF expression levels were downregulated in ANAC092-overexpressing plants. Moreover, yeast one-hybrid and chromatin immunoprecipitation assays confirmed that ANAC092 binds to the promoters of AUXIN RESPONSE FACTOR 8 (ARF8) and PIN-FORMED 4 (PIN4). Furthermore, a dual-luciferase assay indicated that ANAC092 decreases ARF8 and PIN4 promoter activities. We also applied a CRISPR/Cas9 system to mutate ANAC092. The roots of three of the analyzed mutants were longer than normal. Collectively, our findings indicate that ANAC092 negatively affects root development by controlling the auxin pathway.

Tackling Plant Phosphate Starvation by the Roots

Plant responses to phosphate deprivation encompass a wide range of strategies, varying from altering root system architecture, entering symbiotic interactions to excreting root exudates for phosphorous release, and recycling of internal phosphate. These processes are tightly controlled by a complex network of proteins that are specifically upregulated upon phosphate starvation. Although the different effects of phosphate starvation have been intensely studied, the full extent of its contribution to altered root system architecture remains unclear. In this review, we focus on the effect of phosphate starvation on the developmental processes that shape the plant root system and their underlying molecular pathways.

PtrARF2.1 Is Involved in Regulation of Leaf Development and Lignin Biosynthesis in Poplar Trees

Auxin response factors (ARFs) are important regulators modulating the expression of auxin-responsive genes in various biological processes in plants. In the Populus genome, a total of 39 ARF members have been identified, but their detailed functions are still unclear. In this study, six poplar auxin response factor 2 (PtrARF2) members were isolated from P. trichocarpa. Expression pattern analysis showed that PtrARF2.1 is highly expressed in leaf tissues compared with other PtrARF2 genes and significantly repressed by exogenous auxin treatment. PtrARF2.1 is a nuclear-localized protein without transcriptional activation activity. Knockdown of PtrARF2.1 by RNA interference (RNAi) in poplars led to the dwarf plant, altered leaf shape, and reduced size of the leaf blade, while overexpression of PtrARF2.1 resulted in a slight reduction in plant height and the similar leaf phenotype in contrast to the wildtype. Furthermore, histological staining analysis revealed an ectopic deposition of lignin in leaf veins and petioles of PtrARF2.1-RNAi lines. RNA-Seq analysis showed that 74 differential expression genes (DEGs) belonging to 12 transcription factor families, such as NAM, ATAF and CUC (NAC), v-myb avian myeloblastosis viral oncogene homolog (MYB), ethylene response factors (ERF) and basic helix–loop–helix (bHLH), were identified in PtrARF2.1-RNAi leaves and other 24 DEGs were associated with the lignin biosynthetic pathway. Altogether, the data indicate that PtrARF2.1 plays an important role in regulating leaf development and influences the lignin biosynthesis in poplars.

Drought Response in Rice: The miRNA Story

As a semi-aquatic plant, rice requires water for proper growth, development, and orientation of physiological processes. Stress is induced at the cellular and molecular level when rice is exposed to drought or periods of low water availability. Plants have existing defense mechanisms in planta that respond to stress. In this review we examine the role played by miRNAs in the regulation and control of drought stress in rice through a summary of molecular studies conducted on miRNAs with emphasis on their contribution to drought regulatory networks in comparison to other plant systems. The interaction between miRNAs, target genes, transcription factors and their respective roles in drought-induced stresses is elaborated. The cross talk involved in controlling drought stress responses through the up and down regulation of targets encoding regulatory and functional proteins is highlighted. The information contained herein can further be explored to identify targets for crop improvement in the future.

Unravelling the MicroRNA-Mediated Gene Regulation in Developing Pongamia Seeds by High-Throughput Small RNA Profiling

Pongamia (Millettia pinnata syn. Pongamia pinnata) is a multipurpose biofuel tree which can withstand a variety of abiotic stresses. Commercial applications of Pongamia trees may substantially benefit from improvements in their oil-seed productivity, which is governed by complex regulatory mechanisms underlying seed development. MicroRNAs (miRNAs) are important molecular regulators of plant development, while relatively little is known about their roles in seed development, especially for woody plants. In this study, we identified 236 conserved miRNAs within 49 families and 143 novel miRNAs via deep sequencing of Pongamia seeds sampled at three developmental phases. For these miRNAs, 1327 target genes were computationally predicted. Furthermore, 115 differentially expressed miRNAs (DEmiRs) between successive developmental phases were sorted out. The DEmiR-targeted genes were preferentially enriched in the functional categories associated with DNA damage repair and photosynthesis. The combined analyses of expression profiles for DEmiRs and functional annotations for their target genes revealed the involvements of both conserved and novel miRNA-target modules in Pongamia seed development. Quantitative Real-Time PCR validated the expression changes of 15 DEmiRs as well as the opposite expression changes of six targets. These results provide valuable miRNA candidates for further functional characterization and breeding practice in Pongamia and other oilseed plants.

Comparative profile analysis reveals differentially expressed microRNAs regulate anther and pollen development in kenaf cytoplasmic male sterility line

Cytoplasmic male sterility (CMS) is advantageous in extensive crop breeding and represents a perfect model for understanding anther and pollen development research. MicroRNAs (miRNAs) play key roles in regulating various biological processes. However, the miRNA-mediated regulatory network in kenaf CMS occurrence remains largely unknown. In the present study, a comparative deep sequencing approach was used to investigate the miRNAs and their roles in regulating anther and pollen development during CMS occurrence. We identified 283 known and 46 new candidate miRNAs in kenaf anther. A total of 67 differentially expressed miRNAs (DEMs) were discovered between CMS and its maintainer line. Among them, 40 and 27 miRNAs were up- and downregulated, respectively. These 67 DEMs were predicted to target 189 genes. Validation of DEMs and putative target genes were confirmed by using real-time quantitative PCR. In addition, a potential miRNA-mediated regulatory network, which mainly involves the auxin signaling pathway, signal transduction, glycolysis and energy metabolism, gene expression, transmembrane transport, protein modification and metabolism, and floral development, that mediates anther development during CMS occurrence was proposed. Taken together, our findings provide a better understanding of the molecular mechanism of miRNA regulation in pollen development and CMS occurrence in kenaf.

MicroRNA-guided regulation of heat stress response in wheat

Background

With rising global temperature, understanding plants’ adaptation to heat stress has implications in plant breeding. MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs guiding gene expression at the post-transcriptional level. In this study, small RNAs and the degradome (parallel analysis of RNA ends) of leaf tissues collected from control and heat-stressed wheat plants immediately at the end of the stress period and 1 and 4 days later were analysed.

Results

Sequencing of 24 small RNA libraries produced 55.2 M reads while 404 M reads were obtained from the corresponding 24 PARE libraries. From these, 202 miRNAs were ascertained, of which mature miRNA evidence was obtained for 104 and 36 were found to be differentially expressed after heat stress. The PARE analysis identified 589 transcripts targeted by 84 of the ascertained miRNAs. PARE sequencing validated the targets of the conserved members of miRNA156, miR166 and miR393 families as squamosa promoter-binding-like, homeobox leucine-zipper and transport inhibitor responsive proteins, respectively. Heat stress responsive miRNA targeted superoxide dismutases and an array of homeobox leucine-zipper proteins, F-box proteins and protein kinases. Query of miRNA targets to interactome databases revealed a predominant association of stress responses such as signalling, antioxidant activity and ubiquitination to superoxide dismutases, F-box proteins, pentatricopeptide repeat-containing proteins and mitochondrial transcription termination factor-like proteins.

Conclusion

The interlaced data set generated in this study identified and validated heat stress regulated miRNAs and their target genes associated with thermotolerance. Such accurate identification and validation of miRNAs and their target genes are essential to develop novel regulatory gene-based breeding strategies.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5799-6) contains supplementary material, which is available to authorized users.

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.

Sweet cherry fruit miRNAs and effect of high CO2 on the profile associated with ripening

157 known and 55 novel miRNAs were found in sweet cherry fruit. MiRNA target genes involved in fruit ripening and the differentially expressed miRNAs under CO₂ treatment were identified. MicroRNAs (miRNAs) are short non-coding RNAs and play important functions in many biological processes, including fruit ripening and senescence. In the current study, the high-throughput sequencing and bioinformatics methods were implemented to decipher the miRNAs landscape in sweet cherry fruit. A total of 157 known miRNAs belonging to 50 families and 55 putative novel miRNAs were found. Target genes of the miRNAs were predicted and genes involved in fruit ripening were found, including F-box proteins and TFs such as SPL, TCP, NAC, MYB, ARF and AP2/ERF. And these target genes were further confirmed by degradome sequencing. A regulatory network model was constructed to uncover the miRNAs and their targets involved in fruit ripening and senescence. Importantly, elevated carbon dioxide can significantly postpone the ripening and senescence of sweet cherry fruit and the differentially expressed miRNAs exposed to CO₂ were identified. These miRNAs included miR482j, miR6275, miR164, miR166, miR171, miR393, miR858, miR3627a, miR6284, miR6289 and miR7122b, and some of their functions were linked to fruit ripening. This study was the first report to profile miRNAs in sweet cherry fruit and it would provide more information for further study of miRNA roles in the ripening processes and their regulation mechanism underlying the effects of high carbon dioxide treatment on fruit ripening.

Molecular characterization of virus-derived small RNAs in Nicotiana benthamiana plants infected with tobacco curly shoot virus and its β satellite

Tobacco curly shoot virus (TbCSV) is a monopartite DNA virus of the genus Begomovirus, which causes leaf curl symptoms in tobacco and tomato. The β satellite of TbCSV (TbCSB induces more severe symptoms and enhanced virus accumulation when co-infects the host plants with TbCSV. Small interfering RNAs derived from virus(vsiRNAs) induce disease symptoms and promote virus invasion by target and guide the degradation of host transcripts The vsiRNAs derived from TbCSV and TbCSV + TbCSB remained to be explored to elucidate the molecular mechanism of symptoms development in plants. In the present work, two libraries of small RNA from TbCSV-infected and TbCSV + TbCSB-infected N. benthamiana plants were constructed and the vsiRNAs in both samples shared the same characteristics. The size of the vsiRNAs ranged from 18 to 30 nucleotides (nt), with most of them being 21 or 22 nt, which accounted for 29.11% and 23.22% in TbCSV plants and 29.39% and 21.82% in TbCSV + TbCSV plants, respectively. The vsiRNAs with A/U bias at the first site were abundant in both the TbCSV-treated and TbCSV + TbCSB-treated plants. It is discovered that the vsiRNAs continuously, but heterogeneously, distributed through bothe the TbCSV and TbCSB sequences. And the distribution profiles were similar in both the treatments such as mainly in the overlapping region of the AC2/AC3 coding sequences. The host transcripts targeted by vsiRNAs were predicted, and the targeted genes were found to be involved in varied biological processes. It is indicated that the presence of TbCSB does not significantly affect the production of vsiRNAs from TbCSV in plants, the distribution hotsopt of TbCSV vsiRNAs could be useful in designing effective targets for TbCSV resistance exploiting RNA interference.

An Essential Role for miRNA167 in Maternal Control of Embryonic and Seed Development

Maternal cells play a critical role in ensuring the normal development of embryos, endosperms, and seeds. Mutations that disrupt the maternal control of embryogenesis and seed development are difficult to identify. Here, we completely deleted four MICRORNA167 (MIR167) genes in Arabidopsis (Arabidopsis thaliana) using a clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein9 (Cas9) genome-editing technology. We found that plants with a deletion of MIR167A phenocopied plants overexpressing miRNA167-resistant versions of Auxin Response Factor6 (ARF6) or ARF8, two miRNA167 targets. Both the mir167a mutant and the ARF overexpression lines were defective in anther dehiscence and ovule development. Serendipitously, we found that the mir167a (♀) × wild type (♂) crosses failed to produce normal embryos and endosperms, despite the findings that embryos with either mir167a^+/- or mir167a^-/- genotypes developed normally when mir167a^+/- plants were self-pollinated, revealing a central role of MIR167A in maternal control of seed development. The mir167a phenotype is 100% penetrant, providing a great genetic tool for studying the roles of miRNAs and auxin in maternal control. Moreover, we found that mir167a mutants flowered significantly later than wild-type plants, a phenotype that was not observed in the ARF overexpression lines. We show that the reproductive defects of mir167a mutants were suppressed by a decrease of activities of ARF6, ARF8, or both. Our results clearly demonstrate that MIR167A is the predominant MIR167 member in regulating Arabidopsis reproduction and that MIR167A acts as a maternal gene that functions largely through ARF6 and ARF8.

Differential Expression of MicroRNAs During Root Formation in Taxus Chinensis Var. mairei Cultivars

MicroRNAs (miRNAs) have been shown to play key roles in the regulation of plant growth and development by modifying the expression of their target genes. However, the influence of miRNAs on root formation and development in woody plants, such as Taxus chinensis, remains largely unknown. In the current study, we explored the phytohormone-response and nutrition-response miRNA expression profiles during T. chinensis rooting by quantitative real-time PCR (qPCR). We identified six phytohormone-response miRNAs, namely, miR164a, miR165, miR167a, miR171b, miR319, and miR391, and eight nutrition-response miRNAs, namely, miR169b, miR395a, miR399c, miR408, miR826, miR827, miR857, and miR2111a, that were differentially expressed at different rooting phases of T. chinensis. Using northern blot analysis of the putative target genes of these miRNAs, we detected the relative gene expression changes of the target genes. Taken together, our results suggest that miRNAs are involved in root formation of T. chinensis and that miRNAs may play important regulatory roles in primary root, crown root, and root hair formation by targeting phytohormone and/or nutrition response genes in T. chinensis. For the first time, these results expand our understanding of the molecular mechanisms of plant root formation and development in a conifer species.

Salinity-associated microRNAs and their potential roles in mediating salt tolerance in rice colonized by the endophytic root fungus Piriformospora indica

Piriformospora indica (P. indica), an endophytic root fungus, supports the growth and enhanced tolerance of plants to biotic and abiotic stresses. Several recent studies showed the significant role of small RNA (sRNA) molecules including microRNAs (miRNAs) in plant adaption to environmental stress, but little is known concerning the symbiosis-mediated salt stress tolerance regulated at miRNAs level. The overarching goal of this research is to elucidate the impact of miRNAs in regulating the P. indica-mediated salt tolerance in rice. Applying sRNA-seq analysis led to identify a set of 547 differentially abundant miRNAs in response to P. indica inoculation and salt stress. These included 206 rice-specific and 341 previously known miRNAs from other plant species. In silico analysis of miRNAs predictions of the differentially abundant miRNAs led to identifying of 193 putatively target genes, most of which were encoded either genes or transcription factors involved in nutrient uptake, sodium ion transporters, growth regulators, and auxin- responsive proteins. The rice-specific miRNAs targeted the transcription factors involved in the import of potassium ions into the root cells, the export of sodium ions, and plant growth and development. Interestingly, P. indica affected the differential abundance of miRNAs regulated genes and transcription factors linked to salt stress tolerance. Our data helps to understand the molecular basis of salt stress tolerance mediated by symbionts in plant and the potential impact of miRNAs for genetic improvement of rice varieties for tolerance to salt stress.

Micromanagement of Developmental and Stress-Induced Senescence: The Emerging Role of MicroRNAs

MicroRNAs are short (19⁻24-nucleotide-long), non-coding RNA molecules. They downregulate gene expression by triggering the cleavage or translational inhibition of complementary mRNAs. Senescence is a stage of development following growth completion and is dependent on the expression of specific genes. MicroRNAs control the gene expression responsible for plant competence to answer senescence signals. Therefore, they coordinate the juvenile-to-adult phase transition of the whole plant, the growth and senescence phase of each leaf, age-related cellular structure changes during vessel formation, and remobilization of resources occurring during senescence. MicroRNAs are also engaged in the ripening and postharvest senescence of agronomically important fruits. Moreover, the hormonal regulation of senescence requires microRNA contribution. Environmental cues, such as darkness or drought, induce senescence-like processes in which microRNAs also play regulatory roles. In this review, we discuss recent findings concerning the role of microRNAs in the senescence of various plant species.

The interaction between miR160 and miR165/166 in the control of leaf development and drought tolerance in Arabidopsis

MicroRNAs (miRNAs) are a class of non-coding RNAs that play important roles in plant development and abiotic stresses. To date, studies have mainly focused on the roles of individual miRNAs, however, a few have addressed the interactions among multiple miRNAs. In this study, we investigated the interplay and regulatory circuit between miR160 and miR165/166 and its effect on leaf development and drought tolerance in Arabidopsis using Short Tandem Target Mimic (STTM). By crossing STTM160 Arabidopsis with STTM165/166, we successfully generated a double mutant of miR160 and miR165/166. The double mutant plants exhibited a series of compromised phenotypes in leaf development and drought tolerance in comparison to phenotypic alterations in the single STTM lines. RNA-seq and qRT-PCR analyses suggested that the expression levels of auxin and ABA signaling genes in the STTM-directed double mutant were compromised compared to the two single mutants. Our results also suggested that miR160-directed regulation of auxin response factors (ARFs) contribute to leaf development via auxin signaling genes, whereas miR165/166- mediated HD-ZIP IIIs regulation confers drought tolerance through ABA signaling. Our studies further indicated that ARFs and HD-ZIP IIIs may play opposite roles in the regulation of leaf development and drought tolerance that can be further applied to other crops for agronomic traits improvement.

Exploring miRNAs for developing climate-resilient crops: A perspective review

Climate changes and environmental stresses have significant implications on global crop production and necessitate developing crops that can withstand an array of climate changes and environmental perturbations such as irregular water-supplies leading to drought or water-logging, hyper soil-salinity, extreme and variable temperatures, ultraviolet radiations and metal stress. Plants have intricate molecular mechanisms to cope with these dynamic environmental changes, one of the most common and effective being the reprogramming of expression of stress-responsive genes. Plant microRNAs (miRNAs) have emerged as key post-transcriptional and translational regulators of gene-expression for modulation of stress implications. Recent reports are establishing their key roles in epigenetic regulations of stress/adaptive responses as well as in providing plants genome-stability. Several stress responsive miRNAs are being identified from different crop plants and miRNA-driven RNA-interference (RNAi) is turning into a technology of choice for improving crop traits and providing phenotypic plasticity in challenging environments. Here we presents a perspective review on exploration of miRNAs as potent targets for engineering crops that can withstand multi-stress environments via loss-/gain-of-function approaches. This review also shed a light on potential roles plant miRNAs play in genome-stability and their emergence as potent target for genome-editing. Current knowledge on plant miRNAs, their biogenesis, function, their targets, and latest developments in bioinformatics approaches for plant miRNAs are discussed. Though there are recent reviews discussing primarily the individual miRNAs responsive to single stress factors, however, considering practical limitation of this approach, special emphasis is given in this review on miRNAs involved in responses and adaptation of plants to multi-stress environments including at epigenetic and/or epigenomic levels.

Death of female flower microsporocytes progresses independently of meiosis-like process and can be accelerated by specific transcripts in Asparagus officinalis

Asparagus officinalis (garden asparagus) is a dioecious perennial crop, and the dioecy (i.e., sex) of A. officinalis can affect its productivity. In A. officinalis, flower anthers in female plants fail to accumulate callose around microsporocytes, fail to complete meiosis, and degenerate due to cell death. Although 13 genes have been implicated in the anther development of male and female flowers, it is unclear how these genes regulate the cell death in female flower anthers. The aim of this study was to narrow down factors involved in this process. TUNEL staining and Feulgen staining of female flower microsporocytes suggest that female microsporocytes enter a previously undetected meiosis-like process, and that the cell death occurs independently of this meiosis-like process, excluding the possibility that the cell death is caused by the cessation of meiosis. RNA sequencing with individual floral organs (tepals, pistils and stamens) revealed that several genes possibly regulating the cell death, such as metacaspase genes and a Bax inhibitor-1 gene, are differentially regulated between female and male flower anthers, and that genes involved in callose accumulation are up-regulated only in male flower anthers. These genes are likely involved in regulating the cell death in female flower anthers in A. officinalis.

High-throughput sequencing and differential expression analysis of miRNAs in response to Brassinosteroid treatment in Arabidopsis thaliana

Brassinosteroids are a class of phytohormones that play crucial roles in improving stress tolerance in plants. Many biochemical and physiological changes in response to abiotic stress are related to regulation of gene expression and accumulation of associated proteins. MicroRNAs (miRNAs) are class of small non-coding RNAs that regulate gene expression post-transcriptionally. Roles of these regulatory RNAs in brassinosteroid (BR) signalling have however remained elusive. In this study using high-throughput small RNA sequencing method, we present a comprehensive compilation of BR-induced differentially expressed microRNAs in root and shoots of Arabidopsis thaliana seedlings. We identified 229 known miRNAs belonging to 102 families and 27 novel miRNAs that express in response to exogenous BR treatment. Out of 102 families, miRNAs belonging to known 48 families and out of 27 novel miRNAs, 23 were observed to be differentially expressed in response to BR treatment. Among the conserved miRNAs, all members of miR169 were observed to be downregulated in both shoot and root samples. While, auxin-responsive factors were predicted to be direct targets of some novel miRNAs that are upregulated in shoots and suppressed in roots. The BR-responsive tissue-specific miRNome characterized in this study can be used as a starting point by investigators for functional validation studies that will shed light on the underlying molecular mechanism of BR-mediated stress tolerance at the level of post-transcriptional gene regulation.

L-Cysteine desulfhydrase-dependent hydrogen sulfide is required for methane-induced lateral root formation

Methane-triggered lateral root formation is not only a universal event, but also dependent on L-cysteine desulfhydrase-dependent hydrogen sulfide signaling. Whether or how methane (CH₄) triggers lateral root (LR) formation has not been elucidated. In this report, CH₄ induction of lateral rooting and the role of hydrogen sulfide (H₂S) were dissected in tomato and Arabidopsis by using physiological, anatomical, molecular, and genetic approaches. First, we discovered that CH₄ induction of lateral rooting is a universal event. Exogenously applied CH₄ not only triggered tomato lateral rooting, but also increased activities of L-cysteine desulfhydrase (DES; a major synthetic enzyme of H₂S) and induced endogenous H₂S production, and contrasting responses were observed in the presence of hypotaurine (HT; a scavenger of H₂S) or DL-propargylglycine (PAG; an inhibitor of DES) alone. CH₄-triggered lateral rooting were sensitive to the inhibition of endogenous H₂S with HT or PAG. The changes in the transcripts of representative cell cycle regulatory genes, miRNA and its target genes were matched with above phenotypes. In the presence of CH₄, Arabidopsis mutant Atdes1 exhibited defects in lateral rooting, compared with the wild-type. Molecular evidence showed that the transcriptional profiles of representative target genes modulated by CH₄ in wild-type plants were impaired in Atdes1 mutant. Overall, our data demonstrate the main branch of the DES-dependent H₂S signaling cascade in CH₄-triggered LR formation.

Identification of lncRNAs involved in rice ovule development and female gametophyte abortion by genome-wide screening and functional analysis

Background

As important female reproductive tissues, the rice (Oryza sativa L.) ovule and female gametophyte is significant in terms of their fertility. Long noncoding RNAs (lncRNAs) play important and wide-ranging roles in the growth and development of plants and have become a major research focus in recent years. Therefore, we explored the characterization and expression change of lncRNAs during ovule development and female gametophytic abortion.

Results

In our study, whole-transcriptome strand-specific RNA sequencing (ssRNA-seq) was performed in the ovules of a high-frequency female-sterile rice line (fsv1) and a wild-type rice line (Gui99) at the megaspore mother cell meiosis stage (stage 1), functional megaspore mitosis stage (stage 2) and female gametophyte mature stage (stage 3). By comparing two rice lines, we identified 152, 233, and 197 differentially expressed lncRNAs at the three ovule developmental stages. Functional analysis of the coherent target genes of these differentially expressed lncRNAs indicated that many lncRNAs participate in multiple pathways such as hormone and cellular metabolism and signal transduction. Moreover, there were many differentially expressed lncRNAs acting as the precursors of some miRNAs that are involved in the development of ovules and female gametophytes. In addition, we have found that lncRNAs can act as decoys, competing with mRNAs for binding to miRNAs to maintain the normal expression of genes related to ovule and female gametophyte development.

Conclusion

These results provide important clues for elucidating the female gametophyte abortion mechanism in rice. This study also expands our understanding about the biological functions of lncRNAs and the annotation of the rice genome.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5442-6) contains supplementary material, which is available to authorized users.

Target-mimicry based diminution of miRNA167 reinforced flowering-time phenotypes in tobacco via spatial-transcriptional biases of flowering-associated miRNAs

Evolutionarily conserved microRNAs such as miR156, miR159, miR167 and miR172 tightly regulate the extensive array of gene expression during flowering in plants, through instant and long-term alterations in the expression of their target genes. Here we employed a novel target-mimicry approach for the diminution of auxin signalling regulator miRNA167 by developing mimic-transgenic lines in tobacco, to investigate the transcriptional biases of flowering-associated miRNAs in apical and floral meristematic tissues and their phenotypic implications. Recorded morpho-alterations such as uneven flowering-time phenotypes, anomalous floral organ formation, and large variations in the seed forming characteristics permitted us to determine the consequence of the extent of miR167 expression diminution accompanying the transcriptional biases of interrelated miRNAs. We demonstrate that percent diminution of miR167 gene expression is proportionally associated with both early and late flowering-time phenotypes in mimic lines. Also, the associated miRNAs, miR156, miR159, and miR172 showed >90% transcriptional diminution in at least 'early-flowering' miR167 mimic lines. On contrary, low percentages of their respective diminution were recorded in 'late-flowering' lines. Evidently, the misexpression of miR156, miR159, and miR172 led to the over-expression of their respective target genes SPL9, AtMYB33-like and AP2 genes in mimic lines which resulted in assorted phenotypes. We describe the scope of spatial regulation of these microRNAs in floral bud tissues of mimic lines which showed negative- or very low (<25%) misexpression levels in early/late-flowering lines highlighting their roles in the acquisition of flowering mechanism. To our knowledge, this study represents the first characterization of transcriptional biases of flowering associated miRNAs in miR167-mimic lines and certainly augments our understanding of the importance of microRNA-mediated regulation of flowering in plants.

Functional Disruption of the Tomato Putative Ortholog of HAWAIIAN SKIRT Results in Facultative Parthenocarpy, Reduced Fertility and Leaf Morphological Defects

A number of plant microRNAs have been demonstrated to regulate developmental processes by integrating internal and environmental cues. Recently, the Arabidopsis thaliana F-box protein HAWAIIAN SKIRT (HWS) gene has been described for its role in miRNA biogenesis. We have isolated in a forward genetic screen a tomato (Solanum lycopersicum) line mutated in the putative ortholog of HWS. We show that the tomato hws-1 mutant exhibits reduction in leaflet serration, leaflet fusion, some degree of floral organ fusion, and alteration in miRNA levels, similarly to the original A. thaliana hws-1 mutant. We also describe novel phenotypes for hws such as facultative parthenocarpy, reduction in fertility and flowering delay. In slhws-1, the parthenocarpy trait is influenced by temperature, with higher parthenocarpy rate in warmer environmental conditions. Conversely, slhws-1 is able to produce seeds when grown in cooler environment. We show that the reduction in seed production in the mutant is mainly due to a defective male function and that the levels of several miRNAs are increased, in accordance with previous HWS studies, accounting for the abnormal leaf and floral phenotypes as well as the altered flowering and fruit development processes. This is the first study of HWS in fleshy fruit plant, providing new insights in the function of this gene in fruit development.

Deep sequencing identified potential miRNAs involved in defence response, stress and plant growth characteristics of wild genotypes of cardamom

Cardamom has long been used as a food flavouring agent and in ayurvedic medicines for mouth ulcers, digestive problems and even depression. Extensive occurrence of pests and diseases adversely affect its cultivation and result in substantial reductions in total production and productivity. Numerous studies revealed the significant role of miRNAs in plant biotic stress responses. In the current study, miRNA profiling of cultivar and wild cardamom genotypes was performed using an Ion Proton sequencer. We identified 161 potential miRNAs representing 42 families, including monocot/tissue-specific and 14 novel miRNAs in both genotypes. Significant differences in miRNA family abundance between the libraries were observed in read frequencies. A total of 19 miRNAs (from known miRNAs) displayed a twofold difference in expression between wild and cultivar genotypes. We found 1168 unique potential targets for 40 known miRNA families in wild and 1025 potential targets for 42 known miRNA families in cultivar genotypes. The differential expression analysis revealed that most miRNAs identified were highly expressed in cultivars and, furthermore, lower expression of miR169 and higher expression of miR529 in wild cardamom proved evidence that wild genotypes have stronger drought stress tolerance and floral development than cultivars. Potential targets predicted for the newly identified miRNAs from the miRNA libraries of wild and cultivar cardamom genotypes involved in metabolic and developmental processes and in response to various stimuli. qRT-PCR confirmed miRNAs were differentially expressed between wild and cultivar genotypes. Furthermore, four target genes were validated experimentally to confirm miRNA-mRNA target pairing using RNA ligase-mediated 5' Rapid Amplification of cDNA Ends (5'RLM-RACE) PCR.

miRNomes involved in imparting thermotolerance to crop plants

Thermal stress is one of the challenges to crop plants that negatively impacts crop yield. To overcome this ever-growing problem, utilization of regulatory mechanisms, especially microRNAs (miRNAs), that provide efficient and precise regulation in a targeted manner have been found to play determining roles. Besides their roles in plant growth and development, many recent studies have shown differential regulation of several miRNAs during abiotic stresses including heat stress (HS). Thus, understanding the underlying mechanism of miRNA-mediated gene expression during HS will enable researchers to exploit this regulatory mechanism to address HS responses. This review focuses on the miRNAs and regulatory networks that were involved in physiological, metabolic and morphological adaptations during HS in plant, specifically in crops. Illustrated examples including, the miR156-SPL, miR169-NF-YA5, miR395-APS/AST, miR396-WRKY, etc., have been discussed in specific relation to the crop plants. Further, we have also discussed the available plant miRNA databases and bioinformatics tools useful for miRNA identification and study of their regulatory role in response to HS. Finally, we have briefly discussed the future prospects about the miRNA-related mechanisms of HS for improving thermotolerance in crop plants.

Genome-wide identification and profiling of microRNAs in Paulownia tomentosa cambial tissues in response to seasonal changes

MicroRNAs (miRNAs), a group of endogenous small non-coding RNAs, have been shown to play essential roles in the regulation of gene expression at the post-transcriptional level. Although Paulownia tomentosa is an ecologically and economically important timber species due to its rapid growth, few efforts have focused on small RNAs (sRNAs) in the cambial tissues during winter and summer transition. In the present study, we identified 33 known miRNA families and 29 novel miRNAs which include 20 putative novel miRNAs* in P. tomentosa cambial tissues during winter and summer transition. Through differential expression analysis, we showed that 15 known miRNAs and 8 novel miRNAs were preferentially abundant in certain stage of cambial tissues. Based on the P. tomentosa mRNA transcriptome database, 1667 and 78 potential targets were predicted for 29 known and 20 novel miRNAs, respectively and the predicted targets are mostly transcription factors and functional genes. The targets of these miRNAs were enriched in "metabolic process" and "transcription regulation" by using Gene Ontology enrichment analysis. In addition, KEGG pathway analyses revealed the involvement of miRNAs in starch and sucrose metabolism and plant-pathogen interaction metabolism pathways. Noticeably, qRT-PCR expression analysis demonstrated that 9 miRNAs and their targets were existed a negative correlation in P. tomentosa cambial tissues. This study is the first to examine known and novel miRNAs and their potential targets in P. tomentosa cambial tissues during winter and summer transition and identify several candidate genes potentially regulating cambial phase transition, and thus provide a framework for further understanding of miRNAs functions in the regulation of cambial phase transition and wood formation in trees.

Molecular Control by Non-coding RNAs During Fruit Development: From Gynoecium Patterning to Fruit Ripening

Fruits are originated from the transition of a quiescent ovary to a fast-growing young fruit. The evolution of reproductive structures such as ovary and fruit has made seed dispersal easier, which is a key process for reproductive success in flowering plants. The complete fruit development and ripening are characterized by a remarkable phenotypic plasticity which is orchestrated by a myriad of genetic factors. In this context, transcriptional regulation by non-coding small (i.e., microRNAs) and long (lncRNAs) RNAs underlies important mechanisms controlling reproductive organ development. These mechanisms may act together and interact with other pathways (i.e., phytohormones) to regulate cell fate and coordinate reproductive organ development. Functional genomics has shown that non-coding RNAs regulate a diversity of developmental reproductive stages, from carpel formation and ovary development to the softening of the ripe/ripened fruit. This layer of transcriptional control has been associated with ovule, seed, and fruit development as well as fruit ripening, which are crucial developmental processes in breeding programs because of their relevance for crop production. The final ripe fruit is the result of a process under multiple levels of regulation, including mechanisms orchestrated by microRNAs and lncRNAs. Most of the studies we discuss involve work on tomato and Arabidopsis. In this review, we summarize non-coding RNA-controlled mechanisms described in the current literature that act coordinating the main steps of gynoecium development/patterning and fruit ripening.

Tuning of SlARF10A dosage by sly-miR160a is critical for auxin-mediated compound leaf and flower development

miR160 adjusts auxin-mediated development by post-transcriptional regulation of the auxin response factors ARF10/16/17. In tomato, knockdown of miR160 (sly-miR160) suggested that it is required for auxin-driven leaf blade outgrowth, but whether additional developmental events are adjusted by sly-miR160 is not clear. Here, the SlMIR160 genes and the genes of its SlARFs targets were edited by CRISPR/Cas9 resulting in the isolation of loss-of-function mutants. In addition, hypomorphic mutants that accumulate variable reduced levels of sly-miR160a were isolated. We found that the loss-of-function mutants in SlMIR160a (CR-slmir160a-6/7) produced only four wiry leaves, whereas the hypomorphic mutants developed leaves and flowers with graded developmental abnormalities. Phenotypic severity correlated with the upregulation of SlARF10A. Consistent with that, double mutants in SlMIR160a and SlARF10A restored leaf and flower development indicating that over-accumulation of SlARF10A underlay the developmental abnormalities exhibited in the CR-slmir160a mutants. Phenotype severity also correlated with the upregulation of the SHOOT MERISTEMLESS homolog Tomato Knotted 2, which in turn activated the transcription of the cytokinin biosynthesis genes SlIPT2 and SlIPT4. However, no change in Tomato Knotted 2 was detected in the absence of SlARF10A, suggesting that it is upregulated due to auxin signaling suppression by SlARF10A. Knockout of sly-miR160a-targeted SlARFs showed that whereas SlARF10A is indispensable for leaf blade outgrowth and floral organ patterning, the functions of SlARF16A and SlARF17 are redundant. Taken together our results suggest that sly-miR160a promotes blade outgrowth as well as leaf and leaflet initiation and floral organ development through the quantitative regulation of its major target SlARF10A.

Integration of small RNAs and transcriptome sequencing uncovers a complex regulatory network during vernalization and heading stages of orchardgrass (Dactylis glomerata L.)

BACKGROUND

Flowering is a critical reproductive process in higher plants. Timing of optimal flowering depends upon the coordination among seasonal environmental cues. For cool season grasses, such as Dactylis glomerata, vernalization induced by low temperature provides competence to initiate flowering after prolonged cold. We combined analyses of the transcriptome and microRNAs (miRNAs) to generate a comprehensive resource for regulatory miRNAs and their target circuits during vernalization and heading stages.

RESULTS

A total of 3,846 differentially expressed genes (DEGs) and 69 differentially expressed miRNAs were identified across five flowering stages. The expression of miR395, miR530, miR167, miR396, miR528, novel_42, novel_72, novel_107, and novel_123 demonstrated significant variations during vernalization. These miRNA targeted genes were involved in phytohormones, transmembrane transport, and plant morphogenesis in response to vernalization. The expression patterns of DEGs related to plant hormones, stress responses, energy metabolism, and signal transduction changed significantly in the transition from vegetative to reproductive phases.

CONCLUSIONS

Five hub genes, c136110_g1 (BRI1), c131375_g1 (BZR1), c133350_g1 (VRN1), c139830_g1 (VIN3), and c125792_g2 (FT), might play central roles in vernalization response. Our comprehensive analyses have provided a useful platform for investigating consecutive transcriptional and post-transcriptional regulation of critical phases in D. glomerata and provided insights into the genetic engineering of flowering-control in cereal crops.

Systematic identification of long non-coding RNAs during pollen development and fertilization in Brassica rapa

The importance of long non-coding RNAs (lncRNAs) in plant development has been established, but a systematic analysis of lncRNAs expressed during pollen development and fertilization has been elusive. We performed a time series of RNA-seq experiments at five developmental stages during pollen development and three different time points after pollination in Brassica rapa and identified 12 051 putative lncRNAs. A comprehensive view of dynamic lncRNA expression networks underpinning pollen development and fertilization was provided. B. rapa lncRNAs share many common characteristics of lncRNAs: relatively short length, low expression but specific in narrow time windows, and low evolutionary conservation. Gene modules and key lncRNAs regulating reproductive development such as exine formation were uncovered. Forty-seven cis-acting lncRNAs and 451 trans-acting lncRNAs were revealed to be highly coexpressed with their target protein-coding genes. Of particular importance are the discoveries of 14 lncRNAs that were highly coexpressed with 10 function-known pollen-associated coding genes. Fifteen lncRNAs were predicted as endogenous target mimics for 13 miRNAs, and two lncRNAs were proved to be functional target mimics for miR160 after experimental verification and shown to function in pollen development. Our study provides the systematic identification of lncRNAs during pollen development and fertilization in B. rapa and forms the foundation for future genetic, genomic, and evolutionary studies.

Systematic miRNome profiling reveals differential microRNAs in transgenic maize metabolism

BACKGROUND

While some genetically modified organisms (GMOs) are created to produce new double-stranded RNA molecules (dsRNA), in others, such molecules may occur as an unintended effect of the genetic engineering process. Furthermore, GMOs might produce naturally occurring dsRNA molecules in higher or lower quantities than its non-transgenic counterpart. This study is the first to use high-throughput technology to characterize the miRNome of commercialized GM maize events and to investigate potential alterations in miRNA regulatory networks.

RESULTS

Thirteen different conserved miRNAs were found to be dys-regulated in GM samples. The insecticide Bt GM variety had the most distinct miRNome. These miRNAs target a range of endogenous transcripts, such as transcription factors and nucleic acid binding domains, which play key molecular functions in basic genetic regulation. In addition, we have identified 20 potential novel miRNAs with target transcripts involved in lipid metabolism in maize. isomiRs were also found in 96 conserved miRNAs sequences, as well as potential transgenic miRNA sequences, which both can be a source of potential off-target effects in the plant genome. We have also provided information on technical limitations and when to carry on additional in vivo experimental testing.

CONCLUSIONS

These findings do not reveal hazards per se but show that robust and reproducible miRNA profiling technique can strengthen the assessment of risk by detecting any new intended and unintended dsRNA molecules, regardless of the outcome, at any stage of GMO development.

Identification and characterization of differentially expressed Phaseolus vulgaris miRNAs and their targets during mungbean yellow mosaic India virus infection reveals new insight into Phaseolus-MYMIV interaction

Phaseolus vulgaris is an economically important legume in tropical and subtropical regions of Asia, Africa, Latin-America and parts of USA and Europe. However, its production gets severely affected by mungbean yellow mosaic India virus (MYMIV). We aim to identify and characterize differentially expressed miRNAs during MYMIV-infection in P. vulgaris. A total of 422 miRNAs are identified of which 292 are expressed in both MYMIV-treated and mock-treated samples, 109 are expressed only in MYMIV-treated and 21 are expressed only in mock-treated samples. Selected up- and down-regulated miRNAs are validated by RT-qPCR. 3367 target ORFs are identified for 270 miRNAs. Selected targets are validated by 5' RLM-RACE. Differentially expressed miRNAs regulate transcription factors and are involved in improving stress tolerance to MYMIV. These findings will provide an insight into the role of miRNAs during MYMIV infection in P. vulgaris in particular and during any biotic stress conditions in Leguminosae family in general.

Small RNA and degradome sequencing used to elucidate the basis of tolerance to salinity and alkalinity in wheat

BACKGROUND

Soil salinity and/or alkalinity impose a major constraint over crop yield and quality. An understanding of the molecular basis of the plant response to these stresses could inform the breeding of more tolerant varieties. The bread wheat cultivar SR3 exhibits an enhanced level of salinity tolerance, while SR4 is distinguished by its superior tolerance of alkalinity.

RESULTS

The small RNA and degradome sequencing was used to explore the miRNAs and corresponding targets associated with the superior stress tolerance of the SR lines. An examination of the small RNA content of these two closely related lines revealed the presence of 98 known and 219 novel miRNA sequences. Degradome libraries were constructed in order to identify the targets of the miRNAs, leading to the identification of 58 genes targeted by 26 of the known miRNAs and 549 targeted by 65 of the novel ones. The function of two of the stress-responsive miRNAs was explored using virus-induced gene silencing.

CONCLUSIONS

This analysis indicated that regulation mediated by both auxin and epigenetic modification can be important in determining both salinity and alkalinity tolerance, while jasmonate signaling and carbohydrate metabolism are important for salinity tolerance, as is proton transport for alkalinity tolerance.

Molecular improvement of alfalfa for enhanced productivity and adaptability in a changing environment

Due to an expanding world population and increased buying power, the demand for ruminant products such as meat and milk is expected to grow substantially in coming years, and high levels of forage crop production will therefore be a necessity. Unfortunately, urbanization of agricultural land, intensive agricultural practices, and climate change are all predicted to limit crop production in the future, which means that the development of forage cultivars with improved productivity and adaptability will be essential. Because alfalfa (Medicago sativa L.) is one of the most widely cultivated perennial forage crops, it has been the target of much research in this field. In this review, we discuss progress that has been made towards the improvement of productivity, abiotic stress tolerance, and nutrient-use efficiency, as well as disease and pest resistance, in alfalfa using biotechnological techniques. Furthermore, we consider possible future priorities and avenues for attaining further enhancements in this crop as a means of contributing to the realization of food security in a changing environment.

Plant small RNAs: advancement in the understanding of biogenesis and role in plant development

Present review addresses the advances made in the understanding of biogenesis of plant small RNAs and their role in plant development. We discuss the elaborate role of microRNAs (miRNAs) and trans-acting small interfering RNAs (ta-siRNAs) in various aspects of plant growth and development and highlight relevance of small RNA mobility. Small non-coding RNAs regulate various aspects of plant development. Small RNAs (sRNAs) of 21-24 nucleotide length are derived from double-stranded RNAs through the combined activity of several biogenesis and processing components. These sRNAs function by negatively regulating the expression of target genes. miRNAs and ta-siRNAs constitute two important classes of endogenous small RNAs in plants, which play important roles in plant growth and developmental processes like embryogenesis, organ formation and patterning, shoot and root growth, and reproductive development. Biogenesis of miRNAs is a multistep process which includes transcription, processing and modification, and their loading onto RNA-induced silencing complex (RISC). RISC-loaded miRNAs carry out post-transcriptional silencing of their target(s). Recent studies identified orthologues of different biogenesis components of novel and conserved small RNAs from different model plants. Although many small RNAs have been identified from diverse plant species, only a handful of them have been functionally characterized. In this review, we discuss the advances made in understanding the biogenesis, functional conservation/divergence in miRNA-mediated gene regulation, and the developmental role of small RNAs in different plant species.

MicroRNAs regulate the main events in rice drought stress response by manipulating the water supply to shoots

MicroRNAs (miRNAs) are small endogenous regulatory RNAs that are involved in a variety of biological processes related to proliferation, development, and response to biotic and abiotic stresses. miRNA profiles of rice (Oryza sativa L. cv. IR64.) leaves in a partial root zone drying (PRD) system were analysed using a high-throughput sequencing approach to identify miRNAs associated with drought signalling. The treatments performed in this study were as follows: well-watered ("wet" roots, WW), wherein both halves of the pot were watered daily; drought ("dry" roots, DD), wherein water was withheld from both halves of the pot; and well-watered/drought ("wet" and "dry" roots, WD), wherein one half of each pot was watered daily, the same as in WW, and water was withheld from the other part, the same as in DD. High-throughput sequencing enabled us to detect novel miRNAs and study the differential expression of known miRNAs. A total of 209 novel miRNAs were detected in this study. Differential miRNA profiling of the DD, WD and WW conditions showed differential expression of 159 miRNAs, among which 83, 44 and 32 miRNAs showed differential expression under both DD and WD conditions. The detection of putative targets of the differentially expressed miRNAs and investigation of their functions showed that most of these genes encode transcription factors involved in growth and development, leaf morphology, regulation of hormonal homeostasis, and stress response. The most important differences between the DD and WD conditions involved regulation of the levels of hormones such as auxin, cytokinin, abscisic acid, and jasmonic acid and also regulation of phosphor homeostasis. Overall, differentially expressed miRNAs under WD conditions were found to differ from those under DD conditions, with such differences playing a role in adaptation and inducing the normal condition. The mechanisms involved in regulating hormonal homeostasis and involved in energy production and consumption were found to be the most important regulatory pathways distinguishing the DD and WD conditions.