Publications

The current genetic and recombination maps of the cat have fewer than 3,000 markers and a resolution limit greater than 1 Mb. To complement the first-generation domestic cat maps, support higher resolution mapping studies, and aid genome assembly in specific areas as well as in the whole genome, a 15,000(Rad) radiation hybrid (RH) panel for the domestic cat was generated. Fibroblasts from the female Abyssinian cat that was used to generate the cat genomic sequence were fused to a Chinese hamster cell line (A23), producing 150 hybrid lines. The clones were initially characterized using 39 short tandem repeats (STRs) and 1,536 SNP markers. The utility of whole-genome amplification in preserving and extending RH panel DNA was also tested using 10 STR markers; no significant difference in retention was observed. The resolution of the 15,000(Rad) RH panel was established by constructing framework maps across 10 different 1-Mb regions on different feline chromosomes. In these regions, 2-point analysis was used to estimate RH distances, which compared favorably with the estimation of physical distances. The study demonstrates that the 15,000(Rad) RH panel constitutes a powerful tool for constructing high-resolution maps, having an average resolution of 40.1 kb per marker across the ten 1-Mb regions. In addition, the RH panel will complement existing genomic resources for the domestic cat, aid in the accurate re-assemblies of the forthcoming cat genomic sequence, and support cross-species genomic comparisons.

Individual differences in biological ageing (i.e., the rate of physiological response to the passage of time) may be due in part to genotype-specific variation in gene action. However, the sources of heritable variation in human age-related gene expression profiles are largely unknown. We have profiled genome-wide expression in peripheral blood mononuclear cells from 1240 individuals in large families and found 4472 human autosomal transcripts, representing  4349 genes, significantly correlated with age. We identified 623 transcripts that show genotype by age interaction in addition to a main effect of age, defining a large set of novel candidates for characterization of the mechanisms of differential biological ageing. We applied a novel SNP genotype × age interaction test to one of these candidates, the ubiquilin-like gene UBQLNL, and found evidence of joint cis-association and genotype by age interaction as well as trans-genotype by age interaction for UBQLNL expression. Both UBQLNL expression levels at recruitment and cis genotype are associated with longitudinal cancer risk in our study cohort.

BACKGROUND: Imatinib mesylate is currently the drug of choice to treat chronic myeloid leukemia. However, patient resistance and cytotoxicity make secondary lines of treatment, such as omacetaxine mepesuccinate, a necessity. Given that drug cytotoxicity represents a major problem during treatment, it is essential to understand the biological pathways affected to better predict poor drug response and prioritize a treatment regime.

METHODS: We conducted cell viability and gene expression assays to determine heritability and gene expression changes associated with imatinib and omacetaxine treatment of 55 non-cancerous lymphoblastoid cell lines, derived from 17 pedigrees. In total, 48,803 transcripts derived from Illumina Human WG-6 BeadChips were analyzed for each sample using SOLAR, whilst correcting for kinship structure.

RESULTS: Cytotoxicity within cell lines was highly heritable following imatinib treatment (h2 = 0.60-0.73), but not omacetaxine treatment. Cell lines treated with an IC20 dose of imatinib or omacetaxine showed differential gene expression for 956 (1.96%) and 3,892 transcripts (7.97%), respectively; 395 of these (0.8%) were significantly influenced by both imatinib and omacetaxine treatment. k-means clustering and DAVID functional annotation showed expression changes in genes related to kinase binding and vacuole-related functions following imatinib treatment, whilst expression changes in genes related to cell division and apoptosis were evident following treatment with omacetaxine. The enrichment scores for these ontologies were very high (mostly >10).

CONCLUSIONS: Induction of gene expression changes related to different pathways following imatinib and omacetaxine treatment suggests that the cytotoxicity of such drugs may be differentially tolerated by individuals based on their genetic background.

The prevalence of metabolic syndrome (MS) has been rising alarmingly worldwide, including in the United States, but knowledge on specific genetic determinants of MS is very limited. Therefore, we planned to identify the genetic determinants of MS as defined by National Cholesterol Education Program/Adult Treatment Panel III (NCEP/ATPIII) criteria. We performed linkage screen for MS using data from 692 Mexican Americans, who participated in the San Antonio Family Diabetes/Gallbladder Study (SAFDGS). We found strong evidence for linkage of MS on chromosome 7q (LOD = 3.6, empirical P = 6.0 × 10(-5)), between markers D7S2212 and D7S821. In addition, six chromosomal regions exhibited potential evidence for linkage (LOD ≥1.2) with MS. Furthermore, we examined 29 single-nucleotide polymorphisms (SNPs) from the fatty acid translocase (FAT or CD36, 18 SNPs) gene and guanine nucleotide binding protein, α transducing 3 (GNAT3, 11 SNPs) gene, located within the 1-LOD support interval region for their association with MS and its related traits. Several SNPs were associated with MS and its related traits. Remarkably, rs11760281 in GNAT3 and rs1194197 near CD36 exhibited the strongest associations with MS (P = 0.0003, relative risk (RR) = 1.6 and P = 0.004, RR = 1.7, respectively) and several other related traits. These two variants explained  18% of the MS linkage evidence on chromosome 7q21, and together conferred approximately threefold increase in MS risk (RR = 2.7). In conclusion, our linkage and subsequent association studies implicate a region on chromosome 7q21 to influence MS in Mexican Americans.

Genetic variants responsible for susceptibility to obesity and its comorbidities among Hispanic children have not been identified. The VIVA LA FAMILIA Study was designed to genetically map childhood obesity and associated biological processes in the Hispanic population. A genome-wide association study (GWAS) entailed genotyping 1.1 million single nucleotide polymorphisms (SNPs) using the Illumina Infinium technology in 815 children. Measured genotype analysis was performed between genetic markers and obesity-related traits i.e., anthropometry, body composition, growth, metabolites, hormones, inflammation, diet, energy expenditure, substrate utilization and physical activity. Identified genome-wide significant loci: 1) corroborated genes implicated in other studies (MTNR1B, ZNF259/APOA5, XPA/FOXE1 (TTF-2), DARC, CCR3, ABO); 2) localized novel genes in plausible biological pathways (PCSK2, ARHGAP11A, CHRNA3); and 3) revealed novel genes with unknown function in obesity pathogenesis (MATK, COL4A1). Salient findings include a nonsynonymous SNP (rs1056513) in INADL (p = 1.2E-07) for weight; an intronic variant in MTNR1B associated with fasting glucose (p = 3.7E-08); variants in the APOA5-ZNF259 region associated with triglycerides (p = 2.5-4.8E-08); an intronic variant in PCSK2 associated with total antioxidants (p = 7.6E-08); a block of 23 SNPs in XPA/FOXE1 (TTF-2) associated with serum TSH (p = 5.5E-08 to 1.0E-09); a nonsynonymous SNP (p = 1.3E-21), an intronic SNP (p = 3.6E-13) in DARC identified for MCP-1; an intronic variant in ARHGAP11A associated with sleep duration (p = 5.0E-08); and, after adjusting for body weight, variants in MATK for total energy expenditure (p = 2.7E-08) and in CHRNA3 for sleeping energy expenditure (p = 6.0E-08). Unprecedented phenotyping and high-density SNP genotyping enabled localization of novel genetic loci associated with the pathophysiology of childhood obesity.

Obesity is associated with a chronic low inflammatory state characterized by elevated levels of chemokines. Monocyte chemoattractant protein-1 (MCP-1) is a member of the cysteine-cysteine (CC) chemokine family and is increased in obesity. The purpose of this study was to identify loci regulating serum MCP-1 in obese Hispanic children from the Viva La Familia Study. A genome-wide association (GWA) analysis was performed in 815 children, ages 4-19 years, using genotypes assayed with the Illumina HumanOmni1-Quad v1.0 BeadChips. All analyses were performed in SOLAR using a linear regression-based test under an additive model of allelic effect, while accounting for the relatedness of family members via a kinship variance component. The strongest association for MCP-1 levels was found with a non-synonymous single nucleotide polymorphism (SNP), rs12075, resulting in an amino acid substitution (Asp42Gly) in the Duffy antigen receptor for chemokines (DARC) gene product (minor allele frequency=43.6%, p=1.3 × 10(-21)) on chromosome 1. Four other DARC SNPs were also significantly associated with MCP-1 levels (p<10(-16)-10(-6)). The Asp42Gly variant was associated with higher levels of MCP-1 and accounted for approximately 10% of its variability. In addition, MCP-1 levels were significantly associated with SNPs in chemokine receptor 3 (CCR3) and caspase recruitment domain family, member 9 (CARD9). In summary, the association of the DARC Asp42Gly variant with MCP-1 levels replicates previous GWA results substantiating a potential role for DARC in the regulation of pro-inflammatory cytokines.

OBJECTIVE: To identify genetic variation influencing serum bilirubin levels in American Indians, we performed genome-wide screening and association analyses in the Strong Heart Family Study. Bilirubin is an endogenous antioxidant that has demonstrated an inverse relationship with cardiovascular disease. Genetic variation within the promoter region of uridine diphosphate glucuronosyltransferase (UGT1A1) on chromosome 2q has been associated with elevated serum bilirubin levels in European populations. However, no study has investigated the UGT1A1 promoter in American Indians.

METHODS: Statistical analyses were carried out with 3,484 participants aged 14 to 93 years recruited from three geographic areas in the United States; Arizona, Oklahoma, and North and South Dakota.

RESULTS: Variance components linkage analysis detected a quantitative trait locus (QTL) for bilirubin on chromosome 2q in the combined centers (LOD = 6.61, P = 4.24 × 10⁻⁶) and in Oklahoma (LOD = 5.65, P = 4.57 24 × 10⁻⁵). Genetic association of the UGT1A1 promoter polymorphism was significant for all geographic locations. After adjustment using conditional linkage for UGT1A1 promoter variance, the linkage signal dropped to 1.10 in the combined sample and to 3.32 (P = 0.02) in Oklahoma, indicating this polymorphism is not completely responsible for the linkage signal in American Indians. We also detected suggestive linkage signals in the Dakotas on chromosome 10p12 (LOD = 2.18) and in the combined centers (LOD = 2.24) on chromosome 10q21.

CONCLUSIONS: Replication of a serum bilirubin QTL on chromosome 2q in American Indians implicates UGT1A1 but further genotyping is warranted to identify additional causative polymorphisms. Evidence also supports a potential novel locus for bilirubin on chromosome 10.

The genetic architecture of the craniofacial complex has been the subject of intense scrutiny because of the high frequency of congenital malformations. Numerous animal models have been used to document the early development of the craniofacial complex, but few studies have focused directly on the genetic underpinnings of normal variation in the human craniofacial complex. This study examines 80 quantitative traits derived from lateral cephalographs of 981 participants in the Fels Longitudinal Study, Wright State University, Dayton, Ohio. Quantitative genetic analyses were conducted using the Sequential Oligogenic Linkage Analysis Routines analytic platform, a maximum-likelihood variance components method that incorporates all familial information for parameter estimation. Heritability estimates were significant and of moderate to high magnitude for all craniofacial traits. Additionally, significant quantitative trait loci (QTL) were identified for 10 traits from the three developmental components (basicranium, splanchnocranium, and neurocranium) of the craniofacial complex. These QTL were found on chromosomes 3, 6, 11, 12, and 14. This study of the genetic architecture of the craniofacial complex elucidates fundamental information of the genetic architecture of the craniofacial complex in humans.

Variation in anthropometric measurements due to sexual dimorphism can be the result of genotype by sex interactions (G×S). The purpose of this study was to examine the sex-specific genetic architecture in anthropometric measurements in Alaskan Eskimos from the Genetics of Coronary Artery Disease in Alaska Natives (GOCADAN) study. Maximum likelihood-based variance components decomposition methods, implemented in SOLAR, were used for G×S analyses. Anthropometric measurements included BMI, waist circumference (WC), waist/height ratio, percent body fat (%BF), and subscapular and triceps skinfolds. Except for WC, mean values of all phenotypes were significantly different in men and women (P < 0.05). All anthropometric measures were significantly heritable (P < 0.001). In a preliminary analysis not allowing for G×S interaction, evidence of linkage was detected between markers D19S414 and D19S220 on chromosome 19 for WC (logarithm of odds (lod) = 3.5), %BF (lod = 1.7), BMI (lod = 2.4), waist/height ratio (lod = 2.5), subscapular (lod = 2.1), and triceps skinfolds (lod = 1.9). In subsequent analyses which allowed for G×S interaction, linkage was again found between these traits and the same two markers on chromosome 19 with significantly improved lod scores for: WC (lod = 4.5), %BF (lod = 3.8), BMI (lod = 3.5), waist/height ratio (lod = 3.2), subscapular (lod = 3.0), and triceps skinfolds (lod = 2.9). These results support the evidence of a G×S interaction in the expression of genetic effects resulting in sexual dimorphism in anthropometric phenotypes and identify the chromosome 19q12-13 region as important for adiposity-related traits in Alaskan Eskimos.