Publications

2017

Peng, H, F Yeh, J Lin, L G Best, S A Cole, E T Lee, B Howard V, and J Zhao. (2017) 2017. “Plasminogen Activator Inhibitor-1 Is Associated With Leukocyte Telomere Length in American Indians: Findings from the Strong Heart Family Study.”. Journal of Thrombosis and Haemostasis : JTH 15 (6): 1078-85. https://doi.org/10.1111/jth.13689.

UNLABELLED: Essentials Plasminogen activator inhibitor-1 (PAI-1) advanced cellular senescence in experiment studies. No population study exists on the association between PAI-1 and biological aging in American Indians. We found cross-sectional and longitudinal associations between higher PAI-1 and shorter telomere length. Our findings suggest a pathway linking PAI-1 with biological aging beyond metabolic factors.

SUMMARY: Background Plasminogen activator inhibitor-1 (PAI-1) promotes cellular aging both in vitro and in vivo. Telomere length is a marker of biological aging. Objectives To examine the cross-sectional and longitudinal associations between plasma PAI-1 and leukocyte telomere length in a large-scale epidemiological study of American Indians. Methods We measured leukocyte telomere length (LTL) and plasma PAI-1 in 2560 American Indians who were free of overt cardiovascular disease (CVD) and participated in the Strong Heart Family Study (SHFS) clinical examination in 2001-2003. LTL and PAI-1 were repeatedly measured in 475 participants who attended SHFS clinical visits in both 2001-2003 and 1998-1999. A generalized estimating equation model was used to examine the cross-sectional and longitudinal associations between PAI-1 and LTL, adjusting for known risk factors. Results A higher level of plasma PAI-1 was negatively associated with shorter age-adjusted LTL (β = -0.023; 95% CI, -0.034 to -0.013). This association was attenuated (β = -0.015; 95% CI, -0.029 to -0.002) after adjustments for demographics, study site, lifestyle (smoking, drinking and physical activity) and metabolic factors (obesity, blood pressure, fasting glucose, insulin, lipids and kidney function). Further adjustment for hsCRP did not change this association (β = -0.015; 95% CI, -0.029 to -0.001). Longitudinal analysis revealed that change in plasma PAI-1 was also inversely associated with change in LTL after adjusting for demographics, follow-up years, lifestyle factors, changes in metabolic factors, baseline levels of PAI-1 and LTL (β = -0.0005; 95% CI, -0.0009 to -0.0001). Conclusions A higher level of plasma PAI-1 was associated with shorter LTL in American Indians. This finding may suggest a potential role of PAI-1 in biological aging among American Indians.

Sabo, Aniko, Pamela Mishra, Shannon Dugan-Perez, Saroja Voruganti, Jack W Kent, Divya Kalra, Shelley A Cole, et al. (2017) 2017. “Exome Sequencing Reveals Novel Genetic Loci Influencing Obesity-Related Traits in Hispanic Children.”. Obesity (Silver Spring, Md.) 25 (7): 1270-76. https://doi.org/10.1002/oby.21869.

OBJECTIVE: To perform whole exome sequencing in 928 Hispanic children and identify variants and genes associated with childhood obesity.

METHODS: Single-nucleotide variants (SNVs) were identified from Illumina whole exome sequencing data using integrated read mapping, variant calling, and an annotation pipeline (Mercury). Association analyses of 74 obesity-related traits and exonic variants were performed using SeqMeta software. Rare autosomal variants were analyzed using gene-based association analyses, and common autosomal variants were analyzed at the SNV level.

RESULTS: (1) Rare exonic variants in 10 genes and 16 common SNVs in 11 genes that were associated with obesity traits in a cohort of Hispanic children were identified, (2) novel rare variants in peroxisome biogenesis factor 1 (PEX1) associated with several obesity traits (weight, weight z score, BMI, BMI z score, waist circumference, fat mass, trunk fat mass) were discovered, and (3) previously reported SNVs associated with childhood obesity were replicated.

CONCLUSIONS: Convergence of whole exome sequencing, a family-based design, and extensive phenotyping discovered novel rare and common variants associated with childhood obesity. Linking PEX1 to obesity phenotypes poses a novel mechanism of peroxisomal biogenesis and metabolism underlying the development of childhood obesity.

Haring, Bernhard, Adam Omidpanah, Astrid M Suchy-Dicey, Lyle G Best, Steven P Verney, Dean K Shibata, Shelley A Cole, et al. (2017) 2017. “Left Ventricular Mass, Brain Magnetic Resonance Imaging, and Cognitive Performance: Results From the Strong Heart Study.”. Hypertension (Dallas, Tex. : 1979) 70 (5): 964-71. https://doi.org/10.1161/HYPERTENSIONAHA.117.09807.

Left ventricular mass (LVM) has been shown to serve as a measure of target organ damage resulting from chronic exposure to several risk factors. Data on the association of midlife LVM with later cognitive performance are sparse. We studied 721 adults (mean age 56 years at baseline) enrolled in the Strong Heart Study (SHS, 1993-1995) and the ancillary CDCAI (Cerebrovascular Disease and Its Consequences in American Indians) Study (2010-2013), a study population with high prevalence of cardiovascular disease. LVM was assessed with transthoracic echocardiography at baseline in 1993 to 1995. Cranial magnetic resonance imaging and cognitive testing were undertaken between 2010 and 2013. Generalized estimating equations were used to model associations between LVM and later imaging and cognition outcomes. The mean follow-up period was 17 years. A difference of 25 g in higher LVM was associated with marginally lower hippocampal volume (0.01%; 95% confidence interval, 0.02-0.00; P=0.001) and higher white matter grade (0.10; 95% confidence interval, 0.02-0.18; P=0.014). Functionally, participants with higher LVM tended to have slightly lower scores on the modified mini-mental state examination (0.58; 95% confidence interval, 1.08-0.08; P=0.024). The main results persisted after adjusting for blood pressure levels or vascular disease. The small overall effect sizes are partly explained by survival bias because of the high prevalence of cardiovascular disease in our population. Our findings emphasize the role of cardiovascular health in midlife as a target for the prevention of deleterious cognitive and functional outcomes in later life.

Grau-Perez, Maria, Chin-Chi Kuo, Matthew O Gribble, Poojitha Balakrishnan, Miranda Jones Spratlen, Dhananjay Vaidya, Kevin A Francesconi, et al. (2017) 2017. “Association of Low-Moderate Arsenic Exposure and Arsenic Metabolism With Incident Diabetes and Insulin Resistance in the Strong Heart Family Study.”. Environmental Health Perspectives 125 (12): 127004. https://doi.org/10.1289/EHP2566.

BACKGROUND: High arsenic exposure has been related to diabetes, but at low-moderate levels the evidence is mixed. Arsenic metabolism, which is partly genetically controlled and may rely on certain B vitamins, plays a role in arsenic toxicity.

OBJECTIVE: We evaluated the prospective association of arsenic exposure and metabolism with type 2 diabetes and insulin resistance.

METHODS: We included 1,838 American Indian men and women free of diabetes (median age, 36 y). Arsenic exposure was assessed as the sum of inorganic arsenic (iAs), monomethylarsonate (MMA), and dimethylarsinate (DMA) urine concentrations (ΣAs). Arsenic metabolism was evaluated by the proportions of iAs, MMA, and DMA over their sum (iAs%, MMA%, and DMA%). Homeostasis model assessment for insulin resistance (HOMA2-IR) was measured at baseline and follow-up visits. Incident diabetes was evaluated at follow-up.

RESULTS: Median ΣAs, iAs%, MMA%, and DMA% was 4.4 μg/g creatinine, 9.5%, 14.4%, and 75.6%, respectively. Over 10,327 person-years of follow-up, 252 participants developed diabetes. Median HOMA2-IR at baseline was 1.5. The fully adjusted hazard ratio [95% confidence interval (CI)] for incident diabetes per an interquartile range increase in ΣAs was 1.57 (95% CI: 1.18, 2.08) in participants without prediabetes at baseline. Arsenic metabolism was not associated with incident diabetes. ΣAs was positively associated with HOMA2-IR at baseline but negatively with HOMA2-IR at follow-up. Increased MMA% was associated with lower HOMA2-IR when either iAs% or DMA% decreased. The association of arsenic metabolism with HOMA2-IR differed by B-vitamin intake and AS3MT genetics variants.

CONCLUSIONS: Among participants without baseline prediabetes, arsenic exposure was associated with incident diabetes. Low MMA% was cross-sectional and prospectively associated with higher HOMA2-IR. Research is needed to confirm possible interactions of arsenic metabolism with B vitamins and AS3MT variants on diabetes risk. https://doi.org/10.1289/EHP2566.

Fry, Constance L, Thomas C Naugle, Shelley A Cole, Jonathan Gelfond, Geetha Chittoor, Angeline F Mariani, Martin W Goros, Barrett G Haik, and Venkata Saroja Voruganti. (2017) 2017. “The Latino Eyelid: Anthropometric Analysis of a Spectrum of Findings.”. Ophthalmic Plastic and Reconstructive Surgery 33 (6): 440-45. https://doi.org/10.1097/IOP.0000000000000821.

PURPOSE: Published anthropometric measurements of the Latino eyelid are limited. This study describes features spanning the morphologic range from non-Latino whites to East Asians in the spectrum of the Latino eyelid.

METHODS: A cross-sectional study of 68 people (32 Latinos, 18 non-Latino whites, and 18 East Asians, ages 18-39), approved by the Institutional Review Board and HIPAA-compliant, was performed. Saliva samples determined genetic components. Indirect anthropometric measurements were performed with ImageJ software. Eyelid measurements included margin reflex distance, palpebral fissure height, eyelid crease height, orbital height, horizontal fissure length, inner and outer canthal distances, medial and lateral canthal angles, and lateral canthal angle of inclination. Additionally, exophthalmometry and epicanthal folds were recorded.

RESULTS: Analysis of 184 markers from HumanExome Chip data revealed distinct clustering patterns. Genetically, the Asian participants were in 1 group, the whites in another group, and the Latinos spanned the spectrum between these 2 groups. In Latinos, the inner canthal distance and lateral canthal angle of inclination were similar to Asians, whereas the eyelid crease spanned the range from Asians to whites. Half of the Latinos had epicanthal folds.

CONCLUSIONS: Latinos possess a spectrum of eyelid features spanning the morphologic characteristics from those of non-Latino whites to those of East Asians. These normative data on Latinos from Texas and Mexico aid in the diagnoses of Latino eyelid disorders and are a reference for optimizing oculofacial surgery outcomes.

2016

Mamtani, Manju, Hemant Kulkarni, Thomas D Dyer, Harald H H Göring, Jennifer L Neary, Shelley A Cole, Jack W Kent, et al. (2016) 2016. “Genome- and Epigenome-Wide Association Study of Hypertriglyceridemic Waist in Mexican American Families.”. Clinical Epigenetics 8: 6. https://doi.org/10.1186/s13148-016-0173-x.

BACKGROUND: There is growing interest in the hypertriglyceridemic waist (HTGW) phenotype, defined as high waist circumference (≥95 cm in males and ≥80 cm in females) combined with high serum triglyceride concentration (≥2.0 mmol/L in males and ≥1.5 mmol/L in females) as a marker of type 2 diabetes (T2D) and cardiovascular disease. However, the prevalence of this phenotype in high-risk populations, its association with T2D, and the genetic or epigenetic influences on HTGW are not well explored. Using data from large, extended families of Mexican Americans (a high-risk minority population in the USA) we aimed to: (1) estimate the prevalence of this phenotype, (2) test its association with T2D and related traits, and (3) dissect out the genetic and epigenetic associations with this phenotype using genome-wide and epigenome-wide studies, respectively.

RESULTS: Data for this study was from 850 Mexican American participants (representing 39 families) recruited under the ongoing San Antonio Family Heart Study, 26 % of these individuals had HTGW. This phenotype was significantly heritable (h (2) r = 0.52, p = 1.1 × 10(-5)) and independently associated with T2D as well as fasting glucose levels and insulin resistance. We conducted genome-wide association analyses using 759,809 single nucleotide polymorphisms (SNPs) and epigenome-wide association analyses using 457,331 CpG sites. There was no evidence of any SNP associated with HTGW at the genome-wide level but two CpG sites (cg00574958 and cg17058475) in CPT1A and one CpG site (cg06500161) in ABCG1 were significantly associated with HTGW and remained significant after adjusting for the closely related components of metabolic syndrome. CPT1A holds a cardinal position in the metabolism of long-chain fatty acids while ABCG1 plays a role in triglyceride metabolism.

CONCLUSIONS: Our results reemphasize the value of HTGW as a marker of T2D. This phenotype shows association with DNA methylation within CPT1A and ABCG1, genes involved in fatty acid and triglyceride metabolism. Our results underscore the importance of epigenetics in a clinically informative phenotype.

Traurig, Michael, Robert L Hanson, Alejandra Marinelarena, Sayuko Kobes, Paolo Piaggi, Shelley Cole, Joanne E Curran, et al. (2016) 2016. “Analysis of SLC16A11 Variants in 12,811 American Indians: Genotype-Obesity Interaction for Type 2 Diabetes and an Association With RNASEK Expression.”. Diabetes 65 (2): 510-9. https://doi.org/10.2337/db15-0571.

Genetic variants in SLC16A11 were recently reported to be associated with type 2 diabetes in Mexican and other Latin American populations. The diabetes risk haplotype had a frequency of 50% in Native Americans from Mexico but was rare in Europeans and Africans. In the current study, we analyzed SLC16A11 in 12,811 North American Indians and found that the diabetes risk haplotype, tagged by the rs75493593 A allele, was nominally associated with type 2 diabetes (P = 0.001, odds ratio 1.11). However, there was a strong interaction with BMI (P = 5.1 × 10(-7)) such that the diabetes association was stronger in leaner individuals. rs75493593 was also strongly associated with BMI in individuals with type 2 diabetes (P = 3.4 × 10(-15)) but not in individuals without diabetes (P = 0.77). Longitudinal analyses suggest that this is due, in part, to an association of the A allele with greater weight loss following diabetes onset (P = 0.02). Analyses of global gene expression data from adipose tissue, skeletal muscle, and whole blood provide evidence that rs75493593 is associated with expression of the nearby RNASEK gene, suggesting that RNASEK expression may mediate the effect of genotype on diabetes.

Chittoor, Geetha, Jack W Kent, Marcio Almeida, Sobha Puppala, Vidya S Farook, Shelley A Cole, Karin Haack, et al. (2016) 2016. “GWAS and Transcriptional Analysis Prioritize ITPR1 and CNTN4 for a Serum Uric Acid 3p26 QTL in Mexican Americans.”. BMC Genomics 17: 276. https://doi.org/10.1186/s12864-016-2594-5.

BACKGROUND: The variation in serum uric acid concentrations is under significant genetic influence. Elevated SUA concentrations have been linked to increased risk for gout, kidney stones, chronic kidney disease, and cardiovascular disease whereas reduced serum uric acid concentrations have been linked to multiple sclerosis, Parkinson's disease and Alzheimer's disease. Previously, we identified a novel locus on chromosome 3p26 affecting serum uric acid concentrations in Mexican Americans from San Antonio Family Heart Study. As a follow up, we examined genome-wide single nucleotide polymorphism data in an extended cohort of 1281 Mexican Americans from multigenerational families of the San Antonio Family Heart Study and the San Antonio Family Diabetes/Gallbladder Study. We used a linear regression-based joint linkage/association test under an additive model of allelic effect, while accounting for non-independence among family members via a kinship variance component.

RESULTS: Univariate genetic analysis indicated serum uric acid concentrations to be significant heritable (h (2) = 0.50 ± 0.05, p < 4 × 10(-35)), and linkage analysis of serum uric acid concentrations confirmed our previous finding of a novel locus on 3p26 (LOD = 4.9, p < 1 × 10(-5)) in the extended sample. Additionally, we observed strong association of serum uric acid concentrations with variants in following candidate genes in the 3p26 region; inositol 1,4,5-trisphosphate receptor, type 1 (ITPR1), contactin 4 (CNTN4), decapping mRNA 1A (DCP1A); transglutaminase 4 (TGM4) and rho guanine nucleotide exchange factor (GEF) 26 (ARHGEF26) [p < 3 × 10(-7); minor allele frequencies ranged between 0.003 and 0.42] and evidence of cis-regulation for ITPR1 transcripts.

CONCLUSION: Our results confirm the importance of the chromosome 3p26 locus and genetic variants in this region in the regulation of serum uric acid concentrations.

Vazquez-Vidal, Itzel, Geetha Chittoor, Sandra Laston, Sobha Puppala, Zeeba Kayani, Kaizeen Mody, Anthony G Comuzzie, Shelley A Cole, and Saroja Voruganti. (2016) 2016. “Assessment of Cardiovascular Disease Risk Factors in a Genetically Homogenous Population of Parsi Zoroastrians in the United States: A Pilot Study.”. American Journal of Human Biology : The Official Journal of the Human Biology Council 28 (3): 440-3. https://doi.org/10.1002/ajhb.22834.

OBJECTIVE: Genetically isolated and homogenous populations are ideal for detecting genes underlying common complex diseases. The use of isolated populations with reduced disease heterogeneity has led to significant gene discoveries in the past. The aim of this pilot study was to assess the prevalence of cardiovascular disease (CVD) risk phenotypes in a genetically homogenous population of Parsi Zoroastrians in the United States.

METHODS: Anthropometrics, blood pressure, and medical history were collected from 152 men and 186 women participating in a pilot study as part of the Parsi Family Study. The relative pairs used in the study included 60 parent-off springs, 28 siblings, 6 grandparent-grandchild, 7 avuncular, 18 half-siblings, 7 half-avuncular, and one half-first cousin. Estimates of genetic and environmental influence were calculated using a maximum likelihood-based variance components method implemented in SOLAR.

RESULTS: The prevalence of overweight/obesity in adults (62%) was on par with current US prevalence. Hypertension and prehypertension were prevalent in 16% and 46% of the participants, respectively. The quantitative genetic analysis revealed significant heritabilities for all anthropometric phenotypes (P < 0.05). Significant phenotypic correlations were found between blood pressure and anthropometric phenotypes (P < 0.001), whereas significant genetic correlation was found for only diastolic blood pressure and fat free mass (rhoG  = -0.88, P < 0.05).

CONCLUSION: These preliminary data show significant additive genetic effects on CVD-related phenotypes in this population. Our findings represent the first epidemiological data in Parsi Zoroastrians in the United States and offer excellent promise for future genetic studies in this population. Am. J. Hum. Biol. 28:440-443, 2016. © 2016 Wiley Periodicals, Inc.

Peng, Hao, Yun Zhu, Fawn Yeh, Shelley A Cole, Lyle G Best, Jue Lin, Elizabeth Blackburn, et al. (2016) 2016. “Impact of Biological Aging on Arterial Aging in American Indians: Findings from the Strong Heart Family Study.”. Aging 8 (8): 1583-92. https://doi.org/10.18632/aging.101013.

Telomere length, a marker of biological aging, has been associated with cardiovascular disease (CVD). Increased arterial stiffness, an indicator of arterial aging, predicts adverse CVD outcomes. However, the relationship between telomere length and arterial stiffness is less well studied. Here we examined the cross-sectional association between leukocyte telomere length (LTL) and arterial stiffness in 2,165 American Indians in the Strong Heart Family Study (SHFS). LTL was measured by qPCR. Arterial stiffness was assessed by stiffness index β. The association between LTL and arterial stiffness was assessed by generalized estimating equation model, adjusting for sociodemographics (age, sex, education level), study site, metabolic factors (fasting glucose, lipids, systolic blood pressure, and kidney function), lifestyle (BMI, smoking, drinking, and physical activity), and prevalent CVD. Results showed that longer LTL was significantly associated with a decreased arterial stiffness (β=-0.070, P=0.007). This association did not attenuate after further adjustment for hsCRP (β=-0.071, P=0.005) or excluding participants with overt CVD (β=-0.068, P=0.012), diabetes (β=-0.070, P=0.005), or chronic kidney disease (β=-0.090, P=0.001). In summary, shorter LTL was significantly associated with an increased arterial stiffness, independent of known risk factors. This finding may shed light on the potential role of biological aging in arterial aging in American Indians.