Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

Four novel Acinetobacter lwoffii strains isolated from the milk of cows in China with subclinical mastitis

BACKGROUND

Acinetobacter lwoffii (A. lwoffii) is a Gram-negative bacteria common in the environment, and it is the normal flora in human respiratory and digestive tracts. The bacteria is a zoonotic and opportunistic pathogen that causes various infections, including nosocomial infections. The aim of this study was to identify A. lwoffii strains isolated from bovine milk with subclinical mastitis in China and get a better understanding of its antimicrobial susceptibility and resistance profile. This is the first study to analyze the drug resistance spectrum and corresponding mechanisms of A. lwoffii isolated in raw milk.

RESULTS

Four A. lwoffii strains were isolated by PCR method. Genetic evolution analysis using the neighbor-joining method showed that the four strains had a high homology with Acinetobacter lwoffii. The strains were resistant to several antibiotics and carried 17 drug-resistance genes across them. Specifically, among 23 antibiotics, the strains were completely susceptible to 6 antibiotics, including doxycycline, erythromycin, polymyxin, clindamycin, imipenem, and meropenem. In addition, the strains showed variable resistance patterns. A total of 17 resistance genes, including plasmid-mediated resistance genes, were detected across the four strains. These genes mediated resistance to 5 classes of antimicrobials, including beta-lactam, aminoglycosides, fluoroquinolones, tetracycline, sulfonamides, and chloramphenicol.

CONCLUSION

These findings indicated that multi-drug resistant Acinetobacter lwoffii strains exist in raw milk of bovine with subclinical mastitis. Acinetobacter lwoffii are widespread in natural environmental samples, including water, soil, bathtub, soap box, skin, pharynx, conjunctiva, saliva, gastrointestinal tract, and vaginal secretions. The strains carry resistance genes in mobile genetic elements to enhance the spread of these genes. Therefore, more attention should be paid to epidemiological surveillance and drug resistant A. lwoffii.

The Net of Life, a short story: Intricate patterns of gene flows across hundreds of extant genomes, all the way to LUCA

Over the past quarter-century, the field of evolutionary biology has been transformed by the emergence of complete genome sequences and the conceptual framework known as the 'Net of Life.' This paradigm shift challenges traditional notions of evolution as a tree-like process, emphasizing the complex, interconnected network of gene flow that may blur the boundaries between distinct lineages. In this context, gene loss, rather than horizontal gene transfer, is the primary driver of gene content, with vertical inheritance playing a principal role. The 'Net of Life' not only impacts our understanding of genome evolution but also has profound implications for classification systems, the rapid appearance of new traits, and the spread of diseases. Here, we explore the core tenets of the 'Net of Life' and its implications for genome-scale phylogenetic divergence, providing a comprehensive framework for further investigations in evolutionary biology.

Unlocking the genetic potential of Lacticaseibacillus rhamnosus strains: Medical applications of a promising probiotic for human and animal health

Lacticaseibacillus rhamnosus is a group of probiotic strains that have gained popularity for their potential health benefits such as promoting digestive health, boosting the immune system, improving lactose digestion, preventing and treating antibiotic-associated diarrhea, reducing the severity and duration of certain infections, and preventing the formation of dental plaque. In particular, L. rhamnosus strains SD4 and SD11 have promising human and animal health applications due to their ability to inhibit the growth of harmful pathogens. This study presents an in silico genomic analysis of L. rhamnosus strains SD4 and SD11. We analyzed draft genomes and conducted comparative genome analyses against several other probiotic strains, aiming to gain insights into the genomes of the two strains and to compare them to related strains isolated from other sources. We also aimed to clarify the functional mechanisms and adaptation of these strains to specific environments. Comprehensive insights into the genomes of L. rhamnosus SD4 and SD11 could enhance our understanding of their capacity to colonize, adapt, and exhibit probiotic properties after administration. This study holds significance in advancing our understanding of the potential health benefits associated with these strains and in elucidating the underlying mechanisms responsible for their effectiveness in humans and animals.

Long range segmentation of prokaryotic genomes by gene age and functionality

Bacterial and archaeal genomes encompass numerous operons that typically consist of two to five genes. On larger scales, however, gene order is poorly conserved through the evolution of prokaryotes. Nevertheless, non-random localization of different classes of genes on prokaryotic chromosomes could reflect important functional and evolutionary constraints. We explored the patterns of genomic localization of evolutionarily conserved (ancient) and variable (young) genes across the diversity of bacteria and archaea. Nearly all bacterial and archaeal chromosomes were found to encompass large segments of 100-300 kilobases that were significantly enriched in either ancient or young genes. Similar clustering of genes with lethal knockout phenotype (essential genes) was observed as well. Mathematical modeling of genome evolution suggests that this long-range gene clustering in prokaryotic chromosomes reflects perpetual genome rearrangement driven by a combination of selective and neutral processes rather than evolutionary conservation.

Origin, evolution, and maintenance of gene-strand bias in bacteria

Gene-strand bias is a characteristic feature of bacterial genome organization wherein genes are preferentially encoded on the leading strand of replication, promoting co-orientation of replication and transcription. This co-orientation bias has evolved to protect gene essentiality, expression, and genomic stability from the harmful effects of head-on replication-transcription collisions. However, the origin, variation, and maintenance of gene-strand bias remain elusive. Here, we reveal that the frequency of inversions that alter gene orientation exhibits large variation across bacterial populations and negatively correlates with gene-strand bias. The density, distance, and distribution of inverted repeats show a similar negative relationship with gene-strand bias explaining the heterogeneity in inversions. Importantly, these observations are broadly evident across the entire bacterial kingdom uncovering inversions and inverted repeats as primary factors underlying the variation in gene-strand bias and its maintenance. The distinct catalytic subunits of replicative DNA polymerase have co-evolved with gene-strand bias, suggesting a close link between replication and the origin of gene-strand bias. Congruently, inversion frequencies and inverted repeats vary among bacteria with different DNA polymerases. In summary, we propose that the nature of replication determines the fitness cost of replication-transcription collisions, establishing a selection gradient on gene-strand bias by fine-tuning DNA sequence repeats and, thereby, gene inversions.

Unexplored power of CRISPR-Cas9 in neuroscience, a multi-OMICs review

The neuroscience and neurobiology of gene editing to enhance learning and memory is of paramount interest to the scientific community. The advancements of CRISPR system have created avenues to treat neurological disorders by means of versatile modalities varying from expression to suppression of genes and proteins. Neurodegenerative disorders have also been attributed to non-canonical DNA secondary structures by affecting neuron activity through controlling gene expression, nucleosome shape, transcription, translation, replication, and recombination. Changing DNA regulatory elements which could contribute to the fate and function of neurons are thoroughly discussed in this review. This study presents the ability of CRISPR system to boost learning power and memory, treat or cure genetically-based neurological disorders, and alleviate psychiatric diseases by altering the activity and the irritability of the neurons at the synaptic cleft through DNA manipulation, and also, epigenetic modifications using Cas9. We explore and examine how each different OMIC techniques can come useful when altering DNA sequences. Such insight into the underlying relationship between OMICs and cellular behaviors leads us to better neurological and psychiatric therapeutics by intelligently designing and utilizing the CRISPR/Cas9 technology.

Increasing transposase abundance with ocean depth correlates with a particle-associated lifestyle

Transposases are mobile genetic elements that move within and between genomes, promoting genomic plasticity in microorganisms. In marine microbial communities, the abundance of transposases increases with depth, but the reasons behind this trend remain unclear. Our analysis of metagenomes from the Tara Oceans and Malaspina Expeditions suggests that a particle-associated lifestyle is the main covariate for the high occurrence of transposases in the deep ocean, and this trend holds true for individual genomes as well as in a community-wide sense. We observed a strong and depth-independent correlation between transposase abundance and the presence of biofilm-associated genes, as well as the prevalence of secretory enzymes. This suggests that mobile genetic elements readily propagate among microbial communities within crowded biofilms. Furthermore, we show that particle association positively correlates with larger genome size, which is in turn associated with higher transposase abundance. Cassette sequences associated with transposons are enriched with genes related to defense mechanisms, which are more highly expressed in the deep sea. Thus, while transposons spread at the expense of their microbial hosts, they also introduce novel genes and potentially benefit the hosts in helping to compete for limited resources. Overall, our results suggest a new understanding of deep ocean particles as highways for gene sharing among defensively oriented microbial genomes.IMPORTANCEGenes can move within and between microbial genomes via mobile genetic elements, which include transposases and transposons. In the oceans, there is a puzzling increase in transposase abundance in microbial genomes as depth increases. To gain insight into this trend, we conducted an extensive analysis of marine microbial metagenomes and metatranscriptomes. We found a significant correlation between transposase abundance and a particle-associated lifestyle among marine microbes at both the metagenome and genome-resolved levels. We also observed a link between transposase abundance and genes related to defense mechanisms. These results suggest that as microbes become densely packed into crowded particles, mobile genes are more likely to spread and carry genetic material that provides a competitive advantage in crowded habitats. This may enable deep sea microbes to effectively compete in such environments.

A new mathematical model of folate homeostasis in E. coli highlights the potential importance of the folinic acid futile cycle in cell growth

Folate (vitamin B9) plays a central role in one-carbon metabolism in prokaryotes and eukaryotes. This pathway mediates the transfer of one-carbon units, playing a crucial role in nucleotide synthesis, methylation, and amino acid homeostasis. The folinic acid futile cycle adds a layer of intrigue to this pathway, due to its associations with metabolism, cell growth, and dormancy. It also introduces additional complexity to folate metabolism. A logical way to deal with such complexity is to examine it by using mathematical modelling. This work describes the construction and analysis of a model of folate metabolism, which includes the folinic acid futile cycle. This model was tested under three in silico growth conditions. Model simulations revealed: 1) the folate cycle behaved as a stable biochemical system in three growth states (slow, standard, and rapid); 2) the initial concentration of serine had the greatest impact on metabolite concentrations; 3) 5-formyltetrahydrofolate cyclo-ligase (5-FCL) activity had a significant impact on the levels of the 7 products that carry the one-carbon donated from folates, and the redox couple NADP/NADPH; this was particularly evident in the rapid growth state; 4) 5-FCL may be vital to the survival of the cells by maintaining low levels of homocysteine, as high levels can induce toxicity; and 5) the antifolate therapeutic trimethoprim had a greater impact on folate metabolism with higher nutrient availability. These results highlight the important role of 5-FCL in intracellular folate homeostasis and mass generation under different metabolic scenarios.

N-glycosylation in Archaea: Unusual sugars and unique modifications

Archaea are microorganisms that comprise a distinct branch of the universal tree of life and which are best known as extremophiles, residing in a variety of environments characterized by harsh physical conditions. One seemingly universal trait of Archaea is the ability to perform N-glycosylation. At the same time, archaeal N-linked glycans present variety in terms of both composition and architecture not seen in the parallel eukaryal or bacterial processes. In this mini-review, many of the unique and unusual sugars found in archaeal N-linked glycans as identified by nuclear magnetic resonance spectroscopy are described.

Using Machine Learning to Predict Genes Underlying Differentiation of Multipartite and Unipartite Traits in Bacteria

Since the discovery of the second chromosome in the Rhodobacter sphaeroides 2.4.1 by Suwanto and Kaplan in 1989 and the revelation of gene sequences, multipartite genomes have been reported in over three hundred bacterial species under nine different phyla. This phenomenon shattered the dogma of a unipartite genome (a single circular chromosome) in bacteria. Recently, Artificial Intelligence (AI), machine learning (ML), and Deep Learning (DL) have emerged as powerful tools in the investigation of big data in a plethora of disciplines to decipher complex patterns in these data, including the large-scale analysis and interpretation of genomic data. An important inquiry in bacteriology pertains to the genetic factors that underlie the structural evolution of multipartite and unipartite bacterial species. Towards this goal, here we have attempted to leverage machine learning as a means to identify the genetic factors that underlie the differentiation of, in general, bacteria with multipartite genomes and bacteria with unipartite genomes. In this study, deploying ML algorithms yielded two gene lists of interest: one that contains 46 discriminatory genes obtained following an assessment on all gene sets, and another that contains 35 discriminatory genes obtained based on an investigation of genes that are differentially present (or absent) in the genomes of the multipartite bacteria and their respective close relatives. Our study revealed a small pool of genes that discriminate bacteria with multipartite genomes and their close relatives with single-chromosome genomes. Machine learning thus aided in uncovering the genetic factors that underlie the differentiation of bacterial multipartite and unipartite traits.

Studying the Human Microbiota: Advances in Understanding the Fundamentals, Origin, and Evolution of Biological Timekeeping

The recently observed circadian oscillations of the intestinal microbiota underscore the profound nature of the human-microbiome relationship and its importance for health. Together with the discovery of circadian clocks in non-photosynthetic gut bacteria and circadian rhythms in anucleated cells, these findings have indicated the possibility that virtually all microorganisms may possess functional biological clocks. However, they have also raised many essential questions concerning the fundamentals of biological timekeeping, its evolution, and its origin. This narrative review provides a comprehensive overview of the recent literature in molecular chronobiology, aiming to bring together the latest evidence on the structure and mechanisms driving microbial biological clocks while pointing to potential applications of this knowledge in medicine. Moreover, it discusses the latest hypotheses regarding the evolution of timing mechanisms and describes the functions of peroxiredoxins in cells and their contribution to the cellular clockwork. The diversity of biological clocks among various human-associated microorganisms and the role of transcriptional and post-translational timekeeping mechanisms are also addressed. Finally, recent evidence on metabolic oscillators and host-microbiome communication is presented.

Characterization of the pig lower respiratory tract antibiotic resistome

Respiratory diseases and its treatments are highly concerned in both the pig industry and human health. However, the composition and distribution of antibiotic resistance genes (ARGs) in swine lower respiratory tract microbiome remain unknown. The relationships of ARGs with mobile genetic elements (MGEs) and lung health are unclear. Here, we characterize antibiotic resistomes of the swine lower respiratory tract microbiome containing 1228 open reading frames belonging to 372 ARGs using 745 metagenomes from 675 experimental pigs. Twelve ARGs conferring resistance to tetracycline are related to an MGE Tn916 family, and multiple types of ARGs are related to a transposase gene tnpA. Most of the linkage complexes between ARGs and MGEs (the Tn916 family and tnpA) are also observed in pig gut microbiomes and human lung microbiomes, suggesting the high risk of these MGEs mediating ARG transfer to both human and pig health. Gammaproteobacteria are the major ARG carriers, within which Escherichia coli harbored >50 ARGs and >10 MGEs. Although the microbial compositions structure the compositions of ARGs, we identify 73 ARGs whose relative abundances are significantly associated with the severity of lung lesions. Our results provide the first overview of ARG profiles in the swine lower respiratory tract microbiome.

Reconstruction of the last bacterial common ancestor from 183 pangenomes reveals a versatile ancient core genome

BACKGROUND

Cumulative sequencing efforts have yielded enough genomes to construct pangenomes for dozens of bacterial species and elucidate intraspecies gene conservation. Given the diversity of organisms for which this is achievable, similar analyses for ancestral species are feasible through the integration of pangenomics and phylogenetics, promising deeper insights into the nature of ancient life.

RESULTS

We construct pangenomes for 183 bacterial species from 54,085 genomes and identify their core genomes using a novel statistical model to estimate genome-specific error rates and underlying gene frequencies. The core genomes are then integrated into a phylogenetic tree to reconstruct the core genome of the last bacterial common ancestor (LBCA), yielding three main results: First, the gene content of modern and ancestral core genomes are diverse at the level of individual genes but are similarly distributed by functional category and share several poorly characterized genes. Second, the LBCA core genome is distinct from any individual modern core genome but has many fundamental biological systems intact, especially those involving translation machinery and biosynthetic pathways to all major nucleotides and amino acids. Third, despite this metabolic versatility, the LBCA core genome likely requires additional non-core genes for viability, based on comparisons with the minimal organism, JCVI-Syn3A.

CONCLUSIONS

These results suggest that many cellular systems commonly conserved in modern bacteria were not just present in ancient bacteria but were nearly immutable with respect to short-term intraspecies variation. Extending this analysis to other domains of life will likely provide similar insights into more distant ancestral species.

Large-scale phylogenomics of aquatic bacteria reveal molecular mechanisms for adaptation to salinity

The crossing of environmental barriers poses major adaptive challenges. Rareness of freshwater-marine transitions separates the bacterial communities, but how these are related to brackish counterparts remains elusive, as do the molecular adaptations facilitating cross-biome transitions. We conducted large-scale phylogenomic analysis of freshwater, brackish, and marine quality-filtered metagenome-assembled genomes (11,248). Average nucleotide identity analyses showed that bacterial species rarely existed in multiple biomes. In contrast, distinct brackish basins cohosted numerous species, but their intraspecific population structures displayed clear signs of geographic separation. We further identified the most recent cross-biome transitions, which were rare, ancient, and most commonly directed toward the brackish biome. Transitions were accompanied by systematic changes in amino acid composition and isoelectric point distributions of inferred proteomes, which evolved over millions of years, as well as convergent gains or losses of specific gene functions. Therefore, adaptive challenges entailing proteome reorganization and specific changes in gene content constrains the cross-biome transitions, resulting in species-level separation between aquatic biomes.

Comparative genomics of Stutzerimonas balearica (Pseudomonas balearica): diversity, habitats, and biodegradation of aromatic compounds

Stutzerimonas balearica (Pseudomonas balearica) has been found principally in oil-polluted environments. The capability of S. balearica to thrive from the degradation of pollutant compounds makes it a species of interest for potential bioremediation applications. However, little has been reported about the diversity of S. balearica. In this study, genome sequences of S. balearica strains from different origins were analyzed, revealing that it is a diverse species with an open pan-genome that will continue revealing new genes and functionalities as the genomes of more strains are sequenced. The nucleotide signatures and intra- and inter-species variation of the 16S rRNA genes of S. balearica were reevaluated. A strategy of screening 16S rRNA gene sequences in public databases enabled the detection of 158 additional strains, of which only 23% were described as S. balearica. The species was detected from a wide range of environments, although mostly from aquatic and polluted environments, predominantly related to petroleum oil. Genomic and phenotypic analyses confirmed that S. balearica possesses varied inherent capabilities for aromatic compounds degradation. This study increases the knowledge of the biology and diversity of S. balearica and will serve as a basis for future work with the species.

Bacteria and Archaea Regulate Particulate Organic Matter Export in Suspended and Sinking Marine Particle Fractions

The biological carbon pump (BCP) in the Southern Ocean is driven by phytoplankton productivity and is a significant organic matter sink. However, the role of particle-attached (PA) and free-living (FL) prokaryotes (bacteria and archaea) and their diversity in influencing the efficiency of the BCP is still unclear. To investigate this, we analyzed the metagenomes linked to suspended and sinking marine particles from the Sub-Antarctic Southern Ocean Time Series (SOTS) by deploying a Marine Snow Catcher (MSC), obtaining suspended and sinking particulate material, determining organic carbon and nitrogen flux, and constructing metagenome-assembled genomes (MAGs). The suspended and sinking particle-pools were dominated by bacteria with the potential to degrade organic carbon. Bacterial communities associated with the sinking fraction had more genes related to the degradation of complex organic carbon than those in the suspended fraction. Archaea had the potential to drive nitrogen metabolism via nitrite and ammonia oxidation, altering organic nitrogen concentration. The data revealed several pathways for chemoautotrophy and the secretion of recalcitrant dissolved organic carbon (RDOC) from CO₂, with bacteria and archaea potentially sequestering particulate organic matter (POM) via the production of RDOC. These findings provide insights into the diversity and function of prokaryotes in suspended and sinking particles and their role in organic carbon/nitrogen export in the Southern Ocean. IMPORTANCE The biological carbon pump is crucial for the export of particulate organic matter in the ocean. Recent studies on marine microbes have shown the profound influence of bacteria and archaea as regulators of particulate organic matter export. Yet, despite the importance of the Southern Ocean as a carbon sink, we lack comparable insights regarding microbial contributions. This study provides the first insights regarding prokaryotic contributions to particulate organic matter export in the Southern Ocean. We reveal evidence that prokaryotic communities in suspended and sinking particle fractions harbor widespread genomic potential for mediating particulate organic matter export. The results substantially enhance our understanding of the role played by microorganisms in regulating particulate organic matter export in suspended and sinking marine fractions in the Southern Ocean.

Phylogroup-specific variation shapes the clustering of antimicrobial resistance genes and defence systems across regions of genome plasticity in Pseudomonas aeruginosa

BACKGROUND

Pseudomonas aeruginosa is an opportunistic pathogen consisting of three phylogroups (hereafter named A, B, and C). Here, we assessed phylogroup-specific evolutionary dynamics across available and also new P. aeruginosa genomes.

METHODS

In this genomic analysis, we first generated new genome assemblies for 18 strains of the major P. aeruginosa clone type (mPact) panel, comprising a phylogenetically diverse collection of clinical and environmental isolates for this species. Thereafter, we combined these new genomes with 1991 publicly available P. aeruginosa genomes for a phylogenomic and comparative analysis. We specifically explored to what extent antimicrobial resistance (AMR) genes, defence systems, and virulence genes vary in their distribution across regions of genome plasticity (RGPs) and "masked" (RGP-free) genomes, and to what extent this variation differs among the phylogroups.

FINDINGS

We found that members of phylogroup B possess larger genomes, contribute a comparatively larger number of pangenome families, and show lower abundance of CRISPR-Cas systems. Furthermore, AMR and defence systems are pervasive in RGPs and integrative and conjugative/mobilizable elements (ICEs/IMEs) from phylogroups A and B, and the abundance of these cargo genes is often significantly correlated. Moreover, inter- and intra-phylogroup interactions occur at the accessory genome level, suggesting frequent recombination events. Finally, we provide here the mPact panel of diverse P. aeruginosa strains that may serve as a valuable reference for functional analyses.

INTERPRETATION

Altogether, our results highlight distinct pangenome characteristics of the P. aeruginosa phylogroups, which are possibly influenced by variation in the abundance of CRISPR-Cas systems and are shaped by the differential distribution of other defence systems and AMR genes.

FUNDING

German Science Foundation, Max-Planck Society, Leibniz ScienceCampus Evolutionary Medicine of the Lung, BMBF program Medical Infection Genomics, Kiel Life Science Postdoc Award.

Abiotic selection of microbial genome size in the global ocean

Strong purifying selection is considered a major evolutionary force behind small microbial genomes in the resource-poor photic ocean. However, very little is currently known about how the size of prokaryotic genomes evolves in the global ocean and whether patterns reflect shifts in resource availability in the epipelagic and relatively stable deep-sea environmental conditions. Using 364 marine microbial metagenomes, we investigate how the average genome size of uncultured planktonic prokaryotes varies across the tropical and polar oceans to the hadal realm. We find that genome size is highest in the perennially cold polar ocean, reflecting elongation of coding genes and gene dosage effects due to duplications in the interior ocean microbiome. Moreover, the rate of change in genome size due to temperature is 16-fold higher than with depth up to 200 m. Our results demonstrate how environmental factors can influence marine microbial genome size selection and ecological strategies of the microbiome.

Analysis of multipartite bacterial genomes using alignment free and alignment-based pipelines

Since the discovery of second chromosome in Rhodobacter sphaeroides 2.4.1 in 1989, multipartite genomes have been reported in over three hundred bacterial species under nine different phyla. This has shattered the unipartite (single chromosome) genome dogma in bacteria. Since then, many questions on various aspects of multipartite genomes in bacteria have been addressed. However, our understanding of how multipartite genomes emerge and evolve is still lacking. Importantly, the knowledge of genetic factors underlying the differences in multipartite and single-chromosome genomes is lacking. In this work, we have performed comparative evolutionary and functional genomics analyses to identify molecular factors that discriminate multipartite from unipartite bacteria, with the goal to decipher taxon-specific factors, and those that are prevalent across the taxa, underlying these traits. We assessed the roles of evolutionary mechanisms, specifically gene gain, in driving the divergence of bacteria with single and multiple chromosomes. In addition, we performed functional genomic analysis to garner support for our findings from comparative evolutionary analysis. We found genes such as those encoding conserved hypothetical proteins in Deinococcus radiodurans R1, and putative phage phi-C31 gp36 major capsid like and hypothetical proteins in Rhodobacter sphaeroides 2.4.1, which are located on accessory chromosomes in these bacteria but were not found in the inferred ancestral sequences, and on the primary chromosomes, as well as were not found in their closest relatives with single chromosome within the same clade. Our study shines a new light on the potential roles of the secondary chromosomes in helping bacteria with multipartite genomes to adapt to specialized environments or growth conditions.

Transperons: RNA operons as effectors of coordinated gene expression in eukaryotes

Coordinated gene expression allows spatiotemporal control of cellular processes and is achieved by the cotranscription/translation of functionally related genes/proteins. Prokaryotes evolved polycistronic messages (operons) to confer expression from a single promoter to efficiently cotranslate proteins functioning on the same pathway. Yet, despite having far greater diversity (e.g., gene number, distribution, modes of expression), eukaryotic cells employ individual promoters and monocistronic messages. Although gene expression is modular, it does not account for how eukaryotes achieve coordinated localized translation. The RNA operon theory states that mRNAs derived from different chromosomes assemble into ribonucleoprotein particles (RNPs) that act as functional operons to generate protein cohorts upon cotranslation. Work in yeast has now validated this theory and shown that intergenic associations and noncanonical histone functions create pathway-specific RNA operons (transperons) that regulate cell physiology. Herein the involvement of chromatin organization in transperon formation and programmed gene coexpression is discussed.

Patterns of pan-genome occupancy and gene coexpression under water-deficit in Brachypodium distachyon

Natural populations are characterized by abundant genetic diversity driven by a range of different types of mutation. The tractability of sequencing complete genomes has allowed new insights into the variable composition of genomes, summarized as a species pan-genome. These analyses demonstrate that many genes are absent from the first reference genomes, whose analysis dominated the initial years of the genomic era. Our field now turns towards understanding the functional consequence of these highly variable genomes. Here, we analysed weighted gene coexpression networks from leaf transcriptome data for drought response in the purple false brome Brachypodium distachyon and the differential expression of genes putatively involved in adaptation to this stressor. We specifically asked whether genes with variable "occupancy" in the pan-genome - genes which are either present in all studied genotypes or missing in some genotypes - show different distributions among coexpression modules. Coexpression analysis united genes expressed in drought-stressed plants into nine modules covering 72 hub genes (87 hub isoforms), and genes expressed under controlled water conditions into 13 modules, covering 190 hub genes (251 hub isoforms). We find that low occupancy pan-genes are under-represented among several modules, while other modules are over-enriched for low-occupancy pan-genes. We also provide new insight into the regulation of drought response in B. distachyon, specifically identifying one module with an apparent role in primary metabolism that is strongly responsive to drought. Our work shows the power of integrating pan-genomic analysis with transcriptomic data using factorial experiments to understand the functional genomics of environmental response.

In silico approach to identify microsatellite candidate biomarkers to differentiate the biovar of Corynebacterium pseudotuberculosis genomes

Corynebacterium pseudotuberculosis is the causative bacterial agent of the zoonotic disease known as caseous lymphadenitis, and it presents several mechanisms of response to host defenses, including the presence of virulence factors (VFs). The genomes of these bacteria have several polymorphic markers known as microsatellites, or simple sequence repeats (SSRs), that can be used to characterize the genome, to study possible polymorphisms existing among strains, and to verify the effects of such polymorphic markers in coding regions and regions associated with VFs. In this study, several SSRs were identified within coding regions throughout the 54 genomes of this species, revealing possible polymorphisms associated with coding regions that could be used as strain-specific or serotype-specific identifiers of C. pseudotuberculosis. The similarities associated with SSRs amongst the different serum variants of C. pseudotuberculosis, biovars equi and ovis, were also evaluated, and it was possible to identify SSRs located in coding regions responsible for a VF enrolled in pathogenesis known to mediate bacterial adherence (SpaH-type pili virulence factor). Phylogenetic analyses revealed that strains sharing SSR patterns, including the possible polymorphisms identified in the same position of gene-coding regions, were displayed by strains with a common ancestor, corroborating with the Genome Tree Report of the NCBI. Statistical analysis showed that the microsatellite groups belonging to equi and ovis biovars have a significance of 0.006 (p-value) in similarity, thus indicating them as good biomarker candidates for C. pseudotuberculosis.

Nutrition or nature: using elementary flux modes to disentangle the complex forces shaping prokaryote pan-genomes

Background

Microbial pan-genomes are shaped by a complex combination of stochastic and deterministic forces. Even closely related genomes exhibit extensive variation in their gene content. Understanding what drives this variation requires exploring the interactions of gene products with each other and with the organism’s external environment. However, to date, conceptual models of pan-genome dynamics often represent genes as independent units and provide limited information about their mechanistic interactions.

Results

We simulated the stochastic process of gene-loss using the pooled genome-scale metabolic reaction networks of 46 taxonomically diverse bacterial and archaeal families as proxies for their pan-genomes. The frequency by which reactions are retained in functional networks when stochastic gene loss is simulated in diverse environments allowed us to disentangle the metabolic reactions whose presence depends on the metabolite composition of the external environment (constrained by “nutrition”) from those that are independent of the environment (constrained by “nature”). By comparing the frequency of reactions from the first group with their observed frequencies in bacterial and archaeal families, we predicted the metabolic niches that shaped the genomic composition of these lineages. Moreover, we found that the lineages that were shaped by a more diverse metabolic niche also occur in more diverse biomes as assessed by global environmental sequencing datasets.

Conclusion

We introduce a computational framework for analyzing and interpreting pan-reactomes that provides novel insights into the ecological and evolutionary drivers of pan-genome dynamics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02052-3.

In Silico Study of Mercury Resistance Genes Extracted from Pseudomonas spp. Involved in Bioremediation: Understanding the Promoter Regions and Regulatory Elements

Microbial genes and their product were diverse and beneficial for heavy metal bioremediation from the contaminated sites. Screening of genes and gene products plays a significant role in the detoxification of pollutants. Understanding of the promoter region and its regulatory elements is a vital implication of microbial genes. To the best of our knowledge, there is no in silico study reported so far on mer gene families used for heavy metal bioremediation. The motif distribution was observed densely upstream of the TSSs (transcription start sites) between +1 and -350 bp and sparsely distributed above -350 bp, according to the current study. MEME identified the best common candidate motifs of TFs (transcription factors) binding with the lowest e value (7.2 $e$ -033) and is the most statistically significant candidate motif. The EXPREG output of the 11 TFs with varying degrees of function such as activation, repression, transcription, and dual purposes was thoroughly examined. Data revealed that transcriptional gene regulation in terms of activation and repression was observed at 36.4% and 54.56%, respectively. This shows that most TFs are involved in transcription gene repression rather than activation. Likewise, EXPREG output revealed that transcriptional conformational modes, such as monomers, dimers, tetramers, and other factors, were also analyzed. The data indicated that most of the transcriptional conformation mode was dual, which accounts for 96%. CpG island analysis using online and offline tools revealed that the gene body had fewer CpG islands compared to the promoter regions. Understanding the common candidate motifs, transcriptional factors, and regulatory elements of the mer operon gene cluster using a machine learning approach could help us better understand gene expression patterns in heavy metal bioremediation.

Microbial genomic trait evolution is dominated by frequent and rare pulsed evolution

On the macroevolutionary time scale, does trait evolution proceed gradually or by rapid bursts (pulses) separated by prolonged periods of stasis or slow evolution? Although studies have shown that pulsed evolution is prevalent in animals, our knowledge about the tempo and mode of evolution across the tree of life is very limited. This long-standing debate calls for a test in bacteria and archaea, the most ancient and diverse forms of life with unique population genetic properties. Using a likelihood-based framework, we show that pulsed evolution is not only present but also prevalent and predominant in microbial genomic trait evolution. We detected two distinct types of pulsed evolution (small frequent and large rare jumps) that are predicted by the punctuated equilibrium and quantum evolution theories. Our findings suggest that major bacterial lineages could have originated in quick bursts and that pulsed evolution is a common theme across the tree of life.

Realistic Gene Transfer to Gene Duplication Ratios Identify Different Roots in the Bacterial Phylogeny Using a Tree Reconciliation Method

The rooting of phylogenetic trees permits important inferences about ancestral states and the polarity of evolutionary events. Recently, methods that reconcile discordance between gene-trees and species-trees—tree reconciliation methods—are becoming increasingly popular for rooting species trees. Rooting via reconciliation requires values for a particular parameter, the gene transfer to gene duplication ratio (T:D), which in current practice is estimated on the fly from discordances observed in the trees. To date, the accuracy of T:D estimates obtained by reconciliation analyses has not been compared to T:D estimates obtained by independent means, hence the effect of T:D upon inferences of species tree roots is altogether unexplored. Here we investigated the issue in detail by performing tree reconciliations of more than 10,000 gene trees under a variety of T:D ratios for two phylogenetic cases: a bacterial (prokaryotic) tree with 265 species and a fungal-metazoan (eukaryotic) tree with 31 species. We show that the T:D ratios automatically estimated by a current tree reconciliation method, ALE, generate virtually identical T:D ratios across bacterial genes and fungal-metazoan genes. The T:D ratios estimated by ALE differ 10- to 100-fold from robust, ALE-independent estimates from real data. More important is our finding that the root inferences using ALE in both datasets are strongly dependent upon T:D. Using more realistic T:D ratios, the number of roots inferred by ALE consistently increases and, in some cases, clearly incorrect roots are inferred. Furthermore, our analyses reveal that gene duplications have a far greater impact on ALE’s preferences for phylogenetic root placement than gene transfers or gene losses do. Overall, we show that obtaining reliable species tree roots with ALE is only possible when gene duplications are abundant in the data and the number of falsely inferred gene duplications is low. Finding a sufficient sample of true gene duplications for rooting species trees critically depends on the T:D ratios used in the analyses. T:D ratios, while being important parameters of genome evolution in their own right, affect the root inferences with tree reconciliations to an unanticipated degree.

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 to core genes, accessory genes are enriched in genome plasticity functions. We found no evidence for recent horizontal gene transfer 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.​.

Genomic analysis provides novel insights into diversification and taxonomy of Allorhizobium vitis (i.e. Agrobacterium vitis)

BACKGROUND

Allorhizobium vitis (formerly named Agrobacterium vitis or Agrobacterium biovar 3) is the primary causative agent of crown gall disease of grapevine worldwide. We obtained and analyzed whole-genome sequences of diverse All. vitis strains to get insights into their diversification and taxonomy.

RESULTS

Pairwise genome comparisons and phylogenomic analysis of various All. vitis strains clearly indicated that All. vitis is not a single species, but represents a species complex composed of several genomic species. Thus, we emended the description of All. vitis, which now refers to a restricted group of strains within the All. vitis species complex (i.e. All. vitis sensu stricto) and proposed a description of a novel species, All. ampelinum sp. nov. The type strain of All. vitis sensu stricto remains the current type strain of All. vitis, K309T. The type strain of All. ampelinum sp. nov. is S4T. We also identified sets of gene clusters specific to the All. vitis species complex, All. vitis sensu stricto and All. ampelinum, respectively, for which we predicted the biological function and infer the role in ecological diversification of these clades, including some we could experimentally validate. All. vitis species complex-specific genes confer tolerance to different stresses, including exposure to aromatic compounds. Similarly, All. vitis sensu stricto-specific genes confer the ability to degrade 4-hydroxyphenylacetate and a putative compound related to gentisic acid. All. ampelinum-specific genes have putative functions related to polyamine metabolism and nickel assimilation. Congruently with the genome-based classification, All. vitis sensu stricto and All. ampelinum were clearly delineated by MALDI-TOF MS analysis. Moreover, our genome-based analysis indicated that Allorhizobium is clearly separated from other genera of the family Rhizobiaceae.

CONCLUSIONS

Comparative genomics and phylogenomic analysis provided novel insights into the diversification and taxonomy of Allorhizobium vitis species complex, supporting our redefinition of All. vitis sensu stricto and description of All. ampelinum. Our pan-genome analyses suggest that these species have differentiated ecologies, each relying on specialized nutrient consumption or toxic compound degradation to adapt to their respective niche.

Genomic attributes of thermophilic and hyperthermophilic bacteria and archaea

Thermophiles and hyperthermophiles are immensely useful in understanding the evolution of life, besides their utility in environmental and industrial biotechnology. Advancements in sequencing technologies have revolutionized the field of microbial genomics. The massive generation of data enhances the sequencing coverage multi-fold and allows to analyse the entire genomic features of microbes efficiently and accurately. The mandate of a pure isolate can also be bypassed where whole metagenome-assembled genomes and single cell-based sequencing have fulfilled the majority of the criteria to decode various attributes of microbial genomes. A boom has, therefore, been seen in analysing the extremophilic bacteria and archaea using sequence-based approaches. Due to extensive sequence analysis, it becomes easier to understand the gene flow and their evolution among the members of bacteria and archaea. For instance, sequencing unveiled that Thermotoga maritima shares around 24% of genes of archaeal origin. Comparative and functional genomics provide an analytical view to understanding the microbial diversity of thermophilic bacteria and archaea, their interactions with other microbes, their adaptations, gene flow, and evolution over time. In this review, the genomic features of thermophilic bacteria and archaea are dealt with comprehensively.

Genomic Insights into Achromobacter mucicolens IA Antibiotic Resistance

Achromobacter denitrificans is an environmental opportunistic pathogen that is infecting a large number of immunocompromised patients. A more recently identified strain from the historical collection of strains of Achromobacter denitrificans is Achromobacter mucicolens. In hosts with a variety of underlying diseases, Achromobacter spp. can induce a wide spectrum of disorders. Because of the bacterium’s intrinsic genetic constitution and resistance gained over time, antibiotics are challenged to handle A. mucicolens. Due to the fact that A. mucicolens is rare and its taxonomy is not completely understood, it is difficult to define clinical symptoms, acquisition risk factors, and thus the best therapeutic course of action. To help comprehend this intrinsic and acquired resistance, we analyzed the entire genome of the A. mucicolens IA strain and utilized bioinformatics methods to estimate the strain's probable drug resistance profile. In our study, we have isolated and cultured a clinically important A. mucicolens strain and subjected it to antimicrobial susceptibility tests against antibiotics in the Vitek 2 testing system. The strain’s genome sequence as well as an investigation of 27 of its phenotypic traits provides important information regarding this pathogen. The genome of this A. mucicolens IA strain possesses a number of antibiotic resistance genes that code for efflux pump systems and other antibiotic-regulating as well as -modifying enzymes. Our research analysis predicted genes involved in drug resistance, including genes for efflux pump systems, antibiotic efflux, antibiotic inactivation, and antibiotic target alteration. In vitro studies validated the genomic evidence for its ability to exhibit resistance against a wide range of antibiotics. Our investigation paves the way for more research on understanding the functioning of the key discovered genes that contribute toward the pathogenicity of A. mucicolens and hence gives new information and treatment options for this emerging pathogen.

IMPORTANCEAchromobacter species are well-known opportunistic human pathogens that can be found in water and soil and most commonly in hospital settings. They thrive in immunocompromised individuals, producing sporadic cases of pneumonia, septicemia, peritonitis, urinary tract infections, and other illnesses. Achromobacter strains are inherently resistant to a wide spectrum of antibiotics, making them difficult to treat promptly. The strain under study, A. mucicolens, was notably resistant to various antibiotics, and the infection could be controlled only after several rounds of prescription medications at different doses. This consumed a lot of time and put the already immunosuppressed leukemic patient through a great ordeal. The study aimed to raise awareness about the importance of the Achromobacter bacterium’s lethality, and doctors should evaluate the bacterium’s potential for resistance before prescribing antibiotics. Sanitation and other precautions should also be implemented in hospitals and other public places.

Identification of Protein Drug Targets of Biofilm Formation and Quorum Sensing in Multidrug Resistant Enterococcus faecalis

Abstract:

Enterococcus faecalis (E. faecalis) is an opportunistic multidrug-resistant (MDR) pathogen found in the guts of humans and farmed animals. Due to the occurrence of (MDR) strain there is an urgent need to look for an alternative treatment approach. E. faecalis is a Gram-positive bacterium, which is among the most prevalent multidrug resistant hospital pathogens. Its ability to develop quorum sensing (QS) mediated biofilm formation further exacerbates the pathogenicity and triggers lifethreatening infections. Therefore, developing a suitable remedy for curing E. faecalis mediated enterococcal infections is an arduous task. Several putative virulence factors and proteins are involved in the development of biofilms in E. faecalis. Such proteins often play important roles in virulence, disease, and colonization by pathogens. The elucidation of the structure-function relationship of such protein drug targets and the interacting compounds could provide an attractive paradigm towards developing structure-based drugs against E. faecalis. This review provides a comprehensive overview of the current status, enigmas that warrant further studies, and the prospects toward alleviating the antibiotic resistance in E. faecalis. Specifically, the role of biofilm and quorum sensing (QS) in the emergence of MDR strains had been elaborated along with the importance of the protein drug targets involved in both the processes.

Toward a theory of evolution as multilevel learning

We apply the theory of learning to physically renormalizable systems in an attempt to outline a theory of biological evolution, including the origin of life, as multilevel learning. We formulate seven fundamental principles of evolution that appear to be necessary and sufficient to render a universe observable and show that they entail the major features of biological evolution, including replication and natural selection. It is shown that these cornerstone phenomena of biology emerge from the fundamental features of learning dynamics such as the existence of a loss function, which is minimized during learning. We then sketch the theory of evolution using the mathematical framework of neural networks, which provides for detailed analysis of evolutionary phenomena. To demonstrate the potential of the proposed theoretical framework, we derive a generalized version of the Central Dogma of molecular biology by analyzing the flow of information during learning (back propagation) and predicting (forward propagation) the environment by evolving organisms. The more complex evolutionary phenomena, such as major transitions in evolution (in particular, the origin of life), have to be analyzed in the thermodynamic limit, which is described in detail in the paper by Vanchurin et al. [V. Vanchurin, Y. I. Wolf, E. V. Koonin, M. I. Katsnelson, Proc. Natl. Acad. Sci. U.S.A. 119, 10.1073/pnas.2120042119 (2022)].

Ecological selection of bacterial taxa with larger genome sizes in response to polycyclic aromatic hydrocarbons stress

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous priority pollutants that cause great damage to the natural environment and health. Average genome size in a community is critical for shedding light on microbiome's functional response to pollution stress within an environment. Here, microcosms under different concentrations were performed to evaluate the selection of PAHs stress on the average genome size in a community. We found the distinct communities of significantly larger genome size with the increase of PAHs concentration gradients in soils, and consistent trends were discovered in soils at different latitudes. The abundance of Proteobacteria and Deinococcus-Thermus with relatively larger genomes increased along with PAHs stress and well adapted to polluted environments. In contrast, the abundance of Patescibacteria with a highly streamlined and smaller genome decreased, implying complex interactions between environmental selection and functional fitness resulted in bacteria with larger genomes becoming more abundant. Moreover, we confirmed the increased capacity for horizontal transfer of degrading genes between communities by showing an increased connection number per node positively related to the nidA gene along the concentration gradients in the co-occurrence network. Our findings suggest PAHs tend to select bacterial taxa with larger genome sizes, with significant consequences for community stability and potential biodegradation strategies.

Prokaryote pangenomes are dynamic entities

Prokaryote pangenomes are influenced heavily by environmental factors and the opportunity for gene gain and loss events. As the field of pangenome analysis has expanded, so has the need to fully understand the complexity of how eco-evolutionary dynamics shape pangenomes. Here, we describe current models of pangenome evolution and discuss their suitability and accuracy. We suggest that pangenomes are dynamic entities under constant flux, highlighting the influence of two-way interactions between pangenome and environment. New classifications of core and accessory genes are also considered, underscoring the need for continuous evaluation of nomenclature in a fast-moving field. We conclude that future models of pangenome evolution should incorporate eco-evolutionary dynamics to fully encompass their dynamic, changeable nature.

Transposable elements promote the evolution of genome streamlining

Eukaryotes and prokaryotes have distinct genome architectures, with marked differences in genome size, the ratio of coding/non-coding DNA, and the abundance of transposable elements (TEs). As TEs replicate independently of their hosts, the proliferation of TEs is thought to have driven genome expansion in eukaryotes. However, prokaryotes also have TEs in intergenic spaces, so why do prokaryotes have small, streamlined genomes? Using an in silico model describing the genomes of single-celled asexual organisms that coevolve with TEs, we show that TEs acquired from the environment by horizontal gene transfer can promote the evolution of genome streamlining. The process depends on local interactions and is underpinned by rock-paper-scissors dynamics in which populations of cells with streamlined genomes beat TEs, which beat non-streamlined genomes, which beat streamlined genomes, in continuous and repeating cycles. Streamlining is maladaptive to individual cells, but improves lineage viability by hindering the proliferation of TEs. Streamlining does not evolve in sexually reproducing populations because recombination partially frees TEs from the deleterious effects they cause. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Comparative Genomics of Pseudomonas stutzeri Complex: Taxonomic Assignments and Genetic Diversity

Pseudomonas stutzeri is a species complex with extremely broad phenotypic and genotypic diversity. However, very little is known about its diversity, taxonomy and phylogeny at the genomic scale. To address these issues, we systematically and comprehensively defined the taxonomy and nomenclature for this species complex and explored its genetic diversity using hundreds of sequenced genomes. By combining average nucleotide identity (ANI) evaluation and phylogenetic inference approaches, we identified 123 P. stutzeri complex genomes covering at least six well-defined species among all sequenced Pseudomonas genomes; of these, 25 genomes represented novel members of this species complex. ANI values of ≥∼95% and digital DNA-DNA hybridization (dDDH) values of ≥∼60% in combination with phylogenomic analysis consistently and robustly supported the division of these strains into 27 genomovars (most likely species to some extent), comprising 16 known and 11 unknown genomovars. We revealed that 12 strains had mistaken taxonomic assignments, while 16 strains without species names can be assigned to the species level within the species complex. We observed an open pan-genome of the P. stutzeri complex comprising 13,261 gene families, among which approximately 45% gene families do not match any sequence present in the COG database, and a large proportion of accessory genes. The genome contents experienced extensive genetic gain and loss events, which may be one of the major mechanisms driving diversification within this species complex. Surprisingly, we found that the ectoine biosynthesis gene cluster (ect) was present in all genomes of P. stutzeri species complex strains but distributed at very low frequency (43 out of 9548) in other Pseudomonas genomes, suggesting a possible origin of the ancestors of P. stutzeri species complex in high-osmolarity environments. Collectively, our study highlights the potential of using whole-genome sequences to re-evaluate the current definition of the P. stutzeri complex, shedding new light on its genomic diversity and evolutionary history.

Viroids and Viroid-like Circular RNAs: Do They Descend from Primordial Replicators?

Viroids are a unique class of plant pathogens that consist of small circular RNA molecules, between 220 and 450 nucleotides in size. Viroids encode no proteins and are the smallest known infectious agents. Viroids replicate via the rolling circle mechanism, producing multimeric intermediates which are cleaved to unit length either by ribozymes formed from both polarities of the viroid genomic RNA or by coopted host RNAses. Many viroid-like small circular RNAs are satellites of plant RNA viruses. Ribozyviruses, represented by human hepatitis delta virus, are larger viroid-like circular RNAs that additionally encode the viral nucleocapsid protein. It has been proposed that viroids are direct descendants of primordial RNA replicons that were present in the hypothetical RNA world. We argue, however, that much later origin of viroids, possibly, from recently discovered mobile genetic elements known as retrozymes, is a far more parsimonious evolutionary scenario. Nevertheless, viroids and viroid-like circular RNAs are minimal replicators that are likely to be close to the theoretical lower limit of replicator size and arguably comprise the paradigm for replicator emergence. Thus, although viroid-like replicators are unlikely to be direct descendants of primordial RNA replicators, the study of the diversity and evolution of these ultimate genetic parasites can yield insights into the earliest stages of the evolution of life.

Comparative pangenomics: analysis of 12 microbial pathogen pangenomes reveals conserved global structures of genetic and functional diversity

BACKGROUND

With the exponential growth of publicly available genome sequences, pangenome analyses have provided increasingly complete pictures of genetic diversity for many microbial species. However, relatively few studies have scaled beyond single pangenomes to compare global genetic diversity both within and across different species. We present here several methods for "comparative pangenomics" that can be used to contextualize multi-pangenome scale genetic diversity with gene function for multiple species at multiple resolutions: pangenome shape, genes, sequence variants, and positions within variants.

RESULTS

Applied to 12,676 genomes across 12 microbial pathogenic species, we observed several shared resolution-specific patterns of genetic diversity: First, pangenome openness is associated with species' phylogenetic placement. Second, relationships between gene function and frequency are conserved across species, with core genomes enriched for metabolic and ribosomal genes and accessory genomes for trafficking, secretion, and defense-associated genes. Third, genes in core genomes with the highest sequence diversity are functionally diverse. Finally, certain protein domains are consistently mutation enriched across multiple species, especially among aminoacyl-tRNA synthetases where the extent of a domain's mutation enrichment is strongly function-dependent.

CONCLUSIONS

These results illustrate the value of each resolution at uncovering distinct aspects in the relationship between genetic and functional diversity across multiple species. With the continued growth of the number of sequenced genomes, these methods will reveal additional universal patterns of genetic diversity at the pangenome scale.

OUP accepted manuscript

Operons are a hallmark of the genomic and regulatory architecture of prokaryotes. However, the mechanism by which two genes placed far apart gradually come close and form operons remains to be elucidated. Here, we propose a new model of the origin of operons: Mobile genetic elements called insertion sequences can facilitate the formation of operons by consecutive insertion–deletion–excision reactions. This mechanism barely leaves traces of insertion sequences and thus difficult to detect in nature. In this study, as a proof-of-concept, we reproducibly demonstrated operon formation in the laboratory. The insertion sequence IS3 and the insertion sequence excision enhancer are genes found in a broad range of bacterial species. We introduced these genes into insertion sequence-less Escherichia coli and found that, supporting our hypothesis, the activity of the two genes altered the expression of genes surrounding IS3, closed a 2.7 kb gap between a pair of genes, and formed new operons. This study shows how insertion sequences can facilitate the rapid formation of operons through locally increasing the structural mutation rates and highlights how coevolution with mobile elements may shape the organization of prokaryotic genomes and gene regulation.

OUP accepted manuscript

Advances in genome sequencing have produced hundreds of thousands of bacterial genome sequences, many of which have integrated prophages derived from temperate bacteriophages. These prophages play key roles by influencing bacterial metabolism, pathogenicity, antibiotic resistance, and defense against viral attack. However, they vary considerably even among related bacterial strains, and they are challenging to identify computationally and to extract precisely for comparative genomic analyses. Here, we describe DEPhT, a multimodal tool for prophage discovery and extraction. It has three run modes that facilitate rapid screening of large numbers of bacterial genomes, precise extraction of prophage sequences, and prophage annotation. DEPhT uses genomic architectural features that discriminate between phage and bacterial sequences for efficient prophage discovery, and targeted homology searches for precise prophage extraction. DEPhT is designed for prophage discovery in Mycobacterium genomes but can be adapted broadly to other bacteria. We deploy DEPhT to demonstrate that prophages are prevalent in Mycobacterium strains but are absent not only from the few well-characterized Mycobacterium tuberculosis strains, but also are absent from all ∼30 000 sequenced M. tuberculosis strains.

A Large-Scale Multiple Genome Comparison of Acidophilic Archaea (pH ≤ 5.0) Extends Our Understanding of Oxidative Stress Responses in Polyextreme Environments

Acidophilic archaea thrive in anaerobic and aerobic low pH environments (pH < 5) rich in dissolved heavy metals that exacerbate stress caused by the production of reactive oxygen species (ROS) such as hydrogen peroxide (H2O2), hydroxyl radical (OH) and superoxide (O2−). ROS react with lipids, proteins and nucleic acids causing oxidative stress and damage that can lead to cell death. Herein, genes and mechanisms potentially involved in ROS mitigation are predicted in over 200 genomes of acidophilic archaea with sequenced genomes. These organisms are often be subjected to simultaneous multiple stresses such as high temperature, high salinity, low pH and high heavy metal loads. Some of the topics addressed include: (1) the phylogenomic distribution of these genes and what this can tell us about the evolution of these mechanisms in acidophilic archaea; (2) key differences in genes and mechanisms used by acidophilic versus non-acidophilic archaea and between acidophilic archaea and acidophilic bacteria and (3) how comparative genomic analysis predicts novel genes or pathways involved in oxidative stress responses in archaea and likely horizontal gene transfer (HGT) events.

Bacterial cooperation through horizontal gene transfer

Cooperation exists across all scales of biological organization, from genetic elements to complex human societies. Bacteria cooperate by secreting molecules that benefit all individuals in the population (i.e., public goods). Genes associated with cooperation can spread among strains through horizontal gene transfer (HGT). We discuss recent findings on how HGT mediated by mobile genetic elements promotes bacterial cooperation, how cooperation in turn can facilitate more frequent HGT, and how the act of HGT itself may be considered as a form of cooperation. We propose that HGT is an important enforcement mechanism in bacterial populations, thus creating a positive feedback loop that further maintains cooperation. To enforce cooperation, HGT serves as a homogenizing force by transferring the cooperative trait, effectively eliminating cheaters.

Comparative genomic analysis of hyper-ammonia producing Acetoanaerobium sticklandii DSM 519 with purinolytic Gottschalkia acidurici 9a and pathogenic Peptoclostridium difficile 630

Acetoanaerobium sticklandii DSM519 (CST) is a hype-ammonia producing non-pathogenic anaerobe that can use amino acids as important carbon and energy sources through the Stickland reactions. Biochemical aspects of this organism have been extensively studied, but systematic studies addressing its metabolic discrepancy remain scant. In this perspective, we have intensively analyzed its genomic and metabolic characteristics to comprehend the evolutionary conservation of amino acid catabolism by a comparative genomic approach. The whole-genome data indicated that CST has shown a phylogenomic similarity with hyper-ammonia producing, purinolytic, and proteolytic pathogenic Clostridia. CST has shown to common genomic context sharing across the purinolytic Gottschalkia acidurici 9a and pathogenic Peptoclostridium difficile 630. Genome syntenic analysis described that syntenic orthologs might be originated from the recent ancestor at a slow evolution rate and syntenic-out paralogs evolved from either CDF or CAC via α-event and β-event. Collinearity of either gene orders or gene families was adjusted with syntenic out-paralogs across these genomes. The genome-wide metabolic analysis predicted 11 unique putative metabolic subsystems from the CST genome for amino acid catabolism and hydrogen production. The in silico analysis of our study revealed that a characteristic system for amino acid catabolism-directed biofuel synthesis might have slowly evolved and established as a core genomic content of CST.

Genic Selection Within Prokaryotic Pangenomes

Understanding the evolutionary forces shaping prokaryotic pangenome structure is a major goal of microbial evolution research. Recent work has highlighted that a substantial proportion of accessory genes appear to confer niche-specific adaptations. This work has primarily focused on selection acting at the level of individual cells. Herein, we discuss a lower level of selection that also contributes to pangenome variation: genic selection. This refers to cases where genetic elements, rather than individual cells, are the entities under selection. The clearest examples of this form of selection are selfish mobile genetic elements, which are those that have either a neutral or a deleterious effect on host fitness. We review the major classes of these and other mobile elements and discuss the characteristic features of such elements that could be under genic selection. We also discuss how genetic elements that are beneficial to hosts can also be under genic selection, a scenario that may be more prevalent but not widely appreciated, because disentangling the effects of selection at different levels (i.e., organisms vs. genes) is challenging. Nonetheless, an appreciation for the potential action and implications of genic selection is important to better understand the evolution of prokaryotic pangenomes.

Unifying themes and distinct features of carbon and nitrogen assimilation by polysaccharide-degrading bacteria: a summary of four model systems

Our current understanding of enzymatic polysaccharide degradation has come from a huge number of in vitro studies with purified enzymes. While this vast body of work has been invaluable in identifying and characterizing novel mechanisms of action and engineering desirable traits into these enzymes, a comprehensive picture of how these enzymes work as part of a native in vivo system is less clear. Recently, several model bacteria have emerged with genetic systems that allow for a more nuanced study of carbohydrate active enzymes (CAZymes) and how their activity affects bacterial carbon metabolism. With these bacterial model systems, it is now possible to not only study a single nutrient system in isolation (i.e., carbohydrate degradation and carbon metabolism), but also how multiple systems are integrated. Given that most environmental polysaccharides are carbon rich but nitrogen poor (e.g., lignocellulose), the interplay between carbon and nitrogen metabolism in polysaccharide-degrading bacteria can now be studied in a physiologically relevant manner. Therefore, in this review, we have summarized what has been experimentally determined for CAZyme regulation, production, and export in relation to nitrogen metabolism for two Gram-positive (Caldicellulosiruptor bescii and Clostridium thermocellum) and two Gram-negative (Bacteroides thetaiotaomicron and Cellvibrio japonicus) polysaccharide-degrading bacteria. By comparing and contrasting these four bacteria, we have highlighted the shared and unique features of each, with a focus on in vivo studies, in regard to carbon and nitrogen assimilation. We conclude with what we believe are two important questions that can act as guideposts for future work to better understand the integration of carbon and nitrogen metabolism in polysaccharide-degrading bacteria. KEY POINTS: • Regardless of CAZyme deployment system, the generation of a local pool of oligosaccharides is a common strategy among Gram-negative and Gram-positive polysaccharide degraders as a means to maximally recoup the energy expenditure of CAZyme production and export. • Due to the nitrogen deficiency of insoluble polysaccharide-containing substrates, Gram-negative and Gram-positive polysaccharide degraders have a diverse set of strategies for supplementation and assimilation. • Future work needs to precisely characterize the energetic expenditures of CAZyme deployment and bolster our understanding of how carbon and nitrogen metabolism are integrated in both Gram-negative and Gram-positive polysaccharide-degrading bacteria, as both of these will significantly influence a given bacterium's suitability for biotechnology applications.

Harmonization of whole-genome sequencing for outbreak surveillance of Enterobacteriaceae and Enterococci

Whole-genome sequencing (WGS) is becoming the de facto standard for bacterial typing and outbreak surveillance of resistant bacterial pathogens. However, interoperability for WGS of bacterial outbreaks is poorly understood. We hypothesized that harmonization of WGS for outbreak surveillance is achievable through the use of identical protocols for both data generation and data analysis. A set of 30 bacterial isolates, comprising of various species belonging to the Enterobacteriaceae family and Enterococcus genera, were selected and sequenced using the same protocol on the Illumina MiSeq platform in each individual centre. All generated sequencing data were analysed by one centre using BioNumerics (6.7.3) for (i) genotyping origin of replications and antimicrobial resistance genes, (ii) core-genome multi-locus sequence typing (cgMLST) for Escherichia coli and Klebsiella pneumoniae and whole-genome multi-locus sequencing typing (wgMLST) for all species. Additionally, a split k-mer analysis was performed to determine the number of SNPs between samples. A precision of 99.0% and an accuracy of 99.2% was achieved for genotyping. Based on cgMLST, a discrepant allele was called only in 2/27 and 3/15 comparisons between two genomes, for E. coli and K. pneumoniae, respectively. Based on wgMLST, the number of discrepant alleles ranged from 0 to 7 (average 1.6). For SNPs, this ranged from 0 to 11 SNPs (average 3.4). Furthermore, we demonstrate that using different de novo assemblers to analyse the same dataset introduces up to 150 SNPs, which surpasses most thresholds for bacterial outbreaks. This shows the importance of harmonization of data-processing surveillance of bacterial outbreaks. In summary, multi-centre WGS for bacterial surveillance is achievable, but only if protocols are harmonized.

The Role of Gene Duplication in the Divergence of Enzyme Function: A Comparative Approach

Gene duplication is a crucial process involved in the appearance of new genes and functions. It is thought to have played a major role in the growth of enzyme families and the expansion of metabolism at the biosphere's dawn and in recent times. Here, we analyzed paralogous enzyme content within each of the seven enzymatic classes for a representative sample of prokaryotes by a comparative approach. We found a high ratio of paralogs for three enzymatic classes: oxidoreductases, isomerases, and translocases, and within each of them, most of the paralogs belong to only a few subclasses. Our results suggest an intricate scenario for the evolution of prokaryotic enzymes, involving different fates for duplicated enzymes fixed in the genome, where around 20-40% of prokaryotic enzymes have paralogs. Intracellular organisms have a lesser ratio of duplicated enzymes, whereas free-living enzymes show the highest ratios. We also found that phylogenetically close phyla and some unrelated but with the same lifestyle share similar genomic and biochemical traits, which ultimately support the idea that gene duplication is associated with environmental adaptation.