Role of RelA-synthesized (p)ppGpp and ROS-induced mutagenesis in de novo acquisition of antibiotic resistance in E. coli.

The stringent response of bacteria to starvation and stress also fulfills a role in addressing the threat of antibiotics. Within this stringent response, (p)ppGpp, synthesized by RelA or SpoT, functions as a global alarmone. However, the effect of this (p)ppGpp on resistance development is poorly understood. Here, we show that knockout of relA or rpoS curtails resistance development against bactericidal antibiotics. The emergence of mutated genes associated with starvation and (p)ppGpp, among others, indicates the activation of stringent responses. The growth rate is decreased in ΔrelA-resistant strains due to the reduced ability to synthesize (p)ppGpp and the persistence of deacylated tRNA impeding protein synthesis. Sluggish cellular activity causes decreased production of reactive oxygen species (ROS), thereby reducing oxidative damage, leading to weakened DNA mismatch repair, potentially reducing the generation of mutations. These findings offer new targets for mitigating antibiotic resistance development, potentially achieved through inhibiting (p)ppGpp or ROS synthesis.

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.

Metabolic and chromosomal changes in a Bacillus subtilis whiA mutant.

IMPORTANCE: WhiA is a conserved DNA-binding protein that influences cell division in many Gram-positive bacteria and, in B. subtilis, also chromosome segregation. How WhiA works in Bacillus subtilis is unknown. Here, we tested three hypothetical mechanisms using metabolomics, fatty acid analysis, and chromosome confirmation capture experiments. This revealed that WhiA does not influence cell division and chromosome segregation by modulating either central carbon metabolism or fatty acid composition. However, the inactivation of WhiA reduces short-range chromosome interactions. These findings provide new avenues to study the molecular mechanism of WhiA in the future.

Reactive oxygen species accelerate de novo acquisition of antibiotic resistance in E. coli.

Reactive oxygen species (ROS) produced as a secondary effect of bactericidal antibiotics are hypothesized to play a role in killing bacteria. If correct, ROS may play a role in development of de novo resistance. Here we report that single-gene knockout strains with reduced ROS scavenging exhibited enhanced ROS accumulation and more rapid acquisition of resistance when exposed to sublethal levels of bactericidal antibiotics. Consistent with this observation, the ROS scavenger thiourea in the medium decelerated resistance development. Thiourea downregulated the transcriptional level of error-prone DNA polymerase and DNA glycosylase MutM, which counters the incorporation and accumulation of 8-hydroxy-2'-deoxyguanosine (8-HOdG) in the genome. The level of 8-HOdG significantly increased following incubation with bactericidal antibiotics but decreased after treatment with the ROS scavenger thiourea. These observations suggest that in E. coli sublethal levels of ROS stimulate de novo development of resistance, providing a mechanistic basis for hormetic responses induced by antibiotics.

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.

Genome sequences of 10 new carnation mottle virus variants.

Here, we report the genome sequences of 10 Carnation mottle virus variants. Six variants originated from a single proprietary carnation cultivar, and four were derived from four different proprietary cultivars. All variants showed nucleotide differences, but the last four did not show any variation at the amino acid level.

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.

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.

Extreme slow growth as alternative strategy to survive deep starvation in bacteria.

Bacteria can become dormant or form spores when they are starved for nutrients. Here, we find that non-sporulating Bacillus subtilis cells can survive deep starvation conditions for many months. During this period, cells adopt an almost coccoid shape and become tolerant to antibiotics. Unexpectedly, these cells appear to be metabolically active and show a transcriptome profile very different from that of stationary phase cells. We show that these starved cells are not dormant but are growing and dividing, albeit with a doubling time close to 4 days. Very low nutrient levels, comparable to 10,000-fold diluted lysogeny broth (LB), are sufficient to sustain this growth. This extreme slow growth, which we propose to call 'oligotrophic growth state', provides an alternative strategy for B. subtilis to endure nutrient depletion and environmental stresses. Further work is warranted to test whether this state can be found in other bacterial species to survive deep starvation conditions.

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.

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.

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.

Free SepF interferes with recruitment of late cell division proteins.

The conserved cell division protein SepF aligns polymers of FtsZ, the key cell division protein in bacteria, during synthesis of the (Fts)Z-ring at midcell, the first stage in cytokinesis. In addition, SepF acts as a membrane anchor for the Z-ring. Recently, it was shown that SepF overexpression in Mycobacterium smegmatis blocks cell division. Why this is the case is not known. Surprisingly, we found in Bacillus subtilis that SepF overproduction does not interfere with Z-ring assembly, but instead blocks assembly of late division proteins responsible for septum synthesis. Transposon mutagenesis suggested that SepF overproduction suppresses the essential WalRK two-component system, which stimulates expression of ftsZ. Indeed, it emerged that SepF overproduction impairs normal WalK localization. However, transcriptome analysis showed that the WalRK activity was in fact not reduced in SepF overexpressing cells. Further experiments indicated that SepF competes with EzrA and FtsA for binding to FtsZ, and that binding of extra SepF by FtsZ alleviates the cell division defect. This may explain why activation of WalRK in the transposon mutant, which increases ftsZ expression, counteracts the division defect. In conclusion, our data shows that an imbalance in early cell division proteins can interfere with recruitment of late cell division proteins.

Beyond the polymerase-γ theory: Production of ROS as a mode of NRTI-induced mitochondrial toxicity.

Use of some HIV-1 nucleoside reverse transcriptase inhibitors (NRTI) is associated with severe adverse events. However, the exact mechanisms behind their toxicity has not been fully understood. Mitochondrial dysfunction after chronic exposure to specific NRTIs has predominantly been assigned to mitochondrial polymerase-γ inhibition by NRTIs. However, an increasing amount of data suggests that this is not the sole mechanism. Many NRTI induced adverse events have been linked to the incurrence of oxidative stress, although the causality of events leading to reactive oxygen species (ROS) production and their role in toxicity is unclear. In this study we show that short-term effects of first generation NRTIs, which are rarely discussed in the literature, include inhibition of oxygen consumption, decreased ATP levels and increased ROS production. Collectively these events affect fitness and longevity of C. elegans through mitohormetic signalling events. Furthermore, we demonstrate that these effects can be normalized by addition of the anti-oxidant N-acetylcysteine (NAC), which suggests that ROS likely influence the onset and severity of adverse events upon drug exposure.

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.

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.

Confounding Factors in the Transcriptome Analysis of an In-Vivo Exposure Experiment.

CONFOUNDING FACTORS: In transcriptomics experimentation, confounding factors frequently exist alongside the intended experimental factors and can severely influence the outcome of a transcriptome analysis. Confounding factors are regularly discussed in methodological literature, but their actual, practical impact on the outcome and interpretation of transcriptomics experiments is, to our knowledge, not documented. For instance, in-vivo experimental factors; like Individual, Sample-Composition and Time-of-Day are potentially formidable confounding factors. To study these confounding factors, we designed an extensive in-vivo transcriptome experiment (n = 264) with UVR exposure of murine skin containing six consecutive samples from each individual mouse (n = 64). ANALYSIS APPROACH: Evaluation of the confounding factors: Sample-Composition, Time-of-Day, Handling-Stress, and Individual-Mouse resulted in the identification of many genes that were affected by them. These genes sometimes showed over 30-fold expression differences. The most prominent confounding factor was Sample-Composition caused by mouse-dependent skin composition differences, sampling variation and/or influx/efflux of mobile cells. Although we can only evaluate these effects for known cell type specifically expressed genes in our complex heterogeneous samples, it is clear that the observed variations also affect the cumulative expression levels of many other non-cell-type-specific genes. ANOVA: ANOVA analysis can only attempt to neutralize the effects of the well-defined confounding factors, such as Individual-Mouse, on the experimental factors UV-Dose and Recovery-Time. Also, by definition, ANOVA only yields reproducible gene-expression differences, but we found that these differences were very small compared to the fold changes induced by the confounding factors, questioning the biological relevance of these ANOVA-detected differences. Furthermore, it turned out that many of the differentially expressed genes found by ANOVA were also present in the gene clusters associated with the confounding factors. CONCLUSION: Hence our overall conclusion is that confounding factors have a major impact on the outcome of in-vivo transcriptomics experiments. Thus the set-up, analysis, and interpretation of such experiments should be approached with the utmost prudence.