Ren, Yingxue

Abstract Background: Genetic studies have shifted to sequencing-based rare variants discovery after decades of success in identifying common disease variants by Genome-Wide Association Studies using Single Nucleotide Polymorphism chips. Sequencing-based studies require large sample sizes for statistical power but often inadvertently introduce batch effects because samples are typically collected, processed, and sequenced at multiple centers. Conventionally, batch effects are first detected and visualized using Principal Components Analysis and then controlled by including batch covariates in the disease association models. For sequencing-based genetic studies, because all variants included in the association analyses have passed quality control measures, this conventional approach treats every variant as equal and ignores the substantial differences still remaining in variant qualities and characteristics such as genotype quality scores, alternative allele fractions (fraction of reads supporting alternative allele at a variant position) and sequencing depths. Results: In the Alzheimer’s Disease Sequencing Project (ADSP) exome dataset of 9,904 cases and controls, we discovered hidden variant-level differences between sample batches of three sequencing centers and two exome capture kits. Although sequencing centers were included as a covariate in our association models, we observed differences at the variant level in genotype quality and alternative allele fraction between samples processed by different exome capture kits that significantly impacted both the confidence of variant detection and the identification of disease-associated variants. Furthermore, we found that the association signals of a subset of top disease risk variants came exclusively from samples processed by one exome capture kit that was more effective at capturing the alternative alleles compared to the other kit. Conclusions: Our findings highlight the importance of additional variant-level quality control for large sequencing-based genetic studies. More importantly, we demonstrate that automatically filtering out variants with batch differences may lead to false negatives if the batch discordance came largely from quality differences and if the variants from one batch had better quality scores.

Abstract Genetic variants that define two distinct haplotypes at the TMEM106B locus have been implicated in multiple neurodegenerative diseases and in healthy brain ageing. In frontotemporal dementia (FTD), the high expressing TMEM106B risk haplotype was shown to increase susceptibility for FTD with TDP-43 inclusions (FTD-TDP) and to modify disease penetrance in progranulin mutation carriers (FTD-GRN). To elucidate the biological function of TMEM106B and determine whether lowering TMEM106B may be a viable therapeutic strategy, we performed brain transcriptomic analyses in 8-month-old animals from our recently developed Tmem106b−/− mouse model. We included 10 Tmem106b+/+ (wild-type), 10 Tmem106b+/− and 10 Tmem106−/− mice. The most differentially expressed genes (153 downregulated and 60 upregulated) were identified between Tmem106b−/− and wild-type animals, with an enrichment for genes implicated in myelination-related cellular processes including axon ensheathment and oligodendrocyte differentiation. Co-expression analysis also revealed that the most downregulated group of correlated genes was enriched for myelination-related processes. We further detected a significant loss of OLIG2-positive cells in the corpus callosum of Tmem106b−/− mice, which was present already in young animals (21 days) and persisted until old age (23 months), without worsening. Quantitative polymerase chain reaction revealed a reduction of differentiated but not undifferentiated oligodendrocytes cellular markers. While no obvious changes in myelin were observed at the ultrastructure levels in unchallenged animals, treatment with cuprizone revealed that Tmem106b−/− mice are more susceptible to cuprizone-induced demyelination and have a reduced capacity to remyelinate, a finding which we were able to replicate in a newly generated Tmem106b CRISPR/cas9 knock-out mouse model. Finally, using a TMEM106B HeLa knock-out cell line and primary cultured oligodendrocytes, we determined that loss of TMEM106B leads to abnormalities in the distribution of lysosomes and PLP1. Together these findings reveal an important function for TMEM106B in myelination with possible consequences for therapeutic strategies aimed at lowering TMEM106B levels.

APOE4 exacerbates α-synuclein pathology in a synucleinopathy mouse model, and APOE4 genotype associates with increased α-synuclein pathology in humans.

Carrying premature termination codons in 1 allele of the ABCA7 gene is associated with an increased risk for Alzheimer’s disease (AD). While the primary function of ABCA7 is to regulate the transport of phospholipids and cholesterol, ABCA7 is also involved in maintaining homeostasis of the immune system. Since inflammatory pathways causatively or consequently participate in AD pathogenesis, we studied the effects of Abca7 haplodeficiency in mice on brain immune responses under acute and chronic conditions. When acute inflammation was induced through peripheral lipopolysaccharide injection in control or heterozygous Abca7 knockout mice, partial ABCA7 deficiency diminished proinflammatory responses by impairing CD14 expression in the brain. On breeding to App NL-G-F knockin mice, we observed increased amyloid-β (Aβ) accumulation and abnormal endosomal morphology in microglia. Taken together, our results demonstrate that ABCA7 loss of function may contribute to AD pathogenesis by altering proper microglial responses to acute inflammatory challenges and during the development of amyloid pathology, providing insight into disease mechanisms and possible treatment strategies.

Objective: The ε4 allele of the APOE gene ( APOE4 ) is the strongest genetic risk factor for Alzheimer disease when compared with the common ε3 allele. Although there has been significant progress in understanding how apoE4 (apolipoprotein E4) drives amyloid pathology, its effects on amyloid-independent pathways, in particular cerebrovascular integrity and function, are less clear. Approach and Results: Here, we show that brain pericytes, the mural cells of the capillary walls, differentially modulate endothelial cell phenotype in an apoE isoform-dependent manner. Extracellular matrix protein induction, tube-like structure formation, and barrier formation were lower with endothelial cells cocultured with pericytes isolated from apoE4-targeted replacement (TR) mice compared with those from apoE3-TR mice. Importantly, aged apoE4-targeted replacement mice had decreased extracellular matrix protein expression and increased plasma protein leakages compared with apoE3-TR mice. Conclusions: ApoE4 impairs pericyte-mediated basement membrane formation, potentially contributing to the cerebrovascular effects of apoE4.

We previously reported an extremely rare case of follicular dendritic cell sarcoma (FDCS) presented as a thyroid mass. Given the rarity of this disease, there are no personalized and molecularly targeted treatment options due to the lack of knowledge in the genomic makeup of the tumor. A 44-year-old white woman was diagnosed with an extranodal FDCS in thyroid. The patient underwent a total thyroidectomy, central compartment dissection, parathyroid re-implantation, and adjuvant radiation therapy. Tumor DNA sequencing of 236 genes by FoundationOne panel found truncating mutations in PTEN and missense mutations in RET and TP53. However, patient-matched germline DNA was not sequenced which is critical for identification of true somatic mutations. Furthermore, the FoundationOne panel doesn't measure genomic rearrangements which have been shown to be abundant in sarcomas and are associated with sarcoma tumorigenesis and progression. In the current study, we carried out comprehensive genomic sequencing of the tumor, adjacent normal tissues, and patient-matched blood, in an effort to understand the genomic makeup of this rare extranodal FDCS and to identify potential therapeutic targets. Eighty-one somatic point mutations were identified in tumor but not in adjacent normal tissues or blood. A clonal truncating mutation in the CLTCL1 gene, which stabilizes the mitotic spindle, was likely a driver mutation of tumorigenesis and could explain the extensive copy number aberrations (CNAs) and genomic rearrangements in the tumor including a chr15/chr17 local chromothripsis resulted in 6 expressed fusion genes. The fusion gene HDGFRP3→SHC4 led to a 200-fold increase in the expression of oncogene SHC4 which is a potential target of the commercial drug Dasatinib. Missense mutations in ATM and splice-site mutation in VEGFR1 were also detected in addition to the TP53 missense mutation reported by FoundationOne.