ABSTRACT
Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.
Subject(s)
Diabetes Mellitus, Type 2 , Disease Progression , Genetic Predisposition to Disease , Genome-Wide Association Study , Humans , Adipocytes/metabolism , Chromatin/genetics , Chromatin/metabolism , Coronary Artery Disease/complications , Coronary Artery Disease/genetics , Diabetes Mellitus, Type 2/classification , Diabetes Mellitus, Type 2/complications , Diabetes Mellitus, Type 2/genetics , Diabetes Mellitus, Type 2/pathology , Diabetes Mellitus, Type 2/physiopathology , Diabetic Nephropathies/complications , Diabetic Nephropathies/genetics , Endothelial Cells/metabolism , Enteroendocrine Cells , Epigenomics , Genetic Predisposition to Disease/genetics , Islets of Langerhans/metabolism , Multifactorial Inheritance/genetics , Peripheral Arterial Disease/complications , Peripheral Arterial Disease/genetics , Single-Cell AnalysisABSTRACT
Critical illness in COVID-19 is an extreme and clinically homogeneous disease phenotype that we have previously shown1 to be highly efficient for discovery of genetic associations2. Despite the advanced stage of illness at presentation, we have shown that host genetics in patients who are critically ill with COVID-19 can identify immunomodulatory therapies with strong beneficial effects in this group3. Here we analyse 24,202 cases of COVID-19 with critical illness comprising a combination of microarray genotype and whole-genome sequencing data from cases of critical illness in the international GenOMICC (11,440 cases) study, combined with other studies recruiting hospitalized patients with a strong focus on severe and critical disease: ISARIC4C (676 cases) and the SCOURGE consortium (5,934 cases). To put these results in the context of existing work, we conduct a meta-analysis of the new GenOMICC genome-wide association study (GWAS) results with previously published data. We find 49 genome-wide significant associations, of which 16 have not been reported previously. To investigate the therapeutic implications of these findings, we infer the structural consequences of protein-coding variants, and combine our GWAS results with gene expression data using a monocyte transcriptome-wide association study (TWAS) model, as well as gene and protein expression using Mendelian randomization. We identify potentially druggable targets in multiple systems, including inflammatory signalling (JAK1), monocyte-macrophage activation and endothelial permeability (PDE4A), immunometabolism (SLC2A5 and AK5), and host factors required for viral entry and replication (TMPRSS2 and RAB2A).
Subject(s)
COVID-19 , Critical Illness , Genetic Predisposition to Disease , Genetic Variation , Genome-Wide Association Study , Humans , COVID-19/genetics , Genetic Predisposition to Disease/genetics , Genetic Variation/genetics , Genotype , Genotyping Techniques , Monocytes/metabolism , Phenotype , rab GTP-Binding Proteins/genetics , Transcriptome , Whole Genome SequencingABSTRACT
Interpretation of genetic association results is difficult because signals often lack biological context. To generate hypotheses of the functional genetic etiology of complex cardiometabolic traits, we estimated the genetically determined component of gene expression from common variants using PrediXcan (1) and determined genes with differential predicted expression by trait. PrediXcan imputes tissue-specific expression levels from genetic variation using variant-level effect on gene expression in transcriptome data. To explore the value of imputed genetically regulated gene expression (GReX) models across different ancestral populations, we evaluated imputed expression levels for predictive accuracy genome-wide in RNA sequence data in samples drawn from European-ancestry and African-ancestry populations and identified substantial predictive power using European-derived models in a non-European target population. We then tested the association of GReX on 15 cardiometabolic traits including blood lipid levels, body mass index, height, blood pressure, fasting glucose and insulin, RR interval, fibrinogen level, factor VII level and white blood cell and platelet counts in 15 755 individuals across three ancestry groups, resulting in 20 novel gene-phenotype associations reaching experiment-wide significance across ancestries. In addition, we identified 18 significant novel gene-phenotype associations in our ancestry-specific analyses. Top associations were assessed for additional support via query of S-PrediXcan (2) results derived from publicly available genome-wide association studies summary data. Collectively, these findings illustrate the utility of transcriptome-based imputation models for discovery of cardiometabolic effect genes in a diverse dataset.
Subject(s)
Forecasting/methods , Metabolome/genetics , Metabolome/physiology , Adult , Aged , Blood Pressure , Body Mass Index , Chromosome Mapping/methods , Ethnicity/genetics , Female , Genetic Association Studies/methods , Genome-Wide Association Study/methods , Humans , Male , Middle Aged , Multifactorial Inheritance/genetics , Phenotype , Polymorphism, Single Nucleotide/genetics , Transcriptome/genetics , White People/geneticsABSTRACT
Although hundreds of genome-wide association studies-implicated loci have been reported for adult obesity-related traits, less is known about the genetics specific for early-onset obesity and with only a few studies conducted in non-European populations to date. Searching for additional genetic variants associated with childhood obesity, we performed a trans-ancestral meta-analysis of 30 studies consisting of up to 13 005 cases (≥95th percentile of body mass index (BMI) achieved 2-18 years old) and 15 599 controls (consistently <50th percentile of BMI) of European, African, North/South American and East Asian ancestry. Suggestive loci were taken forward for replication in a sample of 1888 cases and 4689 controls from seven cohorts of European and North/South American ancestry. In addition to observing 18 previously implicated BMI or obesity loci, for both early and late onset, we uncovered one completely novel locus in this trans-ancestral analysis (nearest gene, METTL15). The variant was nominally associated with only the European subgroup analysis but had a consistent direction of effect in other ethnicities. We then utilized trans-ancestral Bayesian analysis to narrow down the location of the probable causal variant at each genome-wide significant signal. Of all the fine-mapped loci, we were able to narrow down the causative variant at four known loci to fewer than 10 single nucleotide polymorphisms (SNPs) (FAIM2, GNPDA2, MC4R and SEC16B loci). In conclusion, an ethnically diverse setting has enabled us to both identify an additional pediatric obesity locus and further fine-map existing loci.
Subject(s)
Chromosome Mapping/methods , Genome-Wide Association Study/methods , Pediatric Obesity/genetics , Polymorphism, Single Nucleotide , Wilms Tumor/genetics , Bayes Theorem , Case-Control Studies , Child , Female , Genetic Loci , Genetic Predisposition to Disease , Humans , MaleABSTRACT
BACKGROUND/OBJECTIVES: Mexico has one of the highest prevalence of childhood obesity in the world. Genome-wide association studies (GWAS) for obesity have identified multiple single-nucleotide polymorphisms (SNPs) in populations of European, East Asian, and African descent. The contribution of these loci to obesity in Mexican children is unclear. We assessed the transferability of 98 obesity loci in Mexican children and fine-mapped the association signals. SUBJECTS/METHODS: The study included 405 and 390 Mexican children with normal weight and obesity. Participants were genotyped with a genome-wide dense SNP array designed for Latino populations, allowing for the analysis of GWAS index SNPs as well as fine-mapping SNPs, totaling 750 SNPs covering 98 loci. Two genetic risk scores (GRS) were constructed: a "discovery GRS" and a "best-associated GRS", representing the number of effect alleles at the GWAS index SNPs and at the best-associated SNPs after fine-mapping for each subject. RESULTS: Seventeen obesity loci were significantly associated with obesity, and five had fine-mapping SNPs significantly better associated with obesity than their corresponding GWAS index SNPs in Mexican children. Six obesity-associated SNPs significantly departed from additive to dominant (N = 5) or recessive (N = 1) models, and a significant interaction was found between rs274609 (TNNI3K) and rs1010553 (ITIH4) on childhood obesity risk. The best-associated GRS was significantly more associated with childhood obesity (OR = 1.21 per additional risk allele [95%CI:1.17-1.25], P = 4.8 × 10-25) than the discovery GRS (OR = 1.05 per additional risk allele [95%CI:1.02-1.08], P = 8.0 × 10-4), and was also associated with waist-to-hip ratio, fasting glucose, fasting insulin and triglyceride levels, the association being mediated by obesity. An overall depletion of obesity risk alleles was observed in Mexican children with normal weight when compared to GWAS discovery populations. CONCLUSIONS: Our study indicates a partial transferability of GWAS obesity loci in Mexican children, and supports the pertinence of post-GWAS fine-mapping experiments in the admixed Mexican population.
Subject(s)
Genetic Predisposition to Disease/epidemiology , Genome-Wide Association Study , Obesity/genetics , Polymorphism, Single Nucleotide/genetics , Adolescent , Child , Fasting/blood , Female , Gene Frequency , Genetic Loci , Humans , Life Style , Male , Mexico/epidemiology , Obesity/blood , Obesity/epidemiology , Socioeconomic FactorsABSTRACT
BACKGROUND: Association studies in recently admixed populations are extremely useful to identify the genetic architecture of pigmentation, due to their high genotypic and phenotypic variation. However, to date only four Genome-Wide Association Studies (GWAS) have been carried out in these populations. RESULTS: We present a GWAS of skin pigmentation in an admixed sample from Cuba (N = 762). Additionally, we conducted a meta-analysis including the Cuban sample, and admixed samples from Cape Verde, Puerto Rico and African-Americans from San Francisco. This meta-analysis is one of the largest efforts so far to characterize the genetic basis of skin pigmentation in admixed populations (N = 2,104). We identified five genome-wide significant regions in the meta-analysis, and explored if the markers observed in these regions are associated with the expression of relevant pigmentary genes in human melanocyte cultures. In three of the regions identified in the meta-analysis (SLC24A5, SLC45A2, and GRM5/TYR), the association seems to be driven by non-synonymous variants (rs1426654, rs16891982, and rs1042602, respectively). The rs16891982 polymorphism is strongly associated with the expression of the SLC45A2 gene. In the GRM5/TYR region, in addition to the rs1042602 non-synonymous SNP located on the TYR gene, variants located in the nearby GRM5 gene have an independent effect on pigmentation, possibly through regulation of gene expression of the TYR gene. We also replicated an association recently described near the MFSD12 gene on chromosome 19 (lead variant rs112332856). Additionally, our analyses support the presence of multiple signals in the OCA2/HERC2/APBA2 region on chromosome 15. A clear causal candidate is the HERC2 intronic variant rs12913832, which has a profound influence on OCA2 expression. This variant has pleiotropic effects on eye, hair, and skin pigmentation. However, conditional and haplotype-based analyses indicate the presence of other variants with independent effects on melanin levels in OCA2 and APBA2. Finally, a follow-up of genome-wide signals identified in a recent GWAS for tanning response indicates that there is a substantial overlap in the genetic factors influencing skin pigmentation and tanning response. CONCLUSIONS: Our meta-analysis of skin pigmentation GWAS in recently admixed populations provides new insights about the genetic architecture of this complex trait.
Subject(s)
Genetics, Population , Genome-Wide Association Study , Skin Pigmentation/genetics , Alleles , Genotype , Humans , Linkage Disequilibrium , Polymorphism, Single Nucleotide , Quantitative Trait LociABSTRACT
Like many highly variable human traits, more than a dozen genes are known to contribute to the full range of skin color. However, the historical bias in favor of genetic studies in European and European-derived populations has blinded us to the magnitude of pigmentation's complexity. As deliberate efforts are being made to better characterize diverse global populations and new sequencing technologies, better measurement tools, functional assessments, predictive modeling, and ancient DNA analyses become more widely accessible, we are beginning to appreciate how limited our understanding of the genetic bases of human skin color have been. Novel variants in genes not previously linked to pigmentation have been identified and evidence is mounting that there are hundreds more variants yet to be found. Even for genes that have been exhaustively characterized in European populations like MC1R, OCA2, and SLC24A5, research in previously understudied groups is leading to a new appreciation of the degree to which genetic diversity, epistatic interactions, pleiotropy, admixture, global and local adaptation, and cultural practices operate in population-specific ways to shape the genetic architecture of skin color. Furthermore, we are coming to terms with how factors like tanning response and barrier function may also have influenced selection on skin throughout human history. By examining how our knowledge of pigmentation genetics has shifted in the last decade, we can better appreciate how far we have come in understanding human diversity and the still long road ahead for understanding many complex human traits.
Subject(s)
Biological Evolution , Skin Physiological Phenomena , Skin Pigmentation , Anthropology, Physical , Antiporters/genetics , Genetics, Population , Guanine Nucleotide Exchange Factors/genetics , Humans , Membrane Proteins/genetics , Membrane Transport Proteins/genetics , Receptor, Melanocortin, Type 1/genetics , Skin/metabolism , Skin Pigmentation/genetics , Skin Pigmentation/physiology , Ubiquitin-Protein LigasesABSTRACT
The evidence for a relationship between serum vitamin D levels and nonskeletal health outcomes is inconsistent. The validity of single or predicted measurements of 25-hydroxyvitamin D (25(OH)D) concentration is unknown, as levels of this biomarker are highly seasonally variable. We compared models of 25(OH)D measured at baseline, at multiple time points throughout the year, and averaged over the year among 309 persons in Toronto, Ontario, Canada (43°N latitude) during 2009-2013. Information and blood samples were collected every 2 months. Baseline and average 25(OH)D concentrations were correlated (r = 0.88). Major factors associated with 25(OH)D level were similar across models and included race/ethnicity (concentrations in non-European groups were lower than those in Europeans), vitamin D supplement use of ≥1,000 IU/day (18.9 nmol/L (95% confidence interval (CI): 16.1, 21.8) vs. no supplement use in a full data set with all factors), and the presence of the group-specific component/vitamin D binding protein gene (GC/DBP) rs4588 functional polymorphism (AA vs. CC: -16.7 nmol/L (95% CI: -26.2, -7.1); CA vs. CC: -10.7 nmol/L (95% CI: -14.9, -6.5)). Most factors had similar associations in Europeans and non-Europeans. Genetic factors may play a greater role in average 25(OH)D concentrations. Prediction models for 25(OH)D are challenging and population-specific, but use of genetic factors along with a few common population-relevant, quantifiable nongenetic factors with strong associations may be the most feasible approach to vitamin D assessment over time.
Subject(s)
Skin Pigmentation/genetics , Skin Pigmentation/physiology , Vitamin D/analogs & derivatives , Adolescent , Adult , Asian People/genetics , Black People/genetics , Body Mass Index , Female , Genetic Markers , Humans , Linear Models , Male , Middle Aged , Ontario/epidemiology , Polymorphism, Single Nucleotide , Seasons , Vitamin D/blood , Vitamin D/genetics , White People/genetics , Young AdultABSTRACT
Top signals from genome-wide association studies (GWASs) of type 2 diabetes (T2D) are enriched with expression quantitative trait loci (eQTLs) identified in skeletal muscle and adipose tissue. We therefore hypothesized that such eQTLs might account for a disproportionate share of the heritability estimated from all SNPs interrogated through GWASs. To test this hypothesis, we applied linear mixed models to the Wellcome Trust Case Control Consortium (WTCCC) T2D data set and to data sets representing Mexican Americans from Starr County, TX, and Mexicans from Mexico City. We estimated the proportion of phenotypic variance attributable to the additive effect of all variants interrogated in these GWASs, as well as a much smaller set of variants identified as eQTLs in human adipose tissue, skeletal muscle, and lymphoblastoid cell lines. The narrow-sense heritability explained by all interrogated SNPs in each of these data sets was substantially greater than the heritability accounted for by genome-wide-significant SNPs (â¼10%); GWAS SNPs explained over 50% of phenotypic variance in the WTCCC, Starr County, and Mexico City data sets. The estimate of heritability attributable to cross-tissue eQTLs was greater in the WTCCC data set and among lean Hispanics, whereas adipose eQTLs significantly explained heritability among Hispanics with a body mass index ≥ 30. These results support an important role for regulatory variants in the genetic component of T2D susceptibility, particularly for eQTLs that elicit effects across insulin-responsive peripheral tissues.
Subject(s)
Diabetes Mellitus, Type 2/genetics , Phenotype , Quantitative Trait Loci/genetics , Adipose Tissue/chemistry , Analysis of Variance , Body Mass Index , Genome-Wide Association Study/statistics & numerical data , Humans , Linear Models , Mexican Americans/genetics , Mexico , Muscle, Skeletal/chemistry , Polymorphism, Single Nucleotide/genetics , Principal Component Analysis , TexasABSTRACT
MOTIVATION: Rapid advances in genotyping and genome-wide association studies have enabled the discovery of many new genotype-phenotype associations at the resolution of individual markers. However, these associations explain only a small proportion of theoretically estimated heritability of most diseases. In this work, we propose an integrative mixture model called JBASE: joint Bayesian analysis of subphenotypes and epistasis. JBASE explores two major reasons of missing heritability: interactions between genetic variants, a phenomenon known as epistasis and phenotypic heterogeneity, addressed via subphenotyping. RESULTS: Our extensive simulations in a wide range of scenarios repeatedly demonstrate that JBASE can identify true underlying subphenotypes, including their associated variants and their interactions, with high precision. In the presence of phenotypic heterogeneity, JBASE has higher Power and lower Type 1 Error than five state-of-the-art approaches. We applied our method to a sample of individuals from Mexico with Type 2 diabetes and discovered two novel epistatic modules, including two loci each, that define two subphenotypes characterized by differences in body mass index and waist-to-hip ratio. We successfully replicated these subphenotypes and epistatic modules in an independent dataset from Mexico genotyped with a different platform. AVAILABILITY AND IMPLEMENTATION: JBASE is implemented in C++, supported on Linux and is available at http://www.cs.toronto.edu/â¼goldenberg/JBASE/jbase.tar.gz. The genotype data underlying this study are available upon approval by the ethics review board of the Medical Centre Siglo XXI. Please contact Dr Miguel Cruz at mcruzl@yahoo.com for assistance with the application. CONTACT: anna.goldenberg@utoronto.ca SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Algorithms , Epistasis, Genetic , Phenotype , Bayes Theorem , Body Mass Index , Diabetes Mellitus, Type 2/genetics , Genome-Wide Association Study , Genotype , Genotyping Techniques , Humans , Mexico , Waist-Hip RatioABSTRACT
We carried out an admixture analysis of a sample comprising 1,019 individuals from all the provinces of Cuba. We used a panel of 128 autosomal Ancestry Informative Markers (AIMs) to estimate the admixture proportions. We also characterized a number of haplogroup diagnostic markers in the mtDNA and Y-chromosome in order to evaluate admixture using uniparental markers. Finally, we analyzed the association of 16 single nucleotide polymorphisms (SNPs) with quantitative estimates of skin pigmentation. In the total sample, the average European, African and Native American contributions as estimated from autosomal AIMs were 72%, 20% and 8%, respectively. The Eastern provinces of Cuba showed relatively higher African and Native American contributions than the Western provinces. In particular, the highest proportion of African ancestry was observed in the provinces of Guantánamo (40%) and Santiago de Cuba (39%), and the highest proportion of Native American ancestry in Granma (15%), Holguín (12%) and Las Tunas (12%). We found evidence of substantial population stratification in the current Cuban population, emphasizing the need to control for the effects of population stratification in association studies including individuals from Cuba. The results of the analyses of uniparental markers were concordant with those observed in the autosomes. These geographic patterns in admixture proportions are fully consistent with historical and archaeological information. Additionally, we identified a sex-biased pattern in the process of gene flow, with a substantially higher European contribution from the paternal side, and higher Native American and African contributions from the maternal side. This sex-biased contribution was particularly evident for Native American ancestry. Finally, we observed that SNPs located in the genes SLC24A5 and SLC45A2 are strongly associated with melanin levels in the sample.
Subject(s)
Gene Flow/genetics , Genetics, Population , Haplotypes/genetics , Pigmentation/genetics , Black People/genetics , Chromosomes, Human, Y/genetics , Cuba , DNA, Mitochondrial/genetics , Hispanic or Latino/genetics , Humans , Indians, North American/genetics , Polymorphism, Single Nucleotide/genetics , White People/geneticsABSTRACT
OBJECTIVE: Several algorithms have been proposed to reduce the genotyping effort and cost, while retaining the accuracy of N-acetyltransferase-2 (NAT2) phenotype prediction. Data from the 1000 Genomes (1KG) project and an admixed cohort of Black Brazilians were used to assess the accuracy of NAT2 phenotype prediction using algorithms based on paired single nucleotide polymorphisms (SNPs) (rs1041983 and rs1801280) or a tag SNP (rs1495741). METHODS: NAT2 haplotypes comprising SNPs rs1801279, rs1041983, rs1801280, rs1799929, rs1799930, rs1208 and rs1799931 were assigned according to the arylamine N-acetyltransferases database. Contingency tables were used to visualize the agreement between the NAT2 acetylator phenotypes on the basis of these haplotypes versus phenotypes inferred by the prediction algorithms. RESULTS: The paired and tag SNP algorithms provided more than 96% agreement with the 7-SNP derived phenotypes in Europeans, East Asians, South Asians and Admixed Americans, but discordance of phenotype prediction occurred in 30.2 and 24.8% 1KG Africans and in 14.4 and 18.6% Black Brazilians, respectively. Paired SNP panel misclassification occurs in carriers of NATs haplotypes *13A (282T alone), *12B (282T and 803G), *6B (590A alone) and *14A (191A alone), whereas haplotype *14, defined by the 191A allele, is the major culprit of misclassification by the tag allele. CONCLUSION: Both the paired SNP and the tag SNP algorithms may be used, with economy of scale, to infer NAT2 acetylator phenotypes, including the ultra-slow phenotype, in European, East Asian, South Asian and American populations represented in the 1KG cohort. Both algorithms, however, perform poorly in populations of predominant African descent, including admixed African-Americans, African Caribbeans and Black Brazilians.
Subject(s)
Algorithms , Arylamine N-Acetyltransferase/genetics , Ethnicity/genetics , Polymorphism, Single Nucleotide/genetics , Racial Groups/genetics , Acetylation , Genome-Wide Association Study , Humans , PhenotypeABSTRACT
Although disproportionately affected by increasing rates of type 2 diabetes and dyslipidemias, Hispanic populations are underrepresented in efforts to understand genetic susceptibility to these disorders. Where research has been undertaken, these populations have provided substantial insight into identification of novel risk-associated genes and have aided in the ability to fine map previously described risk loci. Genome-wide analyses in Hispanic and trans-ethnic populations have resulted in identification of more than 40 replicated or novel genes with significant effects for type 2 diabetes or lipid traits. Initial investigations into rare variant effects have identified new risk-associated variants private to Hispanic populations, and preliminary results suggest metagenomic approaches in Hispanic populations, such as characterizing the gut microbiome, will enable the development of new predictive tools and therapeutic targets for type 2 diabetes. Future genome-wide studies in expanded cohorts of Hispanics are likely to result in new insights into the genetic etiology of metabolic health.
Subject(s)
Diabetes Mellitus, Type 2/genetics , Genome-Wide Association Study , Lipids/genetics , Genetic Loci , Genetic Predisposition to Disease , Hispanic or Latino , Humans , Metagenome , Risk FactorsABSTRACT
OBJECTIVES: The main goals of this study are to 1) quantitatively measure skin, hair, and iris pigmentation in a diverse sample of individuals, 2) describe variation within and between these samples, and 3) demonstrate how quantitative measures can facilitate genotype-phenotype association tests. MATERIALS AND METHODS: We quantitatively characterize skin, hair, and iris pigmentation using the Melanin (M) Index (skin) and CIELab values (hair) in 1,450 individuals who self-identify as African American, East Asian, European, Hispanic, or South Asian. We also quantify iris pigmentation in a subset of these individuals using CIELab values from high-resolution iris photographs. We compare mean skin M index and hair and iris CIELab values among populations using ANOVA and MANOVA respectively and test for genotype-phenotype associations in the European sample. RESULTS: All five populations are significantly different for skin (P <2 × 10(-16) ) and hair color (P <2 × 10(-16) ). Our quantitative analysis of iris and hair pigmentation reinforces the continuous, rather than discrete, nature of these traits. We confirm the association of three loci (rs16891982, rs12203592, and rs12913832) with skin pigmentation and four loci (rs12913832, rs12203592, rs12896399, and rs16891982) with hair pigmentation. Interestingly, the derived rs12203592 T allele located within the IRF4 gene is associated with lighter skin but darker hair color. DISCUSSION: The quantitative methods used here provide a fine-scale assessment of pigmentation phenotype and facilitate genotype-phenotype associations, even with relatively small sample sizes. This represents an important expansion of current investigations into pigmentation phenotype and associated genetic variation by including non-European and admixed populations. Am J Phys Anthropol 160:570-581, 2016. © 2015 Wiley Periodicals, Inc.
Subject(s)
Eye Color/genetics , Hair Color/genetics , Racial Groups/genetics , Skin Pigmentation/genetics , Analysis of Variance , Anthropology, Physical , Humans , Linear Models , Polymorphism, Single Nucleotide/geneticsABSTRACT
OBJECTIVES: The main goal of the study was to test the association of 18 polymorphisms located within nine pigmentation candidate genes with quantitative skin pigmentation measures collected in a sample of individuals of East Asian ancestry living in Canada (N = 419). METHODS: The 18 polymorphisms are located within genes that show putative signatures of positive selection in East Asian populations. The genetic markers were selected for genotyping based on an annotation of common East Asian polymorphisms to predict potential functional effects. We restricted our attention to polymorphisms that have an allele frequency difference of at least 30% between East Asian populations and African and European populations, or have alleles that are present in East Asians, but are absent in Africans and Europeans. RESULTS: Two nonsynonymous variants selected within the OCA2 gene, rs1800414 (His615Arg) and rs74653330 (Ala481Thr), were significantly associated with melanin levels in the sample. Both single nucleotide polymorphisms (SNPs) are nonsynonymous polymorphisms located more than 30 kb apart on chromosome 15 and have very different frequencies in the East Asian sample. Additionally, both polymorphisms are predicted to have a deleterious effect on the protein. Linear regression analysis using an additive model indicate that each copy of the derived rs1800414 allele G decreases Melanin Index approximately 0.9 units and each copy of the derived rs74653330 allele A decreases Melanin Index approximately 1.9 units. CONCLUSIONS: Two nonsynonymous OCA2 polymorphisms (rs1800414 and rs74653330) are independently associated with normal skin pigmentation variation in East Asian populations and have very different frequency distributions in East Asia.
Subject(s)
Albinism, Oculocutaneous/epidemiology , Membrane Transport Proteins/genetics , Polymorphism, Single Nucleotide , Skin Pigmentation , Adult , Albinism, Oculocutaneous/genetics , Asia, Eastern/ethnology , Female , Gene Frequency , Genotype , Humans , Male , Membrane Transport Proteins/metabolism , Ontario/epidemiology , Vietnam/ethnology , Young AdultABSTRACT
The inverse correlation between skin pigmentation and latitude observed in human populations is thought to have been shaped by selective pressures favoring lighter skin to facilitate vitamin D synthesis in regions far from the equator. Several candidate genes for skin pigmentation have been shown to exhibit patterns of polymorphism that overlap the geospatial variation in skin color. However, little work has focused on estimating the time frame over which skin pigmentation has changed and on the intensity of selection acting on different pigmentation genes. To provide a temporal framework for the evolution of lighter pigmentation, we used forward Monte Carlo simulations coupled with a rejection sampling algorithm to estimate the time of onset of selective sweeps and selection coefficients at four genes associated with this trait in Europeans: KITLG, TYRP1, SLC24A5, and SLC45A2. Using compound haplotype systems consisting of rapidly evolving microsatellites linked to one single-nucleotide polymorphism in each gene, we estimate that the onset of the sweep shared by Europeans and East Asians at KITLG occurred approximately 30,000 years ago, after the out-of-Africa migration, whereas the selective sweeps for the European-specific alleles at TYRP1, SLC24A5, and SLC45A2 started much later, within the last 11,000-19,000 years, well after the first migrations of modern humans into Europe. We suggest that these patterns were influenced by recent increases in size of human populations, which favored the accumulation of advantageous variants at different loci.
Subject(s)
Evolution, Molecular , Skin Pigmentation/genetics , White People/genetics , Africa , Alleles , Antigens, Neoplasm/genetics , Antigens, Neoplasm/metabolism , Antiporters/genetics , Antiporters/metabolism , Asian People/genetics , DNA/genetics , DNA/isolation & purification , Ethnicity/genetics , Europe , Genetic Loci , Haplotypes , Humans , Membrane Glycoproteins/genetics , Membrane Glycoproteins/metabolism , Membrane Transport Proteins/genetics , Membrane Transport Proteins/metabolism , Microsatellite Repeats , Oxidoreductases/genetics , Oxidoreductases/metabolism , Phenotype , Phylogeny , Polymorphism, Single Nucleotide , Selection, GeneticABSTRACT
BACKGROUND: Currently, there is very limited knowledge about the genes involved in normal pigmentation variation in East Asian populations. We carried out a genome-wide scan of signatures of positive selection using the 1000 Genomes Phase I dataset, in order to identify pigmentation genes showing putative signatures of selective sweeps in East Asia. We applied a broad range of methods to detect signatures of selection including: 1) Tests designed to identify deviations of the Site Frequency Spectrum (SFS) from neutral expectations (Tajima's D, Fay and Wu's H and Fu and Li's D* and F*), 2) Tests focused on the identification of high-frequency haplotypes with extended linkage disequilibrium (iHS and Rsb) and 3) Tests based on genetic differentiation between populations (LSBL). Based on the results obtained from a genome wide analysis of 25 kb windows, we constructed an empirical distribution for each statistic across all windows, and identified pigmentation genes that are outliers in the distribution. RESULTS: Our tests identified twenty genes that are relevant for pigmentation biology. Of these, eight genes (ATRN, EDAR, KLHL7, MITF, OCA2, TH, TMEM33 and TRPM1,) were extreme outliers (top 0.1% of the empirical distribution) for at least one statistic, and twelve genes (ADAM17, BNC2, CTSD, DCT, EGFR, LYST, MC1R, MLPH, OPRM1, PDIA6, PMEL (SILV) and TYRP1) were in the top 1% of the empirical distribution for at least one statistic. Additionally, eight of these genes (BNC2, EGFR, LYST, MC1R, OCA2, OPRM1, PMEL (SILV) and TYRP1) have been associated with pigmentary traits in association studies. CONCLUSIONS: We identified a number of putative pigmentation genes showing extremely unusual patterns of genetic variation in East Asia. Most of these genes are outliers for different tests and/or different populations, and have already been described in previous scans for positive selection, providing strong support to the hypothesis that recent selective sweeps left a signature in these regions. However, it will be necessary to carry out association and functional studies to demonstrate the implication of these genes in normal pigmentation variation.
Subject(s)
Asian People/genetics , Proteins/genetics , Selection, Genetic , Skin Pigmentation/genetics , Asian People/ethnology , China/ethnology , Genetic Variation , Genetics, Population , Haplotypes , Humans , Japan/ethnology , Phenotype , Polymorphism, Single NucleotideABSTRACT
BACKGROUND: Several studies in type 2 diabetes patients have shown significant associations between the SOD2 gene Val16Ala polymorphism and albuminuria, but this association has not been explored in the Mexican population. METHODS: We evaluated the association between the SOD2 gene Val16Ala polymorphism (rs4880) and macroalbuminuria in a sample of 994 unrelated Mexican type 2 diabetes patients. The study included 119 subjects with urinary albumin >300 mg/dL and 875 subjects with urinary albumin ≤ 30 mg/dL. Genotyping of the SOD2 gene Val16Ala SNP was carried out with Real-Time Polymerase Chain Reaction (RT-PCR). RESULTS: The frequency of the TT genotype was 6.7% higher in participants with macroalbuminuria than in the normoalbuminuria group (16.8% vs. 10.1%). Using a logistic regression analysis, we observed that individuals with the CC genotype had significantly lower risks of macroalbuminuria than those with the TT genotype (OR=0.42, p=0.034). We carried out a meta-analysis combining our data with data from four previous studies and estimated an odds ratio (95% CI) for the C allele (with respect to the reference T allele) of 0.65 (0.52-0.80, p<0.001). CONCLUSIONS: A significant association was found between the SOD2 Val16Ala polymorphism and macroalbuminuria in a sample of Mexican type 2 diabetes patients.
Subject(s)
Albuminuria/genetics , Diabetes Mellitus, Type 2/genetics , Superoxide Dismutase/genetics , Aged , Albuminuria/complications , Alleles , Diabetes Mellitus, Type 2/complications , Female , Gene Frequency , Genotype , Humans , Logistic Models , Male , Mexico/ethnology , Middle Aged , Odds Ratio , Polymorphism, Single NucleotideABSTRACT
BACKGROUND: Physiological and lifestyle factors put older adults at an increased risk of vitamin D insufficiency and resulting negative health outcomes. Here we explore the vitamin D status in a sample of community dwelling older adults of diverse ancestry living in the Greater Toronto area (GTA). METHODS: Two hundred and twenty-four (224) adults over 60 years of age were recruited from the Square One Older Adult Centre, in Mississauga, Ontario. Circulating 25-hydroxyvitamin D (25(OH)D) concentrations were measured from dried blood spot cards. Dietary and supplemental intakes of vitamin D were assessed via questionnaires. Skin pigmentation was assessed quantitatively by measuring melanin levels using a reflectometer. RESULTS: The mean 25(OH)D concentration in the total sample was 82.4 nmol/L. There were no statistically significant differences in serum 25(OH)D concentrations, supplemental or dietary vitamin D intakes between the three major ancestral groups (East Asians, Europeans and South Asians). Females had significantly higher 25(OH)D concentrations than males (84.5 nmol/L vs. 72.2 nmol/L, p = 0.012). The proportion of participants with 25(OH)D concentrations below 50 nmol/L and 75 nmol/L were 12.1%, and 38.8%, respectively. The mean daily supplemental intake of vitamin D was 917 IU/day. Vitamin D intake from supplements was the major factor determining 25(OH)D concentrations (p < 0.001). CONCLUSIONS: Mean concentration of 25(OH)D in a sample of older adults of diverse ancestry living in the GTA exceeded 80 nmol/L, and there were no significant differences in 25(OH)D levels between ancestral groups. These results sharply contrast with our recent study focused on young adults of diverse ancestry living in the same geographic area, in which we found substantially lower 25(OH)D concentrations (mean 39.5 nmol/L), low supplemental vitamin D intake (114 IU/day), and significant differences in 25(OH)D levels between ancestral groups. High daily intake of supplemental vitamin D in this sample of older adults likely accounts for such disparate findings with respect to the young adult sample.