Multiple reference panels of a given tissue or multiple tissues often exist, and multiple regression methods could be used for training gene expression imputation models for TWAS. To leverage expression imputation models (i.e., base models) trained with multiple reference panels, regression methods, and tissues, we develop a Stacked Regression based TWAS (SR-TWAS) tool which can obtain optimal linear combinations of base models for a given validation transcriptomic dataset. Both simulation and real studies showed that SR-TWAS improved power due to increased effective training sample sizes and borrowed strength across multiple regression methods and tissues. Leveraging base models across multiple reference panels, tissues, and regression methods, our studies of Alzheimer's disease (AD) dementia and Parkinson's disease (PD) identified respective 11 independent significant risk genes for AD (supplementary motor area tissue) and 12 independent significant risk genes for PD (substantia nigra tissue), including 6 novels for AD and 6 novels for PD.

BACKGROUND AND OBJECTIVES: Develop imputation models using clinical measures to infer Alzheimer's disease neuropathologic changes (AD-NC) in living adults to identify adults at risk for Alzheimer's disease (AD). METHODS: We used clinical and postmortem data of two prospective cohort studies - Memory and Aging Project (MAP) and Religious Orders Study (ROS). We used generalized linear regression models with Elastic-Net penalty to train imputation models of AD-NC traits ( β -Amyloid, tau tangles, global AD pathology, and NIA-Reagan), in MAP decedents using clinical measures collected at last visit as predictors. ROS cohort was used as an independent validation and test data. We validated these models in ROS decedents and applied the models to baseline clinical data of ROS participants to infer baseline AD-NC traits. Baseline clinical data were collected an average of 8 years before last follow-up. We used Cox proportional hazard models to test if inferred baseline AD-NC traits predicted incident AD dementia (ADD). In addition, two-sample t-tests were used to examine if inferred baseline AD-NC traits predicted adults with high risk of pathologic AD profiled at death. RESULTS: By applying imputation models to clinical measures collected at last visit in ROS to validate the imputation models, we obtained predicted R 2 as 0.188 for β -Amyloid, 0.316 for tau tangles, and 0.262 for global AD pathology. Prediction area under the receiver operating characteristic curve (AUC) for the dichotomous NIA-Reagan was 0.765. All four inferred AD-NC traits at last visit strongly discriminated postmortem NIA-Reagan status (p-values < 10 -28 ). The inferred baseline levels of all four AD-NC traits predicted ADD, with higher accuracies for predicting ADD in Year 3 (AUC ranging in 0.861 - 0.919) versus Year 5 (AUC 0.842 - 0.896), and the highest accuracy was obtained using inferred NIA-Reagan status. The inferred baseline levels of all four AD-NC traits significantly discriminate individuals with postmortem pathologic AD (all p-values < 1.5 × 10 -7 ). CONCLUSIONS: Inferred baseline levels of AD-NC traits derived from clinical measures discriminate adults at risk for ADD and pathologic AD profiled at death. Further studies are needed to determine if repeated measures of inferred AD-NC traits can be used to monitor the accumulation of AD-NC traits during the prolonged course of AD.

Background: Alzheimer's disease (AD) is a progressive disorder without a cure. Develop risk prediction models for detecting presymptomatic AD using non-cognitive measures is necessary to enable early interventions. Objective: Examine if non-cognitive metrics alone can be used to construct risk models to identify adults at risk for AD dementia and cognitive impairment. Methods: Clinical data from older adults without dementia from the Memory and Aging Project (MAP, n = 1,179) and Religious Orders Study (ROS, n = 1,103) were analyzed using Cox proportional hazard models to develop risk prediction models for AD dementia and cognitive impairment. Models using only non-cognitive covariates were compared to models that added cognitive covariates. All models were trained in MAP, tested in ROS, and evaluated by the AUC of ROC curve. Results: Models based on non-cognitive covariates alone achieved AUC (0.800,0.785) for predicting AD dementia (3.5) years from baseline. Including additional cognitive covariates improved AUC to (0.916,0.881). A model with a single covariate of composite cognition score achieved AUC (0.905,0.863). Models based on non-cognitive covariates alone achieved AUC (0.717,0.714) for predicting cognitive impairment (3.5) years from baseline. Including additional cognitive covariates improved AUC to (0.783,0.770). A model with a single covariate of composite cognition score achieved AUC (0.754,0.730). Conclusion: Risk models based on non-cognitive metrics predict both AD dementia and cognitive impairment. However, non-cognitive covariates do not provide incremental predictivity for models that include cognitive metrics in predicting AD dementia, but do in models predicting cognitive impairment. Further improved risk prediction models for cognitive impairment are needed.

Most existing TWAS tools require individual-level eQTL reference data and thus are not applicable to summary-level reference eQTL datasets. The development of TWAS methods that can harness summary-level reference data is valuable to enable TWAS in broader settings and enhance power due to increased reference sample size. Thus, we develop a TWAS framework called OTTERS (Omnibus Transcriptome Test using Expression Reference Summary data) that adapts multiple polygenic risk score (PRS) methods to estimate eQTL weights from summary-level eQTL reference data and conducts an omnibus TWAS. We show that OTTERS is a practical and powerful TWAS tool by both simulations and application studies.

Itch, initiated by the activation of sensory neurons, is associated frequently with dermatological diseases. MrgprA3+ sensory neurons have been identified as one of the major itch-sensing neuronal populations. Mounting evidence has demonstrated that peripheral pathological conditions induce physiological regulation of sensory neurons, which is critical for the maintenance of chronic itch sensation. However, the underlying molecular mechanisms are not clear. Here, we performed RNA sequencing of genetically labeled MrgprA3+ neurons under both naïve and allergic contact dermatitis conditions. Our results revealed the unique molecular signature of itch-sensing neurons and the distinct transcriptional profile changes that result in response to dermatitis. We found enrichment of nine Mrgpr family members and two histamine receptors in MrgprA3+ neurons, suggesting that MrgprA3+ neurons are a direct neuronal target for histamine and Mrgpr agonists. In addition, PTPN6 and PCDH12 were identified as highly selective markers of MrgprA3+ neurons. We also discovered that MrgprA3+ neurons respond to skin dermatitis in a way that is unique from other sensory neurons by regulating a combination of transcriptional factors, ion channels, and key molecules involved in synaptic transmission. These results significantly increase our knowledge of itch transmission and uncover potential targets for combating itch.

Existing methods for integrating functional annotations in genome-wide association studies (GWASs) to fine-map and prioritize potential causal variants are limited to using non-overlapped categorical annotations or limited by the computation burden of modeling genome-wide variants. To overcome these limitations, we propose a scalable Bayesian functional GWAS method to account for multivariate quantitative functional annotations (BFGWAS_QUANT), accompanied by a scalable computation algorithm enabling joint modeling of genome-wide variants. Simulation studies validated the performance of BFGWAS_QUANT for accurately quantifying annotation enrichment and improving GWAS power. Applying BFGWAS_QUANT to study five Alzheimer disease (AD)-related phenotypes using individual-level GWAS data (n = ∼1,000), we found that histone modification annotations have higher enrichment than expression quantitative trait locus (eQTL) annotations for all considered phenotypes, with the highest enrichment in H3K27me3 (polycomb regression). We also found that cis-eQTLs in microglia had higher enrichment than eQTLs of bulk brain frontal cortex tissue for all considered phenotypes. A similar enrichment pattern was also identified using the International Genomics of Alzheimer's Project (IGAP) summary-level GWAS data of AD (n = ∼54,000). The strongest known APOE E4 risk allele was identified for all five phenotypes, and the APOE locus was validated using the IGAP data. BFGWAS_QUANT fine-mapped 32 significant variants from 1,073 genome-wide significant variants in the IGAP data. We also demonstrated that the polygenic risk scores (PRSs) using effect size estimates by BFGWAS_QUANT had a similar prediction accuracy as other methods assuming a sparse causal model. Overall, BFGWAS_QUANT is a useful GWAS tool for quantifying annotation enrichment and prioritizing potential causal variants.

Standard transcriptome-wide association study (TWAS) methods first train gene expression prediction models using reference transcriptomic data and then test the association between the predicted genetically regulated gene expression and phenotype of interest. Most existing TWAS tools require cumbersome preparation of genotype input files and extra coding to enable parallel computation. To improve the efficiency of TWAS tools, we developed Transcriptome-Integrated Genetic Association Resource V2 (TIGAR-V2), which directly reads Variant Call Format (VCF) files, enables parallel computation, and reduces up to 90% of computation cost (mainly due to loading genotype data) compared to the original version. TIGAR-V2 can train gene expression imputation models using either nonparametric Bayesian Dirichlet process regression (DPR) or Elastic-Net (as used by PrediXcan), perform TWASs using either individual-level or summary-level genome-wide association study (GWAS) data, and implement both burden and variance-component statistics for gene-based association tests. We trained gene expression prediction models by DPR for 49 tissues using Genotype-Tissue Expression (GTEx) V8 by TIGAR-V2 and illustrated the usefulness of these Bayesian cis-expression quantitative trait locus (eQTL) weights through TWASs of breast and ovarian cancer utilizing public GWAS summary statistics. We identified 88 and 37 risk genes, respectively, for breast and ovarian cancer, most of which are either known or near previously identified GWAS (∼95%) or TWAS (∼40%) risk genes and three novel independent TWAS risk genes with known functions in carcinogenesis. These findings suggest that TWASs can provide biological insight into the transcriptional regulation of complex diseases. The TIGAR-V2 tool, trained Bayesian cis-eQTL weights, and linkage disequilibrium (LD) information from GTEx V8 are publicly available, providing a useful resource for mapping risk genes of complex diseases.

Cross-species translational approaches to human genomic analyses are lacking. The present study uses an integrative framework to investigate how genes associated with nicotine use in model organisms contribute to the genetic architecture of human tobacco consumption. First, we created a model organism geneset by collecting results from five animal models of nicotine exposure (RNA expression changes in brain) and then tested the relevance of these genes and flanking genetic variation using genetic data from human cigarettes per day (UK BioBank N = 123,844; all European Ancestry). We tested three hypotheses: (1) DNA variation in, or around, the ‘model organism geneset’ will contribute to the heritability to human tobacco consumption, (2) that the model organism genes will be enriched for genes associated with human tobacco consumption, and (3) that a polygenic score based off our model organism geneset will predict tobacco consumption in the AddHealth sample (N = 1667; all European Ancestry). Our results suggested that: (1) model organism genes accounted for ~5–36% of the observed SNP-heritability in human tobacco consumption (enrichment: 1.60–31.45), (2) model organism genes, but not negative control genes, were enriched for the gene-based associations (MAGMA, H-MAGMA, SMultiXcan) for human cigarettes per day, and (3) polygenic scores based on our model organism geneset predicted cigarettes per day in an independent sample. Altogether, these findings highlight the advantages of using multiple species evidence to isolate genetic factors to better understand the etiological complexity of tobacco and other nicotine consumption.

Late-life depression is associated with an increased risk for dementia but we have limited knowledge of the molecular mechanisms underlying this association. Here we investigated whether brain microRNAs, important posttranscriptional regulators of gene expression, contribute to this association. Late-life depressive symptoms were assessed annually in 300 participants of the Religious Orders Study and Rush Memory and Aging Project for a mean of 7 years. Participants underwent annual cognitive testing, clinical assessment of cognitive status, and uniform neuropathologic examination after death. microRNAs were profiled from the prefrontal cortex using NanoString platform in the discovery cohort and small RNA sequencing in the replication cohort. A global microRNA association study of late-life depressive symptoms was performed using linear mixed model adjusting for the potential confounding factors. Four brain microRNAs were associated with late-life depressive symptoms at adjusted p < 0.05: miR-484, miR-26b-5p, miR-30d-5p, and miR-197-3p. Lower expression levels of these miRNAs were associated having greater depressive symptoms. Furthermore, lower levels of miR-484 and miR-197-3p were associated with faster decline of cognition over time. Moreover, lower miR-484 level was associated with higher probability of having Alzheimer's dementia. Importantly, the associations between miR-484 and depressive symptoms and Alzheimer's dementia, respectively, were replicated in an independent cohort. Lastly, the predicted targets of miR-484 were enriched in a brain protein co-expression module involving synaptic transmission and regulation of synaptic plasticity. This study identified four brain microRNAs associated with late-life depressive symptoms assessed longitudinally. In addition, we found a molecular connection between late-life depression and dementia through miR-484.

The transcriptome-wide association studies (TWASs) that test for association between the study trait and the imputed gene expression levels from cis-acting expression quantitative trait loci (cis-eQTL) genotypes have successfully enhanced the discovery of genetic risk loci for complex traits. By using the gene expression imputation models fitted from reference datasets that have both genetic and transcriptomic data, TWASs facilitate gene-based tests with GWAS data while accounting for the reference transcriptomic data. The existing TWAS tools like PrediXcan and FUSION use parametric imputation models that have limitations for modeling the complex genetic architecture of transcriptomic data. Therefore, to improve on this, we employ a nonparametric Bayesian method that was originally proposed for genetic prediction of complex traits, which assumes a data-driven nonparametric prior for cis-eQTL effect sizes. The nonparametric Bayesian method is flexible and general because it includes both of the parametric imputation models used by PrediXcan and FUSION as special cases. Our simulation studies showed that the nonparametric Bayesian model improved both imputation R2 for transcriptomic data and the TWAS power over PrediXcan when ≥1% cis-SNPs co-regulate gene expression and gene expression heritability ≤0.2. In real applications, the nonparametric Bayesian method fitted transcriptomic imputation models for 57.8% more genes over PrediXcan, thus improving the power of follow-up TWASs. We implement both parametric PrediXcan and nonparametric Bayesian methods in a convenient software tool “TIGAR” (Transcriptome-Integrated Genetic Association Resource), which imputes transcriptomic data and performs subsequent TWASs using individual-level or summary-level GWAS data.