Adult height is one of the earliest putative examples of polygenic adaptation in humans. However, this conclusion was recently challenged because residual uncorrected stratification from large-scale consortium studies was considered responsible for the previously noted genetic difference. It thus remains an open question whether height loci exhibit signals of polygenic adaptation in any human population. We re-examined this question, focusing on one of the shortest European populations, the Sardinians, in addition to mainland European populations. We utilized height-associated loci from the Biobank Japan (BBJ) dataset to further alleviate concerns of biased ascertainment of GWAS loci and showed that the Sardinians remain significantly shorter than expected under neutrality (∼0.22 standard deviation shorter than Utah residents with ancestry from northern and western Europe [CEU] on the basis of polygenic height scores, p = 3.89 × 10-4). We also found the trajectory of polygenic height scores between the Sardinian and the British populations diverged over at least the last 10,000 years (p = 0.0082), consistent with a signature of polygenic adaptation driven primarily by the Sardinian population. Although the polygenic score-based analysis showed a much subtler signature in mainland European populations, we found a clear and robust adaptive signature in the UK population by using a haplotype-based statistic, the trait singleton density score (tSDS), driven by the height-increasing alleles (p = 9.1 × 10-4). In summary, by ascertaining height loci in a distant East Asian population, we further supported the evidence of polygenic adaptation at height-associated loci among the Sardinians. In mainland Europeans, the adaptive signature was detected in haplotype-based analysis but not in polygenic score-based analysis.

Genetic variants of interleukin-3 (IL-3), a well-studied cytokine, may have a role in the pathophysiology of rheumatoid arthritis (RA); but reports on this association sometimes conflict. A case-control study was designed to investigate association between RA and a single-nucleotide polymorphism (SNP) in the IL-3 promoter region. Comparison of cases of RA versus control individuals yielded a chi(2) value of 14.28 (P=.0002), with a genotype odds ratio of 2.24 (95% confidence interval [95%CI] 1.44-3.49). When female cases with earlier onset were compared with female control individuals, the SNP revealed an even more significant correlation, with chi2=21.75 (P=.000004) and a genotype odds ratio of 7.27 (95%CI 2.80-18.89). The stronger association that we observed in this clinically distinct subgroup (females with early onset), within a region where linkage disequilibrium was not significantly extended, suggested that the genuine RA locus should locate either within or close to the IL-3 gene. Combined genotype data on SNPs on eight other candidate genes were combined with our IL-3 results, to estimate relationships between pairs of loci and RA, by maximum-likelihood analysis. The utility of combining the genotype data in this way to identify possible contributions of various genes to this disease is discussed.

We have investigated the evidence for positive selection in samples of African, European, and East Asian ancestry at 65 loci associated with susceptibility to type 2 diabetes (T2D) previously identified through genome-wide association studies. Selection early in human evolutionary history is predicted to lead to ancestral risk alleles shared between populations, whereas late selection would result in population-specific signals at derived risk alleles. By using a wide variety of tests based on the site frequency spectrum, haplotype structure, and population differentiation, we found no global signal of enrichment for positive selection when we considered all T2D risk loci collectively. However, in a locus-by-locus analysis, we found nominal evidence for positive selection at 14 of the loci. Selection favored the protective and risk alleles in similar proportions, rather than the risk alleles specifically as predicted by the thrifty gene hypothesis, and may not be related to influence on diabetes. Overall, we conclude that past positive selection has not been a powerful influence driving the prevalence of T2D risk alleles.

Circulating cytokines and growth factors are regulators of inflammation and have been implicated in autoimmune and metabolic diseases. In this genome-wide association study (GWAS) of up to 8,293 Finns we identified 27 genome-widely significant loci (p < 1.2 × 10-9) for one or more cytokines. Fifteen of the associated variants had expression quantitative trait loci in whole blood. We provide genetic instruments to clarify the causal roles of cytokine signaling and upstream inflammation in immune-related and other chronic diseases. We further link inflammatory markers with variants previously associated with autoimmune diseases such as Crohn disease, multiple sclerosis, and ulcerative colitis and hereby elucidate the molecular mechanisms underpinning these diseases and suggest potential drug targets.

Platelets play an essential role in hemostasis and thrombosis. We performed a genome-wide association study of platelet count in 12,491 participants of the Hispanic Community Health Study/Study of Latinos by using a mixed-model method that accounts for admixture and family relationships. We discovered and replicated associations with five genes (ACTN1, ETV7, GABBR1-MOG, MEF2C, and ZBTB9-BAK1). Our strongest association was with Amerindian-specific variant rs117672662 (p value = 1.16 × 10(-28)) in ACTN1, a gene implicated in congenital macrothrombocytopenia. rs117672662 exhibited allelic differences in transcriptional activity and protein binding in hematopoietic cells. Our results underscore the value of diverse populations to extend insights into the allelic architecture of complex traits.

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10(-20)) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10(-11)) near ETS1; 3q28 (rs6444305, p = 1.10 × 10(-10)) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10(-10)) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10(-8)) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (pomnibus = 4.20 × 10(-67) to 2.67 × 10(-70)). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR(per-allele)] = 1.44; p = 4.59 × 10(-16)) and rs3130437 in HLA class I (OR(per-allele) = 1.23; p = 8.23 × 10(-9)). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.

For nephronophthisis (NPHP), the primary genetic cause of chronic renal failure in young adults, three loci have been mapped. To identify a new locus for NPHP, we here report on total-genome linkage analysis in seven families with NPHP, in whom we had excluded linkage to all three known NPHP loci. LOD scores >1 were obtained at nine loci, which were then fine mapped at 1-cM intervals. Extensive total-genome haplotype analysis revealed homozygosity in one family, in the region of the PCLN1 gene. Subsequent mutational analysis in this gene revealed PCLN1 mutations, thereby allowing exclusion of this family as a phenocopy. Multipoint linkage analysis for the remaining six families with NPHP together yielded a maximum LOD score (Z(max)) of 8.9 (at D1S253). We thus identified a new locus, NPHP4, for nephronophthisis. Markers D1S2660 and D1S2642 are flanking NPHP4 at a 2.9-cM critical interval. In one family with NPHP4, extensive genealogical studies were conducted, revealing consanguinity during the 17th century. On the basis of haplotype sharing by descent, we obtained a multipoint Z(max) of 5.8 for D1S253 in this kindred alone. In addition, we were able to localize to the NPHP4 locus a new locus for Senior-Løken syndrome, an NPHP variant associated with retinitis pigmentosa.

Age-related hearing impairment (ARHI) is the most common sensory impairment in the aging population; a third of individuals are affected by disabling hearing loss by the age of 65. It causes social isolation and depression and has recently been identified as a risk factor for dementia. The genetic risk factors and underlying pathology of ARHI are largely unknown, meaning that targets for new therapies remain elusive, yet heritability estimates range between 35% and 55%. We performed genome-wide association studies (GWASs) for two self-reported hearing phenotypes, using more than 250,000 UK Biobank (UKBB) volunteers aged between 40 and 69 years. Forty-four independent genome-wide significant loci (p < 5E-08) were identified, considerably increasing the number of established trait loci. Thirty-four loci are novel associations with hearing loss of any form, and only one of the ten known hearing loci has a previously reported association with an ARHI-related trait. Gene sets from these loci are enriched in auditory processes such as synaptic activities, nervous system processes, inner ear morphology, and cognition, while genetic correlation analysis revealed strong positive correlations with multiple personality and psychological traits for the first time. Immunohistochemistry for protein localization in adult mouse cochlea implicate metabolic, sensory, and neuronal functions for NID2, CLRN2, and ARHGEF28. These results provide insight into the genetic landscape underlying ARHI, opening up novel therapeutic targets for further investigation. In a wider context, our study also highlights the viability of using self-report phenotypes for genetic discovery in very large samples when deep phenotyping is unavailable.

Gallbladder disease (GBD) is one of the major digestive diseases. Its risk factors include age, sex, obesity, type 2 diabetes, and metabolic syndrome (MS). The prevalence of GBD is high in minority populations, such as Native and Mexican Americans. Ethnic differences, familial aggregation of GBD, and the identification of susceptibility loci for gallstone disease by use of animal models suggest genetic influences on GBD. However, the major susceptibility loci for GBD in human populations have not been identified. Using ultrasound-based information on GBD occurrence and a 10-cM gene map, we performed multipoint variance-components analysis to localize susceptibility loci for GBD. Phenotypic and genotypic data from 715 individuals in 39 low-income Mexican American families participating in the San Antonio Family Diabetes/Gallbladder Study were used. Two GBD phenotypes were defined for the analyses: (1) clinical or symptomatic GBD, the cases of cholecystectomies due to stones confirmed by ultrasound, and (2) total GBD, the clinical GBD cases plus the stone carriers newly diagnosed by ultrasound. With use of the National Cholesterol Education Program/Adult Treatment Panel III criteria, five MS risk factors were defined: increased waist circumference, hypertriglyceredemia, low high-density lipoprotein cholesterol, hypertension, and high fasting glucose. The MS risk-factor score (range 0-5) for a given individual was used as a single, composite covariate in the genetic analyses. After accounting for the effects of age, sex, and MS risk-factor score, we found stronger linkage signals for the symptomatic GBD phenotype. The highest LOD scores (3.7 and 3.5) occurred on chromosome 1p between markers D1S1597 and D1S407 (1p36.21) and near marker D1S255 (1p34.3), respectively. Other genetic locations (chromosomes 2p, 3q, 4p, 8p, 9p, 10p, and 16q) across the genome exhibited some evidence of linkage (LOD >or=1.2) to symptomatic GBD. Some of these chromosomal regions corresponded with the genetic locations of Lith loci, which influence gallstone formation in mouse models. In conclusion, we found significant evidence of major genetic determinants of symptomatic GBD on chromosome 1p in Mexican Americans.

Takayasu arteritis is a rare inflammatory disease of large arteries. We performed a genetic study in Takayasu arteritis comprising 6,670 individuals (1,226 affected individuals) from five different populations. We discovered HLA risk factors and four non-HLA susceptibility loci in VPS8, SVEP1, CFL2, and chr13q21 and reinforced IL12B, PTK2B, and chr21q22 as robust susceptibility loci shared across ancestries. Functional analysis proposed plausible underlying disease mechanisms and pinpointed ETS2 as a potential causal gene for chr21q22 association. We also identified >60 candidate loci with suggestive association (p < 5 × 10-5) and devised a genetic risk score for Takayasu arteritis. Takayasu arteritis was compared to hundreds of other traits, revealing the closest genetic relatedness to inflammatory bowel disease. Epigenetic patterns within risk loci suggest roles for monocytes and B cells in Takayasu arteritis. This work enhances understanding of the genetic basis and pathophysiology of Takayasu arteritis and provides clues for potential new therapeutic targets.

Identification of the genetic influences on human essential hypertension and other complex diseases has proved difficult, partly because of genetic heterogeneity. In many complex-trait resources, additional phenotypic data have been collected, allowing comorbid intermediary phenotypes to be used to characterize more genetically homogeneous subsets. The traditional approach to analyzing covariate-defined subsets has typically depended on researchers' previous expectations for definition of a comorbid subset and leads to smaller data sets, with a concomitant attrition in power. An alternative is to test for dependence between genetic sharing and covariates across the entire data set. This approach offers the advantage of exploiting the full data set and could be widely applied to complex-trait genome scans. However, existing maximum-likelihood methods can be prohibitively computationally expensive, especially since permutation is often required to determine significance. We developed a less computationally intensive score test and applied it to biometric and biochemical covariate data, from 2,044 sibling pairs with severe hypertension, collected by the British Genetics of Hypertension (BRIGHT) study. We found genomewide-significant evidence for linkage with hypertension and several related covariates. The strongest signals were with leaner-body-mass measures on chromosome 20q (maximum LOD = 4.24) and with parameters of renal function on chromosome 5p (maximum LOD = 3.71). After correction for the multiple traits and genetic locations studied, our global genomewide P value was .046. This is the first identity-by-descent regression analysis of hypertension to our knowledge, and it demonstrates the value of this approach for the incorporation of additional phenotypic information in genetic studies of complex traits.

Complex traits important for humans are often correlated phenotypically and genetically. Joint mapping of quantitative-trait loci (QTLs) for multiple correlated traits plays an important role in unraveling the genetic architecture of complex traits. Compared with single-trait analysis, joint mapping addresses more questions and has advantages for power of QTL detection and precision of parameter estimation. Some statistical methods have been developed to map QTLs underlying multiple traits, most of which are based on maximum-likelihood methods. We develop here a multivariate version of the Bayes methodology for joint mapping of QTLs, using the Markov chain-Monte Carlo (MCMC) algorithm. We adopt a variance-components method to model complex traits in outbred populations (e.g., humans). The method is robust, can deal with an arbitrary number of alleles with arbitrary patterns of gene actions (such as additive and dominant), and allows for multiple phenotype data of various types in the joint analysis (e.g., multiple continuous traits and mixtures of continuous traits and discrete traits). Under a Bayesian framework, parameters--including the number of QTLs--are estimated on the basis of their marginal posterior samples, which are generated through two samplers, the Gibbs sampler and the reversible-jump MCMC. In addition, we calculate the Bayes factor related to each identified QTL, to test coincident linkage versus pleiotropy. The performance of our method is evaluated in simulations with full-sib families. The results show that our proposed Bayesian joint-mapping method performs well for mapping multiple QTLs in situations of either bivariate continuous traits or mixed data types. Compared with the analysis for each trait separately, Bayesian joint mapping improves statistical power, provides stronger evidence of QTL detection, and increases precision in estimation of parameter and QTL position. We also applied the proposed method to a set of real data and detected a coincident linkage responsible for determining bone mineral density and areal bone size of wrist in humans.

To date, 68 loci have been associated with type 2 diabetes (T2D) or glucose homeostasis traits. We report here the results of experiments aimed at functionally characterizing the SNPs replicated for T2D and glucose traits. We sought to determine whether these loci were associated with transcript levels in adipose, muscle, liver, lymphocytes, and pancreatic β-cells. We found an excess of trans, rather than cis, associations among these SNPs in comparison to what was expected in adipose and muscle. Among transcripts differentially expressed (FDR < 0.05) between muscle or adipose cells of insulin-sensitive individuals and those of insulin-resistant individuals (matched on BMI), trans-regulated transcripts, in contrast to the cis-regulated ones, were enriched. The paucity of cis associations with transcripts was confirmed in a study of liver transcriptome and was further supported by an analysis of the most detailed transcriptome map of pancreatic β-cells. Relative to location- and allele-frequency-matched random SNPs, both the 68 loci and top T2D-associated SNPs from two large-scale genome-wide studies were enriched for trans eQTLs in adipose and muscle but not in lymphocytes. Our study suggests that T2D SNPs have broad-reaching and tissue-specific effects that often extend beyond local transcripts and raises the question of whether patterns of cis or trans transcript regulation are a key feature of the architecture of complex traits.

Identifying genomic annotations that differentiate causal from trait-associated variants is essential to fine mapping disease loci. Although many studies have identified non-coding functional annotations that overlap disease-associated variants, these annotations often colocalize, complicating the ability to use these annotations for fine mapping causal variation. We developed a statistical approach (Genomic Annotation Shifter [GoShifter]) to assess whether enriched annotations are able to prioritize causal variation. GoShifter defines the null distribution of an annotation overlapping an allele by locally shifting annotations; this approach is less sensitive to biases arising from local genomic structure than commonly used enrichment methods that depend on SNP matching. Local shifting also allows GoShifter to identify independent causal effects from colocalizing annotations. Using GoShifter, we confirmed that variants in expression quantitative trail loci drive gene-expression changes though DNase-I hypersensitive sites (DHSs) near transcription start sites and independently through 3' UTR regulation. We also showed that (1) 15%-36% of trait-associated loci map to DHSs independently of other annotations; (2) loci associated with breast cancer and rheumatoid arthritis harbor potentially causal variants near the summits of histone marks rather than full peak bodies; (3) variants associated with height are highly enriched in embryonic stem cell DHSs; and (4) we can effectively prioritize causal variation at specific loci.

Blood pressure (BP) is a heritable risk factor for cardiovascular disease. To investigate genetic associations with systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP), we genotyped ~50,000 SNPs in up to 87,736 individuals of European ancestry and combined these in a meta-analysis. We replicated findings in an independent set of 68,368 individuals of European ancestry. Our analyses identified 11 previously undescribed associations in independent loci containing 31 genes including PDE1A, HLA-DQB1, CDK6, PRKAG2, VCL, H19, NUCB2, RELA, HOXC@ complex, FBN1, and NFAT5 at the Bonferroni-corrected array-wide significance threshold (p < 6 × 10(-7)) and confirmed 27 previously reported associations. Bioinformatic analysis of the 11 loci provided support for a putative role in hypertension of several genes, such as CDK6 and NUCB2. Analysis of potential pharmacological targets in databases of small molecules showed that ten of the genes are predicted to be a target for small molecules. In summary, we identified previously unknown loci associated with BP. Our findings extend our understanding of genes involved in BP regulation, which may provide new targets for therapeutic intervention or drug response stratification.

Venous thromboembolism (VTE), the third leading cause of cardiovascular mortality, is a complex thrombotic disorder with environmental and genetic determinants. Although several genetic variants have been found associated with VTE, they explain a minor proportion of VTE risk in cases. We undertook a meta-analysis of genome-wide association studies (GWASs) to identify additional VTE susceptibility genes. Twelve GWASs totaling 7,507 VTE case subjects and 52,632 control subjects formed our discovery stage where 6,751,884 SNPs were tested for association with VTE. Nine loci reached the genome-wide significance level of 5 × 10(-8) including six already known to associate with VTE (ABO, F2, F5, F11, FGG, and PROCR) and three unsuspected loci. SNPs mapping to these latter were selected for replication in three independent case-control studies totaling 3,009 VTE-affected individuals and 2,586 control subjects. This strategy led to the identification and replication of two VTE-associated loci, TSPAN15 and SLC44A2, with lead risk alleles associated with odds ratio for disease of 1.31 (p = 1.67 × 10(-16)) and 1.21 (p = 2.75 × 10(-15)), respectively. The lead SNP at the TSPAN15 locus is the intronic rs78707713 and the lead SLC44A2 SNP is the non-synonymous rs2288904 previously shown to associate with transfusion-related acute lung injury. We further showed that these two variants did not associate with known hemostatic plasma markers. TSPAN15 and SLC44A2 do not belong to conventional pathways for thrombosis and have not been associated to other cardiovascular diseases nor related quantitative biomarkers. Our findings uncovered unexpected actors of VTE etiology and pave the way for novel mechanistic concepts of VTE pathophysiology.

Resolving the molecular processes that mediate genetic risk remains a challenge because most disease-associated variants are non-coding and functional characterization of these signals requires knowledge of the specific tissues and cell-types in which they operate. To address this challenge, we developed a framework for integrating tissue-specific gene expression and epigenomic maps to obtain "tissue-of-action" (TOA) scores for each association signal by systematically partitioning posterior probabilities from Bayesian fine-mapping. We applied this scheme to credible set variants for 380 association signals from a recent GWAS meta-analysis of type 2 diabetes (T2D) in Europeans. The resulting tissue profiles underscored a predominant role for pancreatic islets and, to a lesser extent, adipose and liver, particularly among signals with greater fine-mapping resolution. We incorporated resulting TOA scores into a rule-based classifier and validated the tissue assignments through comparison with data from cis-eQTL enrichment, functional fine-mapping, RNA co-expression, and patterns of physiological association. In addition to implicating signals with a single TOA, we found evidence for signals with shared effects in multiple tissues as well as distinct tissue profiles between independent signals within heterogeneous loci. Lastly, we demonstrated that TOA scores can be directly coupled with eQTL colocalization to further resolve effector transcripts at T2D signals. This framework guides mechanistic inference by directing functional validation studies to the most relevant tissues and can gain power as fine-mapping resolution and cell-specific annotations become richer. This method is generalizable to all complex traits with relevant annotation data and is made available as an R package.

Risk prediction based on genomic profiles has raised a lot of attention recently. However, family history is usually ignored in genetic risk prediction. In this study we proposed a statistical framework for risk prediction given an individual's genotype profile and family history. Genotype information about the relatives can also be incorporated. We allow risk prediction given the current age and follow-up period and consider competing risks of mortality. The framework allows easy extension to any family size and structure. In addition, the predicted risk at any percentile and the risk distribution graphs can be computed analytically. We applied the method to risk prediction for breast and prostate cancers by using known susceptibility loci from genome-wide association studies. For breast cancer, in the population the 10-year risk at age 50 ranged from 1.1% at the 5th percentile to 4.7% at the 95th percentile. If we consider the average 10-year risk at age 50 (2.39%) as the threshold for screening, the screening age ranged from 62 at the 20th percentile to 38 at the 95th percentile (and some never reach the threshold). For women with one affected first-degree relative, the 10-year risks ranged from 2.6% (at the 5th percentile) to 8.1% (at the 95th percentile). For prostate cancer, the corresponding 10-year risks at age 60 varied from 1.8% to 14.9% in the population and from 4.2% to 23.2% in those with an affected first-degree relative. We suggest that for some diseases genetic testing that incorporates family history can stratify people into diverse risk categories and might be useful in targeted prevention and screening.

Epidemiological studies have demonstrated that genetic factors account for at least 50% of the liability for nicotine dependence (ND). Although several linkage studies have been conducted, all samples to date were primarily of European origin. In this study, we conducted a genomewide scan of 1,261 individuals, representing 402 nuclear families, of African American (AA) origin. We examined 385 autosomal microsatellite markers for ND, which was assessed by smoking quantity (SQ), the Heaviness of Smoking Index (HSI), and the Fagerstrom Test for ND (FTND). After performing linkage analyses using various methods implemented in the GENEHUNTER and S.A.G.E. programs, we found a region near marker D10S1432 on chromosome 10q22 that showed a significant linkage to indexed SQ, with a maximum LOD score of 4.17 at 92 cM and suggestive linkage to HSI, SQ, and log-transformed SQ. Additionally, we identified three regions that met the criteria for suggestive linkage to at least one ND measure: on chromosomes 9q31 at marker D9S1825, 11p11 between markers D11S1993 and D11S1344, and 13q13 between markers D13S325 and D13S788. Other locations on chromosomes 15p11, 17q25, and 18q12 exhibited some evidence of linkage for ND (LOD >1.44). The four regions with significant or suggestive linkage were positive for multiple ND measures by multiple statistical methods. Some of these regions have been linked to smoking behavior at nominally significant levels in other studies, which provides independent replication of the regions for ND in different cohorts. In summary, we found significant linkage on chromosome 10q22 and suggestive linkage on chromosomes 9, 11, and 13 for major genetic determinants of ND in an AA sample. Further analysis of these positive regions by fine mapping and/or association analysis is thus warranted. To our knowledge, this study represents the first genomewide linkage scan of ND in an AA sample.

Alternate splicing events can create isoforms that alter gene function, and genetic variants associated with alternate gene isoforms may reveal molecular mechanisms of disease. We used subcutaneous adipose tissue of 426 Finnish men from the METSIM study and identified splice junction quantitative trait loci (sQTLs) for 6,077 splice junctions (FDR < 1%). In the same individuals, we detected expression QTLs (eQTLs) for 59,443 exons and 15,397 genes (FDR < 1%). We identified 595 genes with an sQTL and exon eQTL but no gene eQTL, which could indicate potential isoform differences. Of the significant sQTL signals, 2,114 (39.8%) included at least one proxy variant (linkage disequilibrium r2 > 0.8) located within an intron spanned by the splice junction. We identified 203 sQTLs that colocalized with 141 genome-wide association study (GWAS) signals for cardiometabolic traits, including 25 signals for lipid traits, 24 signals for body mass index (BMI), and 12 signals for waist-hip ratio adjusted for BMI. Among all 141 GWAS signals colocalized with an sQTL, we detected 26 that also colocalized with an exon eQTL for an exon skipped by the sQTL splice junction. At a GWAS signal for high-density lipoprotein cholesterol colocalized with an NR1H3 sQTL splice junction, we show that the alternative splice product encodes an NR1H3 transcription factor that lacks a DNA binding domain and fails to activate transcription. Together, these results detect splicing events and candidate mechanisms that may contribute to gene function at GWAS loci.