1
|
Tseng SP, Wetterer JK, Suarez AV, Lee CY, Yoshimura T, Shoemaker D, Yang CCS. Genetic Diversity and Wolbachia Infection Patterns in a Globally Distributed Invasive Ant. Front Genet 2019; 10:838. [PMID: 31608104 PMCID: PMC6758599 DOI: 10.3389/fgene.2019.00838] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2019] [Accepted: 08/13/2019] [Indexed: 11/23/2022] Open
Abstract
Understanding the phylogeographic history of an invasive species may facilitate reconstructing the history and routes of its invasion. The longhorn crazy ant, Paratrechina longicornis, is a ubiquitous agricultural and household pest throughout much of the tropics and subtropics, but little is known about the history of its spread. Here, we examine worldwide genetic variation in P. longicornis and its associated Wolbachia bacterial symbionts. Analyses of mtDNA sequences of 248 P. longicornis workers (one per colony) from 13 geographic regions reveal two highly diverged mtDNA clades that co-occur in most of the geographic regions. These two mtDNA clades are associated with different Wolbachia infection patterns, but are not congruent with patterns of nDNA (microsatellite) variation. Multilocus sequence typing reveals two distinct Wolbachia strains in P. longicornis, namely, wLonA and wLonF. The evolutionary histories of these two strains differ; wLonA appears to be primarily transmitted maternally, and patterns of mtDNA and nDNA variation and wLonA infection status are consistent with a relatively recent Wolbachia-induced selective sweep. In contrast, the observed patterns of mtDNA variation and wLonF infections suggest frequent horizontal transfer and losses of wLonF infections. The lack of nDNA structure among sampled geographic regions coupled with the finding that numerous mtDNA haplotypes are shared among regions implies that inadvertent long-distance movement through human commerce is common in P. longicornis and has shaped the genetic structure of this invasive ant worldwide.
Collapse
Affiliation(s)
- Shu-Ping Tseng
- Research Institute for Sustainable Humanosphere, Kyoto University, Kyoto, Japan
| | - James K. Wetterer
- Wilkes Honors College, Florida Atlantic University, Jupiter, FL, United States
| | - Andrew V. Suarez
- Department of Evolution, Ecology and Behavior and Department of Entomology, University of Illinois Urbana-Champaign, IL, United States
| | - Chow-Yang Lee
- Department of Entomology, University of California, Riverside, CA, United States
| | - Tsuyoshi Yoshimura
- Research Institute for Sustainable Humanosphere, Kyoto University, Kyoto, Japan
| | - DeWayne Shoemaker
- Entomology and Plant Pathology, University of Tennessee, Knoxville, TN, United States
| | | |
Collapse
|
2
|
Bazykin GA. Changing preferences: deformation of single position amino acid fitness landscapes and evolution of proteins. Biol Lett 2016; 11:rsbl.2015.0315. [PMID: 26445980 DOI: 10.1098/rsbl.2015.0315] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
The fitness landscape-the function that relates genotypes to fitness-and its role in directing evolution are a central object of evolutionary biology. However, its huge dimensionality precludes understanding of even the basic aspects of its shape. One way to approach it is to ask a simpler question: what are the properties of a function that assigns fitness to each possible variant at just one particular site-a single position fitness landscape-and how does it change in the course of evolution? Analyses of genomic data from multiple species and multiple individuals within a species have proved beyond reasonable doubt that fitness functions of positions throughout the genome do themselves change with time, thus shaping protein evolution. Here, I will briefly review the literature that addresses these dynamics, focusing on recent genome-scale analyses of fitness functions of amino acid sites, i.e. vectors of fitnesses of 20 individual amino acid variants at a given position of a protein. The set of amino acids that confer high fitness at a particular position changes with time, and the rate of this change is comparable with the rate at which a position evolves, implying that this process plays a major role in evolutionary dynamics. However, the causes of these changes remain largely unclear.
Collapse
Affiliation(s)
- Georgii A Bazykin
- Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences, Moscow 127051, Russia Faculty of Bioengineering and Bioinformatics and Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow 119234, Russia Pirogov Russian National Research Medical University, Moscow 117997, Russia
| |
Collapse
|
3
|
Vinogradov AE. Accelerated pathway evolution in mouse-like rodents involves cell cycle control. Mamm Genome 2015; 26:609-18. [PMID: 26424469 DOI: 10.1007/s00335-015-9605-8] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2015] [Accepted: 09/23/2015] [Indexed: 10/23/2022]
Abstract
Rodents include both the cancer-susceptible short-lived mouse and the two unrelated cancer-resistant long-lived mole-rats. In this work, their genomes were analyzed with the goal to reveal pathways enriched in genes, which are more similar between the mole-rats than between the mouse and the naked mole-rat. The pathways related to cell cycle control were prominent. They include external signal transduction and all cell cycle stages. There are several stem cell pathways among them. The other enriched pathways involve ubiquitin-dependent protein degradation, immunity, mRNA splicing, and apoptosis. The ubiquitin-dependent protein degradation is a core of network of enriched pathways. However, this phenomenon is not specific for the mouse and the mole-rats. The other muroid species show features similar to the mouse, whereas the non-muroid rodents and the human show features similar to the mole-rats. The higher ratio of non-synonymous to synonymous nucleotide substitutions (dN/dS) indicates the accelerated evolution of revealed pathways in the muroid rodents (except the blind mole-rat). Paradoxically, the dN/dS averaged over the whole genome is lower in the muroids, i.e., the purifying selection is generally stronger in them. In practical sense, these data suggest caveat for using muroid rodents (mouse, rat, and hamsters) as biomedical models of human conditions involving cell cycle and show the network of pathways where muroid genes are most different (compared with non-muroid) from human genes. The guinea pig is emphasized as a more suitable rodent model for biomedical research involving cell cycle.
Collapse
|
4
|
Malyarchuk B, Derenko M, Denisova G. A mitogenomic phylogeny and genetic history of sable (Martes zibellina). Gene 2014; 550:56-67. [PMID: 25110108 DOI: 10.1016/j.gene.2014.08.015] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2014] [Revised: 07/29/2014] [Accepted: 08/06/2014] [Indexed: 12/21/2022]
Abstract
We assessed phylogeny of sable (Martes zibellina, Linnaeus, 1758) by sequence analysis of nearly complete, new mitochondrial genomes in 36 specimens from different localities in northern Eurasia (Primorye, Khabarovsk and Krasnoyarsk regions, the Kamchatka Peninsula, the Kuril Islands and the Urals). Phylogenetic analysis of mtDNA sequences demonstrates that two clades, A and BC, radiated about 200-300 thousandyears ago (kya) according to results of Bayesian molecular clock and RelTime analyses of different mitogenome alignments (nearly complete mtDNA sequences, protein-coding region, and synonymous sites), while the age estimates of clades A, B and C fall within the Late Pleistocene (~50-140 kya). Bayesian skyline plots (BSPs) of sable population size change based on analysis of nearly complete mtDNAs show an expansion around 40 kya in the warm Karganian time, without a decline of population size around the Last Glacial Maximum (21 kya). The BSPs based on synonymous clock rate indicate that M. zibellina experienced demographic expansions later, approximately 22 kya. The A2a clade that colonized Kamchatka ~23-50 kya (depending on the mutation rate used) survived the last glaciation there as demonstrated by the BSP analysis. In addition, we have found evidence of positive selection acting at ND4 and cytochrome b genes, thereby suggesting adaptive evolution of the A2a clade in Kamchatka.
Collapse
Affiliation(s)
- Boris Malyarchuk
- Institute of Biological Problems of the North, Magadan, 685000 Russia.
| | - Miroslava Derenko
- Institute of Biological Problems of the North, Magadan, 685000 Russia
| | - Galina Denisova
- Institute of Biological Problems of the North, Magadan, 685000 Russia
| |
Collapse
|
5
|
Kamaraj B, Kumar A, Purohit R. Evolutionary reconstruction and population genetics analysis of aurora kinases. PLoS One 2013; 8:e75763. [PMID: 24086628 PMCID: PMC3782425 DOI: 10.1371/journal.pone.0075763] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2013] [Accepted: 08/16/2013] [Indexed: 12/03/2022] Open
Abstract
Background Aurora kinases belong to the highly conserved kinase family and play a vital role in cell cycle regulation. The structure and function of these kinases are inter-related and sometimes they also act as substitutes in case of knockdown of other aurora kinases. Method In this work we carried out the evolutionary reconstruction and population genetic studies of aurora kinase proteins. Substitution saturation test, CAI (Codon adaptation index), gene expression and RSCU (Relative synonymous codon usage) values were computed for all the three aurora kinases. Linear regression method was used to check the dependency of gene expression on their CAI values. Results The results suggested that aurora-B and aurora-C has shown convergence in their evolutionary pathway. Moreover, the aurora-A I57V mutation showed high penetrance in human population and exist at very high frequency (84.4%) when compared to the native residue (15.6%). The mutation showed notable range of functional gain and seemed to be promising for the evolution of aurora-A function. Mutant allele might also become a challenging prospect for understanding the pattern of evolution followed by cell cycle kinases. Conclusion The overall result suggested that the aurora-A is currently under the evolutionary transition and to determine the functional significance of the mutation further investigation are required.
Collapse
Affiliation(s)
- Balu Kamaraj
- Bioinformatics Division, School of Bio Sciences and Technology, Vellore Institute of Technology University, Vellore, Tamil Nadu, India
| | - Ambuj Kumar
- Bioinformatics Division, School of Bio Sciences and Technology, Vellore Institute of Technology University, Vellore, Tamil Nadu, India
| | - Rituraj Purohit
- Bioinformatics Division, School of Bio Sciences and Technology, Vellore Institute of Technology University, Vellore, Tamil Nadu, India
- * E-mail:
| |
Collapse
|
6
|
Abstract
The most common models of sequence evolution used to make inferences about adaptation rely on the assumption that selective pressures at a site remain constant through time. Instead, one might plausibly imagine that a change in the environment renders an allele beneficial and that when it fixes, the site is now constrained-until another change in the environment occurs that affects the selective pressures at that site. With this view in mind, we introduce a simple dynamic model for the evolution of coding regions, in which non-synonymous sites alternate between being fixed for the favored allele and being neutral with respect to other alleles. We use the pruning algorithm to derive closed forms for observable patterns of polymorphism and divergence in terms of the model parameters. Using our model, estimates of the fraction of beneficial substitutions α would remain similar to those obtained from existing approaches. In this framework, however, it becomes natural to ask how often adaptive substitutions originate from previously constrained or previously neutral sites, i.e., about the source of adaptive substitutions. We show that counts of coding sites that are both polymorphic in a sample from one species and divergent between two others carry information about this parameter. We also extend the basic model to include the effects of weakly deleterious mutations and discuss the importance of assumptions about the distribution of deleterious mutations among constrained non-synonymous sites. Finally, we derive a likelihood function for the parameters and apply it to a toy example, variation data for coding regions from chromosome 2 of the Drosophila melanogaster subgroup. This modeling work underscores how restrictive assumptions about adaptation have been to date, and how further work in this area will help to reveal unexplored and yet basic characteristics of adaptation.
Collapse
|
7
|
Leffler EM, Bullaughey K, Matute DR, Meyer WK, Ségurel L, Venkat A, Andolfatto P, Przeworski M. Revisiting an old riddle: what determines genetic diversity levels within species? PLoS Biol 2012; 10:e1001388. [PMID: 22984349 PMCID: PMC3439417 DOI: 10.1371/journal.pbio.1001388] [Citation(s) in RCA: 330] [Impact Index Per Article: 27.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
Understanding why some species have more genetic diversity than others is central to the study of ecology and evolution, and carries potentially important implications for conservation biology. Yet not only does this question remain unresolved, it has largely fallen into disregard. With the rapid decrease in sequencing costs, we argue that it is time to revive it.
Collapse
Affiliation(s)
- Ellen M. Leffler
- Department of Human Genetics, University of Chicago, Chicago, Illinois, United States of America
- * E-mail: (EML); (MP)
| | - Kevin Bullaughey
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois, United States of America
| | - Daniel R. Matute
- Department of Human Genetics, University of Chicago, Chicago, Illinois, United States of America
| | - Wynn K. Meyer
- Department of Human Genetics, University of Chicago, Chicago, Illinois, United States of America
| | - Laure Ségurel
- Department of Human Genetics, University of Chicago, Chicago, Illinois, United States of America
- Howard Hughes Medical Institute, University of Chicago, Chicago, Illinois, United States of America
| | - Aarti Venkat
- Department of Human Genetics, University of Chicago, Chicago, Illinois, United States of America
| | - Peter Andolfatto
- Department of Ecology and Evolutionary Biology and the Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, New Jersey, United States of America
| | - Molly Przeworski
- Department of Human Genetics, University of Chicago, Chicago, Illinois, United States of America
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois, United States of America
- Howard Hughes Medical Institute, University of Chicago, Chicago, Illinois, United States of America
- * E-mail: (EML); (MP)
| |
Collapse
|
8
|
Suzuki Y. Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus. Genes Genet Syst 2012; 86:287-94. [PMID: 22362027 DOI: 10.1266/ggs.86.287] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
The number of N-linked glycosylation sites in the globular head of hemagglutinin (HA) has increased during evolution of H3N2 human influenza A virus. Here natural selection operating on the gains of N-linked glycosylation sites was examined by using the single-site analysis and the single-substitution analysis. In the single-site analysis, positive selection was not inferred at the amino acid sites where the substitutions generating N-linked glycosylation sites were observed, but was detected at antigenic sites. In contrast, in the single-substitution analysis, positive selection was detected for the amino acid substitutions generating N-linked glycosylation sites. The single-site analysis and the single-substitution analysis appeared to be suitable for detecting recurrent and episodic natural selection, respectively. The gains of N-linked glycosylation sites were likely to be positively selected for the function of shielding antigenic sites from immune responses. At the antigenic sites, positive selection appeared to have operated not only on the radical substitution but also on the conservative substitution in terms of the charge of amino acids, suggesting that the antigenic drift is not a by-product of the evolution of receptor binding avidity in HA of human H3N2 virus.
Collapse
Affiliation(s)
- Yoshiyuki Suzuki
- Graduate School of Natural Sciences, Nagoya City University, Aichi-ken, Japan.
| |
Collapse
|
9
|
Connallon T, Singh ND, Clark AG. Impact of genetic architecture on the relative rates of X versus autosomal adaptive substitution. Mol Biol Evol 2012; 29:1933-42. [PMID: 22319138 DOI: 10.1093/molbev/mss057] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Molecular evolutionary theory predicts that the ratio of autosomal to X-linked adaptive substitution (K(A)/K(x)) is primarily determined by the average dominance coefficient of beneficial mutations. Although this theory has profoundly influenced analysis and interpretation of comparative genomic data, its predictions are based upon two unverified assumptions about the genetic basis of adaptation. The theory assumes that 1) the rate of adaptively driven molecular evolution is limited by the availability of beneficial mutations, and 2) the scaling of evolutionary parameters between the X and the autosomes (e.g., the beneficial mutation rate, and the fitness effect distribution of beneficial alleles, per X-linked versus autosomal locus) is constant across molecular evolutionary timescales. Here, we show that the genetic architecture underlying bouts of adaptive substitution can influence both assumptions, and consequently, the theoretical relationship between K(A)/K(x) and mean dominance. Quantitative predictions of prior theory apply when 1) many genomically dispersed genes potentially contribute beneficial substitutions during individual steps of adaptive walks, and 2) the population beneficial mutation rate, summed across the set of potentially contributing genes, is sufficiently small to ensure that adaptive substitutions are drawn from new mutations rather than standing genetic variation. Current research into the genetic basis of adaptation suggests that both assumptions are plausibly violated. We find that the qualitative positive relationship between mean dominance and K(A)/K(x) is relatively robust to the specific conditions underlying adaptive substitution, yet the quantitative relationship between dominance and K(A)/K(x) is quite flexible and context dependent. This flexibility may partially account for the puzzlingly variable X versus autosome substitution patterns reported in the empirical evolutionary genomics literature. The new theory unites the previously separate analysis of adaptation using new mutations versus standing genetic variation and makes several useful predictions about the interaction between genetic architecture, evolutionary genetic constraints, and effective population size in determining the ratio of adaptive substitution between autosomal and X-linked genes.
Collapse
Affiliation(s)
- Tim Connallon
- Department of Molecular Biology and Genetics, Cornell University, USA.
| | | | | |
Collapse
|