Towards an advanced testing strategy for genotoxicity using image-based 2D and 3D HepG2 DNA damage response fluorescent protein reporters.

In vitro assessment of mutagenicity is an essential component in the chemical risk assessment. Given the diverse modes of action by which chemicals can induce DNA damage, it is essential that these in vitro assays are carefully evaluated for their possibilities and limitations. In this study, we used a fluorescent protein HepG2 reporter test system in combination with high content imaging. To measure induction of the DNA damage response (DDR), we used three different green fluorescent protein reporters for p53 pathway activation. These allowed for accurate quantification of p53, p21 and BTG2 (BTG anti-proliferation factor 2) protein expression and cell viability parameters at a single cell or spheroid resolution. The reporter lines were cultured as 2D monolayers and as 3D spheroids. Furthermore, liver maturity and cytochrome P450 enzyme expression were increased by culturing in an amino acid-rich (AAGLY) medium. We found that culture conditions that support a sustained proliferative state (2D culturing with normal DMEM medium) give superior sensitivity when genotoxic compounds are tested that do not require metabolisation and of which the mutagenic mode of action is dependent on replication. For compounds, which are metabolically converted to mutagenic metabolites, more differentiated HepG2 DDR reporters (e.g. 3D cultures) showed a higher sensitivity. This study stratifies how different culture methods of HepG2 DDR reporter cells can influence the sensitivity towards diverse genotoxicants and how this provides opportunities for a tiered genotoxicity testing strategy.

The in vitro assessment of the toxicity of volatile, oxidisable, redox-cycling compounds: phenols as an example.

Phenols are regarded as highly toxic chemicals. Their effects are difficult to study in in vitro systems because of their ambiguous fate (degradation, auto-oxidation and volatility). In the course of in vitro studies of a series of redox-cycling phenols, we found evidences of cross-contamination in several in vitro high-throughput test systems, in particular by trimethylbenzene-1, 4-diol/trimethylhydroquinone (TMHQ) and 2,6-di-tertbutyl-4-ethylphenol (DTBEP), and investigated in detail the physicochemical basis for such phenomenon and how to prevent it. TMHQ has fast degradation kinetics followed by significant diffusion rates of the resulting quinone to adjacent wells, other degradation products being able to air-diffuse as well. DTBEP showed lower degradation kinetics, but a higher diffusion rate. In both cases the in vitro toxicity was underestimated because of a decrease in concentration, in addition to cross-contamination to neighbouring wells. We identified four degradation products for TMHQ and five for DTBEP indicating that the current effects measured on cells are not only attributable to the parent phenolic compound. To overcome these drawbacks, we investigated in detail the physicochemical changes occurring in the course of the incubation and made use of gas-permeable and non-permeable plastic seals to prevent it. Diffusion was greatly prevented by the use of both plastic seals, as revealed by GC-MS analysis. Gas non-permeable plastic seals, reduced to a minimum compounds diffusion as well oxidation and did not affect the biological performance of cultured cells. Hence, no toxicological cross-contamination was observed in neighbouring wells, thus allowing a more reliable in vitro assessment of phenol-induced toxicity.

Fluorescent tagging of endogenous Heme oxygenase-1 in human induced pluripotent stem cells for high content imaging of oxidative stress in various differentiated lineages.

Tagging of endogenous stress response genes can provide valuable in vitro models for chemical safety assessment. Here, we present the generation and application of a fluorescent human induced pluripotent stem cell (hiPSC) reporter line for Heme oxygenase-1 (HMOX1), which is considered a sensitive and reliable biomarker for the oxidative stress response. CRISPR/Cas9 technology was used to insert an enhanced green fluorescent protein (eGFP) at the C-terminal end of the endogenous HMOX1 gene. Individual clones were selected and extensively characterized to confirm precise editing and retained stem cell properties. Bardoxolone-methyl (CDDO-Me) induced oxidative stress caused similarly increased expression of both the wild-type and eGFP-tagged HMOX1 at the mRNA and protein level. Fluorescently tagged hiPSC-derived proximal tubule-like, hepatocyte-like, cardiomyocyte-like and neuron-like progenies were treated with CDDO-Me (5.62-1000 nM) or diethyl maleate (5.62-1000 µM) for 24 h and 72 h. Multi-lineage oxidative stress responses were assessed through transcriptomics analysis, and HMOX1-eGFP reporter expression was carefully monitored using live-cell confocal imaging. We found that eGFP intensity increased in a dose-dependent manner with dynamics varying amongst lineages and stressors. Point of departure modelling further captured the specific lineage sensitivities towards oxidative stress. We anticipate that the newly developed HMOX1 hiPSC reporter will become a valuable tool in understanding and quantifying critical target organ cell-specific oxidative stress responses induced by (newly developed) chemical entities.

Comparing in vitro human liver models to in vivo human liver using RNA-Seq.

The liver plays an important role in xenobiotic metabolism and represents a primary target for toxic substances. Many different in vitro cell models have been developed in the past decades. In this study, we used RNA-sequencing (RNA-Seq) to analyze the following human in vitro liver cell models in comparison to human liver tissue: cancer-derived cell lines (HepG2, HepaRG 3D), induced pluripotent stem cell-derived hepatocyte-like cells (iPSC-HLCs), cancerous human liver-derived assays (hPCLiS, human precision cut liver slices), non-cancerous human liver-derived assays (PHH, primary human hepatocytes) and 3D liver microtissues. First, using CellNet, we analyzed whether these liver in vitro cell models were indeed classified as liver, based on their baseline expression profile and gene regulatory networks (GRN). More comprehensive analyses using non-differentially expressed genes (non-DEGs) and differential transcript usage (DTU) were applied to assess the coverage for important liver pathways. Through different analyses, we noticed that 3D liver microtissues exhibited a high similarity with in vivo liver, in terms of CellNet (C/T score: 0.98), non-DEGs (10,363) and pathway coverage (highest for 19 out of 20 liver specific pathways shown) at the beginning of the incubation period (0 h) followed by a decrease during long-term incubation for 168 and 336 h. PHH also showed a high degree of similarity with human liver tissue and allowed stable conditions for a short-term cultivation period of 24 h. Using the same metrics, HepG2 cells illustrated the lowest similarity (C/T: 0.51, non-DEGs: 5623, and pathways coverage: least for 7 out of 20) with human liver tissue. The HepG2 are widely used in hepatotoxicity studies, however, due to their lower similarity, they should be used with caution. HepaRG models, iPSC-HLCs, and hPCLiS ranged clearly behind microtissues and PHH but showed higher similarity to human liver tissue than HepG2 cells. In conclusion, this study offers a resource of RNA-Seq data of several biological replicates of human liver cell models in vitro compared to human liver tissue.

The EU-ToxRisk method documentation, data processing and chemical testing pipeline for the regulatory use of new approach methods.

Hazard assessment, based on new approach methods (NAM), requires the use of batteries of assays, where individual tests may be contributed by different laboratories. A unified strategy for such collaborative testing is presented. It details all procedures required to allow test information to be usable for integrated hazard assessment, strategic project decisions and/or for regulatory purposes. The EU-ToxRisk project developed a strategy to provide regulatorily valid data, and exemplified this using a panel of > 20 assays (with > 50 individual endpoints), each exposed to 19 well-known test compounds (e.g. rotenone, colchicine, mercury, paracetamol, rifampicine, paraquat, taxol). Examples of strategy implementation are provided for all aspects required to ensure data validity: (i) documentation of test methods in a publicly accessible database; (ii) deposition of standard operating procedures (SOP) at the European Union DB-ALM repository; (iii) test readiness scoring accoding to defined criteria; (iv) disclosure of the pipeline for data processing; (v) link of uncertainty measures and metadata to the data; (vi) definition of test chemicals, their handling and their behavior in test media; (vii) specification of the test purpose and overall evaluation plans. Moreover, data generation was exemplified by providing results from 25 reporter assays. A complete evaluation of the entire test battery will be described elsewhere. A major learning from the retrospective analysis of this large testing project was the need for thorough definitions of the above strategy aspects, ideally in form of a study pre-registration, to allow adequate interpretation of the data and to ensure overall scientific/toxicological validity.

The association of diabetes mellitus and insulin treatment with expression of insulin-related proteins in breast tumors.

BACKGROUND: The insulin receptor (INSR) and the insulin growth factor 1 receptor (IGF1R) play important roles in the etiology of both diabetes mellitus and breast cancer. We aimed to evaluate the expression of hormone and insulin-related proteins within or related to the PI3K and MAPK pathway in breast tumors of women with or without diabetes mellitus, treated with or without insulin (analogues). METHODS: Immunohistochemistry was performed on tumor tissue of 312 women with invasive breast cancer, with or without pre-existing diabetes mellitus, diagnosed in 2000-2010, who were randomly selected from a Danish breast cancer cohort. Women with diabetes were 2:1 frequency matched by year of birth and age at breast cancer diagnosis to those without diabetes. Tumor Microarrays were successfully stained for p-ER, EGFR, p-ERK1/2, p-mTOR, and IGF1R, and scored by a breast pathologist. Associations of expression of these proteins with diabetes, insulin treatment (human insulin and insulin analogues) and other diabetes medication were evaluated by multivariable logistic regression adjusting for menopause and BMI; effect modification by menopausal status, BMI, and ER status was assessed using interactions terms. RESULTS: We found no significant differences in expression of any of the proteins in breast tumors of women with (n = 211) and without diabetes (n = 101). Among women with diabetes, insulin use (n = 53) was significantly associated with higher tumor protein expression of IGF1R (OR = 2.36; 95%CI:1.02-5.52; p = 0.04) and p-mTOR (OR = 2.35; 95%CI:1.13-4.88; p = 0.02), especially among women treated with insulin analogues. Menopause seemed to modified the association between insulin and IGF1R expression (p = 0.07); the difference in IGF1R expression was only observed in tumors of premenopausal women (OR = 5.10; 95%CI:1.36-19.14; p = 0.02). We found no associations between other types of diabetes medication, such as metformin, and protein expression of the five proteins evaluated. CONCLUSIONS: In our study, breast tumors of women with pre-existing diabetes did not show an altered expression of selected PI3K/MAPK pathway-related proteins. We observed an association between insulin treatment and increased p-mTOR and IGF1R expression of breast tumors, especially in premenopausal women. This observation, if confirmed, might be clinically relevant since the use of IGF1R and mTOR inhibitors are currently investigated in clinical trials.

Insulin-like growth factor 1 receptor activation promotes mammary gland tumor development by increasing glycolysis and promoting biomass production.

BACKGROUND: The insulin-like growth factor 1 (IGF1) signaling axis plays a major role in tumorigenesis. In a previous experiment, we chronically treated mice with several agonists of the IGF1 receptor (IGF1R). We found that chronic treatment with insulin analogues with high affinity towards the IGF1R (IGF1 and X10) decreased the mammary gland tumor latency time in a p53R270H/+WAPCre mouse model. Frequent injections with insulin analogues that only mildly activated the IGF1R in vivo (glargine and insulin) did not significantly decrease the tumor latency time in this mouse model. METHODS: Here, we performed next-generation RNA sequencing (40 million, 100 bp reads) on 50 mammary gland tumors to unravel the underlying mechanisms of IGF1R-promoted tumorigenesis. Mutational profiling of the individual tumors was performed to screen for treatment-specific mutations. The transcriptomic data were used to construct a support vector machine (SVM) classifier so that the phenotypic characteristics of tumors exposed to the different insulin analogue treatments could be predicted. For translational purposes, we ran the same classifiers on transcriptomic (micro-array) data of insulin analogue-exposed human breast cancer cell lines. Genome-scale metabolic modeling was performed with iMAT. RESULTS: We found that chronic X10 and IGF1 treatment resulted in tumors with an increased and sustained proliferative and invasive transcriptomic profile. Furthermore, a Warburg-like effect with increased glycolysis was observed in tumors of the X10/IGF1 groups and, to a lesser extent, also in glargine-induced tumors. A metabolic flux analysis revealed that this enhanced glycolysis programming in X10/IGF1 tumors was associated with increased biomass production programs. Although none of the treatments induced genetic instability or enhanced mutagenesis, mutations in Ezh2 and Hras were enriched in X10/IGF1 treatment tumors. CONCLUSIONS: Overall, these data suggest that the decreased mammary gland tumor latency time caused by chronic IGF1R activation is related to modulation of tumor progression rather than increased tumor initiation.

Alternative signaling network activation through different insulin receptor family members caused by pro-mitogenic antidiabetic insulin analogues in human mammary epithelial cells.

INTRODUCTION: Insulin analogues are designed to have improved pharmacokinetic parameters compared to regular human insulin. This provides a sustained control of blood glucose levels in diabetic patients. All novel insulin analogues are tested for their mitogenic side effects, however these assays do not take into account the molecular mode of action of different insulin analogues. Insulin analogues can bind the insulin receptor and the insulin-like growth factor 1 receptor with different affinities and consequently will activate different downstream signaling pathways. METHODS: Here we used a panel of MCF7 human breast cancer cell lines that selectively express either one of the isoforms of the INSR or the IGF1R. We applied a transcriptomics approach to assess the differential transcriptional programs activated in these cells by either insulin, IGF1 or X10 treatment. RESULTS: Based on the differentially expressed genes between insulin versus IGF1 and X10 treatment, we retrieved a mitogenic classifier gene set. Validation by RT-qPCR confirmed the robustness of this gene set. The translational potential of these mitogenic classifier genes was examined in primary human mammary cells and in mammary gland tissue of mice in an in vivo model. The predictive power of the classifier genes was evaluated by testing all commercial insulin analogues in the in vitro model and defined X10 and glargine as the most potent mitogenic insulin analogues. CONCLUSIONS: We propose that these mitogenic classifier genes can be used to test the mitogenic potential of novel insulin analogues as well as other alternative molecules with an anticipated affinity for the IGF1R.

Mammary gland tumor promotion by chronic administration of IGF1 and the insulin analogue AspB10 in the p53R270H/⁺WAPCre mouse model.

INTRODUCTION: Insulin analogues are structurally modified molecules with altered pharmaco-kinetic and -dynamic properties compared to regular human insulin used by diabetic patients. While these compounds are tested for undesired mitogenic effects, an epidemiological discussion is ongoing regarding an association between insulin analogue therapy and increased cancer incidence, including breast cancer. Standard in vivo rodent carcinogenesis assays do not pick up this possible increased carcinogenic potential. METHODS: Here we studied the role of insulin analogues in breast cancer development. For this we used the human relevant mammary gland specific p53R270H/⁺WAPCre mouse model. Animals received life long repeated treatment with four different insulin (-like) molecules: normal insulin, insulin glargine, insulin X10 (AspB10) or insulin-like growth factor 1 (IGF1). RESULTS: Insulin-like molecules with strong mitogenic signaling, insulin X10 and IGF1, significantly decreased the time for tumor development. Yet, insulin glargine and normal insulin, did not significantly decrease the latency time for (mammary gland) tumor development. The majority of tumors had an epithelial to mesenchymal transition phenotype (EMT), irrespective of treatment condition. Enhanced extracellular signaling related kinase (Erk) or serine/threonine kinase (Akt) mitogenic signaling was in particular present in tumors from the insulin X10 and IGF1 treatment groups. CONCLUSIONS: These data indicate that insulin-like molecules with enhanced mitogenic signaling increase the risk of breast cancer development. Moreover, the use of a tissue specific cancer model, like the p53R270H/⁺WAPCre mouse model, is relevant to assess the intrinsic pro-carcinogenic potential of mitogenic and non-mitogenic biologicals such as insulin analogues.

Treatment with insulin (analogues) and breast cancer risk in diabetics; a systematic review and meta-analysis of in vitro, animal and human evidence.

INTRODUCTION: Several studies have suggested that anti-diabetic insulin analogue treatment might increase cancer risk. The aim of this study was to review the postulated association between insulin and insulin analogue treatment and breast cancer development, and plausible mechanisms. METHOD: A systematic literature search was performed on breast cell-line, animal and human studies using the key words 'insulin analogue' and 'breast neoplasia' in MEDLINE at PubMed, EMBASE, and ISI Web of Science databases. A quantitative and qualitative review was performed on the epidemiological data; due to a limited number of reported estimates, a meta-analysis was performed for glargine only. A comprehensive overview was composed for in vitro and animal studies. Protein and gene expression was analysed for the cell lines most frequently used in the included in vitro studies. RESULTS: In total 16 in vitro, 5 animal, 2 in vivo human and 29 epidemiological papers were included. Insulin AspB10 showed mitogenic properties in vitro and in animal studies. Glargine was the only clinically available insulin analogue for which an increased proliferative potential was found in breast cancer cell lines. However, the pooled analysis of 13 epidemiological studies did not show evidence for an association between insulin glargine treatment and an increased breast cancer risk (HR 1.04; 95 % CI 0.91-1.17; p=0.49) versus no glargine in patients with diabetes mellitus. It has to be taken into account that the number of animal studies was limited, and epidemiological studies were underpowered and suffered from methodological limitations. CONCLUSION: There is no compelling evidence that any clinically available insulin analogue (Aspart, Determir, Glargine, Glulisine or Lispro), nor human insulin increases breast cancer risk. Overall, the data suggests that insulin treatment is not involved in breast tumour initiation, but might induce breast tumour progression by up regulating mitogenic signalling pathways.

ToxProfiler: A novel human-based reporter assay for in vitro chemical safety assessment.

In vitro chemical safety assessment often relies on simple and general cytotoxicity endpoint measurements and fails to adequately predict human toxicity. To improve the in vitro chemical safety assessment, it is important to understand the underlying mechanisms of toxicity. Here we introduce ToxProfiler, a novel human-based reporter assay that quantifies the chemical-induced stress responses at a single-cell level and reveals the toxicological mode-of-action (MoA) of novel drugs and chemicals. The assay accurately measures the activation of seven major cellular stress response pathways (oxidative stress, cell cycle stress, endoplasmic reticulum stress, ion stress, protein stress, autophagy and inflammation) that play a role in the adaptive responses prior to cellular toxicity. To assess the applicability of the assay in predicting the toxicity MoA of chemicals, we tested a set of 100 chemicals with well-known in vitro and in vivo toxicological profiles. Concentration response modeling and point-of-departure estimation for each reporter protein allowed for chemical potency ranking and revealed the primary toxicological MoA of chemicals. Furthermore, the assay could effectively group chemicals based on their shared toxicity signatures and link them to specific toxicological targets, e.g. mitochondrial toxicity and genotoxicity, and different human pathologies, including liver toxicity and cardiotoxicity. Overall, ToxProfiler is a quantitative in vitro reporter assay that can accurately provide insight into the toxicological MoA of compounds, thereby assisting in the future mechanism-based safety assessment of chemicals.

Mapping the dynamics of Nrf2 antioxidant and NFκB inflammatory responses by soft electrophilic chemicals in human liver cells defines the transition from adaptive to adverse responses.

A comprehensive understanding of the dynamic activation and crosstalk between different cellular stress response pathways that drive cell adversity is crucial in chemical safety assessment. Various chemicals have electrophilic properties that drive cell injury responses in particular oxidative stress signaling and inflammatory signaling. Here we used bacterial artificial chromosome-based GFP cellular stress reporters with live cell confocal imaging, to systematically monitor the differential modulation of the dynamics of stress pathway activation by six different soft electrophiles: sulforaphane, andrographolide, diethyl maleate, CDDO-Me, ethacrynic acid and tert-butyl hydroquinone. The various soft electrophiles showed differential potency and dynamics of Nrf2 activation and nuclear translocation. These differences in Nrf2 dynamics correlated with distinct activation pattern of Nrf2 downstream targets SRNX1 and HMOX1. All soft electrophiles caused a strong dose dependent suppression of a cytokine-induced NFĸB response represented by suppression of NFĸB nuclear oscillation and inhibition of the downstream target gene activation A20 and ICAM1, which followed the potency of Nrf2 modulation but occurred at higher concentration close to saturation of Nrf2 activation. RNAi-based depletion of RelA resulted in a prolonged presence of Nrf2 in the nucleus after soft electrophile treatment; depletion of Nrf2 caused the induction of NFĸB signaling and activation of its downstream targets A20 and ICAM1. A systematic transcriptome analysis confirmed these effects by soft electrophiles on Nrf2 and NFκB signaling crosstalk in human induced-pluripotent stem cell-derived hepatocyte-like cells. Altogether our data indicate that modulation of Nrf2 by soft electrophiles may have consequences for efficient inflammatory signaling.

Systematic transcriptome-based comparison of cellular adaptive stress response activation networks in hepatic stem cell-derived progeny and primary human hepatocytes.

Various adaptive cellular stress response pathways are critical in the pathophysiology of liver disease and drug-induced liver injury. Human-induced pluripotent stem cell (hiPSC)-derived hepatocyte-like cells (HLCs) provide a promising tool to study cellular stress response pathways, but in this context there is limited insight on how HLCs compare to other in vitro liver models. Here, we systematically compared the transcriptomic profiles upon chemical activation in HLCs, hiPSC, primary human hepatocytes (PHH) and HepG2 liver cancer cells. We used targeted RNA-sequencing to map concentration transcriptional response using benchmark concentration modeling for the various stress responses in the different test systems. We found that HLCs are very sensitive towards oxidative stress and inflammation conditions as corresponding genes were activated at over 3 fold lower concentrations of the corresponding pathway inducing compounds as compared to PHH. PHH were the most sensitive model when studying UPR related effects. Due to the non-proliferative nature of PHH and HLCs, these do not pose a good/sensitive model to pick up DNA damage responses, while hiPSC and HepG2 were more sensitive in these conditions. We envision that this study contributes to a better understanding on how HLCs can contribute to the assessment of cell physiological stress response activation to predict hepatotoxic events.

System Microscopy of Stress Response Pathways in Cholestasis Research.

Exposure to oxidative radical species and cytokine-mediated inflammatory stress are established contributors to hepatocyte cell death during cholestasis. Cellular counter measures against those stressors are called adaptive stress response pathways. While in early stages of the disease adaptive stress pathways protect the hepatocytes, in later stages during prolonged stressed conditions they fail. The quantitative imaging-based assessment of cellular stress response pathways using the HepG2 BAC-GFP response reporter platform is a powerful strategy to evaluate the impact of chemical substances and gene knockdown on activation of adaptive stress response pathways, hence allowing systematic screening for positive or negative influences on cholestasis progression. This protocol allows the application of a highly versatile screening tool for a systematic evaluation of the effect of compounds having cholestasis liability and affected genes during cholestatic injury on cellular adaptive stress pathway activation. The approach involves high-throughput live-cell visualization of GFP-tagged key proteins of the oxidative stress response/Nrf2 pathway and inflammatory cytokine signaling. Quantitative image analysis of temporal responses of individual cells is followed by informatics analysis. The overall practical approaches are discussed in this chapter.

Exposure to bacterial signals does not alter pea aphids' survival upon a second challenge or investment in production of winged offspring.

Pea aphids have an obligate nutritional symbiosis with the bacteria Buchneraaphidicola and frequently also harbor one or more facultative symbionts. Aphids are also susceptible to bacterial pathogen infections, and it has been suggested that aphids have a limited immune response towards such pathogen infections compared to other, more well-studied insects. However, aphids do possess at least some of the genes known to be involved in bacterial immune responses in other insects, and immune-competent hemocytes. One possibility is that immune priming with microbial elicitors could stimulate immune protection against subsequent bacterial infections, as has been observed in several other insect systems. To address this hypothesis we challenged aphids with bacterial immune elicitors twenty-four hours prior to live bacterial pathogen infections and then compared their survival rates to aphids that were not pre-exposed to bacterial signals. Using two aphid genotypes, we found no evidence for immune protection conferred by immune priming during infections with either Serratia marcescens or with Escherichia coli. Immune priming was not altered by the presence of facultative, beneficial symbionts in the aphids. In the absence of inducible immune protection, aphids may allocate energy towards other defense traits, including production of offspring with wings that could escape deteriorating conditions. To test this, we monitored the ratio of winged to unwinged offspring produced by adult mothers of a single clone that had been exposed to bacterial immune elicitors, to live E. coli infections or to no challenge. We found no correlation between immune challenge and winged offspring production, suggesting that this mechanism of defense, which functions upon exposure to fungal pathogens, is not central to aphid responses to bacterial infections.

Escherichia coli K-12 pathogenicity in the pea aphid, Acyrthosiphon pisum, reveals reduced antibacterial defense in aphids.

To better understand the molecular basis underlying aphid immune tolerance to beneficial bacteria and immune defense to pathogenic bacteria, we characterized how the pea aphid Acyrthosiphon pisum responds to Escherichia coli K-12 infections. E. coli bacteria, usually cleared in the hemolymph of other insect species, were capable of growing exponentially and killing aphids within a few days. Red fluorescence protein expressing E. coli K-12 laboratory strain multiplied in the aphid hemolymph as well as in the digestive tract, resulting in death of infected aphids. Selected gene deletion mutants of the E. coli K-12 predicted to have reduced virulence during systemic infections showed no difference in either replication or killing rate when compared to the wild type E. coli strain. Of note, however, the XL1-Blue E. coli K-12 strain exhibited a significant lag phase before multiplying and killing aphids. This bacterial strain has recently been shown to be more sensitive to oxidative stress than other E. coli K-12 strains, revealing a potential role for reactive oxygen species-mediated defenses in the otherwise reduced aphid immune system.