201
|
Adams RH, Schield DR, Castoe TA. Recent Advances in the Inference of Gene Flow from Population Genomic Data. ACTA ACUST UNITED AC 2019. [DOI: 10.1007/s40610-019-00120-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
202
|
Pfeifer B, Kapan DD. Estimates of introgression as a function of pairwise distances. BMC Bioinformatics 2019; 20:207. [PMID: 31014244 PMCID: PMC6480520 DOI: 10.1186/s12859-019-2747-z] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Accepted: 03/18/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Research over the last 10 years highlights the increasing importance of hybridization between species as a major force structuring the evolution of genomes and potentially providing raw material for adaptation by natural and/or sexual selection. Fueled by research in a few model systems where phenotypic hybrids are easily identified, research into hybridization and introgression (the flow of genes between species) has exploded with the advent of whole-genome sequencing and emerging methods to detect the signature of hybridization at the whole-genome or chromosome level. Amongst these are a general class of methods that utilize patterns of single-nucleotide polymorphisms (SNPs) across a tree as markers of hybridization. These methods have been applied to a variety of genomic systems ranging from butterflies to Neanderthals to detect introgression, however, when employed at a fine genomic scale these methods do not perform well to quantify introgression in small sample windows. RESULTS We introduce a novel method to detect introgression by combining two widely used statistics: pairwise nucleotide diversity dxy and Patterson's D. The resulting statistic, the distance fraction (df), accounts for genetic distance across possible topologies and is designed to simultaneously detect and quantify introgression. We also relate our new method to the recently published fd and incorporate these statistics into the powerful genomics R-package PopGenome, freely available on GitHub (pievos101/PopGenome) and the Comprehensive R Archive Network (CRAN). The supplemental material contains a wide range of simulation studies and a detailed manual how to perform the statistics within the PopGenome framework. CONCLUSION We present a new distance based statistic df that avoids the pitfalls of Patterson's D when applied to small genomic regions and accurately quantifies the fraction of introgression (f) for a wide range of simulation scenarios.
Collapse
Affiliation(s)
- Bastian Pfeifer
- Institute for Medical Informatics, Statistics and Documentation, Medical University, Graz, Austria
| | - Durrell D. Kapan
- Department of Entomology and Center for Comparative Genomics, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, 55 Music Concourse Dr., San Francisco, USA
| |
Collapse
|
203
|
Pfeifer B, Kapan DD. Estimates of introgression as a function of pairwise distances. BMC Bioinformatics 2019; 20:207. [PMID: 31014244 DOI: 10.1101/154377] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Accepted: 03/18/2019] [Indexed: 05/27/2023] Open
Abstract
BACKGROUND Research over the last 10 years highlights the increasing importance of hybridization between species as a major force structuring the evolution of genomes and potentially providing raw material for adaptation by natural and/or sexual selection. Fueled by research in a few model systems where phenotypic hybrids are easily identified, research into hybridization and introgression (the flow of genes between species) has exploded with the advent of whole-genome sequencing and emerging methods to detect the signature of hybridization at the whole-genome or chromosome level. Amongst these are a general class of methods that utilize patterns of single-nucleotide polymorphisms (SNPs) across a tree as markers of hybridization. These methods have been applied to a variety of genomic systems ranging from butterflies to Neanderthals to detect introgression, however, when employed at a fine genomic scale these methods do not perform well to quantify introgression in small sample windows. RESULTS We introduce a novel method to detect introgression by combining two widely used statistics: pairwise nucleotide diversity dxy and Patterson's D. The resulting statistic, the distance fraction (df), accounts for genetic distance across possible topologies and is designed to simultaneously detect and quantify introgression. We also relate our new method to the recently published fd and incorporate these statistics into the powerful genomics R-package PopGenome, freely available on GitHub (pievos101/PopGenome) and the Comprehensive R Archive Network (CRAN). The supplemental material contains a wide range of simulation studies and a detailed manual how to perform the statistics within the PopGenome framework. CONCLUSION We present a new distance based statistic df that avoids the pitfalls of Patterson's D when applied to small genomic regions and accurately quantifies the fraction of introgression (f) for a wide range of simulation scenarios.
Collapse
Affiliation(s)
- Bastian Pfeifer
- Institute for Medical Informatics, Statistics and Documentation, Medical University, Graz, Austria
| | - Durrell D Kapan
- Department of Entomology and Center for Comparative Genomics, Institute for Biodiversity Science and Sustainability, California Academy of Sciences, 55 Music Concourse Dr., San Francisco, USA.
| |
Collapse
|
204
|
Carter TE, Yared S, Hansel S, Lopez K, Janies D. Sequence-based identification of Anopheles species in eastern Ethiopia. Malar J 2019; 18:135. [PMID: 30992003 PMCID: PMC6469081 DOI: 10.1186/s12936-019-2768-0] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2019] [Accepted: 04/04/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND The recent finding of a typically non-African Anopheles species in eastern Ethiopia emphasizes the need for detailed species identification and characterization for effective malaria vector surveillance. Molecular approaches increase the accuracy and interoperability of vector surveillance data. To develop effective molecular assays for Anopheles identification, it is important to evaluate different genetic loci for the ability to characterize species and population level variation. Here the utility of the internal transcribed spacer 2 (ITS2) and cytochrome oxidase I (COI) loci for detection of Anopheles species from understudied regions of eastern Ethiopia was investigated. METHODS Adult mosquitoes were collected from the Harewe locality (east) and Meki (east central) Ethiopia. PCR and Sanger sequencing were performed for portions of the ITS2 and COI loci. Both NCBI's Basic Local Alignment Search tool (BLAST) and phylogenetic analysis using a maximum-likelihood approach were performed to identify species of Anopheles specimens. RESULTS Two species from the east Ethiopian collection, Anopheles arabiensis and Anopheles pretoriensis were identified. Analyses of ITS2 locus resulted in delineation of both species. In contrast, analysis of COI locus could not be used to delineate An. arabiensis from other taxa in Anopheles gambiae complex, but could distinguish An. pretoriensis sequences from sister taxa. CONCLUSION The lack of clarity from COI sequence analysis highlights potential challenges of species identification within species complexes. These results provide supporting data for the development of molecular assays for delineation of Anopheles in east Ethiopia.
Collapse
Affiliation(s)
- Tamar E Carter
- Department of Biology, Baylor University, Waco, TX, USA.
| | - Solomon Yared
- Department of Biology, Jigjiga University, Jigjiga, Ethiopia
| | - Shantoy Hansel
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC, USA
| | - Karen Lopez
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC, USA
| | - Daniel Janies
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC, USA
| |
Collapse
|
205
|
Feng S, Ru D, Sun Y, Mao K, Milne R, Liu J. Trans-lineage polymorphism and nonbifurcating diversification of the genus Picea. THE NEW PHYTOLOGIST 2019; 222:576-587. [PMID: 30415488 DOI: 10.1111/nph.15590] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/19/2018] [Accepted: 11/02/2018] [Indexed: 06/09/2023]
Abstract
Nonbifurcating divergence caused by introgressive hybridization is continuously reported for groups of closely related species. In this study, we aimed to reconstruct the genome-scale classification of deep lineages of the conifer genus Picea, establish their phylogenetic relationships and test the bifurcating hypothesis between deeply branching lineages based on genomic data. We sequenced the transcriptomes of 35 individuals of 27 taxa covering all main lineages of the genus. Four major lineages, comprising three to 12 taxa each, largely consistent with morphological evidence, were recovered across the coalescent and integrated nuclear phylogeny. However, many of the individual gene trees recovered contradict one another. Moreover, the well-supported coalescent tree inferred here differs from previous studies based on various DNA markers, with respect to topology and inter-lineage relationships. We identified the shared polymorphisms between four major lineages. ABBA-BABA tests confirmed the inter-lineage gene flow and thus violated the bifurcating divergence model. Gene flow occurred more frequently between lineages distributed in the same continent than those disjunct between continents. Our results indicate that introgression and nonbifurcating diversification apply, even between deeply branching lineages of the conifer genus Picea.
Collapse
Affiliation(s)
- Shuo Feng
- State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology & College of Life Sciences, Lanzhou University, Lanzhou, 730000, China
| | - Dafu Ru
- Key Laboratory for Bio-resource and Eco-environment of Ministry of Education, College of Life Science, Sichuan University, Chengdu, 610064, China
| | - Yongshuai Sun
- Key Laboratory for Bio-resource and Eco-environment of Ministry of Education, College of Life Science, Sichuan University, Chengdu, 610064, China
| | - Kangshan Mao
- Key Laboratory for Bio-resource and Eco-environment of Ministry of Education, College of Life Science, Sichuan University, Chengdu, 610064, China
| | - Richard Milne
- Institute of Molecular Plant Sciences, School of Biological Sciences, The University of Edinburgh, Daniel Rutherford Building, The King's Buildings, Mayfield Road, Edinburgh, EH9 3JH, UK
- Royal Botanic Garden Edinburgh, 20A Inverleith Row, Edinburgh, EH3 5LR, UK
| | - Jianquan Liu
- State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology & College of Life Sciences, Lanzhou University, Lanzhou, 730000, China
- Key Laboratory for Bio-resource and Eco-environment of Ministry of Education, College of Life Science, Sichuan University, Chengdu, 610064, China
| |
Collapse
|
206
|
Galtier N. Delineating species in the speciation continuum: A proposal. Evol Appl 2019; 12:657-663. [PMID: 30976300 PMCID: PMC6439491 DOI: 10.1111/eva.12748] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2018] [Accepted: 11/08/2018] [Indexed: 12/18/2022] Open
Abstract
Delineating species is a difficult and seemingly uninteresting issue that is still essential to address. Taxonomic methodology is heterogeneous according to the taxa and scientists involved due to the disparate data quality and quantity and disagreements over the species concept. This has negative impacts on basic and applied research. Genomic data substantially enhance our understanding of the speciation process but do not provide a ubiquitous solution to the species problem. The relevance of comparative approaches in speciation research has nevertheless recently been demonstrated. I suggest moving towards a more unified taxonomic classification through a reference-based decision procedure.
Collapse
Affiliation(s)
- Nicolas Galtier
- UMR 5554 – Institut des Sciences de l'EvolutionCNRS – University of Montpellier – IRD – EPHEMontpellierFrance
| |
Collapse
|
207
|
Abstract
In this perspective, we evaluate the explanatory power of the neutral theory of molecular evolution, 50 years after its introduction by Kimura. We argue that the neutral theory was supported by unreliable theoretical and empirical evidence from the beginning, and that in light of modern, genome-scale data, we can firmly reject its universality. The ubiquity of adaptive variation both within and between species means that a more comprehensive theory of molecular evolution must be sought.
Collapse
Affiliation(s)
- Andrew D Kern
- Department of Genetics, Rutgers University, Piscataway, NJ
| | - Matthew W Hahn
- Department of Biology and Department of Computer Science, Indiana University Bloomington, IN
| |
Collapse
|
208
|
Wang Y, Zhou X, Wang L, Liu X, Yang D, Rokas A. Gene Selection and Evolutionary Modeling Affect Phylogenomic Inference of Neuropterida Based on Transcriptome Data. Int J Mol Sci 2019; 20:E1072. [PMID: 30832228 PMCID: PMC6429444 DOI: 10.3390/ijms20051072] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2018] [Revised: 02/18/2019] [Accepted: 02/26/2019] [Indexed: 11/30/2022] Open
Abstract
Neuropterida is a super order of Holometabola that consists of the orders Megaloptera (dobsonflies, fishflies, and alderflies), Neuroptera (lacewings) and Raphidioptera (snakeflies). Several proposed higher-level relationships within Neuropterida, such as the relationships between the orders or between the families, have been extensively debated. To further understand the evolutionary history of Neuropterida, we conducted phylogenomic analyses of all 13 published transcriptomes of the neuropterid species, as well as of a new transcriptome of the fishfly species Ctenochauliodes similis of Liu and Yang, 2006 (Megaloptera: Corydalidae: Chauliodinae) that we sequenced. Our phylogenomic data matrix contained 1392 ortholog genes from 22 holometabolan species representing six families from Neuroptera, two families from Raphidioptera, and two families from Megaloptera as the ingroup taxa, and nine orders of Holometabola as outgroups. Phylogenetic reconstruction was performed using both concatenation and coalescent-based approaches under a site-homogeneous model as well as under a site-heterogeneous model. Surprisingly, analyses using the site-homogeneous model strongly supported a paraphyletic Neuroptera, with Coniopterygidae assigned as the sister group of all other Neuropterida. In contrast, analyses using the site-heterogeneous model recovered Neuroptera as monophyletic. The monophyly of Neuroptera was also recovered in concatenation and coalescent-based analyses using genes with stronger phylogenetic signals [i.e., higher average bootstrap support (ABS) values and higher relative tree certainty including all conflicting bipartitions (RTCA) values] under the site-homogeneous model. The present study illustrated how both data selection and model selection influence phylogenomic analyses of large-scale data matrices comprehensively.
Collapse
Affiliation(s)
- Yuyu Wang
- College of Plant Protection, Hebei Agricultural University, Baoding 071001, China.
- Department of Entomology, China Agricultural University, Beijing 100193, China.
- Department of Biological Sciences, Vanderbilt University, Nashville, TN 37235, USA.
| | - Xiaofan Zhou
- Department of Biological Sciences, Vanderbilt University, Nashville, TN 37235, USA.
- Guangdong Province Key Laboratory of Microbial Signals and Disease Control, Integrative Microbiology Research Centre, South China Agricultural University, Guangzhou 510642, China.
| | - Liming Wang
- College of Plant Protection, Hebei Agricultural University, Baoding 071001, China.
| | - Xingyue Liu
- Department of Entomology, China Agricultural University, Beijing 100193, China.
| | - Ding Yang
- Department of Entomology, China Agricultural University, Beijing 100193, China.
| | - Antonis Rokas
- Department of Biological Sciences, Vanderbilt University, Nashville, TN 37235, USA.
| |
Collapse
|
209
|
Oakeshott JG, Robin C, Gordon KH. Editorial overview: Revisiting Dobzhansky and the 'modern synthesis' in light of insect evolutionary genomics. CURRENT OPINION IN INSECT SCIENCE 2019; 31:iii-vi. [PMID: 31109682 DOI: 10.1016/j.cois.2019.04.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Affiliation(s)
- John G Oakeshott
- CSIRO, Clunies Ross St, GPO Box 1700, Acton, ACT, 2601, Australia.
| | - Charles Robin
- School of Biosciences, University of Melbourne, Grattan St, Melbourne, VIC, 3010, Australia
| | - Karl Hj Gordon
- CSIRO, Clunies Ross St, GPO Box 1700, Acton, ACT, 2601, Australia
| |
Collapse
|
210
|
Martin SH, Davey JW, Salazar C, Jiggins CD. Recombination rate variation shapes barriers to introgression across butterfly genomes. PLoS Biol 2019; 17:e2006288. [PMID: 30730876 PMCID: PMC6366726 DOI: 10.1371/journal.pbio.2006288] [Citation(s) in RCA: 165] [Impact Index Per Article: 33.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2018] [Accepted: 01/07/2019] [Indexed: 12/30/2022] Open
Abstract
Hybridisation and introgression can dramatically alter the relationships among groups of species, leading to phylogenetic discordance across the genome and between populations. Introgression can also erode species differences over time, but selection against introgression at certain loci acts to maintain postmating species barriers. Theory predicts that species barriers made up of many loci throughout the genome should lead to a broad correlation between introgression and recombination rate, which determines the extent to which selection on deleterious foreign alleles will affect neutral alleles at physically linked loci. Here, we describe the variation in genealogical relationships across the genome among three species of Heliconius butterflies: H. melpomene (mel), H. cydno (cyd), and H. timareta (tim), using whole genomes of 92 individuals, and ask whether this variation can be explained by heterogeneous barriers to introgression. We find that species relationships vary predictably at the chromosomal scale. By quantifying recombination rate and admixture proportions, we then show that rates of introgression are predicted by variation in recombination rate. This implies that species barriers are highly polygenic, with selection acting against introgressed alleles across most of the genome. In addition, long chromosomes, which have lower recombination rates, produce stronger barriers on average than short chromosomes. Finally, we find a consistent difference between two species pairs on either side of the Andes, which suggests differences in the architecture of the species barriers. Our findings illustrate how the combined effects of hybridisation, recombination, and natural selection, acting at multitudes of loci over long periods, can dramatically sculpt the phylogenetic relationships among species.
Collapse
Affiliation(s)
- Simon H. Martin
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - John W. Davey
- Department of Biology, University of York, York, United Kingdom
| | - Camilo Salazar
- Biology Program, Faculty of Natural Sciences and Mathematics, Universidad del Rosario, Bogota, Colombia
| | - Chris D. Jiggins
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
211
|
Wu M, Kostyun JL, Moyle LC. Genome Sequence of Jaltomata Addresses Rapid Reproductive Trait Evolution and Enhances Comparative Genomics in the Hyper-Diverse Solanaceae. Genome Biol Evol 2019; 11:335-349. [PMID: 30608583 PMCID: PMC6368146 DOI: 10.1093/gbe/evy274] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/16/2018] [Indexed: 12/11/2022] Open
Abstract
Within the economically important plant family Solanaceae, Jaltomata is a rapidly evolving genus that has extensive diversity in flower size and shape, as well as fruit and nectar color, among its ∼80 species. Here, we report the whole-genome sequencing, assembly, and annotation, of one representative species (Jaltomata sinuosa) from this genus. Combining PacBio long reads (25×) and Illumina short reads (148×) achieved an assembly of ∼1.45 Gb, spanning ∼96% of the estimated genome. Ninety-six percent of curated single-copy orthologs in plants were detected in the assembly, supporting a high level of completeness of the genome. Similar to other Solanaceous species, repetitive elements made up a large fraction (∼80%) of the genome, with the most recently active element, Gypsy, expanding across the genome in the last 1–2 Myr. Computational gene prediction, in conjunction with a merged transcriptome data set from 11 tissues, identified 34,725 protein-coding genes. Comparative phylogenetic analyses with six other sequenced Solanaceae species determined that Jaltomata is most likely sister to Solanum, although a large fraction of gene trees supported a conflicting bipartition consistent with substantial introgression between Jaltomata and Capsicum after these species split. We also identified gene family dynamics specific to Jaltomata, including expansion of gene families potentially involved in novel reproductive trait development, and loss of gene families that accompanied the loss of self-incompatibility. This high-quality genome will facilitate studies of phenotypic diversification in this rapidly radiating group and provide a new point of comparison for broader analyses of genomic evolution across the Solanaceae.
Collapse
Affiliation(s)
- Meng Wu
- Department of Biology, Indiana University Bloomington
| | - Jamie L Kostyun
- Department of Biology, Indiana University Bloomington.,Department of Plant Biology, University of Vermont
| | | |
Collapse
|
212
|
Flagel L, Brandvain Y, Schrider DR. The Unreasonable Effectiveness of Convolutional Neural Networks in Population Genetic Inference. Mol Biol Evol 2019; 36:220-238. [PMID: 30517664 PMCID: PMC6367976 DOI: 10.1093/molbev/msy224] [Citation(s) in RCA: 95] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Population-scale genomic data sets have given researchers incredible amounts of information from which to infer evolutionary histories. Concomitant with this flood of data, theoretical and methodological advances have sought to extract information from genomic sequences to infer demographic events such as population size changes and gene flow among closely related populations/species, construct recombination maps, and uncover loci underlying recent adaptation. To date, most methods make use of only one or a few summaries of the input sequences and therefore ignore potentially useful information encoded in the data. The most sophisticated of these approaches involve likelihood calculations, which require theoretical advances for each new problem, and often focus on a single aspect of the data (e.g., only allele frequency information) in the interest of mathematical and computational tractability. Directly interrogating the entirety of the input sequence data in a likelihood-free manner would thus offer a fruitful alternative. Here, we accomplish this by representing DNA sequence alignments as images and using a class of deep learning methods called convolutional neural networks (CNNs) to make population genetic inferences from these images. We apply CNNs to a number of evolutionary questions and find that they frequently match or exceed the accuracy of current methods. Importantly, we show that CNNs perform accurate evolutionary model selection and parameter estimation, even on problems that have not received detailed theoretical treatments. Thus, when applied to population genetic alignments, CNNs are capable of outperforming expert-derived statistical methods and offer a new path forward in cases where no likelihood approach exists.
Collapse
Affiliation(s)
- Lex Flagel
- Monsanto Company, Chesterfield, MO
- Department of Plant and Microbial Biology, University of Minnesota, St. Paul, MN
| | - Yaniv Brandvain
- Department of Plant and Microbial Biology, University of Minnesota, St. Paul, MN
| | - Daniel R Schrider
- Department of Genetics, University of North Carolina, Chapel Hill, NC
| |
Collapse
|
213
|
Doellman MM, Feder JL. Genomic transitions during host race and species formation. CURRENT OPINION IN INSECT SCIENCE 2019; 31:84-92. [PMID: 31109679 DOI: 10.1016/j.cois.2018.11.006] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2018] [Revised: 11/19/2018] [Accepted: 11/20/2018] [Indexed: 06/09/2023]
Abstract
Darwin recognized species as discontinuous, yet considered them to be formed by an incremental process of natural selection. Recent theoretical work on 'genome-wide congealing' is bridging this gap between the gradualism of divergent selection and rapid genome-wide divergence, particularly during ecological speciation-with-gene-flow. Host races and species of phytophagous insects, displaying a spectrum of divergence and gene flow among member taxa, provide model systems for testing predicted non-linear transitions from 'genic' divergence at a few uncoupled loci to 'genomic' divergence with genome-wide coupling of selected loci and strong reproductive isolation. Integrating across natural history, genomics, and evolutionary theory, emerging research suggests a tipping point from 'genic' to 'genomic' divergence between host races and species, during both sympatric speciation and secondary contact.
Collapse
Affiliation(s)
- Meredith M Doellman
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN 46556, USA.
| | - Jeffrey L Feder
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN 46556, USA
| |
Collapse
|
214
|
Janzen GM, Wang L, Hufford MB. The extent of adaptive wild introgression in crops. THE NEW PHYTOLOGIST 2019; 221:1279-1288. [PMID: 30368812 DOI: 10.1111/nph.15457] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/08/2018] [Accepted: 08/24/2018] [Indexed: 05/05/2023]
Abstract
The study of crop evolution has focused primarily on the process of initial domestication. Post-domestication adaptation during the expansion of crops from their centers of origin has received considerably less attention. Recent research has revealed that, in at least some instances, crops have received introgression from their wild relatives that has facilitated adaptation to novel conditions encountered during expansion. Such adaptive introgression could have an important impact on the basic study of domestication, affecting estimates of several evolutionary processes of interest (e.g. the strength of the domestication bottleneck, the timing of domestication, the targets of selection during domestication). Identification of haplotypes introgressed from the wild may also help in the identification of alleles that are beneficial under particular environmental conditions. Here we review mounting evidence for substantial adaptive wild introgression in several crops and consider the implications of such gene flow to our understanding of crop histories.
Collapse
Affiliation(s)
- Garrett M Janzen
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, 50011, USA
| | - Li Wang
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, 50011, USA
| | - Matthew B Hufford
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, 50011, USA
| |
Collapse
|
215
|
Taylor SA, Larson EL. Insights from genomes into the evolutionary importance and prevalence of hybridization in nature. Nat Ecol Evol 2019; 3:170-177. [DOI: 10.1038/s41559-018-0777-y] [Citation(s) in RCA: 211] [Impact Index Per Article: 42.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2018] [Accepted: 12/04/2018] [Indexed: 01/27/2023]
|
216
|
The Timing and Direction of Introgression Under the Multispecies Network Coalescent. Genetics 2019; 211:1059-1073. [PMID: 30670542 DOI: 10.1534/genetics.118.301831] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2018] [Accepted: 01/21/2019] [Indexed: 12/26/2022] Open
Abstract
Introgression is a pervasive biological process, and many statistical methods have been developed to infer its presence from genomic data. However, many of the consequences and genomic signatures of introgression remain unexplored from a methodological standpoint. Here, we develop a model for the timing and direction of introgression based on the multispecies network coalescent, and from it suggest new approaches for testing introgression hypotheses. We suggest two new statistics, D 1 and D 2, which can be used in conjunction with other information to test hypotheses relating to the timing and direction of introgression, respectively. D 1 may find use in evaluating cases of homoploid hybrid speciation (HHS), while D 2 provides a four-taxon test for polarizing introgression. Although analytical expectations for our statistics require a number of assumptions to be met, we show how simulations can be used to test hypotheses about introgression when these assumptions are violated. We apply the D 1 statistic to genomic data from the wild yeast Saccharomyces paradoxus-a proposed example of HHS-demonstrating its use as a test of this model. These methods provide new and powerful ways to address questions relating to the timing and direction of introgression.
Collapse
|
217
|
Izquierdo A, Fahrenberger M, Persampieri T, Benedict MQ, Giles T, Catteruccia F, Emes RD, Dottorini T. Evolution of gene expression levels in the male reproductive organs of Anopheles mosquitoes. Life Sci Alliance 2019; 2:e201800191. [PMID: 30623175 PMCID: PMC6315087 DOI: 10.26508/lsa.201800191] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2018] [Revised: 12/21/2018] [Accepted: 12/21/2018] [Indexed: 12/31/2022] Open
Abstract
Modifications in gene expression determine many of the phenotypic differentiations between closely related species. This is particularly evident in reproductive tissues, where evolution of genes is more rapid, facilitating the appearance of distinct reproductive characteristics which may lead to species isolation and phenotypic variation. Large-scale, comparative analyses of transcript expression levels have been limited until recently by lack of inter-species data mining solutions. Here, by combining expression normalisation across lineages, multivariate statistical analysis, evolutionary rate, and protein-protein interaction analysis, we investigate ortholog transcripts in the male accessory glands and testes across five closely related species in the Anopheles gambiae complex. We first demonstrate that the differentiation by transcript expression is consistent with the known Anopheles phylogeny. Then, through clustering, we discover groups of transcripts with tissue-dependent expression patterns conserved across lineages, or lineage-dependent patterns conserved across tissues. The strongest associations with reproductive function, transcriptional regulatory networks, protein-protein subnetworks, and evolutionary rate are found for the groups of transcripts featuring large expression differences in lineage or tissue-conserved patterns.
Collapse
Affiliation(s)
- Abril Izquierdo
- School of Veterinary Medicine and Science, Sutton Bonington Campus, University of Nottingham, Leicestershire, UK
| | - Martin Fahrenberger
- School of Veterinary Medicine and Science, Sutton Bonington Campus, University of Nottingham, Leicestershire, UK
| | - Tania Persampieri
- Department of Experimental Medicine, University of Perugia, Perugia, Italy
| | - Mark Q Benedict
- Centers for Disease Control and Prevention, Division of Parasitic Diseases and Malaria, Entomology Branch, Atlanta, GA, USA
| | - Tom Giles
- Advanced Data Analysis Centre, Sutton Bonington Campus, University of Nottingham, Leicestershire, UK
| | - Flaminia Catteruccia
- Department of Immunology and Infectious Diseases, Harvard T. H. Chan School of Public Health, Boston, MA, USA
| | - Richard D Emes
- School of Veterinary Medicine and Science, Sutton Bonington Campus, University of Nottingham, Leicestershire, UK.,Advanced Data Analysis Centre, Sutton Bonington Campus, University of Nottingham, Leicestershire, UK
| | - Tania Dottorini
- School of Veterinary Medicine and Science, Sutton Bonington Campus, University of Nottingham, Leicestershire, UK
| |
Collapse
|
218
|
Ruzzante L, Reijnders MJ, Waterhouse RM. Of Genes and Genomes: Mosquito Evolution and Diversity. Trends Parasitol 2019; 35:32-51. [DOI: 10.1016/j.pt.2018.10.003] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2018] [Revised: 10/07/2018] [Accepted: 10/08/2018] [Indexed: 12/16/2022]
|
219
|
Advances in Computational Methods for Phylogenetic Networks in the Presence of Hybridization. BIOINFORMATICS AND PHYLOGENETICS 2019. [DOI: 10.1007/978-3-030-10837-3_13] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
220
|
Hanemaaijer MJ, Collier TC, Chang A, Shott CC, Houston PD, Schmidt H, Main BJ, Cornel AJ, Lee Y, Lanzaro GC. The fate of genes that cross species boundaries after a major hybridization event in a natural mosquito population. Mol Ecol 2018; 27:4978-4990. [DOI: 10.1111/mec.14947] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Revised: 07/17/2018] [Accepted: 07/19/2018] [Indexed: 12/30/2022]
Affiliation(s)
- Mark J. Hanemaaijer
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology UC Davis Davis California
| | - Travis C. Collier
- Daniel K. Inouye US Pacific Basin Agricultural Research Center (PBARC), United States Department of Agriculture Agricultural Research Service Hilo Hawaii
| | - Allison Chang
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology UC Davis Davis California
| | - Chloe C. Shott
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology UC Davis Davis California
| | - Parker D. Houston
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology UC Davis Davis California
| | - Hanno Schmidt
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology UC Davis Davis California
| | - Bradley J. Main
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology UC Davis Davis California
| | - Anthony J. Cornel
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology UC Davis Davis California
- Mosquito Control Research Laboratory, Department of Entomology and Nematology University of California Parlier California
- School of Health Systems & Public Health University of Pretoria Pretoria South Africa
| | - Yoosook Lee
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology UC Davis Davis California
| | - Gregory C. Lanzaro
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology UC Davis Davis California
| |
Collapse
|
221
|
Rouard M, Droc G, Martin G, Sardos J, Hueber Y, Guignon V, Cenci A, Geigle B, Hibbins MS, Yahiaoui N, Baurens FC, Berry V, Hahn MW, D’Hont A, Roux N. Three New Genome Assemblies Support a Rapid Radiation in Musa acuminata (Wild Banana). Genome Biol Evol 2018; 10:3129-3140. [PMID: 30321324 PMCID: PMC6282646 DOI: 10.1093/gbe/evy227] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/10/2018] [Indexed: 12/15/2022] Open
Abstract
Edible bananas result from interspecific hybridization between Musa acuminata and Musa balbisiana, as well as among subspecies in M. acuminata. Four particular M. acuminata subspecies have been proposed as the main contributors of edible bananas, all of which radiated in a short period of time in southeastern Asia. Clarifying the evolution of these lineages at a whole-genome scale is therefore an important step toward understanding the domestication and diversification of this crop. This study reports the de novo genome assembly and gene annotation of a representative genotype from three different subspecies of M. acuminata. These data are combined with the previously published genome of the fourth subspecies to investigate phylogenetic relationships. Analyses of shared and unique gene families reveal that the four subspecies are quite homogenous, with a core genome representing at least 50% of all genes and very few M. acuminata species-specific gene families. Multiple alignments indicate high sequence identity between homologous single copy-genes, supporting the close relationships of these lineages. Interestingly, phylogenomic analyses demonstrate high levels of gene tree discordance, due to both incomplete lineage sorting and introgression. This pattern suggests rapid radiation within Musa acuminata subspecies that occurred after the divergence with M. balbisiana. Introgression between M. a. ssp. malaccensis and M. a. ssp. burmannica was detected across the genome, though multiple approaches to resolve the subspecies tree converged on the same topology. To support evolutionary and functional analyses, we introduce the PanMusa database, which enables researchers to exploration of individual gene families and trees.
Collapse
Affiliation(s)
- Mathieu Rouard
- Bioversity International, Parc Scientifique Agropolis II, Montpellier, France
| | - Gaetan Droc
- CIRAD, UMR AGAP, Montpellier, France
- AGAP, Univ Montpellier, CIRAD, INRA, Montpellier SupAgro, France
| | - Guillaume Martin
- CIRAD, UMR AGAP, Montpellier, France
- AGAP, Univ Montpellier, CIRAD, INRA, Montpellier SupAgro, France
| | - Julie Sardos
- Bioversity International, Parc Scientifique Agropolis II, Montpellier, France
| | - Yann Hueber
- Bioversity International, Parc Scientifique Agropolis II, Montpellier, France
| | - Valentin Guignon
- Bioversity International, Parc Scientifique Agropolis II, Montpellier, France
| | - Alberto Cenci
- Bioversity International, Parc Scientifique Agropolis II, Montpellier, France
| | | | - Mark S Hibbins
- Department of Biology, Indiana University
- Department of Computer Science, Indiana University
| | - Nabila Yahiaoui
- CIRAD, UMR AGAP, Montpellier, France
- AGAP, Univ Montpellier, CIRAD, INRA, Montpellier SupAgro, France
| | - Franc-Christophe Baurens
- CIRAD, UMR AGAP, Montpellier, France
- AGAP, Univ Montpellier, CIRAD, INRA, Montpellier SupAgro, France
| | - Vincent Berry
- LIRMM, Université de Montpellier, CNRS, Montpellier, France
| | - Matthew W Hahn
- Department of Biology, Indiana University
- Department of Computer Science, Indiana University
| | - Angelique D’Hont
- CIRAD, UMR AGAP, Montpellier, France
- AGAP, Univ Montpellier, CIRAD, INRA, Montpellier SupAgro, France
| | - Nicolas Roux
- Bioversity International, Parc Scientifique Agropolis II, Montpellier, France
| |
Collapse
|
222
|
Abstract
Convergent evolution provides key evidence for the action of natural selection. The process of convergence is often inferred because the same trait appears in multiple species that are not closely related. However, different parts of the genome can reveal different relationships among species, with some genes or regions uniting lineages that appear unrelated in the species tree. If changes in traits occur in these discordant regions, a false pattern of convergence can be produced (known as “hemiplasy”). Here, we provide a way to quantify the probability that hemiplasy occurs and contrast it with the probability of convergence. We find that hemiplasy is likely to explain many apparent cases of convergent evolution, even when the fraction of discordant regions is low. Convergent evolution—the appearance of the same character state in apparently unrelated organisms—is often inferred when a trait is incongruent with the species tree. However, trait incongruence can also arise from changes that occur on discordant gene trees, a process referred to as hemiplasy. Hemiplasy is rarely taken into account in studies of convergent evolution, despite the fact that phylogenomic studies have revealed rampant discordance. Here, we study the relative probabilities of homoplasy (including convergence and reversal) and hemiplasy for an incongruent trait. We derive expressions for the probabilities of the two events, showing that they depend on many of the same parameters. We find that hemiplasy is as likely—or more likely—than homoplasy for a wide range of conditions, even when levels of discordance are low. We also present a method to calculate the ratio of these two probabilities (the “hemiplasy risk factor”) along the branches of a phylogeny of arbitrary length. Such calculations can be applied to any tree to identify when and where incongruent traits may be due to hemiplasy.
Collapse
|
223
|
Nowling RJ, Emrich SJ. Adjusted likelihood-ratio test for variants with unknown genotypes. J Bioinform Comput Biol 2018; 16:1840020. [DOI: 10.1142/s0219720018400206] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Association tests performed with the Likelihood-Ratio Test (LR Test) can be an alternative to [Formula: see text], which is often used in population genetics to find variants of interest. Because the LR Test has several properties that could make it preferable to [Formula: see text], we propose a novel approach for modeling unknown genotypes in highly-similar species. To show the effectiveness of this LR Test approach, we apply it to single-nucleotide polymorphisms (SNPs) associated with the recent speciation of the malaria vectors Anopheles gambiae and Anopheles coluzzii and compare to [Formula: see text].
Collapse
Affiliation(s)
- Ronald J. Nowling
- Electrical Engineering and Computer Science, Milwaukee School of Engineering, 1025 North Broadway, Milwaukee, WI 53202, USA
| | - Scott J. Emrich
- Electrical Engineering and Computer Science, University of Tennessee, 1520 Middle Drive, Knoxville, TN 37996, USA
| |
Collapse
|
224
|
Deitz KC, Takken W, Slotman MA. The Effect of Hybridization on Dosage Compensation in Member Species of the Anopheles gambiae Species Complex. Genome Biol Evol 2018; 10:1663-1672. [PMID: 29860336 PMCID: PMC6037052 DOI: 10.1093/gbe/evy108] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/30/2018] [Indexed: 12/17/2022] Open
Abstract
Dosage compensation has evolved in concert with Y-chromosome degeneration in many taxa that exhibit heterogametic sex chromosomes. Dosage compensation overcomes the biological challenge of a “half dose” of X chromosome gene transcripts in the heterogametic sex. The need to equalize gene expression of a hemizygous X with that of autosomes arises from the fact that the X chromosomes retain hundreds of functional genes that are actively transcribed in both sexes and interact with genes expressed on the autosomes. Sex determination and heterogametic sex chromosomes have evolved multiple times in Diptera, and in each case the genetic control of dosage compensation is tightly linked to sex determination. In the Anopheles gambiae species complex (Culicidae), maleness is conferred by the Y-chromosome gene Yob, which despite its conserved role between species is polymorphic in its copy number between them. Previous work demonstrated that male An. gambiae s.s. males exhibit complete dosage compensation in pupal and adult stages. In the present study, we have extended this analysis to three sister species in the An. gambiae complex: An. coluzzii, An. arabiensis, and An. quadriannulatus. In addition, we analyzed dosage compensation in bi-directional F1 hybrids between these species to determine if hybridization results in the mis-regulation and disruption of dosage compensation. Our results confirm that dosage compensation operates in the An. gambiae species complex through the hypertranscription of the male X chromosome. Additionally, dosage compensation in hybrid males does not differ from parental males, indicating that hybridization does not result in the mis-regulation of dosage compensation.
Collapse
Affiliation(s)
- Kevin C Deitz
- Department of Entomology, Texas A&M University.,Department of Ecology and Evolutionary Biology and The Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ
| | - Willem Takken
- Laboratory of Entomology, Department of Plant Sciences, Wageningen University, The Netherlands
| | | |
Collapse
|
225
|
Role of sexual imprinting in assortative mating and premating isolation in Darwin's finches. Proc Natl Acad Sci U S A 2018; 115:E10879-E10887. [PMID: 30348758 DOI: 10.1073/pnas.1813662115] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Global biodiversity is being degraded at an unprecedented rate, so it is important to preserve the potential for future speciation. Providing for the future requires understanding speciation as a contemporary ecological process. Phylogenetically young adaptive radiations are a good choice for detailed study because diversification is ongoing. A key question is how incipient species become reproductively isolated from each other. Barriers to gene exchange have been investigated experimentally in the laboratory and in the field, but little information exists from the quantitative study of mating patterns in nature. Although the degree to which genetic variation underlying mate-preference learning is unknown, we provide evidence that two species of Darwin's finches imprint on morphological cues of their parents and mate assortatively. Statistical evidence of presumed imprinting is stronger for sons than for daughters and is stronger for imprinting on fathers than on mothers. In combination, morphology and species-specific song learned from the father constitute a barrier to interbreeding. The barrier becomes stronger the more the species diverge morphologically and ecologically. It occasionally breaks down, and the species hybridize. Hybridization is most likely to happen when species are similar to each other in adaptive morphological traits, e.g., body size and beak size and shape. Hybridization can lead to the formation of a new species reproductively isolated from the parental species as a result of sexual imprinting. Conservation of sufficiently diverse natural habitat is needed to sustain a large sample of extant biota and preserve the potential for future speciation.
Collapse
|
226
|
Ayala D, Zhang S, Chateau M, Fouet C, Morlais I, Costantini C, Hahn MW, Besansky NJ. Association mapping desiccation resistance within chromosomal inversions in the African malaria vector Anopheles gambiae. Mol Ecol 2018; 28:1333-1342. [PMID: 30252170 DOI: 10.1111/mec.14880] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2018] [Revised: 09/08/2018] [Accepted: 09/10/2018] [Indexed: 12/30/2022]
Abstract
Inversion polymorphisms are responsible for many ecologically important phenotypes and are often found under balancing selection. However, the same features that ensure their large role in local adaptation-especially reduced recombination between alternate arrangements-mean that uncovering the precise loci within inversions that control these phenotypes is unachievable using standard mapping approaches. Here, we take advantage of long-term balancing selection on a pair of inversions in the mosquito Anopheles gambiae to map desiccation tolerance via pool-GWAS. Two polymorphic inversions on chromosome 2 of this species (denoted 2La and 2Rb) are associated with arid and hot conditions in Africa and are maintained in spatially and temporally heterogeneous environments. After measuring thousands of wild-caught individuals for survival under desiccation stress, we used phenotypically extreme individuals homozygous for alternative arrangements at the 2La inversion to construct pools for whole-genome sequencing. Genomewide association mapping using these pools revealed dozens of significant SNPs within both 2La and 2Rb, many of which neighboured genes controlling ion channels or related functions. Our results point to the promise of similar approaches in systems with inversions maintained by balancing selection and provide a list of candidate genes underlying the specific phenotypes controlled by the two inversions studied here.
Collapse
Affiliation(s)
- Diego Ayala
- Eck Institute for Global Health and Department of Biological Sciences, University of Notre Dame, Notre Dame, Indiana
| | - Simo Zhang
- Department of Computer Science, Indiana University, Bloomington, Indiana
| | - Mathieu Chateau
- Institut de Recherche pour le Développement, MIVEGEC (IRD, CNRS, Univ. Montpellier), Montpellier, France
| | - Caroline Fouet
- Institut de Recherche pour le Développement, MIVEGEC (IRD, CNRS, Univ. Montpellier), Montpellier, France.,Organisation de Coordination pour la lutte contre les Endémies en Afrique Centrale (OCAEC), Yaoundé, Cameroon
| | - Isabelle Morlais
- Institut de Recherche pour le Développement, MIVEGEC (IRD, CNRS, Univ. Montpellier), Montpellier, France.,Organisation de Coordination pour la lutte contre les Endémies en Afrique Centrale (OCAEC), Yaoundé, Cameroon
| | - Carlo Costantini
- Institut de Recherche pour le Développement, MIVEGEC (IRD, CNRS, Univ. Montpellier), Montpellier, France.,Organisation de Coordination pour la lutte contre les Endémies en Afrique Centrale (OCAEC), Yaoundé, Cameroon
| | - Matthew W Hahn
- Department of Computer Science, Indiana University, Bloomington, Indiana.,Department of Biology, Indiana University, Bloomington, Indiana
| | - Nora J Besansky
- Eck Institute for Global Health and Department of Biological Sciences, University of Notre Dame, Notre Dame, Indiana
| |
Collapse
|
227
|
Cheng C, Kirkpatrick M. Inversions are bigger on the X chromosome. Mol Ecol 2018; 28:1238-1245. [PMID: 30059177 DOI: 10.1111/mec.14819] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2017] [Revised: 03/22/2018] [Accepted: 04/02/2018] [Indexed: 12/22/2022]
Abstract
In many insects, X-linked inversions fix at a higher rate and are much less polymorphic than autosomal inversions. Here, we report that in Drosophila, X-linked inversions also capture 67% more genes. We estimated the number of genes captured through an approximate Bayesian computational analysis of gene orders in nine species of Drosophila. X-linked inversions fixed with a significantly larger gene content. Further, X-linked inversions of intermediate size enjoy highest fixation rate, while the fixation rate of autosomal inversions decreases with size. A less detailed analysis in Anopheles suggests a similar pattern holds in mosquitoes. We develop a population genetic model that assumes the fitness effects of inversions scale with the number of genes captured. We show that the same conditions that lead to a higher fixation rate also produce a larger size for inversions on the X.
Collapse
Affiliation(s)
- Changde Cheng
- Department of Integrative Biology, University of Texas, Austin, Texas
| | - Mark Kirkpatrick
- Department of Integrative Biology, University of Texas, Austin, Texas
| |
Collapse
|
228
|
Elworth RAL, Allen C, Benedict T, Dulworth P, Nakhleh L. ALPHA: a toolkit for Automated Local PHylogenomic Analyses. Bioinformatics 2018; 34:2848-2850. [PMID: 29562324 DOI: 10.1093/bioinformatics/bty173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2017] [Accepted: 03/16/2018] [Indexed: 11/14/2022] Open
Abstract
Summary The evolutionary histories of individual regions across a genomic alignment-called 'local genealogies'-can differ from each other, due to processes such as recombination. Elucidating and analyzing these local genealogies are important for a large number of inference tasks, including those pertaining to species phylogenies, evolutionary processes and trait mapping. In this paper, we present a toolkit for automated local phylogenomic analyses, or ALPHA. The purpose of this toolkit is to provide a wide array of functionalities for automated inference of local genealogies as well as analyses based on these local genealogies. The toolkit uses sliding windows to construct local genealogies and can compute a wide array of local phylogeny based statistics, such as the D-statistic. The toolkit comes with a graphical user interface and several import/export functionalities. Over the last few decades, much emphasis in phylogenomics has been put on developing tools for inferring species phylogenies. This toolkit complements those efforts by emphasizing the 'local' aspect of phylogenomics. Availability and implementation ALPHA is freely available for installation and use, including source code, at https://github.com/chilleo/ALPHA.
Collapse
Affiliation(s)
| | | | | | | | - Luay Nakhleh
- Computer Science, Rice University, Houston, USA.,BioSciences, Rice University, Houston, USA
| |
Collapse
|
229
|
Beckman EJ, Benham PM, Cheviron ZA, Witt C. Detecting introgression despite phylogenetic uncertainty: The case of the South American siskins. Mol Ecol 2018; 27:4350-4367. [DOI: 10.1111/mec.14795] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2017] [Revised: 05/21/2018] [Accepted: 05/23/2018] [Indexed: 12/25/2022]
Affiliation(s)
- Elizabeth J. Beckman
- Division of Biological Sciences University of Montana Missoula Montana
- Department of Biology and Museum of Southwestern Biology University of New Mexico Albuquerque New Mexico
| | - Phred M. Benham
- Division of Biological Sciences University of Montana Missoula Montana
| | | | - Christopher C. Witt
- Department of Biology and Museum of Southwestern Biology University of New Mexico Albuquerque New Mexico
| |
Collapse
|
230
|
Coughlan JM, Willis JH. Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex. Mol Ecol 2018; 28:1343-1357. [PMID: 30028906 DOI: 10.1111/mec.14804] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2018] [Revised: 07/03/2018] [Accepted: 07/06/2018] [Indexed: 01/01/2023]
Abstract
Chromosomal inversions can play an important role in adaptation, but the mechanism of their action in many natural populations remains unclear. An inversion could suppress recombination between locally beneficial alleles, thereby preventing maladaptive reshuffling with less-fit, migrant alleles. The recombination suppression hypothesis has gained much theoretical support but empirical tests are lacking. Here, we evaluated the evolutionary history and phenotypic effects of a chromosomal inversion which differentiates annual and perennial forms of Mimulus guttatus. We found that perennials likely possess the derived orientation of the inversion. In addition, this perennial orientation occurs in a second perennial species, M. decorus, where it is strongly associated with life history differences between co-occurring M. decorus and annual M. guttatus. One prediction of the recombination suppression hypothesis is that loci contributing to local adaptation will predate the inversion. To test whether the loci influencing perenniality pre-date this inversion, we mapped QTLs for life history traits that differ between annual M. guttatus and a more distantly related, collinear perennial species, M. tilingii. Consistent with the recombination suppression hypothesis, we found that this region is associated with life history in the absence of the inversion, and this association can be broken into at least two QTLs. However, the absolute phenotypic effect of the LG8 inversion region on life history is weaker in M. tilingii than in perennials which possess the inversion. Thus, while we find support for the recombination suppression hypothesis, the contribution of this inversion to life history divergence in this group is likely complex.
Collapse
Affiliation(s)
| | - John H Willis
- Biology Department, Duke University, Durham, North Carolina
| |
Collapse
|
231
|
Wiltshire RM, Bergey CM, Kayondo JK, Birungi J, Mukwaya LG, Emrich SJ, Besansky NJ, Collins FH. Reduced-representation sequencing identifies small effective population sizes of Anopheles gambiae in the north-western Lake Victoria basin, Uganda. Malar J 2018; 17:285. [PMID: 30081911 PMCID: PMC6080216 DOI: 10.1186/s12936-018-2432-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Accepted: 07/30/2018] [Indexed: 01/09/2023] Open
Abstract
BACKGROUND Malaria is the leading cause of global paediatric mortality in children below 5 years of age. The number of fatalities has reduced significantly due to an expansion of control interventions but the development of new technologies remains necessary in order to achieve elimination. Recent attention has been focused on the release of genetically modified (GM) mosquitoes into natural vector populations as a mechanism of interrupting parasite transmission but despite successful in vivo laboratory studies, a detailed population genetic assessment, which must first precede any proposed field trial, has yet to be undertaken systematically. Here, the genetic structure of Anopheles gambiae populations in north-western Lake Victoria is explored to assess their suitability as candidates for a pilot field study release of GM mosquitoes. METHODS 478 Anopheles gambiae mosquitoes were collected from six locations and a subset (N = 96) was selected for restriction site-associated DNA sequencing (RADseq). The resulting single nucleotide polymorphism (SNP) marker set was analysed for effective size (Ne), connectivity and population structure (PCA, FST). RESULTS 5175 high-quality genome-wide SNPs were identified. A principal components analysis (PCA) of the collinear genomic regions illustrated that individuals clustered in concordance with geographic origin with some overlap between sites. Genetic differentiation between populations was varied with inter-island comparisons having the highest values (median FST 0.0480-0.0846). Ne estimates were generally small (124.2-1920.3). CONCLUSIONS A reduced-representation SNP marker set for genome-wide An. gambiae genetic analysis in the north-western Lake Victoria basin is reported. Island populations demonstrated low to moderate genetic differentiation and greater structure suggesting some limitation to migration. Smaller estimates of Ne indicate that an introduced effector transgene will be more susceptible to genetic drift but to ensure that it is driven to fixation a robust gene drive mechanism will likely be needed. These findings, together with their favourable location and suitability for frequent monitoring, indicate that the Ssese Islands contain several candidate field locations, which merit further evaluation as potential GM mosquito pilot release sites.
Collapse
Affiliation(s)
- Rachel M Wiltshire
- Eck Institute for Global Health, Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, 46556, USA.
| | - Christina M Bergey
- Departments of Anthropology and Biology, The Pennsylvania State University, University Park, PA, 16802, USA
| | - Jonathan K Kayondo
- Division of Entomology and Vector Biology, Uganda Virus Research Institute, Plot No. 51-59, Nakiwogo Road, Entebbe 49, Uganda
| | - Josephine Birungi
- Division of Entomology and Vector Biology, Uganda Virus Research Institute, Plot No. 51-59, Nakiwogo Road, Entebbe 49, Uganda
| | - Louis G Mukwaya
- Division of Entomology and Vector Biology, Uganda Virus Research Institute, Plot No. 51-59, Nakiwogo Road, Entebbe 49, Uganda
| | - Scott J Emrich
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, TN, 37996, USA
| | - Nora J Besansky
- Eck Institute for Global Health, Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Frank H Collins
- Eck Institute for Global Health, Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, 46556, USA
| |
Collapse
|
232
|
Seixas FA, Boursot P, Melo-Ferreira J. The genomic impact of historical hybridization with massive mitochondrial DNA introgression. Genome Biol 2018; 19:91. [PMID: 30056805 PMCID: PMC6065068 DOI: 10.1186/s13059-018-1471-8] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Accepted: 06/25/2018] [Indexed: 12/21/2022] Open
Abstract
BACKGROUND The extent to which selection determines interspecific patterns of genetic exchange enlightens the role of adaptation in evolution and speciation. Often reported extensive interspecific introgression could be selection-driven, but also result from demographic processes, especially in cases of invasive species replacements, which can promote introgression at their invasion front. Because invasion and selective sweeps similarly mold variation, population genetics evidence for selection can only be gathered in an explicit demographic framework. The Iberian hare, Lepus granatensis, displays in its northern range extensive mitochondrial DNA introgression from L. timidus, an arctic/boreal species that it replaced locally after the last glacial maximum. We use whole-genome sequencing to infer geographic and genomic patterns of nuclear introgression and fit a neutral model of species replacement with hybridization, allowing us to evaluate how selection influenced introgression genome-wide, including for mtDNA. RESULTS Although the average nuclear and mtDNA introgression patterns contrast strongly, they fit a single demographic model of post-glacial invasive replacement of timidus by granatensis. Outliers of elevated introgression include several genes related to immunity, spermatogenesis, and mitochondrial metabolism. Introgression is reduced on the X chromosome and in low recombining regions. CONCLUSIONS General nuclear and mtDNA patterns of introgression can be explained by purely demographic processes. Hybrid incompatibilities and interplay between selection and recombination locally modulate levels of nuclear introgression. Selection promoted introgression of some genes involved in conflicts, either interspecific (parasites) or possibly cytonuclear. In the latter case, nuclear introgression could mitigate the potential negative effects of alien mtDNA on mitochondrial metabolism and male-specific traits.
Collapse
Affiliation(s)
- Fernando A Seixas
- CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Universidade do Porto, Campus Agrário de Vairão, 4485-661, Vairão, Portugal
- Departamento de Biologia, Faculdade de Ciências da Universidade do Porto, Rua Campo Alegre s/n, 4169-007, Porto, Portugal
- Institut des Sciences de l'Évolution, Université de Montpellier, CNRS, IRD, EPHE, Place Eugène Bataillon, 34095, Montpellier, France
| | - Pierre Boursot
- Institut des Sciences de l'Évolution, Université de Montpellier, CNRS, IRD, EPHE, Place Eugène Bataillon, 34095, Montpellier, France.
| | - José Melo-Ferreira
- CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Universidade do Porto, Campus Agrário de Vairão, 4485-661, Vairão, Portugal.
- Departamento de Biologia, Faculdade de Ciências da Universidade do Porto, Rua Campo Alegre s/n, 4169-007, Porto, Portugal.
| |
Collapse
|
233
|
He Z, Li X, Yang M, Wang X, Zhong C, Duke NC, Wu CI, Shi S. Speciation with gene flow via cycles of isolation and migration: insights from multiple mangrove taxa. Natl Sci Rev 2018; 6:275-288. [PMID: 31258952 PMCID: PMC6599600 DOI: 10.1093/nsr/nwy078] [Citation(s) in RCA: 69] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Allopatric speciation requiring an unbroken period of geographical isolation has been the standard model of neo-Darwinism. While doubts have been repeatedly raised, strict allopatry without any gene flow remains a plausible mechanism in most cases. To rigorously reject strict allopatry, genomic sequences superimposed on the geological records of a well-delineated geographical barrier are necessary. The Strait of Malacca, narrowly connecting the Pacific and Indian Ocean coasts, serves at different times either as a geographical barrier or a conduit of gene flow for coastal/marine species. We surveyed 1700 plants from 29 populations of 5 common mangrove species by large-scale DNA sequencing and added several whole-genome assemblies. Speciation between the two oceans is driven by cycles of isolation and gene flow due to the fluctuations in sea level leading to the opening/closing of the Strait to ocean currents. Because the time required for speciation in mangroves is longer than the isolation phases, speciation in these mangroves has proceeded through many cycles of mixing-isolation-mixing, or MIM, cycles. The MIM mechanism, by relaxing the condition of no gene flow, can promote speciation in many more geographical features than strict allopatry can. Finally, the MIM mechanism of speciation is also efficient, potentially yielding mn (m > 1) species after n cycles.
Collapse
Affiliation(s)
- Ziwen He
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, Key Laboratory of Biodiversity Dynamics and Conservation of Guangdong Higher Education Institutes, School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, China
| | - Xinnian Li
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, Key Laboratory of Biodiversity Dynamics and Conservation of Guangdong Higher Education Institutes, School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, China
| | - Ming Yang
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, Key Laboratory of Biodiversity Dynamics and Conservation of Guangdong Higher Education Institutes, School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, China
| | - Xinfeng Wang
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, Key Laboratory of Biodiversity Dynamics and Conservation of Guangdong Higher Education Institutes, School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, China
| | - Cairong Zhong
- Hainan Dongzhai Harbor National Nature Reserve Administration, Haikou 571129, China
| | - Norman C Duke
- Centre for Tropical Water and Aquatic Ecosystem Research, James Cook University, Townsville 4811, Australia
| | - Chung-I Wu
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, Key Laboratory of Biodiversity Dynamics and Conservation of Guangdong Higher Education Institutes, School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, China.,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing 100101, China.,Department of Ecology and Evolution, University of Chicago, Chicago IL 60637, USA
| | - Suhua Shi
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, Key Laboratory of Biodiversity Dynamics and Conservation of Guangdong Higher Education Institutes, School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, China
| |
Collapse
|
234
|
Dolgova O, Lao O. Evolutionary and Medical Consequences of Archaic Introgression into Modern Human Genomes. Genes (Basel) 2018; 9:E358. [PMID: 30022013 PMCID: PMC6070777 DOI: 10.3390/genes9070358] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Revised: 07/07/2018] [Accepted: 07/11/2018] [Indexed: 01/13/2023] Open
Abstract
The demographic history of anatomically modern humans (AMH) involves multiple migration events, population extinctions and genetic adaptations. As genome-wide data from complete genome sequencing becomes increasingly abundant and available even from extinct hominins, new insights of the evolutionary history of our species are discovered. It is currently known that AMH interbred with archaic hominins once they left the African continent. Current non-African human genomes carry fragments of archaic origin. This review focuses on the fitness consequences of archaic interbreeding in current human populations. We discuss new insights and challenges that researchers face when interpreting the potential impact of introgression on fitness and testing hypotheses about the role of selection within the context of health and disease.
Collapse
Affiliation(s)
- Olga Dolgova
- Population Genomics Group, Centre Nacional d'Anàlisi Genòmica, Centre de Regulació Genòmica (CRG-CNAG), Parc Científic de Barcelona, Baldiri Reixac 4, 08028 Barcelona, Catalonia, Spain.
| | - Oscar Lao
- Population Genomics Group, Centre Nacional d'Anàlisi Genòmica, Centre de Regulació Genòmica (CRG-CNAG), Parc Científic de Barcelona, Baldiri Reixac 4, 08028 Barcelona, Catalonia, Spain.
| |
Collapse
|
235
|
Fuller ZL, Leonard CJ, Young RE, Schaeffer SW, Phadnis N. Ancestral polymorphisms explain the role of chromosomal inversions in speciation. PLoS Genet 2018; 14:e1007526. [PMID: 30059505 PMCID: PMC6085072 DOI: 10.1371/journal.pgen.1007526] [Citation(s) in RCA: 48] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Revised: 08/09/2018] [Accepted: 06/29/2018] [Indexed: 01/28/2023] Open
Abstract
Understanding the role of chromosomal inversions in speciation is a fundamental problem in evolutionary genetics. Here, we perform a comprehensive reconstruction of the evolutionary histories of the chromosomal inversions in Drosophila persimilis and D. pseudoobscura. We provide a solution to the puzzling origins of the selfish Sex-Ratio arrangement in D. persimilis and uncover surprising patterns of phylogenetic discordance on this chromosome. These patterns show that, contrary to widely held views, all fixed chromosomal inversions between D. persimilis and D. pseudoobscura were already present in their ancestral population long before the species split. Our results suggest that patterns of higher genomic divergence and an association of reproductive isolation genes with chromosomal inversions may be a direct consequence of incomplete lineage sorting of ancestral polymorphisms. These findings force a reconsideration of the role of chromosomal inversions in speciation, not as protectors of existing hybrid incompatibilities, but as fertile grounds for their formation.
Collapse
Affiliation(s)
- Zachary L. Fuller
- Department of Biology, Erwin W. Mueller Laboratories, The Pennsylvania State University, University Park, PA, United States of America
| | | | - Randee E. Young
- Department of Biology, University of Utah, Salt Lake City, UT, United States of America
| | - Stephen W. Schaeffer
- Department of Biology, Erwin W. Mueller Laboratories, The Pennsylvania State University, University Park, PA, United States of America
| | - Nitin Phadnis
- Department of Biology, University of Utah, Salt Lake City, UT, United States of America
| |
Collapse
|
236
|
Wu M, Kostyun JL, Hahn MW, Moyle LC. Dissecting the basis of novel trait evolution in a radiation with widespread phylogenetic discordance. Mol Ecol 2018; 27:3301-3316. [PMID: 29953708 DOI: 10.1111/mec.14780] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2017] [Revised: 01/15/2018] [Accepted: 01/19/2018] [Indexed: 01/03/2023]
Abstract
Phylogenetic analyses of trait evolution can provide insight into the evolutionary processes that initiate and drive phenotypic diversification. However, recent phylogenomic studies have revealed extensive gene tree-species tree discordance, which can lead to incorrect inferences of trait evolution if only a single species tree is used for analysis. This phenomenon-dubbed "hemiplasy"-is particularly important to consider during analyses of character evolution in rapidly radiating groups, where discordance is widespread. Here, we generate whole-transcriptome data for a phylogenetic analysis of 14 species in the plant genus Jaltomata (the sister clade to Solanum), which has experienced rapid, recent trait evolution, including in fruit and nectar colour, and flower size and shape. Consistent with other radiations, we find evidence for rampant gene tree discordance due to incomplete lineage sorting (ILS) and to introgression events among the well-supported subclades. As both ILS and introgression increase the probability of hemiplasy, we perform several analyses that take discordance into account while identifying genes that might contribute to phenotypic evolution. Despite discordance, the history of fruit colour evolution in Jaltomata can be inferred with high confidence, and we find evidence of de novo adaptive evolution at individual genes associated with fruit colour variation. In contrast, hemiplasy appears to strongly affect inferences about floral character transitions in Jaltomata, and we identify candidate loci that could arise either from multiple lineage-specific substitutions or standing ancestral polymorphisms. Our analysis provides a generalizable example of how to manage discordance when identifying loci associated with trait evolution in a radiating lineage.
Collapse
Affiliation(s)
- Meng Wu
- Department of Biology, Indiana University, Bloomington, Indiana
| | - Jamie L Kostyun
- Department of Biology, Indiana University, Bloomington, Indiana
- Department of Plant Biology, University of Vermont, Burlington, Vermont
| | - Matthew W Hahn
- Department of Biology, Indiana University, Bloomington, Indiana
- Department of Computer Science, Indiana University, Bloomington, Indiana
| | - Leonie C Moyle
- Department of Biology, Indiana University, Bloomington, Indiana
| |
Collapse
|
237
|
Bangs MR, Douglas MR, Mussmann SM, Douglas ME. Unraveling historical introgression and resolving phylogenetic discord within Catostomus (Osteichthys: Catostomidae). BMC Evol Biol 2018; 18:86. [PMID: 29879898 PMCID: PMC5992631 DOI: 10.1186/s12862-018-1197-y] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2018] [Accepted: 05/18/2018] [Indexed: 01/25/2023] Open
Abstract
BACKGROUND Porous species boundaries can be a source of conflicting hypotheses, particularly when coupled with variable data and/or methodological approaches. Their impacts can often be magnified when non-model organisms with complex histories of reticulation are investigated. One such example is the genus Catostomus (Osteichthys, Catostomidae), a freshwater fish clade with conflicting morphological and mitochondrial phylogenies. The former is hypothesized as reflecting the presence of admixed genotypes within morphologically distinct lineages, whereas the latter is interpreted as the presence of distinct morphologies that emerged multiple times through convergent evolution. We tested these hypotheses using multiple methods, to including multispecies coalescent and concatenated approaches. Patterson's D-statistic was applied to resolve potential discord, examine introgression, and test the putative hybrid origin of two species. We also applied naïve binning to explore potential effects of concatenation. RESULTS We employed 14,007 loci generated from ddRAD sequencing of 184 individuals to derive the first highly supported nuclear phylogeny for Catostomus. Our phylogenomic analyses largely agreed with a morphological interpretation,with the exception of the placement of Xyrauchen texanus, which differs from both morphological and mitochondrial phylogenies. Additionally, our evaluation of the putative hybrid species C. columbianus revealed a lack introgression and instead matched the mitochondrial phylogeny. Furthermore, D-statistic tests clarified all discrepancies based solely on mitochondrial data, with agreement among topologies derived from concatenation and multispecies coalescent approaches. Extensive historic introgression was detected across six species-pairs. Potential endemism in the Virgin and Little Colorado Rivers was also apparent, and the former genus Pantosteus was derived as monophyletic, save for C. columbianus. CONCLUSIONS Complex reticulated histories detected herein support the hypothesis that introgression was responsible for conflicts that occurred within the mitochondrial phylogeny, and explains discrepancies found between it and previous morphological phylogenies. Additionally, the hybrid origin of C. columbianus was refuted, but with the caveat that more fine-grain sampling is still needed. Our diverse phylogenomic approaches provided largely concordant results, with naïve binning useful in exploring the single conflict. Considerable diversity was found within Catostomus across southwestern North America, with two drainages [Virgin River (UT) and Little Colorado River (AZ)] reflecting unique composition.
Collapse
Affiliation(s)
- Max R Bangs
- Department of Biological Sciences, University of Arkansas, Fayetteville, AR, 72701, USA. .,School of Fisheries, Aquaculture and Aquatic Sciences, Auburn University, Auburn, AL, 36849, USA.
| | - Marlis R Douglas
- Department of Biological Sciences, University of Arkansas, Fayetteville, AR, 72701, USA
| | - Steven M Mussmann
- Department of Biological Sciences, University of Arkansas, Fayetteville, AR, 72701, USA
| | - Michael E Douglas
- Department of Biological Sciences, University of Arkansas, Fayetteville, AR, 72701, USA
| |
Collapse
|
238
|
James S, Collins FH, Welkhoff PA, Emerson C, Godfray HCJ, Gottlieb M, Greenwood B, Lindsay SW, Mbogo CM, Okumu FO, Quemada H, Savadogo M, Singh JA, Tountas KH, Touré YT. Pathway to Deployment of Gene Drive Mosquitoes as a Potential Biocontrol Tool for Elimination of Malaria in Sub-Saharan Africa: Recommendations of a Scientific Working Group †. Am J Trop Med Hyg 2018; 98:1-49. [PMID: 29882508 PMCID: PMC5993454 DOI: 10.4269/ajtmh.18-0083] [Citation(s) in RCA: 124] [Impact Index Per Article: 20.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Accepted: 04/04/2018] [Indexed: 12/22/2022] Open
Abstract
Gene drive technology offers the promise for a high-impact, cost-effective, and durable method to control malaria transmission that would make a significant contribution to elimination. Gene drive systems, such as those based on clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein, have the potential to spread beneficial traits through interbreeding populations of malaria mosquitoes. However, the characteristics of this technology have raised concerns that necessitate careful consideration of the product development pathway. A multidisciplinary working group considered the implications of low-threshold gene drive systems on the development pathway described in the World Health Organization Guidance Framework for testing genetically modified (GM) mosquitoes, focusing on reduction of malaria transmission by Anopheles gambiae s.l. mosquitoes in Africa as a case study. The group developed recommendations for the safe and ethical testing of gene drive mosquitoes, drawing on prior experience with other vector control tools, GM organisms, and biocontrol agents. These recommendations are organized according to a testing plan that seeks to maximize safety by incrementally increasing the degree of human and environmental exposure to the investigational product. As with biocontrol agents, emphasis is placed on safety evaluation at the end of physically confined laboratory testing as a major decision point for whether to enter field testing. Progression through the testing pathway is based on fulfillment of safety and efficacy criteria, and is subject to regulatory and ethical approvals, as well as social acceptance. The working group identified several resources that were considered important to support responsible field testing of gene drive mosquitoes.
Collapse
Affiliation(s)
- Stephanie James
- Foundation for the National Institutes of Health, Bethesda, Maryland
| | | | | | | | | | - Michael Gottlieb
- Foundation for the National Institutes of Health, Bethesda, Maryland
| | - Brian Greenwood
- London School of Hygiene & Tropical Medicine, London, United Kingdom
| | | | | | - Fredros O. Okumu
- Ifakara Health Institute, Ifakara, Tanzania
- University of Glasgow, Glasgow, Scotland
- University of the Witwatersrand, Johannesburg, South Africa
| | - Hector Quemada
- Donald Danforth Plant Science Center, Saint Louis, Missouri
| | - Moussa Savadogo
- New Partnership for Africa’s Development, Ouagadougou, Burkina Faso
| | - Jerome A. Singh
- Centre for the AIDS Programme of Research in South Africa, Durban, KwaZulu-Natal, South Africa
| | - Karen H. Tountas
- Foundation for the National Institutes of Health, Bethesda, Maryland
| | - Yeya T. Touré
- University of Sciences, Techniques and Technologies of Bamako, Bamako, Mali
| |
Collapse
|
239
|
Buchman AB, Ivy T, Marshall JM, Akbari OS, Hay BA. Engineered Reciprocal Chromosome Translocations Drive High Threshold, Reversible Population Replacement in Drosophila. ACS Synth Biol 2018. [PMID: 29608276 DOI: 10.1101/088393] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Replacement of wild insect populations with transgene-bearing individuals unable to transmit disease or survive under specific environmental conditions using gene drive provides a self-perpetuating method of disease prevention. Mechanisms that require the gene drive element and linked cargo to exceed a high threshold frequency in order for spread to occur are attractive because they offer several points of control: they bring about local, but not global population replacement; and transgenes can be eliminated by reintroducing wildtypes into the population so as to drive the frequency of transgenes below the threshold frequency required for drive. Reciprocal chromosome translocations were proposed as a tool for bringing about high threshold population replacement in 1940 and 1968. However, translocations able to achieve this goal have only been reported once, in the spider mite Tetranychus urticae, a haplo-diploid species in which there is strong selection in haploid males for fit homozygotes. We report the creation of engineered translocation-bearing strains of Drosophila melanogaster, generated through targeted chromosomal breakage and homologous recombination. These strains drive high threshold population replacement in laboratory populations. While it remains to be shown that engineered translocations can bring about population replacement in wild populations, these observations suggest that further exploration of engineered translocations as a tool for controlled population replacement is warranted.
Collapse
Affiliation(s)
- Anna B Buchman
- Division of Biology and Biological Engineering , California Institute of Technology , Pasadena , California 91125 , United States
- Division of Biological Sciences , University of California , San Diego , California 92161 , United States
| | - Tobin Ivy
- Division of Biology and Biological Engineering , California Institute of Technology , Pasadena , California 91125 , United States
| | - John M Marshall
- School of Public Health , University of California , Berkeley , California 94720 , United States
| | - Omar S Akbari
- Division of Biology and Biological Engineering , California Institute of Technology , Pasadena , California 91125 , United States
- Division of Biological Sciences , University of California , San Diego , California 92161 , United States
| | - Bruce A Hay
- Division of Biology and Biological Engineering , California Institute of Technology , Pasadena , California 91125 , United States
- Division of Biological Sciences , University of California , San Diego , California 92161 , United States
| |
Collapse
|
240
|
Buchman AB, Ivy T, Marshall JM, Akbari OS, Hay BA. Engineered Reciprocal Chromosome Translocations Drive High Threshold, Reversible Population Replacement in Drosophila. ACS Synth Biol 2018; 7:1359-1370. [PMID: 29608276 DOI: 10.1021/acssynbio.7b00451] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Replacement of wild insect populations with transgene-bearing individuals unable to transmit disease or survive under specific environmental conditions using gene drive provides a self-perpetuating method of disease prevention. Mechanisms that require the gene drive element and linked cargo to exceed a high threshold frequency in order for spread to occur are attractive because they offer several points of control: they bring about local, but not global population replacement; and transgenes can be eliminated by reintroducing wildtypes into the population so as to drive the frequency of transgenes below the threshold frequency required for drive. Reciprocal chromosome translocations were proposed as a tool for bringing about high threshold population replacement in 1940 and 1968. However, translocations able to achieve this goal have only been reported once, in the spider mite Tetranychus urticae, a haplo-diploid species in which there is strong selection in haploid males for fit homozygotes. We report the creation of engineered translocation-bearing strains of Drosophila melanogaster, generated through targeted chromosomal breakage and homologous recombination. These strains drive high threshold population replacement in laboratory populations. While it remains to be shown that engineered translocations can bring about population replacement in wild populations, these observations suggest that further exploration of engineered translocations as a tool for controlled population replacement is warranted.
Collapse
Affiliation(s)
- Anna B Buchman
- Division of Biology and Biological Engineering , California Institute of Technology , Pasadena , California 91125 , United States
- Division of Biological Sciences , University of California , San Diego , California 92161 , United States
| | - Tobin Ivy
- Division of Biology and Biological Engineering , California Institute of Technology , Pasadena , California 91125 , United States
| | - John M Marshall
- School of Public Health , University of California , Berkeley , California 94720 , United States
| | - Omar S Akbari
- Division of Biology and Biological Engineering , California Institute of Technology , Pasadena , California 91125 , United States
- Division of Biological Sciences , University of California , San Diego , California 92161 , United States
| | - Bruce A Hay
- Division of Biology and Biological Engineering , California Institute of Technology , Pasadena , California 91125 , United States
- Division of Biological Sciences , University of California , San Diego , California 92161 , United States
| |
Collapse
|
241
|
Anselmetti Y, Duchemin W, Tannier E, Chauve C, Bérard S. Phylogenetic signal from rearrangements in 18 Anopheles species by joint scaffolding extant and ancestral genomes. BMC Genomics 2018; 19:96. [PMID: 29764366 PMCID: PMC5954271 DOI: 10.1186/s12864-018-4466-7] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022] Open
Abstract
Background Genomes rearrangements carry valuable information for phylogenetic inference or the elucidation of molecular mechanisms of adaptation. However, the detection of genome rearrangements is often hampered by current deficiencies in data and methods: Genomes obtained from short sequence reads have generally very fragmented assemblies, and comparing multiple gene orders generally leads to computationally intractable algorithmic questions. Results We present a computational method, ADseq, which, by combining ancestral gene order reconstruction, comparative scaffolding and de novo scaffolding methods, overcomes these two caveats. ADseq provides simultaneously improved assemblies and ancestral genomes, with statistical supports on all local features. Compared to previous comparative methods, it runs in polynomial time, it samples solutions in a probabilistic space, and it can handle a significantly larger gene complement from the considered extant genomes, with complex histories including gene duplications and losses. We use ADseq to provide improved assemblies and a genome history made of duplications, losses, gene translocations, rearrangements, of 18 complete Anopheles genomes, including several important malaria vectors. We also provide additional support for a differentiated mode of evolution of the sex chromosome and of the autosomes in these mosquito genomes. Conclusions We demonstrate the method’s ability to improve extant assemblies accurately through a procedure simulating realistic assembly fragmentation. We study a debated issue regarding the phylogeny of the Gambiae complex group of Anopheles genomes in the light of the evolution of chromosomal rearrangements, suggesting that the phylogenetic signal they carry can differ from the phylogenetic signal carried by gene sequences, more prone to introgression. Electronic supplementary material The online version of this article (10.1186/s12864-018-4466-7) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Yoann Anselmetti
- ISEM, Université de Montpellier, CNRS, IRD, EPHE, Montpellier, France.,Univ Lyon, Université Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Evolutive UMR5558, 43 Boulevard du 11 novembre 1918, Villeurbanne cedex, 69622, France
| | - Wandrille Duchemin
- Univ Lyon, Université Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Evolutive UMR5558, 43 Boulevard du 11 novembre 1918, Villeurbanne cedex, 69622, France.,INRIA Grenoble - Rhône-Alpes, 655 Avenue de l'Europe, Montbonnot-Saint-Martin, 38330, France
| | - Eric Tannier
- Univ Lyon, Université Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Evolutive UMR5558, 43 Boulevard du 11 novembre 1918, Villeurbanne cedex, 69622, France.,INRIA Grenoble - Rhône-Alpes, 655 Avenue de l'Europe, Montbonnot-Saint-Martin, 38330, France
| | - Cedric Chauve
- Department of Mathematics, Simon Fraser University, 8888 University Drive, Burnaby, V5A1S6, BC, Canada
| | - Sèverine Bérard
- ISEM, Université de Montpellier, CNRS, IRD, EPHE, Montpellier, France.
| |
Collapse
|
242
|
Analysis of natural female post-mating responses of Anopheles gambiae and Anopheles coluzzii unravels similarities and differences in their reproductive ecology. Sci Rep 2018; 8:6594. [PMID: 29700344 PMCID: PMC5920108 DOI: 10.1038/s41598-018-24923-w] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2017] [Accepted: 03/20/2018] [Indexed: 11/18/2022] Open
Abstract
Anopheles gambiae and An. coluzzii, the two most important malaria vectors in sub-Saharan Africa, are recently radiated sibling species that are reproductively isolated even in areas of sympatry. In females from these species, sexual transfer of male accessory gland products, including the steroid hormone 20-hydroxyecdysone (20E), induces vast behavioral, physiological, and transcriptional changes that profoundly shape their post-mating ecology, and that may have contributed to the insurgence of post-mating, prezygotic reproductive barriers. As these barriers can be detected by studying transcriptional changes induced by mating, we set out to analyze the post-mating response of An. gambiae and An. coluzzii females captured in natural mating swarms in Burkina Faso. While the molecular pathways shaping short- and long-term mating-induced changes are largely conserved in females from the two species, we unravel significant inter-specific differences that suggest divergent regulation of key reproductive processes such as egg development, processing of seminal secretion, and mating behavior, that may have played a role in reproductive isolation. Interestingly, a number of these changes occur in genes previously shown to be regulated by the sexual transfer of 20E and may be due to divergent utilization of this steroid hormone in the two species.
Collapse
|
243
|
Artemov GN, Bondarenko SM, Naumenko AN, Stegniy VN, Sharakhova MV, Sharakhov IV. Partial-arm translocations in evolution of malaria mosquitoes revealed by high-coverage physical mapping of the Anopheles atroparvus genome. BMC Genomics 2018; 19:278. [PMID: 29688842 PMCID: PMC5914054 DOI: 10.1186/s12864-018-4663-4] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2017] [Accepted: 04/12/2018] [Indexed: 02/06/2023] Open
Abstract
Background Malaria mosquitoes have had a remarkable stability in the number of chromosomes in their karyotype (2n = 6) during 100 million years of evolution. Moreover, autosomal arms were assumed to maintain their integrity even if their associations with each other changed via whole-arm translocations. Here we use high-coverage comparative physical genome mapping of three Anopheles species to test the extent of evolutionary conservation of chromosomal arms in malaria mosquitoes. Results In this study, we developed a physical genome map for Anopheles atroparvus, one of the dominant malaria vectors in Europe. Using fluorescence in situ hybridization (FISH) of DNA probes with the ovarian nurse cell polytene chromosomes and synteny comparison, we anchored 56 genomic scaffolds to the An. atroparvus chromosomes. The obtained physical map represents 89.6% of the An. atroparvus genome. This genome has the second highest mapping coverage among Anophelinae assemblies after An. albimanus, which has 98.2% of the genome assigned to its chromosomes. A comparison of the An. atroparvus, An. albimanus, and An. gambiae genomes identified partial-arm translocations between the autosomal arms that break down the integrity of chromosome elements in evolution affecting the structure of the genetic material in the pericentromeric regions. Unlike An. atroparvus and An. albimanus, all chromosome elements of An. gambiae are fully syntenic with chromosome elements of the putative ancestral Anopheles karyotype. We also detected nonrandom distribution of large conserved synteny blocks and confirmed a higher rate of inversion fixation in the X chromosome compared with autosomes. Conclusions Our study demonstrates the power of physical mapping for understanding the genome evolution in malaria mosquitoes. The results indicate that syntenic relationships among chromosome elements of Anopheles species have not been fully preserved because of multiple partial-arm translocations. Electronic supplementary material The online version of this article (10.1186/s12864-018-4663-4) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Gleb N Artemov
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk, 634050, Russia
| | - Semen M Bondarenko
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk, 634050, Russia
| | - Anastasia N Naumenko
- Department of Entomology, Fralin Life Science Institute, Virginia Polytechnic Institute and State University, 360 West Campus Drive, Blacksburg, VA, 24061, USA
| | - Vladimir N Stegniy
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk, 634050, Russia
| | - Maria V Sharakhova
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk, 634050, Russia. .,Department of Entomology, Fralin Life Science Institute, Virginia Polytechnic Institute and State University, 360 West Campus Drive, Blacksburg, VA, 24061, USA.
| | - Igor V Sharakhov
- Laboratory of Ecology, Genetics and Environmental Protection, Tomsk State University, 36 Lenin Avenue, Tomsk, 634050, Russia. .,Department of Entomology, Fralin Life Science Institute, Virginia Polytechnic Institute and State University, 360 West Campus Drive, Blacksburg, VA, 24061, USA.
| |
Collapse
|
244
|
Schrider DR, Ayroles J, Matute DR, Kern AD. Supervised machine learning reveals introgressed loci in the genomes of Drosophila simulans and D. sechellia. PLoS Genet 2018; 14:e1007341. [PMID: 29684059 PMCID: PMC5933812 DOI: 10.1371/journal.pgen.1007341] [Citation(s) in RCA: 69] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2017] [Revised: 05/03/2018] [Accepted: 03/28/2018] [Indexed: 12/30/2022] Open
Abstract
Hybridization and gene flow between species appears to be common. Even though it is clear that hybridization is widespread across all surveyed taxonomic groups, the magnitude and consequences of introgression are still largely unknown. Thus it is crucial to develop the statistical machinery required to uncover which genomic regions have recently acquired haplotypes via introgression from a sister population. We developed a novel machine learning framework, called FILET (Finding Introgressed Loci via Extra-Trees) capable of revealing genomic introgression with far greater power than competing methods. FILET works by combining information from a number of population genetic summary statistics, including several new statistics that we introduce, that capture patterns of variation across two populations. We show that FILET is able to identify loci that have experienced gene flow between related species with high accuracy, and in most situations can correctly infer which population was the donor and which was the recipient. Here we describe a data set of outbred diploid Drosophila sechellia genomes, and combine them with data from D. simulans to examine recent introgression between these species using FILET. Although we find that these populations may have split more recently than previously appreciated, FILET confirms that there has indeed been appreciable recent introgression (some of which might have been adaptive) between these species, and reveals that this gene flow is primarily in the direction of D. simulans to D. sechellia. Understanding the extent to which species or diverged populations hybridize in nature is crucially important if we are to understand the speciation process. Accordingly numerous research groups have developed methodology for finding the genetic evidence of such introgression. In this report we develop a supervised machine learning approach for uncovering loci which have introgressed across species boundaries. We show that our method, FILET, has greater accuracy and power than competing methods in discovering introgression, and in addition can detect the directionality associated with the gene flow between species. Using whole genome sequences from Drosophila simulans and Drosophila sechellia we show that FILET discovers quite extensive introgression between these species that has occurred mostly from D. simulans to D. sechellia. Our work highlights the complex process of speciation even within a well-studied system and points to the growing importance of supervised machine learning in population genetics.
Collapse
Affiliation(s)
- Daniel R. Schrider
- Department of Genetics, Rutgers University, Piscataway, New Jersey, United States of America
- Human Genetics Institute of New Jersey, Rutgers University, Piscataway, New Jersey, United States of America
- * E-mail:
| | - Julien Ayroles
- Ecology and Evolutionary Biology Department, Princeton University, Princeton, New Jersey, United States of America
- Lewis Sigler Institute for Integrative Genomics, Princeton University, Princeton, New Jersey, United States of America
| | - Daniel R. Matute
- Biology Department, University of North Carolina, Chapel Hill, North Carolina, United States of America
| | - Andrew D. Kern
- Department of Genetics, Rutgers University, Piscataway, New Jersey, United States of America
- Human Genetics Institute of New Jersey, Rutgers University, Piscataway, New Jersey, United States of America
| |
Collapse
|
245
|
Pélissié B, Crossley MS, Cohen ZP, Schoville SD. Rapid evolution in insect pests: the importance of space and time in population genomics studies. CURRENT OPINION IN INSECT SCIENCE 2018; 26:8-16. [PMID: 29764665 DOI: 10.1016/j.cois.2017.12.008] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/17/2017] [Revised: 12/28/2017] [Accepted: 12/31/2017] [Indexed: 06/08/2023]
Abstract
Pest species in agroecosystems often exhibit patterns of rapid evolution to environmental and human-imposed selection pressures. Although the role of adaptive processes is well accepted, few insect pests have been studied in detail and most research has focused on selection at insecticide resistance candidate genes. Emerging genomic datasets provide opportunities to detect and quantify selection in insect pest populations, and address long-standing questions about mechanisms underlying rapid evolutionary change. We examine the strengths of recent studies that stratify population samples both in space (along environmental gradients and comparing ancestral vs. derived populations) and in time (using chronological sampling, museum specimens and comparative phylogenomics), resulting in critical insights on evolutionary processes, and providing new directions for studying pests in agroecosystems.
Collapse
Affiliation(s)
- Benjamin Pélissié
- University of Wisconsin-Madison, Department of Entomology, 1630 Linden Drive, 637-643 Russell Labs, Madison, WI 53706, USA.
| | - Michael S Crossley
- University of Wisconsin-Madison, Department of Entomology, 1630 Linden Drive, 637-643 Russell Labs, Madison, WI 53706, USA
| | - Zachary Paul Cohen
- University of Wisconsin-Madison, Department of Entomology, 1630 Linden Drive, 637-643 Russell Labs, Madison, WI 53706, USA
| | - Sean D Schoville
- University of Wisconsin-Madison, Department of Entomology, 1630 Linden Drive, 637-643 Russell Labs, Madison, WI 53706, USA
| |
Collapse
|
246
|
Hanemaaijer MJ, Houston PD, Collier TC, Norris LC, Fofana A, Lanzaro GC, Cornel AJ, Lee Y. Mitochondrial genomes of Anopheles arabiensis, An. gambiae and An. coluzzii show no clear species division. F1000Res 2018; 7:347. [PMID: 31069048 PMCID: PMC6489993 DOI: 10.12688/f1000research.13807.2] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 03/08/2019] [Indexed: 12/05/2022] Open
Abstract
Here we report the complete mitochondrial sequences of 70 individual field collected mosquito specimens from throughout Sub-Saharan Africa. We generated this dataset to identify species specific markers for the following Anopheles species and chromosomal forms: An. arabiensis, An. coluzzii (The Forest and Mopti chromosomal forms) and An. gambiae (The Bamako and Savannah chromosomal forms). The raw Illumina sequencing reads were mapped to the NC_002084 reference mitogenome sequence. A total of 783 single nucleotide polymorphisms (SNPs) were detected on the mitochondrial genome, of which 460 are singletons (58.7%). None of these SNPs are suitable as molecular markers to distinguish among An. arabiensis, An. coluzzii and An. gambiae or any of the chromosomal forms. The lack of species or chromosomal form specific markers is also reflected in the constructed phylogenetic tree, which shows no clear division among the operational taxonomic units considered here.
Collapse
Affiliation(s)
- Mark J. Hanemaaijer
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Parker D. Houston
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Travis C. Collier
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Laura C. Norris
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Abdrahamane Fofana
- Malaria Research and Training Center, University of Bamako, Bamako, E2528, Mali
| | - Gregory C. Lanzaro
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Anthony J. Cornel
- Mosquito Control Research Laboratory, Kearney Agricultural Center, Department of Entomology and Nematology, University of California Davis, Davis, CA, 93648, USA
| | - Yoosook Lee
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| |
Collapse
|
247
|
Hanemaaijer MJ, Houston PD, Collier TC, Norris LC, Fofana A, Lanzaro GC, Cornel AJ, Lee Y. Mitochondrial genomes of Anopheles arabiensis, An. gambiae and An. coluzzii show no clear species division. F1000Res 2018; 7:347. [PMID: 31069048 PMCID: PMC6489993 DOI: 10.12688/f1000research.13807.1] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 02/14/2018] [Indexed: 10/12/2023] Open
Abstract
Here we report the complete mitochondrial sequences of 70 individual field collected mosquito specimens from throughout Sub-Saharan Africa. We generated this dataset to identify species specific markers for the following Anopheles species and chromosomal forms: An. arabiensis, An. coluzzii (The Forest and Mopti chromosomal forms) and An. gambiae (The Bamako and Savannah chromosomal forms). The raw Illumina sequencing reads were mapped to the NC_002084 reference mitogenome sequence. A total of 783 single nucleotide polymorphisms (SNPs) were detected on the mitochondrial genome, of which 460 are singletons (58.7%). None of these SNPs are suitable as molecular markers to distinguish among An. arabiensis, An. coluzzii and An. gambiae or any of the chromosomal forms. The lack of species or chromosomal form specific markers is also reflected in the constructed phylogenetic tree, which shows no clear division among the operational taxonomic units considered here.
Collapse
Affiliation(s)
- Mark J. Hanemaaijer
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Parker D. Houston
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Travis C. Collier
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Laura C. Norris
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Abdrahamane Fofana
- Malaria Research and Training Center, University of Bamako, Bamako, E2528, Mali
| | - Gregory C. Lanzaro
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| | - Anthony J. Cornel
- Mosquito Control Research Laboratory, Kearney Agricultural Center, Department of Entomology and Nematology, University of California Davis, Davis, CA, 93648, USA
| | - Yoosook Lee
- Vector Genetics Laboratory, Department of Pathology, Microbiology and Immunology, University of California Davis , Davis, CA, 95616, USA
| |
Collapse
|
248
|
Sougoufara S, Thiaw O, Cailleau A, Diagne N, Harry M, Bouganali C, Sembène PM, Doucoure S, Sokhna C. The Impact of Periodic Distribution Campaigns of Long-Lasting Insecticidal-Treated Bed Nets on Malaria Vector Dynamics and Human Exposure in Dielmo, Senegal. Am J Trop Med Hyg 2018; 98:1343-1352. [PMID: 29557325 DOI: 10.4269/ajtmh.17-0009] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
The implementation of long-lasting insecticidal-treated bed nets (LLINs) has contributed to halving the mortality rate due to malaria since 2000 in sub-Saharan Africa. These tools are highly effective against indoor-feeding malaria vectors. Thus, to achieve the World Health Assembly's new target to reduce the burden of malaria over the next 15 years by 90%, it is necessary to understand how the spatiotemporal dynamics of malaria vectors and human exposure to bites is modified in the context of scaling up global efforts to control malaria transmission. This study was conducted in Dielmo, a Senegalese village, after the introduction of LLINs and two rounds of LLINs renewals. Data analysis showed that implementation of LLINs correlated with a significant decrease in the biting densities of the main malaria vectors, Anopheles gambiae s.l. and Anopheles funestus, reducing malaria transmission. Other environment factors likely contributed to the decrease in An. funestus, but this trend was enhanced with the introduction of LLINs. The bulk of bites occurred during sleeping hours, but the residual vector populations of An. gambiae s.l. and An. funestus had an increased propensity to bite outdoors, so a risk of infectious bites remained for LLINs users. These results highlight the need to increase the level and correct use of LLINs and to combine this intervention with complementary control measures against residual exposure, such as spatial repellents and larval source management, to achieve the goal of eliminating malaria transmission.
Collapse
Affiliation(s)
- Seynabou Sougoufara
- Aix Marseille University, Institut de Recherche pour le Développement (IDR) (Dakar, Marseille, Papeete), AP-HM, Institut Hospitalo-Universitaire-Méditerranée Infection, UMR Vecteurs-Infections Tropicales et Méditerranéennes (VITROME), Marseille, France.,Département de Biologie Animale, Faculté des Sciences et Techniques, Université Cheikh Anta Diop de Dakar, Dakar Fann, Sénégal
| | - Omar Thiaw
- Aix Marseille University, Institut de Recherche pour le Développement (IDR) (Dakar, Marseille, Papeete), AP-HM, Institut Hospitalo-Universitaire-Méditerranée Infection, UMR Vecteurs-Infections Tropicales et Méditerranéennes (VITROME), Marseille, France
| | - Aurélie Cailleau
- Centre Suisse de Recherches Scientifiques en Cote d'Ivoire (CSRS), Yopougon, Abidjan, Côte d'Ivoire.,Unité d'Entomologie Médicale (UME), Institut Pasteur Dakar, Dakar, Sénégal
| | - Nafissatou Diagne
- Aix Marseille University, Institut de Recherche pour le Développement (IDR) (Dakar, Marseille, Papeete), AP-HM, Institut Hospitalo-Universitaire-Méditerranée Infection, UMR Vecteurs-Infections Tropicales et Méditerranéennes (VITROME), Marseille, France
| | - Myriam Harry
- UMR Évolution, Génomes, Comportement, Écologie (EGCE) Centre National de la Recherche Scientifique (CNRS), IRD-University Paris-Sud, IDEEV, University Paris-Saclay, Gif-sur-Yvette Cedex, France
| | - Charles Bouganali
- Aix Marseille University, Institut de Recherche pour le Développement (IDR) (Dakar, Marseille, Papeete), AP-HM, Institut Hospitalo-Universitaire-Méditerranée Infection, UMR Vecteurs-Infections Tropicales et Méditerranéennes (VITROME), Marseille, France
| | - Pape M Sembène
- Département de Biologie Animale, Faculté des Sciences et Techniques, Université Cheikh Anta Diop de Dakar, Dakar Fann, Sénégal
| | - Souleymane Doucoure
- Aix Marseille University, Institut de Recherche pour le Développement (IDR) (Dakar, Marseille, Papeete), AP-HM, Institut Hospitalo-Universitaire-Méditerranée Infection, UMR Vecteurs-Infections Tropicales et Méditerranéennes (VITROME), Marseille, France
| | - Cheikh Sokhna
- Aix Marseille University, Institut de Recherche pour le Développement (IDR) (Dakar, Marseille, Papeete), AP-HM, Institut Hospitalo-Universitaire-Méditerranée Infection, UMR Vecteurs-Infections Tropicales et Méditerranéennes (VITROME), Marseille, France
| |
Collapse
|
249
|
Kearns AM, Restani M, Szabo I, Schrøder-Nielsen A, Kim JA, Richardson HM, Marzluff JM, Fleischer RC, Johnsen A, Omland KE. Genomic evidence of speciation reversal in ravens. Nat Commun 2018; 9:906. [PMID: 29500409 PMCID: PMC5834606 DOI: 10.1038/s41467-018-03294-w] [Citation(s) in RCA: 85] [Impact Index Per Article: 14.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2016] [Accepted: 02/01/2018] [Indexed: 01/10/2023] Open
Abstract
Many species, including humans, have emerged via complex reticulate processes involving hybridisation. Under certain circumstances, hybridisation can cause distinct lineages to collapse into a single lineage with an admixed mosaic genome. Most known cases of such 'speciation reversal' or 'lineage fusion' involve recently diverged lineages and anthropogenic perturbation. Here, we show that in western North America, Common Ravens (Corvus corax) have admixed mosaic genomes formed by the fusion of non-sister lineages ('California' and 'Holarctic') that diverged ~1.5 million years ago. Phylogenomic analyses and concordant patterns of geographic structuring in mtDNA, genome-wide SNPs and nuclear introns demonstrate long-term admixture and random interbreeding between the non-sister lineages. In contrast, our genomic data support reproductive isolation between Common Ravens and Chihuahuan Ravens (C. cryptoleucus) despite extensive geographic overlap and a sister relationship between Chihuahuan Ravens and the California lineage. These data suggest that the Common Raven genome was formed by secondary lineage fusion and most likely represents a case of ancient speciation reversal that occurred without anthropogenic causes.
Collapse
Affiliation(s)
- Anna M Kearns
- Natural History Museum, University of Oslo, P.O. Box 1172 Blindern, 0318, Oslo, Norway.
- Department of Biological Sciences, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD, 21250, USA.
- Center for Conservation Genomics, Smithsonian Conservation Biology Institute, National Zoological Park, Washington, 20013-7012, DC, USA.
| | - Marco Restani
- Department of Biological Sciences, St. Cloud State University, 720 Fourth Avenue, St. Cloud, MN, 56301-4498, USA
| | - Ildiko Szabo
- Cowan Tetrapod Collection, Beaty Biodiversity Museum, University of British Columbia, 2212 Main Mall, Vancouver, BC, V6T 1Z4, Canada
| | | | - Jin Ah Kim
- Department of Biological Sciences, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD, 21250, USA
| | - Hayley M Richardson
- Department of Biological Sciences, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD, 21250, USA
| | - John M Marzluff
- School of Environmental and Forest Sciences, University of Washington, Box 352100, Seattle, WA, 98195, USA
| | - Robert C Fleischer
- Center for Conservation Genomics, Smithsonian Conservation Biology Institute, National Zoological Park, Washington, 20013-7012, DC, USA
| | - Arild Johnsen
- Natural History Museum, University of Oslo, P.O. Box 1172 Blindern, 0318, Oslo, Norway
| | - Kevin E Omland
- Department of Biological Sciences, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD, 21250, USA
| |
Collapse
|
250
|
Knowles LL, Huang H, Sukumaran J, Smith SA. A matter of phylogenetic scale: Distinguishing incomplete lineage sorting from lateral gene transfer as the cause of gene tree discord in recent versus deep diversification histories. AMERICAN JOURNAL OF BOTANY 2018; 105:376-384. [PMID: 29710372 DOI: 10.1002/ajb2.1064] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2017] [Accepted: 01/30/2018] [Indexed: 06/08/2023]
Abstract
PREMISE OF THE STUDY Discordant gene trees are commonly encountered when sequences from thousands of loci are applied to estimate phylogenetic relationships. Several processes contribute to this discord. Yet, we have no methods that jointly model different sources of conflict when estimating phylogenies. An alternative to analyzing entire genomes or all the sequenced loci is to identify a subset of loci for phylogenetic analysis. If we can identify data partitions that are most likely to reflect descent from a common ancestor (i.e., discordant loci that indeed reflect incomplete lineage sorting [ILS], as opposed to some other process, such as lateral gene transfer [LGT]), we can analyze this subset using powerful coalescent-based species-tree approaches. METHODS Test data sets were simulated where discord among loci could arise from ILS and LGT. Data sets where analyzed using the newly developed program CLASSIPHY (Huang et al., ) to assess whether our ability to distinguish the cause of discord among loci varied when ILS and LGT occurred in the recent versus deep past and whether the accuracy of these inferences were affected by the mutational process. KEY RESULTS We show that accuracy of probabilistic classification of individual loci by the cause of discord differed when ILS and LGT events occurred more recently compared with the distant past and that the signal-to-noise ratio arising from the mutational process contributes to difficulties in inferring LGT data partitions. CONCLUSIONS We discuss our findings in terms of the promise and limitations of identifying subsets of loci for species-tree inference that will not violate the underlying coalescent model (i.e., data partitions in which ILS, and not LGT, contributes to discord). We also discuss the empirical implications of our work given the many recalcitrant nodes in the tree of life (e.g., origins of angiosperms, amniotes, or Neoaves), and recent arguments for concatenating loci.
Collapse
Affiliation(s)
- L Lacey Knowles
- Department of Ecology and Evolutionary Biology, Museum of Zoology, University of Michigan, 1109 Geddes Avenue, Ann Arbor, MI, 48109-1079, USA
| | - Huateng Huang
- Department of Ecology and Evolutionary Biology, Museum of Zoology, University of Michigan, 1109 Geddes Avenue, Ann Arbor, MI, 48109-1079, USA
| | - Jeet Sukumaran
- Department of Ecology and Evolutionary Biology, Museum of Zoology, University of Michigan, 1109 Geddes Avenue, Ann Arbor, MI, 48109-1079, USA
| | - Stephen A Smith
- Department of Ecology and Evolutionary Biology, Museum of Zoology, University of Michigan, 1109 Geddes Avenue, Ann Arbor, MI, 48109-1079, USA
| |
Collapse
|