Cross-ancestry genetic architecture and prediction for cholesterol traits
=========================================================================

* Md. Moksedul Momin
* Xuan Zhou
* Elina Hyppönen
* Beben Benyamin
* S. Hong Lee

## Abstract

While cholesterol is essential for human life, a high level of cholesterol is closely linked with the risk of cardiovascular diseases. Genome-wide association studies (GWASs) have been successful to identify genetic variants associated with cholesterol, which have been conducted mostly in white European populations. Consequently, it remains mostly unknown how genetic effects on cholesterol vary across ancestries. Here, we estimate cross-ancestry genetic correlation to address questions on how genetic effects are shared across ancestries for cholesterol. We find significant genetic heterogeneity between ancestries for total- and LDL-cholesterol. Furthermore, we show that single nucleotide polymorphisms (SNPs), which have concordant effects across ancestries for cholesterol, are more frequently found in the regulatory region, compared to the 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 ([https://www.ebi.ac.uk/gwas/](https://www.ebi.ac.uk/gwas/); details are in web resources section). 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

## Introduction

Cholesterol is a type of lipid that is essential for human life, forming an essential structural component of the cell membrane1-3. While cholesterol is necessary for human body to function, too much cholesterol can harm the body. High cholesterol is linked with a high risk of cardiovascular diseases (CVDs), such as coronary heart disease, stroke, and peripheral vascular disease, which are the leading cause of death worldwide4, accounting for 32% of all deaths in 20195. Specifically, elevated low-density lipoprotein (LDL) and decreased high-density lipoprotein (HDL) cholesterols are associated with increased CVD risk6-9. These cholesterol traits are heritable and known to be polygenic6, 10, 11. Reported heritability estimates for total-, LDL- and HDL-cholesterols are typically in the range of 20 to 60%12.

Over the last two decades, genome-wide association studies (GWASs) have successfully identified several genome-wide significant single nucleotide polymorphisms (SNPs) associated with cholesterol traits4, 13-15. While these findings have provided important insights into the genetics of cholesterol, most GWAS for cholesterol to date have been conducted in populations of white European ancestry16-18. Although the number of GWASs representing non-European populations are gradually increasing, they still remain greatly underrepresented in the efforts of gene discovery16, 19. Consequently, how genetic effects on cholesterol vary across ancestries remain mostly unknown20, 21. It is also not clear to what extent the associated genetic variants discovered in European populations are relevant for other ancestries (e.g., South Asian and African ancestries), and if the polygenic risk prediction of cholesterol can be applied across ancestries22-25.

The genetic effects on most complex traits are likely to vary at least to some extent across different ancestry groups26, 27. Cross-ancestry genetic correlation analyses can dissect the shared genetic architecture between diverse ancestries, also allowing to leverage power from diverse sources of information28. While common causal variants for cholesterol are likely to be shared across ancestries, their per-allele effect sizes may depend on allele frequencies that can differ across ancestries due to different evolutionary force such as selection and genetic drift29. Moreover, each ancestry has a unique genetic background that may affect the magnitude and direction of per-allele effect sizes for complex traits such as cholesterol30. It has been reported that the relationship between allele frequency and per-allele effect size varies across different ancestries, which should be properly accounted for. otherwise, the estimation of cross-ancestry genetic correlation can be biased31, 32.

Cross-ancestry genetic prediction can reduce the potential health disparity for non-European populations that are still underrepresented in public genomic databases including GWAS and polygenic risk scores (PRS)33. It is crucial to understand the source of genetic heterogeneity across ancestries in the genetic prediction. In general, it is not likely that SNP effects estimated from a single ancestry group are always applicable to other ancestries, which has a practical relevance. For example, several studies have reported that the predictive ability of complex traits including cholesterol was poor for Africans, East-Asians, South-Asians and Latinos, when using SNP effects estimated in Europeans19, 34, 35. To obtain more reliable cross-ancestry genetic prediction, it may be important to restrict to functionally homogenous genes or common causal variants across ancestries28, 36. We hypothesize that SNPs in strong linkage disequilibrium (LD) with the functionally homogenous genes have concordant effects, i.e., the same direction of SNP effects, across ancestries.

In this study, we estimate cross-ancestry genetic correlation to address the question about how genetic effects are shared across ancestries for cholesterol traits, accounting for the relationship between allele frequency and per-allele effect size31. In the estimation of cross-ancestry genetic correlation, we also investigate the role of concordant SNPs that are derived from comparing SNP effects between two independent GWAS datasets of UK Biobank and Biobank Japan (BBJ). We evaluate the transferability of genetic prediction across different ancestry groups and suggest a list of SNPs that are suitable for the use in polygenic risk prediction in cross-ancestry analyses.

## Results

### Overview of methods

The total numbers of individuals and SNPs for each ancestry after stringent quality control (QC) (Methods) are shown in **Supplementary Table 1**. From the quality-controlled data of 288,837 white British people, we randomly selected 30,000 individuals to be used in the analyses of cross-ancestry genetic correlations. The remaining 258,792 individuals were used as the discovery dataset in the cross-ancestry genetic risk prediction and in the classification of concordant SNPs (referred to as UKBB discovery). In the cross-ancestry genetic analysis of total-, HDL- and LDL-cholesterol, four ancestry groups were included, i.e. the 30,000 white British ancestry group, 26,457 other European, 6,199 south Asian and 6,179 African ancestry groups (**Supplementary Table 1**). We accounted for the relationship between per-allele effect size and allele frequency31, 37 by using trait-specific and ancestry-specific *α* that was explicitly estimated for each trait and each ancestry, using Akaike information criterion (AIC)31, 32. We used the common SNPs for each pair of ancestries to estimate the cross-ancestry genetic correlation, using the bivariate GREML approach38, accounting for the relationship between allele frequency and per-allele effect size31. We further investigated if the set of concordant SNPs, which were derived by comparing UKBB and Biobank Japan discovery GWAS summary statistics for cholesterol, is enriched in the regulatory region, compared to the other genomic regions. The list of concordant SNPs for total-, HDL- and LDL-cholesterol are now available in GWAS catalogue. Cross-ancestry genetic covariance was partitioned, based on the sets of concordant and discordant SNPs, to see how the genetic heterogeneity is attributed to those SNP sets (see **Supplementary Table 4)**. Finally, cross-ancestry polygenic prediction was performed based on the sets of concordant and discordant SNPs.

### Determining trait-specific and ancestry-specific scale factor (*α*) for each ancestry

The scale factor (*α*) can account for the relationship between allele frequency and per-allele effect size, that is, per-allele effect sizes vary, proportional to [*p* (1 − *p*)] *α*, where *p* is the allele frequency32, 39, 40. It is also reported that the scale factor is not uniformly distributed across ancestries, and there may be an optimal *α* value for each specific ancestry group31. Following the previous approach31, 32, we investigated various *α* values ranging between -1 and 0.5 to determine the ancestry specific *α* value of each ancestry group for total-, LDL- and HDL-cholesterol. To determine optimal *α*, we compared the Akaike Information Criteria (AIC) values across different heritability models with various *α* values for each trait and each ancestry **(Figure 1)**. Detailed values of log-likelihood and AIC are provided in **Supplementary Table 6-9**. As expected, optimal *α* values are not uniformly distributed across traits and across ancestries (Figure 1). These identified *α* values are subsequently used in the estimation of cross-ancestry genetic correlations to dissect the shared genetic architecture and investigate genetic heterogeneity across ancestries for the cholesterol related traits.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/02/02/2023.01.31.23285307/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2023/02/02/2023.01.31.23285307/F1)

Figure 1: Determining the optimal ancestry-specific scaling factors (*α*) for each trait.
The *α* value reflects the relationship between allele frequency and per-allele effect size and can vary across ancestries and traits. ΔAIC values are plotted against scaling factors, *α*, for each ancestry group. The lowest AIC (i.e., ΔAIC=0) indicates the best model. The sample sizes are 30,000, 26,457, 6,199, and 6,179 for white British, other European, South Asian, and African ancestry groups, respectively. TC: total-cholesterol, HDL: high-density lipoprotein cholesterol, LDL: low-density lipoprotein cholesterol.

### Heritability (*h**2*) estimates across ancestries

The estimated SNP-based heritabilities of total-, LDL- and HDL-cholesterol are presented in Figure 1. The estimates are significantly different particularly between European and African ancestries. For total-cholesterol, there is a significant difference in SNP-based heritability estimates between African vs. European (*p-*value=4.26e-03), and African vs white British (*p-* value=1.14e-03). Similarly, the estimate of LDL-cholesterol is significantly lower in white British (*p-*value= 1.11e-03) and other European (*p-*value= 5.19e-03) than African ancestry, which agrees with the previous findings based on twin studies41. We also observed significant heterogeneity of SNP-based heritability for the HDL-cholesterol between South Asian and other Europeans, between South Asian and white British.

### Estimated cross-ancestry genetic correlations

The estimated cross-ancestry genetic correlations (*r**g*) for cholesterol traits are presented in **Figure 3**. For total-cholesterol, we observed a genetic heterogeneity between South Asian vs. white British (*r**g*= 0.399; SE= 0.143; *p*-value= 2.65e-05), South Asian vs. other European (*r**g*= 0.353; SE=0.133; *p*-value= 1.14e-06) and South Asian vs. African ancestry (*r**g* = 0.188; SE=0.197; *p*-value= 3.76e-05). There is also a genetic heterogeneity between African vs. white British (*r**g*= 0.473; SE=0.127; *p*-value= 3.33e-05) and African vs. other European ancestry (*r**g*= 0.315; SE=0.122; *p*-value=1.96e-08). In contrast, white British and other European are genetically homogenous (*r**g* =0.954; SE=0.087; *p*-value= 5.96e-01) (**Figure 3 and Supplementary Table 10**). For LDL-cholesterol, results are similar to total-cholesterol. There is a significant genetic heterogeneity between South Asian vs. white British (*r**g* = 0.296; SE=0.155; *p*-value=5.57e-06), South Asian vs. other European (*r**g* = 0.177; SE=0.138; *p*-value=2.46e-09), South Asian vs. African (*r**g* = 0.110; SE=0.190; *p*-value=2.81e-06), and African vs. other European ancestry (*r**g*= 0.409; SE=0.147; *p*-value=2.81e-06) (**Figure 3 and Supplementary Table 11**). As expected, the cross-ancestry genetic correlation between other European and white British was close to 1 (*r**g*= 1.084; SE=0.128; *p*-value=5.12e-01). We did not observe genetic heterogeneity among the pairs of ancestry groups for HDL-cholesterol (**Figure 3 and Supplementary Table 12)**.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/02/02/2023.01.31.23285307/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2023/02/02/2023.01.31.23285307/F2)

Figure 2: Estimated SNP-based heritability across ancestries for cholesterol traits.
The main bars indicate SNP-based heritability estimates, and the error bars indicate 95% confidence intervals. TC= Total-cholesterol, HDL= high-density lipoprotein cholesterol, LDL= low-density lipoprotein cholesterol.

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/02/02/2023.01.31.23285307/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2023/02/02/2023.01.31.23285307/F3)

Figure 3: Estimated cross-ancestry genetic correlations.
The main bars indicate estimated cross-ancestry genetic correlations, and the error bars indicate 95% confidence intervals of the estimates. WB = White British, OE = Other European, SAS = South Asian, AFR = African.

### Genomic partitioning of cross-ancestry genetic covariance using concordant and discordant SNPs between two diverse ancestries

Some genes are functionally homogeneous across ancestries while the other genes may not be36, 42, 43. It can be hypothesised that the functionally homogenous genes are enriched in the regulatory regions, and they contribute more to phenotypic variation within and between ancestries, compared to the other genes. We obtained a set of concordant SNPs (a proxy of functionally homogenous genes) for total-, HDL-, and LDL-cholesterols, by comparing the direction of SNP effects between two diverse ancestries, using the GWAS summary statistics of UK Biobank and Biobank Japan. For the UK Biobank GWAS, we used 258,792 white British individuals who are not overlapping with anyone in the 4 ancestry groups used in our study (white British, other European, South Asian, and African). For the Biobank Japan, we used GWAS summary statistics that are publicly available. In this concordance/discordance analysis, we considered the same HapMap3 SNPs used in the genetic correlation analyses above. The numbers of concordant and discordant SNPs for each pair of ancestries are presented in **Supplementary Table 4**.

First, we quantified if the concordance SNPs are more frequently found in the regulatory or genic region, compared to the other genomic regions for total-cholesterol. **Figure 4** shows that the number of concordant SNPs in the regulatory region is significantly higher than the non-regulatory region (OR= 1.09, *p*-value=2.2e-26 for *p*-value ≤ 1; OR= 1.21, *p*-value=9.2e-16 for *p*-value ≤ 0.05; OR= 1.18, *p*-value=1.6e-06 for *p*-value ≤ 0.01). When selecting SNPs with a genome-wide association (GWA) *p*-value > 0.05 or 0.01, the odds ratio increases (**Figure 4**). Similarly, the number of concordant SNPs in the genic region is significantly higher than the non-genic region (OR= 1.03, *p*-value=1.8e-16 for *p*-value ≤ 1; OR= 1.17, *p*-value=1.6e-32 for *p*-value ≤ 0.05; OR= 1.18, *p*-value=1.8e-06 for *p*-value ≤ 0.01) (**Figure 4**). Similar results were observed when using the HDL- and LDL-cholesterol traits (**Supplementary Figure 1 and 2**).

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/02/02/2023.01.31.23285307/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2023/02/02/2023.01.31.23285307/F4)

Figure 4: A forest plot with odds ratios indicating that concordant SNPs are more frequently found in the regulatory or genic region.
This analysis is for total-cholesterol phenotypes. Error bar represents 95% confidence intervals. The *p*-value of odds ratio indicates that the odds ratio is significantly different from 1. For regulatory or genic region, a genome-wide association (GWA) *p*-value threshold ≤1, 0.05 or 0.01 was used to select a set of concordant and discordant SNPs using UK Biobank GWAS summary statistics for total-cholesterol.

Subsequently, we partitioned genetic covariance components attributed to the two sets of genomic regions (concordant vs. discordant SNPs). We estimated two genomic relationship matrixes (GRM), using the sets of concordant and discordant SNPs, which were simultaneously fitted in a bivariate multiple random-effects model. When considering the set of concordant SNPs, the estimated genetic covariances between other European (OE) vs. south Asian (SAS), white British (WB) vs. African (AFR) and WB vs. OE were significantly higher than the expectation (the proportion of the concordant SNPs) for total-cholesterol (**Figure 5**). On the other hand, the estimated genetic covariances for these pairs of ancestries were significantly lower than the expectation when using discordant SNPs (**Figure 5**). For HDL-cholesterol (**Supplementary Figure 3**) and LDL-cholesterol (**Supplementary Figure 4**), a similar result was observed that the estimated genetic covariances between OE vs. SAS, WB vs. OE and WB vs. SAS were significantly deviated from the expectation. When using SNPs with genome-wide association *p*-values < 0.05 or 0.01 (**Supplementary Table 4**), the estimated genetic covariances due to concordant and discordant SNPs were more significantly deviated from the expectation in general (**Figure 5, Supplementary Figure 3 and 4**). It is also noted that the estimated genetic covariances for the set of discordant SNPs were not higher than zero (Figure 5), implying that the genetic heterogeneity of cholesterol traits across ancestry might be mostly due to the set of discordant SNPs. This also shows that the set of concordant SNPs may be useful in cross-ancestry polygenic risk predictions. The results are similar when genome-wide association p-values from BBJ are used (**Supplementary Figure 5**).

![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/02/02/2023.01.31.23285307/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2023/02/02/2023.01.31.23285307/F5)

Figure 5: Estimated genetic covariances for concordant and discordant SNPs for total-cholesterol.
Concordant and discordant SNPs were derived from the comparison of SNP effects between two independent GWAS datasets of UK Biobank and BBJ. In this concordant or discordant analysis, a set of SNPs with genome-wide association (GWA) *p*-values < 1, 0.05 or 0.01 was used, where the GWA p-values were from UK Biobank GWAS for total-cholesterol. The main bars represent estimated cross-ancestry genetic covariance using the set of genome-wide SNPs, and the error bars indicate 95% confidence intervals. The horizontal dashed line indicates the expected genetic covariance, assuming all SNPs contribute equally to the genetic covariance, i.e., the expected genetic covariance = the estimated total genetic covariance × the proportion of number of concordant SNPs, where the estimated total genetic covariance is based on all the SNPs including both concordant and discordant SNPs. The value with each bar indicates a *p*-value testing the null hypothesis that the estimated genetic covariance is not significantly different from the expectation. WB = White British, OE = Other European, SAS = South Asian, AFR = African.

We further investigated the impact of concordant SNPs in a cross-ancestry polygenic risk prediction. We used the UKBB discovery dataset, which is independent from the four target datasets including white British, other European, South Asian, and African ancestries, to estimate SNP effects and obtain GWAS summary statistics for cholesterol traits. Using the GWAS summary statistics, we constructed polygenic risk scores for the individuals in the target datasets. The predictive ability (*R*2) of polygenic risk scores for total-cholesterol is significantly higher when using the set of concordant SNPs than when using the set of discordant SNPs for both within- and cross-ancestry predictions (**Figure 6, Supplementary Figure 4**) (*p*-values for the difference between concordant and discordant PRS SNPs are 3.8e-33, 2.2e-25, 1.3e-04 and 5.3e-04 for white British, other European, South Asian and African, respectively). Although not significant, *R*2 is slightly higher when using the set of concordant SNPs, compared to when using the total set of SNPs (**Figure 6**), suggesting that including discordant SNPs may have adverse effects on the cross-ancestry risk predictions. When accounting for the proportion of concordant SNPs, a similar result was observed in that concordant SNPs performed better that discordant SNPs in within- and cross-ancestry risk predictions (**Supplementary Figure 5**). A similar finding was observed when using BBJ discovery GWAS summary statistics, i.e., the cross-ancestry prediction accuracy of the concordant SNPs significantly higher than the discordant SNPs (**Supplementary Figure 4**). Interestingly, the concordant SNPs performs notably better than the total set of SNPs when predicting white British, other European and south Asian ancestries (**Supplementary Figure 4**). Results are invariant when considering LDL- and HDL-cholesterol (Supplementary Figures 6-7).

![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/02/02/2023.01.31.23285307/F6.medium.gif)

[Figure 6:](http://medrxiv.org/content/early/2023/02/02/2023.01.31.23285307/F6)

Figure 6: The predictive ability (*R*2) of polygenic risk scores for total-cholesterol when using the set of concordant, discordant, or total SNPs for cross-ancestry risk predictions.
UK Biobank GWAS was used as the discovery dataset (n= 258,792), while target datasets were other European (n=26,457), south Asian (n=6,199) and African ancestry (n=6,179).

**Left panels:** The main bars represent *R*2 values and error bars correspond to 95% confidence interval.

**Right panels:** Dot points represent the differences between *R*2 values, error bars correspond to 95% confidence intervals of the differences, and *p*-values indicate that the differences of *R*2 are significantly different from zero (null hypothesis). P-values was estimated using an R-package (r2redux)45 based on Wald’s test statistics.

For HDL cholesterol, it is notable that the accuracy of cross-ancestry prediction can be higher than within-ancestry prediction (e.g., South Asian vs. White British in Supplementary Figure 6). We further confirmed this result with a clump-and-threshold (C + T) based PRS method (PRSice)**44** and compared the significance of difference (Supplementary Figure 8). It shows that PGS generated from White British GWAS provides a significantly higher predictive accuracy for South Asian (p-value = 6.67e-16) and African ancestry groups (p-value = 7.35e-04), compared to White British. This may have an important implication in genomic medicine for underrepresented non-European populations.

## Discussion

Cholesterol is an essential structural component of the cell membrane, which is necessary for the body to function1, 2. However, the risk of CVD is associated with a high level of cholesterol that can be determined by genetic risk factors4, 46, 47. Although the genetic study of cholesterol has been conducted, it is not clear how genetic effects on cholesterol vary across different ancestries. In this study, we explicitly estimated cross-ancestry genetic correlations to investigate the shared genetic architecture across ancestries for cholesterol. Importantly, we appropriately accounted for the relationship between allele frequency and per-allele effect size by modelling the ancestry-specific scale factor for cholesterol, which can provide more reliable estimates31.

The reliable estimation of cross-ancestry genetic correlation allows us to understand the shared genetic architecture across ancestries, providing crucial information when for various downstream analyses of complex traits such as cross-ancestry GWAS and cross-ancestry polygenic risk score prediction. Moreover, this may inform best practices for cross-ancestry meta-analysis, multi-ancestry disease mapping, and the transferability of epidemiological findings. Our analysis shows that in general, total- and LDL-cholesterol are both genetically heterogeneous across ancestries, whereas HDL-cholesterol is not48. This finding has important implications for the power of cross-ancestry GWASs and cross-ancestry polygenic risk score prediction, which for HDL-cholesterol may be much higher than that for total- and LDL-cholesterols (**Supplementary Figure 6**).

To identify genetic variants that contribute to the genetic heterogeneity, we investigated concordant and discordant SNP sets that were identified by comparing the direction of SNP effects between UK Biobank and Biobank Japan GWAS summary statistics, noting that the two datasets are independent from the four target ancestry groups used in this study. The concordant SNPs may be associated with genes that are functionally homogeneous across ancestries49, and we show in this study that the concordant SNPs are more often located in the regulatory or genic regions, compared to other genomic regions. We also show that such strong genetic heterogeneity across ancestries for cholesterol can be attributed to the discordant SNPs, but not to the concordant SNPs. We provide evidence that the set of concordant SNPs can be useful in the cross-ancestry polygenic risk predictions, which may improve the transferability of polygenic risk scores to clinical practice16, 50, 51.

There are a number of limitations in this study. For determining optimal *α*, we did not consider the relationship between LD and per-allele effect sizes, i.e., as in LDAK-thin model32 that requires a substantial reduction of the number of SNPs. We also acknowledge that the conclusions from cross-ancestry analyses (cross-ancestry correlation and genomic prediction) in this study are restricted to common variants (MAF ≥ 0.01) and HapMap3 SNPs only; as these are robust and reliable for dissecting cross-ancestry genetic architecture52, 53. A moderate sample size (limited power of the data) was used to estimate optimal scale factors (*α*) for south Asian and African populations. Therefore, the genetic heterogeneity needs to be explored with larger sample size. The concordant SNPs were identified by comparing the direction of SNP effects between white British (UKBB) and East Asian (BBJ) populations, because adequate data was not available from other ancestries. When public genomic databases have sufficient resources across ancestries, we can have a finer set of concordant SNPs by comparing SNP effects across various ancestries.

In conclusion, there is a significant genetic heterogeneity between ancestries for total- and LDL-cholesterol, which is mostly driven by the set of discordant SNPs. Interestingly, the concordant SNPs are more frequently found in the regulatory region as annotated by an independent study54, and restricting to concordant SNPs can provide better accuracy for cross-ancestry polygenic prediction for cholesterol. Our findings contribute to knowledge about the genetic architecture of cholesterol that is shared across ancestries. The proposed cross-ancestry polygenic prediction can be potentially useful in clinical practice. Our analysis protocol can be extended to a wide range of other complex traits and diseases.

## Methods

### Ethical statement

We used publicly available from the UK Biobank ([https://www.ukbiobank.ac.uk/](https://www.ukbiobank.ac.uk/)). Science protocol and operational procedures for the UK Biobank have been reviewed and approved by the North-West Multi-Centre Research Ethics Committee (MREC), National Information Governance Board for Health & Social Care (NIGB), and Community Health Index Advisory Group (CHIAG). The UK Biobank has obtained consent from all participants. The access of the UK Biobank data was approved under the reference number 14575 (“Whole genome approaches for dissecting (shared) genetic architecture and individual risk prediction of complex traits in human populations”). Publicly available GWAS summary statistics of Biobank Japan (BBJ) were used, following BBJ’s guidelines ([http://jenger.riken.jp/en/](http://jenger.riken.jp/en/)). The research ethics approval of this study has been obtained from the University of South Australia Human Research Ethics Committee.

### Participants and stratification of ancestries

Data from the UK Biobank contains 501,748 participants recruited between 2006 and 201055. The participants were recruited from 22 assessment centres in England, Wales, and Scotland, ranging in age from 37 to 73 years old56. All the phenotypic data for cholesterol traits under this study are derived from baseline survey. Principal component analysis was applied to the UK Biobank individuals to stratify participants57 into four different ancestries following previous approach31.

### Genotypic data and quality control

We used the second release of the UK biobank ([https://www.ukbiobank.ac.uk/](https://www.ukbiobank.ac.uk/)) genotype data comprising 488,377 individuals and 92,693,895 imputed autosomal SNPs. The individuals were genotyped by Affymetrix UK BiLEVE Axiom array and Affymetrix UK Biobank Axiom® array. Combination of UK10K and Haplotype Reference Consortium (HRC) data were considered as the reference dataset for the imputation of the UK Biobank genotypic dataset58. In this analysis, to dissect the genetic architecture of disease and complex traits, we retained only HapMap3 SNPs in this analysis52, which are also considered robust and reliable for estimating heritability, genetic correlation52, 59. Stringent quality control (QC) procedure was applied to each ancestry to select high quality individuals and high-quality SNPs. SNPs QC criteria include, SNPs excluded with an INFO score (used to indicate the quality of genotype imputation) <0.660-62, call rate <0.95, a MAF <0.01 and a Hardy–Weinberg equilibrium *p*-value <10−4. We also exclude population outliers (individuals outside ±6SD) and related individuals (--rel-cutoff 0.05) using PLINK63.

Individual level QC criteria include samples with genotype missing rate >0.05, gender mismatch (reported gender does not fit with the genetically assigned sex determined from gene data), poor genotype quality or a sex chromosome aneuploidy was excluded from the main analyses. For the ease of computation, we reduced the number individuals in white British ancestry. The total number of individuals and total number of SNPs after QC shown in **Supplementary Table 1**. The number of common SNPs across different pairs of ancestries presented in **Supplementary Table 2** and the number of common SNP for each genomic region (genomic partitioning) between ancestries presented in **Supplementary Table 3**.

### Functional annotation of the genome

The common SNPs between populations were partitioned into genomic region using genomic annotation reported by Gusev et. al.54, where they partitioned the genome into coding, UTR, promoter, intron, DHS and intergenic regions. For the genomic partitioning analysis, we include promoter, coding, UTR, and DHS regions as regulatory regions64-66, and introns (an integral part of a gene)67, 68 and the intergenic regions as non-genic regions. We also partitioned the whole genome into two predefined functional categories as genic (includes SNPs from promoter, coding, untranslated, intron and DHS region) and non-genic regions (intergenic region).

### Concordant and discordant SNP annotation

To identify concordant and discordant SNPs we compared SNP effects between two independent GWAS datasets of white British from UK Biobank and Biobank Japan (BBJ). The BBJ summary statistics data are publicly available ([http://jenger.riken.jp/en/result](http://jenger.riken.jp/en/result)). We excluded SNPs that were ambiguous or had a strand issue. After excluding these SNPs, there were 4,113,630 SNPs that are common between UKBB and BBJ. To determine concordant and discordant SNPs, we compared the direction of SNP effects between white British from UKBB and BBJ. We used only HapMap3 SNPs from 4,113,630 SNPs for concordant and discordant analysis across different ancestry pairs (**Supplementary Table 4**).

There were four possible combinations of direction of SNP effects (beta):

(+beta, +beta) if the SNP effects are positive in both GWAS.

(+beta, -beta) if the SNP effects are positive and negative in the UKBB and BBJ GWAS.

(-beta, +beta) if the SNP effects are negative and positive in the UKBB and BBJ GWAS.

(-beta, -beta) if the SNP effects are negative in both GWAS.

Each SNP should be in one of four possible combinations and belongs to either concordance or discordance. SNPs belonged to ((+beta, +beta) ∪ (-beta, -beta) were considered concordant, otherwise discordant, i.e. (+beta, -beta) ∪ (-beta, +beta).

### Data analysis

#### Phenotypic adjustment of main traits

Prior to model fitting, all cholesterol traits were adjusted for demographic variables, the UK biobank assessment centre (as factor), genotype measurement batch (as factor) and population structure measured by the first 10 principal components (PCs)64, 69 using linear models in *R*-*software* (4.0.3). Demographic variable includes sex, birth year, education, and Townsend deprivation index (**Supplementary Table 5**). Information of educational qualifications converted to education levels (years) for all the UK Biobank individuals70.

### Determining scale factors for GCTA-*α* model

GCTA model assumes all the SNPs has equal contribution to the genetic variance (has no LD weights), whereas LDAK-thin model32 explicitly considers LD among SNPs. The previously recommended and widely used α are -0.50 and -0.125 for GCTA model71 and LDAK-thin model32, respectively. Here we have used 13 different values of *α* (between -1 and 0.5) following GCTA model (termed as GCTA-*α* model)31. In order to perform a cross-ancestry genetic correlation analysis of cholesterol traits (total cholesterol, HDL cholesterol, and LDL cholesterol), we determined and used optimal *α* based on GCTA models for each trait and ancestry. We did not consider another widely used LDAK-thin model as it will reduce number of common SNPs between ancestry due to LD-pruning.

### Statistical models

#### Univariate Linear Mixed Model

The univariate Linear Mixed Model (LMM) for can be written as, ![Formula][1]</img>  Where **y** is the vector of phenotypic observation, ***b*** is the vector of fixed effects, ***g*** is the vector of additive genetic value and ***e*** is the vector of the residuals. The random effects (***g*** and ***e*** are presumed to be distributed normally with mean zero where **X** and **Z** are incidence matrices

Heritability was estimated using the genetic and residual variances obtained from the univariate LMM, which can be expressed as ![Formula][2]</img>  Here, ![Graphic][3]</img> is the genetic variance and ![Graphic][4]</img> is residual variance. Estimation assumed environmental homogeneity

### Bivariate Linear Mixed Model

The bivariate Linear Mixed Model (LMM) was used to estimate heritability and cross-ancestry genetic correlation using individual level genetic data written as, ![Formula][5]</img>  ![Formula][6]</img>  Where ***y*****1** and ***y*****2** are vector of phenotypic observation, **b****1** and **b****2** are the vector of fixed effects, **g****1** **and g****2** are vector of additive genetic value and **e****1** and **e****2** are the vector of residuals. The random effects (**g****1**, **g****2** and **e****1**, **e****2**) are presumed to be distributed normally with mean zero where **X** and **Z** are incidence matrices i.e. i.e. ![Graphic][7]</img> and ![Graphic][8]</img>.

The variance covariance matrix of observed phenotypes can be written as ![Formula][9]</img>  where, **A** is the genomic relationship matrix (GRM)72-74, which can be estimated based on the genome-wide SNP information, and **I** is an identity matrix which implicitly assumes across individuals of environmental effects and measurement error. The terms, ![Graphic][10]</img> and ![Graphic][11]</img> indicate the genetic and residual variance of the trait for the two-ancestry group, and ![Graphic][12]</img> is the genetic covariances between the two ancestry groups. It is noted that there is no parameter to model residual correlation in **V** because there are no multiple phenotypic measures for any individual, i.e., the phenotypes of the first (second) trait are available only for the first (second) ancestry group.

Cross-ancestry genetic correlation between two random genetic effects can be computed either directly as genetic covariance standardized by the square root of the product of the genetic variances of the two random genetic effects (equation 6) or indirectly by the correlation coefficient of SNP effect sizes38, 75. ![Formula][13]</img>  

### GREML analysis to estimate heritability and cross-ancestry genetic correlations

Bivariate GREML is the cornerstone method to estimate SNP heritability and cross-ancestry genetic correlation using common SNPs across ancestries. The SNPs frequency, heritability model (relationship between heritability and MAF), and the scale factor (*α*) varied across ancestries31. We used a recently proposed approach of estimating GRM31 in combined population, that accounts ancestry specific *α* and ancestry specific allele frequencies for estimating heritability and cross-ancestry genetic correlation. Both estimation of GRM and GREML analysis was implemented in *mtg2*76.

### Genomic prediction

The polygenic score (PGS) is obtained from by aggregating and quantifying single nucleotide polymorphism (SNP) effects. PGS of an individual (*k*) can be defined as cumulative effect of SNP counts with a standard equation as: ![Formula][14]</img>  Here, *β**j* is the SNP effect from discovery GWAS, *m* is the total number of SNPs included in the predictor, *x**jk* is the number of copies (0,1, or 2) of trait associated SNP *j* in the genotype of individual *k*.

## Supporting information

Supplementary information [[supplements/285307_file02.txt]](pending:yes)

## Data Availability

The genotype and phenotype data of the UK Biobank can be accessed through procedures described on its webpage (https://www.ukbiobank.ac.uk/) and summary statistics of BMI and total-, LDL- and HDL-cholesterol from Biobank Japan (BBJ) can be obtained from its website (http://jenger.riken.jp/en/result) MTG2, https://sites.google.com/site/honglee0707/mtg2 PLINK2 version can be downloaded from https://www.cog-genomics.org/plink/ r2redux R-package (https://github.com/mommy003/r2redux from GitHub or from CRAN) The GWAS summary statistics dataset that is generated in this current study and supports the findings have been deposited in the NHGRI-EBI GWAS catalogue with the accession codes GCST90244051, GCST90244052, GCST90244053, GCST90244054, GCST90244055 and GCST90244056; (https://www.ebi.ac.uk/gwas/). GWAS for all SNPs and concordant SNPs for total-, HDL- and LDL-cholesterol can be accessed in following links GWAS of total cholesterol (all SNP) (https://ftp.ebi.ac.uk/pub/databases/gwas/summary\_statistics/GCST90244001-GCST90245000/GCST90244051/) GWAS of total cholesterol (concordant SNP) (https://ftp.ebi.ac.uk/pub/databases/gwas/summary\_statistics/GCST90244001-GCST90245000/GCST90244052/) GWAS of HDL-cholesterol (all SNP) (https://ftp.ebi.ac.uk/pub/databases/gwas/summary\_statistics/GCST90244001-GCST90245000/GCST90244053/) GWAS for HDL-cholesterol (concordant SNP) (https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244054/) GWAS for LDL-cholesterol (all SNP) (https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244055/) GWAS for LDL-cholesterol (concordant SNP) (https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244056/) 

## Web resources and code availability

The genotype and phenotype data of the UK Biobank can be accessed through procedures described on its webpage ([https://www.ukbiobank.ac.uk/](https://www.ukbiobank.ac.uk/)) and summary statistics of BMI and total-, LDL- and HDL-cholesterol from Biobank Japan (BBJ) can be obtained from its website ([http://jenger.riken.jp/en/result](http://jenger.riken.jp/en/result))

MTG2, [https://sites.google.com/site/honglee0707/mtg2](https://sites.google.com/site/honglee0707/mtg2)

PLINK2 version can be downloaded from [https://www.cog-genomics.org/plink/](https://www.cog-genomics.org/plink/)

*r2redux* R-package ([https://github.com/mommy003/r2redux](https://github.com/mommy003/r2redux) from GitHub or from CRAN)

The GWAS summary statistics dataset that is generated in this current study and supports the findings have been deposited in the NHGRI-EBI GWAS catalogue with the accession codes GCST90244051, GCST90244052, GCST90244053, GCST90244054, GCST90244055 and GCST90244056; ([https://www.ebi.ac.uk/gwas/](https://www.ebi.ac.uk/gwas/)). GWAS for all SNPs and concordant SNPs for total-, HDL- and LDL-cholesterol can be accessed in following links

GWAS of total cholesterol (all SNP) ([https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244051/](https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244051/))

GWAS of total cholesterol (concordant SNP) ([https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244052/](https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244052/))

GWAS of HDL-cholesterol (all SNP) ([https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244053/](https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244053/))

GWAS for HDL-cholesterol (concordant SNP) ([https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244054/](https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244054/))

GWAS for LDL-cholesterol (all SNP) ([https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244055/](https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244055/))

GWAS for LDL-cholesterol (concordant SNP) ([https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244056/](https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90244001-GCST90245000/GCST90244056/))

## Declaration of interest

The authors declare that they do not have any competing interests.

## Acknowledgements

This research is supported by the Australian Research Council (DP190100766). We thank the staff and participants of the UK Biobank and Biobank Japan for their important contributions. Our reference number approved by UK Biobank is 14575. The analyses were performed using computational resources provided by the Australian Government through Gadi under the National Computational Merit Allocation Scheme (NCMAS), and HPCs (Statgen server) managed by UniSA IT.

*   Received January 31, 2023.
*   Revision received January 31, 2023.
*   Accepted February 2, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/)

## References

1.  1.Ding, X., et al., The role of cholesterol metabolism in cancer. American journal of cancer research, 2019. 9(2): p. 219.
    
    

2.  2.Yan, S., et al., Bufalin enhances TRAIL-induced apoptosis by redistributing death receptors in lipid rafts in breast cancer cells. Anti-cancer drugs, 2014. 25(6): p. 683–689.
    
    

3.  3.Craig, M.,  S.N.S. Yarrarapu, and  M. Dimri, Biochemistry, cholesterol. 2018.
    
    

4.  4.Musunuru, K. and  S. Kathiresan, Genetics of common, complex coronary artery disease. Cell, 2019. 177(1): p. 132–145.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2019.02.015&link_type=DOI) 

5.  5.WHO, World Health Organization; Cardiovascular diseases (CVDs). 2021.
    
    

6.  6.Trinder, M.,  G.A. Francis, and  L.R. Brunham, Association of monogenic vs polygenic hypercholesterolemia with risk of atherosclerotic cardiovascular disease. JAMA cardiology 2020. 5(4): p. 390–399.
    
    

7.  7.Verbeek, R., et al., Cardiovascular disease risk associated with elevated lipoprotein (a) attenuates at low low-density lipoprotein cholesterol levels in a primary prevention setting. European heart journal, 2018. 39(27): p. 2589–2596.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

8.  8.Go, A.S., et al., Heart disease and stroke statistics—2013 update: a report from the American Heart Association. Circulation, 2013. 127(1): p. e6–e245.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjg6IjEyNy8xL2U2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDIvMDIvMjAyMy4wMS4zMS4yMzI4NTMwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

9.  9.Andaleon, A.,  L.S. Mogil, and  H.E. Wheeler, Gene-based association study for lipid traits in diverse cohorts implicates BACE1 and SIDT2 regulation in triglyceride levels. PearJ 2018. 6: p. e4314.
    
    

10. 10.Trinder, M., et al., Polygenic contribution to low-density lipoprotein cholesterol levels and cardiovascular risk in monogenic familial hypercholesterolemia. Circulation: Genomic Precision Medicine, 2020. 13(5): p. 515–523.
    
    

11. 11.Motazacker, M.M., et al., Evidence of a polygenic origin of extreme high-density lipoprotein cholesterol levels. Arteriosclerosis, Thrombosis and Vascular biology, 2013. 33(7): p. 1521–1528.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYXR2YmFoYSI7czo1OiJyZXNpZCI7czo5OiIzMy83LzE1MjEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMi8wMi8yMDIzLjAxLjMxLjIzMjg1MzA3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

12. 12.Weiss, L.A., et al., The sex-specific genetic architecture of quantitative traits in humans. Nature genetics, 2006. 38(2): p. 218–222.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng1726&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16429159&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234953200018&link_type=ISI) 

13. 13.Ma, L., et al., Genome-wide association analysis of total cholesterol and high-density lipoprotein cholesterol levels using the Framingham heart study data. BMJ medical genetics, 2010. 11(1): p. 1–11.
    
    

14. 14.Klarin, D., et al., Genetics of blood lipids among∼ 300,000 multi-ethnic participants of the Million Veteran Program. Nature Genetics, 2018. 50(11): p. 1514–1523.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0222-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30275531&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

15. 15.Liu, D.J., et al., Exome-wide association study of plasma lipids in> 300,000 individuals. Nature genetics, 2017. 49(12): p. 1758–1766.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3977&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29083408&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

16. 16.Martin, A.R., et al., Clinical use of current polygenic risk scores may exacerbate health disparities. Nature Genetics, 2019. 51(4): p. 584–591.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0379-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30926966&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

17. 17.Bustamante, C.D. and  E.G. Burchard, De la Vega FM. Genomics for the world. Nature, 2011. 475(7355): p. 163–5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/475163a&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21753830&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000292690500024&link_type=ISI) 

18. 18.Oh, S.S., et al., Making precision medicine socially precise. Take a deep breath. American Journal of Respiratory and Critical Care Medicine, 2016.
    
    

19. 19.Duncan, L., et al., Analysis of polygenic risk score usage and performance in diverse human populations. Nature Communications, 2019. 10(1): p. 1–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-13473-y&link_type=DOI) 

20. 20.Morris, A.P., Transethnic meta-analysis of genomewide association studies. Genetic Epidemiology, 2011. 35(8): p. 809–822.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.20630&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22125221&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

21. 21.Okada, Y., et al., Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature, 2014. 506(7488): p. 376–381.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12873&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24390342&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331477800043&link_type=ISI) 

22. 22.Brown, B.C., et al., Transethnic genetic-correlation estimates from summary statistics. The American Journal of Human Genetics, 2016. 99(1): p. 76–88.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2016.05.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27321947&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

23. 23.Galinsky, K.J., et al., Estimating cross-population genetic correlations of causal effect sizes. Genetic Epidemiology, 2019. 43(2): p. 180–188.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

24. 24.Veturi, Y., et al., Modeling heterogeneity in the genetic architecture of ethnically diverse groups using random effect interaction models. Genetics, 2019. 211(4): p. 1395–1407.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6MTA6IjIxMS80LzEzOTUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMi8wMi8yMDIzLjAxLjMxLjIzMjg1MzA3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

25. 25.Takeuchi, F., et al., Interethnic analyses of blood pressure loci in populations of East Asian and European descent. Nature Communications, 2018. 9(1): p. 1–16.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-04053-7&link_type=DOI) 

26. 26.Rosenberg, N.A., et al., Genome-wide association studies in diverse populations. Nature Reviews Genetics, 2010. 11(5): p. 356–366.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg2760&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20395969&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276771400013&link_type=ISI) 

27. 27.McClellan, J. and  M.-C. King, Genetic heterogeneity in human disease. Cell, 2010. 141(2): p. 210–217.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2010.03.032&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20403315&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276738400008&link_type=ISI) 

28. 28.Benyamin, B., et al., Cross-ethnic meta-analysis identifies association of the GPX3-TNIP1 locus with amyotrophic lateral sclerosis. Nature Communications, 2017. 8(1): p. 1–7.
    
    

29. 29.Ding, K. and  I.J. Kullo, Evolutionary genetics of coronary heart disease. Circulation, 2009. 119(3): p. 459–467.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjk6IjExOS8zLzQ1OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAyLzAyLzIwMjMuMDEuMzEuMjMyODUzMDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

30. 30.Shi, H., et al., Population-specific causal disease effect sizes in functionally important regions impacted by selection. Nature communications, 2021. 12(1): p. 1–15.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-24759-5.1&link_type=DOI) 

31. 31.Momin, M.M., et al., A novel method for an unbiased estimate of cross-ancestry genetic correlation using individual-level data. bioRxiv, 2021.
    
    

32. 32.Speed, D., et al., Reevaluation of SNP heritability in complex human traits. Nature Genetics, 2017. 49(7): p. 986–992.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3865&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

33. 33.Martin, A.R., et al., Human demographic history impacts genetic risk prediction across diverse populations. The American Journal of Human Genetics, 2017. 100(4): p. 635–649.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2017.03.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28366442&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

34. 34.Márquez-Luna, C., et al., Multiethnic polygenic risk scores improve risk prediction in diverse populations. Genetic Epidemiology, 2017. 41(8): p. 811–823.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.22083&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

35. 35.Vilhjálmsson, B.J., et al., Modeling linkage disequilibrium increases accuracy of polygenic risk scores. The American Journal of Human Genetics, 2015. 97(4): p. 576–592.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2015.09.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26430803&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

36. 36.Lam, M., et al., Comparative genetic architectures of schizophrenia in East Asian and European populations. Nature Genetics, 2019. 51(12): p. 1670–1678.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/S41588-019-0512-X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

37. 37.Neshat, M., et al., A novel hyper-parameter can increase the prediction accuracy in a single-step genetic evaluation. BioRxiv, 2022: p. 2022.07.03.498620.
    
    

38. 38.Lee, S.H., et al., Estimation of pleiotropy between complex diseases using single-nucleotide polymorphism-derived genomic relationships and restricted maximum likelihood. Bioinformatics, 2012. 28(19): p. 2540–2542.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts474&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22843982&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309687500025&link_type=ISI) 

39. 39.Speed, D., et al., Improved heritability estimation from genome-wide SNPs. The American Journal of Human Genetics, 2012. 91(6): p. 1011–1021.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2012.10.010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23217325&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

40. 40.Speed, D.,  J. Holmes, and  D.J. Balding, Evaluating and improving heritability models using summary statistics. Nature Genetics, 2020. 52(4): p. 458–462.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

41. 41.Iliadou, A., et al., Heritabilities of lipids in young European American and African American twins. Twin Research and Human Genetics, 2005. 8(5): p. 492–498.
    
    

42. 42.Peterson, R.E., et al., Genome-wide association studies in ancestrally diverse populations: opportunities, methods, pitfalls, and recommendations. Cell, 2019. 179(3): p. 589–603.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2019.08.051&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31607513&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

43. 43.Marigorta, U.M. and  A. Navarro, High trans-ethnic replicability of GWAS results implies common causal variants. PLoS Genet, 2013. 9(6): p. e1003566.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1003566&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23785302&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

44. 44.Euesden, J.,  C.M. Lewis, and  P.F. O’Reilly, PRSice: polygenic risk score software. Bioinformatics, 2015. 31(9): p. 1466–1468.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu848&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25550326&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

45. 45.Momin, M.M., et al., Significance tests for R2 of out-of-sample prediction using polygenic scores. The American Journal of Human Genetics, 2023.
    
    

46. 46.Nelson, R.H., Hyperlipidemia as a risk factor for cardiovascular disease. Primary Care: Clinics in Office Practice, 2013. 40(1): p. 195–211.
    
    

47. 47.Tall, A.R., et al., Addressing dyslipidemic risk beyond LDL-cholesterol. The Journal of Clinical Investigation, 2022. 132(1).
    
    

48. 48.Kuchenbaecker, K., et al., The transferability of lipid loci across African, Asian and European cohorts. Nature Communications, 2019. 10(1): p. 1–10.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-13473-y&link_type=DOI) 

49. 49.Cao, C., Analysis of Concordance and Discordance in Genetic Association Studies via Forward-Backward Scoring Scheme, Masters Thesis. 2020, The Ohio State University.
    
    

50. 50.Huang, Q.Q., et al., Transferability of genetic loci and polygenic scores for cardiometabolic traits in British Pakistani and Bangladeshi individuals. Nature communications, 2022. 13(1): p. 1–11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-022-28952-y&link_type=DOI) 

51. 51.Lewis, C.M. and  E. Vassos, Polygenic risk scores: from research tools to clinical instruments. Genomic Medicine, 2020. 12: p. 1–11.
    
    

52. 52.Tropf, F.C., et al., Hidden heritability due to heterogeneity across seven populations. Nature Human Behaviour, 2017. 1(10): p. 757–765.
    
    

53. 53.Bulik-Sullivan, B., et al., ReproGen Consortium Psychiatric Genomics Consortium Genetic Consortium for Anorexia Nervosa of the Wellcome Trust Case Control Consortium 3 An atlas of genetic correlations across human diseases and traits. Nat Genet, 2015. 47(11): p. 1236–1241.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3406&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

54. 54.Gusev, A., et al., Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases. The American Journal of Human Genetics, 2014. 95(5): p. 535–552.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2014.10.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25439723&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

55. 55.Fry, A., et al., Comparison of sociodemographic and health-related characteristics of UK Biobank participants with those of the general population. American Journal of Epidemiology, 2017. 186(9): p. 1026–1034.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwx246&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28641372&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

56. 56.Ollier, W.,  T. Sprosen, and  T. Peakman, UK Biobank: from concept to reality. Future Medicine, 2005.
    
    

57. 57.Novembre, J. and  M. Stephens, Interpreting principal component analyses of spatial population genetic variation. Nature Genetics, 2008. 40(5): p. 646–649.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.139&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18425127&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255366700034&link_type=ISI) 

58. 58.Loh, P.-R., et al., Reference-based phasing using the Haplotype Reference Consortium panel. Nature Genetics, 2016. 48(11): p. 1443.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3679&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27694958&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

59. 59.Bulik-Sullivan, B., et al., ReproGen Consortium Psychiatric Genomics Consortium Genetic Consortium for Anorexia Nervosa of the Wellcome Trust Case Control Consortium 3 An atlas of genetic correlations across human diseases and traits. Nature Genetics, 2015. 47(11): p. 1236–1241.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3406&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

60. 60.Border, R., et al., Imputation of behavioral candidate gene repeat variants in 486,551 publicly-available UK Biobank individuals. European Journal of Human Genetics, 2019. 27(6): p. 963–969.
    
    

61. 61.Lee, S.H.,  W.S.P. Weerasinghe, and  J.H. Van Der Werf, Genotype-environment interaction on human cognitive function conditioned on the status of breastfeeding and maternal smoking around birth. Scientific Reports, 2017. 7(1): p. 1–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-017-14520-8&link_type=DOI) 

62. 62.Peyrot, W.J., et al., Does childhood trauma moderate polygenic risk for depression? A meta-analysis of 5765 subjects from the psychiatric genomics consortium. Biological Psychiatry, 2018. 84(2): p. 138–147.
    
    

63. 63.Purcell, S., et al., PLINK: a tool set for whole-genome association and population-based linkage analyses. The American Journal of Human Genetics, 2007. 81(3): p. 559–575.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/519795&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17701901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

64. 64.Zhou, X.,  H.K. Im, and  S.H. Lee, CORE GREML for estimating covariance between random effects in linear mixed models for complex trait analyses. Nature Communications, 2020. 11(1): p. 1–11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-16113-y&link_type=DOI) 

65. 65.Ni, G., et al., Estimation of genetic correlation via linkage disequilibrium score regression and genomic restricted maximum likelihood. The American Journal of Human Genetics, 2018. 102(6): p. 1185–1194.
    
    

66. 66.Meuleman, W., et al., Index and biological spectrum of human DNase I hypersensitive sites. Nature 2020. 584(7820): p. 244–251.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

67. 67.Yadav, M.L. and  B. Mohapatra, Intergenic regions, also known as spacer DNA. 2018.
    
    

68. 68.Gilbert, W., Genes-in-pieces revisited. Science, 1985. 228: p. 823–825.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czozOiJzY2kiO3M6NToicmVzaWQiO3M6MTI6IjIyOC80NzAxLzgyMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAyLzAyLzIwMjMuMDEuMzEuMjMyODUzMDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

69. 69.Jin, J., et al. Principal components ancestry adjustment for Genetic Analysis Workshop 17 data. in BMC Proceedings. 2011. BioMed Central.
    
    

70. 70.Okbay, A., et al., Genome-wide association study identifies 74 loci associated with educational attainment. Nature, 2016. 533(7604): p. 539–542.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature17671&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27225129&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

71. 71.Yang, J., et al., GCTA: a tool for genome-wide complex trait analysis. The American Journal of Human Genetics, 2011. 88(1): p. 76–82.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2010.11.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21167468&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

72. 72.VanRaden, P.M., Efficient methods to compute genomic predictions. Journal of Dairy Science, 2008. 91(11): p. 4414–4423.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3168/jds.2007-0980&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18946147&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000260277200035&link_type=ISI) 

73. 73.Yang, J., et al., Common SNPs explain a large proportion of the heritability for human height. Nature Genetics, 2010. 42(7): p. 565–569.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.608&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20562875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279242400007&link_type=ISI) 

74. 74.Amin, N.,  C.M. Van Duijn, and  Y.S. Aulchenko, A genomic background based method for association analysis in related individuals. PloS One, 2007. 2(12): p. e1274.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0001274&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18060068&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

75. 75.Bulik-Sullivan, B.K., et al., LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nature genetics, 2015. 47(3): p. 291.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom) 

76. 76.Lee, S.H. and  J.H. Van der Werf, MTG2: an efficient algorithm for multivariate linear mixed model analysis based on genomic information. Bioinformatics, 2016. 32(9): p. 1420–1422.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26755623&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F02%2F02%2F2023.01.31.23285307.atom)

 [1]: /embed/graphic-7.gif
 [2]: /embed/graphic-8.gif
 [3]: /embed/inline-graphic-1.gif
 [4]: /embed/inline-graphic-2.gif
 [5]: /embed/graphic-9.gif
 [6]: /embed/graphic-10.gif
 [7]: /embed/inline-graphic-3.gif
 [8]: /embed/inline-graphic-4.gif
 [9]: /embed/graphic-11.gif
 [10]: /embed/inline-graphic-5.gif
 [11]: /embed/inline-graphic-6.gif
 [12]: /embed/inline-graphic-7.gif
 [13]: /embed/graphic-12.gif
 [14]: /embed/graphic-13.gif