Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape.

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

Factorizing polygenic epistasis improves prediction and uncovers biological pathways in complex traits.

Epistasis is central in many domains of biology, but it has not yet been proven useful for understanding the etiology of complex traits. This is partly because complex-trait epistasis involves polygenic interactions that are poorly captured in current models. To address this gap, we developed a model called Epistasis Factor Analysis (EFA). EFA assumes that polygenic epistasis can be factorized into interactions between a few epistasis factors (EFs), which represent latent polygenic components of the observed complex trait. The statistical goals of EFA are to improve polygenic prediction and to increase power to detect epistasis, while the biological goal is to unravel genetic effects into more-homogeneous units. We mathematically characterize EFA and use simulations to show that EFA outperforms current epistasis models when its assumptions approximately hold. Applied to predicting yeast growth rates, EFA outperforms the additive model for several traits with large epistasis heritability and uniformly outperforms the standard epistasis model. We replicate these prediction improvements in a second dataset. We then apply EFA to four previously characterized traits in the UK Biobank and find statistically significant epistasis in all four, including two that are robust to scale transformation. Moreover, we find that the inferred EFs partly recover pre-defined biological pathways for two of the traits. Our results demonstrate that more realistic models can identify biologically and statistically meaningful epistasis in complex traits, indicating that epistasis has potential for precision medicine and characterizing the biology underlying GWAS results.

Impact of cross-ancestry genetic architecture on GWASs in admixed populations.

Genome-wide association studies (GWASs) have identified thousands of variants for disease risk. These studies have predominantly been conducted in individuals of European ancestries, which raises questions about their transferability to individuals of other ancestries. Of particular interest are admixed populations, usually defined as populations with recent ancestry from two or more continental sources. Admixed genomes contain segments of distinct ancestries that vary in composition across individuals in the population, allowing for the same allele to induce risk for disease on different ancestral backgrounds. This mosaicism raises unique challenges for GWASs in admixed populations, such as the need to correctly adjust for population stratification. In this work we quantify the impact of differences in estimated allelic effect sizes for risk variants between ancestry backgrounds on association statistics. Specifically, while the possibility of estimated allelic effect-size heterogeneity by ancestry (HetLanc) can be modeled when performing a GWAS in admixed populations, the extent of HetLanc needed to overcome the penalty from an additional degree of freedom in the association statistic has not been thoroughly quantified. Using extensive simulations of admixed genotypes and phenotypes, we find that controlling for and conditioning effect sizes on local ancestry can reduce statistical power by up to 72%. This finding is especially pronounced in the presence of allele frequency differentiation. We replicate simulation results using 4,327 African-European admixed genomes from the UK Biobank for 12 traits to find that for most significant SNPs, HetLanc is not large enough for GWASs to benefit from modeling heterogeneity in this way.

Harnessing tissue-specific genetic variation to dissect putative causal pathways between body mass index and cardiometabolic phenotypes.

Body mass index (BMI) is a complex disease risk factor known to be influenced by genes acting via both metabolic pathways and appetite regulation. In this study, we aimed to gain insight into the phenotypic consequences of BMI-associated genetic variants, which may be mediated by their expression in different tissues. First, we harnessed meta-analyzed gene expression datasets derived from subcutaneous adipose (n = 1257) and brain (n = 1194) tissue to identify 86 and 140 loci, respectively, which provided evidence of genetic colocalization with BMI. These two sets of tissue-partitioned loci had differential effects with respect to waist-to-hip ratio, suggesting that the way they influence fat distribution might vary despite their having very similar average magnitudes of effect on BMI itself (adipose = 0.0148 and brain = 0.0149 standard deviation change in BMI per effect allele). For instance, BMI-associated variants colocalized with TBX15 expression in adipose tissue (posterior probability [PPA] = 0.97), but not when we used TBX15 expression data derived from brain tissue (PPA = 0.04) This gene putatively influences BMI via its role in skeletal development. Conversely, there were loci where BMI-associated variants provided evidence of colocalization with gene expression in brain tissue (e.g., NEGR1, PPA = 0.93), but not when we used data derived from adipose tissue, suggesting that these genes might be more likely to influence BMI via energy balance. Leveraging these tissue-partitioned variant sets through a multivariable Mendelian randomization framework provided strong evidence that the brain-tissue-derived variants are predominantly responsible for driving the genetically predicted effects of BMI on cardiovascular-disease endpoints (e.g., coronary artery disease: odds ratio = 1.05, 95% confidence interval = 1.04-1.07, p = 4.67 × 10-14). In contrast, our analyses suggested that the adipose tissue variants might predominantly be responsible for the underlying relationship between BMI and measures of cardiac function, such as left ventricular stroke volume (beta = 0.21, 95% confidence interval = 0.09-0.32, p = 6.43 × 10-4).

Mendelian pathway analysis of laboratory traits reveals distinct roles for ciliary subcompartments in common disease pathogenesis.

Rare monogenic disorders of the primary cilium, termed ciliopathies, are characterized by extreme presentations of otherwise common diseases, such as diabetes, hepatic fibrosis, and kidney failure. However, despite a recent revolution in our understanding of the cilium's role in rare disease pathogenesis, the organelle's contribution to common disease remains largely unknown. Hypothesizing that common genetic variants within Mendelian ciliopathy genes might contribute to common complex diseases pathogenesis, we performed association studies of 16,874 common genetic variants across 122 ciliary genes with 12 quantitative laboratory traits characteristic of ciliopathy syndromes in 452,593 individuals in the UK Biobank. We incorporated tissue-specific gene expression analysis, expression quantitative trait loci, and Mendelian disease phenotype information into our analysis and replicated our findings in meta-analysis. 101 statistically significant associations were identified across 42 of the 122 examined ciliary genes (including eight novel replicating associations). These ciliary genes were widely expressed in tissues relevant to the phenotypes being studied, and eQTL analysis revealed strong evidence for correlation between ciliary gene expression levels and laboratory traits. Perhaps most interestingly, our analysis identified different ciliary subcompartments as being specifically associated with distinct sets of phenotypes. Taken together, our data demonstrate the utility of a Mendelian pathway-based approach to genomic association studies, challenge the widely held belief that the cilium is an organelle important mainly in development and in rare syndromic disease pathogenesis, and provide a framework for the continued integration of common and rare disease genetics to provide insight into the pathophysiology of human diseases of immense public health burden.

PGS-server: accuracy, robustness and transferability of polygenic score methods for biobank scale studies.

Polygenic scores (PGS) are important tools for carrying out genetic prediction of common diseases and disease related complex traits, facilitating the development of precision medicine. Unfortunately, despite the critical importance of PGS and the vast number of PGS methods recently developed, few comprehensive comparison studies have been performed to evaluate the effectiveness of PGS methods. To fill this critical knowledge gap, we performed a comprehensive comparison study on 12 different PGS methods through internal evaluations on 25 quantitative and 25 binary traits within the UK Biobank with sample sizes ranging from 147 408 to 336 573, and through external evaluations via 25 cross-study and 112 cross-ancestry analyses on summary statistics from multiple genome-wide association studies with sample sizes ranging from 1415 to 329 345. We evaluate the prediction accuracy, computational scalability, as well as robustness and transferability of different PGS methods across datasets and/or genetic ancestries, providing important guidelines for practitioners in choosing PGS methods. Besides method comparison, we present a simple aggregation strategy that combines multiple PGS from different methods to take advantage of their distinct benefits to achieve stable and superior prediction performance. To facilitate future applications of PGS, we also develop a PGS webserver (http://www.pgs-server.com/) that allows users to upload summary statistics and choose different PGS methods to fit the data directly. We hope that our results, method and webserver will facilitate the routine application of PGS across different research areas.

Changes in environmental exposures over decades may influence the genetic architecture of severe spermatogenic failure.

STUDY QUESTION: Do the genetic determinants of idiopathic severe spermatogenic failure (SPGF) differ between generations? SUMMARY ANSWER: Our data support that the genetic component of idiopathic SPGF is impacted by dynamic changes in environmental exposures over decades. WHAT IS KNOWN ALREADY: The idiopathic form of SPGF has a multifactorial etiology wherein an interaction between genetic, epigenetic, and environmental factors leads to the disease onset and progression. At the genetic level, genome-wide association studies (GWASs) allow the analysis of millions of genetic variants across the genome in a hypothesis-free manner, as a valuable tool for identifying susceptibility risk loci. However, little is known about the specific role of non-genetic factors and their influence on the genetic determinants in this type of conditions. STUDY DESIGN, SIZE, DURATION: Case-control genetic association analyses were performed including a total of 912 SPGF cases and 1360 unaffected controls. PARTICIPANTS/MATERIALS, SETTING, METHODS: All participants had European ancestry (Iberian and German). SPGF cases were diagnosed during the last decade either with idiopathic non-obstructive azoospermia (n = 547) or with idiopathic non-obstructive oligozoospermia (n = 365). Case-control genetic association analyses were performed by logistic regression models considering the generation as a covariate and by in silico functional characterization of the susceptibility genomic regions. MAIN RESULTS AND THE ROLE OF CHANCE: This analysis revealed 13 novel genetic association signals with SPGF, with eight of them being independent. The observed associations were mostly explained by the interaction between each lead variant and the age-group. Additionally, we established links between these loci and diverse non-genetic factors, such as toxic or dietary habits, respiratory disorders, and autoimmune diseases, which might potentially influence the genetic architecture of idiopathic SPGF. LARGE SCALE DATA: GWAS data are available from the authors upon reasonable request. LIMITATIONS, REASONS FOR CAUTION: Additional independent studies involving large cohorts in ethnically diverse populations are warranted to confirm our findings. WIDER IMPLICATIONS OF THE FINDINGS: Overall, this study proposes an innovative strategy to achieve a more precise understanding of conditions such as SPGF by considering the interactions between a variable exposome through different generations and genetic predisposition to complex diseases. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by the "Plan Andaluz de Investigación, Desarrollo e Innovación (PAIDI 2020)" (ref. PY20_00212, P20_00583), the Spanish Ministry of Economy and Competitiveness through the Spanish National Plan for Scientific and Technical Research and Innovation (ref. PID2020-120157RB-I00 funded by MCIN/ AEI/10.13039/501100011033), and the 'Proyectos I+D+i del Programa Operativo FEDER 2020' (ref. B-CTS-584-UGR20). ToxOmics-Centre for Toxicogenomics and Human Health, Genetics, Oncology and Human Toxicology, is also partially supported by the Portuguese Foundation for Science and Technology (Projects: UIDB/00009/2020; UIDP/00009/2020). The authors declare no competing interests. TRIAL REGISTRATION NUMBER: N/A.

Cis-meQTL for cocaine use-associated DNA methylation in an HIV-positive cohort show pleiotropic effects on multiple traits.

BACKGROUND: Cocaine use (CU) is associated with psychiatric and medical diseases. Little is known about the mechanisms of CU-related comorbidities. Findings from preclinical and clinical studies have suggested that CU is associated with aberrant DNA methylation (DNAm) that may be influenced by genetic variants [i.e., methylation quantitative trait loci (meQTLs)]. In this study, we mapped cis-meQTLs for CU-associated DNAm sites (CpGs) in an HIV-positive cohort (Ntotal = 811) and extended the meQTLs to multiple traits. RESULTS: We conducted cis-meQTL analysis for 224 candidate CpGs selected for their association with CU in blood. We identified 7,101 significant meQTLs [false discovery rate (FDR) < 0.05], which mostly mapped to genes involved in immunological functions and were enriched in immune pathways. We followed up the meQTLs using phenome-wide association study and trait enrichment analyses, which revealed 9 significant traits. We tested for causal effects of CU on these 9 traits using Mendelian Randomization and found evidence that CU plays a causal role in increasing hypertension (p-value = 2.35E-08) and decreasing heel bone mineral density (p-value = 1.92E-19). CONCLUSIONS: These findings suggest that genetic variants for CU-associated DNAm have pleiotropic effects on other relevant traits and provide new insights into the causal relationships between cocaine use and these complex traits.

Complete Sequence of the 22q11.2 Allele in 1,053 Subjects with 22q11.2 Deletion Syndrome Reveals Modifiers of Conotruncal Heart Defects.

The 22q11.2 deletion syndrome (22q11.2DS) results from non-allelic homologous recombination between low-copy repeats termed LCR22. About 60%-70% of individuals with the typical 3 megabase (Mb) deletion from LCR22A-D have congenital heart disease, mostly of the conotruncal type (CTD), whereas others have normal cardiac anatomy. In this study, we tested whether variants in the hemizygous LCR22A-D region are associated with risk for CTDs on the basis of the sequence of the 22q11.2 region from 1,053 22q11.2DS individuals. We found a significant association (FDR p < 0.05) of the CTD subset with 62 common variants in a single linkage disequilibrium (LD) block in a 350 kb interval harboring CRKL. A total of 45 of the 62 variants were associated with increased risk for CTDs (odds ratio [OR) ranges: 1.64-4.75). Associations of four variants were replicated in a meta-analysis of three genome-wide association studies of CTDs in affected individuals without 22q11.2DS. One of the replicated variants, rs178252, is located in an open chromatin region and resides in the double-elite enhancer, GH22J020947, that is predicted to regulate CRKL (CRK-like proto-oncogene, cytoplasmic adaptor) expression. Approximately 23% of patients with nested LCR22C-D deletions have CTDs, and inactivation of Crkl in mice causes CTDs, thus implicating this gene as a modifier. Rs178252 and rs6004160 are expression quantitative trait loci (eQTLs) of CRKL. Furthermore, set-based tests identified an enhancer that is predicted to target CRKL and is significantly associated with CTD risk (GH22J020946, sequence kernal association test (SKAT) p = 7.21 × 10-5) in the 22q11.2DS cohort. These findings suggest that variance in CTD penetrance in the 22q11.2DS population can be explained in part by variants affecting CRKL expression.

A Bayesian method for estimating gene-level polygenicity under the framework of transcriptome-wide association study.

Polygenicity refers to the phenomenon that multiple genetic variants have a nonzero effect on a complex trait. It is defined as the proportion of genetic variants with a nonzero effect on the trait. Evaluation of polygenicity can provide valuable insights into the genetic architecture of the trait. Several recent works have attempted to estimate polygenicity at the single nucleotide polymorphism level. However, evaluating polygenicity at the gene level can be biologically more meaningful. We propose the notion of gene-level polygenicity, defined as the proportion of genes having a nonzero effect on the trait under the framework of a transcriptome-wide association study. We introduce a Bayesian approach genepoly to estimate this quantity for a trait. The method is based on spike and slab prior and simultaneously estimates the subset of non-null genes. Our simulation study shows that genepoly efficiently estimates gene-level polygenicity. The method produces a downward bias for small choices of trait heritability due to a non-null gene, which diminishes rapidly with an increase in the genome-wide association study (GWAS) sample size. While identifying the subset of non-null genes, genepoly offers a high level of specificity and an overall good level of sensitivity-the sensitivity increases as the sample size of the reference panel expression and GWAS data increase. We applied the method to seven phenotypes in the UK Biobank, integrating expression data. We find height to be the most polygenic and asthma to be the least polygenic.

Comparing Within- and Between-Family Polygenic Score Prediction.

Polygenic scores are a popular tool for prediction of complex traits. However, prediction estimates in samples of unrelated participants can include effects of population stratification, assortative mating, and environmentally mediated parental genetic effects, a form of genotype-environment correlation (rGE). Comparing genome-wide polygenic score (GPS) predictions in unrelated individuals with predictions between siblings in a within-family design is a powerful approach to identify these different sources of prediction. Here, we compared within- to between-family GPS predictions of eight outcomes (anthropometric, cognitive, personality, and health) for eight corresponding GPSs. The outcomes were assessed in up to 2,366 dizygotic (DZ) twin pairs from the Twins Early Development Study from age 12 to age 21. To account for family clustering, we used mixed-effects modeling, simultaneously estimating within- and between-family effects for target- and cross-trait GPS prediction of the outcomes. There were three main findings: (1) DZ twin GPS differences predicted DZ differences in height, BMI, intelligence, educational achievement, and ADHD symptoms; (2) target and cross-trait analyses indicated that GPS prediction estimates for cognitive traits (intelligence and educational achievement) were on average 60% greater between families than within families, but this was not the case for non-cognitive traits; and (3) much of this within- and between-family difference for cognitive traits disappeared after controlling for family socio-economic status (SES), suggesting that SES is a major source of between-family prediction through rGE mechanisms. These results provide insights into the patterns by which rGE contributes to GPS prediction, while ruling out confounding due to population stratification and assortative mating.

Phenotypic plasticity in plant height shaped by interaction between genetic loci and diurnal temperature range.

Phenotypic plasticity is observed widely in plants and often studied with reaction norms for adult plant or end-of-season traits. Uncovering genetic, environmental and developmental patterns behind the observed phenotypic variation under natural field conditions is needed. Using a sorghum (Sorghum bicolor) genetic population evaluated for plant height in seven natural field conditions, we investigated the major pattern that differentiated these environments. We then examined the physiological relevance of the identified environmental index by investigating the developmental trajectory of the population with multistage height measurements in four additional environments and conducting crop growth modelling. We found that diurnal temperature range (DTR) during the rapid growth period of sorghum development was an effective environmental index. Three genetic loci (Dw1, Dw3 and qHT7.1) were consistently detected for individual environments, reaction-norm parameters across environments and growth-curve parameters through the season. Their genetic effects changed dynamically along the environmental gradient and the developmental stage. A conceptual model with three-dimensional reaction norms was proposed to showcase the interconnecting components: genotype, environment and development. Beyond genomic and environmental analyses, further integration of development and physiology at the whole-plant and molecular levels into complex trait dissection would enhance our understanding of mechanisms underlying phenotypic variation.

von Willebrand disease and von Willebrand factor.

Progress in both basic and translational research into the molecular mechanisms of VWD can be seen in multiple fields. GENETICS OF VWD: In the past several decades, knowledge of the underlying pathogenesis of von Willebrand disease (VWD) has increased tremendously, thanks in no small part to detailed genetic mapping of the von Willebrand Factor (VWF) gene and advances in genetic and bioinformatic technology. However, these advances do not always easily translate into improved management for patients with VWD and low-VWF levels. VWD AND PREGNANCY: For example, the treatment of pregnant women with VWD both pre- and postpartum can be complicated. While knowledge of the VWF genotype at some amino acid positions can aid in knowledge of who may be at increased risk of thrombocytopenia or insufficient increase in VWF levels during pregnancy, in many cases, VWF levels and bleeding severity is highly heterogeneous, making monitoring recommended during pregnancy to optimize treatment strategies. VWF AND COVID-19: New challenges related to the consequences of dysregulation of hemostasis continue to be discovered. The ongoing COVID-19 pandemic has highlighted that VWF has additional biological roles in the regulation of inflammatory disorders and angiogenesis, disruption of which may contribute to COVID-19 induced vasculopathy. Increased endothelial cell activation and Weibel-Palade body exocytosis in severe COVID-19 lead to markedly increased plasma VWF levels. Coupled with impairment of normal ADAMTS13 multimer regulation, these data suggest a role for VWF in the pathogenesis underlying pulmonary microvascular angiopathy in severe COVID-19. CONCLUSION: With the increased affordability and availability of next-generation sequencing techniques, as well as a push towards a multi-omic approach and personalized medicine in human genetics, there is hope that translational research will improve VWD patient outcomes.

Rethinking the Genetic Etiology of Nonsyndromic Tooth Agenesis.

PURPOSE OF REVIEW: Genetic studies in humans and animal models have improved our understanding of the role of numerous genes in the etiology of nonsyndromic tooth agenesis (TA). The purpose of this review is to discuss recently identified genes potentially contributing to TA. RECENT FINDINGS: Despite research progress, understanding the genetic factors underlying nonsyndromic TA has been challenging given the genetic heterogeneity, variable expressivity, and incomplete penetrance of putatively pathogenic variants often observed associated with the condition. Next-generation sequencing technologies have provided a platform for novel gene and variant discoveries and informed paradigm-shifting concepts in the etiology of TA. This review summarizes the current knowledge on genes and pathways related to nonsyndromic TA with a focus on recently identified genes/variants. Evidence suggesting possible multi-locus variation in TA is also presented.

The adaptive benefit of evolved increases in hemoglobin-O2 affinity is contingent on tissue O2 diffusing capacity in high-altitude deer mice.

BACKGROUND: Complex organismal traits are often the result of multiple interacting genes and sub-organismal phenotypes, but how these interactions shape the evolutionary trajectories of adaptive traits is poorly understood. We examined how functional interactions between cardiorespiratory traits contribute to adaptive increases in the capacity for aerobic thermogenesis (maximal O2 consumption, V̇O2max, during acute cold exposure) in high-altitude deer mice (Peromyscus maniculatus). We crossed highland and lowland deer mice to produce F2 inter-population hybrids, which expressed genetically based variation in hemoglobin (Hb) O2 affinity on a mixed genetic background. We then combined physiological experiments and mathematical modeling of the O2 transport pathway to examine the links between cardiorespiratory traits and V̇O2max. RESULTS: Physiological experiments revealed that increases in Hb-O2 affinity of red blood cells improved blood oxygenation in hypoxia but were not associated with an enhancement in V̇O2max. Sensitivity analyses performed using mathematical modeling showed that the influence of Hb-O2 affinity on V̇O2max in hypoxia was contingent on the capacity for O2 diffusion in active tissues. CONCLUSIONS: These results suggest that increases in Hb-O2 affinity would only have adaptive value in hypoxic conditions if concurrent with or preceded by increases in tissue O2 diffusing capacity. In high-altitude deer mice, the adaptive benefit of increasing Hb-O2 affinity is contingent on the capacity to extract O2 from the blood, which helps resolve controversies about the general role of hemoglobin function in hypoxia tolerance.

The Origin of Additive Genetic Variance Driven by Positive Selection.

Fisher's fundamental theorem of natural selection predicts no additive variance of fitness in a natural population. Consistently, studies in a variety of wild populations show virtually no narrow-sense heritability (h2) for traits important to fitness. However, counterexamples are occasionally reported, calling for a deeper understanding on the evolution of additive variance. In this study, we propose adaptive divergence followed by population admixture as a source of the additive genetic variance of evolutionarily important traits. We experimentally tested the hypothesis by examining a panel of ∼1,000 yeast segregants produced by a hybrid of two yeast strains that experienced adaptive divergence. We measured >400 yeast cell morphological traits and found a strong positive correlation between h2 and evolutionary importance. Because adaptive divergence followed by population admixture could happen constantly, particularly in species with wide geographic distribution and strong migratory capacity (e.g., humans), the finding reconciles the observation of abundant additive variances in evolutionarily important traits with Fisher's fundamental theorem of natural selection. Importantly, the revealed role of positive selection in promoting rather than depleting additive variance suggests a simple explanation for why additive genetic variance can be dominant in a population despite the ubiquitous between-gene epistasis observed in functional assays.

Deciphering the Emerging Complexities of Molecular Mechanisms at GWAS Loci.

Genome-wide association studies (GWASs) have identified thousands of loci associated with hundreds of complex diseases and traits, and progress is being made toward elucidating the causal variants and genes underlying these associations. Functional characterization of mechanisms at GWAS loci is a multi-faceted challenge. Challenges include linkage disequilibrium and allelic heterogeneity at each locus, the noncoding nature of most loci, and the time and cost needed for experimentally evaluating the potential mechanistic contributions of genes and variants. As GWAS sample sizes increase, more loci are identified, and the complexities of individual loci emerge. Loci can consist of multiple association signals, each of which can reflect the influence of multiple variants, inseparable by association analyses. Each signal within a locus can influence the same or different target genes. Experimental studies of genes and variants can differ on the basis of cell type, cellular environment, or other context-specific variables. In this review, we describe the complexity of mechanisms at GWAS loci-including multiple signals, multiple variants, and/or multiple genes-and the implications these complexities hold for experimental study design and interpretation of GWAS mechanisms.

Landscape of Conditional eQTL in Dorsolateral Prefrontal Cortex and Co-localization with Schizophrenia GWAS.

Causal genes and variants within genome-wide association study (GWAS) loci can be identified by integrating GWAS statistics with expression quantitative trait loci (eQTL) and determining which variants underlie both GWAS and eQTL signals. Most analyses, however, consider only the marginal eQTL signal, rather than dissect this signal into multiple conditionally independent signals for each gene. Here we show that analyzing conditional eQTL signatures, which could be important under specific cellular or temporal contexts, leads to improved fine mapping of GWAS associations. Using genotypes and gene expression levels from post-mortem human brain samples (n = 467) reported by the CommonMind Consortium (CMC), we find that conditional eQTL are widespread; 63% of genes with primary eQTL also have conditional eQTL. In addition, genomic features associated with conditional eQTL are consistent with context-specific (e.g., tissue-, cell type-, or developmental time point-specific) regulation of gene expression. Integrating the 2014 Psychiatric Genomics Consortium schizophrenia (SCZ) GWAS and CMC primary and conditional eQTL data reveals 40 loci with strong evidence for co-localization (posterior probability > 0.8), including six loci with co-localization of conditional eQTL. Our co-localization analyses support previously reported genes, identify novel genes associated with schizophrenia risk, and provide specific hypotheses for their functional follow-up.

Functional gene networks reveal distinct mechanisms segregating in migraine families.

Migraine is the most common neurological disorder worldwide and it has been shown to have complex polygenic origins with a heritability of estimated 40-70%. Both common and rare genetic variants are believed to underlie the pathophysiology of the prevalent types of migraine, migraine with typical aura and migraine without aura. However, only common variants have been identified so far. Here we identify for the first time a gene module with rare mutations through a systems genetics approach integrating RNA sequencing data from brain and vascular tissues likely to be involved in migraine pathology in combination with whole genome sequencing of 117 migraine families. We found a gene module in the visual cortex, based on single nuclei RNA sequencing data, that had increased rare mutations in the migraine families and replicated this in a second independent cohort of 1930 patients. This module was mainly expressed by interneurons, pyramidal CA1, and pyramidal SS cells, and pathway analysis showed association with hormonal signalling (thyrotropin-releasing hormone receptor and oxytocin receptor signalling pathways), Alzheimer's disease pathway, serotonin receptor pathway and general heterotrimeric G-protein signalling pathways. Our results demonstrate that rare functional gene variants are strongly implicated in the pathophysiology of migraine. Furthermore, we anticipate that the results can be used to explain the critical mechanisms behind migraine and potentially improving the treatment regime for migraine patients.