Comprehensive multi-center assessment of small RNA-seq methods for quantitative miRNA profiling.

RNA-seq is increasingly used for quantitative profiling of small RNAs (for example, microRNAs, piRNAs and snoRNAs) in diverse sample types, including isolated cells, tissues and cell-free biofluids. The accuracy and reproducibility of the currently used small RNA-seq library preparation methods have not been systematically tested. Here we report results obtained by a consortium of nine labs that independently sequenced reference, 'ground truth' samples of synthetic small RNAs and human plasma-derived RNA. We assessed three commercially available library preparation methods that use adapters of defined sequence and six methods using adapters with degenerate bases. Both protocol- and sequence-specific biases were identified, including biases that reduced the ability of small RNA-seq to accurately measure adenosine-to-inosine editing in microRNAs. We found that these biases were mitigated by library preparation methods that incorporate adapters with degenerate bases. MicroRNA relative quantification between samples using small RNA-seq was accurate and reproducible across laboratories and methods.

Evaluation of commercially available small RNASeq library preparation kits using low input RNA.

BACKGROUND: Evolving interest in comprehensively profiling the full range of small RNAs present in small tissue biopsies and in circulating biofluids, and how the profile differs with disease, has launched small RNA sequencing (RNASeq) into more frequent use. However, known biases associated with small RNASeq, compounded by low RNA inputs, have been both a significant concern and a hurdle to widespread adoption. As RNASeq is becoming a viable choice for the discovery of small RNAs in low input samples and more labs are employing it, there should be benchmark datasets to test and evaluate the performance of new sequencing protocols and operators. In a recent publication from the National Institute of Standards and Technology, Pine et al., 2018, the investigators used a commercially available set of three tissues and tested performance across labs and platforms. RESULTS: In this paper, we further tested the performance of low RNA input in three commonly used and commercially available RNASeq library preparation kits; NEB Next, NEXTFlex, and TruSeq small RNA library preparation. We evaluated the performance of the kits at two different sites, using three different tissues (brain, liver, and placenta) with high (1 µg) and low RNA (10 ng) input from tissue samples, or 5.0, 3.0, 2.0, 1.0, 0.5, and 0.2 ml starting volumes of plasma. As there has been a lack of robust validation platforms for differentially expressed miRNAs, we also compared low input RNASeq data with their expression profiles on three different platforms (Abcam Fireplex, HTG EdgeSeq, and Qiagen miRNome). CONCLUSIONS: The concordance of RNASeq results on these three platforms was dependent on the RNA expression level; the higher the expression, the better the reproducibility. The results provide an extensive analysis of small RNASeq kit performance using low RNA input, and replication of these data on three downstream technologies.

Analysis of extracellular RNA in cerebrospinal fluid.

We examined the extracellular vesicle (EV) and RNA composition of pooled normal cerebrospinal fluid (CSF) samples and CSF from five major neurological disorders: Alzheimer's disease (AD), Parkinson's disease (PD), low-grade glioma (LGG), glioblastoma multiforme (GBM), and subarachnoid haemorrhage (SAH), representing neurodegenerative disease, cancer, and severe acute brain injury. We evaluated: (I) size and quantity of EVs by nanoparticle tracking analysis (NTA) and vesicle flow cytometry (VFC), (II) RNA yield and purity using four RNA isolation kits, (III) replication of RNA yields within and between laboratories, and (IV) composition of total and EV RNAs by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and RNA sequencing (RNASeq). The CSF contained ~106 EVs/µL by NTA and VFC. Brain tumour and SAH CSF contained more EVs and RNA relative to normal, AD, and PD. RT-qPCR and RNASeq identified disease-related populations of microRNAs and messenger RNAs (mRNAs) relative to normal CSF, in both total and EV fractions. This work presents relevant measures selected to inform the design of subsequent replicative CSF studies. The range of neurological diseases highlights variations in total and EV RNA content due to disease or collection site, revealing critical considerations guiding the selection of appropriate approaches and controls for CSF studies.

Total Extracellular Small RNA Profiles from Plasma, Saliva, and Urine of Healthy Subjects.

Interest in circulating RNAs for monitoring and diagnosing human health has grown significantly. There are few datasets describing baseline expression levels for total cell-free circulating RNA from healthy control subjects. In this study, total extracellular RNA (exRNA) was isolated and sequenced from 183 plasma samples, 204 urine samples and 46 saliva samples from 55 male college athletes ages 18-25 years. Many participants provided more than one sample, allowing us to investigate variability in an individual's exRNA expression levels over time. Here we provide a systematic analysis of small exRNAs present in each biofluid, as well as an analysis of exogenous RNAs. The small RNA profile of each biofluid is distinct. We find that a large number of RNA fragments in plasma (63%) and urine (54%) have sequences that are assigned to YRNA and tRNA fragments respectively. Surprisingly, while many miRNAs can be detected, there are few miRNAs that are consistently detected in all samples from a single biofluid, and profiles of miRNA are different for each biofluid. Not unexpectedly, saliva samples have high levels of exogenous sequence that can be traced to bacteria. These data significantly contribute to the current number of sequenced exRNA samples from normal healthy individuals.

A de novo missense mutation in ZMYND11 is associated with global developmental delay, seizures, and hypotonia.

Recently, mutations in the zinc finger MYND-type containing 11 (ZMYND11) gene were identified in patients with autism spectrum disorders, intellectual disability, aggression, and complex neuropsychiatric features, supporting that this gene is implicated in 10p15.3 microdeletion syndrome. We report a novel de novo variant in the ZMYND11 gene (p.Ser421Asn) in a patient with a complex neurodevelopmental phenotype. The patient is a 24-yr-old Caucasian/Filipino female with seizures, global developmental delay, sensorineural hearing loss, hypotonia, dysmorphic features, and other features including a happy disposition and ataxic gait similar to Angelman syndrome. In addition, this patient had uncommon features including eosinophilic esophagitis and multiple, severe allergies not described in similar ZMYND11 cases. This new case further supports the association of ZMYND11 with autistic-like phenotypes and suggests that ZMYND11 should be included in the list of potentially causative candidate genes in cases with complex neurodevelopmental phenotypes.

microRNA changes in liver tissue associated with fibrosis progression in patients with hepatitis C.

BACKGROUND & AIMS: Accumulating evidence indicates that microRNAs play a role in a number of disease processes including the pathogenesis of liver fibrosis in hepatitis C infection. Our goal is to add to the accruing information regarding microRNA deregulation in liver fibrosis to increase our understanding of the underlying mechanisms of pathology and progression. METHODS: We used next generation sequencing to profile all detectable microRNAs in liver tissue and serum from patients with hepatitis C, stages F1-F4 of fibrosis. RESULTS: We found altered expression of several microRNAs, in particular, miR-182, miR199a-5p, miR-200a-5p and miR-183 were found to be significantly upregulated in tissue from liver biopsies of hepatitis C patients with advanced fibrosis, stage F3 and F4, when compared with liver biopsies from patients with early fibrosis, stages F1 and F2. We also found miR-148-5p, miR-1260b, miR-122-3p and miR-378i among the microRNAs most significantly down-regulated from early to advanced fibrosis of the liver. We also sequenced the serum microRNAs; however, we were not able to detect significant changes in circulating microRNAs associated with fibrosis stage after adjusting for multiple tests. CONCLUSIONS: Adding measurements of tissue microRNAs acquired during routine biopsies will continue to increase our knowledge of underlying mechanisms of fibrosis. Our goal is that these data, in combination with studies from other researchers and future long-term studies, could be used to enhance the staging accuracy of liver biopsies and expand the surveillance of patients at increased risk for cancer and progression to advanced fibrosis.

Meeting report: discussions and preliminary findings on extracellular RNA measurement methods from laboratories in the NIH Extracellular RNA Communication Consortium.

Extracellular RNAs (exRNAs) have been identified in all tested biofluids and have been associated with a variety of extracellular vesicles, ribonucleoprotein complexes and lipoprotein complexes. Much of the interest in exRNAs lies in the fact that they may serve as signalling molecules between cells, their potential to serve as biomarkers for prediction and diagnosis of disease and the possibility that exRNAs or the extracellular particles that carry them might be used for therapeutic purposes. Among the most significant bottlenecks to progress in this field is the lack of robust and standardized methods for collection and processing of biofluids, separation of different types of exRNA-containing particles and isolation and analysis of exRNAs. The Sample and Assay Standards Working Group of the Extracellular RNA Communication Consortium is a group of laboratories funded by the U.S. National Institutes of Health to develop such methods. In our first joint endeavour, we held a series of conference calls and in-person meetings to survey the methods used among our members, placed them in the context of the current literature and used our findings to identify areas in which the identification of robust methodologies would promote rapid advancements in the exRNA field.

High-throughput sequencing reveals altered expression of hepatic microRNAs in nonalcoholic fatty liver disease-related fibrosis.

Recent evidence suggests that microRNAs (miRNAs), small, noncoding RNA molecules that regulate gene expression, may play a role in the regulation of metabolic disorders, including nonalcoholic fatty liver disease (NAFLD). To identify miRNAs that mediate NAFLD-related fibrosis, we used high-throughput sequencing to assess miRNAs obtained from liver biopsies of 15 individuals without NAFLD fibrosis (F0) and 15 individuals with severe NAFLD fibrosis or cirrhosis (F3-F4), matched for age, sex, body mass index, type 2 diabetes status, hemoglobin A1c, and use of diabetes medications. We used DESeq2 and Kruskal-Wallis test to identify miRNAs that were differentially expressed between NAFLD patients with or without fibrosis, adjusting for multiple testing using Bonferroni correction. We identified a total of 75 miRNAs showing statistically significant evidence (adjusted P value <0.05) for differential expression between the 2 groups, including 30 upregulated and 45 downregulated miRNAs. Quantitative reverse-transcription polymerase chain reaction analysis of selected miRNAs identified by sequencing validated 9 of 11 of the top differentially expressed miRNAs. We performed functional enrichment analysis of dysregulated miRNAs and identified several potential gene targets related to NAFLD-related fibrosis including hepatic fibrosis, hepatic stellate cell activation, transforming growth factor beta signaling, and apoptosis signaling. We identified forkhead box O3 and F-box WD repeat domain containing 7, E3 ubiquitin protein ligase (FBXW7) as potential targets of miR-182, and found that levels of forkhead box O3, but not FBXW7, were significantly decreased in fibrotic samples. These findings support a role for hepatic miRNAs in the pathogenesis of NAFLD-related fibrosis and yield possible new insight into the molecular mechanisms underlying the initiation and progression of liver fibrosis and cirrhosis.

Profiles of extracellular miRNA in cerebrospinal fluid and serum from patients with Alzheimer's and Parkinson's diseases correlate with disease status and features of pathology.

The discovery and reliable detection of markers for neurodegenerative diseases have been complicated by the inaccessibility of the diseased tissue--such as the inability to biopsy or test tissue from the central nervous system directly. RNAs originating from hard to access tissues, such as neurons within the brain and spinal cord, have the potential to get to the periphery where they can be detected non-invasively. The formation and extracellular release of microvesicles and RNA binding proteins have been found to carry RNA from cells of the central nervous system to the periphery and protect the RNA from degradation. Extracellular miRNAs detectable in peripheral circulation can provide information about cellular changes associated with human health and disease. In order to associate miRNA signals present in cell-free peripheral biofluids with neurodegenerative disease status of patients with Alzheimer's and Parkinson's diseases, we assessed the miRNA content in cerebrospinal fluid and serum from postmortem subjects with full neuropathology evaluations. We profiled the miRNA content from 69 patients with Alzheimer's disease, 67 with Parkinson's disease and 78 neurologically normal controls using next generation small RNA sequencing (NGS). We report the average abundance of each detected miRNA in cerebrospinal fluid and in serum and describe 13 novel miRNAs that were identified. We correlated changes in miRNA expression with aspects of disease severity such as Braak stage, dementia status, plaque and tangle densities, and the presence and severity of Lewy body pathology. Many of the differentially expressed miRNAs detected in peripheral cell-free cerebrospinal fluid and serum were previously reported in the literature to be deregulated in brain tissue from patients with neurodegenerative disease. These data indicate that extracellular miRNAs detectable in the cerebrospinal fluid and serum are reflective of cell-based changes in pathology and can be used to assess disease progression and therapeutic efficacy.

Comparison of Analysis Tools for miRNA High Throughput Sequencing Using Nerve Crush as a Model.

Recent advances in sample preparation and analysis for next generation sequencing have made it possible to profile and discover new miRNAs in a high throughput manner. In the case of neurological disease and injury, these types of experiments have been more limited. Possibly because tissues such as the brain and spinal cord are inaccessible for direct sampling in living patients, and indirect sampling of blood and cerebrospinal fluid are affected by low amounts of RNA. We used a mouse model to examine changes in miRNA expression in response to acute nerve crush. We assayed miRNA from both muscle tissue and blood plasma. We examined how the depth of coverage (the number of mapped reads) changed the number of detectable miRNAs in each sample type. We also found that samples with very low starting amounts of RNA (mouse plasma) made high depth of mature miRNA coverage more difficult to obtain. Each tissue must be assessed independently for the depth of coverage required to adequately power detection of differential expression, weighed against the cost of sequencing that sample to the adequate depth. We explored the changes in total mapped reads and differential expression results generated by three different software packages: miRDeep2, miRNAKey, and miRExpress and two different analysis packages, DESeq and EdgeR. We also examine the accuracy of using miRDeep2 to predict novel miRNAs and subsequently detect them in the samples using qRT-PCR.

Identification of extracellular miRNA in human cerebrospinal fluid by next-generation sequencing.

There has been a growing interest in using next-generation sequencing (NGS) to profile extracellular small RNAs from the blood and cerebrospinal fluid (CSF) of patients with neurological diseases, CNS tumors, or traumatic brain injury for biomarker discovery. Small sample volumes and samples with low RNA abundance create challenges for downstream small RNA sequencing assays. Plasma, serum, and CSF contain low amounts of total RNA, of which small RNAs make up a fraction. The purpose of this study was to maximize RNA isolation from RNA-limited samples and apply these methods to profile the miRNA in human CSF by small RNA deep sequencing. We systematically tested RNA isolation efficiency using ten commercially available kits and compared their performance on human plasma samples. We used RiboGreen to quantify total RNA yield and custom TaqMan assays to determine the efficiency of small RNA isolation for each of the kits. We significantly increased the recovery of small RNA by repeating the aqueous extraction during the phenol-chloroform purification in the top performing kits. We subsequently used the methods with the highest small RNA yield to purify RNA from CSF and serum samples from the same individual. We then prepared small RNA sequencing libraries using Illumina's TruSeq sample preparation kit and sequenced the samples on the HiSeq 2000. Not surprisingly, we found that the miRNA expression profile of CSF is substantially different from that of serum. To our knowledge, this is the first time that the small RNA fraction from CSF has been profiled using next-generation sequencing.