1
|
Ertl HA, Hill MS, Wittkopp PJ. Differential Grainy head binding correlates with variation in chromatin structure and gene expression in Drosophila melanogaster. BMC Genomics 2022; 23:854. [PMID: 36575386 PMCID: PMC9795675 DOI: 10.1186/s12864-022-09082-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Accepted: 12/14/2022] [Indexed: 12/28/2022] Open
Abstract
Phenotypic evolution is often caused by variation in gene expression resulting from altered gene regulatory mechanisms. Genetic variation affecting chromatin remodeling has been identified as a potential source of variable gene expression; however, the roles of specific chromatin remodeling factors remain unclear. Here, we address this knowledge gap by examining the relationship between variation in gene expression, variation in chromatin structure, and variation in binding of the pioneer factor Grainy head between imaginal wing discs of two divergent strains of Drosophila melanogaster and their F1 hybrid. We find that (1) variation in Grainy head binding is mostly due to sequence changes that act in cis but are located outside of the canonical Grainy head binding motif, (2) variation in Grainy head binding correlates with changes in chromatin accessibility, and (3) this variation in chromatin accessibility, coupled with variation in Grainy head binding, correlates with variation in gene expression in some cases but not others. Interactions among these three molecular layers is complex, but these results suggest that genetic variation affecting the binding of pioneer factors contributes to variation in chromatin remodeling and the evolution of gene expression.
Collapse
Affiliation(s)
- Henry A. Ertl
- grid.214458.e0000000086837370Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI 48109 USA
| | - Mark S. Hill
- grid.214458.e0000000086837370Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI 48109 USA ,grid.83440.3b0000000121901201Present address: Cancer Evolution and Genome Instability Laboratory, University College London Cancer Institute and The Francis Crick Institute, London, UK
| | - Patricia J. Wittkopp
- grid.214458.e0000000086837370Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI 48109 USA ,grid.214458.e0000000086837370Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA
| |
Collapse
|
2
|
Li NK, Corander J, Grad YH, Chang HH. Discovering recent selection forces shaping the evolution of dengue viruses based on polymorphism data across geographic scales. Virus Evol 2022; 8:veac108. [PMID: 36601300 PMCID: PMC9789396 DOI: 10.1093/ve/veac108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Revised: 09/23/2022] [Accepted: 11/28/2022] [Indexed: 11/30/2022] Open
Abstract
Incomplete selection makes it challenging to infer selection on genes at short time scales, especially for microorganisms, due to stronger linkage between loci. However, in many cases, the selective force changes with environment, time, or other factors, and it is of great interest to understand selective forces at this level to answer relevant biological questions. We developed a new method that uses the change in dN /dS , instead of the absolute value of dN /dS , to infer the dominating selective force based on sequence data across geographical scales. If a gene was under positive selection, dN /dS was expected to increase through time, whereas if a gene was under negative selection, dN /dS was expected to decrease through time. Assuming that the migration rate decreased and the divergence time between samples increased from between-continent, within-continent different-country, to within-country level, dN /dS of a gene dominated by positive selection was expected to increase with increasing geographical scales, and the opposite trend was expected in the case of negative selection. Motivated by the McDonald-Kreitman (MK) test, we developed a pairwise MK test to assess the statistical significance of detected trends in dN /dS . Application of the method to a global sample of dengue virus genomes identified multiple significant signatures of selection in both the structural and non-structural proteins. Because this method does not require allele frequency estimates and uses synonymous mutations for comparison, it is less prone to sampling error, providing a way to infer selection forces within species using publicly available genomic data from locations over broad geographical scales.
Collapse
Affiliation(s)
- Nien-Kung Li
- Department of Life Science & Institute of Bioinformatics and Structural Biology, National Tsing Hua University, 101, Section 2, Kuang-Fu Road, Hsinchu 300044, Taiwan
| | - Jukka Corander
- Helsinki Institute for Information Technology, Department of Mathematics and Statistics, University of Helsinki, Yliopistonkatu 3, Helsinki 00014, Finland,Department of Biostatistics, University of Oslo, Domus Medica Gaustad Sognsvannsveien 9, Oslo 0372, Norway,Parasites and Microbes, The Wellcome Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, UK
| | | | | |
Collapse
|
3
|
Rodrigues MF, Vibranovski MD, Cogni R. Clinal and seasonal changes are correlated in Drosophila melanogaster natural populations. Evolution 2021; 75:2042-2054. [PMID: 34184262 DOI: 10.1111/evo.14300] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2020] [Revised: 06/08/2021] [Accepted: 06/14/2021] [Indexed: 12/22/2022]
Abstract
Spatial and seasonal variations in the environment are ubiquitous. Environmental heterogeneity can affect natural populations and lead to covariation between environment and allele frequencies. Drosophila melanogaster is known to harbor polymorphisms that change both with latitude and seasons. Identifying the role of selection in driving these changes is not trivial, because nonadaptive processes can cause similar patterns. Given the environment changes in similar ways across seasons and along the latitudinal gradient, one promising approach may be to look for parallelism between clinal and seasonal changes. Here, we test whether there is a genome-wide correlation between clinal and seasonal changes, and whether the pattern is consistent with selection. Allele frequency estimates were obtained from pooled samples from seven different locations along the east coast of the United States, and across seasons within Pennsylvania. We show that there is a genome-wide correlation between clinal and seasonal variations, which cannot be explained by linked selection alone. This pattern is stronger in genomic regions with higher functional content, consistent with natural selection. We derive a way to biologically interpret these correlations and show that around 3.7% of the common, autosomal variants could be under parallel seasonal and spatial selection. Our results highlight the contribution of natural selection in driving fluctuations in allele frequencies in natural fly populations and point to a shared genomic basis to climate adaptation that happens over space and time in D. melanogaster.
Collapse
Affiliation(s)
- Murillo F Rodrigues
- Department of Genetics and Evolutionary Biology, Institute of Biosciences, University of Sao Paulo, Sao Paulo, 05508-090, Brazil.,Current Address: Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, 97403
| | - Maria D Vibranovski
- Department of Genetics and Evolutionary Biology, Institute of Biosciences, University of Sao Paulo, Sao Paulo, 05508-090, Brazil
| | - Rodrigo Cogni
- Department of Ecology, Institute of Biosciences, University of Sao Paulo, Sao Paulo, 05508-090, Brazil
| |
Collapse
|
4
|
Signor S. Transposable elements in individual genotypes of Drosophila simulans. Ecol Evol 2020; 10:3402-3412. [PMID: 32273997 PMCID: PMC7141027 DOI: 10.1002/ece3.6134] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2019] [Revised: 01/25/2020] [Accepted: 02/04/2020] [Indexed: 11/17/2022] Open
Abstract
Transposable elements are abundant, dynamic components of the genome that affect organismal phenotypes and fitness. In Drosophila melanogaster, they have increased in abundance as the species spread out of Africa, and different populations differ in their transposable element content. However, very little is currently known about how transposable elements differ between individual genotypes, and how that relates to the population dynamics of transposable elements overall. The sister species of D. melanogaster, D. simulans, has also recently become cosmopolitan, and panels of inbred genotypes exist from cosmopolitan and African flies. Therefore, we can determine whether the differences in colonizing populations are repeated in D. simulans, what the dynamics of transposable elements are in individual genotypes, and how that compares to wild flies. After estimating copy number in cosmopolitan and African D. simulans, I find that transposable element load is higher in flies from cosmopolitan populations. In addition, transposable element load varies considerably between populations, between genotypes, but not overall between wild and inbred lines. Certain genotypes either contain active transposable elements or are more permissive of transposition and accumulate copies of particular transposable elements. Overall, it is important to quantify genotype-specific transposable element dynamics as well as population averages to understand the dynamics of transposable element accumulation over time.
Collapse
Affiliation(s)
- Sarah Signor
- Department of Biological SciencesNorth Dakota State UniversityFargoNDUSA
| |
Collapse
|
5
|
Mateo L, Rech GE, González J. Genome-wide patterns of local adaptation in Western European Drosophila melanogaster natural populations. Sci Rep 2018; 8:16143. [PMID: 30385770 PMCID: PMC6212444 DOI: 10.1038/s41598-018-34267-0] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Accepted: 10/12/2018] [Indexed: 12/21/2022] Open
Abstract
Signatures of spatially varying selection have been investigated both at the genomic and transcriptomic level in several organisms. In Drosophila melanogaster, the majority of these studies have analyzed North American and Australian populations, leading to the identification of several loci and traits under selection. However, several studies based mainly in North American populations showed evidence of admixture that likely contributed to the observed population differentiation patterns. Thus, disentangling demography from selection might be challenging when analyzing these populations. European populations could help identify loci under spatially varying selection provided that no recent admixture from African populations would have occurred. In this work, we individually sequence the genome of 42 European strains collected in populations from contrasting environments: Stockholm (Sweden) and Castellana Grotte (Southern Italy). We found low levels of population structure and no evidence of recent African admixture in these two populations. We thus look for patterns of spatially varying selection affecting individual genes and gene sets. Besides single nucleotide polymorphisms, we also investigated the role of transposable elements in local adaptation. We concluded that European populations are a good dataset to identify candidate loci under spatially varying selection. The analysis of the two populations sequenced in this work in the context of all the available D. melanogaster data allowed us to pinpoint genes and biological processes likely to be relevant for local adaptation. Identifying and analyzing populations with low levels of population structure and admixture should help to disentangle selective from non-selective forces underlying patterns of population differentiation in other species as well.
Collapse
Affiliation(s)
- Lidia Mateo
- Institute of Evolutionary Biology. CSIC-Universitat Pompeu Fabra. Passeig Maritim de la Barceloneta, 37-49. 08003, Barcelona, Spain
| | - Gabriel E Rech
- Institute of Evolutionary Biology. CSIC-Universitat Pompeu Fabra. Passeig Maritim de la Barceloneta, 37-49. 08003, Barcelona, Spain
| | - Josefa González
- Institute of Evolutionary Biology. CSIC-Universitat Pompeu Fabra. Passeig Maritim de la Barceloneta, 37-49. 08003, Barcelona, Spain.
| |
Collapse
|
6
|
Jackson BC, Campos JL, Haddrill PR, Charlesworth B, Zeng K. Variation in the Intensity of Selection on Codon Bias over Time Causes Contrasting Patterns of Base Composition Evolution in Drosophila. Genome Biol Evol 2017; 9:102-123. [PMID: 28082609 PMCID: PMC5381600 DOI: 10.1093/gbe/evw291] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/07/2016] [Indexed: 12/11/2022] Open
Abstract
Four-fold degenerate coding sites form a major component of the genome, and are often used to make inferences about selection and demography, so that understanding their evolution is important. Despite previous efforts, many questions regarding the causes of base composition changes at these sites in Drosophila remain unanswered. To shed further light on this issue, we obtained a new whole-genome polymorphism data set from D. simulans. We analyzed samples from the putatively ancestral range of D. simulans, as well as an existing polymorphism data set from an African population of D. melanogaster. By using D. yakuba as an outgroup, we found clear evidence for selection on 4-fold sites along both lineages over a substantial period, with the intensity of selection increasing with GC content. Based on an explicit model of base composition evolution, we suggest that the observed AT-biased substitution pattern in both lineages is probably due to an ancestral reduction in selection intensity, and is unlikely to be the result of an increase in mutational bias towards AT alone. By using two polymorphism-based methods for estimating selection coefficients over different timescales, we show that the selection intensity on codon usage has been rather stable in D. simulans in the recent past, but the long-term estimates in D. melanogaster are much higher than the short-term ones, indicating a continuing decline in selection intensity, to such an extent that the short-term estimates suggest that selection is only active in the most GC-rich parts of the genome. Finally, we provide evidence for complex evolutionary patterns in the putatively neutral short introns, which cannot be explained by the standard GC-biased gene conversion model. These results reveal a dynamic picture of base composition evolution.
Collapse
Affiliation(s)
- Benjamin C Jackson
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, United Kingdom
| | - José L Campos
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| | - Penelope R Haddrill
- Centre for Forensic Science, Department of Pure and Applied Chemistry, University of Strathclyde, Glasgow, United Kingdom
| | - Brian Charlesworth
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| | - Kai Zeng
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, United Kingdom
| |
Collapse
|
7
|
Cogni R, Kuczynski K, Koury S, Lavington E, Behrman EL, O’Brien KR, Schmidt PS, Eanes WF. On the Long-term Stability of Clines in Some Metabolic Genes in Drosophila melanogaster. Sci Rep 2017; 7:42766. [PMID: 28220806 PMCID: PMC5318857 DOI: 10.1038/srep42766] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2016] [Accepted: 01/13/2017] [Indexed: 11/16/2022] Open
Abstract
Very little information exists for long-term changes in genetic variation in natural populations. Here we take the unique opportunity to compare a set of data for SNPs in 15 metabolic genes from eastern US collections of Drosophila melanogaster that span a large latitudinal range and represent two collections separated by 12 to 13 years. We also expand this to a 22-year interval for the Adh gene and approximately 30 years for the G6pd and Pgd genes. During these intervals, five genes showed a statistically significant change in average SNP allele frequency corrected for latitude. While much remains unchanged, we see five genes where latitudinal clines have been lost or gained and two where the slope significantly changes. The long-term frequency shift towards a southern favored Adh S allele reported in Australia populations is not observed in the eastern US over a period of 21 years. There is no general pattern of southern-favored or northern-favored alleles increasing in frequency across the genes. This observation points to the fluid nature of some allelic variation over this time period and the action of selective responses or migration that may be more regional than uniformly imposed across the cline.
Collapse
Affiliation(s)
- Rodrigo Cogni
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York, 11794 USA
| | - Kate Kuczynski
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York, 11794 USA
| | - Spencer Koury
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York, 11794 USA
| | - Erik Lavington
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York, 11794 USA
| | - Emily L. Behrman
- Department of Biology, University of Pennsylvania, Philadelphia, PA, USA
| | | | - Paul S. Schmidt
- Department of Biology, University of Pennsylvania, Philadelphia, PA, USA
| | - Walter F. Eanes
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York, 11794 USA
| |
Collapse
|
8
|
Survey of Global Genetic Diversity Within the Drosophila Immune System. Genetics 2016; 205:353-366. [PMID: 27815361 DOI: 10.1534/genetics.116.195016] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2016] [Accepted: 10/28/2016] [Indexed: 11/18/2022] Open
Abstract
Numerous studies across a wide range of taxa have demonstrated that immune genes are routinely among the most rapidly evolving genes in the genome. This observation, however, does not address what proportion of immune genes undergo strong selection during adaptation to novel environments. Here, we determine the extent of very recent divergence in genes with immune function across five populations of Drosophila melanogaster and find that immune genes do not show an overall trend of recent rapid adaptation. Our population-based approach uses a set of carefully matched control genes to account for the effects of demography and local recombination rate, allowing us to identify whether specific immune functions are putative targets of strong selection. We find evidence that viral-defense genes are rapidly evolving in Drosophila at multiple timescales. Local adaptation to bacteria and fungi is less extreme and primarily occurs through changes in recognition and effector genes rather than large-scale changes to the regulation of the immune response. Surprisingly, genes in the Toll pathway, which show a high rate of adaptive substitution between the D. melanogaster and D. simulans lineages, show little population differentiation. Quantifying the flies for resistance to a generalist Gram-positive bacterial pathogen, we found that this genetic pattern of low population differentiation was recapitulated at the phenotypic level. In sum, our results highlight the complexity of immune evolution and suggest that Drosophila immune genes do not follow a uniform trajectory of strong directional selection as flies encounter new environments.
Collapse
|
9
|
Bergland AO, Tobler R, González J, Schmidt P, Petrov D. Secondary contact and local adaptation contribute to genome-wide patterns of clinal variation in Drosophila melanogaster. Mol Ecol 2016; 25:1157-74. [PMID: 26547394 DOI: 10.1111/mec.13455] [Citation(s) in RCA: 93] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2014] [Revised: 10/29/2015] [Accepted: 11/02/2015] [Indexed: 12/12/2022]
Abstract
Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome-wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome-wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an FST outlier approach. Our results provide novel insight into the well-studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.
Collapse
Affiliation(s)
- Alan O Bergland
- Department of Biology, Stanford University, Stanford, CA, 94305-5020, USA
| | - Ray Tobler
- Department of Biology, Stanford University, Stanford, CA, 94305-5020, USA.,Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, Vienna, A-1210, Austria
| | - Josefa González
- Department of Biology, Stanford University, Stanford, CA, 94305-5020, USA.,Institute of Evolutionary Biology (CSIC-Universitat Pompeu Fabra), Passeig Maritim de la Barceloneta 37-49, 0800, 3 Barcelona, Spain
| | - Paul Schmidt
- Department of Biology, The University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Dmitri Petrov
- Department of Biology, Stanford University, Stanford, CA, 94305-5020, USA
| |
Collapse
|
10
|
Edwards SV, Shultz AJ, Campbell-Staton SC. Next-generation sequencing and the expanding domain of phylogeography. FOLIA ZOOLOGICA 2015. [DOI: 10.25225/fozo.v64.i3.a2.2015] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
Affiliation(s)
- Scott V. Edwards
- Department of Organismic and Evolutionary Biology, and Museum of Comparative Zoology, Harvard University, Cambridge, MA 02138, U.S.A.
| | - Allison J. Shultz
- Department of Organismic and Evolutionary Biology, and Museum of Comparative Zoology, Harvard University, Cambridge, MA 02138, U.S.A.
| | - Shane C. Campbell-Staton
- Department of Organismic and Evolutionary Biology, and Museum of Comparative Zoology, Harvard University, Cambridge, MA 02138, U.S.A.
| |
Collapse
|
11
|
Kao JY, Lymer S, Hwang SH, Sung A, Nuzhdin SV. Postmating reproductive barriers contribute to the incipient sexual isolation of the United States and Caribbean Drosophila melanogaster. Ecol Evol 2015; 5:3171-82. [PMID: 26357543 PMCID: PMC4559059 DOI: 10.1002/ece3.1596] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2015] [Revised: 05/22/2015] [Accepted: 05/26/2015] [Indexed: 02/01/2023] Open
Abstract
The nascent stages of speciation start with the emergence of sexual isolation. Understanding the influence of reproductive barriers in this evolutionary process is an ongoing effort. We present a study of Drosophila melanogaster admixed populations from the southeast United States and the Caribbean islands known to be a secondary contact zone of European- and African-derived populations undergoing incipient sexual isolation. The existence of premating reproductive barriers has been previously established, but these types of barriers are not the only source shaping sexual isolation. To assess the influence of postmating barriers, we investigated putative postmating barriers of female remating and egg-laying behavior, as well as hatchability of eggs laid and female longevity after mating. In the central region of our putative hybrid zone of American and Caribbean populations, we observed lower hatchability of eggs laid accompanied by increased resistance to harm after mating to less-related males. These results illustrate that postmating reproductive barriers act alongside premating barriers and genetic admixture such as hybrid incompatibilities and influence early phases of sexual isolation.
Collapse
Affiliation(s)
- Joyce Y Kao
- Section of Molecular and Computational Biology, Department of Biology, University of Southern California Los Angeles, California, 90089 ; Department of Biology, New York University 29 Washington Pl, New York city, New York, 10003
| | - Seana Lymer
- Department of Biology, New York University 29 Washington Pl, New York city, New York, 10003
| | - Sea H Hwang
- Section of Molecular and Computational Biology, Department of Biology, University of Southern California Los Angeles, California, 90089
| | - Albert Sung
- Section of Molecular and Computational Biology, Department of Biology, University of Southern California Los Angeles, California, 90089
| | - Sergey V Nuzhdin
- Section of Molecular and Computational Biology, Department of Biology, University of Southern California Los Angeles, California, 90089 ; St. Petersburg State Polytechnical University St. Petersburg, Russia
| |
Collapse
|
12
|
Coolon JD, Stevenson KR, McManus CJ, Yang B, Graveley BR, Wittkopp PJ. Molecular Mechanisms and Evolutionary Processes Contributing to Accelerated Divergence of Gene Expression on the Drosophila X Chromosome. Mol Biol Evol 2015; 32:2605-15. [PMID: 26041937 DOI: 10.1093/molbev/msv135] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
In species with a heterogametic sex, population genetics theory predicts that DNA sequences on the X chromosome can evolve faster than comparable sequences on autosomes. Both neutral and nonneutral evolutionary processes can generate this pattern. Complex traits like gene expression are not predicted to have accelerated evolution by these theories, yet a "faster-X" pattern of gene expression divergence has recently been reported for both Drosophila and mammals. Here, we test the hypothesis that accelerated adaptive evolution of cis-regulatory sequences on the X chromosome is responsible for this pattern by comparing the relative contributions of cis- and trans-regulatory changes to patterns of faster-X expression divergence observed between strains and species of Drosophila with a range of divergence times. We find support for this hypothesis, especially among male-biased genes, when comparing different species. However, we also find evidence that trans-regulatory differences contribute to a faster-X pattern of expression divergence both within and between species. This contribution is surprising because trans-acting regulators of X-linked genes are generally assumed to be randomly distributed throughout the genome. We found, however, that X-linked transcription factors appear to preferentially regulate expression of X-linked genes, providing a potential mechanistic explanation for this result. The contribution of trans-regulatory variation to faster-X expression divergence was larger within than between species, suggesting that it is more likely to result from neutral processes than positive selection. These data show how accelerated evolution of both coding and noncoding sequences on the X chromosome can lead to accelerated expression divergence on the X chromosome relative to autosomes.
Collapse
Affiliation(s)
- Joseph D Coolon
- Department of Ecology and Evolutionary Biology, University of Michigan
| | - Kraig R Stevenson
- Department of Computational Medicine and Bioinformatics, University of Michigan
| | - C Joel McManus
- Department of Biological Sciences, Carnegie Mellon University Department of Genetics and Developmental Biology, Institute for Systems Genomics, University of Connecticut Health Center
| | - Bing Yang
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan
| | - Brenton R Graveley
- Department of Genetics and Developmental Biology, Institute for Systems Genomics, University of Connecticut Health Center
| | - Patricia J Wittkopp
- Department of Ecology and Evolutionary Biology, University of Michigan Department of Computational Medicine and Bioinformatics, University of Michigan Department of Molecular, Cellular, and Developmental Biology, University of Michigan
| |
Collapse
|
13
|
Kao JY, Zubair A, Salomon MP, Nuzhdin SV, Campo D. Population genomic analysis uncovers African and European admixture inDrosophila melanogasterpopulations from the south-eastern United States and Caribbean Islands. Mol Ecol 2015; 24:1499-509. [DOI: 10.1111/mec.13137] [Citation(s) in RCA: 65] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2014] [Revised: 02/23/2015] [Accepted: 02/25/2015] [Indexed: 01/19/2023]
Affiliation(s)
- Joyce Y. Kao
- Section of Molecular and Computational Biology; Department of Biology; University of Southern California; 1050 Childs Way Los Angeles CA 90089 USA
| | - Asif Zubair
- Section of Molecular and Computational Biology; Department of Biology; University of Southern California; 1050 Childs Way Los Angeles CA 90089 USA
| | - Matthew P. Salomon
- Section of Molecular and Computational Biology; Department of Biology; University of Southern California; 1050 Childs Way Los Angeles CA 90089 USA
| | - Sergey V. Nuzhdin
- Section of Molecular and Computational Biology; Department of Biology; University of Southern California; 1050 Childs Way Los Angeles CA 90089 USA
| | - Daniel Campo
- Section of Molecular and Computational Biology; Department of Biology; University of Southern California; 1050 Childs Way Los Angeles CA 90089 USA
| |
Collapse
|
14
|
Jackson BC, Campos JL, Zeng K. The effects of purifying selection on patterns of genetic differentiation between Drosophila melanogaster populations. Heredity (Edinb) 2014; 114:163-74. [PMID: 25227256 PMCID: PMC4270736 DOI: 10.1038/hdy.2014.80] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2014] [Revised: 06/16/2014] [Accepted: 07/22/2014] [Indexed: 01/21/2023] Open
Abstract
Using the data provided by the Drosophila Population Genomics Project, we investigate factors that affect the genetic differentiation between Rwandan and French populations of D. melanogaster. By examining within-population polymorphisms, we show that sites in long introns (especially those >2000 bp) have significantly lower π (nucleotide diversity) and more low-frequency variants (as measured by Tajima's D, minor allele frequencies, and prevalence of variants that are private to one of the two populations) than short introns, suggesting a positive relationship between intron length and selective constraint. A similar analysis of protein-coding polymorphisms shows that 0-fold (degenerate) sites in more conserved genes are under stronger purifying selection than those in less conserved genes. There is limited evidence that selection on codon bias has an effect on differentiation (as measured by FST) at 4-fold (degenerate) sites, and 4-fold sites and sites in 8–30 bp of short introns ⩽65 bp have comparable FST values. Consistent with the expected effect of purifying selection, sites in long introns and 0-fold sites in conserved genes are less differentiated than those in short introns and less conserved genes, respectively. Genes in non-crossover regions (for example, the fourth chromosome) have very high FST values at both 0-fold and 4-fold degenerate sites, which is probably because of the large reduction in within-population diversity caused by tight linkage between many selected sites. Our analyses also reveal subtle statistical properties of FST, which arise when information from multiple single nucleotide polymorphisms is combined and can lead to the masking of important signals of selection.
Collapse
Affiliation(s)
- B C Jackson
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK
| | - J L Campos
- Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, UK
| | - K Zeng
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK
| |
Collapse
|
15
|
Blumenstiel JP, Chen X, He M, Bergman CM. An age-of-allele test of neutrality for transposable element insertions. Genetics 2014; 196:523-38. [PMID: 24336751 PMCID: PMC3914624 DOI: 10.1534/genetics.113.158147] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2013] [Accepted: 12/06/2013] [Indexed: 01/31/2023] Open
Abstract
How natural selection acts to limit the proliferation of transposable elements (TEs) in genomes has been of interest to evolutionary biologists for many years. To describe TE dynamics in populations, previous studies have used models of transposition-selection equilibrium that assume a constant rate of transposition. However, since TE invasions are known to happen in bursts through time, this assumption may not be reasonable. Here we propose a test of neutrality for TE insertions that does not rely on the assumption of a constant transposition rate. We consider the case of TE insertions that have been ascertained from a single haploid reference genome sequence. By conditioning on the age of an individual TE insertion allele (inferred by the number of unique substitutions that have occurred within the particular TE sequence since insertion), we determine the probability distribution of the insertion allele frequency in a population sample under neutrality. Taking models of varying population size into account, we then evaluate predictions of our model against allele frequency data from 190 retrotransposon insertions sampled from North American and African populations of Drosophila melanogaster. Using this nonequilibrium neutral model, we are able to explain ∼ 80% of the variance in TE insertion allele frequencies based on age alone. Controlling for both nonequilibrium dynamics of transposition and host demography, we provide evidence for negative selection acting against most TEs as well as for positive selection acting on a small subset of TEs. Our work establishes a new framework for the analysis of the evolutionary forces governing large insertion mutations like TEs, gene duplications, or other copy number variants.
Collapse
Affiliation(s)
- Justin P. Blumenstiel
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas 66049
| | - Xi Chen
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas 66049
| | - Miaomiao He
- Faculty of Life Sciences, University of Manchester, Manchester M21 0RG, United Kingdom
| | - Casey M. Bergman
- Faculty of Life Sciences, University of Manchester, Manchester M21 0RG, United Kingdom
| |
Collapse
|
16
|
Cogni R, Kuczynski C, Koury S, Lavington E, Behrman EL, O'Brien KR, Schmidt PS, Eanes WF. THE INTENSITY OF SELECTION ACTING ON THECOUCH POTATOGENE-SPATIAL-TEMPORAL VARIATION IN A DIAPAUSE CLINE. Evolution 2013; 68:538-48. [DOI: 10.1111/evo.12291] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2013] [Accepted: 09/26/2013] [Indexed: 11/27/2022]
Affiliation(s)
- Rodrigo Cogni
- Department of Ecology and Evolution; Stony Brook University; Stony Brook New York
| | - Caitlin Kuczynski
- Department of Ecology and Evolution; Stony Brook University; Stony Brook New York
| | - Spencer Koury
- Department of Ecology and Evolution; Stony Brook University; Stony Brook New York
| | - Erik Lavington
- Department of Ecology and Evolution; Stony Brook University; Stony Brook New York
| | - Emily L. Behrman
- Department of Biology; University of Pennsylvania; Philadelphia Pennsylvania
| | | | - Paul S. Schmidt
- Department of Biology; University of Pennsylvania; Philadelphia Pennsylvania
| | - Walter F. Eanes
- Department of Ecology and Evolution; Stony Brook University; Stony Brook New York
| |
Collapse
|
17
|
Campo D, Lehmann K, Fjeldsted C, Souaiaia T, Kao J, Nuzhdin SV. Whole-genome sequencing of two North American Drosophila melanogaster populations reveals genetic differentiation and positive selection. Mol Ecol 2013; 22:5084-97. [PMID: 24102956 DOI: 10.1111/mec.12468] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2012] [Revised: 07/15/2013] [Accepted: 07/16/2013] [Indexed: 11/29/2022]
Abstract
The prevailing demographic model for Drosophila melanogaster suggests that the colonization of North America occurred very recently from a subset of European flies that rapidly expanded across the continent. This model implies a sudden population growth and range expansion consistent with very low or no population subdivision. As flies adapt to new environments, local adaptation events may be expected. To describe demographic and selective events during North American colonization, we have generated a data set of 35 individual whole-genome sequences from inbred lines of D. melanogaster from a west coast US population (Winters, California, USA) and compared them with a public genome data set from Raleigh (Raleigh, North Carolina, USA). We analysed nuclear and mitochondrial genomes and described levels of variation and divergence within and between these two North American D. melanogaster populations. Both populations exhibit negative values of Tajima's D across the genome, a common signature of demographic expansion. We also detected a low but significant level of genome-wide differentiation between the two populations, as well as multiple allele surfing events, which can be the result of gene drift in local subpopulations on the edge of an expansion wave. In contrast to this genome-wide pattern, we uncovered a 50-kilobase segment in chromosome arm 3L that showed all the hallmarks of a soft selective sweep in both populations. A comparison of allele frequencies within this divergent region among six populations from three continents allowed us to cluster these populations in two differentiated groups, providing evidence for the action of natural selection on a global scale.
Collapse
Affiliation(s)
- D Campo
- Molecular and Computational Biology, University of Southern California, Los Angeles, CA, 90089, USA
| | | | | | | | | | | |
Collapse
|
18
|
Duchen P, Zivkovic D, Hutter S, Stephan W, Laurent S. Demographic inference reveals African and European admixture in the North American Drosophila melanogaster population. Genetics 2013; 193:291-301. [PMID: 23150605 PMCID: PMC3527251 DOI: 10.1534/genetics.112.145912] [Citation(s) in RCA: 117] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2012] [Accepted: 10/30/2012] [Indexed: 11/18/2022] Open
Abstract
Drosophila melanogaster spread from sub-Saharan Africa to the rest of the world colonizing new environments. Here, we modeled the joint demography of African (Zimbabwe), European (The Netherlands), and North American (North Carolina) populations using an approximate Bayesian computation (ABC) approach. By testing different models (including scenarios with continuous migration), we found that admixture between Africa and Europe most likely generated the North American population, with an estimated proportion of African ancestry of 15%. We also revisited the demography of the ancestral population (Africa) and found-in contrast to previous work-that a bottleneck fits the history of the population of Zimbabwe better than expansion. Finally, we compared the site-frequency spectrum of the ancestral population to analytical predictions under the estimated bottleneck model.
Collapse
Affiliation(s)
- Pablo Duchen
- Evolutionary Biology, University of Munich, 82152 Planegg-Martinsried, Germany.
| | | | | | | | | |
Collapse
|
19
|
Fabian DK, Kapun M, Nolte V, Kofler R, Schmidt PS, Schlötterer C, Flatt T. Genome-wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America. Mol Ecol 2012; 21:4748-69. [PMID: 22913798 PMCID: PMC3482935 DOI: 10.1111/j.1365-294x.2012.05731.x] [Citation(s) in RCA: 177] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2012] [Revised: 06/29/2012] [Accepted: 07/04/2012] [Indexed: 01/18/2023]
Abstract
Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome-wide next-generation sequencing of DNA pools ('pool-seq') from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome-wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well-known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.
Collapse
Affiliation(s)
- Daniel K Fabian
- Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, A-1210, Vienna, Austria
| | | | | | | | | | | | | |
Collapse
|
20
|
Fang S, Yukilevich R, Chen Y, Turissini DA, Zeng K, Boussy IA, Wu CI. Incompatibility and competitive exclusion of genomic segments between sibling Drosophila species. PLoS Genet 2012; 8:e1002795. [PMID: 22761593 PMCID: PMC3386244 DOI: 10.1371/journal.pgen.1002795] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2011] [Accepted: 05/15/2012] [Indexed: 11/22/2022] Open
Abstract
The extent and nature of genetic incompatibilities between incipient races and sibling species is of fundamental importance to our view of speciation. However, with the exception of hybrid inviability and sterility factors, little is known about the extent of other, more subtle genetic incompatibilities between incipient species. Here we experimentally demonstrate the prevalence of such genetic incompatibilities between two young allopatric sibling species, Drosophila simulans and D. sechellia. Our experiments took advantage of 12 introgression lines that carried random introgressed D. sechellia segments in different parts of the D. simulans genome. First, we found that these introgression lines did not show any measurable sterility or inviability effects. To study if these sechellia introgressions in a simulans background contained other fitness consequences, we competed and genetically tracked the marked alleles within each introgression against the wild-type alleles for 20 generations. Strikingly, all marked D. sechellia introgression alleles rapidly decreased in frequency in only 6 to 7 generations. We then developed computer simulations to model our competition results. These simulations indicated that selection against D. sechellia introgression alleles was high (average s = 0.43) and that the marker alleles and the incompatible alleles did not separate in 78% of the introgressions. The latter result likely implies that most introgressions contain multiple genetic incompatibilities. Thus, this study reveals that, even at early stages of speciation, many parts of the genome diverge to a point where introducing foreign elements has detrimental fitness consequences, but which cannot be seen using standard sterility and inviability assays. Determining the extent of genomic incompatibilities is a pivotal issue in understanding the process of speciation. A controversial topic that has recently sparked debate is whether there are few isolated genetic regions (so-called “genomic islands of speciation”) or extensive genetic regions (“genomic continents of speciation”) responsible for species divergence. To answer this question, most work has focused on species divergence at the DNA sequence level. Here, we present a new perspective by shifting the focus to the fitness and functional aspects of foreign genomic introgression. To illustrate our point, we performed an introgression experiment on two sibling species, D. sechellia and D. simulans. After introgressing random genomic segments of D. sechellia into D. simulans genetic background, a 20-generation competition experiment revealed that, even at the early stages of speciation, there are virtually always detrimental fitness consequences to introducing random foreign elements from one genome to another. This implies that incipient speciation may be characterized by widespread accumulation of genomic incompatibilities rather than a few isolated genes. This study shows that we should move beyond the sterility and inviability assays in order to understand the full extent of genetic incompatibilities during speciation.
Collapse
Affiliation(s)
- Shu Fang
- Biodiversity Research Center, Academia Sinica, Taipei, Taiwan, Republic of China
- * E-mail: (SF); (RY)
| | - Roman Yukilevich
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois, United States of America
- Biology Department, Union College, Schenectady, New York, United States of America
- * E-mail: (SF); (RY)
| | - Ying Chen
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois, United States of America
- Department of Biology, Loyola University Chicago, Chicago, Illinois, United States of America
| | - David A. Turissini
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois, United States of America
| | - Kai Zeng
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois, United States of America
- State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen (Zhongshan) University, Guangzhou, People's Republic of China
| | - Ian A. Boussy
- Department of Biology, Loyola University Chicago, Chicago, Illinois, United States of America
| | - Chung-I. Wu
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois, United States of America
- State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen (Zhongshan) University, Guangzhou, People's Republic of China
- Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, People's Republic of China
| |
Collapse
|
21
|
Pavlidis P, Jensen JD, Stephan W, Stamatakis A. A critical assessment of storytelling: gene ontology categories and the importance of validating genomic scans. Mol Biol Evol 2012; 29:3237-48. [PMID: 22617950 DOI: 10.1093/molbev/mss136] [Citation(s) in RCA: 159] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
In the age of whole-genome population genetics, so-called genomic scan studies often conclude with a long list of putatively selected loci. These lists are then further scrutinized to annotate these regions by gene function, corresponding biological processes, expression levels, or gene networks. Such annotations are often used to assess and/or verify the validity of the genome scan and the statistical methods that have been used to perform the analyses. Furthermore, these results are frequently considered to validate "true-positives" if the identified regions make biological sense a posteriori. Here, we show that this approach can be potentially misleading. By simulating neutral evolutionary histories, we demonstrate that it is possible not only to obtain an extremely high false-positive rate but also to make biological sense out of the false-positives and construct a sensible biological narrative. Results are compared with a recent polymorphism data set from Drosophila melanogaster.
Collapse
Affiliation(s)
- Pavlos Pavlidis
- The Exelixis Lab, Scientific Computing Group, Heidelberg Institute for Theoretical Studies (HITS gGmbH), Heidelberg, Germany.
| | | | | | | |
Collapse
|