<title>Introduction</title> <p>
<italic>Rickettsia felis</italic> is an arthropod-associated intracellular, gram-negative bacterium with a worldwide distribution ##REF##9284376##[1]##. Molecular and serological evidence suggests <italic>R. felis</italic> is infectious for humans, causing a disease similar to murine typhus ##REF##8027348##[2]##, although direct evidence of horizontal transmission of <italic>R. felis</italic> from arthropods to humans has not been demonstrated. Several studies have detected <italic>R. felis</italic> DNA in a diverse range of flea species of which <italic>Ctenocephalides felis</italic>, the cat flea, is considered the primary vector ##REF##2240368##[3]##. Laboratory studies have confirmed <italic>R. felis</italic> in <italic>C. felis</italic> populations is predominantly maintained by transstadial and transovarial transmission ##REF##1729713##[4]##, ##REF##12125879##[5]##. Stable vertical transmission of <italic>R. felis</italic> has been examined in colonized fleas and the reported prevalence of <italic>R. felis</italic> in commercial and institutional <italic>C. felis</italic> colonies ranged from 43–100% ##REF##9284376##[1]##, ##REF##12125879##[5]##–##REF##16893625##[7]##. Within an isolated colony, such as the <italic>C. felis</italic> colony at Louisiana State University (LSU), prevalence of <italic>R. felis</italic> is variable with studies reporting a 65% and 100% prevalence of <italic>R. felis</italic> in 1999 ##REF##10925806##[6]## and 2007 ##REF##16893625##[7]##, respectively. Vertical transmission of <italic>R. felis</italic> persists in <italic>C. felis</italic> for at least 12 generations without the aid of an <italic>R. felis</italic>-infected bloodmeal, however, over successive generations prevalence wanes to low levels (&lt;10%) ##REF##12125879##[5]##. The specific mechanisms by which <italic>R. felis</italic> prevalence fluctuates within a <italic>C. felis</italic> colony are unknown, but in nature, the prevalence of <italic>R. felis</italic> in a flea population is likely amplified by fleas feeding on <italic>R. felis</italic>-infected mammalian hosts ##REF##1629332##[8]##.</p> <p>Ultrastructural analyses have been used to characterize <italic>R. felis</italic> within the flea host, as well as in cell lines ##REF##2240368##[3]##, ##REF##16885313##[9]##, ##UREF##0##[10]##. Using electron microscopy, the presence of <italic>R. felis</italic> has been demonstrated in <italic>C. felis</italic> midgut epithelial cells, ovaries, salivary glands, and muscle ##REF##2240368##[3]##, ##UREF##0##[10]##. Additionally, molecular detection of <italic>R. felis</italic> has utilized traditional and quantitative real-time PCR (qPCR), typically targeting a portion of the genus-common 17-kDa antigen gene, for field sample diagnostics ##REF##16893625##[7]##, ##REF##15799525##[11]##, ##REF##12183224##[12]##. While qPCR has been utilized to assess the replication kinetics of <italic>R. felis in vitro</italic>
##REF##12183224##[12]##, <italic>in vivo</italic> (arthropod and vertebrate) models of infection, for the purpose of examining the infection kinetics of <italic>R. felis</italic>, have not been examined. Development of such models will further elucidate the mechanisms for maintenance of <italic>R. felis</italic>, specifically in the flea populations that serve as vectors and reservoirs, and broaden our understanding of the ecology of <italic>R. felis</italic> transmission in nature.</p> <p>Recent analysis of rickettsial replication in a <italic>Rickettsia</italic>/tick model (<italic>Rickettsia amblyommii</italic>/<italic>Amblyomma americanum</italic>) utilizing qPCR demonstrated a generalized dissemination of rickettsiae and, notably, a steady state level of rickettsial infection during relatively long (7–14 day) tick feeding periods ##UREF##1##[13]##. While the direct regulators of rickettsial replication are under investigation in the slow-feeding tick host, the kinetics of rickettsial replication within the flea has not yet been examined. After locating a suitable host, cat fleas immediately begin feeding. Within 24 to 36 hr after initiating bloodfeeding, female cat fleas begin laying eggs and reach peak egg production (40 to 50 eggs/day) in 4 to 9 days ##REF##2588462##[14]##–##REF##9017899##[16]##. Feeding adult fleas can live as long as 100 days, however, fleas that have initiated feeding and are removed from the host usually die within 24 to 48 hr ##REF##9017899##[16]##. Similar to other <italic>Rickettsia</italic>/arthropod relationships, <italic>R. felis</italic> is maintained via vertical transmission in the LSU colony of <italic>C. felis</italic>. We hypothesize that within colonized fleas, the intimate relationship between <italic>R. felis</italic> and <italic>C. felis</italic> allows for coordination of rickettsial replication and the metabolically active periods during flea bloodmeal acquisition and oogenesis. In the present study, the replication kinetics of <italic>R. felis</italic> in the cat-fed, LSU <italic>C. felis</italic> colony were investigated by developing a qPCR assay to detect and quantify <italic>R. felis</italic> in actively feeding cat fleas. Likewise, we considered the influence of flea gender and <italic>R. felis</italic> prevalence within the colony on rickettsial load in cat fleas during bloodmeal acquisition. A sensitive and specific method of detection, qPCR, allowed for accurate detection and quantification of <italic>R. felis</italic> in <italic>C. felis</italic>, including samples with low levels of infection.</p>
<title>Materials and Methods</title> <title>Fleas and cats</title> <p>Unfed <italic>C. felis</italic> were obtained from a colony maintained at LSU (Baton Rouge, LA). For the past 20 years, this colony has been persistently infected at a varying prevalence with <italic>R. felis</italic>
##REF##12125879##[5]##. At the Louisiana State University School of Veterinary Medicine (LSU-SVM) unfed newly-emerged, adult <italic>C. felis</italic> (mixed sex) were allowed to feed on 2 year-old, short-hair, flea-naïve, specific pathogen free cats (Harlan, Indianapolis, IN) as previously described by Foil et al. ##REF##9775583##[17]##. All cats were housed individually with 12∶12 light∶dark cycles at 20–22.2°C and 40–60% relative humidity. Pre-study blood collected from all cats was negative for <italic>Rickettsia</italic> by PCR. Briefly, fleas were contained in a capsule secured to the shaved lateral thorax of a cat and were able to feed through the fine mesh of the capsule. Before each trial (Day 0), colony <italic>R. felis</italic>-infection prevalence, infection load, and infection density of a minimum of nine individual fleas were assessed. Animal use in this research was conducted in accordance with protocols approved by the LSU Institutional Animal Care and Use Committee, and an approved protocol is on file in the office of the Division of Laboratory Animal Medicine at LSU-SVM.</p> <title>Experimental design</title> <p>Three separate trials were conducted over the course of one year. Each trial was initiated approximately six months from the start of the previous trial. All fleas were drawn from the same colony, therefore, fleas in subsequent trials were descendants of the previous trial as the colony was not supplemented with outside fleas. The dynamic <italic>R. felis</italic>-infection prevalence in the LSU <italic>C. felis</italic> colony did not allow trials to be exact replicates of one another.</p> <p>In each trial, capsules containing ∼100, unfed adult (mixed sex) fleas were secured to an individual cat. Approximately ten fleas were randomly selected and removed from the capsule for nine consecutive days. After collection, all flea samples were immediately frozen at −20°C until further processing. Every three days, fleas remaining in the study were transferred to a new capsule and returned to the cat.</p> <title>DNA Isolation</title> <p>Fleas recovered daily from the cat host were individually sexed, assigned sample numbers, placed in microcentrifuge tubes and homogenized with sterile plastic pestles. Genomic DNA (gDNA) was extracted using either QIAGEN DNeasy Tissue Kit (Chatsworth, CA) according to the manufacturer's instructions and eluted in 50 µl Buffer AE (Trial 1), or using a modified version of the HotSHOT DNA extraction protocol (Trials 2 and 3) as described next ##REF##10907076##[18]##. Briefly, individual fleas were placed in microcentrifuge tubes partially submerged in liquid nitrogen and pulverized using sterile plastic pestles. Ground samples were incubated at 95°C for 45 min in 30 µl alkaline lysis reagent (25 mM NaOH, 0.2 mM disodium EDTA, pH of 12), cooled to 4°C for 5 min and mixed with 30 µl of neutralizing reagent (40 mM Tris-HCl, pH of 5). All gDNA preparations were stored at −20°C.</p> <title>Rickettsial Detection and Identification by PCR and Sequencing</title> <p>
<italic>R. felis</italic> infection in flea samples was detected by PCR amplification of a small portion of the <italic>R. felis</italic> 17-kDa antigen gene using primers (<italic>Rf</italic>17.153 <named-content content-type="gene">5′-AGGACAGCTTGTCGGAGTAGG-3′</named-content> and <italic>Rf</italic>17.309 <named-content content-type="gene">5′-ACGCCATTCTACGCTAGTGC-3′</named-content>) derived from the complete genome of <italic>R. felis</italic> available on GenBank (Accession #: CP000053). A portion of <italic>C. felis</italic> 18S rDNA, amplified by primers (<italic>Cf</italic>18S <named-content content-type="gene">5′-TGCTCACCGTTTGACTTGG-3′</named-content> and <italic>Cf</italic>18S <named-content content-type="gene">5′-GTTTCTCAGGCTCCCTCTCC-3′</named-content>) derived from an available partial sequence of <italic>C. felis</italic> 18S rDNA in GenBank (Accession #: AF136859), was used as a control to check the integrity of the template DNA. All primers used for standard PCR and qPCR were synthesized by Integrated DNA Technologies, Inc. (Coralville, IA). PCR products were amplified using isolated gDNA (individual flea lysates) or water (negative control) as template, gene-specific primers and 2× PCR Master Mix (Promega, Madison, WI). PCR conditions were as follows: initial denaturation at 94°C for 3 min, followed by 35 cycles of denaturation at 94°C for 30 s, annealing at 60°C for 45 s, extension at 72°C for 1 min, and a final extension at 72°C for 7 min. Amplified PCR products were visualized by electrophoresis on ethidium bromide stained 2.0% agarose gels. To verify identity, representative PCR products from each trial were cloned into pCR4-TOPO vectors (Invitrogen, Carlsbad, CA) and sequenced as previously described by Pornwiroon et al. ##REF##18043659##[19]##.</p> <title>Southern Blot Analysis</title> <p>A Southern blot was conducted to estimate the <italic>C. felis</italic> 18S rDNA copy number. Southern blot was performed as previously described by Horigane et al. ##REF##17366597##[20]##. Briefly, gDNA was extracted from <italic>Rickettsia</italic>-uninfected, unfed adult cat fleas (Heska Corporation, Loveland, CO) using Qiagen DNeasy DNA Extraction Kit and digested with five restriction enzymes (6 µg gDNA/enzyme): <italic>Eag</italic> I, <italic>Pst I</italic>, <italic>Xba</italic> I, <italic>Xho</italic> I (New England BioLabs, Ipswich, MA), and <italic>Eco</italic>R I (Promega, Madison, WI). DNA was separated on a 0.8% TAE gel and transferred to a positively charged nylon membrane (Amersham Biosciences, Piscataway, NJ) with 20× SSC, and fixed by UV. <italic>C. felis</italic> 18S rDNA was identified using the probe, <italic>Cf</italic>18SP, that was constructed based on a portion of the <italic>C. felis</italic> 18S rDNA sequence and amplified with primers <italic>Cf</italic> 18S<sub>for</sub> and <italic>Cf</italic> 18S<sub>rev</sub>. The probe was labeled with PCR DIG Labeling Mix (Roche, Indianapolis, IN) according to the manufacturer's protocol. Anti-digoxigenin-AP (Fab fragment, Roche) and CDP Star (New England Bio Labs) were used for detection.</p> <title>Construction of Internal-Control Plasmid for Quantitative Real-Time PCR</title> <p>To quantify copies of <italic>R. felis</italic> and <italic>C. felis</italic> genes in flea samples, serial dilutions of a plasmid containing single-copy portions of both the <italic>R. felis</italic> 17-kDa antigen gene and <italic>C. felis</italic> 18S rDNA were used to generate a standard curve. The plasmid was constructed by PCR amplification of a 157-base pair (bp) fragment of <italic>R. felis</italic> 17-kDa antigen gene and a 179-bp fragment of <italic>C. felis</italic> 18S rDNA with <italic>Rf</italic>17.153<sub>for</sub>-<italic>Rf</italic>17.309<sub>rev</sub> and <italic>Cf</italic>18S<sub>for</sub>-<italic>Cf</italic>18S<sub>rev</sub> primers, respectively. Each amplicon was cloned into to a pCR4-TOPO vector and sequenced to confirm identity as previously described ##REF##18043659##[19]##. Both amplicons were PCR amplified using a gene-specific primer and either a M13-forward or M13-reverse primer, digested with <italic>Eco</italic>R I and ligated together. The ligation product was amplified by PCR with primer pair <italic>Rf</italic>17.153<sub>for</sub> and <italic>Cf</italic>18S<sub>rev</sub>, cloned and sequenced. The resulting plasmid, pCR4-TOPO-<italic>Rf</italic>17kDa+<italic>Cf</italic>18SrDNA, served as a standard template. The minimum detection limit for pCR4-TOPO-<italic>Rf</italic>17kDa+<italic>Cf</italic>18SrDNA was 10 copies.</p> <title>Quantitative Real-Time PCR</title> <p>For each gene, qPCR components and template that included 2× iTaq SYBR Green Supermix (BioRad, Hercules, CA); 100 nM of each primer; DNase/RNase-free water; and 5 µl of gDNA template (samples), water (negative control), or serial 10-fold dilutions (1×10<sup>8</sup> to 10 copies) of pCR4-TOPO-<italic>Rf</italic>17kDa+<italic>Cf</italic>18SrDNA were pre-mixed in 35 µl volumes in 96-well plates and aliquoted in triplicate 10 µl reactions on 384-well plates. The qPCR was performed with an ABI 7900HT unit (Applied Biosystems, Foster City, CA) at LSU-SVM using conditions previously described ##UREF##1##[13]##. Results were analyzed with ABI 7900HT sequence detection system (SDS v2.3) software. The specificity of the assay was verified and the expected single peak for the internal-control plasmid and positive gDNA samples, but not in the water (negative control) samples, was identified in the dissociation curve. Additionally, representative qPCR products from each trial were verified by gel analysis to confirm the specificity of the reaction and cloned and sequenced to confirm identity (data not shown). Trial <italic>R. felis</italic>-infection prevalence in the <italic>C. felis</italic> colony was quantified as the percent of fleas positive for <italic>R. felis</italic> infection out of the total number examined per trial. <italic>R. felis</italic>-infection load was quantified as the copy number of <italic>Rf</italic>17kDa per individual flea lysate. <italic>R. felis</italic>-infection density was quantified as the ratio of log transformed <italic>Rf</italic>17kDa and log transformed <italic>Cf</italic>18S rDNA copy numbers (<italic>Rf</italic>17kDa/<italic>Cf</italic>18S) per individual flea.</p> <title>Statistical Analysis</title> <p>Rickettsial load in fleas and the ratio of <italic>Rf</italic>17kDa/<italic>Cf</italic>18S were assessed after the logarithmic transformation of the quantity of the genes of interest (<italic>Rf</italic>17kDa and <italic>Cf</italic>18S). Analysis of variance, (SAS statistical package, Version 9.1.3, GLM procedure ANOVA, Cary, NC) was performed to examine potential differences between rickettsial load in fleas and ratio of <italic>Rf</italic>17kDa/<italic>Cf</italic>18S copy number over the study period; when overall significance was found, Tukey's honestly significant difference (HSD) post hoc test was used to examine pairwise differences of means of main effects. Pairwise <italic>t</italic>-tests of least square means were performed to determine any interaction effects between trial, gender, and experimental day for rickettsial infection load and ratio of <italic>Rf</italic>17kDa/<italic>Cf</italic>18S. An <italic>F</italic>-test was used for general comparisons of grouped means. For all comparisons, a <italic>P</italic>-value of &lt;0.05 was considered significantly different.</p>
<title>Results</title> <title>Prevalence of <italic>R. felis</italic> in LSU <italic>C. felis</italic> Colony</title> <p>The specificity of <italic>R. felis</italic>-infection of <italic>C. felis</italic> (LSU colony) was confirmed by sequencing a portion of the 17-kDa antigen gene from a representative subset of fleas positive for <italic>Rickettsia spp</italic>. infection. Obtained sequences were compared to those in GenBank and all samples sequenced demonstrated a 100% identity to <italic>R. felis</italic> (accession number CP000053).</p> <p>The prevalence of <italic>R. felis</italic>-infection in individual <italic>C. felis</italic> lysates was assessed by qPCR using primers that amplified a portion of the 17-kDa antigen gene. Trial prevalence and the daily prevalence range for each trial are listed in ##TAB##0##Table 1##. <italic>R. felis</italic>-infection prevalence was greatest during Trial 1 and decreased 23% and 61% in Trials 2 and 3, respectively. When prevalence was assessed separately for male and female fleas, <italic>R. felis</italic> infection was gender independent.</p> <title>Quantification of <italic>Rf17</italic>kDa and <italic>Cf</italic>18S Gene Copies During Flea Feeding</title> <p>To determine <italic>R. felis</italic>-infection load, fleas were individually titrated and the total number of <italic>R. felis</italic> (17-kDa copy number) in each individual flea sample was determined by qPCR. Serial dilutions of pCR4-TOPO-<italic>Rf</italic>17kDa+<italic>Cf</italic>18SrDNA were used to generate a standard curve and extrapolate <italic>R. felis</italic> quantities per individual flea lysate (##SUPPL##0##Table S1a##). The range of detected <italic>Rf</italic>17kDa was 3.1×10<sup>1</sup> to 3.74×10<sup>5</sup> copies per reaction, corresponding to a 1.3×10<sup>3</sup> to 1.6×10<sup>7</sup> total <italic>Rf</italic>17kDa load per flea lysate. The mean quantities of <italic>R. felis</italic> at each time point for individual flea samples in each trial are presented in ##FIG##0##Figure 1##. Among all trials, fleas had significantly different <italic>R. felis</italic>-infection loads; increasing 4.75-fold during the course of the study between fleas in Trial 1 and Trial 3. Trial 3 fleas had the greatest <italic>R. felis</italic>-infection loads (mean load±SEM across all time points was 7.88×10<sup>6</sup>±5.43×10<sup>5</sup>) followed by Trial 2 (2.54×10<sup>6</sup>±1.40×10<sup>5</sup>) and Trial 1 (1.66×10<sup>6</sup>±1.62×10<sup>5</sup>). Combining data from all trials, when flea gender was considered, female fleas had significantly greater mean <italic>R. felis</italic>-infection loads than males. However, at individual time points within a trial or across trials there were few significant differences in <italic>R. felis</italic>-infection loads between male and female fleas.</p> <p>When the mean daily <italic>R. felis</italic>-infection loads for each trial, and for all trials combined, were compared, significant variability in <italic>R. felis</italic>-infection load was observed by day compared to the starting point at Day 0 within Trials 1 and 2 (##FIG##0##Figure 1##). No significant variability in daily <italic>R. felis</italic>-infection load was observed in Trial 3, when compared to the Day 0 infection load. Analyses of the mean daily <italic>R. felis</italic> infection load combining all Trials demonstrated that there was no single day when <italic>R. felis</italic>-infection load was consistently greater within the study period. Likewise, although variability in <italic>R. felis</italic>-infection load was observed in some trials, there was no consistent pattern of increased or decreased rickettsial infection across all trials. Therefore, a consistent <italic>R. felis</italic>-infection load during flea feeding was observed as there was no significant correlation between <italic>R. felis</italic> replication and flea feeding or oviposition events.</p> <p>To estimate flea sample size, <italic>R. felis</italic>-infected and -uninfected <italic>C. felis</italic> were individually titrated and the <italic>Cf</italic>18S rDNA copy number within a single flea was determined by qPCR. Serial dilutions of pCR4-TOPO-<italic>Rf</italic>17kDa+<italic>Cf</italic>18SrDNA were used to generate a standard curve and extrapolate <italic>Cf</italic>18S rDNA copy number per individual <italic>C. felis</italic> lysate (##SUPPL##0##Table S1b##). Within individual trials, mean <italic>Cf</italic>18S copy numbers did significantly vary during the feeding period on some days, however there was no single day where <italic>Cf</italic>18S copy numbers were consistently different. For both <italic>R. felis</italic>-infected and -uninfected fleas, the mean <italic>Cf</italic>18S copy numbers increased with each subsequent Trial, with a significant 2.6-fold increase in mean <italic>Cf</italic>18S copy numbers between Trials 2 and 3. When data from all three trials were combined, significantly greater <italic>Cf</italic>18S copy numbers were observed in female <italic>C. felis</italic> compared to male <italic>C. felis</italic>. Likewise, uninfected <italic>C. felis</italic> had significantly greater mean <italic>Cf</italic>18S copy numbers compared to <italic>R. felis</italic>-infected <italic>C. felis</italic>.</p> <title>
<italic>R. felis</italic>-infection Load in <italic>C. felis</italic> Expressed as a Ratio of <italic>Rf</italic>17kDa/<italic>Cf</italic>18S</title> <p>
<italic>R. felis</italic>-infection of <italic>C. felis</italic> was further assessed as a ratio of <italic>R. felis</italic> and <italic>C. felis</italic> genes. Examination of the published genome demonstrated that the <italic>R. felis</italic> 17-kDa antigen gene is a single-copy gene ##REF##15984913##[21]##. In this study, the copy number of <italic>C. felis</italic> 18S rDNA was estimated by Southern blot analysis. <italic>C. felis</italic> gDNA was digested with 5 restriction enzymes all resulting in the presence of a single band of variable size (dependant on enzyme used) (##FIG##1##Figure 2##). Although each enzyme produced a single band indicating that the <italic>C. felis</italic> 18S rDNA gene may be single-copy, a tandem arrangement of multiple gene copies may not be discernable using Southern blot analysis alone.</p> <p>For individual flea lysates positive for <italic>R. felis</italic>, <italic>R. felis</italic>-infection density was examined. Assessed by qPCR, copies of <italic>Rf</italic>17kDa and <italic>Cf</italic>18S were determined relative to their positions in the standard curve and extrapolated to the individual flea lysate. To determine <italic>R. felis</italic>-infection density <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios for individual flea samples were calculated after the logarithmic transformation of the quantity of the genes of interest (<italic>Rf</italic>17kDa and <italic>Cf</italic>18S). The mean daily <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios for each trial are presented in ##FIG##2##Figure 3##.</p> <p>The mean (±SEM) trial <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratio increased significantly between subsequent trials from 0.72±0.0064 to 0.75±0.0059 to 0.78±0.0032, in Trials 1–3, respectively. When flea gender was considered, female fleas in Trials 2 and 3 had significantly greater <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios than female fleas in Trial 1, but were not significantly different from one another. Male fleas in Trials 1 and 2 did not have significantly different <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios from one another, but both had significantly lower ratios than male fleas in Trial 3. Comparing female and male fleas within trials, males had significantly greater <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios in Trial 1 than females, but in Trial 2 females had significantly greater <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios than males. Due to the limited occurrence of <italic>R. felis</italic> infection in male fleas in Trial 3, gender-based ratio differences could not be statistically assessed.</p> <p>In individual trials, significant differences in <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios were observed on multiple days, however, among all trials no single day was consistently significantly different from Day 0. The mean daily <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios for all trials were combined and compared, independent of flea gender, and no significant variability in <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratio was observed by day compared to the starting point at Day 0. Analyses of the mean daily <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratio for female fleas, combining and comparing all trials, identified a significantly lower ratio on Day 6 compared to Day 0, a trend similar to that observed for females in Trial 1, but not females in Trials 2 and 3. No significant differences in mean daily <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios for all trials combined were observed for male fleas on any day compared. Additionally, a correlation between the ratio of <italic>Rf</italic>17kDa/<italic>Cf</italic>18S and flea bloodmeal acquisition was not observed. For female fleas, during the time of peak oviposition (Day 6) there was a decrease in <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratio. No distinct relationship was observed between <italic>R. felis</italic>-infection load and <italic>C. felis</italic> 18S copy number. A consistent infection rate of <italic>R. felis</italic>, with minimal significant increases or decreases in infection load is maintained during flea metabolically active periods.</p>
<title>Discussion</title> <p>Both traditional and real-time PCR have been used to further characterize rickettsial infection in fleas. Amplification of a portion of the gene encoding the rickettsial-common 17-kDa antigen allowed for identification and subsequent differentiation by restriction fragment length polymorphism of both <italic>Rickettsia typhi</italic> and <italic>R. felis</italic> in <italic>C. felis</italic>
##REF##9701920##[22]##. In this study, the presence of <italic>R. felis</italic> as the rickettsial species infecting the LSU <italic>C. felis</italic> flea colony was confirmed by traditional PCR and gene sequence analysis. Real-time PCR has also been used for quantitative assessment of <italic>R. felis</italic> replication and as a diagnostic tool to assess <italic>R. felis</italic> infection in arthropod hosts. Replication of four isolates of <italic>R. felis</italic>, including the LSU strain, was measured in XTC-2 cells treated with antibiotics; however, the techniques for establishing rickettsial gene copy numbers were not provided and the qPCR assay utilized is not clear ##REF##12183224##[12]##. More recently, a probe-based qPCR targeting rickettsial <italic>ompB</italic> for detection and differentiation of <italic>R. felis</italic> and <italic>R. typhi</italic> in colonized and wild-caught <italic>C. felis</italic> was described ##REF##16893625##[7]##. It was demonstrated that crude extraction (boiled flea lysate) of DNA was an insensitive procedure, compared to kit-based DNA extraction, resulting in limited detection of <italic>Rickettsia</italic>. In the present report, we utilized qPCR to analyze the rickettsial load in individual <italic>C. felis</italic> during flea feeding and oviposition. Differences in sensitivity between kit-based and HotSHOT DNA recovery techniques utilized in the current study were undetectable. Also, Henry <italic>et al</italic>
##REF##16893625##[7]## utilized serial dilutions of a plasmid containing a portion of the <italic>ompB</italic> target sequence and determined the sensitivity of the assay to be 1 copy/µl; 10 times greater sensitivity than the 17-kDa antigen gene target sequence and SYBR Green assay in the current study. Although there is increased sensitivity in the probe-based assay compared to our SYBR Green assay, the large rickettsial load in LSU fleas does not require detection of low numbers of <italic>R. felis</italic>. Examining gene copy numbers by qPCR is limited in differentiating between live/dead organisms; however, DNA analysis can still provide a reasonable assessment of infection load in hosts that are naturally infected and/or utilized in the vertical transmission of the organism. This study provides the first application of qPCR to examine the kinetics of <italic>R. felis</italic> infection within its anautogenous vector and primary reservoir, the cat flea.</p> <p>Wide dissemination of vertically maintained bacteria is common in arthropods ##REF##10196738##[23]##, ##REF##10759311##[24]##. Within the flea, <italic>R. feli</italic>s infects many types of tissues, including the salivary glands ##REF##2240368##[3]##, ##UREF##0##[10]##; in this study, <italic>R. felis</italic> infection load was examined at the whole individual flea level. Among the <italic>R. felis</italic>-infected fleas, we did identify a mean of 3.9×10<sup>6</sup> rickettsiae during feeding. All trials were significantly different in their mean <italic>R. felis</italic> infection load, with Trial 1 fleas and Trial 3 fleas having the lowest and highest individual flea rickettsial loads, respectively. Regardless of whether infection load was quantified as total load per flea or as a ratio of <italic>R. felis</italic> and <italic>C. felis</italic> genes, rickettsial load was not definitively affected by flea bloodmeal acquisition or oogenesis, nor was there a consistent <italic>R. felis</italic> replication pattern observed across all trials. The significant decrease in the mean trial <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratio on Day 6 is skewed by the results of Trial 1, in which more fleas were positive for <italic>R. felis</italic>-infection. We have utilized qPCR to examine spotted fever group rickettsiae within the tick <italic>A. americanum</italic> and, similar to this study, we identified a steady level of <italic>R. amblyommii</italic> load associated with bloodmeal acquisition when infection was assessed on an individual tissue basis ##UREF##1##[13]##. Therefore, in the two models examined, there appears to be a balance between rickettsial load and host size in actively feeding arthropods. Future studies examining <italic>R. felis</italic> infection within its flea vector and determination of infection density within specific tissues (<italic>e.g.</italic> salivary glands and ovaries) will further elucidate the biological interactions between rickettsiae and flea hosts.</p> <p>In addition to <italic>C. felis</italic>, <italic>R. felis</italic> has been molecularly identified in numerous species of wild-caught fleas, including <italic>C. canis</italic>
##REF##12860622##[25]##, <italic>Xenopsylla cheopis</italic>
##REF##16965716##[26]##, <italic>Archaeopsylla erinacei</italic>
##REF##16606979##[27]##, <italic>Spilopsyllus cuniculi</italic>
##REF##16704850##[28]##, <italic>Echidnophaga gallinacean</italic>
##REF##16704850##[28]##, <italic>Anomiopsylla nudata</italic>
##REF##15799525##[11]##; and, while it is likely that the prevalence of <italic>R. felis</italic> in a flea population is amplified by an infectious bloodmeal, the acquisition mechanisms and stability of transmission are not known. The prevalence of <italic>R. felis</italic> in wild-caught fleas has been reported to be 1–20%, which is typically lower than that observed in colonized fleas ##REF##9284376##[1]##, ##REF##12653300##[29]##–##REF##12023908##[32]##. The efficient replication and vertical transmission of <italic>R. felis</italic> within fleas may minimize the necessity of frequent mammalian infections to maintain <italic>R. felis</italic> in nature; however, experimental evidence indicates the transmission cycle of <italic>R. felis</italic> is more complex. In laboratory reared flea colonies, such as the LSU <italic>C. felis</italic> colony, the presence and dynamic prevalence of <italic>R. felis</italic> make it a valuable tool to examine <italic>Rickettsia</italic>/flea interactions ##REF##18043659##[19]##, ##REF##12023908##[32]##. Wedincamp and Foil ##REF##12125879##[5]## investigated the efficiency of <italic>R. felis</italic> vertical transmission in <italic>C. felis</italic> without the aid of an infectious bloodmeal and demonstrated fluctuating, but decreasing prevalence through twelve generations. During the last 15 years, several independent studies have identified the variable prevalence of <italic>R. felis</italic>, ranging from 43–100% ##REF##16893625##[7]##, ##REF##18043659##[19]##, ##REF##9701920##[22]##, ##REF##8069413##[33]##. The LSU <italic>C. felis</italic> colony is maintained solely on cat hosts, but the role of cat hosts as a source of infectious bloodmeal is unclear. Cats continuously fed on by <italic>R. felis</italic>-infected fleas seroconvert two to four months post-exposure ##REF##10925806##[6]##. <italic>R. felis</italic> DNA has been detected in cat blood ##REF##10925806##[6]##; however, recovery of <italic>R. felis</italic> from these cats has been unsuccessful thus far [K. R. Macaluso, unpublished data]. Each trial utilized a flea-naïve cat host, minimizing the possibility of cat-derived, immune-mediated rickettsial clearance in feeding fleas. Interestingly, even when fleas are fed only on live hosts, there is variance in <italic>R. felis</italic> prevalence between generations as demonstrated in the current study that utilized subsequent generations of fleas for each trial. After a substantial population loss followed by a population expansion in the LSU <italic>C. felis</italic> colony, the prevalence of <italic>R. felis</italic> in the flea population neared 100% ##REF##18043659##[19]##. In individual fleas, 3.1×10<sup>1</sup> to 3.74×10<sup>5</sup> of <italic>Rf</italic>17kDa copies per reaction were observed, with only four fleas having detectable <italic>Rf</italic>17kDa copy numbers (in total lysate) under 1.0×10<sup>4</sup>. Although not proven in this study, the dramatic contrast of <italic>Rf</italic>17kDa copy numbers of these four fleas in comparison to the remaining <italic>R. felis</italic>-infected fleas may suggest a role for low-level horizontal acquisition of <italic>R. felis</italic> in these few fleas from feeding on a shared host or larval cannibalism. Under laboratory conditions with either natural or artificial hosts, the mechanisms of prevalence and infection load fluctuations are intriguing and require further study to assess if carriage of <italic>R. felis</italic> is beneficial to <italic>C. felis</italic>.</p> <p>Alternatively, there may be microbial-dependent influence on <italic>R. felis</italic> prevalence, as we have identified intracellular <italic>Wolbachia</italic> spp. in the fleas ##REF##18043659##[19]##. <italic>Wolbachia</italic> spp. infect many arthropods and are readily able to manipulate their arthropod host (<italic>e.g.</italic> feminization, cytoplasmic incompatibility) with a potential to impact host fitness ##REF##10196738##[23]##, ##REF##15305172##[34]##, ##REF##15012323##[35]##. While <italic>Wolbachia</italic> spp. have been identified in fleas, their impact on flea fitness and relationship with other flea microbiota, such as <italic>R. felis</italic>, have not been examined ##REF##18043659##[19]##, ##REF##12760642##[36]##. Not only can microbiota compete with one another for host resources, they can often manipulate their environment affecting arthropod-host fitness ##REF##15305172##[34]##, ##REF##15012323##[35]##. For example, vertical transmission of the tick-borne human rickettsial pathogens, <italic>R. rickettsii</italic>, <italic>R. parkeri</italic>, and <italic>R. conorii</italic>, are associated with decreased arthropod-host fitness, whereas vertical maintenance of other <italic>Rickettsia</italic> spp. do not impact tick fitness ##UREF##1##[13]##, ##REF##12653135##[37]##, ##REF##12495176##[38]##. Interaction and possible competition of vertically transmitted microbiota and their potential impact on host fitness is likely complex and needs to be scrutinized. Utilization of qPCR and laboratory models of <italic>R. felis</italic> transmission will further elucidate mechanisms of transmission in nature.</p> <p>Although no clear replication pattern was observed in actively bloodfeeding and ovipositing fleas, there was an inverse correlation between colony <italic>R. felis</italic>-infection prevalence and <italic>R. felis</italic>-infection load in individual fleas. Specifically, as the prevalence of <italic>R. felis</italic> decreased in our flea colony from 96% in Trial 1 to 35% in Trial 3, the mean <italic>R. felis</italic>-infection load in individual fleas increased 4.75-fold (##FIG##3##Figure 4##). Quantification of increasing <italic>R. felis</italic>-infection was demonstrated at the whole flea level (counting <italic>Rf</italic>17kDa copy numbers) and verified again when assessing infection as a ratio of <italic>Rickettsia</italic> and flea genes. <italic>C. felis</italic> 18S rDNA copies also increased in both <italic>R. felis</italic>-infected and uninfected fleas across all trials. The ratio of <italic>R. felis</italic> and <italic>C. felis</italic> genes increases significantly from Trial 1 to Trial 3, demonstrating that the larger fleas (greater <italic>Cf</italic>18S quantities) in Trial 3 have an unproportionately greater <italic>R. felis</italic> burden than the smaller fleas in Trial 1. These results verify our <italic>R. felis</italic>-infection load per individual flea lysate findings, supporting that infection load is increasing beyond the expected proportional increase relative to <italic>C. felis</italic> 18S rDNA copy number. Whether or not larger, more densely infected fleas are more competent vectors for <italic>R. felis</italic> will be interesting to explore further.</p> <p>Variable prevalence and infection density in <italic>C. felis</italic> may represent a <italic>R. felis</italic> ecological maintenance strategy, in which waning prevalence signals, by an unknown mechanism, increased infection burdens in individual fleas potentially facilitating more efficient transmission to progeny (vertical transmission) or a reservoir host (horizontal transmission) in order to persist in the flea population. Here we describe the first quantitative assessment of <italic>R. felis</italic> infection in fleas. We also have data to suggest that the previously unrecognized horizontal transmission of <italic>R. felis</italic> occurs among fleas and that prevalence of <italic>R. felis</italic> in the population is correlated to individual flea rickettsial load. The results of this study will help elucidate the epidemiology of <italic>R. felis</italic>, an emerging infectious pathogen, by demonstrating the dynamic prevalence and infection density changes that occur within a flea population.</p>
<p>Conceived and designed the experiments: KER GCH LDF KRM. Performed the experiments: KER RWS GCH KRM. Analyzed the data: KER. Wrote the paper: KER RWS LDF KRM.</p> <title>Background</title> <p>
<italic>Rickettsia felis</italic> is a flea-associated rickettsial pathogen recurrently identified in both colonized and wild-caught cat fleas, <italic>Ctenocephalides felis</italic>. We hypothesized that within colonized fleas, the intimate relationship between <italic>R. felis</italic> and <italic>C. felis</italic> allows for the coordination of rickettsial replication and metabolically active periods during flea bloodmeal acquisition and oogenesis.</p> <title>Methodology/Principal Findings</title> <p>A quantitative real-time PCR assay was developed to quantify <italic>R. felis</italic> in actively feeding <italic>R. felis</italic>-infected fleas. In three separate trials, fleas were allowed to feed on cats, and a mean of 3.9×10<sup>6</sup>
<italic>R. felis</italic> 17-kDa gene copies was detected for each flea. A distinct <italic>R. felis</italic> infection pattern was not observed in fleas during nine consecutive days of bloodfeeding. However, an inverse correlation between the prevalence of <italic>R. felis</italic>-infection, which ranged from 96% in Trial 1 to 35% in Trial 3, and the <italic>R. felis</italic>-infection load in individual fleas was identified. Expression of <italic>R. felis</italic>-infection load as a ratio of <italic>R. felis</italic>/<italic>C. felis</italic> genes confirmed that fleas in Trial 3 had significantly greater rickettsial loads than those in Trial 1.</p> <title>Conclusion/Significance</title> <p>Examining rickettsial infection dynamics in the flea vector will further elucidate the intimate relationship between <italic>R. felis</italic> and <italic>C. felis</italic>, and facilitate a more accurate understanding of the ecology and epidemiology of <italic>R. felis</italic> transmission in nature.</p>
<title>Supporting Information</title>
<p>We thank Mark Guillotte for technical assistance with the LSU flea colony and Mike Kearney for his help with statistical analysis. We also thank Ms. Julie Millard for assistance with monitoring cat vitals.</p>
<fig id="pone-0002805-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002805.g001</object-id><label>Figure 1</label><caption><title>
<italic>R. felis</italic>-infection load in individual <italic>C. felis</italic> lysates.</title><p>For nine consecutive days, 10 fleas were randomly selected and removed from the feeding capsule situated on a cat host. In <italic>Rickettsia</italic> positive fleas, qPCR was used to quantify the mean (±SEM) <italic>R. felis</italic>-infection load. Results are presented as the mean daily <italic>R. felis</italic>-infection load (17-kDa antigen gene copy number) of individual, whole flea lysates during flea bloodmeal acquisition. Within each trial significant differences from Day 0 (unfed fleas) are marked with an asterisk. <sup>Ψ</sup>In Trial 3, Day 3 only one flea was positive for <italic>R. felis</italic> infection.</p></caption></fig> <fig id="pone-0002805-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002805.g002</object-id><label>Figure 2</label><caption><title>Estimation of <italic>C. felis</italic> 18S rDNA copy number using Southern blot.</title><p>
<italic>C. felis</italic> genomic DNA (6 µg/enzyme) was digested with <italic>Eag</italic> I, <italic>Eco</italic>R I, <italic>Pst</italic> I, <italic>Xba</italic> I and <italic>Xho</italic> I. Uncut gDNA served as a negative control and PCR product of a portion of <italic>C. felis</italic> 18S rDNA served as a positive control. Genomic DNA was hybridized with <italic>Cf</italic>18S bp probe to estimate the number of 18S rDNA gene copies in <italic>C. felis</italic>. Digestion with each enzyme results in a single digestion product. (*DNA was not completely digested by <italic>Pst</italic> I. Top band is uncut DNA.)</p></caption></fig> <fig id="pone-0002805-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002805.g003</object-id><label>Figure 3</label><caption><title>
<italic>R. felis</italic>-infection load in <italic>C. felis</italic> expressed as a ratio of <italic>R. felis</italic> and <italic>C. felis</italic> genes.</title><p>For nine consecutive days 10 fleas were randomly selected and removed from the feeding capsule situated on a cat host. In <italic>Rickettsia</italic> positive fleas, <italic>R. felis</italic> infection was determined by quantifying the <italic>R. felis</italic> 17-kDa copy number. For <italic>R. felis</italic>-infected fleas, <italic>C. felis</italic> 18S rDNA copy number was quantified to serve as a comparison point for <italic>R. felis</italic> infection. <italic>R. felis</italic>-infection density was determined by logarithmically transforming and taking the ratio of <italic>R. felis</italic> 17-kDa and <italic>C. felis</italic> 18S rDNA (<italic>Rf</italic>17kDa/<italic>Cf</italic>18S) copy numbers and the mean (±SEM) daily <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratios for each trial were calculated. Within each trial significant differences from Day 0 (unfed fleas) are marked with an asterisk. <sup>Ψ</sup>In Trial 3, Day 3 only one flea was positive for <italic>R. felis</italic> infection. All trials have significantly higher ratios than the previous trial, indicating that rickettsial burdens are increasing in the LSU <italic>C. felis</italic> colony.</p></caption></fig> <fig id="pone-0002805-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002805.g004</object-id><label>Figure 4</label><caption><title>Experimental model depicting the relationship between <italic>R. felis</italic> and <italic>C. felis</italic>.</title><p>Supported by the results of this study, <italic>R. felis</italic> prevalence and individual flea infection load in the LSU <italic>C. felis</italic> colony are inversely correlated. The ratio of <italic>Rf</italic>17kDa/<italic>Cf</italic>18SrDNA also increased significantly between Trial 1 and Trial 3 indicating fleas are infected at a greater density. Trials are situated within the model according to their individual results. As a population, fleas in Trial 1 had the highest prevalence of <italic>R. felis</italic> infection and the lowest mean individual <italic>R. felis</italic>-infection load. Conversely, fleas in Trial 3 had the lowest prevalence of <italic>R. felis</italic>-infection and the greatest mean individual <italic>R. felis</italic>-infection load. Trial 2 fleas represent a median demonstrating the progression of decreasing colony prevalence and increasing infection load from Trial 1 to Trial 3. The ratio of <italic>R. felis</italic> and <italic>C. felis</italic> genes increases significantly from Trial 1 to Trial 3, demonstrating that fleas in Trial 3 have a greater <italic>R. felis</italic> burden (higher <italic>Rf</italic>17kDa/<italic>Cf</italic>18S ratio) than fleas in Trial 1. These results indicate that at increased infection loads, <italic>R. felis</italic> may influence flea fitness to facilitate their own successful transmission to the next generation of fleas or to a susceptible mammalian host.</p></caption></fig>
<table-wrap id="pone-0002805-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002805.t001</object-id><label>Table 1</label><caption><title>Trial <italic>R. felis</italic>-infection prevalence and daily range of cat fleas positive for <italic>R. felis</italic>-infection.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Trial</td><td align="left" rowspan="1" colspan="1">Trial prevalence</td><td align="left" rowspan="1" colspan="1">Daily range of cat fleas positive for <italic>R. felis</italic> infection</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">96%</td><td align="left" rowspan="1" colspan="1">89–100%</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">73%</td><td align="left" rowspan="1" colspan="1">40–90%</td></tr><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">35%</td><td align="left" rowspan="1" colspan="1">10–60%</td></tr><tr><td align="left" rowspan="1" colspan="1">Mean of Trials</td><td align="left" rowspan="1" colspan="1">68%</td><td align="left" rowspan="1" colspan="1">46–83%</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0002805.s001"><label>Table S1</label><caption><p>Raw data for gene copy numbers of (a) <italic>R. felis</italic> 17-kDa and (b) <italic>C. felis</italic> 18S rDNA in individual flea samples. In three independent trials spanning one year, individual fleas were assessed by qPCR for <italic>R. felis</italic> infection prior to the start of the feeding (Day 0) and during nine days of feeding on cat hosts. (a) <italic>R. felis</italic> 17-kDa gene copy number quantities per individual flea lysate were extrapolated based on their position relative to the standard curve and listed in the above table. Individual fleas negative for infection with <italic>R. felis</italic> are represented by ‘-’. (b) <italic>C. felis</italic> 18S rDNA gene copy number quantities per individual flea lysate were also determined by extrapolating their values based on their position relative to the standard curve and are listed in the above table. <italic>C. felis</italic> 18S rDNA counts were performed for all fleas, however, fleas determined to be infected with <italic>R. felis</italic> are shaded in gray. The values in the above tables were log transformed and ratios of log<italic>Rf</italic>17kDa/log<italic>Cf</italic>18SrDNA were generated for individual fleas samples (results presented in main text).</p><p>(0.04 MB PDF)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This research was supported by Louisiana Board of Regents (LEQSF), the National Institutes of Health (NIH; P20RR0201595), and the National Institute of Allergy and Infectious Diseases (R21AI70705). Gretchen Henry was supported by an NIH T35 (RR017504-01) awarded to LSU-SVM.</p></fn></fn-group>
<graphic id="pone-0002805-t001-1" xlink:href="pone.0002805.t001"/> <graphic xlink:href="pone.0002805.g001"/> <graphic xlink:href="pone.0002805.g002"/> <graphic xlink:href="pone.0002805.g003"/> <graphic xlink:href="pone.0002805.g004"/>
<media xlink:href="pone.0002805.s001.pdf"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["10"], "element-citation": ["\n"], "surname": ["Macaluso", "Pornwiroon", "Popov", "Foil"], "given-names": ["KR", "W", "VL", "LD"], "year": ["2008"], "article-title": ["Identification of "], "italic": ["Rickettsia felis"], "source": ["Vector Borne Zoonotic Dis"], "volume": ["8"]}, {"label": ["13"], "element-citation": ["\n"], "surname": ["Zanetti", "Pornwiroon", "Kearney", "Macaluso"], "given-names": ["AS", "W", "MT", "KR"], "year": ["2008"], "article-title": ["Characterization of rickettsial infection in "], "italic": ["Amblyomma americanum"], "source": ["J Med Entomol"], "volume": ["45"]}]
PLoS One. 2008 Jul 30; 3(7):e2805