Maintenance of multipartite genome system and its functional significance in bacteria

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.

Agrobacterium tumefaciens: a Transformative Agent for Fundamental Insights into Host-Microbe Interactions, Genome Biology, Chemical Signaling, and Cell Biology

Agrobacterium tumefaciens incites the formation of readily visible macroscopic structures known as crown galls on plant tissues that it infects. Records from biologists as early as the 17th century noted these unusual plant growths and began examining the basis for their formation. These studies eventually led to isolation of the infectious agent, A. tumefaciens, and decades of study revealed the remarkable mechanisms by which A. tumefaciens causes crown gall through stable horizontal genetic transfer to plants. This fundamental discovery generated a barrage of applications in the genetic manipulation of plants that is still under way. As a consequence of the intense study of A. tumefaciens and its role in plant disease, this pathogen was developed as a model for the study of critical processes that are shared by many bacteria, including host perception during pathogenesis, DNA transfer and toxin secretion, bacterial cell-cell communication, plasmid biology, and more recently, asymmetric cell biology and composite genome coordination and evolution. As such, studies of A. tumefaciens have had an outsized impact on diverse areas within microbiology and plant biology that extend far beyond its remarkable agricultural applications. In this review, we attempt to highlight the colorful history of A. tumefaciens as a study system, as well as current areas that are actively demonstrating its value and utility as a model microorganism.

High-Resolution 3D Genome Map of Brucella Chromosomes in Exponential and Stationary Phases

The three-dimensional (3D) genome structure of an organism or cell is highly relevant to its biological activities, but the availability of 3D genome information for bacteria, especially intracellular pathogens, is still limited. Here, we used Hi-C (high-throughput chromosome conformation capture) technology to determine the 3D chromosome structures of exponential- and stationary-phase Brucella melitensis at a 1-kb resolution. We observed that the contact heat maps of the two B. melitensis chromosomes contain a prominent diagonal and a secondary diagonal. Then, 79 chromatin interaction domains (CIDs) were detected at an optical density at 600 nm (OD600) of 0.4 (exponential phase), with the longest CID being 106 kb and the shortest being 12 kb. Moreover, we obtained 49,363 significant cis-interaction loci and 59,953 significant trans-interaction loci. Meanwhile, 82 CIDs of B. melitensis at an OD600 of 1.5 (stationary phase) were detected, with the longest CID being 94 kb and the shortest being 16 kb. In addition, 25,965 significant cis-interaction loci and 35,938 significant trans-interaction loci were obtained in this phase. Furthermore, we found that as the B. melitensis cells grew from the logarithmic to the plateau phase, the frequency of short-range interactions increased, while that of long-range interactions decreased. Finally, combined analysis of 3D genome and whole-genome transcriptome (RNA-seq) data revealed that the strength of short-range interactions in Chr1 is specifically and strongly correlated with gene expression. Overall, our study provides a global view of the chromatin interactions in the B. melitensis chromosomes, which will serve as a resource for further study of the spatial regulation of gene expression in Brucella. IMPORTANCE The spatial structure of chromatin plays important roles in normal cell functions and in the regulation of gene expression. Three-dimensional genome sequencing has been performed in many mammals and plants, but the availability of such data for bacteria, especially intracellular pathogens, is still limited. Approximately 10% of sequenced bacterial genomes contain more than one replicon. However, how multiple replicons are organized within bacterial cells, how they interact, and whether these interactions help to maintain or segregate these multipartite genomes are unresolved issues. Brucella is a Gram-negative, facultative intracellular, and zoonotic bacterium. Except for Brucella suis biovar 3, Brucella species have two chromosomes. Here, we applied Hi-C technology to determine the 3D genome structures of exponential- and stationary-phase Brucella melitensis chromosomes at a 1-kb resolution. Combined analysis of the 3D genome and RNA-seq data indicated that the strength of short-range interactions in B. melitensis Chr1 is specifically and strongly correlated with gene expression. Our study provides a resource to achieve a deeper understanding of the spatial regulation of gene expression in Brucella.

Bipartite Genomes in Enterobacterales: Independent Origins of Chromids, Elevated Openness and Donors of Horizontally Transferred Genes

Multipartite bacteria have one chromosome and one or more chromid. Chromids are believed to have properties that enhance genomic flexibility, making them a favored integration site for new genes. However, the mechanism by which chromosomes and chromids jointly contribute to this flexibility is not clear. To shed light on this, we analyzed the openness of chromosomes and chromids of the two bacteria, Vibrio and Pseudoalteromonas, both which belong to the Enterobacterales order of Gammaproteobacteria, and compared the genomic openness with that of monopartite genomes in the same order. We applied pangenome analysis, codon usage analysis and the HGTector software to detect horizontally transferred genes. Our findings suggest that the chromids of Vibrio and Pseudoalteromonas originated from two separate plasmid acquisition events. Bipartite genomes were found to be more open compared to monopartite. We found that the shell and cloud pangene categories drive the openness of bipartite genomes in Vibrio and Pseudoalteromonas. Based on this and our two recent studies, we propose a hypothesis that explains how chromids and the chromosome terminus region contribute to the genomic plasticity of bipartite genomes.

Genome-Wide Association Studies across Environmental and Genetic Contexts Reveal Complex Genetic Architecture of Symbiotic Extended Phenotypes

A goal of modern biology is to develop the genotype-phenotype (G→P) map, a predictive understanding of how genomic information generates trait variation that forms the basis of both natural and managed communities. As microbiome research advances, however, it has become clear that many of these traits are symbiotic extended phenotypes, being governed by genetic variation encoded not only by the host's own genome, but also by the genomes of myriad cryptic symbionts. Building a reliable G→P map therefore requires accounting for the multitude of interacting genes and even genomes involved in symbiosis. Here, we use naturally occurring genetic variation in 191 strains of the model microbial symbiont Sinorhizobium meliloti paired with two genotypes of the host Medicago truncatula in four genome-wide association studies (GWAS) to determine the genomic architecture of a key symbiotic extended phenotype-partner quality, or the fitness benefit conferred to a host by a particular symbiont genotype, within and across environmental contexts and host genotypes. We define three novel categories of loci in rhizobium genomes that must be accounted for if we want to build a reliable G→P map of partner quality; namely, (i) loci whose identities depend on the environment, (ii) those that depend on the host genotype with which rhizobia interact, and (iii) universal loci that are likely important in all or most environments. IMPORTANCE Given the rapid rise of research on how microbiomes can be harnessed to improve host health, understanding the contribution of microbial genetic variation to host phenotypic variation is pressing, and will better enable us to predict the evolution of (and select more precisely for) symbiotic extended phenotypes that impact host health. We uncover extensive context-dependency in both the identity and functions of symbiont loci that control host growth, which makes predicting the genes and pathways important for determining symbiotic outcomes under different conditions more challenging. Despite this context-dependency, we also resolve a core set of universal loci that are likely important in all or most environments, and thus, serve as excellent targets both for genetic engineering and future coevolutionary studies of symbiosis.

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.

Centromere Interactions Promote the Maintenance of the Multipartite Genome in Agrobacterium tumefaciens

Many pathogens or symbionts of animals and plants contain multiple replicons, a configuration called a multipartite genome. Multipartite genomes enable those species to replicate their genomes faster and better adapt to new niches. Despite their prevalence, the mechanisms by which multipartite genomes are stably maintained are poorly understood. Agrobacterium tumefaciens is a plant pathogen that contains four replicons: a circular chromosome (Ch1), a linear chromosome (Ch2), and two large plasmids. Recent work indicates that their replication origins are clustered at the cell poles in a manner that depends on their ParB family centromeric proteins: ParB1 for Ch1 and individual RepB paralogs for Ch2 and the plasmids. However, understanding of these interactions and how they contribute to genome maintenance is limited. By combining genome-wide chromosome conformation capture (Hi-C) assays, chromatin-immunoprecipitation sequencing (ChIP-seq), and live cell fluorescence microscopy, we provide evidence here that centromeric clustering is mediated by interactions between these centromeric proteins. We further show that the disruption of centromere clustering results in the loss of replicons. Our data establish the role of centromeric clustering in multipartite genome stability. IMPORTANCE About 10% of sequenced bacteria have multiple replicons, also known as multipartite genomes. How these multipartite genomes are maintained is still poorly understood. Here, we use Agrobacterium tumefaciens as a model and show that the replication origins of the four replicons are clustered through direct interactions between the centromeric proteins; disruption of origin clustering leads to the loss of replicons. Thus, our study provided evidence that centromeric clustering is important for maintaining multipartite genomes.

Conformation and dynamic interactions of the multipartite genome in Agrobacterium tumefaciens

How bacteria with multipartite genomes organize and segregate their DNA is poorly understood. Here, we investigate a prototypical multipartite genome in the plant pathogen Agrobacterium tumefaciens. We identify previously unappreciated interreplicon interactions: the four replicons cluster through interactions at their centromeres, and the two chromosomes, one circular and one linear, interact along their replication arms. Our data suggest that these interreplicon contacts play critical roles in the organization and maintenance of multipartite genomes.

Bacterial species from diverse phyla contain multiple replicons, yet how these multipartite genomes are organized and segregated during the cell cycle remains poorly understood. Agrobacterium tumefaciens has a 2.8-Mb circular chromosome (Ch1), a 2.1-Mb linear chromosome (Ch2), and two large plasmids (pAt and pTi). We used this alpha proteobacterium as a model to investigate the global organization and temporal segregation of a multipartite genome. Using chromosome conformation capture assays, we demonstrate that both the circular and the linear chromosomes, but neither of the plasmids, have their left and right arms juxtaposed from their origins to their termini, generating interarm interactions that require the broadly conserved structural maintenance of chromosomes complex. Moreover, our study revealed two types of interreplicon interactions: “ori-ori clustering” in which the replication origins of all four replicons interact, and “Ch1-Ch2 alignment” in which the arms of Ch1 and Ch2 interact linearly along their lengths. We show that the centromeric proteins (ParB1 for Ch1 and RepB^Ch2 for Ch2) are required for both types of interreplicon contacts. Finally, using fluorescence microscopy, we validated the clustering of the origins and observed their frequent colocalization during segregation. Altogether, our findings provide a high-resolution view of the conformation of a multipartite genome. We hypothesize that intercentromeric contacts promote the organization and maintenance of diverse replicons.

Involvement of Serine / Threonine protein kinases in DNA damage response and cell division in bacteria

The roles of Serine/Threonine protein kinases (STPKs) in bacterial physiology, including bacterial responses to nutritional stresses and under pathogenesis have been well documented. STPKs roles in bacterial cell cycle regulation and DNA damage response have not been much emphasized, possibly because the LexA/RecA type SOS response became the synonym to DNA damage response and cell cycle regulation in bacteria. This review summarizes current knowledge of STPKs genetics, domain organization, and their roles in DNA damage response and cell division regulation in bacteria.

The Pseudoalteromonas multipartite genome: distribution and expression of pangene categories, and a hypothesis for the origin and evolution of the chromid

Bacterial genomes typically consist of one large chromosome, but can also include secondary replicons. These so-called multipartite genomes are scattered on the bacterial tree of life with the majority of cases belonging to Proteobacteria. Within the class gamma-proteobacteria, multipartite genomes are restricted to the two families Vibrionaceae and Pseudoalteromonadaceae. Whereas the genome of vibrios is well studied, information on the Pseudoalteromonadaceae genome is much scarcer. We have studied Pseudoalteromonadaceae with respect to the origin of the chromid, how pangene categories are distributed, how genes are expressed relative to their genomic location, and identified chromid hallmark genes. We calculated the Pseudoalteromonadaceae pangenome based on 25 complete genomes and found that core/softcore are significantly overrepresented in late replicating sectors of the chromid, regardless of how the chromid is replicated. On the chromosome, core/softcore and shell/cloud genes are only weakly overrepresented at the chromosomal replication origin and termination sequences, respectively. Gene expression is trending downwards with increasing distance from the chromosomal oriC, whereas the chromidal expression pattern is more complex. Moreover, we identified 78 chromid hallmark genes, and BLASTp searches suggest that the majority of them were acquired from the ancestral gene pool of Alteromonadales. Finally, our data strongly suggest that the chromid originates from a plasmid that was acquired in a relatively recent event. In summary, this study extends our knowledge on multipartite genomes, and helps us understand how and why secondary replicons are acquired, why they are maintained, and how they are shaped by evolution.

Molecular insights into replication initiation in a multipartite genome harboring bacterium Deinococcus radiodurans

Deinococcus radiodurans harbors a multipartite ploid genome system consisting of two chromosomes and two plasmids present in multiple copies. How these discrete genome elements are maintained and inherited is not well understood. PprA, a pleiotropic protein involved in radioresistance, has been characterized for its roles in DNA repair, genome segregation, and cell division in this bacterium. Here, we show that PprA regulates ploidy of chromosome I and II and inhibits the activity of drDnaA, the initiator protein in D. radiodurans. We found that pprA deletion resulted in an increased genomic content and ploidy of both the chromosomal elements. Expression of PprA in trans rescued the phenotypes of the pprA mutant. To understand the molecular mechanism underlying these phenotypes, we characterized drDnaA and drDnaB. As expected for an initiator protein, recombinant drDnaA showed sequence-specific interactions with the putative oriC sequence in chromosome I (oriCI). Both drDnaA and drDnaB showed ATPase activity, also typical of initiator proteins, but only drDnaB exhibited 5'→3' dsDNA helicase activity in vitro. drDnaA and drDnaB showed homotypic and heterotypic interactions with each other, which were perturbed by PprA. Interestingly, PprA has inhibited the ATPase activity of drDnaA but showed no effect on the activity of drDnaB. Regulation of chromosome copy number and inhibition of the initiator protein functions by PprA strongly suggest that it plays a role as a checkpoint regulator of the DNA replication initiation in D. radiodurans perhaps through its interaction with the replication initiation machinery.

Characterization of ori and parS-like functions in secondary genome replicons in Deinococcus radiodurans

The mechanisms underlying multipartite genome maintenance and its functional significance in extraordinary radioresistance of Deinococcus radiodurans are not well understood. The sequences upstream to parAB operons in chrII (cisII) and MP (cisMP) could stabilize an otherwise, non-replicative colE1 plasmid, in D. radiodurans DnaA and cognate ParB proteins bound specifically with cisII and cisMP elements. The ΔcisII and ΔcisMP cells showed the reduced copy number of cognate replicons and radioresistance as compared with wild type. Fluorescent reporter-operator system inserted in chrI, chrII, and MP in wild type and cisII mutants showed the presence of all three replicons in wild-type cells. Although chrI was present in all the ΔcisII and ΔcisMP cells, nearly half of these cells had chrII and MP, respectively, and the other half had the reduced number of foci representing these replications. These results suggested that cisII and cisMP elements contain both origin of replication and parS-like functions and the secondary genome replicons (chrII and MP) are maintained independent of chrI and have roles in radioresistance of D. radiodurans.

Overcoming the Challenges of Megabase-Sized Plasmid Construction in Escherichia coli

Although Escherichia coli has been a popular tool for plasmid construction, this bacterium was believed to be "unsuitable" for constructing a large plasmid whose size exceeds 500 kilobases. We assumed that traditional plasmid vectors may lack some regulatory DNA elements required for the stable replication and segregation of such a large plasmid. In addition, the use of a few site-specific recombination systems may facilitate cloning of large DNA segments. Here we show two strategies for constructing 1-megabase (1-Mb) secondary chromosomes by using new bacterial artificial chromosome (BAC) vectors. First, the 3-Mb genome of a genome-reduced E. coli strain was split into two chromosomes (2-Mb and 1-Mb), of which the smaller one has the origin of replication and the partitioning locus of the Vibrio tubiashii secondary chromosome. This chromosome fission method (Flp-POP cloning) works via flippase-mediated excision, which coincides with the reassembly of a split chloramphenicol resistance gene, allowing chloramphenicol selection. Next, we developed a new cloning method (oriT-POP cloning) and a fully equipped BAC vector (pMegaBAC1H) for developing a 1-Mb plasmid. Two 0.5-Mb genomic regions were sequentially transferred from two donor strains to a recipient strain via conjugation and captured by pMegaBAC1H in the recipient strain to produce a 1-Mb plasmid. This 1-Mb plasmid was transmissible to another E. coli strain via conjugation. Furthermore, these 1-Mb secondary chromosomes were amplifiable in vitro by using the reconstituted E. coli chromosome replication cycle reaction (RCR). These strategies and technologies would make popular E. coli cells a productive factory for designer chromosome engineering.

Succinylome Analysis Reveals the Involvement of Lysine Succinylation in the Extreme Resistance of Deinococcus radiodurans

Increasing evidence shows that the succinylation of lysine residues mainly regulates enzymes involved in the carbon metabolism pathway, in both prokaryotic and eukaryotic cells. Deinococcus radiodurans is one of the most radioresistant organisms on earth and is famous for its robust resistance. A major goal in the current study of protein succinylation is to explore its function in D. radiodurans. High-resolution LC-MS/MS is used for qualitative proteomics to perform a global succinylation analysis of D. radiodurans and 492 succinylation sites in 270 proteins are identified. These proteins are involved in a variety of biological processes and pathways. It is found that the enzymes involved in nucleic acid binding/processing are enriched in D. radiodurans compared with their previously reported levels in other bacteria. The mutagenesis studies confirm that succinylation regulates the enzymatic activities of species-specific proteins PprI and DdrB, which belong to the radiation-desiccation response regulon. Together, these results provide insight into the role of lysine succinylation in the extreme resistance of D. radiodurans.

ParA proteins of secondary genome elements cross-talk and regulate radioresistance through genome copy number reduction in Deinococcus radiodurans

Deinococcus radiodurans, an extremely radioresistant bacterium has a multipartite genome system and ploidy. Mechanisms underlying such types of bacterial genome maintenance and its role in extraordinary radioresistance are not known in this bacterium. Chromosome I (Chr I), chromosome II (Chr II) and megaplasmid (Mp) encode its own set of genome partitioning proteins. Here, we have characterized P-loop ATPases of Chr II (ParA2) and Mp (ParA3) and their roles in the maintenance of genome copies and extraordinary radioresistance. Purified ParA2 and ParA3 showed nearly similar polymerization kinetics and interaction patterns with DNA. Electron microscopic examination of purified proteins incubated with DNA showed polymerization on nicked circular dsDNA. ParA2 and ParA3 showed both homotypic and heterotypic interactions to each other, but not with ParA1 (ParA of Chr I). Similarly, ParA2 and ParA3 interacted with ParB2 and ParB3 but not with ParB1 in vivo. ParB2 and ParB3 interaction with cis-elements located upstream to the corresponding parAB operon was found to be sequence-specific. Unlike single mutant of parA2 and parA3, their double mutant (ΔparA2ΔParA3) affected copy number of cognate genome elements and resistance to γ-radiation as well as hydrogen peroxide in this bacterium. These results suggested that ParA2 and ParA3 are DNA-binding ATPases producing higher order polymers on DNA and are functionally redundant in the maintenance of secondary genome elements in D. radiodurans. The findings also suggest the involvement of secondary genome elements such as Chr II and Mp in the extraordinary radioresistance of D. radiodurans.