An extension to: Systematic assessment of commercially available low-input miRNA library preparation kits.

High-throughput sequencing has emerged as the favoured method to study microRNA (miRNA) expression, but biases introduced during library preparation have been reported. We recently compared the performance (sensitivity, reliability, titration response and differential expression) of six commercially-available kits on synthetic miRNAs and human RNA, where library preparation was performed by the vendors. We hereby supplement this study with data from two further commonly used kits (NEBNext, NEXTflex) whose manufacturers initially declined to participate. NEXTflex demonstrated the highest sensitivity, which may reflect its use of partially-randomized adapter sequences, but overall performance was lower than the QIAseq and TailorMix kits. NEBNext showed intermediate performance. We reaffirm that biases are kit specific, complicating the comparison of miRNA datasets generated using different kits.

Systematic assessment of commercially available low-input miRNA library preparation kits.

High-throughput sequencing is increasingly favoured to assay the presence and abundance of microRNAs (miRNAs) in biological samples, even from low RNA amounts, and a number of commercial vendors now offer kits that allow miRNA sequencing from sub-nanogram (ng) inputs. Although biases introduced during library preparation have been documented, the relative performance of current reagent kits has not been investigated in detail. Here, six commercial kits capable of handling <100ng total RNA input were used for library preparation, performed by kit manufactures, on synthetic miRNAs of known quantities and human total RNA samples. We compared the performance of miRNA detection sensitivity, reliability, titration response and the ability to detect differentially expressed miRNAs. In addition, we assessed the use of unique molecular identifiers (UMI) sequence tags in one kit. We observed differences in detection sensitivity and ability to identify differentially expressed miRNAs between the kits, but none were able to detect the full repertoire of synthetic miRNAs. The reliability within the replicates of all kits was good, while larger differences were observed between the kits, although none could accurately quantify the relative levels of the majority of miRNAs. UMI tags, at least within the input ranges tested, offered little advantage to improve data utility. In conclusion, biases in miRNA abundance are heavily influenced by the kit used for library preparation, suggesting that comparisons of datasets prepared by different procedures should be made with caution. This article is intended to assist researchers select the most appropriate kit for their experimental conditions.

A comprehensive profile of circulating RNAs in human serum.

Non-coding RNA (ncRNA) molecules have fundamental roles in cells and many are also stable in body fluids as extracellular RNAs. In this study, we used RNA sequencing (RNA-seq) to investigate the profile of small non-coding RNA (sncRNA) in human serum. We analyzed 10 billion Illumina reads from 477 serum samples, included in the Norwegian population-based Janus Serum Bank (JSB). We found that the core serum RNA repertoire includes 258 micro RNAs (miRNA), 441 piwi-interacting RNAs (piRNA), 411 transfer RNAs (tRNA), 24 small nucleolar RNAs (snoRNA), 125 small nuclear RNAs (snRNA) and 123 miscellaneous RNAs (misc-RNA). We also investigated biological and technical variation in expression, and the results suggest that many RNA molecules identified in serum contain signs of biological variation. They are therefore unlikely to be random degradation by-products. In addition, the presence of specific fragments of tRNA, snoRNA, Vault RNA and Y_RNA indicates protection from degradation. Our results suggest that many circulating RNAs in serum can be potential biomarkers.

Microbial Community Composition of Tap Water and Biofilms Treated with or without Copper-Silver Ionization.

Copper-silver ionization (CSI) is an in-house water disinfection method primarily installed to eradicate Legionella bacteria from drinking water distribution systems (DWDS). Its effect on the abundance of culturable Legionella and Legionella infections has been documented in several studies. However, the effect of CSI on other bacteria in DWDS is largely unknown. To investigate these effects, we characterized drinking water and biofilm communities in a hospital using CSI, in a neighboring building without CSI, and in treated drinking water at the local water treatment plant. We used 16S rDNA amplicon sequencing and Legionella culturing. The sequencing results revealed three distinct water groups: (1) cold-water samples (no CSI), (2) warm-water samples at the research institute (no CSI), and (3) warm-water samples at the hospital (after CSI; ANOSIM, p < 0.001). Differences between the biofilm communities exposed and not exposed to CSI were less clear (ANOSIM, p = 0.022). No Legionella were cultured, but limited numbers of Legionella sequences were recovered from all 25 water samples (0.2-1.4% relative abundance). The clustering pattern indicated local selection of Legionella types (Kruskal-Wallis, p < 0.001). Furthermore, one unclassified Betaproteobacteria OTU was highly enriched in CSI-treated warm water samples at the hospital (Kruskal-Wallis, p < 0.001).

Development and application of a whole genome amplicon sequencing method for infectious salmon anemia virus (ISAV).

Infectious salmon anemia (ISA) is an infectious disease primarily affecting farmed Atlantic salmon, Salmo salar, which is caused by the ISA virus (ISAV). ISAV belongs to the Orthomyxoviridae family. The disease is a serious condition resulting in reduced fish welfare and high mortality. In this study, we designed an amplicon-based sequencing protocol for whole genome sequencing of ISAV. The method consists of 80 ISAV-specific primers that cover 92% of the virus genome and was designed to be used on an Illumina MiSeq platform. The sequencing accuracy was investigated by comparing sequences with previously published Sanger sequences. The sequences obtained were nearly identical to those obtained by Sanger sequencing, thus demonstrating that sequences produced by this amplicon sequencing protocol had an acceptable accuracy. The amplicon-based sequencing method was used to obtain the whole genome sequence of 12 different ISAV isolates from a small local epidemic in the northern part of Norway. Analysis of the whole genome sequences revealed that segment reassortment took place between some of the isolates and could identify which segments that had been reassorted.

Multi-omics approach reveals dysregulated genes during hESCs neuronal differentiation exposure to paracetamol.

Prenatal paracetamol exposure has been associated with neurodevelopmental outcomes in childhood. Pharmacoepigenetic studies show differences in cord blood DNA methylation between unexposed and paracetamol-exposed neonates, however, causality and impact of long-term prenatal paracetamol exposure on brain development remain unclear. Using a multi-omics approach, we investigated the effects of paracetamol on an in vitro model of early human neurodevelopment. We exposed human embryonic stem cells undergoing neuronal differentiation with paracetamol concentrations corresponding to maternal therapeutic doses. Single-cell RNA-seq and ATAC-seq integration identified paracetamol-induced chromatin opening changes linked to gene expression. Differentially methylated and/or expressed genes were involved in neurotransmission and cell fate determination trajectories. Some genes involved in neuronal injury and development-specific pathways, such as KCNE3, overlapped with differentially methylated genes previously identified in cord blood associated with prenatal paracetamol exposure. Our data suggest that paracetamol may play a causal role in impaired neurodevelopment.

Robust neuronal differentiation of human embryonic stem cells for neurotoxicology.

Here, we describe a protocol for rapid neuronal differentiation from human embryonic stem cells (hESCs) toward a heterogenous population of telencephalic progenitors, immature and mature neurons, for drug-screening and early-brain differentiation studies. hESC neuronal differentiation depends on adhesion and minimal cell-passaging to avert monolayer cross-connectivity rupture. In this protocol, we detail optimized cell-seeding densities and coating conditions with high cell viability suitable for neurotoxicology and high-resolution single-cell omics studies. Daily media changes reduce compound instability and degradation for optimal screening. For complete details on the use and execution of this protocol, please refer to Samara et al. (2022).

A multi-omics approach to visualize early neuronal differentiation from hESCs in 4D.

Neuronal differentiation of pluripotent stem cells is an established method to study physiology, disease, and medication safety. However, the sequence of events in human neuronal differentiation and the ability of in vitro models to recapitulate early brain development are poorly understood. We developed a protocol optimized for the study of early human brain development and neuropharmacological applications. We comprehensively characterized gene expression and epigenetic profiles at four timepoints, because the cells differentiate from embryonic stem cells towards a heterogeneous population of progenitors, immature and mature neurons bearing telencephalic signatures. A multi-omics roadmap of neuronal differentiation, combined with searchable interactive gene analysis tools, allows for extensive exploration of early neuronal development and the effect of medications.

Exposure to a human relevant mixture of persistent organic pollutants or to perfluorooctane sulfonic acid alone dysregulates the developing cerebellum of chicken embryo.

Prenatal exposure to persistent organic pollutants (POPs) is associated with neurodevelopmental disorders. In the present study, we explored whether a human-relevant POP mixture affects the development of chicken embryo cerebellum. We used a defined mixture of 29 POPs, with chemical composition and concentrations based on blood levels in the Scandinavian population. We also evaluated exposure to a prominent compound in the mixture, perfluorooctane sulfonic acid (PFOS), alone. Embryos (n = 7-9 per exposure group) were exposed by injection directly into the allantois at embryonic day 13 (E13). Cerebella were isolated at E17 and subjected to morphological, RNA-seq and shot-gun proteomics analyses. There was a reduction in thickness of the molecular layer of cerebellar cortex in both exposure scenarios. Exposure to the POP mixture significantly affected expression of 65 of 13,800 transcripts, and 43 of 2,568 proteins, when compared to solvent control. PFOS alone affected expression of 80 of 13,859 transcripts, and 69 of 2,555 proteins. Twenty-five genes and 15 proteins were common for both exposure groups. These findings point to alterations in molecular events linked to retinoid X receptor (RXR) signalling, neuronal cell proliferation and migration, cellular stress responses including unfolded protein response, lipid metabolism, and myelination. Exposure to the POP mixture increased methionine oxidation, whereas PFOS decreased oxidation. Several of the altered genes and proteins are involved in a wide variety of neurological disorders. We conclude that POP exposure can interfere with fundamental aspects of neurodevelopment, altering molecular pathways that are associated with adverse neurocognitive and behavioural outcomes.

Ultralow amounts of DNA from long-term archived serum samples produce high-quality methylomes.

BACKGROUND: Long-term stored serum is considered challenging for epigenomic analyses: as there are no cells, circulating DNA is scarce, and amplification removes epigenetic signals. Additionally, pre-analytical treatments and storage might introduce biases and fragmentation to the DNA. In particular, starting with low-input DNA can result in low-diversity libraries. However, successful whole-genome bisulphite sequencing (WGBS) of such serum samples has the potential to open biobanks for epigenetic analyses and deliver novel prediagnostic biomarkers. Here, we perform WGBS using the Accel-NGS library preparation kit on ultralow amounts of DNA from long-term archived samples with diverse pretreatments from the Janus Serum Bank. RESULTS: Ninety-four of the 96 samples produced satisfactory methylation calls; an average of 578 M reads per sample generated a mean coverage of 17× and mean duplication level of 35%. Failed samples were related to poor bisulphite conversion rather than to sequencing or library preparation. We demonstrate the feasibility of WGBS on ultralow DNA yields from serum samples stored up to 48 years. CONCLUSIONS: Our results show the potential of large serum biobank collections for future epigenomic studies and biomarker discovery.

MicroRNA Expression Differences in Blood-Derived CD19+ B Cells of Methotrexate Treated Rheumatoid Arthritis Patients.

Rheumatoid arthritis (RA) is a complex disease with a wide range of underlying susceptibility factors. Recently, dysregulation of microRNAs (miRNAs) in RA have been reported in several immune cell types from blood. However, B cells have not been studied in detail yet. Given the autoimmune nature of RA with the presence of autoantibodies, CD19+ B cells are a key cell type in RA pathogenesis and alterations in CD19+ B cell subpopulations have been observed in patient blood. Therefore, we aimed to reveal the global miRNA repertoire and to analyze miRNA expression profile differences in homogenous RA patient phenotypes in blood-derived CD19+ B cells. Small RNA sequencing was performed on CD19+ B cells of newly diagnosed untreated RA patients (n=10), successfully methotrexate (MTX) treated RA patients in remission (MTX treated RA patients, n=18) and healthy controls (n=9). The majority of miRNAs was detected across all phenotypes. However, significant expression differences between MTX treated RA patients and controls were observed for 27 miRNAs, while no significant differences were seen between the newly diagnosed patients and controls. Several of the differentially expressed miRNAs were previously found to be dysregulated in RA including miR-223-3p, miR-486-3p and miR-23a-3p. MiRNA target enrichment analysis, using the differentially expressed miRNAs and miRNA-target interactions from miRTarBase as input, revealed enriched target genes known to play important roles in B cell activation, differentiation and B cell receptor signaling, such as STAT3, PRDM1 and PTEN. Interestingly, many of those genes showed a high degree of correlated expression in CD19+ B cells in contrast to other immune cell types. Our results suggest important regulatory functions of miRNAs in blood-derived CD19+ B cells of MTX treated RA patients and motivate for future studies investigating the interactive mechanisms between miRNA and gene targets, as well as the possible predictive power of miRNAs for RA treatment response.