1. Genetics and Genomics

The missing link between genetic association and regulatory function

  1. Noah J Connally  Is a corresponding author
  2. Sumaiya Nazeen
  3. Daniel Lee
  4. Huwenbo Shi
  5. John Stamatoyannopoulos
  6. Sung Chun
  7. Chris Cotsapas  Is a corresponding author
  8. Christopher A Cassa  Is a corresponding author
  9. Shamil R Sunyaev  Is a corresponding author
  1. Department of Biomedical Informatics, Harvard Medical School, United States
  2. Brigham and Women’s Hospital, Division of Genetics, Harvard Medical School, United States
  3. Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, United States
  4. Brigham and Women’s Hospital, Department of Neurology, Harvard Medical School, United States
  5. Department of Epidemiology, Harvard T.H. Chan School of Public Health, United States
  6. Altius Institute, United States
  7. Division of Pulmonary Medicine, Boston Children’s Hospital, United States
  8. Department of Neurology, Yale Medical School, United States
  9. Department of Genetics, Yale Medical School, United States
https://doi.org/10.7554/eLife.74970
Share this article
Cite this article
  1. Noah J Connally
  2. Sumaiya Nazeen
  3. Daniel Lee
  4. Huwenbo Shi
  5. John Stamatoyannopoulos
  6. Sung Chun
  7. Chris Cotsapas
  8. Christopher A Cassa
  9. Shamil R Sunyaev
(2022)
The missing link between genetic association and regulatory function
eLife 11:e74970.
https://doi.org/10.7554/eLife.74970
  2. Download BibTeX
  3. Download .RIS

Abstract

The genetic basis of most traits is highly polygenic and dominated by non-coding alleles. It is widely assumed that such alleles exert small regulatory effects on the expression of cis-linked genes. However, despite the availability of gene expression and epigenomic datasets, few variant-to-gene links have emerged. It is unclear whether these sparse results are due to limitations in available data and methods, or to deficiencies in the underlying assumed model. To better distinguish between these possibilities, we identified 220 gene–trait pairs in which protein-coding variants influence a complex trait or its Mendelian cognate. Despite the presence of expression quantitative trait loci near most GWAS associations, by applying a gene-based approach we found limited evidence that the baseline expression of trait-related genes explains GWAS associations, whether using colocalization methods (8% of genes implicated), transcription-wide association (2% of genes implicated), or a combination of regulatory annotations and distance (4% of genes implicated). These results contradict the hypothesis that most complex trait-associated variants coincide with homeostatic expression QTLs, suggesting that better models are needed. The field must confront this deficit and pursue this ‘missing regulation.’

Editor's evaluation

The findings reported here are important because they address the issue of how complex traits arise from their genetic underpinnings. There is an assumption that genetically mediated variation in transcript abundance, usually detected via analysis of expressed quantitative trait loci, is key to this process, but we lack robust evidence in support of that view. This article finds limited evidence that the baseline expression of trait-related genes explains the associations between complex traits and genetic variants (as identified from genome-wide association studies), leading to the view that the field needs to confront a problem of 'missing regulation.'

https://doi.org/10.7554/eLife.74970.sa0

Introduction

Modern complex trait genetics has uncovered surprises at every turn, including the paucity of associations between traits and coding variants of large effect, and the ‘mystery of missing heritability,’ in which no combination of common and rare variants can explain a large fraction of trait heritability (Manolio et al., 2009). Further work has revealed unexpectedly high polygenicity for most human traits and very small effect sizes for individual variants. Enrichment analyses have demonstrated that a large fraction of heritability resides in regions with gene regulatory potential, predominantly tissue-specific accessible chromatin and enhancer elements, suggesting that trait-associated variants influence gene regulation (Maurano et al., 2012; Trynka et al., 2013; Gusev et al., 2014). Furthermore, genes in trait-associated loci are more likely to have genetic variants that affect their expression levels (expression quantitative trait loci, or eQTLs), and the variants with the strongest trait associations are more likely also to be associated with transcript abundance of at least one proximal gene (Nicolae et al., 2010). Combined, these observations have led to the inference that most trait-associated variants are eQTLs, and their effects arise from altering transcript abundance, rather than protein sequence. Equivalent sQTL (splice QTL) analyses of exon usage data have revealed a more modest overlap with trait-associated alleles, suggesting that a fraction of trait-associated variants influence splicing, and hence the relative abundance of different transcript isoforms, rather than overall expression levels. The genetic variant causing expression changes may lie outside the locus and involve a knock-on effect on gene regulation, with the variant altering transcript abundances for genes elsewhere in the genome (a trans-eQTL), but the consensus view is that trans-eQTLs are typically mediated by the variant influencing a gene in the region (a cis-eQTL) (GTEx Consortium, 2020). Thus, a model has emerged in which most trait-associated variants influence proximal gene regulation.

Here we argue that this unembellished model—in which genome-wide association study (GWAS) peaks are mediated by the effects on the homeostatic expression assayed in tissue samples—is the exception rather than the rule. We highlight the challenges of current strategies linking GWAS variants to genes and call for a reevaluation of the basic model in favor of more complex models possibly involving context-specificity with respect to cell types, developmental stages, cell states, or the constancy of expression effects.

Our argument begins with several observations that challenge the unembellished model. One challenge is the difference between spatial distributions of eQTLs, which are dramatically enriched in close proximity to genes, and GWAS peaks, which are usually farther away (Stranger et al., 2007; Farh et al., 2015; Mostafavi et al., 2022). Another is that expression levels mediate a minority of complex trait heritability (Yao et al., 2020). Finally, many studies have designed tools for colocalization analysis: a test of whether GWAS and eQTL associations are due to the same set of variants, not merely distinct variants in linkage disequilibrium. If the model is correct, most trait associations should also be eQTLs, but across studies, only 5–40% of trait associations colocalize with eQTLs (Giambartolomei et al., 2014; Chun et al., 2017; Giambartolomei et al., 2018; Hormozdiari et al., 2016).

Despite the doubts raised, the fact that most GWAS peaks do not colocalize with eQTLs cannot disprove the predominant, unembellished model. In a sense, negative colocalization results are confusing because their hypothesis is too broad. If we predict merely that GWAS peaks will colocalize with some genes’ expression, it is not clear what is meant by a peak’s failure to colocalize with any individual gene’s expression.

Thus, a narrower, more testable hypothesis requires identifying genes we believe a priori are biologically relevant to the GWAS trait. If these trait-linked genes have nearby GWAS peaks and eQTLs, failure to colocalize would be a meaningful negative result. Earlier studies tested all GWAS peaks; when a peak has no colocalization, the model is inconclusive. But trait-linked genes that fail to colocalize reveal that our method for detecting non-coding variation is, with current data, incompatible with our model for understanding it.

With this distinction in mind, we created a set of trait-associated genes capable of supporting or contradicting the model of non-coding GWAS associations acting as eQTLs. For this purpose, the selection of genes becomes extremely important. Because the model attempts to explain the genetic relationship between traits and gene expression, true positives cannot be selected based on measurements of genetic association to traits (GWAS) or expression (eQTL mapping). With this restriction, one source of true positives is to identify genes that are both in loci associated with a complex trait and are also known to harbor coding mutations tied to a related Mendelian trait or the same complex trait. Using a model not based on expression, Mendelian genes are enriched in common-variant heritability for cognate complex traits (Weiner et al., 2022). The genes and their coding variants may be detected in familial studies of cognate Mendelian disorders or by aggregation in a burden test on the same complex phenotypes as GWAS (Backman et al., 2021).

For genes whose coding variants can cause detectable phenotypic change, the strong expectation is that a variant of small effect influences the gene identified by its rare coding variants. As an example, APOE and LDLR are both low-density lipoprotein receptor genes (Schneider et al., 1981; Goldstein and Brown, 1973). Coding variants in APOE and LDLR can lead to the Mendelian disorder familial hypercholesterolemia (Goldstein and Brown, 1973; Cenarro et al., 2016). Even in the absence of a Mendelian coding variant, experiments in animal models have found that the overexpression of these genes reduces cholesterol levels (Shimano et al., 1992a; Shimano et al., 1992b; Kawashiri et al., 2001). GWAS on human subjects have found significant associations near APOE and LDLR, so it seems reasonable to suspect that any non-coding effects in these loci may be mediated by these genes. This general relationship between Mendelian and complex traits is supported by several lines of evidence summarized in Appendix 1.

Results

To test the model that trait-associated variants influence baseline gene expression, we assembled a list of putatively causative genes for seven polygenic common traits with available large-scale GWAS data, each of which also has an extreme form in which coding variants of large effect alter one or more genes with well-characterized biology (Table 1). Our selection included four common diseases: type II diabetes (T2D) (Mahajan et al., 2018), where early-onset familial forms are caused by rare coding mutations (insulin-independent MODY; neonatal diabetes; maternally inherited diabetes and deafness; familial partial lipodystrophy); ulcerative colitis (UC) and Crohn disease (CD) (Liu et al., 2015; Goyette et al., 2015), which have Mendelian pediatric forms characterized by severity of presentation; and breast cancer (BC) (Zhang et al., 2020), where germline coding mutations (e.g., BRCA1) or somatic tissue (e.g., PIK3CA) are sufficient for disease. We also chose three quantitative traits: low- and high-density lipoprotein levels (LDL and HDL); and height. Between known Mendelian genes and those from Backman et al., 2021, our analysis included 220 unique gene-trait pairs (Figure 1).

Table 1
Putatively causative Mendelian genes.

Each gene includes reference(s) to the known biological role of its coding variants, as established in familial studies, in vitro experiments, and/or animal models. Genes from Backman et al., 2021 are not included here, but can be found in Figure 2.

PhenotypeGenes
Low-density lipoproteinAPOB Soria et al., 1989; Pullinger et al., 1995 
APOC2 Hegele et al., 1991 
APOE de Knijff et al., 1994 
LDLR Brown and Goldstein, 1976 
LPL Heizmann et al., 1991; Clee et al., 2001 
PCSK9 Abifadel et al., 2003
High-density lipoproteinABCA1 Brooks-Wilson et al., 1999; Bodzioch et al., 1999; Rust et al., 1999; Ordovas et al., 1986 
APOA1 Ordovas et al., 1986 
CETP Glueck et al., 1975 
LIPC Isaacs et al., 2004; Grarup et al., 2008; Iijima et al., 2008 
LIPG Yamakawa-Kobayashi et al., 2003 
PLTP Jiang et al., 1999 

SCARB1 Tai et al., 2003; McCarthy et al., 2003
Height ANTXR1 Stránecký et al., 2013; Bayram et al., 2014 
ATR O’Driscoll et al., 2003; Ogi et al., 2012 
BLM Ellis et al., 1995; Foucault et al., 1997 
CDC6 Bicknell et al., 2011a 
CDT1 Bicknell et al., 2011a; Guernsey et al., 2011 
CENPJ AlDosari et al., 2010 
COL1A1 Wallis et al., 1990 
COL1A2 Spotila et al., 1992; De Paepe et al., 1997 
COMP Briggs et al., 1995; Mabuchi et al., 2003 
CREBBP Menke et al., 2016; Menke et al., 2018; Angius et al., 2019 
DNA2 Shaheen et al., 2014 
EP300 Woods et al., 2014; Tsai et al., 2011 
EVC Polymeropoulos et al., 1996; Ruiz-Perez et al., 2003 
EVC2 Ruiz-Perez et al., 2003; Galdzicka et al., 2002 
BN1 Faivre et al., 2003; Le Goff et al., 2011; Horn and Robinson, 2011; Takenouchi et al., 2013 
FGFR3 Hyland et al., 2003; Toydemir et al., 2006; Makrythanasis et al., 2014 
FKBP10 Alanay et al., 2010; Kelley et al., 2011; Barnes et al., 2013 
HR Berg et al., 1993; Woods et al., 1996; Goddard et al., 1995; Ayling et al., 1997 
KRAS Aoki et al., 2005; Schubbert et al., 2006; Carta et al., 2006 
NBN Varon et al., 1998; Tanzanella et al., 2003 
NIPBL Tonkin et al., 2004; Krantz et al., 2004 
ORC1 Bicknell et al., 2011a; Guernsey et al., 2011; Bicknell et al., 2011b 
RC4 Guernsey et al., 2011; Bicknell et al., 2011b 
ORC6L Bicknell et al., 2011a; de Munnik et al., 2012 
PCNT Rauch et al., 2008; Griffith et al., 2008; Piane et al., 2009 
PLOD2 van der Slot et al., 2003; HaVinh et al., 2004; Puig-Hervás et al., 2012 
PTPN11 Tartaglia et al., 2001; Maheshwari et al., 2002; Kosaki et al., 2002 
RAD21 Deardorff et al., 2012; Kruszka et al., 2019; Goel and Parasivam, 2020 
RAF1 Pandit et al., 2007; Razzaque et al., 2007 
RECQL4 Lindor et al., 2000; Beghini et al., 2003; Wang et al., 2003 
RIT1 Aoki et al., 2013; Bertola et al., 2014; Gos et al., 2014 
ROR2 Afzal et al., 2000; van Bokhoven, 2000; Tufan et al., 2005 
SLC26A2 Hästbacka et al., 1993; Rossi and Superti-Furga, 2001; Barreda-Bonis et al., 2018 
SMAD4 Le Goff et al., 2012; Caputo et al., 2012; Lindor et al., 2012 
SRCAP Hood et al., 2012; Le Goff et al., 2013 
WRN Yu et al., 1996; Goto et al., 1997; Yu, 1997
Crohn diseaseATG16L1 Hampe et al., 2007 
CARD9 Rivas et al., 2011 
IL10 Fowler, 2005 
IL10RA Gasche et al., 2003; Mao et al., 2012 
IL10RB Glocker et al., 2009; Begue et al., 2011 
IL23R Duerr et al., 2006; Libioulle et al., 2007; Glas et al., 2007 
IRGM McCarroll et al., 2008; Craddock et al., 2010; Prescott et al., 2010 
NOD2 Ogura et al., 2001; Hugot et al., 2001 
PRDM1 Ellinghaus et al., 2013 
PTPN22 Diaz-Gallo et al., 2011 
Ulcerative colitisATG16L1 Fowler et al., 2008 
CARD9 Rivas et al., 2011 
IL23R Glas et al., 2007; Fisher et al., 2008 
IRGM McCarroll et al., 2008 
PRDM1 Ellinghaus et al., 2013 
PTPN22 Diaz-Gallo et al., 2011 
RNF186 Rivas et al., 2016; Beaudoin et al., 2013 
Type II diabetesABCC8 Reis et al., 2000 
BLK Borowiec et al., 2009 
CEL Bengtsson-Ellmark et al., 2004; Raeder et al., 2006 
EIF2AK3 Harding et al., 2001; Brickwood et al., 2003; Durocher et al., 2006 
GATA4 Shaw-Smith, 2014 
GATA6 Yorifuji et al., 2012; De Franco et al., 2013 
GCK Froguel et al., 1993 
GLIS3 Senée et al., 2006 
HNF1A Yamagata et al., 1996b; Vaxillaire et al., 1997 
HNF1B Horikawa et al., 1997; Lindner et al., 1999 
HNF4A Yamagata et al., 1996a; Stoffel and Duncan, 1997 
IER3IP1 Poulton et al., 2011; Abdel-Salam et al., 2012; Shalev et al., 2014 
INS Støy et al., 2007 
KCNJ11 Hani et al., 1998; Gloyn et al., 2004 j
KLF11 Neve et al., 2005 
LMNA Cao and Hegele, 2000 
NEUROD1 Malecki et al., 1999 
NEUROG3 Gradwohl et al., 2000; Rubio-Cabezas et al., 2011; Pinney et al., 2011 
PAX4 Shimajiri et al., 2001; Mauvais-Jarvis et al., 2004; Plengvidhya et al., 2007 
PDX1 Stoffers et al., 1997; Macfarlane et al., 1999; Hani et al., 1999 
PPARG Deeb et al., 1998; Savage et al., 2002 
PTF1A Sellick et al., 2004 
RFX6 Smith, 2010; Sansbury et al., 2015 
SLC19A2 Labay et al., 1999, Oishi et al., 2002; Shaw-Smith et al., 2012 
SLC2A2 Laukkanen et al., 2005; Sansbury et al., 2012 
WFS1 Strom et al., 1998; Hardy et al., 1999; Khanim et al., 2001 
ZFP57 Mackay et al., 2008; Boonen et al., 2013 
Breast cancer (selected using MutPanning; Dietlein et al., 2020)AKT1
ARID1A
ATM
BRCA1
BRCA2
CBFB
CDH1
CDKN1B
CHEK2
CTCF
ERBB2
ESR1
FGFR2
FOXA1
GATA3
GPS2
HS6ST1
KMT2C
KRAS
LRRC37A3
MAP2K4
MAP3K1
NCOR1
NF1
NUP93
PALB2
PIK3CA
PTEN
RB1
RUNX1
SF3B1
STK11
TBX3
TP53
ZFP36L1
Figure 1 with 2 supplements see all
Download asset Open asset
Putatively causative genes identified by each method category.

The leftmost column in each half of the plot displays the entire group of putatively causative genes for our Mendelian set of genes and our (Backman et al., 2021) set of genes, respectively, as well as noting how many are unique to each set or shared between the two sets. The second column in each half indicates how many genes from each set have a nearby GWAS peak or have both a nearby GWAS peak and an expression QTL (eQTL). The remaining columns indicate how many genes were identified through colocalization, transcriptome-wide association studies (TWAS), or chromatin methods, while noting how many of these genes are unique vs. shared between the Mendelian and Backman sets.

In well-powered GWAS, even relatively rare large-effect coding alleles (mutations in BRCA1 that cause breast cancer, for instance) may be detectable as an association to common variants, which could make the effect of a coding variant appear to be regulatory instead. To account for this possibility, we computed association statistics in each GWAS locus conditional on coding variants. We applied a direct conditional test to datasets with available individual-level genotype data (height, LDL, HDL); for those studies without available genotype data, we computed conditional associations from summary statistics using COJO (Yang et al., 2011; Yang et al., 2012; ‘Materials and methods’). With both methods, the resulting GWAS associations should reflect only non-coding variants.

After controlling for coding variation, we examined whether these genes are more likely than chance to be in close proximity to variants associated with the polygenic form of each trait. In agreement with existing literature (Freund et al., 2018), we observe a significant enrichment for all traits in our combined Mendelian and Backman et al., 2021 gene sets (Figure 1—figure supplement 1).

Of our 220 genes, 147 (67%) fell within 1 Mb of a GWAS locus for the cognate complex trait, over three times as many as the 43 predicted by a random null model (95% confidence interval: 31.5–54.5). Our window of 1 Mb represents roughly the upper bound for distances identified between enhancer–promoter pairs, but most pairs are closer (Nasser et al., 2021), so we would expect enrichment to increase as the window around genes decreases; this proves to be the case. At a distance of 100 kb, we find 104 putatively causative genes (47%), though the null model predicts only 11 (95% CI 4.5–17.0), an order-of-magnitude enrichment (Figure 1—figure supplement 1). Given their known causal roles in the severe forms of each phenotype, these results suggest that the 147 genes near GWAS signals are likely to be the targets of trait-associated non-coding variants. For example, we see a significant GWAS association between breast cancer risk and variants in the estrogen receptor (ESR1) locus even after controlling for coding variation; the baseline expression model would thus predict that non-coding risk alleles alter ESR1 expression to drive breast cancer risk.

We next looked for evidence that the trait-associated variants were also altering the expression of our 147 genes in relevant tissues. Controlling for the number of tests we conducted, 134 (91.1%) of these genes had an eQTL in at least one relevant tissue at a false discovery rate (FDR) of Q < 0.05 (‘Materials and methods’). If these variants act through changes in gene expression, phenotypic associations should be driven by some of the same variants as eQTLs in relevant tissue types. We therefore looked for co-localization between our GWAS signals and eQTLs in relevant tissues (Table 2) drawn from the GTEx Project using three well-documented methods: coloc (Giambartolomei et al., 2014), JLIM (Chun et al., 2017), and eCAVIAR (Hormozdiari et al., 2016). We found support for the colocalization of trait and eQTL association for only 7 genes out of 147 (4.8%) for coloc; 10/147 (6.8%) for JLIM; and 8/147 (5.4%) for eCAVIAR. Accounting for overlap, this represents only 18/220 putatively causative genes (8.2%) or 18/147 (12.2%) putatively causative genes near GWAS peaks, even without full multiple-hypothesis testing correction (‘Materials and methods’), which is not obviously better than random chance. We note that prior estimates of the fraction of GWAS associations colocalizing with eQTLs (25–40%; Giambartolomei et al., 2014; Chun et al., 2017; Hormozdiari et al., 2016; Wen et al., 2017) do not directly evaluate the ability to find causative genes. By contrast, our estimate of the number of putatively causative genes that colocalize with eQTLs tests the consistency of our knowledge, models, and data.

Table 2
Tissue-trait pairs.

Tissues were selected for each trait based on a priori knowledge of disease biology.

Mendelian traitGWAS traitTissues examined
Breast cancerBreast cancerBreast mammary tissue
Crohn diseaseCrohn diseaseSmall intestine terminal ileum
Colon sigmoid
Colon transverse
Ulcerative colitisUlcerative colitisSmall intestine terminal ileum
Colon sigmoid
Colon transverse
Dyslipidemia
Hyperlipidemia
Tangier’s disease		High-density lipoproteinLiver
Adipose (subcutaneous)
Whole blood
Dyslipidemia
Hyperlipidemia		Low-density lipoproteinLiver
Adipose (subcutaneous)
Whole blood
Mendelian short statureHeightSkeletal muscle
Monogenic diabetesType II diabetesPancreas
Skeletal muscle
Adipose (subcutaneous)
Small intestine terminal ileum

A potential weakness of our approach is the restriction of our search to predefined tissues. We believe this is necessary in order to avoid the disadvantages of testing each gene–trait pair in each tissue—either a large number of false positives or a severe multiple-testing correction that may lead to false negatives. However, restricting to the set of tissues with a known biological role and available expression data almost certainly leaves out tissues with relevance in certain contexts. Some of the tissues we do use have smaller sample sizes, limiting their power to detect eQTLs with smaller effects.

To address potential shortcomings from the available sample of tissue contexts, we incorporated the Multivariate Adaptive Shrinkage Method (MASH) (Urbut et al., 2019). MASH is a Bayesian method that takes genetic association summary statistics measured across a variety of conditions and, by determining patterns of similarity across conditions, updates the summary statistics of each individual condition. In our case, if an eQTL is difficult to find in a tissue of interest, incorporating information from other tissues may help us detect it. Unlike meta-analysis, this method generates summary statistics that still correspond to a specific tissue.

We ran MASH on every locus used in our earlier analysis using data from all non-brain GTEx tissues (‘Materials and methods’). Rerunning coloc with these modified statistics increased the number of GWAS-eQTL colocalizations across all genes by 26% (from 389 to 489). However, the 100 new colocalizations identified only four additional putatively causative genes (Figure 1—figure supplement 2). These results indicate that tissue-type selection was not the limiting factor in our analysis.

Transcriptome-wide association studies (TWAS) (Gamazon et al., 2015; Gusev et al., 2016; Mancuso et al., 2017; Barbeira et al., 2018) are another class of methods applied to identify causative genes under GWAS peaks using gene expression. TWAS measures genetic correlation between traits and is not designed to avoid correlations caused by LD, which gives it higher power in the case of allelic heterogeneity or poorly typed causative variants (Wainberg et al., 2019). However, while sensitive, TWAS analyses typically yield expansive result sets that include many false positives and are sensitive to the number of tissue types (Wainberg et al., 2019). Results from the FUSION implementation of TWAS (Mancuso et al., 2017) across all tissues identified our putatively causative genes as likely tied to the GWAS peak in 66/220 loci (30%). However, only 4/220 (1.8%) genes were identified by FUSION when we restricted the analysis to relevant tissues.

Given the paucity of expression-mediated GWAS peaks, we asked whether GWAS variants in our study loci reside in likely regulatory sites. Taking the 128 genes in the Mendelian subset of putatively causative genes, we fine-mapped each nearby GWAS association using the SuSiE algorithm (Wang et al., 2020b). For 37 of these genes, we identified at least one high-confidence fine-mapped variant (defined as a variant with posterior inclusion probability, or PIP, greater than 0.7) within 100 kb of the transcription start site. We tested whether these fine-mapped variants fall within (i) regulatory DNA marked by DNase I hypersensitive sites (DHS) Meuleman et al., 2020; (ii) a narrowly mapped active histone modification feature (H3K27ac, H3K4me1, or H3K4me3; Nasser et al., 2021); or (iii) sites marked as an ‘enhancer’ by ChromHMM (Ernst and Kellis, 2012; Ernst and Kellis, 2017; ‘Materials and methods’). As many as 32/37 (86%) genes identified this way have a fine-mapped variants within a candidate regulatory region across all the tissue types examined, or 25/37 (68%) when restricting to trait-relevant tissues (Figure 2); (Supplementary file 1 and Supplementary file 2). Despite this strong evidence that these GWAS associations arise from effects on regulatory DNA, only 5/25 loci (20%) demonstrably correspond to expression effects in our eQTL analysis.

Figure 2
Download asset Open asset
Genes identified as associated with a complex trait by each method.

Columns 'Mend' and 'Backman' indicate whether a gene is from the Mendelian set of putatively causative genes, the Backman et al. set, or both. Subsequent columns indicate whether a gene was identified as a hit using each of our methods: JLIM, coloc, eCaviar, transcriptome-wide association studies (TWAS), and chromatin analysis.

In order to more analogously compare our regulatory feature analysis to our eQTL analysis, we computed ‘activity-by-distance’ (ABD)—a simplification of the 'activity-by-contact' method that provides a measure connecting a distal regulatory region to a target gene (Nasser et al., 2021; Fulco et al., 2019; 'Materials and methods' Figure 3). Taking each locus’s feature with the highest ABD score, we implicate 5/37 (14%) of our Mendelian subset of genes. As such, even when a GWAS association and trait-relevant gene are in the same locus, they are difficult to link, whether using eQTLs or current approaches to integrating chromatin data with target genes.

Figure 3
Download asset Open asset
Chromatin-based causative gene identification.

Following the fine-mapping of GWAS variants, three parallel methods were used. The first identified fine-mapped variants falling within regions annotated as enhancers by ChromHMM. The second identified variants within histone modification features and evaluated their relevance using an activity-by-distance (ABD) score that combined the strength of the feature (i.e., the strength of the acetylation or methylation peak) with its genomic distance to the gene of interest ('Materials and methods'). The third repeated both of these—checking for fine-mapped variants within a region and calculating the ABD score—for DNase I hypersensitivity sites.

Figure 3—source data 1

Gene-level results for linked expression and traits.

https://cdn.elifesciences.org/articles/74970/elife-74970-fig3-data1-v2.txt

Discussion

Overall, our results are strongly consistent with the idea that complex traits are governed by non-coding genetic variants whose effects on phenotype are mediated by their contribution to the regulation of nearby genes. However, these same results are inconsistent with a model that such effects on gene regulation arise from the genetic influences on baseline gene expression that are captured by eQTLs.

The enrichment of putatively causative genes—selected based on existing biological knowledge—near GWAS peaks supports their role in complex traits. Additionally, the enrichment of fine-mapped GWAS variants in likely regulatory regions marked by DHS and other chromatin features lends support to the general model of GWAS associations arising from effects on gene regulation and expression. However, the inability of varied statistical methods to connect GWAS associations and expression argues against the idea that the causative GWAS variants exert their effects via the homeostatic or blended effects on gene expression that are captured by bulk-tissue eQTLs of the sort discovered by broad expression-data collection projects.

Many explanations have been suggested for the limited success of expression-centric methods to explain the mechanisms of GWAS variants. Undirected, broad approaches—including most GWAS-eQTL linking studies—are designed to be largely independent of a priori biological knowledge and hypotheses. This unconstrained focus is ideal for discovery because it delivers the largest number of positive findings, but it is ill-suited for providing explanations for negative results—when you do not know what you were looking for, it is hard to explain why you did not find it. By testing only loci for which there is a strongly suspected contributing gene, we were better able to distinguish which factors might prevent us from identifying GWAS-to-gene links using expression alone.

As a result, we conclude that a number of explanations often considered when evaluating expression-based variant-to-gene methods are not applicable in the context we examined. These include (i) non-expression-mediated mechanisms, (ii) lack of statistical power for GWAS, (iii) the absence of eQTLs for relevant genes, (iv) and underpowered methods for linking expression to GWAS (Table 3).

Table 3
Proposed explanations for negative results under the unembellished model.

Many explanations have been proposed for GWAS associations that are not explained by cis-QTLs. This table details the explanations inconsistent with our results, which are explained in the left column and addressed on the right. Explanations involving more detailed models of gene regulation can be found in Table 4. Two of the explanations addressed here involve violations of the assumptions of our and other expression-based complex trait studies. If coding and non-coding variants affect fundamentally different biological pathways, or if trait associations rarely depend on cis-eQTLs, our methods of mapping regulation to traits would have nothing to uncover. Even in the presence of eQTL-driven trait associations, insufficient power to detect trait associations, to detect eQTL associations, or to link the two would result in predominantly negative results.

Violated assumptions
Genes implicated via coding variants are irrelevant for non-coding associations

  • Our genes are enriched for GWAS associations even after removing the effects of coding variants

  • Loss-of-function variants, which underlie many Mendelian-trait genes, can be thought of as large-effect eQTLs

  • Genes identified from Backman et al., 2021 are not based on cognate phenotypes, but the same complex phenotypes as GWAS

Regulatory mechanisms other than cis-eQTLs

  • Splice QTLs are consistently found to explain less phenotypic variance than eQTLs, and they cannot explain the many GWAS associations that fall within intergenic regions

  • Trans-eQTLs are believed to rely on their effects as cis-eQTLs for other genes; the few exceptions to this model (e.g., CTCF binding sites) are not broadly applicable

Insufficient power
Lack of GWAS power

  • GWAS have been shown to have sufficient power to identify small effects even in rare variants

  • 2/3 of the genes we used have nearby GWAS associations, reflecting a strong enrichment and indicating that GWAS discovery is not a limiting factor

  • Our analysis is conditioned on the presence of GWAS associations

Lack of eQTL mapping power
Lack of power for colocalization and TWAS methods

  1. eQTL = expression QTL; TWAS = transcriptome-wide association studies.

Our results suggest that there remains to be discovered a major component connecting variants to gene expression, which we term 'missing regulation.' We propose that this currently lacking component could be exposed by considering more nuanced experimental and analytical models of gene expression. One likely model involves context-dependent gene expression—expression whose relevant effects are confined to (i) particular refined cell types or anatomical partitions, or (ii) cell states such as responses to perturbations, whether exogenous (e.g., environmental change) or endogenous (e.g., cell differentiation). A complementary model may incorporate heterogeneity of gene expression or the variance of expression across relatively short time scales. These models and others may depend on or be augmented by thresholding or buffering of variant effects, which may produce changes in gene expression with nonlinear effects on phenotype. A summary of proposed models can be found in Table 4.

Table 4
Explaining negative results with more nuanced models of gene regulation.

To reconcile an expression-based model with our observations requires us to both explain the absence of trait-linked eQTLs as well as explaining away the inconsequence of eQTLs for trait-linked genes. The left-hand side lists additions or changes to the unembellished model, while the right-hand side contains explanations of the models and current relevant research.

Extended models of gene regulation
Context dependency:
a context-specific eQTL, invisible in bulk tissues analyzed to date, replaces or supplements the bulk tissue homeostatic eQTL		Cell type Dobbyn et al., 2018; Zhang et al., 2018; Schmiedel et al., 2018; Glastonbury et al., 2019; Rai et al., 2020; Findley et al., 2021; Neavin et al., 2021; Ota et al., 2021; Patel et al., 2021; Bryois et al., 2021; Arvanitis et al., 2022; Oelen et al., 2022; Perez et al., 2022; Schmiedel et al., 2022; Yazar et al., 2022

  • Only a subset of cell types in the tissue contribute to the GWAS phenotype.

  • An eQTL specific to such a cell type is causative for the phenotype.

  • The eQTL either cannot be detected in bulk tissue because of the cell type’s low prevalence

  • The appropriate cellular or anatomical context has not yet been analyzed

Developmental timing Dobbyn et al., 2018; Strober et al., 2019; Cuomo et al., 2020; Bonder et al., 2021; Jerber et al., 2021; Aygün et al., 2022; Elorbany et al., 2022

  • The GWAS phenotype depends on a specific point in cell/tissue development or differentiation

  • eQTLs present at the correct interval contribute to phenotype, but eQTLs observed at other points do not

Cell state or environment Findley et al., 2021; Ota et al., 2021; Oelen et al., 2022; Schmiedel et al., 2022; Huh and Paulsson, 2011; Knowles et al., 2017; Kim-Hellmuth et al., 2017; Balliu et al., 2021; Mu et al., 2021; Ward et al., 2021; Nathan et al., 2022; Baca et al., 2022

  • The causative eQTL has effects that are undetectable in steady-state expression under normal conditions

  • It may activate only in response to a specific environmental condition, such as immune activation or a metabolic shift

Nonlinear or non-homeostatic:
the relationship between eQTL and genotype is indirect		Nonlinearity Fu et al., 2009; Dori-Bachash et al., 2011; Ghazalpour et al., 2011; Pai et al., 2012; Vogel and Marcotte, 2012; Khan et al., 2013; Wu et al., 2013; McManus et al., 2014; Albert and Kruglyak, 2015; Bader et al., 2015; Battle et al., 2015; Cenik et al., 2015; McManus et al., 2015; Pai et al., 2015; Schafer et al., 2015; Chick et al., 2016; Liu et al., 2016; Schaefke et al., 2018; Buccitelli and Selbach, 2020; Wang et al., 2020a; Kusnadi et al., 2022

  • There may be buffering that prevents a change in expression from producing a change in protein levels

  • Expression below a certain level may not influence phenotype, rendering small eQTLs irrelevant

Steady-state expression may be a poor model Pedraza and Paulsson, 2008; Raj and van Oudenaarden, 2008; Shahrezaei and Swain, 2008; Larson et al., 2009; Raj and van Oudenaarden, 2009; Suter et al., 2011; Dar et al., 2012; Viñuelas et al., 2013; Kumar et al., 2015; Nicolas et al., 2017; Qiu et al., 2019; Wang et al., 2020c

  • Phenotype may depend on the kinetics of expression, which could be cyclical or follow some other pattern

  • Expression may be stochastic, such that only a random subset of cells display the relevant expression pattern at any one time

  1. eQTL = expression QTL.

We propose that finding the 'missing regulation' will require not only identifying novel eQTLs explaining GWAS peaks, but also explaining the phenotypic irrelevance of 'red herring eQTLs'—that is, eQTLs for putatively causative genes that fall near GWAS peaks but do not colocalize with them. Above, we use the example genes APOE and LDLR. Both these genes harbor coding variants causing Mendelian hypercholesterolemia, and both have non-coding variants that GWAS have tied to LDL levels. Both have eQTLs in trait-relevant tissues. For APOE, these points cohere into an explanation: the LDL association is an eQTL for the lipid-binding gene. But for LDLR—and for most genes—the association, the mechanism, and the gene cannot be tied together.

Importantly, our results do not diminish the importance or general utility of eQTLs. Rather, they suggest that current models are deficient in two respects: (i) they fail to unify trait-associated non-coding variants with known trait-associated genes, and (ii) they fail to explain the non-effects of identified 'red herring' eQTLs. These deficiencies highlight a need for new approaches to the role of gene regulation in complex traits.

One long-standing goal of GWAS has been to discover genes contributing to complex traits (Manolio et al., 2009; Farh et al., 2015), but low rates of positive findings for expression-based variant-to-gene methods have constrained this possibility (Chun et al., 2017; Baird et al., 2021). Among other challenges, this has limited the benefit of GWAS and expression data for disease-gene mapping and drug discovery (Baird et al., 2021; Umans et al., 2021). Another practical question raised is the value of current large-scale public datasets. Compared to genotypes, expression data are relatively difficult to collect, particularly from specialized cell contexts. If the most relevant models are shown to depend on effects not observable in homeostatic eQTL mapping, the field may need to consider prioritizing other biological contexts and forms of expression data.

Materials and methods

Gene selection

Request a detailed protocol

By manual literature search, we selected 128 genes harboring large-effect-size coding variants for one of the seven phenotypes (Table 1; specifically, we selected 128 gene–trait pairs, representing 121 unique genes). These genes were identified using familial studies, rare disease exome-sequencing analyses, and, for breast cancer, using the MutPanning method (Dietlein et al., 2020) (citations for each gene are included in Table 1). Review papers, as well as the OMIM database (McKusick-nathans institute of genetic medicine, 2021), were generally used as starting points, but an examination of the original literature was needed to confirm genes’ suitability. For example, though SMC3 is known to cause Cornelia de Lange syndrome, which is characterized in part by short stature, SMC3 mutations lead to a milder form of the syndrome, usually without a marked reduction in stature (Deardorff et al., 2007). Several of these phenotypes—height, HDL, cholesterol, breast cancer, and type II diabetes—were also analyzed in Backman et al., 2021, which, through burden testing, identified a total of 110 genes; after accounting for overlaps, this increased our set of putatively causative genes to 220.

The inclusion of genes from Backman et al. ensures that our results are not dependent on an undetected bias in our selection. The set of genes chosen from familial studies offered the advantage that it was selected based on independent methods and data distinct from the large-scale genotyping studies that have characterized the GWAS era. The tradeoff to this was the impossibility of selecting genes through a fully systematic and non-arbitrary process. Because this work was performed in the UKBB, there is some overlap between their data and ours. However, our work did not use exomes, and most of the variants driving their findings are too rare to influence GWAS results. When this is not the case, our decision to condition on coding variants should make the effects used in our work independent from their findings.

Identifying coding variants

Request a detailed protocol

Because GWAS sample sizes are large enough to detect the low-frequency coding variants used to select some of our genes, it is possible that a coding variant would distort the association signal of nearby eQTLs. To minimize this concern, we removed the effects of coding variants on GWAS. Many variants can fall within coding sequences in rare splice variants, so it is important to remove only those variants that appear commonly as coding. These coding variants were selected based on the pext (proportion of expression across transcripts) data (Cummings et al., 2020). Two filters were used. First, we removed genes whose expression in a trait-relevant tissue was below 50% of their maximum expression across tissues. Second, we removed variants that fell within the coding sequence of less than 25% of splice isoforms in that tissue. The remaining variants were used to correct GWAS signal, as explained below. The code for this analysis, and all other quantitative analyses in this paper, can be accessed at https://github.com/NJC12/missing_link_association_function (NJC12, 2023; copy archived at swh:1:rev:46d9072b7cc13f6532203d1494eec4d0f634e092).

GWAS

Request a detailed protocol

For height, LDL cholesterol, and HDL cholesterol, GWAS were performed using genotypic and phenotypic data from the UKBB. In order to avoid confounding, we restricted our sample to the 337K unrelated individuals with genetically determined British ancestry identified by Bycroft et al., 2017. The GWAS were run using Plink 2.0 (Chang et al., 2015), with the covariates age, sex, body mass index (for LDL and HDL only), 10 principal components, and coding variants.

Conditional analysis

Request a detailed protocol

Because UKBB has limited power for breast cancer, Crohn disease, ulcerative colitis, and type II diabetes, we used publicly available summary statistics. The Conditional and Joint Analysis (COJO) (Yang et al., 2011; Yang et al., 2012) program can condition summary statistics on selected variants—in our case, coding variants—by using an LD reference panel. For this reference, we used TOPMed subjects of European ancestry (Taliun et al., 2021). The ancestry of these subjects was confirmed with FastPCA (Galinsky et al., 2016a; Galinsky et al., 2016b), and the relevant data were extracted using bcftools (Danecek et al., 2021). Our conditional GWAS data are available at doi:10.5061/dryad.612jm644q.

Enrichment analysis

Request a detailed protocol

At each distance, the number of Mendelian and non-Mendelian genes within that window around GWAS peaks are counted. p-Values were calculated using Fisher’s exact test (Figure 1, Figure 1—figure supplement 1). Because Mendelian genes may be unusually important beyond our chosen traits, we conduct a set of controls by measuring the enrichment of non-matching Mendelian and complex traits (CD genes and BC GWAS; BC genes and LDL GWAS; LDL genes and UC GWAS; UC genes and height GWAS; height genes and T2D GWAS; T2D genes and HDL GWAS; HDL genes and CD GWAS).

eQTL detection

Request a detailed protocol

eQTL summary statistics were taken from GTEx v7. Some methods detect colocalization with variants that are individually significant, but would not pass a genome-wide threshold (Chun et al., 2017). Because we tested only a subset of genes, we used the Benjamini–Hochberg method (Benjamini and Hochberg, 1995) to calculate the FDR based on the number of tests we conducted multiplied by a correction factor to account for variants that are tested in combination with a gene but are not reported (a factor of 20 closely matched the genome-wide FDR results for GTEx). With this method, 204/220 (93%) of our genes displayed an eQTL, including 134/147 genes with a nearby GWAS peak (91%). Even using the FDR statistics of the GTEx project—which are based on the assumption of testing every gene in every tissue—107/220 (49%) of our genes and 76/147 (52%) of genes near GWAS peaks had an eQTL at Q < 0.05.

Colocalization

Request a detailed protocol

JLIM (Chun et al., 2017) was run using GWAS summary statistics and GTEx v7 genotypes and phenotypes for the tissues listed in Table 2. Coloc (Giambartolomei et al., 2014) was run using GWAS and GTEx v7 summary statistics for the same tissues. eCAVIAR (Hormozdiari et al., 2016) was run using GWAS and GTEx v7 summary statistics for these tissues, and a reference dataset of LD from UKBB (Weissbrod et al., 2020). MASH was run incorporating data from all non-brain tissues, and coloc was rerun using the adjusted values for the same tissues as before.

MASH

Request a detailed protocol

MASH was applied to all GTEx tissues using the mashr R-package (Urbut et al., 2019). We restricted this model to non-brain tissues—which include all of our trait-selected tissues—due to the known tendency of brain and non-brain tissues to cluster separately in expression analysis (Battle et al., 2017; Park et al., 2017; Gamazon et al., 2019).

Fusion (TWAS)

Request a detailed protocol

We used the FUSION implementation of TWAS, which accounts for the possibility of multiple cis-eQTLs linked to the trait-associated variant by jointly calling sets of genes predicted to include the causative gene, to interrogate our 220 loci (Mancuso et al., 2017). FUSION included our putatively causative genes in the set identified as likely relevant to the GWAS peak in 66/220 loci (30%). However, interpretation of this TWAS result is difficult. For many complex traits, TWAS returns a large number of findings (e.g., >150 for LDL cholesterol and >4800 for height). This is in part due to the multiple genes jointly returned at a locus and can also be a result of the large number of tissues and cell types included in the implementation of FUSION. Most hits are found in tissues without any clear relevance to the trait and absent in relevant tissues—LDL, for example, has more TWAS associations between expression and eQTL in prostate adenocarcinoma (24 genes associated), brain prefrontal cortex (23 genes associated), and transformed fibroblasts (21 genes associated) than it does in adipose (16 genes associated), blood (11 genes associated), or liver (5 genes associated). Individual genes were often identified as hits in multiple tissues, but with an inconsistent direction of effect—that is, increased gene expression correlated with an increase in the quantitative trait or disease risk in some tissues, but a decrease in others, which suggests that the gene in question may not be the one whose expression contributes to the complex trait. Because of this possibility and the known biological role of many of our genes, we restricted our results to tissues with established relevance to our traits.

Fine-mapping GWAS hits

Request a detailed protocol

We fine-mapped the GWAS variants located within ±100 kb of our putatively causative genes by applying the SuSiE algorithm (Wang et al., 2020b) on the unconditional summary statistics from the GWAS of breast cancer, Crohn disease, ulcerative colitis, type II diabetes, height, LDL cholesterol, and HDL cholesterol. An LD reference panel from UKBB subjects of European ancestry was used for this analysis. Fine-mapped variants were annotated using snpEff (v4.3t). Only non-coding variants were kept for further analysis.

Functional genomic annotation of fine-mapped hits

Request a detailed protocol

We projected fine-mapped GWAS variants onto active regions of the genome, identified using three alternative approaches: (i) histone modification features, (ii) DHS, and (iii) ChromHMM enhancers.

First, we looked at three histone modification marks, namely, acetylation of histone H3 lysine 27 residues (H3K27ac), mono-methylation of histone H3 lysine 4 residues (H3K4me1), and tri-methylation of lysine 4 residues (H3K4me3) from the Roadmap Epigenomics Project (Wainberg et al., 2019) to identify functional enhancers, which are key contributors of tissue-specific gene regulation. We downloaded imputed narrowPeak sets for H3K27ac, H3K4me1, and H3K4me3 from the Roadmap Epigenomics Project (Wainberg et al., 2019) ftp site (available here) for 14 different tissue types (Supplementary file 1). For each tissue type, we extracted the narrow peaks that are within ±5 Mb of our putatively causative genes. Then following the approach described in Fulco et al. Barbeira et al., 2018, we extended the 150 bp narrow peaks by 175 bp on both sides to arrive at candidate features of 500 bp in length. All features mapping to blacklisted regions (https://sites.google.com/site/anshul kundaje/projects/blacklists) were removed (Kundaje et al., 2015). The remaining features were re-centered around the peak and overlapping features were merged to give the final set of features per histone modification track. The mean activity/strength of a feature (AF) was calculated by taking the geometric mean of the corresponding peak strengths from H3K27ac, H3K4me1, and H3K4me3 marks. We then combined these activity measurements with the linear distances between the features and the transcription start sites of causative genes to compute ABD scores (a simplified version of ABC scores Barbeira et al., 2018) for gene–feature pairs using the following formula:

ABDscoreF,G=AF×1/dF,G∑allfwithin±5MbofGAf×1/df,G

The ABD score can be thought of as a measure of the contribution of a feature, F to the combined regulatory on gene, G. A high ABD score may serve as a proxy for an increased specificity between a chromatin feature and the gene of interest. We projected the fine-mapped variants onto the chromatin features in different tissue types to assess whether there is an enrichment of likely causal GWAS hits in regulatory features near our putatively causative genes. Both proximity (genomic distance) and specificity (ABD scores) were considered to determine the regulatory contribution of the fine-mapped hits.

Next, we looked at the DHS that are considered to be generic markers of the regulatory DNA and can contain genetic variations associated with traits and diseases (Meuleman et al., 2020). We downloaded the index of human DHS along with biosample metadata from https://www.meuleman.org/research/dhsindex/. The index was in hg38 coordinates, which were converted to hg19 coordinates using the online version of the hgLiftOver package (https://genome.ucsc.edu/cgi-bin/hgLiftOver). We created a DHS index for each tissue type relevant to the traits and diseases we analyzed by including all DHS that are present in at least one biosample from a certain tissue type (Supplementary file 2). We then selected DHS that lie within ±100 kb of the TSS of our putatively causative genes. Since DHS are of variable widths, we recentered the summits in a 350 bp window and merged overlapping sites in the same way as we did for other chromatin marks. We calculated the mean activity (AF) by averaging the strengths of all the merged sites. Next, we calculated the activity by distance score for each DHS and gene pair using the same formula described above. Finally, for each fine-mapped variant, we identified all DHS that fall within ±100 kb of the variant.

Finally, we used in silico chromatin state predictions (chromHMM core 15-state model Wainberg et al., 2019) for relevant tissue types (Supplementary file 1) to identify active enhancer regions in the genome. Tissue-specific chromHMM annotations were downloaded from the Roadmap Epigenomics Project Wainberg et al., 2019 ftp site (available here). We considered a fine-mapped variant to fall in an enhancer region if it was housed within a chromHMM segment described as either enhancer, or bivalent enhancer, or genic enhancer. Since chromHMM annotations are not accompanied by activity measurements, the ABD approach could not be applied here.

Appendix 1

Evidence for the relationship between Mendelian and complex traits

More generally, this expectation is supported by several lines of evidence. Comorbidity between Mendelian and complex traits has been used to identify common variants associated with the complex traits (Blair et al., 2013). Early GWAS found associations near genes identified through familial studies of severe disorders (Voight et al., 2010; Chan et al., 2015), and later implicated some of the same genes in complex and Mendelian forms of cardiovascular (Kathiresan and Srivastava, 2012) and neuropsychiatric (Zhu et al., 2014) traits. More recent analyses have found that GWAS associations are enriched in regions near causative genes for cognate Mendelian traits in blood traits (Vuckovic et al., 2020), lipid traits and diabetes (Weiner et al., 2022), as well as a diverse collection of 62 traits (Freund et al., 2018). Another recent method used transcriptomic, proteomic, and epigenomic data to prioritize genes and found that, in a selection of nine phenotypes, the selected genes were enriched for Mendelian genes causing similar traits (Mountjoy et al., 2021). Together, these suggest that genes causing Mendelian traits also influence cognate complex traits, but not through the same coding mechanism.

Genes can also harbor coding variants tied to less severe forms of a trait. These coding variants are more difficult to identify individually as their effect sizes are much smaller. However, the greater number of variants (in aggregate) and freedom from searching for severe segregating traits allows the use of large population datasets. Backman et al., 2021 used burden testing on UK Biobank data to identify genes whose coding variation affects complex traits, finding many genes not identified through familial studies.

Data availability

Numerical data for results is included in Source Data 1. The dataset generated (GWAS summary statistics conditioned on coding variants) can be found at https://doi.org/10.5061/dryad.612jm644q. Code for this project is available at https://github.com/NJC12/missing_link_association_function (copy archived at swh:1:rev:46d9072b7cc13f6532203d1494eec4d0f634e092). UK Biobank data was used and is available by application to the UK Biobank Data Access Committee https://www.ukbiobank.ac.uk/enable-your-research/apply-for-access. TOPMed Whole Genome Sequencing Project - Freeze 5b, Phases 1 and 2 data was used and can be accessed at https://www.nhlbiwgs.org/topmed-whole-genome-sequencing-project-freeze-5b-phases-1-and-2, full list of TOPMed URLs available in Supplementary file 3.

The following data sets were generated
    1. Connally NJ
    (2022) Dryad Digital Repository
    GWAS results conditioned on coding variants.
    https://doi.org/10.5061/dryad.612jm644q

References

    1. Berg MA
    2. Argente J
    3. Chernausek S
    4. Gracia R
    5. Guevara-Aguirre J
    6. Hopp M
    7. Pérez-Jurado L
    8. Rosenbloom A
    9. Toledo SP
    10. Francke U
    (1993)
    Diverse growth hormone receptor gene mutations in Laron syndrome
    American Journal of Human Genetics 52:998–1005.
    1. Craddock N
    2. Hurles ME
    3. Cardin N
    4. Pearson RD
    5. Plagnol V
    6. Robson S
    7. Vukcevic D
    8. Barnes C
    9. Conrad DF
    10. Giannoulatou E
    11. Holmes C
    12. Marchini JL
    13. Stirrups K
    14. Tobin MD
    15. Wain LV
    16. Yau C
    17. Aerts J
    18. Ahmad T
    19. Andrews TD
    20. Arbury H
    21. Attwood A
    22. Auton A
    23. Ball SG
    24. Balmforth AJ
    25. Barrett JC
    26. Barroso I
    27. Barton A
    28. Bennett AJ
    29. Bhaskar S
    30. Blaszczyk K
    31. Bowes J
    32. Brand OJ
    33. Braund PS
    34. Bredin F
    35. Breen G
    36. Brown MJ
    37. Bruce IN
    38. Bull J
    39. Burren OS
    40. Burton J
    41. Byrnes J
    42. Caesar S
    43. Clee CM
    44. Coffey AJ
    45. Connell JMC
    46. Cooper JD
    47. Dominiczak AF
    48. Downes K
    49. Drummond HE
    50. Dudakia D
    51. Dunham A
    52. Ebbs B
    53. Eccles D
    54. Edkins S
    55. Edwards C
    56. Elliot A
    57. Emery P
    58. Evans DM
    59. Evans G
    60. Eyre S
    61. Farmer A
    62. Ferrier IN
    63. Feuk L
    64. Fitzgerald T
    65. Flynn E
    66. Forbes A
    67. Forty L
    68. Franklyn JA
    69. Freathy RM
    70. Gibbs P
    71. Gilbert P
    72. Gokumen O
    73. Gordon-Smith K
    74. Gray E
    75. Green E
    76. Groves CJ
    77. Grozeva D
    78. Gwilliam R
    79. Hall A
    80. Hammond N
    81. Hardy M
    82. Harrison P
    83. Hassanali N
    84. Hebaishi H
    85. Hines S
    86. Hinks A
    87. Hitman GA
    88. Hocking L
    89. Howard E
    90. Howard P
    91. Howson JMM
    92. Hughes D
    93. Hunt S
    94. Isaacs JD
    95. Jain M
    96. Jewell DP
    97. Johnson T
    98. Jolley JD
    99. Jones IR
    100. Jones LA
    101. Kirov G
    102. Langford CF
    103. Lango-Allen H
    104. Lathrop GM
    105. Lee J
    106. Lee KL
    107. Lees C
    108. Lewis K
    109. Lindgren CM
    110. Maisuria-Armer M
    111. Maller J
    112. Mansfield J
    113. Martin P
    114. Massey DCO
    115. McArdle WL
    116. McGuffin P
    117. McLay KE
    118. Mentzer A
    119. Mimmack ML
    120. Morgan AE
    121. Morris AP
    122. Mowat C
    123. Myers S
    124. Newman W
    125. Nimmo ER
    126. O’Donovan MC
    127. Onipinla A
    128. Onyiah I
    129. Ovington NR
    130. Owen MJ
    131. Palin K
    132. Parnell K
    133. Pernet D
    134. Perry JRB
    135. Phillips A
    136. Pinto D
    137. Prescott NJ
    138. Prokopenko I
    139. Quail MA
    140. Rafelt S
    141. Rayner NW
    142. Redon R
    143. Reid DM
    144. Ring SM
    145. Robertson N
    146. Russell E
    147. St Clair D
    148. Sambrook JG
    149. Sanderson JD
    150. Schuilenburg H
    151. Scott CE
    152. Scott R
    153. Seal S
    154. Shaw-Hawkins S
    155. Shields BM
    156. Simmonds MJ
    157. Smyth DJ
    158. Somaskantharajah E
    159. Spanova K
    160. Steer S
    161. Stephens J
    162. Stevens HE
    163. Stone MA
    164. Su Z
    165. Symmons DPM
    166. Thompson JR
    167. Thomson W
    168. Travers ME
    169. Turnbull C
    170. Valsesia A
    171. Walker M
    172. Walker NM
    173. Wallace C
    174. Warren-Perry M
    175. Watkins NA
    176. Webster J
    177. Weedon MN
    178. Wilson AG
    179. Woodburn M
    180. Wordsworth BP
    181. Young AH
    182. Zeggini E
    183. Carter NP
    184. Frayling TM
    185. Lee C
    186. McVean G
    187. Munroe PB
    188. Palotie A
    189. Sawcer SJ
    190. Scherer SW
    191. Strachan DP
    192. Tyler-Smith C
    193. Brown MA
    194. Burton PR
    195. Caulfield MJ
    196. Compston A
    197. Farrall M
    198. Gough SCL
    199. Hall AS
    200. Hattersley AT
    201. Hill AVS
    202. Mathew CG
    203. Pembrey M
    204. Satsangi J
    205. Stratton MR
    206. Worthington J
    207. Deloukas P
    208. Duncanson A
    209. Kwiatkowski DP
    210. McCarthy MI
    211. Ouwehand W
    212. Parkes M
    213. Rahman N
    214. Todd JA
    215. Samani NJ
    216. Donnelly P
    217. Wellcome Trust Case Control Consortium
    (2010) Genome-wide association study of cnvs in 16,000 cases of eight common diseases and 3,000 shared controls
    Nature 464:713–720.
    https://doi.org/10.1038/nature08979
    1. Glueck CJ
    2. Fallat RW
    3. Millett F
    4. Steiner PM
    (1975)
    Familial hyperalphalipoproteinemia
    Archives of Internal Medicine 135:1025–1028.
    1. Hegele RA
    2. Connelly PW
    3. Maguire GF
    4. Huff MW
    5. Leiter L
    6. Wolfe BM
    7. Evans AJ
    8. Little JA
    (1991)
    An apolipoprotein CII mutation, ciilys19 -- -- ’hr' identified in patients with hyperlipidemia
    Disease Markers 9:73–80.
  176. Book
    1. McKusick-nathans institute of genetic medicine
    (2021)
    Online Mendelian Inheritance in Man, OMIM
    Baltimore, MD: johns hopkins university.
    1. Spotila LD
    2. Sereda L
    3. Prockop DJ
    (1992)
    Partial isodisomy for maternal chromosome 7 and short stature in an individual with a mutation at the COL1A2 locus
    American Journal of Human Genetics 51:1396–1405.
    1. Taliun D
    2. Harris DN
    3. Kessler MD
    4. Carlson J
    5. Szpiech ZA
    6. Torres R
    7. Taliun SAG
    8. Corvelo A
    9. Gogarten SM
    10. Kang HM
    11. Pitsillides AN
    12. LeFaive J
    13. Lee S-B
    14. Tian X
    15. Browning BL
    16. Das S
    17. Emde A-K
    18. Clarke WE
    19. Loesch DP
    20. Shetty AC
    21. Blackwell TW
    22. Smith AV
    23. Wong Q
    24. Liu X
    25. Conomos MP
    26. Bobo DM
    27. Aguet F
    28. Albert C
    29. Alonso A
    30. Ardlie KG
    31. Arking DE
    32. Aslibekyan S
    33. Auer PL
    34. Barnard J
    35. Barr RG
    36. Barwick L
    37. Becker LC
    38. Beer RL
    39. Benjamin EJ
    40. Bielak LF
    41. Blangero J
    42. Boehnke M
    43. Bowden DW
    44. Brody JA
    45. Burchard EG
    46. Cade BE
    47. Casella JF
    48. Chalazan B
    49. Chasman DI
    50. Chen Y-DI
    51. Cho MH
    52. Choi SH
    53. Chung MK
    54. Clish CB
    55. Correa A
    56. Curran JE
    57. Custer B
    58. Darbar D
    59. Daya M
    60. de Andrade M
    61. DeMeo DL
    62. Dutcher SK
    63. Ellinor PT
    64. Emery LS
    65. Eng C
    66. Fatkin D
    67. Fingerlin T
    68. Forer L
    69. Fornage M
    70. Franceschini N
    71. Fuchsberger C
    72. Fullerton SM
    73. Germer S
    74. Gladwin MT
    75. Gottlieb DJ
    76. Guo X
    77. Hall ME
    78. He J
    79. Heard-Costa NL
    80. Heckbert SR
    81. Irvin MR
    82. Johnsen JM
    83. Johnson AD
    84. Kaplan R
    85. Kardia SLR
    86. Kelly T
    87. Kelly S
    88. Kenny EE
    89. Kiel DP
    90. Klemmer R
    91. Konkle BA
    92. Kooperberg C
    93. Köttgen A
    94. Lange LA
    95. Lasky-Su J
    96. Levy D
    97. Lin X
    98. Lin K-H
    99. Liu C
    100. Loos RJF
    101. Garman L
    102. Gerszten R
    103. Lubitz SA
    104. Lunetta KL
    105. Mak ACY
    106. Manichaikul A
    107. Manning AK
    108. Mathias RA
    109. McManus DD
    110. McGarvey ST
    111. Meigs JB
    112. Meyers DA
    113. Mikulla JL
    114. Minear MA
    115. Mitchell BD
    116. Mohanty S
    117. Montasser ME
    118. Montgomery C
    119. Morrison AC
    120. Murabito JM
    121. Natale A
    122. Natarajan P
    123. Nelson SC
    124. North KE
    125. O’Connell JR
    126. Palmer ND
    127. Pankratz N
    128. Peloso GM
    129. Peyser PA
    130. Pleiness J
    131. Post WS
    132. Psaty BM
    133. Rao DC
    134. Redline S
    135. Reiner AP
    136. Roden D
    137. Rotter JI
    138. Ruczinski I
    139. Sarnowski C
    140. Schoenherr S
    141. Schwartz DA
    142. Seo J-S
    143. Seshadri S
    144. Sheehan VA
    145. Sheu WH
    146. Shoemaker MB
    147. Smith NL
    148. Smith JA
    149. Sotoodehnia N
    150. Stilp AM
    151. Tang W
    152. Taylor KD
    153. Telen M
    154. Thornton TA
    155. Tracy RP
    156. Van Den Berg DJ
    157. Vasan RS
    158. Viaud-Martinez KA
    159. Vrieze S
    160. Weeks DE
    161. Weir BS
    162. Weiss ST
    163. Weng L-C
    164. Willer CJ
    165. Zhang Y
    166. Zhao X
    167. Arnett DK
    168. Ashley-Koch AE
    169. Barnes KC
    170. Boerwinkle E
    171. Gabriel S
    172. Gibbs R
    173. Rice KM
    174. Rich SS
    175. Silverman EK
    176. Qasba P
    177. Gan W
    178. Papanicolaou GJ
    179. Nickerson DA
    180. Browning SR
    181. Zody MC
    182. Zöllner S
    183. Wilson JG
    184. Cupples LA
    185. Laurie CC
    186. Jaquish CE
    187. Hernandez RD
    188. O’Connor TD
    189. Abecasis GR
    190. NHLBI Trans-Omics for Precision Medicine Consortium
    (2021) Sequencing of 53,831 diverse genomes from the NHLBI topmed program
    Nature 590:290–299.
    https://doi.org/10.1038/s41586-021-03205-y
    1. Voight BF
    2. Scott LJ
    3. Steinthorsdottir V
    4. Morris AP
    5. Dina C
    6. Welch RP
    7. Zeggini E
    8. Huth C
    9. Aulchenko YS
    10. Thorleifsson G
    11. McCulloch LJ
    12. Ferreira T
    13. Grallert H
    14. Amin N
    15. Wu G
    16. Willer CJ
    17. Raychaudhuri S
    18. McCarroll SA
    19. Langenberg C
    20. Hofmann OM
    21. Dupuis J
    22. Qi L
    23. Segrè AV
    24. van Hoek M
    25. Navarro P
    26. Ardlie K
    27. Balkau B
    28. Benediktsson R
    29. Bennett AJ
    30. Blagieva R
    31. Boerwinkle E
    32. Bonnycastle LL
    33. Bengtsson Boström K
    34. Bravenboer B
    35. Bumpstead S
    36. Burtt NP
    37. Charpentier G
    38. Chines PS
    39. Cornelis M
    40. Couper DJ
    41. Crawford G
    42. Doney ASF
    43. Elliott KS
    44. Elliott AL
    45. Erdos MR
    46. Fox CS
    47. Franklin CS
    48. Ganser M
    49. Gieger C
    50. Grarup N
    51. Green T
    52. Griffin S
    53. Groves CJ
    54. Guiducci C
    55. Hadjadj S
    56. Hassanali N
    57. Herder C
    58. Isomaa B
    59. Jackson AU
    60. Johnson PRV
    61. Jørgensen T
    62. Kao WHL
    63. Klopp N
    64. Kong A
    65. Kraft P
    66. Kuusisto J
    67. Lauritzen T
    68. Li M
    69. Lieverse A
    70. Lindgren CM
    71. Lyssenko V
    72. Marre M
    73. Meitinger T
    74. Midthjell K
    75. Morken MA
    76. Narisu N
    77. Nilsson P
    78. Owen KR
    79. Payne F
    80. Perry JRB
    81. Petersen A-K
    82. Platou C
    83. Proença C
    84. Prokopenko I
    85. Rathmann W
    86. Rayner NW
    87. Robertson NR
    88. Rocheleau G
    89. Roden M
    90. Sampson MJ
    91. Saxena R
    92. Shields BM
    93. Shrader P
    94. Sigurdsson G
    95. Sparsø T
    96. Strassburger K
    97. Stringham HM
    98. Sun Q
    99. Swift AJ
    100. Thorand B
    101. Tichet J
    102. Tuomi T
    103. van Dam RM
    104. van Haeften TW
    105. van Herpt T
    106. van Vliet-Ostaptchouk JV
    107. Walters GB
    108. Weedon MN
    109. Wijmenga C
    110. Witteman J
    111. Bergman RN
    112. Cauchi S
    113. Collins FS
    114. Gloyn AL
    115. Gyllensten U
    116. Hansen T
    117. Hide WA
    118. Hitman GA
    119. Hofman A
    120. Hunter DJ
    121. Hveem K
    122. Laakso M
    123. Mohlke KL
    124. Morris AD
    125. Palmer CNA
    126. Pramstaller PP
    127. Rudan I
    128. Sijbrands E
    129. Stein LD
    130. Tuomilehto J
    131. Uitterlinden A
    132. Walker M
    133. Wareham NJ
    134. Watanabe RM
    135. Abecasis GR
    136. Boehm BO
    137. Campbell H
    138. Daly MJ
    139. Hattersley AT
    140. Hu FB
    141. Meigs JB
    142. Pankow JS
    143. Pedersen O
    144. Wichmann H-E
    145. Barroso I
    146. Florez JC
    147. Frayling TM
    148. Groop L
    149. Sladek R
    150. Thorsteinsdottir U
    151. Wilson JF
    152. Illig T
    153. Froguel P
    154. van Duijn CM
    155. Stefansson K
    156. Altshuler D
    157. Boehnke M
    158. McCarthy MI
    159. MAGIC investigators
    160. GIANT Consortium
    (2010) Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis
    Nature Genetics 42:579–589.
    https://doi.org/10.1038/ng.609
    1. Wallis GA
    2. Starman BJ
    3. Zinn AB
    4. Byers PH
    (1990)
    Variable expression of osteogenesis imperfecta in a nuclear family is explained by somatic mosaicism for a lethal point mutation in the alpha 1 (I) gene (COL1A1) of type I collagen in a parent
    American Journal of Human Genetics 46:1034–1040.
    1. Yu CE
    (1997)
    Mutations in the consensus helicase domains of the werner syndrome gene
    American Journal of Human Genetics 12:1.
    1. Zhang H
    2. Ahearn TU
    3. Lecarpentier J
    4. Barnes D
    5. Beesley J
    6. Qi G
    7. Jiang X
    8. O’Mara TA
    9. Zhao N
    10. Bolla MK
    11. Dunning AM
    12. Dennis J
    13. Wang Q
    14. Ful ZA
    15. Aittomäki K
    16. Andrulis IL
    17. Anton-Culver H
    18. Arndt V
    19. Aronson KJ
    20. Arun BK
    21. Auer PL
    22. Azzollini J
    23. Barrowdale D
    24. Becher H
    25. Beckmann MW
    26. Behrens S
    27. Benitez J
    28. Bermisheva M
    29. Bialkowska K
    30. Blanco A
    31. Blomqvist C
    32. Bogdanova NV
    33. Bojesen SE
    34. Bonanni B
    35. Bondavalli D
    36. Borg A
    37. Brauch H
    38. Brenner H
    39. Briceno I
    40. Broeks A
    41. Brucker SY
    42. Brüning T
    43. Burwinkel B
    44. Buys SS
    45. Byers H
    46. Caldés T
    47. Caligo MA
    48. Calvello M
    49. Campa D
    50. Castelao JE
    51. Chang-Claude J
    52. Chanock SJ
    53. Christiaens M
    54. Christiansen H
    55. Chung WK
    56. Claes KBM
    57. Clarke CL
    58. Cornelissen S
    59. Couch FJ
    60. Cox A
    61. Cross SS
    62. Czene K
    63. Daly MB
    64. Devilee P
    65. Diez O
    66. Domchek SM
    67. Dörk T
    68. Dwek M
    69. Eccles DM
    70. Ekici AB
    71. Evans DG
    72. Fasching PA
    73. Figueroa J
    74. Foretova L
    75. Fostira F
    76. Friedman E
    77. Frost D
    78. Gago-Dominguez M
    79. Gapstur SM
    80. Garber J
    81. García-Sáenz JA
    82. Gaudet MM
    83. Gayther SA
    84. Giles GG
    85. Godwin AK
    86. Goldberg MS
    87. Goldgar DE
    88. González-Neira A
    89. Greene MH
    90. Gronwald J
    91. Guénel P
    92. Häberle L
    93. Hahnen E
    94. Haiman CA
    95. Hake CR
    96. Hall P
    97. Hamann U
    98. Harkness EF
    99. Heemskerk-Gerritsen BAM
    100. Hillemanns P
    101. Hogervorst FBL
    102. Holleczek B
    103. Hollestelle A
    104. Hooning MJ
    105. Hoover RN
    106. Hopper JL
    107. Howell A
    108. Huebner H
    109. Hulick PJ
    110. Imyanitov EN
    111. kConFab Investigators
    112. ABCTB Investigators
    113. Isaacs C
    114. Izatt L
    115. Jager A
    116. Jakimovska M
    117. Jakubowska A
    118. James P
    119. Janavicius R
    120. Janni W
    121. John EM
    122. Jones ME
    123. Jung A
    124. Kaaks R
    125. Kapoor PM
    126. Karlan BY
    127. Keeman R
    128. Khan S
    129. Khusnutdinova E
    130. Kitahara CM
    131. Ko Y-D
    132. Konstantopoulou I
    133. Koppert LB
    134. Koutros S
    135. Kristensen VN
    136. Laenkholm A-V
    137. Lambrechts D
    138. Larsson SC
    139. Laurent-Puig P
    140. Lazaro C
    141. Lazarova E
    142. Lejbkowicz F
    143. Leslie G
    144. Lesueur F
    145. Lindblom A
    146. Lissowska J
    147. Lo W-Y
    148. Loud JT
    149. Lubinski J
    150. Lukomska A
    151. MacInnis RJ
    152. Mannermaa A
    153. Manoochehri M
    154. Manoukian S
    155. Margolin S
    156. Martinez ME
    157. Matricardi L
    158. McGuffog L
    159. McLean C
    160. Mebirouk N
    161. Meindl A
    162. Menon U
    163. Miller A
    164. Mingazheva E
    165. Montagna M
    166. Mulligan AM
    167. Mulot C
    168. Muranen TA
    169. Nathanson KL
    170. Neuhausen SL
    171. Nevanlinna H
    172. Neven P
    173. Newman WG
    174. Nielsen FC
    175. Nikitina-Zake L
    176. Nodora J
    177. Offit K
    178. Olah E
    179. Olopade OI
    180. Olsson H
    181. Orr N
    182. Papi L
    183. Papp J
    184. Park-Simon T-W
    185. Parsons MT
    186. Peissel B
    187. Peixoto A
    188. Peshkin B
    189. Peterlongo P
    190. Peto J
    191. Phillips K-A
    192. Piedmonte M
    193. Plaseska-Karanfilska D
    194. Prajzendanc K
    195. Prentice R
    196. Prokofyeva D
    197. Rack B
    198. Radice P
    199. Ramus SJ
    200. Rantala J
    201. Rashid MU
    202. Rennert G
    203. Rennert HS
    204. Risch HA
    205. Romero A
    206. Rookus MA
    207. Rübner M
    208. Rüdiger T
    209. Saloustros E
    210. Sampson S
    211. Sandler DP
    212. Sawyer EJ
    213. Scheuner MT
    214. Schmutzler RK
    215. Schneeweiss A
    216. Schoemaker MJ
    217. Schöttker B
    218. Schürmann P
    219. Senter L
    220. Sharma P
    221. Sherman ME
    222. Shu X-O
    223. Singer CF
    224. Smichkoska S
    225. Soucy P
    226. Southey MC
    227. Spinelli JJ
    228. Stone J
    229. Stoppa-Lyonnet D
    230. EMBRACE Study
    231. GEMO Study Collaborators
    232. Swerdlow AJ
    233. Szabo CI
    234. Tamimi RM
    235. Tapper WJ
    236. Taylor JA
    237. Teixeira MR
    238. Terry M
    239. Thomassen M
    240. Thull DL
    241. Tischkowitz M
    242. Toland AE
    243. Tollenaar RAEM
    244. Tomlinson I
    245. Torres D
    246. Troester MA
    247. Truong T
    248. Tung N
    249. Untch M
    250. Vachon CM
    251. van den Ouweland AMW
    252. van der Kolk LE
    253. van Veen EM
    254. vanRensburg EJ
    255. Vega A
    256. Wappenschmidt B
    257. Weinberg CR
    258. Weitzel JN
    259. Wildiers H
    260. Winqvist R
    261. Wolk A
    262. Yang XR
    263. Yannoukakos D
    264. Zheng W
    265. Zorn KK
    266. Milne RL
    267. Kraft P
    268. Simard J
    269. Pharoah PDP
    270. Michailidou K
    271. Antoniou AC
    272. Schmidt MK
    273. Chenevix-Trench G
    274. Easton DF
    275. Chatterjee N
    276. García-Closas M
    (2020) Genome-Wide association study identifies 32 novel breast cancer susceptibility loci from overall and subtype-specific analyses
    Nature Genetics 52:572–581.
    https://doi.org/10.1038/s41588-020-0609-2

Article and author information

Author details

  1. Noah J Connally

    1. Department of Biomedical Informatics, Harvard Medical School, Boston, United States
    2. Brigham and Women’s Hospital, Division of Genetics, Harvard Medical School, Boston, United States
    3. Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, United States
    Contribution
    Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing
    For correspondence
    noahconnally@g.harvard.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3818-6739

  2. Sumaiya Nazeen

    1. Department of Biomedical Informatics, Harvard Medical School, Boston, United States
    2. Brigham and Women’s Hospital, Division of Genetics, Harvard Medical School, Boston, United States
    3. Brigham and Women’s Hospital, Department of Neurology, Harvard Medical School, Boston, United States
    Contribution
    Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared

  3. Daniel Lee

    1. Department of Biomedical Informatics, Harvard Medical School, Boston, United States
    2. Brigham and Women’s Hospital, Division of Genetics, Harvard Medical School, Boston, United States
    3. Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, United States
    Contribution
    Data curation, Formal analysis, Investigation, Writing – review and editing
    Competing interests
    No competing interests declared

  4. Huwenbo Shi

    1. Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, United States
    2. Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, United States
    Contribution
    Resources, Data curation, Software, Formal analysis, Investigation, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared

  5. John Stamatoyannopoulos

    Altius Institute, Seattle, United States
    Contribution
    Resources, Funding acquisition, Writing – review and editing
    Competing interests
    No competing interests declared

  6. Sung Chun

    Division of Pulmonary Medicine, Boston Children’s Hospital, Boston, United States
    Contribution
    Conceptualization, Supervision, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared

  7. Chris Cotsapas

    1. Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, United States
    2. Department of Neurology, Yale Medical School, New Haven, United States
    3. Department of Genetics, Yale Medical School, New Haven, United States
    Contribution
    Conceptualization, Supervision, Methodology, Writing – review and editing
    For correspondence
    cotsapas@broadinstitute.org
    Competing interests
    No competing interests declared

  8. Christopher A Cassa

    1. Brigham and Women’s Hospital, Division of Genetics, Harvard Medical School, Boston, United States
    2. Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, United States
    Contribution
    Conceptualization, Resources, Supervision, Methodology, Project administration, Writing – review and editing
    For correspondence
    cassa@mit.edu
    Competing interests
    No competing interests declared

  9. Shamil R Sunyaev

    1. Department of Biomedical Informatics, Harvard Medical School, Boston, United States
    2. Brigham and Women’s Hospital, Division of Genetics, Harvard Medical School, Boston, United States
    3. Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, United States
    Contribution
    Conceptualization, Resources, Supervision, Funding acquisition, Methodology, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    ssunyaev@hms.harvard.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5715-5677

Funding

National Institutes of Health (R35GM127131)

  • Shamil R Sunyaev

National Institutes of Health (R01HG010372)

  • Shamil R Sunyaev

National Institutes of Health (R01MH101244)

  • Shamil R Sunyaev

National Institutes of Health (U01HG012009)

  • Chris Cotsapas

National Institutes of Health (T32GM74897)

  • Shamil R Sunyaev

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Alkes Price, Alex Bloemendal, Benjamin Neale, Bogdan Pasanuic, Sasha (Alexander) Gusev, and Matt Warman for their helpful discussions. This research was supported by NIH grants R35GM127131, R01HG010372, R01MH101244, and U01HG012009. NJC was supported by NIH training grant T32GM74897. UK Biobank was accessed under projects 14048 and 10438. TOPMed data were used under dbGaP project 28674.

Version history

  1. Preprint posted: June 11, 2021 (view preprint)
  2. Received: October 25, 2021
  3. Accepted: December 2, 2022
  4. Accepted Manuscript published: December 14, 2022 (version 1)
  5. Version of Record published: January 16, 2023 (version 2)

Copyright

© 2022, Connally et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 5,829
    Page views
  • 941
    Downloads
  • 32
    Citations

Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Noah J Connally
  2. Sumaiya Nazeen
  3. Daniel Lee
  4. Huwenbo Shi
  5. John Stamatoyannopoulos
  6. Sung Chun
  7. Chris Cotsapas
  8. Christopher A Cassa
  9. Shamil R Sunyaev
(2022)
The missing link between genetic association and regulatory function
eLife 11:e74970.
https://doi.org/10.7554/eLife.74970

Share this article

https://doi.org/10.7554/eLife.74970

Categories and tags

Research organism

Further reading

    1. Genetics and Genomics
    2. Immunology and Inflammation

    Transposable elements regulate thymus development and function

    Jean-David Larouche, Céline M Laumont ... Claude Perreault
    Research Article

    Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/β. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

    1. Genetics and Genomics

    Deciphering the complex relationship between type 2 diabetes and fracture risk with both genetic and observational evidence

    Pianpian Zhao, Zhifeng Sheng ... Hou-Feng Zheng
    Research Article

    The ‘diabetic bone paradox’ suggested that type 2 diabetes (T2D) patients would have higher areal bone mineral density (BMD) but higher fracture risk than individuals without T2D. In this study, we found that the genetically predicted T2D was associated with higher BMD and lower risk of fracture in both weighted genetic risk score (wGRS) and two-sample Mendelian randomization (MR) analyses. We also identified ten genomic loci shared between T2D and fracture, with the top signal at SNP rs4580892 in the intron of gene RSPO3. And the higher expression in adipose subcutaneous and higher protein level in plasma of RSPO3 were associated with increased risk of T2D, but decreased risk of fracture. In the prospective study, T2D was observed to be associated with higher risk of fracture, but BMI mediated 30.2% of the protective effect. However, when stratified by the T2D-related risk factors for fracture, we observed that the effect of T2D on the risk of fracture decreased when the number of T2D-related risk factors decreased, and the association became non-significant if the T2D patients carried none of the risk factors. In conclusion, the genetically determined T2D might not be associated with higher risk of fracture. And the shared genetic architecture between T2D and fracture suggested a top signal around RSPO3 gene. The observed effect size of T2D on fracture risk decreased if the T2D-related risk factors could be eliminated. Therefore, it is important to manage the complications of T2D to prevent the risk of fracture.

    1. Genetics and Genomics

    The BTB-ZF gene Bm-mamo regulates pigmentation in silkworm caterpillars

    Songyuan Wu, Xiaoling Tong ... Fangyin Dai
    Research Article

    The color pattern of insects is one of the most diverse adaptive evolutionary phenotypes. However, the molecular regulation of this color pattern is not fully understood. In this study, we found that the transcription factor Bm-mamo is responsible for black dilute (bd) allele mutations in the silkworm. Bm-mamo belongs to the BTB zinc finger family and is orthologous to mamo in Drosophila melanogaster. This gene has a conserved function in gamete production in Drosophila and silkworms and has evolved a pleiotropic function in the regulation of color patterns in caterpillars. Using RNAi and clustered regularly interspaced short palindromic repeats (CRISPR) technology, we showed that Bm-mamo is a repressor of dark melanin patterns in the larval epidermis. Using in vitro binding assays and gene expression profiling in wild-type and mutant larvae, we also showed that Bm-mamo likely regulates the expression of related pigment synthesis and cuticular protein genes in a coordinated manner to mediate its role in color pattern formation. This mechanism is consistent with the dual role of this transcription factor in regulating both the structure and shape of the cuticle and the pigments that are embedded within it. This study provides new insight into the regulation of color patterns as well as into the construction of more complex epidermal features in some insects.