Publications

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.

For the past five years, genome-wide association studies (GWAS) have identified hundreds of common variants associated with human diseases and traits, including high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglyceride (TG) levels. Approximately 95 loci associated with lipid levels have been identified primarily among populations of European ancestry. The Population Architecture using Genomics and Epidemiology (PAGE) study was established in 2008 to characterize GWAS-identified variants in diverse population-based studies. We genotyped 49 GWAS-identified SNPs associated with one or more lipid traits in at least two PAGE studies and across six racial/ethnic groups. We performed a meta-analysis testing for SNP associations with fasting HDL-C, LDL-C, and ln(TG) levels in self-identified European American ( 20,000), African American ( 9,000), American Indian ( 6,000), Mexican American/Hispanic ( 2,500), Japanese/East Asian ( 690), and Pacific Islander/Native Hawaiian ( 175) adults, regardless of lipid-lowering medication use. We replicated 55 of 60 (92%) SNP associations tested in European Americans at p<0.05. Despite sufficient power, we were unable to replicate ABCA1 rs4149268 and rs1883025, CETP rs1864163, and TTC39B rs471364 previously associated with HDL-C and MAFB rs6102059 previously associated with LDL-C. Based on significance (p<0.05) and consistent direction of effect, a majority of replicated genotype-phentoype associations for HDL-C, LDL-C, and ln(TG) in European Americans generalized to African Americans (48%, 61%, and 57%), American Indians (45%, 64%, and 77%), and Mexican Americans/Hispanics (57%, 56%, and 86%). Overall, 16 associations generalized across all three populations. For the associations that did not generalize, differences in effect sizes, allele frequencies, and linkage disequilibrium offer clues to the next generation of association studies for these traits.

Our objective was to resequence insulin receptor substrate 2 (IRS2) to identify variants associated with obesity- and diabetes-related traits in Hispanic children. Exonic and intronic segments, 5' and 3' flanking regions of IRS2 (∼14.5 kb), were bidirectionally sequenced for single nucleotide polymorphism (SNP) discovery in 934 Hispanic children using 3730XL DNA Sequencers. Additionally, 15 SNPs derived from Illumina HumanOmni1-Quad BeadChips were analyzed. Measured genotype analysis tested associations between SNPs and obesity and diabetes-related traits. Bayesian quantitative trait nucleotide analysis was used to statistically infer the most likely functional polymorphisms. A total of 140 SNPs were identified with minor allele frequencies (MAF) ranging from 0.001 to 0.47. Forty-two of the 70 coding SNPs result in nonsynonymous amino acid substitutions relative to the consensus sequence; 28 SNPs were detected in the promoter, 12 in introns, 28 in the 3'-UTR, and 2 in the 5'-UTR. Two insertion/deletions (indels) were detected. Ten independent rare SNPs (MAF = 0.001-0.009) were associated with obesity-related traits (P = 0.01-0.00002). SNP 10510452_139 in the promoter region was shown to have a high posterior probability (P = 0.77-0.86) of influencing BMI, fat mass, and waist circumference in Hispanic children. SNP 10510452_139 contributed between 2 and 4% of the population variance in body weight and composition. None of the SNPs or indels were associated with diabetes-related traits or accounted for a previously identified quantitative trait locus on chromosome 13 for fasting serum glucose. Rare but not common IRS2 variants may play a role in the regulation of body weight but not an essential role in fasting glucose homeostasis in Hispanic children.

Genetic studies have identified thousands of variants associated with complex traits. However, most association studies are limited to populations of European descent and a single phenotype. The Population Architecture using Genomics and Epidemiology (PAGE) Study was initiated in 2008 by the National Human Genome Research Institute to investigate the epidemiologic architecture of well-replicated genetic variants associated with complex diseases in several large, ethnically diverse population-based studies. Combining DNA samples and hundreds of phenotypes from multiple cohorts, PAGE is well-suited to address generalization of associations and variability of effects in diverse populations; identify genetic and environmental modifiers; evaluate disease subtypes, intermediate phenotypes, and biomarkers; and investigate associations with novel phenotypes. PAGE investigators harmonize phenotypes across studies where possible and perform coordinated cohort-specific analyses and meta-analyses. PAGE researchers are genotyping thousands of genetic variants in up to 121,000 DNA samples from African-American, white, Hispanic/Latino, Asian/Pacific Islander, and American Indian participants. Initial analyses will focus on single nucleotide polymorphisms (SNPs) associated with obesity, lipids, cardiovascular disease, type 2 diabetes, inflammation, various cancers, and related biomarkers. PAGE SNPs are also assessed for pleiotropy using the "phenome-wide association study" approach, testing each SNP for associations with hundreds of phenotypes. PAGE data will be deposited into the National Center for Biotechnology Information's Database of Genotypes and Phenotypes and made available via a custom browser.

OBJECTIVES: The thickness of the brain's cortical gray matter (GM) and the fractional anisotropy (FA) of the cerebral white matter (WM) each follow an inverted U-shape trajectory with age. The two measures are positively correlated and may be modulated by common biological mechanisms. We employed four types of genetic analyses to localize individual genes acting pleiotropically upon these phenotypes.

METHODS: Whole-brain and regional GM thickness and FA values were measured from high-resolution anatomical and diffusion tensor MR images collected from 712, Mexican American participants (438 females, age = 47.9 ± 13.2 years) recruited from 73 (9.7 ± 9.3 individuals/family) large families. The significance of the correlation between two traits was estimated using a bivariate genetic correlation analysis. Localization of chromosomal regions that jointly influenced both traits was performed using whole-genome quantitative trait loci (QTL) analysis. Gene localization was performed using SNP genotyping on Illumina 1M chip and correlation with leukocyte-based gene-expression analyses. The gene-expressions were measured using the Illumina BeadChip. These data were available for 371 subjects.

RESULTS: Significant genetic correlation was observed among GM thickness and FA values. Significant logarithm of odds (LOD ≥ 3.0) QTLs were localized within chromosome 15q22-23. More detailed localization reported no significant association (p < 5·10(-5)) for 1565 SNPs located within the QTLs. Post hoc analysis indicated that 40% of the potentially significant (p ≤ 10(-3)) SNPs were localized to the related orphan receptor alpha (RORA) and NARG2 genes. A potentially significant association was observed for the rs2456930 polymorphism reported as a significant GWAS finding in Alzheimer's disease neuroimaging initiative subjects. The expression levels for RORA and ADAM10 genes were significantly (p < 0.05) correlated with both FA and GM thickness. NARG2 expressions were significantly correlated with GM thickness (p < 0.05) but failed to show a significant correlation (p = 0.09) with FA.

DISCUSSION: This study identified a novel, significant QTL at 15q22-23. SNP correlation with gene-expression analyses indicated that RORA, NARG2, and ADAM10 jointly influence GM thickness and WM-FA values.

BACKGROUND: Genome-wide association studies identified several single nucleotide polymorphisms (SNP) associated with prevalent coronary heart disease (CHD), but less is known of associations with incident CHD. The association of 13 published CHD SNPs was examined in 5 ancestry groups of 4 large US prospective cohorts.

METHODS AND RESULTS: The analyses included incident coronary events over an average 9.1 to 15.7 follow-up person-years in up to 26 617 white individuals (6626 events), 8018 black individuals (914 events), 1903 Hispanic individuals (113 events), 3669 American Indian individuals (595 events), and 885 Asian/Pacific Islander individuals (66 events). We used Cox proportional hazards models (with additive mode of inheritance) adjusted for age, sex, and ancestry (as needed). Nine loci were statistically associated with incident CHD events in white participants: 9p21 (rs10757278; P=4.7 × 10(-41)), 16q23.1 (rs2549513; P=0.0004), 6p24.1 (rs499818; P=0.0002), 2q36.3 (rs2943634; P=6.7 × 10(-6)), MTHFD1L (rs6922269, P=5.1 × 10(-10)), APOE (rs429358; P=2.7×10(-18)), ZNF627 (rs4804611; P=5.0 × 10(-8)), CXCL12 (rs501120; P=1.4 × 10(-6)) and LPL (rs268; P=2.7 × 10(-17)). The 9p21 region showed significant between-study heterogeneity, with larger effects in individuals age 55 years or younger and in women. Inclusion of coronary revascularization procedures among the incident CHD events introduced heterogeneity. The SNPs were not associated with CHD in black participants, and associations varied in other US minorities.

CONCLUSIONS: Prospective analyses of white participants replicated several reported cross-sectional CHD-SNP associations.

Although disrupted in schizophrenia 1 (DISC1) has been implicated in many psychiatric disorders, including schizophrenia, bipolar disorder, schizoaffective disorder and major depression, its biological role in these disorders is unclear. To better understand this gene and its role in psychiatric disease, we conducted transcriptional profiling and genome-wide association analysis in 1232 pedigreed Mexican-American individuals for whom we have neuroanatomic images, neurocognitive assessments and neuropsychiatric diagnoses. SOLAR was used to determine heritability, identify gene expression patterns and perform association analyses on 188 quantitative brain-related phenotypes. We found that the DISC1 transcript is highly heritable (h(2)=0.50; P=1.97 × 10(-22)), and that gene expression is strongly cis-regulated (cis-LOD=3.89) but is also influenced by trans-effects. We identified several DISC1 polymorphisms that were associated with cortical gray matter thickness within the parietal, temporal and frontal lobes. Associated regions affiliated with memory included the entorhinal cortex (rs821639, P=4.11 × 10(-5); rs2356606, P=4.71 × 10(-4)), cingulate cortex (rs16856322, P=2.88 × 10(-4)) and parahippocampal gyrus (rs821639, P=4.95 × 10(-4)); those affiliated with executive and other cognitive processing included the transverse temporal gyrus (rs9661837, P=5.21 × 10(-4); rs17773946, P=6.23 × 10(-4)), anterior cingulate cortex (rs2487453, P=4.79 × 10(-4); rs3738401, P=5.43 × 10(-4)) and medial orbitofrontal cortex (rs9661837; P=7.40 × 10(-4)). Cognitive measures of working memory (rs2793094, P=3.38 × 10(-4)), as well as lifetime history of depression (rs4658966, P=4.33 × 10(-4); rs12137417, P=4.93 × 10(-4)) and panic (rs12137417, P=7.41 × 10(-4)) were associated with DISC1 sequence variation. DISC1 has well-defined genetic regulation and clearly influences important phenotypes related to psychiatric disease.

Genes play an important role in lifelong skeletal health. Genes that influence bone building during childhood have the potential to affect bone health not only throughout childhood but also into adulthood. Given that peak bone mass is a significant predictor of adult fracture risk, it is imperative that the genetic underpinnings of the normal pediatric skeleton are uncovered. In a sample of 600 10-year-old children from 144 families in the Fels Longitudinal Study, we examined radiographic cortical bone measures of the second metacarpal. Morphometic measurements included bone width, medial and lateral cortical thicknesses, and the calculated cortical index representing the amount of cortex relative to bone width. We then conducted genome-wide linkage analysis on these traits in 440 genotyped individuals using the SOLAR analytic platform. Significant quantitative trait loci (QTL) were identified for bone traits on three separate chromosomes. A QTL for medial cortical thickness was localized to chromosome 2p25.2. A QTL for lateral cortical thickness was localized to chromosomal region 3p26.1-3p25.3. Finally, a QTL detected for cortical index was localized to the 17q21.2 chromosomal region. Each region contains plausible candidate genes for pediatric skeletal health, some of which confirm findings from studies of adulthood bone, and for others represent novel candidate genes for skeletal health.

OBJECTIVE: Genome-wide association studies that compare the statistical association between thousands of DNA variations and a human trait have detected 958 loci across 127 different diseases and traits. However, these statistical associations only provide evidence for genomic regions likely to harbor a causal gene(s) and do not directly identify such genes. We combined gene variation and expression data in a human cohort to identify causal genes.

RESEARCH DESIGN AND METHODS: Global gene transcription activity was obtained for each individual in a large human cohort (n = 1,240). These quantitative transcript data were tested for correlation with genotype data generated from the same individuals to identify gene expression patterns influenced by the variants.

RESULTS: Variant rs8050136 lies within intron 1 of the FTO gene on chromosome 16 and marks a locus strongly associated with type 2 diabetes and obesity and widely replicated across many populations. We report that genetic variation at this locus does not influence FTO gene expression levels (P = 0.38), but is strongly correlated with expression of RBL2 (P = 2.7 x 10(-5)), approximately 270,000 base pairs distant to FTO.

CONCLUSIONS: These data suggest that variants at FTO influence RBL2 gene expression at large genetic distances. This observation underscores the complexity of human transcriptional regulation and highlights the utility of large human cohorts in which both genetic variation and global gene expression data are available to identify disease genes. Expedient identification of genes mediating the effects of genome-wide association study-identified loci will enable mechanism-of-action studies and accelerate understanding of human disease processes under genetic influence.