Sulfate Reduction Drives Elevated Methylmercury Formation in the Water Column of a Eutrophic Freshwater Lake

Mercury (Hg) contamination of aquatic food webs is controlled in part by the formation and accumulation of toxic and bioaccumulative methylmercury (MeHg). MeHg production is mediated by metabolically diverse microorganisms carrying the hgcAB gene pair, while the demethylation reaction is mediated by several biotic and abiotic processes. However, the relative importance of these two processes on MeHg accumulation and the environmental factors that influence them are poorly characterized, especially in eutrophic environments. In this study, both Hg methylation and MeHg demethylation in a eutrophic freshwater lake were linked to ambient MeHg concentrations and hgcA abundance and expression. High methylation rate potentials indicated in situ MeHg formation was a key source of MeHg to the water column, driven by high hgcA abundance and transcription. Molybdate treatment decreased methylation rate potentials, highlighting the importance of sulfate reduction in driving MeHg formation. Sulfate-reducing bacteria accounted for over 50% of the hgcA gene transcription, despite representing less than 10% of the hgcA-carrying microbial community. An arsR-like transcriptional regulator preceded many hgcA sequences; these were transcriptionally active and linked to lower hgcA expression. Overall, this study elucidates the microbial and biogeochemical processes that influence the in situ formation of MeHg in understudied eutrophic freshwater environments.

A transcription factor from the cryptic Escherichia coli Rac prophage controls both phage and host operons

Bacterial genomes are shaped by cryptic prophages, which are viral genomes integrated into the bacterial chromosome. Escherichia coli genomes have 10 prophages on average. Though usually inactive, prophage genes can profoundly impact host cell physiology. Among the phage genes in the E. coli chromosome, there are several putative transcription factors (TFs). These prophage TFs are predicted to control only phage promoters; however, their regulatory functions are not well characterized. The cohabitation of prophages and bacteria has led to conditions under which the majority of prophage genes are unexpressed, at least under normal growth conditions. We characterized a Rac prophage TF, YdaT, expression of which is normally inhibited by Rac TFs and, surprisingly, by the host global regulator OxyR. YdaT, when expressed, leads to a toxic phenotype manifested by drastic cell filamentation and cell death. We determined the binding sites and regulatory action for YdaT, finding two sites within the Rac locus, and one upstream of the host rcsA gene, which codes for the global regulator RcsA. The resulting increase in RcsA strongly impacts the bacterial RcsA/B regulon, which includes operons related to motility, capsule biosynthesis, colanic acid production, biofilm formation, and cell division. Our results provide novel insights into the host's genetic network, which appears to integrate YdaT in a complex manner, favoring its maintenance in the silenced state. The fact that the potentially toxic YdaT locus remains unmutated suggests its importance and potential benefits for the host, which may appear under stress conditions that are not yet known.

Unveiling intraspecific diversity and evolutionary dynamics of the foodborne pathogen Bacillus paranthracis through high-quality pan-genome analysis

Understanding the evolutionary dynamics of foodborne pathogens throughout host-associated habitats is of utmost importance. Bacterial pan-genomes, as dynamic entities, are strongly influenced by ecological lifestyles. As a phenotypically diverse species in the Bacillus cereus group, Bacillus paranthracis is recognized as an emerging foodborne pathogen and a probiotic simultaneously. This poorly understood species is a suitable study model for adaptive pan-genome evolution. In this study, we determined the biogeographic distribution, abundance, genetic diversity, and genotypic profiles of key genetic elements of B. paranthracis. Metagenomic read recruitment analyses demonstrated that B. paranthracis members are globally distributed and abundant in host-associated habitats. A high-quality pan-genome of B. paranthracis was subsequently constructed to analyze the evolutionary dynamics involved in ecological adaptation comprehensively. The open pan-genome indicated a flexible gene repertoire with extensive genetic diversity. Significant divergences in the phylogenetic relationships, functional enrichment, and degree of selective pressure between the different components demonstrated different evolutionary dynamics between the core and accessory genomes driven by ecological forces. Purifying selection and gene loss are the main signatures of evolutionary dynamics in B. paranthracis pan-genome. The plasticity of the accessory genome is characterized by horizontal gene transfer (HGT), massive gene losses, and weak purifying or positive selection, which might contribute to niche-specific adaptation. In contrast, although the core genome dominantly undergoes purifying selection, its association with HGT and positively selected mutations indicates its potential role in ecological diversification. Furthermore, host fitness-related dynamics are characterized by the loss of secondary metabolite biosynthesis gene clusters (BGCs) and CAZyme-encoding genes and the acquisition of antimicrobial resistance (AMR) and virulence genes via HGT. This study offers a case study of pan-genome evolution to investigate the ecological adaptations reflected by biogeographical characteristics, thereby advancing the understanding of intraspecific diversity and evolutionary dynamics of foodborne pathogens.

Horizontal Gene Transfer of a key Translation Factor and its Role in Polyproline Proteome Evolution

Prolines cause ribosomes to stall during translation due to their rigid structure. This phenomenon occurs in all domains of life and is exacerbated at polyproline motifs. Such stalling can be eased by the elongation factor P (EF-P) in bacteria. We discovered a potential connection between the loss of ancestral EF-P, the appearance of horizontally transferred EF-P variants, and genomic signs of EF-P dysfunction. Horizontal transfer of the efp gene has occurred several times among bacteria and is associated with the loss of highly conserved polyproline motifs. In this study, we pinpoint cases of horizontal EF-P transfer among a diverse set of bacteria and examine genomic features associated with these events in the phyla Thermotogota and Planctomycetes. In these phyla, horizontal EF-P transfer is also associated with the loss of entire polyproline motif-containing proteins, whose expression is likely dependent on EF-P. In particular, three proteases (Lon, ClpC, and FtsH) and three tRNA synthetases (ValS, IleS1, and IleS2) appear highly sensitive to EF-P transfer. The conserved polyproline motifs within these proteins all reside within close proximity to ATP-binding-regions, some of which are crucial for their function. Our work shows that an ancient EF-P dysfunction has left genomic traces that persist to this day, although it remains unclear whether this dysfunction was strictly due to loss of ancestral EF-P or was related to the appearance of an exogenous variant. The latter possibility would imply that the process of "domesticating" a horizontally transferred efp gene can perturb the overall function of EF-P.

H-NS is a bacterial transposon capture protein

The histone-like nucleoid structuring (H-NS) protein is a DNA binding factor, found in gammaproteobacteria, with functional equivalents in diverse microbes. Universally, such proteins are understood to silence transcription of horizontally acquired genes. Here, we identify transposon capture as a major overlooked function of H-NS. Using genome-scale approaches, we show that H-NS bound regions are transposition "hotspots". Since H-NS often interacts with pathogenicity islands, such targeting creates clinically relevant phenotypic diversity. For example, in Acinetobacter baumannii, we identify altered motility, biofilm formation, and interactions with the human immune system. Transposon capture is mediated by the DNA bridging activity of H-NS and, if absent, more ubiquitous transposition results. Consequently, transcribed and essential genes are disrupted. Hence, H-NS directs transposition to favour evolutionary outcomes useful for the host cell.

Barcoded overexpression screens in gut Bacteroidales identify genes with roles in carbon utilization and stress resistance

A mechanistic understanding of host-microbe interactions in the gut microbiome is hindered by poorly annotated bacterial genomes. While functional genomics can generate large gene-to-phenotype datasets to accelerate functional discovery, their applications to study gut anaerobes have been limited. For instance, most gain-of-function screens of gut-derived genes have been performed in Escherichia coli and assayed in a small number of conditions. To address these challenges, we develop Barcoded Overexpression BActerial shotgun library sequencing (Boba-seq). We demonstrate the power of this approach by assaying genes from diverse gut Bacteroidales overexpressed in Bacteroides thetaiotaomicron. From hundreds of experiments, we identify new functions and phenotypes for 29 genes important for carbohydrate metabolism or tolerance to antibiotics or bile salts. Highlights include the discovery of a D-glucosamine kinase, a raffinose transporter, and several routes that increase tolerance to ceftriaxone and bile salts through lipid biosynthesis. This approach can be readily applied to develop screens in other strains and additional phenotypic assays.

Microbial communities reveal niche partitioning across the slope and bottom zones of the challenger deep

Widespread marine microbiomes exhibit compositional and functional differentiation as a result of adaptation driven by environmental characteristics. We investigated the microbial communities in both seawater and sediments on the slope (7-9 km) and the bottom (9-11 km) of the Challenger Deep of the Mariana Trench to explore community differentiation. Both metagenome-assembled genomes (MAGs) and 16S rRNA amplicon sequence variants (ASVs) showed that the microbial composition in the seawater was similar to that of sediment on the slope, while distinct from that of sediment in the bottom. This scenario suggested a potentially stronger community interaction between seawater and sediment on the slope, which was further confirmed by community assembly and population movement analyses. The metagenomic analysis also indicates a specific stronger potential of nitrate reduction and sulphate assimilation in the bottom seawater, while more versatile nitrogen and sulphur cycling pathways occur on the slope, reflecting functional differentiations among communities in conjunction with environmental features. This work implies that microbial community differentiation occurred in the different hadal niches, and was likely an outcome of microbial adaptation to the extreme hadal trench environment, especially the associated hydrological and geological conditions, which should be considered and measured in situ in future studies.

Measuring the burden of hundreds of BioBricks defines an evolutionary limit on constructability in synthetic biology

Engineered DNA will slow the growth of a host cell if it redirects limiting resources or otherwise interferes with homeostasis. Escape mutants that alleviate this burden can rapidly evolve and take over cell populations, making genetic engineering less reliable and predictable. Synthetic biologists often use genetic parts encoded on plasmids, but their burden is rarely characterized. We measured how 301 BioBrick plasmids affected Escherichia coli growth and found that 59 (19.6%) were burdensome, primarily because they depleted the limited gene expression resources of host cells. Overall, no BioBricks reduced the growth rate of E. coli by >45%, which agreed with a population genetic model that predicts such plasmids should be unclonable. We made this model available online for education ( https://barricklab.org/burden-model ) and added our burden measurements to the iGEM Registry. Our results establish a fundamental limit on what DNA constructs and genetic modifications can be successfully engineered into cells.

Synthetic genomes unveil the effects of synonymous recoding

Engineering the genetic code of an organism provides the basis for (i) making any organism safely resistant to natural viruses and (ii) preventing genetic information flow into and out of genetically modified organisms while (iii) allowing the biosynthesis of genetically encoded unnatural polymers1-4. Achieving these three goals requires the reassignment of multiple of the 64 codons nature uses to encode proteins. However, synonymous codon replacement-recoding-is frequently lethal, and how recoding impacts fitness remains poorly explored. Here, we explore these effects using whole-genome synthesis, multiplexed directed evolution, and genome-transcriptome-translatome-proteome co-profiling on multiple recoded genomes. Using this information, we assemble a synthetic Escherichia coli genome in seven sections using only 57 codons to encode proteins. By discovering the rules responsible for the lethality of synonymous recoding and developing a data-driven multi-omics-based genome construction workflow that troubleshoots synthetic genomes, we overcome the lethal effects of 62,007 synonymous codon swaps and 11,108 additional genomic edits. We show that synonymous recoding induces transcriptional noise including new antisense RNAs, leading to drastic transcriptome and proteome perturbation. As the elimination of select codons from an organism's genetic code results in the widespread appearance of cryptic promoters, we show that synonymous codon choice may naturally evolve to minimize transcriptional noise. Our work provides the first genome-scale description of how synonymous codon changes influence organismal fitness and paves the way for the construction of functional genomes that provide genetic firewalls from natural ecosystems and safely produce biopolymers, drugs, and enzymes with an expanded chemistry.

Timing the evolution of phosphorus-cycling enzymes through geological time using phylogenomics

Phosphorus plays a crucial role in controlling biological productivity, but geological estimates of phosphate concentrations in the Precambrian ocean, during life's origin and early evolution, vary over several orders of magnitude. While reduced phosphorus species may have served as alternative substrates to phosphate, their bioavailability on the early Earth remains unknown. Here, we reconstruct the phylogenomic record of life on Earth and find that phosphate transporting genes (pnas) evolved in the Paleoarchean (ca. 3.6-3.2 Ga) and are consistent with phosphate concentrations above modern levels ( > 3 µM). The first gene optimized for low phosphate levels (pstS; <1 µM) appeared around the same time or in the Mesoarchean depending on the reconstruction method. Most enzymatic pathways for metabolising reduced phosphorus emerged and expanded across the tree of life later. This includes phosphonate-catabolising CP-lyases, phosphite-oxidising pathways and hypophosphite-oxidising pathways. CP-lyases are particularly abundant in dissolved phosphate concentrations below 0.1 µM. Our results thus indicate at least local regions of declining phosphate levels through the Archean, possibly linked to phosphate-scavenging Fe(III), which may have limited productivity. However, reduced phosphorus species did not become widely used until after the Paleoproterozoic Great Oxidation Event (2.3 Ga), possibly linked to expansion of the biosphere at that time.

Comparative genomic analysis and functional investigations for MCs catabolism mechanisms and evolutionary dynamics of MCs-degrading bacteria in ecology

Microcystins (MCs) significantly threaten the ecosystem and public health. Biodegradation has emerged as a promising technology for removing MCs. Many MCs-degrading bacteria have been identified, including an indigenous bacterium Sphingopyxis sp. YF1 that could degrade MC-LR and Adda completely. Herein, we gained insight into the MCs biodegradation mechanisms and evolutionary dynamics of MCs-degrading bacteria, and revealed the toxic risks of the MCs degradation products. The biochemical characteristics and genetic repertoires of strain YF1 were explored. A comparative genomic analysis was performed on strain YF1 and six other MCs-degrading bacteria to investigate their functions. The degradation products were investigated, and the toxicity of the intermediates was analyzed through rigorous theoretical calculation. Strain YF1 might be a novel species that exhibited versatile substrate utilization capabilities. Many common genes and metabolic pathways were identified, shedding light on shared functions and catabolism in the MCs-degrading bacteria. The crucial genes involved in MCs catabolism mechanisms, including mlr and paa gene clusters, were identified successfully. These functional genes might experience horizontal gene transfer events, suggesting the evolutionary dynamics of these MCs-degrading bacteria in ecology. Moreover, the degradation products for MCs and Adda were summarized, and we found most of the intermediates exhibited lower toxicity to different organisms than the parent compound. These findings systematically revealed the MCs catabolism mechanisms and evolutionary dynamics of MCs-degrading bacteria. Consequently, this research contributed to the advancement of green biodegradation technology in aquatic ecology, which might protect human health from MCs.

A global survey of prokaryotic genomes reveals the eco-evolutionary pressures driving horizontal gene transfer

Horizontal gene transfer, the exchange of genetic material through means other than reproduction, is a fundamental force in prokaryotic genome evolution. Genomic persistence of horizontally transferred genes has been shown to be influenced by both ecological and evolutionary factors. However, there is limited availability of ecological information about species other than the habitats from which they were isolated, which has prevented a deeper exploration of ecological contributions to horizontal gene transfer. Here we focus on transfers detected through comparison of individual gene trees to the species tree, assessing the distribution of gene-exchanging prokaryotes across over a million environmental sequencing samples. By analysing detected horizontal gene transfer events, we show distinct functional profiles for recent versus old events. Although most genes transferred are part of the accessory genome, genes transferred earlier in evolution tend to be more ubiquitous within present-day species. We find that co-occurring, interacting and high-abundance species tend to exchange more genes. Finally, we show that host-associated specialist species are most likely to exchange genes with other host-associated specialist species, whereas species found across different habitats have similar gene exchange rates irrespective of their preferred habitat. Our study covers an unprecedented scale of integrated horizontal gene transfer and environmental information, highlighting broad eco-evolutionary trends.

Lethal perturbation of an Escherichia coli regulatory network is triggered by a restriction-modification system's regulator and can be mitigated by excision of the cryptic prophage Rac

Bacterial gene regulatory networks orchestrate responses to environmental challenges. Horizontal gene transfer can bring in genes with regulatory potential, such as new transcription factors (TFs), and this can disrupt existing networks. Serious regulatory perturbations may even result in cell death. Here, we show the impact on Escherichia coli of importing a promiscuous TF that has adventitious transcriptional effects within the cryptic Rac prophage. A cascade of regulatory network perturbations occurred on a global level. The TF, a C regulatory protein, normally controls a Type II restriction-modification system, but in E. coli K-12 interferes with expression of the RacR repressor gene, resulting in de-repression of the normally-silent Rac ydaT gene. YdaT is a prophage-encoded TF with pleiotropic effects on E. coli physiology. In turn, YdaT alters expression of a variety of bacterial regulons normally controlled by the RcsA TF, resulting in deficient lipopolysaccharide biosynthesis and cell division. At the same time, insufficient RacR repressor results in Rac DNA excision, halting Rac gene expression due to loss of the replication-defective Rac prophage. Overall, Rac induction appears to counteract the lethal toxicity of YdaT. We show here that E. coli rewires its regulatory network, so as to minimize the adverse regulatory effects of the imported C TF. This complex set of interactions may reflect the ability of bacteria to protect themselves by having robust mechanisms to maintain their regulatory networks, and/or suggest that regulatory C proteins from mobile operons are under selection to manipulate their host's regulatory networks for their own benefit.

A synthetic biology approach to assemble and reboot clinically relevant Pseudomonas aeruginosa tailed phages

The rise in the frequency of antibiotic resistance has made bacterial infections, specifically Pseudomonas aeruginosa, a cause for greater concern. Phage therapy is a promising solution that uses naturally isolated phages to treat bacterial infections. Ecological limitations, which stipulate a discrete host range and the inevitable evolution of resistance, may be overcome through a better understanding of phage biology and the utilization of engineered phages. In this study, we developed a synthetic biology approach to construct tailed phages that naturally target clinically relevant strains of Pseudomonas aeruginosa. As proof of concept, we successfully cloned and assembled the JG024 and DMS3 phage genomes in yeast using transformation-associated recombination cloning and rebooted these two phage genomes in two different strains of P. aeruginosa. We identified factors that affected phage reboot efficiency like the phage species or the presence of antiviral defense systems in the bacterial strain. We have successfully extended this method to two other phage species and observed that the method enables the reboot of phages that are naturally unable to infect the strain used for reboot. This research represents a critical step toward the construction of clinically relevant, engineered P. aeruginosa phages.IMPORTANCEPseudomonas aeruginosa is a bacterium responsible for severe infections and a common major complication in cystic fibrosis. The use of antibiotics to treat bacterial infections has become increasingly difficult as antibiotic resistance has become more prevalent. Phage therapy is an alternative solution that is already being used in some European countries, but its use is limited by the narrow host range due to the phage receptor specificity, the presence of antiviral defense systems in the bacterial strain, and the possible emergence of phage resistance. In this study, we demonstrate the use of a synthetic biology approach to construct and reboot clinically relevant P. aeruginosa tailed phages. This method enables a significant expansion of possibilities through the construction of engineered phages for therapy applications.

Interaction with a phage gene underlie costs of a β-lactamase

The fitness cost of an antibiotic resistance gene (ARG) can differ across host strains, creating refuges that allow the maintenance of an ARG in the absence of direct selection for its resistance phenotype. Despite the importance of such ARG-host interactions for predicting ARG dynamics, the basis of ARG fitness costs and their variability between hosts are not well understood. We determined the genetic basis of a host-dependent cost of a β-lactamase, blaTEM-116*, that conferred a significant cost in one Escherichia coli strain but was close to neutral in 11 other Escherichia spp. strains. Selection of a blaTEM-116*-encoding plasmid in the strain in which it initially had a high cost resulted in rapid and parallel compensation for that cost through mutations in a P1-like phage gene, relAP1. When the wild-type relAP1 gene was added to a strain in which it was not present and in which blaTEM-116* was neutral, it caused the ARG to become costly. Thus, relAP1 is both necessary and sufficient to explain blaTEM-116* costs in at least some host backgrounds. To our knowledge, these findings represent the first demonstrated case of the cost of an ARG being influenced by a genetic interaction with a phage gene. The interaction between a phage gene and a plasmid-borne ARG highlights the complexity of selective forces determining the maintenance and spread of ARGs and, by extension, encoding phage and plasmids in natural bacterial communities.IMPORTANCEAntibiotic resistance genes (ARGs) play a major role in the increasing problem of antibiotic resistance in clinically relevant bacteria. Selection of these genes occurs in the presence of antibiotics, but their eventual success also depends on the sometimes substantial costs they impose on host bacteria in antibiotic-free environments. We evolved an ARG that confers resistance to penicillin-type antibiotics in one host in which it did confer a cost and in one host in which it did not. We found that costs were rapidly and consistently reduced through parallel genetic changes in a gene encoded by a phage that was infecting the costly host. The unmutated version of this gene was sufficient to cause the ARG to confer a cost in a host in which it was originally neutral, demonstrating an antagonism between the two genetic elements and underlining the range and complexity of pressures determining ARG dynamics in natural populations.

Autotrophic biofilms sustained by deeply sourced groundwater host diverse bacteria implicated in sulfur and hydrogen metabolism

BACKGROUND

Biofilms in sulfide-rich springs present intricate microbial communities that play pivotal roles in biogeochemical cycling. We studied chemoautotrophically based biofilms that host diverse CPR bacteria and grow in sulfide-rich springs to investigate microbial controls on biogeochemical cycling.

RESULTS

Sulfide springs biofilms were investigated using bulk geochemical analysis, genome-resolved metagenomics, and scanning transmission X-ray microscopy (STXM) at room temperature and 87 K. Chemolithotrophic sulfur-oxidizing bacteria, including Thiothrix and Beggiatoa, dominate the biofilms, which also contain CPR Gracilibacteria, Absconditabacteria, Saccharibacteria, Peregrinibacteria, Berkelbacteria, Microgenomates, and Parcubacteria. STXM imaging revealed ultra-small cells near the surfaces of filamentous bacteria that may be CPR bacterial episymbionts. STXM and NEXAFS spectroscopy at carbon K and sulfur L2,3 edges show that filamentous bacteria contain protein-encapsulated spherical elemental sulfur granules, indicating that they are sulfur oxidizers, likely Thiothrix. Berkelbacteria and Moranbacteria in the same biofilm sample are predicted to have a novel electron bifurcating group 3b [NiFe]-hydrogenase, putatively a sulfhydrogenase, potentially linked to sulfur metabolism via redox cofactors. This complex could potentially contribute to symbioses, for example, with sulfur-oxidizing bacteria such as Thiothrix that is based on cryptic sulfur cycling. One Doudnabacteria genome encodes adjacent sulfur dioxygenase and rhodanese genes that may convert thiosulfate to sulfite. We find similar conserved genomic architecture associated with CPR bacteria from other sulfur-rich subsurface ecosystems.

CONCLUSIONS

Our combined metagenomic, geochemical, spectromicroscopic, and structural bioinformatics analyses of biofilms growing in sulfide-rich springs revealed consortia that contain CPR bacteria and sulfur-oxidizing Proteobacteria, including Thiothrix, and bacteria from a new family within Beggiatoales. We infer roles for CPR bacteria in sulfur and hydrogen cycling. Video Abstract.

Invited Review Beyond parasitic convergence: unravelling the evolution of the organellar genomes in holoparasites

BACKGROUND

The molecular evolution of organellar genomes in angiosperms has been studied extensively, with some lineages, such as parasitic ones, displaying unique characteristics. Parasitism has emerged 12 times independently in angiosperm evolution. Holoparasitism is the most severe form of parasitism, and is found in ~10 % of parasitic angiosperms. Although a few holoparasitic species have been examined at the molecular level, most reports involve plastomes instead of mitogenomes. Parasitic plants establish vascular connections with their hosts through haustoria to obtain water and nutrients, which facilitates the exchange of genetic information, making them more susceptible to horizontal gene transfer (HGT). HGT is more prevalent in the mitochondria than in the chloroplast or nuclear compartments.

SCOPE

This review summarizes current knowledge on the plastid and mitochondrial genomes of holoparasitic angiosperms, compares the genomic features across the different lineages, and discusses their convergent evolutionary trajectories and distinctive features. We focused on Balanophoraceae (Santalales), which exhibits extraordinary traits in both their organelles.

CONCLUSIONS

Apart from morphological similarities, plastid genomes of holoparasitic plants also display other convergent features, such as rampant gene loss, biased nucleotide composition and accelerated evolutionary rates. In addition, the plastomes of Balanophoraceae have extremely low GC and gene content, and two unexpected changes in the genetic code. Limited data on the mitochondrial genomes of holoparasitic plants preclude thorough comparisons. Nonetheless, no obvious genomic features distinguish them from the mitochondria of free-living angiosperms, except for a higher incidence of HGT. HGT appears to be predominant in holoparasitic angiosperms with a long-lasting endophytic stage. Among the Balanophoraceae, mitochondrial genomes exhibit disparate evolutionary paths with notable levels of heteroplasmy in Rhopalocnemis and unprecedented levels of HGT in Lophophytum. Despite their differences, these Balanophoraceae share a multichromosomal mitogenome, a feature also found in a few free-living angiosperms.

Four functional profiles for fibre and mucin metabolism in the human gut microbiome

BACKGROUND

With the emergence of metagenomic data, multiple links between the gut microbiome and the host health have been shown. Deciphering these complex interactions require evolved analysis methods focusing on the microbial ecosystem functions. Despite the fact that host or diet-derived fibres are the most abundant nutrients available in the gut, the presence of distinct functional traits regarding fibre and mucin hydrolysis, fermentation and hydrogenotrophic processes has never been investigated.

RESULTS

After manually selecting 91 KEGG orthologies and 33 glycoside hydrolases further aggregated in 101 functional descriptors representative of fibre and mucin degradation pathways in the gut microbiome, we used nonnegative matrix factorization to mine metagenomic datasets. Four distinct metabolic profiles were further identified on a training set of 1153 samples, thoroughly validated on a large database of 2571 unseen samples from 5 external metagenomic cohorts and confirmed with metatranscriptomic data. Profiles 1 and 2 are the main contributors to the fibre-degradation-related metagenome: they present contrasted involvement in fibre degradation and sugar metabolism and are differentially linked to dysbiosis, metabolic disease and inflammation. Profile 1 takes over Profile 2 in healthy samples, and unbalance of these profiles characterize dysbiotic samples. Furthermore, high fibre diet favours a healthy balance between profiles 1 and profile 2. Profile 3 takes over profile 2 during Crohn's disease, inducing functional reorientations towards unusual metabolism such as fucose and H2S degradation or propionate, acetone and butanediol production. Profile 4 gathers under-represented functions, like methanogenesis. Two taxonomic makes up of the profiles were investigated, using either the covariation of 203 prevalent genomes or metagenomic species, both providing consistent results in line with their functional characteristics. This taxonomic characterization showed that profiles 1 and 2 were respectively mainly composed of bacteria from the phyla Bacteroidetes and Firmicutes while profile 3 is representative of Proteobacteria and profile 4 of methanogens.

CONCLUSIONS

Integrating anaerobic microbiology knowledge with statistical learning can narrow down the metagenomic analysis to investigate functional profiles. Applying this approach to fibre degradation in the gut ended with 4 distinct functional profiles that can be easily monitored as markers of diet, dysbiosis, inflammation and disease. Video Abstract.

Metabolically diverse microorganisms mediate methylmercury formation under nitrate-reducing conditions in a dynamic hydroelectric reservoir

Brownlee Reservoir is a mercury (Hg)-impaired hydroelectric reservoir that exhibits dynamic hydrological and geochemical conditions and is located within the Hells Canyon Complex in Idaho, USA. Methylmercury (MeHg) contamination in fish is a concern in the reservoir. While MeHg production has historically been attributed to sulfate-reducing bacteria and methanogenic archaea, microorganisms carrying the hgcA gene are taxonomically and metabolically diverse and the major biogeochemical cycles driving mercury (Hg) methylation are not well understood. In this study, Hg speciation and redox-active compounds were measured throughout Brownlee Reservoir across the stratified period in four consecutive years (2016-2019) to identify the location where and redox conditions under which MeHg is produced. Metagenomic sequencing was performed on a subset of samples to characterize the microbial community with hgcA and identify possible links between biogeochemical cycles and MeHg production. Biogeochemical profiles suggested in situ water column Hg methylation was the major source of MeHg. These profiles, combined with genome-resolved metagenomics focused on hgcA-carrying microbes, indicated that MeHg production occurs in this system under nitrate- or manganese-reducing conditions, which were previously thought to preclude Hg-methylation. Using this multidisciplinary approach, we identified the cascading effects of interannual variability in hydrology on the redox status, microbial metabolic strategies, abundance and metabolic diversity of Hg methylators, and ultimately MeHg concentrations throughout the reservoir. This work expands the known conditions conducive to producing MeHg and suggests that the Hg-methylation mitigation efforts by nitrate or manganese amendment may be unsuccessful in some locations.

Horizontal Gene Transfer and Drug Resistance Involving Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) acquires drug resistance at a rate comparable to that of bacterial pathogens that replicate much faster and have a higher mutation rate. One explanation for this rapid acquisition of drug resistance in Mtb is that drug resistance may evolve in other fast-replicating mycobacteria and then be transferred to Mtb through horizontal gene transfer (HGT). This paper aims to address three questions. First, does HGT occur between Mtb and other mycobacterial species? Second, what genes after HGT tend to survive in the recipient genome? Third, does HGT contribute to antibiotic resistance in Mtb? I present a conceptual framework for detecting HGT and analyze 39 ribosomal protein genes, 23S and 16S ribosomal RNA genes, as well as several genes targeted by antibiotics against Mtb, from 43 genomes representing all major groups within Mycobacterium. I also included mgtC and the insertion sequence IS6110 that were previously reported to be involved in HGT. The insertion sequence IS6110 shows clearly that the Mtb complex participates in HGT. However, the horizontal transferability of genes depends on gene function, as was previously hypothesized. HGT is not observed in functionally important genes such as ribosomal protein genes, rRNA genes, and other genes chosen as drug targets. This pattern can be explained by differential selection against functionally important and unimportant genes after HGT. Functionally unimportant genes such as IS6110 are not strongly selected against, so HGT events involving such genes are visible. For functionally important genes, a horizontally transferred diverged homologue from a different species may not work as well as the native counterpart, so the HGT event involving such genes is strongly selected against and eliminated, rendering them invisible to us. In short, while HGT involving the Mtb complex occurs, antibiotic resistance in the Mtb complex arose from mutations in those drug-targeted genes within the Mtb complex and was not gained through HGT.

Evolutionary barriers to horizontal gene transfer in macrophage-associated Salmonella

Horizontal gene transfer (HGT) is a powerful evolutionary force facilitating bacterial adaptation and emergence of novel phenotypes. Several factors, including environmental ones, are predicted to restrict HGT, but we lack systematic and experimental data supporting these predictions. Here, we address this gap by measuring the relative fitness of 44 genes horizontally transferred from Escherichia coli to Salmonella enterica in infection-relevant environments. We estimated the distribution of fitness effects in each environment and identified that dosage-dependent effects across different environments are a significant barrier to HGT. The majority of genes were found to be deleterious. We also found longer genes had stronger negative fitness consequences than shorter ones, showing that gene length was negatively associated with HGT. Furthermore, fitness effects of transferred genes were found to be environmentally dependent. In summary, a substantial fraction of transferred genes had a significant fitness cost on the recipient, with both gene characteristics and the environment acting as evolutionary barriers to HGT.

The distribution of fitness effects of plasmid pOXA-48 in clinical enterobacteria

Antimicrobial resistance (AMR) in bacteria is a major public health problem. The main route for AMR acquisition in clinically important bacteria is the horizontal transfer of plasmids carrying resistance genes. AMR plasmids allow bacteria to survive antibiotics, but they also entail physiological alterations in the host cell. Multiple studies over the last few years have indicated that these alterations can translate into a fitness cost when antibiotics are absent. However, due to technical limitations, most of these studies are based on analysing new associations between plasmids and bacteria generated in vitro, and we know very little about the effects of plasmids in their native bacterial hosts. In this study, we used a CRISPR-Cas9-tool to selectively cure plasmids from clinical enterobacteria to overcome this limitation. Using this approach, we were able to study the fitness effects of the carbapenem resistance plasmid pOXA-48 in 35 pOXA-48-carrying isolates recovered from hospitalized patients. Our results revealed that pOXA-48 produces variable effects across the collection of wild-type enterobacterial strains naturally carrying the plasmid, ranging from fitness costs to fitness benefits. Importantly, the plasmid was only associated with a significant fitness reduction in four out of 35 clones, and produced no significant changes in fitness in the great majority of isolates. Our results suggest that plasmids produce neutral fitness effects in most native bacterial hosts, helping to explain the great prevalence of plasmids in natural microbial communities.

Evidence of cospeciation between termites and their gut bacteria on a geological time scale

Termites host diverse communities of gut microbes, including many bacterial lineages only found in this habitat. The bacteria endemic to termite guts are transmitted via two routes: a vertical route from parent colonies to daughter colonies and a horizontal route between colonies sometimes belonging to different termite species. The relative importance of both transmission routes in shaping the gut microbiota of termites remains unknown. Using bacterial marker genes derived from the gut metagenomes of 197 termites and one Cryptocercus cockroach, we show that bacteria endemic to termite guts are mostly transferred vertically. We identified 18 lineages of gut bacteria showing cophylogenetic patterns with termites over tens of millions of years. Horizontal transfer rates estimated for 16 bacterial lineages were within the range of those estimated for 15 mitochondrial genes, suggesting that horizontal transfers are uncommon and vertical transfers are the dominant transmission route in these lineages. Some of these associations probably date back more than 150 million years and are an order of magnitude older than the cophylogenetic patterns between mammalian hosts and their gut bacteria. Our results suggest that termites have cospeciated with their gut bacteria since first appearing in the geological record.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Genome-resolved metagenomics reveals depth-related patterns of microbial community structure and functions in a highly stratified, AMD overlaying mine tailings

Acid mine drainage (AMD) is a worldwide environmental problem, yet bioremediation is hampered by a limited knowledge of the reductive microbial processes in the AMD ecosystem. Here, we generate extensive metagenome and geochemical datasets to investigate how microbial populations and metabolic capacities driving major element cycles are structured in a highly stratified, AMD overlaying tailings environment. The results demonstrated an explicit depth-dependent differentiation of microbial community composition and function profiles between the surface and deeper tailings layers, paralleling the dramatic shifts in major physical and geochemical properties. Specifically, key genes involved in sulfur and iron oxidation were significantly enriched in the surface tailings, whereas those associated with reductive nitrogen, sulfur, and iron processes were enriched in the deeper layers. Genome-resolved metagenomics retrieved 406 intermediate or high-quality genomes spanning 26 phyla, including major new groups (e.g., Patescibacteria and DPANN). Metabolic models involving nitrogen, sulfur, iron, and carbon cycles were proposed based on the functional potentials of the abundant microbial genomes, emphasizing syntrophy and the importance of lesser-known taxa in the degradation of complex carbon compounds. These results have implications for in situ AMD bioremediation.

proGenomes3: approaching one million accurately and consistently annotated high-quality prokaryotic genomes

The interpretation of genomic, transcriptomic and other microbial 'omics data is highly dependent on the availability of well-annotated genomes. As the number of publicly available microbial genomes continues to increase exponentially, the need for quality control and consistent annotation is becoming critical. We present proGenomes3, a database of 907 388 high-quality genomes containing 4 billion genes that passed stringent criteria and have been consistently annotated using multiple functional and taxonomic databases including mobile genetic elements and biosynthetic gene clusters. proGenomes3 encompasses 41 171 species-level clusters, defined based on universal single copy marker genes, for which pan-genomes and contextual habitat annotations are provided. The database is available at http://progenomes.embl.de/.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Synthesis and assembly of full-length cyanophage A-4L genome

Artificial cyanophages are considered to be an effective biological method to control harmful cyanobacterial bloom. However, no synthetic cyanophage genome has been constructed and where its obstacles are unclear. Here, we survey a stretch of 16 kb length sequence of cyanophage A-4L that is unclonable in Escherichia coli. We test 12 predicted promoters of cyanophage A-4L which were verified all active in E. coli. Next, we screen for eight ORFs that hindered the assembly of intermediate DNA fragments in E. coli and describe that seven ORFs in the 16 kb sequence could not be separately cloned in E. coli. All of unclonable ORFs in high-copy-number plasmid were successfully cloned using low-copy-number vector, suggesting that these ORFs were copy-number-dependent. We propose a clone strategy abandoned the promotor and the start codon that could be applied for unclonable ORFs. Last, we de novo synthesized and assembled the full-length genome of cyanophage A-4L. This work deepens the understanding of synthetic cyanophages studies.

Comparison of prokaryotes between Mount Everest and the Mariana Trench

BACKGROUND

Mount Everest and the Mariana Trench represent the highest and deepest places on Earth, respectively. They are geographically separated, with distinct extreme environmental parameters that provide unique habitats for prokaryotes. Comparison of prokaryotes between Mount Everest and the Mariana Trench will provide a unique perspective to understanding the composition and distribution of environmental microbiomes on Earth.

RESULTS

Here, we compared prokaryotic communities between Mount Everest and the Mariana Trench based on shotgun metagenomic analysis. Analyzing 25 metagenomes and 1176 metagenome-assembled genomes showed distinct taxonomic compositions between Mount Everest and the Mariana Trench, with little taxa overlap, and significant differences in genome size, GC content, and predicted optimal growth temperature. However, community metabolic capabilities exhibited striking commonality, with > 90% of metabolic modules overlapping among samples of Mount Everest and the Mariana Trench, with the only exception for CO2 fixations (photoautotrophy in Mount Everest but chemoautotrophy in the Mariana Trench). Most metabolic pathways were common but performed by distinct taxa in the two extreme habitats, even including some specialized metabolic pathways, such as the versatile degradation of various refractory organic matters, heavy metal metabolism (e.g., As and Se), stress resistance, and antioxidation. The metabolic commonality indicated the overall consistent roles of prokaryotes in elemental cycling and common adaptation strategies to overcome the distinct stress conditions despite the intuitively huge differences in Mount Everest and the Mariana Trench.

CONCLUSION

Our results, the first comparison between prokaryotes in the highest and the deepest habitats on Earth, may highlight the principles of prokaryotic diversity: although taxa are habitat-specific, primary metabolic functions could be always conserved. Video abstract.

Prokaryotic Pangenomes Act as Evolving Ecosystems

Understanding adaptation to the local environment is a central tenet and a major focus of evolutionary biology. But this is only part of the adaptionist story. In addition to the external environment, one of the main drivers of genome composition is genetic background. In this perspective, I argue that there is a growing body of evidence that intra-genomic selective pressures play a significant part in the composition of prokaryotic genomes and play a significant role in the origin, maintenance and structuring of prokaryotic pangenomes.

Molecular basis for lethal cross-talk between two unrelated bacterial transcription factors - the regulatory protein of a restriction-modification system and the repressor of a defective prophage

Bacterial gene expression depends on the efficient functioning of global transcriptional networks, however their interconnectivity and orchestration rely mainly on the action of individual DNA binding proteins called transcription factors (TFs). TFs interact not only with their specific target sites, but also with secondary (off-target) sites, and vary in their promiscuity. It is not clear yet what mechanisms govern the interactions with secondary sites, and how such rewiring affects the overall regulatory network, but this could clearly constrain horizontal gene transfer. Here, we show the molecular mechanism of one such off-target interaction between two unrelated TFs in Escherichia coli: the C regulatory protein of a Type II restriction-modification system, and the RacR repressor of a defective prophage. We reveal that the C protein interferes with RacR repressor expression, resulting in derepression of the toxic YdaT protein. These results also provide novel insights into regulation of the racR-ydaST operon. We mapped the C regulator interaction to a specific off-target site, and also visualized C protein dynamics, revealing intriguing differences in single molecule dynamics in different genetic contexts. Our results demonstrate an apparent example of horizontal gene transfer leading to adventitious TF cross-talk with negative effects on the recipient's viability. More broadly, this study represents an experimentally-accessible model of a regulatory constraint on horizontal gene transfer.

Phage-Encoded Sigma Factors Alter Bacterial Dormancy

By entering a reversible state of reduced metabolic activity, dormant microorganisms are able to tolerate suboptimal conditions that would otherwise reduce their fitness. Dormancy may also benefit bacteria by serving as a refuge from parasitic infections. Here, we focus on dormancy in the Bacillota, where endospore development is transcriptionally regulated by the expression of sigma factors. A disruption of this process could influence the survivorship or reproduction of phages that infect spore-forming hosts with implications for coevolutionary dynamics. We characterized the distribution of sigma factors in over 4,000 genomes of diverse phages capable of infecting hosts that span the bacterial domain. From this, we identified homologs of sporulation-specific sigma factors in phages that infect spore-forming hosts. Unlike sigma factors required for phage reproduction, we provide evidence that sporulation-like sigma factors are nonessential for lytic infection. However, when expressed in the spore-forming Bacillus subtilis, some of these phage-derived sigma factors can activate the bacterial sporulation gene network and lead to a reduction in spore yield. Our findings suggest that the acquisition of host-like transcriptional regulators may allow phages to manipulate a complex and ancient trait in one of the most abundant cell types on Earth. IMPORTANCE As obligate parasites, phages exert strong top-down pressure on host populations with eco-evolutionary implications for community dynamics and ecosystem functioning. The process of phage infection, however, is constrained by bottom-up processes that influence the energetic and nutritional status of susceptible hosts. Many phages have acquired auxiliary genes from bacteria, which can be used to exploit host metabolism with consequences for phage fitness. In this study, we demonstrate that phages infecting spore-forming bacteria carry homologs of sigma factors, which their hosts use to orchestrate gene expression during spore development. By tapping into regulatory gene networks, phages may manipulate the physiology and survival strategies of nongrowing bacteria in ways that influence host-parasite coevolution.

A Siderophore-Encoding Plasmid Encodes High-Level Virulence in Escherichia coli

A plasmid that harbored the virulence factors highly like those of the virulence plasmid commonly found in clinical hypervirulent Klebsiella pneumoniae strains was detected in a foodborne Escherichia coli strain EC1108 and designated p1108-IncFIB. This virulent-like plasmid was found to be common in E. coli from various sources. To understand the contribution of this plasmid to the virulence of E. coli, plasmid p1108-IncFIB in strain EC1108 was first cured to generate strain EC1108-PC. The virulence plasmid p15WZ-82_Vir in Klebsiella pneumoniae strain 15WZ-82 was then transmitted to EC1108-PC to produce the transconjugant, EC1108-PC-TC to assess the contribution of this virulence plasmid to the virulence level of E. coli. During the process of conjugation, p15WZ-82_Vir was found to be evolved into p15WZ-82_int, which underwent homologous recombination with a plasmid encoding a carbapenemase gene, bla_NDM-1, p1108-NDM, in EC1108-PC. Comparison between the level of virulence in the EC1108, EC1108-PC-TC, and EC1108-PC through serum and macrophage resistance assay, as well as animal experiments, confirmed that plasmid p1108-IncFIB encoded a high level of virulence in E. coli, yet the fusion plasmid derived from p15WZ-82_Vir did not encode virulence but instead imposed a high fitness cost in the E. coli strain EC1108-PC-TC. These findings indicate that E. coli strains carrying the virulence plasmid p1108-IncFIB in multidrug-resistant (MDR) strains may also impose serious public health threats like that of hypervirulent Klebsiella pneumoniae strains harboring the p15WZ-82_Vir plasmid.

IMPORTANCE Acquisition of pLVPK-like virulence plasmid by Klebsiella pneumoniae converts it to hypervirulent K. pneumoniae (HvKP), which has become one of the most important clinical bacterial pathogens. The potential of transmission of this virulence plasmid and its contribution to the virulence of other Enterobacteriaceae, such as E. coli, are not clear yet. In this study, we showed that pLVPK-like virulence plasmid exhibited fitness costs and did not contribute to the virulence in E. coli. However, we identified a novel virulence plasmid, p1108-IncFIB, that encodes similar siderophore genes as those of pLVPK from a foodborne E. coli strain and showed that p1108-IncFIB encoded a high level of virulence in E. coli. BLAST of E. coli genomes from GenBank showed that these siderophore genes were widespread in clinical E. coli strains. Further studies are warranted to understand the impact of this plasmid in the control of clinical infections caused by E. coli.

Mutational Pattern Induced by 5-Fluorouracil and Oxaliplatin in the Gut Microbiome

Chemotherapeutic agents, such as 5-fluorouracil (5-FU) and oxaliplatin (Oxi), can not only kill the cancer cell but also influence the proliferation of gut microbiota; however, the interaction between these drugs and gut microbiota remains poorly understood. In this study, we developed a powerful framework for taxonomy composition and genomic variation analysis to investigate the mutagenesis effect and proliferation influence of chemotherapeutic agents, such as 5-FU and Oxi, on gut microbiota and the interaction between these drugs and gut microbiota during chemotherapy. Using the gut microbiome data, we detected 1.45 million variations among the chemotherapy groups and found the drugs significantly affected mutation signatures of gut microbiota. Oxi notably increased transversion rate, whereas 5-FU reduced the rate. Traits related to cell division and nutrient mobilization showed evidence of strong selection pressure from chemotherapeutic agents. In addition, drug-associated bacteriome shift patterns and functional alterations were found: the metabolism changes in the 5-FU group implied that gut microbiota could provide additional nicotinamide adenine dinucleotide (NAD⁺) to inhibit cancer cell autophagy; in the Oxi group, the ribosome and lysine biosynthesis genes were obviously enriched. Our study provides a blueprint for characterizing the role of microbes and drug–microbe interaction in the gut microbiota response to chemotherapy.

Synthetic Genomics From a Yeast Perspective

Synthetic Genomics focuses on the construction of rationally designed chromosomes and genomes and offers novel approaches to study biology and to construct synthetic cell factories. Currently, progress in Synthetic Genomics is hindered by the inability to synthesize DNA molecules longer than a few hundred base pairs, while the size of the smallest genome of a self-replicating cell is several hundred thousand base pairs. Methods to assemble small fragments of DNA into large molecules are therefore required. Remarkably powerful at assembling DNA molecules, the unicellular eukaryote Saccharomyces cerevisiae has been pivotal in the establishment of Synthetic Genomics. Instrumental in the assembly of entire genomes of various organisms in the past decade, the S. cerevisiae genome foundry has a key role to play in future Synthetic Genomics developments.

Reassembling a cannon in the DNA defense arsenal: Genetics of StySA, a BREX phage exclusion system in Salmonella lab strains

Understanding mechanisms that shape horizontal exchange in prokaryotes is a key problem in biology. A major limit on DNA entry is imposed by restriction-modification (RM) processes that depend on the pattern of DNA modification at host-specified sites. In classical RM, endonucleolytic DNA cleavage follows detection of unprotected sites on entering DNA. Recent investigation has uncovered BREX (BacteRiophage EXclusion) systems. These RM-like activities employ host protection by DNA modification, but immediate replication arrest occurs without evident of nuclease action on unmodified phage DNA. Here we show that the historical stySA RM locus of Salmonella enterica sv Typhimurium is a variant BREX system. A laboratory strain disabled for both the restriction and methylation activity of StySA nevertheless has wild type sequence in pglX, the modification gene homolog. Instead, flanking genes pglZ and brxC each carry multiple mutations (μ) in their C-terminal domains. We further investigate this system in situ, replacing the mutated pglZμ and brxCμ genes with the WT counterpart. PglZ-WT supports methylation in the presence of either BrxCμ or BrxC-WT but not in the presence of a deletion/insertion allele, ΔbrxC::cat. Restriction requires both BrxC-WT and PglZ-WT, implicating the BrxC C-terminus specifically in restriction activity. These results suggests that while BrxC, PglZ and PglX are principal components of the BREX modification activity, BrxL is required for restriction only. Furthermore, we show that a partial disruption of brxL disrupts transcription globally.

Horizontal gene transfer is a major driver of evolution and adaptation in bacteria. Genes from outside may be beneficial or dangerous to the receiving cell. Benefits include new food sources such as sugars, or new homes by adhesion, or new resistances, as to antibiotics. Dangers are posed by bacteriophages--viruses that take over the cell machinery, multiply, and release progeny to kill sister cells. Host-dependent restriction-modification systems enable defense that distinguishes relatives from strangers: using a modification pattern (M) carried by DNA bases added by the host cell to prevent restriction (R). Sisters and cousin cells will have the same protective pattern on DNA, while DNA of foreign origin will have the wrong M pattern and be restricted (R, rejected). Typically, restriction involves nuclease digestion. Here we address the enigmatic StySA RM system, one of the earliest to be genetically characterized. It is a variant of the newly recognized defense mechanism, BREX. BREX systems also track DNA history via modification pattern, but restrict by a novel, uncharacterized mechanism. Like other BREX family systems, StySA-BREX modification requires multiple components. When StySA-BREX transcription is unbalanced, we find global disruption of gene transcription. The disruption pattern does not suggest SOS-inducing damage to DNA.

Recombination resolves the cost of horizontal gene transfer in experimental populations of Helicobacter pylori

Horizontal gene transfer (HGT)—the transfer of DNA between lineages—is responsible for a large proportion of the genetic variation that contributes to evolution in microbial populations. While HGT can bring beneficial genetic innovation, the transfer of DNA from other species or strains can also have deleterious effects. In this study, we evolve populations of the bacteria Helicobacter pylori and use DNA sequencing to identify over 40,000 genetic variants transferred by HGT. We measure the cost of many of these and find that both strongly beneficial mutations and deleterious mutations are genetic variants transferred by natural transformation. Importantly, we also show how recombination that separates linked beneficial and deleterious mutations resolves the cost of HGT.

Horizontal gene transfer (HGT) is important for microbial evolution, yet we know little about the fitness effects and dynamics of horizontally transferred genetic variants. In this study, we evolve laboratory populations of Helicobacter pylori, which take up DNA from their environment by natural transformation, and measure the fitness effects of thousands of transferred genetic variants. We find that natural transformation increases the rate of adaptation but comes at the cost of significant genetic load. We show that this cost is circumvented by recombination, which increases the efficiency of selection by decoupling deleterious and beneficial genetic variants. Our results show that adaptation with HGT, pervasive in natural microbial populations, is shaped by a combination of selection, recombination, and genetic drift not accounted for in existing models of evolution.

Strain-level characterization of broad host range mobile genetic elements transferring antibiotic resistance from the human microbiome

Mobile genetic elements (MGEs) carrying antibiotic resistance genes (ARGs) disseminate ARGs when they mobilise into new bacterial hosts. The nature of such horizontal gene transfer (HGT) events between human gut commensals and pathogens remain poorly characterised. Here, we compare 1354 cultured commensal strains (540 species) to 45,403 pathogen strains (12 species) and find 64,188 MGE-mediated ARG transfer events between the two groups using established methods. Among the 5931 MGEs, we find 15 broad host range elements predicted to have crossed different bacterial phyla while also occurring in animal and environmental microbiomes. We experimentally demonstrate that predicted broad host range MGEs can mobilise from commensals Dorea longicatena and Hungatella hathewayi to pathogen Klebsiella oxytoca, crossing phyla simultaneously. Our work establishes the MGE-mediated ARG dissemination network between human gut commensals and pathogens and highlights broad host range MGEs as targets for future ARG dissemination management.

Gene sharing among plasmids and chromosomes reveals barriers for antibiotic resistance gene transfer

The emergence of antibiotic resistant bacteria is a major threat to modern medicine. Rapid adaptation to antibiotics is often mediated by the acquisition of plasmids carrying antibiotic resistance (ABR) genes. Nonetheless, the determinants of plasmid-mediated ABR gene transfer remain debated. Here, we show that the propensity of ABR gene transfer via plasmids is higher for accessory chromosomal ABR genes in comparison with core chromosomal ABR genes, regardless of the resistance mechanism. Analysing the pattern of ABR gene occurrence in the genomes of 2635 Enterobacteriaceae isolates, we find that 33% of the 416 ABR genes are shared between chromosomes and plasmids. Phylogenetic reconstruction of ABR genes occurring on both plasmids and chromosomes supports their evolution by lateral gene transfer. Furthermore, accessory ABR genes (encoded in less than 10% of the chromosomes) occur more abundantly in plasmids in comparison with core ABR genes (encoded in greater than or equal to 90% of the chromosomes). The pattern of ABR gene occurrence in plasmids and chromosomes is similar to that in the total Escherichia genome. Our results thus indicate that the previously recognized barriers for gene acquisition by lateral gene transfer apply also to ABR genes. We propose that the functional complexity of the underlying ABR mechanism is an important determinant of ABR gene transferability. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Cargo Genes of Tn7-Like Transposons Comprise an Enormous Diversity of Defense Systems, Mobile Genetic Elements, and Antibiotic Resistance Genes

Transposition is a major mechanism of horizontal gene mobility in prokaryotes. However, exploration of the genes mobilized by transposons (cargo) is hampered by the difficulty in delineating integrated transposons from their surrounding genetic context. Here, we present a computational approach that allowed us to identify the boundaries of 6,549 Tn7-like transposons. We found that 96% of these transposons carry at least one cargo gene. Delineation of distinct communities in a gene-sharing network demonstrates how transposons function as a conduit of genes between phylogenetically distant hosts. Comparative analysis of the cargo genes reveals significant enrichment of mobile genetic elements (MGEs) nested within Tn7-like transposons, such as insertion sequences and toxin-antitoxin modules, and of genes involved in recombination, anti-MGE defense, and antibiotic resistance. More unexpectedly, cargo also includes genes encoding central carbon metabolism enzymes. Twenty-two Tn7-like transposons carry both an anti-MGE defense system and antibiotic resistance genes, illustrating how bacteria can overcome these combined pressures upon acquisition of a single transposon. This work substantially expands the distribution of Tn7-like transposons, defines their evolutionary relationships, and provides a large-scale functional classification of prokaryotic genes mobilized by transposition.

Deeplasmid: deep learning accurately separates plasmids from bacterial chromosomes

Plasmids are mobile genetic elements that play a key role in microbial ecology and evolution by mediating horizontal transfer of important genes, such as antimicrobial resistance genes. Many microbial genomes have been sequenced by short read sequencers and have resulted in a mix of contigs that derive from plasmids or chromosomes. New tools that accurately identify plasmids are needed to elucidate new plasmid-borne genes of high biological importance. We have developed Deeplasmid, a deep learning tool for distinguishing plasmids from bacterial chromosomes based on the DNA sequence and its encoded biological data. It requires as input only assembled sequences generated by any sequencing platform and assembly algorithm and its runtime scales linearly with the number of assembled sequences. Deeplasmid achieves an AUC–ROC of over 89%, and it was more accurate than five other plasmid classification methods. Finally, as a proof of concept, we used Deeplasmid to predict new plasmids in the fish pathogen Yersinia ruckeri ATCC 29473 that has no annotated plasmids. Deeplasmid predicted with high reliability that a long assembled contig is part of a plasmid. Using long read sequencing we indeed validated the existence of a 102 kb long plasmid, demonstrating Deeplasmid's ability to detect novel plasmids.

A CRISPR/dCas9-assisted system to clone toxic genes in Escherichia coli

BACKGROUND

The cloning of toxic genes in E. coli requires strict regulation of the target genes' leaky expression. Many methods facilitating successful gene cloning of toxic genes are commonly exploited, but the applicability is severely limited.

METHODS

A CRISPR/dCas9-assisted system was used to clone toxic genes in E. coli. The plasmid-based and genome-integrated systems were designed in this study. And the green fluorescent protein characterization system was used to test the repression efficiency of the two systems.

RESULTS

We optimized the plasmid-based CRISPR/dCas9-assisted repression system via testing different sgRNAs targeting the Ptrc promoter and achieved inhibition efficiency up to 64.8%. The genome-integrated system represented 35.9% decreased GFP expression and was successfully employed to cloned four toxic genes from Corynebacterium glutamicum in E. coli.

CONCLUSIONS

Using this method, we successfully cloned four C. glutamicum-derived toxic genes that had been failed to clone in conventional ways. The CRISPR/dCas9-assisted gene cloning method was a promising tool to facilitate precise gene cloning of different origins in E. coli.

GENERAL SIGNIFICANCE

This system will be useful for cloning toxic genes from different origins in E. coli, and can accelerate the related research of gene characterization and heterologous expression in the metagenomic era.

mOTUs: Profiling Taxonomic Composition, Transcriptional Activity and Strain Populations of Microbial Communities

The mOTU profiler, or mOTUs for short, is a software tool that enables the profiling of microbial communities in terms of their taxonomic composition, relative abundance of metabolically active members, and diversity of strain populations. To this end, it maintains a database of single-copy phylogenetic marker gene sequences, which are used as a reference to which short read metagenomic and metatranscriptomic reads are mapped for the identification and quantification of microbial taxa. Here, we describe the most common use cases of the mOTU profiler in two basic protocols. Additional supporting protocols provide information on its installation and in-depth guidance on adjusting its settings for increasing or decreasing the stringency with which taxa are detected and quantified, as well as for customizing the output file format. Guidelines for understanding the profiling results are provided, along with additional information on unique features, methodological details, and the development history of the tool. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Metagenomic and metatranscriptomic mOTU profiling Basic Protocol 2: Metagenomic SNV profiling Support Protocol 1: Installing mOTUs Support Protocol 2: Profiling pipeline-step by step Support Protocol 3: The mOTUs profiling routine using advanced parameters Support Protocol 4: Metagenomic SNV calling: advanced parameters.

An efficient, non-phylogenetic method for detecting genes sharing evolutionary signals in phylogenomic datasets

Assessing the compatibility between gene family phylogenies is a crucial and often computationally demanding step in many phylogenomic analyses. Here, we describe the Evolutionary Similarity Index (IES), a means to assess shared evolution between gene families using a weighted orthogonal distance regression model applied to sequence distances. The utilization of pairwise distance matrices circumvents comparisons between gene tree topologies, which are inherently uncertain and sensitive to evolutionary model choice, phylogenetic reconstruction artifacts, and other sources of error. Furthermore, IES enables the many-to-many pairing of multiple copies between similarly evolving gene families. This is done by selecting non-overlapping pairs of copies, one from each assessed family, and yielding the least sum of squared residuals. Analyses of simulated gene family data sets show that IES’s accuracy is on par with popular tree-based methods while also less susceptible to noise introduced by sequence alignment and evolutionary model fitting. Applying IES to an empirical data set of 1,322 genes from 42 archaeal genomes identified eight major clusters of gene families with compatible evolutionary trends. The most cohesive cluster consisted of 62 genes with compatible evolutionary signal, which occur as both single-copy and multiple homologs per genome; phylogenetic analysis of concatenated alignments from this cluster produced a tree closely matching previously published species trees for Archaea. Four other clusters are mainly composed of accessory genes with limited distribution among Archaea and enriched toward specific metabolic functions. Pairwise evolutionary distances obtained from these accessory gene clusters suggest patterns of interphyla horizontal gene transfer. An IES implementation is available at https://github.com/lthiberiol/evolSimIndex.

Essential gene acquisition destabilizes plasmid inheritance

Extra-chromosomal genetic elements are important drivers of evolutionary transformations and ecological adaptations in prokaryotes with their evolutionary success often depending on their ‘utility’ to the host. Examples are plasmids encoding antibiotic resistance genes, which are known to proliferate in the presence of antibiotics. Plasmids carrying an essential host function are recognized as permanent residents in their host. Essential plasmids have been reported in several taxa where they often encode essential metabolic functions; nonetheless, their evolution remains poorly understood. Here we show that essential genes are rarely encoded on plasmids; evolving essential plasmids in Escherichia coli we further find that acquisition of an essential chromosomal gene by a plasmid can lead to plasmid extinction. A comparative genomics analysis of Escherichia isolates reveals few plasmid-encoded essential genes, yet these are often integrated into plasmid-related functions; an example is the GroEL/GroES chaperonin. Experimental evolution of a chaperonin-encoding plasmid shows that the acquisition of an essential gene reduces plasmid fitness regardless of the stability of plasmid inheritance. Our results suggest that essential plasmid emergence leads to a dose effect caused by gene redundancy. The detrimental effect of essential gene acquisition on plasmid inheritance constitutes a barrier for plasmid-mediated lateral gene transfer and supplies a mechanistic understanding for the rarity of essential genes in extra-chromosomal genetic elements.

Mobile genetic elements have been extensively studied due to their role as agents of genetic innovation and rapid adaptation in prokaryotes. Specifically, prokaryotic plasmids have been the focus of investigation in the context of bacterial survival under growth limiting conditions with the prime example of resistance to antibiotics and heavy metals. In contrast, plasmids that encode for functions that are essential to their host viability are rarely described. We investigate the evolution of plasmids that encode for genes previously identified as essential for bacterial life. Our analysis of Escherichia isolates reveals only few plasmid-encoded essential genes, which likely function in the plasmid rather than the host life cycle. Following the evolution of plasmids encoding an essential gene in Escherichia coli in real time, we further find that the acquisition of a chromosomal essential gene may lead to plasmid loss. Our study supplies data and a mechanistic understanding on the rarity of essential genes in mobile genetic elements. We conclude that prokaryotic plasmids are rarely essential for their bacterial host.

Nucleoid-associated Rok differentially affects chromosomal transformation on Bacillus subtilis recombination-deficient cells

Rok, a Bacillus subtilis nucleoid-associated protein (NAP), negatively regulates competence development and silences xenogeneic genes. We show that rok inactivation increases rpoB482 natural intraspecies chromosomal transformation (CT) and plasmid transformation to a different extent. In ΔaddAB, ΔrecO, recF15, ΔrecU, ΔruvAB or rec⁺ cells intraspecies CT significantly increases, but the ΔrecD2 mutation reduces, and the ΔrecX, ΔradA or ΔdprA mutation further decreases CT in the Δrok context when compared to rok⁺ cells. These observations support the idea that rok inactivation, by altering the topology of the recipient DNA, differentially affects the integration of homologous DNA in rec-deficient strains, and in minor extent the competent subpopulation size. The impairment of other NAP (Hbsu or LrpC) also increased intra- and interspecies CT (nonself-DNA, ~8% nucleotide sequence divergence) in rec⁺ cells, but differentially reduced both types of CTs in certain rec-deficient strains. We describe that rok inactivation significantly stimulates intra and interspecies CT but differentially reduces them in transformation-deficient cells, perhaps by altering the nucleoid architecture. We extend the observation to other NAPs (Hbsu, LrpC).

Helarchaeota and co-occurring sulfate-reducing bacteria in subseafloor sediments from the Costa Rica Margin

Deep sediments host many archaeal lineages, including the Asgard superphylum which contains lineages predicted to require syntrophic partnerships. Our knowledge about sedimentary archaeal diversity and their metabolic pathways and syntrophic partners is still very limited. We present here new genomes of Helarchaeota and the co-occurring sulfate-reducing bacteria (SRB) recovered from organic-rich sediments off Costa Rica Margin. Phylogenetic analyses revealed three new metagenome-assembled genomes (MAGs) affiliating with Helarchaeota, each of which has three variants of the methyl-CoM reductase-like (MCR-like) complex that may enable them to oxidize short-chain alkanes anaerobically. These Helarchaeota have no multi-heme cytochromes but have Group 3b and Group 3c [NiFe] hydrogenases, and formate dehydrogenase, and therefore have the capacity to transfer the reducing equivalents (in the forms of hydrogen and formate) generated from alkane oxidation to external partners. We also recovered five MAGs of SRB affiliated with the class of Desulfobacteria, two of which showed relative abundances (represented by genome coverages) positively correlated with those of the three Helarchaeota. Genome analysis suggested that these SRB bacteria have the capacity of H2 and formate utilization and could facilitate electron transfers from other organisms by means of these reduced substances. Their co-occurrence and metabolic features suggest that Helarchaeota may metabolize synergistically with some SRB, and together exert an important influence on the carbon cycle by mitigating the hydrocarbon emission from sediments to the overlying ocean.

Genetic Context Significantly Influences the Maintenance and Evolution of Degenerate Pathways

Understanding the evolution of novel physiological traits is highly relevant for expanding the characterization and manipulation of biological systems. Acquisition of new traits can be achieved through horizontal gene transfer (HGT). Here, we investigate drivers that promote or deter the maintenance of HGT-driven degeneracy, occurring when processes accomplish identical functions through nonidentical components. Subsequent evolution can optimize newly acquired functions; for example, beneficial alleles identified in an engineered Methylorubrum extorquens strain allowed it to utilize a “Foreign” formaldehyde oxidation pathway substituted for its Native pathway for methylotrophic growth. We examined the fitness consequences of interactions between these alleles when they were combined with the Native pathway or both (Dual) pathways. Unlike the Foreign pathway context where they evolved, these alleles were often neutral or deleterious when moved into these alternative genetic backgrounds. However, there were instances where combinations of multiple alleles resulted in higher fitness outcomes than individual allelic substitutions could provide. Importantly, the genetic context accompanying these allelic substitutions significantly altered the fitness landscape, shifting local fitness peaks and restricting the set of accessible evolutionary trajectories. These findings highlight how genetic context can negatively impact the probability of maintaining native and HGT-introduced functions together, making it difficult for degeneracy to evolve. However, in cases where the cost of maintaining degeneracy was mitigated by adding evolved alleles impacting the function of these pathways, we observed rare opportunities for pathway coevolution to occur. Together, our results highlight the importance of genetic context and resulting epistasis in retaining or losing HGT-acquired degenerate functions.

Fasting alters the gut microbiome reducing blood pressure and body weight in metabolic syndrome patients

Periods of fasting and refeeding may reduce cardiometabolic risk elevated by Western diet. Here we show in the substudy of NCT02099968, investigating the clinical parameters, the immunome and gut microbiome exploratory endpoints, that in hypertensive metabolic syndrome patients, a 5-day fast followed by a modified Dietary Approach to Stop Hypertension diet reduces systolic blood pressure, need for antihypertensive medications, body-mass index at three months post intervention compared to a modified Dietary Approach to Stop Hypertension diet alone. Fasting alters the gut microbiome, impacting bacterial taxa and gene modules associated with short-chain fatty acid production. Cross-system analyses reveal a positive correlation of circulating mucosa-associated invariant T cells, non-classical monocytes and CD4⁺ effector T cells with systolic blood pressure. Furthermore, regulatory T cells positively correlate with body-mass index and weight. Machine learning analysis of baseline immunome or microbiome data predicts sustained systolic blood pressure response within the fasting group, identifying CD8⁺ effector T cells, Th17 cells and regulatory T cells or Desulfovibrionaceae, Hydrogenoanaerobacterium, Akkermansia, and Ruminococcaceae as important contributors to the model. Here we report that the high-resolution multi-omics data highlight fasting as a promising non-pharmacological intervention for the treatment of high blood pressure in metabolic syndrome patients.

Nutritional modification including fasting has been shown to reduce cardiometabolic risk linked to western diet. Here the authors show implementation of fasting resulted in alterations to the intestinal microbiota, and circulating immune cells, improving blood pressure and body weight in patients with metabolic syndrome.

Evidence for Horizontal and Vertical Transmission of Mtr-Mediated Extracellular Electron Transfer among the Bacteria

Some bacteria and archaea have evolved the means to use extracellular electron donors and acceptors for energy metabolism, a phenomenon broadly known as extracellular electron transfer (EET). One such EET mechanism is the transmembrane electron conduit MtrCAB, which has been shown to transfer electrons derived from metabolic substrates to electron acceptors, like Fe(III) and Mn(IV) oxides, outside the cell. Although most studies of MtrCAB-mediated EET have been conducted in Shewanella oneidensis MR-1, recent investigations in Vibrio and Aeromonas species have revealed that the electron-donating proteins that support MtrCAB in Shewanella are not as representative as previously thought. This begs the question of how widespread the capacity for MtrCAB-mediated EET is, the changes it has accrued in different lineages, and where these lineages persist today. Here, we employed a phylogenetic and comparative genomics approach to identify the MtrCAB system across all domains of life. We found mtrCAB in the genomes of numerous diverse Bacteria from a wide range of environments, and the patterns therein strongly suggest that mtrCAB was distributed through both horizontal and subsequent vertical transmission, and with some cases indicating downstream modular diversification of both its core and accessory components. Our data point to an emerging evolutionary story about metal-oxidizing and -reducing metabolism, demonstrates that this capacity for EET has broad relevance to a diversity of taxa and the biogeochemical cycles they drive, and lays the foundation for further studies to shed light on how this mechanism may have coevolved with Earth's redox landscape. IMPORTANCE While many metabolisms make use of soluble, cell-permeable substrates like oxygen or hydrogen, there are other substrates, like iron or manganese, that cannot be brought into the cell. Some bacteria and archaea have evolved the means to directly "plug in" to such environmental electron reservoirs in a process known as extracellular electron transfer (EET), making them powerful agents of biogeochemical change and promising vehicles for bioremediation and alternative energy. Yet the diversity, distribution, and evolution of EET mechanisms are poorly constrained. Here, we present findings showing that the genes encoding one such EET system (mtrCAB) are present in a broad diversity of bacteria found in a wide range of environments, emphasizing the ubiquity and potential impact of EET in our biosphere. Our results suggest that these genes have been disseminated largely through horizontal transfer, and the changes they have accrued in these lineages potentially reflect adaptations to changing environments.