Genome Sequence of the Endosymbiont Rickettsia peacockii and Comparison with Virulent Rickettsia rickettsii: Identification of Virulence Factors

Rickettsia peacockii, also known as the East Side Agent, is a non-pathogenic obligate intracellular bacterium found as an endosymbiont in Dermacentor andersoni ticks in the western USA and Canada. Its presence in ticks is correlated with reduced prevalence of Rickettsia rickettsii, the agent of Rocky Mountain Spotted Fever. It has been proposed that a virulent SFG rickettsia underwent changes to become the East Side Agent. We determined the genome sequence of R. peacockii and provide a comparison to a closely related virulent R. rickettsii. The presence of 42 chromosomal copies of the ISRpe1 transposon in the genome of R. peacockii is associated with a lack of synteny with the genome of R. rickettsii and numerous deletions via recombination between transposon copies. The plasmid contains a number of genes from distantly related organisms, such as part of the glycosylation island of Pseudomonas aeruginosa. Genes deleted or mutated in R. peacockii which may relate to loss of virulence include those coding for an ankyrin repeat containing protein, DsbA, RickA, protease II, OmpA, ScaI, and a putative phosphoethanolamine transferase. The gene coding for the ankyrin repeat containing protein is especially implicated as it is mutated in R. rickettsii strain Iowa, which has attenuated virulence. Presence of numerous copies of the ISRpe1 transposon, likely acquired by lateral transfer from a Cardinium species, are associated with extensive genomic reorganization and deletions. The deletion and mutation of genes possibly involved in loss of virulence have been identified by this genomic comparison. It also illustrates that the introduction of a transposon into the genome can have varied effects; either correlating with an increase in pathogenicity as in Francisella tularensis or a loss of pathogenicity as in R. peacockii and the recombination enabled by multiple transposon copies can cause significant deletions in some genomes while not in others.     

Lethal Effect of Rickettsia rickettsii on Its Tick Vector (Dermacentor andersoni)

Rickettsia rickettsii, the causative agent of Rocky Mountain spotted fever, was lethal for the majority of experimentally and transovarially infected Rocky Mountain wood ticks (Dermacentor andersoni). Overall, 94.1% of nymphs infected as larvae by feeding on rickettsemic guinea pigs died during the molt into adults and 88.3% of adult female ticks infected as nymphs died prior to feeding. In contrast, only 2.8% of uninfected larvae failed to develop into adults over two generations. Infected female ticks incubated at 4°C had a lower mortality (80.9%) than did those held at 21°C (96.8%). Rickettsiae were vertically transmitted to 39.0% of offspring, and significantly fewer larvae developed from infected ticks. The lethal effect of R. rickettsii may explain the low prevalence of infected ticks in nature and affect its enzootic maintenance.     

Comparative genomics of four Liliales families inferred from the complete chloroplast genome sequence of Veratrum patulum O. Loes. (Melanthiaceae)

The sequence of the chloroplast genome, which is inherited maternally, contains useful information for many scientific fields such as plant systematics, biogeography and biotechnology because its characteristics are highly conserved among species. There is an increase in chloroplast genomes of angiosperms that have been sequenced in recent years. In this study, the nucleotide sequence of the chloroplast genome (cpDNA) of Veratrum patulum Loes. (Melanthiaceae, Liliales) was analyzed completely. The circular double-stranded DNA of 153,699 bp consists of two inverted repeat (IR) regions of 26,360 bp each, a large single copy of 83,372 bp, and a small single copy of 17,607 bp. This plastome contains 81 protein-coding genes, 30 distinct tRNA and four genes of rRNA. In addition, there are six hypothetical coding regions (ycf1, ycf2, ycf3, ycf4, ycf15 and ycf68) and two open reading frames (ORF42 and ORF56), which are also found in the chloroplast genomes of the other species. The gene orders and gene contents of the V. patulum plastid genome are similar to that of Smilax china, Lilium longiflorum and Alstroemeria aurea, members of the Smilacaceae, Liliaceae and Alstroemeriaceae (Liliales), respectively. However, the loss rps16 exon 2 in V. patulum results in the difference in the large single copy regions in comparison with other species. The base substitution rate is quite similar among genes of these species. Additionally, the base substitution rate of inverted repeat region was smaller than that of single copy regions in all observed species of Liliales. The IR regions were expanded to trnH_GUG in V. patulum, a part of rps19 in L. longiflorum and A. aurea, and whole sequence of rps19 in S. china. Furthermore, the IGS lengths of rbcL-accD-psaI region were variable among Liliales species, suggesting that this region might be a hotspot of indel events and the informative site for phylogenetic studies in Liliales. In general, the whole chloroplast genome of V. patulum, a potential medicinal plant, will contribute to research on the genetic applications of this genus.     

New chloroplast microsatellite markers suitable for assessing genetic diversity of Lolium perenne and other related grass species

Background and Aims

Lolium perenne (perennial ryegrass) is the most important forage grass species of temperate regions. We have previously released the chloroplast genome sequence of L. perenne ‘Cashel’. Here nine chloroplast microsatellite markers are published, which were designed based on knowledge about genetically variable regions within the L. perenne chloroplast genome. These markers were successfully used for characterizing the genetic diversity in Lolium and different grass species.

Methods

Chloroplast genomes of 14 Poaceae taxa were screened for mononucleotide microsatellite repeat regions and primers designed for their amplification from nine loci. The potential of these markers to assess genetic diversity was evaluated on a set of 16 Irish and 15 European L. perenne ecotypes, nine L. perenne cultivars, other Lolium taxa and other grass species.

Key Results

All analysed Poaceae chloroplast genomes contained more than 200 mononucleotide repeats (chloroplast simple sequence repeats, cpSSRs) of at least 7 bp in length, concentrated mainly in the large single copy region of the genome. Nucleotide composition varied considerably among subfamilies (with Pooideae biased towards poly A repeats). The nine new markers distinguish L. perenne from all non-Lolium taxa. TeaCpSSR28 was able to distinguish between all Lolium species and Lolium multiflorum due to an elongation of an A₈ mononucleotide repeat in L. multiflorum. TeaCpSSR31 detected a considerable degree of microsatellite length variation and single nucleotide polymorphism. TeaCpSSR27 revealed variation within some L. perenne accessions due to a 44-bp indel and was hence readily detected by simple agarose gel electrophoresis. Smaller insertion/deletion events or single nucleotide polymorphisms detected by these new markers could be visualized by polyacrylamide gel electrophoresis or DNA sequencing, respectively.

Conclusions

The new markers are a valuable tool for plant breeding companies, seed testing agencies and the wider scientific community due to their ability to monitor genetic diversity within breeding pools, to trace maternal inheritance and to distinguish closely related species.     

Ixodid fauna and zoonotic agents in ticks from dogs: first report of Rickettsia rickettsii in Rhipicephalus sanguineus in the state of Mato Grosso do Sul, mid-western Brazil

Ticks from 148 dogs from the urban area of the municipality of Campo Grande, state of Mato Grosso do Sul, Brazil, were collected, classified and analyzed using polymerase chain reaction (PCR) for the identification of Rickettsia spp., Trypanosoma cruzi and Leishmania spp. A total of 2015 ticks were collected. The species Rhipicephalus sanguineus (98.9 %) and Amblyomma cajennense (1.1 %) were identified. Molecular analysis revealed that no tick samples were infected by T. cruzi. Regarding Leishmania spp., tick samples from 36 dogs spread across all regions of the municipality were positive for L. chagasi. One tick sample was positive for Rickettsia spp. (gltA gene) in the PCR reaction. This sample was submitted to further PCR based on the ompA gene and the amplicon was sequenced. Identity of 100 % was found with homologous sequences of R. rickettsii available in GenBank. This paper is the first to report the natural infection of R. sanguineus by R. rickettsii in the municipality of Campo Grande, state of Mato Grosso do Sul, mid-western Brazil.     

Transformation Frequency of a mariner-Based Transposon in Rickettsia rickettsii

Transformation frequencies of a mariner-based transposon system in Rickettsia rickettsii were determined using a plaque assay system for enumeration and isolation of mutants. Sequence analysis of insertion sites in both R. rickettsii and R. prowazekii indicated that insertions were random. Transposon mutagenesis provides a useful tool for rickettsial research.     

Genetic Variation in the Plasmodium falciparum Circumsporozoite Protein in India and Its Relevance to RTS,S Malaria Vaccine

RTS,S is the most advanced malaria vaccine candidate, currently under phase-III clinical trials in Africa. This Plasmodium falciparum vaccine contains part of the central repeat region and the complete C-terminal T cell epitope region (Th2R and Th3R) of the circumsporozoite protein (CSP). Since naturally occurring polymorphisms at the vaccine candidate loci are critical determinants of the protective efficacy of the vaccines, it is imperative to investigate these polymorphisms in field isolates. In this study we have investigated the genetic diversity at the central repeat, C-terminal T cell epitope (Th2R and Th3R) and N-terminal T cell epitope regions of the CSP, in P. falciparum isolates from Madhya Pradesh state of India. These isolates were collected through a 5-year prospective study aimed to develop a well-characterized field-site for the future evaluation of malaria vaccine in India. Our results revealed that the central repeat (63 haplotypes, n = 161) and C-terminal Th2R/Th3R epitope (24 haplotypes, n = 179) regions were highly polymorphic, whereas N-terminal non-repeat region was less polymorphic (5 haplotypes, n = 161) in this population. We did not find any evidence of the role of positive natural selection in maintaining the genetic diversity at the Th2R/Th3R regions of CSP. Comparative analysis of the Th2R/Th3R sequences from this study to the global isolates (n = 1160) retrieved from the GenBank database revealed two important points. First, the majority of the sequences (∼61%, n = 179) from this study were identical to the Dd2/Indochina type, which is also the predominant Th2R/Th3R haplotype in Asia (∼59%, n = 974). Second, the Th2R/Th3R sequences in Asia, South America and Africa are geographically distinct with little allele sharing between continents. In conclusion, this study provides an insight on the existing polymorphisms in the CSP in a parasite population from India that could potentially influence the efficacy of RTS,S vaccine in this region.     

ISSR and DAMD markers revealed high genetic variability within Flavoparmelia caperata in Western Himalaya (India)

Flavoparmelia caperata (L.) Hale is medicinally very important and possesses antifungal and antibacterial activities. F. caperata is the only species found in India. Inter simple sequence repeat (ISSR) and Directed amplification of minisatellite DNA (DAMD) methods were used to analyze the genetic variability within F. caperata from the Western Himalayan region of India. Eleven ISSR and 10 DAMD primers produced 139 and 117 polymorphic bands, and detected 91.44 and 82.34 % polymorphisms, respectively. Cumulative band data generated for ISSR and DAMD markers resulted in 86.86 % polymorphism across all the accessions of F. caperata. The average Polymorphic information content (PIC) value obtained with ISSR, DAMD, and cumulative band data were 0.28, 0.27, and 0.27, respectively. The clustering of the F. caperata accessions in the UPGMA dendrogram showed that these accessions are intermingled with each other in different subclusters irrespective of their geographical affiliations. The pattern of genetic variations within F. caperata accessions could be due to free exchange of spores that might have taken place among these accessions in the wild. ISSR and DAMD markers efficiently and reliably resulted in discrete banding patterns and polymorphic profiles. These markers despite targeting different regions of genome, revealed almost similar levels of polymorphism across all the accessions. The wide range of genetic distance and high level of polymorphism detected by ISSR and DAMD reflected a high genetic variability among the different accessions of F. caperata.     

Sequence polymorphisms in ribosomal RNA genes and variations in chromosomal loci of Oenothera odorata and O. laciniata

Numerous studies on Oenothera species have been investigated for the physiological and ecological characteristics. However, no such an information based on molecular cytogenetic has yet been introduced, in turn, is very essential for identifying sequence polymorphisms of rRNA genes with their loci on mitotic phases for further biological researches. In this study, sequence variations of rRNA genes in Oenothera odorata and O. laciniata were examined to identify informative factors as unique or repeat sequences in intra- and inter-specific variations. Intra-specific variation revealed that the sequences of the rRNA genes including the spacer regions were highly conserved revealing only a few variations. From the inter-specific variation, spacer regions of species were completely different as (1) non-homologous sequences in NTS and (2) different type repeat sequences in ITS 1, 2 and 5.8S rRNA, whereas the remaining coding regions were highly conserved. FISH was carried out on mitotic phases using the 5S rDNA of the analyzed sequences. From the interphase and metaphase chromosomes of the examined species, two loci of 5S rDNA in O. odorata and four loci in O. laciniata were confirmed on the telomeric region of the short arm. Due to the small size and unclear centromere of the chromosomes, karyotype could not be completed. However, we confirmed that the chromosomes are organized by meta- and acrocentric chromosomes and the chromosomes with identified loci were assumed to be paired by the location of loci at the telomeric region on the short arm with relative lengths.     

Intraspecific independent evolution of floral spur length in response to local flower visitor size in Japanese Aquilegia in different mountain regions

Geographic differences in floral traits may reflect geographic differences in effective pollinator assemblages. Independent local adaptation to pollinator assemblages in multiple regions would be expected to cause parallel floral trait evolution, although sufficient evidence for this is still lacking. Knowing the intraspecific evolutionary history of floral traits will reveal events that occur in the early stages of trait diversification. In this study, we investigated the relationship between flower spur length and pollinator size in 16 populations of Aquilegia buergeriana var. buergeriana distributed in four mountain regions in the Japanese Alps. We also examined the genetic relationship between yellow‐ and red‐flowered individuals, to see if color differences caused genetic differentiation by pollinator isolation. Genetic relationships among 16 populations were analyzed based on genome‐wide single‐nucleotide polymorphisms. Even among populations within the same mountain region, pollinator size varied widely, and the average spur length of A. buergeriana var. buergeriana in each population was strongly related to the average visitor size of that population. Genetic relatedness between populations was not related to the similarity of spur length between populations; rather, it was related to the geographic proximity of populations in each mountain region. Our results indicate that spur length in each population evolved independently of the population genetic structure but in parallel in response to local flower visitor size in different mountain regions. Further, yellow‐ and red‐flowered individuals of A. buergeriana var. buergeriana were not genetically differentiated. Unlike other Aquilegia species in Europe and America visited by hummingbirds and hawkmoths, the Japanese Aquilegia species is consistently visited by bumblebees. As a result, genetic isolation by flower color may not have occurred.     

Nucleotide Sequence of Polymerase Chain Reaction Product Amplified from Rickettsia japonica DNA Using Rickettsia rickettsii 190-Kilodalton Surface Antigen Gene Primers

The PCR product amplified from Rickettsia japonica with the primer pair Rr 190.70p and Rr 190.602n of R. rickettsii 190-kDa antigen gene was cloned into M13mp19 RF DNA at the EcoRI site and sequenced by chemiluminescent DNA sequencing. The sequence revealed a molecular size of 533 base pairs (bp). The primer-flanking region of 491 bp, an open reading frame, was compared with the corresponding region of R. rickettsii, demonstrating 35 nucleotide substitutions in R. japonica. The sequence of primer portions in R. japonica DNA was also analyzed, revealing one nucleotide substitution in the Rr 190.70p and two in the Rr 190.602n portion. The homology in the overall sequence of PCR-amplified regions between R. japonica and R. rickettsii was 93% in nucleotide and 85% in putative amino acid structure. The sequence contains no cleavage site for the restriction endonuclease AfaI but two PstI sites giving three fragments of 121, 159, and 253 bp, which differentiated R. japonica from other spotted fever group rickettsiae in addition to R. rickettsii. The cleavage sites for endonucleases AluI, HinfI, and MunI that disappeared or appeared in the sequence by nucleotide substitution differentiated R. japonica from others, as did PstI. The estimation of molecular size of DNA fragments on polyacrylamide gel electrophoresis is discussed.     

Tandem repeat copy-number variation in protein-coding regions of human genes

Background

Tandem repeat variation in protein-coding regions will alter protein length and may introduce frameshifts. Tandem repeat variants are associated with variation in pathogenicity in bacteria and with human disease. We characterized tandem repeat polymorphism in human proteins, using the UniGene database, and tested whether these were associated with host defense roles.

Results

Protein-coding tandem repeat copy-number polymorphisms were detected in 249 tandem repeats found in 218 UniGene clusters; observed length differences ranged from 2 to 144 nucleotides, with unit copy lengths ranging from 2 to 57. This corresponded to 1.59% (218/13,749) of proteins investigated carrying detectable polymorphisms in the copy-number of protein-coding tandem repeats. We found no evidence that tandem repeat copy-number polymorphism was significantly elevated in defense-response proteins (p = 0.882). An association with the Gene Ontology term 'protein-binding' remained significant after covariate adjustment and correction for multiple testing. Combining this analysis with previous experimental evaluations of tandem repeat polymorphism, we estimate the approximate mean frequency of tandem repeat polymorphisms in human proteins to be 6%. Because 13.9% of the polymorphisms were not a multiple of three nucleotides, up to 1% of proteins may contain frameshifting tandem repeat polymorphisms.

Conclusion

Around 1 in 20 human proteins are likely to contain tandem repeat copy-number polymorphisms within coding regions. Such polymorphisms are not more frequent among defense-response proteins; their prevalence among protein-binding proteins may reflect lower selective constraints on their structural modification. The impact of frameshifting and longer copy-number variants on protein function and disease merits further investigation.     

Vectorial competence of Amblyomma tonelliae to transmit Rickettsia rickettsii

The aim of this work was to test the vectorial competence of Amblyomma tonelliae (Ixodida: Ixodidae) to transmit Rickettsia rickettsii (Rickettsiales: Rickettsiaceae), the agent of Rocky Mountain spotted fever (RMSF). All parasitic stages of A. tonelliae were exposed to R. rickettsii by allowing each stage to feed on hosts inoculated with this pathogen. Thereafter, ticks were fed on uninfected hosts. All stages of A. tonelliae were able to acquire the R. rickettsii infection and maintain it by transstadial and transovarial transmission. When infected ticks fed on uninfected hosts, the hosts developed rickettsiosis disease. This study demonstrates the vectorial competence of A. tonelliae to transmit R. rickettsii. These results have epidemiological relevance because A. tonelliae is one of the tick species most likely to infest humans in Argentina, including in areas in which RMSF has been reported.     

Association of Bitter Taste Receptor T2R38 Polymorphisms,Oral Microbiota,and Rheumatoid Arthritis

The association of taste genetics and the oral microbiome in autoimmune diseases such as rheumatoid arthritis (RA) has not been reported. We explored a novel oral mucosal innate immune pathway involving the bitter taste G protein-coupled receptor T2R38. This case–control study aimed to evaluate whether T2R38 polymorphisms associate with the buccal microbial composition in RA. Genomic DNA was obtained from buccal swabs of 35 RA patients and 64 non-RA controls. TAS2R38 genotypes were determined by Sanger sequencing. The buccal microbiome was assessed by Illumina MiSeq sequencing of the V4-16S rRNA gene. Bacterial community differences were analyzed with alpha and beta diversity measures. Linear discriminant analysis effect size identified taxa discriminating between RA versus non-RA and across TAS2R38 genotypes. TAS2R38 genotype frequency was similar between RA and non-RA controls (PAV/PAV; PAV/AVI; AVI/AVI: RA 42.9%; 45.7%; 11.4% versus controls 32.8%; 48.4%; 18.8%, chi-square (2, N = 99) = 2.1, p = 0.35). The relative abundance of Porphyromonas, among others, differed between RA and non-RA controls. The relative abundance of several bacterial species also differed across TAS2R38 genotypes. These findings suggest an association between T2R38 polymorphisms and RA buccal microbial composition. However, further research is needed to understand the impact of T2R38 in oral health and RA development.     

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein‐coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6‐trnE‐Cytb‐D‐loop‐trnL2 (CUN)‐ND5‐D‐loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.     

Genetic polymorphisms in the glutamate-rich protein of Plasmodium falciparum field isolates from a malaria-endemic area of Brazil

The genetic diversity displayed by Plasmodium falciparum, the most deadly Plasmodium species, is a significant obstacle for effective malaria vaccine development. In this study, we identified genetic polymorphisms in P. falciparum glutamate-rich protein (GLURP), which is currently being tested in clinical trials as a malaria vaccine candidate, from isolates found circulating in the Brazilian Amazon at variable transmission levels. The study was performed using samples collected in 1993 and 2008 from rural villages situated near Porto Velho, in the state of Rondônia. DNA was extracted from 126 P. falciparum-positive thick blood smears using the phenol-chloroform method and subjected to a nested polymerase chain reaction protocol with specific primers against two immunodominant regions of GLURP, R0 and R2. Only one R0 fragment and four variants of the R2 fragment were detected. No differences were observed between the two time points with regard to the frequencies of the fragment variants. Mixed infections were uncommon. Our results demonstrate conservation of GLURP-R0 and limited polymorphic variation of GLURP-R2 in P. falciparum isolates from individuals living in Porto Velho. This is an important finding, as genetic polymorphisms in B and T-cell epitopes could have implications for the immunological properties of the antigen.     

Extensive Chromosome Homoeology among Brassiceae Species Were Revealed by Comparative Genetic Mapping with High-Density EST-Based SNP Markers in Radish (Raphanus sativus L.)

A linkage map of expressed sequence tag (EST)-based markers in radish (Raphanus sativus L.) was constructed using a low-cost and high-efficiency single-nucleotide polymorphism (SNP) genotyping method named multiplex polymerase chain reaction–mixed probe dot-blot analysis developed in this study. Seven hundred and forty-six SNP markers derived from EST sequences of R. sativus were assigned to nine linkage groups with a total length of 806.7 cM. By BLASTN, 726 markers were found to have homologous genes in Arabidopsis thaliana, and 72 syntenic regions, which have great potential for utilizing genomic information of the model species A. thaliana in basic and applied genetics of R. sativus, were identified. By construction and analysis of the genome structures of R. sativus based on the 24 genomic blocks within the Brassicaceae ancestral karyotype, 23 of the 24 genomic blocks were detected in the genome of R. sativus, and half of them were found to be triplicated. Comparison of the genome structure of R. sativus with those of the A, B, and C genomes of Brassica species and that of Sinapis alba L. revealed extensive chromosome homoeology among Brassiceae species, which would facilitate transfer of the genomic information from one Brassiceae species to another.     

Intergenomic single nucleotide polymorphisms as a tool for bacterial artificial chromosome contig building of homoeologous Brassica napus regions

Background

Homoeologous sequences pose a particular challenge if bacterial artificial chromosome (BAC) contigs shall be established for specific regions of an allopolyploid genome. Single nucleotide polymorphisms (SNPs) differentiating between homoeologous genomes (intergenomic SNPs) may represent a suitable screening tool for such purposes, since they do not only identify homoeologous sequences but also differentiate between them.

Results

Sequence alignments between Brassica rapa (AA) and Brassica oleracea (CC) sequences mapping to corresponding regions on chromosomes A1 and C1, respectively were used to identify single nucleotide polymorphisms between the A and C genomes. A large fraction of these polymorphisms was also present in Brassica napus (AACC), an allopolyploid species that originated from hybridisation of A and C genome species. Intergenomic SNPs mapping throughout homoeologous chromosome segments spanning approximately one Mbp each were included in Illumina’s GoldenGate® Genotyping Assay and used to screen multidimensional pools of a Brassica napus bacterial artificial chromosome library with tenfold genome coverage. Based on the results of 50 SNP assays, a BAC contig for the Brassica napus A subgenome was established that spanned the entire region of interest. The C subgenome region was represented in three BAC contigs.

Conclusions

This proof-of-concept study shows that sequence resources of diploid progenitor genomes can be used to deduce intergenomic SNPs suitable for multiplex polymerase chain reaction (PCR)-based screening of multidimensional BAC pools of a polyploid organism. Owing to their high abundance and ease of identification, intergenomic SNPs represent a versatile tool to establish BAC contigs for homoeologous regions of a polyploid genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-560) contains supplementary material, which is available to authorized users.     

Experimental infection of the rabbit tick,<Emphasis Type="Italic"> Haemaphysalis leporispalustris</Emphasis>, with the bacterium<Emphasis Type="Italic"> Rickettsia rickettsii</Emphasis>, and comparative biology of infected and uninfected tick lineages

The present study consisted of two experiments that evaluated experimental infections of Haemaphysalis leporispalustris ticks by a Brazilian strain of Rickettsia rickettsii, and their effect on tick biology. In experiment I, ticks were exposed to R. rickettsii during the larval, nymphal or adult stages by feeding on rabbits (Oryctolagus cuniculus) needle-inoculated with R. rickettsii, and thereafter reared on uninfected rabbits for the entire next tick generation. Regardless of the tick stage that acquired the infection, all subsequent tick stages were shown to be infected by PCR (infection rates varying from 1.3 to 41.7%), and were able to transmit R. rickettsii to uninfected rabbits, as demonstrated by rabbit seroconversion, guinea pig inoculation with rabbit blood, and PCR on rabbit blood. In Experiment II, ticks were exposed to R. rickettsii during the larval stage by feeding on rabbits co-infested with R. rickettsii-infected adult ticks, and thereafter reared on uninfected rabbits until the next generation of larvae. Again, all subsequent tick stages were shown to be infected by PCR (infection rates varying from 3.0 to 40.0%), and were able to transmit R. rickettsii to uninfected rabbits. Thus, it was demonstrated that larvae, nymphs, and adults of H. leporispalustris were able to acquire and maintain the R. rickettsii infection by transstadial and transovarial transmissions within the tick population, with active transmission of the bacterium to susceptible rabbits by all parasitic stages. Analyses of biological parameters of uninfected and R. rickettsii-infected tick lineages were performed in order to evaluate possible deleterious effects of R. rickettsii to the infected tick lineages. Surprisingly, all but one of the four R. rickettsii-experimental groups of the present study showed overall better biological performance than their sibling uninfected control ticks. Results of the present study showed that H. leporispalustris could support infection by a high virulent strain of R. rickettsii for at least two generations, in which infected tick lineages tended to have better performance than uninfected ticks. Our results support a possible role of H. leporispalustris in the enzootic maintenance of R. rickettsii in Latin America, as previously suggested by earlier works.