MicroRNA target prediction by expression analysis of host genes

Gene regulatory networks in disease and ageing

The precise control of gene expression is required for the maintenance of cellular homeostasis and proper cellular function, and the declining control of gene expression with age is considered a major contributor to age-associated changes in cellular physiology and disease. The coordination of gene expression can be represented through models of the molecular interactions that govern gene expression levels, so-called gene regulatory networks. Gene regulatory networks can represent interactions that occur through signal transduction, those that involve regulatory transcription factors, or statistical models of gene-gene relationships based on the premise that certain sets of genes tend to be coexpressed across a range of conditions and cell types. Advances in experimental and computational technologies have enabled the inference of these networks on an unprecedented scale and at unprecedented precision. Here, we delineate different types of gene regulatory networks and their cell-biological interpretation. We describe methods for inferring such networks from large-scale, multi-omics datasets and present applications that have aided our understanding of cellular ageing and disease mechanisms.

A computational modeling of pri-miRNA expression

MicroRNAs (miRNAs) play crucial roles in gene regulation. Most studies focus on mature miRNAs, which leaves many unknowns about primary miRNAs (pri-miRNAs). To fill the gap, we attempted to model the expression of pri-miRNAs in 1829 primary cell types, cell lines, and tissues in this study. We demonstrated that the expression of pri-miRNAs can be modeled well by the expression of specific sets of mRNAs, which we termed their associated mRNAs. These associated mRNAs differ from their corresponding target mRNAs and are enriched with specific functions. Most associated mRNAs of a miRNA are shared across conditions, while on average, about one-fifth of the associated mRNAs are condition-specific. Our study shed new light on understanding miRNA biogenesis and general gene transcriptional regulation.

Exosomal MicroRNA and Protein Profiles of Hepatitis B Virus-Related Hepatocellular Carcinoma Cells

Infection with hepatitis B virus (HBV) is a main risk factor for hepatocellular carcinoma (HCC). Extracellular vesicles, such as exosomes, play an important role in tumor development and metastasis, including regulation of HBV-related HCC. In this study, we have characterized exosome microRNA and proteins released in vitro from hepatitis B virus (HBV)-related HCC cell lines SNU-423 and SNU-182 and immortalized normal hepatocyte cell lines (THLE2 and THLE3) using microRNA sequencing and mass spectrometry. Bioinformatics, including functional enrichment and network analysis, combined with survival analysis using data related to HCC in The Cancer Genome Atlas (TCGA) database, were applied to examine the prognostic significance of the results. More than 40 microRNAs and 200 proteins were significantly dysregulated (p < 0.05) in the exosomes released from HCC cells in comparison with the normal liver cells. The functional analysis of the differentially expressed exosomal miRNAs (i.e., mir-483, mir-133a, mir-34a, mir-155, mir-183, mir-182), their predicted targets, and exosomal differentially expressed proteins (i.e., POSTN, STAM, EXOC8, SNX9, COL1A2, IDH1, FN1) showed correlation with pathways associated with HBV, virus activity and invasion, exosome formation and adhesion, and exogenous protein binding. The results from this study may help in our understanding of the role of HBV infection in the development of HCC and in the development of new targets for treatment or non-invasive predictive biomarkers of HCC.

Vitamin D, a Secosteroid Hormone and Its Multifunctional Receptor, Vitamin D Receptor, in Alzheimer's Type Neurodegeneration

Vitamin D is a secosteroid hormone exerting neurosteroid-like properties. Its well-known nuclear hormone receptor, and recently proposed as a mitochondrial transcription factor, vitamin D receptor, acts for its primary functions. The second receptor is an endoplasmic reticulum protein, protein disulfide isomerase A3 (PDIA3), suggested to act as a rapid response. Vitamin D has effects on various systems, particularly through calcium metabolism. Among them, the nervous system has an important place in the context of our subject. Recent studies have shown that vitamin D and its receptors have numerous effects on the nervous system. Neurodegeneration is a long-term process. Throughout a human life span, so is vitamin D deficiency. Our previous studies and others have suggested that the out-come of long-term vitamin D deficiency (hypovitaminosis D or inefficient utilization of vitamin D), may lead neurons to be vulnerable to aging and neurodegeneration. We suggest that keeping vitamin D levels at adequate levels at all stages of life, considering new approaches such as agonists that can activate vitamin D receptors, and utilizing other derivatives produced in the synthesis process with UVB are crucial when considering vitamin D-based intervention studies. Given most aspects of vitamin D, this review outlines how vitamin D and its receptors work and are involved in neurodegeneration, emphasizing Alzheimer's disease.

microRNAs: Critical Players during Helminth Infections

microRNAs (miRNAs) are a group of small non-coding RNAs that regulate gene expression post-transcriptionally through their interaction with the 3' untranslated regions (3' UTR) of target mRNAs, affecting their stability and/or translation. Therefore, miRNAs regulate biological processes such as signal transduction, cell death, autophagy, metabolism, development, cellular proliferation, and differentiation. Dysregulated expression of microRNAs is associated with infectious diseases, where miRNAs modulate important aspects of the parasite-host interaction. Helminths are parasitic worms that cause various neglected tropical diseases affecting millions worldwide. These parasites have sophisticated mechanisms that give them a surprising immunomodulatory capacity favoring parasite persistence and establishment of infection. In this review, we analyze miRNAs in infections caused by helminths, emphasizing their role in immune regulation and its implication in diagnosis, prognosis, and the development of therapeutic strategies.

Ability of a Combined FIB4/miRNA181a Score to Predict Significant Liver Fibrosis in NAFLD Patients

Liver biopsy is the gold standard for assessing fibrosis, but there is a need to seek non-invasive biomarkers for this purpose. The aim of this study was to evaluate the correlation between the serum levels of the microRNAs miR-21, miR-29a, miR-122, miR-155 and miR-181a and the phenotypic expression of NAFLD. A cross-sectional study was carried out on 108 NAFLD patients diagnosed by liver biopsy. FIB-4 and NAFLD fibrosis scores were calculated. The comparison between the distributions of microRNA values according to the presence or absence of histological fibrosis (F2–F4) was performed. A multivariate logistic regression analysis was performed to build a score for predicting fibrosis using FIB-4 and Ln (miR-181a) as independent variables. Only miR-181a showed a statistical difference between patients with significant liver fibrosis (>F2) and those without (F0–F1) (p = 0.017). FIB-4 revealed an AUC on the ROC curve of 0.667 to predict clinically significant fibrosis (F2–F4). When assessed using the score in association with Ln (miR-181a), there was an improvement in the ROC curve, with an AUC of 0.71. miR-181a can be used as a non-invasive method of predicting fibrosis in NAFLD, and an association with FIB-4 has the potential to increase the accuracy of each method alone.

Pan-Cancer Analysis Reveals Common and Specific Relationships between Intragenic miRNAs and Their Host Genes

MicroRNAs (miRNAs) are small endogenous non-coding RNAs that play important roles in regulating gene expression. Most miRNAs are located within or close to genes (host). miRNAs and their host genes have either coordinated or independent transcription. We performed a comprehensive investigation on co-transcriptional patterns of miRNAs and host genes based on 4707 patients across 21 cancer types. We found that only 11.6% of miRNA-host pairs were co-transcribed consistently and strongly across cancer types. Most miRNA-host pairs showed a strong coexpression only in some specific cancer types, demonstrating a high heterogenous pattern. For two particular types of intergenic miRNAs, readthrough and divergent miRNAs, readthrough miRNAs showed higher coexpression with their host genes than divergent ones. miRNAs located within non-coding genes had tighter co-transcription with their hosts than those located within protein-coding genes, especially exonic and junction miRNAs. A few precursor miRNAs changed their dominate form between 5' and 3' strands in different cancer types, including miR-486, miR-99b, let-7e, miR-125a, let-7g, miR-339, miR-26a, miR-16, and miR-218, whereas only two miRNAs with multiple host genes switched their co-transcriptional partner in different cancer types (miR-219a-1 with SLC39A7/HSD17B8 and miR-3615 with RAB37/SLC9A3R1). miRNAs generated from distinct precursors (such as miR-125b from miR-125b-1 or miR-125b-2) were more likely to have cancer-dependent main contributors. miRNAs and hosts were less co-expressed in KIRC than other cancer types, possibly due to its frequent VHL mutations. Our findings shed new light on miRNA biogenesis and cancer diagnosis and treatments.

TFEB Signalling-Related MicroRNAs and Autophagy

The oncogenic Transcription Factor EB (TFEB), a member of MITF-TFE family, is known to be the most important regulator of the transcription of genes responsible for the control of lysosomal biogenesis and functions, autophagy, and vesicles flux. TFEB activation occurs in response to stress factors such as nutrient and growth factor deficiency, hypoxia, lysosomal stress, and mitochondrial damage. To reach the final functional status, TFEB is regulated in multimodal ways, including transcriptional rate, post-transcriptional regulation, and post-translational modifications. Post-transcriptional regulation is in part mediated by miRNAs. miRNAs have been linked to many cellular processes involved both in physiology and pathology, such as cell migration, proliferation, differentiation, and apoptosis. miRNAs also play a significant role in autophagy, which exerts a crucial role in cell behaviour during stress or survival responses. In particular, several miRNAs directly recognise TFEB transcript or indirectly regulate its function by targeting accessory molecules or enzymes involved in its post-translational modifications. Moreover, the transcriptional programs triggered by TFEB may be influenced by the miRNA-mediated regulation of TFEB targets. Finally, recent important studies indicate that the transcription of many miRNAs is regulated by TFEB itself. In this review, we describe the interplay between miRNAs with TFEB and focus on how these types of crosstalk affect TFEB activation and cellular functions.

The interaction of Akkermansia muciniphila with host-derived substances, bacteria and diets

Abstract

Trillions of microbes inhabit the human gut and build extremely complex communities. Gut microbes contribute to host metabolisms for better or worse and are widely studied and associated with health and disease. Akkermansia muciniphila is a gut microbiota member, which uses mucin as both carbon and nitrogen sources. Many studies on A. muciniphila have been conducted since this unique bacterium was first described in 2004. A. muciniphila can play an important role in our health because of its beneficial effects, such as improving type II diabetes and obesity and anti-inflammation. A. muciniphila establishes its position as a next-generation probiotic. Besides the effect of A. muciniphila on host health, a technique for boosting has been investigated. In this review, we show what factors can modulate the abundance of A. muciniphila focusing on the interaction with host-derived substances, other bacteria and diets. This review also refers to the possibility of the interaction between medicine and A. muciniphila; this will open up future treatment strategies that can increase A. muciniphila abundance in the gut.

Key points

• Host-derived substances such as bile, microRNA and melatonin as well as mucin have beneficial effects on A. muciniphila.

• Gut and probiotic bacteria and diet ingredients such as carbohydrates and phytochemicals could boost the abundance of A. muciniphila.

• Several medicines could affect the growth of A. muciniphila.

Telomerase reverse transcriptase activates transcription of miR500A to inhibit Hedgehog signalling and promote cell invasiveness

Telomerase reverse transcriptase (TERT) maintains telomere homeostasis, thus ensuring chromosome stability and cell proliferation. In addition, several telomere-independent functions of human TERT have been described. In this study, we report that TERT binds directly to the TCF binding elements located upstream of the oncomiR miR500A, and induces its transcription. This function was independent of the telomerase activity, as shown with experiments using catalytically inactive TERT and inhibitors of TERT and the TERT RNA component. miR500A was in turn found to target three key components of the Hedgehog signalling pathway: Patched 1; Gli family zinc finger 3; and Cullin 3, thereby promoting cell invasion. Our results point to the crucial role of the TERT-miR500A-Hedgehog axis in tumour aggressiveness and highlight the therapeutic potential of targeting noncanonical TERT functions in cancer.

microRNA-seq of cartilage reveals an overabundance of miR-140-3p which contains functional isomiRs

miR-140 is selectively expressed in cartilage. Deletion of the entire Mir140 locus in mice results in growth retardation and early-onset osteoarthritis-like pathology; however, the relative contribution of miR-140-5p or miR-140-3p to the phenotype remains to be determined. An unbiased small RNA sequencing approach identified miR-140-3p as significantly more abundant (>10-fold) than miR-140-5p in human cartilage. Analysis of these data identified multiple miR-140-3p isomiRs differing from the miRBase annotation at both the 5' and 3' end, with >99% having one of two seed sequences (5' bases 2-8). Canonical (miR-140-3p.2) and shifted (miR-140-3p.1) seed isomiRs were overexpressed in chondrocytes and transcriptomics performed to identify targets. miR-140-3p.1 and miR-140-3p.2 significantly down-regulated 694 and 238 genes, respectively, of which only 162 genes were commonly down-regulated. IsomiR targets were validated using 3'UTR luciferase assays. miR-140-3p.1 targets were enriched within up-regulated genes in rib chondrocytes of Mir140-null mice and within down-regulated genes during human chondrogenesis. Finally, through imputing the expression of miR-140 from the expression of the host gene WWP2 in 124 previously published data sets, an inverse correlation with miR-140-3p.1 predicted targets was identified. Together these data suggest the novel seed containing isomiR miR-140-3p.1 is more functional than original consensus miR-140-3p seed containing isomiR.

Pten and Dicer1 loss in the mouse uterus causes poorly differentiated endometrial adenocarcinoma

Endometrial cancer remains the most common gynecological malignancy in the United States. While the loss of the tumor suppressor, PTEN (phosphatase and tensin homolog), is well studied in endometrial cancer, recent studies suggest that DICER1, the endoribonuclease responsible for miRNA genesis, also plays a significant role in endometrial adenocarcinoma. Conditional uterine deletion of Dicer1 and Pten in mice resulted in poorly differentiated endometrial adenocarcinomas, which expressed Napsin A and HNF1B (hepatocyte nuclear factor 1 homeobox B), markers of clear-cell adenocarcinoma. Adenocarcinomas were hormone-independent. Treatment with progesterone did not mitigate poorly differentiated adenocarcinoma, nor did it affect adnexal metastasis. Transcriptomic analyses of DICER1 deleted uteri or Ishikawa cells revealed unique transcriptomic profiles and global miRNA downregulation. Computational integration of miRNA with mRNA targets revealed deregulated let-7 and miR-16 target genes, similar to published human DICER1-mutant endometrial cancers from TCGA (The Cancer Genome Atlas). Similar to human endometrial cancers, tumors exhibited dysregulation of ephrin-receptor signaling and transforming growth factor-beta signaling pathways. LIM kinase 2 (LIMK2), an essential molecule in p21 signal transduction, was significantly upregulated and represents a novel mechanism for hormone-independent pathogenesis of endometrial adenocarcinoma. This preclinical mouse model represents the first genetically engineered mouse model of poorly differentiated endometrial adenocarcinoma.

Computational prediction of miRNA/mRNA duplexomes at the whole human genome scale reveals functional subnetworks of interacting genes with embedded miRNA annealing motifs

Perfect annealing between microRNAs (miRNAs) and messenger RNAs (mRNAs) was computationally searched at a broad scale in the human genome to determine whether theoretical pairing is restrictively represented in functional subnetworks or is randomly distributed. Massive RNA interference (RNAi) pairing motifs in genes constitute a remarkable subnetwork that displays highly genetically and biochemically interconnected genes. These analyses show unexpected repertoires of genes defined by their congruence in comatching with miRNAs at numerous sites and by their interconnection based on protein/protein interactions or proteins regulating the activity of others. This offers insights into the putatively coregulated homeostasis of large networks of genes by RNAi, whereas other networks seem to be independent of this regulatory mode. Genes accordingly defined by theoretical RNAi pairing cluster mainly in subnetworks related to cellular, metabolic and developmental processes and their regulation. Indeed, genes harboring numerous potential sites of hybridization with miRNAs are highly enriched with GO terms depicting the abovementioned processes and are grouped in a subnetwork of genes that are significantly more highly connected than they would be according to a random distribution. The significant number of interacting genes that present numerous potential comatches with miRNAs suggests that they may be under the control of the integrative and concerted action of multiple miRNAs.

miRNA target identification and prediction as a function of time in gene expression data

The understanding of miRNA target interactions is still limited due to conflicting data and the fact that high-quality validation of targets is a time-consuming process. Faster methods like high-throughput screens and bioinformatics predictions are employed but suffer from several problems. One of these, namely the potential occurrence of downstream (i.e. secondary) effects in high-throughput screens has been only little discussed so far. However, such effects limit usage for both the identification of interactions and for the training of bioinformatics tools. In order to analyse this problem more closely, we performed time-dependent microarray screening experiments overexpressing human miR-517a-3p, and, together with published time-dependent datasets of human miR-17-5p, miR-135b and miR-124 overexpression, we analysed the dynamics of deregulated genes. We show that the number of deregulated targets increases over time, whereas seed sequence content and performance of several miRNA target prediction algorithms actually decrease over time. Bioinformatics recognition success of validated miR-17 targets was comparable to that of data gained only 12 h post-transfection. We therefore argue that the timing of microarray experiments is of critical importance for detecting direct targets with high confidence and for the usability of these data for the training of bioinformatics prediction tools.

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.

Identification of key microRNAs involved in tumorigenesis and prognostic microRNAs in breast cancer

Breast cancer is a commonly diagnosed cancer in women, and one of the leading causes of cancer-related death among female patients However, the key microRNAs involved in its tumorigenesis and microRNAs of prognostic values have not been fully understood. In the present study, we aimed to perform a systematic analysis of microRNA expression profiles to identify some key microRNAs associated with tumor initiation and prognosis. Using TCGA breast cancer datasets, we identified 110 differentially expressed microRNAs. The functional enrichment analysis of the upregulated microRNAs revealed signaling transduction pathways, such as Notch and Wnt signaling pathway, and metabolism-related pathways such as sugar and nucleotide sugar metabolism, and oxidative stress response. Moreover, multivariable Cox model based on three variables of hsa-mir-130a, hsa-mir-3677, and hsa-mir-1247 stratified patients into high-risk and low-risk groups, which showed significant prognostic difference. In addition, we also tested the performance of this model in patient cohorts of any specific breast cancer subtypes or different TNM stages. The high performance in risk prediction was also observed in all of breast cancer subtypes and TNM stages. We also observed that there were highly possible interactions between hsa-mir-130a and seven target genes. Among these target genes, VAV3 and ESR1 were predicted as the target genes of hsa-mir-130a, suggesting that hsa-mir-130a may function by regulating the expression of VAV3 and ESR1 in breast cancer. In conclusion, the stratification based on the multivariable Cox model showed high performance in risk prediction. The dysregulated microRNAs and prognostic microRNAs greatly improved our understanding of the microRNA-related molecular mechanism underlying breast cancer.

AAV-miR-204 Protects from Retinal Degeneration by Attenuation of Microglia Activation and Photoreceptor Cell Death

Inherited retinal diseases (IRDs) represent a frequent cause of genetic blindness. Their high genetic heterogeneity hinders the application of gene-specific therapies to the vast majority of patients. We recently demonstrated that the microRNA miR-204 is essential for retinal function, although the underlying molecular mechanisms remain poorly understood. Here, we investigated the therapeutic potential of miR-204 in IRDs. We subretinally delivered an adeno-associated viral (AAV) vector carrying the miR-204 precursor to two genetically different IRD mouse models. The administration of AAV-miR-204 preserved retinal function in a mouse model for a dominant form of retinitis pigmentosa (RHO-P347S). This was associated with a reduction of apoptotic photoreceptor cells and with a better preservation of photoreceptor marker expression. Transcriptome analysis showed that miR-204 shifts expression profiles of transgenic retinas toward those of healthy retinas by the downregulation of microglia activation and photoreceptor cell death. Delivery of miR-204 exerted neuroprotective effects also in a mouse model of Leber congenital amaurosis, due to mutations of the Aipl1 gene. Our study highlights the mutation-independent therapeutic potential of AAV-miR204 in slowing down retinal degeneration in IRDs and unveils the previously unreported role of this miRNA in attenuating microglia activation and photoreceptor cell death.

Modulation and Evolution of Animal Development through microRNA Regulation of Gene Expression

microRNAs regulate gene expression by blocking the translation of mRNAs and/or promoting their degradation. They, therefore, play important roles in gene regulatory networks (GRNs) by modulating the expression levels of specific genes and can tune GRN outputs more broadly as part of feedback loops. These roles for microRNAs provide developmental buffering on one hand but can facilitate evolution of development on the other. Here we review how microRNAs can modulate GRNs during animal development as part of feedback loops and through their individual or combinatorial targeting of multiple different genes in the same network. We then explore how changes in the expression of microRNAs and consequently targets can facilitate changes in GRNs that alter development and lead to phenotypic evolution. The reviewed studies exemplify the key roles played by microRNAs in the regulation and evolution of gene expression during developmental processes in animals.

Long noncoding RNA: a crosslink in biological regulatory network

Long noncoding RNAs (lncRNAs) had been defined as a novel class of functional RNAs longer than 200 nucleotides around a decade ago. It is widely acknowledged that lncRNAs play a significant role in regulation of gene expression, but the biological and molecular mechanisms are diverse and complex, and remain to be determined. Especially, the regulatory network of lncRNAs associated with other biological molecules is still a controversial matter, thus becoming a new frontier of the studies on transcriptome. Recent advance in high-throughput sequencing technologies and bioinformatics approaches may be an accelerator to lift the mysterious veil. In this review, we will outline well-known associations between lncRNAs and other biological molecules, demonstrate the diverse bioinformatics approaches applied in prediction and analysis of lncRNA interaction and perform a case study for lncRNA linc00460 to concretely decipher the lncRNA regulatory network.

Pan-cancer analysis on microRNA-associated gene activation

Background

While microRNAs (miRNAs) were widely considered to repress target genes at mRNA and/or protein levels, emerging evidence from in vitro experiments has shown that miRNAs can also activate gene expression in particular contexts. However, this counterintuitive observation has rarely been reported or interpreted in in vivo conditions.

Methods

We systematically explored the positive correlation between miRNA and gene expressions and its potential implications in tumorigenesis, based on 8375 patient samples across 31 major human cancers from The Cancer Genome Atlas (TCGA).

Findings

We found that positive miRNA-gene correlations are surprisingly prevalent and consistent across cancer types, and show distinct patterns than negative correlations. The top-ranked positive correlations are significantly involved in the immune cell differentiation and cell membrane signaling related processes, and display strong power in stratifying patients in terms of survival rate. Although intragenic miRNAs generally tend to co-express with their host genes, a substantial portion of miRNAs shows no obvious correlation with their host gene plausibly due to non-conservation. A miRNA can upregulate a gene by inhibiting its upstream suppressor, or shares transcription factors with that gene, both leading to positive correlation. The miRNA/gene sites associated with the top-ranked positive correlations are more likely to form super-enhancers compared to randomly chosen pairs. Wet-lab experiments revealed that positive correlations partially remain in in vitro condition.

Interpretation

Our study brings new insights into the critical role of miRNA in gene regulation and the complex mechanisms underlying miRNA functions, and reveals both biological and clinical significance of miRNA-associated gene activation.

Integrative microRNA-mRNA Analysis of Muscle Tissues in Qianhua Mutton Merino and Small Tail Han Sheep Reveals Key Roles for oar-miR-655-3p and oar-miR-381-5p

The Qianhua Mutton Merino (QHMM) is a new variety of sheep (Ovis aries) with improved meat performance compared with the traditional Small Tail Han (STH) sheep variety. We recently reported the transcriptome profiling of longissimus muscle tissues between QHMM and STH sheep. In the present study, we aimed to evaluate key micro (mi)RNA-mRNA networks associated with sheep muscle growth and development. We used miRNA sequencing to obtain longissimus muscle miRNA profiles from QHMM and STH sheep. We identified a total of 153 known sheep miRNAs, of which 4 were differentially expressed (DE) between the 2 sheep varieties. We combined these results with mRNA library data to build an miRNA-mRNA network, including 26 target genes of the 4 DE miRNAs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that 26 target genes were significantly enriched in 86 biological processes, including muscle organogenesis, myoblast migration, cell proliferation, and adipose tissue development, and in 9 metabolic pathways, including carbohydrate, nucleotide, and amino acid metabolic pathways. oar-miR-655-3p and its target gene ACSM3 and oar-miR-381-5p and its target gene ABAT were selected for subsequent analysis based on GO and KEGG analyses. The binding sites of oar-miR-655-3p with ACSM3 and oar-miR-381-5p with ABAT were validated by a dual-luciferase reporter gene detection system. This represents the first integrative analysis of miRNA-mRNA networks in QHMM and STH muscles and suggests that DE miRNAs, especially oar-miR-655-3p and oar-miR-381-5p, play crucial roles in muscle growth and development.

Measurement of Urinary Level of a Specific Competing endogenous RNA network (FOS and RCAN mRNA/ miR-324-5p, miR-4738-3p, /lncRNA miR-497-HG) Enables Diagnosis of Bladder Cancer

OBJECTIVE

To assist in the diagnosis, treatment, and prognostic prediction of bladder cancer, the molecular patterns associated with it should be elucidated. Competing endogenous RNA network: MicroRNA (miRNA), long noncoding RNA (lncRNA), and their target autophagy genes have been strongly implicated in tumor development and metastasis.

PATIENTS AND METHODS

Bioinformatics analysis was performed to retrieve a ceRNA: lncRNA-miRNA-mRNA network linked to autophagy and relevant to bladder cancer. Expression of selected noncoding human RNAs (miR-324-5p, miR-4738-3p, and lncRNA miR-497-HG) and their target genes (RCAN1 mRNA and FOSB mRNA) was examined by qPCR in bladder tissues and urine samples obtained from 196 individuals (98 patients with bladder cancer, 48 patients with benign lesions, and 50 healthy controls).

RESULTS

Expression levels of the selected genes in urine samples in the bladder cancer group were significantly different from those in the control group (P < 0.001). Expression in bladder cancer tissue samples correlated with that in urine samples. Urinary expression levels of all biomarkers had high accuracy to distinguish patients with and without bladder cancer, with FOSB mRNA and RCAN1 mRNA having the highest accuracy (99% for RCAN1 mRNA or FOSB mRNA, 87.8% for miR-324-5p, 84.7% for miR-4738-3p, and 90.5% for lncRNA miR-497-HG). FOSB mRNA and RCAN1 mRNA expression showed also a higher accuracy than cytology (77.6%).

CONCLUSION

The significant differential expression of the ceRNA network: lncRNA-miRNA-mRNA network in bladder cancer as compared to noncancerous controls has revealed the superior accuracy of the chosen biomarkers to cytology, especially FOSB mRNA and RCAN1 mRNA, suggesting their involvement in bladder cancer pathogenesis and promising role for future diagnosis, and targeted therapy.

Single-cell microRNA-mRNA co-sequencing reveals non-genetic heterogeneity and mechanisms of microRNA regulation

Measuring multiple omics profiles from the same single cell opens up the opportunity to decode molecular regulation that underlies intercellular heterogeneity in development and disease. Here, we present co-sequencing of microRNAs and mRNAs in the same single cell using a half-cell genomics approach. This method demonstrates good robustness (~95% success rate) and reproducibility (R2 = 0.93 for both microRNAs and mRNAs), yielding paired half-cell microRNA and mRNA profiles, which we can independently validate. By linking the level of microRNAs to the expression of predicted target mRNAs across 19 single cells that are phenotypically identical, we observe that the predicted targets are significantly anti-correlated with the variation of abundantly expressed microRNAs. This suggests that microRNA expression variability alone may lead to non-genetic cell-to-cell heterogeneity. Genome-scale analysis of paired microRNA-mRNA co-profiles further allows us to derive and validate regulatory relationships of cellular pathways controlling microRNA expression and intercellular variability.

A miRNA Host Response Signature Accurately Discriminates Acute Respiratory Infection Etiologies

Background: Acute respiratory infections (ARIs) are the leading indication for antibacterial prescriptions despite a viral etiology in the majority of cases. The lack of available diagnostics to discriminate viral and bacterial etiologies contributes to this discordance. Recent efforts have focused on the host response as a source for novel diagnostic targets although none have explored the ability of host-derived microRNAs (miRNA) to discriminate between these etiologies. Methods: In this study, we compared host-derived miRNAs and mRNAs from human H3N2 influenza challenge subjects to those from patients with Streptococcus pneumoniae pneumonia. Sparse logistic regression models were used to generate miRNA signatures diagnostic of ARI etiologies. Generalized linear modeling of mRNAs to identify differentially expressed (DE) genes allowed analysis of potential miRNA:mRNA relationships. High likelihood miRNA:mRNA interactions were examined using binding target prediction and negative correlation to further explore potential changes in pathway regulation in response to infection. Results: The resultant miRNA signatures were highly accurate in discriminating ARI etiologies. Mean accuracy was 100% [88.8-100; 95% Confidence Interval (CI)] in discriminating the healthy state from S. pneumoniae pneumonia and 91.3% (72.0-98.9; 95% CI) in discriminating S. pneumoniae pneumonia from influenza infection. Subsequent differential mRNA gene expression analysis revealed alterations in regulatory networks consistent with known biology including immune cell activation and host response to viral infection. Negative correlation network analysis of miRNA:mRNA interactions revealed connections to pathways with known immunobiology such as interferon regulation and MAP kinase signaling. Conclusion: We have developed novel human host-response miRNA signatures for bacterial and viral ARI etiologies. miRNA host response signatures reveal accurate discrimination between S. pneumoniae pneumonia and influenza etiologies for ARI and integrated analyses of the host-pathogen interface are consistent with expected biology. These results highlight the differential miRNA host response to bacterial and viral etiologies of ARI, offering new opportunities to distinguish these entities.

Amyloid Beta 1-42 Alters the Expression of miRNAs in Cortical Neurons

Recently, Aβ1-42 was demonstrated to have the potential to translocate into the nucleus and to be involved in the transcriptional regulation of certain neurodegeneration-related genes. This data raises the question of whether Aβ-induced neurodegeneration might include the expression of miRNAs. Thus, our aim in this study was to investigate the effects of Aβ1-42 on certain miRNAs which are related with vitamin D metabolism, neuronal differentiation, development, and memory. This question was investigated in primary cortical neurons that were treated with 10 μM Aβ and/or 10^-8 M 1,25-dihydroxyvitamin D3 at different time points by expression analysis of let-7a-5p, miR-26b-5p, miR-27b-3p, miR-31a-5p, miR-125b-5p, and miR-192-5p with qRT-PCR. Our data indicate that amyloid pathology has effects on the expression of miRNAs. Furthermore, some of these miRNAs simultaneously regulate the proteins or the enzymes involved in neuronal metabolism. The experimental setup that we used and the data we acquired supply valuable information about the miRNAs that play a part in the Aβ pathology and suggested Aβ as a counterpart of vitamin D at the crossroads of neuronal differentiation, development, and memory.

Microrna-26b attenuates monocrotaline-induced pulmonary vascular remodeling via targeting connective tissue growth factor (CTGF) and cyclin D1 (CCND1)

MicroRNAs are involved in the control of cell growth, and deregulated pulmonary artery smooth muscle cell proliferation plays an essential role in the development of pulmonary hypertension. The objective of this study was to identify differentially expressed microRNA(s) and explore its therapeutic role in treatment of the disease. MicroRNA expression profile analysis showed microRNA-26b was differentially expressed in pulmonary artery smooth muscle cells harvested from monocrotaline-treated rats, and we validated microRNA-26b targets, in vitro and in vivo, CTGF and CCND1, both of which have been shown, in our previous work, to be involved in the pathogenesis of pulmonary hypertension. In vivo experiments demonstrated monocrotaline-induced pulmonary artery remodeling could be almost completely abolished by administration of microRNA-26b, while CTGF or CCND1 shRNA significantly, but only partially, attenuated the remodeling by silencing the designed target. Additionally, exogenous expression of the microRNA-26b substantially downregulated CTGF and CCND1 in human pulmonary artery smooth muscle cells. MicroRNA-26b might be a potent therapeutic tool to treat pulmonary hypertension.

Integrated micro/messenger RNA regulatory networks in essential thrombocytosis

Essential thrombocytosis (ET) is a chronic myeloproliferative disorder with an unregulated surplus of platelets. Complications of ET include stroke, heart attack, and formation of blood clots. Although platelet-enhancing mutations have been identified in ET cohorts, genetic networks causally implicated in thrombotic risk remain unestablished. In this study, we aim to identify novel ET-related miRNA-mRNA regulatory networks through comparisons of transcriptomes between healthy controls and ET patients. Four network discovery algorithms have been employed, including (a) Pearson correlation network, (b) sparse supervised canonical correlation analysis (sSCCA), (c) sparse partial correlation network analysis (SPACE), and, (d) (sparse) Bayesian network analysis-all through a combined data-driven and knowledge-based analysis. The result predicts a close relationship between an 8-miRNA set (miR-9, miR-490-5p, miR-490-3p, miR-182, miR-34a, miR-196b, miR-34b*, miR-181a-2*) and a 9-mRNA set (CAV2, LAPTM4B, TIMP1, PKIG, WASF1, MMP1, ERVH-4, NME4, HSD17B12). The majority of the identified variables have been linked to hematologic functions by a number of studies. Furthermore, it is observed that the selected mRNAs are highly relevant to ET disease, and provide an initial framework for dissecting both platelet-enhancing and functional consequences of dysregulated platelet production.

A distinct microRNA expression profile is associated with α[11C]-methyl-L-tryptophan (AMT) PET uptake in epileptogenic cortical tubers resected from patients with tuberous sclerosis complex

Tuberous sclerosis complex (TSC) is characterized by hamartomatous lesions in various organs and arises due to mutations in the TSC1 or TSC2 genes. TSC mutations lead to a range of neurological manifestations including epilepsy, cognitive impairment, autism spectrum disorders (ASD), and brain lesions that include cortical tubers. There is evidence that seizures arise at or near cortical tubers, but it is unknown why some tubers are epileptogenic while others are not. We have previously reported increased tryptophan metabolism measured with α[11C]-methyl-l-tryptophan (AMT) positron emission tomography (PET) in epileptogenic tubers in approximately two-thirds of patients with tuberous sclerosis and intractable epilepsy. However, the underlying mechanisms leading to seizure onset in TSC remain poorly characterized. MicroRNAs are enriched in the brain and play important roles in neurodevelopment and brain function. Recent reports have shown aberrant microRNA expression in epilepsy and TSC. In this study, we performed microRNA expression profiling in brain specimens obtained from TSC patients undergoing epilepsy surgery for intractable epilepsy. Typically, in these resections several non-seizure onset tubers are resected together with the seizure-onset tubers because of their proximity. We directly compared seizure onset tubers, with and without increased tryptophan metabolism measured with PET, and non-onset tubers to assess the role of microRNAs in epileptogenesis associated with these lesions. Whether a particular tuber was epileptogenic or non-epileptogenic was determined with intracranial electrocorticography, and tryptophan metabolism was measured with AMT PET. We identified a set of five microRNAs (miR-142-3p, 142-5p, 223-3p, 200b-3p and 32-5p) that collectively distinguish among the three primary groups of tubers: non-onset/AMT-cold (NC), onset/AMT-cold (OC), and onset/AMT-hot (OH). These microRNAs were significantly upregulated in OH tubers compared to the other two groups, and microRNA expression was most significantly associated with AMT-PET uptake. The microRNAs target a group of genes enriched for synaptic signaling and epilepsy risk, including SLC12A5, SYT1, GRIN2A, GRIN2B, KCNB1, SCN2A, TSC1, and MEF2C. We confirmed the interaction between miR-32-5p and SLC12A5 using a luciferase reporter assay. Our findings provide a new avenue for subsequent mechanistic studies of tuber epileptogenesis in TSC.

A Looking-Glass of Non-coding RNAs in oral cancer

Oral cancer is a multifactorial pathology and is characterized by the lack of efficient treatment and accurate diagnostic tools. This is mainly due the late diagnosis; therefore, reliable biomarkers for the timely detection of the disease and patient stratification are required. Non-coding RNAs (ncRNAs) are key elements in the physiological and pathological processes of various cancers, which is also reflected in oral cancer development and progression. A better understanding of their role could give a more thorough perspective on the future treatment options for this cancer type. This review offers a glimpse into the ncRNA involvement in oral cancer, which can help the medical community tap into the world of ncRNAs and lay the ground for more powerful diagnostic, prognostic and treatment tools for oral cancer that will ultimately help build a brighter future for these patients.

Stem Cell-Derived Exosomes, Autophagy, Extracellular Matrix Turnover, and miRNAs in Cardiac Regeneration during Stem Cell Therapy

Stem cell therapy (SCT) raises the hope for cardiac regeneration in ischemic hearts. However, underlying molecular mechanisms for repair of dead myocardium by SCT in the ischemic heart is poorly understood. Growing evidences suggest that cardiac matrix stiffness and differential expressions of miRNAs play a crucial role in stem cell survival and differentiation. However, their roles on transplanted stem cells, for myocardial repair of the ischemic heart, remain unclear. Transplanted stem cells may act in an autocrine and/or paracrine manner to regenerate the dead myocardium. Paracrine mediators such as stem cell-derived exosomes are emerging as a novel therapeutic strategy to overcome some of the limitations of SCT. These exosomes carry microRNAs (miRNAs) that may regulate stem cell differentiation into a specific lineage. MicroRNAs may also contribute to stiffness of surrounding matrix by regulating extracellular matrix (ECM) turnover. The survival of transplanted stem cell depends on its autophagic process that maintains cellular homeostasis. Therefore, exosomes, miRNAs, extracellular matrix turnover, and autophagy may have an integral role in improving the efficacy of SCT. This review elaborates the specific roles of these regulatory components on cardiac regeneration in the ischemic heart during SCT.

MicroRNAs and regulated interaction networks reveal differences between adult and pediatric acute myeloid leukemia

OBJECTIVES

The purpose of this study was to identify featured microRNAs and their regulated network between adult and pediatric acute myeloid leukemia (AML) and find potential utility as biomarkers for diagnosis and treatment of pediatric AML.

METHODS

We downloaded the microRNA expression dataset GSE35320 from Gene Expression Omnibus database and selected expression chips from bone marrow of 71 pediatric AML samples and 6 adulthood AML samples. Differentially expressed microRNAs were identified by Wilcox test. The target genes of these microRNAs were predicted using an integrative method and their functional enrichment analysis was performed using DAVID. Finally, STRING database and Cytoscape software was used to construct and analyze the interaction network.

RESULTS

A total of 7 differentially expressed microRNAs were identified and the remarkably up-regulated and down-regulated microRNAs were miR-16 and miR-142-5p which included 323 and 22 predicted target genes, respectively. The target genes of 7 microRNAs were most associated with regulation of cell cycle, p53 signaling pathway, Wnt signaling pathway and neurotrophin signaling pathway. The interaction network of miR-16 target genes was constructed among 354 high confidence interaction pairs. The core genes of the network, such as TP53, BCL2, VEGFA, had a role in prognosis of children with AML.

CONCLUSIONS

The featured microRNAs and their target genes are significant in the occurrence and development of pediatric AML, which is likely to be important for the identification of therapeutic targets and biomarkers for these patients.

Hypoxia-induced microRNA-26b inhibition contributes to hypoxic pulmonary hypertension via CTGF

The objective of this study was to explore the role of miRNAs in the control of HPH as well as molecular mechanism underlying. Computational analysis and luciferase assay were carried out to search the target gene of miR-26b. Luciferase assay, RT-PCR and western-blot analysis was performed to test interaction among hypoxia, miR-26b, SRF and CTGF. MiR-26b was significantly downregulated; meanwhile, CTGF and SRF were significantly upregulated in HPH rat model. Using computational analysis, CTGF was found to be a virtual target gene of miR-26b, and only cell transfected with vectors containing wild-type CTGF 3'UTR and miR-26b showed a lower luciferase activity than scramble control. Hypoxia significantly inhibited miR-26b promoter, and promoted SRF promoter. Meanwhile, hypoxia had no effect on CTGF promoter. In addition, SRF promoted the promoter of CTGF. MiR-26b was significantly downregulated; meanwhile, CTGF and SRF were upregulated in PASMCs exposed to hypoxia. In addition, miR-26b and SRF siRNA, but not CTGF siRNA, significantly inhibited SRF expression. Meanwhile, miR-26b, SRF siRNA, and CTGF siRNA significantly inhibited CTGF expression in hypoxia-treated cell. PASMCs treated with hypoxia showed higher cell viability and higher percentage cells in S phase than the control, which could be reversed by miR-26b, SRF siRNA, and CTGF siRNA transfection. These findings suggested that hypoxia induced miR-26b inhibition and SRF and CTGF upregulation in HPH rat model. CTGF mediated hypoxia-induced regulation of miR-26b and SRF in proliferation of PASMCs, which indicated that hypoxia-induced miR-26b inhibition contributed to the pathogenesis of HPH via CTGF.

Biomarker MicroRNAs for Diagnosis of Oral Squamous Cell Carcinoma Identified Based on Gene Expression Data and MicroRNA-mRNA Network Analysis

Oral squamous cell carcinoma is one of the most malignant tumors with high mortality rate worldwide. Biomarker discovery is critical for early diagnosis and precision treatment of this disease. MicroRNAs are small noncoding RNA molecules which often regulate essential biological processes and are good candidates for biomarkers. By integrative analysis of both the cancer-associated gene expression data and microRNA-mRNA network, miR-148b-3p, miR-629-3p, miR-27a-3p, and miR-142-3p were screened as novel diagnostic biomarkers for oral squamous cell carcinoma based on their unique regulatory abilities in the network structure of the conditional microRNA-mRNA network and their important functions. These findings were confirmed by literature verification and functional enrichment analysis. Future experimental validation is expected for the further investigation of their molecular mechanisms.

The Tyr113His T/C rs1051740 and 'very slow' phenotype of the EPHX1 gene alters miR-26b-5p and miR-1207-5p expression in pregnancy

BACKGROUND

Environmental insults and microsomal epoxide hydrolase 1 (EPHX1) single nucleotide polymorphisms (SNPs), Tyr113His T/C rs1051740 and His139Arg A/G rs2234922, aberrantly alters microRNA (miR) expression and are linked to low birthweights (LBW).

OBJECTIVES

To investigate the interplay between pollution, EPHX1 SNPs and miRs during pregnancy and associated LBW outcomes.

METHODS

South African pregnant women (n=241) were recruited in the MACE birth cohort study in Durban, a city with high levels of industry and traffic related pollutants. EPHX1 SNPs were genotyped using PCR-RFLP and grouped into their respective phenotypes, i.e. normal (N), slow (S), very slow (VS) and fast (F). EPHX1, miR-26b-5p, miR-193b-3p and miR-1207-5p expression were determined using quantitative PCR.

RESULTS

Mothers with the Tyr113His SNP had low iron levels [TT vs. TC+CC: mean difference (MD)=0.67g/dl; p=0.0167], LBW [TT vs. TC+CC: MD=189.30g; p=0.0067], and low EPHX1 expression; p<0.0001. miR-26b-5p and miR-1207-5p expression were significantly higher in the CC genotypes compared to TT+TC groups; p<0.0001. The opposite trend occurred for miR-193b-3p; p=0.0045. Mothers with the VS phenotype had low iron levels [N vs. VS and VS vs. F: MD=2.03 and -1.96g/dl; p=0.0021, respectively], decreased gestational age [VS vs. F: MD=-2.14weeks; p=0.0051, respectively], and LBW [N vs. VS, VS vs. F and S vs. VS: MD=1000, -940.30 and 968.80g; p<0.0001, respectively]; F phenotype had the highest EPHX1 expression [N vs. F, VS vs. F and S vs. F: MD=-1.067, -1.854 and -1.379; p=0.0002, respectively]; and N phenotype had low miR-26b-5p [N vs. VS: MD=-0.6100; p=0.0159] and miR-1207-5p [N vs. VS and VS vs. F: MD=-0.834 and 1.103; p=0.0007, respectively] expression. miR-193b-3p expression between phenotypes remained unchanged.

CONCLUSION

The Tyr113His T/C variant of rs1051740 and VS phenotype alters EPHX1, miR-26b-5p and miR-1207-5p expression, and contributes towards low blood iron levels and LBW.

Targeted regulation of MiR-98 on E2F1 increases chemosensitivity of leukemia cells K562/A02

Background

miRNA is a microRNA that negatively regulates protein expression at post-transcriptional or translational level. It is widely involved in the pathogenesis of tumors. miR-98 belongs to the let-7 family, and its overexpression can increase the sensitivity to drugs in solid cancer cells. However, the function of miR-98 in leukemia is still unclear. In this study, the effect of miR-98 on drug resistance and proliferation of leukemia cells were investigated.

Methods

Real-time quantitative polymerase chain reaction analyzed the expression difference between miR-98 and E2F1 in leukemia cell lines, K562 and K562/A02. The downstream target gene of miR-98 was predicted by TargetScan; K562/A02 was transiently transfected with miR-98 mimic to upregulate the expression of miR-98; real-time quantitative polymerase chain reaction and Western blot were used to analyze the expression alterations of E2F1; cell counting kit-8 was used to evaluate the influence on K562/A02 proliferation and sensitivity to chemotherapeutic drugs; meanwhile, Western blot was used to analyze the expression of p21, Bax, matrix metalloproteinase 9 and ABCG2 proteins.

Results

E2F1 is one of the target genes of miR-98 proved by bioinformatics. Compared with the K562, the level of miRNA-98 expression was decreased in K562/A02, but the level of E2F1 expression was upregulated. Leukemia cell line K562/A02 was transfected with miR-98 mimic to upregulate the expression of miR-98, the expression of E2F1 was significantly decreased. After upregulating the miR-98 expression in K562/A02, the proliferation was weakened, and the sensitivity to chemotherapy was increased. Western blot showed that upregulated miR-98 expression increased the levels of p21 and BAX proteins in K562/A02 cells, and decreased the levels of matrix metalloprotease 9 and ABCG2 proteins, which were significantly different compared with those before miR-98 mimic transfection.

Conclusion

In the leukemia drug-resistant cell line K562/A02, the targeted upregulated expression of miR-98 could decrease the proliferation of leukemia cells and improve the sensitivity to chemotherapeutics by inhibiting E2F1 expression. miR-98 might be a potential target for overcoming leukemia multidrug resistance.

miRNAs target databases: developmental methods and target identification techniques with functional annotations

PURPOSE

microRNA (miRNA) regulates diverse biological mechanisms and metabolisms in plants and animals. Thus, the discoveries of miRNA has revolutionized the life sciences and medical research.The miRNA represses and cleaves the targeted mRNA by binding perfect or near perfect or imperfect complementary base pairs by RNA-induced silencing complex (RISC) formation during biogenesis process. One miRNA interacts with one or more mRNA genes and vice versa, hence takes part in causing various diseases. In this paper, the different microRNA target databases and their functional annotations developed by various researchers have been reviewed. The concurrent research review aims at comprehending the significance of miRNA and presenting the existing status of annotated miRNA target resources built by researchers henceforth discovering the knowledge for diagnosis and prognosis.

METHODS AND RESULTS

This review discusses the applications and developmental methodologies for constructing target database as well as the utility of user interface design. An integrated architecture is drawn and a graphically comparative study of present status of miRNA targets in diverse diseases and various biological processes is performed. These databases comprise of information such as miRNA target-associated disease, transcription factor binding sites (TFBSs) in miRNA genomic locations, polymorphism in miRNA target, A-to-I edited target, Gene Ontology (GO), genome annotations, KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, target expression analysis, TF-miRNA and miRNA-mRNA interaction networks, drugs-targets interactions, etc.

CONCLUSION

miRNA target databases contain diverse experimentally and computationally predicted target through various algorithms. The comparison of various miRNA target database has been performed on various parameters. The computationally predicted target databases suffer from false positive information as there is no common theory for prediction of miRNA targets. The review conclusion emphasizes the need of more intelligent computational improvement for the miRNA target identification, their functional annotations and datasbase development.

Gene-microRNA network module analysis for ovarian cancer

Background

MicroRNAs (miRNAs) are involved in many biological processes by regulating post-transcriptional gene expression. The alterations of the regulatory pathways can cause different diseases including cancer. Although many works have been done to study the gene-miRNA regulatory network, the intertwined relationship is far from being fully understood. The objective of this study is to integrate both gene expression data and miRNA data so as to explore the complex relationships among them.

Methods

By integrating the networks consisting of gene coexpression, miRNA coexpression, gene-miRNA coexpression, and the known gene-miRNA interactions, we aim to find the most connected network modules so as to study their functions and properties. In this paper, we proposed an optimization model for identification of the modules in the integrated networks. This model tries to find both the modules in the gene-gene and miRNA-miRNA coexpression networks and the densely connected gene-miRNA subneworks. An approximation computational method was developed to solve the optimization problem.

Results

We applied the method to 556 human ovarian cancer samples with both gene expression data and miRNA expression data. The identified modules are significantly enriched by miRNA clusters, GO-BPs, and KEGG pathways. We compared our method with some existing methods and showed the better performance of our method. We also showed that the miRNAs and genes in our identified modules are associated with cancers, especially ovarian cancer.

Conclusions

This study provides strong support that the subnetworks consisting of genes and miRNAs with close interactions contribute the cancers. The proposed computational method can be applied to other studies that are related to different types of networks.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0357-1) contains supplementary material, which is available to authorized users.

YM500v3: a database for small RNA sequencing in human cancer research

We previously presented the YM500 database, which contains >8000 small RNA sequencing (smRNA-seq) data sets and integrated analysis results for various cancer miRNome studies. In the updated YM500v3 database (http://ngs.ym.edu.tw/ym500/) presented herein, we not only focus on miRNAs but also on other functional small non-coding RNAs (sncRNAs), such as PIWI-interacting RNAs (piRNAs), tRNA-derived fragments (tRFs), small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs). There is growing knowledge of the role of sncRNAs in gene regulation and tumorigenesis. We have also incorporated >10 000 cancer-related RNA-seq and >3000 more smRNA-seq data sets into the YM500v3 database. Furthermore, there are two main new sections, ‘Survival' and ‘Cancer', in this updated version. The ‘Survival’ section provides the survival analysis results in all cancer types or in a user-defined group of samples for a specific sncRNA. The ‘Cancer’ section provides the results of differential expression analyses, miRNA–gene interactions and cancer miRNA-related pathways. In the ‘Expression’ section, sncRNA expression profiles across cancer and sample types are newly provided. Cancer-related sncRNAs hold potential for both biotech applications and basic research.

Regulation of Connexins Expression Levels by MicroRNAs, an Update

Control of cell-cell coordination and communication is regulated by several factors, including paracrine and autocrine release of biomolecules, and direct exchange of soluble factors between cells through gap junction channels. Additionally, hemichannels also participate in cell-cell coordination through the release of signaling molecules, such as ATP and glutamate. A family of transmembrane proteins named connexins forms both gap junction channels and hemichannels. Because of their importance in cell and tissue coordination, connexins are controlled both by post-translational and post-transcriptional modifications. In recent years, non-coding RNAs have garnered research interest due to their ability to exert post-transcriptional regulation of gene expression. One of the most recent, well-documented control mechanisms of protein synthesis is found through the action of small, single-stranded RNA, called micro RNAs (miRNAs or miRs). Put simply, miRNAs are negative regulators of the expression of a myriad proteins involved in many physiological and pathological processes. This mini review will briefly summarize what is currently known about the action of miRNAs over Cxs expression/function in different organs under some relevant physiological and pathological conditions.

Regulatory network analysis of genes and microRNAs in human hepatoblastoma

Hepatoblastoma (HB) is a common type of primary tumor in children. Previous studies have examined the expression of genes, including transcription factors (TFs), target genes, host genes and microRNAs (miRNAs or miRs) associated with HB. However, the regulatory pathways of miRNAs and genes remain unclear. In the present study, a novel perspective is proposed, which focuses on HB and the associated regulatory pathways, to construct three networks at various levels, including a differentially expressed network, an associated network and a global network. Genes and miRNAs are considered as key factors in the network. In the three networks, the associations between each pair of factors, including TFs that regulate miRNAs, miRNAs that interact with target genes and miRNAs that are located at host genes, were analyzed. The differentially expressed network is considered to be the most crucial of the three networks. All factors in the differentially expressed network were mutated or differentially expressed, which indicated that the majority of the factors were cancerogenic factors that may lead to HB. In addition, the network contained numerous abnormal linkages that may trigger HB. If the expression of each factor was corrected to a normal level, HB may be successfully treated. The associated network included more HB-associated genes and miRNAs, and was useful for analyzing the pathogenesis of HB. By analyzing these close associations, the first and the last factor of the regulatory pathways were revealed to have important roles in HB. For example, v-myc avian myelocytomatosis viral oncogene neuroblastoma derived homolog (MYCN) was observed to regulate Homo sapiens (hsa)-miR-221, hsa-miR-18a and hsa-miR-17-5p, but no miRNAs targeted MYCN. In conclusion, the pathways and mechanisms underlying HB were expounded in the present study, which proposed a fundamental hypothesis for additional studies.

Discovering miRNA Regulatory Networks in Holt-Oram Syndrome Using a Zebrafish Model

MicroRNAs (miRNAs) are small non-coding RNAs that play an important role in the post-transcriptional regulation of gene expression. miRNAs are involved in the regulation of many biological processes such as differentiation, apoptosis, and cell proliferation. miRNAs are expressed in embryonic, postnatal, and adult hearts, and they have a key role in the regulation of gene expression during cardiovascular development and disease. Aberrant expression of miRNAs is associated with abnormal cardiac cell differentiation and dysfunction. Tbx5 is a member of the T-box gene family, which acts as transcription factor involved in the vertebrate heart development. Alteration of Tbx5 level affects the expression of hundreds of genes. Haploinsufficiency and gene duplication of Tbx5 are at the basis of the cardiac abnormalities associated with Holt–Oram syndrome (HOS). Recent data indicate that miRNAs might be an important part of the regulatory circuit through which Tbx5 controls heart development. Using high-throughput technologies, we characterized genome-widely the miRNA and mRNA expression profiles in WT- and Tbx5-depleted zebrafish embryos at two crucial developmental time points, 24 and 48 h post fertilization (hpf). We found that several miRNAs, which are potential effectors of Tbx5, are differentially expressed; some of them are already known to be involved in cardiac development and functions, such as miR-30, miR-34, miR-190, and miR-21. We performed an integrated analysis of miRNA expression data with gene expression profiles to refine computational target prediction approaches by means of the inversely correlation of miRNA–mRNA expressions, and we highlighted targets, which have roles in cardiac contractility, cardiomyocyte proliferation/apoptosis, and morphogenesis, crucial functions regulated by Tbx5. This approach allowed to discover complex regulatory circuits involving novel miRNAs and protein coding genes not considered before in the HOS such as miR-34a and miR-30 and their targets.

Tiny giants of gene regulation: experimental strategies for microRNA functional studies

The discovery over two decades ago of short regulatory microRNAs (miRNAs) has led to the inception of a vast biomedical research field dedicated to understanding these powerful orchestrators of gene expression. Here we aim to provide a comprehensive overview of the methods and techniques underpinning the experimental pipeline employed for exploratory miRNA studies in animals. Some of the greatest challenges in this field have been uncovering the identity of miRNA-target interactions and deciphering their significance with regard to particular physiological or pathological processes. These endeavors relied almost exclusively on the development of powerful research tools encompassing novel bioinformatics pipelines, high-throughput target identification platforms, and functional target validation methodologies. Thus, in an unparalleled manner, the biomedical technology revolution unceasingly enhanced and refined our ability to dissect miRNA regulatory networks and understand their roles in vivo in the context of cells and organisms. Recurring motifs of target recognition have led to the creation of a large number of multifactorial bioinformatics analysis platforms, which have proved instrumental in guiding experimental miRNA studies. Subsequently, the need for discovery of miRNA-target binding events in vivo drove the emergence of a slew of high-throughput multiplex strategies, which now provide a viable prospect for elucidating genome-wide miRNA-target binding maps in a variety of cell types and tissues. Finally, deciphering the functional relevance of miRNA post-transcriptional gene silencing under physiological conditions, prompted the evolution of a host of technologies enabling systemic manipulation of miRNA homeostasis as well as high-precision interference with their direct, endogenous targets. For further resources related to this article, please visit the WIREs website.

De novo assembly of the blunt snout bream (Megalobrama amblycephala) gill transcriptome to identify ammonia exposure associated microRNAs and their targets

De novo transcriptome sequencing is a robust method for microRNA (miRNA) target gene prediction, especially for organisms without reference genomes. Following exposure of Megalobrama amblycephala to ammonia (0.1 or 20 mg L⁻¹ ), two cDNA libraries were constructed from the fish gills and sequenced using Illumina HiSeq 2000. Over 90 million reads were generated and de novo assembled into 46, 615 unigenes, which were then extensively annotated by comparing to different protein databases, followed by biochemical pathway prediction. The expression of 2666 unigenes significantly differed; 1961 were up-regulated, while 975 were down-regulated. Among these, 250 unigenes were identified as the targets for 10 conserved and 4 putative novel miRNA families by miRNA target computational prediction. We examined expression of ssa-miRNA-21 and its target genes by real-time quantitative PCR and found agreement with the sequencing data. This study demonstrates the feasibility of identifying miRNA targets by transcriptome analysis. The transcriptome assembly data represent a substantial increase in the genomic resources available for Megalobrama amblycephala and will be useful for gene expression profile analysis and miRNA functional annotation.

•

Several differentially expressed genes and miRNA were identified.

•

Target prediction indicated that a number of miRNAs involved in immunity.

•

The study provides an understanding of molecular mechanisms of ammonia-induced toxicology in fish.

Comparison of Prognostic MicroRNA Biomarkers in Blood and Tissues for Gastric Cancer

Gastric cancer (GC) still keeps up high mortality worldwide with poor prognosis. Efficient and non-invasive prognostic biomarkers are urgently needed. MicroRNAs are non-coding RNAs playing roles in post-transcriptional gene regulation, which contribute to various biological processes such as development, differentiation and carcinogenesis. MicroRNA expression profiles have been associated with the prognosis and outcome in GC. MicroRNA prognostic biomarkers have been identified from blood or tissues samples, but with different prognostic features. Understanding the various roles of microRNAs in different sample sources of GC will provide deep insights into GC progression. In this review, we highlight the distinct prognostic roles of microRNAs biomarkers in blood and tissue according to their relationships with prognostic parameters, survival rates and target pathways. This will be useful for non-invasive biomarker development and selection in prognosis of GC.

MicroRNAs involved in the browning process of adipocytes

The present review focuses on the role of miRNAs in the control of white adipose tissue browning, a process which describes the recruitment of adipocytes showing features of brown adipocytes in white adipose tissue. MicroRNAs (miRNAs) are a class of short non-coding RNAs (19-22 nucleotides) involved in gene regulation. Although the main effect of miRNAs is the inhibition of the translational machinery, thereby preventing the production of the protein product, the activation of protein translation has also been described in the literature. In addition to modifying translation, miRNAs binding to its target mRNAs also trigger the recruitment and association of mRNA decay factors, leading to mRNA destabilization, degradation, and thus to the decrease in expression levels. Although a great number of miRNAs have been reported to potentially regulate genes that play important roles in the browning process, only a reduced number of studies have demonstrated experimentally an effect on this process associated to changes in miRNA expressions, so far. These studies have shown, by using either primary adipocyte cultures or experimental models of mice (KO mice, mice overexpressing a specific miRNA), that miR-196a, miR-26, and miR-30 are needed for browning process development. By contrast, miR-155, miR-133, miR-27b, and miR-34 act as negative regulators of this process [corrected]. Further studies are needed to fully describe the miRNA network-involved white adipose tissue browning regulation.

CancerNet: a database for decoding multilevel molecular interactions across diverse cancer types

Protein–protein interactions (PPIs) and microRNA (miRNA)–target interactions are important for deciphering the mechanisms of tumorigenesis. However, current PPI databases do not support cancer-specific analysis. Also, no available databases can be used to retrieve cancer-associated miRNA–target interactions. As the pathogenesis of human cancers is affected by several miRNAs rather than a single miRNA, it is needed to uncover miRNA synergism in a systems level. Here for each cancer type, we constructed a miRNA–miRNA functionally synergistic network based on the functions of miRNA targets and their topological features in that cancer PPI network. And for the first time, we report the cancer-specific database CancerNet (http://bis.zju.edu.cn/CancerNet), which contains information about PPIs, miRNA–target interactions and functionally synergistic miRNA–miRNA pairs across 33 human cancer types. In addition, PPI information across 33 main normal tissues and cell types are included. Flexible query methods are allowed to retrieve cancer molecular interactions. Network viewer can be used to visualize interactions that users are interested in. Enrichment analysis tool was designed to detect significantly overrepresented Gene Ontology categories of miRNA targets. Thus, CancerNet serves as a comprehensive platform for assessing the roles of proteins and miRNAs, as well as their interactions across human cancers.

Integrated analysis of miRNA and mRNA paired expression profiling of prenatal skeletal muscle development in three genotype pigs

MicroRNAs (miRNAs) play a vital role in muscle development by binding to messenger RNAs (mRNAs). Based on prenatal skeletal muscle at 33, 65 and 90 days post-coitus (dpc) from Landrace, Tongcheng and Wuzhishan pigs, we carried out integrated analysis of miRNA and mRNA expression profiling. We identified 33, 18 and 67 differentially expressed miRNAs and 290, 91 and 502 mRNA targets in Landrace, Tongcheng and Wuzhishan pigs, respectively. Subsequently, 12 mRNAs and 3 miRNAs differentially expressed were validated using quantitative real-time PCR (qPCR), and 5 predicted miRNA targets were confirmed via dual luciferase reporter or western blot assays. We identified a set of miRNAs and mRNA genes differentially expressed in muscle development. Gene ontology (GO) enrichment analysis suggests that the miRNA targets are primarily involved in muscle contraction, muscle development and negative regulation of cell proliferation. Our data indicated that more mRNAs are regulated by miRNAs at earlier stages than at later stages of muscle development. Landrace and Tongcheng pigs also had longer phases of myoblast proliferation than Wuzhishan pigs. This study will be helpful to further explore miRNA-mRNA interactions in myogenesis and aid to uncover the molecular mechanisms of muscle development and phenotype variance in pigs.