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.

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.

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 11^th, 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.

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.

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

•

Z. mobilis ER79ag and ER79ap ethanol mutants were obtained by adaptive evolution.

•

ER79ap had a better cell viability than the WT and ER79ap under ethanol stress.

•

Mutants shared SNVs in clpP and spoT/relA, in addition ER79ap has a SNP in clpB.

•

Mutant allele spoT/relA of ER79ap seems to be more important to ethanol tolerance.

•

Glucose consumption and ethanol production were not affected in mutant strains.

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.

The Salmonella enterica serovar Typhi ltrR-ompR-ompC-ompF genes are involved in resistance to the bile salt sodium deoxycholate and in bacterial transformation

A characterization of the LtrR regulator, an S. Typhi protein belonging to the LysR family is presented. Proteomics, outer membrane protein profiles and transcriptional analyses demonstrated that LtrR is required for the synthesis of OmpR, OmpC and OmpF. DNA-protein interaction analysis showed that LtrR binds to the regulatory region of ompR and then OmpR interacts with the ompC and ompF promoters inducing porin synthesis. LtrR-dependent and independent ompR promoters were identified, and both promoters are involved in the synthesis of OmpR for OmpC and OmpF production. To define the functional role of the ltrR-ompR-ompC-ompF genetic network, mutants in each gene were obtained. We found that ltrR, ompR, ompC and ompF were involved in the control of bacterial transformation, while the two regulators and ompC are necessary for the optimal growth of S. Typhi in the presence of one of the major bile salts found in the gut, sodium deoxycholate. The data presented establish the pivotal role of LtrR in the regulatory network of porin synthesis and reveal new genetic strategies of survival and cellular adaptation to the environment used by Salmonella.

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.

Prokaryotic regulatory systems biology: Common principles governing the functional architectures of Bacillus subtilis and Escherichia coli unveiled by the natural decomposition approach

Escherichia coli and Bacillus subtilis are two of the best-studied prokaryotic model organisms. Previous analyses of their transcriptional regulatory networks have shown that they exhibit high plasticity during evolution and suggested that both converge to scale-free-like structures. Nevertheless, beyond this suggestion, no analyses have been carried out to identify the common systems-level components and principles governing these organisms. Here we show that these two phylogenetically distant organisms follow a set of common novel biologically consistent systems principles revealed by the mathematically and biologically founded natural decomposition approach. The discovered common functional architecture is a diamond-shaped, matryoshka-like, three-layer (coordination, processing, and integration) hierarchy exhibiting feedback, which is shaped by four systems-level components: global transcription factors (global TFs), locally autonomous modules, basal machinery and intermodular genes. The first mathematical criterion to identify global TFs, the κ-value, was reassessed on B. subtilis and confirmed its high predictive power by identifying all the previously reported, plus three potential, master regulators and eight sigma factors. The functionally conserved cores of modules, basal cell machinery, and a set of non-orthologous common physiological global responses were identified via both orthologous genes and non-orthologous conserved functions. This study reveals novel common systems principles maintained between two phylogenetically distant organisms and provides a comparison of their lifestyle adaptations. Our results shed new light on the systems-level principles and the fundamental functions required by bacteria to sustain life.

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.

Regulatory network of Escherichia coli: consistency between literature knowledge and microarray profiles

The transcriptional network of Escherichia coli may well be the most complete experimentally characterized network of a single cell. A rule-based approach was built to assess the degree of consistency between whole-genome microarray experiments in different experimental conditions and the accumulated knowledge in the literature compiled in RegulonDB, a data base of transcriptional regulation and operon organization in E. coli. We observed a high and statistical significant level of consistency, ranging from 70%-87%. When effector metabolites of regulatory proteins are not considered in the prediction of the active or inactive state of the regulators, consistency falls by up to 40%. Similarly, consistency decreases when rules for multiple regulatory interactions are altered or when "on" and "off" entries were assigned randomly. We modified the initial state of regulators and evaluated the propagation of errors in the network that do not correlate linearly with the connectivity of regulators. We interpret this deviation mainly as a result of the existence of redundant regulatory interactions. Consistency evaluation opens a new space of dialogue between theory and experiment, as the consequences of different assumptions can be evaluated and compared.

Environmental conditions and transcriptional regulation inEscherichia coli: a physiological integrative approach

Bacteria develop a number of devices for sensing, responding, and adapting to different environmental conditions. Understanding within a genomic perspective how the transcriptional machinery of bacteria is modulated, as a response for changing conditions, is a major challenge for biologists. Knowledge of which genes are turned on or turned off under specific conditions is essential for our understanding of cell behavior. In this study we describe how the information pertaining to gene expression and associated growth conditions (even with very little knowledge of the associated regulatory mechanisms) is gathered from the literature and incorporated into RegulonDB, a database on transcriptional regulation and operon organization in E. coli. The link between growth conditions, signal transduction, and transcriptional regulation is modeled in the database in a simple format that highlights biological relevant information. As far as we know, there is no other database that explicitly clarifies the effect of environmental conditions on gene transcription. We discuss how this knowledge constitutes a benchmark that will impact future research aimed at integration of regulatory responses in the cell; for instance, analysis of microarrays, predicting culture behavior in biotechnological processes, and comprehension of dynamics of regulatory networks. This integrated knowledge will contribute to the future goal of modeling the behavior of E. coli as an entire cell. The RegulonDB database can be accessed on the web at the URL: http://www.cifn.unam.mx/Computational_Biology/regulondb/.