Antiquity of microRNAs and their targets in land plants

Hairpin in a haystack: In silico identification and characterization of plant-conserved microRNA in Rafflesiaceae

Rafflesiaceae is a family of endangered plants whose members are solely parasitic to the tropical grape vine Tetrastigma (Vitaceae). Currently, the genetics of their crosstalk with the host remains unexplored. In this study, we use homology-based in silico approaches to characterize micro-RNAs (miRNAs) expressed by Sapria himalayana and Rafflesia cantleyi from published omics data. Derived from secondary structures or hairpins, miRNAs are small regulators of gene expression. We found that some plant-conserved miRNA still exists in Rafflesiaceae. Out of 9 highly conserved miRNA families in plants, 7 families (156/157, 159/319, 160, 165/166, 171, 172, 390) were identified with a total of 22 variants across Rafflesiaceae. Some miRNAs were missing endogenous targets and may have evolved to target host miRNA, though this requires experimental verification. Rafflesiaceae miRNA promoters are mostly inducible by ethylene that mediates stress response in the host but could be perceived by the parasites as a signal for growth. This study provides evidence that certain miRNAs with ancient origins in land plants still exist in Rafflesiaceae, though some may have been coopted by parasites to target host genes.

Differential Stress Responses to Rice Blast Fungal Infection Associated with the Vegetative Growth Phase in Rice

During vegetative growth, plants undergo various morphological and physiological changes in the transition from the juvenile phase to the adult phase. In terms of stress resistance, it has been suggested that plants gain or reinforce disease resistance during the process of maturation, which is recognized as adult plant resistance or age-related resistance. While much knowledge has been obtained about changes in disease resistance as growth stages progress, knowledge about changes in plant responses to pathogens with progressing age in plants is limited. In this study, we experimentally compared rice blast resistance in rice leaves sampled from plants at different growth phases. The results indicate differential infection progression and fungal status depending on growth stage. Transcriptome analysis following blast fungus infection revealed that several genes involved in the defense response were upregulated in both the juvenile and intermediate stage, but the expression changes of many genes were growth phase-specific. These findings highlight differences in rice leaf stress responses to blast infection at different growth stages.

microRNA858 represses the transcription factor gene SbMYB47 and regulates flavonoid biosynthesis in Scutellaria baicalensis

MicroRNAs (miRNAs) are noncoding endogenous single-stranded RNAs that regulate target gene expression by reducing their transcription and translation. Several miRNAs in plants function in secondary metabolism. The dried root of Scutellaria baicalensis Georgi is a traditional Chinese medicine that contains flavonoids (baicalin, wogonoside, and baicalein) as its main active ingredients. Although the S. baicalensis genome sequence has been published, information regarding its miRNAs is lacking. In this study, 12 small RNA libraries of different S. baicalensis tissues were compiled, including roots, stems, leaves, and flowers. A total of 129 miRNAs were identified, including 99 miRNAs from 27 miRNA families and 30 predicted miRNAs. Furthermore, 46 reliable target genes of 15 miRNA families were revealed using psRNATarget and confirmed by degradome sequencing. It was speculated that the microRNA858 (miR858)-SbMYB47 module might be involved in flavonoid biosynthesis. Transient assays in Nicotiana benthamiana leaves indicated that miR858 targets SbMYB47 and suppresses its expression. Artificial miRNA-mediated knockdown of miR858 and overexpression of SbMYB47 significantly increased the flavonoid content in S. baicalensis hairy roots, while SbMYB47 knockdown inhibited flavonoid accumulation. Yeast one-hybrid and dual-luciferase assays indicated that SbMYB47 directly binds to and activates the S. baicalensis phenylalanine ammonia-lyase 3 (SbPAL-3) and flavone synthase II (SbFNSⅡ-2) promoters. Our findings reveal the link between the miR858-SbMYB47 module and flavonoid biosynthesis, providing a potential strategy for the production of flavonoids with important pharmacological activities through metabolic engineering.

Advantage looping: Gene regulatory circuits between microRNAs and their target transcription factors in plants

The 20 to 24 nucleotide microRNAs (miRNAs) and their target transcription factors (TF) have emerged as key regulators of diverse processes in plants, including organ development and environmental resilience. In several instances, the mature miRNAs degrade the TF-encoding transcripts, while their protein products in turn bind to the promoters of the respective miRNA-encoding genes and regulate their expression, thus forming feedback loops (FBLs) or feedforward loops (FFLs). Computational analysis suggested that such miRNA-TF loops are recurrent motifs in gene regulatory networks (GRNs) in plants as well as animals. In recent years, modeling and experimental studies have suggested that plant miRNA-TF loops in GRNs play critical roles in driving organ development and abiotic stress responses. Here, we discuss the miRNA-TF FBLs and FFLs that have been identified and studied in plants over the past decade. We then provide some insights into the possible roles of such motifs within GRNs. Lastly, we provide perspectives on future directions for dissecting the functions of miRNA-centric GRNs in plants.

INDETERMINATE DOMAIN Transcription Factors in Crops: Plant Architecture, Disease Resistance, Stress Response, Flowering, and More

INDETERMINATE DOMAIN (IDD) genes encode plant-specific transcription factors containing a conserved IDD domain with four zinc finger motifs. Previous studies on Arabidopsis IDDs (AtIDDs) have demonstrated that these genes play roles in diverse physiological and developmental processes, including plant architecture, seed and root development, flowering, stress responses, and hormone signaling. Recent studies have revealed important functions of IDDs from rice and maize, especially in regulating leaf differentiation, which is related to the evolution of C4 leaves from C3 leaves. Moreover, IDDs in crops are involved in the regulation of agriculturally important traits, including disease and stress resistance, seed development, and flowering. Thus, IDDs are valuable targets for breeding manipulation. This review explores the role of IDDs in plant development, environmental responses, and evolution, which provides idea for agricultural application.

Differential expression of microRNAs in response to Papaya ringspot virus infection in differentially responding genotypes of papaya (Carica papaya L.) and its wild relative

Papaya ringspot virus (PRSV) is one of the most devastating viruses of papaya that has significantly hampered papaya production across the globe. Although PRSV resistance is known in some of its wild relatives, such as Vasconcellea cauliflora and in some of the improved papaya genotypes, the molecular basis of this resistance mechanism has not been studied and understood. Plant microRNAs are an important class of small RNAs that regulate the gene expression in several plant species against the invading plant pathogens. These miRNAs are known to manifest the expression of genes involved in resistance against plant pathogens, through modulation of the plant's biochemistry and physiology. In this study we made an attempt to study the overall expression pattern of small RNAs and more specifically the miRNAs in different papaya genotypes from India, that exhibit varying levels of tolerance or resistance to PRSV. Our study found that the expression of some of the miRNAs was differentially regulated in these papaya genotypes and they had entirely different miRNA expression profile in healthy and PRSV infected symptomatic plants. This data may help in improvement of papaya cultivars for resistance against PRSV through new breeding initiatives or biotechnological approaches such as genome editing.

Review: Plant microRNAs in pathogen defense: A panacea or a piece of the puzzle?

Plant microRNAs (miRNAs) control key agronomic traits that are associated with their conserved role(s) in development. However, despite a multitude of studies, the utility of miRNAs in plant-pathogen resistance remains less certain. Reviewing the literature identifies three general classes of miRNAs regarding plant pathogen defense. Firstly, a number of evolutionary dynamic 22 nucleotide miRNA families that repress large numbers of plant immunity genes, either directly, or through triggering the biogenesis of secondary siRNAs. However, understanding of their role in defense and of their manipulation to enhance pathogen resistance are still lacking. Secondly, highly conserved miRNAs that indirectly impact disease resistance through their targets that are primarily regulating development or hormone signaling. Any alteration of these miRNAs usually results in pleiotropic impacts, which may alter disease resistance in some plant species, and against some pathogens. Thirdly, are the comparatively diverse and evolutionary dynamic set of non-conserved miRNAs, some of which contribute to pathogen resistance, but whose narrow evolutionary presence will likely restrict their utility. Therefore, reflecting the diverse and evolving nature of plant-pathogen interactions, a complex interplay of plant miRNAs with pathogen responses exists. Any miRNA-based solution for pathogen resistance will likely be highly specific, rather than a general panacea.

Coordination of floral and fiber development in cotton (Gossypium) by hormone- and flavonoid-signalling associated regulatory miRNAs

Various plant development activities and stress responses are tightly regulated by various microRNAs (miRNA) and their target genes, or transcription factors in a spatiotemporal manner. Here, to exemplify how flowering-associated regulatory miRNAs synchronize their expression dynamics during floral and fiber development in cotton, constitutive expression diminution transgenic lines of auxin-signaling regulatory Gh-miR167 (35S-MIM167) were developed through target mimicry approach. 'Moderate' (58% to 80%)- and 'high' (> 80%)-Gh-miR167 diminution mimic lines showed dosage-dependent developmental deformities in anther development, pollen maturation, and fruit (= boll) formation. Cross pollination of 'moderate' 35S-MIM167 mimic lines with wild type (WT) plant partially restored boll formation and emergence of fiber initials on the ovule surface. Gh-miR167 diminution favored organ-specific transcription biases in miR159, miR166 as well as miR160, miR164, and miR172 along with their target genes during anther and petal development, respectively. Similarly, accumulative effect of percent Gh-miR167 diminution, cross regulation of its target ARF6/8 genes, and temporal mis-expression of hormone signaling- and flavonoid biosynthesis-associated regulatory miRNAs at early fiber initiation stage caused irregular fiber formation. Spatial and temporal transcription proportions of regulatory miRNAs were also found crucial for the execution of hormone- and flavonoid-dependent progression of floral and fiber development. These observations discover how assorted regulatory genetic circuits get organized in response to Gh-miR167 diminution and converge upon ensuing episodes of floral and fiber development in cotton.

Identification of conserved miRNAs and their targets in Jatropha curcas: an in silico approach

BACKGROUND

MicroRNAs (miRNAs) are small endogenous RNAs with an approximate length of 18-22 nucleotides and involved in the regulation of gene expression in transcriptional or post-transcriptional levels. They were found to be associated with leaf morphogenesis, flowering time, vegetative phase change, and response to environmental cues in plants, where they act as a critical regulatory factor. The nature of high conservancy of plant miRNAs within the plant species made it possible to detect the conserved miRNAs by computational approaches. Expressed Sequence Tags (EST) based comparative genomic approaches provide advantages over wet lab approaches as it is convenient, easy to carry out and less time consuming. EST-based in silico approach can unravel new conserved miRNAs in plants, even when the complete genome sequence is not available.

RESULTS

To identify the novel miRNAs, a total of 46,865 ESTs from Jatropha curcas were searched for homology to all available 6746 mature miRNAs of plant eudicotyledons. Finally, we ended up with 12 novel miRNAs in Jatropha that range from 18 to 19 nucleotides where their respective precursor miRNAs had 54.11-71.76% (A + U) content. The putative miRNAs belong to 12 individual miRNA family and most of them have higher (A + U) content ranging from 47.36 to 77.77% than their respective miRNA homologs. Many of the target genes by the newly identified miRNAs were associated with plant growth and development, stress response, defense and hormone signaling, and oil synthesis pathways.

CONCLUSION

These findings have the potential to speed up miRNA identification and expand our understanding of miRNA functions in J. curcas.

MicroRNA162 regulates stomatal conductance in response to low night temperature stress via abscisic acid signaling pathway in tomato

MicroRNAs (miRNAs) mediate the degradation of target mRNA and inhibit mRNA translation to regulate gene expression at the transcriptional and post-transcriptional levels in response to environmental stress in plants. We characterized the post-transcriptional mechanism by deep sequencing small RNA (sRNA) to examine how miRNAs were involved in low night temperature (LNT) stress in tomato and whether the molecular mechanism depended on the abscisic acid (ABA) signaling pathway. We annotated conserved miRNAs and novel miRNAs with four sRNA libraries composed of wild-type (WT) tomato plants and ABA-deficient mutant (sit) plants under normal growth and LNT stress conditions. Reverse genetics analysis suggested that miR162 participated in LNT resistance and the ABA-dependent signaling pathway in tomato. miR162-overexpressing (pRI-miR162) and miR162-silenced (pRNAi-miR162) transgenic tomato plants were generated to evaluate miR162 functions in response to LNT stress. miR162 deficiency exhibited high photosynthetic capacity and regulated stomatal opening, suggesting negative regulation of miR162 in the ABA-dependent signaling pathway in response to LNT stress. As feedback regulation, miR162 positively regulated ABA to maintain homeostasis of tomato under diverse abiotic stresses. The mRNA of DICER-LIKE1 (DCL1) was targeted by miR162, and miR162 inhibited DCL1 cleavage in LNT response, including the regulation of miRNA160/164/171a and their targets. The DCL1-deficient mutants (dcl1) with CRISPR/Cas9 prevented stomatal opening to influence photosynthesis in the ABA signaling pathway under LNT stress. Finally, we established the regulatory mechanism of ABA-miR162-DCL1, which systematically mediated cold tolerance in tomato. This study suggests that post-transcriptional modulators acted as systemic signal responders via the stress hormone signaling pathway, and the model at the post-transcriptional level presents a new direction for research in plant abiotic stress resistance.

Oligo-Not Only for Silencing: Overlooked Potential for Multidirectional Action in Plants

Oligo technology is a low-cost and easy-to-implement method for direct manipulation of gene activity. The major advantage of this method is that gene expression can be changed without requiring stable transformation. Oligo technology is mainly used for animal cells. However, the use of oligos in plants seems to be even easier. The oligo effect could be similar to that induced by endogenous miRNAs. In general, the action of exogenously introduced nucleic acids (Oligo) can be divided into a direct interaction with nucleic acids (genomic DNA, hnRNA, transcript) and an indirect interaction via the induction of processes regulating gene expression (at the transcriptional and translational levels) involving regulatory proteins using endogenous cellular mechanisms. Presumed mechanisms of oligonucleotides' action in plant cells (including differences from animal cells) are described in this review. Basic principles of oligo action in plants that allow bidirectional changes in gene activity and even those that lead to heritable epigenetic changes in gene expression are presented. The effect of oligos is related to the target sequence at which they are directed. This paper also compares different delivery methods and provides a quick guide to using IT tools to help design oligonucleotides.

MIR390 Is Involved in Regulating Anthracnose Resistance in Apple

As an important cash crop in China, apple has a good flavor and is rich in nutrients. Fungal attacks have become a major obstacle in apple cultivation. Colletotrichum gloeosporioides is one of the most devastating fungal pathogens in apple. Thus, discovering resistance genes in response to C. gloeosporioides may aid in designing safer control strategies and facilitate the development of apple resistance breeding. A previous study reported that 'Hanfu' autotetraploid apple displayed higher C. gloeosporioides resistance than 'Hanfu' apple, and the expression level of mdm-MIR390b was significantly upregulated in autotetraploid plants compared to that in 'Hanfu' plants, as demonstrated by digital gene expression (DGE) analysis. It is still unclear, however, whether mdm-MIR390b regulates apple anthracnose resistance. Apple MIR390b was transformed into apple 'GL-3' plants to identify the functions of mdm-MIR390b in anthracnose resistance. C. gloeosporioides treatment analysis indicated that the overexpression of mdm-MIR390b reduced fungal damage to apple leaves and fruit. Physiology analysis showed that mdm-MIR390b increased C. gloeosporioides resistance by improving superoxide dismutase (SOD) and peroxidase (POD) activity to alleviate the damage caused by O₂^- and H₂O₂. Our results demonstrate that mdm-MIR390b can improve apple plants' anthracnose resistance.

Integrated High-Throughput Sequencing, Microarray Hybridization and Degradome Analysis Uncovers MicroRNA-Mediated Resistance Responses of Maize to Pathogen Curvularia lunata

Curvularia lunata (Wakker) Boed, the causal agent of leaf spot in maize, is prone to mutation, making it difficult to control. RNAi technology has proven to be an important tool of genetic engineering and functional genomics aimed for crop improvement. MicroRNAs (miRNAs), which act as post-transcriptional regulators, often cause translational repression and gene silencing. In this article, four small RNA (sRNA) libraries were generated from two maize genotypes inoculated by C. lunata; among these, ltR1 and ltR2 were from the susceptible variety Huangzao 4 (HZ), ltR3 and ltR4, from the resistant variety Luyuan (LY), and 2286, 2145, 1556 and 2504 reads were annotated as miRNA in these four sRNA libraries, respectively. Through the combined analysis of high-throughput sequencing, microarray hybridization and degradome, 48 miRNAs were identified as being related to maize resistance to C. lunata. Among these, PC-732 and PC-169, two new maize miRNAs discovered, were predicted to cleave mRNAs of metacaspase 1 (AMC1) and thioredoxin family protein (Trx), respectively, possibly playing crucial roles in the resistance of maize to C. lunata. To further confirm the role of PC-732 in the interaction of maize and C. lunata, the miRNA was silenced through STTM (short tandem target mimic) technology, and we found that knocking down PC-732 decreased the susceptibility of maize to C. lunata. Precisely speaking, the target gene of PC-732 might inhibit the expression of disease resistance-related genes during the interaction between maize and C. lunata. Overall, the findings of this study indicated the existence of miRNAs involved in the resistance of maize to C. lunata and will contribute to rapidly clarify the resistant mechanism of maize to C. lunata.

Gma-miR408 Enhances Soybean Cyst Nematode Susceptibility by Suppressing Reactive Oxygen Species Accumulation

Soybean cyst nematode (SCN, Heterodera glycine) is a serious damaging disease in soybean worldwide, thus resulting in severe yield losses. MicroRNA408 (miR408) is an ancient and highly conserved miRNA involved in regulating plant growth, development, biotic and abiotic stress response. Here, we analyzed the evolution of miR408 in plants and verified four miR408 members in Glycine max. In the current research, highly upregulated gma-miR408 expressing was detected during nematode migration and syncytium formation response to soybean cyst nematode infection. Overexpressing and silencing miR408 vectors were transformed to soybean to confirm its potential role in plant and nematode interaction. Significant variations were observed in the MAPK signaling pathway with low OXI1, PR1, and wounding of the overexpressing lines. Overexpressing miR408 could negatively regulate soybean resistance to SCN by suppressing reactive oxygen species accumulation. Conversely, silencing miR408 positively regulates soybean resistance to SCN. Overall, gma-miR408 enhances soybean cyst nematode susceptibility by suppressing reactive oxygen species accumulation.

Analysis of Homologous Regions of Small RNAs MIR397 and MIR408 Reveals the Conservation of Microsynteny among Rice Crop-Wild Relatives

MIRNAs are small non-coding RNAs that play important roles in a wide range of biological processes in plant growth and development. MIR397 (involved in drought, low temperature, and nitrogen and copper (Cu) starvation) and MIR408 (differentially expressed in response to environmental stresses such as copper, light, mechanical stress, dehydration, cold, reactive oxygen species, and drought) belong to conserved MIRNA families that either negatively or positively regulate their target genes. In the present study, we identified the homologs of MIR397 and MIR408 in Oryza sativa and its six wild progenitors, three non-Oryza species, and one dicot species. We analyzed the 100 kb segments harboring MIRNA homologs from 11 genomes to obtain a comprehensive view of their community evolution around these loci in the farthest (distant) relatives of rice. Our study showed that mature MIR397 and MIR408 were highly conserved among all Oryza species. Comparative genomics analyses also revealed that the microsynteny of the 100 kb region surrounding MIRNAs was only conserved in Oryza spp.; disrupted in Sorghum, maize, and wheat; and completely lost in Arabidopsis. There were deletions, rearrangements, and translocations within the 100 kb segments in Oryza spp., but the overall microsynteny of the region was maintained. The phylogenetic analyses of the precursor regions of all MIRNAs under study revealed a bimodal clade of common origin. This comparative analysis of miRNA involved in abiotic stress tolerance in plants provides a powerful tool for future Oryza research. Crop wild relatives (CWRs) offer multiple traits with potential to decrease the amount of yield loss owing to biotic and abiotic stresses. Using a comparative genomics approach, the exploration of CWRs as a source of tolerance to these stresses by understanding their evolution can be further used to leverage their yield potential.

Identification of microRNAs That Provide a Low Light Stress Tolerance-Mediated Signaling Pathway during Vegetative Growth in Rice

Low light intensity affects several physiological parameters during the different growth stages in rice. Plants have various regulatory mechanisms to cope with stresses. One of them is the differential and temporal expression of genes, which is governed by post-transcriptional gene expression regulation through endogenous miRNAs. To decipher low light stress-responsive miRNAs in rice, miRNA expression profiling was carried out using next-generation sequencing of low-light-tolerant (Swarnaprabha) and -sensitive (IR8) rice genotypes through Illumina sequencing. Swarnaprabha and IR8 were subjected to 25% low light treatment for one day, three days, and five days at the active tillering stage. More than 43 million raw reads and 9 million clean reads were identified in Swarnaprabha, while more than 41 million raw reads and 8.5 million clean reads were identified in IR8 after NGS. Importantly, 513 new miRNAs in rice were identified, whose targets were mostly regulated by the genes involved in photosynthesis and metabolic pathways. Additionally, 114 known miRNAs were also identified. Five novel (osa-novmiR1, osa-novmiR2, osa-novmiR3, osa-novmiR4, and osa-novmiR5) and three known (osa-miR166c-3p, osa-miR2102-3p, and osa-miR530-3p) miRNAs were selected for their expression validation through miRNA-specific qRT-PCR. The expression analyses of most of the predicted targets of corresponding miRNAs show negative regulation. Hence, miRNAs modulated the expression of genes providing tolerance/susceptibility to low light stress. This information might be useful in the improvement of crop productivity under low light stress.

Evolutionary conservation of nested MIR159 structural microRNA genes and their promoter characterization in Arabidopsis thaliana

MicroRNAs (miRNAs) are endogenous small RNAs, that are vital for gene expression regulation in eukaryotes. Whenever a pri-miRNA precursor includes another miRNA precursor, and both of these precursors may generate independent, non-overlapping mature miRNAs, we named them nested miRNAs. However, the extent of nested miR159 structural evolutionary conservation and its promoter characterization remains unknown. In this study, the sequence alignment and phylogenetic analysis reveal that the MIR159 family is ancient, and its nested miR159 structures are evolutionary conserved in different plant species. The overexpression of ath-MIR159a, including the 1.2 kb downstream region, has no effect on rescuing the mir159ab phenotype. The promoter truncation results revealed that the 1.0 kb promoter of ath-MIR159a is sufficient for rescuing the mir159ab phenotype. The cis-regulatory elements in the ath-miR159a promoters indicated functions related to different phytohormones, abiotic stresses, and transcriptional activation. While the MybSt1 motif-containing region is not responsible for activating the regulation of the miR159a promoter. The qRT-PCR results showed that overexpression of ath-MIR159a led to high expression levels of miR159a.1-5 and miR159a.1-3 and complemented the growth defect of mir159ab via downregulation of MYB33 and MYB65. Furthermore, continuously higher expression of the miR159a.2 duplex in transgenic lines with the curly leaf phenotype indicates that miR159a.2 is functional in Arabidopsis and suggests that it is possible for a miRNA precursor to encode several regulatory small RNAs in plants. Taken together, our study demonstrates that the nested miR159 structure is evolutionary conserved and miRNA-mediated gene regulation is more complex than previously thought.

Identification of genes differentially expressed between prostrate shoots and erect shoots in the lycophyte Selaginella nipponica using an RNA-seq approach

Lycophytes are the earliest vascular plants and Selaginella is the most studied genus among them. Prostrate shoots are produced during early growth and erect shoots emerge later in S. nipponica, thus providing an opportunity for exploring the evolution of the mechanism underlying the transition between growth phases. Six libraries were sequenced for the prostrate and the erect shoots, and a total of 206 768 genes were identified. Some genes were differentially expressed in prostate and erect shoot, with relatively high expression in the prostate shoots being related to hormone responses and defence reactions, while higher expression in the erect shoots was related to spore formation and shoot development. Some SPL genes possessed a miR156 binding site and were highly expressed in the erect shoots, while AP2-like genes were more highly expressed in the prostrate shoots but simultaneously lacked any miR172 binding site. MiR156 was detected at a higher concentration in the prostrate shoots. Thus, the mechanism for the vegetative to reproductive transition of sporophytes probably originated in the common ancestor of vascular plants and must have experienced stepwise development during evolution.

Dosage-sensitive miRNAs trigger modulation of gene expression during genomic imbalance in maize

The genomic imbalance caused by varying the dosage of individual chromosomes or chromosomal segments (aneuploidy) has more detrimental effects than altering the dosage of complete chromosome sets (ploidy). Previous analysis of maize (Zea mays) aneuploids revealed global modulation of gene expression both on the varied chromosome (cis) and the remainder of the genome (trans). However, little is known regarding the role of microRNAs (miRNAs) under genomic imbalance. Here, we report the impact of aneuploidy and polyploidy on the expression of miRNAs. In general, cis miRNAs in aneuploids present a predominant gene-dosage effect, whereas trans miRNAs trend toward the inverse level, although other types of responses including dosage compensation, increased effect, and decreased effect also occur. By contrast, polyploids show less differential miRNA expression than aneuploids. Significant correlations between expression levels of miRNAs and their targets are identified in aneuploids, indicating the regulatory role of miRNAs on gene expression triggered by genomic imbalance.

Regulatory role of microRNAs (miRNAs) in the recent development of abiotic stress tolerance of plants

MicroRNAs (miRNAs) are a distinct groups of single-stranded non-coding, tiny regulatory RNAs approximately 20-24 nucleotides in length. miRNAs negatively influence gene expression at the post-transcriptional level and have evolved considerably in the development of abiotic stress tolerance in a number of model plants and economically important crop species. The present review aims to deliver the information on miRNA-mediated regulation of the expression of major genes or Transcription Factors (TFs), as well as genetic and regulatory pathways. Also, the information on adaptive mechanisms involved in plant abiotic stress responses, prediction, and validation of targets, computational tools, and databases available for plant miRNAs, specifically focus on their exploration for engineering abiotic stress tolerance in plants. The regulatory function of miRNAs in plant growth, development, and abiotic stresses consider in this review, which uses high-throughput sequencing (HTS) technologies to generate large-scale libraries of small RNAs (sRNAs) for conventional screening of known and novel abiotic stress-responsive miRNAs adds complexity to regulatory networks in plants. The discoveries of miRNA-mediated tolerance to multiple abiotic stresses, including salinity, drought, cold, heat stress, nutritional deficiency, UV-radiation, oxidative stress, hypoxia, and heavy metal toxicity, are highlighted and discussed in this review.

Transcriptome and Small RNA Profiling of Potato Virus Y Infected Potato Cultivars, Including Systemically Infected Russet Burbank

Potatoes are the world’s most produced non-grain crops and an important food source for billions of people. Potatoes are susceptible to numerous pathogens that reduce yield, including Potato virus Y (PVY). Genetic resistance to PVY is a sustainable way to limit yield and quality losses due to PVY infection. Potato cultivars vary in their susceptibility to PVY and include susceptible varieties such as Russet Burbank, and resistant varieties such as Payette Russet. Although the loci and genes associated with PVY-resistance have been identified, the genes and mechanisms involved in limiting PVY during the development of systemic infections have yet to be fully elucidated. To increase our understanding of PVY infection, potato antiviral responses, and resistance, we utilized RNA sequencing to characterize the transcriptomes of two potato cultivars. Since transcriptional responses associated with the extreme resistance response occur soon after PVY contact, we analyzed the transcriptome and small RNA profile of both the PVY-resistant Payette Russet cultivar and PVY-susceptible Russet Burbank cultivar 24 h post-inoculation. While hundreds of genes, including terpene synthase and protein kinase encoding genes, exhibited increased expression, the majority, including numerous genes involved in plant pathogen interactions, were downregulated. To gain greater understanding of the transcriptional changes that occur during the development of systemic PVY-infection, we analyzed Russet Burbank leaf samples one week and four weeks post-inoculation and identified similarities and differences, including higher expression of genes involved in chloroplast function, photosynthesis, and secondary metabolite production, and lower expression of defense response genes at those time points. Small RNA sequencing identified different populations of 21- and 24-nucleotide RNAs and revealed that the miRNA profiles in PVY-infected Russet Burbank plants were similar to those observed in other PVY-tolerant cultivars and that during systemic infection ~32% of the NLR-type disease resistance genes were targeted by 21-nt small RNAs. Analysis of alternative splicing in PVY-infected potato plants identified splice variants of several hundred genes, including isoforms that were more dominant in PVY-infected plants. The description of the PVYN-Wi-associated transcriptome and small RNA profiles of two potato cultivars described herein adds to the body of knowledge regarding differential outcomes of infection for specific PVY strain and host cultivar pairs, which will help further understanding of the mechanisms governing genetic resistance and/or virus-limiting responses in potato plants.

MiRNA fine tuning for crop improvement: using advance computational models and biotechnological tools

MiRNAs modulate target genes expression at post-transcriptional levels, by reducing spatial abundance of mRNAs. MiRNAs regulats plant metabolism, and emerged as regulators of plant stress responses. Which make miRNAs promising candidates for fine tuning to affectively alter crop stress tolerance and other important traits. With recent advancements in the computational biology and biotechnology miRNAs structure and target prediction is possible resulting in pin point editing; miRNA modulation can be done by up or down regulating miRNAs using recently available biotechnological tools (CRISPR Cas9, TALENS and RNAi). In this review we have focused on miRNA biogenesis, miRNA roles in plant development, plant stress responses and roles in signaling pathways. Additionally we have discussed latest computational prediction models for miRNA to target gene interaction and biotechnological systems used recently for miRNA modulation. We have also highlighted setbacks and limitations in the way of miRNA modulation; providing entirely a new direction for improvement in plant genomics primarily focusing miRNAs.

ASRmiRNA: Abiotic Stress-Responsive miRNA Prediction in Plants by Using Machine Learning Algorithms with Pseudo K-Tuple Nucleotide Compositional Features

MicroRNAs (miRNAs) play a significant role in plant response to different abiotic stresses. Thus, identification of abiotic stress-responsive miRNAs holds immense importance in crop breeding programmes to develop cultivars resistant to abiotic stresses. In this study, we developed a machine learning-based computational method for prediction of miRNAs associated with abiotic stresses. Three types of datasets were used for prediction, i.e., miRNA, Pre-miRNA, and Pre-miRNA + miRNA. The pseudo K-tuple nucleotide compositional features were generated for each sequence to transform the sequence data into numeric feature vectors. Support vector machine (SVM) was employed for prediction. The area under receiver operating characteristics curve (auROC) of 70.21, 69.71, 77.94 and area under precision-recall curve (auPRC) of 69.96, 65.64, 77.32 percentages were obtained for miRNA, Pre-miRNA, and Pre-miRNA + miRNA datasets, respectively. Overall prediction accuracies for the independent test set were 62.33, 64.85, 69.21 percentages, respectively, for the three datasets. The SVM also achieved higher accuracy than other learning methods such as random forest, extreme gradient boosting, and adaptive boosting. To implement our method with ease, an online prediction server “ASRmiRNA” has been developed. The proposed approach is believed to supplement the existing effort for identification of abiotic stress-responsive miRNAs and Pre-miRNAs.

The Evolution and Functional Roles of miR408 and Its Targets in Plants

MicroRNA408 (miR408) is an ancient and highly conserved miRNA, which is involved in the regulation of plant growth, development and stress response. However, previous research results on the evolution and functional roles of miR408 and its targets are relatively scattered, and there is a lack of a systematic comparison and comprehensive summary of the detailed evolutionary pathways and regulatory mechanisms of miR408 and its targets in plants. Here, we analyzed the evolutionary pathway of miR408 in plants, and summarized the functions of miR408 and its targets in regulating plant growth and development and plant responses to various abiotic and biotic stresses. The evolutionary analysis shows that miR408 is an ancient and highly conserved microRNA, which is widely distributed in different plants. miR408 regulates the growth and development of different plants by down-regulating its targets, encoding blue copper (Cu) proteins, and by transporting Cu to plastocyanin (PC), which affects photosynthesis and ultimately promotes grain yield. In addition, miR408 improves tolerance to stress by down-regulating target genes and enhancing cellular antioxidants, thereby increasing the antioxidant capacity of plants. This review expands and promotes an in-depth understanding of the evolutionary and regulatory roles of miR408 and its targets in plants.

Viral suppressor of RNA silencing in vascular plants also interferes with the development of the bryophyte Physcomitrella patens

Plant viruses are important pathogens able to overcome plant defense mechanisms using their viral suppressors of RNA silencing (VSR). Small RNA pathways of bryophytes and vascular plants have significant similarities, but little is known about how viruses interact with mosses. This study elucidated the responses of Physcomitrella patens to two different VSRs. We transformed P. patens plants to express VSR P19 from tomato bushy stunt virus and VSR 2b from cucumber mosaic virus, respectively. RNA sequencing and quantitative PCR were used to detect the effects of VSRs on gene expression. Small RNA (sRNA) sequencing was used to estimate the influences of VSRs on the sRNA pool of P. patens. Expression of either VSR-encoding gene caused developmental disorders in P. patens. The transcripts of four different transcription factors (AP2/erf, EREB-11 and two MYBs) accumulated in the P19 lines. sRNA sequencing revealed that VSR P19 significantly changed the microRNA pool in P. patens. Our results suggest that VSR P19 is functional in P. patens and affects the abundance of specific microRNAs interfering with gene expression. The results open new opportunities for using Physcomitrella as an alternative system to study plant-virus interactions.

Plant RNA-mediated gene regulatory network

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

Exploring diverse roles of micro RNAs in banana: Current status and future prospective

Micro RNAs (miRNAs) are 20-24 nucleotides long non-coding RNA sequences identified and characterized in multiple plant and animal systems. miRNAs play multifarious roles ranging from plant development to stress tolerance by synchronizing physiological processes at the level of transcription and translation. Banana is a major horticultural crop with colossal production worldwide. Despite the recent encouraging developments, the information on functions of miRNAs in banana plants is still in its infancy. The available literature pertaining to miRNAs in banana plants hints towards their contribution as master regulators in crucial physiological processes for instance abiotic stress responses, pathogenic defence response, fruit ripening and so on. This review is focused on biogenesis of miRNAs, their identification and deciphering their respective roles in banana plants with special emphasis on abiotic stress responses, plant immune responses, fruit ripening and storage. Based on the prior reports, we identified a few miRNAs with prospective roles in stress tolerance and illustrated the potential applications of miRNAs in banana crop improvement utilizing recent biotechnological tools such as CRISPR cas9, RNAi and the nano particle based foliar spray of miRNAs. The review briefly explained the future directions in banana research with a special emphasis on miRNA regulatory networks and agronomic traits improvement. Finally, future domains in miRNA research in plants and their possible applications towards crop improvement in agriculture are described briefly.

The APETALA2 homolog CaFFN regulates flowering time in pepper

Flowering time is an important agronomic trait that contributes to fitness in plants. However, the genetic basis of flowering time has not been extensively studied in pepper. To understand the genetics underlying flowering time, we constructed an F₂ population by crossing a spontaneous early flowering mutant and a late-flowering pepper line. Using bulked segregant RNA-seq, a major locus controlling flowering time in this population was mapped to the end of chromosome 2. An APETALA2 (AP2) homolog (CaFFN) cosegregated with flowering time in 297 individuals of the F₂ population. A comparison between the parents revealed a naturally occurring rare SNP (SNP2T > C) that resulted in the loss of a start codon in CaFFN in the early flowering mutant. Transgenic Nicotiana benthamiana plants with high CaFFN expression exhibited a delay in flowering time and floral patterning defects. On the other hand, pepper plants with CaFFN silencing flowered early. Therefore, the CaFFN gene acts as a flowering repressor in pepper. CaFFN may function as a transcriptional activator to activate the expression of CaAGL15 and miR156e and as a transcriptional repressor to repress the expression of CaAG, CaAP1, CaSEP3, CaSOC1, and miR172b based on a qRT-PCR assay. Direct activation of CaAGL15 by CaFFN was detected using yeast one-hybrid and dual-luciferase reporter assays, consistent with the hypothesis that CaFFN regulates flowering time. Moreover, the CaFFN gene association analysis revealed a significant association with flowering time in a natural pepper population, indicating that the CaFFN gene has a broad effect on flowering time in pepper. Finally, the phylogeny, evolutionary expansion and expression patterns of CaFFN/AP2 homologs were analyzed to provide valuable insight into CaFFN. This study increases our understanding of the involvement of CaFFN in controlling flowering time in pepper, thus making CaFFN a target gene for breeding early maturing pepper.

Identification and expression analysis of miRNAs in germination and seedling growth of Tibetan hulless barley

Germination and seedling growth are crucial for plant development and agricultural production. While, the regulatory mechanisms during these processes in Tibetan hulless barley (Hordeum vulgare L. var. nudum) are not well understood. Given the regulatory roles of microRNAs (miRNAs) in crop plants and the irreplaceability of barley in the highland area of China, we herein presented a genome-wide survey of miRNAs to reveal a potential regulatory network in the early developmental stages of two Tibetan hulless barleys, from which a total of 156 miRNAs was identified including 35 known and 121 novel ones. Six of the identified novel miRNAs were further experimentally validated. According to the evolutionary analysis, miR156, miR166, miR168, and miR171 were conserved across Tibetan hulless barleys and eight other seed plants. Expression profiles of ten known miRNAs showed that they were involved in phytohormone signaling, carbohydrate and lipid metabolism, as well as juvenile-adult transition during barley development. Moreover, a total of 1280 genes targeted by 101 miRNAs were predicted from both barley libraries. Three genes (PLN03212, MATE eukaryotic, and GRAS) were validated via the RNA ligase-mediated 5'-rapid amplification of cDNA ends (RLM-5' RACE) to be the targets of hvu-miR159a, hvu-miR166a, and hvu-miR171-3p, respectively. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of putative targets, the most abundant pathways were related to "metabolism". These results revealed that miRNA-target pairs participating in the regulation of multigene expression and the embryonic development of Tibetan hulless barleys were controlled by complex mechanisms involving the concordant expression of different miRNAs and feedback loops among miRNAs as well as their targets. The study provides insight into the regulatory network of barley miRNAs for better understanding of miRNA functions during germination and seedling growth.

Global expression dynamics and miRNA evolution profile govern floral/fiber architecture in the modern cotton (Gossypium)

Majority of differentially expressed miRNAs with functional attributes have been recruited independently and parallelly during allopolyploidy followed by the millennia of human selection of both domesticated G. hirsutum and G. barbadense. The genus Gossypium is a marvelous evolutionary model for studying allopolyploidy and morpho-evolution of long-spinnable fibers from the ancestral wild-fuzz. Many genes, transcription factors, and notably, the regulatory miRNAs essentially govern such remarkable modern fiber phenotypes. To comprehend the impact of allopolyploidy on the evolutionary selection of transcriptional dynamicity of key miRNAs, comparative transcriptome profiling of vegetative and fiber tissues of domesticated diploid G. arboreum (A₂) and allopolyploid cotton species G. hirsutum (AD₁), and G. barbadense (AD₂) identified > 300 differentially expressed miRNAs (DEmiRs) within or between corresponding tissues of A₂, AD₁ and AD₂ species. Up to 49% and 32% DEmiRs were up- and down-regulated at fiber initiation stage of AD₁ and AD₂ species, respectively, whereas 50% and 18% DEmiRs were up- and down-regulated at fiber elongation stage of both the allopolyploid species. Interestingly, A-subgenome-specific DEmiRs exhibit expression dominance in the allopolyploid genetic backgrounds. Comparative spatio-temporal expression analyses of AD₁ and AD₂ species discovered that a majority of DEmiRs were recruited independently under millennia of human selection during domestication. Functional annotations of these DEmiRs revealed selection of associated molecular functions such as hormone-signaling, calcium-signaling and reactive oxygen species (ROS) signaling during fiber initiation and elongation. To validate the functional attributes of annotated DEmiRs, we demonstrated for the first time that the target-mimicry-based constitutive diminution of auxin-signaling associated miR167 directly affected the differentiation of floral and fiber tissues of transgenic cotton. These results strongly suggested that the evolutionarily favored DEmiRs including miR167 were involved in the transcriptional regulation of numerous genes during cotton evolution for enhanced fiber-associated agronomic traits.

The Evolution of euAPETALA2 Genes in Vascular Plants: From Plesiomorphic Roles in Sporangia to Acquired Functions in Ovules and Fruits

The field of evolutionary developmental biology can help address how morphological novelties evolve, a key question in evolutionary biology. In Arabidopsis thaliana, APETALA2 (AP2) plays a role in the development of key plant innovations including seeds, flowers, and fruits. AP2 belongs to the AP2/ETHYLENE RESPONSIVE ELEMENT BINDING FACTOR family which has members in all viridiplantae, making it one of the oldest and most diverse gene lineages. One key subclade, present across vascular plants is the euAPETALA2 (euAP2) clade, whose founding member is AP2. We reconstructed the evolution of the euAP2 gene lineage in vascular plants to better understand its impact on the morphological evolution of plants, identifying seven major duplication events. We also performed spatiotemporal expression analyses of euAP2/TOE3 genes focusing on less explored vascular plant lineages, including ferns, gymnosperms, early diverging angiosperms and early diverging eudicots. Altogether, our data suggest that euAP2 genes originally contributed to spore and sporangium development, and were subsequently recruited to ovule, fruit and floral organ development. Finally, euAP2 protein sequences are highly conserved; therefore, changes in the role of euAP2 homologs during development are most likely due to changes in regulatory regions.

Endophytic Microbiome Responses to Sulfur Availability in Beta vulgaris (L.)

Sulfur is an essential plant macronutrient, and its adequate supply allows an efficient root storage and sugar extractability in sugar beets (Beta vulgaris L.). In this study, we investigated the effect of changes in sulfur availability on the endophytic community structure of sugar beets. Plants were hydroponically grown in a complete nutrient solution (S-supplied), a nutrient solution without MgSO₄ (S-deprived), and a nutrient solution without MgSO₄ for six days and resupplied with 100 μM MgSO₄ for 48 h (S-resupplied). The sulfur status was monitored by inductively coupled plasma ICP–OES, and combustion analysis together with the evaluation of microRNA395 as a biomarker for sulfate status. Metabarcoding of the bacterial 16S rRNA gene was carried out in order to determine leaf endophytic community structure. The Shannon diversity index significantly differed (p < 0.05) between sulfate-supplied and sulfate-deprived seedlings. Validation by Real-Time PCR showed a significant increase (p < 0.05) of Burkholderia spp. in sulfate-deprived plants as compared to sulfate-supplied ones. The study sheds new light on the effects of nutrient deficiency on the microbiome of sugar beet plants.

Molecular Cloning, Transcriptional Profiling, Subcellular Localization, and miRNA-Binding Site Analysis of Six SCL9 Genes in Poplar

The SCL9 subfamily is a key member of the GRAS family that regulates plant development and stress responses. Nevertheless, the functional role of these genes in the growth and development of poplar still unclear. Here, we reported the six SCL9 genes, which were found to be differentially expressed during poplar adventitious root formation. The full-length sequences of PeSCL9 genes of ‘Nanlin895’ poplar (Populus deltoids × Populus euramericana) were cloned by the RACE technique All PeSCL9 genes lacked introns. RT-qPCR revealed that PeSCL9 genes displayed a dynamic expression pattern in the adventitious root of poplar, according to RT-qPCR data. A series of comprehensive genes characteristics analysis were carried out for six genes by bioinformation. Meanwhile, transient expression analysis of the Populus protoplasts showed that all the PeSCL9 proteins were localized in the nucleus. In addition, the degradome and sRNA of ‘Nanlin895’ poplar in combination were used to predict miRNAs that regulate PeSCL9. It was found that miR396a and miR396c may affect PeSCL9 expression via cleavage, which was further verified by a transient expression experiment in Populus protoplasts. Overall, the development of poplar adventitious root and other tissues was closely related to these six SCL9 genes, and they serve as a starting point for further research into the mechanisms regulating poplar growth and development.

The microRNA476a-RFL module regulates adventitious root formation through a mitochondria-dependent pathway in Populus

For woody plants, clonal propagation efficiency is largely determined by adventitious root (AR) formation at the bases of stem cuttings. However, our understanding of the molecular mechanisms contributing to AR morphogenesis in trees remains limited, despite the importance of vegetative propagation, currently the most common practice for tree breeding and commercialization. Here, we identified Populus-specific miR476a as a regulator of wound-induced adventitious rooting that acts by orchestrating mitochondrial homeostasis. MiR476a exhibited inducible expression during AR formation and directly targeted several Restorer of Fertility like (RFL) genes encoding mitochondrion-localized pentatricopeptide repeat proteins. Genetic modification of miR476a-RFL expression revealed that miR476a/RFL-mediated dynamic regulation of mitochondrial homeostasis influences AR formation in poplar. Mitochondrial perturbation via exogenous application of a chemical inhibitor indicated that miR476a/RFL-directed AR formation depends on mitochondrial regulation that acts via auxin signaling. Our results thus establish a microRNA-directed mitochondrion-auxin signaling cascade required for AR development, providing insights into the role of mitochondrial regulation in the developmental plasticity of plants.

Increased hippocampal excitability in miR-324-null mice

MicroRNAs are non-coding RNAs that act to downregulate the expression of target genes by translational repression and degradation of messenger RNA molecules. Individual microRNAs have the ability to specifically target a wide array of gene transcripts, therefore allowing each microRNA to play key roles in multiple biological pathways. miR-324 is a microRNA predicted to target thousands of RNA transcripts and is expressed far more highly in the brain than in any other tissue, suggesting that it may play a role in one or multiple neurological pathways. Here we present data from the first global miR-324-null mice, in which increased excitability and interictal discharges were identified in vitro in the hippocampus. RNA sequencing was used to identify differentially expressed genes in miR-324-null mice which may contribute to this increased hippocampal excitability, and 3'UTR luciferase assays and western blotting revealed that two of these, Suox and Cd300lf, are novel direct targets of miR-324. Characterisation of microRNAs that produce an effect on neurological activity, such as miR-324, and identification of the pathways they regulate will allow a better understanding of the processes involved in normal neurological function and in turn may present novel pharmaceutical targets in treating neurological disease.

Pan-Genome miRNomics in Brachypodium

Pan-genomes are efficient tools for the identification of conserved and varying genomic sequences within lineages of a species. Investigating genetic variations might lead to the discovery of genes present in a subset of lineages, which might contribute into beneficial agronomic traits such as stress resistance or yield. The content of varying genomic regions in the pan-genome could include protein-coding genes as well as microRNA(miRNAs), small non-coding RNAs playing key roles in the regulation of gene expression. In this study, we performed in silico miRNA identification from the genomic sequences of 54 lineages of Brachypodium distachyon, aiming to explore varying miRNA contents and their functional interactions. A total of 115 miRNA families were identified in 54 lineages, 56 of which were found to be present in all lineages. The miRNA families were classified based on their conservation among lineages and potential mRNA targets were identified. Obtaining information about regulatory mechanisms stemming from these miRNAs offers strong potential to provide a better insight into the complex traits that were potentially present in some lineages. Future work could lead us to introduce these traits to different lineages or other economically important plant species in order to promote their survival in different environmental conditions.

Next-Generation Sequencing Identification and Characterization of MicroRNAs in Dwarfed Citrus Trees Infected With Citrus Dwarfing Viroid in High-Density Plantings

Citrus dwarfing viroid (CDVd) induces stunting on sweet orange trees [Citrus sinensis (L.) Osbeck], propagated on trifoliate orange rootstock [Citrus trifoliata (L.), syn. Poncirus trifoliata (L.) Raf.]. MicroRNAs (miRNAs) are a class of non-coding small RNAs (sRNAs) that play important roles in the regulation of tree gene expression. To identify miRNAs in dwarfed citrus trees, grown in high-density plantings, and their response to CDVd infection, sRNA next-generation sequencing was performed on CDVd-infected and non-infected controls. A total of 1,290 and 628 miRNAs were identified in stem and root tissues, respectively, and among those, 60 were conserved in each of these two tissue types. Three conserved miRNAs (csi-miR479, csi-miR171b, and csi-miR156) were significantly downregulated (adjusted p-value < 0.05) in the stems of CDVd-infected trees compared to the non-infected controls. The three stem downregulated miRNAs are known to be involved in various physiological and developmental processes some of which may be related to the characteristic dwarfed phenotype displayed by CDVd-infected C. sinensis on C. trifoliata rootstock field trees. Only one miRNA (csi-miR535) was significantly downregulated in CDVd-infected roots and it was predicted to target genes controlling a wide range of cellular functions. Reverse transcription quantitative polymerase chain reaction analysis performed on selected miRNA targets validated the negative correlation between the expression levels of these targets and their corresponding miRNAs in CDVd-infected trees. Our results indicate that CDVd-responsive plant miRNAs play a role in regulating important citrus growth and developmental processes that may participate in the cellular changes leading to the observed citrus dwarf phenotype.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Small RNA and degradome deep sequencing reveals important roles of microRNAs in cotton (Gossypium hirsutum L.) response to root-knot nematode Meloidogyne incognita infection

Investigation of cotton response to nematode infection will allow us to better understand the cotton immune defense mechanism and design a better biotechnological approach for efficiently managing pest nematodes in cotton. In this study, we firstly treated cotton by root knot nematode (RKN, Meloidogyne incognita) infections, then we employed the high throughput deep sequencing technology to sequence and genome-widely identify all miRNAs in cotton; finally, we analyzed the functions of these miRNAs in cotton response to RKN infections. A total of 266 miRNAs, including 193 known and 73 novel miRNAs, were identified by deep sequencing technology, which belong to 67 conserved and 66 novel miRNA families, respectively. A majority of identified miRNA families only contain one miRNA; however, miR482 family contains 14 members and some others contain 2-13 members. Certain miRNAs were specifically expressed in RKN-infected cotton roots and others were completely inhibited by RKN infection. A total of 50 miRNAs were differentially expressed after RKN infection, in which 28 miRNAs were up-regulated and 22 were inhibited by RKN treatment. Based on degradome sequencing, 87 gene targets were identified to be targeted by 57 miRNAs. These miRNA-targeted genes are involved in the interaction of cotton plants and nematode infection. Based on GO (gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, 466 genes from all 636 miRNA targets were mapped to 6340 GO terms, 181 genes from 228 targets of differentially expressed miRNAs were mapped to 1588 GO terms. The GO terms were then categorized into the three main GO classes: biological processes, cellular components, and molecular functions. The targets of differentially expressed miRNAs were enriched in 43 GO terms, including 22 biological processes, 10 cellular components, and 11 molecular functions (p < 0.05). Many identified processes were associated with organism responses to the environmental stresses, including regulation of nematode larval development, response to nematode, and response to flooding. Our results will enhance the study and application of developing new cotton cultivars for nematode resistance.

Conservation and turnover of miRNAs and their highly complementary targets in early branching animals

MicroRNAs (miRNAs) are crucial post-transcriptional regulators that have been extensively studied in Bilateria, a group comprising the majority of extant animals, where more than 30 conserved miRNA families have been identified. By contrast, bilaterian miRNA targets are largely not conserved. Cnidaria is the sister group to Bilateria and thus provides a unique opportunity for comparative studies. Strikingly, like their plant counterparts, cnidarian miRNAs have been shown to predominantly have highly complementary targets leading to transcript cleavage by Argonaute proteins. Here, we assess the conservation of miRNAs and their targets by small RNA sequencing followed by miRNA target prediction in eight species of Anthozoa (sea anemones and corals), the earliest-branching cnidarian class. We uncover dozens of novel miRNAs but only a few conserved ones. Further, given their high complementarity, we were able to computationally identify miRNA targets in each species. Besides evidence for conservation of specific miRNA target sites, which are maintained between sea anemones and stony corals across 500 Myr of evolution, we also find indications for convergent evolution of target regulation by different miRNAs. Our data indicate that cnidarians have only few conserved miRNAs and corresponding targets, despite their high complementarity, suggesting a high evolutionary turnover.

Structural and Functional Characteristics of miRNAs in Five Strategic Millet Species and Their Utility in Drought Tolerance

Millets are the strategic food crops in arid and drought-prone ecologies. Millets, by virtue of nature, are very well-adapted to drought conditions and able to produce sustainable yield. Millets have important nutrients that can help prevent micro-nutrient malnutrition. As a result of the adverse effect of climate change and widespread malnutrition, millets have attained a strategic position to sustain food and nutritional security. Although millets can adapt well to the drought ecologies where other cereals fail completely, the yield level is very low under stress. There is a tremendous opportunity to increase the genetic potential of millet crops in dry lands when the genetics of the drought-tolerance mechanism is fully explained. MicroRNAs (miRNAs) are the class of small RNAs that control trait expression. They are part of the gene regulation but little studied in millets. In the present study, novel miRNAs and gene targets were identified from the genomic resources of pearl millet, sorghum, foxtail millet, finger millet, and proso millet through in silico approaches. A total of 1,002 miRNAs from 280 families regulating 23,158 targets were identified using different filtration criteria in five millet species. The unique as well as conserved structural features and functional characteristics of miRNA across millets were explained. About 84 miRNAs were conserved across millets in different species combinations, which explained the evolutionary relationship of the millets. Further, 215 miRNAs controlling 155 unique major drought-responsive genes, transcription factors, and protein families revealed the genetics of drought tolerance that are accumulated in the millet genomes. The miRNAs regulating the drought stress through specific targets or multiple targets showed through a network analysis. The identified genes regulated by miRNA genes could be useful in developing functional markers and used for yield improvement under drought in millets as well as in other crops.

miR-339-3p regulated acute pancreatitis induced by caerulein through targeting TNF receptor-associated factor 3 in AR42J cells

Acute pancreatitis (AP) is an inflammatory disease with high morbidity and mortality. The regulation mechanism of miRNA is involved in the production and development of various diseases, but the regulation mechanism of miRNA in AP is still not fully elucidated. The expression of miR-339-3p was detected using quantitative real-time PCR. The levels of TNF-α, IL-1β, and IL-6 were detected using enzyme-linked immunosorbent assay. Cell apoptosis was measured using flow cytometry. The protein expressions of TNF receptor-associated factor 3 (TRAF3), Bcl-2, C-caspase 3, Bax, p-p38, and p38 were measured using western blot. Luciferase reporter assay and RNA immunoprecipitation assay were applied to ensure that miR-399-3p targeted TRAF3. Caerulein promoted the expression of TNF-α, IL-1β, and IL-6, enhanced the expression of C-caspase 3 and Bax while inhibited Bcl-2 protein expression. Meanwhile, caerulein also reduced the expression of miR-339-3p and induced the expression of TRAF3 in rat pancreatic acinar cells. miR-399-3p transfection inhibited the levels of TNF-α, IL-1β, and IL-6 and C-caspase 3 and Bax protein expression as well as suppressed cell apoptosis, while increased Bcl-2 protein expression in caerulein-induced AP. TRAF3 has been verified as a target of miR-339-3p. Interestingly, the reduction of miR-399-3p inhibited the p38 pathway, which was impaired by the upregulation of TRAF3. In addition, the suppression effects of miR-339-3p on cell inflammation and apoptosis in caerulein-induced AP were reversed by enhancing TRAF3 expression. In this study, in vitro model of AP was characterized by strong inflammation and cell apoptosis. We have first demonstrated the regulatory network of miR-339-3p and TRAF3. Overexpression of miR-339-3p inhibited cell inflammation and cell apoptosis in caerulein-induced AP through modulating TRAF3 expression via the p38 pathway, providing a new therapeutic target in the treatment of AP.

CRISPR-based assessment of genomic structure in the conserved SQUAMOSA promoter-binding-like gene clusters in rice

Although SQUAMOSA promoter-binding-like (SPL) transcription factors are important regulators of development in rice (Oryza sativa), prior assessments of the SPL family have been limited to single genes. A functional comparison across the full gene family in standardized genetic backgrounds has not been reported previously. Here, we demonstrate that the SPL gene family in rice is enriched due to the most recent whole genome duplication (WGD). Notably, 10 of 19 rice SPL genes (52%) cluster in four units that have persisted for at least 50 million years. We show that SPL gene grouping and retention following WGD is widespread in angiosperms, suggesting the conservatism and importance of this gene arrangement. We used Cas9 editing to generate transformation lines for all 19 SPL genes in a common set of backgrounds, and found that knockouts of 14 SPL genes exhibited defects in plant height, 10 exhibited defects in panicle size, and nine had altered grain lengths. We observed subfunctionalization of genes in the paleoduplicated pairs, but little evidence of neofunctionalization. Expression of OsSPL3 was negatively correlated with that of its closest neighbor in its synteny group, OsSPL4, and its sister paired gene, OsSPL12, in the opposing group. Nucleotide diversity was lower in eight of the nine singleton genes in domesticated rice, relative to wild rice, whereas the reverse was true for the paired genes. Together, these results provide functional information on eight previously unexamined OsSPL family members and suggest that paleoduplicate pair redundancy benefits plant survival and innovation.

Non-coding RNA polymerases that silence transposable elements and reprogram gene expression in plants

Multisubunit RNA polymerase (Pol) complexes are the core machinery for gene expression in eukaryotes. The enzymes Pol I, Pol II and Pol III transcribe distinct subsets of nuclear genes. This family of nuclear RNA polymerases expanded in terrestrial plants by the duplication of Pol II subunit genes. Two Pol II-related enzymes, Pol IV and Pol V, are highly specialized in the production of regulatory, non-coding RNAs. Pol IV and Pol V are the central players of RNA-directed DNA methylation (RdDM), an RNA interference pathway that represses transposable elements (TEs) and selected genes. Genetic and biochemical analyses of Pol IV/V subunits are now revealing how these enzymes evolved from ancestral Pol II to sustain non-coding RNA biogenesis in silent chromatin. Intriguingly, Pol IV-RdDM regulates genes that influence flowering time, reproductive development, stress responses and plant–pathogen interactions. Pol IV target genes vary among closely related taxa, indicating that these regulatory circuits are often species-specific. Data from crops like maize, rice, tomato and Brassicarapa suggest that dynamic repositioning of TEs, accompanied by Pol IV targeting to TE-proximal genes, leads to the reprogramming of plant gene expression over short evolutionary timescales.

Aux/IAA14 Regulates microRNA-Mediated Cold Stress Response in Arabidopsis Roots

The phytohormone auxin and microRNA-mediated regulation of gene expressions are key regulators of plant growth and development at both optimal and under low-temperature stress conditions. However, the mechanistic link between microRNA and auxin in regulating plant cold stress response remains elusive. To better understand the role of microRNA (miR) in the crosstalk between auxin and cold stress responses, we took advantage of the mutants of Arabidopsis thaliana with altered response to auxin transport and signal. Screening of the mutants for root growth recovery after cold stress at 4 °C revealed that the auxin signaling mutant, solitary root 1 (slr1; mutation in Aux/IAA14), shows a hypersensitive response to cold stress. Genome-wide expression analysis of miRs in the wild-type and slr1 mutant roots using next-generation sequencing revealed 180 known and 71 novel cold-responsive microRNAs. Cold stress also increased the abundance of 26-31 nt small RNA population in slr1 compared with wild type. Comparative analysis of microRNA expression shows significant differential expression of 13 known and 7 novel miRs in slr1 at 4 °C compared with wild type. Target gene expression analysis of the members from one potential candidate miR, miR169, revealed the possible involvement of miR169/NF-YA module in the Aux/IAA14-mediated cold stress response. Taken together, these results indicate that SLR/IAA14, a transcriptional repressor of auxin signaling, plays a crucial role in integrating miRs in auxin and cold responses.

The trehalose 6-phosphate pathway impacts vegetative phase change in Arabidopsis thaliana

The vegetative phase change marks the beginning of the adult phase in the life cycle of plants and is associated with a gradual decline in the microRNA miR156, in response to sucrose status. Trehalose 6-phosphate (T6P) is a sugar molecule with signaling function reporting the current sucrose state. To elucidate the role of T6P signaling in vegetative phase change, molecular, genetic, and metabolic analyses were performed using Arabidopsis thaliana loss-of-function lines in TREHALOSE PHOSPHATE SYNTHASE1 (TPS1), a gene coding for an enzyme that catalyzes the production of T6P. These lines show a significant delay in vegetative phase change, under both short and long day conditions. Induced expression of TPS1 complements this delay in the TPS1 knockout mutant (tps1-2 GVG::TPS1). Further analyses indicate that the T6P pathway promotes vegetative phase transition by suppressing miR156 expression and thereby modulating the levels of its target transcripts, the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes. TPS1 knockdown plants, with a delayed vegetative phase change phenotype, accumulate significantly more sucrose than wild-type plants as a result of a feedback mechanism. In summary, we conclude that the T6P pathway forms an integral part of an endogenous mechanism that influences phase transitions dependent on the metabolic state.

Stress induced dynamic adjustment of conserved miR164:NAC module

Aims including the rationale

Salinity and drought are the two major stresses limiting the productivity of economically important crops such as Glycine max (soybean). The incidence of these stresses during the pod development stages affects the quality and quantity of seeds, which compromise the yield of soybean. The miR164:NAC module has been shown to play a critical role in regulating the response to salt and drought stress in several plant species. However, biological role of miR164:NAC module in salt stress in soybean is not fully understood.

Methods

In this study, we identified 215 salt responsive miRNAs, using miScript miRNA array with a sensitive and a tolerant soybean genotype, William82 and INCASoy36, respectively. The targets of these salt regulated miRNAs were searched in the degradome datasets.

Key results

It was found that four salt stress deregulated miRNAs targeted the NAC transcription factor and among these miR164k and miR408d showed antagonistic expression in the two soybean genotypes. The expression of miR164k was higher in salt tolerant INCASoy36 as compared to salt sensitive William82, under unstressed conditions. However under salt stress, miR164k was downregulated in INCASoy36 (-2.65 fold), whereas it was upregulated in William82 (4.68 fold). A transient co-expression assay validated that gma-miR164k directs the cleavage of GmNAC1 transcript. Bioinformatics analysis revealed that the regulation of NAC transcription factor family by members of miR164 family is conserved across many species. The dynamic expression profiles of miR164 and NAC-TFs were captured in different tissues of rice, tobacco, and two soybean genotypes under drought and salt stress conditions.

Main conclusion

Collectively, our results suggest that genetically determined dynamic modulation of the conserved miR164:NAC-TF module may play an important role in determining the adaptive response of plants to stress.

A role for the Auxin Response Factors ARF6 and ARF8 homologs in petal spur elongation and nectary maturation in Aquilegia

The petal spur of the basal eudicot Aquilegia is a key innovation associated with the adaptive radiation of the genus. Previous studies have shown that diversification of Aquilegia spur length can be predominantly attributed to variation in cell elongation. However, the genetic pathways that control the development of petal spurs are still being investigated. Here, we focus on a pair of closely related homologs of the AUXIN RESPONSE FACTOR family, AqARF6 and AqARF8, to explore their roles in Aquileiga coerulea petal spur development. Expression analyses of the two genes show that they are broadly expressed in vegetative and floral organs, but have relatively higher expression in petal spurs, particularly at later stages. Knockdown of the two AqARF6 and AqARF8 transcripts using virus-induced gene silencing resulted in largely petal-specific defects, including a significant reduction in spur length due to a decrease in cell elongation. These spurs also exhibited an absence of nectar production, which was correlated with downregulation of STYLISH homologs that have previously been shown to control nectary development. This study provides the first evidence of ARF6/8 homolog-mediated petal development outside the core eudicots. The genes appear to be specifically required for cell elongation and nectary maturation in the Aquilegia petal spur.

Plant-pathogen interactions: MicroRNA-mediated trans-kingdom gene regulation in fungi and their host plants

MicroRNAs (miRNAs) have been prevalently studied in plants, animals, and viruses. However, recent studies show evidences of miRNA-like RNAs (milRNAs) in fungi as well. It is known that after successful infection, pathogens hijack the host machinery and use it for their own growth and multiplication. Alternatively, resistant plants can overcome the pathogen attack by a variety of mechanisms. Based on this prior knowledge, we computationally predicted milRNAs from 13 fungi, and identified their targets in transcriptomes of the respective fungi as well as their host plants. The expressions of the milRNAs and targets were confirmed using qRT-PCR. We found that plant miRNAs targeted fungal virulence genes, while fungal milRNAs targeted plant resistance genes; corroborating miRNA-mediated trans-kingdom gene regulation and the roles of miRNAs in plant-pathogen interactions. Transgenic plants with miRNAs targeting fungal virulence genes, or anti-sense of fungal milRNAs, would be expected to be highly resistant to the fungal pathogens.