Evolution of asexuality via different mechanisms in grass thrips (thysanoptera: Aptinothrips).

Asexual lineages can derive from sexual ancestors via different mechanisms and at variable rates, which affects the diversity of the asexual population and thereby its ecological success. We investigated the variation and evolution of reproductive systems in Aptinothrips, a genus of grass thrips comprising four species. Extensive population surveys and breeding experiments indicated sexual reproduction in A. elegans, asexuality in A. stylifer and A. karnyi, and both sexual and asexual lineages in A. rufus. Asexuality in A. stylifer and A. rufus coincides with a worldwide distribution, with sexual A. rufus lineages confined to a limited area. Inference of molecular phylogenies and antibiotic treatment revealed different causes of asexuality in different species. Asexuality in A. stylifer and A. karnyi has most likely genetic causes, while it is induced by endosymbionts in A. rufus. Endosymbiont-community characterization revealed presence of Wolbachia, and lack of other bacteria known to manipulate host reproduction. However, only 69% asexual A. rufus females are Wolbachia-infected, indicating that either an undescribed endosymbiont causes asexuality in this species or that Wolbachia was lost in several lineages that remained asexual. These results open new perspectives for studies on the maintenance of mixed sexual and asexual reproduction in natural populations.

Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

Effect of immune pressure on hepatitis C virus evolution: insights from a single-source outbreak.

The host's immune response to hepatitis C virus (HCV) can result in the selection of characteristic mutations (adaptations) that enable the virus to escape this response. The ability of the virus to mutate at these sites is dependent on the incoming virus, the fitness cost incurred by the mutation, and the benefit to the virus in escaping the response. Studies examining viral adaptation in chronic HCV infection have shown that these characteristic immune escape mutations can be observed at the population level as human leukocyte antigen (HLA)-specific viral polymorphisms. We examined 63 individuals with chronic HCV infection who were infected from a single HCV genotype 1b source. Our aim was to determine the extent to which the host's immune pressure affects HCV diversity and the ways in which the sequence of the incoming virus, including preexisting escape mutations, can influence subsequent mutations in recipients and infection outcomes. Conclusion: HCV sequences from these individuals revealed 29 significant associations between specific HLA types within the new hosts and variations within their viruses, which likely represent new viral adaptations. These associations did not overlap with previously reported adaptations for genotypes 1a and 3a and possibly reflected a combination of constraint due to the incoming virus and genetic distance between the strains. However, these sites accounted for only a portion of the sites in which viral diversity was observed in the new hosts. Furthermore, preexisting viral adaptations in the incoming (source) virus likely influenced the outcomes in the new hosts.

Perturbation expansions of multilocus fixation probabilities for frequency-dependent selection with applications to the Hill-Robertson effect and to the joint evolution of helping and punishment.

Natural populations are of finite size and organisms carry multilocus genotypes. There are, nevertheless, few results on multilocus models when both random genetic drift and natural selection affect the evolutionary dynamics. In this paper we describe a formalism to calculate systematic perturbation expansions of moments of allelic states around neutrality in populations of constant size. This allows us to evaluate multilocus fixation probabilities (long-term limits of the moments) under arbitrary strength of selection and gene action. We show that such fixation probabilities can be expressed in terms of selection coefficients weighted by mean first passages times of ancestral gene lineages within a single ancestor. These passage times extend the coalescence times that weight selection coefficients in one-locus perturbation formulas for fixation probabilities. We then apply these results to investigate the Hill-Robertson effect and the coevolution of helping and punishment. Finally, we discuss limitations and strengths of the perturbation approach. In particular, it provides accurate approximations for fixation probabilities for weak selection regimes only (Ns < or = 1), but it provides generally good prediction for the direction of selection under frequency-dependent selection.

Evolution and expression of vitellogenin genes.

A major event during the growth period of oocytes in nonmammalian animals is the accumulation of yolk. The genes coding for the yolk protein precursor, known as vitallogenin, are well characterized in a few vertebrate and invertebrate species. Studies on the evolution of these genes and on the regulatory mechanisms involved in their time, tissue- and hormone-specific expression are presented and discussed in this review.

Patterns of evolution of host proteins involved in retroviral pathogenesis.

BACKGROUND: Evolutionary analysis may serve as a useful approach to identify and characterize host defense and viral proteins involved in genetic conflicts. We analyzed patterns of coding sequence evolution of genes with known (TRIM5alpha and APOBEC3G) or suspected (TRIM19/PML) roles in virus restriction, or in viral pathogenesis (PPIA, encoding Cyclophilin A), in the same set of human and non-human primate species. RESULTS AND CONCLUSION: This analysis revealed previously unidentified clusters of positively selected sites in APOBEC3G and TRIM5alpha that may delineate new virus-interaction domains. In contrast, our evolutionary analyses suggest that PPIA is not under diversifying selection in primates, consistent with the interaction of Cyclophilin A being limited to the HIV-1M/SIVcpz lineage. The strong sequence conservation of the TRIM19/PML sequences among primates suggests that this gene does not play a role in antiretroviral defense.

Birth and adaptive evolution of a hominoid gene that supports high neurotransmitter flux.

The enzyme glutamate dehydrogenase (GDH) is important for recycling the chief excitatory neurotransmitter, glutamate, during neurotransmission. Human GDH exists in housekeeping and brain-specific isotypes encoded by the genes GLUD1 and GLUD2, respectively. Here we show that GLUD2 originated by retroposition from GLUD1 in the hominoid ancestor less than 23 million years ago. The amino acid changes responsible for the unique brain-specific properties of the enzyme derived from GLUD2 occurred during a period of positive selection after the duplication event.

Assessing rapid evolution in a changing environment.

Climate change poses a serious threat to species persistence. Effective modelling of evolutionary responses to rapid climate change is therefore essential. In this review we examine recent advances in phylogenetic comparative methods, techniques normally used to study adaptation over long periods, which allow them to be applied to the study of adaptation over shorter time scales. This increased applicability is largely due to the emergence of more flexible models of character evolution and the parallel development of molecular technologies that can be used to assess adaptive variation at loci scattered across the genome. The merging of phylogenetic and population genetic approaches to the study of adaptation has significant potential to advance our understanding of rapid responses to environmental change.

Molecular phylogenetics of shrews (Mammalia: Soricidae) reveal timing of transcontinental colonizations.

We sequenced 2167 base pairs (bp) of mitochondrial DNA cytochrome b and 16S, and 1390 bp of nuclear genes BRCA1 and ApoB in shrews taxa (Eulipotyphla, family Soricidae). The aim was to study the relationships at higher taxonomic levels within this family, and in particular the position of difficult clades such as Anourosorex and Myosorex. The data confirmed two monophyletic subfamilies, Soricinae and Crocidurinae. In the former, the tribes Anourosoricini, Blarinini, Nectogalini, Notiosoricini, and Soricini were supported. The latter was formed by the tribes Myosoricini and Crocidurini. The genus Suncus appeared to be paraphyletic and included Sylvisorex. We further suggest a biogeographical hypothesis, which shows that North America was colonized by three independent lineages of Soricinae during middle Miocene. Our hypothesis is congruent with the first fossil records for these taxa. Using molecular dating, the first exchanges between Africa and Eurasia occurred during the middle Miocene. The last one took place in the Late Miocene, with the dispersion of the genus Crocidura through the old world.

Splicing and the evolution of proteins in mammals.

It is often supposed that a protein's rate of evolution and its amino acid content are determined by the function and anatomy of the protein. Here we examine an alternative possibility, namely that the requirement to specify in the unprocessed RNA, in the vicinity of intron-exon boundaries, information necessary for removal of introns (e.g., exonic splice enhancers) affects both amino acid usage and rates of protein evolution. We find that the majority of amino acids show skewed usage near intron-exon boundaries, and that differences in the trends for the 2-fold and 4-fold blocks of both arginine and leucine show this to be owing to effects mediated at the nucleotide level. More specifically, there is a robust relationship between the extent to which an amino acid is preferred/avoided near boundaries and its enrichment/paucity in splice enhancers. As might then be expected, the rate of evolution is lowest near intron-exon boundaries, at least in part owing to splice enhancers, such that domains flanking intron-exon junctions evolve on average at under half the rate of exon centres from the same gene. In contrast, the rate of evolution of intronless retrogenes is highest near the domains where intron-exon junctions previously resided. The proportion of sequence near intron-exon boundaries is one of the stronger predictors of a protein's rate of evolution in mammals yet described. We conclude that after intron insertion selection favours modification of amino acid content near intron-exon junctions, so as to enable efficient intron removal, these changes then being subject to strong purifying selection even if nonoptimal for protein function. Thus there exists a strong force operating on protein evolution in mammals that is not explained directly in terms of the biology of the protein.

Gene organization of the mating type regions in the ectomycorrhizal fungus Laccaria bicolor reveals distinct evolution between the two mating type loci

In natural conditions, basidiomycete ectomycorrhizal fungi such as Laccaria bicolor are typically in the dikaryotic state when forming symbioses with trees, meaning that two genetically different individuals have to fuse or 'mate'. Nevertheless, nothing is known about the molecular mechanisms of mating in these ecologically important fungi. Here, advantage was taken of the first sequenced genome of the ectomycorrhizal fungus, Laccaria bicolor, to determine the genes that govern the establishment of cell-type identity and orchestrate mating. The L. bicolor mating type loci were identified through genomic screening. The evolutionary history of the genomic regions that contained them was determined by genome-wide comparison of L. bicolor sequences with those of known tetrapolar and bipolar basidiomycete species, and by phylogenetic reconstruction of gene family history. It is shown that the genes of the two mating type loci, A and B, are conserved across the Agaricales, but they are contained in regions of the genome with different evolutionary histories. The A locus is in a region where the gene order is under strong selection across the Agaricales. By contrast, the B locus is in a region where the gene order is likely under a low selection pressure but where gene duplication, translocation and transposon insertion are frequent.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Prenatal presentation and postnatal evolution of a patient with Jansen metaphyseal dysplasia with a novel missense mutation in PTH1R.

Wave-shaped ribs were detected at prenatal ultrasound in a 20(+1) week female fetus. At birth, skeletal radiographs showed marked hypomineralization and suggested hypophosphatasia. However, elevated blood calcium and alkaline phosphatase excluded hypophosphatasia and raised the possibility of Jansen metaphyseal dysplasia. Molecular analysis of the PTH/PTHrP receptor gene (PTH1R) showed heterozygosity for a previously undescribed transversion variant (c.1373T>A), which predicts p.Ile458Lys. In vitro evaluation of wild type and mutant PTH/PTHrP receptors supported the pathogenic role of the p.Ile458Lys substitution, and confirmed the diagnosis of Jansen metaphyseal dysplasia. This disorder may present prenatally with wavy ribs and in the newborn with hypomineralization, and may therefore be confused with hypophosphatasia. The mottled metaphyseal lesions typically associated with this disease appear only in childhood.

Evidence for the evolution of multiple genomes in arbuscular mycorrhizal fungi.

Ancient asexuals directly contradict the evolutionary theories that explain why organisms should evolve a sexual life history. The mutualistic, arbuscular mycorrhizal fungi are thought to have been asexual for approximately 400 million years. In the absence of sex, highly divergent descendants of formerly allelic nucleotide sequences are thought to evolve in a genome. In mycorrhizal fungi, where individual offspring receive hundreds of nuclei from the parent, it has been hypothesized that a population of genetically different nuclei should evolve within one individual. Here we use DNA-DNA fluorescent in situ hybridization to show that genetically different nuclei co-exist in individual arbuscular mycorrhizal fungi. We also show that the population genetics techniques used in other organisms are unsuitable for detecting recombination because the assumptions and underlying processes do not fit the fungal genomic structure shown here. Instead we used a phylogenetic approach to show that the within-individual genetic variation that occurs in arbuscular mycorrhizal fungi probably evolved through accumulation of mutations in an essentially clonal genome, with some infrequent recombination events. We conclude that mycorrhizal fungi have evolved to be multi-genomic.

Evolution at Two Levels in Fire Ants: The Relationship between Patterns of Gene Expression and Protein Sequence Evolution.

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.