Publications

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.

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.

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.

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.

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.

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.

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.

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.

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.

Glucose-dependent insulinotropic polypeptide (GIP) has important actions on whole body metabolic function. GIP and its receptor are also present in the central nervous system and have been linked to neurotrophic actions. Metabolic effects of central nervous system GIP signaling have not been reported. We investigated whether centrally administered GIP could increase peripheral plasma GIP concentrations and influence the metabolic response to a mixed macronutrient meal in nonhuman primates. An infusion and sampling system was developed to enable continuous intracerebroventricular (ICV) infusions with serial venous sampling in conscious nonhuman primates. Male baboons (Papio sp.) that were healthy and had normal body weights (28.9 ± 2.1 kg) were studied (n = 3). Animals were randomized to receive continuous ICV infusions of GIP (20 pmol·kg-1·h-1) or vehicle before and over the course of a 300-min mixed meal test (15 kcal/kg, 1.5g glucose/kg) on two occasions. A significant increase in plasma GIP concentration was observed under ICV GIP infusion (66.5 ± 8.0 vs. 680.6 ± 412.8 pg/ml, P = 0.04) before administration of the mixed meal. Increases in postprandial, but not fasted, insulin (P = 0.01) and pancreatic polypeptide (P = 0.04) were also observed under ICV GIP. Effects of ICV GIP on fasted or postprandial glucagon, glucose, triglyceride, and free fatty acids were not observed. Our data demonstrate that central GIP signaling can promote increased plasma GIP concentrations independent of nutrient stimulation and increase insulin and pancreatic polypeptide responses to a mixed meal.