Genome and sequence determinants governing the expression of horizontally acquired DNA in bacteria

Engineering the D-lactic acid responsive promoter/repressor system as dynamic metabolic engineering tool in Lactobacillus delbrueckii subsp. bulgaricus for controlled D-lactic acid biosynthesis

Dynamic metabolic engineering integrates synthetic logic circuits into cellular systems, optimizing metabolic fluxes and augmenting biosynthesis of target metabolites. This study evaluated a D-lactic acid (DLA)-responsive promoter-repressor system from Pseudomonas fluorescens A506, re-engineered for heightened sensitivity and functional efficacy in Lactobacillus delbrueckii subsp. bulgaricus VI104. The codon-optimized regulatory architecture exhibited peak performance at DLA inducer concentration range of 60-100 mM, validated by fluorometry and microscopy. As an application, overexpression of D-lactate dehydrogenase (dldh) downstream of the engineered promoter repressor system enabled finely tuned modulation of DLA biosynthesis, autonomously regulating the transition between growth and production phases, thereby attenuating overall metabolic load. Cross-species compatibility was confirmed by excising regulatory elements from the promoter-repressor system and functionally assessing them in recombinant L. bulgaricus. Molecular docking elucidated critical noncovalent interactions between D-LldR repressor and operator nucleotide sequence in absence of inducer DLA. The engineered promoter construct with high efficiency effectively elevated DLA biosynthesis by 2.15-folds and expanded the overall fermentation time relative to constitutive systems, attaining maximum DLA titre of 9.02 g L⁻1 in bioreactor setup. These results substantially broaden the molecular cloning toolkit available for L. bulgaricus, fostering potential future applications in biotherapeutics and probiotics.

Profiling lateral gene transfer events in the human microbiome using WAAFLE

Lateral gene transfer (LGT), also known as horizontal gene transfer, facilitates genomic diversification in microbial populations. While previous work has surveyed LGT in human-associated microbial isolate genomes, the landscape of LGT arising in personal microbiomes is not well understood, as there are no widely adopted methods to characterize LGT from complex communities. Here we developed, benchmarked and validated a computational algorithm (WAAFLE or Workflow to Annotate Assemblies and Find LGT Events) to profile LGT from assembled metagenomes. WAAFLE prioritizes specificity while maintaining high sensitivity for intergenus LGT. Applying WAAFLE to >2,000 human metagenomes from diverse body sites, we identified >100,000 high-confidence previously uncharacterized LGT (~2 per microbial genome-equivalent). These were enriched for mobile elements, as well as restriction-modification functions associated with the destruction of foreign DNA. LGT frequency was influenced by biogeography, phylogenetic similarity of involved pairs (for example, Fusobacterium periodonticum and F. nucleatum) and donor abundance. These forces manifest as networks in which hub taxa donate unequally with phylogenetic neighbours. Our findings suggest that human microbiome LGT may be more ubiquitous than previously described.

The presence of genes encoding carbapenem-hydrolyzing oxacillinase and lack of carbapenem resistance in non-baumannii Acinetobacter misidentified as Acinetobacter baumannii causing bloodstream infections in a tertiary hospital over a 3-year period

OBJECTIVE

This study aims to analyze the genomic and clinical characteristics of Non-baumannii Acinetobacter strains misidentified as A. baumannii, causing bloodstream infections (BSIs) in our hospital.

MATERIALS AND METHODS

Whole genome sequencing was performed and average nucleotide identity (ANI) was analyzed. Susceptibility testing was conducted using micro-broth methods. The distribution of antimicrobial resistance genes (ARGs) and mobile genetic elements (MGEs) was examined using online software tools. The prevalence of virulence factors (VFs) was investigated through nucleotide coding sequence comparisons. Genetic structures of blaOXA genes were analyzed by Gcluster software. Clinical information was collected from electronic medical records for patient characterization.

RESULTS

ANI analysis identified five strains as Acinetobacter pittii, with the remaining four identified as A. geminorum, A. nosocomialis, A. soli and A. bereziniae. The GC content of all isolates was less than 38.9 % except for A. soli 16,294. All Non-baumannii Acinetobacter strains were relatively susceptible to antibiotics, except for one A. pittii isolate. Nine blaOXA variants were identified in seven isolates, with two isolates co-carrying 2 different types of blaOXA. Twenty-four insertion sequences (ISs) were identified, with ISAba and IS17 being the primary ISs. Five A. pittii isolates shared the same genetic structures around blaOXA. Genes related to adherence, immune modulation, and nutritional/metabolic factors were the most frequent. Few VFs were detected in A. soli 16,294 and A.bereziniae 14,325.

CONCLUSIONS

The presence of carbapenem hydrolyzing oxacillinase encoding genes did not confer carbapenem resistance, possibly due to the lack of ISs in the blaOXA flanking sequences. Different blaOXA variants within distinct strains shared the same genetic structures, suggesting potential for multidrug resistance development, which warrants our attention.

Genomic, Phylogenetic and Physiological Characterization of the PAH-Degrading Strain Gordonia polyisoprenivorans 135

The strain Gordonia polyisoprenivorans 135 is able to utilize a wide range of aromatic compounds. The aim of this work was to study the features of genetic organization and biotechnological potential of the strain G. polyisoprenivorans 135 as a degrader of aromatic compounds. The study of the genome of the strain 135 and the pangenome of the G. polyisoprenivorans species revealed that some genes, presumably involved in PAH catabolism, are atypical for Gordonia and belong to the pangenome of Actinobacteria. Analyzing the intergenic regions of strain 135 alongside the "panIGRome" of G. polyisoprenivorans showed that some intergenic regions in strain 135 also differ from those located between the same pairs of genes in related strains. The strain G. polyisoprenivorans 135 in our work utilized naphthalene (degradation degree 39.43%) and grew actively on salicylate. At present, this is the only known strain of G. polyisoprenivorans with experimentally confirmed ability to utilize these compounds.

Exploring the host range for genetic transfer of magnetic organelle biosynthesis

Magnetosomes produced by magnetotactic bacteria have great potential for application in biotechnology and medicine due to their unique physicochemical properties and high biocompatibility. Attempts to transfer the genes for magnetosome biosynthesis into non-magnetic organisms have had mixed results. Here we report on a systematic study to identify key components needed for magnetosome biosynthesis after gene transfer. We transfer magnetosome genes to 25 proteobacterial hosts, generating seven new magnetosome-producing strains. We characterize the recombinant magnetosomes produced by these strains and demonstrate that denitrification and anaerobic photosynthesis are linked to the ability to synthesize magnetosomes upon the gene transfer. In addition, we show that the number of magnetosomes synthesized by a foreign host negatively correlates with the guanine-cytosine content difference between the host and the gene donor. Our findings have profound implications for the generation of magnetized living cells and the potential for transgenic biogenic magnetic nanoparticle production.

r/K selection of GC content in prokaryotes

The guanine/cytosine (GC) content of prokaryotic genomes is species-specific, taking values from 16% to 77%. This diversity of selection for GC content remains contentious. We analyse the correlations between GC content and a range of phenotypic and genotypic data in thousands of prokaryotes. GC content integrates well with these traits into r/K selection theory when phenotypic plasticity is considered. High GC-content prokaryotes are r-strategists with cheaper descendants thanks to a lower average amino acid metabolic cost, colonize unstable environments thanks to flagella and a bacillus form and are generalists in terms of resource opportunism and their defence mechanisms. Low GC content prokaryotes are K-strategists specialized for stable environments that maintain homeostasis via a high-cost outer cell membrane and endospore formation as a response to nutrient deprivation, and attain a higher nutrient-to-biomass yield. The lower proteome cost of high GC content prokaryotes is driven by the association between GC-rich codons and cheaper amino acids in the genetic code, while the correlation between GC content and genome size may be partly due to functional diversity driven by r/K selection. In all, molecular diversity in the GC content of prokaryotes may be a consequence of ecological r/K selection.

Comprehensive essentiality analysis of the Mycobacterium kansasii genome by saturation transposon mutagenesis and deep sequencing

Mycobacterium kansasii (Mk) is an opportunistic pathogen that is frequently isolated from urban water systems, posing a health risk to susceptible individuals. Despite its ability to cause tuberculosis-like pulmonary disease, very few studies have probed the genetics of this opportunistic pathogen. Here, we report a comprehensive essentiality analysis of the Mk genome. Deep sequencing of a high-density library of Mk Himar1 transposon mutants revealed that 86.8% of the chromosomal thymine-adenine (TA) dinucleotide target sites were permissive to insertion, leaving 13.2% TA sites unoccupied. Our analysis identified 394 of the 5,350 annotated open reading frames (ORFs) as essential. The majority of these essential ORFs (84.8%) share essential mutual orthologs with Mycobacterium tuberculosis (Mtb). A comparative genomics analysis identified 139 Mk essential ORFs that share essential orthologs in four other species of mycobacteria. Thirteen Mk essential ORFs share orthologs in all four species that were identified as being not essential, while only two Mk essential ORFs are absent in all species compared. We used the essentiality data and a comparative genomics analysis reported here to highlight differences in essentiality between candidate Mtb drug targets and the corresponding Mk orthologs. Our findings suggest that the Mk genome encodes redundant or additional pathways that may confound validation of potential Mtb drugs and drug target candidates against the opportunistic pathogen. Additionally, we identified 57 intergenic regions containing four or more consecutive unoccupied TA sites. A disproportionally large number of these regions were located upstream of pe/ppe genes. Finally, we present an essentiality and orthology analysis of the Mk pRAW-like plasmid, pMK1248. IMPORTANCE Mk is one of the most common nontuberculous mycobacterial pathogens associated with tuberculosis-like pulmonary disease. Drug resistance emergence is a threat to the control of Mk infections, which already requires long-term, multidrug courses. A comprehensive understanding of Mk biology is critical to facilitate the development of new and more efficacious therapeutics against Mk. We combined transposon-based mutagenesis with analysis of insertion site identification data to uncover genes and other genomic regions required for Mk growth. We also compared the gene essentiality data set of Mk to those available for several other mycobacteria. This analysis highlighted key similarities and differences in the biology of Mk compared to these other species. Altogether, the genome-wide essentiality information generated and the results of the cross-species comparative genomics analysis represent valuable resources to assist the process of identifying and prioritizing potential Mk drug target candidates and to guide future studies on Mk biology.

Profiling novel lateral gene transfer events in the human microbiome

Lateral gene transfer (LGT) is an important mechanism for genome diversification in microbial populations, including the human microbiome. While prior work has surveyed LGT events in human-associated microbial isolate genomes, the scope and dynamics of novel LGT events arising in personal microbiomes are not well understood, as there are no widely adopted computational methods to detect, quantify, and characterize LGT from complex microbial communities. We addressed this by developing, benchmarking, and experimentally validating a computational method (WAAFLE) to profile novel LGT events from assembled metagenomes. Applying WAAFLE to >2K human metagenomes from diverse body sites, we identified >100K putative high-confidence but previously uncharacterized LGT events (~2 per assembled microbial genome-equivalent). These events were enriched for mobile elements (as expected), as well as restriction-modification and transport functions typically associated with the destruction of foreign DNA. LGT frequency was quantifiably influenced by biogeography, the phylogenetic similarity of the involved taxa, and the ecological abundance of the donor taxon. These forces manifest as LGT networks in which hub species abundant in a community type donate unequally with their close phylogenetic neighbors. Our findings suggest that LGT may be a more ubiquitous process in the human microbiome than previously described. The open-source WAAFLE implementation, documentation, and data from this work are available at http://huttenhower.sph.harvard.edu/waafle.

Influence of class 2 integron integrase concentration on gene cassette insertion and excision in vivo

Integron can capture and express antimicrobial resistance gene cassettes and plays important roles in horizontal gene transfer. The establishment of a complete in vitro reaction system will help to reveal integron integrase mediated site-specific recombination process and regulation mechanism. As an enzymatic reaction, the concentration of integrase is assumed to have a great influence on the reaction rate. To determine the influence of different concentrations of integrase on the reaction rate and to find the best range of enzyme concentration were essential to optimizing the in vitro reaction system. In this study, plasmids with gradient transcription levels of class 2 integron integrase gene intI2 under different promoters were constructed. Among plasmids pI2W16, pINTI2N, pI2W, and pI2NW, intI2 transcription levels ranged from about 0.61-fold to 49.65-fold of that in pINTI2N. And the frequencies of gene cassette sat2 integration and excision catalyzed by IntI2 were positively correlated with the transcription levels of intI2 within this range. Western blotting results indicated high expression of IntI2 partly existed in the form of an inclusion body. When compared with Pc of class 1 integron, the spacer sequence of PintI2 can increase the strength of PcW but decrease the strength of PcS. In conclusion, the frequencies of gene cassette integration and excision were positively correlated with the concentration of IntI2. intI2 driving by PcW with PintI2 spacer sequence can obtain the optimum IntI2 concentration required to achieve the maximum recombination efficiency in vivo in this study.

Comprehensive Genomic Characterization of Marine Bacteria Thalassospira spp. Provides Insights into Their Ecological Roles in Aromatic Hydrocarbon-Exposed Environments

The marine bacterial genus Thalassospira has often been identified as an abundant member of polycyclic aromatic hydrocarbon (PAH)-exposed microbial communities. However, despite their potential usability for biotechnological applications, functional genes that are conserved in their genomes have barely been investigated. Thus, the goal of this study was to comprehensively examine the functional genes that were potentially responsible for aromatic hydrocarbon biodegradation in the Thalassospira genomes available from databases, including a new isolate of Thalassospira, strain GO-4, isolated from a phenanthrene-enriched marine bacterial consortium. Strain GO-4 was used in this study as a model organism to link the genomic data and their metabolic functions. Strain GO-4 growth assays confirmed that it utilized a common phenanthrene biodegradation intermediate 2-carboxybenzaldehyde (CBA) as the sole source of carbon and energy, but did not utilize phenanthrene. Consistently, strain GO-4 was found to possess homologous genes of phdK, pht, and pca that encode enzymes for biodegradation of CBA, phthalic acid, and protocatechuic acid, respectively. Further comprehensive genomic analyses for 33 Thalassospira genomes from databases showed that a gene cluster that consisted of phdK and pht homologs was conserved in 13 of the 33 strains. pca gene homologs were found in all examined genomes; however, homologs of the known PAH-degrading genes, such as the pah, phn, or nah genes, were not found. Possibly Thalassospira spp. co-occupy niches with other PAH-degrading bacteria that provide them with PAH degradation intermediates and facilitated their inhabitation in PAH-exposed microbial ecosystems. IMPORTANCE Comprehensive investigation of multiple genomic data sets from targeted microbial taxa deposited in databases may provide substantial information to predict metabolic capabilities and ecological roles in different environments. This study is the first report that details the functional profiling of Thalassospira spp. that have repeatedly been found in polycyclic aromatic hydrocarbon (PAH)-exposed marine bacterial communities by using genomic data from a new isolate, Thalassospira strain GO-4, and other strains in databases. Through screening of functional genes potentially involved in aromatic hydrocarbon biodegradation across 33 Thalassospira genomes and growth assays for strain GO-4, it was suggested that Thalassospira spp. unexceptionally conserved the ability to metabolize single-ring, PAH biodegradation intermediates, while being incapable of utilizing PAHs. This expanded our understanding of this potentially important contributing member to PAH-degrading microbial ecosystems; such species are considered to be specialized in driving downstream reactions of PAH biodegradation.

A positive correlation between GC content and growth temperature in prokaryotes

BACKGROUND

GC pairs are generally more stable than AT pairs; GC-rich genomes were proposed to be more adapted to high temperatures than AT-rich genomes. Previous studies consistently showed positive correlations between growth temperature and the GC contents of structural RNA genes. However, for the whole genome sequences and the silent sites of the codons in protein-coding genes, the relationship between GC content and growth temperature is in a long-lasting debate.

RESULTS

With a dataset much larger than previous studies (681 bacteria and 155 archaea with completely assembled genomes), our phylogenetic comparative analyses showed positive correlations between optimal growth temperature (Topt) and GC content both in bacterial and archaeal structural RNA genes and in bacterial whole genome sequences, chromosomal sequences, plasmid sequences, core genes, and accessory genes. However, in the 155 archaea, we did not observe a significant positive correlation of Topt with whole-genome GC content (GCw) or GC content at four-fold degenerate sites. We randomly drew 155 samples from the 681 bacteria for 1000 rounds. In most cases (> 95%), the positive correlations between Topt and genomic GC contents became statistically nonsignificant (P > 0.05). This result suggested that the small sample sizes might account for the lack of positive correlations between growth temperature and genomic GC content in the 155 archaea and the bacterial samples of previous studies. Comparing the GC content among four categories (psychrophiles/psychrotrophiles, mesophiles, thermophiles, and hyperthermophiles) also revealed a positive correlation between GCw and growth temperature in bacteria. By including the GCw of incompletely assembled genomes, we expanded the sample size of archaea to 303. Positive correlations between GCw and Topt appear especially after excluding the halophilic archaea whose GC contents might be strongly shaped by intense UV radiation.

CONCLUSIONS

This study explains the previous contradictory observations and ends a long debate. Prokaryotes growing in high temperatures have higher GC contents. Thermal adaptation is one possible explanation for the positive association. Meanwhile, we propose that the elevated efficiency of DNA repair in response to heat mutagenesis might have the by-product of increasing GC content like that happens in intracellular symbionts and marine bacterioplankton.

High-Throughput Transcriptional Characterization of Regulatory Sequences from Bacterial Biosynthetic Gene Clusters

Recent efforts to sequence, survey, and functionally characterize the diverse biosynthetic capabilities of bacteria have identified numerous Biosynthetic Gene Clusters (BGCs). Genes found within BGCs are typically transcriptionally silent, suggesting their expression is tightly regulated. To better elucidate the underlying mechanisms and principles that govern BGC regulation on a DNA sequence level, we employed high-throughput DNA synthesis and multiplexed reporter assays to build and to characterize a library of BGC-derived regulatory sequences. Regulatory sequence transcription levels were measured in the Actinobacteria Streptomyces albidoflavus J1074, a popular model strain from a genus rich in BGC diversity. Transcriptional activities varied over 1000-fold in range and were used to identify key features associated with expression, including GC content, transcription start sites, and sequence motifs. Furthermore, we demonstrated that transcription levels could be modulated through coexpression of global regulatory proteins. Lastly, we developed and optimized a S. albidoflavus cell-free expression system for rapid characterization of regulatory sequences. This work helps to elucidate the regulatory landscape of BGCs and provides a diverse library of characterized regulatory sequences for rational engineering and activation of cryptic BGCs.