High-Throughput Sequencing of Plasma MicroRNA in Chronic Fatigue Syndrome/Myalgic Encephalomyelitis

Background

MicroRNAs (miRNAs) are known to regulate many biological processes and their dysregulation has been associated with a variety of diseases including Chronic Fatigue Syndrome/Myalgic Encephalomyelitis (CFS/ME). The recent discovery of stable and reproducible miRNA in plasma has raised the possibility that circulating miRNAs may serve as novel diagnostic markers. The objective of this study was to determine the role of plasma miRNA in CFS/ME.

Results

Using Illumina high-throughput sequencing we identified 19 miRNAs that were differentially expressed in the plasma of CFS/ME patients in comparison to non-fatigued controls. Following RT-qPCR analysis, we were able to confirm the significant up-regulation of three miRNAs (hsa-miR-127-3p, hsa-miR-142-5p and hsa-miR-143-3p) in the CFS/ME patients.

Conclusion

Our study is the first to identify circulating miRNAs from CFS/ME patients and also to confirm three differentially expressed circulating miRNAs in CFS/ME patients, providing a basis for further study to find useful CFS/ME biomarkers.     

Site-Specific MicroRNA Expression May Lead to Different Subtypes in Ulcerative Colitis

Background and Aim

Ulcerative Colitis (UC) is a type of inflammatory bowel disease, considered as an important disease of gastrointestinal tract having a huge impact on the health of the patient. Prolonged inflammation of colon in UC patients increases the risk of developing colorectal cancer. MiRNA are reported as a connecting link between inflammation and cancer. Differential miRNA expression is reported in Crohn’s disease (CD) patients involving various regions of the gastrointestinal tract. The current study was performed to dissect out the site specific miRNA expression in the colon biopsy samples of UC patients from Northern India.

Methods

Biopsy samples were collected from UC patients and healthy controls from Rectosigmoid Area (RS) and Ascending Colon (AC). MiRNA expression was compared between patients with RS and AC using a microarray platform. Differential expression was further validated by Real Time PCR analysis. Demographic and pathological data of UC -associated CRC patients was collected from the hospital database and analyzed for assessing the site of cancer.

Results

Upon analysis of data generated on a microarray platform and qRT PCR revealed that the expression of six miRNAs hsa-miR-146b-5p, hsa-miR-335-3p, hsa-miR-342-3p, hsa-miR-644b-3p, hsa-miR-491-3p, hsa-miR-4732-3p were downregulated in patients where RS was involved as compared to AC. The expression of hsa-miR-141-3p was upregulated in patients where RS region was involved as compared to AC. Analysis of the registered UC patient’s database from the hospital revealed that the site of CRC was predomimnantly the rectosigmoid region of the colon in most of the cases.

Conclusion

This is the first study to show the differential expression of miRNA involving different sites of colon in UC patients. Taking our data and previous reports into consideration, we propose that differential miRNA expression during UC perhaps contribute in the development of UC-associated CRC at the rectosigmoid area.     

Differential MicroRNA Expression Profile between Stimulated PBMCs from HIV-1 Infected Elite Controllers and Viremic Progressors

Background

The emerging relationship between microRNAs (miRNA) and viral-control is a topic of interest in the field of HIV. Host-genome might play an important role in the control of viremia. The aim of this study was to assess the specific miRNA profile that could contribute to the control of HIV replication in Elite Controllers

Results

After adequate normalization, expression profile of 286 human miRNAs (hsa-miR) was evaluated in phytohaemagglutinin-stimulated PBMCs from 29 individuals classified in 4 groups: 8 elite controllers (EC; viral load <50 cp/ml without treatment), 8 viremic progressors (VP; VL>5000 cp/ml without treatment), 8 patients under antiretroviral treatment (ART; VL<200 cp/ml) and 5 uninfected individuals (HIV-) through TaqMan Array Human microRNA Cards v3.0. A differential expression pattern consisting of 23 miRNAs became significantly different when comparing EC and VP. Profiling analysis segregated the population in two different blocks: while EC and HIV- clustered together in the same block (EC/HIV-_block 1), VP and ART individuals clustered together in a second block (VP/ART_block 2). Two inversely expressed miRNA patterns were determined within those two blocks: a set of 4 miRNAs (hsa-miR-221, -27a, -27b and -29b) was up-expressed in EC/HIV-_block and down-expressed in VP/ART_block while 19 miRNAs were down-expressed in block 1 and up-expressed in block 2. Differential miRNAs were successfully validated through individual RT-qPCR assays.

Conclusions

Profile in EC resembled HIV- and differentially clusters with VP and ART. Therefore, differential clustering does not rely on undetectable viremia.     

Salivary MicroRNA in Pancreatic Cancer Patients

Background

Pancreatic cancer is the fourth leading cause of cancer death in Western countries, with the lowest 1-year survival rate among commonly diagnosed cancers. Reliable biomarkers for pancreatic cancer diagnosis are lacking and are urgently needed to allow for curative surgery. As microRNA (miRNA) recently emerged as candidate biomarkers for this disease, we explored in the present pilot study the differences in salivary microRNA profiles between patients with pancreatic tumors that are not eligible for surgery, precancerous lesions, inflammatory disease or cancer-free patients as a potential early diagnostic tool.

Methods

Whole saliva samples from patients with pancreatic cancer (n = 7), pancreatitis (n = 4), IPMN (n = 2), or healthy controls (n = 4) were obtained during endoscopic examination. After total RNA isolation, expression of 94 candidate miRNAs was screened by q(RT)PCR using Biomark Fluidgm. Human-derived pancreatic cancer cells were xenografted in athymic mice as an experimental model of pancreatic cancer.

Results

We identified hsa-miR-21, hsa-miR-23a, hsa-miR-23b and miR-29c as being significantly upregulated in saliva of pancreatic cancer patients compared to control, showing sensitivities of 71.4%, 85.7%, 85,7% and 57%, respectively and excellent specificity (100%). Interestingly, hsa-miR-23a and hsa-miR23b are overexpressed in the saliva of patients with pancreatic cancer precursor lesions. We found that hsa-miR-210 and let-7c are overexpressed in the saliva of patients with pancreatitis as compared to the control group, with sensitivity of 100% and 75%, and specificity of 100% and 80%, respectively. Last hsa-miR-216 was upregulated in cancer patients as compared to patients diagnosed with pancreatitis, with sensitivity of 50% and specificity of 100%. In experimental models of PDAC, salivary microRNA detection precedes systemic detection of cancer cells markers.

Conclusions

Our novel findings indicate that salivary miRNA are discriminatory in pancreatic cancer patients that are not eligible for surgery. In addition, we demonstrate in experimental models that salivary miRNA detection precedes systemic detection of cancer cells markers. This study stems for the use of salivary miRNA as biomarker for the early diagnosis of patients with unresectable pancreatic cancer.     

MiRNA Analysis by Quantitative PCR in Preterm Human Breast Milk Reveals Daily Fluctuations of hsa-miR-16-5p

Background and Aims

Human breast milk is an extremely dynamic fluid containing many biologically-active components which change throughout the feeding period and throughout the day. We designed a miRNA assay on minimized amounts of raw milk obtained from mothers of preterm infants. We investigated changes in miRNA expression within month 2 of lactation and then over the course of 24 hours.

Materials and Methods

Analyses were performed on pooled breast milk, made by combining samples collected at different clock times from the same mother donor, along with time series collected over 24 hours from four unsynchronized mothers. Whole milk, lipids or skim milk fractions were processed and analyzed by qPCR. We measured hsa-miR-16-5p, hsa-miR-21-5p, hsa-miR-146-5p, and hsa-let-7a, d and g (all -5p). Stability of miRNA endogenous controls was evaluated using RefFinder, a web tool integrating geNorm, Normfinder, BestKeeper and the comparative ΔΔCt method.

Results

MiR-21 and miR-16 were stably expressed in whole milk collected within month 2 of lactation from four mothers. Analysis of lipids and skim milk revealed that miR-146b and let-7d were better references in both fractions. Time series (5H-23H) allowed the identification of a set of three endogenous reference genes (hsa-let-7d, hsa-let-7g and miR-146b) to normalize raw quantification cycle (Cq) data. We identified a daily oscillation of miR-16-5p.

Perspectives

Our assay allows exploring miRNA levels of breast milk from mother with preterm baby collected in time series over 48–72 hours.     

MicroRNA Signature Characterizes Primary Tumors That Metastasize in an Esophageal Adenocarcinoma Rat Model

Objective

To establish a miRNA signature for metastasis in an animal model of esophageal adenocarcinoma (EAC).

Background

The incidence of esophageal adenocarcinoma (EAC) has dramatically increased and esophageal cancer is now the sixth leading cause of cancer deaths worldwide. Mortality rates remain high among patients with advanced stage disease and esophagectomy is associated with high complication rates. Hence, early identification of potentially metastatic disease would better guide treatment strategies.

Methods

The modified Levrat’s surgery was performed to induce EAC in Sprague-Dawley rats. Primary EAC and distant metastatic sites were confirmed via histology and immunofluorescence. miRNA profiling was performed on primary tumors with or without metastasis. A unique subset of miRNAs expressed in primary tumors and metastases was identified with Ingenuity Pathway Analysis (IPA) along with upstream and downstream targets. miRNA-linked gene expression analysis was performed on a secondary cohort of metastasis positive (n=5) and metastasis negative (n=28) primary tumors.

Results

The epithelial origin of distant metastasis was established by IF using villin (VIL1) and mucin 5AC (MUC5AC) antibodies. miRNome analysis identified four down-regulated miRNAs in metastasis positive primary tumors compared to metastasis negative tumors: miR-92a-3p (p=0.0001), miR-141-3p (p=0.0022), miR-451-1a (p=0.0181) and miR133a-3p (p=0.0304). Six target genes identified in the top scoring networks by IPA were validated as significantly, differentially expressed in metastasis positive primary tumors: Ago2, Akt1, Kras, Bcl2L11, CDKN1B and Zeb2.

Conclusion

In vivo metastasis was confirmed in the modified Levrat’s model. Analysis of the primary tumor identified a distinctive miRNA signature for primary tumors that metastasized.     

Integrated miRNA and mRNA Expression Profiling in Inflamed Colon of Patients with Ulcerative Colitis

Background

Ulcerative colitis (UC) is associated with differential colonic expression of genes involved in immune response (e.g. IL8) and barrier integrity (e.g. cadherins). MicroRNAs (miRNAs) are regulators of gene expression and are involved in various immune-related diseases. In this study, we investigated (1) if miRNA expression in UC mucosa is altered and (2) if any of these changes correlate with mucosal mRNA expression. Integration of mRNA and miRNA expression profiling may allow the identification of functional links between dysregulated miRNAs and their target mRNA.

Methodology

Colonic mucosal biopsies were obtained from 17 UC (10 active and 7 inactive) patients and 10 normal controls. Total RNA was used to analyze miRNA and mRNA expression via Affymetrix miRNA 2.0 and Affymetrix Human Gene 1.0ST arrays, respectively. Both miRNA and gene expression profiles were integrated by correlation analysis to identify dysregulated miRNAs with their corresponding predicted target mRNA. Microarray data were validated with qRT-PCR. Regulation of IL8 and CDH11 expression by hsa-miR-200c-3p was determined by luciferase reporter assays.

Results

When comparing active UC patients vs. controls, 51 miRNAs and 1543 gene probe sets gave significantly different signals. In contrast, in inactive UC vs. controls, no significant miRNA expression differences were found while 155 gene probe sets had significantly different signals. We then identified potential target genes of the significantly dysregulated miRNAs and genes in active UC vs. controls and found a highly significant inverse correlation between hsa-miR-200c-3p and IL8, an inflammatory marker, and between hsa-miR-200c-3p and CDH11, a gene related to intestinal epithelial barrier function. We could demonstrate that hsa-miR-200c-3p directly regulates IL8 and CDH11 expression.

Conclusion

Differential expression of immune- and barrier-related genes in inflamed UC mucosa may be influenced by altered expression of miRNAs. Integrated analysis of miRNA and mRNA expression profiles revealed hsa-miR-200c-3p for use of miRNA mimics as therapeutics.     

Towards Clinical Applications of Blood-Borne miRNA Signatures: The Influence of the Anticoagulant EDTA on miRNA Abundance

Background

Circulating microRNAs (miRNAs) from blood are increasingly recognized as biomarker candidates for human diseases. Clinical routine settings frequently include blood sampling in tubes with EDTA as anticoagulant without considering the influence of phlebotomy on the overall miRNA expression pattern. We collected blood samples from six healthy individuals each in an EDTA blood collection tube. Subsequently, the blood was transferred into PAXgene^TM tubes at three different time points, i.e. directly (0 min), 10 min, and 2 h after phlebotomy. As control blood was also directly collected in PAXgene^TM blood RNA tubes that contain a reagent to directly lyse blood cells and stabilize their content. For all six blood donors at the four conditions (24 samples) we analyzed the abundance of 1,205 miRNAs by human Agilent miRNA V16 microarrays.

Results

While we found generally a homogenous pattern of the miRNA abundance in all 24 samples, the duration of the EDTA treatment appears to influence the miRNA abundance of specific miRNAs. The most significant changes are observed after longer EDTA exposition. Overall, the impact of the different blood sample conditions on the miRNA pattern was substantially lower than intra-individual variations. While samples belonging to one of the six individuals mostly cluster together, there was no comparable clustering for any of the four tested blood sampling conditions. The most affected miRNA was miR-769-3p that was not detected in any of the six PAXgene blood samples, but in all EDTA 2h samples. Accordingly, hsa-miR-769-3p was also the only miRNA that showed a significantly different abundance between the 4 blood sample conditions by an ANOVA analysis (Benjamini-Hochberg adjusted p-value of 0.003). Validation by qRT-PCR confirmed this finding.

Conclusion

The pattern of blood-borne miRNA abundance is rather homogenous between the four tested blood sample conditions of six blood donors. There was a clustering between the miRNA profiles that belong to a specific blood donor, but not between any of the four tested blood sampling conditions. The results show a limited overall impact of the blood sampling conditions on the miRNA pattern. Notwithstanding, the abundance of single miRNAs can be significantly altered by different blood sampling conditions.     

Reproducibility of Circulating MicroRNAs in Stored Plasma Samples

Background

Most studies of microRNA (miRNA) and disease have examined tissue-specific expression in limited numbers of samples. The presence of circulating miRNAs in plasma samples provides the opportunity to examine prospective associations between miRNA expression and disease in initially healthy individuals. However, little data exist on the reproducibility of miRNAs in stored plasma.

Methods

We used Real-Time PCR to measure 61 pre-selected microRNA candidates in stored plasma. Coefficients of variation (CVs) were used to assess inter-assay reliability (n = 15) and within-person stability over one year (n = 80). Intraclass correlation coefficients (ICCs) and polychoric correlation coefficients were used to assess within-person stability and delayed processing reproducibility (whole blood stored at 4°C for 0, 24 and 48 hours; n = 12 samples).

Results

Of 61 selected miRNAs, 23 were detected in at least 50% of samples and had average CVs below 20% for inter-assay reproducibility and 31 for delayed processing reproducibility. Ten miRNAs were detected in at least 50% of samples, had average CVs below 20% and had ICCs above 0.4 for within-person stability over 1–2 years, six of which satisfied criteria for both interassay reproducibility and short-term within-person stability (miR-17-5p, -191-5p, -26a-5p, -27b-3p, -320a, and -375) and two all three types of reproducibility (miR-27b-3p and -26a-5p). However, many miRNAs with acceptable average CVs had high maximum CVs, most had low expression levels, and several had low ICCs with delayed processing.

Conclusions

About a tenth of miRNAs plausibly related to chronic disease were reliably detected in stored samples of healthy adults.     

MicroRNA expression signatures of bladder cancer revealed by deep sequencing

Background

MicroRNAs (miRNAs) are a class of small noncoding RNAs that regulate gene expression. They are aberrantly expressed in many types of cancers. In this study, we determined the genome-wide miRNA profiles in bladder urothelial carcinoma by deep sequencing.

Methodology/Principal Findings

We detected 656 differentially expressed known human miRNAs and miRNA antisense sequences (miRNA*s) in nine bladder urothelial carcinoma patients by deep sequencing. Many miRNAs and miRNA*s were significantly upregulated or downregulated in bladder urothelial carcinoma compared to matched histologically normal urothelium. hsa-miR-96 was the most significantly upregulated miRNA and hsa-miR-490-5p was the most significantly downregulated one. Upregulated miRNAs were more common than downregulated ones. The hsa-miR-183, hsa-miR-200b∼429, hsa-miR-200c∼141 and hsa-miR-17∼92 clusters were significantly upregulated. The hsa-miR-143∼145 cluster was significantly downregulated. hsa-miR-182, hsa-miR-183, hsa-miR-200a, hsa-miR-143 and hsa-miR-195 were evaluated by Real-Time qPCR in a total of fifty-one bladder urothelial carcinoma patients. They were aberrantly expressed in bladder urothelial carcinoma compared to matched histologically normal urothelium (p<0.001 for each miRNA).

Conclusions/Significance

To date, this is the first study to determine genome-wide miRNA expression patterns in human bladder urothelial carcinoma by deep sequencing. We found that a collection of miRNAs were aberrantly expressed in bladder urothelial carcinoma compared to matched histologically normal urothelium, suggesting that they might play roles as oncogenes or tumor suppressors in the development and/or progression of this cancer. Our data provide novel insights into cancer biology.     

Reorganization of metastamiRs in the evolution of metastatic aggressive neuroblastoma cells

Background

MetastamiRs have momentous clinical relevance and have been correlated with disease progression in many tumors. In this study, we identified neuroblastoma metastamiRs exploiting unique mouse models of favorable and high-risk metastatic human neuroblastoma. Further, we related their deregulation to the modulation of target proteins and established their association with clinical outcomes.

Results

Whole genome miRNA microarray analysis identified 74 metastamiRs across the manifold of metastatic tumors. RT-qPCR on select miRNAs validated profile expression. Results from bio-informatics across the ingenuity pathway, miRCancer, and literature data-mining endorsed the expression of these miRNAs in multiple tumor systems and showed their role in metastasis, identifying them as metastamiRs. Immunoblotting and TMA-IHC analyses revealed alterations in the expression/phosphorylation of metastamiRs’ targets, including ADAMTS-1, AKT1/2/3, ASK1, AURKβ, Birc1, Birc2, Bric5, β-CATENIN, CASP8, CD54, CDK4, CREB, CTGF, CXCR4, CYCLIN-D1, EGFR, ELK1, ESR1, CFOS, FOSB, FRA, GRB10, GSK3β, IL1α, JUND, kRAS, KRTAP1, MCP1, MEGF10, MMP2, MMP3, MMP9, MMP10, MTA2, MYB, cMYC, NF2, NOS3, P21, pP38, PTPN3, CLEAVED PARP, PKC, SDF-1β, SEMA3D, SELE, STAT3, TLR3, TNFα, TNFR1, and VEGF in aggressive cells ex vivo and in a manifold of metastatic tumors in vivo. miRNA mimic (hsa-miR-125b, hsa-miR-27b, hsa-miR-93, hsa-miR-20a) and inhibitor (hsa-miR-1224-3p, hsa-miR-1260) approach for select miRNAs revealed the direct influence of the altered metastamiRs in the regulation of identified protein targets. Clinical outcome association analysis with the validated metastamiRs’ targets corresponded strongly with poor overall and relapse-free survival.

Conclusions

For the first time, these results identified a comprehensive list of neuroblastoma metastamiRs, related their deregulation to altered expression of protein targets, and established their association with poor clinical outcomes. The identified set of distinctive neuroblastoma metastamiRs could serve as potential candidates for diagnostic markers for the switch from favorable to high-risk metastatic disease.

Electronic supplementary material

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