Publications

2013

Taylor, Kira C, Cara L Carty, Logan Dumitrescu, Petra Bůžkova, Shelley A Cole, Lucia Hindorff, Fred R Schumacher, et al. (2013) 2013. “Investigation of Gene-by-Sex Interactions for Lipid Traits in Diverse Populations from the Population Architecture Using Genomics and Epidemiology Study.”. BMC Genetics 14: 33. https://doi.org/10.1186/1471-2156-14-33.

BACKGROUND: High-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglyceride (TG) levels are influenced by both genes and the environment. Genome-wide association studies (GWAS) have identified  100 common genetic variants associated with HDL-C, LDL-C, and/or TG levels, mostly in populations of European descent, but little is known about the modifiers of these associations. Here, we investigated whether GWAS-identified SNPs for lipid traits exhibited heterogeneity by sex in the Population Architecture using Genomics and Epidemiology (PAGE) study.

RESULTS: A sex-stratified meta-analysis was performed for 49 GWAS-identified SNPs for fasting HDL-C, LDL-C, and ln(TG) levels among adults self-identified as European American (25,013). Heterogeneity by sex was established when phet < 0.001. There was evidence for heterogeneity by sex for two SNPs for ln(TG) in the APOA1/C3/A4/A5/BUD13 gene cluster: rs28927680 (p(het) = 7.4 x 10(-7)) and rs3135506 (p(het) = 4.3 x 10(-4)one SNP in PLTP for HDL levels (rs7679; p(het) = 9.9 x 10(-4)), and one in HMGCR for LDL levels (rs12654264; p(het) = 3.1 x 10(-5)). We replicated heterogeneity by sex in five of seventeen loci previously reported by genome-wide studies (binomial p = 0.0009). We also present results for other racial/ethnic groups in the supplementary materials, to provide a resource for future meta-analyses.

CONCLUSIONS: We provide further evidence for sex-specific effects of SNPs in the APOA1/C3/A4/A5/BUD13 gene cluster, PLTP, and HMGCR on fasting triglyceride levels in European Americans from the PAGE study. Our findings emphasize the need for considering context-specific effects when interpreting genetic associations emerging from GWAS, and also highlight the difficulties in replicating interaction effects across studies and across racial/ethnic groups.

Rodríguez-Sánchez, I P, M L Garza-Rodríguez, K Mohamed-Noriega, V S Voruganti, M E Tejero, I Delgado-Enciso, D C Pérez-Ibave, et al. (2013) 2013. “Olfactomedin-Like 3 (OLFML3) Gene Expression in Baboon and Human Ocular Tissues: Cornea, Lens, Uvea, and Retina.”. Journal of Medical Primatology 42 (3): 105-11. https://doi.org/10.1111/jmp.12037.

BACKGROUND: Olfactomedin-like is a family of polyfunctional polymeric glycoproteins. This family has at least four members. One member of this family is OLFML3, which is preferentially expressed in placenta but is also detected in other adult tissues including the liver and heart. However, its orthologous rat gene is expressed in the iris, sclera, trabecular meshwork, retina, and optic nerve.

METHODS: OLFML3 messenger amplification was performed by RT-PCR from human and baboon ocular tissues. The products were cloned and sequenced.

RESULTS: We report OLFML3 expression in human and baboon eye. The full coding DNA sequence has 1221 bp, from which an open reading frame of 406 amino acid was obtained. The baboon OLFML3 gene nucleotidic sequence has 98% and amino acidic 99% similarity with humans.

CONCLUSIONS: OLFML3 gene expression in human and baboon ocular tissues and its high similarity make the baboon a powerful model to deduce the physiological and/or metabolic function of this protein in the eye.

Zhang, Lili, Petra Bůžkova, Christina L Wassel, Mary J Roman, Kari E North, Dana C Crawford, Jonathan Boston, et al. (2013) 2013. “Lack of Associations of Ten Candidate Coronary Heart Disease Risk Genetic Variants and Subclinical Atherosclerosis in Four US Populations: The Population Architecture Using Genomics and Epidemiology (PAGE) Study.”. Atherosclerosis 228 (2): 390-9. https://doi.org/10.1016/j.atherosclerosis.2013.02.038.

BACKGROUND: A number of genetic variants have been discovered by recent genome-wide association studies for their associations with clinical coronary heart disease (CHD). However, it is unclear whether these variants are also associated with the development of CHD as measured by subclinical atherosclerosis phenotypes, ankle brachial index (ABI), carotid artery intima-media thickness (cIMT) and carotid plaque.

METHODS: Ten CHD risk single nucleotide polymorphisms (SNPs) were genotyped in individuals of European American (EA), African American (AA), American Indian (AI), and Mexican American (MA) ancestry in the Population Architecture using Genomics and Epidemiology (PAGE) study. In each individual study, we performed linear or logistic regression to examine population-specific associations between SNPs and ABI, common and internal cIMT, and plaque. The results from individual studies were meta-analyzed using a fixed effect inverse variance weighted model.

RESULTS: None of the ten SNPs was significantly associated with ABI and common or internal cIMT, after Bonferroni correction. In the sample of 13,337 EA, 3809 AA, and 5353 AI individuals with carotid plaque measurement, the GCKR SNP rs780094 was significantly associated with the presence of plaque in AI only (OR = 1.32, 95% confidence interval: 1.17, 1.49, P = 1.08 × 10(-5)), but not in the other populations (P = 0.90 in EA and P = 0.99 in AA). A 9p21 region SNP, rs1333049, was nominally associated with plaque in EA (OR = 1.07, P = 0.02) and in AI (OR = 1.10, P = 0.05).

CONCLUSIONS: We identified a significant association between rs780094 and plaque in AI populations, which needs to be replicated in future studies. There was little evidence that the index CHD risk variants identified through genome-wide association studies in EA influence the development of CHD through subclinical atherosclerosis as assessed by cIMT and ABI across ancestries.

Zhu, Yun, Saroja Voruganti, Jue Lin, Tet Matsuguchi, Elizabeth Blackburn, Lyle G Best, Elisa T Lee, Jean W MacCluer, Shelley A Cole, and Jinying Zhao. (2013) 2013. “QTL Mapping of Leukocyte Telomere Length in American Indians: The Strong Heart Family Study.”. Aging 5 (9): 704-16.

Telomeres play a central role in cellular senescence and are associated with a variety of age-related disorders such as dementia, Alzheimer's disease and atherosclerosis. Telomere length varies greatly among individuals of the same age, and is heritable. Here we performed a genome-wide linkage scan to identify quantitative trait loci (QTL) influencing leukocyte telomere length (LTL) measured by quantitative PCR in 3,665 American Indians (aged 14-93 years) from 94 large, multi-generational families. All participants were recruited by the Strong Heart Family Study (SHFS), a prospective study to identify genetic factors for cardiovascular disease and its risk factors in American Indians residing in Oklahoma, Arizona and Dakota. LTL heritability was estimated to be between 51% and 62%, suggesting a strong genetic predisposition to interindividual variation of LTL in this population. Significant QTLs were localized to chromosome 13 (Logarithm of odds score (LOD)=3.9) at 13q12.11, to 18q22.2 (LOD=3.2) and to 3p14.1 (LOD=3.0) for Oklahoma. This is the first study to identify susceptibility loci influencing leukocyte telomere variation in American Indians, a minority group suffering from a disproportionately high rate of type 2 diabetes and other age-related disorders.

Dumitrescu, Logan, Cara L Carty, Nora Franceschini, Lucia A Hindorff, Shelley A Cole, Petra Bůžkova, Fredrick R Schumacher, et al. (2013) 2013. “Post-Genome-Wide Association Study Challenges for Lipid Traits: Describing Age As a Modifier of Gene-Lipid Associations in the Population Architecture Using Genomics and Epidemiology (PAGE) Study.”. Annals of Human Genetics 77 (5): 416-25. https://doi.org/10.1111/ahg.12027.

Numerous common genetic variants that influence plasma high-density lipoprotein cholesterol, low-density lipoprotein cholesterol (LDL-C), and triglyceride distributions have been identified via genome-wide association studies (GWAS). However, whether or not these associations are age-dependent has largely been overlooked. We conducted an association study and meta-analysis in more than 22,000 European Americans between 49 previously identified GWAS variants and the three lipid traits, stratified by age (males: <50 or ≥50 years of age; females: pre- or postmenopausal). For each variant, a test of heterogeneity was performed between the two age strata and significant Phet values were used as evidence of age-specific genetic effects. We identified seven associations in females and eight in males that displayed suggestive heterogeneity by age (Phet < 0.05). The association between rs174547 (FADS1) and LDL-C in males displayed the most evidence for heterogeneity between age groups (Phet = 1.74E-03, I(2) = 89.8), with a significant association in older males (P = 1.39E-06) but not younger males (P = 0.99). However, none of the suggestive modifying effects survived adjustment for multiple testing, highlighting the challenges of identifying modifiers of modest SNP-trait associations despite large sample sizes.

Farook, Vidya S, Dawn K Coletta, Sobha Puppala, Jennifer Schneider, Geetha Chittoor, Shirley L Hu, Deidre A Winnier, et al. (2013) 2013. “Linkage of Type 2 Diabetes on Chromosome 9p24 in Mexican Americans: Additional Evidence from the Veterans Administration Genetic Epidemiology Study (VAGES).”. Human Heredity 76 (1): 36-46. https://doi.org/10.1159/000354849.

OBJECTIVE: Type 2 diabetes (T2DM) is a complex metabolic disease and is more prevalent in certain ethnic groups such as the Mexican Americans. The goal of our study was to perform a genome-wide linkage (GWL) analysis to localize T2DM susceptibility loci in Mexican Americans.

METHODS: We used the phenotypic and genotypic data from 1,122 Mexican-American individuals (307 families) who participated in the Veterans Administration Genetic Epidemiology Study (VAGES). GWL analysis was performed using the variance components approach. Data from 2 additional Mexican-American family studies, the San Antonio Family Heart Study (SAFHS) and the San Antonio Family Diabetes/Gallbladder Study (SAFDGS), were combined with the VAGES data to test for improved linkage evidence.

RESULTS: After adjusting for covariate effects, T2DM was found to be under significant genetic influences (h2 = 0.62, p = 2.7 × 10(-6)). The strongest evidence for linkage of T2DM occurred between markers D9S1871 and D9S2169 on chromosome 9p24.2-p24.1 (LOD = 1.8). Given that we previously reported suggestive evidence for linkage of T2DM at this region also in SAFDGS, we found the significant and increased linkage evidence (LOD = 4.3, empirical p = 1.0 × 10(-5), genome-wide p = 1.6 × 10(-3)) for T2DM at the same chromosomal region, when we performed a GWL analysis of the VAGES data combined with the SAFHS and SAFDGS data.

CONCLUSION: Significant T2DM linkage evidence was found on chromosome 9p24 in Mexican Americans. Importantly, the chromosomal region of interest in this study overlaps with several recent genome-wide association studies involving T2DM-related traits. Given its overlap with such findings and our own initial T2DM association findings in the 9p24 chromosomal region, high throughput sequencing of the linked chromosomal region could identify the potential causal T2DM genes.

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.