Inferring stabilizing mutations from protein phylogenies : application to influenza hemagglutinin

One selection pressure shaping sequence evolution is the requirement that a protein fold with sufficient stability to perform its biological functions. We present a conceptual framework that explains how this requirement causes the probability that a particular amino acid mutation is fixed during evolution to depend on its effect on protein stability. We mathematically formalize this framework to develop a Bayesian approach for inferring the stability effects of individual mutations from homologous protein sequences of known phylogeny. This approach is able to predict published experimentally measured mutational stability effects (ΔΔG values) with an accuracy that exceeds both a state-of-the-art physicochemical modeling program and the sequence-based consensus approach. As a further test, we use our phylogenetic inference approach to predict stabilizing mutations to influenza hemagglutinin. We introduce these mutations into a temperature-sensitive influenza virus with a defect in its hemagglutinin gene and experimentally demonstrate that some of the mutations allow the virus to grow at higher temperatures. Our work therefore describes a powerful new approach for predicting stabilizing mutations that can be successfully applied even to large, complex proteins such as hemagglutinin. This approach also makes a mathematical link between phylogenetics and experimentally measurable protein properties, potentially paving the way for more accurate analyses of molecular evolution

Diversifying Selection Underlies the Origin of Allozyme Polymorphism at the Phosphoglucose Isomerase Locus in Tigriopus californicus

The marine copepod Tigriopus californicus lives in intertidal rock pools along the Pacific coast, where it exhibits strong, temporally stable population genetic structure. Previous allozyme surveys have found high frequency private alleles among neighboring subpopulations, indicating that there is limited genetic exchange between populations. Here we evaluate the factors responsible for the diversification and maintenance of alleles at the phosphoglucose isomerase (Pgi) locus by evaluating patterns of nucleotide variation underlying previously identified allozyme polymorphism. Copepods were sampled from eleven sites throughout California and Baja California, revealing deep genetic structure among populations as well as genetic variability within populations. Evidence of recombination is limited to the sample from Pescadero and there is no support for linkage disequilibrium across the Pgi locus. Neutrality tests and codon-based models of substitution suggest the action of natural selection due to elevated non-synonymous substitutions at a small number of sites in Pgi. Two sites are identified as the charge-changing residues underlying allozyme polymorphisms in T. californicus. A reanalysis of allozyme variation at several focal populations, spanning a period of 26 years and over 200 generations, shows that Pgi alleles are maintained without notable frequency changes. Our data suggest that diversifying selection accounted for the origin of Pgi allozymes, while McDonald-Kreitman tests and the temporal stability of private allozyme alleles suggests that balancing selection may be involved in the maintenance of amino acid polymorphisms within populations

Molecular Adaptation of rbcL in the Heterophyllous Aquatic Plant Potamogeton

Heterophyllous aquatic plants show marked phenotypic plasticity. They adapt to environmental changes by producing different leaf types: submerged, floating and terrestrial leaves. By contrast, homophyllous plants produce only submerged leaves and grow entirely underwater. Heterophylly and submerged homophylly evolved under selective pressure modifying the species-specific optima for photosynthesis, but little is known about the evolutionary outcome of habit. Recent evolutionary analyses suggested that rbcL, a chloroplast gene that encodes a catalytic subunit of RuBisCO, evolves under positive selection in most land plant lineages. To examine the adaptive evolutionary process linked to heterophylly or homophylly, we analyzed positive selection in the rbcL sequences of ecologically diverse aquatic plants, Japanese Potamogeton.Phylogenetic and maximum likelihood analyses of codon substitution models indicated that Potamogeton rbcL has evolved under positive Darwinian selection. The positive selection has operated specifically in heterophyllous lineages but not in homophyllous ones in the branch-site models. This suggests that the selective pressure on this chloroplast gene was higher for heterophyllous lineages than for homophyllous lineages. The replacement of 12 amino acids occurred at structurally important sites in the quaternary structure of RbcL, two of which (residue 225 and 281) were identified as potentially under positive selection.Our analysis did not show an exact relationship between the amino acid replacements and heterophylly or homophylly but revealed that lineage-specific positive selection acted on the Potamogeton rbcL. The contrasting ecological conditions between heterophyllous and homophyllous plants have imposed different selective pressures on the photosynthetic system. The increased amino acid replacement in RbcL may reflect the continuous fine-tuning of RuBisCO under varying ecological conditions

Adaptive Evolution in Zinc Finger Transcription Factors

The majority of human genes are conserved among mammals, but some gene families have undergone extensive expansion in particular lineages. Here, we present an evolutionary analysis of one such gene family, the poly–zinc-finger (poly-ZF) genes. The human genome encodes approximately 700 members of the poly-ZF family of putative transcriptional repressors, many of which have associated KRAB, SCAN, or BTB domains. Analysis of the gene family across the tree of life indicates that the gene family arose from a small ancestral group of eukaryotic zinc-finger transcription factors through many repeated gene duplications accompanied by functional divergence. The ancestral gene family has probably expanded independently in several lineages, including mammals and some fishes. Investigation of adaptive evolution among recent paralogs using dN/dS analysis indicates that a major component of the selective pressure acting on these genes has been positive selection to change their DNA-binding specificity. These results suggest that the poly-ZF genes are a major source of new transcriptional repression activity in humans and other primates

Evolutionary patterns of two major reproduction candidate genes (Zp2 and Zp3) reveal no contribution to reproductive isolation between bovine species

<p>Abstract</p> <p>Background</p> <p>It has been established that mammalian egg zona pellucida (ZP) glycoproteins are responsible for species-restricted binding of sperm to unfertilized eggs, inducing the sperm acrosome reaction, and preventing polyspermy. In mammals, ZP apparently represents a barrier to heterospecific fertilization and thus probably contributes to reproductive isolation between species. The evolutionary relationships between some members of the tribe Bovini are complex and highly debatable, particularly, those involving <it>Bos </it>and <it>Bison </it>species for which interspecific hybridization is extensively documented. Because reproductive isolation is known to be a major precursor of species divergence, testing evolutionary patterns of ZP glycoproteins may shed some light into the speciation process of these species. To this end, we have examined intraspecific and interspecific genetic variation of two ZP genes (<it>Zp2 </it>and <it>Zp3</it>) for seven representative species (111 individuals) from the Bovini tribe, including five species from <it>Bos </it>and <it>Bison</it>, and two species each from genera <it>Bubalus </it>and <it>Syncerus</it>.</p> <p>Results</p> <p>A pattern of low levels of intraspecific polymorphism and interspecific divergence was detected for the two sequenced fragments each for <it>Zp2 </it>and <it>Zp3</it>. At intraspecific level, none of neutrality tests detected deviations from neutral equilibrium expectations for the two genes. Several haplotypes in both genes were shared by multiple species from <it>Bos </it>and <it>Bison</it>.</p> <p>Conclusions</p> <p>Here we argue that neither ancestral polymorphism nor introgressive hybridization alone can fully account for haplotype sharing among species from <it>Bos </it>and <it>Bison</it>, and that both scenarios have contributed to such a pattern of haplotype sharing observed here. Additionally, codon-based tests revealed strong evidence for purifying selection in the <it>Zp3 </it>coding haplotype sequences and weak evidence for purifying selection in the <it>Zp2 </it>coding haplotype sequences. Contrary to a general genetic pattern that genes or genomic regions contributing to reproductive isolation between species often evolve rapidly and show little or no gene flow between species, these results demonstrate that, particularly, those sequenced exons of the <it>Zp2 </it>and the <it>Zp3 </it>did not show any contribution to reproductive isolation between the bovine species studied here.</p

Evidence for maintenance of sex determinants but not of sexual stages in red yeasts, a group of early diverged basidiomycetes

<p>Abstract</p> <p>Background</p> <p>The red yeasts are an early diverged group of basidiomycetes comprising sexual and asexual species. Sexuality is based on two compatible mating types and sexual identity is determined by <it>MAT </it>loci that encode homeodomain transcription factors, peptide pheromones and their receptors. The objective of the present study was to investigate the presence and integrity of <it>MAT </it>genes throughout the phylogenetic diversity of red yeasts belonging to the order Sporidiobolales.</p> <p>Results</p> <p>We surveyed 18 sexual heterothallic and self-fertile species and 16 asexual species. Functional pheromone receptor homologues (<it>STE3.A1 </it>and <it>STE3.A2</it>) were found in multiple isolates of most of the sexual and asexual species. For each of the two mating types, sequence comparisons with whole-genome data indicated that synteny tended to be conserved along the pheromone receptor region. For the homeodomain transcription factor, likelihood methods suggested that diversifying selection acting on the self/non-self recognition region promotes diversity in sexual species, while rapid evolution seems to be due to relaxed selection in asexual strains.</p> <p>Conclusions</p> <p>The majority of both sexual and asexual species of red yeasts have functional pheromone receptors and homeodomain homologues. This and the frequent existence of asexual strains within sexual species, makes the separation between sexual and asexual species imprecise. Events of loss of sexuality seem to be recent and frequent, but not uniformly distributed within the Sporidiobolales. Loss of sex could promote speciation by fostering the emergence of asexual lineages from an ancestral sexual stock, but does not seem to contribute to the generation of exclusively asexual lineages that persist for a long time.</p

HIV Evolution in Early Infection: Selection Pressures, Patterns of Insertion and Deletion, and the Impact of APOBEC

The pattern of viral diversification in newly infected individuals provides information about the host environment and immune responses typically experienced by the newly transmitted virus. For example, sites that tend to evolve rapidly across multiple early-infection patients could be involved in enabling escape from common early immune responses, could represent adaptation for rapid growth in a newly infected host, or could represent reversion from less fit forms of the virus that were selected for immune escape in previous hosts. Here we investigated the diversification of HIV-1 env coding sequences in 81 very early B subtype infections previously shown to have resulted from transmission or expansion of single viruses (n = 78) or two closely related viruses (n = 3). In these cases, the sequence of the infecting virus can be estimated accurately, enabling inference of both the direction of substitutions as well as distinction between insertion and deletion events. By integrating information across multiple acutely infected hosts, we find evidence of adaptive evolution of HIV-1 env and identify a subset of codon sites that diversified more rapidly than can be explained by a model of neutral evolution. Of 24 such rapidly diversifying sites, 14 were either i) clustered and embedded in CTL epitopes that were verified experimentally or predicted based on the individual's HLA or ii) in a nucleotide context indicative of APOBEC-mediated G-to-A substitutions, despite having excluded heavily hypermutated sequences prior to the analysis. In several cases, a rapidly evolving site was embedded both in an APOBEC motif and in a CTL epitope, suggesting that APOBEC may facilitate early immune escape. Ten rapidly diversifying sites could not be explained by CTL escape or APOBEC hypermutation, including the most frequently mutated site, in the fusion peptide of gp41. We also examined the distribution, extent, and sequence context of insertions and deletions, and we provide evidence that the length variation seen in hypervariable loop regions of the envelope glycoprotein is a consequence of selection and not of mutational hotspots. Our results provide a detailed view of the process of diversification of HIV-1 following transmission, highlighting the role of CTL escape and hypermutation in shaping viral evolution during the establishment of new infections

Advanced paternal age effects in neurodevelopmental disorders?review of potential underlying mechanisms

Multiple epidemiological studies suggest a relationship between advanced paternal age (APA) at conception and adverse neurodevelopmental outcomes in offspring, particularly with regard to increased risk for autism and schizophrenia. Conclusive evidence about how age-related changes in paternal gametes, or age-independent behavioral traits affect neural development is still lacking. Recent evidence suggests that the origins of APA effects are likely to be multidimensional, involving both inherited predisposition and de novo events. Here we provide a review of the epidemiological and molecular findings to date. Focusing on the latter, we present the evidence for genetic and epigenetic mechanisms underpinning the association between late fatherhood and disorder in offspring. We also discuss the limitations of the APA literature. We propose that different hypotheses relating to the origins of the APA effects are not mutually exclusive. Instead, multiple mechanisms likely contribute, reflecting the etiological complexity of neurodevelopmental disorders