1
|
Rakotoarivelo AR, Rambuda T, Taron UH, Stalder G, O'Donoghue P, Robovský J, Hartmann S, Hofreiter M, Moodley Y. Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini). Mol Phylogenet Evol 2024; 198:108131. [PMID: 38909875 DOI: 10.1016/j.ympev.2024.108131] [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: 12/26/2023] [Revised: 05/19/2024] [Accepted: 06/15/2024] [Indexed: 06/25/2024]
Abstract
The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.
Collapse
Affiliation(s)
- Andrinajoro R Rakotoarivelo
- Department of Biological Sciences, University of Venda, Private Bag X5050, Thohoyandou 0950, Republic of South Africa; Department of Zoology and Entomology, University of the Free State: QwaQwa Campus, Private Bag X13, Phuthaditjhaba 9866, Republic of South Africa
| | - Thabelo Rambuda
- Department of Biological Sciences, University of Venda, Private Bag X5050, Thohoyandou 0950, Republic of South Africa; Department of Genetics, University of Pretoria, Private Bag X20, Hatfield 0028, Republic of South Africa
| | - Ulrike H Taron
- Evolutionary Adaptive Genomics, Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht-Strasse 24-25, 14476 Potsdam, Germany
| | - Gabrielle Stalder
- Department of Interdisciplinary Life Sciences, Research Institute of Wildlife Ecology, University of Veterinary Medicine, Savoyenstraße 1, A-1160 Wien, Austria
| | | | - Jan Robovský
- Department of Zoology, Faculty of Science, University of South Bohemia, Branišovská 1760, 37005 České Budějovice, Czech Republic
| | - Stefanie Hartmann
- Evolutionary Adaptive Genomics, Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht-Strasse 24-25, 14476 Potsdam, Germany
| | - Michael Hofreiter
- Evolutionary Adaptive Genomics, Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht-Strasse 24-25, 14476 Potsdam, Germany
| | - Yoshan Moodley
- Department of Biological Sciences, University of Venda, Private Bag X5050, Thohoyandou 0950, Republic of South Africa.
| |
Collapse
|
2
|
Ma XG, Ren YB, Sun H. Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China. Cladistics 2024; 40:357-373. [PMID: 38197450 DOI: 10.1111/cla.12570] [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: 04/17/2023] [Revised: 11/27/2023] [Accepted: 12/09/2023] [Indexed: 01/11/2024] Open
Abstract
Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.
Collapse
Affiliation(s)
- Xiang-Guang Ma
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, 650201, China
| | - Yue-Bo Ren
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, 650201, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Hang Sun
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, 650201, China
| |
Collapse
|
3
|
Adams K, Roux O. No sexual pheromones in Anopheles mosquitoes? CURRENT OPINION IN INSECT SCIENCE 2024; 64:101227. [PMID: 38936474 DOI: 10.1016/j.cois.2024.101227] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/24/2024] [Revised: 05/16/2024] [Accepted: 06/21/2024] [Indexed: 06/29/2024]
Abstract
Swarming behavior is the cornerstone of the anopheline mating system. At dusk, males congregate in monospecific swarms in which females come to find a mate once in their lives. Although many Anopheles species coexist in sympatry, hybrids are infrequent, suggesting the existence of strong premating reproductive barriers. Chemical cues, particularly pheromones, often play a crucial role in bringing sexes together in a species-specific manner among insects. While the existence of sexual pheromones in Anopheles species has been postulated, only a few studies developed experimental designs to investigate their presence. Here, we discuss the contrasting and debatable findings regarding both long-range and contact sex pheromones in the context of swarm ecology in Anopheles species.
Collapse
Affiliation(s)
- Kelsey Adams
- Immunology & Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, USA; Howard Hughes Medical Institute, Chevy Chase, USA
| | - Olivier Roux
- MIVEGEC, Univ. Montpellier, IRD, CNRS, Montpellier, France.
| |
Collapse
|
4
|
Zhang D, She H, Wang S, Wang H, Li S, Cheng Y, Song G, Jia C, Qu Y, Rheindt FE, Olsson U, Alström P, Lei F. Phylogenetic Conflict Between Species Tree and Maternally Inherited Gene Trees in a Clade of Emberiza Buntings (Aves: Emberizidae). Syst Biol 2024; 73:279-289. [PMID: 38157277 DOI: 10.1093/sysbio/syad078] [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/04/2022] [Revised: 12/06/2023] [Accepted: 12/27/2023] [Indexed: 01/03/2024] Open
Abstract
Different genomic regions may reflect conflicting phylogenetic topologies primarily due to incomplete lineage sorting and/or gene flow. Genomic data are necessary to reconstruct the true species tree and explore potential causes of phylogenetic conflict. Here, we investigate the phylogenetic relationships of 4 Emberiza species (Aves: Emberizidae) and discuss the potential causes of the observed mitochondrial non-monophyly of Emberiza godlewskii (Godlewski's bunting) using phylogenomic analyses based on whole genome resequencing data from 41 birds. Analyses based on both the whole mitochondrial genome and ~39 kilobases from the non-recombining W chromosome reveal sister relationships between each the northern and southern populations of E. godlewskii with E. cioides and E. cia, respectively. In contrast, the monophyly of E. godlewskii is reflected by the phylogenetic signal of autosomal and Z chromosomal sequence data as well as demographic inference analyses, which-in combination-support the following tree topology: ([{E. godlewskii, E. cia}, E. cioides], E. jankowskii). Using D-statistics, we detected multiple gene flow events among different lineages, indicating pervasive introgressive hybridization within this clade. Introgression from an unsampled lineage that is sister to E. cioides or introgression from an unsampled mitochondrial + W chromosomal lineage of E. cioides into northern E. godlewskii may explain the phylogenetic conflict between the species tree estimated from genome-wide data versus mtDNA/W tree topologies. These results underscore the importance of using genomic data for phylogenetic reconstruction and species delimitation.
Collapse
Affiliation(s)
- Dezhi Zhang
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
| | - Huishang She
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
| | - Shangyu Wang
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
- College of Life Sciences, University of Chinese Academy of Sciences, 100049 Beijing, China
| | - Haitao Wang
- School of Life Sciences, Jilin Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun 130024, China
| | - Shi Li
- College of Animal Science and Technology, Jilin Agricultural University, Changchun 130118, China
| | - Yalin Cheng
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
| | - Gang Song
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
| | - Chenxi Jia
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
| | - Yanhua Qu
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
| | - Frank E Rheindt
- Department of Biological Sciences, National University of Singapore, Singapore 117543, Republic of Singapore
| | - Urban Olsson
- Department of Biology and Environmental Science, University of Gothenburg, Box 463, SE-405 30 Gothenburg, Sweden
- Gothenburg Global Biodiversity Centre, Box 461, SE-405 30 Gothenburg, Sweden
| | - Per Alström
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
- Animal Ecology, Department of Ecology and Genetics, Evolutionary Biology Centre, Uppsala University, Norbyvägen 18 D, SE-752 36 Uppsala, Sweden
| | - Fumin Lei
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, 100101 Beijing, China
- College of Life Sciences, University of Chinese Academy of Sciences, 100049 Beijing, China
| |
Collapse
|
5
|
Leal JL, Milesi P, Hodková E, Zhou Q, James J, Eklund DM, Pyhäjärvi T, Salojärvi J, Lascoux M. Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles. Syst Biol 2024; 73:392-418. [PMID: 38613229 PMCID: PMC11282369 DOI: 10.1093/sysbio/syae012] [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/17/2023] [Revised: 12/18/2023] [Accepted: 03/28/2024] [Indexed: 04/14/2024] Open
Abstract
Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids, it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.
Collapse
Affiliation(s)
- J Luis Leal
- Plant Ecology and Evolution, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| | - Pascal Milesi
- Plant Ecology and Evolution, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
- Science for Life Laboratory (SciLifeLab), Uppsala University, 75237 Uppsala, Sweden
| | - Eva Hodková
- Plant Ecology and Evolution, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
- Faculty of Environmental Sciences, Czech University of Life Sciences Prague, Kamýcká 129, 16521 Prague, Czech Republic
| | - Qiujie Zhou
- Plant Ecology and Evolution, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| | - Jennifer James
- Plant Ecology and Evolution, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| | - D Magnus Eklund
- Physiology and Environmental Toxicology, Department of Organismal Biology, Uppsala University, Norbyvägen 18A, 75236 Uppsala, Sweden
| | - Tanja Pyhäjärvi
- Organismal and Evolutionary Biology Research Program, Faculty of Biological and Environmental Sciences, and Viikki Plant Science Centre, University of Helsinki, P.O. Box 65 (Viikinkaari 1), 00014 Helsinki, Finland
- Department of Forest Sciences, University of Helsinki, 00014 Helsinki, Finland
| | - Jarkko Salojärvi
- School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551, Singapore
- Organismal and Evolutionary Biology Research Program, Faculty of Biological and Environmental Sciences, and Viikki Plant Science Centre, University of Helsinki, P.O. Box 65 (Viikinkaari 1), 00014 Helsinki, Finland
| | - Martin Lascoux
- Plant Ecology and Evolution, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
- Science for Life Laboratory (SciLifeLab), Uppsala University, 75237 Uppsala, Sweden
| |
Collapse
|
6
|
Reyna-Blanco CS, Caduff M, Galimberti M, Leuenberger C, Wegmann D. Inference of Locus-Specific Population Mixtures from Linked Genome-Wide Allele Frequencies. Mol Biol Evol 2024; 41:msae137. [PMID: 38958167 PMCID: PMC11255385 DOI: 10.1093/molbev/msae137] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 06/26/2024] [Accepted: 06/27/2024] [Indexed: 07/04/2024] Open
Abstract
Admixture between populations and species is common in nature. Since the influx of new genetic material might be either facilitated or hindered by selection, variation in mixture proportions along the genome is expected in organisms undergoing recombination. Various graph-based models have been developed to better understand these evolutionary dynamics of population splits and mixtures. However, current models assume a single mixture rate for the entire genome and do not explicitly account for linkage. Here, we introduce TreeSwirl, a novel method for inferring branch lengths and locus-specific mixture proportions by using genome-wide allele frequency data, assuming that the admixture graph is known or has been inferred. TreeSwirl builds upon TreeMix that uses Gaussian processes to estimate the presence of gene flow between diverged populations. However, in contrast to TreeMix, our model infers locus-specific mixture proportions employing a hidden Markov model that accounts for linkage. Through simulated data, we demonstrate that TreeSwirl can accurately estimate locus-specific mixture proportions and handle complex demographic scenarios. It also outperforms related D- and f-statistics in terms of accuracy and sensitivity to detect introgressed loci.
Collapse
Affiliation(s)
- Carlos S Reyna-Blanco
- Department of Biology, University of Fribourg, Fribourg 1700, Switzerland
- Swiss Institute of Bioinformatics, Fribourg 1700, Switzerland
| | - Madleina Caduff
- Department of Biology, University of Fribourg, Fribourg 1700, Switzerland
- Swiss Institute of Bioinformatics, Fribourg 1700, Switzerland
| | - Marco Galimberti
- Department of Biology, University of Fribourg, Fribourg 1700, Switzerland
- Swiss Institute of Bioinformatics, Fribourg 1700, Switzerland
- Department of Psychiatry, Yale University School of Medicine, New Haven, CT, USA
- Veterans Affairs Connecticut Healthcare System, West Haven, CT, USA
| | | | - Daniel Wegmann
- Department of Biology, University of Fribourg, Fribourg 1700, Switzerland
- Swiss Institute of Bioinformatics, Fribourg 1700, Switzerland
| |
Collapse
|
7
|
Hancock PA, Ochomo E, Messenger LA. Genetic surveillance of insecticide resistance in African Anopheles populations to inform malaria vector control. Trends Parasitol 2024; 40:604-618. [PMID: 38760258 DOI: 10.1016/j.pt.2024.04.016] [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/28/2024] [Revised: 04/24/2024] [Accepted: 04/26/2024] [Indexed: 05/19/2024]
Abstract
Insecticide resistance in malaria vector populations poses a major threat to malaria control, which relies largely on insecticidal interventions. Contemporary vector-control strategies focus on combatting resistance using multiple insecticides with differing modes of action within the mosquito. However, diverse genetic resistance mechanisms are present in vector populations, and continue to evolve. Knowledge of the spatial distribution of these genetic mechanisms, and how they impact the efficacy of different insecticidal products, is critical to inform intervention deployment decisions. We developed a catalogue of genetic-resistance mechanisms in African malaria vectors that could guide molecular surveillance. We highlight situations where intervention deployment has led to resistance evolution and spread, and identify challenges in understanding and mitigating the epidemiological impacts of resistance.
Collapse
Affiliation(s)
- Penelope A Hancock
- Department of Infectious Disease Epidemiology, Imperial College London, London, UK.
| | - Eric Ochomo
- Centre for Global Health Research, Kenya Medical Research Institute, Kisumu, Kenya; Vector Group, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool, UK
| | - Louisa A Messenger
- Department of Environmental and Occupational Health, School of Public Health, University of Nevada, Las Vegas, USA; Parasitology and Vector Biology (PARAVEC) Laboratory, School of Public Health, University of Nevada, Las Vegas, USA
| |
Collapse
|
8
|
Manthey JD, Spellman GM. Recombination rate variation shapes genomic variability of phylogeographic structure in a widespread North American songbird (Aves: Certhia americana). Mol Phylogenet Evol 2024; 196:108088. [PMID: 38697377 DOI: 10.1016/j.ympev.2024.108088] [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: 01/29/2024] [Revised: 04/11/2024] [Accepted: 04/24/2024] [Indexed: 05/05/2024]
Abstract
The nonrandom distribution of chromosomal characteristics and functional elements-genomic architecture-impacts the relative strengths and impacts of population genetic processes across the genome. Due to this relationship, genomic architecture has the potential to shape variation in population genetic structure across the genome. Population genetic structure has been shown to vary across the genome in a variety of taxa, but this body of work has largely focused on pairwise population genomic comparisons between closely related taxa. Here, we used whole genome sequencing of seven phylogeographically structured populations of a North American songbird, the Brown Creeper (Certhia americana), to determine the impacts of genomic architecture on phylogeographic structure variation across the genome. Using multiple methods to infer phylogeographic structure-ordination, clustering, and phylogenetic methods-we found that recombination rate variation explained a large proportion of phylogeographic structure variation. Genomic regions with low recombination showed phylogeographic structure consistent with the genome-wide pattern. In regions with high recombination, we found strong phylogeographic structure, but with discordant patterns relative to the genome-wide pattern. In regions with high recombination rate, we found that populations with small effective population sizes evolve relatively more rapidly than larger populations, leading to discordant signatures of phylogeographic structure. These results suggest that the interplay between recombination rate variation and effective population sizes shape the relative impacts of selection and genetic drift in different parts of the genome. Overall, the combined interactions of population genetic processes, genomic architecture, and effective population sizes shape patterns of variability in phylogeographic structure across the genome of the Brown Creeper.
Collapse
Affiliation(s)
- Joseph D Manthey
- Department of Biological Sciences, Texas Tech University. Lubbock, TX, USA.
| | - Garth M Spellman
- Department of Zoology, Denver Museum of Nature & Science, Denver, CO, USA
| |
Collapse
|
9
|
Shi YJ, Mi JX, Huang JL, Tian FF, He F, Zhong Y, Yang HB, Wang F, Xiao Y, Yang LK, Zhang F, Chen LH, Wan XQ. A new species of Populus and the extensive hybrid speciation arising from it on the Qinghai-Tibet Plateau. Mol Phylogenet Evol 2024; 196:108072. [PMID: 38615706 DOI: 10.1016/j.ympev.2024.108072] [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/10/2023] [Revised: 03/25/2024] [Accepted: 04/11/2024] [Indexed: 04/16/2024]
Abstract
While the diversity of species formation is broadly acknowledged, significant debate exists regarding the universal nature of hybrid species formation. Through an 18-year comprehensive study of all Populus species on the Qinghai-Tibet Plateau, 23 previously recorded species and 8 new species were identified. Based on morphological characteristics, these can be classified into three groups: species in section Leucoides, species with large leaves, and species with small leaves in section Tacamahaca. By conducting whole-genome re-sequencing of 150 genotypes from these 31 species, 2.28 million single nucleotide polymorphisms (SNPs) were identified. Phylogenetic analysis utilizing these SNPs not only revealed a highly intricate evolutionary network within the large-leaf species of section Tacamahaca but also confirmed that a new species, P. curviserrata, naturally hybridized with P. cathayana, P. szechuanica, and P. ciliata, resulting in 11 hybrid species. These findings indicate the widespread occurrence of hybrid species formation within this genus, with hybridization serving as a key evolutionary mechanism for Populus on the plateau. A novel hypothesis, "Hybrid Species Exterminating Their Ancestral Species (HSEAS)," is introduced to explain the mechanisms of hybrid species formation at three different scales: the entire plateau, the southeastern mountain region, and individual river valleys.
Collapse
Affiliation(s)
- Yu-Jie Shi
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China; Zhejiang Provincial Key Laboratory of Plant Evolutionary Ecology and Conservation, Taizhou University, Taizhou 318000, China
| | - Jia-Xuan Mi
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China
| | - Jin-Liang Huang
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China
| | - Fei-Fei Tian
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China
| | - Fang He
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China
| | - Yu Zhong
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China; Sichuan Province Key Laboratory of Ecological Forestry Engineering On the Upper Reaches of the Yangtze River, China
| | - Han-Bo Yang
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China; Sichuan Province Key Laboratory of Ecological Forestry Engineering On the Upper Reaches of the Yangtze River, China
| | - Fang Wang
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China
| | - Yu Xiao
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China
| | - Lin-Kai Yang
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China
| | - Fan Zhang
- College of Landscape Architecture, Sichuan Agricultural University, Chengdu 611130, China
| | - Liang-Hua Chen
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China; Sichuan Province Key Laboratory of Ecological Forestry Engineering On the Upper Reaches of the Yangtze River, China
| | - Xue Qin Wan
- College of Forestry, Sichuan Agricultural University, Chengdu 611130, China; Sichuan Province Key Laboratory of Ecological Forestry Engineering On the Upper Reaches of the Yangtze River, China.
| |
Collapse
|
10
|
Mathieu-Bégné E, Kincaid-Smith J, Chaparro C, Allienne JF, Rey O, Boissier J, Toulza E. Schistosoma haematobium and Schistosoma bovis first generation hybrids undergo gene expressions changes consistent with species compatibility and heterosis. PLoS Negl Trop Dis 2024; 18:e0012267. [PMID: 38954732 PMCID: PMC11249247 DOI: 10.1371/journal.pntd.0012267] [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: 05/22/2023] [Revised: 07/15/2024] [Accepted: 06/03/2024] [Indexed: 07/04/2024] Open
Abstract
When two species hybridize, the two parental genomes are brought together and some alleles might interact for the first time. To date, the extent of the transcriptomic changes in first hybrid generations, along with their functional outcome constitute an important knowledge gap, especially in parasite species. Here we explored the molecular and functional outcomes of hybridization in first-generation hybrids between the blood fluke parasites Schistosoma haematobium and S. bovis. Through a transcriptomic approach, we measured gene expression in both parental species and hybrids. We described and quantified expression profiles encountered in hybrids along with the main biological processes impacted. Up to 7,100 genes fell into a particular hybrid expression profile (intermediate between the parental expression levels, over-expressed, under-expressed, or expressed like one of the parental lines). Most of these genes were different depending on the direction of the parental cross (S. bovis mother and S. haematobium father or the reverse) and depending on the sex. For a given sex and cross direction, the vast majority of genes were hence unassigned to a hybrid expression profile: either they were differentially expressed genes but not typical of any hybrid expression profiles or they were not differentially expressed neither between hybrids and parental lines nor between parental lines. The most prevalent profile of gene expression in hybrids was the intermediate one (24% of investigated genes). These results suggest that transcriptomic compatibility between S. haematobium and S. bovis remains quite high. We also found support for an over-dominance model (over- and under-expressed genes in hybrids compared to parental lines) potentially associated with heterosis. In females in particular, processes such as reproductive processes, metabolism and cell interactions as well as signaling pathways were indeed affected. Our study hence provides new insight on the biology of Schistosoma hybrids with evidences supporting compatibility and heterosis.
Collapse
Affiliation(s)
| | - Julien Kincaid-Smith
- IHPE, Université de Montpellier, CNRS, Ifremer, Université Perpignan Via Domitia, Perpignan, France
- CBGP, IRD, CIRAD, INRAE, Institut Agro, Université de Montpellier, Montpellier, France
| | - Cristian Chaparro
- IHPE, Université de Montpellier, CNRS, Ifremer, Université Perpignan Via Domitia, Perpignan, France
| | - Jean-François Allienne
- IHPE, Université de Montpellier, CNRS, Ifremer, Université Perpignan Via Domitia, Perpignan, France
| | - Olivier Rey
- IHPE, Université de Montpellier, CNRS, Ifremer, Université Perpignan Via Domitia, Perpignan, France
| | - Jérôme Boissier
- IHPE, Université de Montpellier, CNRS, Ifremer, Université Perpignan Via Domitia, Perpignan, France
| | - Eve Toulza
- IHPE, Université de Montpellier, CNRS, Ifremer, Université Perpignan Via Domitia, Perpignan, France
| |
Collapse
|
11
|
Knief U, Müller IA, Stryjewski KF, Metzler D, Sorenson MD, Wolf JBW. Evolution of Chromosomal Inversions across an Avian Radiation. Mol Biol Evol 2024; 41:msae092. [PMID: 38743589 PMCID: PMC11152452 DOI: 10.1093/molbev/msae092] [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/08/2024] [Revised: 04/05/2024] [Accepted: 05/03/2024] [Indexed: 05/16/2024] Open
Abstract
Chromosomal inversions are structural mutations that can play a prominent role in adaptation and speciation. Inversions segregating across species boundaries (trans-species inversions) are often taken as evidence for ancient balancing selection or adaptive introgression, but can also be due to incomplete lineage sorting. Using whole-genome resequencing data from 18 populations of 11 recognized munia species in the genus Lonchura (N = 176 individuals), we identify four large para- and pericentric inversions ranging in size from 4 to 20 Mb. All four inversions cosegregate across multiple species and predate the numerous speciation events associated with the rapid radiation of this clade across the prehistoric Sahul (Australia, New Guinea) and Bismarck Archipelago. Using coalescent theory, we infer that trans-specificity is improbable for neutrally segregating variation despite substantial incomplete lineage sorting characterizing this young radiation. Instead, the maintenance of all three autosomal inversions (chr1, chr5, and chr6) is best explained by selection acting along ecogeographic clines not observed for the collinear parts of the genome. In addition, the sex chromosome inversion largely aligns with species boundaries and shows signatures of repeated positive selection for both alleles. This study provides evidence for trans-species inversion polymorphisms involved in both adaptation and speciation. It further highlights the importance of informing selection inference using a null model of neutral evolution derived from the collinear part of the genome.
Collapse
Affiliation(s)
- Ulrich Knief
- Division of Evolutionary Biology, Faculty of Biology, LMU Munich, 82152 Planegg-Martinsried, Germany
- Evolutionary Biology & Ecology, Faculty of Biology, University of Freiburg, 79104 Freiburg, Germany
| | - Ingo A Müller
- Division of Evolutionary Biology, Faculty of Biology, LMU Munich, 82152 Planegg-Martinsried, Germany
- Department of Bioinformatics and Genetics, Swedish Museum of Natural History, 11418 Stockholm, Sweden
- Division of Systematics and Evolution, Department of Zoology, Stockholm University, 11418 Stockholm, Sweden
| | | | - Dirk Metzler
- Division of Evolutionary Biology, Faculty of Biology, LMU Munich, 82152 Planegg-Martinsried, Germany
| | | | - Jochen B W Wolf
- Division of Evolutionary Biology, Faculty of Biology, LMU Munich, 82152 Planegg-Martinsried, Germany
| |
Collapse
|
12
|
Bouafou L, Makanga BK, Rahola N, Boddé M, Ngangué MF, Daron J, Berger A, Mouillaud T, Makunin A, Korlević P, Nwezeobi J, Kengne P, Paupy C, Lawniczak MKN, Ayala D. Host preference patterns in domestic and wild settings: Insights into Anopheles feeding behavior. Evol Appl 2024; 17:e13693. [PMID: 38828055 PMCID: PMC11143308 DOI: 10.1111/eva.13693] [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/28/2023] [Revised: 03/28/2024] [Accepted: 04/02/2024] [Indexed: 06/05/2024] Open
Abstract
The adaptation of Anopheles malaria vectors to domestic settings is directly linked to their ability to feed on humans. The strength of this species-habitat association is unequal across the species within the genus, with the major vectors being particularly dependent on humans. However, our understanding of how blood-feeding behavior interacts with and adapts to environmental settings, including the presence of humans, remains limited. Using a field-based approach, we first investigated Anopheles community structure and feeding behavior patterns in domestic and sylvatic settings in La Lopé National Park in Gabon, Central Africa. We characterized the preference indices using a dual-host choice sampling approach across mosquito species, habitats, and seasons. We then quantified the plastic biting behavior of mosquito species in each habitat. We collected individuals from 16 Anopheles species that exhibited significant differences in species composition and abundance between sylvatic and domestic settings. The host-seeking behavior also varied among the seven most abundant species. The general attractiveness to each host, human or animal, remained relatively constant for each species, but with significant variations between habitats across species. These variations, to more generalist and to more anthropophilic behavior, were related to seasonal changes and distance from the village, respectively. Finally, we pointed out that the host choice of major malaria vectors changed in the absence of humans, revealing a plastic feeding behavior of these species. This study highlights the effect of humans on Anopheles distribution and feeding evolution. The characterization of feeding behavior in wild and domestic settings provides opportunities to better understand the interplay between genetic determinants of host preference and ecological factors. Our findings suggest that protected areas may offer alternative thriving conditions to major malaria vectors.
Collapse
Affiliation(s)
- Lemonde Bouafou
- UMR MIVEGEC, University of Montpellier, CNRS, IRDMontpellierFrance
- CIRMFFrancevilleGabon
| | | | - Nil Rahola
- UMR MIVEGEC, University of Montpellier, CNRS, IRDMontpellierFrance
| | | | | | - Josquin Daron
- UMR MIVEGEC, University of Montpellier, CNRS, IRDMontpellierFrance
| | - Audric Berger
- UMR MIVEGEC, University of Montpellier, CNRS, IRDMontpellierFrance
| | - Theo Mouillaud
- UMR MIVEGEC, University of Montpellier, CNRS, IRDMontpellierFrance
| | | | | | | | - Pierre Kengne
- UMR MIVEGEC, University of Montpellier, CNRS, IRDMontpellierFrance
- CIRMFFrancevilleGabon
| | - Christophe Paupy
- UMR MIVEGEC, University of Montpellier, CNRS, IRDMontpellierFrance
| | | | - Diego Ayala
- UMR MIVEGEC, University of Montpellier, CNRS, IRDMontpellierFrance
- Medical Entomology UnitInstitut Pasteur de MadagascarAntananarivoMadagascar
| |
Collapse
|
13
|
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
|
14
|
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
|
15
|
Musher LJ, Del-Rio G, Marcondes RS, Brumfield RT, Bravo GA, Thom G. Geogenomic Predictors of Genetree Heterogeneity Explain Phylogeographic and Introgression History: A Case Study in an Amazonian Bird (Thamnophilus aethiops). Syst Biol 2024; 73:36-52. [PMID: 37804132 DOI: 10.1093/sysbio/syad061] [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/20/2022] [Revised: 09/14/2023] [Accepted: 10/04/2023] [Indexed: 10/08/2023] Open
Abstract
Can knowledge about genome architecture inform biogeographic and phylogenetic inference? Selection, drift, recombination, and gene flow interact to produce a genomic landscape of divergence wherein patterns of differentiation and genealogy vary nonrandomly across the genomes of diverging populations. For instance, genealogical patterns that arise due to gene flow should be more likely to occur on smaller chromosomes, which experience high recombination, whereas those tracking histories of geographic isolation (reduced gene flow caused by a barrier) and divergence should be more likely to occur on larger and sex chromosomes. In Amazonia, populations of many bird species diverge and introgress across rivers, resulting in reticulated genomic signals. Herein, we used reduced representation genomic data to disentangle the evolutionary history of 4 populations of an Amazonian antbird, Thamnophilus aethiops, whose biogeographic history was associated with the dynamic evolution of the Madeira River Basin. Specifically, we evaluate whether a large river capture event ca. 200 Ka, gave rise to reticulated genealogies in the genome by making spatially explicit predictions about isolation and gene flow based on knowledge about genomic processes. We first estimated chromosome-level phylogenies and recovered 2 primary topologies across the genome. The first topology (T1) was most consistent with predictions about population divergence and was recovered for the Z-chromosome. The second (T2), was consistent with predictions about gene flow upon secondary contact. To evaluate support for these topologies, we trained a convolutional neural network to classify our data into alternative diversification models and estimate demographic parameters. The best-fit model was concordant with T1 and included gene flow between non-sister taxa. Finally, we modeled levels of divergence and introgression as functions of chromosome length and found that smaller chromosomes experienced higher gene flow. Given that (1) genetrees supporting T2 were more likely to occur on smaller chromosomes and (2) we found lower levels of introgression on larger chromosomes (and especially the Z-chromosome), we argue that T1 represents the history of population divergence across rivers and T2 the history of secondary contact due to barrier loss. Our results suggest that a significant portion of genomic heterogeneity arises due to extrinsic biogeographic processes such as river capture interacting with intrinsic processes associated with genome architecture. Future phylogeographic studies would benefit from accounting for genomic processes, as different parts of the genome reveal contrasting, albeit complementary histories, all of which are relevant for disentangling the intricate geogenomic mechanisms of biotic diversification. [Amazonia; biogeography; demographic modeling; gene flow; gene tree; genome architecture; geogenomics; introgression; linked selection; neural network; phylogenomic; phylogeography; reproductive isolation; speciation; species tree.].
Collapse
Affiliation(s)
- Lukas J Musher
- Department of Ornithology, The Academy of Natural Sciences of Drexel University, Philadelphia, PA 19103, USA
- Department of Ornithology, American Museum of Natural History, New York, NY 10024, USA
| | - Glaucia Del-Rio
- Cornell Laboratory of Ornithology and Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY 14853, USA
- Florida Museum of Natural History, University of Florida, Gainesville, FL 32611, USA
| | - Rafael S Marcondes
- Department of Biology and Museum of Natural Science, Louisiana State University, Baton Rouge, LA 70803, USA
- Department of BioSciences, Rice University, Houston, TX 77005, USA
| | - Robb T Brumfield
- Department of Biology and Museum of Natural Science, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Gustavo A Bravo
- Sección de Ornitología, Colecciones Biológicas, Instituto de Investigación de Recursos Biológicos Alexander von Humboldt, Claustro de San Agustín, Villa de Leyva, Boyacá 111311, Colombia
- Museum of Comparative Zoology and Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
| | - Gregory Thom
- Department of Biology and Museum of Natural Science, Louisiana State University, Baton Rouge, LA 70803, USA
| |
Collapse
|
16
|
Daron J, Bouafou L, Tennessen JA, Rahola N, Makanga B, Akone-Ella O, Ngangue MF, Longo Pendy NM, Paupy C, Neafsey DE, Fontaine MC, Ayala D. Genomic Signatures of Microgeographic Adaptation in Anopheles coluzzii Along an Anthropogenic Gradient in Gabon. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.16.594472. [PMID: 38798379 PMCID: PMC11118577 DOI: 10.1101/2024.05.16.594472] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
Species distributed across heterogeneous environments often evolve locally adapted populations, but understanding how these persist in the presence of homogenizing gene flow remains puzzling. In Gabon, Anopheles coluzzii, a major African malaria mosquito is found along an ecological gradient, including a sylvatic population, away of any human presence. This study identifies into the genomic signatures of local adaptation in populations from distinct environments including the urban area of Libreville, and two proximate sites 10km apart in the La Lopé National Park (LLP), a village and its sylvatic neighborhood. Whole genome re-sequencing of 96 mosquitoes unveiled ∼ 5.7millions high-quality single nucleotide polymorphisms. Coalescent-based demographic analyses suggest an ∼ 8,000-year-old divergence between Libreville and La Lopé populations, followed by a secondary contact ( ∼ 4,000 ybp) resulting in asymmetric effective gene flow. The urban population displayed reduced effective size, evidence of inbreeding, and strong selection pressures for adaptation to urban settings, as suggested by the hard selective sweeps associated with genes involved in detoxification and insecticide resistance. In contrast, the two geographically proximate LLP populations showed larger effective sizes, and distinctive genomic differences in selective signals, notably soft-selective sweeps on the standing genetic variation. Although neutral loci and chromosomal inversions failed to discriminate between LLP populations, our findings support that microgeographic adaptation can swiftly emerge through selection on standing genetic variation despite high gene flow. This study contributes to the growing understanding of evolution of populations in heterogeneous environments amid ongoing gene flow and how major malaria mosquitoes adapt to human. Significance Anopheles coluzzii , a major African malaria vector, thrives from humid rainforests to dry savannahs and coastal areas. This ecological success is linked to its close association with domestic settings, with human playing significant roles in driving the recent urban evolution of this mosquito. Our research explores the assumption that these mosquitoes are strictly dependent on human habitats, by conducting whole-genome sequencing on An. coluzzii specimens from urban, rural, and sylvatic sites in Gabon. We found that urban mosquitoes show de novo genetic signatures of human-driven vector control, while rural and sylvatic mosquitoes exhibit distinctive genetic evidence of local adaptations derived from standing genetic variation. Understanding adaptation mechanisms of this mosquito is therefore crucial to predict evolution of vector control strategies.
Collapse
|
17
|
van der Gulik PTS, Hoff WD, Speijer D. The contours of evolution: In defence of Darwin's tree of life paradigm. Bioessays 2024; 46:e2400012. [PMID: 38436469 DOI: 10.1002/bies.202400012] [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: 01/15/2024] [Revised: 02/12/2024] [Accepted: 02/15/2024] [Indexed: 03/05/2024]
Abstract
Both the concept of a Darwinian tree of life (TOL) and the possibility of its accurate reconstruction have been much criticized. Criticisms mostly revolve around the extensive occurrence of lateral gene transfer (LGT), instances of uptake of complete organisms to become organelles (with the associated subsequent gene transfer to the nucleus), as well as the implications of more subtle aspects of the biological species concept. Here we argue that none of these criticisms are sufficient to abandon the valuable TOL concept and the biological realities it captures. Especially important is the need to conceptually distinguish between organismal trees and gene trees, which necessitates incorporating insights of widely occurring LGT into modern evolutionary theory. We demonstrate that all criticisms, while based on important new findings, do not invalidate the TOL. After considering the implications of these new insights, we find that the contours of evolution are best represented by a TOL.
Collapse
Affiliation(s)
| | - Wouter D Hoff
- Department of Microbiology and Molecular Genetics, Oklahoma State University, Stillwater, Oklahoma, USA
| | - Dave Speijer
- Department of Medical Biochemistry, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
| |
Collapse
|
18
|
Kamau L, Bennett KL, Ochomo E, Herren J, Agumba S, Otieno S, Omoke D, Matoke-Muhia D, Mburu D, Mwangangi J, Ramaita E, Juma EO, Mbogo C, Barasa S, Miles A. The Anopheles coluzzii range extends into Kenya: detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa. Malar J 2024; 23:122. [PMID: 38671462 PMCID: PMC11046809 DOI: 10.1186/s12936-024-04950-x] [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/13/2024] [Accepted: 04/15/2024] [Indexed: 04/28/2024] Open
Abstract
BACKGROUND Anopheles coluzzii is a primary vector of malaria found in West and Central Africa, but its presence has hitherto never been documented in Kenya. A thorough understanding of vector bionomics is important as it enables the implementation of targeted and effective vector control interventions. Malaria vector surveillance efforts in the country have tended to focus on historically known primary vectors. The current study sought to determine the taxonomic status of samples collected from five different malaria epidemiological zones in Kenya as well as describe the population genetic structure and insecticide resistance profiles in relation to other An. coluzzii populations. METHODS Mosquitoes were sampled as larvae from Busia, Kwale, Turkana, Kirinyaga and Kiambu counties, representing the range of malaria endemicities in Kenya, in 2019 and 2021 and emergent adults analysed using Whole Genome Sequencing (WGS) data processed in accordance with the Anopheles gambiae 1000 Genomes Project phase 3. Where available, historical samples from the same sites were included for WGS. Comparisons were made with An. coluzzii cohorts from West and Central Africa. RESULTS This study reports the detection of An. coluzzii for the first time in Kenya. The species was detected in Turkana County across all three time points from which samples were analyzed and its presence confirmed through taxonomic analysis. Additionally, there was a lack of strong population genetic differentiation between An. coluzzii from Kenya and those from the more northerly regions of West and Central Africa, suggesting they represent a connected extension to the known species range. Mutations associated with target-site resistance to DDT and pyrethroids and metabolic resistance to DDT were found at high frequencies up to 64%. The profile and frequencies of the variants observed were similar to An. coluzzii from West and Central Africa but the ace-1 mutation linked to organophosphate and carbamate resistance present in An. coluzzii from coastal West Africa was absent in Kenya. CONCLUSIONS These findings emphasize the need for the incorporation of genomics in comprehensive and routine vector surveillance to inform on the range of malaria vector species, and their insecticide resistance status to inform the choice of effective vector control approaches.
Collapse
Affiliation(s)
- Luna Kamau
- Centre for Biotechnology Research and Development (CBRD), Kenya Medical Research Institute, PO Box 54840, Nairobi, 00200, Kenya.
| | - Kelly L Bennett
- Malaria Vector Genomic Surveillance, Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK
| | - Eric Ochomo
- Centre for Global Health Research (CGHR), Kenya Medical Research Institute, Nairobi, Kenya
- Liverpool School of Tropical Medicine, Liverpool, United Kingdom
| | - Jeremy Herren
- International Center for Insect Physiology and Ecology (Icipe), Nairobi, Kenya
| | - Silas Agumba
- Centre for Global Health Research (CGHR), Kenya Medical Research Institute, Nairobi, Kenya
| | - Samson Otieno
- Centre for Global Health Research (CGHR), Kenya Medical Research Institute, Nairobi, Kenya
| | - Diana Omoke
- Centre for Global Health Research (CGHR), Kenya Medical Research Institute, Nairobi, Kenya
| | - Damaris Matoke-Muhia
- Centre for Biotechnology Research and Development (CBRD), Kenya Medical Research Institute, PO Box 54840, Nairobi, 00200, Kenya
- Pan African Mosquito Control Association (PAMCA), Nairobi, Kenya
| | - David Mburu
- Pwani University Biosciences Research Centre (PUBReC), Kilifi, Kenya
| | - Joseph Mwangangi
- Centre for Geographic Medicine Research-Coast (CGMR-C), Kenya Medical Research Institute, Nairobi, Kenya
| | - Edith Ramaita
- Ministry of Health-National Malaria Control Programme (NMCP), Kenya, Nairobi, Kenya
| | - Elijah O Juma
- Pan African Mosquito Control Association (PAMCA), Nairobi, Kenya
| | - Charles Mbogo
- Pan African Mosquito Control Association (PAMCA), Nairobi, Kenya
| | - Sonia Barasa
- Malaria Vector Genomic Surveillance, Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK
- Pan African Mosquito Control Association (PAMCA), Nairobi, Kenya
| | - Alistair Miles
- Malaria Vector Genomic Surveillance, Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK
| |
Collapse
|
19
|
Soboleva ES, Kirilenko KM, Fedorova VS, Kokhanenko AA, Artemov GN, Sharakhov IV. Two Nested Inversions in the X Chromosome Differentiate the Dominant Malaria Vectors in Europe, Anopheles atroparvus and Anopheles messeae. INSECTS 2024; 15:312. [PMID: 38786868 PMCID: PMC11122324 DOI: 10.3390/insects15050312] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Revised: 04/22/2024] [Accepted: 04/23/2024] [Indexed: 05/25/2024]
Abstract
The Maculipennis subgroup of malaria mosquitoes includes both dominant malaria vectors and non-vectors in Eurasia. Understanding the genetic factors, particularly chromosomal inversions, that differentiate Anopheles species can provide valuable insights for vector control strategies. Although autosomal inversions between the species in this subgroup have been characterized based on the chromosomal banding patterns, the number and positions of rearrangements in the X chromosome remain unclear due to the divergent banding patterns. Here, we identified two large X chromosomal inversions, approximately 13 Mb and 10 Mb in size, using fluorescence in situ hybridization. The inversion breakpoint regions were mapped by hybridizing 53 gene markers with polytene chromosomes of An. messeae. The DNA probes were designed based on gene sequences from the annotated An. atroparvus genome. The two nested inversions resulted in five syntenic blocks. Only two small syntenic blocks, which encompass 181 annotated genes in the An. atroparvus genome, changed their position and orientation in the X chromosome. The analysis of the An. atroparvus genome revealed an enrichment of gene ontology terms associated with immune system and mating behavior in the rearranged syntenic blocks. Additionally, the enrichment of DNA transposons was found in sequences homologous to three of the four breakpoint regions. This study demonstrates the successful application of the physical genome mapping approach to identify rearrangements that differentiate species in insects with polytene chromosomes.
Collapse
Affiliation(s)
- Evgenia S. Soboleva
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk 634050, Russia
| | - Kirill M. Kirilenko
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk 634050, Russia
| | - Valentina S. Fedorova
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk 634050, Russia
| | - Alina A. Kokhanenko
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk 634050, Russia
| | - Gleb N. Artemov
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk 634050, Russia
| | - Igor V. Sharakhov
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk 634050, Russia
- Department of Entomology, the Fralin Life Sciences Institute, Virginia Polytechnic Institute and State University, 360 West Campus Drive, Blacksburg, VA 24061, USA
| |
Collapse
|
20
|
Dietz L, Mayer C, Stolle E, Eberle J, Misof B, Podsiadlowski L, Niehuis O, Ahrens D. Metazoa-level USCOs as markers in species delimitation and classification. Mol Ecol Resour 2024; 24:e13921. [PMID: 38146909 DOI: 10.1111/1755-0998.13921] [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/15/2023] [Revised: 12/06/2023] [Accepted: 12/13/2023] [Indexed: 12/27/2023]
Abstract
Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.
Collapse
Affiliation(s)
- Lars Dietz
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change, Bonn, Germany
| | - Christoph Mayer
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change, Bonn, Germany
| | - Eckart Stolle
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change, Bonn, Germany
| | - Jonas Eberle
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change, Bonn, Germany
- Paris-Lodron-University, Salzburg, Austria
| | - Bernhard Misof
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change, Bonn, Germany
- Rheinische Friedrich-Wilhelms-Universität Bonn, Bonn, Germany
| | - Lars Podsiadlowski
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change, Bonn, Germany
| | - Oliver Niehuis
- Abt. Evolutionsbiologie und Ökologie, Institut für Biologie I, Albert-Ludwigs-Universität Freiburg, Freiburg, Germany
| | - Dirk Ahrens
- Museum A. Koenig, Leibniz Institute for the Analysis of Biodiversity Change, Bonn, Germany
| |
Collapse
|
21
|
Bossert S, Pauly A, Danforth BN, Orr MC, Murray EA. Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae). Mol Ecol Resour 2024; 24:e13925. [PMID: 38183389 DOI: 10.1111/1755-0998.13925] [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/21/2023] [Revised: 12/08/2023] [Accepted: 12/21/2023] [Indexed: 01/08/2024]
Abstract
Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.
Collapse
Affiliation(s)
- Silas Bossert
- Department of Entomology, Washington State University, Pullman, Washington, USA
- Department of Entomology, National Museum of Natural History, Smithsonian Institution, Washington, DC, USA
| | - Alain Pauly
- Royal Belgian Institute of Natural Sciences, O.D. Taxonomy and Phylogeny, Brussels, Belgium
| | - Bryan N Danforth
- Department of Entomology, Cornell University, Ithaca, New York, USA
| | - Michael C Orr
- Entomologie, Staatliches Museum für Naturkunde Stuttgart, Stuttgart, Germany
| | - Elizabeth A Murray
- Department of Entomology, Washington State University, Pullman, Washington, USA
| |
Collapse
|
22
|
Caputo B, De Marco CM, Pichler V, Bottà G, Bennett KL, Clarkson CS, Tennessen JA, Weetman D, Miles A, Torre AD. Speciation within the Anopheles gambiae complex: high-throughput whole genome sequencing reveals evidence of a putative new cryptic taxon in 'far-west' Africa. RESEARCH SQUARE 2024:rs.3.rs-3914444. [PMID: 38562903 PMCID: PMC10984024 DOI: 10.21203/rs.3.rs-3914444/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae - are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as "intermediate" between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1,190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identify a novel putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, these populations lack insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.
Collapse
Affiliation(s)
- B Caputo
- Dipartimento di Sanità Pubblica e Malattie Infettive, Istituto Pasteur Italia-Fondazione Cenci-Bolognetti, Università di Roma "Sapienza", Rome Italy
| | - C M De Marco
- Dipartimento di Sanità Pubblica e Malattie Infettive, Istituto Pasteur Italia-Fondazione Cenci-Bolognetti, Università di Roma "Sapienza", Rome Italy
| | - V Pichler
- Dipartimento di Sanità Pubblica e Malattie Infettive, Istituto Pasteur Italia-Fondazione Cenci-Bolognetti, Università di Roma "Sapienza", Rome Italy
| | - G Bottà
- Dipartimento di Sanità Pubblica e Malattie Infettive, Istituto Pasteur Italia-Fondazione Cenci-Bolognetti, Università di Roma "Sapienza", Rome Italy
| | - K L Bennett
- Wellcome Sanger Genomic Surveillance Unit, Wellcome Sanger Institute, Hinxton, Cambridge, UK
| | - C S Clarkson
- Wellcome Sanger Genomic Surveillance Unit, Wellcome Sanger Institute, Hinxton, Cambridge, UK
| | - J A Tennessen
- Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Broad Institute, Cambridge, MA, USA
| | - D Weetman
- Liverpool School of Tropical Medicine, Liverpool, UK
| | - A Miles
- Wellcome Sanger Genomic Surveillance Unit, Wellcome Sanger Institute, Hinxton, Cambridge, UK
| | - A Della Torre
- Dipartimento di Sanità Pubblica e Malattie Infettive, Istituto Pasteur Italia-Fondazione Cenci-Bolognetti, Università di Roma "Sapienza", Rome Italy
| |
Collapse
|
23
|
Bernabeu M, Manzano-Morales S, Gabaldón T. On the impact of incomplete taxon sampling on the relative timing of gene transfer events. PLoS Biol 2024; 22:e3002460. [PMID: 38498548 PMCID: PMC10947667 DOI: 10.1371/journal.pbio.3002460] [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: 07/24/2023] [Accepted: 12/04/2023] [Indexed: 03/20/2024] Open
Abstract
A recent study questioned the use of branch length methods to assess the relative timing of horizontal gene transfers because of the effects of so-called "ghost" lineages. This Formal Comment discusses key considerations regarding the potential effect of missing lineages when assessing relative timing of evolutionary events.
Collapse
Affiliation(s)
- Moisès Bernabeu
- Barcelona Supercomputing Center (BSC), Barcelona, Spain
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Saioa Manzano-Morales
- Barcelona Supercomputing Center (BSC), Barcelona, Spain
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Toni Gabaldón
- Barcelona Supercomputing Center (BSC), Barcelona, Spain
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Barcelona, Spain
- Catalan Institution for Research and Advanced Studies (ICREA), Barcelona, Spain
- Centro de Investigación Biomédica En Red de Enfermedades Infecciosas (CIBERINFEC), Barcelona, Spain
| |
Collapse
|
24
|
Jiang Z, Zang W, Ericson PGP, Song G, Wu S, Feng S, Drovetski SV, Liu G, Zhang D, Saitoh T, Alström P, Edwards SV, Lei F, Qu Y. Gene flow and an anomaly zone complicate phylogenomic inference in a rapidly radiated avian family (Prunellidae). BMC Biol 2024; 22:49. [PMID: 38413944 PMCID: PMC10900574 DOI: 10.1186/s12915-024-01848-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: 11/11/2023] [Accepted: 02/15/2024] [Indexed: 02/29/2024] Open
Abstract
BACKGROUND Resolving the phylogeny of rapidly radiating lineages presents a challenge when building the Tree of Life. An Old World avian family Prunellidae (Accentors) comprises twelve species that rapidly diversified at the Pliocene-Pleistocene boundary. RESULTS Here we investigate the phylogenetic relationships of all species of Prunellidae using a chromosome-level de novo assembly of Prunella strophiata and 36 high-coverage resequenced genomes. We use homologous alignments of thousands of exonic and intronic loci to build the coalescent and concatenated phylogenies and recover four different species trees. Topology tests show a large degree of gene tree-species tree discordance but only 40-54% of intronic gene trees and 36-75% of exonic genic trees can be explained by incomplete lineage sorting and gene tree estimation errors. Estimated branch lengths for three successive internal branches in the inferred species trees suggest the existence of an empirical anomaly zone. The most common topology recovered for species in this anomaly zone was not similar to any coalescent or concatenated inference phylogenies, suggesting presence of anomalous gene trees. However, this interpretation is complicated by the presence of gene flow because extensive introgression was detected among these species. When exploring tree topology distributions, introgression, and regional variation in recombination rate, we find that many autosomal regions contain signatures of introgression and thus may mislead phylogenetic inference. Conversely, the phylogenetic signal is concentrated to regions with low-recombination rate, such as the Z chromosome, which are also more resistant to interspecific introgression. CONCLUSIONS Collectively, our results suggest that phylogenomic inference should consider the underlying genomic architecture to maximize the consistency of phylogenomic signal.
Collapse
Affiliation(s)
- Zhiyong Jiang
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
| | - Wenqing Zang
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
| | - Per G P Ericson
- Department of Bioinformatics and Genetics, Swedish Museum of Natural History, PO Box 50007, Stockholm, SE-104 05, Sweden
| | - Gang Song
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| | - Shaoyuan Wu
- Jiangsu International Joint Center of Genomics, Jiangsu Key Laboratory of Phylogenomics & Comparative Genomics, School of Life Sciences, Jiangsu Normal University, Xuzhou, 221116, Jiangsu, China
| | - Shaohong Feng
- Center for Evolutionary & Organismal Biology, Zhejiang University School of Medicine, Hangzhou, 310058, China
- Liangzhu Laboratory, Zhejiang University, 1369 West Wenyi Road, Hangzhou, 311121, China
- Innovation Center of Yangtze River Delta, Zhejiang University, Jiashan, 314102, China
| | - Sergei V Drovetski
- National Museum of Natural History, Smithsonian Institution, Washington, DC, 20004, USA
- Present address: U.S. Geological Survey, Eastern Ecological Science Center at Patuxent Research Refuge, Laurel, MD, 20708, USA
| | - Gang Liu
- Chinese Academy of Forestry, Institute of Ecological Conservation and Restoration, Beijing, 100091, China
| | - Dezhi Zhang
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| | - Takema Saitoh
- Yamashina Institute for Ornithology, Abiko, Chiba, Japan
| | - Per Alström
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- Animal Ecology, Department of Ecology and Genetics, Evolutionary Biology Centre, Uppsala University, Norbyvägen 18 D, 752 36, Uppsala, Sweden
| | - Scott V Edwards
- Museum of Comparative Zoology and Department of Organismic & Evolutionary Biology, Harvard University, 26 Oxford Street, Cambridge, MA, 02138, USA
| | - Fumin Lei
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
| | - Yanhua Qu
- Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China.
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China.
- Department of Bioinformatics and Genetics, Swedish Museum of Natural History, PO Box 50007, Stockholm, SE-104 05, Sweden.
| |
Collapse
|
25
|
Kamau L, Bennett KL, Ochomo E, Herren J, Agumba S, Otieno S, Omoke D, Matoke-Muhia D, Mburu D, Mwangangi J, Ramaita E, Juma EO, Mbogo C, Barasa S, Miles A. The Anopheles coluzzii range extends into Kenya: Detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa. RESEARCH SQUARE 2024:rs.3.rs-3953608. [PMID: 38410447 PMCID: PMC10896386 DOI: 10.21203/rs.3.rs-3953608/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/28/2024]
Abstract
Background Anopheles coluzzii is a primary vector of malaria found in West and Central Africa, but its presence has hitherto never been documented in Kenya. A thorough understanding of vector bionomics is important as it enables the implementation of targeted and effective vector control interventions. Malaria vector surveillance efforts in the country have tended to focus on historically known primary vectors. In the current study, we sought to determine the taxonomic status of samples collected from five different malaria epidemiological zones in Kenya as well asdescribe the population genetic structure and insecticide resistance profiles in relation to other An. coluzzi populations. Methods Mosquitoes were sampled as larvae from Busia, Kwale, Turkana, Kirinyaga and Kiambu counties, representing the range of malaria endemicities in Kenya, in 2019 and 2021 and emergent adults analysed using Whole Genome Sequencing data processed in accordance with the Anopheles gambiae 1000 Genomes Project phase 3. Where available, historical samples from the same sites were included for WGS. Results This study reports the detection of Anopheles coluzzii for the first time in Kenya. The species was detected in Turkana County across all three time points sampled and its presence confirmed through taxonomic analysis. Additionally, we found a lack of strong population genetic differentiation between An. coluzzii from Kenya and those from the more northerly regions of West and Central Africa, suggesting they represent a connected extension to the known species range. Mutations associated with target-site resistance to DDT and pyrethroids and metabolic resistance to DDT were found at high frequencies of ~60%. The profile and frequencies of the variants observed were similar to An. coluzzii from West and Central Africa but the ace-1 mutation linked to organophosphate and carbamate resistance present in An. coluzzii from coastal West Africa was absent in Kenya. Conclusions These findings emphasise the need for the incorporation of genomics in comprehensive and routine vector surveillance to inform on the range of malaria vector species, and their insecticide resistance status to inform the choice of effective vector control approaches.
Collapse
Affiliation(s)
- Luna Kamau
- Centre for Biotechnology Research and Development (CBRD), Kenya Medical Research Institute
| | - Kelly L Bennett
- Malaria Vector Genomic Surveillance, Wellcome Trust Sanger Institute
| | - Eric Ochomo
- Centre for Global Health Research (CGHR), Kenya Medical Research Institute
| | - Jeremy Herren
- International Centre of Insect Physiology and Ecology
| | - Silas Agumba
- Centre for Global Health Research (CGHR), Kenya Medical Research Institute
| | - Samson Otieno
- Centre for Global Health Research (CGHR), Kenya Medical Research Institute
| | - Diana Omoke
- Centre for Biotechnology Research and Development (CBRD), Kenya Medical Research Institute
| | - Damaris Matoke-Muhia
- Centre for Biotechnology Research and Development (CBRD), Kenya Medical Research Institute
| | | | - Joseph Mwangangi
- Centre for Geographic Medicine Research-Coast (CGMR-C), Kenya Medical Research Institute
| | - Edith Ramaita
- Ministry of Health - National Malaria Control Programme (NMCP)
| | | | | | - Sonia Barasa
- Malaria Vector Genomic Surveillance, Wellcome Trust Sanger Institute
| | - Alistair Miles
- Malaria Vector Genomic Surveillance, Wellcome Trust Sanger Institute
| |
Collapse
|
26
|
Foley NM, Harris AJ, Bredemeyer KR, Ruedi M, Puechmaille SJ, Teeling EC, Criscitiello MF, Murphy WJ. Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation. CELL GENOMICS 2024; 4:100482. [PMID: 38237599 PMCID: PMC10879000 DOI: 10.1016/j.xgen.2023.100482] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/05/2023] [Revised: 09/17/2023] [Accepted: 12/10/2023] [Indexed: 02/17/2024]
Abstract
The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.
Collapse
Affiliation(s)
- Nicole M Foley
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX, USA.
| | - Andrew J Harris
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX, USA; Interdisciplinary Program in Genetics & Genomics, Texas A&M University, College Station, TX, USA
| | - Kevin R Bredemeyer
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX, USA; Interdisciplinary Program in Genetics & Genomics, Texas A&M University, College Station, TX, USA
| | - Manuel Ruedi
- Department of Mammalogy and Ornithology, Natural History Museum of Geneva, Route de Malagnou 1, BP 6434, 1211 Geneva 6, Switzerland
| | - Sebastien J Puechmaille
- Institut des Sciences de l'Évolution, Montpellier (ISEM), Université de Montpellier, CNRS, EPHE, IRD, Montpellier, France; Institut Universitaire de France, Paris, France
| | - Emma C Teeling
- School of Biology and Environmental, Science, Science Centre West, University College Dublin, Belfield, Ireland
| | - Michael F Criscitiello
- Interdisciplinary Program in Genetics & Genomics, Texas A&M University, College Station, TX, USA; Department of Veterinary Pathobiology, Texas A&M University, College Station, TX 77843, USA
| | - 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.
| |
Collapse
|
27
|
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
|
28
|
Cribbie EP, Doerr D, Chauve C. AGO, a Framework for the Reconstruction of Ancestral Syntenies and Gene Orders. Methods Mol Biol 2024; 2802:247-265. [PMID: 38819563 DOI: 10.1007/978-1-0716-3838-5_10] [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
Reconstructing ancestral gene orders from the genome data of extant species is an important problem in comparative and evolutionary genomics. In a phylogenomics setting that accounts for gene family evolution through gene duplication and gene loss, the reconstruction of ancestral gene orders involves several steps, including multiple sequence alignment, the inference of reconciled gene trees, and the inference of ancestral syntenies and gene adjacencies. For each of the steps of such a process, several methods can be used and implemented using a growing corpus of, often parameterized, tools; in practice, interfacing such tools into an ancestral gene order reconstruction pipeline is far from trivial. This chapter introduces AGO, a Python-based framework aimed at creating ancestral gene order reconstruction pipelines allowing to interface and parameterize different bioinformatics tools. The authors illustrate the features of AGO by reconstructing ancestral gene orders for the X chromosome of three ancestral Anopheles species using three different pipelines. AGO is freely available at https://github.com/cchauve/AGO-pipeline .
Collapse
Affiliation(s)
- Evan P Cribbie
- Department of Mathematics, Simon Fraser University, Burnaby, BC, Canada
| | - Daniel Doerr
- Department for Endocrinology and Diabetology, Medical Faculty and University Hospital Düsseldorf, German Diabetes Center (DDZ), Leibniz Institute for Diabetes Research, and Center for Digital Medicine, Heinrich Heine University, Düsseldorf, Germany
| | - Cedric Chauve
- Department of Mathematics, Simon Fraser University, Burnaby, BC, Canada.
| |
Collapse
|
29
|
Thawornwattana Y, Seixas F, Yang Z, Mallet J. Major patterns in the introgression history of Heliconius butterflies. eLife 2023; 12:RP90656. [PMID: 38108819 PMCID: PMC10727504 DOI: 10.7554/elife.90656] [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] [Indexed: 12/19/2023] Open
Abstract
Gene flow between species, although usually deleterious, is an important evolutionary process that can facilitate adaptation and lead to species diversification. It also makes estimation of species relationships difficult. Here, we use the full-likelihood multispecies coalescent (MSC) approach to estimate species phylogeny and major introgression events in Heliconius butterflies from whole-genome sequence data. We obtain a robust estimate of species branching order among major clades in the genus, including the 'melpomene-silvaniform' group, which shows extensive historical and ongoing gene flow. We obtain chromosome-level estimates of key parameters in the species phylogeny, including species divergence times, present-day and ancestral population sizes, as well as the direction, timing, and intensity of gene flow. Our analysis leads to a phylogeny with introgression events that differ from those obtained in previous studies. We find that Heliconius aoede most likely represents the earliest-branching lineage of the genus and that 'silvaniform' species are paraphyletic within the melpomene-silvaniform group. Our phylogeny provides new, parsimonious histories for the origins of key traits in Heliconius, including pollen feeding and an inversion involved in wing pattern mimicry. Our results demonstrate the power and feasibility of the full-likelihood MSC approach for estimating species phylogeny and key population parameters despite extensive gene flow. The methods used here should be useful for analysis of other difficult species groups with high rates of introgression.
Collapse
Affiliation(s)
| | - Fernando Seixas
- Department of Organismic and Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College LondonLondonUnited Kingdom
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard UniversityCambridgeUnited States
| |
Collapse
|
30
|
Tan X, Qi J, Liu Z, Fan P, Liu G, Zhang L, Shen Y, Li J, Roos C, Zhou X, Li M. Phylogenomics Reveals High Levels of Incomplete Lineage Sorting at the Ancestral Nodes of the Macaque Radiation. Mol Biol Evol 2023; 40:msad229. [PMID: 37823401 PMCID: PMC10638670 DOI: 10.1093/molbev/msad229] [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/23/2023] [Revised: 09/06/2023] [Accepted: 10/08/2023] [Indexed: 10/13/2023] Open
Abstract
The genus Macaca includes 23 species assigned into 4 to 7 groups. It exhibits the largest geographic range and represents the most successful example of adaptive radiation of nonhuman primates. However, intrageneric phylogenetic relationships among species remain controversial and have not been resolved so far. In this study, we conducted a phylogenomic analysis on 16 newly generated and 8 published macaque genomes. We found strong evidence supporting the division of this genus into 7 species groups. Incomplete lineage sorting (ILS) was the primary factor contributing to the discordance observed among gene trees; however, we also found evidence of hybridization events, specifically between the ancestral arctoides/sinica and silenus/nigra lineages that resulted in the hybrid formation of the fascicularis/mulatta group. Combined with fossil data, our phylogenomic data were used to establish a scenario for macaque radiation. These findings provide insights into ILS and potential ancient introgression events that were involved in the radiation of macaques, which will lead to a better understanding of the rapid speciation occurring in nonhuman primates.
Collapse
Affiliation(s)
- Xinxin Tan
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
- Geneplus-Beijing Institute, Beijing 102206, China
| | - Jiwei Qi
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Zhijin Liu
- College of Life Sciences, Capital Normal University, Beijing 100049, China
| | - Pengfei Fan
- School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, China
| | - Gaoming Liu
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Liye Zhang
- Primate Genetics Laboratory, German Primate Center, Leibniz Institute for Primate Research, Göttingen 37077, Germany
| | - Ying Shen
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Jing Li
- Key Laboratory of Bio-resources and Eco-environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu 610064, China
| | - Christian Roos
- Primate Genetics Laboratory, German Primate Center, Leibniz Institute for Primate Research, Göttingen 37077, Germany
- Gene Bank of Primates, German Primate Center, Leibniz Institute for Primate Research, Göttingen 37077, Germany
| | - Xuming Zhou
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Ming Li
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| |
Collapse
|
31
|
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
|
32
|
Kormos A, Dimopoulos G, Bier E, Lanzaro GC, Marshall JM, James AA. Conceptual risk assessment of mosquito population modification gene-drive systems to control malaria transmission: preliminary hazards list workshops. Front Bioeng Biotechnol 2023; 11:1261123. [PMID: 37965050 PMCID: PMC10641379 DOI: 10.3389/fbioe.2023.1261123] [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: 07/18/2023] [Accepted: 10/09/2023] [Indexed: 11/16/2023] Open
Abstract
The field-testing and eventual adoption of genetically-engineered mosquitoes (GEMs) to control vector-borne pathogen transmission will require them meeting safety criteria specified by regulatory authorities in regions where the technology is being considered for use and other locales that might be impacted. Preliminary risk considerations by researchers and developers may be useful for planning the baseline data collection and field research used to address the anticipated safety concerns. Part of this process is to identify potential hazards (defined as the inherent ability of an entity to cause harm) and their harms, and then chart the pathways to harm and evaluate their probability as part of a risk assessment. The University of California Malaria Initiative (UCMI) participated in a series of workshops held to identify potential hazards specific to mosquito population modification strains carrying gene-drive systems coupled to anti-parasite effector genes and their use in a hypothetical island field trial. The hazards identified were placed within the broader context of previous efforts discussed in the scientific literature. Five risk areas were considered i) pathogens, infections and diseases, and the impacts of GEMs on human and animal health, ii) invasiveness and persistence of GEMs, and interactions of GEMs with target organisms, iii) interactions of GEMs with non-target organisms including horizontal gene transfer, iv) impacts of techniques used for the management of GEMs and v) evolutionary and stability considerations. A preliminary hazards list (PHL) was developed and is made available here. This PHL is useful for internal project risk evaluation and is available to regulators at prospective field sites. UCMI project scientists affirm that the subsequent processes associated with the comprehensive risk assessment for the application of this technology should be driven by the stakeholders at the proposed field site and areas that could be affected by this intervention strategy.
Collapse
Affiliation(s)
- Ana Kormos
- Vector Genetics Laboratory, University of California, Davis, Davis, CA, United States
| | - George Dimopoulos
- W. Harry Feinstone Department of Molecular Microbiology and Immunology, Bloomberg School of Public Health, Malaria Research Institute, Johns Hopkins University, Baltimore, MD, United States
| | - Ethan Bier
- Department of Cell and Developmental Biology, University of California, San Diego, San Diego, CA, United States
| | - Gregory C. Lanzaro
- Vector Genetics Laboratory, University of California, Davis, Davis, CA, United States
| | - John M. Marshall
- Divisions of Epidemiology and Biostatistics, School of Public Health, University of California, Berkeley, Berkeley, CA, United States
| | - Anthony A. James
- Departments of Microbiology and Molecular Genetics and Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
| |
Collapse
|
33
|
Acford-Palmer H, Campos M, Bandibabone J, N'Do S, Bantuzeko C, Zawadi B, Walker T, Phelan JE, Messenger LA, Clark TG, Campino S. Detection of insecticide resistance markers in Anopheles funestus from the Democratic Republic of the Congo using a targeted amplicon sequencing panel. Sci Rep 2023; 13:17363. [PMID: 37833354 PMCID: PMC10575962 DOI: 10.1038/s41598-023-44457-0] [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/18/2023] [Accepted: 10/09/2023] [Indexed: 10/15/2023] Open
Abstract
Vector control strategies have been successful in reducing the number of malaria cases and deaths globally, but the spread of insecticide resistance represents a significant threat to disease control. Insecticide resistance has been reported across Anopheles (An.) vector populations, including species within the An. funestus group. These mosquitoes are responsible for intense malaria transmission across sub-Saharan Africa, including in the Democratic Republic of the Congo (DRC), a country contributing > 12% of global malaria infections and mortality events. To support the continuous efficacy of vector control strategies, it is essential to monitor insecticide resistance using molecular surveillance tools. In this study, we developed an amplicon sequencing ("Amp-seq") approach targeting An. funestus, and using multiplex PCR, dual index barcoding, and next-generation sequencing for high throughput and low-cost applications. Using our Amp-seq approach, we screened 80 An. funestus field isolates from the DRC across a panel of nine genes with mutations linked to insecticide resistance (ace-1, CYP6P4, CYP6P9a, GSTe2, vgsc, and rdl) and mosquito speciation (cox-1, mtND5, and ITS2). Amongst the 18 non-synonymous mutations detected, was N485I, in the ace-1 gene associated with carbamate resistance. Overall, our panel represents an extendable and much-needed method for the molecular surveillance of insecticide resistance in An. funestus populations.
Collapse
Affiliation(s)
- Holly Acford-Palmer
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK
| | - Monica Campos
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK
| | - Janvier Bandibabone
- Centre de Recherche en Sciences Naturelles de Lwiro, Sud-Kivu, Democratic Republic of the Congo
| | - Sévérin N'Do
- Médecins Sans Frontières (MSF) OCBA, Barcelona, Spain
- Institut de Recherche en Sciences de La Santé (IRSS), Bobo-Dioulasso, Burkina Faso
| | - Chimanuka Bantuzeko
- Centre de Recherche en Sciences Naturelles de Lwiro, Sud-Kivu, Democratic Republic of the Congo
- Université Officielle de Bukavu (UOB), Bukavu, Democratic Republic of the Congo
| | - Bertin Zawadi
- Centre de Recherche en Sciences Naturelles de Lwiro, Sud-Kivu, Democratic Republic of the Congo
| | - Thomas Walker
- School of Life Sciences, Gibbet Hill Campus, University of Warwick, Coventry, CV4 7AL, UK
| | - Jody E Phelan
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK
| | - Louisa A Messenger
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK
- Department of Environmental and Occupational Health, School of Public Health, University of Nevada, Las Vegas, Las Vegas, USA
| | - Taane G Clark
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK
- Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, UK
| | - Susana Campino
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK.
| |
Collapse
|
34
|
Soghigian J, Sither C, Justi SA, Morinaga G, Cassel BK, Vitek CJ, Livdahl T, Xia S, Gloria-Soria A, Powell JR, Zavortink T, Hardy CM, Burkett-Cadena ND, Reeves LE, Wilkerson RC, Dunn RR, Yeates DK, Sallum MA, Byrd BD, Trautwein MD, Linton YM, Reiskind MH, Wiegmann BM. Phylogenomics reveals the history of host use in mosquitoes. Nat Commun 2023; 14:6252. [PMID: 37803007 PMCID: PMC10558525 DOI: 10.1038/s41467-023-41764-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2023] [Accepted: 09/08/2023] [Indexed: 10/08/2023] Open
Abstract
Mosquitoes have profoundly affected human history and continue to threaten human health through the transmission of a diverse array of pathogens. The phylogeny of mosquitoes has remained poorly characterized due to difficulty in taxonomic sampling and limited availability of genomic data beyond the most important vector species. Here, we used phylogenomic analysis of 709 single copy ortholog groups from 256 mosquito species to produce a strongly supported phylogeny that resolves the position of the major disease vector species and the major mosquito lineages. Our analyses support an origin of mosquitoes in the early Triassic (217 MYA [highest posterior density region: 188-250 MYA]), considerably older than previous estimates. Moreover, we utilize an extensive database of host associations for mosquitoes to show that mosquitoes have shifted to feeding upon the blood of mammals numerous times, and that mosquito diversification and host-use patterns within major lineages appear to coincide in earth history both with major continental drift events and with the diversification of vertebrate classes.
Collapse
Affiliation(s)
- John Soghigian
- Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC, USA
- Faculty of Veterinary Medicine, University of Calgary, Calgary, AB, Canada
| | - Charles Sither
- Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC, USA
| | - Silvia Andrade Justi
- Walter Reed Biosystematics Unit, Smithsonian Institution Museum Support Center, Suitland, MD, USA
- One Health Branch, Walter Reed Army Institute of Research, Silver Spring, MD, USA
- Department of Entomology, Smithsonian Institution National Museum of Natural History, Washington, DC, USA
| | - Gen Morinaga
- Faculty of Veterinary Medicine, University of Calgary, Calgary, AB, Canada
| | - Brian K Cassel
- Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC, USA
| | - Christopher J Vitek
- Center for Vector-Borne Diseases, University of Texas Rio Grande Valley, Edinburg, TX, USA
| | - Todd Livdahl
- Department of Biology, Clark University, Worcester, MA, USA
| | - Siyang Xia
- Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT, USA
| | - Andrea Gloria-Soria
- Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT, USA
- Department of Entomology, Center for Vector Biology & Zoonotic Diseases, The Connecticut Agricultural Experiment Station, New Haven, CT, USA
| | - Jeffrey R Powell
- Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT, USA
| | - Thomas Zavortink
- Bohart Museum of Entomology, University of California, Davis, CA, USA
| | | | - Nathan D Burkett-Cadena
- Florida Medical Entomology Laboratory, Institute of Food and Agricultural Sciences, University of Florida, Vero Beach, FL, USA
| | - Lawrence E Reeves
- Florida Medical Entomology Laboratory, Institute of Food and Agricultural Sciences, University of Florida, Vero Beach, FL, USA
| | - Richard C Wilkerson
- Walter Reed Biosystematics Unit, Smithsonian Institution Museum Support Center, Suitland, MD, USA
| | - Robert R Dunn
- Department of Applied Ecology, North Carolina State University, Raleigh, NC, USA
| | - David K Yeates
- Australian National Insect Collection, CSIRO National Collections and Marine Infrastructure, Canberra, ACT, Australia
| | - Maria Anice Sallum
- Departamento de Epidemiologia, Faculdade de Saude Publica, Universidade de Sao Paulo, Sao Paulo, Brazil
| | - Brian D Byrd
- College of Health and Human Sciences, School of Health Sciences, Western Carolina University, Cullowhee, NC, USA
| | - Michelle D Trautwein
- Entomology Department, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, San Francisco, CA, USA
| | - Yvonne-Marie Linton
- Walter Reed Biosystematics Unit, Smithsonian Institution Museum Support Center, Suitland, MD, USA
- One Health Branch, Walter Reed Army Institute of Research, Silver Spring, MD, USA
- Department of Entomology, Smithsonian Institution National Museum of Natural History, Washington, DC, USA
| | - Michael H Reiskind
- Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC, USA
| | - Brian M Wiegmann
- Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC, USA.
| |
Collapse
|
35
|
Kamimura Y, Nishikawa M, Yamasako J. DNA barcoding of Japanese earwig species (Insecta, Dermaptera), with sequence diversity analyses of three species of Anisolabididae. Biodivers Data J 2023; 11:e107001. [PMID: 38318508 PMCID: PMC10840520 DOI: 10.3897/bdj.11.e107001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Accepted: 09/19/2023] [Indexed: 02/07/2024] Open
Abstract
Dermaptera is a polyneopteran insect order that includes more than 2,000 described species, commonly known as earwigs, that mainly inhabit tropical, subtropical and warm temperate regions. Although 40 species have been found in Japan, their distribution and habitat preferences have remained ambiguous due to sample misidentification, particularly amongst immature specimens. To overcome this problem, we sequenced and analysed the DNA barcoding region of the mitochondrial cytochrome oxidase I gene (cox1) of dermapteran species recorded from Japan. Including publicly available data, 72.5% of known Japanese dermapteran species were subjected to molecular identification. We extensively sampled three wingless species of subfamily Anisolabidinae (Anisolabididae): Anisolabismaritima, Anisolabellamarginalis and Euborelliapallipes. Although these species exhibit similar habitat preferences as semi-synanthropes, A.maritima, a cosmopolitan species with the highest affinity to seashore, had significantly higher sequence diversity than the latter two species, which are considered endemic to East Asia. A similar trend was observed for (at least partly) winged cosmopolitan species of other families. Introgression with the congener Anisolabisseirokui is also suggested for A.maritima. Possible causes of the varying levels of sequence diversity are discussed.
Collapse
Affiliation(s)
| | - Masaru Nishikawa
- Ehime University, Matsuyama, JapanEhime UniversityMatsuyamaJapan
| | - Junsuke Yamasako
- National Agriculture and Food Research Organization, Tsukuba, JapanNational Agriculture and Food Research OrganizationTsukubaJapan
| |
Collapse
|
36
|
Cicconardi F, Milanetti E, Pinheiro de Castro EC, Mazo-Vargas A, Van Belleghem SM, Ruggieri AA, Rastas P, Hanly J, Evans E, Jiggins CD, Owen McMillan W, Papa R, Di Marino D, Martin A, Montgomery SH. Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies. Nat Commun 2023; 14:5620. [PMID: 37699868 PMCID: PMC10497600 DOI: 10.1038/s41467-023-41412-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Accepted: 08/30/2023] [Indexed: 09/14/2023] Open
Abstract
Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.
Collapse
Affiliation(s)
- Francesco Cicconardi
- School of Biological Sciences, Bristol University, Bristol, United Kingdom.
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom.
| | - Edoardo Milanetti
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185, Rome, Italy
- Center for Life Nano- & Neuro-Science, Italian Institute of Technology, Viale Regina Elena 291, 00161, Rome, Italy
| | | | - Anyi Mazo-Vargas
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, 14853, USA
| | - Steven M Van Belleghem
- Department of Biology, University of Puerto Rico, Rio Piedras, PR, Puerto Rico
- Ecology, Evolution and Conservation Biology, Biology Department, KU Leuven, Leuven, Belgium
| | | | - Pasi Rastas
- Institute of Biotechnology, University of Helsinki, Helsinki, Finland
| | - Joseph Hanly
- Department of Biological Sciences, The George Washington University, Washington DC, WA, 20052, USA
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Elizabeth Evans
- Department of Biology, University of Puerto Rico, Rio Piedras, PR, Puerto Rico
| | - Chris D Jiggins
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - W Owen McMillan
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Riccardo Papa
- Department of Biology, University of Puerto Rico, Rio Piedras, PR, Puerto Rico
- Molecular Sciences and Research Center, University of Puerto Rico, San Juan, PR, Puerto Rico
- Comprehensive Cancer Center, University of Puerto Rico, San Juan, PR, Puerto Rico
| | - Daniele Di Marino
- Department of Life and Environmental Sciences, New York-Marche Structural Biology Center (NY-MaSBiC), Polytechnic University of Marche, Via Brecce Bianche, 60131, Ancona, Italy
- Neuronal Death and Neuroprotection Unit, Department of Neuroscience, Mario Negri Institute for Pharmacological Research-IRCCS, Via Mario Negri 2, 20156, Milano, Italy
- National Biodiversity Future Center (NBFC), Palermo, Italy
| | - Arnaud Martin
- Department of Biological Sciences, The George Washington University, Washington DC, WA, 20052, USA
| | - Stephen H Montgomery
- School of Biological Sciences, Bristol University, Bristol, United Kingdom.
- Smithsonian Tropical Research Institute, Panama City, Panama.
| |
Collapse
|
37
|
Thomas GWC, Hughes JJ, Kumon T, Berv JS, Nordgren CE, Lampson M, Levine M, Searle JB, Good JM. The genomic landscape, causes, and consequences of extensive phylogenomic discordance in Old World mice and rats. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.28.555178. [PMID: 37693498 PMCID: PMC10491188 DOI: 10.1101/2023.08.28.555178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/12/2023]
Abstract
A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is often widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents (Old World mice and rats) - a large and ecologically diverse group that gave rise to the mouse and rat model systems. Combining new linked-read genome assemblies for seven murine species with eleven published rodent genomes, we first used ultra-conserved elements (UCEs) to infer a robust species tree. We then used whole genomes to examine finer-scale patterns of discordance and found that phylogenies built from proximate chromosomal regions had similar phylogenies. However, there was no relationship between tree similarity and local recombination rates in house mice, suggesting that genetic linkage influences phylogenetic patterns over deeper timescales. This signal may be independent of contemporary recombination landscapes. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to less discordance, while genes experiencing positive selection showed more discordant and variable phylogenetic signals. Finally, we show that assuming a single species tree can result in high error rates when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.
Collapse
Affiliation(s)
- Gregg W. C. Thomas
- Division of Biological Sciences, University of Montana, Missoula, MT, 59801
- Informatics Group, Harvard University, Cambridge, MA, 02138
| | - Jonathan J. Hughes
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, 14853
- Department of Evolution, Ecology, and Organismal Biology, University of California Riverside, Riverside, CA, 92521
| | - Tomohiro Kumon
- Department of Biology, University of Pennsylvania, Philadelphia, PA, 19104
| | - Jacob S. Berv
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, 14853
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, 48109
| | - C. Erik Nordgren
- Department of Biology, University of Pennsylvania, Philadelphia, PA, 19104
| | - Michael Lampson
- Department of Biology, University of Pennsylvania, Philadelphia, PA, 19104
| | - Mia Levine
- Department of Biology, University of Pennsylvania, Philadelphia, PA, 19104
| | - Jeremy B. Searle
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, 14853
| | - Jeffrey M. Good
- Division of Biological Sciences, University of Montana, Missoula, MT, 59801
| |
Collapse
|
38
|
Giesen A, Blanckenhorn WU, Schäfer MA, Shimizu KK, Shimizu-Inatsugi R, Misof B, Podsiadlowski L, Niehuis O, Lischer HEL, Aeschbacher S, Kapun M. Geographic Variation in Genomic Signals of Admixture Between Two Closely Related European Sepsid Fly Species. Evol Biol 2023; 50:395-412. [PMID: 37854269 PMCID: PMC10579158 DOI: 10.1007/s11692-023-09612-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Accepted: 07/28/2023] [Indexed: 10/20/2023]
Abstract
The extent of interspecific gene flow and its consequences for the initiation, maintenance, and breakdown of species barriers in natural systems remain poorly understood. Interspecific gene flow by hybridization may weaken adaptive divergence, but can be overcome by selection against hybrids, which may ultimately promote reinforcement. An informative step towards understanding the role of gene flow during speciation is to describe patterns of past gene flow among extant species. We investigate signals of admixture between allopatric and sympatric populations of the two closely related European dung fly species Sepsis cynipsea and S. neocynipsea (Diptera: Sepsidae). Based on microsatellite genotypes, we first inferred a baseline demographic history using Approximate Bayesian Computation. We then used genomic data from pooled DNA of natural and laboratory populations to test for past interspecific gene flow based on allelic configurations discordant with the inferred population tree (ABBA-BABA test with D-statistic). Comparing the detected signals of gene flow with the contemporary geographic relationship among interspecific pairs of populations (sympatric vs. allopatric), we made two contrasting observations. At one site in the French Cevennes, we detected an excess of past interspecific gene flow, while at two sites in Switzerland we observed lower signals of past microsatellite genotypes gene flow among populations in sympatry compared to allopatric populations. These results suggest that the species boundaries between these two species depend on the past and/or present eco-geographic context in Europe, which indicates that there is no uniform link between contemporary geographic proximity and past interspecific gene flow in natural populations. Supplementary Information The online version contains supplementary material available at 10.1007/s11692-023-09612-5.
Collapse
Affiliation(s)
- Athene Giesen
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Wolf U. Blanckenhorn
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Martin A. Schäfer
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Kentaro K. Shimizu
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Rie Shimizu-Inatsugi
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Bernhard Misof
- Zoological Research Museum Alexander Koenig, Bonn, Germany
| | | | - Oliver Niehuis
- Department of Evolutionary Biology and Ecology, Institute of Biology I (Zoology), Albert Ludwig University, Freiburg, Germany
| | - Heidi E. L. Lischer
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
- Interfaculty Bioinformatics Unit, University of Bern, Bern, Switzerland
| | - Simon Aeschbacher
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Martin Kapun
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
- Division of Cell & Developmental Biology, Medical University of Vienna, Vienna, Austria
- Natural History Museum Vienna, Burgring 7, 1010 Vienna, Austria
| |
Collapse
|
39
|
Knyshov A, Gordon ERL, Masonick PK, Castillo S, Forero D, Hoey-Chamberlain R, Hwang WS, Johnson KP, Lemmon AR, Moriarty Lemmon E, Standring S, Zhang J, Weirauch C. Chromosome-Aware Phylogenomics of Assassin Bugs (Hemiptera: Reduvioidea) Elucidates Ancient Gene Conflict. Mol Biol Evol 2023; 40:msad168. [PMID: 37494292 PMCID: PMC10411492 DOI: 10.1093/molbev/msad168] [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: 07/15/2023] [Accepted: 07/19/2023] [Indexed: 07/28/2023] Open
Abstract
Though the phylogenetic signal of loci on sex chromosomes can differ from those on autosomes, chromosomal-level genome assemblies for nonvertebrates are still relatively scarce and conservation of chromosomal gene content across deep phylogenetic scales has therefore remained largely unexplored. We here assemble a uniquely large and diverse set of samples (17 anchored hybrid enrichment, 24 RNA-seq, and 70 whole-genome sequencing samples of variable depth) for the medically important assassin bugs (Reduvioidea). We assess the performance of genes based on multiple features (e.g., nucleotide vs. amino acid, nuclear vs. mitochondrial, and autosomal vs. X chromosomal) and employ different methods (concatenation and coalescence analyses) to reconstruct the unresolved phylogeny of this diverse (∼7,000 spp.) and old (>180 Ma) group. Our results show that genes on the X chromosome are more likely to have discordant phylogenies than those on autosomes. We find that the X chromosome conflict is driven by high gene substitution rates that impact the accuracy of phylogenetic inference. However, gene tree clustering showed strong conflict even after discounting variable third codon positions. Alternative topologies were not particularly enriched for sex chromosome loci, but spread across the genome. We conclude that binning genes to autosomal or sex chromosomes may result in a more accurate picture of the complex evolutionary history of a clade.
Collapse
Affiliation(s)
- Alexander Knyshov
- Department of Entomology, University of California, Riverside, CA, USA
| | - Eric R L Gordon
- Ecology and Evolutionary Biology Department, University of Connecticut, Storrs, CT, USA
| | - Paul K Masonick
- Department of Entomology, University of California, Riverside, CA, USA
| | | | - Dimitri Forero
- Instituto de Ciencias Naturales, Universidad Nacional de Colombia, Bogota, Colombia
| | | | - Wei Song Hwang
- Lee Kong Chian Natural History Museum, National University of Singapore, Queenstown, Singapore
| | - Kevin P Johnson
- Illinois Natural History Survey, Prairie Research Institute, University of Illinois, Champaign, IL, USA
| | - Alan R Lemmon
- Department of Scientific Computing, Florida State University, Tallahassee, FL, USA
| | | | | | - Junxia Zhang
- Key Laboratory of Zoological Systematics and Application of Hebei Province, Institute of Life Science and Green Development, College of Life Sciences, Hebei University, Baoding, Hebei, China
| | | |
Collapse
|
40
|
Carballar-Lejarazú R, Dong Y, Pham TB, Tushar T, Corder RM, Mondal A, Sánchez C. HM, Lee HF, Marshall JM, Dimopoulos G, James AA. Dual effector population modification gene-drive strains of the African malaria mosquitoes, Anopheles gambiae and Anopheles coluzzii. Proc Natl Acad Sci U S A 2023; 120:e2221118120. [PMID: 37428915 PMCID: PMC10629562 DOI: 10.1073/pnas.2221118120] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Accepted: 06/05/2023] [Indexed: 07/12/2023] Open
Abstract
Proposed genetic approaches for reducing human malaria include population modification, which introduces genes into vector mosquitoes to reduce or prevent parasite transmission. We demonstrate the potential of Cas9/guide RNA (gRNA)-based gene-drive systems linked to dual antiparasite effector genes to spread rapidly through mosquito populations. Two strains have an autonomous gene-drive system coupled to dual anti-Plasmodium falciparum effector genes comprising single-chain variable fragment monoclonal antibodies targeting parasite ookinetes and sporozoites in the African malaria mosquitoes Anopheles gambiae (AgTP13) and Anopheles coluzzii (AcTP13). The gene-drive systems achieved full introduction within 3 to 6 mo after release in small cage trials. Life-table analyses revealed no fitness loads affecting AcTP13 gene-drive dynamics but AgTP13 males were less competitive than wild types. The effector molecules reduced significantly both parasite prevalence and infection intensities. These data supported transmission modeling of conceptual field releases in an island setting that shows meaningful epidemiological impacts at different sporozoite threshold levels (2.5 to 10 k) for human infection by reducing malaria incidence in optimal simulations by 50 to 90% within as few as 1 to 2 mo after a series of releases, and by ≥90% within 3 mo. Modeling outcomes for low sporozoite thresholds are sensitive to gene-drive system fitness loads, gametocytemia infection intensities during parasite challenges, and the formation of potentially drive-resistant genome target sites, extending the predicted times to achieve reduced incidence. TP13-based strains could be effective for malaria control strategies following validation of sporozoite transmission threshold numbers and testing field-derived parasite strains. These or similar strains are viable candidates for future field trials in a malaria-endemic region.
Collapse
Affiliation(s)
| | - Yuemei Dong
- W. Harry Feinstone Department of Molecular Microbiology and Immunology, Bloomberg School of Public Health, Malaria Research Institute, Johns Hopkins University, Baltimore, MD21205
| | - Thai Binh Pham
- Department of Microbiology & Molecular Genetics, University of California, Irvine, CA92697-4025
| | - Taylor Tushar
- Department of Microbiology & Molecular Genetics, University of California, Irvine, CA92697-4025
| | - Rodrigo M. Corder
- Divisions of Epidemiology and Biostatistics, School of Public Health, University of California, Berkeley, CA94720
| | - Agastya Mondal
- Divisions of Epidemiology and Biostatistics, School of Public Health, University of California, Berkeley, CA94720
| | - Héctor M. Sánchez C.
- Divisions of Epidemiology and Biostatistics, School of Public Health, University of California, Berkeley, CA94720
| | - Hsu-Feng Lee
- Department of Microbiology & Molecular Genetics, University of California, Irvine, CA92697-4025
| | - John M. Marshall
- Divisions of Epidemiology and Biostatistics, School of Public Health, University of California, Berkeley, CA94720
| | - George Dimopoulos
- W. Harry Feinstone Department of Molecular Microbiology and Immunology, Bloomberg School of Public Health, Malaria Research Institute, Johns Hopkins University, Baltimore, MD21205
| | - Anthony A. James
- Department of Microbiology & Molecular Genetics, University of California, Irvine, CA92697-4025
- Department of Molecular Biology & Biochemistry, University of California, Irvine, CA92697-3900
| |
Collapse
|
41
|
Zhang L, Abhari N, Colijn C, Wu Y. A fast and scalable method for inferring phylogenetic networks from trees by aligning lineage taxon strings. Genome Res 2023; 33:1053-1060. [PMID: 37217252 PMCID: PMC10538497 DOI: 10.1101/gr.277669.123] [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: 01/06/2023] [Accepted: 05/16/2023] [Indexed: 05/24/2023]
Abstract
The reconstruction of phylogenetic networks is an important but challenging problem in phylogenetics and genome evolution, as the space of phylogenetic networks is vast and cannot be sampled well. One approach to the problem is to solve the minimum phylogenetic network problem, in which phylogenetic trees are first inferred, and then the smallest phylogenetic network that displays all the trees is computed. The approach takes advantage of the fact that the theory of phylogenetic trees is mature, and there are excellent tools available for inferring phylogenetic trees from a large number of biomolecular sequences. A tree-child network is a phylogenetic network satisfying the condition that every nonleaf node has at least one child that is of indegree one. Here, we develop a new method that infers the minimum tree-child network by aligning lineage taxon strings in the phylogenetic trees. This algorithmic innovation enables us to get around the limitations of the existing programs for phylogenetic network inference. Our new program, named ALTS, is fast enough to infer a tree-child network with a large number of reticulations for a set of up to 50 phylogenetic trees with 50 taxa that have only trivial common clusters in about a quarter of an hour on average.
Collapse
Affiliation(s)
- Louxin Zhang
- Department of Mathematics and Centre for Data Science and Machine Learning, National University of Singapore, Singapore 119076, Singapore;
| | - Niloufar Abhari
- Department of Mathematics, Simon Fraser University, Burnaby, British Columbia V5A 1S6, Canada
| | - Caroline Colijn
- Department of Mathematics, Simon Fraser University, Burnaby, British Columbia V5A 1S6, Canada
| | - Yufeng Wu
- Department of Computer Science and Engineering and Institute for Systems Genomics, University of Connecticut, Storrs, Connecticut 06269, USA
| |
Collapse
|
42
|
Dagilis AJ, Matute DR. The fitness of an introgressing haplotype changes over the course of divergence and depends on its size and genomic location. PLoS Biol 2023; 21:e3002185. [PMID: 37459351 PMCID: PMC10374083 DOI: 10.1371/journal.pbio.3002185] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Revised: 07/27/2023] [Accepted: 06/06/2023] [Indexed: 07/28/2023] Open
Abstract
The genomic era has made clear that introgression, or the movement of genetic material between species, is a common feature of evolution. Examples of both adaptive and deleterious introgression exist in a variety of systems. What is unclear is how the fitness of an introgressing haplotype changes as species diverge or as the size of the introgressing haplotype changes. In a simple model, we show that introgression may more easily occur into parts of the genome which have not diverged heavily from a common ancestor. The key insight is that alleles from a shared genetic background are likely to have positive epistatic interactions, increasing the fitness of a larger introgressing block. In regions of the genome where few existing substitutions are disrupted, this positive epistasis can be larger than incompatibilities with the recipient genome. Further, we show that early in the process of divergence, introgression of large haplotypes can be favored more than introgression of individual alleles. This model is consistent with observations of a positive relationship between recombination rate and introgression frequency across the genome; however, it generates several novel predictions. First, the model suggests that the relationship between recombination rate and introgression may not exist, or may be negative, in recently diverged species pairs. Furthermore, the model suggests that introgression that replaces existing derived variation will be more deleterious than introgression at sites carrying ancestral variants. These predictions are tested in an example of introgression in Drosophila melanogaster, with some support for both. Finally, the model provides a potential alternative explanation to asymmetry in the direction of introgression, with expectations of higher introgression from rapidly diverged populations into slowly evolving ones.
Collapse
Affiliation(s)
- Andrius J Dagilis
- Biology Department, University of North Carolina, Chapel Hill, North Carolina, United States of America
| | - Daniel R Matute
- Biology Department, University of North Carolina, Chapel Hill, North Carolina, United States of America
| |
Collapse
|
43
|
Leal JL, Milesi P, Salojärvi J, Lascoux M. Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences. Syst Biol 2023; 72:372-390. [PMID: 36932679 PMCID: PMC10275558 DOI: 10.1093/sysbio/syad009] [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: 01/05/2022] [Revised: 10/14/2022] [Accepted: 03/10/2023] [Indexed: 03/19/2023] Open
Abstract
Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].
Collapse
Affiliation(s)
- J Luis Leal
- Plant Ecology and Evolution, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| | - Pascal Milesi
- Plant Ecology and Evolution, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
- Science for Life Laboratory (SciLifeLab), Uppsala University, 75237 Uppsala, Sweden
| | - Jarkko Salojärvi
- Organismal and Evolutionary Biology Research Program, Faculty of Biological and Environmental Sciences, and Viikki Plant Science Centre, University of Helsinki, P.O. Box 65 (Viikinkaari 1), 00014 Helsinki, Finland
- School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551, Singapore
| | - Martin Lascoux
- Plant Ecology and Evolution, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
- Science for Life Laboratory (SciLifeLab), Uppsala University, 75237 Uppsala, Sweden
| |
Collapse
|
44
|
Zhang BL, Chen W, Wang Z, Pang W, Luo MT, Wang S, Shao Y, He WQ, Deng Y, Zhou L, Chen J, Yang MM, Wu Y, Wang L, Fernández-Bellon H, Molloy S, Meunier H, Wanert F, Kuderna L, Marques-Bonet T, Roos C, Qi XG, Li M, Liu Z, Schierup MH, Cooper DN, Liu J, Zheng YT, Zhang G, Wu DD. Comparative genomics reveals the hybrid origin of a macaque group. SCIENCE ADVANCES 2023; 9:eadd3580. [PMID: 37262187 PMCID: PMC10413639 DOI: 10.1126/sciadv.add3580] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Accepted: 01/25/2023] [Indexed: 06/03/2023]
Abstract
Although species can arise through hybridization, compelling evidence for hybrid speciation has been reported only rarely in animals. Here, we present phylogenomic analyses on genomes from 12 macaque species and show that the fascicularis group originated from an ancient hybridization between the sinica and silenus groups ~3.45 to 3.56 million years ago. The X chromosomes and low-recombination regions exhibited equal contributions from each parental lineage, suggesting that they were less affected by subsequent backcrossing and hence could have played an important role in maintaining hybrid integrity. We identified many reproduction-associated genes that could have contributed to the development of the mixed sexual phenotypes characteristic of the fascicularis group. The phylogeny within the silenus group was also resolved, and functional experimentation confirmed that all extant Western silenus species are susceptible to HIV-1 infection. Our study provides novel insights into macaque evolution and reveals a hybrid speciation event that has occurred only very rarely in primates.
Collapse
Affiliation(s)
- Bao-Lin Zhang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
| | - Wu Chen
- Guangzhou Zoo and Guangzhou Wildlife Research Center, Guangzhou 510070, China
| | - Zefu Wang
- Key Laboratory for Bio-resource and Eco-environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu 610065, China
- Co-Innovation Center for Sustainable Forestry in Southern China, College of Biology and the Environment, Nanjing Forestry University, Nanjing 210037, China
| | - Wei Pang
- Key Laboratory of Animal Models and Human Disease Mechanisms of the Chinese Academy of Sciences, KIZ-CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
| | - Meng-Ting Luo
- Key Laboratory of Animal Models and Human Disease Mechanisms of the Chinese Academy of Sciences, KIZ-CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
| | - Sheng Wang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
| | - Yong Shao
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
| | - Wen-Qiang He
- Key Laboratory of Animal Models and Human Disease Mechanisms of the Chinese Academy of Sciences, KIZ-CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
| | - Yuan Deng
- BGI-Shenzhen, Shenzhen 518083, China
- Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen DK-2100, Denmark
| | - Long Zhou
- Center for Evolutionary and Organismal Biology and Women’s Hospital at Zhejiang University School of Medicine, Hangzhou 310058, China
| | | | - Min-Min Yang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
| | - Yajiang Wu
- Guangzhou Zoo and Guangzhou Wildlife Research Center, Guangzhou 510070, China
| | - Lu Wang
- Shaanxi Key Laboratory for Animal Conservation, College of Life Sciences, Northwest University, Xi’an, China
| | | | | | - Hélène Meunier
- Centre de Primatologie, de l'Université de Strasbourg, Niederhausbergen, France
- Laboratoire de Neurosciences Cognitives et Adaptatives, UMR 7364, Université de Strasbourg, Strasbourg, France
| | - Fanélie Wanert
- Plateforme SILABE, Université de Strasbourg, Niederhausbergen, France
| | - Lukas Kuderna
- Genome Interpretation Department, Illumina Inc., Foster City, CA, USA
| | - Tomas Marques-Bonet
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, Dr. Aiguader 88, Barcelona 08003, Spain
- Catalan Institution of Research and Advanced Studies (ICREA), Passeig de Lluís Companys, 23, Barcelona 08010, Spain
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, Barcelona 08028, Spain
- Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Edifici ICTA-ICP, c/Columnes s/n, 08193 Cerdanyola del Vallès, Barcelona, Spain
| | - Christian Roos
- Primate Genetics Laboratory, German Primate Center, Göttingen, Germany
- Gene Bank of Primates, German Primate Center, Göttingen, Germany
| | - Xiao-Guang Qi
- Shaanxi Key Laboratory for Animal Conservation, College of Life Sciences, Northwest University, Xi’an, China
| | - Ming Li
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Zhijin Liu
- College of Life Sciences, Capital Normal University, Beijing 100048, China
| | | | - David N. Cooper
- Institute of Medical Genetics, School of Medicine, Cardiff University, Cardiff CF14 4XN, UK
| | - Jianquan Liu
- Key Laboratory for Bio-resource and Eco-environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu 610065, China
- State Key Laboratory of Grassland Agro-ecosystem, Institute of Innovation Ecology and College of Life Sciences, Lanzhou University, Lanzhou 730000, China
| | - Yong-Tang Zheng
- Key Laboratory of Animal Models and Human Disease Mechanisms of the Chinese Academy of Sciences, KIZ-CUHK Joint Laboratory of Bioresources and Molecular Research in Common Diseases, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
- National Resource Center for Non-Human Primates, Kunming Primate Research Center and National Research Facility for Phenotypic and Genetic Analysis of Model Animals (Primate Facility), Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650107, China
| | - Guojie Zhang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
- Center for Evolutionary and Organismal Biology and Women’s Hospital at Zhejiang University School of Medicine, Hangzhou 310058, China
- Liangzhu Laboratory, Zhejiang University Medical Center, 1369 West Wenyi Road, Hangzhou 311121, China
- Villum Center for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen 2100, Denmark
| | - Dong-Dong Wu
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming 650223, China
- National Resource Center for Non-Human Primates, Kunming Primate Research Center and National Research Facility for Phenotypic and Genetic Analysis of Model Animals (Primate Facility), Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650107, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming 650223, China
- Kunming Natural History Museum of Zoology, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650223, China
| |
Collapse
|
45
|
Hibbins MS, Breithaupt LC, Hahn MW. Phylogenomic comparative methods: Accurate evolutionary inferences in the presence of gene tree discordance. Proc Natl Acad Sci U S A 2023; 120:e2220389120. [PMID: 37216509 PMCID: PMC10235958 DOI: 10.1073/pnas.2220389120] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Accepted: 04/24/2023] [Indexed: 05/24/2023] Open
Abstract
Phylogenetic comparative methods have long been a mainstay of evolutionary biology, allowing for the study of trait evolution across species while accounting for their common ancestry. These analyses typically assume a single, bifurcating phylogenetic tree describing the shared history among species. However, modern phylogenomic analyses have shown that genomes are often composed of mosaic histories that can disagree both with the species tree and with each other-so-called discordant gene trees. These gene trees describe shared histories that are not captured by the species tree, and therefore that are unaccounted for in classic comparative approaches. The application of standard comparative methods to species histories containing discordance leads to incorrect inferences about the timing, direction, and rate of evolution. Here, we develop two approaches for incorporating gene tree histories into comparative methods: one that constructs an updated phylogenetic variance-covariance matrix from gene trees, and another that applies Felsenstein's pruning algorithm over a set of gene trees to calculate trait histories and likelihoods. Using simulation, we demonstrate that our approaches generate much more accurate estimates of tree-wide rates of trait evolution than standard methods. We apply our methods to two clades of the wild tomato genus Solanum with varying rates of discordance, demonstrating the contribution of gene tree discordance to variation in a set of floral traits. Our approaches have the potential to be applied to a broad range of classic inference problems in phylogenetics, including ancestral state reconstruction and the inference of lineage-specific rate shifts.
Collapse
Affiliation(s)
- Mark S. Hibbins
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ONM5S 3B2, Canada
- Department of Biology, Indiana University, Bloomington, IN47405
| | - Lara C. Breithaupt
- Department of Biology, Indiana University, Bloomington, IN47405
- Department of Computer Science, Duke University, Durham, NC27710
| | - Matthew W. Hahn
- Department of Biology, Indiana University, Bloomington, IN47405
- Department of Computer Science, Indiana University, Bloomington, IN47405
| |
Collapse
|
46
|
Foley NM, Mason VC, Harris AJ, Bredemeyer KR, Damas J, Lewin HA, Eizirik E, Gatesy J, Karlsson EK, Lindblad-Toh K, Springer MS, Murphy WJ, Andrews G, Armstrong JC, Bianchi M, Birren BW, Bredemeyer KR, Breit AM, Christmas MJ, Clawson H, Damas J, Di Palma F, Diekhans M, Dong MX, Eizirik E, Fan K, Fanter C, Foley NM, Forsberg-Nilsson K, Garcia CJ, Gatesy J, Gazal S, Genereux DP, Goodman L, Grimshaw J, Halsey MK, Harris AJ, Hickey G, Hiller M, Hindle AG, Hubley RM, Hughes GM, Johnson J, Juan D, Kaplow IM, Karlsson EK, Keough KC, Kirilenko B, Koepfli KP, Korstian JM, Kowalczyk A, Kozyrev SV, Lawler AJ, Lawless C, Lehmann T, Levesque DL, Lewin HA, Li X, Lind A, Lindblad-Toh K, Mackay-Smith A, Marinescu VD, Marques-Bonet T, Mason VC, Meadows JRS, Meyer WK, Moore JE, Moreira LR, Moreno-Santillan DD, Morrill KM, Muntané G, Murphy WJ, Navarro A, Nweeia M, Ortmann S, Osmanski A, Paten B, Paulat NS, Pfenning AR, Phan BN, Pollard KS, Pratt HE, Ray DA, Reilly SK, Rosen JR, Ruf I, Ryan L, Ryder OA, Sabeti PC, Schäffer DE, Serres A, Shapiro B, Smit AFA, Springer M, Srinivasan C, Steiner C, Storer JM, Sullivan KAM, Sullivan PF, Sundström E, Supple MA, Swofford R, Talbot JE, Teeling E, Turner-Maier J, Valenzuela A, Wagner F, Wallerman O, Wang C, Wang J, Weng Z, Wilder AP, Wirthlin ME, Xue JR, Zhang X. A genomic timescale for placental mammal evolution. Science 2023; 380:eabl8189. [PMID: 37104581 DOI: 10.1126/science.abl8189] [Citation(s) in RCA: 27] [Impact Index Per Article: 27.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.
Collapse
Affiliation(s)
- Nicole M Foley
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX, USA
| | - Victor C Mason
- Institute of Cell Biology, University of Bern, Bern, Switzerland
| | - Andrew J Harris
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX, USA
- Interdisciplinary Program in Genetics and Genomics, Texas A&M University, College Station, TX, USA
| | - Kevin R Bredemeyer
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX, USA
- Interdisciplinary Program in Genetics and Genomics, Texas A&M University, College Station, TX, USA
| | - Joana Damas
- The Genome Center, University of California, Davis, CA, USA
| | - Harris A Lewin
- The Genome Center, University of California, Davis, CA, USA
- Department of Evolution and Ecology, University of California, Davis, CA, USA
| | - Eduardo Eizirik
- School of Health and Life Sciences, Pontifical Catholic University of Rio Grande do Sul, Porto Alegre, Brazil
| | - John Gatesy
- Division of Vertebrate Zoology, American Museum of Natural History, New York, NY, USA
| | - Elinor K Karlsson
- Program in Bioinformatics and Integrative Biology, UMass Chan Medical School, Worcester, MA 01605, USA
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
- Program in Molecular Medicine, University of Massachussetts Chan Medical School, Worcester, MA 01605, USA
| | - Kerstin Lindblad-Toh
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
- Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, 751 32 Uppsala, Sweden
| | - Mark S Springer
- Department of Evolution, Ecology, and Organismal Biology, University of California, Riverside, CA, USA
| | - William J Murphy
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX, USA
- Interdisciplinary Program in Genetics and Genomics, Texas A&M University, College Station, TX, USA
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
47
|
Ravazi A, de Oliveira J, Madeira FF, Nunes GM, dos Reis YV, de Oliveira ABB, Azevedo LMS, Galvão C, de Azeredo-Oliveira MTV, da Rosa JA, Alevi KCC. Climate and Environmental Changes and Their Potential Effects on the Dynamics of Chagas Disease: Hybridization in Rhodniini (Hemiptera, Triatominae). INSECTS 2023; 14:378. [PMID: 37103193 PMCID: PMC10143345 DOI: 10.3390/insects14040378] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Revised: 02/24/2023] [Accepted: 03/01/2023] [Indexed: 06/19/2023]
Abstract
Chagas disease affects about eight million people. In view of the issues related to the influence of anthropogenic changes in the dynamics of the distribution and reproductive interaction of triatomines, we performed experimental crosses between species of the Rhodniini tribe in order to evaluate interspecific reproductive interactions and hybrid production capacity. Reciprocal crossing experiments were conducted among Rhodnius brethesi × R. pictipes, R. colombiensis × R. ecuadoriensis, R. neivai × R. prolixus, R. robustus × R. prolixus, R. montenegrensis × R. marabaensis; R. montenegrensis × R. robustus, R. prolixus × R. nasutus and R. neglectus × R. milesi. With the exception of crosses between R. pictipes ♀ × R. brethesi ♂, R. ecuadoriensis ♀ × R. colombiensis ♂ and R. prolixus ♀ × R. neivai ♂, all experimental crosses resulted in hybrids. Our results demonstrate that both allopatric and sympatric species produce hybrids, which can generate concern for public health agencies in the face of current anthropogenic events. Thus, we demonstrate that species of the Rhodniini tribe are capable of producing hybrids under laboratory conditions. These results are of great epidemiological importance and raise an important discussion about the influence of climatic and environmental interactions on Chagas disease dynamics.
Collapse
Affiliation(s)
- Amanda Ravazi
- Instituto de Biociências de Botucatu, Universidade Estadual Paulista “Júlio de Mesquita Filho” (UNESP), Rua Dr. Antônio Celso Wagner Zanin, 250, Distrito de Rubião Junior, Botucatu 18618-689, SP, Brazil
| | - Jader de Oliveira
- Laboratório de Entomologia em Saúde Pública, Faculdade de Saúde Pública, Universidade de São Paulo (USP), Av. Dr. Arnaldo 715, São Paulo 01246-904, SP, Brazil
| | - Fernanda Fernandez Madeira
- Laboratório de Biologia Celular, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista “Júlio de Mesquita Filho” (UNESP), Rua Cristóvão Colombo 2265, São José do Rio Preto 15054-000, SP, Brazil
| | - Giovana Menezes Nunes
- Laboratório de Biologia Celular, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista “Júlio de Mesquita Filho” (UNESP), Rua Cristóvão Colombo 2265, São José do Rio Preto 15054-000, SP, Brazil
| | - Yago Visinho dos Reis
- Instituto de Biociências de Botucatu, Universidade Estadual Paulista “Júlio de Mesquita Filho” (UNESP), Rua Dr. Antônio Celso Wagner Zanin, 250, Distrito de Rubião Junior, Botucatu 18618-689, SP, Brazil
| | - Ana Beatriz Bortolozo de Oliveira
- Laboratório de Biologia Celular, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista “Júlio de Mesquita Filho” (UNESP), Rua Cristóvão Colombo 2265, São José do Rio Preto 15054-000, SP, Brazil
| | - Luísa Martins Sensato Azevedo
- Laboratório de Biologia Celular, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista “Júlio de Mesquita Filho” (UNESP), Rua Cristóvão Colombo 2265, São José do Rio Preto 15054-000, SP, Brazil
| | - Cleber Galvão
- Laboratório Nacional e Internacional de Referência em Taxonomia de Triatomíneos, Instituto Oswaldo Cruz (FIOCRUZ), Av. Brazil 4365, Pavilhão Rocha Lima, Sala 505, Rio de Janeiro 21040-360, RJ, Brazil
| | - Maria Tercília Vilela de Azeredo-Oliveira
- Laboratório de Biologia Celular, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista “Júlio de Mesquita Filho” (UNESP), Rua Cristóvão Colombo 2265, São José do Rio Preto 15054-000, SP, Brazil
| | - João Aristeu da Rosa
- Laboratório de Parasitologia, Faculdade de Ciências Farmacêuticas, Universidade Estadual Paulista “Júlio de Mesquita Filho” (UNESP), Rodovia Araraquara-Jaú km 1, Araraquara 14801-902, SP, Brazil
| | - Kaio Cesar Chaboli Alevi
- Instituto de Biociências de Botucatu, Universidade Estadual Paulista “Júlio de Mesquita Filho” (UNESP), Rua Dr. Antônio Celso Wagner Zanin, 250, Distrito de Rubião Junior, Botucatu 18618-689, SP, Brazil
- Laboratório de Entomologia em Saúde Pública, Faculdade de Saúde Pública, Universidade de São Paulo (USP), Av. Dr. Arnaldo 715, São Paulo 01246-904, SP, Brazil
- Laboratório Nacional e Internacional de Referência em Taxonomia de Triatomíneos, Instituto Oswaldo Cruz (FIOCRUZ), Av. Brazil 4365, Pavilhão Rocha Lima, Sala 505, Rio de Janeiro 21040-360, RJ, Brazil
| |
Collapse
|
48
|
Yurchenko AA, Naumenko AN, Artemov GN, Karagodin DA, Hodge JM, Velichevskaya AI, Kokhanenko AA, Bondarenko SM, Abai MR, Kamali M, Gordeev MI, Moskaev AV, Caputo B, Aghayan SA, Baricheva EM, Stegniy VN, Sharakhova MV, Sharakhov IV. Phylogenomics revealed migration routes and adaptive radiation timing of Holarctic malaria mosquito species of the Maculipennis Group. BMC Biol 2023; 21:63. [PMID: 37032389 PMCID: PMC10084679 DOI: 10.1186/s12915-023-01538-w] [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: 08/10/2022] [Accepted: 02/08/2023] [Indexed: 04/11/2023] Open
Abstract
BACKGROUND Phylogenetic analyses of closely related species of mosquitoes are important for better understanding the evolution of traits contributing to transmission of vector-borne diseases. Six out of 41 dominant malaria vectors of the genus Anopheles in the world belong to the Maculipennis Group, which is subdivided into two Nearctic subgroups (Freeborni and Quadrimaculatus) and one Palearctic (Maculipennis) subgroup. Although previous studies considered the Nearctic subgroups as ancestral, details about their relationship with the Palearctic subgroup, and their migration times and routes from North America to Eurasia remain controversial. The Palearctic species An. beklemishevi is currently included in the Nearctic Quadrimaculatus subgroup adding to the uncertainties in mosquito systematics. RESULTS To reconstruct historic relationships in the Maculipennis Group, we conducted a phylogenomic analysis of 11 Palearctic and 2 Nearctic species based on sequences of 1271 orthologous genes. The analysis indicated that the Palearctic species An. beklemishevi clusters together with other Eurasian species and represents a basal lineage among them. Also, An. beklemishevi is related more closely to An. freeborni, which inhabits the Western United States, rather than to An. quadrimaculatus, a species from the Eastern United States. The time-calibrated tree suggests a migration of mosquitoes in the Maculipennis Group from North America to Eurasia about 20-25 million years ago through the Bering Land Bridge. A Hybridcheck analysis demonstrated highly significant signatures of introgression events between allopatric species An. labranchiae and An. beklemishevi. The analysis also identified ancestral introgression events between An. sacharovi and its Nearctic relative An. freeborni despite their current geographic isolation. The reconstructed phylogeny suggests that vector competence and the ability to enter complete diapause during winter evolved independently in different lineages of the Maculipennis Group. CONCLUSIONS Our phylogenomic analyses reveal migration routes and adaptive radiation timing of Holarctic malaria vectors and strongly support the inclusion of An. beklemishevi into the Maculipennis Subgroup. Detailed knowledge of the evolutionary history of the Maculipennis Subgroup provides a framework for examining the genomic changes related to ecological adaptation and susceptibility to human pathogens. These genomic variations may inform researchers about similar changes in the future providing insights into the patterns of disease transmission in Eurasia.
Collapse
Affiliation(s)
- Andrey A Yurchenko
- Department of Entomology, the Fralin Life Sciences Institute, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA
- Kurchatov Genomics Center, the Federal Research Center, Institute of Cytology and Genetics, Novosibirsk, Russia
- Current Address: INSERM U981, Gustave Roussy Institute, Université Paris-Saclay, Villejuif, France
| | - Anastasia N Naumenko
- Department of Entomology, the Fralin Life Sciences Institute, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA
| | - Gleb N Artemov
- Department of Genetics and Cell Biology and the Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, Tomsk, Russia
| | - Dmitry A Karagodin
- Laboratory of Cell Differentiation Mechanisms, the Federal Research Center, Institute of Cytology and Genetics, Novosibirsk, Russia
| | - James M Hodge
- Department of Entomology, the Fralin Life Sciences Institute, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA
| | - Alena I Velichevskaya
- Department of Genetics and Cell Biology and the Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, Tomsk, Russia
| | - Alina A Kokhanenko
- Department of Genetics and Cell Biology and the Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, Tomsk, Russia
| | - Semen M Bondarenko
- Department of Entomology, the Fralin Life Sciences Institute, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA
- Department of Genetics and Cell Biology and the Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, Tomsk, Russia
| | - Mohammad R Abai
- Department of Medical Entomology and Vector Control, Tehran University of Medical Sciences, Tehran, Iran
| | - Maryam Kamali
- Department of Medical Entomology and Parasitology, Tarbiat Modares University, Tehran, Iran
| | - Mikhail I Gordeev
- Department of General Biology and Ecology, State University of Education, Mytishchi, Russia
| | - Anton V Moskaev
- Department of General Biology and Ecology, State University of Education, Mytishchi, Russia
| | - Beniamino Caputo
- Dipartimento Di Sanità Pubblica E Malattie Infettive, Università Sapienza, Rome, Italy
| | - Sargis A Aghayan
- Scientific Center of Zoology and Hydroecology, National Academy of Sciences of the Republic of Armenia, Yerevan, Armenia
- Department of Zoology, Yerevan State University, Yerevan, Armenia
| | - Elina M Baricheva
- Laboratory of Cell Differentiation Mechanisms, the Federal Research Center, Institute of Cytology and Genetics, Novosibirsk, Russia
| | - Vladimir N Stegniy
- Department of Genetics and Cell Biology and the Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, Tomsk, Russia
| | - Maria V Sharakhova
- Department of Entomology, the Fralin Life Sciences Institute, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA.
- Laboratory of Cell Differentiation Mechanisms, the Federal Research Center, Institute of Cytology and Genetics, Novosibirsk, Russia.
| | - Igor V Sharakhov
- Department of Entomology, the Fralin Life Sciences Institute, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA.
- Department of Genetics and Cell Biology and the Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, Tomsk, Russia.
| |
Collapse
|
49
|
Nguyen VT, Collier TC, Seok S, Wang X, Mburu MM, Simubali L, Gebhardt ME, Norris DE, Lee Y. The first genome sequence of Anopheles squamous from Macha, Zambia. F1000Res 2023; 12:330. [PMID: 37842340 PMCID: PMC10568211 DOI: 10.12688/f1000research.130734.1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 03/10/2023] [Indexed: 10/17/2023] Open
Abstract
Despite efforts to minimize the impacts of malaria and reduce the number of primary vectors, malaria has yet to be eliminated in Zambia. Understudied vector species may perpetuate malaria transmission in pre-elimination settings. Anopheles squamosus is one of the most abundantly caught mosquito species in southern Zambia and has previously been found with Plasmodium falciparum sporozoites, a causal agent of human malaria. This species may be a critical vector of malaria transmission, however, there is a lack of genetic information available for An. squamosus. We report the first genome data and the first complete mitogenome (Mt) sequence of An. squamosus. The sequence was extracted from one individual mosquito from the Chidakwa area in Macha, Zambia. The raw reads were obtained using Illumina Novaseq 6000 and assembled through NOVOplasty alignment with related species. The length of the An. squamosus Mt was 15,351 bp, with 77.9 % AT content. The closest match to the whole mitochondrial genome in the phylogenetic tree is the African malaria mosquito, Anopheles gambiae. Its genome data is available through National Center for Biotechnology Information (NCBI) Sequencing Reads Archive (SRA) with accession number SRR22114392. The mitochondrial genome was deposited in NCBI GenBank with the accession number OP776919. The ITS2 containing contig sequence was deposited in GenBank with the accession number OQ241725. Mitogenome annotation and a phylogenetic tree with related Anopheles mosquito species are provided.
Collapse
Affiliation(s)
- Valerie T. Nguyen
- Florida Medical Entomology Laboratory, Department of Entomology and Nematology, Institute of Food and Agricultural Sciences, University of Florida, Gainesville, Florida, 32962, USA
| | | | - Sangwoo Seok
- Florida Medical Entomology Laboratory, Department of Entomology and Nematology, Institute of Food and Agricultural Sciences, University of Florida, Gainesville, Florida, 32962, USA
| | - Xiaodi Wang
- Florida Medical Entomology Laboratory, Department of Entomology and Nematology, Institute of Food and Agricultural Sciences, University of Florida, Gainesville, Florida, 32962, USA
| | | | | | - Mary E. Gebhardt
- The W. Harry Feinstone Department of Molecular Microbiology and Immunology, The Johns Hopkins Malaria Research Institute, Johns Hopkins Bloomberg School of Public Health, Johns Hopkins University, Baltimore, Maryland, 21205, USA
| | - Douglas E. Norris
- The W. Harry Feinstone Department of Molecular Microbiology and Immunology, The Johns Hopkins Malaria Research Institute, Johns Hopkins Bloomberg School of Public Health, Johns Hopkins University, Baltimore, Maryland, 21205, USA
| | - Yoosook Lee
- Florida Medical Entomology Laboratory, Department of Entomology and Nematology, Institute of Food and Agricultural Sciences, University of Florida, Gainesville, Florida, 32962, USA
| |
Collapse
|
50
|
Small ST, Costantini C, Sagnon N, Guelbeogo MW, Emrich SJ, Kern AD, Fontaine MC, Besansky NJ. Standing genetic variation and chromosome differences drove rapid ecotype formation in a major malaria mosquito. Proc Natl Acad Sci U S A 2023; 120:e2219835120. [PMID: 36881629 PMCID: PMC10089221 DOI: 10.1073/pnas.2219835120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2022] [Accepted: 02/09/2023] [Indexed: 03/08/2023] Open
Abstract
Species distributed across heterogeneous environments often evolve locally adapted ecotypes, but understanding of the genetic mechanisms involved in their formation and maintenance in the face of gene flow is incomplete. In Burkina Faso, the major African malaria mosquito Anopheles funestus comprises two strictly sympatric and morphologically indistinguishable yet karyotypically differentiated forms reported to differ in ecology and behavior. However, knowledge of the genetic basis and environmental determinants of An. funestus diversification was impeded by lack of modern genomic resources. Here, we applied deep whole-genome sequencing and analysis to test the hypothesis that these two forms are ecotypes differentially adapted to breeding in natural swamps versus irrigated rice fields. We demonstrate genome-wide differentiation despite extensive microsympatry, synchronicity, and ongoing hybridization. Demographic inference supports a split only ~1,300 y ago, closely following the massive expansion of domesticated African rice cultivation ~1,850 y ago. Regions of highest divergence, concentrated in chromosomal inversions, were under selection during lineage splitting, consistent with local adaptation. The origin of nearly all variations implicated in adaptation, including chromosomal inversions, substantially predates the ecotype split, suggesting that rapid adaptation was fueled mainly by standing genetic variation. Sharp inversion frequency differences likely facilitated adaptive divergence between ecotypes by suppressing recombination between opposing chromosomal orientations of the two ecotypes, while permitting free recombination within the structurally monomorphic rice ecotype. Our results align with growing evidence from diverse taxa that rapid ecological diversification can arise from evolutionarily old structural genetic variants that modify genetic recombination.
Collapse
Affiliation(s)
- Scott T. Small
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN46556
- Eck Institute for Global Health, University of Notre Dame, Notre Dame, IN46556
- Institute for Ecology and Evolution, University of Oregon, Eugene, OR97403
| | - Carlo Costantini
- Centre National de Recherche et Formation sur le Paludisme, Ouagadougou01 BP 2208, Burkina Faso
- Infectious Diseases and Vectors: Ecology, Genetics, Evolution and Control (MIVEGEC), Université de Montpellier, CNRS 5290, Institute of Research for Development (IRD) 224, F-34394Montpellier, France
| | - N’Fale Sagnon
- Centre National de Recherche et Formation sur le Paludisme, Ouagadougou01 BP 2208, Burkina Faso
| | - Moussa W. Guelbeogo
- Centre National de Recherche et Formation sur le Paludisme, Ouagadougou01 BP 2208, Burkina Faso
| | - Scott J. Emrich
- Eck Institute for Global Health, University of Notre Dame, Notre Dame, IN46556
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, IN46556
| | - Andrew D. Kern
- Institute for Ecology and Evolution, University of Oregon, Eugene, OR97403
| | - Michael C. Fontaine
- Infectious Diseases and Vectors: Ecology, Genetics, Evolution and Control (MIVEGEC), Université de Montpellier, CNRS 5290, Institute of Research for Development (IRD) 224, F-34394Montpellier, France
- Groningen Institute for Evolutionary Life Sciences, University of Groningen, 9747 AGGroningen, The Netherlands
| | - Nora J. Besansky
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN46556
- Eck Institute for Global Health, University of Notre Dame, Notre Dame, IN46556
| |
Collapse
|