1
|
Jin Y, Du X, Jiang C, Ji W, Yang P. Disentangling sources of gene tree discordance for Hordeum species via target-enriched sequencing assays. Mol Phylogenet Evol 2024; 199:108160. [PMID: 39019201 DOI: 10.1016/j.ympev.2024.108160] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2024] [Revised: 07/04/2024] [Accepted: 07/14/2024] [Indexed: 07/19/2024]
Abstract
Hordeum is an economically and evolutionarily important genus within the Triticeae tribe of the family Poaceae, and contains 33 widely distributed and diverse species which cytologically represent four subgenomes (H, Xa, Xu and I). These wild species (except Hordeum spontaneum, which is the primary gene pool of barley) are secondary or tertiary gene-pool germplasms for barley and wheat improvement, and uncovering their complicated evolutionary relationships would benefit for future breeding programs. Here, we developed a complexity-reduced pipeline via capturing genome-wide distributed fragments via two novel target-enriched assays (HorCap v1.0 and BarPlex v1.0) in conjugation with high-throughput sequencing of the enrichments. Both assays were tested for genotyping 40 species from three genera (Hordeum, Triticum, and Aegilops) containing 82 samples 67 accessions. Either of both assays worked efficiently in genotyping, while integration of both assays can significantly improve the robustness and resolution of the Hordeum phylogenetic trees. Interestingly, the incomplete lineage sorting (ILS) was inferred for the first time as the major factor causing phylogenetic discordance among the four subgenomes, whereas in New World species (carrying I genome) post-speciation introgression events were revealed. Through revising the evolutionary relationships of the Hordeum species based on an ancestral state reconstruction for the diploids and parental donor inference for the polyploids, our results raised new queries about the Hordeum phylogeny. Moreover, both newly-developed assays are applicable in genotyping and phylogenetic analysis of Hordeum and other Triticeae wild species.
Collapse
Affiliation(s)
- Yanlong Jin
- State Key Laboratory of Crop Gene Resources and Breeding, Key Laboratory of Grain Crop Genetic Resources Evaluation and Utilization (MARA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China; State Key Laboratory of Crop Stress Biology for Arid Areas, College of Agronomy, Northwest AandF University, Yangling 712100, China
| | - Xin Du
- State Key Laboratory of Crop Stress Biology for Arid Areas, College of Agronomy, Northwest AandF University, Yangling 712100, China
| | - Congcong Jiang
- State Key Laboratory of Crop Gene Resources and Breeding, Key Laboratory of Grain Crop Genetic Resources Evaluation and Utilization (MARA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Wanquan Ji
- State Key Laboratory of Crop Stress Biology for Arid Areas, College of Agronomy, Northwest AandF University, Yangling 712100, China
| | - Ping Yang
- State Key Laboratory of Crop Gene Resources and Breeding, Key Laboratory of Grain Crop Genetic Resources Evaluation and Utilization (MARA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China.
| |
Collapse
|
2
|
Gutiérrez-Larruscain D, Vargas P, Fernández-Mazuecos M, Pausas JG. Phylogenomic analysis reveals the evolutionary history of Paleartic needle-leaved junipers. Mol Phylogenet Evol 2024; 199:108162. [PMID: 39067655 DOI: 10.1016/j.ympev.2024.108162] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Revised: 07/10/2024] [Accepted: 07/24/2024] [Indexed: 07/30/2024]
Abstract
Needle-leaved junipers (Juniperus sect. Juniperus, Cupressaceae) are coniferous trees and shrubs with red or blue fleshy cones. They are distributed across Asia, Macaronesia and the Mediterranean Basin, with one species (J. communis) having a circumboreal distribution. Here we aim to resolve the phylogeny of this clade to infer its intricate evolutionary history. To do so, we built a comprehensive, time-calibrated phylogeny using genotyping-by-sequencing (GBS) and combine it with species occurrence using phylogeographic tools. Our results provide solid phylogenetic resolution to propose a new taxonomic classification and a biogeographical history of the section. Specifically, we confirm the monophyly of two groups within J. sect. Juniperus: the Asian (blue-cone) species including the circumboreal J. communis, and the Mediterranean-Macaronesian (red-cone) species. In addition, we provide strong phylogenetic evidence for three distinct species (J. badia, J. conferta, J. lutchuensis) previously considered subspecies or varieties, as well as for the differentiation between the eastern and western Mediterranean lineages of J. macrocarpa. Our findings suggest that the Mediterranean basin was the primary center of diversification for Juniperus sect. Juniperus, followed by an East Asian-Tethyan disjunction resulting from uplifts of the Qinghai-Tibetan Plateau and climatic shifts. The colonization history of Macaronesia by red-cone junipers from the western Mediterranean appears to have taken place independently in two different geological periods: the Miocene (Azores) and the Pliocene (Madeira-Canary Islands). Overall, genomic data and phylogenetic analysis are key to consider a new taxonomic proposal and reconstruct the biogeographical history of the iconic needle-leaved junipers across the Paleartic.
Collapse
Affiliation(s)
- David Gutiérrez-Larruscain
- Department of Ecology and Global Change, Desertification Research Centre (CIDE: CSIC-UV-GVA), Valencia, Spain.
| | - Pablo Vargas
- Department of Biodiversity and Conservation, Real Jardín Botánico (RJB: CSIC), Madrid, Spain
| | - Mario Fernández-Mazuecos
- Department of Biology (Botany), Faculty of Sciences, Universidad Autónoma de Madrid (UAM), Madrid, Spain
| | - Juli G Pausas
- Department of Ecology and Global Change, Desertification Research Centre (CIDE: CSIC-UV-GVA), Valencia, Spain
| |
Collapse
|
3
|
Campbell MA, Hammer MP, Adams M, Raadik TA, Unmack PJ. Evolutionary relationships and fine-scale geographic structuring in the temperate percichthyid genus Gadopsis (blackfishes) to support fisheries and conservation management. Mol Phylogenet Evol 2024; 199:108159. [PMID: 39029548 DOI: 10.1016/j.ympev.2024.108159] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2024] [Revised: 07/04/2024] [Accepted: 07/14/2024] [Indexed: 07/21/2024]
Abstract
Gadopsis (Percichthyidae) is a freshwater genus distributed in south-eastern Australia, including Tasmania, and comprises two recognized species. Previous molecular phylogenetic investigations of the genus, mostly conducted in the pre-genomics era and reflecting a range of geographic and molecular sampling intensities, have supported the recognition of up to seven candidate species. Here we analyze a genome-wide SNP dataset that provides comprehensive geographic and genomic coverage of Gadopsis to produce a robust hypothesis of species boundaries and evolutionary relationships. We then leverage the SNP dataset to characterize relationships within candidate species that lack clear intraspecific phylogenetic relationships. We find further support for the seven previously identified candidate species of Gadopsis and evidence that the Bass Strait centered candidate species (SBA) originated from ancient hybridization. The SNP dataset permits a high degree of intraspecific resolution, providing improvements over previous studies, with numerous candidate species showing intraspecific divisions in phylogenetic analysis. Further population genetic analysis of the Murray-Darling candidate species (NMD) and SBA finds support for K = 6 and K = 7 genetic clusters, respectively. The SNP data generated for this study have diverse applications in natural resource management for these fishes of conservation concern.
Collapse
Affiliation(s)
- Matthew A Campbell
- The University of California Davis, Davis, California, USA; University of Alaska Museum of the North, Fairbanks, Alaska, USA.
| | - Michael P Hammer
- Museum and Art Gallery of the Northern Territory, Darwin, Northern Territory, Australia
| | - Mark Adams
- South Australian Museum, Adelaide, South Australia, Australia; School of Biological Sciences, University of Adelaide, Adelaide, South Australia, Australia
| | - Tarmo A Raadik
- Arthur Rylah Institute for Environmental Research, Department of Energy, Environment and Climate Action, Heidelberg, Victoria, Australia
| | - Peter J Unmack
- University of Canberra, Canberra, Australian Capital Territory, Australia
| |
Collapse
|
4
|
Bjornson S, Verbruggen H, Upham NS, Steenwyk JL. Reticulate evolution: Detection and utility in the phylogenomics era. Mol Phylogenet Evol 2024; 201:108197. [PMID: 39270765 DOI: 10.1016/j.ympev.2024.108197] [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: 05/12/2024] [Revised: 08/13/2024] [Accepted: 09/08/2024] [Indexed: 09/15/2024]
Abstract
Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.
Collapse
Affiliation(s)
- Saelin Bjornson
- School of BioSciences, University of Melbourne, Victoria, Australia
| | - Heroen Verbruggen
- School of BioSciences, University of Melbourne, Victoria, Australia; CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Campus de Vairão, Universidade do Porto, 4485-661 Vairão, Portugal
| | - Nathan S Upham
- School of Life Sciences, Arizona State University, Tempe, AZ, USA.
| | - Jacob L Steenwyk
- Howards Hughes Medical Institute and the Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, CA, USA.
| |
Collapse
|
5
|
Pezzi PH, Wheeler LC, Freitas LB, Smith SD. Incomplete lineage sorting and hybridization underlie tree discordance in Petunia and related genera (Petunieae, Solanaceae). Mol Phylogenet Evol 2024; 198:108136. [PMID: 38909873 DOI: 10.1016/j.ympev.2024.108136] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2024] [Revised: 06/06/2024] [Accepted: 06/17/2024] [Indexed: 06/25/2024]
Abstract
Despite the overarching history of species divergence, phylogenetic studies often reveal distinct topologies across regions of the genome. The sources of these gene tree discordances are variable, but incomplete lineage sorting (ILS) and hybridization are among those with the most biological importance. Petunia serves as a classic system for studying hybridization in the wild. While field studies suggest that hybridization is frequent, the extent of reticulation within Petunia and its closely related genera has never been examined from a phylogenetic perspective. In this study, we used transcriptomic data from 11 Petunia, 16 Calibrachoa, and 10 Fabiana species to illuminate the relationships between these species and investigate whether hybridization played a significant role in the diversification of the clade. We inferred that gene tree discordance within genera is linked to hybridization events along with high levels of ILS due to their rapid diversification. Moreover, network analyses estimated deeper hybridization events between Petunia and Calibrachoa, genera that have different chromosome numbers. Although these genera cannot hybridize at the present time, ancestral hybridization could have played a role in their parallel radiations, as they share the same habitat and life history.
Collapse
Affiliation(s)
- Pedro H Pezzi
- Department of Genetics, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil.
| | - Lucas C Wheeler
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, USA
| | - Loreta B Freitas
- Department of Genetics, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Stacey D Smith
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, USA
| |
Collapse
|
6
|
Wang Y, Wu X, Chen Y, Xu C, Wang Y, Wang Q. Phylogenomic analyses revealed widely occurring hybridization events across Elsholtzieae (Lamiaceae). Mol Phylogenet Evol 2024; 198:108112. [PMID: 38806075 DOI: 10.1016/j.ympev.2024.108112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2024] [Revised: 05/14/2024] [Accepted: 05/22/2024] [Indexed: 05/30/2024]
Abstract
Obtaining a robust phylogeny proves challenging due to the intricate evolutionary history of species, where processes such as hybridization and incomplete lineage sorting can introduce conflicting signals, thereby complicating phylogenetic inference. In this study, we conducted comprehensive sampling of Elsholtzieae, with a particular focus on its largest genus, Elsholtzia. We utilized 503 nuclear loci and complete plastome sequences obtained from 99 whole-genome sequencing datasets to elucidate the interspecific relationships within the Elsholtzieae. Additionally, we explored various sources of conflicts between gene trees and species trees. Fully supported backbone phylogenies were recovered, and the monophyly of Elsholtzia and Keiskea was not supported. Significant gene tree heterogeneity was observed at numerous nodes, particularly regarding the placement of Vuhuangia and the E. densa clade. Further investigations into potential causes of this discordance revealed that incomplete lineage sorting (ILS), coupled with hybridization events, has given rise to substantial gene tree discordance. Several species, represented by multiple samples, exhibited a closer association with geographical distribution rather than following a strictly monophyletic pattern in plastid trees, suggesting chloroplast capture within Elsholtzieae and providing evidence of hybridization. In conclusion, this study provides phylogenomic insights to untangle taxonomic problems in the tribe Elsholtzieae, especially the genus Elsholtzia.
Collapse
Affiliation(s)
- Yan Wang
- State Key Laboratory of Plant Diversity and Specialty Crops, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; National Botanical Garden, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xuexue Wu
- State Key Laboratory of Plant Diversity and Specialty Crops, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; National Botanical Garden, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yanyi Chen
- State Key Laboratory of Plant Diversity and Specialty Crops, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; National Botanical Garden, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Chao Xu
- State Key Laboratory of Plant Diversity and Specialty Crops, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; National Botanical Garden, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
| | - Yinghui Wang
- State Key Laboratory of Plant Diversity and Specialty Crops, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; National Botanical Garden, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Qiang Wang
- State Key Laboratory of Plant Diversity and Specialty Crops, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; National Botanical Garden, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China.
| |
Collapse
|
7
|
Döcker J, Linz S, Wicke K. Bounding the Softwired Parsimony Score of a Phylogenetic Network. Bull Math Biol 2024; 86:121. [PMID: 39174812 PMCID: PMC11341636 DOI: 10.1007/s11538-024-01350-9] [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: 05/29/2024] [Accepted: 08/16/2024] [Indexed: 08/24/2024]
Abstract
In comparison to phylogenetic trees, phylogenetic networks are more suitable to represent complex evolutionary histories of species whose past includes reticulation such as hybridisation or lateral gene transfer. However, the reconstruction of phylogenetic networks remains challenging and computationally expensive due to their intricate structural properties. For example, the small parsimony problem that is solvable in polynomial time for phylogenetic trees, becomes NP-hard on phylogenetic networks under softwired and parental parsimony, even for a single binary character and structurally constrained networks. To calculate the parsimony score of a phylogenetic network N, these two parsimony notions consider different exponential-size sets of phylogenetic trees that can be extracted from N and infer the minimum parsimony score over all trees in the set. In this paper, we ask: What is the maximum difference between the parsimony score of any phylogenetic tree that is contained in the set of considered trees and a phylogenetic tree whose parsimony score equates to the parsimony score of N? Given a gap-free sequence alignment of multi-state characters and a rooted binary level-k phylogenetic network, we use the novel concept of an informative blob to show that this difference is bounded by k + 1 times the softwired parsimony score of N. In particular, the difference is independent of the alignment length and the number of character states. We show that an analogous bound can be obtained for the softwired parsimony score of semi-directed networks, while under parental parsimony on the other hand, such a bound does not hold.
Collapse
Affiliation(s)
- Janosch Döcker
- School of Computer Science, University of Auckland, Auckland, New Zealand
| | - Simone Linz
- School of Computer Science, University of Auckland, Auckland, New Zealand.
| | - Kristina Wicke
- Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, USA
| |
Collapse
|
8
|
Du K, Ricci JMB, Lu Y, Garcia-Olazabal M, Walter RB, Warren WC, Dodge TO, Schumer M, Park H, Meyer A, Schartl M. Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation. Nat Commun 2024; 15:6609. [PMID: 39098897 PMCID: PMC11298535 DOI: 10.1038/s41467-024-50852-6] [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: 12/31/2023] [Accepted: 07/16/2024] [Indexed: 08/06/2024] Open
Abstract
Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.
Collapse
Affiliation(s)
- Kang Du
- The Xiphophorus Genetic Stock Center, Texas State University, San Marcos, Texas, TX, USA
| | | | - Yuan Lu
- The Xiphophorus Genetic Stock Center, Texas State University, San Marcos, Texas, TX, USA
| | - Mateo Garcia-Olazabal
- The Xiphophorus Genetic Stock Center, Texas State University, San Marcos, Texas, TX, USA
| | - Ronald B Walter
- The Xiphophorus Genetic Stock Center, Texas State University, San Marcos, Texas, TX, USA
| | - Wesley C Warren
- Department of Animal Sciences, Department of Surgery, Institute for Data Science and Informatics, University of Missouri, Bond Life Sciences Center, Columbia, MI, USA
| | - Tristram O Dodge
- Department of Biology & Howard Hughes Medical Institute, Stanford University, Stanford, CA, USA
| | - Molly Schumer
- Department of Biology & Howard Hughes Medical Institute, Stanford University, Stanford, CA, USA
| | - Hyun Park
- Division of Biotechnology, College of Life Science and Biotechnology, Korea University, Seoul, Republic of Korea
| | - Axel Meyer
- Department of Biology, University of Konstanz, Konstanz, Germany.
| | - Manfred Schartl
- The Xiphophorus Genetic Stock Center, Texas State University, San Marcos, Texas, TX, USA.
- Developmental Biochemistry, Biocenter, University of Wuerzburg, Am Hubland, Wuerzburg, Germany.
- Research Department for Limnology, University of Innsbruck, Mondsee, Austria.
| |
Collapse
|
9
|
Weinell JL, Burbrink FT, Das S, Brown RM. Novel phylogenomic inference and 'Out of Asia' biogeography of cobras, coral snakes and their allies. ROYAL SOCIETY OPEN SCIENCE 2024; 11:240064. [PMID: 39113776 PMCID: PMC11303032 DOI: 10.1098/rsos.240064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Revised: 05/27/2024] [Accepted: 05/31/2024] [Indexed: 08/10/2024]
Abstract
Estimation of evolutionary relationships among lineages that rapidly diversified can be challenging, and, in such instances, inaccurate or unresolved phylogenetic estimates can lead to erroneous conclusions regarding historical geographical ranges of lineages. One example underscoring this issue has been the historical challenge posed by untangling the biogeographic origin of elapoid snakes, which includes numerous dangerously venomous species as well as species not known to be dangerous to humans. The worldwide distribution of this lineage makes it an ideal group for testing hypotheses related to historical faunal exchanges among the many continents and other landmasses occupied by contemporary elapoid species. We developed a novel suite of genomic resources, included worldwide sampling, and inferred a robust estimate of evolutionary relationships, which we leveraged to quantitatively estimate geographical range evolution through the deep-time history of this remarkable radiation. Our phylogenetic and biogeographical estimates of historical ranges definitively reject a lingering former 'Out of Africa' hypothesis and support an 'Out of Asia' scenario involving multiple faunal exchanges between Asia, Africa, Australasia, the Americas and Europe.
Collapse
Affiliation(s)
- Jeffrey L. Weinell
- Department of Ecology and Evolutionary Biology and Biodiversity Institute, University of Kansas, 1345 Jayhawk Blvd, Lawrence, KS66045, USA
- Department of Herpetology, American Museum of Natural History, 200 Central Park West, New York, NY10024, USA
| | - Frank T. Burbrink
- Department of Herpetology, American Museum of Natural History, 200 Central Park West, New York, NY10024, USA
| | - Sunandan Das
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki00014, Finland
| | - Rafe M. Brown
- Department of Ecology and Evolutionary Biology and Biodiversity Institute, University of Kansas, 1345 Jayhawk Blvd, Lawrence, KS66045, USA
| |
Collapse
|
10
|
Wang M, Zhu M, Qian J, Yang Z, Shang F, Egan AN, Li P, Liu L. Phylogenomics of mulberries (Morus, Moraceae) inferred from plastomes and single copy nuclear genes. Mol Phylogenet Evol 2024; 197:108093. [PMID: 38740145 DOI: 10.1016/j.ympev.2024.108093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2023] [Revised: 04/30/2024] [Accepted: 05/06/2024] [Indexed: 05/16/2024]
Abstract
Mulberries (genus Morus), belonging to the order Rosales, family Moraceae, are important woody plants due to their economic values in sericulture, as well as for nutritional benefits and medicinal values. However, the taxonomy and phylogeny of Morus, especially for the Asian species, remains challenging due to its wide geographical distribution, morphological plasticity, and interspecific hybridization. To better understand the evolutionary history of Morus, we combined plastomes and a large-scale nuclear gene analyses to investigate their phylogenetic relationships. We assembled the plastomes and screened 211 single-copy nuclear genes from 13 Morus species and related taxa. The plastomes of Morus species were relatively conserved in terms of genome size, gene content, synteny, IR boundary and codon usage. Using nuclear data, our results elucidated identical topologies based on coalescent and concatenation methods. The genus Morus was supported as monophyletic, with M. notabilis as the first diverging lineage and the two North American Morus species, M. microphylla and M. rubra, as sister to the other Asian species. In the Asian Morus species, interspecific relationships were completely resolved. However, cyto-nuclear discordances and gene tree-species tree conflicts were detected in the phylogenies of Morus, with multiple evidences supporting hybridization/introgression as the main cause of discordances between nuclear and plastid phylogenies, while gene tree-species tree conflicts were mainly caused by ILS.
Collapse
Affiliation(s)
- Meizhen Wang
- College of Life Sciences, Henan Normal University, Xinxiang 453000, China; Systematic & Evolutionary Botany and Biodiversity Group, MOE Key Laboratory of Biosystems Homeostasis and Protection, College of Life Sciences, Zhejiang University, Hangzhou 310058, China
| | - Mengmeng Zhu
- Laboratory of Plant Germplasm and Genetic Engineering, School of Life Sciences, Henan University, Kaifeng 475001, China
| | - Jiayi Qian
- College of Life Sciences, Henan Agricultural University, Zhengzhou 450002, China
| | - Zhaoping Yang
- College of Life Sciences and Technologies, Tarim University, Alar 843300, China
| | - Fude Shang
- Laboratory of Plant Germplasm and Genetic Engineering, School of Life Sciences, Henan University, Kaifeng 475001, China; College of Life Sciences, Henan Agricultural University, Zhengzhou 450002, China.
| | - Ashley N Egan
- Department of Biology, Utah Valley University, Orem, UT 84058, United States.
| | - Pan Li
- Systematic & Evolutionary Botany and Biodiversity Group, MOE Key Laboratory of Biosystems Homeostasis and Protection, College of Life Sciences, Zhejiang University, Hangzhou 310058, China.
| | - Luxian Liu
- College of Life Sciences, Henan Normal University, Xinxiang 453000, China; Laboratory of Plant Germplasm and Genetic Engineering, School of Life Sciences, Henan University, Kaifeng 475001, China.
| |
Collapse
|
11
|
Allman ES, Baños H, Garrote-Lopez M, Rhodes JA. Identifiability of Level-1 Species Networks from Gene Tree Quartets. Bull Math Biol 2024; 86:110. [PMID: 39052074 PMCID: PMC11272829 DOI: 10.1007/s11538-024-01339-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/22/2024] [Accepted: 07/01/2024] [Indexed: 07/27/2024]
Abstract
When hybridization or other forms of lateral gene transfer have occurred, evolutionary relationships of species are better represented by phylogenetic networks than by trees. While inference of such networks remains challenging, several recently proposed methods are based on quartet concordance factors-the probabilities that a tree relating a gene sampled from the species displays the possible 4-taxon relationships. Building on earlier results, we investigate what level-1 network features are identifiable from concordance factors under the network multispecies coalescent model. We obtain results on both topological features of the network, and numerical parameters, uncovering a number of failures of identifiability related to 3-cycles in the network. Addressing these identifiability issues is essential for designing statistically consistent inference methods.
Collapse
Affiliation(s)
- Elizabeth S Allman
- Department of Mathematics and Statistics, University of Alaska, Fairbanks, AK, USA.
| | - Hector Baños
- Department of Mathematics, California State University San Bernadino, San Bernadino, CA, USA
| | | | - John A Rhodes
- Department of Mathematics and Statistics, University of Alaska, Fairbanks, AK, USA
| |
Collapse
|
12
|
Valdés-Florido A, González-Toral C, Maguilla E, Cires E, Díaz-Lifante Z, Andrés-Camacho C, Nieto Feliner G, Arroyo J, Escudero M. Polyploidy and hybridization in the Mediterranean: unravelling the evolutionary history of Centaurium (Gentianaceae). ANNALS OF BOTANY 2024; 134:247-262. [PMID: 38687133 PMCID: PMC11232519 DOI: 10.1093/aob/mcae066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/16/2024] [Accepted: 04/29/2024] [Indexed: 05/02/2024]
Abstract
BACKGROUND AND AIMS Polyploidy is considered one of the main mechanisms of plant evolution and speciation. In the Mediterranean Basin, polyploidy has contributed to making this region a biodiversity hotspot, along with its geological and climatic history and other ecological and biogeographical factors. The Mediterranean genus Centaurium (Gentianaceae) comprises ~25 species, of which 60 % are polyploids, including tetraploids and hexaploids. To date, the evolutionary history of centauries has been studied using Sanger sequencing phylogenies, which have been insufficient to fully understand the phylogenetic relationships in this lineage. The goal of this study is to gain a better understanding of the evolutionary history of Centaurium by exploring the mechanisms that have driven its diversification, specifically hybridization and polyploidy. We aim to identify the parentage of hybrid species, at the species or clade level, as well as assessing whether morphological traits are associated with particular ploidy levels. METHODS We sequenced RADseq markers from 42 samples of 28 Centaurium taxa, and performed phylogenomic analyses using maximum likelihood, summary coalescent SVDquartets and Neighbor-Net approaches. To identify hybrid taxa, we used PhyloNetworks and the fastSTRUCTURE algorithm. To infer the putative parental species of the allopolyploids, we employed genomic analyses (SNIPloid). The association between different traits and particular ploidy levels was explored with non-metric multidimensional scaling. KEY RESULTS Our phylogenetic analyses confirmed the long-suspected occurrence of recurrent hybridization. The allopolyploid origin of the tetraploid C. serpentinicola and the hexaploids C. mairei, C. malzacianum and C. centaurioides was also confirmed, unlike that of C. discolor. We inferred additional signatures of hybridization events within the genus and identified morphological traits differentially distributed in different ploidy levels. CONCLUSIONS This study highlights the important role that hybridization has played in the evolution of a Mediterranean genus such as Centaurium, leading to a polyploid complex, which facilitated its diversification and may exemplify that of other Mediterranean groups.
Collapse
Affiliation(s)
- Ana Valdés-Florido
- Department of Plant Biology and Ecology, Faculty of Biology, University of Seville, Seville, 41012, Spain
| | | | - Enrique Maguilla
- Department of Molecular Biology and Biochemical Engineering, Pablo de Olavide University, Seville, 41013, Spain
| | - Eduardo Cires
- Department of Organisms and Systems Biology, University of Oviedo, Oviedo, 33071, Spain
- Institute of Natural Resources and Territorial Planning (INDUROT), Campus de Mieres, Mieres, 33600, Spain
| | - Zoila Díaz-Lifante
- Department of Plant Biology and Ecology, Faculty of Biology, University of Seville, Seville, 41012, Spain
| | - Cristina Andrés-Camacho
- Department of Plant Biology and Ecology, Faculty of Biology, University of Seville, Seville, 41012, Spain
| | | | - Juan Arroyo
- Department of Plant Biology and Ecology, Faculty of Biology, University of Seville, Seville, 41012, Spain
| | - Marcial Escudero
- Department of Plant Biology and Ecology, Faculty of Biology, University of Seville, Seville, 41012, Spain
| |
Collapse
|
13
|
Rose JP, Kriebel R, Sytsma KJ, Drew BT. Phylogenomic perspectives on speciation and reproductive isolation in a North American biodiversity hotspot: an example using California sages (Salvia subgenus Audibertia: Lamiaceae). ANNALS OF BOTANY 2024; 134:295-310. [PMID: 38733329 PMCID: PMC11232522 DOI: 10.1093/aob/mcae073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/16/2024] [Accepted: 05/07/2024] [Indexed: 05/13/2024]
Abstract
BACKGROUND AND AIMS The California Floristic Province (CA-FP) is the most species-rich region of North America north of Mexico. One of several proposed hypotheses explaining the exceptional diversity of the region is that the CA-FP harbours myriad recently diverged lineages with nascent reproductive barriers. Salvia subgenus Audibertia is a conspicuous element of the CA-FP, with multiple sympatric and compatible species. METHODS Using 305 nuclear loci and both organellar genomes, we reconstruct species trees, examine genomic discordance, conduct divergence-time estimation, and analyse contemporaneous patterns of gene flow and mechanical reproductive isolation. KEY RESULTS Despite strong genomic discordance, an underlying bifurcating tree is supported. Organellar genomes capture additional introgression events not detected in the nuclear genome. Most interfertility is found within clades, indicating that reproductive barriers arise with increasing genetic divergence. Species are generally not mechanically isolated, suggesting that it is unlikely to be the primary factor leading to reproductive isolation. CONCLUSIONS Rapid, recent speciation with some interspecific gene flow in conjunction with the onset of a Mediterranean-like climate is the underlying cause of extant diversity in Salvia subgenus Audibertia. Speciation has largely not been facilitated by gene flow. Its signal in the nuclear genome seems to mostly be erased by backcrossing, but organellar genomes each capture different instances of historical gene flow, probably characteristic of many CA-FP lineages. Mechanical reproductive isolation appears to be only part of a mosaic of factors limiting gene flow.
Collapse
Affiliation(s)
- Jeffrey P Rose
- Department of Biology, University of Nebraska at Kearney, Kearney, NE 68849, USA
- Department of Botany, University of Wisconsin-Madison, 430 Lincoln Drive, Madison, WI 53706, USA
| | - Ricardo Kriebel
- Department of Botany, University of Wisconsin-Madison, 430 Lincoln Drive, Madison, WI 53706, USA
- California Academy of Sciences, San Francisco, CA 94118, USA
| | - Kenneth J Sytsma
- Department of Botany, University of Wisconsin-Madison, 430 Lincoln Drive, Madison, WI 53706, USA
| | - Bryan T Drew
- Department of Biology, University of Nebraska at Kearney, Kearney, NE 68849, USA
| |
Collapse
|
14
|
Przelomska NAS, Diaz RA, Ávila FA, Ballen GA, Cortés-B R, Kistler L, Chitwood DH, Charitonidou M, Renner SS, Pérez-Escobar OA, Antonelli A. Morphometrics and Phylogenomics of Coca (Erythroxylum spp.) Illuminate Its Reticulate Evolution, With Implications for Taxonomy. Mol Biol Evol 2024; 41:msae114. [PMID: 38982580 PMCID: PMC11233275 DOI: 10.1093/molbev/msae114] [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: 08/08/2023] [Revised: 05/01/2024] [Accepted: 05/10/2024] [Indexed: 07/11/2024] Open
Abstract
South American coca (Erythroxylum coca and E. novogranatense) has been a keystone crop for many Andean and Amazonian communities for at least 8,000 years. However, over the last half-century, global demand for its alkaloid cocaine has driven intensive agriculture of this plant and placed it in the center of armed conflict and deforestation. To monitor the changing landscape of coca plantations, the United Nations Office on Drugs and Crime collects annual data on their areas of cultivation. However, attempts to delineate areas in which different varieties are grown have failed due to limitations around identification. In the absence of flowers, identification relies on leaf morphology, yet the extent to which this is reflected in taxonomy is uncertain. Here, we analyze the consistency of the current naming system of coca and its four closest wild relatives (the "coca clade"), using morphometrics, phylogenomics, molecular clocks, and population genomics. We include name-bearing type specimens of coca's closest wild relatives E. gracilipes and E. cataractarum. Morphometrics of 342 digitized herbarium specimens show that leaf shape and size fail to reliably discriminate between species and varieties. However, the statistical analyses illuminate that rounder and more obovate leaves of certain varieties could be associated with the subtle domestication syndrome of coca. Our phylogenomic data indicate extensive gene flow involving E. gracilipes which, combined with morphometrics, supports E. gracilipes being retained as a single species. Establishing a robust evolutionary-taxonomic framework for the coca clade will facilitate the development of cost-effective genotyping methods to support reliable identification.
Collapse
Affiliation(s)
- Natalia A S Przelomska
- School of Biological Sciences, University of Portsmouth, Portsmouth PO1 2DY, UK
- Royal Botanic Gardens, Kew, Richmond, Surrey TW9 3AE, UK
- Department of Anthropology, National Museum of Natural History, Smithsonian Institution, Washington DC 20560, USA
| | - Rudy A Diaz
- Royal Botanic Gardens, Kew, Richmond, Surrey TW9 3AE, UK
| | | | - Gustavo A Ballen
- Instituto de Biociências, Universidade Estadual Paulista, Botucatu, São Paulo, Brazil
- School of Biological and Behavioural Sciences, Queen Mary University of London, London E1 4NS, UK
| | - Rocío Cortés-B
- Herbario Forestal Universidad Distrital, Campus El Vivero, CR 5E 15-82 Bogotá, Colombia
| | - Logan Kistler
- Department of Anthropology, National Museum of Natural History, Smithsonian Institution, Washington DC 20560, USA
| | - Daniel H Chitwood
- Department of Horticulture, Michigan State University, East Lansing, MI 48824, USA
- Department of Computational Mathematics, Science & Engineering, Michigan State University, East Lansing, MI 48824, USA
| | - Martha Charitonidou
- Department of Biological Applications and Technology, University of Ioannina, 45110 Ioannina, Greece
| | - Susanne S Renner
- Department of Biology, Washington University, Saint Louis, MO 63130, USA
| | | | - Alexandre Antonelli
- Royal Botanic Gardens, Kew, Richmond, Surrey TW9 3AE, UK
- Gothenburg Global Biodiversity Centre, Department of Biological and Environmental Sciences, University of Gothenburg, SE 41319 Göteborg, Sweden
- Department of Biology, University of Oxford, Oxford OX1 3RB, UK
| |
Collapse
|
15
|
Bacon CD, Hill A. Hybridization in palms (Arecaceae). Ecol Evol 2024; 14:e70014. [PMID: 39011137 PMCID: PMC11246834 DOI: 10.1002/ece3.70014] [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: 09/13/2023] [Revised: 07/01/2024] [Accepted: 07/01/2024] [Indexed: 07/17/2024] Open
Abstract
Hybridization has significant evolutionary consequences across the Tree of Life. The process of hybridization has played a major role in plant evolution and has contributed to species richness and trait variation. Since morphological traits are partially a product of their environment, there may be a link between hybridization and ecology. Plant hybrid species richness is noted to be higher in harsh environments, and we explore this hypothesis with a keystone tropical plant lineage, palms (Arecaceae). Leveraging a recent literature review of naturally occurring palm hybrids, we developed a method to calculate hybrid frequency, and then tested if there is phylogenetic signal of hybrids using a phylogeny of all palms. Further, we used phylogenetic comparative methods to examine the interaction between hybrid frequency and presence in dry environments, on islands, and the species richness of genera. Phylogenetic generalized least squares models had stronger support than models of random association, indicating phylogenetic signal for the presence of hybrids in dry and island environments. However, all p-values were >.05 and therefore the correlation was poor between hybridization and the trait frequencies examined. Presence in particular environments are not strongly correlated to hybrid frequency, but phylogenetic signal suggests a role in its distribution in different habitats. Hybridization in palms is not evenly distributed across subfamilies, tribes, subtribes yet plays an important role in palm diversity, nonetheless. Increasing our understanding hybridization in this economically and culturally important plant family is essential, particularly since rates are projected to increase with climate change, reconfiguring the dynamics and distribution of biodiversity.
Collapse
Affiliation(s)
- Christine D. Bacon
- Department of Biological and Environmental SciencesUniversity of GothenburgGothenburgSweden
- Gothenburg Global Biodiversity CentreGothenburgSweden
| | - Adrian Hill
- Department of Biological and Environmental SciencesUniversity of GothenburgGothenburgSweden
- Gothenburg Global Biodiversity CentreGothenburgSweden
| |
Collapse
|
16
|
Mo YK, Hahn MW, Smith ML. Applications of machine learning in phylogenetics. Mol Phylogenet Evol 2024; 196:108066. [PMID: 38565358 DOI: 10.1016/j.ympev.2024.108066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2023] [Revised: 02/16/2024] [Accepted: 03/21/2024] [Indexed: 04/04/2024]
Abstract
Machine learning has increasingly been applied to a wide range of questions in phylogenetic inference. Supervised machine learning approaches that rely on simulated training data have been used to infer tree topologies and branch lengths, to select substitution models, and to perform downstream inferences of introgression and diversification. Here, we review how researchers have used several promising machine learning approaches to make phylogenetic inferences. Despite the promise of these methods, several barriers prevent supervised machine learning from reaching its full potential in phylogenetics. We discuss these barriers and potential paths forward. In the future, we expect that the application of careful network designs and data encodings will allow supervised machine learning to accommodate the complex processes that continue to confound traditional phylogenetic methods.
Collapse
Affiliation(s)
- Yu K Mo
- Department of Computer Science, Indiana University, Bloomington, IN 47405, USA
| | - Matthew W Hahn
- Department of Computer Science, Indiana University, Bloomington, IN 47405, USA; Department of Biology, Indiana University, Bloomington, IN 47405, USA
| | - Megan L Smith
- Department of Biological Sciences, Mississippi State University, Starkville, MS 39762, USA.
| |
Collapse
|
17
|
Franco FF, Amaral DT, Bonatelli IAS, Meek JB, Moraes EM, Zappi DC, Taylor NP, Eaton DAR. A historical stepping-stone path for an island-colonizing cactus across a submerged "bridge" archipelago. Heredity (Edinb) 2024; 132:296-308. [PMID: 38637723 PMCID: PMC11166651 DOI: 10.1038/s41437-024-00683-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: 07/29/2023] [Revised: 03/26/2024] [Accepted: 03/27/2024] [Indexed: 04/20/2024] Open
Abstract
Here we use population genomic data (ddRAD-Seq) and ecological niche modeling to test biogeographic hypotheses for the divergence of the island-endemic cactus species Cereus insularis Hemsl. (Cereeae; Cactaceae) from its sister species C. fernambucensis Lem. The Cereus insularis grows in the Fernando de Noronha Islands (FNI), a Neotropical archipelago located 350 km off the Brazilian Atlantic Forest (BAF) coast. Phylogeographic reconstructions support a northward expansion by the common ancestor of C. insularis and C. fernambucensis along the mainland BAF coast, with C. insularis diverging from the widespread mainland taxon C. fernambucensis after colonizing FNI in the late Pleistocene. The morphologically distinct C. insularis is monophyletic and nested within C. fernambucensis, as expected from a progenitor-derivative speciation model. We tested alternative biogeographic and demographic hypotheses for the colonization of the FNI using Approximate Bayesian Computation. We found the greatest support for a stepping-stone path that emerged during periods of decreased sea level (the "bridge" hypothesis), in congruence with historical ecological niche modeling that shows highly suitable habitats on stepping-stone islands during glacial periods. The outlier analyses reveal signatures of selection in C. insularis, suggesting a putative role of adaptation driving rapid anagenic differentiation of this species in FNI.
Collapse
Affiliation(s)
- Fernando Faria Franco
- Departamento de Biologia. Centro de Ciências Humanas e Biológicas, Universidade Federal de São Carlos (UFSCar), Sorocaba, Brazil.
| | - Danilo Trabuco Amaral
- Centro de Ciências Naturais e Humanas, Universidade Federal do ABC (UFABC), Santo André, São Paulo, Brazil
- Programa de Pós Graduação em Biologia Comparada. Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto, Universidade de São Paulo (USP), Ribeirão Preto, Brazil
| | - Isabel A S Bonatelli
- Instituto de Ciências Ambientais, Químicas e Farmacêuticas. Departamento de Ecologia e Biologia Evolutiva, Universidade Federal de São Paulo, Diadema, SP, Brazil
| | - Jared B Meek
- Department of Ecology, Evolution and Environmental Biology, Columbia University, New York, NY, 10027, USA
| | - Evandro Marsola Moraes
- Departamento de Biologia. Centro de Ciências Humanas e Biológicas, Universidade Federal de São Carlos (UFSCar), Sorocaba, Brazil
| | - Daniela Cristina Zappi
- Programa de Pós Graduação em Botânica, Instituto de Ciências Biológicas, Universidade de Brasília, PO Box 04457, Brasília, DF, 70910970, Brazil
| | - Nigel Paul Taylor
- Departamento de Biologia. Centro de Ciências Humanas e Biológicas, Universidade Federal de São Carlos (UFSCar), Sorocaba, Brazil
| | - Deren A R Eaton
- Department of Ecology, Evolution and Environmental Biology, Columbia University, New York, NY, 10027, USA
| |
Collapse
|
18
|
Liu G, Pan Q, Dai Y, Wang X, Li M, Zhu P, Zhou X. Phylogenomics of Afrotherian mammals and improved resolution of extant Paenungulata. Mol Phylogenet Evol 2024; 195:108047. [PMID: 38460890 DOI: 10.1016/j.ympev.2024.108047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Revised: 02/19/2024] [Accepted: 03/06/2024] [Indexed: 03/11/2024]
Abstract
Molecular investigations have gathered a diverse set of mammals-predominantly African natives like elephants, hyraxes, and aardvarks-into a clade known as Afrotheria. Nevertheless, the precise phylogenetic relationships among these species remain contentious. Here, we sourced orthologous markers and ultraconserved elements to discern the interordinal connections among Afrotherian mammals. Our phylogenetic analyses bolster the common origin of Afroinsectiphilia and Paenungulata, and propose Afrosoricida as the closer relative to Macroscelidea rather than Tubulidentata, while also challenging the notion of Sirenia and Hyracoidea as sister taxa. The approximately unbiased test and the gene concordance factor uniformly recognized the alliance of Proboscidea with Hyracoidea as the dominant topology within Paenungulata. Investigation into sites with extremly high phylogenetic signal unveiled their potential to intensify conflicts in the Paenungulata topology. Subsequent exploration suggested that incomplete lineage sorting was predominantly responsible for the observed contentious relationships, whereas introgression exerted a subsidiary influence. The divergence times estimated in our study hint at the Cretaceous-Paleogene (K-Pg) extinction event as a catalyst for Afrotherian diversification. Overall, our findings deliver a tentative but insightful overview of Afrotheria phylogeny and divergence, elucidating these relationships through the lens of phylogenomics.
Collapse
Affiliation(s)
- Gaoming Liu
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Qi Pan
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yichen Dai
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Xiao Wang
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Meng Li
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Pingfen Zhu
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Xuming Zhou
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China.
| |
Collapse
|
19
|
Rancilhac L, Enbody ED, Harris R, Saitoh T, Irestedt M, Liu Y, Lei F, Andersson L, Alström P. Introgression Underlies Phylogenetic Uncertainty But Not Parallel Plumage Evolution in a Recent Songbird Radiation. Syst Biol 2024; 73:12-25. [PMID: 37801684 PMCID: PMC11129591 DOI: 10.1093/sysbio/syad062] [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: 04/26/2023] [Revised: 09/11/2023] [Accepted: 10/05/2023] [Indexed: 10/08/2023] Open
Abstract
Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails," a group of 5 species in the songbird genus Motacilla: 1 species, Motacilla alba, shows wide intra-specific plumage variation, while the 4r others form 2 pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the 2 species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome to quantify the extent of discordance across genomic regions and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investing genome-wide patterns of gene tree heterogeneity to help understand the mechanisms underlying phenotypic evolution. [Gene tree heterogeneity; incomplete lineage sorting; introgression; parallel evolution; phylogenomics; plumage evolution; wagtails.].
Collapse
Affiliation(s)
- Loïs Rancilhac
- Animal Ecology, Department of Ecology and Genetics, Evolutionary Biology Centre, Uppsala University, Norbyvägen 18 D, 752 36 Uppsala, Sweden
| | - Erik D Enbody
- Department of Medical Biochemistry and Microbiology, Uppsala University, 751 23 Uppsala, Sweden
- Biomolecular Engineering, University of California, 95064 Santa Cruz, CA, USA
| | - Rebecca Harris
- Department of Biology, University of Washington, Seattle, WA 98105, USA
| | - Takema Saitoh
- Yamashina Institute for Ornithology, 115 Konoyama, Abiko, Chiba 270-1145, Japan
| | - Martin Irestedt
- Department of Bioinformatics and Genetics, Swedish Museum of Natural History, P.O. Box 50007, 104 05 Stockholm, Sweden
| | - Yang Liu
- State Key Laboratory of Biocontrol, School of Ecology, Sun Yat-sen University, Shenzhen 518107, China
| | - Fumin Lei
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
| | - Leif Andersson
- Department of Medical Biochemistry and Microbiology, Uppsala University, 751 23 Uppsala, Sweden
- Department of Veterinary Integrative Biosciences, Texas A&M University, College Station, TX 77843, USA
| | - Per Alström
- Animal Ecology, Department of Ecology and Genetics, Evolutionary Biology Centre, Uppsala University, Norbyvägen 18 D, 752 36 Uppsala, Sweden
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
| |
Collapse
|
20
|
Pang XX, Zhang DY. Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information. Syst Biol 2024; 73:207-222. [PMID: 38224495 PMCID: PMC11129598 DOI: 10.1093/sysbio/syad077] [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: 05/01/2023] [Revised: 12/17/2023] [Accepted: 12/27/2023] [Indexed: 01/17/2024] Open
Abstract
In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.
Collapse
Affiliation(s)
- Xiao-Xu Pang
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, China
| | - Da-Yong Zhang
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
21
|
Jewett EM. CORRECTING MODEL MISSPECIFICATION IN RELATIONSHIP ESTIMATES. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.13.594005. [PMID: 38868169 PMCID: PMC11167690 DOI: 10.1101/2024.05.13.594005] [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: 06/14/2024]
Abstract
The datasets of large genotyping biobanks and direct-to-consumer genetic testing companies contain many related individuals. Until now, it has been widely accepted that the most distant relationships that can be detected are around fifteen degrees (approximately 8 th cousins) and that practical relationship estimates have a ceiling around ten degrees (approximately 5 th cousins). However, we show that these assumptions are incorrect and that they are due to a misapplication of relationship estimators. In particular, relationship estimators are applied almost exclusively to putative relatives who have been identified because they share detectable tracts of DNA identically by descent (IBD). However, no existing relationship estimator conditions on the event that two individuals share at least one detectable segment of IBD anywhere in the genome. As a result, the relationship estimates obtained using existing estimators are dramatically biased for distant relationships, inferring all sufficiently distant relationships to be around ten degrees regardless of the depth of the true relationship. Moreover, existing relationship estimators are derived under a model that assumes that each pair of related individuals shares a single common ancestor (or mating pair of ancestors). This model breaks down for relationships beyond 10 generations in the past because individuals share many thousands of cryptic common ancestors due to pedigree collapse. We first derive a corrected likelihood that conditions on the event that at least one segment is observed between a pair of putative relatives and we demonstrate that the corrected likelihood largely eliminates the bias in estimates of pairwise relationships and provides a more accurate characterization of the uncertainty in these estimates. We then reformulate the relationship inference problem to account for the fact that individuals share many common ancestors, not just one. We demonstrate that the most distant relationship that can be inferred may be forty degrees or more, rather than ten, extending the time-to-common ancestor from approximately 200 years in the past to approximately 600 years in the past or more. This dramatic increase in the range of relationship estimators makes it possible to infer relationships whose common ancestors lived before historical events such as European settlement of the Americas and the Transatlantic Slave Trade, and possibly much earlier.
Collapse
|
22
|
Morales-Saldaña S, Hipp AL, Valencia-Ávalos S, Hahn M, González-Elizondo MS, Gernandt DS, Pham KK, Oyama K, González-Rodríguez A. Divergence and reticulation in the Mexican white oaks: ecological and phylogenomic evidence on species limits and phylogenetic networks in the Quercus laeta complex (Fagaceae). ANNALS OF BOTANY 2024; 133:1007-1024. [PMID: 38428030 PMCID: PMC11089265 DOI: 10.1093/aob/mcae030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 02/28/2024] [Indexed: 03/03/2024]
Abstract
BACKGROUND AND AIMS Introgressive hybridization poses a challenge to taxonomic and phylogenetic understanding of taxa, particularly when there are high numbers of co-occurring, intercrossable species. The genus Quercus exemplifies this situation. Oaks are highly diverse in sympatry and cross freely, creating syngameons of interfertile species. Although a well-resolved, dated phylogeny is available for the American oak clade, evolutionary relationships within many of the more recently derived clades remain to be defined, particularly for the young and exceptionally diverse Mexican white oak clade. Here, we adopted an approach bridging micro- and macroevolutionary scales to resolve evolutionary relationships in a rapidly diversifying clade endemic to Mexico. METHODS Ecological data and sequences of 155 low-copy nuclear genes were used to identify distinct lineages within the Quercus laeta complex. Concatenated and coalescent approaches were used to assess the phylogenetic placement of these lineages relative to the Mexican white oak clade. Phylogenetic network methods were applied to evaluate the timing and genomic significance of recent or historical introgression among lineages. KEY RESULTS The Q. laeta complex comprises six well-supported lineages, each restricted geographically and with mostly divergent climatic niches. Species trees corroborated that the different lineages are more closely related to other species of Mexican white oaks than to each other, suggesting that this complex is polyphyletic. Phylogenetic networks estimated events of ancient introgression that involved the ancestors of three present-day Q. laeta lineages. CONCLUSIONS The Q. laeta complex is a morphologically and ecologically related group of species rather than a clade. Currently, oak phylogenetics is at a turning point, at which it is necessary to integrate phylogenetics and ecology in broad regional samples to figure out species boundaries. Our study illuminates one of the more complicated of the Mexican white oak groups and lays groundwork for further taxonomic study.
Collapse
Affiliation(s)
- Saddan Morales-Saldaña
- Instituto de Investigaciones en Ecosistemas y Sustentabilidad, Universidad Nacional Autónoma de México (UNAM), Antigua Carretera a Pátzcuaro No. 8701, Col. Ex-Hacienda de San José de la Huerta, Morelia, 58190, Michoacán, México
| | - Andrew L Hipp
- The Morton Arboretum, Lisle, IL 60532-1293, USA
- The Field Museum, Chicago, IL 60605, USA
| | - Susana Valencia-Ávalos
- Herbario de la Facultad de Ciencias, Departamento de Biología Comparada, Universidad Nacional Autónoma de México (UNAM), 04510, Ciudad de México, México
| | | | | | - David S Gernandt
- Departamento de Botánica, Instituto de Biología, Universidad Nacional Autónoma de México (UNAM), 04510, Ciudad de México, México
| | - Kasey K Pham
- Department of Biology, University of Florida, Gainesville, FL 32611, USA
| | - Ken Oyama
- Escuela Nacional de Estudios Superiores Unidad Morelia, Universidad Nacional Autónoma de México (UNAM), Antigua Carretera a Pátzcuaro No. 8701, Col. Ex‐Hacienda de San José de la Huerta, Morelia, 58190, Michoacán, México
| | - Antonio González-Rodríguez
- Instituto de Investigaciones en Ecosistemas y Sustentabilidad, Universidad Nacional Autónoma de México (UNAM), Antigua Carretera a Pátzcuaro No. 8701, Col. Ex-Hacienda de San José de la Huerta, Morelia, 58190, Michoacán, México
| |
Collapse
|
23
|
Myers EA, Rautsaw RM, Borja M, Jones J, Grünwald CI, Holding ML, Grazziotin F, Parkinson CL. Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus). Syst Biol 2024:syae018. [PMID: 38695290 DOI: 10.1093/sysbio/syae018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Indexed: 08/11/2024] Open
Abstract
Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.
Collapse
Affiliation(s)
- Edward A Myers
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA
- Department of Herpetology, California Academy of Sciences, San Francisco, CA 94118, USA
| | - Rhett M Rautsaw
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA
| | - Miguel Borja
- Facultad de Ciencias Biológicas, Universdad Juárez del Estado de Durango, Av. Universidad s/n. Fracc. Filadelfia, Gómez Palacio, Durango., Mex
| | - Jason Jones
- Herp.mx A.C. C.P. 28989, Villa de Álvarez, Colima, Mexico
| | - Christoph I Grünwald
- Herp.mx A.C. C.P. 28989, Villa de Álvarez, Colima, Mexico
- Biodiversa A.C., Avenida de la Ribera #203, C.P. 45900, Chapala, Jalisco, Mexico
| | - Matthew L Holding
- Department of Biological Sciences, Clemson University, Clemson, SC 29634, USA
- Life Sciences Institute, University of Michigan, Ann Arbor, MI 48109, USA
| | - Felipe Grazziotin
- Laboratório Especial de Coleções Zoológicas, Instituto Butantan, Avenida Vital Brasil, São Paulo, 05503-900, Brazil
| | | |
Collapse
|
24
|
Rurik I, Melichárková A, Gbúrová Štubová E, Kučera J, Kochjarová J, Paun O, Vďačný P, Slovák M. Homoplastic versus xenoplastic evolution: exploring the emergence of key intrinsic and extrinsic traits in the montane genus Soldanella (Primulaceae). THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2024; 118:753-765. [PMID: 38217489 DOI: 10.1111/tpj.16630] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 12/02/2023] [Accepted: 12/27/2023] [Indexed: 01/15/2024]
Abstract
Specific ecological conditions in the high mountain environment exert a selective pressure that often leads to convergent trait evolution. Reticulations induced by incomplete lineage sorting and introgression can lead to discordant trait patterns among gene and species trees (hemiplasy/xenoplasy), providing a false illusion that the traits under study are homoplastic. Using phylogenetic species networks, we explored the effect of gene exchange on trait evolution in Soldanella, a genus profoundly influenced by historical introgression. At least three features evolved independently multiple times: the single-flowered dwarf phenotype, dysploid cytotype, and ecological generalism. The present analyses also indicated that the recurring occurrence of stoloniferous growth might have been prompted by an introgression event between an ancestral lineage and a still extant species, although its emergence via convergent evolution cannot be completely ruled out. Phylogenetic regression suggested that the independent evolution of larger genomes in snowbells is most likely a result of the interplay between hybridization events of dysploid and euploid taxa and hostile environments at the range margins of the genus. The emergence of key intrinsic and extrinsic traits in snowbells has been significantly impacted not only by convergent evolution but also by historical and recent introgression events.
Collapse
Affiliation(s)
- Ivan Rurik
- Department of Zoology, Comenius University Bratislava, Ilkovičova 6, 842 15, Bratislava, Slovak Republic
| | - Andrea Melichárková
- Institute of Botany, Plant Science and Biodiversity Centre, Slovak Academy of Sciences, Dúbravská cesta 9, 845 23, Bratislava, Slovak Republic
| | - Eliška Gbúrová Štubová
- Institute of Botany, Plant Science and Biodiversity Centre, Slovak Academy of Sciences, Dúbravská cesta 9, 845 23, Bratislava, Slovak Republic
- Slovak National Museum, Natural History Museum, Vajanského nábrežie 2, 810 06, Bratislava, Slovak Republic
| | - Jaromír Kučera
- Institute of Botany, Plant Science and Biodiversity Centre, Slovak Academy of Sciences, Dúbravská cesta 9, 845 23, Bratislava, Slovak Republic
| | - Judita Kochjarová
- Department of Phytology, Faculty of Forestry, Technical University Zvolen, Masarykova 24, 960 53, Zvolen, Slovak Republic
| | - Ovidiu Paun
- Department of Botany and Biodiversity Research, University of Vienna, Rennweg 14, 1030, Vienna, Austria
| | - Peter Vďačný
- Department of Zoology, Comenius University Bratislava, Ilkovičova 6, 842 15, Bratislava, Slovak Republic
| | - Marek Slovák
- Institute of Botany, Plant Science and Biodiversity Centre, Slovak Academy of Sciences, Dúbravská cesta 9, 845 23, Bratislava, Slovak Republic
- Department of Botany, Charles University, Benátská 2, 128 01, Prague, Czech Republic
| |
Collapse
|
25
|
Das S, Greenbaum E, Brecko J, Pauwels OSG, Ruane S, Pirro S, Merilä J. Phylogenomics of Psammodynastes and Buhoma (Elapoidea: Serpentes), with the description of a new Asian snake family. Sci Rep 2024; 14:9489. [PMID: 38664489 PMCID: PMC11045840 DOI: 10.1038/s41598-024-60215-2] [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: 07/23/2023] [Accepted: 04/19/2024] [Indexed: 04/28/2024] Open
Abstract
Asian mock vipers of the genus Psammodynastes and African forest snakes of the genus Buhoma are two genera belonging to the snake superfamily Elapoidea. The phylogenetic placements of Psammodynastes and Buhoma within Elapoidea has been extremely unstable which has resulted in their uncertain and debated taxonomy. We used ultraconserved elements and traditional nuclear and mitochondrial markers to infer the phylogenetic relationships of these two genera with other elapoids. Psammodynastes, for which a reference genome has been sequenced, were found, with strong branch support, to be a relatively early diverging split within Elapoidea that is sister to a clade consisting of Elapidae, Micrelapidae and Lamprophiidae. Hence, we allocate Psammodynastes to its own family, Psammodynastidae new family. However, the phylogenetic position of Buhoma could not be resolved with a high degree of confidence. Attempts to identify the possible sources of conflict in the rapid radiation of elapoid snakes suggest that both hybridisation/introgression during the rapid diversification, including possible ghost introgression, as well as incomplete lineage sorting likely have had a confounding role. The usual practice of combining mitochondrial loci with nuclear genomic data appears to mislead phylogeny reconstructions in rapid radiation scenarios, especially in the absence of genome scale data.
Collapse
Affiliation(s)
- Sunandan Das
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, 00014, Helsinki, Finland.
| | - Eli Greenbaum
- Department of Biological Sciences, University of Texas at El Paso, 500 W. University Avenue, El Paso, TX, 79968, USA
| | - Jonathan Brecko
- Royal Belgian Institute of Natural Sciences, Rue Vautier 29, 1000, Brussels, Belgium
- Royal Museum for Central Africa, Tervuren, Belgium
| | - Olivier S G Pauwels
- Royal Belgian Institute of Natural Sciences, Rue Vautier 29, 1000, Brussels, Belgium
| | - Sara Ruane
- Life Sciences Section, Negaunee Integrative Research Center, Field Museum, Chicago, IL, USA
| | - Stacy Pirro
- Iridian Genomes Inc., Bethesda, MD, 20817, USA
| | - Juha Merilä
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, 00014, Helsinki, Finland
- Area of Ecology and Biodiversity, School of Biological Sciences, Kadoorie Biological Sciences Building, The University of Hong Kong, Pokfulam Road, Hong Kong, SAR, China
| |
Collapse
|
26
|
Allman ES, Baños H, Mitchell JD, Rhodes JA. TINNiK: Inference of the Tree of Blobs of a Species Network Under the Coalescent. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.04.20.590418. [PMID: 38712257 PMCID: PMC11071406 DOI: 10.1101/2024.04.20.590418] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2024]
Abstract
The tree of blobs of a species network shows only the tree-like aspects of relationships of taxa on a network, omitting information on network substructures where hybridization or other types of lateral transfer of genetic information occur. By isolating such regions of a network, inference of the tree of blobs can serve as a starting point for a more detailed investigation, or indicate the limit of what may be inferrable without additional assumptions. Building on our theoretical work on the identifiability of the tree of blobs from gene quartet distributions under the Network Multispecies Coalescent model, we develop an algorithm, TINNiK, for statistically consistent tree of blobs inference. We provide examples of its application to both simulated and empirical datasets, utilizing an implementation in the MSCquartets 2.0 R package. MSC Classification 92D15, 92D20.
Collapse
|
27
|
Roycroft E, Ford F, Ramm T, Schembri R, Breed WG, Burns PA, Rowe KC, Moritz C. Speciation across biomes: Rapid diversification with reproductive isolation in the Australian delicate mice. Mol Ecol 2024; 33:e17301. [PMID: 38385302 DOI: 10.1111/mec.17301] [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/17/2023] [Revised: 01/25/2024] [Accepted: 01/30/2024] [Indexed: 02/23/2024]
Abstract
Phylogeographic studies of continental clades, especially when combined with palaeoclimate modelling, provide powerful insight into how environment drives speciation across climatic contexts. Australia, a continent characterized by disparate modern biomes and dynamic climate change, provides diverse opportunity to reconstruct the impact of past and present environments on diversification. Here, we use genomic-scale data (1310 exons and whole mitogenomes from 111 samples) to investigate Pleistocene diversification, cryptic diversity, and secondary contact in the Australian delicate mice (Hydromyini: Pseudomys), a recent radiation spanning almost all Australian environments. Across northern Australia, we find no evidence for introgression between cryptic lineages within Pseudomys delicatulus sensu lato, with palaeoclimate models supporting contraction and expansion of suitable habitat since the last glacial maximum. Despite multiple contact zones, we also find little evidence of introgression at a continental scale, with the exception of a potential hybrid zone in the mesic biome. In the arid zone, combined insights from genetic data and palaeomodels support a recent expansion in the arid specialist P. hermannsburgensis and contraction in the semi-arid P. bolami. In the face of repeated secondary contact, differences in sperm morphology and chromosomal rearrangements are potential mechanisms that maintain species boundaries in these recently diverged species. Additionally, we describe the western delicate mouse as a new species and recommend taxonomic reinstatement of the eastern delicate mouse. Overall, we show that speciation in an evolutionarily young and widespread clade has been driven by environmental change, and potentially maintained by divergence in reproductive morphology and chromosome rearrangements.
Collapse
Affiliation(s)
- Emily Roycroft
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Acton, Australian Capital Territory, Australia
| | - Fred Ford
- Biodiversity Conservation and Science, New South Wales Department of Planning and Environment, Queanbeyan, New South Wales, Australia
- Australian National Wildlife Collection, National Research Collections Australia, CSIRO, Canberra, Australian Capital Territory, Australia
| | - Till Ramm
- Zoo Leipzig, Leipzig, Germany
- Museum für Naturkunde Berlin, Berlin, Germany
| | - Rhiannon Schembri
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Acton, Australian Capital Territory, Australia
- School of Natural Sciences, Macquarie University, Macquarie Park, New South Wales, Australia
| | - William G Breed
- School of Biological Sciences, University of Adelaide, Adelaide, South Australia, Australia
| | - Phoebe A Burns
- Wildlife Conservation and Science, Zoos Victoria, Parkville, Victoria, Australia
| | - Kevin C Rowe
- Sciences Department, Museums Victoria, Melbourne, Victoria, Australia
- School of Biosciences, The University of Melbourne, Parkville, Victoria, Australia
| | - Craig Moritz
- Division of Ecology and Evolution, Research School of Biology, The Australian National University, Acton, Australian Capital Territory, Australia
| |
Collapse
|
28
|
Wang MZ, Wu J, Zhang SL, Mao LM, Ohi-Toma T, Takano A, Zhang YH, Cameron KM, Li P. Species delimitation in Amana (Liliaceae): transcriptomes battle with evolutionary complexity. Cladistics 2024; 40:135-156. [PMID: 37983640 DOI: 10.1111/cla.12565] [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: 12/14/2022] [Revised: 10/10/2023] [Accepted: 10/18/2023] [Indexed: 11/22/2023] Open
Abstract
Species delimitation has long been a subject of controversy, and there are many alternative concepts and approaches used to define species in plants. The genus Amana (Liliaceae), known as "East Asian tulips" has a number of cryptic species and a huge genome size (1C = 21.48-57.35 pg). It also is intriguing how such a spring ephemeral genus thrives in subtropical areas. However, phylogenetic relationships and species delimitation within Amana are challenging. Here we included all species and 84 populations of Amana, which are collected throughout its distribution range. A variety of methods were used to clarify its species relationships based on a combination of morphological, ecological, genetic, evolutionary and phylogenetic species concepts. This evidence supports the recognition of at least 12 species in Amana. Moreover, we explored the complex evolutionary history within the genus and detected several historical hybridization and introgression events based on phylogenetic trees (transcriptomic and plastid), phylonetworks, admixture and ABBA-BABA analyses. Morphological traits have undergone parallel evolution in the genus. This spring ephemeral genus might have originated from a temperate region, yet finally thrives in subtropical areas, and three hypotheses about its adaptive evolution are proposed for future testing. In addition, we propose a new species, Amana polymorpha, from eastern Zhejiang Province, China. This research also demonstrates that molecular evidence at the genome level (such as transcriptomes) has greatly improved the accuracy and reasonability of species delimitation and taxon classification.
Collapse
Affiliation(s)
- Mei-Zhen Wang
- Laboratory of Systematic & Evolutionary Botany and Biodiversity, College of Life Sciences, Zhejiang University, Hangzhou, 310058, Zhejiang, China
| | - Jing Wu
- Laboratory of Systematic & Evolutionary Botany and Biodiversity, College of Life Sciences, Zhejiang University, Hangzhou, 310058, Zhejiang, China
| | - Sheng-Lu Zhang
- Plant Quarantine Station of Lin'an District, Hangzhou, 311300, Zhejiang, China
| | - Li-Mi Mao
- Nanjing Institute of Geology and Palaeontology, Chinese Academy of Sciences, Nanjing, 210008, Jiangsu, China
| | - Tetsuo Ohi-Toma
- Nature Fieldwork Center, Okayama University of Science, Okayama, 700-0005, Japan
| | - Atsuko Takano
- Museum of Nature and Human Activities, Hyogo 6 chome, Yayoigaoka, Sanda, Hyogo, 669-1546, Japan
| | - Yong-Hua Zhang
- College of Life and Environmental Science, Wenzhou University, Wenzhou, 325035, China
| | - Kenneth M Cameron
- Department of Botany, University of Wisconsin, Madison, WI, 53706, USA
| | - Pan Li
- Laboratory of Systematic & Evolutionary Botany and Biodiversity, College of Life Sciences, Zhejiang University, Hangzhou, 310058, Zhejiang, China
| |
Collapse
|
29
|
Bedoya AM. Botany and geogenomics: Constraining geological hypotheses in the neotropics with large-scale genetic data derived from plants. AMERICAN JOURNAL OF BOTANY 2024; 111:e16306. [PMID: 38557829 DOI: 10.1002/ajb2.16306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2024] [Accepted: 01/24/2024] [Indexed: 04/04/2024]
Abstract
Decades of empirical research have revealed how the geological history of our planet shaped plant evolution by establishing well-known patterns (e.g., how mountain uplift resulted in high rates of diversification and replicate radiations in montane plant taxa). This follows a traditional approach where botanical data are interpreted in light of geological events. In this synthesis, I instead describe how by integrating natural history, phylogenetics, and population genetics, botanical research can be applied alongside geology and paleontology to inform our understanding of past geological and climatic processes. This conceptual shift aligns with the goals of the emerging field of geogenomics. In the neotropics, plant geogenomics is a powerful tool for the reciprocal exploration of two long standing questions in biology and geology: how the dynamic landscape of the region came to be and how it shaped the evolution of the richest flora. Current challenges that are specific to analytical approaches for plant geogenomics are discussed. I describe the scale at which various geological questions can be addressed from biological data and what makes some groups of plants excellent model systems for geogenomics research. Although plant geogenomics is discussed with reference to the neotropics, the recommendations given here for approaches to plant geogenomics can and should be expanded to exploring long-standing questions on how the earth evolved with the use of plant DNA.
Collapse
Affiliation(s)
- Ana M Bedoya
- Department of Biological Sciences, Louisiana State University, Baton Rouge, 70803, LA, USA
| |
Collapse
|
30
|
Ané C, Fogg J, Allman ES, Baños H, Rhodes JA. Anomalous networks under the multispecies coalescent: theory and prevalence. J Math Biol 2024; 88:29. [PMID: 38372830 DOI: 10.1007/s00285-024-02050-7] [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: 07/15/2023] [Revised: 01/18/2024] [Accepted: 01/21/2024] [Indexed: 02/20/2024]
Abstract
Reticulations in a phylogenetic network represent processes such as gene flow, admixture, recombination and hybrid speciation. Extending definitions from the tree setting, an anomalous network is one in which some unrooted tree topology displayed in the network appears in gene trees with a lower frequency than a tree not displayed in the network. We investigate anomalous networks under the Network Multispecies Coalescent Model with possible correlated inheritance at reticulations. Focusing on subsets of 4 taxa, we describe a new algorithm to calculate quartet concordance factors on networks of any level, faster than previous algorithms because of its focus on 4 taxa. We then study topological properties required for a 4-taxon network to be anomalous, uncovering the key role of [Formula: see text]-cycles: cycles of 3 edges parent to a sister group of 2 taxa. Under the model of common inheritance, that is, when each gene tree coalesces within a species tree displayed in the network, we prove that 4-taxon networks are never anomalous. Under independent and various levels of correlated inheritance, we use simulations under realistic parameters to quantify the prevalence of anomalous 4-taxon networks, finding that truly anomalous networks are rare. At the same time, however, we find a significant fraction of networks close enough to the anomaly zone to appear anomalous, when considering the quartet concordance factors observed from a few hundred genes. These apparent anomalies may challenge network inference methods.
Collapse
Affiliation(s)
- Cécile Ané
- Department of Statistics, University of Wisconsin - Madison, Madison, WI, 53706, USA.
- Department of Botany, University of Wisconsin - Madison, Madison, WI, 53706, USA.
| | - John Fogg
- Department of Statistics, University of Wisconsin - Madison, Madison, WI, 53706, USA
| | - Elizabeth S Allman
- Department of Mathematics and Statistics, University of Alaska Fairbanks, Fairbanks, AK, 99775-6660, USA
| | - Hector Baños
- Department of Biochemistry & Molecular Biology, Dalhousie University, Halifax, NS, Canada
- Department of Mathematics and Statistics, Dalhousie University, Halifax, NS, Canada
| | - John A Rhodes
- Department of Mathematics and Statistics, University of Alaska Fairbanks, Fairbanks, AK, 99775-6660, USA
| |
Collapse
|
31
|
Wu Z, Solís-Lemus C. Ultrafast learning of four-node hybridization cycles in phylogenetic networks using algebraic invariants. BIOINFORMATICS ADVANCES 2024; 4:vbae014. [PMID: 38384862 PMCID: PMC10879748 DOI: 10.1093/bioadv/vbae014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Revised: 12/23/2023] [Accepted: 02/06/2024] [Indexed: 02/23/2024]
Abstract
Motivation The abundance of gene flow in the Tree of Life challenges the notion that evolution can be represented with a fully bifurcating process which cannot capture important biological realities like hybridization, introgression, or horizontal gene transfer. Coalescent-based network methods are increasingly popular, yet not scalable for big data, because they need to perform a heuristic search in the space of networks as well as numerical optimization that can be NP-hard. Here, we introduce a novel method to reconstruct phylogenetic networks based on algebraic invariants. While there is a long tradition of using algebraic invariants in phylogenetics, our work is the first to define phylogenetic invariants on concordance factors (frequencies of four-taxon splits in the input gene trees) to identify level-1 phylogenetic networks under the multispecies coalescent model. Results Our novel hybrid detection methodology is optimization-free as it only requires the evaluation of polynomial equations, and as such, it bypasses the traversal of network space, yielding a computational speed at least 10 times faster than the fastest-to-date network methods. We illustrate our method's performance on simulated and real data from the genus Canis. Availability and implementation We present an open-source publicly available Julia package PhyloDiamond.jl available at https://github.com/solislemuslab/PhyloDiamond.jl with broad applicability within the evolutionary community.
Collapse
Affiliation(s)
- Zhaoxing Wu
- Department of Statistics, Wisconsin Institute for Discovery, University of Wisconsin-Madison, Madison, WI 53706, United States
| | - Claudia Solís-Lemus
- Department of Plant Pathology, Wisconsin Institute for Discovery, University of Wisconsin-Madison, Madison, WI 53706, United States
| |
Collapse
|
32
|
Kolb AW, Chau VQ, Miller DL, Yannuzzi NA, Brandt CR. Phylogenetic and Recombination Analysis of Clinical Vitreous Humor-Derived Adenovirus Isolates Reveals Discordance Between Serotype and Phylogeny. Invest Ophthalmol Vis Sci 2024; 65:12. [PMID: 38319669 PMCID: PMC10854415 DOI: 10.1167/iovs.65.2.12] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Accepted: 01/19/2024] [Indexed: 02/07/2024] Open
Abstract
Purpose To sequence, identify, and perform phylogenetic and recombination analysis on three clinical adenovirus samples taken from the vitreous humor at the Bascom Palmer Eye Institute. Methods The PacBio Sequel II was used to sequence the genomes of the three clinical adenovirus isolates. To identify the isolates, a full genome-based multiple sequence alignment (MSA) of 722 mastadenoviruses was generated using multiple alignment using fast Fourier transform (MAFFT). MAFFT was also used to generate genome-based human adenovirus B (HAdV-B) MSAs, as well as HAdV-B fiber, hexon, and penton protein-based MSAs. To examine recombination within HAdV-B, RF-Net 2 and Bootscan software programs were used. Results In the course of classifying three new atypical ocular adenovirus samples, taken from the vitreous humor, we found that all three isolates were HAdV-B species. The three Bascom Palmer HAdV-B genomes were then combined with over 300 HAdV-B genome sequences, including nine ocular HAdV-B genome sequences. Attempts to categorize the penton, hexon, and fiber serotypes using phylogeny of the three Bascom Palmer samples were inconclusive due to incongruence between serotype and phylogeny in the dataset. Recombination analysis using a subset of HAdV-B strains to generate a hybridization network detected recombination between nonhuman primate and human-derived strains, recombination between one HAdV-B strain and the HAdV-E outgroup, and limited recombination between the B1 and B2 clades. Conclusions The discordance between serotype and phylogeny detected in this study suggests that the current classification system does not accurately describe the natural history and phylogenetic relationships among adenoviruses.
Collapse
Affiliation(s)
- Aaron W. Kolb
- Department of Ophthalmology and Visual Sciences, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin, United States
| | - Viet Q. Chau
- Bascom Palmer Eye Institute, University of Miami, Miami, Florida, United States
| | - Darlene L. Miller
- Bascom Palmer Eye Institute, University of Miami, Miami, Florida, United States
| | - Nicolas A. Yannuzzi
- Bascom Palmer Eye Institute, University of Miami, Miami, Florida, United States
| | - Curtis R. Brandt
- Department of Ophthalmology and Visual Sciences, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin, United States
- McPherson Eye Research Institute, University of Wisconsin-Madison, Madison, Wisconsin, United States
- Department of Medical Microbiology and Immunology, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin, United States
| |
Collapse
|
33
|
Nevado B, Atchison GW, Bridges EL, Orzell S, Filatov D, Hughes CE. Pleistocene diversification of unifoliolate-leaved Lupinus (Leguminosae: Papilionoideae) in Florida. Mol Ecol 2024; 33:e17232. [PMID: 38205900 DOI: 10.1111/mec.17232] [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: 08/02/2023] [Revised: 11/16/2023] [Accepted: 11/21/2023] [Indexed: 01/12/2024]
Abstract
The importance and prevalence of recent ice-age and post-glacial speciation and species diversification during the Pleistocene across many organismal groups and physiographic settings are well established. However, the extent to which Pleistocene diversification can be attributed to climatic oscillations and their effects on distribution ranges and population structure remains debatable. In this study, we use morphologic, geographic and genetic (RADseq) data to document Pleistocene speciation and intra-specific diversification of the unifoliolate-leaved clade of Florida Lupinus, a small group of species largely restricted to inland and coastal sand ridges across the Florida peninsula and panhandle. Phylogenetic and demographic analyses alongside morphological and geographic evidence suggest that recent speciation and intra-specific divergence within this clade were driven by a combination of non-adaptive allopatric divergence caused by edaphic niche conservatism and opportunities presented by the emergence of new post-glacial sand ridge habitats. These results highlight the central importance of even modest geographic isolation and short periods of allopatric divergence following range expansion in the emergence of new taxa and add to the growing evidence that Pleistocene climatic oscillations may contribute to rapid diversification in a myriad of physiographic settings. Furthermore, our results shed new light on long-standing taxonomic debate surrounding the number of species in the Florida unifoliate Lupinus clade providing support for recognition of five species and a set of intra-specific variants. The important conservation implications for the narrowly restricted, highly endangered species Lupinus aridorum, which we show to be genetically distinct from its sister species Lupinus westianus, are discussed.
Collapse
Affiliation(s)
- Bruno Nevado
- Faculty of Sciences, cE3c - Centre for Ecology, Evolution and Environmental Changes & CHANGE - Global Change and Sustainability Institute, University of Lisbon, Lisbon, Portugal
- Department of Animal Biology, Faculty of Sciences, University of Lisbon, Lisbon, Portugal
| | - Guy W Atchison
- Department of Systematic and Evolutionary Botany, University of Zurich, Zurich, Switzerland
| | - Edwin L Bridges
- Botanical and Ecological Consultant, Gig Harbor, Washington, USA
| | - Steve Orzell
- Avon Park Air Force Range, Avon Park, Florida, USA
| | | | - Colin E Hughes
- Department of Systematic and Evolutionary Botany, University of Zurich, Zurich, Switzerland
| |
Collapse
|
34
|
Allman ES, Baños H, Garrote-Lopez M, Rhodes JA. IDENTIFIABILITY OF LEVEL-1 SPECIES NETWORKS FROM GENE TREE QUARTETS. ARXIV 2024:arXiv:2401.06290v1. [PMID: 38259350 PMCID: PMC10802673] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Subscribe] [Scholar Register] [Indexed: 01/24/2024]
Abstract
When hybridization or other forms of lateral gene transfer have occurred, evolutionary relationships of species are better represented by phylogenetic networks than by trees. While inference of such networks remains challenging, several recently proposed methods are based on quartet concordance factors - the probabilities that a tree relating a gene sampled from the species displays the possible 4-taxon relationships. Building on earlier results, we investigate what level-1 network features are identifiable from concordance factors under the network multispecies coalescent model. We obtain results on both topological features of the network, and numerical parameters, uncovering a number of failures of identifiability related to 3-cycles in the network.
Collapse
|
35
|
Thom G, Moreira LR, Batista R, Gehara M, Aleixo A, Smith BT. Genomic Architecture Predicts Tree Topology, Population Structuring, and Demographic History in Amazonian Birds. Genome Biol Evol 2024; 16:evae002. [PMID: 38236173 PMCID: PMC10823491 DOI: 10.1093/gbe/evae002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 10/26/2023] [Accepted: 12/12/2023] [Indexed: 01/19/2024] Open
Abstract
Geographic barriers are frequently invoked to explain genetic structuring across the landscape. However, inferences on the spatial and temporal origins of population variation have been largely limited to evolutionary neutral models, ignoring the potential role of natural selection and intrinsic genomic processes known as genomic architecture in producing heterogeneity in differentiation across the genome. To test how variation in genomic characteristics (e.g. recombination rate) impacts our ability to reconstruct general patterns of differentiation between species that cooccur across geographic barriers, we sequenced the whole genomes of multiple bird populations that are distributed across rivers in southeastern Amazonia. We found that phylogenetic relationships within species and demographic parameters varied across the genome in predictable ways. Genetic diversity was positively associated with recombination rate and negatively associated with species tree support. Gene flow was less pervasive in genomic regions of low recombination, making these windows more likely to retain patterns of population structuring that matched the species tree. We further found that approximately a third of the genome showed evidence of selective sweeps and linked selection, skewing genome-wide estimates of effective population sizes and gene flow between populations toward lower values. In sum, we showed that the effects of intrinsic genomic characteristics and selection can be disentangled from neutral processes to elucidate spatial patterns of population differentiation.
Collapse
Affiliation(s)
- Gregory Thom
- Department of Ornithology, American Museum of Natural History, New York, NY, USA
- Museum of Natural Science, Louisiana State University, Baton Rouge, LA, USA
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, USA
| | - Lucas Rocha Moreira
- Program in Bioinformatics and Integrative Biology, University of Massachusetts Chan Medical School, Worcester, MA, USA
- Department of Vertebrate Genomics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Romina Batista
- Programa de Coleções Biológicas, Instituto Nacional de Pesquisas da Amazônia, Manaus, Brazil
- School of Science, Engineering and Environment, University of Salford, Manchester, UK
| | - Marcelo Gehara
- Department of Earth and Environmental Sciences, Rutgers University, Newark, NJ, USA
| | - Alexandre Aleixo
- Finnish Museum of Natural History, University of Helsinki, Helsinki, Finland
- Department of Environmental Genomics, Instituto Tecnológico Vale, Belém, Brazil
| | - Brian Tilston Smith
- Department of Ornithology, American Museum of Natural History, New York, NY, USA
| |
Collapse
|
36
|
Rodríguez-Machado S, Elías DJ, McMahan CD, Gruszkiewicz-Tolli A, Piller KR, Chakrabarty P. Disentangling historical relationships within Poeciliidae (Teleostei: Cyprinodontiformes) using ultraconserved elements. Mol Phylogenet Evol 2024; 190:107965. [PMID: 37977500 DOI: 10.1016/j.ympev.2023.107965] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2023] [Revised: 10/18/2023] [Accepted: 11/12/2023] [Indexed: 11/19/2023]
Abstract
Poeciliids (Cyprinodontiformes: Poeciliidae), commonly known as livebearers, are popular fishes in the aquarium trade (e.g., guppies, mollies, swordtails) that are widely distributed in the Americas, with 274 valid species in 27 genera. This group has undergone various taxonomic changes recently, spurred by investigations using traditional genetic markers. Here we used over 1,000 ultraconserved loci to infer the relationships within Poeciliidae in the first attempt at understanding their diversification based on genome-scale data. We explore gene tree discordance and investigate potential incongruence between concatenation and coalescent inference methods. Our aim is to examine the influence of incomplete lineage sorting and reticulate evolution on the poeciliids' evolutionary history and how these factors contribute to the observed gene tree discordace. Our concatenated and coalescent phylogenomic inferences recovered four major clades within Poeciliidae. Most supra-generic level relationships we inferred were congruent with previous molecular studies, but we found some disagreements; the Middle American taxa Phallichthys and Poecilia (Mollienesia) were recovered as non-monophyletic, and unlike other recent molecular studies, we recovered Brachyrhaphis as monophyletic. Our study is the first to provide signatures of reticulate evolution in Poeciliidae at the family level; however, continued finer-scale investigations are needed to understand the complex evolutionary history of the family along with a much-needed taxonomic re-evaluation.
Collapse
Affiliation(s)
- Sheila Rodríguez-Machado
- Museum of Natural Science, Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, United States.
| | - Diego J Elías
- Museum of Natural Science, Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, United States; Field Museum of Natural History, Chicago, IL 60605, United States
| | - Caleb D McMahan
- Field Museum of Natural History, Chicago, IL 60605, United States
| | - Anna Gruszkiewicz-Tolli
- Department of Biological Sciences, Southeastern Louisiana University, Hammond, LA 70402, United States
| | - Kyle R Piller
- Department of Biological Sciences, Southeastern Louisiana University, Hammond, LA 70402, United States
| | - Prosanta Chakrabarty
- Museum of Natural Science, Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, United States
| |
Collapse
|
37
|
Patané JSL, Martins J, Setubal JC. A Guide to Phylogenomic Inference. Methods Mol Biol 2024; 2802:267-345. [PMID: 38819564 DOI: 10.1007/978-1-0716-3838-5_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.
Collapse
Affiliation(s)
- José S L Patané
- Laboratório de Genética e Cardiologia Molecular, Instituto do Coração/Heart Institute Hospital das Clínicas - Faculdade de Medicina da Universidade de São Paulo São Paulo, São Paulo, SP, Brazil
| | - Joaquim Martins
- Integrative Omics group, Biorenewables National Laboratory, Brazilian Center for Research in Energy and Materials, Campinas, SP, Brazil
| | - João Carlos Setubal
- Departmento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo, SP, Brazil.
| |
Collapse
|
38
|
Qin YQ, Zhang MH, Yang CY, Nie ZL, Wen J, Meng Y. Phylogenomics and divergence pattern of Polygonatum (Asparagaceae: Polygonateae) in the north temperate region. Mol Phylogenet Evol 2024; 190:107962. [PMID: 37926394 DOI: 10.1016/j.ympev.2023.107962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Revised: 10/23/2023] [Accepted: 11/01/2023] [Indexed: 11/07/2023]
Abstract
Polygonatum is the largest genus of tribe Polygonateae (Asparagaceae) and is widely distributed in the temperate Northern Hemisphere, especially well diversified in southwestern China to northeastern Asia. Phylogenetic relationships of many species are still controversial. Hence it is necessary to clarify their phylogenetic relationships and infer possible reticulate relationships for the genus. In this study, genome-wide data of 43 species from Polygonatum and its closely related taxa were obtained by Hyb-Seq sequencing. The phylogenetic trees constructed from genome-wide nuclear and chloroplast sequences strongly supported the monophyly of Polygonatum with division into three major clades. A high level of incongruence was detected between nuclear and chloroplast trees as well as among gene trees within the genus, but all occurred within each major clade. However, introgression tests and reticulate evolution analyses revealed low level of gene flow and weak introgression events in the genus, suggesting hybridization and introgression were not dominant during the evolutionary diversification of Polygonatum in the Northern Hemisphere. This study provides important insights into reconstructing evolutionary relationships and speciation pattern of taxa from the north temperate flora.
Collapse
Affiliation(s)
- Yu-Qian Qin
- College of Biology and Environmental Sciences, Jishou University, Jishou, Hunan 416000, China
| | - Meng-Hua Zhang
- College of Biology and Environmental Sciences, Jishou University, Jishou, Hunan 416000, China
| | - Chu-Yun Yang
- College of Biology and Environmental Sciences, Jishou University, Jishou, Hunan 416000, China
| | - Ze-Long Nie
- College of Biology and Environmental Sciences, Jishou University, Jishou, Hunan 416000, China
| | - Jun Wen
- Department of Botany, National Museum of Natural History, Smithsonian Institution, Washington, DC 20013-7012, USA
| | - Ying Meng
- College of Biology and Environmental Sciences, Jishou University, Jishou, Hunan 416000, China.
| |
Collapse
|
39
|
Frankel LE, Ané C. Summary Tests of Introgression Are Highly Sensitive to Rate Variation Across Lineages. Syst Biol 2023; 72:1357-1369. [PMID: 37698548 DOI: 10.1093/sysbio/syad056] [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: 01/24/2023] [Revised: 07/07/2023] [Accepted: 08/29/2023] [Indexed: 09/13/2023] Open
Abstract
The evolutionary implications and frequency of hybridization and introgression are increasingly being recognized across the tree of life. To detect hybridization from multi-locus and genome-wide sequence data, a popular class of methods are based on summary statistics from subsets of 3 or 4 taxa. However, these methods often carry the assumption of a constant substitution rate across lineages and genes, which is commonly violated in many groups. In this work, we quantify the effects of rate variation on the D test (also known as ABBA-BABA test), the D3 test, and HyDe. All 3 tests are used widely across a range of taxonomic groups, in part because they are very fast to compute. We consider rate variation across species lineages, across genes, their lineage-by-gene interaction, and rate variation across gene-tree edges. We simulated species networks according to a birth-death-hybridization process, so as to capture a range of realistic species phylogenies. For all 3 methods tested, we found a marked increase in the false discovery of reticulation (type-1 error rate) when there is rate variation across species lineages. The D3 test was the most sensitive, with around 80% type-1 error, such that D3 appears to more sensitive to a departure from the clock than to the presence of reticulation. For all 3 tests, the power to detect hybridization events decreased as the number of hybridization events increased, indicating that multiple hybridization events can obscure one another if they occur within a small subset of taxa. Our study highlights the need to consider rate variation when using site-based summary statistics, and points to the advantages of methods that do not require assumptions on evolutionary rates across lineages or across genes.
Collapse
Affiliation(s)
- Lauren E Frankel
- Department of Botany, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Cécile Ané
- Department of Botany, University of Wisconsin-Madison, Madison, WI 53706, USA
- Department of Statistics, University of Wisconsin-Madison, Madison, WI 53706, USA
| |
Collapse
|
40
|
Cai ZY, Niu ZY, Zhang YY, Tong YH, Vu TC, Goh WL, Sungkaew S, Teerawatananon A, Xia NH. Phylogenomic analyses reveal reticulate evolution between Neomicrocalamus and Temochloa (Poaceae: Bambusoideae). FRONTIERS IN PLANT SCIENCE 2023; 14:1274337. [PMID: 38111884 PMCID: PMC10726129 DOI: 10.3389/fpls.2023.1274337] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Accepted: 11/14/2023] [Indexed: 12/20/2023]
Abstract
Neomicrocalamus and Temochloa are closely related to bamboo genera. However, when considered with newly discovered and morphologically similar material from China and Vietnam, the phylogenetic relationship among these three groups was ambiguous in the analyses based on DNA regions. Here, as a means of investigating the relationships among the three bamboo groups and exploring potential sources of genomic conflicts, we present a phylogenomic examination based on the whole plastome, single-nucleotide polymorphism (SNP), and single-copy nuclear (SCN) gene datasets. Three different phylogenetic hypotheses were found. The inconsistency is attributed to the combination of incomplete lineage sorting and introgression. The origin of newly discovered bamboos is from introgressive hybridization between Temochloa liliana (which contributed 80.7% of the genome) and Neomicrocalamus prainii (19.3%), indicating that the newly discovered bamboos are closer to T. liliana in genetics. The more similar morphology and closer distribution elevation also imply a closer relationship between Temochloa and newly discovered bamboos.
Collapse
Affiliation(s)
- Zhuo-Yu Cai
- Key Laboratory of Plant Resources, Conservation and Sustainable Utilization/Guangdong Provincial Key Laboratory of Digital Botanical Garden, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
- South China National Botanical Garden, Guangzhou, China
| | - Zheng-Yang Niu
- Key Laboratory of Plant Resources, Conservation and Sustainable Utilization/Guangdong Provincial Key Laboratory of Digital Botanical Garden, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
- South China National Botanical Garden, Guangzhou, China
| | - You-Yuan Zhang
- Key Laboratory of Plant Resources, Conservation and Sustainable Utilization/Guangdong Provincial Key Laboratory of Digital Botanical Garden, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
- Quality Management Office, Guiyang Vocational and Technical College, Guiyang, China
| | - Yi-Hua Tong
- Key Laboratory of Plant Resources, Conservation and Sustainable Utilization/Guangdong Provincial Key Laboratory of Digital Botanical Garden, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
- South China National Botanical Garden, Guangzhou, China
- State Key Laboratory of Plant Diversity and Specialty Crops, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
| | - Tien Chinh Vu
- Vietnam National Museum of Nature, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Wei Lim Goh
- Faculty of Science, Universiti Tunku Abdul Rahman, Jalan Universiti, Kampar, Perak, Malaysia
| | - Sarawood Sungkaew
- Department of Forest Biology, Faculty of Forestry, Kasetsart University, Bangkok, Thailand
| | | | - Nian-He Xia
- Key Laboratory of Plant Resources, Conservation and Sustainable Utilization/Guangdong Provincial Key Laboratory of Digital Botanical Garden, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
- South China National Botanical Garden, Guangzhou, China
| |
Collapse
|
41
|
Jin ZT, Hodel RGJ, Ma DK, Wang H, Liu GN, Ren C, Ge BJ, Fan Q, Jin SH, Xu C, Wu J, Liu BB. Nightmare or delight: Taxonomic circumscription meets reticulate evolution in the phylogenomic era. Mol Phylogenet Evol 2023; 189:107914. [PMID: 37666378 DOI: 10.1016/j.ympev.2023.107914] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2023] [Revised: 08/25/2023] [Accepted: 08/29/2023] [Indexed: 09/06/2023]
Abstract
Phylogenetic studies in the phylogenomics era have demonstrated that reticulate evolution greatly impedes the accuracy of phylogenetic inference, and consequently can obscure taxonomic treatments. However, the systematics community lacks a broadly applicable strategy for taxonomic delimitation in groups characterized by pervasive reticulate evolution. The red-fruit genus, Stranvaesia, provides an ideal model to examine the influence of reticulation on generic circumscription, particularly where hybridization and allopolyploidy dominate the evolutionary history. In this study, we conducted phylogenomic analyses integrating data from hundreds of single-copy nuclear (SCN) genes and plastomes, and interrogated nuclear paralogs to clarify the inter/intra-generic relationship of Stranvaesia and its allies in the framework of Maleae. Analyses of phylogenomic discord and phylogenetic networks showed that allopolyploidization and introgression promoted the origin and diversification of the Stranvaesia clade, a conclusion further bolstered by cytonuclear and gene tree discordance. With a well-inferred phylogenetic backbone, we propose an updated generic delimitation of Stranvaesia and introduce a new genus, Weniomeles. This new genus is distinguished by its purple-black fruits, thorns trunk and/or branches, and a distinctive fruit core anatomy characterized by multilocular separated by a layer of sclereids and a cluster of sclereids at the top of the locules. Through this study, we highlight a broadly-applicable workflow that underscores the significance of reticulate evolution analyses in shaping taxonomic revisions from phylogenomic data.
Collapse
Affiliation(s)
- Ze-Tao Jin
- College of Horticulture, State Key Laboratory of Crop Genetics & Germplasm Enhancement and Utilization, Nanjing Agricultural University, Nanjing, Jiangsu 210095, China; State Key Laboratory of Plant Diversity and Specialty Crops / State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; China National Botanical Garden, Beijing 100093, China
| | - Richard G J Hodel
- Department of Botany, National Museum of Natural History, Smithsonian Institution, Washington, DC 20013-7012, USA
| | - Dai-Kun Ma
- State Key Laboratory of Plant Diversity and Specialty Crops / State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; China National Botanical Garden, Beijing 100093, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hui Wang
- State Key Laboratory of Plant Diversity and Specialty Crops / State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; China National Botanical Garden, Beijing 100093, China; College of Forestry and Biotechnology, Zhejiang Agriculture and Forestry University, Hangzhou, Zhejiang 311300, China
| | | | - Chen Ren
- Key Laboratory of Plant Resources Conservation and Sustainable Utilization & Guangdong Provincial Key Laboratory of Applied Botany, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, Guangdong 510650, China
| | - Bin-Jie Ge
- Eastern China Conservation Center for Wild Endangered Plant Resources, Shanghai Chenshan Botanical Garden, Shanghai 201602, China
| | - Qiang Fan
- State Key Laboratory of Biocontrol and Guangdong Provincial Key Laboratory of Plant Resources, School of Life Sciences, Sun Yat-sen University, Guangzhou, Guangdong 510275, China
| | - Shui-Hu Jin
- College of Forestry and Biotechnology, Zhejiang Agriculture and Forestry University, Hangzhou, Zhejiang 311300, China
| | - Chao Xu
- State Key Laboratory of Plant Diversity and Specialty Crops / State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; China National Botanical Garden, Beijing 100093, China
| | - Jun Wu
- College of Horticulture, State Key Laboratory of Crop Genetics & Germplasm Enhancement and Utilization, Nanjing Agricultural University, Nanjing, Jiangsu 210095, China
| | - Bin-Bin Liu
- State Key Laboratory of Plant Diversity and Specialty Crops / State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China; China National Botanical Garden, Beijing 100093, China.
| |
Collapse
|
42
|
Jensen A, Swift F, de Vries D, Beck RMD, Kuderna LFK, Knauf S, Chuma IS, Keyyu JD, Kitchener AC, Farh K, Rogers J, Marques-Bonet T, Detwiler KM, Roos C, Guschanski K. Complex Evolutionary History With Extensive Ancestral Gene Flow in an African Primate Radiation. Mol Biol Evol 2023; 40:msad247. [PMID: 37987553 PMCID: PMC10691879 DOI: 10.1093/molbev/msad247] [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: 07/05/2023] [Revised: 10/17/2023] [Accepted: 11/09/2023] [Indexed: 11/22/2023] Open
Abstract
Understanding the drivers of speciation is fundamental in evolutionary biology, and recent studies highlight hybridization as an important evolutionary force. Using whole-genome sequencing data from 22 species of guenons (tribe Cercopithecini), one of the world's largest primate radiations, we show that rampant gene flow characterizes their evolutionary history and identify ancient hybridization across deeply divergent lineages that differ in ecology, morphology, and karyotypes. Some hybridization events resulted in mitochondrial introgression between distant lineages, likely facilitated by cointrogression of coadapted nuclear variants. Although the genomic landscapes of introgression were largely lineage specific, we found that genes with immune functions were overrepresented in introgressing regions, in line with adaptive introgression, whereas genes involved in pigmentation and morphology may contribute to reproductive isolation. In line with reports from other systems that hybridization might facilitate diversification, we find that some of the most species-rich guenon clades are of admixed origin. This study provides important insights into the prevalence, role, and outcomes of ancestral hybridization in a large mammalian radiation.
Collapse
Affiliation(s)
- Axel Jensen
- Department of Ecology and Genetics, Animal Ecology, Uppsala University, Uppsala SE-75236, Sweden
| | - Frances Swift
- School of Biological Sciences, Institute of Ecology and Evolution, University of Edinburgh, Edinburgh, UK
| | - Dorien de Vries
- School of Science, Engineering & Environment, University of Salford, Salford M5 4WT, UK
| | - Robin M D Beck
- School of Science, Engineering & Environment, University of Salford, Salford M5 4WT, UK
| | - Lukas F K Kuderna
- Illumina Artificial Intelligence Laboratory, Illumina Inc., Foster City, CA 94404, USA
| | - Sascha Knauf
- Institute of International Animal Health/One Health, Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Greifswald – Insel Riems 17493, Germany
| | | | - Julius D Keyyu
- Tanzania Wildlife Research Institute (TAWIRI), Arusha, Tanzania
| | - Andrew C Kitchener
- Department of Natural Sciences, National Museums Scotland, Edinburgh EH1 1JF, UK
- School of Geosciences, University of Edinburgh, Edinburgh EH8 9XP, UK
| | - Kyle Farh
- Illumina Artificial Intelligence Laboratory, Illumina Inc., Foster City, CA 94404, USA
| | - Jeffrey Rogers
- Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Tomas Marques-Bonet
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, Barcelona 08003, Spain
- Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Barcelona, Spain
- Catalan Institution of Research and Advanced Studies (ICREA), Barcelona, Spain
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Barcelona 08028, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA) and Universitat Pompeu Fabra, Barcelona 08010, Spain
| | - Kate M Detwiler
- Department of Biological Sciences, Florida Atlantic University, Boca Raton, FL, USA
| | - Christian Roos
- Gene Bank of Primates and Primate Genetics Laboratory, German Primate Center, Leibniz Institute for Primate Research, Göttingen 37077, Germany
| | - Katerina Guschanski
- Department of Ecology and Genetics, Animal Ecology, Uppsala University, Uppsala SE-75236, Sweden
- School of Biological Sciences, Institute of Ecology and Evolution, University of Edinburgh, Edinburgh, UK
| |
Collapse
|
43
|
Lescroart J, Bonilla-Sánchez A, Napolitano C, Buitrago-Torres DL, Ramírez-Chaves HE, Pulido-Santacruz P, Murphy WJ, Svardal H, Eizirik E. Extensive Phylogenomic Discordance and the Complex Evolutionary History of the Neotropical Cat Genus Leopardus. Mol Biol Evol 2023; 40:msad255. [PMID: 37987559 PMCID: PMC10701098 DOI: 10.1093/molbev/msad255] [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: 08/31/2023] [Revised: 11/07/2023] [Accepted: 11/13/2023] [Indexed: 11/22/2023] Open
Abstract
Even in the genomics era, the phylogeny of Neotropical small felids comprised in the genus Leopardus remains contentious. We used whole-genome resequencing data to construct a time-calibrated consensus phylogeny of this group, quantify phylogenomic discordance, test for interspecies introgression, and assess patterns of genetic diversity and demographic history. We infer that the Leopardus radiation started in the Early Pliocene as an initial speciation burst, followed by another in its subgenus Oncifelis during the Early Pleistocene. Our findings challenge the long-held notion that ocelot (Leopardus pardalis) and margay (L. wiedii) are sister species and instead indicate that margay is most closely related to the enigmatic Andean cat (L. jacobita), whose whole-genome data are reported here for the first time. In addition, we found that the newly sampled Andean tiger cat (L. tigrinus pardinoides) population from Colombia associates closely with Central American tiger cats (L. tigrinus oncilla). Genealogical discordance was largely attributable to incomplete lineage sorting, yet was augmented by strong gene flow between ocelot and the ancestral branch of Oncifelis, as well as between Geoffroy's cat (L. geoffroyi) and southern tiger cat (L. guttulus). Contrasting demographic trajectories have led to disparate levels of current genomic diversity, with a nearly tenfold difference in heterozygosity between Andean cat and ocelot, spanning the entire range of variability found in extant felids. Our analyses improved our understanding of the speciation history and diversity patterns in this felid radiation, and highlight the benefits to phylogenomic inference of embracing the many heterogeneous signals scattered across the genome.
Collapse
Affiliation(s)
- Jonas Lescroart
- Department of Biology, University of Antwerp, Antwerp, Belgium
- School of Health and Life Sciences, Pontifical Catholic University of Rio Grande do Sul, Porto Alegre, Brazil
| | - Alejandra Bonilla-Sánchez
- School of Health and Life Sciences, Pontifical Catholic University of Rio Grande do Sul, Porto Alegre, Brazil
- Faculty of Exact and Natural Sciences, University of Antioquia, Medellín, Colombia
| | - Constanza Napolitano
- Department of Biological Sciences and Biodiversity, University of Los Lagos, Osorno, Chile
- Institute of Ecology and Biodiversity, Concepción, Chile
- Cape Horn International Center, Puerto Williams, Chile
- Andean Cat Alliance, Villa Carlos Paz, Argentina
| | - Diana L Buitrago-Torres
- School of Health and Life Sciences, Pontifical Catholic University of Rio Grande do Sul, Porto Alegre, Brazil
| | - Héctor E Ramírez-Chaves
- Department of Biological Sciences, University of Caldas, Manizales, Colombia
- Centro de Museos, Museo de Historia Natural, University of Caldas, Manizales, Colombia
| | | | - William J Murphy
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX, USA
- Interdisciplinary Program in Genetics & Genomics, Texas A&M University, College Station, TX, USA
| | - Hannes Svardal
- Department of Biology, University of Antwerp, Antwerp, Belgium
- Naturalis Biodiversity Center, Leiden, Netherlands
| | - Eduardo Eizirik
- School of Health and Life Sciences, Pontifical Catholic University of Rio Grande do Sul, Porto Alegre, Brazil
- Instituto Pró-Carnívoros, Atibaia, Brazil
| |
Collapse
|
44
|
Cornwell A, Zhang Y, Thondamal M, Johnson DW, Thakar J, Samuelson AV. The C. elegans Myc-family of transcription factors coordinate a dynamic adaptive response to dietary restriction. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.22.568222. [PMID: 38045350 PMCID: PMC10690244 DOI: 10.1101/2023.11.22.568222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/05/2023]
Abstract
Dietary restriction (DR), the process of decreasing overall food consumption over an extended period of time, has been shown to increase longevity across evolutionarily diverse species and delay the onset of age-associated diseases in humans. In Caenorhabditis elegans, the Myc-family transcription factors (TFs) MXL-2 (Mlx) and MML-1 (MondoA/ChREBP), which function as obligate heterodimers, and PHA-4 (orthologous to forkhead box transcription factor A) are both necessary for the full physiological benefits of DR. However, the adaptive transcriptional response to DR and the role of MML-1::MXL-2 and PHA-4 remains elusive. We identified the transcriptional signature of C. elegans DR, using the eat-2 genetic model, and demonstrate broad changes in metabolic gene expression in eat-2 DR animals, which requires both mxl-2 and pha-4. While the requirement for these factors in DR gene expression overlaps, we found many of the DR genes exhibit an opposing change in relative gene expression in eat-2;mxl-2 animals compared to wild-type, which was not observed in eat-2 animals with pha-4 loss. We further show functional deficiencies of the mxl-2 loss in DR outside of lifespan, as eat-2;mxl-2 animals exhibit substantially smaller brood sizes and lay a proportion of dead eggs, indicating that MML-1::MXL-2 has a role in maintaining the balance between resource allocation to the soma and to reproduction under conditions of chronic food scarcity. While eat-2 animals do not show a significantly different metabolic rate compared to wild-type, we also find that loss of mxl-2 in DR does not affect the rate of oxygen consumption in young animals. The gene expression signature of eat-2 mutant animals is consistent with optimization of energy utilization and resource allocation, rather than induction of canonical gene expression changes associated with acute metabolic stress -such as induction of autophagy after TORC1 inhibition. Consistently, eat-2 animals are not substantially resistant to stress, providing further support to the idea that chronic DR may benefit healthspan and lifespan through efficient use of limited resources rather than broad upregulation of stress responses, and also indicates that MML-1::MXL-2 and PHA-4 may have different roles in promotion of benefits in response to different pro-longevity stimuli.
Collapse
Affiliation(s)
- Adam Cornwell
- Department of Biomedical Genetics, University of Rochester Medical Center, 601 Elmwood Avenue, Rochester, NY 14642, USA
| | - Yun Zhang
- Department of Biomedical Genetics, University of Rochester Medical Center, 601 Elmwood Avenue, Rochester, NY 14642, USA
| | - Manjunatha Thondamal
- Department of Biomedical Genetics, University of Rochester Medical Center, 601 Elmwood Avenue, Rochester, NY 14642, USA
- Department of Biological Sciences, GITAM University, Andhra Pradesh, India
| | - David W Johnson
- Department of Biomedical Genetics, University of Rochester Medical Center, 601 Elmwood Avenue, Rochester, NY 14642, USA
- Department of Math and Science, Genesee Community College, One College Rd Batavia, NY 14020, USA
| | - Juilee Thakar
- Department of Biomedical Genetics, University of Rochester Medical Center, 601 Elmwood Avenue, Rochester, NY 14642, USA
- Department of Biostatistics and Computational Biology, University of Rochester Medical Center, 601 Elmwood Avenue, Rochester, NY 14642, USA
- Department of Microbiology and Immunology, University of Rochester Medical Center, 601 Elmwood Avenue, Rochester, NY 14642, USA
| | - Andrew V Samuelson
- Department of Biomedical Genetics, University of Rochester Medical Center, 601 Elmwood Avenue, Rochester, NY 14642, USA
| |
Collapse
|
45
|
Zhang Q, Folk RA, Mo ZQ, Ye H, Zhang ZY, Peng H, Zhao JL, Yang SX, Yu XQ. Phylotranscriptomic analyses reveal deep gene tree discordance in Camellia (Theaceae). Mol Phylogenet Evol 2023; 188:107912. [PMID: 37648181 DOI: 10.1016/j.ympev.2023.107912] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Revised: 08/09/2023] [Accepted: 08/27/2023] [Indexed: 09/01/2023]
Abstract
Gene tree discordance is a significant legacy of biological evolution. Multiple factors can result in incongruence among genes, such as introgression, incomplete lineage sorting (ILS), gene duplication or loss. Resolving the background of gene tree discordance is a critical way to uncover the process of species diversification. Camellia, the largest genus in Theaceae, has controversial taxonomy and systematics due in part to a complex evolutionary history. We used 60 transcriptomes of 55 species, which represented 15 sections of Camellia to investigate its phylogeny and the possible causes of gene tree discordance. We conducted gene tree discordance analysis based on 1,617 orthologous low-copy nuclear genes, primarily using coalescent species trees and polytomy tests to distinguish hard and soft conflict. A selective pressure analysis was also performed to assess the impact of selection on phylogenetic topology reconstruction. Our results detected different levels of gene tree discordance in the backbone of Camellia, and recovered rapid diversification as one of the possible causes of gene tree discordance. Furthermore, we confirmed that none of the currently proposed sections of Camellia was monophyletic. Comparisons among datasets partitioned under different selective pressure regimes showed that integrating all orthologous genes provided the best phylogenetic resolution of the species tree of Camellia. The findings of this study reveal rapid diversification as a major source of gene tree discordance in Camellia and will facilitate future investigation of reticulate relationships at the species level in this important plant genus.
Collapse
Affiliation(s)
- Qiong Zhang
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, China; College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Ryan A Folk
- Department of Biological Sciences, Mississippi State University, MS 39762, United States
| | - Zhi-Qiong Mo
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, China
| | - Hang Ye
- Guangxi Key Laboratory of Special Non-wood Forest Cultivation and Utilization, Guangxi Forestry Research Institute, Nanning 530002, Guangxi, China
| | - Zhao-Yuan Zhang
- Guangxi Key Laboratory of Special Non-wood Forest Cultivation and Utilization, Guangxi Forestry Research Institute, Nanning 530002, Guangxi, China
| | - Hua Peng
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, China
| | - Jian-Li Zhao
- Yunnan Key Laboratory of Plant Reproductive Adaptation and Evolutionary Ecology and Institute of Biodiversity, School of Ecology and Environmental Science, Yunnan University, Kunming 650091, China.
| | - Shi-Xiong Yang
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, China.
| | - Xiang-Qin Yu
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, China.
| |
Collapse
|
46
|
Skopalíková J, Leong-Škorničková J, Šída O, Newman M, Chumová Z, Zeisek V, Jarolímová V, Poulsen AD, Dantas-Queiroz MV, Fér T, Záveská E. Ancient hybridization in Curcuma (Zingiberaceae)-Accelerator or brake in lineage diversifications? THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2023; 116:773-785. [PMID: 37537754 DOI: 10.1111/tpj.16408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 07/13/2023] [Accepted: 07/19/2023] [Indexed: 08/05/2023]
Abstract
Hybridization is a widespread phenomenon in the evolution of plants and exploring its role is crucial to understanding diversification processes of many taxonomic groups. Recently, more attention is focused on the role of ancient hybridization that has repeatedly been shown as triggers of evolutionary radiation, although in some cases, it can prevent further diversification. The causes, frequency, and consequences of ancient hybridization remain to be explored. Here, we present an account of several events of ancient hybridization in turmeric, the economically important plant genus Curcuma (Zingiberaceae), which harbors about 130 known species. We analyzed 1094 targeted low-copy genes and plastomes obtained by next-generation sequencing of 37 species of Curcuma, representing the known genetic diversity and spanning the geographical distribution of the genus. Using phylogenetic network analysis, we show that the entire genus Curcuma as well as its most speciose lineage arose via introgression from the genus Pyrgophyllum and one of the extinct lineages, respectively. We also document a single event of ancient hybridization, with C. vamana as a product, that represents an evolutionary dead end. We further discuss distinct circumstances of those hybridization events that deal mainly with (in)congruence in chromosome counts of the parental lineages.
Collapse
Affiliation(s)
- Jana Skopalíková
- Department of Botany, Charles University, Prague, Czech Republic
- Czech Academy of Sciences, Institute of Botany, Průhonice, Czech Republic
| | - Jana Leong-Škorničková
- The Herbarium, Singapore Botanic Gardens, Singapore
- Department of Biological Sciences, National University of Singapore, Singapore
| | - Otakar Šída
- Department of Botany, National Museum in Prague, Prague, Czech Republic
| | - Mark Newman
- Royal Botanic Garden Edinburgh, Edinburgh, Scotland, UK
| | - Zuzana Chumová
- Czech Academy of Sciences, Institute of Botany, Průhonice, Czech Republic
| | - Vojtěch Zeisek
- Department of Botany, Charles University, Prague, Czech Republic
- Czech Academy of Sciences, Institute of Botany, Průhonice, Czech Republic
| | - Vlasta Jarolímová
- Czech Academy of Sciences, Institute of Botany, Průhonice, Czech Republic
| | | | | | - Tomáš Fér
- Department of Botany, Charles University, Prague, Czech Republic
| | - Eliška Záveská
- Czech Academy of Sciences, Institute of Botany, Průhonice, Czech Republic
| |
Collapse
|
47
|
Fogg J, Allman ES, Ané C. PhyloCoalSimulations: A Simulator for Network Multispecies Coalescent Models, Including a New Extension for the Inheritance of Gene Flow. Syst Biol 2023; 72:1171-1179. [PMID: 37254872 DOI: 10.1093/sysbio/syad030] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2023] [Revised: 05/03/2023] [Accepted: 05/15/2023] [Indexed: 06/01/2023] Open
Abstract
We consider the evolution of phylogenetic gene trees along phylogenetic species networks, according to the network multispecies coalescent process, and introduce a new network coalescent model with correlated inheritance of gene flow. This model generalizes two traditional versions of the network coalescent: with independent or common inheritance. At each reticulation, multiple lineages of a given locus are inherited from parental populations chosen at random, either independently across lineages or with positive correlation according to a Dirichlet process. This process may account for locus-specific probabilities of inheritance, for example. We implemented the simulation of gene trees under these network coalescent models in the Julia package PhyloCoalSimulations, which depends on PhyloNetworks and its powerful network manipulation tools. Input species phylogenies can be read in extended Newick format, either in numbers of generations or in coalescent units. Simulated gene trees can be written in Newick format, and in a way that preserves information about their embedding within the species network. This embedding can be used for downstream purposes, such as to simulate species-specific processes like rate variation across species, or for other scenarios as illustrated in this note. This package should be useful for simulation studies and simulation-based inference methods. The software is available open source with documentation and a tutorial at https://github.com/cecileane/PhyloCoalSimulations.jl.
Collapse
Affiliation(s)
- John Fogg
- Department of Statistics, University of Wisconsin - Madison, WI, 53706, USA
| | - Elizabeth S Allman
- Department of Mathematics and Statistics, University of Alaska - Fairbanks, AK, 99775, USA
| | - Cécile Ané
- Department of Statistics, University of Wisconsin - Madison, WI, 53706, USA
- Department of Botany, University of Wisconsin - Madison, WI, 53706, USA
| |
Collapse
|
48
|
San Jose M, Doorenweerd C, Geib S, Barr N, Dupuis JR, Leblanc L, Kauwe A, Morris KY, Rubinoff D. Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex. Mol Phylogenet Evol 2023; 188:107892. [PMID: 37524217 DOI: 10.1016/j.ympev.2023.107892] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Revised: 07/07/2023] [Accepted: 07/28/2023] [Indexed: 08/02/2023]
Abstract
As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.
Collapse
Affiliation(s)
- Michael San Jose
- University of Hawaii, College of Tropical Agriculture and Human Resources, Department of Plant and Environmental Protection Sciences, Entomology Section, 3050 Maile Way, Honolulu, HI, 96822-2231, USA.
| | - Camiel Doorenweerd
- University of Hawaii, College of Tropical Agriculture and Human Resources, Department of Plant and Environmental Protection Sciences, Entomology Section, 3050 Maile Way, Honolulu, HI, 96822-2231, USA
| | - Scott Geib
- Tropical Crop and Commodity Protection Research Unit, Daniel K Inouye U.S. Pacific Basin Agricultural Center, USDA Agricultural Research Services, Hilo, HI, USA
| | - Norman Barr
- United States Department of Agriculture, Animal and Plant Health Inspection Service, Plant Protection and Quarantine, Science & Technology, Insect Management and Molecular Diagnostics Laboratory, 22675 N. Moorefield Road, Edinburg, TX 78541, USA
| | - Julian R Dupuis
- University of Kentucky, Department of Entomology, S225 Ag Science Center North, 1100 South Limestone, Lexington, KY, 40546-0091, USA
| | - Luc Leblanc
- University of Idaho, Department of Entomology, Plant Pathology and Nematology, 875 Perimeter Drive, MS2329, Moscow, ID, 83844-2329, USA
| | - Angela Kauwe
- Tropical Crop and Commodity Protection Research Unit, Daniel K Inouye U.S. Pacific Basin Agricultural Center, USDA Agricultural Research Services, Hilo, HI, USA
| | - Kimberley Y Morris
- University of Hawaii, College of Tropical Agriculture and Human Resources, Department of Plant and Environmental Protection Sciences, Entomology Section, 3050 Maile Way, Honolulu, HI, 96822-2231, USA; Tropical Crop and Commodity Protection Research Unit, Daniel K Inouye U.S. Pacific Basin Agricultural Center, USDA Agricultural Research Services, Hilo, HI, USA
| | - Daniel Rubinoff
- University of Hawaii, College of Tropical Agriculture and Human Resources, Department of Plant and Environmental Protection Sciences, Entomology Section, 3050 Maile Way, Honolulu, HI, 96822-2231, USA
| |
Collapse
|
49
|
Flouri T, Jiao X, Huang J, Rannala B, Yang Z. Efficient Bayesian inference under the multispecies coalescent with migration. Proc Natl Acad Sci U S A 2023; 120:e2310708120. [PMID: 37871206 PMCID: PMC10622872 DOI: 10.1073/pnas.2310708120] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2023] [Accepted: 08/15/2023] [Indexed: 10/25/2023] Open
Abstract
Analyses of genome sequence data have revealed pervasive interspecific gene flow and enriched our understanding of the role of gene flow in speciation and adaptation. Inference of gene flow using genomic data requires powerful statistical methods. Yet current likelihood-based methods involve heavy computation and are feasible for small datasets only. Here, we implement the multispecies-coalescent-with-migration model in the Bayesian program bpp, which can be used to test for gene flow and estimate migration rates, as well as species divergence times and population sizes. We develop Markov chain Monte Carlo algorithms for efficient sampling from the posterior, enabling the analysis of genome-scale datasets with thousands of loci. Implementation of both introgression and migration models in the same program allows us to test whether gene flow occurred continuously over time or in pulses. Analyses of genomic data from Anopheles mosquitoes demonstrate rich information in typical genomic datasets about the mode and rate of gene flow.
Collapse
Affiliation(s)
- Tomáš Flouri
- Department of Genetics, Evolution, and Environment, University College London, LondonWC1E 6BT, United Kingdom
| | - Xiyun Jiao
- Department of Statistics and Data Science, China Southern University of Science and Technology, Shenzhen518055, China
| | - Jun Huang
- Department of Intelligent Medical Engineering, School of Biomedical Engineering, Capital Medical University, Beijing100069, China
| | - Bruce Rannala
- Department of Evolution and Ecology, University of California, Davis, CA95616
| | - Ziheng Yang
- Department of Genetics, Evolution, and Environment, University College London, LondonWC1E 6BT, United Kingdom
| |
Collapse
|
50
|
Zhang C, Meng R, Meng Y, Guo BL, Liu QR, Nie ZL. Parallel evolution, atavism, and extensive introgression explain the radiation of Epimedium sect. Diphyllon (Berberidaceae) in southern East Asia. FRONTIERS IN PLANT SCIENCE 2023; 14:1234148. [PMID: 37915504 PMCID: PMC10616310 DOI: 10.3389/fpls.2023.1234148] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/03/2023] [Accepted: 09/29/2023] [Indexed: 11/03/2023]
Abstract
East Asia is the richest region of plant biodiversity in the northern temperate zone, and its radiation provides key insights for understanding rapid speciation, including evolutionary patterns and processes. However, it is challenging to investigate the recent evolutionary radiation among plants because of the lack of genetic divergence, phenotypic convergence, and interspecific gene flow. Epimedium sect. Diphyllon is a rarely studied plant lineage endemic to East Asia, especially highly diversified in its southern part. In this study, we report a robust phylogenomic analysis based on genotyping-by-sequencing data of this lineage. The results revealed a clear biogeographic pattern for Epimedium sect. Diphyllon with recognition into two major clades corresponding to the Sino-Himalayan and Sino-Japanese subkingdoms of East Asian Flora and rapid diversification of the extant species dated to the Pleistocene. Evolutionary radiation of Epimedium sect. Diphyllon is characterized by recent and predominant parallel evolution and atavism between the two subkingdom regions, with extensive reticulating hybridization within each region during the course of diversification in southern East Asia. A parallel-atavism-introgression hypothesis is referred to in explaining the radiation of plant diversity in southern East Asia, which represents a potential model for the rapid diversification of plants under global climate cooling in the late Tertiary. Our study advances our understanding of the evolutionary processes of plant radiation in East Asia as well as in other biodiversity hotspot regions.
Collapse
Affiliation(s)
- Cheng Zhang
- Key Laboratory of Biodiversity Science and Ecological Engineering of Ministry of Education, College of Life Sciences, Beijing Normal University, Beijing, China
- College of Biological Resources and Environmental Sciences, Jishou University, Jishou, Hunan, China
| | - Ran Meng
- College of Biological Resources and Environmental Sciences, Jishou University, Jishou, Hunan, China
| | - Ying Meng
- College of Biological Resources and Environmental Sciences, Jishou University, Jishou, Hunan, China
| | - Bao-Lin Guo
- Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, China
| | - Quan-Ru Liu
- Key Laboratory of Biodiversity Science and Ecological Engineering of Ministry of Education, College of Life Sciences, Beijing Normal University, Beijing, China
| | - Ze-Long Nie
- College of Biological Resources and Environmental Sciences, Jishou University, Jishou, Hunan, China
| |
Collapse
|