1
|
Chen Y, Branicki W, Walsh S, Nothnagel M, Kayser M, Liu F. The impact of correlations between pigmentation phenotypes and underlying genotypes on genetic prediction of pigmentation traits. Forensic Sci Int Genet 2020; 50:102395. [PMID: 33070049 DOI: 10.1016/j.fsigen.2020.102395] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2020] [Revised: 08/25/2020] [Accepted: 09/15/2020] [Indexed: 12/31/2022]
Abstract
Predicting appearance phenotypes from genotypes is relevant for various areas of human genetic research and applications such as genetic epidemiology, human history, anthropology, and particularly in forensics. Many appearance phenotypes, and thus their underlying genotypes, are highly correlated, with pigmentation traits serving as primary examples. However, all available genetic prediction models, including those for pigmentation traits currently used in forensic DNA phenotyping, ignore phenotype correlations. Here, we investigated the impact of appearance phenotype correlations on genetic appearance prediction in the exemplary case of three pigmentation traits. We used data for categorical eye, hair and skin colour as well as 41 DNA markers utilized in the recently established HIrisPlex-S system from 762 individuals with complete phenotype and genotype information. Based on these data, we performed genetic prediction modelling of eye, hair and skin colour via three different strategies, namely the established approach of predicting phenotypes solely based on genotypes while not considering phenotype correlations, and two novel approaches that considered phenotype correlations, either incorporating truly observed correlated phenotypes or DNA-predicted correlated phenotypes in addition to the DNA predictors. We found that using truly observed correlated pigmentation phenotypes as additional predictors increased the DNA-based prediction accuracies for almost all eye, hair and skin colour categories, with the largest increase for intermediate eye colour, brown hair colour, dark to black skin colour, and particularly for dark skin colour. Outcomes of dedicated computer simulations suggest that this prediction accuracy increase is due to the additional genetic information that is implicitly provided by the truly observed correlated pigmentation phenotypes used, yet not covered by the DNA predictors applied. In contrast, considering DNA-predicted correlated pigmentation phenotypes as additional predictors did not improve the performance of the genetic prediction of eye, hair and skin colour, which was in line with the results from our computer simulations. Hence, in practical applications of DNA-based appearance prediction where no phenotype knowledge is available, such as in forensic DNA phenotyping, it is not advised to use DNA-predicted correlated phenotypes as predictors in addition to the DNA predictors. In the very least, this is not recommended for the pigmentation traits and the established pigmentation DNA predictors tested here.
Collapse
Affiliation(s)
- Yan Chen
- Department of Genetic Identification, Erasmus MC University Medical Center Rotterdam, Rotterdam, the Netherlands; Key Laboratory of Genomic and Precision Medicine, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, China; University of Chinese Academy of Sciences, Beijing, China
| | - Wojciech Branicki
- Malopolska Centre of Biotechnology, Jagiellonian University, Krakow, Poland
| | - Susan Walsh
- Department of Biology, Indiana University Purdue University Indianapolis (IUPUI), Indianapolis, IN, USA
| | - Michael Nothnagel
- Cologne Center for Genomics, University of Cologne, Cologne, Germany; University Hospital Cologne, Cologne, Germany
| | - Manfred Kayser
- Department of Genetic Identification, Erasmus MC University Medical Center Rotterdam, Rotterdam, the Netherlands
| | - Fan Liu
- Department of Genetic Identification, Erasmus MC University Medical Center Rotterdam, Rotterdam, the Netherlands; Key Laboratory of Genomic and Precision Medicine, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, China; University of Chinese Academy of Sciences, Beijing, China.
| | | |
Collapse
|
2
|
Zeisel SH. Precision (Personalized) Nutrition: Understanding Metabolic Heterogeneity. Annu Rev Food Sci Technol 2020; 11:71-92. [DOI: 10.1146/annurev-food-032519-051736] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
People differ in their requirements for and responses to nutrients and bioactive molecules in the diet. Many inputs contribute to metabolic heterogeneity (including variations in genetics, epigenetics, microbiome, lifestyle, diet intake, and environmental exposure). Precision nutrition is not about developing unique prescriptions for individual people but rather about stratifying people into different subgroups of the population on the basis of biomarkers of the above-listed sources of metabolic variation and then using this stratification to better estimate the different subgroups’ dietary requirements, thereby enabling better dietary recommendations and interventions. The hope is that we will be able to subcategorize people into ever-smaller groups that can be targeted in terms of recommendations, but we will never achieve this at the individual level, thus, the choice of precision nutrition rather than personalized nutrition to designate this new field. This review focuses mainly on genetically related sources of metabolic heterogeneity and identifies challenges that need to be overcome to achieve a full understanding of the complex interactions between the many sources of metabolic heterogeneity that make people differ from one another in their requirements for and responses to foods. It also discusses the commercial applications of precision nutrition.
Collapse
Affiliation(s)
- Steven H. Zeisel
- Nutrition Research Institute, Department of Nutrition, University of North Carolina, Kannapolis, North Carolina 28081, USA
| |
Collapse
|
3
|
Abstract
Classically, phenotype is what is observed, and genotype is the genetic makeup. Statistical studies aim to project phenotypic likelihoods of genotypic patterns. The traditional genotype-to-phenotype theory embraces the view that the encoded protein shape together with gene expression level largely determines the resulting phenotypic trait. Here, we point out that the molecular biology revolution at the turn of the century explained that the gene encodes not one but ensembles of conformations, which in turn spell all possible gene-associated phenotypes. The significance of a dynamic ensemble view is in understanding the linkage between genetic change and the gained observable physical or biochemical characteristics. Thus, despite the transformative shift in our understanding of the basis of protein structure and function, the literature still commonly relates to the classical genotype-phenotype paradigm. This is important because an ensemble view clarifies how even seemingly small genetic alterations can lead to pleiotropic traits in adaptive evolution and in disease, why cellular pathways can be modified in monogenic and polygenic traits, and how the environment may tweak protein function.
Collapse
Affiliation(s)
- Ruth Nussinov
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, Maryland, United States of America
- Sackler Institute of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Chung-Jung Tsai
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, Maryland, United States of America
| | - Hyunbum Jang
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, Maryland, United States of America
| |
Collapse
|
4
|
Jiang T, Huang M, Jiang T, Gu Y, Wang Y, Wu Y, Ma H, Jin G, Dai J, Hu Z. Genome-wide compound heterozygosity analysis highlighted 4 novel susceptibility loci for congenital heart disease in Chinese population. Clin Genet 2018; 94:296-302. [PMID: 29774522 DOI: 10.1111/cge.13384] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2018] [Revised: 04/25/2018] [Accepted: 05/15/2018] [Indexed: 01/25/2023]
Abstract
Genome-wide association studies (GWASs) have achieved great success in deciphering the genetic cause of congenital heart disease (CHD). However, the heritability of CHD remains to be clarified, and numerous genetic factors responsible for occurrence of CHD are yet unclear. In this study, we performed a genome-wide search for relaxed forms of compound heterozygosity (CH) in association with CHD using our existing GWAS data including 2265 individuals (957 CHD cases and 1308 controls). CollapsABEL was used to iteratively test the association between the CH genotype and the CHD phenotype in a sliding window manner. We highlighted 17 genetic loci showing suggestive CH-like associations with CHD (P < 5 × 10-8 ), among which 4 genetic loci had expression quantitative trait loci (eQTL) effects in blood (PeQTL < 0.01). After conditional association analysis, each loci had only 1 independently effective signal reaching the significance threshold (rs2071477/rs3129299 at 6p21.32, P = 2.47 × 10-10 ; rs10773097/rs2880921 at 12q24.31, P = 3.30 × 10-8 ; rs73032040/rs7259476 at 19q13.11, P = 1.14 × 10-8 ; rs10416386/rs4239517 at 19q13.31, P = 1.15 × 10-9 ), together explained 7.83% of the CHD variance. Among these 4 associated loci, outstanding candidates for CHD-associated genes included UBC, CFM2, ZNF302, LYPD3 and CADM4. Although replication studies with larger sample size are warranted, the first CH GWAS of CHD may extend our current knowledge of the genetic contributions to CHD in the Han Chinese population.
Collapse
Affiliation(s)
- T Jiang
- State Key Laboratory of Reproductive Medicine, Nanjing Medical University, Nanjing, China.,Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - M Huang
- Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - T Jiang
- State Key Laboratory of Reproductive Medicine, Nanjing Medical University, Nanjing, China.,Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - Y Gu
- State Key Laboratory of Reproductive Medicine, Nanjing Medical University, Nanjing, China.,Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - Y Wang
- State Key Laboratory of Reproductive Medicine, Nanjing Medical University, Nanjing, China.,Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - Y Wu
- Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - H Ma
- Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - G Jin
- State Key Laboratory of Reproductive Medicine, Nanjing Medical University, Nanjing, China.,Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - J Dai
- Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| | - Z Hu
- State Key Laboratory of Reproductive Medicine, Nanjing Medical University, Nanjing, China.,Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
| |
Collapse
|
5
|
Abstract
Relatedness within a sample can be of ancient (population stratification) or recent (familial structure) origin, and can either be known (pedigree data) or unknown (cryptic relatedness). All of these forms of familial relatedness have the potential to confound the results of genome-wide association studies. This chapter reviews the major methods available to researchers to adjust for the biases introduced by relatedness and maximize power to detect associations. The advantages and disadvantages of different methods are presented with reference to elements of study design, population characteristics, and computational requirements.
Collapse
Affiliation(s)
- Russell Thomson
- Centre for Research in Mathematics, School of Computing, Engineering and Mathematics, Western Sydney University, Parramatta, Australia.
| | - Rebekah McWhirter
- Menzies Institute for Medical Research, University of Tasmania, Hobart, TAS, Australia
| |
Collapse
|
6
|
Zhong K, Zhu G, Jing X, Hendriks AEJ, Drop SLS, Ikram MA, Gordon S, Zeng C, Uitterlinden AG, Martin NG, Liu F, Kayser M. Genome-wide compound heterozygote analysis highlights alleles associated with adult height in Europeans. Hum Genet 2017; 136:1407-1417. [PMID: 28921393 PMCID: PMC5702380 DOI: 10.1007/s00439-017-1842-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2017] [Accepted: 08/26/2017] [Indexed: 01/08/2023]
Abstract
Adult height is the most widely genetically studied common trait in humans; however, the trait variance explainable by currently known height-associated single nucleotide polymorphisms (SNPs) identified from the previous genome-wide association studies (GWAS) is yet far from complete given the high heritability of this complex trait. To exam if compound heterozygotes (CH) may explain extra height variance, we conducted a genome-wide analysis to screen for CH in association with adult height in 10,631 Dutch Europeans enriched with extremely tall people, using our recently developed method implemented in the software package CollapsABEL. The analysis identified six regions (3q23, 5q35.1, 6p21.31, 6p21.33, 7q21.2, and 9p24.3), where multiple pairs of SNPs as CH showed genome-wide significant association with height (P < 1.67 × 10−10). Of those, 9p24.3 represents a novel region influencing adult height, whereas the others have been highlighted in the previous GWAS on height based on analysis of individual SNPs. A replication analysis in 4080 Australians of European ancestry confirmed the significant CH-like association at 9p24.3 (P < 0.05). Together, the collapsed genotypes at these six loci explained 2.51% of the height variance (after adjusting for sex and age), compared with 3.23% explained by the 14 top-associated SNPs at 14 loci identified by traditional GWAS in the same data set (P < 5 × 10−8). Overall, our study empirically demonstrates that CH plays an important role in adult height and may explain a proportion of its “missing heritability”. Moreover, our findings raise promising expectations for other highly polygenic complex traits to explain missing heritability identifiable through CH-like associations.
Collapse
Affiliation(s)
- Kaiyin Zhong
- Department of Genetic Identification, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands
| | - Gu Zhu
- Queensland Institute of Medical Research, Brisbane, 4029, Australia
| | - Xiaoxi Jing
- Key Laboratory of Genomic and Precision Medicine, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, China.,University of Chinese Academy of Sciences, Beijing, China
| | - A Emile J Hendriks
- Division of Endocrinology, Department of Pediatrics, Sophia Children's Hospital, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands.,Department of Pediatrics, University of Cambridge, Cambridge, UK
| | - Sten L S Drop
- Division of Endocrinology, Department of Pediatrics, Sophia Children's Hospital, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands
| | - M Arfan Ikram
- Department of Internal Medicine, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands
| | - Scott Gordon
- Queensland Institute of Medical Research, Brisbane, 4029, Australia
| | - Changqing Zeng
- Key Laboratory of Genomic and Precision Medicine, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, China
| | - Andre G Uitterlinden
- Department of Internal Medicine, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands.,Department of Epidemiology, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands
| | | | - Fan Liu
- Department of Genetic Identification, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands. .,Key Laboratory of Genomic and Precision Medicine, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, China. .,University of Chinese Academy of Sciences, Beijing, China.
| | - Manfred Kayser
- Department of Genetic Identification, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands.
| |
Collapse
|
7
|
Zhong K, Verkouteren JA, Jacobs LC, Uitterlinden AG, Hofman A, Liu F, Nijsten T, Kayser M. Pigmentation-Independent Susceptibility Loci for Actinic Keratosis Highlighted by Compound Heterozygosity Analysis. J Invest Dermatol 2017; 137:77-84. [DOI: 10.1016/j.jid.2016.09.007] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2016] [Revised: 08/29/2016] [Accepted: 09/02/2016] [Indexed: 10/21/2022]
|
8
|
Liu F, Hamer MA, Deelen J, Lall JS, Jacobs L, van Heemst D, Murray PG, Wollstein A, de Craen AJM, Uh HW, Zeng C, Hofman A, Uitterlinden AG, Houwing-Duistermaat JJ, Pardo LM, Beekman M, Slagboom PE, Nijsten T, Kayser M, Gunn DA. The MC1R Gene and Youthful Looks. Curr Biol 2016; 26:1213-20. [PMID: 27133870 DOI: 10.1016/j.cub.2016.03.008] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2015] [Revised: 02/12/2016] [Accepted: 03/01/2016] [Indexed: 11/28/2022]
Abstract
Looking young for one's age has been a desire since time immemorial. This desire is attributable to the belief that appearance reflects health and fecundity. Indeed, perceived age predicts survival [1] and associates with molecular markers of aging such as telomere length [2]. Understanding the underlying molecular biology of perceived age is vital for identifying new aging therapies among other purposes, but studies are lacking thus far. As a first attempt, we performed genome-wide association studies (GWASs) of perceived facial age and wrinkling estimated from digital facial images by analyzing over eight million SNPs in 2,693 elderly Dutch Europeans from the Rotterdam Study. The strongest genetic associations with perceived facial age were found for multiple SNPs in the MC1R gene (p < 1 × 10(-7)). This effect was enhanced for a compound heterozygosity marker constructed from four pre-selected functional MC1R SNPs (p = 2.69 × 10(-12)), which was replicated in 599 Dutch Europeans from the Leiden Longevity Study (p = 0.042) and in 1,173 Europeans of the TwinsUK Study (p = 3 × 10(-3)). Individuals carrying the homozygote MC1R risk haplotype looked on average up to 2 years older than non-carriers. This association was independent of age, sex, skin color, and sun damage (wrinkling, pigmented spots) and persisted through different sun-exposure levels. Hence, a role for MC1R in youthful looks independent of its known melanin synthesis function is suggested. Our study uncovers the first genetic evidence explaining why some people look older for their age and provides new leads for further investigating the biological basis of how old or young people look.
Collapse
Affiliation(s)
- Fan Liu
- Key Laboratory of Genomic and Precision Medicine, China Gastrointestinal Cancer Research Center, Beijing Institute of Genomics, Chinese Academy of Sciences, No.1 Beichen West Road, Chaoyang District, Beijing 100101, China; Department of Genetic Identification, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands
| | - Merel A Hamer
- Department of Dermatology, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands
| | - Joris Deelen
- Department of Molecular Epidemiology, Leiden University Medical Center, P.O. Box 9600, 2300 RC Leiden, the Netherlands
| | - Japal S Lall
- Unilever R&D, Colworth Science Park, Sharnbrook, Bedfordshire MK44 1LQ, UK
| | - Leonie Jacobs
- Department of Dermatology, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands
| | - Diana van Heemst
- Department of Gerontology and Geriatrics, Leiden University Medical Center, P.O. Box 9600, 2300 RC Leiden, the Netherlands
| | - Peter G Murray
- Unilever R&D, Colworth Science Park, Sharnbrook, Bedfordshire MK44 1LQ, UK
| | - Andreas Wollstein
- Department of Genetic Identification, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands; Section of Evolutionary Biology, Department of Biology II, Ludwig Maximilians University Munich, Großhaderner Str. 2, 82152 Planegg-Martinsried, Germany
| | - Anton J M de Craen
- Department of Gerontology and Geriatrics, Leiden University Medical Center, P.O. Box 9600, 2300 RC Leiden, the Netherlands
| | - Hae-Won Uh
- Department of Medical Statistics and Bioinformatics, Leiden University Medical Center, P.O. Box 9600, 2300 RC Leiden, the Netherlands
| | - Changqing Zeng
- Key Laboratory of Genomic and Precision Medicine, China Gastrointestinal Cancer Research Center, Beijing Institute of Genomics, Chinese Academy of Sciences, No.1 Beichen West Road, Chaoyang District, Beijing 100101, China
| | - Albert Hofman
- Department of Epidemiology, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands
| | - André G Uitterlinden
- Department of Epidemiology, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands; Department of Internal Medicine, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands
| | - Jeanine J Houwing-Duistermaat
- Department of Medical Statistics and Bioinformatics, Leiden University Medical Center, P.O. Box 9600, 2300 RC Leiden, the Netherlands; Department of Statistics, School of Mathematics, University of Leeds, Leeds LS2 9JT, UK
| | - Luba M Pardo
- Department of Dermatology, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands
| | - Marian Beekman
- Department of Molecular Epidemiology, Leiden University Medical Center, P.O. Box 9600, 2300 RC Leiden, the Netherlands
| | - P Eline Slagboom
- Department of Molecular Epidemiology, Leiden University Medical Center, P.O. Box 9600, 2300 RC Leiden, the Netherlands
| | - Tamar Nijsten
- Department of Dermatology, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands
| | - Manfred Kayser
- Department of Genetic Identification, Erasmus MC University Medical Center Rotterdam, P.O. Box 2040, 3000 CA Rotterdam, the Netherlands.
| | - David A Gunn
- Unilever R&D, Colworth Science Park, Sharnbrook, Bedfordshire MK44 1LQ, UK.
| |
Collapse
|
9
|
Zhong K, Karssen LC, Kayser M, Liu F. CollapsABEL: an R library for detecting compound heterozygote alleles in genome-wide association studies. BMC Bioinformatics 2016; 17:156. [PMID: 27059780 PMCID: PMC4826552 DOI: 10.1186/s12859-016-1006-9] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2015] [Accepted: 03/31/2016] [Indexed: 11/10/2022] Open
Abstract
Background Compound Heterozygosity (CH) in classical genetics is the presence of two different recessive mutations at a particular gene locus. A relaxed form of CH alleles may account for an essential proportion of the missing heritability, i.e. heritability of phenotypes so far not accounted for by single genetic variants. Methods to detect CH-like effects in genome-wide association studies (GWAS) may facilitate explaining the missing heritability, but to our knowledge no viable software tools for this purpose are currently available. Results In this work we present the Generalized Compound Double Heterozygosity (GCDH) test and its implementation in the R package CollapsABEL. Time-consuming procedures are optimized for computational efficiency using Java or C++. Intermediate results are stored either in an SQL database or in a so-called big.matrix file to achieve reasonable memory footprint. Our large scale simulation studies show that GCDH is capable of discovering genetic associations due to CH-like interactions with much higher power than a conventional single-SNP approach under various settings, whether the causal genetic variations are available or not. CollapsABEL provides a user-friendly pipeline for genotype collapsing, statistical testing, power estimation, type I error control and graphics generation in the R language. Conclusions CollapsABEL provides a computationally efficient solution for screening general forms of CH alleles in densely imputed microarray or whole genome sequencing datasets. The GCDH test provides an improved power over single-SNP based methods in detecting the prevalence of CH in human complex phenotypes, offering an opportunity for tackling the missing heritability problem. Binary and source packages of CollapsABEL are available on CRAN (https://cran.r-project.org/web/packages/CollapsABEL) and the website of the GenABEL project (http://www.genabel.org/packages). Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-1006-9) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Kaiyin Zhong
- Department of Genetic Identification, Erasmus University Medical Center Rotterdam, Rotterdam, The Netherlands
| | | | - Manfred Kayser
- Department of Genetic Identification, Erasmus University Medical Center Rotterdam, Rotterdam, The Netherlands
| | - Fan Liu
- Department of Genetic Identification, Erasmus University Medical Center Rotterdam, Rotterdam, The Netherlands. .,Key Laboratory of Genomic and Precision Medicine, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, China.
| |
Collapse
|
10
|
Liu F, Visser M, Duffy DL, Hysi PG, Jacobs LC, Lao O, Zhong K, Walsh S, Chaitanya L, Wollstein A, Zhu G, Montgomery GW, Henders AK, Mangino M, Glass D, Bataille V, Sturm RA, Rivadeneira F, Hofman A, van IJcken WFJ, Uitterlinden AG, Palstra RJTS, Spector TD, Martin NG, Nijsten TEC, Kayser M. Genetics of skin color variation in Europeans: genome-wide association studies with functional follow-up. Hum Genet 2015; 134:823-35. [PMID: 25963972 PMCID: PMC4495261 DOI: 10.1007/s00439-015-1559-0] [Citation(s) in RCA: 99] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2015] [Accepted: 04/20/2015] [Indexed: 02/05/2023]
Abstract
In the International Visible Trait Genetics (VisiGen) Consortium, we investigated the genetics of human skin color by combining a series of genome-wide association studies (GWAS) in a total of 17,262 Europeans with functional follow-up of discovered loci. Our GWAS provide the first genome-wide significant evidence for chromosome 20q11.22 harboring the ASIP gene being explicitly associated with skin color in Europeans. In addition, genomic loci at 5p13.2 (SLC45A2), 6p25.3 (IRF4), 15q13.1 (HERC2/OCA2), and 16q24.3 (MC1R) were confirmed to be involved in skin coloration in Europeans. In follow-up gene expression and regulation studies of 22 genes in 20q11.22, we highlighted two novel genes EIF2S2 and GSS, serving as competing functional candidates in this region and providing future research lines. A genetically inferred skin color score obtained from the 9 top-associated SNPs from 9 genes in 940 worldwide samples (HGDP-CEPH) showed a clear gradual pattern in Western Eurasians similar to the distribution of physical skin color, suggesting the used 9 SNPs as suitable markers for DNA prediction of skin color in Europeans and neighboring populations, relevant in future forensic and anthropological
investigations.
Collapse
Affiliation(s)
- Fan Liu
- Department of Forensic Molecular Biology, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands,
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
11
|
Jacobs LC, Liu F, Pardo LM, Hofman A, Uitterlinden AG, Kayser M, Nijsten T. IRF4, MC1R and TYR genes are risk factors for actinic keratosis independent of skin color. Hum Mol Genet 2015; 24:3296-303. [PMID: 25724930 DOI: 10.1093/hmg/ddv076] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2014] [Accepted: 02/24/2015] [Indexed: 12/13/2022] Open
Abstract
Actinic keratosis (AK) is a pre-malignant skin disease, highly prevalent in elderly Europeans. This study investigates genetic susceptibility to AK with a genome-wide association study (GWAS). A full body skin examination was performed in 3194 elderly individuals from the Rotterdam Study (RS) of exclusive north-western European origin (aged 51-99 years, 45% male). Physicians graded the number of AK into four severity levels: none (76%), 1-3 (14%), 4-9 (6%) and ≥10 (5%), and skin color was quantified using a spectrophotometer on sun-unexposed skin. A GWAS for AK severity was conducted, where promising signals at IRF4 and MC1R (P < 4.2 × 10(-7)) were successfully replicated in an additional cohort of 623 RS individuals (IRF4, rs12203592, Pcombined = 6.5 × 10(-13) and MC1R, rs139810560, Pcombined = 4.1 × 10(-9)). Further, in an analysis of ten additional well-known human pigmentation genes, TYR also showed significant association with AK (rs1393350, P = 5.3 × 10(-4)) after correction for multiple testing. Interestingly, the strength and significance of above-mentioned associations retained largely the same level after skin color adjustment. Overall, our data strongly suggest that IRF4, MC1R and TYR genes likely have pleiotropic effects, a combination of pigmentation and oncogenic functions, resulting in an increased risk of AK.
Collapse
Affiliation(s)
| | - Fan Liu
- Department of Forensic Molecular Biology
| | | | | | - André G Uitterlinden
- Department of Epidemiology and Department of Internal Medicine, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands
| | | | | |
Collapse
|
12
|
A Genome-Wide Association Study Identifies the Skin Color Genes IRF4, MC1R, ASIP, and BNC2 Influencing Facial Pigmented Spots. J Invest Dermatol 2015; 135:1735-1742. [PMID: 25705849 DOI: 10.1038/jid.2015.62] [Citation(s) in RCA: 84] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2014] [Revised: 01/27/2015] [Accepted: 02/09/2015] [Indexed: 02/02/2023]
Abstract
Facial pigmented spots are a common skin aging feature, but genetic predisposition has yet to be thoroughly investigated. We conducted a genome-wide association study for pigmented spots in 2,844 Dutch Europeans from the Rotterdam Study (mean age: 66.9±8.0 years; 47% male). Using semi-automated image analysis of high-resolution digital facial photographs, facial pigmented spots were quantified as the percentage of affected skin area (mean women: 2.0% ±0.9, men: 0.9% ±0.6). We identified genome-wide significant association with pigmented spots at three genetic loci: IRF4 (rs12203592, P=1.8 × 10(-27)), MC1R (compound heterozygosity score, P=2.3 × 10(-24)), and RALY/ASIP (rs6059655, P=1.9 × 10(-9)). In addition, after adjustment for the other three top-associated loci the BNC2 locus demonstrated significant association (rs62543565, P=2.3 × 10(-8)). The association signals observed at all four loci were successfully replicated (P<0.05) in an independent Dutch cohort (Leiden Longevity Study n=599). Although the four genes have previously been associated with skin color variation and skin cancer risk, all association signals remained highly significant (P<2 × 10(-8)) when conditioning the association analyses on skin color. We conclude that genetic variations in IRF4, MC1R, RALY/ASIP, and BNC2 contribute to the acquired amount of facial pigmented spots during aging, through pathways independent of the basal melanin production.
Collapse
|
13
|
Walsh S, Chaitanya L, Clarisse L, Wirken L, Draus-Barini J, Kovatsi L, Maeda H, Ishikawa T, Sijen T, de Knijff P, Branicki W, Liu F, Kayser M. Developmental validation of the HIrisPlex system: DNA-based eye and hair colour prediction for forensic and anthropological usage. Forensic Sci Int Genet 2013; 9:150-61. [PMID: 24528593 DOI: 10.1016/j.fsigen.2013.12.006] [Citation(s) in RCA: 119] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2013] [Revised: 12/13/2013] [Accepted: 12/17/2013] [Indexed: 12/19/2022]
Abstract
Forensic DNA Phenotyping or 'DNA intelligence' tools are expected to aid police investigations and find unknown individuals by providing information on externally visible characteristics of unknown suspects, perpetrators and missing persons from biological samples. This is especially useful in cases where conventional DNA profiling or other means remain non-informative. Recently, we introduced the HIrisPlex system, capable of predicting both eye and hair colour from DNA. In the present developmental validation study, we demonstrate that the HIrisPlex assay performs in full agreement with the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines providing an essential prerequisite for future HIrisPlex applications to forensic casework. The HIrisPlex assay produces complete profiles down to only 63 pg of DNA. Species testing revealed human specificity for a complete HIrisPlex profile, while only non-human primates showed the closest full profile at 20 out of the 24 DNA markers, in all animals tested. Rigorous testing of simulated forensic casework samples such as blood, semen, saliva stains, hairs with roots as well as extremely low quantity touch (trace) DNA samples, produced complete profiles in 88% of cases. Concordance testing performed between five independent forensic laboratories displayed consistent reproducible results on varying types of DNA samples. Due to its design, the assay caters for degraded samples, underlined here by results from artificially degraded DNA and from simulated casework samples of degraded DNA. This aspect was also demonstrated previously on DNA samples from human remains up to several hundreds of years old. With this paper, we also introduce enhanced eye and hair colour prediction models based on enlarged underlying databases of HIrisPlex genotypes and eye/hair colour phenotypes (eye colour: N = 9188 and hair colour: N = 1601). Furthermore, we present an online web-based system for individual eye and hair colour prediction from full and partial HIrisPlex DNA profiles. By demonstrating that the HIrisPlex assay is fully compatible with the SWGDAM guidelines, we provide the first forensically validated DNA test system for parallel eye and hair colour prediction now available to forensic laboratories for immediate casework application, including missing person cases. Given the robustness and sensitivity described here and in previous work, the HIrisPlex system is also suitable for analysing old and ancient DNA in anthropological and evolutionary studies.
Collapse
Affiliation(s)
- Susan Walsh
- Department of Forensic Molecular Biology, Erasmus MC University Medical Centre Rotterdam, Rotterdam, The Netherlands
| | - Lakshmi Chaitanya
- Department of Forensic Molecular Biology, Erasmus MC University Medical Centre Rotterdam, Rotterdam, The Netherlands
| | - Lindy Clarisse
- Department of Human Biological Traces, Netherlands Forensic Institute, The Hague, The Netherlands
| | - Laura Wirken
- Forensic Laboratory for DNA Research, Department of Human Genetics, Leiden University Medical Center, Leiden, The Netherlands
| | | | - Leda Kovatsi
- Laboratory of Forensic Medicine & Toxicology, School of Medicine, Aristotle University of Thessaloniki, Thessaloniki, Greece
| | - Hitoshi Maeda
- Department of Legal Medicine, Osaka City University, Medical School, Osaka, Japan
| | - Takaki Ishikawa
- Department of Legal Medicine, Osaka City University, Medical School, Osaka, Japan; Division of Legal Medicine, Faculty of Medicine, Tottori University, 86 Nichicho Yonago, Japan
| | - Titia Sijen
- Department of Human Biological Traces, Netherlands Forensic Institute, The Hague, The Netherlands
| | - Peter de Knijff
- Forensic Laboratory for DNA Research, Department of Human Genetics, Leiden University Medical Center, Leiden, The Netherlands
| | - Wojciech Branicki
- Section of Forensic Genetics, Institute of Forensic Research, Kraków, Poland; Department of Genetics and Evolution, Institute of Zoology, Jagiellonian University, Kraków, Poland
| | - Fan Liu
- Department of Forensic Molecular Biology, Erasmus MC University Medical Centre Rotterdam, Rotterdam, The Netherlands
| | - Manfred Kayser
- Department of Forensic Molecular Biology, Erasmus MC University Medical Centre Rotterdam, Rotterdam, The Netherlands.
| |
Collapse
|
14
|
Dissecting genome-wide association signals for loss-of-function phenotypes in sorghum flavonoid pigmentation traits. G3-GENES GENOMES GENETICS 2013; 3:2085-94. [PMID: 24048646 PMCID: PMC3815067 DOI: 10.1534/g3.113.008417] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Genome-wide association studies are a powerful method to dissect the genetic basis of traits, although in practice the effects of complex genetic architecture and population structure remain poorly understood. To compare mapping strategies we dissected the genetic control of flavonoid pigmentation traits in the cereal grass sorghum by using high-resolution genotyping-by-sequencing single-nucleotide polymorphism markers. Studying the grain tannin trait, we find that general linear models (GLMs) are not able to precisely map tan1-a, a known loss-of-function allele of the Tannin1 gene, with either a small panel (n = 142) or large association panel (n = 336), and that indirect associations limit the mapping of the Tannin1 locus to Mb-resolution. A GLM that accounts for population structure (Q) or standard mixed linear model that accounts for kinship (K) can identify tan1-a, whereas a compressed mixed linear model performs worse than the naive GLM. Interestingly, a simple loss-of-function genome scan, for genotype-phenotype covariation only in the putative loss-of-function allele, is able to precisely identify the Tannin1 gene without considering relatedness. We also find that the tan1-a allele can be mapped with gene resolution in a biparental recombinant inbred line family (n = 263) using genotyping-by-sequencing markers but lower precision in the mapping of vegetative pigmentation traits suggest that consistent gene-level resolution will likely require larger families or multiple recombinant inbred lines. These findings highlight that complex association signals can emerge from even the simplest traits given epistasis and structured alleles, but that gene-resolution mapping of these traits is possible with high marker density and appropriate models.
Collapse
|
15
|
Bacanu SA. Testing for modes of inheritance involving compound heterozygotes. Genet Epidemiol 2013; 37:522-8. [PMID: 23633151 DOI: 10.1002/gepi.21732] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2012] [Revised: 03/26/2013] [Accepted: 04/01/2013] [Indexed: 11/09/2022]
Abstract
Functional variants change the protein product or the expression of genes. Due to the latest advances in sequencing technology, most known functional variants can now be assayed in a cost-effective manner. However, to fully use the information from functional variants, researchers need to model the joint effect of these variants. In this article, we propose methods that model the action/interaction of loss-of-function (LOF) mutations, i.e., those mutations that eliminate the protein product of a gene. When multiple LOFs occur in the same causal gene/region, their effect on a phenotype might depend on whether these mutations lie on the same DNA strand/haplotype. When compared to LOFs occurring on the same strand, if these mutations lie on different strands, both copies of the gene are impaired and the impact on the relevant phenotypes is likely to be more severe. To use the information from LOF strand colocalization, we propose three methods that utilize the information from the estimated number of affected strands. We compare the performance of the proposed and competing methods by using simulations of common and rare LOF variants. Two of the proposed methods exhibited desirable power profiles, the first for both common and rare LOFs and the second only for common LOFs. One of the existing methods, collapsed double heterozygosity, exhibits good power to detect compound models for rare variants, especially when no haplotype harbors two or more rare alleles. Consequently, we recommend these three methods to be used for the analysis of functional variants coming from sequencing studies.
Collapse
Affiliation(s)
- Silviu-Alin Bacanu
- Virginia Institute for Psychiatric and Behavioral Genetics BIOTECH I, Richmond, Virginia, USA.
| |
Collapse
|
16
|
Liu F, Wen B, Kayser M. Colorful DNA polymorphisms in humans. Semin Cell Dev Biol 2013; 24:562-75. [PMID: 23587773 DOI: 10.1016/j.semcdb.2013.03.013] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2013] [Accepted: 03/26/2013] [Indexed: 10/26/2022]
Abstract
In this review article we summarize current knowledge on how variation on the DNA level influences human pigmentation including color variation of iris, hair, and skin. We review recent progress in the field of human pigmentation genetics by focusing on the genes and DNA polymorphisms discovered to be involved in determining human pigmentation traits, their association with diseases particularly skin cancers, and their power to predict human eye, hair, and skin colors with potential utilization in forensic investigations.
Collapse
Affiliation(s)
- Fan Liu
- Department of Forensic Molecular Biology, Erasmus MC University Medical Center Rotterdam, Rotterdam, The Netherlands.
| | | | | |
Collapse
|
17
|
Comprehensive candidate gene study highlights UGT1A and BNC2 as new genes determining continuous skin color variation in Europeans. Hum Genet 2012; 132:147-58. [PMID: 23052946 DOI: 10.1007/s00439-012-1232-9] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2012] [Accepted: 09/18/2012] [Indexed: 01/09/2023]
Abstract
Natural variation in human skin pigmentation is primarily due to genetic causes rooted in recent evolutionary history. Genetic variants associated with human skin pigmentation confer risk of skin cancer and may provide useful information in forensic investigations. Almost all previous gene-mapping studies of human skin pigmentation were based on categorical skin color information known to oversimplify the continuous nature of human skin coloration. We digitally quantified skin color into hue and saturation dimensions for 5,860 Dutch Europeans based on high-resolution skin photographs. We then tested an extensive list of 14,185 single nucleotide polymorphisms in 281 candidate genes potentially involved in human skin pigmentation for association with quantitative skin color phenotypes. Confirmatory association was revealed for several known skin color genes including HERC2, MC1R, IRF4, TYR, OCA2, and ASIP. We identified two new skin color genes: genetic variants in UGT1A were significantly associated with hue and variants in BNC2 were significantly associated with saturation. Overall, digital quantification of human skin color allowed detecting new skin color genes. The variants identified in this study may also contribute to the risk of skin cancer. Our findings are also important for predicting skin color in forensic investigations.
Collapse
|
18
|
Walsh S, Liu F, Wollstein A, Kovatsi L, Ralf A, Kosiniak-Kamysz A, Branicki W, Kayser M. The HIrisPlex system for simultaneous prediction of hair and eye colour from DNA. Forensic Sci Int Genet 2012; 7:98-115. [PMID: 22917817 DOI: 10.1016/j.fsigen.2012.07.005] [Citation(s) in RCA: 249] [Impact Index Per Article: 20.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2012] [Revised: 06/25/2012] [Accepted: 07/23/2012] [Indexed: 02/03/2023]
Abstract
Recently, the field of predicting phenotypes of externally visible characteristics (EVCs) from DNA genotypes with the final aim of concentrating police investigations to find persons completely unknown to investigating authorities, also referred to as Forensic DNA Phenotyping (FDP), has started to become established in forensic biology. We previously developed and forensically validated the IrisPlex system for accurate prediction of blue and brown eye colour from DNA, and recently showed that all major hair colour categories are predictable from carefully selected DNA markers. Here, we introduce the newly developed HIrisPlex system, which is capable of simultaneously predicting both hair and eye colour from DNA. HIrisPlex consists of a single multiplex assay targeting 24 eye and hair colour predictive DNA variants including all 6 IrisPlex SNPs, as well as two prediction models, a newly developed model for hair colour categories and shade, and the previously developed IrisPlex model for eye colour. The HIrisPlex assay was designed to cope with low amounts of template DNA, as well as degraded DNA, and preliminary sensitivity testing revealed full DNA profiles down to 63pg input DNA. The power of the HIrisPlex system to predict hair colour was assessed in 1551 individuals from three different parts of Europe showing different hair colour frequencies. Using a 20% subset of individuals, while 80% were used for model building, the individual-based prediction accuracies employing a prediction-guided approach were 69.5% for blond, 78.5% for brown, 80% for red and 87.5% for black hair colour on average. Results from HIrisPlex analysis on worldwide DNA samples imply that HIrisPlex hair colour prediction is reliable independent of bio-geographic ancestry (similar to previous IrisPlex findings for eye colour). We furthermore demonstrate that it is possible to infer with a prediction accuracy of >86% if a brown-eyed, black-haired individual is of non-European (excluding regions nearby Europe) versus European (including nearby regions) bio-geographic origin solely from the strength of HIrisPlex eye and hair colour probabilities, which can provide extra intelligence for future forensic applications. The HIrisPlex system introduced here, including a single multiplex test assay, an interactive tool and prediction guide, and recommendations for reporting final outcomes, represents the first tool for simultaneously establishing categorical eye and hair colour of a person from DNA. The practical forensic application of the HIrisPlex system is expected to benefit cases where other avenues of investigation, including STR profiling, provide no leads on who the unknown crime scene sample donor or the unknown missing person might be.
Collapse
Affiliation(s)
- Susan Walsh
- Department of Forensic Molecular Biology, Erasmus MC University Medical Centre Rotterdam, The Netherlands
| | | | | | | | | | | | | | | |
Collapse
|