Publications

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD.

Although DNA methylation is now recognized as an important mediator of complex diseases, the extent to which the genetic basis of such diseases is accounted for by DNA methylation is unknown. In the setting of large, extended families representing a minority, high-risk population of the USA, we aimed to characterize the role of epigenome-wide DNA methylation in type 2 diabetes (T2D). Using Illumina HumanMethylation450 BeadChip arrays, we tested for association of DNA methylation at 446 356 sites with age, sex and phenotypic traits related to T2D in 850 pedigreed Mexican-American individuals. Robust statistical analyses showed that (i) 15% of the methylome is significantly heritable, with a median heritability of 0.14; (ii) DNA methylation at 14% of CpG sites is associated with nearby sequence variants; (iii) 22% and 3% of the autosomal CpG sites are associated with age and sex, respectively; (iv) 53 CpG sites were significantly associated with liability to T2D, fasting blood glucose and insulin resistance; (v) DNA methylation levels at five CpG sites, mapping to three well-characterized genes (TXNIP, ABCG1 and SAMD12) independently explained 7.8% of the heritability of T2D (vi) methylation at these five sites was unlikely to be influenced by neighboring DNA sequence variation. Our study has identified novel epigenetic indicators of T2D risk in Mexican Americans who have increased risk for this disease. These results provide new insights into potential treatment targets of T2D.

Arsenic toxicokinetics are important for disease risks in exposed populations, but genetic determinants are not fully understood. We examined urine arsenic species patterns measured by HPLC-ICPMS among 2189 Strong Heart Study participants 18 years of age and older with data on  400 genome-wide microsatellite markers spaced  10 cM and arsenic speciation (683 participants from Arizona, 684 from Oklahoma, and 822 from North and South Dakota). We logit-transformed % arsenic species (% inorganic arsenic, %MMA, and %DMA) and also conducted principal component analyses of the logit % arsenic species. We used inverse-normalized residuals from multivariable-adjusted polygenic heritability analysis for multipoint variance components linkage analysis. We also examined the contribution of polymorphisms in the arsenic metabolism gene AS3MT via conditional linkage analysis. We localized a quantitative trait locus (QTL) on chromosome 10 (LOD 4.12 for %MMA, 4.65 for %DMA, and 4.84 for the first principal component of logit % arsenic species). This peak was partially but not fully explained by measured AS3MT variants. We also localized a QTL for the second principal component of logit % arsenic species on chromosome 5 (LOD 4.21) that was not evident from considering % arsenic species individually. Some other loci were suggestive or significant for 1 geographical area but not overall across all areas, indicating possible locus heterogeneity. This genome-wide linkage scan suggests genetic determinants of arsenic toxicokinetics to be identified by future fine-mapping, and illustrates the utility of principal component analysis as a novel approach that considers % arsenic species jointly.

Populations and individuals differ in susceptibility to infections because of a number of factors, including host genetic variation. We previously demonstrated that differences in antibody titer, which reflect infection history, are significantly heritable. Here we attempt to identify the genetic factors influencing variation in these serological phenotypes. Blood samples from >1300 Mexican Americans were quantified for IgG antibody level against 12 common infections, selected on the basis of their reported role in cardiovascular disease risk: Chlamydia pneumoniae; Helicobacter pylori; Toxoplasma gondii; cytomegalovirus; herpes simplex I virus; herpes simplex II virus; human herpesvirus 6 (HHV6); human herpesvirus 8 (HHV8); varicella zoster virus; hepatitis A virus (HAV); influenza A virus; and influenza B virus. Pathogen-specific quantitative antibody levels were analyzed, as were three measures of pathogen burden. Genome-wide linkage and joint linkage and association analyses were performed using  1 million SNPs. Significant linkage (lod scores >3.0) was obtained for HHV6 (on chromosome 7), HHV8 (on chromosome 6), and HAV (on chromosome 13). SNP rs4812712 on chromosome 20 was significantly associated with C. pneumoniae (P=5.3 × 10(-8)). However, no genome-wide significant loci were obtained for the other investigated antibodies. We conclude that it is possible to localize host genetic factors influencing some of these antibody traits, but that further larger-scale investigations will be required to elucidate the genetic mechanisms contributing to variation in antibody levels.

Brain-derived neurotrophic factor (BDNF) plays a key role in energy balance. In population studies, SNPs of the BDNF locus have been linked to obesity, but the mechanism by which these variants cause weight gain is unknown. Here, we examined human hypothalamic BDNF expression in association with 44 BDNF SNPs. We observed that the minor C allele of rs12291063 is associated with lower human ventromedial hypothalamic BDNF expression (p < 0.001) and greater adiposity in both adult and pediatric cohorts (p values < 0.05). We further demonstrated that the major T allele for rs12291063 possesses a binding capacity for the transcriptional regulator, heterogeneous nuclear ribonucleoprotein D0B, knockdown of which disrupts transactivation by the T allele. Binding and transactivation functions are both disrupted by substituting C for T. These findings provide a rationale for BDNF augmentation as a targeted treatment for obesity in individuals who have the rs12291063 CC genotype.

BACKGROUND AND OBJECTIVES: CKD disproportionally affects American Indians, who similar to other populations, show genetic susceptibility to kidney outcomes. Recent studies have identified several loci associated with kidney traits, but their relevance in American Indians is unknown.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS: This study used data from a large, family-based genetic study of American Indians (the Strong Heart Family Study), which includes 94 multigenerational families enrolled from communities located in Oklahoma, the Dakotas, and Arizona. Individuals were recruited from the Strong Heart Study, a population-based study of cardiovascular disease in American Indians. This study selected 25 single nucleotide polymorphisms in 23 loci identified from recently published kidney-related genome-wide association studies in individuals of European ancestry to evaluate their associations with kidney function (estimated GFR; individuals 18 years or older, up to 3282 individuals) and albuminuria (urinary albumin to creatinine ratio; n=3552) in the Strong Heart Family Study. This study also examined the association of single nucleotide polymorphisms in the APOL1 region with estimated GFR in 1121 Strong Heart Family Study participants. GFR was estimated using the abbreviated Modification of Diet in Renal Disease Equation. Additive genetic models adjusted for age and sex were used.

RESULTS: This study identified significant associations of single nucleotide polymorphisms with estimated GFR in or nearby PRKAG2, SLC6A13, UBE2Q2, PIP5K1B, and WDR72 (P<2.1 × 10(-3) to account for multiple testing). Single nucleotide polymorphisms in these loci explained 2.2% of the estimated GFR total variance and 2.9% of its heritability. An intronic variant of BCAS3 was significantly associated with urinary albumin to creatinine ratio. APOL1 single nucleotide polymorphisms were not associated with estimated GFR in a single variant test or haplotype analyses, and the at-risk variants identified in individuals with African ancestry were not detected in DNA sequencing of American Indians.

CONCLUSION: This study extends the genetic associations of loci affecting kidney function to American Indians, a population at high risk of kidney disease, and provides additional support for a potential biologic relevance of these loci across ancestries.

Telomeres play a central role in cellular aging, and shorter telomere length has been associated with age-related disorders including diabetes. However, a causal link between telomere shortening and diabetes risk has not been established. In a well-characterized longitudinal cohort of American Indians participating in the Strong Heart Family Study, we examined whether leukocyte telomere length (LTL) at baseline predicts incident diabetes independent of known diabetes risk factors. Among 2,328 participants free of diabetes at baseline, 292 subjects developed diabetes during an average 5.5 years of follow-up. Compared with subjects in the highest quartile (longest) of LTL, those in the lowest quartile (shortest) had an almost twofold increased risk of incident diabetes (hazard ratio [HR] 1.83 [95% CI 1.26-2.66]), whereas the risk for those in the second (HR 0.87 [95% CI 0.59-1.29]) and the third (HR 0.95 [95% CI 0.65-1.38]) quartiles was statistically nonsignificant. These findings suggest a nonlinear association between LTL and incident diabetes and indicate that LTL could serve as a predictive marker for diabetes development in American Indians, who suffer from disproportionately high rates of diabetes.

Inorganic arsenic is methylated in the body by arsenic (III) methyltransferase (AS3MT). Arsenic methylation is thought to play a role in arsenic-related epigenetic phenomena, including aberrant DNA and histone methylation. However, it is unclear whether the promoter of the AS3MT gene, which codes for AS3MT, is differentially methylated as a function of arsenic exposure. In this study, we evaluated AS3MT promoter methylation according to exposure, assessed by urinary arsenic excretion in a stratified random sample of 48 participants from the Strong Heart Study who had urine arsenic measured at baseline and DNA available from 1989 to 1991 and 1998-1999. For this study, all data are from the 1989-1991 visit. We measured AS3MT promoter methylation at its 48 CpG loci by bisulphite sequencing. We compared mean  % methylation at each CpG locus by arsenic exposure group using linear regression adjusted for study centre, age and sex. A hypomethylated region in the AS3MT promoter was associated with higher arsenic exposure. In vitro, arsenic induced AS3MT promoter hypomethylation, and it increased AS3MT expression in human peripheral blood mononuclear cells. These findings may suggest that arsenic exposure influences the epigenetic regulation of a major arsenic metabolism gene.

OBJECTIVE: Gene × diet interaction plays an important role in atherosclerosis, an inflammatory disorder. Leukotrienes are the most potent inflammatory mediators, and genetic variants encoding leukotriene genes have been implicated in atherosclerosis. This study tests nutrigenetic interaction of a previously defined leukotriene haplotype on carotid artery hypertrophy and atherosclerosis in American Indians.

METHODS: This study included 3402 American Indians participating in the Strong Heart Family Study (SHFS). Carotid artery measurements, including intima-media thickness (IMT), vascular mass, and plaque, were assessed using ultrasound. Eleven tagSNPs in the leukotriene A4 hydrolase (LTA4H) gene were genotyped in all subjects. Main haplotype effect and haplotype × diet interaction were examined by generalized estimating equation, adjusting for known risk factors.

RESULTS: There was no significant main effect of haplotype or diet on any of the carotid artery measures. However, a previously defined LTA4H haplotype, called HapE, significantly interacted with dietary intake of n-3 and n-6 fatty acids on both IMT (P(HapE × n3) = 0.018, P(HapE × n6) = 0.040) and vascular mass (P(HapE × n3) = 0.012, P(HapE × n6) = 0.018), but not plaque. The direction of this nutrigenetic interaction on IMT was consistent with that reported in a recent study of Caucasian twins.

CONCLUSION: Dietary intake of polyunsaturated fatty acids significantly modifies the effect of a leukotriene haplotype on carotid artery hypertrophy but not atherosclerosis in American Indians, independent of established cardiovascular risk factors. Replication of nutrigenetic interaction in two distinct ethnic groups suggests the robustness and generalizability of our findings to diverse populations.

BACKGROUND: C-reactive protein (CRP) is a biomarker of inflammation. Genome-wide association studies (GWAS) have identified single-nucleotide polymorphisms (SNPs) associated with CRP concentrations and inflammation-related traits such as cardiovascular disease, type 2 diabetes mellitus, and obesity. We aimed to replicate previous CRP-SNP associations, assess whether these associations generalize to additional race/ethnicity groups, and evaluate inflammation-related SNPs for a potentially pleiotropic association with CRP.

METHODS AND RESULTS: We selected and analyzed 16 CRP-associated and 250 inflammation-related GWAS SNPs among 40 473 African American, American Indian, Asian/Pacific Islander, European American, and Hispanic participants from 7 studies collaborating in the Population Architecture using Genomics and Epidemiology (PAGE) study. Fixed-effect meta-analyses combined study-specific race/ethnicity-stratified linear regression estimates to evaluate the association between each SNP and high-sensitivity CRP. Overall, 18 SNPs in 8 loci were significantly associated with CRP (Bonferroni-corrected P<3.1×10(-3) for replication, P<2.0×10(-4) for pleiotropy): Seven of these were specific to European Americans, while 9 additionally generalized to African Americans (1), Hispanics (5), or both (3); 1 SNP was seen only in African Americans and Hispanics. Two SNPs in the CELSR2/PSRC1/SORT1 locus showed a potentially novel association with CRP: rs599839 (P=2.0×10(-6)) and rs646776 (P=3.1×10(-5)).

CONCLUSIONS: We replicated 16 SNP-CRP associations, 10 of which generalized to African Americans and/or Hispanics. We also identified potentially novel pleiotropic associations with CRP for two SNPs previously associated with coronary artery disease and/or low-density lipoprotein-cholesterol. These findings demonstrate the benefit of evaluating genotype-phenotype associations in multiple race/ethnicity groups and looking for pleiotropic relationships among SNPs previously associated with related phenotypes.