Abstract
Rickettsia typhi, transmitted by rat fleas, causes most human flea-borne rickettsioses worldwide. Another rickettsia, Rickettsia felis, found in cat fleas, Ctenocephalides felis, has also been implicated as a potential human pathogen. In the continental United States, human cases of flea-borne rickettsioses are reported primarily from the southern regions of Texas and California where the cat flea is considered the principal vector. In California, more than 90% of locally acquired human cases are reported from suburban communities within Los Angeles and Orange counties despite the almost ubiquitous presence of cat fleas and their hosts throughout the state. The objective of this study is to assess the presence and infection rate of Rickettsia species in cat fleas from selected endemic and nonendemic regions of California. Cat fleas were collected from cats in Los Angeles County (endemic region) and Sacramento and Contra Costa counties (nonendemic region). Sequencing of 17 amplicons confirmed the presence of R. felis in both the endemic and non-endemic regions with a calculated maximum likelihood estimation of 131 and 234 per 1000 fleas, respectively. R. typhi was not detected in any flea pools. Two R. felis-like genotypes were also detected in fleas from Los Angeles County; Genotype 1 was detected in 1 flea pool and Genotype 2 was found in 10 flea pools. Genotype 1 was also detected in a single flea pool from Sacramento County. Results from this study show that R. felis is widespread in cat flea populations in both flea-borne rickettsioses endemic and nonendemic regions of California, suggesting that a high prevalence of this bacterium in cat fleas does not predispose to increased risk of human infection. Further studies are needed to elucidate the role of R. felis and the two R. felis-like organisms as etiologic agents of human flea-borne rickettsioses in California.
Introduction
R
In the continental United States, human cases of flea-borne rickettsioses are reported from the southern portions of Texas and California where cat fleas and opossums have been implicated as the principal vector and potential reservoirs, respectively (Adams et al. 1970, Civen and Ngo 2008). A limited number of epidemiological studies have been performed to assess the overall prevalence of R. felis and R. typhi in cat fleas in Southern California. Rickettsial DNA was identified by polymerase chain reaction (PCR) in 52 (31.1%) of 167 cat fleas collected from opossums in Los Angeles County, 12 of which were further characterized as R. typhi (9 fleas) and R. felis (3 fleas) (Williams et al. 1992). Eremeeva et al. (2012) detected R. felis DNA in 882 (47.1%), R. typhi DNA in 25 (1.3%), and DNA of both organisms in 32 (1.7%) of 1873 cat flea pools recovered from cats, roof rats, and opossums trapped within the vicinity of households of human cases in Los Angeles and Orange counties. Rickettsia DNA was not recovered from the blood of any of the host animals nor were animals seropositive for Rickettsia species. R. felis was detected, however, in other tissues of four opossums. In the neighboring San Bernardino County, only R. felis DNA was detected in 152 (26.7%) of 570 cat flea pools and in 26 (22.6%) of 115 individual fleas recovered from opossums (Abramowicz et al. 2012).
In California, human cases of flea-borne rickettsioses are reported to public health authorities as “Rickettsial Diseases” (non-Rocky Mountain spotted fever), including typhus or typhus-like illness. Between 2010 and 2014, an average of 80 such cases was reported annually. Cat fleas and their hosts are ubiquitous in California, yet locally acquired human cases of flea-borne rickettsioses are reported almost exclusively from suburban communities within Los Angeles and Orange counties (Civen and Ngo 2008, Eremeeva et al. 2012). The reasons for the apparent regional endemicity of disease in California are not known, but could be due to differences in the presence, infection rate, and/or species of rickettsia in populations of cat fleas from other areas of the state. Surveillance bias due to physician awareness of the disease could also be a contributing factor. Documenting the distribution of R. typhi and R. felis in the presumed cat flea vector in California is necessary to elucidate the epidemiology of flea-borne rickettsioses and the apparent regional distribution of human disease. This study provides an initial comparison of the presence and species composition of rickettsiae in cat fleas from three counties of the state, comprising both endemic and nonendemic regions.
Materials and Methods
Flea collection
Cat fleas were collected from domestic cats in three counties of California. Collection sites included two animal shelters in Los Angeles County (flea-borne rickettsioses endemic region) and two trap/neuter/release (TNR) programs, one in Sacramento County and one in Contra Costa County (flea-borne rickettsioses nonendemic region). One to ten fleas were collected from each cat using a fine-tooth comb and were placed directly into a microcentrifuge tube (one per cat) containing 70% ethanol (Carolina Biological Supply Co.). All fleas in each tube made up a flea pool.
Molecular analyses
Pooled fleas from each cat were prepared and underwent molecular analyses. Fleas were washed in molecular-grade water to remove alcohol and opened midsection using disposable 18G to 20G needles. Flea tissues were vortex mixed with 180 μL of lysis buffer and 20 μL of proteinase K and DNA extracted following the manufacturer's instructions (QIAmp DNA Mini Kit; Qiagen, Inc.). DNA was eluted in 100 μL of elution buffer. A sample containing only elution buffer was extracted with each batch of flea samples and used as a negative extraction control.
The absence of PCR inhibitors was determined by the amplification of a fragment of ∼200 bp of the mitochondrial 16S rDNA using primers Sen_mt16S: 5′-TACATAACACGAGAAGACC-3′ and Rev_mt16S: 5′-GTGATTGCGCTGTTATCC-3′, as previously described (Vobis et al. 2004), with modifications. Briefly, all assays were performed in 15 μL of total volume using 1 μL of extracted DNA as template. The reaction mixture contained 7.5 pmol of each primer and a 1× “ready to use” premix (SYBR Premix Ex Taq Tli RNAse H Plus, Takara Bio) containing an optimized buffer, nucleotides, a heat-activated DNA polymerase, and SYBR Green I for detection of amplifications using an intercalator-based real-time PCR technology. A 1× of ROX reference dye was used to standardize fluorescence detection across wells. All reactions were run on a StepOne Plus Real-Time PCR System (Life Technologies). Amplifications were performed under the following conditions: one hot-start cycle at 95°C for 1 min followed by 45 cycles of 94°C for 10 s, 56°C for 20 s, and 72°C for 20 s. Amplification was detected by SYBR Green I fluorescence with emission at 522 nm. Melting curves were obtained for each amplification using temperature increments of 0.5°C/s. Positive amplification were characterized by a melting temperature peak between 80.9°C and 82.5°C. Negative extraction controls were used in each batch tested. Flea DNA was amplified from all samples.
Rickettsia genus-specific primers (ompB-F1: 5′-CAGGATTGGTAACTGCTTC-3′; ompB-R1: 5′-ACCGTTACCTATATCACCGG-3′) were chosen to amplify a 1560 bp region spanning approximately the first one-third of the ompB gene, based on the R. felis Marseille-URRWXCal2 strain (GenBank ID NC_007109). PCRs were carried out in a 50 μL volume reaction mix (0.30 μM of each primer, 0.3 mM dNTPs, 1 mM MgCl2, 1× PCR buffer, and 0.05 U/μL Taq DNA polymerase) with 100 ng of DNA. Reactions were performed in a Veriti™ 96-well thermal cycler (Applied Biosystems) using the following conditions: 95°C for 2 min, 40 cycles of 95°C for 30 s, 54.8°C for 45 s, and 72°C for 30 s, followed by a final extension of 7 min at 72°C. Positive (Rickettsia rickettsii or Rickettsia amblyommii) and negative controls were used in each run. Amplicons were visualized on 1.5% agarose gels stained with ethidium bromide.
All positive PCR products were restriction enzyme digested with PstI to distinguish between R. felis (expected fragment sizes: 885, 597, 78 bp) and R. typhi (expected fragment sizes: 957, 300, 294 bp). Seventeen samples with the expected R. felis restriction pattern were chosen for nucleotide sequence analysis. In addition, restriction enzyme digestion of PCR products producing fragments distinctly different from R. felis or R. typhi was further characterized. Amplicons from PCR-positive samples were spin column purified (QIAquick PCR Purification Kit; Qiagen) and quantified photometrically (NanoVue Spectrophotometer; Fisher Scientific). The purified amplicons were sequenced bidirectionally by a commercial facility (Eurofins MWG Operon), using the same primer pairs designed for amplification plus two internal primers, ompB-F2 (5′-GCTGCGCCTTCTACATTA-3′) and ompB-R2 (5′-CAACAAATCCTGCAGCTC-3′). Sample sequences were compared to the R. felis Marseille-URRWXCal2 strain (GenBank ID NC_007109) and R. typhi Wilmington strain (GenBank ID NC_006142) reference sequences using CLC Main Workbench 7.0.3 (CLC Bio) and NCBI/BLAST (
Statistical analysis
The pooled infection rate of fleas with Rickettsia species was calculated using maximum likelihood estimation (MLE) methods for unequal pool sizes. The bias-corrected MLE method for point estimations and the skewness-corrected function to compute asymptotic 95% confidence limits were calculated using PooledInfRate version 4.0 (Brad Biggerstaff, CDC; available from
Results
All fleas in the study were identified as C. felis using a standard taxonomic key (Lewis et al. 1998). A total of 130 cats were sampled and therefore 130 flea pools collected: 76 from the flea-borne rickettsiosis endemic region of Los Angeles County [Animal Shelter 1 (n = 41) and Animal Shelter 2 (n = 35)] and 54 from the nonendemic region of Sacramento County (n = 29) and Contra Costa County (n = 25). A total of 520 fleas were collected, with a median of five fleas per cat. Rickettsia DNA was detected in 77 flea pools by PCR analysis (Table 1); the Rickettsia species in three samples from Southern California could not be determined due to insufficient DNA available for restriction enzyme digestion or sequencing. These samples were considered as positives, however, when analyzing the overall infection rate. The pooled infection rate by MLE of Rickettsia species was not significantly higher in fleas tested from the endemic region (42 pools out of 281 fleas tested; pooled infection rate of 213 per 1000 fleas) when compared to the nonendemic region (35 pools out of 239 fleas tested; pooled infection rate of 234 per 1000 fleas) (p = 0.283). When comparing the pooled infection rate of R. felis only, there was a significant difference between endemic and nonendemic areas: 29 pools out of 277 fleas tested (pooled infection rate of 131 per 1000 fleas) were R. felis positive from the endemic region and 35 pools out of 239 fleas tested (pooled infection rate of 234 per 1000 fleas) were R. felis positive from the nonendemic region (p < 0.0001).
TNR—trap, neuter, release program for stray and feral cats.
Coinfection of R. felis with Rickettsia Genotype 1.
Coinfection of R. felis and Rickettsia Genotype 2.
All amplicons but one from fleas collected in the nonendemic region had restriction patterns consistent with R. felis; none had restriction patterns consistent with R. typhi. Two restriction patterns, distinct from R. typhi and R. felis, were seen in amplicons from fleas collected in the endemic region and in a single flea pool from the nonendemic region. Genotype 1 (GenBank ID KP398499) was identified in one flea pool from Animal Shelter 1 in the endemic region. In addition, one flea pool from the nonendemic region (Sacramento County TNR program) presented with restriction patterns, suggesting a coinfection of Genotype 1 with R. felis. Genotype 2 (GenBank ID KP398500) was identified in 10 pools from the endemic region; 3 flea pools from Animal Shelter 1, 6 flea pools from Animal Shelter 2, and 1 flea pool from Animal Shelter 1, which had restriction patterns suggesting the presence of a coinfection of Genotype 2 with R. felis. Although the Rickettsia species could not be determined for three samples, the remaining 28 infected cat flea pools from the endemic region were shown to harbor R. felis DNA only (Table 1).
Sequence analysis of 17 amplicons (7 from Los Angeles County, 5 from Sacramento County, and 5 from Contra Costa County) having the expected restriction pattern for R. felis confirmed a 99–100% identity to the R. felis Marseille-URRWXCal2 strain (GenBank ID NC_007109). The nucleotide sequence of Genotype 1 was 99% identical to a proposed new species, Candidatus Rickettsia asemboensis (GenBank ID JN315972.2), and 94% identical to R. felis. Sequencing analysis of a flea pool from the nonendemic region demonstrated a coinfection of R. felis and Genotype 1. Genotype 2 was 99% identical to the newly reported Candidatus R. senegalensis identified in fleas from Senegal (GenBank ID KF666470.1) and shared 96% identity with R. felis. A phylogenetic tree (Fig. 1) was generated based on 1431-bp sequences of the ompB of the two R. felis-like genotypes from the flea pools and closely related organisms by using the maximum-likelihood method on the basis of the Kimura 2-parameter model (Kimura 1980), with bootstrap support based on 1000 replicates (Molecular Evolutionary Genetics Analysis [MEGA] software, version 6) (Tamura et al. 2013).

Phylogenetic tree based on 1431 bp from the outer membrane protein B (ompB) sequences of Rickettsia felis-like genotypes obtained from flea pools in Southern California (in boldface) and closely related organisms constructed by using the maximum-likelihood method on the basis of the Kimura 2-parameter model. Each bacterial name is followed by geographic origin data and the GenBank accession number is provided in parentheses. The numbers at the nodes indicate percentages of bootstrap support based on 1000 replicates. Percentages corresponding to partitions reproduced in fewer than 50% of bootstrap replicates are collapsed. The scale bar indicates 0.02 substitutions per nucleotide position.
Sequencing primers were designed to differentiate between R. felis and Genotype 2 when sequencing an amplicon suspected of containing both sequences, based on the restriction digest pattern. Primer ompB-15bp-R (5′-TACTGCATTAGCAGGACC-3′) binds to the wild-type R. felis sequence in a region that contains a 15 bp deletion in Genotype 2. Primer ompB-15bpdel-R (5′- CCATTATTAGTACCTATTACTGC-3′) binds to the Genotype 2 sequence but not to R. felis. The presence of a mixed infection of R. felis and Genotype 2 was confirmed in one flea pool collected from a cat at Los Angeles County Animal Shelter 1.
Discussion
With the exception of R. typhi, the global epidemiology and disease ecology of flea-borne Rickettsia species are still poorly understood. Such a conundrum exists in California where the majority of reported human cases (over 90%) occur in residents of Los Angeles and Orange counties, but rarely elsewhere in the state despite statewide prevalence of the implicated cat flea vector and its potential hosts. In the 1970s, a cycle involving R. typhi, cat fleas, and opossums was proposed to explain this Southern California phenomenon (Adams et al. 1970), and subsequent studies suggested that, rather than R. typhi, a different Rickettsia species named R. felis might be the etiologic agent of human illness in this region (Reif and Macaluso 2009, Eremeeva et al. 2012). The etiologic agent causing human flea-borne rickettsioses in California, however, is yet to be confirmed by molecular methods such as the detection of rickettsial DNA using PCR analysis or isolation of the infectious agent from clinical specimens. Rickettsiae are often found sequestered within endothelial cells reducing the number of organisms circulating in the bloodstream (Labruna and Walker 2014). Because of this, it can be extremely difficult to detect or isolate these organisms in a patient's blood sample even before antibiotic therapy has begun. In addition, the antibody-based testing protocols currently used for diagnosis use R. typhi antigens, which may result in cross-reactivity between antigenically related rickettsiae (Civen and Ngo 2008, Reif and Macaluso 2009) and thus prevent confirmation of the Rickettsia species involved.
In this study, the infection rate in cat fleas with R. felis and the absence of R. typhi were consistent with findings from past studies in Southern California, which demonstrated a minimum infection prevalence of 13–32% R. felis and 0–2% R. typhi (Abramowicz et al. 2012, Eremeeva et al. 2012). R. felis has also been reported from regions of Texas and Mexico endemic for flea-borne rickettsioses. In Corpus Christi, Texas, 14 (2.6%) of 529 cat fleas recovered from opossums were PCR positive for the presence of Rickettsia DNA and sequencing analysis demonstrated R. felis and R. typhi presence in 11 and 3 cat flea specimens, respectively (Boostrom et al. 2002). In Yucatán, Mexico, 11 (20%) of 54 pools of cat fleas collected from domestic dogs were positive for Rickettsia DNA and sequences from 2 positive flea pools were 100% identical to R. felis (Zavala-Velázquez et al. 2002). The wide distribution of R. felis in cat fleas from flea-borne rickettsioses endemic and nonendemic regions of North America suggests that a high prevalence of this bacterium in cat fleas does not predispose to increase risk of human infection. Human flea-borne rickettsioses cases are uncommon in California with a total of 217 confirmed cases reported between 2010 and 2014 (
A comparison of cat flea pools positive for R. felis among the three California counties sampled in this study revealed that the infection rate of this rickettsia species in regions considered flea-borne rickettsioses nonendemic (i.e., Sacramento and Contra Costa counties) was higher than in Los Angeles County, where ∼70% of the human cases are reported. Thus, differential risk of acquiring flea-borne rickettsioses between endemic and nonendemic areas of California cannot be ascribed to differential exposure to R. felis from infected cat fleas. Other factors may be responsible for the nonendemic status of regions in California such as lack of identification of the disease agent by physicians when patients present with potentially nonspecific symptoms of flea-borne rickettsioses (Civen and Ngo 2008) or involvement of other hematophagous arthropods in R. felis transmission. R. felis DNA has been reported previously in other flea species, ticks, and mosquitoes, suggesting the role of other ectoparasites as potential vectors of R. felis (Socolovschi et al. 2012, Abarca et al. 2013, Ramírez-Hernández et al. 2013, Dieme et al. 2015, Peniche-Lara et al. 2015).
The identification of two R. felis-like genotypes in cat fleas raises questions as to whether other Rickettsia species may be causative agents of flea-borne rickettsioses. Genotype 1, detected in two flea pools, is closely related to Candidatus R. asemboensis, a newly recognized rickettsia detected in five species of fleas (Ctenocephalides canis, C. felis, Echidnophaga gallinacea, X. cheopis, and Pulex irritans) from Asembo, Kenya (Jiang et al. 2013). In that study, four C. canis and 71 (94.7%) of 75 C. felis pools were infected with Candidatus R. asemboensis, whereas R. felis was detected in only three C. felis pools. One cat flea pool was coinfected with both agents (Jiang et al. 2013). Genotype 2, detected in 10 flea pools in this study, was most similar to Candidatus R. senegalensis recently detected in cat fleas from Senegal (Mediannikov et al. 2014). The distribution of these rickettsiae in cat flea populations and their potential to cause human disease should be considered for further investigation.
Why human cases of flea-borne rickettsioses cases in California are primarily reported from Los Angeles and Orange counties remains unclear. The results from this study demonstrate that R. felis is present in cat fleas in both endemic and nonendemic regions of the state and suggest that at least two additional R. felis-like genotypes may be circulating in populations of these fleas. However, because only a limited number of fleas and locations were examined in this study, a more comprehensive survey might shed more light on the prevalence of these Rickettsia species and R. typhi. Other considerations in future research should include flea pool size and PCR targets. Comparisons across uneven numbers of pooled fleas, combined with a PCR targeting a relatively large DNA fragment, have the potential to produce a dilution effect that may decrease the efficiency of rickettsia detection. At the same time, the current study design required amplification of a large DNA fragment, when used in conjunction with the restriction enzyme digestion, to adequately differentiate between R. felis and R. typhi. Many questions regarding the basic biology and ecology of flea-borne rickettsioses in California are yet to be answered, including confirmation of the etiologic agent(s).
Footnotes
Acknowledgments
The authors thank Crystal Perreira and Leslie Foss for their assistance with flea collections from the Contra Costa Co. TNR program, Jessica Chinison and Dr. Dominique Griffon for lab support, and Dr. Mark Novak for his assistance with flea collections from the Sacramento Co. TNR program. They thank the veterinarians and staff from each collection site for providing the opportunity to collect flea specimens. The authors acknowledge Dr. Vicki Kramer and Dr. Charsey Cole Porse for providing critical review of the article.
Author Disclosure Statement
No competing financial interests exist.
