Recent work has shown that much of the missing heritability of complex traits can be resolved by estimates of heritability explained by all genotyped SNPs. However, it is currently unknown how much heritability is missing due to poor tagging or additional causal variants at known GWAS loci. Here, we use variance components to quantify the heritability explained by all SNPs at known GWAS loci in nine diseases from WTCCC1 and WTCCC2. After accounting for expectation, we observed all SNPs at known GWAS loci to explain 1.29 x more heritability than GWAS-associated SNPs on average (P=3.3 x 10⁻⁵). For some diseases, this increase was individually significant: 2.07 x for Multiple Sclerosis (MS) (P=6.5 x 10⁻⁹) and 1.48 x for Crohns Disease (CD) (P = 1.3 x 10⁻³); all analyses of autoimmune diseases excluded the well-studied MHC region. Additionally, we found that GWAS loci from other related traits also explained significant heritability. The union of all autoimmune disease loci explained 7.15 x more MS heritability than known MS SNPs (P < 1.0 x 10⁻¹⁶ and 2.20 x more CD heritability than known CD SNPs (P = 6.1 x 10⁻⁹), with an analogous increase for all autoimmune diseases analyzed. We also observed significant increases in an analysis of > 20,000 Rheumatoid Arthritis (RA) samples typed on ImmunoChip, with 2.37 x more heritability from all SNPs at GWAS loci (P = 2.3 x 10⁻⁶) and 5.33 x more heritability from all autoimmune disease loci (P < 1 x 10⁻¹⁶ compared to known RA SNPs (including those identified in this cohort). Our methods adjust for LD between SNPs, which can bias standard estimates of heritability from SNPs even if all causal variants are typed. By comparing adjusted estimates, we hypothesize that the genome-wide distribution of causal variants is enriched for low-frequency alleles, but that causal variants at known GWAS loci are skewed towards common alleles. These findings have important ramifications for fine-mapping study design and our understanding of complex disease architecture.

The pharmacokinetic and pharmacodynamic disciplines address pharmacological traits, including efficacy and adverse events. Pharmacogenomics studies have identified pervasive genetic effects on treatment outcomes, resulting in the development of genetic biomarkers for optimization of drug therapy. Pharmacogenomics-based tests are already being applied in clinical decision making. However, despite substantial progress in identifying the genetic etiology of pharmacological response, current biomarker panels still largely rely on single gene tests with a large portion of the genetic effects remaining to be discovered. Future research must account for the combined effects of multiple genetic variants, incorporate pathway-based approaches, explore gene-gene interactions and nonprotein coding functional genetic variants, extend studies across ancestral populations, and prioritize laboratory characterization of molecular mechanisms. Because genetic factors can play a key role in drug response, accurate biomarker tests capturing the main genetic factors determining treatment outcomes have substantial potential for improving individual clinical care.

Genome-wide association studies (GWASs) have focused primarily on populations of European descent, but it is essential that diverse populations become better represented. Increasing diversity among study participants will advance our understanding of genetic architecture in all populations and ensure that genetic research is broadly applicable. To facilitate and promote research in multi-ancestry and admixed cohorts, we outline key methodological considerations and highlight opportunities, challenges, solutions, and areas in need of development. Despite the perception that analyzing genetic data from diverse populations is difficult, it is scientifically and ethically imperative, and there is an expanding analytical toolbox to do it well.

Transcriptomic imputation approaches combine eQTL reference panels with large-scale genotype data in order to test associations between disease and gene expression. These genic associations could elucidate signals in complex genome-wide association study (GWAS) loci and may disentangle the role of different tissues in disease development. We used the largest eQTL reference panel for the dorso-lateral prefrontal cortex (DLPFC) to create a set of gene expression predictors and demonstrate their utility. We applied DLPFC and 12 GTEx-brain predictors to 40,299 schizophrenia cases and 65,264 matched controls for a large transcriptomic imputation study of schizophrenia. We identified 413 genic associations across 13 brain regions. Stepwise conditioning identified 67 non-MHC genes, of which 14 did not fall within previous GWAS loci. We identified 36 significantly enriched pathways, including hexosaminidase-A deficiency, and multiple porphyric disorder pathways. We investigated developmental expression patterns among the 67 non-MHC genes and identified specific groups of pre- and postnatal expression.

In the HTML version of the article originally published, the author group 'The Schizophrenia Working Group of the Psychiatric Genomics Consortium' was displayed incorrectly. The error has been corrected in the HTML version of the article.

Genetic variants underlying complex traits, including disease susceptibility, are enriched within the transcriptional regulatory elements, promoters and enhancers. There is emerging evidence that regulatory elements associated with particular traits or diseases share similar patterns of transcriptional activity. Accordingly, shared transcriptional activity (coexpression) may help prioritise loci associated with a given trait, and help to identify underlying biological processes. Using cap analysis of gene expression (CAGE) profiles of promoter- and enhancer-derived RNAs across 1824 human samples, we have analysed coexpression of RNAs originating from trait-associated regulatory regions using a novel quantitative method (network density analysis; NDA). For most traits studied, phenotype-associated variants in regulatory regions were linked to tightly-coexpressed networks that are likely to share important functional characteristics. Coexpression provides a new signal, independent of phenotype association, to enable fine mapping of causative variants. The NDA coexpression approach identifies new genetic variants associated with specific traits, including an association between the regulation of the OCT1 cation transporter and genetic variants underlying circulating cholesterol levels. NDA strongly implicates particular cell types and tissues in disease pathogenesis. For example, distinct groupings of disease-associated regulatory regions implicate two distinct biological processes in the pathogenesis of ulcerative colitis; a further two separate processes are implicated in Crohn's disease. Thus, our functional analysis of genetic predisposition to disease defines new distinct disease endotypes. We predict that patients with a preponderance of susceptibility variants in each group are likely to respond differently to pharmacological therapy. Together, these findings enable a deeper biological understanding of the causal basis of complex traits.

Background

Treatment strategies blocking tumour necrosis factor (anti-TNF) have proven very successful in patients with rheumatoid arthritis (RA). However, a significant subset of patients does not respond for unknown reasons. Currently, there are no means of identifying these patients before treatment. This study was aimed at identifying genetic factors predicting anti-TNF treatment outcome in patients with RA using a genome-wide association approach.

Methods

We conducted a multistage, genome-wide association study with a primary analysis of 2 557 253 single-nucleotide polymorphisms (SNPs) in 882 patients with RA receiving anti-TNF therapy included through the Dutch Rheumatoid Arthritis Monitoring (DREAM) registry and the database of Apotheekzorg. Linear regression analysis of changes in the Disease Activity Score in 28 joints after 14 weeks of treatment was performed using an additive model. Markers with p<10(-3) were selected for replication in 1821 patients from three independent cohorts. Pathway analysis including all SNPs with p<10(-3) was performed using Ingenuity.

Results

772 markers showed evidence of association with treatment outcome in the initial stage. Eight genetic loci showed improved p value in the overall meta-analysis compared with the first stage, three of which (rs1568885, rs1813443 and rs4411591) showed directional consistency over all four cohorts studied. We were unable to replicate markers previously reported to be associated with anti-TNF outcome. Network analysis indicated strong involvement of biological processes underlying inflammatory response and cell morphology.

Conclusions

Using a multistage strategy, we have identified eight genetic loci associated with response to anti-TNF treatment. Further studies are required to validate these findings in additional patient collections.

Genome-wide association studies have identified over 150 loci associated with lipid traits, however, no large-scale studies exist for Hispanics and other minority populations. Additionally, the genetic architecture of lipid-influencing loci remains largely unknown. We performed one of the most racially/ethnically diverse fine-mapping genetic studies of HDL-C, LDL-C, and triglycerides to-date using SNPs on the MetaboChip array on 54,119 individuals: 21,304 African Americans, 19,829 Hispanic Americans, 12,456 Asians, and 530 American Indians. The majority of signals found in these groups generalize to European Americans. While we uncovered signals unique to racial/ethnic populations, we also observed systematically consistent lipid associations across these groups. In African Americans, we identified three novel signals associated with HDL-C (LPL, APOA5, LCAT) and two associated with LDL-C (ABCG8, DHODH). In addition, using this population, we refined the location for 16 out of the 58 known MetaboChip lipid loci. These results can guide tailored screening efforts, reveal population-specific responses to lipid-lowering medications, and aid in the development of new targeted drug therapies.

Susceptibility to schizophrenia is inversely correlated with general cognitive ability at both the phenotypic and the genetic level. Paradoxically, a modest but consistent positive genetic correlation has been reported between schizophrenia and educational attainment, despite the strong positive genetic correlation between cognitive ability and educational attainment. Here we leverage published genome-wide association studies (GWASs) in cognitive ability, education, and schizophrenia to parse biological mechanisms underlying these results. Association analysis based on subsets (ASSET), a pleiotropic meta-analytic technique, allowed jointly associated loci to be identified and characterized. Specifically, we identified subsets of variants associated in the expected ("concordant") direction across all three phenotypes (i.e., greater risk for schizophrenia, lower cognitive ability, and lower educational attainment); these were contrasted with variants that demonstrated the counterintuitive ("discordant") relationship between education and schizophrenia (i.e., greater risk for schizophrenia and higher educational attainment). ASSET analysis revealed 235 independent loci associated with cognitive ability, education, and/or schizophrenia at p < 5 × 10^-8. Pleiotropic analysis successfully identified more than 100 loci that were not significant in the input GWASs. Many of these have been validated by larger, more recent single-phenotype GWASs. Leveraging the joint genetic correlations of cognitive ability, education, and schizophrenia, we were able to dissociate two distinct biological mechanisms-early neurodevelopmental pathways that characterize concordant allelic variation and adulthood synaptic pruning pathways-that were linked to the paradoxical positive genetic association between education and schizophrenia. Furthermore, genetic correlation analyses revealed that these mechanisms contribute not only to the etiopathogenesis of schizophrenia but also to the broader biological dimensions implicated in both general health outcomes and psychiatric illness.

To reveal post-traumatic stress disorder (PTSD) genetic risk influences on tissue-specific gene expression, we use brain and non-brain transcriptomic imputation. We impute genetically regulated gene expression (GReX) in 29,539 PTSD cases and 166,145 controls from 70 ancestry-specific cohorts and identify 18 significant GReX-PTSD associations corresponding to specific tissue-gene pairs. The results suggest substantial genetic heterogeneity based on ancestry, cohort type (military versus civilian), and sex. Two study-wide significant PTSD associations are identified in European and military European cohorts; ZNF140 is predicted to be upregulated in whole blood, and SNRNP35 is predicted to be downregulated in dorsolateral prefrontal cortex, respectively. In peripheral leukocytes from 175 marines, the observed PTSD differential gene expression correlates with the predicted differences for these individuals, and deployment stress produces glucocorticoid-regulated expression changes that include downregulation of both ZNF140 and SNRNP35. SNRNP35 knockdown in cells validates its functional role in U12-intron splicing. Finally, exogenous glucocorticoids in mice downregulate prefrontal Snrnp35 expression.