Interpreting meta-analyses of genome-wide association studies

Genome-wide association study for lung cancer in 6531 African Americans reveals new susceptibility loci

Despite lung cancer affecting all races and ethnicities, disparities are observed in incidence and mortality rates among different ethnic groups in the United States. Non-Hispanic African Americans had a high incidence rate of lung cancer at 55.8 per 100 000 people, as well as the highest death rate at 37.2 per 100 000 people from 2016 to 2020. While previous genome-wide association studies (GWAS) have identified over 45 susceptibility risk loci that influence lung cancer development, few GWAS have investigated the etiology of lung cancer in African Americans. To address this gap in knowledge, we conducted GWAS of lung cancer focused on studying African Americans, comprising 2267 lung cancer cases and 4264 controls. We identified three loci associated with lung cancer, one with lung adenocarcinoma, and four with lung squamous cell carcinoma in this population at the genomic-wide significance level. Among them, three novel loci were identified near VWF at 12p13.31 for overall lung cancer and GACAT3 at 2p24.3 and LMAN1L at 15q24.1 for lung squamous cell carcinoma. In addition, we confirmed previously reported risk loci with known or new lead variants near CHRNA5 at 15q25.1 and CYP2A6 at 19q13.2 associated with lung cancer and TRIP13 at 5p15.33 and ERC1 at 12p13.33 associated with lung squamous cell carcinoma. Further multi-step functional analyses shed light on biological mechanisms underlying these associations of lung cancer in this population. Our study highlights the importance of ancestry-specific studies for the potential alleviation of lung cancer burden in African Americans.

Integrating large-scale meta-GWAS and PigGTEx resources to decipher the genetic basis of 232 complex traits in pigs

Understanding the molecular and cellular mechanisms underlying complex traits in pigs is crucial for enhancing genetic gain via artificial selection and utilizing pigs as models for human disease and biology. Here, we conducted comprehensive genome-wide association studies (GWAS) followed by a cross-breed meta-analysis for 232 complex traits and a within-breed meta-analysis for 12 traits, using 28.3 million imputed sequence variants in 70 328 animals across 14 pig breeds. We identified 6878 quantitative trait loci (QTL) for 139 complex traits. Leveraging the Pig Genotype-Tissue Expression resource, we systematically investigated the biological context and regulatory mechanisms behind these trait-QTLs, ultimately prioritizing 14 829 variant-gene-tissue-trait regulatory circuits. For instance, rs344053754 regulates UGT2B31 expression in the liver and intestines, potentially by modulating enhancer activity, ultimately influencing litter weight at weaning in pigs. Furthermore, we observed conservation of certain genetic and regulatory mechanisms underlying complex traits between humans and pigs. Overall, our cross-breed meta-GWAS in pigs provides invaluable resources and novel insights into the genetic regulatory and evolutionary mechanisms of complex traits in mammals.

Detection of positive selection on depression-associated genes

Although depression significantly impacts fitness, some hypotheses suggest that it may offer a survival benefit. However, there has been limited systematic investigation into the selection pressures acting on genes associated with depression at the genomic level. Here, we conducted comparative genomic analyses and computational molecular evolutionary analyses on 320 depression-associated genes at two levels, i.e., across the primate phylogeny (long timescale selection) and in modern human populations (recent selection). We identified seven genes under positive selection in the human lineage, and 46 genes under positive selection in modern human populations. Most positively selected variants in modern human populations were at UTR regions and non-coding exons, indicating the importance of gene expression regulation in the evolution of depression-associated genes. Positively selected genes are not only related to immune responses, but also function in reproduction and dietary adaptation. Notably, the proportion of depression-associated genes under positive selection was significantly higher than the positively selected genes at the genome-wide average level in African, East Asian, and South Asian populations. We also identified two positively selected loci that happened to be associated with depression in the South Asian population. Our study revealed that depression-associated genes are subject to varying selection pressures across different populations. We suggest that, in precision medicine-particularly in gene therapy-it is crucial to consider the specific functions of genes within distinct populations.

Systematic dissection of pleiotropic loci and critical regulons in excitatory neurons and microglia relevant to neuropsychiatric and ocular diseases

Advancements in single-cell multimodal techniques have greatly enhanced our understanding of disease-relevant loci identified through genome-wide association studies (GWASs). To investigate the biological connections between the eye and brain, we integrated bulk and single-cell multiomic profiles with GWAS summary statistics for eight neuropsychiatric and five ocular diseases. Our analysis uncovered five latent factors explaining 61.7% of the genetic variance across these 13 diseases, revealing diverse correlational patterns among them. We identified 45 pleiotropic loci with 91 candidate genes that contribute to disease risk. By integrating GWAS and single-cell profiles, we implicated excitatory neurons and microglia as key contributors in the eye-brain connections. Polygenic enrichment analysis further identified 15 pleiotropic regulons in excitatory neurons and 16 in microglia that were linked to comorbid conditions. Functionally, excitatory neuron-specific regulons were involved in axon guidance and synaptic activity, while microglia-specific regulons were associated with immune response and cell activation. In sum, these findings underscore the genetic link between psychiatric disorders and ocular diseases.

Multi-ancestry genome-wide association analyses: a comparison of meta- and mega-analyses in the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) study

BACKGROUND

There is increasing need for effective incorporation of high-dimensional genetics data from individuals with varied ancestry in genome-wide association (GWAS) analyses. Classically, multi-ancestry GWAS analyses are performed using statistical meta-analysis to combine results conducted within homogeneous ancestry groups. The emergence of cosmopolitan reference panels makes collective preprocessing of GWAS data possible, but impact on downstream GWAS results in a mega-analysis framework merits investigation. We utilized GWAS data from the multi-national Hyperglycemia and Adverse Pregnancy Outcome Study to investigate differences in GWAS findings using a homogeneous ancestry meta-analysis versus a heterogeneous ancestry mega-analysis pipeline. Maternal fasting and 1-hr glucose and metabolomics measured during a 2-hr 75-gram oral glucose tolerance test during early third trimester pregnancy were evaluated as phenotypes.

RESULTS

For the homogeneous ancestry meta-analysis pipeline, variant data were prepared by identifying sets of individuals with similar ancestry and imputing to ancestry-specific reference panels. GWAS was conducted within each ancestry group and results were combined using random-effects meta-analysis. For the heterogeneous ancestry mega-analysis pipeline, data for all individuals were collectively imputed to the Trans-Omics for Precision Medicine (TOPMed) cosmopolitan reference panel, and GWAS was conducted using a unified mega-analysis. The meta-analysis pipeline identified genome-wide significant associations for 15 variants in a region close to GCK on chromosome 7 with maternal fasting glucose and no significant findings for 1-hr glucose. Associations in this same region were identified using the mega-analysis pipeline, along with a well-documented association at MTNR1B on chromosome 11 with both fasting and 1-hr maternal glucose. For metabolomics analyses, the number of significant findings in the heterogeneous ancestry mega-analysis far exceeded those from the homogeneous ancestry meta-analysis and confirmed many previously documented associations, but genomic inflation factors were much more variable.

CONCLUSIONS

For multi-ancestry GWAS, heterogeneous ancestry mega-analysis generates a rich set of variants for analysis using a cosmopolitan reference panel and results in vastly more significant, biologically credible and previously documented associations than a homogeneous ancestry meta-analysis approach. Genomic inflation factors do indicate that findings from the mega-analysis pipeline may merit cautious interpretation and further follow-up.

The Impact of Ancestry on Genome-Wide Association Studies

Genome-wide association studies (GWAS) are an important tool for the study of complex disease genetics. Decisions regarding the quality control (QC) procedures employed as part of a GWAS can have important implications on the results and their biological interpretation. Many GWAS have been conducted predominantly in cohorts of European ancestry, but many initiatives aim to increase the representation of diverse ancestries in genetic studies. The question of how these data should be combined and the consequences that genetic variation across ancestry groups might have on GWAS results warrants further investigation. In this study, we focus on several commonly used methods for combining genetic data across diverse ancestry groups and the impact these decisions have on the outcome of GWAS summary statistics. We ran GWAS on two binary phenotypes using ancestry-specific, multi-ancestry mega-analysis, and meta-analysis approaches. We found that while multi-ancestry mega-analysis and meta-analysis approaches can aid in identifying signals shared across ancestries, they can diminish the signal of ancestry-specific associations and modify their effect sizes. These results demonstrate the potential impact on downstream post-GWAS analyses and follow-up studies. Decisions regarding how the genetic data are combined has the potential to mask important findings that might serve individuals of ancestries that have been historically underrepresented in genetic studies. New methods that consider ancestry-specific variants in conjunction with the shared variants need to be developed.

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

Hearing difficulty (HD) is a major health burden in older adults. While ageing-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. We analysed a large-scale HD genome-wide association study (GWAS; ntotal = 501 825, 56% females) and GWAS data related to 3935 brain imaging-derived phenotypes (IDPs) assessed in up to 33 224 individuals (52% females) using multiple MRI modalities. To investigate HD pleiotropy with brain structure and function, we conducted genetic correlation, latent causal variable, Mendelian randomization and multivariable generalized linear regression analyses. Additionally, we performed local genetic correlation and multi-trait co-localization analyses to identify genomic regions and loci implicated in the pleiotropic mechanisms shared between HD and brain IDPs. We observed a widespread genetic correlation of HD with 120 IDPs in females, 89 in males and 171 in the sex-combined analysis. The latent causal variable analysis showed that some of these genetic correlations could be due to cause-effect relationships. For seven of them, the causal effects were also confirmed by the Mendelian randomization approach: vessel volume→HD in the sex-combined analysis; hippocampus volume→HD, cerebellum grey matter volume→HD, primary visual cortex volume→HD and HD→fluctuation amplitudes of node 46 in resting-state functional MRI dimensionality 100 in females; global mean thickness→HD and HD→mean orientation dispersion index in superior corona radiata in males. The local genetic correlation analysis identified 13 pleiotropic regions between HD and these seven IDPs. We also observed a co-localization signal for the rs13026575 variant between HD, primary visual cortex volume and SPTBN1 transcriptomic regulation in females. Brain structure and function may have a role in the sex differences in HD predisposition via possible cause-effect relationships and shared regulatory mechanisms.

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.

Evaluating the Contribution of Cell Type-Specific Alternative Splicing to Variation in Lipid Levels

BACKGROUND

Genome-wide association studies have identified hundreds of loci associated with lipid levels. However, the genetic mechanisms underlying most of these loci are not well-understood. Recent work indicates that changes in the abundance of alternatively spliced transcripts contribute to complex trait variation. Consequently, identifying genetic loci that associate with alternative splicing in disease-relevant cell types and determining the degree to which these loci are informative for lipid biology is of broad interest.

METHODS

We analyze gene splicing in 83 sample-matched induced pluripotent stem cell (iPSC) and hepatocyte-like cell lines (n=166), as well as in an independent collection of primary liver tissues (n=96) to perform discovery of splicing quantitative trait loci (sQTLs).

RESULTS

We observe that transcript splicing is highly cell type specific, and the genes that are differentially spliced between iPSCs and hepatocyte-like cells are enriched for metabolism pathway annotations. We identify 1384 hepatocyte-like cell sQTLs and 1455 iPSC sQTLs at a false discovery rate of <5% and find that sQTLs are often shared across cell types. To evaluate the contribution of sQTLs to variation in lipid levels, we conduct colocalization analysis using lipid genome-wide association data. We identify 19 lipid-associated loci that colocalize either with an hepatocyte-like cell expression quantitative trait locus or sQTL. Only 2 loci colocalize with both a sQTL and expression quantitative trait locus, indicating that sQTLs contribute information about genome-wide association studies loci that cannot be obtained by analysis of steady-state gene expression alone.

CONCLUSIONS

These results provide an important foundation for future efforts that use iPSC and iPSC-derived cells to evaluate genetic mechanisms influencing both cardiovascular disease risk and complex traits in general.

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.

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.

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.

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