Genome-wide Significance Thresholds for Admixture Mapping Studies

Adjusting for principal components can induce spurious associations in genome-wide association studies in admixed populations

Principal component analysis (PCA) is widely used to control for population structure in genome-wide association studies (GWAS). Top principal components (PCs) typically reflect population structure, but challenges arise in deciding how many PCs are needed and ensuring that PCs do not capture other artifacts such as regions with atypical linkage disequilibrium (LD). In response to the latter, many groups suggest performing LD pruning or excluding known high LD regions prior to PCA. However, these suggestions are not universally implemented and the implications for GWAS are not fully understood, especially in the context of admixed populations. In this paper, we investigate the impact of pre-processing and the number of PCs included in GWAS models in African American samples from the Women's Women's Health Initiative SNP Health Association Resource and two Trans-Omics for Precision Medicine Whole Genome Sequencing Project contributing studies (Jackson Heart Study and Genetic Epidemiology of Chronic Obstructive Pulmonary Disease Study). In all three samples, we find the first PC is highly correlated with genome-wide ancestry whereas later PCs often capture local genomic features. The pattern of which, and how many, genetic variants are highly correlated with individual PCs differs from what has been observed in prior studies focused on European populations and leads to distinct downstream consequences: adjusting for such PCs yields biased effect size estimates and elevated rates of spurious associations due to the phenomenon of collider bias. Excluding high LD regions identified in previous studies does not resolve these issues. LD pruning proves more effective, but the optimal choice of thresholds varies across datasets. Altogether, our work highlights unique issues that arise when using PCA to control for ancestral heterogeneity in admixed populations and demonstrates the importance of careful pre-processing and diagnostics to ensure that PCs capturing multiple local genomic features are not included in GWAS models.

Estimating heritability explained by local ancestry and evaluating stratification bias in admixture mapping from summary statistics

The heritability explained by local ancestry markers in an admixed population (hγ2) provides crucial insight into the genetic architecture of a complex disease or trait. Estimation of hγ2 can be susceptible to biases due to population structure in ancestral populations. Here, we present heritability estimation from admixture mapping summary statistics (HAMSTA), an approach that uses summary statistics from admixture mapping to infer heritability explained by local ancestry while adjusting for biases due to ancestral stratification. Through extensive simulations, we demonstrate that HAMSTA hγ2 estimates are approximately unbiased and are robust to ancestral stratification compared to existing approaches. In the presence of ancestral stratification, we show a HAMSTA-derived sampling scheme provides a calibrated family-wise error rate (FWER) of ∼5% for admixture mapping, unlike existing FWER estimation approaches. We apply HAMSTA to 20 quantitative phenotypes of up to 15,988 self-reported African American individuals in the Population Architecture using Genomics and Epidemiology (PAGE) study. We observe hˆγ2 in the 20 phenotypes range from 0.0025 to 0.033 (mean hˆγ2 = 0.012 ± 9.2 × 10-4), which translates to hˆ2 ranging from 0.062 to 0.85 (mean hˆ2 = 0.30 ± 0.023). Across these phenotypes we find little evidence of inflation due to ancestral population stratification in current admixture mapping studies (mean inflation factor of 0.99 ± 0.001). Overall, HAMSTA provides a fast and powerful approach to estimate genome-wide heritability and evaluate biases in test statistics of admixture mapping studies.

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.

Admixture mapping of peripheral artery disease in a Dominican population reveals a putative risk locus on 2q35

Peripheral artery disease (PAD) is a form of atherosclerotic cardiovascular disease, affecting ∼8 million Americans, and is known to have racial and ethnic disparities. PAD has been reported to have a significantly higher prevalence in African Americans (AAs) compared to non-Hispanic European Americans (EAs). Hispanic/Latinos (HLs) have been reported to have lower or similar rates of PAD compared to EAs, despite having a paradoxically high burden of PAD risk factors; however, recent work suggests prevalence may differ between sub-groups. Here, we examined a large cohort of diverse adults in the BioMe biobank in New York City. We observed the prevalence of PAD at 1.7% in EAs vs. 8.5% and 9.4% in AAs and HLs, respectively, and among HL sub-groups, the prevalence was found at 11.4% and 11.5% in Puerto Rican and Dominican populations, respectively. Follow-up analysis that adjusted for common risk factors demonstrated that Dominicans had the highest increased risk for PAD relative to EAs [OR = 3.15 (95% CI 2.33-4.25), p < 6.44 × 10-14]. To investigate whether genetic factors may explain this increased risk, we performed admixture mapping by testing the association between local ancestry and PAD in Dominican BioMe participants (N = 1,813) separately from European, African, and Native American (NAT) continental ancestry tracts. The top association with PAD was an NAT ancestry tract at chromosome 2q35 [OR = 1.96 (SE = 0.16), p < 2.75 × 10-05) with 22.6% vs. 12.9% PAD prevalence in heterozygous NAT tract carriers versus non-carriers, respectively. Fine-mapping at this locus implicated tag SNP rs78529201 located within a long intergenic non-coding RNA (lincRNA) LINC00607, a gene expression regulator of key genes related to thrombosis and extracellular remodeling of endothelial cells, suggesting a putative link of the 2q35 locus to PAD etiology. Efforts to reproduce the signal in other Hispanic cohorts were unsuccessful. In summary, we showed how leveraging health system data helped understand nuances of PAD risk across HL sub-groups and admixture mapping approaches elucidated a putative risk locus in a Dominican population.

Admixture mapping identifies novel Alzheimer's disease risk regions in African Americans

BACKGROUND

This study used admixture mapping to prioritize the genetic regions associated with Alzheimer's disease (AD) in African American (AA) individuals, followed by ancestry-aware regression analysis to fine-map the prioritized regions.

METHODS

We analyzed 10,271 individuals from 17 different AA datasets. We performed admixture mapping and meta-analyzed the results. We then used regression analysis, adjusting for local ancestry main effects and interactions with genotype, to refine the regions identified from admixture mapping. Finally, we leveraged in silico annotation and differential gene expression data to prioritize AD-related variants and genes.

RESULTS

Admixture mapping identified two genome-wide significant loci on chromosomes 17p13.2 (p = 2.2 × 10-5 ) and 18q21.33 (p = 1.2 × 10-5 ). Our fine mapping of the chromosome 17p13.2 and 18q21.33 regions revealed several interesting genes such as the MINK1, KIF1C, and BCL2.

DISCUSSION

Our ancestry-aware regression approach showed that AA individuals have a lower risk of AD if they inherited African ancestry admixture block at the 17p13.2 locus.

HIGHLIGHTS

We identified two genome-wide significant admixture mapping signals: on chromosomes 17p13.2 and 18q21.33, which are novel in African American (AA) populations. Our ancestry-aware regression approach showed that AA individuals have a lower risk of Alzheimer's disease (AD) if they inherited African ancestry admixture block at the 17p13.2 locus. We found that the overall proportion of African ancestry does not differ between the cases and controls that suggest African genetic ancestry alone is not likely to explain the AD prevalence difference between AA and non-Hispanic White populations.

Genome-wide admixture and association analysis identifies African ancestry-specific risk loci of eosinophilic esophagitis in African Americans

BACKGROUND

Eosinophilic esophagitis (EoE), a chronic allergic inflammatory disease, is linked to multiple genetic risk factors, but studies have focused on populations of European ancestry. Few studies have assessed Black or African American (AA) populations for loci involved in EoE susceptibility.

OBJECTIVE

We performed admixture mapping (AM) and genome-wide association study (GWAS) of EoE using participants from AA populations.

METHODS

We conducted AM and GWAS of EoE using 137 EoE cases and 1465 healthy controls from the AA population. Samples were genotyped using molecular evolutionary genetics analysis (MEGA). Genotype imputation was carried out with the Consortium on Asthma Among African-Ancestry Populations in the Americas (CAAPA) reference panel using the Michigan Imputation Server. Global and local ancestry inference was carried out, followed by fine mapping and RNA sequencing. After quality control filtering, over 6,000,000 variants were tested by logistic regression adjusted for sex, age, and global ancestry.

RESULTS

The global African ancestry proportion was found to be significantly lower among cases than controls (0.751 vs 0.786, P = .012). Case-only AM identified 3 significant loci (9p13.3, 12q24.22-23, and 15q11.2) associated with EoE, of which 12q24.22-23 and 9p13.3 were further replicated in the case-control analysis, with associations observed with African ancestry. Fine mapping and multiomic functional annotations prioritized the variants rs11068264 (FBXW8) and rs7307331 (VSIG10) at 12q24.23 and rs2297879 (ARHGEF39) at 9p13.3. GWAS identified 1 genome-wide significant locus at chromosome 1p22.3 (rs17131726, DDAH1) and 10 other suggestive loci. Most GWAS variants were low-frequency African ancestry-specific variants. RNA sequencing revealed that esophageal DDAH1 and VSIG10 were downregulated and ARHGEF39 upregulated among EoE cases.

CONCLUSIONS

GWAS and AM for EoE in AA revealed that African ancestry-specific genetic susceptibility loci exist at 1p22.3, 9p13.3, and 12q24.23, providing evidence of ancestry-specific inheritance of EoE. More independent genetic studies of different ancestries for EoE are needed.

Estimating heritability explained by local ancestry and evaluating stratification bias in admixture mapping from summary statistics

The heritability explained by local ancestry markers in an admixed population h γ 2 provides crucial insight into the genetic architecture of a complex disease or trait. Estimation of h γ 2 can be susceptible to biases due to population structure in ancestral populations. Here, we present a novel approach, Heritability estimation from Admixture Mapping Summary STAtistics (HAMSTA), which uses summary statistics from admixture mapping to infer heritability explained by local ancestry while adjusting for biases due to ancestral stratification. Through extensive simulations, we demonstrate that HAMSTA h γ 2 estimates are approximately unbiased and are robust to ancestral stratification compared to existing approaches. In the presence of ancestral stratification, we show a HAMSTA-derived sampling scheme provides a calibrated family-wise error rate (FWER) of ~5% for admixture mapping, unlike existing FWER estimation approaches. We apply HAMSTA to 20 quantitative phenotypes of up to 15,988 self-reported African American individuals in the Population Architecture using Genomics and Epidemiology (PAGE) study. We observe h ˆ γ 2 in the 20 phenotypes range from 0.0025 to 0.033 (mean h ˆ γ 2 = 0.012 + / - 9.2 × 10 - 4 ), which translates to h ˆ 2 ranging from 0.062 to 0.85 (mean h ˆ 2 = 0.30 + / - 0.023 ). Across these phenotypes we find little evidence of inflation due to ancestral population stratification in current admixture mapping studies (mean inflation factor of 0.99 +/- 0.001). Overall, HAMSTA provides a fast and powerful approach to estimate genome-wide heritability and evaluate biases in test statistics of admixture mapping studies.

Admixture Mapping of Peripheral Artery Disease in a Dominican Population Reveals a Novel Risk Locus on 2q35

Peripheral artery disease (PAD) is a form of atherosclerotic cardiovascular disease, affecting ∼8 million Americans, and is known to have racial and ethnic disparities. PAD has been reported to have significantly higher prevalence in African Americans (AAs) compared to non-Hispanic European Americans (EAs). Hispanic/Latinos (HLs) have been reported to have lower or similar rates of PAD compared to EAs, despite having a paradoxically high burden of PAD risk factors, however recent work suggests prevalence may differ between sub-groups. Here we examined a large cohort of diverse adults in the Bio Me biobank in New York City (NYC). We observed the prevalence of PAD at 1.7% in EAs vs 8.5% and 9.4% in AAs and HLs, respectively; and among HL sub-groups, at 11.4% and 11.5% in Puerto Rican and Dominican populations, respectively. Follow-up analysis that adjusted for common risk factors demonstrated that Dominicans had the highest increased risk for PAD relative to EAs (OR=3.15 (95% CI 2.33-4.25), P <6.44×10 -14 ). To investigate whether genetic factors may explain this increased risk, we performed admixture mapping by testing the association between local ancestry (LA) and PAD in Dominican Bio Me participants (N=1,940) separately for European (EUR), African (AFR) and Native American (NAT) continental ancestry tracts. We identified a NAT ancestry tract at chromosome 2q35 that was significantly associated with PAD (OR=2.05 (95% CI 1.51-2.78), P <4.06×10 -6 ) with 22.5% vs 12.5% PAD prevalence in heterozygous NAT tract carriers versus non-carriers, respectively. Fine-mapping at this locus implicated tag SNP rs78529201 located within a long intergenic non-coding RNA (lincRNA) LINC00607 , a gene expression regulator of key genes related to thrombosis and extracellular remodeling of endothelial cells, suggesting a putative link of the 2q35 locus to PAD etiology. In summary, we showed how leveraging health systems data helped understand nuances of PAD risk across HL sub-groups and admixture mapping approaches elucidated a novel risk locus in a Dominican population.

Genetic ancestry, differential gene expression, and survival in pediatric B-cell acute lymphoblastic leukemia

BACKGROUND

Black children have lower incidence yet worse survival than White and Latinx children with B-cell acute lymphoblastic leukemia (B-ALL). It is unclear how reported race/ethnicity (RRE) is associated with death in B-ALL after accounting for differentially expressed genes associated with genetic ancestry.

METHODS

Using Phase 1 and 2 NCI TARGET B-ALL cases (N = 273; RRE-Black = 21, RRE-White = 162, RRE-Latinx = 69, RRE-Other = 9, RRE-Unknown = 12), we estimated proportions of African (AFR), European (EUR), and Amerindian (AMR) genetic ancestry. We estimated hazard ratios (HR) and 95% confidence intervals (95% CI) between ancestry and death while adjusting for RRE and clinical measures. We identified genes associated with genetic ancestry and adjusted for them in RRE and death associations.

RESULTS

Genetic ancestry varied within RRE (RRE-Black, AFR proportion: Mean: 78.5%, Range: 38.2%-93.6%; RRE-White, EUR proportion: Mean: 94%, Range: 1.6%-99.9%; RRE-Latinx, AMR proportion: Mean: 52.0%, Range: 1.2%-98.7%). We identified 10, 1, and 6 differentially expressed genes (p_adjusted  <0.05) associated with AFR, AMR, and EUR ancestry proportion, respectively. We found AMR and AFR ancestry were statistically significantly associated with death (AMR each 10% HR: 1.05, 95% CI: 1.03-1.17, AFR each 10% increase HR: 1.03, 95% CI:1.01-1.19). RRE differences in the risk of death were larger in magnitude upon adjustment for genes associated with genetic ancestry for RRE-Black, but not RRE-Latinx children (RRE-Black HR: 3.35, 95% CI: 1.31, 8.53; RRE-Latinx HR: 1.47, 0.88-2.45).

CONCLUSIONS

Our work highlights B-ALL survival differences by RRE after adjusting for ancestry differentially expressed genes suggesting other factors impacting survival are important.

Localized variation in ancestral admixture identifies pilocytic astrocytoma risk loci among Latino children

Background

Pilocytic astrocytoma (PA) is the most common pediatric brain tumor. PA has at least a 50% higher incidence in populations of European ancestry compared to other ancestral groups, which may be due in part to genetic differences.

Methods

We first compared the global proportions of European, African, and Amerindian ancestries in 301 PA cases and 1185 controls of self-identified Latino ethnicity from the California Biobank. We then conducted admixture mapping analysis to assess PA risk with local ancestry.

Results

We found PA cases had a significantly higher proportion of global European ancestry than controls (case median = 0.55, control median = 0.51, P value = 3.5x10^-3). Admixture mapping identified 13 SNPs in the 6q14.3 region (SNX14) contributing to risk, as well as three other peaks approaching significance on chromosomes 7, 10 and 13. Downstream fine mapping in these regions revealed several SNPs potentially contributing to childhood PA risk.

Conclusions

There is a significant difference in genomic ancestry associated with Latino PA risk and several genomic loci potentially mediating this risk.

Childhood brain tumors are among the most prevalent and lethal childhood cancers. Despite this, the epidemiology as well as genetic risks are not well defined. For example, children of European ancestry have a higher risk of contracting pilocytic astrocytoma (PA) compared to other ancestries, but the genetic or environmental basis for this is unknown. Latino children are a mixture of multiple ancestries including European, African, and Native American. Using a group of Californian Latino children, we show that the risk of PA increases when a Latino child has a higher proportion of European ancestry. This global ancestry difference shows that germline genetic risk alleles contribute to a higher PA risk in children of European descendent. Moreover, this ancestral risk is localized to specific regions of the genome, especially in Chromosome 6 near the SNX14 gene, which is associated with cancer-related growth signaling pathway described by MAPK/ERK. This result brings us one step closer to understanding the etiology of this common childhood brain tumor.

Genome-Wide Admixture Mapping of Estimated Glomerular Filtration Rate and Chronic Kidney Disease Identifies European and African Ancestry-of-Origin Loci in Hispanic and Latino Individuals in the United States

BACKGROUND

Admixture mapping is a powerful approach for gene mapping of complex traits that leverages the diverse genetic ancestry in populations with recent admixture, such as Hispanic or Latino individuals in the United States. These individuals have an increased risk of CKD.

METHODS

We performed genome-wide admixture mapping for both CKD and eGFR in a sample of 12,601 participants from the Hispanic Community Health Study/Study of Latinos, with validation in a sample of 8191 Black participants from the Women's Health Initiative (WHI). We also compared the findings with those from a conventional genome-wide association study.

RESULTS

Three novel ancestry-of-origin loci were identified on chromosomes 2, 14, and 15 for CKD and eGFR. The chromosome 2 locus comprises two European ancestry regions encompassing the FSHR and NRXN1 genes, with European ancestry at this locus associated with increased CKD risk. The chromosome 14 locus, found within the DLK1-DIO3 imprinted domain, was associated with lower eGFR and driven by European ancestry. The eGFR-associated locus on chromosome 15 included intronic variants of RYR3 and was within an African-specific genomic region associated with higher eGFR. The genome-wide association study failed to identify significant associations in these regions. We validated the chromosome 14 and 15 loci associated with eGFR in the WHI Black participants.

CONCLUSIONS

This study provides evidence of shared ancestry-specific genomic regions influencing eGFR in Hispanic or Latino individuals and Black individuals and illustrates the potential for leveraging genetic ancestry in recently admixed populations for the discovery of novel candidate loci for kidney phenotypes.

Local Ancestry Adjusted Allelic Association Analysis Robustly Captures Tuberculosis Susceptibility Loci

Pulmonary tuberculosis (TB), caused by Mycobacterium tuberculosis, is a complex disease. The risk of developing active TB is in part determined by host genetic factors. Most genetic studies investigating TB susceptibility fail to replicate association signals particularly across diverse populations. South African populations arose because of multi-wave genetic admixture from the indigenous KhoeSan, Bantu-speaking Africans, Europeans, Southeast Asian-and East Asian populations. This has led to complex genetic admixture with heterogenous patterns of linkage disequilibrium and associated traits. As a result, precise estimation of both global and local ancestry is required to prevent both false positive and false-negative associations. Here, 820 individuals from South Africa were genotyped on the SNP-dense Illumina Multi-Ethnic Genotyping Array (∼1.7M SNPs) followed by local and global ancestry inference using RFMix. Local ancestry adjusted allelic association (LAAA) models were utilized owing to the extensive genetic heterogeneity present in this population. Hence, an interaction term, comprising the identification of the minor allele that corresponds to the ancestry present at the specific locus under investigation, was included as a covariate. One SNP (rs28647531) located on chromosome 4q22 was significantly associated with TB susceptibility and displayed a SNP minor allelic effect (G allele, frequency = 0.204) whilst correcting for local ancestry for Bantu-speaking African ancestry (p-value = 5.518 × 10^-7; OR = 3.065; SE = 0.224). Although no other variants passed the significant threshold, clear differences were observed between the lead variants identified for each ancestry. Furthermore, the LAAA model robustly captured the source of association signals in multi-way admixed individuals from South Africa and allowed the identification of ancestry-specific disease risk alleles associated with TB susceptibility that have previously been missed.

Admixture mapping analysis reveals differential genetic ancestry associated with Chagas disease susceptibility in the Colombian population

Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP) and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the Major Histocompatibility Complex region (OR = 0.74, 95%CI = 0.66-0.83, lowest p-value = 4.53x10-8). The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95%CI = 0.90-0.97, p-value = 3.54x10-4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the Human Leukocyte Antigen (HLA) genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils, and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease.

Admixture mapping reveals the association between Native American ancestry at 3q13.11 and reduced risk of Alzheimer's disease in Caribbean Hispanics

BACKGROUND

Genetic studies have primarily been conducted in European ancestry populations, identifying dozens of loci associated with late-onset Alzheimer's disease (AD). However, much of AD's heritability remains unexplained; as the prevalence of AD varies across populations, the genetic architecture of the disease may also vary by population with the presence of novel variants or loci.

METHODS

We conducted genome-wide analyses of AD in a sample of 2565 Caribbean Hispanics to better understand the genetic contribution to AD in this population. Statistical analysis included both admixture mapping and association testing. Evidence for differential gene expression within regions of interest was collected from independent transcriptomic studies comparing AD cases and controls in samples with primarily European ancestry.

RESULTS

Our genome-wide association study of AD identified no loci reaching genome-wide significance. However, a genome-wide admixture mapping analysis that tests for association between a haplotype's ancestral origin and AD status detected a genome-wide significant association with chromosome 3q13.11 (103.7-107.7Mb, P = 8.76E-07), driven by a protective effect conferred by the Native American ancestry (OR = 0.58, 95%CI = 0.47-0.73). Within this region, two variants were significantly associated with AD after accounting for the number of independent tests (rs12494162, P = 2.33E-06; rs1731642, P = 6.36E-05). The significant admixture mapping signal is composed of 15 haplotype blocks spanning 5 protein-coding genes (ALCAM, BBX, CBLB, CCDC54, CD47) and four brain-derived topologically associated domains, and includes markers significantly associated with the expression of ALCAM, BBX, CBLB, and CD47 in the brain. ALCAM and BBX were also significantly differentially expressed in the brain between AD cases and controls with European ancestry.

CONCLUSION

These results provide multiethnic evidence for a relationship between AD and multiple genes at 3q13.11 and illustrate the utility of leveraging genetic ancestry diversity via admixture mapping for new insights into AD.

The impact of global and local Polynesian genetic ancestry on complex traits in Native Hawaiians

Epidemiological studies of obesity, Type-2 diabetes (T2D), cardiovascular diseases and several common cancers have revealed an increased risk in Native Hawaiians compared to European- or Asian-Americans living in the Hawaiian islands. However, there remains a gap in our understanding of the genetic factors that affect the health of Native Hawaiians. To fill this gap, we studied the genetic risk factors at both the chromosomal and sub-chromosomal scales using genome-wide SNP array data on ~4,000 Native Hawaiians from the Multiethnic Cohort. We estimated the genomic proportion of Native Hawaiian ancestry ("global ancestry," which we presumed to be Polynesian in origin), as well as this ancestral component along each chromosome ("local ancestry") and tested their respective association with binary and quantitative cardiometabolic traits. After attempting to adjust for non-genetic covariates evaluated through questionnaires, we found that per 10% increase in global Polynesian genetic ancestry, there is a respective 8.6%, and 11.0% increase in the odds of being diabetic (P = 1.65×10-4) and having heart failure (P = 2.18×10-4), as well as a 0.059 s.d. increase in BMI (P = 1.04×10-10). When testing the association of local Polynesian ancestry with risk of disease or biomarkers, we identified a chr6 region associated with T2D. This association was driven by an uniquely prevalent variant in Polynesian ancestry individuals. However, we could not replicate this finding in an independent Polynesian cohort from Samoa due to the small sample size of the replication cohort. In conclusion, we showed that Polynesian ancestry, which likely capture both genetic and lifestyle risk factors, is associated with an increased risk of obesity, Type-2 diabetes, and heart failure, and that larger cohorts of Polynesian ancestry individuals will be needed to replicate the putative association on chr6 with T2D.

Whole genome sequence analyses of eGFR in 23,732 people representing multiple ancestries in the NHLBI trans-omics for precision medicine (TOPMed) consortium

Background

Genetic factors that influence kidney traits have been understudied for low frequency and ancestry-specific variants.

Methods

We combined whole genome sequencing (WGS) data from 23,732 participants from 10 NHLBI Trans-Omics for Precision Medicine (TOPMed) Program multi-ethnic studies to identify novel loci for estimated glomerular filtration rate (eGFR). Participants included European, African, East Asian, and Hispanic ancestries. We applied linear mixed models using a genetic relationship matrix estimated from the WGS data and adjusted for age, sex, study, and ethnicity.

Findings

When testing single variants, we identified three novel loci driven by low frequency variants more commonly observed in non-European ancestry (PRKAA2, rs180996919, minor allele frequency [MAF] 0.04%, P = 6.1 × 10⁻¹¹; METTL8, rs116951054, MAF 0.09%, P = 4.5 × 10⁻⁹; and MATK, rs539182790, MAF 0.05%, P = 3.4 × 10⁻⁹). We also replicated two known loci for common variants (rs2461702, MAF=0.49, P = 1.2 × 10⁻⁹, nearest gene GATM, and rs71147340, MAF=0.34, P = 3.3 × 10⁻⁹, CDK12). Testing aggregated variants within a gene identified the MAF gene. A statistical approach based on local ancestry helped to identify replication samples for ancestry-specific variants.

Interpretation

This study highlights challenges in studying variants influencing kidney traits that are low frequency in populations and more common in non-European ancestry.

Lack of Association of IL6 polymorphism with the susceptibility to Chagas disease in Latin American populations

The aim of the present study was to analyze IL6 rs1800795 genetic variant in the susceptibility to Trypanosoma cruzi infection and in the development of chronic Chagas cardiomyopathy (CCC), in five independent Latin American cohorts. A total of 3,087 individuals from Latin American countries (Argentina, Bolivia, Peru, and two cohorts from Colombia) were studied. In all cohorts, patients were classified as seropositive for T. cruzi antigens (n= 1,963) and seronegative (n= 1,124). Based on clinical evaluation, the seropositive patients, were classified as CCC (n= 900) and asymptomatic (n= 1,063). No statistically significant differences in the frequency of IL6 rs1800795 between seropositive and seronegative, or between CCC and asymptomatic patients, were found. Furthermore, after the meta-analysis no statistically significant differences were observed. Our results do not support a contribution of IL6 rs1800795 genetic variant in the susceptibility to the infection and the development of chronic Chagas cardiomyopathy in the studied populations.

Selection and hybridization shaped the rapid spread of African honey bee ancestry in the Americas

Recent biological invasions offer 'natural' laboratories to understand the genetics and ecology of adaptation, hybridization, and range limits. One of the most impressive and well-documented biological invasions of the 20th century began in 1957 when Apis mellifera scutellata honey bees swarmed out of managed experimental colonies in Brazil. This newly-imported subspecies, native to southern and eastern Africa, both hybridized with and out-competed previously-introduced European honey bee subspecies. Populations of scutellata-European hybrid honey bees rapidly expanded and spread across much of the Americas in less than 50 years. We use broad geographic sampling and whole genome sequencing of over 300 bees to map the distribution of scutellata ancestry where the northern and southern invasions have presently stalled, forming replicated hybrid zones with European bee populations in California and Argentina. California is much farther from Brazil, yet these hybrid zones occur at very similar latitudes, consistent with the invasion having reached a climate barrier. At these range limits, we observe genome-wide clines for scutellata ancestry, and parallel clines for wing length that span hundreds of kilometers, supporting a smooth transition from climates favoring scutellata-European hybrid bees to climates where they cannot survive winter. We find no large effect loci maintaining exceptionally steep ancestry transitions. Instead, we find most individual loci have concordant ancestry clines across South America, with a build-up of somewhat steeper clines in regions of the genome with low recombination rates, consistent with many loci of small effect contributing to climate-associated fitness trade-offs. Additionally, we find no substantial reductions in genetic diversity associated with rapid expansions nor complete dropout of scutellata ancestry at any individual loci on either continent, which suggests that the competitive fitness advantage of scutellata ancestry at lower latitudes has a polygenic basis and that scutellata-European hybrid bees maintained large population sizes during their invasion. To test for parallel selection across continents, we develop a null model that accounts for drift in ancestry frequencies during the rapid expansion. We identify several peaks within a larger genomic region where selection has pushed scutellata ancestry to high frequency hundreds of kilometers past the present cline centers in both North and South America and that may underlie high-fitness traits driving the invasion.

Coagulation factor VIII: Relationship to cardiovascular disease risk and whole genome sequence and epigenome-wide analysis in African Americans

BACKGROUND

Prospective studies have suggested higher factor VIII (FVIII) levels are an independent risk factor for coronary heart disease (CHD) and stroke. However, limited information, including on genetic and epigenetic contributors to FVIII variation, is available specifically among African Americans (AAs), who have higher FVIII levels than Europeans.

OBJECTIVES

We measured FVIII levels in ~3400 AAs from the community-based Jackson Heart Study and assessed genetic, epigenetic, and epidemiological correlates of FVIII, as well as incident cardiovascular disease (CVD) associations.

METHODS

We assessed cross-sectional associations of FVIII with CVD risk factors as well as incident CHD, stroke, heart failure, and mortality associations. We additionally assessed associations with TOPMed whole genome sequencing data and an epigenome-wide methylation array.

RESULTS

Our results confirmed associations between FVIII and risk of incident CHD events and total mortality in AAs; mortality associations were largely independent of traditional risk factors. We also demonstrate an association of FVIII with incident heart failure, independent of B-type natriuretic peptide. Two genomic regions were strongly associated with FVIII (ABO and VWF). The index variant at VWF is specific to individuals of African descent and is distinct from the previously reported European VWF association signal. Epigenome-wide association analysis showed significant FVIII associations with several CpG sites in the ABO region. However, after adjusting for ABO genetic variants, ABO CpG sites were not significant.

CONCLUSIONS

Larger sample sizes of AAs will be required to discover additional genetic and epigenetic contributors to FVIII phenotypic variation, which may have consequences for CVD health disparities.

Association of IL18 genetic polymorphisms with Chagas disease in Latin American populations

Host genetic factors have been suggested to play an important role in the susceptibility to Chagas disease. Given the influence of interleukin 18 (IL-18) in the development of the disease, in the present study, we analyzed three IL18 genetic variants (rs2043055, rs1946518, rs360719) regarding the predisposition to Trypanosoma cruzi infection and the development of chronic Chagas cardiomyopathy (CCC), in different Latin America populations. Genetic data of 3,608 patients from Colombia, Bolivia, Argentina, and Brazil were meta-analyzed to validate previous findings with increased statistical power. Seropositive and seronegative individuals were compared for T. cruzi infection susceptibility. In the Colombian cohort, the allelic frequencies of the three variants showed a significant association, with adjustment for sex and age, and also after applying multiple testing adjustments. Among the Colombian and Argentinean cohorts, rs360719 showed a significant genetic effect in a fixed-effects meta-analysis after a Bonferroni correction (OR: 0.76, CI: 0.66-0.89, P = 0.001). For CCC, the rs2043055 showed an association with protection from cardiomyopathy in the Colombian cohort (OR: 0.79, CI: 0.64-0.99, P = 0.037), with adjustment for sex and age, and after applying multiple testing adjustments. The meta-analysis of the CCC vs. asymptomatic patients from the four cohorts showed no evidence of association. In conclusion, our results validated the association found previously in the Colombian cohort suggesting that IL18 rs360719 plays an important role in the susceptibility to T. cruzi infection and no evidence of association was found between the IL18 genetic variants and CCC in the Latin American population studied.