Systematic tissue-specific functional annotation of the human genome highlights immune-related DNA elements for late-onset Alzheimer's disease

Brain Cell-based Genetic Subtyping and Drug Repositioning for Alzheimer Disease

Alzheimer's Disease (AD) is characterized by its complex and heterogeneous etiology and gradual progression, leading to high drug failure rates in late-stage clinical trials. In order to better stratify individuals at risk for AD and discern potential therapeutic targets we employed a novel procedure utilizing cell-based co-regulated gene networks and polygenic risk scores (cbPRSs). After defining genetic subtypes using extremes of cbPRS distributions, we evaluated correlations of the genetic subtypes with previously defined AD subtypes defined on the basis of domain-specific cognitive functioning and neuroimaging biomarkers. Employing a PageRank algorithm, we identified priority gene targets for the genetic subtypes. Pathway analysis of priority genes demonstrated associations with neurodegeneration and suggested candidate drugs currently utilized in diabetes, hypertension, and epilepsy for repositioning in AD. Experimental validation utilizing human induced pluripotent stem cell (hiPSC)-derived astrocytes demonstrated the modifying effects of estradiol, levetiracetam, and pioglitazone on expression of APOE and complement C4 genes, suggesting potential repositioning for AD.

NCAD v1.0: a database for non-coding variant annotation and interpretation

The application of whole genome sequencing is expanding in clinical diagnostics across various genetic disorders, and the significance of non-coding variants in penetrant diseases is increasingly being demonstrated. Therefore, it is urgent to improve the diagnostic yield by exploring the pathogenic mechanisms of variants in non-coding regions. However, the interpretation of non-coding variants remains a significant challenge, due to the complex functional regulatory mechanisms of non-coding regions and the current limitations of available databases and tools. Hence, we develop the non-coding variant annotation database (NCAD, http://www.ncawdb.net/), encompassing comprehensive insights into 665,679,194 variants, regulatory elements, and element interaction details. Integrating data from 96 sources, spanning both GRCh37 and GRCh38 versions, NCAD v1.0 provides vital information to support the genetic diagnosis of non-coding variants, including allele frequencies of 12 diverse populations, with a particular focus on the population frequency information for 230,235,698 variants in 20,964 Chinese individuals. Moreover, it offers prediction scores for variant functionality, five categories of regulatory elements, and four types of non-coding RNAs. With its rich data and comprehensive coverage, NCAD serves as a valuable platform, empowering researchers and clinicians with profound insights into non-coding regulatory mechanisms while facilitating the interpretation of non-coding variants.

eQTL studies: from bulk tissues to single cells

An expression quantitative trait locus (eQTL) is a chromosomal region where genetic variants are associated with the expression levels of specific genes that can be both nearby or distant. The identifications of eQTLs for different tissues, cell types, and contexts have led to a better understanding of the dynamic regulations of gene expressions and implications of functional genes and variants for complex traits and diseases. Although most eQTL studies have been performed on data collected from bulk tissues, recent studies have demonstrated the importance of cell-type-specific and context-dependent gene regulations in biological processes and disease mechanisms. In this review, we discuss statistical methods that have been developed to enable the detection of cell-type-specific and context-dependent eQTLs from bulk tissues, purified cell types, and single cells. We also discuss the limitations of the current methods and future research opportunities.

multi-GPA-Tree: Statistical approach for pleiotropy informed and functional annotation tree guided prioritization of GWAS results

Genome-wide association studies (GWAS) have successfully identified over two hundred thousand genotype-trait associations. Yet some challenges remain. First, complex traits are often associated with many single nucleotide polymorphisms (SNPs), most with small or moderate effect sizes, making them difficult to detect. Second, many complex traits share a common genetic basis due to 'pleiotropy' and and though few methods consider it, leveraging pleiotropy can improve statistical power to detect genotype-trait associations with weaker effect sizes. Third, currently available statistical methods are limited in explaining the functional mechanisms through which genetic variants are associated with specific or multiple traits. We propose multi-GPA-Tree to address these challenges. The multi-GPA-Tree approach can identify risk SNPs associated with single as well as multiple traits while also identifying the combinations of functional annotations that can explain the mechanisms through which risk-associated SNPs are linked with the traits. First, we implemented simulation studies to evaluate the proposed multi-GPA-Tree method and compared its performance with existing statistical approaches. The results indicate that multi-GPA-Tree outperforms existing statistical approaches in detecting risk-associated SNPs for multiple traits. Second, we applied multi-GPA-Tree to a systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA), and to a Crohn's disease (CD) and ulcertive colitis (UC) GWAS, and functional annotation data including GenoSkyline and GenoSkylinePlus. Our results demonstrate that multi-GPA-Tree can be a powerful tool that improves association mapping while facilitating understanding of the underlying genetic architecture of complex traits and potential mechanisms linking risk-associated SNPs with complex traits.

Evaluating 17 methods incorporating biological function with GWAS summary statistics to accelerate discovery demonstrates a tradeoff between high sensitivity and high positive predictive value

Where sufficiently large genome-wide association study (GWAS) samples are not currently available or feasible, methods that leverage increasing knowledge of the biological function of variants may illuminate discoveries without increasing sample size. We comprehensively evaluated 17 functional weighting methods for identifying novel associations. We assessed the performance of these methods using published results from multiple GWAS waves across each of five complex traits. Although no method achieved both high sensitivity and positive predictive value (PPV) for any trait, a subset of methods utilizing pleiotropy and expression quantitative trait loci nominated variants with high PPV (>75%) for multiple traits. Application of functionally weighting methods to enhance GWAS power for locus discovery is unlikely to circumvent the need for larger sample sizes in truly underpowered GWAS, but these results suggest that applying functional weighting to GWAS can accurately nominate additional novel loci from available samples for follow-up studies.

Integrative genome-wide analyses identify novel loci associated with kidney stones and provide insights into its genetic architecture

Kidney stone disease (KSD) is a complex disorder with high heritability and prevalence. We performed a large genome-wide association study (GWAS) meta-analysis for KSD to date, including 720,199 individuals with 17,969 cases in European population. We identified 44 susceptibility loci, including 28 novel loci. Cell type-specific analysis pinpointed the proximal tubule as the most relevant cells where susceptibility variants might act through a tissue-specific fashion. By integrating kidney-specific omics data, we prioritized 223 genes which strengthened the importance of ion homeostasis, including calcium and magnesium in stone formation, and suggested potential target drugs for the treatment. The genitourinary and digestive diseases showed stronger genetic correlations with KSD. In this study, we generate an atlas of candidate genes, tissue and cell types involved in the formation of KSD. In addition, we provide potential drug targets for KSD treatment and insights into shared regulation with other diseases.

Alzheimer's disease heterogeneity explained by polygenic risk scores derived from brain transcriptomic profiles

INTRODUCTION

Alzheimer's disease (AD) is heterogeneous, both clinically and neuropathologically. We investigated whether polygenic risk scores (PRSs) integrated with transcriptome profiles from AD brains can explain AD clinical heterogeneity.

METHODS

We conducted co-expression network analysis and identified gene sets (modules) that were preserved in three AD transcriptome datasets and associated with AD-related neuropathological traits including neuritic plaques (NPs) and neurofibrillary tangles (NFTs). We computed the module-based PRSs (mbPRSs) for each module and tested associations with mbPRSs for cognitive test scores, cognitively defined AD subgroups, and brain imaging data.

RESULTS

Of the modules significantly associated with NPs and/or NFTs, the mbPRSs from two modules (M6 and M9) showed distinct associations with language and visuospatial functioning, respectively. They matched clinical subtypes and brain atrophy at specific regions.

DISCUSSION

Our findings demonstrate that polygenic profiling based on co-expressed gene sets can explain heterogeneity in AD patients, enabling genetically informed patient stratification and precision medicine in AD.

HIGHLIGHTS

Co-expression gene-network analysis in Alzheimer's disease (AD) brains identified gene sets (modules) associated with AD heterogeneity. AD-associated modules were selected when genes in each module were enriched for neuritic plaques and neurofibrillary tangles. Polygenic risk scores from two selected modules were linked to the matching cognitively defined AD subgroups (language and visuospatial subgroups). Polygenic risk scores from the two modules were associated with cognitive performance in language and visuospatial domains and the associations were confirmed in regional-specific brain atrophy data.

Epigenome-wide association study of peripheral immune cell populations in Parkinson's disease

Understanding the contribution of immune mechanisms to Parkinson's disease pathogenesis is an important challenge, potentially of major therapeutic implications. To further elucidate the involvement of peripheral immune cells, we studied epigenome-wide DNA methylation in isolated populations of CD14+ monocytes, CD19+ B cells, CD4+ T cells, and CD8+ T cells from Parkinson's disease patients and healthy control participants. We included 25 patients with a maximum five years of disease duration and 25 controls, and isolated four immune cell populations from each fresh blood sample. Epigenome-wide DNA methylation profiles were generated from 186 samples using the Illumina MethylationEpic array and association with disease status was tested using linear regression models. We identified six differentially methylated CpGs in CD14+ monocytes and one in CD8 + T cells. Four differentially methylated regions were identified in monocytes, including a region upstream of RAB32, a gene that has been linked to LRRK2. Methylation upstream of RAB32 correlated negatively with mRNA expression, and RAB32 expression was upregulated in Parkinson's disease both in our samples and in summary statistics from a previous study. Our epigenome-wide association study of early Parkinson's disease provides evidence for methylation changes across different peripheral immune cell types, highlighting monocytes and the RAB32 locus. The findings were predominantly cell-type-specific, demonstrating the value of isolating purified cell populations for genomic studies.

XMAP: Cross-population fine-mapping by leveraging genetic diversity and accounting for confounding bias

Fine-mapping prioritizes risk variants identified by genome-wide association studies (GWASs), serving as a critical step to uncover biological mechanisms underlying complex traits. However, several major challenges still remain for existing fine-mapping methods. First, the strong linkage disequilibrium among variants can limit the statistical power and resolution of fine-mapping. Second, it is computationally expensive to simultaneously search for multiple causal variants. Third, the confounding bias hidden in GWAS summary statistics can produce spurious signals. To address these challenges, we develop a statistical method for cross-population fine-mapping (XMAP) by leveraging genetic diversity and accounting for confounding bias. By using cross-population GWAS summary statistics from global biobanks and genomic consortia, we show that XMAP can achieve greater statistical power, better control of false positive rate, and substantially higher computational efficiency for identifying multiple causal signals, compared to existing methods. Importantly, we show that the output of XMAP can be integrated with single-cell datasets, which greatly improves the interpretation of putative causal variants in their cellular context at single-cell resolution.

Beyond the amyloid cascade: An update of Alzheimer's disease pathophysiology

Alzheimer's disease (AD) is a multi-etiology disease. The biological system of AD is associated with multidomain genetic, molecular, cellular, and network brain dysfunctions, interacting with central and peripheral immunity. These dysfunctions have been primarily conceptualized according to the assumption that amyloid deposition in the brain, whether from a stochastic or a genetic accident, is the upstream pathological change. However, the arborescence of AD pathological changes suggests that a single amyloid pathway might be too restrictive or inconsistent with a cascading effect. In this review, we discuss the recent human studies of late-onset AD pathophysiology in an attempt to establish a general updated view focusing on the early stages. Several factors highlight heterogenous multi-cellular pathological changes in AD, which seem to work in a self-amplifying manner with amyloid and tau pathologies. Neuroinflammation has an increasing importance as a major pathological driver, and perhaps as a convergent biological basis of aging, genetic, lifestyle and environmental risk factors.

Neurogenetic mechanisms of risk for ADHD: Examining associations of polygenic scores and brain volumes in a population cohort

BACKGROUND

ADHD polygenic scores (PGSs) have been previously shown to predict ADHD outcomes in several studies. However, ADHD PGSs are typically correlated with ADHD but not necessarily reflective of causal mechanisms. More research is needed to elucidate the neurobiological mechanisms underlying ADHD. We leveraged functional annotation information into an ADHD PGS to (1) improve the prediction performance over a non-annotated ADHD PGS and (2) test whether volumetric variation in brain regions putatively associated with ADHD mediate the association between PGSs and ADHD outcomes.

METHODS

Data were from the Philadelphia Neurodevelopmental Cohort (N = 555). Multiple mediation models were tested to examine the indirect effects of two ADHD PGSs-one using a traditional computation involving clumping and thresholding and another using a functionally annotated approach (i.e., AnnoPred)-on ADHD inattention (IA) and hyperactivity-impulsivity (HI) symptoms, via gray matter volumes in the cingulate gyrus, angular gyrus, caudate, dorsolateral prefrontal cortex (DLPFC), and inferior temporal lobe.

RESULTS

A direct effect was detected between the AnnoPred ADHD PGS and IA symptoms in adolescents. No indirect effects via brain volumes were detected for either IA or HI symptoms. However, both ADHD PGSs were negatively associated with the DLPFC.

CONCLUSIONS

The AnnoPred ADHD PGS was a more developmentally specific predictor of adolescent IA symptoms compared to the traditional ADHD PGS. However, brain volumes did not mediate the effects of either a traditional or AnnoPred ADHD PGS on ADHD symptoms, suggesting that we may still be underpowered in clarifying brain-based biomarkers for ADHD using genetic measures.

A genome-wide association study of frailty identifies significant genetic correlation with neuropsychiatric, cardiovascular, and inflammation pathways

Frailty is an aging-related clinical phenotype defined as a state in which there is an increase in a person's vulnerability for dependency and/or mortality when exposed to a stressor. While underlying mechanisms leading to the occurrence of frailty are complex, the importance of genetic factors has not been fully investigated. We conducted a large-scale genome-wide association study (GWAS) of frailty, as defined by the five criteria (weight loss, exhaustion, physical activity, walking speed, and grip strength) captured in the Fried Frailty Score (FFS), in 386,565 European descent participants enrolled in the UK Biobank (mean age 57 [SD 8] years, 208,481 [54%] females). We identified 37 independent, novel loci associated with the FFS (p < 5 × 10-8), including seven loci without prior described associations with other traits. The variants associated with FFS were significantly enriched in brain tissues as well as aging-related pathways. Our post-GWAS bioinformatic analyses revealed significant genetic correlations between FFS and cardiovascular-, neurological-, and inflammation-related diseases/traits, and subsequent Mendelian Randomization analyses identified causal associations with chronic pain, obesity, diabetes, education-related traits, joint disorders, and depressive/neurological, metabolic, and respiratory diseases. The GWAS signals were replicated in the Health and Retirement Study (HRS, n = 9,720, mean age 73 [SD 7], 5,582 [57%] females), where the polygenic risk score built from UKB GWAS was significantly associated with the FFS in HRS individuals (OR per SD of the score 1.27, 95% CI 1.22-1.31, p = 1.3 × 10-11). These results provide new insight into the biology of frailty by comprehensively evaluating its genetic architecture.

Integrative epigenomic and functional characterization assay based annotation of regulatory activity across diverse human cell types

We introduce ChromActivity, a computational framework for predicting and annotating regulatory activity across the genome through integration of multiple epigenomic maps and various functional characterization datasets. ChromActivity generates genomewide predictions of regulatory activity associated with each functional characterization dataset across many cell types based on available epigenomic data. It then for each cell type produces (1) ChromScoreHMM genome annotations based on the combinatorial and spatial patterns within these predictions and (2) ChromScore tracks of overall predicted regulatory activity. ChromActivity provides a resource for analyzing and interpreting the human regulatory genome across diverse cell types.

graph-GPA 2.0: improving multi-disease genetic analysis with integration of functional annotation data

Genome-wide association studies (GWAS) have successfully identified a large number of genetic variants associated with traits and diseases. However, it still remains challenging to fully understand the functional mechanisms underlying many associated variants. This is especially the case when we are interested in variants shared across multiple phenotypes. To address this challenge, we propose graph-GPA 2.0 (GGPA 2.0), a statistical framework to integrate GWAS datasets for multiple phenotypes and incorporate functional annotations within a unified framework. Our simulation studies showed that incorporating functional annotation data using GGPA 2.0 not only improves the detection of disease-associated variants, but also provides a more accurate estimation of relationships among diseases. Next, we analyzed five autoimmune diseases and five psychiatric disorders with the functional annotations derived from GenoSkyline and GenoSkyline-Plus, along with the prior disease graph generated by biomedical literature mining. For autoimmune diseases, GGPA 2.0 identified enrichment for blood-related epigenetic marks, especially B cells and regulatory T cells, across multiple diseases. Psychiatric disorders were enriched for brain-related epigenetic marks, especially the prefrontal cortex and the inferior temporal lobe for bipolar disorder and schizophrenia, respectively. In addition, the pleiotropy between bipolar disorder and schizophrenia was also detected. Finally, we found that GGPA 2.0 is robust to the use of irrelevant and/or incorrect functional annotations. These results demonstrate that GGPA 2.0 can be a powerful tool to identify genetic variants associated with each phenotype or those shared across multiple phenotypes, while also promoting an understanding of functional mechanisms underlying the associated variants.

A statistical framework to identify cell types whose genetically regulated proportions are associated with complex diseases

Finding disease-relevant tissues and cell types can facilitate the identification and investigation of functional genes and variants. In particular, cell type proportions can serve as potential disease predictive biomarkers. In this manuscript, we introduce a novel statistical framework, cell-type Wide Association Study (cWAS), that integrates genetic data with transcriptomics data to identify cell types whose genetically regulated proportions (GRPs) are disease/trait-associated. On simulated and real GWAS data, cWAS showed good statistical power with newly identified significant GRP associations in disease-associated tissues. More specifically, GRPs of endothelial and myofibroblasts in lung tissue were associated with Idiopathic Pulmonary Fibrosis and Chronic Obstructive Pulmonary Disease, respectively. For breast cancer, the GRP of blood CD8+ T cells was negatively associated with breast cancer (BC) risk as well as survival. Overall, cWAS is a powerful tool to reveal cell types associated with complex diseases mediated by GRPs.

Multi-trait genome-wide association analyses leveraging alcohol use disorder findings identify novel loci for smoking behaviors in the Million Veteran Program

Smoking behaviors and alcohol use disorder (AUD), both moderately heritable traits, commonly co-occur in the general population. Single-trait genome-wide association studies (GWAS) have identified multiple loci for smoking and AUD. However, GWASs that have aimed to identify loci contributing to co-occurring smoking and AUD have used small samples and thus have not been highly informative. Applying multi-trait analysis of GWASs (MTAG), we conducted a joint GWAS of smoking and AUD with data from the Million Veteran Program (N = 318,694). By leveraging GWAS summary statistics for AUD, MTAG identified 21 genome-wide significant (GWS) loci associated with smoking initiation and 17 loci associated with smoking cessation compared to 16 and 8 loci, respectively, identified by single-trait GWAS. The novel loci for smoking behaviors identified by MTAG included those previously associated with psychiatric or substance use traits. Colocalization analysis identified 10 loci shared by AUD and smoking status traits, all of which achieved GWS in MTAG, including variants on SIX3, NCAM1, and near DRD2. Functional annotation of the MTAG variants highlighted biologically important regions on ZBTB20, DRD2, PPP6C, and GCKR that contribute to smoking behaviors. In contrast, MTAG of smoking behaviors and alcohol consumption (AC) did not enhance discovery compared with single-trait GWAS for smoking behaviors. We conclude that using MTAG to augment the power of GWAS enables the identification of novel genetic variants for commonly co-occuring phenotypes, providing new insights into their pleiotropic effects on smoking behavior and AUD.

Investigation of convergent and divergent genetic influences underlying schizophrenia and alcohol use disorder

BACKGROUND

Alcohol use disorder (AUD) and schizophrenia (SCZ) frequently co-occur, and large-scale genome-wide association studies (GWAS) have identified significant genetic correlations between these disorders.

METHODS

We used the largest published GWAS for AUD (total cases = 77 822) and SCZ (total cases = 46 827) to identify genetic variants that influence both disorders (with either the same or opposite direction of effect) and those that are disorder specific.

RESULTS

We identified 55 independent genome-wide significant single nucleotide polymorphisms with the same direction of effect on AUD and SCZ, 8 with robust effects in opposite directions, and 98 with disorder-specific effects. We also found evidence for 12 genes whose pleiotropic associations with AUD and SCZ are consistent with mediation via gene expression in the prefrontal cortex. The genetic covariance between AUD and SCZ was concentrated in genomic regions functional in brain tissues (p = 0.001).

CONCLUSIONS

Our findings provide further evidence that SCZ shares meaningful genetic overlap with AUD.

eQTL Studies: from Bulk Tissues to Single Cells

An expression quantitative trait locus (eQTL) is a chromosomal region where genetic variants are associated with the expression levels of certain genes that can be both nearby or distant. The identifications of eQTLs for different tissues, cell types, and contexts have led to better understanding of the dynamic regulations of gene expressions and implications of functional genes and variants for complex traits and diseases. Although most eQTL studies to date have been performed on data collected from bulk tissues, recent studies have demonstrated the importance of cell-type-specific and context-dependent gene regulations in biological processes and disease mechanisms. In this review, we discuss statistical methods that have been developed to enable the detections of cell-type-specific and context-dependent eQTLs from bulk tissues, purified cell types, and single cells. We also discuss the limitations of the current methods and future research opportunities.

PALM: a powerful and adaptive latent model for prioritizing risk variants with functional annotations

MOTIVATION

The findings from genome-wide association studies (GWASs) have greatly helped us to understand the genetic basis of human complex traits and diseases. Despite the tremendous progress, much effects are still needed to address several major challenges arising in GWAS. First, most GWAS hits are located in the non-coding region of human genome, and thus their biological functions largely remain unknown. Second, due to the polygenicity of human complex traits and diseases, many genetic risk variants with weak or moderate effects have not been identified yet.

RESULTS

To address the above challenges, we propose a powerful and adaptive latent model (PALM) to integrate cell-type/tissue-specific functional annotations with GWAS summary statistics. Unlike existing methods, which are mainly based on linear models, PALM leverages a tree ensemble to adaptively characterize non-linear relationship between functional annotations and the association status of genetic variants. To make PALM scalable to millions of variants and hundreds of functional annotations, we develop a functional gradient-based expectation-maximization algorithm, to fit the tree-based non-linear model in a stable manner. Through comprehensive simulation studies, we show that PALM not only controls false discovery rate well, but also improves statistical power of identifying risk variants. We also apply PALM to integrate summary statistics of 30 GWASs with 127 cell type/tissue-specific functional annotations. The results indicate that PALM can identify more risk variants as well as rank the importance of functional annotations, yielding better interpretation of GWAS results.

AVAILABILITY AND IMPLEMENTATION

The source code is available at https://github.com/YangLabHKUST/PALM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Sparse prediction informed by genetic annotations using the logit normal prior for Bayesian regression tree ensembles

Using high-dimensional genetic variants such as single nucleotide polymorphisms (SNP) to predict complex diseases and traits has important applications in basic research and other clinical settings. For example, predicting gene expression is a necessary first step to identify (putative) causal genes in transcriptome-wide association studies. Due to weak signals, high-dimensionality, and linkage disequilibrium (correlation) among SNPs, building such a prediction model is challenging. However, functional annotations at the SNP level (e.g., as epigenomic data across multiple cell- or tissue-types) are available and could be used to inform predictor importance and aid in outcome prediction. Existing approaches to incorporate annotations have been based mainly on (generalized) linear models. Bayesian additive regression trees (BART), in contrast, is a reliable method to obtain high-quality nonlinear out of sample predictions without overfitting. Unfortunately, the default prior from BART may be too inflexible to handle sparse situations where the number of predictors approaches or surpasses the number of observations. Motivated by our real data application, this article proposes an alternative prior based on the logit normal distribution because it provides a framework that is adaptive to sparsity and can model informative functional annotations. It also provides a framework to incorporate prior information about the between SNP correlations. Computational details for carrying out inference are presented along with the results from a simulation study and a genome-wide prediction analysis of the Alzheimer's Disease Neuroimaging Initiative data.

Liver-Specific Polygenic Risk Score Is Associated with Alzheimer's Disease Diagnosis

BACKGROUND

Our understanding of the pathophysiology underlying Alzheimer's disease (AD) has benefited from genomic analyses, including those that leverage polygenic risk score (PRS) models of disease. The use of functional annotation has been able to improve the power of genomic models.

OBJECTIVE

We sought to leverage genomic functional annotations to build tissue-specific AD PRS models and study their relationship with AD and its biomarkers.

METHODS

We built 13 tissue-specific AD PRS and studied the scores' relationships with AD diagnosis, cerebrospinal fluid (CSF) amyloid, CSF tau, and other CSF biomarkers in two longitudinal cohort studies of AD.

RESULTS

The AD PRS model that was most predictive of AD diagnosis (even without APOE) was the liver AD PRS: n = 1,115; odds ratio = 2.15 (1.67-2.78), p = 3.62×10-9. The liver AD PRS was also statistically significantly associated with cerebrospinal fluid biomarker evidence of amyloid-β (Aβ42:Aβ40 ratio, p = 3.53×10-6) and the phosphorylated tau:amyloid-β ratio (p = 1.45×10-5).

CONCLUSION

These findings provide further evidence of the role of the liver-functional genome in AD and the benefits of incorporating functional annotation into genomic research.

Evaluation of Efficacy and Safety Using Low Dose Radiation Therapy with Alzheimer's Disease: A Protocol for Multicenter Phase II Clinical Trial

BACKGROUND

Alzheimer's disease (AD), the most common cause of dementia, is a neurodegenerative disease resulting from extracellular and intracellular deposits of amyloid-β (Aβ) and neurofibrillary tangles in the brain. Although many clinical studies evaluating pharmacological approaches have been conducted, most have shown disappointing results; thus, innovative strategies other than drugs have been actively attempted.

OBJECTIVE

This study aims to explore low-dose radiation therapy (LDRT) for the treatment of patients with AD based on preclinical evidence, case reports, and a small pilot trial in humans.

METHODS

This study is a phase II, multicenter, prospective, single-blinded, randomized controlled trial that will evaluate the efficacy and safety of LDRT to the whole brain using a linear accelerator in patients with mild AD. Sixty participants will be randomly assigned to three groups: experimental I (24 cGy/6 fractions), experimental II (300 cGy/6 fractions), or sham RT group (0 cGy/6 fractions). During LDRT and follow-up visits after LDRT, possible adverse events will be assessed by the physician's interview and neurological examinations. Furthermore, the effectiveness of LDRT will be measured using neurocognitive function tests and imaging tools at 6 and 12 months after LDRT. We will also monitor the alterations in cytokines, Aβ42/Aβ40 ratio, and tau levels in plasma. Our primary endpoint is the change in cognitive function test scores estimated by the Alzheimer's Disease Assessment Scale-Korea compared to baseline after 6 months of LDRT.

CONCLUSIONS

This study is registered at ClinicalTrials.gov [NCT05635968] and is currently recruiting patients. This study will provide evidence that LDRT is a new treatment strategy for AD.

Microglial efferocytosis: Diving into the Alzheimer's disease gene pool

Genome-wide association studies and functional genomics studies have linked specific cell types, genes, and pathways to Alzheimer's disease (AD) risk. In particular, AD risk alleles primarily affect the abundance or structure, and thus the activity, of genes expressed in macrophages, strongly implicating microglia (the brain-resident macrophages) in the etiology of AD. These genes converge on pathways (endocytosis/phagocytosis, cholesterol metabolism, and immune response) with critical roles in core macrophage functions such as efferocytosis. Here, we review these pathways, highlighting relevant genes identified in the latest AD genetics and genomics studies, and describe how they may contribute to AD pathogenesis. Investigating the functional impact of AD-associated variants and genes in microglia is essential for elucidating disease risk mechanisms and developing effective therapeutic approaches.

Unfolding the genotype-to-phenotype black box of cardiovascular diseases through cross-scale modeling

Complex traits such as cardiovascular diseases (CVD) are the results of complicated processes jointly affected by genetic and environmental factors. Genome-wide association studies (GWAS) identified genetic variants associated with diseases but usually did not reveal the underlying mechanisms. There could be many intermediate steps at epigenetic, transcriptomic, and cellular scales inside the black box of genotype-phenotype associations. In this article, we present a machine-learning-based cross-scale framework GRPath to decipher putative causal paths (pcPaths) from genetic variants to disease phenotypes by integrating multiple omics data. Applying GRPath on CVD, we identified 646 and 549 pcPaths linking putative causal regions, variants, and gene expressions in specific cell types for two types of heart failure, respectively. The findings suggest new understandings of coronary heart disease. Our work promoted the modeling of tissue- and cell type-specific cross-scale regulation to uncover mechanisms behind disease-associated variants, and provided new findings on the molecular mechanisms of CVD.

•

We defined one type of cross-scale genotype-to-phenotype regulation path

•

We designed a framework GRPath to uncover putative regulation paths for diseases

•

GRPath helped uncover molecular mechanisms for two major types of heart failure

Gene-lifestyle interactions in the genomics of human complex traits

The role and biological significance of gene-environment interactions in human traits and diseases remain poorly understood. To address these questions, the CHARGE Gene-Lifestyle Interactions Working Group conducted series of genome-wide interaction studies (GWIS) involving up to 610,475 individuals across four ancestries for three lipids and four blood pressure traits, while accounting for interaction effects with drinking and smoking exposures. Here we used GWIS summary statistics from these studies to decipher potential differences in genetic associations and G×E interactions across phenotype-exposure-ancestry combinations, and to derive insights on the potential mechanistic underlying G×E through in-silico functional analyses. Our analyses show first that interaction effects likely contribute to the commonly reported ancestry-specific genetic effect in complex traits, and second, that some phenotype-exposures pairs are more likely to benefit from a greater detection power when accounting for interactions. It also highlighted modest correlation between marginal and interaction effects, providing material for future methodological development and biological discussions. We also estimated contributions to phenotypic variance, including in particular the genetic heritability conditional on the exposure, and heritability partitioned across a range of functional annotations and cell types. In these analyses, we found multiple instances of potential heterogeneity of functional partitions between exposed and unexposed individuals, providing new evidence for likely exposure-specific genetic pathways. Finally, along this work, we identified potential biases in methods used to jointly meta-analyze genetic and interaction effects. We performed simulations to characterize these limitations and to provide the community with guidelines for future G×E studies.

A general framework for predicting the transcriptomic consequences of non-coding variation and small molecules

Genome wide association studies (GWASs) for complex traits have implicated thousands of genetic loci. Most GWAS-nominated variants lie in noncoding regions, complicating the systematic translation of these findings into functional understanding. Here, we leverage convolutional neural networks to assist in this challenge. Our computational framework, peaBrain, models the transcriptional machinery of a tissue as a two-stage process: first, predicting the mean tissue specific abundance of all genes and second, incorporating the transcriptomic consequences of genotype variation to predict individual abundance on a subject-by-subject basis. We demonstrate that peaBrain accounts for the majority (>50%) of variance observed in mean transcript abundance across most tissues and outperforms regularized linear models in predicting the consequences of individual genotype variation. We highlight the validity of the peaBrain model by calculating non-coding impact scores that correlate with nucleotide evolutionary constraint that are also predictive of disease-associated variation and allele-specific transcription factor binding. We further show how these tissue-specific peaBrain scores can be leveraged to pinpoint functional tissues underlying complex traits, outperforming methods that depend on colocalization of eQTL and GWAS signals. We subsequently: (a) derive continuous dense embeddings of genes for downstream applications; (b) highlight the utility of the model in predicting transcriptomic impact of small molecules and shRNA (on par with in vitro experimental replication of external test sets); (c) explore how peaBrain can be used to model difficult-to-study processes (such as neural induction); and (d) identify putatively functional eQTLs that are missed by high-throughput experimental approaches.

GPA-Tree: statistical approach for functional-annotation-tree-guided prioritization of GWAS results

MOTIVATION

In spite of great success of genome-wide association studies (GWAS), multiple challenges still remain. First, complex traits are often associated with many single nucleotide polymorphisms (SNPs), each with small or moderate effect sizes. Second, our understanding of the functional mechanisms through which genetic variants are associated with complex traits is still limited. To address these challenges, we propose GPA-Tree and it simultaneously implements association mapping and identifies key combinations of functional annotations related to risk-associated SNPs by combining a decision tree algorithm with a hierarchical modeling framework.

RESULTS

First, we implemented simulation studies to evaluate the proposed GPA-Tree method and compared its performance with existing statistical approaches. The results indicate that GPA-Tree outperforms existing statistical approaches in detecting risk-associated SNPs and identifying the true combinations of functional annotations with high accuracy. Second, we applied GPA-Tree to a systemic lupus erythematosus (SLE) GWAS and functional annotation data including GenoSkyline and GenoSkylinePlus. The results from GPA-Tree highlight the dysregulation of blood immune cells, including but not limited to primary B, memory helper T, regulatory T, neutrophils and CD8+ memory T cells in SLE. These results demonstrate that GPA-Tree can be a powerful tool that improves association mapping while facilitating understanding of the underlying genetic architecture of complex traits and potential mechanisms linking risk-associated SNPs with complex traits.

AVAILABILITY AND IMPLEMENTATION

The GPATree software is available at https://dongjunchung.github.io/GPATree/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Leveraging auxiliary data from arbitrary distributions to boost GWAS discovery with Flexible cFDR

Genome-wide association studies (GWAS) have identified thousands of genetic variants that are associated with complex traits. However, a stringent significance threshold is required to identify robust genetic associations. Leveraging relevant auxiliary covariates has the potential to boost statistical power to exceed the significance threshold. Particularly, abundant pleiotropy and the non-random distribution of SNPs across various functional categories suggests that leveraging GWAS test statistics from related traits and/or functional genomic data may boost GWAS discovery. While type 1 error rate control has become standard in GWAS, control of the false discovery rate can be a more powerful approach. The conditional false discovery rate (cFDR) extends the standard FDR framework by conditioning on auxiliary data to call significant associations, but current implementations are restricted to auxiliary data satisfying specific parametric distributions, typically GWAS p-values for related traits. We relax these distributional assumptions, enabling an extension of the cFDR framework that supports auxiliary covariates from arbitrary continuous distributions ("Flexible cFDR"). Our method can be applied iteratively, thereby supporting multi-dimensional covariate data. Through simulations we show that Flexible cFDR increases sensitivity whilst controlling FDR after one or several iterations. We further demonstrate its practical potential through application to an asthma GWAS, leveraging various functional genomic data to find additional genetic associations for asthma, which we validate in the larger, independent, UK Biobank data resource.

VannoPortal: multiscale functional annotation of human genetic variants for interrogating molecular mechanism of traits and diseases

Interpreting the molecular mechanism of genomic variations and their causal relationship with diseases/traits are important and challenging problems in the human genetic study. To provide comprehensive and context-specific variant annotations for biologists and clinicians, here, by systematically integrating over 4TB genomic/epigenomic profiles and frequently-used annotation databases from various biological domains, we develop a variant annotation database, called VannoPortal. In general, the database has following major features: (i) systematically integrates 40 genome-wide variant annotations and prediction scores regarding allele frequency, linkage disequilibrium, evolutionary signature, disease/trait association, tissue/cell type-specific epigenome, base-wise functional prediction, allelic imbalance and pathogenicity; (ii) equips with our recent novel index system and parallel random-sweep searching algorithms for efficient management of backend databases and information extraction; (iii) greatly expands context-dependent variant annotation to incorporate large-scale epigenomic maps and regulatory profiles (such as EpiMap) across over 33 tissue/cell types; (iv) compiles many genome-scale base-wise prediction scores for regulatory/pathogenic variant classification beyond protein-coding region; (v) enables fast retrieval and direct comparison of functional evidence among linked variants using highly interactive web panel in addition to plain table; (vi) introduces many visualization functions for more efficient identification and interpretation of functional variants in single web page. VannoPortal is freely available at http://mulinlab.org/vportal.

Novel characterization of the multivariate genetic architecture of internalizing psychopathology and alcohol use

Genetic correlations suggest that the genetic relationship of alcohol use with internalizing psychopathology depends on the measure of alcohol use. Problematic alcohol use (PAU) is positively genetically correlated with internalizing psychopathology, whereas alcohol consumption ranges from not significantly correlated to moderately negatively correlated with internalizing psychopathology. To explore these different genetic relationships of internalizing psychopathology with alcohol use, we performed a multivariate genome-wide association study of four correlated factors (internalizing psychopathology, PAU, quantity of alcohol consumption, and frequency of alcohol consumption) and then assessed genome-wide and local genetic covariance between these factors. We identified 14 significant regions of local, largely positive, genetic covariance between PAU and internalizing psychopathology and 12 regions of significant local genetic covariance (including both positive and negative genetic covariance) between consumption factors and internalizing psychopathology. Partitioned genetic covariance among functional annotations suggested that brain tissues contribute significantly to positive genetic covariance between internalizing psychopathology and PAU but not to the genetic covariance between internalizing psychopathology and quantity or frequency of alcohol consumption. We hypothesize that genome-wide genetic correlations between alcohol use and psychiatric traits may not capture the more complex shared or divergent genetic architectures at the locus or tissue specific level. This study highlights the complexity of genetic architectures of alcohol use and internalizing psychopathology, and the differing shared genetics of internalizing disorders with PAU compared to consumption.

Integrative analysis of scRNA-seq and GWAS data pinpoints periportal hepatocytes as the relevant liver cell types for blood lipids

Liver, a heterogeneous tissue consisting of various cell types, is known to be relevant for blood lipid traits. By integrating summary statistics from genome-wide association studies (GWAS) of lipid traits and single-cell transcriptome data of the liver, we sought to identify specific cell types in the liver that were most relevant for blood lipid levels. We conducted differential expression analyses for 40 cell types from human and mouse livers in order to construct the cell-type specifically expressed gene sets, which we refer to as construction of the liver cell-type specifically expressed gene sets (CT-SEGS). Under the assumption that CT-SEGS represented specific functions of each cell type, we applied stratified linkage disequilibrium score regression to determine cell types that were most relevant for complex traits and diseases. We first confirmed the validity of this method (of delineating functionally relevant cell types) by identifying the immune cell types as relevant for autoimmune diseases. We further showed that lipid GWAS signals were enriched in the human and mouse periportal hepatocytes. Our results provide important information to facilitate future cellular studies of the metabolic mechanism affecting blood lipid levels.

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

Hematopoietic mosaic chromosomal alterations increase the risk for diverse types of infection

Age is the dominant risk factor for infectious diseases, but the mechanisms linking age to infectious disease risk are incompletely understood. Age-related mosaic chromosomal alterations (mCAs) detected from genotyping of blood-derived DNA, are structural somatic variants indicative of clonal hematopoiesis, and are associated with aberrant leukocyte cell counts, hematological malignancy, and mortality. Here, we show that mCAs predispose to diverse types of infections. We analyzed mCAs from 768,762 individuals without hematological cancer at the time of DNA acquisition across five biobanks. Expanded autosomal mCAs were associated with diverse incident infections (hazard ratio (HR) 1.25; 95% confidence interval (CI) = 1.15-1.36; P = 1.8 × 10-7), including sepsis (HR 2.68; 95% CI = 2.25-3.19; P = 3.1 × 10-28), pneumonia (HR 1.76; 95% CI = 1.53-2.03; P = 2.3 × 10-15), digestive system infections (HR 1.51; 95% CI = 1.32-1.73; P = 2.2 × 10-9) and genitourinary infections (HR 1.25; 95% CI = 1.11-1.41; P = 3.7 × 10-4). A genome-wide association study of expanded mCAs identified 63 loci, which were enriched at transcriptional regulatory sites for immune cells. These results suggest that mCAs are a marker of impaired immunity and confer increased predisposition to infections.

Detecting local genetic correlations with scan statistics

Genetic correlation analysis has quickly gained popularity in the past few years and provided insights into the genetic etiology of numerous complex diseases. However, existing approaches oversimplify the shared genetic architecture between different phenotypes and cannot effectively identify precise genetic regions contributing to the genetic correlation. In this work, we introduce LOGODetect, a powerful and efficient statistical method to identify small genome segments harboring local genetic correlation signals. LOGODetect automatically identifies genetic regions showing consistent associations with multiple phenotypes through a scan statistic approach. It uses summary association statistics from genome-wide association studies (GWAS) as input and is robust to sample overlap between studies. Applied to seven phenotypically distinct but genetically correlated neuropsychiatric traits, we identify 227 non-overlapping genome regions associated with multiple traits, including multiple hub regions showing concordant effects on five or more traits. Our method addresses critical limitations in existing analytic strategies and may have wide applications in post-GWAS analysis.

Heritability Enrichment Implicates Microglia in Parkinson's Disease Pathogenesis

Objective:

Understanding how different parts of the immune system contribute to pathogenesis in Parkinson’s disease is a burning challenge with important therapeutic implications. We studied enrichment of common variant heritability for Parkinson’s disease stratified by immune and brain cell types.

Methods:

We used summary statistics from the most recent meta-analysis of genomewide association studies in Parkinson’s disease and partitioned heritability using linkage disequilibrium score regression, stratified for specific cell types, as defined by open chromatin regions. We also validated enrichment results using a polygenic risk score approach and intersected disease-associated variants with epigenetic data and expression quantitative loci to nominate and explore a putative microglial locus.

Results:

We found significant enrichment of Parkinson’s disease risk heritability in open chromatin regions of microglia and monocytes. Genomic annotations overlapped substantially between these 2 cell types, and only the enrichment signal for microglia remained significant in a joint model. We present evidence suggesting P2RY12, a key microglial gene and target for the antithrombotic agent clopidogrel, as the likely driver of a significant Parkinson’s disease association signal on chromosome 3.

Interpretation:

Our results provide further support for the importance of immune mechanisms in Parkinson’s disease pathogenesis, highlight microglial dysregulation as a contributing etiological factor, and nominate a targetable microglial gene candidate as a pathogenic player. Immune processes can be modulated by therapy, with potentially important clinical implications for future treatment in Parkinson’s disease.

Genome-wide meta-analysis, fine-mapping and integrative prioritization implicate new Alzheimer's disease risk genes

Genome-wide association studies have discovered numerous genomic loci associated with Alzheimer's disease (AD); yet the causal genes and variants are incompletely identified. We performed an updated genome-wide AD meta-analysis, which identified 37 risk loci, including new associations near CCDC6, TSPAN14, NCK2 and SPRED2. Using three SNP-level fine-mapping methods, we identified 21 SNPs with >50% probability each of being causally involved in AD risk and others strongly suggested by functional annotation. We followed this with colocalization analyses across 109 gene expression quantitative trait loci datasets and prioritization of genes by using protein interaction networks and tissue-specific expression. Combining this information into a quantitative score, we found that evidence converged on likely causal genes, including the above four genes, and those at previously discovered AD loci, including BIN1, APH1B, PTK2B, PILRA and CASS4.

A Review of Statistical Methods for Identifying Trait-Relevant Tissues and Cell Types

Genome-wide association studies (GWASs) have identified and replicated many genetic variants that are associated with diseases and disease-related complex traits. However, the biological mechanisms underlying these identified associations remain largely elusive. Exploring the biological mechanisms underlying these associations requires identifying trait-relevant tissues and cell types, as genetic variants likely influence complex traits in a tissue- and cell type-specific manner. Recently, several statistical methods have been developed to integrate genomic data with GWASs for identifying trait-relevant tissues and cell types. These methods often rely on different genomic information and use different statistical models for trait-tissue relevance inference. Here, we present a comprehensive technical review to summarize ten existing methods for trait-tissue relevance inference. These methods make use of different genomic information that include functional annotation information, expression quantitative trait loci information, genetically regulated gene expression information, as well as gene co-expression network information. These methods also use different statistical models that range from linear mixed models to covariance network models. We hope that this review can serve as a useful reference both for methodologists who develop methods and for applied analysts who apply these methods for identifying trait relevant tissues and cell types.

Pinpointing miRNA and genes enrichment over trait-relevant tissue network in Genome-Wide Association Studies

BACKGROUND

Understanding gene regulation is important but difficult. Elucidating tissue-specific gene regulation mechanism is even more challenging and requires gene co-expression network assembled from protein-protein interaction, transcription factor and gene binding, and post-transcriptional regulation (e.g., miRNA targeting) information. The miRNA binding affinity could therefore be changed by SNP(s) located at the 3' untranslated regions (3'UTR) of the target messenger RNA (mRNA) which miRNA(s) interacts with. Genome-wide association study (GWAS) has reported significant numbers of loci hosting SNPs associated with many traits. The goal of this study is to pinpoint GWAS functional variants located in 3'UTRs and elucidate if the genes harboring these variants along with their targeting miRNAs are associated with genetic traits relevant to certain tissues.

METHODS

By applying MIGWAS, CoCoNet, ANNOVAR, and DAVID bioinformatics software and utilizing the gene expression database (e.g. GTEx data) to study GWAS summary statistics for 43 traits from 28 GWAS studies, we have identified a list of miRNAs and targeted genes harboring 3'UTR variants, which could contribute to trait-relevant tissue over miRNA-target gene network.

RESULTS

Our result demonstrated that strong association between traits and tissues exists, and in particular, the Primary Biliary Cirrhosis (PBC) trait has the most significant p-value for all 180 tissues among all 43 traits used for this study. We reported SNPs located in 3'UTR regions of genes (SFMBT2, ZC3HAV1, and UGT3A1) targeted by miRNAs for PBC trait and its tissue association network. After employing Gene Ontology (GO) analysis for PBC trait, we have also identified a very important miRNA targeted gene over miRNA-target gene network, PFKL, which encodes the liver subunit of an enzyme.

CONCLUSIONS

The non-coding variants identified from GWAS studies are casually assumed to be not critical to translated protein product. However, 3' untranslated regions (3'UTRs) of genes harbor variants can often change the binding affinity of targeting miRNAs playing important roles in protein translation degree. Our study has shown that GWAS variants could play important roles on miRNA-target gene networks by contributing the association between traits and tissues. Our analysis expands our knowledge on trait-relevant tissue network and paves way for future human disease studies.

Hematopoietic mosaic chromosomal alterations and risk for infection among 767,891 individuals without blood cancer

Age is the dominant risk factor for infectious diseases, but the mechanisms linking the two are incompletely understood1,2. Age-related mosaic chromosomal alterations (mCAs) detected from blood-derived DNA genotyping, are structural somatic variants associated with aberrant leukocyte cell counts, hematological malignancy, and mortality3-11. Whether mCAs represent independent risk factors for infection is unknown. Here we use genome-wide genotyping of blood DNA to show that mCAs predispose to diverse infectious diseases. We analyzed mCAs from 767,891 individuals without hematological cancer at DNA acquisition across four countries. Expanded mCA (cell fraction >10%) prevalence approached 4% by 60 years of age and was associated with diverse incident infections, including sepsis, pneumonia, and coronavirus disease 2019 (COVID-19) hospitalization. A genome-wide association study of expanded mCAs identified 63 significant loci. Germline genetic alleles associated with expanded mCAs were enriched at transcriptional regulatory sites for immune cells. Our results link mCAs with impaired immunity and predisposition to infections. Furthermore, these findings may also have important implications for the ongoing COVID-19 pandemic, particularly in prioritizing individual preventive strategies and evaluating immunization responses.

Hematopoietic mosaic chromosomal alterations and risk for infection among 767,891 individuals without blood cancer

Age is the dominant risk factor for infectious diseases, but the mechanisms linking the two are incompletely understood1,2. Age-related mosaic chromosomal alterations (mCAs) detected from blood-derived DNA genotyping, are structural somatic variants associated with aberrant leukocyte cell counts, hematological malignancy, and mortality3-11. Whether mCAs represent independent risk factors for infection is unknown. Here we use genome-wide genotyping of blood DNA to show that mCAs predispose to diverse infectious diseases. We analyzed mCAs from 767,891 individuals without hematological cancer at DNA acquisition across four countries. Expanded mCA (cell fraction >10%) prevalence approached 4% by 60 years of age and was associated with diverse incident infections, including sepsis, pneumonia, and coronavirus disease 2019 (COVID-19) hospitalization. A genome-wide association study of expanded mCAs identified 63 significant loci. Germline genetic alleles associated with expanded mCAs were enriched at transcriptional regulatory sites for immune cells. Our results link mCAs with impaired immunity and predisposition to infections. Furthermore, these findings may also have important implications for the ongoing COVID-19 pandemic, particularly in prioritizing individual preventive strategies and evaluating immunization responses.

Leveraging functional annotation to identify genes associated with complex diseases

To increase statistical power to identify genes associated with complex traits, a number of transcriptome-wide association study (TWAS) methods have been proposed using gene expression as a mediating trait linking genetic variations and diseases. These methods first predict expression levels based on inferred expression quantitative trait loci (eQTLs) and then identify expression-mediated genetic effects on diseases by associating phenotypes with predicted expression levels. The success of these methods critically depends on the identification of eQTLs, which may not be functional in the corresponding tissue, due to linkage disequilibrium (LD) and the correlation of gene expression between tissues. Here, we introduce a new method called T-GEN (Transcriptome-mediated identification of disease-associated Genes with Epigenetic aNnotation) to identify disease-associated genes leveraging epigenetic information. Through prioritizing SNPs with tissue-specific epigenetic annotation, T-GEN can better identify SNPs that are both statistically predictive and biologically functional. We found that a significantly higher percentage (an increase of 18.7% to 47.2%) of eQTLs identified by T-GEN are inferred to be functional by ChromHMM and more are deleterious based on their Combined Annotation Dependent Depletion (CADD) scores. Applying T-GEN to 207 complex traits, we were able to identify more trait-associated genes (ranging from 7.7% to 102%) than those from existing methods. Among the identified genes associated with these traits, T-GEN can better identify genes with high (>0.99) pLI scores compared to other methods. When T-GEN was applied to late-onset Alzheimer's disease, we identified 96 genes located at 15 loci, including two novel loci not implicated in previous GWAS. We further replicated 50 genes in an independent GWAS, including one of the two novel loci.

Genome-wide association study of smoking trajectory and meta-analysis of smoking status in 842,000 individuals

Here we report a large genome-wide association study (GWAS) for longitudinal smoking phenotypes in 286,118 individuals from the Million Veteran Program (MVP) where we identified 18 loci for smoking trajectory of current versus never in European Americans, one locus in African Americans, and one in Hispanic Americans. Functional annotations prioritized several dozen genes where significant loci co-localized with either expression quantitative trait loci or chromatin interactions. The smoking trajectories were genetically correlated with 209 complex traits, for 33 of which smoking was either a causal or a consequential factor. We also performed European-ancestry meta-analyses for smoking status in the MVP and GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN) (N_total = 842,717) and identified 99 loci for smoking initiation and 13 loci for smoking cessation. Overall, this large GWAS of longitudinal smoking phenotype in multiple populations, combined with a meta-GWAS for smoking status, adds new insights into the genetic vulnerability for smoking behavior.

Genome-wide association studies (GWASs) for cigarette smoking have identified several hundred loci that account for a small proportion of the overall genetic risk. Here, the authors report a large GWAS for smoking trajectories and meta-analysis for smoking status, finding multiple plausible loci.

Ultrafast and scalable variant annotation and prioritization with big functional genomics data

The advances of large-scale genomics studies have enabled compilation of cell type–specific, genome-wide DNA functional elements at high resolution. With the growing volume of functional annotation data and sequencing variants, existing variant annotation algorithms lack the efficiency and scalability to process big genomic data, particularly when annotating whole-genome sequencing variants against a huge database with billions of genomic features. Here, we develop VarNote to rapidly annotate genome-scale variants in large and complex functional annotation resources. Equipped with a novel index system and a parallel random-sweep searching algorithm, VarNote shows substantial performance improvements (two to three orders of magnitude) over existing algorithms at different scales. It supports both region-based and allele-specific annotations and introduces advanced functions for the flexible extraction of annotations. By integrating massive base-wise and context-dependent annotations in the VarNote framework, we introduce three efficient and accurate pipelines to prioritize the causal regulatory variants for common diseases, Mendelian disorders, and cancers.

LPM: a latent probit model to characterize the relationship among complex traits using summary statistics from multiple GWASs and functional annotations

MOTIVATION

Much effort has been made toward understanding the genetic architecture of complex traits and diseases. In the past decade, fruitful GWAS findings have highlighted the important role of regulatory variants and pervasive pleiotropy. Because of the accumulation of GWAS data on a wide range of phenotypes and high-quality functional annotations in different cell types, it is timely to develop a statistical framework to explore the genetic architecture of human complex traits by integrating rich data resources.

RESULTS

In this study, we propose a unified statistical approach, aiming to characterize relationship among complex traits, and prioritize risk variants by leveraging regulatory information collected in functional annotations. Specifically, we consider a latent probit model (LPM) to integrate summary-level GWAS data and functional annotations. The developed computational framework not only makes LPM scalable to hundreds of annotations and phenotypes but also ensures its statistically guaranteed accuracy. Through comprehensive simulation studies, we evaluated LPM's performance and compared it with related methods. Then, we applied it to analyze 44 GWASs with 9 genic category annotations and 127 cell-type specific functional annotations. The results demonstrate the benefits of LPM and gain insights of genetic architecture of complex traits.

AVAILABILITY AND IMPLEMENTATION

The LPM package, all simulation codes and real datasets in this study are available at https://github.com/mingjingsi/LPM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genetic architecture of Alzheimer's disease

Advances in genetic and genomic technologies over the last thirty years have greatly enhanced our knowledge concerning the genetic architecture of Alzheimer's disease (AD). Several genes including APP, PSEN1, PSEN2, and APOE have been shown to exhibit large effects on disease susceptibility, with the remaining risk loci having much smaller effects on AD risk. Notably, common genetic variants impacting AD are not randomly distributed across the genome. Instead, these variants are enriched within regulatory elements active in human myeloid cells, and to a lesser extent liver cells, implicating these cell and tissue types as critical to disease etiology. Integrative approaches are emerging as highly effective for identifying the specific target genes through which AD risk variants act and will likely yield important insights related to potential therapeutic targets in the coming years. In the future, additional consideration of sex- and ethnicity-specific contributions to risk as well as the contribution of complex gene-gene and gene-environment interactions will likely be necessary to further improve our understanding of AD genetic architecture.

A method for scoring the cell type-specific impacts of noncoding variants in personal genomes

Here we use the expression and accessibility data from a diverse set of cell types to learn a model for the dependence of the accessibility of a regulatory element on its DNA sequence and TF expression. Using GTEx samples with WGS data, we show that the noncoding variants predicted to affect accessibility are more strongly associated with the expression of nearby genes. To interpret a personal genome, we combine the sequence information with context-specific TF expression to prioritize variants and regulatory elements in any genomic region of interest. This approach should be helpful in the study of risk loci previously identified by GWAS. Results from analysis of height and WGS data from the GTEx project support this hypothesis.

A person’s genome typically contains millions of variants which represent the differences between this personal genome and the reference human genome. The interpretation of these variants, i.e., the assessment of their potential impact on a person’s phenotype, is currently of great interest in human genetics and medicine. We have developed a prioritization tool called OpenCausal which takes as inputs 1) a personal genome and 2) a reference context-specific TF expression profile and returns a list of noncoding variants prioritized according to their impact on chromatin accessibility for any given genomic region of interest. We applied OpenCausal to 6,430 samples across 18 tissues derived from the GTEx project and found that the variants prioritized by OpenCausal are highly enriched for eQTLs and caQTLs. We further propose a strategy to integrate the predicted open scores with genome-wide association studies (GWAS) data to prioritize putative causal variants and regulatory elements for a given risk locus (i.e., fine-mapping analysis). As an initial example, we applied this method to a GWAS dataset of human height and found that the prioritized putative variants and elements are correlated with the phenotype (i.e., heights of individuals) better than others.

Statistical Methods in Genome-Wide Association Studies

Since the initial success of genome-wide association studies (GWAS) in 2005, tens of thousands of genetic variants have been identified for hundreds of human diseases and traits. In a GWAS, genotype information at up to millions of genetic markers is collected from up to hundreds of thousands of individuals, together with their phenotype information. Several scientific goals can be accomplished through the analysis of GWAS data, including the identification of variants, genes, and pathways associated with diseases and traits of interest; the inference of the genetic architecture of these traits; and the development of genetic risk prediction models. In this review, we provide an overview of the statistical challenges in achieving these goals and recent progress in statistical methodology to address these challenges.

Low-Dose Ionizing Radiation Modulates Microglia Phenotypes in the Models of Alzheimer's Disease

Alzheimer’s disease (AD) is the most common type of dementia. AD involves major pathologies such as amyloid-β (Aβ) plaques and neurofibrillary tangles in the brain. During the progression of AD, microglia can be polarized from anti-inflammatory M2 to pro-inflammatory M1 phenotype. The activation of triggering receptor expressed on myeloid cells 2 (TREM2) may result in microglia phenotype switching from M1 to M2, which finally attenuated Aβ deposition and memory loss in AD. Low-dose ionizing radiation (LDIR) is known to ameliorate Aβ pathology and cognitive deficits in AD; however, the therapeutic mechanisms of LDIR against AD-related pathology have been little studied. First, we reconfirm that LDIR (two Gy per fraction for five times)-treated six-month 5XFAD mice exhibited (1) the reduction of Aβ deposition, as reflected by thioflavins S staining, and (2) the improvement of cognitive deficits, as revealed by Morris water maze test, compared to sham-exposed 5XFAD mice. To elucidate the mechanisms of LDIR-induced inhibition of Aβ accumulation and memory loss in AD, we examined whether LDIR regulates the microglial phenotype through the examination of levels of M1 and M2 cytokines in 5XFAD mice. In addition, we investigated the direct effects of LDIR on lipopolysaccharide (LPS)-induced production and secretion of M1/M2 cytokines in the BV-2 microglial cells. In the LPS- and LDIR-treated BV-2 cells, the M2 phenotypic marker CD206 was significantly increased, compared with LPS- and sham-treated BV-2 cells. Finally, the effect of LDIR on M2 polarization was confirmed by detection of increased expression of TREM2 in LPS-induced BV2 cells. These results suggest that LDIR directly induced phenotype switching from M1 to M2 in the brain with AD. Taken together, our results indicated that LDIR modulates LPS- and Aβ-induced neuroinflammation by promoting M2 polarization via TREM2 expression, and has beneficial effects in the AD-related pathology such as Aβ deposition and memory loss.

Rare variant association testing in the non-coding genome

The development of next-generation sequencing technologies has opened-up some new possibilities to explore the contribution of genetic variants to human diseases and in particular that of rare variants. Statistical methods have been developed to test for association with rare variants that require the definition of testing units and, in these testing units, the selection of qualifying variants to include in the test. In the coding regions of the genome, testing units are usually the different genes and qualifying variants are selected based on their functional effects on the encoded proteins. Extending these tests to the non-coding regions of the genome is challenging. Testing units are difficult to define as the non-coding genome organisation is still rather unknown. Qualifying variants are difficult to select as the functional impact of non-coding variants on gene expression is hard to predict. These difficulties could explain why very few investigators so far have analysed the non-coding parts of their whole genome sequencing data. These non-coding parts yet represent the vast majority of the genome and some studies suggest that they could play a major role in disease susceptibility. In this review, we discuss recent experimental and statistical developments to gain knowledge on the non-coding genome and how this knowledge could be used to include rare non-coding variants in association tests. We describe the few studies that have considered variants from the non-coding genome in association tests and how they managed to define testing units and select qualifying variants.