Evolutionary implications of horizontal gene transfer

Gene expression, nucleotide composition and codon usage bias of genes associated with human Y chromosome

Analysis of codon usage pattern is important to understand the genetic and evolutionary characteristics of genomes. We have used bioinformatic approaches to analyze the codon usage bias (CUB) of the genes located in human Y chromosome. Codon bias index (CBI) indicated that the overall extent of codon usage bias was low. The relative synonymous codon usage (RSCU) analysis suggested that approximately half of the codons out of 59 synonymous codons were most frequently used, and possessed a T or G at the third codon position. The codon usage pattern was different in different genes as revealed from correspondence analysis (COA). A significant correlation between effective number of codons (ENC) and various GC contents suggests that both mutation pressure and natural selection affect the codon usage pattern of genes located in human Y chromosome. In addition, Y-linked genes have significant difference in GC contents at the second and third codon positions, expression level, and codon usage pattern of some codons like the SPANX genes in X chromosome.

Ecological networks to unravel the routes to horizontal transposon transfers

Transposable elements (TEs) represent the single largest component of numerous eukaryotic genomes, and their activity and dispersal constitute an important force fostering evolutionary innovation. The horizontal transfer of TEs (HTT) between eukaryotic species is a common and widespread phenomenon that has had a profound impact on TE dynamics and, consequently, on the evolutionary trajectory of many species' lineages. However, the mechanisms promoting HTT remain largely unknown. In this article, we argue that network theory combined with functional ecology provides a robust conceptual framework and tools to delineate how complex interactions between diverse organisms may act in synergy to promote HTTs.

Clusterflock: a flocking algorithm for isolating congruent phylogenomic datasets

BACKGROUND

Collective animal behavior, such as the flocking of birds or the shoaling of fish, has inspired a class of algorithms designed to optimize distance-based clusters in various applications, including document analysis and DNA microarrays. In a flocking model, individual agents respond only to their immediate environment and move according to a few simple rules. After several iterations the agents self-organize, and clusters emerge without the need for partitional seeds. In addition to its unsupervised nature, flocking offers several computational advantages, including the potential to reduce the number of required comparisons.

FINDINGS

In the tool presented here, Clusterflock, we have implemented a flocking algorithm designed to locate groups (flocks) of orthologous gene families (OGFs) that share an evolutionary history. Pairwise distances that measure phylogenetic incongruence between OGFs guide flock formation. We tested this approach on several simulated datasets by varying the number of underlying topologies, the proportion of missing data, and evolutionary rates, and show that in datasets containing high levels of missing data and rate heterogeneity, Clusterflock outperforms other well-established clustering techniques. We also verified its utility on a known, large-scale recombination event in Staphylococcus aureus. By isolating sets of OGFs with divergent phylogenetic signals, we were able to pinpoint the recombined region without forcing a pre-determined number of groupings or defining a pre-determined incongruence threshold.

CONCLUSIONS

Clusterflock is an open-source tool that can be used to discover horizontally transferred genes, recombined areas of chromosomes, and the phylogenetic 'core' of a genome. Although we used it here in an evolutionary context, it is generalizable to any clustering problem. Users can write extensions to calculate any distance metric on the unit interval, and can use these distances to 'flock' any type of data.

Dioxygenases of Chlorobiphenyl-Degrading Species Rhodococcus wratislaviensis G10 and Chlorophenol-Degrading Species Rhodococcus opacus 1CP Induced in Benzoate-Grown Cells and Genes Potentially Involved in These Processes

Dioxygenases induced during benzoate degradation by the actinobacterium Rhodococcus wratislaviensis G10 strain degrading haloaromatic compounds were studied. Rhodococcus wratislaviensis G10 completely degraded 2 g/liter benzoate during 30 h and 10 g/liter during 200 h. Washed cells grown on benzoate retained respiration activity for more than 90 days, and a high activity of benzoate dioxygenase was recorded for 10 days. Compared to the enzyme activities with benzoate, the activity of benzoate dioxygenases was 10-30% with 13 of 35 substituted benzoate analogs. Two dioxygenases capable of cleaving the aromatic ring were isolated and characterized: protocatechuate 3,4-dioxygenase and catechol 1,2-dioxygenase. Catechol inhibited the activity of protocatechuate 3,4-dioxygenase. Protocatechuate did not affect the activity of catechol 1,2-dioxygenase. A high degree of identity was shown by MALDI-TOF mass spectrometry for protein peaks of the R. wratislaviensis G10 and Rhodococcus opacus 1CP cells grown on benzoate or LB. DNA from the R. wratislaviensis G10 strain was specifically amplified using specific primers to variable regions of genes coding α- and β-subunits of protocatechuate 3,4-dioxygenase and to two genes of the R. opacus 1CP coding catechol 1,2-dioxygenase. The products were 99% identical with the corresponding regions of the R. opacus 1CP genes. This high identity (99%) between the genes coding degradation of aromatic compounds in the R. wratislaviensis G10 and R. opacus 1CP strains isolated from sites of remote location (1400 km) and at different time (20-year difference) indicates a common origin of biodegradation genes of these strains and a wide distribution of these genes among rhodococci.

Advances in biodiversity: metagenomics and the unveiling of biological dark matter

BACKGROUND

Efforts to harmonize genomic data standards used by the biodiversity and metagenomic research communities have shown that prokaryotic data cannot be understood or represented in a traditional, classical biological context for conceptual reasons, not technical ones.

RESULTS

Biology, like physics, has a fundamental duality-the classical macroscale eukaryotic realm vs. the quantum microscale microbial realm-with the two realms differing profoundly, and counter-intuitively, from one another. Just as classical physics is emergent from and cannot explain the microscale realm of quantum physics, so classical biology is emergent from and cannot explain the microscale realm of prokaryotic life. Classical biology describes the familiar, macroscale realm of multi-cellular eukaryotic organisms, which constitute a highly derived and constrained evolutionary subset of the biosphere, unrepresentative of the vast, mostly unseen, microbial world of prokaryotic life that comprises at least half of the planet's biomass and most of its genetic diversity. The two realms occupy fundamentally different mega-niches: eukaryotes interact primarily mechanically with the environment, prokaryotes primarily physiologically. Further, many foundational tenets of classical biology simply do not apply to prokaryotic biology.

CONCLUSIONS

Classical genetics one held that genes, arranged on chromosomes like beads on a string, were the fundamental units of mutation, recombination, and heredity. Then, molecular analysis showed that there were no fundamental units, no beads, no string. Similarly, classical biology asserts that individual organisms and species are fundamental units of ecology, evolution, and biodiversity, composing an evolutionary history of objectively real, lineage-defined groups in a single-rooted tree of life. Now, metagenomic tools are forcing a recognition that there are no completely objective individuals, no unique lineages, and no one true tree. The newly revealed biosphere of microbial dark matter cannot be understood merely by extending the concepts and methods of eukaryotic macrobiology. The unveiling of biological dark matter is allowing us to see, for the first time, the diversity of the entire biosphere and, to paraphrase Darwin, is providing a new view of life. Advancing and understanding that view will require major revisions to some of the most fundamental concepts and theories in biology.

Genomic Recoding Broadly Obstructs the Propagation of Horizontally Transferred Genetic Elements

Horizontally transferred genetic elements such as viruses and conjugative plasmids move DNA between organisms, increasing genetic diversity but destabilizing engineered biological systems. Here, we used a genomically recoded Escherichia coli strain lacking UAG stop codons and the recognition protein release factor 1 to study how an alternative genetic code influences horizontally transferred genetic element propagation. The alternative genetic code conferred resistance to multiple viruses (λ, M13, P1, MS2) at titers up to 10(11) PFU/ml and impaired conjugative plasmids (F and RK2) up to 10(5)-fold. By recoding UAG codons to UAA in viruses and plasmids, we restored viral infectivity and conjugative function. Propagating viruses on a mixed community of cells with standard and alternative genetic codes reduced viral titer, and over time viruses adapted to the alternative genetic code. This work demonstrates that altering the genetic code broadly obstructs the propagation of horizontally transferred genetic elements and supports the use of genomic recoding as a strategy to stabilize engineered biological systems.

Transfer of DNA from Bacteria to Eukaryotes

Historically, the members of the Agrobacterium genus have been considered the only bacterial species naturally able to transfer and integrate DNA into the genomes of their eukaryotic hosts. Yet, increasing evidence suggests that this ability to genetically transform eukaryotic host cells might be more widespread in the bacterial world. Indeed, analyses of accumulating genomic data reveal cases of horizontal gene transfer from bacteria to eukaryotes and suggest that it represents a significant force in adaptive evolution of eukaryotic species. Specifically, recent reports indicate that bacteria other than Agrobacterium, such as Bartonella henselae (a zoonotic pathogen), Rhizobium etli (a plant-symbiotic bacterium related to Agrobacterium), or even Escherichia coli, have the ability to genetically transform their host cells under laboratory conditions. This DNA transfer relies on type IV secretion systems (T4SSs), the molecular machines that transport macromolecules during conjugative plasmid transfer and also during transport of proteins and/or DNA to the eukaryotic recipient cells. In this review article, we explore the extent of possible transfer of genetic information from bacteria to eukaryotic cells as well as the evolutionary implications and potential applications of this transfer.

CDI Systems Are Stably Maintained by a Cell-Contact Mediated Surveillance Mechanism

Contact-dependent growth inhibition (CDI) systems are widespread amongst Gram-negative bacteria where they play important roles in inter-cellular competition and biofilm formation. CDI+ bacteria use cell-surface CdiA proteins to bind neighboring bacteria and deliver C-terminal toxin domains. CDI+ cells also express CdiI immunity proteins that specifically neutralize toxins delivered from adjacent siblings. Genomic analyses indicate that cdi loci are commonly found on plasmids and genomic islands, suggesting that these Type 5 secretion systems are spread through horizontal gene transfer. Here, we examine whether CDI toxin and immunity activities serve to stabilize mobile genetic elements using a minimal F plasmid that fails to partition properly during cell division. This F plasmid is lost from Escherichia coli populations within 50 cell generations, but is maintained in ~60% of the cells after 100 generations when the plasmid carries the cdi gene cluster from E. coli strain EC93. By contrast, the ccdAB "plasmid addiction" module normally found on F exerts only a modest stabilizing effect. cdi-dependent plasmid stabilization requires the BamA receptor for CdiA, suggesting that plasmid-free daughter cells are inhibited by siblings that retain the CDI+ plasmid. In support of this model, the CDI+ F plasmid is lost rapidly from cells that carry an additional cdiI immunity gene on a separate plasmid. These results indicate that plasmid stabilization occurs through elimination of non-immune cells arising in the population via plasmid loss. Thus, genetic stabilization reflects a strong selection for immunity to CDI. After long-term passage for more than 300 generations, CDI+ plasmids acquire mutations that increase copy number and result in 100% carriage in the population. Together, these results show that CDI stabilizes genetic elements through a toxin-mediated surveillance mechanism in which cells that lose the CDI system are detected and eliminated by their siblings.

Phylogeny Reconstruction with Alignment-Free Method That Corrects for Horizontal Gene Transfer

Advances in sequencing have generated a large number of complete genomes. Traditionally, phylogenetic analysis relies on alignments of orthologs, but defining orthologs and separating them from paralogs is a complex task that may not always be suited to the large datasets of the future. An alternative to traditional, alignment-based approaches are whole-genome, alignment-free methods. These methods are scalable and require minimal manual intervention. We developed SlopeTree, a new alignment-free method that estimates evolutionary distances by measuring the decay of exact substring matches as a function of match length. SlopeTree corrects for horizontal gene transfer, for composition variation and low complexity sequences, and for branch-length nonlinearity caused by multiple mutations at the same site. We tested SlopeTree on 495 bacteria, 73 archaea, and 72 strains of Escherichia coli and Shigella. We compared our trees to the NCBI taxonomy, to trees based on concatenated alignments, and to trees produced by other alignment-free methods. The results were consistent with current knowledge about prokaryotic evolution. We assessed differences in tree topology over different methods and settings and found that the majority of bacteria and archaea have a core set of proteins that evolves by descent. In trees built from complete genomes rather than sets of core genes, we observed some grouping by phenotype rather than phylogeny, for instance with a cluster of sulfur-reducing thermophilic bacteria coming together irrespective of their phyla. The source-code for SlopeTree is available at: http://prodata.swmed.edu/download/pub/slopetree_v1/slopetree.tar.gz.

Due to their lack of distinct morphological features, bacteria and archaea were extremely difficult to classify until technology was developed to obtain their DNA sequences; these sequences could then be compared to estimate evolutionary relationships. Now, due to technological advances, there is a flood of available sequences from a wide variety of organisms. These advances have spurred the development of algorithms which can estimate evolutionary relationships using whole genomes, in contrast to the more traditional methods which used single genes earlier and now typically use groups of conserved genes. However, there are many challenges when attempting to infer evolutionary relationships, in particular horizontal gene transfer, where DNA is transferred from one organism to another, resulting in an organism’s genome containing DNA that does not reflect its evolution by descent. We developed a new whole-genome method for estimating evolutionary distances which identifies and corrects for horizontal transfer. We found that for SlopeTree and all other whole-genome methods we applied, horizontal transfer causes some evolutionary distances to be grossly underestimated, and that our correction corrects for this.

A functional difference between native and horizontally acquired genes in bdelloid rotifers

The form of RNA processing known as SL trans-splicing involves the transfer of a short conserved sequence, the spliced leader (SL), from a noncoding SL RNA to the 5' ends of mRNA molecules. SL trans-splicing occurs in several animal taxa, including bdelloid rotifers (Rotifera, Bdelloidea). One striking feature of these aquatic microinvertebrates is the large proportion of foreign genes, i.e. those acquired by horizontal gene transfer from other organisms, in their genomes. However, whether such foreign genes behave similarly to native genes has not been tested in bdelloids or any other animal. We therefore used a combination of experimental and computational methods to examine whether transcripts of foreign genes in bdelloids were SL trans-spliced, like their native counterparts. We found that many foreign transcripts contain SLs, use similar splice acceptor sequences to native genes, and are able to undergo alternative trans-splicing. However, a significantly lower proportion of foreign mRNAs contains SL sequences than native transcripts. This demonstrates a novel functional difference between foreign and native genes in bdelloids and suggests that SL trans-splicing is not essential for the expression of foreign genes, but is acquired during their domestication.

Impact of genomics on the understanding of microbial evolution and classification: the importance of Darwin's views on classification

Analyses of genome sequences, by some approaches, suggest that the widespread occurrence of horizontal gene transfers (HGTs) in prokaryotes disguises their evolutionary relationships and have led to questioning of the Darwinian model of evolution for prokaryotes. These inferences are critically examined in the light of comparative genome analysis, characteristic synapomorphies, phylogenetic trees and Darwin's views on examining evolutionary relationships. Genome sequences are enabling discovery of numerous molecular markers (synapomorphies) such as conserved signature indels (CSIs) and conserved signature proteins (CSPs), which are distinctive characteristics of different prokaryotic taxa. Based on these molecular markers, exhibiting high degree of specificity and predictive ability, numerous prokaryotic taxa of different ranks, currently identified based on the 16S rRNA gene trees, can now be reliably demarcated in molecular terms. Within all studied groups, multiple CSIs and CSPs have been identified for successive nested clades providing reliable information regarding their hierarchical relationships and these inferences are not affected by HGTs. These results strongly support Darwin's views on evolution and classification and supplement the current phylogenetic framework based on 16S rRNA in important respects. The identified molecular markers provide important means for developing novel diagnostics, therapeutics and for functional studies providing important insights regarding prokaryotic taxa.

Nothing in Evolution Makes Sense Except in the Light of Genomics: Read-Write Genome Evolution as an Active Biological Process

The 21st century genomics-based analysis of evolutionary variation reveals a number of novel features impossible to predict when Dobzhansky and other evolutionary biologists formulated the neo-Darwinian Modern Synthesis in the middle of the last century. These include three distinct realms of cell evolution; symbiogenetic fusions forming eukaryotic cells with multiple genome compartments; horizontal organelle, virus and DNA transfers; functional organization of proteins as systems of interacting domains subject to rapid evolution by exon shuffling and exonization; distributed genome networks integrated by mobile repetitive regulatory signals; and regulation of multicellular development by non-coding lncRNAs containing repetitive sequence components. Rather than single gene traits, all phenotypes involve coordinated activity by multiple interacting cell molecules. Genomes contain abundant and functional repetitive components in addition to the unique coding sequences envisaged in the early days of molecular biology. Combinatorial coding, plus the biochemical abilities cells possess to rearrange DNA molecules, constitute a powerful toolbox for adaptive genome rewriting. That is, cells possess "Read-Write Genomes" they alter by numerous biochemical processes capable of rapidly restructuring cellular DNA molecules. Rather than viewing genome evolution as a series of accidental modifications, we can now study it as a complex biological process of active self-modification.

Evolutionary change and phylogenetic relationships in light of horizontal gene transfer

Horizontal gene transfer has, over the past 25 years, become a part of evolutionary thinking. In the present paper I discuss horizontal gene transfer (HGT) in relation to contingency, natural selection, evolutionary change speed and the Tree-of-Life endeavour, with the aim of contributing to the understanding of the role of HGT in evolutionary processes. In addition, the challenges that HGT imposes on the current view of evolution are emphasized.

Comparative Genome Analysis Provides Insights into the Pathogenicity of Flavobacterium psychrophilum

Flavobacterium psychrophilum is a fish pathogen in salmonid aquaculture worldwide that causes cold water disease (CWD) and rainbow trout fry syndrome (RTFS). Comparative genome analyses of 11 F. psychrophilum isolates representing temporally and geographically distant populations were used to describe the F. psychrophilum pan-genome and to examine virulence factors, prophages, CRISPR arrays, and genomic islands present in the genomes. Analysis of the genomic DNA sequences were complemented with selected phenotypic characteristics of the strains. The pan genome analysis showed that F. psychrophilum could hold at least 3373 genes, while the core genome contained 1743 genes. On average, 67 new genes were detected for every new genome added to the analysis, indicating that F. psychrophilum possesses an open pan genome. The putative virulence factors were equally distributed among isolates, independent of geographic location, year of isolation and source of isolates. Only one prophage-related sequence was found which corresponded to the previously described prophage 6H, and appeared in 5 out of 11 isolates. CRISPR array analysis revealed two different loci with dissimilar spacer content, which only matched one sequence in the database, the temperate bacteriophage 6H. Genomic Islands (GIs) were identified in F. psychrophilum isolates 950106-1/1 and CSF 259–93, associated with toxins and antibiotic resistance. Finally, phenotypic characterization revealed a high degree of similarity among the strains with respect to biofilm formation and secretion of extracellular enzymes. Global scale dispersion of virulence factors in the genomes and the abilities for biofilm formation, hemolytic activity and secretion of extracellular enzymes among the strains suggested that F. psychrophilum isolates have a similar mode of action on adhesion, colonization and destruction of fish tissues across large spatial and temporal scales of occurrence. Overall, the genomic characterization and phenotypic properties may provide new insights to the mechanisms of pathogenicity in F. psychrophilum.

Abiotic Gene Transfer: Rare or Rampant?

Phylogenetic studies reveal that horizontal gene transfer (HGT) plays a prominent role in evolution and genetic variability of life. Five biotic mechanisms of HGT among prokaryotic organisms have been extensively characterized: conjugation, competence, transduction, gene transfer agent particles, and transitory fusion with recombination, but it is not known whether they can account for all natural HGT. It is even less clear how HGT could have occurred before any of these mechanisms had developed. Here, we consider contemporary conditions and experiments on microorganisms to estimate possible roles of abiotic HGT-currently and throughout evolution. Candidate mechanisms include freeze-and-thaw, microbeads-agitation, and electroporation-based transformation, and we posit that these laboratory techniques have analogues in nature acting as mechanisms of abiotic HGT: freeze-and-thaw cycles in polar waters, agitation by sand at foreshores and riverbeds, and lightning-triggered electroporation in near-surface aqueous habitats. We derive conservative order-of-magnitude estimates for rates of microorganisms subjected to freeze-and-thaw cycles, sand agitation, and lightning-triggered electroporation, at 10²⁴, 10¹⁹, and 10¹⁷ per year, respectively. Considering the yield of viable transformants, which is by far the highest in electroporation, we argue this may still favor lightning-triggered transformation over the other two mechanisms. Electroporation-based gene transfer also appears to be the most general of these abiotic candidates, and perhaps even of all known HGT mechanisms. Future studies should provide improved estimates of gene transfer rates and cell viability, currently and in the past, but to assess the importance of abiotic HGT in nature will likely require substantial progress-also in knowledge of biotic HGT.

Molecular Taxonomic Evidence for Two Distinct Genotypes of Mycobacterium yongonense via Genome-Based Phylogenetic Analysis

Recently, we introduced a distinct Mycobacterium intracellulare INT-5 genotype, distantly related to other genotypes of M. intracellulare (INT-1 to -4). The aim of this study is to determine the exact taxonomic status of the M. intracellulare INT-5 genotype via genome-based phylogenetic analysis. To this end, genome sequences of the two INT-5 strains, MOTT-H4Y and MOTT-36Y were compared with M. intracellulare ATCC 13950^T and Mycobacterium yongonense DSM 45126^T. Our phylogenetic analysis based on complete genome sequences, multi-locus sequence typing (MLST) of 35 target genes, and single nucleotide polymorphism (SNP) analysis indicated that the two INT-5 strains were more closely related to M. yongonense DSM 45126^T than the M. intracellulare strains. These results suggest their taxonomic transfer from M. intracellulare into M. yongonense. Finally, we selected 5 target genes (argH, dnaA, deaD, hsp65, and recF) and used SNPs for the identification of M. yongonese strains from other M. avium complex (MAC) strains. The application of the SNP analysis to 14 MAC clinical isolates enabled the selective identification of 4 M. yongonense clinical isolates from the other MACs. In conclusion, our genome-based phylogenetic analysis showed that the taxonomic status of two INT-5 strains, MOTT-H4Y and MOTT-36Y should be revised into M. yongonense. Our results also suggest that M. yongonense could be divided into 2 distinct genotypes (the Type I genotype with the M. parascrofulaceum rpoB gene and the Type II genotype with the M. intracellulare rpoB gene) depending on the presence of the lateral gene transfer of rpoB from M. parascrofulaceum.

Detection of Prokaryotic Genes in the Amphimedon queenslandica Genome

Horizontal gene transfer (HGT) is common between prokaryotes and phagotrophic eukaryotes. In metazoans, the scale and significance of HGT remains largely unexplored but is usually linked to a close association with parasites and endosymbionts. Marine sponges (Porifera), which host many microorganisms in their tissues and lack an isolated germ line, are potential carriers of genes transferred from prokaryotes. In this study, we identified a number of potential horizontally transferred genes within the genome of the sponge, Amphimedon queenslandica. We further identified homologs of some of these genes in other sponges. The transferred genes, most of which possess catalytic activity for carbohydrate or protein metabolism, have assimilated host genome characteristics and are actively expressed. The diversity of functions contributed by the horizontally transferred genes is likely an important factor in the adaptation and evolution of A. queenslandica. These findings highlight the potential importance of HGT on the success of sponges in diverse ecological niches.

Inferring Orthologs: Open Questions and Perspectives

With the increasing number of sequenced genomes and their comparisons, the detection of orthologs is crucial for reliable functional annotation and evolutionary analyses of genes and species. Yet, the dynamic remodeling of genome content through gain, loss, transfer of genes, and segmental and whole-genome duplication hinders reliable orthology detection. Moreover, the lack of direct functional evidence and the questionable quality of some available genome sequences and annotations present additional difficulties to assess orthology. This article reviews the existing computational methods and their potential accuracy in the high-throughput era of genome sequencing and anticipates open questions in terms of methodology, reliability, and computation. Appropriate taxon sampling together with combination of methods based on similarity, phylogeny, synteny, and evolutionary knowledge that may help detecting speciation events appears to be the most accurate strategy. This review also raises perspectives on the potential determination of orthology throughout the whole species phylogeny.

Continuous Influx of Genetic Material from Host to Virus Populations

Many genes of large double-stranded DNA viruses have a cellular origin, suggesting that host-to-virus horizontal transfer (HT) of DNA is recurrent. Yet, the frequency of these transfers has never been assessed in viral populations. Here we used ultra-deep DNA sequencing of 21 baculovirus populations extracted from two moth species to show that a large diversity of moth DNA sequences (n = 86) can integrate into viral genomes during the course of a viral infection. The majority of the 86 different moth DNA sequences are transposable elements (TEs, n = 69) belonging to 10 superfamilies of DNA transposons and three superfamilies of retrotransposons. The remaining 17 sequences are moth sequences of unknown nature. In addition to bona fide DNA transposition, we uncover microhomology-mediated recombination as a mechanism explaining integration of moth sequences into viral genomes. Many sequences integrated multiple times at multiple positions along the viral genome. We detected a total of 27,504 insertions of moth sequences in the 21 viral populations and we calculate that on average, 4.8% of viruses harbor at least one moth sequence in these populations. Despite this substantial proportion, no insertion of moth DNA was maintained in any viral population after 10 successive infection cycles. Hence, there is a constant turnover of host DNA inserted into viral genomes each time the virus infects a moth. Finally, we found that at least 21 of the moth TEs integrated into viral genomes underwent repeated horizontal transfers between various insect species, including some lepidopterans susceptible to baculoviruses. Our results identify host DNA influx as a potent source of genetic diversity in viral populations. They also support a role for baculoviruses as vectors of DNA HT between insects, and call for an evaluation of possible gene or TE spread when using viruses as biopesticides or gene delivery vectors.

While gene exchange is known to occur between viruses and their hosts, this phenomenon has never been studied at the level of the viral population. Here we report that each time a virus from the Baculoviridae family infects a moth, a large number (dozens to hundreds) and high diversity of moth DNA sequences (86 different sequences) can integrate into replicating viral genomes. These findings show that viral populations carry a measurable load of host DNA sequences, further supporting the role of viruses as vectors of horizontal transfer of DNA between insect species. The potential uncontrolled gene spread associated with the use of viruses produced in insect cells as gene delivery vectors and/or biopesticides should therefore be evaluated.

The Immune Syntax Revisited: Opening New Windows on Language Evolution

Recent research has added new dimensions to our understanding of classical evolution, according to which evolutionary novelties result from gene mutations inherited from parents to offspring. Language is surely one such novelty. Together with specific changes in our genome and epigenome, we suggest that two other (related) mechanisms may have contributed to the brain rewiring underlying human cognitive evolution and, specifically, the changes in brain connectivity that prompted the emergence of our species-specific linguistic abilities: the horizontal transfer of genetic material by viral and non-viral vectors and the brain/immune system crosstalk (more generally, the dialogue between the microbiota, the immune system, and the brain).

Evolutionary Ancestry of Eukaryotic Protein Kinases and Choline Kinases

The reversible phosphorylation of proteins catalyzed by protein kinases in eukaryotes supports an important role for eukaryotic protein kinases (ePKs) in the emergence of nucleated cells in the third superkingdom of life. Choline kinases (ChKs) could also be critical in the early evolution of eukaryotes, because of their function in the biosynthesis of phosphatidylcholine, which is unique to eukaryotic membranes. However, the genomic origins of ePKs and ChKs are unclear. The high degeneracy of protein sequences and broad expansion of ePK families have made this fundamental question difficult to answer. In this study, we identified two class-I aminoacyl-tRNA synthetases with high similarities to consensus amino acid sequences of human protein-serine/threonine kinases. Comparisons of primary and tertiary structures supported that ePKs and ChKs evolved from a common ancestor related to glutaminyl aminoacyl-tRNA synthetases, which may have been one of the key factors in the successful of emergence of ancient eukaryotic cells from bacterial colonies.

The population genetics of drug resistance evolution in natural populations of viral, bacterial and eukaryotic pathogens

Drug resistance is a costly consequence of pathogen evolution and a major concern in public health. In this review, we show how population genetics can be used to study the evolution of drug resistance and also how drug resistance evolution is informative as an evolutionary model system. We highlight five examples from diverse organisms with particular focus on: (i) identifying drug resistance loci in the malaria parasite Plasmodium falciparum using the genomic signatures of selective sweeps, (ii) determining the role of epistasis in drug resistance evolution in influenza, (iii) quantifying the role of standing genetic variation in the evolution of drug resistance in HIV, (iv) using drug resistance mutations to study clonal interference dynamics in tuberculosis and (v) analysing the population structure of the core and accessory genome of Staphylococcus aureus to understand the spread of methicillin resistance. Throughout this review, we discuss the uses of sequence data and population genetic theory in studying the evolution of drug resistance.

Insect-induced effects on plants and possible effectors used by galling and leaf-mining insects to manipulate their host-plant

Gall-inducing insects are iconic examples in the manipulation and reprogramming of plant development, inducing spectacular morphological and physiological changes of host-plant tissues within which the insect feeds and grows. Despite decades of research, effectors involved in gall induction and basic mechanisms of gall formation remain unknown. Recent research suggests that some aspects of the plant manipulation shown by gall-inducers may be shared with other insect herbivorous life histories. Here, we illustrate similarities and contrasts by reviewing current knowledge of metabolic and morphological effects induced on plants by gall-inducing and leaf-mining insects, and ask whether leaf-miners can also be considered to be plant reprogrammers. We review key plant functions targeted by various plant reprogrammers, including plant-manipulating insects and nematodes, and functionally characterize insect herbivore-derived effectors to provide a broader understanding of possible mechanisms used in host-plant manipulation. Consequences of plant reprogramming in terms of ecology, coevolution and diversification of plant-manipulating insects are also discussed.

New Structural Templates for Clinically Validated and Novel Targets in Antimicrobial Drug Research and Development

The development of bacterial resistance against current antibiotic drugs necessitates a continuous renewal of the arsenal of efficacious drugs. This imperative has not been met by the output of antibiotic research and development of the past decades for various reasons, including the declining efforts of large pharma companies in this area. Moreover, the majority of novel antibiotics are chemical derivatives of existing structures that represent mostly step innovations, implying that the available chemical space may be exhausted. This review negates this impression by showcasing recent achievements in lead finding and optimization of antibiotics that have novel or unexplored chemical structures. Not surprisingly, many of the novel structural templates like teixobactins, lysocin, griselimycin, or the albicidin/cystobactamid pair were discovered from natural sources. Additional compounds were obtained from the screening of synthetic libraries and chemical synthesis, including the gyrase-inhibiting NTBI's and spiropyrimidinetrione, the tarocin and targocil inhibitors of wall teichoic acid synthesis, or the boronates and diazabicyclo[3.2.1]octane as novel β-lactamase inhibitors. A motif that is common to most clinically validated antibiotics is that they address hotspots in complex biosynthetic machineries, whose functioning is essential for the bacterial cell. Therefore, an introduction to the biological targets-cell wall synthesis, topoisomerases, the DNA sliding clamp, and membrane-bound electron transport-is given for each of the leads presented here.

Deep sequencing revealed molecular signature of horizontal gene transfer of plant like transcripts in the mosquito Anopheles culicifacies: an evolutionary puzzle

In prokaryotes, horizontal gene transfer (HGT) has been regarded as an important evolutionary drive to acquire and retain beneficial genes for their survival in diverse ecologies. However, in eukaryotes, the functional role of HGTs remains questionable, although current genomic tools are providing increased evidence of acquisition of novel traits within non-mating metazoan species. Here, we provide another transcriptomic evidence for the acquisition of massive plant genes in the mosquito, Anopheles culicifacies. Our multiple experimental validations including genomic PCR, RT-PCR, real-time PCR, immuno-blotting and immuno-florescence microscopy, confirmed that plant like transcripts (PLTs) are of mosquito origin and may encode functional proteins. A comprehensive molecular analysis of the PLTs and ongoing metagenomic analysis of salivary microbiome provide initial clues that mosquitoes may have survival benefits through the acquisition of nuclear as well as chloroplast encoded plant genes. Our findings of PLTs further support the similar questionable observation of HGTs in other higher organisms, which is still a controversial and debatable issue in the community of evolutionists. We believe future understanding of the underlying mechanism of the feeding associated molecular responses may shed new insights in the functional role of PLTs in the mosquito.

What's in your next-generation sequence data? An exploration of unmapped DNA and RNA sequence reads from the bovine reference individual

BACKGROUND

Next-generation sequencing projects commonly commence by aligning reads to a reference genome assembly. While improvements in alignment algorithms and computational hardware have greatly enhanced the efficiency and accuracy of alignments, a significant percentage of reads often remain unmapped.

RESULTS

We generated de novo assemblies of unmapped reads from the DNA and RNA sequencing of the Bos taurus reference individual and identified the closest matching sequence to each contig by alignment to the NCBI non-redundant nucleotide database using BLAST. As expected, many of these contigs represent vertebrate sequence that is absent, incomplete, or misassembled in the UMD3.1 reference assembly. However, numerous additional contigs represent invertebrate species. Most prominent were several species of Spirurid nematodes and a blood-borne parasite, Babesia bigemina. These species are either not present in the US or are not known to infect taurine cattle and the reference animal appears to have been host to unsequenced sister species.

CONCLUSIONS

We demonstrate the importance of exploring unmapped reads to ascertain sequences that are either absent or misassembled in the reference assembly and for detecting sequences indicative of parasitic or commensal organisms.

Horizontal functional gene transfer from bacteria to fishes

Invertebrates can acquire functional genes via horizontal gene transfer (HGT) from bacteria but fishes are not known to do so. We provide the first reliable evidence of one HGT event from marine bacteria to fishes. The HGT appears to have occurred after emergence of the teleosts. The transferred gene is expressed and regulated developmentally. Its successful integration and expression may change the genetic and metabolic repertoire of fishes. In addition, this gene contains conserved domains and similar tertiary structures in fishes and their putative donor bacteria. Thus, it may function similarly in both groups. Evolutionary analyses indicate that it evolved under purifying selection, further indicating its conserved function. We document the first likely case of HGT of functional gene from prokaryote to fishes. This discovery certifies that HGT can influence vertebrate evolution.

Evolutionary formation of gene clusters by reorganization: the meleagrin/roquefortine paradigm in different fungi

The biosynthesis of secondary metabolites in fungi is catalyzed by enzymes encoded by genes linked in clusters that are frequently co-regulated at the transcriptional level. Formation of gene clusters may take place by de novo assembly of genes recruited from other cellular functions, but also novel gene clusters are formed by reorganization of progenitor clusters and are distributed by horizontal gene transfer. This article reviews (i) the published information on the roquefortine/meleagrin/neoxaline gene clusters of Penicillium chrysogenum (Penicillium rubens) and the short roquefortine cluster of Penicillium roqueforti, and (ii) the correlation of the genes present in those clusters with the enzymes and metabolites derived from these pathways. The P. chrysogenum roq/mel cluster consists of seven genes and includes a gene (roqT) encoding a 12-TMS transporter protein of the MFS family. Interestingly, the orthologous P. roquefortine gene cluster has only four genes and the roqT gene is present as a residual pseudogene that encodes only small peptides. Two of the genes present in the central region of the P. chrysogenum roq/mel cluster have been lost during the evolutionary formation of the short cluster and the order of the structural genes in the cluster has been rearranged. The two lost genes encode a N1 atom hydroxylase (nox) and a roquefortine scaffold-reorganizing oxygenase (sro). As a consequence P. roqueforti has lost the ability to convert the roquefortine-type carbon skeleton to the glandicoline/meleagrin-type scaffold and is unable to produce glandicoline B, meleagrin and neoxaline. The loss of this genetic information is not recent and occurred probably millions of years ago when a progenitor Penicillium strain got adapted to life in a few rich habitats such as cheese, fermented cereal grains or silage. P. roqueforti may be considered as a "domesticated" variant of a progenitor common to contemporary P. chrysogenum and related Penicillia.

Community rescue in experimental metacommunities

The conditions that allow biodiversity to recover following severe environmental degradation are poorly understood. We studied community rescue, the recovery of a viable community through the evolutionary rescue of many populations within an evolving community, in metacommunities of soil microbes adapting to a herbicide. The metacommunities occupied a landscape of crossed spatial gradients of the herbicide (Dalapon) and a resource (glucose), whereas their constituent communities were either isolated or connected by dispersal. The spread of adapted communities across the landscape and the persistence of communities when that landscape was degraded were strongly promoted by dispersal, and the capacity to adapt to lethal stress was also related to community size and initial diversity. After abrupt and lethal stress, community rescue was most frequent in communities that had previously experienced sublethal levels of stress and had been connected by dispersal. Community rescue occurred through the evolutionary rescue of both initially common taxa, which remained common, and of initially rare taxa, which grew to dominate the evolved community. Community rescue may allow productivity and biodiversity to recover from severe environmental degradation.

Population genomics supports baculoviruses as vectors of horizontal transfer of insect transposons

Horizontal transfer (HT) of DNA is an important factor shaping eukaryote evolution. Although several hundreds of eukaryote-to-eukaryote HTs of transposable elements (TEs) have been reported, the vectors underlying these transfers remain elusive. Here, we show that multiple copies of two TEs from the cabbage looper (Trichoplusia ni) transposed in vivo into genomes of the baculovirus Autographa californica multiple nucleopolyhedrovirus (AcMNPV) during caterpillar infection. We further demonstrate that both TEs underwent recent HT between several sympatric moth species (T. ni, Manduca sexta, Helicoverpa spp.) showing different degrees of susceptibility to AcMNPV. Based on two independent population genomics data sets (reaching a total coverage >330,000X), we report a frequency of one moth TE in ~8,500 AcMNPV genomes. Together, our results provide strong support for the role of viruses as vectors of TE HT between animals, and they call for a systematic evaluation of the frequency and impact of virus-mediated HT on the evolution of host genomes.

Horizontal transfer of DNA is common among eukaryotes but the vectors involved remain elusive. Here, Gilbert et al. show high frequency of in vivo transposition from the cabbage looper moth into genomes of a baculovirus, suggesting that viruses can act as vectors of horizontal transfer between animals.

Multiple recent horizontal transfers of a large genomic region in cheese making fungi

While the extent and impact of horizontal transfers in prokaryotes are widely acknowledged, their importance to the eukaryotic kingdom is unclear and thought by many to be anecdotal. Here we report multiple recent transfers of a huge genomic island between Penicillium spp. found in the food environment. Sequencing of the two leading filamentous fungi used in cheese making, P. roqueforti and P. camemberti, and comparison with the penicillin producer P. rubens reveals a 575 kb long genomic island in P. roqueforti—called Wallaby—present as identical fragments at non-homologous loci in P. camemberti and P. rubens. Wallaby is detected in Penicillium collections exclusively in strains from food environments. Wallaby encompasses about 250 predicted genes, some of which are probably involved in competition with microorganisms. The occurrence of multiple recent eukaryotic transfers in the food environment provides strong evidence for the importance of this understudied and probably underestimated phenomenon in eukaryotes.

Horizontal gene transfers are known to play an important role in prokaryote evolution but their impact and prevalence in eukaryotes is less clear. Here, the authors sequence the genomes of cheese making fungi P. roqueforti and P. camemberti, and provide evidence for recent horizontal transfers of a large genomic region.

Synthesis of phylogeny and taxonomy into a comprehensive tree of life

Reconstructing the phylogenetic relationships that unite all lineages (the tree of life) is a grand challenge. The paucity of homologous character data across disparately related lineages currently renders direct phylogenetic inference untenable. To reconstruct a comprehensive tree of life, we therefore synthesized published phylogenies, together with taxonomic classifications for taxa never incorporated into a phylogeny. We present a draft tree containing 2.3 million tips-the Open Tree of Life. Realization of this tree required the assembly of two additional community resources: (i) a comprehensive global reference taxonomy and (ii) a database of published phylogenetic trees mapped to this taxonomy. Our open source framework facilitates community comment and contribution, enabling the tree to be continuously updated when new phylogenetic and taxonomic data become digitally available. Although data coverage and phylogenetic conflict across the Open Tree of Life illuminate gaps in both the underlying data available for phylogenetic reconstruction and the publication of trees as digital objects, the tree provides a compelling starting point for community contribution. This comprehensive tree will fuel fundamental research on the nature of biological diversity, ultimately providing up-to-date phylogenies for downstream applications in comparative biology, ecology, conservation biology, climate change, agriculture, and genomics.

Recurrent Domestication by Lepidoptera of Genes from Their Parasites Mediated by Bracoviruses

Bracoviruses are symbiotic viruses associated with tens of thousands of species of parasitic wasps that develop within the body of lepidopteran hosts and that collectively parasitize caterpillars of virtually every lepidopteran species. Viral particles are produced in the wasp ovaries and injected into host larvae with the wasp eggs. Once in the host body, the viral DNA circles enclosed in the particles integrate into lepidopteran host cell DNA. Here we show that bracovirus DNA sequences have been inserted repeatedly into lepidopteran genomes, indicating this viral DNA can also enter germline cells. The original mode of Horizontal Gene Transfer (HGT) unveiled here is based on the integrative properties of an endogenous virus that has evolved as a gene transfer agent within parasitic wasp genomes for ≈100 million years. Among the bracovirus genes thus transferred, a phylogenetic analysis indicated that those encoding C-type-lectins most likely originated from the wasp gene set, showing that a bracovirus-mediated gene flux exists between the 2 insect orders Hymenoptera and Lepidoptera. Furthermore, the acquisition of bracovirus sequences that can be expressed by Lepidoptera has resulted in the domestication of several genes that could result in adaptive advantages for the host. Indeed, functional analyses suggest that two of the acquired genes could have a protective role against a common pathogen in the field, baculovirus. From these results, we hypothesize that bracovirus-mediated HGT has played an important role in the evolutionary arms race between Lepidoptera and their pathogens.

Eukaryotes are generally thought to evolve mainly through the modification of existing genetic information. However, evidence of horizontal gene transfer (HGT) in eukaryotes-the accidental acquisition of a novel gene from another species, allowing acquisition of novel traits—is now recognized as an important factor in their evolution. We show here that in several lineages, lepidopteran genomes have acquired genes from a bracovirus that is symbiotically used by parasitic wasps to inhibit caterpillar host immune defences. Integration of parts of the viral genome into host caterpillar DNA strongly suggests that integration can sporadically occur in the germline, leading to the production of lepidopteran lineages that harbor bracovirus sequences. Moreover, some of the transferred bracovirus genes reported here originate from the wasp genome, demonstrating that a gene flux exists between the two insect orders Hymenoptera and Lepidoptera that diverged ≈300 MYA. As bracovirus gene organisation has evolved to allow expression in Lepidoptera, these transferred genes can be readily domesticated. Additionally, we present functional analyses suggesting that some of the acquired genes confer to caterpillars a protection toward baculovirus, a very common pathogen in the field. This phenomenon may have implications for understanding how caterpillars acquire resistance against baculoviruses used in biological control.

Microbial Speciation

What are species? How do they arise? These questions are not easy to answer and have been particularly controversial in microbiology. Yet, for those microbiologists studying environmental questions or dealing with clinical issues, the ability to name and recognize species, widely considered the fundamental units of ecology, can be practically useful. On a more fundamental level, the speciation problem, the focus here, is more mechanistic and conceptual. What is the origin of microbial species, and what evolutionary and ecological mechanisms keep them separate once they begin to diverge? To what extent are these mechanisms universal across diverse types of microbes, and more broadly across the entire the tree of life? Here, we propose that microbial speciation must be viewed in light of gene flow, which defines units of genetic similarity, and of natural selection, which defines units of phenotype and ecological function. We discuss to what extent ecological and genetic units overlap to form cohesive populations in the wild, based on recent evolutionary modeling and population genomics studies. These studies suggest a continuous "speciation spectrum," which microbial populations traverse in different ways depending on their balance of gene flow and natural selection.

Distribution of putative xenogeneic silencers in prokaryote genomes

Gene silencing is an important function as it keeps newly acquired foreign DNA repressed, thereby avoiding possible deleterious effects in the host organism. Known transcriptional regulators associated with this process are called xenogeneic silencers (XS) and belong to either the H-NS, Lsr2, MvaT or Rok families. In the work described here we looked for XS-like regulators and their distribution in prokaryotic organisms was evaluated. Our analysis showed that putative XS regulators similar to H-NS, Lsr2, MvaT or Rok are present only in bacteria (31.7%). This does not exclude the existence of alternative XS in the rest of the organisms analyzed. Additionally, of the four XS groups evaluated in this work, those from the H-NS family have diversified more than the other groups. In order to compare the distribution of these putative XS regulators we also searched for other nucleoid-associated proteins (NAPs) not included in this group such as Fis, EbfC/YbaB, HU/IHF and Alba. Results showed that NAPs from the Fis, EbfC/YbaB, HU/IHF and Alba families are widely (94%) distributed among prokaryotes. These NAPs were found in multiple combinations with or without XS-like proteins. In regard with XS regulators, results showed that only XS proteins from one family were found in those organisms containing them. This suggests specificity for this type of regulators and their corresponding genomes.

The search for a promising cell factory system for production of edible vaccine

Despite worldwide vaccination against devastating diseases for decades, millions of children in remote and impoverished regions of the globe die every year from vaccine-preventable infectious diseases. The reasons for incomplete coverage of vaccination programs are based in part on the relatively high costs of conventional vaccinations, including mass production, refrigeration, transportation, and training as well as funding personnel for their administration. Plant-based edible vaccines (PEVs) have been introduced as a revolutionary cost-effective vaccination modality. However, they suffer from major deficiencies that have restricted their application to bench-scale. This article discusses the deficiencies of PEVs and also provides concise overview on the health-promoting, biological and biotechnological features of spirulina (Arthrospira). In short, we envision that spirulina could be considered as a potential alternative biofactory system to the plants toward the production of edible vaccines in high-yield with low-costs that other hosts cannot yet offer.

Differential gene retention as an evolutionary mechanism to generate biodiversity and adaptation in yeasts

The evolutionary history of the characters underlying the adaptation of microorganisms to food and biotechnological uses is poorly understood. We undertook comparative genomics to investigate evolutionary relationships of the dairy yeast Geotrichum candidum within Saccharomycotina. Surprisingly, a remarkable proportion of genes showed discordant phylogenies, clustering with the filamentous fungus subphylum (Pezizomycotina), rather than the yeast subphylum (Saccharomycotina), of the Ascomycota. These genes appear not to be the result of Horizontal Gene Transfer (HGT), but to have been specifically retained by G. candidum after the filamentous fungi-yeasts split concomitant with the yeasts' genome contraction. We refer to these genes as SRAGs (Specifically Retained Ancestral Genes), having been lost by all or nearly all other yeasts, and thus contributing to the phenotypic specificity of lineages. SRAG functions include lipases consistent with a role in cheese making and novel endoglucanases associated with degradation of plant material. Similar gene retention was observed in three other distantly related yeasts representative of this ecologically diverse subphylum. The phenomenon thus appears to be widespread in the Saccharomycotina and argues that, alongside neo-functionalization following gene duplication and HGT, specific gene retention must be recognized as an important mechanism for generation of biodiversity and adaptation in yeasts.

Frequent, independent transfers of a catabolic gene from bacteria to contrasted filamentous eukaryotes

Even genetically distant prokaryotes can exchange genes between them, and these horizontal gene transfer events play a central role in adaptation and evolution. While this was long thought to be restricted to prokaryotes, certain eukaryotes have acquired genes of bacterial origin. However, gene acquisitions in eukaryotes are thought to be much less important in magnitude than in prokaryotes. Here, we describe the complex evolutionary history of a bacterial catabolic gene that has been transferred repeatedly from different bacterial phyla to stramenopiles and fungi. Indeed, phylogenomic analysis pointed to multiple acquisitions of the gene in these filamentous eukaryotes-as many as 15 different events for 65 microeukaryotes. Furthermore, once transferred, this gene acquired introns and was found expressed in mRNA databases for most recipients. Our results show that effective inter-domain transfers and subsequent adaptation of a prokaryotic gene in eukaryotic cells can happen at an unprecedented magnitude.

The Ty1-copia LTR retroelement family PARTC is highly conserved in conifers over 200 MY of evolution

Long Terminal Repeat retroelements (LTR-RTs) are a major component of many plant genomes. Although well studied and described in angiosperms, their features and dynamics are poorly understood in gymnosperms. Representative complete copies of a Ty1-copia element isolate in Picea abies and named PARTC were identified in six other conifer species (Picea glauca, Pinus sylvestris, Pinus taeda, Abies sibirica, Taxus baccata and Juniperus communis) covering more than 200 million years of evolution. Here we characterized the structure of this element, assessed its abundance across conifers, studied the modes and timing of its amplification, and evaluated the degree of conservation of its extant copies at nucleotide level over distant species. We demonstrated that the element is ancient, abundant, widespread and its paralogous copies are present in the genera Picea, Pinus and Abies as an LTR-RT family. The amplification leading to the extant copies of PARTC occurred over long evolutionary times spanning 10s of MY and mostly took place after the speciation of the conifers analyzed. The level of conservation of PARTC is striking and may be explained by low substitution rates and limited removal mechanisms for LTR-RTs. These PARTC features and dynamics are representative of a more general scenario for LTR-RTs in gymnosperms quite different from that characterizing the vast majority of LTR-RT elements in angiosperms.

Rational design of antibiotic treatment plans: a treatment strategy for managing evolution and reversing resistance

The development of reliable methods for restoring susceptibility after antibiotic resistance arises has proven elusive. A greater understanding of the relationship between antibiotic administration and the evolution of resistance is key to overcoming this challenge. Here we present a data-driven mathematical approach for developing antibiotic treatment plans that can reverse the evolution of antibiotic resistance determinants. We have generated adaptive landscapes for 16 genotypes of the TEM β-lactamase that vary from the wild type genotype “TEM-1” through all combinations of four amino acid substitutions. We determined the growth rate of each genotype when treated with each of 15 β-lactam antibiotics. By using growth rates as a measure of fitness, we computed the probability of each amino acid substitution in each β-lactam treatment using two different models named the Correlated Probability Model (CPM) and the Equal Probability Model (EPM). We then performed an exhaustive search through the 15 treatments for substitution paths leading from each of the 16 genotypes back to the wild type TEM-1. We identified optimized treatment paths that returned the highest probabilities of selecting for reversions of amino acid substitutions and returning TEM to the wild type state. For the CPM model, the optimized probabilities ranged between 0.6 and 1.0. For the EPM model, the optimized probabilities ranged between 0.38 and 1.0. For cyclical CPM treatment plans in which the starting and ending genotype was the wild type, the probabilities were between 0.62 and 0.7. Overall this study shows that there is promise for reversing the evolution of resistance through antibiotic treatment plans.

Vegetative hyphal fusion and subsequent nuclear behavior in Epichloë grass endophytes

Epichloë species (including the former genus Neotyphodium) are fungal symbionts of many agronomically important forage grasses, and provide their grass hosts with protection from a wide range of biotic and abiotic stresses. Epichloë species include many interspecific hybrids with allodiploid-like genomes, which may provide the potential for combined traits or recombination to generate new traits. Though circumstantial evidence suggests that such interspecific hybrids might have arisen from nuclear fusion events following vegetative hyphal fusion between different Epichloë strains, this hypothesis has not been addressed empirically. Here, we investigated vegetative hyphal fusion and subsequent nuclear behavior in Epichloë species. A majority of Epichloë strains, especially those having a sexual stage, underwent self vegetative hyphal fusion. Vegetative fusion also occurred between two hyphae from different Epichloë strains. Though Epichloë spp. are uninucleate fungi, hyphal fusion resulted in two nuclei stably sharing the same cytoplasm, which might ultimately lead to nuclear fusion. In addition, protoplast fusion experiments gave rise to uninucleate putative hybrids, which apparently had two markers, one from each parent within the same nucleus. These results are consistent with the notion that interspecific hybrids arise from vegetative hyphal fusion. However, we also discuss additional factors, such as post-hybridization selection, that may be important to explain the recognized prevalence of hybrids in Epichloë species.

Expression of multiple horizontally acquired genes is a hallmark of both vertebrate and invertebrate genomes

Background

A fundamental concept in biology is that heritable material, DNA, is passed from parent to offspring, a process called vertical gene transfer. An alternative mechanism of gene acquisition is through horizontal gene transfer (HGT), which involves movement of genetic material between different species. HGT is well-known in single-celled organisms such as bacteria, but its existence in higher organisms, including animals, is less well established, and is controversial in humans.

Results

We have taken advantage of the recent availability of a sufficient number of high-quality genomes and associated transcriptomes to carry out a detailed examination of HGT in 26 animal species (10 primates, 12 flies and four nematodes) and a simplified analysis in a further 14 vertebrates. Genome-wide comparative and phylogenetic analyses show that HGT in animals typically gives rise to tens or hundreds of active ‘foreign’ genes, largely concerned with metabolism. Our analyses suggest that while fruit flies and nematodes have continued to acquire foreign genes throughout their evolution, humans and other primates have gained relatively few since their common ancestor. We also resolve the controversy surrounding previous evidence of HGT in humans and provide at least 33 new examples of horizontally acquired genes.

Conclusions

We argue that HGT has occurred, and continues to occur, on a previously unsuspected scale in metazoans and is likely to have contributed to biochemical diversification during animal evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0607-3) contains supplementary material, which is available to authorized users.

Roles for Synonymous Codon Usage in Protein Biogenesis

Owing to the degeneracy of the genetic code, a protein sequence can be encoded by many different synonymous mRNA coding sequences. Synonymous codon usage was once thought to be functionally neutral, but evidence now indicates it is shaped by evolutionary selection and affects other aspects of protein biogenesis beyond specifying the amino acid sequence of the protein. Synonymous rare codons, once thought to have only negative impacts on the speed and accuracy of translation, are now known to play an important role in diverse functions, including regulation of cotranslational folding, covalent modifications, secretion, and expression level. Mutations altering synonymous codon usage are linked to human diseases. However, much remains unknown about the molecular mechanisms connecting synonymous codon usage to efficient protein biogenesis and proper cell physiology. Here we review recent literature on the functional effects of codon usage, including bioinformatics approaches aimed at identifying general roles for synonymous codon usage.

Comparative analysis of transposable elements highlights mobilome diversity and evolution in vertebrates

Transposable elements (TEs) are major components of vertebrate genomes, with major roles in genome architecture and evolution. In order to characterize both common patterns and lineage-specific differences in TE content and TE evolution, we have compared the mobilomes of 23 vertebrate genomes, including 10 actinopterygian fish, 11 sarcopterygians, and 2 nonbony vertebrates. We found important variations in TE content (from 6% in the pufferfish tetraodon to 55% in zebrafish), with a more important relative contribution of TEs to genome size in fish than in mammals. Some TE superfamilies were found to be widespread in vertebrates, but most elements showed a more patchy distribution, indicative of multiple events of loss or gain. Interestingly, loss of major TE families was observed during the evolution of the sarcopterygian lineage, with a particularly strong reduction in TE diversity in birds and mammals. Phylogenetic trends in TE composition and activity were detected: Teleost fish genomes are dominated by DNA transposons and contain few ancient TE copies, while mammalian genomes have been predominantly shaped by nonlong terminal repeat retrotransposons, along with the persistence of older sequences. Differences were also found within lineages: The medaka fish genome underwent more recent TE amplification than the related platyfish, as observed for LINE retrotransposons in the mouse compared with the human genome. This study allows the identification of putative cases of horizontal transfer of TEs, and to tentatively infer the composition of the ancestral vertebrate mobilome. Taken together, the results obtained highlight the importance of TEs in the structure and evolution of vertebrate genomes, and demonstrate their major impact on genome diversity both between and within lineages.

Alignment-free genetic sequence comparisons: a review of recent approaches by word analysis

Modern sequencing and genome assembly technologies have provided a wealth of data, which will soon require an analysis by comparison for discovery. Sequence alignment, a fundamental task in bioinformatics research, may be used but with some caveats. Seminal techniques and methods from dynamic programming are proving ineffective for this work owing to their inherent computational expense when processing large amounts of sequence data. These methods are prone to giving misleading information because of genetic recombination, genetic shuffling and other inherent biological events. New approaches from information theory, frequency analysis and data compression are available and provide powerful alternatives to dynamic programming. These new methods are often preferred, as their algorithms are simpler and are not affected by synteny-related problems. In this review, we provide a detailed discussion of computational tools, which stem from alignment-free methods based on statistical analysis from word frequencies. We provide several clear examples to demonstrate applications and the interpretations over several different areas of alignment-free analysis such as base-base correlations, feature frequency profiles, compositional vectors, an improved string composition and the D2 statistic metric. Additionally, we provide detailed discussion and an example of analysis by Lempel-Ziv techniques from data compression.

Detection of atypical genes in virus families using a one-class SVM

BACKGROUND

The diversity of viruses, the absence of universally common genes in them, and their ability to act as carriers of genetic material make assessment of evolutionary paths of viral genes very difficult. One important factor contributing to this complexity is horizontal gene transfer.

RESULTS

We explore the possibility for the systematic identification of atypical genes within virus families, including viruses whose genome is not encoded by a double-stranded DNA. Our method is based on gene statistical features that differ in genes that were subject of recent horizontal gene transfer from those of the genome in which they are observed. We employ a one-class SVM approach to detect atypical genes within a virus family basing of their statistical signatures and without explicit knowledge of the source species. The simplicity of the statistical features used makes the method applicable to various viruses irrespective of their genome size or type.

CONCLUSIONS

On simulated data, the method can robustly identify alien genes irrespective of the coding nucleic acid found in a virus. It also compares well to results obtained in related studies for double-stranded DNA viruses. Its value in practice is confirmed by the identification of isolated examples of horizontal gene transfer events that have already been described in the literature. A Python package implementing the method and the results for the analyzed virus families are available at http://svm-agp.bioinf.mpi-inf.mpg.de.