Statistical Genetic Approaches to Investigate Genotype-by-Environment Interaction: Review and Novel Extension of Models.

Statistical genetic models of genotype-by-environment (G×E) interaction can be divided into two general classes, one on G×E interaction in response to dichotomous environments (e.g., sex, disease-affection status, or presence/absence of an exposure) and the other in response to continuous environments (e.g., physical activity, nutritional measurements, or continuous socioeconomic measures). Here we develop a novel model to jointly account for dichotomous and continuous environments. We develop the model in terms of a joint genotype-by-sex (for the dichotomous environment) and genotype-by-social determinants of health (SDoH; for the continuous environment). Using this model, we show how a depression variable, as measured by the Beck Depression Inventory-II survey instrument, is not only underlain by genetic effects (as has been reported elsewhere) but is also significantly determined by joint G×Sex and G×SDoH interaction effects. This model has numerous applications leading to potentially transformative research on the genetic and environmental determinants underlying complex diseases.

Gene-by-Environment Interaction in Non-Alcoholic Fatty Liver Disease and Depression: The Role of Hepatic Transaminases.

Non-alcoholic fatty liver disease (NAFLD) encompasses a range of liver conditions, from benign fatty accumulation to severe fibrosis. The global prevalence of NAFLD has risen to 25-30%, with variations across ethnic groups. NAFLD may advance to hepatocellular carcinoma, increases cardiovascular risk, is associated with chronic kidney disease, and is an independent metabolic disease risk factor. Assessment methods for liver health include liver biopsy, magnetic resonance imaging, ultrasound, and vibration-controlled transient elastography (VCTE by FibroScan). Hepatic transaminases are cost-effective and minimally invasive liver health assessment methods options. This study focuses on the interaction between genetic factors underlying the traits (hepatic transaminases and the FibroScan results) on the one hand and the environment (depression) on the other. We examined 525 individuals at risk for metabolic disorders. We utilized variance components models and likelihood-based statistical inference to examine potential GxE interactions in markers of NAFLD, including aspartate aminotransferase (AST), alanine aminotransferase (ALT), and the AST/ALT ratio, and Vibration-Controlled Transient Elastography (VCTE by FibroScan). We calculated the Fibroscan-AST (FAST) score (a score that identifies the risk of progressive non-alcoholic steatohepatitis (NASH) and screened for depression using the Beck Depression Inventory-II (BDI-II). We identified significant G × E interactions for AST/ALT ratio × BDI-II, but not AST, ALT, or the FAST score. Our findings support that genetic factors play a role in hepatic transaminases, especially the AST/ALT ratio, with depression influencing this relationship. These insights contribute to understanding the complex interplay of genetics, environment, and liver health, potentially guiding future personalized interventions.

Non-alcoholic Fatty Liver Disease and Depression: Evidence for Genotype × Environment Interaction in Mexican Americans.

This study examines the impact of G × E interaction effects on non-alcoholic fatty liver disease (NAFLD) among Mexican Americans in the Rio Grande Valley (RGV) of South Texas. We examined potential G × E interaction using variance components models and likelihood-based statistical inference in the phenotypic expression of NAFLD, including hepatic steatosis and hepatic fibrosis (identified using vibration controlled transient elastography and controlled attenuation parameter measured by the FibroScan Device). We screened for depression using the Beck Depression Inventory-II (BDI-II). We identified significant G × E interactions for hepatic fibrosis × BDI-II. These findings provide evidence that genetic factors interact with depression to influence the expression of hepatic fibrosis.

Genetic variants in ALDH1L1 and GLDC influence the serine-to-glycine ratio in Hispanic children.

BACKGROUND: Glycine is a proteogenic amino acid that is required for numerous metabolic pathways, including purine, creatine, heme, and glutathione biosynthesis. Glycine formation from serine, catalyzed by serine hydroxy methyltransferase, is the major source of this amino acid in humans. Our previous studies in a mouse model have shown a crucial role for the 10-formyltetrahydrofolate dehydrogenase enzyme in serine-to-glycine conversion. OBJECTIVES: We sought to determine the genomic influence on the serine-glycine ratio in 803 Hispanic children from 319 families of the Viva La Familia cohort. METHODS: We performed a genome-wide association analysis for plasma serine, glycine, and the serine-glycine ratio in Sequential Oligogenic Linkage Analysis Routines while accounting for relationships among family members. RESULTS: All 3 parameters were significantly heritable (h2 = 0.22-0.78; P < 0.004). The strongest associations for the serine-glycine ratio were with single nucleotide polymorphisms (SNPs) in aldehyde dehydrogenase 1 family member L1 (ALDH1L1) and glycine decarboxylase (GLDC) and for glycine with GLDC (P < 3.5 × 10-8; effect sizes, 0.03-0.07). No significant associations were found for serine. We also conducted a targeted genetic analysis with ALDH1L1 exonic SNPs and found significant associations between the serine-glycine ratio and rs2886059 (ß = 0.68; SE, 0.25; P = 0.006) and rs3796191 (ß = 0.25; SE, 0.08; P = 0.003) and between glycine and rs3796191 (ß = -0.08; SE, 0.02; P = 0.0004). These exonic SNPs were further associated with metabolic disease risk factors, mainly adiposity measures (P < 0.006). Significant genetic and phenotypic correlations were found for glycine and the serine-glycine ratio with metabolic disease risk factors, including adiposity, insulin sensitivity, and inflammation-related phenotypes [estimate of genetic correlation = -0.37 to 0.35 (P < 0.03); estimate of phenotypic correlation = -0.19 to 0.13 (P < 0.006)]. The significant genetic correlations indicate shared genetic effects among glycine, the serine-glycine ratio, and adiposity and insulin sensitivity phenotypes. CONCLUSIONS: Our study suggests that ALDH1L1 and GLDC SNPs influence the serine-to-glycine ratio and metabolic disease risk.

The Prevalence of Glaucoma in the Jirel Ethnic Group of Nepal.

Glaucoma is one of the leading causes of blindness worldwide with individuals in Asia disproportionately affected. Using a cross-sectional study design as part of the Jiri Eye Study, we assessed the prevalence of glaucoma in the Jirel population of Nepal and provide new information on the occurrence of glaucoma in south central Asia. Over a four-year period, 2,042 members of the Jirel population, aged 18 years and older, underwent a detailed ocular examination. Glaucoma was diagnosed using the International Society of Geographical and Epidemiological Ophthalmology criteria. The mean (SD) age at exam was 42.3 (16.7) years and 54.1% of the sample was female. In the total sample, the mean (SD) intraocular pressure (IOP) and vertical cup-to-disc ratio (VCDR) was 14.55 (2.42) mmHg and 0.31 (0.15), respectively. The 97.5th and 99.5th percentile for IOP and VCDR was 20 mmHg and 22 mmHg, and 0.7 and 0.8, respectively. The overall prevalence of glaucoma in the population was 2.30% (n = 47). Of these 47 individuals, 37 (78.7%) had primary open angle glaucoma, 6 (12.8%) had primary angle closure glaucoma, and 4 (8.5%) had secondary glaucoma. There was a significant (p = 5.86×10-6) increase in the prevalence of glaucoma with increasing age overall and across glaucoma subtypes. Six individuals with glaucoma (12.8%) were blind in at least one eye. Of the individuals with glaucoma, 93.6% were previously undiagnosed. In individuals aged 40 years or older (n = 1057, 51.4% female), the mean (SD) IOP and VCDR was 14.39 (2.63) mmHg and 0.34 (0.16), respectively, and glaucoma prevalence was 4.16% (n = 44). The prevalence of glaucoma and undiagnosed disease is high in the Jirel population of Nepal. This study will inform strategies to minimize glaucoma-associated burden in Nepal.

Genetic variants and physical activity interact to affect bone density in Hispanic children.

BACKGROUND: Our aim was to investigate if moderate to vigorous physical activity (MVPA), calcium intake interacts with bone mineral density (BMD)-related single nucleotide polymorphisms (SNPs) to influence BMD in 750 Hispanic children (4-19y) of the cross-sectional Viva La Familia Study. METHODS: Physical activity and dietary intake were measured by accelerometers and multiple-pass 24 h dietary recalls, respectively. Total body and lumbar spine BMD were measured by dual energy X-ray absorptiometry. A polygenic risk score (PRS) was computed based on SNPs identified in published literature. Regression analysis was conducted with PRSs, MVPA and calcium intake with total body and lumbar spine BMD. RESULTS: We found evidence of statistically significant interaction effects between the PRS and MVPA on total body BMD and lumbar spine BMD (p < 0.05). Higher PRS was associated with a lower total body BMD (ß = - 0.040 ± 0.009, p = 1.1 × 10- 5) and lumbar spine BMD (ß = - 0.042 ± 0.013, p = 0.0016) in low MVPA group, as compared to high MVPA group (ß = - 0.015 ± 0.006, p = 0.02; ß = 0.008 ± 0.01, p = 0.4, respectively). DISCUSSION: The study indicated that calcium intake does not modify the relationship between genetic variants and BMD, while it implied physical activity interacts with genetic variants to affect BMD in Hispanic children. Due to limited sample size of our study, future research on gene by environment interaction on bone health and functional studies to provide biological insights are needed. CONCLUSIONS: Bone health in Hispanic children with high genetic risk for low BMD is benefitted more by MVPA than children with low genetic risk. Our results may be useful to predict disease risk and tailor dietary and physical activity advice delivery to people, especially children.

Genetic variants affecting bone mineral density and bone mineral content at multiple skeletal sites in Hispanic children.

CONTEXT: Osteoporosis is a major public health burden with significant economic costs. However, the correlates of bone health in Hispanic children are understudied. OBJECTIVE: We aimed to identify genetic variants associated with bone mineral density (BMD) and bone mineral content (BMC) at multiple skeletal sites in Hispanic children. METHODS: We conducted a cross-sectional genome-wide linkage analysis, genome-wide and exome-wide association analysis of BMD and BMC. The Viva La Familia Study is a family-based cohort with a total of 1030 Hispanic children (4-19 years old at baseline) conducted in Houston, TX. BMD and BMC were measured by Dual-energy X-ray absorptiometry. RESULTS: Significant heritability were observed for BMC and BMD at multiple skeletal sites ranging between 44 and 68% (P < 2.8 × 10-9). Significant evidence for linkage was found for BMD of pelvis and left leg on chromosome 7p14, lumbar spine on 20q13 and left rib on 6p21, and BMC of pelvis on chromosome 20q12 and total body on 14q22-23 (logarithm of odds score > 3). We found genome-wide significant association between BMC of right arm and rs762920 at PVALB (P = 4.6 × 10-8), and between pelvis BMD and rs7000615 at PTK2B (P = 7.4 × 10-8). Exome-wide association analysis revealed novel association of variants at MEGF10 and ABRAXAS2 with left arm and lumber spine BMC, respectively (P < 9 × 10-7). CONCLUSIONS: We identified novel loci associated with BMC and BMD in Hispanic children, with strongest evidence for PTK2B. These findings provide better understanding of bone genetics and shed light on biological mechanisms underlying BMD and BMC variation.

Fine mapping and identification of serum urate loci in American Indians: The Strong Heart Family Study.

While studies have reported genetic loci affecting serum urate (SU) concentrations, few studies have been conducted in minority populations. Our objective for this study was to identify genetic loci regulating SU in a multigenerational family-based cohort of American Indians, the Strong Heart Family Study (SHFS). We genotyped 162,718 single nucleotide polymorphisms (SNPs) in 2000 SHFS participants using an Illumina MetaboChip array. A genome-wide association analysis of SU was conducted using measured genotype analysis approach accounting for kinships in SOLAR, and meta-analysis in METAL. Our results showed strong association of SU with rs4481233, rs9998811, rs7696092 and rs13145758 (minor allele frequency (MAF) = 25-44%; P < 3 × 10-14) of solute carrier family 2, member 9 (SLC2A9) and rs41481455, rs2231142 and rs1481012 (MAF = 29%; p < 3 × 10-9) of ATP-binding cassette protein, subfamily G, member 2 (ABCG2). Carriers of G alleles of rs9998811, rs4148155 and rs1481012 and A alleles of rs4481233, rs7696092 and rs13145758 and rs2231142 had lower SU concentrations as compared to non-carriers. Genetic analysis of SU conditional on significant SLC2A9 and ABCG2 SNPs revealed new loci, nucleobindin 1 (NUCB1) and neuronal PAS domain protein 4 (NPAS4) (p <6× 10-6). To identify American Indian-specific SNPs, we conducted targeted sequencing of key regions of SLC2A9. A total of 233 SNPs were identified of which 89 were strongly associated with SU (p < 7.1 × 10-10) and 117 were American Indian specific. Analysis of key SNPs in cohorts of Mexican-mestizos, European, Indian and East Asian ancestries showed replication of common SNPs, including our lead SNPs. Our results demonstrate the association of SU with uric acid transporters in a minority population of American Indians and potential novel associations of SU with neuronal-related genes which warrant further investigation.

Glycated Serum Protein Genetics and Pleiotropy with Cardiometabolic Risk Factors.

Measurements of fasting glucose (FG) or glycated hemoglobin A1c (HbA1c) are two clinically approved approaches commonly used to determine glycemia, both of which are influenced by genetic factors. Obtaining accurate measurements of FG or HbA1c is not without its challenges, though. Measuring glycated serum protein (GSP) offers an alternative approach for assessing glycemia. The aim of this study was to estimate the heritability of GSP and GSP expressed as a percentage of total serum albumin (%GA) using a variance component approach and localize genomic regions (QTLs) that harbor genes likely to influence GSP and %GA trait variation in a large extended multigenerational pedigree from Jiri, Nepal (n = 1,800). We also performed quantitative bivariate analyses to assess the relationship between GSP or %GA and several cardiometabolic traits. Additive genetic effects significantly influence variation in GSP and %GA levels (p values: 1.15 × 10-5 and 3.39 × 10-5, respectively). We localized a significant (LOD score = 3.18) and novel GSP QTL on chromosome 11q, which has been previously linked to type 2 diabetes. Two common (MAF > 0.4) SNPs within the chromosome 11 QTL were associated with GSP (adjusted pvalue < 5.87 × 10-5): an intronic variant (rs10790184) in the DSCAML1 gene and a 3'UTR variant (rs8258) in the CEP164 gene. Significant positive correlations were observed between GSP or %GA and blood pressure, and lipid traits (p values: 0.0062 to 1.78 × 10-9). A significant negative correlation was observed between %GA and HDL cholesterol (p = 1.12 × 10-5). GSP is influenced by genetic factors and can be used to assess glycemia and diabetes risk. Thus, GSP measurements can facilitate glycemic studies when accurate FG and/or HbA1c measurements are difficult to obtain. GSP can also be measured from frozen blood (serum) samples, which allows the prospect of retrospective glycemic studies using archived samples.

Genetic Variants Related to Cardiometabolic Traits Are Associated to B Cell Function, Insulin Resistance, and Diabetes Among AmeriCan Indians: The Strong Heart Family Study.

Background: Genetic research may inform underlying mechanisms for disparities in the burden of type 2 diabetes mellitus among American Indians. Our objective was to assess the association of genetic variants in cardiometabolic candidate genes with B cell dysfunction via HOMA-B, insulin resistance via HOMA-IR, and type 2 diabetes mellitus in the Strong Heart Family Study (SHFS). Methods and Results: We examined the association of variants, previously associated with cardiometabolic traits (∼200,000 from Illumina Cardio MetaboChip), using mixed models of HOMA-B residuals corrected for HOMA-IR (cHOMA-B), log transformed HOMA-IR, and incident diabetes, adjusted for age, sex, population stratification, and familial relatedness. Center-specific estimates were combined using fixed effect meta-analyses. We used Bonferroni correction to account for multiple testing (P < 4.13 × 10-7). We also assessed the association between variants in candidate diabetes genes with these metabolic traits. We explored the top SNPs in an independent, replication sample from Southwestern Arizona. We identified significant associations with cHOMA-B for common variants at 26 loci of which 8 were novel (PRSS7, FCRL5, PEL1, LRP12, IGLL1, ARHGEF10, PARVA, FLJ16686). The most significant variant association with cHOMA-B was observed on chromosome 5 for an intergenic variant near PARP8 (rs2961831, P = 6.39 × 10-9). In the replication study, we found a signal at rs4607517 near GCK/YKT6 (P = 0.01). Variants near candidate diabetes genes (especially GCK and KCNQ1) were also nominally associated with HOMA-IR and cHOMA-B. Conclusion: We identified variants at novel loci and confirmed those at known candidate diabetes loci associations for cHOMA-B. This study also provided evidence for association of variants at KCNQ2, CTNAA2, and KCNQ1with cHOMA-B among American Indians. Further studies are needed to account for the high heritability of diabetes among the American Indian participants of the SHFS cohort.

Genetic variation underlying renal uric acid excretion in Hispanic children: the Viva La Familia Study.

BACKGROUND: Reduced renal excretion of uric acid plays a significant role in the development of hyperuricemia and gout in adults. Hyperuricemia has been associated with chronic kidney disease and cardiovascular disease in children and adults. There are limited genome-wide association studies associating genetic polymorphisms with renal urate excretion measures. Therefore, we investigated the genetic factors that influence the excretion of uric acid and related indices in 768 Hispanic children of the Viva La Familia Study. METHODS: We performed a genome-wide association analysis for 24-h urinary excretion measures such as urinary uric acid/urinary creatinine ratio, uric acid clearance, fractional excretion of uric acid, and glomerular load of uric acid in SOLAR, while accounting for non-independence among family members. RESULTS: All renal urate excretion measures were significantly heritable (p <2 × 10-6) and ranged from 0.41 to 0.74. Empirical threshold for genome-wide significance was set at p <1 × 10-7. We observed a strong association (p < 8 × 10-8) of uric acid clearance with a single nucleotide polymorphism (SNP) in zinc finger protein 446 (ZNF446) (rs2033711 (A/G), MAF: 0.30). The minor allele (G) was associated with increased uric acid clearance. Also, we found suggestive associations of uric acid clearance with SNPs in ZNF324, ZNF584, and ZNF132 (in a 72 kb region of 19q13; p <1 × 10-6, MAFs: 0.28-0.31). CONCLUSION: For the first time, we showed the importance of 19q13 region in the regulation of renal urate excretion in Hispanic children. Our findings indicate differences in inherent genetic architecture and shared environmental risk factors between our cohort and other pediatric and adult populations.

Association of Cardiometabolic Genes with Arsenic Metabolism Biomarkers in American Indian Communities: The Strong Heart Family Study (SHFS).

BACKGROUND: Metabolism of inorganic arsenic (iAs) is subject to inter-individual variability, which is explained partly by genetic determinants. OBJECTIVES: We investigated the association of genetic variants with arsenic species and principal components of arsenic species in the Strong Heart Family Study (SHFS). METHODS: We examined variants previously associated with cardiometabolic traits (~ 200,000 from Illumina Cardio MetaboChip) or arsenic metabolism and toxicity (670) among 2,428 American Indian participants in the SHFS. Urine arsenic species were measured by high performance liquid chromatography-inductively coupled plasma mass spectrometry (HPLC-ICP-MS), and percent arsenic species [iAs, monomethylarsonate (MMA), and dimethylarsinate (DMA), divided by their sum × 100] were logit transformed. We created two orthogonal principal components that summarized iAs, MMA, and DMA and were also phenotypes for genetic analyses. Linear regression was performed for each phenotype, dependent on allele dosage of the variant. Models accounted for familial relatedness and were adjusted for age, sex, total arsenic levels, and population stratification. Single nucleotide polymorphism (SNP) associations were stratified by study site and were meta-analyzed. Bonferroni correction was used to account for multiple testing. RESULTS: Variants at 10q24 were statistically significant for all percent arsenic species and principal components of arsenic species. The index SNP for iAs%, MMA%, and DMA% (rs12768205) and for the principal components (rs3740394, rs3740393) were located near AS3MT, whose gene product catalyzes methylation of iAs to MMA and DMA. Among the candidate arsenic variant associations, functional SNPs in AS3MT and 10q24 were most significant (p < 9.33 × 10-5). CONCLUSIONS: This hypothesis-driven association study supports the role of common variants in arsenic metabolism, particularly AS3MT and 10q24. Citation: Balakrishnan P, Vaidya D, Franceschini N, Voruganti VS, Gribble MO, Haack K, Laston S, Umans JG, Francesconi KA, Goessler W, North KE, Lee E, Yracheta J, Best LG, MacCluer JW, Kent J Jr., Cole SA, Navas-Acien A. 2017. Association of cardiometabolic genes with arsenic metabolism biomarkers in American Indian communities: the Strong Heart Family Study (SHFS). Environ Health Perspect 125:15-22; http://dx.doi.org/10.1289/EHP251.

Assessment of cardiovascular disease risk factors in a genetically homogenous population of Parsi Zoroastrians in the United States: A pilot study.

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.

Fatty acids linked to cardiovascular mortality are associated with risk factors.

BACKGROUND: Although saturated fatty acids (FAs) have been linked to cardiovascular mortality, it is not clear whether this outcome is attributable solely to their effects on low-density lipoprotein cholesterol (LDL-C) or whether other risk factors are also associated with FAs. The Western Alaskan Native population, with its rapidly changing lifestyles, shift in diet from unsaturated to saturated fatty acids and dramatic increase in cardiovascular disease (CVD), presents an opportunity to elucidate any associations between specific FAs and known CVD risk factors. OBJECTIVE: We tested the hypothesis that the specific FAs previously identified as related to CVD mortality are also associated with individual CVD risk factors. METHODS: In this community-based, cross-sectional study, relative proportions of FAs in plasma and red blood cell membranes were compared with CVD risk factors in a sample of 758 men and women aged ≥35 years. Linear regression analyses were used to analyze relations between specific FAs and CVD risk factors (LDL-C, high-density lipoprotein cholesterol, triglycerides, C-reactive protein, systolic blood pressure, diastolic blood pressure, heart rate, body mass index, fasting glucose and fasting insulin, 2-hour glucose and 2-hour insulin). RESULTS: The specific saturated FAs previously identified as related to CVD mortality, the palmitic and myristic acids, were adversely associated with most CVD risk factors, whereas unsaturated linoleic acid (18:2n-6) and the marine n-3 FAs were not associated or were beneficially associated with CVD risk factors. CONCLUSIONS: The results suggest that CVD risk factors are more extensively affected by individual FAs than hitherto recognized, and that risk for CVD, MI and stroke can be reduced by reducing the intake of palmitate, myristic acid and simple carbohydrates and improved by greater intake of linoleic acid and marine n-3 FAs.

Linkage Analysis of Urine Arsenic Species Patterns in the Strong Heart Family Study.

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.

Serum uric acid concentrations and SLC2A9 genetic variation in Hispanic children: the Viva La Familia Study.

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.

Genetics of kidney disease and related cardiometabolic phenotypes in Zuni Indians: the Zuni Kidney Project.

The objective of this study is to identify genetic factors associated with chronic kidney disease (CKD) and related cardiometabolic phenotypes among participants of the Genetics of Kidney Disease in Zuni Indians study. The study was conducted as a community-based participatory research project in the Zuni Indians, a small endogamous tribe in rural New Mexico. We recruited 998 members from 28 extended multigenerational families, ascertained through probands with CKD who had at least one sibling with CKD. We used the Illumina Infinium Human1M-Duo version 3.0 BeadChips to type 1.1 million single nucleotide polymorphisms (SNPs). Prevalence estimates for CKD, hyperuricemia, diabetes, and hypertension were 24%, 30%, 17% and 34%, respectively. We found a significant (p < 1.58 × 10(-7)) association for a SNP in a novel gene for serum creatinine (PTPLAD2). We replicated significant associations for genes with serum uric acid (SLC2A9), triglyceride levels (APOA1, BUD13, ZNF259), and total cholesterol (PVRL2). We found novel suggestive associations (p < 1.58 × 10(-6)) for SNPs in genes with systolic (OLFML2B), and diastolic blood pressure (NFIA). We identified a series of genes associated with CKD and related cardiometabolic phenotypes among Zuni Indians, a population with a high prevalence of kidney disease. Illuminating genetic variations that modulate the risk for these disorders may ultimately provide a basis for novel preventive strategies and therapeutic interventions.

Mapping of a blood pressure QTL on chromosome 17 in American Indians of the strong heart family study.

BACKGROUND: Blood pressure (BP) is a complex trait, with a heritability of 30 to 40%. Several genome wide associated BP loci explain only a small fraction of the phenotypic variation. Family studies can provide an important tool for gene discovery by utilizing trait and genetic transmission information among relative-pairs. We have previously described a quantitative trait locus at chromosome 17q25.3 influencing systolic BP in American Indians of the Strong Heart Family Study (SHFS). This locus has been reported to associate with variation in BP traits in family studies of Europeans, African Americans and Hispanics. METHODS: To follow-up persuasive linkage findings at this locus, we performed comprehensive genotyping in the 1-LOD unit support interval region surrounding this QTL using a multi-step strategy. We first genotyped 1,334 single nucleotide polymorphisms (SNPs) in 928 individuals from families that showed evidence of linkage for BP. We then genotyped a second panel of 306 SNPs in all SHFS participants (N = 3,807) for genes that displayed the strongest evidence of association in the region, and, in a third step, included additional genotyping to better cover the genes of interest and to interrogate plausible candidate genes in the region. RESULTS: Three genes had multiple SNPs marginally associated with systolic BP (TBC1D16, HRNBP3 and AZI1). In BQTN analysis, used to estimate the posterior probability that any variant in each gene had an effect on the phenotype, AZI1 showed the most prominent findings (posterior probability of 0.66). Importantly, upon correction for multiple testing, none of our study findings could be distinguished from chance. CONCLUSION: Our findings demonstrate the difficulty of follow-up studies of linkage studies for complex traits, particularly in the context of low powered studies and rare variants underlying linkage peaks.

Evidence of heterogeneity by race/ethnicity in genetic determinants of QT interval.

BACKGROUND: QT interval (QT) prolongation is an established risk factor for ventricular tachyarrhythmia and sudden cardiac death. Previous genome-wide association studies in populations of the European descent have identified multiple genetic loci that influence QT, but few have examined these loci in ethnically diverse populations. METHODS: Here, we examine the direction, magnitude, and precision of effect sizes for 21 previously reported SNPs from 12 QT loci, in populations of European (n = 16,398), African (n = 5,437), American Indian (n = 5,032), Hispanic (n = 1,143), and Asian (n = 932) descent as part of the Population Architecture using Genomics and Epidemiology (PAGE) study. Estimates obtained from linear regression models stratified by race/ethnicity were combined using inverse-variance weighted meta-analysis. Heterogeneity was evaluated using Cochran's Q test. RESULTS: Of 21 SNPs, 7 showed consistent direction of effect across all 5 populations, and an additional 9 had estimated effects that were consistent across 4 populations. Despite consistent direction of effect, 9 of 16 SNPs had evidence (P < 0.05) of heterogeneity by race/ethnicity. For these 9 SNPs, linkage disequilibrium plots often indicated substantial variation in linkage disequilibrium patterns among the various racial/ethnic groups, as well as possible allelic heterogeneity. CONCLUSIONS: These results emphasize the importance of analyzing racial/ethnic groups separately in genetic studies. Furthermore, they underscore the possible utility of trans-ethnic studies to pinpoint underlying casual variants influencing heritable traits such as QT.

Generalization of associations of kidney-related genetic loci to American Indians.

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.