Correction: Identification of Candidate Adherent-Invasive E. coli Signature Transcripts by Genomic/Transcriptomic Analysis

The following information is missing from the Funding section: Dr. Tarr is also supported by the National Institutes of Health [5P30 DK052574 (DDRCC)]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Copyright: © 2015 Zhang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


Introduction
Crohn's disease (CD) is a form of inflammatory bowel disease (IBD) that is characterized by skip lesions of transmural inflammation, and can occur at multiple sites in the digestive tract.Inflammation can be found anywhere in the gastrointestinal tract from the mouth to the anus, but in most (60-80%) CD patients, the distal small intestine is frequently involved [1,2].Factors implicated in the pathogenesis of IBD include host genetic predisposition, and continual activation of the mucosal immune system by luminal bacteria and their products [3,4].From 16S ribosomal RNA gene sequence data, several laboratories have demonstrated imbalances in the gut microbial composition of CD patients, particularly those with ileal involvement when compared to unaffected individuals [5][6][7][8][9][10][11][12][13][14][15][16].A consistent feature is a reduction in the relative frequency of Faecalibacterium prausnitzii [8] and an increase in Proteobacteria, particularly Escherichia coli [5].A greater relative abundance of E. coli has been associated with CD, and particularly in active disease compared to patients in remission [17].Mucosa-associated E. coli in particular are more abundant in CD [18] and in several small studies were isolated from inflamed tissue that include areas with ulcers and granulomas [19,20].In addition E. coli from the neoterminal ileum in post-surgical CD patients are linked to early recurrence of the disease [2].
Adherent invasive E. coli (AIEC) are considered to be pathobionts [21][22][23] and are isolated from the intestinal mucosa in humans with a higher prevalence in CD patients than in healthy subjects [2,24,25].The AIEC phenotype requires adherence and invasion of intestinal epithelial cells and survival and replication within macrophages [26,27].Only a few commensal E. coli have been tested for this phenotype [28].Using these methods, AIEC strains are detected in 22-52% of ileal CD patients and in 6-18% of non-IBD subjects [2,18,[29][30][31].However, these studies differ with respect to the number of biopsies analyzed, the anatomical location of the biopsies, and disease activity.
The design of a culture independent assay is hindered by the fact that although AIEC usually belong to the B2 or D groups, they are phylogenetically heterogeneous [18,32].Jensen et al [33] reported a quantitative real-time PCR (RT-qPCR) to determine the proportion of E. coli LF82 in DNA from human intestinal biopsies using spiked samples, but has not reported the results of this assay using clinical samples.Furthermore the genomic target of this assay, the pMT1-like plasmid, is not conserved among AIEC.Dogan et al, [34] reported that genes encoding processes responsible for propanediol utilization (pdu operon) and iron acquisition (yersiniabactin, chu operon) are overrepresented in human and dog AIEC genomes and might represent AIEC virulence factors.
To gain insight into biological pathways that contribute to AIEC pathogenicity we conducted a comparative transcriptomic analysis of the reference AIEC strain LF82 and the noninvasive commensal strain HS, grown in pure cultures.Furthermore, the genomic sequences of 11 non-invasive E. coli strains, including MG1655 [35] and HS [36], and a panel of 13 AIEC strains [34,[37][38][39][40][41] were compared to identify coding regions that could potentially serve as AIEC probes.Five of these gene targets and the previously described gene pduC, were tested by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) following extraction of RNA from fecal and ileal biopsy samples from 53 patients with and without IBD.

Bacterial RNA isolation, sequencing and alignment to genomes
The reference AIEC strain LF82, originally isolated by Dr. Darfeuille-Michaud, was provided as a gift by Dr. Phillip Sherman (University of Toronto) and its identity was confirmed by multi-locus sequence typing [43].The non-invasive HS strain was purchased from American Type Culture Collection (ATCC 700891).Triplicate Luria broth cultures (37°C) of LF82 and HS were grown with continuous shaking for 2 hours (exponential phase) and 24 h without shaking (stationary phase).Total RNA was extracted from the cells using the RiboPure Bacteria kit (Life Technologies Corp. Carlsbad, CA), following the manufacturer's protocol.The average RNA Integrity Number (RIN) over all samples was 7. Two micrograms of RNA was depleted of ribosomal RNA using the RiboMinus Transcriptome Isolation Kit (Life Technologies Corp. Carlsbad, CA).These samples were then used as a template for strand-specific cDNA synthesis and subjected to single-end 150 bp Illumina sequencing.The RNA-Seq libraries were prepared and sequenced at the New York Genome Center (NYGC).Raw sequences were filtered to remove human sequence contamination, remove short reads (< 50 bp), depleted of duplicate reads, and quality trimmed using Trimmomatic (v 0.32) [44].rRNA sequences were identified and culled using SortMe RNA (v1.9)[45].Raw sequence reads for LF82 and HS were mapped to NCBI reference genomes NC_011993 and NC_009800, respectively [37] using the Burroughs Wheeler aligner (BWA) [46].Counts for each annotated genomic loci were determined by HTseq-count (version 0.6.1)[47].The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE69020.
Differentially Expressed Genes (DEGs) in LF82 compared to HS Two DEG algorithms were employed, edgeR [48] and DESeq [49].The raw counts produced by HTseq-count provided the input variables for the DESeq and edgeR packages.DEGs were defined as 2 fold change and FDR < 0.05 and LF82 and HS transcripts were compared at 2h or 24h, independently.DEGs resulting from edgeR were the input variables for knowledge based biological functions using the Gene Ontology (GO) plugin BiNGO [50] and the custom ontology and annotation files found on the Gene Ontology website [51,52].DEGs resulting from DESeq were the input variables for knowledge based pathways/modules defined either by the Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/)[53] or a set of modules obtained through clustering a network of high quality functional interactions predicted for E. coli [54].The up-regulated and down-regulated output from DESeq for each time point were entered to identify the perturbed pathways regardless of the overall polarity.

Ethics Statement
This study was approved by the Institutional Review Board (IRB) at Stony Brook University Hospital.Pediatric (age 7 years) and adult patients are recruited in a consecutive fashion by the Stony Brook Digestive Diseases Research Tissue Procurement Facility and provide verbal and written consent for chart abstraction, blood, stool, tissue biopsies and/or surgical waste collection with analysis for research purposes and for their information to be stored in the hospital database.For children between 7-17 years old participating in this study, both oral and written parent/legal guardian permission and a separate oral and written assent from the child was obtained.The IRB at Stony Brook University Hospital approved this consent procedure.

Enrollment of patients and collection of samples
After receiving IRB approval, participants previously scheduled to undergo colonoscopy or intestinal resection, were identified and consented.Pediatric (ages 7 years) and adult patients were recruited in a consecutive fashion by the Stony Brook Digestive Diseases Research Tissue Procurement Facility.The period of enrollment was between March 2011 and June 2014.
Patients with a confirmed diagnosis of IBD were phenotyped based on endoscopic and radiographic studies as previously described [55].Tissue specimens were collected and immediately placed into RNAlater (Life Technologies, Carlsbad, CA).

DNA isolation from bacteria
Nine bacterial strains were processed for DNA isolation: LF82, MG1655, HS, and 6 MS AIEC strains.Following overnight culture, a single colony of each bacterial strain was placed in 5 ml of tryptic soy broth and incubated overnight at 37°C with shaking.Total bacterial DNA was extracted using the QIAamp DNA Mini Kit and according to the manufacturer's protocol and stored at -20°C until batch analysis.

PCR and electrophoresis
The forward and reverse primers for the Cas genes (strains LF82_088, LF82_091, LF82_092 and LF82_093) were designed using the NCBI primer designing tool Primer-BLAST [56].The E. coli 16S rRNA forward and reverse primers were previously validated [57].The predicted PCR products were 340 bp for E. coli 16S rRNA, 107 bp for LF82_088, 109 bp for LF82_091, 97 bp for LF82_092, and 125 bp for LF82_093.Amplification was performed in a 15 μL reaction volume and consisting of 1.5 μL 10X PCR buffer (Qiagen), 3 μL Q solution, nuclease free water, 0.5 μM forward and reverse primers, 0.1 uL Qiagen Taq DNA polymerase, and 1μL template.PCR was performed using an Eppendorf Mastercycler EPGradient S. The following thermal cycling conditions were used: 5 min at 94°C and 36 cycles of amplification consisting of 30 seconds at 95°C, 30 seconds at 56°C, and 1 min at 72°C, with 5 min at 72°C for the final extension.PCR product bands were analyzed after electrophoresis in a 1% agarose gel in 1X TBE containing ethidium bromide and digital imaging using The ChemiDoc MP system (Biorad, Hercules, CA).

RNA isolation from stool and bacteria
Total bacterial RNA was extracted from each stool sample using a fecal RNA isolation kit (Zymo Research Corporation, Irvine, CA) according to the manufacturer's protocol.RNA from strains LF82, MG1655, and HS was extracted using the same kit after culture for 2 and 24 hours.RNA was archived at -80 0 C until batch analysis.

RNA isolation from ileal biopsies
Fresh frozen ileal biopsies were homogenized individually in 2 ml of Trizol solution (Life Technologies) with the PowerGen125 homogenizer (Fisher Scientific) and 1 ml aliquots placed into 1.5 mL microcentrifuge tubes.RNA was subsequently extracted using phenol/chloroform extraction methods as previously described [58].The RNA was reconstituted in 50ul of RNA Storing Solution (Life Technologies) and stored at -80 0 C until batch analysis.

Reverse transcription quantitative polymerase chain reaction (RT-qPCR) of E. coli transcripts
For cDNA production, 500 nanograms of RNA was added to a 20 μL reaction using the Super-Script VILO cDNA Synthesis Kit (Life Technologies, Carlsbad, CA).Quantitative PCR was conducted in triplicate on 1:2 dilutions of cDNA from fecal samples and 1:2, 1:4 and 1:8 dilutions of cDNA from pure E. coli cultures and using 1 μL volumes.Amplification was performed in a 20 μL reaction volume and consisting of 10 μl of 2x SYBR Green Master Mix, 1 μl each of 10uM forward and reverse primers, 1 μL of cDNA, and 7 μL of nuclease free water.The thermal cycling conditions were: 10 min at 95°C and 40 cycles of amplification consisting of 30 seconds at 95°C and 60 seconds at 60°C using a Mastercycler EPGradient S (Eppendorf).Primers included Total bacteria and E. coli 16S rRNA forward and reverse primers as previously validated [57] and the pduC gene as previously described [34].Primers were designed for 5 candidate genes LF82_088, LF82_091, LF82_092, LF82_093, and LF82_095, using an online primer design tool [56].The sequences of all primers are listed in S2 Table.

Statistical analysis
All analyses were performed using the GraphPad Prism 5 software suite (GraphPad, San Diego, CA).For each RT-qPCR assay, the average cycle threshold (Ct) of Functional profiling of genes was accomplished using the Gene Ontology (GO) plugin BiNGO [50] and the custom ontology and annotation files on the Gene Ontology website (http://www.geneontology.org).This analysis revealed that multiple functional categories have overlapping datasets as shown in Tables A-D in S2 File.Examples include "siderophore metabolic process like enterobactin", which are up-regulated at both time points, and"glycogen metabolic process" and "oxidation-reduction process", which are down regulated at both time points (Table 3).Analysis using alternative pathways/modules gene sets [53,54] facilitated visualization of patterns of gene expression against a very complex background.For example, the functional category chemotaxis is up-regulated (FDR = 0.008) in LF82 at 2h, but bacterialtype flagellum-dependent cell motility is down regulated (FDR = 1.4 x 10 −6 ) at 24h.However as shown in Fig 1, the polarity of the DEGs are preserved at both time points.These networkbased results draw attention to modules that do not overlap with the GO categories (e.g.modules 24 and 79 in Fig 1).

Selection of candidate AIEC signature transcripts
To identify candidate AIEC signature transcripts, we took a subtractive approach to identify coding DNA sequences that were present in the genome of the reference AIEC strain LF82 but not homologous to sequences in 11 non-invasive E. coli strains.In addition to the HS and MG1655 strains, we included 9 strains from patients with and without IBD and phenotyped with respect to their inability to invade epithelial cells and survive within macrophages.
To test the in silico results and validate the PCR primers we amplified DNA for each of the 4 candidate genes (S1 Table ).Agarose gel electrophoresis of PCR reactions verified amplification products of the expected sizes (see methods) for candidate genes LF82_091, LF82_092, LF82_093 and LF82_088 in strains LF82, MS145-7, and MS57-2 (S2 Table ).All other strains, including MG1655, HS and the 4 MS AIEC strains without homologous Cas genes, exhibited no PCR amplification with these primers.All samples produced the expected band at 340 base pairs for the E. coli 16S rRNA gene product (S1 & S2 Figs).

Screening Candidate Gene Transcripts in Human Clinical Specimens
RNA was isolated from fecal samples collected from 53 individuals at Stony Brook University.Within this collection, 43 (81.1%)stool samples were acquired from children (Table 3).Twenty-two were IBD patients and 31 individuals were non-IBD controls.Non-IBD patients included subjects with functional GI disorders, Celiac disease, lactose intolerance and one patient with juvenile polyps.The number of male patients was significantly higher in both IBD cohorts compared to controls, p = 0.009 and 0.029 for CD and UC respectively.CD patients were significantly older (p = 0.024).The median ages for CD, UC/IC and controls were 20, 16 and 15 years, respectively.The IBD patients included 14 patients with CD, 6 patients with UC and 2 with indeterminate colitis (IC).Three of the CD patients were diagnosed at enrollment.Parallel ileal biopsies were available for 10 CD patients, 3 UC patients and 23 non-IBD controls.Table 5 displays the characteristics of all subjects.For IBD patients, age of diagnosis, disease location and disease behavior (CD) are as defined by the Montreal classification [60].Also included are disease duration, body mass index (BMI), smoking, surgical management of IBD, and IBD medications.
To compare the relative abundance of E. coli between clinical specimens, we performed RT-qPCR with E. coli-specific 16S rRNA gene primers and normalized results to total bacterial 16S rRNA gene expression (Tables 6 & 7).The median ΔCT values (Total-E.coli Ct) among CD, UC/IC and non-IBD fecal samples were -14.40, -7.14, and -13.56, respectively.There was no statistically significant difference in E. coli abundance compared to non-IBD controls.Among ileal biopsy specimens, the mean ΔCT values for CD, UC/IC and non-IBD samples were -9.94, -10.50, and -11.84, respectively.There was no statistically significant elevation in E. coli abundance in IBD specimens compared to controls.
The threshold of detection of transcripts corresponding to excisionase (LF82_095), pduC (LF82_332) and four Cas homologous genes (LF82_088, LF82_091, LF82_092, and LF82_093 was set at Ct 35.The negative Ct values ranged between 39 and 40.A higher proportion of CD fecal (Table 6) and ileal (Table 7) cDNA samples were positive for LF82_091 and  LF82_092 transcripts than non-IBD fecal and ileal RNA samples (p <0.05)A higher proportion of CD fecal samples were positive for LF82_088 in CD vs. non-IBD samples and a higher proportion of CD ileal samples were positive for LF82_093 AND LF82_095 in CD vs. non-IBD.
The median ΔCT values (Total-E.coli Ct) among CD, UC/IC and non-IBD fecal samples were -14.40, -7.14, and -13.56, respectively.There was no statistically significant difference in E. coli abundance when compared to non-IBD controls.Among ileal biopsy specimens, the mean ΔCT values for CD, UC/IC and non-IBD samples were -9.94, -10.50, and -11.84, respectively.There was no statistically significant elevation in E. coli abundance in IBD specimens compared to controls.

Discussion
Although a higher proportion of CD patients harbor AIEC, such organisms can also be recovered from non-IBD patients.Conversely, NIEC strains are recovered from IBD patients (Table 1).The pathogenic potential of AIEC may vary depending on host susceptibility.Host factors such as IBD risk alleles and Paneth cell function have been linked to alterations in ileal mucosa-associated microbial composition and the Escherichia/Shigella genus [12, 14, 57, 61,  62].In-vitro analysis has not been performed for many human commensal E. coli strains.In this study the complete genomes for 13 AIEC and 11NIEC, all with prior in-vitro phenotypic analysis were compared.
Multiple studies have demonstrated that CD patients, particularly those with ileal disease, have altered intestinal microbial biodiversity and composition.Because most of these studies are based on 16S rRNA sequence analysis, they do not address alterations in microbial function, or in subgroups within identified species.Shotgun bacterial DNA metagenomics and bacterial metatranscriptomics measure alterations in microbial function more directly than does 16S rRNA sequence analysis.The advantage of bacterial transcriptomic data over shotgun metagenomics data is that the former provides information on which bacterial genes are actually transcribed.In this study we compared the transcriptomes of a reference AIEC strain, LF82 to a control strain HS to identify genes associated with the AIEC phenotype.We selected HS as the control strain which was previously demonstrated to be non-invasive [28].
A comparative analysis of genes shared between the LF82 and HS genomes indicated that many of the DEG had a relatively low fold change (~2-4 fold) making them less suited for clinical assays.Up-regulated genes in LF82 are involved in many key pathways including iron metabolism, supporting the recent report that AIEC strains are enriched for genes involved in iron utilization [37], a feature of many B2 phylotype members.Comparison of the transcriptional profiles revealed a significant effect of growth conditions (see Fig 1).We identified six genes with no detectable expression in HS (Table 2) at both growth conditions.Four of the genes code for identical proteins in the enteropathogenic bacteria Salmonella and Shigella.Further characterization of these proteins in AIEC and non-invasive E. coli strains is necessary to determine if they are a component of the AIEC phenotype.
In the comparative analysis of RNA-seq data, 747 CDS that mapped to the LF82 genome did not share homology with CDS in HS (S1 Table ).We extended our comparative analysis to 13 E. coli strains with the AIEC pathotype and 11 NIEC (Table 1).Using a subtractive genomics approach, we found that the 166 CDS present only in LF82 were not homologous in all 11 NIEC (S1 Table ).However, none of the 166 CDS were present in the 13 AIEC strains surveyed.This observation supports the concept that the AIEC pathovar is formed by a heterogeneous collection of serogroups and serotypes.As shown in Table 3, AIEC genomes are enriched in genes belonging to the pdu operon, the ibe operon, and the type VI secretion system [34,37,38].The pdu operon is a component of a metabolic pathway required for fucose utilization [63], and is present in enterpathogenic bacteria and offers a competitive advantage for energy production under anaerobic conditions [34,63].The ibeA gene (invasion of brain endothelium) encodes an invasion protein found in several extraintestinal pathogenic E. coli (ExPEC) strains [64].This gene may also play a role in E. coli resistance to H 2 O s stress [65].IbeA is a necessary component for invasion of IECs and absence or mutation of this gene limits survival of AIEC within macrophage [66].The type VI secretion system has been implicated in targeting other bacterial and eukaryotic cells [67].We found homologous CDS for chuA and yersiniabactin, in 6 of 11 NIEC.These iron uptake genes are enriched among AIEC strains [37] and other pathogenic E. coli including ExPEC and EHEC.However, it remains to be determined whether these genes are expressed in the noninvasive strains.This analysis is limited by the fact that growth in pure cultures represents a very different environment than within the human intestine, and thus does not take into consideration complex microbe-microbe and host-microbe interactions.In addition, our subtractive genomics approach was limited to CDS expressed in the reference AIEC strain LF82.NRG857C has a genome that is highly similar to LF82 however, CDS in NRG857C but absent in LF82 were present in as many as six of the 13 other AIEC strains.Additional CDS that were homologous among three or more AIEC except LF82 and absent in the 11 NIEC are listed in S3 Table.Nonetheless, the results of this analysis provide a useful baseline repertoire of E. coli transcriptional patterns that may aid in the analysis of complex patient based metatranscriptomic data.
Among the 166 CDS mapping to the LF82 genome, we identified four potential signature transcripts belonging to CRISPR-associated (Cas) genes.These genes map to a region of the LF82 genome that is highly specific [37] and in our analysis these CDS were conserved in 4 of 6 AIEC strains.We did not find homologous CDS in DEC, although they are homologous to CDS in three ExPEC.Among the strains with these specific Cas genes, four of the AIEC strains and the three ExPEC are of the B2 phylotype.AIEC of the B2 phylotype are described to be among the most abundant and the most virulent [68].CRISPR-Cas forms the adaptive immunity system [69][70][71].Bacterial strains express Cas proteins that recognize foreign genetic elements in plasmids and phages and insert fragments of the exogenous DNA into their own genomes.Most E. coli harbor CRISPR-Cas systems that belong to subtype I-E [72].LF82 has the I-F system which has 3 CRISPR arrays and an operon of 6 cas-F genes (cas6f, csy3, csy2, csy1, cas2, cas3, and cas1) [72].This system is also found in Yersinia pestis an enterotoxigenic E. coli (strain B7A) and a subset of B2 phylotype E. coli [72].Toro et al, [73] examined the relationship between CRISPR-Cas systems and virulence in Shiga toxin-producing E. coli (STEC) and observed conservation of CRIPR spacer contents among strains of the same serotype and that the highly virulent STEC strains had fewer spacers within CRISPR arrays.Two other groups have recently identified CRISPR-Cas gene loci for the development of serotype-specific PCR assays of STEC [74,75] and Salmonella enterica serotypes Typhi and Paratyphi A [76].
We analyzed 53 fecal samples (Table 6) using 4 Cas gene assays and 35.7% of CD compared to 6.2% of non-IBD control samples (p = 0.02) revealed positive assays for 3 of the 4 assays.Using the pduC primers described in Dogan et al [34], expression of the pduC gene was detected in 21.4% CD compared to 25% of non-IBD controls (p = 1.0).For the excisionase gene 64% of CD compared to 44% of non-IBD control samples (p = 0.34) had positive assays.We also analyzed 38 parallel ileal biopsy samples (Table 7) and 50% of CD compared to 13% of non-IBD control samples (p = 0.04) had positive assays.Expression of the pduC gene was detected in 17% CD compared to 4% of non-IBD controls (p = 0.27).For the excisionase gene 33% of CD compared to 0% of non-IBD control samples (p = 0.0095) had positive assays.All 4 excisionase positive samples were correspondingly positive for Cas genes.The p-values for the Cas assays did not reach significance after applying the Bonferroni correction for multiple comparisons (p <0.01).Nevertheless we observed a similar trend in fecal and/or ileal biopsies for all four of the Cas genes tested.Our data suggests the Cas genes may serve as promising AIEC biomarkers; this will need to be confirmed in a larger set of patient samples.We did not detect a significant difference in E. coli 16S rRNA gene expression (ΔCT) relative to total bacteria in cases compared to non-IBD controls.
Altogether our sample sizes were small and pduC expression was less discriminating for AIEC infected samples.However, it may be a useful target for therapeutic intervention as previously described.It is also possible that other pdu operon genes are more specific and could serve as better targets.Our study is consistent with other reports that no single gene is able to distinguish AIEC from NIEC.Furthermore, it remains to be demonstrated whether any candidate AIEC signature transcripts with utility as a microbial biomarker, has a functional role in pathogenicity.
In summary, these results identify potential candidate AIEC signature transcripts, which may be more prevalent among CD patients than non-IBD patients and serve as proof of principle for our comparative genomic/transcriptomic analysis of AIEC and NIEC.

Fig 1 .
Fig 1.Functional modules differentially expressed in strain LF82 compared to non-invasive strain HS at 2h and 24h time points.As mentioned above, the differences can not be related to the duration of culture vs. the agitation.Functional modules were obtained from a high quality set of proteinprotein interaction previously defined for E. coli[36].Node size indicates fold change (log2) of differential expression (DE) based on the DESeq algorithm [33].Node color indicates direction of regulation.Flagella and chemotaxis module: (a) most chemotaxis genes are up-regulated in LF82 relative to HS at 2h (DE > 3), while at 24h (b), half of the flagella genes are down-regulated in LF82 (DE > 2).(c) Most nitrate reductase genes are down-regulated in LF82 at 2h (DE > 2).(d) At 2 hours, iron-sulphur cluster assembly genes in LF82 are down-regulated (DE > = 2.7), and enterobactin synthesis genes (involved in iron transport) are up-regulated (DE > = 2.3).doi:10.1371/journal.pone.0130902.g001 , Mesalamine, Olsalazine, Sulfasalazine c Immunomodulators: Imuran, Methotrexate d Biologics: Adalimumab, Certolizumab, Infliximab doi:10.1371/journal.pone.0130902.t005

Table 2 .
These 6 CDS are homologous in LF82 and HS but transcripts are detected only in LF82 at both 2h and 24 h.The LF82 NCBI Locus Tags and the bacterial gene names (if available) are shown.The mean normalized RPKM at 2h and 24 h is shown.

Table 3 .
Selected common up-regulated and down-regulated biological pathways in LF82 at 2h and 24h time points.For more comprehensive lists of up-regulated and down-regulated pathways at 2h and 24h please see S2A and S2B Table.The false discovery rate (FDR) is indicated for both the 2h and 24h cultures. doi:10.1371/journal.pone.0130902.t003

Table 4 .
LF82 transcripts that share homology with at least 4 other AIEC genomes but none of the 11 NIEC genomes.(See also S3 Table).Putative protein function is based on sequence homology as listed in NCBI GENE.The mean RPKM are show in the 2h and 24 h LF82 cultures.The number of AIEC strains (total of 13) with CDS sharing >85% sequence homology is listed.The genes selected for exploratory RT-qPCR analysis of patient samples are in bold.

Table 5 .
Clinical characteristics of CD, UC/IC and non-IBD patients.

Table 6 .
Fecal RT-qPCR results for candidate AIEC transcripts.The number of positive fecal stool samples are shown for each candidate AIEC transcript.Transcript is defined by LF82 locus tag and hypothetical function.Fisher's exact tests were used to compare the frequencies of positive results."*" represents P values of <0.05.The median ΔCt E. coli (range) is shown.The nonparametric Mann-Whitney test was used to compare values for IBD compared to non-IBD and CD compared to non-IBD for fecal and ileal biopsy samples, respectively.

Table 7 .
Ileal RT qPCR results for candidate AIEC transcripts.The number of positive ileal biopsy samples are shown for each candidate AIEC transcript.Transcript is defined by LF82 locus tag and hypothetical function.Fisher's exact tests were used to compare the frequencies of positive results."*" represents P values of <0.05.The median ΔCt E. coli (interquartile range) is shown.The nonparametric Mann-Whitney test was used to compare values for IBD compared to non-IBD and CD compared to non-IBD for fecal and ileal biopsy samples, respectively.