The Effect of the Stringent Response and Oxidative Stress Response on Fitness Costs of De Novo Acquisition of Antibiotic Resistance.

Resistance evolution during exposure to non-lethal levels of antibiotics is influenced by various stress responses of bacteria which are known to affect growth rate. Here, we aim to disentangle how the interplay between resistance development and associated fitness costs is affected by stress responses. We performed de novo resistance evolution of wild-type strains and single-gene knockout strains in stress response pathways using four different antibiotics. Throughout resistance development, the increase in minimum inhibitory concentration (MIC) is accompanied by a gradual decrease in growth rate, most pronounced in amoxicillin or kanamycin. By measuring biomass yield on glucose and whole-genome sequences at intermediate and final time points, we identified two patterns of how the stress responses affect the correlation between MIC and growth rate. First, single-gene knockout E. coli strains associated with reactive oxygen species (ROS) acquire resistance faster, and mutations related to antibiotic permeability and pumping out occur earlier. This increases the metabolic burden of resistant bacteria. Second, the ΔrelA knockout strain, which has reduced (p)ppGpp synthesis, is restricted in its stringent response, leading to diminished growth rates. The ROS-related mutagenesis and the stringent response increase metabolic burdens during resistance development, causing lower growth rates and higher fitness costs.

The influence of oxygen and oxidative stress on de novo acquisition of antibiotic resistance in E. coli and Lactobacillus lactis.

BACKGROUND: Bacteria can acquire resistance through DNA mutations in response to exposure to sub-lethal concentrations of antibiotics. According to the radical-based theory, reactive oxygen species (ROS), a byproduct of the respiratory pathway, and oxidative stress caused by reactive metabolic byproducts, play a role in cell death as secondary killing mechanism. In this study we address the question whether ROS also affects development of resistance, in the conditions that the cells is not killed by the antibiotic. RESULTS: To investigate whether oxygen and ROS affect de novo acquisition of antibiotic resistance, evolution of resistance due to exposure to non-lethal levels of antimicrobials was compared in E. coli wildtype and ΔoxyR strains under aerobic and anaerobic conditions. Since Lactococcus lactis (L. lactis) does not have an active electron transport chain (ETC) even in the presence of oxygen, and thus forms much less ROS, resistance development in L. lactis was used to distinguish between oxygen and ROS. The resistance acquisition in E. coli wildtype under aerobic and anaerobic conditions did not differ much. However, the aerobically grown ΔoxyR strain gained resistance faster than the wildtype or anaerobic ΔoxyR. Inducing an ETC by adding heme increased the rate at which L. lactis acquired resistance. Whole genome sequencing identified specific mutations involved in the acquisition of resistance. These mutations were specific for each antibiotic. The lexA mutation in ΔoxyR strain under aerobic conditions indicated that the SOS response was involved in resistance acquisition. CONCLUSIONS: The concept of hormesis can explain the beneficial effects of low levels of ROS and reactive metabolic byproducts, while high levels are lethal. DNA repair and mutagenesis may therefore expedite development of resistance. Taken together, the results suggest that oxygen as such barely affects resistance development. Nevertheless, non-lethal levels of ROS stimulate de novo acquisition of antibiotic resistance.

The protein modifier SUMO is critical for integrity of the Arabidopsis shoot apex at warm ambient temperatures.

SUMO is a protein modification whose conjugate levels peak during acute heat stress. We find that SUMO is also critical for plant longevity when Arabidopsis experiences a prolonged non-damaging period of only 28 degrees Celsius. Remarkably, this thermo-lethality at 28 degrees was not seen with any other mutant of the SUMO pathway tested. Autoimmunity due to low SUMO1/2 expression levels was not causal for this thermo-lethality. The role of SUMO for thermo-resilience was also distinct from its requirement for thermomorphogenesis - a growth response triggered by the same warm temperature, as only the latter response was dependent on the SUMO ligase SIZ1 as well. Thermo-resilience at 28 degrees Celsius and (acquired) thermotolerance (a response that allows plants to recover and acclimate to brief extreme temperatures) both depend on the HEAT SHOCK TRANSCRIPTION FACTOR A1 (HSFA1). Acquired thermotolerance was, however, normal in the sumo1/2 knockdown mutant. Thus, SUMO-dependent thermo-resilience is potentially controlled in a different way than the protein damage pathway that underpins thermotolerance. Close inspection of shoot apices revealed that the cell patterning and tissue integrity of the shoot apex of the SUMO1/2 knockdown mutant was lost at 28, but not 22 degrees Celsius. We thus describe a novel SUMO-dependent phenotype.

Increased comparability between RNA-Seq and microarray data by utilization of gene sets.

The field of transcriptomics uses and measures mRNA as a proxy of gene expression. There are currently two major platforms in use for quantifying mRNA, microarray and RNA-Seq. Many comparative studies have shown that their results are not always consistent. In this study we aim to find a robust method to increase comparability of both platforms enabling data analysis of merged data from both platforms. We transformed high dimensional transcriptomics data from two different platforms into a lower dimensional, and biologically relevant dataset by calculating enrichment scores based on gene set collections for all samples. We compared the similarity between data from both platforms based on the raw data and on the enrichment scores. We show that the performed data transforms the data in a biologically relevant way and filters out noise which leads to increased platform concordance. We validate the procedure using predictive models built with microarray based enrichment scores to predict subtypes of breast cancer using enrichment scores based on sequenced data. Although microarray and RNA-Seq expression levels might appear different, transforming them into biologically relevant gene set enrichment scores significantly increases their correlation, which is a step forward in data integration of the two platforms. The gene set collections were shown to contain biologically relevant gene sets. More in-depth investigation on the effect of the composition, size, and number of gene sets that are used for the transformation is suggested for future research.

The Arabidopsis SUMO E3 ligase SIZ1 mediates the temperature dependent trade-off between plant immunity and growth.

Increased ambient temperature is inhibitory to plant immunity including auto-immunity. SNC1-dependent auto-immunity is, for example, fully suppressed at 28°C. We found that the Arabidopsis sumoylation mutant siz1 displays SNC1-dependent auto-immunity at 22°C but also at 28°C, which was EDS1 dependent at both temperatures. This siz1 auto-immune phenotype provided enhanced resistance to Pseudomonas at both temperatures. Moreover, the rosette size of siz1 recovered only weakly at 28°C, while this temperature fully rescues the growth defects of other SNC1-dependent auto-immune mutants. This thermo-insensitivity of siz1 correlated with a compromised thermosensory growth response, which was independent of the immune regulators PAD4 or SNC1. Our data reveal that this high temperature induced growth response strongly depends on COP1, while SIZ1 controls the amplitude of this growth response. This latter notion is supported by transcriptomics data, i.e. SIZ1 controls the amplitude and timing of high temperature transcriptional changes including a subset of the PIF4/BZR1 gene targets. Combined our data signify that SIZ1 suppresses an SNC1-dependent resistance response at both normal and high temperatures. At the same time, SIZ1 amplifies the dark and high temperature growth response, likely via COP1 and upstream of gene regulation by PIF4 and BRZ1.

Effects of a previously selected antibiotic resistance on mutations acquired during development of a second resistance in Escherichia coli.

BACKGROUND: The effect of mutations conferring antibiotic resistance can depend on the genetic background. To determine if a previously de novo acquired antibiotic resistance influences the adaptation to a second antibiotic, antibiotic resistance was selected for by exposure to stepwise increasing sublethal levels of amoxicillin, enrofloxacin, kanamycin, or tetracycline. E. coli populations adapted to either a single or two antibiotics sequentially were characterized using whole genome population sequencing and MIC measurements. RESULTS: In a wild-type background, adaptation to any of the antibiotics resulted in the appearance of well-known mutations, as well as a number of mutated genes not known to be associated with antibiotic resistance. Development of a second resistance in a strain with an earlier acquired resistance to a different antibiotic did not always result in the appearance of all mutations associated with resistance in a wild-type background. In general, a more varied set of mutations was acquired during secondary adaptation. The ability of E. coli to maintain the first resistance during this process depended on the combination of antibiotics used. The maintenance of mutations associated with resistance to the first antibiotic did not always predict the residual MIC for that compound. CONCLUSIONS: In general, the data presented here indicate that adaptation to each antibiotic is unique and independent. The mutational trajectories available in already resistant cells appear more varied than in wild-type cells, indicating that the genetic background of E. coli influences resistance development. The observed mutations cannot always fully explain the resistance pattern observed, indicating a crucial role for adaptation on the gene expression level in de novo acquisition of antibiotic resistance.

Expression of distinct maternal and somatic 5.8S, 18S, and 28S rRNA types during zebrafish development.

There is mounting evidence that the ribosome is not a static translation machinery, but a cell-specific, adaptive system. Ribosomal variations have mostly been studied at the protein level, even though the essential transcriptional functions are primarily performed by rRNAs. At the RNA level, oocyte-specific 5S rRNAs are long known for Xenopus. Recently, we described for zebrafish a similar system in which the sole maternal-type 5S rRNA present in eggs is replaced completely during embryonic development by a somatic-type. Here, we report the discovery of an analogous system for the 45S rDNA elements: 5.8S, 18S, and 28S. The maternal-type 5.8S, 18S, and 28S rRNA sequences differ substantially from those of the somatic-type, plus the maternal-type rRNAs are also replaced by the somatic-type rRNAs during embryogenesis. We discuss the structural and functional implications of the observed sequence differences with respect to the translational functions of the 5.8S, 18S, and 28S rRNA elements. Finally, in silico evidence suggests that expansion segments (ES) in 18S rRNA, previously implicated in ribosome-mRNA interaction, may have a preference for interacting with specific mRNA genes. Taken together, our findings indicate that two distinct types of ribosomes exist in zebrafish during development, each likely conducting the translation machinery in a unique way.

Linking maternal and somatic 5S rRNA types with different sequence-specific non-LTR retrotransposons.

5S rRNA is a ribosomal core component, transcribed from many gene copies organized in genomic repeats. Some eukaryotic species have two 5S rRNA types defined by their predominant expression in oogenesis or adult tissue. Our next-generation sequencing study on zebrafish egg, embryo, and adult tissue identified maternal-type 5S rRNA that is exclusively accumulated during oogenesis, replaced throughout the embryogenesis by a somatic-type, and thus virtually absent in adult somatic tissue. The maternal-type 5S rDNA contains several thousands of gene copies on chromosome 4 in tandem repeats with small intergenic regions, whereas the somatic-type is present in only 12 gene copies on chromosome 18 with large intergenic regions. The nine-nucleotide variation between the two 5S rRNA types likely affects TFIII binding and riboprotein L5 binding, probably leading to storage of maternal-type rRNA. Remarkably, these sequence differences are located exactly at the sequence-specific target site for genome integration by the 5S rRNA-specific Mutsu retrotransposon family. Thus, we could define maternal- and somatic-type MutsuDr subfamilies. Furthermore, we identified four additional maternal-type and two new somatic-type MutsuDr subfamilies, each with their own target sequence. This target-site specificity, frequently intact maternal-type retrotransposon elements, plus specific presence of Mutsu retrotransposon RNA and piRNA in egg and adult tissue, suggest an involvement of retrotransposons in achieving the differential copy number of the two types of 5S rDNA loci.

The Conserved DNA Binding Protein WhiA Influences Chromosome Segregation in Bacillus subtilis.

The DNA binding protein WhiA is conserved in Gram-positive bacteria and is present in the genetically simple cell wall-lacking mycoplasmas. The protein shows homology to eukaryotic homing endonucleases but lacks nuclease activity. WhiA was first characterized in streptomycetes, where it regulates the expression of key differentiation genes, including the cell division gene ftsZ, which is essential for sporulation. For Bacillus subtilis, it was shown that WhiA is essential when certain cell division genes are deleted. However, in B. subtilis, WhiA is not required for sporulation, and it does not seem to function as a transcription factor, despite its DNA binding activity. The exact function of B. subtilis WhiA remains elusive. We noticed that whiA mutants show an increased space between their nucleoids, and here, we describe the results of fluorescence microscopy, genetic, and transcriptional experiments to further investigate this phenomenon. It appeared that the deletion of whiA is synthetic lethal when either the DNA replication and segregation regulator ParB or the DNA replication inhibitor YabA is absent. However, WhiA does not seem to affect replication initiation. We found that a ΔwhiA mutant is highly sensitive for DNA-damaging agents. Further tests revealed that the deletion of parAB induces the SOS response, including the cell division inhibitor YneA. When yneA was inactivated, the viability of the synthetic lethal ΔwhiA ΔparAB mutant was restored. However, the nucleoid segregation phenotype remained. These findings underline the importance of WhiA for cell division and indicate that the protein also plays a role in DNA segregation.IMPORTANCE The conserved WhiA protein family can be found in most Gram-positive bacteria, including the genetically simple cell wall-lacking mycoplasmas, and these proteins play a role in cell division. WhiA has some homology with eukaryotic homing endonucleases but lacks nuclease activity. Because of its DNA binding activity, it is assumed that the protein functions as a transcription factor, but this is not the case in the model system B. subtilis The function of this protein in B. subtilis remains unclear. We noticed that a whiA mutant has a mild chromosome segregation defect. Further studies of this phenomenon provided new support for a functional role of WhiA in cell division and indicated that the protein is required for normal chromosome segregation.

Genome rearrangements in Escherichia coli during de novo acquisition of resistance to a single antibiotic or two antibiotics successively.

BACKGROUND: The ability of bacteria to acquire resistance to antibiotics relies to a large extent on their capacity for genome modification. Prokaryotic genomes are highly plastic and can utilize horizontal gene transfer, point mutations, and gene deletions or amplifications to realize genome expansion and rearrangements. The contribution of point mutations to de novo acquisition of antibiotic resistance is well-established. In this study, the internal genome rearrangement of Escherichia coli during to de novo acquisition of antibiotic resistance was investigated using whole-genome sequencing. RESULTS: Cells were made resistant to one of the four antibiotics and subsequently to one of the three remaining. This way the initial genetic rearrangements could be documented together with the effects of an altered genetic background on subsequent development of resistance. A DNA fragment including ampC was amplified by a factor sometimes exceeding 100 as a result of exposure to amoxicillin. Excision of prophage e14 was observed in many samples with a double exposure history, but not in cells exposed to a single antibiotic, indicating that the activation of the SOS stress response alone, normally the trigger for excision, was not sufficient to cause excision of prophage e14. Partial deletion of clpS and clpA occurred in strains exposed to enrofloxacin and tetracycline. Other deletions were observed in some strains, but not in replicates with the exact same exposure history. Various insertion sequence transpositions correlated with exposure to specific antibiotics. CONCLUSIONS: Many of the genome rearrangements have not been reported before to occur during resistance development. The observed correlation between genome rearrangements and specific antibiotic pressure, as well as their presence in independent replicates indicates that these events do not occur randomly. Taken together, the observed genome rearrangements illustrate the plasticity of the E. coli genome when exposed to antibiotic stress.

Identifying small RNAs derived from maternal- and somatic-type rRNAs in zebrafish development.

rRNAs are non-coding RNAs present in all prokaryotes and eukaryotes. In eukaryotes there are four rRNAs: 18S, 5.8S, 28S, originating from a common precursor (45S), and 5S. We have recently discovered the existence of two distinct developmental types of rRNA: a maternal-type, present in eggs and a somatic-type, expressed in adult tissues. Lately, next-generation sequencing has allowed the discovery of new small-RNAs deriving from longer non-coding RNAs, including small-RNAs from rRNAs (srRNAs). Here, we systemically investigated srRNAs of maternal- or somatic-type 18S, 5.8S, 28S, with small-RNAseq from many zebrafish developmental stages. We identified new srRNAs for each rRNA. For 5.8S, we found srRNA consisting of the 5' or 3' halves, with only the latter having different sequence for the maternal- and somatic-types. For 18S, we discovered 21 nt srRNA from the 5' end of the 18S rRNA with a striking resemblance to microRNAs; as it is likely processed from a stem-loop precursor and present in human and mouse Argonaute-complexed small-RNA. For 28S, an abundant 80 nt srRNA from the 3' end of the 28S rRNA was found. The expression levels during embryogenesis of these srRNA indicate they are not generated from rRNA degradation and might have a role in the zebrafish development.

A Novel Feedback Loop That Controls Bimodal Expression of Genetic Competence.

Gene expression can be highly heterogeneous in isogenic cell populations. An extreme type of heterogeneity is the so-called bistable or bimodal expression, whereby a cell can differentiate into two alternative expression states. Stochastic fluctuations of protein levels, also referred to as noise, provide the necessary source of heterogeneity that must be amplified by specific genetic circuits in order to obtain a bimodal response. A classical model of bimodal differentiation is the activation of genetic competence in Bacillus subtilis. The competence transcription factor ComK activates transcription of its own gene, and an intricate regulatory network controls the switch to competence and ensures its reversibility. However, it is noise in ComK expression that determines which cells activate the ComK autostimulatory loop and become competent for genetic transformation. Despite its important role in bimodal gene expression, noise remains difficult to investigate due to its inherent stochastic nature. We adapted an artificial autostimulatory loop that bypasses all known ComK regulators to screen for possible factors that affect noise. This led to the identification of a novel protein Kre (YkyB) that controls the bimodal regulation of ComK. Interestingly, Kre appears to modulate the induction of ComK by affecting the stability of comK mRNA. The protein influences the expression of many genes, however, Kre is only found in bacteria that contain a ComK homologue and, importantly, kre expression itself is downregulated by ComK. The evolutionary significance of this new feedback loop for the reduction of transcriptional noise in comK expression is discussed. Our findings show the importance of mRNA stability in bimodal regulation, a factor that requires more attention when studying and modelling this non-deterministic developmental mechanism.

Improving small RNA-seq by using a synthetic spike-in set for size-range quality control together with a set for data normalization.

There is an increasing interest in complementing RNA-seq experiments with small-RNA (sRNA) expression data to obtain a comprehensive view of a transcriptome. Currently, two main experimental challenges concerning sRNA-seq exist: how to check the size distribution of isolated sRNAs, given the sensitive size-selection steps in the protocol; and how to normalize data between samples, given the low complexity of sRNA types. We here present two separate sets of synthetic RNA spike-ins for monitoring size-selection and for performing data normalization in sRNA-seq. The size-range quality control (SRQC) spike-in set, consisting of 11 oligoribonucleotides (10-70 nucleotides), was tested by intentionally altering the size-selection protocol and verified via several comparative experiments. We demonstrate that the SRQC set is useful to reproducibly track down biases in the size-selection in sRNA-seq. The external reference for data-normalization (ERDN) spike-in set, consisting of 19 oligoribonucleotides, was developed for sample-to-sample normalization in differential-expression analysis of sRNA-seq data. Testing and applying the ERDN set showed that it can reproducibly detect differential expression over a dynamic range of 2(18). Hence, biological variation in sRNA composition and content between samples is preserved while technical variation is effectively minimized. Together, both spike-in sets can significantly improve the technical reproducibility of sRNA-seq.

Dynamics of Mutations during Development of Resistance by Pseudomonas aeruginosa against Five Antibiotics.

Pseudomonas aeruginosa is an opportunistic pathogen that causes considerable morbidity and mortality, specifically during intensive care. Antibiotic-resistant variants of this organism are more difficult to treat and cause substantial extra costs compared to susceptible strains. In the laboratory, P. aeruginosa rapidly developed resistance to five medically relevant antibiotics upon exposure to stepwise increasing concentrations. At several time points during the acquisition of resistance, samples were taken for whole-genome sequencing. The increase in the MIC of ciprofloxacin was linked to specific mutations in gyrA, parC, and gyrB, appearing sequentially. In the case of tobramycin, mutations in fusA, HP02880, rplB, and capD were induced. The MICs of the beta-lactam compounds meropenem and ceftazidime and the combination of piperacillin and tazobactam correlated linearly with beta-lactamase activity but not always with individual mutations. The genes that were mutated during the development of beta-lactam resistance differed for each antibiotic. A quantitative relationship between the frequency of mutations and the increase in resistance could not be established for any of the antibiotics. When the adapted strains are grown in the absence of the antibiotic, some mutations remained and others were reversed, but this reversal did not necessarily lower the MIC. The increased MIC came at the cost of moderately reduced cellular functions or a somewhat lower growth rate. In all cases except ciprofloxacin, the increase in resistance seems to be the result of complex interactions among several cellular systems rather than individual mutations.

Absence/presence calling in microarray-based CGH experiments with non-model organisms.

Structural variations in genomes are commonly studied by (micro)array-based comparative genomic hybridization. The data analysis methods to infer copy number variation in model organisms (human, mouse) are established. In principle, the procedures are based on signal ratios between test and reference samples and the order of the probe targets in the genome. These procedures are less applicable to experiments with non-model organisms, which frequently comprise non-sequenced genomes with an unknown order of probe targets. We therefore present an additional analysis approach, which does not depend on the structural information of a reference genome, and quantifies the presence or absence of a probe target in an unknown genome. The principle is that intensity values of target probes are compared with the intensities of negative-control probes and positive-control probes from a control hybridization, to determine if a probe target is absent or present. In a test, analyzing the genome content of a known bacterial strain: Staphylococcus aureus MRSA252, this approach proved to be successful, demonstrated by receiver operating characteristic area under the curve values larger than 0.9995. We show its usability in various applications, such as comparing genome content and validating next-generation sequencing reads from eukaryotic non-model organisms.

Differences in gene expression profiles between human preimplantation embryos cultured in two different IVF culture media.

STUDY QUESTION: Is gene expression in human preimplantation embryos affected by the medium used for embryo culture in vitro during an IVF treatment? SUMMARY ANSWER: Six days of in vitro culture of human preimplantation embryos resulted in medium-dependent differences in expression level of genes involved in apoptosis, protein degradation, metabolism and cell-cycle regulation. WHAT IS KNOWN ALREADY: Several human studies have shown an effect of culture medium on embryo development, pregnancy outcome and birthweight. However, the underlying mechanisms in human embryos are still unknown. In animal models of human development, it has been demonstrated that culture of preimplantation embryos in vitro affects gene expression. In humans, it has been found that culture medium affects gene expression of cryopreserved embryos that, after thawing, were cultured in two different media for 2 more days. STUDY DESIGN, SIZE, DURATION: In a multicenter trial, women were randomly assigned to two culture medium groups [G5 and human tubal fluid (HTF)]. Data on embryonic development were collected for all embryos. In one center, embryos originating from two pronuclei (2PN) zygotes that were not selected for transfer or cryopreservation on Day 2 or 3 because of lower morphological quality, were cultured until Day 6 and used in this study, if couples consented. PARTICIPANTS/MATERIALS, SETTING, METHODS: Ten blastocysts each from the G5 and HTF study groups, matched for fertilization method, maternal age and blastocyst quality, were selected and their mRNA was isolated and amplified. Embryos were examined individually for genome-wide gene expression using Agilent microarrays and PathVisio was used to identify the pathways that showed a culture medium-dependent activity. MAIN RESULTS AND THE ROLE OF CHANCE: Expression of 951 genes differed significantly (P < 0.01) between the G5 and HTF groups. Eighteen pathways, involved in apoptosis, metabolism, protein processing and cell-cycle regulation, showed a significant overrepresentation of differentially expressed genes. The DNA replication, G1 to S cell-cycle control and oxidative phosphorylation pathways were up-regulated in the G5 group compared with the HTF group. This is in agreement with the morphological assessment of the 1527 embryos (originating from 2PN zygotes), which showed that embryos consisted of more cells on Day 2 (3.73 ± 1.30 versus 3.40 ± 1.35, P < 0.001) and Day 3 (7.00 ± 2.41 versus 5.84 ± 2.36, P < 0.001) in the G5 group when compared with the HTF group. Furthermore, the implantation rate was significantly higher in the G5 group compared with the HTF group (26.7% versus 14.7%, P = 0.002) after transfer on the second or the third day after fertilization. LIMITATIONS, REASONS FOR CAUTION: Despite careful matching of the embryos, it cannot be excluded that the differences observed between the study groups are caused by factors that we did not investigate. Extrapolation of these results to embryos used for transfer demands caution as in the present study embryos that were not selected for either embryo transfer or cryopreservation have been used for the culture experiment until Day 6. WIDER IMPLICATIONS OF THE FINDINGS: This study shows that gene expression in human preimplantation embryos is altered by the culture medium used during IVF treatment and provides insight into the biological pathways that are affected. Whether these changes in gene expression have any long-term effects on children born after IVF remains unknown. However, it is possible that early adaptations of the preimplantation embryo to its environment persist during fetal and post-natal development. STUDY FUNDING/COMPETING INTERESTS: No funding and no competing interests declared. TRIAL REGISTRATION NUMBER: Not applicable.

A novel toxicogenomics-based approach to categorize (non-)genotoxic carcinogens.

Alternative methods to detect non-genotoxic carcinogens are urgently needed, as this class of carcinogens goes undetected in the current testing strategy for carcinogenicity under REACH. A complicating factor is that non-genotoxic carcinogens act through several distinctive modes of action, which makes prediction of their carcinogenic property difficult. We have recently demonstrated that gene expression profiling in primary mouse hepatocytes is a useful approach to categorize non-genotoxic carcinogens according to their modes of action. In the current study, we improved the methods used for analysis and added mouse embryonic stem cells as a second in vitro test system, because of their features complementary to hepatocytes. Our approach involved an unsupervised analysis based on the 30 most significantly up- and down-regulated genes per chemical. Mouse embryonic stem cells and primary mouse hepatocytes were exposed to a selected set of chemicals and subsequently subjected to gene expression profiling. We focused on non-genotoxic carcinogens, but also included genotoxic carcinogens and non-carcinogens to test the robustness of this approach. Application of the optimized comparison approach resulted in improved categorization of non-genotoxic carcinogens. Mouse embryonic stem cells were a useful addition, especially for genotoxic substances, but also for detection of non-genotoxic carcinogens that went undetected by primary hepatocytes. The approach presented here is an important step forward to categorize chemicals, especially those that are carcinogenic.

A day and night difference in the response of the hepatic transcriptome to cyclophosphamide treatment.

Application of omics-based technologies is a widely used approach in research aiming to improve testing strategies for human health risk assessment. In most of these studies, however, temporal variations in gene expression caused by the circadian clock are a commonly neglected pitfall. In the present study, we investigated the impact of the circadian clock on the response of the hepatic transcriptome after exposure of mice to the chemotherapeutic agent cyclophosphamide (CP). Analysis of the data without considering clock progression revealed common responses in terms of regulated pathways between light and dark phase exposure, including DNA damage, oxidative stress, and a general immune response. The overall response, however, was stronger in mice exposed during the day. Use of time-matched controls, thereby eliminating non-CP-responsive circadian clock-controlled genes, showed that this difference in response was actually even more pronounced: CP-related responses were only identified in mice exposed during the day. Only minor differences were found in acute toxicity pathways, namely lymphocyte counts and kidney weights, indicating that gene expression is subject to time of day effects. This study is the first to highlight the impact of the circadian clock on the identification of toxic responses by omics approaches.

Staphylococcus aureus ST398 gene expression profiling during ex vivo colonization of porcine nasal epithelium.

BACKGROUND: Staphylococcus aureus is a common human and animal opportunistic pathogen. In humans nasal carriage of S. aureus is a risk factor for various infections. Methicillin-resistant S. aureus ST398 is highly prevalent in pigs in Europe and North America. The mechanism of successful pig colonization by MRSA ST398 is poorly understood. Previously, we developed a nasal colonization model of porcine nasal mucosa explants to identify molecular traits involved in nasal MRSA colonization of pigs. RESULTS: We report the analysis of changes in the transcription of MRSA ST398 strain S0462 during colonization on the explant epithelium. Major regulated genes were encoding metabolic processes and regulation of these genes may represent metabolic adaptation to nasal mucosa explants. Colonization was not accompanied by significant changes in transcripts of the main virulence associated genes or known human colonization factors. Here, we documented regulation of two genes which have potential influence on S. aureus colonization; cysteine extracellular proteinase (scpA) and von Willebrand factor-binding protein (vWbp, encoded on SaPIbov5). Colonization with isogenic-deletion strains (Δvwbp and ΔscpA) did not alter the ex vivo nasal S. aureus colonization compared to wild type. CONCLUSIONS: Our results suggest that nasal colonization with MRSA ST398 is a complex event that is accompanied with changes in bacterial gene expression regulation and metabolic adaptation.

In vivo murine hepatic microRNA and mRNA expression signatures predicting the (non-)genotoxic carcinogenic potential of chemicals.

There is a high need to improve the assessment of, especially non-genotoxic, carcinogenic features of chemicals. We therefore explored a toxicogenomics-based approach using genome-wide microRNA and mRNA expression profiles upon short-term exposure in mice. For this, wild-type mice were exposed for seven days to three different classes of chemicals, i.e., four genotoxic carcinogens (GTXC), seven non-genotoxic carcinogens (NGTXC), and five toxic non-carcinogens. Hepatic expression patterns of mRNA and microRNA transcripts were determined after exposure and used to assess the discriminative power of the in vivo transcriptome for GTXC and NGTXC. A final classifier set, discriminative for GTXC and NGTXC, was generated from the transcriptomic data using a tiered approach. This appeared to be a valid approach, since the predictive power of the final classifier set in three different classifier algorithms was very high for the original training set of chemicals. Subsequent validation in an additional set of chemicals revealed that the predictive power for GTXC remained high, in contrast to NGTXC, which appeared to be more troublesome. Our study demonstrated that the in vivo microRNA-ome has less discriminative power to correctly identify (non-)genotoxic carcinogen classes. The results generally indicate that single mRNA transcripts do have the potential to be applied in risk assessment, but that additional (genomic) strategies are necessary to correctly predict the non-genotoxic carcinogenic potential of a chemical.