Impact of genome build on RNA-seq interpretation and diagnostics.

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network and Genomics Research to Elucidate the Genetics of Rare Disease Consortium. Across six routinely collected biospecimens, 61% of quantified genes were not influenced by genome build. However, we identified 1,492 genes with build-dependent quantification, 3,377 genes with build-exclusive expression, and 9,077 genes with annotation-specific expression across six routinely collected biospecimens, including 566 clinically relevant and 512 known OMIM genes. Further, we demonstrate that between builds for a given gene, a larger difference in quantification is well correlated with a larger change in expression outlier calling. Combined, we provide a database of genes impacted by build choice and recommend that transcriptomics-guided analyses and diagnoses are cross referenced with these data for robustness.

Integration of transcriptomics and long-read genomics prioritizes structural variants in rare disease.

Rare structural variants (SVs) - insertions, deletions, and complex rearrangements - can cause Mendelian disease, yet they remain difficult to accurately detect and interpret. We sequenced and analyzed Oxford Nanopore long-read genomes of 68 individuals from the Undiagnosed Disease Network (UDN) with no previously identified diagnostic mutations from short-read sequencing. Using our optimized SV detection pipelines and 571 control long-read genomes, we detected 716 long-read rare (MAF < 0.01) SV alleles per genome on average, achieving a 2.4x increase from short-reads. To characterize the functional effects of rare SVs, we assessed their relationship with gene expression from blood or fibroblasts from the same individuals, and found that rare SVs overlapping enhancers were enriched (LOR = 0.46) near expression outliers. We also evaluated tandem repeat expansions (TREs) and found 14 rare TREs per genome; notably these TREs were also enriched near overexpression outliers. To prioritize candidate functional SVs, we developed Watershed-SV, a probabilistic model that integrates expression data with SV-specific genomic annotations, which significantly outperforms baseline models that don't incorporate expression data. Watershed-SV identified a median of eight high-confidence functional SVs per UDN genome. Notably, this included compound heterozygous deletions in FAM177A1 shared by two siblings, which were likely causal for a rare neurodevelopmental disorder. Our observations demonstrate the promise of integrating long-read sequencing with gene expression towards improving the prioritization of functional SVs and TREs in rare disease patients.

Nanopore sequencing of 1000 Genomes Project samples to build a comprehensive catalog of human genetic variation.

Less than half of individuals with a suspected Mendelian condition receive a precise molecular diagnosis after comprehensive clinical genetic testing. Improvements in data quality and costs have heightened interest in using long-read sequencing (LRS) to streamline clinical genomic testing, but the absence of control datasets for variant filtering and prioritization has made tertiary analysis of LRS data challenging. To address this, the 1000 Genomes Project ONT Sequencing Consortium aims to generate LRS data from at least 800 of the 1000 Genomes Project samples. Our goal is to use LRS to identify a broader spectrum of variation so we may improve our understanding of normal patterns of human variation. Here, we present data from analysis of the first 100 samples, representing all 5 superpopulations and 19 subpopulations. These samples, sequenced to an average depth of coverage of 37x and sequence read N50 of 54 kbp, have high concordance with previous studies for identifying single nucleotide and indel variants outside of homopolymer regions. Using multiple structural variant (SV) callers, we identify an average of 24,543 high-confidence SVs per genome, including shared and private SVs likely to disrupt gene function as well as pathogenic expansions within disease-associated repeats that were not detected using short reads. Evaluation of methylation signatures revealed expected patterns at known imprinted loci, samples with skewed X-inactivation patterns, and novel differentially methylated regions. All raw sequencing data, processed data, and summary statistics are publicly available, providing a valuable resource for the clinical genetics community to discover pathogenic SVs.

APOE loss-of-function variants: Compatible with longevity and associated with resistance to Alzheimer's disease pathology.

The ε4 allele of apolipoprotein E (APOE) is the strongest genetic risk factor for sporadic Alzheimer's disease (AD). Knockdown of ε4 may provide a therapeutic strategy for AD, but the effect of APOE loss of function (LoF) on AD pathogenesis is unknown. We searched for APOE LoF variants in a large cohort of controls and patients with AD and identified seven heterozygote carriers of APOE LoF variants. Five carriers were controls (aged 71-90 years), one carrier was affected by progressive supranuclear palsy, and one carrier was affected by AD with an unremarkable age at onset of 75 years. Two APOE ε3/ε4 controls carried a stop-gain affecting ε4: one was cognitively normal at 90 years and had no neuritic plaques at autopsy; the other was cognitively healthy at 79 years, and lumbar puncture at 76 years showed normal levels of amyloid. These results suggest that ε4 drives AD risk through the gain of abnormal function and support ε4 knockdown as a viable therapeutic option.

Impact of genome build on RNA-seq interpretation and diagnostics.

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network (UDN) and Genomics Research to Elucidate the Genetics of Rare Disease (GREGoR) Consortium. We identified 2,800 genes with build-dependent quantification across six routinely-collected biospecimens, including 1,391 protein-coding genes and 341 known rare disease genes. We further observed multiple genes that only have detectable expression in a subset of genome builds. Finally, we characterized how genome build impacts the detection of outlier transcriptomic events. Combined, we provide a database of genes impacted by build choice, and recommend that transcriptomics-guided analyses and diagnoses are cross-referenced with these data for robustness.

Culture shock: microglial heterogeneity, activation, and disrupted single-cell microglial networks in vitro.

BACKGROUND: Microglia, the resident immune cells of the brain, play a critical role in numerous diseases, but are a minority cell type and difficult to genetically manipulate in vivo with viral vectors and other approaches. Primary cultures allow a more controlled setting to investigate these cells, but morphological and transcriptional changes upon removal from their normal brain environment raise many caveats from in vitro studies. METHODS: To investigate whether cultured microglia recapitulate in vivo microglial signatures, we used single-cell RNA sequencing (scRNAseq) to compare microglia freshly isolated from the brain to primary microglial cultures. We performed cell population discovery, differential expression analysis, and gene co-expression module analysis to compare signatures between in vitro and in vivo microglia. We constructed causal predictive network models of transcriptional regulators from the scRNAseq data and identified a set of potential key drivers of the cultured phenotype. To validate this network analysis, we knocked down two of these key drivers, C1qc and Prdx1, in primary cultured microglia and quantified changes in microglial activation markers. RESULTS: We found that, although often assumed to be a relatively homogenous population of cells in culture, in vitro microglia are a highly heterogeneous population consisting of distinct subpopulations of cells with transcriptional profiles reminiscent of macrophages and monocytes, and are marked by transcriptional programs active in neurodegeneration and other disease states. We found that microglia in vitro presented transcriptional activation of a set of "culture shock genes" not found in freshly isolated microglia, characterized by strong upregulation of disease-associated genes including Apoe, Lyz2, and Spp1, and downregulation of homeostatic microglial markers, including Cx3cr1, P2ry12, and Tmem119. Finally, we found that cultured microglia prominently alter their transcriptional machinery modulated by key drivers from the homeostatic to activated phenotype. Knockdown of one of these drivers, C1qc, resulted in downregulation of microglial activation genes Lpl, Lyz2, and Ccl4. CONCLUSIONS: Overall, our data suggest that when removed from their in vivo home environment, microglia suffer a severe case of "culture shock", drastically modulating their transcriptional regulatory network state from homeostatic to activated through upregulation of modules of culture-specific genes. Consequently, cultured microglia behave as a disparate cell type that does not recapitulate the homeostatic signatures of microglia in vivo. Finally, our predictive network model discovered potential key drivers that may convert activated microglia back to their homeostatic state, allowing for more accurate representation of in vivo states in culture. Knockdown of key driver C1qc partially attenuated microglial activation in vitro, despite C1qc being only weakly upregulated in culture. This suggests that even genes that are not strongly differentially expressed across treatments or preparations may drive downstream transcriptional changes in culture.

Accelerated identification of disease-causing variants with ultra-rapid nanopore genome sequencing.

Whole-genome sequencing (WGS) can identify variants that cause genetic disease, but the time required for sequencing and analysis has been a barrier to its use in acutely ill patients. In the present study, we develop an approach for ultra-rapid nanopore WGS that combines an optimized sample preparation protocol, distributing sequencing over 48 flow cells, near real-time base calling and alignment, accelerated variant calling and fast variant filtration for efficient manual review. Application to two example clinical cases identified a candidate variant in <8 h from sample preparation to variant identification. We show that this framework provides accurate variant calls and efficient prioritization, and accelerates diagnostic clinical genome sequencing twofold compared with previous approaches.

Widespread choroid plexus contamination in sampling and profiling of brain tissue.

The choroid plexus, a tissue responsible for producing cerebrospinal fluid, is found predominantly in the lateral and fourth ventricles of the brain. This highly vascularized and ciliated tissue is made up of specialized epithelial cells and capillary networks surrounded by connective tissue. Given the complex structure of the choroid plexus, this can potentially result in contamination during routine tissue dissection. Bulk and single-cell RNA sequencing studies, as well as genome-wide in situ hybridization experiments (Allen Brain Atlas), have identified several canonical markers of choroid plexus such as Ttr, Folr1, and Prlr. We used the Ttr gene as a marker to query the Gene Expression Omnibus database for transcriptome studies of brain tissue and identified at least some level of likely choroid contamination in numerous studies that could have potentially confounded data analysis and interpretation. We also analyzed transcriptomic datasets from human samples from Allen Brain Atlas and the Genotype-Tissue Expression (GTEx) database and found abundant choroid contamination, with regions in closer proximity to choroid more likely to be impacted such as hippocampus, cervical spinal cord, substantia nigra, hypothalamus, and amygdala. In addition, analysis of both the Allen Brain Atlas and GTEx datasets for differentially expressed genes between likely "high contamination" and "low contamination" groups revealed a clear enrichment of choroid plexus marker genes and gene ontology pathways characteristic of these ciliated choroid cells. Inclusion of these contaminated samples could result in biological misinterpretation or simply add to the statistical noise and mask true effects. We cannot assert that Ttr or other genes/proteins queried in targeted assays are artifacts from choroid contamination as some of these differentials may be due to true biological effects. However, for studies that have an unequal distribution of choroid contamination among groups, investigators may wish to remove contaminated samples from analyses or incorporate choroid marker gene expression into their statistical modeling. In addition, we suggest that a simple RT-qPCR or western blot for choroid markers would mitigate unintended choroid contamination for any experiment, but particularly for samples intended for more costly omic profiling. This study highlights an unexpected problem for neuroscientists, but it is also quite possible that unintended contamination of adjacent structures occurs during dissections for other tissues but has not been widely recognized.

Clusterin ameliorates tau pathology in vivo by inhibiting fibril formation.

The molecular chaperone Clusterin (CLU) impacts the amyloid pathway in Alzheimer's disease (AD) but its role in tau pathology is unknown. We observed CLU co-localization with tau aggregates in AD and primary tauopathies and CLU levels were upregulated in response to tau accumulation. To further elucidate the effect of CLU on tau pathology, we utilized a gene delivery approach in CLU knock-out (CLU KO) mice to drive expression of tau bearing the P301L mutation. We found that loss of CLU was associated with exacerbated tau pathology and anxiety-like behaviors in our mouse model of tauopathy. Additionally, we found that CLU dramatically inhibited tau fibrilization using an in vitro assay. Together, these results demonstrate that CLU plays a major role in both amyloid and tau pathologies in AD.

Systematic analysis of dark and camouflaged genes reveals disease-relevant genes hiding in plain sight.

BACKGROUND: The human genome contains "dark" gene regions that cannot be adequately assembled or aligned using standard short-read sequencing technologies, preventing researchers from identifying mutations within these gene regions that may be relevant to human disease. Here, we identify regions with few mappable reads that we call dark by depth, and others that have ambiguous alignment, called camouflaged. We assess how well long-read or linked-read technologies resolve these regions. RESULTS: Based on standard whole-genome Illumina sequencing data, we identify 36,794 dark regions in 6054 gene bodies from pathways important to human health, development, and reproduction. Of these gene bodies, 8.7% are completely dark and 35.2% are ≥ 5% dark. We identify dark regions that are present in protein-coding exons across 748 genes. Linked-read or long-read sequencing technologies from 10x Genomics, PacBio, and Oxford Nanopore Technologies reduce dark protein-coding regions to approximately 50.5%, 35.6%, and 9.6%, respectively. We present an algorithm to resolve most camouflaged regions and apply it to the Alzheimer's Disease Sequencing Project. We rescue a rare ten-nucleotide frameshift deletion in CR1, a top Alzheimer's disease gene, found in disease cases but not in controls. CONCLUSIONS: While we could not formally assess the association of the CR1 frameshift mutation with Alzheimer's disease due to insufficient sample-size, we believe it merits investigating in a larger cohort. There remain thousands of potentially important genomic regions overlooked by short-read sequencing that are largely resolved by long-read technologies.