ALBERT

Description: The genomes of related species contain valuable information on the history of the considered taxa. Great apes in particular exhibit variation of evolutionary patterns along their genomes. However, the great ape data also bring new challenges, such as the presence of incomplete lineage sorting and ancestral shared polymorphisms. Previous methods for genome-scale analysis are restricted to very few individuals or cannot disentangle the contribution of mutation rates and fixation biases. This represents a limitation both for the understanding of these forces as well as for the detection of regions affected by selection. Here, we present a new model designed to estimate mutation rates and fixation biases from genetic variation within and between species. We relax the assumption of instantaneous substitutions, modeling substitutions as mutational events followed by a gradual fixation. Hence, we straightforwardly account for shared ancestral polymorphisms and incomplete lineage sorting. We analyze genome-wide synonymous site alignments of human, chimpanzee, and two orangutan species. From each taxon, we include data from several individuals. We estimate mutation rates and GC-biased gene conversion intensity. We find that both mutation rates and biased gene conversion vary with GC content. We also find lineage-specific differences, with weaker fixation biases in orangutan species, suggesting a reduced historical effective population size. Finally, our results are consistent with directional selection acting on coding sequences in relation to exonic splicing enhancers.

Description: Genetic control of male or female gonad development displays between different groups of organisms a remarkable diversity of "master sex-determining genes" at the top of the genetic hierarchies, whereas downstream components surprisingly appear to be evolutionarily more conserved. Without much further studies, conservation of sequence has been equalized to conservation of function. We have used the medaka fish to investigate the generality of this paradigm. In medaka, the master male sex-determining gene is dmrt1bY , a highly conserved downstream regulator of sex determination in vertebrates. To understand its function in orchestrating the complex gene regulatory network, we have identified targets genes and regulated pathways of Dmrt1bY. Monitoring gene expression and interactions by transgenic fluorescent reporter fish lines, in vivo tissue-chromatin immunoprecipitation and in vitro gene regulation assays revealed concordance but also major discrepancies between mammals and medaka, notably amongst spatial, temporal expression patterns and regulations of the canonical Hedgehog and R-spondin/Wnt/Follistatin signaling pathways. Examination of Foxl2 protein distribution in the medaka ovary defined a new subpopulation of theca cells, where ovarian-type aromatase transcriptional regulation appears to be independent of Foxl2. In summary, these data show that the regulation of the downstream regulatory network of sex determination is less conserved than previously thought.

Description: Synonymous codon usage patterns are shaped by a balance between mutation, drift, and natural selection. To date, detection of translational selection in vertebrates has proven to be a challenging task, obscured by small long-term effective population sizes in larger animals and the existence of isochores in some species. The consensus is that, in such species, natural selection is either completely ineffective at overcoming mutational pressures and genetic drift or perhaps is effective but so weak that it is not detectable. The aim of this research is to understand the interplay between mutation, selection, and genetic drift in vertebrates. We observe that although variation in mutational bias is undoubtedly the dominant force influencing codon usage, translational selection acts as a weak additional factor influencing synonymous codon usage. These observations indicate that translational selection is a widespread phenomenon in vertebrates and is not limited to a few species.

Description: Class IV homeodomain leucine zipper (C4HDZ) genes are plant-specific transcription factors that, based on phenotypes in Arabidopsis thaliana , play an important role in epidermal development. In this study, we sampled all major extant lineages and their closest algal relatives for C4HDZ homologs and phylogenetic analyses result in a gene tree that mirrors land plant evolution with evidence for gene duplications in many lineages, but minimal evidence for gene losses. Our analysis suggests an ancestral C4HDZ gene originated in an algal ancestor of land plants and a single ancestral gene was present in the last common ancestor of land plants. Independent gene duplications are evident within several lineages including mosses, lycophytes, euphyllophytes, seed plants, and, most notably, angiosperms. In recently evolved angiosperm paralogs, we find evidence of pseudogenization via mutations in both coding and regulatory sequences. The increasing complexity of the C4HDZ gene family through the diversification of land plants correlates to increasing complexity in epidermal characters.

Description: Resolving difficult nodes for any part of the vertebrate tree of life often requires analyzing a large number of loci. Developing molecular markers that are workable for the groups of interest is often a bottleneck in phylogenetic research. Here, on the basis of a nested polymerase chain reaction (PCR) strategy, we present a universal toolkit including 102 nuclear protein-coding locus (NPCL) markers for vertebrate phylogenomics. The 102 NPCL markers have a broad range of evolutionary rates, which makes them useful for a wide range of time depths. The new NPCL toolkit has three important advantages compared with all previously developed NPCL sets: 1) the kit is universally applicable across vertebrates, with a PCR success rate of 94.6% in 16 widely divergent tested vertebrate species; 2) more than 90% of PCR reactions produce strong and single bands of the expected sizes that can be directly sequenced; and 3) all cleanup PCR reactions can be sequenced with only two specific universal primers. To test its actual phylogenetic utility, 30 NPCLs from this toolkit were used to address the higher level relationships of living salamanders. Of the 639 target PCR reactions performed on 19 salamanders and several outgroup species, 632 (98.9%) were successful, and 602 (94.1%) were directly sequenced. Concatenation and species-tree analyses on this 30-locus data set produced a fully resolved phylogeny and showed that Cryptobranchoidea (Cryptobranchidae + Hynobiidae) branches first within the salamander tree, followed by Sirenidae. Our experimental tests and our demonstration for a particular case show that our NPCL toolkit is a highly reliable, fast, and cost-effective approach for vertebrate phylogenomic studies and thus has the potential to accelerate the completion of many parts of the vertebrate tree of life.

Description: Genetic incompatibilities are commonly observed between hybridizing species. Although this type of isolating mechanism has received considerable attention, we have few examples describing how genetic incompatibilities evolve. We investigated the evolution of two loci involved in a classic example of a Bateson–Dobzhansky–Muller (BDM) incompatibility in Xiphophorus , a genus of freshwater fishes from northern Central America. Hybrids develop a lethal melanoma due to the interaction of two loci, an oncogene and its repressor. We cloned and sequenced the putative repressor locus in 25 Xiphophorus species and an outgroup species, and determined the status of the oncogene in those species from the literature. Using phylogenetic analyses, we find evidence that a repeat region in the proximal promoter of the repressor is coevolving with the oncogene. The data support a hypothesis that departs from the standard BDM model: it appears the alleles that cause the incompatibilities have coevolved simultaneously within lineages, rather than in allopatric or temporal isolation.

Description: Incompatibility systems in which individuals bearing identical alleles reject each other favor the maintenance of a diversity of alleles. Mushroom mating type loci ( MAT ) encode for dozens or hundreds of incompatibility alleles whose loss from the population is greatly restricted through negative frequency selection, leading to a system of alleles with highly divergent sequences. Here, we use DNA sequences of homeodomain (HD) encoding genes at the MAT locus of five closely related species of the root rot basidiomycete Heterobasidion annosum sensu lato to show that the extended coalescence time of MAT alleles greatly predates speciation in the group, contrasting loci outside of MAT that show allele divergences largely consistent with the species phylogeny with those of MAT , which show rampant trans-species polymorphism. We observe a roughly 6-fold greater genealogical depth and polymorphism of MAT compared with non- MAT that argues for the maintenance of balanced polymorphism for a minimum duration of 24 My based on a molecular-clock calibrated species phylogeny. As with other basidiomycete HD genes, balancing selection appears to be concentrated at the specificity-determining region in the N-terminus of the protein based on identification of codons under selection and the absence of recombination within the region. However, the elevated polymorphism extends into the nonspecificity determining regions as well as a neighboring non- MAT gene, the mitochondrial intermediate peptidase ( MIP ). In doing so, increased divergence should decrease recombination among alleles and as a by-product create incompatibilities in the functional domains not involved in allele recognition but in regulating sexual development.

Description: Here we present computational machinery to efficiently and accurately identify transposable element (TE) insertions in 146 next-generation sequenced inbred strains of Drosophila melanogaster . The panel of lines we use in our study is composed of strains from a pair of genetic mapping resources: the Drosophila Genetic Reference Panel (DGRP) and the Drosophila Synthetic Population Resource (DSPR). We identified 23,087 TE insertions in these lines, of which 83.3% are found in only one line. There are marked differences in the distribution of elements over the genome, with TEs found at higher densities on the X chromosome, and in regions of low recombination. We also identified many more TEs per base pair of intronic sequence and fewer TEs per base pair of exonic sequence than expected if TEs are located at random locations in the euchromatic genome. There was substantial variation in TE load across genes. For example, the paralogs derailed and derailed-2 show a significant difference in the number of TE insertions, potentially reflecting differences in the selection acting on these loci. When considering TE families, we find a very weak effect of gene family size on TE insertions per gene, indicating that as gene family size increases the number of TE insertions in a given gene within that family also increases. TEs are known to be associated with certain phenotypes, and our data will allow investigators using the DGRP and DSPR to assess the functional role of TE insertions in complex trait variation more generally. Notably, because most TEs are very rare and often private to a single line, causative TEs resulting in phenotypic differences among individuals may typically fail to replicate across mapping panels since individual elements are unlikely to segregate in both panels. Our data suggest that "burden tests" that test for the effect of TEs as a class may be more fruitful.

Description: Sex determination mechanisms are highly variable across teleost fishes and sexual development is often plastic. Nevertheless, downstream factors establishing the two sexes are presumably conserved. Here, we study sequence evolution and gene expression of core genes of sexual development in a prime model system in evolutionary biology, the East African cichlid fishes. Using the available five cichlid genomes, we test for signs of positive selection in 28 genes including duplicates from the teleost whole-genome duplication, and examine the expression of these candidate genes in three cichlid species. We then focus on a particularly striking case, the A- and B-copies of the aromatase cyp19a1 , and detect different evolutionary trajectories: cyp19a1A evolved under strong positive selection, whereas cyp19a1B remained conserved at the protein level, yet is subject to regulatory changes at its transcription start sites. Importantly, we find shifts in gene expression in both copies. Cyp19a1 is considered the most conserved ovary-factor in vertebrates, and in all teleosts investigated so far, cyp19a1A and cyp19a1B are expressed in ovaries and the brain, respectively. This is not the case in cichlids, where we find new expression patterns in two derived lineages: the A-copy gained a novel testis-function in the Ectodine lineage, whereas the B-copy is overexpressed in the testis of the speciest-richest cichlid group, the Haplochromini. This suggests that even key factors of sexual development, including the sex steroid pathway, are not conserved in fish, supporting the idea that flexibility in sexual determination and differentiation may be a driving force of speciation.

Description: 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.

Description: The ability to survey polymorphism on a genomic scale has enabled genome-wide scans for the targets of natural selection. Theory that connects patterns of genetic variation to evidence of natural selection most often assumes a diallelic locus and no recurrent mutation. Although these assumptions are suitable to selection that targets single nucleotide variants, fundamentally different types of mutation generate abundant polymorphism in genomes. Moreover, recent empirical results suggest that mutationally complex, multiallelic loci including microsatellites and copy number variants are sometimes targeted by natural selection. Given their abundance, the lack of inference methods tailored to the mutational peculiarities of these types of loci represents a notable gap in our ability to interrogate genomes for signatures of natural selection. Previous theoretical investigations of mutation-selection balance at multiallelic loci include assumptions that limit their application to inference from empirical data. Focusing on microsatellites, we assess the dynamics and population-level consequences of selection targeting mutationally complex variants. We develop general models of a multiallelic fitness surface, a realistic model of microsatellite mutation, and an efficient simulation algorithm. Using these tools, we explore mutation-selection-drift equilibrium at microsatellites and investigate the mutational history and selective regime of the microsatellite that causes Friedreich’s ataxia. We characterize microsatellite selective events by their duration and cost, note similarities to sweeps from standing point variation, and conclude that it is premature to label microsatellites as ubiquitous agents of efficient adaptive change. Together, our models and simulation algorithm provide a powerful framework for statistical inference, which can be used to test the neutrality of microsatellites and other multiallelic variants.

Description: Although endogenous retroviruses are common across vertebrate genomes, the koala retrovirus (KoRV) is the only retrovirus known to be currently invading the germ line of its host. KoRV is believed to have first infected koalas in northern Australia less than two centuries ago. We examined KoRV in 28 koala museum skins collected in the late 19th and 20th centuries and deep sequenced the complete proviral envelope region from five northern Australian specimens. Strikingly, KoRV env sequences were conserved among koalas collected over the span of a century, and two functional motifs that affect viral infectivity were fixed across the museum koala specimens. We detected only 20 env polymorphisms among the koalas, likely representing derived mutations subject to purifying selection. Among northern Australian koalas, KoRV was already ubiquitous by the late 19th century, suggesting that KoRV evolved and spread among koala populations more slowly than previously believed. Given that museum and modern koalas share nearly identical KoRV sequences, it is likely that koala populations, for more than a century, have experienced increased susceptibility to diseases caused by viral pathogenesis.

Description: Protein interaction networks play central roles in biological systems, from simple metabolic pathways through complex programs permitting the development of organisms. Multicellularity could only have arisen from a careful orchestration of cellular and molecular roles and responsibilities, all properly controlled and regulated. Disease reflects a breakdown of this organismal homeostasis. To better understand the evolution of interactions whose dysfunction may be contributing factors to disease, we derived the human protein coevolution network using our MatrixMatchMaker algorithm and using the Orthologous MAtrix project (OMA) database as a source for protein orthologs from 103 eukaryotic genomes. We annotated the coevolution network using protein–protein interaction data, many functional data sources, and we explored the evolutionary rates and dates of emergence of the proteins in our data set. Strikingly, clustering based only on the topology of the coevolution network partitions it into two subnetworks, one generally representing ancient eukaryotic functions and the other functions more recently acquired during animal evolution. That latter subnetwork is enriched for proteins with roles in cell–cell communication, the control of cell division, and related multicellular functions. Further annotation using data from genetic disease databases and cancer genome sequences strongly implicates these proteins in both ciliopathies and cancer. The enrichment for such disease markers in the animal network suggests a functional link between these coevolving proteins. Genetic validation corroborates the recruitment of ancient cilia in the evolution of multicellularity.

Description: Here, we sequenced the 5,419,609 bp circular genome of an Enterobacter aerogenes clinical isolate that killed a patient and was resistant to almost all current antibiotics (except gentamicin) commonly used to treat Enterobacterial infections, including colistin. Genomic and phylogenetic analyses explain the discrepancies of this bacterium and show that its core genome originates from another genus, Klebsiella. Atypical characteristics of this bacterium (i.e., motility, presence of ornithine decarboxylase, and lack of urease activity) are attributed to genomic mosaicism, by acquisition of additional genes, such as the complete 60,582 bp flagellar assembly operon acquired "en bloc" from the genus Serratia. The genealogic tree of the 162,202 bp multidrug-resistant conjugative plasmid shows that it is a chimera of transposons and integrative conjugative elements from various bacterial origins, resembling a rhizome. Moreover, we demonstrate biologically that a G53S mutation in the pmrA gene results in colistin resistance. E. aerogenes has a large RNA population comprising 8 rRNA operons and 87 cognate tRNAs that have the ability to translate transferred genes that use different codons, as exemplified by the significantly different codon usage between genes from the core genome and the "mobilome." On the basis of our findings, the evolution of this bacterium to become a "killer bug" with new genomic repertoires was from three criteria that are "opportunity, power, and usage" to indicate a sympatric lifestyle: "opportunity" to meet other bacteria and exchange foreign sequences since this bacteria was similar to sympatric bacteria; "power" to integrate these foreign sequences such as the acquisition of several mobile genetic elements (plasmids, integrative conjugative element, prophages, transposons, flagellar assembly system, etc.) found in his genome; and "usage" to have the ability to translate these sequences including those from rare codons to serve as a translator of foreign languages.

Description: Diatoms are the most species-rich group of microalgae, and their contribution to marine primary production is important on a global scale. Diatoms can form dense blooms through rapid asexual reproduction; mutations acquired and propagated during blooms likely provide the genetic, and thus phenotypic, variability upon which natural selection may act. Positive selection was tested using genome and transcriptome-wide pair-wise comparisons of homologs in three genera of diatoms ( Pseudo-nitzschia, Ditylum, and Thalassiosira ) that represent decreasing phylogenetic distances. The signal of positive selection was greatest between two strains of Thalassiosira pseudonana . Further testing among seven strains of T. pseudonana yielded 809 candidate genes of positive selection, which are 7% of the protein-coding genes. Orphan genes and genes encoding protein-binding domains and transcriptional regulators were enriched within the set of positively selected genes relative to the genome as a whole. Positively selected genes were linked to the potential selective pressures of nutrient limitation and sea surface temperature based on analysis of gene expression profiles and identification of positively selected genes in subsets of strains from locations with similar environmental conditions. The identification of positively selected genes presents an opportunity to test new hypotheses in natural populations and the laboratory that integrate selected genotypes in T. pseudonana with their associated phenotypes and selective forces.

Description: Antisense transcription, or transcription on the opposite strand of the same genomic locus as another transcript, has been observed in many organisms, including yeast. Several antisense transcripts are known to be conserved across various species of yeast, and a few antisense transcripts are associated with functional regulation of the sense transcript. We detect antisense transcription from approximately 90% of protein-coding genes, and antisense transcription is generally associated with histone modifications indicative of a transcriptionally active state. The pattern of genome-wide antisense transcription in two species of budding yeast, Saccharomyces cerevisiae and S. paradoxus , is widely evolutionarily conserved: Antisense transcripts exhibit conserved expression levels and localization with respect to gene annotations. Regions of genes exhibiting conserved antisense transcription also show less sequence divergence than regions of genes without antisense transcription. These findings provide further support that widespread antisense transcription is functional in yeast, and expand the catalog of putative functional antisense transcripts to include nonpolyadenylated transcripts. Because antisense transcripts are less divergent in sequence than expected, they likely contain sequences important to their function.

Description: Reconstruction of the past is an important task of evolutionary biology. It takes place at different points in a hierarchy of molecular variation, including genes, individuals, populations, and species. Statistical inference about population histories has recently received considerable attention, following the development of computational tools to provide tractable approaches to this very challenging problem. Here, we introduce a likelihood-based approach which generalizes a recently developed model for random fluctuations in allele frequencies based on an approximation to the neutral Wright–Fisher diffusion. Our new framework approximates the infinite alleles Wright–Fisher model and uses an implementation with an adaptive Markov chain Monte Carlo algorithm. The method is especially well suited to data sets harboring large population samples and relatively few loci for which other likelihood-based models are currently computationally intractable. Using our model, we reconstruct the global population history of a major human pathogen, Streptococcus pneumoniae . The results illustrate the potential to reach important biological insights to an evolutionary process by a population genetics approach, which can appropriately accommodate very large population samples.

Description: Humans’ ability for rapid dispersal and adaptation has allowed us to colonize diverse geographic and climatic regions of the planet, creating a complex evolutionary history. This complexity can be understood, at least partially, by modeling the underlying demographic parameters in the evolutionary process. In this study, we analyze a model of human evolution in which population size, gene flow (GF), and time are varied. Specifically, we simulate mitochondrial DNA for 42 demographic scenarios, represented by 42 parameter combinations, to describe the initial dispersal of modern humans out of Africa. The analyses include three values for colonization size (CS; 1%, 10%, and 30% of the African population), seven values for rate of GF (10 –6 –0.5), and two values for time of colonization (50,000 and 100,000 years ago). We then estimate summary statistics for the simulated data sets to calculate the percent of explained variation by each parameter and to identify which parameter combinations generate distinct differences in genetic variation, that is, which demographic scenarios can be distinguished from each other. On the basis of these results, we make recommendations about which summary statistics to use according to the parameter of interest. Our results show that CS, GF, and their interaction have the largest effect on genetic variation under our model of human evolution. Comparison with empirical data suggests that 1% of the existing African mitochondrial genetic variation left and colonized the rest of the world (i.e., CS = 1%) and bidirectional GF continued at a level of ~10 individuals per generation (i.e., GF = 10 –3 ) after the initial colonization. Our study serves as a model to bridge the gap between the use of simulations for theoretical population genetics and empirical data analysis such as approximate Bayesian computation approaches and is, thus, applicable to the study of molecular evolution in any organism.

Description: In recent years, ancient DNA has increasingly been used for estimating molecular timescales, particularly in studies of substitution rates and demographic histories. Molecular clocks can be calibrated using temporal information from ancient DNA sequences. This information comes from the ages of the ancient samples, which can be estimated by radiocarbon dating the source material or by dating the layers in which the material was deposited. Both methods involve sources of uncertainty. The performance of Bayesian phylogenetic inference depends on the information content of the data set, which includes variation in the DNA sequences and the structure of the sample ages. Various sources of estimation error can reduce our ability to estimate rates and timescales accurately and precisely. We investigated the impact of sample-dating uncertainties on the estimation of evolutionary timescale parameters using the software BEAST. Our analyses involved 11 published data sets and focused on estimates of substitution rate and root age. We show that, provided that samples have been accurately dated and have a broad temporal span, it might be unnecessary to account for sample-dating uncertainty in Bayesian phylogenetic analyses of ancient DNA. We also investigated the sample size and temporal span of the ancient DNA sequences needed to estimate phylogenetic timescales reliably. Our results show that the range of sample ages plays a crucial role in determining the quality of the results but that accurate and precise phylogenetic estimates of timescales can be made even with only a few ancient sequences. These findings have important practical consequences for studies of molecular rates, timescales, and population dynamics.

Description: Genetic exchange by conjugation is responsible for the spread of resistance, virulence, and social traits among prokaryotes. Recent works unraveled the functioning of the underlying type IV secretion systems (T4SS) and its distribution and recruitment for other biological processes (exaptation), notably pathogenesis. We analyzed the phylogeny of key conjugation proteins to infer the evolutionary history of conjugation and T4SS. We show that single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) conjugation, while both based on a key AAA + ATPase, diverged before the last common ancestor of bacteria. The two key ATPases of ssDNA conjugation are monophyletic, having diverged at an early stage from dsDNA translocases. Our data suggest that ssDNA conjugation arose first in diderm bacteria, possibly Proteobacteria, and then spread to other bacterial phyla, including bacterial monoderms and Archaea. Identifiable T4SS fall within the eight monophyletic groups, determined by both taxonomy and structure of the cell envelope. Transfer to monoderms might have occurred only once, but followed diverse adaptive paths. Remarkably, some Firmicutes developed a new conjugation system based on an atypical relaxase and an ATPase derived from a dsDNA translocase. The observed evolutionary rates and patterns of presence/absence of specific T4SS proteins show that conjugation systems are often and independently exapted for other functions. This work brings a natural basis for the classification of all kinds of conjugative systems, thus tackling a problem that is growing as fast as genomic databases. Our analysis provides the first global picture of the evolution of conjugation and shows how a self-transferrable complex multiprotein system has adapted to different taxa and often been recruited by the host. As conjugation systems became specific to certain clades and cell envelopes, they may have biased the rate and direction of gene transfer by conjugation within prokaryotes.

Description: Hydrothermal vents are typically located in midocean ridges and back-arc basins and are usually generated by the movement of tectonic plates. Life thrives in these environments despite the extreme conditions. In addition to chemoautotrophic bacteria, decapod crustaceans are dominant in many of the hydrothermal vents discovered to date. Contrary to the hypothesis that these species are remnants of relic fauna, increasing evidence supports the notion that hydrothermal vent decapods have diversified in more recent times with previous research attributing the origin of alvinocarid shrimps to the Miocene. This study investigated seven representative decapod species from four hydrothermal vents throughout the Western Pacific and Indian Oceans. A partitioned mix-model phylogenomic analysis of mitochondrial DNA produced a consistent phylogenetic topology of these vent-endemic species. Additionally, molecular dating analysis calibrated using multiple fossils suggested that both bythograeid crabs and alvinocarid shrimps originated in the late Mesozoic and early Cenozoic. Although of limited sampling, our estimates support the extinction/repopulation hypothesis, which postulates recent diversification times for most hydrothermal vent species due to their mass extinction by global deep-water anoxic/dysoxic events during the Late Cretaceous and Early Tertiary. The continental-derived property of the West Pacific province is compatible with the possibility that vent decapods diversified from ancestors from shallow-water regions such as cold seeps. Our results move us a step closer toward understanding the evolutionary origin of hydrothermal vent species and their distribution in the Western Pacific–Indian Ocean Region.

Description: Proteins in the superfamily of voltage-gated ion channels mediate behavior across the tree of life. These proteins regulate the movement of ions across cell membranes by opening and closing a central pore that controls ion flow. The best-known members of this superfamily are the voltage-gated potassium, calcium (Ca v ), and sodium (Na v ) channels, which underlie impulse conduction in nerve and muscle. Not all members of this family are opened by changes in voltage, however. NALCN (NA + leak channel nonselective) channels, which encode a voltage-insensitive "sodium leak" channel, have garnered a growing interest. This study examines the phylogenetic relationship among Na v /Ca v voltage-gated and voltage-insensitive channels in the eukaryotic group Opisthokonta, which includes animals, fungi, and their unicellular relatives. We show that NALCN channels diverged from voltage-gated channels before the divergence of fungi and animals and that the closest relatives of NALCN channels are fungal calcium channels, which they functionally resemble.

Description: Newer parts of sex chromosomes, neo-sex chromosomes, offer unique possibilities for studying gene degeneration and sequence evolution in response to loss of recombination and population size decrease. We have recently described a neo-sex chromosome system in Sylvioidea passerines that has resulted from a fusion between the first half (10 Mb) of chromosome 4a and the ancestral sex chromosomes. In this study, we report the results of molecular analyses of neo-Z and neo-W gametologs and intronic parts of neo-Z and autosomal genes on the second half of chromosome 4a in three species within different Sylvioidea lineages (Acrocephalidea, Timaliidae, and Alaudidae). In line with hypotheses of neo-sex chromosome evolution, we observe 1) lower genetic diversity of neo-Z genes compared with autosomal genes, 2) moderate synonymous and weak nonsynonymous sequence divergence between neo-Z and neo-W gametologs, and 3) lower GC content on neo-W than neo-Z gametologs. Phylogenetic reconstruction of eight neo-Z and neo-W gametologs suggests that recombination continued after the split of Alaudidae from the rest of the Sylvioidea lineages (i.e., after ~42.2 Ma) and with some exceptions also after the split of Acrocephalidea and Timaliidae (i.e., after ~39.4 Ma). The Sylvioidea neo-sex chromosome shares classical evolutionary features with the ancestral sex chromosomes but, as expected from its more recent origin, shows weaker divergence between gametologs.

Description: In bacteria, physiological change may be effected by a single gene acquisition, producing ecological differentiation without genetic isolation. Natural selection acting on such differences can reduce the frequency of genotypes that arise from recombination at these loci. However, gene acquisition can only account for recombination interference in the fraction of the genome that is tightly linked to the integration site. To identify additional loci that contribute to adaptive differences, we examined orthologous genes in species of Enterobacteriaceae to identify significant differences in the degree of codon selection. Significance was assessed using the Adaptive Codon Enrichment metric, which accounts for the variation in codon usage bias that is expected to arise from mutation and drift; large differences in codon usage bias were identified in more genes than would be expected to arise from stochastic processes alone. Genes in the same operon showed parallel differences in codon usage bias, suggesting that changes in the overall levels of gene expression led to changes in the degree of adaptive codon usage. Most significant differences between orthologous operons were found among those involved with specific environmental adaptations, whereas "housekeeping" genes rarely showed significant changes. When considered together, the loci experiencing significant changes in codon selection outnumber potentially adaptive gene acquisition events. The identity of genes under strong codon selection seems to be influenced by the habitat from which the bacteria were isolated. We propose a two-stage model for how adaptation to different selective regimes can drive bacterial speciation. Initially, gene acquisitions catalyze rapid ecological differentiation, which modifies the utilization of genes, thereby changing the strength of codon selection on them. Alleles develop fitness variation by substitution, producing recombination interference at these loci in addition to those flanking acquired genes, allowing sequences to diverge across the entire genome and establishing genetic isolation (i.e., protection from frequent homologous recombination).

Description: Gene expression levels correlate with multiple aspects of gene sequence and gene structure in phylogenetically diverse taxa, suggesting an important role of gene expression levels in the evolution of protein-coding genes. Here we present results of a genome-wide study of the influence of gene expression on synonymous codon usage, amino acid composition, and gene structure in the red flour beetle, Tribolium castaneum . Consistent with the action of translational selection, we find that synonymous codon usage bias increases with gene expression. However, the correspondence between tRNA gene copy number and optimal codons is weak. At the amino acid level, translational selection is suggested by the positive correlation between tRNA gene numbers and amino acid usage, which is stronger for highly expressed genes. In addition, there is a clear trend for increased use of metabolically cheaper, less complex amino acids as gene expression increases. tRNA gene numbers also correlate negatively with amino acid size/complexity (S/C) score indicating the coupling between translational selection and selection to minimize the use of large/complex amino acids. Interestingly, the analysis of 10 additional genomes suggests that the correlation between tRNA gene numbers and amino acid S/C score is widespread and might be explained by selection against negative consequences of protein misfolding. At the level of gene structure, three major trends are detected: 1) complete coding region length increases across low and intermediate expression levels but decreases in highly expressed genes; 2) the average intron size shows the opposite trend, first decreasing with expression, followed by a slight increase in highly expressed genes; and 3) intron density remains nearly constant across all expression levels. These changes in gene architecture are only in partial agreement with selection favoring reduced cost of biosynthesis.

Description: New protein-coding genes can originate either through modification of existing genes or de novo. Recently, the importance of de novo origination has been recognized in eukaryotes, although eukaryotic genes originated de novo are relatively rare and difficult to identify. In contrast, viruses contain many de novo genes , namely those in which an existing gene has been "overprinted" by a new open reading frame, a process that generates a new protein-coding gene overlapping the ancestral gene. We analyzed the evolution of 12 experimentally validated viral genes that originated de novo and estimated their relative ages. We found that young de novo genes have a different codon usage from the rest of the genome. They evolve rapidly and are under positive or weak purifying selection. Thus, young de novo genes might have strain-specific functions, or no function, and would be difficult to detect using current genome annotation methods that rely on the sequence signature of purifying selection. In contrast to young de novo genes, older de novo genes have a codon usage that is similar to the rest of the genome. They evolve slowly and are under stronger purifying selection. Some of the oldest de novo genes evolve under stronger selection pressure than the ancestral gene they overlap, suggesting an evolutionary tug of war between the ancestral and the de novo gene.

Description: Here, we present a study of the molecular evolution of the pheromone receptor genes ( pre-1 and pre-2 ) in Neurospora taxa with different mating systems. We focus on comparisons between heterothallic and homothallic taxa, reproducing sexually by outcrossing and by intrahaploid selfing, respectively. Our general aim was to use a phylogenetic framework to investigate whether the evolutionary trajectory of the pheromone and receptor genes in Neurospora differs between heterothallic and homothallic taxa, and among the homothallic lineages/clades previously indicated to represent independent switches from heterothallism to homothallism in the evolutionary history of the genus. We complemented molecular evolution analyses with an expression study of the pre genes and their upstream regulators, the mating-type ( mat ) genes, in homothallic taxa. Our analyses suggest that the pheromone receptor gene pre-1 is functionally conserved in both heterothallic and homothallic taxa. Moreover, we found evidence of positive selection for a small fraction of codons in the cytoplasmic signal-transducing C-terminal region of the protein PRE-1. Distribution of positively selected codons differs between heterothallic and homothallic groups, suggesting functional divergence associated with mating system. The gene pre-2 was shown to evolve under high selective constraints, with no strong evidence for positive selection. Although our data suggest that both pre-1 and pre-2 are overall functional in homothallic taxa, individual taxa display frame-shift mutations causing premature stop codons, which might indicate loss of function. Transcriptional patterns of pre and mat genes in six homothallic taxa, selected to represent six different switches from heterothallism to homothallism, do not support a universal pattern of regulation of these genes during reproductive tissue development. Taken together, our analyses suggest that the pheromone receptor genes pre-1 and pre-2 are in general functional in homothallic Neurospora taxa, in contrast with the situation for the mat genes that are generally degenerate in these taxa.

Description: Comparative genome biology has unveiled the polyploid origin of all angiosperms and the role of recurrent polyploidization in the amplification of gene families and the structuring of genomes. Which species share certain ancient polyploidy events, and which do not, is ill defined because of the limited number of sequenced genomes and transcriptomes and their uneven phylogenetic distribution. Previously, it has been suggested that most, but probably not all, of the eudicots have shared an ancient hexaploidy event, referred to as the gamma triplication. In this study, detailed phylogenies of subfamilies of MADS-box genes suggest that the gamma triplication has occurred before the divergence of Gunnerales but after the divergence of Buxales and Trochodendrales. Large-scale phylogenetic and K S - based approaches on the inflorescence transcriptomes of Gunnera manicata (Gunnerales) and Pachysandra terminalis (Buxales) provide further support for this placement, enabling us to position the gamma triplication in the stem lineage of the core eudicots. This triplication likely initiated the functional diversification of key regulators of reproductive development in the core eudicots, comprising 75% of flowering plants. Although it is possible that the gamma event triggered early core eudicot diversification, our dating estimates suggest that the event occurred early in the stem lineage, well before the rapid speciation of the earliest core eudicot lineages. The evolutionary significance of this paleopolyploidy event may thus rather lie in establishing a species lineage that was resilient to extinction, but with the genomic potential for later diversification. We consider that the traits generated from this potential characterize extant core eudicots both chemically and morphologically.

Description: Genome reduction in obligately intracellular bacteria is one of the most well-established patterns in the field of molecular evolution. In the extreme, many sap-feeding insects harbor nutritional symbionts with genomes that are so reduced that it is not clear how they perform basic cellular functions. For example, the primary symbiont of psyllids ( Carsonella ) maintains one of the smallest and most AT-rich bacterial genomes ever identified and has surprisingly lost many genes that are thought to be essential for its role in provisioning its host with amino acids. However, our understanding of this extreme case of genome reduction is limited, as genomic data for Carsonella are available from only a single host species, and little is known about the functional role of "secondary" bacterial symbionts in psyllids. To address these limitations, we analyzed complete Carsonella genomes from pairs of congeneric hosts in three divergent genera within the Psyllidae ( Ctenarytaina , Heteropsylla , and Pachypsylla ) as well as complete secondary symbiont genomes from two of these host species ( Ctenarytaina eucalypti and Heteropsylla cubana ). Although the Carsonella genomes are generally conserved in size, structure, and GC content and exhibit genome-wide signatures of purifying selection, we found that gene loss has remained active since the divergence of the host species and had a particularly large impact on the amino acid biosynthesis pathways that define the symbiotic role of Carsonella . In some cases, the presence of additional bacterial symbionts may compensate for gene loss in Carsonella , as functional gene content indicates a high degree of metabolic complementarity between co-occurring symbionts. The genomes of the secondary symbionts also show signatures of long-term evolution as vertically transmitted, intracellular bacteria, including more extensive genome reduction than typically observed in facultative symbionts. Therefore, a history of co-evolution with secondary bacterial symbionts can partially explain the ongoing genome reduction in Carsonella . However, the absence of these secondary symbionts in other host lineages indicates that the relationships are dynamic and that other mechanisms, such as changes in host diet or functional coordination with the host genome, must also be at play.

Description: Gene duplication is a major driver of organismal adaptation and evolution and plays an important role in multiple human diseases. Whole-genome analyses have shown similar and high rates of gene duplication across a variety of eukaryotic species. Most of these studies, however, did not address the possible impact of interlocus gene conversion (IGC) on the evolution of gene duplicates. Because IGC homogenizes pairs of duplicates, widespread conversion would cause gene duplication events that happened long ago to appear more recent, resulting in artificially high estimates of duplication rates. Although the majority of genome-wide studies (including in the budding yeast Saccharomyces cerevisiae [Scer]) point to levels of IGC between paralogs ranging from 2% to 18%, Gao and Innan (Gao LZ, Innan H. 2004. Very low gene duplication rate in the yeast genome. Science 306:1367–1370.) found that gene conversion in yeast affected 〉80% of paralog pairs. If conversion rates really are this high, it would imply that the rate of gene duplication in eukaryotes is much lower than previously reported. In this work, we apply four different methodologies—including one approach that closely mirrors Gao and Innan’s method—to estimate the level of IGC in Scer. Our analyses point to a maximum conversion level of 13% between paralogs in this species, in close agreement with most estimates of IGC in eukaryotes. We also show that the exceedingly high levels of conversion found previously derive from application of an accurate method to an inappropriate data set. In conclusion, our work provides the most striking evidence to date supporting the reduced incidence of IGC among Scer paralogs and sets up a framework for future analyses in other eukaryotes.

Description: Cyclin B3 evolution has the unique peculiarity of an abrupt 3-fold increase of the protein size in the mammalian lineage due to the extension of a single exon. We have analyzed the evolution of the gene to define the modalities of this event and the possible consequences on the function of the protein. Database searches can trace the appearance of the gene to the origin of metazoans. Most introns were already present in early metazoans, and the intron–exon structure as well as the protein size were fairly conserved in invertebrates and nonmammalian vertebrates. Although intron gains are considered as rare events, we identified two cases, one at the prochordate–chordate transition and one in murids, resulting from different mechanisms. At the emergence of mammals, the gene was relocated from chromosome 6 of platypus to the X chromosome in marsupials, but the exon extension occurred only in placental mammals. A repetitive structure of 18 amino acids, of uncertain origin, is detectable in the 3,000-nt mammalian exon-encoded sequence, suggesting an extension by multiple internal duplications, some of which are still detectable in the primate lineage. Structure prediction programs suggest that the repetitive structure has no associated three-dimensional structure but rather a tendency for disorder. Splice variant isoforms were detected in several mammalian species but without conserved pattern, notably excluding the constant coexistence of premammalian-like transcripts, without the extension. The yeast two-hybrid method revealed that, in human, the extension allowed new interactions with ten unrelated proteins, most of them with specific three-dimensional structures involved in protein–protein interactions, and some highly expressed in testis, as is cyclin B3. The interactions with activator of cAMP-responsive element modulator in testis (ACT), germ cell-less homolog 1, and chromosome 1 open reading frame 14 remain to be verified in vivo since they may not be expressed at the same stages of spermatogenesis as cyclin B3.

Description: The X chromosome has a large effect on hybrid dysfunction, particularly on hybrid male sterility. Although the evidence for this so-called large-X effect is clear, its molecular causes are not yet fully understood. One possibility is that, under certain conditions, evolution proceeds faster in X-linked than in autosomal loci (i.e., faster-X effect) due to both natural selection and their hemizygosity in males, an effect that is expected to be greatest in genes with male-biased expression. Here, I study genome-wide variation in transcript abundance between Drosophila yakuba and D. santomea , within these species and in their hybrid males to evaluate both the faster-X and large-X effects at the level of expression. I find that in X-linked male-biased genes (MBGs) expression evolves faster than in their autosomal counterparts, an effect that is accompanied by a unique reduction in expression polymorphism. This suggests that Darwinian selection is driving expression differences between species, likely enhanced by the hemizygosity of the X chromosome in males. Despite the recent split of the two sister species under study, abundant changes in both cis - and trans -regulatory elements underlie expression divergence in the majority of the genes analyzed, with significant differences in allelic ratios of transcript abundance between the two reciprocal F 1 hybrid males. Cis–trans coevolution at molecular level, evolved shortly after populations become isolated, may therefore contribute to explain the breakdown of the regulation of gene expression in hybrid males. Additionally, the X chromosome plays a large role in this hybrid male misexpression, which affects not only MBG but also, to a lesser degree, nonsex-biased genes. Interestingly, hybrid male misexpression is concentrated mostly in autosomal genes, likely facilitated by the rapid evolution of sex-linked trans -acting factors. I suggest that the faster evolution of X-linked MBGs, at both protein and expression levels, contributes to explain the large effect of the X chromosome on hybrid male sterility, likely mediating widespread autosomal misexpression through the preferential recognition of cis -regulatory elements by conspecific trans -acting factors (i.e., cis–trans conspecific recognition).

Description: The sex chromosomes of the tropical crop papaya ( Carica papaya ) are evolutionarily young and consequently allow for the examination of evolutionary mechanisms that drive early sex chromosome divergence. We conducted a molecular population genetic analysis of four X/Y gene pairs from a collection of 45 wild papaya accessions. These population genetic analyses reveal striking differences in the patterns of polymorphism between the X and Y chromosomes that distinguish them from other sex chromosome systems. In most sex chromosome systems, the Y chromosome displays significantly reduced polymorphism levels, whereas the X chromosome maintains a level of polymorphism that is comparable to autosomal loci. However, the four papaya sex-linked loci that we examined display diversity patterns that are opposite this trend: the papaya X alleles exhibit significantly reduced polymorphism levels, whereas the papaya Y alleles maintain greater than expected levels of diversity. Our analyses suggest that selective sweeps in the regions of the X have contributed to this pattern while also revealing geographically restricted haplogroups on the Y. We discuss the possible role sexual selection and/or genomic conflict have played in shaping the contrasting patterns of polymorphism found for the papaya X and Y chromosomes.

Description: The HLA region shows diversity concerning the number and content of DRB genes present per haplotype. Similar observations are made for the equivalent regions in other primate species. To elucidate the evolutionary history of the various HLA-DRB genes, a large panel of intron sequences obtained from humans, chimpanzees, rhesus macaques, and common marmosets has been subjected to phylogenetic analyses. Special attention was paid to the presence and absence of particular transposable elements and/or to their segments. The sharing of different parts of the same long interspersed nuclear element-2 (LINE2, L2) and various Alu insertions by the species studied demonstrates that one precursor gene must have been duplicated several times before the Old World monkey (OWM) and hominid (HOM) divergence. At least four ancestral DRB gene families appear to have been present before the radiation of OWM and HOM, and one of these even predates the speciation of Old and New World primates. Two of these families represent the pseudogenes DRB6/DRB2 and DRB7 , which have been locked in the genomes of various primate species over long evolutionary time spans. Furthermore, all phylogenies of different intron segments show consistently that, apart from the pseudogenes, only DRB5 genes are shared by OWM and HOM, and they demonstrate the common history of certain DRB genes/lineages of humans and chimpanzees. In contrast, the evolutionary history of some other DRB loci is difficult to decipher, thus illustrating the complex history of the evolution of DRB genes due to a combination of mutations and recombination-like events. The selected approach allowed us to shed light on the ancestral DRB gene pool in primates and on the evolutionary relationship of the various HLA-DRB genes.

Description: Chemosensory genes are frequently the target of positive selection and are often present in large gene families, but little is known about heterogeneity of selection in these cases and its relation to function. Here, we use the vomeronasal-1 receptor ( V1R ) repertoire of mouse lemurs ( Microcebus spp.) as a model system to study patterns of selection of chemosensory genes at several different levels. Mouse lemurs are small nocturnal strepsirrhine primates and have a large (~200 loci) repertoire of V1R loci that are likely important for intraspecific pheromonal communication and interspecific interactions, for example, recognition of predator cues. We investigated signals and patterns of positive selection among the 105 identified full length V1R loci in the gray mouse lemur and within 7 V1R loci amplified across multiple mouse lemur species. Phylogenetic reconstructions of published sequences revealed at least nine monophyletic clusters of V1R s in gray mouse lemurs that have diversified since the split between lemurs and lorisoid primates. A large majority of clusters evolved under significant positive selection. Similar results were found in V1R s of closely related greater galagos. Comparison with function of related V1R clusters in mice suggested a potential relationship between receptor function and strength of selection. Interestingly, most codons identified as being under positive selection are located in the extracellular domains of the receptors and hence likely indicate the position of residues involved in ligand binding. Positive selection was also detected within five V1R loci (=71% of analyzed loci) sequenced from 6 to 10 mouse lemur species, indicating ongoing selection within the genus, which may be related to sexual selection and, potentially, speciation processes. Variation in strength of positive selection on V1R s showed no simple relationship to cluster size. The diversity of V1R loci in mouse lemurs reflects their adaptive evolution and is most likely related to the fundamental relevance of olfactory communication and predator recognition in these primates. Overall, adaptive evolution is the predominant mode of evolution of V1R loci at all levels, and the substantial heterogeneity in the strength of selection may be related to receptor function.

Description: The large number of sexually transmitted diseases and ocular trachoma cases that are caused globally each year by Chlamydia trachomatis has made this organism a World Health Organization priority for vaccine development. However, there is no gene transfer system for Chlamydia to help identify potential vaccine targets. To accelerate discoveries toward this goal, here we analyzed the broadest diversity of C. trachomatis genomes to date, including 25 geographically dispersed clinical and seven reference strains representing 14 of the 19 known serotypes. Strikingly, all 32 genomes were found to have evidence of DNA acquisition by homologous recombination in their history. Four distinct clades were identified, which correspond to all C. trachomatis disease phenotypes: lymphogranuloma venereum (LGV; Clade 1); noninvasive urogenital infections (Clade 2); ocular trachoma (Clade 3); and protocolitis (Clade 4; also includes some noninvasive urogenital infections). Although the ancestral relationship between clades varied, most strains acted as donor and recipient of recombination with no evidence for barriers to genetic exchange. The niche-specific LGV and trachoma clades have undergone less recombination, although the opportunity for mixing with strains from other clades that infect the rectal and ocular mucosa, respectively, is evident. Furthermore, there are numerous occasions for gene conversion events through sequential infections at the same anatomic sites. The size of recombinant segments is relatively small (~357 bp) compared with in vitro experiments of various C. trachomatis strains but is consistent with in vitro estimates for other bacterial species including Escherichia coli and Helicobacter pylori . Selection has also played a crucial role during the diversification of the organism. Clade 2 had the lowest nonsynonymous to synonymous ratio (d N /d S ) but the highest effect of recombination, which is consistent with the widespread occurrence of synonymous substitutions in recombined genomic segments. The trachoma Clade 3 had the highest d N /d S estimates, which may be caused by an increased effect of genetic drift from niche specialization and a reduced effective population size. The degree of drift, selection, and recombination in C. trachomatis suggests that the challenge will remain to identify genomic regions that are stable and cross protective for the development of an efficacious vaccine.

Description: We surveyed genetic variation in alr2 , an allodeterminant of the colonial hydroid Hydractinia symbiolongicarpus . We generated cDNA from a sample of 239 Hydractinia colonies collected at Lighthouse Point, Connecticut, and identified 473 alr2 alleles, 198 of which were unique. Rarefaction analysis suggested that the sample was near saturation. Most alleles were rare, with 86% occurring at frequencies of 1% or less. Alleles were highly variable, diverging on average by 18% of the amino acids in a predicted extracellular domain of the molecule. Analysis of 152 full-length alleles confirmed the existence of two structural types, defined by exons 4–8 of the gene. Several residues of the predicted immunoglobulin superfamily-like domains display signatures of positive selection. We also identified 77 unique alr2 pseudogene sequences from 85 colonies. Twenty-seven of these sequences matched expressed alr2 sequences from other colonies. This observation is consistent with pseudogenes contributing to alr2 diversification through sequence donation. A more limited collection of animals was made from a distant, relict population of H. symbiolongicarpus . Sixty percent of the unique sequences identified in this sample were found to match sequences from the Lighthouse Point population. The large number of alr2 alleles, their degree of divergence, the predominance of rare alleles in the population, their persistence over broad spatial and temporal scales, and the signatures of positive selection in multiple residues of the putative recognition domain paint a consistent picture of negative-frequency-dependent selection operating in this system. The genetic diversity observed at alr2 is comparable to that of the most highly polymorphic genetic systems known to date.

Description: The mechanosensory lateral line, found only in fishes and amphibians, is an important sense organ associated with aquatic life. Lateral line patterns differ among teleost, the most diverse vertebrate taxa, hypothetically in response to selective pressures from different aquatic habitats. In this article, we conduct evolutionary genomic analyses of 34 genes associated with lateral line system development in teleosts to elucidate the significance of contrasting evolutionary rates and changes in the protein coding sequences. We find that duplicated copies of these genes are preferentially retained in the teleost genomes and that episodic events of positive selection have occurred in 22 of the 30 postduplication branches. In general, teleost genes evolved at a faster rate relative to their tetrapod counterparts, and the mutation rates of 26 of the 34 genes differed among teleosts and tetrapods. We conclude that following whole genome duplication, evolutionary rates and episodic events of positive selection on the lateral line system development genes might have been one of the factors favoring the subsequent adaptive radiation of teleosts into diverse habitats. These results provide the foundation for further detailed explorations into lateral line system genes and the evolution of diverse phenotypes and adaptations.

Description: In any comparative studies striving to understand the similarities and differences of the living organisms at the molecular genetic level, the crucial first step is to establish the homology (orthology and paralogy) of genes between different organisms. Determination of the homology of genes becomes complicated when the genes have undergone a rapid divergence in sequence or when the involved genes are members of a gene family that has experienced a differential gain or loss of its constituents in different taxonomic groups. Organisms with duplicated genomes such as teleost fishes might have been especially prone to these problems because the functional redundancies provided by the duplicate copies of genes would have allowed a rapid divergence or loss of genes during evolution. In this study, we will demonstrate that much of the ambiguities in the determination of the homology between fish and tetrapod genes resulting from the problems like these can be eliminated by complementing the sequence-based phylogenies with nonsequence information, such as the exon–intron structure of a gene or the composition of a gene’s genomic neighbors. We will use the Tbx6/16 subfamily genes of zebrafish ( tbx6 , tbx16 , tbx24 , and mga genes), which have been well known for the ambiguity of their evolutionary relationships to the Tbx6/16 subfamily genes of tetrapods, as an illustrative example. We will show that, despite the similarity of sequence and expression to the tetrapod Tbx6 genes, zebrafish tbx6 gene is actually a novel T-box gene more closely related to the tetrapod Tbx16 genes, whereas the zebrafish tbx24 gene, hitherto considered to be a novel gene due to the high level of sequence divergence, is actually an ortholog of tetrapod Tbx6 genes. We will also show that, after their initial appearance by the multiplication of a common ancestral gene at the beginning of vertebrate evolution, the Tbx6/16 subfamily of vertebrate T-box genes might have experienced differential losses of member genes in different vertebrate groups and gradual pooling of member gene’s functions in surviving members, which might have prevented the revelation of the true identity of member genes by way of the comparison of sequence and function.

Description: Cyanobacteria are among the most ancient organisms known to have circadian rhythms. The cpmA gene is involved in controlling the circadian output signal. We studied polymorphism and divergence of this gene in six populations of a stress-tolerant cyanobacterium, Chroococcidiopsis sp., sampled in extreme habitats across the globe. Despite high haplotype diversity (0.774), nucleotide diversity of cpmA is very low ( = 0.0034): the gene appears to be even more conserved than housekeeping genes. Even though the populations were sampled thousands kilometers apart, they manifested virtually no genetic differentiation at this locus ( F ST = 0.0228). Using various tests for neutrality, we determined that evolution of cpmA significantly departures from the neutral model and is governed by episodic positive selection.

Description: Environmental drivers of biodiversity can be identified by relating patterns of community similarity to ecological factors. Community variation has traditionally been assessed by considering changes in species composition and more recently by incorporating phylogenetic information to account for the relative similarity of taxa. Here, we describe how an important class of measures including Bray–Curtis, Canberra, and UniFrac can be extended to allow community variation to be computed on a phylogenetic network. We focus on phylogenetic split systems, networks that are produced by the widely used median network and neighbor-net methods, which can represent incongruence in the evolutionary history of a set of taxa. Calculating β diversity over a split system provides a measure of community similarity averaged over uncertainty or conflict in the available phylogenetic signal. Our freely available software, Network Diversity, provides 11 qualitative (presence–absence, unweighted) and 14 quantitative (weighted) network-based measures of community similarity that model different aspects of community richness and evenness. We demonstrate the broad applicability of network-based diversity approaches by applying them to three distinct data sets: pneumococcal isolates from distinct geographic regions, human mitochondrial DNA data from the Indonesian island of Nias, and proteorhodopsin sequences from the Sargasso and Mediterranean Seas. Our results show that major expected patterns of variation for these data sets are recovered using network-based measures, which indicates that these patterns are robust to phylogenetic uncertainty and conflict. Nonetheless, network-based measures of community similarity can differ substantially from measures ignoring phylogenetic relationships or from tree-based measures when incongruent signals are present in the underlying data. Network-based measures provide a methodology for assessing the robustness of β-diversity results in light of incongruent phylogenetic signal and allow β diversity to be calculated over widely used network structures such as median networks.

Description: We implement an isolation with migration model for three species, with migration occurring between two closely related species while an out-group species is used to provide further information concerning gene trees and model parameters. The model is implemented in the likelihood framework for analyzing multilocus genomic sequence alignments, with one sequence sampled from each of the three species. The prior distribution of gene tree topology and branch lengths at every locus is calculated using a Markov chain characterization of the genealogical process of coalescent and migration, which integrates over the histories of migration events analytically. The likelihood function is calculated by integrating over branch lengths in the gene trees (coalescent times) numerically. We analyze the model to study the gene tree-species tree mismatch probability and the time to the most recent common ancestor at a locus. The model is used to construct a likelihood ratio test (LRT) of speciation with gene flow. We conduct computer simulations to evaluate the LRT and found that the test is in general conservative, with the false positive rate well below the significance level. For the test to have substantial power, hundreds of loci are needed. Application of the test to a human–chimpanzee–gorilla genomic data set suggests gene flow around the time of speciation of the human and the chimpanzee.

Description: Although both genotypes with elevated mutation rate (mutators) and mobilization of insertion sequence (IS) elements have substantial impact on genome diversification, their potential interactions are unknown. Moreover, the evolutionary forces driving gradual accumulation of these elements are unclear: Do these elements spread in an initially transposon-free bacterial genome as they enable rapid adaptive evolution? To address these issues, we inserted an active IS 1 element into a reduced Escherichia coli genome devoid of all other mobile DNA. Evolutionary laboratory experiments revealed that IS elements increase mutational supply and occasionally generate variants with especially large phenotypic effects. However, their impact on adaptive evolution is small compared with mismatch repair mutator alleles, and hence, the latter impede the spread of IS-carrying strains. Given their ubiquity in natural populations, such mutator alleles could limit early phase of IS element evolution in a new bacterial host. More generally, our work demonstrates the existence of an evolutionary conflict between mutation-promoting mechanisms.

Description: Muscles are composed of multinucleated muscle fibers with different contractile and physiological properties, which result from specific slow or fast gene expression programs in the differentiated muscle cells. In the zebra fish embryo, the slow program is under the control of Hedgehog signaling from the notochord and floor plate. This pathway activates the expression of the conserved transcriptional repressor, Prdm1 (Blimp1), which in turn represses the fast program and promotes the slow program in adaxial cells of the somite and their descendants. In the mouse embryo, myogenesis is also initiated in the myotomal compartment of the somite, but the slow muscle program is not confined to a specific subset of cells. We now show that Prdm1 is expressed in the first differentiated myocytes of the early myotome from embryonic day (E)9.5–E11.5. During this period, muscle formation depends on the myogenic regulatory factors, Myf5 and Mrf4. In their absence, Prdm1 is not activated, in apparent contrast to zebra fish where Prdm1 is expressed in the absence of Myf5 and MyoD that drive myogenesis in adaxial cells. However, as in zebra fish, Prdm1 expression in the mouse myotome does not occur in the absence of Hedgehog signaling. Analysis of the muscle phenotype of Prdm1 mutant embryos shows that myogenesis appears to proceed normally. Notably, there is no requirement for Prdm1 activation of the slow muscle program in the mouse myotome. Furthermore, the gene for the transcriptional repressor, Sox6, which is repressed by Prdm1 to permit slow muscle differentiation in zebra fish, is not expressed in the mouse myotome. We propose that the lack of functional conservation for mouse Prdm1, that can nevertheless partially rescue the adaxial cells of zebra fish Prdm1 mutants, reflects differences in the evolution of the role of key regulators such as Prdm1 or Sox6, in initiating the onset of the slow muscle program, between teleosts and mammals.

Description: The specific recognition of antigen by T cells is critical to the generation of adaptive immune responses in vertebrates. T cells recognize antigen using a somatically diversified T-cell receptor (TCR). All jawed vertebrates use four TCR chains called α, β, , and , which are expressed as either a αβ or heterodimer. Nonplacental mammals (monotremes and marsupials) are unusual in that their genomes encode a fifth TCR chain, called TCRµ, whose function is not known but is also somatically diversified like the conventional chains. The origins of TCRµ are also unclear, although it appears distantly related to TCR. Recent analysis of avian and amphibian genomes has provided insight into a model for understanding the evolution of the TCR genes in tetrapods that was not evident from humans, mice, or other commonly studied placental (eutherian) mammals. An analysis of the genes encoding the TCR chains in the duckbill platypus revealed the presence of a highly divergent variable (V) gene, indistinguishable from immunoglobulin heavy (IgH) chain V genes (VH) and related to V genes used in TCRµ. They are expressed as part of TCR repertoire (VH) and similar to what has been found in frogs and birds. This, however, is the first time a VH has been found in a mammal and provides a critical link in reconstructing the evolutionary history of TCRµ. The current structure of TCR and TCRµ genes in tetrapods suggests ancient and possibly recurring translocations of gene segments between the IgH and TCR genes, as well as translocations of TCR genes out of the TCRα/ locus early in mammals, creating the TCRµ locus.

Description: Starch is the most widespread and abundant storage carbohydrate in crops and its production is critical to both crop yield and quality. In regard to the starch content in the seeds of crop plants, there is a distinct difference between grasses (Poaceae) and dicots. However, few studies have described the evolutionary pattern of genes in the starch biosynthetic pathway in these two groups of plants. In this study, therefore, an attempt was made to compare evolutionary rate, gene duplication, and selective pattern of the key genes involved in this pathway between the two groups, using five grasses and five dicots as materials. The results showed 1) distinct differences in patterns of gene duplication and loss between grasses and dicots; duplication in grasses mainly occurred before the divergence of grasses, whereas duplication mostly occurred in individual species within the dicots; there is less gene loss in grasses than in dicots, 2) a considerably higher evolutionary rate in grasses than in dicots in most gene families analyzed, and 3) evidence of a different selective pattern between grasses and dicots; positive selection may have occurred asymmetrically in grasses in some gene families, for example, ADP-glucose pyrophosphorylase small subunit. Therefore, we deduced that gene duplication contributes to, and a higher evolutionary rate is associated with, the higher starch content in grasses. In addition, two novel aspects of the evolution of the starch biosynthetic pathway were observed.

Description: From studies investigating the differences in evolutionary rates between genes, gene compactness and gene expression level have been identified as important determinants of gene-level protein evolutionary rate, as represented by nonsynonymous to synonymous substitution rate ( d N / d S ) ratio. However, the causes of exon-level variances in d N / d S are less understood. Here, we use principal component regression to examine to what extent 13 exon features explain the variance in d N , d S , and the d N / d S ratio of human–rhesus macaque or human–mouse orthologous exons. The exon features were grouped into six functional categories: expression features, mRNA splicing features, structural–functional features, compactness features, exon duplicability, and other features, including G + C content and exon length. Although expression features are important for determining d N and d N / d S between exons of different genes, structural–functional features and splicing features explained more of the variance for exons of the same genes. Furthermore, we show that compactness features can explain only a relatively small percentage of variance in exon-level d N or d N / d S in either between-gene or within-gene comparison. By contrast, d S yielded inconsistent results in the human–mouse comparison and the human–rhesus macaque comparison. This inconsistency may suggest rapid evolutionary changes of the mutation landscape in mammals. Our results suggest that between-gene and within-gene variation in d N / d S (and d N ) are driven by different evolutionary forces and that the role of mRNA splicing in causing the variation in evolutionary rates of coding sequences may be underappreciated.

Description: Extensive synonymous codon modification of viral genomes appears to be an effective way of attenuating strains for use as live vaccines. An assumption of this method is that codon changes have individually small effects, such that codon-attenuated viruses will be slow to evolve back to high fitness (and thus to high virulence). The major capsid gene of the bacterial virus T7 was modified to have varying levels of suboptimal synonymous codons in different constructs, and fitnesses declined linearly with the number of changes. Adaptation of the most extreme design, with 182 codon changes, resulted in a slow fitness recovery by standards of previous experimental evolution with this virus, although fitness effects of substitutions were higher than expected from the average effect of an engineered codon modification. Molecular evolution during recovery was modest, and changes evolved both within the modified gene and outside it. Some changes within the modified gene evolved in parallel across replicates, but with no obvious explanation. Overall, the study supports the premise that codon-modified viruses recover fitness slowly, although the evolution is substantially more rapid than expected from the design principle.

Description: In the fungal kingdom, the evolution of mating systems is highly dynamic, varying even among closely related species. Rearrangements in the mating-type ( mat ) locus, which contains the major regulators of sexual development, are expected to underlie the transitions between self-sterility (heterothallism) and self-fertility (homothallism). However, both the genetic mechanisms and the direction of evolutionary transitions in fungal mating systems are under debate. Here, we present new sequences of the mat locus of four homothallic and one heterothallic species of the model genus Neurospora (Ascomycota). By examining the patterns of synteny among these sequences and previously published data, we show that the locus is conserved among heterothallic species belonging to distinct phylogenetic clades, while different gene arrangements characterize the four homothallic species. These results allowed us to ascertain a heterothallic ancestor for the genus, confirming the prediction of the dead-end theory on unidirectional transitions toward selfing. We show that at least four shifts from heterothallism to homothallism have occurred in Neurospora , three of which involve the acquisition of sequences of both mating types into the same haploid genome. We present evidence for two genetic mechanisms allowing these shifts: translocation and unequal crossover. Finally, we identified two novel retrotransposons and suggest that these have played a major role in mating-system transitions, by facilitating multiple rearrangements of the mat locus.

Description: Estimates of the proportion of amino acid substitutions that have been fixed by selection ( α ) vary widely among taxa, ranging from zero in humans to over 50% in Drosophila . This wide range may reflect differences in the efficacy of selection due to differences in the effective population size ( N e ). However, most comparisons have been made among distantly related organisms that differ not only in N e but also in many other aspects of their biology. Here, we estimate α in three closely related lineages of house mice that have a similar ecology but differ widely in N e : Mus musculus musculus ( N e ~ 25,000–120,000), M. m. domesticus ( N e ~ 58,000–200,000), and M. m. castaneus ( N e ~ 200,000–733,000). Mice were genotyped using a high-density single nucleotide polymorphism array, and the proportions of replacement and silent mutations within subspecies were compared with those fixed between each subspecies and an outgroup, Mus spretus . There was significant evidence of positive selection in M. m. castaneus , the lineage with the largest N e , with α estimated to be approximately 40%. In contrast, estimates of α for M. m. domesticus ( α = 13%) and for M. m. musculus ( α = 12 %) were much smaller. Interestingly, the higher estimate of α for M. m. castaneus appears to reflect not only more adaptive fixations but also more effective purifying selection. These results support the hypothesis that differences in N e contribute to differences among species in the efficacy of selection.

Description: Neanderthals have been shown to share more genetic variants with present-day non-Africans than Africans. Recent admixture between Neanderthals and modern humans outside of Africa was proposed as the most parsimonious explanation for this observation. However, the hypothesis of ancient population structure within Africa could not be ruled out as an alternative explanation. We use simulations to test whether the site frequency spectrum, conditioned on a derived Neanderthal and an ancestral Yoruba (African) nucleotide (the doubly conditioned site frequency spectrum [ dcfs ]), can distinguish between models that assume recent admixture or ancient population structure. We compare the simulations to the dcfs calculated from data taken from populations of European, Chinese, and Japanese descent in the Complete Genomics Diversity Panel. Simulations under a variety of plausible demographic parameters were used to examine the shape of the dcfs for both models. The observed shape of the dcfs cannot be explained by any set of parameter values used in the simulations of the ancient structure model. The dcfs simulations for the recent admixture model provide a good fit to the observed dcfs for non-Africans, thereby supporting the hypothesis that recent admixture with Neanderthals accounts for the greater similarity of Neanderthals to non-Africans than Africans.

Description: The Trinidadian pike cichlid ( Crenicichla frenata ) is a major predator of the guppy ( Poecilia reticulata ), a model system for visual ecology research, and visual predation by the pike cichlid is known to select for male guppies with reduced short-wavelength reflectance. However, an early study of the pike cichlid’s visual system suggested a lack of short-wavelength–sensitive cone photoreceptors, a surprising finding as many African cichlids have highly developed short-wavelength vision. In this study, we found evidence for only four expressed cone opsins (LWS, RH2a, SWS2a, and SWS2b), plus one pseudogene (RH2b). Taken together with our microspectrophotometry data, which revealed the presence of three types of cone photoreceptor, including one sensitive to short-wavelength light, this would indicate a broader spectral capacity than previously believed from earlier visual studies of this fish. Relative to the highly diverse African cichlids, however, this Neotropical cichlid appears to have a greatly reduced opsin complement, reflecting both gene loss along the Neotropical lineage (lacking functional RH2b and, possibly, SWS1 opsins) and gene duplication within the African clade (which possesses paralogous RH2aα and RH2aβ opsins). Molecular evolutionary analyses show that positive selection has shaped the SWS2b and RH1 opsins along the Neotropical lineage, which may be indicative of adaptive evolution to alter nonspectral aspects of opsin biology. These results represent the first molecular evolutionary study of visual pigments in a Neotropical cichlid and thus provide a foundation for further study of a morphologically and ecologically diverse clade that has been understudied with respect to the link between visual ecology and diversification.

Description: On an evolutionary time scale, polymorphic alleles are believed to have a short life, persisting at most tens of millions of years even under long-term balancing selection. Here, we report highly diverged trans -species dimorphism of the proteasome subunit beta type 8 ( PSMB8 ) gene, which encodes a catalytic subunit of the immunoproteasome responsible for the generation of peptides presented by major histocompatibility complex (MHC) class I molecules, in lower teleosts including Cypriniformes (zebrafish and loach) and Salmoniformes (trout and salmon), whose last common ancestor dates to 300 Ma. Moreover, phylogenetic analyses indicated that these dimorphic alleles share lineages with two shark paralogous genes, suggesting that these two lineages have been maintained for more than 500 My either as alleles or as paralogs, and that conversion between alleles and paralogs has occurred at least once during vertebrate evolution. Two lineages termed PSMB8A and PSMB8F show an A 31 F substitution that would probably affect their cleaving specificity, and whereas the PSMB8A lineage has been retained by all analyzed jawed vertebrates, the PSMB8F lineage has been lost by most jawed vertebrates except for cartilaginous fish and basal teleosts. However, a possible functional equivalent of the PSMB8F lineage has been revived as alleles within the PSMB8A lineage at least twice during vertebrate evolution in the amphibian Xenopus and teleostean Oryzias species. Dynamic evolution of the PSMB8 polymorphism through long-term persistence, loss, and regaining of dimorphism and conversion between alleles and paralogs implies the presence of strong selective pressure for functional polymorphism of this gene.

Description: Adaptive immune systems are present only in vertebrates. How do all the remaining animals withstand continuous attacks of permanently evolving pathogens? Even in the absence of adaptive immunity, every organism must be able to unambiguously distinguish "self" cells from any imaginable "nonself." Here, we analyzed the function of highly polymorphic gene vCRL1 , which is expressed in follicle and blood cells of Ciona intestinalis , pointing to possible recognition roles either during fertilization or in immune reactions. By using segregation analysis, we demonstrate that vCRL1 locus is not involved in the control of self-sterility. Interestingly, genetic knockdown of vCRL1 in all tissues or specifically in hemocytes results in a drastic developmental arrest during metamorphosis exactly when blood system formation in Ciona normally occurs. Our data demonstrate that vCRL1 gene might be essential for the establishment of a functional blood system in Ciona . Presumably, presence of the vCRL1 receptor on the surface of blood cells renders them as self, whereas any cell lacking it is referred to as nonself and will be consequently destroyed. We propose that individual-specific receptor vCRL1 might be utilized to facilitate somatic self/nonself discrimination.

Description: 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 .

Description: Synonymous codons are widely selected for various biological mechanisms in both prokaryotes and eukaryotes. Recent evidence suggests that microRNA (miRNA) function may affect synonymous codon choices near miRNA target sites. To better understand this, we perform genome-wide analysis on synonymous codon usage around miRNA target sites in four plant genomes. We observed a general trend of increased site accessibility around miRNA target sites in plants. Guanine-cytosine (GC)–poor codons are preferred in the flank region of miRNA target sites. Within-genome analyses show significant variation among miRNA targets in species. GC content of the target gene can partly explain the variation of site accessibility among miRNA targets. miRNA targets in GC-rich genes show stronger selection signals than those in GC-poor genes. Gene’s codon usage bias and the conservation level of miRNA and its target also have some effects on site accessibility, but the expression level of miRNA or its target and the mechanism of miRNA activity do not contribute to site accessibility differences among miRNA targets. We suggest that synonymous codons near miRNA targets are selected for efficient miRNA binding and proper miRNA function. Our results present a new dimension of natural selection on synonymous codons near miRNA target sites in plants, which will have important implications of coding sequence evolution.

Description: Chloroplasts have evolved from a cyanobacterial endosymbiont and their continuity has been maintained by chloroplast division, which is performed by the constriction of a ring-like division complex at the division site. It is believed that the synchronization of the endosymbiotic and host cell division events was a critical step in establishing a permanent endosymbiotic relationship, such as is commonly seen in existing algae. In the majority of algal species, chloroplasts divide once per specific period of the host cell division cycle. In order to understand both the regulation of the timing of chloroplast division in algal cells and how the system evolved, we examined the expression of chloroplast division genes and proteins in the cell cycle of algae containing chloroplasts of cyanobacterial primary endosymbiotic origin (glaucophyte, red, green, and streptophyte algae). The results show that the nucleus-encoded chloroplast division genes and proteins of both cyanobacterial and eukaryotic host origin are expressed specifically during the S phase, except for FtsZ in one graucophyte alga. In this glaucophyte alga, FtsZ is persistently expressed throughout the cell cycle, whereas the expression of the nucleus-encoded MinD and MinE as well as FtsZ ring formation are regulated by the phases of the cell cycle. In contrast to the nucleus-encoded division genes, it has been shown that the expression of chloroplast-encoded division genes is not regulated by the host cell cycle. The endosymbiotic gene transfer of minE and minD from the chloroplast to the nuclear genome occurred independently on multiple occasions in distinct lineages, whereas the expression of nucleus-encoded MIND and MINE is regulated by the cell cycle in all lineages examined in this study. These results suggest that the timing of chloroplast division in algal cell cycle is restricted by the cell cycle-regulated expression of some but not all of the chloroplast division genes. In addition, it is suggested that the regulation of each division-related gene was established shortly after the endosymbiotic gene transfer, and this event occurred multiple times independently in distinct genes and in distinct lineages.

Description: During allopolyploid speciation, two divergent nuclear genomes merge, yet only one (usually the maternal) of the two sets of progenitor organellar genomes is maintained. Rubisco (1,5-bisphosphate carboxylase/oxygenase) is composed of nuclear-encoded small subunits (SSUs) and plastome-encoded large subunits (LSUs), providing an ideal system to explore the evolutionary process of cytonuclear accommodation. Here, we take initial steps in this direction, using Gossypium allopolyploids as our model. SSU copies from divergent (5–10 My) progenitor diploids ("A" and "D" genomes) were combined at the time of polyploid formation 1–2 Ma, with the LSU encoded by the maternal A-genome parent. LSU genes from A- and D-genome diploids and AD-genome allopolyploids were sequenced, revealing several nonsynonymous substitutions and suggesting the possibility of differential selection on the nuclear-encoded rbcS partner following allopolyploid formation. Sequence data for the rbcS gene family revealed nonreciprocal homoeologous recombination between A- and D- rbcS homoeologs in all polyploid species but not in a synthetic intergenomic F1 hybrid, demonstrating "gene conversion" during allopolyploid evolution. All progenitor rbcS genes are retained and expressed in the five extant allopolyploid species, but analysis of the leaf transcriptome showed that A-homoeologs are preferentially expressed in both the allopolyploid and hybrid, consistent with the maternal origin of rbcL . Although rbcS genes from both progenitor genomes are expressed, some appear to have experienced mutations that may represent cytonuclear coevolution.

Description: In the genome of Artiodactyla (cow, sheep, pigs, camels, and whales), a major retroposon group originated from a presumable horizontal transfer of BovB, a retrotransposon-like element retroposon, between 52 and 70 million years ago. Since then, BovB retroposons have proliferated and today occupy a quarter of the cow’s genome sequence. The BovB-related short interspersed elements (SINEs) were used for resolving the phylogeny of Bovinae (cows, spiral-horned antelopes, and nilgais) and their relatives. In silico screening of 55,000 intronic retroposon insertions in the cow genome and experimental validation of 126 introns resulted in 29 informative retroposon markers for resolving bovine evolutionary relationships. A transposition-in-transposition analysis identifies three different phases of SINE activity and show how BovB elements have expanded in the cattle genome.

Description: The use of codon substitution models to compare synonymous and nonsynonymous substitution rates is a widely used approach to detecting positive Darwinian selection affecting protein evolution. However, in several recent papers, Hughes and colleagues claim that codon-based likelihood-ratio tests (LRTs) are logically flawed as they lack prior hypotheses and fail to accommodate random fluctuations in synonymous and nonsynonymous substitutions Friedman and Hughes (2007) also used site-based LRTs to analyze 605 gene families consisting of human and mouse paralogues. They found that the outcome of the tests was largely determined by irrelevant factors such as the GC content at the third codon positions and the synonymous rate d S , but not by the nonsynonymous rate d N or the d N / d S ratio, factors that should be related to selection. Here, we reanalyze those data. Contra Friedman and Hughes, we found that the test results are related to sequence length and the average d N / d S ratio. We examine the criticisms of Hughes and suggest that they are based on misunderstandings of the codon models and on statistical errors. Our analyses suggest that codon-based tests are useful tools for comparative analysis of genomic data sets.

Description: Giardia intestinalis is a major cause of waterborne enteric disease in humans. The species is divided into eight assemblages suggested to represent separate Giardia species based on host specificities and the genetic divergence of marker genes. We have investigated whether genome-wide recombination occurs between assemblages using the three available G. intestinalis genomes. First, the relative nonsynonymous substitution rates of the homologs were compared for 4,009 positional homologs. The vast majority of these comparisons indicate genetic isolation without interassemblage recombinations. Only a region of 6 kbp suggests genetic exchange between assemblages A and E, followed by gene conversion events. Second, recombination-detecting software fails to identify within-gene recombination between the different assemblages for most of the homologs. Our results indicate very low frequency of recombination between the syntenic core genes, suggesting that G. intestinalis assemblages are genetically isolated lineages and thus should be viewed as separated Giardia species.

Description: During the early stages of speciation, interspecific gene flow may be impeded by deleterious epistatic interactions in hybrids, which maintain parental allelic combinations at the speciation genes. The resulting semipermeable nature of the barrier to interspecific gene flow provides a valuable framework to identify the genes involved in hybrid mortality or sterility, as well as the evolutionary mechanisms that initially caused their divergence. The two Atlantic eels Anguilla anguilla and A. rostrata are partially isolated sister species that naturally hybridize, but whose genetic basis of postzygotic isolation remains unknown. We collected high-throughput sequencing data from the transcriptomes of 58 individuals and discovered 94 genes showing differentially fixed mutations between species. Evidence for positive selection at nuclear diagnostic genes was obtained using multilocus extensions of the McDonald–Kreitman test with polymorphism data from each species. In contrast, mitochondrial protein-coding genes experienced strong purifying selection and mostly diverged at synonymous sites, except for the mt-atp6 gene, which showed an atypically high nonsynonymous to synonymous rate ratio. Nuclear-encoded protein interactors of the mt-atp6 gene in the ATP synthase complex were significantly overrepresented in the list of nuclear diagnostic genes. Further analysis of resequencing data showed that positive selection has operated at both the mt-atp6 gene and its nuclear interactor atp5c1 . These findings suggest that a cytonuclear incompatibility caused by a disruption of normal ATP synthase function in hybrids contributes to partial reproductive isolation between European and American eels.

Description: Gibbons (Hylobatidae) are small, arboreal apes indigenous to Southeast Asia that diverged from other apes ~15–18 Ma. Extant lineages radiated rapidly 6–10 Ma and are organized into four genera ( Hylobates , Hoolock , Symphalangus , and Nomascus ) consisting of 12–19 species. The use of short interspersed elements (SINEs) as phylogenetic markers has seen recent popularity due to several desirable characteristics: the ancestral state of a locus is known to be the absence of an element, rare potentially homoplasious events are relatively easy to resolve, and samples can be quickly and inexpensively genotyped. During radiation of primates, one particular family of SINEs, the Alu family, has proliferated in primate genomes. Nomascus leucogenys (northern white-cheeked gibbon) sequences were analyzed for repetitive content with RepeatMasker using a custom library. The sequences containing Alu elements identified as members of a gibbon-specific subfamily were then compared with orthologous positions in other primate genomes. A primate phylogenetic panel consisting of 18 primate species, including 13 gibbon species representing all four extant genera, was assayed for all loci, and a total of 125 gibbon-specific Alu insertions were identified. The resulting amplification patterns were used to generate a phylogenetic tree. We demonstrate significant support for Symphalangus as the most basal lineage within the family. Our findings also place Nomascus as a derived lineage, sister to Hoolock , with the Nomascus–Hoolock clade sister to Hylobates . Further, our analysis groups N. leucogenys and Nomascus siki as sister taxa to the exclusion of the other Nomascus species assayed. This study represents the first use of SINEs to determine the genus level phylogenetic relationships within the family Hylobatidae. These relationships have been resolved with robust support at most internal nodes, demonstrating the utility of SINE-based phylogenetic analysis. We postulate that hybridization and rapid radiation may have contributed to the complex and contradictory findings of the previous studies. Our findings will aid in the conservation of these threatened primates and inform future studies of the biogeographical history and distribution of modern gibbon species.

Description: Translation termination is accomplished by proteins of the Class I release factor family (RF) that recognize stop codons and catalyze the ribosomal release of the newly synthesized peptide. Bacteria have two canonical RFs: RF1 recognizes UAA and UAG, RF2 recognizes UAA and UGA. Despite that these two release factor proteins are sufficient for de facto translation termination, the eukaryotic organellar RF protein family, which has evolved from bacterial release factors, has expanded considerably, comprising multiple subfamilies, most of which have not been functionally characterized or formally classified. Here, we integrate multiple sources of information to analyze the remarkable differentiation of the RF family among organelles. We document the origin, phylogenetic distribution and sequence structure features of the mitochondrial and plastidial release factors: mtRF1a, mtRF1, mtRF2a, mtRF2b, mtRF2c, ICT1, C12orf65, pRF1, and pRF2, and review published relevant experimental data. The canonical release factors (mtRF1a, mtRF2a, pRF1, and pRF2) and ICT1 are derived from bacterial ancestors, whereas the others have resulted from gene duplications of another release factor. These new RF family members have all lost one or more specific motifs relevant for bona fide release factor function but are mostly targeted to the same organelle as their ancestor. We also characterize the subset of canonical release factor proteins that bear nonclassical PxT/SPF tripeptide motifs and provide a molecular-model-based rationale for their retained ability to recognize stop codons. Finally, we analyze the coevolution of canonical RFs with the organellar genetic code. Although the RF presence in an organelle and its stop codon usage tend to coevolve, we find three taxa that encode an RF2 without using UGA stop codons, and one reverse scenario, where mamiellales green algae use UGA stop codons in their mitochondria without having a mitochondrial type RF2. For the latter, we put forward a "stop-codon reinvention" hypothesis that involves the retargeting of the plastid release factor to the mitochondrion.

Description: The population history of the indigenous populations in island Southeast Asia is generally accepted to have been shaped by two major migrations: the ancient "Out of Africa" migration ~50,000 years before present (YBP) and the relatively recent "Out of Taiwan" expansion of Austronesian agriculturalists approximately 5,000 YBP. The Negritos are believed to have originated from the ancient migration, whereas the majority of island Southeast Asians are associated with the Austronesian expansion. We determined 86 mitochondrial DNA (mtDNA) complete genome sequences in four indigenous Malaysian populations, together with a reanalysis of published autosomal single-nucleotide polymorphism (SNP) data of Southeast Asians to test the plausibility and impact of those migration models. The three Austronesian groups (Bidayuh, Selatar, and Temuan) showed high frequencies of mtDNA haplogroups, which originated from the Asian mainland ~30,000–10,000 YBP, but low frequencies of "Out of Taiwan" markers. Principal component analysis and phylogenetic analysis using autosomal SNP data indicate a dichotomy between continental and island Austronesian groups. We argue that both the mtDNA and autosomal data suggest an "Early Train" migration originating from Indochina or South China around the late-Pleistocene to early-Holocene period, which predates, but may not necessarily exclude, the Austronesian expansion.

Description: Haploid genomes greater than 25,000 Mb are rare, within the animals only the lungfish and some of the salamanders and crustaceans are known to have genomes this large. There is very little data on the structure of genomes this size. It is known, however, that for animal genomes up to 3,000 Mb, there is in general a good correlation between genome size and the percent of the genome composed of repetitive sequence and that this repetitive component is highly dynamic. In this study, we sampled the Australian lungfish genome using three mini-genomic libraries and found that with very little sequence, the results converged on an estimate of 40% of the genome being composed of recognizable transposable elements (TEs), chiefly from the CR1 and L2 long interspersed nuclear element clades. We further characterized the CR1 and L2 elements in the lungfish genome and show that although most CR1 elements probably represent recent amplifications, the L2 elements are more diverse and are more likely the result of a series of amplifications. We suggest that our sampling method has probably underestimated the recognizable TE content. However, on the basis of the most likely sources of error, we suggest that this very large genome is not largely composed of recently amplified, undetected TEs but may instead include a large component of older degenerate TEs. Based on these estimates, and on Thomson’s (Thomson K. 1972. An attempt to reconstruct evolutionary changes in the cellular DNA content of lungfish. J Exp Zool. 180:363–372) inference that in the lineage leading to the extant Australian lungfish, there was massive increase in genome size between 350 and 200 mya, after which the size of the genome changed little, we speculate that the very large Australian lungfish genome may be the result of a massive amplification of TEs followed by a long period with a very low rate of sequence removal and some ongoing TE activity.

Description: L1 elements are mammalian non–long terminal repeat retrotransposons, or long interspersed elements (LINEs), that significantly influence the dynamics and fluidity of the genome. A series of observations suggest that plant L1-clade LINEs, just as mammalian L1s, mobilize both short interspersed elements (SINEs) and certain messenger RNA by recognizing the 3'-poly(A) tail of RNA. However, one L1 lineage in monocots was shown to possess a conserved 3'-end sequence with a solid RNA structure also observed in maize and sorghum SINEs. This strongly suggests that plant LINEs require a particular 3'-end sequence during initiation of reverse transcription. As one L1-clade LINE was also found to share the 3'-end sequence with a SINE in a green algal genome, I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized the specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Description: Plant artificial micro-RNAs (amiRs) have been engineered to target viral genomes and induce their degradation. However, the exceptional evolutionary plasticity of RNA viruses threatens the durability of the resistance conferred by these amiRs. It has recently been shown that viral populations not experiencing strong selective pressure from an antiviral amiR may already contain enough genetic variability in the target sequence to escape plant resistance in an almost deterministic manner. Furthermore, it has also been shown that viral populations exposed to subinhibitory concentrations of the antiviral amiR speed up this process. In this article, we have characterized the molecular evolutionary dynamics of an amiR target sequence in a viral genome under both conditions. The use of Illumina ultradeep sequencing has allowed us to identify virus sequence variants at frequencies as low as 2 x 10 –6 and to track their variation in time before and after the viral population was able of successfully infecting plants fully resistant to the ancestral virus. We found that every site in the amiR-target sequence of the viral genome presented variation and that the variant that eventually broke resistance was sampled among the many coexisting ones. In this system, viral evolution in fully susceptible plants results from an equilibrium between mutation and genetic drift, whereas evolution in partially resistant plants originates from more complex dynamics involving mutation, selection, and drift.

Description: Molecular evolutionary rate estimates have been shown to depend on the time period over which they are estimated. Factors such as demographic processes, calibration errors, purifying selection, and the heterogeneity of substitution rates among sites (RHAS) are known to affect the accuracy with which rates of evolution are estimated. We use mathematical modeling and Bayesian analyses of simulated sequence alignments to explore how mutational hotspots can lead to time-dependent rate estimates. Mathematical modeling shows that underestimation of molecular rates over increasing time scales is inevitable when RHAS is ignored. Although a gamma distribution is commonly used to model RHAS, we show that when the actual RHAS deviates from a gamma-like distribution, rates can either be under- or overestimated in a time-dependent manner. Simulations performed under different scenarios of RHAS confirm the mathematical modeling and demonstrate the impacts of time-dependent rates on estimates of divergence times. Most notably, erroneous rate estimates can have narrow credibility intervals, leading to false confidence in biased estimates of rates, and node ages. Surprisingly, large errors in estimates of overall molecular rate do not necessarily generate large errors in divergence time estimates. Finally, we illustrate the correlation between time-dependent rate patterns and differential saturation between quickly and slowly evolving sites. Our results suggest that data partitioning or simple nonparametric mixture models of RHAS significantly improve the accuracy with which node ages and substitution rates can be estimated.

Description: Multiple visual pigments, prerequisites for color vision, are found in arthropods, but the evolutionary origin of their diversity remains obscure. In this study, we explore the opsin genes in five distantly related species of Onychophora, using deep transcriptome sequencing and screening approaches. Surprisingly, our data reveal the presence of only one opsin gene (onychopsin) in each onychophoran species, and our behavioral experiments indicate a maximum sensitivity of onychopsin to blue–green light. In our phylogenetic analyses, the onychopsins represent the sister group to the monophyletic clade of visual r-opsins of arthropods. These results concur with phylogenomic support for the sister-group status of the Onychophora and Arthropoda and provide evidence for monochromatic vision in velvet worms and in the last common ancestor of Onychophora and Arthropoda. We conclude that the diversification of visual pigments and color vision evolved in arthropods, along with the evolution of compound eyes—one of the most sophisticated visual systems known.

Description: Cis- regulatory DNA has been suspected to play a preeminent role in adaptive evolution, but understanding the role of cis -regulatory mutations in gene expression divergence first requires an accurate analysis of the functional differences associated with these regions. We analyzed allele-specific expression (ASE) in leaf and floral tissues of F1 interspecific hybrids generated between the two closely related species Arabidopsis thaliana and A. lyrata with a whole-genome SNP (single nucleotide polymorphism) tiling array. We observed 2,205 genes showing ASE pattern in at least one tissue. Nearly 90% of genes displaying ASE preferentially expressed the allele of A. lyrata . Genome-wide comparison of sequence divergence revealed that genes displaying ASE had a higher ratio of nonsynonymous to synonymous substitutions in coding regions. We further observe that the epigenetic landscape of histone methylation in A. thaliana genome associate with ASE. The asymmetry in the direction of allele-specific expression suggests interspecific differences in the efficiency of gene silencing in F1 hybrids.

Description: The olfactory receptor (OR) gene family is the largest gene family found in mammalian genomes. It is known to evolve through a birth-and-death process. Here, we characterized the sequences of 16 segregating OR pseudogenes in the samples of the wolf and the Chinese village dog (CVD) and compared them with the sequences from dogs of different breeds. Our results show that the segregating OR pseudogenes in breed dogs are under strong purifying selection, while evolving neutrally in the CVD, and show a more complicated pattern in the wolf. In the wolf, we found a trend to remove deleterious polymorphisms and accumulate nondeleterious polymorphisms. On the basis of protein structure of the ORs, we found that the distribution of different types of polymorphisms (synonymous, nonsynonymous, tolerated, and untolerated) varied greatly between the wolf and the breed dogs. In summary, our results suggest that different forms of selection have acted on the segregating OR pseudogenes in the CVD since domestication, breed dogs after breed formation, and ancestral wolf population, which has driven the evolution of these genes in different directions.

Description: Although gene duplications occur at a higher rate, only a small fraction of these are retained. The position of a gene’s encoded product in the protein–protein interaction network has recently emerged as a determining factor of gene duplicability. However, the direction of the relationship between network centrality and duplicability is not universal: In Escherichia coli , yeast, fly, and worm, duplicated genes more often act at the periphery of the network, whereas in humans, such genes tend to occupy the most central positions. Herein, we have inferred duplication events that took place in the different branches of the primate phylogeny. In agreement with previous observations, we found that duplications generally affected the most central network genes, which is presumably the process that has most influenced the trend in humans. However, the opposite trend—that is, duplication being more common in genes whose encoded products are peripheral in the network—is observed for three recent branches, including, quite counterintuitively, the external branch leading to humans. This indicates a shift in the relationship between centrality and duplicability during primate evolution. Furthermore, we found that genes encoding interacting proteins exhibit phylogenetic tree topologies that are more similar than expected for random pairs and that genes duplicated in a given branch of the phylogeny tend to interact with those that duplicated in the same lineage. These results indicate that duplication of a gene increases the likelihood of duplication of its interacting partners. Our observations indicate that the structure of the primate protein–protein interaction network affects gene duplicability in previously unrecognized ways.

Description: An absolute timescale for evolution is essential if we are to associate evolutionary phenomena, such as adaptation or speciation, with potential causes, such as geological activity or climatic change. Timescales in most phylogenetic studies use geologically dated fossils or phylogeographic events as calibration points, but more recently, it has also become possible to use experimentally derived estimates of the mutation rate as a proxy for substitution rates. The large radiation of drosophilid taxa endemic to the Hawaiian islands has provided multiple calibration points for the Drosophila phylogeny, thanks to the "conveyor belt" process by which this archipelago forms and is colonized by species. However, published date estimates for key nodes in the Drosophila phylogeny vary widely, and many are based on simplistic models of colonization and coalescence or on estimates of island age that are not current. In this study, we use new sequence data from seven species of Hawaiian Drosophila to examine a range of explicit coalescent models and estimate substitution rates. We use these rates, along with a published experimentally determined mutation rate, to date key events in drosophilid evolution. Surprisingly, our estimate for the date for the most recent common ancestor of the genus Drosophila based on mutation rate (25–40 Ma) is closer to being compatible with independent fossil-derived dates (20–50 Ma) than are most of the Hawaiian-calibration models and also has smaller uncertainty. We find that Hawaiian-calibrated dates are extremely sensitive to model choice and give rise to point estimates that range between 26 and 192 Ma, depending on the details of the model. Potential problems with the Hawaiian calibration may arise from systematic variation in the molecular clock due to the long generation time of Hawaiian Drosophila compared with other Drosophila and/or uncertainty in linking island formation dates with colonization dates. As either source of error will bias estimates of divergence time, we suggest mutation rate estimates be used until better models are available.

Description: We used next-generation sequencing to characterize the genomes of nine species of Orobanchaceae of known phylogenetic relationships, different life forms, and including a polyploid species. The study species are the autotrophic, nonparasitic Lindenbergia philippensis , the hemiparasitic Schwalbea americana , and seven nonphotosynthetic parasitic species of Orobanche ( Orobanche crenata , Orobanche cumana , Orobanche gracilis (tetraploid), and Orobanche pancicii ) and Phelipanche ( Phelipanche lavandulacea , Phelipanche purpurea , and Phelipanche ramosa ). Ty3/Gypsy elements comprise 1.93%–28.34% of the nine genomes and Ty1/Copia elements comprise 8.09%–22.83%. When compared with L. philippensis and S. americana , the nonphotosynthetic species contain higher proportions of repetitive DNA sequences, perhaps reflecting relaxed selection on genome size in parasitic organisms. Among the parasitic species, those in the genus Orobanche have smaller genomes but higher proportions of repetitive DNA than those in Phelipanche , mostly due to a diversification of repeats and an accumulation of Ty3/Gypsy elements. Genome downsizing in the tetraploid O. gracilis probably led to sequence loss across most repeat types.

Description: Burkina Faso is located in the heart of West Africa and is a representative of the local structured patterns of human variability. Here, different cultures and languages are found in a geographic contiguity, as a result of several waves of migration and the succession of long- and short-term empires. However, historical documentation for this area is only partial, focusing predominantly on the recent empires, and linguistic surveys lack the power to fully elucidate the social context of the contact-induced changes. In this paper, we report Y-chromosomal data and complete mtDNA genome sequences for ten populations from Burkina Faso whose languages belong to two very distantly related branches of the Niger–Congo phylum, the Gur and Mande language families. In addition, two further populations, the Mande-speaking Mandenka from Senegal and the Yoruba from Nigeria, were included for regional comparison. We focus on the different historical trajectories undergone by the maternal and paternal lineages. Our results reveal a striking structure in the paternal line, which matches the linguistic affiliation of the ethnolinguistic groups, in contrast to the near-complete homogeneity of the populations in the maternal line. However, while the ancient structure along the linguistic lines is apparent in the Y-chromosomal haplogroup affiliation, this has clearly been overlain by more recent migrations, as shown by significant correlations between the genetic distances based on Y chromosome short tandem repeats and geographic distances between the populations, as well as by the patterns of shared haplotypes. Using the complete mtDNA sequences, we are able to reconstruct population size variation in the past, showing a strong sign of expansion in the concomitance with the Holocene Climate Optimum approximately 12,000–10,000 years ago, which has been suggested as the cause of the spread of the Niger–Congo phylum in the area. However, subsequent climatic fluctuations do not appear to have had an impact on the demography of the inhabitants of West Africa, probably reflecting the adaptive advantages of cultural innovations, such as pastoralism and agriculture.

Description: In computational evolutionary biology, verification and benchmarking is a challenging task because the evolutionary history of studied biological entities is usually not known. Computer programs for simulating sequence evolution in silico have shown to be viable test beds for the verification of newly developed methods and to compare different algorithms. However, current simulation packages tend to focus either on gene-level aspects of genome evolution such as character substitutions and insertions and deletions (indels) or on genome-level aspects such as genome rearrangement and speciation events. Here, we introduce Artificial Life Framework (ALF), which aims at simulating the entire range of evolutionary forces that act on genomes: nucleotide, codon, or amino acid substitution (under simple or mixture models), indels, GC-content amelioration, gene duplication, gene loss, gene fusion, gene fission, genome rearrangement, lateral gene transfer (LGT), or speciation. The other distinctive feature of ALF is its user-friendly yet powerful web interface. We illustrate the utility of ALF with two possible applications: 1) we reanalyze data from a study of selection after globin gene duplication and test the statistical significance of the original conclusions and 2) we demonstrate that LGT can dramatically decrease the accuracy of two well-established orthology inference methods. ALF is available as a stand-alone application or via a web interface at http://www.cbrg.ethz.ch/alf .

Description: Although human bitter taste perception is hypothesized to be a dietary adaptation, little is known about genetic signatures of selection and patterns of bitter taste perception variability in ethnically diverse populations with different diets, particularly from Africa. To better understand the genetic basis and evolutionary history of bitter taste sensitivity, we sequenced a 2,975 bp region encompassing TAS2R38 , a bitter taste receptor gene, in 611 Africans from 57 populations in West Central and East Africa with diverse subsistence patterns, as well as in a comparative sample of 132 non-Africans. We also examined the association between genetic variability at this locus and threshold levels of phenylthiocarbamide (PTC) bitterness in 463 Africans from the above populations to determine how variation influences bitter taste perception. Here, we report striking patterns of variation at TAS2R38 , including a significant excess of novel rare nonsynonymous polymorphisms that recently arose only in Africa, high frequencies of haplotypes in Africa associated with intermediate bitter taste sensitivity, a remarkably similar frequency of common haplotypes across genetically and culturally distinct Africans, and an ancient coalescence time of common variation in global populations. Additionally, several of the rare nonsynonymous substitutions significantly modified levels of PTC bitter taste sensitivity in diverse Africans. While ancient balancing selection likely maintained common haplotype variation across global populations, we suggest that recent selection pressures may have also resulted in the unusually high level of rare nonsynonymous variants in Africa, implying a complex model of selection at the TAS2R38 locus in African populations. Furthermore, the distribution of common haplotypes in Africa is not correlated with diet, raising the possibility that common variation may be under selection due to their role in nondietary biological processes. In addition, our data indicate that novel rare mutations contribute to the phenotypic variance of PTC sensitivity, illustrating the influence of rare variation on a common trait, as well as the relatively recent evolution of functionally diverse alleles at this locus.

Description: Evolution of floral monosymmetry is thought to be a major driving force of angiosperm radiation, making angiosperms the most successful land plant group in terms of species richness. Monosymmetry evolved from a polysymmetric ancestor repeatedly in different angiosperm lineages, where it likely facilitated diversification through the interaction with insects. Most monosymmetric taxa are thus dominated by monosymmetric members. However, in the Brassicaceae, only few members develop a monosymmetric corolla with two petal pairs of unequal size, making them an ideal system to study the evolution of molecular mechanisms enhancing flower complexity. Monosymmetry is controlled by the TCP transcription factors that belong to the CYC2 clade in distantly related taxa. In Iberis amara , the first crucifer analyzed in terms of monosymmetry development, unequal corolla formation is due to a stronger CYC2 clade gene expression in the smaller adaxial petals compared with the larger abaxial ones. Phylogenetic reconstruction of the crucifer family reveals that the monosymmetric genera Iberis , Calepina , and Teesdalia belong to one major crucifer lineage. Monosymmetry is most pronounced in Iberis and less so in Calepina and Teesdalia , with a positive dosage-dependent correlation between the strength of a CYC2 expression difference and the extent of monosymmetry formation. An early adaxial CYC2 expression in floral meristems, observed in many distantly related taxa, might have facilitated the repeated evolution of CYC2 -controlled monosymmetry. Comparison of early and late CYC2 expression in monosymmetric and polysymmetric crucifers representative for the four major crucifer lineages reveals that an adaxial CYC2 expression in floral meristems is likely ancestral for the Brassicaceae. However, it got lost in all analyzed monosymmetric members and is, as such, not a prerequisite for the establishment of corolla monosymmetry in crucifers. Here, monosymmetry evolved via a heterochronic CYC2 expression shift from an ancestral early adaxial expression in floral meristems to an adaxial CYC2 transcript accumulation later in petal development. This study emphasizes the potential of regulatory changes in the evolution of morphological novelties, like corolla monosymmetry in the Brassicaceae. In combination with a corymboid inflorescence, monosymmetry might have served as a key invention driving diversification in the genus Iberis comprising more than 20 monosymmetric species.

Description: Plant-specific TCP domain transcription factors have been shown to regulate morphological novelties during plant evolution, including the complex architecture of the Asteraceae inflorescence that involves different types of flowers. We conducted comparative analysis of the CYCLOIDEA/TEOSINTE BRANCHED1 ( CYC/TB1 ) gene family in Gerbera hybrida (gerbera) and Helianthus annuus (sunflower), two species that represent distant tribes within Asteraceae. Our data confirm that the CYC/TB1 gene family has expanded in Asteraceae, a condition that appears to be connected with the increased developmental complexity and evolutionary success of this large plant family. Phylogenetic analysis of the CYC/TB1 gene family revealed both shared and lineage-specific duplications in gerbera and sunflower, corresponding to the three gene lineages previously identified as specific to core eudicots: CYC1, CYC2, and CYC3. Expression analyses of early stages of flower primordia development indicated that especially within the CYC2 clade, with the greatest number of secondary gene duplications, gene expression patterns are conserved between the species and associated with flower and inflorescence development. All sunflower and gerbera CYC2 clade genes showed differential expression between developing flower types, being upregulated in marginal ray (and trans) flowers. One gene in gerbera ( GhCYC3 ) and two in sunflower ( HaCYC2d and HaCYC2c ) were indicated to be strong candidates as regulators of ray flower identity, a function that is specific for Asteraceae. Our data further showed that other CYC2 clade genes are likely to have more specialized functions at the level of single flowers, including the late functions in floral reproductive organs that may be more conserved across plant families. The expression patterns of CYC1 and CYC3 clade genes showed more differences between the two species but still pointed to possible conserved functions during vegetative plant development. Pairwise protein–protein interaction assays gave the first molecular evidence that CYC/TB1-like proteins function in complexes. Compared with sunflower, the gerbera proteins showed higher capacity for dimerization, between as well as within CYC clades. Our data from two distant species within the Asteraceae suggest that the expansion and the apparent conservation of especially the CYC2 clade CYC/TB1 -like genes are associated with the evolution of the increased complexity of the Asteraceae inflorescence architecture.

Description: The plant domestication process is associated with considerable modifications of plant phenotype. The identification of the genetic basis of this adaptation is of great interest for evolutionary biology. One of the methods used to identify such genes is the detection of signatures of selection. However, domestication is generally associated with major demographic effects. It is therefore crucial to disentangle the effects of demography and selection on diversity. In this study, we investigated selection in a flowering time pathway during domestication of pearl millet. We first used a random set of 20 genes to model pearl millet domestication using approximate Bayesian computation. This analysis showed that a model with exponential growth and wild–cultivated gene flow was well supported by our data set. Under this model, the domestication date of pearl millet is estimated at around 4,800 years ago. We assessed selection in 15 pearl millet DNA sequences homologous to flowering time genes and showed that these genes underwent selection more frequently than expected. We highlighted significant signatures of selection in six pearl millet flowering time genes associated with domestication or improvement of pearl millet. Moreover, higher deviations from neutrality were found for circadian clock–associated genes. Our study provides new insights into the domestication process of pearl millet and shows that a category of genes of the flowering pathway were preferentially selected during pearl millet domestication.

Description: Plasmids are vessels of genetic exchange in microbial communities. They are known to transfer between different host organisms and acquire diverse genetic elements from chromosomes and/or other plasmids. Therefore, they constitute an important element in microbial evolution by rapidly disseminating various genetic properties among different communities. A paradigmatic example of this is the dissemination of antibiotic resistance (AR) genes that has resulted in the emergence of multiresistant pathogenic bacterial strains. To globally analyze the evolutionary dynamics of plasmids, we built a large graph in which 2,343 plasmids (nodes) are connected according to the proteins shared by each other. The analysis of this gene-sharing network revealed an overall coherence between network clustering and the phylogenetic classes of the corresponding microorganisms, likely resulting from genetic barriers to horizontal gene transfer between distant phylogenetic groups. Habitat was not a crucial factor in clustering as plasmids from organisms inhabiting different environments were often found embedded in the same cluster. Analyses of network metrics revealed a statistically significant correlation between plasmid mobility and their centrality within the network, providing support to the observation that mobile plasmids are particularly important in spreading genes in microbial communities. Finally, our study reveals an extensive (and previously undescribed) sharing of AR genes between Actinobacteria and Gammaproteobacteria, suggesting that the former might represent an important reservoir of AR genes for the latter.

Description: Although it is thought that changes in gene expression play an important role in adaptation, the identification of gene-regulatory sequences that have been targets of positive selection has proved difficult. Here, we identify a cis -regulatory element of the Drosophila melanogaster CG9509 gene that is associated with a selective sweep in a derived non-African population of the species. Expression analyses indicate that CG9509 consistently shows greater expression in non-African than in African strains of D . melanogaster . We find that a 1.8 kb region located just upstream of the CG9509 coding region is devoid of DNA sequence polymorphism in a European population sample and that this is best explained by the recent action of positive selection (within the past 4,000–10,000 years). Using a reporter gene construct and phiC31-mediated site-specific integration, we show that the European version of the CG9509 upstream region drives 2–3 times greater expression than the African version in an otherwise identical genetic background. This expression difference corresponds well to that of the native gene and indicates that sequence variation within the CG9509 upstream region can completely account for its high expression in the European population. Selection appears to have favored a quantitative increase in gene expression in the Malphigian tubule, the tissue where CG9509 is predominantly expressed.

Description: In accordance with the hypothesis that cancer formation is a process of somatic evolution driven by natural selection, signature of positive selection has been detected on a number of cancer-related nuclear genes. It remains, however, controversial whether a similar selective pressure has also acted on mitochondrial DNA (mtDNA), a small molecule in mitochondrion that may play an important role in tumorigenesis by altering oxidative phosphorylation. To better understand the mutational pattern on cancerous mtDNA and decipher the genetic signature left by natural selection, a total of 186 entire mitochondrial genomes of cancerous and adjacent normal tissues from 93 esophageal cancer patients were obtained and extensively studied. Our results revealed that the observed mutational pattern on the cancerous mtDNAs might be best explained as relaxation of negative selection. Taking into account an additional 1,235 cancerous (nearly) complete mtDNA sequences retrieved from the literature, our results suggested that the relaxed selective pressure was the most likely explanation for the accumulation of mtDNA variation in different types of cancer. This notion is in good agreement with the observation that aerobic glycolysis, instead of mitochondrial respiration, plays the key role in generating energy in cancer cells. Furthermore, our study provided solid evidence demonstrating that problems in some of the published cancerous mtDNA data adequately explained the previously contradictory conclusions about the selective pressure on cancer mtDNA, thus serving as a paradigm emphasizing the importance of data quality in affecting our understanding on the role of mtDNA in tumorigenesis.

Description: Pathogen recognition is the first step of immune reactions. In animals and plants, direct or indirect pathogen recognition is often mediated by a wealth of fast-evolving receptors, many of which contain ligand-binding and signal transduction domains, such as leucine-rich or tetratricopeptide repeat (LRR/TPR) and NB-ARC domains, respectively. In order to identify candidates potentially involved in algal defense, we mined the genome of the brown alga Ectocarpus siliculosus for homologues of these genes and assessed the evolutionary pressures acting upon them. We thus annotated all Ectocarpus LRR-containing genes, in particular an original group of LRR-containing GTPases of the ROCO family, and 24 NB-ARC–TPR proteins. They exhibit high birth and death rates, while a diversifying selection is acting on their LRR (respectively TPR) domain, probably affecting the ligand-binding specificities. Remarkably, each repeat is encoded by an exon, and the intense exon shuffling underpins the variability of LRR and TPR domains. We conclude that the Ectocarpus ROCO and NB-ARC–TPR families are excellent candidates for being involved in recognition/transduction events linked to immunity. We further hypothesize that brown algae may generate their immune repertoire via controlled somatic recombination, so far only known from the vertebrate adaptive immune systems.

Description: By exploiting the large body of genome data and the considerable progress in phylogenetic methodology, recent phylogenomic studies have provided new insights into the relationships among major eukaryotic groups. However, confident placement of the eukaryotic root remains a major challenge. This is due to the large evolutionary distance separating eukaryotes from their closest relatives, the Archaea, implying a weak phylogenetic signal and strong long-branch attraction artifacts. Here, we apply a new approach to the rooting of the eukaryotic tree by using a subset of genomic information with more recent evolutionary origin—mitochondrial sequences, whose closest relatives are α-Proteobacteria. For this, we identified and assembled a data set of 42 mitochondrial proteins (mainly encoded by the nuclear genome) and performed Bayesian and maximum likelihood analyses. Taxon sampling includes the recently sequenced Thecamonas trahens , a member of the phylogenetically elusive Apusozoa. This data set confirms the relationships of several eukaryotic supergroups seen before and places the eukaryotic root between the monophyletic "unikonts" and "bikonts." We further show that T. trahens branches sister to Opisthokonta with significant statistical support and question the bikont/excavate affiliation of Malawimonas species. The mitochondrial data set developed here (to be expanded in the future) constitutes a unique alternative means in resolving deep eukaryotic relationships.

Description: Oligoadenylate synthetases (OAS) are interferon-induced enzymes that participate in the first line of defense against a wide range of viral infection in animals. Upon activation by viral double-stranded RNA, OAS synthesizes (2–5) oligoadenylates, which activate RNase L, leading to the nonspecific degradation of cellular and viral RNA. Some association studies in humans suggest that variation at one of the OAS genes, OAS1 , could be influencing host susceptibility to viral infection. We assessed the diversity of OAS1 in hominoid primates with a focus on chimpanzees. We found that the OAS1 gene is extremely polymorphic in Central African chimpanzee and exhibits levels of silent and replacement diversity much higher than neutral regions of the chimpanzee genome. This level of variation strongly suggests that balancing selection is acting on OAS1 , and indeed, this conclusion was validated by several tests of neutrality. We further demonstrated that balancing selection has been acting at this locus since the split between chimpanzees, humans, and gorillas (~8.6 Ma) and caused the persistence of two deeply divergent allelic lineages in Central African chimpanzees. These two groups of OAS1 alleles differ by a large number of amino acids (a.a.), including several a.a. putatively involved in RNA binding. It is therefore very likely that variation at the OAS1 locus affects the innate immune response of individuals to specific viral infection. Our data strongly suggest that interactions between viral RNA and OAS1 are responsible for the maintenance of ancestral polymorphisms at this locus for at least 13.2 My.