Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression

Elimination of ligation dependent artifacts in T4 RNA ligase to achieve high efficiency and low bias microRNA capture

Adapter ligation is a critical first step in many microRNA analysis methods including microarray, qPCR, and sequencing. Previous studies have shown that ligation bias can have dramatic effects on both the fidelity of expression profiles and reproducibility across samples. We have developed a method for high efficiency and low bias microRNA capture by 3′ adapter ligation using T4 RNA ligase that does not require pooled adapters. Using a panel of 20 microRNA, we investigated the effects of ligase type, PEG concentration, ligase amount, adapter concentration, incubation time, incubation temperature, and adapter design on capture efficiency and bias. Of these factors, high PEG% was found to be critical in suppressing ligation bias. We obtained high average capture efficiency and low CV across the 20 microRNA panel, both in idealized buffer conditions (86%±10%) and total RNA spiking conditions (64%±17%). We demonstrate that this method is reliable across microRNA species that previous studies have had difficulty capturing and that our adapter design performs significantly better than the common adapter designs. Further, we demonstrate that the optimization methodology must be specifically designed for minimizing bias in order to obtain the ideal reaction parameters.

A detailed analysis of next generation sequencing reads of microRNA expression in Barrett's esophagus: absolute versus relative quantification

BACKGROUND

Next generation sequencing (NGS) is a state of the art technology for microRNA (miRNA) analysis. The quantitative interpretation of the primary output of NGS i.e. the read counts for a miRNA sequence that can vary by several orders of magnitude (1 to 107) remains incompletely understood.

FINDINGS

NGS (SOLiD 3 technology) was performed on biopsies from 6 Barrett's esophagus (BE) and 5 Gastroesophageal Reflux Disease (GERD) patients. Read sequences were aligned to miRBase 18.0. Differential expression analysis was adjusted for false discovery rate of 5%. Quantitative real-time polymerase chain reaction (qRT-PCR) was performed for 36 miRNA in a validation cohort of 47 patients (27 BE and 20 GERD). Correlation coefficients, accuracy, precision and recall of NGS compared to qRT-PCR were calculated. Increase in NGS reads was associated with progressively lower Cq values, p < 0.05. Although absolute quantification between NGS reads and Cq values correlated modestly: -0.38, p = 0.01 for BE and -0.32, p = 0.05 for GERD, relative quantification (fold changes) of miRNA expression between BE &GERD by NGS correlated highly with qRT-PCR 0.86, p = 2.45E-11. Fold change correlations were unaffected when different thresholds of NGS read counts were compared (>1000 vs. <1000, >500 vs. <500 and >100 vs. <100). The accuracy, precision and recall of NGS to label a miRNA as differentially expressed were 0.71, 0.88 and 0.74 respectively.

CONCLUSION

Absolute NGS reads correlated modestly with qRT-PCR but fold changes correlated highly. NGS is robust at relative but not absolute quantification of miRNA levels and accurate for high-throughput identification of differentially expressed miRNA.

Circulating microRNA biomarkers for glioma and predicting response to therapy

The need for glioma biomarkers with improved sensitivity and specificity has sparked research into short non-coding RNA known as microRNA (miRNA). Altered miRNA biogenesis and expression in glioma plays a vital role in important signaling pathways associated with a range of tumor characteristics including gliomagenesis, invasion, and malignancy. This review will discuss current research into the role of miRNA in glioma and altered miRNA expression in biofluids as candidate biomarkers with a particular focus on glioblastoma, the most malignant form of glioma. The isolation and characterization of miRNA using cellular and molecular biology techniques from the circulation of glioma patients could potentially be used for improved diagnosis, prognosis, and treatment decisions. We aim to highlight the links between research into miRNA function, their use as biomarkers, and how these biomarkers can be used to predict response to therapy. Furthermore, increased understanding of miRNA in glioma biology through biomarker research has led to the development of miRNA therapeutics which could restore normal miRNA expression and function and improve the prognosis of glioma patients. A panel of important miRNA biomarkers for glioma in various biofluids discovered to date has been summarized here. There is still a need, however, to standardize techniques for biomarker characterization to bring us closer to clinically relevant miRNA-based diagnostic and therapeutic signatures. A clinically validated biomarker panel has potential to improve time to diagnosis, predicting response to treatment and ultimately the prognosis of glioma patients.

miRNA dysregulation in cancer: towards a mechanistic understanding

It is now well known that gene expression is intricately regulated inside each cell especially in mammals. There are multiple layers of gene regulation active inside a cell at a given point of time. Gene expression is regulated post-transcriptionally by microRNAs and other factors. Mechanistically, microRNAs are known to bind to the 3’ UTR of mRNAs and cause repression of gene expression and the number of known microRNAs continues to increase every day. Dysregulated microRNA signatures in different types of cancer are being uncovered consistently implying their importance in cellular homeostasis. However when studied in isolation in mouse models, clear-cut cellular and molecular mechanisms have been described only for a select few microRNAs. What is the reason behind this discrepancy? Are microRNAs small players in gene regulation helping only to fine tune gene expression? Or are their roles tissue and cell type-specific with single-cell level effects on mRNA expression and microRNA threshold levels? Or does it all come down to the technical limitations of high-throughput techniques, resulting in false positive results? In this review, we will assess the challenges facing the field and potential avenues for resolving the cellular and molecular mechanisms of these small but important regulators of gene expression.

Analysis options for high-throughput sequencing in miRNA expression profiling

Background

Recently high-throughput sequencing (HTS) using next generation sequencing techniques became useful in digital gene expression profiling.

Our study introduces analysis options for HTS data based on mapping to miRBase or counting and grouping of identical sequence reads. Those approaches allow a hypothesis free detection of miRNA differential expression.

Methods

We compare our results to microarray and qPCR data from one set of RNA samples. We use Illumina platforms for microarray analysis and miRNA sequencing of 20 samples from benign follicular thyroid adenoma and malignant follicular thyroid carcinoma. Furthermore, we use three strategies for HTS data analysis to evaluate miRNA biomarkers for malignant versus benign follicular thyroid tumors.

Results

High correlation of qPCR and HTS data was observed for the proposed analysis methods. However, qPCR is limited in the differential detection of miRNA isoforms. Moreover, we illustrate a much broader dynamic range of HTS compared to microarrays for small RNA studies. Finally, our data confirm hsa-miR-197-3p, hsa-miR-221-3p, hsa-miR-222-3p and both hsa-miR-144-3p and hsa-miR-144-5p as potential follicular thyroid cancer biomarkers.

Conclusions

Compared to microarrays HTS provides a global profile of miRNA expression with higher specificity and in more detail. Summarizing of HTS reads as isoform groups (analysis pipeline B) or according to functional criteria (seed analysis pipeline C), which better correlates to results of qPCR are promising new options for HTS analysis. Finally, data opens future miRNA research perspectives for HTS and indicates that qPCR might be limited in validating HTS data in detail.

Label-free quantification of microRNAs using ligase-assisted sandwich hybridization on a DNA microarray

MicroRNAs (miRNAs) can be used as biomarkers for cancer and other human diseases; therefore, high-throughput and reliable miRNA-quantification methods are required to exploit these markers for diagnostic testing. In this report, we describe the construction of a platform for miRNA-quantification using ligase-assisted sandwich hybridization (LASH) without miRNA-labeling. T4 DNA ligase was used to compensate for the low affinity between miRNAs and two short complementary DNA probes, and it improved the hybridization yield ∼50,000 times. The LASH assay enabled synthesized miR-143 to be quantified at concentrations ranging from 30 fM to 30 pM. The LASH assay could also quantify endogenous miR-143 released from cultured cells as well as some miRNAs in total RNAs derived from blood. Furthermore, multi-color detection enabled us to distinguish between the highly homologous miR-141 and miR-200a. This simple label-free quantification technique is an easy-to-use approach that can be applied to disease diagnosis.

MicroRNA dysregulation in spinal cord injury: causes, consequences and therapeutics

Trauma to the spinal cord causes permanent disability to more than 180,000 people every year worldwide. The initial mechanical damage triggers a complex set of secondary events involving the neural, vascular, and immune systems that largely determine the functional outcome of the spinal cord injury (SCI). Cellular and biochemical mechanisms responsible for this secondary injury largely depend on activation and inactivation of specific gene programs. Recent studies indicate that microRNAs function as gene expression switches in key processes of the SCI. Microarray data from rodent contusion models reveal that SCI induces changes in the global microRNA expression patterns. Variations in microRNA abundance largely result from alterations in the expression of the cells at the damaged spinal cord. However, microRNA expression levels after SCI are also influenced by the infiltration of immune cells to the injury site and the death and migration of specific neural cells after injury. Evidences on the role of microRNAs in the SCI pathophysiology have come from different sources. Bioinformatic analysis of microarray data has been used to identify specific variations in microRNA expression underlying transcriptional changes in target genes, which are involved in key processes in the SCI. Direct evidences on the role of microRNAs in SCI are scarcer, although recent studies have identified several microRNAs (miR-21, miR-486, miR-20) involved in key mechanisms of the SCI such as cell death or astrogliosis, among others. From a clinical perspective, different evidences make clear that microRNAs can be potent therapeutic tools to manipulate cell state and molecular processes in order to enhance functional recovery. The present article reviews the actual knowledge on how injury affects microRNA expression and the meaning of these changes in the SCI pathophysiology, to finally explore the clinical potential of microRNAs in the SCI.

Profiling of estrogen-regulated microRNAs in breast cancer cells

Estrogen plays vital roles in mammary gland development and breast cancer progression. It mediates its function by binding to and activating the estrogen receptors (ERs), ERα, and ERβ. ERα is frequently upregulated in breast cancer and drives the proliferation of breast cancer cells. The ERs function as transcription factors and regulate gene expression. Whereas ERα's regulation of protein-coding genes is well established, its regulation of noncoding microRNA (miRNA) is less explored. miRNAs play a major role in the post-transcriptional regulation of genes, inhibiting their translation or degrading their mRNA. miRNAs can function as oncogenes or tumor suppressors and are also promising biomarkers. Among the miRNA assays available, microarray and quantitative real-time polymerase chain reaction (qPCR) have been extensively used to detect and quantify miRNA levels. To identify miRNAs regulated by estrogen signaling in breast cancer, their expression in ERα-positive breast cancer cell lines were compared before and after estrogen-activation using both the µParaflo-microfluidic microarrays and Dual Labeled Probes-low density arrays. Results were validated using specific qPCR assays, applying both Cyanine dye-based and Dual Labeled Probes-based chemistry. Furthermore, a time-point assay was used to identify regulations over time. Advantages of the miRNA assay approach used in this study is that it enables a fast screening of mature miRNA regulations in numerous samples, even with limited sample amounts. The layout, including the specific conditions for cell culture and estrogen treatment, biological and technical replicates, and large-scale screening followed by in-depth confirmations using separate techniques, ensures a robust detection of miRNA regulations, and eliminates false positives and other artifacts. However, mutated or unknown miRNAs, or regulations at the primary and precursor transcript level, will not be detected. The method presented here represents a thorough investigation of estrogen-mediated miRNA regulation.

Methodological challenges in utilizing miRNAs as circulating biomarkers

MicroRNAs (miRNAs) have emerged as important regulators in the post-transcriptional control of gene expression. The discovery of their presence not only in tissues but also in extratissular fluids, including blood, urine and cerebro-spinal fluid, together with their changes in expression in various pathological conditions, has implicated these extracellular miRNAs as informative biomarkers of disease. However, exploiting miRNAs in this capacity requires methodological rigour. Here, we report several key procedural aspects of miRNA isolation from plasma and serum, as exemplified by research in cardiovascular and pulmonary diseases. We also highlight the advantages and disadvantages of various profiling methods to determine the expression levels of plasma- and serum-derived miRNAs. Attention to such methodological details is critical, as circulating miRNAs become diagnostic tools for various human diseases.

Convergent microRNA actions coordinate neocortical development

Neocortical development is a complex process that, at the cellular level, involves tight control of self-renewal, cell fate commitment, survival, differentiation and delamination/migration. These processes require, at the molecular level, the precise regulation of intrinsic signaling pathways and extrinsic factors with coordinated action in a spatially and temporally specific manner. Transcriptional regulation plays an important role during corticogenesis; however, microRNAs (miRNAs) are emerging as important post-transcriptional regulators of various aspects of central nervous system development. miRNAs are a class of small, single-stranded noncoding RNA molecules that control the expression of the majority of protein coding genes (i.e., targets). How do different miRNAs achieve precise control of gene networks during neocortical development? Here, we critically review all the miRNA–target interactions validated in vivo, with relevance to the generation and migration of pyramidal-projection glutamatergic neurons, and for the initial formation of cortical layers in the embryonic development of rodent neocortex. In particular, we focus on convergent miRNA actions, which are still a poorly understood layer of complexity in miRNA signaling, but potentially one of the keys to disclosing how miRNAs achieve the precise coordination of complex biological processes such as neocortical development.

MicroRNA responses to focal cerebral ischemia in male and female mouse brain

Stroke occurs with greater frequency in men than in women across diverse ethnic backgrounds and nationalities. Work from our lab and others have revealed a sex-specific sensitivity to cerebral ischemia whereby males exhibit a larger extent of brain damage resulting from an ischemic event compared to females. Previous studies revealed that microRNA (miRNA) expression is regulated by cerebral ischemia in males; however, no studies to date have examined the effect of ischemia on miRNA responses in females. Thus, we examined miRNA responses in male and female brain in response to cerebral ischemia using miRNA arrays. These studies revealed that in male and female brains, ischemia leads to both a universal miRNA response as well as a sexually distinct response to challenge. Target prediction analysis of the miRNAs increased in male or female ischemic brain reveal sex-specific differences in gene targets and protein pathways. These data support that the mechanisms underlying sexually dimorphic responses to cerebral ischemia includes distinct changes in miRNAs in male and female brain, in addition to a miRNA signature response to ischemia that is common to both.

The role of microRNAs in melanoma

Melanoma is the most dangerous form of skin cancer, being largely resistant to conventional therapies at advanced stages. Understanding the molecular mechanisms behind this disease might be the key for the development of novel therapeutic strategies. MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally control gene expression, thereby regulating various cellular signaling pathways involved in the initiation and progression of different cancer types, including melanoma. In this review, we summarize approaches for the identification of candidate miRNAs and their target genes and review the functions of miRNAs in melanoma. Finally, we highlight the recent progress in pre-clinical use of miRNAs as prognostic markers and therapeutic targets.

Differential expression of miRNAs in Brassica napus root following infection with Plasmodiophora brassicae

Canola (oilseed rape, Brassica napus L.) is susceptible to infection by the biotrophic protist Plasmodiophora brassicae, the causal agent of clubroot. To understand the roles of microRNAs (miRNAs) during the post-transcriptional regulation of disease initiation and progression, we have characterized the changes in miRNA expression profiles in canola roots during clubroot disease development and have compared these to uninfected roots. Two different stages of clubroot development were targeted in this miRNA profiling study: an early time of 10-dpi for disease initiation and a later 20-dpi, by which time the pathogen had colonized the roots (as evident by visible gall formation and histological observations). P. brassicae responsive miRNAs were identified and validated by qRT-PCR of miRNAs and the subsequent validation of the target mRNAs through starBase degradome analysis, and through 5' RLM-RACE. This study identifies putative miRNA-regulated genes with roles during clubroot disease initiation and development. Putative target genes identified in this study included: transcription factors (TFs), hormone-related genes, as well as genes associated with plant stress response regulation such as cytokinin, auxin/ethylene response elements. The results of our study may assist in elucidating the role of miRNAs in post-transcriptional regulation of target genes during disease development and may contribute to the development of strategies to engineer durable resistance to this important phytopathogen.

Expression of miR-15/107 family microRNAs in human tissues and cultured rat brain cells

The miR-15/107 family comprises a group of 10 paralogous microRNAs (miRNAs), sharing a 5′ AGCAGC sequence. These miRNAs have overlapping targets. In order to characterize the expression of miR-15/107 family miRNAs, we employed customized TaqMan Low-Density micro-fluid PCR-array to investigate the expression of miR-15/107 family members, and other selected miRNAs, in 11 human tissues obtained at autopsy including the cerebral cortex, frontal cortex, primary visual cortex, thalamus, heart, lung, liver, kidney, spleen, stomach and skeletal muscle. miR-103, miR-195 and miR-497 were expressed at similar levels across various tissues, whereas miR-107 is enriched in brain samples. We also examined the expression patterns of evolutionarily conserved miR-15/107 miRNAs in three distinct primary rat brain cell preparations (enriched for cortical neurons, astrocytes and microglia, respectively). In primary cultures of rat brain cells, several members of the miR-15/107 family are enriched in neurons compared to other cell types in the central nervous system (CNS). In addition to mature miRNAs, we also examined the expression of precursors (pri-miRNAs). Our data suggested a generally poor correlation between the expression of mature miRNAs and their precursors. In summary, we provide a detailed study of the tissue and cell type-specific expression profile of this highly expressed and phylogenetically conserved family of miRNA genes.

Transient Maternal Hypothyroxinemia Potentiates the Transcriptional Response to Exogenous Thyroid Hormone in the Fetal Cerebral Cortex Before the Onset of Fetal Thyroid Function: A Messenger and MicroRNA Profiling Study

Thyroid hormone (TH) is essential for brain development both before and after birth. We have used gene expression microarrays to identify TH-regulated genes in the fetal cerebral cortex prior to the onset of fetal thyroid function to better understand the role of TH in early cortical development. TH levels were transiently manipulated in pregnant mice by treatment with goitrogens from gestational day (GD) 13-16 and/or by injection of TH 12 h before sacrifice on GD 16. The transcriptional response to exogenous TH in the GD 16 fetal cortex was potentiated by transient goitrogen treatment, suggesting that the hypothyroxinemic brain is a different substrate upon which TH can act, or that robust compensatory mechanisms are induced by transient hypothyroxinemia. Several known TH-responsive genes were identified including Klf9, and several novel TH-responsive genes such as Appbp2, Ppap2b, and Fgfr1op2 were identified in which TH response elements were confirmed. We also identified specific microRNAs whose expression in the fetal cortex was affected by TH treatment, and determined that Ppap2b and Klf9 are the target genes of miR-16 and miR-106, respectively. Thus, a complex redundant functional network appears to coordinate TH-mediated gene expression in the developing brain.

Methods and matrices: approaches to identifying miRNAs for nasopharyngeal carcinoma

Background

Nasopharyngeal carcinoma (NPC) is a solid tumor of the head and neck. Multimodal therapy is highly effective when NPC is detected early. However, due to the location of the tumor and the absence of clinical signs, early detection is difficult, making a biomarker for the early detection of NPC a priority. The dysregulation of small non-coding RNAs (miRNAs) during carcinogenesis is the focus of much current biomarker research. Herein, we examine several miRNA discovery methods using two sample matrices to identify circulating miRNAs (c-miRNAs) associated with NPC.

Methods

We tested two miRNA discovery workflows on two sample sources for miRNAs associated with NPC. In the first workflow, we assumed that NPC tumor tissue would be enriched for miRNAs, so we compared miRNA expression in FFPE from NPC cases and controls using microarray and RNA-Seq technologies. Candidate miRNAs from both technologies were verified by qPCR in FFPE and sera from an independent NPC sample set. In a second workflow, we directly interrogated NPC case and control sera by RNA-Seq for c-miRNAs associated with NPC, with candidate c-miRNAs verified by qPCR in the sera from the same independent NPC sample set.

Results

Both microarray and RNA-Seq narrowed the miRNA signature to 1-5% of the known mature human miRNAs. Moreover, these two methods produced similar results when applied to the same sample type (FFPE), with RNA-Seq additionally indicating “unknown” miRNAs associated with NPC. However, we found different miRNA profiles in NPC sera compared to FFPE using RNA-Seq, with the few overlapping miRNAs found to be significantly up-regulated in FFPE significantly down-regulated in sera (and vice versa). Despite the different miRNA profiles found in FFPE and sera, both profiles strongly associated with NPC, providing two potential sources for biomarker signatures for NPC.

Conclusions

We determined that the direct interrogation of sera by RNA-Seq was the most informative method for identifying a c-miRNA signature associated with NPC. We also showed that there are different miRNA expression profiles associated with NPC for tumor tissue and sera. These results reflect on the methods and meaning of miRNA biomarkers for NPC in tissue and peripheral blood.

Identification and expression profiling of Vigna mungo microRNAs from leaf small RNA transcriptome by deep sequencing

MicroRNAs (miRNAs) represent a class of small non-coding RNA molecules that play a crucial role in post-transcriptional gene regulation. Several conserved and species-specific miRNAs have been characterized to date, predominantly from the plant species whose genome is well characterized. However, information on the variability of these regulatory RNAs in economically important but genetically less characterized crop species are limited. Vigna mungo is an important grain legume, which is grown primarily for its protein-rich edible seeds. miRNAs from this species have not been identified to date due to lack of genome sequence information. To identify miRNAs from V. mungo, a small RNA library was constructed from young leaves. High-throughput Illumina sequencing technology and bioinformatic analysis of the small RNA reads led to the identification of 66 miRNA loci represented by 45 conserved miRNAs belonging to 19 families and eight non-conserved miRNAs belonging to seven families. Besides, 13 novel miRNA candidates in V. mungo were also identified. Expression patterns of selected conserved, non-conserved, and novel miRNA candidates have been demonstrated in leaf, stem, and root tissues by quantitative polymerase chain reaction, and potential target genes were predicted for most of the conserved miRNAs. This information offers genomic resources for better understanding of miRNA mediated post-transcriptional gene regulation.

Analysis of microRNA expression using machine learning

The systematic analysis of miRNA expression and its potential mRNA targets constitutes a basal objective in miRNA research in addition to miRNA gene detection and miRNA target prediction. In this chapter we address methodical issues of miRNA expression analysis using self-organizing maps (SOM), a neural network machine learning algorithm with strong visualization and second-level analysis capabilities widely used to categorize large-scale, high-dimensional data. We shortly review selected experimental and theoretical aspects of miRNA expression analysis. Then, the protocol of our SOM method is outlined with special emphasis on miRNA/mRNA coexpression. The method allows extracting differentially expressed RNA transcripts, their functional context, and also characterization of global properties of expression states and profiles. In addition to the separate study of miRNA and mRNA expression landscapes, we propose the combined analysis of both entities using a covariance SOM.

High-throughput approaches for microRNA expression analysis

Profiling microRNA (miRNA) expression is of widespread interest due to their critical roles in diverse biological processes, including development, cell proliferation, differentiation, and apoptosis. Profiling can be achieved via three major methods: amplification-based (real-time quantitative PCR, qRT-PCR), hybridization-based (microarrays), and sequencing-based (next-generation sequencing (NGS)) technologies. The gold standard is qRT-PCR and serves as a platform for single reverse PCR amplification experiments and for a large number of miRNAs in parallel, both by multiplexing and plate based arrays. Currently, qRT-PCR is used for the validation of miRNA profiling results from other platforms. Hybridization based miRNA profiling by microarrays has become a widely used method especially for biomarker and therapeutic target identification. The data obtained from microarrays also enables functional prediction of miRNAs by correlating miRNA expression patterns to corresponding mRNA and protein profiles. Additionally, miRNA profiling strategies based on deep sequencing allow both the identification of novel miRNAs and relative quantification of miRNAs. Each miRNA profiling strategy has specific strengths and challenges that have to be considered depending on the nature of the research context.In this chapter the high-throughput approaches that can be applied to microRNA profiling are discussed starting from small-scale qRT-PCR technology to a wider one, NGS.

MicroRNA methodology: advances in miRNA technologies

There is an emerging trend in microRNA research and thus substantial progress in microRNA technologies. In this chapter we provide insights into the main microRNA specific methodologies and critical steps of microRNA expression profiling, target gene identification, and functional confirmation of microRNA effects up to in vivo application of antagomirs.

Conserved miRNAs are candidate post-transcriptional regulators of developmental arrest in free-living and parasitic nematodes

Animal development is complex yet surprisingly robust. Animals may develop alternative phenotypes conditional on environmental changes. Under unfavorable conditions, Caenorhabditis elegans larvae enter the dauer stage, a developmentally arrested, long-lived, and stress-resistant state. Dauer larvae of free-living nematodes and infective larvae of parasitic nematodes share many traits including a conserved endocrine signaling module (DA/DAF-12), which is essential for the formation of dauer and infective larvae. We speculated that conserved post-transcriptional regulatory mechanism might also be involved in executing the dauer and infective larvae fate. We used an unbiased sequencing strategy to characterize the microRNA (miRNA) gene complement in C. elegans, Pristionchus pacificus, and Strongyloides ratti. Our study raised the number of described miRNA genes to 257 for C. elegans, tripled the known gene set for P. pacificus to 362 miRNAs, and is the first to describe miRNAs in a Strongyloides parasite. Moreover, we found a limited core set of 24 conserved miRNA families in all three species. Interestingly, our estimated expression fold changes between dauer versus nondauer stages and infective larvae versus free-living stages reveal that despite the speed of miRNA gene set evolution in nematodes, homologous gene families with conserved “dauer-infective” expression signatures are present. These findings suggest that common post-transcriptional regulatory mechanisms are at work and that the same miRNA families play important roles in developmental arrest and long-term survival in free-living and parasitic nematodes.

Assessing an improved protocol for plasma microRNA extraction

The first step in biomarkers discovery is to identify the best protocols for their purification and analysis. This issue is critical when considering peripheral blood samples (plasma and serum) that are clinically interesting but meet several methodological problems, mainly complexity and low biomarker concentration. Analysis of small molecules, such as circulating microRNAs, should overcome these disadvantages. The present study describes an optimal RNA extraction method of microRNAs from human plasma samples. Different reagents and commercially available kits have been analyzed, identifying also the best pre-analytical conditions for plasma isolation. Between all of them, the column-based approaches were shown to be the most effective. In this context, miRNeasy Serum/Plasma Kit (from Qiagen) rendered more concentrated RNA, that was better suited for microarrays studies and did not require extra purification steps for sample concentration and purification than phenol based extraction methods. We also present evidences that the addition of low doses of an RNA carrier before starting the extraction process improves microRNA purification while an already published carrier dose can result in significant bias over microRNA profiles. Quality controls for best protocol selection were developed by spectrophotometry measurement of contaminants and microfluidics electrophoresis (Agilent 2100 Bioanalyzer) for RNA integrity. Selected donor and patient plasma samples and matched biopsies were tested by Affymetrix microarray technology to compare differentially expressed microRNAs. In summary, this study defines an optimized protocol for microRNA purification from human blood samples, increasing the performance of assays and shedding light over the best way to discover and use these biomarkers in clinical practice.

A new microarray platform for whole-genome expression profiling of Mycobacterium tuberculosis

Microarrays have allowed gene expression profiling to progress from the gene level to the genome level, and oligonucleotide microarrays have become the platform of choice for large-scale, targeted gene expression studies. cDNA arrays and spotted oligonucleotide arrays have gradually given way to in situ synthesized oligonucleotide-based DNA microarrays for whole-genome expression profiling. With the identification of new coding and regulatory sequences, it is imperative that microarrays be updated to enable more complete expression profiling of genomes. We report here a new in situ synthesized oligonucleotide-based microarray platform for Mycobacterium tuberculosis that has been updated for the latest genome information and incorporates hitherto unannotated genes with described biological functions. This microarray has greater coverage of mycobacterial genes than any other array reported to date. We have also evaluated different labeled-target preparation methods and hybridization conditions for this new microarray to obtain high quality data and reproducible results. The new design has been rigorously validated for its specificity and performance using samples isolated from mycobacteria grown under different environment conditions. Further, the quality of the generated data has been compared with published data and is superior to that obtained using spotted oligonucleotide microarrays.

Comprehensive annotation of microRNA expression profiles

BACKGROUND

MicroRNAs (miRNAs) regulate many biological processes by post-translational gene silencing. Analysis of miRNA expression profiles is a reliable method for investigating particular biological processes due to the stability of miRNA and the development of advanced sequencing methods. However, this approach is limited by the broad specificity of miRNAs, which may target several mRNAs.

RESULT

In this study, we developed a method for comprehensive annotation of miRNA array or deep sequencing data for investigation of cellular biological effects. Using this method, the specific pathways and biological processes involved in Alzheimer's disease were predicted with high correlation in four independent samples. Furthermore, this method was validated for evaluation of cadmium telluride (CdTe) nanomaterial cytotoxicity. As a result, apoptosis pathways were selected as the top pathways associated with CdTe nanoparticle exposure, which is consistent with previous studies.

CONCLUSIONS

Our findings contribute to the validation of miRNA microarray or deep sequencing results for early diagnosis of disease and evaluation of the biological safety of new materials and drugs.

MicroRNA expression signatures of stage, grade, and progression in clear cell RCC

Clear cell RCC is the most common, and more likely to metastasize, of the three main histological types of RCC. Pathologic stage is the most important prognostic indicator and nuclear grade can predict outcome within stages of localized RCC. Epithelial tumors are thought to accumulate a series of genetic and epigenetic changes as they progress through well-defined clinical and histopathological changes. MicroRNAs (miRNAs) are involved in the regulation of mRNA expression from many human genes and miRNA expression is dysregulated in cancer. To better understand the contribution of dysregulated miRNA expression to the progression and biology of ccRCC, we examined the differences in expression levels of 723 human miRNAs through a series of analyses by stage, grade, and disease progression status in a large series of 94 ccRCC. We found a consistent signature that included significant upregulation of miR-21-5p, 142-3p, let-7g-5p, let-7i-5p and 424-5p, as well as downregulation of miR-204-5p, to be associated with ccRCC of high stage, or high grade, or progression. Discrete signatures associated with each of stage, grade, or progression were also identified. The let-7 family was significantly downregulated in ccRCC compared with normal renal parenchyma. Expression of the 6 most significantly differentially expressed miRNAs between ccRCC was verified by stem-loop qRT-PCR. Pathways predicted as targets of the most significantly dysregulated miRNAs included signaling, epithelial cancers, metabolism, and epithelial to mesenchymal transition. Our studies help to further elucidate the biology underlying the progression of ccRCC and identify miRNAs for potential translational application.

Dysregulated microRNAs in amyotrophic lateral sclerosis microglia modulate genes linked to neuroinflammation

MicroRNAs (miRNAs) regulate gene expression at post-transcriptional level and are key modulators of immune system, whose dysfunction contributes to the progression of neuroinflammatory diseaseas such as amyotrophic lateral sclerosis (ALS), the most widespread motor neuron disorder. ALS is a non-cell-autonomous disease targeting motor neurons and neighboring glia, with microgliosis directly contributing to neurodegeneration. As limited information exists on miRNAs dysregulations in ALS, we examined this topic in primary microglia from superoxide dismutase 1-G93A mouse model. We compared miRNAs transcriptional profiling of non-transgenic and ALS microglia in resting conditions and after inflammatory activation by P2X7 receptor agonist. We identified upregulation of selected immune-enriched miRNAs, recognizing miR-22, miR-155, miR-125b and miR-146b among the most highly modulated. We proved that miR-365 and miR-125b interfere, respectively, with the interleukin-6 and STAT3 pathway determining increased tumor necrosis factor alpha (TNFα) transcription. As TNFα directly upregulated miR-125b, and inhibitors of miR-365/miR-125b reduced TNFα transcription, we recognized the induction of miR-365 and miR-125b as a vicious gateway culminating in abnormal TNFα release. These results strengthen the impact of miRNAs in modulating inflammatory genes linked to ALS and identify specific miRNAs as pathogenetic mechanisms in the disease.

Analysis of organ-enriched microRNAs in plasma as an approach to development of Universal Screening Test: feasibility study

Background

Early disease detection with a minimally invasive screening test will significantly increase effectiveness and decrease the cost of treatment. Here we propose a framework of a novel approach – Universal Screening Test (UST) for the detection of pathological processes in a particular organ system, organ, or tissue by RT-qPCR analysis of circulating cell-free miRNAs in plasma. As the first step towards assessing the feasibility of this concept, the present study was designed to analyze whether the same microRNAs (miRNAs) can detect various diseases of a particular organ system.

Methods

RNA was extracted from plasma using Trizol treatment and silica binding. Levels of miRNAs were measured by single target RT-qPCR. The following innovations have been tested and proven effective: (i) the use of organ system/organ/tissue-enriched miRNAs; (ii) the use of miRNAs associated with broad disease categories, such as cancer and inflammation, in combination with the organ-enriched miRNAs; and (iii) the use of “miRNA pairs” for selecting miRNA combinations with the highest sensitivity and specificity.

Results

Here we report biomarker miRNA pairs effectively differentiating (i) patients with pulmonary system diseases (asthma, pneumonia and non-small cell lung cancer) and gastrointestinal (GI) system diseases (Crohn’s disease, stages I/II esophageal, gastric and colon cancers) from controls, each with 95% accuracy; (ii) patients with a pathology of the pulmonary system from patients with a pathology of the GI system with 94% accuracy; and (iii) cancer patients (stages I/II esophageal, gastric, colon cancers, or non-small cell lung cancer) from patients with inflammatory diseases (asthma, pneumonia, or Crohn’s disease) with 93%-95% accuracy.

Conclusions

The results obtained in the present study, along with the data reported by us and others previously, are encouraging and lay the ground for further investigation of the described approach for UST development.

Expression of herpes simplex virus 1 microRNAs in cell culture models of quiescent and latent infection

To facilitate studies of herpes simplex virus 1 latency, cell culture models of quiescent or latent infection have been developed. Using deep sequencing, we analyzed the expression of viral microRNAs (miRNAs) in two models employing human fibroblasts and one using rat neurons. In all cases, the expression patterns differed from that in productively infected cells, with the rat neuron pattern most closely resembling that found in latently infected human or mouse ganglia in vivo.

Individual and combined effects of DNA methylation and copy number alterations on miRNA expression in breast tumors

Background

The global effect of copy number and epigenetic alterations on miRNA expression in cancer is poorly understood. In the present study, we integrate genome-wide DNA methylation, copy number and miRNA expression and identify genetic mechanisms underlying miRNA dysregulation in breast cancer.

Results

We identify 70 miRNAs whose expression was associated with alterations in copy number or methylation, or both. Among these, five miRNA families are represented. Interestingly, the members of these families are encoded on different chromosomes and are complementarily altered by gain or hypomethylation across the patients. In an independent breast cancer cohort of 123 patients, 41 of the 70 miRNAs were confirmed with respect to aberration pattern and association to expression. In vitro functional experiments were performed in breast cancer cell lines with miRNA mimics to evaluate the phenotype of the replicated miRNAs. let-7e-3p, which in tumors is found associated with hypermethylation, is shown to induce apoptosis and reduce cell viability, and low let-7e-3p expression is associated with poorer prognosis. The overexpression of three other miRNAs associated with copy number gain, miR-21-3p, miR-148b-3p and miR-151a-5p, increases proliferation of breast cancer cell lines. In addition, miR-151a-5p enhances the levels of phosphorylated AKT protein.

Conclusions

Our data provide novel evidence of the mechanisms behind miRNA dysregulation in breast cancer. The study contributes to the understanding of how methylation and copy number alterations influence miRNA expression, emphasizing miRNA functionality through redundant encoding, and suggests novel miRNAs important in breast cancer.

A cross-platform comparison of affymetrix and Agilent microarrays reveals discordant miRNA expression in lung tumors of c-Raf transgenic mice

Non-coding RNAs play major roles in the translational control of gene expression. In order to identify disease-associated miRNAs in precursor lesions of lung cancer, RNA extracts from lungs of either c-Raf transgenic or wild-type (WT) mice were hybridized to the Agilent and Affymetrix miRNA microarray platforms, respectively. This resulted in the detection of a range of miRNAs varying between 111 and 267, depending on the presence or absence of the transgene, on the gender, and on the platform used. Importantly, when the two platforms were compared, only 11–16% of the 586 overlapping genes were commonly detected. With the Agilent microarray, seven miRNAs were identified as significantly regulated, of which three were selectively up-regulated in male transgenic mice. Much to our surprise, when the same samples were analyzed with the Affymetrix platform, only two miRNAs were identified as significantly regulated. Quantitative PCR performed with lung RNA extracts from WT and transgenic mice confirmed only partially the differential expression of significant regulated miRNAs and established that the Agilent platform failed to detect miR-433. Finally, bioinformatic analyses predicted a total of 152 mouse genes as targets of the regulated miRNAs of which 4 and 11 genes were significantly regulated at the mRNA level, respectively in laser micro-dissected lung dysplasia and lung adenocarcinomas of c-Raf transgenic mice. Furthermore, for many of the predicted mouse target genes expression of the coded protein was also repressed in human lung cancer when the publically available database of the Human Protein Atlas was analyzed, thus supporting the clinical significance of our findings. In conclusion, a significant difference in a cross-platform comparison was observed that will have important implications for research into miRNAs.

Determination of reliable reference genes for multi-generational gene expression analysis on C. elegans exposed to abused drug nicotine

RATIONAL

More research has recently been focused on multigenerational toxicogenomics impacts. Such studies rely on behavioral as well as genetic and epigenetic analyses using various biotechniques. Of these technologies, quantitative reverse transcriptase PCR is considered as a mature discovery and validation tool. Nevertheless, the interpretation of the resulting gene expression necessitates the establishment of reliable internal controls for normalization. No study has been performed to identify reliable reference genes in multigenerational settings.

OBJECTIVES

The primary aim was to evaluate the stability of 16 reference gene candidates in Caenorhabditis elegans exposed to nicotine and their two subsequent generations for determining the most reliable reference genes for multigenerational study.

METHODS

We exposed C. elegans to nicotine in the F0 generation and investigated the relative stabilities of 16 housekeeping genes in L4 larvae across three generations (F0, F1, and F2) using five statistical approaches (geNorm, ∆Ct method, NormFinder, BestKeeper, and RefFinder).

RESULTS

geNorm shows that CDC-42 and Y45F10D.4 were the most stable reference genes. Based on NormFinder, TBA-1, EIF3.C, ARP-6, CDC-42, and MDH2 may serve as the top reliable reference genes. Comparative ∆Ct method ranked TBA-1, CDC-42, EIF3.C, ARP-6, and Y45F10D.4 as the most stable reference genes. BestKeeper shows that Y45F10D.4, F35G12.2, TBA-1, CDC-42, and CSQ-1were better reference genes. Overall, TBA-1, CDC-42, EIF3.C, ARP-6, and Y45F10D.4 were the most reliable reference genes for mutigenerational nicotine-exposed study.

CONCLUSIONS

Of the 16 tested gene candidates, TBA-1 and CDC-42 were the two most stable reference genes for performing reliable gene expression normalization in the multigenerational impact of nicotine exposure.

Do transgenesis and marker-assisted backcross breeding produce substantially equivalent plants? A comparative study of transgenic and backcross rice carrying bacterial blight resistant gene Xa21

Background

The potential impact of genetically modified (GM) plants on human health has attracted much attention worldwide, and the issue remains controversial. This is in sharp contrast to the broad acceptance of plants produced by breeding through Marker Assisted Backcrossing (MAB).

Results

Focusing on transcriptome variation and perturbation to signaling pathways, we assessed the molecular and biological aspects of substantial equivalence, a general principle for food safety endorsed by the Food and Agricultural Organization and the World Health Organization, between a transgenic crop and a plant from MAB breeding. We compared a transgenic rice line (DXT) and a MAB rice line (DXB), both of which contain the gene Xa21 providing resistance to bacterial leaf blight. By using Next-Generation sequencing data of DXT, DXB and their parental line (D62B), we compared the transcriptome variation of DXT and DXB. Remarkably, DXT had 43% fewer differentially expressed genes (DEGs) than DXB. The genes exclusively expressed in DXT and in DXB have pathogen and stress defense functions. Functional categories of DEGs in DXT were comparable to that in DXB, and seven of the eleven pathways significantly affected by transgenesis were also perturbed by MAB breeding.

Conclusions

These results indicated that the transgenic rice and rice from MAB breeding are substantial equivalent at the transcriptome level, and paved a way for further study of transgenic rice, e.g., understanding the chemical and nutritional properties of the DEGs identified in the current study.

A study of small RNAs from cerebral neocortex of pathology-verified Alzheimer's disease, dementia with lewy bodies, hippocampal sclerosis, frontotemporal lobar dementia, and non-demented human controls

MicroRNAs (miRNAs) are small (20-22 nucleotides) regulatory non-coding RNAs that strongly influence gene expression. Most prior studies addressing the role of miRNAs in neurodegenerative diseases (NDs) have focused on individual diseases such as Alzheimer's disease (AD), making disease-to-disease comparisons impossible. Using RNA deep sequencing, we sought to analyze in detail the small RNAs (including miRNAs) in the temporal neocortex gray matter from non-demented controls (n = 2), AD (n = 5), dementia with Lewy bodies (n = 4), hippocampal sclerosis of aging (n = 4), and frontotemporal lobar dementia (FTLD) (n = 5) cases, together accounting for the most prevalent ND subtypes. All cases had short postmortem intervals, relatively high-quality RNA, and state-of-the-art neuropathological diagnoses. The resulting data (over 113 million reads in total, averaging 5.6 million reads per sample) and secondary expression analyses constitute an unprecedented look into the human cerebral cortical miRNome at a nucleotide resolution. While we find no apparent changes in isomiR or miRNA editing patterns in correlation with ND pathology, our results validate and extend previous miRNA profiling studies with regard to quantitative changes in NDs. In agreement with this idea, we provide independent cohort validation for changes in miR-132 expression levels in AD (n = 8) and FTLD (n = 14) cases when compared to controls (n = 8). The identification of common and ND-specific putative novel brain miRNAs and/or short-hairpin molecules is also presented. The challenge now is to better understand the impact of these and other alterations on neuronal gene expression networks and neuropathologies.

High-throughput sequencing reveals differential expression of miRNAs in intestine from sea cucumber during aestivation

The regulatory role of miRNA in gene expression is an emerging hot new topic in the control of hypometabolism. Sea cucumber aestivation is a complicated physiological process that includes obvious hypometabolism as evidenced by a decrease in the rates of oxygen consumption and ammonia nitrogen excretion, as well as a serious degeneration of the intestine into a very tiny filament. To determine whether miRNAs play regulatory roles in this process, the present study analyzed profiles of miRNA expression in the intestine of the sea cucumber (Apostichopus japonicus), using Solexa deep sequencing technology. We identified 308 sea cucumber miRNAs, including 18 novel miRNAs specific to sea cucumber. Animals sampled during deep aestivation (DA) after at least 15 days of continuous torpor, were compared with animals from a non-aestivation (NA) state (animals that had passed through aestivation and returned to the active state). We identified 42 differentially expressed miRNAs [RPM (reads per million) >10, |FC| (|fold change|) ≥ 1, FDR (false discovery rate) <0.01] during aestivation, which were validated by two other miRNA profiling methods: miRNA microarray and real-time PCR. Among the most prominent miRNA species, miR-200-3p, miR-2004, miR-2010, miR-22, miR-252a, miR-252a-3p and miR-92 were significantly over-expressed during deep aestivation compared with non-aestivation animals. Preliminary analyses of their putative target genes and GO analysis suggest that these miRNAs could play important roles in global transcriptional depression and cell differentiation during aestivation. High-throughput sequencing data and microarray data have been submitted to GEO database.

Evaluation of internal reference genes for quantitative expression analysis by real-time reverse transcription-PCR in somatic cells from goat milk

Reverse transcription (RT) quantitative real-time PCR (qPCR) is the most accurate and easy-to-perform technique to measure the expression level of a selected gene of interest by quantifying mRNA transcripts. The use of reference genes is commonly accepted as the most reliable approach to normalize RT-qPCR data and reduce possible errors generated in the quantification of gene expression. The optimal number and choice of reference genes are experimentally validated for specific tissues or cell types and experimental designs. To date, data on qPCR normalization in goats are scarce and the most suitable reference genes in this species have been identified for only a limited number of tissues. The aim of this study was to determine an optimal combination of stably expressed reference genes in caprine milk somatic cells (MSC) from healthy and infected mammary glands. For the purpose, we performed RT-qPCR for 10 commonly used reference genes from various functional classes and then determined their expression level in MSC from goats intramammary challenged with Staphylococcus aureus and in MSC from healthy controls, with a view to select genes whose stability would be unaffected under infection conditions. The geNorm and NormFinder algorithms were used for validating the reference genes. Furthermore, to demonstrate the importance of normalization of gene expression with appropriate reference genes, we tested the effect of using a combination of the least stable genes for expression analysis evaluation. On the basis of our evaluation, we recommend the use of a panel of reference genes that should include G6PD, YWHAZ, and ACTB for caprine MSC gene expression profiling. The expression of the 2 genes of interest, pentraxin-related protein (PTX3) and secreted phosphoprotein 1 (SPP1), was evaluated by RT-qPCR in all samples collected pre- and postinfection, and the recommended reference genes were used to normalize the data. Our study provides a validated panel of optimal reference genes for the identification of genes differentially expressed by qRT-PCR in caprine MSC. Moreover, we provided a set of intron-spanning primer sequences that could be suitable for gene expression experiments using SYBR Green chemistry on other caprine tissues and cells.

Performance comparison of digital microRNA profiling technologies applied on human breast cancer cell lines

MicroRNA profiling represents an important first-step in deducting individual RNA-based regulatory function in a cell, tissue, or at a specific developmental stage. Currently there are several different platforms to choose from in order to make the initial miRNA profiles. In this study we investigate recently developed digital microRNA high-throughput technologies. Four different platforms were compared including next generation SOLiD ligation sequencing and Illumina HiSeq sequencing, hybridization-based NanoString nCounter, and miRCURY locked nucleic acid RT-qPCR. For all four technologies, full microRNA profiles were generated from human cell lines that represent noninvasive and invasive tumorigenic breast cancer. This study reports the correlation between platforms, as well as a more extensive analysis of the accuracy and sensitivity of data generated when using different platforms and important consideration when verifying results by the use of additional technologies. We found all the platforms to be highly capable for microRNA analysis. Furthermore, the two NGS platforms and RT-qPCR all have equally high sensitivity, and the fold change accuracy is independent of individual miRNA concentration for NGS and RT-qPCR. Based on these findings we propose new guidelines and considerations when performing microRNA profiling.

Candidate microRNA biomarkers of pancreatic ductal adenocarcinoma: meta-analysis, experimental validation and clinical significance

BACKGROUND

The diagnostic and prognostic value of microRNA (miRNA) expression aberrations in pancreatic ductal adenocarcinoma (PDAC) has been studied extensively in recent years. However, differences in measurement platforms and lab protocols as well as small sample sizes can render gene expression levels incomparable.

METHODS

A comprehensive meta-review of published studies in PDAC that compared the miRNA expression profiles of PDAC tissues and paired neighbouring noncancerous pancreatic tissues was performed to determine candidate miRNA biomarkers for PDAC. Both a miRNA vote-counting strategy and a recently published Robust Rank Aggregation method were employed. In this review, a total of 538 tumour and 206 noncancerous control samples were included.

RESULTS

We identified a statistically significant miRNA meta-signature of seven up- and three down-regulated miRNAs. The experimental validation results showed that the miRNA expression levels were in accordance with the meta-signature. The results from the vote-counting strategy were consistent with those from the Robust Rank Aggregation method. The experimental validation confirmed that the statistically unique profiles identified by the meta-review approach could discriminate PDAC tissues from paired nonmalignant pancreatic tissues. In a cohort of 70 patients, the high expression of miR-21 (p=0.018, HR=2.610; 95% CI=1.179-5.777) and miR-31 (p=0.039, HR=2.735; 95% CI=1.317-6.426), the low expression of miR-375 (p=0.022, HR=2.337; 95% CI=1.431-5.066) were associated with poor overall survival following resection, independent of clinical covariates.

CONCLUSIONS

The identified miRNAs may be used to develop a panel of diagnostic and prognostic biomarkers for PDAC with sufficient sensitivity and specificity for use in a clinical setting.

Role of microRNAs and long-non-coding RNAs in CD4(+) T-cell differentiation

CD4(+) T lymphocytes orchestrate adaptive immune responses by differentiating into various subsets of effector T cells such as T-helper 1 (Th1), Th2, Th17, and regulatory T cells. These subsets have been generally described by master transcription factors that dictate the expression of cytokines and receptors, which ultimately define lymphocyte effector functions. However, the view of T-lymphocyte subsets as stable and terminally differentiated lineages has been challenged by increasing evidence of functional plasticity within CD4(+) T-cell subsets, which implies flexible programming of effector functions depending on time and space of T-cell activation. An outstanding question with broad basic and traslational implications relates to the mechanisms, besides transcriptional regulation, which define the plasticity of effector functions. In this study, we discuss the emerging role of regulatory non-coding RNAs in T-cell differentiation and plasticity. Not only microRNAs have been proven to be important for CD4(+) T-cell differentiation, but it is also likely that the overall T-cell functioning is the result of a multilayered network composed by coding RNAs as well as by short and long non-coding RNAs. The integrated study of all the nodes of this network will provide a comprehensive view of the molecular mechanisms underlying T-cell functions in health and disease.

Convergent transcription induces dynamic DNA methylation at disiRNA loci

Cytosine methylation of DNA is an important epigenetic gene silencing mechanism in plants, fungi, and animals. In the filamentous fungus Neurospora crassa, nearly all known DNA methylations occur in transposon relics and repetitive sequences, and DNA methylation does not depend on the canonical RNAi pathway. disiRNAs are Dicer-independent small non-coding RNAs that arise from gene-rich part of the Neurospora genome. Here we describe a new type of DNA methylation that is associated with the disiRNA loci. Unlike the known DNA methylation in Neurospora, disiRNA loci DNA methylation (DLDM) is highly dynamic and is regulated by an on/off mechanism. Some disiRNA production appears to rely on pol II directed transcription. Importantly, DLDM is triggered by convergent transcription and enriched in promoter regions. Together, our results establish a new mechanism that triggers DNA methylation.

DNA methylation in eukayrotes refers to the modification of cytidines at 5^th position with methyl group (5mC). Though absent in some species, DNA methylation is conserved across fungi, plants and animals and plays a critical role in X chromosome inactivation, genomic imprinting, transposon silencing etc. In addition, DNA methylation also occurs at the promoter sequence to regulate gene expression. Filamentous fungus Neurospora crassa has a well-known mechanism of DNA methylation for genomic defense. During sexual stage repetitive sequences (e.g. transposons) are recognized and point mutations are introduced. During vegetative stage these mutations serve as signals for establishing static DNA methylation to silence all copies of the sequences. In this study, we report a new type of DNA methylation in Neurospora. It is tightly linked to a type of non-coding small RNA termed dicer-independent siRNA (disiRNA) and therefore was termed disiRNA loci DNA methylation (DLDM). DLDM is dynamic regulated and shows an on/off pattern, i.e. most alleles contain no 5mC but some are densely methylated. Interestingly, DLDM can be triggered by convergent transcription and is accumulated at promoter regions. In summary, our findings demonstrate a new type of dynamic DNA methylation.

Direct quantification of mRNA and miRNA from cell lysates using reverse transcription real time PCR: a multidimensional analysis of the performance of reagents and workflows

Substantial efforts have been devoted to in vitro testing of candidate chemotherapeutics by profiling transcriptional changes across the collection of NCI-60 cell-lines. A work-flow with reagents that enable the direct quantification of RNA of different molecular sizes simultaneously in the same sample without laborious total RNA isolation will invariably increase the throughput and accuracy of the study. MicroRNAs (miRNAs) are known to regulate most cellular functions, acting post-transcriptionally by repressing numerous eukaryotic mRNAs. Recent findings on the remarkable stability of miRNA prompted us to investigate the feasibility of quantifying the expression levels of both mRNA and miRNA directly from cell lysates (cell-to-Ct). Multidimensional analyses of the expressions of mRNA and miRNA across seven NCI-60 cell lines and multiple reagents were conducted to assess the performances of these reagents and workflows for cell-to-Ct measurements using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Quantification of RNA species using lysates prepared from an in-house and one of the commercial reagents demonstrated comparable performance to those prepared by the more laborious and conventional method of using guanidinium-phenol-chloroform. Additionally, miRNA was found to be highly stable in the cell lysates when incubated at room temperature for prolonged period of time and subjected to multiple freeze-thaw cycles. In summary, this study demonstrated significant differences in pre-analytical performance of a variety of commercially available reagents and described a cost-effective reagent useful for rapid, scalable, and high-throughput workflow for the detection of mRNA and miRNA from the same biological sample.

Trametes versicolor protein YZP activates regulatory B lymphocytes - gene identification through de novo assembly and function analysis in a murine acute colitis model

Background

Trametes versicolor (Yun-Zhi) is a medicinal fungus used as a chemotherapy co-treatment to enhance anti-tumor immunity. Although the efficacies of T. versicolor extracts have been documented, the active ingredients and mechanisms underlying the actions of these extracts remain uncharacterized.

Results

We purified a new protein, YZP, from the fruiting bodies of T. versicolor and identified the gene encoding YZP using RNA-seq and de novo assembly technologies. YZP is a 12-kDa non-glycosylated protein comprising 139 amino acids, including an 18-amino acids signal peptide. YZP induced a greater than 60-fold increase in IL-10 secretion in mice B lymphocytes; moreover, YZP specifically triggered the differentiation of CD1d⁺ B cells into IL-10-producing regulatory B cells (Bregs) and enhanced the expression of CD1d. YZP-induced B cells suppressed approximately 40% of the LPS-activated macrophage production of inflammatory cytokines in a mixed leukocyte reaction and significantly alleviated the disease activity and colonic inflammation in a DSS-induced acute colitis murine model. Furthermore, YZP activated Breg function via interaction with TLR2 and TLR4 and up-regulation of the TLR-mediated signaling pathway.

Conclusions

We purified a novel Breg-stimulating protein, YZP, from T. versicolor and developed an advanced approach combining RNA-seq and de novo assembly technologies.to clone its gene. We demonstrated that YZP activated CD1d⁺ Breg differentiation through TLR2/4-mediated signaling pathway, and the YZP-stimulated B cells exhibited anti-inflammatory efficacies in vitro and in murine acute colitis models.

Universal digital high-resolution melt: a novel approach to broad-based profiling of heterogeneous biological samples

Comprehensive profiling of nucleic acids in genetically heterogeneous samples is important for clinical and basic research applications. Universal digital high-resolution melt (U-dHRM) is a new approach to broad-based PCR diagnostics and profiling technologies that can overcome issues of poor sensitivity due to contaminating nucleic acids and poor specificity due to primer or probe hybridization inaccuracies for single nucleotide variations. The U-dHRM approach uses broad-based primers or ligated adapter sequences to universally amplify all nucleic acid molecules in a heterogeneous sample, which have been partitioned, as in digital PCR. Extensive assay optimization enables direct sequence identification by algorithm-based matching of melt curve shape and Tm to a database of known sequence-specific melt curves. We show that single-molecule detection and single nucleotide sensitivity is possible. The feasibility and utility of U-dHRM is demonstrated through detection of bacteria associated with polymicrobial blood infection and microRNAs (miRNAs) associated with host response to infection. U-dHRM using broad-based 16S rRNA gene primers demonstrates universal single cell detection of bacterial pathogens, even in the presence of larger amounts of contaminating bacteria; U-dHRM using universally adapted Lethal-7 miRNAs in a heterogeneous mixture showcases the single copy sensitivity and single nucleotide specificity of this approach.

Profiling of circulating microRNAs after a bout of acute resistance exercise in humans

Recent studies have revealed a new aspect of physiological regulation in which microRNAs (miRNAs) play fundamental roles in diverse biological and pathological processes. Furthermore, it was recently discovered that miRNAs are stably secreted into blood and that circulating miRNAs may play important roles in cell-cell communication. Here, we examined whether the circulating miRNA profile is affected by acute resistance exercise. Twelve males performed a resistance exercise session (bench press and leg press), consisting of five sets of 10 repetitions at 70% of maximum strength, with a 1 min rest between sets. Blood samples were taken before exercise, and at 0 and 60 min, 1 day, and 3 days after exercise. The circulating miRNA profile was determined by microarray analysis. Quantitative real-time PCR confirmed that the miR-149* level increased three days after resistance exercise. In contrast, the miR-146a and miR-221 levels decreased three days after resistance exercise. Our findings suggest that circulating miRNA levels change in response to acute resistance exercise, and miRNAs may play important roles in resistance-exercise-induced adaptation.

Unraveling genomic variation from next generation sequencing data

Elucidating the content of a DNA sequence is critical to deeper understand and decode the genetic information for any biological system. As next generation sequencing (NGS) techniques have become cheaper and more advanced in throughput over time, great innovations and breakthrough conclusions have been generated in various biological areas. Few of these areas, which get shaped by the new technological advances, involve evolution of species, microbial mapping, population genetics, genome-wide association studies (GWAs), comparative genomics, variant analysis, gene expression, gene regulation, epigenetics and personalized medicine. While NGS techniques stand as key players in modern biological research, the analysis and the interpretation of the vast amount of data that gets produced is a not an easy or a trivial task and still remains a great challenge in the field of bioinformatics. Therefore, efficient tools to cope with information overload, tackle the high complexity and provide meaningful visualizations to make the knowledge extraction easier are essential. In this article, we briefly refer to the sequencing methodologies and the available equipment to serve these analyses and we describe the data formats of the files which get produced by them. We conclude with a thorough review of tools developed to efficiently store, analyze and visualize such data with emphasis in structural variation analysis and comparative genomics. We finally comment on their functionality, strengths and weaknesses and we discuss how future applications could further develop in this field.

Transcriptome-wide identification and characterization of microRNAs from castor bean (Ricinus communis L.)

BACKGROUND

MicroRNAs (miRNAs) are endogenously encoded small RNAs that post-transcriptionally regulate gene expression and play essential roles in numerous developmental and physiological processes. Currently, little information on the transcriptome and tissue-specific expression of miRNAs is available in the model non-edible oilseed crop castor bean (Ricinus communis L.), one of the most important non-edible oilseed crops cultivated worldwide. Recent advances in sequencing technologies have allowed the identification of conserved and novel miRNAs in many plant species. Here, we used high-throughput sequencing technologies to identify and characterize the miRNAs in castor bean.

RESULTS

Five small RNA libraries were constructed for deep sequencing from root tips, leaves, developing seeds (at the initial stage, seed1; and at the fast oil accumulation stage, seed2) and endosperms in castor bean. High-throughput sequencing generated a large number of sequence reads of small RNAs in this study. In total, 86 conserved miRNAs were identified, including 63 known and 23 newly identified. Sixteen miRNA isoform variants in length were found from the conserved miRNAs of castor bean. MiRNAs displayed diverse organ-specific expression levels among five libraries. Combined with criteria for miRNA annotation and a RT-PCR approach, 72 novel miRNAs and their potential precursors were annotated and 20 miRNAs newly identified were validated. In addition, new target candidates for miRNAs newly identified in this study were proposed.

CONCLUSIONS

The current study presents the first high-throughput small RNA sequencing study performed in castor bean to identify its miRNA population. It characterizes and increases the number of miRNAs and their isoforms identified in castor bean. The miRNA expression analysis provides a foundation for understanding castor bean miRNA organ-specific expression patterns. The present study offers an expanded picture of miRNAs for castor bean and other members in the family Euphorbiaceae.