miRSeq: a user-friendly standalone toolkit for sequencing quality evaluation and miRNA profiling

Metformin inhibits gastric cancer cell proliferation by regulation of a novel Loc100506691-CHAC1 axis

Long noncoding RNAs (lncRNAs) are a group of nonprotein coding transcripts that play a critical role in cancer progression. However, the role of lncRNA in metformin-induced inhibition of cell growth and its biological function in gastric cancer remain largely unknown. In this study, we identified an oncogenic lncRNA, Loc100506691, the expression of which was decreased in gastric cancer cells with metformin treatment. Moreover, Loc100506691 was significantly overexpressed in gastric cancer compared with adjacent normal tissues (p < 0.001), and high Loc100506691 expression was significantly correlated with poor survival of patients with gastric cancer. Additionally, Loc100506691 knockdown could significantly suppress gastric cancer cell growth in vitro, and ectopic Loc100506691 expression accelerated tumor growth in an in vivo mouse model. Analysis of the cell cycle revealed that Loc100506691 knockdown induced cell cycle arrest at the G2/M phase by impairing cell entry from the G2/M to G1 phase. Loc100506691 negatively regulated CHAC1 expression by modulating miR-26a-5p/miR-330-5p expression, and CHAC1 knockdown markedly attenuated Loc100506691 knockdown-induced gastric cancer cell growth and motility suppression. We concluded that anti-proliferative effects of metformin in gastric cancer may be partially caused by suppression of the Loc100506691-miR-26a-5p/miR-330-5p-CHAC1 axis.

LOC550643, a Long Non-coding RNA, Acts as Novel Oncogene in Regulating Breast Cancer Growth and Metastasis

Metastatic disease is responsible for over 90% of death in patients with breast cancer. Therefore, identifying the molecular mechanisms that regulate metastasis and developing useful therapies are crucial tasks. Long non-coding RNAs (lncRNAs), which are non-coding transcripts with >200 nucleotides, have recently been identified as critical molecules for monitoring cancer progression. This study examined the novel lncRNAs involved in the regulation of tumor progression in breast cancer. This study identified 73 metastasis-related lncRNA candidates from comparison of paired isogenic high and low human metastatic breast cancer cell lines, and their expression levels were verified in clinical tumor samples by using The Cancer Genome Atlas. Among the cell lines, a novel lncRNA, LOC550643, was highly expressed in breast cancer cells. Furthermore, the high expression of LOC550643 was significantly correlated with the poor prognosis of breast cancer patients, especially those with triple-negative breast cancer. Knockdown of LOC550643 inhibited cell proliferation of breast cancer cells by blocking cell cycle progression at S phase. LOC550643 promoted important in vitro metastatic traits such as cell migration and invasion. Furthermore, LOC550643 could inhibit miR-125b-2-3p expression to promote breast cancer cell growth and invasiveness. In addition, by using a xenograft mouse model, we demonstrated that depletion of LOC550643 suppressed the lung metastatic potential of breast cancer cells. Overall, our study shows that LOC550643 plays an important role in breast cancer cell metastasis and growth, and LOC550643 could be a potential diagnosis biomarker and therapeutic target for breast cancer.

Identifying Circulating MicroRNA in Kawasaki Disease by Next-Generation Sequencing Approach

Kawasaki disease (KD) typically occurs in children aged under 5 years and can cause coronary artery lesions (CALs). Early diagnosis and treatment with intravenous immunoglobulin can reduce the occurrence of CALs; therefore, identifying a good biomarker for diagnosing KD is essential. Here, using next-generation sequencing in patients with recurrent KD, those with viral infection, and healthy controls, we identified dysregulated circulating microRNAs as diagnostic biomarkers for KD. Pathway enrichment analysis illustrated the putative role of these miRNAs in KD progression. Their expression levels were validated using real-time polymerase chain reaction (qPCR). Fifteen dysregulated circulating miRNAs (fold changes >2 and <0.5) were differentially expressed in the recurrent KD group compared with the viral infection and control groups. These miRNAs were significantly involved in the transforming growth factor-β, epithelial-mesenchymal transition, and cell apoptosis signaling pathways. Notably, their expression levels were frequently restored after intravenous immunoglobulin treatment. Among the candidates, miR-24-3p expression level was significantly higher in patients with recurrent KD compared with healthy controls or viral infection controls (p < 0.001). Receiver operating characteristic analysis revealed that high miR-24-3p expression levels may be a potential biomarker for KD diagnosis. In conclusion, we identified miR-24-3p significantly higher in KD patients, which may be a potential diagnostic biomarker for KD.

Mutual Interaction of Clinical Factors and Specific microRNAs to Predict Mild Cognitive Impairment in Patients Receiving Hemodialysis

Cognitive impairment (CI) is not uncommon in dialysis patients. Various factors have been implicated. This study aims to examine mutual interaction of various clinical factors for CI in patients receiving hemodialysis. A total of 48 hemodialysis patients in outpatient clinic were recruited from 2015 to 2017. Demographics, circulating uremic toxin concentrations, miRNA concentrations, and nerve injury protein concentrations were collected. Clinical dementia rating (CDR) scores were used to stratify the functional scores of the patients. Receiver operating characteristic (ROC) analysis was used to evaluate diagnostic test performance for predicting dichotomous results, and cumulative ROC analysis was used to examine the combined contribution of clinical factors. CDR scale 0 included 15 patients (mean age, 59.1 years); CDR > 0.5 included 33 patients (mean age, 64.0 years). On cumulative ROC analysis, the major predictors of mild CI were hemoglobin, age, sex, homocysteine, neuron-specific enolase (NSE), and miR-486. The cumulative area under the curve (AUC) on combining hemoglobin, age, and miR-486 was the highest (0.897, 95% confidence interval 0.806–0.988). Two dichotomized variables reached 81.82% sensitivity and 86.67% specificity, with the likelihood ratio for positive and negative results being 6.14 and 0.21, respectively. In conclusion, hemoglobin, age, and miR-486 display high-degree combined effects on mild CI in patients receiving hemodialysis.

Upregulation of miR-941 in Circulating CD14+ Monocytes Enhances Osteoclast Activation via WNT16 Inhibition in Patients with Psoriatic Arthritis

Psoriatic arthritis (PsA) is a destructive joint disease mediated by osteoclasts. MicroRNAs (miRNAs) regulate several important pathways in osteoclastogenesis. We profiled the expression of miRNAs in CD14+ monocytes from PsA patients and investigated how candidate microRNAs regulate the pathophysiology in osteoclastogenesis. The RNA from circulatory CD14+ monocytes was isolated from PsA patients, psoriasis patients without arthritis (PsO), and healthy controls (HCs). The miRNAs were initially profiled by next-generation sequencing (NGS). The candidate miRNAs revealed by NGS were validated by PCR in 40 PsA patients, 40 PsO patients, and 40 HCs. The osteoclast differentiation and its functional resorption activity were measured with or without RNA interference against the candidate miRNA. The microRNA-941 was selectively upregulated in CD14+ monocytes from PsA patients. Osteoclast development and resorption ability were increased in CD14+ monocytes from PsA patients. Inhibition of miR-941 abrogated the osteoclast development and function while increased the expression of WNT16. After successful treatment, the increased miR-941 expression in CD14+ monocytes from PsA patients was revoked. The expression of miR-941 in CD14+ monocytes is associated with PsA disease activity. MiR-941 enhances osteoclastogenesis in PsA via WNT16 repression. The miR-941 could be a potential biomarker and treatment target for PsA.

Effect of Low-Fat Diet in Obese Mice Lacking Toll-like Receptors

Background: This study aimed at assessing the effect of a low-fat diet (LFD) in obese mice lacking toll–like receptors (Tlr) and understanding the expression and regulation of microRNAs during weight reduction. Methods: C57BL/6, Tlr5^−/−, Tlr2^−/− and Tlr4^−/− mice were used in this study. A group of mice were fed with a high-fat diet (HFD) (58% kcal) for 12 weeks to induce obesity (diet-induced obesity, DIO). Another group that had been fed with HFD for eight weeks (obese mice) were switched to a low-fat diet (LFD) (10.5% kcal) for the next four weeks to reduce their body weight. The control mice were fed with a standard AIN-76A diet for the entire 12 weeks. The body weight of the mice was measured weekly. At the end of the experiment, epididymal fat weight and adipocyte size were measured. The differentially expressed miRNAs in the fat tissue was determined by next-generation sequencing with real-time quantitative reverse transcription polymerase chain reaction (RT–qPCR). Target prediction and functional annotation of miRNAs were performed using miRSystem database. Results: Switching to LFD significantly reduced the body weight and epididymal fat mass in the HFD-fed C57BL/6 and Tlr5^−/− mice but not in Tlr2^−/− and Tlr4^−/− mice. Weight reduction significantly decreased the size of adipocytes in C57BL/6 but not in the Tlr knockout mice. In Tlr2^−/− and Tlr4^−/− mice, feeding with HFD and the subsequent weight reduction resulted in an aberrant miRNA expression in the epididymal fat tissue unlike in C57BL/6 and Tlr5^−/−. However, target prediction and functional annotation by miRSystem database revealed that all the top 10 Kyoto Encyclopedia of Genes and Genomes (KEGG) database pathways of the dysregulated miRNAs during weight reduction in the C57BL/6 mice were also found in the regulated pathways of Tlr5^−/−, Tlr2^−/−, and Tlr4^−/− strains. However, among these pathways, gene sets involved in arginine and proline metabolism and glutathione metabolism were mainly involved in the Tlr knockout mice but not in the C57BL/6 mice. Conclusions: In this study, we demonstrated that feeding of LFD leads to significant body weight reduction in C57BL/6 and Tlr5^−/− mice, but not in Tlr2^−/− and Tlr4^−/− mice. Significant reduction in the size of adipocytes of epididymal fat was only found in C57BL/6, but not in Tlr5^−/−, Tlr2^−/−, and Tlr4^−/− mice. The dysregulated miRNAs in Tlr2^−/− and Tlr4^−/− mice were different from those in C57BL/6 and Tlr5^−/− strains. Among those miRNA-regulated pathways, arginine and proline metabolism as well as glutathione metabolism may have important roles in the Tlr knockout mice rather than in C57BL/6 mice.

Next generation sequencing identifies miRNA-based biomarker panel for lupus nephritis

The symptomatology of lupus nephritis (LN) consists of foamy urine and lower leg edema, as well as such systemic manifestations as oral ulcers, arthralgia/arthritis, and lymphadenopathy. However, these symptoms may appear mild and non-specific. If these symptoms are unrecognized, thus delaying treatment, approximately 10% of LN patients will develop permanent kidney damage and end-stage kidney disease. Therefore, the purpose of this study is to identify a surrogate biomarker for the early detection of LN. In this study, we first adopted next generation sequencing (NGS) in order to screen differential expression levels of microRNA between SLE patients with and without LN. The results of both the NGS and the literature review confirmed the potential of 15 microRNAs through real-time qPCR. We further considered clinical laboratory data for additional analysis. In total, 41 microRNAs demonstrated significant differences through NGS screening. We then verified eight microRNAs from NGS and seven microRNAs from the literature review using the real-time qPCR method in peripheral mononuclear cells. Ultimately, mir-125a-5p, miR-146a-5p, and mir-221-3p were found to be statistically significant not only in the screening study but also in the real-time qPCR verification studies. miR-146a-5p was observed to have a significant correlation with clinical biochemistry markers, as well as to be a surrogate biomarker for the early detection of lupus nephritis. This study is the first to show that the intracellular biomarker miR-146a-5p may serve as a useful specific biomarker for the detection of lupus nephritis among lupus patients in the future, regardless of serum albumin levels and spot urine protein/creatinine ratio.

Next-generation sequencing identifies articular cartilage and subchondral bone miRNAs after ESWT on early osteoarthritis knee

Extracorporeal shockwave therapy (ESWT) has shown chondroprotective effects on the initiation of the osteoarthritis (OA) changes of the rat knee. This study evaluated 69 significant expressed profiles of microRNA (miRNA) in the articular cartilage and subchondral bone after ESWT. There were 118 target genes identified for miRNAs of interest in articular cartilage and 214 target genes in subchondral bone by next generation sequencing (NGS). In principal component analysis (PCA), the relationships of miRNA expression in bone and cartilage were improved after ESWT. Global functional annotation showed that predicted targets were involved in cartilage development, inflammatory and immune response, ion binding, angiogenesis, cell adhesion, cell cycle, transcription and translation, gene expression, NTP binding, signal transduction, collagen fibril organization, apoptotic process, chondrocyte differentiation, cell differentiation, bone development as well as cell proliferation. The miRNAs profile and the target genes were comprehensively surveyed and compared in articular cartilage and subchondral bone of early OA knee before and after ESWT. Our study represents the direct assessment to date of miRNA expression profiling in early OA articular cartilage and subchondral bone. The results provide insights that could contribute to the development of new biomarkers and therapeutic strategies for OA changes and the treatment with ESWT.

Blood-Bourne MicroRNA Biomarker Evaluation in Attention-Deficit/Hyperactivity Disorder of Han Chinese Individuals: An Exploratory Study

Background: Attention-deficit/hyperactivity disorder (ADHD) is a highly genetic neurodevelopmental disorder, and its dysregulation of gene expression involves microRNAs (miRNAs). The purpose of this study was to identify potential miRNAs biomarkers and then use these biomarkers to establish a diagnostic panel for ADHD. Design and methods: RNA samples from white blood cells (WBCs) of five ADHD patients and five healthy controls were combined to create one pooled patient library and one control library. We identified 20 candidate miRNAs with the next-generation sequencing (NGS) technique (Illumina). Blood samples were then collected from a Training Set (68 patients and 54 controls) and a Testing Set (20 patients and 20 controls) to identify the expression profiles of these miRNAs with real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR). We used receiver operating characteristic (ROC) curves and the area under the curve (AUC) to evaluate both the specificity and sensitivity of the probability score yielded by the support vector machine (SVM) model. Results: We identified 13 miRNAs as potential ADHD biomarkers. The ΔCt values of these miRNAs in the Training Set were integrated to create a biomarker model using the SVM algorithm, which demonstrated good validity in differentiating ADHD patients from control subjects (sensitivity: 86.8%, specificity: 88.9%, AUC: 0.94, p < 0.001). The results of the blind testing showed that 85% of the subjects in the Testing Set were correctly classified using the SVM model alignment (AUC: 0.91, p < 0.001). The discriminative validity is not influenced by patients' age or gender, indicating both the robustness and the reliability of the SVM classification model. Conclusion: As measured in peripheral blood, miRNA-based biomarkers can aid in the differentiation of ADHD in clinical settings. Additional studies are needed in the future to clarify the ADHD-associated gene functions and biological mechanisms modulated by miRNAs.

Highly specific real-time quantification of diverse microRNAs in human samples using universal primer set frame

In this study, one group of universal primer set frame, composed by one reverse transcription (RT) primer frame and a pair of quantitative real-time polymerase chain reaction (qRT-PCR) primer frames, was elaborately screened and designed by homebuilt software for sensitive and specific quantification of diverse miRNAs. The universal primer set frame can be applied for multiplex miRNAs detection by simply changing the RT-X part of RT primer frame and RP-Y part of qRT-PCR reverse primer frame based on target sequence. The maximum similarity of RT-Y, RT-Z and qRT-PCR forward primer to the human genome and human transcriptome is less than 76%, ensuring the high specificity in human sample detection. The high sensitivity and broad dynamic linear range of the developed approaches by using designed primer set frame were demonstrated on the in vitro detection of miR-21 and miR-155, with dynamic range of 10 fM to 10 nM and detection limit of 3.74 × 10^-15 M and 5.81 × 10^-15 M for miR-21 and miR-155, respectively. In particular, the developed assays also have high sequence specificity which could clearly discriminate a single base difference in miRNA sequence. The contents of miR-21 and miR-155 in tissue and serum samples have been successfully detected using the developed assays. Results indicated that miR-21 and miR-155 were elevated in cancer tissue and serum specimens than that of normal samples, implying the developed assays hold a great promise for further application in biomedical research and early clinical diagnosis. More importantly, the primer set frame can be universally used in any miRNA investigation.

MiR-30a-5p Inhibits Epithelial-to-Mesenchymal Transition and Upregulates Expression of Tight Junction Protein Claudin-5 in Human Upper Tract Urothelial Carcinoma Cells

The involvement of microRNAs (miRNAs) in cancer development and their potential as prognostic biomarkers are becoming increasingly known. However, the signature of miRNAs and their regulatory roles in tumorigenesis of upper tract urothelial carcinoma (UTUC) remain to be elucidated. This study aimed to profile the miRNA expression pattern in UTUC tumor tissues and identify candidate miRNAs with prognostic and/or therapeutic functions. Methods and Results: We collected 22 UTUC tissue and adjacent normal tissues samples from patients who underwent nephroureterectomy. The miRNAs signatures of three selected UTUC samples using next-generation sequencing showed that miR-30a-5p was significantly downregulated in UTUC tumors compared to adjacent normal tissues. The differentially-expressed miRNAs were specifically validated by quantitative real-time polymerase chain reaction. In addition, the miRNA expression signatures were analyzed with the transcriptome profile characterized by microarray. Further in vitro studies indicated that overexpression of miR-30a-5p significantly suppressed proliferation, migration, and epithelial-to-mesenchymal transition (EMT) in cultured BFTC-909 UTUC cells. As a potential target gene of miR-30a-5p in the tight junction pathway suggested by the pathway enrichment analysis, the reduced expression of tight junction protein claudin-5 in UTUC cells was demonstrated to be upregulated by miR-30a-5p genetic delivery. Conclusions: Taken together, our findings demonstrated that miR-30a-5p inhibits proliferation, metastasis, and EMT, and upregulates the expression of tight junction claudin-5 in UTUC cells. Thus, miR-30a-5p may provide a promising therapeutic strategy for UTUC treatment.

Major methylation alterations on the CpG markers of inflammatory immune associated genes after IVIG treatment in Kawasaki disease

Background

Kawasaki disease (KD) is an autoimmune disease preferentially attacking children younger than five years worldwide. So far, the principal treatment to KD is the administration of Intravenous immunoglobulin (IVIG). Although DNA methylation plays important regulation roles in diseases, few studies investigated the regulation roles of DNA methylation in KD.

Methods

In this study, we focused not only on the DNA methylation alterations resulted from KD onset but also on DNA methylation alterations resulted from IVIG administration. To do so, we investigated the effects of KD’s onset and IVIG administration on CpG marker’s methylation alterations.

Results

We first found that DNA methylation alterations reflecting disease onset or IVIG administration are contributed mainly by the CpG markers on autosomes. In addition, we also demonstrated that some CpG markers carry methylation alteration among samples, forcing the expression abundance of the downstream genes to be also altered and negatively correlated with methylation profile. Finally, compared with KD’s onset, IVIG administration brings stronger impact on methylation alteration. And, such alterations were conducted mainly by hyper-methylating CpG markers, forcing the corresponding genes to keep lower expression levels. Moreover, the genes regulated by the altered CpG markers with IVIG administration are enriched in the pathways associated with inflammatory immune response.

Conclusions

In summary, our result provides researchers with another way into the regulation mechanism through which IVIG represses excessive inflammatory responses.

Early and late effects of prenatal corticosteroid treatment on the microRNA profiles of lung tissue in rats

Glucocorticoids have been administered to mothers at risk of premature delivery to induce maturation of preterm fetal lungs and prevent the development of respiratory distress syndrome. Micro (mi)RNAs serve various crucial functions in cell proliferation, differentiation and organ development; however, few studies have demonstrated an association between miRNAs and lung development. The aim of the present study was to investigate alterations in the miRNA profiles of rat lung tissue following prenatal glucocorticoid therapy for fetal lung development. The differences in miRNA expression profiles were compared between postnatal days 7 (D7) and 120 (D120) rat lung tissues, followed by validation using reverse transcription-quantitative polymerase chain reaction. The miRNA profiles of rat lung tissues following prenatal dexamethasone (DEX) therapy were also investigated. miRNAs with 2-fold changes were selected for further analysis. At D120, 6 upregulated and 6 downregulated miRNAs were detected, compared with D7. Among these differentially expressed miRNAs, miR-101-3p and miR-99b-5p were associated with the lowest and highest expressions of miRNA at D7, respectively. A limited impact on the miRNA profiles of rat lung tissues was observed following prenatal DEX treatment, which may help to further clarify the mechanisms underlying normal lung development. However, the results of the present study cannot entirely elucidate the effects of prenatal DEX treatment on the lung development of premature infants, and further studies investigating the impact of prenatal corticosteroids on fetal lung miRNA profiles are required.

Evaluation and application of the strand-specific protocol for next-generation sequencing

Next-generation sequencing (NGS) has become a powerful sequencing tool, applied in a wide range of biological studies. However, the traditional sample preparation protocol for NGS is non-strand-specific (NSS), leading to biased estimates of expression for transcripts overlapped at the antisense strand. Strand-specific (SS) protocols have recently been developed. In this study, we prepared the same RNA sample by using the SS and NSS protocols, followed by sequencing with Illumina HiSeq platform. Using real-time quantitative PCR as a standard, we first proved that the SS protocol more precisely estimates gene expressions compared with the NSS protocol, particularly for those overlapped at the antisense strand. In addition, we also showed that the sequence reads from the SS protocol are comparable with those from conventional NSS protocols in many aspects. Finally, we also mapped a fraction of sequence reads back to the antisense strand of the known genes, originally without annotated genes located. Using sequence assembly and PCR validation, we succeeded in identifying and characterizing the novel antisense genes. Our results show that the SS protocol performs more accurately than the traditional NSS protocol and can be applied in future studies.