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

Characterization and expression analysis of heat shock proteins HfHspc1 and HfHspc4 in Haemaphysalis flava ticks

In this study, the cDNAs of heat shock cognate 90 protein 1 (HfHspc1) and heat shock cognate 90 protein 4 (HfHspc4) from Haemaphysalis flava (Acari: Ixodidae) were obtained using the rapid amplification of cDNA ends (RACE) approach, and the expression patterns of HfHspc1 and HfHspc4 in different developmental stages, engorgement stages and tick organs were analyzed by qPCR. The full length of HfHspc1 was 2411 bp, and its open reading frame (ORF) was 2196 bp, encoding a protein of 732 aa, containing five HSPC family signatures, with MEEVD motif at its extreme C-terminal. The full length of HfHspc4 was 2800 bp, and its ORF was 2364 bp, encoding a protein of 789 aa, containing a signal peptide and five family signatures, with HEEL motif at its extreme C-terminal. The expression of HfHspc1 and HfHspc4 was the highest in males, while it was significantly the highest in the ovaries of fully engorged females, potentially implying the roles of these proteins in the successful digestion of blood and development of eggs.

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.

Egg protein profile and dynamics during embryogenesis in Haemaphysalis flava ticks

Tick eggs contain all essential proteins for embryogenesis, and egg proteins are a potential reservoir of tick-protective antigens. However, the protein profile and dynamics during embryonic development remain unknown. This study aimed to depict the protein profile and dynamics in tick embryogenesis, further providing protein candidates for targeted interventions. Eggs from Haemaphysalis flava ticks were incubated at 28 °C and 85% relative humidity. On days 0 (newly laid eggs without incubation), 7, 14 and 21, eggs were collected, dewaxed and subject to protein extraction. Extracted proteins were digested by filter-aided sample preparation and analyzed by liquid chromatography-tandem mass spectrometry (LC/MS-MS). MS data were searched against an in-house H. flava protein database for tick-derived protein identification. Abundances of 40 selected high-confidence proteins were further quantified by LC-parallel reaction monitoring (PRM)/MS analysis throughout egg incubation. A total of 93 high-confidence proteins were identified in eggs on 0-day incubation. Identified proteins belonged to seven functional categories: transporters, enzymes, proteinase inhibitors, immunity-related proteins, cytoskeletal proteins, heat shock proteins and uncharacterized proteins. The enzyme category contained the most types of proteins. Neutrophil elastase inhibitors represented the most abundant proteins in terms of intensity-based absolute-protein-quantification. LC-PRM/MS revealed that the abundances of 20 proteins increased including enolase, calreticulin, actin, GAPDH et cetera, and the abundances of 11 proteins decreased including vitellogenins, neutrophil elastase inhibitor, carboxypeptidase Q, et cetera from 0- to 21-day incubation. This study provides the most comprehensive egg protein profile and dynamics during tick embryogenesis. Further investigations are needed to test the tick-control efficacy by targeting the egg proteins.

Cloning, expression, and function of ferritins in the tick Haemaphysalis flava

The full-length cDNA of two ferritins of Haemaphysalis flava were cloned after which recombinant Hf-FER1 and Hf-FER2 were expressed and their function was analyzed. In addition, RNA interference (RNAi) based on the injection of Hf-fer1 or Hf-fer2 dsRNA into fully engorged female ticks was performed. The cDNA encoding Hf-FER1 is 834 bp in length. It contains an iron-responsive element in the 5' untranslated region and encodes 174 amino acid residues. The full-length cDNA of Hf-FER2 contains 696 bp and encodes 199 amino acids, including a putative signal peptide sequence. Hf-FER1 and Hf-FER2 both have the ferroxidase iron center and the ferrihydrite nucleation center. The evolutionary relationship of Hf-FER1 and Hf-FER2 was established, and the predicted quaternary structures were assembled as typical spherical shells composed of 24 subunits which was demonstrated by nature PAGE. Real-time PCR showed that Hf-fer1 and Hf-fer2 were expressed in all developmental stages, with the highest expression in fully engorged females. The expression of Hf-fer1 and Hf-fer2 were relatively high in unfed larvae. Hf-fer1 was expressed in all tissues and was especially abundant in the salivary glands of fully engorged females. In contrast, the highest levels of Hf-fer2 were found in the midgut of fully engorged females, and no expression was found in the salivary glands of this life stage. Both recombinant Hf-FER1 and Hf-FER2 had iron-binding capabilities. Silencing of both Hf-fer1 and Hf-fer2 affected fecundity. Compared to the control, the percentage of ticks that laid eggs in the Hf-fer1 and Hf-fer2 RNAi groups was 73.3% and 66.7%, respectively. The silenced ticks that laid eggs had lower egg weight to body weight ratios, and the eggs had abnormal morphologies. The hatchability of eggs with normal morphology in the Hf-fer1 and Hf-fer2 silenced groups was 47.8% and 22.8%, respectively, which was significantly different from the control group (P < 0.005). These findings indicate that Hf-FER1 and Hf-FER2 play important roles in the iron storage of H. flava.