Publications

2013

Pandey, Gouri Shankar, Chen Yanover, Lisa M Miller-Jenkins, Susan Garfield, Shelley A Cole, Joanne E Curran, Eric K Moses, et al. (2013) 2013. “Endogenous Factor VIII Synthesis from the Intron 22-Inverted F8 Locus May Modulate the Immunogenicity of Replacement Therapy for Hemophilia A.”. Nature Medicine 19 (10): 1318-24. https://doi.org/10.1038/nm.3270.

Neutralizing antibodies (inhibitors) to replacement factor VIII (FVIII, either plasma derived or recombinant) impair the effective management of hemophilia A. Individuals with hemophilia A due to major deletions of the FVIII gene (F8) lack antigenically cross-reactive material in their plasma ("CRM-negative"), and the prevalence of inhibitors in these individuals may be as high as 90%. Conversely, individuals with hemophilia A caused by F8 missense mutations are CRM-positive, and their overall prevalence of inhibitors is <10% (ref. 2). Individuals with the F8 intron 22 inversion (found in ∼50% of individuals with severe hemophilia A) have been grouped with the former on the basis of their genetic defect and CRM-negative status. However, only ∼20% of these individuals develop inhibitors. Here we demonstrate that the levels of F8 mRNA and intracellular FVIII protein in B lymphoblastoid cells and liver biopsies from individuals with the intron 22 inversion are comparable to those in healthy controls. These results support the hypothesis that most individuals with the intron 22 inversion are tolerized to FVIII and thus do not develop inhibitors. Furthermore, we developed a new pharmacogenetic algorithm that permits the stratification of inhibitor risk for individuals and subpopulations by predicting the immunogenicity of replacement FVIII using, as input, the number of putative T cell epitopes in the infused protein and the competence of major histocompatibility complex class II molecules to present such epitopes. This algorithm showed statistically significant accuracy in predicting the presence of inhibitors in 25 unrelated individuals with the intron 22 inversion.

Franceschini, N, K Haack, H H H Goring, V S Voruganti, S Laston, L Almasy, E T Lee, et al. (2013) 2013. “Epidemiology and Genetic Determinants of Progressive Deterioration of Glycaemia in American Indians: The Strong Heart Family Study.”. Diabetologia 56 (10): 2194-202. https://doi.org/10.1007/s00125-013-2988-8.

AIMS/HYPOTHESIS: Type 2 diabetes is a chronic, heterogeneous disease and a major risk factor for cardiovascular diseases. The underlying mechanisms leading to progression to type 2 diabetes are not fully understood and genetic tools may help to identify important pathways of glycaemic deterioration.

METHODS: Using prospective data on American Indians from the Strong Heart Family Study, we identified 373 individuals defined as progressors (diabetes incident cases), 566 individuals with transitory impaired fasting glucose (IFG) and 1,011 controls (normal fasting glycaemia at all visits). We estimated the heritability (h(2)) of the traits and the evidence for association with 16 known variants identified in type 2 diabetes genome-wide association studies.

RESULTS: We noted high h(2) for diabetes progression (h(2) = 0.65 ± 0.16, p = 2.7 × 10(-6)) but little contribution of genetic factors to transitory IFG (h(2) = 0.09 ± 0.10, p = 0.19) for models adjusted for multiple risk factors. At least three variants (in WFS1, TSPAN8 and THADA) were nominally associated with diabetes progression in age- and sex-adjusted analyses with estimates showing the same direction of effects as reported in the discovery European ancestry studies.

CONCLUSIONS/INTERPRETATION: Our findings do not exclude these loci for diabetes susceptibility in American Indians and suggest phenotypic heterogeneity of the IFG trait, which may have implications for genetic studies when diagnosis is based on a single time-point measure.

Chittoor, G, V S Farook, S Puppala, S P Fowler, J Schneider, T D Dyer, S A Cole, et al. (2013) 2013. “Localization of a Major Susceptibility Locus Influencing Preterm Birth.”. Molecular Human Reproduction 19 (10): 687-96. https://doi.org/10.1093/molehr/gat036.

Preterm birth (PTB) is a complex trait, but little is known regarding its major genetic determinants. The objective of this study is to localize genes that influence susceptibility to PTB in Mexican Americans (MAs), a minority population in the USA, using predominantly microfilmed birth certificate-based data obtained from the San Antonio Family Birth Weight Study. Only 1302 singleton births from 288 families with information on PTB and significant covariates were considered for genetic analysis. PTB is defined as a childbirth that occurs at <37 completed weeks of gestation, and the prevalence of PTB in this sample was 6.4%. An ∼10 cM genetic map was used to conduct a genome-wide linkage analysis using the program SOLAR. The heritability of PTB was high (h(2) ± SE: 0.75 ± 0.20) and significant (P = 4.5 × 10(-5)), after adjusting for the significant effects of birthweight and birth order. We found significant evidence for linkage of PTB (LOD = 3.6; nominal P = 2.3 × 10(-5); empirical P = 1.0 × 10(-5)) on chromosome 18q between markers D18S1364 and D18S541. Several other chromosomal regions (2q, 9p, 16q and 20q) were also potentially linked with PTB. A strong positional candidate gene in the 18q linked region is SERPINB2 or PAI-2, a member of the plasminogen activator system that is associated with various reproductive processes. In conclusion, to our knowledge, perhaps for the first time in MAs or US populations, we have localized a major susceptibility locus for PTB on chromosome 18q21.33-q23.

Bella, Jonathan N, Shelley A Cole, Sandy Laston, Laura Almasy, Anthony Comuzzie, Elisa T Lee, Lyle G Best, et al. (2013) 2013. “Genome-Wide Linkage Analysis of Carotid Artery Lumen Diameter: The Strong Heart Family Study.”. International Journal of Cardiology 168 (4): 3902-8. https://doi.org/10.1016/j.ijcard.2013.06.048.

BACKGROUND: A significant proportion of the variability in carotid artery lumen diameter is attributable to genetic factors.

METHODS: Carotid ultrasonography and genotyping were performed in the 3300 American Indian participants in the Strong Heart Family Study (SHFS) to identify chromosomal regions harboring novel genes associated with inter-individual variation in carotid artery lumen diameter. Genome-wide linkage analysis was conducted using standard variance component linkage methods, implemented in SOLAR, based on multipoint identity-by-descent matrices.

RESULTS: Genome-wide linkage analysis revealed a significant evidence for linkage for a locus for left carotid artery diastolic and systolic lumen diameters in Arizona SHFS participants on chromosome 7 at 120 cM (lod = 4.85 and 3.77, respectively, after sex and age adjustment, and lod = 3.12 and 2.72, respectively, after adjustment for sex, age, height, weight, systolic and diastolic blood pressure, diabetes mellitus and current smoking). Other regions with suggestive evidence of linkage for left carotid artery diastolic and systolic lumen diameter were found on chromosome 12 at 153 cM (lod = 2.20 and 2.60, respectively, after sex and age adjustment, and lod = 2.44 and 2.16, respectively, after full covariate adjustment) in Oklahoma SHFS participants; suggestive linkage for right carotid artery diastolic and systolic lumen diameter was found on chromosome 9 at 154 cM (lod = 2.72 and 3.19, respectively after sex and age adjustment, and lod = 2.36 and 2.21, respectively, after full covariate adjustment) in Oklahoma SHFS participants.

CONCLUSION: We found significant evidence for loci influencing carotid artery lumen diameter on chromosome 7q and suggestive linkage on chromosomes 12q and 9q.

Rubicz, Rohina, Jianhui Zhu, Sandra Laston, Shelley A Cole, Saroja Voruganti, Sven O E Ebbesson, Barbara Howard V, et al. (2013) 2013. “Statistical Genetic Analysis of Serological Measures of Common, Chronic Infections in Alaska Native Participants in the GOCADAN Study.”. Genetic Epidemiology 37 (7): 751-7. https://doi.org/10.1002/gepi.21745.

This paper describes genetic investigations of seroreactivity to five common infectious pathogens in the Genetics of Coronary Artery Disease in Alaska Natives (GOCADAN) study. Antibody titers and seroprevalence were available for 495 to 782 (depending on the phenotype) family members at two time points, approximately 15 years apart, for Chlamydophila pneumoniae, Helicobacter pylori, cytomegalovirus (CMV), herpes simplex virus 1 (HSV-1), and herpes simplex virus 2 (HSV-2). Seroprevalence rates indicate that infections with most of these pathogens are common (≥20% for all of them, >80% for H. pylori, CMV, and HSV-1). Seropositive individuals typically remain seropositive over time, with seroreversion rates of <1% to 10% over ∼15 years. Antibody titers were significantly heritable for most pathogens, with the highest estimate being 0.61 for C. pneumoniae. Significant genome-wide linkage evidence was obtained for C. pneumoniae on chromosome 15 (logarithm of odds, LOD score of 3.13). These results demonstrate that individual host genetic differences influence antibody measures of common infections in this population, and further investigation may elucidate the underlying immunological processes and genes involved.

Gribble, Matthew O, Venkata Saroja Voruganti, Cheryl D Cropp, Kevin A Francesconi, Walter Goessler, Jason G Umans, Ellen K Silbergeld, et al. (2013) 2013. “SLCO1B1 Variants and Urine Arsenic Metabolites in the Strong Heart Family Study.”. Toxicological Sciences : An Official Journal of the Society of Toxicology 136 (1): 19-25. https://doi.org/10.1093/toxsci/kft181.

Arsenic species patterns in urine are associated with risk for cancer and cardiovascular diseases. The organic anion transporter coded by the gene SLCO1B1 may transport arsenic species, but its association with arsenic metabolites in human urine has not yet been studied. The objective of this study is to evaluate associations of urine arsenic metabolites with variants in the candidate gene SLCO1B1 in adults from the Strong Heart Family Study. We estimated associations between % arsenic species biomarker traits and 5 single-nucleotide polymorphisms (SNPs) in the SLCO1B1 gene in 157 participants, assuming additive genetics. Linear regression models for each SNP accounted for kinships and were adjusted for sex, body mass index, and study center. The minor allele of rs1564370 was associated with lower %MMA (p = .0003) and higher %DMA (p = .0002), accounting for 8% of the variance for %MMA and 9% for %DMA. The rs1564370 minor allele homozygote frequency was 17% and the heterozygote frequency was 43%. The minor allele of rs2291075 was associated with lower %MMA (p = .0006) and higher %DMA (p = .0014), accounting for 7% of the variance for %MMA and 5% for %DMA. The frequency of rs2291075 minor allele homozygotes was 1% and of heterozygotes was 15%. Common variants in SLCO1B1 were associated with differences in arsenic metabolites in a preliminary candidate gene study. Replication of this finding in other populations and analyses with respect to disease outcomes are needed to determine whether this novel candidate gene is important for arsenic-associated disease risks.

Kochunov, P, J Charlesworth, A Winkler, L E Hong, T E Nichols, J E Curran, E Sprooten, et al. (2013) 2013. “Transcriptomics of Cortical Gray Matter Thickness Decline During Normal Aging.”. NeuroImage 82: 273-83. https://doi.org/10.1016/j.neuroimage.2013.05.066.

INTRODUCTION: We performed a whole-transcriptome correlation analysis, followed by the pathway enrichment and testing of innate immune response pathway analyses to evaluate the hypothesis that transcriptional activity can predict cortical gray matter thickness (GMT) variability during normal cerebral aging.

METHODS: Transcriptome and GMT data were available for 379 individuals (age range=28-85) community-dwelling members of large extended Mexican American families. Collection of transcriptome data preceded that of neuroimaging data by 17 years. Genome-wide gene transcriptome data consisted of 20,413 heritable lymphocytes-based transcripts. GMT measurements were performed from high-resolution (isotropic 800 μm) T1-weighted MRI. Transcriptome-wide and pathway enrichment analysis was used to classify genes correlated with GMT. Transcripts for sixty genes from seven innate immune pathways were tested as specific predictors of GMT variability.

RESULTS: Transcripts for eight genes (IGFBP3, LRRN3, CRIP2, SCD, IDS, TCF4, GATA3, and HN1) passed the transcriptome-wide significance threshold. Four orthogonal factors extracted from this set predicted 31.9% of the variability in the whole-brain and between 23.4 and 35% of regional GMT measurements. Pathway enrichment analysis identified six functional categories including cellular proliferation, aggregation, differentiation, viral infection, and metabolism. The integrin signaling pathway was significantly (p<10(-6)) enriched with GMT. Finally, three innate immune pathways (complement signaling, toll-receptors and scavenger and immunoglobulins) were significantly associated with GMT.

CONCLUSION: Expression activity for the genes that regulate cellular proliferation, adhesion, differentiation and inflammation can explain a significant proportion of individual variability in cortical GMT. Our findings suggest that normal cerebral aging is the product of a progressive decline in regenerative capacity and increased neuroinflammation.

Dumitrescu, Logan, Cara L Carty, Nora Franceschini, Lucia A Hindorff, Shelley A Cole, Petra Bůžkova, Fredrick R Schumacher, et al. (2013) 2013. “No Evidence of Interaction Between Known Lipid-Associated Genetic Variants and Smoking in the Multi-Ethnic PAGE Population.”. Human Genetics 132 (12): 1427-31. https://doi.org/10.1007/s00439-013-1375-3.

Genome-wide association studies (GWAS) have identified many variants that influence high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, and/or triglycerides. However, environmental modifiers, such as smoking, of these known genotype-phenotype associations are just recently emerging in the literature. We have tested for interactions between smoking and 49 GWAS-identified variants in over 41,000 racially/ethnically diverse samples with lipid levels from the Population Architecture Using Genomics and Epidemiology (PAGE) study. Despite their biological plausibility, we were unable to detect significant SNP × smoking interactions.

Voruganti, Venkata Saroja, Jack W Kent, Subrata Debnath, Shelley A Cole, Karin Haack, Harald H H Göring, Melanie A Carless, et al. (2013) 2013. “Genome-Wide Association Analysis Confirms and Extends the Association of SLC2A9 With Serum Uric Acid Levels to Mexican Americans.”. Frontiers in Genetics 4: 279. https://doi.org/10.3389/fgene.2013.00279.

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.

2012

Lee, M, A C Choh, K D Williams, V Schroeder, T D Dyer, J Blangero, S A Cole, et al. (2012) 2012. “Genome-Wide Linkage Scan for Quantitative Trait Loci Underlying Normal Variation in Heel Bone Ultrasound Measures.”. The Journal of Nutrition, Health & Aging 16 (1): 8-13.

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.