miRmapper: A Tool for Interpretation of miRNA⁻mRNA Interaction Networks

Hsa_circular RNA_0045474 Facilitates Osteoarthritis Via Modulating microRNA-485-3p and Augmenting Transcription Factor 4

Circular RNA (circRNA) influences on the pathological process of osteoarthritis (OA) and may be a potential marker for disease diagnosis. The study was to scrutinize the association of circ_0045474 with OA. Clinical samples of OA patients were collected, and 12 circRNAs derived from KPNA2 gene were examined. CHON-001 cells were stimulated with IL-1β to construct an OA chondrocyte model. miR-485-3p, transcription factor 4 (TCF4) and circ_0045474, type II procollagen (COL2A1), and human collagenase-3 (MMP13) were tested. Furthermore, cell activities were analyzed. The relationship between miR-485-3p, TCF4, and circ_0045474 was determined. The role of circ_0045474 in vivo was further confirmed by constructing an OA mouse model by anterior cruciate ligament transection. circ_0045474 expression was elevated in OA patients. Suppressing circ_0045474 restrained IL-1β-stimulated extracellular matrix degradation, inflammatory cytokine secretion, and chondrocyte apoptosis. Circ_0045474 competitively combined with miR-485-3p, while TCF4 was the target of miR-485-3p. Circ_0045474 modulated IL-1β-stimulated extracellular matrix degradation, inflammatory cytokine secretion, and chondrocyte apoptosis via miR-485-3p/TCF4 axis. Suppressing circ 0045474 was effective to alleviate OA in mice. Silenced circ_0045474 suppresses OA progression in vitro and vivo via miR-485-3p/TCF4 axis. In short, circ_0045474 can be considered a novel therapeutic target for OA.

Transcriptome-wide analysis uncovers regulatory elements of the antennal transcriptome repertoire of bumblebee at different life stages

Bumblebees are crucial pollinators, providing essential ecosystem services and global food production. The success of pollination services relies on the interaction between sensory organs and the environment. The antenna functions as a versatile multi-sensory organ, pivotal in mediating chemosensory/olfactory information, and governs adaptive responses to environmental changes. Despite an increasing number of RNA-sequencing studies on insect antenna, comprehensive antennal transcriptome studies at the different life stages were not elucidated systematically. Here, we quantified the expression profile and dynamics of coding/microRNA genes of larval head and antennal tissues from early- and late-stage pupa to the adult of Bombus terrestris as suitable model organism among pollinators. We further performed Pearson correlation analyses on the gene expression profiles of the antennal transcriptome from larval head tissue to adult stages, exploring both positive and negative expression trends. The positively correlated coding genes were primarily enriched in sensory perception of chemical stimuli, ion transport, transmembrane transport processes and olfactory receptor activity. Negatively correlated genes were mainly enriched in organic substance biosynthesis and regulatory mechanisms underlying larval body patterning and the formation of juvenile antennal structures. As post-transcriptional regulators, miR-1000-5p, miR-13b-3p, miR-263-5p and miR-252-5p showed positive correlations, whereas miR-315-5p, miR-92b-3p, miR-137-3p, miR-11-3p and miR-10-3p exhibited negative correlations in antennal tissue. Notably, based on the inverse expression relationship, positively and negatively correlated microRNA (miRNA)-mRNA target pairs revealed that differentially expressed miRNAs predictively targeted genes involved in antennal development, shaping antennal structures and regulating antenna-specific functions. Our data serve as a foundation for understanding stage-specific antennal transcriptomes and large-scale comparative analysis of transcriptomes in different insects.

Curcumin modulates cell type-specific miRNA networks to induce cytotoxicity in ovarian cancer cells

AIM

To understand the epigenetic role of curcumin, a natural polyphenolic compound extracted from the spice Curcuma longa in inducing cytotoxicity in two molecularly distinct ovarian cancer cell lines: PA1 and A2780.

MATERIALS AND METHODS

An integrated mRNA-miRNA sequence analysis was performed to determine the curcumin-induced mRNA-miRNA regulatory networks in the induction of cytotoxicity. The miRNA-mRNA pathways, the miRNAs and their targets implicated in apoptosis, autophagy, DNA damage, and stemness markers were validated. Gene/miRNA expressions were validated using qPCR and protein expressions by western blotting. Curcumin-induced oncogenic /tumor-suppressor miRNAs were profiled utilising the oncomiRdb database. Similarly, the expressions of oncogenes/tumor suppressor genes were profiled and correlated with the TCGA ovarian cancer dataset. A dual luciferase assay was performed to investigate the interaction of miR-199a-5p to its direct target, DDR1.

KEY FINDINGS

The expression of several miRNAs demonstrated an inverse correlation with their respective direct targets. In curcumin-treated PA1 cells, miR-335-5p target ATG5 (autophagic), and OCT4 (pluripotent gene) were downregulated, miR-32a target PTEN (tumor suppressor) was upregulated, miR-1285 target P53 (tumor suppressor) was upregulated, and both miR-182-5p and miR-503-3p target BCL2, were down-regulated. Contrastingly, in curcumin-treated A2780 cells, miR-181a-3p target ATG5, miR-30a-5p, and miR-216a target BECN1 (autophagic) were upregulated, and miR-129a-5p target BCL2 were downregulated. The reversal of the oncomiR/TSmiR profile revealed suppression of oncogenic processes by curcumin. Curcumin treatment induced a moderate cisplatin-sensitisation effect and impaired epithelial-to-mesenchymal transition (EMT) characteristics. Curcumin also regulated the miR-199a-5p/DDR1 axis with a decrease in collagen deposition.

SIGNIFICANCE

The activity of curcumin is cell-type specific. Distinct miRNA regulatory networks were activated to induce multiple modes of cellular cytotoxicity in these ovarian cancer cells. This study further highlights the molecular mechanism of curcumin action in ovarian cancers establishing its candidacy as a promising drug candidate.

Blood miRNA levels associated with ADHD traits in children across six European birth cohorts

BACKGROUND

Attention-deficit/hyperactivity disorder (ADHD) is a prevalent and highly heritable neurodevelopmental disorder of major societal concern. Diagnosis can be challenging and there are large knowledge gaps regarding its etiology, though studies suggest an interplay of genetic and environmental factors involving epigenetic mechanisms. MicroRNAs (miRNAs) show promise as biomarkers of human pathology and novel therapies, and here we aimed to identify blood miRNAs associated with traits of ADHD as possible biomarker candidates and further explore their biological relevance.

METHODS

Our study population consisted of 1126 children (aged 5-12 years, 46% female) from the Human Early Life Exposome study, a study spanning six ongoing population-based European birth cohorts. Expression profiles of miRNAs in whole blood samples were quantified by microarray and tested for association with ADHD-related measures of behavior and neuropsychological functions from questionnaires (Conner's Rating Scale and Child Behavior Checklist) and computer-based tests (the N-back task and Attention Network Test).

RESULTS

We identified 29 miRNAs significantly associated (false discovery rate < .05) with the Conner's questionnaire-rated trait hyperactivity, 15 of which have been linked to ADHD in previous studies. Investigation into their biological relevance revealed involvement in several pathways related to neurodevelopment and function, as well as being linked with other neurodevelopmental or psychiatric disorders known to overlap with ADHD both in symptomology, genetic risk, and co-occurrence, such as autism spectrum disorder or schizophrenia. An additional three miRNAs were significantly associated with Conner's-rated inattention. No associations were found with questionnaire-rated total ADHD index or with computer-based tests.

CONCLUSIONS

The large overlap of our hyperactivity-associated miRNAs with previous studies on ADHD is intriguing and warrant further investigation. Though this study should be considered explorative and preliminary, these findings contribute towards identifying a set of miRNAs for use as blood-based biomarkers to aid in earlier and easier ADHD diagnosis.

Dysregulation of miR-1-3p: An Early Event in Colitis-Associated Dysplasia

Detection of colorectal dysplasia during surveillance colonoscopy remains the best method of determining risk for colitis-associated colorectal cancer (CAC). miRNAs (miRs) show great promise as tissue-specific biomarkers of neoplasia. The goal of this study was to explore the miR expression profile of precancerous dysplastic lesions in the AOM/DSS mouse model and identify early molecular changes associated with CAC. Epithelial cells were laser-microdissected from the colonic mucosa (inflamed versus dysplastic) of mice with AOM/DSS-induced colitis. A miR signature that can distinguish inflamed non-neoplastic mucosa from dysplasia was identified. Bioinformatic analyses led to the discovery of associated miR gene targets and enriched pathways and supported the construction of a network interaction map. miR-1a-3p was one of the miRs with the highest number of predicted targets, including Cdk6. Interestingly, miR-1a-3p and Cdk6 were down- and up-regulated in dysplastic lesions, respectively. Transfection of HCT116 and RKO cells with miR-1a-3p mimics induced apoptosis and cell cycle arrest in G1, suggesting its biological function. A slight reduction in the level of CDK6 transcripts was also observed in cells transfected with miR-1. These data provide novel insight into the early molecular alterations that accompany the development of CAC and identify a miR signature that represents a promising biomarker for the early detection of colitis-associated dysplasia.

Curcumin alters distinct molecular pathways in breast cancer subtypes revealed by integrated miRNA/mRNA expression analysis

BACKGROUND

Curcumin is well known for its anticancer properties. Its cytotoxic activity has been documented in several cancer cell lines, including breast cancer. The pleiotropic activity of curcumin as an antioxidant, an antiangiogenic, antiproliferative, and pro-apoptotic, is due to its diverse targets, such as signaling pathways, protein/enzyme, or noncoding gene.

AIM

This study aimed to identify key miRNAs and mRNAs induced by curcumin in breast cancer cells MCF7, T47D (hormone positive), versus MDA-MB231 (hormone negative) using comparative analysis of global gene expression profiles.

METHODS

RNA was isolated and subjected to mRNA and miRNA library sequencing to study the global gene expression profile of curcumin-treated breast cancer cells. The differential expression of gene and miRNA was performed using the DESeq R package. The enriched pathways were studied using cluster profileR, and integrated miRNA-mRNA analysis was carried out using miRtarvis and miRmapper tools.

RESULTS

Curcumin treatment led to upregulation of 59% TSGs in MCF7, 21% in MDA-MB-231 cells, and 36% TSGs in T47D, and downregulation of 57% oncogenes in MCF7, 76% in MDA-MB-231, and 91% in T47D. Similarly, curcumin treatment led to upregulation of 32% TSmiRs in MCF7, 37.5% in MDA-MB231, and 62.5% in T47D, and downregulation of 77% oncomiRs in MCF7, 50% in MDA-MB231 and 28.6% in T47D. Integrated analysis of miRNA-mRNA led to the identification of a common NFKB pathway altered by curcumin in all three cell lines. Analysis of uniquely enriched pathway revealed non-integrin membrane-ECM interactions and laminin interactions in MCF7; extracellular matrix organization and degradation in MDA-MB-231 and cell cycle arrest and G2/M transition in T47D.

CONCLUSION

Curcumin regulates miRNA and mRNA in a cell type-specific manner. The integrative analysis led to the detection of miRNAs and mRNAs pairs, which can be used as biomarkers associated with carcinogenesis, diagnostic, and treatment response in breast cancer.

Feline microRNAome in ovary and testis: Exploration of in-silico miRNA-mRNA networks involved in gonadal function and cellular stress response

The aim of the study was to perform the first in-depth analysis of miRNAs in ovarian and testicular tissues of the domestic cat, a critical biomedical model. Specifically, potential miRNA involvement was explored in gonadal function, testis development, and cellular stress response to preservation protocols. We performed miRNA-sequencing on 20 ovarian and 20 testicular samples from 15 cats, including different ages and tissue treatments. Using fresh tissues (n = 15), we confirmed gonadal expression of 183 miRNA precursors and discovered additional 52 novel feline candidate precursors. We integrated the mRNA data from our previous study on the same age and treatment groups to create in-silico miRNA-mRNA networks and their functional enrichment, which allows comprehensive exploration into possible miRNA functions in cat gonads. Clusters of miRNAs united by shared differentially expressed mRNA targets are potentially involved in testicular development and spermatogenesis. MicroRNAs could play a significant role in ovarian tissue response to stress from microwave-assisted dehydration, with smaller roles in cellular response to vitrification in both ovary and testis. This new list of miRNAs with potential function in cat gonads is a major step towards understanding the gonadal biology, as well as optimizing fertility preservation protocols.

ST08 Altered NF-κB Pathway in Breast Cancer Cells In Vitro as Revealed by miRNA-mRNA Analysis and Enhanced the Effect of Cisplatin on Tumour Reduction in EAC Mouse Model

ST08 is a novel curcumin derivative that exhibited apoptotic and anti-migratory activity in MDA-MB-231, triple-negative breast cancer cells reported earlier. In this study, we further explored the anticancer properties of ST08. ST08 reduced tumor burden in vivo and induced apoptosis through the mitochondrial pathway both in vitro and in vivo. ST08 potentiated the effect of cisplatin in vitro and in vivo in mouse EAC breast cancer models with minimal toxicity. ST08 induced alterations in the gene expression were studied by parallel analysis of miRNA and mRNA. 74 differentially expressed miRNA regulated 114 mRNA in triple-negative (MDA-MB-231) cancer cells. Pathway related to the ECM was altered in mesenchymal MDA-MB-231 cells. We constructed a unique miRNA-mRNA interaction network, and one of the pathways regulated by miRNA was NF-κB. Targets of NF-κB like MMP1, PTX3, and MMP2 were downregulated in MDA-MB-231 in response to ST08 treatment. PMA induced cell proliferation was abrogated by ST08 treatment, and no additional cell cytotoxicity was observed when used in combination with IKK-16 indicating ST08 regulation of NF-κB pathway in MDA-MB-231 cells.

Early downregulation of hsa-miR-144-3p in serum from drug-naïve Parkinson's disease patients

Advanced age represents one of the major risk factors for Parkinson's Disease. Recent biomedical studies posit a role for microRNAs, also known to be remodelled during ageing. However, the relationship between microRNA remodelling and ageing in Parkinson's Disease, has not been fully elucidated. Therefore, the aim of the present study is to unravel the relevance of microRNAs as biomarkers of Parkinson's Disease within the ageing framework. We employed Next Generation Sequencing to profile serum microRNAs from samples informative for Parkinson's Disease (recently diagnosed, drug-naïve) and healthy ageing (centenarians) plus healthy controls, age-matched with Parkinson's Disease patients. Potential microRNA candidates markers, emerging from the combination of differential expression and network analyses, were further validated in an independent cohort including both drug-naïve and advanced Parkinson's Disease patients, and healthy siblings of Parkinson's Disease patients at higher genetic risk for developing the disease. While we did not find evidences of microRNAs co-regulated in Parkinson's Disease and ageing, we report that hsa-miR-144-3p is consistently down-regulated in early Parkinson's Disease patients. Moreover, interestingly, functional analysis revealed that hsa-miR-144-3p is involved in the regulation of coagulation, a process known to be altered in Parkinson's Disease. Our results consistently show the down-regulation of hsa-mir144-3p in early Parkinson's Disease, robustly confirmed across a variety of analytical and experimental analyses. These promising results ask for further research to unveil the functional details of the involvement of hsa-mir144-3p in Parkinson's Disease.

miRNA and lncRNA Expression Networks Modulate Cell Cycle and DNA Repair Inhibition in Senescent Prostate Cells

Cellular senescence is a state of permanent growth arrest that arises once cells reach the limit of their proliferative capacity. It creates an inflammatory microenvironment favouring the initiation and progression of various age-related diseases, including prostate cancer. Non-coding RNAs (ncRNAs) have emerged as important regulators of cellular gene expression. Nonetheless, very little is known about the interplay of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) and how deregulation of ncRNA networks promotes cellular senescence. To investigate this, human prostate epithelial cells were cultured through different passages until senescent, and their RNA was extracted and sequenced using RNA sequencing (RNAseq) and microRNA sequencing (miRNA-seq) miRNAseq. Differential expression (DE) gene analysis was performed to compare senescent and proliferating cells with Limma, miRNA-target interactions with multiMiR, lncRNA-target interactions using TCGA data and network evaluation with miRmapper. We found that miR-335-3p, miR-543 and the lncRNAs H19 and SMIM10L2A all play central roles in the regulation of cell cycle and DNA repair processes. Expression of most genes belonging to these pathways were down-regulated by senescence. Using the concept of network centrality, we determined the top 10 miRNAs and lncRNAs, with miR-335-3p and H19 identified as the biggest hubs for miRNAs and lncRNA respectively. These ncRNAs regulate key genes belonging to pathways involved in cell senescence and prostate cancer demonstrating their central role in these processes and opening the possibility for their use as biomarkers or therapeutic targets to mitigate against prostate ageing and carcinogenesis.

MicroRNA Targeting

MicroRNAs (miRNAs) are small noncoding elements that play essential roles in the posttranscriptional regulation of biochemical processes. miRNAs recognize and target multiple mRNAs; therefore, investigating miRNA dysregulation is an indispensable strategy to understand pathological conditions and to design innovative drugs. Targeting miRNAs in diseases improve outcomes of several therapeutic strategies thus, this present study highlights miRNA targeting methods through experimental assays and bioinformatics tools. The first part of this review focuses on experimental miRNA targeting approaches for elucidating key biochemical pathways. A growing body of evidence about the miRNA world reveals the fact that it is not possible to uncover these molecules' structural and functional characteristics related to the biological processes with a deterministic approach. Instead, a systemic point of view is needed to truly understand the facts behind the natural complexity of interactions and regulations that miRNA regulations present. This task heavily depends both on computational and experimental capabilities. Fortunately, several miRNA bioinformatics tools catering to nonexperts are available as complementary wet-lab approaches. For this purpose, this work provides recent research and information about computational tools for miRNA targeting research.

MicroRNA Databases and Tools

In this era of big data, sets of methodologies and strategies are designed to extract knowledge from huge volumes of data. However, the cost of where and how to get this information accurately and quickly is extremely important, given the diversity of genomes and the different ways of representing that information. Among the huge set of information and relationships that the genome carries, there are sequences called miRNAs (microRNAs). These sequences were described in the 1990s and are mainly involved in mechanisms of regulation and gene expression. Having this in mind, this chapter focuses on exploring the available literature and providing useful and practical guidance on the miRNA database and tools topic. For that, we organized and present this text in two ways: (a) the update reviews and articles, which best summarize and discuss the theme; and (b) our update investigation on miRNA literature and portals about databases and tools. Finally, we present the main challenge and a possible solution to improve resources and tools.

Curcumin derivative ST09 modulates the miR-199a-5p/DDR1 axis and regulates proliferation and migration in ovarian cancer cells

Ovarian cancers are among the fatal malignancies affecting women globally, mainly due to their metastatic and chemoresistant nature. In this study, we report a potent curcumin derivative ST09 effective against ovarian cancers. Prior in-vitro studies with ST09 drug showed cytotoxicity in tumorigenic cells compared to normal cells and in-vivo, significant tumor reduction was observed with least systemic toxicity. ST09 induced cytotoxicity in the ovarian cancer cells triggering mitochondria-mediated intrinsic apoptotic pathway. Delving deeper to understand the underlying molecular mechanisms involved in ovarian cancer pathogenesis, we identified an inverse correlation of miR-199a-5p with DDR1, a collagen receptor with receptor tyrosine kinase activity. The ST09 treatment in ovarian cancer cell lines resulted in the deregulation of the miR-199a-5p/DDR1 axis, conferring tumor-suppressive functions. We established DDR1 to be a direct target of miR-199a-5p and that ST09-induced DDR1 loss in these ovarian cancer cells resulted in the inactivation of its downstream MMP activation, migration, EMT, and prosurvival NF-κB pathway. Overall this study demonstrates ST09, a potent drug candidate for ovarian cancer treatment which exhibits anti-invasive and migrastatic properties.

Long non-coding RNAs and their potential impact on diagnosis, prognosis, and therapy in prostate cancer: racial, ethnic, and geographical considerations

INTRODUCTION

Advances in high-throughput sequencing have greatly advanced our understanding of long non-coding RNAs (lncRNAs) in a relatively short period of time. This has expanded our knowledge of cancer, particularly how lncRNAs drive many important cancer phenotypes via their regulation of gene expression.

AREAS COVERED

Men of African descent are disproportionately affected by PC in terms of incidence, morbidity, and mortality. LncRNAs could serve as biomarkers to differentiate low-risk from high-risk diseases. Additionally, they may represent therapeutic targets for advanced and castrate-resistant cancer. We review current research surrounding lncRNAs and their association with PC. We discuss how lncRNAs can provide new insights and diagnostic biomarkers for African American men. Finally, we review advances in computational approaches that predict the regulatory effects of lncRNAs in cancer.

EXPERT OPINION

PC diagnostic biomarkers that offer high specificity and sensitivity are urgently needed. PC specific lncRNAs are compelling as diagnostic biomarkers owing to their high tissue and tumor specificity and presence in bodily fluids. Recent studies indicate that PCA3 clinical utility might be restricted to men of European descent. Further work is required to develop lncRNA biomarkers tailored for men of African descent.

Network based multifactorial modelling of miRNA-target interactions

Competing endogenous RNA (ceRNA) regulations and crosstalk between various types of non-coding RNA in humans is an important and under-explored subject. Several studies have pointed out that an alteration in miRNA:target interaction can result in unexpected changes due to indirect and complex interactions. In this article, we defined a new network-based model that incorporates miRNA:ceRNA interactions with expression values. Our approach calculates network-wide effects of perturbations in the expression level of one or more nodes in the presence or absence of miRNA interaction factors such as seed type, binding energy. We carried out the analysis of large-scale miRNA:target networks from breast cancer patients. Highly perturbing genes identified by our approach coincide with breast cancer-associated genes and miRNAs. Our network-based approach takes the sponge effect into account and helps to unveil the crosstalk between nodes in miRNA:target network. The model has potential to reveal unforeseen regulations that are only evident in the network context. Our tool is scalable and can be plugged in with emerging miRNA effectors such as circRNAs, lncRNAs, and available as R package ceRNAnetsim: https://www.bioconductor.org/packages/release/bioc/html/ceRNAnetsim.html.

Frequency of MicroRNA Response Elements Identifies Pathologically Relevant Signaling Pathways in Triple-Negative Breast Cancer

Complex interactions between mRNAs and microRNAs influence cellular functions. The mRNA-microRNA interactions also determine the post-transcriptional availability of mRNAs and unbound microRNAs. MicroRNAs binds to one or more microRNA response elements (MREs) located on the 3′UTR of mRNAs. In this study, we leveraged MREs and their frequencies in cancer and matched normal tissues to obtain insights into disease-specific interactions between mRNAs and microRNAs. We developed a bioinformatics method “ReMIx” that utilizes RNA sequencing (RNA-Seq) data to quantify MRE frequencies across the transcriptome. We applied ReMIx to triple-negative (TN) breast cancer tumor-normal adjacent pairs and identified MREs specific to TN tumors. ReMIx identified candidate mRNAs and microRNAs in the MAPK signaling cascade. Further analysis of MAPK gene regulatory networks revealed microRNA partners that influence and modulate MAPK signaling. In conclusion, we demonstrate a novel method of using MREs in the identification of functionally relevant mRNA-microRNA interactions in TN breast cancer.

•

Bioinformatics method ReMIx identify differential microRNA response rlements (MRE)

•

Tumor-specific MREs frequency observed in triple-negative breast cancer (TNBC)

•

MRE analysis identify MAPK signaling genes as therapeutic target for TNBC

•

MREs frequency can be used to identify pathologically relevant pathways

Computational Network Analysis Identifies Evolutionarily Conserved miRNA Gene Interactions Potentially Regulating Immune Response in Bovine Trypanosomosis

Bovine trypanosomosis is a devastating disease that causes huge economic loss to the global cattle industry on a yearly basis. Selection of accurate biomarkers are important in early disease diagnosis and treatment. Of late, micro-RNAs (miRNAs) are becoming the most useful biomarkers for both infectious and non-infectious diseases in humans, but this is not the case in animals. miRNAs are non-coding RNAs that regulate gene expression through binding to the 3'-, 5'-untranslated regions (UTR) or coding sequence (CDS) region of one or more target genes. The molecular identification of miRNAs that regulates the expression of immune genes responding to bovine trypanosomosis is poorly defined, as is the possibility that these miRNAs could serve as potential biomarkers for disease diagnosis and treatment currently unknown. To this end, we utilized in silico tools to elucidate conserved miRNAs regulating immune response genes during infection, in addition to cataloging significant genes. Based on the p value of 1.77E-32, we selected 25 significantly expressed immune genes. Using prediction analysis, we identified a total of 4,251 bovine miRNAs targeting these selected genes across the 3'UTR, 5'UTR and CDS regions. Thereafter, we identified candidate miRNAs based on the number of gene targets and their abundance at the three regions. In all, we found the top 13 miRNAs that are significantly conserved targeting 7 innate immune response genes, including bta-mir-2460, bta-mir-193a, bta-mir-2316, and bta-mir-2456. Our gene ontology analysis suggests that these miRNAs are involved in gene silencing, cellular protein modification process, RNA-induced silencing complex, regulation of humoral immune response mediated by circulating immunoglobulin and negative regulation of chronic inflammatory response, among others. In conclusion, this study identifies specific miRNAs that may be involved in the regulation of gene expression during bovine trypanosomosis. These miRNAs have the potential to be used as biomarkers in the animal and veterinary research community to facilitate the development of tools for early disease diagnosis/detection, drug targeting, and the rational design of drugs to facilitate disease treatment.

Clustering analysis of microRNA and mRNA expression data from TCGA using maximum edge-weighted matching algorithms

Background

microRNA (miRNA) is a short RNA (~ 22 nt) that regulates gene expression at the posttranscriptional level. Aberration of miRNA expressions could affect their targeting mRNAs involved in cancer-related signaling pathways. We conduct clustering analysis of miRNA and mRNA using expression data from the Cancer Genome Atlas (TCGA). We combine the Hungarian algorithm and blossom algorithm in graph theory. Data analysis is done using programming language R and Python.

Methods

We first quantify edge-weights of the miRNA-mRNA pairs by combining their expression correlation coefficient in tumor (T_CC) and correlation coefficient in normal (N_CC). We thereby introduce a bipartite graph partition procedure to identify cluster candidates. Specifically, we propose six weight formulas to quantify the change of miRNA-mRNA expression T_CC relative to N_CC, and apply the traditional hierarchical clustering to subjectively evaluate the different weight formulas of miRNA-mRNA pairs. Among these six different weight formulas, we choose the optimal one, which we define as the integrated mean value weights, to represent the connections between miRNA and mRNAs. Then the Hungarian algorithm and the blossom algorithm are employed on the miRNA-mRNA bipartite graph to passively determine the clusters. The combination of Hungarian and the blossom algorithms is dubbed maximum weighted merger method (MWMM).

Results

MWMM identifies clusters of different sizes that meet the mathematical criterion that internal connections inside a cluster are relatively denser than external connections outside the cluster and biological criterion that the intra-cluster Gene Ontology (GO) term similarities are larger than the inter-cluster GO term similarities. MWMM is developed using breast invasive carcinoma (BRCA) as training data set, but can also applies to other cancer type data sets. MWMM shows advantage in GO term similarity in most cancer types, when compared to other algorithms.

Conclusions

miRNAs and mRNAs that are likely to be affected by common underlying causal factors in cancer can be clustered by MWMM approach and potentially be used as candidate biomarkers for different cancer types and provide clues for targets of precision medicine in cancer treatment.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0562-z) contains supplementary material, which is available to authorized users.