A natural barrier to lateral gene transfer from prokaryotes to eukaryotes revealed from genomes: the 70 % rule

Glycosaminoglycan lyase: A new competition between bacteria and the pacific white shrimp Litopenaeus vannamei

Horizontal gene transfer (HGT) is an important evolutionary force in the formation of prokaryotic and eukaryotic genomes. In recent years, many HGT genes horizontally transferred from prokaryotes to eukaryotes have been reported, and most of them are present in arthropods. The Pacific white shrimp Litopenaeus vannamei, an important economic species of arthropod, has close relationships with bacteria, providing a platform for horizontal gene transfer (HGT). In this study, we analyzed bacteria-derived HGT based on a high-quality genome of L. vannamei via a homology search and phylogenetic analysis, and six HGT genes were identified. Among these six horizontally transferred genes, we found one gene (LOC113799989) that contains a bacterial chondroitinase AC structural domain and encodes an unknown glycosaminoglycan (GAG) lyase in L. vannamei. The real-time quantitative PCR results showed that the mRNA expression level of LOC113799989 was highest in the hepatopancreas and heart, and after stimulation by Vibrio parahaemolyticus, its mRNA expression level was rapidly up-regulated within 12 h. Furthermore, after injecting si-RNA and stimulation by V. parahaemolyticus, we found that the experimental group had a higher cumulative mortality rate in 48 h than the control group, indicating that the bacteria-derived GAG lyase can reduce the mortality of shrimp with respect to infection by V. parahaemolyticus and might be related to the resistance of shrimp to bacterial diseases. Our findings contribute to the study of the function of GAGs and provide new insights into GAG-related microbial pathogenesis and host defense mechanisms in arthropods.

Are Viruses Taxonomic Units? A Protein Domain and Loop-Centric Phylogenomic Assessment

Virus taxonomy uses a Linnaean-like subsumption hierarchy to classify viruses into taxonomic units at species and higher rank levels. Virus species are considered monophyletic groups of mobile genetic elements (MGEs) often delimited by the phylogenetic analysis of aligned genomic or metagenomic sequences. Taxonomic units are assumed to be independent organizational, functional and evolutionary units that follow a 'natural history' rationale. Here, I use phylogenomic and other arguments to show that viruses are not self-standing genetically-driven systems acting as evolutionary units. Instead, they are crucial components of holobionts, which are units of biological organization that dynamically integrate the genetics, epigenetic, physiological and functional properties of their co-evolving members. Remarkably, phylogenomic analyses show that viruses share protein domains and loops with cells throughout history via massive processes of reticulate evolution, helping spread evolutionary innovations across a wider taxonomic spectrum. Thus, viruses are not merely MGEs or microbes. Instead, their genomes and proteomes conduct cellularly integrated processes akin to those cataloged by the GO Consortium. This prompts the generation of compositional hierarchies that replace the 'is-a-kind-of' by a 'is-a-part-of' logic to better describe the mereology of integrated cellular and viral makeup. My analysis demands a new paradigm that integrates virus taxonomy into a modern evolutionarily centered taxonomy of organisms.

Database Bias in the Detection of Interdomain Horizontal Gene Transfer Events in Pezizomycotina

Horizontal gene transfer (HGT) is a widely acknowledged phenomenon in prokaryotes for generating genetic diversity. However, the impact of this process in eukaryotes, particularly interdomain HGT, is a topic of debate. Although there have been observed biases in interdomain HGT detection, little exploration has been conducted on the effects of imbalanced databases. In our study, we conducted experiments to assess how different databases affect the detection of interdomain HGT using proteomes from the Pezizomycotina fungal subphylum as our focus group. Our objective was to simulate the database imbalance commonly found in public biological databases, where bacterial and eukaryotic sequences are unevenly represented, and demonstrate that an increase in uploaded eukaryotic sequences leads to a decrease in predicted HGTs. For our experiments, four databases with varying proportions of eukaryotic sequences but consistent proportions of bacterial sequences were utilized. We observed a significant reduction in detected interdomain HGT candidates as the proportion of eukaryotes increased within the database. Our data suggest that the imbalance in databases bias the interdomain HGT detection and highlights challenges associated with confirming the presence of interdomain HGT among Pezizomycotina fungi and potentially other groups within Eukarya.

Modeling the mosaic structure of bacterial genomes to infer their evolutionary history

The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.

Dynamics of nucleic acid mobility

Advances in sequencing technologies and bioinformatic analyses are accelerating the quantity and quality of data from all domains of life. This rich resource has the potential to reveal a number of important incidences with respect to possible exchange of nucleic acids. Ancient events have impacted species evolution and adaptation to new ecological niches. However, we still lack a full picture of processes ongoing within and between somatic cells, gametes, and different organisms. We propose that events linked to acceptance of alien nucleic acids grossly could be divided into 2 main routes in plants: one, when plants are exposed to extreme challenges and, the second level, a more everyday or season-related stress incited by biotic or abiotic factors. Here, many events seem to comprise somatic cells. Are the transport and acceptance processes of alien sequences random or are there specific regulatory systems not yet fully understood? Following entrance into a new cell, a number of intracellular processes leading to chromosomal integration and function are required. Modification of nucleic acids and possibly exchange of sequences within a cell may also occur. Such fine-tune events are most likely very common. There are multiple questions that we will discuss concerning different types of vesicles and their roles in nucleic acid transport and possible intracellular sequence exchange between species.

Omics of an Enigmatic Marine Amoeba Uncovers Unprecedented Gene Trafficking from Giant Viruses and Provides Insights into Its Complex Life Cycle

Amoebozoa include lineages of diverse ecology, behavior, and morphology. They are assumed to encompass members with the largest genome sizes of all living things, yet genomic studies in the group are limited. Trichosphaerium, a polymorphic, multinucleate, marine amoeba with a complicated life cycle, has puzzled experts for over a century. In an effort to explore the genomic diversity and investigate extraordinary behavior observed among the Amoebozoa, we used integrated omics approaches to study this enigmatic marine amoeba. Omics data, including single-cell transcriptomics and cytological data, demonstrate that Trichosphaerium sp. possesses the complete meiosis toolkit genes. These genes are expressed in life stages of the amoeba including medium and large cells. The life cycle of Trichosphaerium sp. involves asexual processes via binary fission and multiple fragmentation of giant cells, as well as sexual-like processes involving genes implicated in sexual reproduction and polyploidization. These findings are in stark contrast to a life cycle previously reported for this amoeba. Despite the extreme morphological plasticity observed in Trichosphaerium, our genomic data showed that populations maintain a species-level intragenomic variation. A draft genome of Trichosphaerium indicates elevated lateral gene transfer (LGT) from bacteria and giant viruses. Gene trafficking in Trichosphaerium is the highest within Amoebozoa and among the highest in microbial eukaryotes.

Extremophilic red algae as models for understanding adaptation to hostile environments and the evolution of eukaryotic life on the early earth

Extremophiles have always garnered great interest because of their exotic lifestyles and ability to thrive at the physical limits of life. In hot springs environments, the Cyanidiophyceae red algae are the only photosynthetic eukaryotes able to live under extremely low pH (0-5) and relatively high temperature (35ºC to 63ºC). These extremophiles live as biofilms in the springs, inhabit acid soils near the hot springs, and form endolithic populations in the surrounding rocks. Cyanidiophyceae represent a remarkable source of knowledge about the evolution of extremophilic lifestyles and their genomes encode specialized enzymes that have applied uses. Here we review the evolutionary origin, taxonomy, genome biology, industrial applications, and use of Cyanidiophyceae as genetic models. Currently, Cyanidiophyceae comprise a single order (Cyanidiales), three families, four genera, and nine species, including the well-known Cyanidioschyzon merolae and Galdieria sulphuraria. These algae have small, gene-rich genomes that are analogous to those of prokaryotes they live and compete with. There are few spliceosomal introns and evidence exists for horizontal gene transfer as a driver of local adaptation to gain access to external fixed carbon and to extrude toxic metals. Cyanidiophyceae offer a variety of commercial opportunities such as phytoremediation to detoxify contaminated soils or waters and exploitation of their mixotrophic lifestyles to support the efficient production of bioproducts such as phycocyanin and floridosides. In terms of exobiology, Cyanidiophyceae are an ideal model system for understanding the evolutionary effects of foreign gene acquisition and the interactions between different organisms inhabiting the same harsh environment on the early Earth. Finally, we describe ongoing research with C. merolae genetics and summarize the unique insights they offer to the understanding of algal biology and evolution.

Metazoan tryptophan indole-lyase: Are they still active?

Tryptophan indole-lyase (TIL), also known as tryptophanase, is a pyridoxal-5'-phosphate dependent bacterial enzyme that catalyzes the reversible hydrolytic cleavage of l-tryptophan (l-Trp) to indole and ammonium pyruvate. TIL is also found in some metazoans, and they may have been acquired by horizontal gene transfer. In this study, two metazoans, Nematostella vectensis (starlet sea anemone) and Bradysia coprophila (fungus gnat) TILs were bacterially expressed and characterized. The kcat values of metazoan TILs were low, < 1/200 of the kcat of Escherichia coli TIL. By contrast, metazoan TILs showed lower Km values than the TILs of common bacteria, indicating that their affinity for l-Trp is higher than that of bacterial TILs. Analysis of a series of chimeric enzymes based on B. coprophila and bacterial TILs revealed that the low Km value of B. coprophila TIL is not accidental due to the substitution of a single residue, but is due to the cooperative effect of multiple residues. This suggests that high affinity for l-Trp was positively selected during the molecular evolution of metazoan TIL. This is the first report that metazoan TILs have low but obvious activity.

Antibiotic resistance gene profile in aerobic granular reactor under antibiotic stress: Can eukaryotic microalgae act as inhibiting factor?

Antibiotic resistance gene (ARG) pollution is critical environmental problem, and horizontal gene transfer acts as a driving evolutionary force. In theory, due to the phylogenetic distance between eukaryotes and prokaryotes, eukaryotic microalgae can be a natural barrier that plays a negative role in ARG transfer among the symbiotic bacteria to decrease ARG abundance in sludge during wastewater treatment. However, this hypothesis is far from proven and needs to be tested experimentally, so this study investigated the influence of eukaryote microalgae (Scenedesmus) on the ARG profile of symbiotic bacteria based on aerobic granular reactor. The results indicated that Scenedesmus symbiosis could affect ARG diversity of bacteria, and the detected numbers of ARG in aerobic granular sludge (AG) group and algae-bacteria granular consortia (AAG) group were 45-53 and 44-47, respectively. In terms of relative abundance, after target microalgae symbiosis, the total abundance of ARGs significantly decreased from 1.17 × 10°, 2.69 × 10° and 1.36 × 10^-1 to 6.53 × 10^-1, 9.64 × 10^-1 and 1.04 × 10^-1 in the systems with the addition of streptomycin, azithromycin and vancomycin, respectively (P < 0.05), yet there was no significant difference between AG and AAG under the stress of ampicillin, sulfamethazine and tetracycline (P > 0.05). Redundancy analysis showed that the eukaryotic microalgae were significant factor explaining the change in ARG relative abundance (P < 0.05), which contributed 15.3% of ARG variation. Furthermore, the results show that, except for the tetracycline treatment system, the total relative abundances of MGEs in the AAG under the stress of the other five antibiotics were 3.54 × 10^-2-7.13 × 10^-1, which were all significantly lower than those in the AG (8.38 × 10^-2-1.59 × 10°). There was a more significant positive correlation relationship between ARGs and mobile genetic elements (MGEs) than that between ARGs and dominated bacteria.

Comparative Genomics of Microsporidia

The microsporidia are a phylum of intracellular parasites that represent the eukaryotic cell in a state of extreme reduction, with genomes and metabolic capabilities embodying eukaryotic cells in arguably their most streamlined state. Over the past 20 years, microsporidian genomics has become a rapidly expanding field starting with sequencing of the genome of Encephalitozoon cuniculi, one of the first ever sequenced eukaryotes, to the current situation where we have access to the data from over 30 genomes across 20+ genera. Reaching back further in evolutionary history, to the point where microsporidia diverged from other eukaryotic lineages, we now also have genomic data for some of the closest known relatives of the microsporidia such as Rozella allomycis, Metchnikovella spp. and Amphiamblys sp. Data for these organisms allow us to better understand the genomic processes that shaped the emergence of the microsporidia as a group. These intensive genomic efforts have revealed some of the processes that have shaped microsporidian cells and genomes including patterns of genome expansions and contractions through gene gain and loss, whole genome duplication, differential patterns of invasion and purging of transposable elements. All these processes have been shown to occur across short and longer time scales to give rise to a phylum of parasites with dynamic genomes with a diversity of sizes and organisations.

Reconciling Asgardarchaeota Phylogenetic Proximity to Eukaryotes and Planctomycetes Cellular Features in the Evolution of Life

The relationship between the three domains of life—Archaea, Bacteria, and Eukarya—is one of Biology’s greatest mysteries. Current favored models imply two ancestral domains, Bacteria and Archaea, with eukaryotes originating within Archaea. This type of models has been supported by the recent description of the Asgardarchaeota, the closest prokaryotic relatives of eukaryotes. However, there are many problems associated with any scenarios implying that eukaryotes originated from within the Archaea, including genome mosaicism, phylogenies, the cellular organization of the Archaea, and their ancestral character. By contrast, all models of eukaryogenesis fail to consider two relevant discoveries: the detection of membrane coat proteins, and of phagocytosis-related processes in Planctomycetes, which are among the bacteria with the most developed endomembrane system.

Consideration of these often overlooked features and others found in Planctomycetes and related bacteria suggest an evolutionary model based on a single ancestral domain. In this model, the proximity of Asgard and eukaryotes is not rejected but instead, Asgard are considered as diverging away from a common ancestor instead of on the way toward the eukaryotic ancestor. This model based on a single ancestral domain solves most of the ambiguities associated with the ones based on two ancestral domains. The single-domain model is better suited to explain the origin and evolution of all three domains of life, blurring the distinctions between them. Support for this model as well as the opportunities that it presents not only for reinterpreting previous results, but also for planning future experiments, are explored.

Gene Duplications Trace Mitochondria to the Onset of Eukaryote Complexity

The last eukaryote common ancestor (LECA) possessed mitochondria and all key traits that make eukaryotic cells more complex than their prokaryotic ancestors, yet the timing of mitochondrial acquisition and the role of mitochondria in the origin of eukaryote complexity remain debated. Here, we report evidence from gene duplications in LECA indicating an early origin of mitochondria. Among 163,545 duplications in 24,571 gene trees spanning 150 sequenced eukaryotic genomes, we identify 713 gene duplication events that occurred in LECA. LECA's bacterial-derived genes include numerous mitochondrial functions and were duplicated significantly more often than archaeal-derived and eukaryote-specific genes. The surplus of bacterial-derived duplications in LECA most likely reflects the serial copying of genes from the mitochondrial endosymbiont to the archaeal host's chromosomes. Clustering, phylogenies and likelihood ratio tests for 22.4 million genes from 5,655 prokaryotic and 150 eukaryotic genomes reveal no evidence for lineage-specific gene acquisitions in eukaryotes, except from the plastid in the plant lineage. That finding, and the functions of bacterial genes duplicated in LECA, suggests that the bacterial genes in eukaryotes are acquisitions from the mitochondrion, followed by vertical gene evolution and differential loss across eukaryotic lineages, flanked by concomitant lateral gene transfer among prokaryotes. Overall, the data indicate that recurrent gene transfer via the copying of genes from a resident mitochondrial endosymbiont to archaeal host chromosomes preceded the onset of eukaryotic cellular complexity, favoring mitochondria-early over mitochondria-late hypotheses for eukaryote origin.

The past, present and future of the tree of life

The advent of comparative genomics in the late 1990s led to the discovery of extensive lateral gene transfer in prokaryotes. The resulting debate over whether life as a whole is best represented as a tree or a network has since given way to a general consensus in which trees and networks co-exist rather than stand in opposition. Embracing this consensus allows us to move beyond the question of which is true or false. The future of the tree of life debate lies in asking what trees and networks can, and should, do for science.

The macronuclear genome of anaerobic ciliate Entodinium caudatum reveals its biological features adapted to the distinct rumen environment

Entodinium caudatum is an anaerobic binucleated ciliate representing the most dominant protozoal species in the rumen. However, its biological features are largely unknown due to the inability to establish an axenic culture. In this study, we primally sequenced its macronucleus (MAC) genome to aid the understanding of its metabolism, physiology, ecology. We isolated the MAC of E. caudatum strain MZG-1 and sequenced the MAC genome using Illumina MiSeq, MinION, and PacBio RSII systems. De novo assembly of the MiSeq sequence reads followed with subsequent scaffolding with MinION and PacBio reads resulted in a draft MAC genome about 117 Mbp. A large number of carbohydrate-active enzymes were likely acquired through horizontal gene transfer. About 8.74% of the E. caudatum predicted proteome was predicted as proteases. The MAC genome of E. caudatum will help better understand its important roles in rumen carbohydrate metabolism, and interaction with other members of the rumen microbiome.

Acquisition of pcnB [poly(A) polymerase I] genes via horizontal transfer from the β, γ-Proteobacteria

Poly(A) polymerases (PAPs) and tRNA nucleotidyltransferases belong to a superfamily of nucleotidyltransferases and modify RNA 3'-ends. The product of the pcnB gene, PAP I, has been characterized in a few β-, γ- and δ-Proteobacteria. Using the PAP I signature sequence, putative PAPs were identified in bacterial species from the α- and ε-Proteobacteria and from four other bacterial phyla (Firmicutes, Actinobacteria, Bacteroidetes and Aquificae). Phylogenetic analysis, alien index and G+C content calculations strongly suggest that the PAPs in the species identified in this study arose by horizontal gene transfer from the β- and γ-Proteobacteria.

Horizontally acquired antibacterial genes associated with adaptive radiation of ladybird beetles

BACKGROUND

Horizontal gene transfer (HGT) has been documented in many herbivorous insects, conferring the ability to digest plant material and promoting their remarkable ecological diversification. Previous reports suggest HGT of antibacterial enzymes may have contributed to the insect immune response and limit bacterial growth. Carnivorous insects also display many evolutionary successful lineages, but in contrast to the plant feeders, the potential role of HGTs has been less well-studied.

RESULTS

Using genomic and transcriptomic data from 38 species of ladybird beetles, we identified a set of bacterial cell wall hydrolase (cwh) genes acquired by this group of beetles. Infection with Bacillus subtilis led to upregulated expression of these ladybird cwh genes, and their recombinantly produced proteins limited bacterial proliferation. Moreover, RNAi-mediated cwh knockdown led to downregulation of other antibacterial genes, indicating a role in antibacterial immune defense. cwh genes are rare in eukaryotes, but have been maintained in all tested Coccinellinae species, suggesting that this putative immune-related HGT event played a role in the evolution of this speciose subfamily of predominant predatory ladybirds.

CONCLUSION

Our work demonstrates that, in a manner analogous to HGT-facilitated plant feeding, enhanced immunity through HGT might have played a key role in the prey adaptation and niche expansion that promoted the diversification of carnivorous beetle lineages. We believe that this represents the first example of immune-related HGT in carnivorous insects with an association with a subsequent successful species radiation.

Comparative Genomics Reveals Novel Target Genes towards Specific Control of Plant-Parasitic Nematodes

Plant-parasitic nematodes cause extensive annual yield losses to worldwide agricultural production. Most cultivated plants have no known resistance against nematodes and the few bearing a resistance gene can be overcome by certain species. Chemical methods that have been deployed to control nematodes have largely been banned from use due to their poor specificity and high toxicity. Hence, there is an urgent need for the development of cleaner and more specific control methods. Recent advances in nematode genomics, including in phytoparasitic species, provide an unprecedented opportunity to identify genes and functions specific to these pests. Using phylogenomics, we compared 61 nematode genomes, including 16 for plant-parasitic species and identified more than 24,000 protein families specific to these parasites. In the genome of Meloidogyne incognita, one of the most devastating plant parasites, we found ca. 10,000 proteins with orthologs restricted only to phytoparasitic species and no further homology in protein databases. Among these phytoparasite-specific proteins, ca. 1000 shared the same properties as known secreted effectors involved in essential parasitic functions. Of these, 68 were novel and showed strong expression during the endophytic phase of the nematode life cycle, based on both RNA-seq and RT-qPCR analyses. Besides effector candidates, transcription-related and neuro-perception functions were enriched in phytoparasite-specific proteins, revealing interesting targets for nematode control methods. This phylogenomics analysis constitutes a unique resource for the further understanding of the genetic basis of nematode adaptation to phytoparasitism and for the development of more efficient control methods.

A spectrum of verticality across genes

Lateral gene transfer (LGT) has impacted prokaryotic genome evolution, yet the extent to which LGT compromises vertical evolution across individual genes and individual phyla is unknown, as are the factors that govern LGT frequency across genes. Estimating LGT frequency from tree comparisons is problematic when thousands of genomes are compared, because LGT becomes difficult to distinguish from phylogenetic artefacts. Here we report quantitative estimates for verticality across all genes and genomes, leveraging a well-known property of phylogenetic inference: phylogeny works best at the tips of trees. From terminal (tip) phylum level relationships, we calculate the verticality for 19,050,992 genes from 101,422 clusters in 5,655 prokaryotic genomes and rank them by their verticality. Among functional classes, translation, followed by nucleotide and cofactor biosynthesis, and DNA replication and repair are the most vertical. The most vertically evolving lineages are those rich in ecological specialists such as Acidithiobacilli, Chlamydiae, Chlorobi and Methanococcales. Lineages most affected by LGT are the α-, β-, γ-, and δ- classes of Proteobacteria and the Firmicutes. The 2,587 eukaryotic clusters in our sample having prokaryotic homologues fail to reject eukaryotic monophyly using the likelihood ratio test. The low verticality of α-proteobacterial and cyanobacterial genomes requires only three partners-an archaeal host, a mitochondrial symbiont, and a plastid ancestor-each with mosaic chromosomes, to directly account for the prokaryotic origin of eukaryotic genes. In terms of phylogeny, the 100 most vertically evolving prokaryotic genes are neither representative nor predictive for the remaining 97% of an average genome. In search of factors that govern LGT frequency, we find a simple but natural principle: Verticality correlates strongly with gene distribution density, LGT being least likely for intruding genes that must replace a preexisting homologue in recipient chromosomes. LGT is most likely for novel genetic material, intruding genes that encounter no competing copy.

Genome structure and content of the rice root-knot nematode (Meloidogyne graminicola)

Discovered in the 1960s, Meloidogyne graminicola is a root-knot nematode species considered as a major threat to rice production. Yet, its origin, genomic structure, and intraspecific diversity are poorly understood. So far, such studies have been limited by the unavailability of a sufficiently complete and well-assembled genome. In this study, using a combination of Oxford Nanopore Technologies and Illumina sequencing data, we generated a highly contiguous reference genome (283 scaffolds with an N50 length of 294 kb, totaling 41.5 Mb). The completeness scores of our assembly are among the highest currently published for Meloidogyne genomes. We predicted 10,284 protein-coding genes spanning 75.5% of the genome. Among them, 67 are identified as possibly originating from horizontal gene transfers (mostly from bacteria), which supposedly contribute to nematode infection, nutrient processing, and plant defense manipulation. Besides, we detected 575 canonical transposable elements (TEs) belonging to seven orders and spanning 2.61% of the genome. These TEs might promote genomic plasticity putatively related to the evolution of M. graminicola parasitism. This high-quality genome assembly constitutes a major improvement regarding previously available versions and represents a valuable molecular resource for future phylogenomic studies of Meloidogyne species. In particular, this will foster comparative genomic studies to trace back the evolutionary history of M. graminicola and its closest relatives.

Whitefly genomes contain ribotoxin coding genes acquired from plants

Ribosome inactivating proteins (RIPs) are RNA N-glycosidases that depurinate a specific adenine residue in the conserved sarcin/ricin loop of 28S rRNA. These enzymes are widely distributed among plants and bacteria. Previously, we have described for the first time RIP genes in mosquitoes belonging to the Culicidae family. We showed that these genes are derived from a single event of horizontal gene transfer (HGT) from a prokaryotic donor. Mosquito RIP genes are evolving under purifying selection, strongly suggesting that these toxins have acquired a functional role. In this work, we show the existence of two RIP encoding genes in the genome of the whitefly Bemisia tabaci, a hemiptera species belonging to the Aleyrodidae family distantly related to mosquitoes. Contamination artifacts were ruled out analyzing three independent B. tabaci genome databases. In contrast to mosquito RIPs, whitefly genes harbor introns and according to transcriptomic evidence are transcribed and spliced. Phylogeny and the taxonomic distribution strongly support that whitefly RIP genes are derived from an independent HGT event from a plant source. These results, along with our previous description of RIPs in Diptera, suggest that the acquired genes are functional in these insects and confer some fitness advantage.

The invisible life inside plants: Deciphering the riddles of endophytic bacterial diversity

Endophytic bacteria often promote plant growth and protect their host plant against pathogens, herbivores, and abiotic stresses including drought, increased salinity or pollution. Current agricultural practices are being challenged in terms of climate change and the ever-increasing demand for food. Therefore, the rational exploitation of bacterial endophytes to increase the productivity and resistance of crops appears to be very promising. However, the efficient and larger-scale use of bacterial endophytes for more effective and sustainable agriculture is hindered by very little knowledge on molecular aspects of plant-endophyte interactions and mechanisms driving bacterial communities in planta. In addition, since most of the information on bacterial endophytes has been obtained through culture-dependent techniques, endophytic bacterial diversity and its full biotechnological potential still remain highly unexplored. In this study, we discuss the diversity and role of endophytic populations as well as complex interactions that the endophytes have with the plant and vice versa, including the interactions leading to plant colonization. A description of biotic and abiotic factors influencing endophytic bacterial communities is provided, along with a summary of different methodologies suitable for determining the diversity of bacterial endophytes, mechanisms governing the assembly and structure of bacterial communities in the endosphere, and potential biotechnological applications of endophytes in the future.

Testing the “(Neo-)Darwinian” Principles against Reticulate Evolution: How Variation, Adaptation, Heredity and Fitness, Constraints and Affordances, Speciation, and Extinction Surpass Organisms and Species

Variation, adaptation, heredity and fitness, constraints and affordances, speciation, and extinction form the building blocks of the (Neo-)Darwinian research program, and several of these have been called “Darwinian principles”. Here, we suggest that caution should be taken in calling these principles Darwinian because of the important role played by reticulate evolutionary mechanisms and processes in also bringing about these phenomena. Reticulate mechanisms and processes include symbiosis, symbiogenesis, lateral gene transfer, infective heredity mediated by genetic and organismal mobility, and hybridization. Because the “Darwinian principles” are brought about by both vertical and reticulate evolutionary mechanisms and processes, they should be understood as foundational for a more pluralistic theory of evolution, one that surpasses the classic scope of the Modern and the Neo-Darwinian Synthesis. Reticulate evolution moreover demonstrates that what conventional (Neo-)Darwinian theories treat as intra-species features of evolution frequently involve reticulate interactions between organisms from very different taxonomic categories. Variation, adaptation, heredity and fitness, constraints and affordances, speciation, and extinction therefore cannot be understood as “traits” or “properties” of genes, organisms, species, or ecosystems because the phenomena are irreducible to specific units and levels of an evolutionary hierarchy. Instead, these general principles of evolution need to be understood as common goods that come about through interactions between different units and levels of evolutionary hierarchies, and they are exherent rather than inherent properties of individuals.

Desiccation does not drastically increase the accessibility of exogenous DNA to nuclear genomes: evidence from the frequency of endosymbiotic DNA transfer

Background

Although horizontal gene transfer (HGT) is a widely accepted force in the evolution of prokaryotic genomes, its role in the evolution of eukaryotic genomes remains hotly debated. Some bdelloid rotifers that are resistant to extreme desiccation and radiation undergo a very high level of HGT, whereas in another desiccation-resistant invertebrate, the tardigrade, the pattern does not exist. Overall, the DNA double-strand breaks (DSBs) induced by prolonged desiccation have been postulated to open a gateway to the nuclear genome for exogenous DNA integration and thus to facilitate the HGT process, thereby enhancing the rate of endosymbiotic DNA transfer (EDT).

Results

We first surveyed the abundance of nuclear mitochondrial DNAs (NUMTs) and nuclear plastid DNAs (NUPTs) in five eukaryotes that are highly resistant to desiccation: the bdelloid rotifers Adineta vaga and Adineta ricciae, the tardigrade Ramazzottius varieornatus, and the resurrection plants Dorcoceras hygrometricum and Selaginella tamariscina. Excessive NUMTs or NUPTs were not detected. Furthermore, we compared 24 groups of desiccation-tolerant organisms with their relatively less desiccation-tolerant relatives but did not find a significant difference in NUMT/NUPT contents.

Conclusions

Desiccation may induce DSBs, but it is unlikely to dramatically increase the frequency of exogenous sequence integration in most eukaryotes. The capture of exogenous DNA sequences is possible only when DSBs are repaired through a subtype of non-homologous end joining, named alternative end joining (alt-EJ). Due to the deleterious effects of the resulting insertion mutations, alt-EJ is less frequently initiated than other mechanisms.

Exploration of space to achieve scientific breakthroughs

Living organisms adapt to changing environments using their amazing flexibility to remodel themselves by a process called evolution. Environmental stress causes selective pressure and is associated with genetic and phenotypic shifts for better modifications, maintenance, and functioning of organismal systems. The natural evolution process can be used in complement to rational strain engineering for the development of desired traits or phenotypes as well as for the production of novel biomaterials through the imposition of one or more selective pressures. Space provides a unique environment of stressors (e.g., weightlessness and high radiation) that organisms have never experienced on Earth. Cells in the outer space reorganize and develop or activate a range of molecular responses that lead to changes in cellular properties. Exposure of cells to the outer space will lead to the development of novel variants more efficiently than on Earth. For instance, natural crop varieties can be generated with higher nutrition value, yield, and improved features, such as resistance against high and low temperatures, salt stress, and microbial and pest attacks. The review summarizes the literature on the parameters of outer space that affect the growth and behavior of cells and organisms as well as complex colloidal systems. We illustrate an understanding of gravity-related basic biological mechanisms and enlighten the possibility to explore the outer space environment for application-oriented aspects. This will stimulate biological research in the pursuit of innovative approaches for the future of agriculture and health on Earth.

Characterization of horizontally acquired ribotoxin encoding genes and their transcripts in Aedes aegypti

Ribosome Inactivating Proteins (RIPs) are RNA N-glycosidases that depurinate a specific adenine residue in the conserved sarcin/ricin loop of the 28S rRNA. The occurrence of RIP genes has been described in a wide range of plant taxa, as well as in several species of bacteria and fungi. A remarkable case is the presence of these genes in metazoans belonging to the Culicinae subfamily. We reported that these genes are derived from a single horizontal gene transfer event, most likely from a bacterial donor species. Moreover, we have shown evidence that mosquito RIP genes are evolving under purifying selection, suggesting that these toxins have acquired a functional role in these organisms. In the present work, we characterized the intra-specific sequence variability of Aedes aegypti RIP genes (RIPAe1, RIPAe2, and RIPAe3) and tested their expression at the mRNA level. Our results show that RIPAe2 and RIPAe3 are transcribed and polyadenylated, and their expression levels are modulated across the developmental stages. Varibility among genes was observed, including the existence of null alleles for RIPAe1 and RIPAe2, with variants showing partial deletions. These results further support the existence of a physiological function for these foreign genes in mosquitoes. The possible nature of this functionality is discussed.

Phytoplankton pangenome reveals extensive prokaryotic horizontal gene transfer of diverse functions

The extent and role of horizontal gene transfer (HGT) in phytoplankton and, more broadly, eukaryotic evolution remain controversial topics. Recent studies substantiate the importance of HGT in modifying or expanding functions such as metal or reactive species detoxification and buttressing halotolerance. Yet, the potential of HGT to significantly alter the fate of species in a major eukaryotic assemblage remains to be established. We provide such an example for the ecologically important lineages encompassed by cryptophytes, rhizarians, alveolates, stramenopiles, and haptophytes ("CRASH" taxa). We describe robust evidence of prokaryotic HGTs in these taxa affecting functions such as polysaccharide biosynthesis. Numbers of HGTs range from 0.16 to 1.44% of CRASH species gene inventories, comparable to the ca. 1% prokaryote-derived HGTs found in the genomes of extremophilic red algae. Our results substantially expand the impact of HGT in eukaryotes and define a set of general principles for prokaryotic gene fixation in phytoplankton genomes.

Horizontal Transfer of Promiscuous Activity from Nonphotosynthetic Bacteria Contributed to Evolution of Chlorophyll Degradation Pathway

The relationship between enzymes and substrates does not perfectly match the "lock and key" model, because enzymes act on molecules other than their true substrate in different catalytic reactions. Such biologically nonfunctional reactions are called "promiscuous activities." Promiscuous activities are apparently useless, but they can be an important starting point for enzyme evolution. It has been hypothesized that enzymes with low promiscuous activity will show enhanced promiscuous activity under selection pressure and become new specialists through gene duplication. Although this is the prevailing scenario, there are two major problems: 1) it would not apply to prokaryotes because horizontal gene transfer is more significant than gene duplication and 2) there is no direct evidence that promiscuous activity is low without selection pressure. We propose a new scenario including various levels of promiscuous activity throughout a clade and horizontal gene transfer. STAY-GREEN (SGR), a chlorophyll a-Mg dechelating enzyme, has homologous genes in bacteria lacking chlorophyll. We found that some bacterial SGR homologs have much higher Mg-dechelating activities than those of green plant SGRs, while others have no activity, indicating that the level of promiscuous activity varies. A phylogenetic analysis suggests that a bacterial SGR homolog with high dechelating activity was horizontally transferred to a photosynthetic eukaryote. Some SGR homologs acted on various chlorophyll molecules that are not used as substrates by green plant SGRs, indicating that SGR acquired substrate specificity after transfer to eukaryotes. We propose that horizontal transfer of high promiscuous activity is one process of new enzyme acquisition.

Plant Evolution: Assembling Land Plants

Traditional evolutionary scenarios posit that land plants emerged from land plant-like relatives, the charophytes. New phylogenies suggest a closer affinity to simpler pond scum relatives, and evidence the gradual assembly of the land plant genome, revealing a phenotypic simplification from the complex ancestors envisaged by traditional scenarios.

On the evolutionary significance of horizontal gene transfers in plants

Horizontal gene transfer (HGT) has long been seen as a crucial process in the evolution of prokaryotic species, but until recently it was thought to have little, if any, effect on the evolution of eukaryotic life forms. Detecting and describing HGT events in eukaryotes is difficult, making this phenomenon at times controversial. However, modern advances in genomics and bioinformatics have radically altered our view of HGT in eukaryotes, especially in plants. It now appears that HGT to and from plant lineages is more common than previously suspected. Importantly, the transfer of functional nuclear genes with adaptive significance has been reported in numerous taxa. Here we review several recent studies that have found evidence of the horizontal transfer of nuclear genes, and argue that HGT has undoubtedly had profound impacts on plant evolution as a whole.

Genomic mosaicism underlies the adaptation of marine Synechococcus ecotypes to distinct oceanic iron niches

Phytoplankton are limited by iron (Fe) in ~40% of the world's oceans including high-nutrient low-chlorophyll (HNLC) regions. While low-Fe adaptation has been well-studied in large eukaryotic diatoms, less is known for small, prokaryotic marine picocyanobacteria. This study reveals key physiological and genomic differences underlying Fe adaptation in marine picocyanobacteria. HNLC ecotype CRD1 strains have greater physiological tolerance to low Fe congruent with their expanded repertoire of Fe transporter, storage and regulatory genes compared to other ecotypes. From metagenomic analysis, genes encoding ferritin, flavodoxin, Fe transporters and siderophore uptake genes were more abundant in low-Fe waters, mirroring paradigms of low-Fe adaptation in diatoms. Distinct Fe-related gene repertories of HNLC ecotypes CRD1 and CRD2 also highlight how coexisting ecotypes have evolved independent approaches to life in low-Fe habitats. Synechococcus and Prochlorococcus HNLC ecotypes likewise exhibit independent, genome-wide reductions of predicted Fe-requiring genes. HNLC ecotype CRD1 interestingly was most similar to coastal ecotype I in Fe physiology and Fe-related gene content, suggesting populations from these different biomes experience similar Fe-selective conditions. This work supports an improved perspective that phytoplankton are shaped by more nuanced Fe niches in the oceans than previously implied from mostly binary comparisons of low- versus high-Fe habitats and populations.

PhySpeTree: an automated pipeline for reconstructing phylogenetic species trees

Background

Phylogenetic species trees are widely used in inferring evolutionary relationships. Existing software and algorithms mainly focus on phylogenetic inference. However, less attention has been paid to intermediate steps, such as processing extremely large sequences and preparing configure files to connect multiple software. When the species number is large, the intermediate steps become a bottleneck that may seriously affect the efficiency of tree building.

Results

Here, we present an easy-to-use pipeline named PhySpeTree to facilitate the reconstruction of species trees across bacterial, archaeal, and eukaryotic organisms. Users need only to input the abbreviations of species names; PhySpeTree prepares complex configure files for different software, then automatically downloads genomic data, cleans sequences, and builds trees. PhySpeTree allows users to perform critical steps such as sequence alignment and tree construction by adjusting advanced options. PhySpeTree provides two parallel pipelines based on concatenated highly conserved proteins and small subunit ribosomal RNA sequences, respectively. Accessory modules, such as those for inserting new species, generating visualization configurations, and combining trees, are distributed along with PhySpeTree.

Conclusions

Together with accessory modules, PhySpeTree significantly simplifies tree reconstruction. PhySpeTree is implemented in Python running on modern operating systems (Linux, macOS, and Windows). The source code is freely available with detailed documentation (https://github.com/yangfangs/physpetools).

Evolutionary entanglement of mobile genetic elements and host defence systems: guns for hire

All cellular life forms are afflicted by diverse genetic parasites, including viruses and other types of mobile genetic elements (MGEs), and have evolved multiple, diverse defence systems that protect them from MGE assault via different mechanisms. Here, we provide our perspectives on how recent evidence points to tight evolutionary connections between MGEs and defence systems that reach far beyond the proverbial arms race. Defence systems incur a fitness cost for the hosts; therefore, at least in prokaryotes, horizontal mobility of defence systems, mediated primarily by MGEs, is essential for their persistence. Moreover, defence systems themselves possess certain features of selfish elements. Common components of MGEs, such as site-specific nucleases, are 'guns for hire' that can also function as parts of defence mechanisms and are often shuttled between MGEs and defence systems. Thus, evolutionary and molecular factors converge to mould the multifaceted, inextricable connection between MGEs and anti-MGE defence systems.

Putative group I introns in the eukaryote nuclear internal transcribed spacers

Group I introns are mobile genetic elements that interrupt genes encoding proteins and RNAs. In the rRNA operon, introns can insert in the small subunit (SSU) and large subunit (LSU) of a wide variety of protists and various prokaryotes, but they were never found in the ITS region. In this study, unusually long ITS regions of fungi and closely related unicellular organisms (Polychytrium aggregatum, Mitosporidium daphniae, Amoeboaphelidium occidentale and Nuclearia simplex) were analysed. While the insertion of repeats is responsible for long ITS in other eukaryotes, the increased size of the sequences analysed herein seems rather due to the presence of introns in ITS-1 or ITS-2. The identified insertions can be folded in secondary structures according to group I intron models, and they cluster within introns in conserved core-based phylogeny. In addition, for Mitosporidium, Amoeboaphelidium and Nuclearia, more conventional ITS-2 structures can be deduced once spacer introns are removed. Sequences of five shark species were also analysed for their structure and included in phylogeny because of unpublished work reporting introns in their ITS, obtaining congruent results. Overall, the data presented herein indicate that spacer regions may contain introns.

Phylogeny and Evolution of RNA 3'-Nucleotidyltransferases in Bacteria

The tRNA nucleotidyltransferases and poly(A) polymerases belong to a superfamily of nucleotidyltransferases. The amino acid sequences of a number of bacterial tRNA nucleotidyltransferases and poly(A) polymerases have been used to construct a rooted, neighbor-joining phylogenetic tree. Using information gleaned from that analysis, along with data from the rRNA-based phylogenetic tree, structural data available on a number of members of the superfamily and other biochemical information on the superfamily, it is possible to suggest a scheme for the evolution of the bacterial tRNA nucleotidyltransferases and poly(A) polymerases from ancestral species. Elements of that scheme are discussed along with questions arising from the scheme which can be explored experimentally.

Adapting for life in the extreme

Red algae have adapted to extreme environments by acquiring genes from bacteria and archaea.

Protist Evolution: Stealing Genes to Gut It Out

The amount and evolutionary impact of horizontal gene transfer in eukaryotes remain contentious issues. A new phylogenomic study suggests that gene transfer from prokaryotes has contributed significantly to the adaptation and metabolic evolution of Blastocystis, the most widespread human gut eukaryotic parasite.

The genomes of polyextremophilic cyanidiales contain 1% horizontally transferred genes with diverse adaptive functions

The role and extent of horizontal gene transfer (HGT) in eukaryotes are hotly disputed topics that impact our understanding of the origin of metabolic processes and the role of organelles in cellular evolution. We addressed this issue by analyzing 10 novel Cyanidiales genomes and determined that 1% of their gene inventory is HGT-derived. Numerous HGT candidates share a close phylogenetic relationship with prokaryotes that live in similar habitats as the Cyanidiales and encode functions related to polyextremophily. HGT candidates differ from native genes in GC-content, number of splice sites, and gene expression. HGT candidates are more prone to loss, which may explain the absence of a eukaryotic pan-genome. Therefore, the lack of a pan-genome and cumulative effects fail to provide substantive arguments against our hypothesis of recurring HGT followed by differential loss in eukaryotes. The maintenance of 1% HGTs, even under selection for genome reduction, underlines the importance of non-endosymbiosis related foreign gene acquisition.

Tree reconciliation combined with subsampling improves large scale inference of orthologous group hierarchies

BACKGROUND

An orthologous group (OG) comprises a set of orthologous and paralogous genes that share a last common ancestor (LCA). OGs are defined with respect to a chosen taxonomic level, which delimits the position of the LCA in time to a specified speciation event. A hierarchy of OGs expands on this notion, connecting more general OGs, distant in time, to more recent, fine-grained OGs, thereby spanning multiple levels of the tree of life. Large scale inference of OG hierarchies with independently computed taxonomic levels can suffer from inconsistencies between successive levels, such as the position in time of a duplication event. This can be due to confounding genetic signal or algorithmic limitations. Importantly, inconsistencies limit the potential use of OGs for functional annotation and third-party applications.

RESULTS

Here we present a new methodology to ensure hierarchical consistency of OGs across taxonomic levels. To resolve an inconsistency, we subsample the protein space of the OG members and perform gene tree-species tree reconciliation for each sampling. Differently from previous approaches, by subsampling the protein space, we avoid the notoriously difficult task of accurately building and reconciling very large phylogenies. We implement the method into a high-throughput pipeline and apply it to the eggNOG database. We use independent protein domain definitions to validate its performance.

CONCLUSION

The presented consistency pipeline shows that, contrary to previous limitations, tree reconciliation can be a useful instrument for the construction of OG hierarchies. The key lies in the combination of sampling smaller trees and aggregating their reconciliations for robustness. Results show comparable or greater performance to previous pipelines. The code is available on Github at: https://github.com/meringlab/og_consistency_pipeline .

Pan-genome analyses of model fungal species

The concept of the species ‘pan-genome’, the union of ‘core’ conserved genes and all ‘accessory’ non-conserved genes across all strains of a species, was first proposed in prokaryotes to account for intraspecific variability. Species pan-genomes have been extensively studied in prokaryotes, but evidence of species pan-genomes has also been demonstrated in eukaryotes such as plants and fungi. Using a previously published methodology based on sequence homology and conserved microsynteny, in addition to bespoke pipelines, we have investigated the pan-genomes of four model fungal species: Saccharomyces cerevisiae, Candida albicans, Cryptococcus neoformans var. grubii and Aspergillus fumigatus. Between 80 and 90 % of gene models per strain in each of these species are core genes that are highly conserved across all strains of that species, many of which are involved in housekeeping and conserved survival processes. In many of these species, the remaining ‘accessory’ gene models are clustered within subterminal regions and may be involved in pathogenesis and antimicrobial resistance. Analysis of the ancestry of species core and accessory genomes suggests that fungal pan-genomes evolve by strain-level innovations such as gene duplication as opposed to wide-scale horizontal gene transfer. Our findings lend further supporting evidence to the existence of species pan-genomes in eukaryote taxa.

Genomes of trombidid mites reveal novel predicted allergens and laterally transferred genes associated with secondary metabolism

Background

Trombidid mites have a unique life cycle in which only the larval stage is ectoparasitic. In the superfamily Trombiculoidea ("chiggers"), the larvae feed preferentially on vertebrates, including humans. Species in the genus Leptotrombidium are vectors of a potentially fatal bacterial infection, scrub typhus, that affects 1 million people annually. Moreover, chiggers can cause pruritic dermatitis (trombiculiasis) in humans and domesticated animals. In the Trombidioidea (velvet mites), the larvae feed on other arthropods and are potential biological control agents for agricultural pests. Here, we present the first trombidid mites genomes, obtained both for a chigger, Leptotrombidium deliense, and for a velvet mite, Dinothrombium tinctorium.

Results

Sequencing was performed using Illumina technology. A 180 Mb draft assembly for D. tinctorium was generated from two paired-end and one mate-pair library using a single adult specimen. For L. deliense, a lower-coverage draft assembly (117 Mb) was obtained using pooled, engorged larvae with a single paired-end library. Remarkably, both genomes exhibited evidence of ancient lateral gene transfer from soil-derived bacteria or fungi. The transferred genes confer functions that are rare in animals, including terpene and carotenoid synthesis. Thirty-seven allergenic protein families were predicted in the L. deliense genome, of which nine were unique. Preliminary proteomic analyses identified several of these putative allergens in larvae.

Conclusions

Trombidid mite genomes appear to be more dynamic than those of other acariform mites. A priority for future research is to determine the biological function of terpene synthesis in this taxon and its potential for exploitation in disease control.

Failure to Recover Major Events of Gene Flux in Real Biological Data Due to Method Misapplication

In prokaryotes, known mechanisms of lateral gene transfer (transformation, transduction, conjugation, and gene transfer agents) generate new combinations of genes among chromosomes during evolution. In eukaryotes, whose host lineage is descended from archaea, lateral gene transfer from organelles to the nucleus occurs at endosymbiotic events. Recent genome analyses studying gene distributions have uncovered evidence for sporadic, discontinuous events of gene transfer from bacteria to archaea during evolution. Other studies have used traditional models designed to investigate gene family size evolution (Count) to support claims that gene transfer to archaea was continuous during evolution, rather than involving occasional periodic mass gene influx events. Here, we show that the methodology used in analyses favoring continuous gene transfers to archaea was misapplied in other studies and does not recover known events of single simultaneous origin for many genes followed by differential loss in real data: plastid genomes. Using the same software and the same settings, we reanalyzed presence/absence pattern data for proteins encoded in plastid genomes and for eukaryotic protein families acquired from plastids. Contrary to expectations under a plastid origin model, we found that the methodology employed inferred that gene acquisitions occurred uniformly across the plant tree. Sometimes as many as nine different acquisitions by plastid DNA were inferred for the same protein family. That is, the methodology that recovered gradual and continuous lateral gene transfer among lineages for archaea obtains the same result for plastids, even though it is known that massive gains followed by gradual differential loss is the true evolutionary process that generated plastid gene distribution data. Our findings caution against the use of models designed to study gene family size evolution for investigating gene transfer processes, especially when transfers involving more than one gene per event are possible.

The last universal common ancestor between ancient Earth chemistry and the onset of genetics

All known life forms trace back to a last universal common ancestor (LUCA) that witnessed the onset of Darwinian evolution. One can ask questions about LUCA in various ways, the most common way being to look for traits that are common to all cells, like ribosomes or the genetic code. With the availability of genomes, we can, however, also ask what genes are ancient by virtue of their phylogeny rather than by virtue of being universal. That approach, undertaken recently, leads to a different view of LUCA than we have had in the past, one that fits well with the harsh geochemical setting of early Earth and resembles the biology of prokaryotes that today inhabit the Earth's crust.

A Novel Eukaryotic Denitrification Pathway in Foraminifera

Benthic foraminifera are unicellular eukaryotes inhabiting sediments of aquatic environments. Several species were shown to store and use nitrate for complete denitrification, a unique energy metabolism among eukaryotes. The population of benthic foraminifera reaches high densities in oxygen-depleted marine habitats, where they play a key role in the marine nitrogen cycle. However, the mechanisms of denitrification in foraminifera are still unknown, and the possibility of a contribution of associated bacteria is debated. Here, we present evidence for a novel eukaryotic denitrification pathway that is encoded in foraminiferal genomes. Large-scale genome and transcriptomes analyses reveal the presence of a denitrification pathway in foraminifera species of the genus Globobulimina. This includes the enzymes nitrite reductase (NirK) and nitric oxide reductase (Nor) as well as a wide range of nitrate transporters (Nrt). A phylogenetic reconstruction of the enzymes' evolutionary history uncovers evidence for an ancient acquisition of the foraminiferal denitrification pathway from prokaryotes. We propose a model for denitrification in foraminifera, where a common electron transport chain is used for anaerobic and aerobic respiration. The evolution of hybrid respiration in foraminifera likely contributed to their ecological success, which is well documented in palaeontological records since the Cambrian period.

A massive incorporation of microbial genes into the genome of Tetranychus urticae, a polyphagous arthropod herbivore

A number of horizontal gene transfers (HGTs) have been identified in the spider mite Tetranychus urticae, a chelicerate herbivore. However, the genome of this mite species has at present not been thoroughly mined for the presence of HGT genes. Here, we performed a systematic screen for HGT genes in the T. urticae genome using the h-index metric. Our results not only validated previously identified HGT genes but also uncovered 25 novel HGT genes. In addition to HGT genes with a predicted biochemical function in carbohydrate, lipid and folate metabolism, we also identified the horizontal transfer of a ketopantoate hydroxymethyltransferase and a pantoate β-alanine ligase gene. In plants and bacteria, both genes are essential for vitamin B5 biosynthesis and their presence in the mite genome strongly suggests that spider mites, similar to Bemisia tabaci and nematodes, can synthesize their own vitamin B5. We further show that HGT genes were physically embedded within the mite genome and were expressed in different life stages. By screening chelicerate genomes and transcriptomes, we were able to estimate the evolutionary histories of these HGTs during chelicerate evolution. Our study suggests that HGT has made a significant and underestimated impact on the metabolic repertoire of plant-feeding spider mites.

Comparative analyses of whole-genome protein sequences from multiple organisms

Phylogenies based on entire genomes are a powerful tool for reconstructing the Tree of Life. Several methods have been proposed, most of which employ an alignment-free strategy. Average sequence similarity methods are different than most other whole-genome methods, because they are based on local alignments. However, previous average similarity methods fail to reconstruct a correct phylogeny when compared against other whole-genome trees. In this study, we developed a novel average sequence similarity method. Our method correctly reconstructs the phylogenetic tree of in silico evolved E. coli proteomes. We applied the method to reconstruct a whole-proteome phylogeny of 1,087 species from all three domains of life, Bacteria, Archaea, and Eucarya. Our tree was automatically reconstructed without any human decisions, such as the selection of organisms. The tree exhibits a concentric circle-like structure, indicating that all the organisms have similar total branch lengths from their common ancestor. Branching patterns of the members of each phylum of Bacteria and Archaea are largely consistent with previous reports. The topologies are largely consistent with those reconstructed by other methods. These results strongly suggest that this approach has sufficient taxonomic resolution and reliability to infer phylogeny, from phylum to strain, of a wide range of organisms.