Comprehensive comparative analysis of the periodontal pathogen Porphyromonas gingivalis: exploring the pan-genome, the reconstruction of the gene regulatory network and genome-scale metabolic network.

Porphyromonas gingivalis is a nonmotile, obligate anaerobic, Gram-negative bacterium known for its association with periodontal disease and its involvement in systemic diseases such as atherosclerosis, cardiovascular disease, colon cancer, and Alzheimer's disease. This bacterium produces several virulence factors, including capsules, fimbriae, lipopolysaccharides, proteolytic enzymes, and hemagglutinins. A comparative genomic analysis revealed the open pangenome of P. gingivalis and identified complete type IV secretion systems in strain KCOM2805 and almost complete type VI secretion systems in strains KCOM2798 and ATCC49417, which is a new discovery as previous studies did not find the proteins involved in secretion systems IV and VI. Conservation of some virulence factors between different strains was observed, regardless of their genetic diversity and origin. In addition, we performed for the first time a reconstruction analysis of the gene regulatory network, identifying transcription factors and proteins involved in the regulatory mechanisms of bacterial pathogenesis. In particular, QseB regulates the expression of hemagglutinin and arginine deaminase, while Rex may suppress the release of gingipain through interactions with PorV and the formatum/nitrate transporter. Our study highlights the central role of conserved virulence factors and regulatory pathways, particularly QseB and Rex, in P. gingivalis and provides insights into potential therapeutic targets.

Screening and Structural Characterization of Heat Shock Response Elements (HSEs) in Entamoeba histolytica Promoters.

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.

Differential Impact of CD43 and CD28 on T-Cell Differentiation Depending on the Order of Engagement with the TCR.

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.

Machine learning and statistics shape a novel path in archaeal promoter annotation.

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.

Characterization and distribution of GHRH, PACAP, TRH, SST and IGF1 mRNAs in the green iguana.

The somatotropic axis (SA) regulates numerous aspects of vertebrate physiology such as development, growth, and metabolism and has influence on several tissues including neural, immune, reproductive and gastric tract. Growth hormone (GH) is a key component of SA, it is synthesized and released mainly by pituitary somatotrophs, although now it is known that virtually all tissues can express GH, which, in addition to its well-described endocrine roles, also has autocrine/paracrine/intracrine actions. In the pituitary, GH expression is regulated by several hypothalamic neuropeptides including GHRH, PACAP, TRH and SST. GH, in turn, regulates IGF1 synthesis in several target tissues, adding complexity to the system since GH effects can be exerted either directly or mediated by IGF1. In reptiles, little is known about the SA components and their functional interactions. The aim of this work was to characterize the mRNAs of the principal SA components in the green iguana and to develop the tools that allow the study of the structural and functional evolution of this system in reptiles. By employing RT-PCR and RACE, the cDNAs encoding for GHRH, PACAP, TRH, SST and IGF1 were amplified and sequenced. Results showed that these cDNAs coded for the corresponding protein precursors of 154, 170, 243, 113, and 131 amino acids, respectively. Of these, GHRH, PACAP, SST and IGF1 precursors exhibited a high structural conservation with respect to its counterparts in other vertebrates. On the other hand, iguana's TRH precursor showed 7 functional copies of mature TRH (pyr-QHP-NH2), as compared to 4 and 6 copies of TRH in avian and mammalian proTRH sequences, respectively. It was found that in addition to its primary production site (brain for GHRH, PACAP, TRH and SST, and liver for IGF1), they were also expressed in other peripheral tissues, i.e. testes and ovaries expressed all the studied mRNAs, whereas TRH and IGF1 mRNAs were observed ubiquitously in all tissues considered. These results show that the main SA components in reptiles of the Squamata Order maintain a good structural conservation among vertebrate phylogeny, and suggest important physiological interactions (endocrine, autocrine and/or paracrine) between them due to their wide peripheral tissue expression.

RegulonDB version 9.0: high-level integration of gene regulation, coexpression, motif clustering and beyond.

RegulonDB (http://regulondb.ccg.unam.mx) is one of the most useful and important resources on bacterial gene regulation,as it integrates the scattered scientific knowledge of the best-characterized organism, Escherichia coli K-12, in a database that organizes large amounts of data. Its electronic format enables researchers to compare their results with the legacy of previous knowledge and supports bioinformatics tools and model building. Here, we summarize our progress with RegulonDB since our last Nucleic Acids Research publication describing RegulonDB, in 2013. In addition to maintaining curation up-to-date, we report a collection of 232 interactions with small RNAs affecting 192 genes, and the complete repertoire of 189 Elementary Genetic Sensory-Response units (GENSOR units), integrating the signal, regulatory interactions, and metabolic pathways they govern. These additions represent major progress to a higher level of understanding of regulated processes. We have updated the computationally predicted transcription factors, which total 304 (184 with experimental evidence and 120 from computational predictions); we updated our position-weight matrices and have included tools for clustering them in evolutionary families. We describe our semiautomatic strategy to accelerate curation, including datasets from high-throughput experiments, a novel coexpression distance to search for 'neighborhood' genes to known operons and regulons, and computational developments.

Genome misclassification of Klebsiella variicola and Klebsiella quasipneumoniae isolated from plants, animals and humans.

OBJECTIVE: Due to the fact that K. variicola, K. quasipneumoniae and K. pneumoniae are closely related bacterial species, misclassification can occur due to mistakes either in normal biochemical tests or during submission to public databases. The objective of this work was to identify K. variicola and K. quasipneumoniae genomes misclassified in GenBank database. MATERIALS AND METHODS: Both rpoB phylogenies and average nucleotide identity (ANI) were used to identify a significant number of misclassified Klebsiella spp. genomes. RESULTS: Here we report an update of K. variicola and K. Quasipneumoniae genomes correctly classified and a list of isolated genomes obtained from humans, plants, animals and insects, described originally as K. pneumoniae or K. variicola, but known now to be misclassified. CONCLUSIONS: This work contributes to recognize the extensive presence of K. variicola and K. quasipneumoniae isolates in diverse sites and samples.

Dissecting the protein architecture of DNA-binding transcription factors in bacteria and archaea.

Gene regulation at the transcriptional level is a central process in all organisms where DNA-binding transcription factors play a fundamental role. This class of proteins binds specifically at DNA sequences, activating or repressing gene expression as a function of the cell's metabolic status, operator context and ligand-binding status, among other factors, through the DNA-binding domain (DBD). In addition, TFs may contain partner domains (PaDos), which are involved in ligand binding and protein-protein interactions. In this work, we systematically evaluated the distribution, abundance and domain organization of DNA-binding TFs in 799 non-redundant bacterial and archaeal genomes. We found that the distributions of the DBDs and their corresponding PaDos correlated with the size of the genome. We also identified specific combinations between the DBDs and their corresponding PaDos. Within each class of DBDs there are differences in the actual angle formed at the dimerization interface, responding to the presence/absence of ligands and/or crystallization conditions, setting the orientation of the resulting helices and wings facing the DNA. Our results highlight the importance of PaDos as central elements that enhance the diversity of regulatory functions in all bacterial and archaeal organisms, and our results also demonstrate the role of PaDos in sensing diverse signal compounds. The highly specific interactions between DBDs and PaDos observed in this work, together with our structural analysis highlighting the difficulty in predicting both inter-domain geometry and quaternary structure, suggest that these systems appeared once and evolved with diverse duplication events in all the analysed organisms.

Identification of beta-Lactamases and beta-Lactam-Related Proteins in Human Pathogenic Bacteria using a Computational Search Approach.

A systematic analysis of beta-lactamases based on comparative proteomics has not been performed thus far. In this report, we searched for the presence of beta-lactam-related proteins in 591 bacterial proteomes belonging to 52 species that are pathogenic to humans. The amino acid sequences for 19 different types of beta-lactamases (ACT, CARB, CifA, CMY, CTX, FOX, GES, GOB, IMP, IND, KPC, LEN, OKP, OXA, OXY, SHV, TEM, NDM, and VIM) were obtained from the ARG-ANNOT database and were used to construct 19 HMM profiles, which were used to identify potential beta-lactamases in the completely sequenced bacterial proteomes. A total of 2877 matches that included the word "beta-lactamase" and/or "penicillin" in the functional annotation and/or in any of its regions were obtained. These enzymes were mainly described as "penicillin-binding proteins," "beta-lactamases," and "metallo-beta-lactamases" and were observed in 47 of the 52 species studied. In addition, proteins classified as "beta-lactamases" were observed in 39 of the species included. A positive correlation between the number of beta-lactam-related proteins per species and the proteome size was observed (R 0.78, P < 0.00001). This correlation partially explains the high presence of beta-lactam-related proteins in large proteomes, such as Nocardia brasiliensis, Bacillus anthracis, and Mycobacterium tuberculosis, along with their absence in small proteomes, such as Chlamydia spp. and Mycoplasma spp. We detected only five types of beta-lactamases (TEM, SHV, CTX, IMP, and OXA) and other related proteins in particular species that corresponded with those reported in the literature. We additionally detected other potential species-specific beta-lactamases that have not yet been reported. In the future, better results will be achieved due to more accurate sequence annotations and a greater number of sequenced genomes.

Transcriptomic analysis of the process of biofilm formation in Rhizobium etli CFN42.

Organisms belonging to the genus Rhizobium colonize leguminous plant roots and establish a mutually beneficial symbiosis. Biofilms are structured ecosystems in which microbes are embedded in a matrix of extracellular polymeric substances, and their development is a multistep process. The biofilm formation processes of R. etli CFN42 were analyzed at an early (24-h incubation) and mature stage (72 h), comparing cells in the biofilm with cells remaining in the planktonic stage. A genome-wide microarray analysis identified 498 differentially regulated genes, implying that expression of ~8.3 % of the total R. etli gene content was altered during biofilm formation. In biofilms-attached cells, genes encoding proteins with diverse functions were overexpressed including genes involved in membrane synthesis, transport and chemotaxis, repression of flagellin synthesis, as well as surface components (particularly exopolysaccharides and lipopolysaccharides), in combination with the presence of activators or stimulators of N-acyl-homoserine lactone synthesis This suggests that R. etli is able to sense surrounding environmental conditions and accordingly regulate the transition from planktonic and biofilm growth. In contrast, planktonic cells differentially expressed genes associated with transport, motility (flagellar and twitching) and inhibition of exopolysaccharide synthesis. To our knowledge, this is the first report of nodulation and nitrogen assimilation-related genes being involved in biofilm formation in R. etli. These results contribute to the understanding of the physiological changes involved in biofilm formation by bacteria.

Identification of DNA Methyltransferase Genes in Human Pathogenic Bacteria by Comparative Genomics.

DNA methylation plays an important role in gene expression and virulence in some pathogenic bacteria. In this report, we describe DNA methyltransferases (MTases) present in human pathogenic bacteria and compared them with related species, which are not pathogenic or less pathogenic, based in comparative genomics. We performed a search in the KEGG database of the KEGG database orthology groups associated with adenine and cytosine DNA MTase activities (EC: 2.1.1.37, EC: 2.1.1.113 and EC: 2.1.1.72) in 37 human pathogenic species and 18 non/less pathogenic relatives and performed comparisons of the number of these MTases sequences according to their genome size, the DNA MTase type and with their non-less pathogenic relatives. We observed that Helicobacter pylori and Neisseria spp. presented the highest number of MTases while ten different species did not present a predicted DNA MTase. We also detected a significant increase of adenine MTases over cytosine MTases (2.19 vs. 1.06, respectively, p < 0.001). Adenine MTases were the only MTases associated with restriction modification systems and DNA MTases associated with type I restriction modification systems were more numerous than those associated with type III restriction modification systems (0.84 vs. 0.17, p < 0.001); additionally, there was no correlation with the genome size and the total number of DNA MTases, indicating that the number of DNA MTases is related to the particular evolution and lifestyle of specific species, regulating the expression of virulence genes in some pathogenic bacteria.

The lifestyle of prokaryotic organisms influences the repertoire of promiscuous enzymes.

The metabolism of microbial organisms and its diversity are partly the result of an adaptation process to the characteristics of the environments that they inhabit. In this work, we analyze the influence of lifestyle on the content of promiscuous enzymes in 761 nonredundant bacterial and archaeal genomes. Promiscuous enzymes were defined as those proteins whose catalytic activities are defined by two or more different Enzyme Commission (E.C.) numbers. The genomes analyzed were categorized into four lifestyles for their exhaustive comparisons: free-living, extremophiles, pathogens, and intracellular. From these analyses we found that free-living organisms have larger genomes and an enrichment of promiscuous enzymes. In contrast, intracellular organisms showed smaller genomes and the lesser proportion of promiscuous enzymes. On the basis of our data, we show that the proportion of promiscuous enzymes in an organism is mainly influenced by the lifestyle, where fluctuating environments promote its emergence. Finally, we evidenced that duplication processes occur preferentially in metabolism of free-living and extremophiles species.

Electrostatic analysis of bacterial expansins.

Expansins are a family of proteins with plant cell wall remodeling-activity, which bind cell wall components through hydrophobic and electrostatic interactions. A shallow area on the surface of the protein serves as the polysaccharide binding site (PBS) and it is composed of conserved residues. However, electric charge differences on the opposite face of the PBS produce basic, neutral, or acidic proteins. An analysis of forty-four bacterial expansins, homologues of BsEXLX1, revealed two main groups defined by: (a) the presence or absence of disulfide bonds; and (b) by the proteins isoelectric point (pI). We determined the location of the residues responsible for the pI on the structure of representative expansins. Our results suggest that the electric charge at the opposite site of the PBS may help in substrate differentiation among expansins from different species; in addition, electrostatic polarization between the front and the back of the molecule could affect expansin activity on cellulose.

The alignment of enzymatic steps reveals similar metabolic pathways and probable recruitment events in Gammaproteobacteria.

BACKGROUND: It is generally accepted that gene duplication followed by functional divergence is one of the main sources of metabolic diversity. In this regard, there is an increasing interest in the development of methods that allow the systematic identification of these evolutionary events in metabolism. Here, we used a method not based on biomolecular sequence analysis to compare and identify common and variable routes in the metabolism of 40 Gammaproteobacteria species. METHOD: The metabolic maps deposited in the KEGG database were transformed into linear Enzymatic Step Sequences (ESS) by using the breadth-first search algorithm. These ESS represent subsequent enzymes linked to each other, where their catalytic activities are encoded in the Enzyme Commission numbers. The ESS were compared in an all-against-all (pairwise comparisons) approach by using a dynamic programming algorithm, leaving only a set of significant pairs. RESULTS AND CONCLUSION: From these comparisons, we identified a set of functionally conserved enzymatic steps in different metabolic maps, in which cell wall components and fatty acid and lysine biosynthesis were included. In addition, we found that pathways associated with biosynthesis share a higher proportion of similar ESS than degradation pathways and secondary metabolism pathways. Also, maps associated with the metabolism of similar compounds contain a high proportion of similar ESS, such as those maps from nucleotide metabolism pathways, in particular the inosine monophosphate pathway. Furthermore, diverse ESS associated with the low part of the glycolysis pathway were identified as functionally similar to multiple metabolic pathways. In summary, our comparisons may help to identify similar reactions in different metabolic pathways and could reinforce the patchwork model in the evolution of metabolism in Gammaproteobacteria.

Comparative genomics of nucleotide metabolism: a tour to the past of the three cellular domains of life.

BACKGROUND: Nucleotide metabolism is central to all biological systems, due to their essential role in genetic information and energy transfer, which in turn suggests its possible presence in the last common ancestor (LCA) of Bacteria, Archaea and Eukarya. In this context, elucidation of the contribution of the origin and diversification of de novo and salvage pathways of nucleotide metabolism will allow us to understand the links between the enzymatic steps associated with the LCA and the emergence of the first metabolic pathways. RESULTS: In this work, the taxonomical distribution of the enzymes associated with nucleotide metabolism was evaluated in 1,606 complete genomes. 151 sequence profiles associated with 120 enzymatic reactions were used. The evaluation was based on profile comparisons, using RPS-Blast. Organisms were clustered based on their taxonomical classifications, in order to obtain a normalized measure of the taxonomical distribution of enzymes according to the average of presence/absence of enzymes per genus, which in turn was used for the second step, to calculate the average presence/absence of enzymes per Clade. CONCLUSION: From these analyses, it was suggested that divergence at the enzymatic level correlates with environmental changes and related modifications of the cell wall and membranes that took place during cell evolution. Specifically, the divergence of the 5-(carboxyamino) imidazole ribonucleotide mutase to phosphoribosylaminoimidazole carboxylase could be related to the emergence of multicellularity in eukaryotic cells. In addition, segments of salvage and de novo pathways were probably complementary in the LCA to the synthesis of purines and pyrimidines. We also suggest that a large portion of the pathway to inosine 5'-monophosphate (IMP) in purines could have been involved in thiamine synthesis or its derivatives in early stages of cellular evolution, correlating with the fact that these molecules may have played an active role in the protein-RNA world. The analysis presented here provides general observations concerning the adaptation of the enzymatic steps in the early stages of the emergence of life and the LCA.

Characterization of pituitary growth hormone and its receptor in the green iguana (Iguana iguana).

Pituitary growth hormone (GH) has been studied in most vertebrate groups; however, only a few studies have been carried out in reptiles. Little is known about pituitary hormones in the order Squamata, to which the green iguana (gi) belongs. In this work, we characterized the hypophysis of Iguana iguana morphologically. The somatotrophs (round cells of 7.6-10 µm containing 250- to 300-nm secretory granules where the giGH is stored) were found, by immunohistochemistry and in situ hybridization, exclusively in the caudal lobe of the pars distalis, whereas the lactotrophs were distributed only in the rostral lobe. A pituitary giGH-like protein was obtained by immuno-affinity chromatography employing a heterologous antibody against chicken GH. giGH showed molecular heterogeneity (22, 44, and 88 kDa by SDS-PAGE/Western blot under non-reducing conditions and at least four charge variants (pIs 6.2, 6.5, 6.9, 7.4) by isoelectric focusing. The pituitary giGH cDNA (1016 bp), amplified by PCR and RACE, encodes a pre-hormone of 218 aa, of which 190 aa correspond to the mature protein and 28 aa to the signal peptide. The giGH receptor cDNA was also partially sequenced. Phylogenetic analyses of the amino acid sequences of giGH and giGHR homologs in vertebrates suggest a parallel evolution and functional relationship between the GH and its receptor.

DeepReg: a deep learning hybrid model for predicting transcription factors in eukaryotic and prokaryotic genomes.

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.

The interaction of InvF-RNAP is mediated by the chaperone SicA in Salmonella sp: an in silico prediction.

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.

Comparative genomics sheds light on transcription factor-mediated regulation in the extreme acidophilic Acidithiobacillia representatives.

Extreme acidophiles thrive in acidic environments, confront a multitude of challenges, and demonstrate remarkable adaptability in their metabolism to cope with the ever-changing environmental fluctuations, which encompass variations in temperature, pH levels, and the availability of electron acceptors and donors. The survival and proliferation of members within the Acidithiobacillia class rely on the deployment of transcriptional regulatory systems linked to essential physiological traits. The study of these transcriptional regulatory systems provides valuable insights into critical processes, such as energy metabolism and nutrient assimilation, and how they integrate into major genetic-metabolic circuits. In this study, we examined the transcriptional regulatory repertoires and potential interactions of forty-three Acidithiobacillia complete and draft genomes, encompassing nine species. To investigate the function and diversity of Transcription Factors (TFs) and their DNA Binding Sites (DBSs), we conducted a genome-wide comparative analysis, which allowed us to identify these regulatory elements in representatives of Acidithiobacillia. We classified TFs into gene families and compared their occurrence among all representatives, revealing conservation patterns across the class. The results identified conserved regulators for several pathways, including iron and sulfur oxidation, the main pathways for energy acquisition, providing new evidence for viable regulatory interactions and branch-specific conservation in Acidithiobacillia. The identification of TFs and DBSs not only corroborates existing experimental information for selected species, but also introduces novel candidates for experimental validation. Moreover, these promising candidates have the potential for further extension to new representatives within the class.

Identification of modules and key genes associated with breast cancer subtypes through network analysis.

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.