<title>Background</title> <p>Arthropod-borne intracellular organisms that parasitize the cells of mammalian hosts must be able to manipulate a diversity of host cells to support their own growth and life cycle. Revealing how they accomplish this will illuminate not only pathogenesis but also cell biology.<italic>Anaplasma phagocytophilum </italic>(<italic>Ap</italic>) is a gram-negative obligate intracellular bacterium, the agent of human granulocytic anaplasmosis (HGA), an emerging tick-borne disease. <italic>Ap </italic>has a 1.47 million base pair genome with 1411 annotated features [##REF##16482227##1##]. Clinically, membrane bound <italic>Ap </italic>colonies, called morulae, are seen in peripheral blood neutrophils. The white-footed mouse (<italic>Peromyscus leucopus</italic>) is considered to be the primary reservoir for the <italic>Ap </italic>variant responsible for HGA, but other mammals are also susceptible [##REF##16482227##1##, ####REF##8748270##2##, ##REF##9157141##3##, ##REF##10830537##4####10830537##4##]. Ticks do not pass <italic>Ap </italic>to their offspring, but to mammals they feed upon, which transmit it back to ticks, and so the organism cycles between tick and mammalian hosts.</p> <p>HGA is a potentially severe illness with symptoms, including pancytopenia and limb edema, that suggest other cells or tissues, beside neutrophils, are infected [##REF##15113857##5##, ####REF##8980758##6##, ##REF##3558086##7####3558086##7##] In mice, <italic>Ap </italic>infects endothelial cells [##REF##16481545##8##] and human bone marrow cells support infection <italic>in vivo </italic>and <italic>in vitro </italic>[##REF##15113857##5##,##REF##9359749##9##]. The specific cells infected in ticks have not been unambiguously identified, however evidence indicates they reside within midgut and salivary gland tissues [##REF##14742506##10##, ####REF##15200875##11##, ##REF##16717118##12####16717118##12##]. Tick cell lines have been developed that support <italic>Ap </italic>replication, including ISE6, which was isolated from <italic>Ixodes scapularis</italic>, the primary vector of HGA in North America [##REF##10405394##13##]. Susceptible human cell lines include HL-60, a promyelocytic leukemia cell line that serves as a model for neutrophils, and the microvascular endothelial cell line HMEC-1 [##REF##1361507##14##]. <italic>Ap </italic>produces distinct infection phenotypes and growth kinetics in these cell lines, suggesting, along with its broad host range, that the organism adapts to each host by shifting its gene expression.</p> <p>The obligate intracellular lifestyle of <italic>Ap </italic>makes direct biochemical, genetic, and observational study approaches inherently difficult. Transformation of <italic>Ap </italic>with fluorescent reporters has recently been achieved and should improve visualization of live bacteria, and open avenues for directed genetic research [##REF##17076894##15##]. Nevertheless, methods for functional genomic analysis, for example, specific gene knockout, are still lacking. Gene transcription and expression analyses in animal models are largely impractical because <italic>Ap </italic>levels in tick and mammal tissues are too low for recovery of sufficient bacterial RNA or protein. In vitro studies have focused on characterization of the immunodominant <italic>p44/msp2 </italic>genes, which encode a large family of major surface proteins whose expression varies according to whether the organisms were derived from tick or mammalian host cells [##REF##11709787##16##]. In addition, genes encoding the type IV secretion system of <italic>Ap </italic>have been identified, transcriptionally analyzed, and described [##REF##16469062##17##,##REF##11895979##18##], but their function and regulation remain undefined. DNA microarrays have been used to measure changes in host cell gene transcription during infection, with an aim to infer the mechanisms and strategies applied by <italic>Ap </italic>[##REF##16876385##19##, ####REF##12471136##20##, ##REF##15760455##21##, ##REF##15879137##22##, ##REF##16005178##23##, ##REF##16299303##24####16299303##24##], but no microarray studies that directly measure <italic>Ap </italic>transcription have been published.</p> <p>The release of an annotated <italic>Ap </italic>genome sequence [##REF##16482227##1##], and development of maskless, photolithographic, digital light processor technology (DLP) [##UREF##0##25##] have made it feasible to characterize global transcript levels in <italic>Ap </italic>using tiling microarrays [##REF##17967045##26##,##REF##15607417##27##]. With these technologies entire genomes can be probed instead of sampling only selected sequences. The continuous data generated can be plotted in genomic order as a line graph, with transcribed genes appearing as peaks rising from a baseline of non-transcribed or intergenic sequence, and peak height corresponding to relative transcript abundance. A direct alignment of this to a parallel, annotated map of the genome can provide a visually striking and intuitive way to assess the data. Through Affymetrix (Santa Clara, CA) and NimbleGen Systems, Inc. (Madison, WI), we designed a tiling microarray for the entire genome of <italic>Ap </italic>(1.47 Mbp) and characterized <italic>Ap </italic>gene transcription in three cell lines representative of its life cycle (ISE6 tick) and pathogenesis in humans (HL-60 and HMEC-1).</p>
<title>Methods</title> <title>Cell lines, <italic>Ap </italic>strain, and growth conditions</title> <p>Sterile and <italic>Ap</italic>-infected HL-60 cells (American Type Culture Collection, Manassas, VA, USA; ATCC CCL-240) were maintained in RPMI 1640 medium supplemented with 10% fetal bovine serum (FBS) and 25 mM HEPES. Cultures infected with <italic>Ap </italic>isolate HZ were subcultured weekly by 1:50 (v/v) dilution of &gt; 90% infected cells into sterile HL-60 cultures [##REF##8531996##28##]. The HMEC-1 cell line was received from the Centers for Disease Control (Atlanta, GA), and both sterile and infected cells likewise cultured in RPMI 1640 medium with 10% FBS and 25 mM HEPES [##REF##15201033##29##]. Infected HMEC-1 cultures were fed daily and <italic>Ap </italic>subcultured 1:50 bi-weekly when &gt; 80% of cells were infected. HL-60 and HMEC-1 cultures were kept at 37°C in a humidified atmosphere of 5% CO<sub>2 </sub>in air. ISE6 cells were propagated in L15B300 medium with 5% tryptose phosphate broth (BD, Sparks MD, USA), 5% FBS, and 0.1% lipoprotein concentrate (MPBiomedical, Irvine CA, USA) at 34°C [##REF##10405394##13##]. <italic>Ap</italic>-infected ISE6 cultures were fed twice weekly with medium buffered to pH 7.6 using 0.25% NaHCO<sub>3 </sub>and 25 mM HEPES, and subcultured 1:50 bi-weekly [##REF##10405394##13##].</p> <p><italic>Ap </italic>strain HZ was cultured from the blood of a New York state patient by co-culture with HL-60 [Goodman et al. unpublished; [##REF##8531996##28##] ] HZ-<italic>Ap</italic>-infected HL-60 cells (passage 8) were simultaneously inoculated into the three cell lines. These infected parallel cultures were continuously subcultured and served as the source of infected cell samples for tiling array analysis. All samples from each cell line were from <italic>Ap </italic>cultures between passages 21 and 34.</p> <title>Tiling array design and manufacture</title> <p>Through consultation with Affymetrix (Santa Clara, CA), a library of 258,480 complimentary (perfect match) 25-mer oligonucleotide probes covering both DNA strands of the <italic>Ap </italic>genome (isolate HZ) [##REF##16482227##1##] was designed. Each probe overlapped its neighbor by 11 bases for a probe resolution of 14 bases, the distance from the center of one probe to the next. Probes were "hard pruned" – ridden of highly repetitive sequence elements thought to be irrelevant using an algorithm (Affymetrix) to identify, somewhat subjectively, long repeat sequences. Probes for these were not included, though probes for many "shorter" repeating sequences were. Pruned sequences can be viewed easily in the Artemis graphs. They are characterized by successive data points with the same or similar value that together produce large blunt peaks. For examples see additional file ##SUPPL##0##1## coordinates 665858–666184, 1025792–1026289, and 645698–646032. NimbleGen Systems, Inc. (Madison, WI) synthesized the oligonucleotide probes <italic>in situ </italic>using a photo-mediated, maskless process in which the synthesis of each probe is directed by a digital light processor [##UREF##0##25##].</p> <title>Isolation of RNA</title> <p><italic>Ap </italic>genomic transcription was measured in each of the three cell lines when cultures were approximately 95% infected. Typically, cells contained hundreds of bacteria (Figure ##FIG##0##1##: Microscopic images of Giemsa stained cells infected with <italic>Ap</italic>). RNA was extracted from three <italic>Ap</italic>-infected and three uninfected samples of each cell line (18 samples total). Each sample was from a separate culture and consisted of approximately 10<sup>7 </sup>infected cells or uninfected control cells. Cells were suspended by pipetting (HL-60 and ISE6) or with a cell scraper (HMEC-1) and immediately centrifuged at 300 × g for 2 minutes. The supernatant was aspirated and discarded; cell pellets were loosened by flicking and immediately dissolved in TRI REAGENT™ (Sigma, Saint Louis, MO, USA). All steps were performed at room temperature. Total RNA was then isolated according to the TRI REAGENT™ product instructions. In brief, samples in TRI REAGENT™ were extracted with chloroform and centrifuged at 12,000 × g for 15 minutes at 4°C. RNA in the aqueous, upper phase was precipitated in isopropanol, collected by centrifugation at 12,000 × g for 10 minutes at 4°C, and washed twice in cold 75% ethanol. RNA pellets were dissolved in 100 μL RNase-free water, quantified by spectrophotometry, and processed for array analysis.</p> <title>Preparation of tiling array "target"</title> <p>Total RNA from <italic>Ap</italic>-infected or sterile control cells was processed for <italic>Ap </italic>transcript measurement according to the Affymetrix "Prokaryotic Target Preparation" protocol using random priming of total RNA to synthesize a single strand of cDNA. The cDNA was recovered by column purification, fragmented with DNase I, and end labeled with biotin. These biotinylated cDNA fragment "targets" were hybridized to the "probes" contained on the tiling arrays, labeled with a streptavidin-phycoerythrin conjugate, and probe hybridization was quantified by laser scanning. The detailed protocol was as follows.</p> <title>cDNA synthesis</title> <p>In a volume of 30 μL, 10 μg of total RNA (from <italic>Ap </italic>infected or sterile control cells) was combined with random primers (25 ng/μL final concentration) (Invitrogen, Carlsbad, CA), and, in a thermocycler, incubated 10 minutes at 70°C followed by 10 minutes at 25°C, then chilled to 4°C. To this reaction mixture was added 30 μL of the following master mix: 12 μL 5× 1<sup>st </sup>Strand Buffer, 6 μL 100 mM DTT, 3 μL 10 mM dNTPs, 1.5 μL SUPERaseIn™ (20U/μL) (Ambion, Austin TX, USA), 7.5 μL SuperScript II (200 U/μL) (Invitrogen). Samples (60 μL) were incubated in a thermocycler 10 minutes at 25°C, 60 minutes at 37°C, 60 minutes at 42°C, 10 minutes at 70°C, and chilled to 4°C.</p> <title>cDNA isolation and fragmentation</title> <p>To degrade RNA, 20 μL of 1N NaOH was added to each sample, incubated at 65°C for 30 minutes, and neutralized by addition of 20 μL 1N HCl. MiniElute PCR Purification Columns (Qiagen, Valencia CA, USA) were used according to product instructions to purify cDNA from the samples. Typical cDNA yields were 3–4 μg. cDNA in 10 μL was combined with 2 μL 10× One-Phor-All Buffer (Amersham Biosciences, Piscataway, NJ), 0.6 U DNase I/μg cDNA (Amersham Biosciences), plus sufficient water for 20 μL total volume, and incubated 10 minutes at 37°C. DNase I was inactivated by heating to 98°C for 10 minutes. cDNA fragments produced were 50–200 bases in length.</p> <title>Biotinylation of 3' termini of cDNA fragments</title> <p>The GeneChip<sup>® </sup>DNA labeling kit (Affymetrix) was used as follows: 20 μL fragmented cDNA was combined with 10 μL 5× reaction buffer, 2 μL 7.5 mM GeneChip DNA labeling reagent, 2 μL terminal deoxynucleotidyl transferase, and 16 μL water and incubated at 37°C for 60 minutes. The reaction was stopped with 2 μL of 0.5 M EDTA and then frozen at -20°C until it was applied to an array.</p> <title>Tiling array hybridization and scanning</title> <p>Samples were hybridized to tiling arrays and scanned at the BioMedical Genomics Center at the University of Minnesota using the Affymetrix Fluidics Station 400. Arrays were scanned using an Affymetrix Genechip 3000 scanner according to standard Affymetrix protocols.</p> <title>Tiling array data analysis</title> <p>"Cel" files generated by the University of Minnesota's microarray facility were joined to Affymetrix BPMAP files specific to the tiling array using Affymetrix<sup>® </sup>Tiling Analysis Software (TAS). TAS generated a list of signal intensities and arranged them in order of genomic location and DNA strand. The data are available at the NCBI Gene Expression Omnibus (GEO) database (study #GSE11487 <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/projects/geo/query/acc.cgi?acc=GSE11487"/>).</p> <p>Graphical representation of these data along with their annotations was accomplished with the JAVA based program "Artemis" <ext-link ext-link-type="uri" xlink:href="http://www.sanger.ac.uk/"/>. Using a script developed internally, the intensity plots were reformatted and imported into Artemis along with an annotation feature list <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov"/>. The resulting graphics give a visual overview of transcription as it relates to genomic organization, and provide clues to operon structure (see additional file ##SUPPL##0##1##: Artemis transcription graph of the entire, annotated <italic>Ap </italic>genome during infection of HL-60, HMEC-1, and ISE6 cell lines). The complete genome coverage provided by the overlapping probes on the tiling array translates into 90 spot intensities generated for a 1000 base open reading frame (ORF). This large number of intensities, coupled with the quality of data suggested that creating a linear graph, and measuring the area under the peaks in regions corresponding to annotated open reading frames – ORF transcription areas – would be a simple and useful method to quantify transcripts for each ORF. To compute these ORF transcription areas, the intensities were normalized via quantiles [##REF##16939796##30##] and imported into the IgorPro data analysis program (WaveMetrics Lake Oswego OR, USA) along with the ORF and structural RNA annotations available from <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov"/>. A script was written to index a trapezoidal integration algorithm of the intensity list with the start and end genomic positions indicated on the annotation. This script operation generated a list of 1411 transcription areas.</p> <title>Statistical evaluation of area differences, T values &amp; Fold change</title> <p>ORF transcription areas computed from the quantile normalized data (3 each for HL-60, HMEC-1 and ISE6) and paired 2 tail Students t-test, were performed on: HL-60 vs. ISE6, HMEC-1 vs. ISE6, and HL-60 vs. HMEC-1. ORF transcription area comparisons with p values ≤ 0.05 were considered significant for determination of the number and identity of genes transcribed. Determination of differentially expressed genes utilized the additional requirement that the mean ORF transcription area be at least twice the mean ORF transcription area of the same gene of the compared cell line.</p> <p>The number of expressed ORFs was determined by T-test comparison between the ORF transcription areas from infected cell monolayers, and those of uninfected control cell monolayers. The signal intensity of these arrays was baseline corrected using the signal intensities of twelve manually selected intergenic regions devoid of obvious signal from across the span of the genome. ORF transcription area comparisons with p values ≤ 0.05 were considered significant for determination of the number and identity of genes transcribed.</p> <title>Validation of tiling array data by quantitative reverse transcription-PCR (qRT-PCR)</title> <p>Five <italic>Ap </italic>genes with known products were assayed for relative transcript abundance by qRT-PCR. Tiling data indicated that four of the genes had differential transcription patterns between the human and tick cells: major outer membrane protein (<italic>omp-1A</italic>; <italic>APH_1359</italic>), outer membrane efflux protein (<italic>APH_1110</italic>), major surface protein 4 (<italic>msp4</italic>; <italic>APH_1240</italic>), and the 60 kDa chaperonin (<italic>APH_0240</italic>). The fifth gene, which codes for succinyl-CoA synthetase beta subunit (<italic>APH_1052</italic>), was transcribed equally in all three cell lines (Figure ##FIG##1##2##; Artemis transcription profiles for five genes chosen for assay by qRT-PCR).</p> <p>Total RNA (from portions of samples prepared for array analysis), from three separate cultures of each <italic>Ap</italic>-infected cell line (9 samples), were assayed in triplicate by qRT-PCR. To eliminate any DNA contamination, samples were DNase I treated using DNA-free™ (Ambion). DNase I was inactivated and RNA purified using RNeasy mini columns (Qiagen, Valentia, CA). mRNA from each of the five genes was reverse transcribed and amplified quantitatively with primers designed using MacVector (Cary, NC) and Netprimer (Palo Alto, CA) (see additional file ##SUPPL##1##2##: qRT-PCR primers). The primers were tested by conventional PCR on a Stratagene (La Jolla CA, USA) Robocycler with temperature gradient capability, using <italic>Ap </italic>strain HZ DNA as target. Formation of appropriate product sizes was verified and a single annealing temperature (60°C) and primer concentration (150 nM) suitable for all five primer pairs were determined, allowing RNA from each of the cell lines to be qRT-PCR-amplified together for best determination of relative transcript levels. Reverse transcription and subsequent quantitative PCR were performed on 100 ng of each RNA sample in 96-well plates using the Brilliant II SYBR Green 1-step qRT-PCR kit (Stratagene), and Stratagene's Mx3005P thermal cycler. To initiate the qRT-PCR, reverse transcription was allowed to proceed for 30 minutes at 50°C, followed by heat treatment for 10 minutes at 95°C to activate DNA polymerase and deactivate reverse transcriptase. cDNA was then amplified during 40 cycles of 30 seconds at 95°C, 1 minute at 60°C, and 1 minute at 72°C.</p>
<title>Results</title> <title>Percentage of <italic>Ap </italic>genes measured as transcribed in each cell line</title> <p>Of the 1411 annotated features [##REF##16482227##1##] in the <italic>Ap </italic>genome, 983 (69.6%) were significantly transcribed (p-value ≤ 0.05) in HL-60, 620 (43.9%) in HMEC-1, and 974 (69.0%) in ISE6, compared to negative control samples (RNA from uninfected cells).</p> <title>Differential <italic>Ap </italic>gene transcription between cell lines</title> <p>Between HL-60 and HMEC-1, 71 <italic>Ap </italic>ORFs (5%) were differentially (p-value ≤ 0.05) transcribed (see additional file ##SUPPL##2##3##: <italic>Ap-</italic>HL-60 vs. <italic>Ap</italic>-HMEC-1 differential transcription). Between HL-60 and ISE6, 585 <italic>Ap </italic>ORFs (41.5%) were differentially transcribed. Between HMEC-1 and ISE6, 304 <italic>Ap </italic>ORFs (21.5%) were differentially transcribed. Adding a fold change criterion of 2 or greater, only one <italic>Ap </italic>gene between <italic>Ap </italic>from HL-60 (<italic>Ap</italic>-HL-60) and <italic>Ap </italic>from HMEC-1 (<italic>Ap</italic>-HMEC-1) passed: <italic>APH_1342</italic>, one of the <italic>p44</italic>/<italic>msp2 </italic>paralogs. Between <italic>Ap</italic>-HL-60 and <italic>Ap </italic>from ISE6 (<italic>Ap</italic>-ISE6), 117 ORFs (8.5%), and between <italic>Ap</italic>-HMEC-1 and <italic>Ap</italic>-ISE6, 61 (4.3%) were at least 2-fold different (Table ##TAB##0##1##). The relatively low percentage of ORFs measured as transcribed in <italic>Ap</italic>-HMEC-1 (43.9%) was probably due to lower average signal intensity from those samples (850 vs. 2120 in HL-60). We determined this to be the result of suboptimal biotin labeling after using a particular batch of terminal transferase. A new aliquot of terminal transferase used in the preparation of one of the samples of <italic>Ap</italic>-HL-60 produced a particularly bright signal, resulting in a higher signal to noise ratio for the <italic>Ap-</italic>HL-60 data. Because of this, and because differential transcription was low between the human cell lines (compared to that between the human and tick cell lines), subsequent descriptions of differential transcription in <italic>Ap</italic>-ISE6 are based on comparisons to <italic>Ap-</italic>HL-60.</p> <p>Of the 117 <italic>Ap </italic>ORFs differentially transcribed (p ≤ 0.05, ≥ two-fold difference) between the HL-60 and ISE6 cells, 76 had higher levels in HL-60 and 41 had higher levels in ISE6. The 76 <italic>Ap</italic>-HL-60 ORFs comprise 35 known and 41 hypothetical proteins (54%). All but three of the ORFs that were up-regulated in <italic>Ap</italic>-ISE6 are annotated as hypothetical (93%) (see Table ##TAB##1##2##: Genes differentially transcribed between human (HL-60) and tick (ISE6) cells). By comparison, 40% of all <italic>Ap </italic>genes are annotated as hypothetical.</p> <p>The amino acid sequences derived from the 117 differentially transcribed genes were analyzed using the secretomeP CBS prediction server [##REF##16212653##31##] and the CELLO subcellular localization predictor [##REF##16752418##32##] to determine the probable cellular location of each of the gene products – periplasm, inner or outer membrane, extracellular (secreted), or cytoplasmic. While 25% of all <italic>Ap </italic>genes products are membrane associated (non cytoplasmic), 43% of the 76 genes differentially transcribed in HL-60 cells and 46% of the 41 genes differentially transcribed in ISE6 cells code for non cytoplasmic proteins. As illustrated in Table ##TAB##1##2##, the greater a gene's differential transcription, the more likely it was to encode a membrane associated protein (i.e. differentially transcribed genes were over-represented by membrane associated proteins; see Table ##TAB##1##2##: Summary of <italic>Ap</italic>-HL-60 vs. <italic>Ap</italic>-ISE6 differential gene transcription).</p> <p>As illustrated in the Artemis transcript level graphs (see additional file ##SUPPL##0##1##: materials for graphing transcript level data in Artemis), when the data are displayed as linear graphs alongside a map of the annotated genome, numerous transcription behaviors are revealed. Transcribed sequences are seen to rise from the over-all flat baseline and generally correspond well to annotated ORFs. However, there are examples of transcript signal extending beyond ORF boundaries (<italic>APH_ </italic>numbers <italic>0005, 0406, 0793, 0808, 0811, 0859, 0906</italic>, and <italic>1151</italic>), transcription apparently not associated with an ORF (coordinates 46672–46738, 944100–944549, 692299–692983, and 1306128–1306875), and transcribed unannotated ORFs (875684–876751, 1445252–1445797 and 1241148–1241727). The ORF identified between coordinates 1241148 and 1241727 is another <italic>p44/msp2 </italic>paralog, bringing the total number of <italic>p44 </italic>loci now identified to 114 (113 were originally annotated; [##REF##16482227##1##]. Peaks and plateaus of varying profile representing gene transcription are clearly discernible. Often they slope downward from 5' to 3', but sometimes they are flat (Figure ##FIG##2##3##: Examples of flat and sloped transcription peaks). There are also numerous ORFs and operons that showed no significant transcription in any of the cell lines (see additional file ##SUPPL##3##4##: genes and operons with no detected transcripts).</p> <p>Paralogs of the <italic>p44/msp2 </italic>family of outer membrane proteins form a characteristic hybridization pattern that is somewhat perplexing. Since <italic>p44 </italic>is abundantly expressed in <italic>Ap</italic>, transcripts with sequences that correspond to the conserved ends of the gene should bind to all the probes on the array that are complimentary – i.e., those of over 100 genes. Signals associated with the conserved ends of the <italic>p44 </italic>paralogs do rise sharply, while those that correspond to the hypervariable region (HVR) in between are generally near baseline. This produces a double horn shaped signature. Most paralogs are not expressed within a population of bacteria [##REF##18285495##33##] therefore those that display bridged horns – representing transcript hybridization to the HVR – are likely to be specifically transcribed. In HL-60, <italic>APH_1152 </italic>(similar to p44-47) and <italic>APH_1351 </italic>(similar to <italic>p44-35</italic>), and in HMEC-1, <italic>APH_1253 </italic>(similar to <italic>p44-39</italic>), <italic>APH_1342 </italic>(similar to <italic>p44-31</italic>), and <italic>APH_1350 </italic>(similar to <italic>p44-51</italic>) had strong signals associated with their HVRs, suggesting those paralogs were expressed. <italic>Ap</italic>-ISE6 produced no significant hybridization to any of the <italic>p44 </italic>HVRs, however along with <italic>Ap</italic>-HL-60 and <italic>Ap-</italic>HMEC-1, <italic>Ap</italic>-ISE6 produced strong signals to the conserved <italic>p44 </italic>sequences. In all three cell lines, signals to the conserved <italic>p44 </italic>sequences were greater than those from the HVRs – of the expressing paralogs noted in <italic>Ap</italic>-HL-60 and <italic>Ap</italic>-HMEC-1. In addition, this pattern of excessive hybridization to the conserved ends of the <italic>p44 </italic>ORFs, is "reflected" in the non-coding DNA strand. Probes to sequences opposite conserved <italic>p44 </italic>sense sequences are hybridized significantly in the human cell samples, and as strongly in the tick cells as the sense probes, such that the horned profile appears reflected in the opposite DNA strand. (Figure ##FIG##3##4##: <italic>p44 </italic>transcription phenomena: horns, reflecting, and HVR associated signal)</p> <p>Exceptions are <italic>p44-70, p44-71, p44-72, and p44-79</italic>, which have "conserved" ends that differ significantly from the other <italic>p44s</italic>; they produced no horns or reflections (see additional file ##SUPPL##0##1##, coordinates 680648–684696 and 1418814–1420199). Subtler reflecting was also seen in several non-<italic>p44 </italic>ORFs, such as <italic>APH_1387</italic>, which codes for outer membrane protein HGE2 [##REF##16482227##1##], and the hypothetical <italic>APH_0536 </italic>(Figure ##FIG##4##5##: Reflecting).</p> <p>Like conserved <italic>p44 </italic>sequences, repeat sequences, which are common throughout the genome, generally displayed strong signals on both DNA strands (see additional file ##SUPPL##4##5##: Repeat-sequence-based sense and anti sense signal).</p> <p>At the <italic>p44 </italic>expression locus (<italic>APH_1221</italic>) both <italic>Ap</italic>-HL-60 and <italic>Ap-</italic>HMEC-1 showed strong transcription beginning near base 1289280, just before the start of the <italic>omp-1N </italic>gene, and continuing through the <italic>p44 </italic>expression site, while <italic>Ap</italic>-ISE6 did not. The <italic>p44 </italic>"horns" seen in <italic>Ap</italic>-ISE6 within the expression locus, are likely examples of the generalized hybridization to conserved <italic>p44 </italic>sequence noted above. The <italic>tr1 </italic>gene (<italic>APH_1218</italic>) upstream of the <italic>p44 </italic>expression locus, which encodes a putative transcription regulator [##REF##17905983##34##], is well transcribed by <italic>Ap-</italic>ISE6 but not by <italic>Ap</italic>-HL-60 or <italic>Ap</italic>-HMEC-1. The DNA binding protein <italic>ApxR </italic>(<italic>APH_0515</italic>; [##REF##17905983##34##] was weakly transcribed in the human cell lines but not at all in the tick cell line (Figure ##FIG##5##6##: Artemis transcription plots of the <italic>p44 </italic>expression site, and <italic>ApxR</italic>, a putative <italic>p44 </italic>transcription regulator).</p> <p>The type IV secretion system genes identified by Hotopp et al. [##REF##16482227##1##] consistently showed little activity in any of the host cells, while <italic>sodB </italic>(<italic>APH_0371</italic>), an iron superoxide dismutase shown to be co-transcribed with components of the type IV secretion system of <italic>E. chaffeensis </italic>and <italic>Ap </italic>[##REF##11895979##18##], was moderately transcribed by <italic>Ap </italic>in all three cell lines. <italic>Ank </italic>(<italic>APH_0740</italic>) was strongly transcribed in <italic>Ap</italic>-HMEC-1, somewhat less so in <italic>Ap-</italic>HL-60, and only marginally in <italic>Ap</italic>-ISE6. This <italic>Ap </italic>gene encodes a protein that is translocated to the nucleus of infected HL-60 cells [##REF##15236641##35##,##REF##10948155##36##] and phosphorylated there within minutes [##REF##17250594##37##], presumably as an effector molecule delivered via the <italic>Ap </italic>type IV secretion system [##REF##17587335##38##]. Located between genome coordinates 1194300 and 1203600 are eight paralogs of the <italic>TrbC/VirB2 </italic>gene family (pfam04956), six of which showed measurable transcript levels either only in the tick cell line (<italic>APH_1131 – APH_1134</italic>), or the human cell lines (<italic>APH_1144 </italic>and <italic>APH_1145</italic>). The relationship by amino acid sequence of these eight paralogs is illustrated in Figure ##FIG##6##7## (Phylogenetic tree of eight <italic>virB2 </italic>paralogs by amino acid sequence), and indicates those transcribed in ISE6 are more closely related to each other than those transcribed in HL-60 and HMEC-1. Amino acid sequence alignments for the eight <italic>virB2 </italic>paralogs of <italic>Ap </italic>(see additional file ##SUPPL##5##6##) show identities that rank from a high of 93% between tick cell expressed paralogs <italic>APH_1133 </italic>and <italic>APH_1134</italic>, and a low of 22% between non-expressed <italic>APH_1136 </italic>and human cell expressed <italic>APH_1145</italic>. Multiple alignment showed higher identity and similarity between the C termini of paralogs, which contain the functional portion of the proteins.</p> <p>Two apparent tick-cell-specific operons were identified. ORFs between coordinates 1448342 and 1445170, which include locus tags <italic>APH_1386 </italic>through <italic>APH_1382</italic>, were transcribed only in the tick cell line (see additional file ##SUPPL##6##7##: Tick- and human-specific <italic>Ap </italic>operons). Locus tag <italic>APH_1380 </italic>appears to be part of the operon and as such was transcribed in the tick cell line, and, at a lower level in the human cell lines. The functions of the hypothetical proteins of these six ORFs are not known. However, a BLAST homology search produced E values of 9e-18 to 4e-9, indicating the six ORFs are related. The transcription profile around <italic>APH_1380 </italic>and sequence characteristics just up-stream, suggest that the ORF actually begins with the methionine at coordinate 1445107. In support of this, there is a ribosomal binding site at coordinate 1445120. This upstream area shows significant amino acid sequence homology with the N-termini of the other ORF members of this putative operon, also suggesting the sequence is part of that ORF. Between coordinates 1445252 and 1445797 an un-annotated ORF appears to be transcribed only in the tick cell line, and also shows significant homology to the other ORFs in this putative operon. If this is a true ORF, and the start of <italic>APH_1380 </italic>is extended to coordinate 1445107, the two putative ORFs <italic>APH_1381 </italic>and <italic>APH_1382 </italic>on the positive DNA strand may not be true ORFs, since they are situated opposite coding sequences in the operon and showed no transcription signal (see additional file ##SUPPL##6##7## panel A). The other apparent tick specific operon includes locus tags <italic>APH_0726 </italic>through <italic>APH_0720 </italic>(see additional file ##SUPPL##6##7## panel B). All but the small locus tags <italic>APH_0721 </italic>and <italic>APH_0722 </italic>were transcribed. Although these genes are also annotated as encoding hypothetical proteins, searches using SignalP [##REF##15223320##39##] and TMHMM [##REF##9783223##40##] prediction servers indicated they all have transmembrane domains. There was also a group of <italic>Ap </italic>genes transcribed only in the human cells: <italic>APH_0837</italic>, <italic>APH_0838</italic>, <italic>APH_0839</italic>, and <italic>APH_0842 </italic>(see additional file ##SUPPL##6##7## panel C). All encode hypothetical proteins and all are related by amino acid sequence, especially <italic>APH_0838</italic>, <italic>APH_0839</italic>, and <italic>APH_0842</italic>.</p> <title>qRT-PCR</title> <p>Relative transcript levels for the five selected <italic>Ap </italic>genes, within and between cell lines, confirm those indicated by the array data (Figure ##FIG##7##8##; Tiling vs. qRT-PCR graphs).</p>
<title>Discussion</title> <p>Total RNA from <italic>Ap </italic>infected human and tick cells was used to establish host cell specific <italic>Ap </italic>transcription profiles by hybridization to complementary oligonucleotides representing the entire genome of <italic>Ap </italic>on tiling arrays. The high percentages of genes measured as transcribed (69.6% in HL-60, 43.9% in HMEC-1, and 69.0% in ISE6), and the low levels of hybridization produced by the uninfected control samples, demonstrate that the method and array design produced sensitive, consistent, and specific transcription measurements. This is encouraging since efforts to fractionate or amplify RNA samples inevitably skew results. However, the culture samples analyzed were heavily infected and therefore optimal for such a direct approach. The three cell lines – HL-60 human promyelocytic, HMEC-1 human microvascular endothelial, and ISE6 tick – each produced bacteria with distinct transcription profiles, suggesting that <italic>Ap </italic>gene expression is closely dependent on the phenotype and genotype (species origin) of its host cell. The bacteria assayed were not synchronized, they were the result of 1:50 inoculations, and therefore the transcription profiles generated were an average, perhaps with a "late stage" bias, of the infection process in each cell line.</p> <p>Transcription profiles between the two human cell lines appeared similar, however with better and more consistent biotin labeling the percentage of <italic>Ap </italic>ORFs transcribed in HMEC-1 (43.9%) is predicted to be closer to that seen in HL-60 and ISE6 (~70%), and differences in transcription profiles between <italic>Ap</italic>-HL-60 and <italic>Ap</italic>-HMEC-1 would be magnified to reveal additional essential characteristics of <italic>Ap </italic>transcription in the human promyelocytic versus endothelial cells. Transcription differences between the human and tick cells were extensive; there were many genes and apparent operons transcribed in the tick cells but not in the human cells, and vice versa. The fact that the vast majority of tick cell specific transcripts are for hypothetical genes is tantalizing, and likely reflects our ignorance of the molecular patho-physiology of ticks and their associated bacteria.</p> <p>The observation that in all three cell lines some <italic>Ap </italic>genes and operons remained inactive, is either an indication that there are genetic capabilities not called for by these in vitro infection conditions – the particular intracellular environments of each cell line and the laboratory growth conditions – or the failure of this method to measure the transcription of those genes. Genes and operons that were truly silent may, among other possibilities, encode products specific to earlier stages of infection, to colonization of ticks following blood-meal uptake, or to parasitism of different hosts. Given the distinct transcription profiles produced between the human and tick cells, and the diversity of animal hosts and cell types infected within each, all are possible explanations.</p> <p>The <italic>virB2 </italic>paralogs of the type IV secretion system (T4SS) identified as differentially transcribed (6 of 8) between the human and tick cells (<italic>APH_1144 </italic>and <italic>APH_1145</italic>, and <italic>APH_1131 – APH_1134</italic>, respectively) represent host cell specific usage of type IV secretion system components. <italic>VirB2 </italic>is the major protein that makes up the T4SS pilus, and has been shown to be necessary for full virulence in <italic>Brucella abortus </italic>[##REF##15322008##41##]. In <italic>Ap</italic>, seven of the eight <italic>virB2 </italic>paralogs are annotated as being <italic>TrbC/VirB2 </italic>(<italic>pfam04956</italic>) family members on the Entrez Protein entries for each individual protein. <italic>APH_1145</italic>, although not annotated as <italic>virB2</italic>, shares homology with and is located next to the other seven. Several other bacteria within the family <italic>Anaplasmataceae </italic>also possess multiple paralogs of <italic>virB2</italic>, which is unusual, as the majority of bacteria with type IV secretion systems have only one or two <italic>virB2 </italic>genes. A blast search done with <italic>APH_1133 </italic>shows, for example, that <italic>Anaplasma marginale</italic>, as well as <italic>Ehrlichia </italic>and <italic>Wolbachia </italic>species, also have multiple loci annotated as <italic>TrbC/VirB2 </italic>family members (see additional file ##SUPPL##7##8##: Examples of other <italic>Anaplasmataceae </italic>bacteria with multiple <italic>virB2 </italic>loci). These bacteria might also express specific <italic>virB2 </italic>paralogs in a host cell dependent manner.</p> <p>The absence of <italic>p44 </italic>transcription in ISE6 at the <italic>p44 </italic>expression locus and clear transcription in HL-60 and HMEC-1, is consistent with the observation that the tick cell samples produced little or no hybridization to <italic>p44 </italic>HVRs, while the human samples did, and indicates that in ISE6 little if any transcript was generated from any of the 22 full-length <italic>p44 </italic>genes. The lack of <italic>ApxR </italic>transcript in the tick cells is consistent with the findings of Wang et al., who performed quantitative reverse transcription PCR on <italic>Ap</italic>-infected ISE6 cells and tick salivary glands and found that <italic>ApxR </italic>is not transcribed [##REF##17905983##34##]. It was suggested that <italic>ApxR </italic>generally regulates transcription in mammalian host cells and specifically regulates <italic>p44 </italic>transcription by binding to the <italic>tr1 </italic>promoter. The strong transcription of <italic>tr1 </italic>in the tick cells in this study may be due to a lack of suppression by <italic>ApxR</italic>, which is not transcribed in the tick cells. The function of <italic>tr1</italic>, therefore, is unclear.</p> <p>The apparent over-representation of transcript from conserved <italic>p44 </italic>sequences, along with its reflecting behavior in the anti-sense strand, is unexpected. It may be the result of transcriptional "read-through" followed by the formation of stable double stranded, conserved sequence RNA. Bacteria are known to have poor control over transcription termination, and transcription of anti-sense sequence has been identified in <italic>Mycoplasma genitalium </italic>[##REF##17660438##42##]. Since <italic>p44 </italic>paralogs are scattered throughout the genome on both DNA strands, any adjacent gene transcription that continues into sense or anti-sense <italic>p44 </italic>sequences will create "false transcripts," the conserved sequences of which are complementary. Conserved anti-sense false transcript may anneal to conserved sense "true" and false transcript to form double stranded conserved sequence RNA, which is relatively stable compared to single stranded RNA and thus would accumulate in the bacteria (Figure ##FIG##8##9##: Diagram of possible mechanism to explain the over-representation of <italic>p44 </italic>conserved sequence transcripts and their anti-sense counterparts). Sense and anti-sense <italic>p44 </italic>false transcripts could come from many of the numerous <italic>p44 </italic>paralogs, but a possible source of anti-sense <italic>p44 </italic>transcript in the tick cells is via read-through from the <italic>msp4 </italic>gene (see additional file ##SUPPL##8##9##: <italic>msp4 </italic>transcription), which is opposite and just downstream of <italic>p44-15b </italic>and <italic>p44-13</italic>, strongly transcribed in the tick cells, not transcribed in the human cells, and has no obvious transcription terminator.</p> <p>It is possible that the anti-sense transcription noted in some genes, along with the prominent <italic>p44 </italic>transcription phenomena, function to regulate gene expression. In prokaryotes, <italic>cis</italic>- and <italic>trans</italic>-encoded anti-sense transcripts regulate coding sequence lying directly opposite or elsewhere in the genome, respectively [##REF##17387036##43##]. Although anti-sense mediated expression regulation mechanisms are poorly understood, some possible modes have been discussed and include: imprinting through DNA methylation, RNA processing interference, and ribosome interference [##REF##17139297##44##,##REF##17073772##45##]. In the case of <italic>p44</italic>, anti-sense transcripts may serve to silence leaky expression occurring from any of the 22 identified full-length <italic>p44 </italic>paralogs [##REF##16482227##1##], which are apparently capable of being expressed independently from the <italic>p44 </italic>expression locus [##REF##10364227##46##]. <italic>P44 </italic>silencing may be especially important in tick cells and account for the particular abundance of anomalous <italic>p44 </italic>conserved sequence transcripts in <italic>Ap</italic>-ISE6, which showed no <italic>p44 </italic>HVR transcription. Sense and anti-sense RNA homologous to the conserved ends of the <italic>p44 </italic>genes may even facilitate the process of non-reciprocal recombination by which <italic>p44 </italic>paralogs move into and out of the expression locus [##REF##16552034##47##]. If they are not purposeful, it is likely that these gene transcription phenomena are the result of poorly controlled transcription or are artifacts of the tiling arrays. The repeat sequence associated sense and anti-sense "transcripts" do appear to be an artifact of the tiling arrays, as they are consistently seen wherever repeat sequences occur, whether inside or outside of coding sequences. However, the transcription behavior of <italic>p44 </italic>is unique in the genome, and most genes do not display anti-sense transcription, therefore the over-representation of transcript to conserved <italic>p44 </italic>sequences and its reflection, and the anti-sense transcription of some genes, are intriguing and merit further investigation.</p>
<title>Conclusion</title> <p>Obligate intracellular pathogens like <italic>Ap </italic>control the cells they parasitize – to prevent immune attacks, divert cellular resources, and prevent host cells from apoptosing. Our understanding of tick genes is poor so it is not surprising that the up-regulated <italic>Ap </italic>genes in tick cells are nearly all "hypothetical." Matched with our limited understanding of <italic>Ap </italic>genes, the tick cell data are particularly difficult to interpret. Conversely, it makes sense that the most differentially active <italic>Ap </italic>genes in HL-60 cells are better characterized, since human cell lines have mainly been used to study the biology of <italic>Ap</italic>, and, perhaps, <italic>Ap </italic>genes evolved to interact in human cells would tend to be related to characterized effectors. It also makes sense that the differentially transcribed <italic>Ap </italic>genes in HL-60 and ISE6 are over-represented by membrane associated gene products, since survival in such disparate host cells would seem to require substantial specialization at the interface of the organism with its host cell: the bacterial membrane. The fact that a majority of <italic>Ap </italic>genes have no known function poses the greatest challenge to interpreting these data. However, some things are clear: 1. Genes differentially transcribed between the human and tick cells disproportionately represent surface proteins (~45% compared to ~25% of all proteins) (Table ##TAB##1##2##). 2. There are genes, paralogs, and operons exclusively transcribed in the tick and the human cells, some of which may encode excellent vaccine candidates. 3. The particular paralogs of the <italic>p44 </italic>family of membrane proteins (114) expressed in a population of <italic>Ap </italic>may be identified by the elevated signal produced within the HVR of each as compared to silent paralogs. 4. Whole RNA isolated from <italic>Ap </italic>infected host cells can be used to reveal details of bacterial gene transcription, including that from anti-sense sequences. 5. Global transcription profiles can likely be generated for <italic>Ap </italic>in any host cells, and for all aspects of the cell infection cycle – cell binding, entry, growth, and escape – although some enrichment for bacteria or bacterial mRNA may be necessary. Coupling <italic>Ap </italic>transcription data with that of infected host cells will facilitate the discovery of <italic>Ap </italic>and host cell gene functions.</p> <p>Having transcription data for all of an organism's DNA sequence allows a line graph display for both DNA strands parallel to an annotated map of the genome. This way one can readily see transcriptional behavior that may be less accessible through other analysis tools. For example, anti-sense transcription, and the variation in transcription profiles of genes – sloped, flat, horned, and reflected – may lead to important insights into <italic>Ap </italic>gene regulation, as well as for other intracellular organisms that subvert host cell processes for their own benefit.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p><italic>Anaplasma phagocytophilum </italic>(<italic>Ap</italic>) is an obligate intracellular bacterium and the agent of human granulocytic anaplasmosis, an emerging tick-borne disease. <italic>Ap </italic>alternately infects ticks and mammals and a variety of cell types within each. Understanding the biology behind such versatile cellular parasitism may be derived through the use of tiling microarrays to establish high resolution, genome-wide transcription profiles of the organism as it infects cell lines representative of its life cycle (tick; ISE6) and pathogenesis (human; HL-60 and HMEC-1).</p> <title>Results</title> <p>Detailed, host cell specific transcriptional behavior was revealed. There was extensive differential <italic>Ap </italic>gene transcription between the tick (ISE6) and the human (HL-60 and HMEC-1) cell lines, with far fewer differentially transcribed genes between the human cell lines, and all disproportionately represented by membrane or surface proteins. There were <italic>Ap </italic>genes exclusively transcribed in each cell line, apparent human- and tick-specific operons and paralogs, and anti-sense transcripts that suggest novel expression regulation processes. Seven <italic>virB2 </italic>paralogs (of the bacterial type IV secretion system) showed human or tick cell dependent transcription. Previously unrecognized genes and coding sequences were identified, as were the expressed <italic>p44/msp2 </italic>(major surface proteins) paralogs (of 114 total), through elevated signal produced to the unique hypervariable region of each – 2/114 in HL-60, 3/114 in HMEC-1, and none in ISE6.</p> <title>Conclusion</title> <p>Using these methods, whole genome transcription profiles can likely be generated for <italic>Ap</italic>, as well as other obligate intracellular organisms, in any host cells and for all stages of the cell infection process. Visual representation of comprehensive transcription data alongside an annotated map of the genome renders complex transcription into discernable patterns.</p>
<title>Authors' contributions</title> <p>CN conceived of and carried out the tiling array study, led the array design effort, participated in data analysis, and drafted the manuscript. MH participated in the design of the study and the array, transformed the raw array data into interpretable formats, including adapting it to the Artemis annotation software, analyzed the data, and assisted in drafting the manuscript. RF participated in designing the study and the array, and assisted in the analysis and interpretation of the data, and in the preparation of the manuscript. BS performed the quantitative RT-PCR experiments. SG assisted with data analysis. AOC performed the <italic>virB2 </italic>sequence comparisons. TK and UM assisted in designing the study, the array, and in the preparation of the manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank: Wayne Xu (University of Minnesota, Super Computing Institute, University budget funded), for helping with bioinformatics issues and data interpretation; Arkady Khodursky (Biochemistry, Molecular Biology, and Biophysics, University of Minnesota), for his advice on data normalization and interpretation; and Gerald Baldridge (NIH grant Nr. AIO42792) for his help with editing the manuscript. The work presented herein was funded by a grant from the National Research Fund for Tick-Borne Diseases, Inc. to TJK, and a grant from NIH, Nr. AIO42792, to UGM.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Microscopic images of Giemsa stained cells infected with <italic>Ap</italic>. (A) <italic>Ap</italic>-infected HMEC-1 (B) <italic>Ap</italic>-infected ISE6 (C) <italic>Ap-</italic>infected HL-60. Cell nuclei are labeled "N" and arrows point to <italic>Ap </italic>morulae. Scale bar = 10 μm</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Artemis profiles depicting the relative transcription levels of five <italic>Ap </italic>genes during infection of HL-60 (red), HMEC-1 (green), and ISE6 (blue) cells. Plots were "smoothed" by setting the sliding window average to 5. (A) Major outer membrane protein gene (<italic>omp-1A</italic>; <italic>APH_1359</italic>) transcription greater in the human cell lines compared to the tick cell line. (B) Outer membrane efflux protein (<italic>APH_1110</italic>) greater in the tick cell line compared to the human cell lines. (C) Transcription of the major surface protein 4 gene (<italic>msp4</italic>; <italic>APH_1240</italic>) only in the tick cell line. (D) Transcription of the 60 kDa chaperonin gene (<italic>groL</italic>; <italic>APH_0240</italic>) was greatest in HL-60, significantly lower in HMEC-1, and least in ISE6. (E) Equal transcription of the succinyl-CoA synthetase beta subunit gene (<italic>sucC</italic>; <italic>APH_1052</italic>) in all three cell lines.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Artemis transcription plots showing examples of flat and sloped gene transcription profiles (Red: <italic>Ap</italic>-HL-60, Green: <italic>Ap</italic>-HMEC-1, Blue: <italic>Ap</italic>-ISE6; plots were "smoothed" by setting the sliding window average to 5). (A) Polynucleotide phosphorylase gene (<italic>pnp</italic>) with an over-all flat transcription profile in all three cell lines. (B) Two examples of genes – <italic>APH_0756 </italic>(hypothetical) and <italic>rpoH </italic>(heat shock sigma factor sigma 32) – with transcription profiles that slope downward from 5' to 3'.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p>Artemis transcription plots of characteristic <italic>p44 </italic>transcription profiles (Red: <italic>Ap</italic>-HL-60, Green: <italic>Ap</italic>-HMEC-1, Blue: <italic>Ap</italic>-ISE6; plots were "smoothed" by setting the sliding window average to 5). Arrows in panel B indicate <italic>p44 </italic>conserved sequence "horns" on the coding (minus) strand, and "reflected" horns (panel A) in the anti-sense (plus) strand. A strong signal (green) associated with the HVR in <italic>APH_1342 </italic>(*), likely indicates expression of the corresponding <italic>p44 </italic>paralog (<italic>p44-31</italic>) in HMEC-1. The lack of HVR associated signal in <italic>APH_1343</italic>, but strong conserved sequence associated signals (horns), is typical of most <italic>p44 </italic>paralogs. An unannotated segment of <italic>p44 </italic>conserved sequence lies between <italic>APH_1343 </italic>and <italic>APH_1344 </italic>(yellow) on the minus strand. It also showed strong sense (B) and anti-sense (A) signals. <italic>APH_1344 </italic>and <italic>APH_1345 </italic>show typical transcription profiles: signal on the sense strand (B) but not on the anti-sense strand (A).</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p>Artemis transcription plots of two genes showing "reflecting" transcription patterns on the anti-sense strands (Red: <italic>Ap</italic>-HL-60, Green: <italic>Ap</italic>-HMEC-1, Blue: <italic>Ap</italic>-ISE6; plots were "smoothed" by setting the sliding window average to 5). (A) HGE2 protein <italic>APH_1387</italic>. (B) Hypothetical protein APH 0546. Note that <italic>Ap </italic>in all three cell lines produced sense and anti-sense transcript for <italic>APH_1387 </italic>(panel A), while in the case of <italic>APH_0546 </italic>(panel B) only <italic>Ap</italic>-ISE6 produced sense and anti-sense transcript.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p>Artemis plots illustrating transcription activity at the <italic>p44 </italic>expression site, and at <italic>ApxR</italic>, a putative <italic>p44 </italic>transcription regulator. (Red: <italic>Ap</italic>-HL-60, Green: <italic>Ap</italic>-HMEC-1, Blue: <italic>Ap</italic>-ISE6; plots were "smoothed" by setting the sliding window average to 5). (A) In the human cell lines, <italic>Ap </italic>shows transcription beginning upstream of <italic>omp-1N </italic>(and <italic>p44-18ES</italic>, the <italic>p44 </italic>expression locus) near coordinate 1289280, but there is no specific transcription in the tick cell line. Transcription regulator <italic>tr1 </italic>(<italic>APH_1218</italic>) is not transcribed in the human cell lines but is in the tick cell line. (B) <italic>ApxR </italic>(<italic>APH_0515</italic>), a putative regulator of <italic>p44 </italic>transcription – through binding to and inhibiting the <italic>tr1 </italic>promoter – shows low-level transcription in the human cell lines but none in the tick cell line.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p>Phylogenetic tree showing the relationship, based on amino acid sequence, of eight <italic>virB2 </italic>paralogs in the <italic>Ap </italic>genome. Four were transcribed only in ISE6 (<italic>APH_1131 – APH_1134</italic>), and two only in HL-60 and HMEC-1 (<italic>APH_1144 </italic>and <italic>APH_1145</italic>). No transcript from <italic>APH_1130 </italic>or <italic>APH_1136 </italic>was measured. The tick cell line associated paralogs are closely related to each other, while those transcribed in the human cell lines form a separate group and are less related to each other. The tree was constructed with PAUP 4.0 using neighbor-joining: absolute variation. Values shown in branches correspond to 2000 bootstraps analysis.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p>Tiling array (area under gene curve) vs. qRT-PCR (40 minus threshold cycle) measurements of transcript levels of five <italic>Ap </italic>genes (key to bars indicated) during growth in HL-60, HMEC-1, and ISE6 cells. Relative transcript levels for the five selected <italic>Ap </italic>genes, within and between cell lines, confirm those indicated by the array data. qRT-PCR data was converted by subtracting the Ct (threshold cycle) from forty (total PCR cycles), since lower threshold cycles correspond to higher transcript levels.</p></caption></fig> <fig position="float" id="F9"><label>Figure 9</label><caption><p>A proposed model for generation of the observed anomalous <italic>p44 </italic>conserved sequence transcripts (sense and anti-sense). "Read through" transcription of genes lying just upstream of anti-sense <italic>p44 </italic>sequence (e.g. "<italic>Gene 1</italic>" and <italic>msp4</italic>) may produce anti-sense <italic>p44 </italic>transcript, which, together with <italic>p44 </italic>sense transcript, forms double stranded RNA (dsRNA). Because the HVR sequences are not complementary they do not form dsRNA and are therefore degraded. However, the conserved, complementary sequences do form dsRNA so are stabilized, accumulate in the bacteria, and are measured as over-abundant by the arrays.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Summary of differential <italic>Ap </italic>gene transcription between HL-60, HMEC-1, and ISE6</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="right"># ORFs<break/>Differentially <break/>transcribed<break/>(p ≤ 0.05)</td><td align="right">% of total <break/>ORFs<break/><break/></td><td align="right"># ORFs ≥ 2-fold   <break/>(p ≤ 0.05)<break/>differentially<break/>transcribed</td><td align="right">% of total <break/>ORFs<break/><break/></td></tr></thead><tbody><tr><td align="left"><italic>Ap</italic>-HL-60 vs. <italic>Ap</italic>-HMEC-1</td><td align="right"><bold>71</bold></td><td align="right">5.0</td><td align="right"><bold>1</bold></td><td align="right">0.07</td></tr><tr><td align="left"><italic>Ap</italic>-HL-60 vs. <italic>Ap</italic>-ISE6</td><td align="right"><bold>585</bold></td><td align="right">41.5</td><td align="right"><bold>117</bold></td><td align="right">8.5</td></tr><tr><td align="left"><italic>Ap</italic>-HMEC-1 vs. <italic>Ap</italic>-ISE6</td><td align="right"><bold>304</bold></td><td align="right">21.5</td><td align="right"><bold>61</bold></td><td align="right">4.3</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Summary of <italic>Ap</italic>-HL-60 vs. <italic>Ap</italic>-ISE6 differential gene transcription</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left"><bold>Gene Product</bold></td><td align="center"><bold>Locus</bold></td><td align="left"><bold>Predicted Cellular Location</bold></td><td align="center"><bold>Fold Change Ap-HL-60/Ap-ISE6</bold></td></tr></thead><tbody><tr><td align="right">1</td><td align="left">DNA-binding protein</td><td align="center">APH_1100</td><td align="left">Cytoplasmic</td><td align="center">4.9</td></tr><tr><td align="right">2</td><td align="left">HGE-14 protein</td><td align="center">APH_0387</td><td align="left"><bold>Extracellular</bold></td><td align="center">4.8</td></tr><tr><td align="right">3</td><td align="left">hypothetical protein</td><td align="center">APH_1412</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">4.2</td></tr><tr><td align="right">4</td><td align="left">hypothetical protein</td><td align="center">APH_0915</td><td align="left"><bold>Outer Membrane, Extracellular</bold></td><td align="center">4.1</td></tr><tr><td align="right">5</td><td align="left">hypothetical protein</td><td align="center">APH_0906</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">4.1</td></tr><tr><td align="right">6</td><td align="left">major outer membrane protein OMP-1A</td><td align="center">APH_1359</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">4.0</td></tr><tr><td align="right">7</td><td align="left">hypothetical protein</td><td align="center">APH_1378</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">3.8</td></tr><tr><td align="right">8</td><td align="left">hypothetical protein</td><td align="center">APH_0842</td><td align="left">Cytoplasmic</td><td align="center">3.7</td></tr><tr><td align="right">9</td><td align="left">hypothetical protein</td><td align="center">APH_0838</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">3.6</td></tr><tr><td align="right">10</td><td align="left">hypothetical protein</td><td align="center">APH_0388</td><td align="left">Cytoplasmic</td><td align="center">3.6</td></tr><tr><td align="right">11</td><td align="left">hypothetical protein</td><td align="center">APH_1145</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">3.5</td></tr><tr><td align="right">12</td><td align="left">OmpA family protein</td><td align="center">APH_0338</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">3.4</td></tr><tr><td align="right">13</td><td align="left">DNA-binding response regulator</td><td align="center">APH_1099</td><td align="left">Cytoplasmic</td><td align="center">3.4</td></tr><tr><td align="right">14</td><td align="left">hypothetical protein</td><td align="center">APH_1144</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">3.4</td></tr><tr><td align="right">15</td><td align="left">hypothetical protein</td><td align="center">APH_0837</td><td align="left">Cytoplasmic</td><td align="center">3.4</td></tr><tr><td align="right">16</td><td align="left">HGE-14 protein</td><td align="center">APH_0382</td><td align="left"><bold>Extracellular</bold></td><td align="center">3.3</td></tr><tr><td align="right">17</td><td align="left">hypothetical protein</td><td align="center">APH_0005</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">3.3</td></tr><tr><td align="right">18</td><td align="left">hypothetical protein</td><td align="center">APH_0756</td><td align="left"><bold>Inner Membrane, Cytoplasmic</bold></td><td align="center">3.2</td></tr><tr><td align="right">19</td><td align="left">10 kDa chaperonin</td><td align="center">APH_0241</td><td align="left"><bold>Periplasmic</bold></td><td align="center">3.2</td></tr><tr><td align="right">20</td><td align="left">hypothetical protein</td><td align="center">APH_0032</td><td align="left"><bold>Outer Membrane, Extracellular</bold></td><td align="center">3.1</td></tr><tr><td align="right">21</td><td align="left">hypothetical protein</td><td align="center">APH_0874</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">3.1</td></tr><tr><td align="right">22</td><td align="left">hypothetical protein</td><td align="center">APH_0233</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">3.1</td></tr><tr><td align="right">23</td><td align="left">HGE-14 protein</td><td align="center">APH_0385</td><td align="left">Cytoplasmic</td><td align="center">3.1</td></tr><tr><td align="right">24</td><td align="left">signal peptidase II</td><td align="center">APH_1160</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">3.0</td></tr><tr><td align="right">25</td><td align="left">hypothetical protein</td><td align="center">APH_1156</td><td align="left">Cytoplasmic</td><td align="center">3.0</td></tr><tr><td align="right">26</td><td align="left">hypothetical protein</td><td align="center">APH_0793</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">2.9</td></tr><tr><td align="right">27</td><td align="left">HGE-14 protein</td><td align="center">APH_0455</td><td align="left"><bold>Extracellular</bold></td><td align="center">2.9</td></tr><tr><td align="right">28</td><td align="left">Omp-1N</td><td align="center">APH_1220</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">2.8</td></tr><tr><td align="right">29</td><td align="left">hypothetical protein</td><td align="center">APH_0949</td><td align="left"><bold>Inner Membrane, Cytoplasmic</bold></td><td align="center">2.7</td></tr><tr><td align="right">30</td><td align="left">hypothetical protein</td><td align="center">APH_0033</td><td align="left">Cytoplasmic</td><td align="center">2.7</td></tr><tr><td align="right">31</td><td align="left">hypothetical protein</td><td align="center">APH_1307</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">2.7</td></tr><tr><td align="right">32</td><td align="left">hypothetical protein</td><td align="center">APH_1157</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">2.7</td></tr><tr><td align="right">33</td><td align="left">hypothetical protein</td><td align="center">APH_1151</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">2.6</td></tr><tr><td align="right">34</td><td align="left">antioxidant AhpC/Tsa family</td><td align="center">APH_0795</td><td align="left">Cytoplasmic</td><td align="center">2.6</td></tr><tr><td align="right">35</td><td align="left">RNA polymerase sigma-32 factor</td><td align="center">APH_0759</td><td align="left">Cytoplasmic</td><td align="center">2.6</td></tr><tr><td align="right">36</td><td align="left">60 kDa chaperonin</td><td align="center">APH_0240</td><td align="left">Cytoplasmic</td><td align="center">2.6</td></tr><tr><td align="right">37</td><td align="left">hypothetical protein</td><td align="center">APH_1235</td><td align="left">Cytoplasmic</td><td align="center">2.5</td></tr><tr><td align="right">38</td><td align="left">hypothetical protein</td><td align="center">APH_0922</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">2.5</td></tr><tr><td align="right">39</td><td align="left">hypothetical protein</td><td align="center">APH_1262</td><td align="left">Cytoplasmic</td><td align="center">2.5</td></tr><tr><td align="right">40</td><td align="left">hypothetical protein</td><td align="center">APH_0757</td><td align="left">Cytoplasmic</td><td align="center">2.5</td></tr><tr><td align="right">41</td><td align="left">chaperone protein DnaK</td><td align="center">APH_0346</td><td align="left">Cytoplasmic</td><td align="center">2.4</td></tr><tr><td align="right">42</td><td align="left">hypothetical protein</td><td align="center">APH_1236</td><td align="left">Cytoplasmic</td><td align="center">2.4</td></tr><tr><td align="right">43</td><td align="left">hypothetical protein</td><td align="center">APH_0363</td><td align="left">Cytoplasmic</td><td align="center">2.4</td></tr><tr><td align="right">44</td><td align="left">translation initiation factor IF-3</td><td align="center">APH_1263</td><td align="left">Cytoplasmic</td><td align="center">2.4</td></tr><tr><td align="right">45</td><td align="left">glyceraldehyde-3-phosphate dehydrogenase type I</td><td align="center">APH_1349</td><td align="left">Cytoplasmic</td><td align="center">2.4</td></tr><tr><td align="right">46</td><td align="left">hypothetical protein</td><td align="center">APH_0873</td><td align="left">Cytoplasmic</td><td align="center">2.4</td></tr><tr><td align="right">47</td><td align="left">hypothetical protein</td><td align="center">APH_0919</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">2.4</td></tr><tr><td align="right">48</td><td align="left">hypothetical protein</td><td align="center">APH_1072</td><td align="left">Cytoplasmic</td><td align="center">2.4</td></tr><tr><td align="right">49</td><td align="left">hypothetical protein</td><td align="center">APH_1320</td><td align="left">Cytoplasmic</td><td align="center">2.3</td></tr><tr><td align="right">50</td><td align="left">HGE-14 protein</td><td align="center">APH_0453</td><td align="left">Cytoplasmic</td><td align="center">2.3</td></tr><tr><td align="right">51</td><td align="left">outer membrane protein MSP2 family</td><td align="center">APH_1325</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">2.2</td></tr><tr><td align="right">52</td><td align="left">hypothetical protein</td><td align="center">APH_0643</td><td align="left">Cytoplasmic</td><td align="center">2.2</td></tr><tr><td align="right">53</td><td align="left">hypothetical protein</td><td align="center">APH_0839</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">2.2</td></tr><tr><td align="right">54</td><td align="left">putative acyl carrier protein</td><td align="center">APH_0929</td><td align="left">Cytoplasmic</td><td align="center">2.2</td></tr><tr><td align="right">55</td><td align="left">Es1 family protein</td><td align="center">APH_0006</td><td align="left">Cytoplasmic</td><td align="center">2.2</td></tr><tr><td align="right">56</td><td align="left">hypothetical protein</td><td align="center">APH_0179</td><td align="left">Cytoplasmic</td><td align="center">2.2</td></tr><tr><td align="right">57</td><td align="left">iron-sulfur cluster assembly accessory protein</td><td align="center">APH_0676</td><td align="left">Cytoplasmic</td><td align="center">2.2</td></tr><tr><td align="right">58</td><td align="left">putative ATP synthase F0 B' subunit</td><td align="center">APH_1190</td><td align="left">Cytoplasmic</td><td align="center">2.1</td></tr><tr><td align="right">59</td><td align="left">hypothetical protein</td><td align="center">APH_0719</td><td align="left">Cytoplasmic</td><td align="center">2.1</td></tr><tr><td align="right">60</td><td align="left">hypothetical protein</td><td align="center">APH_0991</td><td align="left">Cytoplasmic</td><td align="center">2.1</td></tr><tr><td align="right">61</td><td align="left">succinate dehydrogenase cytochrome b556 subunit</td><td align="center">APH_0999</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">2.1</td></tr><tr><td align="right">62</td><td align="left">pyruvate phosphate dikinase</td><td align="center">APH_0185</td><td align="left">Cytoplasmic</td><td align="center">2.1</td></tr><tr><td align="right">63</td><td align="left">iron-binding protein</td><td align="center">APH_0051</td><td align="left">Cytoplasmic</td><td align="center">2.1</td></tr><tr><td align="right">64</td><td align="left">nucleoside diphosphate kinase</td><td align="center">APH_1217</td><td align="left">Cytoplasmic</td><td align="center">2.1</td></tr><tr><td align="right">65</td><td align="left">malonyl CoA-acyl carrier protein transacylase</td><td align="center">APH_0092</td><td align="left">Cytoplasmic</td><td align="center">2.0</td></tr><tr><td align="right">66</td><td align="left">hypothetical protein</td><td align="center">APH_0786</td><td align="left">Cytoplasmic</td><td align="center">2.0</td></tr><tr><td align="right">67</td><td align="left">co-chaperone GrpE</td><td align="center">APH_0036</td><td align="left">Cytoplasmic</td><td align="center">2.0</td></tr><tr><td align="right">68</td><td align="left">hypothetical protein</td><td align="center">APH_0771</td><td align="left">Cytoplasmic</td><td align="center">2.0</td></tr><tr><td align="right">69</td><td align="left">hypothetical protein</td><td align="center">APH_0585</td><td align="left">Cytoplasmic</td><td align="center">2.0</td></tr><tr><td align="right">70</td><td align="left">hypothetical protein</td><td align="center">APH_0655</td><td align="left">Cytoplasmic</td><td align="center">2.0</td></tr><tr><td align="right">71</td><td align="left">ribonucleoside-diphosphate reductase alpha subunit</td><td align="center">APH_0331</td><td align="left">Cytoplasmic</td><td align="center">2.0</td></tr><tr><td align="right">72</td><td align="left">P44-45 outer membrane protein</td><td align="center">APH_0171</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">2.0</td></tr><tr><td align="right">73</td><td align="left">adenylosuccinate lyase</td><td align="center">APH_0867</td><td align="left">Cytoplasmic</td><td align="center">2.0</td></tr><tr><td align="right">74</td><td align="left">P44-36 outer membrane protein</td><td align="center">APH_1168</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">2.0</td></tr><tr><td align="right">75</td><td align="left">aspartate aminotransferase</td><td align="center">APH_0660</td><td align="left">Cytoplasmic</td><td align="center">2.0</td></tr><tr><td align="right">76</td><td align="left">cytochrome C membrane-bound</td><td align="center">APH_0180</td><td align="left"><bold>Periplasmic</bold></td><td align="center">2.0</td></tr><tr><td/><td align="center" colspan="2"><bold>35/76 named genes ≥ 2-fold up in <italic>Ap</italic>-HL-60 = 54% hypothetical</bold></td><td align="center" colspan="2"><bold>33/76 genes membrane associated = 43%</bold></td></tr><tr><td align="right">77</td><td align="left">hypothetical protein</td><td align="center">APH_0197</td><td align="left"><bold>Periplasmic</bold></td><td align="center">0.5</td></tr><tr><td align="right">78</td><td align="left">hypothetical protein</td><td align="center">APH_0369</td><td align="left">Cytoplasmic</td><td align="center">0.5</td></tr><tr><td align="right">79</td><td align="left">hypothetical protein</td><td align="center">APH_0497</td><td align="left">Cytoplasmic</td><td align="center">0.5</td></tr><tr><td align="right">80</td><td align="left">hypothetical protein</td><td align="center">APH_0425</td><td align="left">Cytoplasmic</td><td align="center">0.5</td></tr><tr><td align="right">81</td><td align="left">hypothetical protein</td><td align="center">APH_0587</td><td align="left">Cytoplasmic</td><td align="center">0.5</td></tr><tr><td align="right">82</td><td align="left">hypothetical protein</td><td align="center">APH_0963</td><td align="left">Cytoplasmic</td><td align="center">0.5</td></tr><tr><td align="right">83</td><td align="left">hypothetical protein</td><td align="center">APH_1130</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">0.5</td></tr><tr><td align="right">84</td><td align="left">hypothetical protein</td><td align="center">APH_0467</td><td align="left">Cytoplasmic</td><td align="center">0.5</td></tr><tr><td align="right">85</td><td align="left">thiamine biosynthesis protein ThiC truncation</td><td align="center">APH_0586</td><td align="left">Cytoplasmic</td><td align="center">0.5</td></tr><tr><td align="right">86</td><td align="left">hypothetical protein</td><td align="center">APH_0806</td><td align="left"><bold>Periplasmic</bold></td><td align="center">0.5</td></tr><tr><td align="right">87</td><td align="left">hypothetical protein</td><td align="center">APH_0599</td><td align="left">Cytoplasmic</td><td align="center">0.5</td></tr><tr><td align="right">88</td><td align="left">hypothetical protein</td><td align="center">APH_0827</td><td align="left">Cytoplasmic</td><td align="center">0.4</td></tr><tr><td align="right">89</td><td align="left">outer membrane efflux protein</td><td align="center">APH_1110</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">0.4</td></tr><tr><td align="right">90</td><td align="left">hypothetical protein</td><td align="center">APH_1131</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">0.4</td></tr><tr><td align="right">91</td><td align="left">hypothetical protein</td><td align="center">APH_0829</td><td align="left">Cytoplasmic</td><td align="center">0.4</td></tr><tr><td align="right">92</td><td align="left">hypothetical protein</td><td align="center">APH_0818</td><td align="left">Cytoplasmic</td><td align="center">0.4</td></tr><tr><td align="right">93</td><td align="left">hypothetical protein</td><td align="center">APH_0841</td><td align="left">Cytoplasmic</td><td align="center">0.4</td></tr><tr><td align="right">94</td><td align="left">hypothetical protein</td><td align="center">APH_1382</td><td align="left">Cytoplasmic</td><td align="center">0.4</td></tr><tr><td align="right">95</td><td align="left">hypothetical protein</td><td align="center">APH_0550</td><td align="left">Cytoplasmic</td><td align="center">0.4</td></tr><tr><td align="right">96</td><td align="left">hypothetical protein</td><td align="center">APH_0485</td><td align="left">Cytoplasmic</td><td align="center">0.4</td></tr><tr><td align="right">97</td><td align="left">hypothetical protein</td><td align="center">APH_0355</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">0.4</td></tr><tr><td align="right">98</td><td align="left">hypothetical protein</td><td align="center">APH_1132</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">0.3</td></tr><tr><td align="right">99</td><td align="left">hypothetical protein</td><td align="center">APH_0720</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">0.3</td></tr><tr><td align="right">100</td><td align="left">hypothetical protein</td><td align="center">APH_1384</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">0.3</td></tr><tr><td align="right">101</td><td align="left">hypothetical protein</td><td align="center">APH_1380</td><td align="left">Cytoplasmic</td><td align="center">0.3</td></tr><tr><td align="right">102</td><td align="left">hypothetical protein</td><td align="center">APH_1370</td><td align="left">Cytoplasmic</td><td align="center">0.3</td></tr><tr><td align="right">103</td><td align="left">hypothetical protein</td><td align="center">APH_0320</td><td align="left">Cytoplasmic</td><td align="center">0.3</td></tr><tr><td align="right">104</td><td align="left">hypothetical protein</td><td align="center">APH_0726</td><td align="left"><bold>Membrane</bold></td><td align="center">0.3</td></tr><tr><td align="right">105</td><td align="left">hypothetical protein</td><td align="center">APH_1369</td><td align="left">Cytoplasmic</td><td align="center">0.3</td></tr><tr><td align="right">106</td><td align="left">hypothetical protein</td><td align="center">APH_1368</td><td align="left">Cytoplasmic</td><td align="center">0.3</td></tr><tr><td align="right">107</td><td align="left">hypothetical protein</td><td align="center">APH_1385</td><td align="left">Cytoplasmic</td><td align="center">0.2</td></tr><tr><td align="right">108</td><td align="left">hypothetical protein</td><td align="center">APH_0724</td><td align="left"><bold>Membrane</bold></td><td align="center">0.2</td></tr><tr><td align="right">109</td><td align="left">hypothetical protein</td><td align="center">APH_0805</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">0.2</td></tr><tr><td align="right">110</td><td align="left">hypothetical protein</td><td align="center">APH_0723</td><td align="left"><bold>Membrane</bold></td><td align="center">0.2</td></tr><tr><td align="right">111</td><td align="left">hypothetical protein</td><td align="center">APH_0487</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">0.2</td></tr><tr><td align="right">112</td><td align="left">hypothetical protein</td><td align="center">APH_1386</td><td align="left">Cytoplasmic</td><td align="center">0.2</td></tr><tr><td align="right">113</td><td align="left">hypothetical protein</td><td align="center">APH_0177</td><td align="left"><bold>Extracellular</bold></td><td align="center">0.1</td></tr><tr><td align="right">114</td><td align="left">hypothetical protein</td><td align="center">APH_0546</td><td align="left"><bold>Extracellular</bold></td><td align="center">0.1</td></tr><tr><td align="right">115</td><td align="left">major surface protein 4</td><td align="center">APH_1240</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">0.1</td></tr><tr><td align="right">116</td><td align="left">hypothetical protein</td><td align="center">APH_0916</td><td align="left"><bold>Inner Membrane</bold></td><td align="center">0.1</td></tr><tr><td align="right">117</td><td align="left">hypothetical protein</td><td align="center">APH_0406</td><td align="left"><bold>Outer Membrane</bold></td><td align="center">0.1</td></tr><tr><td/><td align="center" colspan="2"><bold>3/41 named genes ≥ 2-fold up in <italic>Ap</italic>-ISE6 = 93% hypothetical</bold></td><td align="center" colspan="2"><bold>19/41 genes membrane associated = 46%</bold></td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Materials and instructions for graphing transcript level data in Artemis</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>qRT-PCR primers</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p><italic>Ap-</italic>HL-60 vs. <italic>Ap</italic>-HMEC-1 differential transcription. 71 ORFs differentially transcribed (p ≤ 0.05) between <italic>Ap</italic>-HL-60 and <italic>Ap</italic>-HMEC-1</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Genes and operons with no detected transcripts. Artemis transcription plots showing examples of genes and an operon with no detectable transcript signal in any of the cell lines (Red: <italic>Ap</italic>-HL-60, Green: <italic>Ap</italic>-HMEC-1, Blue: <italic>Ap</italic>-ISE6; "smoothed" using a sliding window average of 5). (A) Genes <italic>radC </italic>(DNA repair) and <italic>lipB </italic>(lipoyl (octanoyl)-acyl carrier protein B). (B) An operon including loci <italic>APH_0778 – APH_0783</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>Repeat-sequence-based sense and anti sense signal. Artemis plots showing sense and anti-sense transcript signal (shaded) within repeat sequences in two ORFs (Red: <italic>Ap</italic>-HL-60, Green: <italic>Ap</italic>-HMEC-1, Blue: <italic>Ap</italic>-ISE6; plots were "smoothed" by setting the sliding window average to 5). (A) <italic>APH_0377 </italic>(hypothetical). (B) <italic>APH_0455 </italic>(hypothetical).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p>Amino acid sequence alignments for the eight <italic>virB2 </italic>paralogs of <italic>Ap</italic>. (A) Pair-wise alignment of the eight <italic>virB2 </italic>amino acid sequences showing the percentages of identity and similarity between each. Identities rank from a high of 93% between tick cell expressed paralogs <italic>APH_1133 </italic>and <italic>APH_1134</italic>, and a low of 22% between non-expressed <italic>APH_1136 </italic>and human cell expressed <italic>APH_1145</italic>. (B) Multiple alignment showing higher identity and similarity between the C termini of paralogs, which contain the functional portion of the proteins. * Identical amino acids. · Conservative substitution.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p>Tick- and human-specific <italic>Ap </italic>operons. Artemis transcription plots showing <italic>Ap </italic>operons specific to the tick cell line (A, B) and the human cell lines (C). (Red: <italic>Ap</italic>-HL-60, Green: <italic>Ap</italic>-HMEC-1, Blue: <italic>Ap</italic>-ISE6.) (A) A tick operon (<italic>APH_1380 </italic>to <italic>APH_1386</italic>) that appears to include an unannotated ORF between <italic>APH_1380 </italic>and <italic>APH_1383 </italic>(in the same reading frame as <italic>APH_1380</italic>). Annotated loci <italic>APH_1381 </italic>and <italic>APH_1382 </italic>appear to be false ORFs. (B) Another tick-specific operon that includes loci <italic>APH_0720 </italic>to <italic>APH_0726</italic>. Loci <italic>APH_0721 </italic>and <italic>APH_0722 </italic>showed no transcription. (C) <italic>Ap </italic>genes transcribed only in the human cells: <italic>APH_0837</italic>, <italic>APH_0838</italic>, <italic>APH_0839</italic>, and <italic>APH_0842</italic>. Plots on panel A were smoothed using a sliding average of 5, and on panels B and C using a sliding average of 10.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional file 8</title><p>Examples of other <italic>Anaplasmataceae </italic>bacteria with multiple <italic>virB2 </italic>loci. Examples of <italic>Anaplasmataceae </italic>bacteria with multiple <italic>virB2 </italic>loci.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional file 9</title><p><italic>msp4 </italic>transcription. Artemis transcription plots showing the position of the <italic>msp4 </italic>gene just downstream of anti-sense <italic>p44 </italic>sequences and its strong transcription by <italic>Ap</italic>-ISE6 (Red: <italic>Ap</italic>-HL-60, Green: <italic>Ap</italic>-HMEC-1, Blue: <italic>Ap</italic>-ISE6; plots were "smoothed" by setting the sliding window average to 5). "Read through" transcription of <italic>msp4 </italic>may be a source of anti-sense <italic>p44 </italic>transcripts in <italic>Ap</italic>-ISE6.</p></caption></supplementary-material>
<table-wrap-foot><p>Summary of differential <italic>Ap </italic>gene transcription between HL-60, HMEC-1, and ISE6.</p></table-wrap-foot> <table-wrap-foot><p>Genes differentially transcribed (p ≤ 0.05, ≥ two-fold difference) between human (HL-60) and tick (ISE6) cells. Gene products (117) are listed in descending order of their transcript abundance in HL-60, with their fold change indicated in the right hand column. Gene products 1–76 were those more highly transcribed in HL-60, and 77–117 were those more highly transcribed in ISE6. Gene products predicted to be membrane associated (in bold) tended to be those most differentially transcribed – at the top (most abundant in <italic>Ap</italic>-HL-60) and bottom (most abundant in <italic>Ap</italic>-ISE6) of the list. The percentages of hypothetical genes, and the percentages of gene products that are membrane associated are also indicated. (~25% of all <italic>Ap </italic>genes code for proteins that are membrane associated.)</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-364-1"/> <graphic xlink:href="1471-2164-9-364-2"/> <graphic xlink:href="1471-2164-9-364-3"/> <graphic xlink:href="1471-2164-9-364-4"/> <graphic xlink:href="1471-2164-9-364-5"/> <graphic xlink:href="1471-2164-9-364-6"/> <graphic xlink:href="1471-2164-9-364-7"/> <graphic xlink:href="1471-2164-9-364-8"/> <graphic xlink:href="1471-2164-9-364-9"/>
<media xlink:href="1471-2164-9-364-S1.zip" mimetype="application" mime-subtype="x-zip-compressed"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-364-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-364-S3.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-364-S4.png" mimetype="image" mime-subtype="png"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-364-S5.png" mimetype="image" mime-subtype="png"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-364-S6.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-364-S7.png" mimetype="image" mime-subtype="png"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-364-S8.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-364-S9.png" mimetype="image" mime-subtype="png"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Cerrina", "Blattnerb", "Huanga", "Huea", "Greenc", "Singh-Gassonb", "Sussmanb"], "given-names": ["Fa", "F", "W", "Y", "R", "S", "M"], "suffix": ["1"], "article-title": ["Biological lithography: development of a maskless microarray synthesizer for DNA chips"], "source": ["Microelectronic Engineering"], "year": ["2002"], "volume": ["61\u201362"], "fpage": ["33"], "lpage": ["40"], "pub-id": ["10.1016/S0167-9317(02)00566-X"]}]
BMC Genomics. 2008 Jul 31; 9:364