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

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.

The diurnal salivary glands transcriptome of Dermacentor nuttalli from the first four days of blood feeding

The ixodid tick Dermacentor nuttalli is distributed from southern Siberia to North China and is a vector of many pathogens. This species can have severe impacts on animal husbandry and human health. To date, the control of D. nuttalli is limited to the use of acaricides such as organophosphorus, synthetic pyrethroids and amidine pesticides. There are no environmentally friendly or reliable prevention and control measures, and little is known regarding key antigens involved in blood feeding. Salivary glands are major tissues involved in the blood feeding and pathogen transmission of ticks. Therefore, this study focused on salivary glands tissue to identify the dominant antigens of D. nuttalli involved in tick feeding. For this, high-throughput RNA sequencing (RNA-seq) was used for analysis. The transcriptome of female D. nuttalli ticks was assembled and characterized, and differentially expressed genes (DEGs) were identified in the salivary glands of ticks that had not fed (0 h) and of ticks after 24, 48, 72 and 96 h of feeding. There were 22,802,784, 22,275,013, 26,629,453, 24,982,389, and 22,596,230 high-quality clean reads obtained from salivary glands tissues at the five different blood feeding time points. The total number of annotated unigenes was 100,347. The differences in gene expression between different time points were compared, and functional enrichment was performed. Quantitative reverse transcription PCR (RT‒qPCR) was used to validate the RNA-seq results, the results of which showed that the differences in expressed transcripts presented similar trends. Among the identified DEGs, the most numerous were those with catalytic and binding activities and those involved in diverse metabolic pathways and cellular processes. The expression patterns of homologous and family-member proteins throughout the blood feeding period exhibited significant differences, strongly suggesting that the transcriptome composition is highly dynamic and likely subjected to important variation throughout the life cycle. Studies of gene sequences in D. nuttalli will greatly increase the information on tick protective antigens, which could potentially function as effective vaccine candidates or drug targets for the development of environmentally friendly acaricides.

De novo assembled transcriptomics assisted label-free quantitative proteomics analysis reveals sex-specific proteins in the intestinal tissue of Haemaphysalis qinghaiensis

The hard tick Haemaphysalis qinghaiensis is the vector of a wide variety of infectious agents, such as spirochetes and other bacteria as well as viruses in the western plateau of China. Tick midgut is the key tissue involved in the host-pathogen-vector interface. Multiple midgut proteins are related to key functions in blood digestion, tick survival, and tick-borne pathogen transmission. However, information on the sex-specific proteins expressed in the midgut tissue of H. qinghaiensis for which the genome has not been sequenced is limited. Hence, we assembled and characterized the transcriptome of the H. qinghaiensis midgut and identified the differentially expressed genes (DEGs) in female and male ticks. The sequencing of the mRNA for this nonmodel species is essential for producing a protein database for mass spectrometry-based identification. Here, we combined high-throughput parallel sequencing and label-free quantitative proteomics analysis to extensively characterize the tick midgut using massive RNA sequencing and mass spectrometry, which allowed the detection of genes and proteins. A total of 279,186 transcripts were annotated into 125,790 coding sequences (CDSs), which were manually curated into 96 different gene families. A total of 12,837 DEGs between the two sexes were found by RNA-seq analysis. Of these, 5401 were upregulated genes, while 7436 were downregulated genes. The most common molecular functions were those related to the endocrine system, translation, signal transduction, transport, and catabolism. Meanwhile, the most common biological processes were related to cellular processes, metabolic processes, cellular anatomical entities, and cargo receptor activities. An analysis of the label-free protein quantitation dataset showed 272 upregulated proteins and 46 downregulated proteins when the fold-change was >2.0 (LC-MS/MS). Association analysis of the transcriptome and proteome with GO functional enrichment showed that the majority of the genes (proteins) were those related to catalytic activity, binding, cellular processes, metabolic processes, and responses to stimuli. This study aims to elucidate the digestive physiology of H. qinghaiensis as well as its physiological sexual dimorphism. This will allow the identification of protein candidates with physiological importance that could be used as targets to control the vector as well as the transmission of tick-borne pathogens to humans and animals.

Comparative analysis of the anticoagulant activities and immunogenicity of HSC70 and HSC70TKD of Haemaphysalis flava

BACKGROUND

Haemaphysalis flava is a hematophagous ectoparasite that acquires the nutrition needed for development and reproduction by sucking blood and digesting the blood meal. During blood-sucking and blood-meal digestion, the prevention of blood coagulation is important for this tick. Previous studies have shown that heat shock cognate 70 (HSC70) protein has certain anticoagulant activities, but its immunogenicity remains unclear. Also, whether the mutation of individual bases of the TKD-like peptide of HSC70 through the overlap extension method can change its anticoagulant activities and immunogenicity remains to be investigated.

METHODS

The gene encoding the HSC70 protein was cloned from a complementary DNA library synthesized from H. flava. The coding gene of the TKD-like peptide of HSC70 was mutated into a TKD peptide coding gene (HSC70^TKD) using the overlap extension method. Escherichia coli prokaryotic expression plasmids were constructed to obtain the recombinant proteins of HSC70 (rHSC70) and HSC70^TKD (rHSC70^TKD). The purified rHSC70 and rHSC70^TKD were evaluated at different concentrations for anticoagulant activities using four in vitro clotting assays. Emulsifying recombinant proteins with complete and incomplete Freund's adjuvants were subcutaneously immunized in Sprague Dawley rats. The serum antibody titers and serum concentrations of interferon-gamma (IFN-γ) and interleukin-4 (IL-4) were detected using an indirect enzyme-linked immunosorbent assay to assess the immunogenicity of rHSC70 and rHSC70^TKD.

RESULTS

The open reading frame of HSC70 was successfully amplified and found to have a length of 1958 bp. The gene encoding the TKD-like peptide of HSC70 was artificially mutated, with the 1373-position adenine (A) of the original sequence mutated into guanine (G), the 1385-position cytosine (C) mutated into G and the 1386-position G mutated into C. rHSC70 and rHSC70^TKD that fused with His-tag were obtained using the expression plasmids pET-28a-HSC70 and pET-28a-HSC70^TKD, respectively. rHSC70 and rHSC70^TKD prolonged the thrombin time (TT) and reduced the fibrinogen (FIB) content in the plasma, but did not affect the prothrombin time (PT) or activated partial thromboplastin time (APTT) when compared to the negative control. Interestingly, the ability of rHSC70^TKD to prolong the TT and reduce the FIB content in the plasma was better than that of rHSC70. The specific antibody titers of both rHSC70 and rHSC70^TKD in rat serum reached 1:124,000 14 days after the third immunization. The serum concentration of IFN-γ in the rHSC70^TKD group was higher than that in the rHSC70 group. The rHSC70 group has the highest serum concentration of IL-4, and the serum concentration of IL-4 in the rHSC70^TKD group was higher than that in the negative group.

CONCLUSIONS

rHSC70 and rHSC70^TKD exhibited anticoagulant activities by prolonging the TT and reducing the FIB content in vitro. rHSC70^TKD had better anticoagulant activities than rHSC70. Both rHSC70 and rHSC70^TKD had good immunogenicity and induced humoral and cellular immunity.

Contemporary diagnostics for medically relevant fastidious microorganisms belonging to the genera Anaplasma, Bartonella, Coxiella, Orientia, and Rickettsia

Many of the human infectious pathogens—especially the zoonotic or vector-borne bacteria—are fastidious organisms that are difficult to cultivate because of their strong adaption to the infected host culminating in their near-complete physiological dependence on this environment. These bacterial species exhibit reduced multiplication rates once they are removed from their optimal ecological niche. This fact complicates the laboratory diagnosis of the disease and hinders the detection and further characterization of the underlying organisms, e.g. at the level of their resistance to antibiotics due to their slow growth. Here, we describe the current state of microbiological diagnostics for five genera of human pathogens with a fastidious laboratory lifestyle. For Anaplasma spp., Bartonella spp., Coxiella burnetii, Orientia spp. and Rickettsia spp., we will summarize the existing diagnostic protocols, the specific limitations for implementation of novel diagnostic approaches and the need for further optimization or expansion of the diagnostic armamentarium. We will reflect upon the diagnostic opportunities provided by new technologies including mass spectrometry and next-generation nucleic acid sequencing. Finally, we will review the (im)possibilities of rapidly developing new in vitro diagnostic tools for diseases of which the causative agents are fastidiously growing and therefore hard to detect.

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.

Quantitative Visions of Reality at the Tick-Host Interface: Biochemistry, Genomics, Proteomics, and Transcriptomics as Measures of Complete Inventories of the Tick Sialoverse

Species have definitive genomes. Even so, the transcriptional and translational products of the genome are dynamic and subject to change over time. This is especially true for the proteins secreted by ticks at the tick-host feeding interface that represent a complex system known as the sialoverse. The sialoverse represent all of the proteins derived from tick salivary glands for all tick species that may be involved in tick-host interaction and the modulation of the host's defense mechanisms. The current study contemplates the advances made over time to understand and describe the complexity present in the sialoverse. Technological advances at given periods in time allowed detection of functions, genes, and proteins enabling a deeper insight into the complexity of the sialoverse and a concomitant expansion in complexity with as yet, no end in sight. The importance of systematic classification of the sialoverse is highlighted with the realization that our coverage of transcriptome and proteome space remains incomplete, but that complete descriptions may be possible in the future. Even so, analysis and integration of the sialoverse into a comprehensive understanding of tick-host interactions may require further technological advances given the high level of expected complexity that remains to be uncovered.