Integrative common and rare variant analyses provide insights into the genetic architecture of liver cirrhosis.

We report a multi-ancestry genome-wide association study on liver cirrhosis and its associated endophenotypes, alanine aminotransferase (ALT) and γ-glutamyl transferase. Using data from 12 cohorts, including 18,265 cases with cirrhosis, 1,782,047 controls, up to 1 million individuals with liver function tests and a validation cohort of 21,689 cases and 617,729 controls, we identify and validate 14 risk associations for cirrhosis. Many variants are located near genes involved in hepatic lipid metabolism. One of these, PNPLA3 p.Ile148Met, interacts with alcohol intake, obesity and diabetes on the risk of cirrhosis and hepatocellular carcinoma (HCC). We develop a polygenic risk score that associates with the progression from cirrhosis to HCC. By focusing on prioritized genes from common variant analyses, we find that rare coding variants in GPAM associate with lower ALT, supporting GPAM as a potential target for therapeutic inhibition. In conclusion, this study provides insights into the genetic underpinnings of cirrhosis.

Genetics and epidemiology of mutational barcode-defined clonal hematopoiesis.

Clonal hematopoiesis (CH) arises when a substantial proportion of mature blood cells is derived from a single hematopoietic stem cell lineage. Using whole-genome sequencing of 45,510 Icelandic and 130,709 UK Biobank participants combined with a mutational barcode method, we identified 16,306 people with CH. Prevalence approaches 50% in elderly participants. Smoking demonstrates a dosage-dependent impact on risk of CH. CH associates with several smoking-related diseases. Contrary to published claims, we find no evidence that CH is associated with cardiovascular disease. We provide evidence that CH is driven by genes that are commonly mutated in myeloid neoplasia and implicate several new driver genes. The presence and nature of a driver mutation alters the risk profile for hematological disorders. Nevertheless, most CH cases have no known driver mutations. A CH genome-wide association study identified 25 loci, including 19 not implicated previously in CH. Splicing, protein and expression quantitative trait loci were identified for CD164 and TCL1A.

Large-scale plasma proteomics comparisons through genetics and disease associations.

High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people2, for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity.

Complex effects of sequence variants on lipid levels and coronary artery disease.

Many sequence variants have additive effects on blood lipid levels and, through that, on the risk of coronary artery disease (CAD). We show that variants also have non-additive effects and interact to affect lipid levels as well as affecting variance and correlations. Variance and correlation effects are often signatures of epistasis or gene-environmental interactions. These complex effects can translate into CAD risk. For example, Trp154Ter in FUT2 protects against CAD among subjects with the A1 blood group, whereas it associates with greater risk of CAD in others. His48Arg in ADH1B interacts with alcohol consumption to affect lipid levels and CAD. The effect of variants in TM6SF2 on blood lipids is greatest among those who never eat oily fish but absent from those who often do. This work demonstrates that variants that affect variance of quantitative traits can allow for the discovery of epistasis and interactions of variants with the environment.

Histopathology and levels of proteins in plasma associate with survival after colorectal cancer diagnosis.

BACKGROUND: The TNM system is used to assess prognosis after colorectal cancer (CRC) diagnosis. Other prognostic factors reported include histopathological assessments of the tumour, tumour mutations and proteins in the blood. As some of these factors are strongly correlated, it is important to evaluate the independent effects they may have on survival. METHODS: Tumour samples from 2162 CRC patients were visually assessed for amount of tumour stroma, severity of lymphocytic infiltrate at the tumour margins and the presence of lymphoid follicles. Somatic mutations in the tumour were assessed for 2134 individuals. Pre-surgical levels of 4963 plasma proteins were measured in 128 individuals. The associations between these features and prognosis were inspected by a Cox Proportional Hazards Model (CPH). RESULTS: Levels of stroma, lymphocytic infiltration and presence of lymphoid follicles all associate with prognosis, along with high tumour mutation burden, high microsatellite instability and TP53 and BRAF mutations. The somatic mutations are correlated with the histopathology and none of the somatic mutations associate with survival in a multivariate analysis. Amount of stroma and lymphocytic infiltration associate with local invasion of tumours. Elevated levels of two plasma proteins, CA-125 and PPP1R1A, associate with a worse prognosis. CONCLUSIONS: Tumour stroma and lymphocytic infiltration variables are strongly associated with prognosis of CRC and capture the prognostic effects of tumour mutation status. CA-125 and PPP1R1A may be useful prognostic biomarkers in CRC.

Evaluation of Large-Scale Proteomics for Prediction of Cardiovascular Events.

Importance: Whether protein risk scores derived from a single plasma sample could be useful for risk assessment for atherosclerotic cardiovascular disease (ASCVD), in conjunction with clinical risk factors and polygenic risk scores, is uncertain. Objective: To develop protein risk scores for ASCVD risk prediction and compare them to clinical risk factors and polygenic risk scores in primary and secondary event populations. Design, Setting, and Participants: The primary analysis was a retrospective study of primary events among 13 540 individuals in Iceland (aged 40-75 years) with proteomics data and no history of major ASCVD events at recruitment (study duration, August 23, 2000 until October 26, 2006; follow-up through 2018). We also analyzed a secondary event population from a randomized, double-blind lipid-lowering clinical trial (2013-2016), consisting of individuals with stable ASCVD receiving statin therapy and for whom proteomic data were available for 6791 individuals. Exposures: Protein risk scores (based on 4963 plasma protein levels and developed in a training set in the primary event population); polygenic risk scores for coronary artery disease and stroke; and clinical risk factors that included age, sex, statin use, hypertension treatment, type 2 diabetes, body mass index, and smoking status at the time of plasma sampling. Main Outcomes and Measures: Outcomes were composites of myocardial infarction, stroke, and coronary heart disease death or cardiovascular death. Performance was evaluated using Cox survival models and measures of discrimination and reclassification that accounted for the competing risk of non-ASCVD death. Results: In the primary event population test set (4018 individuals [59.0% women]; 465 events; median follow-up, 15.8 years), the protein risk score had a hazard ratio (HR) of 1.93 per SD (95% CI, 1.75 to 2.13). Addition of protein risk score and polygenic risk scores significantly increased the C index when added to a clinical risk factor model (C index change, 0.022 [95% CI, 0.007 to 0.038]). Addition of the protein risk score alone to a clinical risk factor model also led to a significantly increased C index (difference, 0.014 [95% CI, 0.002 to 0.028]). Among White individuals in the secondary event population (6307 participants; 432 events; median follow-up, 2.2 years), the protein risk score had an HR of 1.62 per SD (95% CI, 1.48 to 1.79) and significantly increased C index when added to a clinical risk factor model (C index change, 0.026 [95% CI, 0.011 to 0.042]). The protein risk score was significantly associated with major adverse cardiovascular events among individuals of African and Asian ancestries in the secondary event population. Conclusions and Relevance: A protein risk score was significantly associated with ASCVD events in primary and secondary event populations. When added to clinical risk factors, the protein risk score and polygenic risk score both provided statistically significant but modest improvement in discrimination.

Sequence variants affecting voice pitch in humans.

The genetic basis of the human vocal system is largely unknown, as are the sequence variants that give rise to individual differences in voice and speech. Here, we couple data on diversity in the sequence of the genome with voice and vowel acoustics in speech recordings from 12,901 Icelanders. We show how voice pitch and vowel acoustics vary across the life span and correlate with anthropometric, physiological, and cognitive traits. We found that voice pitch and vowel acoustics have a heritable component and discovered correlated common variants in ABCC9 that associate with voice pitch. The ABCC9 variants also associate with adrenal gene expression and cardiovascular traits. By showing that voice and vowel acoustics are influenced by genetics, we have taken important steps toward understanding the genetics and evolution of the human vocal system.

Thirty novel sequence variants impacting human intracranial volume.

Intracranial volume, measured through magnetic resonance imaging and/or estimated from head circumference, is heritable and correlates with cognitive traits and several neurological disorders. We performed a genome-wide association study meta-analysis of intracranial volume (n = 79 174) and found 64 associating sequence variants explaining 5.0% of its variance. We used coding variation, transcript and protein levels, to uncover 12 genes likely mediating the effect of these variants, including GLI3 and CDK6 that affect cranial synostosis and microcephaly, respectively. Intracranial volume correlates genetically with volumes of cortical and sub-cortical regions, cognition, learning, neonatal and neurological traits. Parkinson's disease cases have greater and attention deficit hyperactivity disorder cases smaller intracranial volume than controls. Our Mendelian randomization studies indicate that intracranial volume associated variants either increase the risk of Parkinson's disease and decrease the risk of attention deficit hyperactivity disorder and neuroticism or correlate closely with a confounder.

Multiomics study of nonalcoholic fatty liver disease.

Nonalcoholic fatty liver (NAFL) and its sequelae are growing health problems. We performed a genome-wide association study of NAFL, cirrhosis and hepatocellular carcinoma, and integrated the findings with expression and proteomic data. For NAFL, we utilized 9,491 clinical cases and proton density fat fraction extracted from 36,116 liver magnetic resonance images. We identified 18 sequence variants associated with NAFL and 4 with cirrhosis, and found rare, protective, predicted loss-of-function variants in MTARC1 and GPAM, underscoring them as potential drug targets. We leveraged messenger RNA expression, splicing and predicted coding effects to identify 16 putative causal genes, of which many are implicated in lipid metabolism. We analyzed levels of 4,907 plasma proteins in 35,559 Icelanders and 1,459 proteins in 47,151 UK Biobank participants, identifying multiple proteins involved in disease pathogenesis. We show that proteomics can discriminate between NAFL and cirrhosis. The present study provides insights into the development of noninvasive evaluation of NAFL and new therapeutic options.

The sequences of 150,119 genomes in the UK Biobank.

Detailed knowledge of how diversity in the sequence of the human genome affects phenotypic diversity depends on a comprehensive and reliable characterization of both sequences and phenotypic variation. Over the past decade, insights into this relationship have been obtained from whole-exome sequencing or whole-genome sequencing of large cohorts with rich phenotypic data1,2. Here we describe the analysis of whole-genome sequencing of 150,119 individuals from the UK Biobank3. This constitutes a set of high-quality variants, including 585,040,410 single-nucleotide polymorphisms, representing 7.0% of all possible human single-nucleotide polymorphisms, and 58,707,036 indels. This large set of variants allows us to characterize selection based on sequence variation within a population through a depletion rank score of windows along the genome. Depletion rank analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UK Biobank: a large British Irish cohort, a smaller African cohort and a South Asian cohort. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large-scale whole-genome sequencing studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on whole-exome sequencing and/or imputation.

Neuroimaging Findings in Neurodevelopmental Copy Number Variants: Identifying Molecular Pathways to Convergent Phenotypes.

Genomic copy number variants (CNVs) are associated with a high risk of neurodevelopmental disorders. A growing body of genetic studies suggests that these high-risk genetic variants converge in common molecular pathways and that common pathways also exist across clinically distinct disorders, such as schizophrenia and autism spectrum disorder. A key question is how common molecular mechanisms converge into similar clinical outcomes. We review emerging evidence for convergent cognitive and brain phenotypes across distinct CNVs. Multiple CNVs were shown to have similar effects on core sensory, cognitive, and motor traits. Emerging data from multisite neuroimaging studies have provided valuable information on how these CNVs affect brain structure and function. However, most of these studies examined one CNV at a time, making it difficult to fully understand the proportion of shared brain effects. Recent studies have started to combine neuroimaging data from multiple CNV carriers and identified similar brain effects across CNVs. Some early findings also support convergence in CNV animal models. Systems biology, through integration of multilevel data, provides new insights into convergent molecular mechanisms across genetic risk variants (e.g., altered synaptic activity). However, the link between such key molecular mechanisms and convergent psychiatric phenotypes is still unknown. To better understand this link, we need new approaches that integrate human molecular data with neuroimaging, cognitive, and animal model data, while taking into account critical developmental time points. Identifying risk mechanisms across genetic loci can elucidate the pathophysiology of neurodevelopmental disorders and identify new therapeutic targets for cross-disorder applications.

Large-scale integration of the plasma proteome with genetics and disease.

The plasma proteome can help bridge the gap between the genome and diseases. Here we describe genome-wide association studies (GWASs) of plasma protein levels measured with 4,907 aptamers in 35,559 Icelanders. We found 18,084 associations between sequence variants and levels of proteins in plasma (protein quantitative trait loci; pQTL), of which 19% were with rare variants (minor allele frequency (MAF) < 1%). We tested plasma protein levels for association with 373 diseases and other traits and identified 257,490 associations. We integrated pQTL and genetic associations with diseases and other traits and found that 12% of 45,334 lead associations in the GWAS Catalog are with variants in high linkage disequilibrium with pQTL. We identified 938 genes encoding potential drug targets with variants that influence levels of possible biomarkers. Combining proteomics, genomics and transcriptomics, we provide a valuable resource that can be used to improve understanding of disease pathogenesis and to assist with drug discovery and development.

Predicting the probability of death using proteomics.

Predicting all-cause mortality risk is challenging and requires extensive medical data. Recently, large-scale proteomics datasets have proven useful for predicting health-related outcomes. Here, we use measurements of levels of 4,684 plasma proteins in 22,913 Icelanders to develop all-cause mortality predictors both for short- and long-term risk. The participants were 18-101 years old with a mean follow up of 13.7 (sd. 4.7) years. During the study period, 7,061 participants died. Our proposed predictor outperformed, in survival prediction, a predictor based on conventional mortality risk factors. We could identify the 5% at highest risk in a group of 60-80 years old, where 88% died within ten years and 5% at the lowest risk where only 1% died. Furthermore, the predicted risk of death correlates with measures of frailty in an independent dataset. Our results show that the plasma proteome can be used to assess general health and estimate the risk of death.

Analysis of Diffusion Tensor Imaging Data From the UK Biobank Confirms Dosage Effect of 15q11.2 Copy Number Variation on White Matter and Shows Association With Cognition.

BACKGROUND: Copy number variations at the 15q11.2 BP1-BP2 locus are present in 0.5%-1.0% of the population, and the deletion is associated with several neurodevelopmental disorders. Previously, we showed a reciprocal effect of 15q11.2 copy number variation on fractional anisotropy, with widespread increases in deletion carriers. We aim to expand these findings using a larger sample of participants (N = 29,166) and higher resolution imaging and by examining the implications for cognitive performance. METHODS: Diffusion tensor imaging measures from participants with no neurological or psychiatric diagnoses were obtained from the UK Biobank database. We compared 15q11.2 BP1-BP2 deletion (n = 102) and duplication (n = 113) carriers to a large cohort of control individuals with no neuropsychiatric copy number variants (n = 28,951). Additionally, we assessed how changes in white matter mediated the association between carrier status and cognitive performance. RESULTS: Deletion carriers showed increases in fractional anisotropy in the internal capsule and cingulum and decreases in the posterior thalamic radiation compared with both duplication carriers and control subjects (who had intermediate values). Compared with control subjects, deletion carriers had lower scores across cognitive tasks, which were partly influenced by white matter. Reduced fractional anisotropy in the posterior thalamic radiation partially contributed to worse cognitive performance in deletion carriers. CONCLUSIONS: These results, together with our previous findings, provide convergent evidence for an effect of 15q11.2 BP1-BP2 on white matter microstructure, this being more pronounced in deletion carriers. Additionally, changes in white matter were found to partially mediate cognitive ability in deletion carriers, providing a link between white matter changes in 15q11.2 BP1-BP2 carriers and cognitive function.

Sequence Variants in TAAR5 and Other Loci Affect Human Odor Perception and Naming.

Olfactory receptor (OR) genes in humans form a special class characterized by unusually high DNA sequence diversity, which should give rise to differences in perception and behavior. In the largest genome-wide association study to date based on olfactory testing, we investigated odor perception and naming with smell tasks performed by 9,122 Icelanders, with replication in a separate sample of 2,204 individuals. We discovered an association between a low-frequency missense variant in TAAR5 and reduced intensity rating of fish odor containing trimethylamine (p.Ser95Pro, pcombined = 5.6 × 10-15). We demonstrate that TAAR5 genotype affects aversion to fish odor, reflected by linguistic descriptions of the odor and pleasantness ratings. We also discovered common sequence variants in two canonical olfactory receptor loci that associate with increased intensity and naming of licorice odor (trans-anethole: lead variant p.Lys233Asn in OR6C70, pcombined = 8.8 × 10-16 and pcombined = 1.4 × 10-9) and enhanced naming of cinnamon (trans-cinnamaldehyde; intergenic variant rs317787-T, pcombined = 5.0 × 10-17). Together, our results show that TAAR5 genotype variation influences human odor responses and highlight that sequence diversity in canonical OR genes can lead to enhanced olfactory ability, in contrast to the view that greater tolerance for mutations in the human OR repertoire leads to diminished function.

Reciprocal White Matter Changes Associated With Copy Number Variation at 15q11.2 BP1-BP2: A Diffusion Tensor Imaging Study.

BACKGROUND: The 15q11.2 BP1-BP2 cytogenetic region has been associated with learning and motor delays, autism, and schizophrenia. This region includes a gene that codes for the cytoplasmic FMR1 interacting protein 1 (CYFIP1). The CYFIP1 protein is involved in actin cytoskeletal dynamics and interacts with the fragile X mental retardation protein. Absence of fragile X mental retardation protein causes fragile X syndrome. Because abnormal white matter microstructure has been reported in both fragile X syndrome and psychiatric disorders, we looked at the impact of 15q11.2 BP1-BP2 dosage on white matter microstructure. METHODS: Combining a brain-wide voxel-based approach and a regional-based analysis, we analyzed diffusion tensor imaging data from healthy individuals with the deletion (n = 30), healthy individuals with the reciprocal duplication (n = 27), and IQ-matched control subjects with no large copy number variants (n = 19), recruited from a large genotyped population sample. RESULTS: We found global mirror effects (deletion > control > duplication) on fractional anisotropy. The deletion group showed widespread increased fractional anisotropy when compared with duplication. Regional analyses revealed a greater effect size in the posterior limb of the internal capsule and a tendency for decreased fractional anisotropy in duplication. CONCLUSIONS: These results show a reciprocal effect of 15q11.2 BP1-BP2 on white matter microstructure, suggesting that reciprocal chromosomal imbalances may lead to opposite changes in brain structure. Findings in the deletion overlap with previous white matter differences reported in fragile X syndrome patients, suggesting common pathogenic mechanisms derived from disruptions of cytoplasmic CYFIP1-fragile X mental retardation protein complexes. Our data begin to identify specific components of the 15q11.2 BP1-BP2 phenotype and neurobiological mechanisms of potential relevance to the increased risk for disorder.

MAP1B mutations cause intellectual disability and extensive white matter deficit.

Discovery of coding variants in genes that confer risk of neurodevelopmental disorders is an important step towards understanding the pathophysiology of these disorders. Whole-genome sequencing of 31,463 Icelanders uncovers a frameshift variant (E712KfsTer10) in microtubule-associated protein 1B (MAP1B) that associates with ID/low IQ in a large pedigree (genome-wide corrected P = 0.022). Additional stop-gain variants in MAP1B (E1032Ter and R1664Ter) validate the association with ID and IQ. Carriers have 24% less white matter (WM) volume (ß = -2.1SD, P = 5.1 × 10-8), 47% less corpus callosum (CC) volume (ß = -2.4SD, P = 5.5 × 10-10) and lower brain-wide fractional anisotropy (P = 6.7 × 10-4). In summary, we show that loss of MAP1B function affects general cognitive ability through a profound, brain-wide WM deficit with likely disordered or compromised axons.

Neuronal event detection in fMRI time series using iterative deconvolution techniques.

An iterative estimation algorithm for deconvolution of neuronal activity from Blood Oxygen Level Dependent (BOLD) time series data is presented. The algorithm requires knowledge of the hemodynamic impulse response function but does not require knowledge of the stimulation function. The method uses majorization-minimization of a cost function to find an optimal solution to the inverse problem. The cost function includes penalties for the l(1) norm, total variation and negativity. The algorithm is able to identify the occurrence of neuronal activity bursts from BOLD time series accurately. The accuracy of the algorithm was tested in simulations and experimental fMRI data using blocked and event-related designs. The simulations revealed that the algorithm is most sensitive to contrast-to-noise ratio levels and to errors in the assumed hemodynamic model and least sensitive to autocorrelation in the noise. Within normal fMRI conditions, the method is effective for event detection.