Genome organization of major tandem repeats in the hard tick, Ixodes scapularis

Insights into the Karyotype Evolution of Charinidae, the Early-Diverging Clade of Whip Spiders (Arachnida: Amblypygi)

Whip spiders (Amblypygi) represent an ancient order of tetrapulmonate arachnids with a low diversity. Their cytogenetic data are confined to only a few reports. Here, we analyzed the family Charinidae, a lineage almost at the base of the amblypygids, providing an insight into the ancestral traits and basic trajectories of amblypygid karyotype evolution. We performed Giemsa staining, selected banding techniques, and detected 18S ribosomal DNA and telomeric repeats by fluorescence in situ hybridization in four Charinus and five Sarax species. Both genera exhibit a wide range of diploid chromosome numbers (2n = 42-76 and 22-74 for Charinus and Sarax, respectively). The 2n reduction was accompanied by an increase of proportion of biarmed elements. We further revealed a single NOR site (probably an ancestral condition for charinids), the presence of a (TTAGG)n telomeric motif localized mostly at the chromosome ends, and an absence of heteromorphic sex chromosomes. Our data collectively suggest a high pace of karyotype repatterning in amblypygids, with probably a high ancestral 2n and its subsequent gradual reduction by fusions, and the action of pericentric inversions, similarly to what has been proposed for neoamblypygids. The possible contribution of fissions to charinid karyotype repatterning, however, cannot be fully ruled out.

Karyotype changes in long-term cultured tick cell lines

Tick cell lines are an easy-to-handle system for the study of viral and bacterial infections and other aspects of tick cellular processes. Tick cell cultures are often continuously cultivated, as freezing can affect their viability. However, the long-term cultivation of tick cells can influence their genome stability. In the present study, we investigated karyotype and genome size of tick cell lines. Though 16S rDNA sequencing showed the similarity between Ixodes spp. cell lines at different passages, their karyotypes differed from 2n = 28 chromosomes for parental Ixodes spp. ticks, and both increase and decrease in chromosome numbers were observed. For example, the highly passaged Ixodes scapularis cell line ISE18 and Ixodes ricinus cell lines IRE/CTVM19 and IRE/CTVM20 had modal chromosome numbers 48, 23 and 48, respectively. Also, the Ornithodoros moubata cell line OME/CTVM22 had the modal chromosome number 33 instead of 2n = 20 chromosomes for Ornithodoros spp. ticks. All studied tick cell lines had a larger genome size in comparison to the genomes of the parental ticks. Thus, highly passaged tick cell lines can be used for research purposes, but possible differences in encoded genetic information and downstream cellular processes, between different cell populations, should be taken into account.

Ixodes persulcatus/pavlovskyi natural hybrids in Siberia: Occurrence in sympatric areas and infection by a wide range of tick-transmitted agents

Ixodes persulcatus and Ixodes pavlovskyi ticks, two closely related species of the I. ricinus - I. persulcatus group, are widely distributed in the southern part of Western Siberia. Recently, the existence of natural hybrids of I. persulcatus and I. pavlovskyi ticks has been demonstrated. The aim of this study was to evaluate the abundance of I. persulcatus/pavlovskyi hybrids in several locations with different ratios of parental tick species and to investigate the prevalence and genetic variability of a wide range of infectious agents in these hybrids compared to the parental tick species. Natural hybrids of I. persulcatus and I. pavlovskyi ticks were identified in all examined locations in Altai and Novosibirsk, Western Siberia, Russia. The abundance of hybrids varied from 7% to 40% in different locations and was maximal in a location with similar proportions of I. persulcatus and I. pavlovskyi ticks. For the first time, it was shown that hybrids can be infected with the same agents as their parental tick species: tick-borne encephalitis and Kemerovo viruses, Borrelia afzelii, Borrelia bavariensis, Borrelia garinii, Borrelia miyamotoi, Rickettsia helvetica, Rickettsia raoultii, Rickettsia sibirica, "Candidatus Rickettsia tarasevichiae", Anaplasma phagocytophilum, Ehrlichia muris, "Candidatus Neoehrlichia mikurensis", and Babesia microti. The prevalence of most bacterial agents in hybrids was intermediate compared to their parental tick species. Most genetic variants of the identified agents have been previously found in the parental tick species. Wide distribution of I. persulcatus/pavlovskyi natural hybrids implies that I. persulcatus, I. pavlovskyi and their hybrids coexist in all I. persulcatus - I. pavlovskyi sympatric areas.

A Roadmap for Tick-Borne Flavivirus Research in the "Omics" Era

Tick-borne flaviviruses (TBFs) affect human health globally. Human vaccines provide protection against some TBFs, and antivirals are available, yet TBF-specific control strategies are limited. Advances in genomics offer hope to understand the viral complement transmitted by ticks, and to develop disruptive, data-driven technologies for virus detection, treatment, and control. The genome assemblies of Ixodes scapularis, the North American tick vector of the TBF, Powassan virus, and other tick vectors, are providing insights into tick biology and pathogen transmission and serve as nucleation points for expanded genomic research. Systems biology has yielded insights to the response of tick cells to viral infection at the transcript and protein level, and new protein targets for vaccines to limit virus transmission. Reverse vaccinology approaches have moved candidate tick antigenic epitopes into vaccine development pipelines. Traditional drug and in silico screening have identified candidate antivirals, and target-based approaches have been developed to identify novel acaricides. Yet, additional genomic resources are required to expand TBF research. Priorities include genome assemblies for tick vectors, “omic” studies involving high consequence pathogens and vectors, and emphasizing viral metagenomics, tick-virus metabolomics, and structural genomics of TBF and tick proteins. Also required are resources for forward genetics, including the development of tick strains with quantifiable traits, genetic markers and linkage maps. Here we review the current state of genomic research on ticks and tick-borne viruses with an emphasis on TBFs. We outline an ambitious 10-year roadmap for research in the “omics era,” and explore key milestones needed to accomplish the goal of delivering three new vaccines, antivirals and acaricides for TBF control by 2030.

Genome scaffolding and annotation for the pathogen vector Ixodes ricinus by ultra-long single molecule sequencing

Background

Global warming and other ecological changes have facilitated the expansion of Ixodes ricinus tick populations. Ixodes ricinus is the most important carrier of vector-borne pathogens in Europe, transmitting viruses, protozoa and bacteria, in particular Borrelia burgdorferi (sensu lato), the causative agent of Lyme borreliosis, the most prevalent vector-borne disease in humans in the Northern hemisphere. To faster control this disease vector, a better understanding of the I. ricinus tick is necessary. To facilitate such studies, we recently published the first reference genome of this highly prevalent pathogen vector. Here, we further extend these studies by scaffolding and annotating the first reference genome by using ultra-long sequencing reads from third generation single molecule sequencing. In addition, we present the first genome size estimation for I. ricinus ticks and the embryo-derived cell line IRE/CTVM19.

Results

235,953 contigs were integrated into 204,904 scaffolds, extending the currently known genome lengths by more than 30% from 393 to 516 Mb and the N50 contig value by 87% from 1643 bp to a N50 scaffold value of 3067 bp. In addition, 25,263 sequences were annotated by comparison to the tick’s North American relative Ixodes scapularis. After (conserved) hypothetical proteins, zinc finger proteins, secreted proteins and P450 coding proteins were the most prevalent protein categories annotated. Interestingly, more than 50% of the amino acid sequences matching the homology threshold had 95–100% identity to the corresponding I. scapularis gene models. The sequence information was complemented by the first genome size estimation for this species. Flow cytometry-based genome size analysis revealed a haploid genome size of 2.65Gb for I. ricinus ticks and 3.80 Gb for the cell line.

Conclusions

We present a first draft sequence map of the I. ricinus genome based on a PacBio-Illumina assembly. The I. ricinus genome was shown to be 26% (500 Mb) larger than the genome of its American relative I. scapularis. Based on the genome size of 2.65 Gb we estimated that we covered about 67% of the non-repetitive sequences. Genome annotation will facilitate screening for specific molecular pathways in I. ricinus cells and provides an overview of characteristics and functions.

Electronic supplementary material

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

Paralog analyses reveal gene duplication events and genes under positive selection in Ixodes scapularis and other ixodid ticks

Background

Hard ticks (family Ixodidae) are obligatory hematophagous ectoparasites of worldwide medical and veterinary importance. The haploid genomes of multiple species of ixodid ticks exceed 1 Gbp, prompting questions regarding gene, segmental and whole genome duplication in this phyletic group. The availability of the genome assembly for the black legged tick, Ixodes scapularis, and transcriptome datasets for multiple species of ticks offers an opportunity to assess the contribution of gene duplication to the genome. Here we developed a bioinformatics pipeline to identify and analyze duplicated genes (paralogs) using gene models from the prostriate tick, I. scapularis IscaW1.1 annotation and expressed sequence tags (ESTs) from I. scapularis and the metastriate ticks, Rhipicephalus microplus (southern cattle tick), R. appendiculatus (brown ear tick) and Amblyomma variegatum (tropical bont tick).

Results

Approximately 1-2 % of I. scapularis gene models and 2-14 % of ESTs from the four species represent duplicated genes. The ratio of non-synonymous to synonymous nucleotide substitution rates suggests ~ 25 % of duplicated genes are under positive selection pressure in each species. Analyses of synonymous substitution rates provide evidence for two duplication events in I. scapularis and R. microplus involving several hundred genes. Conservative molecular clock estimates based on synonymous substitution rates for species of Anopheles mosquitoes and the fruit fly, Drosophila melanogaster, suggest these events occurred within the last 50 MYA. Mapping of paralogs to the I. scapularis genome assembly supports tandem, or possibly segmental duplication events.

Conclusions

The present study marks the first genome-level analyses of gene duplication for the Ixodidae and provides insights into mechanisms shaping genome evolution in this group. At least two duplication events involving hundreds of genes may have occurred independently in the lineages prostriata and metastriata, with tandem and segmental duplication the most likely mechanisms for paralog generation. Duplicated genes under positive selection pressure may be linked to emerging functions in the tick and represent important candidates for further study.

Electronic supplementary material

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

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.

Natural hybridization between Ixodes ricinus and Ixodes persulcatus ticks evidenced by molecular genetics methods

The recently shown phenomenon of natural hybridization between Ixodes persulcatus and Ixodes pavlovskyi ticks (Kovalev et al., 2015) stimulated similar studies in the sympatric zones of other tick species. In the present paper, 265 Ixodes ricinus and I. persulcatus ticks from Estonia were subjected to a search for interspecific hybrids based on nuclear (ITS2) and mitochondrial (cox1) markers as well as morphological features. Surprisingly, only 72.1% of ticks morphologically identified as I. ricinus actually were I. ricinus both at nuclear and mitochondrial markers, while the accuracy of morphological species identification for I. persulcatus was 99.3%. Among ticks morphologically identified as I. ricinus, 24.6% turned out to be interspecific hybrids and 3.3% were I. persulcatus. Generally, about 11% of the individuals studied were shown to be interspecific hybrids with different levels of nuclear DNA introgression. The analysis of hybrid populations proved the mating pair female I. ricinus×male I. persulcatus to form hybrids more efficiently, then female I. persulcatus×male I. ricinus. The same trend can be observed for backcrosses preferentially mating with I. ricinus. Hybridization between I. ricinus and I. persulcatus proved the existing view about their reproductive isolation to be untenable. Interspecific hybridization occurring between both closely (I. persulcatus and I. pavlovskyi) and more distantly (I. ricinus and I. persulcatus) related Ixodes species could introduce novel alleles that modify vector competence, host use or the ability to exploit diverse microhabitats.

Heterochromatin, histone modifications, and nuclear architecture in disease vectors

Interactions between a pathogen and a vector are plastic and dynamic. Such interactions can be more rapidly accommodated by epigenetic changes than by genetic mutations. Gene expression can be affected by the proximity to the heterochromatin, by local histone modifications, and by the three-dimensional position within the nucleus. Recent studies of disease vectors indicate that gene regulation by these factors can be important for susceptibility to pathogens, reproduction, immunity, development, and longevity. Knowledge about heterochromatin, histone modifications, and nuclear architecture will help our understanding of epigenetic mechanisms that control gene function at traits related to vectorial capacity.

An Ixodes scapularis cell line with a predominantly neuron-like phenotype

The Ixodes scapularis embryo-derived cell line ISE6 is the most widely utilized tick-derived cell line due to its susceptibility to a wide variety of tick- and non-tick-vectored pathogens. Little is known about its tissue origin or biological background. Protein expression of ISE6 cells was compared with that of another I. scapularis-derived cell line, IDE12, and dissected tick synganglia. Results demonstrated the presence of a neuronal marker protein, type 3 β-tubulin, in all three samples, as well as other shared and unique neuronal and immune response-associated proteins. Of neuronal proteins shared between the two cell lines, ISE6 expressed several in significantly greater quantities than IDE12. Stimulation of ISE6 cells by in vivo exposure to the hemocoel environment in unfed larval and molting nymphal ticks, but not unfed nymphal ticks, resulted in the development of neuron-like morphologic characteristics in the implanted cells.

Proteomics characterization of tick-host-pathogen interactions

Ticks are blood-feeding arthropod ectoparasites of wild and domestic animals that transmit disease-causing pathogens to humans and animals worldwide and a good model for the characterization of tick-host-pathogen interactions. Tick-host-pathogen interactions consist of dynamic processes involving genetic traits of hosts, pathogens, and ticks that mediate their development and survival. Proteomics provides information on the protein content of cells and tissues that may differ from results at the transcriptomics level and may be relevant for basic biological studies and vaccine antigen discovery. In this chapter, we describe various methods for protein extraction and for proteomics analysis in ticks based on one-dimensional gel electrophoresis to characterize tick-host-pathogen interactions. Particularly relevant for this characterization is the use of blood-fed ticks. Therefore, we put special emphasis on working with replete ticks collected after feeding on vertebrate hosts.

Origin and distribution of AT-rich repetitive DNA families in Triatoma infestans (Heteroptera)

Triatoma infestans, one of the most important vectors of Trypanosoma cruzi, is very interesting model, because it shows large interpopulation variation in the amount and distribution of heterochromatin. This polymorphism involved the three large pairs up to almost all autosomal pairs, including the sex chromosomes. To understand the dynamics of heterochromatin variation in T. infestans, we isolated the AT-rich satDNA portion of this insect using reassociation kinetics (C0t), followed by cloning, sequencing and FISH. After chromosome localization, immunolabeling with anti-5-methylcytosine, anti-H4K5ac and anti-H3K9me2 antibodies was performed to determine the functional characteristics of heterochromatin. The results allowed us to reorganize the karyotype of T. infestans in accordance with the distribution of the families of repetitive DNA using seven different markers. We found that two arrays with lengths of 79 and 33bp have a strong relationship with transposable element sequences, suggesting that these two families of satDNA probably originated from Polintons. The results also allowed us to identify at least four chromosome rearrangements involved in the amplification/dispersion of AT-rich satDNA of T. infestans. These data should be very useful in new studies including those examining the cytogenomic and population aspects of this very important species of insect.

Comparative proteomics for the characterization of the most relevant Amblyomma tick species as vectors of zoonotic pathogens worldwide

Ticks transmit zoonotic pathogens worldwide. Nevertheless, very little information is available on their genome, transcriptome and proteome. Herein, we characterized the proteome of Amblyomma americanum adults and nymphs because of their role in pathogen transmission and compared the proteome of A. americanum, A. cajennense and A. variegatum adult ticks. We also used de novo sequencing proteomics data for the analysis of the phylogenetic relationships between the three Amblyomma spp. in a proof of concept for phyloproteomics. The results showed that host and tick proteins involved in blood digestion, heme detoxification, development and innate immunity were differentially represented between adults and nymphs. Although these ticks were unfed, over-represented host proteins may supply nutrients during off-host periods. Tick proteins involved in tick attachment, feeding, heat shock response, protease inhibition and heme detoxification were differentially represented between Amblyomma spp., suggesting adaptation processes to biotic and abiotic factors. These results suggested that phyloproteomics might be a useful tool for the phylogenetic analysis of tick species in which sequence data is a limiting factor and demonstrate the possibilities of proteomics studies for the characterization of relevant tick vector species and provide new relevant information to understand the physiology, development and evolution of these tick species.

BIOLOGICAL SIGNIFICANCE

This is the first report on the proteome of the most important Amblyomma tick species for their relevance as vectors of zoonotic pathogens worldwide. Nevertheless, very little information is available on the genome, transcriptome and proteome of these vector ectoparasites. The results reported herein provide new relevant information to understand the physiology, development and evolution of these tick species. Phyloproteomics using de novo protein sequencing was assayed as a new approach for the phylogenetic analysis of tick species in which sequence data is a limiting factor.This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Development of genomic resources for the tick Ixodes ricinus: isolation and characterization of single nucleotide polymorphisms

Assessing the genetic variability of the tick Ixodes ricinus-an important vector of pathogens in Europe-is an essential step for setting up antitick control methods. Here, we report the first identification of a set of SNPs isolated from the genome of I. ricinus, by applying a reduction in genomic complexity, pyrosequencing and new bioinformatics tools. Almost 1.4 million of reads (average length: 528 nt) were generated with a full Roche 454 GS FLX run on two reduced representation libraries of I. ricinus. A newly developed bioinformatics tool (DiscoSnp), which isolates SNPs without requiring any reference genome, was used to obtain 321 088 putative SNPs. Stringent selection criteria were applied in a bioinformatics pipeline to select 1768 SNPs for the development of specific primers. Among 384 randomly SNPs tested by Fluidigm genotyping technology on 464 individuals ticks, 368 SNPs loci (96%) exhibited the presence of the two expected alleles. Hardy-Weinberg equilibrium tests conducted on six natural populations of ticks have shown that from 26 to 46 of the 384 loci exhibited significant heterozygote deficiency.

Host specialization in ticks and transmission of tick-borne diseases: a review

Determining patterns of host use, and the frequency at which these patterns change, are of key importance if we are to understand tick population dynamics, the evolution of tick biodiversity, and the circulation and evolution of associated pathogens. The question of whether ticks are typically host specialists or host generalists has been subject to much debate over the last half-century. Indeed, early research proposed that morphological diversity in ticks was linked to host specific adaptations and that most ticks were specialists. Later work disputed this idea and suggested that ticks are largely limited by biogeographic conditions and tend to use all locally available host species. The work presented in this review suggests that the actual answer likely lies somewhere between these two extremes. Although recent observational studies support the view that phylogenetically diverse host species share ticks when found on similar ecological ranges, theory on host range evolution predicts that host specialization should evolve in ticks given their life history characteristics. Contemporary work employing population genetic tools to examine host-associated population structure in several tick systems support this prediction and show that simple species records are not enough to determine whether a parasite is a true host generalist; host specialization does evolve in ticks at local scales, but may not always lead to speciation. Ticks therefore seem to follow a pattern of being global generalists, local specialists. Given this, the notion of host range needs to be modified from an evolutionary perspective, where one simply counts the number of hosts used across the geographic distribution, to a more ecological view, where one considers host use at a local scale, if we are to better understand the circulation of tick-borne pathogens and exposure risks for humans and livestock.

Sex chromosome pairing and extensive NOR polymorphism in Wadicosa fidelis (Araneae: Lycosidae)

In terms of cytogenetics, entelegyne araneomorphs are the best studied clade of spiders. The typical karyotype of entelegyne males consists of acrocentric chromosomes, including 2 non-homologous X chromosomes. The present study is focused on the karyotype, nucleolus organising regions (NORs) and sex chromosome behaviour during meiosis of the entelegyne Wadicosa fidelis (Lycosidae). Preparations stained by Giemsa were used to study karyotype and meiosis. NORs were visualised by silver staining and fluorescence in situ hybridisation with 18S rDNA probe. The male karyotype consists of 28 acrocentric elements, including 2 X chromosomes. In contrast to the majority of other spiders, the male sex chromosomes pair during the major part of meiosis. Following an initial period of parallel pairing, the attachment of male sex chromosomes is restricted to centromeric areas and continues until metaphase II. Our study revealed an enormous number of NORs in the population from Galilee and indicates a considerable variability of NOR numbers in this population. The distal regions of 9 or 10 autosomal pairs contain NORs. The obtained data indicate the rapid spread of NORs in the karyotype of W. fidelis, which was presumably caused by ectopic recombinations and subsequent hybridisations of individuals with different NOR genotypes that produced heterozygotes.

Why so many unknown genes? Partitioning orphans from a representative transcriptome of the lone star tick Amblyomma americanum

Background

Genomic resources within the phylum Arthropoda are largely limited to the true insects but are beginning to include unexplored subphyla, such as the Crustacea and Chelicerata. Investigations of these understudied taxa uncover high frequencies of orphan genes, which lack detectable sequence homology to genes in pre-existing databases. The ticks (Acari: Chelicerata) are one such understudied taxon for which genomic resources are urgently needed. Ticks are obligate blood-feeders that vector major diseases of humans, domesticated animals, and wildlife. In analyzing a transcriptome of the lone star tick Amblyomma americanum, one of the most abundant disease vectors in the United States, we find a high representation of unannotated sequences. We apply a general framework for quantifying the origin and true representation of unannotated sequences in a dataset and for evaluating the biological significance of orphan genes.

Results

Expressed sequence tags (ESTs) were derived from different life stages and populations of A. americanum and combined with ESTs available from GenBank to produce 14,310 ESTs, over twice the number previously available. The vast majority (71%) has no sequence homology to proteins archived in UniProtKB. We show that poor sequence or assembly quality is not a major contributor to this high representation by orphan genes. Moreover, most unannotated sequences are functional: a microarray experiment demonstrates that 59% of functional ESTs are unannotated. Lastly, we attempt to further annotate our EST dataset using genomic datasets from other members of the Acari, including Ixodes scapularis, four other tick species and the mite Tetranychus urticae. We find low homology with these species, consistent with significant divergence within this subclass.

Conclusions

We conclude that the abundance of orphan genes in A. americanum likely results from 1) taxonomic isolation stemming from divergence within the tick lineage and limited genomic resources for ticks and 2) lineage-specific genes needing functional genomic studies to evaluate their association with the unique biology of ticks. The EST sequences described here will contribute substantially to the development of tick genomics. Moreover, the framework provided for the evaluation of orphan genes can guide analyses of future transcriptome sequencing projects.

Molecular and pharmacological characterization of two D(1)-like dopamine receptors in the Lyme disease vector, Ixodes scapularis

Advancements in tick neurobiology may impact the development of acaricides to control those species that transmit human and animal diseases. Here, we report the first cloning and pharmacological characterization of two neurotransmitter binding G protein-coupled receptors in the Lyme disease (blacklegged) tick, Ixodes scapularis. The genes IscaGPRdop1 and IscaGPRdop2 were identified in the I. scapularis genome assembly and predicted as orthologs of previously characterized D(1)-like dopamine receptors in the fruit fly Drosophila melanogaster and honeybee Apis mellifera. Heterologous expression in HEK 293 cells demonstrated that each receptor functioned as a D(1)-like dopamine receptor because significant increases in levels of intracellular cyclic adenosine monophosphate (cAMP) were detected following dopamine treatment. Importantly, the receptors were distinct in their pharmacological properties regarding concentration-dependent response to dopamine, constitutive activity, and response to other biogenic amines. Exposure to a variety of dopamine receptor agonists and antagonists further demonstrated a D(1)-like pharmacology of these dopamine receptors and highlighted their differential activities in vitro.

Differential transcription of two highly divergent gut-expressed Bm86 antigen gene homologues in the tick Rhipicephalus appendiculatus (Acari: Ixodida)

The transcriptional control of gene expression is not well documented in the Arthropoda. We describe transcriptional analysis of two exceptionally divergent homologues (Ra86) of the Bm86 gut antigen from Rhipicephalus appendiculatus. Bm86 forms the basis of a commercial vaccine for the control of Rhipicephalus (Boophilus) microplus. The R. appendiculatus Ra86 proteins contain 654 and 693 amino acids, with only 80% amino acid sequence identity. Reverse-transcription PCR of gut cDNA showed transcription of only one genotype in individual female ticks. PCR amplification of 3' untranslated sequences from genomic DNA indicated that both variants could be encoded within a single genome. When both variants were present, one of the two Ra86 genotypes was transcriptionally dominant.