While the genetics of MS risk susceptibility are well-described, and recent progress has been made on the genetics of disease severity, the genetics of disease progression remain elusive. We therefore investigated the genetic determinants of MS progression on longitudinal brain MRI: change in brain volume (BV) and change in T2 lesion volume (T2LV), reflecting progressive tissue loss and increasing disease burden, respectively. We performed genome-wide association studies of change in BV (N = 3401) and change in T2LV (N = 3513) across six randomized clinical trials from Biogen and Roche/Genentech: ADVANCE, ASCEND, DECIDE, OPERA I & II, and ORATORIO. Analyses were adjusted for randomized treatment arm, age, sex, and ancestry. Results were pooled in a meta-analysis, and were evaluated for enrichment of MS risk variants. Variant colocalization and cell-specific expression analyses were performed using published cohorts. The strongest peaks were in PTPRD (rs77321193-C/A, p = 3.9 × 10–7) for BV change, and NEDD4L (rs11398377-GC/G, p = 9.3 × 10–8) for T2LV change. Evidence of colocalization was observed for NEDD4L, and both genes showed increased expression in neuronal and/or glial populations. No association between MS risk variants and MRI outcomes was observed. In this unique, precompetitive industry partnership, we report putative regions of interest in the neurodevelopmental gene PTPRD, and the ubiquitin ligase gene NEDD4L. These findings are distinct from known MS risk genetics, indicating an added role for genetic progression analyses and informing drug discovery.

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored 1 , 2 . FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10 –11 ) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

With advances in NGS technologies, transcriptional profiling of human tissue across many diseases is becoming more routine, leading to the generation of petabytes of data deposited in public repositories. There is a need for bench scientists with little computational expertise to be able to access and mine this data to understand disease pathology, identify robust biomarkers of disease and the effect of interventions (in vivo or in vitro). To this end we release an open source analytics and visualization platform for expression data called OmicsView, http://omicsview.org. This platform comes preloaded with 1000 s of samples across many disease areas and normal tissue, including the GTEx database, all processed with a harmonized pipeline. We demonstrate the power and ease-of-use of the platform by means of a Crohn's disease data mining exercise where we can quickly uncover disease pathology and identify strong biomarkers of disease and response to treatment.

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 lifespan 1 . 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.

ABSTRACT Population isolates such as Finland provide benefits in genetic studies because the allelic spectrum of damaging alleles in any gene is often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%), which survived the founding bottleneck, as opposed to being distributed over a much larger number of ultra--rare variants. While this advantage is well-- established in Mendelian genetics, its value in common disease genetics has been less explored. FinnGen aims to study the genome and national health register data of 500,000 Finns, already reaching 224,737 genotyped and phenotyped participants. Given the relatively high median age of participants (63 years) and dominance of hospital-based recruitment, FinnGen is enriched for many disease endpoints often underrepresented in population-based studies (e.g., rarer immune-mediated diseases and late onset degenerative and ophthalmologic endpoints). We report here a genome-wide association study (GWAS) of 1,932 clinical endpoints defined from nationwide health registries. We identify genome--wide significant associations at 2,491 independent loci. Among these, finemapping implicates 148 putatively causal coding variants associated with 202 endpoints, 104 with low allele frequency (AF<10%) of which 62 were over two-fold enriched in Finland. We studied a benchmark set of 15 diseases that had previously been investigated in large genome-wide association studies. FinnGen discovery analyses were meta-analysed in Estonian and UK biobanks. We identify 30 novel associations, primarily low-frequency variants strongly enriched, in or specific to, the Finnish population and Uralic language family neighbors in Estonia and Russia. These findings demonstrate the power of bottlenecked populations to find unique entry points into the biology of common diseases through low-frequency, high impact variants. Such high impact variants have a potential to contribute to medical translation including drug discovery.

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 association studies (GWAS) have demonstrated the contribution of common genetic variation to cognitive function. However, little is known about the impact of rare variants on cognitive function at the population level. Here, we leveraged exome sequencing from UK Biobank to investigate the impact of rare variants on cognition.

Understanding the genetic basis of Parkinson’s disease （PD） has dramatically progressed since the identification of the first genetic mutation in alpha-synuclein, to become a highly complex model of oligogenic disease （Lubbe et al., 2016; Keogh et al., 2018）. There are more than 90 loci in the genome significantly associated with PD （Nalls et al., 2019）, but more interestingly, most of them are with incomplete penetrance and with only small to moderate risk levels,

SUMMARY Genetic interactions (GIs), the joint impact of different genes or variants on a phenotype, are foundational to the genetic architecture of complex traits. However, identifying GIs through human genetics is challenging since it necessitates very large population sizes, while findings from model systems not always translate to humans. Here, we combined exome-sequencing and genotyping in the UK Biobank with combinatorial RNA-interference (coRNAi) screening to systematically test for pairwise GIs between 30 lipid GWAS genes. Gene-based protein-truncating variant (PTV) burden analyses from 240,970 exomes revealed additive GIs for APOB with PCSK9 and LPL , respectively. Both, genetics and coRNAi identified additive GIs for 12 additional gene pairs. Overlapping non-additive GIs were detected only for TOMM40 at the APOE locus with SORT1 and NCAN . Our study identifies distinct gene pairs that modulate both, plasma and cellular lipid levels via additive and non-additive effects and nominates drug target pairs for improved lipid-lowering combination therapies.

Identifying the biological correlates of late life cognitive function is important if we are to ascertain biomarkers for, and develop treatments to help reduce, age-related cognitive decline. This study investigated the associations between plasma levels of 91 neurology-related proteins (Olink® Proteomics) and general fluid cognitive ability in the Lothian Birth Cohort 1936 (LBC1936, N=798), the Lothian Birth Cohort 1921 (LBC1921, N=165), and the INTERVAL BioResource, (N=4,451). In LBC1936, we also examined mediation of protein-cognitive ability associations by MRI-derived indices of brain structure. In the LBC1936, 22 of the proteins and the first principal component (PC) created from a PC analysis of the 91 proteins, were associated with general fluid cognitive ability (β between −0.11 and −0.17, p<0.0029). Total brain volume partially mediated the association between 10 of these proteins and general fluid cognitive ability. Effect sizes for the 22 proteins, although smaller, were all in the same direction as in LBC1936 in an age-matched subsample of INTERVAL. Similar effect sizes were found for the majority of these 22 proteins in the older LBC1921. The associations were not replicated in a younger subset of INTERVAL. In conclusion, we identified plasma levels of a number of neurology-related proteins that were associated with general fluid cognitive ability in later life, some of which were mediated by brain volume.

Background: GBA variants are the most common genetic risk factors for Parkinson's disease (PD) world-wide, and can be found in up to 20% of Ashkenazi PD patients. The E326K variant, which is not considered a Gaucher's disease causing mutation, was recently shown to increase the risk for PD. Since E326K is a common variant among Europeans, Finnish and Ashkenazi (2.4, 8.6 and 1.2% carrier rate, respectively), we aimed to refine its involvement in PD. Methods: 1200 consecutively recruited PD patients of a full Ashkenazi origin were genotyped for 10 GBA variants, the LRRK2-G2019S and the SMPD1-L302P. Alleles' frequencies were compared to controls, composed of 378 elderly healthy individuals and the non-neuro gnomAD Ashkenazi database. Odds-Ratio (OR) and age-at-motor-symptom-onset (AAO) were also calculated for all genotypes. Results: All allelic variations tested had significant allelic ORs, demonstrating a wide range (1.86-12.84). The lowest allelic OR was observed for E326K (p = .013). Forty-five patients (of 1200, 3.75%) had at least two mutations (of the 12 tested), compared to 2 (0.53%) among 378 controls (p = .0013). Of the E326K carrier patients, 37% (10/27) carried additional mutations and the genotypic OR for individuals who carried only the E326K variant was 1.07. It did not reach statistical significance even when simulating the expected carrier frequency of E326K in 100,000 Ashkenazi controls (p = .39). In addition, an additive effect was demonstrated for risk in carriers of two mutations, the LRRK2-G2019S and a mild-GBA mutation (N370S or R496H), compared to carriers of only one mutation in one of these genes (simulated OR 11.79 compared to 7.58 and 2.49, respectively). An additive effect was also suggested for earlier AAO (5.0 years earlier than in non-carriers, compared to 3.1 and 2.2 years, respectively). Conclusions: Compared to previous studies, we demonstrate here a higher frequency of PD patients that carry two mutations. The GBA-E326K is more likely to affect PD risk when accompanied by another mutation, and an additive effect on risk and earlier AAO was proposed for carriers of LRRK2/mild-GBA double mutations. Altogether, these data support an oligogenic approach to PD genetics.

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P< 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.

BACKGROUND:Individuals with an extremely rare inherited condition, termed Congenital Insensitivity to Pain (CIP), do not feel pain in response to noxious stimuli. Variants in SCN9A, encoding the transmembrane voltage-gated sodium channel Nav1.7, have previously been reported in subjects with CIP accompanied by anosmia, which are typically transmitted in a recessive pattern. Functional characterisations of some of these SCN9A mutations show that they result in complete loss-of-function of Nav1.7.METHODS:In a consanguineous family we performed whole exome sequencing of three members who have a diagnosis of CIP and one unaffected family member. The functional effects of the segregating variant in SCN9A were determined using patch clamp electrophysiology in human embryonic kidney (HEK) 293 cells transfected with the variant.RESULTS:We found that each CIP subject was homozygous for a putatively nonsense variant, R1488*, in SCN9A. This variant was reported elsewhere in a subject with CIP, though the functional effect was not determined. Using electrophysiology, we confirm that this variant results in a complete loss-of-function of Nav1.7.CONCLUSIONS:We confirm through electrophysiological analysis that this R1488* variant in SCN9A results in complete loss-of-function of Nav1.7, which is consistent with reports on other variants in this gene in subjects with CIP.

Background: Cytokines are critical to human disease and are attractive therapeutic targets given their widespread influence on gene regulation and transcription. Defining the downstream regulatory mechanisms influenced by cytokines is central to defining drug and disease mechanisms. One promising strategy is to use interactions between expression quantitative trait loci (eQTLs) and cytokine levels to define target genes and mechanisms. Results: In a clinical trial for anti-IL-6 in patients with systemic lupus erythematosus, we measure interferon (IFN) status, anti-IL-6 drug exposure, and whole blood genome-wide gene expression at three time points. We show that repeat transcriptomic measurements increases the number of cis eQTLs identified compared to using a single time point We observe a statistically significant enrichment of in vivo eQTL interactions with IFN status and anti-IL-6 drug exposure and find many novel interactions that have not been previously described. Finally, we find transcription factor binding motifs interrupted by eQTL interaction SNPs, which point to key regulatory mediators of these environmental stimuli and therefore potential therapeutic targets for autoimmune diseases. In particular, genes with IFN interactions are enriched for ISRE binding site motifs, while those with anti-IL-6 interactions are enriched for IRF4 motifs. Conclusions: This study highlights the potential to exploit clinical trial data to discover in vivo eQTL interactions with therapeutically relevant environmental variables.

To test imputed HLA alleles, AA and peptide binding grooves in 2,957 MS cases and 11,828 controls of European (EUR) ancestry (N=14,785), 227 cases and 908 controls of Latino (LAT) ancestry (N=1,135), and 172 cases and 1,376 controls of African (AFR) ancestry (N=1,548) (total N=17,468).

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We hypothesised that these "resistant smokers" may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the "resistant smokers" and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34 × 10(-4)) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.