Chromids Aid Genome Expansion and Functional Diversification in the Family Burkholderiaceae

Colony morphotype variation in Burkholderia: implications for success of applications and therapeutics

The Burkholderia genus includes both environmental and pathogenic isolates known for their phenotypic plasticity and adaptability. Burkholderia spp. are intrinsically resistant to many antibiotics, often requiring prolonged therapies during infection. A key feature of Burkholderia spp. is colony morphotype variation (CMV), which allows for rapid adaptation to environmental changes and influences virulence, antibiotic resistance, and pathogenicity by impacting the expression of key virulence factors such as lipopolysaccharides, extracellular DNA, efflux pumps, and flagella. While alternative treatments, such as vaccines and phage therapies, hold promise, CMV has the potential to undermine their efficacy by modifying essential therapeutic targets. Despite its importance, the prevalence and underlying mechanisms of CMV remain poorly understood, leaving critical gaps in our knowledge that may hinder the development of sustainable solutions for managing Burkholderia infections. Addressing these gaps is crucial not only for improving infection management but also for enabling the safe reuse of Burkholderia in biotechnology, where their plant growth-promoting and bioremediation properties are highly valuable. Our goal is to raise awareness within the scientific community about the significance of CMV in Burkholderia, highlighting the urgent need to uncover the mechanisms driving CMV. A deeper understanding of CMV's role in virulence and resistance is essential to developing robust, long-term therapeutic strategies.

Physiology, Heavy Metal Resistance, and Genome Analysis of Two Cupriavidus gilardii Strains Isolated from the Naica Mine (Mexico)

Here, we report the characterization of two Cupriavidus strains, NOV2-1 and OV2-1, isolated from an iron-oxide deposit in an underground tunnel of the Naica mine in Mexico. This unique biotope, characterized by its high temperature (≈50 °C) and the presence of heavy metals, is no longer available for sampling at this time. The genomes of NOV2-1 and OV2-1 comprised two replicons: a chromosome of 3.58 and 3.53 Mb, respectively, and a chromid of 2.1 Mb in both strains. No plasmids were found. The average nucleotide identity and the core genome phylogeny showed that NOV2-1 and OV2-1 belonged to the Cupriavidus gilardii species. NOV2-1 and OV2-1 grew up to 48 °C, with an optimal temperature of 42 °C. Discrete differences were observed between C. gilardii CCUG38401T, NOV2-1, and OV2-1 in the biochemical tests. NOV2-1 and OV2-1 presented resistance to zinc, lead, copper, cadmium, nickel, and cobalt. Several complete and incomplete gene clusters related to the resistance to these heavy metals (ars, czc, cop 1, sil-cop 2, cup, mmf, and mer) were detected in the genome of these strains. Although further studies are needed to determine the origin and role of the detected gene clusters, it is suggested that the czc system may have been mobilized by horizontal gene transfer. This study expands the extreme biotopes where Cupriavidus strains can be retrieved.

Taxonomic distribution of SbmA/BacA and BacA-like antimicrobial peptide transporters suggests independent recruitment and convergent evolution in host-microbe interactions

Antimicrobial peptides (AMPs) are often produced by eukaryotes to control bacterial populations in both pathogenic and mutualistic symbioses. Several pathogens and nitrogen-fixing legume symbionts depend on transporters called SbmA (or BacA) or BclA (BacA-like) to survive exposure to AMPs. However, how broadly these transporters are distributed amongst bacteria, and their evolutionary history, is poorly understood. We used computational approaches, including phylogenetic and sequence similarity analyses, to examine the distribution of SbmA/BacA and BclA proteins across 1,255 species spanning the domain Bacteria, leading to the identification of 71 and 177 SbmA/BacA and BclA proteins, respectively. In vitro sensitivity assays using legume AMPs and several BclA proteins confirmed that AMP transport is a common feature of BclA homologues. Our analyses indicated that SbmA/BacA homologues are encoded only by species in the phylum Pseudomonadota and are primarily found in just two orders: Hyphomicrobiales and Enterobacterales. BclA homologues are somewhat more broadly distributed and were found in clusters across four phyla. These included several orders of the phyla Pseudomonadota and Cyanobacteriota, the order Mycobacteriales (phylum Actinomycetota) and the class Negativicutes (phylum Bacillota). Many of the clades enriched for species encoding SbmA/BacA or BclA homologues are rich in species that interact with eukaryotic hosts in mutualistic or pathogenic interactions. These observations suggest that SbmA/BacA and BclA proteins have been repeatedly co-opted to facilitate associations with eukaryotic hosts by allowing bacteria to cope with host-encoded AMPs.

The Complex Epigenetic Panorama in the Multipartite Genome of the Nitrogen-Fixing Bacterium Sinorhizobium meliloti

In prokaryotes, DNA methylation plays roles in DNA repair, gene expression, cell cycle progression, and immune recognition of foreign DNA. Genome-wide methylation patterns can vary between strains, influencing phenotype, and gene transfer. However, broader evolutionary studies on bacterial epigenomic variation remain limited. In this study, we conducted an epigenomic analysis using single-molecule real-time sequencing on 21 strains of Sinorhizobium meliloti, a facultative plant nitrogen-fixing alphaproteobacterium. This species is notable for its multipartite genome structure, consisting of a chromosome, chromid, and megaplasmid, leading to significant genomic and phenotypic diversity. We identified 16 palindromic and nonpalindromic methylated DNA motifs, including N4-methylcytosine and N6-methyladenine modifications, and analyzed their associated methyltransferases. Some motifs were methylated across all strains, forming a core set of epigenomic signatures, while others exhibited variable methylation frequencies, indicating a dispensable (shell) epigenome. Additionally, we observed differences in methylation frequency between replicons and within coding sequences versus regulatory regions, suggesting that methylation patterns may reflect multipartite genome evolution and influence gene regulation. Overall, our findings reveal extensive epigenomic diversity in S. meliloti, with complex epigenomic signatures varying across replicons and genomic regions. These results enhance our understanding of multipartite genome evolution and highlight the potential role of epigenomic diversity in phenotypic variation.

Transcriptomic Response of Rhizobium leguminosarum to Acidic Stress and Nutrient Limitation Is Versatile and Substantially Influenced by Extrachromosomal Gene Pool

Multipartite genomes are thought to confer evolutionary advantages to bacteria by providing greater metabolic flexibility in fluctuating environments and enabling rapid adaptation to new ecological niches and stress conditions. This genome architecture is commonly found in plant symbionts, including nitrogen-fixing rhizobia, such as Rhizobium leguminosarum bv. trifolii TA1 (RtTA1), whose genome comprises a chromosome and four extrachromosomal replicons (ECRs). In this study, the transcriptomic responses of RtTA1 to partial nutrient limitation and low acidic pH were analyzed using high-throughput RNA sequencing. RtTA1 growth under these conditions resulted in the differential expression of 1035 to 1700 genes (DEGs), which were assigned to functional categories primarily related to amino acid and carbohydrate metabolism, ribosome and cell envelope biogenesis, signal transduction, and transcription. These results highlight the complexity of the bacterial response to stress. Notably, the distribution of DEGs among the replicons indicated that ECRs played a significant role in the stress response. The transcriptomic data align with the Rhizobium pangenome analysis, which revealed an over-representation of functional categories related to transport, metabolism, and regulatory functions on ECRs. These findings confirm that ECRs contribute substantially to the ability of rhizobia to adapt to challenging environmental conditions.

Ancient and remote quartzite caves as a novel source of culturable microbes with biotechnological potential

Quartzite caves located on table-top mountains (tepuis) in the Guyana Shield, are ancient, remote, and pristine subterranean environments where microbes have evolved peculiar metabolic strategies to thrive in silica-rich, slightly acidic and oligotrophic conditions. In this study, we explored the culturable fraction of the microbiota inhabiting the (ortho)quartzite cave systems in Venezuelan tepui (remote table-top mountains) and we investigated their metabolic and enzymatic activities in relation with silica solubilization and extracellular hydrolytic activities as well as the capacity to produce antimicrobial compounds. Eighty microbial strains were isolated with a range of different enzymatic capabilities. More than half of the isolated strains performed at least three enzymatic activities and four bacterial strains displayed antimicrobial activities. The antimicrobial producers Paraburkholderia bryophila CMB_CA002 and Sphingomonas sp. MEM_CA187, were further analyzed by conducting chemotaxonomy, phylogenomics, and phenomics. While the isolate MEM_CA187 represents a novel species of the genus Sphingomonas, for which the name Sphingomonas imawarii sp. nov. is proposed, P. bryophila CMB_CA002 is affiliated with a few strains of the same species that are antimicrobial producers. Chemical analyses demonstrated that CMB_CA002 produces ditropolonyl sulfide that has a broad range of activity and a possibly novel siderophore. Although the antimicrobial compounds produced by MEM_CA187 could not be identified through HPLC-MS analysis due to the absence of reference compounds, it represents the first soil-associated Sphingomonas strain with the capacity to produce antimicrobials. This work provides first insights into the metabolic potential present in quartzite cave systems pointing out that these environments are a novel and still understudied source of microbial strains with biotechnological potential.

Exploring the Role of Inulin in Targeting the Gut Microbiota: An Innovative Strategy for Alleviating Colonic Fibrosis Induced By Irradiation

The use of radiation therapy to treat pelvic and abdominal cancers can lead to the development of either acute or chronic radiation enteropathy. Radiation-induced chronic colonic fibrosis is a common gastrointestinal disorder resulting from the above radiation therapy. In this study, we establish the efficacy of inulin supplements in safeguarding against colonic fibrosis caused by irradiation therapy. Studies have demonstrated that inulin supplements enhance the proliferation of bacteria responsible to produce short-chain fatty acids (SCFAs) and elevate the levels of SCFAs in feces. In a mouse model of chronic radiation enteropathy, the transplantation of gut microbiota and its metabolites from feces of inulin-treated mice were found to reduce colonic fibrosis in validation experiments. Administering inulin-derived metabolites from gut microbiota led to a notable decrease in the expression of genes linked to fibrosis and collagen production in mouse embryonic fibroblast cell line NIH/3T3. In the cell line, inulin-derived metabolites also suppressed the expression of genes linked to the extracellular matrix synthesis pathway. The results indicate a novel and practical approach to safeguarding against chronic radiation-induced colonic fibrosis.

Integrating global microbiome data into antibiotic resistance assessment in large rivers

Rivers are important in spreading antimicrobial resistance (AMR). Assessing AMR risk in large rivers is challenged by large spatial scale and numerous contamination sources. Integrating river resistome data into a global framework may help addressing this difficulty. Here, we conducted an omics-based assessment of AMR in a large river (i.e. the Pearl River in China) with global microbiome data. Results showed that antibiotic resistome in river water and sediment was more diversified than that in other rivers, with contamination levels in some river reaches higher than global baselines. Discharge of WWTP effluent and landfill waste drove AMR prevalence in the river, and the resistome level was highly associated with human and animal sources. Detection of 54 risk rank I ARGs and emerging mobilizable mcr and tet(X) highlighted AMR risk in the river reaches with high human population density and livestock pollution. Florfenicol-resistant floR therein deserved priority concerns due to its high detection frequency, dissimilar phylogenetic distance, mobilizable potential, and presence in multiple pathogens. Co-sharing of ARGs across taxonomic ranks implied their transfer potentials in the community. By comparing with global genomic data, we found that Burkholderiaceae, Enterobacteriaceae, Moraxellaceae and Pseudomonadaceae were important potential ARG-carrying bacteria in the river, and WHO priority carbapenem-resistant Enterobacteriaceae, A. baumannii and P. aeruginosa should be included in future surveillance. Collectively, the findings from this study provide an omics-benchmarked assessment strategy for public risk associated with AMR in large rivers.

The benzoyl-CoA pathway serves as a genomic marker to identify the oxygen requirements in the degradation of aromatic hydrocarbons

The first step of anaerobic benzoate degradation is the formation of benzoyl-coenzyme A by benzoate-coenzyme A ligase (BCL). The anaerobic route is steered by benzoyl-CoA reductase, which promotes benzoyl-CoA breakdown, which is subsequently oxidized. In certain bacteria at low oxygen conditions, the aerobic metabolism of monoaromatic hydrocarbons occurs through the degradation Box pathway. These pathways have undergone experimental scrutiny in Alphaproteobacteria and Betaproteobacteria and have also been explored bioinformatically in representative Betaproteobacteria. However, there is a gap in our knowledge regarding the distribution of the benzoyl-CoA pathway and the evolutionary forces propelling its adaptation beyond that of representative bacteria. To address these questions, we used bioinformatic procedures to identify the BCLs and the lower pathways that transform benzoyl-CoA. These procedures included the identification of conserved motifs. As a result, we identified two motifs exclusive to BCLs, describing some of the catalytic properties of this enzyme. These motifs helped to discern BCLs from other aryl-CoA ligases effectively. The predicted BCLs and the enzymes of lower pathways were used as genomic markers for identifying aerobic, anaerobic, or hybrid catabolism, which we found widely distributed in Betaproteobacteria. Despite these enhancements, our approach failed to distinguish orthologs from a small cluster of paralogs exhibiting all the specified features to predict an ortholog. Nonetheless, the conducted phylogenetic analysis and the properties identified in the genomic context aided in formulating hypotheses about how this redundancy contributes to refining the catabolic strategy employed by these bacteria to degrade the substrates.

Transposon sequencing reveals the essential gene set and genes enabling gut symbiosis in the insect symbiont Caballeronia insecticola

Caballeronia insecticola is a bacterium belonging to the Burkholderia genus sensu lato, which is able to colonize multiple environments like soils and the gut of the bean bug Riptortus pedestris. We constructed a saturated Himar1 mariner transposon library and revealed by transposon-sequencing that 498 protein-coding genes constitute the essential genome of Caballeronia insecticola for growth in free-living conditions. By comparing essential gene sets of Caballeronia insecticola and seven related Burkholderia s.l. strains, only 120 common genes were identified, indicating that a large part of the essential genome is strain-specific. In order to reproduce specific nutritional conditions that are present in the gut of Riptortus pedestris, we grew the mutant library in minimal media supplemented with candidate gut nutrients and identified several condition-dependent fitness-defect genes by transposon-sequencing. To validate the robustness of the approach, insertion mutants in six fitness genes were constructed and their growth deficiency in media supplemented with the corresponding nutrient was confirmed. The mutants were further tested for their efficiency in Riptortus pedestris gut colonization, confirming that gluconeogenic carbon sources, taurine and inositol, are nutrients consumed by the symbiont in the gut. Thus, our study provides insights about specific contributions provided by the insect host to the bacterial symbiont.

Cell cycle-coordinated maintenance of the Vibrio bipartite genome

To preserve the integrity of their genome, bacteria rely on several genome maintenance mechanisms that are co-ordinated with the cell cycle. All members of the Vibrio family have a bipartite genome consisting of a primary chromosome (Chr1) homologous to the single chromosome of other bacteria such as Escherichia coli and a secondary chromosome (Chr2) acquired by a common ancestor as a plasmid. In this review, we present our current understanding of genome maintenance in Vibrio cholerae, which is the best-studied model for bacteria with multi-partite genomes. After a brief overview on the diversity of Vibrio genomic architecture, we describe the specific, common, and co-ordinated mechanisms that control the replication and segregation of the two chromosomes of V. cholerae. Particular attention is given to the unique checkpoint mechanism that synchronizes Chr1 and Chr2 replication.

Identification of the lipodepsipeptide selethramide encoded in a giant nonribosomal peptide synthetase from a Burkholderia bacterium

Bacterial natural products have found many important industrial applications. Yet traditional discovery pipelines often prioritize individual natural product families despite the presence of multiple natural product biosynthetic gene clusters in each bacterial genome. Systematic characterization of talented strains is a means to expand the known natural product space. Here, we report genomics, epigenomics, and metabolomics studies of Burkholderia sp. FERM BP-3421, a soil isolate and known producer of antitumor spliceostatins. Its genome is composed of two chromosomes and two plasmids encoding at least 29 natural product families. Metabolomics studies showed that FERM BP-3421 also produces antifungal aminopyrrolnitrin and approved anticancer romidepsin. From the orphan metabolome features, we connected a lipopeptide of 1,928 Da to an 18-module nonribosomal peptide synthetase encoded as a single gene in chromosome 1. Isolation and structure elucidation led to the identification of selethramide which contains a repeating pattern of serine and leucine and is cyclized at the side chain oxygen of the one threonine residue at position 13. A (R)-3-hydroxybutyric acid moiety decorates the N-terminal serine. Initial attempts to obtain deletion mutants to probe the role of selethramide failed. After acquiring epigenome (methylome) data for FERM BP-3421, we employed a mimicry by methylation strategy that improved DNA transfer efficiency. Mutants defective in selethramide biosynthesis showed reduced surfactant activity and impaired swarming motility that could be chemically complemented with selethramide. This work unveils a lipopeptide that promotes surface motility, establishes improved DNA transfer efficiency, and sets the stage for continued natural product identification from a prolific strain.

Independent origins and evolution of the secondary replicons of the class Gammaproteobacteria

Multipartite genomes, consisting of more than one replicon, have been found in approximately 10 % of bacteria, many of which belong to the phylum Proteobacteria. Many aspects of their origin and evolution, and the possible advantages related to this type of genome structure, remain to be elucidated. Here, we performed a systematic analysis of the presence and distribution of multipartite genomes in the class Gammaproteobacteria, which includes several genera with diverse lifestyles. Within this class, multipartite genomes are mainly found in the order Alteromonadales (mostly in the genus Pseudoalteromonas) and in the family Vibrionaceae. Our data suggest that the emergence of secondary replicons in Gammaproteobacteria is rare and that they derive from plasmids. Despite their multiple origins, we highlighted the presence of evolutionary trends such as the inverse proportionality of the genome to chromosome size ratio, which appears to be a general feature of bacteria with multipartite genomes irrespective of taxonomic group. We also highlighted some functional trends. The core gene set of the secondary replicons is extremely small, probably limited to essential genes or genes that favour their maintenance in the genome, while the other genes are less conserved. This hypothesis agrees with the idea that the primary advantage of secondary replicons could be to facilitate gene acquisition through horizontal gene transfer, resulting in replicons enriched in genes associated with adaptation to different ecological niches. Indeed, secondary replicons are enriched both in genes that could promote adaptation to harsh environments, such as those involved in antibiotic, biocide and metal resistance, and in functional categories related to the exploitation of environmental resources (e.g. carbohydrates), which can complement chromosomal functions.

Selection in coral mitogenomes, with insights into adaptations in the deep sea

Corals are a dominant benthic fauna that occur across a vast range of depths from just below the ocean's surface to the abyssopelagic zone. However, little is known about the evolutionary mechanisms that enable them to inhabit such a wide range of environments. The mitochondrial (mt) genome, which is involved in energetic pathways, may be subject to selection pressures at greater depths to meet the metabolic demands of that environment. Here, we use a phylogenomic framework combined with codon-based models to evaluate whether mt protein-coding genes (PCGs) associated with cellular energy functions are under positive selection across depth in three groups of corals: Octocorallia, Scleractinia, and Antipatharia. The results demonstrated that mt PCGs of deep- and shallow-water species of all three groups were primarily under strong purifying selection (0.0474 < ω < 0.3123), with the exception of positive selection in atp6 (ω = 1.3263) of deep-sea antipatharians. We also found evidence for positive selection at fifteen sites across cox1, mtMutS, and nad1 in deep-sea octocorals and nad3 of deep-sea antipatharians. These results contribute to our limited understanding of mt adaptations as a function of depth and provide insight into the molecular response of corals to the extreme deep-sea environment.

Genomics of the "tumorigenes" clade of the family Rhizobiaceae and description of Rhizobium rhododendri sp. nov

Tumorigenic members of the family Rhizobiaceae, known as agrobacteria, are responsible for crown and cane gall diseases of various crops worldwide. Tumorigenic agrobacteria are commonly found in the genera Agrobacterium, Allorhizobium, and Rhizobium. In this study, we analyzed a distinct "tumorigenes" clade of the genus Rhizobium, which includes the tumorigenic species Rhizobium tumorigenes, as well as strains causing crown gall disease on rhododendron. Here, high-quality, closed genomes of representatives of the "tumorigenes" clade were generated, followed by comparative genomic and phylogenomic analyses. Additionally, the phenotypic characteristics of representatives of the "tumorigenes" clade were analyzed. Our results showed that the tumorigenic strains isolated from rhododendron represent a novel species of the genus Rhizobium for which the name Rhizobium rhododendri sp. nov. is proposed. This species also includes additional strains originating from blueberry and Himalayan blackberry in the United States, whose genome sequences were retrieved from GenBank. Both R. tumorigenes and R. rhododendri contain multipartite genomes, including a chromosome, putative chromids, and megaplasmids. Synteny and phylogenetic analyses indicated that a large putative chromid of R. rhododendri resulted from the cointegration of an ancestral megaplasmid and two putative chromids, following its divergence from R. tumorigenes. Moreover, gene clusters specific for both species of the "tumorigenes" clade were identified, and their biological functions and roles in the ecological diversification of R. rhododendri and R. tumorigenes were predicted and discussed.

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-resolved analyses show an extensive diversification in key aerobic hydrocarbon-degrading enzymes across bacteria and archaea

Background

Hydrocarbons (HCs) are organic compounds composed solely of carbon and hydrogen that are mainly accumulated in oil reservoirs. As the introduction of all classes of hydrocarbons including crude oil and oil products into the environment has increased significantly, oil pollution has become a global ecological problem. However, our perception of pathways for biotic degradation of major HCs and key enzymes in these bioconversion processes has mainly been based on cultured microbes and is biased by uneven taxonomic representation. Here we used Annotree to provide a gene-centric view of the aerobic degradation ability of aliphatic and aromatic HCs in 23,446 genomes from 123 bacterial and 14 archaeal phyla. 

Results

Apart from the widespread genetic potential for HC degradation in Proteobacteria, Actinobacteriota, Bacteroidota, and Firmicutes, genomes from an additional 18 bacterial and 3 archaeal phyla also hosted key HC degrading enzymes. Among these, such degradation potential has not been previously reported for representatives in the phyla UBA8248, Tectomicrobia, SAR324, and Eremiobacterota. Genomes containing whole pathways for complete degradation of HCs were only detected in Proteobacteria and Actinobacteriota. Except for several members of Crenarchaeota, Halobacterota, and Nanoarchaeota that have tmoA, ladA, and alkB/M key genes, respectively, representatives of archaeal genomes made a small contribution to HC degradation. None of the screened archaeal genomes coded for complete HC degradation pathways studied here; however, they contribute significantly to peripheral routes of HC degradation with bacteria.

Conclusion

Phylogeny reconstruction showed that the reservoir of key aerobic hydrocarbon-degrading enzymes in Bacteria and Archaea undergoes extensive diversification via gene duplication and horizontal gene transfer. This diversification could potentially enable microbes to rapidly adapt to novel and manufactured HCs that reach the environment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08906-w.

Minimal overall divergence of the gut microbiome in an adaptive radiation of Cyprinodon pupfishes despite potential adaptive enrichment for scale-eating

Adaptive radiations offer an excellent opportunity to understand the eco-evolutionary dynamics of gut microbiota and host niche specialization. In a laboratory common garden, we compared the gut microbiota of two novel derived trophic specialist pupfishes, a scale-eater and a molluscivore, to closely related and distant outgroup generalist populations, spanning both rapid trophic evolution within 10 kya and stable generalist diets persisting over 11 Mya. We predicted an adaptive and highly divergent microbiome composition in the trophic specialists reflecting their rapid rates of craniofacial and behavioral diversification. We sequenced 16S rRNA amplicons of gut microbiomes from lab-reared adult pupfishes raised under identical conditions and fed the same high protein diet. In contrast to our predictions, gut microbiota largely reflected phylogenetic distance among species, rather than generalist or specialist life history, in support of phylosymbiosis. However, we did find significant enrichment of Burkholderiaceae bacteria in replicated lab-reared scale-eater populations. These bacteria sometimes digest collagen, the major component of fish scales, supporting an adaptive shift. We also found some enrichment of Rhodobacteraceae and Planctomycetia in lab-reared molluscivore populations, but these bacteria target cellulose. Overall phylogenetic conservation of microbiome composition contrasts with predictions of adaptive radiation theory and observations of rapid diversification in all other trophic traits in these hosts, including craniofacial morphology, foraging behavior, aggression, and gene expression, suggesting that the functional role of these minor shifts in microbiota will be important for understanding the role of the microbiome in trophic diversification.

Discovery of Chlorophyll d: Isolation and Characterization of a Far-Red Cyanobacterium from the Original Site of Manning and Strain (1943) at Moss Beach, California

We have isolated a chlorophyll-d-containing cyanobacterium from the intertidal field site at Moss Beach, on the coast of Central California, USA, where Manning and Strain (1943) originally discovered this far-red chlorophyll. Here, we present the cyanobacterium’s environmental description, culturing procedure, pigment composition, ultrastructure, and full genome sequence. Among cultures of far-red cyanobacteria obtained from red algae from the same site, this strain was an epiphyte on a brown macroalgae. Its Q_yin vivo absorbance peak is centered at 704–705 nm, the shortest wavelength observed thus far among the various known Acaryochloris strains. Its Chl a/Chl d ratio was 0.01, with Chl d accounting for 99% of the total Chl d and Chl a mass. TEM imagery indicates the absence of phycobilisomes, corroborated by both pigment spectra and genome analysis. The Moss Beach strain codes for only a single set of genes for producing allophycocyanin. Genomic sequencing yielded a 7.25 Mbp circular chromosome and 10 circular plasmids ranging from 16 kbp to 394 kbp. We have determined that this strain shares high similarity with strain S15, an epiphyte of red algae, while its distinct gene complement and ecological niche suggest that this strain could be the closest known relative to the original Chl d source of Manning and Strain (1943). The Moss Beach strain is designated Acaryochloris sp. (marina) strain Moss Beach.

Natural Chromosome-Chromid Fusion across rRNA Operons in a Burkholderiaceae Bacterium

Chromids (secondary chromosomes) in bacterial genomes that are present in addition to the main chromosome appear to be evolutionarily conserved in some specific bacterial groups. In rare cases among these groups, a small number of strains from Rhizobiales and Vibrionales were shown to possess a naturally fused single chromosome that was reported to have been generated through intragenomic homologous recombination between repeated sequences on the chromosome and chromid. Similar examples have never been reported in the family Burkholderiaceae, a well-documented group that conserves chromids. Here, an in-depth genomic characterization was performed on a Burkholderiaceae bacterium that was isolated from a soil bacterial consortium maintained on diesel fuel and mutagenic benzo[a]pyrene. This organism, Cupriavidus necator strain KK10, was revealed to carry a single chromosome with unexpectedly large size (>6.6 Mb), and results of comparative genomics with the genome of C. necator N-1T indicated that the single chromosome of KK10 was generated through fusion of the prototypical chromosome and chromid at the rRNA operons. This fusion hypothetically occurred through homologous recombination with a crossover between repeated rRNA operons on the chromosome and chromid. Some metabolic functions that were likely expressed from genes on the prototypical chromid region were indicated to be retained. If this phenomenon-the bacterial chromosome-chromid fusion across the rRNA operons through homologous recombination-occurs universally in prokaryotes, the multiple rRNA operons in bacterial genomes may not only contribute to the robustness of ribosome function, but also provide more opportunities for genomic rearrangements through frequent recombination. IMPORTANCE A bacterial chromosome that was naturally fused with the secondary chromosome, or "chromid," and presented as an unexpectedly large single replicon was discovered in the genome of Cupriavidus necator strain KK10, a biotechnologically useful member of the family Burkholderiaceae. Although Burkholderiaceae is a well-documented group that conserves chromids in their genomes, this chromosomal fusion event has not been previously reported for this family. This fusion has hypothetically occurred through intragenomic homologous recombination between repeated rRNA operons and, if so, provides novel insight into the potential of multiple rRNA operons in bacterial genomes to lead to chromosome-chromid fusion. The harsh conditions under which strain KK10 was maintained-a genotoxic hydrocarbon-enriched milieu-may have provided this genotype with a niche in which to survive.

Mining of Cyanobacterial Genomes Indicates Natural Product Biosynthetic Gene Clusters Located in Conjugative Plasmids

Microbial natural products are compounds with unique chemical structures and diverse biological activities. Cyanobacteria commonly possess a wide range of biosynthetic gene clusters (BGCs) to produce natural products. Although natural product BGCs have been found in almost all cyanobacterial genomes, little attention has been given in cyanobacterial research to the partitioning of these biosynthetic pathways in chromosomes and plasmids. Cyanobacterial plasmids are believed to disperse several natural product BGCs, such as toxins, by plasmids through horizontal gene transfer. Therefore, plasmids may confer the ability to produce toxins and may play a role in the evolution of diverse natural product BGCs from cyanobacteria. Here, we performed an analysis of the distribution of natural product BGCs in 185 genomes and mapped the presence of genes involved in the conjugation in plasmids. The 185 analyzed genomes revealed 1817 natural products BGCs. Individual genomes contained 1–42 biosynthetic pathways (mean 8), 95% of which were present in chromosomes and the remaining 5% in plasmids. Of the 424 analyzed cyanobacterial plasmids, 12% contained homologs of genes involved in conjugation and natural product biosynthetic pathways. Among the biosynthetic pathways in plasmids, manual curation identified those to produce aeruginosin, anabaenopeptin, ambiguine, cryptophycin, hassallidin, geosmin, and microcystin. These compounds are known toxins, protease inhibitors, odorous compounds, antimicrobials, and antitumorals. The present study provides in silico evidence using genome mining that plasmids may be involved in the distribution of natural product BGCs in cyanobacteria. Consequently, cyanobacterial plasmids have importance in the context of biotechnology, water management, and public health risk assessment. Future research should explore in vivo conjugation and the end products of natural product BGCs in plasmids via chemical analyses.

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.

Changes in Growth Performance and Ileal Microbiota Composition by Xylanase Supplementation in Broilers Fed Wheat-Based Diets

Xylanase exerts key roles in improving growth performance and intestinal health of broilers fed wheat-based diets. However, knowledge is limited regarding effects of xylanase supplementation on ileal microbiota in broilers. A total of 128 one-day-old broilers (initial BW 48.03 ± 0.33 g) were selected to investigate effects of xylanase (AT-xynA) on growth performance, ileal morphology, microbiota composition, immune response, antioxidant capacity, and endocrine peptide levels in broilers. Broilers were randomly allotted into two dietary treatments (n = 8), namely, a wheat-soybean basal diet and a basal diet with 4,000 U/kg AT-xynA (XY). On days 7, 14, 21, and 42, broilers were weighted and ileal tissues were sampled. Ileal digesta samples were collected for analyzing microbiota composition on days 21 and 42. The results showed that AT-xynA could improve average daily weight gain and average daily feed intake, and there were interactions between diet and age of broilers (p < 0.05). On days 21 and 42, xylanase supplementation decreased ileal microbiota α-diversity, and the relative abundance of potentially pathogenic microbiota, such as phylum Proteobacteria, family Moraxellaceae and Staphylococcaceae, genus Staphylococcus, Pseudomonas, Streptococcus, and Enterococcus, increased the abundance of Lactobacillus (p < 0.05). Moreover, the reduction in acetate concentration and abundance of short-chain fatty acid-producing bacteria was also observed in broilers from XY group (p < 0.05). AT-xynA increased ileal villus height, glucagon-like peptide-1, and insulin-like growth factor-1 concentrations and decreased interleukin-1β, interleukin-6, tumor necrosis factor-α, and malondialdehyde content in broilers, and these positive effects on intestinal health were greater in young broilers. In conclusion, xylanase supplementation to wheat-based diets could improve ileal intestinal morphology and immune function, and alleviate excess fermentation of bacteria, which may be related to changes of intestinal microbiota. In addition, the positive effects of xylanase on intestinal health were more pronounced in young broilers, thus contributing to subsequent improvement in growth performance of broilers.

Presence of the Hmq System and Production of 4-Hydroxy-3-Methyl-2-Alkylquinolines Are Heterogeneously Distributed between Burkholderia cepacia Complex Species and More Prevalent among Environmental than Clinical Isolates

The Burkholderia cepacia complex (Bcc) comprises several species of closely related, versatile bacteria. Some Bcc strains produce 4-hydroxy-3-methyl-2-alkylquinolines (HMAQs), analogous to the 4-hydroxy-2-alkylquinolines of Pseudomonas aeruginosa. Using in silico analyses, we previously estimated that the hmqABCDEFG operon, which encodes enzymes involved in the biosynthesis of HMAQs, is carried by about one-third of Bcc strains, with considerable inter- and intraspecies variability. In the present study, we investigated by PCR, using consensus primers, the distribution of hmqABCDEFG in a collection of 312 Bcc strains (222 of clinical and 90 of environmental origins) belonging to 18 Bcc species. We confirmed that this operon is not distributed evenly among Bcc species. Among the 30% of strains bearing the hmqABCDEFG operon, we found that 92% of environmental isolates and 82% of clinically isolated Bcc strains produce levels of HMAQs detectable by liquid chromatography-mass spectrometry in at least one of the tested culture conditions. Among the hmqABCDEFG-positive but HMAQ-negative strains, none expressed the hmqA gene under the specified culture conditions. Interestingly, the hmqABCDEFG operon is more prevalent among plant root environment species (e.g., Burkholderia ambifaria and Burkholderia cepacia) and absent in species commonly found in chronically colonized individuals with cystic fibrosis (e.g., Burkholderia cenocepacia and Burkholderia multivorans), suggesting a role for the Hmq system in niche adaptation. We investigated the impact of the Hmq system on plant growth promotion and found that Pisum sativum root development by B. ambifaria required a functional HMAQ system.

IMPORTANCE Environmental bacteria belonging to the various closely related species forming the Burkholderia cepacia complex (Bcc) can infect plants and animals, including humans. Their pathogenicity is regulated by intercellular communication, or quorum sensing, allowing them to collaborate instead of acting individually. Bcc organisms generally exploit interacting quorum sensing systems based on N-acyl-homoserine lactones as signaling molecules. Several Bcc strains also carry an hmqABCDEFG operon responsible for the biosynthesis of 4-hydroxy-3-methyl-2-alkylquinolines (HMAQs), molecules analogous to the Pseudomonas quinolone signal (PQS) system of P. aeruginosa. Our finding that the prevalences of the Hmq system and HMAQ production are very different between various Bcc species suggests a key role in niche adaptation or pathogenicity. This is supported by a significant reduction in plant growth promotion in the absence of HMAQ production for a beneficial Bcc strain.

The Location of Substitutions and Bacterial Genome Arrangements

Increasing evidence supports the notion that different regions of a genome have unique rates of molecular change. This variation is particularly evident in bacterial genomes where previous studies have reported gene expression and essentiality tend to decrease, whereas substitution rates usually increase with increasing distance from the origin of replication. Genomic reorganization such as rearrangements occur frequently in bacteria and allow for the introduction and restructuring of genetic content, creating gradients of molecular traits along genomes. Here, we explore the interplay of these phenomena by mapping substitutions to the genomes of Escherichia coli, Bacillus subtilis, Streptomyces, and Sinorhizobium meliloti, quantifying how many substitutions have occurred at each position in the genome. Preceding work indicates that substitution rate significantly increases with distance from the origin. Using a larger sample size and accounting for genome rearrangements through ancestral reconstruction, our analysis demonstrates that the correlation between the number of substitutions and the distance from the origin of replication is significant but small and inconsistent in direction. Some replicons had a significantly decreasing trend (E. coli and the chromosome of S. meliloti), whereas others showed the opposite significant trend (B. subtilis, Streptomyces, pSymA and pSymB in S. meliloti). dN, dS, and ω were examined across all genes and there was no significant correlation between those values and distance from the origin. This study highlights the impact that genomic rearrangements and location have on molecular trends in some bacteria, illustrating the importance of considering spatial trends in molecular evolutionary analysis. Assuming that molecular trends are exclusively in one direction can be problematic.

Plasmids Related to the Symbiotic Nitrogen Fixation Are Not Only Cooperated Functionally but Also May Have Evolved over a Time Span in Family Rhizobiaceae

Rhizobia are soil bacteria capable of forming symbiotic nitrogen-fixing nodules associated with leguminous plants. In fast-growing legume-nodulating rhizobia, such as the species in the family Rhizobiaceae, the symbiotic plasmid is the main genetic basis for nitrogen-fixing symbiosis, and is susceptible to horizontal gene transfer. To further understand the symbioses evolution in Rhizobiaceae, we analyzed the pan-genome of this family based on 92 genomes of type/reference strains and reconstructed its phylogeny using a phylogenomics approach. Intriguingly, although the genetic expansion that occurred in chromosomal regions was the main reason for the high proportion of low-frequency flexible gene families in the pan-genome, gene gain events associated with accessory plasmids introduced more genes into the genomes of nitrogen-fixing species. For symbiotic plasmids, although horizontal gene transfer frequently occurred, transfer may be impeded by, such as, the host’s physical isolation and soil conditions, even among phylogenetically close species. During coevolution with leguminous hosts, the plasmid system, including accessory and symbiotic plasmids, may have evolved over a time span, and provided rhizobial species with the ability to adapt to various environmental conditions and helped them achieve nitrogen fixation. These findings provide new insights into the phylogeny of Rhizobiaceae and advance our understanding of the evolution of symbiotic nitrogen fixation.

Rhizobium indicum sp. nov., isolated from root nodules of pea (Pisum sativum) cultivated in the Indian trans-Himalayas

Three strains of rhizobia isolated from effective root nodules of pea (Pisum sativum L.) collected from the Indian trans-Himalayas were characterized using 16S rRNA, atpD and recA genes. Phylogeny of the 16S rRNA genes revealed that the newly isolated strains were members of the genus Rhizobium with ≥99.9% sequence similarity to the members within the "Rhizobium leguminosarum" group. Phylogenetic analyses based on the concatenated sequences of atpD and recA gene, and 92 core genes extracted from the genome sequences indicated that strains JKLM 12A2T and JKLM 13E are grouped as a separate clade closely related to R. laguerreae FB206T. In contrast, the strain JKLM 19E was placed with "R. hidalgonense" FH14T. Whole-genome average nucleotide identity (ANI) values were 97.6% within strains JKLM 12A2T and JKLM 13E, and less than 94% with closely related species. The digital DNA-DNA hybridization (dDDH) values were 81.45 within the two strains and less than 54.8% to closely related species. The major cellular fatty acids were C18:1w7c in summed feature 8, C14:0 3OH/C16:1 iso I in summed feature 2, and C18:0. The DNA G+C content of JKLM 12A2T and JKLM 13E was 60.8mol%. The data on genomic, chemotaxonomic, and phenotypic characteristics indicates that the strains JKLM 12A2T and JKLM 13E represent a novel species, Rhizobium indicum sp. nov. The type strain is JKLM 12A2T (= MCC 3961T=KACC 21380T=JCM 33658T). However, the strain JKLM 19E represents a member of "R. hidalgonense" and the symbiovar viciae.

Mycoavidus sp. Strain B2-EB: Comparative Genomics Reveals Minimal Genomic Features Required by a Cultivable Burkholderiaceae-Related Endofungal Bacterium

Obligate bacterial endosymbionts are critical to the existence of many eukaryotes. Such endobacteria are usually characterized by reduced genomes and metabolic dependence on the host, which may cause difficulty in isolating them in pure cultures. Family Burkholderiaceae-related endofungal bacteria affiliated with the Mycoavidus-Glomeribacter clade can be associated with the fungal subphyla Mortierellomycotina and Glomeromycotina. In this study, a cultivable endosymbiotic bacterium, Mycoavidus sp. strain B2-EB, present in the fungal host Mortierella parvispora was obtained successfully. The B2-EB genome (1.88 Mb) represents the smallest genome among the endofungal bacterium Mycoavidus cysteinexigens (2.64-2.80 Mb) of Mortierella elongata and the uncultured endosymbiont "Candidatus Glomeribacter gigasporarum" (1.37 to 2.36 Mb) of arbuscular mycorrhizal fungi. Despite a reduction in genome size, strain B2-EB displays a high genome completeness, suggesting a nondegenerative reduction in the B2-EB genome. Compared with a large proportion of transposable elements (TEs) in other known Mycoavidus genomes (7.2 to 11.5% of the total genome length), TEs accounted for only 2.4% of the B2-EB genome. This pattern, together with a high proportion of single-copy genes in the B2-EB genome, suggests that the B2-EB genome reached a state of relative evolutionary stability. These results represent the most streamlined structure among the cultivable endofungal bacteria and suggest the minimal genome features required by both an endofungal lifestyle and artificial culture. This study allows us to understand the genome evolution of Burkholderiaceae-related endosymbionts and to elucidate microbiological interactions.IMPORTANCE This study attempted the isolation of a novel endobacterium, Mycoavidus sp. B2-EB (JCM 33615), harbored in the fungal host Mortierella parvispora E1425 (JCM 39028). We report the complete genome sequence of this strain, which possesses a reduced genome size with relatively high genome completeness and a streamlined genome structure. The information indicates the minimal genomic features required by both the endofungal lifestyle and artificial cultivation, which furthers our understanding of genome reduction in fungal endosymbionts and extends the culture resources for biotechnological development on engineering synthetic microbiomes.

Analysis of the biodegradative and adaptive potential of the novel polychlorinated biphenyl degrader Rhodococcus sp. WAY2 revealed by its complete genome sequence

The complete genome sequence of Rhodococcus sp. WAY2 (WAY2) consists of a circular chromosome, three linear replicons and a small circular plasmid. The linear replicons contain typical actinobacterial invertron-type telomeres with the central CGTXCGC motif. Comparative phylogenetic analysis of the 16S rRNA gene along with phylogenomic analysis based on the genome-to-genome blast distance phylogeny (GBDP) algorithm and digital DNA–DNA hybridization (dDDH) with other Rhodococcus type strains resulted in a clear differentiation of WAY2, which is likely a new species. The genome of WAY2 contains five distinct clusters of bph, etb and nah genes, putatively involved in the degradation of several aromatic compounds. These clusters are distributed throughout the linear plasmids. The high sequence homology of the ring-hydroxylating subunits of these systems with other known enzymes has allowed us to model the range of aromatic substrates they could degrade. Further functional characterization revealed that WAY2 was able to grow with biphenyl, naphthalene and xylene as sole carbon and energy sources, and could oxidize multiple aromatic compounds, including ethylbenzene, phenanthrene, dibenzofuran and toluene. In addition, WAY2 was able to co-metabolize 23 polychlorinated biphenyl congeners, consistent with the five different ring-hydroxylating systems encoded by its genome. WAY2 could also use n-alkanes of various chain-lengths as a sole carbon source, probably due to the presence of alkB and ladA gene copies, which are only found in its chromosome. These results show that WAY2 has a potential to be used for the biodegradation of multiple organic compounds.