Multi-ancestry genome-wide association study of major depression aids locus discovery, fine mapping, gene prioritization and causal inference

Most genome-wide association studies (GWAS) of major depression (MD) have been conducted in samples of European ancestry. Here we report a multi-ancestry GWAS of MD, adding data from 21 cohorts with 88,316 MD cases and 902,757 controls to previously reported data. This analysis used a range of measures to define MD and included samples of African (36% of effective sample size), East Asian (26%) and South Asian (6%) ancestry and Hispanic/Latin American participants (32%). The multi-ancestry GWAS identified 53 significantly associated novel loci. For loci from GWAS in European ancestry samples, fewer than expected were transferable to other ancestry groups. Fine mapping benefited from additional sample diversity. A transcriptome-wide association study identified 205 significantly associated novel genes. These findings suggest that, for MD, increasing ancestral and global diversity in genetic studies may be particularly important to ensure discovery of core genes and inform about transferability of findings.

Genome-wide association studies and cross-population meta-analyses investigating short and long sleep duration

Sleep duration has been linked to a wide range of negative health outcomes and to reduced life expectancy. We present genome-wide association studies of short ( ≤ 5 h) and long ( ≥ 10 h) sleep duration in adults of European (N = 445,966), African (N = 27,785), East Asian (N = 3141), and admixed-American (N = 16,250) ancestry from UK Biobank and the Million Veteran Programme. In a cross-population meta-analysis, we identify 84 independent loci for short sleep and 1 for long sleep. We estimate SNP-based heritability for both sleep traits in each ancestry based on population derived linkage disequilibrium (LD) scores using cov-LDSC. We identify positive genetic correlation between short and long sleep traits (rg = 0.16 ± 0.04; p = 0.0002), as well as similar patterns of genetic correlation with other psychiatric and cardiometabolic phenotypes. Mendelian randomisation reveals a directional causal relationship between short sleep and depression, and a bidirectional causal relationship between long sleep and depression.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Discovery and systematic characterization of risk variants and genes for coronary artery disease in over a million participants

The discovery of genetic loci associated with complex diseases has outpaced the elucidation of mechanisms of disease pathogenesis. Here we conducted a genome-wide association study (GWAS) for coronary artery disease (CAD) comprising 181,522 cases among 1,165,690 participants of predominantly European ancestry. We detected 241 associations, including 30 new loci. Cross-ancestry meta-analysis with a Japanese GWAS yielded 38 additional new loci. We prioritized likely causal variants using functionally informed fine-mapping, yielding 42 associations with less than five variants in the 95% credible set. Similarity-based clustering suggested roles for early developmental processes, cell cycle signaling and vascular cell migration and proliferation in the pathogenesis of CAD. We prioritized 220 candidate causal genes, combining eight complementary approaches, including 123 supported by three or more approaches. Using CRISPR-Cas9, we experimentally validated the effect of an enhancer in MYO9B, which appears to mediate CAD risk by regulating vascular cell motility. Our analysis identifies and systematically characterizes >250 risk loci for CAD to inform experimental interrogation of putative causal mechanisms for CAD.

A saturated map of common genetic variants associated with human height

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Adolescent Verbal Memory as a Psychosis Endophenotype: A Genome-Wide Association Study in an Ancestrally Diverse Sample

Verbal memory impairment is one of the most prominent cognitive deficits in psychosis. However, few studies have investigated the genetic basis of verbal memory in a neurodevelopmental context, and most genome-wide association studies (GWASs) have been conducted in European-ancestry populations. We conducted a GWAS on verbal memory in a maximum of 11,017 participants aged 8.9 to 11.1 years in the Adolescent Brain Cognitive Development Study®, recruited from a diverse population in the United States. Verbal memory was assessed by the Rey Auditory Verbal Learning Test, which included three measures of verbal memory: immediate recall, short-delay recall, and long-delay recall. We adopted a mixed-model approach to perform a joint GWAS of all participants, adjusting for ancestral background and familial relatedness. The inclusion of participants from all ancestries increased the power of the GWAS. Two novel genome-wide significant associations were found for short-delay and long-delay recall verbal memory. In particular, one locus (rs9896243) associated with long-delay recall was mapped to the NSF (N-Ethylmaleimide Sensitive Factor, Vesicle Fusing ATPase) gene, indicating the role of membrane fusion in adolescent verbal memory. Based on the GWAS in the European subset, we estimated the SNP-heritability to be 15% to 29% for the three verbal memory traits. We found that verbal memory was genetically correlated with schizophrenia, providing further evidence supporting verbal memory as an endophenotype for psychosis.

Type 2 diabetes risks and determinants in second-generation migrants and mixed ethnicity people of South Asian and African Caribbean descent in the UK

AIMS/HYPOTHESIS

Excess risks of type 2 diabetes in UK South Asians (SA) and African Caribbeans (AC) compared with Europeans remain unexplained. We studied risks and determinants of type 2 diabetes in first- and second-generation (born in the UK) migrants, and in those of mixed ethnicity.

METHODS

Data from the UK Biobank, a population-based cohort of ~500,000 participants aged 40-69 at recruitment, were used. Type 2 diabetes was assigned using self-report and HbA1c. Ethnicity was both self-reported and genetically assigned using admixture level scores. European, mixed European/South Asian (MixESA), mixed European/African Caribbean (MixEAC), SA and AC groups were analysed, matched for age and sex to enable comparison. In the frames of this cross-sectional study, we compared type 2 diabetes in second- vs first-generation migrants, and mixed ethnicity vs non-mixed groups. Risks and explanations were analysed using logistic regression and mediation analysis, respectively.

RESULTS

Type 2 diabetes prevalence was markedly elevated in SA (599/3317 = 18%) and AC (534/4180 = 13%) compared with Europeans (140/3324 = 4%). Prevalence was lower in second- vs first-generation SA (124/1115 = 11% vs 155/1115 = 14%) and AC (163/2200 = 7% vs 227/2200 = 10%). Favourable adiposity (i.e. lower waist/hip ratio or BMI) contributed to lower risk in second-generation migrants. Type 2 diabetes in mixed populations (MixESA: 52/831 = 6%, MixEAC: 70/1045 = 7%) was lower than in comparator ethnic groups (SA: 18%, AC: 13%) and higher than in Europeans (4%). Greater socioeconomic deprivation accounted for 17% and 42% of the excess type 2 diabetes risk in MixESA and MixEAC compared with Europeans, respectively. Replacing self-reported with genetically assigned ethnicity corroborated the mixed ethnicity analysis.

CONCLUSIONS/INTERPRETATION

Type 2 diabetes risks in second-generation SA and AC migrants are a fifth lower than in first-generation migrants. Mixed ethnicity risks were markedly lower than SA and AC groups, though remaining higher than in Europeans. Distribution of environmental risk factors, largely obesity and socioeconomic status, appears to play a key role in accounting for ethnic differences in type 2 diabetes risk.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

The Genetic Architecture of Depression in Individuals of East Asian Ancestry: A Genome-Wide Association Study

Importance

Most previous genome-wide association studies (GWAS) of depression have used data from individuals of European descent. This limits the understanding of the underlying biology of depression and raises questions about the transferability of findings between populations.

Objective

To investigate the genetics of depression among individuals of East Asian and European descent living in different geographic locations, and with different outcome definitions for depression.

Design, Setting, and Participants

Genome-wide association analyses followed by meta-analysis, which included data from 9 cohort and case-control data sets comprising individuals with depression and control individuals of East Asian descent. This study was conducted between January 2019 and May 2021.

Exposures

Associations of genetic variants with depression risk were assessed using generalized linear mixed models and logistic regression. The results were combined across studies using fixed-effects meta-analyses. These were subsequently also meta-analyzed with the largest published GWAS for depression among individuals of European descent. Additional meta-analyses were carried out separately by outcome definition (clinical depression vs symptom-based depression) and region (East Asian countries vs Western countries) for East Asian ancestry cohorts.

Main Outcomes and Measures

Depression status was defined based on health records and self-report questionnaires.

Results

There were a total of 194 548 study participants (approximate mean age, 51.3 years; 62.8% women). Participants included 15 771 individuals with depression and 178 777 control individuals of East Asian descent. Five novel associations were identified, including 1 in the meta-analysis for broad depression among those of East Asian descent: rs4656484 (β = -0.018, SE = 0.003, P = 4.43x10-8) at 1q24.1. Another locus at 7p21.2 was associated in a meta-analysis restricted to geographically East Asian studies (β = 0.028, SE = 0.005, P = 6.48x10-9 for rs10240457). The lead variants of these 2 novel loci were not associated with depression risk in European ancestry cohorts (β = -0.003, SE = 0.005, P = .53 for rs4656484 and β = -0.005, SE = 0.004, P = .28 for rs10240457). Only 11% of depression loci previously identified in individuals of European descent reached nominal significance levels in the individuals of East Asian descent. The transancestry genetic correlation between cohorts of East Asian and European descent for clinical depression was r = 0.413 (SE = 0.159). Clinical depression risk was negatively genetically correlated with body mass index in individuals of East Asian descent (r = -0.212, SE = 0.084), contrary to findings for individuals of European descent.

Conclusions and Relevance

These results support caution against generalizing findings about depression risk factors across populations and highlight the need to increase the ancestral and geographic diversity of samples with consistent phenotyping.

Performance of Breast Cancer Polygenic Risk Scores in 760 Female CHEK2 Germline Mutation Carriers

BACKGROUND

Genome-wide association studies suggest that the combined effects of breast cancer (BC)-associated single nucleotide polymorphisms (SNPs) can improve BC risk stratification using polygenic risk scores (PRSs). The performance of PRSs in genome-wide association studies-independent clinical cohorts is poorly studied in individuals carrying mutations in moderately penetrant BC predisposition genes such as CHEK2.

METHODS

A total of 760 female CHEK2 mutation carriers were included; 561 women were affected with BC, of whom 74 developed metachronous contralateral BC (mCBC). For PRS calculations, 2 SNP sets covering 77 (SNP set 1, developed for BC risk stratification in women unselected for their BRCA1/2 germline mutation status) and 88 (SNP set 2, developed for BC risk stratification in female BRCA1/2 mutation carriers) BC-associated SNPs were used. All statistical tests were 2-sided.

RESULTS

Both SNP sets provided concordant PRS results at the individual level (r = 0.91, P < 2.20 × 10-16). Weighted cohort Cox regression analyses revealed statistically significant associations of PRSs with the risk for first BC. For SNP set 1, a hazard ratio of 1.71 per SD of the PRS was observed (95% confidence interval = 1.36 to 2.15, P = 3.87 × 10-6). PRSs identify a subgroup of CHEK2 mutation carriers with a predicted lifetime risk for first BC that exceeds the surveillance thresholds defined by international guidelines. Association of PRS with mCBC was examined via Cox regression analysis (SNP set 1 hazard ratio = 1.23, 95% confidence interval = 0.86 to 1.78, P = .26).

CONCLUSIONS

PRSs may be used to personalize risk-adapted preventive measures for women with CHEK2 mutations. Larger studies are required to assess the role of PRSs in mCBC predisposition.

Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals

Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to ~1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequency ≤ 0.01) variant BP associations (P < 5 × 10-8), of which 32 were in new BP-associated loci and 55 were independent BP-associated single-nucleotide variants within known BP-associated regions. Average effects of rare variants (44% coding) were ~8 times larger than common variant effects and indicate potential candidate causal genes at new and known loci (for example, GATA5 and PLCB3). BP-associated variants (including rare and common) were enriched in regions of active chromatin in fetal tissues, potentially linking fetal development with BP regulation in later life. Multivariable Mendelian randomization suggested possible inverse effects of elevated systolic and diastolic BP on large artery stroke. Our study demonstrates the utility of rare-variant analyses for identifying candidate genes and the results highlight potential therapeutic targets.

Trans-ancestry meta-analyses identify rare and common variants associated with blood pressure and hypertension

High blood pressure is a major risk factor for cardiovascular disease and premature death. However, there is limited knowledge on specific causal genes and pathways. To better understand the genetics of blood pressure, we genotyped 242,296 rare, low-frequency and common genetic variants in up to 192,763 individuals and used ∼155,063 samples for independent replication. We identified 30 new blood pressure- or hypertension-associated genetic regions in the general population, including 3 rare missense variants in RBM47, COL21A1 and RRAS with larger effects (>1.5 mm Hg/allele) than common variants. Multiple rare nonsense and missense variant associations were found in A2ML1, and a low-frequency nonsense variant in ENPEP was identified. Our data extend the spectrum of allelic variation underlying blood pressure traits and hypertension, provide new insights into the pathophysiology of hypertension and indicate new targets for clinical intervention.

Genomic variation in the MAP3K5 gene is associated with β-thalassemia disease severity and hydroxyurea treatment efficacy

AIM

In this study we explored the association between genetic variations in MAP3K5 and PDE7B genes, residing on chromosome 6q23, and disease severity in β-hemoglobinopathy patients, as well as the association between these variants with response to hydroxyurea (HU) treatment. Furthermore, we examined MAP3K5 expression in the context of high fetal hemoglobin (HbF) and upon HU treatment in erythroid progenitor cells from healthy and KLF1 haploinsufficient individuals.

MATERIALS & METHODS

For this purpose, we genotyped β-thalassemia intermedia and major patients and healthy controls, as well as a cohort of compound heterozygous sickle cell disease/β-thalassemia patients receiving HU as HbF augmentation treatment. Furthermore, we examined MAP3K5 expression in the context of high HbF and upon HU treatment in erythroid progenitor cells from healthy and KLF1 haploinsufficient individuals.

RESULTS

A short tandem repeat in the MAP3K5 promoter and two intronic MAP3K5 gene variants, as well as a PDE7B variant, are associated with low HbF levels and a severe disease phenotype. Moreover, MAP3K5 mRNA expression levels are altered in the context of high HbF and are affected by the presence of HU. Lastly, the abovementioned MAP3K5 variants are associated with HU treatment efficacy.

CONCLUSION

Our data suggest that these MAP3K5 variants are indicative of β-thalassemia disease severity and response to HU treatment.

Population-specific documentation of pharmacogenomic markers and their allelic frequencies in FINDbase

Aims: Population and ethnic group-specific allele frequencies of pharmacogenomic markers are poorly documented and not systematically collected in structured data repositories. We developed the Frequency of Inherited Disorders Pharmacogenomics database (FINDbase-PGx), a separate module of the FINDbase, aiming to systematically document pharmacogenomic allele frequencies in various populations and ethnic groups worldwide. Materials & methods: We critically collected and curated 214 scientific articles reporting pharmacogenomic markers allele frequencies in various populations and ethnic groups worldwide. Subsequently, in order to host the curated data, support data visualization and data mining, we developed a website application, utilizing Microsoft™ PivotViewer software. Results: Curated allelic frequency data pertaining to 144 pharmacogenomic markers across 14 genes, representing approximately 87,000 individuals from 150 populations worldwide, are currently included in FINDbase-PGx. A user-friendly query interface allows for easy data querying, based on numerous content criteria, such as population, ethnic group, geographical region, gene, drug and rare allele frequency. Conclusion: FINDbase-PGx is a comprehensive database, which, unlike other pharmacogenomic knowledgebases, fulfills the much needed requirement to systematically document pharmacogenomic allelic frequencies in various populations and ethnic groups worldwide.