Pyrosequencing and characterization of immune response genes from the American dog tick, Dermacentor variabilis (L.)

Saliva proteome of partially- and fully-engorged adult female Haemaphysalis flava ticks

Tick saliva is a reservoir of bioactive proteins. Saliva protein compositions change dynamically during blood-feeding. Decipherment of protein profiles in different blood-feeding stages may bring deeper insight into tick feeding physiology and provide targets for immunologic control alternatives. However, having the infancy of tick genome sequencing, assembly, annotation, and limited knowledge of tick salivary proteins restrain the data interpretation. Here, we aimed to depict the saliva protein profile in partially- (PE) and fully-engorged (FE) Haemaphysalis flava ticks, with a special focus on the analysis of those uncharacterized proteins. Saliva was collected from PE and FE adult female H. flava ticks. Saliva proteins were analyzed by high-performance liquid chromatography-tandem mass spectrometry (HPLC/MS-MS). MS data were searched against an in-house salivary gland transcriptome library for identification of tick-derived proteins. Abundances of proteins were compared between PE and FE ticks. The uncharacterized proteins detected in saliva were further bioinformatically analyzed. In total, 614 proteins were identified including 94 host proteins and 520 tick-derived proteins. The 226 tick-derived high-confidence proteins were classified into 10 categories: transporters, enzymes, protease inhibitors, immunity-related proteins, lipocalins, glycine-rich proteins, muscle proteins, secreted proteins, uncharacterized proteins and others. A total of 98 proteins were shared in both PE and FE with 74 only in PE and 54 only in FE. Abundances of 24 shared proteins were significantly higher in PE. The profile of top 15 most abundant proteins was also different between PE and FE ticks. The 65 uncharacterized proteins detected in tick saliva were branched into subclusters 1 A, 1B, 2, 3 A, 3B and 3 C based on particular motifs like RGD, LRR, indicating their diverse predicted functions like anti-coagulation, regulation of innate immune, or other functions. This study provides and compares saliva proteomes of H. flava ticks in two feeding stages with special cluster analysis on the uncharacterized proteins. Further investigations are needed to confirm the roles of these uncharacterized proteins in ticks.

Early Transcriptional Changes in the Midgut of Ornithodoros moubata after Feeding and Infection with Borrelia duttonii

Studies on tick-pathogen-host interactions are helping to identify candidates for vaccines against ticks and tick-borne diseases and to discover potent bioactive tick molecules. The tick midgut is the main tissue involved in blood feeding and, moreover, the first organ to have contact with pathogens ingested through the blood meal. As little is known about the molecular biology of feeding and tick defence mechanisms against microorganisms, but important for understanding vector-pathogen interactions, we explored the early transcriptional changes in the midgut of Ornithodoros moubata after feeding and in response to challenge with the relapsing-fever spirochete Borrelia duttonii using the Ion S5XL platform. Besides transcripts with metabolic function and immune-related transcripts we discovered numerous putative and uncharacterized protein sequences. Overall, our analyses support previous studies and provides a valuable reference database for further functional proteomic analysis of midgut proteins of O. moubata.

Cracking the Challenge of Antimicrobial Drug Resistance with CRISPR/Cas9, Nanotechnology and Other Strategies in ESKAPE Pathogens

Antimicrobial resistance is mushrooming as a silent pandemic. It is considered among the most common priority areas identified by both national and international agencies. The global development of multidrug-resistant strains now threatens public health care improvement by introducing antibiotics against infectious agents. These strains are the product of both continuous evolution and unchecked antimicrobial usage (AMU). The ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) are the leading cause of nosocomial infections throughout the world. Most of them are now multidrug-resistant, which pose significant challenges in clinical practice. Understanding these bacteria’s resistance mechanisms is crucial for developing novel antimicrobial agents or other alternative tools to fight against these pathogens. A mechanistic understanding of resistance in these pathogens would also help predict underlying or even unknown mechanisms of resistance of other emerging multidrug-resistant pathogens. Research and development to find better antibacterial drugs and research on tools like CRISPER-Cas9, vaccines, and nanoparticles for treatment of infections that can be further explored in the clinical practice health sector have recognized these alternatives as essential and highly effective tools to mitigate antimicrobial resistance. This review summarizes the known antimicrobial resistance mechanisms of ESKAPE pathogens and strategies for overcoming this resistance with an extensive overview of efforts made in this research area.

Comparative Proteomic Analysis of Rhipicephalus sanguineus sensu lato (Acari: Ixodidae) Tropical and Temperate Lineages: Uncovering Differences During Ehrlichia canis Infection

The tick vector Rhipicephalus sanguineus is established as a complex of closely related species with high veterinary-medical significance, in which the presence of different genetic, morphological, and biological traits has resulted in the recognition of different lineages within taxa. One of the most striking differences in the "temperate" and "tropical" lineages of R. sanguineus (s.l.) is the vector competence to Ehrlichia canis, suggesting that these ticks tolerate and react differently to pathogen infection. The present study addresses the SG and MG proteome of the R. sanguineus tropical and temperate lineages and compares their proteomic profile during E. canis infection. Batches of nymphs from the two lineages were allowed to feed on naïve and experimentally E. canis infected dogs and after molting, adults were dissected, and salivary glands and midgut tissues separated. Samples were screened for the presence of E. canis before proteomic analyses. The representation of the proteins identified in infected and non-infected tissues of each lineage was compared and gene ontology used for protein classification. Results highlight important differences in those proteomic profiles that added to previous reported genetic, biological, behavioral, and morphological differences, strengthening the hypothesis of the existence of two different species. Comparing infected and non-infected tissues, the results show that, while in midgut tissues the response to E. canis infection is similar in the salivary glands, the two lineages show a different pattern of protein representation. Focusing on the proteins found only in the infected condition, the data suggests that the cement cone produced during tick feeding may be implicated in pathogen infection. This study adds useful information to the debate on the controversial R. sanguineus systematic status, to the discussion related with the different vectorial competence occurring between the two lineages and identifies potential targets for efficient tick and tick-borne disease control.

Transcriptome analysis of the innate immune system of Hyalomma asiaticum

Ticks are considered to be the second most important vectors of human infectious diseases. The innate immune system is the key factor that affects its vector competence. Hyalomma asiaticum is the primary vector of Crimean-Congo hemorrhagic fever virus (CCHFV). However, the immune system of H. asiaticum remains virtually unknown. Here, a high throughput full-length mRNA sequencing method was adopted to define the immunotranscriptome of H. asiaticum infected with the fungal pathogen Beauveria bassiana and gram-negative bacterium Enterobacter cloacae. The analysis yielded 22,300 isoforms with an average length of 3233 bps. In total, 68 potential immunity-related genes were identified based on similarity to the homologs known to be involved in immunity. These included most members of the Toll and JAK/STAT signaling pathways, but not the IMD signaling pathway. Moreover, two copies of Dicer-2 and five copies of Argonaute-2 were detected. These genes are postulated to be involved in the RNA interference (RNAi) pathway, which is an important defense against RNA viruses. Overall, this study provides the foundation for understanding the immune response of H. asiaticum to CCHFV.

Transcriptome profiling analysis of the intoxication response in midgut tissue of Agrotis ipsilon larvae to Bacillus thuringiensis Vip3Aa protoxin

Vip insecticidal proteins are produced by Bacillus thuringiensis (Bt) during its vegetative growth phase. In the present study, Vip3Aa11 and Vip3Aa39 proteins were investigated. These two proteins present 39 amino acid differential sites and they shared 95.06% amino acid sequence similarity. They are effective against some Lepidoptera insect larvae. In a previous study, using artificial diet bioassays, we estimated the LC₅₀ of Vip3Aa11 and Vip3Aa39 strains against Agrotis ipsilon larvae were 73.41 μg/mL (with 95% confidence interval of 2.34-11.19) and 5.43 μg/mL (with 95% confidence interval of 43.20-115.03), respectively. To investigate the response of Agrotis ipsilon transcriptome in defending against Vip3Aa11 and Vip3Aa39 toxins, we performed high-throughput RNA-sequencing on cDNA generated from the midguts of Agrotis ipsilon larvae that consumed a control diet (CK-M-A), Vip3Aa11 (Vip3Aa11-M-A) and Vip3Aa39 (Vip3Aa39-M-A) proteins. We generated about 98.87 Gb bases in total on BGISEQ-500 sequencing platform. After assembling all samples together and filtering the abundance, we got 51,887 unigenes, the total length, average length, N50 and GC content of unigenes are 64,523,651 bp, 1243 bp, 2330 bp and 41.81% respectively. We revealed 558 midgut genes differential expressed in Vip3Aa11-M-A and 65 midgut genes differentially expressed in Vip3Aa39-M-A. The differentially expressed genes were enriched for serine proteases and potential Bt Vip toxin midgut receptor genes. Eleven serine proteases related genes and 13 Bt toxin potential receptor genes with differential expression were found. Based on transcriptome profiling, we focused on validation the sensitivity of these two Vip3Aa proteins to trypsin and their binding properties to Agrotis ipsilon midgut BBMV (Brush Border Membrane Vesicles). The results show that the sensitivity of the two proteins to trypsin is similar. Binding experiments revealed that both proteins can bind to Agrotis ipsilon midgut BBMV, and there is a competitive binding between them. This transcriptome dataset provided a comprehensive sequence resource of Agrotis ipsilon and provides a foundation for comparative studies with other species of insects.

Comparative proteomics of the vector Dermacentor reticulatus revealed differentially regulated proteins associated with pathogen transmission in response to laboratory infection with Rickettsia slovaca

BACKGROUND

Tick-borne rickettsial diseases are caused by pathogens acquired from hard ticks. In particular, Rickettsia slovaca, a zoonotic infectious bacterium causing tick-borne lymphadenopathy (TIBOLA), is transmitted by the vectors Dermacentor spp. that can be found all over Europe. Although recent studies point out the extreme complexity of bacteria-induced effects in these blood-feeding vectors, the knowledge of individual molecules involved in the preservation and transmission of the pathogen is still limited. System biology tools, including proteomics, may contribute greatly to the understanding of pathogen-tick-host interactions.

METHODS

Herein, we performed a comparative proteomics study of the tick vector Dermacentor reticulatus that was experimentally infected with the endosymbiotic bacterium R. slovaca. Rickettsia-free ticks, collected in the southern region of Slovakia, were infected with the bacterium by a capillary tube-feeding system, and the dynamics of infection was assessed by quantitative PCR method after 5, 10, 15 and 27 days.

RESULTS

At the stage of controlled proliferation (at 27 dpi), 33 (from 481 profiled) differentially abundant protein spots were detected on a two-dimensional gel. From the aforementioned protein spots, 21 were successfully identified by tandem mass spectrometry.

CONCLUSIONS

Although a few discovered proteins were described as having structural or housekeeping functions, the vast majority of the affected proteins were suggested to be essential for tick attachment and feeding on the host, host immune system evasion and defensive response modulation to ensure successful pathogen transmission.

Transcriptome profiling of Diachasmimorpha longicaudata towards useful molecular tools for population management

Background

Diachasmimorpha longicaudata (Hymenoptera: Braconidae) is a solitary parasitoid of Tephritidae (Diptera) fruit flies of economic importance currently being mass-reared in bio-factories and successfully used worldwide. A peculiar biological aspect of Hymenoptera is its haplo-diploid life cycle, where females (diploid) develop from fertilized eggs and males (haploid) from unfertilized eggs. Diploid males were described in many species and recently evidenced in D. longicaudata by mean of inbreeding studies. Sex determination in this parasitoid is based on the Complementary Sex Determination (CSD) system, with alleles from at least one locus involved in early steps of this pathway. Since limited information is available about genetics of this parasitoid species, a deeper analysis on D. longicaudata’s genomics is required to provide molecular tools for achieving a more cost effective production under artificial rearing conditions.

Results

We report here the first transcriptome analysis of male-larvae, adult females and adult males of D. longicaudata using 454-pyrosequencing. A total of 469766 reads were analyzed and 8483 high-quality isotigs were assembled. After functional annotation, a total of 51686 unigenes were produced, from which, 7021 isotigs and 20227 singletons had at least one BLAST hit against the NCBI non-redundant protein database. A preliminary comparison of adult female and male evidenced that 98 transcripts showed differential expression profiles, with at least a 10-fold difference. Among the functionally annotated transcripts we detected four sequences potentially involved in sex determination and three homologues to two known genes involved in the sex determination cascade. Finally, a total of 4674SimpleSequence Repeats (SSRs) were in silico identified and characterized.

Conclusion

The information obtained here will significantly contribute to the development of D. longicaudata functional genomics, genetics and population-based genome studies. Thousands of new microsatellite markers were identified as toolkits for population genetics analysis. The transcriptome characterized here is the starting point to elucidate the molecular bases of the sex determination mechanism in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2759-2) contains supplementary material, which is available to authorized users.

Pyrosequencing the Midgut Transcriptome of the Banana Weevil Cosmopolites sordidus (Germar) (Coleoptera: Curculionidae) Reveals Multiple Protease-Like Transcripts

The banana weevil Cosmopolites sordidus is an important and serious insect pest in most banana and plantain-growing areas of the world. In spite of the economic importance of this insect pest very little genomic and transcriptomic information exists for this species. In the present study, we characterized the midgut transcriptome of C. sordidus using massive 454-pyrosequencing. We generated over 590,000 sequencing reads that assembled into 30,840 contigs with more than 400 bp, representing a significant expansion of existing sequences available for this insect pest. Among them, 16,427 contigs contained one or more GO terms. In addition, 15,263 contigs were assigned an EC number. In-depth transcriptome analysis identified genes potentially involved in insecticide resistance, peritrophic membrane biosynthesis, immunity-related function and defense against pathogens, and Bacillus thuringiensis toxins binding proteins as well as multiple enzymes involved with protein digestion. This transcriptome will provide a valuable resource for understanding larval physiology and for identifying novel target sites and management approaches for this important insect pest.

Sialomes and Mialomes: A Systems-Biology View of Tick Tissues and Tick-Host Interactions

Tick saliva facilitates tick feeding and infection of the host. Gene expression analysis of tick salivary glands and other tissues involved in host-pathogen interactions has revealed a wide range of bioactive tick proteins. Transcriptomic analysis has been a milestone in the field and has recently been enhanced by next-generation sequencing (NGS). Furthermore, the application of quantitative proteomics to ticks with unknown genomes has provided deeper insights into the molecular mechanisms underlying tick hematophagy, pathogen transmission, and tick-host-pathogen interactions. We review current knowledge on the transcriptomics and proteomics of tick tissues from a systems-biology perspective and discuss future challenges in the field.

Prediction of Leymus arenarius (L.) antimicrobial peptides based on de novo transcriptome assembly

Leymus arenarius is a unique wild growing Poaceae plant exhibiting extreme tolerance to environmental conditions. In this study we for the first time performed whole-transcriptome sequencing of lymegrass seedlings using Illumina platform followed by de novo transcriptome assembly and functional annotation. Our goal was to identify transcripts encoding antimicrobial peptides (AMPs), one of the key components of plant innate immunity. Using the custom software developed for this study that predicted AMPs and classified them into families, we revealed more than 160 putative AMPs in lymegrass seedlings. We classified them into 7 families based on their cysteine motifs and sequence similarity. The families included defensins, thionins, hevein-like peptides, snakins, cyclotide, alfa-hairpinins and LTPs. This is the first communication about the presence of almost all known AMP families in trascriptomic data of a single plant species. Additionally, cysteine-rich peptides that potentially represent novel families of AMPs were revealed. We have confirmed by RT-PCR validation the presence of 30 transcripts encoding selected AMPs in lymegrass seedlings. In summary, the presented method of pAMP prediction developed by us can be applied for relatively fast and simple screening of novel components of plant immunity system and is well suited for whole-transcriptome or genome analysis of uncharacterized plants.

Deep Sequencing Analysis of the Ixodes ricinus Haemocytome

BACKGROUND

Ixodes ricinus is the main tick vector of the microbes that cause Lyme disease and tick-borne encephalitis in Europe. Pathogens transmitted by ticks have to overcome innate immunity barriers present in tick tissues, including midgut, salivary glands epithelia and the hemocoel. Molecularly, invertebrate immunity is initiated when pathogen recognition molecules trigger serum or cellular signalling cascades leading to the production of antimicrobials, pathogen opsonization and phagocytosis. We presently aimed at identifying hemocyte transcripts from semi-engorged female I. ricinus ticks by mass sequencing a hemocyte cDNA library and annotating immune-related transcripts based on their hemocyte abundance as well as their ubiquitous distribution.

METHODOLOGY/PRINCIPAL FINDINGS

De novo assembly of 926,596 pyrosequence reads plus 49,328,982 Illumina reads (148 nt length) from a hemocyte library, together with over 189 million Illumina reads from salivary gland and midgut libraries, generated 15,716 extracted coding sequences (CDS); these are displayed in an annotated hyperlinked spreadsheet format. Read mapping allowed the identification and annotation of tissue-enriched transcripts. A total of 327 transcripts were found significantly over expressed in the hemocyte libraries, including those coding for scavenger receptors, antimicrobial peptides, pathogen recognition proteins, proteases and protease inhibitors. Vitellogenin and lipid metabolism transcription enrichment suggests fat body components. We additionally annotated ubiquitously distributed transcripts associated with immune function, including immune-associated signal transduction proteins and transcription factors, including the STAT transcription factor.

CONCLUSIONS/SIGNIFICANCE

This is the first systems biology approach to describe the genes expressed in the haemocytes of this neglected disease vector. A total of 2,860 coding sequences were deposited to GenBank, increasing to 27,547 the number so far deposited by our previous transcriptome studies that serves as a discovery platform for studies with I. ricinus biochemistry and physiology.

Thioester-containing proteins of the tick Ixodes ricinus: gene expression, response to microbial challenge and their role in phagocytosis of the yeast Candida albicans

The ability of ticks to act as vectors for a wide range of serious human and animal infectious diseases is apparently linked to the insufficiency of the tick immune system to effectively eliminate pathogens they transmit. At the tick-pathogen interface, an important role is presumably played by components of an ancient complement system that includes a repertoire of thioester-containing proteins (TEPs), which in Ixodes sp. comprises three α2-macroglobulins (A2M), three C3 complement component-related molecules (C3), two macroglobulin complement-related (Mcr) and one insect-type TEPs (Tep). In order to assess the function of TEPs in tick immunity, a quantitative real-time PCR expression analysis of tick TEPs was performed at various developmental stages of Ixodes ricinus, and in tissues dissected from adult females. Expression of TEP genes was mostly tissue specific; IrA2M1, IrC3-1, IrC3-3 were found to be expressed in cells of tick fat body adjacent to the tracheal trunks, IrA2M2 in hemocytes, IrTep in ovaries, IrMcr1 in salivary glands and only IrA2M3, IrC3-2 and IrMcr2 mRNAs were present in multiple organs. Expression of tick TEPs was further examined in response to injection of model microbes representing Gram-negative, Gram-positive bacteria and yeast. The greatest expression induction was observed for IrA2M1 and IrC3-1 after challenge with the yeast Candida albicans. Phagocytosis of the yeast was strongly dependent on an active thioester bond and the subsequent silencing of individual tick TEPs by RNA interference demonstrated the involvement of IrC3-1 and IrMcr2. This result suggests the existence of a distinct complement-like pathway, different from that leading to phagocytosis of Gram-negative bacteria. Understanding of the tick immune response against model microbes should provide new concepts for investigating interactions between ticks and relevant tick-borne pathogens.

Identification of host proteins involved in rickettsial invasion of tick cells

Tick-borne spotted fever group (SFG) Rickettsia species are obligate intracellular bacteria capable of infecting both vertebrate and invertebrate host cells, an essential process for subsequent bacterial survival in distinct hosts. The host cell signaling molecules involved in the uptake of Rickettsia into mammalian and Drosophila cells have been identified; however, invasion into tick cells is understudied. Considering the movement of SFG Rickettsia between vertebrate and invertebrate hosts, the hypothesis is that conserved mechanisms are utilized for host cell invasion. The current study employed biochemical inhibition assays to determine the tick proteins involved in Rickettsia montanensis infection of tick-derived cells from a natural host, Dermacentor variabilis. The results revealed several tick proteins important for rickettsial invasion, including actin filaments, actin-related protein 2/3 complex, phosphatidylinositol-3'-kinase, protein tyrosine kinases (PTKs), Src family PTK, focal adhesion kinase, Rho GTPase Rac1, and neural Wiskott-Aldrich syndrome protein. Delineating the molecular mechanisms of rickettsial infection is critical to a thorough understanding of rickettsial transmission in tick populations and the ecology of tick-borne rickettsial diseases.

De novo transcriptome of the Hemimetabolous German cockroach (Blattella germanica)

Background

The German cockroach, Blattella germanica, is an important insect pest that transmits various pathogens mechanically and causes severe allergic diseases. This insect has long served as a model system for studies of insect biology, physiology and ecology. However, the lack of genome or transcriptome information heavily hinder our further understanding about the German cockroach in every aspect at a molecular level and on a genome-wide scale. To explore the transcriptome and identify unique sequences of interest, we subjected the B. germanica transcriptome to massively parallel pyrosequencing and generated the first reference transcriptome for B. germanica.

Methodology/Principal Findings

A total of 1,365,609 raw reads with an average length of 529 bp were generated via pyrosequencing the mixed cDNA library from different life stages of German cockroach including maturing oothecae, nymphs, adult females and males. The raw reads were de novo assembled to 48,800 contigs and 3,961 singletons with high-quality unique sequences. These sequences were annotated and classified functionally in terms of BLAST, GO and KEGG, and the genes putatively coding detoxification enzyme systems, insecticide targets, key components in systematic RNA interference, immunity and chemoreception pathways were identified. A total of 3,601 SSRs (Simple Sequence Repeats) loci were also predicted.

Conclusions/Significance

The whole transcriptome pyrosequencing data from this study provides a usable genetic resource for future identification of potential functional genes involved in various biological processes.

Immunity-related genes in Ixodes scapularis--perspectives from genome information

Ixodes scapularis, commonly known as the deer tick, transmits a wide array of human and animal pathogens including Borrelia burgdorferi. Despite substantial advances in our understanding of immunity in model arthropods, including other disease vectors, precisely how I. scapularis immunity functions and influences persistence of invading pathogens remains largely unknown. This review provides a comprehensive analysis of the recently sequenced I. scapularis genome for the occurrence of immune-related genes and related pathways. We will also discuss the potential influence of immunity-related genes on the persistence of tick-borne pathogens with an emphasis on the Lyme disease pathogen B. burgdorferi. Further enhancement of our knowledge of tick immune responses is critical to understanding the molecular basis of the persistence of tick-borne pathogens and development of novel interventions against the relevant infections.

IrSPI, a tick serine protease inhibitor involved in tick feeding and Bartonella henselae infection

Ixodes ricinus is the most widespread and abundant tick in Europe, frequently bites humans, and is the vector of several pathogens including those responsible for Lyme disease, Tick-Borne Encephalitis, anaplasmosis, babesiosis and bartonellosis. These tick-borne pathogens are transmitted to vertebrate hosts via tick saliva during blood feeding, and tick salivary gland (SG) factors are likely implicated in transmission. In order to identify such tick factors, we characterized the transcriptome of female I. ricinus SGs using next generation sequencing techniques, and compared transcriptomes between Bartonella henselae-infected and non-infected ticks. High-throughput sequencing of I. ricinus SG transcriptomes led to the generation of 24,539 isotigs. Among them, 829 and 517 transcripts were either significantly up- or down-regulated respectively, in response to bacterial infection. Searches based on sequence identity showed that among the differentially expressed transcripts, 161 transcripts corresponded to nine groups of previously annotated tick SG gene families, while the others corresponded to genes of unknown function. Expression patterns of five selected genes belonging to the BPTI/Kunitz family of serine protease inhibitors, the tick salivary peptide group 1 protein, the salp15 super-family, and the arthropod defensin family, were validated by qRT-PCR. IrSPI, a member of the BPTI/Kunitz family of serine protease inhibitors, showed the highest up-regulation in SGs in response to Bartonella infection. IrSPI silencing impaired tick feeding, as well as resulted in reduced bacterial load in tick SGs. This study provides a comprehensive analysis of I. ricinus SG transcriptome and contributes significant genomic information about this important disease vector. This in-depth knowledge will enable a better understanding of the molecular interactions between ticks and tick-borne pathogens, and identifies IrSPI, a candidate to study now in detail to estimate its potentialities as vaccine against the ticks and the pathogens they transmit.

Identification of Rhipicephalus microplus genes that modulate the infection rate of the rickettsia Anaplasma marginale

Arthropod vectors transmit a diversity of animal and human pathogens, ranging from RNA viruses to protozoal parasites. Chemotherapeutic control of pathogens has classically focused either on insecticides that kill the vector itself or antimicrobials for infected patients. The limitation of the former is that it targets both infected and uninfected vectors and selects for resistant populations while the latter requires prompt and accurate diagnosis. An alternative strategy is to target vector molecules that permit the pathogen to establish itself, replicate, and/or develop within the vector. Using the rickettsial pathogen Anaplasma marginale and its tropical tick vector, Rhipicephalus microplus, as a model, we tested whether silencing specific gene targets would affect tick infection rates (the % of fed ticks that are infected with the pathogen) and pathogen levels within infected ticks. Silencing of three R. microplus genes, CK187220, CV437619 and TC18492, significantly decreased the A. marginale infection rate in salivary glands, whereas gene silencing of TC22382, TC17129 and TC16059 significantly increased the infection rate in salivary glands. However in all cases of significant difference in the infection rate, the pathogen levels in the ticks that did become infected, were not significantly different. These results are consistent with the targeted genes affecting the pathogen at early steps in infection of the vector rather than in replication efficiency. Identifying vector genes and subsequent determination of the encoded functions are initial steps in discovery of new targets for inhibiting pathogen development and subsequent transmission.

Strategies and molecular tools to fight antimicrobial resistance: resistome, transcriptome, and antimicrobial peptides

The increasing number of antibiotic resistant bacteria motivates prospective research toward discovery of new antimicrobial active substances. There are, however, controversies concerning the cost-effectiveness of such research with regards to the description of new substances with novel cellular interactions, or description of new uses of existing substances to overcome resistance. Although examination of bacteria isolated from remote locations with limited exposure to humans has revealed an absence of antibiotic resistance genes, it is accepted that these genes were both abundant and diverse in ancient living organisms, as detected in DNA recovered from Pleistocene deposits (30,000 years ago). Indeed, even before the first clinical use of antibiotics more than 60 years ago, resistant organisms had been isolated. Bacteria can exhibit different strategies for resistance against antibiotics. New genetic information may lead to the modification of protein structure affecting the antibiotic carriage into the cell, enzymatic inactivation of drugs, or even modification of cellular structure interfering in the drug-bacteria interaction. There are still plenty of new genes out there in the environment that can be appropriated by putative pathogenic bacteria to resist antimicrobial agents. On the other hand, there are several natural compounds with antibiotic activity that may be used to oppose them. Antimicrobial peptides (AMPs) are molecules which are wide-spread in all forms of life, from multi-cellular organisms to bacterial cells used to interfere with microbial growth. Several AMPs have been shown to be effective against multi-drug resistant bacteria and have low propensity to resistance development, probably due to their unique mode of action, different from well-known antimicrobial drugs. These substances may interact in different ways with bacterial cell membrane, protein synthesis, protein modulation, and protein folding. The analysis of bacterial transcriptome may contribute to the understanding of microbial strategies under different environmental stresses and allows the understanding of their interaction with novel AMPs.

The ovarian transcriptome of the cattle tick, Rhipicephalus (Boophilus) microplus, feeding upon a bovine host infected with Babesia bovis

BACKGROUND

Cattle babesiosis is a tick-borne disease of cattle with the most severe form of the disease caused by the apicomplexan, Babesia bovis. Babesiosis is transmitted to cattle through the bite of infected cattle ticks of the genus Rhipicephalus. The most prevalent species is Rhipicephalus (Boophilus) microplus, which is distributed throughout the tropical and subtropical countries of the world. The transmission of B. bovis is transovarian and a previous study of the R. microplus ovarian proteome identified several R. microplus proteins that were differentially expressed in response to infection. Through various approaches, we studied the reaction of the R. microplus ovarian transcriptome in response to infection by B. bovis.

METHODS

A group of ticks were allowed to feed on a B. bovis-infected splenectomized calf while a second group fed on an uninfected splenectomized control calf. RNA was purified from dissected adult female ovaries of both infected and uninfected ticks and a subtracted B. bovis-infected cDNA library was synthesized, subtracting with the uninfected ovarian RNA. Four thousand ESTs were sequenced from the ovary subtracted library and annotated.

RESULTS

The subtracted library dataset assembled into 727 unique contigs and 2,161 singletons for a total of 2,888 unigenes, Microarray experiments designed to detect B. bovis-induced gene expression changes indicated at least 15 transcripts were expressed at a higher level in ovaries from ticks feeding upon the B. bovis-infected calf as compared with ovaries from ticks feeding on an uninfected calf. We did not detect any transcripts from these microarray experiments that were expressed at a lower level in the infected ovaries compared with the uninfected ovaries. Using the technique called serial analysis of gene expression, 41 ovarian transcripts from infected ticks were differentially expressed when compared with transcripts of controls.

CONCLUSION

Collectively, our experimental approaches provide the first comprehensive profile of the R. microplus ovarian transcriptome responding to infection by B. bovis. This dataset should prove useful in molecular studies of host-pathogen interactions between this tick and its apicomplexan parasite.

Interaction of the tick immune system with transmitted pathogens

Ticks are hematophagous arachnids transmitting a wide variety of pathogens including viruses, bacteria, and protozoans to their vertebrate hosts. The tick vector competence has to be intimately linked to the ability of transmitted pathogens to evade tick defense mechanisms encountered on their route through the tick body comprising midgut, hemolymph, salivary glands or ovaries. Tick innate immunity is, like in other invertebrates, based on an orchestrated action of humoral and cellular immune responses. The direct antimicrobial defense in ticks is accomplished by a variety of small molecules such as defensins, lysozymes or by tick-specific antimicrobial compounds such as microplusin/hebraein or 5.3-kDa family proteins. Phagocytosis of the invading microbes by tick hemocytes is likely mediated by the primordial complement-like system composed of thioester-containing proteins, fibrinogen-related lectins and convertase-like factors. Moreover, an important role in survival of the ingested microbes seems to be played by host proteins and redox balance maintenance in the tick midgut. Here, we summarize recent knowledge about the major components of tick immune system and focus on their interaction with the relevant tick-transmitted pathogens, represented by spirochetes (Borrelia), rickettsiae (Anaplasma), and protozoans (Babesia). Availability of the tick genomic database and feasibility of functional genomics based on RNA interference greatly contribute to the understanding of molecular and cellular interplay at the tick-pathogen interface and may provide new targets for blocking the transmission of tick pathogens.

Gut transcriptome of replete adult female cattle ticks, Rhipicephalus (Boophilus) microplus, feeding upon a Babesia bovis-infected bovine host

As it feeds upon cattle, Rhipicephalus (Boophilus) microplus is capable of transmitting a number of pathogenic organisms, including the apicomplexan hemoparasite Babesia bovis, a causative agent of bovine babesiosis. The R. microplus female gut transcriptome was studied for two cohorts: adult females feeding on a bovine host infected with B. bovis and adult females feeding on an uninfected bovine. RNA was purified and used to generate a subtracted cDNA library from B. bovis-infected female gut, and 4,077 expressed sequence tags (ESTs) were sequenced. Gene expression was also measured by a microarray designed from the publicly available R. microplus gene index: BmiGI Version 2. We compared gene expression in the tick gut from females feeding upon an uninfected bovine to gene expression in tick gut from females feeding upon a splenectomized bovine infected with B. bovis. Thirty-three ESTs represented on the microarray were expressed at a higher level in female gut samples from the ticks feeding upon a B. bovis-infected calf compared to expression levels in female gut samples from ticks feeding on an uninfected calf. Forty-three transcripts were expressed at a lower level in the ticks feeding upon B. bovis-infected female guts compared with expression in female gut samples from ticks feeding on the uninfected calf. These array data were used as initial characterization of gene expression associated with the infection of R. microplus by B. bovis.

Nutrigenomics in Arma chinensis: transcriptome analysis of Arma chinensis fed on artificial diet and Chinese oak silk moth Antheraea pernyi pupae

Background

The insect predator, Arma chinensis, is capable of effectively controlling many pests, such as Colorado potato beetle, cotton bollworm, and mirid bugs. Our previous study demonstrated several life history parameters were diminished for A. chinensis reared on an artificial diet compared to a natural food source like the Chinese oak silk moth pupae. The molecular mechanisms underlying the nutritive impact of the artificial diet on A. chinensis health are unclear. So we utilized transcriptome information to better understand the impact of the artificial diet on A. chinensis at the molecular level.

Methodology/Principal Findings

Illumina HiSeq2000 was used to sequence 4.79 and 4.70 Gb of the transcriptome from pupae-fed and artificial diet-fed A. chinensis libraries, respectively, and a de novo transcriptome assembly was performed (Trinity short read assembler). This resulted in 112,029 and 98,724 contigs, clustered into 54,083 and 54,169 unigenes for pupae-fed and diet-fed A. chinensis, respectively. Unigenes from each sample’s assembly underwent sequence splicing and redundancy removal to acquire non-redundant unigenes. We obtained 55,189 unigenes of A. chinensis, including 12,046 distinct clusters and 43,143 distinct singletons. Unigene sequences were aligned by BLASTx to nr, Swiss-Prot, KEGG and COG (E-value <10⁻⁵), and further aligned by BLASTn to nt (E-value <10⁻⁵), retrieving proteins of highest sequence similarity with the given unigenes along with their protein functional annotations. Totally, 22,964, 7,898, 18,069, 15,416, 8,066 and 5,341 unigenes were annotated in nr, nt, Swiss-Prot, KEGG, COG and GO, respectively. We compared gene expression variations and found thousands of genes were differentially expressed between pupae-fed and diet-fed A. chinensis.

Conclusions/Significance

Our study provides abundant genomic data and offers comprehensive sequence information for studying A. chinensis. Additionally, the physiological roles of the differentially expressed genes enable us to predict effects of some dietary ingredients and subsequently propose formulation improvements to artificial diets.

Ixodes scapularis JAK-STAT pathway regulates tick antimicrobial peptides, thereby controlling the agent of human granulocytic anaplasmosis

Ixodes scapularis transmits the agent of human granulocytic anaplasmosis, among other pathogens. The mechanisms used by the tick to control Anaplasma phagocytophilum are not known. We demonstrate that the I. scapularis Janus kinase (JAK)-signaling transducer activator of transcription (STAT) pathway plays a critical role in A. phagocytophilum infection of ticks. The A. phagocytophilum burden increases in salivary glands and hemolymph when the JAK-STAT pathway is suppressed by RNA interference. The JAK-STAT pathway exerts its anti-Anaplasma activity presumably through STAT-regulated effectors. A salivary gland gene family encoding 5.3-kDa antimicrobial peptides is highly induced upon A. phagocytophilum infection of tick salivary glands. Gene expression and electrophoretic mobility shift assays showed that the 5.3-kDa antimicrobial peptide-encoding genes are regulated by tick STAT. Silencing of these genes increased A. phagocytophilum infection of tick salivary glands and transmission to mammalian host. These data suggest that the JAK-STAT signaling pathway plays a key role in controlling A. phagocytophilum infection in ticks by regulating the expression of antimicrobial peptides.

Pyrosequencing the Bemisia tabaci transcriptome reveals a highly diverse bacterial community and a robust system for insecticide resistance

Background

Bemisia tabaci (Gennadius) is a phloem-feeding insect poised to become one of the major insect pests in open field and greenhouse production systems throughout the world. The high level of resistance to insecticides is a main factor that hinders continued use of insecticides for suppression of B. tabaci. Despite its prevalence, little is known about B. tabaci at the genome level. To fill this gap, an invasive B. tabaci B biotype was subjected to pyrosequencing-based transcriptome analysis to identify genes and gene networks putatively involved in various physiological and toxicological processes.

Methodology and Principal Findings

Using Roche 454 pyrosequencing, 857,205 reads containing approximately 340 megabases were obtained from the B. tabaci transcriptome. De novo assembly generated 178,669 unigenes including 30,980 from insects, 17,881 from bacteria, and 129,808 from the nohit. A total of 50,835 (28.45%) unigenes showed similarity to the non-redundant database in GenBank with a cut-off E-value of 10–5. Among them, 40,611 unigenes were assigned to one or more GO terms and 6,917 unigenes were assigned to 288 known pathways. De novo metatranscriptome analysis revealed highly diverse bacterial symbionts in B. tabaci, and demonstrated the host-symbiont cooperation in amino acid production. In-depth transcriptome analysis indentified putative molecular markers, and genes potentially involved in insecticide resistance and nutrient digestion. The utility of this transcriptome was validated by a thiamethoxam resistance study, in which annotated cytochrome P450 genes were significantly overexpressed in the resistant B. tabaci in comparison to its susceptible counterparts.

Conclusions

This transcriptome/metatranscriptome analysis sheds light on the molecular understanding of symbiosis and insecticide resistance in an agriculturally important phloem-feeding insect pest, and lays the foundation for future functional genomics research of the B. tabaci complex. Moreover, current pyrosequencing effort greatly enriched the existing whitefly EST database, and makes RNAseq a viable option for future genomic analysis.

Functional characterization of two defensin isoforms of the hard tick Ixodes ricinus

BACKGROUND

The immune system of ticks is stimulated to produce many pharmacologically active molecules during feeding and especially during pathogen invasion. The family of cationic peptides - defensins - represents a specific group of antimicrobial compounds with six conserved cysteine residues in a molecule.

RESULTS

Two isoforms of the defensin gene (def1 and def2) were identified in the European tick Ixodes ricinus. Expression of both genes was induced in different tick organs by a blood feeding or pathogen injection. We have tested the ability of synthetic peptides def1 and def2 to inhibit the growth or directly kill several pathogens. The antimicrobial activities (expressed as minimal inhibition concentration and minimal bactericidal concentration values) against Gram positive bacteria were confirmed, while Gram negative bacteria, yeast, Tick Borne Encephalitis and West Nile Viruses were shown to be insensitive. In addition to antimicrobial activities, the hemolysis effect of def1 and def2 on human erythrocytes was also established.

CONCLUSIONS

Although there is nothing known about the realistic concentration of defensins in I. ricinus tick body, these results suggest that defensins play an important role in defence against different pathogens. Moreover this is a first report of a one amino acid substitution in a defensins molecule and its impact on antimicrobial activity.

Transcriptomics of the bed bug (Cimex lectularius)

BACKGROUND

Bed bugs (Cimex lectularius) are blood-feeding insects poised to become one of the major pests in households throughout the United States. Resistance of C. lectularius to insecticides/pesticides is one factor thought to be involved in its sudden resurgence. Despite its high-impact status, scant knowledge exists at the genomic level for C. lectularius. Hence, we subjected the C. lectularius transcriptome to 454 pyrosequencing in order to identify potential genes involved in pesticide resistance.

METHODOLOGY AND PRINCIPAL FINDINGS

Using 454 pyrosequencing, we obtained a total of 216,419 reads with 79,596,412 bp, which were assembled into 35,646 expressed sequence tags (3902 contigs and 31744 singletons). Nearly 85.9% of the C. lectularius sequences showed similarity to insect sequences, but 44.8% of the deduced proteins of C. lectularius did not show similarity with sequences in the GenBank non-redundant database. KEGG analysis revealed putative members of several detoxification pathways involved in pesticide resistance. Lamprin domains, Protein Kinase domains, Protein Tyrosine Kinase domains and cytochrome P450 domains were among the top Pfam domains predicted for the C. lectularius sequences. An initial assessment of putative defense genes, including a cytochrome P450 and a glutathione-S-transferase (GST), revealed high transcript levels for the cytochrome P450 (CYP9) in pesticide-exposed versus pesticide-susceptible C. lectularius populations. A significant number of single nucleotide polymorphisms (296) and microsatellite loci (370) were predicted in the C. lectularius sequences. Furthermore, 59 putative sequences of Wolbachia were retrieved from the database.

CONCLUSIONS

To our knowledge this is the first study to elucidate the genetic makeup of C. lectularius. This pyrosequencing effort provides clues to the identification of potential detoxification genes involved in pesticide resistance of C. lectularius and lays the foundation for future functional genomics studies.