Publications

Increased serum uric acid (SUA) is a risk factor for gout and renal and cardiovascular disease (CVD). The purpose of this study was to identify genetic factors that affect the variation in SUA in 632 Mexican Americans participants of the San Antonio Family Heart Study (SAFHS). A genome-wide association (GWA) analysis was performed using the Illumina Human Hap 550K single nucleotide polymorphism (SNP) microarray. We used a linear regression-based association test under an additive model of allelic effect, while accounting for non-independence among family members via a kinship variance component. All analyses were performed in the software package SOLAR. SNPs rs6832439, rs13131257, and rs737267 in solute carrier protein 2 family, member 9 (SLC2A9) were associated with SUA at genome-wide significance (p < 1.3 × 10(-7)). The minor alleles of these SNPs had frequencies of 36.2, 36.2, and 38.2%, respectively, and were associated with decreasing SUA levels. All of these SNPs were located in introns 3-7 of SLC2A9, the location of the previously reported associations in European populations. When analyzed for association with cardiovascular-renal disease risk factors, conditional on SLC2A9 SNPs strongly associated with SUA, significant associations were found for SLC2A9 SNPs with BMI, body weight, and waist circumference (p < 1.4 × 10(-3)) and suggestive associations with albumin-creatinine ratio and total antioxidant status (TAS). The SLC2A9 gene encodes an urate transporter that has considerable influence on variation in SUA. In addition to the primary association locus, suggestive evidence (p < 1.9 × 10(-6)) for joint linkage/association (JLA) was found at a previously-reported urate quantitative trait locus (Logarithm of odds score = 3.6) on 3p26.3. In summary, our GWAS extends and confirms the association of SLC2A9 with SUA for the first time in a Mexican American cohort and also shows for the first time its association with cardiovascular-renal disease risk factors.

Quantitative ultrasound (QUS) traits are correlated with bone mineral density (BMD), but predict risk for future fracture independent of BMD. Only a few studies, however, have sought to identify specific genes influencing calcaneal QUS measures. The aim of this study was to conduct a genome-wide linkage scan to identify quantitative trait loci (QTL) influencing normal variation in QUS traits. QUS measures were collected from a total of 719 individuals (336 males and 383 females) from the Fels Longitudinal Study who have been genotyped and have at least one set of QUS measurements. Participants ranged in age from 18.0 to 96.6 years and were distributed across 110 nuclear and extended families. Using the Sahara ® bone sonometer, broadband ultrasound attenuation (BUA), speed of sound (SOS) and stiffness index (QUI) were collected from the right heel. Variance components based linkage analysis was performed on the three traits using 400 polymorphic short tandem repeat (STR) markers spaced approximately 10 cM apart across the autosomes to identify QTL influencing the QUS traits. Age, sex, and other significant covariates were simultaneously adjusted. Heritability estimates (h²) for the QUS traits ranged from 0.42 to 0.57. Significant evidence for a QTL influencing BUA was found on chromosome 11p15 near marker D11S902 (LOD = 3.11). Our results provide additional evidence for a QTL on chromosome 11p that harbors a potential candidate gene(s) related to BUA and bone metabolism.

A theory of aging holds that senescence is caused by a dysregulated nuclear factor kappa B (NF-κB) signal transduction network (STN). We adopted a systems genetics approach in our study of the NF-κB STN. Ingenuity Pathways Analysis (IPA) was used to identify gene/gene product interactions between NF-κB and the genes in our transcriptional profiling array. Principal components factor analysis (PCFA) was performed on a sub-network of 19 genes, including two initiators of the toll-like receptor (TLR) pathway, myeloid differentiation primary response gene (88) (MyD88) and TIR (Toll/interleukin-1 receptor)-domain-containing adapter-inducing interferon-β (TRIF). TLR pathways are either MyD88-dependent or TRIF-dependent. Therefore, we also performed PCFA on a subset excluding the MyD88 transcript, and on another subset excluding two TRIF transcripts. Using linkage analysis we found that each set gave rise to at least one factor with a logarithm of the odds (LOD) score greater than 3, two on chromosome 15 at 15q12 and 15q22.2, and another two on chromosome 17 at 17p13.3 and 17q25.3. We also found several suggestive signals (2<LOD score<3) at 1q32.1, 1q41, 2q34, 3q23, and 7p15.3. We are currently examining potential associations with single nucleotide polymorphisms within the 1-LOD intervals of our linkage signals.

BACKGROUND: Several genome scans have explored the linkage of chronic kidney disease phenotypes to chromosomic regions with disparate results. Genome scan meta-analysis (GSMA) is a quantitative method to synthesize linkage results from independent studies and assess their concordance.

METHODS: We searched PubMed to identify genome linkage analyses of renal function traits in humans, such as estimated glomerular filtration rate (GFR), albuminuria, serum creatinine concentration and creatinine clearance. We contacted authors for numerical data and extracted information from individual studies. We applied the GSMA nonparametric approach to combine results across 14 linkage studies for GFR, 11 linkage studies for albumin creatinine ratio, 11 linkage studies for serum creatinine and 4 linkage studies for creatinine clearance.

RESULTS: No chromosomal region reached genome-wide statistical significance in the main analysis which included all scans under each phenotype; however, regions on Chromosomes 7, 10 and 16 reached suggestive significance for linkage to two or more phenotypes. Subgroup analyses by disease status or ethnicity did not yield additional information.

CONCLUSIONS: While heterogeneity across populations, methodologies and study designs likely explain this lack of agreement, it is possible that linkage scan methodologies lack the resolution for investigating complex traits. Combining family-based linkage studies with genome-wide association studies may be a powerful approach to detect private mutations contributing to complex renal phenotypes.

BACKGROUND: Despite overwhelming evidence that major depression is highly heritable, recent studies have localized only a single depression-related locus reaching genome-wide significance and have yet to identify a causal gene. Focusing on family-based studies of quantitative intermediate phenotypes or endophenotypes, in tandem with studies of unrelated individuals using categorical diagnoses, should improve the likelihood of identifying major depression genes. However, there is currently no empirically derived statistically rigorous method for selecting optimal endophentypes for mental illnesses. Here, we describe the endophenotype ranking value, a new objective index of the genetic utility of endophenotypes for any heritable illness.

METHODS: Applying endophenotype ranking value analysis to a high-dimensional set of over 11,000 traits drawn from behavioral/neurocognitive, neuroanatomic, and transcriptomic phenotypic domains, we identified a set of objective endophenotypes for recurrent major depression in a sample of Mexican American individuals (n = 1122) from large randomly selected extended pedigrees.

RESULTS: Top-ranked endophenotypes included the Beck Depression Inventory, bilateral ventral diencephalon volume, and expression levels of the RNF123 transcript. To illustrate the utility of endophentypes in this context, each of these traits were utlized along with disease status in bivariate linkage analysis. A genome-wide significant quantitative trait locus was localized on chromsome 4p15 (logarithm of odds = 3.5) exhibiting pleiotropic effects on both the endophenotype (lymphocyte-derived expression levels of the RNF123 gene) and disease risk.

CONCLUSIONS: The wider use of quantitative endophenotypes, combined with unbiased methods for selecting among these measures, should spur new insights into the biological mechanisms that influence mental illnesses like major depression.

Serum sodium concentration is the clinical index of systemic water balance. Although disordered water balance is common and morbid, little is known about genetic effects on serum sodium concentration at the population level. Prior studies addressed only participants of European descent and either failed to demonstrate significant heritability or showed only modest effect. We investigated heritability of serum sodium concentration in large cohorts reflecting a range of races/ethnicities, including the Framingham Heart Study (FHS, non-Hispanic Caucasian), the Heredity and Phenotype Intervention Heart Study (HAPI, Amish Caucasian), the Jackson Heart Study (JHS, African American), the Strong Heart Family Study (SHFS, American Indian), and the Genetics of Kidney Disease in Zuni Indians Study (GKDZI, American Indian). Serum sodium was transformed for the osmotic effect of glucose, and participants with markedly elevated glucose or reduced estimated glomerular filtration rate (eGFR) were excluded. Using a standard variance components method, incorporating covariates of age, glucose, and eGFR, we found heritability to be high in African American and American Indian populations and much more modest in non-Hispanic Caucasian populations. Estimates among females increased after stratification on sex and were suggestive among female participants in FHS (0.18 ± 0.12, P = 0.057) and male participants in JHS (0.24 ± 0.16, P = 0.067) and statistically significant among female participants in JHS (0.44 ± 0.09, P = 1 × 10 ⁻⁷), SHFS (0.59 ± 0.05, P = 9.4 × 10⁻⁴⁶), and GKDZI (0.46 ± 0.15, P = 1.7 × 10⁻⁴), and male participants in HAPI (0.18 ± 0.12, P = 0.03) and SHFS (0.67 ± 0.07, P = 5.4 × 10⁻²⁶). Exclusion of diuretic users increased heritability among females and was significant in all cohorts where data were available. In aggregate, these data strongly support the heritability of systemic water balance and underscore sex and ethnicity-specific effects.

BACKGROUND: Genome-wide association studies (GWAS) have identified loci associated with ischemic stroke (IS) and cardiovascular disease (CVD) in European-descent individuals, but their replication in different populations has been largely unexplored.

METHODS AND RESULTS: Nine single nucleotide polymorphisms (SNPs) selected from GWAS and meta-analyses of stroke, and 86 SNPs previously associated with myocardial infarction and CVD risk factors, including blood lipids (high density lipoprotein [HDL], low density lipoprotein [LDL], and triglycerides), type 2 diabetes, and body mass index (BMI), were investigated for associations with incident IS in European Americans (EA) N=26 276, African-Americans (AA) N=8970, and American Indians (AI) N=3570 from the Population Architecture using Genomics and Epidemiology Study. Ancestry-specific fixed effects meta-analysis with inverse variance weighting was used to combine study-specific log hazard ratios from Cox proportional hazards models. Two of 9 stroke SNPs (rs783396 and rs1804689) were significantly associated with [corrected] IS hazard in AA; none were significant in this large EA cohort. Of 73 CVD risk factor SNPs tested in EA, 2 (HDL and triglycerides SNPs) were associated with IS. In AA, SNPs associated with LDL, HDL, and BMI were significantly associated with IS (3 of 86 SNPs tested). Out of 58 SNPs tested in AI, 1 LDL SNP was significantly associated with IS.

CONCLUSIONS: Our analyses showing lack of replication in spite of reasonable power for many stroke SNPs and differing results by ancestry highlight the need to follow up on GWAS findings and conduct genetic association studies in diverse populations. We found modest IS associations with BMI and lipids SNPs, though these findings require confirmation.

The high-density-lipoprotein-(HDL-) associated esterase paraoxonase 1 (PON1) is a likely contributor to the antioxidant and antiatherosclerotic capabilities of HDL. Two nonsynonymous mutations in the structural gene, PON1, have been associated with variation in activity levels, but substantial interindividual differences remain unexplained and are greatest for substrates other than the eponymous paraoxon. PON1 activity levels were measured for three substrates-organophosphate paraoxon, arylester phenyl acetate, and lactone dihydrocoumarin-in 767 Mexican American individuals from San Antonio, Texas. Genetic influences on activity levels for each substrate were evaluated by association with approximately one million single nucleotide polymorphism (SNPs) while conditioning on PON1 genotypes. Significant associations were detected at five loci including regions on chromosomes 4 and 17 known to be associated with atherosclerosis and lipoprotein regulation and loci on chromosome 3 that regulate ubiquitous transcription factors. These loci explain 7.8% of variation in PON1 activity with lactone as a substrate, 5.6% with the arylester, and 3.0% with paraoxon. In light of the potential importance of PON1 in preventing cardiovascular disease/events, these novel loci merit further investigation.

BACKGROUND AND PURPOSE: We hypothesized that the P-selectin (SELP) gene, localized to a region on chromosome 1q24, pleiotropically contributes to increased blood pressure and cerebral atrophy. We tested this hypothesis by performing genetic correlation analyses for 13 mRNA gene expression measures from P-selectin and 11 other genes located in 1q24 region and three magnetic resonance imaging derived indices of cerebral integrity.

METHODS: The subject pool consisted of 369 (219F; aged 28-85, average = 47.1 ± 12.7 years) normally aging, community-dwelling members of large extended Mexican-American families. Genetic correlation analysis decomposed phenotypic correlation coefficients into genetic and environmental components among 13 leukocyte-based mRNA gene expressions and three whole-brain and regional measurements of cerebral integrity: cortical gray matter thickness, fractional anisotropy of cerebral white matter, and the volume of hyperintensive WM lesions.

RESULTS: From the 13 gene expressions, significant phenotypic correlations were only found for the P- and L-selectin expression levels. Increases in P-selectin expression levels tracked with decline in cerebral integrity while the opposite trend was observed for L-selectin expression. The correlations for the P-selectin expression were driven by shared genetic factors, while the correlations with L-selectin expression were due to shared environmental effects.

CONCLUSION: This study demonstrated that P-selectin expression shared a significant variance with measurements of cerebral integrity and posits elevated P-selectin expression levels as a potential risk factor of hypertension-related cerebral atrophy.

The delta-5 and delta-6 desaturases (D5D and D6D), encoded by fatty acid desaturase 1 (FADS1) and 2 (FADS2) genes, respectively, are rate-limiting enzymes in the metabolism of ω-3 and ω-6 fatty acids. The objective of this study was to identify genes influencing variation in estimated D5D and D6D activities in plasma and erythrocytes in Alaskan Eskimos (n = 761) participating in the genetics of coronary artery disease in Alaska Natives (GOCADAN) study. Desaturase activity was estimated by product: precursor ratio of polyunsaturated fatty acids. We found evidence of linkage for estimated erythrocyte D5D (eD5D) on chromosome 11q12-q13 (logarithm of odds score = 3.5). The confidence interval contains candidate genes FADS1, FADS2, 7-dehydrocholesterol reductase (DHCR7), and carnitine palmitoyl transferase 1A, liver (CPT1A). Measured genotype analysis found association between CPT1A, FADS1, and FADS2 single-nucleotide polymorphisms (SNPs) and estimated eD5D activity (p-values between 10(-28) and 10(-5)). A Bayesian quantitative trait nucleotide analysis showed that rs3019594 in CPT1A, rs174541 in FADS1, and rs174568 in FADS2 had posterior probabilities > 0.8, thereby demonstrating significant statistical support for a functional effect on eD5D activity. Highly significant associations of FADS1, FADS2, and CPT1A transcripts with their respective SNPs (p-values between 10(-75) and 10(-7)) in Mexican Americans of the San Antonio Family Heart Study corroborated our results. These findings strongly suggest a functional role for FADS1, FADS2, and CPT1A SNPs in the variation in eD5D activity.