Genome-wide interaction analysis of pathological hallmarks in Alzheimer's disease

Novel early-onset Alzheimer-associated genes influence risk through dysregulation of glutamate, immune activation, and intracell signaling pathways

Alzheimer Disease (AD) is a highly polygenic disease that presents with relatively earlier onset (≤70yo; EOAD) in about 5% of cases. Around 90% of these EOAD cases remain unexplained by pathogenic mutations. Using data from EOAD cases and controls, we performed a genome-wide association study (GWAS) and trans-ancestry meta-analysis on non-Hispanic Whites (NHW, NCase=6,282, NControl=13,386), African Americans (AA NCase=782, NControl=3,663) and East Asians (NCase=375, NControl=838 CO). We identified eight novel significant loci: six in the ancestry-specific analyses and two in the trans-ancestry analysis. By integrating gene-based analysis, eQTL, pQTL and functional annotations, we nominate four novel genes that are involved in microglia activation, glutamate production, and signaling pathways. These results indicate that EOAD, although sharing many genes with LOAD, harbors unique genes and pathways that could be used to create better prediction models or target identification for this type of AD.

A Genome-Wide Interaction Study of Erythrocyte ω-3 Polyunsaturated Fatty Acid Species and Memory in the Framingham Heart Study Offspring Cohort

BACKGROUND

Cognitive decline, and more specifically Alzheimer's disease, continues to increase in prevalence globally, with few, if any, adequate preventative approaches. Several tests of cognition are utilized in the diagnosis of cognitive decline that assess executive function, short- and long-term memory, cognitive flexibility, and speech and motor control. Recent studies have separately investigated the genetic component of both cognitive health, using these measures, and circulating fatty acids.

OBJECTIVES

We aimed to examine the potential moderating effect of main species of ω-3 polyunsaturated fatty acids (PUFAs) on an individual's genetically conferred risk of cognitive decline.

METHODS

The Offspring cohort from the Framingham Heart Study was cross-sectionally analyzed in this genome-wide interaction study (GWIS). Our sample included all individuals with red blood cell ω-3 PUFA, genetic, cognitive testing (via Trail Making Tests [TMTs]), and covariate data (N = 1620). We used linear mixed effects models to predict each of the 3 cognitive measures (TMT A, TMT B, and TMT D) by each ω-3 PUFA, single nucleotide polymorphism (SNP) (0, 1, or 2 minor alleles), ω-3 PUFA by SNP interaction term, and adjusting for sex, age, education, APOE ε4 genotype status, and kinship (relatedness).

RESULTS

Our analysis identified 31 unique SNPs from 24 genes reaching an exploratory significance threshold of 1×10-5. Fourteen of the 24 genes have been previously associated with the brain/cognition, and 5 genes have been previously associated with circulating lipids. Importantly, 8 of the genes we identified, DAB1, SORCS2, SERINC5, OSBPL3, CPA6, DLG2, MUC19, and RGMA, have been associated with both cognition and circulating lipids. We identified 22 unique SNPs for which individuals with the minor alleles benefit substantially from increased ω-3 fatty acid concentrations and 9 unique SNPs for which the common homozygote benefits.

CONCLUSIONS

In this GWIS of ω-3 PUFA species on cognitive outcomes, we identified 8 unique genes with plausible biology suggesting individuals with specific polymorphisms may have greater potential to benefit from increased ω-3 PUFA intake. Additional replication in prospective settings with more diverse samples is needed.

Associations of Heavy Metals with Cognitive Function: An Epigenome-Wide View of DNA Methylation and Mediation Analysis

OBJECTIVE

Exposure to heavy metals has been reported to be associated with impaired cognitive function, but the underlying mechanisms remain unclear. This pilot study aimed to identify key heavy metal elements associated with cognitive function and further explore the potential mediating role of metal-related DNA methylation.

METHODS

Blood levels of arsenic, cadmium, lead, copper, manganese, and zinc and genome-wide DNA methylations were separately detected in peripheral blood in 155 older adults. Cognitive function was evaluated using the Mini-Mental State Examination (MMSE). Least absolute shrinkage and selection operator penalized regression and Bayesian kernel machine regression were used to identify metals associated with cognitive function. An epigenome-wide association study examined the DNA methylation profile of the identified metal, and mediation analysis investigated its mediating role.

RESULTS

The MMSE scores showed a significant decrease of 1.61 (95% confidence interval [CI]: -2.64, -0.59) with each 1 standard deviation increase in ln-transformed arsenic level; this association was significant in multiple-metal models and dominated the overall negative effect of 6 heavy metal mixture on cognitive function. Seventy-three differentially methylated positions were associated with blood arsenic (p < 1.0 × 10-5). The methylation levels at cg05226051 (annotated to TDRD3) and cg18886932 (annotated to GAL3ST3) mediated 24.8% and 25.5% of the association between blood arsenic and cognitive function, respectively (all p < 0.05).

INTERPRETATION

Blood arsenic levels displayed a negative association with the cognitive function of older adults. This finding shows that arsenic-related DNA methylation alterations are critical partial mediators that may serve as potential biomarkers for further mechanism-related studies. ANN NEUROL 2024.

Loci on chromosome 20 interact with rs16969968 to influence cigarettes per day in European ancestry individuals

BACKGROUND

The understanding of the molecular genetic contributions to smoking is largely limited to the additive effects of individual single nucleotide polymorphisms (SNPs), but the underlying genetic risk is likely to also include dominance, epistatic, and gene-environment interactions.

METHODS

To begin to address this complexity, we attempted to identify genetic interactions between rs16969968, the most replicated SNP associated with smoking quantity, and all SNPs and genes across the genome.

RESULTS

Using the UK Biobank European subsample, we found one SNP, rs1892967, and two genes, PCNA and TMEM230, that showed a significant genome-wide interaction with rs16969968 for log10 CPD and raw CPD, respectively, in a sample of 116 442 individuals who self-reported currently or previously smoking. We extended these analyses to individuals of South Asian descent and meta-analyzed the combined sample of 117 212 individuals of European and South Asian ancestry. We replicated the gene findings in a meta-analysis of five Finnish samples (N=40 140): FinHealth, FINRISK, Finnish Twin Cohort, GeneRISK, and Health-2000-2011.

CONCLUSIONS

To our knowledge, this represents the first reliable epistatic association between single nucleotide polymorphisms for smoking behaviors and provides a novel direction for possible future functional studies related to this interaction. Furthermore, this work demonstrates the feasibility of these analyses by pooling multiple datasets across various ancestries, which may be applied to other top SNPs for smoking and/or other phenotypes.

Genetic factors associated with suicidal behaviors and alcohol use disorders in an American Indian population

American Indians (AI) demonstrate the highest rates of both suicidal behaviors (SB) and alcohol use disorders (AUD) among all ethnic groups in the US. Rates of suicide and AUD vary substantially between tribal groups and across different geographical regions, underscoring a need to delineate more specific risk and resilience factors. Using data from over 740 AI living within eight contiguous reservations, we assessed genetic risk factors for SB by investigating: (1) possible genetic overlap with AUD, and (2) impacts of rare and low-frequency genomic variants. Suicidal behaviors included lifetime history of suicidal thoughts and acts, including verified suicide deaths, scored using a ranking variable for the SB phenotype (range 0-4). We identified five loci significantly associated with SB and AUD, two of which are intergenic and three intronic on genes AACSP1, ANK1, and FBXO11. Nonsynonymous rare and low-frequency mutations in four genes including SERPINF1 (PEDF), ZNF30, CD34, and SLC5A9, and non-intronic rare and low-frequency mutations in genes OPRD1, HSD17B3 and one lincRNA were significantly associated with SB. One identified pathway related to hypoxia-inducible factor (HIF) regulation, whose 83 nonsynonymous rare and low-frequency variants on 10 genes were significantly linked to SB as well. Four additional genes, and two pathways related to vasopressin-regulated water metabolism and cellular hexose transport, also were strongly associated with SB. This study represents the first investigation of genetic factors for SB in an American Indian population that has high risk for suicide. Our study suggests that bivariate association analysis between comorbid disorders can increase statistical power; and rare and low-frequency variant analysis in a high-risk population enabled by whole-genome sequencing has the potential to identify novel genetic factors. Although such findings may be population specific, rare functional mutations relating to PEDF and HIF regulation align with past reports and suggest a biological mechanism for suicide risk and a potential therapeutic target for intervention.

Proteome-wide association studies using summary proteomic data identified 23 risk genes of Alzheimer's disease

Characterizing the genetic mechanisms underlying Alzheimer's disease (AD) dementia is crucial for developing new therapeutics. Proteome-wide association study (PWAS) integrating proteomics data with genome-wide association study (GWAS) summary data was shown as a powerful tool for detecting risk genes. The identified PWAS risk genes can be interpretated as having genetic effects mediated through the genetically regulated protein abundances. Existing PWAS analyses of AD often rely on the availability of individual-level proteomics and genetics data of a reference cohort. Leveraging summary-level protein quantitative trait loci (pQTL) reference data of multiple relevant tissues is expected to improve PWAS findings for studying AD. Here, we applied our recently developed OTTERS tool to conduct PWAS of AD dementia, by leveraging summary-level pQTL data of brain, cerebrospinal fluid (CSF), and plasma tissues, and multiple statistical methods. For each target protein, imputation models of the protein abundance with genetic predictors were trained from summary-level pQTL data, estimating a set of pQTL weights for considered genetic predictors. PWAS p-values were obtained by integrating GWAS summary data of AD dementia with estimated pQTL weights. PWAS p-values from multiple statistical methods were combined by the aggregated Cauchy association test to yield one omnibus PWAS p-value for the target protein. We identified significant PWAS risk genes through omnibus PWAS p-values and analyzed their protein-protein interactions using STRING. Their potential causal effects were assessed by the probabilistic Mendelian randomization (PMR-Egger). As a result, we identified a total of 23 significant PWAS risk genes for AD dementia in brain, CSF, and plasma tissues, including 7 novel findings. We showed that 15 of these risk genes were interconnected within a protein-protein interaction network involving the well-known AD risk gene of APOE and 5 novel findings, and enriched in immune functions and lipids pathways including positive regulation of immune system process, positive regulation of macrophage proliferation, humoral immune response, and high-density lipoprotein particle clearance. Existing biological evidence was found to relate our novel findings with AD. We validated the mediated causal effects of 14 risk genes (60.8%). In conclusion, we identified both known and novel PWAS risk genes, providing novel insights into the genetic mechanisms in brain, CSF, and plasma tissues, and targeted therapeutics development of AD dementia. Our study also demonstrated the effectiveness of integrating public available summary-level pQTL data with GWAS summary data for mapping risk genes of complex human diseases.

PDE4 Gene Family Variants Are Associated with Response to Apremilast Treatment in Psoriasis

Moderate-to-severe psoriasis (Ps) treatment includes systemic drugs and biological agents. Apremilast, a small molecule primarily metabolized by cytochrome CYP3A4, modulates the immune system by specifically inhibiting phosphodiesterase type 4 (PDE4) isoforms and is currently used for the treatment of Ps and psoriatic arthritis (PsA). Clinical trials and real-world data showed variable efficacy in response among Ps patients underlying the need for personalized therapy. This study implements a candidate-gene and a network-based approach to identify genetic markers associated with apremilast response in forty-nine Greek Ps patients. Our data revealed an association of sixty-four SNPs within or near PDE4 and CYP3A4 genes, four SNPs in ncRNAs ANRIL, LINC00941 and miR4706, which influence the abundance or function of PDE4s, and thirty-three SNPs within fourteen genes whose protein products either interact directly with PDE4 proteins or constitute components of the cAMP signaling pathway which is modulated by PDE4s. Notably, fifty-six of the aforementioned SNPs constitute eQTLs for the respective genes in relevant to psoriasis tissues/cells implying that these variants could be causal. Our analysis provides a number of novel genetic variants that, upon validation in larger cohorts, could be utilized as predictive markers regarding the response of Ps patients to apremilast treatment.

Characterization of brain resilience in Alzheimer's disease using polygenic risk scores and further improvement by integrating mitochondria-associated loci

INTRODUCTION

Identification of high-risk people for Alzheimer's disease (AD) is critical for prognosis and early management. Longitudinal epidemiologic studies have observed heterogeneity in the brain and cognitive aging. Brain resilience was described as above-expected cognitive function. The "resilience" framework has been shown to correlate with individual characteristics such as genetic factors and age. Besides, accumulative evidence has confirmed the association of mitochondria with the pathogenesis of AD. However, it is challenging to assess resilience through genetic metrics, in particular incorporating mitochondria-associated loci.

OBJECTIVES

In this paper, we first demonstrated that polygenic risk scores (PRS) could characterize individuals' resilience levels. Then, we indicated that mitochondria-associated loci could improve the performance of PRSs, providing more reliable measurements for the prevention and diagnosis of AD.

METHODS

The discovery (N = 1,550) and independent validation samples (N = 2,090) were used to construct nine types of PRSs containing mitochondria-related loci (PRSMT) from both biological and statistical aspects and combined them with known AD risk loci derived from genome-wide association studies (GWAS).Individuals' levels of brain resilience were comprehensively measured by linear regression models using eight pathological characteristics.

RESULTS

It was found that PRSs could characterize brain resilience levels (e.g., Pearson correlation test Pmin = 7.96×10-9). Moreover, the performance of PRS models could be efficiently improved by incorporating a small number of mitochondria-related loci (e.g., Pearson correlation test P improved from 1.41×10-3 to 6.09×10-6). PRSs' ability to characterize brain resilience was validated. More importantly, by incorporating some mitochondria-related loci, the performance of PRSs in measuring brain resilience could be significantly improved.

CONCLUSION

Our findings imply that mitochondria may play an important role in brain resilience, and targeting mitochondria may open a new door to AD prevention and therapy.

Serum proteomics reveals APOE dependent and independent protein signatures in Alzheimer's disease

The current demand for early intervention, prevention, and treatment of late onset Alzheimer's disease (LOAD) warrants deeper understanding of the underlying molecular processes which could contribute to biomarker and drug target discovery. Utilizing high-throughput proteomic measurements in serum from a prospective population-based cohort of older adults (n = 5,294), we identified 303 unique proteins associated with incident LOAD (median follow-up 12.8 years). Over 40% of these proteins were associated with LOAD independently of APOE-ε4 carrier status. These proteins were implicated in neuronal processes and overlapped with protein signatures of LOAD in brain and cerebrospinal fluid. We found 17 proteins which LOAD-association was strongly dependent on APOE-ε4 carrier status. Most of them showed consistent associations with LOAD in cerebrospinal fluid and a third had brain-specific gene expression. Remarkably, four proteins in this group (TBCA, ARL2, S100A13 and IRF6) were downregulated by APOE-ε4 yet upregulated as a consequence of LOAD as determined in a bi-directional Mendelian randomization analysis, reflecting a potential response to the disease onset. Accordingly, the direct association of these proteins to LOAD was reversed upon APOE-ε4 genotype adjustment, a finding which we replicate in an external cohort (n = 719). Our findings provide an insight into the dysregulated pathways that may lead to the development and early detection of LOAD, including those both independent and dependent on APOE-ε4. Importantly, many of the LOAD-associated proteins we find in the circulation have been found to be expressed - and have a direct link with AD - in brain tissue. Thus, the proteins identified here, and their upstream modulating pathways, provide a new source of circulating biomarker and therapeutic target candidates for LOAD.

DNA methylation at birth and lateral ventricular volume in childhood: a neuroimaging epigenetics study

BACKGROUND

Lateral ventricular volume (LVV) enlargement has been repeatedly linked to schizophrenia; yet, what biological factors shape LVV during early development remain unclear. DNA methylation (DNAm), an essential process for neurodevelopment that is altered in schizophrenia, is a key molecular system of interest.

METHODS

In this study, we conducted the first epigenome-wide association study of neonatal DNAm in cord blood with LVV in childhood (measured using T1-weighted brain scans at 10 years), based on data from a large population-based birth cohort, the Generation R Study (N = 840). Employing both probe-level and methylation profile score (MPS) approaches, we further examined whether epigenetic modifications identified at birth in cord blood are: (a) also observed cross-sectionally in childhood using peripheral blood DNAm at age of 10 years (Generation R, N = 370) and (b) prospectively associated with LVV measured in young adulthood in an all-male sample from the Avon Longitudinal Study of Parents and Children (ALSPAC, N = 114).

RESULTS

At birth, DNAm levels at four CpGs (annotated to potassium channel tetramerization domain containing 3, KCTD3; SHH signaling and ciliogenesis regulator, SDCCAG8; glutaredoxin, GLRX) prospectively associated with childhood LVV after genome-wide correction; these genes have been implicated in brain development and psychiatric traits including schizophrenia. An MPS capturing a broader epigenetic profile of LVV - but not individual top hits - showed significant cross-sectional associations with LVV in childhood in Generation R and prospectively associated with LVV in early adulthood within ALSPAC.

CONCLUSIONS

This study finds suggestive evidence that DNAm at birth prospectively associates with LVV at different life stages, albeit with small effect sizes. The prediction of MPS on LVV in a childhood sample and an independent male adult sample further underscores the stability and reproducibility of DNAm as a potential marker for LVV. Future studies with larger samples and comparable time points across development are needed to further elucidate how DNAm associates with this clinically relevant brain structure and risk for neuropsychiatric disorders, and what factors explain the identified DNAm profile of LVV at birth.

Genetic association study of Preterm birth and Gestational age in a population-based case-control study in Peru

Preterm birth (PTB) is an adverse pregnancy outcome affecting ∼15 million pregnancies worldwide. Genetic studies have identified several candidate loci for PTB, but results remain inconclusive and limited to European populations. Thus, we conducted a genome-wide association study (GWAS) of PTB and gestational age at delivery (GA) among 2,212 Peruvian women. PTB cases delivered ≥ 20 weeks' but < 37 weeks' gestation, while controls delivered at term (≥ 37 weeks but < 42 weeks). After imputation (TOPMED) and quality control, we assessed the association of ∼6 million SNPs with PTB and GA using multivariable regression models adjusted for maternal age and the first two genetic principal components. In silico functional analysis (FUMA-GWAS) was conducted among top signals detected with an arbitrary P < 1.0×10 -5 in each GWAS. We sought to replicate genetic associations with PTB and GA identified in Europeans, and we developed a genetic risk score for GA based on European markers. Mean GA was 30 ± 4 weeks in PTB cases (N=933) and 39 ± 1 in the controls (N=1,279). PTB cases were slightly older and had higher C-sections and vaginal bleeding than controls. No association was identified at genome-wide level. Top suggestive ( P < 1.0×10 -5 ) signals were seen at rs13151645 ( LINC01182 ) for PTB, and at rs72824565 ( CTNNA2 ) for GA. Top PTB variants were enriched for biological pathways associated with polyketide, progesterone, steroid hormones, and glycosyl metabolism. Top GA variants were enriched in intronic regions and cancer pathways, and these genes were upregulated in the brain and subcutaneous adipose tissue. In combination with non-genetic risk factors, top SNPs explained 14% and 15% of the phenotypic variance of PTB and GA in our sample, but these results need to be interpreted with caution. Variants in WNT4 associated with GA in Europeans were replicated in our study. The genetic risk score based in European markers, was associated with a 2-day longer GA (R 2 =0.003, P =0.002) per standard deviation increase in the score in our sample. This genetic association study identified various signals suggestively associated with PTB and GA in a non- European population; they were linked to relevant biological pathways related to the metabolism of progesterone, prostanoid, and steroid hormones, and genes associated with GA were significantly upregulated in relevant tissues for the pathophysiology of PTB based on the in- silico functional analysis. None of these top variants overlapped with signals previously identified for PTB or GA in Europeans.

Serum proteomics reveals APOE dependent and independent protein signatures in Alzheimer's disease

The current demand for early intervention, prevention, and treatment of late onset Alzheimer's disease (LOAD) warrants deeper understanding of the underlying molecular processes which could contribute to biomarker and drug target discovery. Utilizing high-throughput proteomic measurements in serum from a prospective population-based cohort of older adults (n=5,294), we identified 303 unique proteins associated with incident LOAD (median follow-up 12.8 years). Over 40% of these proteins were associated with LOAD independently of APOE-ε4 carrier status. These proteins were implicated in neuronal processes and overlapped with protein signatures of LOAD in brain and cerebrospinal fluid. We found 17 proteins which LOAD-association was strongly dependent on APOE-ε4 carrier status. Most of them showed consistent associations with LOAD in cerebrospinal fluid and a third had brain-specific gene expression. Remarkably, four proteins in this group (TBCA, ARL2, S100A13 and IRF6) were downregulated by APOE-ε4 yet upregulated as a consequence of LOAD as determined in a bi-directional Mendelian randomization analysis, reflecting a potential response to the disease onset. Accordingly, the direct association of these proteins to LOAD was reversed upon APOE-ε4 genotype adjustment, a finding which we replicate in an external cohort (n=719). Our findings provide an insight into the dysregulated pathways that may lead to the development and early detection of LOAD, including those both independent and dependent on APOE-ε4. Importantly, many of the LOAD-associated proteins we find in the circulation have been found to be expressed - and have a direct link with AD - in brain tissue. Thus, the proteins identified here, and their upstream modulating pathways, provide a new source of circulating biomarker and therapeutic target candidates for LOAD.

Structural Variation Detection and Association Analysis of Whole-Genome-Sequence Data from 16,905 Alzheimer's Diseases Sequencing Project Subjects

Structural variations (SVs) are important contributors to the genetics of human diseases. However, their role in Alzheimer's disease (AD) remains largely unstudied due to challenges in accurately detecting SVs. We analyzed whole-genome sequencing data from the Alzheimer's Disease Sequencing Project (N = 16,905) and identified 400,234 (168,223 high-quality) SVs. Laboratory validation yielded a sensitivity of 82% (85% for high-quality). We found a significant burden of deletions and duplications in AD cases, particularly for singletons and homozygous events. On AD genes, we observed the ultra-rare SVs associated with the disease, including protein-altering SVs in ABCA7, APP, PLCG2, and SORL1. Twenty-one SVs are in linkage disequilibrium (LD) with known AD-risk variants, exemplified by a 5k deletion in complete LD with rs143080277 in NCK2. We also identified 16 SVs associated with AD and 13 SVs linked to AD-related pathological/cognitive endophenotypes. This study highlights the pivotal role of SVs in shaping our understanding of AD genetics.

DPYSL2/CRMP2 isoform B knockout in human iPSC-derived glutamatergic neurons confirms its role in mTOR signaling and neurodevelopmental disorders

The DPYSL2/CRMP2 gene encodes a microtubule-stabilizing protein crucial for neurogenesis and is associated with numerous psychiatric and neurodegenerative disorders including schizophrenia, bipolar disorder, and Alzheimer's disease. DPYSL2 generates multiple RNA and protein isoforms, but few studies have differentiated between them. We previously reported an association of a functional variant in the DPYSL2-B isoform with schizophrenia (SCZ) and demonstrated in HEK293 cells that this variant reduced the length of cellular projections and created transcriptomic changes that captured schizophrenia etiology by disrupting mTOR signaling-mediated regulation. In the present study, we follow up on these results by creating, to our knowledge, the first models of endogenous DPYSL2-B knockout in human induced pluripotent stem cells (iPSCs) and neurons. CRISPR/Cas9-faciliated knockout of DPYSL2-B in iPSCs followed by Ngn2-induced differentiation to glutamatergic neurons showed a reduction in DPYSL2-B/CRMP2-B RNA and protein with no observable impact on DPYSL2-A/CRMP2-A. The average length of dendrites in knockout neurons was reduced up to 58% compared to controls. Transcriptome analysis revealed disruptions in pathways highly relevant to psychiatric disease including mTOR signaling, cytoskeletal dynamics, immune function, calcium signaling, and cholesterol biosynthesis. We also observed a significant enrichment of the differentially expressed genes in SCZ-associated loci from genome-wide association studies (GWAS). Our findings expand our previous results to neuronal cells, clarify the functions of the human DPYSL2-B isoform and confirm its involvement in molecular pathologies shared between many psychiatric diseases.

Disease-specific selective vulnerability and neuroimmune pathways in dementia revealed by single cell genomics

The development of successful therapeutics for dementias requires an understanding of their shared and distinct molecular features in the human brain. We performed single-nuclear RNAseq and ATACseq in Alzheimer disease (AD), Frontotemporal degeneration (FTD), and Progressive Supranuclear Palsy (PSP), analyzing 40 participants, yielding over 1.4M cells from three brain regions ranging in vulnerability and pathological burden. We identify 35 shared disease-associated cell types and 14 that are disease-specific, replicating those previously identified in AD. Disease - specific cell states represent molecular features of disease-specific glial-immune mechanisms and neuronal vulnerability in each disorder, layer 4/5 intra-telencephalic neurons in AD, layer 2/3 intra-telencephalic neurons in FTD, and layer 5/6 near-projection neurons in PSP. We infer intrinsic disease-associated gene regulatory networks, which we empirically validate by chromatin footprinting. We find that causal genetic risk acts in specific neuronal and glial cells that differ across disorders, primarily non-neuronal cells in AD and specific neuronal subtypes in FTD and PSP. These data illustrate the heterogeneous spectrum of glial and neuronal composition and gene expression alterations in different dementias and identify new therapeutic targets by revealing shared and disease-specific cell states.

Genome-wide association studies and cross-population meta-analyses investigating short and long sleep duration

Sleep duration has been linked to a wide range of negative health outcomes and to reduced life expectancy. We present genome-wide association studies of short ( ≤ 5 h) and long ( ≥ 10 h) sleep duration in adults of European (N = 445,966), African (N = 27,785), East Asian (N = 3141), and admixed-American (N = 16,250) ancestry from UK Biobank and the Million Veteran Programme. In a cross-population meta-analysis, we identify 84 independent loci for short sleep and 1 for long sleep. We estimate SNP-based heritability for both sleep traits in each ancestry based on population derived linkage disequilibrium (LD) scores using cov-LDSC. We identify positive genetic correlation between short and long sleep traits (rg = 0.16 ± 0.04; p = 0.0002), as well as similar patterns of genetic correlation with other psychiatric and cardiometabolic phenotypes. Mendelian randomisation reveals a directional causal relationship between short sleep and depression, and a bidirectional causal relationship between long sleep and depression.

Structural Variation Detection and Association Analysis of Whole-Genome-Sequence Data from 16,905 Alzheimer's Diseases Sequencing Project Subjects

Structural variations (SVs) are important contributors to the genetics of numerous human diseases. However, their role in Alzheimer's disease (AD) remains largely unstudied due to challenges in accurately detecting SVs. Here, we analyzed whole-genome sequencing data from the Alzheimer's Disease Sequencing Project (ADSP, N=16,905 subjects) and identified 400,234 (168,223 high-quality) SVs. We found a significant burden of deletions and duplications in AD cases (OR=1.05, P=0.03), particularly for singletons (OR=1.12, P=0.0002) and homozygous events (OR=1.10, P<0.0004). On AD genes, the ultra-rare SVs, including protein-altering SVs in ABCA7, APP, PLCG2, and SORL1, were associated with AD (SKAT-O P=0.004). Twenty-one SVs are in linkage disequilibrium (LD) with known AD-risk variants, e.g., a deletion (chr2:105731359-105736864) in complete LD (R2=0.99) with rs143080277 (chr2:105749599) in NCK2. We also identified 16 SVs associated with AD and 13 SVs associated with AD-related pathological/cognitive endophenotypes. Our findings demonstrate the broad impact of SVs on AD genetics.

Integrative Bayesian tensor regression for imaging genetics applications

Identifying biomarkers for Alzheimer's disease with a goal of early detection is a fundamental problem in clinical research. Both medical imaging and genetics have contributed informative biomarkers in literature. To further improve the performance, recently, there is an increasing interest in developing analytic approaches that combine data across modalities such as imaging and genetics. However, there are limited methods in literature that are able to systematically combine high-dimensional voxel-level imaging and genetic data for accurate prediction of clinical outcomes of interest. Existing prediction models that integrate imaging and genetic features often use region level imaging summaries, and they typically do not consider the spatial configurations of the voxels in the image or incorporate the dependence between genes that may compromise prediction ability. We propose a novel integrative Bayesian scalar-on-image regression model for predicting cognitive outcomes based on high-dimensional spatially distributed voxel-level imaging data, along with correlated transcriptomic features. We account for the spatial dependencies in the imaging voxels via a tensor approach that also enables massive dimension reduction to address the curse of dimensionality, and models the dependencies between the transcriptomic features via a Graph-Laplacian prior. We implement this approach via an efficient Markov chain Monte Carlo (MCMC) computation strategy. We apply the proposed method to the analysis of longitudinal ADNI data for predicting cognitive scores at different visits by integrating voxel-level cortical thickness measurements derived from T1w-MRI scans and transcriptomics data. We illustrate that the proposed imaging transcriptomics approach has significant improvements in prediction compared to prediction using a subset of features from only one modality (imaging or genetics), as well as when using imaging and transcriptomics features but ignoring the inherent dependencies between the features. Our analysis is one of the first to conclusively demonstrate the advantages of prediction based on combining voxel-level cortical thickness measurements along with transcriptomics features, while accounting for inherent structural information.

Genetic Variants of Adrenoceptors

Adrenoceptors are class A G-protein-coupled receptors grouped into three families (α1-, α2-, and β-adrenoceptors), each one including three members. All nine corresponding adrenoceptor genes display genetic variation in their coding and adjacent non-coding genomic region. Coding variants, i.e., nucleotide exchanges within the transcribed and translated receptor sequence, may result in a difference in amino acid sequence thus altering receptor function and signaling. Such variants have been intensely studied in vitro in overexpression systems and addressed in candidate-gene studies for distinct clinical parameters. In recent years, large cohorts were analyzed in genome-wide association studies (GWAS), where variants are detected as significant in context with specific traits. These studies identified two of the in-depth characterized 18 coding variants in adrenoceptors as repeatedly statistically significant genetic risk factors - p.Arg389Gly in the β1- and p.Thr164Ile in the β2-adrenoceptor, along with 56 variants in the non-coding regions adjacent to the adrenoceptor gene loci, the functional role of which is largely unknown at present. This chapter summarizes current knowledge on the two coding variants in adrenoceptors that have been consistently validated in GWAS and provides a prospective overview on the numerous non-coding variants more recently attributed to adrenoceptor gene loci.

Genetic Factors Associated with Suicidal Behaviors and Alcohol Use Disorders in an American Indian Population

American Indians (AI) demonstrate the highest rates of both suicidal behaviors (SB) and alcohol use disorders (AUD) among all ethnic groups in the US. Rates of suicide and AUD vary substantially between tribal groups and across different geographical regions, underscoring a need to delineate more specific risk and resilience factors. Using data from over 740 AI living within eight contiguous reservations, we assessed genetic risk factors for SB by investigating: (1) possible genetic overlap with AUD, and (2) impacts of rare and low frequency genomic variants. Suicidal behaviors included lifetime history of suicidal thoughts and acts, including verified suicide deaths, scored using a ranking variable for the SB phenotype (range 0-4). We identified five loci significantly associated with SB and AUD, two of which are intergenic and three intronic on genes AACSP1, ANK1, and FBXO11. Nonsynonymous rare mutations in four genes including SERPINF1 (PEDF), ZNF30, CD34, and SLC5A9, and non-intronic rare mutations in genes OPRD1, HSD17B3 and one lincRNA were significantly associated with SB. One identified pathway related to hypoxia-inducible factor (HIF) regulation, whose 83 nonsynonymous rare variants on 10 genes were significantly linked to SB as well. Four additional genes, and two pathways related to vasopressin-regulated water metabolism and cellular hexose transport, also were strongly associated with SB. This study represents the first investigation of genetic factors for SB in an American Indian population that has high risk for suicide. Our study suggests that bivariate association analysis between comorbid disorders can increase statistical power; and rare variant analysis in a high-risk population enabled by whole-genome sequencing has the potential to identify novel genetic factors. Although such findings may be population specific, rare functional mutations relating to PEDF and HIF regulation align with past reports and suggest a biological mechanism for suicide risk and a potential therapeutic target for intervention.

Whole-exome identifies germline variants in families with obstructive sleep apnea syndrome

Background: Obstructive sleep apnea syndrome (OSAS) (OMIM #107650) is characterized by complete or partial obstruction of the upper airways, resulting in periods of sleep associated apnea. OSAS increases morbidity and mortality risk from cardiovascular and cerebrovascular diseases. While heritability of OSAS is estimated at ∼40%, the precise underlying genes remain elusive. Brazilian families with OSAS that follows as seemingly autosomal dominant inheritance pattern were recruited. Methods: The study included nine individuals from two Brazilian families displaying a seemingly autosomal dominant inheritance pattern of OSAS. Whole exome sequencing of germline DNA were analyzed using Mendel, MD software. Variants selected were analyzed using Varstation^® with subsequent analyses that included validation by Sanger sequencing, pathogenic score assessment by ACMG criteria, co-segregation analyses (when possible) allele frequency, tissue expression patterns, pathway analyses, effect on protein folding modeling using Swiss-Model and RaptorX. Results: Two families (six affected patients and three unaffected controls) were analyzed. A comprehensive multistep analysis yielded variants in COX20 (rs946982087) (family A), PTPDC1 (rs61743388) and TMOD4 (rs141507115) (family B) that seemed to be strong candidate genes for being OSAS associated genes in these families. Conclusion: Sequence variants in COX20, PTPDC1 and TMOD4 seemingly are associated with OSAS phenotype in these families. Further studies in more, ethnically diverse families and non-familial OSAS cases are needed to better define the role of these variants as contributors to OSAS phenotype.

The genomics of visuospatial neurocognition in obsessive-compulsive disorder: A preliminary GWAS

BACKGROUND

The study of Obsessive-Compulsive Disorder (OCD) genomics has primarily been tackled by Genome-wide association studies (GWAS), which have encountered troubles in identifying replicable single nucleotide polymorphisms (SNPs). Endophenotypes have emerged as a promising avenue of study in trying to elucidate the genomic bases of complex traits such as OCD.

METHODS

We analyzed the association of SNPs across the whole genome with the construction of visuospatial information and executive performance through four neurocognitive variables assessed by the Rey-Osterrieth Complex Figure Test (ROCFT) in a sample of 133 OCD probands. Analyses were performed at SNP- and gene-level.

RESULTS

No SNP reached genome-wide significance, although there was one SNP almost reaching significant association with copy organization (rs60360940; P = 9.98E-08). Suggestive signals were found for the four variables at both SNP- (P < 1E-05) and gene-levels (P < 1E-04). Most of the suggestive signals pointed to genes and genomic regions previously associated with neurological function and neuropsychological traits.

LIMITATIONS

Our main limitations were the sample size, which was limited to identify associated signals at a genome-wide level, and the composition of the sample, more representative of rather severe OCD cases than a population-based OCD sample with a broad severity spectrum.

CONCLUSIONS

Our results suggest that studying neurocognitive variables in GWAS would be more informative on the genetic basis of OCD than the classical case/control GWAS, facilitating the genetic characterization of OCD and its different clinical profiles, the development of individualized treatment approaches, and the improvement of prognosis and treatment response.

Large-scale genome sequencing redefines the genetic footprints of high-altitude adaptation in Tibetans

BACKGROUND

Tibetans are genetically adapted to high-altitude environments. Though many studies have been conducted, the genetic basis of the adaptation remains elusive due to the poor reproducibility for detecting selective signatures in the Tibetan genomes.

RESULTS

Here, we present whole-genome sequencing (WGS) data of 1001 indigenous Tibetans, covering the major populated areas of the Qinghai-Tibetan Plateau in China. We identify 35 million variants, and more than one-third of them are novel variants. Utilizing the large-scale WGS data, we construct a comprehensive map of allele frequency and linkage disequilibrium and provide a population-specific genome reference panel, referred to as 1KTGP. Moreover, with the use of a combined approach, we redefine the signatures of Darwinian-positive selection in the Tibetan genomes, and we characterize a high-confidence list of 4320 variants and 192 genes that have undergone selection in Tibetans. In particular, we discover four new genes, TMEM132C, ATP13A3, SANBR, and KHDRBS2, with strong signals of selection, and they may account for the adaptation of cardio-pulmonary functions in Tibetans. Functional annotation and enrichment analysis indicate that the 192 genes with selective signatures are likely involved in multiple organs and physiological systems, suggesting polygenic and pleiotropic effects.

CONCLUSIONS

Overall, the large-scale Tibetan WGS data and the identified adaptive variants/genes can serve as a valuable resource for future genetic and medical studies of high-altitude populations.

Deep Learning-Based Feature Extraction with MRI Data in Neuroimaging Genetics for Alzheimer's Disease

The prognosis and treatment of patients suffering from Alzheimer's disease (AD) have been among the most important and challenging problems over the last few decades. To better understand the mechanism of AD, it is of great interest to identify genetic variants associated with brain atrophy. Commonly, in these analyses, neuroimaging features are extracted based on one of many possible brain atlases with FreeSurf and other popular software; this, however, may cause the loss of important information due to our incomplete knowledge about brain function embedded in these suboptimal atlases. To address the issue, we propose convolutional neural network (CNN) models applied to three-dimensional MRI data for the whole brain or multiple, divided brain regions to perform completely data-driven and automatic feature extraction. These image-derived features are then used as endophenotypes in genome-wide association studies (GWASs) to identify associated genetic variants. When we applied this method to ADNI data, we identified several associated SNPs that have been previously shown to be related to several neurodegenerative/mental disorders, such as AD, depression, and schizophrenia.

The role of structural variations in Alzheimer's disease and other neurodegenerative diseases

Dozens of single nucleotide polymorphisms (SNPs) related to Alzheimer's disease (AD) have been discovered by large scale genome-wide association studies (GWASs). However, only a small portion of the genetic component of AD can be explained by SNPs observed from GWAS. Structural variation (SV) can be a major contributor to the missing heritability of AD; while SV in AD remains largely unexplored as the accurate detection of SVs from the widely used array-based and short-read technology are still far from perfect. Here, we briefly summarized the strengths and weaknesses of available SV detection methods. We reviewed the current landscape of SV analysis in AD and SVs that have been found associated with AD. Particularly, the importance of currently less explored SVs, including insertions, inversions, short tandem repeats, and transposable elements in neurodegenerative diseases were highlighted.

Genome-wide association study of brain biochemical phenotypes reveals distinct genetic architecture of Alzheimer's disease related proteins

BACKGROUND

Alzheimer's disease (AD) is neuropathologically characterized by amyloid-beta (Aβ) plaques and neurofibrillary tangles. The main protein components of these hallmarks include Aβ40, Aβ42, tau, phosphor-tau, and APOE. We hypothesize that genetic variants influence the levels and solubility of these AD-related proteins in the brain; identifying these may provide key insights into disease pathogenesis.

METHODS

Genome-wide genotypes were collected from 441 AD cases, imputed to the haplotype reference consortium (HRC) panel, and filtered for quality and frequency. Temporal cortex levels of five AD-related proteins from three fractions, buffer-soluble (TBS), detergent-soluble (Triton-X = TX), and insoluble (Formic acid = FA), were available for these same individuals. Variants were tested for association with each quantitative biochemical measure using linear regression, and GSA-SNP2 was used to identify enriched Gene Ontology (GO) terms. Implicated variants and genes were further assessed for association with other relevant variables.

RESULTS

We identified genome-wide significant associations at seven novel loci and the APOE locus. Genes and variants at these loci also associate with multiple AD-related measures, regulate gene expression, have cell-type specific enrichment, and roles in brain health and other neuropsychiatric diseases. Pathway analysis identified significant enrichment of shared and distinct biological pathways.

CONCLUSIONS

Although all biochemical measures tested reflect proteins core to AD pathology, our results strongly suggest that each have unique genetic architecture and biological pathways that influence their specific biochemical states in the brain. Our novel approach of deep brain biochemical endophenotype GWAS has implications for pathophysiology of proteostasis in AD that can guide therapeutic discovery efforts focused on these proteins.

CRISPRi screening reveals regulators of tau pathology shared between exosomal and vesicle-free tau

The aggregation of the microtubule-associated protein tau is a defining feature of Alzheimer's disease and other tauopathies. Tau pathology is believed to be driven by free tau aggregates and tau carried within exosome-like extracellular vesicles, both of which propagate trans-synaptically and induce tau pathology in recipient neurons by a corrupting process of seeding. Here, we performed a genome-wide CRISPRi screen in tau biosensor cells and identified cellular regulators shared by both mechanisms of tau seeding. We identified ANKLE2, BANF1, NUSAP1, EIF1AD, and VPS18 as the top validated regulators that restrict tau aggregation initiated by both exosomal and vesicle-free tau seeds. None of our validated hits affected the uptake of either form of tau seeds, supporting the notion that they operate through a cell-autonomous mechanism downstream of the seed uptake. Lastly, validation studies with human brain tissue also revealed that several of the identified protein hits are down-regulated in the brains of Alzheimer's patients, suggesting that their decreased activity may be required for the emergence or progression of tau pathology in the human brain.

Loss of SORCS2 is Associated with Neuronal DNA Double-Strand Breaks

SORCS2 is one of five proteins that constitute the Vps10p-domain receptor family. Members of this family play important roles in cellular processes linked to neuronal survival, differentiation and function. Genetic and functional studies implicate SORCS2 in cognitive function, as well as in neurodegenerative and psychiatric disorders. DNA damage and DNA repair deficits are linked to ageing and neurodegeneration, and transient neuronal DNA double-strand breaks (DSBs) also occur as a result of neuronal activity. Here, we report a novel role for SORCS2 in DSB formation. We show that SorCS2 loss is associated with elevated DSB levels in the mouse dentate gyrus and that knocking out SORCS2 in a human neuronal cell line increased Topoisomerase IIβ-dependent DSB formation and reduced neuronal viability. Neuronal stimulation had no impact on levels of DNA breaks in vitro, suggesting that the observed differences may not be the result of aberrant neuronal activity in these cells. Our findings are consistent with studies linking the VPS10 receptors and DNA damage to neurodegenerative conditions.

MED13 and glycolysis are conserved modifiers of α-synuclein-associated neurodegeneration

α-Synuclein (α-syn) is important in synucleinopathies such as Parkinson's disease (PD). While genome-wide association studies (GWASs) of synucleinopathies have identified many risk loci, the underlying genes have not been shown for most loci. Using Drosophila, we screened 3,471 mutant chromosomes for genetic modifiers of α-synuclein and identified 12 genes. Eleven modifiers have human orthologs associated with diseases, including MED13 and CDC27, which lie within PD GWAS loci. Drosophila Skd/Med13 and glycolytic enzymes are co-upregulated by α-syn-associated neurodegeneration. While elevated α-syn compromises mitochondrial function, co-expressing skd/Med13 RNAi and α-syn synergistically increase the ratio of oxidized-to-reduced glutathione. The resulting neurodegeneration can be suppressed by overexpressing a glycolytic enzyme or treatment with deferoxamine, suggesting that compensatory glycolysis is neuroprotective. In addition, the functional relationship between α-synuclein, MED13, and glycolytic enzymes is conserved between flies and mice. We propose that hypoxia-inducible factor and MED13 are part of a druggable pathway for PD.

Single Nucleus Transcriptome Data from Alzheimer’s Disease Mouse Models Yield New Insight into Pathophysiology

Background: 5XFAD humanized mutant mice and Trem2 knockout (T2KO) mice are two mouse models relevant to the study of Alzheimer’s disease (AD)-related pathology. Objective: To determine hippocampal transcriptomic and polyadenylation site usage alterations caused by genetic mutations engineered in 5XFAD and T2KO mice. Methods: Employing a publicly available single-nucleus RNA sequencing dataset, we used Seurat and Sierra analytic programs to identify differentially expressed genes (DEGs) and differential transcript usage (DTU), respectively, in hippocampal cell types from each of the two mouse models. We analyzed cell type-specific DEGs further using Ingenuity Pathway Analysis (IPA). Results: We identified several DEGs in both neuronal and glial cell subtypes in comparisons of wild type (WT) versus 5XFAD and WT versus T2KO mice, including Ttr, Fth1, Pcsk1n, Malat1, Rpl37, Rtn1, Sepw1, Uba52, Mbp, Arl6ip5, Gm26917, Vwa1, and Pgrmc1. We also observed DTU in common between the two comparisons in neuronal and glial subtypes, specifically in the genes Prnp, Rbm4b, Pnisr, Opcml, Cpne7, Adgrb1, Gabarapl2, Ubb, Ndfip1, Car11, and Stmn4. IPA identified three statistically significant canonical pathways that appeared in multiple cell types and that overlapped between 5XFAD and T2KO comparisons to WT, including ‘FXR/RXR Activation’, ‘LXR/RXR Activation’, and ‘Acute Phase Response Signaling’. Conclusion: DEG, DTU, and IPA findings, derived from two different mouse models of AD, highlight the importance of energy imbalance and inflammatory processes in specific hippocampal cell types, including subtypes of neurons and glial cells, in the development of AD-related pathology. Additional studies are needed to further characterize these findings.

Genome-wide analysis identified abundant genetic modulators of contributions of the apolipoprotein E alleles to Alzheimer's disease risk

INTRODUCTION

The apolipoprotein E (APOE) ε2 and ε4 alleles have beneficial and adverse impacts on Alzheimer's disease (AD), respectively, with incomplete penetrance, which may be modulated by other genetic variants.

METHODS

We examined whether the associations of the APOE alleles with other polymorphisms in the genome can be sensitive to AD-affection status.

RESULTS

We identified associations of the ε2 and ε4 alleles with 314 and 232 polymorphisms, respectively. Of them, 35 and 31 polymorphisms had significantly different effects in AD-affected and -unaffected groups, suggesting their potential involvement in the AD pathogenesis by modulating the effects of the ε2 and ε4 alleles, respectively. Our survival-type analysis of the AD risk supported modulating roles of multiple group-specific polymorphisms. Our functional analysis identified gene enrichment in multiple immune-related biological processes, for example, B cell function.

DISCUSSION

These findings suggest involvement of local and inter-chromosomal modulators of the effects of the APOE alleles on the AD risk.

Transcriptional profiles in olfactory pathway-associated brain regions of African green monkeys: Associations with age and Alzheimer's disease neuropathology

Introduction

Olfactory impairment in older individuals is associated with an increased risk of Alzheimer's disease (AD). Characterization of age versus neuropathology‐associated changes in the brain olfactory pathway may elucidate processes underlying early AD pathogenesis. Here, we report age versus AD neuropathology–associated differential transcription in four brain regions in the olfactory pathway of 10 female African green monkeys (vervet, Chlorocebus aethiops sabaeus), a well‐described model of early AD‐like neuropathology.

Methods

Transcriptional profiles were determined by microarray in the olfactory bulb (OB), piriform cortex (PC), temporal lobe white matter (WM), and inferior temporal cortex (ITC). Amyloid beta (Aβ) plaque load in parietal and temporal cortex was determined by immunohistochemistry, and concentrations of Aβ42, Aβ40, and norepinephrine in ITC were determined by enzyme‐linked immuosorbent assay (ELISA). Transcriptional profiles were compared between middle‐aged and old animals, and associations with AD‐relevant neuropathological measures were determined.

Results

Transcriptional profiles varied by brain region and age group. Expression levels of TRO and RNU4‐1 were significantly lower in all four regions in the older group. An additional 29 genes were differentially expressed by age in three of four regions. Analyses of a combined expression data set of all four regions identified 77 differentially expressed genes (DEGs) by age group. Among these DEGs, older subjects had elevated levels of CTSB, EBAG9, LAMTOR3, and MRPL17, and lower levels of COMMD10 and TYW1B. A subset of these DEGs was associated with neuropathology biomarkers. Notably, CTSB was positively correlated with Aβ plaque counts, Aβ42:Aβ40 ratios, and norepinephrine levels in all brain regions.

Discussion

These data demonstrate age differences in gene expression in olfaction‐associated brain regions. Biological processes exhibiting age‐related enrichment included the regulation of cell death, vascular function, mitochondrial function, and proteostasis. A subset of DEGs was specifically associated with AD phenotypes. These may represent promising targets for future mechanistic investigations and perhaps therapeutic intervention.

Epigenome-wide association study of global cortical volumes in generation Scotland: Scottish family health study

A complex interplay of genetic and environmental risk factors influence global brain structural alterations associated with brain health and disease. Epigenome-wide association studies (EWAS) of global brain imaging phenotypes have the potential to reveal the mechanisms of brain health and disease and can lead to better predictive analytics through the development of risk scores.We perform an EWAS of global brain volumes in Generation Scotland using peripherally measured whole blood DNA methylation (DNAm) from two assessments, (i) at baseline recruitment, ~6 years prior to MRI assessment (N = 672) and (ii) concurrent with MRI assessment (N=565). Four CpGs at baseline were associated with global cerebral white matter, total grey matter, and whole-brain volume (Bonferroni p≤7.41×10-8, βrange = -1.46x10-6 to 9.59 × 10-7). These CpGs were annotated to genes implicated in brain-related traits, including psychiatric disorders, development, and ageing. We did not find significant associations in the meta-analysis of the EWAS of the two sets concurrent with imaging at the corrected level.These findings reveal global brain structural changes associated with DNAm measured ~6 years previously, indicating a potential role of early DNAm modifications in brain structure. Although concurrent DNAm was not associated with global brain structure, the nominally significant findings identified here present a rationale for future investigation of associations between DNA methylation and structural brain phenotypes in larger population-based samples.

Deep multiview learning to identify imaging-driven subtypes in mild cognitive impairment

BACKGROUND

In Alzheimer's Diseases (AD) research, multimodal imaging analysis can unveil complementary information from multiple imaging modalities and further our understanding of the disease. One application is to discover disease subtypes using unsupervised clustering. However, existing clustering methods are often applied to input features directly, and could suffer from the curse of dimensionality with high-dimensional multimodal data. The purpose of our study is to identify multimodal imaging-driven subtypes in Mild Cognitive Impairment (MCI) participants using a multiview learning framework based on Deep Generalized Canonical Correlation Analysis (DGCCA), to learn shared latent representation with low dimensions from 3 neuroimaging modalities.

RESULTS

DGCCA applies non-linear transformation to input views using neural networks and is able to learn correlated embeddings with low dimensions that capture more variance than its linear counterpart, generalized CCA (GCCA). We designed experiments to compare DGCCA embeddings with single modality features and GCCA embeddings by generating 2 subtypes from each feature set using unsupervised clustering. In our validation studies, we found that amyloid PET imaging has the most discriminative features compared with structural MRI and FDG PET which DGCCA learns from but not GCCA. DGCCA subtypes show differential measures in 5 cognitive assessments, 6 brain volume measures, and conversion to AD patterns. In addition, DGCCA MCI subtypes confirmed AD genetic markers with strong signals that existing late MCI group did not identify.

CONCLUSION

Overall, DGCCA is able to learn effective low dimensional embeddings from multimodal data by learning non-linear projections. MCI subtypes generated from DGCCA embeddings are different from existing early and late MCI groups and show most similarity with those identified by amyloid PET features. In our validation studies, DGCCA subtypes show distinct patterns in cognitive measures, brain volumes, and are able to identify AD genetic markers. These findings indicate the promise of the imaging-driven subtypes and their power in revealing disease structures beyond early and late stage MCI.

Systems Genetic Identification of Mitochondrion-Associated Alzheimer's Disease Genes and Implications for Disease Risk Prediction

Cumulative evidence has revealed the association between mitochondrial dysfunction and Alzheimer’s disease (AD). Because the number of mitochondrial genes is very limited, the mitochondrial pathogenesis of AD must involve certain nuclear genes. In this study, we employed systems genetic methods to identify mitochondrion-associated nuclear genes that may participate in the pathogenesis of AD. First, we performed a mitochondrial genome-wide association study (MiWAS, n = 809) to identify mitochondrial single-nucleotide polymorphisms (MT-SNPs) associated with AD. Then, epistasis analysis was performed to examine interacting SNPs between the mitochondrial and nuclear genomes. Weighted co-expression network analysis (WGCNA) was applied to transcriptomic data from the same sample (n = 743) to identify AD-related gene modules, which were further enriched by mitochondrion-associated genes. Using hub genes derived from these modules, random forest models were constructed to predict AD risk in four independent datasets (n = 743, n = 542, n = 161, and n = 540). In total, 9 potentially significant MT-SNPs and 14,340 nominally significant MT-nuclear interactive SNPs were identified for AD, which were validated by functional analysis. A total of 6 mitochondrion-related modules involved in AD pathogenesis were found by WGCNA, from which 91 hub genes were screened and used to build AD risk prediction models. For the four independent datasets, these models perform better than those derived from AD genes identified by genome-wide association studies (GWASs) or differential expression analysis (DeLong’s test, p  <  0.05). Overall, through systems genetics analyses, mitochondrion-associated SNPs/genes with potential roles in AD pathogenesis were identified and preliminarily validated, illustrating the power of mitochondrial genetics in AD pathogenesis elucidation and risk prediction.

Inter- and intra-chromosomal modulators of the APOE ɛ2 and ɛ4 effects on the Alzheimer's disease risk

The mechanisms of incomplete penetrance of risk-modifying impacts of apolipoprotein E (APOE) ε2 and ε4 alleles on Alzheimer's disease (AD) have not been fully understood. We performed genome-wide analysis of differences in linkage disequilibrium (LD) patterns between 6,136 AD-affected and 10,555 AD-unaffected subjects from five independent studies to explore whether the association of the APOE ε2 allele (encoded by rs7412 polymorphism) and ε4 allele (encoded by rs429358 polymorphism) with AD was modulated by autosomal polymorphisms. The LD analysis identified 24 (mostly inter-chromosomal) and 57 (primarily intra-chromosomal) autosomal polymorphisms with significant differences in LD with either rs7412 or rs429358, respectively, between AD-affected and AD-unaffected subjects, indicating their potential modulatory roles. Our Cox regression analysis showed that minor alleles of four inter-chromosomal and ten intra-chromosomal polymorphisms exerted significant modulating effects on the ε2- and ε4-associated AD risks, respectively, and identified ε2-independent (rs2884183 polymorphism, 11q22.3) and ε4-independent (rs483082 polymorphism, 19q13.32) associations with AD. Our functional analysis highlighted ε2- and/or ε4-linked processes affecting the lipid and lipoprotein metabolism and cell junction organization which may contribute to AD pathogenesis. These findings provide insights into the ε2- and ε4-associated mechanisms of AD pathogenesis, underlying their incomplete penetrance.

Integration of multidimensional splicing data and GWAS summary statistics for risk gene discovery

A common strategy for the functional interpretation of genome-wide association study (GWAS) findings has been the integrative analysis of GWAS and expression data. Using this strategy, many association methods (e.g., PrediXcan and FUSION) have been successful in identifying trait-associated genes via mediating effects on RNA expression. However, these approaches often ignore the effects of splicing, which can carry as much disease risk as expression. Compared to expression data, one challenge to detect associations using splicing data is the large multiple testing burden due to multidimensional splicing events within genes. Here, we introduce a multidimensional splicing gene (MSG) approach, which consists of two stages: 1) we use sparse canonical correlation analysis (sCCA) to construct latent canonical vectors (CVs) by identifying sparse linear combinations of genetic variants and splicing events that are maximally correlated with each other; and 2) we test for the association between the genetically regulated splicing CVs and the trait of interest using GWAS summary statistics. Simulations show that MSG has proper type I error control and substantial power gains over existing multidimensional expression analysis methods (i.e., S-MultiXcan, UTMOST, and sCCA+ACAT) under diverse scenarios. When applied to the Genotype-Tissue Expression Project data and GWAS summary statistics of 14 complex human traits, MSG identified on average 83%, 115%, and 223% more significant genes than sCCA+ACAT, S-MultiXcan, and UTMOST, respectively. We highlight MSG's applications to Alzheimer's disease, low-density lipoprotein cholesterol, and schizophrenia, and found that the majority of MSG-identified genes would have been missed from expression-based analyses. Our results demonstrate that aggregating splicing data through MSG can improve power in identifying gene-trait associations and help better understand the genetic risk of complex traits.

Replication of gene polymorphisms associated with periodontitis-related traits in an elderly cohort: the Washington Heights/Inwood Community Aging Project Ancillary Study of Oral Health

AIM

We sought to replicate findings from published genome-wide association studies (GWAS), linking specific candidate gene loci with periodontitis-related clinical/microbial traits.

MATERIALS AND METHODS

In the published GWAS, a total of 2196 single nucleotide polymorphisms associated with periodontitis-related traits at a p ≤ 5 × 10^-6 and mapped to 136 gene loci. The replication cohort included 1124 individuals, 65-98 years old (67% female, 45% Hispanic, 30% Black, 23% White) with available genome-wide genotypes and full-mouth periodontal status. Microbial profiles using checkerboard DNA-DNA hybridization and 16SrRNA sequencing were available from 912 and 739 participants, respectively.

RESULTS

Using gene-specific p-values after linkage disequilibrium pruning, the following gene/phenotype associations replicated successfully: CLEC19A with edentulism and %teeth with pocket depth (PD) ≥4 mm; IL37, HPVC1, TRPS1, ABHD12B, LDLRAD4 (C180rF1), TGM3, and GRK5 with %teeth with PD ≥4 mm; DAB2IP with presence of PD ≥6 mm; KIAA1715(LNPK), ROBO2, RAB28, LINC01017, NELL1, LDLRAD4(C18orF1), and CRYBB2P1 with %teeth with clinical attachment level (CAL) ≥3 mm; RUNX2 and LAMA2 with %teeth with CAL ≥5 mm; and KIAA1715(LNPK) with high colonization by Aggregatibacter actinomycetemcomitans. In addition, CLEC19A, IQSEC1, and EMR1 associated with microbial abundance based on checkerboard data, LBP and NCR2 with abundance based on sequencing data, and NCR2 with microbial diversity based on sequencing data.

CONCLUSIONS

Several gene loci identified in published GWAS as associated with periodontitis-related phenotypes replicated successfully in an elderly cohort.

Using the Alzheimer's Disease Neuroimaging Initiative to improve early detection, diagnosis, and treatment of Alzheimer's disease

INTRODUCTION

The Alzheimer's Disease Neuroimaging Initiative (ADNI) has accumulated 15 years of clinical, neuroimaging, cognitive, biofluid biomarker and genetic data, and biofluid samples available to researchers, resulting in more than 3500 publications. This review covers studies from 2018 to 2020.

METHODS

We identified 1442 publications using ADNI data by conventional search methods and selected impactful studies for inclusion.

RESULTS

Disease progression studies supported pivotal roles for regional amyloid beta (Aβ) and tau deposition, and identified underlying genetic contributions to Alzheimer's disease (AD). Vascular disease, immune response, inflammation, resilience, and sex modulated disease course. Biologically coherent subgroups were identified at all clinical stages. Practical algorithms and methodological changes improved determination of Aβ status. Plasma Aβ, phosphorylated tau181, and neurofilament light were promising noninvasive biomarkers. Prognostic and diagnostic models were externally validated in ADNI but studies are limited by lack of ethnocultural cohort diversity.

DISCUSSION

ADNI has had a profound impact in improving clinical trials for AD.

Macrophage mitochondrial bioenergetics and tissue invasion are boosted by an Atossa-Porthos axis in Drosophila

Cellular metabolism must adapt to changing demands to enable homeostasis. During immune responses or cancer metastasis, cells leading migration into challenging environments require an energy boost, but what controls this capacity is unclear. Here, we study a previously uncharacterized nuclear protein, Atossa (encoded by CG9005), which supports macrophage invasion into the germband of Drosophila by controlling cellular metabolism. First, nuclear Atossa increases mRNA levels of Porthos, a DEAD‐box protein, and of two metabolic enzymes, lysine‐α‐ketoglutarate reductase (LKR/SDH) and NADPH glyoxylate reductase (GR/HPR), thus enhancing mitochondrial bioenergetics. Then Porthos supports ribosome assembly and thereby raises the translational efficiency of a subset of mRNAs, including those affecting mitochondrial functions, the electron transport chain, and metabolism. Mitochondrial respiration measurements, metabolomics, and live imaging indicate that Atossa and Porthos power up OxPhos and energy production to promote the forging of a path into tissues by leading macrophages. Since many crucial physiological responses require increases in mitochondrial energy output, this previously undescribed genetic program may modulate a wide range of cellular behaviors.

Association between WWOX/MAF variants and dementia-related neuropathologic endophenotypes

The genetic locus containing the WWOX and MAF genes was implicated as a clinical Alzheimer's disease (AD) risk locus in two recent large meta-analytic genome wide association studies (GWAS). In a prior GWAS, we identified a variant in WWOX as a suggestive risk allele for hippocampal sclerosis. We hypothesized that the WWOX/MAF locus may be preferentially associated with non-plaque- and non-tau-related neuropathological changes (NC). Data from research participants with GWAS and autopsy measures from the National Alzheimer's Coordinating Center and the Religious Orders Study and the Rush Memory and Aging Project were meta-analyzed. Notably, no variants in the locus were significantly associated with ADNC. However, several WWOX/MAF variants had significant adjusted associations with limbic-predominant age-related TDP-43 encephalopathy NC (LATE-NC), HS, and brain arteriolosclerosis. These associations remained largely unchanged after adjustment for ADNC (operationalized with standard semiquantitative staging), suggesting that these associations are independent of ADNC. Thus, WWOX genetic variants were associated pathologically with LATE-NC, not ADNC.

Methylome-wide association study of antidepressant use in Generation Scotland and the Netherlands Twin Register implicates the innate immune system

Antidepressants are an effective treatment for major depressive disorder (MDD), although individual response is unpredictable and highly variable. Whilst the mode of action of antidepressants is incompletely understood, many medications are associated with changes in DNA methylation in genes that are plausibly linked to their mechanisms. Studies of DNA methylation may therefore reveal the biological processes underpinning the efficacy and side effects of antidepressants. We performed a methylome-wide association study (MWAS) of self-reported antidepressant use accounting for lifestyle factors and MDD in Generation Scotland (GS:SFHS, N = 6428, EPIC array) and the Netherlands Twin Register (NTR, N = 2449, 450 K array) and ran a meta-analysis of antidepressant use across these two cohorts. We found ten CpG sites significantly associated with self-reported antidepressant use in GS:SFHS, with the top CpG located within a gene previously associated with mental health disorders, ATP6V1B2 (β = -0.055, pcorrected = 0.005). Other top loci were annotated to genes including CASP10, TMBIM1, MAPKAPK3, and HEBP2, which have previously been implicated in the innate immune response. Next, using penalised regression, we trained a methylation-based score of self-reported antidepressant use in a subset of 3799 GS:SFHS individuals that predicted antidepressant use in a second subset of GS:SFHS (N = 3360, β = 0.377, p = 3.12 × 10-11, R2 = 2.12%). In an MWAS analysis of prescribed selective serotonin reuptake inhibitors, we showed convergent findings with those based on self-report. In NTR, we did not find any CpGs significantly associated with antidepressant use. The meta-analysis identified the two CpGs of the ten above that were common to the two arrays used as being significantly associated with antidepressant use, although the effect was in the opposite direction for one of them. Antidepressants were associated with epigenetic alterations in loci previously associated with mental health disorders and the innate immune system. These changes predicted self-reported antidepressant use in a subset of GS:SFHS and identified processes that may be relevant to our mechanistic understanding of clinically relevant antidepressant drug actions and side effects.

Past, present and future of therapeutic strategies against amyloid-β peptides in Alzheimer's disease: a systematic review

Alzheimer's disease (AD) is the most prevalent neurodegenerative disease in ageing, affecting around 46 million people worldwide but few treatments are currently available. The etiology of AD is still puzzling, and new drugs development and clinical trials have high failure rates. Urgent outline of an integral (multi-target) and effective treatment of AD is needed. Accumulation of amyloid-β (Aβ) peptides is considered one of the fundamental neuropathological pillars of the disease, and its dyshomeostasis has shown a crucial role in AD onset. Therefore, many amyloid-targeted therapies have been investigated. Here, we will systematically review recent (from 2014) investigational, follow-up and review studies focused on anti-amyloid strategies to summarize and analyze their current clinical potential. Combination of anti-Aβ therapies with new developing early detection biomarkers and other therapeutic agents acting on early functional AD changes will be highlighted in this review. Near-term approval seems likely for several drugs acting against Aβ, with recent FDA approval of a monoclonal anti-Aβ oligomers antibody -aducanumab- raising hopes and controversies. We conclude that, development of oligomer-epitope specific Aβ treatment and implementation of multiple improved biomarkers and risk prediction methods allowing early detection, together with therapies acting on other factors such as hyperexcitability in early AD, could be the key to slowing this global pandemic.

The Fly Homologue of MFSD11 Is Possibly Linked to Nutrient Homeostasis and Has a Potential Role in Locomotion: A First Characterization of the Atypical Solute Carrier CG18549 in Drosophila Melanogaster

Simple Summary

The body is dependent on nutrients and ions to work normally. Within the hu-man body there is a group of proteins named transporters or solute carriers. These transporters are vital for the transport of nutrients such as glucose, amino acids, and fats, as well as ions such as sodium, calcium, and potassium. Despite being vital for normal physiology, as well as pathophysiology, a large number (approximately one third) of the transporters are orphans, where information about their expression and function is missing. Here, we aimed to begin to unravel the expression and function of one of these orphan transporters, MFSD11, by studying its orthologue in fruit flies (CG18549). We found that the fly orthologue is expressed in the brain of fruit flies and that it is possibly involved in metabolism and/or locomotion of the flies. The exact mechanism behind the observed behaviors is not fully understood, but our study provides new insights into the expression and function of CG18549. Clearly, these results, among others about the orphan transporters, provide a strong example as to why it is vital to fully characterize them and through that gain knowledge about the body during normal condition and disease.

Abstract

Cellular transport and function are dependent on substrate influx and efflux of various compounds. In humans, the largest superfamily of transporters is the SoLute Carriers (SLCs). Many transporters are orphans and little to nothing is known about their expression and/or function, yet they have been assigned to a cluster called atypical SLCs. One of these atypical SLCs is MFSD11. Here we present a first in-depth characterization of the MFSD11, CG18549. By gene expression and behavior analysis on ubiquitous and brain-specific knockdown flies. CG18549 knockdown flies were found to have altered adipokinetic hormone and adipokinteic hormone receptor expression as well as reduced vesicular monoamine transporter expression; to exhibit an altered locomotor behavior, and to have an altered reaction to stress stimuli. Furthermore, the gene expression of CG18549 in the brain was visualized and abundant expression in both the larvae and adult brain was observed, a result that is coherent with the FlyAtlas Anatomy microarray. The exact mechanism behind the observed behaviors is not fully understood, but this study provides new insights into the expression and function of CG18549. Clearly, these results provide a strong example as to why it is vital to fully characterize orphan transporters and through that gain knowledge about the body during normal condition and disease.

The heterogeneity among subgroups of haplogroup J influencing Alzheimer's disease risk

Introduction

The impact of mitochondrial haplogroups on Alzheimer's disease (AD) risk has not been fully elucidated and warrants further investigation at the subgroup level.

Objectives

The aim of this research is to evaluate the association between mitochondrial haplogroups and AD risk in subgroups level.

Methods

In total, 809 AD Neuroimaging Initiative subjects were assessed using mtDNA sequencing, the AD Assessment Scale-Cognitive Subscale (ADAS-cog), hippocampal volume measurements, the hypometabolic convergence index (HCI), and MCI-to-AD conversion proportion measurements.

Results

The frequency of haplogroup J was significantly higher than that of other haplogroups in the AD group (p = 0.013). According to the correlation between haplogroup J-specific SNPs and ADAS-cog, haplogroup J was divided into four subgroups harboring exacerbating SNPs, protective SNPs, both exacerbating and protective SNPs, or irrelevant SNPs. The subgroups harboring exacerbating SNPs exhibited higher AD risk represented by the levels of ADAS-cog, hippocampal volume, HCI, and MCI-to-AD conversion proportion than other subgroups.

Conclusion

Heterogeneity existed among the subgroups of haplogroup J, which suggested that different subgroups exhibited different levels of AD risk. This study provides novel insights into the correlation between mitochondrial haplogroups and AD risk.

RP11-362K2.2:RP11-767I20.1 Genetic Variation Is Associated with Post-Reperfusion Therapy Parenchymal Hematoma. A GWAS Meta-Analysis

Stroke is one of the most common causes of death and disability. Reperfusion therapies are the only treatment available during the acute phase of stroke. Due to recent clinical trials, these therapies may increase their frequency of use by extending the time-window administration, which may lead to an increase in complications such as hemorrhagic transformation, with parenchymal hematoma (PH) being the more severe subtype, associated with higher mortality and disability rates. Our aim was to find genetic risk factors associated with PH, as that could provide molecular targets/pathways for their prevention/treatment and study its genetic correlations to find traits sharing genetic background. We performed a GWAS and meta-analysis, following standard quality controls and association analysis (fastGWAS), adjusting age, NIHSS, and principal components. FUMA was used to annotate, prioritize, visualize, and interpret the meta-analysis results. The total number of patients in the meta-analysis was 2034 (216 cases and 1818 controls). We found rs79770152 having a genome-wide significant association (beta 0.09, p-value 3.90 × 10⁻⁸) located in the RP11-362K2.2:RP11-767I20.1 gene and a suggestive variant (rs13297983: beta 0.07, p-value 6.10 × 10⁻⁸) located in PCSK5 associated with PH occurrence. The genetic correlation showed a shared genetic background of PH with Alzheimer’s disease and white matter hyperintensities. In addition, genes containing the ten most significant associations have been related to aggregated amyloid-β, tau protein, white matter microstructure, inflammation, and matrix metalloproteinases.

Interplay between stress-related genes may influence Alzheimer's disease development: The results of genetic interaction analyses of human data

Emerging evidence from experimental and clinical research suggests that stress-related genes may play key roles in AD development. The fact that genome-wide association studies were not able to detect a contribution of such genes to AD indicates the possibility that these genes may influence AD non-linearly, through interactions of their products. In this paper, we selected two stress-related genes (GCN2/EIF2AK4 and APP) based on recent findings from experimental studies which suggest that the interplay between these genes might influence AD in humans. To test this hypothesis, we evaluated the effects of interactions between SNPs in these two genes on AD occurrence, using the Health and Retirement Study data on white indidividuals. We found several interacting SNP-pairs whose associations with AD remained statistically significant after correction for multiple testing. These findings emphasize the importance of nonlinear mechanisms of polygenic AD regulation that cannot be detected in traditional association studies. To estimate collective effects of multiple interacting SNP-pairs on AD, we constructed a new composite index, called Interaction Polygenic Risk Score, and showed that its association with AD is highly statistically significant. These results open a new avenue in the analyses of mechanisms of complex multigenic AD regulation.

The Emerging Physiological Role of AGMO 10 Years after Its Gene Identification

The gene encoding alkylglycerol monooxygenase (AGMO) was assigned 10 years ago. So far, AGMO is the only known enzyme capable of catalysing the breakdown of alkylglycerols and lyso-alkylglycerophospholipids. With the knowledge of the genetic information, it was possible to relate a potential contribution for mutations in the AGMO locus to human diseases by genome-wide association studies. A possible role for AGMO was implicated by genetic analyses in a variety of human pathologies such as type 2 diabetes, neurodevelopmental disorders, cancer, and immune defence. Deficient catabolism of stored lipids carrying an alkyl bond by an absence of AGMO was shown to impact on the overall lipid composition also outside the ether lipid pool. This review focuses on the current evidence of AGMO in human diseases and summarises experimental evidence for its role in immunity, energy homeostasis, and development in humans and several model organisms. With the progress in lipidomics platform and genetic identification of enzymes involved in ether lipid metabolism such as AGMO, it is now possible to study the consequence of gene ablation on the global lipid pool and further on certain signalling cascades in a variety of model organisms in more detail.

Combined change of behavioral traits for domestication and gene-networks in mice selectively bred for active tameness

Tameness is a major element of animal domestication and involves two components: motivation to approach humans (active tameness) and reluctance to avoid humans (passive tameness). To understand the behavioral and genetic mechanisms of active tameness in mice, we had previously conducted selective breeding for long durations of contact and heading toward human hands in an active tameness test using a wild-derived heterogeneous stock. Although the study showed a significant increase in contacting and heading with the 12th generation of breeding, the effect on other behavioral indices related to tameness and change of gene expression levels underlying selective breeding was unclear. Here, we analyzed nine tameness-related traits at a later stage of selective breeding and analyzed how gene expression levels were changed by the selective breeding. We found that five traits, including contacting and heading, showed behavioral change in the selective groups comparing to the control through the generations. Furthermore, we conducted cluster analyses to evaluate the relationships among the nine traits and found that contacting and heading combined in an independent cluster in the selected groups, but not in the control groups. RNA-Seq of hippocampal tissue revealed differential expression of 136 genes between the selection and control groups, while the pathway analysis identified the networks associated with these genes. These results suggest that active tameness was hidden in the control groups but became apparent in the selected populations by selective breeding, potentially driven by changes in gene expression networks.