Accounting for linkage disequilibrium in association analysis of diverse populations.

The National Human Genome Research Institute's catalog of published genome-wide association studies (GWAS) lists over 10,000 genetic variants collectively associated with over 800 human diseases or traits. Most of these GWAS have been conducted in European-ancestry populations. Findings gleaned from these studies have led to identification of disease-associated loci and biologic pathways involved in disease etiology. In multiple instances, these genomic findings have led to the development of novel medical therapies or evidence for prescribing a given drug as the appropriate treatment for a given individual beyond phenotypic appearances or socially defined constructs of race or ethnicity. Such findings have implications for populations throughout the globe and GWAS are increasingly being conducted in more diverse populations. A major challenge for investigators seeking to follow up genomic findings between diverse populations is discordant patterns of linkage disequilibrium (LD). We provide an overview of common measures of LD and opportunities for their use in novel methods designed to address challenges associated with following up GWAS conducted in European-ancestry populations in African-ancestry populations or, more generally, between populations with discordant LD patterns. We detail the strengths and weaknesses associated with different approaches. We also describe application of these strategies in follow-up studies of populations with concordant LD patterns (replication) or discordant LD patterns (transferability) as well as fine-mapping studies. We review application of these methods to a variety of traits and diseases.

Evolutionary context for the association of γ-globin, serum uric acid, and hypertension in African Americans.

BACKGROUND: Hyperuricemia and associated cardio-metabolic disorders are more prevalent in African Americans than in European Americans. We used genome-wide admixture mapping and association testing to identify loci with ancestry effects on serum uric acid levels. METHODS: We analyzed 1,976 African Americans from Washington, D.C, including 1,322 individuals from 328 pedigrees and 654 unrelated individuals, enrolled in the Howard University Family Study. We performed admixture mapping and genome-wide association testing using ~800 k autosomal single-nucleotide polymorphisms (SNPs). We performed fine mapping by dense genotyping. We assessed functionality by a combination of bioinformatic annotation, reporter gene assays, and gel shift experiments. We also analyzed 12,641 individuals enrolled in the National Health and Nutrition Examination Survey. RESULTS: We detected a genome-wide significant locus on chromosome 11p15.4 at which serum uric acid levels increased with increasing African ancestry, independent of kidney function. Fine-mapping identified two independent signals in the ß-globin locus. The ancestral allele at SNP rs2855126, located upstream of the hemoglobin, gamma A gene HBG1, was associated with increased serum uric acid levels and higher expression of a reporter gene relative to the derived allele. Hyperuricemia was associated with increased risk of hypertension in 3,767 African Americans (Odds Ratio = 2.48, p = 2.71 × 10(-19)). CONCLUSIONS: Given that increased expression of γ-globin leads to increased levels of fetal hemoglobin which confers protection against malaria, we hypothesize that evolution in Africa of protection against malaria may have occurred at the cost of increased serum uric acid levels, contributing to the high rates of hyperuricemia and associated cardio-metabolic disorders observed in African Americans.

Variants of the adenosine A(2A) receptor gene are protective against proliferative diabetic retinopathy in patients with type 1 diabetes.

AIMS: The adenosine A(2A) receptor (ADORA(2A)) may ameliorate deleterious physiologic effects associated with tissue injury in individuals with diabetes. We explored associations between variants of the ADORA(2A) gene and proliferative diabetic retinopathy (PDR) in a cohort of patients with type 1 diabetes (T1D). METHODS: The participants were from the Pittsburgh Epidemiology of Diabetes Complications prospective study of childhood-onset T1D. Stereoscopic photographs of the retinal fundus taken at baseline, then biennially, for 10 years were used to define PDR according to the modified Airlie House system. Two tagging single nucleotide polymorphisms (tSNPs; rs2236624-C/T and rs4822489-G/T) in the ADORA(2A) gene were selected using the HapMap (haplotype map) reference database. RESULTS: A significant association was observed between SNP rs2236624 and PDR in the recessive genetic model. Participants homozygous for the T allele displayed a decreased risk of developing prevalent PDR (odds ratio, OR = 0.36; p = 0.04) and incident PDR (hazard ratio = 0.156; p = 0.009), and for all cases of PDR combined (OR = 0.23; p = 0.001). The protective effect of T allele homozygosity remained after adjusting for covariates. Similarly, for SNP rs4822489, an association between PDR and T allele homozygosity was observed following covariate adjustment (OR = 0.55; 95% CI: 0.31-0.92; p = 0.04). CONCLUSION: Genetic variants of ADORA(2A) offer statistically significant protection against PDR development in patients with T1D.

Analyses of genome wide association data, cytokines, and gene expression in African-Americans with benign ethnic neutropenia.

BACKGROUND: Benign ethnic neutropenia (BEN) is a hematologic condition associated with people of African ancestry and specific Middle Eastern ethnic groups. Prior genetic association studies in large population showed that rs2814778 in Duffy Antigen Receptor for Chemokines (DARC) gene, specifically DARC null red cell phenotype, was associated with BEN. However, the mechanism of this red cell phenotype leading to low white cell count remained elusive. METHODS: We conducted an extreme phenotype design genome-wide association study (GWAS), analyzed ~16 million single nucleotide polymorphisms (SNP) in 1,178 African-Americans individuals from the Reasons for Geographic and Racial Differences in Stroke (REGARDS) study and replicated from 819 African-American participants in the Atherosclerosis Risk in Communities (ARIC) study. Conditional analyses on rs2814778 were performed to identify additional association signals on chromosome 1q22. In a separate cohort of healthy individuals with and without BEN, whole genome gene expression from peripheral blood neutrophils were analyzed for DARC. RESULTS: We confirmed that rs2814778 in DARC was associated with BEN (p = 4.09×10-53). Conditioning on rs2814778 abolished other significant chromosome 1 associations. Inflammatory cytokines (IL-2, 6, and 10) in participants in the Howard University Family Study (HUFS) and Multi-Ethnic Study in Atherosclerosis (MESA) showed similar levels in individuals homozygous for the rs2814778 allele compared to others, indicating cytokine sink hypothesis played a minor role in leukocyte homeostasis. Gene expression in neutrophils of individuals with and without BEN was also similar except for low DARC expression in BEN, suggesting normal function. BEN neutrophils had slightly activated profiles in leukocyte migration and hematopoietic stem cell mobilization pathways (expression fold change <2). CONCLUSIONS: These results in humans support the notion of DARC null erythroid progenitors preferentially differentiating to myeloid cells, leading to activated DARC null neutrophils egressing from circulation to the spleen, and causing relative neutropenia. Collectively, these human data sufficiently explained the mechanism DARC null red cell phenotype causing BEN and further provided a biologic basis that BEN is clinically benign.

UGT1A1 is a major locus influencing bilirubin levels in African Americans.

Total serum bilirubin is associated with several clinical outcomes, including cardiovascular disease, diabetes and drug metabolism. We conducted a genome-wide association study in 619 healthy unrelated African Americans in an attempt to replicate reported findings in Europeans and Asians and to identify novel loci influencing total serum bilirubin levels. We analyzed a dense panel of over two million genotyped and imputed SNPs in additive genetic models adjusting for age, sex, and the first two significant principal components from the sample covariance matrix of genotypes. Thirty-nine SNPs spanning a 78 kb region within the UGT1A1 displayed P-values <5 × 10(-8). The lowest P-value was 1.7 × 10(-22) for SNP rs887829. None of SNPs in the UGT1A1 remained statistically significant in conditional association analyses that adjusted for rs887829. In addition, SNP rs10929302 located in phenobarbital response enhancer module was significantly associated with bilirubin level with a P-value of 1.37 × 10(-11); this enhancer module is believed to have a critical role in phenobarbital treatment of hyperbilirubinemia. Interestingly, the lead SNP, rs887829, is in strong linkage disequilibrium (LD) (r(2)≥0.74) with rs10929302. Taking advantage of the lower LD and shorter haplotypes in African-ancestry populations, we identified rs887829 as a more refined proxy for the causative variant influencing bilirubin levels. Also, we replicated the reported association between variants in SEMA3C and bilirubin levels. In summary, UGT1A1 is a major locus influencing bilirubin levels and the results of this study promise to contribute to understanding of the etiology and treatment of hyperbilirubinaemia in African-ancestry populations.

The roles of IL-6, IL-10, and IL-1RA in obesity and insulin resistance in African-Americans.

OBJECTIVE: The aim of the study was to investigate the associations between IL-1 receptor antagonist (IL-1RA), IL-6, IL-10, measures of obesity, and insulin resistance in African-Americans. RESEARCH DESIGN AND METHODS: Nondiabetic participants (n = 1025) of the Howard University Family Study were investigated for associations between serum IL (IL-1RA, IL-6, IL-10), measures of obesity, and insulin resistance, with adjustment for age and sex. Measures of obesity included body mass index, waist circumference, hip circumference, waist-to-hip ratio, and percent fat mass. Insulin resistance was assessed using the homeostasis model assessment of insulin resistance (HOMA-IR). Data were analyzed with R statistical software using linear regression and likelihood ratio tests. RESULTS: IL-1RA and IL-6 were associated with measures of obesity and insulin resistance, explaining 4-12.7% of the variance observed (P values < 0.001). IL-1RA was bimodally distributed and therefore was analyzed based on grouping those with low vs. high IL-1RA levels. High IL-1RA explained up to 20 and 12% of the variance in measures of obesity and HOMA-IR, respectively. Among the IL, only high IL-1RA improved the fit of models regressing HOMA-IR on measures of obesity. In contrast, all measures of obesity improved the fit of models regressing HOMA-IR on IL. IL-10 was not associated with obesity measures or HOMA-IR. CONCLUSIONS: High IL-1RA levels and obesity measures are associated with HOMA-IR in this population-based sample of African-Americans. The results suggest that obesity and increased levels of IL-1RA both contribute to the development of insulin resistance.

A genome-wide association study of serum uric acid in African Americans.

BACKGROUND: Uric acid is the primary byproduct of purine metabolism. Hyperuricemia is associated with body mass index (BMI), sex, and multiple complex diseases including gout, hypertension (HTN), renal disease, and type 2 diabetes (T2D). Multiple genome-wide association studies (GWAS) in individuals of European ancestry (EA) have reported associations between serum uric acid levels (SUAL) and specific genomic loci. The purposes of this study were: 1) to replicate major signals reported in EA populations; and 2) to use the weak LD pattern in African ancestry population to better localize (fine-map) reported loci and 3) to explore the identification of novel findings cognizant of the moderate sample size. METHODS: African American (AA) participants (n = 1,017) from the Howard University Family Study were included in this study. Genotyping was performed using the Affymetrix® Genome-wide Human SNP Array 6.0. Imputation was performed using MACH and the HapMap reference panels for CEU and YRI. A total of 2,400,542 single nucleotide polymorphisms (SNPs) were assessed for association with serum uric acid under the additive genetic model with adjustment for age, sex, BMI, glomerular filtration rate, HTN, T2D, and the top two principal components identified in the assessment of admixture and population stratification. RESULTS: Four variants in the gene SLC2A9 achieved genome-wide significance for association with SUAL (p-values ranging from 8.88 × 10(-9) to 1.38 × 10(-9)). Fine-mapping of the SLC2A9 signals identified a 263 kb interval of linkage disequilibrium in the HapMap CEU sample. This interval was reduced to 37 kb in our AA and the HapMap YRI samples. CONCLUSIONS: The most strongly associated locus for SUAL in EA populations was also the most strongly associated locus in this AA sample. This finding provides evidence for the role of SLC2A9 in uric acid metabolism across human populations. Additionally, our findings demonstrate the utility of following-up EA populations GWAS signals in African-ancestry populations with weaker linkage disequilibrium.