Identification of miRNAs and Their Target Genes in Peach (Prunus persica L.) Using High-Throughput Sequencing and Degradome Analysis

MicroRNAs play critical roles in various biological and metabolic processes. The function of miRNAs has been widely studied in model plants such as Arabidopsis and rice. However, the number of identified miRNAs and related miRNA targets in peach (Prunus persica) is limited. To understand further the relationship between miRNAs and their target genes during tissue development in peach, a small RNA library and three degradome libraries were constructed from three tissues for deep sequencing. We identified 117 conserved miRNAs and 186 novel miRNA candidates in peach by deep sequencing and 19 conserved miRNAs and 13 novel miRNAs were further evaluated for their expression by RT-qPCR. The number of gene targets that were identified for 26 conserved miRNA families and 38 novel miRNA candidates, were 172 and 87, respectively. Some of the identified miRNA targets were abundantly represented as conserved miRNA targets in plant. However, some of them were first identified and showed important roles in peach development. Our study provides information concerning the regulatory network of miRNAs in peach and advances our understanding of miRNA functions during tissue development.     

MicroRNA-499 Expression Distinctively Correlates to Target Genes sox6 and rod1 Profiles to Resolve the Skeletal Muscle Phenotype in Nile Tilapia

A class of small non-coding RNAs, the microRNAs (miRNAs), has been shown to be essential for the regulation of specific cell pathways, including skeletal muscle development, maintenance and homeostasis in vertebrates. However, the relative contribution of miRNAs for determining the red and white muscle cell phenotypes is far from being fully comprehended. To better characterize the role of miRNA in skeletal muscle cell biology, we investigated muscle-specific miRNA (myomiR) signatures in Nile tilapia fish. Quantitative (RT-qPCR) and spatial (FISH) expression analyses revealed a highly differential expression (forty-four-fold) of miR-499 in red skeletal muscle compared to white skeletal muscle, whereas the remaining known myomiRs were equally expressed in both muscle cell types. Detailed examination of the miR-499 targets through bioinformatics led us to the sox6 and rod1 genes, which had low expression in red muscle cells according to RT-qPCR, FISH, and protein immunofluorescence profiling experiments. Interestingly, we verified that the high expression of miR-499 perfectly correlates with a low expression of sox6 and rod1 target genes, as verified by a distinctive predominance of mRNA destabilization and protein translational decay to these genes, respectively. Through a genome-wide comparative analysis of SOX6 and ROD1 protein domains and through an in silico gene regulatory network, we also demonstrate that both proteins are essentially similar in vertebrate genomes, suggesting their gene regulatory network may also be widely conserved. Overall, our data shed light on the potential regulation of targets by miR-499 associated with the slow-twitch muscle fiber type phenotype. Additionally the results provide novel insights into the evolutionary dynamics of miRNA and target genes enrolled in a putative constrained molecular pathway in the skeletal muscle cells of vertebrates.     

Selection of candidate reference genes for real-time PCR studies in lettuce under abiotic stresses

The process of selection and validation of reference genes is the first step in studies of gene expression by real-time quantitative polymerase chain reaction (RT-qPCR). The genome of lettuce, the most popular leaf vegetable cultivated worldwide, has recently been sequenced; therefore, suitable reference genes for reliable results in RT-qPCR analyses are required. In the present study, 17 candidate reference genes were selected, and their expression stability in lettuce leaves under drought, salt, heavy metal, and UV-C irradiation conditions and under the application of abscisic acid (ABA) was evaluated using geNorm and NormFinder software. The candidate reference genes included protein-coding traditional and novel reference genes and microRNAs (miRNAs). The results indicate that the expression stability is dependent on the experimental conditions. The novel protein-coding reference genes were more suitable than the traditional reference genes under drought, UV-C irradiation, and heavy metal conditions and under the application of ABA. Only under salinity conditions were the traditional protein-coding reference genes more stable than the novel genes. In addition, the miRNAs, mainly MIR169, MIR171/170 and MIR172, were stably expressed under the abiotic stresses evaluated, representing a suitable alternative approach for gene expression data normalization. The expression of phenylalanine ammonia lyase (PAL) and 4-hydroxyphenylpyruvate dioxygenase (HPPD) was used to further confirm the validated protein-coding reference genes, and the expression of MIR172 and MIR398 was used to confirm the validated miRNA genes, showing that the use of an inappropriate reference gene induces erroneous results. This work is the first survey of the stability of reference genes in lettuce and provides guidelines to obtain more accurate RT-qPCR results in lettuce studies.     

Dysregulated genes and miRNAs in the apoptosis pathway in colorectal cancer patients

Apoptosis is genetically regulated and involves intrinsic and extrinsic pathways. We examined 133 genes within these pathways to identify whether they are expressed differently in colorectal carcinoma (CRC) and normal tissue (N?=?217) and if they are associated with similar differential miRNA expression. Gene expression data (RNA-Seq) and miRNA expression data (Agilent Human miRNA Microarray V19.0) were generated. We focused on dysregulated genes with a fold change (FC) of >?1.50 or <?0.67, that were significant after adjustment for multiple comparisons. miRNA:mRNA seed-region matches were determined. Twenty-three genes were significantly downregulated (FC?<?0.67) and 18 were significantly upregulated (FC?>?1.50). Of these 41 genes, 11 were significantly associated with miRNA differential expression. BIRC5 had the greatest number of miRNA associations (14) and the most miRNAs with a seed-region match (10). Four of these matches, miR-145-5p, miR-150-5p, miR-195-5p, and miR-650, had a negative beta coefficient. CSF2RB was associated with ten total miRNAs (five with a seed-region match, and one miRNA, miR-92a-3p, with a negative beta coefficient). Of the three miRNAs associated with CTSS, miR-20b-5p, and miR-501-3p, had a seed-region match and a negative beta coefficient between miRNA:mRNA pairs. Several miRNAs that were associated with dysregulated gene expression, seed-region matches, and negative beta coefficients also were associated with CRC-specific survival. Our data suggest that miRNAs could influence several apoptosis-related genes. BIRC5, CTSS, and CSF2R all had seed-region matches with miRNAs that would favor apoptosis. Our study identifies several miRNA associated with apoptosis-related genes, that if validated, could be important therapeutic targets.     

The p53-signaling pathway and colorectal cancer: Interactions between downstream p53 target genes and miRNAs

IntroductionWe examined expression of genes in the p53-signaling pathway. We determine if genes that have significantly different expression in carcinoma tissue compared to normal mucosa also have significantly differentially expressed miRNAs. We utilize a sample of 217 CRC cases.MethodsWe focused on fold change (FC) > 1.50 or <0.67 for genes and miRNAs, that were statistically significant after adjustment for multiple comparisons. We evaluated the linear association between the differential expression of miRNA and mRNA. miRNA:mRNA seed-region matches also were determined.ResultsEleven dysregulated genes were associated with 37 dysregulated miRNAs; all were down-stream from the TP53 gene. MiR-150-5p (HR = 0.82) and miR-196b-5p (HR 0.73) significantly reduced the likelihood of dying from CRC when miRNA expression increased in rectal tumors.ConclusionsOur data suggest that activation of p53 from cellular stress, could target downstream genes that in turn could influence cell cycle arrest, apoptosis, and angiogenesis through mRNA:miRNA interactions.     

The use of miRNAs as reference genes for miRNA expression normalization during <Emphasis Type="Italic">Lilium</Emphasis> somatic embryogenesis by real-time reverse transcription PCR analysis

Expression profiling of miRNAs has the ability to reveal the essence of somatic embryogenesis (SE). qRT-PCR is one of the most commonly used techniques for dynamic miRNA detection but requires optimal reference genes for data reliability. This is the first report on reference gene validation for miRNA expression normalization in Lilium (Lilium pumilum DC. Fisch. and Lilium davidii var. unicolor). In this study, seventeen miRNAs together with two snRNAs (U4, U6), one rRNA (5S rRNA) and three protein-coding genes (FP, ACT, GAPDH) were selected as reference candidates, and their expression stability was validated by qRT-PCR among eleven developing SE cultures in two lilies. Four normalization algorithms, including geNorm, BestKeeper, NormFinder and RefFinder, were also used to evaluate the stability of the reference candidates. For Lilium pumilum DC. Fisch., lpu-miR159a was the optimal reference gene during SE, followed by lpu-miR408b, while U6 was the least stable reference candidate. For Lilium davidii var. unicolor, FP presented greater stability than did half of the miRNA candidates, but the best reference gene was lda-miR162, followed by lda-miR159a. Further analysis of the expression level of miR156 and miR529 was used to evaluate the validity of the reference genes in both lilies. In general, miRNAs are superior to common protein-coding genes and snRNAs / rRNAs as reference genes for miRNA expression normalization during Lilium SE, and the most suitable reference miRNA is different between two species in the same Lilium genus. This is a pioneer study using suitable miRNAs as reference genes in Lilium and constitutes a small but essential step for the further exploration of miRNA function in Lilium, thus offering valuable references for other plants.     

Temperature during early development has long-term effects on microRNA expression in Atlantic cod

Background

Environmental temperature has serious implications in life cycle of aquatic ectotherms. Understanding the molecular mechanisms of temperature acclimation and adaptation of marine organisms is of the uttermost importance for ecology, fisheries, and aquaculture, as it allows modeling the effects of global warming on population dynamics. Regulatory molecules are major modulators of acclimation and adaptation; among them, microRNAs (miRNAs) are versatile and substantial contributors to regulatory networks of development and adaptive plasticity. However, their role in thermal plasticity is poorly known. We have asked whether the temperature and its shift during the early ontogeny (embryonic and larval development) affect the miRNA repertoire of Atlantic cod (Gadus morhua), and if thermal experience has long-term consequences in the miRNA profile.

Results

We characterized miRNA during different developmental stages and in juvenile tissues using next generation sequencing. We identified 389 putative miRNA precursor loci, 120 novel precursor miRNAs, and 281 mature miRNAs. Some miRNAs showed stage- or tissue-enriched expression and miRNAs, such as the miR-17 ~ 92 cluster, myomiRs (miR-206), neuromiRs (miR-9, miR-124), miR-130b, and miR-430 showed differential expression in different temperature regimes. Long-term effect of embryonic incubation temperature was revealed on expression of some miRNAs in juvenile pituitary (miR-449), gonad (miR-27c, miR-30c, and miR-200a), and liver (let-7 h, miR-7a, miR-22, miR-34c, miR-132a, miR-192, miR-221, miR-451, miR-2188, and miR-7550), but not in brain. Some of differentially expressed miRNAs in the liver were confirmed using LNA-based rt-qPCR. The effect of temperature on methylation status of selected miRNA promoter regions was mostly inconclusive.

Conclusions

Temperature elevation by several degrees during embryonic and larval developmental stages significantly alters the miRNA profile, both short-term and long-term. Our results suggest that a further rise in seas temperature might affect life history of Atlantic cod.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1503-7) contains supplementary material, which is available to authorized users.     

Deep Sequencing of Small RNA Repertoires in Mice Reveals Metabolic Disorders-Associated Hepatic miRNAs

Obesity and associated metabolic disorders contribute importantly to the metabolic syndrome. On the other hand, microRNAs (miRNAs) are a class of small non-coding RNAs that repress target gene expression by inducing mRNA degradation and/or translation repression. Dysregulation of specific miRNAs in obesity may influence energy metabolism and cause insulin resistance, which leads to dyslipidemia, steatosis hepatis and type 2 diabetes. In the present study, we comprehensively analyzed and validated dysregulated miRNAs in ob/ob mouse liver, as well as miRNA groups based on miRNA gene cluster and gene family by using deep sequencing miRNA datasets. We found that over 13.8% of the total analyzed miRNAs were dysregulated, of which 37 miRNA species showed significantly differential expression. Further RT-qPCR analysis in some selected miRNAs validated the similar expression patterns observed in deep sequencing. Interestingly, we found that miRNA gene cluster and family always showed consistent dysregulation patterns in ob/ob mouse liver, although they had various enrichment levels. Functional enrichment analysis revealed the versatile physiological roles (over six signal pathways and five human diseases) of these miRNAs. Biological studies indicated that overexpression of miR-126 or inhibition of miR-24 in AML-12 cells attenuated free fatty acids-induced fat accumulation. Taken together, our data strongly suggest that obesity and metabolic disturbance are tightly associated with functional miRNAs. We also identified hepatic miRNA candidates serving as potential biomarkers for the diagnose of the metabolic syndrome.