Publications

Domestic cats have a unique breeding history and can be used as models for human hereditary and infectious diseases. In the current era of genome-wide association studies, insights regarding linkage disequilibrium (LD) are essential for efficient association studies. The objective of this study is to investigate the extent of LD in the domestic cat, Felis silvestris catus, particularly within its breeds. A custom illumina GoldenGate Assay consisting of 1536 single nucleotide polymorphisms (SNPs) equally divided over ten 1 Mb chromosomal regions was developed, and genotyped across 18 globally recognized cat breeds and two distinct random bred populations. The pair-wise LD descriptive measure (r(2)) was calculated between the SNPs in each region and within each population independently. LD decay was estimated by determining the non-linear least-squares of all pair-wise estimates as a function of distance using established models. The point of 50% decay of r(2) was used to compare the extent of LD between breeds. The longest extent of LD was observed in the Burmese breed, where the distance at which r(2) ≈ 0.25 was ∼380 kb, comparable to several horse and dog breeds. The shortest extent of LD was found in the Siberian breed, with an r(2) ≈ 0.25 at approximately 17 kb, comparable to random bred cats and human populations. A comprehensive haplotype analysis was also conducted. The haplotype structure of each region within each breed mirrored the LD estimates. The LD of cat breeds largely reflects the breeds' population history and breeding strategies. Understanding LD in diverse populations will contribute to an efficient use of the newly developed SNP array for the cat in the design of genome-wide association studies, as well as to the interpretation of results for the fine mapping of disease and phenotypic traits.

Infection with Epstein-Barr virus (EBV) is highly prevalent worldwide, and it has been associated with infectious mononucleosis and severe diseases including Burkitt lymphoma, Hodgkin lymphoma, nasopharyngeal lymphoma, and lymphoproliferative disorders. Although EBV has been the focus of extensive research, much still remains unknown concerning what makes some individuals more sensitive to infection and to adverse outcomes as a result of infection. Here we use an integrative genomics approach in order to localize genetic factors influencing levels of Epstein Barr virus (EBV) nuclear antigen-1 (EBNA-1) IgG antibodies, as a measure of history of infection with this pathogen, in large Mexican American families. Genome-wide evidence of both significant linkage and association was obtained on chromosome 6 in the human leukocyte antigen (HLA) region and replicated in an independent Mexican American sample of large families (minimum p-value in combined analysis of both datasets is 1.4×10(-15) for SNPs rs477515 and rs2516049). Conditional association analyses indicate the presence of at least two separate loci within MHC class II, and along with lymphocyte expression data suggest genes HLA-DRB1 and HLA-DQB1 as the best candidates. The association signals are specific to EBV and are not found with IgG antibodies to 12 other pathogens examined, and therefore do not simply reveal a general HLA effect. We investigated whether SNPs significantly associated with diseases in which EBV is known or suspected to play a role (namely nasopharyngeal lymphoma, Hodgkin lymphoma, systemic lupus erythematosus, and multiple sclerosis) also show evidence of associated with EBNA-1 antibody levels, finding an overlap only for the HLA locus, but none elsewhere in the genome. The significance of this work is that a major locus related to EBV infection has been identified, which may ultimately reveal the underlying mechanisms by which the immune system regulates infection with this pathogen.

BACKGROUND: Atherosclerosis is the underlying cause of cardiovascular disease, the leading cause of morbidity and mortality in all American populations, including American Indians. Genetic factors play an important role in the pathogenesis of atherosclerosis. Although a single-nucleotide polymorphism (SNP) may explain only a small portion of variability in disease, the joint effect of multiple variants in a pathway on disease susceptibility could be large.

METHODS AND RESULTS: Using a gene-family analysis, we investigated the joint associations of 61 tag SNPs in 7 nicotinic acetylcholine receptor genes with subclinical atherosclerosis, as measured by carotid intima-media thickness and plaque score, in 3665 American Indians from 94 families recruited by the Strong Heart Family Study (SHFS). Although multiple SNPs showed marginal association with intima-media thickness and plaque score individually, only a few survived adjustments for multiple testing. However, simultaneously modeling of the joint effect of all 61 SNPs in 7 nicotinic acetylcholine receptor genes revealed significant association of the nicotinic acetylcholine receptor gene family with both intima-media thickness and plaque score independent of known coronary risk factors.

CONCLUSIONS: Genetic variants in the nicotinic acetylcholine receptor gene family jointly contribute to subclinical atherosclerosis in American Indians who participated in the SHFS. These variants may influence the susceptibility of atherosclerosis through pathways other than cigarette smoking per se.

BACKGROUND: Arsenic (III) methyltransferase (AS3MT) has been related to urine arsenic metabolites in association studies. Other genes might also play roles in arsenic metabolism and excretion.

OBJECTIVE: We evaluated genetic determinants of urine arsenic metabolites in American Indian adults from the Strong Heart Study (SHS).

METHODS: We evaluated heritability of urine arsenic metabolites [percent inorganic arsenic (%iAs), percent monomethylarsonate (%MMA), and percent dimethylarsinate (%DMA)] in 2,907 SHS participants with urine arsenic measurements and at least one relative within the cohort. We conducted a preliminary linkage analysis in a subset of 487 participants with available genotypes on approximately 400 short tandem repeat markers using a general pedigree variance component approach for localizing quantitative trait loci (QTL).

RESULTS: The medians (interquartile ranges) for %iAs, %MMA, and %DMA were 7.7% (5.4-10.7%), 13.6% (10.5-17.1%), and 78.4% (72.5-83.1%), respectively. The estimated heritability was 53% for %iAs, 50% for %MMA, and 59% for %DMA. After adjustment for sex, age, smoking, body mass index, alcohol consumption, region, and total urine arsenic concentrations, LOD [logarithm (to the base of 10) of the odds] scores indicated suggestive evidence for genetic linkage with QTLs influencing urine arsenic metabolites on chromosomes 5 (LOD = 2.03 for %iAs), 9 (LOD = 2.05 for %iAs and 2.10 for %MMA), and 11 (LOD = 1.94 for %iAs). A peak for %DMA on chromosome 10 within 2 Mb of AS3MT had an LOD of 1.80.

CONCLUSIONS: This population-based family study in American Indian communities supports a genetic contribution to variation in the distribution of arsenic metabolites in urine and, potentially, the involvement of genes other than AS3MT.

PURPOSE: The role of v-ATPases in cancer biology is being increasingly recognized. Yeast studies indicate that the tyrosine kinase inhibitor imatinib may interact with the v-ATPase genes and alter the course of cancer progression. Data from humans in this regard are lacking.

METHODS: We constructed 55 lymphoblastoid cell lines from pedigreed, cancer-free human subjects and treated them with IC20 concentration of imatinib mesylate. Using these cell lines, we (i) estimated the heritability and differential expression of 19 genes encoding several subunits of the v-ATPase protein in response to imatinib treatment; (ii) estimated the genetic similarity among these genes; and (iii) conducted a high-density scan to find cis-regulating genetic variation associated with differential expression of these genes.

RESULTS: We found that the imatinib response of the genes encoding v-ATPase subunits is significantly heritable and can be clustered to identify novel drug targets in imatinib therapy. Further, five of these genes were significantly cis-regulated and together represented nearly half-log fold change in response to imatinib (p = 0.0107) that was homogenous (p = 0.2598).

CONCLUSIONS: Our results proffer support to the growing view that personalized regimens using proton pump inhibitors or v-ATPase inhibitors may improve outcomes of imatinib therapy in various cancers.

OBJECTIVE: Several genome-wide association studies (GWAS) have demonstrated that common genetic variants contribute to obesity. However, studies of this complex trait have focused on ancestrally European populations, despite the high prevalence of obesity in some minority groups.

DESIGN AND METHODS: As part of the "Population Architecture using Genomics and Epidemiology (PAGE)" Consortium, we investigated the association between 13 GWAS-identified single-nucleotide polymorphisms (SNPs) and BMI and obesity in 69,775 subjects, including 6,149 American Indians, 15,415 African-Americans, 2,438 East Asians, 7,346 Hispanics, 604 Pacific Islanders, and 37,823 European Americans. For the BMI-increasing allele of each SNP, we calculated β coefficients using linear regression (for BMI) and risk estimates using logistic regression (for obesity defined as BMI ≥ 30) followed by fixed-effects meta-analysis to combine results across PAGE sites. Analyses stratified by racial/ethnic group assumed an additive genetic model and were adjusted for age, sex, and current smoking. We defined "replicating SNPs" (in European Americans) and "generalizing SNPs" (in other racial/ethnic groups) as those associated with an allele frequency-specific increase in BMI.

RESULTS: By this definition, we replicated 9/13 SNP associations (5 out of 8 loci) in European Americans. We also generalized 8/13 SNP associations (5/8 loci) in East Asians, 7/13 (5/8 loci) in African Americans, 6/13 (4/8 loci) in Hispanics, 5/8 in Pacific Islanders (5/8 loci), and 5/9 (4/8 loci) in American Indians.

CONCLUSION: Linkage disequilibrium patterns suggest that tagSNPs selected for European Americans may not adequately tag causal variants in other ancestry groups. Accordingly, fine-mapping in large samples is needed to comprehensively explore these loci in diverse populations.

A loss-of-function mutation (Q141K, rs2231142) in the ATP-binding cassette, subfamily G, member 2 gene (ABCG2) has been shown to be associated with serum uric acid levels and gout in Asians, Europeans, and European and African Americans; however, less is known about these associations in other populations. Rs2231142 was genotyped in 22,734 European Americans, 9,720 African Americans, 3,849 Mexican Americans, and 3,550 American Indians in the Population Architecture using Genomics and Epidemiology (PAGE) Study (2008-2012). Rs2231142 was significantly associated with serum uric acid levels (P = 2.37 × 10(-67), P = 3.98 × 10(-5), P = 6.97 × 10(-9), and P = 5.33 × 10(-4) in European Americans, African Americans, Mexican Americans, and American Indians, respectively) and gout (P = 2.83 × 10(-10), P = 0.01, and P = 0.01 in European Americans, African Americans, and Mexican Americans, respectively). Overall, the T allele was associated with a 0.24-mg/dL increase in serum uric acid level (P = 1.37 × 10(-80)) and a 1.75-fold increase in the odds of gout (P = 1.09 × 10(-12)). The association between rs2231142 and serum uric acid was significantly stronger in men, postmenopausal women, and hormone therapy users compared with their counterparts. The association with gout was also significantly stronger in men than in women. These results highlight a possible role of sex hormones in the regulation of ABCG2 urate transporter and its potential implications for the prevention, diagnosis, and treatment of hyperuricemia and gout.

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.

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.

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.