De novo assembly and analysis of midgut transcriptome of Haemaphysalis flava and identification of genes involved in blood digestion, feeding and defending from pathogens

Functional annotation and analysis of the hard tick Dermacentor nuttalli midgut genes

Ticks are hematophagous vectors that transmit a variety of pathogens, posing significant threats to the health of both humans and animals. Tick midgut proteins play essential roles in blood digestion, feeding, toxic waste processing, and pathogen transmission. Dermacentor nuttalli is the primary vector of tick-borne pathogens, including rickettsioses in the Qinghai-Tibet Plateau. However, there is a lack of genomic, transcriptomic, and proteomic information regarding the biology of D. nuttalli. In this study, we assembled and compared the midgut transcriptomes of female D. nuttalli ticks at 0, 24, 48, 72, and 96 h during blood feeding, identifying the genes with differentially regulated expression following feeding. The obtained data were compiled and annotated in multiple databases including Nr, NT, PFAM, KOG, KEGG, and GO. The high-quality clean readings of midgut tissue at the different blood-feeding times were recorded as 22,524,912, 23,752,325, 20,377,718, 21,300,710, and 20,378,658, respectively. The transcripts were classified into eight large categories, including immunogenic proteases (8.37%), protease inhibitors (0.85%), transporters (3.96%), ligand binding proteins (1.98%), ribosomal function proteins (0.94%), heat shock proteins (0.30%), other proteases and miscellaneous proteins (57.61%), and unknown proteins (26.00%). Significant differences were observed in the genes obtained at 0, 24, 48, 72, and 96 h during blood feeding. The differentially expressed genes include catalytic proteins that play an important role in accelerating biochemical reactions, binding activity proteins which are involved in various molecular interactions, and proteins that actively participate in multiple metabolic pathways and cellular processes. Notably, the gene expression in the midgut of D. nuttalli shows dynamic changes every 24 h throughout the blood-feeding process. This change may represent an equivalent strategy of antigenic variation for ticks, designed to protect their essential feeding function against the host's immune system. The tick antigens identified in this study may serve as promising candidates for the development of effective vaccines or as drug targets for acaricides.

The proteomic content of Varroa destructor gut varies according to the developmental stage of its host

The nutritional physiology of parasites is often overlooked although it is at the basis of host-parasite interactions. In the case of Varroa destructor, one of the major pests of the Western honey bee Apis mellifera, the nature of molecules and tissues ingested by the parasite is still not completely understood. Here, the V. destructor feeding biology was explored through artificial feeding, dissection of the mite's gut and proteomic analyses. More specifically, the proteome of guts extracted from starved mites and honey bee-fed mites was compared to highlight both the parasite proteins likely involved in food processing and the honey bee proteins actually ingested by the mite. We could identify 25 V. destructor candidate proteins likely involved in the parasite digestion. As the host developmental stages infested by the mite are diverse, we also focused on the identity and on the origin of honey bee proteins ingested by the mite when it feeds on larvae, pupae or adults. We highlighted profiles of consumed honey bee proteins and their variations throughout the V. destructor life cycle. These variations matched the ones observed in the honey bee hemolymph, showing that this tissue is an important part of the mite's diet. Based on the variations of abundance of the most consumed honey bee proteins and on their functions, the potential implication of these key candidate nutrients in V. destructor reproduction is also discussed.

A multi-omics approach for understanding blood digestion dynamics in Ixodes scapularis and identification of anti-tick vaccine targets

Ixodes scapularis, the black-legged tick, is a major arthropod vector that transmits the causative agents of Lyme disease and several other pathogens of human significance. The tick midgut is the main tissue involved in blood acquisition and digestion and the first organ to have contact with pathogens ingested through the blood meal. Gene expression in the midgut before, during, and after a blood meal may vary in response to the physiological changes due to blood feeding. A systems biology approach based on RNA and protein sequencing was used to gain insight into the changes in tick midgut transcripts and proteins during blood ingestion (unfed and partially fed) and digestion (1-, 2-, 7-, and 14 days post detachment from the host) by the Ixodes scapularis female ticks. A total of 2,726 differentially expressed transcripts, and 449 proteins were identified across the time points. Genes involved in detoxification of xenobiotics, proteases, protease inhibitors, metabolism, and immunity were differentially expressed in response to blood feeding. Similarly, proteins corresponding to the same groups were also differentially expressed. Nine genes from major gene categories were chosen as potential vaccine candidates, and, using RNA interference, the effect of these gene knockdowns on tick biology was investigated. Knockdown of these genes had variable negative impacts on tick physiology, such as the inability to engorge fully and to produce eggs and increased mortality. These and additional gene targets provide opportunities to explore novel tick control strategies.

Transcriptome analysis of Haemaphysalis flava female using Illumina HiSeq 4000 sequencing: de novo assembly, functional annotation and discovery of SSR markers

BACKGROUND

Ticks are ectoparasites capable of directly damaging their hosts and transmitting vector-borne diseases. The ixodid tick Haemaphysalis flava has a broad distribution that extends from East to South Asia. This tick is a reservoir of severe fever with thrombocytopenia syndrome virus (SFTSV) that causes severe hemorrhagic disease, with cases reported from China, Japan and South Korea. Recently, the distribution of H. flava in South Korea was found to overlap with the occurrence of SFTSV.

METHODS

This study was undertaken to discover the molecular resources of H. flava female ticks using the Illumina HiSeq 4000 system, the Trinity de novo sequence assembler and annotation against public databases. The locally curated Protostome database (PANM-DB) was used to screen the putative adaptation-related transcripts classified to gene families, such as angiotensin-converting enzyme, aquaporin, adenylate cyclase, AMP-activated protein kinase, glutamate receptors, heat shock proteins, molecular chaperones, insulin receptor, mitogen-activated protein kinase and solute carrier family proteins. Also, the repeats and simple sequence repeats (SSRs) were screened from the unigenes using RepeatMasker (v4.0.6) and MISA (v1.0) software tools, followed by the designing of SSRs flanking primers using BatchPrimer 3 (v1.0) software.

RESULTS

The transcriptome produced a total of 69,822 unigenes, of which 46,175 annotated to the homologous proteins in the PANM-DB. The unigenes were also mapped to the EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) specializations. Promiscuous presence of protein kinase, zinc finger (C2H2-type), reverse transcriptase, and RNA recognition motif domains was observed in the unigenes. A total of 3480 SSRs were screened, of which 1907 and 1274 were found as tri- and dinucleotide repeats, respectively. A list of primer sequences flanking the SSR motifs was detailed for validation of polymorphism in H. flava and the related tick species.

CONCLUSIONS

The reference transcriptome information on H. flava female ticks will be useful for an enriched understanding of tick biology, its competency to act as a vector and the study of species diversity related to disease transmission.

De novo assembly and annotation of the Amblyomma hebraeum tick midgut transcriptome response to Ehrlichia ruminantium infection

The South African bont tick Amblyomma hebraeum is a hematophagous vector for the heartwater disease pathogen Ehrlichia ruminantium in southern Africa. During feeding, the tick's enterocytes express proteins that perform vital functions in blood digestion, including proteins that may be involved in E. ruminantium acquisition, colonization or immunity. To delineate the molecular mechanism of midgut response to E. ruminantium infection, we performed comparative analyses of midgut transcriptomes of E. ruminantium infected engorged A. hebraeum nymphs, and infected adult male and female ticks with their corresponding matched uninfected controls, before and during feeding. A total of 102,036 unigenes were annotated in public databases and their expression levels analyzed for engorged nymphs as well as unfed and partly-fed adult ticks. There were 2,025 differentially expressed genes (DEGs) in midguts, of which 1,225 unigenes were up-regulated and 800 unigenes were down-regulated in the midguts of infected ticks. Annotation of DEGs revealed an increase in metabolic and cellular processes among E. ruminantium infected ticks. Notably, among the infected ticks, there was up-regulation in the expression of genes involved in tick immunity, histone proteins and oxidative stress responses. We also observed up-regulation of glycoproteins that E. ruminantium could potentially use as docking sites for host cell entry. Insights uncovered in this study offer a platform for further investigations into the molecular interaction between E. ruminantium and A. hebraeum.

Integrative analysis highlights molecular and immune responses of tick Amblyomma americanum to Escherichia coli challenge

Ticks are ectoparasites that can transmit various pathogens capable of causing life-threatening illnesses in people and animals, making them a severe public health threat. Understanding how ticks respond to bacterial infection is crucial for deciphering their immune defense mechanisms and identifying potential targets for controlling tick-borne diseases. In this study, an in-depth transcriptome analysis was used to investigate the molecular and immune responses of Amblyomma americanum to infection caused by the microinjection of Escherichia coli. With an abundance of differentially expressed genes discovered at different times, the analysis demonstrated significant changes in gene expression profiles in response to E. coli challenge. Notably, we found alterations in crucial immune markers, including the antimicrobial peptides defensin and microplusin, suggesting they may play an essential role in the innate immune response. Furthermore, KEGG analysis showed that following E. coli exposure, a number of key enzymes, including lysosomal alpha-glucosidase, fibroblast growth factor, legumain, apoptotic protease-activating factor, etc., were altered, impacting the activity of the lysosome, mitogen-activated protein kinase, antigen processing and presentation, bacterial invasion, apoptosis, and the Toll and immune deficiency pathways. In addition to the transcriptome analysis, we constructed protein interaction networks to elucidate the molecular interactions underlying the tick's response to E. coli challenge. Hub genes were identified, and their functional enrichment provided insights into the regulation of cytoskeleton rearrangement, apoptotic processes, and kinase activity that may occur in infected cells. Collectively, the findings shed light on the potential immune responses in A. americanum that control E. coli infection.

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.

Effect of food source availability in the salivary gland transcriptome of the unique burying beetle Nicrophorus pustulatus (Coleoptera: Silphidae)

Nicrophorus is a genus of beetles that bury and transform small vertebrate carcasses into a brood ball coated with their oral and anal secretions to prevent decay and that will serve as a food source for their young. Nicrophorus pustulatus is an unusual species with the ability to overtake brood of other burying beetles and whose secretions, unlike other Nicrophorus species, has been reported not to exhibit antimicrobial properties. This work aims to better understand how the presence or absence of a food source influences the expression of genes involved in the feeding process of N. pustulatus. To achieve that, total RNA was extracted from pooled samples of salivary gland tissue from N. pustulatus and sequenced using an Illumina platform. The resulting reads were used to assemble a de novo transcriptome using Trinity. Duplicates with more than 95% similarity were removed to obtain a "unigene" set. Annotation of the unigene set was done using the Trinotate pipeline. Transcript abundance was determined using Kallisto and differential gene expression analysis was performed using edgeR. A total of 651 genes were found to be differentially expressed, including 390 upregulated and 261 downregulated genes in fed insects compared to starved. Several genes upregulated in fed beetles are associated with the insect immune response and detoxification processes with only one transcript encoding for the antimicrobial peptide (AMP) defensin. These results confirm that N. pustulatus does not upregulate the production of genes encoding AMPs during feeding. This study provides a snapshot of the changes in gene expression in the salivary glands of N. pustulatus following feeding while providing a well described transcriptome for the further analysis of this unique burying beetle.

Characterization of AV422 from Haemaphysalis flava ticks in vitro

Ticks are hematophagous ectoparasites and cause a major public health threat worldwide. Development of anti-tick vaccines is regarded to be an optimal alternative for tick control. AV422, a unique protein in ticks, is secreted into hosts during blood-feeding, but its roles are not confirmed in Haemaphysalis flava ticks. We retrieved a gene fragment encoding AV422 from a transcriptome dataset of H. flava, and based on it, we reconstructed the full length of AV422 from H. flava (Hf-AV422) by rapid amplification of cDNA ends. Expression profiles of Hf-AV422 in whole ticks and organs of different engorgement levels were determined by qPCR. Then its opening reading frame (ORF) was expressed in Escherichia coli strain BL21 (DE3). The prothrombin time (PT), activated partial thromboplastin time (APTT) and thrombin time (TT) assays were conducted to test anticoagulant activities of the purified recombinant protein (rHf-AV422). The full length of AV422 was 1152 bp. Hf-AV422 showed to be conserved as indicated by multiple sequence alignment. Expression of Hf-AV422 was significantly higher in salivary glands and cuticles than in ovaries. Its expression in whole ticks decreased during engorgement with the highest levels in 1/4 engorged ticks. rHf-AV422 prolonged PT, APTT and TT when incubated with rabbit plasma. Our data demonstrated that Hf-AV422 is a conserved salivary protein with anticoagulant activity. Further studies are needed to test in detail its functional properties to ensure it an adequate antigen candidate for the development of broad-spectrum vaccines against ticks.

De novo transcriptome sequencing and comparative profiling of the ovary in partially engorged and fully engorged Haemaphysalis flava ticks

Haemaphysalis flava is the vector of several pathogens and has medical and veterinary importance. Transcriptome information of the ovary of H. flava is unavailable and limits understanding of its molecular basis of reproduction. We studied the ovary transcriptome of partially engorged and fully engorged H. flava using high-throughput RNA sequencing technology. A total of 53,025,360 and 57,942,890 clean reads were obtained with 7.95 GB and 8.69 GB clean bases in partially engorged ticks (PETs) and fully engorged ticks (FETs), respectively. The clean reads were assembled into 138,711 unigenes. A total of 72,043 unigenes (51.93%) were annotated and 66,668 unigenes (48.07%) were unknown. A total of 38,487 differentially expressed genes (DEGs) were found between PET and FET with 19,031 upregulated genes and 19,456 downregulated genes. The RNA-seq results were validated by qRT-PCR, including six upregulated genes and three downregulated genes. Some unigenes coding for nutrient transporters, proteases, and protease inhibitors were found and analyzed. This study was the first time to perform the transcriptome sequences of the ovary of partially engorged and fully engorged H. flava. The results can benefit the understanding of the molecular basis of ovary maturation and oogenesis of the H. flava and boost the development of the strategies for control of H. flava.

Borrelia burgdorferi infection modifies protein content in saliva of Ixodes scapularis nymphs

BACKGROUND

Lyme disease (LD) caused by Borrelia burgdorferi is the most prevalent tick-borne disease. There is evidence that vaccines based on tick proteins that promote tick transmission of B. burgdorferi could prevent LD. As Ixodes scapularis nymph tick bites are responsible for most LD cases, this study sought to identify nymph tick saliva proteins associated with B. burgdorferi transmission using LC-MS/MS. Tick saliva was collected using a non-invasive method of stimulating ticks (uninfected and infected: unfed, and every 12 h during feeding through 72 h, and fully-fed) to salivate into 2% pilocarpine-PBS for protein identification using LC-MS/MS.

RESULTS

We identified a combined 747 tick saliva proteins of uninfected and B. burgdorferi infected ticks that were classified into 25 functional categories: housekeeping-like (48%), unknown function (18%), protease inhibitors (9%), immune-related (6%), proteases (8%), extracellular matrix (7%), and small categories that account for <5% each. Notably, B. burgdorferi infected ticks secreted high number of saliva proteins (n=645) than uninfected ticks (n=376). Counter-intuitively, antimicrobial peptides, which function to block bacterial infection at tick feeding site were suppressed 23-85 folds in B. burgdorferi infected ticks. Similar to glycolysis enzymes being enhanced in mammalian cells exposed to B. burgdorferi : eight of the 10-glycolysis pathway enzymes were secreted at high abundance by B. burgdorferi infected ticks. Of significance, rabbits exposed to B. burgdorferi infected ticks acquired potent immunity that caused 40-60% mortality of B. burgdorferi infected ticks during the second infestation compared to 15-28% for the uninfected. This might be explained by ELISA data that show that high expression levels of immunogenic proteins in B. burgdorferi infected ticks.

CONCLUSION

Data here suggest that B. burgdorferi infection modified protein content in tick saliva to promote its survival at the tick feeding site. For instance, enzymes; copper/zinc superoxide dismutase that led to production of H2O2 that is toxic to B. burgdorferi were suppressed, while, catalase and thioredoxin that neutralize H2O2, and pyruvate kinase which yields pyruvate that protects Bb from H2O2 killing were enhanced. We conclude data here is an important resource for discovery of effective antigens for a vaccine to prevent LD.

De novo assembly and analysis of the transcriptome of the Dermacentor marginatus genes differentially expressed after blood-feeding and long-term starvation

Background

The ixodid tick Dermacentor marginatus is a vector of many pathogens wide spread in Eurasia. Studies of gene sequence on many tick species have greatly increased the information on tick protective antigen which might have the potential to function as effective vaccine candidates or drug targets for eco-friendly acaricide development. In the current study, RNA-seq was applied to identify D. marginatus sequences and analyze differentially expressed unigenes.

Methods

To obtain a broader picture of gene sequences and changes in expression level, RNA-seq was performed to obtain the whole-body transcriptome data of D. marginatus adult female ticks after engorgement and long-term starvation. Subsequently, the real-time quantitative PCR (RT-qPCR) was applied to validate the RNA-seq data.

Results

RNA-seq produced 30,251 unigenes, of which 32% were annotated. Gene expression was compared among groups that differed by status as newly molted, starved and engorged female adult ticks. Nearly one third of the unigenes in each group were differentially expressed compared to the other two groups, and the most numerous were genes encoding proteins involved in catalytic and binding activities and apoptosis. Selected up-regulated differentially expressed genes in each group were associated to protein, lipids, carbohydrate and chitin metabolism. Blood-feeding and long-term starvation also caused genes differentially expressed in the defense response and antioxidant response. RT-qPCR results indicated 6 differentially expressed transcripts showed similar trends in expression changes with RNA-seq results confirming that the gene expression profiles in transcriptome data is in consistent with RT-qPCR validation.

Conclusions

Obtaining the sequence information of D. marginatus and characterizing the expression pattern of the genes involved in blood-feeding and during starvation would be helpful in understanding molecular physiology of D. marginatus and provides data for anti-tick vaccine and drug development for controlling the tick.

Cloning of four HSPA multigene family members in Haemaphysalis flava ticks

The heat shock protein 70 (HSPA) family and their genes have been studied in ticks and are considered as possible antigen candidates for the development of anti-tick vaccines. However, knowledge about their members, structure and function in ticks is incomplete. Based on our transcriptomic data, the full length of four HSPA genes in Haemaphysalis flava (Acari: Ixodidae) was cloned via rapid amplification of cDNA ends. The open reading frame of HSPA2A, HSPA2B, HSPA5 and HSPA9 was 1920, 1911, 1983 and 2088 bp in length, respectively. Three family signatures and one localization motif were in the encoding proteins. HSPA2A and HSPA2B were predicted to be located at cytoplasm/nucleus, whereas HSPA5 and HSPA9 were at endoplasmic reticulum and mitochondria, respectively. In silico simulation demonstrated that those proteins had distinct numbers of α-helixes, extended strands and coils, and different antigenic epitopes. Expression of HSPA5 and HSPA9 in the salivary gland was significantly higher in partially-engorged female adult ticks than the fully-engorged (P < 0.01) as shown by a quantitative polymerase chain reaction. Our data indicated that H. flava ticks had at least four HSPA genes encoding proteins with different cellular locations, structures and expression profiles, suggesting their diverse roles in tick biology.

A survey of proteins in midgut contents of the tick, Haemaphysalis flava, by proteome and transcriptome analysis

Tick blood meals are stored and digested in their midguts. Blood digestion is complex, and many proteins are involved. Study of the tick-derived proteins in the midgut content may aid in the discovery of active molecules that would be useful for anti-tick vaccines. We analyzed the midgut content proteomes of partially engorged female Haemaphysalis flava, fully engorged female H. flava, and hedgehog serum using liquid chromatography tandem-mass spectrometry and label-free quantitation. In this study, high-confidence protein profiling of tick midgut content was determined. Based on the search against our in-house transcriptome database, the 28 high-confidence proteins were identified. Of these, 17 were identified as tick-derived, and the rest were of unspecified origin (proteins that could not be differentiated as host-derived or tick-derived proteins). The function of these midgut content proteins identified here may involve nutrient transportation, anti-coagulation, erythrocyte lysis, detoxification, lipid metabolism, and immunization. The presence of hemoglobin suggested that the red blood cells were lysed in the gut lumen. The midgut contents contain a large amount of fibrinogen and it has the ability to clot immediately. The midgut contained mostly host-derived proteins, and these host proteins provide rich nutrients for tick development and reproduction. However, some intracellular proteins were also identified, suggesting the possibility of shedding of the midgut epithelium and ingestion of saliva during feeding. This finding advances our understanding of the digestive mechanism and will be useful in the screening of vaccine antigens.

Gene expression during the early stages of host perception and attachment in adult female Rhipicephalus microplus ticks

The cattle tick, Rhipicephalus microplus, is a serious pest of cattle, with significant economic consequences to the livestock industries of tropical and semitropical countries. Rhipicephalus microplus belongs to the Metastriata group of the Ixodidae family known as hard ticks. When adult hard ticks feed, mating has not yet occurred and an initial host attachment phase of 1-2 days is followed by a slow feeding phase that can last several days. Once mating occurs, feeding concludes with a rapid engorgement phase that is completed in 12-36 h. Our group's interest in mining the genome and transcriptome of R. microplus for novel targets for development of tick control technologies led us to investigate the early transcriptional events occurring upon tick attachment and subsequent feeding. We placed newly molted unfed adult R. microplus females upon a bovine host and harvested the attached ticks after 3, 6, 12, and 24 h. We also placed a group of these ticks in a gas-permeable tube taped onto the side of the bovine host. These ticks were able to sense the host but unable to penetrate the tube to begin attachment and were ultimately harvested after 3 h. This study produced a comprehensive transcriptome from newly molted adult ticks and will provide a useful resource for studies of tick feeding and host perception and also assist genome annotation refinements.

De novo transcriptome sequencing and comparative analysis of Haemaphysalis flava Neumann, 1897 at larvae and nymph stages

Haemaphysalis flava Neumann, 1897 (H. flava) is of public health significance due to its capability of transmitting several pathogens such as Rickettsia, Ehrlichia, Bartonella and Francisella tularensis. However, lack of complete genome, transcriptome and proteome information limits the understanding of the biology of H. flava. Here, the total RNA of H. flava was collected separately at larvae and nymph stages and analyzed with high-throughput RNA sequencing technology. The obtained data were assembled and annotated based on the near origin species in the Nr database. The functions of the unigenes were annotated and classified by seven databases, including Nr, Nt, Pfam, KOG, Swiss-Prot, GO and KEGG. A total of 61,850,967 and 79,579,368 clean reads were obtained with a data bulk of 9.28 G and 11.94 G in larvae and nymph stages, respectively. The number of unigenes was 440,896, with 48.6% of them being matched to the Nr database and 51.4% remaining unknown. Additionally, 1,776,404 SNPs were identified in the unigenes. Differential analysis revealed 80 differentially expressed genes (DEGs), including 56 up-regulated genes and 24 down-regulated genes in the nymph versus larvae. qPCR confirmed 4 of the 56 up-regulated genes and 4 of the down-regulated genes. KEGG analysis of the DEGs showed that aldehyde dehydrogenase and sorbitol dehydrogenase, two up-regulated unigenes in nymph versus larvae, were both matched to the top three enriched pathways: "chloroalkane and chloroalkene degradation", "fatty acid degradation" and "glycolysis and gluconeogenesis". This is the first report on the whole transcriptome of H. flava at larvae and nymph stages. This study contributes to the understanding of H. flava at the gene expression level in different developmental stages and provides a theoretical basis for the development of vaccines against H. flava.

Cathepsin L-a novel cysteine protease from Haemaphysalis flava Neumann, 1897

Ixodid ticks are ectoparasites responsible for the transmission of a large number of bacterial, viral, and protozoan pathogens to animals and humans. As long-term blood-pool feeders, the digestion of host blood is critical to their development as well as to the establishment of the sexual cycle of hemoparasites such as Babesia parasites, the agents of human and animal babesiosis. Previous studies have demonstrated that cysteine proteases are involved in blood digestion, embryogenesis, and pathogen transmission in other species of ticks, but their characteristics and functions are still unidentified in Haemaphysalis flava. Here, we describe the characterization of a cysteine protease HfCL from H. flava. We show that HfCL belongs to the L-like papain family of proteases, exhibits high expression in nymphs and adults, and localizes to both the midgut and salivary glands. Biochemical assays using purified recombinant enzyme reveal that rHfCL can hydrolyze the fluorogenic substrate Z-phe-Arg-MCA with optimal activity detected at pH 6. Furthermore, the short-term growth assay indicates that rHfCL can inhibit the intraerythrocytic development of Babesia microti and Babesia gibsoni in vitro.

De novo assembly and analysis of midgut transcriptome of the argasid tick Ornithodoros erraticus and identification of genes differentially expressed after blood feeding

Ticks are hematophagous vectors of great medical and veterinary importance because they transmit numerous pathogenic microorganisms to humans and animals. The argasid Ornithodoros erraticus is the main vector of tick-borne human relapsing fever and African swine fever in the Mediterranean Basin. Tick enterocytes express bioactive molecules that perform key functions in blood digestion, feeding, toxic waste processing and pathogen transmission. To explore new strategies for tick control, in this work we have obtained and compared the midgut transcriptomes of O. erraticus female ticks before and after a blood meal and identified the genes whose expression is differentially regulated after feeding. The transcript sequences were annotated, functionally and structurally characterised and their expression levels compared between both physiological conditions (unfed females and fed females at 2 days post-engorgement). Up to 29,025 transcripts were assembled, and 9290 of them corresponded to differentially expressed genes (DEGs) after feeding. Of these, 4656 genes were upregulated and nearly the same number of genes was downregulated in fed females compared to unfed females. BLASTN and BLASTX analyses of the 29,025 transcripts allowed the annotation of 9072 transcripts/proteins. Among them, the most numerous were those with catalytic and binding activities and those involved in diverse metabolic pathways and cellular processes. The analyses of functional groups of upregulated DEGs potentially related to the digestion of proteins, carbohydrates and lipids, and the genes involved in the defence response and response to oxidative stress, confirm that these processes are narrowly regulated in ticks, highlighting their complexity and importance in tick biology. The expression patterns of six genes throughout the blood digestion period revealed significant differences between these patterns, strongly suggesting that the transcriptome composition is highly dynamic and subjected to important variation along the trophogonic cycle. This may guide future studies aimed at improving the understanding of the molecular physiology of tick digestion and digestion-related processes. The current work provides a more robust and comprehensive understanding of the argasid tick digestive system.

The Mexican bean beetle (Epilachna varivestis) regurgitome and insights into beetle-borne virus specificity

For nearly 400 million years, insects and plants have been embattled in an evolutionary arms race. Insects have developed diverse feeding strategies and behaviors in an effort to circumvent and overcome an extensive collection of plant defense tactics. Sap-sucking insects often inject saliva into hosts plants, which contains a suite of effector proteins and even microbial communities that can alter the plant's defenses. Lacking salivary glands, leaf-feeding beetles represent an interesting group of phytophagous insects. Feeding beetles regurgitate onto leaf surfaces and it is thought that these oral secretions influence insect-plant interactions and even play a role in virus-vector specificity. Since the molecular and biological makeup of the regurgitant is virtually unknown, we carried out RNA sequencing and 16S rDNA analysis on a major soybean pest, Epilachna varivestis, to generate the first ever beetle "regurgitome" and characterize its microbiome. Interestingly, the regurgitant is comprised of a rich molecular assortment of genes encoding putative extracellular proteins involved in digestion, molting, immune defense, and detoxification. By carrying out plant inoculation assays, we reinforced the fundamental role of the regurgitant in beetle-borne virus specificity. Ultimately, these studies begin to characterize the importance of regurgitant in virus transmission and beetle-plant interactions.

Characterization of a glycine-rich protein from Rhipicephalus microplus: tissue expression, gene silencing and immune recognition

Salivary molecules, as glycine-rich proteins (GRPs), are essential to tick attachment and feeding on the host and are suggested to be involved in the host's immune system evasion, therefore representing natural candidates in the search for protective vaccine antigens. This work shows the molecular characterization of a GRP from Rhipicephalus microplus (RmGRP). The cDNA and putative amino acid sequences were analysed, as well as the transcription level in tick tissues/developmental stages, showing the highest levels of gene expression in 1-day-old larvae and salivary glands of fully engorged females. RmGRP gene silencing resulted in a lower hatching rate of larvae from treated females. In addition, recombinant RmGRP (rRmGRP) was recognized by sera from naturally and experimentally infested bovines, displaying considerable differences among the individuals tested. rRmGRP was recognized by anti-saliva and anti-salivary glands sera, while anti-rRmGRP serum recognized RmGRP in saliva and salivary glands, indicating its secretion into the host. The data collected indicate that RmGRP may present roles other than in the tick–host relationship, especially in embryo development. In addition, the high expression in adult females, antigenicity and presence of shared characteristics with other tick protective GRPs turns RmGRP a potential candidate to compose an anti-tick vaccine cocktail.

A proteomic insight into the midgut proteome of Ornithodoros moubata females reveals novel information on blood digestion in argasid ticks

Background

The argasid tick Ornithodoros moubata is the main African vector of the human relapsing fever agent Borrelia duttoni and the African swine fever virus. Together with saliva, the tick midgut forms part of the host-tick-pathogen interface, and numerous midgut proteins play key functions in the blood digestion-related process and the infection and transmission of pathogens. This work explores the composition of the midgut proteome of unfed and fed O. moubata females with the aim to complete the biological information already obtained from the midgut transcriptome and provide a more robust and comprehensive perspective of this biological system.

Methods

Midgut tissues taken from females before feeding and 48 h after feeding were subjected to LC/MS-MS analysis. After functional characterization and classification of the proteins identified, the differences in the proteome between unfed and fed females were analysed and discussed. Additionally, a detailed analysis of particular groups of proteins that are involved in the processes of nutrient digestion and responses to the oxidative stress was carried out.

Results

1491 non-redundant tick proteins were identified: 1132 of them in the midgut of unfed ticks, 1138 in the midgut of fed ticks, and up to 779 shared by both physiological conditions. Overall, the comparative analysis of the midgut proteomes of O. moubata females before and after feeding did not reveal great differences in the number or class of proteins expressed, enzymatic composition or functional classification.

Conclusions

The hemoglobinolytic system in ixodids and argasids is very similar in spite of the fact that they display very different feeding and reproductive strategies. Although the main source of nutrients in ticks are proteins, lipids and carbohydrates also constitute significant nutritional sources and play an important part in the process of blood digestion. The genes and proteins involved in intracellular transport mechanisms, defensive responses, detoxifying responses and stress responses seem to be closely regulated, highlighting the complexity and importance of these processes in tick biology, which in turn assigns them a great interest as targets for therapeutic and immunological interventions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2300-8) contains supplementary material, which is available to authorized users.

New approaches and omics tools for mining of vaccine candidates against vector-borne diseases

Vector-borne diseases (VBDs) present a major threat to human and animal health, as well as place a substantial burden on livestock production. As a way of sustainable VBD control, focus is set on vaccine development. Advances in genomics and other "omics" over the past two decades have given rise to a "third generation" of vaccines based on technologies such as reverse vaccinology, functional genomics, immunomics, structural vaccinology and the systems biology approach. The application of omics approaches is shortening the time required to develop the vaccines and increasing the probability of discovery of potential vaccine candidates. Herein, we review the development of new generation vaccines for VBDs, and discuss technological advancement and overall challenges in the vaccine development pipeline. Special emphasis is placed on the development of anti-tick vaccines that can quell both vectors and pathogens.

Gut transcriptome analysis on females of Ornithodoros mimon (Acari: Argasidae) and phylogenetic inference of ticks

Abstract Ornithodoros mimon is an argasid tick that parasitizes bats, birds and opossums and is also harmful to humans. Knowledge of the transcripts present in the tick gut helps in understanding the role of vital molecules in the digestion process and parasite-host relationship, while also providing information about the evolution of arthropod hematophagy. Thus, the present study aimed to know and ascertain the main molecules expressed in the gut of argasid after their blood meal, through analysis on the gut transcriptome of engorged females of O. mimon using 454-based RNA sequencing. The gut transcriptome analysis reveals several transcripts associated with hemoglobin digestion, such as serine, cysteine, aspartic proteases and metalloenzymes. The phylogenetic analysis on the peptidases confirmed that most of them are clustered with other tick genes. We recorded the presence a cathepsin O peptidase-coding transcript in ticks. The topology of the phylogenetic inferences, based on transcripts of inferred families of homologues, was similar to that of previous reports based on mitochondrial genome and nuclear rRNA sequences. We deposited 2,213 sequence of O. mimon to the public databases. Our findings may help towards better understanding of important argasid metabolic processes, such as digestion, nutrition and immunity.

PCR denaturing gradient gel electrophoresis as a useful method to identify of intestinal bacteria flora in Haemaphysalis flava ticks

Ticks are blood-sucking ectoparasites of great medical and veterinary significance that can transmit bacteria, protozoa, fungi and viruses, causing a variety of human and animal diseases worldwide. In the present study, the intestinal bacterial flora associated with Haemaphysalis flava ticks in different developmental stages were analyzed using polymerase chain reaction (PCR) and denaturing gradient gel electrophoresis (DGGE). Eleven distinct DGGE bands were found using PCR-DGGE method. Sequences analyses indicated that they belonged to Bacillus cereus, Candidatus rickettsia, Erwinia sp., Klebsiella pneumoniae, Pectobacterium carotovorum, Pseudomonas aeruginosa, Rickettsia peacockii, Rickettsia helvetica, Rickettsia slovaca, Staphylococcus simulans and Uncultured bacterium clone. Our find that the K. pneumoniae and P. aeruginosa isolates were presented in all H. flava ticks in different developmental stages. The present results indicated that zoonotic pathogens are present in H. flava ticks in Henan province, China. To our knowledge, this is the first report on intestinal bacterial flora associated with H. flava ticks in China.

Functional annotation and analysis of the Ornithodoros moubata midgut genes differentially expressed after blood feeding

The argasid tick Ornithodoros moubata is the main vector of the African swine fever and the human relapsing fever in Africa. As part of the host-parasite-pathogen interface, the tick midgut expresses key proteins for tick survival and tick-borne pathogen transmission. Accordingly, midgut proteins are potential targets for the development of new drugs and vaccines aimed at tick control, and obtaining proteomic and transcriptomic data from the O. moubata midgut would facilitate the identification of such target candidates. With this aim, we have assembled and characterized the midgut transcriptome of O. moubata females before and 48h after a blood meal, and identified the genes that are differentially expressed in the midgut after feeding. Overall, 23,863 transcripts were obtained, and of them, 9,164 were identified and annotated. The most represented molecular functions were catalytic and binding activities, and the most represented biological processes were metabolic, cellular and single-organism processes. KEGG analysis of the annotated sequences assigned up to 3,053 of them to 130 active pathways, among which, the top 30 pathways were mostly metabolic routes. Differential expression analysis between unfed and fed ticks detected 8,026 Differentially Expressed Genes (DEGs), 4,093 up-regulated and 3,933 down-regulated, respectively. The biological significance of these DEGs was further investigated using the KEEG, Pfam and GO databases. The functional groups of the genes/proteins predicted to be involved in the processes of blood digestion, nutrient transport and metabolism, and in responses related to defence and oxidative stress are discussed in more detail. This work reports the first midgut transcriptome analysis of an argasid tick species, and provides a wealth of novel molecular information about the argasid machinery involved in blood feeding and digestion. This information represents a starting point for the development of alternative strategies for tick control.

Cloning and expression pattern of a heat shock cognate protein 70 gene in ticks (Haemaphysalis flava)

Ticks and tick-borne-diseases have serious public health implications, and screening feasible protein candidates for vaccines development is identified to be an effective alternative to control of tick infestations. In current study, we focused on cloning the full-length gene encoding a heat shock cognate protein 70 (Hsc70), a molecular chaperone of critical functional roles belonging to heat shock protein 70 (HSP70) family, in salivary glands of Haemaphysalis flava, namely Hf-Hsc70, and analyzing the expression of Hf-Hsc70 in different life phases, organs and ambient temperatures. Rapid amplification of cDNA ends (RACE) was performed to amplify the 5' and 3' ends of Hf-Hsc70. The expression profiles of Hf-Hsc70 were studied by semi-quantitative real-time PCR (RT-PCR). The full-length of Hf-Hsc70 was 2363 bp, and contained an ORF of 1965 bp encoding a protein of 648 amino acids. The expression levels of Hf-Hsc70 at different life phases were in the order of female larvae < female fully engorged nymphs < male adult ticks < female full engorged adult ticks < female half engorged adult ticks. The relative expression of Hf-Hsc70 in salivary glands was steadily higher than that in midguts (p < 0.05) regardless of feeding status. A 3-h of heat stress did not significantly induce the up-regulation of Hf-Hsc70 transcription. These results indicated that Hf-Hsc70 was a constitutive form of HSP70 family, and its expression pattern in different life phases and organs suggested a possible role in blood feeding, which would further make Hsc70 a potential candidate for the development of vaccines against ticks.

De novo transcriptome sequencing and analysis of the juvenile and adult stages of Fasciola gigantica

Fasciola gigantica is regarded as the major liver fluke causing fasciolosis in livestock in tropical countries. Despite the significant economic and public health impacts of F. gigantica there are few studies on the pathogenesis of this parasite and our understanding is further limited by the lack of genome and transcriptome information. In this study, de novo Illumina RNA sequencing (RNA-seq) was performed to obtain a comprehensive transcriptome profile of the juvenile (42days post infection) and adult stages of F. gigantica. A total of 49,720 unigenes were produced from juvenile and adult stages of F. gigantica, with an average length of 1286 nucleotides (nt) and N50 of 2076nt. A total of 27,862 (56.03%) unigenes were annotated by BLAST similarity searches against the NCBI non-redundant protein database. Because F. gigantica needs to feed and/or digest host tissues, some proteases (including cysteine proteases and aspartic proteases), which play a role in the degradation of host tissues (protein), have been paid more attention in the present study. A total of 6511 distinct genes were found differentially expressed between juveniles and adults, of which 3993 genes were up-regulated and 2518 genes were down-regulated in adults versus juveniles, respectively. Moreover, stage-specific differentially expressed genes were identified in juvenile (17,009) and adult (6517) F. gigantica. The significantly divergent pathways of differentially expressed genes included cAMP signaling pathway (226; 4.12%), proteoglycans in cancer (256; 4.67%) and focal adhesion (199; 3.63%). The transcription pattern also revealed two egg-laying-associated pathways: cGMP-PKG signaling pathway and TGF-β signaling pathway. This study provides the first comparative transcriptomic data concerning juvenile and adult stages of F. gigantica that will be of great value for future research efforts into understanding parasite pathogenesis and developing vaccines against this important parasite.

RNA-seq analyses of the midgut from blood- and serum-fed Ixodes ricinus ticks

Adult females of the genus Ixodes imbibe blood meals exceeding about 100 times their own weight within 7‒9 days. During this period, ticks internalise components of host blood by endocytic digest cells that line the tick midgut epithelium. Using RNA-seq, we aimed to characterise the midgut transcriptome composition in adult Ixodes ricinus females during early and late phase of engorgement. To address specific adaptations to the haemoglobin-rich diet, we compared the midgut transcriptomes of genetically homogenous female siblings fed either bovine blood or haemoglobin-depleted serum. We noted that tick gut transcriptomes are subject to substantial temporal-dependent expression changes between day 3 and day 8 of feeding. In contrast, the number of transcripts significantly affected by the presence or absence of host red blood cells was low. Transcripts relevant to the processes associated with blood-meal digestion were analysed and involvement of selected encoded proteins in the tick midgut physiology discussed. A total of 7215 novel sequences from I. ricinus were deposited in public databases as an additional outcome of this study. Our results broaden the current knowledge of tick digestive system and may lead to the discovery of potential molecular targets for efficient tick control.