RESUMEN
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.
Asunto(s)
Genoma Humano , RNA-Seq , Humanos , RNA-Seq/métodos , Genómica/métodos , Transcriptoma , Enfermedades Raras/genética , Enfermedades Raras/diagnóstico , Perfilación de la Expresión Génica/métodosRESUMEN
PURPOSE: The function of FAM177A1 and its relationship to human disease is largely unknown. Recent studies have demonstrated FAM177A1 to be a critical immune-associated gene. One previous case study has linked FAM177A1 to a neurodevelopmental disorder in 4 siblings. METHODS: We identified 5 individuals from 3 unrelated families with biallelic variants in FAM177A1. The physiological function of FAM177A1 was studied in a zebrafish model organism and human cell lines with loss-of-function variants similar to the affected cohort. RESULTS: These individuals share a characteristic phenotype defined by macrocephaly, global developmental delay, intellectual disability, seizures, behavioral abnormalities, hypotonia, and gait disturbance. We show that FAM177A1 localizes to the Golgi complex in mammalian and zebrafish cells. Intersection of the RNA sequencing and metabolomic data sets from FAM177A1-deficient human fibroblasts and whole zebrafish larvae demonstrated dysregulation of pathways associated with apoptosis, inflammation, and negative regulation of cell proliferation. CONCLUSION: Our data shed light on the emerging function of FAM177A1 and defines FAM177A1-related neurodevelopmental disorder as a new clinical entity.
RESUMEN
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.
Asunto(s)
Encéfalo , Plexo Coroideo , Biomarcadores/metabolismo , Encéfalo/metabolismo , Plexo Coroideo/metabolismo , Receptor 1 de Folato/metabolismo , Hipocampo/metabolismo , Humanos , Transcriptoma/genéticaAsunto(s)
Cuidados Críticos , Secuenciación de Nanoporos/métodos , Trastornos del Neurodesarrollo/diagnóstico , Adolescente , Preescolar , Femenino , Humanos , Lactante , Recién Nacido , Masculino , Persona de Mediana Edad , Mutación , Secuenciación de Nanoporos/economía , Trastornos del Neurodesarrollo/genética , Análisis de Secuencia de ADN/métodos , Estado Epiléptico/genéticaRESUMEN
CSF1R-related disorder (CSF1R-RD) is a neurodegenerative condition that predominantly affects white matter due to genetic alterations in the CSF1R gene, which is expressed by microglia. We studied an elderly man with a hereditary, progressive dementing disorder of unclear etiology. Standard genetic testing for leukodystrophy and other neurodegenerative conditions was negative. Brain autopsy revealed classic features of adult-onset leukoencephalopathy with axonal spheroids and pigmented glia (ALSP), including confluent white matter degeneration with axonal spheroids and pigmented glial cells in the affected white matter, consistent with CSF1R-RD. Subsequent long-read sequencing identified a novel deletion in CSF1R that was not detectable with short-read exome sequencing. To gain insight into potential mechanisms underlying white matter degeneration in CSF1R-RD, we studied multiple brain regions exhibiting varying degrees of white matter pathology. We found decreased CSF1R transcript and protein across brain regions, including intact white matter. Single nuclear RNA sequencing (snRNAseq) identified two disease-associated microglial cell states: lipid-laden microglia (expressing GPNMB, ATG7, LGALS1, LGALS3) and inflammatory microglia (expressing IL2RA, ATP2C1, FCGBP, VSIR, SESN3), along with a small population of CD44+ peripheral monocyte-derived macrophages exhibiting migratory and phagocytic signatures. GPNMB+ lipid-laden microglia with ameboid morphology represented the end-stage disease microglia state. Disease-associated oligodendrocytes exhibited cell stress signatures and dysregulated apoptosis-related genes. Disease-associated oligodendrocyte precursor cells (OPCs) displayed a failure in their differentiation into mature myelin-forming oligodendrocytes, as evidenced by upregulated LRP1, PDGFRA, SOX5, NFIA, and downregulated NKX2-2, NKX6.2, SOX4, SOX8, TCF7L2, YY1, ZNF488. Overall, our findings highlight microglia-oligodendroglia crosstalk in demyelination, with CSF1R dysfunction promoting phagocytic and inflammatory microglia states, an arrest in OPC differentiation, and oligodendrocyte depletion.
Asunto(s)
Neuroglía , Receptores de Factor Estimulante de Colonias de Granulocitos y Macrófagos , Humanos , Masculino , Receptores de Factor Estimulante de Colonias de Granulocitos y Macrófagos/genética , Receptores de Factor Estimulante de Colonias de Granulocitos y Macrófagos/metabolismo , Neuroglía/patología , Neuroglía/metabolismo , Leucoencefalopatías/genética , Leucoencefalopatías/patología , Leucoencefalopatías/metabolismo , Anciano , Microglía/patología , Microglía/metabolismo , Perfilación de la Expresión Génica , Transcriptoma , Sustancia Blanca/patología , Sustancia Blanca/metabolismo , Encéfalo/patología , Encéfalo/metabolismo , Receptor de Factor Estimulante de Colonias de MacrófagosRESUMEN
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.
RESUMEN
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.
RESUMEN
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.
Asunto(s)
Enfermedad de Alzheimer , Humanos , Alelos , Enfermedad de Alzheimer/genética , Enfermedad de Alzheimer/patología , Apolipoproteína E4/genética , Apolipoproteínas E/genética , Genotipo , Longevidad/genéticaRESUMEN
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.
RESUMEN
Reading is an essential skill that requires focused attention. However, much reading is done in non-optimal environments. These days, reading is often done on digital devices or with a digital device nearby. These devices often introduce momentary distractions during reading, interrupting with alerts, notifications, and pop-ups. In two eye-tracking experiments, we investigated how such momentary distractions affect reading. Participants read paragraphs while their eye movements were monitored. During half of the paragraphs, distractions appeared periodically on the screen that required a response from the participants. In Experiment 1, the distractions were arrows that the participant had to respond to and then could immediately forget. In Experiment 2, the participants performed a 1-back task that required them to remember the identity of the last distractor. Compared with the no-distraction condition, the respond-and-forget distractors of Experiment 1 had minimal impact on reading behaviour and comprehension, but the working-memory-load distractors of Experiment 2 led to increased rereading and decreased reading comprehension. It seems a simple pop-up does not disrupt reading, but a message you must remember will.
Asunto(s)
Comprensión , Lectura , Humanos , Comprensión/fisiología , Movimientos Oculares , Atención/fisiología , Recuerdo MentalRESUMEN
Background and Objectives: Single nucleotide variants near TMEM106B associate with risk of frontotemporal lobar dementia with TDP-43 inclusions (FTLD-TDP) and Alzheimer's disease (AD) in genome-wide association studies (GWAS), but the causal variant at this locus remains unclear. Here we asked whether a novel structural variant on TMEM106B is the causal variant. Methods: An exploratory analysis identified structural variants on neurodegeneration-related genes. Subsequent analyses focused on an Alu element insertion on the 3'UTR of TMEM106B. This study included data from longitudinal aging and neurogenerative disease cohorts at Stanford University, case-control cohorts in the Alzheimer's Disease Sequencing Project (ADSP), and expression and proteomics data from Washington University in St. Louis (WUSTL). 432 individuals from two Stanford aging cohorts were whole-genome long-read and short-read sequenced. 16,906 samples from ADSP were short-read sequenced. Genotypes, transcriptomics, and proteomics data were available in 1,979 participants from an aging and dementia cohort at WUSTL. Selection criteria were specific to each cohort. In primary analyses, the linkage disequilibrium between the TMEM106B locus variants in the FTLD-TDP GWAS and the 3'UTR insertion was estimated. We then estimated linkage by ancestry in the ADSP and evaluated the effect of the TMEM106B lead variant on mRNA and protein levels. Results: The primary analysis included 432 participants (52.5% females, age range 45-92 years old). We identified a 316 bp Alu insertion overlapping the TMEM106B 3'UTR tightly linked with top GWAS variants rs3173615(C) and rs1990622(A). In ADSP European-ancestry participants, this insertion is in equivalent linkage with rs1990622(A) (R2=0.962, D'=0.998) and rs3173615(C) (R2=0.960, D'=0.996). In African-ancestry participants, the insertion is in stronger linkage with rs1990622(A) (R2=0.992, D'=0.998) than with rs3173615(C) (R2=0.811, D'=0.994). In public datasets, rs1990622 was consistently associated with TMEM106B protein levels but not with mRNA expression. In the WUSTL dataset, rs1990622 is associated with TMEM106B protein levels in plasma and cerebrospinal fluid, but not with TMEM106B mRNA expression. Discussion: We identified a novel Alu element insertion in the 3'UTR of TMEM106B in tight linkage with the lead FTLD-TDP risk variant. The lead variant is associated with TMEM106B protein levels, but not expression. The 3'UTR insertion is a lead candidate for the causal variant at this complex locus, pending confirmation with functional studies.
RESUMEN
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.
Asunto(s)
Encéfalo , Microglía , Encéfalo/metabolismo , Regulación hacia Abajo , Homeostasis , Macrófagos , Microglía/metabolismoRESUMEN
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.
Asunto(s)
Secuenciación de Nanoporos , Nanoporos , Mapeo Cromosómico , Secuenciación de Nucleótidos de Alto Rendimiento/métodos , Humanos , Secuenciación Completa del Genoma/métodosRESUMEN
Bacterial antibiotic resistance is becoming a significant health threat, and rapid identification of antibiotic-resistant bacteria is essential to save lives and reduce the spread of antibiotic resistance. This paper analyzes the ability of machine learning algorithms (MLAs) to process data from a novel spectroscopic diagnostic device to identify antibiotic-resistant genes and bacterial species by comparison to available bacterial DNA sequences. Simulation results show that the algorithms attain from 92% accuracy (for genes) up to 99% accuracy (for species). This novel approach identifies genes and species by optically reading the percentage of A, C, G, T bases in 1000s of short 10-base DNA oligomers instead of relying on conventional DNA sequencing in which the sequence of bases in long oligomers provides genetic information. The identification algorithms are robust in the presence of simulated random genetic mutations and simulated random experimental errors. Thus, these algorithms can be used to identify bacterial species, to reveal antibiotic resistance genes, and to perform other genomic analyses. Some MLAs evaluated here are shown to be better than others at accurate gene identification and avoidance of false negative identification of antibiotic resistance.
RESUMEN
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.
Asunto(s)
Clusterina/genética , Clusterina/metabolismo , Agregación Patológica de Proteínas/genética , Tauopatías/genética , Proteínas tau/metabolismo , Anciano , Enfermedad de Alzheimer/genética , Enfermedad de Alzheimer/metabolismo , Enfermedad de Alzheimer/patología , Enfermedad de Alzheimer/fisiopatología , Animales , Ansiedad/fisiopatología , Humanos , Técnicas In Vitro , Ratones , Ratones Noqueados , Persona de Mediana Edad , Enfermedad de Pick/genética , Enfermedad de Pick/metabolismo , Enfermedad de Pick/patología , Enfermedad de Pick/fisiopatología , Agregación Patológica de Proteínas/metabolismo , Agregación Patológica de Proteínas/patología , Agregación Patológica de Proteínas/fisiopatología , Tauopatías/metabolismo , Tauopatías/patología , Tauopatías/fisiopatologíaRESUMEN
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.