Showing papers in "Genome Biology and Evolution in 2014"

A Phylogenomic View of Ecological Specialization in the Lachnospiraceae, a Family of Digestive Tract-Associated Bacteria

Abstract: Several bacterial families are known to be highly abundant within the human microbiome, but their ecological roles and evolutionary histories have yet to be investigated in depth. One such family, Lachnospiraceae (phylum Firmicutes, class Clostridia) is abundant in the digestive tracts of many mammals and relatively rare elsewhere. Members of this family have been linked to obesity and protection from colon cancer in humans, mainly due to the association of many species within the group with the production of butyric acid, a substance that is important for both microbial and host epithelial cell growth. We examined the genomes of 30 Lachnospiraceae isolates to better understand the origin of butyric acid capabilities and other ecological adaptations within this group. Butyric acid production-related genes were detected in fewer than half of the examined genomes with the distribution of this function likely arising in part from lateral gene transfer (LGT). An investigation of environment-specific functional signatures indicated that human gut-associated Lachnospiraceae possess genes for endospore formation, whereas other members of this family lack key sporulation-associated genes, an observation supported by analysis of metagenomes from the human gut, oral cavity, and bovine rumen. Our analysis demonstrates that adaptation to an ecological niche and acquisition of defining functional roles within a microbiome can arise through a combination of both habitat-specific gene loss and LGT.

An expanded genomic representation of the phylum cyanobacteria.

Abstract: Molecular surveys of aphotic habitats have indicated the presence of major uncultured lineages phylogenetically classified as members of the Cyanobacteria. One of these lineages has recently been proposed as a nonphotosynthetic sister phylum to the Cyanobacteria, the Melainabacteria, based on recovery of population genomes from human gut and groundwater samples. Here, we expand the phylogenomic representation of the Melainabacteria through sequencing of six diverse population genomes from gut and bioreactor samples supporting the inference that this lineage is nonphotosynthetic, but not the assertion that they are strictly fermentative. We propose that the Melainabacteria is a class within the phylogenetically defined Cyanobacteria based on robust monophyly and shared ancestral traits with photosynthetic representatives. Our findings are consistent with theories that photosynthesis occurred late in the Cyanobacteria and involved extensive lateral gene transfer and extends the recognized functionality of members of this phylum.

The Genomic Diversification of the Whole Acinetobacter Genus: Origins, Mechanisms, and Consequences

Abstract: Bacterial genomics has greatly expanded our understanding of microdiversification patterns within a species, but analyses at higher taxonomical levels are necessary to understand and predict the independent rise of pathogens in a genus We have sampled, sequenced, and assessed the diversity of genomes of validly named and tentative species of the Acinetobacter genus, a clade including major nosocomial pathogens and biotechnologically important species We inferred a robust global phylogeny and delimited several new putative species The genus is very ancient and extremely diverse: Genomes of highly divergent species share more orthologs than certain strains within a species We systematically characterized elements and mechanisms driving genome diversification, such as conjugative elements, insertion sequences, and natural transformation We found many error-prone polymerases that may play a role in resistance to toxins, antibiotics, and in the generation of genetic variation Surprisingly, temperate phages, poorly studied in Acinetobacter, were found to account for a significant fraction of most genomes Accordingly, many genomes encode clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems with some of the largest CRISPR-arrays found so far in bacteria Integrons are strongly overrepresented in Acinetobacter baumannii, which correlates with its frequent resistance to antibiotics Our data suggest that A baumannii arose from an ancient population bottleneck followed by population expansion under strong purifying selection The outstanding diversification of the species occurred largely by horizontal transfer, including some allelic recombination, at specific hotspots preferentially located close to the replication terminus Our work sets a quantitative basis to understand the diversification of Acinetobacter into emerging resistant and versatile pathogens

Evolutionary origins and dynamics of octoploid strawberry subgenomes revealed by dense targeted capture linkage maps.

Abstract: Whole-genome duplications are radical evolutionary events that have driven speciation and adaptation in many taxa. Higher-order polyploids have complex histories often including interspecific hybridization and dynamic genomic changes. This chromosomal reshuffling is poorly understood for most polyploid species, despite their evolutionary and agricultural importance, due to the challenge of distinguishing homologous sequences from each other. Here, we use dense linkage maps generated with targeted sequence capture to improve the diploid strawberry (Fragaria vesca) reference genome and to disentangle the subgenomes of the wild octoploid progenitors of cultivated strawberry, Fragaria virginiana and Fragaria chiloensis. Our novel approach, POLiMAPS (Phylogenetics Of Linkage-Map-Anchored Polyploid Subgenomes), leverages sequence reads to associate informative interhomeolog phylogenetic markers with linkage groups and reference genome positions. In contrast to a widely accepted model, we find that one of the four subgenomes originates with the diploid cytoplasm donor F. vesca, one with the diploid Fragaria iinumae, and two with an unknown ancestor close to F. iinumae. Extensive unidirectional introgression has converted F. iinumae-like subgenomes to be more F. vesca-like, but never the reverse, due either to homoploid hybridization in the F. iinumae-like diploid ancestors or else strong selection spreading F. vesca-like sequence among subgenomes through homeologous exchange. In addition, divergence between homeologous chromosomes has been substantially augmented by interchromosomal rearrangements. Our phylogenetic approach reveals novel aspects of the complicated web of genetic exchanges that occur during polyploid evolution and suggests a path forward for unraveling other agriculturally and ecologically important polyploid genomes.

Genome degeneration and adaptation in a nascent stage of symbiosis

Abstract: Symbiotic associations between animals and microbes are ubiquitous in nature, with an estimated 15% of all insect species harboring intracellular bacterial symbionts. Most bacterial symbionts share many genomic features including small genomes, nucleotide composition bias, high coding density, and a paucity of mobile DNA, consistent with long-term host association. In this study, we focus on the early stages of genome degeneration in a recently derived insect-bacterial mutualistic intracellular association. We present the complete genome sequence and annotation of Sitophilus oryzae primary endosymbiont (SOPE). We also present the finished genome sequence and annotation of strain HS, a close free-living relative of SOPE and other insect symbionts of the Sodalis-allied clade, whose gene inventory is expected to closely resemble the putative ancestor of this group. Structural, functional, and evolutionary analyses indicate that SOPE has undergone extensive adaptation toward an insect-associated lifestyle in a very short time period. The genome of SOPE is large in size when compared with many ancient bacterial symbionts; however, almost half of the protein-coding genes in SOPE are pseudogenes. There is also evidence for relaxed selection on the remaining intact protein-coding genes. Comparative analyses of the whole-genome sequence of strain HS and SOPE highlight numerous genomic rearrangements, duplications, and deletions facilitated by a recent expansion of insertions sequence elements, some of which appear to have catalyzed adaptive changes. Functional metabolic predictions suggest that SOPE has lost the ability to synthesize several essential amino acids and vitamins. Analyses of the bacterial cell envelope and genes encoding secretion systems suggest that these structures and elements have become simplified in the transition to a mutualistic association.

High Variability of Mitochondrial Gene Order among Fungi

Abstract: From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Fungal cytochrome P450 monooxygenases: Their distribution, structure, functions, family expansion, and evolutionary origin

Abstract: Cytochrome P450 (CYP) monooxygenase superfamily contributes a broad array of biological functions in living organisms. In fungi, CYPs play diverse and pivotal roles in versatile metabolism and fungal adaptation to specific ecological niches. In this report, CYPomes in the 47 genomes of fungi belong to the phyla Ascomycota, Basidiomycota, Chytridiomycota, and Zygomycota have been studied. The comparison of fungal CYPomes suggests that generally fungi possess abundant CYPs belonging to a variety of families with the two global families CYP51 and CYP61, indicating individuation of CYPomes during the evolution of fungi. Fungal CYPs show highly conserved characteristic motifs, but very low overall sequence similarities. The characteristic motifs of fungal CYPs are distinguishable from those of CYPs in animals, plants, and especially archaea and bacteria. The four representative motifs contribute to the general function of CYPs. Fungal CYP51s and CYP61s can be used as the models for the substrate recognition sites analysis. The CYP proteins are clustered into 15 clades and the phylogenetic analyses suggest that the wide variety of fungal CYPs has mainly arisen from gene duplication. Two large duplication events might have been associated with the booming of Ascomycota and Basidiomycota. In addition, horizontal gene transfer also contributes to the diversification of fungal CYPs. Finally, a possible evolutionary scenario for fungal CYPs along with fungal divergences is proposed. Our results provide the fundamental information for a better understanding of CYP distribution, structure and function, and new insights into the evolutionary events of fungal CYPs along with the evolution of fungi.

microRNA in teleost fish

Abstract: MicroRNAs (miRNAs) are transcriptional and posttranscriptional regulators involved in nearly all known biological processes in distant eukaryotic clades. Their discovery and functional characterization have broadened our understanding of biological regulatory mechanisms in animals and plants. They show both evolutionary conserved and unique features across Metazoa. Here, we present the current status of the knowledge about the role of miRNA in development, growth, and physiology of teleost fishes, in comparison to other vertebrates. Infraclass Teleostei is the most abundant group among vertebrate lineage. Fish are an important component of aquatic ecosystems and human life, being the prolific source of animal proteins worldwide and a vertebrate model for biomedical research. We review miRNA biogenesis, regulation, modifications, and mechanisms of action. Specific sections are devoted to the role of miRNA in teleost development, organogenesis, tissue differentiation, growth, regeneration, reproduction, endocrine system, and responses to environmental stimuli. Each section discusses gaps in the current knowledge and pinpoints the future directions of research on miRNA in teleosts.

Genome Evolution and Plasticity of Serratia marcescens, an Important Multidrug-Resistant Nosocomial Pathogen

Abstract: Serratia marcescens is an important nosocomial pathogen that can cause an array of infections, most notably of the urinary tract and bloodstream. Naturally, it is found in many environmental niches, and is capable of infecting plants and animals. The emergence and spread of multidrug-resistant strains producing extended-spectrum or metallo beta-lactamases now pose a threat to public health worldwide. Here we report the complete genome sequences of two carefully selected S. marcescens strains, a multidrug-resistant clinical isolate (strain SM39) and an insect isolate (strain Db11). Our comparative analyses reveal the core genome of S. marcescens and define the potential metabolic capacity, virulence, and multidrug resistance of this species. We show a remarkable intraspecies genetic diversity, both at the sequence level and with regards genome flexibility, which may reflect the diversity of niches inhabited by members of this species. A broader analysis with other Serratia species identifies a set of approximately 3,000 genes that characterize the genus. Within this apparent genetic diversity, we identified many genes implicated in the high virulence potential and antibiotic resistance of SM39, including the metallo beta-lactamase and multiple other drug resistance determinants carried on plasmid pSMC1. We further show that pSMC1 is most closely related to plasmids circulating in Pseudomonas species. Our data will provide a valuable basis for future studies on S. marcescens and new insights into the genetic mechanisms that underlie the emergence of pathogens highly resistant to multiple antimicrobial agents.

Evolutionary History of the Clostridium difficile Pathogenicity Locus

Abstract: The symptoms of Clostridium difficile infection are caused by toxins expressed from its 19 kb pathogenicity locus (PaLoc). Stable integration of the PaLoc is suggested by its single chromosomal location and the clade specificity of its different genetic variants. However, the PaLoc is variably present, even among closely related strains, and thus resembles a mobile genetic element. Our aim was to explain these apparently conflicting observations by reconstructing the evolutionary history of the PaLoc. Phylogenetic analyses and annotation of the regions spanning the PaLoc were performed using C. difficile population-representative genomes chosen from a collection of 1,693 toxigenic (PaLoc present) and nontoxigenic (PaLoc absent) isolates. Comparison of the core genome and PaLoc phylogenies demonstrated an eventful evolutionary history, with distinct PaLoc variants acquired clade specifically after divergence. In particular, our data suggest a relatively recent PaLoc acquisition in clade 4. Exchanges and losses of the PaLoc DNA have also occurred, via long homologous recombination events involving flanking chromosomal sequences. The most recent loss event occurred ∼30 years ago within a clade 1 genotype. The genetic organization of the clade 3 PaLoc was unique in containing a stably integrated novel transposon (designated Tn6218), variants of which were found at multiple chromosomal locations. Tn6218 elements were Tn916-related but nonconjugative and occasionally contained genes conferring resistance to clinically relevant antibiotics. The evolutionary histories of two contrasting but clinically important genetic elements were thus characterized: the PaLoc, mobilized rarely via homologous recombination, and Tn6218, mobilized frequently through transposition.

The Evolution of the GPCR Signaling System in Eukaryotes: Modularity, Conservation, and the Transition to Metazoan Multicellularity

Abstract: The G-protein-coupled receptor (GPCR) signaling system is one of the main signaling pathways in eukaryotes. Here, we analyze the evolutionary history of all its components, from receptors to regulators, to gain a broad picture of its system-level evolution. Using eukaryotic genomes covering most lineages sampled to date, we find that the various components of the GPCR signaling pathway evolved independently, highlighting the modular nature of this system. Our data show that some GPCR families, G proteins, and regulators of G proteins diversified through lineage-specific diversifications and recurrent domain shuffling. Moreover, most of the gene families involved in the GPCR signaling system were already present in the last common ancestor of eukaryotes. Furthermore, we show that the unicellular ancestor of Metazoa already had most of the cytoplasmic components of the GPCR signaling system, including, remarkably, all the G protein alpha subunits, which are typical of metazoans. Thus, we show how the transition to multicellularity involved conservation of the signaling transduction machinery, as well as a burst of receptor diversification to cope with the new multicellular necessities.

Dissecting vancomycin-intermediate resistance in staphylococcus aureus using genome-wide association.

Abstract: Vancomycin-intermediate Staphylococcus aureus (VISA) is currently defined as having minimal inhibitory concentration (MIC) of 4–8 µg/ml. VISA evolves through changes in multiple genetic loci with at least 16 candidate genes identified in clinical and in vitro-selected VISA strains. We report a whole-genome comparative analysis of 49 vancomycin-sensitive S. aureus and 26 VISA strains. Resistance to vancomycin was determined by broth microdilution, Etest, and population analysis profile-area under the curve (PAP-AUC). Genome-wide association studies (GWAS) of 55,977 single-nucleotide polymorphisms identified in one or more strains found one highly significant association (P = 8.78E-08) between a nonsynonymous mutation at codon 481 (H481) of the rpoB gene and increased vancomycin MIC. Additionally, we used a database of public S. aureus genome sequences to identify rare mutations in candidate genes associated with VISA. On the basis of these data, we proposed a preliminary model called ECM+RMCG for the VISA phenotype as a benchmark for future efforts. The model predicted VISA based on the presence of a rare mutation in a set of candidate genes (walKR, vraSR, graSR, and agrA) and/or three previously experimentally verified mutations (including the rpoB H481 locus) with an accuracy of 81% and a sensitivity of 73%. Further, the level of resistance measured by both Etest and PAP-AUC regressed positively with the number of mutations present in a strain. This study demonstrated 1) the power of GWAS for identifying common genetic variants associated with antibiotic resistance in bacteria and 2) that rare mutations in candidate gene, identified using large genomic data sets, can also be associated with resistance phenotypes.

Phylogenomics and the dynamic genome evolution of the genus Streptococcus.

Abstract: The genus Streptococcus comprises important pathogens that have a severe impact on human health and are responsible for substantial economic losses to agriculture. Here, we utilize 46 Streptococcus genome sequences (44 species), including eight species sequenced here, to provide the first genomic level insight into the evolutionary history and genetic basis underlying the functional diversity of all major groups of this genus. Gene gain/loss analysis revealed a dynamic pattern of genome evolution characterized by an initial period of gene gain followed by a period of loss, as the major groups within the genus diversified. This was followed by a period of genome expansion associated with the origins of the present extant species. The pattern is concordant with an emerging view that genomes evolve through a dynamic process of expansion and streamlining. A large proportion of the pan-genome has experienced lateral gene transfer (LGT) with causative factors, such as relatedness and shared environment, operating over different evolutionary scales. Multiple gene ontology terms were significantly enriched for each group, and mapping terms onto the phylogeny showed that those corresponding to genes born on branches leading to the major groups represented approximately one-fifth of those enriched. Furthermore, despite the extensive LGT, several biochemical characteristics have been retained since group formation, suggesting genomic cohesiveness through time, and that these characteristics may be fundamental to each group. For example, proteolysis: mitis group; urea metabolism: salivarius group; carbohydrate metabolism: pyogenic group; and transcription regulation: bovis group.

Restriction and Recruitment—Gene Duplication and the Origin and Evolution of Snake Venom Toxins

Abstract: Snake venom has been hypothesized to have originated and diversified through a process that involves duplication of genes encoding body proteins with subsequent recruitment of the copy to the venom gland, where natural selection acts to develop or increase toxicity. However, gene duplication is known to be a rare event in vertebrate genomes, and the recruitment of duplicated genes to a novel expression domain (neofunctionalization) is an even rarer process that requires the evolution of novel combinations of transcription factor binding sites in upstream regulatory regions. Therefore, although this hypothesis concerning the evolution of snake venom is very unlikely and should be regarded with caution, it is nonetheless often assumed to be established fact, hindering research into the true origins of snake venom toxins. To critically evaluate this hypothesis, we have generated transcriptomic data for body tissues and salivary and venom glands from five species of venomous and nonvenomous reptiles. Our comparative transcriptomic analysis of these data reveals that snake venom does not evolve through the hypothesized process of duplication and recruitment of genes encoding body proteins. Indeed, our results show that many proposed venom toxins are in fact expressed in a wide variety of body tissues, including the salivary gland of nonvenomous reptiles and that these genes have therefore been restricted to the venom gland following duplication, not recruited. Thus, snake venom evolves through the duplication and subfunctionalization of genes encoding existing salivary proteins. These results highlight the danger of the elegant and intuitive “just-so story” in evolutionary biology.

Evolution and Classification of Myosins, a Paneukaryotic Whole-Genome Approach

Abstract: Myosins are key components of the eukaryotic cytoskeleton, providing motility for a broad diversity of cargoes. Therefore, understanding the origin and evolutionary history of myosin classes is crucial to address the evolution of eukaryote cell biology. Here, we revise the classification of myosins using an updated taxon sampling that includes newly or recently sequenced genomes and transcriptomes from key taxa. We performed a survey of eukaryotic genomes and phylogenetic analyses of the myosin gene family, reconstructing the myosin toolkit at different key nodes in the eukaryotic tree of life. We also identified the phylogenetic distribution of myosin diversity in terms of number of genes, associated protein domains and number of classes in each taxa. Our analyses show that new classes (i.e., paralogs) and domain architectures were continuously generated throughout eukaryote evolution, with a significant expansion of myosin abundance and domain architectural diversity at the stem of Holozoa, predating the origin of animal multicellularity. Indeed, single-celled holozoans have the most complex myosin complement among eukaryotes, with paralogs of most myosins previously considered animal specific. We recover a dynamic evolutionary history, with several lineagespecific expansions (e.g., the myosin III-like gene family diversification in choanoflagellates), convergence in protein domain architectures (e.g., fungal and animal chitin synthase myosins), and important secondary losses. Overall, our evolutionary scheme demonstrates that the ancestral eukaryote likely had a complex myosin repertoire that included six genes with different protein domain architectures. Finally, we provide an integrative and robust classification, useful for future genomic and functional studies on this crucial eukaryotic gene family.

Extensive Horizontal Gene Transfer during Staphylococcus aureus Co-colonization In Vivo

Abstract: Staphylococcus aureus is a commensal and major pathogen of humans and animals. Comparative genomics of S. aureus populations suggests that colonization of different host species is associated with carriage of mobile genetic elements (MGE), particularly bacteriophages and plasmids capable of encoding virulence, resistance, and immune evasion pathways. Antimicrobial-resistant S. aureus of livestock are a potential zoonotic threat to human health if they adapt to colonize humans efficiently. We utilized the technique of experimental evolution and co-colonized gnotobiotic piglets with both human- and pig-associated variants of the lineage clonal complex 398, and investigated growth and genetic changes over 16 days using whole genome sequencing. The human isolate survived co-colonization on piglets more efficiently than in vitro. During co-colonization, transfer of MGE from the pig to the human isolate was detected within 4 h. Extensive and repeated transfer of two bacteriophages and three plasmids resulted in colonization with isolates carrying a wide variety of mobilomes. Whole genome sequencing of progeny bacteria revealed no acquisition of core genome polymorphisms, highlighting the importance of MGE. Staphylococcus aureus bacteriophage recombination and integration into novel sites was detected experimentally for the first time. During colonization, clones coexisted and diversified rather than a single variant dominating. Unexpectedly, each piglet carried unique populations of bacterial variants, suggesting limited transmission of bacteria between piglets once colonized. Our data show that horizontal gene transfer occurs at very high frequency in vivo and significantly higher than that detectable in vitro.

Comparative Genomic Analysis and Virulence Differences in Closely Related Salmonella enterica Serotype Heidelberg Isolates from Humans, Retail Meats, and Animals

Abstract: Salmonella enterica subsp. enterica serovar Heidelberg (S. Heidelberg) is one of the top serovars causing human salmonellosis. Recently, an antibiotic-resistant strain of this serovar was implicated in a large 2011 multistate outbreak resulting from consumption of contaminated ground turkey that involved 136 confirmed cases, with one death. In this study, we assessed the evolutionary diversity of 44 S. Heidelberg isolates using whole-genome sequencing (WGS) generated by the 454 GS FLX (Roche) platform. The isolates, including 30 with nearly indistinguishable (one band difference) Xbal pulsed-field gel electrophoresis patterns (JF6X01.0032, JF6X01.0058), were collected from various sources between 1982 and 2011 and included nine isolates associated with the 2011 outbreak. Additionally, we determined the complete sequence for the chromosome and three plasmids from a clinical isolate associated with the 2011 outbreak using the Pacific Biosciences (PacBio) system. Using single-nucleotide polymorphism (SNP) analyses, we were able to distinguish highly clonal isolates, including strains isolated at different times in the same year. The isolates from the recent 2011 outbreak clustered together with a mean SNP variation of only 17 SNPs. The S. Heidelberg isolates carried a variety of phages, such as prophage P22, P4, lambda-like prophage Gifsy-2, and the P2-like phage which carries the sopE1 gene, virulence genes including 62 pathogenicity, and 13 fimbrial markers and resistance plasmids of the incompatibility (Inc)I1, IncA/C, and IncHI2 groups. Twenty-one strains contained an IncX plasmid carrying a type IV secretion system. On the basis of the recent and historical isolates used in this study, our results demonstrated that, in addition to providing detailed genetic information for the isolates, WGS can identify SNP targets that can be utilized for differentiating highly clonal S. Heidelberg isolates.

Six Subgroups and Extensive Recent Duplications Characterize the Evolution of the Eukaryotic Tubulin Protein Family

Abstract: Tubulins belong to the most abundant proteins in eukaryotes providing the backbone for many cellular substructures like the mitotic and meiotic spindles, the intracellular cytoskeletal network, and the axonemes of cilia and flagella. Homologs have even been reported for archaea and bacteria. However, a taxonomically broad and whole-genome-based analysis of the tubulin protein family has never been performed, and thus, the number of subfamilies, their taxonomic distribution, and the exact grouping of the supposed archaeal and bacterial homologs are unknown. Here, we present the analysis of 3,524 tubulins from 504 species. The tubulins formed six major subfamilies, a to z. Species of all major kingdoms of the eukaryotes encode members of these subfamilies implying that they must have already been present in the last common eukaryotic ancestor. The proposed archaeal homologs grouped together with the bacterial TubZ proteins as sister clade to the FtsZ proteins indicating that tubulins are unique to eukaryotes. Most species containeda- and/orb-tubulin gene duplicates resulting from recent branch- and species-specific duplication events. This shows that tubulins cannot be used for constructing species phylogenies without resolving their ortholog–paralog relationships. The many gene duplicates and also the independent loss of the d-, "-, or z-tubulins, which have been shown to be part of the triplet microtubules in basal bodies, suggest that tubulins can functionally substitute each other.

Genetic convergence in the adaptation of dogs and humans to the high altitude environment of the Tibetan plateau

Abstract: The high-altitude hypoxic environment represents one of the most extreme challenges for mammals. Previous studies of humans on the Tibetan plateau and in the Andes Mountains have identified statistical signatures of selection in different sets of loci. Here, we first measured the hemoglobin levels in village dogs from Tibet and those from Chinese lowlands. We found that the hemoglobin levels are very similar between the two groups, suggesting that Tibetan dogs might share similar adaptive strategies as the Tibetan people. Through a whole-genome sequencing approach, we have identified EPAS1 and HBB as candidate genes for the hypoxic adaptation on the Tibetan plateau. The population genetic analysis shows a significant convergence between humans and dogs in Tibet. The similarities in the sets of loci that exhibit putative signatures of selection and the hemoglobin levels between humans and dogs of the same environment, but not between human populations in different regions, suggests an extraordinary landscape of convergent evolution between human beings and their best friend on the Tibetan plateau.

Epigenetics and the Evolution of Darwin's Finches

Abstract: The prevailing theory for the molecular basis of evolution involves genetic mutations that ultimately generate the heritable phenotypic variation on which natural selection acts. However, epigenetic transgenerational inheritance of phenotypic variation may also play an important role in evolutionary change. A growing number of studies have demonstrated the presence of epigenetic inheritance in a variety of different organisms that can persist for hundreds of generations. The possibility that epigenetic changes can accumulate over longer periods of evolutionary time has seldom been tested empirically. This study was designed to compare epigenetic changes among several closely related species of Darwin's finches, a well-known example of adaptive radiation. Erythrocyte DNA was obtained from five species of sympatric Darwin's finches that vary in phylogenetic relatedness. Genome-wide alterations in genetic mutations using copy number variation (CNV) were compared with epigenetic alterations associated with differential DNA methylation regions (epimutations). Epimutations were more common than genetic CNV mutations among the five species; furthermore, the number of epimutations increased monotonically with phylogenetic distance. Interestingly, the number of genetic CNV mutations did not consistently increase with phylogenetic distance. The number, chromosomal locations, regional clustering, and lack of overlap of epimutations and genetic mutations suggest that epigenetic changes are distinct and that they correlate with the evolutionary history of Darwin's finches. The potential functional significance of the epimutations was explored by comparing their locations on the genome to the location of evolutionarily important genes and cellular pathways in birds. Specific epimutations were associated with genes related to the bone morphogenic protein, toll receptor, and melanogenesis signaling pathways. Species-specific epimutations were significantly overrepresented in these pathways. As environmental factors are known to result in heritable changes in the epigenome, it is possible that epigenetic changes contribute to the molecular basis of the evolution of Darwin's finches.

Evidence for Soft Selective Sweeps in the Evolution of Pneumococcal Multidrug Resistance and Vaccine Escape

Abstract: The multidrug-resistant Streptococcus pneumoniae Taiwan 19F -14, or PMEN14, clone was first observed with a 19F serotype, which is targeted by the heptavalent polysaccharide conjugate vaccine (PCV7). However, “vaccine escape” PMEN14 isolates with a 19A serotype became an increasingly important cause of disease post-PCV7. Whole genome sequencing was used to characterize the recent evolution of 173 pneumococci of, or related to, PMEN14. This suggested that PMEN14 is a single lineage that originated in the late 1980s in parallel with the acquisition of multiple resistances by close relatives. One of the four detected serotype switches to 19A generated representatives of the sequence type (ST) 320 isolates that have been highly successful post-PCV7. A second produced an ST236 19A genotype with reduced resistance to b-lactams owing to alteration of pbp1a and pbp2x sequences through the same recombination that caused the change in serotype. A third, which generated a mosaic capsule biosynthesis locus, resulted in serotype 19A ST271 isolates. The rapid diversification through homologous recombination seen in the global collection was similarly observed in the absence of vaccination in a set of isolates from the Maela refugee camp in Thailand, a collection that also allowed variation to be observed within carriage through longitudinal sampling. This suggests that some pneumococcal genotypes generate a pool of standing variation that is sufficiently extensive to result in “soft” selective sweeps: The emergence of multiple mutants in parallel

Population-genetic inference from pooled-sequencing data.

Abstract: Although pooled-population sequencing has become a widely used approach for estimating allele frequencies, most work has proceeded in the absence of a proper statistical framework. We introduce a self-sufficient, closed-form, maximum-likelihood estimator for allele frequencies that accounts for errors associated with sequencing, and a likelihood-ratio test statistic that provides a simple means for evaluating the null hypothesis of monomorphism. Unbiased estimates of allele frequencies (where N is the number of individuals sampled) appear to be unachievable, and near-certain identification of a polymorphism requires a minor-allele frequency . A framework is provided for testing for significant differences in allele frequencies between populations, taking into account sampling at the levels of individuals within populations and sequences within pooled samples. Analyses that fail to account for the two tiers of sampling suffer from very large false-positive rates and can become increasingly misleading with increasing depths of sequence coverage. The power to detect significant allele-frequency differences between two populations is very limited unless both the number of sampled individuals and depth of sequencing coverage exceed 100.

Massive Losses of Taste Receptor Genes in Toothed and Baleen Whales

Abstract: Taste receptor genes are functionally important in animals, with a surprising exception in the bottlenose dolphin, which shows extensive losses of sweet, umami, and bitter taste receptor genes. To examine the generality of taste gene loss, we examined seven toothed whales and five baleen whales and sequenced the complete repertoire of three sweet/umami (T1Rs) and ten bitter (T2Rs) taste receptor genes. We found all amplified T1Rs and T2Rs to be pseudogenes in all 12 whales, with a shared premature stop codon in 10 of the 13 genes, which demonstrated massive losses of taste receptor genes in the common ancestor of whales. Furthermore, we analyzed three genome sequences from two toothed whales and one baleen whale and found that the sour taste marker gene Pkd2l1 is a pseudogene, whereas the candidate salty taste receptor genes are intact and putatively functional. Additionally, we examined three genes that are responsible for taste signal transduction and found the relaxation of functional constraints on taste signaling pathways along the ancestral branch leading to whales. Together, our results strongly suggest extensive losses of sweet, umami, bitter, and sour tastes in whales, and the relaxation of taste function most likely arose in the common ancestor of whales between 36 and 53 Ma. Therefore, whales represent the first animal group to lack four of five primary tastes, probably driven by the marine environment with high concentration of sodium, the feeding behavior of swallowing prey whole, and the dietary switch from plants to meat in the whale ancestor.

Gene Loss Rather Than Gene Gain Is Associated with a Host Jump from Monocots to Dicots in the Smut Fungus Melanopsichium pennsylvanicum

Abstract: Smut fungi are well-suited to investigate the ecology and evolution of plant pathogens, as they are strictly biotrophic, yet cultivable on media. Here we report the genome sequence of Melanopsichium pennsylvanicum, closely related to Ustilago maydis and other Poaceae-infecting smuts, but parasitic to a dicot plant. To explore the evolutionary patterns resulting from host adaptation after this huge host jump, the genome of Me. pennsylvanicum was sequenced and compared with the genomes of U. maydis, Sporisorium reilianum ,a ndU. hordei. Although all four genomes had a similar completeness in CEGMA (Core Eukaryotic Genes Mapping Approach) analysis, gene absence was highest in Me. pennsylvanicum, and most pronounced in putative secreted proteins, which are often considered as effector candidates. In contrast, the amount of private genes was similar among the species, highlighting that gene loss rather than gene gain is the hallmark of adaptation after the host jump to the dicot host. Our analyses revealed a trend of putative effectors to be next to another putative effector, but the majority of these are not in clusters and thus the focus on pathogenicity clusters might not be appropriate for all smut genomes. Positive selection studies revealed that Me. pennsylvanicum has the highest number and proportion of genes under positive selection. In general, putative effectors showed a higher proportion of positively selected genes than noneffector candidates. The 248 putative secreted effectors found in all four smut genomes might constitute a core set needed for pathogenicity, whereas those 92 that are found in all grass-parasitic smuts but have no ortholog in Me. pennsylvanicum might constitute a set of effectors important for successful colonization of grass hosts.

Chromera velia, endosymbioses and the rhodoplex hypothesis--plastid evolution in cryptophytes, alveolates, stramenopiles, and haptophytes (CASH lineages).

Abstract: The discovery of Chromera velia, a free-living photosynthetic relative of apicomplexan pathogens, has provided an unexpected opportunity to study the algal ancestry of malaria parasites. In this work, we compared the molecular footprints of a eukaryote-to-eukaryote endosymbiosis in C. velia to their equivalents in peridinin-containing dinoflagellates (PCD) to reevaluate recent claims in favor of a common ancestry of their plastids. To this end, we established the draft genome and a set of full-length cDNA sequences from C. velia via next-generation sequencing. We documented the presence of a single coxI gene in the mitochondrial genome, which thus represents the genetically most reduced aerobic organelle identified so far, but focused our analyses on five “lucky genes” of the Calvin cycle. These were selected because of their known support for a common origin of complex plastids from cryptophytes, alveolates (represented by PCDs), stramenopiles, and haptophytes (CASH) via a single secondary endosymbiosis with a red alga. As expected, our broadly sampled phylogenies of the nuclear-encoded Calvin cycle markers support a rhodophycean origin for the complex plastid of Chromera. However, they also suggest an independent origin of apicomplexan and dinophycean (PCD) plastids via two eukaryote-to-eukaryote endosymbioses. Although at odds with the current view of a common photosynthetic ancestry for alveolates, this conclusion is nonetheless in line with the deviant plastome architecture in dinoflagellates and the morphological paradox of four versus three plastid membranes in the respective lineages. Further support for independent endosymbioses is provided by analysis of five additional markers, four of them involved in the plastid protein import machinery. Finally, we introduce the “rhodoplex hypothesis” as a convenient way to designate evolutionary scenarios where CASH plastids are ultimately the product of a single secondary endosymbiosis with a red alga but were subsequently horizontally spread via higher-order eukaryote-to-eukaryote endosymbioses.

Genome Comparison of Candida orthopsilosis Clinical Strains Reveals the Existence of Hybrids between Two Distinct Subspecies

Abstract: The Candida parapsilosis species complex comprises a group of emerging human pathogens of varying virulence. This complex was recently subdivided into three different species: C. parapsilosis sensu stricto, C. metapsilosis, and C. orthopsilosis. Within the latter, at least two clearly distinct subspecies seem to be present among clinical isolates (Type 1 and Type 2). To gain insight into the genomic differences between these subspecies, we undertook the sequencing of a clinical isolate classified as Type 1 and compared it with the available sequence of a Type 2 clinical strain. Unexpectedly, the analysis of the newly sequenced strain revealed a highly heterozygous genome, which we show to be the consequence of a hybridization event between both identified subspecies. This implicitly suggests that C. orthopsilosis is able to mate, a so-far unanswered question. The resulting hybrid shows a chimeric genome that maintains a similar gene dosage from both parental lineages and displays ongoing loss of heterozygosity. Several of the differences found between the gene content in both strains relate to virulent-related families, with the hybrid strain presenting a higher copy number of genes coding for efflux pumps or secreted lipases. Remarkably, two clinical strains isolated from distant geographical locations (Texas and Singapore) are descendants of the same hybrid line, raising the intriguing possibility of a relationship between the hybridization event and the global spread of a virulent clone.

Expression Profiling of the Wheat Pathogen Zymoseptoria tritici Reveals Genomic Patterns of Transcription and Host-Specific Regulatory Programs

Abstract: Host specialization by pathogens requires a repertoire of virulence factors as well as fine-tuned regulation of gene expression. The fungal wheat pathogen Zymoseptoria tritici (synonym Mycosphaerella graminicola) is a powerful model system for the discovery of genetic elements that underlie virulence and host specialization. We transcriptionally profiled the early stages of Z. tritici infection of a compatible host (wheat) and a noncompatible host (Brachypodium distachyon). The results revealed infection regulatory programs common to both hosts and genes with striking wheat-specific expression, with many of the latter showing sequence signatures of positive selection along the Z. tritici lineage. Genes specifically regulated during infection of wheat populated two large clusters of coregulated genes that may represent candidate pathogenicity islands. On evolutionarily labile, repeat-rich accessory chromosomes (ACs), we identified hundreds of highly expressed genes with signatures of evolutionary constraint and putative biological function. Phylogenetic analyses suggested that gene duplication events on these ACs were rare and largely preceded the diversification of Zymoseptoria species. Together, our data highlight the likely relevance for fungal growth and virulence of hundreds of Z. tritici genes, deepening the annotation and functional inference of the genes of this model pathogen.

Pangenome Evidence for Extensive Interdomain Horizontal Transfer Affecting Lineage Core and Shell Genes in Uncultured Planktonic Thaumarchaeota and Euryarchaeota

Abstract: Horizontal gene transfer (HGT) is an important force in evolution, which may lead, among other things, to the adaptation to new environments by the import of new metabolic functions. Recent studies based on phylogenetic analyses of a few genome fragments containing archaeal 16S rRNA genes and fosmid-end sequences from deep-sea metagenomic libraries have suggested that marine planktonic archaea could be affected by high HGT frequency. Likewise, a composite genome of an uncultured marine euryarchaeote showed high levels of gene sequence similarity to bacterial genes. In this work, we ask whether HGT is frequent and widespread in genomes of these marine archaea, and whether HGT is an ancient and/or recurrent phenomenon. To answer these questions, we sequenced 997 fosmid archaeal clones from metagenomic libraries of deep-Mediterranean waters (1,000 and 3,000 m depth) and built comprehensive pangenomes for planktonic Thaumarchaeota (Group I archaea) and Euryarchaeota belonging to the uncultured Groups II and III Euryarchaeota (GII/III-Euryarchaeota). Comparison with available reference genomes of Thaumarchaeota and a composite marine surface euryarchaeote genome allowed us to define sets of core, lineage-specific core, and shell gene ortholog clusters for the two archaeal lineages. Molecular phylogenetic analyses of all gene clusters showed that 23.9% of marine Thaumarchaeota genes and 29.7% of GII/III-Euryarchaeota genes had been horizontally acquired from bacteria. HGT is not only extensive and directional but also ongoing, with high HGT levels in lineage-specific core (ancient transfers) and shell (recent transfers) genes. Many of the acquired genes are related to metabolism and membrane biogenesis, suggesting an adaptive value for life in cold, oligotrophic oceans. We hypothesize that the acquisition of an important amount of foreign genes by the ancestors of these archaeal groups significantly contributed to their divergence and ecological success.

Archaeal “Dark Matter” and the Origin of Eukaryotes

Abstract: Current hypotheses about the history of cellular life are mainly based on analyses of cultivated organisms, but these represent only a small fraction of extant biodiversity. The sequencing of new environmental lineages therefore provides an opportunity to test, revise, or reject existing ideas about the tree of life and the origin of eukaryotes. According to the textbook three domains hypothesis, the eukaryotes emerge as the sister group to a monophyletic Archaea. However, recent analyses incorporating better phylogenetic models and an improved sampling of the archaeal domain have generally supported the competing eocyte hypothesis, in which core genes of eukaryotic cells originated from within the Archaea, with important implications for eukaryogenesis. Given this trend, it was surprising that a recent analysis incorporating new genomes from uncultivated Archaea recovered a strongly supported three domains tree. Here, we show that this result was due in part to the use of a poorly fitting phylogenetic model and also to the inclusion by an automated pipeline of genes of putative bacterial origin rather than nucleocytosolic versions for some of the eukaryotes analyzed. When these issues were resolved, analyses including the new archaeal lineages placed core eukaryotic genes within the Archaea. These results are consistent with a number of recent studies in which improved archaeal sampling and better phylogenetic models agree in supporting the eocyte tree over the three domains hypothesis.

Settling down: the genome of Serratia symbiotica from the aphid Cinara tujafilina zooms in on the process of accommodation to a cooperative intracellular life.

Abstract: Particularly interesting cases of mutualistic endosymbioses come from the establishment of co-obligate associations of more than one species of endosymbiotic bacteria. Throughout symbiotic accommodation from a free-living bacterium, passing through a facultative stage and ending as an obligate intracellular one, the symbiont experiences massive genomic losses and phenotypic adjustments. Here, we scrutinized the changes in the coevolution of Serratia symbiotica and Buchnera aphidicola endosymbionts in aphids, paying particular attention to the transformations undergone by S. symbiotica to become an obligate endosymbiont. Although it is already known that S. symbiotica is facultative in Acyrthosiphon pisum, in Cinara cedri it has established a co-obligate endosymbiotic consortium along with B. aphidicola to fulfill the aphid’s nutritional requirements. The state of this association in C. tujafilina, an aphid belonging to the same subfamily (Lachninae) that C. cedri, remained unknown. Here, we report the genome of S. symbiotica strain SCt-VLC from the aphid C. tujafilina. While being phylogenetically and genomically very closely related to the facultative endosymbiont S. symbiotica from the aphid A. pisum, it shows a variety of metabolic, genetic, and architectural features, which point toward this endosymbiont being one step closer to an obligate intracellular one. We also describe in depth the process of genome rearrangements suffered by S. symbiotica and the role mobile elements play in gene inactivations. Finally, we postulate the supply to the host of the essential riboflavin (vitamin B2) as key to the establishment of S. symbiotica as a co-obligate endosymbiont in the aphids belonging to the subfamily Lachninane.