Abstract Background The stable fly, Stomoxys calcitrans , is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the USA alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies. Results This study examines stable fly biology by utilizing a combination of high-quality genome sequencing and RNA-Seq analyses targeting multiple developmental stages and tissues. In conjunction, 1600 genes were manually curated to characterize genetic features related to stable fly reproduction, vector host interactions, host-microbe dynamics, and putative targets for control. Most notable was characterization of genes associated with reproduction and identification of expanded gene families with functional associations to vision, chemosensation, immunity, and metabolic detoxification pathways. Conclusions The combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-Seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and new data will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship of Stomoxys to other blood-feeding (horn flies and Glossina ) and non-blood-feeding flies (house flies, medflies, Drosophila ) will facilitate understanding of the evolutionary processes associated with development of blood feeding among the Cyclorrhapha.
Abstract Background The stable fly, Stomoxys calcitrans , is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the United States alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies. Results This study examines stable fly biology by utilizing a combination of high-quality genome sequencing, microbiome analyses, and RNA-seq analyses targeting multiple developmental stages and tissues. In conjunction, manual curation of over 1600 genes was used to examine gene content related to stable fly reproduction, interactions with their host, host-microbe dynamics, and putative routes for control. Most notable was establishment of reproduction-associated genes and identification of expanded vision, chemosensation, immune repertoire, and metabolic detoxification pathway gene families. Conclusions The combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and knowledge will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship of Stomoxys to other blood-feeding ( Glossina ) and non-blood-feeding flies (medflies, Drosophila , house flies) will allow for understanding the evolution of blood feeding among Cyclorrhapha flies.
Acaricide resistant Rhipicephalus microplus populations have become a major problem for many cattle producing areas of the world. Pyrethroid resistance in arthropods is typically associated with mutations in domains I, II, III, and IV of voltage-gated sodium channel genes. In R. microplus, known resistance mutations include a domain II change (C190A) in populations from Australia, Africa, and South America and a domain III mutation (T2134A) that only occurs in Mexico and the U.S. We investigated pyrethroid resistance in cattle fever ticks from Texas and Mexico by estimating resistance levels in field-collected ticks using larval packet discriminating dose (DD) assays and identifying single nucleotide polymorphisms (SNPs) in the para-sodium channel gene that associated with resistance. We then developed qPCR assays for three SNPs and screened a larger set of 1,488 R. microplus ticks, representing 77 field collections and four laboratory strains, for SNP frequency. We detected resistance SNPs in 21 of 68 U.S. field collections and six of nine Mexico field collections. We expected to identify the domain III SNP (T2134A) at a high frequency; however, we only found it in three U.S. collections. A much more common SNP in the U.S. (detected in 19 of 21 field collections) was the C190A domain II mutation, which has never before been reported from North America. We also discovered a novel domain II SNP (T170C) in ten U.S. and two Mexico field collections. The T170C transition mutation has previously been associated with extreme levels of resistance (super-knockdown resistance) in insects. We found a significant correlation (r = 0.81) between the proportion of individuals in field collections that carried any two resistance SNPs and the percent survivorship of F1 larvae from these collections in DD assays. This relationship is accurately predicted by a simple linear regression model (R2 = 0.6635). These findings demonstrate that multiple mutations in the para-sodium channel gene independently associate with pyrethroid resistance in R. microplus ticks, which is likely a consequence of human-induced selection.
Acaricide resistant Rhipicephalus microplus populations have become a major problem for many cattle producing areas of the world. Pyrethroid resistance in arthropods is typically associated with mutations in domains I, II, III, and IV of voltage-gated sodium channel genes. In R. microplus, known resistance mutations include a domain II change (C190A) in populations from Australia, Africa, and South America and a domain III mutation (T2134A) that only occurs in Mexico and the U.S. We investigated pyrethroid resistance in cattle fever ticks from Texas and Mexico by estimating resistance levels in field-collected ticks using larval packet discriminating dose (DD) assays and identifying single nucleotide polymorphisms (SNPs) in the para-sodium channel gene that associated with resistance. We then developed qPCR assays for three SNPs and screened a larger set of 1,488 R. microplus ticks, representing 77 field collections and four laboratory strains, for SNP frequency. We detected resistance SNPs in 21 of 68 U.S. field collections and six of nine Mexico field collections. We expected to identify the domain III SNP (T2134A) at a high frequency; however, we only found it in three U.S. collections. A much more common SNP in the U.S. (detected in 19 of 21 field collections) was the C190A domain II mutation, which has never before been reported from North America. We also discovered a novel domain II SNP (T170C) in ten U.S. and two Mexico field collections. The T170C transition mutation has previously been associated with extreme levels of resistance (super-knockdown resistance) in insects. We found a significant correlation (r = 0.81) between the proportion of individuals in field collections that carried any two resistance SNPs and the percent survivorship of F1 larvae from these collections in DD assays. This relationship is accurately predicted by a simple linear regression model (R2 = 0.6635). These findings demonstrate that multiple mutations in the para-sodium channel gene independently associate with pyrethroid resistance in R. microplus ticks, which is likely a consequence of human-induced selection.
Abstract Background Rhipicephalus (Boophilus) microplus causes significant cattle production losses worldwide because it transmits Babesia bovis and B. bigemina causative agents of bovine babesiosis. Control of these ticks primarily has relied on treatment of cattle with chemical acaricides, but frequent use, exacerbated by the one-host life cycle of these ticks, has led to high-level resistance to multiple classes of acaricides. Consequently, new approaches for control, such as anti-tick vaccines, are critically important. Key to this approach is targeting highly conserved antigenic epitopes to reduce the risk of vaccine escape in heterologous tick populations. Methods We evaluated amino acid conservation within 14 tick proteins across 167 R. microplus collected from geographically diverse locations in the Americas and Pakistan using PCR amplicon sequencing and in silico translation of exons. Results We found that amino acid conservation varied considerably across these proteins. Only one target, the voltage-dependent anion channel, was fully conserved in all 167 R. microplus samples (protein similarity 1.0). Five other proteins were highly conserved: the aquaporin RmAQP1 (0.989), vitellogenin receptor (0.985), serpin-1 (0.985), and subolesin (0.981). In contrast, the glycoprotein protease Bm86 was one of the least conserved (0.889). The Bm86 sequence used in the original Australian TickGARD vaccine carried many amino acid replacements compared to the R. microplus populations examined here, supporting the hypothesis that this vaccine target is not optimal for use in the Americas. By mapping amino acid replacements onto predicted 3D protein models, we also identified amino acid changes within several small peptide vaccines targeting portions of the aquaporin RmAQP2, chitinase, and Bm86. Conclusions These findings emphasize the importance of thoroughly analyzing protein variation within anti-tick vaccine targets across diverse tick populations before selecting candidate vaccine antigens. When considering protein conservation alone, RmAQP1, vitellogenin receptor, serpin-1, subolesin, and especially the voltage-dependent anion channel rank as high priority anti-tick vaccine candidates for use in the Americas and perhaps globally.
