Although a considerable proportion of serum lipids loci identified in European ancestry individuals (EA) replicate in African Americans (AA), interethnic differences in the distribution of serum lipids suggest that some genetic determinants differ by ethnicity. We conducted a comprehensive evaluation of five lipid candidate genes to identify variants with ethnicity-specific effects. We sequenced ABCA1, LCAT, LPL, PON1, and SERPINE1 in 48 AA individuals with extreme serum lipid concentrations (high HDLC/low TG or low HDLC/high TG). Identified variants were genotyped in the full population-based sample of AA (n = 1694) and tested for an association with serum lipids. rs328 (LPL) and correlated variants were associated with higher HDLC and lower TG. Interestingly, a stronger effect was observed on a “European” vs. “African” genetic background at this locus. To investigate this effect, we evaluated the region among West Africans (WA). For TG, the effect size among WA was the same in AA with only African local ancestry (2–3% lower TG), while the larger association among AA with local European ancestry matched previous reports in EA (10%). For HDLC, there was no association with rs328 in AA with only African local ancestry or in WA, while the association among AA with European local ancestry was much greater than what has been observed for EA (15 vs. ∼5 mg/dl), suggesting an interaction with an environmental or genetic factor that differs by ethnicity. Beyond this ancestry effect, the importance of African ancestry-focused, sequence-based work was also highlighted by serum lipid associations of variants that were in higher frequency (or present only) among those of African ancestry. By beginning our study with the sequence variation present in AA individuals, investigating local ancestry effects, and seeking replication in WA, we were able to comprehensively evaluate the role of a set of candidate genes in serum lipids in AA.
Most of the work on the genetic epidemiology of serum lipids in African Americans (AA) has focused on replicating findings that were identified in European ancestry individuals. While this can be very informative about the generalizability of lipids loci across populations, African ancestry-specific variation will be missed using this approach. Our aim was to comprehensively evaluate five lipid candidate genes in an AA population, from the identification of variants of interest to population-level analysis of high-density lipoprotein cholesterol (HDLC) and triglycerides (TG). We sequenced five genes in individuals with extreme lipids (n = 48) drawn from a population-based study of AA. The variants identified were genotyped in 1,694 AA and analyzed. Notable among the findings were the observation of ancestry specific effect for several variants in the LPL gene among these admixed individuals, with a greater effect observed among those with European ancestry in this region. These associations were further elucidated by replication in West Africans. By beginning with the sequence variation present among AA, investigating ancestry effects, and seeking replication in West Africans, we were able to comprehensively evaluate these candidate genes with a focus on African ancestry individuals.
Research in industrialized countries has demonstrated that a key factor limiting the control of hypertension is poor patient adherence and that the most successful interventions for long-term adherence employ multiple strategies. Very little data exist on this question in low-income countries, where medication-taking behavior may be less well developed. We conducted a treatment adherence trial of 544 subjects (mean age ~63 years, mean BP ~168/92 mmHg) with previously untreated hypertension in urban and rural Nigeria. Eligible participants were randomized to one of two arms: clinic management only, or clinic management + home visits. Both interventions included three elements: a community based, nurse-led treatment program with physician backup; facilitation of clinic visits and health education; and the use of diuretics plus a beta blocker as needed. After initial diagnosis, the management protocol was implemented by a nurse with physician backup. Participants were evaluated monthly for 6 months. Medication adherence was assessed with pill count and urine testing. Dropout rates, by treatment group, ranged from 12% to 28%. Among participants who completed the 6 month trial, overall adherence was high (~77% of participants took > 98% of prescribed pills). Adherence did not differ by treatment arm, but was better at the rural than the urban site and among those with higher baseline BP. Hypertension control (BP < 140/90 mmHg) was achieved in ~66% of participants at 6 months. This community-based intervention confirms relatively modest default rates compared to industrialized societies, and suggests that medication adherence can be high in developing world settings in clinic attenders.
anti-hypertensive; adherence; compliance trial
Although an increasing number of hypertension-associated genetic variants is being reported, replication of these findings in independent studies has been challenging. Several genes in a human chromosome 1q linkage region have been reported to be associated with hypertension. We examined polymorphisms in three of these genes (ATP1B1, RGS5 and SELE) in relation to hypertension and blood pressure in a cohort of African–Americans.
We genotyped 87 single nucleotide polymorphisms (SNPs) from the ATP1B1, RGS5 and SELE genes in a well characterized cohort of 968 African–Americans and performed a case–control study to identify susceptibility alleles for hypertension and blood pressure regulation. Single SNP and haplotype association testing was done under an additive genetic model with adjustment for age, sex, BMI and ancestry-by-genotype (principal components).
A total of 12 SNPs showed nominal association with hypertension and/or blood pressure. The strongest signal for hypertension was for rs2815272 in the RGS5 gene (P = 9.3 × 10–3). For SBP, rs3917420 in the SELE gene (P = 9.0 × 10–4) and rs4657251 in the RGS5 gene (P = 9.7 × 10–3) were the top hits. Effect size for each of these variants was approximately 2–3 mmHg. A five-SNP haplotype in the SELE gene also showed significant association with SBP after correction for multiple testing (P < 0.01).
These findings provide additional support for the genetic role of ATP1B1, RGS5 and SELE in hypertension and blood pressure regulation.
African–Americans; candidate gene; haplotype; hypertension; single nucleotide polymorphism
Genome-wide association studies (GWAS) have identified 36 loci associated with body mass index (BMI), predominantly in populations of European ancestry. We conducted a meta-analysis to examine the association of >3.2 million SNPs with BMI in 39,144 men and women of African ancestry, and followed up the most significant associations in an additional 32,268 individuals of African ancestry. We identified one novel locus at 5q33 (GALNT10, rs7708584, p=3.4×10−11) and another at 7p15 when combined with data from the Giant consortium (MIR148A/NFE2L3, rs10261878, p=1.2×10−10). We also found suggestive evidence of an association at a third locus at 6q16 in the African ancestry sample (KLHL32, rs974417, p=6.9×10−8). Thirty-two of the 36 previously established BMI variants displayed directionally consistent effect estimates in our GWAS (binomial p=9.7×10−7), of which five reached genome-wide significance. These findings provide strong support for shared BMI loci across populations as well as for the utility of studying ancestrally diverse populations.
Insulin resistance (IR) is a key determinant of type 2 diabetes (T2D) and other metabolic disorders. This genome-wide association study (GWAS) was designed to shed light on the genetic basis of fasting insulin (FI) and IR in 927 non-diabetic African Americans. 5 396 838 single-nucleotide polymorphisms (SNPs) were tested for associations with FI or IR with adjustments for age, sex, body mass index, hypertension status and first two principal components. Genotyped SNPs (n = 12) with P < 5 × 10−6 in African Americans were carried forward for de novo genotyping in 570 non-diabetic West Africans. We replicated SNPs in or near SC4MOL and TCERG1L in West Africans. The meta-analysis of 1497 African Americans and West Africans yielded genome-wide significant associations for SNPs in the SC4MOL gene: rs17046216 (P = 1.7 × 10−8 and 2.9 × 10−8 for FI and IR, respectively); and near the TCERG1L gene with rs7077836 as the top scoring (P = 7.5 × 10−9 and 4.9 × 10−10 for FI and IR, respectively). In silico replication in the MAGIC study (n = 37 037) showed weak but significant association (adjusted P-value of 0.0097) for rs34602777 in the MYO5A gene. In addition, we replicated previous GWAS findings for IR and FI in Europeans for GCKR, and for variants in four T2D loci (FTO, IRS1, KLF14 and PPARG) which exert their action via IR. In summary, variants in/near SC4MOL, and TCERG1L were associated with FI and IR in this cohort of African Americans and were replicated in West Africans. SC4MOL is under-expressed in an animal model of T2D and plays a key role in lipid biosynthesis, with implications for the regulation of energy metabolism, obesity and dyslipidemia. TCERG1L is associated with plasma adiponectin, a key modulator of obesity, inflammation, IR and diabetes.
Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10−6 were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10−8 for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10−8 for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10−8; RREB1: p = 5.7×10−8). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity.
Central obesity is a marker of body fat distribution and is known to have a genetic underpinning. Few studies have reported genome-wide association study (GWAS) results among individuals of predominantly African ancestry (AA). We performed a collaborative meta-analysis in order to identify genetic loci associated with body fat distribution in AA individuals using waist circumference (WC) and waist to hip ratio (WHR) as measures of fat distribution, with and without adjustment for body mass index (BMI). We uncovered 2 genetic loci potentially associated with fat distribution: LHX2 in association with WC-adjusted-for-BMI and at RREB1 for WHR-adjusted-for-BMI. Six of fourteen previously reported loci for waist in EA populations were significant in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). These findings reinforce the concept that there are loci for body fat distribution that are independent of generalized adiposity.
C-reactive protein (CRP) largely has been studied in white non-Hispanic cohorts. There is limited information on CRP’s range of values, heritability and relation to cardiovascular disease (CVD) risk factors in African Americans. We sought to evaluate the distribution, clinical correlates, heritability and genetic linkage of log-transformed CRP in participants of the middle-aged to elderly African American community-based Jackson Heart Study. The distribution and correlates of CRP were analyzed for the entire study cohort who underwent the first examination (2001–2004). Heritability was estimated for the family cohort nested within the larger Jackson Heart Study (246 families, n=1,317). The relation between CRP and CVD risk factors were tested with multivariable stepwise regression analyses. Heritability was estimated using a variance components method. Linkage analysis was performed using the multipoint variance components approach. The study sample consisted of 4,919 participants (mean age 55±13 years, 63% women); median CRP concentration was 2.7 mg/L. In stepwise models traditional risk factors explained 23.8% of CRP’s variability, with body mass index (BMI, partial R2=13.6%) explaining 57.1% of the variability of CRP due to traditional risk factors. The heritability of CRP (adjusted for age, sex and BMI) was 0.45. The strongest linkage evidence for CRP was observed on chromosome 11 (11p13–11p11.2) with a logarithm of odds score of 2.72. In conclusion, in this large population-based cohort of African Americans, circulating CRP concentration was heritable and associated with several traditional cardiovascular risk factors, particularly BMI.
C-reactive protein; risk factors; genetics; heritability; blood pressure; cholesterol; body mass index; African Americans
C-reactive protein (CRP) is an acute phase reactant protein produced primarily by the liver. Circulating CRP levels are influenced by genetic and non-genetic factors, including infection and obesity. Genome-wide association studies (GWAS) provide an unbiased approach towards identifying loci influencing CRP levels. None of the six GWAS for CRP levels has been conducted in an African ancestry population. The present study aims to: (i) identify genetic variants that influence serum CRP in African Americans (AA) using a genome-wide association approach and replicate these findings in West Africans (WA), (ii) assess transferability of major signals for CRP reported in European ancestry populations (EA) to AA and (iii) use the weak linkage disequilibrium (LD) structure characteristic of African ancestry populations to fine-map the previously reported CRP locus. The discovery cohort comprised 837 unrelated AA, with the replication of significant single-nucleotide polymorphisms (SNPs) assessed in 486 WA. The association analysis was conducted with 2 366 856 genotyped and imputed SNPs under an additive genetic model with adjustment for appropriate covariates. Genome-wide and replication significances were set at P < 5 × 10−8 and P < 0.05, respectively. Ten SNPs in (CRP pseudogene-1) CRPP1 and CRP genes were associated with serum CRP (P = 2.4 × 10−09 to 4.3 × 10−11). All but one of the top-scoring SNPs associated with CRP in AA were successfully replicated in WA. CRP signals previously identified in EA samples were transferable to AAs, and we were able to fine-map this signal, reducing the region of interest from the 25 kb of LD around the locus in the HapMap CEU sample to only 8 kb in our AA sample.
High risk HPV (hrHPV) infection is a necessary cause of cervical cancer but the host genetic determinants of infection are poorly understood. We enrolled 267 women who presented to our cervical cancer screening program in Abuja, Nigeria between April 2012 and August 2012. We collected information on demographic characteristics, risk factors of cervical cancer and obtained samples of blood and cervical exfoliated cells from all participants. We used Roche Linear Array HPV Genotyping Test® to characterize the prevalent HPV according to manufacturer's instruction; Sequenom Mass Array to test 21 SNPs in genes/regions previously associated with hrHPV and regression models to examine independent factors associated with HPV infection. We considered a p<0.05 as significant because this is a replication study. There were 65 women with and 202 women without hrHPV infection. Under the allelic model, we found significant association between two SNPs, rs2305809 on RPS19 and rs2342700 on TYMS, and prevalent hrHPV infection. Multivariate analysis of hrHPV risk adjusted for age, body mass index, smoking, age of menarche, age at sexual debut, lifetime total number of sexual partners and the total number of pregnancies as covariates, yielded a p-value of 0.071 and 0.010 for rs2305809 and rs2342700, respectively. Our findings in this unique population suggest that a number of genetic risk variants for hrHPV are shared with other population groups. Definitive studies with larger sample sizes and using genome wide approaches are needed to understand the genetic architecture of hrHPV risk in multiple populations.
Rapid advances in high throughput genomic technologies and next generation sequencing are making medical genomic research more readily accessible and affordable, including the sequencing of patient and control whole genomes and exomes in order to elucidate genetic factors underlying disease. Over the next five years, the Human Heredity and Health in Africa (H3Africa) Initiative, funded by the Wellcome Trust (United Kingdom) and the National Institutes of Health (United States of America), will contribute greatly towards sequencing of numerous African samples for biomedical research.
Funding agencies and journals often require submission of genomic data from research participants to databases that allow open or controlled data access for all investigators. Access to such genotype-phenotype and pedigree data, however, needs careful control in order to prevent identification of individuals or families. This is particularly the case in Africa, where many researchers and their patients are inexperienced in the ethical issues accompanying whole genome and exome research; and where an historical unidirectional flow of samples and data out of Africa has created a sense of exploitation and distrust. In the current study, we analysed the implications of the anticipated surge of next generation sequencing data in Africa and the subsequent data sharing concepts on the protection of privacy of research subjects. We performed a retrospective analysis of the informed consent process for the continent and the rest-of-the-world and examined relevant legislation, both current and proposed. We investigated the following issues: (i) informed consent, including guidelines for performing culturally-sensitive next generation sequencing research in Africa and availability of suitable informed consent documents; (ii) data security and subject privacy whilst practicing data sharing; (iii) conveying the implications of such concepts to research participants in resource limited settings.
We conclude that, in order to meet the unique requirements of performing next generation sequencing-related research in African populations, novel approaches to the informed consent process are required. This will help to avoid infringement of privacy of individual subjects as well as to ensure that informed consent adheres to acceptable data protection levels with regard to use and transfer of such information.
African populations; Ethical; legal; and societal issues; Next generation sequencing; Whole genome and whole exome sequencing
Interleukins (ILs) are key mediators of the immune response and inflammatory process. Plasma levels of IL-10, IL-1Ra, and IL-6 are associated with metabolic conditions, show large inter-individual variations, and are under strong genetic control. Therefore, elucidation of the genetic variants that influence levels of these ILs provides useful insights into mechanisms of immune response and pathogenesis of diseases. We conducted a genome-wide association study (GWAS) of IL-10, IL-1Ra, and IL-6 levels in 707 non-diabetic African Americans using 5,396,780 imputed and directly genotyped single nucleotide polymorphisms (SNPs) with adjustment for gender, age, and body mass index. IL-10 levels showed genome-wide significant associations (p<5×10−8) with eight SNPs, the most significant of which was rs5743185 in thePMS1 gene (p=2.30×10−10). We tested replication of SNPs that showed genome-wide significance in 425 non-diabetic individuals from West Africa, and successfully replicated SNP rs17365948 in the YWHAZ gene (p=0.02). IL-1Ra levels showed suggestive associations with two SNPs in the ASB3 gene (p=2.55×10−7), 10 SNPs in the IL-1 gene family (IL1F5, IL1F8, IL1F10, and IL1Ra, p=1.04×10−6 to 1.75×10−6), and 23 SNPs near the IL1A gene (p=1.22×10−6 to 1.63×10−6). We also successfully replicated rs4251961 (p=0.009); this SNP was reported to be associated with IL-1Ra levels in a candidate gene study of Europeans. IL-6 levels showed genome-wide significant association with one SNP (RP11-314E23.1; chr6:133397598; p=8.63×10−9). To our knowledge, this is the first GWAS on IL-10, IL-1Ra, and IL-6 levels. Follow-up of these findings may provide valuable insight into the pathobiology of IL actions and dysregulations in inflammation and human diseases.
interleukin; interleukin-10; interleukin-1Ra; interleukin-6; genome-wide association study; African American
Systemic lupus erythematosus (SLE) is a chronic systemic autoimmune disease with variable clinical presentation frequently affecting the skin, joints, haemopoietic system, kidneys, lungs and central nervous system. It can be life threatening when major organs are involved. The full pathological and genetic mechanisms of this complex disease are yet to be elucidated; although roles have been described for environmental triggers such as sunlight, drugs and chemicals, and infectious agents. Cellular processes such as inefficient clearing of apoptotic DNA fragments and generation of autoantibodies have been implicated in disease progression. A diverse array of disease-associated genes and microRNA regulatory molecules that are dysregulated through polymorphism and copy number variation have also been identified; and an effect of ethnicity on susceptibility has been described.
Systemic lupus erythematosus; Autoimmunity; Genetic susceptibility; Apoptosis; dsDNA; Disease genes
It has been suggested that adiponectin may offer protection against the adverse health effects of obesity. In this study, we determined the prevalence of paradoxically high adiponectin or paradoxical hyperadiponectinemia (PHA) among obese African Americans and investigated its relationship with the metabolically healthy obese (MHO) phenotype.
Total adiponectin and metabolic markers including fasting glucose, insulin, serum lipids and obesity measures were determined in 822 unrelated participants from the Howard University Family Study (HUFS). Logistic regression models were used to evaluate the association between MHO phenotype and PHA while adjusting for relevant covariates.
Overall, men had significantly lower adiponectin levels than women. However, adiponectin level was associated with obesity measures, glucose, insulin and insulin resistance index in both men and women. Equal proportion of the obese male and female subjects (19.2%; 66/343) had PHA; these obese individuals with PHA had a healthier metabolic profile including higher HDL-cholesterol, lower insulin levels and smaller waist circumference and insulin levels compared to those without PHA. Also, 28% (96/343) of the study participants met the criteria of MHO phenotype. Interestingly, 42% (28/66) of the obese individuals with PHA also had the MHO phenotype. Finally, the MHO phenotype was associated with PHA in both men and women.
These findings confirm the presence of MHO in African Americans and demonstrate the association of PHA with the MHO phenotype. In all, our findings along with other published results provide evidence for a more systematic investigation of the mechanisms underlying the protective function of adiponectin and its potential therapeutic applications in human metabolic disorders.
Obesity; paradoxical hyperadiponectinemia (PHA); metabolically healthy obese phenotype (MHO)
Habitual levels of dietary sodium and potassium are correlated with age-related increases in blood pressure (BP) and likely play a role in this phenomenon. Although extensive published evidence exists from randomized trials, relatively few large-scale community surveys with multiple 24-hour urine collections have been reported. We obtained three 24-hour samples on 2,704 individuals from Nigeria, Jamaica and the US to evaluate patterns of intake and within-person relationships to blood pressure. The average (±s.d.) age and weight of participants across all three sites were 39.9±8.6 years and 76.1±21.2 kg, respectively, and 55% of the total participants were females. Sodium excretion increased across the East-West gradient (e.g., 123.9±54.6, 134.1±48.8, 176.6±71.0 (±s.d.) mmol, Nigeria, Jamaica and US, respectively), while potassium was essentially unchanged (e.g., 46.3±22.9, 40.7±16.1, 44.7±16.4 (±s.d.) mmol, respectively). In multivariate analyses both sodium (positively) and potassium (negatively) were strongly correlated with blood pressure (p < 0.001); quantitatively the association was stronger, and more consistent in each site individually, for potassium. Within-population day-to-day variation was also greater for sodium than for potassium. Among each population group a significant correlation was observed between sodium and urine volume, supporting the prior finding of sodium as a determinant of fluid intake in free-living individuals. These data confirm the consistency with the possible role of dietary electrolytes as hypertension risk factors, reinforcing the relevance of potassium in these populations.
blood pressure; sodium excretion; potassium excretion; African Diaspora
Malaria transmission requires the production of male and female gametocytes in the human host followed by fertilization and sporogonic development in the mosquito midgut. Although essential for the spread of malaria through the population, little is known about the initiation of gametocytogenesis in vitro or in vivo. Using a gametocyte-defective parasite line and genetic complementation, we show that Plasmodium falciparum
gametocyte development 1 gene (Pfgdv1), encoding a peri-nuclear protein, is critical for early sexual differentiation. Transcriptional analysis of Pfgdv1 negative and positive parasite lines identified a set of gametocytogenesis early genes (Pfge) that were significantly down-regulated (>10 fold) in the absence of Pfgdv1 and expression was restored after Pfgdv1 complementation. Progressive accumulation of Pfge transcripts during successive rounds of asexual replication in synchronized cultures suggests that gametocytes are induced continuously during asexual growth. Comparison of Pfge gene transcriptional profiles in patient samples divided the genes into two groups differing in their expression in mature circulating gametocytes and providing candidates to evaluate gametocyte induction and maturation separately in vivo. The expression profile of one of the early gametocyte specific genes, Pfge1, correlated significantly with asexual parasitemia, which is consistent with the ongoing induction of gametocytogenesis during asexual growth observed in vitro and reinforces the need for sustained transmission-blocking strategies to eliminate malaria.
As malaria control efforts move toward eradication it becomes increasingly important to develop interventions that block transmission. Consequently, advances are needed in our understanding of the production of gametocytes, which are required to transmit the disease. This report provides a first view of the initial stages of gametocytogenesis in vitro and in vivo and demonstrates that during each asexual replication cycle a subpopulation of parasites convert to gametocyte development providing a long transmission window. We also identify a gene that is critical for gametocyte production, P. falciparum
gametocyte development 1 (Pfgdv1) and a set of genes specifically expressed during early gametocytogenesis in P. falciparum (Pfge genes). The expression profile and peri-nuclear location of Pfgdv1 in a subpopulation of schizonts is consistent with a role in an early step in gametocytogenesis. The RNA levels of Pfgdv1 and the Pfge genes accumulated gradually over several asexual cycles in vitro suggesting ongoing gametocyte formation during asexual growth. The further evaluation of these genes in a cohort of malaria infected patients indicated they are good candidates for markers to distinguish ring stage parasites committed to gametocyte production from circulating mature gametocytes, allowing direct analysis of the initiation of sexual differentiation in vivo.
Serum urate concentrations are highly heritable and elevated serum urate is a key risk factor for gout. Genome-wide association studies (GWAS) of serum urate in African American (AA) populations are lacking. We conducted a meta-analysis of GWAS of serum urate levels and gout among 5820 AA and a large candidate gene study among 6890 AA and 21 708 participants of European ancestry (EA) within the Candidate Gene Association Resource Consortium. Findings were tested for replication among 1996 independent AA individuals, and evaluated for their association among 28 283 EA participants of the CHARGE Consortium. Functional studies were conducted using 14C-urate transport assays in mammalian Chinese hamster ovary cells. In the discovery GWAS of serum urate, three loci achieved genome-wide significance (P< 5.0 × 10−8): a novel locus near SGK1/SLC2A12 on chromosome 6 (rs9321453, P= 1.0 × 10−9), and two loci previously identified in EA participants, SLC2A9 (P= 3.8 × 10−32) and SLC22A12 (P= 2.1 × 10−10). A novel rare non-synonymous variant of large effect size in SLC22A12, rs12800450 (minor allele frequency 0.01, G65W), was identified and replicated (beta −1.19 mg/dl, P= 2.7 × 10−16). 14C-urate transport assays showed reduced urate transport for the G65W URAT1 mutant. Finally, in analyses of 11 loci previously associated with serum urate in EA individuals, 10 of 11 lead single-nucleotide polymorphisms showed direction-consistent association with urate among AA. In summary, we identified and replicated one novel locus in association with serum urate levels and experimentally characterize the novel G65W variant in URAT1 as a functional allele. Our data support the importance of multi-ethnic GWAS in the identification of novel risk loci as well as functional variants.
Podoconiosis is a tropical lymphedema resulting from long-term barefoot exposure to red-clay soil derived from volcanic rock. The World Health Organization recently designated it as a neglected tropical disease. Podoconiosis develops in only a subgroup of exposed people, and studies have shown familial clustering with high heritability (63%).
We conducted a genomewide association study of 194 case patients and 203 controls from southern Ethiopia. Findings were validated by means of family-based association testing in 202 family trios and HLA typing in 94 case patients and 94 controls.
We found a genomewide significant association of podoconiosis with the single-nucleotide polymorphism (SNP) rs17612858, located 5.8 kb from the HLA-DQA1 locus (in the allelic model: odds ratio, 2.44; 95% confidence interval [CI], 1.82 to 3.26; P = 1.42×10−9; and in the additive model: odds ratio, 2.19; 95% CI, 1.66 to 2.90; P = 3.44×10−8), and suggestive associations (P<1.0×10−5) with seven other SNPs in or near HLA-DQB1, HLA-DQA1, and HLA-DRB1. We confirmed these associations using family-based association testing. HLA typing showed the alleles HLA-DRB1*0701 (odds ratio, 2.00), DQA1*0201 (odds ratio, 1.91), and DQB1*0202 (odds ratio, 1.79) and the HLA-DRB1*0701–DQB1*0202 haplotype (odds ratio, 1.92) were risk variants for podoconiosis.
Association between variants in HLA class II loci with podoconiosis (a noncommuni-cable disease) suggests that the condition may be a T-cell–mediated inflammatory disease and is a model for gene–environment interactions that may be relevant to other complex genetic disorders. (Funded by the Wellcome Trust and others.)
A recent, large genome-wide association study (GWAS) of European ancestry individuals has identified multiple genetic variants influencing serum lipids. Studies of the transferability of these associations to African Americans remain few, an important limitation given interethnic differences in serum lipids and the disproportionate burden of lipid-associated metabolic diseases among African Americans.
We attempted to evaluate the transferability of 95 lipid-associated loci recently identified in European ancestry individuals to 887 non-diabetic, unrelated African Americans from a population-based sample in the Washington, DC area. Additionally, we took advantage of the generally reduced linkage disequilibrium among African ancestry populations in comparison to European ancestry populations to fine-map replicated GWAS signals.
We successfully replicated reported associations for 10 loci (CILP2/SF4, STARD3, LPL, CYP7A1, DOCK7/ANGPTL3, APOE, SORT1, IRS1, CETP, and UBASH3B). Through trans-ethnic fine-mapping, we were able to reduce associated regions around 75% of the loci that replicated.
Between this study and previous work in African Americans, 40 of the 95 loci reported in a large GWAS of European ancestry individuals also influence lipid levels in African Americans. While there is now evidence that the lipid-influencing role of a number of genetic variants is observed in both European and African ancestry populations, the still considerable lack of concordance highlights the importance of continued ancestry-specific studies to elucidate the genetic underpinnings of these traits.
Lipids; Genetics; African Americans; Genome-wide association study; Ethnicity
The incidence of chronic kidney disease varies by ethnic group in the USA, with African Americans displaying a two-fold higher rate than European Americans. One of the two defining variables underlying staging of chronic kidney disease is the glomerular filtration rate. Meta-analysis in individuals of European ancestry has identified 23 genetic loci associated with the estimated glomerular filtration rate (eGFR). We conducted a follow-up study of these 23 genetic loci using a population-based sample of 1,018 unrelated admixed African Americans. We included in our follow-up study two variants in APOL1 associated with end-stage kidney disease discovered by admixture mapping in admixed African Americans. To address confounding due to admixture, we estimated local ancestry at each marker and global ancestry. We performed regression analysis stratified by local ancestry and combined the resulting regression estimates across ancestry strata using an inverse variance-weighted fixed effects model. We found that 11 of the 24 loci were significantly associated with eGFR in our sample. The effect size estimates were not significantly different between the subgroups of individuals with two copies of African ancestry vs. two copies of European ancestry for any of the 11 loci. In contrast, allele frequencies were significantly different at 10 of the 11 loci. Collectively, the 11 loci, including four secondary signals revealed by conditional analyses, explained 14.2% of the phenotypic variance in eGFR, in contrast to the 1.4% explained by the 24 loci in individuals of European ancestry. Our findings provide insight into the genetic basis of variation in renal function among admixed African Americans.
Low levels of high-density cholesterol (HDLc) accompany chronic kidney disease, but the association between HDLc and the estimated glomerular filtration rate (eGFR) in the general population is unclear. We investigated the HDLc-eGFR association in nondiabetic Han Chinese (HC, n = 1100), West Africans (WA, n = 1497), and African Americans (AA, n = 1539).
There were significant differences by ancestry: HDLc was positively associated with eGFR in HC (β = 0.13, P < 0.0001), but negatively associated among African ancestry populations (WA: −0.19, P < 0.0001; AA: −0.09, P = 0.02). These differences were also seen in nationally-representative NHANES data (among European Americans: 0.09, P = 0.005; among African Americans −0.14, P = 0.03). To further explore the findings in African ancestry populations, we investigated the role of an African ancestry-specific nephropathy risk variant, rs73885319, in the gene encoding HDL-associated APOL1. Among AA, an inverse HDLc-eGFR association was observed only with the risk genotype (−0.38 versus 0.001; P = 0.03). This interaction was not seen in WA.
In summary, counter to expectation, an inverse HDLc-eGFR association was observed among those of African ancestry. Given the APOL1 × HDLc interaction among AA, genetic factors may contribute to this paradoxical association. Notably, these findings suggest that the unexplained mechanism by which APOL1 affects kidney-disease risk may involve HDLc.
Advances in technology and reduced costs are facilitating large-scale sequencing of genes and exomes as well as entire genomes. Recently, we described an approach based on haplotypes called SCARVA1 that enables the simultaneous analysis of the association between rare and common variants in disease etiology. Here, we describe an extension of SCARVA that evaluates individual markers instead of haplotypes. This modified method (SCARVAsnp) is implemented in four stages. First, all common variants in a pre-specified region (eg, gene) are evaluated individually. Second, a union procedure is used to combined all rare variants (RVs) in the index region, and the ratio of the log likelihood with one RV excluded to the log likelihood of a model with all the collapsed RVs is calculated. On the basis of previously-reported simulation studies,1 a likelihood ratio ≥1.3 is considered statistically significant. Third, the direction of the association of the removed RV is determined by evaluating the change in λ values with the inclusion and exclusion of that RV. Lastly, significant common and rare variants, along with covariates, are included in a final regression model to evaluate the association between the trait and variants in that region. We apply simulated and real data sets to show that the method is simple to use, computationally effcient, and that it can accurately identify both common and rare risk variants. This method overcomes several limitations of existing methods. For example, SCARVAsnp limits loss of statistical power by not including variants that are not associated with the trait of interest in the final model. Also, SCARVAsnp takes into consideration the direction of association by effectively modelling positively and negatively associated variants.
complex traits; rare and common variants
Deficiency of prolyl 3-hydroxylase 1, encoded by LEPRE1, causes recessive osteogenesis imperfecta. We previously identified a LEPRE1 mutation, exclusively in African Americans and contemporary West Africans. We hypothesized that this allele originated in West Africa and was introduced to the Americas with the Atlantic slave trade. We aimed to determine the frequency of carriers for this mutation among African Americans and West Africans, and the mutation origin and age.
Genomic DNA was screened for the mutation using PCR and restriction digestion, and a custom TaqMan genomic SNP assay. The mutation age was estimated using microsatellites and short tandem repeats spanning 4.2 Mb surrounding LEPRE1 in probands and carriers.
Approximately 0.4% of Mid-Atlantic African Americans carry this mutation, estimating recessive OI in 1/260,000 births in this population. In Nigeria and Ghana, 1.48% of unrelated individuals are heterozygous carriers, predicting 1/18,260 births will be affected with recessive OI, equal to the incidence of de novo dominant OI. The mutation was not detected in Africans from surrounding countries. All carriers shared a haplotype of 63-770 Kb, consistent with a single founder for this mutation. Using linkage disequilibrium analysis, the mutation was estimated to have originated between 650 and 900 years before present (1100-1350 C.E.).
We identified a West African founder mutation for recessive OI in LEPRE1. Nearly 1.5% of Ghanians and Nigerians are carriers. The age of this allele is consistent with introduction to North America via the Atlantic slave trade (1501 – 1867 C.E).
LEPRE1; osteogenesis imperfecta; founder mutation; West Africa