Detecting past positive selection through ongoing negative selection

Genetic Diversity and Wolbachia Infection Patterns in a Globally Distributed Invasive Ant

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.

Changing preferences: deformation of single position amino acid fitness landscapes and evolution of proteins

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.

Accelerated pathway evolution in mouse-like rodents involves cell cycle control

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.

A mitogenomic phylogeny and genetic history of sable (Martes zibellina)

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.

Evolutionary reconstruction and population genetics analysis of aurora kinases

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.

Modeling the effect of changing selective pressures on polymorphism and divergence

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.

Revisiting an old riddle: what determines genetic diversity levels within species?

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.

Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus

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.

Impact of genetic architecture on the relative rates of X versus autosomal adaptive substitution

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.