Metagenomic analysis of carbohydrate-active enzymes and their contribution to marine sediment biodiversity.

Marine sediments constitute the world's most substantial long-term carbon repository. The microorganisms dwelling in these sediments mediate the transformation of fixed oceanic carbon, but their contribution to the carbon cycle is not fully understood. Previous culture-independent investigations into sedimentary microorganisms have underscored the significance of carbohydrates in the carbon cycle. In this study, we employ a metagenomic methodology to investigate the distribution and abundance of carbohydrate-active enzymes (CAZymes) in 37 marine sediments sites. These sediments exhibit varying oxygen availability and were isolated in diverse regions worldwide. Our comparative analysis is based on the metabolic potential for oxygen utilisation, derived from genes present in both oxic and anoxic environments. We found that extracellular CAZyme modules targeting the degradation of plant and algal detritus, necromass, and host glycans were abundant across all metagenomic samples. The analysis of these results indicates that the oxic/anoxic conditions not only influence the taxonomic composition of the microbial communities, but also affect the occurrence of CAZyme modules involved in the transformation of necromass, algae and plant detritus. To gain insight into the sediment microbial taxa, we reconstructed metagenome assembled genomes (MAG) and examined the presence of primary extracellular carbohydrate active enzyme (CAZyme) modules. Our findings reveal that the primary CAZyme modules and the CAZyme gene clusters discovered in our metagenomes were prevalent in the Bacteroidia, Gammaproteobacteria, and Alphaproteobacteria classes. We compared those MAGs to organisms from the same taxonomic classes found in soil, and we found that they were similar in its CAZyme repertoire, but the soil MAG contained a more abundant and diverse CAZyme content. Furthermore, the data indicate that abundant classes in our metagenomic samples, namely Alphaproteobacteria, Bacteroidia and Gammaproteobacteria, play a pivotal role in carbohydrate transformation within the initial few metres of the sediments.

Glucose consumption rate-dependent transcriptome profiling of Escherichia coli provides insight on performance as microbial factories.

BACKGROUND: The modification of glucose import capacity is an engineering strategy that has been shown to improve the characteristics of Escherichia coli as a microbial factory. A reduction in glucose import capacity can have a positive effect on production strain performance, however, this is not always the case. In this study, E. coli W3110 and a group of four isogenic derivative strains, harboring single or multiple deletions of genes encoding phosphoenolpyruvate:sugar phosphotransferase system (PTS)-dependent transporters as well as non-PTS transporters were characterized by determining their transcriptomic response to reduced glucose import capacity. RESULTS: These strains were grown in bioreactors with M9 mineral salts medium containing 20 g/L of glucose, where they displayed specific growth rates ranging from 0.67 to 0.27 h-1, and specific glucose consumption rates (qs) ranging from 1.78 to 0.37 g/g h. RNA-seq analysis revealed a transcriptional response consistent with carbon source limitation among all the mutant strains, involving functions related to transport and metabolism of alternate carbon sources and characterized by a decrease in genes encoding glycolytic enzymes and an increase in gluconeogenic functions. A total of 107 and 185 genes displayed positive and negative correlations with qs, respectively. Functions displaying positive correlation included energy generation, amino acid biosynthesis, and sugar import. CONCLUSION: Changes in gene expression of E. coli strains with impaired glucose import capacity could be correlated with qs values and this allowed an inference of the physiological state of each mutant. In strains with lower qs values, a gene expression pattern is consistent with energy limitation and entry into the stationary phase. This physiological state could explain why these strains display a lower capacity to produce recombinant protein, even when they show very low rates of acetate production. The comparison of the transcriptomes of the engineered strains employed as microbial factories is an effective approach for identifying favorable phenotypes with the potential to improve the synthesis of biotechnological products.

Global transcriptomic analysis of an engineered Escherichia coli strain lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system during shikimic acid production in rich culture medium.

BACKGROUND: Efficient production of SA in Escherichia coli has been achieved by modifying key genes of the central carbon metabolism and SA pathway, resulting in overproducing strains grown in batch- or fed-batch-fermentor cultures using a complex broth including glucose and YE. In this study, we performed a GTA to identify those genes significantly upregulated in an engineered E. coli strain, PB12.SA22, in mid EXP (5 h), early STA (STA1, 9 h), and late STA (STA2, 44 h) phases, grown in complex fermentation broth in batch-fermentor cultures. RESULTS: Growth of E. coli PB12.SA22 in complex fermentation broth for SA production resulted in an EXP growth during the first 9 h of cultivation depending of supernatant available aromatic amino acids provided by YE because, when tryptophan was totally consumed, cells entered into a second, low-growth phase (even in the presence of glucose) until 26 h of cultivation. At this point, glucose was completely consumed but SA production continued until the end of the fermentation (50 h) achieving the highest accumulation (7.63 g/L of SA). GTA between EXP/STA1, EXP/STA2 and STA1/STA2 comparisons showed no significant differences in the regulation of genes encoding enzymes of central carbon metabolism as in SA pathway, but those genes encoding enzymes involved in sugar, amino acid, nucleotide/nucleoside, iron and sulfur transport; amino acid catabolism and biosynthesis; nucleotide/nucleoside salvage; acid stress response and modification of IM and OM were upregulated between comparisons. CONCLUSIONS: GTA during SA production in batch-fermentor cultures of strain PB12.SA22 grown in complex fermentation broth during the EXP, STA1 and STA2 phases was studied. Significantly, upregulated genes during the EXP and STA1 phases were associated with transport, amino acid catabolism, biosynthesis, and nucleotide/nucleoside salvage. In STA2, upregulation of genes encoding transporters and enzymes involved in the synthesis and catabolism of Arg suggests that this amino acid could have a key role in the fuelling of carbon toward SA synthesis, whereas upregulation of genes involved in pH stress response, such as membrane modifications, suggests a possible response to environmental conditions imposed on the cell at the end of the fermentation.

Glycosyltransferases Expression Changes in Leuconostoc mesenteroides subsp. mesenteroides ATCC 8293 Grown on Different Carbon Sources.

Leuconostoc mesenteroides strains are common contributors in fermented foods producing a wide variety of polysaccharides from sucrose through glycosyltransferases (GTFs). These polymers have been proposed as protective barriers against acidity, dehydration, heat, and oxidative stress. Despite its presence in many traditional fermented products and their association with food functional properties, regulation of GTFs expression in Ln. mesenteroides is still poorly understood. The strain Ln. mesenteroides ATCC 8293 contains three glucansucrases genes not found in operons, and three fructansucrases genes arranged in two operons, levLX and levC-scrB, a Glycoside-hydrolase. We described the first differential gene expression analysis of this strain when cultivated in different carbon sources. We observed that while GTFs are expressed in the presence of most sugars, they are down-regulated in xylose. We ruled out the regulatory effect of CcpA over GTFs and did not find regulatory elements with a direct effect on glucansucrases in the condition assayed. Our findings suggest that only operon levLX is repressed in xylose by LexA and that both fructansucrases operons can be regulated by the VicK/VicR system and PerR. It is essential to further explore the effect of environmental conditions in Ln. mesenteroides bacteria to better understand GTFs regulation and polymer function.

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

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

Omicron-BA.1 Dispersion Rates in Mexico Varied According to the Regional Epidemic Patterns and the Diversity of Local Delta Subvariants.

PURPOSE: The Omicron subvariant BA.1 of SARS-CoV-2 was first detected in November 2021 and quickly spread worldwide, displacing the Delta variant. In this work, a characterization of the spread of this variant in Mexico is presented. METHODS: The time to fixation of BA.1, the diversity of Delta sublineages, the population density, and the level of virus circulation during the inter-wave interval were determined to analyze differences in BA.1 spread. RESULTS: BA.1 began spreading during the first week of December 2021 and became dominant in the next three weeks, causing the fourth COVID-19 epidemiological surge in Mexico. Unlike previous variants, BA.1 did not exhibit a geographically distinct circulation pattern. However, a regional difference in the speed of the replacement of the Delta variant was observed. CONCLUSIONS: Viral diversity and the relative abundance of the virus in a particular area around the time of the introduction of a new lineage seem to have influenced the spread dynamics, in addition to population density. Nonetheless, if there is a significant difference in the fitness of the variants, or if the time allowed for the competition is sufficiently long, it seems the fitter virus will eventually become dominant, as observed in the eventual dominance of the BA.1.x variant in Mexico.

Two-year follow-up of the COVID-19 pandemic in Mexico.

Background: After the initial outbreak in China (December 2019), the World Health Organization declared COVID-19 a pandemic on March 11th, 2020. This paper aims to describe the first 2 years of the pandemic in Mexico. Design and methods: This is a population-based longitudinal study. We analyzed data from the national COVID-19 registry to describe the evolution of the pandemic in terms of the number of confirmed cases, hospitalizations, deaths and reported symptoms in relation to health policies and circulating variants. We also carried out logistic regression to investigate the major risk factors for disease severity. Results: From March 2020 to March 2022, the coronavirus disease 2019 (COVID-19) pandemic in Mexico underwent four epidemic waves. Out of 5,702,143 confirmed cases, 680,063 were hospitalized (11.9%), and 324,436 (5.7%) died. Even if there was no difference in susceptibility by gender, males had a higher risk of death (CFP: 7.3 vs. 4.2%) and hospital admission risk (HP: 14.4 vs. 9.5%). Severity increased with age. With respect to younger ages (0-17 years), the 60+ years or older group reached adjusted odds ratios of 9.63 in the case of admission and 53.05 (95% CI: 27.94-118.62) in the case of death. The presence of any comorbidity more than doubled the odds ratio, with hypertension-diabetes as the riskiest combination. While the wave peaks increased over time, the odds ratios for developing severe disease (waves 2, 3, and 4 to wave 1) decreased to 0.15 (95% CI: 0.12-0.18) in the fourth wave. Conclusion: The health policy promoted by the Mexican government decreased hospitalizations and deaths, particularly among older adults with the highest risk of admission and death. Comorbidities augment the risk of developing severe illness, which is shown to rise by double in the Mexican population, particularly for those reported with hypertension-diabetes. Factors such as the decrease in the severity of the SARS-CoV2 variants, changes in symptomatology, and advances in the management of patients, vaccination, and treatments influenced the decrease in mortality and hospitalizations.

Dominance of Three Sublineages of the SARS-CoV-2 Delta Variant in Mexico.

In this study, we analyzed the sequences of SARS-CoV-2 isolates of the Delta variant in Mexico, which has completely replaced other previously circulating variants in the country due to its transmission advantage. Among all the Delta sublineages that were detected, 81.5 % were classified as AY.20, AY.26, and AY.100. According to publicly available data, these only reached a world prevalence of less than 1%, suggesting a possible Mexican origin. The signature mutations of these sublineages are described herein, and phylogenetic analyses and haplotype networks are used to track their spread across the country. Other frequently detected sublineages include AY.3, AY.62, AY.103, and AY.113. Over time, the main sublineages showed different geographical distributions, with AY.20 predominant in Central Mexico, AY.26 in the North, and AY.100 in the Northwest and South/Southeast. This work describes the circulation, from May to November 2021, of the primary sublineages of the Delta variant associated with the third wave of the COVID-19 pandemic in Mexico and highlights the importance of SARS-CoV-2 genomic surveillance for the timely identification of emerging variants that may impact public health.

The Alpha Variant (B.1.1.7) of SARS-CoV-2 Failed to Become Dominant in Mexico.

During the coronavirus disease 2019 (COVID-19) pandemic, the emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in the United Kingdom in September 2020, was well documented in different areas of the world and became a global public health concern because of its increased transmissibility. The B.1.1.7 lineage was first detected in Mexico during December 2020, showing a slow progressive increase in its circulation frequency, which reached its maximum in May 2021 but never became predominant. In this work, we analyzed the patterns of diversity and distribution of this lineage in Mexico using phylogenetic and haplotype network analyses. Despite the reported increase in transmissibility of the B.1.1.7 lineage, in most Mexican states, it did not displace cocirculating lineages, such as B.1.1.519, which dominated the country from February to May 2021. Our results show that the states with the highest prevalence of B.1.1.7 were those at the Mexico-U.S. border. An apparent pattern of dispersion of this lineage from the northern states of Mexico toward the center or the southeast was observed in the largest transmission chains, indicating possible independent introduction events from the United States. However, other entry points cannot be excluded, as shown by multiple introduction events. Local transmission led to a few successful haplotypes with a localized distribution and specific mutations indicating sustained community transmission. IMPORTANCE The emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) throughout the world were due to its increased transmissibility. However, it did not displace cocirculating lineages in most of Mexico, particularly B.1.1.519, which dominated the country from February to May 2021. In this work, we analyzed the distribution of B.1.1.7 in Mexico using phylogenetic and haplotype network analyses. Our results show that the states with the highest prevalence of B.1.1.7 (around 30%) were those at the Mexico-U.S. border, which also exhibited the highest lineage diversity, indicating possible introduction events from the United States. Also, several haplotypes were identified with a localized distribution and specific mutations, indicating that sustained community transmission occurred in the country.

Transcriptional profiling of fetal hypothalamic TRH neurons.

BACKGROUND: During murine hypothalamic development, different neuroendocrine cell phenotypes are generated in overlapping periods; this suggests that cell-type specific developmental programs operate to achieve complete maturation. A balance between programs that include cell proliferation, cell cycle withdrawal as well as epigenetic regulation of gene expression characterizes neurogenesis. Thyrotropin releasing hormone (TRH) is a peptide that regulates energy homeostasis and autonomic responses. To better understand the molecular mechanisms underlying TRH neuron development, we performed a genome wide study of its transcriptome during fetal hypothalamic development. RESULTS: In primary cultures, TRH cells constitute 2% of the total fetal hypothalamic cell population. To purify these cells, we took advantage of the fact that the segment spanning -774 to +84 bp of the Trh gene regulatory region confers specific expression of the green fluorescent protein (GFP) in the TRH cells. Transfected TRH cells were purified by fluorescence activated cell sorting, various cell preparations pooled, and their transcriptome compared to that of GFP- hypothalamic cells. TRH cells undergoing the terminal phase of differentiation, expressed genes implicated in protein biosynthesis, intracellular signaling and transcriptional control. Among the transcription-associated transcripts, we identified the transcription factors Klf4, Klf10 and Atf3, which were previously uncharacterized within the hypothalamus. CONCLUSION: To our knowledge, this is one of the first reports identifying transcripts with a potentially important role during the development of a specific hypothalamic neuronal phenotype. This genome-scale study forms a rational foundation for identifying genes that might participate in the development and function of hypothalamic TRH neurons.

Distribution and preservation of the components of the engulfment. What is beyond representative genomes?

Engulfment requires the coordinated, targeted synthesis and degradation of peptidoglycan at the leading edge of the engulfing membrane to allow the mother cell to completely engulf the forespore. Proteins such as the DMP and Q:AH complexes in Bacillus subtilis are essential for engulfment, as are a set of accessory proteins including GerM and SpoIIB, among others. Experimental and bioinformatic studies of these proteins in bacteria distinct from Bacillus subtilis indicate that fundamental differences exist regarding the organization and mechanisms used to successfully perform engulfment. As a consequence, the distribution and prevalence of the proteins involved in engulfment and other proteins that participate in different sporulation stages have been studied using bioinformatic approaches. These works are based on the prediction of orthologs in the genomes of representative Firmicutes and have been helpful in tracing hypotheses about the origin and evolution of sporulation genes, some of which have been postulated as sporulation signatures. To date, an extensive study of these signatures outside of the representative Firmicutes is not available. Here, we asked whether phyletic profiles of proteins involved in engulfment can be used as signatures able to describe the sporulation phenotype. We tested this hypothesis in a set of 954 Firmicutes, finding preserved phyletic profiles defining signatures at the genus level. Finally, a phylogenetic reconstruction based on non-redundant phyletic profiles at the family level shows the non-monophyletic origin of these proteins due to gain/loss events along the phylum Firmicutes.

Definition of the Metagenomic Profile of Ocean Water Samples From the Gulf of Mexico Based on Comparison With Reference Samples From Sites Worldwide.

Computational and statistical analysis of shotgun metagenomes can predict gene abundance and is helpful for elucidating the functional and taxonomic compositions of environmental samples. Gene products are compared against physicochemical conditions or perturbations to shed light on the functions performed by the microbial community of an environmental sample; however, this information is not always available. The present study proposes a method for inferring the metabolic potential of metagenome samples by constructing a reference based on determining the probability distribution of the counts of each enzyme annotated. To test the methodology, we used marine water samples distributed worldwide as references. Then, the references were utilized to compare the annotated enzymes of two different water samples extracted from the Gulf of Mexico (GoM) to distinguish those enzymes with atypical behavior. The enzymes whose annotation counts presented frequencies significantly different from those of the reference were used to perform metabolic reconstruction, which naturally identified pathways. We found that several of the enzymes were involved in the biodegradation of petroleum, which is consistent with the impact of human hydrocarbon extraction activity and its ubiquitous presence in the GoM. The examination of other reconstructed pathways revealed significant enzymes indicating the presence of microbial communities characterizing each ocean depth and ocean cycle, providing a fingerprint of each sampled site.

Characterization of Enterobacter cloacae BAGM01 Producing a Thermostable and Alkaline-Tolerant Rhamnolipid Biosurfactant from the Gulf of Mexico.

The search for novel biosurfactants (Bs) requires the isolation of microorganisms from different environments. The Gulf of Mexico (GoM) is a geographical area active in the exploration and exploitation of hydrocarbons. Recent metagenomic and microbiologic studies in this area suggested a potential richness for novel Bs microbial producers. In this work, nineteen bacterial consortia from the GoM were isolated at different depths of the water column and marine sediments. Bs production from four bacterial consortia was detected by the CTAB test and their capacity to reduce surface tension (ST), emulsion index (EI24), and hemolytic activity. These bacterial consortia produced Bs in media supplemented with kerosene, diesel, or sucrose. Cultivable bacteria from these consortia were isolated and identified by bacterial polyphasic characterization. In some consortia, Enterobacter cloacae was the predominant specie. E. cloacae BAGM01 presented Bs activity in minimal medium and was selected to improve its Bs production using a Taguchi and Box-Behnken experimental design; this strain was able to grow and presented Bs activity at 35 g L-1 of NaCl. This Bs decreased ST to around 34.5 ± 0.56 mNm-1 and presented an EI24 of 71 ± 1.27%. Other properties of this Bs were thermal stability, stability in alkaline conditions, and stability at high salinity, conferring important and desirable characteristics in multiple industries. The analysis of the genome of E. cloacae BAGM01 showed the presence of rhlAB genes that have been reported in the synthesis of rhamnolipids, and alkAB genes that are related to the degradation of alkanes. The bioactive molecule was identified as a rhamnolipid after HPLC derivatization, 1H NMR, and UPLC-QTOF-MS analysis.

New insights into the regulatory networks of paralogous genes in bacteria.

Extensive genomic studies on gene duplication in model organisms such as Escherichia coli and Saccharomyces cerevisiae have recently been undertaken. In these models, it is commonly considered that a duplication event may include a transcription factor (TF), a target gene, or both. Following a gene duplication episode, varying scenarios have been postulated to describe the evolution of the regulatory network. However, in most of these, the TFs have emerged as the most important and in some cases the only factor shaping the regulatory network as the organism responds to a natural selection process, in order to fulfil its metabolic needs. Recent findings concerning the regulatory role played by elements other than TFs have indicated the need to reassess these early models. Thus, we performed an exhaustive review of paralogous gene regulation in E. coli and Bacillus subtilis based on published information, available in the NCBI PubMed database and in well-established regulatory databases. Our survey reinforces the notion that despite TFs being the most prominent components shaping the regulatory networks, other elements are also important. These include small RNAs, riboswitches, RNA-binding proteins, sigma factors, protein-protein interactions and DNA supercoiling, which modulate the expression of genes involved in particular metabolic processes or induce a more complex response in terms of the regulatory networks of paralogous genes in an integrated interplay with TFs.

Structure and Evolution of Acinetobacter baumannii Plasmids.

Acinetobacter baumannii is an emergent bacterial pathogen that provokes many types of infections in hospitals around the world. The genome of this organism consists of a chromosome and plasmids. These plasmids vary over a wide size range and many of them have been linked to the acquisition of antibiotic-resistance genes. Our bioinformatic analyses indicate that A. baumannii plasmids belong to a small number of plasmid lineages. The general structure of these lineages seems to be very stable and consists not only of genes involved in plasmid maintenance functions but of gene sets encoding poorly characterized proteins, not obviously linked to survival in the hospital setting, and opening the possibility that they improve the parasitic properties of plasmids. An analysis of genes involved in replication, suggests that members of the same plasmid lineage are part of the same plasmid incompatibility group. The same analysis showed the necessity of classifying the Rep proteins in ten new groups, under the scheme proposed by Bertini et al. (2010). Also, we show that some plasmid lineages have the potential capacity to replicate in many bacterial genera including those embracing human pathogen species, while others seem to replicate only within the limits of the Acinetobacter genus. Moreover, some plasmid lineages are widely distributed along the A. baumannii phylogenetic tree. Despite this, a number of them lack genes involved in conjugation or mobilization functions. Interestingly, only 34.6% of the plasmids analyzed here possess antibiotic resistance genes and most of them belong to fourteen plasmid lineages of the twenty one described here. Gene flux between plasmid lineages appears primarily limited to transposable elements, which sometimes carry antibiotic resistance genes. In most plasmid lineages transposable elements and antibiotic resistance genes are secondary acquisitions. Finally, broad host-range plasmids appear to have played a crucial role.

Author Correction: Analysis of sequencing strategies and tools for taxonomic annotation: Defining standards for progressive metagenomics.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Metagenomic Profiling and Microbial Metabolic Potential of Perdido Fold Belt (NW) and Campeche Knolls (SE) in the Gulf of Mexico.

The Gulf of Mexico (GoM) is a particular environment that is continuously exposed to hydrocarbon compounds that may influence the microbial community composition. We carried out a metagenomic assessment of the bacterial community to get an overall view of this geographical zone. We analyzed both taxonomic and metabolic markers profiles to explain how the indigenous GoM microorganims participate in the biogeochemical cycling. Two geographically distant regions in the GoM, one in the north-west (NW) and one in the south-east (SE) of the GoM were analyzed and showed differences in their microbial composition and metabolic potential. These differences provide evidence the delicate equilibrium that sustains microbial communities and biogeochemical cycles. Based on the taxonomy and gene groups, the NW are more oxic sediments than SE ones, which have anaerobic conditions. Both water and sediments show the expected sulfur, nitrogen, and hydrocarbon metabolism genes, with particularly high diversity of the hydrocarbon-degrading ones. Accordingly, many of the assigned genera were associated with hydrocarbon degradation processes, Nitrospira and Sva0081 were the most abundant in sediments, while Vibrio, Alteromonas, and Alcanivorax were mostly detected in water samples. This basal-state analysis presents the GoM as a potential source of aerobic and anaerobic hydrocarbon degradation genes important for the ecological dynamics of hydrocarbons and the potential use for water and sediment bioremediation processes.

Identification of network topological units coordinating the global expression response to glucose in Bacillus subtilis and its comparison to Escherichia coli.

BACKGROUND: Glucose is the preferred carbon and energy source for Bacillus subtilis and Escherichia coli. A complex regulatory network coordinates gene expression, transport and enzymatic activities, in response to the presence of this sugar. We present a comparison of the cellular response to glucose in these two model organisms, using an approach combining global transcriptome and regulatory network analyses. RESULTS: Transcriptome data from strains grown in Luria-Bertani medium (LB) or LB+glucose (LB+G) were analyzed, in order to identify differentially transcribed genes in B. subtilis. We detected 503 genes in B. subtilis that change their relative transcript levels in the presence of glucose. A similar previous study identified 380 genes in E. coli, which respond to glucose. Catabolic repression was detected in the case of transport and metabolic interconversion activities for both bacteria in LB+G. We detected an increased capacity for de novo synthesis of nucleotides, amino acids and proteins. A comparison between orthologous genes revealed that global regulatory functions such as transcription, translation, replication and genes relating to the central carbon metabolism, presented similar changes in their levels of expression. An analysis of the regulatory network of a subset of genes in both organisms revealed that the set of regulatory proteins responsible for similar physiological responses observed in the transcriptome analysis are not orthologous. An example of this observation is that of transcription factors mediating catabolic repression for most of the genes that displayed reduced transcript levels in the case of both organisms. In terms of topological functional units in both these bacteria, we found interconnected modules that cluster together genes relating to heat shock, respiratory functions, carbon and peroxide metabolism. Interestingly, B. subtilis functions not found in E. coli, such as sporulation and competence were shown to be interconnected, forming modules subject to catabolic repression at the level of transcription. CONCLUSION: Our results demonstrate that the response to glucose is partially conserved in model organisms E. coli and B. subtilis, including genes encoding basic functions such as transcription, translation, replication and genes involved in the central carbon metabolism.

Phenotypic and genomic analysis of Zymomonas mobilis ZM4 mutants with enhanced ethanol tolerance.

Zymomonas mobilis ZM4 is an ethanol-producing microbe that is constitutively tolerant to this solvent. For a better understanding of the ethanol tolerance phenomenon we obtained and characterized two ZM4 mutants (ER79ap and ER79ag) with higher ethanol tolerance than the wild-type. Mutants were evaluated in different ethanol concentrations and this analysis showed that mutant ER79ap was more tolerant and had a better performance in terms of cell viability, than the wild-type strain and ER79ag mutant. Genotyping of the mutant strains showed that both carry non-synonymous mutations in clpP and spoT/relA genes. A third non-synonymous mutation was found only in strain ER79ap, in the clpB gene. Considering that ER79ap has the best tolerance to added ethanol, the mutant alleles of this strain were evaluated in ZM4 and here we show that while all of them contribute to ethanol tolerance, mutation within spoT/relA gene seems to be the most important.

Prediction of protein architectures involved in the signaling-pathway initiating sporulation in Firmicutes.

OBJECTIVES: Like many other proteins, those belonging to the signal transduction cascade initiating sporulation (Spo0 pathway) have conserved protein domains (Capra and Laub in Annu Rev Microbiol 66:325-47, 2012). Improvements in bioinformatics applications to discover proteins involved in the initiation of the sporulating cascade in newly sequenced genomes is an important task that requires rigorous comparative genomic methods and manual curation to identify endospore-forming bacteria. This note aims to present a collection of predicted proteins involved in the Spo0 pathway found in the proteomes of fully sequenced and manually curated endospore-forming Firmicutes species. This collection may serve as a guide to conduct future experiments in endospore formers in genomic and metagenomic projects. DATA DESCRIPTION: Similar to the report of Davidson et al. (PLoS Genet 14:1-33, 2018), we used Pfam profiles (El-Gebali et al. in Nucleic Acids Res 47:D427-32, 2019) defining each protein and the genomic context surrounding the query gene to predict probable orthologs of the Spo0 pathway in Firmicutes. We present in this note a collection of 325 Firmicutes species organized by phylogenetic class and classified as spore formers, non-spore formers or unknown spore phenotype based on published literature, for which we predicted probable orthologs defining the signal transduction pathway initiating sporulation.