Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution

We show that epigenome- and transcriptome-wide association studies (EWAS and TWAS) are prone to significant inflation and bias of test statistics, an unrecognized phenomenon introducing spurious findings if left unaddressed. Neither GWAS-based methodology nor state-of-the-art confounder adjustment methods completely remove bias and inflation. We propose a Bayesian method to control bias and inflation in EWAS and TWAS based on estimation of the empirical null distribution. Using simulations and real data, we demonstrate that our method maximizes power while properly controlling the false positive rate. We illustrate the utility of our method in large-scale EWAS and TWAS meta-analyses of age and smoking

Discovery and fine-mapping of glycaemic and obesity-related trait loci using high-density imputation

Reference panels from the 1000 Genomes (1000G) Project Consortium provide near complete coverage of common and low-frequency genetic variation with minor allele frequency ≥0.5% across European ancestry populations. Within the European Network for Genetic and Genomic Epidemiology (ENGAGE) Consortium, we have undertaken the first large-scale meta-analysis of genome-wide association studies (GWAS), supplemented by 1000G imputation, for four quantitative glycaemic and obesity-related traits, in up to 87,048 individuals of European ancestry. We identified two loci for body mass index (BMI) at genome-wide significance, and two for fasting glucose (FG), none of which has been previously reported in larger meta-analysis efforts to combine GWAS of European ancestry. Through conditional analysis, we also detected multiple distinct signals of association mapping to established loci for waist-hip ratio adjusted for BMI (RSPO3) and FG (GCK and G6PC2). The index variant for one association signal at the G6PC2 locus is a low-frequency coding allele, H177Y, which has recently been demonstrated to have a functional role in glucose regulation. Fine-mapping analyses revealed that the non-coding variants most likely to drive association signals at established and novel loci were enriched for overlap with enhancer elements, which for FG mapped to promoter and transcription factor binding sites in pancreatic islets, in particular. Our study demonstrates that 1000G imputation and genetic fine-mapping of common and low-frequency variant association signals at GWAS loci, integrated with genomic annotation in relevant tissues, can provide insight into the functional and regulatory mechanisms through which their effects on glycaemic and obesity-related traits are mediated

The impact of low-frequency and rare variants on lipid levels

Using a genome-wide screen of 9.6 million genetic variants achieved through 1000 Genomes Project imputation in 62,166 samples, we identify association to lipid traits in 93 loci, including 79 previously identified loci with new lead SNPs and 10 new loci, 15 loci with a low-frequency lead SNP and 10 loci with a missense lead SNP, and 2 loci with an accumulation of rare variants. In six loci, SNPs with established function in lipid genetics (CELSR2, GCKR, LIPC and APOE) or candidate missense mutations with predicted damaging function (CD300LG and TM6SF2) explained the locus associations. The low-frequency variants increased the proportion of variance explained, particularly for low-density lipoprotein cholesterol and total cholesterol. Altogether, our results highlight the impact of low-frequency variants in complex traits and show that imputation offers a cost-effective alternative to resequencing

Genomewide meta-analysis identifies loci associated with IGF-I and IGFBP-3 levels with impact on age-related traits

The growth hormone/insulin-like growth factor (IGF) axis can be manipulated in animal models to promote longevity, and IGF-related proteins including IGF-I and IGF-binding protein-3 (IGFBP-3) have also been implicated in risk of human diseases including cardiovascular diseases, diabetes, and cancer. Throug

Discovery and Fine-Mapping of Glycaemic and Obesity-Related Trait Loci Using High-Density Imputation

Reference panels from the 1000 Genomes (1000G) Project Consortium provide near complete coverage of common and low-frequency genetic variation with minor allele frequency ≥0.5% across European ancestry populations. Within the European Network for Genetic and Genomic Epidemiology (ENGAGE) Consortium, we have undertaken the fi

Meta-analysis of 49 549 individuals imputed with the 1000 Genomes Project reveals an exonic damaging variant in ANGPTL4 determining fasting TG levels

Background So far, more than 170 loci have been associated with circulating lipid levels through genomewide association studies (GWAS). These associations are largely driven by common variants, their function is often not known, and many are likely to be markers for the causal variants. In this study we aimed to identify more new rare and low-frequency functional variants associated with circulating lipid levels. Methods We used the 1000 Genomes Project as a reference panel for the imputations of GWAS data from ~60 000 individuals in the discovery stage and ~90 000 samples in the replication stage. Results Our study resu

WGS-based telomere length analysis in Dutch family trios implicates stronger maternal inheritance and a role for RRM1 gene

Telomere length (TL) regulation is an important factor in ageing, reproduction and cancer development. Genetic, hereditary and environmental factors regulating TL are currently widely investigated, however, their relative contribution to TL variability is still understudied. We have used whole genome sequencing data of 250 family trios from the Genome of the Netherlands project to perform computational measurement of TL and a series of regression and genome-wide association analyses to reveal TL inheritance patterns and associated genetic factors. Our results confirm that TL is a largely heritable trait, primarily with mother’s, and, to a lesser extent, with father’s TL having the strongest influence on the offspring. In this cohort, mother’s, but not father’s age at conception was positively linked to offspring TL. Age-related TL attrition of 40 bp/year had relatively small influence on TL variability. Finally, we have identified TL-associated variations in ribonuclease reductase catalytic subunit M1 (RRM1 gene), which is known to regulate telomere maintenance in yeast. We also highlight the importance of multivariate approach and the limitations of existing tools for the analysis of TL as a polygenic heritable quantitative trait

Schizotypy-related magnetization of cortex in healthy adolescence is colocated with expression of Schizophrenia-related genes

Background: Genetic risk is thought to drive clinical variation on a spectrum of schizophrenia-like traits, but the underlying changes in brain structure that mechanistically link genomic variation to schizotypal experience and behavior are unclear. Methods: We assessed schizotypy using a self-reported questionnaire and measured magnetization transfer as a putative microstructural magnetic resonance imaging marker of intracortical myelination in 68 brain regions in 248 healthy young people (14–25 years of age). We used normative adult brain gene expression data and partial least squares analysis to find the weighted gene expression pattern that was most colocated with the cortical map of schizotypy-related magnetization. Results: Magnetization was significantly correlated with schizotypy in the bilateral posterior cingulate cortex and precuneus (and for disorganized schizotypy, also in medial prefrontal cortex; all false discovery rate–corrected ps < .05), which are regions of the default mode network specialized for social and memory functions. The genes most positively weighted on the whole-genome expression map colocated with schizotypy-related magnetization were enriched for genes that were significantly downregulated in two prior case-control histological studies of brain gene expression in schizophrenia. Conversely, the most negatively weighted genes were enriched for genes that were transcriptionally upregulated in schizophrenia. Positively weighted (downregulated) genes were enriched for neuronal, specifically interneuronal, affiliations and coded a network of proteins comprising a few highly interactive “hubs” such as parvalbumin and calmodulin. Conclusions: Microstructural magnetic resonance imaging maps of intracortical magnetization can be linked to both the behavioral traits of schizotypy and prior histological data on dysregulated gene expression in schizophrenia

Conservative and disruptive modes of adolescent change in human brain functional connectivity

Adolescent changes in human brain function are not entirely understood. Here, we used multiecho functional MRI (fMRI) to measure developmental change in functional connectivity (FC) of resting-state oscillations between pairs of 330 cortical regions and 16 subcortical regions in 298 healthy adolescents scanned 520 times. Participants were aged 14 to 26 y and were scanned on 1 to 3 occasions at least 6 mo apart. We found 2 distinct modes of age-related change in FC: “conservative” and “disruptive.” Conservative development was characteristic of primary cortex, which was strongly connected at 14 y and became even more connected in the period from 14 to 26 y. Disruptive development was characteristic of association cortex and subcortical regions, where connectivity was remodeled: connections that were weak at 14 y became stronger during adolescence, and connections that were strong at 14 y became weaker. These modes of development were quantified using the maturational index (MI), estimated as Spearman’s correlation between edgewise baseline FC (at 14 y, FC14) and adolescent change in FC ( ), at each region. Disruptive systems (with negative MI) were activated by social cognition and autobiographical memory tasks in prior fMRI data and significantly colocated with prior maps of aerobic glycolysis (AG), AG-related gene expression, postnatal cortical surface expansion, and adolescent shrinkage of cortical thickness. The presence of these 2 modes of development was robust to numerous sensitivity analyses. We conclude that human brain organization is disrupted during adolescence by remodeling of FC between association cortical and subcortical areas

A Germline Variant at 8q24 Contributes to Familial Clustering of Prostate Cancer in Men of African Ancestry

Although men of African ancestry have a high risk of prostate cancer (PCa), no genes or mutations have been identified that contribute to familial clustering of PCa in this population. We investigated whether the African ancestry–specific PCa risk variant at 8q24, rs72725854, is enriched in men with a PCa family history in 9052 cases, 143 cases from high-risk families, and 8595 controls of African ancestry. We found the risk allele to be significantly associated with earlier age at diagnosis, more aggressive disease, and enriched in men with a PCa family history (32% of high-risk familial cases carried the variant vs 23% of cases without a family history and 12% of controls). For cases with two or more first-degree relatives with PCa who had at least one family member diagnosed at age <60 yr, the odds ratios for TA heterozygotes and TT homozygotes were 3.92 (95% confidence interval [CI] = 2.13–7.22) and 33.41 (95% CI = 10.86–102.84), respectively. Among men with a PCa family history, the absolute risk by age 60 yr reached 21% (95% CI = 17–25%) for TA heterozygotes and 38% (95% CI = 13–65%) for TT homozygotes. We estimate that in men of African ancestry, rs72725854 accounts for 32% of the total familial risk explained by all known PCa risk variants. Patient summary: We found that rs72725854, an African ancestry–specific risk variant, is more common in men with a family history of prostate cancer and in those diagnosed with prostate cancer at younger ages. Men of African ancestry may benefit from the knowledge of their carrier status for this genetic risk variant to guide decisions about prostate cancer screening. © 2020 The AuthorsThe African ancestry–specific prostate cancer risk variant at 8q24, rs72725854, is enriched in men diagnosed at younger ages and men with a prostate cancer family history. Carriers of this risk allele would benefit from regular and earlier prostate cancer screening