Assessing characteristics of RNA amplification methods for single cell RNA sequencing.

BACKGROUND: Recently, measurement of RNA at single cell resolution has yielded surprising insights. Methods for single-cell RNA sequencing (scRNA-seq) have received considerable attention, but the broad reliability of single cell methods and the factors governing their performance are still poorly known. RESULTS: Here, we conducted a large-scale control experiment to assess the transfer function of three scRNA-seq methods and factors modulating the function. All three methods detected greater than 70% of the expected number of genes and had a 50% probability of detecting genes with abundance greater than 2 to 4 molecules. Despite the small number of molecules, sequencing depth significantly affected gene detection. While biases in detection and quantification were qualitatively similar across methods, the degree of bias differed, consistent with differences in molecular protocol. Measurement reliability increased with expression level for all methods and we conservatively estimate measurements to be quantitative at an expression level greater than ~5-10 molecules. CONCLUSIONS: Based on these extensive control studies, we propose that RNA-seq of single cells has come of age, yielding quantitative biological information.

Sporadic ALS has compartment-specific aberrant exon splicing and altered cell-matrix adhesion biology.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterized by progressive weakness from loss of motor neurons. The fundamental pathogenic mechanisms are unknown and recent evidence is implicating a significant role for abnormal exon splicing and RNA processing. Using new comprehensive genomic technologies, we studied exon splicing directly in 12 sporadic ALS and 10 control lumbar spinal cords acquired by a rapid autopsy system that processed nervous systems specifically for genomic studies. ALS patients had rostral onset and caudally advancing disease and abundant residual motor neurons in this region. We created two RNA pools, one from motor neurons collected by laser capture microdissection and one from the surrounding anterior horns. From each, we isolated RNA, amplified mRNA, profiled whole-genome exon splicing, and applied advanced bioinformatics. We employed rigorous quality control measures at all steps and validated findings by qPCR. In the motor neuron enriched mRNA pool, we found two distinct cohorts of mRNA signals, most of which were up-regulated: 148 differentially expressed genes (P

TDP-43 and ubiquitinated cytoplasmic aggregates in sporadic ALS are low frequency and widely distributed in the lower motor neuron columns independent of disease spread.

Ubiqitinated and TDP-43 immunoreactive cytoplasmic aggregates are hallmark features of ALS molecular pathology. Since clinically most ALS begins focally and advances contiguously, it is important to characterize their distribution. Our objective was to determine the extent and distribution of TDP-43 immunoreactive aggregates in the lower motor neuron columns as a function of disease onset, and to correlate ubiquitinated with TDP-43 aggregates in the lumbar region. We examined TDP-43 cytoplasmic aggregates at four separate neuraxis levels - hypoglossal nucleus and cervical, thoracic, and lumbar anterior horns - in five controls and 20 sporadic ALS nervous systems from patients whose disease began in various sites, i.e. five bulbar, five arm, five trunk, and five leg onsets. We correlated ubiquitinated to TDP-43 aggregates on adjacent histological sections for the lumbar regions. We found that TDP-43 cytoplasmic aggregates are seen in about 8% of motor neurons but there is marked variability between nervous systems, ranging from 0.4% to 20.6%. The aggregates are uniformly distributed within individual nervous systems. There is no obvious correlation between site of disease onset and rate of spread. Almost all ubiquitinated aggregates correlate to TDP-43 aggregates. Thus, TDP-43 immunoreactive cytoplasmic aggregates have a low overall average frequency that does not correlate with either disease course or clinical spread and is the prime ubiquitinated protein.

No evidence of HIV pol gene in spinal cord tissues in sporadic ALS by real-time RT-PCR.

An association of retroviruses and especially HIV has been reported in SALS, but CNS tissues have not been directly tested. To be valid, direct testing should be performed in light of the highly variable distribution of neuropathology from one patient to another and the degradable nature of RNA viruses. We aimed to seek direct evidence of retrovirus and especially HIV in SALS by real-time RT-PCR. We tested 20 SALS spinal cords in our repository specifically created for RNA preservation. We extracted RNA from anterior horns and surrounding tissue scrape from each sample and validated RNA quality. We indexed our search to neuropathologic degeneration and tested multiple areas in many of the nervous systems. We used real-time RT-PCR to a highly conserved region in the highly conserved pol gene and established the lowest limit of detection using a positive control. HIV-1 pol gene sequence was not detected in any of the 20 SALS nervous systems. Our assay was able to detect down to 60 viral copy numbers. In conclusion, we could not find direct evidence of retroviral/HIV activity in SALS CNS to support a relationship. However, our method could miss novel retroviruses with non-homology in the pol gene.

Robust RNA-Seq of aRNA-amplified single cell material collected by patch clamp.

Most single cell RNA sequencing protocols start with single cells dispersed from intact tissue. High-throughput processing of the separated cells is enabled using microfluidics platforms. However, dissociation of tissue results in loss of information about cell location and morphology and potentially alters the transcriptome. An alternative approach for collecting RNA from single cells is to re-purpose the electrophysiological technique of patch clamp recording. A hollow patch pipette is attached to individual cells, enabling the recording of electrical activity, after which the cytoplasm may be extracted for single cell RNA-Seq ("Patch-Seq"). Since the tissue is not disaggregated, the location of cells is readily determined, and the morphology of the cells is maintained, making possible the correlation of single cell transcriptomes with cell location, morphology and electrophysiology. Recent Patch-Seq studies utilizes PCR amplification to increase amount of nucleic acid material to the level required for current sequencing technologies. PCR is prone to create biased libraries - especially with the extremely high degrees of exponential amplification required for single cell amounts of RNA. We compared a PCR-based approach with linear amplifications and demonstrate that aRNA amplification (in vitro transcription, IVT) is more sensitive and robust for single cell RNA collected by a patch clamp pipette.

Gene Expression in Patient-Derived Neural Progenitors Implicates WNT5A Signaling in the Etiology of Schizophrenia.

BACKGROUND: Genome-wide association studies of schizophrenia have demonstrated that variations in noncoding regions are responsible for most of the common variation heritability of the disease. It is hypothesized that these risk variants alter gene expression. Therefore, studying alterations in gene expression in schizophrenia may provide a direct approach to understanding the etiology of the disease. In this study we use cultured neural progenitor cells derived from olfactory neuroepithelium (CNON cells) as a genetically unaltered cellular model to elucidate the neurodevelopmental aspects of schizophrenia. METHODS: We performed a gene expression study using RNA sequencing of CNON cells from 111 control subjects and 144 individuals with schizophrenia. Differentially expressed genes were identified with DESeq2 software, using covariates to correct for sex, age, library batches, and 1 surrogate variable component. RESULTS: A total of 80 genes were differentially expressed (false discovery rate < 10%), showing enrichment in cell migration, cell adhesion, developmental process, synapse assembly, cell proliferation, and related Gene Ontology categories. Cadherin and Wnt signaling pathways were positive in overrepresentation test, and, in addition, many genes were specifically involved in WNT5A signaling. The differentially expressed genes were modestly, but significantly, enriched in the genes overlapping single nucleotide polymorphisms with genome-wide significant association from the Psychiatric Genomics Consortium genome-wide association study of schizophrenia. We also found substantial overlap with genes associated with other psychiatric disorders or brain development, enrichment in the same Gene Ontology categories as genes with mutations de novo in schizophrenia, and studies of induced pluripotent stem cell-derived neural progenitor cells. CONCLUSIONS: CNON cells are a good model of the neurodevelopmental aspects of schizophrenia and can be used to elucidate the etiology of the disorder.

Dorsal Amygdala Neurotrophin-3 Decreases Anxious Temperament in Primates.

BACKGROUND: An early-life anxious temperament (AT) is a risk factor for the development of anxiety, depression, and comorbid substance abuse. We validated a nonhuman primate model of early-life AT and identified the dorsal amygdala as a core component of AT's neural circuit. Here, we combine RNA sequencing, viral-vector gene manipulation, functional brain imaging, and behavioral phenotyping to uncover AT's molecular substrates. METHODS: In response to potential threat, AT and brain metabolism were assessed in 46 young rhesus monkeys. We identified AT-related transcripts using RNA-sequencing data from dorsal amygdala tissue (including central nucleus of the amygdala [Ce] and dorsal regions of the basal nucleus). Based on the results, we overexpressed the neurotrophin-3 gene, NTF3, in the dorsal amygdala using intraoperative magnetic resonance imaging-guided surgery (n = 5 per group). RESULTS: This discovery-based approach identified AT-related alterations in the expression of well-established and novel genes, including an inverse association between NTRK3 expression and AT. NTRK3 is an interesting target because it is a relatively unexplored neurotrophic factor that modulates intracellular neuroplasticity pathways. Overexpression of the transcript for NTRK3's endogenous ligand, NTF3, in the dorsal amygdala resulted in reduced AT and altered function in AT's neural circuit. CONCLUSIONS: Together, these data implicate neurotrophin-3/NTRK3 signaling in the dorsal amygdala in mediating primate anxiety. More generally, this approach provides an important step toward understanding the molecular underpinnings of early-life AT and will be useful in guiding the development of treatments to prevent the development of stress-related psychopathology.