Search

The chromosomal position effect can significantly affect the transgene expression, which may provide an efficient strategy for the inauguration of alien genes in new hosts, but has been less explored rationally. The bacterium Myxococcus xanthus harbors a large circular high-GC genome, and the position effect in this chassis may result in a thousand-fold expression variation of alien natural products. In this study, we conducted transposon insertion at TA sites on the M. xanthus genome, and used enrichment and dilution indexes to respectively appraise high and low expression potentials of alien genes at insertion sites. The enrichment sites are characteristically distributed along the genome, and the dilution sites are overlapped well with the horizontal transfer genes. We experimentally demonstrated the enrichment sites as high expression integration sites (HEISs), and the dilution sites unsuitable for gene integration expression. This work highlights that HEISs are the plug-and-play sites for efficient expression of integrated genes.

Abstract Background Functional redundancy (FR) is widely present, but there is no consensus on its formation process and influencing factors. Taxonomically distinct microorganisms possessing genes for the same function in a community lead to within-community FR, and distinct assemblies of microorganisms in different communities playing the same functional roles are termed between-community FR. We proposed two formulas to respectively quantify the degree of functional redundancy within and between communities and analyzed the FR degrees of carbohydrate degradation functions in global environment samples using the genetic information of glycoside hydrolases (GHs) encoded by prokaryotes. Results Our results revealed that GHs are each encoded by multiple taxonomically distinct prokaryotes within a community, and the enzyme-encoding prokaryotes are further distinct between almost any community pairs. The within- and between-FR degrees are primarily affected by the alpha and beta community diversities, respectively, and are also affected by environmental factors (e.g., pH, temperature, and salinity). The FR degree of the prokaryotic community is determined by deterministic factors. Conclusions We conclude that the functional redundancy of GHs is a stabilized community characteristic. This study helps to determine the FR formation process and influencing factors and provides new insights into the relationships between prokaryotic community biodiversity and ecosystem functions. Video Abstract

ABSTRACT Gene duplication events happen prevalently during evolution, and the mechanisms governing the loss or retention of duplicated genes are mostly elusive. Our genome scanning analysis revealed that trigger factor (TF), the one and only bacterial ribosome-associated molecular chaperone, is singly copied in virtually every bacterium except for a very few that possess two or more copies. However, even in these exceptions, only one complete TF copy exists, while other homologs lack the N-terminal domain that contains the conserved ribosome binding site (RBS) motif. Consistently, we demonstrated that the overproduction of the N-terminal complete TF proteins is detrimental to the cell, which can be rescued by removing the N-terminal domain. Our findings also indicated that TF overproduction leads to a decrease in protein productivity and profile changes in proteome due to its characteristic ribosome binding and holdase activities. Additionally, these N-terminal deficient TF homologs in bacteria with multiple TF homologs partition the function of TF via subfunctionalization. Our results revealed that TF is subjected to a dosage constraint that originates from its own intrinsic functions, which may drive the evolution and fates of duplicated TFs in bacteria.IMPORTANCEGene duplication events presumably occur in tig, which encodes the ribosome-associated molecular chaperone trigger factor (TF). However, TF is singly copied in virtually every bacterium, and these exceptions with multiple TF homologs always retain only one complete copy while other homologs lack the N-terminal domain. Here, we reveal the manner and mechanism underlying the evolution and fates of TF duplicates in bacteria. We discovered that the mutation-to-loss or retention-to-sub/neofunctionalization of TF duplicates is associated with the dosage constraint of N-terminal complete TF. The dosage constraint of TF is attributed to its characteristic ribosome binding and substrate-holding activities, causing a decrease in protein productivity and profile changes in cellular proteome.

ABSTRACT ClpXP is a protease complex that plays important roles in protein quality control and cell cycle regulation, but the functions of multiple ClpXs and multiple ClpPs in M. xanthus remain unknown. The genome of Myxococcus xanthus DK1622 contains two clpPs and three clpXs. The clpP1 and clpX1 genes are cotranscribed and are both essential, while the other copies are isolated in the genome and are deletable. The deletion of clpX2 caused the mutant to be deficient in fruiting body development, while the clpX3 gene is involved in resistance to thermal stress. Both ClpPs possess catalytic active sites, but only ClpP1 shows in vitro peptidase activity on the typical substrate Suc-LY-AMC. All of these clpP and clpX genes exhibit strong transcriptional upregulation in the stationary phase, and the transcription of the three clpX genes appears to be coordinated. Our results demonstrated that multiple ClpPs and multiple ClpXs are functionally divergent and may assist in the environmental adaptation and functional diversification of M. xanthus.IMPORTANCEClpXP is an important protease complex of bacteria and is involved in various physiological processes. Myxococcus xanthus DK1622 possesses two ClpPs and three ClpXs with unclear functions. We investigated the functions of these genes and demonstrated the essential roles of clpP1 and clpX1. Only ClpP1 has in vitro peptidase activity on Suc-LY-AMC, and the isolated clpX copies participate in distinct cellular processes. All of these genes exhibited significant transcriptional upregulation in the stationary phase. Divergent functions appear in multiple ClpPs and multiple ClpXs in M. xanthus DK1622.

Myxobacteria are well known for multicellular social behaviors and valued for biosynthesis of natural products. Myxobacteria social behaviors such as clumping growth severely hamper strain cultivation and genetic manipulation. Using Myxococcus xanthus DK1622, we engineered Hu04, which is deficient in multicellular behavior and pigmentation. Hu04, while maintaining nutritional growth and a similar metabolic background, exhibits improved dispersed growth, streamlining operational procedures. It achieves high cell densities in culture and is promising for synthetic biology applications.

ABSTRACT The clustered regularly interspaced short palindromic repeats and their associated proteins (CRISPR-Cas) system widely occurs in prokaryotic organisms to recognize and destruct genetic invaders. Systematic collation and characterization of endogenous CRISPR-Cas systems are conducive to our understanding and potential utilization of this natural genetic machinery. In this study, we screened 39 complete and 692 incomplete genomes of myxobacteria using a combined strategy to dispose of the abridged genome information and revealed at least 19 CRISPR-Cas subtypes, which were distributed with a taxonomic difference and often lost stochastically in intraspecies strains. The cas genes in each subtype were evolutionarily clustered but deeply separated, while most of the CRISPRs were divided into four types based on the motif characteristics of repeat sequences. The spacers recorded in myxobacterial CRISPRs were in high G+C content, matching lots of phages, tiny amounts of plasmids, and, surprisingly, massive organismic genomes. We experimentally demonstrated the immune and self-target immune activities of three endogenous systems in Myxococcus xanthus DK1622 against artificial genetic invaders and revealed the microhomology-mediated end-joining mechanism for the immunity-induced DNA repair but not homology-directed repair. The panoramic view and immune activities imply potential omnipotent immune functions and applications of the endogenous CRISPR-Cas machinery.IMPORTANCEServing as an adaptive immune system, clustered regularly interspaced short palindromic repeats and their associated proteins (CRISPR-Cas) empower prokaryotes to fend off the intrusion of external genetic materials. Myxobacteria are a collective of swarming Gram-stain-negative predatory bacteria distinguished by intricate multicellular social behavior. An in-depth analysis of their intrinsic CRISPR-Cas systems is beneficial for our understanding of the survival strategies employed by host cells within their environmental niches. Moreover, the experimental findings presented in this study not only suggest the robust immune functions of CRISPR-Cas in myxobacteria but also their potential applications.

ABSTRACTThe chaperone 70 kDa heat shock protein (Hsp70) is important for cells from bacteria to humans to maintain proteostasis, and all eukaryotes and several prokaryotes encode Hsp70 paralogs. Although the mechanisms of Hsp70 function have been clearly illuminated, the function and evolution of Hsp70 paralogs is not well studied. DnaK is a highly conserved bacterial Hsp70 family. Here, we show that dnaK is present in 98.9% of bacterial genomes, and 6.4% of them possess two or more DnaK paralogs. We found that the duplication of dnaK is positively correlated with an increase in proteomic complexity (proteome size, number of domains). We identified the interactomes of the two DnaK paralogs of Myxococcus xanthus DK1622 (MxDnaKs), which revealed that they are mostly nonoverlapping, although both prefer α and β domain proteins. Consistent with the entire M. xanthus proteome, MxDnaK substrates have both significantly more multi-domain proteins and a higher isoelectric point than that of Escherichia coli, which encodes a single DnaK homolog. MxDnaK1 is transcriptionally upregulated in response to heat shock and prefers to bind cytosolic proteins, while MxDnaK2 is downregulated by heat shock and is more associated with membrane proteins. Using domain swapping, we show that the nucleotide-binding domain and the substrate-binding β domain are responsible for the significant differences in DnaK interactomes, and the nucleotide binding domain also determines the dimerization of MxDnaK2, but not MxDnaK1. Our work suggests that bacterial DnaK has been duplicated in order to deal with a more complex proteome, and that this allows evolution of distinct domains to deal with different subsets of target proteins.IMPORTANCEAll eukaryotic and ~40% of prokaryotic species encode multiple 70 kDa heat shock protein (Hsp70) homologs with similar but diversified functions. Here, we show that duplication of canonical Hsp70 (DnaK in prokaryotes) correlates with increasing proteomic complexity and evolution of particular regions of the protein. Using the Myxococcus xanthus DnaK duplicates as a case, we found that their substrate spectrums are mostly nonoverlapping, and are both consistent to that of Escherichia coli DnaK in structural and molecular characteristics, but show differential enrichment of membrane proteins. Domain/region swapping demonstrated that the nucleotide-binding domain and the β substrate-binding domain (SBDβ), but not the SBDα or disordered C-terminal tail region, are responsible for this functional divergence. This work provides the first direct evidence for regional evolution of DnaK paralogs.

The Ebola virus (EBOV) is a lethal pathogen causing hemorrhagic fever syndrome which remains a global health challenge. In the EBOV, two multifunctional proteins, VP35 and VP40, have significant roles in replication, virion assembly, and budding from the cell and have been identified as druggable targets. In this study, we employed in silico methods comprising molecular docking, molecular dynamic simulations, and pharmacological properties to identify prospective drugs for inhibiting VP35 and VP40 proteins from the myxobacterial bioactive natural product repertoire. Cystobactamid 934-2, Cystobactamid 919-1, and Cittilin A bound firmly to VP35. Meanwhile, 2-Hydroxysorangiadenosine, Enhypyrazinone B, and Sorangiadenosine showed strong binding to the matrix protein VP40. Molecular dynamic simulations revealed that, among these compounds, Cystobactamid 919-1 and 2-Hydroxysorangiadenosine had stable interactions with their respective targets. Similarly, molecular mechanics Poisson–Boltzmann surface area (MMPBSA) calculations indicated close-fitting receptor binding with VP35 or VP40. These two compounds also exhibited good pharmacological properties. In conclusion, we identified Cystobactamid 919-1 and 2-Hydroxysorangiadenosine as potential ligands for EBOV that target VP35 and VP40 proteins. These findings signify an essential step in vitro and in vivo to validate their potential for EBOV inhibition.

Microorganisms are important sources of lipolytic enzymes with characteristics for wide promising usages in the specific industrial biotechnology. The cellulolytic myxobacterium Sorangium cellulosum is rich of lipolytic enzymes in the genome, but little has been investigated. Here, we discerned 406 potential lipolytic enzymes in 13 sequenced S. cellulosum genomes. These lipolytic enzymes belonged to 12 families, and most are novel with low identities (14–37%) to those reported. We characterized a new carboxylesterase, LipB, from the alkaline-adaptive So0157-2. This enzyme, belonging to family VIII, hydrolyzed glyceryl tributyrate and p-nitrophenyl esters with short chain fatty acids (≤C12), and exhibited the highest activity against p-nitrophenyl butyrate. It retained over 50% of the activities in a broad temperature range (from 20°C to 60°C), alkaline conditions (pH 8.0–9.5), and the enzymatic activity was stable with methanol, ethanol and isopropanol, and stimulated significantly in the presence of 5 mM Ni2+. LipB also exhibited β-lactamase activity on nitrocefin, but not ampicillin, cefotaxime and imipenem. The bioinformatic analysis and specific enzymatic characteristics indicate that S. cellulosum is a promising resource to explore lipolytic enzymes for industrial adaptations.

Myxobacteria are famous for their capacity for social behavior and natural product biosynthesis. The unique sociality of myxobacteria is not only an intriguing scientific topic but also the main limiting factor for their manipulation. After more than half a century of research, a series of genetic techniques for myxobacteria have been developed, rendering these mysterious bacteria manipulable. Here, we review the advances in genetic manipulation of myxobacteria, with a particular focus on the exploitation of secondary metabolism. We emphasize the necessity and urgency of constructing the myxobacterial chassis for synthetic biology research and the exploitation of untapped secondary metabolism.

ABSTRACT Different habitats harbor different microbial communities with elusive assembly mechanisms. This study comprehensively investigated the global assembly mechanisms of microbial communities and effects of community-internal influencing factors using the Earth Microbiome Project (EMP) data set. We found that deterministic and stochastic processes contribute approximately equally to global microbial community assembly, and, specifically, deterministic processes generally play a major role in free-living and plant-associated (but not plant corpus) environments, while stochastic processes are the major contributor in animal-associated environments. In contrast with the assembly of microorganisms, the assembly of functional genes, predicted from PICRUSt, is mainly attributed to deterministic processes in all microbial communities. The sink and source microbial communities are normally assembled using similar mechanisms, and the core microorganisms are specific to different environment types. On a global scale, deterministic processes are positively related to the community alpha diversity, microbial interaction degree and bacterial predatory-specific gene abundance. Our analysis provides a panoramic picture and regularities of global and environment-typical microbial community assemblies. IMPORTANCE With the development of sequencing technologies, the research topic of microbial ecology has evolved from the analysis of community composition to community assembly, including the relative contribution of deterministic and stochastic processes for the formation and maintenance of community diversity. Many studies have reported the microbial assembly mechanisms in various habitats, but the assembly regularities of global microbial communities remain unknown. In this study, we analyzed the EMP data set using a combined pipeline to explore the assembly mechanisms of global microbial communities, microbial sources to construct communities, core microbes in different environment types, and community-internal factors influencing assembly. The results provide a panoramic picture and rules of global and environment-typical microbial community assemblies, which enhances our understandings of the mechanisms globally controlling community diversity and species coexistence.

Zhang et al. identify three hot regions of macrolide glycosyltransferase BsGT-1, that can play a role in epothilone B and sugar donor specificity. Further mutagenesis of BsGT-1 showed that specificity could be engineered and when the mutations were introduced into homologous GTs, promiscuity was improved, setting the stage for further rational glycodiversification of macrolides.

Abstract Background Polyketide synthases (PKSs) include ketone synthase (KS), acyltransferase (AT) and acyl carrier protein (ACP) domains to catalyse the elongation of polyketide chains. Some PKSs also contain ketoreductase (KR), dehydratase (DH) and enoylreductase (ER) domains as modification domains. Insertion, deletion or substitution of the catalytic domains may lead to the production of novel polyketide derivatives or to the accumulation of desired products. Epothilones are 16-membered macrolides that have been used as anticancer drugs. The substrate promiscuity of the module 4 AT domain of the epothilone PKS (EPOAT4) results in production of epothilone mixtures; substitution of this domain may change the ratios of epothilones. In addition, there are two dormant domains in module 9 of the epothilone PKS. Removing these redundant domains to generate a simpler and more efficient assembly line is a desirable goal. Results The substitution of module 4 drastically diminished the activity of epothilone PKS. However, with careful design of the KS-AT linker and the post-AT linker, replacing EPOAT4 with EPOAT2, EPOAT6, EPOAT7 or EPOAT8 (specifically incorporating methylmalonyl-CoA (MMCoA)) significantly increased the ratio of epothilone D (4) to epothilone C (3) (the highest ratio of 4:3 = 4.6:1), whereas the ratio of 4:3 in the parental strain Schlegelella brevitalea 104-1 was 1.4:1. We also obtained three strains by swapping EPOAT4 with EPOAT3, EPOAT5, or EPOAT9, which specifically incorporate malonyl-CoA (MCoA). These strains produced only epothilone C, and the yield was increased by a factor of 1.8 compared to that of parental strain 104-1. Furthermore, mutations of five residues in the AT domain identified Ser310 as the critical factor for MMCoA recognition in EPOAT4. Then, the mutation of His308 to valine or tyrosine combined with the mutation of Phe310 to serine further altered the product ratios. At the same time, we successfully deleted the inactive module 9 DH and ER domains and fused the ΨKR domain with the KR domain through an ~ 25-residue linker to generate a productive and simplified epothilone PKS. Conclusions These results suggested that the substitution and deletion of catalytic domains effectively produces desirable compounds and that selection of the linkers between domains is crucial for maintaining intact PKS catalytic activity.

Glycoside hydrolases (GHs) are greatly diverse in sequences and functions, but systematic studies of GH relationships based on structural information are lacking. Here, we report that GHs have multiple evolutionary origins and are structurally derived from 27 homologous superfamilies and 16 folds, but GHs are highly biased to distribute in a few superfamilies and folds. Six of these superfamilies are widely encoded by archaea, bacteria, and eukaryotes, indicating that they may be the most ancient in origin. Most superfamilies vary in enzyme function, and some, such as the superfamilies of (β/α)8-barrel and (α/α)6-barrel structures, exhibit extreme functional diversity; this is highly positively correlated with sequence diversity. More than one-third of glycosidase activities show a phenomenon of convergent evolution, especially the degradation functions of GHs on polysaccharides. The GHs of most superfamilies have relatively narrow environmental distributions, normally with the highest abundance in host-associated environments and a distribution preference for moderate low-temperature and acidic environments. Overall, our expanded analysis facilitates an understanding of complex GH sequence–structure–function relationships and may guide our screening and engineering of GHs.

ABSTRACT Differential transcription of functionally divergent duplicate genes is critical for bacterial cells to properly and competitively function in the environment, but the transcriptional regulation mechanisms remain in mystery. Myxococcus xanthus DK1622 possesses two duplicate groELs with divergent functions. Here, we report that MXAN_4468, an orphan gene located upstream of groEL2, encodes a response regulator (RR) and is responsible for the differential expression regulation of duplicate groELs. This RR protein realizes its negative regulatory role via a novel dual-mode functioning manner: binding to the transcription repressor HrcA to enhance its transcriptional inhibition of duplicate groELs and binding to the 3′ end of the MXAN_4468 sequence to specifically decrease the transcription of the following groEL2. Phosphorylation at the conserved 61st aspartic acid is required to trigger the regulatory functions of MXAN_4468. Pull-down experiment and mutation demonstrated that two noncognate CheA proteins, respectively belonging to the Che8 and Che7 chemosensory pathways, are involved in the protein phosphorylation. A transcriptome analysis, as well as the pull-down experiment, suggested that MXAN_4468 plays a global negative regulatory role in M. xanthus. This study elucidates, for the first time, the regulatory mechanism of differential transcription of bacterial duplicate groELs and suggests a global regulatory role of a dual-functional orphan RR. IMPORTANCE Multiply copied groELs require precise regulation of transcriptions for their divergent cellular functions. Here, we reported that an orphan response regulator (RR) tunes the transcriptional discrepancy of the duplicate groELs in Myxococcus xanthus DK1622 in a dual-functional mode. This RR protein has a conserved phosphorylation site, and the phosphorylation is required for the regulatory functions. Transcriptomic analysis, as well as a pull-down experiment, suggests that the RR plays a global regulatory role in M. xanthus. This study highlights that the dual-functional orphan RR might be involved in conducting the transcriptional symphony to stabilize the complex biological functions in cells.

Toxic effectors secreted by the type VI secretion system (T6SS) facilitate interbacterial warfare, as well as pathogenesis toward humans, animals and plants. However, systematically predicting T6SS effectors remains challenging due to their sequence and functional diversity. In this study, we systematically identified putative T6SS toxic effectors in prokaryotic genomes on the basis of the observation that genes encoding adaptor proteins and genes encoding cognate effector proteins are generally adjacent in the genome. Adaptor proteins are mediators that help to load their cognate effectors onto the T6SS spike complex. The contextual genes of the known adaptor proteins (DUF1795, DUF2169 or DUF4123) all exhibited a high proportion of encoding T6SS spike complex protein (VgrG or PAAR) and effector proteins. On the basis of the genomic context, we found that PRK06147 might be a novel adaptor protein. These four adaptors are widely distributed among the bacterial genomes. From neighbors of 5297 adaptor genes, we identified 1356 putative effector genes from 92 different families, and two-thirds were currently annotated as hypothetical proteins or as having unknown functions. Our results indicate that each class of adaptors can be used as an effective marker to identify T6SS toxic effectors, moreover, this approach can promote the discovery of new effectors.

ABSTRACT Protein toxins secreted by prokaryotes have been found to affect the pathogenicity of pathogens or directly mediate antagonistic interactions between prokaryotes. PAAR proteins are important carriers of toxic effectors and are located at the forefront of either the type VI secretion system (T6SS) or the extracellular contractile injection system (eCIS). This study systematically investigated PAAR homologues and related toxic effectors. We found that PAAR homologues were divided into 8 types and 16 subtypes and distributed in 23.1% of bacterial genomes and 7.8% of archaeal genomes. PAAR proteins of all types fold into a highly similar conical structure, even from relatively diverse underlying sequences. PAAR homologues associated with different secretion systems display a mixed phylogenetic relationship, indicating that PAAR proteins from such a subtype can be assembled on either a T6SS or an eCIS. More than 1,300 PAAR-related toxic effector genes were identified; one PAAR subtype can be associated with toxins of over 40 families, and toxins from one family can be associated with more than 10 PAAR subtypes. A large-scale comparison of Earth Microbiome Project data and prokaryotic genomes revealed that prokaryotes encoding PAAR genes are widely present in diverse environments worldwide, and taxa encoding multiple PAAR gene copies exhibit a wider distribution in environments than other taxa. Overall, our studies highlighted that PAAR proteins are versatile clips loaded with antimicrobial toxin bullets for secretion weapons (T6SS and eCIS), greatly enriching the weapon arsenal of prokaryotes, which, often together with VgrG, help prokaryotes fight for survival advantages in crowded environments. IMPORTANCE Infectious diseases caused by microbial pathogens are severe threats to human health and economic development. To respond to these threats, it is necessary to understand how microorganisms survive in and adapt to complex environments. Microorganic toxins, which are widely distributed in nature, are the key weapons in life domain interactions. PAAR proteins are important carriers of prokaryotic toxic effectors. We reveal the versatility of PAAR proteins between secretory systems and the massive diversity of toxic effectors carried by PAAR proteins, which helps prokaryotes enrich their arsenal and expand their ability to attack their neighbors. A large number of PAAR homologues and related toxic effectors enhance the survival competitiveness of prokaryotic populations. In conclusion, our work provides an example for large-scale analysis of the global distribution and ecological functions of prokaryotic functional genes.

The biosynthetic genes for secondary metabolites are often clustered into giant operons with no transcription terminator before the end. The long transcripts are frangible and the transcription efficiency declines along with the process. Internal promoters might occur in operons to coordinate the transcription of individual genes, but their effects on the transcription of operon genes and the yield of metabolites have been less investigated. Epothilones are a kind of antitumor polyketides synthesized by seven multifunctional enzymes encoded by a 56-kb operon. In this study, we identified multiple internal promoters in the epothilone operon. We performed CRISPR-dCas9–mediated transcription activation of internal promoters, combined activation of different promoters, and activation in different epothilone-producing M. xanthus strains. We found that activation of internal promoters in the operon was able to promote the gene transcription, but the activation efficiency was distinct from the activation of separate promoters. The transcription of genes in the operon was influenced by not only the starting promoter but also internal promoters of the operon; internal promoters affected the transcription of the following and neighboring upstream/downstream genes. Multiple interferences between internal promoters thus changed the transcriptional profile of operon genes and the production of epothilones. Better activation efficiency for the gene transcription and the epothilone production was obtained in the low epothilone-producing strains. Our results highlight that interactions between promoters in the operon are critical for the gene transcription and the metabolite production efficiency.

Hsp70 proteins are highly conserved chaperones that occur in all kingdoms of life. Multiple copies of Hsp70s are often present in genome-sequenced prokaryotes, especially taxa with complex life cycles, such as myxobacteria.

Summary Glycosylation of natural products can influence their pharmacological properties, and efficient glycosyltransferases (GTs) are critical for this purpose. The polyketide epothilones are potent anti‐tumour compounds, and YjiC is the only reported GT for the glycosylation of epothilone. In this study, we phylogenetically analysed 8261 GTs deposited in CAZy database and revealed that YjiC locates in a subbranch of the Macrolide I group, forming the YjiC‐subbranch with 160 GT sequences. We demonstrated that the YjiC‐subbranch GTs are normally efficient in epothilone glycosylation, but some showed low glycosylation activities. Sequence alignment of YjiC‐subbranch showed that the 66th and 77th amino acid residues, which were close to the catalytic cavity in molecular docking model, were conserved in five high‐active GTs (Q66 and P77) but changed in two low‐efficient GTs. Site‐directed residues swapping at the two positions in the two low‐active GTs (BssGT and BamGT) and the high‐active GT BsGT‐1 demonstrated that the two amino acid residues played an important role in the catalytic efficiency of epothilone glycosylation. This study highlights that the potent GTs for appointed compounds are phylogenetically grouped with conserved residues for the catalytic efficiency.

ABSTRACT Vegetation represents probably the most crucial step for the ecosystem functions of wetlands, but it is unclear how microbial populations and functions shift along with vegetation. In this study, we found that the richness and diversity of soil bacteria increased with vegetation levels and that the community composition was distinctly shifted from bare to vegetative places. The bare land displayed an extremely high abundance of Cyanobacteria as a monospecies genus, while a Gemmatimonadetes genus was predominant as multiple species in all the vegetative wetlands, suggesting their important ecosystem functions and potential mechanisms. Expression of the genes related to photosynthesis was enriched exclusively in bare land. Genes involved in biological organic carbon metabolism and the cycling of main elements (C, N, S, and P) were highly expressed in vegetative wetlands and were mostly included in the metagenome-assembled genome (MAG) of Gemmatimonadetes. Some compounds identified from soil metabolomic results also corresponded to pathways involving these key active genes. Cyanobacteria is thus responsible for the carbon sink in early infertile wetlands, and Gemmatimonadetes plays a crucial role in ecosystem functions in vegetative wetlands. Our results highlight that the soil microbial populations execute ecosystem functions for wetlands and that vegetation is the determinant for the population and functional shifts in the coastal estuarine wetland of the Yellow River Delta. IMPORTANCE Vegetation probably represents the most crucial step for the ecosystem functions of wetlands, but it is unclear how microbial populations and functions shift in pace with the colonization and succession of vegetation. In this study, we found that a Cyanobacteria monospecies genus and a Gemmatimonadetes multispecies genus are fastidiously predominant in the bare and vegetative wetlands of the Yellow River Delta, respectively. Consistently, photosynthesis genes were enriched exclusively in bare land, while genes involved in biological organic carbon metabolism and the cycling of main elements were highly expressed in vegetative wetlands, were mostly included in the MAG of Gemmatimonadetes, and were consistent with soil metabolomic results. Our results provide insight into the adaptive succession of predominant bacterial species and their ecosystem functions in response to the presence of vegetation.

Myxococcus xanthus DK1622 has two RecA genes, recA1 (MXAN_1441) and recA2 (MXAN_1388), with unknown functional differentiation. Herein, we showed that both recA genes were induced by ultraviolet (UV) irradiation but that the induction of recA1 was more delayed than that of recA2. Deletion of recA1 did not affect the growth but significantly decreased the UV-radiation survival, homologous recombination (HR) ability, and induction of LexA-dependent SOS genes. In contrast, the deletion of recA2 markedly prolonged the lag phase of bacterial growth and increased the sensitivity to DNA damage caused by hydrogen peroxide but did not change the UV-radiation resistance or SOS gene inducibility. Protein activity analysis demonstrated that RecA1, but not RecA2, catalyzed DNA strand exchange (DSE) and LexA autocleavage in vitro. Transcriptomic analysis indicated that RecA2 has evolved mainly to regulate gene expression for cellular transportation and antioxidation. This is the first report of functional divergence of duplicated bacterial recA genes. The results highlight the evolutionary strategy of M. xanthus cells for DNA HR and genome sophistication.

Abstract Background The CRISPR/dCas9 system is a powerful tool to activate the transcription of target genes in eukaryotic or prokaryotic cells, but lacks assays in complex conditions, such as the biosynthesis of secondary metabolites. Results In this study, to improve the transcription of the heterologously expressed biosynthetic genes for the production of epothilones, we established the CRISPR/dCas9-mediated activation technique in Myxococcus xanthus and analyzed some key factors involving in the CRISPR/dCas9 activation. We firstly optimized the cas9 codon to fit the M. xanthus cells, mutated the gene to inactivate the nuclease activity, and constructed the dCas9-activator system in an epothilone producer. We compared the improvement efficiency of different sgRNAs on the production of epothilones and the expression of the biosynthetic genes. We also compared the improvement effects of different activator proteins, the ω and α subunits of RNA polymerase, and the sigma factors σ54 and CarQ. By using a copper-inducible promoter, we determined that higher expressions of dCas9-activator improved the activation effects. Conclusions Our results showed that the CRISPR/dCas-mediated transcription activation is a simple and broadly applicable technique to improve the transcriptional efficiency for the production of secondary metabolites in microorganisms. This is the first time to construct the CRISPR/dCas9 activation system in myxobacteria and the first time to assay the CRISPR/dCas9 activations for the biosynthesis of microbial secondary metabolites.

Abstract Background The CRISPR/Cas9 system is a powerful tool for genome editing, in which the sgRNA binds and guides the Cas9 protein for the sequence-specific cleavage. The protocol is employable in different organisms, but is often limited by cell damage due to the endonuclease activity of the introduced Cas9 and the potential off-target DNA cleavage from incorrect guide by the 20 nt spacer. Results In this study, after resolving some critical limits, we have established an efficient CRISPR/Cas9 system for the deletion of large genome fragments related to the biosynthesis of secondary metabolites in Myxococcus xanthus cells. We revealed that the high expression of a codon-optimized cas9 gene in M. xanthus was cytotoxic, and developed a temporally high expression strategy to reduce the cell damage from high expressions of Cas9. We optimized the deletion protocol by using the tRNA–sgRNA–tRNA chimeric structure to ensure correct sgRNA sequence. We found that, in addition to the position-dependent nucleotide preference, the free energy of a 20 nt spacer was a key factor for the deletion efficiency. Conclusions By using the developed protocol, we achieved the CRISPR/Cas9-induced deletion of large biosynthetic gene clusters for secondary metabolites in M. xanthus DK1622 and its epothilone-producing mutant. The findings and the proposals described in this paper were suggested to be workable in other organisms, for example, other Gram negative bacteria with high GC content.

Abstract Background Microbial secondary metabolism is regulated by a complex and mostly-unknown network of global and pathway-specific regulators. A dozen biosynthetic gene clusters for secondary metabolites have been reported in myxobacteria, but a few regulation factors have been identified. Results We identified a transcription regulator Esi for the biosynthesis of epothilones. Inactivation of esi promoted the epothilone production, while overexpression of the gene suppressed the production. The regulation was determined to be resulted from the transcriptional changes of epothilone genes. Esi was able to bind, probably via the N-terminus of the protein, to an inverted repeat sequence in the promoter of the epothilone biosynthetic gene cluster. The Esi-homologous sequences retrieved from the RefSeq database are all of the Proteobacteria. However, the Esi regulation is not universal in myxobacteria, because the esi gene exists only in a few myxobacterial genomes. Conclusions Esi binds to the epothilone promoter and down-regulates the transcriptional level of the whole gene cluster to affect the biosynthesis of epothilone. This is the first transcription regulator identified for epothilone biosynthesis.

Although plasmids provide additional functions for cellular adaptation to the environment, they also create a metabolic burden, which causes the host cells to be less competitive with their siblings. Low-copy-number plasmids have thus evolved several mechanisms for their long-term maintenance in host cells. pMF1, discovered in Myxococcus fulvus 124B02, is the only endogenous autonomously replicated plasmid yet found in myxobacteria. Here we report that a post-segregational killing system, encoded by a co-transcriptional gene pair of pMF1.19 and pMF1.20, is involved in maintaining the pMF1 plasmid in its host cells. We demonstrate that the protein encoded by pMF1.20 is a new kind of nuclease, which is able to cleave DNA in vitro. The nuclease activity can be neutralized by the protein encoded by pMF1.19 through protein–protein interaction, suggesting that the protein is an immune protein for nuclease cleavage. We propose that the post-segregational killing mechanism of the nuclease toxin and immune protein pair encoded by pMF1.20 and pMF1.19 is helpful for the stable maintenance of pMF1 in M. fulvus cells.

Due to the high similarity in their requirements for space and food, close bacterial relatives may be each other's strongest competitors. Close bacterial relatives often form visible boundaries to separate their swarming colonies, a phenomenon termed colony-merger incompatibility. While bacterial species are known to have many incompatible strains, it is largely unclear which traits lead to multiple incompatibilities and the interactions between multiple incompatible siblings. To investigate the competitive interactions of closely related incompatible strains, we mutated Myxococcus xanthus DK1622, a predatory bacterium with complex social behavior. From 3392 random transposon mutations, we obtained 11 self-identification (SI) deficient mutants that formed unmerged colony boundaries with the ancestral strain. The mutations were at nine loci with unknown functions and formed nine independent SI mutants. Compared with their ancestral strain, most of the SI mutants showed reduced growth, swarming and development abilities, but some remained unchanged from their monocultures. When pairwise mixed with their ancestral strain for co-cultivation, these mutants exhibited improved, reduced or unchanged competitive abilities compared with the ancestral strain. The sporulation efficiencies were affected by the DK1622 partner, ranging from almost complete inhibition to 360% stimulation. The differences in competitive growth between the SI mutants and DK1622 were highly correlated with the differences in their sporulation efficiencies. However, the competitive efficiencies of the mutants in mixture were inconsistent with their growth or sporulation abilities in monocultures. We propose that the colony-merger incompatibility in M. xanthus is associated with multiple independent genetic loci, and the incompatible strains hold competitive interaction abilities, which probably determine the complex relationships between multiple incompatible M. xanthus strains and their co-existence strategies.

Chaperonin GroEL (Cpn60) requires cofactor GroES (Cpn10) for protein refolding in bacteria that possess single groEL and groES genes in a bicistronic groESL operon. Among 4,861 completely-sequenced prokaryotic genomes, 884 possess duplicate groEL genes and 770 possess groEL genes with no neighboring groES. It is unclear whether stand-alone groEL requires groES in order to function and, if required, how duplicate groEL genes and unequal groES genes balance their expressions. In Myxococcus xanthus DK1622, we determined that, while duplicate groELs were alternatively deletable, the single groES that clusters with groEL1 was essential for cell survival. Either GroEL1 or GroEL2 required interactions with GroES for in vitro and in vivo functions. Deletion of groEL1 or groEL2 resulted in decreased expressions of both groEL and groES; and ectopic complementation of groEL recovered not only the groEL but also groES expressions. The addition of an extra groES gene upstream groEL2 to form a bicistronic operon had almost no influence on groES expression and the cell survival rate, whereas over-expression of groES using a self-replicating plasmid simultaneously increased the groEL expressions. The results indicated that M. xanthus DK1622 cells coordinate expressions of the duplicate groEL and single groES genes for synergistic functions of GroELs and GroES. We proposed a potential regulation mechanism for the expression coordination.

Two unrecognizable strains of the same bacterial species form a distinct colony boundary. During growth as colonies, Myxococcus xanthus uses multiple factors to establish cooperation between recognized strains and prevent interactions with unrecognized strains of the same species. Here, ΔMXAN_0049 is a mutant strain deficient in immunity for the paired nuclease gene, MXAN_0050, that has a function in the colony-merger incompatibility of Myxococcus xanthus DK1622. With the aim to investigate the factors involved in boundary formation, a proteome and metabolome study was employed. Visualization of the boundary between DK1622 and ΔMXAN_0049 was done scanning electron microscope (SEM), which displayed the presence of many damaged cells in the boundary. Proteome analysis of the DK1622- boundary disclosed many possible proteins, such as cold shock proteins, cell shape-determining protein MreC, along with a few pathways, such as RNA degradation, phenylalanine, tyrosine and tryptophan biosynthesis, and Type VI secretion system (T6SS), which may play major roles in the boundary formation. Metabolomics studies revealed various secondary metabolites that were significantly produced during boundary formation. Overall, the results concluded that multiple factors participated in the boundary formation in M. xanthus, leading to cellular damage that is helpful in solving the mystery of the boundary formation mechanism.

Genome sequencing of cellulolytic myxobacterium Sorangium cellulosum reveals many open-reading frames (ORFs) encoding various degradation enzymes with low sequence similarity to those reported, but none of them has been characterized. In this paper, a predicted lipase gene (lipA) was cloned from S. cellulosum strain So0157-2 and characterized. lipA is 981-bp in size, encoding a polypeptide of 326 amino acids that contains the pentapeptide (GHSMG) and catalytic triad residues (Ser114, Asp250 and His284). Searching in the GenBank database shows that the LipA protein has only the 30% maximal identity to a human monoglyceride lipase. The novel lipA gene was expressed in Escherichia coli BL21 and the recombinant protein (r-LipA) was purified using Ni-NTA affinity chromatography. The enzyme hydrolyzed the p-nitrophenyl (pNP) esters of short or medium chain fatty acids (≤C10), and the maximal activity was on pNP acetate.The r-LipA is a cold-adapted lipase, with high enzymatic activity in a wide range of temperature and pH values. At 4 °C and 30 °C, the Km values of r-LipA on pNP acetate are 0.037 ± 0.001 and 0.174 ± 0.006 mM, respectively. Higher pH and temperature conditions promoted hydrolytic activity toward the pNP esters with longer chain fatty acids. Remarkably, this lipase retained much of its activity in the presence of commercial detergents and organic solvents. The results suggest that the r-LipA protein has some new characteristics potentially promising for industrial applications and S. cellulosum is an intriguing resource for lipase screening.

Myxococcus xanthus DK1622 is a rich source of novel secondary metabolites, and it is often used as an expression host of exogenous biosynthetic gene clusters. However, the frequency of obtaining large genome-deletion variants by using traditional strategies is low, and progenies generated by homologous recombination contain irregular deletions. The present study aims to develop an efficient genome-engineering system for this bacterium based on the Cre/loxP system. We first verified the functionality of the native cre system that was integrated into the chromosome with an inducible promoter PcuoA. Then we assayed the deletion frequency of 8-bp-spacer-sequence mutants in loxP by Cre recombinase which was expressed by suicide vector pBJ113 or self-replicative vector pZJY41. It was found that higher guanine content in a spacer sequence had higher deletion frequency, and the self-replicative vector was more suitable for the Cre/loxP system, probably due to the leaky expression of inducible promoter PcuoA. We also inspected the effects of different antibiotics and the native or synthetic cre gene. Polymerase chain reaction (PCR) and sequencing of new genome joints confirmed that the Cre/loxP system was able to delete a 466 kb fragment in M. xanthus. This Cre/loxP-mediated recombination could serve as an alternative genetic manipulation method.

The gene encoding the GroEL chaperonin is duplicated in nearly 30% of bacterial genomes; and although duplicated groEL genes have been comprehensively determined to have distinct physiological functions in different species, the mechanisms involved have not been characterized to date. Myxococcus xanthus DK1622 has two copies of the groEL gene, each of which can be deleted without affecting cell viability; however, the deletion of either gene does result in distinct defects in the cellular heat-shock response, predation, and development. In this study, we show that, from the expression levels of different groELs, the distinct functions of groEL1 and groEL2 in predation and development are probably the result of the substrate selectivity of the paralogous GroEL chaperonins, whereas the lethal effect of heat shock due to the deletion of groEL1 is caused by a decrease in the total groEL expression level. Following a bioinformatics analysis of the composition characteristics of GroELs from different bacteria, we performed region-swapping assays in M. xanthus, demonstrating that the differences in the apical and the C-terminal equatorial regions determine the substrate specificity of the two GroELs. Site-directed mutagenesis experiments indicated that the GGM repeat sequence at the C-terminus of GroEL1 plays an important role in functional divergence. Divergent functions of duplicated GroELs, which have similar patterns of variation in different bacterial species, have thus evolved mainly via alteration of the apical and the C-terminal equatorial regions. We identified the specific substrates of strain DK1622's GroEL1 and GroEL2 using immunoprecipitation and mass spectrometry techniques. Although 68 proteins bound to both GroEL1 and GroEL2, 83 and 46 proteins bound exclusively to GroEL1 or GroEL2, respectively. The GroEL-specific substrates exhibited distinct molecular sizes and secondary structures, providing an encouraging indication for GroEL evolution for functional divergence.

The bacterial type IV pilus (T4P) is the strongest biological motor known to date as its retraction can generate forces well over 100 pN. Myxococcus xanthus, a δ-proteobacterium, provides a good model for T4P investigations because its social (S) gliding motility is powered by T4P. In this study, the interactions among M. xanthus T4P proteins were investigated using genetics and the yeast two-hybrid (Y2H) system. Our genetic analysis suggests that there is an integrated T4P structure that crosses the inner membrane (IM), periplasm and the outer membrane (OM). Moreover, this structure exists in the absence of the pilus filament. A systematic Y2H survey provided evidence for direct interactions among IM and OM proteins exposed to the periplasm. For example, the IM lipoprotein PilP interacted with its cognate OM protein PilQ. In addition, interactions among T4P proteins from the thermophile Thermus thermophilus were investigated by Y2H. The results indicated similar protein-protein interactions in the T4P system of this non-proteobacterium despite significant sequence divergence between T4P proteins in T. thermophilus and M. xanthus. The observations here support the model of an integrated T4P structure in the absence of a pilus in diverse bacterial species.

Migrating cells employ sophisticated signal transduction systems to respond to their environment and polarize towards attractant sources. Bacterial cells also regulate their polarity dynamically to reverse their direction of movement. In Myxococcus xanthus, a GTP-bound Ras-like G-protein, MglA, activates the motility machineries at the leading cell pole. Reversals are provoked by pole-to-pole switching of MglA, which is under the control of a chemosensory-like signal transduction cascade (Frz). It was previously known that the asymmetric localization of MglA at one cell pole is regulated by MglB, a GTPase Activating Protein (GAP). In this process, MglB specifically localizes at the opposite lagging cell pole and blocks MglA localization at that pole. However, how MglA is targeted to the leading pole and how Frz activity switches the localizations of MglA and MglB synchronously remained unknown. Here, we show that MglA requires RomR, a previously known response regulator protein, to localize to the leading cell pole efficiently. Specifically, RomR-MglA and RomR-MglB complexes are formed and act complementarily to establish the polarity axis, segregating MglA and MglB to opposite cell poles. Finally, we present evidence that Frz signaling may regulate MglA localization through RomR, suggesting that RomR constitutes a link between the Frz-signaling and MglAB polarity modules. Thus, in Myxococcus xanthus, a response regulator protein governs the localization of a small G-protein, adding further insight to the polarization mechanism and suggesting that motility regulation evolved by recruiting and combining existing signaling modules of diverse origins.

pMF1 is the only autonomously replicating plasmid that has been recently identified in myxobacteria. This study characterized the partitioning (par) system of this plasmid. The fragment that significantly increased the retaining stability of plasmids in Myxococcus cells in the absence of selective antibiotics contained three open reading frames (ORFs) pMF1.21-pMF1.23 (parCAB). The pMF1.22 ORF (parA) is homologous to members of the parA ATPase family, with the highest similarity (56%) to the Sphingobium japonicum ParA-like protein, while the other two ORFs had no homologs in GenBank. DNase I footprinting and electrophoretic mobility shift assays showed that the pMF1.23 (parB) product is a DNA-binding protein of iteron DNA sequences, while the product of pMF1.21 (parC) has no binding activity but is able to enhance the DNA-binding activity of ParB to iterons. The ParB protein autogenously repressed the expression of the par genes, consistent with the type Ib par pattern, while the ParC protein has less repressive activity. The ParB-binding iteron sequences are distributed not only near the partitioning gene loci but also along pMF1. These results indicate that the pMF1 par system has novel structural and functional characteristics.

In the title compound, [Sn2(CH3)6(C6H2O4S)(H2O)]n, each of the two crystallographically independent Sn atoms exhibits a distorted trigonal–bipyramidal coordination geometry formed by two O and three C atoms. The coordinated water molecule plays an important role in crystal packing consolidation via O—H...O hydrogen bonding.