Interpreting meta-analyses of genome-wide association studies

Genome-Wide Association Study Points to Novel Locus for Gilles de la Tourette Syndrome

BACKGROUND

Tourette syndrome (TS) is a childhood-onset neurodevelopmental disorder of complex genetic architecture and is characterized by multiple motor tics and at least one vocal tic persisting for more than 1 year.

METHODS

We performed a genome-wide meta-analysis integrating a novel TS cohort with previously published data, resulting in a sample size of 6133 individuals with TS and 13,565 ancestry-matched control participants.

RESULTS

We identified a genome-wide significant locus on chromosome 5q15. Integration of expression quantitative trait locus, Hi-C (high-throughput chromosome conformation capture), and genome-wide association study data implicated the NR2F1 gene and associated long noncoding RNAs within the 5q15 locus. Heritability partitioning identified statistically significant enrichment in brain tissue histone marks, while polygenic risk scoring of brain volume data identified statistically significant associations with right and left thalamus volumes and right putamen volume.

CONCLUSIONS

Our work presents novel insights into the neurobiology of TS, thereby opening up new directions for future studies.

Functional Genomics in Psoriasis

Psoriasis is an autoimmune cutaneous condition that significantly impacts quality of life and represents a burden on society due to its prevalence. Genome-wide association studies (GWASs) have pinpointed several psoriasis-related risk loci, underlining the disease's complexity. Functional genomics is paramount to unveiling the role of such loci in psoriasis and disentangling its complex nature. In this review, we aim to elucidate the main findings in this field and integrate our discussion with gold-standard techniques in molecular biology-i.e., Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-and high-throughput technologies. These tools are vital to understanding how disease risk loci affect gene expression in psoriasis, which is crucial in identifying new targets for personalized treatments in advanced precision medicine.

Shared molecular mechanisms and transdiagnostic potential of neurodevelopmental disorders and immune disorders

The co-occurrence and familial clustering of neurodevelopmental disorders and immune disorders suggest shared genetic risk factors. Based on genome-wide association summary statistics from five neurodevelopmental disorders and four immune disorders, we conducted genome-wide, local genetic correlation and polygenic overlap analysis. We further performed a cross-trait GWAS meta-analysis. Pleotropic loci shared between the two categories of diseases were mapped to candidate genes using multiple algorithms and approaches. Significant genetic correlations were observed between neurodevelopmental disorders and immune disorders, including both positive and negative correlations. Neurodevelopmental disorders exhibited higher polygenicity compared to immune disorders. Around 50%-90% of genetic variants of the immune disorders were shared with neurodevelopmental disorders. The cross-trait meta-analysis revealed 154 genome-wide significant loci, including 8 novel pleiotropic loci. Significant associations were observed for 30 loci with both types of diseases. Pathway analysis on the candidate genes at these loci revealed common pathways shared by the two types of diseases, including neural signaling, inflammatory response, and PI3K-Akt signaling pathway. In addition, 26 of the 30 lead SNPs were associated with blood cell traits. Neurodevelopmental disorders exhibit complex polygenic architecture, with a subset of individuals being at a heightened genetic risk for both neurodevelopmental and immune disorders. The identification of pleiotropic loci has important implications for exploring opportunities for drug repurposing, enabling more accurate patient stratification, and advancing genomics-informed precision in the medical field of neurodevelopmental disorders.

Genome-wide association study identifies high-impact susceptibility loci for HCC in North America

BACKGROUND AND AIMS

Despite the substantial impact of environmental factors, individuals with a family history of liver cancer have an increased risk for HCC. However, genetic factors have not been studied systematically by genome-wide approaches in large numbers of individuals from European descent populations (EDP).

APPROACH AND RESULTS

We conducted a 2-stage genome-wide association study (GWAS) on HCC not affected by HBV infections. A total of 1872 HCC cases and 2907 controls were included in the discovery stage, and 1200 HCC cases and 1832 controls in the validation. We analyzed the discovery and validation samples separately and then conducted a meta-analysis. All analyses were conducted in the presence and absence of HCV. The liability-scale heritability was 24.4% for overall HCC. Five regions with significant ORs (95% CI) were identified for nonviral HCC: 3p22.1, MOBP , rs9842969, (0.51, [0.40-0.65]); 5p15.33, TERT , rs2242652, (0.70, (0.62-0.79]); 19q13.11, TM6SF2 , rs58542926, (1.49, [1.29-1.72]); 19p13.11 MAU2 , rs58489806, (1.53, (1.33-1.75]); and 22q13.31, PNPLA3 , rs738409, (1.66, [1.51-1.83]). One region was identified for HCV-induced HCC: 6p21.31, human leukocyte antigen DQ beta 1, rs9275224, (0.79, [0.74-0.84]). A combination of homozygous variants of PNPLA3 and TERT showing a 6.5-fold higher risk for nonviral-related HCC compared to individuals lacking these genotypes. This observation suggests that gene-gene interactions may identify individuals at elevated risk for developing HCC.

CONCLUSIONS

Our GWAS highlights novel genetic susceptibility of nonviral HCC among European descent populations from North America with substantial heritability. Selected genetic influences were observed for HCV-positive HCC. Our findings indicate the importance of genetic susceptibility to HCC development.

Complementation testing identifies genes mediating effects at quantitative trait loci underlying fear-related behavior

Knowing the genes involved in quantitative traits provides an entry point to understanding the biological bases of behavior, but there are very few examples where the pathway from genetic locus to behavioral change is known. To explore the role of specific genes in fear behavior, we mapped three fear-related traits, tested fourteen genes at six quantitative trait loci (QTLs) by quantitative complementation, and identified six genes. Four genes, Lamp, Ptprd, Nptx2, and Sh3gl, have known roles in synapse function; the fifth, Psip1, was not previously implicated in behavior; and the sixth is a long non-coding RNA, 4933413L06Rik, of unknown function. Variation in transcriptome and epigenetic modalities occurred preferentially in excitatory neurons, suggesting that genetic variation is more permissible in excitatory than inhibitory neuronal circuits. Our results relieve a bottleneck in using genetic mapping of QTLs to uncover biology underlying behavior and prompt a reconsideration of expected relationships between genetic and functional variation.

Integrative cross-omics and cross-context analysis elucidates molecular links underlying genetic effects on complex traits

Genetic effects on functionally related 'omic' traits often co-occur in relevant cellular contexts, such as tissues. Motivated by the multi-tissue methylation quantitative trait loci (mQTLs) and expression QTLs (eQTLs) analysis, we propose X-ING (Cross-INtegrative Genomics) for cross-omics and cross-context integrative analysis. X-ING takes as input multiple matrices of association statistics, each obtained from different omics data types across multiple cellular contexts. It models the latent binary association status of each statistic, captures the major association patterns among omics data types and contexts, and outputs the posterior mean and probability for each input statistic. X-ING enables the integration of effects from different omics data with varying effect distributions. In the multi-tissue cis-association analysis, X-ING shows improved detection and replication of mQTLs by integrating eQTL maps. In the trans-association analysis, X-ING reveals an enrichment of trans-associations in many disease/trait-relevant tissues.

GenomeMUSter mouse genetic variation service enables multitrait, multipopulation data integration and analysis

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Complementation testing identifies causal genes at quantitative trait loci underlying fear related behavior

Knowing the genes involved in quantitative traits provides a critical entry point to understanding the biological bases of behavior, but there are very few examples where the pathway from genetic locus to behavioral change is known. Here we address a key step towards that goal by deploying a test that directly queries whether a gene mediates the effect of a quantitative trait locus (QTL). To explore the role of specific genes in fear behavior, we mapped three fear-related traits, tested fourteen genes at six QTLs, and identified six genes. Four genes, Lsamp, Ptprd, Nptx2 and Sh3gl, have known roles in synapse function; the fifth gene, Psip1, is a transcriptional co-activator not previously implicated in behavior; the sixth is a long non-coding RNA 4933413L06Rik with no known function. Single nucleus transcriptomic and epigenetic analyses implicated excitatory neurons as likely mediating the genetic effects. Surprisingly, variation in transcriptome and epigenetic modalities between inbred strains occurred preferentially in excitatory neurons, suggesting that genetic variation is more permissible in excitatory than inhibitory neuronal circuits. Our results open a bottleneck in using genetic mapping of QTLs to find novel biology underlying behavior and prompt a reconsideration of expected relationships between genetic and functional variation.

A compendium of genetic regulatory effects across pig tissues

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

Genome-wide association study meta-analysis of blood pressure traits and hypertension in sub-Saharan African populations: an AWI-Gen study

Most hypertension-related genome-wide association studies (GWASs) focus on non-African populations, despite hypertension (a major risk factor for cardiovascular disease) being highly prevalent in Africa. The AWI-Gen study GWAS meta-analysis for blood pressure (BP)-related traits (systolic and diastolic BP, pulse pressure, mean-arterial pressure and hypertension) from three sub-Saharan African geographic regions (N = 10,775), identifies two novel genome-wide significant signals (p < 5E-08): systolic BP near P2RY1 (rs77846204; intergenic variant, p = 4.95E-08) and pulse pressure near LINC01256 (rs80141533; intergenic variant, p = 1.76E-08). No genome-wide signals are detected for the AWI-Gen GWAS meta-analysis with previous African-ancestry GWASs (UK Biobank (African), Uganda Genome Resource). Suggestive signals (p < 5E-06) are observed for all traits, with 29 SNPs associating with more than one trait and several replicating known associations. Polygenic risk scores (PRSs) developed from studies on different ancestries have limited transferability, with multi-ancestry PRS providing better prediction. This study provides insights into the genetics of BP variation in African populations.

Mouse phenome database: curated data repository with interactive multi-population and multi-trait analyses

The Mouse Phenome Database continues to serve as a curated repository and analysis suite for measured attributes of members of diverse mouse populations. The repository includes annotation to community standard ontologies and guidelines, a database of allelic states for 657 mouse strains, a collection of protocols, and analysis tools for flexible, interactive, user directed analyses that increasingly integrates data across traits and populations. The database has grown from its initial focus on a standard set of inbred strains to include heterogeneous mouse populations such as the Diversity Outbred and mapping crosses and well as Collaborative Cross, Hybrid Mouse Diversity Panel, and recombinant inbred strains. Most recently the system has expanded to include data from the International Mouse Phenotyping Consortium. Collectively these data are accessible by API and provided with an interactive tool suite that enables users' persistent selection, storage, and operation on collections of measures. The tool suite allows basic analyses, advanced functions with dynamic visualization including multi-population meta-analysis, multivariate outlier detection, trait pattern matching, correlation analyses and other functions. The data resources and analysis suite provide users a flexible environment in which to explore the basis of phenotypic variation in health and disease across the lifespan.

Polygenic risk scores identify heterogeneity in asthma and chronic obstructive pulmonary disease

BACKGROUND

Asthma and chronic obstructive pulmonary disease (COPD) have distinct and overlapping genetic and clinical features.

OBJECTIVE

We sought to test the hypothesis that polygenic risk scores (PRSs) for asthma (PRSAsthma) and spirometry (FEV1 and FEV1/forced vital capacity; PRSspiro) would demonstrate differential associations with asthma, COPD, and asthma-COPD overlap (ACO).

METHODS

We developed and tested 2 asthma PRSs and applied the higher performing PRSAsthma and a previously published PRSspiro to research (Genetic Epidemiology of COPD study and Childhood Asthma Management Program, with spirometry) and electronic health record-based (Mass General Brigham Biobank and Genetic Epidemiology Research on Adult Health and Aging [GERA]) studies. We assessed the association of PRSs with COPD and asthma using modified random-effects and binary-effects meta-analyses, and ACO and asthma exacerbations in specific cohorts. Models were adjusted for confounders and genetic ancestry.

RESULTS

In meta-analyses of 102,477 participants, the PRSAsthma (odds ratio [OR] per SD, 1.16 [95% CI, 1.14-1.19]) and PRSspiro (OR per SD, 1.19 [95% CI, 1.17-1.22]) both predicted asthma, whereas the PRSspiro predicted COPD (OR per SD, 1.25 [95% CI, 1.21-1.30]). However, results differed by cohort. The PRSspiro was not associated with COPD in GERA and Mass General Brigham Biobank. In the Genetic Epidemiology of COPD study, the PRSAsthma (OR per SD: Whites, 1.3; African Americans, 1.2) and PRSspiro (OR per SD: Whites, 2.2; African Americans, 1.6) were both associated with ACO. In GERA, the PRSAsthma was associated with asthma exacerbations (OR, 1.18) in Whites; the PRSspiro was associated with asthma exacerbations in White, LatinX, and East Asian participants.

CONCLUSIONS

PRSs for asthma and spirometry are both associated with ACO and asthma exacerbations. Genetic prediction performance differs in research versus electronic health record-based cohorts.

Exploring the cross-cancer effect of smoking and its fingerprints in blood DNA methylation on multiple cancers: A Mendelian randomization study

Aberrant smoking-related DNA methylation has been widely investigated as a carcinogenesis mechanism, but whether the cross-cancer epigenetic pathways exist remains unclear. We conducted two-sample Mendelian randomization (MR) analyses respectively on smoking behaviors (age of smoking initiation, smoking initiation, smoking cessation, and lifetime smoking index [LSI]) and smoking-related DNA methylation to investigate their effect on 15 site-specific cancers, based on a genome-wide association study (GWAS) of 1.2 million European individuals and an epigenome-WAS (EWAS) of 5907 blood samples of Europeans for smoking and 15 GWASs of European ancestry for multiple site-specific cancers. Significantly identified CpG sites were further used for colocalization analysis, and those with cross-cancer effect were validated by overlapping with tissue-specific eQTLs. In the genomic MR, smoking measurements of smoking initiation, smoking cessation and LSI were suggested to be casually associated with risk of seven types of site-specific cancers, among which cancers at lung, cervix and colorectum were provided with strong evidence. In the epigenetic MR, methylation at 75 CpG sites were reported to be significantly associated with increased risks of multiple cancers. Eight out of 75 CpG sites were observed with cross-cancer effect, among which cg06639488 (EFNA1), cg12101586 (CYP1A1) and cg14142171 (HLA-L) were validated by eQTLs at specific cancer sites, and cg07932199 (ATXN2) had strong evidence to be associated with cancers of lung (coefficient, 0.65, 95% confidence interval [CI], 0.31-1.00), colorectum (0.90 [0.61, 1.18]), breast (0.31 [0.20, 0.43]) and endometrium (0.98 [0.68, 1.27]). These findings highlight the potential practices targeting DNA methylation-involved cross-cancer pathways.

Statistical assessment of biomarker replicability using MAJAR method

In the era of precision medicine, many biomarkers have been discovered to be associated with drug efficacy and safety responses, which can be used for patient stratification and drug response prediction. Due to the small sample size and limited power of randomized clinical studies, meta-analysis is usually conducted to aggregate all available studies to maximize the power for identifying prognostic and predictive biomarkers. However, it is often challenging to find an independent study to replicate the discoveries from the meta-analysis (e.g. meta-analysis of pharmacogenomics genome-wide association studies (PGx GWAS)), which seriously limits the potential impacts of the discovered biomarkers. To overcome this challenge, we develop a novel statistical framework, MAJAR (meta-analysis of joint effect associations for biomarker replicability assessment), to jointly test prognostic and predictive effects and assess the replicability of identified biomarkers by implementing an enhanced expectation-maximization algorithm and calculating their posterior-probability-of-replicabilities and Bayesian false discovery rates (Fdr). Extensive simulation studies were conducted to compare the performance of MAJAR and existing methods in terms of Fdr, power, and computational efficiency. The simulation results showed improved statistical power with well-controlled Fdr of MAJAR over existing methods and robustness to outliers under different data generation processes. We further demonstrated the advantages of MAJAR over existing methods by applying MAJAR to the PGx GWAS summary statistics data from a large cardiovascular randomized clinical trial. Compared to testing main effects only, MAJAR identified 12 novel variants associated with the treatment-related low-density lipoprotein cholesterol reduction from baseline.

PheWAS and cross-disorder analysis reveal genetic architecture, pleiotropic loci and phenotypic correlations across 11 autoimmune disorders

Introduction

Autoimmune disorders (ADs) are a group of about 80 disorders that occur when self-attacking autoantibodies are produced due to failure in the self-tolerance mechanisms. ADs are polygenic disorders and associations with genes both in the human leukocyte antigen (HLA) region and outside of it have been described. Previous studies have shown that they are highly comorbid with shared genetic risk factors, while epidemiological studies revealed associations between various lifestyle and health-related phenotypes and ADs.

Methods

Here, for the first time, we performed a comparative polygenic risk score (PRS) - Phenome Wide Association Study (PheWAS) for 11 different ADs (Juvenile Idiopathic Arthritis, Primary Sclerosing Cholangitis, Celiac Disease, Multiple Sclerosis, Rheumatoid Arthritis, Psoriasis, Myasthenia Gravis, Type 1 Diabetes, Systemic Lupus Erythematosus, Vitiligo Late Onset, Vitiligo Early Onset) and 3,254 phenotypes available in the UK Biobank that include a wide range of socio-demographic, lifestyle and health-related outcomes. Additionally, we investigated the genetic relationships of the studied ADs, calculating their genetic correlation and conducting cross-disorder GWAS meta-analyses for the observed AD clusters.

Results

In total, we identified 508 phenotypes significantly associated with at least one AD PRS. 272 phenotypes were significantly associated after excluding variants in the HLA region from the PRS estimation. Through genetic correlation and genetic factor analyses, we identified four genetic factors that run across studied ADs. Cross-trait meta-analyses within each factor revealed pleiotropic genome-wide significant loci.

Discussion

Overall, our study confirms the association of different factors with genetic susceptibility for ADs and reveals novel observations that need to be further explored.

An RNA-informed dosage sensitivity map reflects the intrinsic functional nature of genes

Understanding dosage sensitivity or why Mendelian diseases have dominant vs. recessive modes of inheritance is crucial for uncovering the etiology of human disease. Previous knowledge of dosage sensitivity is mainly based on observations of rare loss-of-function mutations or copy number changes, which are underpowered due to ultra rareness of such variants. Thus, the functional underpinnings of dosage constraint remain elusive. In this study, we aim to systematically quantify dosage perturbations from cis-regulatory variants in the general population to yield a tissue-specific dosage constraint map of genes and further explore their underlying functional logic. We reveal an inherent divergence of dosage constraints in genes by functional categories with signaling genes (transcription factors, protein kinases, ion channels, and cellular machinery) being dosage sensitive, while effector genes (transporters, metabolic enzymes, cytokines, and receptors) are generally dosage resilient. Instead of being a metric of functional dispensability, we show that dosage constraint reflects underlying homeostatic constraints arising from negative feedback. Finally, we employ machine learning to integrate DNA and RNA metrics to generate a comprehensive, tissue-specific map of dosage sensitivity (MoDs) for autosomal genes.

Candidate pathway analysis of surfactant proteins identifies CTSH and SFTA2 that influences lung cancer risk

Pulmonary surfactant is a lipoprotein synthesized and secreted by alveolar type II cells in lung. We evaluated the associations between 200,139 single nucleotide polymorphisms (SNPs) of 40 surfactant-related genes and lung cancer risk using genotyped data from two independent lung cancer genome-wide association studies. Discovery data included 18,082 cases and 13,780 controls of European ancestry. Replication data included 1,914 cases and 3,065 controls of European descent. Using multivariate logistic regression, we found novel SNPs in surfactant-related genes CTSH [rs34577742 C > T, odds ratio (OR) = 0.90, 95% confidence interval (CI) = 0.89-0.93, P = 7.64 × 10-9] and SFTA2 (rs3095153 G > A, OR = 1.16, 95% CI = 1.10-1.21, P = 1.27 × 10-9) associated with overall lung cancer in the discovery data and validated in an independent replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.80-0.96, P = 5.76 × 10-3) and SFTA2 (rs3095153 G > A, OR = 1.14, 95% CI = 1.01-1.28, P = 3.25 × 10-2). Among ever smokers, we found SNPs in CTSH (rs34577742 C > T, OR = 0.89, 95% CI = 0.85-0.92, P = 1.94 × 10-7) and SFTA2 (rs3095152 G > A, OR = 1.20, 95% CI = 1.14-1.27, P = 4.25 × 10-11) associated with overall lung cancer in the discovery data and validated in the replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.79-0.97, P = 1.64 × 10-2) and SFTA2 (rs3095152 G > A, OR = 1.15, 95% CI = 1.01-1.30, P = 3.81 × 10-2). Subsequent transcriptome-wide association study using expression weights from a lung expression quantitative trait loci study revealed genes most strongly associated with lung cancer are CTSH (PTWAS = 2.44 × 10-4) and SFTA2 (PTWAS = 2.32 × 10-6).

GWAS meta-analysis of over 29,000 people with epilepsy identifies 26 risk loci and subtype-specific genetic architecture

Epilepsy is a highly heritable disorder affecting over 50 million people worldwide, of which about one-third are resistant to current treatments. Here we report a multi-ancestry genome-wide association study including 29,944 cases, stratified into three broad categories and seven subtypes of epilepsy, and 52,538 controls. We identify 26 genome-wide significant loci, 19 of which are specific to genetic generalized epilepsy (GGE). We implicate 29 likely causal genes underlying these 26 loci. SNP-based heritability analyses show that common variants explain between 39.6% and 90% of genetic risk for GGE and its subtypes. Subtype analysis revealed markedly different genetic architectures between focal and generalized epilepsies. Gene-set analyses of GGE signals implicate synaptic processes in both excitatory and inhibitory neurons in the brain. Prioritized candidate genes overlap with monogenic epilepsy genes and with targets of current antiseizure medications. Finally, we leverage our results to identify alternate drugs with predicted efficacy if repurposed for epilepsy treatment.

Sex differences in the pleiotropy of hearing difficulty with imaging-derived phenotypes: a brain-wide investigation

Background

Hearing difficulty (HD) is one of the major health burdens in older adults. While aging-related changes in the peripheral auditory system play an important role, genetic variation associated with brain structure and function could also be involved in HD predisposition.

Methods

We analyzed a large-scale HD genome-wide association study (GWAS; N total = 501,825, 56% females) and GWAS data related to 3,935 brain imaging-derived phenotypes (IDPs) assessed in up to 33,224 individuals (52% females) using multiple magnetic resonance imaging (MRI) modalities. To investigate HD pleiotropy with brain structure and function, we conducted genetic correlation, latent causal variable (LCV), Mendelian randomization (MR), and multivariable generalized linear regression analyses. Additionally, we performed local genetic correlation and multi-trait colocalization analyses to identify genomic regions and loci implicated in the pleiotropic mechanisms shared between HD and brain IDPs.

Results

We observed a widespread genetic correlation of HD with 120 IDPs in females, 89 IDPs in males, and 171 IDPs in the sex-combined analysis. The LCV analyses showed that some of these genetic correlations could be due to cause-effect relationships. For seven correlations, the causal effects were also confirmed by the MR approach: vessel volume→HD in the sex-combined analysis; hippocampus volume→HD, cerebellum grey matter volume→HD, primary visual cortex volume→HD, and HD→rfMRI-ICA100 node 46 in females; global mean thickness→HD and HD→mean orientation dispersion index in superior corona radiata in males. The local genetic correlation analyses identified 13 pleiotropic regions between HD and these seven IDPs. We also observed a colocalization signal for the rs13026575 variant between HD, primary visual cortex volume, and SPTBN1 transcriptomic regulation in females.

Conclusion

Brain structure and function may have a role in the sex differences in HD predisposition via possible cause-effect relationships and shared regulatory mechanisms.

Cross-trait analyses identify shared genetics between migraine, headache, and glycemic traits, and a causal relationship with fasting proinsulin

The co-occurrence of migraine and glycemic traits has long been reported in observational epidemiological studies, but it has remained unknown how they are linked genetically. We used large-scale GWAS summary statistics on migraine, headache, and nine glycemic traits in European populations to perform cross-trait analyses to estimate genetic correlation, identify shared genomic regions, loci, genes, and pathways, and test for causal relationships. Out of the nine glycemic traits, significant genetic correlation was observed for fasting insulin (FI) and glycated haemoglobin (HbA1c) with both migraine and headache, while 2-h glucose was genetically correlated only with migraine. Among 1703 linkage disequilibrium (LD) independent regions of the genome, we found pleiotropic regions between migraine and FI, fasting glucose (FG), and HbA1c, and pleiotropic regions between headache and glucose, FI, HbA1c, and fasting proinsulin. Cross-trait GWAS meta-analysis with glycemic traits, identified six novel genome-wide significant lead SNPs with migraine, and six novel lead SNPs with headache (Pmeta < 5.0 × 10-8 and Psingle-trait < 1 × 10-4), all of which were LD-independent. Genes with a nominal gene-based association (Pgene ≤ 0.05) were significantly enriched (overlapping) across the migraine, headache, and glycemic traits. Mendelian randomisation analyses produced intriguing, but inconsistent, evidence for a causal relationship between migraine and headache with multiple glycemic traits; and consistent evidence suggesting increased fasting proinsulin levels may causally decrease the risk of headache. Our findings indicate that migraine, headache, and glycemic traits share a common genetic etiology and provide genetic insights into the molecular mechanisms contributing to their comorbid relationship.

The Impact of Stability Considerations on Genetic Fine-Mapping

Fine-mapping methods, which aim to identify genetic variants responsible for complex traits following genetic association studies, typically assume that sufficient adjustments for confounding within the association study cohort have been made, e.g., through regressing out the top principal components (i.e., residualization). Despite its widespread use, however, residualization may not completely remove all sources of confounding. Here, we propose a complementary stability-guided approach that does not rely on residualization, which identifies consistently fine-mapped variants across different genetic backgrounds or environments. We demonstrate the utility of this approach by applying it to fine-map eQTLs in the GEUVADIS data. Using 378 different functional annotations of the human genome, including recent deep learning-based annotations (e.g., Enformer), we compare enrichments of these annotations among variants for which the stability and traditional residualization-based fine-mapping approaches agree against those for which they disagree, and find that the stability approach enhances the power of traditional fine-mapping methods in identifying variants with functional impact. Finally, in cases where the two approaches report distinct variants, our approach identifies variants comparably enriched for functional annotations. Our findings suggest that the stability principle, as a conceptually simple device, complements existing approaches to fine-mapping, reinforcing recent advocacy of evaluating cross-population and cross-environment portability of biological findings. To support visualization and interpretation of our results, we provide a Shiny app, available at: https://alan-aw.shinyapps.io/stability_v0/.

Genome-wide Association Study Meta-analysis of Blood Pressure Traits and Hypertension in Sub-Saharan African Populations: An AWI-Gen Study

Most hypertension-related genome-wide association studies (GWAS) focus on non-African populations, despite hypertension (a major risk factor for cardiovascular disease) being highly prevalent in Africa. The AWI-Gen study GWAS meta-analysis for blood pressure-related traits (systolic and diastolic blood pressure, pulse pressure, mean-arterial pressure and hypertension) from three sub-Saharan African geographic regions (N=10,775), identified two genome-wide significant signals (p<5E-08): systolic blood pressure near P2RY1 (rs77846204; intergenic variant, p=4.25E-08) and pulse pressure near Linc01256 (rs80141533; intergenic variant, p=4.25E-08). No genome-wide signals were detected for the AWI-Gen GWAS meta-analysis with previous African-ancestry GWASs (UK Biobank (African), Uganda Genome Resource). Suggestive signals (p<5E-06) were observed for all traits, with 29 displaying pleiotropic effects and several replicating known associations. Polygenic risk scores developed from studies on different ancestries had limited transferability, with multi-ancestry models providing better prediction. This study provides insights into the genetics and physiology of blood pressure variation in African populations.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes

BACKGROUND

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes.

METHODS

We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome.

RESULTS

We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs.

CONCLUSIONS

Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

From SNP to pathway-based GWAS meta-analysis: do current meta-analysis approaches resolve power and replication in genetic association studies?

Genome-wide association studies (GWAS) have benefited greatly from enhanced high-throughput technology in recent decades. GWAS meta-analysis has become increasingly popular to highlight the genetic architecture of complex traits, informing about the replicability and variability of effect estimations across human ancestries. A wealth of GWAS meta-analysis methodologies have been developed depending on the input data and the outcome information of interest. We present a survey of current approaches from SNP to pathway-based meta-analysis by acknowledging the range of resources and methodologies in the field, and we provide a comprehensive review of different categories of Genome-Wide Meta-analysis methods employed. These methods highlight different levels at which GWAS meta-analysis may be done, including Single Nucleotide Polymorphisms, Genes and Pathways, for which we describe their framework outline. We also discuss the strengths and pitfalls of each approach and make suggestions regarding each of them.

A genetic correlation and bivariate genome-wide association study of grip strength and depression

Grip strength is an important biomarker reflecting muscle strength, and depression is a psychiatric disorder all over the world. Several studies found a significant inverse association between grip strength and depression, and there is also evidence for common physiological mechanisms between them. We used twin data from Qingdao, China to calculate genetic correlations, and we performed a bivariate GWAS to explore potential SNPs, genes, and pathways in common between grip strength and depression. 139 pairs of Dizygotic twins were used for bivariate GWAS. VEAGSE2 and PASCAL software were used for gene-based analysis and pathway enrichment analysis, respectively. And the resulting SNPs were subjected to eQTL analysis and pleiotropy analysis. The genetic correlation coefficient between grip strength and depression was -0.41 (-0.96, -0.15). In SNP-based analysis, 7 SNPs exceeded the genome-wide significance level (P<5×10-8) and a total of 336 SNPs reached the level of suggestive significance (P<1×10-5). Gene-based analysis and pathway-based analysis identified genes and pathways related to muscle strength and the nervous system. The results of eQTL analysis were mainly enriched in tissues such as the brain, thyroid, and skeletal muscle. Pleiotropy analysis shows that 9 of the 15 top SNPs were associated with both grip strength and depression. In conclusion, this bivariate GWAS identified potentially common pleiotropic SNPs, genes, and pathways in grip strength and depression.

Clinical and genetic analysis of familial neuromyelitis optica spectrum disorder in Chinese: associated with ubiquitin-specific peptidase USP18 gene variants

BACKGROUND

Familial clustering of neuromyelitis optica spectrum disorder (NMOSD) was present in Chinese. This study was to investigate the clinical characteristics and genetic background of familial NMOSD.

METHODS

Through questionnaires in four medical centres in 2016-2020, we identified 10 families with NMOSD aggregation. The statistical differences of clinical characteristics between familial and sporadic NMOSD (22 cases and 459 cases) were summarised. The whole-exome sequencing (WES) for seven families (13 cases and 13 controls) was analysed, compared with our previous WES data for sporadic NMOSD (228 cases and 1 400 controls). The family-based and population-based association and linkage analysis were conducted to identify the pathogenetic genes, the variant impacts were predicted.

RESULTS

The familial occurrence was 0.87% in Chinese. Familial patients had higher expanded disability status scale score than sporadic patients (p=0.03). The single-nucleotide polymorphism (SNP) rs2252257 in the promoter and enhancer of ubiquitin-specific peptidase USP18 was linked to familial NMOSD (p=7.8E-05, logarithm of the odds (LOD)=3.1), SNPs rs361553, rs2252257 and rs5746523 were related to sporadic NMOSD (p=1.29E-10, 3.45E-07 and 2.01E-09, respectively). Patients with the SNP rs361553 T/T genotype had higher recurrence rate than C/T or C/C genotype (1.22±0.85 vs 0.69±0.57 and 0.81±0.65, p=0.003 and 0.001, respectively). SNPs rs361553 and rs2252257 altered USP18 expression in brain and nerve tissues.

CONCLUSION

Most clinical characteristics of familial NMOSD were indistinguishable from sporadic NMOSD except for the worst episodes severity. USP18 with impaired intronic regulatory function contributed to the pathogenesis of NMOSD.

Genetic heterogeneity: Challenges, impacts, and methods through an associative lens

Genetic heterogeneity describes the occurrence of the same or similar phenotypes through different genetic mechanisms in different individuals. Robustly characterizing and accounting for genetic heterogeneity is crucial to pursuing the goals of precision medicine, for discovering novel disease biomarkers, and for identifying targets for treatments. Failure to account for genetic heterogeneity may lead to missed associations and incorrect inferences. Thus, it is critical to review the impact of genetic heterogeneity on the design and analysis of population level genetic studies, aspects that are often overlooked in the literature. In this review, we first contextualize our approach to genetic heterogeneity by proposing a high-level categorization of heterogeneity into "feature," "outcome," and "associative" heterogeneity, drawing on perspectives from epidemiology and machine learning to illustrate distinctions between them. We highlight the unique nature of genetic heterogeneity as a heterogeneous pattern of association that warrants specific methodological considerations. We then focus on the challenges that preclude effective detection and characterization of genetic heterogeneity across a variety of epidemiological contexts. Finally, we discuss systems heterogeneity as an integrated approach to using genetic and other high-dimensional multi-omic data in complex disease research.

H3AGWAS: a portable workflow for genome wide association studies

Background

Genome-wide association studies (GWAS) are a powerful method to detect associations between variants and phenotypes. A GWAS requires several complex computations with large data sets, and many steps may need to be repeated with varying parameters. Manual running of these analyses can be tedious, error-prone and hard to reproduce.

Results

The H3AGWAS workflow from the Pan-African Bioinformatics Network for H3Africa is a powerful, scalable and portable workflow implementing pre-association analysis, implementation of various association testing methods and post-association analysis of results.

Conclusions

The workflow is scalable—laptop to cluster to cloud (e.g., SLURM, AWS Batch, Azure). All required software is containerised and can run under Docker or Singularity.

Mouse Phenome Database: towards a more FAIR-compliant and TRUST-worthy data repository and tool suite for phenotypes and genotypes

The Mouse Phenome Database (MPD; https://phenome.jax.org; RRID:SCR_003212), supported by the US National Institutes of Health, is a Biomedical Data Repository listed in the Trans-NIH Biomedical Informatics Coordinating Committee registry. As an increasingly FAIR-compliant and TRUST-worthy data repository, MPD accepts phenotype and genotype data from mouse experiments and curates, organizes, integrates, archives, and distributes those data using community standards. Data are accompanied by rich metadata, including widely used ontologies and detailed protocols. Data are from all over the world and represent genetic, behavioral, morphological, and physiological disease-related characteristics in mice at baseline or those exposed to drugs or other treatments. MPD houses data from over 6000 strains and populations, representing many reproducible strain types and heterogenous populations such as the Diversity Outbred where each mouse is unique but can be genotyped throughout the genome. A suite of analysis tools is available to aggregate, visualize, and analyze these data within and across studies and populations in an increasingly traceable and reproducible manner. We have refined existing resources and developed new tools to continue to provide users with access to consistent, high-quality data that has translational relevance in a modernized infrastructure that enables interaction with a suite of bioinformatics analytic and data services.

Leveraging pleiotropy for joint analysis of genome-wide association studies with per trait interpretations

We introduce pleiotropic association test (PAT) for joint analysis of multiple traits using genome-wide association study (GWAS) summary statistics. The method utilizes the decomposition of phenotypic covariation into genetic and environmental components to create a likelihood ratio test statistic for each genetic variant. Though PAT does not directly interpret which trait(s) drive the association, a per trait interpretation of the omnibus p-value is provided through an extension to the meta-analysis framework, m-values. In simulations, we show PAT controls the false positive rate, increases statistical power, and is robust to model misspecifications of genetic effect. Additionally, simulations comparing PAT to three multi-trait methods, HIPO, MTAG, and ASSET, show PAT identified 15.3% more omnibus associations over the next best method. When these associations were interpreted on a per trait level using m-values, PAT had 37.5% more true per trait interpretations with a 0.92% false positive assignment rate. When analyzing four traits from the UK Biobank, PAT discovered 22,095 novel variants. Through the m-values interpretation framework, the number of per trait associations for two traits were almost tripled and were nearly doubled for another trait relative to the original single trait GWAS.

Genetic Overlap Analysis Identifies a Shared Etiology between Migraine and Headache with Type 2 Diabetes

Migraine and headache frequently co-occur with type 2 diabetes (T2D), suggesting a shared aetiology between the two conditions. We used genome-wide association study (GWAS) data to investigate the genetic overlap and causal relationship between migraine and headache with T2D. Using linkage disequilibrium score regression (LDSC), we found a significant genetic correlation between migraine and T2D (rg = 0.06, p = 1.37 × 10−5) and between headache and T2D (rg = 0.07, p = 3.0 × 10−4). Using pairwise GWAS (GWAS-PW) analysis, we identified 11 pleiotropic regions between migraine and T2D and 5 pleiotropic regions between headache and T2D. Cross-trait SNP meta-analysis identified 23 novel SNP loci (Pmeta < 5 × 10−8) associated with migraine and T2D, and three novel SNP loci associated with headache and T2D. Cross-trait gene-based overlap analysis identified 33 genes significantly associated (Pgene-based < 3.85 × 10−6) with migraine and T2D, and 11 genes associated with headache and T2D, with 7 genes (EHMT2, SLC44A4, PLEKHA1, CFDP1, TMEM170A, CHST6, and BCAR1) common between them. There was also a significant overlap of genes nominally associated (Pgene-based < 0.05) with both migraine and T2D (Pbinomial-test = 2.83 × 10−46) and headache and T2D (Pbinomial-test = 4.08 × 10−29). Mendelian randomisation (MR) analyses did not provide consistent evidence for a causal relationship between migraine and T2D. However, we found headache was causally associated (inverse-variance weighted, ORIVW = 0.90, Pivw = 7 × 10−3) with T2D. Our findings robustly confirm the comorbidity of migraine and headache with T2D, with shared genetically controlled biological mechanisms contributing to their co-occurrence, and evidence for a causal relationship between headache and T2D.

Multinational Genome-Wide Association Study and Functional Genomics Analysis Implicates Decreased SIRT3 Expression Underlying Intracranial Aneurysm Risk

BACKGROUND

The genetic mechanisms regulating intracranial aneurysm (IA) formation and rupture are largely unknown. To identify germline-genetic risk factors for IA, we perform a multinational genome-wide association study (GWAS) of individuals from the United Kingdom, Finland, and Japan.

OBJECTIVE

To identify a shared, multinational genetic basis of IA.

METHODS

Using GWAS summary statistics from UK Biobank, FinnGen, and Biobank Japan, we perform a meta-analysis of IA, containing ruptured and unruptured IA cases. Logistic regression was used to identify IA-associated single-nucleotide polymorphisms. Effect size was calculated using the coefficient r , estimating the contribution of the single-nucleotide polymorphism to the genetic variance of the trait. Genome-wide significance was set at 5.0 × 10 -8 . Expression quantitative trait loci mapping and functional genomics approaches were used to infer mechanistic consequences of implicated variants.

RESULTS

Our cohort contained 155 154 individuals (3132 IA cases and 152 022 controls). We identified 4 genetic loci reaching genome-wide: rs73392700 ( SIRT3 , effect size = 0.28, P = 4.3 × 10 -12 ), rs58721068 ( EDNRA , effect size = -0.20, P = 4.8 × 10 -12 ), rs4977574 ( AL359922.1 , effect size = 0.18, P = 7.9 × 10 -12 ), and rs11105337 ( ATP2B1 , effect size = -0.15, P = 3.4 × 10 -8 ). Expression quantitative trait loci mapping suggests that rs73392700 has a large effect size on SIRT3 gene expression in arterial and muscle, but not neurological, tissues. Functional genomics analysis suggests that rs73392700 causes decreased SIRT3 gene expression.

CONCLUSION

We perform a multinational GWAS of IA and identify 4 genetic risk loci, including 2 novel IA risk loci ( SIRT3 and AL359922.1 ). Identification of high-risk genetic loci across ancestries will enable population-genetic screening approaches to identify patients with IA.

NLRP3 Susceptible Gene Polymorphisms in Patients with Primary Gouty Arthritis and Hyperuricemia

Polymorphisms have been identified to predispose to primary gouty arthritis (GA) and hyperuricemia (HUA). Here, we accessed the five polymorphisms of rs10754558, rs35829419, rs3738448, rs3806268, and rs7525979 in NLRP3 on GA and HUA susceptibility. We collected 1198 samples (314 GA, 377 HUA, and 507 controls) for this case-control study. Our data detected that the rs3806268 (GA vs. AA: OR = 0.65, p = 0.012) was significantly associated with the susceptibility to GA. The rs3738448 (TT vs. GG: OR = 2.05, p = 0.024) and rs7525979 (TT vs. CC: OR = 1.96, p = 0.037) were significantly associated with the susceptibility to HUA. The rs3806268 AG genotype presented decreased risk of GA among the hypertension (OR = 0.54, p = 0.0093), smoking (OR = 0.59, p = 0.018), and no obesity (OR = 0.60, p = 0.0097) subjects compared to the GG genotype group. The rs3738448 TT genotype demonstrated increased risk of HUA among the hypertension (OR = 4.10, p = 0.0056) and no drinking population (OR = 3.56, p = 0.016) compared to the GG genotype group. The rs7525979 TT genotype demonstrated increased risk of HUA among the hypertension (OR = 4.01, p = 0.0064) and no drinking population (OR = 3.24, p = 0.034) compared to the CC genotype group. Furthermore, a significant haplotype effect of rs10754558/C-rs35829419/C-rs3738448/G-rs3806268/A-rs7525979/C was found (OR = 1.60, p = 0.0046) compared with GCGAC haplotype. Bioinformatics analyses indicated that rs3738448, rs3806268, and rs7525979 might influence the gene regulation, while the T-allele of rs3738448 increased the stability of NLRP3-mRNA. Collectively, our case-control study confirms NLRP3 polymorphisms might participate in regulating immune and inflammation responses in GA and HUA.

Cross-ancestry genome-wide meta-analysis of 61,047 cases and 947,237 controls identifies new susceptibility loci contributing to lung cancer

To identify new susceptibility loci to lung cancer among diverse populations, we performed cross-ancestry genome-wide association studies in European, East Asian and African populations and discovered five loci that have not been previously reported. We replicated 26 signals and identified 10 new lead associations from previously reported loci. Rare-variant associations tended to be specific to populations, but even common-variant associations influencing smoking behavior, such as those with CHRNA5 and CYP2A6, showed population specificity. Fine-mapping and expression quantitative trait locus colocalization nominated several candidate variants and susceptibility genes such as IRF4 and FUBP1. DNA damage assays of prioritized genes in lung fibroblasts indicated that a subset of these genes, including the pleiotropic gene IRF4, potentially exert effects by promoting endogenous DNA damage.

Genome-wide association meta-analysis of 88,250 individuals highlights pleiotropic mechanisms of five ocular diseases in UK Biobank

BACKGROUND

Ocular diseases may exhibit common clinical symptoms and epidemiological comorbidity. However, the extent of pleiotropic mechanisms across ocular diseases remains unclear. We aim to examine shared genetic etiology in age-related macular degeneration (AMD), diabetic retinopathy (DR), glaucoma, retinal detachment (RD), and myopia.

METHODS

We analyzed genome-wide association analyses for the five ocular diseases in 43,877 cases and 44,373 controls of European ancestry from UK Biobank, estimated their genetic relationships (LDSC, GNOVA, and Genomic SEM), and identified pleiotropic loci (ASSET and METASOFT).

FINDINGS

The genetic correlation of common SNPs revealed a meaningful genetic structure within these diseases, identifying genetic correlations between AMD, DR, and glaucoma. Cross-trait meta-analysis identified 23 pleiotropic loci associated with at least two ocular diseases and 14 loci unique to individual disorders (non-pleiotropic). We found that the genes associated with these shared genetic loci are involved in neuron differentiation (P = 8.80 × 10^-6) and eye development systems (P = 3.86 × 10^-5), and single cell RNA sequencing data reveals their heightened gene expression from multipotent progenitors to other differentiated retinal cells during retina developmental process.

INTERPRETATION

These results highlighted the potential common genetic architectures among these ocular diseases and can deepen the understanding of the molecular mechanisms underlying the related diseases.

FUNDING

The National Natural Science Foundation of China (61871294), Zhejiang Provincial Natural Science Foundation of China (LR19C060001), and the Scientific Research Foundation for Talents of Wenzhou Medical University (QTJ18023).

Genetics of the human microglia regulome refines Alzheimer's disease risk loci

Microglia are brain myeloid cells that play a critical role in neuroimmunity and the etiology of Alzheimer's disease (AD), yet our understanding of how the genetic regulatory landscape controls microglial function and contributes to AD is limited. Here, we performed transcriptome and chromatin accessibility profiling in primary human microglia from 150 donors to identify genetically driven variation and cell-specific enhancer-promoter (E-P) interactions. Integrative fine-mapping analysis identified putative regulatory mechanisms for 21 AD risk loci, of which 18 were refined to a single gene, including 3 new candidate risk genes (KCNN4, FIBP and LRRC25). Transcription factor regulatory networks captured AD risk variation and identified SPI1 as a key putative regulator of microglia expression and AD risk. This comprehensive resource capturing variation in the human microglia regulome provides insights into the etiology of neurodegenerative disease.

Multivariate genome-wide association study of depression, cognition, and memory phenotypes and validation analysis identify 12 cross-ethnic variants

To date, little is known about the pleiotropic genetic variants among depression, cognition, and memory. The current research aimed to identify the potential pleiotropic single nucleotide polymorphisms (SNPs), genes, and pathways of the three phenotypes by conducting a multivariate genome-wide association study and an additional pleiotropy analysis among Chinese individuals and further validate the top variants in the UK Biobank (UKB). In the discovery phase, the participants were 139 pairs of dizygotic twins from the Qingdao Twins Registry. The genome-wide efficient mixed-model analysis identified 164 SNPs reaching suggestive significance (P < 1 × 10-5). Among them, rs3967317 (P = 1.21 × 10-8) exceeded the genome-wide significance level (P < 5 × 10-8) and was also demonstrated to be associated with depression and memory in pleiotropy analysis, followed by rs9863698, rs3967316, and rs9261381 (P = 7.80 × 10-8-5.68 × 10-7), which were associated with all three phenotypes. After imputation, a total of 457 SNPs reached suggestive significance. The top SNP chr6:24597173 was located in the KIAA0319 gene, which had biased expression in brain tissues. Genes and pathways related to metabolism, immunity, and neuronal systems demonstrated nominal significance (P < 0.05) in gene-based and pathway enrichment analyses. In the validation phase, 12 of the abovementioned SNPs reached the nominal significance level (P < 0.05) in the UKB. Among them, three SNPs were located in the KIAA0319 gene, and four SNPs were identified as significant expression quantitative trait loci in brain tissues. These findings may provide evidence for pleiotropic variants among depression, cognition, and memory and clues for further exploring the shared genetic pathogenesis of depression with Alzheimer's disease.

Epigenomic and transcriptomic analyses define core cell types, genes and targetable mechanisms for kidney disease

More than 800 million people suffer from kidney disease, yet the mechanism of kidney dysfunction is poorly understood. In the present study, we define the genetic association with kidney function in 1.5 million individuals and identify 878 (126 new) loci. We map the genotype effect on the methylome in 443 kidneys, transcriptome in 686 samples and single-cell open chromatin in 57,229 kidney cells. Heritability analysis reveals that methylation variation explains a larger fraction of heritability than gene expression. We present a multi-stage prioritization strategy and prioritize target genes for 87% of kidney function loci. We highlight key roles of proximal tubules and metabolism in kidney function regulation. Furthermore, the causal role of SLC47A1 in kidney disease is defined in mice with genetic loss of Slc47a1 and in human individuals carrying loss-of-function variants. Our findings emphasize the key role of bulk and single-cell epigenomic information in translating genome-wide association studies into identifying causal genes, cellular origins and mechanisms of complex traits.

Toxicological Mechanism of Individual Susceptibility to Formaldehyde-Induced Respiratory Effects

Understanding the mechanisms of individual susceptibility to exposure to environmental pollutants has been a challenge in health risk assessment. Here, an integrated approach combining a CRISPR screen in human cells and epidemiological analysis was developed to identify the individual susceptibility to the adverse health effects of air pollutants by taking formaldehyde (FA) and the associated chronic obstructive pulmonary disease (COPD) as a case study. Among the primary hits of CRISPR screening of FA in human A549 cells, HTR4 was the only gene genetically associated with COPD susceptibility in global populations. However, the association between HTR4 and FA-induced respiratory toxicity is unknown in the literature. Adverse outcome pathway (AOP) network analysis of CRISPR screen hits provided a potential mechanistic link between activation of HTR4 (molecular initiating event) and FA-induced lung injury (adverse outcome). Systematic toxicology tests (in vitro and animal experiments) were conducted to reveal the HTR4-involved biological mechanisms underlying the susceptibility to adverse health effects of FA. Functionality and enhanced expression of HTR4 were required for susceptibility to FA-induced lung injury, and FA-induced epigenetic changes could result in enhanced expression of HTR4. Specific epigenetic and genetic characteristics of HTR4 were associated with the progression and prevalence of COPD, respectively, and these genetic risk factors for COPD could be potential biomarkers of individual susceptibility to adverse respiratory effects of FA. These biomarkers could be of great significance for defining subpopulations susceptible to exposure to FA and reducing uncertainty in the next-generation health risk assessment of air pollutants. Our study delineated a novel toxicological pathway mediated by HTR4 in FA-induced lung injury, which could provide a mechanistic understanding of the potential biomarkers of individual susceptibility to adverse respiratory effects of FA.

Genetic characterization of outbred Sprague Dawley rats and utility for genome-wide association studies

Sprague Dawley (SD) rats are among the most widely used outbred laboratory rat populations. Despite this, the genetic characteristics of SD rats have not been clearly described, and SD rats are rarely used for experiments aimed at exploring genotype-phenotype relationships. In order to use SD rats to perform a genome-wide association study (GWAS), we collected behavioral data from 4,625 SD rats that were predominantly obtained from two commercial vendors, Charles River Laboratories and Harlan Sprague Dawley Inc. Using double-digest genotyping-by-sequencing (ddGBS), we obtained dense, high-quality genotypes at 291,438 SNPs across 4,061 rats. This genetic data allowed us to characterize the variation present in Charles River vs. Harlan SD rats. We found that the two populations are highly diverged (F_ST > 0.4). Furthermore, even for rats obtained from the same vendor, there was strong population structure across breeding facilities and even between rooms at the same facility. We performed multiple separate GWAS by fitting a linear mixed model that accounted for population structure and using meta-analysis to jointly analyze all cohorts. Our study examined Pavlovian conditioned approach (PavCA) behavior, which assesses the propensity for rats to attribute incentive salience to reward-associated cues. We identified 46 significant associations for the various metrics used to define PavCA. The surprising degree of population structure among SD rats from different sources has important implications for their use in both genetic and non-genetic studies.

Outbred Sprague Dawley rats are among the most commonly used rats for neuroscience, physiology and pharmacological research; in the year 2020, 4,188 publications contained the keyword “Sprague Dawley”. Rats identified as “Sprague Dawley” are sold by several commercial vendors, including Charles River Laboratories and Harlan Sprague Dawley Inc. (now Envigo). Despite their widespread use, little is known about the genetic diversity of SD. We genotyped more than 4,000 SD rats, which we used for a genome-wide association study (GWAS) and to characterize genetic differences between SD rats from Charles River Laboratories and Harlan. Our analysis revealed extensive population structure both between and within vendors. The GWAS for Pavlovian conditioned approach (PavCA) identified a number of genome-wide significant loci for that complex behavioral trait. Our results demonstrate that, despite sharing an identical name, SD rats that are obtained from different vendors are very different. Future studies should carefully define the exact source of SD rats being used and may exploit their genetic diversity for genetic studies of complex traits.

Investigating the genetic architecture of eye colour in a Canadian cohort

Eye color is highly variable in populations with European ancestry, ranging from low to high quantities of melanin in the iris. Polymorphisms in the HERC2/OCA2 locus have the largest effect on eye color in these populations, although other genomic regions also influence eye color. We performed genome-wide association studies of eye color in a Canadian cohort of European ancestry (N = 5,641) and investigated candidate causal variants. We uncovered several candidate causal signals in the HERC2/OCA2 region, whereas other loci likely harbor a single causal signal. We observed colocalization of eye color signals with the expression or methylation profiles of cultured primary melanocytes. Genetic correlations of eye and hair color suggest high genome-wide pleiotropy, but locus-level differences in the genetic architecture of both traits. Overall, we provide a better picture of the polymorphisms underpinning eye color variation, which may be a consequence of specific molecular processes in the iris melanocytes.

•

Genome-wide association studies of eye color in 5,641 participants

•

Multiple independent candidate causal variants were identified across HERC2/OCA2

•

Single candidate causal variants observed on or near IRF4, SLC24A4, TYR, and TYRP1

•

Colocalization of eye color signals with expression and methylation profiles

UACA locus is associated with breast cancer chemoresistance and survival

Few germline genetic variants have been robustly linked with breast cancer outcomes. We conducted trans-ethnic meta genome-wide association study (GWAS) of overall survival (OS) in 3973 breast cancer patients from the Pathways Study, one of the largest prospective breast cancer survivor cohorts. A locus spanning the UACA gene, a key regulator of tumor suppressor Par-4, was associated with OS in patients taking Par-4 dependent chemotherapies, including anthracyclines and anti-HER2 therapy, at a genome-wide significance level (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P = 1.27 \times 10^{ - 9}$$\end{document}P=1.27×10−9). This association was confirmed in meta-analysis across four independent prospective breast cancer cohorts (combined hazard ratio = 1.84, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P = 1.28 \times 10^{ - 11}$$\end{document}P=1.28×10−11). Transcriptome-wide association study revealed higher UACA gene expression was significantly associated with worse OS (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P = 4.68 \times 10^{ - 7}$$\end{document}P=4.68×10−7). Our study identified the UACA locus as a genetic predictor of patient outcome following treatment with anthracyclines and/or anti-HER2 therapy, which may have clinical utility in formulating appropriate treatment strategies for breast cancer patients based on their genetic makeup.

Identification of genetic variants influencing methylation in brain with pleiotropic effects on psychiatric disorders

Psychiatric disorders affect 29% of the global population at least once in the lifespan, and genetic studies have proved a shared genetic basis among them, although the underlying molecular mechanisms remain largely unknown. DNA methylation plays an important role in complex disorders and, remarkably, enrichment of common genetic variants influencing allele-specific methylation (ASM) has been reported among variants associated with specific psychiatric disorders. In the present study we assessed the contribution of ASM to a set of eight psychiatric disorders by combining genetic, epigenetic and expression data. We interrogated a list of 3896 ASM tagSNPs in the brain in the summary statistics of a cross-disorder GWAS meta-analysis of eight psychiatric disorders from the Psychiatric Genomics Consortium, including more than 162,000 cases and 276,000 controls. We identified 80 SNPs with pleiotropic effects on psychiatric disorders that show an opposite directional effect on methylation and gene expression. These SNPs converge on eight candidate genes: ZSCAN29, ZSCAN31, BTN3A2, DDAH2, HAPLN4, ARTN, FAM109B and NAGA. ZSCAN29 shows the broadest pleiotropic effects, showing associations with five out of eight psychiatric disorders considered, followed by ZSCAN31 and BTN3A2, associated with three disorders. All these genes overlap with CNVs related to cognitive phenotypes and psychiatric traits, they are expressed in the brain, and seven of them have previously been associated with specific psychiatric disorders, supporting our results. To sum up, our integrative functional genomics analysis identified eight psychiatric disease risk genes that impact a broad list of disorders and highlight an etiologic role of SNPs that influence DNA methylation and gene expression in the brain.

Genetic overlap analysis of endometriosis and asthma identifies shared loci implicating sex hormones and thyroid signalling pathways

STUDY QUESTION

Is there a shared genetic or causal association of endometriosis with asthma or what biological mechanisms may underlie their potential relationships?

SUMMARY ANSWER

Our results confirm a significant but non-causal association of endometriosis with asthma implicating shared genetic susceptibility and biological pathways in the mechanisms of the disorders, and potentially, their co-occurrence.

WHAT IS KNOWN ALREADY

Some observational studies have reported a pattern of co-occurring relationship between endometriosis and asthma; however, there is conflicting evidence and the aetiology, as well as the underlying mechanisms of the relationship, remain unclear.

STUDY DESIGN, SIZE, DURATION

We applied multiple statistical genetic approaches in the analysis of well-powered, genome-wide association study (GWAS) summary data to comprehensively assess the relationship of endometriosis with asthma. Endometriosis GWAS from the International Endogene Consortium (IEC, 17 054 cases and 191 858 controls) and asthma GWAS from the United Kingdom Biobank (UKB, 26 332 cases and 375 505 controls) were analysed. Additional asthma data from the Trans-National Asthma Genetic Consortium (TAGC, 19 954 cases and 107 715 controls) were utilized for replication testing.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We assessed single-nucleotide polymorphism (SNP)-level genetic overlap and correlation between endometriosis and asthma using SNP effect concordance analysis (SECA) and linkage disequilibrium score regression analysis (LDSC) methods, respectively. GWAS meta-analysis, colocalization (GWAS-PW), gene-based and pathway-based functional enrichment analysis methods were applied, respectively, to identify SNP loci, genomic regions, genes and biological pathways shared by endometriosis and asthma. Potential causal associations between endometriosis and asthma were assessed using Mendelian randomization (MR) methods.

MAIN RESULTS AND THE ROLE OF CHANCE

SECA revealed significant concordance of SNP risk effects across the IEC endometriosis and the UKB asthma GWAS. Also, LDSC analysis found a positive and significant genetic correlation (rG = 0.16, P = 2.01 × 10-6) between the two traits. GWAS meta-analysis of the IEC endometriosis and UKB asthma GWAS identified 14 genome-wide significant (Pmeta-analysis < 5.0 × 10-8) independent loci, five of which are putatively novel. Three of these loci were consistently replicated using TAGC asthma GWAS and reinforced in colocalization and gene-based analyses. Additional shared genomic regions were identified in the colocalization analysis. MR found no evidence of a significant causal association between endometriosis and asthma. However, combining gene-based association results across the GWAS for endometriosis and asthma, we identified 17 shared genes with a genome-wide significant Fisher's combined P-value (FCPgene) <2.73 × 10-6. Additional analyses (independent gene-based analysis) replicated evidence of gene-level genetic overlap between endometriosis and asthma. Biological mechanisms including 'thyroid hormone signalling', 'abnormality of immune system physiology', 'androgen biosynthetic process' and 'brain-derived neurotrophic factor signalling pathway', among others, were significantly enriched for endometriosis and asthma in a pathway-based analysis.

LARGE SCALE DATA

The GWAS for endometriosis data were sourced from the International Endogen Consortium (IEC) and can be accessed by contacting the consortium. The GWAS data for asthma are freely available online at Lee Lab (https://www.leelabsg.org/resources) and from the Trans-National Asthma Genetic Consortium (TAGC).

LIMITATIONS, REASONS FOR CAUTION

Given we analysed GWAS datasets from mainly European populations, our results may not be generalizable to other ancestries.

WIDER IMPLICATIONS OF THE FINDINGS

This study provides novel insights into mechanisms underpinning endometriosis and asthma, and potentially their observed relationship. Findings support a co-occurring relationship of endometriosis with asthma largely due to shared genetic components. Agents targeting 'selective androgen receptor modulators' may be therapeutically relevant in both disorders. Moreover, SNPs, loci, genes and biological pathways identified in our study provide potential targets for further investigation in endometriosis and asthma.

STUDY FUNDING/COMPETING INTEREST(S)

National Health and Medical Research Council (NHMRC) of Australia (241,944, 339,462, 389,927, 389,875, 389,891, 389,892, 389,938, 443,036, 442,915, 442,981, 496,610, 496,739, 552,485, 552,498, 1,026,033 and 1,050,208), Wellcome Trust (awards 076113 and 085475) and the Lundbeck Foundation (R102-A9118 and R155-2014-1724). All researchers had full independence from the funders. Authors do not have any conflict of interest.

Identification and Validation of Candidate Genes from Genome-Wide Association Studies

Exploiting the statistical associations coming out from a GWAS experiment to identify and validate candidate genes may be potentially difficult and time consuming. To fill the gap between the identification of candidate genes toward their functional validation onto the trait performance, the prioritization of variants underlying the GWAS-associated regions is necessary. In parallel, recent developments in genomics and statistical methods have been achieved notably in human genetic and they are accordingly being adopted in plant breeding toward the study of the genetic architecture of traits to sustain genetic gains. In this chapter, we aim at providing both theoretical and practical aspects underlying three main options including (1) the MetaGWAS analysis, (2) the statistical fine mapping and (3) the integration of functional data toward the identification and validation of candidate genes from a GWAS experiment.

Structural variants are a major source of gene expression differences in humans and often affect multiple nearby genes

Structural variants (SVs) are an important source of human genome diversity, but their functional effects are poorly understood. We mapped 61,668 SVs in 613 individuals from the GTEx project and measured their effects on gene expression. We estimate that common SVs are causal at 2.66% of eQTLs, a 10.5-fold enrichment relative to their abundance in the genome. Duplications and deletions were the most impactful variant types, whereas the contribution of mobile element insertions was small (0.12% of eQTLs, 1.9-fold enriched). Multitissue analysis of eQTLs revealed that gene-altering SVs show more constitutive effects than other variant types, with 62.09% of coding SV-eQTLs active in all tissues with eQTL activity compared with 23.08% of coding SNV- and indel-eQTLs. Noncoding SVs, SNVs and indels show broadly similar patterns. We also identified 539 rare SVs associated with nearby gene expression outliers. Of these, 62.34% are noncoding SVs that affect gene expression but have modest enrichment at regulatory elements, showing that rare noncoding SVs are a major source of gene expression differences but remain difficult to predict from current annotations. Both common and rare SVs often affect the expression of multiple genes: SV-eQTLs affect an average of 1.82 nearby genes, whereas SNV- and indel-eQTLs affect an average of 1.09 genes, and 21.34% of rare expression-altering SVs show effects on two to nine different genes. We also observe significant effects on rare gene expression changes extending 1 Mb from the SV. This provides a mechanism by which individual SVs may have strong or pleiotropic effects on phenotypic variation.

Reviewing and assessing existing meta-analysis models and tools

Over the past few years, meta-analysis has become popular among biomedical researchers for detecting biomarkers across multiple cohort studies with increased predictive power. Combining datasets from different sources increases sample size, thus overcoming the issue related to limited sample size from each individual study and boosting the predictive power. This leads to an increased likelihood of more accurately predicting differentially expressed genes/proteins or significant biomarkers underlying the biological condition of interest. Currently, several meta-analysis methods and tools exist, each having its own strengths and limitations. In this paper, we survey existing meta-analysis methods, and assess the performance of different methods based on results from different datasets as well as assessment from prior knowledge of each method. This provides a reference summary of meta-analysis models and tools, which helps to guide end-users on the choice of appropriate models or tools for given types of datasets and enables developers to consider current advances when planning the development of new meta-analysis models and more practical integrative tools.

A large Canadian cohort provides insights into the genetic architecture of human hair colour

Hair colour is a polygenic phenotype that results from differences in the amount and ratio of melanins located in the hair bulb. Genome-wide association studies (GWAS) have identified many loci involved in the pigmentation pathway affecting hair colour. However, most of the associated loci overlap non-protein coding regions and many of the molecular mechanisms underlying pigmentation variation are still not understood. Here, we conduct GWAS meta-analyses of hair colour in a Canadian cohort of 12,741 individuals of European ancestry. By performing fine-mapping analyses we identify candidate causal variants in pigmentation loci associated with blonde, red and brown hair colour. Additionally, we observe colocalization of several GWAS hits with expression and methylation quantitative trait loci (QTLs) of cultured melanocytes. Finally, transcriptome-wide association studies (TWAS) further nominate the expression of EDNRB and CDK10 as significantly associated with hair colour. Our results provide insights on the mechanisms regulating pigmentation biology in humans.

Performing post-genome-wide association study analysis: overview, challenges and recommendations

Genome-wide association studies (GWAS) provide  huge information on statistically significant single-nucleotide polymorphisms (SNPs) associated with various human complex traits and diseases. By performing GWAS studies, scientists have successfully identified the association of hundreds of thousands to  millions of SNPs to a single phenotype. Moreover, the association of some SNPs with rare diseases has been intensively tested. However, classic GWAS studies have not yet provided solid, knowledgeable insight into functional and biological mechanisms underlying phenotypes or mechanisms of diseases. Therefore, several post-GWAS (pGWAS) methods have been recommended. Currently, there is no simple scientific document to provide a quick guide for performing pGWAS analysis. pGWAS is a crucial step for a better understanding of the biological machinery beyond the SNPs. Here, we provide an overview to performing pGWAS analysis and demonstrate the challenges behind each method. Furthermore, we direct readers to key articles for each pGWAS method and present the overall issues in pGWAS analysis.  Finally, we include a custom pGWAS pipeline to guide new users when performing their research.