Comparison of Hepatocellular Carcinoma miRNA Expression Profiling as Evaluated by Next Generation Sequencing and Microarray

MicroRNA (miRNA) expression profiling has proven useful in diagnosing and understanding the development and progression of several diseases. Microarray is the standard method for analyzing miRNA expression profiles; however, it has several disadvantages, including its limited detection of miRNAs. In recent years, advances in genome sequencing have led to the development of next-generation sequencing (NGS) technologies, which significantly advance genome sequencing speed and discovery. In this study, we compared the expression profiles obtained by next generation sequencing (NGS) with the profiles created using microarray to assess if NGS could produce a more accurate and complete miRNA profile. Total RNA from 14 hepatocellular carcinoma tumors (HCC) and 6 matched non-tumor control tissues were sequenced with Illumina MiSeq 50-bp single-end reads. Micro RNA expression profiles were estimated using miRDeep2 software. As a comparison, miRNA expression profiles for 11 out of 14 HCCs were also established by microarray (Agilent human microRNA microarray). The average total sequencing exceeded 2.2 million reads per sample and of those reads, approximately 57% mapped to the human genome. The average correlation for miRNA expression between microarray and NGS and subtraction were 0.613 and 0.587, respectively, while miRNA expression between technical replicates was 0.976. The diagnostic accuracy of HCC, p-value, and AUC were 90.0%, 7.22×10⁻⁴, and 0.92, respectively. In summary, NGS created an miRNA expression profile that was reproducible and comparable to that produced by microarray. Moreover, NGS discovered novel miRNAs that were otherwise undetectable by microarray. We believe that miRNA expression profiling by NGS can be a useful diagnostic tool applicable to multiple fields of medicine.     

miRNA 206 and miRNA 574-5p are highly expression in coronary artery disease

Coronary artery disease (CAD) is the leading cause of human morbidity and mortality worldwide. Innovative diagnostic biomarkers are a pressing need for this disease. miRNAs profiling is an innovative method of identifying biomarkers for many diseases and could be proven as a powerful tool in the diagnosis and treatment of CAD. We performed miRNA microarray analysis from the plasma of three CAD patients and three healthy controls. Subsequently, we performed quantitative real-time PCR (qRT-PCR) analysis of miRNA expression in plasma of another 67 CAD patients and 67 healthy controls. We identified two miRNAs (miR-206 and miR-574-5p) that were significantly up-regulated in CAD patients as compared with healthy controls (P<0.05). The receiver operating characteristic (ROC) curves indicated these two miRNAs had great potential to provide sensitive and specific diagnostic value for CAD.     

Comparative microarray analysis of microRNA expression profiles in primary cutaneous malignant melanoma, cutaneous malignant melanoma metastases, and benign melanocytic nevi

Perturbations in microRNA (miRNA) expression profiles have been reported for cutaneous malignant melanoma (CMM) predominantly when examined in cell lines. Despite the rapidly growing number of newly discovered human miRNA sequences, the availability of up-to-date miRNA expression profiles for clinical samples of primary cutaneous malignant melanoma (PCMM), cutaneous malignant melanoma metastases (CMMM), and benign melanocytic nevi (BMN) is limited. Specimens excised from the center of tumors (lesional) from patients with PCMM (n=9), CMMM (n=4), or BMN (n=8) were obtained during surgery. An exploratory microarray analysis was performed by miRNA expression profiling based on Agilent platform screening for 1205 human miRNAs. The results from the microarray analysis were validated by TaqMan quantitative real-time polymerase chain reaction. In addition to several miRNAs previously known to be associated with CMM, 19 unidentified miRNA candidates were found to be dysregulated in CMM patient samples. Among the 19 novel miRNA candidates, the genes hsa-miR-22, hsa-miR-130b, hsa-miR-146b-5p, hsa-miR-223, hsa-miR-301a, hsa-miR-484, hsa-miR-663, hsa-miR-720, hsa-miR-1260, hsa-miR-1274a, hsa-miR-1274b, hsa-miR-3663-3p, hsa-miR-4281, and hsa-miR-4286 were upregulated, and the genes hsa-miR-24-1*, hsa-miR-26a, hsa-miR-4291, hsa-miR-4317, and hsa-miR-4324 were downregulated. The results of this study partially confirm previous CMM miRNA profiling studies identifying miRNAs that are dysregulated in CMM. However, we report several novel miRNA candidates in CMM tumors; these miRNA sequences require further validation and functional analysis to evaluate whether they play a role in the pathogenesis of CMM.     

The use of miRNA microarrays for the analysis of cancer samples with global miRNA decrease

Recent studies have established that mutations or deletions in microRNA (miRNA) processing enzymes resulting in a global decrease of miRNA expression are frequent across cancers and can be associated with a poorer prognosis. While very popular in miRNA profiling studies, it remains unclear whether miRNA microarrays are suited or not to accurately detecting global miRNA decreases seen in cancers. In this work, we analyzed the miRNA profiles of samples with global miRNA decreases using Affymetrix miRNA microarrays following the inducible genetic deletion of Dicer1. Surprisingly, up to a third of deregulated miRNAs identified upon Dicer1 depletion were found to be up-regulated following standard robust multichip average (RMA) background correction and quantile normalization, indicative of normalization bias. Our comparisons of five preprocess steps performed at the probe level demonstrated that the use of cyclic loess relying on non-miRNA small RNAs present on the Affymetrix platform significantly improved specificity and sensitivity of detection of decreased miRNAs. These findings were validated in samples from patients with prostate cancer, where conjugation of robust normal-exponential background correction with cyclic loess normalization and array weights correctly identified the greatest number of decreased miRNAs, and the lowest amount of false-positive up-regulated miRNAs. These findings highlight the importance of miRNA microarray normalization for the detection of miRNAs that are truly differentially expressed and suggest that the use of cyclic loess based on non-miRNA small RNAs can help to improve the sensitivity and specificity of miRNA profiling in cancer samples with global miRNA decrease.     

Circulating micro-RNAs as potential blood-based markers for early stage breast cancer detection

Introduction

MicroRNAs (miRNAs, miRs) are a class of small, non-coding RNA molecules with relevance as regulators of gene expression thereby affecting crucial processes in cancer development. MiRNAs offer great potential as biomarkers for cancer detection due to their remarkable stability in blood and their characteristic expression in many different diseases. We investigated whether microarray-based miRNA profiling on whole blood could discriminate between early stage breast cancer patients and healthy controls.

Methods

We performed microarray-based miRNA profiling on whole blood of 48 early stage breast cancer patients at diagnosis along with 57 healthy individuals as controls. This was followed by a real-time semi-quantitative Polymerase Chain Reaction (RT-qPCR) validation in a separate cohort of 24 early stage breast cancer patients from a breast cancer screening unit and 24 age matched controls using two differentially expressed miRNAs (miR-202, miR-718).

Results

Using the significance level of p<0.05, we found that 59 miRNAs were differentially expressed in whole blood of early stage breast cancer patients compared to healthy controls. 13 significantly up-regulated miRNAs and 46 significantly down-regulated miRNAs in our microarray panel of 1100 miRNAs and miRNA star sequences could be detected. A set of 240 miRNAs that was evaluated by radial basis function kernel support vector machines and 10-fold cross validation yielded a specificity of 78.8%, and a sensitivity of 92.5%, as well as an accuracy of 85.6%. Two miRNAs were validated by RT-qPCR in an independent cohort. The relative fold changes of the RT-qPCR validation were in line with the microarray data for both miRNAs, and statistically significant differences in miRNA-expression were found for miR-202.

Conclusions

MiRNA profiling in whole blood has potential as a novel method for early stage breast cancer detection, but there are still challenges that need to be addressed to establish these new biomarkers in clinical use.     

miR-1207-5p and miR-1266 suppress gastric cancer growth and invasion by targeting telomerase reverse transcriptase

hTERT is the catalytic subunit of the telomerase complex. Elevated expression of hTERT is associated with the expansion and metastasis of gastric tumor. In this study, we aimed to identify novel tumor suppressor miRNAs that restrain hTERT expression. We began our screen for hTERT-targeting miRNAs with a miRNA microarray. miRNA candidates were further filtered by bioinformatic analysis, general expression pattern in different cell lines, gain-of-function effects on hTERT protein and the potential of these effects to suppress hTERT 3′ untranslated region (3′UTR) luciferase activity. The clinical relevance of two miRNAs (miR-1207-5p and miR-1266) was evaluated by real-time RT-PCR. The effects of these miRNAs on cell growth, cell cycle and invasion of gastric cancer cells were measured with CCK-8, flow cytometry and transwell assays. Finally, the ability of these miRNAs to suppress the transplanted tumors was also investigated. Fourteen miRNAs were identified using a combination of bioinformatics and miRNA microarray analysis. Of these fourteen miRNAs, nine were expressed at significantly lower levels in hTERT-positive cell lines compared with hTERT-negative cell lines and five could downregulate hTERT protein expression. Only miR-1207-5p and miR-1266 interacted with the 3′ UTR of hTERT and the expression levels of these two miRNAs were significantly decreased in gastric cancer tissues. These two miRNAs also inhibited gastric tumor growth in vitro and in vivo. Altogether, miR-1207-5p and miR-1266 were determined to be hTERT suppressors in gastric cancer, and the delivery of these two miRNAs represents a novel therapeutic strategy for gastric cancer treatment.     

Highly sensitive and specific microRNA expression profiling using BeadArray technology

We have developed a highly sensitive, specific and reproducible method for microRNA (miRNA) expression profiling, using the BeadArray™ technology. This method incorporates an enzyme-assisted specificity step, a solid-phase primer extension to distinguish between members of miRNA families. In addition, a universal PCR is used to amplify all targets prior to array hybridization. Currently, assay probes are designed to simultaneously analyse 735 well-annotated human miRNAs. Using this method, highly reproducible miRNA expression profiles were generated with 100–200 ng total RNA input. Furthermore, very similar expression profiles were obtained with total RNA and enriched small RNA species (R² ≥ 0.97). The method has a 3.5–4 log (10⁵–10⁹ molecules) dynamic range and is able to detect 1.2- to 1.3-fold-differences between samples. Expression profiles generated by this method are highly comparable to those obtained with RT–PCR (R² = 0.85–0.90) and direct sequencing (R = 0.87–0.89). This method, in conjunction with the 96-sample array matrix should prove useful for high-throughput expression profiling of miRNAs in large numbers of tissue samples.     

Labeled microRNA pull-down assay system: an experimental approach for high-throughput identification of microRNA-target mRNAs

We developed a simple, direct and cost-effective approach to search for the most likely target genes of a known microRNA (miRNA) in vitro. We term this method ‘labeled miRNA pull-down (LAMP)’ assay system. Briefly, the pre-miRNA is labeled with digoxigenin (DIG), mixed with cell extracts and immunoprecipitated by anti-DIG antiserum. When the DIG-labeled miRNA and bound mRNA complex are obtained, the total cDNAs are then subcloned and sequenced, or RT–PCR-amplified, to search for the putative target genes of a known miRNA. After successfully identifying the known target genes of Caenorhabditis elegans miRNAs lin-4 and let-7 and zebrafish let-7, we applied LAMP to find the unknown target gene of zebrafish miR-1, which resulted in the identification of hand2. We then confirmed hand2 as a novel target gene of miR-1 by whole-mount in situ hybridization and luciferase reporter gene assay. We further validated this target gene by microarray analysis, and the results showed that hand2 is the top-scoring among 302 predicted putative target genes. We concluded that LAMP is an experimental approach for high-throughput identification of the target gene of known miRNAs from both C. elegans and zebrafish, yielding fewer false positive results than those produced by using only the bioinformatics approach.     

Extremely Low-Frequency Electromagnetic Fields Affect the miRNA-Mediated Regulation of Signaling Pathways in the GC-2 Cell Line

Extremely low-frequency electromagnetic fields (ELF-EMFs) can affect male reproductive function, but the underlying mechanism of this effect remains unknown. miRNA-mediated regulation has been implicated as an important epigenetic mechanism for regulatory pathways. Herein, we profiled miRNA expression in response to ELF-EMFs in vitro. Mouse spermatocyte-derived GC–2 cells were intermittently exposed to a 50 Hz ELF-EMF for 72 h (5 min on/10 min off) at magnetic field intensities of 1 mT, 2 mT and 3 mT. Cell viability was assessed using the CCK–8 assay. Apoptosis and the cell cycle were analyzed with flow cytometry. miRNA expression was profiled using Affymetrix Mouse Genechip miRNA 3.0 arrays. Our data showed that the growth, apoptosis or cell cycle arrest of GC–2 cells exposed to the 50 Hz ELF-EMF did not significantly change. However, we identified a total of 55 miRNAs whose expression significantly changed compared with the sham group, including 19 differentially expressed miRNAs (7 miRNAs were upregulated, and 12 were downregulated) in the 1 mT exposure group and 36 (9 miRNAs were upregulated, and 27 were downregulated) in the 3 mT exposure group. The changes in the expression of 15 selected miRNAs measured by real-time PCR were consistent with the microarray results. A network analysis was used to predict core miRNAs and target genes, including miR-30e-5p, miR-210-5p, miR-196b-5p, miR-504-3p, miR-669c-5p and miR-455-3p. We found that these miRNAs were differentially expressed in response to different magnetic field intensities of ELF-EMFs. GO term and KEGG pathway annotation based on the miRNA expression profiling results showed that miRNAs may regulate circadian rhythms, cytokine-cytokine receptor interactions and the p53 signaling pathway. These results suggested that miRNAs could serve as potential biomarkers, and the miRNA-mediated regulation of signaling pathways might play significant roles in the biological effects of ELF-EMFs.     

MicroRNA Expression Differences in Human Hematopoietic Cell Lineages Enable Regulated Transgene Expression

Blood microRNA (miRNA) levels have been associated with and shown to participate in disease pathophysiology. However, the hematopoietic cell of origin of blood miRNAs and the individual blood cell miRNA profiles are poorly understood. We report the miRNA content of highly purified normal hematopoietic cells from the same individuals. Although T-cells, B-cells and granulocytes had the highest miRNA content per cell, erythrocytes contributed more cellular miRNA to the blood, followed by granulocytes and platelets. miRNA profiling revealed different patterns and different expression levels of miRNA specific for each lineage. miR-30c-5p was determined to be an appropriate reference normalizer for cross-cell qRT-PCR comparisons. miRNA profiling of 5 hematopoietic cell lines revealed differential expression of miR-125a-5p. We demonstrated endogenous levels of miR-125a-5p regulate reporter gene expression in Meg-01 and Jurkat cells by (1) constructs containing binding sites for miR-125a-5p or (2) over-expressing or inhibiting miR-125a-5p. This quantitative analysis of the miRNA profiles of peripheral blood cells identifies the circulating hematopoietic cellular miRNAs, supports the use of miRNA profiles for distinguishing different hematopoietic lineages and suggests that endogenously expressed miRNAs can be exploited to regulate transgene expression in a cell-specific manner.     

A novel microRNA-based signature predicts prognosis among nasopharyngeal cancer patients

Nasopharyngeal cancer is one of the most common malignant tumors in the head and neck. Identification of promising miRNA biomarkers might benefit a lot to the detection of nasopharyngeal carcinoma. miRNA expression profile and clinical information were obtained from two microarray profiling data sets from the Gene Expression Omnibus (GEO) database. miRNA signature model was constructed via univariate Cox survival analysis, multivariate Cox survival analysis, and least absolute shrinkage and selection operator Cox regression analysis. Kaplan–Meier curve, area under the curve (AUC), decision curve analysis, Box plot, and nomogram were used to evaluate the prognosis of the model to patients. 67 up-regulated and 93 down-regulated miRNAs were identified from GEO microarray data sets (P < 0.05). A three-miRNA signature (has-miR-142-3p, has-miR-29c, and has-miR-30e) was obviously associated with the overall survival of nasopharyngeal carcinoma patients (P  < 0.001). The AUCs for the signature were 0.74, 0.7 for the training set and external validation set. The AUC of disease free survival and distant metastasis-free survival were also high. The model has better clinical independence and has better clinical prediction effect when combined with clinical characteristics (P < 0.0001). Compared with the published models, our model had a higher AUC. Our results revealed that a three-miRNA signature was a potential novel prognostic biomarker for nasopharyngeal carcinoma.Impact statementNasopharyngeal cancer is one of the most common malignant tumors in the head and neck. Identification of promising miRNA biomarkers might benefit a lot to the detection of nasopharyngeal carcinoma. A three-miRNA signature (has-miR-142-3p, has-miR-29c, and has-miR-30e) was obviously associated with the overall survival of nasopharyngeal carcinoma patients. The model has better clinical independence and has better clinical prediction effect when combined with clinical characteristics. Our results revealed that a three-miRNA signature was a potential novel prognostic biomarker for nasopharyngeal carcinoma.     

A Genome-Wide Investigation of MicroRNA Expression Identifies Biologically-Meaningful MicroRNAs That Distinguish between High-Risk and Low-Risk Intraductal Papillary Mucinous Neoplasms of the Pancreas

Background

Intraductal papillary mucinous neoplasms (IPMNs) are pancreatic ductal adenocarcinoma (PDAC) precursors. Differentiating between high-risk IPMNs that warrant surgical resection and low-risk IPMNs that can be monitored is a significant clinical problem, and we sought to discover a panel of mi(cro)RNAs that accurately classify IPMN risk status.

Methodology/Principal Findings

In a discovery phase, genome-wide miRNA expression profiling was performed on 28 surgically-resected, pathologically-confirmed IPMNs (19 high-risk, 9 low-risk) using Taqman MicroRNA Arrays. A validation phase was performed in 21 independent IPMNs (13 high-risk, 8 low-risk). We also explored associations between miRNA expression level and various clinical and pathological factors and examined genes and pathways regulated by the identified miRNAs by integrating data from bioinformatic analyses and microarray analysis of miRNA gene targets. Six miRNAs (miR-100, miR-99b, miR-99a, miR-342-3p, miR-126, miR-130a) were down-regulated in high-risk versus low-risk IPMNs and distinguished between groups (P<10⁻³, area underneath the curve (AUC) = 87%). The same trend was observed in the validation phase (AUC = 74%). Low miR-99b expression was associated with main pancreatic duct involvement (P = 0.021), and serum albumin levels were positively correlated with miR-99a (r = 0.52, P = 0.004) and miR-100 expression (r = 0.49, P = 0.008). Literature, validated miRNA:target gene interactions, and pathway enrichment analysis supported the candidate miRNAs as tumor suppressors and regulators of PDAC development. Microarray analysis revealed that oncogenic targets of miR-130a (ATG2B, MEOX2), miR-342-3p (DNMT1), and miR-126 (IRS-1) were up-regulated in high- versus low-risk IPMNs (P<0.10).

Conclusions

This pilot study highlights miRNAs that may aid in preoperative risk stratification of IPMNs and provides novel insights into miRNA-mediated progression to pancreatic malignancy. The miRNAs identified here and in other recent investigations warrant evaluation in biofluids in a well-powered prospective cohort of individuals newly-diagnosed with IPMNs and other pancreatic cysts and those at increased genetic risk for these lesions.     

Identification and Validation of miRNAs as Endogenous Controls for RQ-PCR in Blood Specimens for Breast Cancer Studies

Introduction

A prerequisite to accurate interpretation of RQ-PCR data is robust data normalization. A commonly used method is to compare the cycle threshold (C_T) of target miRNAs with those of a stably expressed endogenous (EC) miRNA(s) from the same sample. Despite the large number of studies reporting miRNA expression patterns, comparatively few appropriate ECs have been reported thus far. The purpose of this study was to identify stably expressed miRNAs with which to normalize RQ-PCR data derived from human blood specimens.

Methods

MiRNA profiling of approximately 380 miRNAs was performed on RNA derived from blood specimens from 10 women with breast cancer and 10 matched controls. Analysis of mean expression values across the dataset (GME) identified stably expressed candidates. Additional candidates were selected from the literature and analyzed by the geNorm algorithm. Further validation of three candidate ECs by RQ-PCR was performed in a larger cohort (n = 40 cancer, n = 20 control) was performed, including analysis by geNorm and NormFinder algorithms.

Results

Microarray screening identified 10 candidate ECs with expression patterns closest to the global mean. Geometric averaging of candidate ECs from the literature using geNorm identified miR-425 as the most stably expressed miRNA. MiR-425 and miR-16 were the best combination, achieving the lowest V-value of 0.185. Further validation by RQ-PCR confirmed that miR-16 and miR-425 were the most stably expressed ECs overall. Their combined use to normalize expression data enabled the detection of altered target miRNA expression that reliably differentiated between cancers and controls in human blood specimens.

Conclusion

This study identified that the combined use of 2 miRNAs, (miR-16 and miR-425) to normalize RQ-PCR data generated more reliable results than using either miRNA alone, or use of U6. Further investigation into suitable ECs for use in miRNA RQ-PCR studies is warranted.