RESUMO

The evolution, survival, and adaptation of microbes are consequences of gene duplication, acquisition, and divergence in response to environmental challenges. In this context, enzymes play a central role in the evolution of organisms, because they are fundamental in cell metabolism. Here, we analyzed the enzymatic repertoire in 6,467 microbial genomes, including their abundances, and their associations with metabolic maps. We found that the enzymes follow a power-law distribution, in relation to the genome sizes. Therefore, we evaluated the total proportion enzymatic classes in relation to the genomes, identifying a descending-order proportion: transferases (EC:2.-), hydrolases (EC:3.-), oxidoreductases (EC:1.-), ligases (EC:6.-), lyases (EC:4.-), isomerases (EC:5.-), and translocases (EC:7-.). In addition, we identified a preferential use of enzymatic classes in metabolism pathways for xenobiotics, cofactors and vitamins, carbohydrates, amino acids, glycans, and energy. Therefore, this analysis provides clues about the functional constraints associated with the enzymatic repertoire of functions in Bacteria and Archaea.

RESUMO

Breast cancer is the most common malignancy in women around the world. Intratumor and intertumoral heterogeneity persist in mammary tumors. Therefore, the identification of biomarkers is essential for the treatment of this malignancy. This study analyzed 28,143 genes expressed in 49 breast cancer cell lines using a Weighted Gene Co-expression Network Analysis to determine specific target proteins for Basal A, Basal B, Luminal A, Luminal B, and HER2 ampl breast cancer subtypes. Sixty-five modules were identified, of which five were characterized as having a high correlation with breast cancer subtypes. Genes overexpressed in the tumor were found to participate in the following mechanisms: regulation of the apoptotic process, transcriptional regulation, angiogenesis, signaling, and cellular survival. In particular, we identified the following genes, considered as hubs: IFIT3, an inhibitor of viral and cellular processes; ETS1, a transcription factor involved in cell death and tumorigenesis; ENSG00000259723 lncRNA, expressed in cancers; AL033519.3, a hypothetical gene; and TMEM86A, important for regulating keratinocyte membrane properties, considered as a key in Basal A, Basal B, Luminal A, Luminal B, and HER2 ampl breast cancer subtypes, respectively. The modules and genes identified in this work can be used to identify possible biomarkers or therapeutic targets in different breast cancer subtypes.

Assuntos

RESUMO

Deep learning models (DLMs) have gained importance in predicting, detecting, translating, and classifying a diversity of inputs. In bioinformatics, DLMs have been used to predict protein structures, transcription factor-binding sites, and promoters. In this work, we propose a hybrid model to identify transcription factors (TFs) among prokaryotic and eukaryotic protein sequences, named Deep Regulation (DeepReg) model. Two architectures were used in the DL model: a convolutional neural network (CNN), and a bidirectional long-short-term memory (BiLSTM). DeepReg reached a precision of 0.99, a recall of 0.97, and an F1-score of 0.98. The quality of our predictions, the bias-variance trade-off approach, and the characterization of new TF predictions were evaluated and compared against those produced by DeepTFactor, as well as against experimental data from three model organisms. Predictions based on our DLM tended to exhibit less variance and bias than those from DeepTFactor, thus increasing reliability and decreasing overfitting.

Assuntos

RESUMO

In this work we carried out an in silico analysis to understand the interaction between InvF-SicA and RNAP in the bacterium Salmonella Typhimurium strain LT2. Structural analysis of InvF allowed the identification of three possible potential cavities for interaction with SicA. This interaction could occur with the structural motif known as tetratricopeptide repeat (TPR) 1 and 2 in the two cavities located in the interface of the InvF and α-CTD of RNAP. Indeed, molecular dynamics simulations showed that SicA stabilizes the Helix-turn-Helix DNA-binding motifs, i.e., maintaining their proper conformation, mainly in the DNA Binding Domain (DBD). Finally, to evaluate the role of amino acids that contribute to protein-protein affinity, an alanine scanning mutagenesis approach, indicated that R177 and R181, located in the DBD motif, caused the greatest changes in binding affinity with α-CTD, suggesting a central role in the stabilization of the complex. However, it seems that the N-terminal region also plays a key role in the protein-protein interaction, especially the amino acid R40, since we observed conformational flexibility in this region allowing it to interact with interface residues. We consider that this analysis opens the possibility to validate experimentally the amino acids involved in protein-protein interactions and explore other regulatory complexes where chaperones are involved.

Assuntos

RESUMO

The combination of signals from the T-cell receptor (TCR) and co-stimulatory molecules triggers transcriptional programs that lead to proliferation, cytokine secretion, and effector functions. We compared the impact of engaging the TCR with CD28 and/or CD43 at different time points relative to TCR engagement on T-cell function. TCR and CD43 simultaneous engagement resulted in higher CD69 and PD-1 expression levels than in TCR and CD28-stimulated cells, with a cytokine signature of mostly effector, inflammatory, and regulatory cytokines, while TCR and CD28-activated cells secreted all categories of cytokines, including stimulatory cytokines. Furthermore, the timing of CD43 engagement relative to TCR ligation, and to a lesser degree that of CD28, resulted in distinct patterns of expression of cytokines, chemokines, and growth factors. Complete cell activation was observed when CD28 or CD43 were engaged simultaneously with or before the TCR, but ligating the TCR before CD43 or CD28 failed to complete a cell activation program regarding cytokine secretion. As the order in which CD43 or CD28 and the TCR were engaged resulted in different combinations of cytokines that shape distinct T-cell immune programs, we analyzed their upstream sequences to assess whether the combinations of cytokines were associated with different sets of regulatory elements. We found that the order in which the TCR and CD28 or CD43 are engaged predicts the recruitment of specific sets of chromatin remodelers and TFSS, which ultimately regulate T-cell polarization and plasticity. Our data underscore that the combination of co-stimulatory molecules and the time when they are engaged relative to the TCR can change the cell differentiation program.

Assuntos

RESUMO

Salmonella enterica serovar Typhimurium causes gastroenteritis and systemic infections in humans. For this bacterium the expression of a type III secretion system (T3SS) and effector proteins encoded in the Salmonella pathogenicity island-1 (SPI-1), is keystone for the virulence of this bacterium. Expression of these is controlled by a regulatory cascade starting with the transcriptional regulators HilD, HilC and RtsA that induce the expression of HilA, which then activates expression of the regulator InvF, a transcriptional regulator of the AraC/XylS family. InvF needs to interact with the chaperone SicA to activate transcription of SPI-1 genes including sicA, sopB, sptP, sopE, sopE2, and STM1239. InvF very likely acts as a classical activator; however, whether InvF interacts with the RNA polymerase alpha subunit RpoA has not been determined. Results from this study confirm the interaction between InvF with SicA and reveal that both proteins interact with the RNAP alpha subunit. Thus, our study further supports that the InvF/SicA complex acts as a classical activator. Additionally, we showed for the first time an interaction between a chaperone of T3SS effectors (SicA) and the RNAP.

Assuntos

RESUMO

Entamoeba histolytica (E. histolytica) exhibits a remarkable capacity to respond to thermal shock stress through a sophisticated genetic regulation mechanism. This process is carried out via Heat Shock Response Elements (HSEs), which are recognized by Heat Shock Transcription Factors (EhHSTFs), enabling fine and precise control of gene expression. Our study focused on screening for HSEs in the promoters of the E. histolytica genome, specifically analyzing six HSEs, including Ehpgp5, EhrabB1, EhrabB4, EhrabB5, Ehmlbp, and Ehhsp100. We discovered 2578 HSEs, with 1412 in promoters of hypothetical genes and 1166 in coding genes. We observed that a single promoter could contain anywhere from one to five HSEs. Gene ontology analysis revealed the presence of HSEs in essential genes for the amoeba, including cysteine proteinases, ribosomal genes, Myb family DNA-binding proteins, and Rab GTPases, among others. Complementarily, our molecular docking analyses indicate that these HSEs are potentially recognized by EhHSTF5, EhHSTF6, and EhHSTF7 factors in their trimeric conformation. These findings suggest that E. histolytica has the capability to regulate a wide range of critical genes via HSE-EhHSTFs, not only for thermal stress response but also for vital functions of the parasite. This is the first comprehensive study of HSEs in the genome of E. histolytica, significantly contributing to the understanding of its genetic regulation and highlighting the complexity and precision of this mechanism in the parasite's survival.

Assuntos

RESUMO

An easy and fast strategy to compare functionally the metabolic maps is described. The KEGG metabolic maps are transformed into linear Enzymatic Step Sequences (ESS) using the Breadth First Search (BFS) algorithm. To do this, the KGML files are retrieved, and directed graph representations are created; where the nodes represent enzymes or enzymatic complexes, and the edges represent a compound, that is the 'product' from one reaction and a 'substrate' for the next. Then, a set of initialization nodes are selected, and used as the root for the construction of the BFS tree. This tree is used as a guide to the construction of the ESS. From each leaf (terminal node), the path is traced backwards until it reaches the root metabolic map and with two or fewer neighbors in the graph. In a second step, the ESS are compared with a Dynamic Programing algorithm, considering an "ad hoc" substitution matrix, and minimizing the global score. The dissimilarity values between two EC numbers ranged from 0 to 1, where 0 indicates similar EC numbers, and 1 indicates different EC numbers. Finally, the alignment is evaluated by using the normalized entropy-based function, considering a threshold of ≤ 0.27 as significant.•The KEGG metabolic maps are transformed into linear Enzymatic Step Sequences (ESS) using the Breadth First Search (BFS) algorithm.•Nodes represent enzymes or enzymatic complexes, and the edges represent a compound, that is 'product' from one reaction and a 'substrate' for the next.•The ESS are compared with a Dynamic Programing algorithm, considering an "ad hoc" substitution matrix, and minimizing the global score.

RESUMO

Biological systems respond to environmental perturbations and a large diversity of compounds through gene interactions, and these genetic factors comprise complex networks. Experimental information from transcriptomic studies has allowed the identification of gene networks that contribute to our understanding of microbial adaptations. In this study, we analyzed the gene co-expression networks of three Bifidobacterium species in response to different types of human milk oligosaccharides (HMO) using weighted gene co-expression analysis (WGCNA). RNA-seq data obtained from Geo Datasets were obtained for Bifidobacterium longum subsp. Infantis, Bifidobacterium bifidum and Bifidobacterium longum subsp. Longum. Between 10 and 20 co-expressing modules were obtained for each dataset. HMO-associated genes appeared in the modules with more genes for B. infantis and B. bifidum, in contrast with B. longum. Hub genes were identified in each module, and in general they participated in conserved essential processes. Certain modules were differentially enriched with LacI-like transcription factors, and others with certain metabolic pathways such as the biosynthesis of secondary metabolites. The three Bifidobacterium transcriptomes showed distinct regulation patterns for HMO utilization. HMO-associated genes in B. infantis co-expressed in two modules according to their participation in galactose or N-Acetylglucosamine utilization. Instead, B. bifidum showed a less structured co-expression of genes participating in HMO utilization. Finally, this category of genes in B. longum clustered in a small module, indicating a lack of co-expression with main cell processes and suggesting a recent acquisition. This study highlights distinct co-expression architectures in these bifidobacterial genomes during HMO consumption, and contributes to understanding gene regulation and co-expression in these species of the gut microbiome.

RESUMO

Archaea are a vast and unexplored cellular domain that thrive in a high diversity of environments, having central roles in processes mediating global carbon and nutrient fluxes. For these organisms to balance their metabolism, the appropriate regulation of their gene expression is essential. A key momentum in regulating genes responsible for the life maintenance of archaea is when transcription factor proteins bind to the promoter element. This DNA segment is conserved, which enables its exploration by machine learning techniques. Here, we trained and tested a support vector machine with 3935 known archaeal promoter sequences. All promoter sequences were coded into DNA Duplex Stability. After, we performed a model interpretation task to map the decision pattern of the classification procedure. We also used a dataset of known-promoter sequences for validation. Our results showed that an AT rich region around position - 27 upstream (relative to the start TSS) is the most conserved in the analyzed organisms. In addition, we were able to identify the BRE element (- 33), the PPE (at - 10) and a position at + 3, that provides a more understandable picture of how promoters are organized in all the archaeal organisms. Finally, we used the interpreted model to identify potential promoter sequences of 135 unannotated organisms, delivering regulatory regions annotation of archaea in a scale never accomplished before ( https://pcyt.unam.mx/gene-regulation/ ). We consider that this approach will be useful to understand how gene regulation is achieved in other organisms apart from the already established transcription factor binding sites.

Assuntos

RESUMO

Introduction: Biological systems respond to environmental disturbances and a wide range of compounds through complex gene interaction networks. The enormous growth of experimental information obtained using large-scale genomic techniques such as microarrays and RNA sequencing led to the construction of a wide variety of gene co-expression networks in recent years. These networks allow the discovery of clusters of co-expressed genes that potentially work in the same process linking them to biological processes often of interest to industrial, medicinal, and academic research. Methods: In this study, we built the gene co-expression network of Ustilago maydis from the gene expression data of 168 samples belonging to 19 series, which correspond to the GPL3681 platform deposited in the NCBI using WGCNA software. This network was analyzed to identify clusters of co-expressed genes, gene hubs and Gene Ontology terms. Additionally, we identified relevant modules through a hypergeometric approach based on a predicted set of transcription factors and virulence genes. Results and Discussion: We identified 13 modules in the gene co-expression network of U. maydis. The TFs enriched in the modules of interest belong to the superfamilies of Nucleic acid-binding proteins, Winged helix DNA-binding, and Zn2/Cys6 DNA-binding. On the other hand, the modules enriched with virulence genes were classified into diseases related to corn smut, Invasive candidiasis, among others. Finally, a large number of hypothetical, a large number of hypothetical genes were identified as highly co-expressed with virulence genes, making them possible experimental targets.

RESUMO

Gene regulation is a key process for all microorganisms, as it allows them to adapt to different environmental stimuli. However, despite the relevance of gene expression control, for only a handful of organisms is there related information about genome regulation. In this work, we inferred the gene regulatory networks (GRNs) of bacterial and archaeal genomes by comparisons with six organisms with well-known regulatory interactions. The references we used are: Escherichia coli K-12 MG1655, Bacillus subtilis 168, Mycobacterium tuberculosis, Pseudomonas aeruginosa PAO1, Salmonella enterica subsp. enterica serovar typhimurium LT2, and Staphylococcus aureus N315. To this end, the inferences were achieved in two steps. First, the six model organisms were contrasted in an all-vs-all comparison of known interactions based on Transcription Factor (TF)-Target Gene (TG) orthology relationships and Transcription Unit (TU) assignments. In the second step, we used a guilt-by-association approach to infer the GRNs for 12,230 bacterial and 649 archaeal genomes based on TF-TG orthology relationships of the six bacterial models determined in the first step. Finally, we discuss examples to show the most relevant results obtained from these inferences. A web server with all the predicted GRNs is available at https://regulatorynetworks.unam.mx/ or http://132.247.46.6/.

RESUMO

DNA-binding transcription factors (TFs) play a central role in the gene expression of all organisms, from viruses to humans, including bacteria and archaea. The role of these proteins is the fate of gene expression in the context of environmental challenges. Because thousands of genomes have been sequenced to date, predictions of the encoded proteins are validated through the use of bioinformatics tools to obtain the necessary experimental, posterior knowledge. In this chapter, we describe three approaches to identify TFs in protein sequences. The first approach integrates the results of sequence comparisons and PFAM assignments, using as reference a manually curated collection of TFs. The second approach considers the prediction of DNA-binding structures, such as the classical helix-turn-helix (HTH); and the third approach considers a deep learning model. We suggest that all approaches must be considered together to increase the possibility of identifying new TFs in bacterial and archaeal genomes.

Assuntos

RESUMO

In this work, we inferred the gene regulatory network (GRN) of the fungus Fusarium oxysporum by using the regulatory networks of Aspergillus nidulans FGSC A4, Neurospora crassa OR74A, Saccharomyces cerevisiae S288c, and Fusarium graminearum PH-1 as templates for sequence comparisons. Topological properties to infer the role of transcription factors (TFs) and to identify functional modules were calculated in the GRN. From these analyzes, five TFs were identified as hubs, including FOXG_04688 and FOXG_05432, which regulate 2,404 and 1,864 target genes, respectively. In addition, 16 communities were identified in the GRN, where the largest contains 1,923 genes and the smallest contains 227 genes. Finally, the genes associated with virulence were extracted from the GRN and exhaustively analyzed, and we identified a giant module with ten TFs and 273 target genes, where the most highly connected node corresponds to the transcription factor FOXG_05265, homologous to the putative bZip transcription factor CPTF1 of Claviceps purpurea, which is involved in ergotism disease that affects cereal crops and grasses. The results described in this work can be used for the study of gene regulation in this organism and open the possibility to explore putative genes associated with virulence against their host.

RESUMO

BACKGROUND: Archaea are a vast and unexplored domain. Bioinformatic techniques might enlighten the path to a higher quality genome annotation in varied organisms. Promoter sequences of archaea have the action of a plethora of proteins upon it. The conservation found in a structural level of the binding site of proteins such as TBP, TFB, and TFE aids RNAP-DNA stabilization and makes the archaeal promoter prone to be explored by statistical and machine learning techniques. RESULTS AND DISCUSSIONS: In this study, experimentally verified promoter sequences of the organisms Haloferax volcanii, Sulfolobus solfataricus, and Thermococcus kodakarensis were converted into DNA duplex stability attributes (i.e. numerical variables) and were classified through Artificial Neural Networks and an in-house statistical method of classification, being tested with three forms of controls. The recognition of these promoters enabled its use to validate unannotated promoter sequences in other organisms. As a result, the binding site of basal transcription factors was located through a DNA duplex stability codification. Additionally, the classification presented satisfactory results (above 90%) among varied levels of control. CONCLUDING REMARKS: The classification models were employed to perform genomic annotation into the archaea Aciduliprofundum boonei and Thermofilum pendens, from which potential promoters have been identified and uploaded into public repositories.

Assuntos

RESUMO

Penicillium echinulatum 2HH is an ascomycete well known for its production of cellulolytic enzymes. Understanding lignocellulolytic and sugar uptake systems is essential to obtain efficient fungi strains for the production of bioethanol. In this study we performed a genome-wide functional annotation of carbohydrate-active enzymes and sugar transporters involved in the lignocellulolytic system of P. echinulatum 2HH and S1M29 strains (wildtype and mutant, respectively) and eleven related fungi. Additionally, signal peptide and orthology prediction were carried out. We encountered a diverse assortment of cellulolytic enzymes in P. echinulatum, especially in terms of ß-glucosidases and endoglucanases. Other enzymes required for the breakdown of cellulosic biomass were also found, including cellobiohydrolases, lytic cellulose monooxygenases and cellobiose dehydrogenases. The S1M29 mutant, which is known to produce an increased cellulase activity, and the 2HH wild type strain of P. echinulatum did not show significant differences between their enzymatic repertoire. Nevertheless, we unveiled an amino acid substitution for a predicted intracellular ß-glucosidase of the mutant, which might contribute to hyperexpression of cellulases through a cellodextrin induction pathway. Most of the P. echinulatum enzymes presented orthologs in P. oxalicum 114-2, supporting the presence of highly similar cellulolytic mechanisms and a close phylogenetic relationship between these fungi. A phylogenetic analysis of intracellular ß-glucosidases and sugar transporters allowed us to identify several proteins potentially involved in the accumulation of intracellular cellodextrins. These may prove valuable targets in the genetic engineering of P. echinulatum focused on industrial cellulases production. Our study marks an important step in characterizing and understanding the molecular mechanisms employed by P. echinulatum in the enzymatic hydrolysis of lignocellulosic biomass.

Assuntos

RESUMO

Bioinformatics is a very important informatics tool for health and biological sciences, focusing on biological data management. The objective of this work was to perform a bibliometric analysis regarding the development of Mexican bioinformatics. An exhaustive revision of the literature associated with Mexican bioinformatics in a period of 25-years was performed. Bibliometric tools, such as performance analysis and science mapping were included in the analysis. We identified the main actors as well as the structure and dynamics of Mexican bioinformatics. Some of the main findings were as follows: the thematic structure in the field is defined by the research lines of outstanding authors; the outstanding collaborations of Mexican institutions with foreign countries and institutions are influenced by the geographic proximity and binational agreements, as well as philanthropic and academic programs that promote collaborations, and there is an inclination for health issues promoted by public health financing and philanthropic organizations. It is identified that publications had an explosion since 2012, we consider that this growth may be influenced by the democratization of data, derived from the mass sequencing of biological molecules stored in public databases.