1
|
Zheng Y, Lin C, Wang WJ, Wang L, Qian Y, Mao L, Li B, Lou L, Mao Y, Li N, Zheng J, Jiang N, He C, Wang Q, Zhou Q, Chen F, Jin F. Post-implantation analysis of genomic variations in the progeny from developing fetus to birth. Hum Genomics 2024; 18:79. [PMID: 39010135 PMCID: PMC11247737 DOI: 10.1186/s40246-024-00634-4] [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: 02/06/2024] [Accepted: 06/06/2024] [Indexed: 07/17/2024] Open
Abstract
The analysis of genomic variations in offspring after implantation has been infrequently studied. In this study, we aim to investigate the extent of de novo mutations in humans from developing fetus to birth. Using high-depth whole-genome sequencing, 443 parent-offspring trios were studied to compare the results of de novo mutations (DNMs) between different groups. The focus was on fetuses and newborns, with DNA samples obtained from the families' blood and the aspirated embryonic tissues subjected to deep sequencing. It was observed that the average number of total DNMs in the newborns group was 56.26 (54.17-58.35), which appeared to be lower than that the multifetal reduction group, which was 76.05 (69.70-82.40) (F = 2.42, P = 0.12). However, after adjusting for parental age and maternal pre-pregnancy body mass index (BMI), significant differences were found between the two groups. The analysis was further divided into single nucleotide variants (SNVs) and insertion/deletion of a small number of bases (indels), and it was discovered that the average number of de novo SNVs associated with the multifetal reduction group and the newborn group was 49.89 (45.59-54.20) and 51.09 (49.22-52.96), respectively. No significant differences were noted between the groups (F = 1.01, P = 0.32). However, a significant difference was observed for de novo indels, with a higher average number found in the multifetal reduction group compared to the newborn group (F = 194.17, P < 0.001). The average number of de novo indels among the multifetal reduction group and the newborn group was 26.26 (23.27-29.05) and 5.17 (4.82-5.52), respectively. To conclude, it has been observed that the quantity of de novo indels in the newborns experiences a significant decrease when compared to that in the aspirated embryonic tissues (7-9 weeks). This phenomenon is evident across all genomic regions, highlighting the adverse effects of de novo indels on the fetus and emphasizing the significance of embryonic implantation and intrauterine growth in human genetic selection mechanisms.
Collapse
Affiliation(s)
- Yingming Zheng
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
| | - Chuanping Lin
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
- Reproductive Medical Center, the Second Affiliated Hospital of Wenzhou Medical College and Yuying Children's hospital, Wenzhou, Zhejiang, 325027, China
| | | | - Liya Wang
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
| | - Yeqing Qian
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
| | - Luna Mao
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
| | - Baohua Li
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
| | - Lijun Lou
- Affiliated Dongyang Hospital of Wenzhou Medical University, Dongyang, Zhejiang, 322100, China
| | - Yuchan Mao
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
| | - Na Li
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
| | - Jiayong Zheng
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
| | - Nan Jiang
- Reproductive Medical Center, the First Affiliated Hospital, School of Medicine, Zhejiang University, Hangzhou, Zhejiang, 310003, China
| | - Chaying He
- Hangzhou Women's Hospital (Hangzhou Maternity and Child Health Care Hospital), Hangzhou, Zhejiang, 310008, China
| | - Qijing Wang
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China
| | - Qing Zhou
- BGI Research, Shenzhen, Guangdong, 518083, China
| | - Fang Chen
- BGI Research, Shenzhen, Guangdong, 518083, China
| | - Fan Jin
- Department of Reproductive Endocrinology, Key Laboratory of Reproductive Genetics of National Ministry of Education, Women's Reproductive Health Laboratory of Zhejiang Province, Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Hangzhou, Zhejiang, 310006, China.
| |
Collapse
|
2
|
Imamura M, Nakai R, Ohnuki M, Hamazaki Y, Tanabe H, Sato M, Harishima Y, Horikawa M, Watanabe M, Oota H, Nakagawa M, Suzuki S, Enard W. Generation of chimpanzee induced pluripotent stem cell lines for cross-species comparisons. In Vitro Cell Dev Biol Anim 2024; 60:544-554. [PMID: 38386235 DOI: 10.1007/s11626-024-00853-y] [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: 11/27/2023] [Accepted: 01/04/2024] [Indexed: 02/23/2024]
Abstract
As humans' closest living relatives, chimpanzees offer valuable insights into human evolution. However, technical and ethical limitations hinder investigations into the molecular and cellular foundations that distinguish chimpanzee and human traits. Recently, induced pluripotent stem cells (iPSCs) have emerged as a novel model for functional comparative studies and provided a non-invasive alternative for studying embryonic phenomena. In this study, we generated five new chimpanzee iPSC lines from peripheral blood cells and skin fibroblasts with SeV vectors carrying four reprogramming factors (human OCT3/4, SOX2, KLF4, and L-MYC) and characterized their pluripotency and differentiation potential. We also examined the expression of a human-specific non-coding RNA, HSTR1, which is predicted to be involved in human brain development. Our results show that the chimpanzee iPSCs possess pluripotent characteristics and can differentiate into various cell lineages. Moreover, we found that HSTR1 is expressed in human iPSCs and their neural derivatives but not in chimpanzee counterparts, supporting its possible role in human-specific brain development. As iPSCs are inherently variable due to genetic and epigenetic differences in donor cells or reprogramming procedures, it is essential to expand the number of chimpanzee iPSC lines to comprehensively capture the molecular and cellular properties representative of chimpanzees. Hence, our cells provide a valuable resource for investigating the function and regulation of human-specific transcripts such as HSTR1 and for understanding human evolution more generally.
Collapse
Affiliation(s)
- Masanori Imamura
- Molecular Biology Section, Center for the Evolutionary Origins of Human Behavior, Kyoto University, Inuyama, Aichi, 484-8506, Japan.
| | - Risako Nakai
- Molecular Biology Section, Center for the Evolutionary Origins of Human Behavior, Kyoto University, Inuyama, Aichi, 484-8506, Japan
- iPSC-Based Drug Discovery and Development Team, RIKEN BioResource Research Center, Soraku, Kyoto, 619-0237, Japan
- Center for iPS Cell Research and Application (CiRA), Kyoto University, Kyoto, 606-8507, Japan
| | - Mari Ohnuki
- Center for iPS Cell Research and Application (CiRA), Kyoto University, Kyoto, 606-8507, Japan
- Institute for the Advanced Study of Human Biology (WPI-ASHBi), Kyoto University, Kyoto, 606-8501, Japan
- Hakubi Center, Kyoto University, Kyoto, 606-8501, Japan
- Anthropology and Human Genomics, Faculty of Biology, Ludwig-Maximilians-Universität München, München, Germany
| | - Yusuke Hamazaki
- Molecular Biology Section, Center for the Evolutionary Origins of Human Behavior, Kyoto University, Inuyama, Aichi, 484-8506, Japan
| | - Hideyuki Tanabe
- Research Center for Integrative Evolutionary Science, SOKENDAI (The Graduate University for Advanced Studies), Hayama, 240-0193, Japan
| | - Momoka Sato
- Department of Agricultural and Life Sciences, Faculty of Agriculture, Shinshu University, Kami-Ina, Nagano, 399-4598, Japan
| | - Yu Harishima
- Department of Bioengineering, University of California, Berkeley, CA, 94704, USA
| | - Musashi Horikawa
- Department of Biological Sciences, Graduate School of Science, University of Tokyo, Tokyo, 113-0033, Japan
| | - Mao Watanabe
- Department of Agricultural and Life Sciences, Faculty of Agriculture, Shinshu University, Kami-Ina, Nagano, 399-4598, Japan
| | - Hiroki Oota
- Department of Biological Sciences, Graduate School of Science, University of Tokyo, Tokyo, 113-0033, Japan
| | - Masato Nakagawa
- Center for iPS Cell Research and Application (CiRA), Kyoto University, Kyoto, 606-8507, Japan
| | - Shunsuke Suzuki
- Department of Agricultural and Life Sciences, Faculty of Agriculture, Shinshu University, Kami-Ina, Nagano, 399-4598, Japan
| | - Wolfgang Enard
- Anthropology and Human Genomics, Faculty of Biology, Ludwig-Maximilians-Universität München, München, Germany
| |
Collapse
|
3
|
Waneka G, Pate B, Monroe JG, Sloan DB. Investigating low frequency somatic mutations in Arabidopsis with Duplex Sequencing. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.31.578196. [PMID: 38352550 PMCID: PMC10862904 DOI: 10.1101/2024.01.31.578196] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/22/2024]
Abstract
Mutations are the source of novel genetic diversity but can also lead to disease and maladaptation. The conventional view is that mutations occur randomly with respect to their environment-specific fitness consequences. However, intragenomic mutation rates can vary dramatically due to transcription coupled repair and based on local epigenomic modifications, which are non-uniformly distributed across genomes. One sequence feature associated with decreased mutation is higher expression level, which can vary depending on environmental cues. To understand whether the association between expression level and mutation rate creates a systematic relationship with environment-specific fitness effects, we perturbed expression through a heat treatment in Arabidopsis thaliana. We quantified gene expression to identify differentially expressed genes, which we then targeted for mutation detection using Duplex Sequencing. This approach provided a highly accurate measurement of the frequency of rare somatic mutations in vegetative plant tissues, which has been a recent source of uncertainty in plant mutation research. We included mutant lines lacking mismatch repair (MMR) and base excision repair (BER) capabilities to understand how repair mechanisms may drive biased mutation accumulation. We found wild type (WT) and BER mutant mutation frequencies to be very low (mean variant frequency 1.8×10-8 and 2.6×10-8, respectively), while MMR mutant frequencies were significantly elevated (1.13×10-6). These results show that somatic variant frequencies are extremely low in WT plants, indicating that larger datasets will be needed to address the fundamental evolutionary question as to whether environmental change leads to gene-specific changes in mutation rate.
Collapse
Affiliation(s)
- Gus Waneka
- Department of Biology, Colorado State University, Fort Collins, Colorado, USA
| | - Braden Pate
- Department of Biology, Colorado State University, Fort Collins, Colorado, USA
| | - J Grey Monroe
- Department of Plant Sciences, University of California, Davis, Davis, CA USA
| | - Daniel B Sloan
- Department of Biology, Colorado State University, Fort Collins, Colorado, USA
| |
Collapse
|
4
|
Lin Y, Darolti I, van der Bijl W, Morris J, Mank JE. Extensive variation in germline de novo mutations in Poecilia reticulata. Genome Res 2023; 33:1317-1324. [PMID: 37442578 PMCID: PMC10547258 DOI: 10.1101/gr.277936.123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2023] [Accepted: 07/07/2023] [Indexed: 07/15/2023]
Abstract
The rate of germline mutation is fundamental to evolutionary processes, as it generates the variation upon which selection acts. The guppy, Poecilia reticulata, is a model of rapid adaptation, however the relative contribution of standing genetic variation versus de novo mutation (DNM) to evolution in this species remains unclear. Here, we use pedigree-based approaches to quantify and characterize germline DNMs in three large guppy families. Our results suggest germline mutation rate in the guppy varies substantially across individuals and families. Most DNMs are shared across multiple siblings, suggesting they arose during early embryonic development. DNMs are randomly distributed throughout the genome, and male-biased mutation rate is low, as would be expected from the short guppy generation time. Overall, our study shows remarkable variation in germline mutation rate and provides insights into rapid evolution of guppies.
Collapse
Affiliation(s)
- Yuying Lin
- Department of Zoology and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada;
| | - Iulia Darolti
- Department of Ecology and Evolution, University of Lausanne, CH-1015 Lausanne, Switzerland
| | - Wouter van der Bijl
- Department of Zoology and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada
| | - Jake Morris
- School of Biological Science, University of Bristol, Bristol BS8 1TQ, United Kingdom
| | - Judith E Mank
- Department of Zoology and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada
| |
Collapse
|
5
|
Abstract
Mutations in the germline are an essential component of evolution. Current DNA sequencing systems allow quantification of these mutations through full genome sequencing of father/mother/child trios. A recent paper reports such a large-scale study spanning 68 vertebrate species, and establishes important features of the mutation process.
Collapse
Affiliation(s)
- Bertrand Jordan
- Biologiste, généticien et immunologiste, Président d'Aprogène (Association pour la promotion de la Génomique), 13007 Marseille, France
| |
Collapse
|
6
|
Bergeron LA, Besenbacher S, Zheng J, Li P, Bertelsen MF, Quintard B, Hoffman JI, Li Z, St Leger J, Shao C, Stiller J, Gilbert MTP, Schierup MH, Zhang G. Evolution of the germline mutation rate across vertebrates. Nature 2023; 615:285-291. [PMID: 36859541 PMCID: PMC9995274 DOI: 10.1038/s41586-023-05752-y] [Citation(s) in RCA: 63] [Impact Index Per Article: 63.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Accepted: 01/23/2023] [Indexed: 03/03/2023]
Abstract
The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself1. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies2. Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis3. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.
Collapse
Affiliation(s)
- Lucie A Bergeron
- Villum Centre for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Søren Besenbacher
- Department of Molecular Medicine, Aarhus University, Aarhus, Denmark
| | - Jiao Zheng
- BGI-Shenzhen, Shenzhen, China
- BGI Education Center, University of Chinese Academy of Sciences, Shenzhen, China
| | | | | | | | - Joseph I Hoffman
- Department of Animal Behaviour, Bielefeld University, Bielefeld, Germany
- British Antarctic Survey, High Cross, Cambridge, UK
| | - Zhipeng Li
- College of Animal Science and Technology, Jilin Agricultural University, Changchun, China
| | - Judy St Leger
- Department of Biomedical Sciences, Cornell University, Ithaca, NY, USA
| | - Changwei Shao
- Key Lab of Sustainable Development of Marine Fisheries, Ministry of Agriculture and Rural Affairs, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao, China
| | - Josefin Stiller
- Villum Centre for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - M Thomas P Gilbert
- Center for Evolutionary Hologenomics, The GLOBE Institute, University of Copenhagen, Copenhagen, Denmark
- University Museum, NTNU, Trondheim, Norway
| | | | - Guojie Zhang
- Villum Centre for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
- Centre for Evolutionary & Organismal Biology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China.
- Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, China.
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China.
| |
Collapse
|
7
|
Ren X, Liu Y, Zhao Y, Li B, Bai D, Bou G, Zhang X, Du M, Wang X, Bou T, Shen Y, Dugarjaviin M. Analysis of the Whole-Genome Sequences from an Equus Parent-Offspring Trio Provides Insight into the Genomic Incompatibilities in the Hybrid Mule. Genes (Basel) 2022; 13:genes13122188. [PMID: 36553455 PMCID: PMC9778318 DOI: 10.3390/genes13122188] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Revised: 11/07/2022] [Accepted: 11/21/2022] [Indexed: 11/25/2022] Open
Abstract
Interspecific hybridization often shows negative effects on hybrids. However, only a few multicellular species, limited to a handful of plants and animals, have shown partial genetic mechanisms by which hybridization leads to low fitness in hybrids. Here, to explore the outcome of combining the two genomes of a horse and donkey, we analyzed the whole-genome sequences from an Equus parent-offspring trio using Illumina platforms. We generated 41.39× and 46.21× coverage sequences for the horse and mule, respectively. For the donkey, a 40.38× coverage sequence was generated and stored in our laboratory. Approximately 24.86 million alleles were discovered that varied from the reference genome. Single nucleotide polymorphisms were used as polymorphic markers for assigning alleles to their parental genomic inheritance. We identified 25,703 Mendelian inheritance error single nucleotide polymorphisms in the mule genome that were not inherited from the parents through Mendelian inheritance. A total of 555 de novo single nucleotide polymorphisms were also identified. The rate of de novo single nucleotide polymorphisms was 2.21 × 10-7 in the mule from the Equus parent-offspring trio. This rate is obviously higher than the natural mutation rate for Equus, which is also consistent with the previous hypothesis that interracial crosses may have a high mutation rate. The genes associated with these single nucleotide polymorphisms are mainly involved in immune processes, DNA repair, and cancer processes. The results of the analysis of three genomes from an Equus parent-offspring trio improved our knowledge of the consequences of the integration of parental genomes in mules.
Collapse
|
8
|
de Manuel M, Wu FL, Przeworski M. A paternal bias in germline mutation is widespread in amniotes and can arise independently of cell divisions. eLife 2022; 11:80008. [PMID: 35916372 PMCID: PMC9439683 DOI: 10.7554/elife.80008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2022] [Accepted: 08/01/2022] [Indexed: 11/13/2022] Open
Abstract
In humans and other mammals, germline mutations are more likely to arise in fathers than in mothers. Although this sex bias has long been attributed to DNA replication errors in spermatogenesis, recent evidence from humans points to the importance of mutagenic processes that do not depend on cell division, calling into question our understanding of this basic phenomenon. Here, we infer the ratio of paternal-to-maternal mutations, α, in 42 species of amniotes, from putatively neutral substitution rates of sex chromosomes and autosomes. Despite marked differences in gametogenesis, physiologies and environments across species, fathers consistently contribute more mutations than mothers in all the species examined, including mammals, birds, and reptiles. In mammals, α is as high as 4 and correlates with generation times; in birds and snakes, α appears more stable around 2. These observations are consistent with a simple model, in which mutations accrue at equal rates in both sexes during early development and at a higher rate in the male germline after sexual differentiation, with a conserved paternal-to-maternal ratio across species. Thus, α may reflect the relative contributions of two or more developmental phases to total germline mutations, and is expected to depend on generation time even if mutations do not track cell divisions.
Collapse
Affiliation(s)
- Marc de Manuel
- Department of Biological Sciences, Columbia University, New York, United States
| | - Felix L Wu
- Department of Biological Sciences, Columbia University, New York, United States
| | - Molly Przeworski
- Department of Systems Biology, Columbia University, New York, United States
| |
Collapse
|
9
|
Bergeron LA, Besenbacher S, Turner T, Versoza CJ, Wang RJ, Price AL, Armstrong E, Riera M, Carlson J, Chen HY, Hahn MW, Harris K, Kleppe AS, López-Nandam EH, Moorjani P, Pfeifer SP, Tiley GP, Yoder AD, Zhang G, Schierup MH. The mutationathon highlights the importance of reaching standardization in estimates of pedigree-based germline mutation rates. eLife 2022; 11:73577. [PMID: 35018888 PMCID: PMC8830884 DOI: 10.7554/elife.73577] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Accepted: 01/11/2022] [Indexed: 11/13/2022] Open
Abstract
In the past decade, several studies have estimated the human per-generation germline mutation rate using large pedigrees. More recently, estimates for various nonhuman species have been published. However, methodological differences among studies in detecting germline mutations and estimating mutation rates make direct comparisons difficult. Here, we describe the many different steps involved in estimating pedigree-based mutation rates, including sampling, sequencing, mapping, variant calling, filtering, and appropriately accounting for false-positive and false-negative rates. For each step, we review the different methods and parameter choices that have been used in the recent literature. Additionally, we present the results from a ‘Mutationathon,’ a competition organized among five research labs to compare germline mutation rate estimates for a single pedigree of rhesus macaques. We report almost a twofold variation in the final estimated rate among groups using different post-alignment processing, calling, and filtering criteria, and provide details into the sources of variation across studies. Though the difference among estimates is not statistically significant, this discrepancy emphasizes the need for standardized methods in mutation rate estimations and the difficulty in comparing rates from different studies. Finally, this work aims to provide guidelines for computational and statistical benchmarks for future studies interested in identifying germline mutations from pedigrees.
Collapse
Affiliation(s)
- Lucie A Bergeron
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Søren Besenbacher
- Department of Molecular Medicine (MOMA), Aarhus University, Aarhus N, Denmark
| | - Tychele Turner
- Department of Genetics, Washington University in St. Louis, Saint Louis, United States
| | - Cyril J Versoza
- Center for Evolution and Medicine, Arizona State University, Tempe, United States
| | - Richard J Wang
- Department of Biology, Indiana University, Bloomington, United States
| | - Alivia Lee Price
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Ellie Armstrong
- Department of Biology, Stanford University, Stanford, United States
| | - Meritxell Riera
- Bioinformatics Research Centre, Aarhus University, Aarhus, Denmark
| | - Jedidiah Carlson
- Department of Genome Sciences, University of Washington, Seattle, United States
| | - Hwei-Yen Chen
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Matthew W Hahn
- Department of Biology, Indiana University, Bloomington, United States
| | - Kelley Harris
- Department of Genome Sciences, University of Washington, Seattle, United States
| | | | | | - Priya Moorjani
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, United States
| | - Susanne P Pfeifer
- School of Life Sciences, Arizona State University, Tempe, United States
| | - George P Tiley
- Department of Biology, Duke University, Durham, United States
| | - Anne D Yoder
- Department of Biology, Duke University, Durham, United States
| | - Guojie Zhang
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | | |
Collapse
|
10
|
Rashid I, Campos M, Collier T, Crepeau M, Weakley A, Gripkey H, Lee Y, Schmidt H, Lanzaro GC. Spontaneous mutation rate estimates for the principal malaria vectors Anopheles coluzzii and Anopheles stephensi. Sci Rep 2022; 12:226. [PMID: 34996998 PMCID: PMC8742016 DOI: 10.1038/s41598-021-03943-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2021] [Accepted: 12/07/2021] [Indexed: 11/17/2022] Open
Abstract
Using high-depth whole genome sequencing of F0 mating pairs and multiple individual F1 offspring, we estimated the nuclear mutation rate per generation in the malaria vectors Anopheles coluzzii and Anopheles stephensi by detecting de novo genetic mutations. A purpose-built computer program was employed to filter actual mutations from a deep background of superficially similar artifacts resulting from read misalignment. Performance of filtering parameters was determined using software-simulated mutations, and the resulting estimate of false negative rate was used to correct final mutation rate estimates. Spontaneous mutation rates by base substitution were estimated at 1.00 × 10−9 (95% confidence interval, 2.06 × 10−10—2.91 × 10−9) and 1.36 × 10−9 (95% confidence interval, 4.42 × 10−10—3.18 × 10−9) per site per generation in A. coluzzii and A. stephensi respectively. Although similar studies have been performed on other insect species including dipterans, this is the first study to empirically measure mutation rates in the important genus Anopheles, and thus provides an estimate of µ that will be of utility for comparative evolutionary genomics, as well as for population genetic analysis of malaria vector mosquito species.
Collapse
Affiliation(s)
- Iliyas Rashid
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, UC Davis, 1089 Veterinary Medicine Dr, 4225 VM3B, Davis, CA, 95616, USA.,Section of Cell and Developmental Biology, University of California, San Diego, La Jolla, CA, USA.,Tata Institute for Genetics and Society, Center at inStem, Bangalore, Karnataka, 560065, India
| | - Melina Campos
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, UC Davis, 1089 Veterinary Medicine Dr, 4225 VM3B, Davis, CA, 95616, USA
| | - Travis Collier
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, UC Davis, 1089 Veterinary Medicine Dr, 4225 VM3B, Davis, CA, 95616, USA
| | - Marc Crepeau
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, UC Davis, 1089 Veterinary Medicine Dr, 4225 VM3B, Davis, CA, 95616, USA
| | - Allison Weakley
- Department of ChEM-H Operations, Stanford University, 450 Serra Mall, Stanford, CA, 94305, USA
| | - Hans Gripkey
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, UC Davis, 1089 Veterinary Medicine Dr, 4225 VM3B, Davis, CA, 95616, USA
| | - Yoosook Lee
- Florida Medical Entomology Laboratory, University of Florida, 200 9th St SE, Vero Beach, FL, 32962, USA
| | - Hanno Schmidt
- Anthropology, Institute of Organismic and Molecular Evolution (iomE), Johannes Gutenberg University of Mainz, Saarstraße 21, 55122, Mainz, Germany
| | - Gregory C Lanzaro
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, UC Davis, 1089 Veterinary Medicine Dr, 4225 VM3B, Davis, CA, 95616, USA.
| |
Collapse
|
11
|
Warmuth VM, Burgess MD, Laaksonen T, Manica A, Mägi M, Nord A, Primmer CR, Sætre GP, Winkel W, Ellegren H. Major population splits coincide with episodes of rapid climate change in a forest-dependent bird. Proc Biol Sci 2021; 288:20211066. [PMID: 34727712 PMCID: PMC8564624 DOI: 10.1098/rspb.2021.1066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Accepted: 10/07/2021] [Indexed: 11/17/2022] Open
Abstract
Climate change influences population demography by altering patterns of gene flow and reproductive isolation. Direct mutation rates offer the possibility for accurate dating on the within-species level but are currently only available for a handful of vertebrate species. Here, we use the first directly estimated mutation rate in birds to study the evolutionary history of pied flycatchers (Ficedula hypoleuca). Using a combination of demographic inference and species distribution modelling, we show that all major population splits in this forest-dependent system occurred during periods of increased climate instability and rapid global temperature change. We show that the divergent Spanish subspecies originated during the Eemian-Weichselian transition 115-104 thousand years ago (kya), and not during the last glacial maximum (26.5-19 kya), as previously suggested. The magnitude and rates of climate change during the glacial-interglacial transitions that preceded population splits in pied flycatchers were similar to, or exceeded, those predicted to occur in the course of the current, human-induced climate crisis. As such, our results provide a timely reminder of the strong impact that episodes of climate instability and rapid temperature changes can have on species' evolutionary trajectories, with important implications for the natural world in the Anthropocene.
Collapse
Affiliation(s)
- Vera M. Warmuth
- Department of Evolutionary Biology, Biozentrum Martinsried, Ludwig-Maximilians Universität München, Planegg-Martinsried, Germany
- Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden
| | - Malcolm D. Burgess
- Centre for Animal Behaviour, University of Exeter, Exeter, UK
- RSPB Centre for Conservation Science, Sandy, UK
| | - Toni Laaksonen
- Department of Biology, University of Turku, Turku, Finland
| | - Andrea Manica
- Department of Zoology, University of Cambridge, Cambridge, UK
| | - Marko Mägi
- Institute of Ecology and Earth Sciences, University of Tartu, Tartu, Estonia
| | - Andreas Nord
- Department of Biology, Section for Evolutionary Ecology, Lund University, Lund, Sweden
| | - Craig R. Primmer
- Organismal and Evolutionary Biology Research Program, University of Helsinki, Finland
- Institute of Biotechnology, Helsinki Institute of Life Sciences (HiLIFE), University of Helsinki, Finland
| | - Glenn-Peter Sætre
- Centre for Ecological and Evolutionary Synthesis, University of Oslo, Oslo, Norway
| | - Wolfgang Winkel
- Institute of Avian Research, ‘Vogelwarte Helgoland’, Wilhelmshaven, Germany
| | - Hans Ellegren
- Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden
| |
Collapse
|
12
|
Bergeron LA, Besenbacher S, Schierup MH, Zhang G. Studying mutation rate evolution in primates-a need for systematic comparison of computational pipelines. Gigascience 2021; 10:giab072. [PMID: 34673928 PMCID: PMC8529962 DOI: 10.1093/gigascience/giab072] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2021] [Revised: 10/04/2021] [Accepted: 10/04/2021] [Indexed: 11/14/2022] Open
Abstract
The lack of consensus methods to estimate germline mutation rates from pedigrees has led to substantial differences in computational pipelines in the published literature. Here, we answer Susanne Pfeifer's opinion piece discussing the pipeline choices of our recent article estimating the germline mutation rate of rhesus macaques (Macaca mulatta). We acknowledge the differences between the method that we applied and the one preferred by Pfeifer. Yet, we advocate for full transparency and justification of choices as long as rigorous comparison of pipelines remains absent because it is the only way to conclude on best practices for the field.
Collapse
Affiliation(s)
- Lucie A Bergeron
- Section for Ecology and Evolution, Department of Biology, University of
Copenhagen, 2100, Copenhagen, Denmark
| | - Søren Besenbacher
- Department of Molecular Medicine, Aarhus University, 8200,
Aarhus, Denmark
| | - Mikkel H Schierup
- Bioinformatics Research Centre, Aarhus University, 8000,
Aarhus, Denmark
| | - Guojie Zhang
- Section for Ecology and Evolution, Department of Biology, University of
Copenhagen, 2100, Copenhagen, Denmark
- BGI-Shenzhen, Shenzhen, Guangdong 518083, China
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute
of Zoology, Chinese Academy of Sciences, Kunming 650223, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of
Sciences, Kunming 650223, China
| |
Collapse
|
13
|
Pfeifer SP. Studying mutation rate evolution in primates-the effects of computational pipelines and parameter choices. Gigascience 2021; 10:6407243. [PMID: 34673929 PMCID: PMC8529961 DOI: 10.1093/gigascience/giab069] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2021] [Accepted: 09/13/2021] [Indexed: 11/12/2022] Open
Abstract
This commentary investigates the important role of computational pipeline and parameter choices in performing mutation rate estimation, using the recent article published in this journal by Bergeron et al. entitled "The germline mutational process in rhesus macaque and its implications for phylogenetic dating" as an illustrative example.
Collapse
Affiliation(s)
- Susanne P Pfeifer
- Correspondence address. Susanne P. Pfeifer, School of Life Sciences, Arizona
State University, Tempe, AZ 85281, USA. E-mail:
| |
Collapse
|
14
|
Campbell CR, Tiley GP, Poelstra JW, Hunnicutt KE, Larsen PA, Lee HJ, Thorne JL, Dos Reis M, Yoder AD. Pedigree-based and phylogenetic methods support surprising patterns of mutation rate and spectrum in the gray mouse lemur. Heredity (Edinb) 2021; 127:233-244. [PMID: 34272504 PMCID: PMC8322134 DOI: 10.1038/s41437-021-00446-5] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2020] [Revised: 05/25/2021] [Accepted: 05/26/2021] [Indexed: 02/06/2023] Open
Abstract
Mutations are the raw material on which evolution acts, and knowledge of their frequency and genomic distribution is crucial for understanding how evolution operates at both long and short timescales. At present, the rate and spectrum of de novo mutations have been directly characterized in relatively few lineages. Our study provides the first direct mutation-rate estimate for a strepsirrhine (i.e., the lemurs and lorises), which comprises nearly half of the primate clade. Using high-coverage linked-read sequencing for a focal quartet of gray mouse lemurs (Microcebus murinus), we estimated the mutation rate to be among the highest calculated for a mammal at 1.52 × 10-8 (95% credible interval: 1.28 × 10-8-1.78 × 10-8) mutations/site/generation. Further, we found an unexpectedly low count of paternal mutations, and only a modest overrepresentation of mutations at CpG sites. Despite the surprising nature of these results, we found both the rate and spectrum to be robust to the manipulation of a wide range of computational filtering criteria. We also sequenced a technical replicate to estimate a false-negative and false-positive rate for our data and show that any point estimate of a de novo mutation rate should be considered with a large degree of uncertainty. For validation, we conducted an independent analysis of context-dependent substitution types for gray mouse lemur and five additional primate species for which de novo mutation rates have also been estimated. These comparisons revealed general consistency of the mutation spectrum between the pedigree-based and the substitution-rate analyses for all species compared.
Collapse
Affiliation(s)
- C Ryan Campbell
- Department of Biology, Duke University, Durham, NC, USA
- Department of Evolutionary Anthropology, Duke University, Durham, NC, USA
| | | | | | - Kelsie E Hunnicutt
- Department of Biology, Duke University, Durham, NC, USA
- Department of Biological Sciences, University of Denver, Denver, CO, USA
| | - Peter A Larsen
- Department of Biology, Duke University, Durham, NC, USA
- Department of Veterinary and Biomedical Sciences, University of Minnesota, St. Paul, MN, USA
| | - Hui-Jie Lee
- Department of Biostatistics and Bioinformatics, Duke University, Durham, NC, USA
| | - Jeffrey L Thorne
- Bioinformatics Research Center, North Carolina State University, Raleigh, NC, USA
| | - Mario Dos Reis
- School of Biological and Chemical Sciences, Queen Mary University of London, London, UK
| | - Anne D Yoder
- Department of Biology, Duke University, Durham, NC, USA.
| |
Collapse
|
15
|
Yoder AD, Tiley GP. The challenge and promise of estimating the de novo mutation rate from whole-genome comparisons among closely related individuals. Mol Ecol 2021; 30:6087-6100. [PMID: 34062029 DOI: 10.1111/mec.16007] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 04/22/2021] [Accepted: 05/26/2021] [Indexed: 12/20/2022]
Abstract
Germline mutations are the raw material for natural selection, driving species evolution and the generation of earth's biodiversity. Without this driver of genetic diversity, life on earth would stagnate. Yet, it is a double-edged sword. An excess of mutations can have devastating effects on fitness and population viability. It is therefore one of the great challenges of molecular ecology to determine the rate and mechanisms by which these mutations accrue across the tree of life. Advances in high-throughput sequencing technologies are providing new opportunities for characterizing the rates and mutational spectra within species and populations thus informing essential evolutionary parameters such as the timing of speciation events, the intricacies of historical demography, and the degree to which lineages are subject to the burdens of mutational load. Here, we will focus on both the challenge and promise of whole-genome comparisons among parents and their offspring from known pedigrees for the detection of germline mutations as they arise in a single generation. The potential of these studies is high, but the field is still in its infancy and much uncertainty remains. Namely, the technical challenges are daunting given that pedigree-based genome comparisons are essentially searching for needles in a haystack given the very low signal to noise ratio. Despite the challenges, we predict that rapidly developing methods for whole-genome comparisons hold great promise for integrating empirically derived estimates of de novo mutation rates and mutation spectra across many molecular ecological applications.
Collapse
Affiliation(s)
- Anne D Yoder
- Department of Biology, Duke University, Durham, NC, USA
| | | |
Collapse
|
16
|
Stability across the Whole Nuclear Genome in the Presence and Absence of DNA Mismatch Repair. Cells 2021; 10:cells10051224. [PMID: 34067668 PMCID: PMC8156620 DOI: 10.3390/cells10051224] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Revised: 05/13/2021] [Accepted: 05/14/2021] [Indexed: 01/06/2023] Open
Abstract
We describe the contribution of DNA mismatch repair (MMR) to the stability of the eukaryotic nuclear genome as determined by whole-genome sequencing. To date, wild-type nuclear genome mutation rates are known for over 40 eukaryotic species, while measurements in mismatch repair-defective organisms are fewer in number and are concentrated on Saccharomyces cerevisiae and human tumors. Well-studied organisms include Drosophila melanogaster and Mus musculus, while less genetically tractable species include great apes and long-lived trees. A variety of techniques have been developed to gather mutation rates, either per generation or per cell division. Generational rates are described through whole-organism mutation accumulation experiments and through offspring–parent sequencing, or they have been identified by descent. Rates per somatic cell division have been estimated from cell line mutation accumulation experiments, from systemic variant allele frequencies, and from widely spaced samples with known cell divisions per unit of tissue growth. The latter methods are also used to estimate generational mutation rates for large organisms that lack dedicated germlines, such as trees and hyphal fungi. Mechanistic studies involving genetic manipulation of MMR genes prior to mutation rate determination are thus far confined to yeast, Arabidopsis thaliana, Caenorhabditis elegans, and one chicken cell line. A great deal of work in wild-type organisms has begun to establish a sound baseline, but far more work is needed to uncover the variety of MMR across eukaryotes. Nonetheless, the few MMR studies reported to date indicate that MMR contributes 100-fold or more to genome stability, and they have uncovered insights that would have been impossible to obtain using reporter gene assays.
Collapse
|
17
|
Bergeron LA, Besenbacher S, Bakker J, Zheng J, Li P, Pacheco G, Sinding MHS, Kamilari M, Gilbert MTP, Schierup MH, Zhang G. The germline mutational process in rhesus macaque and its implications for phylogenetic dating. Gigascience 2021; 10:giab029. [PMID: 33954793 PMCID: PMC8099771 DOI: 10.1093/gigascience/giab029] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2020] [Revised: 01/05/2021] [Accepted: 03/29/2021] [Indexed: 01/13/2023] Open
Abstract
BACKGROUND Understanding the rate and pattern of germline mutations is of fundamental importance for understanding evolutionary processes. RESULTS Here we analyzed 19 parent-offspring trios of rhesus macaques (Macaca mulatta) at high sequencing coverage of ∼76× per individual and estimated a mean rate of 0.77 × 10-8de novo mutations per site per generation (95% CI: 0.69 × 10-8 to 0.85 × 10-8). By phasing 50% of the mutations to parental origins, we found that the mutation rate is positively correlated with the paternal age. The paternal lineage contributed a mean of 81% of the de novo mutations, with a trend of an increasing male contribution for older fathers. Approximately 3.5% of de novo mutations were shared between siblings, with no parental bias, suggesting that they arose from early development (postzygotic) stages. Finally, the divergence times between closely related primates calculated on the basis of the yearly mutation rate of rhesus macaque generally reconcile with divergence estimated with molecular clock methods, except for the Cercopithecoidea/Hominoidea molecular divergence dated at 58 Mya using our new estimate of the yearly mutation rate. CONCLUSIONS When compared to the traditional molecular clock methods, new estimated rates from pedigree samples can provide insights into the evolution of well-studied groups such as primates.
Collapse
Affiliation(s)
- Lucie A Bergeron
- Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Universitetsparken 15, 2100 Copenhagen Ø, Denmark
| | - Søren Besenbacher
- Department of Molecular Medicine, Aarhus University, Brendstrupgårdsvej 21A, 8200 Aarhus N, Denmark
| | - Jaco Bakker
- Animal Science Department, Biomedical Primate Research Centre, Lange Kleiweg 161, 2288 GJ Rijswijk, Netherlands
| | - Jiao Zheng
- BGI-Shenzhen, Shenzhen 518083, Guangdong, China
- BGI Education Center, University of Chinese Academy of Sciences, Shenzhen 518083, Guangdong, China
| | - Panyi Li
- BGI-Shenzhen, Shenzhen 518083, Guangdong, China
| | - George Pacheco
- Section for Evolutionary Genomics, The GLOBE Institute, University of Copenhagen, Oester Voldgade 5-7, 1350 Copenhagen K, Denmark
| | - Mikkel-Holger S Sinding
- Department of genetics, Trinity College Dublin, 2 college green, D02 VF25, Dublin, Ireland
- Greenland Institute of Natural Resources, Kivioq 2, 3900 Nuuk, Greenland
| | - Maria Kamilari
- Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Universitetsparken 15, 2100 Copenhagen Ø, Denmark
| | - M Thomas P Gilbert
- Section for Evolutionary Genomics, The GLOBE Institute, University of Copenhagen, Oester Voldgade 5-7, 1350 Copenhagen K, Denmark
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology (NTNU), NO-7491 Trondheim, Norway
| | - Mikkel H Schierup
- Bioinformatics Research Centre, Aarhus University, C.F.Møllers Allé 8, 8000, Aarhus C, Denmark
| | - Guojie Zhang
- Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Universitetsparken 15, 2100 Copenhagen Ø, Denmark
- BGI-Shenzhen, Shenzhen 518083, Guangdong, China
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming 650223, China
| |
Collapse
|
18
|
Ho EKH, Macrae F, Latta LC, McIlroy P, Ebert D, Fields PD, Benner MJ, Schaack S. High and Highly Variable Spontaneous Mutation Rates in Daphnia. Mol Biol Evol 2021; 37:3258-3266. [PMID: 32520985 PMCID: PMC7820357 DOI: 10.1093/molbev/msaa142] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The rate and spectrum of spontaneous mutations are critical parameters in basic and applied biology because they dictate the pace and character of genetic variation introduced into populations, which is a prerequisite for evolution. We use a mutation–accumulation approach to estimate mutation parameters from whole-genome sequence data from multiple genotypes from multiple populations of Daphnia magna, an ecological and evolutionary model system. We report extremely high base substitution mutation rates (µ-n,bs = 8.96 × 10−9/bp/generation [95% CI: 6.66–11.97 × 10−9/bp/generation] in the nuclear genome and µ-m,bs = 8.7 × 10−7/bp/generation [95% CI: 4.40–15.12 × 10−7/bp/generation] in the mtDNA), the highest of any eukaryote examined using this approach. Levels of intraspecific variation based on the range of estimates from the nine genotypes collected from three populations (Finland, Germany, and Israel) span 1 and 3 orders of magnitude, respectively, resulting in up to a ∼300-fold difference in rates among genomic partitions within the same lineage. In contrast, mutation spectra exhibit very consistent patterns across genotypes and populations, suggesting the mechanisms underlying the mutational process may be similar, even when the rates at which they occur differ. We discuss the implications of high levels of intraspecific variation in rates, the importance of estimating gene conversion rates using a mutation–accumulation approach, and the interacting factors influencing the evolution of mutation parameters. Our findings deepen our knowledge about mutation and provide both challenges to and support for current theories aimed at explaining the evolution of the mutation rate, as a trait, across taxa.
Collapse
Affiliation(s)
- Eddie K H Ho
- Department of Biology, Reed College, Portland, OR
| | | | - Leigh C Latta
- Department of Biology, Reed College, Portland, OR.,Division of Natural Sciences and Mathematics, Lewis-Clark State College, Lewiston, ID
| | | | - Dieter Ebert
- Department of Environmental Sciences, Zoology, University of Basel, Basel, Switzerland
| | - Peter D Fields
- Department of Environmental Sciences, Zoology, University of Basel, Basel, Switzerland
| | | | | |
Collapse
|
19
|
Finke K, Kourakos M, Brown G, Dang HT, Tan SJS, Simons YB, Ramdas S, Schäffer AA, Kember RL, Bućan M, Mathieson S. Ancestral haplotype reconstruction in endogamous populations using identity-by-descent. PLoS Comput Biol 2021; 17:e1008638. [PMID: 33635861 PMCID: PMC7946327 DOI: 10.1371/journal.pcbi.1008638] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2020] [Revised: 03/10/2021] [Accepted: 12/15/2020] [Indexed: 12/24/2022] Open
Abstract
In this work we develop a novel algorithm for reconstructing the genomes of ancestral individuals, given genotype or sequence data from contemporary individuals and an extended pedigree of family relationships. A pedigree with complete genomes for every individual enables the study of allele frequency dynamics and haplotype diversity across generations, including deviations from neutrality such as transmission distortion. When studying heritable diseases, ancestral haplotypes can be used to augment genome-wide association studies and track disease inheritance patterns. The building blocks of our reconstruction algorithm are segments of Identity-By-Descent (IBD) shared between two or more genotyped individuals. The method alternates between identifying a source for each IBD segment and assembling IBD segments placed within each ancestral individual. Unlike previous approaches, our method is able to accommodate complex pedigree structures with hundreds of individuals genotyped at millions of SNPs. We apply our method to an Old Order Amish pedigree from Lancaster, Pennsylvania, whose founders came to North America from Europe during the early 18th century. The pedigree includes 1338 individuals from the past 12 generations, 394 with genotype data. The motivation for reconstruction is to understand the genetic basis of diseases segregating in the family through tracking haplotype transmission over time. Using our algorithm thread, we are able to reconstruct an average of 224 ancestral individuals per chromosome. For these ancestral individuals, on average we reconstruct 79% of their haplotypes. We also identify a region on chromosome 16 that is difficult to reconstruct—we find that this region harbors a short Amish-specific copy number variation and the gene HYDIN. thread was developed for endogamous populations, but can be applied to any extensive pedigree with the recent generations genotyped. We anticipate that this type of practical ancestral reconstruction will become more common and necessary to understand rare and complex heritable diseases in extended families. When analyzing complex heritable traits, genomic data from many generations of an extended family increases the amount of information available for statistical inference. However, typically only genomic data from the recent generations of a pedigree are available, as ancestral individuals are deceased. In this work we present an algorithm, called thread, for reconstructing the genomes of ancestral individuals, given a complex pedigree and genomic data from the recent generations. Previous approaches have not been able to accommodate large datasets (both in terms of sites and individuals), made simplifying assumptions about pedigree structure, or did not tie reconstructed sequences back to specific individuals. We apply thread to a complex Old Order Amish pedigree of 1338 individuals, 394 with genotype data.
Collapse
Affiliation(s)
- Kelly Finke
- Department of Computer Science, Swarthmore College, Swarthmore, Pennsylvania, United States of America
- Department of Biology, Swarthmore College, Swarthmore, Pennsylvania, United States of America
| | - Michael Kourakos
- Department of Computer Science, Swarthmore College, Swarthmore, Pennsylvania, United States of America
| | - Gabriela Brown
- Department of Computer Science, Swarthmore College, Swarthmore, Pennsylvania, United States of America
| | - Huyen Trang Dang
- Department of Computer Science, Bryn Mawr College, Bryn Mawr, Pennsylvania, United States of America
| | - Shi Jie Samuel Tan
- Department of Computer Science, Haverford College, Haverford, Pennsylvania, United States of America
| | - Yuval B. Simons
- Department of Genetics, Stanford University, Stanford, California, United States of America
| | - Shweta Ramdas
- Department of Genetics, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Alejandro A. Schäffer
- Cancer Data Science Laboratory, National Cancer Institute, National Institutes of Health, Bethesda, Maryland, United States of America
| | - Rachel L. Kember
- Department of Psychiatry, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Maja Bućan
- Department of Genetics, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Sara Mathieson
- Department of Computer Science, Haverford College, Haverford, Pennsylvania, United States of America
- * E-mail:
| |
Collapse
|
20
|
Changes in life history and population size can explain the relative neutral diversity levels on X and autosomes in extant human populations. Proc Natl Acad Sci U S A 2020; 117:20063-20069. [PMID: 32747577 DOI: 10.1073/pnas.1915664117] [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] [Indexed: 01/17/2023] Open
Abstract
In human populations, the relative levels of neutral diversity on the X and autosomes differ markedly from each other and from the naïve theoretical expectation of 3/4. Here we propose an explanation for these differences based on new theory about the effects of sex-specific life history and given pedigree-based estimates of the dependence of human mutation rates on sex and age. We demonstrate that life history effects, particularly longer generation times in males than in females, are expected to have had multiple effects on human X-to-autosome (X:A) diversity ratios, as a result of male-biased mutation rates, the equilibrium X:A ratio of effective population sizes, and the differential responses to changes in population size. We also show that the standard approach of using divergence between species to correct for male mutation bias results in biased estimates of X:A effective population size ratios. We obtain alternative estimates using pedigree-based estimates of the male mutation bias, which reveal that X:A ratios of effective population sizes are considerably greater than previously appreciated. Finally, we find that the joint effects of historical changes in life history and population size can explain the observed X:A diversity ratios in extant human populations. Our results suggest that ancestral human populations were highly polygynous, that non-African populations experienced a substantial reduction in polygyny and/or increase in the male-to-female ratio of generation times around the Out-of-Africa bottleneck, and that current diversity levels were affected by fairly recent changes in sex-specific life history.
Collapse
|
21
|
Wu FL, Strand AI, Cox LA, Ober C, Wall JD, Moorjani P, Przeworski M. A comparison of humans and baboons suggests germline mutation rates do not track cell divisions. PLoS Biol 2020; 18:e3000838. [PMID: 32804933 PMCID: PMC7467331 DOI: 10.1371/journal.pbio.3000838] [Citation(s) in RCA: 40] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2019] [Revised: 09/02/2020] [Accepted: 07/28/2020] [Indexed: 12/19/2022] Open
Abstract
In humans, most germline mutations are inherited from the father. This observation has been widely interpreted as reflecting the replication errors that accrue during spermatogenesis. If so, the male bias in mutation should be substantially lower in a closely related species with similar rates of spermatogonial stem cell divisions but a shorter mean age of reproduction. To test this hypothesis, we resequenced two 3-4 generation nuclear families (totaling 29 individuals) of olive baboons (Papio anubis), who reproduce at approximately 10 years of age on average, and analyzed the data in parallel with three 3-generation human pedigrees (26 individuals). We estimated a mutation rate per generation in baboons of 0.57×10-8 per base pair, approximately half that of humans. Strikingly, however, the degree of male bias in germline mutations is approximately 4:1, similar to that of humans-indeed, a similar male bias is seen across mammals that reproduce months, years, or decades after birth. These results mirror the finding in humans that the male mutation bias is stable with parental ages and cast further doubt on the assumption that germline mutations track cell divisions. Our mutation rate estimates for baboons raise a further puzzle, suggesting a divergence time between apes and Old World monkeys of 65 million years, too old to be consistent with the fossil record; reconciling them now requires not only a slowdown of the mutation rate per generation in humans but also in baboons.
Collapse
Affiliation(s)
- Felix L. Wu
- Department of Systems Biology, Columbia University, New York, New York, United States of America
- Integrated Program in Cellular, Molecular, and Biomedical Studies, Columbia University, New York, New York, United States of America
| | - Alva I. Strand
- Department of Biological Sciences, Columbia University, New York, New York, United States of America
| | - Laura A. Cox
- Center for Precision Medicine, Department of Internal Medicine, Section of Molecular Medicine, Wake Forest School of Medicine, Winston-Salem, North Carolina, United States of America
- Southwest National Primate Research Center, Texas Biomedical Research Institute, San Antonio, Texas, United States of America
| | - Carole Ober
- Department of Human Genetics, The University of Chicago, Chicago, Illinois, United States of America
| | - Jeffrey D. Wall
- Institute for Human Genetics, Department of Epidemiology & Biostatistics, University of California, San Francisco, San Francisco, California, United States of America
| | - Priya Moorjani
- Department of Biological Sciences, Columbia University, New York, New York, United States of America
| | - Molly Przeworski
- Department of Systems Biology, Columbia University, New York, New York, United States of America
- Department of Biological Sciences, Columbia University, New York, New York, United States of America
| |
Collapse
|
22
|
Yuzon JD, Travadon R, Malar C M, Tripathy S, Rank N, Mehl HK, Rizzo DM, Cobb R, Small C, Tang T, McCown HE, Garbelotto M, Kasuga T. Asexual Evolution and Forest Conditions Drive Genetic Parallelism in Phytophthora ramorum. Microorganisms 2020; 8:E940. [PMID: 32580470 PMCID: PMC7357085 DOI: 10.3390/microorganisms8060940] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2020] [Revised: 06/15/2020] [Accepted: 06/18/2020] [Indexed: 11/16/2022] Open
Abstract
It is commonly assumed that asexual lineages are short-lived evolutionarily, yet many asexual organisms can generate genetic and phenotypic variation, providing an avenue for further evolution. Previous work on the asexual plant pathogen Phytophthora ramorum NA1 revealed considerable genetic variation in the form of Structural Variants (SVs). To better understand how SVs arise and their significance to the California NA1 population, we studied the evolutionary histories of SVs and the forest conditions associated with their emergence. Ancestral state reconstruction suggests that SVs arose by somatic mutations among multiple independent lineages, rather than by recombination. We asked if this unusual phenomenon of parallel evolution between isolated populations is transmitted to extant lineages and found that SVs persist longer in a population if their genetic background had a lower mutation load. Genetic parallelism was also found in geographically distant demes where forest conditions such as host density, solar radiation, and temperature, were similar. Parallel SVs overlap with genes involved in pathogenicity such as RXLRs and have the potential to change the course of an epidemic. By combining genomics and environmental data, we identified an unexpected pattern of repeated evolution in an asexual population and identified environmental factors potentially driving this phenomenon.
Collapse
Affiliation(s)
- Jennifer David Yuzon
- Department of Plant Pathology, University of California, Davis, CA 95616, USA; (R.T.); (H.K.M.); (D.M.R.); (C.S.); (T.T.); (H.E.M.)
| | - Renaud Travadon
- Department of Plant Pathology, University of California, Davis, CA 95616, USA; (R.T.); (H.K.M.); (D.M.R.); (C.S.); (T.T.); (H.E.M.)
| | - Mathu Malar C
- CSIR Indian Institute of Chemical Biology, Kolkata 700032, India; (M.M.C.); (S.T.)
| | - Sucheta Tripathy
- CSIR Indian Institute of Chemical Biology, Kolkata 700032, India; (M.M.C.); (S.T.)
| | - Nathan Rank
- Department of Biology, Sonoma State University, Rohnert Park, CA 94928, USA;
| | - Heather K. Mehl
- Department of Plant Pathology, University of California, Davis, CA 95616, USA; (R.T.); (H.K.M.); (D.M.R.); (C.S.); (T.T.); (H.E.M.)
| | - David M. Rizzo
- Department of Plant Pathology, University of California, Davis, CA 95616, USA; (R.T.); (H.K.M.); (D.M.R.); (C.S.); (T.T.); (H.E.M.)
| | - Richard Cobb
- Department of Natural Resources and Environmental Science, California Polytechnic State University, San Luis Obispo, CA 93407, USA;
| | - Corinn Small
- Department of Plant Pathology, University of California, Davis, CA 95616, USA; (R.T.); (H.K.M.); (D.M.R.); (C.S.); (T.T.); (H.E.M.)
| | - Tiffany Tang
- Department of Plant Pathology, University of California, Davis, CA 95616, USA; (R.T.); (H.K.M.); (D.M.R.); (C.S.); (T.T.); (H.E.M.)
| | - Haley E. McCown
- Department of Plant Pathology, University of California, Davis, CA 95616, USA; (R.T.); (H.K.M.); (D.M.R.); (C.S.); (T.T.); (H.E.M.)
| | - Matteo Garbelotto
- Department of Environmental Science, Policy and Management, University of California, Berkeley, CA 94720, USA;
| | - Takao Kasuga
- Crops Pathology and Genetics Research Unit, USDA Agricultural Research Service, Davis, CA 95616, USA
| |
Collapse
|
23
|
|
24
|
Katju V, Bergthorsson U. Old Trade, New Tricks: Insights into the Spontaneous Mutation Process from the Partnering of Classical Mutation Accumulation Experiments with High-Throughput Genomic Approaches. Genome Biol Evol 2019; 11:136-165. [PMID: 30476040 PMCID: PMC6330053 DOI: 10.1093/gbe/evy252] [Citation(s) in RCA: 70] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/22/2018] [Indexed: 12/17/2022] Open
Abstract
Mutations spawn genetic variation which, in turn, fuels evolution. Hence, experimental investigations into the rate and fitness effects of spontaneous mutations are central to the study of evolution. Mutation accumulation (MA) experiments have served as a cornerstone for furthering our understanding of spontaneous mutations for four decades. In the pregenomic era, phenotypic measurements of fitness-related traits in MA lines were used to indirectly estimate key mutational parameters, such as the genomic mutation rate, new mutational variance per generation, and the average fitness effect of mutations. Rapidly emerging next-generating sequencing technology has supplanted this phenotype-dependent approach, enabling direct empirical estimates of the mutation rate and a more nuanced understanding of the relative contributions of different classes of mutations to the standing genetic variation. Whole-genome sequencing of MA lines bears immense potential to provide a unified account of the evolutionary process at multiple levels-the genetic basis of variation, and the evolutionary dynamics of mutations under the forces of selection and drift. In this review, we have attempted to synthesize key insights into the spontaneous mutation process that are rapidly emerging from the partnering of classical MA experiments with high-throughput sequencing, with particular emphasis on the spontaneous rates and molecular properties of different mutational classes in nuclear and mitochondrial genomes of diverse taxa, the contribution of mutations to the evolution of gene expression, and the rate and stability of transgenerational epigenetic modifications. Future advances in sequencing technologies will enable greater species representation to further refine our understanding of mutational parameters and their functional consequences.
Collapse
Affiliation(s)
- Vaishali Katju
- Department of Veterinary Integrative Biosciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX 77843-4458
| | - Ulfar Bergthorsson
- Department of Veterinary Integrative Biosciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX 77843-4458
| |
Collapse
|
25
|
Thomas GWC, Wang RJ, Puri A, Harris RA, Raveendran M, Hughes DST, Murali SC, Williams LE, Doddapaneni H, Muzny DM, Gibbs RA, Abee CR, Galinski MR, Worley KC, Rogers J, Radivojac P, Hahn MW. Reproductive Longevity Predicts Mutation Rates in Primates. Curr Biol 2018; 28:3193-3197.e5. [PMID: 30270182 DOI: 10.1016/j.cub.2018.08.050] [Citation(s) in RCA: 57] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2018] [Revised: 07/26/2018] [Accepted: 08/22/2018] [Indexed: 12/30/2022]
Abstract
Mutation rates vary between species across several orders of magnitude, with larger organisms having the highest per-generation mutation rates. Hypotheses for this pattern typically invoke physiological or population-genetic constraints imposed on the molecular machinery preventing mutations [1]. However, continuing germline cell division in multicellular eukaryotes means that organisms with longer generation times and of larger size will leave more mutations to their offspring simply as a byproduct of their increased lifespan [2, 3]. Here, we deeply sequence the genomes of 30 owl monkeys (Aotus nancymaae) from six multi-generation pedigrees to demonstrate that paternal age is the major factor determining the number of de novo mutations in this species. We find that owl monkeys have an average mutation rate of 0.81 × 10-8 per site per generation, roughly 32% lower than the estimate in humans. Based on a simple model of reproductive longevity that does not require any changes to the mutational machinery, we show that this is the expected mutation rate in owl monkeys. We further demonstrate that our model predicts species-specific mutation rates in other primates, including study-specific mutation rates in humans based on the average paternal age. Our results suggest that variation in life history traits alone can explain variation in the per-generation mutation rate among primates, and perhaps among a wide range of multicellular organisms.
Collapse
Affiliation(s)
- Gregg W C Thomas
- Department of Biology, Indiana University, 107 S. Indiana Avenue, Bloomington, IN 47405, USA; Department of Computer Science, Indiana University, 107 S. Indiana Avenue, Bloomington, IN 47405, USA.
| | - Richard J Wang
- Department of Biology, Indiana University, 107 S. Indiana Avenue, Bloomington, IN 47405, USA
| | - Arthi Puri
- Department of Computer Science, Indiana University, 107 S. Indiana Avenue, Bloomington, IN 47405, USA
| | - R Alan Harris
- Human Genome Sequencing Center, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA
| | - Muthuswamy Raveendran
- Human Genome Sequencing Center, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA
| | - Daniel S T Hughes
- Human Genome Sequencing Center, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA
| | - Shwetha C Murali
- Human Genome Sequencing Center, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA
| | - Lawrence E Williams
- Keeling Center for Comparative Medicine and Research, University of Texas, MD Anderson Cancer Center, 650 Cool Water Drive, Bastrop, TX 78602, USA
| | - Harsha Doddapaneni
- Human Genome Sequencing Center, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA
| | - Donna M Muzny
- Human Genome Sequencing Center, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA
| | - Richard A Gibbs
- Human Genome Sequencing Center, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA
| | - Christian R Abee
- Keeling Center for Comparative Medicine and Research, University of Texas, MD Anderson Cancer Center, 650 Cool Water Drive, Bastrop, TX 78602, USA
| | - Mary R Galinski
- Emory Vaccine Center, Yerkes National Primate Research Center, Emory University, 201 Dowman Drive, Atlanta, GA, USA; Division of Infectious Diseases, Department of Medicine, Emory University, 201 Dowman Drive, Atlanta, GA, USA
| | - Kim C Worley
- Human Genome Sequencing Center, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA
| | - Jeffrey Rogers
- Human Genome Sequencing Center, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA; Department of Molecular and Human Genetics, Baylor College of Medicine, 1 Baylor Plaza, Houston, TX 77030, USA
| | - Predrag Radivojac
- Department of Computer Science, Indiana University, 107 S. Indiana Avenue, Bloomington, IN 47405, USA
| | - Matthew W Hahn
- Department of Biology, Indiana University, 107 S. Indiana Avenue, Bloomington, IN 47405, USA; Department of Computer Science, Indiana University, 107 S. Indiana Avenue, Bloomington, IN 47405, USA.
| |
Collapse
|