A computational-based update on microRNAs and their targets in barley (Hordeum vulgare L.)

Integrated transcriptome and miRNA sequencing analyses reveal that hypoxia stress induces immune and metabolic disorders in gill of genetically improved farmed tilapia (GIFT, Oreochromis niloticus)

The survival and growth of fish are significantly impacted by a hypoxic environment (low dissolved oxygen). In this study, we compared tissue structure, physiological changes, and mRNA/miRNA transcriptome, in gills of genetically improved farmed tilapia (GIFT, Oreochromis niloticus) between the hypoxic group (DO: 0.55 mg/L, HG) and the control group (DO: 5 mg/L, CG). The results showed that the gill filaments in the hypoxic group showed curling, engorgement, and apoptotic cells increased, and that exposure for 96 h resulted in a reduction in the antioxidant capacity. We constructed and sequenced miRNA and mRNA libraries from gill tissues of GIFT at 96 h of hypoxia stress. Between the HG and CG, a total of 14 differentially expressed (DE) miRNAs and 1557 DE genes were obtained. GO and KEGG enrichment showed that DE genes were mainly enriched in immune and metabolic pathways such as natural killer cell mediated cytotoxicity, steroid biosynthesis, primary immunodeficiency, and synthesis and degradation of ketone bodies. Based on the results of mRNA sequencing and screening for miRNA-mRNA pairs, we selected and verified six DE miRNAs and their probable target genes. The sequencing results were consistent with the qRT-PCR validation results. The result showed that under hypoxia stress, the innate immune response was up-regulated, and the adaptive immune response was down-regulated in the gill of GIFT. The synthesis of cholesterol in gill cells is reduced, which is conducive to the absorption of solvent oxygen. These findings offer fresh information about the processes of fish adaptation to hypoxic stress.

Integrated miRNA, target mRNA, and metabolome profiling of Tinospora cordifolia with reference to berberine biosynthesis

MicroRNAs play a central role in gene regulation and emerge as novel targets for secondary metabolites improvement in plants. The crops thus can be improved through knowledge obtained by the study of miRNAs because of their conserved nature in gene regulation. The present study has been carried out on Tinospora cordifolia (T. cordifolia), because of its illimitable application for the treatment of various diseases. This plant has tremendous medicinal properties, yet unexplored at the molecular level, and has not received much recognition in the scientific field. Thus, here computational analysis was performed to identify T. cordifolia miRNAs using EST database. Using these miRNAs, we predicted their targets which were found to be associated with the regulation of diverse gene networks including 433 berberine biosynthesis genes in T. cordifolia. Further, selected miRNAs were validated and their expression was detected in different T. cordifolia tissues followed by expression analysis of their target mRNAs. These data were then compared with the metabolic profile of T. cordifolia with an emphasis on therapeutically important compound berberine. In this study, we did simultaneous miRNA/target gene expression and metabolome analysis which opens a new way for initiating new proposition and prioritization of miRNAs/genes/metabolites for targeted follow‑up metabolic engineering experimentations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03342-9.

Genome-wide identification and characterization of caffeoyl-coenzyme A O-methyltransferase genes related to the Fusarium head blight response in wheat

Background

Lignin is one of the main components of the cell wall and is directly associated with plant development and defence mechanisms in plants, especially in response to Fusarium graminearum (Fg) infection. Caffeoyl-coenzyme A O-methyltransferase (CCoAOMT) is the main regulator determining the efficiency of lignin synthesis and composition. Although it has been characterized in many plants, to date, the importance of the CCoAOMT family in wheat is not well understood.

Results

Here, a total of 21 wheat CCoAOMT genes (TaCCoAOMT) were identified through an in silico genome search method and they were classified into four groups based on phylogenetic analysis, with the members of the same group sharing similar gene structures and conserved motif compositions. Furthermore, the expression patterns and co-expression network in which TaCCoAOMT is involved were comprehensively investigated using 48 RNA-seq samples from Fg infected and mock samples of 4 wheat genotypes. Combined with qRT-PCR validation of 11 Fg-responsive TaCCoAOMT genes, potential candidates involved in the FHB response and their regulation modules were preliminarily suggested. Additionally, we investigated the genetic diversity and main haplotypes of these CCoAOMT genes in bread wheat and its relative populations based on resequencing data.

Conclusions

This study identified and characterized the CCoAOMT family in wheat, which not only provided potential targets for further functional analysis, but also contributed to uncovering the mechanism of lignin biosynthesis and its role in FHB tolerance in wheat and beyond.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07849-y.

Current status of microRNA-mediated regulation of drought stress responses in cereals

Drought is one of the most important abiotic stress factors impeding crop productivity. With the uncovering of their role as potential regulators of gene expression, microRNAs (miRNAs) have been recognized as new targets for developing stress resistance. MicroRNAs are small noncoding RNAs whose abundance is significantly altered under stress conditions. Interestingly, plant miRNAs predominantly targets transcription factors (TFs), and some of which are also the most critical drought-responsive genes that in turn could regulate the expression of numerous loci with drought-adaptive potential. The phytohormone ABA plays important roles in regulating stomatal conductance and in initiating an adaptive response to drought stress. miRNAs are implicated in regulating ABA-(abscisic acid) and non-ABA-mediated drought resistance pathways. For instance, miR159-MYB module and miR169-NFYA module participates in an ABA-dependent pathway, whereas several other ABA-independent miRNA-target modules (miR156-SPL; miR393-TIR1; miR160-ARF10, ARF16, ARF17; miR167-ARF6 and ARF8; miR390/TAS3siRNA-ARF2, ARF3, ARF4) collectively regulate drought responses in plants. Overall, miRNA-mediated drought response manifests diverse molecular, biochemical and physiological processes. Because of their immense role in controlling gene expression, miRNA manipulation has significant potential to augment plant tolerance to drought stress. This review compiles the current understanding of drought-responsive miRNAs in major cereals. Also, potential miRNA manipulation strategies currently in use along with the challenges and future perspectives are discussed.

Genetic regulation of the traits contributing to wheat nitrogen use efficiency

High nitrogen application aimed at increasing crop yield is offset by higher production costs and negative environmental consequences. For wheat, only one third of the applied nitrogen is utilized, which indicates there is scope for increasing Nitrogen Use Efficiency (NUE). However, achieving greater NUE is challenged by the complexity of the trait, which comprises processes associated with nitrogen uptake, transport, reduction, assimilation, translocation and remobilization. Thus, knowledge of the genetic regulation of these processes is critical in increasing NUE. Although primary nitrogen uptake and metabolism-related genes have been well studied, the relative influence of each towards NUE is not fully understood. Recent attention has focused on engineering transcription factors and identification of miRNAs acting on expression of specific genes related to NUE. Knowledge obtained from model species needs to be translated into wheat using recently-released whole genome sequences, and by exploring genetic variations of NUE-related traits in wild relatives and ancient germplasm. Recent findings indicate the genetic basis of NUE is complex. Pyramiding various genes will be the most effective approach to achieve a satisfactory level of NUE in the field.

Responses of functional miRNA-mRNA regulatory modules to a high-fat diet in the liver of hybrid yellow catfish (Pelteobagrus fulvidraco × P. vachelli)

Fatty liver disease is common in cultured yellow catfish as a result of high fat contents in feeds. However, little is known about the mechanism by which the excessive deposition of liver fat causes fatty liver disease. Hybrid yellow catfish (Pelteobagrus fulvidraco♀ × P. vachelli♂) were fed a high-fat diet (HFD) or a normal-fat diet (NFD) for 60 days. Compared with the NFD group, the HFD group showed lower growth performance, higher hepatosomatic and viscerosomatic indexes, increased hepatic triglyceride and cholesterol contents, and more and larger lipid droplets in liver tissue. Whole transcriptome mRNA libraries and microRNA libraries from fish in the NFD and HFD groups were constructed by high-throughput sequencing. Twelve miRNAs were differentially expressed (DE) between the HFD and NFD groups. Seven negatively correlated DE miRNA-DE mRNA pairs were selected, and the expression patterns of both were confirmed using qRT-PCR. Hybrid yellow catfish showed mediated oxidative degradation of liver glucose and fatty acid peroxidation, regulation of antioxidant enzyme activity, and various immune and inflammatory responses to fat deposition and stress. These findings have important biological significance for protecting the liver against stress, as well as economic significance for establishing healthy aquaculture conditions.

Identification of Known and Novel MicroRNAs in Raspberry Organs Through High-Throughput Sequencing

MicroRNAs (miRNAs) are a class of small endogenous RNAs that play important regulatory roles in plants by negatively affecting gene expression. Studies on the identification of miRNAs and their functions in various plant species and organs have significantly contributed to plant development research. In the current study, we utilized high-throughput sequencing to detect the miRNAs in the root, stem, and leaf tissues of raspberry (Rubus idaeus). A total of more than 35 million small RNA reads ranging in size from 18 to 35 nucleotides were obtained, with 147 known miRNAs and 542 novel miRNAs identified among the three organs. Sequence verification and the relative expression profiles of the six known miRNAs were investigated by stem-loop quantitative real-time PCR. Furthermore, the potential target genes of the known and novel miRNAs were predicted and subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway annotation. Enrichment analysis of the GO-associated biological processes and molecular functions revealed that these target genes were potentially involved in a wide range of metabolic pathways and developmental processes. Moreover, the miRNA target prediction revealed that most of the targets predicted as transcription factor-coding genes are involved in cellular and metabolic processes. This report is the first to identify miRNAs in raspberry. The detected miRNAs were analyzed by cluster analysis according to their expression, which revealed that these conservative miRNAs are necessary for plant functioning. The results add novel miRNAs to the raspberry transcriptome, providing a useful resource for the further elucidation of the functional roles of miRNAs in raspberry growth and development.

Identification, characterization and expression analysis of passion fruit (Passiflora edulis) microRNAs

microRNAs (miRNAs) are highly conserved, short (~ 21-nucleotide), endogenous, non-coding RNA molecules that play major roles in post-transcriptional silencing by guiding target mRNA cleavage or translational inhibition. In this study, applying high-stringent genome-wide computational-based approaches, a total of 28 putative miRNAs belonging to 17 miRNA families were identified from an antioxidant-rich medicinal plant passion fruit (Passiflora edulis). Inter-tissue (leaves and fruits) and inter-varietal (yellow and purple fruit varieties) quantitative study of six putative passion fruit miRNAs (ped-miR160, ped-miR164, ped-miR166, ped-miR393, ped-miR394, and ped-miR398) showed differential expression. Using psRNATarget tool, a total of 25 potential target proteins of the characterized passion fruit miRNAs were also identified. Most of the target proteins identified in this study, including SQUAMOSA promoter binding, Class III HD-Zip, NAC, Scarecrow, APETALA2, Auxin response factor, MYB, and superoxide dismutase, were found to be involved in development, metabolism, and defense/stress response signaling.

The Use of Nitrogen and Its Regulation in Cereals: Structural Genes, Transcription Factors, and the Role of miRNAs

Cereals and, especially, rice, maize, and wheat, are essential commodities, on which human nutrition is based. Expanding population and food demand have required higher production which has been achieved by increasing fertilization, and especially nitrogen supply to cereal crops. In fact, nitrogen is a crucial nutrient for the plant, but excessive use poses serious environmental and health issues. Therefore, increasing nitrogen use efficiency in cereals is of pivotal importance for sustainable agriculture. The main steps in the use of nitrogen are uptake and transport, reduction and assimilation, and translocation and remobilization. Many studies have been carried out on the genes involved in these phases, and on transcription factors regulating these genes. Lately, increasing attention has been paid to miRNAs responding to abiotic stress, including nutrient deficiency. Many miRNAs have been found to regulate transcription factors acting on the expression of specific genes for nitrogen uptake or remobilization. Recent studies on gene regulatory networks have also demonstrated that miRNAs can interact with several nodes in the network, functioning as key regulators in nitrogen metabolism.

MicroRNA Sequencing Revealed Citrus Adaptation to Long-Term Boron Toxicity through Modulation of Root Development by miR319 and miR171

Boron (B) toxicity in Citrus is a common physiological disorder leading to reductions in both productivity and quality. Studies on how Citrus roots evade B toxicity may provide new insight into plant tolerance to B toxicity. Here, using Illumina sequencing, differentially expressed microRNAs (miRNAs) were identified in B toxicity-treated Citrus sinensis (tolerant) and C. grandis (intolerant) roots. The results showed that 37 miRNAs in C. grandis and 11 miRNAs in C. sinensis were differentially expressed when exposed to B toxicity. Among them, miR319, miR171, and miR396g-5p were confirmed via 5'-RACE and qRT-PCR to target a myeloblastosis (MYB) transcription factor gene, a SCARECROW-like protein gene, and a cation transporting ATPase gene, respectively. Maintenance of SCARECROW expression in B treated Citrus roots might fulfill stem cell maintenance, quiescent center, and endodermis specification, thus allowing regular root elongation under B-toxic stress. Down-regulation of MYB due to up-regulation of miR319 in B toxicity-treated C. grandis roots might decrease the number of root tips, thereby dramatically changing root system architecture. Our findings suggested that miR319 and miR171 play a pivotal role in Citrus adaptation to long-term B toxicity by targeting MYB and SCARECROW, respectively, both of which are responsible for root growth and development.

miR172 downregulates the translation of cleistogamy 1 in barley

BACKGROUND AND AIMS

Floret opening in barley is induced by the swelling of the lodicule, a trait under the control of the cleistogamy1 (cly1) gene. The product of cly1 is a member of the APETALA2 (AP2) transcription factor family, which inhibits lodicule development. A sequence polymorphism at the miR172 target site within cly1 has been associated with variation in lodicule development and hence with the cleistogamous phenotype. It was unclear whether miR172 actually functions in cly1 regulation and, if it does, which miR172 gene contributes to cleistogamy. It was also interesting to explore whether miR172-mediated cly1 regulation occurs at transcriptional level or at translational level.

METHODS

Deep sequencing of small RNA identified the miR172 sequences expressed in barley immature spikes. miR172 genes were confirmed by computational and expression analysis. miR172 and cly1 expression profiles were determined by in situ hybridization and quantitative expression analysis. Immunoblot analysis provided the CLY1 protein quantifications. Definitive evidence of the role of miR172 in cleistogamy was provided by a transposon Ds-induced mutant of Hv-miR172a.

KEY RESULTS

A small RNA analysis of the immature barley spike revealed three isomers, miR172a, b and c, of which miR172a was the most abundant. In situ hybridization analysis showed that miR172 and cly1 co-localize in the lodicule primordium, suggesting that these two molecules potentially interact with one another. Immunoblot analysis showed that the sequence polymorphism at the miR172 target site within cly1 reduced the abundance of the CLY1 protein, but not that of its transcript. In a Ds-induced mutant of Hv-miR172a, which generates no mature miR172a, the lodicules fail to grow, resulting in a very small lodicule.

CONCLUSIONS

Direct evidence is presented to show that miR172a acts to reduce the abundance of the CLY1 protein, which enables open flowering in barley.

The use of high-throughput small RNA sequencing reveals differentially expressed microRNAs in response to aster yellows phytoplasma-infection in Vitis vinifera cv. 'Chardonnay'

Phytoplasmas are cell wall-less plant pathogenic bacteria responsible for major crop losses throughout the world. In grapevine they cause grapevine yellows, a detrimental disease associated with a variety of symptoms. The high economic impact of this disease has sparked considerable interest among researchers to understand molecular mechanisms related to pathogenesis. Increasing evidence exist that a class of small non-coding endogenous RNAs, known as microRNAs (miRNAs), play an important role in post-transcriptional gene regulation during plant development and responses to biotic and abiotic stresses. Thus, we aimed to dissect complex high-throughput small RNA sequencing data for the genome-wide identification of known and novel differentially expressed miRNAs, using read libraries constructed from healthy and phytoplasma-infected Chardonnay leaf material. Furthermore, we utilised computational resources to predict putative miRNA targets to explore the involvement of possible pathogen response pathways. We identified multiple known miRNA sequence variants (isomiRs), likely generated through post-transcriptional modifications. Sequences of 13 known, canonical miRNAs were shown to be differentially expressed. A total of 175 novel miRNA precursor sequences, each derived from a unique genomic location, were predicted, of which 23 were differentially expressed. A homology search revealed that some of these novel miRNAs shared high sequence similarity with conserved miRNAs from other plant species, as well as known grapevine miRNAs. The relative expression of randomly selected known and novel miRNAs was determined with real-time RT-qPCR analysis, thereby validating the trend of expression seen in the normalised small RNA sequencing read count data. Among the putative miRNA targets, we identified genes involved in plant morphology, hormone signalling, nutrient homeostasis, as well as plant stress. Our results may assist in understanding the role that miRNA pathways play during plant pathogenesis, and may be crucial in understanding disease symptom development in aster yellows phytoplasma-infected grapevines.

Cross-Kingdom Regulation of Putative miRNAs Derived from Happy Tree in Cancer Pathway: A Systems Biology Approach

MicroRNAs (miRNAs) are well-known key regulators of gene expression primarily at the post-transcriptional level. Plant-derived miRNAs may pass through the gastrointestinal tract, entering into the body fluid and regulate the expression of endogenous mRNAs. Camptotheca acuminata, a highly important medicinal plant known for its anti-cancer potential was selected to investigate cross-kingdom regulatory mechanism and involvement of miRNAs derived from this plant in cancer-associated pathways through in silico systems biology approach. In this study, total 33 highly stable putative novel miRNAs were predicted from the publically available 53,294 ESTs of C. acuminata, out of which 14 miRNAs were found to be regulating 152 target genes in human. Functional enrichment, gene-disease associations and network analysis of these target genes were carried out and the results revealed their association with prominent types of cancers like breast cancer, leukemia and lung cancer. Pathways like focal adhesion, regulation of lipolysis in adipocytes and mTOR signaling pathways were found significantly associated with the target genes. The regulatory network analysis showed the association of some important hub proteins like GSK3B, NUMB, PEG3, ITGA2 and DLG2 with cancer-associated pathways. Based on the analysis results, it can be suggested that the ingestion of the C. acuminata miRNAs may have a functional impact on tumorigenesis in a cross-kingdom way and may affect the physiological condition at genetic level. Thus, the predicted miRNAs seem to hold potentially significant role in cancer pathway regulation and therefore, may be further validated using in vivo experiments for a better insight into their mechanism of epigenetic action of miRNA.

Dehydration-responsive miRNAs in foxtail millet: genome-wide identification, characterization and expression profiling

A set of novel and known dehydration-responsive miRNAs have been identified in foxtail millet. These findings provide new insights into understanding the functional role of miRNAs and their respective targets in regulating plant response to dehydration stress. MicroRNAs perform significant regulatory roles in growth, development and stress response of plants. Though the miRNA-mediated gene regulatory networks under dehydration stress remain largely unexplored in plant including foxtail millet (Setaria italica), which is a natural abiotic stress tolerant crop. To find out the dehydration-responsive miRNAs at the global level, four small RNA libraries were constructed from control and dehydration stress treated seedlings of two foxtail millet cultivars showing contrasting tolerance behavior towards dehydration stress. Using Illumina sequencing technology, 55 known and 136 novel miRNAs were identified, representing 22 and 48 miRNA families, respectively. Eighteen known and 33 novel miRNAs were differentially expressed during dehydration stress. After the stress treatment, 32 dehydration-responsive miRNAs were up-regulated in tolerant cultivar and 22 miRNAs were down-regulated in sensitive cultivar, suggesting that miRNA-mediated molecular regulation might play important roles in providing contrasting characteristics to these cultivars. Predicted targets of identified miRNAs were found to encode various transcription factors and functional enzymes, indicating their involvement in broad spectrum regulatory functions and biological processes. Further, differential expression patterns of seven known miRNAs were validated by northern blot and expression of ten novel dehydration-responsive miRNAs were confirmed by SL-qRT PCR. Differential expression behavior of five miRNA-target genes was verified under dehydration stress treatment and two of them also validated by RLM RACE. Overall, the present study highlights the importance of dehydration stress-associated post-transcriptional regulation governed by miRNAs and their targets in a naturally stress-tolerant model crop.

Transcriptome-wide identification of Rauvolfia serpentina microRNAs and prediction of their potential targets

MicroRNAs (miRNAs) are small non-coding RNAs of ∼ 19-24 nucleotides (nt) in length and considered as potent regulators of gene expression at transcriptional and post-transcriptional levels. Here we report the identification and characterization of 15 conserved miRNAs belonging to 13 families from Rauvolfia serpentina through in silico analysis of available nucleotide dataset. The identified mature R. serpentina miRNAs (rse-miRNAs) ranged between 20 and 22nt in length, and the average minimal folding free energy index (MFEI) value of rse-miRNA precursor sequences was found to be -0.815 kcal/mol. Using the identified rse-miRNAs as query, their potential targets were predicted in R. serpentina and other plant species. Gene Ontology (GO) annotation showed that predicted targets of rse-miRNAs include transcription factors as well as genes involved in diverse biological processes such as primary and secondary metabolism, stress response, disease resistance, growth, and development. Few rse-miRNAs were predicted to target genes of pharmaceutically important secondary metabolic pathways such as alkaloids and anthocyanin biosynthesis. Phylogenetic analysis showed the evolutionary relationship of rse-miRNAs and their precursor sequences to homologous pre-miRNA sequences from other plant species. The findings under present study besides giving first hand information about R. serpentina miRNAs and their targets, also contributes towards the better understanding of miRNA-mediated gene regulatory processes in plants.

Bioinformatic identification and expression analysis of banana microRNAs and their targets

MicroRNAs (miRNAs) represent a class of endogenous non-coding small RNAs that play important roles in multiple biological processes by degrading targeted mRNAs or repressing mRNA translation. Thousands of miRNAs have been identified in many plant species, whereas only a limited number of miRNAs have been predicted in M. acuminata (A genome) and M. balbisiana (B genome). Here, previously known plant miRNAs were BLASTed against the Expressed Sequence Tag (EST) and Genomic Survey Sequence (GSS), a database of banana genes. A total of 32 potential miRNAs belonging to 13 miRNAs families were detected using a range of filtering criteria. 244 miRNA:target pairs were subsequently predicted, most of which encode transcription factors or enzymes that participate in the regulation of development, growth, metabolism, and other physiological processes. In order to validate the predicted miRNAs and the mutual relationship between miRNAs and their target genes, qRT-PCR was applied to detect the tissue-specific expression levels of 12 putative miRNAs and 6 target genes in roots, leaves, flowers, and fruits. This study provides some important information about banana pre-miRNAs, mature miRNAs, and miRNA target genes and these findings can be applied to future research of miRNA functions.

Role of microRNAs in plant drought tolerance

Drought is a normal and recurring climate feature in most parts of the world and plays a major role in limiting crop productivity. However, plants have their own defence systems to cope with adverse climatic conditions. One of these defence mechanisms is the reprogramming of gene expression by microRNAs (miRNAs). miRNAs are small noncoding RNAs of approximately 22 nucleotides length, which have emerged as important regulators of genes at post-transcriptional levels in a range of organisms. Some miRNAs are functionally conserved across plant species and are regulated by drought stress. These properties suggest that miRNA-based genetic modifications have the potential to enhance drought tolerance in cereal crops. This review summarizes the current understanding of the regulatory mechanisms of plant miRNAs, involvement of plant miRNAs in drought stress responses in barley (Hordeum vulgare L.), wheat (Triticum spp.) and other plant species, and the involvement of miRNAs in plant-adaptive mechanisms under drought stress. Potential strategies and directions for future miRNA research and the utilization of miRNAs in the improvement of cereal crops for drought tolerance are also discussed.

Post-transcriptional and post-translational regulations of drought and heat response in plants: a spider's web of mechanisms

Drought and heat tolerance are complex quantitative traits. Moreover, the adaptive significance of some stress-related traits is more related to plant survival than to agronomic performance. A web of regulatory mechanisms fine-tunes the expression of stress-related traits and integrates both environmental and developmental signals. Both post-transcriptional and post-translational modifications contribute substantially to this network with a pivotal regulatory function of the transcriptional changes related to cellular and plant stress response. Alternative splicing and RNA-mediated silencing control the amount of specific transcripts, while ubiquitin and SUMO modify activity, sub-cellular localization and half-life of proteins. Interactions across these modification mechanisms ensure temporally and spatially appropriate patterns of downstream-gene expression. For key molecular components of these regulatory mechanisms, natural genetic diversity exists among genotypes with different behavior in terms of stress tolerance, with effects upon the expression of adaptive morphological and/or physiological target traits.

Identification of conserved and novel microRNAs in Catharanthus roseus by deep sequencing and computational prediction of their potential targets

MicroRNAs are small endogenous non-coding RNAs of ~19-24 nucleotides and perform regulatory roles in many plant processes. To identify miRNAs involved in regulatory networks controlling diverse biological processes including secondary metabolism in Catharanthus roseus, an important medicinal plant, we employed deep sequencing of small RNA from leaf tissue. A total of 88 potential miRNAs comprising of 81 conserved miRNAs belonging to 35 families and seven novel miRNAs were identified. Precursors for 16 conserved and seven novel cro-miRNAs were identified, and their stem-loop hairpin structures were predicted. Selected cro-miRNAs were analyzed by stem-loop qRT-PCR and differential expression patterns were observed in different vegetative tissues of C. roseus. Targets were predicted for conserved and novel cro-miRNAs, which were found to be involved in diverse biological role(s) including secondary metabolism. Our study enriches available resources and information regarding miRNAs and their potential targets for better understanding of miRNA-mediated gene regulation in plants.

MicroRNAs responding to southern rice black-streaked dwarf virus infection and their target genes associated with symptom development in rice

Southern rice black-streaked dwarf virus (SRBSDV) is a recently emerged rice virus that has spread across Asia. This devastating virus causes rice plants to produce a variety of symptoms during different growth stages. MicroRNAs (miRNAs) comprise a large group of 21-24-nt RNA molecules that are important regulators of plant development processes and stress responses. In this study, we used microarray profiling to investigate rice miRNAs responding to SRBSDV infection at 3, 9, 15, and 20 days post-inoculation (dpi). Expression levels of 56 miRNAs were altered in SRBSDV-infected rice plants, with these changes classified into eight different regulation patterns according to their temporal expression dynamics. Fourteen miRNAs belonging to six families (miR164, R396, R530, R1846, R1858, and R2097) were significantly regulated at 20 dpi. We used RT-qPCR to search for expression level correlations between members of these families and their putative targets at 3, 9, and 15 dpi. Some members of the miR164, R396, R530, and R1846 families were found to be positively or negatively correlated with their respective targets during 3-15 days after SRBSDV infection, whereas in more cases the rice miRNAs were not in correlation with their targets along the post-inoculation period, suggesting that some additional factors may be involved in rice miRNA-target interactions. The reported functions of rice genes targeted by the miR164, R396, R530, R1846, and R1858 families indicated that these genes are associated with symptom development. These results provide insights into miRNA-mediated SRBSDV-rice interactions.

Comparison of the transcriptomes of ginger (Zingiber officinale Rosc.) and mango ginger (Curcuma amada Roxb.) in response to the bacterial wilt infection

Bacterial wilt in ginger (Zingiber officinale Rosc.) caused by Ralstonia solanacearum is one of the most important production constraints in tropical, sub-tropical and warm temperature regions of the world. Lack of resistant genotype adds constraints to the crop management. However, mango ginger (Curcuma amada Roxb.), which is resistant to R. solanacearum, is a potential donor, if the exact mechanism of resistance is understood. To identify genes involved in resistance to R. solanacearum, we have sequenced the transcriptome from wilt-sensitive ginger and wilt-resistant mango ginger using Illumina sequencing technology. A total of 26387032 and 22268804 paired-end reads were obtained after quality filtering for C. amada and Z. officinale, respectively. A total of 36359 and 32312 assembled transcript sequences were obtained from both the species. The functions of the unigenes cover a diverse set of molecular functions and biological processes, among which we identified a large number of genes associated with resistance to stresses and response to biotic stimuli. Large scale expression profiling showed that many of the disease resistance related genes were expressed more in C. amada. Comparative analysis also identified genes belonging to different pathways of plant defense against biotic stresses that are differentially expressed in either ginger or mango ginger. The identification of many defense related genes differentially expressed provides many insights to the resistance mechanism to R. solanacearum and for studying potential pathways involved in responses to pathogen. Also, several candidate genes that may underline the difference in resistance to R. solanacearum between ginger and mango ginger were identified. Finally, we have developed a web resource, ginger transcriptome database, which provides public access to the data. Our study is among the first to demonstrate the use of Illumina short read sequencing for de novo transcriptome assembly and comparison in non-model species of Zingiberaceae.

Identification of boron-deficiency-responsive microRNAs in Citrus sinensis roots by Illumina sequencing

BACKGROUND

Boron (B)-deficiency is a widespread problem in many crops, including Citrus. MicroRNAs (miRNAs) play important roles in nutrient deficiencies. However, little is known on B-deficiency-responsive miRNAs in plants. In this study, we first identified miRNAs and their expression pattern in B-deficient Citrus sinensis roots by Illumina sequencing in order to identify miRNAs that might be involved in the tolerance of plants to B-deficiency.

RESULTS

We isolated 52 (40 known and 12 novel) up-regulated and 82 (72 known and 10 novel) down-regulated miRNAs from B-deficient roots, demonstrating remarkable metabolic flexibility of roots, which might contribute to the tolerance of plants to B-deficiency. A model for the possible roles of miRNAs in the tolerance of roots to B-deficiency was proposed. miRNAs might regulate the adaptations of roots to B-deficiency through following several aspects: (a) inactivating reactive oxygen species (ROS) signaling and scavenging through up-regulating miR474 and down-regulating miR782 and miR843; (b) increasing lateral root number by lowering miR5023 expression and maintaining a certain phenotype favorable for B-deficiency-tolerance by increasing miR394 expression; (c) enhancing cell transport by decreasing the transcripts of miR830, miR5266 and miR3465; (d) improving osmoprotection (miR474) and regulating other metabolic reactions (miR5023 and miR821). Other miRNAs such as miR472 and miR2118 in roots increased in response to B-deficiency, thus decreasing the expression of their target genes, which are involved in disease resistance, and hence, the disease resistance of roots.

CONCLUSIONS

Our work demonstrates the possible roles of miRNAs and related mechanisms in the response of plant roots to B-deficiency.

Transcriptome analysis of the Portunus trituberculatus: de novo assembly, growth-related gene identification and marker discovery

Background

The swimming crab, Portunus trituberculatus, is an important farmed species in China, has been attracting extensive studies, which require more and more genome background knowledge. To date, the sequencing of its whole genome is unavailable and transcriptomic information is also scarce for this species. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset for major tissues of Portunus trituberculatus by the Illumina paired-end sequencing technology.

Results

Total RNA was isolated from eyestalk, gill, heart, hepatopancreas and muscle. Equal quantities of RNA from each tissue were pooled to construct a cDNA library. Using the Illumina paired-end sequencing technology, we generated a total of 120,137 transcripts with an average length of 1037 bp. Further assembly analysis showed that all contigs contributed to 87,100 unigenes, of these, 16,029 unigenes (18.40% of the total) can be matched in the GenBank non-redundant database. Potential genes and their functions were predicted by GO, KEGG pathway mapping and COG analysis. Based on our sequence analysis and published literature, many putative genes with fundamental roles in growth and muscle development, including actin, myosin, tropomyosin, troponin and other potentially important candidate genes were identified for the first time in this specie. Furthermore, 22,673 SSRs and 66,191 high-confidence SNPs were identified in this EST dataset.

Conclusion

The transcriptome provides an invaluable new data for a functional genomics resource and future biological research in Portunus trituberculatus. The data will also instruct future functional studies to manipulate or select for genes influencing growth that should find practical applications in aquaculture breeding programs. The molecular markers identified in this study will provide a material basis for future genetic linkage and quantitative trait loci analyses, and will be essential for accelerating aquaculture breeding programs with this species.

Identification and characterization of the microRNA transcriptome of a moth orchid Phalaenopsis aphrodite

Orchids display unique phenotypes, functional characteristics and ecological adaptations that are not found in model plants. In this study, we aimed to characterize the microRNA (miRNA) transcriptome and identify species- and tissue-specific miRNAs in Phalaenopsis aphrodite. After data filtering and cleanup, a total of 59,387,374 reads, representing 1,649,996 unique reads, were obtained from four P. aphrodite small RNA libraries. A systematic bioinformatics analysis pipeline was developed that can be used for miRNA and precursor mining, and target gene prediction in non-model plants. A total of 3,251 unique reads for 181 known plant miRNAs (belonging to 88 miRNA families), 23 new miRNAs and 91 precursors were identified. All the miRNA star sequences (miRNA*), the complementary strands of miRNA that from miRNA/miRNA* duplexes, of the predicted new miRNAs were detected in our small RNA libraries, providing additional evidence for their existence as new miRNAs in P. aphrodite. Furthermore, 240 potential miRNA-targets that appear to be involved in many different biological activities and molecular functions, especially transcription factors, were identified, suggesting that miRNAs can impact multiple processes in P. aphrodite. We also verified the cleavage sites for six targets using RNA ligase-mediated rapid amplification of 5' ends assay. The results provide valuable information about the composition, expression and function of miRNA in P. aphrodite, and will aid functional genomics studies of orchids.

Identification of potential miRNAs and their targets in Vriesea carinata (Poales, Bromeliaceae)

The miRNAs play important roles in regulation of gene expression at the post-transcriptional level. A small RNA and RNA-seq of libraries were constructed to identify miRNAs in Vriesea carinata, a native bromeliad species from Brazilian Atlantic Rainforest. Illumina technology was used to perform high throughput sequencing and data was analyzed using bioinformatics tools. We obtained 2,191,509 mature miRNAs sequences representing 54 conserved families in plant species. Further analysis allowed the prediction of secondary structures for 19 conserved and 16 novel miRNAs. Potential targets were predicted from pre-miRNAs by sequence homology and validated using RTqPCR approach. This study provides the first identification of miRNAs and their potential targets of a bromeliad species.

Evaluation of three RT-qPCR-based miRNA detection methods using seven rice miRNAs

Three frequently-used reverse transcription-quantitative polymerase chain reaction (RT-qPCR)-based miRNA detection methods, stem-loop RT-qPCR, poly(A)-tailing RT-qPCR, and miQPCR, were evaluated using seven selected rice miRNAs. The results revealed that miRNA abundance and sequence characteristics can affect capability of detection. The stem-loop amplification technique detected highly and moderately abundant miRNAs. The poly(A)-tailing method detected both highly abundant and sparsely present miRNAs, but failed to detect miRNAs with a hairpin structure. Only a few miRNAs were detectable by the miQPCR method. We suggest that a combination of methods should be used for reliable quantitative investigation of miRNAs.

Transcriptome analysis reveals ginsenosides biosynthetic genes, microRNAs and simple sequence repeats in Panax ginseng C. A. Meyer

BACKGROUND

Panax ginseng C. A. Meyer is one of the most widely used medicinal plants. Complete genome information for this species remains unavailable due to its large genome size. At present, analysis of expressed sequence tags is still the most powerful tool for large-scale gene discovery. The global expressed sequence tags from P. ginseng tissues, especially those isolated from stems, leaves and flowers, are still limited, hindering in-depth study of P. ginseng.

RESULTS

Two 454 pyrosequencing runs generated a total of 2,423,076 reads from P. ginseng roots, stems, leaves and flowers. The high-quality reads from each of the tissues were independently assembled into separate and shared contigs. In the separately assembled database, 45,849, 6,172, 4,041 and 3,273 unigenes were only found in the roots, stems, leaves and flowers database, respectively. In the jointly assembled database, 178,145 unigenes were observed, including 86,609 contigs and 91,536 singletons. Among the 178,145 unigenes, 105,522 were identified for the first time, of which 65.6% were identified in the stem, leaf or flower cDNA libraries of P. ginseng. After annotation, we discovered 223 unigenes involved in ginsenoside backbone biosynthesis. Additionally, a total of 326 potential cytochrome P450 and 129 potential UDP-glycosyltransferase sequences were predicted based on the annotation results, some of which may encode enzymes responsible for ginsenoside backbone modification. A BLAST search of the obtained high-quality reads identified 14 potential microRNAs in P. ginseng, which were estimated to target 100 protein-coding genes, including transcription factors, transporters and DNA binding proteins, among others. In addition, a total of 13,044 simple sequence repeats were identified from the 178,145 unigenes.

CONCLUSIONS

This study provides global expressed sequence tags for P. ginseng, which will contribute significantly to further genome-wide research and analyses in this species. The novel unigenes identified here enlarge the available P. ginseng gene pool and will facilitate gene discovery. In addition, the identification of microRNAs and the prediction of targets from this study will provide information on gene transcriptional regulation in P. ginseng. Finally, the analysis of simple sequence repeats will provide genetic makers for molecular breeding and genetic applications in this species.

Boron stress responsive microRNAs and their targets in barley

Boron stress is an environmental factor affecting plant development and production. Recently, microRNAs (miRNAs) have been found to be involved in several plant processes such as growth regulation and stress responses. In this study, miRNAs associated with boron stress were identified and characterized in barley. miRNA profiles were also comparatively analyzed between root and leave samples. A total of 31 known and 3 new miRNAs were identified in barley; 25 of them were found to respond to boron treatment. Several miRNAs were expressed in a tissue specific manner; for example, miR156d, miR171a, miR397, and miR444a were only detected in leaves. Additionally, a total of 934 barley transcripts were found to be specifically targeted and degraded by miRNAs. In silico analysis of miRNA target genes demonstrated that many miRNA targets are conserved transcription factors such as Squamosa promoter-binding protein, Auxin response factor (ARF), and the MYB transcription factor family. A majority of these targets were responsible for plant growth and response to environmental changes. We also propose that some of the miRNAs in barley such as miRNA408 might play critical roles against boron exposure. In conclusion, barley may use several pathways and cellular processes targeted by miRNAs to cope with boron stress.

Identification of microRNAs from Eugenia uniflora by high-throughput sequencing and bioinformatics analysis

Background

microRNAs or miRNAs are small non-coding regulatory RNAs that play important functions in the regulation of gene expression at the post-transcriptional level by targeting mRNAs for degradation or inhibiting protein translation. Eugenia uniflora is a plant native to tropical America with pharmacological and ecological importance, and there have been no previous studies concerning its gene expression and regulation. To date, no miRNAs have been reported in Myrtaceae species.

Results

Small RNA and RNA-seq libraries were constructed to identify miRNAs and pre-miRNAs in Eugenia uniflora. Solexa technology was used to perform high throughput sequencing of the library, and the data obtained were analyzed using bioinformatics tools. From 14,489,131 small RNA clean reads, we obtained 1,852,722 mature miRNA sequences representing 45 conserved families that have been identified in other plant species. Further analysis using contigs assembled from RNA-seq allowed the prediction of secondary structures of 25 known and 17 novel pre-miRNAs. The expression of twenty-seven identified miRNAs was also validated using RT-PCR assays. Potential targets were predicted for the most abundant mature miRNAs in the identified pre-miRNAs based on sequence homology.

Conclusions

This study is the first large scale identification of miRNAs and their potential targets from a species of the Myrtaceae family without genomic sequence resources. Our study provides more information about the evolutionary conservation of the regulatory network of miRNAs in plants and highlights species-specific miRNAs.

Identification of common carp (Cyprinus carpio) microRNAs and microRNA-related SNPs

BACKGROUND

MicroRNAs (miRNAs) exist pervasively across viruses, plants and animals and play important roles in the post-transcriptional regulation of genes. In the common carp, miRNA targets have not been investigated. In model species, single-nucleotide polymorphisms (SNPs) have been reported to impair or enhance miRNA regulation as well as to alter miRNA biogenesis. SNPs are often associated with diseases or traits. To date, no studies into the effects of SNPs on miRNA biogenesis and regulation in the common carp have been reported.

RESULTS

Using homology-based prediction combined with small RNA sequencing, we have identified 113 common carp mature miRNAs, including 92 conserved miRNAs and 21 common carp specific miRNAs. The conserved miRNAs had significantly higher expression levels than the specific miRNAs. The miRNAs were clustered into three phylogenetic groups. Totally 394 potential miRNA binding sites in 206 target mRNAs were predicted for 83 miRNAs. We identified 13 SNPs in the miRNA precursors. Among them, nine SNPs had the potential to either increase or decrease the energy of the predicted secondary structures of the precursors. Further, two SNPs in the 3' untranslated regions of target genes were predicted to either disturb or create miRNA-target interactions.

CONCLUSIONS

The common carp miRNAs and their target genes reported here will help further our understanding of the role of miRNAs in gene regulation. The analysis of the miRNA-related SNPs and their effects provided insights into the effects of SNPs on miRNA biogenesis and function. The resource data generated in this study will help advance the study of miRNA function and phenotype-associated miRNA identification.

Deep sequencing-based transcriptome analysis of chicken spleen in response to avian pathogenic Escherichia coli (APEC) infection

Avian pathogenic Escherichia coli (APEC) leads to economic losses in poultry production and is also a threat to human health. The goal of this study was to characterize the chicken spleen transcriptome and to identify candidate genes for response and resistance to APEC infection using Solexa sequencing. We obtained 14422935, 14104324, and 14954692 Solexa read pairs for non-challenged (NC), challenged-mild pathology (MD), and challenged-severe pathology (SV), respectively. A total of 148197 contigs and 98461 unigenes were assembled, of which 134949 contigs and 91890 unigenes match the chicken genome. In total, 12272 annotated unigenes take part in biological processes (11664), cellular components (11927), and molecular functions (11963). Summing three specific contrasts, 13650 significantly differentially expressed unigenes were found in NC Vs. MD (6844), NC Vs. SV (7764), and MD Vs. SV (2320). Some unigenes (e.g. CD148, CD45 and LCK) were involved in crucial pathways, such as the T cell receptor (TCR) signaling pathway and microbial metabolism in diverse environments. This study facilitates understanding of the genetic architecture of the chicken spleen transcriptome, and has identified candidate genes for host response to APEC infection.

miRNA regulation in the early development of barley seed

BACKGROUND

During the early stages of seed development many genes are under dynamic regulation to ensure the proper differentiation and establishment of the tissue that will constitute the mature grain. To investigate how miRNA regulation contributes to this process in barley, a combination of small RNA and mRNA degradome analyses were used to identify miRNAs and their targets.

RESULTS

Our analysis identified 84 known miRNAs and 7 new miRNAs together with 96 putative miRNA target genes regulated through a slicing mechanism in grain tissues during the first 15 days post anthesis. We also identified many potential miRNAs including several belonging to known miRNA families. Our data gave us evidence for an increase in miRNA-mediated regulation during the transition between pre-storage and storage phases. Potential miRNA targets were found in various signalling pathways including components of four phytohormone pathways (ABA, GA, auxin, ethylene) and the defence response to powdery mildew infection. Among the putative miRNA targets we identified were two essential genes controlling the GA response, a GA3oxidase1 and a homolog of the receptor GID1, and a homolog of the ACC oxidase which catalyses the last step of ethylene biosynthesis. We found that two MLA genes are potentially miRNA regulated, establishing a direct link between miRNAs and the R gene response.

CONCLUSION

Our dataset provides a useful source of information on miRNA regulation during the early development of cereal grains and our analysis suggests that miRNAs contribute to the control of development of the cereal grain, notably through the regulation of phytohormone response pathways.

Identification and characterization of microRNAs from barley (Hordeum vulgare L.) by high-throughput sequencing

MicroRNAs (miRNAs) are a class of endogenous RNAs that regulates the gene expression involved in various biological and metabolic processes. Barley is one of the most important cereal crops worldwide and is a model organism for genetic and genomic studies in Triticeae species. However, the miRNA research in barley has lagged behind other model species in grass family. To obtain more information of miRNA genes in barley, we sequenced a small RNA library created from a pool of equal amounts of RNA from four different tissues using Solexa sequencing. In addition to 126 conserved miRNAs (58 families), 133 novel miRNAs belonging to 50 families were identified from this sequence data set. The miRNA* sequences of 15 novel miRNAs were also discovered, suggesting the additional evidence for existence of these miRNAs. qRT-PCR was used to examine the expression pattern of six randomly selected miRNAs. Some miRNAs involved in drought and salt stress response were also identified. Furthermore, the potential targets of these putative miRNAs were predicted using the psRNATarget tools. Our results significantly increased the number of novel miRNAs in barley, which should be useful for further investigation into the biological functions and evolution of miRNAs in barley and other species.

Analysis of conserved microRNAs in floral tissues of sexual and apomictic Boechera species

Background

Apomixis or asexual seed formation represents a potentially important agronomic trait whose introduction into crop plants could be an effective way to fix and perpetuate a desirable genotype through successive seed generations. However, the gene regulatory pathways underlying apomixis remain unknown. In particular, the potential function of microRNAs, which are known to play crucial roles in many aspects of plant growth and development, remains to be determined with regards to the switch from sexual to apomictic reproduction.

Results

Using bioinformatics and microarray validation procedures, 51 miRNA families conserved among angiosperms were identified in Boechera. Microarray assay confirmed 15 of the miRNA families that were identified by bioinformatics techniques. 30 cDNA sequences representing 26 miRNAs could fold back into stable pre-miRNAs. 19 of these pre-miRNAs had miRNAs with Boechera-specific nucleotide substitutions (NSs). Analysis of the Gibbs free energy (ΔG) of these pre-miRNA stem-loops with NSs showed that the Boechera-specific miRNA NSs significantly (p ≤ 0.05) enhance the stability of stem-loops. Furthermore, six transcription factors, the Squamosa promoter binding protein like SPL6, SPL11 and SPL15, Myb domain protein 120 (MYB120), RELATED TO AP2.7 DNA binding (RAP2.7, TOE1 RAP2.7) and TCP family transcription factor 10 (TCP10) were found to be expressed in sexual or apomictic ovules. However, only SPL11 showed differential expression with significant (p ≤ 0.05) up-regulation at the megaspore mother cell (MMC) stage of ovule development in apomictic genotypes.

Conclusions

This study constitutes the first extensive insight into the conservation and expression of microRNAs in Boechera sexual and apomictic species. The miR156/157 target squamosa promoter binding protein-like 11 (SPL11) was found differentially expressed with significant (p ≤ 0.05) up-regulation at the MMC stage of ovule development in apomictic genotypes. The results also demonstrate that nucleotide changes in mature miRNAs significantly (p ≤ 0.05) enhance the thermodynamic stability of pre-miRNA stem-loops.

Discovery of barley miRNAs through deep sequencing of short reads

Background

MicroRNAs are important components of the regulatory network of biological systems and thousands have been discovered in both animals and plants. Systematic investigations performed in species with sequenced genomes such as Arabidopsis, rice, poplar and Brachypodium have provided insights into the evolutionary relationships of this class of small RNAs among plants. However, miRNAs from barley, one of the most important cereal crops, remain unknown.

Results

We performed a large scale study of barley miRNAs through deep sequencing of small RNAs extracted from leaves of two barley cultivars. By using the presence of miRNA precursor sequences in related genomes as one of a number of supporting criteria, we identified up to 100 miRNAs in barley. Of these only 56 have orthologs in wheat, rice or Brachypodium that are known to be expressed, while up to 44 appear to be specifically expressed in barley.

Conclusions

Our study, the first large scale investigation of small RNAs in barley, has identified up to 100 miRNAs. We demonstrate that reliable identification of miRNAs via deep sequencing in a species whose genome has not been sequenced requires a more careful analysis of sequencing errors than is commonly performed. We devised a read filtering procedure for dealing with errors. In addition, we found that the use of a large dataset of almost 35 million reads permits the use of read abundance distributions along putative precursor sequences as a practical tool for isolating miRNAs in a large background of reads originating from other non-coding and coding RNAs. This study therefore provides a generic approach for discovering novel miRNAs where no genome sequence is available.