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.

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.

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.

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.

A streamlined sequencing process enables identification of disease-causing variants in the clinic within 8 hours.

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.

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.

Genetic diversity of the attachment (H) protein gene of current field isolates of canine distemper virus

To characterize the variability of recent field isolates of canine distemper virus (CDV) from different hosts and geographical areas, we conducted nucleotide sequence analysis of the gene encoding the haemagglutinin (H), the attachment protein of this virus. Pronounced differences between field isolates were revealed in comparison to the Convac and Onderstepoort vaccine strains. The diversity of CDV appeared to exceed that determined for measles virus. Phylogenetic analysis also separated the field isolates of CDV from the vaccine strains and provided evidence for the existence of different contemporary genotypes of CDV. Isolates from a Greenlandic sledge dog and a Siberian seal formed a distinct lineage. The remaining isolates formed a group. This group contained two European isolates from mink and ferret, a single lineage comprising three European dog isolates, and another separate lineage of North American isolates from dog, javelina, raccoon and captive leopards.

Comparative analysis of the attachment protein gene (H) of dolphin morbillivirus

DMV, dolphin morbillivirus, a paramyxovirus of uncertain origin recently emerged in Mediterranean dolphins. This study presents the complete nucleotide sequence of the hemagglutinin (H) gene including the gene boundaries. The single open reading frame of the DMV H gene encodes a protein of 604 residues which exhibits overall sequence characteristics similar to the H genes of other morbilliviruses. When compared to its closest homologues, measles virus (MV) and rinderpest virus (RPV), DMV has, respectively, 44 and 46% of amino acid residues in identical positions. The primary sequence of the DMV H protein is markedly less conserved than that of the fusion protein. The comparative data at the genomic level correspond with cross-neutralization studies with different morbilliviruses. Retrospective serogical studies dating back to 1983 indicate DMV-like infections in whales of the eastern Atlantic. The presented data support and extend previous studies suggesting that this novel morbillivirus is one of the phylogenetically oldest morbilliviruses known to circulate today. The relationship of DMV and established morbilliviruses to the newly emerged candidate morbillivirus infecting horse and man is discussed.

Secretion of pancreastatins from the porcine digestive tract

Pancreastatin, a 49-amino acid peptide with a COOH-terminal glycine amide, was originally isolated from porcine pancreas, but pancreastatin immunoreactivity has been found in several neuroendocrine tissues. There are strong indications that pancreastatin is derived from chromogranin A, since the amino acid sequence 240-288 in porcine chromogranin A corresponds to pancreastatin flanked by typical signals for proteolytic processing. We have studied the effect of electric stimulation of the nervous supply to perfused porcine pancreas, antrum, nonantral stomach, and small intestine on the release of immunoreactive pancreastatin, and we characterized the molecular nature of the secreted immunoreactivity by using a radioimmunoassay specific for the COOH-terminal glycine amide of porcine pancreastatin in combination with chromatography. In all tissues nerve stimulation significantly increased the release of immunoreactive pancreastatin. The secreted immunoreactive pancreastatin was heterogeneous, consisting of pancreastatin itself, a COOH-terminal pancreastatin fragment, and NH2-terminally extended pancreastatin forms. Pancreastatin predominated in the perfusate from pancreas and antrum, whereas mainly NH2-terminally extended molecular forms were secreted from the antrectomized stomach and small intestine. The different molecular forms of pancreastatin were secreted from the perfused organs in the same molar ratio as they occur in extracts of the corresponding tissues. Thus, pancreastatin and other chromogranin A-derived peptides in organ-specific proportions regularly accompany the secretion of the peptide hormones from the gastrointestinal tissues on appropriate stimulation.

Secretion of pancreastatins from isolated, perfused, porcine adrenal glands

Pancreastatin is a 49 amino acid peptide with a C-terminal glycine amide originally isolated from porcine pancreas. There are strong indications that pancreastatin is derived from chromogranin A, since the amino acid sequence 240-288 in porcine chromogranin A contains pancreastatin flanked by typical signals for proteolytic processing. Several molecular forms of immunoreactive pancreastatin have been reported to be present in porcine adrenal medulla, but no information on the secretion of the peptides is available. We studied stimuli for the release of immunoreactive pancreastatin from adrenal glands as well as the molecular nature of the released immunoreactivity using isolated, perfused, porcine adrenal glands with intact splanchnic nerve supply and a radioimmunoassay specific for the C-terminal glycine amide of porcine pancreastatin in combination with chromatography. Stimulation of the splanchnic nerves greatly enhanced the release of immunoreactive pancreastatin by a mechanism that seems to involve cholinergic nicotinic transmission. The pancreastatin immunoreactivity was due to large amounts of a N-terminally extended pancreastatin form, possibly corresponding to the chromogranin A(1-288) fragment and small amounts of pancreastatin and a C-terminal pancreastatin fragment. Adrenal extracts also contained unprocessed chromogranin A. We conclude that in the porcine adrenals at least 25% of chromogranin A is processed to smaller molecular forms with the pancreastatin sequence forming their C-terminus. These forms appear to be released in parallel with the catecholamines.

Binding of C-reactive protein to Aspergillus fumigatus fractions

Calcium-dependent binding of C-reactive protein (CRP) to Aspergillus fumigatus was determined by enzyme-linked immunosorbent assay. A homogenate of young hyphae was fractionated by hydrophobic interaction chromatography followed by gel filtration. High CRP-binding activity was found in a fraction of mol. wt c. 500,000 which was characterised by strong binding to the hydrophobic column. Three fractions of less conspicuous CRP-binding activity were identified (c. 500 000, 150 000 and 150 000-50 000 mol. wt respectively). In these four fractions, phosphorylcholine was detected by an anti-phosphorylcholine mouse hybridoma antibody. Some CRP-binding activity in fractions with low affinity for the hydrophobic column did not correspond closely with the presence of phosphorylcholine. It is suggested that C-reactive substance in A. fumigatus is heterogeneous. The C-reactive substances did not correspond with fractions containing major antigens (470 000 and 250 000 mol. wt respectively) which elicit a strong immune response in man.