Large-scale investigation of species-specific orphan genes in the human gut microbiome elucidates their evolutionary origins.

Species-specific genes, also known as orphans, are ubiquitous across life's domains. In prokaryotes, species-specific orphan genes (SSOGs) are mostly thought to originate in external elements such as viruses followed by horizontal gene transfer, whereas the scenario of native origination, through rapid divergence or de novo, is mostly dismissed. However, quantitative evidence supporting either scenario is lacking. Here, we systematically analyzed genomes from 4644 human gut microbiome species and identified more than 600,000 unique SSOGs, representing an average of 2.6% of a given species' pangenome. These sequences are mostly rare within each species yet show signs of purifying selection. Overall, SSOGs use optimal codons less frequently, and their proteins are more disordered than those of conserved genes (i.e., non-SSOGs). Importantly, across species, the GC content of SSOGs closely matches that of conserved ones. In contrast, the ∼5% of SSOGs that share similarity to known viral sequences have distinct characteristics, including lower GC content. Thus, SSOGs with similarity to viruses differ from the remaining SSOGs, contrasting an external origination scenario for most of them. By examining the orthologous genomic region in closely related species, we show that a small subset of SSOGs likely evolved natively de novo and find that these genes also differ in their properties from the remaining SSOGs. Our results challenge the notion that external elements are the dominant source of prokaryotic genetic novelty and will enable future studies into the biological role and relevance of species-specific genes in the human gut.

Active in vivo translocation of the Methanosarcina mazei Gö1 Casposon.

Casposons are transposable elements containing the CRISPR associated gene Cas1solo. Identified in many archaeal genomes, casposons are discussed as the origin of CRISPR-Cas systems due to their proposed Cas1solo-dependent translocation. However, apart from bioinformatic approaches and the demonstration of Cas1solo integrase and endonuclease activity in vitro, casposon transposition has not yet been shown in vivo. Here, we report on active casposon translocations in Methanosarcina mazei Gö1 using two independent experimental approaches. First, mini-casposons, consisting of a R6Kγ origin and two antibiotic resistance cassettes, flanked by target site duplications (TSDs) and terminal inverted repeats (TIRs), were generated, and shown to actively translocate from a suicide plasmid and integrate into the chromosomal MetMaz-C1 TSD IS1a. Second, casposon excision activity was confirmed in a long-term evolution experiment using a Cas1solo overexpression strain in comparison to an empty vector control under four different treatments (native, high temperature, high salt, mitomycin C) to study stress-induced translocation. Analysis of genomic DNA using a nested qPCR approach provided clear evidence of casposon activity in single cells and revealed significantly different casposon excision frequencies between treatments and strains. Our results, providing the first experimental evidence for in vivo casposon activity are summarized in a modified hypothetical translocation model.

Tackling the Pangenome Dilemma Requires the Concerted Analysis of Multiple Population Genetic Processes.

The pangenome is the set of all genes present in a prokaryotic population. Most pangenomes contain many accessory genes of low and intermediate frequencies. Different population genetics processes contribute to the shape of these pangenomes, namely selection and fitness-independent processes such as gene transfer, gene loss, and migration. However, their relative importance is unknown and highly debated. Here, we argue that the debate around prokaryotic pangenomes arose due to the imprecise application of population genetics models. Most importantly, two different processes of horizontal gene transfer act on prokaryotic populations, which are frequently confused, despite their fundamentally different behavior. Genes acquired from distantly related organisms (termed here acquiring gene transfer) are most comparable to mutation in nucleotide sequences. In contrast, gene gain within the population (termed here spreading gene transfer) has an effect on gene frequencies that is identical to the effect of positive selection on single genes. We thus show that selection and fitness-independent population genetic processes affecting pangenomes are indistinguishable at the level of single gene dynamics. Nevertheless, population genetics processes are fundamentally different when considering the joint distribution of all accessory genes across individuals of a population. We propose that, to understand to which degree the different processes shaped pangenome diversity, the development of comprehensive models and simulation tools is mandatory. Furthermore, we need to identify summary statistics and measurable features that can distinguish between the processes, where considering the joint distribution of accessory genes across individuals of a population will be particularly relevant.

Highly efficient reduction of ammonia emissions from livestock waste by the synergy of novel manure acidification and inhibition of ureolytic bacteria.

The global livestock system is one of the largest sources of ammonia emissions and there is an urgent need for ammonia mitigation. Here, we designed and constructed a novel strategy to abate ammonia emissions via livestock manure acidification based on a synthetic lactic acid bacteria community (LAB SynCom). The LAB SynCom possessed a wide carbon source spectrum and pH profile, high adaptability to the manure environment, and a high capability of generating lactic acid. The mitigation strategy was optimized based on the test and performance by adjusting the LAB SynCom inoculation ratio and the adding frequency of carbon source, which contributed to a total ammonia reduction efficiency of 95.5 %. Furthermore, 16S rDNA amplicon sequencing analysis revealed that the LAB SynCom treatment reshaped the manure microbial community structure. Importantly, 22 manure ureolytic microbial genera and urea hydrolysis were notably inhibited by the LAB SynCom treatment during the treatment process. These findings provide new insight into manure acidification that the conversion from ammonia to ammonium ions and the inhibition of ureolytic bacteria exerted a synergistic effect on ammonia mitigation. This work systematically developed a novel strategy to mitigate ammonia emissions from livestock waste, which is a crucial step forward from traditional manure acidification to novel and environmental-friendly acidification.

Arabidopsis latent virus 1, a comovirus widely spread in Arabidopsis thaliana collections.

Transcriptome studies of Illumina RNA-Seq datasets of different Arabidopsis thaliana natural accessions and T-DNA mutants revealed the presence of two virus-like RNA sequences which showed the typical two-segmented genome characteristics of a comovirus. This comovirus did not induce any visible symptoms in infected A. thaliana plants cultivated under standard laboratory conditions. Hence it was named Arabidopsis latent virus 1 (ArLV1). Virus infectivity in A. thaliana plants was confirmed by quantitative reverse transcription polymerase chain reaction, transmission electron microscopy and mechanical inoculation. Arabidopsis latent virus 1 can also mechanically infect Nicotiana benthamiana, causing distinct mosaic symptoms. A bioinformatics investigation of A. thaliana RNA-Seq repositories, including nearly 6500 Sequence Read Archives (SRAs) in the NCBI SRA database, revealed the presence of ArLV1 in 25% of all archived natural A. thaliana accessions and in 8.5% of all analyzed SRAs. Arabidopsis latent virus 1 could also be detected in A. thaliana plants collected from the wild. Arabidopsis latent virus 1 is highly seed-transmissible with up to 40% incidence on the progeny derived from infected A. thaliana plants. This has probably led to a worldwide distribution in the model plant A. thaliana with as yet unknown effects on plant performance in a substantial number of studies.

Fixation dynamics of beneficial alleles in prokaryotic polyploid chromosomes and plasmids.

Theoretical population genetics has been mostly developed for sexually reproducing diploid and for monoploid (haploid) organisms, focusing on eukaryotes. The evolution of bacteria and archaea is often studied by models for the allele dynamics in monoploid populations. However, many prokaryotic organisms harbor multicopy replicons-chromosomes and plasmids-and theory for the allele dynamics in populations of polyploid prokaryotes remains lacking. Here, we present a population genetics model for replicons with multiple copies in the cell. Using this model, we characterize the fixation process of a dominant beneficial mutation at 2 levels: the phenotype and the genotype. Our results show that depending on the mode of replication and segregation, the fixation of the mutant phenotype may precede genotypic fixation by many generations; we term this time interval the heterozygosity window. We furthermore derive concise analytical expressions for the occurrence and length of the heterozygosity window, showing that it emerges if the copy number is high and selection strong. Within the heterozygosity window, the population is phenotypically adapted, while both alleles persist in the population. Replicon ploidy thus allows for the maintenance of genetic variation following phenotypic adaptation and consequently for reversibility in adaptation to fluctuating environmental conditions.

Pangenome Evolution in Environmentally Transmitted Symbionts of Deep-Sea Mussels Is Governed by Vertical Inheritance.

Microbial pangenomes vary across species; their size and structure are determined by genetic diversity within the population and by gene loss and horizontal gene transfer (HGT). Many bacteria are associated with eukaryotic hosts where the host colonization dynamics may impact bacterial genome evolution. Host-associated lifestyle has been recognized as a barrier to HGT in parentally transmitted bacteria. However, pangenome evolution of environmentally acquired symbionts remains understudied, often due to limitations in symbiont cultivation. Using high-resolution metagenomics, here we study pangenome evolution of two co-occurring endosymbionts inhabiting Bathymodiolus brooksi mussels from a single cold seep. The symbionts, sulfur-oxidizing (SOX) and methane-oxidizing (MOX) gamma-proteobacteria, are environmentally acquired at an early developmental stage and individual mussels may harbor multiple strains of each symbiont species. We found differences in the accessory gene content of both symbionts across individual mussels, which are reflected by differences in symbiont strain composition. Compared with core genes, accessory genes are enriched in genome plasticity functions. We found no evidence for recent HGT between both symbionts. A comparison between the symbiont pangenomes revealed that the MOX population is less diverged and contains fewer accessory genes, supporting that the MOX association with B. brooksi is more recent in comparison to that of SOX. Our results show that the pangenomes of both symbionts evolved mainly by vertical inheritance. We conclude that genome evolution of environmentally transmitted symbionts that associate with individual hosts over their lifetime is affected by a narrow symbiosis where the frequency of HGT is constrained.

Disentangling the genetic basis of rhizosphere microbiome assembly in tomato.

Microbiomes play a pivotal role in plant growth and health, but the genetic factors involved in microbiome assembly remain largely elusive. Here, we map the molecular features of the rhizosphere microbiome as quantitative traits of a diverse hybrid population of wild and domesticated tomato. Gene content analysis of prioritized tomato quantitative trait loci suggests a genetic basis for differential recruitment of various rhizobacterial lineages, including a Streptomyces-associated 6.31 Mbp region harboring tomato domestication sweeps and encoding, among others, the iron regulator FIT and the water channel aquaporin SlTIP2.3. Within metagenome-assembled genomes of root-associated Streptomyces and Cellvibrio, we identify bacterial genes involved in metabolism of plant polysaccharides, iron, sulfur, trehalose, and vitamins, whose genetic variation associates with specific tomato QTLs. By integrating 'microbiomics' and quantitative plant genetics, we pinpoint putative plant and reciprocal rhizobacterial traits underlying microbiome assembly, thereby providing a first step towards plant-microbiome breeding programs.

Co-transfer of functionally interdependent genes contributes to genome mosaicism in lambdoid phages.

Lambdoid (or Lambda-like) phages are a group of related temperate phages that can infect Escherichia coli and other gut bacteria. A key characteristic of these phages is their mosaic genome structure, which served as the basis for the 'modular genome hypothesis'. Accordingly, lambdoid phages evolve by transferring genomic regions, each of which constitutes a functional unit. Nevertheless, it is unknown which genes are preferentially transferred together and what drives such co-transfer events. Here we aim to characterize genome modularity by studying co-transfer of genes among 95 distantly related lambdoid (pro-)phages. Based on gene content, we observed that the genomes cluster into 12 groups, which are characterized by a highly similar gene content within the groups and highly divergent gene content across groups. Highly similar proteins can occur in genomes of different groups, indicating that they have been transferred. About 26 % of homologous protein clusters in the four known operons (i.e. the early left, early right, immunity and late operon) engage in gene transfer, which affects all operons to a similar extent. We identified pairs of genes that are frequently co-transferred and observed that these pairs tend to be near one another on the genome. We find that frequently co-transferred genes are involved in related functions and highlight interesting examples involving structural proteins, the cI repressor and Cro regulator, proteins interacting with DNA, and membrane-interacting proteins. We conclude that epistatic effects, where the functioning of one protein depends on the presence of another, play an important role in the evolution of the modular structure of these genomes.

Characterization of Blf4, an Archaeal Lytic Virus Targeting a Member of the Methanomicrobiales.

Today, the number of known viruses infecting methanogenic archaea is limited. Here, we report on a novel lytic virus, designated Blf4, and its host strain Methanoculleus bourgensis E02.3, a methanogenic archaeon belonging to the Methanomicrobiales, both isolated from a commercial biogas plant in Germany. The virus consists of an icosahedral head 60 nm in diameter and a long non-contractile tail of 125 nm in length, which is consistent with the new isolate belonging to the Siphoviridae family. Electron microscopy revealed that Blf4 attaches to the vegetative cells of M. bourgensis E02.3 as well as to cellular appendages. Apart from M. bourgensis E02.3, none of the tested Methanoculleus strains were lysed by Blf4, indicating a narrow host range. The complete 37 kb dsDNA genome of Blf4 contains 63 open reading frames (ORFs), all organized in the same transcriptional direction. For most of the ORFs, potential functions were predicted. In addition, the genome of the host M. bourgensis E02.3 was sequenced and assembled, resulting in a 2.6 Mbp draft genome consisting of nine contigs. All genes required for a hydrogenotrophic lifestyle were predicted. A CRISPR/Cas system (type I-U) was identified with six spacers directed against Blf4, indicating that this defense system might not be very efficient in fending off invading Blf4 virus.

Mass burial genomics reveals outbreak of enteric paratyphoid fever in the Late Medieval trade city Lübeck.

Medieval Europe was repeatedly affected by outbreaks of infectious diseases, some of which reached epidemic proportions. A Late Medieval mass burial next to the Heiligen-Geist-Hospital in Lübeck (present-day Germany) contained the skeletal remains of more than 800 individuals who had presumably died from infectious disease. From 92 individuals, we screened the ancient DNA extracts for the presence of pathogens to determine the cause of death. Metagenomic analysis revealed evidence of Salmonella enterica subsp. enterica serovar Paratyphi C, suggesting an outbreak of enteric paratyphoid fever. Three reconstructed S. Paratyphi C genomes showed close similarity to a strain from Norway (1200 CE). Radiocarbon dates placed the disease outbreak in Lübeck between 1270 and 1400 cal CE, with historical records indicating 1367 CE as the most probable year. The deceased were of northern and eastern European descent, confirming Lübeck as an important trading center of the Hanseatic League in the Baltic region.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak.

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

Functional diversity enables multiple symbiont strains to coexist in deep-sea mussels.

Genetic diversity of closely related free-living microorganisms is widespread and underpins ecosystem functioning, but most evolutionary theories predict that it destabilizes intimate mutualisms. Accordingly, strain diversity is assumed to be highly restricted in intracellular bacteria associated with animals. Here, we sequenced metagenomes and metatranscriptomes of 18 Bathymodiolus mussel individuals from four species, covering their known distribution range at deep-sea hydrothermal vents in the Atlantic. We show that as many as 16 strains of intracellular, sulfur-oxidizing symbionts coexist in individual Bathymodiolus mussels. Co-occurring symbiont strains differed extensively in key functions, such as the use of energy and nutrient sources, electron acceptors and viral defence mechanisms. Most strain-specific genes were expressed, highlighting their potential to affect fitness. We show that fine-scale diversity is pervasive in Bathymodiolus sulfur-oxidizing symbionts, and hypothesize that it may be widespread in low-cost symbioses where the environment, rather than the host, feeds the symbionts.

A Phage Protein Aids Bacterial Symbionts in Eukaryote Immune Evasion.

Phages are increasingly recognized as important members of host-associated microbiomes, with a vast genomic diversity. The new frontier is to understand how phages may affect higher order processes, such as in the context of host-microbe interactions. Here, we use marine sponges as a model to investigate the interplay between phages, bacterial symbionts, and eukaryotic hosts. Using viral metagenomics, we find that sponges, although massively filtering seawater, harbor species-specific and even individually unique viral signatures that are taxonomically distinct from other environments. We further discover a symbiont phage-encoded ankyrin-domain-containing protein, which is widely spread in phages of many host-associated contexts including human. We confirm in macrophage infection assays that the ankyrin protein (ANKp) modulates the eukaryotic host immune response against bacteria. We predict that the role of ANKp in nature is to facilitate coexistence in the tripartite interplay between phages, symbionts, and sponges and possibly many other host-microbe associations.

Horizontally transmitted symbiont populations in deep-sea mussels are genetically isolated.

Eukaryotes are habitats for bacterial organisms where the host colonization and dispersal among individual hosts have consequences for the bacterial ecology and evolution. Vertical symbiont transmission leads to geographic isolation of the microbial population and consequently to genetic isolation of microbiotas from individual hosts. In contrast, the extent of geographic and genetic isolation of horizontally transmitted microbiota is poorly characterized. Here we show that chemosynthetic symbionts of individual Bathymodiolus brooksi mussels constitute genetically isolated subpopulations. The reconstruction of core genome-wide strains from high-resolution metagenomes revealed distinct phylogenetic clades. Nucleotide diversity and strain composition vary along the mussel life span and individual hosts show a high degree of genetic isolation. Our results suggest that the uptake of environmental bacteria is a restricted process in B. brooksi, where self-infection of the gill tissue results in serial founder effects during symbiont evolution. We conclude that bacterial colonization dynamics over the host life cycle is thus an important determinant of population structure and genome evolution of horizontally transmitted symbionts.

Rates of Molecular Evolution in a Marine Synechococcus Phage Lineage.

Cyanophages are characterized by vast genomic diversity and the formation of stable ecotypes over time. The evolution of phage diversity includes vertical processes, such as mutation, and horizontal processes, such as recombination and gene transfer. Here, we study the contribution of vertical and horizontal processes to short-term evolution of marine cyanophages. Analyzing time series data of Synechococcus-infecting Myoviridae ecotypes spanning up to 17 years, we found a high contribution of recombination relative to mutation (r/m) in all ecotypes. Additionally, we found a molecular clock of substitution and recombination in one ecotype, RIM8. The estimated RIM8 evolutionary rates are 2.2 genome-wide substitutions per year (1.275 × 10-5 substitutions/site/year) and 29 genome-wide nucleotide alterations due to recombination per year. We found 26 variable protein families, of which only two families have a predicted functional annotation, suggesting that they are auxiliary metabolic genes with bacterial homologs. A comparison of our rate estimates to other phage evolutionary rate estimates in the literature reveals a negative correlation of phage substitution rates with their genome size. A comparison to evolutionary rates in bacterial organisms further shows that phages have high rates of mutation and recombination compared to their bacterial hosts. We conclude that the increased recombination rate in phages likely contributes to their vast genomic diversity.

Characterization of the lytic archaeal virus Drs3 infecting Methanobacterium formicicum.

Viruses are ubiquitous in the biosphere and greatly affect the hosts they infect. It is generally accepted that members of every microbial taxon are susceptible to at least one virus, and a plethora of bacterial viruses are known. In contrast, knowledge of the archaeal virosphere is still limited. Here, a novel lytic archaeal virus is described, designated "Drs3", as well as its host, Methanobacterium formicicum strain Khl10. This hydrogenotrophic methanogenic archaeon and its virus were isolated from the anaerobic digester of an experimental biogas plant in Germany. The tailed virus has an icosahedral head with a diameter of approximately 60 nm and a long non-contractile tail of approximately 230 nm. These structural observations suggest that the new isolate belongs to the family Siphoviridae, but it could not be assigned to an existing genus. Lysis of the host Khl10 was observed 40-44 h after infection. Lysis of the type strain Methanobacterium formicicum DSMZ 1535 was not observed in the presence of Drs3, pointing towards resistance in the type strain or a rather narrow host range of this newly isolated archaeal virus. The complete 37-kb linear dsDNA genome of Drs3 contains 39 open reading frames, only 12 of which show similarity to genes with predicted functions.

Segregational Drift and the Interplay between Plasmid Copy Number and Evolvability.

The ubiquity of plasmids in all prokaryotic phyla and habitats and their ability to transfer between cells marks them as prominent constituents of prokaryotic genomes. Many plasmids are found in their host cell in multiple copies. This leads to an increased mutational supply of plasmid-encoded genes and genetically heterogeneous plasmid genomes. Nonetheless, the segregation of plasmid copies into daughter cells during cell division is considered to occur in the absence of selection on the plasmid alleles. We investigate the implications of random genetic drift of multicopy plasmids during cell division-termed here "segregational drift"-to plasmid evolution. Performing experimental evolution of low- and high-copy non-mobile plasmids in Escherichia coli, we find that the evolutionary rate of multicopy plasmids does not reflect the increased mutational supply expected according to their copy number. In addition, simulated evolution of multicopy plasmid alleles demonstrates that segregational drift leads to increased loss frequency and extended fixation time of plasmid mutations in comparison to haploid chromosomes. Furthermore, an examination of the experimentally evolved hosts reveals a significant impact of the plasmid type on the host chromosome evolution. Our study demonstrates that segregational drift of multicopy plasmids interferes with the retention and fixation of novel plasmid variants. Depending on the selection pressure on newly emerging variants, plasmid genomes may evolve slower than haploid chromosomes, regardless of their higher mutational supply. We suggest that plasmid copy number is an important determinant of plasmid evolvability due to the manifestation of segregational drift.

Recombination Signal in Mycobacterium tuberculosis Stems from Reference-guided Assemblies and Alignment Artefacts.

DNA acquisition via genetic recombination is considered advantageous as it has the potential to bring together beneficial mutations that emerge independently within a population. Furthermore, recombination is considered to contribute to the maintenance of genome stability by purging slightly deleterious mutations. The prevalence of recombination differs among prokaryotic species and depends on the accessibility of DNA transfer mechanisms. An exceptional example is the human pathogen Mycobacterium tuberculosis (MTB) where no clear transfer mechanisms have been so far characterized and the presence of recombination is questioned. Here, we analyze completely assembled MTB genomes in search for evidence of recombination. We find that putative recombination events are enriched in strains reconstructed by reference-guided assembly and in regions with unreliable alignments. In addition, assembly and alignment artefacts introduce phylogenetic signals that are conflicting the established MTB phylogeny. Our results reveal that the so far reported recombination events in MTB are likely to stem from methodological artefacts. We conclude that no reliable signal of recombination is observed in the currently available MTB genomes. Moreover, our study demonstrates the limitations of reference-guided genome assembly for phylogenetic reconstructions. Rigorously de novo assembled genomes of high quality are mandatory in order to distinguish true evolutionary signal from noise, in particular for low diversity species such as MTB.

Bioinformatics Meets Virology: The European Virus Bioinformatics Center's Second Annual Meeting.

The Second Annual Meeting of the European Virus Bioinformatics Center (EVBC), held in Utrecht, Netherlands, focused on computational approaches in virology, with topics including (but not limited to) virus discovery, diagnostics, (meta-)genomics, modeling, epidemiology, molecular structure, evolution, and viral ecology. The goals of the Second Annual Meeting were threefold: (i) to bring together virologists and bioinformaticians from across the academic, industrial, professional, and training sectors to share best practice; (ii) to provide a meaningful and interactive scientific environment to promote discussion and collaboration between students, postdoctoral fellows, and both new and established investigators; (iii) to inspire and suggest new research directions and questions. Approximately 120 researchers from around the world attended the Second Annual Meeting of the EVBC this year, including 15 renowned international speakers. This report presents an overview of new developments and novel research findings that emerged during the meeting.