Publications

2015

Voruganti, Saroja, Sandra Laston, Karin Haack, Nitesh R Mehta, Shelley A Cole, Nancy F Butte, and Anthony G Comuzzie. (2015) 2015. “Serum Uric Acid Concentrations and SLC2A9 Genetic Variation in Hispanic Children: The Viva La Familia Study.”. The American Journal of Clinical Nutrition 101 (4): 725-32. https://doi.org/10.3945/ajcn.114.095364.

BACKGROUND: Elevated concentrations of serum uric acid are associated with increased risk of gout and renal and cardiovascular diseases. Genetic studies in adults have consistently identified associations of solute carrier family 2, member 9 (SLC2A9), polymorphisms with variation in serum uric acid. However, it is not known whether the association of serum uric acid with SLC2A9 polymorphisms manifests in children.

OBJECTIVE: The aim was to investigate whether variation in serum uric acid is under genetic influence and whether the association with SLC2A9 polymorphisms generalizes to Hispanic children of the Viva La Familia Study.

DESIGN: We conducted a genomewide association study with 1.1 million genetic markers in 815 children.

RESULTS: We found serum uric acid to be significantly heritable [h(2) ± SD = 0.45 ± 0.08, P = 5.8 × 10(-11)] and associated with SLC2A9 variants (P values between 10(-16) and 10(-7)). Several of the significantly associated polymorphisms were previously identified in studies in adults. We also found positive genetic correlations between serum uric acid and BMI z score (ρG = 0.45, P = 0.002), percentage of body fat (ρG = 0.28, P = 0.04), fat mass (ρG = 0.34, P = 0.02), waist circumference (ρG = 0.42, P = 0.003), and waist-to-height ratio (ρG = 0.46, P = 0.001).

CONCLUSIONS: Our results show that variation in serum uric acid in Hispanic children is under considerable genetic influence and is associated with obesity-related phenotypes. As in adults, genetic variation in SLC2A9 is associated with serum uric acid concentrations, an important biomarker of renal and cardiovascular disease risk, in Hispanic children.

Zhou, Yi, Hung-I H Chen, A L Lin, H Dang, Karin Haack, Shelley A Cole, Yufei Huang, Haiyang Yu, Yidong Chen, and Chih-Ko Yeh. (2015) 2015. “Early Gene Expression in Salivary Gland After Isoproterenol Treatment.”. Journal of Cellular Biochemistry 116 (3): 431-7. https://doi.org/10.1002/jcb.24995.

Isoproterenol, a β-adrenergic agonist, has been shown to induce salivary gland hyperplasia. However, the mechanism involved in this pharmacological phenomenon is not well understood. To gain a better understanding of the underlying changes, including genes, networks and pathways altered by isoproterenol, microarray-based gene expression analysis was conducted on rat parotid glands at 10, 30, and 60 min after isoproterenol injection. After isoproterenol treatment, the number of differentially expressed genes was increased in a time-dependent manner. Pathway analysis showed that cell hyperplasia, p38(MAPK), and IGF-1 were the most altered function, network and pathway, respectively. The balanced regulation of up- and down-expression of genes related to cell proliferation/survival may provide a better understanding of the mechanism of isoproterenol-induced parotid gland enlargement without tumor transformation.

González-Alvarez, Rafael, María de Lourdes Garza-Rodríguez, Iván Delgado-Enciso, Víctor Manuel Treviño-Alvarado, Ricardo Canales-Del-Castillo, Laura Elia Martínez-De-Villarreal, Ángel Lugo-Trampe, et al. (2015) 2015. “Molecular Evolution and Expression Profile of the Chemerine Encoding Gene RARRES2 in Baboon and Chimpanzee.”. Biological Research 48 (1): 31. https://doi.org/10.1186/s40659-015-0020-0.

BACKGROUND: Chemerin, encoded by the retinoic acid receptor responder 2 (RARRES2) gene is an adipocytesecreted protein with autocrine/paracrine functions in adipose tissue, metabolism and inflammation with a recently described function in vascular tone regulation, liver, steatosis, etc. This molecule is believed to represent a critical endocrine signal linking obesity to diabetes. There are no data available regarding evolution of RARRES2 in non-human primates and great apes. Expression profile and orthology in RARRES2 genes are unknown aspects in the biology of this multigene family in primates. Thus; we attempt to describe expression profile and phylogenetic relationship as complementary knowledge in the function of this gene in primates. To do that, we performed A RT-PCR from different tissues obtained during necropsies. Also we tested the hypotheses of positive evolution, purifying selection, and neutrality. And finally a phylogenetic analysis was made between primates RARRES2 protein.

RESULTS: RARRES2 transcripts were present in liver, lung, adipose tissue, ovary, pancreas, heart, hypothalamus and pituitary tissues. Expression in kidney and leukocytes were not detectable in either species. It was determined that the studied genes are orthologous.

CONCLUSIONS: RARRES2 evolution fits the hypothesis of purifying selection. Expression profiles of the RARRES2 gene are similar in baboons and chimpanzees and are also phylogenetically related.

Butte, Nancy F, Yan Liu, Issa F Zakeri, Robert P Mohney, Nitesh Mehta, Saroja Voruganti, Harald Göring, Shelley A Cole, and Anthony G Comuzzie. (2015) 2015. “Global Metabolomic Profiling Targeting Childhood Obesity in the Hispanic Population.”. The American Journal of Clinical Nutrition 102 (2): 256-67. https://doi.org/10.3945/ajcn.115.111872.

BACKGROUND: Metabolomics may unravel important biological pathways involved in the pathophysiology of childhood obesity.

OBJECTIVES: We aimed to 1) identify metabolites that differ significantly between nonobese and obese Hispanic children; 2) collapse metabolites into principal components (PCs) associated with obesity and metabolic risk, specifically hyperinsulinemia, hypertriglyceridemia, hyperleptinemia, and hyperuricemia; and 3) identify metabolites associated with energy expenditure and fat oxidation.

DESIGN: This trial was a cross-sectional observational study of metabolomics by using gas chromatography-mass spectrometry and ultrahigh-performance liquid chromatography-tandem mass spectrometry analyses performed on fasting plasma samples from 353 nonobese and 450 obese Hispanic children.

RESULTS: Branched-chained amino acids (BCAAs) (Leu, Ile, and Val) and their catabolites, propionylcarnitine and butyrylcarnitine, were significantly elevated in obese children. Strikingly lower lysolipids and dicarboxylated fatty acids were seen in obese children. Steroid derivatives were markedly higher in obese children as were markers of inflammation and oxidative stress. PC6 (BCAAs and aromatic AAs) and PC10 (asparagine, glycine, and serine) made the largest contributions to body mass index, and PC10 and PC12 (acylcarnitines) made the largest contributions to adiposity. Metabolic risk factors and total energy expenditure were associated with PC6, PC9 (AA and tricarboxylic acid cycle metabolites), and PC10. Fat oxidation was inversely related to PC8 (lysolipids) and positively related to PC16 (acylcarnitines).

CONCLUSIONS: Global metabolomic profiling in nonobese and obese children replicates the increased BCAA and acylcarnitine catabolism and changes in nucleotides, lysolipids, and inflammation markers seen in obese adults; however, a strong signature of reduced fatty acid catabolism and increased steroid derivatives may be unique to obese children. Metabolic flexibility in fuel use observed in obese children may occur through the activation of alternative intermediary pathways. Insulin resistance, hyperleptinemia, hypertriglyceridemia, hyperuricemia, and oxidative stress and inflammation evident in obese children are associated with distinct metabolomic profiles.

Iyengar, Sudha K, John R Sedor, Barry I Freedman, W H Linda Kao, Matthias Kretzler, Benjamin J Keller, Hanna E Abboud, et al. (2015) 2015. “Genome-Wide Association and Trans-Ethnic Meta-Analysis for Advanced Diabetic Kidney Disease: Family Investigation of Nephropathy and Diabetes (FIND).”. PLoS Genetics 11 (8): e1005352. https://doi.org/10.1371/journal.pgen.1005352.

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.

Kulkarni, Hemant, Mark Z Kos, Jennifer Neary, Thomas D Dyer, Jack W Kent, Harald H H Göring, Shelley A Cole, et al. (2015) 2015. “Novel Epigenetic Determinants of Type 2 Diabetes in Mexican-American Families.”. Human Molecular Genetics 24 (18): 5330-44. https://doi.org/10.1093/hmg/ddv232.

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.

Gribble, Matthew O, Venkata Saroja Voruganti, Shelley A Cole, Karin Haack, Poojitha Balakrishnan, Sandra L Laston, Maria Tellez-Plaza, et al. (2015) 2015. “Linkage Analysis of Urine Arsenic Species Patterns in the Strong Heart Family Study.”. Toxicological Sciences : An Official Journal of the Society of Toxicology 148 (1): 89-100. https://doi.org/10.1093/toxsci/kfv164.

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.

Rubicz, Rohina, Robert Yolken, Eugene Drigalenko, Melanie A Carless, Thomas D Dyer, Jack Kent, Joanne E Curran, et al. (2015) 2015. “Genome-Wide Genetic Investigation of Serological Measures of Common Infections.”. European Journal of Human Genetics : EJHG 23 (11): 1544-8. https://doi.org/10.1038/ejhg.2015.24.

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.

Mou, Zongyang, Thomas M Hyde, Barbara K Lipska, Keri Martinowich, Peter Wei, Chiew-Jen Ong, Lindsay A Hunter, et al. (2015) 2015. “Human Obesity Associated With an Intronic SNP in the Brain-Derived Neurotrophic Factor Locus.”. Cell Reports 13 (6): 1073-80. https://doi.org/10.1016/j.celrep.2015.09.065.

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.

2014

Franceschini, Nora, Karin Haack, Laura Almasy, Sandra Laston, Elisa T Lee, Lyle G Best, Richard R Fabsitz, et al. (2014) 2014. “Generalization of Associations of Kidney-Related Genetic Loci to American Indians.”. Clinical Journal of the American Society of Nephrology : CJASN 9 (1): 150-8. https://doi.org/10.2215/CJN.02300213.

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.