Operon Structure Research Articles

Improved BioCyc Operon Prediction: Revisiting the Operon Prediction Problem.

Operon prediction is a valuable component of microbial-genome annotation because operon organization can yield inferences about gene function, and because knowledge of operon structure can aid the interpretation of gene expression data. We present a number of improvements to the existing Pathway Tools operon predictor based mostly on 7 new features that we hypothesized would increase its performance. The new features include shared Gene Ontology biological process terms, similarity of codon usage and GC content, correlated gene expression, and shared protein complex. We evaluated the proposed 7 new features and found that the addition of 6 of them improved the performance of the operon predictor from 79.55% to 83.49%, a decrease in error rate of 19.3%. When gene expression data was not included, the accuracy decreased to 82.547, still an improvement of 14.7%. One of the proposed features as well as a previously used feature had no effect and were removed. Although some of the new features had strong predictive value individually, when combined with the other features they did not have a large impact on predictive accuracy, suggesting that they were not independent from the other features.

De novo biosynthesis of 3-hydroxy-3-methylbutyrate as anti-catabolic supplement by metabolically engineered Escherichia coli

3-Hydroxy-3-methylbutyrate (HMB) is a five-carbon branch-chain hydroxy acid currently used as a dietary supplement to treat sarcopenia and exercise training. However, its current production relies on conventional chemical processes which require toxic substances and are generally non-sustainable. While bio-based syntheses of HMB have been developed, they are dependent on biotransformation of its direct precursors which are generally costly. Therefore, in this work, we developed a synthetic de novo HMB biosynthetic pathway that enables HMB production from renewable resources. This novel HMB biosynthesis employs heterologous enzymes from mevalonate pathway and myxobacterial iso-fatty acid pathway for converting acetyl-CoA to HMB-CoA. Subsequently, HMB-CoA is hydrolyzed by a thioesterase to yield HMB. Upon expression of this pathway, our initial Escherichia coli strain produced 660 mg/L of HMB from glucose in 48 hours. Through optimization of coenzyme A removal from HMB-CoA and genetic operon structure, our final strain achieved HMB production titer of 17.7 g/L in glucose minimal media using a bench-top bioreactor. This engineered strain was further demonstrated to produce HMB from other renewable carbon sources such as xylose, glycerol, and acetate. The results from this work provided a flexible and environmentally benign method for producing HMB.

Azithromycin resistance in Escherichia coli and Salmonella from food-producing animals and meat in Europe.

To characterize the genetic basis of azithromycin resistance in Escherichia coli and Salmonella collected within the EU harmonized antimicrobial resistance (AMR) surveillance programme in 2014-18 and the Danish AMR surveillance programme in 2016-19. WGS data of 1007 E. coli [165 azithromycin resistant (MIC > 16 mg/L)] and 269 Salmonella [29 azithromycin resistant (MIC > 16 mg/L)] were screened for acquired macrolide resistance genes and mutations in rplDV, 23S rRNA and acrB genes using ResFinder v4.0, AMRFinder Plus and custom scripts. Genotype-phenotype concordance was determined for all isolates. Transferability of mef(C)-mph(G)-carrying plasmids was assessed by conjugation experiments. mph(A), mph(B), mef(B), erm(B) and mef(C)-mph(G) were detected in E. coli and Salmonella, whereas erm(C), erm(42), ere(A) and mph(E)-msr(E) were detected in E. coli only. The presence of macrolide resistance genes, alone or in combination, was concordant with the azithromycin-resistant phenotype in 69% of isolates. Distinct mph(A) operon structures were observed in azithromycin-susceptible (n = 50) and -resistant (n = 136) isolates. mef(C)-mph(G) were detected in porcine and bovine E. coli and in porcine Salmonella enterica serovar Derby and Salmonella enterica 1,4, [5],12:i:-, flanked downstream by ISCR2 or TnAs1 and associated with IncIγ and IncFII plasmids. Diverse azithromycin resistance genes were detected in E. coli and Salmonella from food-producing animals and meat in Europe. Azithromycin resistance genes mef(C)-mph(G) and erm(42) appear to be emerging primarily in porcine E. coli isolates. The identification of distinct mph(A) operon structures in susceptible and resistant isolates increases the predictive power of WGS-based methods for in silico detection of azithromycin resistance in Enterobacterales.

SAFPred: synteny-aware gene function prediction for bacteria using protein embeddings.

Today, we know the function of only a small fraction of the protein sequences predicted from genomic data. This problem is even more salient for bacteria, which represent some of the most phylogenetically and metabolically diverse taxa on Earth. This low rate of bacterial gene annotation is compounded by the fact that most function prediction algorithms have focused on eukaryotes, and conventional annotation approaches rely on the presence of similar sequences in existing databases. However, often there are no such sequences for novel bacterial proteins. Thus, we need improved gene function prediction methods tailored for bacteria. Recently, transformer-based language models-adopted from the natural language processing field-have been used to obtain new representations of proteins, to replace amino acid sequences. These representations, referred to as protein embeddings, have shown promise for improving annotation of eukaryotes, but there have been only limited applications on bacterial genomes. To predict gene functions in bacteria, we developed SAFPred, a novel synteny-aware gene function prediction tool based on protein embeddings from state-of-the-art protein language models. SAFpred also leverages the unique operon structure of bacteria through conserved synteny. SAFPred outperformed both conventional sequence-based annotation methods and state-of-the-art methods on multiple bacterial species, including for distant homolog detection, where the sequence similarity to the proteins in the training set was as low as 40%. Using SAFPred to identify gene functions across diverse enterococci, of which some species are major clinical threats, we identified 11 previously unrecognized putative novel toxins, with potential significance to human and animal health. https://github.com/AbeelLab/safpred.

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species.

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

Optimization of degradation conditions and analysis of degradation mechanism for nitrite by Bacillus aryabhattai 47

Excessive nitrite levels cause significant damage to aquaculture, making it crucial to explore green and reliable nitrite removal technologies. In this study, A Bacillus aryabhattai (designated as the strain 47) isolated from aquaculture wastewater was used as the experimental strain. The nitrite degradation conditions of the strain 47 were optimized, and the optimal conditions are: glucose was 12.74 g/L, fermented special soybean meal was 21.27 g/L, MgCl2 369 mg/L, pH 7.0, incubated at 30 °C with the inoculum size of 2 % and the rotation speed of 170 rpm. Under the optimal conditions, the nitrite concentration of the culture solution was 200 mg/L, and the nitrite removal rate reached 91.4 %. Meanwhile, the mechanism by which Mg2+ enhanced the nitrite degradation ability of the strain 47 was investigated by transcriptomics. An operon structure directed cellular trafficking of Mg2+, and then, the Mg2+-mediated catalytic reaction of multiple enzymes enhanced and improved cellular metabolic processes (e.g. the transport and metabolism of nitrite, central carbohydrate metabolism oxidative phosphorylation). At the same time, with the progress of cell metabolism, cells secreted a series of enzymes related to nitrite transport and metabolism to promote the metabolism of nitrite. And the process of the assimilated nitrate reduction pathway of nitrite degradation in the strain 47 was elaborated at the transcriptome level. This study provided a new insight into nitrite treatment mediated by microbial organisms.

Cross-evaluation of E. coli’s operon structures via a whole-cell model suggests alternative cellular benefits for low- versus high-expressing operons

Many bacteria use operons to coregulate genes, but it remains unclear how operons benefit bacteria. We integrated E.coli's 788 polycistronic operons and 1,231 transcription units into an existing whole-cell model and found inconsistencies between the proposed operon structures and the RNA-seq read counts that the model was parameterized from. We resolved these inconsistencies through iterative, model-guided corrections to both datasets, including the correction of RNA-seq counts of short genes that were misreported as zero by existing alignment algorithms. The resulting model suggested two main modes by which operons benefit bacteria. For 86% of low-expression operons, adding operons increased the co-expression probabilities of their constituent proteins, whereas for 92% of high-expression operons, adding operons resulted in more stable expression ratios between the proteins. These simulations underscored the need for further experimental work on how operons reduce noise and synchronize both the expression timing and the quantity of constituent genes. A record of this paper's transparent peer review process is included in the supplemental information.

Retracing the evolution of a modern periplasmic binding protein.

Investigating the evolution of structural features in modern multidomain proteins helps to understand their immense diversity and functional versatility. The class of periplasmic binding proteins (PBPs) offers an opportunity to interrogate one of the main processes driving diversification: the duplication and fusion of protein sequences to generate new architectures. The symmetry of their two-lobed topology, their mechanism of binding, and the organization of their operon structure led to the hypothesis that PBPs arose through a duplication and fusion event of a single common ancestor. To investigate this claim, we set out to reverse the evolutionary process and recreate the structural equivalent of a single-lobed progenitor using ribose-binding protein (RBP) as our model. We found that this modern PBP can be deconstructed into its lobes, producing two proteins that represent possible progenitor halves. The isolated halves of RBP are well folded and monomeric proteins, albeit with a lower thermostability, and do not retain the original binding function. However, the two entities readily form a heterodimer in vitro and in-cell. The x-ray structure of the heterodimer closely resembles the parental protein. Moreover, the binding function is fully regained upon formation of the heterodimer with a ligand affinity similar to that observed in the modern RBP. This highlights how a duplication event could have given rise to a stable and functional PBP-like fold and provides insights into how more complex functional structures can evolve from simpler molecular components.

The Broken Chloroplast Gene Clusters in Gymnosperms Exhibit Elevated Substitution Rates

Plant chloroplast (cp) gene clusters consist of genes arranged closely together on the cp genome. These genes are organized in operon structures and participate in cotranscription, typically exhibiting conservation. Broken gene clusters have been observed in gymnosperms. In order to investigate whether the substitution rates and selection pressure of associated genes are affected following the disruption of gene clusters, the cp genomes of 80 species (78 gymnosperms and 2 outgroups) were analyzed. A phylogenetic analysis was conducted using 58 shared genes to examine the evolutionary rates and selection pressure of genes associated with gene clusters and protein-coding genes in Sciadopitys verticillata. The results demonstrate that S. verticillata exhibited the highest number of rearrangements compared to the Cycas revoluta genome. Four gene clusters (rps2, psbB, rpoB, and petL clusters) in S. verticillata were disrupted, while rps2 in Callitris rhomboidea experienced disruption. Significantly increased evolutionary rates were observed in 12 out of 18 gene cluster-related genes in S. verticillata. Following disruption, S. verticillata and C. rhomboidea exhibited an increase in gene cluster-related genes, particularly rps2, and higher selection pressure on both rps2 and atpA genes compared to other species. Furthermore, among the 58 genes shared by S. verticillata, the evolutionary rates of 36 genes increased, and the selection pressure on 13 genes exceeded that of other species. These results indicate an increased substitution rate of gene clusters in S. verticillata and C. rhomboidea. The large-scale rearrangement and elevated substitution rates of the cp genome in S. verticillata were revealed. This study sheds light on the heterogeneity of cp genome evolution in gymnosperms.

Phylogenetic investigation of Gammaproteobacteria proteins involved in exogenous long-chain fatty acid acquisition and assimilation

BackgroundThe incorporation of exogenous fatty acids into the cell membrane yields structural modifications that directly influence membrane phospholipid composition and indirectly contribute to virulence. FadL and FadD are responsible for importing and activating exogenous fatty acids, while acyltransferases (PlsB, PlsC, PlsX, PlsY) incorporate fatty acids into the cell membrane. Many Gammaproteobacteria species possess multiple homologs of these proteins involved in exogenous fatty acid metabolism, suggesting the evolutionary acquisition and maintenance of this transport pathway. MethodsThis study developed phylogenetic trees based on amino acid and nucleotide sequences of homologs of FadL, FadD, PlsB, PlsC, PlsX, and PlsY via Mr. Bayes and RAxML algorithms. We also explored the operon arrangement of genes encoding for FadL. Additionally, FadL homologs were modeled via SWISS-MODEL, validated and refined by SAVES, Galaxy Refine, and GROMACS, and docked with fatty acids via AutoDock Vina. Resulting affinities were analyzed by 2-way ANOVA test and Tukey's post-hoc test. ResultsOur phylogenetic trees revealed grouping based on operon structure, original homolog blasted from, and order of the homolog, suggesting a more ancestral origin of the multiple homolog phenomena. Our molecular docking simulations indicated a similar binding pattern for the fatty acids between the different FadL homologs. General significanceOur study is the first to illustrate the phylogeny of these proteins and to investigate the binding of various FadL homologs across orders with fatty acids. This study helps unravel the mystery surrounding these proteins and presents topics for future research.

Host range of strand-biased circularizing integrative elements: a new class of mobile DNA elements nesting in Gammaproteobacteria

BackgroundThe strand-biased circularizing integrative elements (SEs) are putatively non-mobilizable integrative elements for transmitting antimicrobial resistance genes. The transposition mode and the prevalence of SEs in prokaryotes remain vague.ResultsTo corroborate the transposition mode and the prevalence of SEs, hypothetical transposition intermediates of an SE were searched for in genomic DNA fractions of an SE host. Then, the SE core genes were defined based on gene knockout experiments, and the synteny blocks of their distant homologs were searched for in the RefSeq complete genome sequence database using PSI-BLAST. A genomic DNA fractionation experiment revealed that SE copies are present in a double-stranded nicked circular form in vivo. Operonic structure of three conserved coding sequences (intA, tfp, intB) and srap located at the left end of SEs were identified as essential for attL × attR recombination. The synteny blocks of tfp and srap homologs were detected in 3.6% of the replicons of Gammaproteobacteria but not in other taxa, implying that SE movement is host-dependent. SEs have been discovered most frequently in the orders Vibrionales (19% of replicons), Pseudomonadales (18%), Alteromonadales (17%), and Aeromonadales (12%). Genomic comparisons revealed 35 new SE members with identifiable termini. SEs are present at 1 to 2 copies per replicon and have a median length of 15.7 kb. Three newly identified SE members carry antimicrobial resistance genes, like tmexCD-toprJ, mcr-9, and blaGMA-1. Further experiments validated that three new SE members possess the strand-biased attL × attR recombination activity.ConclusionsThis study suggested that transposition intermediates of SEs are double-stranded circular DNA. The main hosts of SEs are a subset of free-living Gammaproteobacteria; this represents a rather narrow host range compared to those of mobile DNA element groups discovered to date. As the host range, genetic organization, and movements are unique among the mobile DNA elements, SEs provide a new model system for host-mobile DNA element coevolution studies.

Deep Profiling of the Proteome Dynamics of Pseudomonas aeruginosa Reference Strain PAO1 under Different Growth Conditions.

As one of the most common bacterial pathogens causing nosocomial infections, Pseudomonas aeruginosa is highly adaptable to survive under various conditions. Here, we profiled the abundance dynamics of 3489 proteins across different growth stages in the P. aeruginosa reference strain PAO1 using data-independent acquisition-based quantitative proteomics. The proteins differentially expressed during the planktonic growth exhibit several distinct patterns of expression profiles and are relevant to various biological processes, highlighting the continuous adaptation of the PAO1 proteome during the transition from the acceleration phase to the stationary phase. By contrasting the protein expressions in a biofilm to planktonic cells, the known roles of T6SS, phenazine biosynthesis, quorum sensing, and c-di-GMP signaling in the biofilm formation process were confirmed. Additionally, we also discovered several new functional proteins that may play roles in the biofilm formation process. Lastly, we demonstrated the general concordance of protein expressions within operons across various growth states, which permits the study of coexpression protein units, and reversely, the study of regulatory components in the operon structure. Taken together, we present a high-quality and valuable resource on the proteomic dynamics of the P. aeruginosa reference strain PAO1, with the potential of advancing our understanding of the overall physiology of Pseudomonas bacteria.

Full Copper Resistance in Cupriavidus metallidurans Requires the Interplay of Many Resistance Systems.

The metal-resistant bacterium Cupriavidus metallidurans uses its copper resistance components to survive the synergistic toxicity of copper ions and gold complexes in auriferous soils. The cup, cop, cus, and gig determinants encode as central component the Cu(I)-exporting PIB1-type ATPase CupA, the periplasmic Cu(I)-oxidase CopA, the transenvelope efflux system CusCBA, and the Gig system with unknown function, respectively. The interplay of these systems with each other and with glutathione (GSH) was analyzed. Copper resistance in single and multiple mutants up to the quintuple mutant was characterized in dose-response curves, Live/Dead-staining, and atomic copper and glutathione content of the cells. The regulation of the cus and gig determinants was studied using reporter gene fusions and in case of gig also RT-PCR studies, which verified the operon structure of gigPABT. All five systems contributed to copper resistance in the order of importance: Cup, Cop, Cus, GSH, and Gig. Only Cup was able to increase copper resistance of the Δcop Δcup Δcus Δgig ΔgshA quintuple mutant but the other systems were required to increase copper resistance of the Δcop Δcus Δgig ΔgshA quadruple mutant to the parent level. Removal of the Cop system resulted in a clear decrease of copper resistance in most strain backgrounds. Cus cooperated with and partially substituted Cop. Gig and GSH cooperated with Cop, Cus, and Cup. Copper resistance is thus the result of an interplay of many systems. IMPORTANCE The ability of bacteria to maintain homeostasis of the essential-but-toxic "Janus"-faced element copper is important for their survival in many natural environments but also in case of pathogenic bacteria in their respective host. The most important contributors to copper homeostasis have been identified in the last decades and comprise PIB1-type ATPases, periplasmic copper- and oxygen-dependent copper oxidases, transenvelope efflux systems, and glutathione; however, it is not known how all these players interact. This publication investigates this interplay and describes copper homeostasis as a trait emerging from a network of interacting resistance systems.

Molecular evolution and population genetics of glutamate decarboxylase acid resistance pathway in lactic acid bacteria.

Glutamate decarboxylase (GAD) pathway (GDP) is a major acid resistance mechanism enabling microorganisms' survival in low pH environments. We aimed to study the molecular evolution and population genetics of GDP in Lactic Acid Bacteria (LAB) to understand evolutionary processes shaping adaptation to acidic environments comparing species where the GDP genes are organized in an operon structure (Levilactobacillus brevis) versus lack of an operon structure (Lactiplantibacillus plantarum). Within species molecular population genetic analyses of GDP genes in L. brevis and L. plantarum sampled from diverse fermented food and other environments showed abundant synonymous and non-synonymous nucleotide diversity, mostly driven by low frequency changes, distributed throughout the coding regions for all genes in both species. GAD genes showed higher level of replacement polymorphism compared to transporter genes (gadC and YjeM) for both species, and GAD genes that are outside of an operon structure showed even higher level of replacement polymorphism. Population genetic tests suggest negative selection against replacement changes in all genes. Molecular structure and amino acid characteristics analyses showed that in none of the GDP genes replacement changes alter 3D structure or charge distribution supporting negative selection against non-conservative amino acid changes. Phylogenetic and between species divergence analyses suggested adaptive protein evolution on GDP genes comparing phylogenetically distant species, but conservative evolution comparing closely related species. GDP genes within an operon structure showed slower molecular evolution and higher conservation. All GAD and transporter genes showed high codon usage bias in examined LAB species suggesting high expression and utilization of acid resistance genes. Substantial discordances between species, GAD, and transporter gene tree topologies were observed suggesting molecular evolution of GDP genes do not follow speciation events. Distribution of operon structure on the species tree suggested multiple independent gain or loss of operon structure in LABs. In conclusion, GDP genes in LABs exhibit a dynamic molecular evolutionary history shaped by gene loss, gene transfer, negative and positive selection to maintain its active role in acid resistance mechanism, and enable organisms to thrive in acidic environments.

Encapsulin cargo loading: progress and potential.

Encapsulins are a recently discovered class of prokaryotic self-assembling icosahedral protein nanocompartments measuring between 24 and 42 nm in diameter, capable of selectively encapsulating dedicated cargo proteins in vivo. They have been classified into four families based on sequence identity and operon structure, and thousands of encapsulin systems have recently been computationally identified across a wide range of bacterial and archaeal phyla. Cargo encapsulation is mediated by the presence of specific targeting motifs found in all native cargo proteins that interact with the interior surface of the encapsulin shell during self-assembly. Short C-terminal targeting peptides (TPs) are well documented in Family 1 encapsulins, while more recently, larger N-terminal targeting domains (TDs) have been discovered in Family 2. The modular nature of TPs and their facile genetic fusion to non-native cargo proteins of interest has made cargo encapsulation, both in vivo and in vitro, readily exploitable and has therefore resulted in a range of rationally engineered nano-compartmentalization systems. This review summarizes current knowledge on cargo protein encapsulation within encapsulins and highlights select studies that utilize TP fusions to non-native cargo in creative and useful ways.

Genome-based analyses to learn from and about Paenibacillus sonchi genomovar Riograndensis SBR5T.

Paenibacillus sonchi genomovar Riograndensis SBR5T is a plant growth-promoting rhizobacterium (PGPR) isolated in the Brazilian state of Rio Grande do Sul from the rhizosphere of Triticum aestivum. It fixes nitrogen, produces siderophores as well as the phytohormone indole-3-acetic acid, solubilizes phosphate and displays antagonist activity against Listeria monocytogenes and Pectobacterium carotovorum. Comprehensive omics analysis and the development of genetic tools are key to characterizing and engineering such non-model microorganisms. Therefore, the complete genome of SBR5T was sequenced, and shown to encode 6,705 proteins, 87 tRNAs, and 27 rRNAs and it enabled a landscape transcriptome analysis that unveiled conserved transcriptional and translational patterns and characterized operon structures and riboswitches. The pangenome of P. sonchi species is open with a stable core pangenome. At the same time, the analysis of genes coding for nitrogenases revealed that the trait of nitrogen fixation is sparse within the Paenibacillaceae family and the presence of Fe-only nitrogenase in the P. sonchi group was exclusive to SBR5T. The development of genetic tools for SBR5T enabled genetic transformation, plasmid construction for constitutive and inducible gene expression, and gene repression using the CRISPRi system. Altogether, the work with P. sonchi can guide the study of non-model bacteria with economic potential.

Structural and functional mapping of ars gene cluster in Deinococcus indicus DR1

Deinococcus indicus DR1 is a novel Gram-negative bacterium, isolated from the Dadri wetlands in Uttar Pradesh, India. In addition to being radiation-resistant, the rod-shaped, red-pigmented organism shows extraordinary resistance to arsenic. The proteins of the corresponding ars gene cluster involved in arsenic extrusion in D. indicus DR1 have not yet been characterized. Additionally, how these proteins regulate each other providing arsenic resistance is still unclear. Here, we present a computational model of the operonic structure and the corresponding characterization of the six proteins of the ars gene cluster in D. indicus DR1. Additionally, we show the expression of the genes in the presence of arsenic using qRT-PCR. The ars gene cluster consists of two transcriptional regulators (ArsR1, ArsR2), two arsenate reductases (ArsC2, ArsC3), one metallophosphatase family protein (MPase), and a transmembrane arsenite efflux pump (ArsB). The transcriptional regulators are trans-acting repressors, and the reductases reduce arsenate (As5+) ions to arsenite (As3+) ions for favourable extrusion. The proteins modelled using RoseTTAFold, and their conformationally stable coordinates obtained after MD simulation indicate their various functional roles with respect to arsenic. Excluding ArsB, all the proteins belong to the α + β class of proteins. ArsB, being a membrane protein, is fully α-helical, with 12 transmembrane helices. The results show the degree of similarity or divergence of the mechanism utilized by these proteins of ars gene cluster in D. indicus DR1 to confer high levels of arsenic tolerance. This structural characterization study of the ars genes will enable new and deeper insights of arsenic tolerance.

ClbP Gene, a Potential New Member of the β-Lactamase Family.

The colibactin island (pks) of Escherichia coli formed by 19 genes (55-Kb), encodes non-ribosomal peptide (NRP) and polyketide (PK) synthases, which allow the synthesis of colibactin, a suspected hybrid PK-NRP compound that causes damage to DNA in eukaryotic cells. The clbP, an unusual essential gene, is found in the operon structure with the clbS gene in the pks-encoded machinery. Interestingly, the clbP gene has been annotated as a β-lactamase but no previous study has reported its β-lactamase characteristics. In this study, we (i) investigated the β-lactamase properties of the clbP gene in silico by analysing its phylogenetic relationship with bacterial β-lactamase and peptidase enzymes, (ii) compared its three-dimensional (3D) protein structure with those of bacterial β-lactamase proteins using the Phyr2 database and PyMOL software, and (iii) evaluated in vitro its putative enzymatic activities, including β-lactamase, nuclease, and ribonuclease using protein expression and purification from an E. coli BL21 strain. In this study, we reveal a structural configuration of toxin/antitoxin systems in this island. Thus, similar to the toxin/antitoxin systems, the role of the clbP gene within the pks-island gene group appears as an antitoxin, insofar as it is responsible for the activation of the toxin, which is colibactin. In silico, our analyses revealed that ClbP belonged to the superfamily of β-lactamase, class C. Furthermore, in vitro we were unable to demonstrate its β-lactamase activity, likely due to the fact that the clbP gene requires co-expression with other genes, such as the genes present in the pks-island (19 genes). More research is needed to better understand its actions, particularly with regards to antibiotics, and to discover whether it has any additional functions due to the importance of this gene and its toxicity.

A putative bifunctional CPD/ (6-4) photolyase from the cyanobacteria Synechococcus sp. PCC 7335 is encoded by a UV-B inducible operon: New insights into the evolution of photolyases.

Photosynthetic organisms are continuously exposed to solar ultraviolet radiation-B (UV-B) because of their autotrophic lifestyle. UV-B provokes DNA damage, such as cyclobutane pyrimidine dimers (CPD) or pyrimidine (6-4) pyrimidone photoproducts (6-4 PPs). The cryptochrome/photolyase family (CPF) comprises flavoproteins that can bind damaged or undamaged DNA. Photolyases (PHRs) are enzymes that repair either CPDs or 6-4 PPs. A natural bifunctional CPD/(6-4)- PHR (PhrSph98) was recently isolated from the UV-resistant bacteria Sphingomonas sp. UV9. In this work, phylogenetic studies of bifunctional CPD/(6-4)- photolyases and their evolutionary relationship with other CPF members were performed. Amino acids involved in electron transfer and binding to FAD cofactor and DNA lesions were conserved in proteins from proteobacteria, planctomycete, bacteroidete, acidobacteria and cyanobacteria clades. Genome analysis revealed that the cyanobacteria Synechococcus sp. PCC 7335 encodes a two-gene assembly operon coding for a PHR and a bifunctional CPD/(6-4) PHR- like. Operon structure was validated by RT-qPCR analysis and the polycistronic transcript accumulated after 15 min of UV-B irradiation. Conservation of structure and evolution is discussed. This study provides evidence for a UV-B inducible PHR operon that encodes a CPD/(6-4)- photolyase homolog with a putative bifunctional role in the repair of CPDs and 6-4 PPs damages in oxygenic photosynthetic organisms.

Expression of divergent methyl/alkyl coenzyme M reductases from uncultured archaea

Methanogens and anaerobic methane-oxidizing archaea (ANME) are important players in the global carbon cycle. Methyl-coenzyme M reductase (MCR) is a key enzyme in methane metabolism, catalyzing the last step in methanogenesis and the first step in anaerobic methane oxidation. Divergent mcr and mcr-like genes have recently been identified in uncultured archaeal lineages. However, the assembly and biochemistry of MCRs from uncultured archaea remain largely unknown. Here we present an approach to study MCRs from uncultured archaea by heterologous expression in a methanogen, Methanococcus maripaludis. Promoter, operon structure, and temperature were important determinants for MCR production. Both recombinant methanococcal and ANME-2 MCR assembled with the host MCR forming hybrid complexes, whereas tested ANME-1 MCR and ethyl-coenzyme M reductase only formed homogenous complexes. Together with structural modeling, this suggests that ANME-2 and methanogen MCRs are structurally similar and their reaction directions are likely regulated by thermodynamics rather than intrinsic structural differences.