High trans-ethnic replicability of GWAS results implies common causal variants

Massively parallel variant-to-function mapping determines functional regulatory variants of non-small cell lung cancer

Genome-wide association studies have identified thousands of genetic variants associated with non-small cell lung cancer (NSCLC), however, it is still challenging to determine the causal variants and to improve disease risk prediction. Here, we applied massively parallel reporter assays to perform NSCLC variant-to-function mapping at scale. A total of 1249 candidate variants were evaluated, and 30 potential causal variants within 12 loci were identified. Accordingly, we proposed three genetic architectures underlying NSCLC susceptibility: multiple causal variants in a single haplotype block (e.g. 4q22.1), multiple causal variants in multiple haplotype blocks (e.g. 5p15.33), and a single causal variant (e.g. 20q11.23). We developed a modified polygenic risk score using the potential causal variants from Chinese populations, improving the performance of risk prediction in 450,821 Europeans from the UK Biobank. Our findings not only augment the understanding of the genetic architecture underlying NSCLC susceptibility but also provide strategy to advance NSCLC risk stratification.

Population molecular genetics in Brazil: From genomic databases and research to the implementation of precision medicine

Precision medicine (PM) stands on the brink of revolutionizing medical practice throughout the world, holding significant potential for enhancing patient outcomes. However, its practical implementation, particularly in resource-limited countries, is not without challenges. The success of PM largely hinges on the availability of extensive datasets, including genetic and genomic information. This paper delves into the PM landscape and the current state of genetic and genomic testing in Brazil. We also shed light on the unique challenges posed by the country's diverse population and discuss ongoing initiatives to tackle these obstacles.

GWAS Meta-analysis of Kidney Function Traits in Japanese Populations

BACKGROUND

Genetic epidemiological evidence for the kidney function traits in East Asian populations, including Japanese, remain still relatively unclarified. Especially, the number of genome-wide association studies (GWASs) for kidney traits reported still remains limited, and the sample size of each independent study is relatively small. Given the genetic variability between ancestries/ethnicities, implementation of GWAS with sufficiently large sample sizes in specific population of Japanese is considered meaningful.

METHODS

We conducted the GWAS meta-analyses of kidney traits by leveraging the GWAS summary data of the representative large genome cohort studies with about 200,000 Japanese participants (n = 202,406 for estimated glomerular filtration rate [eGFR] and n = 200,845 for serum creatinine [SCr]).

RESULTS

In the present GWAS meta-analysis, we identified 110 loci with 169 variants significantly associated with eGFR (on chromosomes 1-13 and 15-22; P < 5 × 10-8), whereas we also identified 112 loci with 176 variants significantly associated with SCr (on chromosomes 1-22; P < 5 × 10-8), of which one locus (more than 1 Mb distant from known loci) with one variant (CD36 rs146148222 on chromosome 7) for SCr was considered as the truly novel finding.

CONCLUSION

The present GWAS meta-analysis of the largest genome cohort studies in Japanese subjects provided some original genomic loci associated with kidney function, which may contribute to the possible development of personalized prevention of kidney diseases based on genomic information in the near future.

Expression of interferon-stimulated genes, but not polymorphisms in the interferon α/β receptor 2 gene, is associated with coronavirus disease 2019 mortality

Excessive inflammatory response is a hallmark of severe COVID-19. This study investigated the associations between interferon-stimulated genes (ISGs) expression, genetic variation in the interferon α/β receptor 2 (IFNAR2) gene, and COVID-19 mortality. We investigated 67 patients with moderate-to-severe COVID-19. Of them, 22 patients (32.8 %) died because of COVID-19. We examined the expression of ISGs in total RNA of peripheral whole blood. We observed a significant increase in the expression of all ISGs examined in non-surviving patients, indicating a heightened interferon type I signaling activation in non-survived patients. Subsequently, we analyzed whether the increase in ISGs expression was correlated with polymorphism within the IFNAR2 gene. Intriguingly, no significant association was observed between IFNAR2 gene polymorphism and COVID-19 mortality. Similarly, no association was noted between the IFNAR2 and ISGs expression levels. Overall, our data showed that higher ISGs expression, which presumably indicates heightened interferon type I activation, is associated with COVID-19 mortality.

Identification of susceptibility loci and relevant cell type for IgA nephropathy in Han Chinese by integrative genome-wide analysis

Although many susceptibility loci for IgA nephropathy (IgAN) have been identified, they only account for 11.0% of the overall IgAN variance. We performed a large genome-wide meta-analysis of IgAN in Han Chinese with 3616 cases and 10 417 controls to identify additional genetic loci of IgAN. Considering that inflammatory bowel disease (IBD) and asthma might share an etiology of dysregulated mucosal immunity with IgAN, we performed cross-trait integrative analysis by leveraging functional annotations of relevant cell type and the pleiotropic information from IBD and asthma. Among 8 669 456 imputed variants, we identified a novel locus at 4p14 containing the long noncoding RNA LOC101060498. Cell type enrichment analysis based on annotations suggested that PMA-I-stimulated CD4+CD25-IL17+ Th17 cell was the most relevant cell type for IgAN, which highlights the essential role of Th17 pathway in the pathogenesis of IgAN. Furthermore, we identified six more novel loci associated with IgAN, which included three loci showing pleiotropic effects with IBD or asthma (2q35/PNKD, 6q25.2/SCAF8, and 22q11.21/UBE2L3) and three loci specific to IgAN (14q32.32/TRAF3, 16q22.2/TXNL4B, and 21q21.3/LINC00113) in the pleiotropic analysis. Our findings support the involvement of mucosal immunity, especially T cell immune response and IL-17 signal pathway, in the development of IgAN and shed light on further investigation of IgAN.

A comprehensive framework for trans-ancestry pathway analysis using GWAS summary data from diverse populations

As more multi-ancestry GWAS summary data become available, we have developed a comprehensive trans-ancestry pathway analysis framework that effectively utilizes this diverse genetic information. Within this framework, we evaluated various strategies for integrating genetic data at different levels-SNP, gene, and pathway-from multiple ancestry groups. Through extensive simulation studies, we have identified robust strategies that demonstrate superior performance across diverse scenarios. Applying these methods, we analyzed 6,970 pathways for their association with schizophrenia, incorporating data from African, East Asian, and European populations. Our analysis identified over 200 pathways significantly associated with schizophrenia, even after excluding genes near genome-wide significant loci. This approach substantially enhances detection efficiency compared to traditional single-ancestry pathway analysis and the conventional approach that amalgamates single-ancestry pathway analysis results across different ancestry groups. Our framework provides a flexible and effective tool for leveraging the expanding pool of multi-ancestry GWAS summary data, thereby improving our ability to identify biologically relevant pathways that contribute to disease susceptibility.

Hierarchical joint analysis of marginal summary statistics-Part I: Multipopulation fine mapping and credible set construction

Recent advancement in genome-wide association studies (GWAS) comes from not only increasingly larger sample sizes but also the shift in focus towards underrepresented populations. Multipopulation GWAS increase power to detect novel risk variants and improve fine-mapping resolution by leveraging evidence and differences in linkage disequilibrium (LD) from diverse populations. Here, we expand upon our previous approach for single-population fine-mapping through Joint Analysis of Marginal SNP Effects (JAM) to a multipopulation analysis (mJAM). Under the assumption that true causal variants are common across studies, we implement a hierarchical model framework that conditions on multiple SNPs while explicitly incorporating the different LD structures across populations. The mJAM framework can be used to first select index variants using the mJAM likelihood with different feature selection approaches. In addition, we present a novel approach leveraging the ideas of mediation to construct credible sets for these index variants. Construction of such credible sets can be performed given any existing index variants. We illustrate the implementation of the mJAM likelihood through two implementations: mJAM-SuSiE (a Bayesian approach) and mJAM-Forward selection. Through simulation studies based on realistic effect sizes and levels of LD, we demonstrated that mJAM performs well for constructing concise credible sets that include the underlying causal variants. In real data examples taken from the most recent multipopulation prostate cancer GWAS, we showed several practical advantages of mJAM over other existing multipopulation methods.

Cross-ancestry genetic architecture and prediction for cholesterol traits

While cholesterol is essential, a high level of cholesterol is associated with the risk of cardiovascular diseases. Genome-wide association studies (GWASs) have proven successful in identifying genetic variants that are linked to cholesterol levels, predominantly in white European populations. However, the extent to which genetic effects on cholesterol vary across different ancestries remains largely unexplored. Here, we estimate cross-ancestry genetic correlation to address questions on how genetic effects are shared across ancestries. We find significant genetic heterogeneity between ancestries for cholesterol traits. Furthermore, we demonstrate that single nucleotide polymorphisms (SNPs) with concordant effects across ancestries for cholesterol are more frequently found in regulatory regions compared to other genomic regions. Indeed, the positive genetic covariance between ancestries is mostly driven by the effects of the concordant SNPs, whereas the genetic heterogeneity is attributed to the discordant SNPs. We also show that the predictive ability of the concordant SNPs is significantly higher than the discordant SNPs in the cross-ancestry polygenic prediction. The list of concordant SNPs for cholesterol is available in GWAS Catalog. These findings have relevance for the understanding of shared genetic architecture across ancestries, contributing to the development of clinical strategies for polygenic prediction of cholesterol in cross-ancestral settings.

Positive selection in the genomes of two Papua New Guinean populations at distinct altitude levels

Highlanders and lowlanders of Papua New Guinea have faced distinct environmental stress, such as hypoxia and environment-specific pathogen exposure, respectively. In this study, we explored the top genomics regions and the candidate driver SNPs for selection in these two populations using newly sequenced whole-genomes of 54 highlanders and 74 lowlanders. We identified two candidate SNPs under selection - one in highlanders, associated with red blood cell traits and another in lowlanders, which is associated with white blood cell count - both potentially influencing the heart rate of Papua New Guineans in opposite directions. We also observed four candidate driver SNPs that exhibit linkage disequilibrium with an introgressed haplotype, highlighting the need to explore the possibility of adaptive introgression within these populations. This study reveals that the signatures of positive selection in highlanders and lowlanders of Papua New Guinea align closely with the challenges they face, which are specific to their environments.

Genetic control of DNA methylation is largely shared across European and East Asian populations

DNA methylation is an ideal trait to study the extent of the shared genetic control across ancestries, effectively providing hundreds of thousands of model molecular traits with large QTL effect sizes. We investigate cis DNAm QTLs in three European (n = 3701) and two East Asian (n = 2099) cohorts to quantify the similarities and differences in the genetic architecture across populations. We observe 80,394 associated mQTLs (62.2% of DNAm probes with significant mQTL) to be significant in both ancestries, while 28,925 mQTLs (22.4%) are identified in only a single ancestry. mQTL effect sizes are highly conserved across populations, with differences in mQTL discovery likely due to differences in allele frequency of associated variants and differing linkage disequilibrium between causal variants and assayed SNPs. This study highlights the overall similarity of genetic control across ancestries and the value of ancestral diversity in increasing the power to detect associations and enhancing fine mapping resolution.

Incorporating genetic similarity of auxiliary samples into eGene identification under the transfer learning framework

BACKGROUND

The term eGene has been applied to define a gene whose expression level is affected by at least one independent expression quantitative trait locus (eQTL). It is both theoretically and empirically important to identify eQTLs and eGenes in genomic studies. However, standard eGene detection methods generally focus on individual cis-variants and cannot efficiently leverage useful knowledge acquired from auxiliary samples into target studies.

METHODS

We propose a multilocus-based eGene identification method called TLegene by integrating shared genetic similarity information available from auxiliary studies under the statistical framework of transfer learning. We apply TLegene to eGene identification in ten TCGA cancers which have an explicit relevant tissue in the GTEx project, and learn genetic effect of variant in TCGA from GTEx. We also adopt TLegene to the Geuvadis project to evaluate its usefulness in non-cancer studies.

RESULTS

We observed substantial genetic effect correlation of cis-variants between TCGA and GTEx for a larger number of genes. Furthermore, consistent with the results of our simulations, we found that TLegene was more powerful than existing methods and thus identified 169 distinct candidate eGenes, which was much larger than the approach that did not consider knowledge transfer across target and auxiliary studies. Previous studies and functional enrichment analyses provided empirical evidence supporting the associations of discovered eGenes, and it also showed evidence of allelic heterogeneity of gene expression. Furthermore, TLegene identified more eGenes in Geuvadis and revealed that these eGenes were mainly enriched in cells EBV transformed lymphocytes tissue.

CONCLUSION

Overall, TLegene represents a flexible and powerful statistical method for eGene identification through transfer learning of genetic similarity shared across auxiliary and target studies.

Genetic Architecture of Neurological Disorders and Their Endophenotypes: Insights from Genetic Association Studies

Population-scale genetic association studies of complex neurologic diseases have identified the underlying genetic architecture as multifactorial. Despite the study sample sizes reaching the millions, the identified disease-related genes explain only a small fraction of the phenotypic variance. Notable advancements in statistical methods now enable researchers to gain insights even from genomic regions where genotype-phenotype associations do not reach statistical significance. Such studies confirm a highly interconnected molecular network comprising a core group of genes directly involved in the disease process, alongside an expanded peripheral network, each contributing a small but potentially important (modulatory) effect. Additionally, causal inference methods, utilizing genetic instruments, have shed light on putative causal links between risk factors and clinical endpoints. In light of the pervasive genetic overlap or pleiotropy, however, caution is warranted in interpreting causal relationships inferred from these analyses. In this chapter, I will introduce the genetic association model, provide insights into the current state of genetic association studies, and discuss potential future directions.

Identifying QTLs involved in hybrid performance and heterotic group complementarity: new GWAS models applied to factorial and admixed diallel maize hybrid panels

KEY MESSAGE

An original GWAS model integrating the ancestry of alleles was proposed and allowed the detection of background specific additive and dominance QTLs involved in heterotic group complementarity and hybrid performance. Maize genetic diversity is structured into genetic groups selected and improved relative to each other. This process increases group complementarity and differentiation over time and ensures that the hybrids produced from inter-group crosses exhibit high performances and heterosis. To identify loci involved in hybrid performance and heterotic group complementarity, we introduced an original association study model that disentangles allelic effects from the heterotic group origin of the alleles and compared it with a conventional additive/dominance model. This new model was applied on a factorial between Dent and Flint lines and a diallel between Dent-Flint admixed lines with two different layers of analysis: within each environment and in a multiple-environment context. We identified several strong additive QTLs for all traits, including some well-known additive QTLs for flowering time (in the region of Vgt1/2 on chromosome 8). Yield trait displayed significant non-additive effects in the diallel panel. Most of the detected Yield QTLs exhibited overdominance or, more likely, pseudo-overdominance effects. Apparent overdominance at these QTLs contributed to a part of the genetic group complementarity. The comparison between environments revealed a higher stability of additive QTL effects than non-additive ones. Several QTLs showed variations of effects according to the local heterotic group origin. We also revealed large chromosomic regions that display genetic group origin effects. Altogether, our results illustrate how admixed panels combined with dedicated GWAS modeling allow the identification of new QTLs that could not be revealed by a classical hybrid panel analyzed with traditional modeling.

The phenotype-genotype reference map: Improving biobank data science through replication

Population-scale biobanks linked to electronic health record data provide vast opportunities to extend our knowledge of human genetics and discover new phenotype-genotype associations. Given their dense phenotype data, biobanks can also facilitate replication studies on a phenome-wide scale. Here, we introduce the phenotype-genotype reference map (PGRM), a set of 5,879 genetic associations from 523 GWAS publications that can be used for high-throughput replication experiments. PGRM phenotypes are standardized as phecodes, ensuring interoperability between biobanks. We applied the PGRM to five ancestry-specific cohorts from four independent biobanks and found evidence of robust replications across a wide array of phenotypes. We show how the PGRM can be used to detect data corruption and to empirically assess parameters for phenome-wide studies. Finally, we use the PGRM to explore factors associated with replicability of GWAS results.

Strategies for the Genomic Analysis of Admixed Populations

Admixed populations constitute a large portion of global human genetic diversity, yet they are often left out of genomics analyses. This exclusion is problematic, as it leads to disparities in the understanding of the genetic structure and history of diverse cohorts and the performance of genomic medicine across populations. Admixed populations have particular statistical challenges, as they inherit genomic segments from multiple source populations-the primary reason they have historically been excluded from genetic studies. In recent years, however, an increasing number of statistical methods and software tools have been developed to account for and leverage admixture in the context of genomics analyses. Here, we provide a survey of such computational strategies for the informed consideration of admixture to allow for the well-calibrated inclusion of mixed ancestry populations in large-scale genomics studies, and we detail persisting gaps in existing tools.

Genetic Risk Factors Related to Coronary Artery Disease and Role of Transforming Growth Factor Beta 1 Polymorphisms

Coronary artery disease (CAD) is one of the leading causes of mortality globally and has long been known to be heritable; however, the specific genetic factors involved have yet to be identified. Recent advances have started to unravel the genetic architecture of this disease and set high expectations about the future use of novel susceptibility variants for its prevention, diagnosis, and treatment. In the past decade, there has been major progress in this area. New tools, like common variant association studies, genome-wide association studies, meta-analyses, and genetic risk scores, allow a better understanding of the genetic risk factors driving CAD. In recent years, researchers have conducted further studies that confirmed the role of numerous genetic factors in the development of CAD. These include genes that affect lipid and carbohydrate metabolism, regulate the function of the endothelium and vascular smooth muscles, influence the coagulation system, or affect the immune system. Many CAD-associated single-nucleotide polymorphisms have been identified, although many of their functions are largely unknown. The inflammatory process that occurs in the coronary vessels is very important in the development of CAD. One important mediator of inflammation is TGFβ1. TGFβ1 plays an important role in the processes leading to CAD, such as by stimulating macrophage and fibroblast chemotaxis, as well as increasing extracellular matrix synthesis. This review discusses the genetic risk factors related to the development of CAD, with a particular focus on polymorphisms of the transforming growth factor β (TGFβ) gene and its receptor.

Fine-mapping across diverse ancestries drives the discovery of putative causal variants underlying human complex traits and diseases

Genome-wide association studies (GWAS) of human complex traits or diseases often implicate genetic loci that span hundreds or thousands of genetic variants, many of which have similar statistical significance. While statistical fine-mapping in individuals of European ancestries has made important discoveries, cross-population fine-mapping has the potential to improve power and resolution by capitalizing on the genomic diversity across ancestries. Here we present SuSiEx, an accurate and computationally efficient method for cross-population fine-mapping, which builds on the single-population fine-mapping framework, Sum of Single Effects (SuSiE). SuSiEx integrates data from an arbitrary number of ancestries, explicitly models population-specific allele frequencies and LD patterns, accounts for multiple causal variants in a genomic region, and can be applied to GWAS summary statistics. We comprehensively evaluated SuSiEx using simulations, a range of quantitative traits measured in both UK Biobank and Taiwan Biobank, and schizophrenia GWAS across East Asian and European ancestries. In all evaluations, SuSiEx fine-mapped more association signals, produced smaller credible sets and higher posterior inclusion probability (PIP) for putative causal variants, and captured population-specific causal variants.

Evaluating significance of European-associated index SNPs in the East Asian population for 31 complex phenotypes

BACKGROUND

Genome-wide association studies (GWASs) have identified many single-nucleotide polymorphisms (SNPs) associated with complex phenotypes in the European (EUR) population; however, the extent to which EUR-associated SNPs can be generalized to other populations such as East Asian (EAS) is not clear.

RESULTS

By leveraging summary statistics of 31 phenotypes in the EUR and EAS populations, we first evaluated the difference in heritability between the two populations and calculated the trans-ethnic genetic correlation. We observed the heritability estimates of some phenotypes varied substantially across populations and 53.3% of trans-ethnic genetic correlations were significantly smaller than one. Next, we examined whether EUR-associated SNPs of these phenotypes could be identified in EAS using the trans-ethnic false discovery rate method while accounting for winner's curse for SNP effect in EUR and difference of sample sizes in EAS. We found on average 54.5% of EUR-associated SNPs were also significant in EAS. Furthermore, we discovered non-significant SNPs had higher effect heterogeneity, and significant SNPs showed more consistent linkage disequilibrium and allele frequency patterns between the two populations. We also demonstrated non-significant SNPs were more likely to undergo natural selection.

CONCLUSIONS

Our study revealed the extent to which EUR-associated SNPs could be significant in the EAS population and offered deep insights into the similarity and diversity of genetic architectures underlying phenotypes in distinct ancestral groups.

Similarity and diversity of genetic architecture for complex traits between East Asian and European populations

BACKGROUND

Genome-wide association studies have detected a large number of single-nucleotide polymorphisms (SNPs) associated with complex traits in diverse ancestral groups. However, the trans-ethnic similarity and diversity of genetic architecture is not well understood currently.

RESULTS

By leveraging summary statistics of 37 traits from East Asian (N_max=254,373) or European (N_max=693,529) populations, we first evaluated the trans-ethnic genetic correlation (ρ_g) and found substantial evidence of shared genetic overlap underlying these traits between the two populations, with [Formula: see text] ranging from 0.53 (se = 0.11) for adult-onset asthma to 0.98 (se = 0.17) for hemoglobin A1c. However, 88.9% of the genetic correlation estimates were significantly less than one, indicating potential heterogeneity in genetic effect across populations. We next identified common associated SNPs using the conjunction conditional false discovery rate method and observed 21.7% of trait-associated SNPs can be identified simultaneously in both populations. Among these shared associated SNPs, 20.8% showed heterogeneous influence on traits between the two ancestral populations. Moreover, we demonstrated that population-common associated SNPs often exhibited more consistent linkage disequilibrium and allele frequency pattern across ancestral groups compared to population-specific or null ones. We also revealed population-specific associated SNPs were much likely to undergo natural selection compared to population-common associated SNPs.

CONCLUSIONS

Our study provides an in-depth understanding of similarity and diversity regarding genetic architecture for complex traits across diverse populations, and can assist in trans-ethnic association analysis, genetic risk prediction, and causal variant fine mapping.

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

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

Novel Genetic Variants Associated with Chronic Kidney Disease Progression

SIGNIFICANCE STATEMENT

eGFR slope has been used as a surrogate outcome for progression of CKD. However, genetic markers associated with eGFR slope among patients with CKD were unknown. We aimed to identify genetic susceptibility loci associated with eGFR slope. A two-phase genome-wide association study identified single nucleotide polymorphisms (SNPs) in TPPP and FAT1-LINC02374 , and 22 of them were used to derive polygenic risk scores that mark the decline of eGFR by disrupting binding of nearby transcription factors. This work is the first to identify the impact of TPPP and FAT1-LINC02374 on CKD progression, providing predictive markers for the decline of eGFR in patients with CKD.

BACKGROUND

The incidence of CKD is associated with genetic factors. However, genetic markers associated with the progression of CKD have not been fully elucidated.

METHODS

We conducted a genome-wide association study among 1738 patients with CKD, mainly from the KoreaN cohort study for Outcomes in patients With CKD. The outcome was eGFR slope. We performed a replication study for discovered single nucleotide polymorphisms (SNPs) with P <10 -6 in 2498 patients with CKD from the Chronic Renal Insufficiency Cohort study. Several expression quantitative trait loci (eQTL) studies, pathway enrichment analyses, exploration of epigenetic architecture, and predicting disruption of transcription factor (TF) binding sites explored potential biological implications of the loci. We developed and evaluated the effect of polygenic risk scores (PRS) on incident CKD outcomes.

RESULTS

SNPs in two novel loci, TPPP and FAT1-LINC02374 , were replicated (rs59402340 in TPPP , Pdiscovery =7.11×10 -7 , PCRIC =8.13×10 -4 , Pmeta =7.23×10 -8 ; rs28629773 in FAT1-LINC02374 , Pdiscovery =6.08×10 -7 , PCRIC =4.33×10 -2 , Pmeta =1.87×10 -7 ). The eQTL studies revealed that the replicated SNPs regulated the expression level of nearby genes associated with kidney function. Furthermore, these SNPs were near gene enhancer regions and predicted to disrupt the binding of TFs. PRS based on the independently significant top 22 SNPs were significantly associated with CKD outcomes.

CONCLUSIONS

This study demonstrates that SNP markers in the TPPP and FAT1-LINC02374 loci could be predictive markers for the decline of eGFR in patients with CKD.

Causal effects on complex traits are similar for common variants across segments of different continental ancestries within admixed individuals

Individuals of admixed ancestries (for example, African Americans) inherit a mosaic of ancestry segments (local ancestry) originating from multiple continental ancestral populations. This offers the unique opportunity of investigating the similarity of genetic effects on traits across ancestries within the same population. Here we introduce an approach to estimate correlation of causal genetic effects (radmix) across local ancestries and analyze 38 complex traits in African-European admixed individuals (N = 53,001) to observe very high correlations (meta-analysis radmix = 0.95, 95% credible interval 0.93-0.97), much higher than correlation of causal effects across continental ancestries. We replicate our results using regression-based methods from marginal genome-wide association study summary statistics. We also report realistic scenarios where regression-based methods yield inflated heterogeneity-by-ancestry due to ancestry-specific tagging of causal effects, and/or polygenicity. Our results motivate genetic analyses that assume minimal heterogeneity in causal effects by ancestry, with implications for the inclusion of ancestry-diverse individuals in studies.

Cross-Ancestry DNA Methylation Marks of Insulin Resistance in Pregnancy: An Integrative Epigenome-Wide Association Study

Although there are some epigenome-wide association studies (EWAS) of insulin resistance, for most of them authors did not replicate their findings, and most are focused on populations of European ancestry, limiting the generalizability. In the Epigenetics in Pregnancy (EPIPREG; n = 294 Europeans and 162 South Asians) study, we conducted an EWAS of insulin resistance in maternal peripheral blood leukocytes, with replication in the Born in Bradford (n = 879; n = 430 Europeans and 449 South Asians), Methyl Epigenome Network Association (MENA) (n = 320), and Botnia (n = 56) cohorts. In EPIPREG, we identified six CpG sites inversely associated with insulin resistance across ancestry, of which five were replicated in independent cohorts (cg02988288, cg19693031, and cg26974062 in TXNIP; cg06690548 in SLC7A11; and cg04861640 in ZSCAN26). From methylation quantitative trait loci analysis in EPIPREG, we identified gene variants related to all five replicated cross-ancestry CpG sites, which were associated with several cardiometabolic phenotypes. Mediation analyses suggested that the gene variants regulate insulin resistance through DNA methylation. To conclude, our cross-ancestry EWAS identified five CpG sites related to lower insulin resistance.

The genetic and evolutionary basis of gene expression variation in East Africans

BACKGROUND

Mapping of quantitative trait loci (QTL) associated with molecular phenotypes is a powerful approach for identifying the genes and molecular mechanisms underlying human traits and diseases, though most studies have focused on individuals of European descent. While important progress has been made to study a greater diversity of human populations, many groups remain unstudied, particularly among indigenous populations within Africa. To better understand the genetics of gene regulation in East Africans, we perform expression and splicing QTL mapping in whole blood from a cohort of 162 diverse Africans from Ethiopia and Tanzania. We assess replication of these QTLs in cohorts of predominantly European ancestry and identify candidate genes under selection in human populations.

RESULTS

We find the gene regulatory architecture of African and non-African populations is broadly shared, though there is a considerable amount of variation at individual loci across populations. Comparing our analyses to an equivalently sized cohort of European Americans, we find that QTL mapping in Africans improves the detection of expression QTLs and fine-mapping of causal variation. Integrating our QTL scans with signatures of natural selection, we find several genes related to immunity and metabolism that are highly differentiated between Africans and non-Africans, as well as a gene associated with pigmentation.

CONCLUSION

Extending QTL mapping studies beyond European ancestry, particularly to diverse indigenous populations, is vital for a complete understanding of the genetic architecture of human traits and can reveal novel functional variation underlying human traits and disease.

An empirical Bayes approach to improving population-specific genetic association estimation by leveraging cross-population data

Populations of non-European ancestry are substantially underrepresented in genome-wide association studies (GWAS). As genetic effects can differ between ancestries due to possibly different causal variants or linkage disequilibrium patterns, a meta-analysis that includes GWAS of all populations yields biased estimation in each of the populations and the bias disproportionately impacts non-European ancestry populations. This is because meta-analysis combines study-specific estimates with inverse variance as the weights, which causes biases towards studies with the largest sample size, typical of the European ancestry population. In this paper, we propose two empirical Bayes (EB) estimators to borrow the strength of information across populations although accounting for between-population heterogeneity. Extensive simulation studies show that the proposed EB estimators are largely unbiased and improve efficiency compared to the population-specific estimator. In contrast, even though the meta-analysis estimator has a much smaller variance, it yields significant bias when the genetic effect is heterogeneous across populations. We apply the proposed EB estimators to a large-scale trans-ancestry GWAS of stroke and demonstrate that the EB estimators reduce the variance of the population-specific estimator substantially, with the effect estimates close to the population-specific estimates.

Multi-ancestry transcriptome-wide association analyses yield insights into tobacco use biology and drug repurposing

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

Impact of cross-ancestry genetic architecture on GWAS in admixed populations

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

Genome-wide association study for systemic lupus erythematosus in an egyptian population

Systemic lupus erythematosus (SLE) susceptibility has a strong genetic component. Genome-wide association studies (GWAS) across trans-ancestral populations show both common and distinct genetic variants of susceptibility across European and Asian ancestries, while many other ethnic populations remain underexplored. We conducted the first SLE GWAS on Egyptians-an admixed North African/Middle Eastern population-using 537 patients and 883 controls. To identify novel susceptibility loci and replicate previously known loci, we performed imputation-based association analysis with 6,382,276 SNPs while accounting for individual admixture. We validated the association analysis using adaptive permutation tests (n = 109). We identified a novel genome-wide significant locus near IRS1/miR-5702 (Pcorrected = 1.98 × 10-8) and eight novel suggestive loci (Pcorrected < 1.0 × 10-5). We also replicated (Pperm < 0.01) 97 previously known loci with at least one associated nearby SNP, with ITGAM, DEF6-PPARD and IRF5 the top three replicated loci. SNPs correlated (r 2 > 0.8) with lead SNPs from four suggestive loci (ARMC9, DIAPH3, IFLDT1, and ENTPD3) were associated with differential gene expression (3.5 × 10-95 < p < 1.0 × 10-2) across diverse tissues. These loci are involved in cellular proliferation and invasion-pathways prominent in lupus and nephritis. Our study highlights the utility of GWAS in an admixed Egyptian population for delineating new genetic associations and for understanding SLE pathogenesis.

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.

Detecting associated genes for complex traits shared across East Asian and European populations under the framework of composite null hypothesis testing

Background

Detecting trans-ethnic common associated genetic loci can offer important insights into shared genetic components underlying complex diseases/traits across diverse continental populations. However, effective statistical methods for such a goal are currently lacking.

Methods

By leveraging summary statistics available from global-scale genome-wide association studies, we herein proposed a novel genetic overlap detection method called CONTO (COmposite Null hypothesis test for Trans-ethnic genetic Overlap) from the perspective of high-dimensional composite null hypothesis testing. Unlike previous studies which generally analyzed individual genetic variants, CONTO is a gene-centric method which focuses on a set of genetic variants located within a gene simultaneously and assesses their joint significance with the trait of interest. By borrowing the similar principle of joint significance test (JST), CONTO takes the maximum P value of multiple associations as the significance measurement.

Results

Compared to JST which is often overly conservative, CONTO is improved in two aspects, including the construction of three-component mixture null distribution and the adjustment of trans-ethnic genetic correlation. Consequently, CONTO corrects the conservativeness of JST with well-calibrated P values and is much more powerful validated by extensive simulation studies. We applied CONTO to discover common associated genes for 31 complex diseases/traits between the East Asian and European populations, and identified many shared trait-associated genes that had otherwise been missed by JST. We further revealed that population-common genes were generally more evolutionarily conserved than population-specific or null ones.

Conclusion

Overall, CONTO represents a powerful method for detecting common associated genes across diverse ancestral groups; our results provide important implications on the transferability of GWAS discoveries in one population to others.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03637-8.

Importance of Including Non-European Populations in Large Human Genetic Studies to Enhance Precision Medicine

One goal of genomic medicine is to uncover an individual's genetic risk for disease, which generally requires data connecting genotype to phenotype, as done in genome-wide association studies (GWAS). While there may be clinical promise to employing prediction tools such as polygenic risk scores (PRS), it currently stands that individuals of non-European ancestry may not reap the benefits of genomic medicine because of underrepresentation in large-scale genetics studies. Here, we discuss why this inequity poses a problem for genomic medicine and the reasons for the low transferability of PRS across populations. We also survey the ancestry representation of published GWAS and investigate how estimates of ancestry diversity in GWASparticipants might be biased. We highlight the importance of expanding genetic research in Africa, one of the most underrepresented regions in human genomics research, and discuss issues of ethics, resources, and technology for equitable advancement of genomic medicine.

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.

DETECTING MULTIPLE REPLICATING SIGNALS USING ADAPTIVE FILTERING PROCEDURES

Replicability is a fundamental quality of scientific discoveries: we are interested in those signals that are detectable in different laboratories, different populations, across time etc. Unlike meta-analysis which accounts for experimental variability but does not guarantee replicability, testing a partial conjunction (PC) null aims specifically to identify the signals that are discovered in multiple studies. In many contemporary applications, for example, comparing multiple high-throughput genetic experiments, a large number M of PC nulls need to be tested simultaneously, calling for a multiple comparisons correction. However, standard multiple testing adjustments on the M PC p -values can be severely conservative, especially when M is large and the signals are sparse. We introduce AdaFilter, a new multiple testing procedure that increases power by adaptively filtering out unlikely candidates of PC nulls. We prove that AdaFilter can control FWER and FDR as long as data across studies are independent, and has much higher power than other existing methods. We illustrate the application of AdaFilter with three examples: microarray studies of Duchenne muscular dystrophy, single-cell RNA sequencing of T cells in lung cancer tumors and GWAS for metabolomics.

A Method for Bridging Population-Specific Genotypes to Detect Gene Modules Associated with Alzheimer's Disease

BACKGROUND

Genome-wide association studies have successfully identified variants associated with multiple conditions. However, generalizing discoveries across diverse populations remains challenging due to large variations in genetic composition. Methods that perform gene expression imputation have attempted to address the transferability of gene discoveries across populations, but with limited success.

METHODS

Here, we introduce a pipeline that combines gene expression imputation with gene module discovery, including a dense gene module search and a gene set variation analysis, to address the transferability issue. Our method feeds association probabilities of imputed gene expression with a selected phenotype into tissue-specific gene-module discovery over protein interaction networks to create higher-level gene modules.

RESULTS

We demonstrate our method's utility in three case-control studies of Alzheimer's disease (AD) for three different race/ethnic populations (Whites, African descent and Hispanics). We discovered 182 AD-associated genes from gene modules shared between these populations, highlighting new gene modules associated with AD.

CONCLUSIONS

Our innovative framework has the potential to identify robust discoveries across populations based on gene modules, as demonstrated in AD.

Mendelian randomization analysis of 37 clinical factors and coronary artery disease in East Asian and European populations

BACKGROUND

Coronary artery disease (CAD) remains the leading cause of mortality worldwide despite enormous efforts devoted to its prevention and treatment. While many genetic loci have been identified to associate with CAD, the intermediate causal risk factors and etiology have not been fully understood. This study assesses the causal effects of 37 heritable clinical factors on CAD in East Asian and European populations.

METHODS

We collected genome-wide association summary statistics of 37 clinical factors from the Biobank Japan (42,793 to 191,764 participants) and the UK Biobank (314,658 to 442,817 participants), paired with summary statistics of CAD from East Asians (29,319 cases and 183,134 controls) and Europeans (91,753 cases and 311,344 controls). These clinical factors covered 12 cardiometabolic traits, 13 hematological indices, 7 hepatological and 3 renal function indices, and 2 serum electrolyte indices. We performed univariable and multivariable Mendelian randomization (MR) analyses in East Asians and Europeans separately, followed by meta-analysis.

RESULTS

Univariable MR analyses identified reliable causal evidence (P < 0.05/37) of 10 cardiometabolic traits (height, body mass index [BMI], blood pressure, glycemic and lipid traits) and 4 other clinical factors related to red blood cells (red blood cell count [RBC], hemoglobin, hematocrit) and uric acid (UA). Interestingly, while generally consistent, we identified population heterogeneity in the causal effects of BMI and UA, with higher effect sizes in East Asians than those in Europeans. After adjusting for cardiometabolic factors in multivariable MR analysis, red blood cell traits (RBC, meta-analysis odds ratio 1.07 per standard deviation increase, 95% confidence interval 1.02-1.13; hemoglobin, 1.10, 1.03-1.16; hematocrit, 1.10, 1.04-1.17) remained significant (P < 0.05), while UA showed an independent causal effect in East Asians only (1.12, 1.06-1.19, P = 3.26×10-5).

CONCLUSIONS

We confirmed the causal effects of 10 cardiometabolic traits on CAD and identified causal risk effects of RBC, hemoglobin, hematocrit, and UA independent of traditional cardiometabolic factors. We found no causal effects for 23 clinical factors, despite their reported epidemiological associations. Our findings suggest the physiology of red blood cells and the level of UA as potential intervention targets for the prevention of CAD.

Rheumatoid arthritis and the risk of major cardiometabolic diseases: a Mendelian randomization study

OBJECTIVE

Rheumatoid arthritis (RA) is suggested to be implicated in the development of cardiometabolic diseases. We conducted a Mendelian randomization (MR) study to assess potential causality for associations of RA with the risk of cardiometabolic diseases, including type 2 diabetes (T2D), coronary artery disease (CAD), and ischaemic stroke.

METHOD

Seventy independent single-nucleotide polymorphisms (SNPs) associated with RA were identified as instrumental variables from a genome-wide association study (GWAS) of 58 284 European subjects. Summary-level data for the associations of the 70 genetic variants with T2D, CAD, and ischaemic stroke were taken from three GWASs with a total of 1 529 131 participants. Inverse-variance weighted (IVW) MR was used in the main analyses.

RESULTS

The main IVW MR analysis showed that genetically determined RA was associated with higher risks of T2D [odds ratio (OR): 1.04, 95% confidence interval (CI) 1.02-1.05; p < 0.001] and CAD (OR: 1.02, 95% CI 1.00-1.03; p = 0.012), but not ischaemic stroke (OR: 1.00, 95% CI 0.99-1.02; p = 0.961). Sensitivity analyses with multiple MR methods confirmed these associations. MR-Egger regression showed no evidence of pleiotropy in the association between genetically determined RA and the risk of T2D, CAD, and ischaemic stroke. Leave-one-out sensitivity analysis showed that the association between genetically determined RA and the risk of T2D, CAD, and ischaemic stroke was not driven by any individual SNP.

CONCLUSION

Genetically determined RA was associated with increased risks of T2D and CAD, suggesting that RA plays a crucial role in the pathogenesis of T2D and CAD.

A multi-population phenome-wide association study of genetically-predicted height in the Million Veteran Program

BACKGROUND

Height has been associated with many clinical traits but whether such associations are causal versus secondary to confounding remains unclear in many cases. To systematically examine this question, we performed a Mendelian Randomization-Phenome-wide association study (MR-PheWAS) using clinical and genetic data from a national healthcare system biobank.

METHODS AND FINDINGS

Analyses were performed using data from the US Veterans Affairs (VA) Million Veteran Program in non-Hispanic White (EA, n = 222,300) and non-Hispanic Black (AA, n = 58,151) adults in the US. We estimated height genetic risk based on 3290 height-associated variants from a recent European-ancestry genome-wide meta-analysis. We compared associations of measured and genetically-predicted height with phenome-wide traits derived from the VA electronic health record, adjusting for age, sex, and genetic principal components. We found 345 clinical traits associated with measured height in EA and an additional 17 in AA. Of these, 127 were associated with genetically-predicted height at phenome-wide significance in EA and 2 in AA. These associations were largely independent from body mass index. We confirmed several previously described MR associations between height and cardiovascular disease traits such as hypertension, hyperlipidemia, coronary heart disease (CHD), and atrial fibrillation, and further uncovered MR associations with venous circulatory disorders and peripheral neuropathy in the presence and absence of diabetes. As a number of traits associated with genetically-predicted height frequently co-occur with CHD, we evaluated effect modification by CHD status of genetically-predicted height associations with risk factors for and complications of CHD. We found modification of effects of MR associations by CHD status for atrial fibrillation/flutter but not for hypertension, hyperlipidemia, or venous circulatory disorders.

CONCLUSIONS

We conclude that height may be an unrecognized but biologically plausible risk factor for several common conditions in adults. However, more studies are needed to reliably exclude horizontal pleiotropy as a driving force behind at least some of the MR associations observed in this study.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV₁] and its ratio to forced vital capacity [FEV₁/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV₁ and FEV₁/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10^-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Whole-genome sequencing of 1,171 elderly admixed individuals from São Paulo, Brazil

As whole-genome sequencing (WGS) becomes the gold standard tool for studying population genomics and medical applications, data on diverse non-European and admixed individuals are still scarce. Here, we present a high-coverage WGS dataset of 1,171 highly admixed elderly Brazilians from a census-based cohort, providing over 76 million variants, of which ~2 million are absent from large public databases. WGS enables identification of ~2,000 previously undescribed mobile element insertions without previous description, nearly 5 Mb of genomic segments absent from the human genome reference, and over 140 alleles from HLA genes absent from public resources. We reclassify and curate pathogenicity assertions for nearly four hundred variants in genes associated with dominantly-inherited Mendelian disorders and calculate the incidence for selected recessive disorders, demonstrating the clinical usefulness of the present study. Finally, we observe that whole-genome and HLA imputation could be significantly improved compared to available datasets since rare variation represents the largest proportion of input from WGS. These results demonstrate that even smaller sample sizes of underrepresented populations bring relevant data for genomic studies, especially when exploring analyses allowed only by WGS.

Whole genome sequencing (WGS) data on non-European and admixed individuals remains scarce. Here, the authors analyse WGS data from 1,171 admixed elderly Brazilians from a census cohort, characterising population-specific genetic variation and exploring the clinical utility of this expanded dataset.

Transethnic analysis of psoriasis susceptibility in South Asians and Europeans enhances fine-mapping in the MHC and genomewide

Because transethnic analysis may facilitate prioritization of causal genetic variants, we performed a genomewide association study (GWAS) of psoriasis in South Asians (SAS), consisting of 2,590 cases and 1,720 controls. Comparison with our existing European-origin (EUR) GWAS showed that effect sizes of known psoriasis signals were highly correlated in SAS and EUR (Spearman ρ = 0.78; p < 2 × 10-14). Transethnic meta-analysis identified two non-MHC psoriasis loci (1p36.22 and 1q24.2) not previously identified in EUR, which may have regulatory roles. For these two loci, the transethnic GWAS provided higher genetic resolution and reduced the number of potential causal variants compared to using the EUR sample alone. We then explored multiple strategies to develop reference panels for accurately imputing MHC genotypes in both SAS and EUR populations and conducted a fine-mapping of MHC psoriasis associations in SAS and the largest such effort for EUR. HLA-C*06 was the top-ranking MHC locus in both populations but was even more prominent in SAS based on odds ratio, disease liability, model fit and predictive power. Transethnic modeling also substantially boosted the probability that the HLA-C*06 protein variant is causal. Secondary MHC signals included coding variants of HLA-C and HLA-B, but also potential regulatory variants of these two genes as well as HLA-A and several HLA class II genes, with effects on both chromatin accessibility and gene expression. This study highlights the shared genetic basis of psoriasis in SAS and EUR populations and the value of transethnic meta-analysis for discovery and fine-mapping of susceptibility loci.

The correlation of substitution effects across populations and generations in the presence of nonadditive functional gene action

Allele substitution effects at quantitative trait loci (QTL) are part of the basis of quantitative genetics theory and applications such as association analysis and genomic prediction. In the presence of nonadditive functional gene action, substitution effects are not constant across populations. We develop an original approach to model the difference in substitution effects across populations as a first order Taylor series expansion from a "focal" population. This expansion involves the difference in allele frequencies and second-order statistical effects (additive by additive and dominance). The change in allele frequencies is a function of relationships (or genetic distances) across populations. As a result, it is possible to estimate the correlation of substitution effects across two populations using three elements: magnitudes of additive, dominance, and additive by additive variances; relationships (Nei's minimum distances or Fst indexes); and assumed heterozygosities. Similarly, the theory applies as well to distinct generations in a population, in which case the distance across generations is a function of increase of inbreeding. Simulation results confirmed our derivations. Slight biases were observed, depending on the nonadditive mechanism and the reference allele. Our derivations are useful to understand and forecast the possibility of prediction across populations and the similarity of GWAS effects.

Advanced genomics and clinical phenotypes in psoriatic arthritis

Psoriatic Arthritis (PsA) is a complex polygenic inflammatory disease showing a variable musculoskeletal involvement in patients with skin psoriasis. PsA coexist in 25-40 % of patients with the dermatological manifestations, but PsA may also predate the appearance of psoriasis. Nonetheless, the immunopathogenesis of psoriasis and PsA manifest significant similarities, with a major role of the individual susceptibility in both cases. Genome wide association studies (GWAS) identified several genes/loci associated with the risk to develop PsA, both dependent and independent of psoriasis. The major challenge is thus represented by the need to translate the identification of functional polymorphisms and other genetics findings into biological mechanisms along with the identification of novel putative drug targets. A functional genomics approach aims to increase GWAS power and recent evidence supports the use of a multilayer process, including eQTL, methylome, chromatin conformation analysis and genome editing to discover novel genes that can be affected by disease-associated variants, such as PsA. The available data have considered PsA as a unique homogeneous clinical entity while the clinical experience supports a wide variability of skin and joint manifestations coexisting in diverse patients with different mechanisms underlying the musculoskeletal and dermatological domains. A better discrimination of the patient features is encouraged by the limited data on functional genomics. We provide herein a review of the latest findings on PsA functional genomics highlighting the exciting developments in the field and how these might lead to a better understanding of gene regulation underpinning disease mechanisms and ultimately refine clinical phenotyping.

Bench Research Informed by GWAS Results

Scientifically interesting as well as practically important phenotypes often belong to the realm of complex traits. To the extent that these traits are hereditary, they are usually 'highly polygenic'. The study of such traits presents a challenge for researchers, as the complex genetic architecture of such traits makes it nearly impossible to utilise many of the usual methods of reverse genetics, which often focus on specific genes. In recent years, thousands of genome-wide association studies (GWAS) were undertaken to explore the relationships between complex traits and a large number of genetic factors, most of which are characterised by tiny effects. In this review, we aim to familiarise 'wet biologists' with approaches for the interpretation of GWAS results, to clarify some issues that may seem counterintuitive and to assess the possibility of using GWAS results in experiments on various complex traits.

Identifying pleiotropic genes for complex phenotypes with summary statistics from a perspective of composite null hypothesis testing

Pleiotropy has important implication on genetic connection among complex phenotypes and facilitates our understanding of disease etiology. Genome-wide association studies provide an unprecedented opportunity to detect pleiotropic associations; however, efficient pleiotropy test methods are still lacking. We here consider pleiotropy identification from a methodological perspective of high-dimensional composite null hypothesis and propose a powerful gene-based method called MAIUP. MAIUP is constructed based on the traditional intersection-union test with two sets of independent P-values as input and follows a novel idea that was originally proposed under the high-dimensional mediation analysis framework. The key improvement of MAIUP is that it takes the composite null nature of pleiotropy test into account by fitting a three-component mixture null distribution, which can ultimately generate well-calibrated P-values for effective control of family-wise error rate and false discover rate. Another attractive advantage of MAIUP is its ability to effectively address the issue of overlapping subjects commonly encountered in association studies. Simulation studies demonstrate that compared with other methods, only MAIUP can maintain correct type I error control and has higher power across a wide range of scenarios. We apply MAIUP to detect shared associated genes among 14 psychiatric disorders with summary statistics and discover many new pleiotropic genes that are otherwise not identified if failing to account for the issue of composite null hypothesis testing. Functional and enrichment analyses offer additional evidence supporting the validity of these identified pleiotropic genes associated with psychiatric disorders. Overall, MAIUP represents an efficient method for pleiotropy identification.

Identifying causal variants by fine mapping across multiple studies

Increasingly large Genome-Wide Association Studies (GWAS) have yielded numerous variants associated with many complex traits, motivating the development of "fine mapping" methods to identify which of the associated variants are causal. Additionally, GWAS of the same trait for different populations are increasingly available, raising the possibility of refining fine mapping results further by leveraging different linkage disequilibrium (LD) structures across studies. Here, we introduce multiple study causal variants identification in associated regions (MsCAVIAR), a method that extends the popular CAVIAR fine mapping framework to a multiple study setting using a random effects model. MsCAVIAR only requires summary statistics and LD as input, accounts for uncertainty in association statistics using a multivariate normal model, allows for multiple causal variants at a locus, and explicitly models the possibility of different SNP effect sizes in different populations. We demonstrate the efficacy of MsCAVIAR in both a simulation study and a trans-ethnic, trans-biobank fine mapping analysis of High Density Lipoprotein (HDL).

The Meaning of "Cause" in Genetics

Causation has multiple distinct meanings in genetics. One reason for this is meaning slippage between two concepts of the gene: Mendelian and molecular. Another reason is that a variety of genetic methods address different kinds of causal relationships. Some genetic studies address causes of traits in individuals, which can only be assessed when single genes follow predictable inheritance patterns that reliably cause a trait. A second sense concerns the causes of trait differences within a population. Whereas some single genes can be said to cause population-level differences, most often these claims concern the effects of many genes. Polygenic traits can be understood using heritability estimates, which estimate the relative influences of genetic and environmental differences to trait differences within a population. Attempts to understand the molecular mechanisms underlying polygenic traits have been developed, although causal inference based on these results remains controversial. Genetic variation has also recently been leveraged as a randomizing factor to identify environmental causes of trait differences. This technique-Mendelian randomization-offers some solutions to traditional epidemiological challenges, although it is limited to the study of environments with known genetic influences.

Genome-Wide Analyses Reveal the Genetic Architecture and Candidate Genes of Indicine, Taurine, Synthetic Crossbreds, and Locally Adapted Cattle in Brazil

Cattle population history, breeding systems, and geographic subdivision may be reflected in runs of homozygosity (ROH), effective population size (N_e), and linkage disequilibrium (LD) patterns. Thus, the assessment of this information has become essential to the implementation of genomic selection on purebred and crossbred cattle breeding programs. In this way, we assessed the genotype of 19 cattle breeds raised in Brazil belonging to taurine, indicine, synthetic crossbreds, and Iberian-derived locally adapted ancestries to evaluate the overall LD decay patterns, N_e, ROH, and breed composition. We were able to obtain a general overview of the genomic architecture of cattle breeds currently raised in Brazil and other tropical countries. We found that, among the evaluated breeds, different marker densities should be used to improve the genomic prediction accuracy and power of genome-wide association studies. Breeds showing low N_e values indicate a recent inbreeding, also reflected by the occurrence of longer ROH, which demand special attention in the matting schemes to avoid extensive inbreeding. Candidate genes (e.g., ABCA7, PENK, SPP1, IFNAR1, IFNAR2, SPEF2, PRLR, LRRTM1, and LRRTM4) located in the identified ROH islands were evaluated, highlighting biological processes involved with milk production, behavior, rusticity, and fertility. Furthermore, we were successful in obtaining the breed composition regarding the taurine and indicine composition using single-nucleotide polymorphism (SNP) data. Our results were able to observe in detail the genomic backgrounds that are present in each breed and allowed to better understand the various contributions of ancestor breeds to the modern breed composition to the Brazilian cattle.

Causal associations of serum matrix metalloproteinase-8 level with ischaemic stroke and ischaemic stroke subtypes: a Mendelian randomization study

BACKGROUND AND PURPOSE

Elevated serum matrix metalloproteinase-8 (MMP-8) concentrations are associated with high risk of vascular disease, but the causality remains unclear. A two-sample Mendelian randomization (MR) study was performed to examine the causal effect of serum MMP-8 concentrations on the risk of ischaemic stroke, ischaemic stroke subtypes and coronary artery disease.

METHODS

Ten independent single-nucleotide polymorphisms related to serum MMP-8 concentrations were identified as instrumental variables from a genome-wide association study of 6049 European subjects. Genetic association estimates for ischaemic stroke were obtained from the Multiancestry Genome-wide Association Study of Stroke consortium with 446,696 European individuals. The inverse-variance weighted method was applied to assess the causal associations of serum MMP-8 with ischaemic stroke and its subtypes in the main analysis.

RESULTS

No significant causal association was observed for MMP-8 levels with total ischaemic stroke, large artery stroke or cardioembolic stroke. Genetically determined 1 - unit higher log-transformed serum MMP-8 concentration was associated with an increased risk of small vessel stroke (odds ratio 1.25; 95% confidence interval 1.12-1.39; p < 0.001). In secondary analysis, a similar adverse impact was reported for MMP-8 on coronary artery disease (odds ratio 1.05; 95% confidence interval 1.01-1.10; p = 0.017). Sensitivity analyses further confirmed the relationship between serum MMP-8 level and small vessel stroke and coronary artery disease. Mendelian randomization Egger regression showed no evidence of pleiotropic bias.

CONCLUSIONS

High serum MMP-8 concentrations were causally associated with increased risks of small vessel stroke and coronary artery disease. The mechanism underlying the effect of serum MMP-8 on the vascular system requires further investigation.

A Genome-Wide Association Study of Novel Genetic Variants Associated With Anthropometric Traits in Koreans

Most previous genome-wide association studies (GWAS) have identified genetic variants associated with anthropometric traits. However, most of the evidence were reported in European populations. Anthropometric traits such as height and body fat distribution are significantly affected by gender and genetic factors. Here we performed GWAS involving 64,193 Koreans to identify the genetic factors associated with anthropometric phenotypes including height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip ratio. We found nine novel single-nucleotide polymorphisms (SNPs) and 59 independent genetic signals in genomic regions that were reported previously. Of the 19 SNPs reported previously, eight genetic variants at RP11-513I15.6 and one genetic variant at the RP11-977G19.10 region and six Asian-specific genetic variants were newly found. We compared our findings with those of previous studies in other populations. Five overlapping genetic regions (PAN2, ANKRD52, RNF41, HGMA1, and C6orf106) had been reported previously but none of the SNPs were independently identified in the current study. Seven of the nine newly found novel loci associated with height in women revealed a statistically significant skeletal expression of quantitative trait loci. Our study provides additional insight into the genetic effects of anthropometric phenotypes in East Asians.