1
|
Fisher SA, Patrick K, Hoang T, Marcq E, Behrouzfar K, Young S, Miller TJ, Robinson BWS, Bueno R, Nowak AK, Lesterhuis WJ, Morahan G, Lake RA. The MexTAg collaborative cross: host genetics affects asbestos related disease latency, but has little influence once tumours develop. FRONTIERS IN TOXICOLOGY 2024; 6:1373003. [PMID: 38694815 PMCID: PMC11061428 DOI: 10.3389/ftox.2024.1373003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Accepted: 04/02/2024] [Indexed: 05/04/2024] Open
Abstract
Objectives: This study combines two innovative mouse models in a major gene discovery project to assess the influence of host genetics on asbestos related disease (ARD). Conventional genetics studies provided evidence that some susceptibility to mesothelioma is genetic. However, the identification of host modifier genes, the roles they may play, and whether they contribute to disease susceptibility remain unknown. Here we report a study designed to rapidly identify genes associated with mesothelioma susceptibility by combining the Collaborative Cross (CC) resource with the well-characterised MexTAg mesothelioma mouse model. Methods: The CC is a powerful mouse resource that harnesses over 90% of common genetic variation in the mouse species, allowing rapid identification of genes mediating complex traits. MexTAg mice rapidly, uniformly, and predictably develop mesothelioma, but only after asbestos exposure. To assess the influence of host genetics on ARD, we crossed 72 genetically distinct CC mouse strains with MexTAg mice and exposed the resulting CC-MexTAg (CCMT) progeny to asbestos and monitored them for traits including overall survival, the time to ARD onset (latency), the time between ARD onset and euthanasia (disease progression) and ascites volume. We identified phenotype-specific modifier genes associated with these traits and we validated the role of human orthologues in asbestos-induced carcinogenesis using human mesothelioma datasets. Results: We generated 72 genetically distinct CCMT strains and exposed their progeny (2,562 in total) to asbestos. Reflecting the genetic diversity of the CC, there was considerable variation in overall survival and disease latency. Surprisingly, however, there was no variation in disease progression, demonstrating that host genetic factors do have a significant influence during disease latency but have a limited role once disease is established. Quantitative trait loci (QTL) affecting ARD survival/latency were identified on chromosomes 6, 12 and X. Of the 97-protein coding candidate modifier genes that spanned these QTL, eight genes (CPED1, ORS1, NDUFA1, HS1BP3, IL13RA1, LSM8, TES and TSPAN12) were found to significantly affect outcome in both CCMT and human mesothelioma datasets. Conclusion: Host genetic factors affect susceptibility to development of asbestos associated disease. However, following mesothelioma establishment, genetic variation in molecular or immunological mechanisms did not affect disease progression. Identification of multiple candidate modifier genes and their human homologues with known associations in other advanced stage or metastatic cancers highlights the complexity of ARD and may provide a pathway to identify novel therapeutic targets.
Collapse
Affiliation(s)
- Scott A. Fisher
- National Centre for Asbestos Related Diseases (NCARD), Perth, WA, Australia
- School of Biomedical Sciences, The University of Western Australia, Perth, WA, Australia
- Institute for Respiratory Health, University of Western Australia, Perth, WA, Australia
| | - Kimberley Patrick
- National Centre for Asbestos Related Diseases (NCARD), Perth, WA, Australia
- School of Biomedical Sciences, The University of Western Australia, Perth, WA, Australia
- Institute for Respiratory Health, University of Western Australia, Perth, WA, Australia
| | - Tracy Hoang
- National Centre for Asbestos Related Diseases (NCARD), Perth, WA, Australia
- Institute for Respiratory Health, University of Western Australia, Perth, WA, Australia
| | - Elly Marcq
- Center for Oncological Research (CORE), University of Antwerp, Antwerp, Belgium
- Lab of Dendritic Cell Biology and Cancer Immunotherapy, VIB Center for Inflammation Research, Brussels, Belgium
- Brussels Center for Immunology, Vrije Universiteit Brussel, Brussels, Belgium
| | - Kiarash Behrouzfar
- National Centre for Asbestos Related Diseases (NCARD), Perth, WA, Australia
- School of Biomedical Sciences, The University of Western Australia, Perth, WA, Australia
- Institute for Respiratory Health, University of Western Australia, Perth, WA, Australia
| | - Sylvia Young
- Centre for Diabetes Research, Harry Perkins Institute of Medical Research, Perth, WA, Australia
| | - Timothy J. Miller
- Medical School, The University of Western Australia, Perth, WA, Australia
| | - Bruce W. S. Robinson
- National Centre for Asbestos Related Diseases (NCARD), Perth, WA, Australia
- Institute for Respiratory Health, University of Western Australia, Perth, WA, Australia
- Medical School, The University of Western Australia, Perth, WA, Australia
| | - Raphael Bueno
- Division of Thoracic Surgery, The Lung Center and the International Mesothelioma Program, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA, United States
| | - Anna K. Nowak
- National Centre for Asbestos Related Diseases (NCARD), Perth, WA, Australia
- Institute for Respiratory Health, University of Western Australia, Perth, WA, Australia
- Medical School, The University of Western Australia, Perth, WA, Australia
| | | | - Grant Morahan
- Centre for Diabetes Research, Harry Perkins Institute of Medical Research, Perth, WA, Australia
| | - Richard A. Lake
- National Centre for Asbestos Related Diseases (NCARD), Perth, WA, Australia
- Institute for Respiratory Health, University of Western Australia, Perth, WA, Australia
| |
Collapse
|
2
|
Garin V, Diallo C, Tékété ML, Théra K, Guitton B, Dagno K, Diallo AG, Kouressy M, Leiser W, Rattunde F, Sissoko I, Touré A, Nébié B, Samaké M, Kholovà J, Berger A, Frouin J, Pot D, Vaksmann M, Weltzien E, Témé N, Rami JF. Characterization of adaptation mechanisms in sorghum using a multireference back-cross nested association mapping design and envirotyping. Genetics 2024; 226:iyae003. [PMID: 38381593 PMCID: PMC10990433 DOI: 10.1093/genetics/iyae003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2023] [Accepted: 12/20/2023] [Indexed: 02/23/2024] Open
Abstract
Identifying the genetic factors impacting the adaptation of crops to environmental conditions is of key interest for conservation and selection purposes. It can be achieved using population genomics, and evolutionary or quantitative genetics. Here we present a sorghum multireference back-cross nested association mapping population composed of 3,901 lines produced by crossing 24 diverse parents to 3 elite parents from West and Central Africa-back-cross nested association mapping. The population was phenotyped in environments characterized by differences in photoperiod, rainfall pattern, temperature levels, and soil fertility. To integrate the multiparental and multi-environmental dimension of our data we proposed a new approach for quantitative trait loci (QTL) detection and parental effect estimation. We extended our model to estimate QTL effect sensitivity to environmental covariates, which facilitated the integration of envirotyping data. Our models allowed spatial projections of the QTL effects in agro-ecologies of interest. We utilized this strategy to analyze the genetic architecture of flowering time and plant height, which represents key adaptation mechanisms in environments like West Africa. Our results allowed a better characterization of well-known genomic regions influencing flowering time concerning their response to photoperiod with Ma6 and Ma1 being photoperiod-sensitive and the region of possible candidate gene Elf3 being photoperiod-insensitive. We also accessed a better understanding of plant height genetic determinism with the combined effects of phenology-dependent (Ma6) and independent (qHT7.1 and Dw3) genomic regions. Therefore, we argue that the West and Central Africa-back-cross nested association mapping and the presented analytical approach constitute unique resources to better understand adaptation in sorghum with direct application to develop climate-smart varieties.
Collapse
Affiliation(s)
- Vincent Garin
- Crop Physiology Laboratory, International Crops Research Institute for the Semi-Arid Tropics, Patancheru, 502 324, India
- CIRAD, UMR AGAP Institut, Montpellier, F-34398, France
- UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, F-34398, France
| | - Chiaka Diallo
- Sorghum Program, International Crops Research Institute for the Semi-Arid Tropics, Bamako, BP 320, Mali
- Département d’Enseignement et de Recherche des Sciences et Techniques Agricoles, Institut polytechnique rural de formation et de recherche appliquée de Katibougou, Koulikoro, BP 06, Mali
| | - Mohamed Lamine Tékété
- Institut d’Economie Rurale, Bamako, BP 262, Mali
- Faculté des Sciences et Techniques, Université des Sciences des Techniques et des Technologies de Bamako, Bamako, BP E 3206, Mali
| | | | - Baptiste Guitton
- CIRAD, UMR AGAP Institut, Montpellier, F-34398, France
- UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, F-34398, France
| | - Karim Dagno
- Institut d’Economie Rurale, Bamako, BP 262, Mali
| | | | | | - Willmar Leiser
- Sorghum Program, International Crops Research Institute for the Semi-Arid Tropics, Bamako, BP 320, Mali
| | - Fred Rattunde
- Agronomy Department, University of Wisconsin, Madison, WI 53705, WI, USA
| | - Ibrahima Sissoko
- Sorghum Program, International Crops Research Institute for the Semi-Arid Tropics, Bamako, BP 320, Mali
| | - Aboubacar Touré
- Sorghum Program, International Crops Research Institute for the Semi-Arid Tropics, Bamako, BP 320, Mali
| | - Baloua Nébié
- Dryland Crops Program, International Maize and Wheat Improvement Center (CIMMYT-Senegal) U/C CERAAS, Thiès, Po Box 3320, Senegal
| | - Moussa Samaké
- Faculté des Sciences et Techniques, Université des Sciences des Techniques et des Technologies de Bamako, Bamako, BP E 3206, Mali
| | - Jana Kholovà
- Crop Physiology Laboratory, International Crops Research Institute for the Semi-Arid Tropics, Patancheru, 502 324, India
- Department of Information Technologies, Faculty of Economics and Management, Czech University of Life Sciences, Prague, 165 00, Czech Republic
| | - Angélique Berger
- CIRAD, UMR AGAP Institut, Montpellier, F-34398, France
- UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, F-34398, France
| | - Julien Frouin
- CIRAD, UMR AGAP Institut, Montpellier, F-34398, France
- UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, F-34398, France
| | - David Pot
- CIRAD, UMR AGAP Institut, Montpellier, F-34398, France
- UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, F-34398, France
| | - Michel Vaksmann
- CIRAD, UMR AGAP Institut, Montpellier, F-34398, France
- UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, F-34398, France
| | - Eva Weltzien
- Sorghum Program, International Crops Research Institute for the Semi-Arid Tropics, Bamako, BP 320, Mali
- Agronomy Department, University of Wisconsin, Madison, WI 53705, WI, USA
| | - Niaba Témé
- Institut d’Economie Rurale, Bamako, BP 262, Mali
| | - Jean-François Rami
- CIRAD, UMR AGAP Institut, Montpellier, F-34398, France
- UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, F-34398, France
| |
Collapse
|
3
|
Ye H, Xu Z, Bello SF, Zhu Q, Kong S, Zheng M, Fang X, Jia X, Xu H, Zhang X, Nie Q. Haplotype analysis of genomic prediction by incorporating genomic pathway information based on high-density SNP marker in Chinese yellow-feathered chicken. Poult Sci 2023; 102:102549. [PMID: 36907129 PMCID: PMC10024239 DOI: 10.1016/j.psj.2023.102549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2022] [Revised: 01/16/2023] [Accepted: 01/27/2023] [Indexed: 02/09/2023] Open
Abstract
Genomic selection using single nucleotide polymorphism (SNP) markers is now intensively investigated in breeding and has been widely utilized for genetic improvement. Currently, several studies have used haplotype (consisting of multiallelic SNPs) for genomic prediction and revealed its performance advantage. In this study, we comprehensively evaluated the performance of haplotype models for genomic prediction in 15 traits, including 6 growth, 5 carcass, and 4 feeding traits in a Chinese yellow-feathered chicken population. We adopted 3 methods to define haplotypes from high-density SNP panels, and our strategy included combining Kyoto Encyclopedia of Genes and Genomes pathway information and considering linkage disequilibrium (LD) information. Our results showed an increase in prediction accuracy due to haplotypes ranging from -0.04∼27.16% in all traits, where the significant improvements were found in 12 traits. The estimates of haplotype epistasis heritability were strongly correlated with the accuracy increase by haplotype models. In addition, incorporating genomic annotation information could further increase the accuracy of the haplotype model, where the further increase in accuracy is significantly relative to the increase of relative haplotype epistasis heritability. The genomic prediction using LD information for constructing haplotypes has the best prediction performance among the 4 traits. These results uncovered that haplotype methods were beneficial for genomic prediction, and the accuracy could be further increased by incorporating genomic annotation information. Moreover, using LD information would potentially improve the performance of genomic prediction.
Collapse
Affiliation(s)
- Haoqiang Ye
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Zhenqiang Xu
- Wen's Nanfang Poultry Breeding Co. Ltd, Guangdong Province, Yunfu 527400, China
| | - Semiu Folaniyi Bello
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Qianghui Zhu
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China
| | - Shaofen Kong
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Ming Zheng
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Xiang Fang
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Xinzheng Jia
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, Foshan University, Foshan, 528225 China
| | - Haiping Xu
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Xiquan Zhang
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Qinghua Nie
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China.
| |
Collapse
|
4
|
Hong-Le T, Crouse WL, Keele GR, Holl K, Seshie O, Tschannen M, Craddock A, Das SK, Szalanczy AM, McDonald B, Grzybowski M, Klotz J, Sharma NK, Geurts AM, Key CCC, Hawkins G, Valdar W, Mott R, Solberg Woods LC. Genetic Mapping of Multiple Traits Identifies Novel Genes for Adiposity, Lipids, and Insulin Secretory Capacity in Outbred Rats. Diabetes 2023; 72:135-148. [PMID: 36219827 PMCID: PMC9797320 DOI: 10.2337/db22-0252] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 10/04/2022] [Indexed: 01/21/2023]
Abstract
Despite the successes of human genome-wide association studies, the causal genes underlying most metabolic traits remain unclear. We used outbred heterogeneous stock (HS) rats, coupled with expression data and mediation analysis, to identify quantitative trait loci (QTLs) and candidate gene mediators for adiposity, glucose tolerance, serum lipids, and other metabolic traits. Physiological traits were measured in 1,519 male HS rats, with liver and adipose transcriptomes measured in >410 rats. Genotypes were imputed from low-coverage whole-genome sequencing. Linear mixed models were used to detect physiological and expression QTLs (pQTLs and eQTLs, respectively), using both single nucleotide polymorphism (SNP)- and haplotype-based models for pQTL mapping. Genes with cis-eQTLs that overlapped pQTLs were assessed as causal candidates through mediation analysis. We identified 14 SNP-based pQTLs and 19 haplotype-based pQTLs, of which 10 were in common. Using mediation, we identified the following genes as candidate mediators of pQTLs: Grk5 for fat pad weight and serum triglyceride pQTLs on Chr1, Krtcap3 for fat pad weight and serum triglyceride pQTLs on Chr6, Ilrun for a fat pad weight pQTL on Chr20, and Rfx6 for a whole pancreatic insulin content pQTL on Chr20. Furthermore, we verified Grk5 and Ktrcap3 using gene knockdown/out models, thereby shedding light on novel regulators of obesity.
Collapse
Affiliation(s)
- Thu Hong-Le
- Genetics Institute, University College London, London, U.K
| | - Wesley L. Crouse
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC
| | | | - Katie Holl
- Medical College of Wisconsin, Milwaukee, WI
| | - Osborne Seshie
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | | | - Ann Craddock
- Department of Biochemistry, Wake Forest University School of Medicine, Winston-Salem, NC
| | - Swapan K. Das
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | - Alexandria M. Szalanczy
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | - Bailey McDonald
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | | | | | - Neeraj K. Sharma
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | | | - Chia-Chi Chuang Key
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | - Gregory Hawkins
- Department of Biochemistry, Wake Forest University School of Medicine, Winston-Salem, NC
| | - William Valdar
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC
| | - Richard Mott
- Genetics Institute, University College London, London, U.K
| | - Leah C. Solberg Woods
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| |
Collapse
|
5
|
Morón-García O, Garzón-Martínez GA, Martínez-Martín MJP, Brook J, Corke FMK, Doonan JH, Camargo Rodríguez AV. Genetic architecture of variation in Arabidopsis thaliana rosettes. PLoS One 2022; 17:e0263985. [PMID: 35171969 PMCID: PMC8849614 DOI: 10.1371/journal.pone.0263985] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Accepted: 02/01/2022] [Indexed: 12/04/2022] Open
Abstract
Rosette morphology across Arabidopsis accessions exhibits considerable variation. Here we report a high-throughput phenotyping approach based on automatic image analysis to quantify rosette shape and dissect the underlying genetic architecture. Shape measurements of the rosettes in a core set of Recombinant Inbred Lines from an advanced mapping population (Multiparent Advanced Generation Inter-Cross or MAGIC) derived from inter-crossing 19 natural accessions. Image acquisition and analysis was scaled to extract geometric descriptors from time stamped images of growing rosettes. Shape analyses revealed heritable morphological variation at early juvenile stages and QTL mapping resulted in over 116 chromosomal regions associated with trait variation within the population. Many QTL linked to variation in shape were located near genes related to hormonal signalling and signal transduction pathways while others are involved in shade avoidance and transition to flowering. Our results suggest rosette shape arises from modular integration of sub-organ morphologies and can be considered a functional trait subjected to selective pressures of subsequent morphological traits. On an applied aspect, QTLs found will be candidates for further research on plant architecture.
Collapse
Affiliation(s)
- Odín Morón-García
- The National Plant Phenomics Centre, Institute of Biological, Rural and Environmental Sciences (IBERS), Aberystwyth University, Aberystwyth, United Kingdom
| | - Gina A. Garzón-Martínez
- The National Plant Phenomics Centre, Institute of Biological, Rural and Environmental Sciences (IBERS), Aberystwyth University, Aberystwyth, United Kingdom
| | - M. J. Pilar Martínez-Martín
- The National Plant Phenomics Centre, Institute of Biological, Rural and Environmental Sciences (IBERS), Aberystwyth University, Aberystwyth, United Kingdom
| | - Jason Brook
- The National Plant Phenomics Centre, Institute of Biological, Rural and Environmental Sciences (IBERS), Aberystwyth University, Aberystwyth, United Kingdom
| | - Fiona M. K. Corke
- The National Plant Phenomics Centre, Institute of Biological, Rural and Environmental Sciences (IBERS), Aberystwyth University, Aberystwyth, United Kingdom
| | - John H. Doonan
- The National Plant Phenomics Centre, Institute of Biological, Rural and Environmental Sciences (IBERS), Aberystwyth University, Aberystwyth, United Kingdom
- * E-mail: (AVCR); (JHD)
| | - Anyela V. Camargo Rodríguez
- The National Plant Phenomics Centre, Institute of Biological, Rural and Environmental Sciences (IBERS), Aberystwyth University, Aberystwyth, United Kingdom
- * E-mail: (AVCR); (JHD)
| |
Collapse
|
6
|
Amadeu RR, Muñoz PR, Zheng C, Endelman JB. QTL mapping in outbred tetraploid (and diploid) diallel populations. Genetics 2021; 219:iyab124. [PMID: 34740237 PMCID: PMC8570786 DOI: 10.1093/genetics/iyab124] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Accepted: 07/24/2021] [Indexed: 11/14/2022] Open
Abstract
Over the last decade, multiparental populations have become a mainstay of genetics research in diploid species. Our goal was to extend this paradigm to autotetraploids by developing software for quantitative trait locus (QTL) mapping in connected F1 populations derived from a set of shared parents. For QTL discovery, phenotypes are regressed on the dosage of parental haplotypes to estimate additive effects. Statistical properties of the model were explored by simulating half-diallel diploid and tetraploid populations with different population sizes and numbers of parents. Across scenarios, the number of progeny per parental haplotype (pph) largely determined the statistical power for QTL detection and accuracy of the estimated haplotype effects. Multiallelic QTL with heritability 0.2 were detected with 90% probability at 25 pph and genome-wide significance level 0.05, and the additive haplotype effects were estimated with over 90% accuracy. Following QTL discovery, the software enables a comparison of models with multiple QTL and nonadditive effects. To illustrate, we analyzed potato tuber shape in a half-diallel population with three tetraploid parents. A well-known QTL on chromosome 10 was detected, for which the inclusion of digenic dominance lowered the Deviance Information Criterion (DIC) by 17 points compared to the additive model. The final model also contained a minor QTL on chromosome 1, but higher-order dominance and epistatic effects were excluded based on the DIC. In terms of practical impacts, the software is already being used to select offspring based on the effect and dosage of particular haplotypes in breeding programs.
Collapse
Affiliation(s)
- Rodrigo R Amadeu
- Horticultural Sciences Department, University of Florida, Gainesville, FL 32611, USA
| | - Patricio R Muñoz
- Horticultural Sciences Department, University of Florida, Gainesville, FL 32611, USA
| | - Chaozhi Zheng
- Biometris, Wageningen University and Research, 6708 PB Wageningen, The Netherlands
| | - Jeffrey B Endelman
- Department of Horticulture, University of Wisconsin, Madison, WI 53706, USA
| |
Collapse
|
7
|
Mosedale M, Cai Y, Eaddy JS, Kirby PJ, Wolenski FS, Dragan Y, Valdar W. Human-relevant mechanisms and risk factors for TAK-875-Induced liver injury identified via a gene pathway-based approach in Collaborative Cross mice. Toxicology 2021; 461:152902. [PMID: 34418498 PMCID: PMC8936092 DOI: 10.1016/j.tox.2021.152902] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Revised: 08/05/2021] [Accepted: 08/16/2021] [Indexed: 10/20/2022]
Abstract
Development of TAK-875 was discontinued when a small number of serious drug-induced liver injury (DILI) cases were observed in Phase 3 clinical trials. Subsequent studies have identified hepatocellular oxidative stress, mitochondrial dysfunction, altered bile acid homeostasis, and immune response as mechanisms of TAK-875 DILI and the contribution of genetic risk factors in oxidative response and mitochondrial pathways to the toxicity susceptibility observed in patients. We tested the hypothesis that a novel preclinical approach based on gene pathway analysis in the livers of Collaborative Cross mice could be used to identify human-relevant mechanisms of toxicity and genetic risk factors at the level of the hepatocyte as reported in a human genome-wide association study. Eight (8) male mice (4 matched pairs) from each of 45 Collaborative Cross lines were treated with a single oral (gavage) dose of either vehicle or 600 mg/kg TAK-875. As expected, liver injury was not detected histologically and few changes in plasma biomarkers of hepatotoxicity were observed. However, gene expression profiling in the liver identified hundreds of transcripts responsive to TAK-875 treatment across all strains reflecting alterations in immune response and bile acid homeostasis and the interaction of treatment and strain reflecting oxidative stress and mitochondrial dysfunction. Fold-change expression values were then used to develop pathway-based phenotypes for genetic mapping which identified candidate risk factor genes for TAK-875 toxicity susceptibility at the level of the hepatocyte. Taken together, these findings support our hypothesis that a gene pathway-based approach using Collaborative Cross mice could inform sensitive strains, human-relevant mechanisms of toxicity, and genetic risk factors for TAK-875 DILI. This novel preclinical approach may be helpful in understanding, predicting, and ultimately preventing clinical DILI for other drugs.
Collapse
Affiliation(s)
- Merrie Mosedale
- Division of Pharmacotherapy and Experimental Therapeutics, UNC Eshelman School of Pharmacy, Chapel Hill, NC, 27599, United States.
| | - Yanwei Cai
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, United States.
| | - J Scott Eaddy
- Division of Pharmacotherapy and Experimental Therapeutics, UNC Eshelman School of Pharmacy, Chapel Hill, NC, 27599, United States.
| | - Patrick J Kirby
- Takeda Pharmaceuticals International Co., Cambridge, MA, 02139, United States.
| | - Francis S Wolenski
- Takeda Pharmaceuticals International Co., Cambridge, MA, 02139, United States.
| | - Yvonne Dragan
- Takeda Pharmaceuticals International Co., Cambridge, MA, 02139, United States.
| | - William Valdar
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, United States; Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, United States.
| |
Collapse
|
8
|
Holcomb D, Hamasaki-Katagiri N, Laurie K, Katneni U, Kames J, Alexaki A, Bar H, Kimchi-Sarfaty C. New approaches to predict the effect of co-occurring variants on protein characteristics. Am J Hum Genet 2021; 108:1502-1511. [PMID: 34256028 DOI: 10.1016/j.ajhg.2021.06.011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Accepted: 06/14/2021] [Indexed: 12/14/2022] Open
Abstract
Predicting the effect of a mutated gene before the onset of symptoms of genetic diseases would greatly facilitate diagnosis and potentiate early intervention. There have been myriad attempts to predict the effects of single-nucleotide variants. However, the applicability of these efforts does not scale to co-occurring variants. Furthermore, an increasing number of protein therapeutics contain co-occurring nucleotide variations, adding uncertainty during development to the safety and efficiency of these drugs. Co-occurring nucleotide variants may often have synergistic, additive, or antagonistic effects on protein attributes, further complicating the task of outcome prediction. We tested four models based on the cooperative and antagonistic effects of co-occurring variants to predict pathogenicity and effectiveness of protein therapeutics. A total of 30 attributes, including amino acid and nucleotide features, as well as existing single-variant effect prediction tools, were considered on the basis of previous studies on single-nucleotide variants. Importantly, the effects of synonymous variants, often seen in protein therapeutics, were also included in our models. We used 12 datasets of people with monogenic diseases and controls with co-occurring genetic variants to evaluate the accuracy of our models, accomplishing a degree of accuracy comparable to that of prediction tools for single-nucleotide variants. More importantly, our framework is generalizable to new, well-curated datasets of monogenic diseases and new variant scoring tools. This approach successfully assists in addressing the challenging task of predicting the effect of co-occurring variants on pathogenicity and protein effectiveness and is applicable for a wide range of protein therapeutics and genetic diseases.
Collapse
|
9
|
Keele GR, Prokop JW, He H, Holl K, Littrell J, Deal AW, Kim Y, Kyle PB, Attipoe E, Johnson AC, Uhl KL, Sirpilla OL, Jahanbakhsh S, Robinson M, Levy S, Valdar W, Garrett MR, Solberg Woods LC. Sept8/SEPTIN8 involvement in cellular structure and kidney damage is identified by genetic mapping and a novel human tubule hypoxic model. Sci Rep 2021; 11:2071. [PMID: 33483609 PMCID: PMC7822875 DOI: 10.1038/s41598-021-81550-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Accepted: 01/05/2021] [Indexed: 01/29/2023] Open
Abstract
Chronic kidney disease (CKD), which can ultimately progress to kidney failure, is influenced by genetics and the environment. Genes identified in human genome wide association studies (GWAS) explain only a small proportion of the heritable variation and lack functional validation, indicating the need for additional model systems. Outbred heterogeneous stock (HS) rats have been used for genetic fine-mapping of complex traits, but have not previously been used for CKD traits. We performed GWAS for urinary protein excretion (UPE) and CKD related serum biochemistries in 245 male HS rats. Quantitative trait loci (QTL) were identified using a linear mixed effect model that tested for association with imputed genotypes. Candidate genes were identified using bioinformatics tools and targeted RNAseq followed by testing in a novel in vitro model of human tubule, hypoxia-induced damage. We identified two QTL for UPE and five for serum biochemistries. Protein modeling identified a missense variant within Septin 8 (Sept8) as a candidate for UPE. Sept8/SEPTIN8 expression increased in HS rats with elevated UPE and tubulointerstitial injury and in the in vitro hypoxia model. SEPTIN8 is detected within proximal tubule cells in human kidney samples and localizes with acetyl-alpha tubulin in the culture system. After hypoxia, SEPTIN8 staining becomes diffuse and appears to relocalize with actin. These data suggest a role of SEPTIN8 in cellular organization and structure in response to environmental stress. This study demonstrates that integration of a rat genetic model with an environmentally induced tubule damage system identifies Sept8/SEPTIN8 and informs novel aspects of the complex gene by environmental interactions contributing to CKD risk.
Collapse
Affiliation(s)
| | - Jeremy W Prokop
- HudsonAlpha Institute, Huntsville, AL, USA
- Department of Pediatrics and Human Development, Department of Pharmacology, Michigan State University, Grand Rapids, MI, USA
| | - Hong He
- Departments of Pediatrics and Physiology, Medical College of Wisconsin, Milwaukee, WI, USA
| | - Katie Holl
- Departments of Pediatrics and Physiology, Medical College of Wisconsin, Milwaukee, WI, USA
| | - John Littrell
- Departments of Pediatrics and Physiology, Medical College of Wisconsin, Milwaukee, WI, USA
| | - Aaron W Deal
- Department of Internal Medicine, Wake Forest School of Medicine, Winston-Salem, NC, USA
| | - Yunjung Kim
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Patrick B Kyle
- Department of Pharmacology, Medicine (Nephrology), Pediatrics (Genetics), University of Mississippi Medical Center, Jackson, MS, USA
- Department of Pathology, University of Mississippi Medical Center, Jackson, MS, USA
| | - Esinam Attipoe
- Department of Pharmacology, Medicine (Nephrology), Pediatrics (Genetics), University of Mississippi Medical Center, Jackson, MS, USA
| | - Ashley C Johnson
- Department of Pharmacology, Medicine (Nephrology), Pediatrics (Genetics), University of Mississippi Medical Center, Jackson, MS, USA
| | - Katie L Uhl
- Department of Pediatrics and Human Development, Department of Pharmacology, Michigan State University, Grand Rapids, MI, USA
| | - Olivia L Sirpilla
- Department of Pediatrics and Human Development, Department of Pharmacology, Michigan State University, Grand Rapids, MI, USA
| | - Seyedehameneh Jahanbakhsh
- Department of Pediatrics and Human Development, Department of Pharmacology, Michigan State University, Grand Rapids, MI, USA
| | | | - Shawn Levy
- HudsonAlpha Institute, Huntsville, AL, USA
| | - William Valdar
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
- Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
| | - Michael R Garrett
- Department of Pharmacology, Medicine (Nephrology), Pediatrics (Genetics), University of Mississippi Medical Center, Jackson, MS, USA
| | - Leah C Solberg Woods
- Department of Internal Medicine, Wake Forest School of Medicine, Winston-Salem, NC, USA.
| |
Collapse
|
10
|
Crouse WL, Kelada SNP, Valdar W. Inferring the Allelic Series at QTL in Multiparental Populations. Genetics 2020; 216:957-983. [PMID: 33082282 PMCID: PMC7768242 DOI: 10.1534/genetics.120.303393] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2020] [Accepted: 10/12/2020] [Indexed: 12/25/2022] Open
Abstract
Multiparental populations (MPPs) are experimental populations in which the genome of every individual is a mosaic of known founder haplotypes. These populations are useful for detecting quantitative trait loci (QTL) because tests of association can leverage inferred founder haplotype descent. It is difficult, however, to determine how haplotypes at a locus group into distinct functional alleles, termed the allelic series. The allelic series is important because it provides information about the number of causal variants at a QTL and their combined effects. In this study, we introduce a fully Bayesian model selection framework for inferring the allelic series. This framework accounts for sources of uncertainty found in typical MPPs, including the number and composition of functional alleles. Our prior distribution for the allelic series is based on the Chinese restaurant process, a relative of the Dirichlet process, and we leverage its connection to the coalescent to introduce additional prior information about haplotype relatedness via a phylogenetic tree. We evaluate our approach via simulation and apply it to QTL from two MPPs: the Collaborative Cross (CC) and the Drosophila Synthetic Population Resource (DSPR). We find that, although posterior inference of the exact allelic series is often uncertain, we are able to distinguish biallelic QTL from more complex multiallelic cases. Additionally, our allele-based approach improves haplotype effect estimation when the true number of functional alleles is small. Our method, Tree-Based Inference of Multiallelism via Bayesian Regression (TIMBR), provides new insight into the genetic architecture of QTL in MPPs.
Collapse
Affiliation(s)
- Wesley L Crouse
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina, Chapel Hill, North Carolina 27599
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina 27599
| | - Samir N P Kelada
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina 27599
- Marsico Lung Institute, University of North Carolina, Chapel Hill, North Carolina 27599
| | - William Valdar
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina 27599
- Lineberger Comprehensive Cancer Center, University of North Carolina, Chapel Hill, North Carolina 27599
| |
Collapse
|
11
|
Hsiao K, Noble C, Pitman W, Yadav N, Kumar S, Keele GR, Terceros A, Kanke M, Conniff T, Cheleuitte-Nieves C, Tolwani R, Sethupathy P, Rajasethupathy P. A Thalamic Orphan Receptor Drives Variability in Short-Term Memory. Cell 2020; 183:522-536.e19. [PMID: 32997977 DOI: 10.1016/j.cell.2020.09.011] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Revised: 06/09/2020] [Accepted: 09/01/2020] [Indexed: 12/31/2022]
Abstract
Working memory is a form of short-term memory that involves maintaining and updating task-relevant information toward goal-directed pursuits. Classical models posit persistent activity in prefrontal cortex (PFC) as a primary neural correlate, but emerging views suggest additional mechanisms may exist. We screened ∼200 genetically diverse mice on a working memory task and identified a genetic locus on chromosome 5 that contributes to a substantial proportion (17%) of the phenotypic variance. Within the locus, we identified a gene encoding an orphan G-protein-coupled receptor, Gpr12, which is sufficient to drive substantial and bidirectional changes in working memory. Molecular, cellular, and imaging studies revealed that Gpr12 enables high thalamus-PFC synchrony to support memory maintenance and choice accuracy. These findings identify an orphan receptor as a potent modifier of short-term memory and supplement classical PFC-based models with an emerging thalamus-centric framework for the mechanistic understanding of working memory.
Collapse
Affiliation(s)
- Kuangfu Hsiao
- Laboratory of Neural Dynamics & Cognition, The Rockefeller University, New York, NY 10065, USA
| | - Chelsea Noble
- Laboratory of Neural Dynamics & Cognition, The Rockefeller University, New York, NY 10065, USA
| | - Wendy Pitman
- Department of Biomedical Sciences, Cornell University, Ithaca, NY 14853, USA
| | - Nakul Yadav
- Laboratory of Neural Dynamics & Cognition, The Rockefeller University, New York, NY 10065, USA
| | - Suraj Kumar
- Laboratory of Neural Dynamics & Cognition, The Rockefeller University, New York, NY 10065, USA
| | | | - Andrea Terceros
- Laboratory of Neural Dynamics & Cognition, The Rockefeller University, New York, NY 10065, USA
| | - Matt Kanke
- Department of Biomedical Sciences, Cornell University, Ithaca, NY 14853, USA
| | - Tara Conniff
- Laboratory of Neural Dynamics & Cognition, The Rockefeller University, New York, NY 10065, USA
| | | | - Ravi Tolwani
- Comparative Bioscience Center, The Rockefeller University, New York, NY 10065, USA
| | - Praveen Sethupathy
- Department of Biomedical Sciences, Cornell University, Ithaca, NY 14853, USA.
| | | |
Collapse
|
12
|
Mosedale M, Cai Y, Eaddy JS, Corty RW, Nautiyal M, Watkins PB, Valdar W. Identification of Candidate Risk Factor Genes for Human Idelalisib Toxicity Using a Collaborative Cross Approach. Toxicol Sci 2020; 172:265-278. [PMID: 31501888 DOI: 10.1093/toxsci/kfz199] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Idelalisib is a phosphatidylinositol 3-kinase inhibitor highly selective for the delta isoform that has shown good efficacy in treating chronic lymphocytic leukemia and follicular lymphoma. In clinical trials, however, idelalisib was associated with rare, but potentially serious liver and lung toxicities. In this study, we used the Collaborative Cross (CC) mouse population to identify genetic factors associated with the drug response that may inform risk management strategies for idelalisib in humans. Eight male mice (4 matched pairs) from 50 CC lines were treated once daily for 14 days by oral gavage with either vehicle or idelalisib at a dose selected to achieve clinically relevant peak plasma concentrations (150 mg/kg/day). The drug was well tolerated across all CC lines, and there were no observations of overt liver injury. Differences across CC lines were seen in drug concentration in plasma samples collected at the approximate Tmax on study Days 1, 7, and 14. There were also small but statistically significant treatment-induced alterations in plasma total bile acids and microRNA-122, and these may indicate early hepatocellular stress required for immune-mediated hepatotoxicity in humans. Idelalisib treatment further induced significant elevations in the total cell count of terminal bronchoalveolar lavage fluid, which may be analogous to pneumonitis observed in the clinic. Genetic mapping identified loci associated with interim plasma idelalisib concentration and the other 3 treatment-related endpoints. Thirteen priority candidate quantitative trait genes identified in CC mice may now guide interrogation of risk factors for adverse drug responses associated with idelalisib in humans.
Collapse
Affiliation(s)
- Merrie Mosedale
- Institute for Drug Safety Sciences, University of North Carolina at Chapel Hill, Research Triangle Park, North Carolina 27709.,Division of Pharmacotherapy and Experimental Therapeutics, UNC Eshelman School of Pharmacy, Chapel Hill, North Carolina 27599
| | - Yanwei Cai
- Institute for Drug Safety Sciences, University of North Carolina at Chapel Hill, Research Triangle Park, North Carolina 27709.,Department of Genetics
| | - John Scott Eaddy
- Institute for Drug Safety Sciences, University of North Carolina at Chapel Hill, Research Triangle Park, North Carolina 27709
| | | | - Manisha Nautiyal
- Institute for Drug Safety Sciences, University of North Carolina at Chapel Hill, Research Triangle Park, North Carolina 27709
| | - Paul B Watkins
- Institute for Drug Safety Sciences, University of North Carolina at Chapel Hill, Research Triangle Park, North Carolina 27709.,Division of Pharmacotherapy and Experimental Therapeutics, UNC Eshelman School of Pharmacy, Chapel Hill, North Carolina 27599
| | - William Valdar
- Department of Genetics.,Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599
| |
Collapse
|
13
|
Keele GR, Quach BC, Israel JW, Chappell GA, Lewis L, Safi A, Simon JM, Cotney P, Crawford GE, Valdar W, Rusyn I, Furey TS. Integrative QTL analysis of gene expression and chromatin accessibility identifies multi-tissue patterns of genetic regulation. PLoS Genet 2020; 16:e1008537. [PMID: 31961859 PMCID: PMC7010298 DOI: 10.1371/journal.pgen.1008537] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2019] [Revised: 02/10/2020] [Accepted: 11/23/2019] [Indexed: 01/08/2023] Open
Abstract
Gene transcription profiles across tissues are largely defined by the activity of regulatory elements, most of which correspond to regions of accessible chromatin. Regulatory element activity is in turn modulated by genetic variation, resulting in variable transcription rates across individuals. The interplay of these factors, however, is poorly understood. Here we characterize expression and chromatin state dynamics across three tissues-liver, lung, and kidney-in 47 strains of the Collaborative Cross (CC) mouse population, examining the regulation of these dynamics by expression quantitative trait loci (eQTL) and chromatin QTL (cQTL). QTL whose allelic effects were consistent across tissues were detected for 1,101 genes and 133 chromatin regions. Also detected were eQTL and cQTL whose allelic effects differed across tissues, including local-eQTL for Pik3c2g detected in all three tissues but with distinct allelic effects. Leveraging overlapping measurements of gene expression and chromatin accessibility on the same mice from multiple tissues, we used mediation analysis to identify chromatin and gene expression intermediates of eQTL effects. Based on QTL and mediation analyses over multiple tissues, we propose a causal model for the distal genetic regulation of Akr1e1, a gene involved in glycogen metabolism, through the zinc finger transcription factor Zfp985 and chromatin intermediates. This analysis demonstrates the complexity of transcriptional and chromatin dynamics and their regulation over multiple tissues, as well as the value of the CC and related genetic resource populations for identifying specific regulatory mechanisms within cells and tissues.
Collapse
Affiliation(s)
- Gregory R. Keele
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- The Jackson Laboratory, Bar Harbor, Maine, United States of America
| | - Bryan C. Quach
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- Center for Omics Discovery and Epidemiology, Research Triangle Institute (RTI) International, Research Triangle Park, North Carolina, United States of America
| | - Jennifer W. Israel
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Grace A. Chappell
- Department of Veterinary Integrative Biosciences, Texas A&M University, College Station, Texas, United States of America
| | - Lauren Lewis
- Department of Veterinary Integrative Biosciences, Texas A&M University, College Station, Texas, United States of America
| | - Alexias Safi
- Department of Pediatrics, Duke University, Durham, North Carolina, United States of America
- Center for Genomic and Computational Biology, Duke University, Durham, North Carolina, United States of America
| | - Jeremy M. Simon
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Paul Cotney
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Gregory E. Crawford
- Department of Pediatrics, Duke University, Durham, North Carolina, United States of America
- Center for Genomic and Computational Biology, Duke University, Durham, North Carolina, United States of America
| | - William Valdar
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Ivan Rusyn
- Department of Veterinary Integrative Biosciences, Texas A&M University, College Station, Texas, United States of America
| | - Terrence S. Furey
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- Department of Biology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| |
Collapse
|
14
|
Abstract
In this chapter we will review both the rationale and experimental design for using Heterogeneous Stock (HS) populations for fine-mapping of complex traits in mice and rats. We define an HS as an outbred population derived from an intercross between two or more inbred strains. HS have been used to perform genome-wide association studies (GWAS) for multiple behavioral, physiological, and gene expression traits. GWAS using HS require four key steps, which we review: selection of an appropriate HS population, phenotyping, genotyping, and statistical analysis. We provide advice on the selection of an HS, comment on key issues related to phenotyping, discuss genotyping methods relevant to these populations, and describe statistical genetic analyses that are applicable to genetic analyses that use HS.
Collapse
|
15
|
He T, Hill CB, Angessa TT, Zhang XQ, Chen K, Moody D, Telfer P, Westcott S, Li C. Gene-set association and epistatic analyses reveal complex gene interaction networks affecting flowering time in a worldwide barley collection. JOURNAL OF EXPERIMENTAL BOTANY 2019; 70:5603-5616. [PMID: 31504706 PMCID: PMC6812734 DOI: 10.1093/jxb/erz332] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/05/2019] [Accepted: 08/13/2019] [Indexed: 05/10/2023]
Abstract
Single-marker genome-wide association studies (GWAS) have successfully detected associations between single nucleotide polymorphisms (SNPs) and agronomic traits such as flowering time and grain yield in barley. However, the analysis of individual SNPs can only account for a small proportion of genetic variation, and can only provide limited knowledge on gene network interactions. Gene-based GWAS approaches provide enormous opportunity both to combine genetic information and to examine interactions among genetic variants. Here, we revisited a previously published phenotypic and genotypic data set of 895 barley varieties grown in two years at four different field locations in Australia. We employed statistical models to examine gene-phenotype associations, as well as two-way epistasis analyses to increase the capability to find novel genes that have significant roles in controlling flowering time in barley. Genetic associations were tested between flowering time and corresponding genotypes of 174 putative flowering time-related genes. Gene-phenotype association analysis detected 113 genes associated with flowering time in barley, demonstrating the unprecedented power of gene-based analysis. Subsequent two-way epistasis analysis revealed 19 pairs of gene×gene interactions involved in controlling flowering time. Our study demonstrates that gene-based association approaches can provide higher capacity for future crop improvement to increase crop performance and adaptation to different environments.
Collapse
Affiliation(s)
- Tianhua He
- Western Barley Genetics Alliance, Western Australian State Agricultural Biotechnology Centre, College of Science, Health, Engineering and Education, Murdoch University, Murdoch, WA, Australia
| | - Camilla Beate Hill
- Western Barley Genetics Alliance, Western Australian State Agricultural Biotechnology Centre, College of Science, Health, Engineering and Education, Murdoch University, Murdoch, WA, Australia
| | - Tefera Tolera Angessa
- Western Barley Genetics Alliance, Western Australian State Agricultural Biotechnology Centre, College of Science, Health, Engineering and Education, Murdoch University, Murdoch, WA, Australia
| | - Xiao-Qi Zhang
- Western Barley Genetics Alliance, Western Australian State Agricultural Biotechnology Centre, College of Science, Health, Engineering and Education, Murdoch University, Murdoch, WA, Australia
| | - Kefei Chen
- SAGI-WEST, Faculty of Science and Engineering, Curtin University, Bentley, WA, Australia
| | | | - Paul Telfer
- Australian Grain Technologies Pty Ltd (AGT), SA, Australia
| | - Sharon Westcott
- Agriculture and Food, Department of Primary Industries and Regional Development, South Perth, WA, Australia
| | - Chengdao Li
- Western Barley Genetics Alliance, Western Australian State Agricultural Biotechnology Centre, College of Science, Health, Engineering and Education, Murdoch University, Murdoch, WA, Australia
- Agriculture and Food, Department of Primary Industries and Regional Development, South Perth, WA, Australia
- Hubei Collaborative Innovation Centre for Grain Industry, Yangtze University, Hubei Jingzhou, China
- Correspondence:
| |
Collapse
|
16
|
Abstract
The Collaborative Cross (CC) is a mouse genetic reference population whose range of applications includes quantitative trait loci (QTL) mapping. The design of a CC QTL mapping study involves multiple decisions, including which and how many strains to use, and how many replicates per strain to phenotype, all viewed within the context of hypothesized QTL architecture. Until now, these decisions have been informed largely by early power analyses that were based on simulated, hypothetical CC genomes. Now that more than 50 CC strains are available and more than 70 CC genomes have been observed, it is possible to characterize power based on realized CC genomes. We report power analyses from extensive simulations and examine several key considerations: 1) the number of strains and biological replicates, 2) the QTL effect size, 3) the presence of population structure, and 4) the distribution of functionally distinct alleles among the founder strains at the QTL. We also provide general power estimates to aide in the design of future experiments. All analyses were conducted with our R package, SPARCC (Simulated Power Analysis in the Realized Collaborative Cross), developed for performing either large scale power analyses or those tailored to particular CC experiments.
Collapse
|
17
|
Mosedale M. Mouse Population-Based Approaches to Investigate Adverse Drug Reactions. Drug Metab Dispos 2018; 46:1787-1795. [DOI: 10.1124/dmd.118.082834] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Accepted: 07/06/2018] [Indexed: 01/19/2023] Open
|
18
|
Haplotype-Based Genome-Wide Prediction Models Exploit Local Epistatic Interactions Among Markers. G3-GENES GENOMES GENETICS 2018; 8:1687-1699. [PMID: 29549092 PMCID: PMC5940160 DOI: 10.1534/g3.117.300548] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Genome-wide prediction approaches represent versatile tools for the analysis and prediction of complex traits. Mostly they rely on marker-based information, but scenarios have been reported in which models capitalizing on closely-linked markers that were combined into haplotypes outperformed marker-based models. Detailed comparisons were undertaken to reveal under which circumstances haplotype-based genome-wide prediction models are superior to marker-based models. Specifically, it was of interest to analyze whether and how haplotype-based models may take local epistatic effects between markers into account. Assuming that populations consisted of fully homozygous individuals, a marker-based model in which local epistatic effects inside haplotype blocks were exploited (LEGBLUP) was linearly transformable into a haplotype-based model (HGBLUP). This theoretical derivation formally revealed that haplotype-based genome-wide prediction models capitalize on local epistatic effects among markers. Simulation studies corroborated this finding. Due to its computational efficiency the HGBLUP model promises to be an interesting tool for studies in which ultra-high-density SNP data sets are studied. Applying the HGBLUP model to empirical data sets revealed higher prediction accuracies than for marker-based models for both traits studied using a mouse panel. In contrast, only a small subset of the traits analyzed in crop populations showed such a benefit. Cases in which higher prediction accuracies are observed for HGBLUP than for marker-based models are expected to be of immediate relevance for breeders, due to the tight linkage a beneficial haplotype will be preserved for many generations. In this respect the inheritance of local epistatic effects very much resembles the one of additive effects.
Collapse
|
19
|
Joint Analysis of Strain and Parent-of-Origin Effects for Recombinant Inbred Intercrosses Generated from Multiparent Populations with the Collaborative Cross as an Example. G3-GENES GENOMES GENETICS 2018; 8:599-605. [PMID: 29255115 PMCID: PMC5919741 DOI: 10.1534/g3.117.300483] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Multiparent populations (MPP) have become popular resources for complex trait mapping because of their wider allelic diversity and larger population size compared with traditional two-way recombinant inbred (RI) strains. In mice, the collaborative cross (CC) is one of the most popular MPP and is derived from eight genetically diverse inbred founder strains. The strategy of generating RI intercrosses (RIX) from MPP in general and from the CC in particular can produce a large number of completely reproducible heterozygote genomes that better represent the (outbred) human population. Since both maternal and paternal haplotypes of each RIX are readily available, RIX is a powerful resource for studying both standing genetic and epigenetic variations of complex traits, in particular, the parent-of-origin (PoO) effects, which are important contributors to many complex traits. Furthermore, most complex traits are affected by >1 genes, where multiple quantitative trait locus mapping could be more advantageous. In this paper, for MPP-RIX data but taking CC-RIX as a working example, we propose a general Bayesian variable selection procedure to simultaneously search for multiple genes with founder allelic effects and PoO effects. The proposed model respects the complex relationship among RIX samples, and the performance of the proposed method is examined by extensive simulations.
Collapse
|
20
|
Bayesian Diallel Analysis Reveals Mx1-Dependent and Mx1-Independent Effects on Response to Influenza A Virus in Mice. G3-GENES GENOMES GENETICS 2018; 8:427-445. [PMID: 29187420 PMCID: PMC5919740 DOI: 10.1534/g3.117.300438] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Influenza A virus (IAV) is a respiratory pathogen that causes substantial morbidity and mortality during both seasonal and pandemic outbreaks. Infection outcomes in unexposed populations are affected by host genetics, but the host genetic architecture is not well understood. Here, we obtain a broad view of how heritable factors affect a mouse model of response to IAV infection using an 8 × 8 diallel of the eight inbred founder strains of the Collaborative Cross (CC). Expanding on a prior statistical framework for modeling treatment response in diallels, we explore how a range of heritable effects modify acute host response to IAV through 4 d postinfection. Heritable effects in aggregate explained ∼57% of the variance in IAV-induced weight loss. Much of this was attributable to a pattern of additive effects that became more prominent through day 4 postinfection and was consistent with previous reports of antiinfluenza myxovirus resistance 1 (Mx1) polymorphisms segregating between these strains; these additive effects largely recapitulated haplotype effects observed at the Mx1 locus in a previous study of the incipient CC, and are also replicated here in a CC recombinant intercross population. Genetic dominance of protective Mx1 haplotypes was observed to differ by subspecies of origin: relative to the domesticus null Mx1 allele, musculus acts dominantly whereas castaneus acts additively. After controlling for Mx1, heritable effects, though less distinct, accounted for ∼34% of the phenotypic variance. Implications for future mapping studies are discussed.
Collapse
|
21
|
Mosedale M, Kim Y, Brock WJ, Roth SE, Wiltshire T, Eaddy JS, Keele GR, Corty RW, Xie Y, Valdar W, Watkins PB. Editor's Highlight: Candidate Risk Factors and Mechanisms for Tolvaptan-Induced Liver Injury Are Identified Using a Collaborative Cross Approach. Toxicol Sci 2018; 156:438-454. [PMID: 28115652 DOI: 10.1093/toxsci/kfw269] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
Clinical trials of tolvaptan showed it to be a promising candidate for the treatment of Autosomal Dominant Polycystic Kidney Disease (ADPKD) but also revealed potential for idiosyncratic drug-induced liver injury (DILI) in this patient population. To identify risk factors and mechanisms underlying tolvaptan DILI, 8 mice in each of 45 strains of the genetically diverse Collaborative Cross (CC) mouse population were treated with a single oral dose of either tolvaptan or vehicle. Significant elevations in plasma alanine aminotransferase (ALT) were observed in tolvaptan-treated animals in 3 of the 45 strains. Genetic mapping coupled with transcriptomic analysis in the liver was used to identify several candidate susceptibility genes including epoxide hydrolase 2, interferon regulatory factor 3, and mitochondrial fission factor. Gene pathway analysis revealed that oxidative stress and immune response pathways were activated in response to tolvaptan treatment across all strains, but genes involved in regulation of bile acid homeostasis were most associated with tolvaptan-induced elevations in ALT. Secretory leukocyte peptidase inhibitor (Slpi) mRNA was also induced in the susceptible strains and was associated with increased plasma levels of Slpi protein, suggesting a potential serum marker for DILI susceptibility. In summary, tolvaptan induced signs of oxidative stress, mitochondrial dysfunction, and innate immune response in all strains, but variation in bile acid homeostasis was most associated with susceptibility to the liver response. This CC study has indicated potential mechanisms underlying tolvaptan DILI and biomarkers of susceptibility that may be useful in managing the risk of DILI in ADPKD patients.
Collapse
Affiliation(s)
- Merrie Mosedale
- Institute for Drug Safety Sciences, University of North Carolina at Chapel Hill, Research Triangle Park, North Carolina 27709.,Division of Pharmacotherapy and Experimental Therapeutics, UNC Eshelman School of Pharmacy, Chapel Hill, North Carolina 27599
| | - Yunjung Kim
- Institute for Drug Safety Sciences, University of North Carolina at Chapel Hill, Research Triangle Park, North Carolina 27709.,Department of Genetics, UNC School of Medicine, Chapel Hill, North Carolina 27599
| | - William J Brock
- Otsuka Pharmaceutical Development and Commercialization, Inc., Rockville, Maryland 20850.,Brock Scientific Consulting, Montgomery Village, Maryland 20886
| | - Sharin E Roth
- Otsuka Pharmaceutical Development and Commercialization, Inc., Rockville, Maryland 20850
| | - Tim Wiltshire
- Division of Pharmacotherapy and Experimental Therapeutics, UNC Eshelman School of Pharmacy, Chapel Hill, North Carolina 27599.,Department of Genetics, UNC School of Medicine, Chapel Hill, North Carolina 27599
| | - J Scott Eaddy
- Institute for Drug Safety Sciences, University of North Carolina at Chapel Hill, Research Triangle Park, North Carolina 27709.,Division of Pharmacotherapy and Experimental Therapeutics, UNC Eshelman School of Pharmacy, Chapel Hill, North Carolina 27599
| | - Gregory R Keele
- Department of Genetics, UNC School of Medicine, Chapel Hill, North Carolina 27599
| | - Robert W Corty
- Department of Genetics, UNC School of Medicine, Chapel Hill, North Carolina 27599
| | - Yuying Xie
- Department of Genetics, UNC School of Medicine, Chapel Hill, North Carolina 27599
| | - William Valdar
- Department of Genetics, UNC School of Medicine, Chapel Hill, North Carolina 27599.,Lineberger Comprehensive Cancer Center, Chapel Hill, North Carolina 27599
| | - Paul B Watkins
- Institute for Drug Safety Sciences, University of North Carolina at Chapel Hill, Research Triangle Park, North Carolina 27709.,Division of Pharmacotherapy and Experimental Therapeutics, UNC Eshelman School of Pharmacy, Chapel Hill, North Carolina 27599
| |
Collapse
|
22
|
Keele GR, Prokop JW, He H, Holl K, Littrell J, Deal A, Francic S, Cui L, Gatti DM, Broman KW, Tschannen M, Tsaih SW, Zagloul M, Kim Y, Baur B, Fox J, Robinson M, Levy S, Flister MJ, Mott R, Valdar W, Solberg Woods LC. Genetic Fine-Mapping and Identification of Candidate Genes and Variants for Adiposity Traits in Outbred Rats. Obesity (Silver Spring) 2018; 26:213-222. [PMID: 29193816 PMCID: PMC5740008 DOI: 10.1002/oby.22075] [Citation(s) in RCA: 49] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Revised: 10/20/2017] [Accepted: 10/21/2017] [Indexed: 12/20/2022]
Abstract
OBJECTIVE Obesity is a major risk factor for multiple diseases and is in part heritable, yet the majority of causative genetic variants that drive excessive adiposity remain unknown. Here, outbred heterogeneous stock (HS) rats were used in controlled environmental conditions to fine-map novel genetic modifiers of adiposity. METHODS Body weight and visceral fat pad weights were measured in male HS rats that were also genotyped genome-wide. Quantitative trait loci (QTL) were identified by genome-wide association of imputed single-nucleotide polymorphism (SNP) genotypes using a linear mixed effect model that accounts for unequal relatedness between the HS rats. Candidate genes were assessed by protein modeling and mediation analysis of expression for coding and noncoding variants, respectively. RESULTS HS rats exhibited large variation in adiposity traits, which were highly heritable and correlated with metabolic health. Fine-mapping of fat pad weight and body weight revealed three QTL and prioritized five candidate genes. Fat pad weight was associated with missense SNPs in Adcy3 and Prlhr and altered expression of Krtcap3 and Slc30a3, whereas Grid2 was identified as a candidate within the body weight locus. CONCLUSIONS These data demonstrate the power of HS rats for identification of known and novel heritable mediators of obesity traits.
Collapse
Affiliation(s)
- Gregory R. Keele
- Department of Genetics, University of North Carolina at Chapel Hill, NC
| | | | - Hong He
- Medical College of Wisconsin, Departments of Pediatrics and Physiology, Milwaukee, WI
| | - Katie Holl
- Medical College of Wisconsin, Departments of Pediatrics and Physiology, Milwaukee, WI
| | - John Littrell
- Medical College of Wisconsin, Departments of Pediatrics and Physiology, Milwaukee, WI
| | - Aaron Deal
- Wake Forest School of Medicine, Department of Internal Medicine, Winston Salem, NC
| | - Sanja Francic
- University College London Genetics Institute, London, UK
| | - Leilei Cui
- University College London Genetics Institute, London, UK
| | | | - Karl W. Broman
- Department of Biostatistics and Medical Informatics, University of Wisconsin, Madison, WI
| | - Michael Tschannen
- Medical College of Wisconsin, Departments of Pediatrics and Physiology, Milwaukee, WI
| | - Shirng-Wern Tsaih
- Medical College of Wisconsin, Departments of Pediatrics and Physiology, Milwaukee, WI
| | - Maie Zagloul
- Department of Genetics, University of North Carolina at Chapel Hill, NC
| | - Yunjung Kim
- Department of Genetics, University of North Carolina at Chapel Hill, NC
| | - Brittany Baur
- Medical College of Wisconsin, Departments of Pediatrics and Physiology, Milwaukee, WI
| | - Joseph Fox
- Medical College of Wisconsin, Departments of Pediatrics and Physiology, Milwaukee, WI
| | | | | | - Michael J. Flister
- Medical College of Wisconsin, Departments of Pediatrics and Physiology, Milwaukee, WI
| | - Richard Mott
- University College London Genetics Institute, London, UK
| | - William Valdar
- Department of Genetics, University of North Carolina at Chapel Hill, NC
| | - Leah C Solberg Woods
- Wake Forest School of Medicine, Department of Internal Medicine, Winston Salem, NC
| |
Collapse
|
23
|
Oreper D, Cai Y, Tarantino LM, de Villena FPM, Valdar W. Inbred Strain Variant Database (ISVdb): A Repository for Probabilistically Informed Sequence Differences Among the Collaborative Cross Strains and Their Founders. G3 (BETHESDA, MD.) 2017; 7:1623-1630. [PMID: 28592645 PMCID: PMC5473744 DOI: 10.1534/g3.117.041491] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/05/2017] [Accepted: 03/20/2017] [Indexed: 02/07/2023]
Abstract
The Collaborative Cross (CC) is a panel of recently established multiparental recombinant inbred mouse strains. For the CC, as for any multiparental population (MPP), effective experimental design and analysis benefit from detailed knowledge of the genetic differences between strains. Such differences can be directly determined by sequencing, but until now whole-genome sequencing was not publicly available for individual CC strains. An alternative and complementary approach is to infer genetic differences by combining two pieces of information: probabilistic estimates of the CC haplotype mosaic from a custom genotyping array, and probabilistic variant calls from sequencing of the CC founders. The computation for this inference, especially when performed genome-wide, can be intricate and time-consuming, requiring the researcher to generate nontrivial and potentially error-prone scripts. To provide standardized, easy-to-access CC sequence information, we have developed the Inbred Strain Variant Database (ISVdb). The ISVdb provides, for all the exonic variants from the Sanger Institute mouse sequencing dataset, direct sequence information for CC founders and, critically, the imputed sequence information for CC strains. Notably, the ISVdb also: (1) provides predicted variant consequence metadata; (2) allows rapid simulation of F1 populations; and (3) preserves imputation uncertainty, which will allow imputed data to be refined in the future as additional sequencing and genotyping data are collected. The ISVdb information is housed in an SQL database and is easily accessible through a custom online interface (http://isvdb.unc.edu), reducing the analytic burden on any researcher using the CC.
Collapse
Affiliation(s)
- Daniel Oreper
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina, Chapel Hill, North Carolina 27599-7265
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina 27599-7265
| | - Yanwei Cai
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina, Chapel Hill, North Carolina 27599-7265
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina 27599-7265
| | - Lisa M Tarantino
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina 27599-7265
- Division of Pharmacotherapy and Experimental Therapeutics, Eshelman School of Pharmacy University of North Carolina, Chapel Hill, North Carolina 27599-7265
| | - Fernando Pardo-Manuel de Villena
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina 27599-7265
- Lineberger Comprehensive Cancer Center
| | - William Valdar
- Department of Genetics, University of North Carolina, Chapel Hill, North Carolina 27599-7265
- Lineberger Comprehensive Cancer Center
| |
Collapse
|
24
|
Schughart K, Williams RW. The Collaborative Cross Resource for Systems Genetics Research of Infectious Diseases. Methods Mol Biol 2017; 1488:579-596. [PMID: 27933545 PMCID: PMC7120135 DOI: 10.1007/978-1-4939-6427-7_28] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
An increasing body of evidence highlights the role of host genetic variation in driving susceptibility to severe disease following pathogen infection. In order to fully appreciate the importance of host genetics on infection susceptibility and resulting disease, genetically variable experimental model systems should be employed. These systems allow for the identification, characterization, and mechanistic dissection of genetic variants that cause differential disease responses. Herein we discuss application of the Collaborative Cross (CC) panel of recombinant inbred strains to study viral pathogenesis, focusing on practical considerations for experimental design, assessment and analysis of disease responses within the CC, as well as some of the resources developed for the CC. Although the focus of this chapter is on viral pathogenesis, many of the methods presented within are applicable to studies of other pathogens, as well as to case-control designs in genetically diverse populations.
Collapse
Affiliation(s)
- Klaus Schughart
- Department of Infection Genetics, Helmholtz Centre for Infection Research & University of Veterinary Medicine Hannover, Braunschweig, Niedersachsen Germany
| | - Robert W. Williams
- Department of Microbiology, Immunology and Biochemistry, University of Tennessee Health Science Center, Memphis, Tennessee USA
| |
Collapse
|
25
|
A Random-Model Approach to QTL Mapping in Multiparent Advanced Generation Intercross (MAGIC) Populations. Genetics 2015; 202:471-86. [PMID: 26715662 DOI: 10.1534/genetics.115.179945] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2015] [Accepted: 12/15/2015] [Indexed: 11/18/2022] Open
Abstract
Most standard QTL mapping procedures apply to populations derived from the cross of two parents. QTL detected from such biparental populations are rarely relevant to breeding programs because of the narrow genetic basis: only two alleles are involved per locus. To improve the generality and applicability of mapping results, QTL should be detected using populations initiated from multiple parents, such as the multiparent advanced generation intercross (MAGIC) populations. The greatest challenges of QTL mapping in MAGIC populations come from multiple founder alleles and control of the genetic background information. We developed a random-model methodology by treating the founder effects of each locus as random effects following a normal distribution with a locus-specific variance. We also fit a polygenic effect to the model to control the genetic background. To improve the statistical power for a scanned marker, we release the marker effect absorbed by the polygene back to the model. In contrast to the fixed-model approach, we estimate and test the variance of each locus and scan the entire genome one locus at a time using likelihood-ratio test statistics. Simulation studies showed that this method can increase statistical power and reduce type I error compared with composite interval mapping (CIM) and multiparent whole-genome average interval mapping (MPWGAIM). We demonstrated the method using a public Arabidopsis thaliana MAGIC population and a mouse MAGIC population.
Collapse
|
26
|
Abstract
The Collaborative Cross (CC) was designed to facilitate rapid gene mapping and consists of hundreds of recombinant inbred lines descended from eight diverse inbred founder strains. A decade in production, it can now be applied to mapping projects. Here, we provide a proof of principle for rapid identification of major-effect genes using the CC. To do so, we chose coat color traits since the location and identity of many relevant genes are known. We ascertained in 110 CC lines six different coat phenotypes: albino, agouti, black, cinnamon, and chocolate coat colors and the white-belly trait. We developed a pipeline employing modifications of existing mapping tools suitable for analyzing the complex genetic architecture of the CC. Together with analysis of the founders’ genome sequences, mapping was successfully achieved with sufficient resolution to identify the causative genes for five traits. Anticipating the application of the CC to complex traits, we also developed strategies to detect interacting genes, testing joint effects of three loci. Our results illustrate the power of the CC and provide confidence that this resource can be applied to complex traits for detection of both qualitative and quantitative trait loci.
Collapse
|
27
|
Morgan AP, Welsh CE. Informatics resources for the Collaborative Cross and related mouse populations. Mamm Genome 2015; 26:521-39. [PMID: 26135136 DOI: 10.1007/s00335-015-9581-z] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2015] [Accepted: 06/23/2015] [Indexed: 02/05/2023]
Affiliation(s)
- Andrew P Morgan
- Department of Genetics, University of North Carolina, Chapel Hill, NC, USA
| | - Catherine E Welsh
- Department of Mathematics & Computer Science, Rhodes College, Memphis, TN, USA.
| |
Collapse
|
28
|
Quantitative trait locus mapping methods for diversity outbred mice. G3-GENES GENOMES GENETICS 2014; 4:1623-33. [PMID: 25237114 PMCID: PMC4169154 DOI: 10.1534/g3.114.013748] [Citation(s) in RCA: 150] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
Genetic mapping studies in the mouse and other model organisms are used to search for genes underlying complex phenotypes. Traditional genetic mapping studies that employ single-generation crosses have poor mapping resolution and limit discovery to loci that are polymorphic between the two parental strains. Multiparent outbreeding populations address these shortcomings by increasing the density of recombination events and introducing allelic variants from multiple founder strains. However, multiparent crosses present new analytical challenges and require specialized software to take full advantage of these benefits. Each animal in an outbreeding population is genetically unique and must be genotyped using a high-density marker set; regression models for mapping must accommodate multiple founder alleles, and complex breeding designs give rise to polygenic covariance among related animals that must be accounted for in mapping analysis. The Diversity Outbred (DO) mice combine the genetic diversity of eight founder strains in a multigenerational breeding design that has been maintained for >16 generations. The large population size and randomized mating ensure the long-term genetic stability of this population. We present a complete analytical pipeline for genetic mapping in DO mice, including algorithms for probabilistic reconstruction of founder haplotypes from genotyping array intensity data, and mapping methods that accommodate multiple founder haplotypes and account for relatedness among animals. Power analysis suggests that studies with as few as 200 DO mice can detect loci with large effects, but loci that account for <5% of trait variance may require a sample size of up to 1000 animals. The methods described here are implemented in the freely available R package DOQTL.
Collapse
|