Identifying drug targets for neurological and psychiatric disease via genetics and the brain transcriptome.

Discovering drugs that efficiently treat brain diseases has been challenging. Genetic variants that modulate the expression of potential drug targets can be utilized to assess the efficacy of therapeutic interventions. We therefore employed Mendelian Randomization (MR) on gene expression measured in brain tissue to identify drug targets involved in neurological and psychiatric diseases. We conducted a two-sample MR using cis-acting brain-derived expression quantitative trait loci (eQTLs) from the Accelerating Medicines Partnership for Alzheimer's Disease consortium (AMP-AD) and the CommonMind Consortium (CMC) meta-analysis study (n = 1,286) as genetic instruments to predict the effects of 7,137 genes on 12 neurological and psychiatric disorders. We conducted Bayesian colocalization analysis on the top MR findings (using P<6x10-7 as evidence threshold, Bonferroni-corrected for 80,557 MR tests) to confirm sharing of the same causal variants between gene expression and trait in each genomic region. We then intersected the colocalized genes with known monogenic disease genes recorded in Online Mendelian Inheritance in Man (OMIM) and with genes annotated as drug targets in the Open Targets platform to identify promising drug targets. 80 eQTLs showed MR evidence of a causal effect, from which we prioritised 47 genes based on colocalization with the trait. We causally linked the expression of 23 genes with schizophrenia and a single gene each with anorexia, bipolar disorder and major depressive disorder within the psychiatric diseases and 9 genes with Alzheimer's disease, 6 genes with Parkinson's disease, 4 genes with multiple sclerosis and two genes with amyotrophic lateral sclerosis within the neurological diseases we tested. From these we identified five genes (ACE, GPNMB, KCNQ5, RERE and SUOX) as attractive drug targets that may warrant follow-up in functional studies and clinical trials, demonstrating the value of this study design for discovering drug targets in neuropsychiatric diseases.

Amino acid residues in five separate HLA genes can explain most of the known associations between the MHC and primary biliary cholangitis.

Primary Biliary Cholangitis (PBC) is a chronic autoimmune liver disease characterised by progressive destruction of intrahepatic bile ducts. The strongest genetic association is with HLA-DQA1*04:01, but at least three additional independent HLA haplotypes contribute to susceptibility. We used dense single nucleotide polymorphism (SNP) data in 2861 PBC cases and 8514 controls to impute classical HLA alleles and amino acid polymorphisms using state-of-the-art methodologies. We then demonstrated through stepwise regression that association in the HLA region can be largely explained by variation at five separate amino acid positions. Three-dimensional modelling of protein structures and calculation of electrostatic potentials for the implicated HLA alleles/amino acid substitutions demonstrated a correlation between the electrostatic potential of pocket P6 in HLA-DP molecules and the HLA-DPB1 alleles/amino acid substitutions conferring PBC susceptibility/protection, highlighting potential new avenues for future functional investigation.

Genetic association analysis identifies variants associated with disease progression in primary sclerosing cholangitis.

OBJECTIVE: Primary sclerosing cholangitis (PSC) is a genetically complex, inflammatory bile duct disease of largely unknown aetiology often leading to liver transplantation or death. Little is known about the genetic contribution to the severity and progression of PSC. The aim of this study is to identify genetic variants associated with PSC disease progression and development of complications. DESIGN: We collected standardised PSC subphenotypes in a large cohort of 3402 patients with PSC. After quality control, we combined 130 422 single nucleotide polymorphisms of all patients-obtained using the Illumina immunochip-with their disease subphenotypes. Using logistic regression and Cox proportional hazards models, we identified genetic variants associated with binary and time-to-event PSC subphenotypes. RESULTS: We identified genetic variant rs853974 to be associated with liver transplant-free survival (p=6.07×10-9). Kaplan-Meier survival analysis showed a 50.9% (95% CI 41.5% to 59.5%) transplant-free survival for homozygous AA allele carriers of rs853974 compared with 72.8% (95% CI 69.6% to 75.7%) for GG carriers at 10 years after PSC diagnosis. For the candidate gene in the region, RSPO3, we demonstrated expression in key liver-resident effector cells, such as human and murine cholangiocytes and human hepatic stellate cells. CONCLUSION: We present a large international PSC cohort, and report genetic loci associated with PSC disease progression. For liver transplant-free survival, we identified a genome-wide significant signal and demonstrated expression of the candidate gene RSPO3 in key liver-resident effector cells. This warrants further assessments of the role of this potential key PSC modifier gene.

Generation of primary human intestinal T cell transcriptomes reveals differential expression at genetic risk loci for immune-mediated disease.

OBJECTIVE: Genome-wide association studies (GWAS) have identified genetic variants within multiple risk loci as predisposing to intestinal inflammatory diseases, including Crohn's disease, ulcerative colitis and coeliac disease. Most risk variants affect regulation of transcription, but a critical challenge is to identify which genes and which cell types these variants affect. We aimed to characterise whole transcriptomes for each common T lymphocyte subset resident within the gut mucosa, and use these to infer biological insights and highlight candidate genes of interest within GWAS risk loci. DESIGN: We isolated the four major intestinal T cell populations from pinch biopsies from healthy subjects and generated transcriptomes for each. We computationally integrated these transcriptomes with GWAS data from immune-related diseases. RESULTS: Robust, high quality transcriptomic data were generated from 1 ng of RNA from precisely sorted cell subsets. Gene expression patterns clearly differentiated intestinal T cells from counterparts in peripheral blood and revealed distinct signalling pathways for each intestinal T cell subset. Intestinal-specific T cell transcripts were enriched in GWAS risk loci for Crohn's disease, ulcerative colitis and coeliac disease, but also specific extraintestinal immune-mediated diseases, allowing prediction of novel candidate genes. CONCLUSIONS: This is the first report of transcriptomes for minimally manipulated intestinal T lymphocyte subsets in humans. We have demonstrated that careful processing of mucosal biopsies allows the generation of transcriptomes from as few as 1000 highly purified cells with minimal interindividual variation. Bioinformatic integration of transcriptomic data with recent GWAS data identified specific candidate genes and cell types for inflammatory pathologies.

GWAS findings for human iris patterns: associations with variants in genes that influence normal neuronal pattern development.

Human iris patterns are highly variable. The origins of this variation are of interest in the study of iris-related eye diseases and forensics, as well as from an embryological developmental perspective, with regard to their possible relationship to fundamental processes of neurodevelopment. We have performed genome-wide association scans on four iris characteristics (crypt frequency, furrow contractions, presence of peripupillary pigmented ring, and number of nevi) in three Australian samples of European descent. Both the discovery (n = 2121) and replication (n = 499 and 73) samples showed evidence for association between (1) crypt frequency and variants in the axonal guidance gene SEMA3A (p = 6.6 × 10(-11)), (2) furrow contractions and variants within the cytoskeleton gene TRAF3IP1 (p = 2.3 × 10(-12)), and (3) the pigmented ring and variants in the well-known pigmentation gene SLC24A4 (p = 7.6 × 10(-21)). These replicated findings individually accounted for around 1.5%-3% of the variance for these iris characteristics. Because both SEMA3A and TRAFIP1 are implicated in pathways that control neurogenesis, neural migration, and synaptogenesis, we also examined the evidence of enhancement among such genes, finding enrichment for crypts and furrows. These findings suggest that genes involved in normal neuronal pattern development may also influence tissue structures in the human iris.

Whole-exome sequencing in UK Biobank reveals rare genetic architecture for depression.

Nearly two hundred common-variant depression risk loci have been identified by genome-wide association studies (GWAS). However, the impact of rare coding variants on depression remains poorly understood. Here, we present whole-exome sequencing analyses of depression with seven different definitions based on survey, questionnaire, and electronic health records in 320,356 UK Biobank participants. We showed that the burden of rare damaging coding variants in loss-of-function intolerant genes is significantly associated with risk of depression with various definitions. We compared the rare and common genetic architecture across depression definitions by genetic correlation and showed different genetic relationships between definitions across common and rare variants. In addition, we demonstrated that the effects of rare damaging coding variant burden and polygenic risk score on depression risk are additive. The gene set burden analyses revealed overlapping rare genetic variant components with developmental disorder, autism, and schizophrenia. Our study provides insights into the contribution of rare coding variants, separately and in conjunction with common variants, on depression with various definitions and their genetic relationships with neurodevelopmental disorders.

A versatile gene-based test for genome-wide association studies.

We have derived a versatile gene-based test for genome-wide association studies (GWAS). Our approach, called VEGAS (versatile gene-based association study), is applicable to all GWAS designs, including family-based GWAS, meta-analyses of GWAS on the basis of summary data, and DNA-pooling-based GWAS, where existing approaches based on permutation are not possible, as well as singleton data, where they are. The test incorporates information from a full set of markers (or a defined subset) within a gene and accounts for linkage disequilibrium between markers by using simulations from the multivariate normal distribution. We show that for an association study using singletons, our approach produces results equivalent to those obtained via permutation in a fraction of the computation time. We demonstrate proof-of-principle by using the gene-based test to replicate several genes known to be associated on the basis of results from a family-based GWAS for height in 11,536 individuals and a DNA-pooling-based GWAS for melanoma in approximately 1300 cases and controls. Our method has the potential to identify novel associated genes; provide a basis for selecting SNPs for replication; and be directly used in network (pathway) approaches that require per-gene association test statistics. We have implemented the approach in both an easy-to-use web interface, which only requires the uploading of markers with their association p-values, and a separate downloadable application.

The impact of rare protein coding genetic variation on adult cognitive function.

Compelling evidence suggests that human cognitive function is strongly influenced by genetics. Here, we conduct a large-scale exome study to examine whether rare protein-coding variants impact cognitive function in the adult population (n = 485,930). We identify eight genes (ADGRB2, KDM5B, GIGYF1, ANKRD12, SLC8A1, RC3H2, CACNA1A and BCAS3) that are associated with adult cognitive function through rare coding variants with large effects. Rare genetic architecture for cognitive function partially overlaps with that of neurodevelopmental disorders. In the case of KDM5B we show how the genetic dosage of one of these genes may determine the variability of cognitive, behavioral and molecular traits in mice and humans. We further provide evidence that rare and common variants overlap in association signals and contribute additively to cognitive function. Our study introduces the relevance of rare coding variants for cognitive function and unveils high-impact monogenic contributions to how cognitive function is distributed in the normal adult population.

No association of candidate genes with cannabis use in a large sample of Australian twin families.

While there is solid evidence that cannabis use is heritable, attempts to identify genetic influences at the molecular level have yielded mixed results. Here, a large twin family sample (n = 7452) was used to test for association between 10 previously reported candidate genes and lifetime frequency of cannabis use using a gene-based association test. None of the candidate genes reached even nominal significance (P < 0.05). The lack of replication may point to our limited understanding of the neurobiology of cannabis involvement and also to potential publication bias and false-positive findings in previous studies.

The burden of rare protein-truncating genetic variants on human lifespan.

Genetic predisposition has been shown to contribute substantially to the age at which we die. Genome-wide association studies (GWASs) have linked more than 20 loci to phenotypes related to human lifespan1. However, little is known about how lifespan is impacted by gene loss of function. Through whole-exome sequencing of 352,338 UK Biobank participants of European ancestry, we assessed the relevance of protein-truncating variant (PTV) gene burden on individual and parental survival. We identified four exome-wide significant (P < 4.2 × 10-7) human lifespan genes, BRCA1, BRCA2, ATM and TET2. Gene and gene-set, PTV-burden, phenome-wide association studies support known roles of these genes in cancer to impact lifespan at the population level. The TET2 PTV burden was associated with a lifespan through somatic mutation events presumably due to clonal hematopoiesis. The overlap between PTV burden and common variant-based lifespan GWASs was modest, underscoring the value of exome sequencing in well-powered biobank cohorts to complement GWASs for identifying genes underlying complex traits.

Pairwise effects between lipid GWAS genes modulate lipid plasma levels and cellular uptake.

Complex traits are characterized by multiple genes and variants acting simultaneously on a phenotype. However, studying the contribution of individual pairs of genes to complex traits has been challenging since human genetics necessitates very large population sizes, while findings from model systems do not always translate to humans. Here, we combine genetics with combinatorial RNAi (coRNAi) to systematically test for pairwise additive effects (AEs) and genetic interactions (GIs) between 30 lipid genome-wide association studies (GWAS) genes. Gene-based burden tests from 240,970 exomes show that in carriers with truncating mutations in both, APOB and either PCSK9 or LPL ("human double knock-outs") plasma lipid levels change additively. Genetics and coRNAi identify overlapping AEs for 12 additional gene pairs. Overlapping GIs are observed for TOMM40/APOE with SORT1 and NCAN. Our study identifies distinct gene pairs that modulate plasma and cellular lipid levels primarily via AEs and nominates putative drug target pairs for improved lipid-lowering combination therapies.

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.

Genome-wide association study of height and body mass index in Australian twin families.

Human height and body mass index are influenced by a large number of genes, each with small effects, along with environment. To identify common genetic variants associated with these traits, we performed genome-wide association studies in 11,536 individuals composed of Australian twins, family members, and unrelated individuals at approximately 550,000 genotyped SNPs. We identified a single genome-wide significant variant for height (P value=1.06x10(-9)) located in HHIP, a well-replicated height-associated gene. Suggestive levels of association were found for other known genes associated with height (P values<1x10(-6)): ADAMTSL3, EFEMP1, GPR126, and HMGA2; and BMI (P values<1x10(-4)): FTO and MC4R. Together, these variants explain less than 2% of total phenotypic variation for height and 0.5% for BMI.

Phenome-wide Mendelian randomization mapping the influence of the plasma proteome on complex diseases.

The human proteome is a major source of therapeutic targets. Recent genetic association analyses of the plasma proteome enable systematic evaluation of the causal consequences of variation in plasma protein levels. Here we estimated the effects of 1,002 proteins on 225 phenotypes using two-sample Mendelian randomization (MR) and colocalization. Of 413 associations supported by evidence from MR, 130 (31.5%) were not supported by results of colocalization analyses, suggesting that genetic confounding due to linkage disequilibrium is widespread in naïve phenome-wide association studies of proteins. Combining MR and colocalization evidence in cis-only analyses, we identified 111 putatively causal effects between 65 proteins and 52 disease-related phenotypes ( https://www.epigraphdb.org/pqtl/ ). Evaluation of data from historic drug development programs showed that target-indication pairs with MR and colocalization support were more likely to be approved, evidencing the value of this approach in identifying and prioritizing potential therapeutic targets.

Association Between Common Variants in RBFOX1, an RNA-Binding Protein, and Brain Amyloidosis in Early and Preclinical Alzheimer Disease.

Importance: Genetic studies of Alzheimer disease have focused on the clinical or pathologic diagnosis as the primary outcome, but little is known about the genetic basis of the preclinical phase of the disease. Objective: To examine the underlying genetic basis for brain amyloidosis in the preclinical phase of Alzheimer disease. Design, Setting, and Participants: In the first stage of this genetic association study, a meta-analysis was conducted using genetic and imaging data acquired from 6 multicenter cohort studies of healthy older individuals between 1994 and 2019: the Anti-Amyloid Treatment in Asymptomatic Alzheimer Disease Study, the Berkeley Aging Cohort Study, the Wisconsin Registry for Alzheimer's Prevention, the Biomarkers of Cognitive Decline Among Normal Individuals cohort, the Baltimore Longitudinal Study of Aging, and the Alzheimer Disease Neuroimaging Initiative, which included Alzheimer disease and mild cognitive impairment. The second stage was designed to validate genetic observations using pathologic and clinical data from the Religious Orders Study and Rush Memory and Aging Project. Participants older than 50 years with amyloid positron emission tomographic (PET) imaging data and DNA from the 6 cohorts were included. The largest cohort, the Anti-Amyloid Treatment in Asymptomatic Alzheimer Disease Study (n = 3154), was the PET screening cohort used for a secondary prevention trial designed to slow cognitive decline associated with brain amyloidosis. Six smaller, longitudinal cohort studies (n = 1160) provided additional amyloid PET imaging data with existing genetic data. The present study was conducted from March 29, 2019, to February 19, 2020. Main Outcomes and Measures: A genome-wide association study of PET imaging amyloid levels. Results: From the 4314 analyzed participants (age, 52-96 years; 2478 participants [57%] were women), a novel locus for amyloidosis was noted within RBFOX1 (ß = 0.61, P = 3 × 10-9) in addition to APOE. The RBFOX1 protein localized around plaques, and reduced expression of RBFOX1 was correlated with higher amyloid-ß burden (ß = -0.008, P = .002) and worse cognition (ß = 0.007, P = .006) during life in the Religious Orders Study and Rush Memory and Aging Project cohort. Conclusions and Relevance: RBFOX1 encodes a neuronal RNA-binding protein known to be expressed in neuronal tissues and may play a role in neuronal development. The findings of this study suggest that RBFOX1 is a novel locus that may be involved in the pathogenesis of Alzheimer disease.

Case-control association mapping by proxy using family history of disease.

Collecting cases for case-control genetic association studies can be time-consuming and expensive. In some situations (such as studies of late-onset or rapidly lethal diseases), it may be more practical to identify family members of cases. In randomly ascertained cohorts, replacing cases with their first-degree relatives enables studies of diseases that are absent (or nearly absent) in the cohort. We refer to this approach as genome-wide association study by proxy (GWAX) and apply it to 12 common diseases in 116,196 individuals from the UK Biobank. Meta-analysis with published genome-wide association study summary statistics replicated established risk loci and yielded four newly associated loci for Alzheimer's disease, eight for coronary artery disease and five for type 2 diabetes. In addition to informing disease biology, our results demonstrate the utility of association mapping without directly observing cases. We anticipate that GWAX will prove useful in future genetic studies of complex traits in large population cohorts.

Genome-wide association study of primary sclerosing cholangitis identifies new risk loci and quantifies the genetic relationship with inflammatory bowel disease.

Primary sclerosing cholangitis (PSC) is a rare progressive disorder leading to bile duct destruction; ∼75% of patients have comorbid inflammatory bowel disease (IBD). We undertook the largest genome-wide association study of PSC (4,796 cases and 19,955 population controls) and identified four new genome-wide significant loci. The most associated SNP at one locus affects splicing and expression of UBASH3A, with the protective allele (C) predicted to cause nonstop-mediated mRNA decay and lower expression of UBASH3A. Further analyses based on common variants suggested that the genome-wide genetic correlation (rG) between PSC and ulcerative colitis (UC) (rG = 0.29) was significantly greater than that between PSC and Crohn's disease (CD) (rG = 0.04) (P = 2.55 × 10-15). UC and CD were genetically more similar to each other (rG = 0.56) than either was to PSC (P < 1.0 × 10-15). Our study represents a substantial advance in understanding of the genetics of PSC.

Detection and interpretation of shared genetic influences on 42 human traits.

We performed a scan for genetic variants associated with multiple phenotypes by comparing large genome-wide association studies (GWAS) of 42 traits or diseases. We identified 341 loci (at a false discovery rate of 10%) associated with multiple traits. Several loci are associated with multiple phenotypes; for example, a nonsynonymous variant in the zinc transporter SLC39A8 influences seven of the traits, including risk of schizophrenia (rs13107325: log-transformed odds ratio (log OR) = 0.15, P = 2 × 10(-12)) and Parkinson disease (log OR = -0.15, P = 1.6 × 10(-7)), among others. Second, we used these loci to identify traits that have multiple genetic causes in common. For example, variants associated with increased risk of schizophrenia also tended to be associated with increased risk of inflammatory bowel disease. Finally, we developed a method to identify pairs of traits that show evidence of a causal relationship. For example, we show evidence that increased body mass index causally increases triglyceride levels.

Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses.

Very few genetic variants have been associated with depression and neuroticism, likely because of limitations on sample size in previous studies. Subjective well-being, a phenotype that is genetically correlated with both of these traits, has not yet been studied with genome-wide data. We conducted genome-wide association studies of three phenotypes: subjective well-being (n = 298,420), depressive symptoms (n = 161,460), and neuroticism (n = 170,911). We identify 3 variants associated with subjective well-being, 2 variants associated with depressive symptoms, and 11 variants associated with neuroticism, including 2 inversion polymorphisms. The two loci associated with depressive symptoms replicate in an independent depression sample. Joint analyses that exploit the high genetic correlations between the phenotypes (|ρ^| ≈ 0.8) strengthen the overall credibility of the findings and allow us to identify additional variants. Across our phenotypes, loci regulating expression in central nervous system and adrenal or pancreas tissues are strongly enriched for association.