A chromosome-level assembly of the cat flea genome uncovers rampant gene duplication and genome size plasticity

BACKGROUND

Fleas (Insecta: Siphonaptera) are small flightless parasites of birds and mammals; their blood-feeding can transmit many serious pathogens (i.e., the etiological agents of bubonic plague, endemic and murine typhus). The lack of flea genome assemblies has hindered research, especially comparisons to other disease vectors. Accordingly, we sequenced the genome of the cat flea, Ctenocephalides felis, an insect with substantial human health and veterinary importance across the globe.

RESULTS

By combining Illumina and PacBio sequencing of DNA derived from multiple inbred female fleas with Hi-C scaffolding techniques, we generated a chromosome-level genome assembly for C. felis. Unexpectedly, our assembly revealed extensive gene duplication across the entire genome, exemplified by ~ 38% of protein-coding genes with two or more copies and over 4000 tRNA genes. A broad range of genome size determinations (433-551 Mb) for individual fleas sampled across different populations supports the widespread presence of fluctuating copy number variation (CNV) in C. felis. Similarly, broad genome sizes were also calculated for individuals of Xenopsylla cheopis (Oriental rat flea), indicating that this remarkable "genome-in-flux" phenomenon could be a siphonapteran-wide trait. Finally, from the C. felis sequence reads, we also generated closed genomes for two novel strains of Wolbachia, one parasitic and one symbiotic, found to co-infect individual fleas.

CONCLUSION

Rampant CNV in C. felis has dire implications for gene-targeting pest control measures and stands to complicate standard normalization procedures utilized in comparative transcriptomics analysis. Coupled with co-infection by novel Wolbachia endosymbionts-potential tools for blocking pathogen transmission-these oddities highlight a unique and underappreciated disease vector.

The Rickettsia type IV secretion system: unrealized complexity mired by gene family expansion

Many prokaryotes utilize type IV secretion systems (T4SSs) to translocate substrates (e.g. nucleoprotein, DNA, protein) across the cell envelope, and/or to elaborate surface structures (i.e. pili or adhesins). Among eight distinct T4SS classes, P-T4SSs are typified by the Agrobacterium tumefaciens vir T4SS, which is comprised of 12 scaffold components (VirB1-VirB11, VirD4). While most P-T4SSs include all 12 Vir proteins, some differ from the vir archetype by either containing additional scaffold components not analogous to Vir proteins or lacking one or more of the Vir proteins. In a special case, the Rickettsiales vir homolog (rvh) P-T4SS comprises unprecedented gene family expansion. rvh contains three families of gene duplications (rvhB9, rvhB8, rvhB4): RvhB9,8,4-I are conserved relative to equivalents in other P-T4SSs, while RvhB9,8,4-II have evolved atypical features that deviate substantially from other homologs. Furthermore, rvh contains five VirB6-like genes (rvhB6a-e), which are tandemly arrayed and contain large N- and C-terminal extensions. Our work herein focuses on the complexity underpinned by rvh gene family expansion. Furthermore, we describe an RvhB10 insertion, which occurs in a region that forms the T4SS pore. The significance of these curious properties to rvh structure and function is evaluated, shedding light on a highly complex T4SS.

Which Way In? The RalF Arf-GEF Orchestrates Rickettsia Host Cell Invasion

Bacterial Sec7-domain-containing proteins (RalF) are known only from species of Legionella and Rickettsia, which have facultative and obligate intracellular lifestyles, respectively. L. pneumophila RalF, a type IV secretion system (T4SS) effector, is a guanine nucleotide exchange factor (GEF) of ADP-ribosylation factors (Arfs), activating and recruiting host Arf1 to the Legionella-containing vacuole. In contrast, previous in vitro studies showed R. prowazekii (Typhus Group) RalF is a functional Arf-GEF that localizes to the host plasma membrane and interacts with the actin cytoskeleton via a unique C-terminal domain. As RalF is differentially encoded across Rickettsia species (e.g., pseudogenized in all Spotted Fever Group species), it may function in lineage-specific biology and pathogenicity. Herein, we demonstrate RalF of R. typhi (Typhus Group) interacts with the Rickettsia T4SS coupling protein (RvhD4) via its proximal C-terminal sequence. RalF is expressed early during infection, with its inactivation via antibody blocking significantly reducing R. typhi host cell invasion. For R. typhi and R. felis (Transitional Group), RalF ectopic expression revealed subcellular localization with the host plasma membrane and actin cytoskeleton. Remarkably, R. bellii (Ancestral Group) RalF showed perinuclear localization reminiscent of ectopically expressed Legionella RalF, for which it shares several structural features. For R. typhi, RalF co-localization with Arf6 and PI(4,5)P₂ at entry foci on the host plasma membrane was determined to be critical for invasion. Thus, we propose recruitment of PI(4,5)P₂ at entry foci, mediated by RalF activation of Arf6, initiates actin remodeling and ultimately facilitates bacterial invasion. Collectively, our characterization of RalF as an invasin suggests that, despite carrying a similar Arf-GEF unknown from other bacteria, different intracellular lifestyles across Rickettsia and Legionella species have driven divergent roles for RalF during infection. Furthermore, our identification of lineage-specific Arf-GEF utilization across some rickettsial species illustrates different pathogenicity factors that define diverse agents of rickettsial diseases.

Phylogenomics analysis indicates divergent mechanisms for host cell invasion across diverse species of obligate intracellular Rickettsia. For instance, only some Rickettsia species carry RalF, the rare bacterial Arf-GEF effector utilized by Legionella pneumophila to facilitate fusion of ER-derived membranes with its host-derived vacuole. For R. prowazekii (Typhus Group, TG), prior in vitro studies suggested the Arf-GEF activity of RalF, which is absent from Spotted Fever Group species, might be spatially regulated at the host plasma membrane. Herein, we demonstrate RalF of R. typhi (TG) and R. felis (Transitional Group) localizes to the host plasma membrane, yet R. bellii (Ancestral Group) RalF shows perinuclear localization reminiscent of RalF-mediated recruitment of Arf1 by L. pneumophila to its vacuole. For R. typhi, RalF expression occurs early during infection, with RalF inactivation significantly reducing host cell invasion. Furthermore, RalF co-localization with Arf6 and the phosphoinositide PI(4,5)P₂ at the host plasma membrane was determined to be critical for R. typhi invasion. Thus, our work illustrates that different intracellular lifestyles across species of Rickettsia and Legionella have driven divergent roles for RalF during host cell infection. Collectively, we identify lineage-specific Arf-GEF utilization across diverse rickettsial species, previously unappreciated mechanisms for host cell invasion and infection.

Novel identification of Dermacentor variabilis Arp2/3 complex and its role in rickettsial infection of the arthropod vector

Tick-borne spotted fever group (SFG) Rickettsia species must be able to infect both vertebrate and arthropod host cells. The host actin-related protein 2/3 (Arp2/3) complex is important in the invasion process and actin-based motility for several intracellular bacteria, including SFG Rickettsia in Drosophila and mammalian cells. To investigate the role of the tick Arp2/3 complex in tick-Rickettsia interactions, open reading frames of all subunits of the protein including Arp2, Arp3, ARPC1, ARPC2, ARPC3, ARPC4, and ARPC5 were identified from Dermacentor variabilis. Amino acid sequence analysis showed variation (ranging from 25–88%) in percent identity compared to the corresponding subunits of the complex from Drosophila melanogaster, Mus musculus, Homo sapiens, and Saccharomyces cerevisiae. Potential ATP binding sites were identified in D. variabilis (Dv) Arp2 and Arp3 subunits as well as five putative WD (Trp-Asp) motifs which were observed in DvARPC1. Transcriptional profiles of all subunits of the DvArp2/3 complex revealed greater mRNA expression in both Rickettsia-infected and -uninfected ovary compared to midgut and salivary glands. In response to R. montanensis infection of the tick ovary, the mRNA level of only DvARPC4 was significantly upregulated compared to uninfected tissues. Arp2/3 complex inhibition bioassays resulted in a decrease in the ability of R. montanensis to invade tick tissues with a significant difference in the tick ovary, indicating a role for the Arp2/3 complex in rickettsial invasion of tick cells. Characterization of tick-derived molecules associated with rickettsial infection is imperative in order to better comprehend the ecology of tick-borne rickettsial diseases.

Molecular and functional characterization of vacuolar-ATPase from the American dog tick Dermacentor variabilis

Vacuolar (V)-ATPase is a proton-translocating enzyme that acidifies cellular compartments for various functions such as receptor-mediated endocytosis, intracellular trafficking and protein degradation. Previous studies in Dermacentor variabilis chronically infected with Rickettsia montanensis have identified V-ATPase as one of the tick-derived molecules transcribed in response to rickettsial infection. To examine the role of the tick V-ATPase in tick-Rickettsia interactions, a full-length 2887-bp cDNA (2532-bp open reading frame) clone corresponding to the transcript of the V0 domain subunit a of D. variabilis V-ATPase (DvVATPaseV0a) gene encoding an 843 amino acid protein with an estimated molecular weight of ~96 kDa was isolated from D. variabilis. Amino acid sequence analysis of DvVATPaseV0a showed the highest similarity to VATPaseV0a from Ixodes scapularis. A potential N-glycosylation site and eight putative transmembrane segments were identified in the sequence. Western blot analysis of tick tissues probed with polyclonal antibody raised against recombinant DvVATPaseV0a revealed the expression of V-ATPase in the tick ovary. Transcriptional profiles of DvVATPaseV0a demonstrated a greater mRNA expression in the tick ovary, compared with the midgut and salivary glands; however, the mRNA level in each of these tick tissues remained unchanged after infection with R. montanensis for 1 h. V-ATPase inhibition bioassays resulted in a significant decrease in the ability of R. montanensis to invade tick cells in vitro, suggesting a role of V-ATPase in rickettsial infection of tick cells. Characterization of tick-derived molecules involved in rickettsial infection is essential for a thorough understanding of rickettsial transmission within tick populations and the ecology of tick-borne rickettsial diseases.