Neuropeptide signaling sequences identified by pyrosequencing of the American dog tick synganglion transcriptome during blood feeding and reproduction

Developmental Regulation of Corazonin, Eclosion Hormone, and Bursicon Messages and RNAi Suppression of Corazonin in Adult, Female American Dog Ticks, Dermacentor variabilis

The insect molting process is critical to growth and development and is regulated in part by the neuropeptides corazonin, eclosion hormone, and α and β bursicon. We found messages in a synganglion transcriptome from adult, female American dog ticks, Dermacentor variabilis (that do not molt), with a high similarity to the larval insect neuropeptides that control molting. The phylogenetic analysis of the tick putative neuropeptides compared to other arthropods is discussed in detail. The relative gene expression of these peptides was determined by quantitative PCR during the following adult developmental stages: (i) virgin, unfed 0-24 h after entering the adult stage (non-host-seeking), (ii) host-seeking, unfed, and not mated (3 d after emergence), (iii) part-fed (unmated, attached to host; 1st and 3rd day after emergence), (iv) mated (females are part-fed; allowed to mate for ≤1 day, 7th day after emergence), (v) mated repletes (completion of blood feeding but still attached to host), and (vi) post-drop-off (from host) with egg laying starting within 1 d of detachment. Eclosion hormone transcript levels peaked at mating and at drop-off. Bursicon α levels were highest just after molting into adults, with a second smaller peak in replete females. Bursicon β levels were highest (32-fold) post-drop-off. Corazonin message levels peaked in part-feds and were much higher (40-fold) in repletes compared to 0-24 h after emergence. RNAi suppression of the corazonin message by injection in newly molted ticks reduced oviposition and the number of vitellogenic eggs in the ovaries at drop-off but had no apparent effect on host-seeking, partial feeding, mating, feeding to repletion, and drop-off. The possible roles of these transcripts in adult, female tick development are discussed.

Behavioral Manipulation of Ixodes scapularis by Ehrlichia muris eauclairensis: Implications for Tick-Borne Disease Transmission

Tick-borne diseases pose significant risks to both animals and humans, with emerging pathogens like Ehrlichia muris eauclairensis (EME) underscoring the need for a deeper understanding of pathogen-vector interactions and tick fitness. This study investigates the impact of EME on Ixodes scapularis nymphs, revealing significant behavioral changes in EME-positive ticks. These ticks exhibited increased movement speed, faster bite site-seeking for attachment, and prolonged feeding durations compared to control ticks. Proteomic analyses of the tick synganglion during resting and feeding phases identified 196 differentially expressed proteins in EME-positive ticks, including multiple proteins associated with nicotinic acetylcholine signaling pathways. Our findings indicated altered neuropeptide expression related to stimulus response and activity, suggesting changes in neurophysiology. This research provides the first evidence of behavioral manipulation by an Ehrlichia species, indicating that the tick nervous system is a site of bacterial influence and a potential target for interventions. These findings offer new insights into pathogen-vector dynamics that could lead to the development of transmission-blocking therapies, significantly impacting tick fitness and disease transmission.

Genome sequences of four Ixodes species expands understanding of tick evolution

BACKGROUND

Ticks, hematophagous Acari, pose a significant threat by transmitting various pathogens to their vertebrate hosts during feeding. Despite advances in tick genomics, high-quality genomes were lacking until recently, particularly in the genus Ixodes, which includes the main vectors of Lyme disease.

RESULTS

Here, we present the genome sequences of four tick species, derived from a single female individual, with a particular focus on the European species Ixodes ricinus, achieving a chromosome-level assembly. Additionally, draft assemblies were generated for the three other Ixodes species, I. persulcatus, I. pacificus, and I. hexagonus. The quality of the four genomes and extensive annotation of several important gene families have allowed us to study the evolution of gene repertoires at the level of the genus Ixodes and of the tick group. We have determined gene families that have undergone major amplifications during the evolution of ticks, while an expression atlas obtained for I. ricinus reveals striking patterns of specialization both between and within gene families. Notably, several gene family amplifications are associated with a proliferation of single-exon genes-most strikingly for fatty acid elongases and sulfotransferases.

CONCLUSIONS

The integration of our data with existing genomes establishes a solid framework for the study of gene evolution, improving our understanding of tick biology. In addition, our work lays the foundations for applied research and innovative control targeting these organisms.

Identification and characterization of ecdysis-related neuropeptides in the lone star tick Amblyomma americanum

Introduction

The lone star tick, Amblyomma americanum, is an important ectoparasite known for transmitting diseases to humans and animals. Ecdysis-related neuropeptides (ERNs) control behaviors crucial for arthropods to shed exoskeletons. However, ERN identification and characterization in A. americanum remain incomplete.

Methods

We investigated ERNs in A. americanum, assessing their evolutionary relationships, protein properties, and functions. Phylogeny, sequence alignment, and domain structures of ERNs were analyzed. ERN functionality was explored using enrichment analysis, and developmental and tissue-specific ERN expression profiles were examined using qPCR and RNAi experiments.

Results and discussion

The study shows that ERN catalogs (i.e., eclosion hormone, corazonin, and bursicon) are found in most arachnids, and these ERNs in A. americanum have high evolutionary relatedness with other tick species. Protein modeling analysis indicates that ERNs primarily consist of secondary structures and protein stabilizing forces (i.e., hydrophobic clusters, hydrogen bond networks, and salt bridges). Gene functional analysis shows that ENRs are involved in many ecdysis-related functions, including ecdysis-triggering hormone activity, neuropeptide signaling pathway, and corazonin receptor binding. Bursicon proteins have functions in chitin binding and G protein-coupled receptor activity and strong interactions with leucine-rich repeat-containing G-protein coupled receptor 5. ERNs were expressed in higher levels in newly molted adults and synganglia. RNAi-mediated knockdown of burs α and burs β expression led to a significant decrease in the expression of an antimicrobial peptide, defensin, suggesting they might act in signaling or regulatory pathways that control the expression of immune-related genes. Arthropods are vulnerable immediately after molting because new cuticles are soft and susceptible to injury and pathogen infections. Bursicon homodimers act in prophylactic immunity during this vulnerable period by increasing the synthesis of transcripts encoding antimicrobial peptides to protect them from microbial invasion. Collectively, the expression pattern and characterization of ERNs in this study contribute to a deeper understanding of the physiological processes in A. americanum.

Putative target sites in synganglion for novel ixodid tick control strategies

Acaricide resistance is a global problem that has impacts worldwide. Tick populations with broad resistance to all commercially available acaricides have been reported. Since resistance selection in ticks and their role in pathogen transmission to animals and humans result in important economic and public health burden, it is essential to develop new strategies for their control (i.e., novel chemical compounds, vaccines, biological control). The synganglion is the tick central nervous system and it is responsible for synthesizing and releasing signaling molecules with different physiological functions. Synganglion proteins are the targets of the majority of available acaricides. In this review we provide an overview of the mode-of-action and resistance mechanisms against neurotoxic acaricides in ticks, as well as putative target sites in synganglion, as a supporting tool to identify new target proteins and to develop new strategies for tick control.

Mating, Sperm Transfer and Oviposition in Soft Ticks (Acari: Argasidae), a Review

This review addresses the physiology and behavioral events involved in the reproduction of soft ticks (family Argasidae), with special attention to the events of their adult life: mating, sperm transfer and egg-laying. Many of these aspects are held in common with hard ticks, but the repeated short duration of feeding bouts in soft ticks, in contrast to the extended single engorgements of hard ticks, has consequences peculiar to soft tick reproduction. Reviewed are the dramatic external mechanism of sperm transfer, the unusual maturation and unique morphology and motility of the spermatozoa, the mechanism of oogenesis and its hormonal control, the mystery of fertilization, the involvement of pheromones in mating, the control of reproductive arrests and the vertical transmission of symbiotes in reproduction. Jumping-off points for further investigation are discussed throughout.

A draft of the genome of the Gulf Coast tick, Amblyomma maculatum

The Gulf Coast tick, Amblyomma maculatum, inhabits the Southeastern states of the USA bordering the Gulf of Mexico, Mexico, and other Central and South American countries. More recently, its U.S. range has extended West to Arizona and Northeast to New York state and Connecticut. It is a vector of Rickettsia parkeri and Hepatozoon americanum. This tick species has become a model to study tick/Rickettsia interactions. To increase our knowledge of the basic biology of A. maculatum we report here a draft genome of this tick and an extensive functional classification of its proteome. The DNA from a single male tick was used as a genomic source, and a 10X genomics protocol determined 28,460 scaffolds having equal or more than 10 Kb, totaling 1.98 Gb. The N50 scaffold size was 19,849 Kb. The BRAKER pipeline was used to find the protein-coding gene boundaries on the assembled A. maculatum genome, discovering 237,921 CDS. After trimming and classifying the transposable elements, bacterial contaminants, and truncated genes, a set of 25,702 were annotated and classified as the core gene products. A BUSCO analysis revealed 83.4% complete BUSCOs. A hyperlinked spreadsheet is provided, allowing browsing of the individual gene products and their matches to several databases.

Pyrokinin receptor silencing in females of the southern cattle tick Rhipicephalus (Boophilus) microplus is associated with a reproductive fitness cost

BACKGROUND

Rhipicephalus microplus is the vector of deadly cattle pathogens, especially Babesia spp., for which a recombinant vaccine is not available. Therefore, disease control depends on tick vector control. However, R. microplus populations worldwide have developed resistance to available acaricides, prompting the search for novel acaricide targets. G protein-coupled receptors (GPCRs) are involved in the regulation of many physiological processes and have been suggested as druggable targets for the control of arthropod vectors. Arthropod-specific signaling systems of small neuropeptides are being investigated for this purpose. The pyrokinin receptor (PKR) is a GPCR previously characterized in ticks. Myotropic activity of pyrokinins in feeding-related tissues of Rhipicephalus sanguineus and Ixodes scapularis was recently reported.

METHODS

The R. microplus pyrokinin receptor (Rhimi-PKR) was silenced through RNA interference (RNAi) in female ticks. To optimize RNAi, a dual-luciferase assay was applied to determine the silencing efficiency of two Rhimi-PKR double-stranded RNAs (dsRNA) prior to injecting dsRNA in ticks to be placed on cattle. Phenotypic variables of female ticks obtained at the endpoint of the RNAi experiment were compared to those of control female ticks (non-injected and beta-lactamase dsRNA-injected). Rhimi-PKR silencing was verified by quantitative reverse-transcriptase PCR in whole females and dissected tissues.

RESULTS

The Rhimi-PKR transcript was expressed in all developmental stages. Rhimi-PKR silencing was confirmed in whole ticks 4 days after injection, and in the tick carcass, ovary and synganglion 6 days after injection. Rhimi-PKR silencing was associated with an increased mortality and decreased weight of both surviving females and egg masses (P < 0.05). Delays in repletion, pre-oviposition and incubation periods were observed (P < 0.05).

CONCLUSIONS

Rhimi-PKR silencing negatively affected female reproductive fitness. The PKR appears to be directly or indirectly associated with the regulation of female feeding and/or reproductive output in R. microplus. Antagonists of the pyrokinin signaling system could be explored for tick control.

Transcriptome of the synganglion in the tick Ixodes ricinus and evolution of the cys-loop ligand-gated ion channel family in ticks

BACKGROUND

Ticks represent a major health issue for humans and domesticated animals. Exploring the expression landscape of the tick's central nervous system (CNS), known as the synganglion, would be an important step in understanding tick physiology and in managing tick-borne diseases, but studies on that topic are still relatively scarce. Neuron-specific genes like the cys-loop ligand-gated ion channels (cys-loop LGICs, or cysLGICs) are important pharmacological targets of acaricides. To date their sequence have not been well catalogued for ticks, and their phylogeny has not been fully studied.

RESULTS

We carried out the sequencing of transcriptomes of the I. ricinus synganglion, for adult ticks in different conditions (unfed males, unfed females, and partially-fed females). The de novo assembly of these transcriptomes allowed us to obtain a large collection of cys-loop LGICs sequences. A reference meta-transcriptome based on synganglion and whole body transcriptomes was then produced, showing high completeness and allowing differential expression analyses between synganglion and whole body. Many of the genes upregulated in the synganglion were associated with neurotransmission and/or localized in neurons or the synaptic membrane. As the first step of a functional study of cysLGICs, we cloned the predicted sequence of the resistance to dieldrin (RDL) subunit homolog, and functionally reconstituted the first GABA-gated receptor of Ixodes ricinus. A phylogenetic study was performed for the nicotinic acetylcholine receptors (nAChRs) and other cys-loop LGICs respectively, revealing tick-specific expansions of some types of receptors (especially for Histamine-like subunits and GluCls).

CONCLUSIONS

We established a large catalogue of genes preferentially expressed in the tick CNS, including the cysLGICs. We discovered tick-specific gene family expansion of some types of cysLGIC receptors, and a case of intragenic duplication, suggesting a complex pattern of gene expression among different copies or different alternative transcripts of tick neuro-receptors.

Neuropeptides in Rhipicephalus microplus and other hard ticks

The synganglion is the central nervous system of ticks and, as such, controls tick physiology. It does so through the production and release of signaling molecules, many of which are neuropeptides. These peptides can function as neurotransmitters, neuromodulators and/or neurohormones, although in most cases their functions remain to be established. We identified and performed in silico characterization of neuropeptides present in different life stages and organs of Rhipicephalus microplus, generating transcriptomes from ovary, salivary glands, fat body, midgut and embryo. Annotation of synganglion transcripts led to the identification of 32 functional categories of proteins, of which the most abundant were: secreted, energetic metabolism and oxidant metabolism/detoxification. Neuropeptide precursors are among the sequences over-represented in R. microplus synganglion, with at least 5-fold higher transcription compared with other stages/organs. A total of 52 neuropeptide precursors were identified: ACP, achatin, allatostatins A, CC and CCC, allatotropin, bursicon A/B, calcitonin A and B, CCAP, CCHamide, CCRFamide, CCH/ITP, corazonin, DH31, DH44, eclosion hormone, EFLamide, EFLGGPamide, elevenin, ETH, FMRFamide myosuppressin-like, glycoprotein A2/B5, gonadulin, IGF, inotocin, insulin-like peptides, iPTH, leucokinin, myoinhibitory peptide, NPF 1 and 2, orcokinin, proctolin, pyrokinin/periviscerokinin, relaxin, RYamide, SIFamide, sNPF, sulfakinin, tachykinin and trissin. Several of these neuropeptides have not been previously reported in ticks, as the presence of ETH that was first clearly identified in Parasitiformes, which include ticks and mites. Prediction of the mature neuropeptides from precursor sequences was performed using available information about these peptides from other species, conserved domains and motifs. Almost all neuropeptides identified are also present in other tick species. Characterizing the role of neuropeptides and their respective receptors in tick physiology can aid the evaluation of their potential as drug targets.

Prevention of tick-borne diseases: challenge to recent medicine

Abstract

Ticks represent important vectors and reservoirs of pathogens, causing a number of diseases in humans and animals, and significant damage to livestock every year. Modern research into protection against ticks and tick-borne diseases focuses mainly on the feeding stage, i.e. the period when ticks take their blood meal from their hosts during which pathogens are transmitted. Physiological functions in ticks, such as food intake, saliva production, reproduction, development, and others are under control of neuropeptides and peptide hormones which may be involved in pathogen transmission that cause Lyme borreliosis or tick-borne encephalitis. According to current knowledge, ticks are not reservoirs or vectors for the spread of COVID-19 disease. The search for new vaccination methods to protect against ticks and their transmissible pathogens is a challenge for current science in view of global changes, including the increasing migration of the human population.

Highlights

• Tick-borne diseases have an increasing incidence due to climate change and increased human migration

• To date, there is no evidence of transmission of coronavirus COVID-19 by tick as a vector

• To date, there are only a few modern, effective, and actively- used vaccines against ticks or tick-borne diseases

• Neuropeptides and their receptors expressed in ticks may be potentially used for vaccine design

Functional characterization of the insulin signaling pathway in the hard tick Ixodes ricinus

Ticks are blood-feeding arachnids transmitting a variety of pathogens to humans and animals. A unique trait in tick physiology is their ability to engorge and digest large amounts of host blood, ensuring their high reproductive potential. Activation of the blood digestive machinery in the tick gut, as well as processes controlling maturation of ovaries, are triggered upon blood meal uptake by still largely unknown mechanisms. Sensing of the nutritional status in metazoan organisms is facilitated by the evolutionarily conserved Insulin Signaling Pathway (ISP) and the interlinked Target of Rapamycin (TOR) pathway. Recently, we have identified three components of these pathways in the hard tick Ixodes ricinus midgut transcriptome, namely a putative insulin receptor (InR), and the downstream intracellular serine/threonine kinases AKT and TOR. In this study, we primarily focus on the molecular and functional characterization of the I. ricinus insulin receptor (IrInR), the first InR characterized in Chelicerates. A phylogenetic analysis across the major Arthropod lineages demonstrated that ticks possess only one gene encoding an InR-related molecule. Tissue expression profiling by quantitative PCR in semi-engorged I. ricinus females revealed that the IrInR, as well as AKT (IrAKT) and TOR (IrTOR) are expressed in various organs, with the highest expression being detected in ovaries. We have further evaluated the impact of RNAi-mediated knock-down (KD) of IrInR, IrAKT, and IrTOR on tick blood-feeding and reproductive capacity. Weights of engorged IrInR KD females and laid egg clutches were reduced compared to the control group, and these quantitative parameters clearly correlated with the efficiency of RNAi-KD achieved in individual ticks. The most striking phenotype was observed for IrAKT KD that impaired tick feeding and completely aborted egg production. A recombinant extracellular fragment of the IrInR α-subunit was used to produce antibodies in experimental rabbits to assess its potential as a protective antigen against tick feeding and reproduction. Our data clearly indicate the functionality of the ISP in ticks and demonstrate the need for further investigation of specific roles played by the endogenous insulin-like peptides in tick physiological processes.

Dynamics of Insulin Signaling in the Black-Legged Tick, Ixodes scapularis

Insulin-like peptides (ILPs) have been identified in several invertebrates, particularly insects, and work on these ILPs has revealed many roles including regulation of energy homeostasis, growth, development, and lifespan to name a few. However, information on arthropod ILPs outside of insects is sparse. Studies of Ixodid tick ILPs are particularly scarce, despite their importance as vectors of infectious agents, most notably Lyme disease. The recent publication of the genome of the black-legged tick, Ixodes scapularis, has advanced opportunities to study this organism from a molecular standpoint, a resource sorely needed for an organism with challenging life history requirements for study in the laboratory, such as a long life cycle and obligate, prolonged, blood-feeding at each life stage. Through bioinformatics searches of the tick genome and other available I. scapularis databases, we identified four putative ILP sequences. Full-length sequences of these ILP transcripts were confirmed, and quantitative RT-PCR was used to examine expression levels of these ILPs in different life stages, feeding states, and adult tissues. This work serves as an initial characterization of ILP expression in ticks and provides the foundation for further exploration of the roles of ILPs in these important arthropod vectors.

Tick Haller's Organ, a New Paradigm for Arthropod Olfaction: How Ticks Differ from Insects

Ticks are the vector of many human and animal diseases; and host detection is critical to this process. Ticks have a unique sensory structure located exclusively on the 1st pairs of legs; the fore-tarsal Haller’s organ, not found in any other animals, presumed to function like the insect antennae in chemosensation but morphologically very different. The mechanism of tick chemoreception is unknown. Utilizing next-generation sequencing and comparative transcriptomics between the 1st and 4th legs (the latter without the Haller’s organ), we characterized 1st leg specific and putative Haller’s organ specific transcripts from adult American dog ticks, Dermacentor variabilis. The analysis suggested that the Haller’s organ is involved in olfaction, not gustation. No known odorant binding proteins like those found in insects, chemosensory lipocalins or typical insect olfactory mechanisms were identified; with the transcriptomic data only supporting a possible olfactory G-protein coupled receptor (GPCR) signal cascade unique to the Haller’s organ. Each component of the olfactory GPCR signal cascade was identified and characterized. The expression of GPCR, G_αo and β-arrestin transcripts identified exclusively in the 1st leg transcriptome, and putatively Haller’s organ specific, were examined in unfed and blood-fed adult female and male D. variabilis. Blood feeding to repletion in adult females down-regulated the expression of all three chemosensory transcripts in females but not in males; consistent with differences in post-feeding tick behavior between sexes and an expected reduced chemosensory function in females as they leave the host. Data are presented for the first time of the potential hormonal regulation of tick chemosensation; behavioral assays confirmed the role of the Haller’s organ in N,N-diethyl-meta-toluamide (DEET) repellency but showed no role for the Haller’s organ in host attachment. Further research is needed to understand the potential role of the GPCR cascade in olfaction.

A repertoire of protease inhibitor families in Amblyomma americanum and other tick species: inter-species comparative analyses

BACKGROUND

Protease inhibitors (PIs) are important regulators of physiology and represent anti-parasitic druggable and vaccine targets. We conducted bioinformatic analyses of genome and transcriptome data to determine the protease inhibitor (PI) repertoire in Amblyomma americanum and in 25 other ixodid tick species. For A. americanum, we compared the PI repertoires in fed and unfed, male and female A. americanum ticks. We also analyzed PI repertoires of female 48, 96 and 120 h-fed midgut (MG) and salivary gland (SG) tissues.

RESULTS

We found 1,595 putative non-redundant PI sequences across 26 ixodid tick species. Ticks express PIs from at least 18 different families: I1, I2, I4, I8, I21, I25, I29, I31, I32, I35, I39, I43, I51, I53, I63, I68, I72 and I74 (MEROPS). The largest PI families were I2, I4 and I8 and lowest in I21, I31, I32, I35 and I68. The majority (75%) of tick PIs putatively inhibit serine proteases, with ~11 and 9% putatively regulating cysteine or metalloprotease-mediated pathways, respectively, and ~4% putatively regulating multiple/mixed protease types. In A. americanum, we found 370 PIs in female and 354 in male ticks. In A. americanum we found 231 and 442 in unfed and fed ticks, respectively. In females, we found 206 and 164 PIs in SG and MG, respectively. The majority of highly cross-tick species conserved PIs were in families I1, I2, I8, I21, I25, I29, I39 and I43.

CONCLUSIONS

Ticks appear to express large and diverse repertoires of PIs that primarily target serine protease-mediated pathways. We speculate that PI families with the highest repertoires may contain functionally redundant members while those with the lowest repertoires are functionally non-redundant PIs. We found some highly conserved PIs in the latter category, which we propose as potential candidates for broad-spectrum anti-tick vaccine candidates or druggable targets in tick control.

Transcription Factors as a Target for Vaccination Against Ticks and Mites

Ticks and mites are well-known ectoparasites as potential vectors for numerous bacteria, viruses, and parasites. Many being blood feeders add to physiological deterioration, morbidity, and mortality of their vertebrate hosts. To control them, transcription factors have been identified and studied in their role to sustain such arthropod pests. This paper summarizes some of the work done on those factors involved during blood feeding, reproduction, or when interacting with their pathogens and symbiont populations. Any transcription factor supporting the equilibrium developed by the ticks/mites could become a potential target for new control methods to prevent some of their key physiological functions.

Neuropeptide Evolution: Chelicerate Neurohormone and Neuropeptide Genes may reflect one or more whole genome duplications

Four genomes and two transcriptomes from six Chelicerate species were analyzed for the presence of neuropeptide and neurohormone precursors and their GPCRs. The genome from the spider Stegodyphus mimosarum yielded 87 neuropeptide precursors and 101 neuropeptide GPCRs. High neuropeptide transcripts were also found in the trancriptomes of three other spiders, Latrodectus hesperus, Parasteatoda tepidariorum and Acanthoscurria geniculata. For the scorpion Mesobuthus martensii the numbers are 79 and 74 respectively. The very small genome of the house dust mite, Dermatophagoides farinae, on the other hand contains a much smaller number of such genes. A few new putative Arthropod neuropeptide genes were discovered. Thus, both spiders and the scorpion have an achatin gene and in spiders there are two different genes encoding myosuppressin-like peptides while spiders also have two genes encoding novel LGamides. Another finding is the presence of trissin in spiders and scorpions, while neuropeptide genes that seem to be orthologs of Lottia LFRYamide and Platynereis CCRFamide were also found. Such genes were also found in various insect species, but seem to be lacking from the Holometabola. The Chelicerate neuropeptide and neuropeptide GPCR genes often have paralogs. As the large majority of these are probably not due to local gene duplications, is not impossible that they reflect the effects of one or more ancient whole genome duplications.

Neuropeptide evolution: Chelicerate neurohormone and neuropeptide genes may reflect one or more whole genome duplications

Four genomes and two transcriptomes from six Chelicerate species were analyzed for the presence of neuropeptide and neurohormone precursors and their GPCRs. The genome from the spider Stegodyphus mimosarum yielded 87 neuropeptide precursors and 120 neuropeptide GPCRs. Many neuropeptide transcripts were also found in the transcriptomes of three other spiders, Latrodectus hesperus, Parasteatoda tepidariorum and Acanthoscurria geniculata. For the scorpion Mesobuthus martensii the numbers are 79 and 93 respectively. The very small genome of the house dust mite, Dermatophagoides farinae, on the other hand contains a much smaller number of such genes. A few new putative Arthropod neuropeptide genes were discovered. Thus, both spiders and the scorpion have an achatin gene and in spiders there are two different genes encoding myosuppressin-like peptides while spiders also have two genes encoding novel LGamides. Another finding is the presence of trissin in spiders and scorpions, while neuropeptide genes that seem to be orthologs of Lottia LFRYamide and Platynereis CCRFamide were also found. Such genes were also found in various insect species, but seem to be lacking from the Holometabola. The Chelicerate neuropeptide and neuropeptide GPCR genes often have paralogs. As the large majority of these are probably not due to local gene duplications, is plausible that they reflect the effects of one or more ancient whole genome duplications.

Mevalonate-Farnesal Biosynthesis in Ticks: Comparative Synganglion Transcriptomics and a New Perspective

Juvenile hormone (JH) controls the growth, development, metamorphosis, and reproduction of insects. For many years, the general assumption has been that JH regulates tick and other acarine development and reproduction the same as in insects. Although researchers have not been able to find the common insect JHs in hard and soft tick species and JH applications appear to have no effect on tick development, it is difficult to prove the negative or to determine whether precursors to JH are made in ticks. The tick synganglion contains regions which are homologous to the corpora allata, the biosynthetic source for JH in insects. Next-gen sequencing of the tick synganglion transcriptome was conducted separately in adults of the American dog tick, Dermacentor variabilis, the deer tick, Ixodes scapularis, and the relapsing fever tick, Ornithodoros turicata as a new approach to determine whether ticks can make JH or a JH precursor. All of the enzymes that make up the mevalonate pathway from acetyl-CoA to farnesyl diphosphate (acetoacetyl-CoA thiolase, HMG-S, HMG-R, mevalonate kinase, phosphomevalonate kinase, diphosphomevalonate decarboxylase, and farnesyl diphosphate synthase) were found in at least one of the ticks studied but most were found in all three species. Sequence analysis of the last enzyme in the mevalonate pathway, farnesyl diphosphate synthase, demonstrated conservation of the seven prenyltransferase regions and the aspartate rich motifs within those regions typical of this enzyme. In the JH branch from farnesyl diphosphate to JH III, we found a putative farnesol oxidase used for the conversion of farnesol to farnesal in the synganglion transcriptome of I. scapularis and D. variabilis. Methyltransferases (MTs) that add a methyl group to farnesoic acid to make methyl farnesoate were present in all of the ticks studied with similarities as high as 36% at the amino acid level to insect JH acid methyltransferase (JHAMT). However, when the tick MTs were compared to the known insect JHAMTs from several insect species at the amino acid level, the former lacked the farnesoic acid binding motif typical in insects. The P450s shown in insects to add the C10,11 epoxide to methyl farnesoate, are in the CYP15 family; this family was absent in our tick transcriptomes and in the I. scapularis genome, the only tick genome available. These data suggest that ticks do not synthesize JH III but have the mevalonate pathway and may produce a JH III precursor.

Prediction of G protein-coupled receptor encoding sequences from the synganglion transcriptome of the cattle tick, Rhipicephalus microplus

The cattle tick, Rhipicephalus (Boophilus) microplus, is a pest which causes multiple health complications in cattle. The G protein-coupled receptor (GPCR) super-family presents a candidate target for developing novel tick control methods. However, GPCRs share limited sequence similarity among orthologous family members, and there is no reference genome available for R. microplus. This limits the effectiveness of alignment-dependent methods such as BLAST and Pfam for identifying GPCRs from R. microplus. However, GPCRs share a common structure consisting of seven transmembrane helices. We present an analysis of the R. microplus synganglion transcriptome using a combination of structurally-based and alignment-free methods which supplement the identification of GPCRs by sequence similarity. TMHMM predicts the number of transmembrane helices in a protein sequence. GPCRpred is a support vector machine-based method developed to predict and classify GPCRs using the dipeptide composition of a query amino acid sequence. These two bioinformatic tools were applied to our transcriptome assembly of the cattle tick synganglion. Together, BLAST and Pfam identified 85 unique contigs as encoding partial or full length candidate cattle tick GPCRs. Collectively, TMHMM and GPCRpred identified 27 additional GPCR candidates that BLAST and Pfam missed. This demonstrates that the addition of structurally-based and alignment-free bioinformatic approaches to transcriptome annotation and analysis produces a greater collection of prospective GPCRs than an analysis based solely upon methodologies dependent upon sequence alignment and similarity.

Genomic insights into the Ixodes scapularis tick vector of Lyme disease

Ticks transmit more pathogens to humans and animals than any other arthropod. We describe the 2.1 Gbp nuclear genome of the tick, Ixodes scapularis (Say), which vectors pathogens that cause Lyme disease, human granulocytic anaplasmosis, babesiosis and other diseases. The large genome reflects accumulation of repetitive DNA, new lineages of retro-transposons, and gene architecture patterns resembling ancient metazoans rather than pancrustaceans. Annotation of scaffolds representing ∼57% of the genome, reveals 20,486 protein-coding genes and expansions of gene families associated with tick-host interactions. We report insights from genome analyses into parasitic processes unique to ticks, including host 'questing', prolonged feeding, cuticle synthesis, blood meal concentration, novel methods of haemoglobin digestion, haem detoxification, vitellogenesis and prolonged off-host survival. We identify proteins associated with the agent of human granulocytic anaplasmosis, an emerging disease, and the encephalitis-causing Langat virus, and a population structure correlated to life-history traits and transmission of the Lyme disease agent.

IDENTIFICATION OF THE FGL-AMIDE ALLATOSTATIN GENE OF THE PRIMITIVE TERMITE Mastotermes darwiniensis AND THE WOODROACH Cryptocercus darwini

Allatostatins with the C-terminal ending Tyr/Phe-Xaa-Phe-Gly-Leu/Ile-amide (FGLa/ASTs) are widespread neuropeptides with multiple functions. The gene encoding the FGLa/AST polypeptide precursor was first isolated from cockroaches and since then could be identified in many insects and crustaceans. With its strictly conserved regions in combination with variable regions the gene seems to be a good candidate for phylogenetic analyses between closely and distantly related species. Here, the structure of the FGLa/AST gene of the most primitive termite, the giant northern termite Mastotermes darwiniensis Froggatt, was identified. The FGLa/AST gene of the woodroach Cryptocercus darwini was also determined. Precursor sequences of both species possess the general organization of dictyopteran FGLa/AST precursors containing 14 putative FGLa/AST peptides. In M. darwiniensis, only 11 out of the 14 FGLa/AST-like peptides possess the C-terminal conserved region Y/FXFGL/I/V/M and four of the putative peptide structures are not followed by a Gly residue that would lead to nonamidated peptides. Phylogenetic analyses show the high degree of similarity of dictyopteran FGLa/AST sequences. The position of termites, nested within the Blattaria, confirms that termites have evolved from primitive cockroaches.

Comparison of synganglion neuropeptides, neuropeptide receptors and neurotransmitter receptors and their gene expression in response to feeding in Ixodes scapularis (Ixodidae) vs. Ornithodoros turicata (Argasidae)

Illumina GAII high-throughput sequencing was used to compare expressed genes for female synganglion neuropeptides, neuropeptide receptors and neurotransmitter receptors of the soft tick Ornithodoros turicata with the hard tick Ixodes scapularis. Gene ontology molecular level three mapping revealed no significant differences amongst the same categories represented in O. turicata and I. scapularis. Transcripts predicting 22 neuropeptides or their receptors in the O. turicata synganglion were similar to annotations for 23 neuropeptides or receptors previously identified from I scapularis, with minor exceptions. A transcript predicting ecdysis triggering hormone receptor was identified in O. turicata; transcripts encoding for proprotein convertase and glycoprotein B were identified in both species. Transcripts predicting the same neurotransmitter receptors were found in the synganglion of both species. Gene expression of the transcripts showed numerous differences in response to feeding. Major differences were observed in expression of genes believed important in regulating slow vs. rapid feeding, blood water elimination, cuticle synthesis plasticity and in signalling reproductive activity. Although the glutamate receptor was strongly upregulated in both species, the gamma aminobutyric acid receptor, which inhibits glutamate, was upregulated significantly only in I. scapularis. These differences are consistent with the slow vs. rapid action of the pharyngeal pump in the two species.

Neuropeptide discovery in Symphylella vulgaris (Myriapoda, Symphyla): In silico prediction of the first myriapod peptidome

Arthropods have contributed greatly to our understanding of peptidergic control of physiology and behavior, and being the largest and most diverse animal phylum, represent a model for investigating peptide hormone evolution. Surprisingly, one arthropod subphylum, the Myriapoda, is uninvestigated in terms of its peptide hormones. The public deposition of a transcriptome for Symphylella vulgaris, a pseudocentipede, provides a means for peptide discovery in myriapods. Here, in silico transcriptome mining was used to identify 47 S. vulgaris neuropeptide-encoding transcripts within this dataset. The identified transcripts allowed for the deduction of 31 unique pre/preprohormone sequences, with 97 distinct mature peptides predicted from the deduced proteins. The predicted S. vulgaris peptidome includes members of the adipokinetic hormone/red pigment concentrating hormone, adipokinetic hormone-corazonin-like peptide, allatostatin A, allatostatin C (AST-C), allatotropin, CCHamide, crustacean cardioactive peptide, GSEFLamide, insulin-like peptide, intocin, proctolin, pyrokinin, short neuropeptide F, SIFamide and sulfakinin families. This is the first, and thus far only, peptidome predicted for a myriapod. Of particular note were a modified AST-C, TYWKQCAFNAVSRFamide, that lacks one of two cysteine residues (i.e. one at position 13) stereotypically present in members of this peptide family (and hence is missing the disulfide bridge that spans these residues) and a SIFamide, PPFNGSIFamide, that is truncated due to a lysine for arginine substitution in the dibasic residue pair commonly located at positions 3 and 4 of stereotypical full-length isoforms (e.g. the crustacean peptide GYRKPPFNGSIFamide). The peptides predicted here represent the only extant resource for initiating investigations of native peptidergic signaling in the Myriapoda.

Molecular and pharmacological characterization of the Chelicerata pyrokinin receptor from the southern cattle tick, Rhipicephalus (Boophilus) microplus

We identified the first pyrokinin receptor (Rhimi-PKR) in Chelicerata and analyzed structure-activity relationships of cognate ligand neuropeptides and their analogs. Based on comparative and phylogenetic analyses, this receptor, which we cloned from larvae of the cattle tick Rhipicephalus microplus (Acari: Ixodidae), is the ortholog of the insect pyrokinin (PK)/pheromone biosynthesis activating neuropeptide (PBAN)/diapause hormone (DH) neuropeptide family receptor. Rhimi-PKR functional analyses using calcium bioluminescence were performed with a developed stable recombinant CHO-K1 cell line. Rhimi-PKR was activated by four endogenous PKs from the Lyme disease vector, the tick Ixodes scapularis (EC50s range: 85.4 nM-546 nM), and weakly by another tick PRX-amide peptide, periviscerokinin (PVK) (EC50 = 24.5 μM). PK analogs with substitutions of leucine, isoleucine or valine at the C-terminus for three tick PK peptides, Ixosc-PK1, Ixosc-PK2, and Ixosc-PK3, retained their potency on Rhimi-PKR. Therefore, Rhimi-PKR is less selective and substantially more tolerant than insect PK receptors of C-terminal substitutions of leucine to isoleucine or valine, a key structural feature that serves to distinguish insect PK from PVK/CAP2b receptors. In females, ovary and synganglion had the highest Rhimi-PKR relative transcript abundance followed by the rectal sac, salivary glands, Malpighian tubules, and midgut. This is the first pharmacological analysis of a PK/PBAN/DH-like receptor from the Chelicerata, which will now permit the discovery of the endocrinological roles of this neuropeptide family in vectors of vertebrate pathogens.

Orcokinin-like immunoreactivity in central neurons innervating the salivary glands and hindgut of ixodid ticks

Orcokinins are conserved neuropeptides within the Arthropoda but their cellular distribution and functions in ticks are unknown. We use an antibody against the highly conserved N-terminal (NFDEIDR) of mature orcokinin peptides to examine their distribution in six ixodid species: Amblyomma variegatum, Dermacentor reticulatus, Hyalomma anatolicum, Ixodes scapularis, Ixodes ricinus and Rhipicephalus appendiculatus. Numerous immunoreactive neurons (~100) were detected in various regions of the synganglion (central nervous system) in all examined tick species. Immunoreactive projections of two prominent groups of efferent neurons in the post-oesophageal region were examined in detail: (1) neurons innervating the salivary glands; (2) neurons innervating the hindgut. Using matrix-assisted laser desorption/ionisation-time-of-flight (MALDI-TOF), we detected orcokinin peaks in extracts of the synganglia and hindguts but not in the salivary glands of I. scapularis females. Our data provide further evidence of the presence of orcokinin in ixodid ticks and establish a morphological basis for functional studies of identified peptidergic neuronal networks.

In silico characterization of the neuropeptidome of the Western black widow spider Latrodectus hesperus

Technological advancements in high-throughput sequencing have resulted in the production/public deposition of an ever-growing number of arthropod transcriptomes. While most sequencing projects have focused on hexapods, transcriptomes have also been generated for members of the Chelicerata. One chelicerate for which a large transcriptome has recently been released is the Western black widow Latrodectus hesperus, a member of the Araneae (true spiders). Here, a neuropeptidome for L. hesperus was predicted using this resource. Thirty-eight peptide-encoding transcripts were mined from the L. hesperus transcriptome, with 216 distinct peptides predicted from the deduced pre/preprohormones. The identified peptides included members of the allatostatin A, allatostatin B, allatostatin C, allatotropin, bursicon α, bursicon β, CAPA/periviscerokinin/pyrokinin, CCHamide, corazonin, crustacean cardioactive peptide, crustacean hyperglycemic hormone/ion transport peptide, diuretic hormone 31, diuretic hormone 44, FMRFamide-like peptide (FLP), GSEFLamide, insulin-like peptide, neuropeptide F (NPF), orcokinin, proctolin, short neuropeptide F, SIFamide, sulfakinin and tachykinin-related peptide (TRP) families. Of particular note were the identifications of a carboxyl (C)-terminally extended corazonin, FLPs possessing -IMRFamide, -MMYFamide, and -MIHFamide C-termini, a NPF and a sulfakinin each ending in -RYamide rather than -RFamide, a precursor whose orcokinins include C-terminally amidated isoforms, and a collection of TRPs possessing -FXPXLamide rather than the stereotypical -FXGXLamide C-termini. The L. hesperus peptidome is by far the largest thus far published for any member of the Chelicerata. Taken collectively, these data serve as a reference for future neuropeptide discovery in the Araneae and provide a foundation for future studies of peptidergic control in L. hesperus and other spiders.

Transcriptome of the female synganglion of the black-legged tick Ixodes scapularis (Acari: Ixodidae) with comparison between Illumina and 454 systems

Illumina and 454 pyrosequencing were used to characterize genes from the synganglion of female Ixodes scapularis. GO term searching success for biological processes was similar for samples sequenced by both methods. However, for molecular processes, it was more successful for the Illumina samples than for 454 samples. Functional assignments of transcripts predicting neuropeptides, neuropeptide receptors, neurotransmitter receptors and other genes of interest was done, supported by strong e-values (<-6), and high consensus sequence alignments. Transcripts predicting 15 putative neuropeptide prepropeptides ((allatostatin, allatotropin, bursicon α, corticotropin releasing factor (CRF), CRF-binding protein, eclosion hormone, FMRFamide, glycoprotein A, insulin-like peptide, ion transport peptide, myoinhibitory peptide, inotocin ( =  neurophysin-oxytocin), Neuropeptide F, sulfakinin and SIFamide)) and transcripts predicting receptors for 14 neuropeptides (allatostatin, calcitonin, cardioacceleratory peptide, corazonin, CRF, eclosion hormone, gonadotropin-releasing hormone/AKH-like, insulin-like peptide, neuropeptide F, proctolin, pyrokinin, SIFamide, sulfakinin and tachykinin) are reported. Similar to Dermacentor variabilis, we found transcripts matching pro-protein convertase, essential for converting neuropeptide hormones to their mature form. Additionally, transcripts predicting 6 neurotransmitter/neuromodulator receptors (acetylcholine, GABA, dopamine, glutamate, octopamine and serotonin) and 3 neurotransmitter transporters (GABA transporter, noradrenalin-norepinephrine transporter and Na+-neurotransmitter/symporter) are described. Further, we found transcripts predicting genes for pheromone odorant receptor, gustatory receptor, novel GPCR messages, ecdysone nuclear receptor, JH esterase binding protein, steroidogenic activating protein, chitin synthase, chitinase, and other genes of interest. Also found were transcripts predicting genes for spermatogenesis-associated protein, major sperm protein, spermidine oxidase and spermidine synthase, genes not normally expressed in the female CNS of other invertebrates. The diversity of messages predicting important genes identified in this study offers a valuable resource useful for understanding how the tick synganglion regulates important physiological functions.

Neuropeptidergic control of the hindgut in the black-legged tick Ixodes scapularis

The hindgut, as a part of the tick excretory system, plays an important physiological role in maintaining homoeostases and waste elimination. Immunoreactive projections from the synganglion to the hindgut were found using antibodies against four different neuropeptides: FGLamide related allatostatin, myoinhibitory peptide, SIFamide, and orcokinin. The presence of FGLamide related allatostatin, myoinhibitory peptide and SIFamide in both synganglia (source) and hindgut (target organ) extracts was confirmed by MALDI-TOF. Tissue-specific PCR revealed the expression of four putative FGLamide related allatostatin receptors and an SIFamide receptor in the hindgut. An antibody against Ixodes scapularis SIFamide receptor detected immunoreactive spots in epithelial cells as well as the visceral muscles surrounding the rectal sac, while staining with the antibody against myoinhibitory peptide receptor 1 revealed that the immunoreactivity was only associated with the visceral muscles. In hindgut motility assays, SIFamide activated hindgut motility in a dose-dependent manner. None of other three neuropeptides (FGLamide related allatostatin, myoinhibitory peptide and orcokinin) activated hindgut motility when tested alone. Myoinhibitory peptide antagonised the SIFamide-stimulated hindgut mobility when it was tested in combination with SIFamide.

Examining the role of foraging and malvolio in host-finding behavior in the honey bee parasite, Varroa destructor (Anderson & Trueman)

When a female varroa mite, Varroa destructor (Anderson & Trueman), invades a honey bee brood cell, the physiology rapidly changes from feeding phoretic to reproductive. Changes in foraging and malvolio transcript levels in the brain have been associated with modulated intra-specific food searching behaviors in insects and other invertebrates. Transcription profiles for both genes were examined during and immediately following brood cell invasion to assess their role as potential control elements. Vdfor and Vdmvl transcripts were found in all organs of varroa mites with the highest Vdfor transcript levels in ovary-lyrate organs and the highest Vdmvl in Malpighian tubules. Changes in transcript levels of Vdfor and Vdmvl in synganglia were not associated with the cell invasion process, remaining comparable between early reproductive mites (collected from the pre-capping brood cells) and phoretic mites. However, Vdfor and Vdmvl transcript levels were lowered by 37 and 53%, respectively, in synganglia from reproductive mites compared to early reproductive mites, but not significantly different to levels in synganglia from phoretic mites. On the other hand, in whole body preparations the Vdfor and Vdmvl had significantly higher levels of transcript in reproductive mites compared to phoretic and early reproductive, mainly due to the presence of both transcripts accumulating in the eggs carried by the ovipositing mite. Varroa mites are a critical component for honey bee population decline and finding varroa mite genes associated with brood cell invasion, reproduction, ion balance and other physiological processes will facilitate development of novel control avenues for this honey bee parasite.

Molecular and functional characterization of the first tick CAP2b (periviscerokinin) receptor from Rhipicephalus (Boophilus) microplus (Acari: Ixodidae)

The cDNA of the receptor for CAP(2b)/periviscerokinin (PVK) neuropeptides, designated Rhimi-CAP(2b)-R, was cloned from synganglia of tick Rhipicephalus (Boophilus) microplus. This receptor is the ortholog of the insect CAP(2b)/PVK receptor, as concluded from analyses of the predicted protein sequence, phylogenetics and functional expression. Expression analyses of synganglion, salivary gland, Malpighian tubule, and ovary revealed Rhimi-CAP(2b)-R transcripts. The expression in mammalian cells of the open reading frame of Rhimi-CAP(2b)-R cDNA fused with a hemagglutinin tag at the receptor N-terminus was confirmed by immunocytochemistry. In a calcium bioluminescence assay the recombinant receptor was activated by the tick Ixodes scapularis CAP(2b)/PVK and a PVK analog with EC₅₀s of 64 nM and 249 nM, respectively. Tick pyrokinins were not active. This is the first report on the functional characterization of the CAP(2b)/PVK receptor from any tick species which will now permit the discovery of the physiological roles of these neuropeptides in ticks, as neurohormones, neuromodulators and/or neurotransmitters.

Genomic organization and reproductive regulation of a large lipid transfer protein in the varroa mite, Varroa destructor (Anderson & Trueman)

The complete genomic region and corresponding transcript of the most abundant protein in phoretic varroa mites, Varroa destructor (Anderson & Trueman), were sequenced and have homology with acarine hemelipoglycoproteins and the large lipid transfer protein (LLTP) super family. The genomic sequence of VdLLTP included 14 introns and the mature transcript coded for a predicted polypeptide of 1575 amino acid residues. VdLLTP shared a minimum of 25% sequence identity with acarine LLTPs. Phylogenetic assessment showed VdLLTP was most closely related to Metaseiulus occidentalis vitellogenin and LLTP proteins of ticks; however, no heme binding by VdLLTP was detected. Analysis of lipids associated with VdLLTP showed that it was a carrier for free and esterified C12 -C22 fatty acids from triglycerides, diacylglycerides and monoacylglycerides. Additionally, cholesterol and β-sitosterol were found as cholesterol esters linked to common fatty acids. Transcript levels of VdLLTP were 42 and 310 times higher in phoretic female mites when compared with males and quiescent deutonymphs, respectively. Coincident with initiation of the reproductive phase, VdLLTP transcript levels declined to a third of those in phoretic female mites. VdLLTP functions as an important lipid transporter and should provide a significant RNA interference target for assessing the control of varroa mites.

Molecular characterization of tick salivary gland glutaminyl cyclase

Glutaminyl cyclase (QC) catalyzes the cyclization of N-terminal glutamine residues into pyroglutamate. This post-translational modification extends the half-life of peptides and, in some cases, is essential in binding to their cognate receptor. Due to its potential role in the post-translational modification of tick neuropeptides, we report the molecular, biochemical and physiological characterization of salivary gland QC during the prolonged blood feeding of the black-legged tick (Ixodes scapularis) and the gulf-coast tick (Amblyomma maculatum). QC sequences from I. scapularis and A. maculatum showed a high degree of amino acid identity to each other and other arthropods and residues critical for zinc binding/catalysis (D159, E202, and H330) or intermediate stabilization (E201, W207, D248, D305, F325, and W329) are conserved. Analysis of QC transcriptional gene expression kinetics depicts an upregulation during the bloodmeal of adult female ticks prior to fast-feeding phases in both I. scapularis and A. maculatum suggesting a functional link with bloodmeal uptake. QC enzymatic activity was detected in saliva and extracts of tick salivary glands and midguts. Recombinant QC was shown to be catalytically active. Furthermore, knockdown of QC transcript by RNA interference resulted in lower enzymatic activity, and small, unviable egg masses in both studied tick species as well as lower engorged tick weights for I. scapularis. These results suggest that the post-translational modification of neurotransmitters and other bioactive peptides by QC is critical to oviposition and potentially other physiological processes. Moreover, these data suggest that tick-specific QC-modified neurotransmitters/hormones or other relevant parts of this system could potentially be used as novel physiological targets for tick control.

In silico cloning of genes encoding neuropeptides, neurohormones and their putative G-protein coupled receptors in a spider mite

The genome of the spider mite was prospected for the presence of genes coding neuropeptides, neurohormones and their putative G-protein coupled receptors. Fifty one candidate genes were found to encode neuropeptides or neurohormones. These include all known insect neuropeptides and neurohormones, with the exception of sulfakinin, corazonin, neuroparsin and PTTH. True orthologs of adipokinetic hormone (AKH) were neither found, but there are three genes encoding peptides similar in structure to both AKH and the AKH-corazonin-related peptide. We were also unable to identify the precursors for pigment dispersing factor (PDF) or the recently discovered trissin. However, the spider mite probably does have such genes, as we found their putative receptors. A novel arthropod neuropeptide gene was identified that shows similarity to previously described molluscan neuropeptide genes and was called EFLamide. A total of 65 putative neuropeptide GPCR genes were also identified, of these 58 belong to the A-family and 7 to the B-family. Phylogenetic analysis showed that 50 of them are closely related to insect GPCRs, which allowed the identification of their putative ligand in 39 cases with varying degrees of certainty. Other spider mite GPCRs however have no identifiable orthologs in the genomes of the four holometabolous insect species best analyzed. Whereas some of the latter have orthologs in hemimetabolous insect species, crustaceans or ticks, for others such arthropod homologs are currently unknown.

An evolutionary comparison of leucine-rich repeat containing G protein-coupled receptors reveals a novel LGR subtype

Leucine-rich repeat containing G protein-coupled receptors or LGRs are receptors with important functions in development and reproduction. Belonging to this evolutionarily conserved group of receptors are the well-studied glycoprotein hormone receptors and relaxin receptors in mammals, as well as the bursicon receptor, which triggers cuticle hardening and tanning in freshly enclosed insects. In this study, the numerous LGR sequences in different animal phyla are analyzed and compared. Based on these data a phylogenetic tree was generated. This information sheds new light on structural and evolutionary aspects regarding this receptor group. Apart from vertebrates and insects, LGRs are also present in early chordates (Urochordata, Cephalochordata and Hyperoartia) and other arthropods (Arachnida and Branchiopoda) as well as in Mollusca, Echinodermata, Hemichordata, Nematoda, and even in ancient animal life forms, such as Cnidaria and Placozoa. Three distinct types of LGR exist, distinguishable by their number of leucine-rich repeats (LRRs), their type-specific hinge region and the presence or absence of an LDLa motif. Type C LGRs containing only one LDLa (C1 subtype) appear to be present in nearly all animal phyla. We here describe a second subtype, C2, containing multiple LDLa motifs, which was discovered in echinoderms, mollusks and in one insect species (Pediculus humanis corporis). In addition, eight putative LGRs can be predicted from the genome data of the placozoan species Trichoplax adhaerens. They may represent an ancient form of the LGRs, however, more genomic data will be required to confirm this hypothesis.

Akt is an essential player in regulating cell/organ growth at the adult stage in the hard tick Haemaphysalis longicornis

Ticks grow rapidly during blood feeding, and their body weight may ultimately increase 100-fold more than that before feeding. The molecular mechanisms controlling growth during blood feeding in ticks remain largely unknown. The conserved insulin/PI3K/Akt signaling pathway regulates growth and metabolism in eukaryotes. Here, we show evidence for the involvement of Akt in growth during blood feeding in the parthenogenetic strain of the hard tick Haemaphysalis longicornis. We identified a homolog of the Ser/Thr kinase Akt (HlAkt) from the EST database of the H. longicornis embryo. HlAkt cDNA had a 1,590 bp ORF that encodes 529 amino acids with a predicted molecular weight of 60 kDa. HlAkt possesses a PH domain, a Ser/Thr kinase domain, a hydrophobic motif, and dual phosphorylation residues (Thr 338 and Ser 503) that are essential for kinase activation. Knockdown of HlAkt by RNA interference caused inhibition of blood feeding in female ticks. Histological observation demonstrated that HlAkt knockdown led to the arrest of growth in internal organs. HlAkt knockdown also affected the expressions of blood meal-induced genes that are essential for blood digestion, development, and reproduction in the female tick. These results strongly indicate that HlAkt is essential to complete the blood feeding process accompanied by the growth of internal organs in adult ticks. This is the first report of identification and characterization of Akt in Chelicerata, including ticks.

First transcriptome of the testis-vas deferens-male accessory gland and proteome of the spermatophore from Dermacentor variabilis (Acari: Ixodidae)

Ticks are important vectors of numerous human diseases and animal diseases. Feeding stimulates spermatogenesis, mating and insemination of male factors that trigger female reproduction. The physiology of male reproduction and its regulation of female development are essentially a black box. Several transcriptomes have catalogued expression of tick genes in the salivary glands, synganglion and midgut but no comprehensive investigation has addressed male reproduction and mating. Consequently, a new global approach using transcriptomics, proteomics, and quantitative gene expression is needed to understand male reproduction and stimulation of female reproduction.

This first transcriptome to the reproductive biology of fed male ticks, Dermacentor variabilis, was obtained by 454 pyrosequencing (563,093 reads, 12,804 contigs). Gene Ontology (Biological Processes level III) recognized 3,866 transcripts in 73 different categories; spermiogenesis; spermatogenesis; peptidases, lipases and hydrolases; oxidative and environmental stress; immune defense; and protein binding. Reproduction-associated genes included serine/threonine kinase, metalloendoproteinases, ferritins, serine proteases, trypsin, cysteine proteases, serpins, a cystatin, GPCR and others. qRT-PCR showed significant upregulation from unfed versus fed adult male reproductive organs of zinc metalloprotease, astacin metalloprotease and serine protease, enzymes important in spermiogenesis and mating activity in insects, as well as a GPCR with the greatest similarity to a SIFamide receptor known to be important in regulating courtship behavior in Drosophila. Proteomics on these organs and the spermatophore by tryptic digestion/Liquid chromatography/Mass spectrometry/Mass spectrometry (LC/MS/MS) demonstrated expression of many of the same messages found by 454 sequencing, supporting their identification, and revealed differences in protein distribution in the reproductive system versus the spermatophore. We found Efα but no EF β in the transcriptome and neither of these proteins in the spermatophore. Thus, the previously described model for male regulation of female reproduction may not apply to other ticks. A new paradigm is needed to explain male stimulation of female tick reproduction.

Synganglion transcriptome and developmental global gene expression in adult females of the American dog tick, Dermacentor variabilis (Acari: Ixodidae)

454 Pyrosequencing was used to characterize the expressed genes from the synganglion and associated neurosecretory organs of unfed and partially fed virgin and mated replete females of the American dog tick, Dermacentor variabilis. A total of 14,881 contiguous sequences (contigs) was assembled, with an average size of 229 bp. Gene ontology terms for Level 2 biological processes were assigned to 4366 contigs. Seven acetylcholinesterases, a muscarinic acetylcholine (ACh) receptor, two nicotinic ACh receptor β-subunits, two ACh unc-18 regulators, two dopamine receptors, two gamma aminobutyric acid (GABA) receptors, two GABA transporters, two norepinephrine transporters and an octopamine receptor are described. Microarrays were conducted to examine global gene expression and quantitative real-time polymerase chain reaction was used to verify expression of selected neuropeptides. Hierarchical clustering of all differentially expressed transcripts grouped part-fed and replete ticks as being more similar in terms of differentially expressed genes with unfed ticks as the outgroup. Nine putative neuropeptides (allatostatin, bursicon-β, preprocorazonin, glycoprotein hormone α, insulin-like peptide, three orcokinins, preprosulphakinin) and a gonadotropin releasing hormone receptor were differentially expressed, and their developmental expression and role in reproduction was investigated. The presence of eclosion hormone, corazonin and bursicon in the synganglion, which in insects regulate behaviour and cuticle development associated with moulting, suggest that this system may be used in ticks to regulate blood feeding, cuticle expansion and development related to female reproduction; adult ticks do not moult.

The fat body transcriptomes of the yellow fever mosquito Aedes aegypti, pre- and post- blood meal

BACKGROUND

The fat body is the main organ of intermediary metabolism in insects and the principal source of hemolymph proteins. As part of our ongoing efforts to understand mosquito fat body physiology and to identify novel targets for insect control, we have conducted a transcriptome analysis of the fat body of Aedes aegypti before and in response to blood feeding.

RESULTS

We created two fat body non-normalized EST libraries, one from mosquito fat bodies non-blood fed (NBF) and another from mosquitoes 24 hrs post-blood meal (PBM). 454 pyrosequencing of the non-normalized libraries resulted in 204,578 useable reads from the NBF sample and 323,474 useable reads from the PBM sample. Alignment of reads to the existing reference Ae. aegypti transcript libraries for analysis of differential expression between NBF and PBM samples revealed 116,912 and 115,051 matches, respectively. De novo assembly of the reads from the NBF sample resulted in 15,456 contigs, and assembly of the reads from the PBM sample resulted in 15,010 contigs. Collectively, 123 novel transcripts were identified within these contigs. Prominently expressed transcripts in the NBF fat body library were represented by transcripts encoding ribosomal proteins. Thirty-five point four percent of all reads in the PBM library were represented by transcripts that encode yolk proteins. The most highly expressed were transcripts encoding members of the cathepsin b, vitellogenin, vitellogenic carboxypeptidase, and vitelline membrane protein families.

CONCLUSION

The two fat body transcriptomes were considerably different from each other in terms of transcript expression in terms of abundances of transcripts and genes expressed. They reflect the physiological shift of the pre-feeding fat body from a resting state to vitellogenic gene expression after feeding.

Bioinformatic prediction of arthropod/nematode-like peptides in non-arthropod, non-nematode members of the Ecdysozoa

The Onychophora, Priapulida and Tardigrada, along with the Arthropoda, Nematoda and several other small phyla, form the superphylum Ecdysozoa. Numerous peptidomic studies have been undertaken for both the arthropods and nematodes, resulting in the identification of many peptides from each group. In contrast, little is known about the peptides used as paracrines/hormones by species from the other ecdysozoan taxa. Here, transcriptome mining and bioinformatic peptide prediction were used to identify peptides in members of the Onychophora, Priapulida and Tardigrada, the only non-arthropod, non-nematode members of the Ecdysozoa for which there are publicly accessible expressed sequence tags (ESTs). The extant ESTs for each phylum were queried using 106 arthropod/nematode peptide precursors. Transcripts encoding calcitonin-like diuretic hormone and pigment-dispersing hormone (PDH) were identified for the onychophoran Peripatopsis sedgwicki, with transcripts encoding C-type allatostatin (C-AST) and FMRFamide-like peptide identified for the priapulid Priapulus caudatus. For the Tardigrada, transcripts encoding members of the A-type allatostatin, C-AST, insect kinin, orcokinin, PDH and tachykinin-related peptide families were identified, all but one from Hypsibius dujardini (the exception being a Milnesium tardigradum orcokinin-encoding transcript). The proteins deduced from these ESTs resulted in the prediction of 48 novel peptides, six onychophoran, eight priapulid and 34 tardigrade, which are the first described from these phyla.

New approach for the study of mite reproduction: The first transcriptome analysis of a mite, Phytoseiulus persimilis (Acari: Phytoseiidae)

Many species of mites and ticks are of agricultural and medical importance. Much can be learned from the study of transcriptomes of acarines which can generate DNA-sequence information of potential target genes for the control of acarine pests. High throughput transcriptome sequencing can also yield sequences of genes critical during physiological processes poorly understood in acarines, i.e., the regulation of female reproduction in mites. The predatory mite, Phytoseiulus persimilis, was selected to conduct a transcriptome analysis using 454 pyrosequencing. The objective of this project was to obtain DNA-sequence information of expressed genes from P. persimilis with special interest in sequences corresponding to vitellogenin (Vg) and the vitellogenin receptor (VgR). These genes are critical to the understanding of vitellogenesis, and they will facilitate the study of the regulation of mite female reproduction. A total of 12,556 contiguous sequences (contigs) were assembled with an average size of 935bp. From these sequences, the putative translated peptides of 11 contigs were similar in amino acid sequences to other arthropod Vgs, while 6 were similar to VgRs. We selected some of these sequences to conduct stage-specific expression studies to further determine their function.

Identification of chelicerate neuropeptides using bioinformatics of publicly accessible expressed sequence tags

While numerous investigations have focused on the identification of neuropeptides in arthropods, most have been conducted on members of the Hexapoda or Crustacea, and little is currently known about those in the Chelicerata. Here, publicly accessible expressed sequence tags (ESTs) were mined for putative chelicerate neuropeptide-encoding transcripts; the peptides encoded by the ESTs were deduced using on-line peptide prediction programs and homology to known isoforms. Fifty-eight ESTs representing eight peptide families/subfamilies were identified using this strategy. Of note was the prediction of the first authentic chelicerate C-type allatostatin, pQIRYHQCYFNPISCF, from the mite Tetranychus urticae, as well as the prediction a novel allatostatin CC peptide, GEGKMFWRCYFNAVSCF, from both the tick Amblyomma variegatum and the scorpion Mesobuthus gibbosus. Also identified from T. urticae were authentic crustacean cardioactive peptide (CCAP), several peptides belonging to the crustacean hyperglycemic hormone/ion transport peptide superfamily, members of the calcitonin-like diuretic hormone/diuretic hormone 31 family, and several FMRFamide-like peptides, specifically members of the neuropeptide F (NPF) and short neuropeptide F subfamilies. To the best of our knowledge the identifications of CCAP and NPF in T. urticae are the first for the Chelicerata. In addition, several novel orcokinins were identified from the scorpion Scorpiops jendeki and the spider Loxosceles laeta; in S. jendeki previously unknown isoforms of SIFamide, ESRNPPLNGSMFamide and ESKNPPLNGSMFamide, were also predicted. Taken collectively, the data presented in our study expand the catalog of known chelicerate neuropeptides and provide a foundation for future physiological studies of them in these animals.