Comparative GO: a web application for comparative gene ontology and gene ontology-based gene selection in bacteria

Pungency related gene network in Allium sativum L., response to sulfur treatments

Pungency of garlic (Allium sativum L.) is generated from breakdown of the alk(en)yl cysteine sulphoxide (CSO), alliin and its subsequent breakdown to allicin under the activity of alliinase (All). Based on recent evidence, two other important genes including Sulfite reductase (SiR) and Superoxide dismutase (SOD) are thought to be related to sulfur metabolism. These three gene functions are in sulfate assimilation pathway. However, whether it is involved in stress response in crops is largely unknown. In this research, the order and priority of simultaneous expression of three genes including All, SiR and SOD were measured on some garlic ecotypes of Iran, collected from Zanjan, Hamedan and Gilan, provinces under sulfur concentrations (0, 6, 12, 24 and 60 g/ per experimental unit: pot) using real-time quantitative PCR (RT-qPCR) analysis. For understanding the network interactions between studied genes and other related genes, in silico gene network analysis was constructed to investigate various mechanisms underlying stimulation of A. sativum L. to cope with imposed sulfur. Complicated network including TF-TF, miRNA-TF, and miRNA-TF-gene, was split into sub-networks to have a deeper insight. Analysis of q-RT-PCR data revealed the highest expression in All and SiR genes respectively. To distinguish and select significant pathways in sulfur metabolism, RESNET Plant database of Pathway Studio software v.10 (Elsevier), and other relative data such as chemical reactions, TFs, miRNAs, enzymes, and small molecules were extracted. Complex sub-network exhibited plenty of routes between stress response and sulfate assimilation pathway. Even though Alliinase did not display any connectivity with other stress response genes, it showed binding relation with lectin functional class, as a result of which connected to leucine zipper, exocellulase, peroxidase and ARF functional class indirectly. Integration network of these genes revealed their involvement in various biological processes such as, RNA splicing, stress response, gene silencing by miRNAs, and epigenetic. The findings of this research can be used to extend further research on the garlic metabolic engineering, garlic stress related genes, and also reducing or enhancing the activity of the responsible genes for garlic pungency for health benefits and industry demands.

Quantitative Proteomic Analysis of Cyanide and Mercury Detoxification by Pseudomonas pseudoalcaligenes CECT 5344

The cyanide-degrading bacterium Pseudomonas pseudoalcaligenes CECT 5344 uses cyanide and different metal-cyanide complexes as the sole nitrogen source. Under cyanotrophic conditions, this strain was able to grow with up to 100 μM mercury, which was accumulated intracellularly. A quantitative proteomic analysis by liquid chromatography-tandem mass spectrometry (LC-MS/MS) has been applied to unravel the molecular basis of the detoxification of both cyanide and mercury by the strain CECT 5344, highlighting the relevance of the cyanide-insensitive alternative oxidase CioAB and the nitrilase NitC in the tolerance and assimilation of cyanide, independently of the presence or absence of mercury. Proteins overrepresented in the presence of cyanide and mercury included mercury transporters, mercuric reductase MerA, transcriptional regulator MerD, arsenate reductase and arsenical resistance proteins, thioredoxin reductase, glutathione S-transferase, proteins related to aliphatic sulfonates metabolism and sulfate transport, hemin import transporter, and phosphate starvation induced protein PhoH, among others. A transcriptional study revealed that from the six putative merR genes present in the genome of the strain CECT 5344 that could be involved in the regulation of mercury resistance/detoxification, only the merR2 gene was significantly induced by mercury under cyanotrophic conditions. A bioinformatic analysis allowed the identification of putative MerR2 binding sites in the promoter regions of the regulatory genes merR5, merR6, arsR, and phoR, and also upstream from the structural genes encoding glutathione S-transferase (fosA and yghU), dithiol oxidoreductase (dsbA), metal resistance chaperone (cpxP), and amino acid/peptide extruder involved in quorum sensing (virD), among others. IMPORTANCE Cyanide, mercury, and arsenic are considered very toxic chemicals that are present in nature as cocontaminants in the liquid residues generated by different industrial activities like mining. Considering the huge amounts of toxic cyanide- and mercury-containing wastes generated at a large scale and the high biotechnological potential of P. pseudoalcaligenes CECT 5344 in the detoxification of cyanide present in these industrial wastes, in this work, proteomic, transcriptional, and bioinformatic approaches were used to characterize the molecular response of this bacterium to cyanide and mercury, highlighting the mechanisms involved in the simultaneous detoxification of both compounds. The results generated could be applied for developing bioremediation strategies to detoxify wastes cocontaminated with cyanide, mercury, and arsenic, such as those generated at a large scale in the mining industry.

New RoxS sRNA Targets Identified in Bacillus subtilis by Pulsed SILAC

Non-coding RNAs (sRNA) play a key role in controlling gene expression in bacteria, typically by base-pairing with ribosome binding sites to block translation. The modification of ribosome traffic along the mRNA generally affects its stability. However, a few cases have been described in bacteria where sRNAs can affect translation without a major impact on mRNA stability. To identify new sRNA targets in Bacillus subtilis potentially belonging to this class of mRNAs, we used pulsed-SILAC (stable isotope labeling by amino acids in cell culture) to label newly synthesized proteins after short expression of the RoxS sRNA, the best characterized sRNA in this bacterium. RoxS sRNA was previously shown to interfere with the expression of genes involved in central metabolism, permitting control of the NAD+/NADH ratio in B. subtilis. In this study, we confirmed most of the known targets of RoxS, showing the efficiency of the method. We further expanded the number of mRNA targets encoding enzymes of the TCA cycle and identified new targets. One of these is YcsA, a tartrate dehydrogenase that uses NAD+ as co-factor, in excellent agreement with the proposed role of RoxS in management of NAD+/NADH ratio in Firmicutes. IMPORTANCE Non-coding RNAs (sRNA) play an important role in bacterial adaptation and virulence. The identification of the most complete set of targets for these regulatory RNAs is key to fully identifying the perimeter of its function(s). Most sRNAs modify both the translation (directly) and mRNA stability (indirectly) of their targets. However, sRNAs can also influence the translation efficiency of the target primarily, with little or no impact on mRNA stability. The characterization of these targets is challenging. We describe here the application of the pulsed SILAC method to identify such targets and obtain the most complete list of targets for a defined sRNA.

Proteomic Analysis of Arsenic Resistance during Cyanide Assimilation by Pseudomonas pseudoalcaligenes CECT 5344

Wastewater from mining and other industries usually contains arsenic and cyanide, two highly toxic pollutants, thereby creating the need to develop bioremediation strategies. Here, molecular mechanisms triggered by the simultaneous presence of cyanide and arsenite were analyzed by quantitative proteomics, complemented with qRT-PCR analysis and determination of analytes in the cyanide-assimilating bacterium Pseudomonas pseudoalcaligenes CECT 5344. Several proteins encoded by two ars gene clusters and other Ars-related proteins were up-regulated by arsenite, even during cyanide assimilation. Although some proteins encoded by the cio gene cluster responsible for cyanide-insensitive respiration decreased in the presence of arsenite, the nitrilase NitC required for cyanide assimilation was unaffected, thus allowing bacterial growth with cyanide and arsenic. Two complementary As-resistance mechanisms were developed in this bacterium, the extrusion of As(III) and its extracellular sequestration in biofilm, whose synthesis increased in the presence of arsenite, and the formation of organoarsenicals such as arseno-phosphoglycerate and methyl-As. Tetrahydrofolate metabolism was also stimulated by arsenite. In addition, the ArsH2 protein increased in the presence of arsenite or cyanide, suggesting its role in the protection from oxidative stress caused by both toxics. These results could be useful for the development of bioremediation strategies for industrial wastes co-contaminated with cyanide and arsenic.

Knockdown Experiment Reveals an Essential GTPase CgtA's Involvement in Growth, Viability, Motility, Morphology, and Persister Phenotypes in Vibrio cholerae

CgtA is an essential bacterial GTPase consisting of a highly conserved N-terminal Spo0B-associated GTP-binding protein (Obg) domain, a central GTPase domain, and a variable C-terminal domain (CTD). This study reports global changes in the proteome and transcriptome of wild-type (Wt) versus full-length CgtA-depleted Vibrio cholerae in minimal media. Comparative transcriptome sequencing (RNA-Seq), followed by comparative proteomic analyses, revealed that the knockdown of cgtA significantly altered expressions of 311 proteins involved in diverse cellular activities, many of which are associated with the survival of V. cholerae. Various intracellular functional roles of CgtA in growth, viability, motility, morphology, and persister phenotypes in V. cholerae are revealed based on subsequent confirmatory experiments. Furthermore, a more sustained mRNA expression pattern of cgtA in a minimal medium than in a rich medium was also observed for Wt V. cholerae, where the highest level of mRNA expression of cgtA was observed during the logarithmic growth phase. Thereby, we propose that minimal medium-associated reduced growth rate coupled with cgtA depletion aggravates the intracellular stress in V. cholerae, interrupting vital cellular processes. The functional role of the CTD in V. cholerae is not fully understood. Hence, to specifically investigate the intracellular role of the 57-amino-acid-long CTD of CgtAVC, the CTD-only portion of CgtA was deleted. Subsequent proteomics studies revealed an altered expression of 240 proteins in the CgtA(ΔCTD) mutant, having major overlap with the full-length cgtA-deleted condition. Overall, our study reveals an alternative facet of the survival mechanism of V. cholerae during nutritional downshift as per the concomitant consequences of cgtA depletion. IMPORTANCE It is very important that we must find new drug target proteins from multidrug-resistant human-pathogenic organisms like V. cholerae. CgtA is among such potential candidates, and here, we are reporting about some newly identified cellular roles of it that are important for the survival of V. cholerae. Briefly, we knocked down the full-length cgtA gene, as well as did a partial deletion of this gene from the V. cholerae genome followed by RNA-Seq and proteomics studies. Results from our study revealed up- and downregulation of several known and unknown genes and proteins as the effect of the cgtA knockdown experiment. Also, we have presented some interesting observations that are linked with cgtA for growth, viability, motility, morphology, and persister phenotypes in V. cholerae. Our study enhances the importance of CgtA and paves the way for further exploration based on our provided data.

Holistic view of biological nitrogen fixation and phosphorus mobilization in Azotobacter chroococcum NCIMB 8003

Nitrogen (N) and phosphorus (P) deficiencies are two of the most agronomic problems that cause significant decrease in crop yield and quality. N and P chemical fertilizers are widely used in current agriculture, causing environmental problems and increasing production costs. Therefore, the development of alternative strategies to reduce the use of chemical fertilizers while maintaining N and P inputs are being investigated. Although dinitrogen is an abundant gas in the atmosphere, it requires biological nitrogen fixation (BNF) to be transformed into ammonium, a nitrogen source assimilable by living organisms. This process is bioenergetically expensive and, therefore, highly regulated. Factors like availability of other essential elements, as phosphorus, strongly influence BNF. However, the molecular mechanisms of these interactions are unclear. In this work, a physiological characterization of BNF and phosphorus mobilization (PM) from an insoluble form (Ca₃(PO₄)₂) in Azotobacter chroococcum NCIMB 8003 was carried out. These processes were analyzed by quantitative proteomics in order to detect their molecular requirements and interactions. BNF led to a metabolic change beyond the proteins strictly necessary to carry out the process, including the metabolism related to other elements, like phosphorus. Also, changes in cell mobility, heme group synthesis and oxidative stress responses were observed. This study also revealed two phosphatases that seem to have the main role in PM, an exopolyphosphatase and a non-specific alkaline phosphatase PhoX. When both BNF and PM processes take place simultaneously, the synthesis of nitrogenous bases and L-methionine were also affected. Thus, although the interdependence is still unknown, possible biotechnological applications of these processes should take into account the indicated factors.

The NtrYX Two-Component System of Paracoccus denitrificans Is Required for the Maintenance of Cellular Iron Homeostasis and for a Complete Denitrification under Iron-Limited Conditions

Denitrification consists of the sequential reduction of nitrate to nitrite, nitric oxide, nitrous oxide, and dinitrogen. Nitrous oxide escapes to the atmosphere, depending on copper availability and other environmental factors. Iron is also a key element because many proteins involved in denitrification contain iron-sulfur or heme centers. The NtrYX two-component regulatory system mediates the responses in a variety of metabolic processes, including denitrification. A quantitative proteomic analysis of a Paracoccus denitrificans NtrY mutant grown under denitrifying conditions revealed the induction of different TonB-dependent siderophore transporters and proteins related to iron homeostasis. This mutant showed lower intracellular iron content than the wild-type strain, and a reduced growth under denitrifying conditions in iron-limited media. Under iron-rich conditions, it releases higher concentrations of siderophores and displayes lower nitrous oxide reductase (NosZ) activity than the wild-type, thus leading to nitrous oxide emission. Bioinformatic and qRT-PCR analyses revealed that NtrYX is a global transcriptional regulatory system that responds to iron starvation and, in turn, controls expression of the iron-responsive regulators fur, rirA, and iscR, the denitrification regulators fnrP and narR, the nitric oxide-responsive regulator nnrS, and a wide set of genes, including the cd₁-nitrite reductase NirS, nitrate/nitrite transporters and energy electron transport proteins.

The Exoproteome of Staphylococcus pasteuri Isolated from Cervical Mucus during the Estrus Phase in Water Buffalo (Bubalus bubalis)

Bacterial extracellular proteins participate in the host cell communication by virtue of the modulation of pathogenicity, commensalism and mutualism. Studies on the microbiome of cervical mucus of the water buffalo (Bubalus bubalis) have shown the occurrence of Staphylococcus pasteuri and that the presence of this bacterium is indicative of various physiological and reproductive states in the host. Recently, S. pasteuri has been isolated from the cervical mucus of the buffalo during the different phases of estrous cycle, and has proved to be much more pronounced during the estrus phase. The basis underlying the availability of a significantly increased S. pasteuri population, specifically during the estrus phase, is not known. Consequently, it is important to determine the significance of the specific abundance of S. pasteuri during the estrus phase of the buffalo host, particularly from the perspective of whether this bacterial species is capable of contributing to sexual communication via its extracellular proteins and volatiles. Therefore, the relevance of S. pasteuri exoproteome in the buffalo cervical mucus during the estrus phase was analyzed using LC-MS/MS. As many as 219 proteins were identified, among which elongation factor Tu (EF-Tu), 60-kDa chaperonin (Cpn60), enolase, fructose-bisphosphate aldolase class 1 (FBP aldolase), enoyl-[acyl-carrier-protein] reductase [NADPH] (ENR) and lipoprotein (Lpp) were the functionally important candidates. Most of the proteins present in the exoproteome of S. pasteuri were those involved in cellular–metabolic functions, as well as catalytic- and binding activities. Moreover, computational studies of Lpp have shown enhanced interaction with volatiles such as acetic-, butanoic-, isovaleric- and valeric acids, which were identified in the cervical mucus S. pasteuri culture supernatant. The present findings suggest that S. pasteuri extracellular proteins may play an important role in buffalo sexual communication during the estrus phase.

Splice-disrupt genomic variants in prostate cancer

Background

Splice-disrupt genomic variants are one of the causes of cancer-causing errors in gene expression. Little is known about splice-disrupt genomic variants.

Methods and results

Here, pattern of splice-disrupt variants was investigated using 21,842,764 genomic variants in different types of prostate cancer. A particular attention was paid to genomic locations of splice-disrupt variants on target genes. HLA-A in prostate cancer, MSR1 in familial prostate cancer, and EGFR in both castration-resistant prostate cancer and metastatic castration-resistant had the highest allele frequencies of splice-disrupt variations. Some splice-disrupt variants, located on coding sequences of NCOR2, PTPRC, and CRP, were solely present in the advanced metastatic castration-resistant prostate cancer. High-risk splice-disrupt variants were identified based on computationally calculated Polymorphism Phenotyping (PolyPhen), Sorting Intolerant From Tolerant (SIFT), and Genomic Evolutionary Rate Profiling (GERP) + + scores as well as the recorded clinical significance in dbSNP database of NCBI. Functional annotation of damaging splice-disrupt variants highlighted important cancer-associated functions, including endocrine resistance, lipid metabolic process, steroid metabolic process, regulation of mitotic cell cycle, and regulation of metabolic process. This is the first study that profiles the splice-disrupt genomic variants and their target genes in prostate cancer. Literature mining based variant analysis highlighted the importance of rs1800716 variant, located on the CYP2D6 gene, involved in a range of important functions, such as RNA spicing, drug interaction, death, and urotoxicity.

Conclusions

This is the first study that profiles the splice-disrupt genomic variants and their target genes in different types of prostate cancer. Unravelling alternative splicing opens a new avenue towards the establishment of new diagnostic and prognostic markers for prostate cancer progression and metastasis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-022-07257-9.

MOET: a web-based gene set enrichment tool at the Rat Genome Database for multiontology and multispecies analyses

Biological interpretation of a large amount of gene or protein data is complex. Ontology analysis tools are imperative in finding functional similarities through overrepresentation or enrichment of terms associated with the input gene or protein lists. However, most tools are limited by their ability to do ontology-specific and species-limited analyses. Furthermore, some enrichment tools are not updated frequently with recent information from databases, thus giving users inaccurate, outdated or uninformative data. Here, we present MOET or the Multi-Ontology Enrichment Tool (v.1 released in April 2019 and v.2 released in May 2021), an ontology analysis tool leveraging data that the Rat Genome Database (RGD) integrated from in-house expert curation and external databases including the National Center for Biotechnology Information (NCBI), Mouse Genome Informatics (MGI), The Kyoto Encyclopedia of Genes and Genomes (KEGG), The Gene Ontology Resource, UniProt-GOA, and others. Given a gene or protein list, MOET analysis identifies significantly overrepresented ontology terms using a hypergeometric test and provides nominal and Bonferroni corrected P-values and odds ratios for the overrepresented terms. The results are shown as a downloadable list of terms with and without Bonferroni correction, and a graph of the P-values and number of annotated genes for each term in the list. MOET can be accessed freely from https://rgd.mcw.edu/rgdweb/enrichment/start.html.

Alternative Pathway for 3-Cyanoalanine Assimilation in Pseudomonas pseudoalcaligenes CECT5344 under Noncyanotrophic Conditions

3-Cyanoalanine and cyanohydrins are intermediate nitriles produced in cyanide degradation pathways in plants and bacteria. 3-Cyanoalanine is generated from cyanide by the 3-cyanoalanine synthase, an enzyme mainly characterized in cyanogenic plants. NIT4-type nitrilases use 3-cyanoalanine as a substrate, forming ammonium and aspartate. In some organisms, this enzyme also generates asparagine through an additional nitrile hydratase activity. The alkaliphilic bacterium Pseudomonas pseudoalcaligenes CECT5344 assimilates cyanide through an intermediate cyanohydrin, which is further converted into ammonium by the nitrilase NitC. This bacterium also contains three additional nitrilases, including Nit4. In this work, a proteomic analysis of P. pseudoalcaligenes CECT5344 cells grown with 3-cyanoalanine as the sole nitrogen source has revealed the overproduction of different proteins involved in nitrogen metabolism, including the nitrilase NitC. In contrast, the nitrilase Nit4 was not induced by 3-cyanoalanine, and it was only overproduced in cells grown with a cyanide-containing jewelry-manufacturing residue. Phenotypes of single and double mutant strains defective in nit4 or/and nitC revealed the implication of the nitrilase NitC in the assimilation of 3-cyanoalanine and suggest that the 3-cyanoalanine assimilation pathway in P. pseudoalcaligenes CECT5344 depends on the presence or absence of cyanide. When cyanide is present, 3-cyanoalanine is assimilated via Nit4, but in the absence of cyanide, a novel pathway for 3-cyanoalanine assimilation, in which the nitrilase NitC uses the nitrile generated after deamination of the α-amino group from 3-cyanoalanine, is proposed. IMPORTANCE Nitriles are organic cyanides with important industrial applications, but they are also found in nature. 3-Cyanoalanine is synthesized by plants and some bacteria to detoxify cyanide from endogenous or exogenous sources, but this nitrile may be also involved in other processes such as stress tolerance, nitrogen and sulfur metabolism, and signaling. The cyanide-degrading bacterium Pseudomonas pseudoalcaligenes CECT5344 grows with 3-cyanoalanine as the sole nitrogen source, but it does not use this nitrile as an intermediate in the cyanide assimilation pathway. In this work, a quantitative proteomic analysis by liquid chromatography-tandem mass spectrometry (LC-MS/MS) was performed to study, for the first time, the response to 3-cyanoalanine at the proteomic level. Proteomic data, together with phenotypes of different nitrilase-defective mutants of P. pseudoalcaligenes CECT5344, provide evidence that in the absence of cyanide, the nitrilase Nit4 is not involved in 3-cyanoalanine assimilation, and instead, the nitrilase NitC participates in a novel alternative 3-cyanoalanine assimilation pathway.

Metabolic Remodeling during Nitrogen Fixation in Zymomonas mobilis

Zymomonas mobilis is an ethanologenic bacterium currently being developed for production of advanced biofuels. Recent studies have shown that Z. mobilis can fix dinitrogen gas (N2) as a sole nitrogen source. During N2 fixation, Z. mobilis exhibits increased biomass-specific rates of ethanol production. In order to better understand the physiology of Z. mobilis during N2 fixation and during changes in ammonium (NH4+) availability, we performed liquid chromatography-mass spectrometry (LC-MS)-based targeted metabolomics and shotgun proteomics under three regimes of nitrogen availability: continuous N2 fixation, gradual NH4+ depletion, and acute NH4+ addition to N2-fixing cells. We report dynamic changes in abundance of proteins and metabolites related to nitrogen fixation, motility, ammonium assimilation, amino acid biosynthesis, nucleotide biosynthesis, isoprenoid biosynthesis, and Entner-Doudoroff (ED) glycolysis, providing insight into the regulatory mechanisms that control these processes in Z. mobilis. Our analysis identified potential physiological mechanisms that may contribute to increased specific ethanol production during N2 fixation, including decreased activity of biosynthetic pathways, increased protein abundance of alcohol dehydrogenase (ADHI), and increased thermodynamic favorability of the ED pathway. Of particular relevance to advanced biofuel production, we found that intermediates in the methylerythritol phosphate (MEP) pathway for isoprenoid biosynthesis were depleted during N2 fixation, coinciding with decreased protein abundance of deoxyxylulose 5-phosphate synthase (DXS), the first enzyme in the pathway. This implies that DXS protein abundance serves as a native control point in regulating MEP pathway activity in Z. mobilis. The results of this study will inform metabolic engineering to further develop Z. mobilis as a platform organism for biofuel production. IMPORTANCE Biofuels and bioproducts have the potential to serve as environmentally sustainable replacements for petroleum-derived fuels and commodity molecules. Advanced fuels such as higher alcohols and isoprenoids are more suitable gasoline replacements than bioethanol. Developing microbial systems to generate advanced biofuels requires metabolic engineering to reroute carbon away from ethanol and other native products and toward desired pathways, such as the MEP pathway for isoprenoid biosynthesis. However, rational engineering of microbial metabolism relies on understanding metabolic control points, in terms of both enzyme activity and thermodynamic favorability. In Z. mobilis, the factors that control glycolytic rates, ethanol production, and isoprenoid production are still not fully understood. In this study, we performed metabolomic, proteomic, and thermodynamic analysis of Z. mobilis during N2 fixation. This analysis identified key changes in metabolite levels, enzyme abundance, and glycolytic thermodynamic favorability that occurred during changes in NH4+ availability, helping to inform future efforts in metabolic engineering.

Approaching In Vivo Models of Pneumococcus-Host Interaction: Insights into Surface Proteins, Capsule Production, and Extracellular Vesicles

Infections caused by the Gram-positive bacterium Streptococcus pneumoniae have become a major health problem worldwide because of their high morbidity and mortality rates, especially in developing countries. This microorganism colonizes the human upper respiratory tract and becomes pathogenic under certain circumstances, which are not well known. In the interaction with the host, bacterial surface structures and proteins play major roles. To gain knowledge into gradual changes and adaptive mechanisms that this pathogen undergoes from when it enters the host, we mimicked several in vivo situations representing interaction with epithelial and macrophage cells, as well as a condition of presence in blood. Then, we analyzed, in four pneumococcal strains, two major surface structures, the capsule and extracellular vesicles produced by the pneumococci, as well as surface proteins by proteomics, using the “shaving” approach, followed by LC-MS/MS. We found important differences in both surface ultrastructures and proteins among the culture conditions and strains used. Thus, this work provides insights into physiological adaptations of the pneumococcus when it interacts with the host, which may be useful for the design of strategies to combat infections caused by this pathogen.

Effect of pH on the denitrification proteome of the soil bacterium Paracoccus denitrificans PD1222

Denitrification is a respiratory process by which nitrate is reduced to dinitrogen. Incomplete denitrification results in the emission of the greenhouse gas nitrous oxide and this is potentiated in acidic soils, which display reduced denitrification rates and high N2O/N2 ratios compared to alkaline soils. In this work, impact of pH on the proteome of the soil denitrifying bacterium Paracoccus denitrificans PD1222 was analysed with nitrate as sole energy and nitrogen source under anaerobic conditions at pH ranging from 6.5 to 7.5. Quantitative proteomic analysis revealed that the highest difference in protein representation was observed when the proteome at pH 6.5 was compared to the reference proteome at pH 7.2. However, this difference in the extracellular pH was not enough to produce modification of intracellular pH, which was maintained at 6.5 ± 0.1. The biosynthetic pathways of several cofactors relevant for denitrification and nitrogen assimilation like cobalamin, riboflavin, molybdopterin and nicotinamide were negatively affected at pH 6.5. In addition, peptide representation of reductases involved in nitrate assimilation and denitrification were reduced at pH 6.5. Data highlight the strong negative impact of pH on NosZ synthesis and intracellular copper content, thus impairing active NosZ assembly and, in turn, leading to elevated nitrous oxide emissions.

Synthetase of the methyl donor S-adenosylmethionine from nitrogen-fixing α-rhizobia can bind functionally diverse RNA species

Function of bacterial small non-coding RNAs (sRNAs) and overall RNA metabolism is largely shaped by a vast diversity of RNA-protein interactions. However, in non-model bacteria with defined non-coding transcriptomes the sRNA interactome remains almost unexplored. We used affinity chromatography to capture proteins associated in vivo with MS2-tagged trans-sRNAs that regulate nutrient uptake (AbcR2 and NfeR1) and cell cycle (EcpR1) mRNAs by antisense-based translational inhibition in the nitrogen-fixing α-rhizobia Sinorhizobium meliloti. The three proteomes were rather distinct, with that of EcpR1 particularly enriched in cell cycle-related enzymes, whilst sharing several transcription/translation-related proteins recurrently identified associated with sRNAs. Strikingly, MetK, the synthetase of the major methyl donor S-adenosylmethionine, was reliably recovered as a binding partner of the three sRNAs, which reciprocally co-immunoprecipitated with a FLAG-tagged MetK variant. Induced (over)expression of the trans-sRNAs and MetK depletion did not influence canonical riboregulatory traits, `for example, protein titration or sRNA stability, respectively. An in vitro filter assay confirmed binding of AbcR2, NfeR1 and EcpR1 to MetK and further revealed interaction of the protein with other non-coding and coding transcripts but not with the 5S rRNA. These findings uncover a broad specificity for RNA binding as an unprecedented feature of this housekeeping prokaryotic enzyme.

MegaGO: A Fast Yet Powerful Approach to Assess Functional Gene Ontology Similarity across Meta-Omics Data Sets

The study of microbiomes has gained in importance over the past few years and has led to the emergence of the fields of metagenomics, metatranscriptomics, and metaproteomics. While initially focused on the study of biodiversity within these communities, the emphasis has increasingly shifted to the study of (changes in) the complete set of functions available in these communities. A key tool to study this functional complement of a microbiome is Gene Ontology (GO) term analysis. However, comparing large sets of GO terms is not an easy task due to the deeply branched nature of GO, which limits the utility of exact term matching. To solve this problem, we here present MegaGO, a user-friendly tool that relies on semantic similarity between GO terms to compute the functional similarity between multiple data sets. MegaGO is high performing: Each set can contain thousands of GO terms, and results are calculated in a matter of seconds. MegaGO is available as a web application at https://megago.ugent.be and is installable via pip as a standalone command line tool and reusable software library. All code is open source under the MIT license and is available at https://github.com/MEGA-GO/.

Global Transcriptional Response of Methylorubrum extorquens to Formaldehyde Stress Expands the Role of EfgA and Is Distinct from Antibiotic Translational Inhibition

The potency and indiscriminate nature of formaldehyde reactivity upon biological molecules make it a universal stressor. However, some organisms such as Methylorubrum extorquens possess means to rapidly and effectively mitigate formaldehyde-induced damage. EfgA is a recently identified formaldehyde sensor predicted to halt translation in response to elevated formaldehyde as a means to protect cells. Herein, we investigate growth and changes in gene expression to understand how M. extorquens responds to formaldehyde with and without the EfgA-formaldehyde-mediated translational response, and how this mechanism compares to antibiotic-mediated translation inhibition. These distinct mechanisms of translation inhibition have notable differences: they each involve different specific players and in addition, formaldehyde also acts as a general, multi-target stressor and a potential carbon source. We present findings demonstrating that in addition to its characterized impact on translation, functional EfgA allows for a rapid and robust transcriptional response to formaldehyde and that removal of EfgA leads to heightened proteotoxic and genotoxic stress in the presence of increased formaldehyde levels. We also found that many downstream consequences of translation inhibition were shared by EfgA-formaldehyde- and kanamycin-mediated translation inhibition. Our work uncovered additional layers of regulatory control enacted by functional EfgA upon experiencing formaldehyde stress, and further demonstrated the importance this protein plays at both transcriptional and translational levels in this model methylotroph.

The identification of co-expressed gene modules in Streptococcus pneumonia from colonization to infection to predict novel potential virulence genes

Background

Streptococcus pneumonia (pneumococcus) is a human bacterial pathogen causing a range of mild to severe infections. The complicated transcriptome patterns of pneumococci during the colonization to infection process in the human body are usually determined by measuring the expression of essential virulence genes and the comparison of pathogenic with non-pathogenic bacteria through microarray analyses. As systems biology studies have demonstrated, critical co-expressing modules and genes may serve as key players in biological processes. Generally, Sample Progression Discovery (SPD) is a computational approach traditionally used to decipher biological progression trends and their corresponding gene modules (clusters) in different clinical samples underlying a microarray dataset. The present study aimed to investigate the bacterial gene expression pattern from colonization to severe infection periods (specimens isolated from the nasopharynx, lung, blood, and brain) to find new genes/gene modules associated with the infection progression. This strategy may lead to finding novel gene candidates for vaccines or drug design.

Results

The results included essential genes whose expression patterns varied in different bacterial conditions and have not been investigated in similar studies.

Conclusions

In conclusion, the SPD algorithm, along with differentially expressed genes detection, can offer new ways of discovering new therapeutic or vaccine targeted gene products.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-020-02059-0.

The role of the msaABCR operon in implant-associated chronic osteomyelitis in Staphylococcus aureus USA300 LAC

BACKGROUND

The msaABCR operon regulates several staphylococcal phenotypes such as biofilm formation, capsule production, protease production, pigmentation, antibiotic resistance, and persister cells formation. The msaABCR operon is required for maintaining the cell wall integrity via affecting peptidoglycan cross-linking. The msaABCR operon also plays a role in oxidative stress defense mechanism, which is required to facilitate persistent and recurrent staphylococcal infections. Staphylococcus aureus is the most frequent cause of chronic implant-associated osteomyelitis (OM). The CA-MRSA USA300 strains are predominant in the United States and cause severe infections, including bone and joint infections.

RESULTS

The USA300 LAC strain caused significant bone damage, as evidenced by the presence of severe bone necrosis with multiple foci of sequestra and large numbers of multinucleated osteoclasts. Intraosseous survival and biofilm formation on the K-wires by USA300 LAC strains was pronounced. However, the msaABCR deletion mutant was attenuated. We observed minimal bone necrosis, with no evidence of intramedullary abscess and/or fibrosis, along reduced intraosseous bacterial population and significantly less biofilm formation on the K-wires by the msaABCR mutant. microCT analysis of infected bone showed significant bone loss and damage in the USA300 LAC and complemented strain, whereas the msaABCR mutant's effect was reduced. In addition, we observed increased osteoblasts response and new bone formation around the K-wires in the bone infected by the msaABCR mutant. Whole-cell proteomics analysis of msaABCR mutant cells showed significant downregulation of proteins, cell adhesion factors, and virulence factors that interact with osteoblasts and are associated with chronic OM caused by S. aureus.

CONCLUSION

This study showed that deletion of msaABCR operon in USA300 LAC strain lead to defective biofilm in K-wire implants, decreased intraosseous survival, and reduced cortical bone destruction. Thus, msaABCR plays a role in implant-associated chronic osteomyelitis by regulating extracellular proteases, cell adhesions factors and virulence factors. However additional studies are required to further define the contribution of msaABCR-regulated molecules in osteomyelitis pathogenesis.

Role of the Dihydrodipicolinate Synthase DapA1 on Iron Homeostasis During Cyanide Assimilation by the Alkaliphilic Bacterium Pseudomonas pseudoalcaligenes CECT5344

Cyanide is a toxic compound widely used in mining and jewelry industries, as well as in the synthesis of many different chemicals. Cyanide toxicity derives from its high affinity for metals, which causes inhibition of relevant metalloenzymes. However, some cyanide-degrading microorganisms like the alkaliphilic bacterium Pseudomonas pseudoalcaligenes CECT5344 may detoxify hazardous industrial wastewaters that contain elevated cyanide and metal concentrations. Considering that iron availability is strongly reduced in the presence of cyanide, mechanisms for iron homeostasis should be required for cyanide biodegradation. Previous omic studies revealed that in the presence of a cyanide-containing jewelry residue the strain CECT5344 overproduced the dihydrodipicolinate synthase DapA1, a protein involved in lysine metabolism that also participates in the synthesis of dipicolinates, which are excellent metal chelators. In this work, a dapA1 ^- mutant of P. pseudoalcaligenes CECT5344 has been generated and characterized. This mutant showed reduced growth and cyanide consumption in media with the cyanide-containing wastewater. Intracellular levels of metals like iron, copper and zinc were increased in the dapA1 ^- mutant, especially in cells grown with the jewelry residue. In addition, a differential quantitative proteomic analysis by LC-MS/MS was carried out between the wild-type and the dapA1 ^- mutant strains in media with jewelry residue. The mutation in the dapA1 gene altered the expression of several proteins related to urea cycle and metabolism of arginine and other amino acids. Additionally, the dapA1 ^- mutant showed increased levels of the global nitrogen regulator PII and the glutamine synthetase. This proteomic study has also highlighted that the DapA1 protein is relevant for cyanide resistance, oxidative stress and iron homeostasis response, which is mediated by the ferric uptake regulator Fur. DapA1 is required to produce dipicolinates that could act as iron chelators, conferring protection against oxidative stress and allowing the regeneration of Fe-S centers to reactivate cyanide-damaged metalloproteins.

Transcriptomic Analysis of the Brucella melitensis Rev.1 Vaccine Strain in an Acidic Environment: Insights Into Virulence Attenuation

The live attenuated Brucella melitensis Rev.1 (Elberg-originated) vaccine strain is widely used to control the zoonotic infection brucellosis in small ruminants, but the molecular mechanisms underlying the attenuation of this strain have not been fully characterized. Following their uptake by the host cell, Brucella replicate inside a membrane-bound compartment—the Brucella-containing vacuole—whose acidification is essential for the survival of the pathogen. Therefore, identifying the genes that contribute to the survival of Brucella in acidic environments will greatly assist our understanding of its molecular pathogenic mechanisms and of the attenuated virulence of the Rev.1 strain. Here, we conducted a comprehensive comparative transcriptome analysis of the Rev.1 vaccine strain against the virulent reference strain 16M in cultures grown under either normal or acidic conditions. We found 403 genes that respond differently to acidic conditions in the two strains (FDR < 0.05, fold change ≥ 2). These genes are involved in crucial cellular processes, including metabolic, biosynthetic, and transport processes. Among the highly enriched genes that were downregulated in Rev.1 under acidic conditions were acetyl-CoA synthetase, aldehyde dehydrogenase, cell division proteins, a cold-shock protein, GroEL, and VirB3. The downregulation of these genes may explain the attenuated virulence of Rev.1 and provide new insights into the virulence mechanisms of Brucella.

Non-coding and coding genomic variants distinguish prostate cancer, castration-resistant prostate cancer, familial prostate cancer, and metastatic castration-resistant prostate cancer from each other

A considerable number of deposited variants has provided new possibilities for knowledge discovery in different types of prostate cancer. Here, we analyzed variants located on 3'UTR, 5'UTR, CDs, Intergenic, and Intronic regions in castration-resistant prostate cancer (8496 variants), familial prostate cancer (3241 variants), metastatic castration-resistant prostate cancer (3693 variants), and prostate cancer (16599 variants). Chromosome regions 10p15-p14 and 2p13 were highly enriched (P < 0.00001) for variants located in 3'UTR, 5'UTR, CDs, intergenic, and intronic regions in castration-resistant prostate cancer. In contrast, 10p15-p14, 10q23.3, 12q13.11, 13q12.3, 1q25, and 8p22 regions were enriched (P < 0.001) in familial prostate cancer. In metastatic castration-resistant prostate cancer, 10p15-p14, 10q23.3, 11q22-q23, 14q21.1, and 14q32.13 were highly variant regions (P < 0.001). Chromosome 2 and chromosome 1 hosted many enriched variant regions. AKR1C3, BRCA1, BRCA2, CHGA, CYP19A1, HOXB13, KLK3, and PTEN contained the highest number of 3'UTR, 5'UTR, CDs, Intergenic, and Intronic variants. Network analysis showed that these genes are upstream of important functions including prostate gland development, tumor recurrence, prostate cancer-specific survival, tumor progression, cancer mortality, long-term survival, cancer recurrence, angiogenesis, and AR. Interestingly, all of EGFR, JAK2, NR3C1, PDZD2, and SEMA3C genes had single nucleotide polymorphisms (SNP) in castration-resistant prostate cancer, consistent with high selection pressure on these genes during drug treatment and consequent resistance. High occurrence of variants in 3'UTRs suggests the importance of regulatory variants in different types of prostate cancer; an area that has been neglected compared with coding variants. This study provides a comprehensive overview of genomic regions contributing to different types of prostate cancer.

Identification of Small RNA-Protein Partners in Plant Symbiotic Bacteria

The identification of the protein partners of bacterial small noncoding RNAs (sRNAs) is essential to understand the mechanistic principles and functions of riboregulation in prokaryotic cells. Here, we describe an optimized affinity chromatography protocol that enables purification of in vivo formed sRNA-protein complexes in Sinorhizobium meliloti, a genetically tractable nitrogen-fixing plant symbiotic bacterium. The procedure requires the tagging of the desired sRNA with the MS2 aptamer, which is affinity-captured by the MS2-MBP protein conjugated to an amylose resin. As proof of principle, we show recovery of the RNA chaperone Hfq associated to the strictly Hfq-dependent AbcR2 trans-sRNA. This method can be applied for the investigation of sRNA-protein interactions on a broad range of genetically tractable α-proteobacteria.

Unified Transcriptomic Signature of Arbuscular Mycorrhiza Colonization in Roots of Medicago truncatula by Integration of Machine Learning, Promoter Analysis, and Direct Merging Meta-Analysis

Plant root symbiosis with Arbuscular mycorrhizal (AM) fungi improves uptake of water and mineral nutrients, improving plant development under stressful conditions. Unraveling the unified transcriptomic signature of a successful colonization provides a better understanding of symbiosis. We developed a framework for finding the transcriptomic signature of Arbuscular mycorrhiza colonization and its regulating transcription factors in roots of Medicago truncatula. Expression profiles of roots in response to AM species were collected from four separate studies and were combined by direct merging meta-analysis. Batch effect, the major concern in expression meta-analysis, was reduced by three normalization steps: Robust Multi-array Average algorithm, Z-standardization, and quartiling normalization. Then, expression profile of 33685 genes in 18 root samples of Medicago as numerical features, as well as study ID and Arbuscular mycorrhiza type as categorical features, were mined by seven models: RELIEF, UNCERTAINTY, GINI INDEX, Chi Squared, RULE, INFO GAIN, and INFO GAIN RATIO. In total, 73 genes selected by machine learning models were up-regulated in response to AM (Z-value difference > 0.5). Feature weighting models also documented that this signature is independent from study (batch) effect. The AM inoculation signature obtained was able to differentiate efficiently between AM inoculated and non-inoculated samples. The AP2 domain class transcription factor, GRAS family transcription factors, and cyclin-dependent kinase were among the highly expressed meta-genes identified in the signature. We found high correspondence between the AM colonization signature obtained in this study and independent RNA-seq experiments on AM colonization, validating the repeatability of the colonization signature. Promoter analysis of upregulated genes in the transcriptomic signature led to the key regulators of AM colonization, including the essential transcription factors for endosymbiosis establishment and development such as NF-YA factors. The approach developed in this study offers three distinct novel features: (I) it improves direct merging meta-analysis by integrating supervised machine learning models and normalization steps to reduce study-specific batch effects; (II) seven attribute weighting models assessed the suitability of each gene for the transcriptomic signature which contributes to robustness of the signature (III) the approach is justifiable, easy to apply, and useful in practice. Our integrative framework of meta-analysis, promoter analysis, and machine learning provides a foundation to reveal the transcriptomic signature and regulatory circuits governing Arbuscular mycorrhizal symbiosis and is transferable to the other biological settings.

Cross-Species Meta-Analysis of Transcriptomic Data in Combination With Supervised Machine Learning Models Identifies the Common Gene Signature of Lactation Process

Lactation, a physiologically complex process, takes place in mammary gland after parturition. The expression profile of the effective genes in lactation has not comprehensively been elucidated. Herein, meta-analysis, using publicly available microarray data, was conducted identify the differentially expressed genes (DEGs) between pre- and post-peak milk production. Three microarray datasets of Rat, Bos Taurus, and Tammar wallaby were used. Samples related to pre-peak (n = 85) and post-peak (n = 24) milk production were selected. Meta-analysis revealed 31 DEGs across the studied species. Interestingly, 10 genes, including MRPS18B, SF1, UQCRC1, NUCB1, RNF126, ADSL, TNNC1, FIS1, HES5 and THTPA, were not detected in original studies that highlights meta-analysis power in biosignature discovery. Common target and regulator analysis highlighted the high connectivity of CTNNB1, CDD4 and LPL as gene network hubs. As data originally came from three different species, to check the effects of heterogeneous data sources on DEGs, 10 attribute weighting (machine learning) algorithms were applied. Attribute weighting results showed that the type of organism had no or little effect on the selected gene list. Systems biology analysis suggested that these DEGs affect the milk production by improving the immune system performance and mammary cell growth. This is the first study employing both meta-analysis and machine learning approaches for comparative analysis of gene expression pattern of mammary glands in two important time points of lactation process. The finding may pave the way to use of publically available to elucidate the underlying molecular mechanisms of physiologically complex traits such as lactation in mammals.

Exploring the Denitrification Proteome of Paracoccus denitrificans PD1222

Denitrification is a respiratory process that produces nitrous oxide as an intermediate, which may escape to the atmosphere before its reduction to dinitrogen through the nitrous oxide reductase NosZ. In this work, the denitrification process carried out by Paracoccus denitrificans PD1222 has been explored through a quantitative proteomic analysis. Under anaerobic conditions, with nitrate as sole nitrogen source, the synthesis of all the enzymes involved in denitrification, the respiratory nitrate, nitrite, nitric oxide, and nitrous oxide reductases, was increased. However, the periplasmic and assimilatory nitrate reductases decreased. Synthesis of transporters for alcohols, D-methionine, sulfate and copper, most of the enzymes involved in the tricarboxylic acid cycle, and proteins involved in other metabolic processes like lysine catabolism, fatty acids degradation and acetyl-CoA synthesis, was increased during denitrification in P. denitrificans PD1222. As consequence, an enhanced production of the central metabolite acetyl-CoA was observed. After establishing the key features of the denitrification proteome, its changes by the influence of a competitive electron acceptor, oxygen, or competitive nitrogen source, ammonium, were evaluated.

First gene-ontology enrichment analysis based on bacterial coregenome variants: insights into adaptations of Salmonella serovars to mammalian- and avian-hosts

Background

Many of the bacterial genomic studies exploring evolution processes of the host adaptation focus on the accessory genome describing how the gains and losses of genes can explain the colonization of new habitats. Consequently, we developed a new approach focusing on the coregenome in order to describe the host adaptation of Salmonella serovars.

Methods

In the present work, we propose bioinformatic tools allowing (i) robust phylogenetic inference based on SNPs and recombination events, (ii) identification of fixed SNPs and InDels distinguishing homoplastic and non-homoplastic coregenome variants, and (iii) gene-ontology enrichment analyses to describe metabolic processes involved in adaptation of Salmonella enterica subsp. enterica to mammalian- (S. Dublin), multi- (S. Enteritidis), and avian- (S. Pullorum and S. Gallinarum) hosts.

Results

The ‘VARCall’ workflow produced a robust phylogenetic inference confirming that the monophyletic clade S. Dublin diverged from the polyphyletic clade S. Enteritidis which includes the divergent clades S. Pullorum and S. Gallinarum (i). The scripts ‘phyloFixedVar’ and ‘FixedVar’ detected non-synonymous and non-homoplastic fixed variants supporting the phylogenetic reconstruction (ii). The scripts ‘GetGOxML’ and ‘EveryGO’ identified representative metabolic pathways related to host adaptation using the first gene-ontology enrichment analysis based on bacterial coregenome variants (iii).

Conclusions

We propose in the present manuscript a new coregenome approach coupling identification of fixed SNPs and InDels with regards to inferred phylogenetic clades, and gene-ontology enrichment analysis in order to describe the adaptation of Salmonella serovars Dublin (i.e. mammalian-hosts), Enteritidis (i.e. multi-hosts), Pullorum (i.e. avian-hosts) and Gallinarum (i.e. avian-hosts) at the coregenome scale. All these polyvalent Bioinformatic tools can be applied on other bacterial genus without additional developments.

Electronic supplementary material

The online version of this article (10.1186/s12866-017-1132-1) contains supplementary material, which is available to authorized users.

msaABCR operon is involved in persister cell formation in Staphylococcus aureus

Background

Persister cells comprise a phenotypic variant that shows extreme antibiotic tolerance resulting in treatment failures of bacterial infections. While this phenomenon has posed a great threat in public health, mechanisms underlying their formation in Staphylococcus aureus remain largely unknown. Increasing evidences of the presence of persister cells in recalcitrant infections underscores the great urgency to unravel the mechanism by which these cells develop. Previously, we characterized msaABCR operon that plays roles in regulation of virulence, biofilm development and antibiotic resistance. We also characterized the function of MsaB protein and showed that MsaB is a putative transcription factor that binds target DNA in response to nutrients availability.

Results

In this study, we compared the number of persister cell in wild type, msaABCR deletion mutant and the complemented strain in two backgrounds USA300 LAC and Mu50. Herein, we report that msaABCR deletion mutant forms significantly less number of persister cells relative to wild type after challenge with various antibiotics in planktonic and biofilm growth conditions. Complementation of the msaABCR operon restored wild type phenotype. Combined antibiotic therapy along with msaABCR deletion significantly improves the killing kinetics of stationary phase and biofilm S. aureus cells. Transcriptomics analysis showed that msaABCR regulates several metabolic genes, transcription factors, transporters and enzymes that may play role in persister cells formation, which we seek to define in the future.

Conclusions

This study presented a new regulator, msaABCR operon, that is involved in the persister cells formation, which is a poorly understood in S. aureus. Indeed, we showed that msaABCR deletion significantly reduces the persister cells formation in all growth phases tested. Although, we have not yet defined the mechanism, we have shown that msaABCR regulates several metabolic, transporters, and extracellular proteases genes that have been previously linked with persister cells formation in other bacterial systems. Taken together, this study showed that inactivation of the msaABCR operon enhances the effectiveness of antibiotics for the treatment of S. aureus infections, especially in context of persister cells.

Electronic supplementary material

The online version of this article (10.1186/s12866-017-1129-9) contains supplementary material, which is available to authorized users.

A novel hypothesis-unbiased method for Gene Ontology enrichment based on transcriptome data

Gene Ontology (GO) classification of statistically significantly differentially expressed genes is commonly used to interpret transcriptomics data as a part of functional genomic analysis. In this approach, all significantly expressed genes contribute equally to the final GO classification regardless of their actual expression levels. Gene expression levels can significantly affect protein production and hence should be reflected in GO term enrichment. Genes with low expression levels can also participate in GO term enrichment through cumulative effects. In this report, we have introduced a new GO enrichment method that is suitable for multiple samples and time series experiments that uses a statistical outlier test to detect GO categories with special patterns of variation that can potentially identify candidate biological mechanisms. To demonstrate the value of our approach, we have performed two case studies. Whole transcriptome expression profiles of Salmonella enteritidis and Alzheimer's disease (AD) were analysed in order to determine GO term enrichment across the entire transcriptome instead of a subset of differentially expressed genes used in traditional GO analysis. Our result highlights the key role of inflammation related functional groups in AD pathology as granulocyte colony-stimulating factor receptor binding, neuromedin U binding, and interleukin were remarkably upregulated in AD brain when all using all of the gene expression data in the transcriptome. Mitochondrial components and the molybdopterin synthase complex were identified as potential key cellular components involved in AD pathology.

The Zebrafish Equivalent of Alzheimer's Disease-Associated PRESENILIN Isoform PS2V Regulates Inflammatory and Other Responses to Hypoxic Stress

Dominant mutations in the PRESENILIN genes PSEN1 and PSEN2 cause familial Alzheimer's disease (fAD) that usually shows onset before 65 years of age. In contrast, genetic variation at the PSEN1 and PSEN2 loci does not appear to contribute to risk for the sporadic, late onset form of the disease (sAD), leading to doubts that these genes play a role in the majority of AD cases. However, a truncated isoform of PSEN2, PS2V, is upregulated in sAD brains and is induced by hypoxia and high cholesterol intake. PS2V can increase γ-secretase activity and suppress the unfolded protein response (UPR), but detailed analysis of its function has been hindered by lack of a suitable, genetically manipulable animal model since mice and rats lack this PRESENILIN isoform. We recently showed that zebrafish possess an isoform, PS1IV, that is cognate to human PS2V. Using an antisense morpholino oligonucleotide, we can block specifically the induction of PS1IV that normally occurs under hypoxia. Here, we exploit this ability to identify gene regulatory networks that are modulated by PS1IV. When PS1IV is absent under hypoxia-like conditions, we observe changes in expression of genes controlling inflammation (particularly sAD-associated IL1B and CCR5), vascular development, the UPR, protein synthesis, calcium homeostasis, catecholamine biosynthesis, TOR signaling, and cell proliferation. Our results imply an important role for PS2V in sAD as a component of a pathological mechanism that includes hypoxia/oxidative stress and support investigation of the role of PS2V in other diseases, including schizophrenia, when these are implicated in the pathology.

Identification of OtpR regulated sRNAs in Brucella melitensis expressed under acidic stress and their roles in pathogenesis and metabolism

Small RNAs (sRNAs) are the small regulatory molecules that regulate various biological processes in bacteria. Though the functions of sRNAs are well documented, very little information is available on the sRNAs of Brucella spp. The otpR is the response regulator of a two-component regulatory system, which plays a significant role in Brucella virulence. In this study, we identified the sRNAs expressed in B. melitensis 16M and its otpR mutant under acidic stress from the RNAseq dataset. We identified 94 trans-encoded and 948 cis-encoded sRNAs in B. melitensis 16M. In B. melitensis 16M ΔotpR under acidic stress 99 trans-encoded and 877 cis-encoded sRNAs were identified. Among these, 12 trans-encoded and 43 cis-encoded sRNAs were upregulated in B. melitensis 16M ΔotpR, with an adjusted P-value≤0.05. The mRNA targets of these sRNAs were predicted. Further, the levels of mRNA targets were examined, and the sRNA-mediated regulatory network was predicted. Functional classification and pathway analysis of mRNA targets provided evidence that sRNAs are involved in different metabolic pathways including carbohydrates, amino acids, lipids, nucleotides transport and metabolism, cell membrane biogenesis and intracellular trafficking of Brucella. We also found that eight transcriptional regulators including a quorum sensing regulator, vjbR are down-regulated by sRNAs. These transcriptional regulators might be responsible for the regulation of several other genes in the otpR mutant. The trans-encoded BsnR88 and cis-encoded BsnR980, BsnR998, BsnR881, BsnR1001, BsnR891, BsnR883, BsnR905 are regulating virB operon genes coding for type IV secretion system (T4SS), which is the major virulence factor of Brucella.

Isolation, Characterization, Cryopreservation of Human Amniotic Stem Cells and Differentiation to Osteogenic and Adipogenic Cells

Human stem cells and progenitor cells can be used to treat cancer and replace dysfunctional cells within a tissue or organ. The objective of this study was to identify the appropriate cells type in regenerative medicine and targeted therapy. As an alternative to embryonic and bone marrow stem cells, we examined human amniotic fluid stem cells (hAFSCs), one of the potential source of multipotent stem cells isolated from both cell pellet (using single-stage method), and supernatant of human amniotic fluid. Source of isolation and unique property of the cells emphasize that these cells are one of the promising new tools in therapeutic field. Double sources for isolation and availability of the left over samples in diagnostic laboratory at the same time have less legal and ethical concerns compared with embryonic stem cell studies. Cells were isolated, cultured for 18^th passage for 6 months and characterized using qPCR and flow cytometry. Cells showed good proliferative ability in culture condition. The cells successfully differentiated into the adipogenic and osteogenic lineages. Based on these findings, amniotic fluid can be considered as an appropriate and convenient source of human amniotic fluid stem cells. These cells provide potential tools for therapeutic applications in the field of regenerative medicine. To get a better understanding of crosstalk between Oct4/NANOG with osteogenesis and adipogenesis, we used network analysis based on Common Targets algorithm and Common Regulators algorithm as well as subnetwork discovery based on gene set enrichment. Network analysis highlighted the possible role of MIR 302A and MIR let-7g. We demonstrated the high expression of MIR 302A and low expression of MIR let7g in hAFSCs by qPCR.

Identification of the key regulating genes of diminished ovarian reserve (DOR) by network and gene ontology analysis

Diminished ovarian reserve (DOR) is one of the reasons for infertility that not only affects both older and young women. Ovarian reserve assessment can be used as a new prognostic tool for infertility treatment decision making. Here, up- and down-regulated gene expression profiles of granulosa cells were analysed to generate a putative interaction map of the involved genes. In addition, gene ontology (GO) analysis was used to get insight intol the biological processes and molecular functions of involved proteins in DOR. Eleven up-regulated genes and nine down-regulated genes were identified and assessed by constructing interaction networks based on their biological processes. PTGS2, CTGF, LHCGR, CITED, SOCS2, STAR and FSTL3 were the key nodes in the up-regulated networks, while the IGF2, AMH, GREM, and FOXC1 proteins were key in the down-regulated networks. MIRN101-1, MIRN153-1 and MIRN194-1 inhibited the expression of SOCS2, while CSH1 and BMP2 positively regulated IGF1 and IGF2. Ossification, ovarian follicle development, vasculogenesis, sequence-specific DNA binding transcription factor activity, and golgi apparatus are the major differential groups between up-regulated and down-regulated genes in DOR. Meta-analysis of publicly available transcriptomic data highlighted the high coexpression of CTGF, connective tissue growth factor, with the other key regulators of DOR. CTGF is involved in organ senescence and focal adhesion pathway according to GO analysis. These findings provide a comprehensive system biology based insight into the aetiology of DOR through network and gene ontology analyses.

DrugMiner: comparative analysis of machine learning algorithms for prediction of potential druggable proteins

Application of computational methods in drug discovery has received increased attention in recent years as a way to accelerate drug target prediction. Based on 443 sequence-derived protein features, we applied the most commonly used machine learning methods to predict whether a protein is druggable as well as to opt for superior algorithm in this task. In addition, feature selection procedures were used to provide the best performance of each classifier according to the optimum number of features. When run on all features, Neural Network was the best classifier, with 89.98% accuracy, based on a k-fold cross-validation test. Among all the algorithms applied, the optimum number of most-relevant features was 130, according to the Support Vector Machine-Feature Selection (SVM-FS) algorithm. This study resulted in the discovery of new drug target which potentially can be employed in cell signaling pathways, gene expression, and signal transduction. The DrugMiner web tool was developed based on the findings of this study to provide researchers with the ability to predict druggable proteins. DrugMiner is freely available at www.DrugMiner.org.

Characterization of Pneumococcal Genes Involved in Bloodstream Invasion in a Mouse Model

Streptococcus pneumoniae (the pneumococcus) continues to account for significant morbidity and mortality worldwide, causing life-threatening diseases such as pneumonia, bacteremia and meningitis, as well as less serious infections such as sinusitis, conjunctivitis and otitis media. Current polysaccharide vaccines are strictly serotype-specific and also drive the emergence of non-vaccine serotype strains. In this study, we used microarray analysis to compare gene expression patterns of either serotype 4 or serotype 6A pneumococci in the nasopharynx and blood of mice, as a model to identify genes involved in invasion of blood in the context of occult bacteremia in humans. In this manner, we identified 26 genes that were significantly up-regulated in the nasopharynx and 36 genes that were significantly up-regulated in the blood that were common to both strains. Gene Ontology classification revealed that transporter and DNA binding (transcription factor) activities constitute the significantly different molecular functional categories for genes up-regulated in the nasopharynx and blood. Targeted mutagenesis of selected genes from both niches and subsequent virulence and pathogenesis studies identified the manganese-dependent superoxide dismutase (SodA) as most likely to be essential for colonization, and the cell wall-associated serine protease (PrtA) as important for invasion of blood. This work extends our previous analyses and suggests that both PrtA and SodA warrant examination in future studies aimed at prevention and/or control of pneumococcal disease.

PanCoreGen - Profiling, detecting, annotating protein-coding genes in microbial genomes

A large amount of genomic data, especially from multiple isolates of a single species, has opened new vistas for microbial genomics analysis. Analyzing the pan-genome (i.e. the sum of genetic repertoire) of microbial species is crucial in understanding the dynamics of molecular evolution, where virulence evolution is of major interest. Here we present PanCoreGen - a standalone application for pan- and core-genomic profiling of microbial protein-coding genes. PanCoreGen overcomes key limitations of the existing pan-genomic analysis tools, and develops an integrated annotation-structure for a species-specific pan-genomic profile. It provides important new features for annotating draft genomes/contigs and detecting unidentified genes in annotated genomes. It also generates user-defined group-specific datasets within the pan-genome. Interestingly, analyzing an example-set of Salmonella genomes, we detect potential footprints of adaptive convergence of horizontally transferred genes in two human-restricted pathogenic serovars - Typhi and Paratyphi A. Overall, PanCoreGen represents a state-of-the-art tool for microbial phylogenomics and pathogenomics study.

CompGO: an R package for comparing and visualizing Gene Ontology enrichment differences between DNA binding experiments

Background

Gene ontology (GO) enrichment is commonly used for inferring biological meaning from systems biology experiments. However, determining differential GO and pathway enrichment between DNA-binding experiments or using the GO structure to classify experiments has received little attention.

Results

Herein, we present a bioinformatics tool, CompGO, for identifying Differentially Enriched Gene Ontologies, called DiEGOs, and pathways, through the use of a z-score derivation of log odds ratios, and visualizing these differences at GO and pathway level. Through public experimental data focused on the cardiac transcription factor NKX2-5, we illustrate the problems associated with comparing GO enrichments between experiments using a simple overlap approach.

Conclusions

We have developed an R/Bioconductor package, CompGO, which implements a new statistic normally used in epidemiological studies for performing comparative GO analyses and visualizing comparisons from .BED data containing genomic coordinates as well as gene lists as inputs. We justify the statistic through inclusion of experimental data and compare to the commonly used overlap method. CompGO is freely available as a R/Bioconductor package enabling easy integration into existing pipelines and is available at: http://www.bioconductor.org/packages/release/bioc/html/CompGO.html packages/release/bioc/html/CompGO.html

Transcriptome Profiling of Wild-Type and pga-Knockout Mutant Strains Reveal the Role of Exopolysaccharide in Aggregatibacter actinomycetemcomitans

Exopolysaccharides have a diverse set of functions in most bacteria including a mechanistic role in protecting bacteria against environmental stresses. Among the many functions attributed to the exopolysaccharides, biofilm formation, antibiotic resistance, immune evasion and colonization have been studied most extensively. The exopolysaccharide produced by many Gram positive as well as Gram negative bacteria including the oral pathogen Aggregatibacter actinomycetemcomitans is the homopolymer of β(1,6)-linked N-acetylglucosamine. Recently, we reported that the PGA-deficient mutant of A. actinomycetemcomitans failed to colonize or induce bone resorption in a rat model of periodontal disease, and the colonization genes, apiA and aae, were significantly down regulated in the mutant strain. To understand the role of exopolysaccharide and the pga locus in the global expression of A. actinomycetemcomitans, we have used comparative transcriptome profiling to identify differentially expressed genes in the wild-type strain in relation to the PGA-deficient strain. Transcriptome analysis revealed that about 50% of the genes are differently expressed (P < 0.05 and fold change >1.5). Our study demonstrated that the absence of the pga locus affects the genes involved in peptidoglycan recycling, glycogen storage, and virulence. Further, using confocal microscopy and plating assays, we show that the viability of pga mutant strain is significantly reduced during biofilm growth. Thus, this study highlights the importance of pga genes and the exopolysaccharide in the virulence of A. actinomycetemcomitans.

Differential expression of seven conserved microRNAs in response to abiotic stress and their regulatory network in Helianthus annuus

Biotic and abiotic stresses affect plant development and production through alternation of the gene expression pattern. Gene expression itself is under the control of different regulators such as miRNAs and transcription factors (TFs). MiRNAs are known to play important roles in regulation of stress responses via interacting with their target mRNAs. Here, for the first time, seven conserved miRNAs, associated with drought, heat, salt and cadmium stresses were characterized in sunflower. The expression profiles of miRNAs and their targets were comparatively analyzed between leaves and roots of plants grown under the mentioned stress conditions. Gene ontology analysis of target genes revealed that they are involved in several important pathways such as auxin and ethylene signaling, RNA mediated silencing and DNA methylation processes. Gene regulatory network highlighted the existence of cross-talks between these stress-responsive miRNAs and the other stress responsive genes in sunflower. Based on network analysis, we suggest that some of these miRNAs in sunflower such as miR172 and miR403 may play critical roles in epigenetic responses to stress. It seems that depending on the stress type, theses miRNAs target several pathways and cellular processes to help sunflower to cope with drought, heat, salt and cadmium stress conditions in a tissue-associated manner.

A functional genomics catalogue of activated transcription factors during pathogenesis of pneumococcal disease

Background

Streptococcus pneumoniae (the pneumococcus) is the world’s foremost microbial pathogen, killing more people each year than HIV, TB or malaria. The capacity to penetrate deeper host tissues contributes substantially to the ability of this organism to cause disease. Here we investigated, for the first time, functional genomics modulation of 3 pneumococcal strains (serotype 2 [D39], serotype 4 [WCH43] and serotype 6A [WCH16]) during transition from the nasopharynx to lungs to blood and to brain of mice at both promoter and domain activation levels.

Results

We found 7 highly activated transcription factors (TFs) [argR, codY, hup, rpoD, rr02, scrR and smrC] capable of binding to a large number of up-regulated genes, potentially constituting the regulatory backbone of pneumococcal pathogenesis. Strain D39 showed a distinct profile in employing a large number of TFs during blood infection. Interestingly, the same highly activated TFs used by D39 in blood are also used by WCH16 and WCH43 during brain infection. This indicates that different pneumococcal strains might activate a similar set of TFs and regulatory elements depending on the final site of infection. Hierarchical clustering analysis showed that all the highly activated TFs, except rpoD, clustered together with a high level of similarity in all 3 strains, which might suggest redundancy in the regulatory roles of these TFs during infection. Discriminant function analysis of the TFs in various niches highlights differential regulatory backgrounds of the 3 strains, and pathogenesis data confirms codY as the most significant predictor discriminating between these strains in various niches, particularly in the blood. Moreover, the predicted TF and domain activation profiles of the 3 strains correspond with their distinct pathogenicity characteristics.

Conclusions

Our findings suggest that the pneumococcus changes the short binding sites in the promoter regions of genes in a niche-specific manner to enhance its ability to disseminate from one host niche to another. This study provides a framework for an improved understanding of the dynamics of pneumococcal pathogenesis, and opens a new avenue into similar investigations in other pathogenic bacteria.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-769) contains supplementary material, which is available to authorized users.

Comparison of hematopoietic cancer stem cells with normal stem cells leads to discovery of novel differentially expressed SSRs

Tandem repeat expansion in the transcriptomics level has been considered as one of the underlying causes of different cancers. Cancer stem cells are a small portion of cancer cells within the main neoplasm and can remain alive during chemotherapy and re-induce tumor growth. The EST-SSR background of cancer stem cells and possible roles of expressed SSRs in altering normal stem cells to cancer ones have not been investigated yet. Here, SSR distributions in hematopoietic normal and cancer stem cells were compared based on the expressed EST-SSR. One hundred eighty nine and 223 EST-SSRs were identified in cancer and normal stem cells, respectively. The EST-SSR expression pattern was significantly different between normal and cancer stem cells. The frequencies of AC/GT and TA/TA EST-SSRs were about 10% higher in cancer than normal stem cells. Remarkably, the number of triplets in cancer stem cells was 1.5 times higher than that in normal stem cells. GAT EST-SSR was frequent in cancer stem cells, but, conversely, normal stem cells did not express GAT EST-SSR. We suggest this EST-SSR as a novel triplet in cancer stem cell induction. Translating EST-SSRs to amino acids demonstrated that Asp and Ile were more abundant in cancer stem cells compared to normal stem cells. Finally, Gene Ontology (GO) enrichment analysis was carried out on genes containing triplet SSRs and showed that SSRs intentionally visit some specific GO classes. Interestingly, a NF-kappa (nuclear factor-kB) binding transcription factor was significantly hit by SSR instability which is a hallmark for leukemia stem cells. NF-kappa is an over represented transcription factor during cancer progression. It seems that there is a crosstalk between the NF-kB transcription factor and expressed GAT tandem repeat which negatively regulate apoptosis. In addition to better understanding of tumorigenesis, the findings of this study offer new DNA markers for diagnostic purposes and identifying at risk populations. In addition, a new approach for gene discovery in cancer by target analysis of differentially expressed EST-SSRs between cancer and normal stem cells is presented here.