SVARAP and aSVARAP: simple tools for quantitative analysis of nucleotide and amino acid variability and primer selection for clinical microbiology

Contribution of RNA-RNA Interactions Mediated by the Genome Packaging Signals for the Selective Genome Packaging of Influenza A Virus

The influenza A virus genome is composed of eight single-stranded negative-sense viral RNA segments (vRNAs). The eight vRNAs are selectively packaged into each progeny virion. This process likely involves specific interactions between the vRNAs via segment-specific packaging signals located in both the 3'- and 5'-terminal regions of the respective vRNAs. To assess the importance of vRNA-vRNA interactions via packaging signals for selective genome packaging, we generated mutant viruses possessing silent mutations in the packaging signal region of the hemagglutinin (HA) vRNA. A mutant virus possessing silent mutations in nucleotides (nt) 1664 to 1676 resulted in defects in HA vRNA incorporation and showed a reduction in viral growth. After serial passage, the mutant virus acquired additional mutations in the 5'-terminal packaging signal regions of both the HA and polymerase basic 2 (PB2) vRNAs. These mutations contributed to the recovery of viral growth and HA vRNA packaging efficiency. In addition, an RNA-RNA interaction between the 5' ends of HA and PB2 vRNAs was confirmed in vitro, and this interaction was disrupted following the introduction of silent mutations in the HA vRNA. Thus, our results demonstrated that RNA-RNA interactions between the packaging signal regions of HA vRNA and PB2 vRNA are important for selective genome packaging. IMPORTANCE While numerous viral genomes comprise a single genome segment, the influenza A virus possesses eight segmented genomes. Influenza A virus can benefit from having a segmented genome because the segments can reassort with other strains of the influenza virus to create new genetically distinct strains. The influenza A virus efficiently incorporates one copy of each of its eight genomic segments per viral particle. However, the mechanism by which each segment is specifically selected is poorly understood. The genome segments contain RNA signals that facilitate the incorporation of segments into virus particles. These regions may facilitate specific interactions between the genome segments, creating an eight-segment complex, which can then be packaged into individual particles. In this study, we provide evidence that RNA signals contribute to specific interactions between two of the influenza virus genome segments.

Fluorescent Hybridization of Mycobacterium leprae in Skin Samples Collected in Burkina Faso

Leprosy is caused by Mycobacterium leprae, and it remains underdiagnosed in Burkina Faso. We investigated the use of fluorescent in situ hybridization (FISH) for detecting M. leprae in 27 skin samples (skin biopsy samples, slit skin samples, and skin lesion swabs) collected from 21 patients from Burkina Faso and three from Côte d'Ivoire who were suspected of having cutaneous leprosy. In all seven Ziehl-Neelsen-positive skin samples (four skin biopsy samples and three skin swabs collected from the same patient), FISH specifically identified M. leprae, including one FISH-positive skin biopsy sample that remained negative after testing with PCR targeting the rpoB gene and with the GenoType LepraeDR assay. Twenty other skin samples and three negative controls all remained negative for Ziehl-Neelsen staining, FISH, and rpoB PCR. These data indicate the usefulness of a microscopic examination of skin samples after FISH for first-line diagnosis of cutaneous leprosy. Accordingly, FISH represents a potentially useful point-of-care test for the diagnosis of cutaneous leprosy.

Core gene-based molecular detection and identification of Acanthamoeba species

Acanthamoeba spp. are predominant free-living amoebae of water and soil. They have been used as tools for the isolation and culture of microbes that resist after their phagocytosis, such as Legionella-like bacteria, and, more recently giant viruses for which differences in permissiveness have been reported. However, problems have been reported regarding their identification at the species level. The present work implemented specific PCR systems for the detection and identification of Acanthamoeba species through comparison of sequences and phylogenetic analyses. Thirty-three Acanthamoeba isolates were studied, including 20 reference strains and 13 isolates retrieved from water, soil or clinical samples. Previous delineation of a core genome encompassing 826 genes based on draft genome sequences from 14 Acanthamoeba species allowed designing PCR systems for one of these core genes that encodes an alanine-tRNA ligase. These primers allowed an efficient and specific screening to detect Acanthamoeba presence. In addition, they identified all 20 reference strains, while partial and complete sequences coding for 18S ribosomal RNA identified only 11 (55%). We found that four isolates may be considered as new Acanthamoeba species. Consistent with previous studies, we demonstrated that some Acanthamoeba isolates were incorrectly assigned to species using the 18S rDNA sequences. Our implemented tool may help determining which Acanthamoeba strains are the most efficient for the isolation of associated microorganisms.

Genome-wide analysis of influenza viral RNA and nucleoprotein association

Influenza A virus (IAV) genomes are composed of eight single-stranded RNA segments that are coated by viral nucleoprotein (NP) molecules. Classically, the interaction between NP and viral RNA (vRNA) is depicted as a uniform pattern of ‘beads on a string’. Using high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP), we identified the vRNA binding profiles of NP for two H1N1 IAV strains in virions. Contrary to the prevailing model for vRNA packaging, NP does not bind vRNA uniformly in the A/WSN/1933 and A/California/07/2009 strains, but instead each vRNA segment exhibits a unique binding profile, containing sites that are enriched or poor in NP association. Intriguingly, both H1N1 strains have similar yet distinct NP binding profiles despite extensive sequence conservation. Peaks identified by HITS-CLIP were verified as true NP binding sites based on insensitivity to DNA antisense oligonucleotide-mediated RNase H digestion. Moreover, nucleotide content analysis of NP peaks revealed that these sites are relatively G-rich and U-poor compared to the genome-wide nucleotide content, indicating an as-yet unidentified sequence bias for NP association in vivo. Taken together, our genome-wide study of NP–vRNA interaction has implications for the understanding of influenza vRNA architecture and genome packaging.

Hepatitis C virus of subtype 2l in Marseille, southeastern France

The rate of eradication of chronic hepatitis C considerably increases with direct-acting antiviral agents, particularly hepatitis C virus (HCV) polymerase inhibitors. While implementing full-length HCV NS5B polymerase sequencing in our clinical microbiology laboratory, we identified atypical HCV sequences, classified as subtype 2l, from 2 patients. HCV-2l NS5B polymerase sequences were detected from 5 and 14 additional patients by screening our laboratory hepatitis virus sequence database and the NCBI GenBank sequence database. Phylogenetic analyses show unambiguously that all HCV-2l sequences are clustered apart from HCV 2 non-l sequences, which compose a second cluster. Mean (±SD) nucleotide identity between near full-length NS5B fragments of subtype 2l was 93.4 ± 0.8% (range: 92.4-95.1). Of note, all HCV-2l sequences obtained in our laboratory and in other centers were from serum samples collected in France. Analysis of the HCV-2l NS5B polymerase amino acid sequences at 30 positions critical for interaction with or resistance to HCV polymerase inhibitors showed specific patterns.

Serotype-associated polymorphisms in a partial rpoB gene sequence of Salmonella enterica

Salmonella enterica is a zoonotic bacterium with more than 2500 serotypes, which affect a wide range of hosts and produce diverse clinical outcomes. Strain identification usually involves costly and time-demanding procedures. This paper describes the sequencing of a rpoB hypervariable gene segment (847 bp) that allows identification of serotypes in S. enterica strains isolated from several hosts. The nucleotide similarity values among S. enterica serotypes ranged from 98.23% to 99.88%, with potential usefulness for devising a simple one-step sequencing as a first approach for identification of S. enterica strains. In conclusion, the analysis of polymorphisms in the partial rpoB sequence can discriminate S. enterica strains at the subspecies level.

Real-time PCR systems targeting giant viruses of amoebae and their virophages

Giant viruses that infect amoebae, including mimiviruses and marseilleviruses, were first described in 2003. Virophages were subsequently described that infect mimiviruses. Culture isolation with Acanthamoeba spp. and metagenomic studies have shown that these giant viruses are common inhabitants of our biosphere and have enabled the recent detection of these viruses in human samples. However, the genomes of these viruses display substantial genetic diversity, making it a challenge to examine their presence in environmental and clinical samples using conventional and real-time PCR. We designed and evaluated the performance of PCR systems capable of detecting all currently isolated mimiviruses, marseilleviruses and virophages to assess their prevalence in various samples. Our real-time PCR assays accurately detected all or most of the members of the currently delineated lineages of giant viruses infecting acanthamoebae as well as the mimivirus virophages, and enabled accurate classification of the mimiviruses of amoebae in lineages A, B or C. We were able to detect four new mimiviruses directly from environmental samples and correctly classified these viruses within mimivirus lineage C. This was subsequently confirmed by culture on amoebae followed by partial Sanger sequencing. PCR systems such as those implemented here may contribute to an improved understanding of the prevalence of mimiviruses, their virophages and marseilleviruses in humans.

rpoB gene as a novel molecular marker to infer phylogeny in Planctomycetales

The 16S rRNA gene has been used in the last decades as a gold standard for determining the phylogenetic position of bacteria and their taxonomy. It is a well conserved gene, with some variations, present in all bacteria and allows the reconstruction of genealogies of microorganisms. Nevertheless, this gene has its limitations when inferring phylogenetic relationships between closely related isolates. To overcome this problem, DNA-DNA hybridization appeared as a solution to clarify interspecies relationships when the sequence similarity of the 16S rRNA gene is above 97 %. However, this technique is time consuming, expensive and laborious and so, researchers developed other molecular markers such as sequencing of housekeeping or functional genes for accurate determination of bacterial phylogeny. One of these genes that have been used successfully, particularly in clinical microbiology, codes for the beta subunit of the RNA polymerase (rpoB). The rpoB gene is sufficiently conserved to be used as a molecular clock, it is present in all bacteria and it is a mono-copy gene. In this study, rpoB gene sequencing was applied to the phylum Planctomycetes. Based on the genomes of 19 planctomycetes it was possible to determine the correlation between the rpoB gene sequence and the phylogenetic position of the organisms at a 95-96 % sequence similarity threshold for a novel species. A 1200-bp fragment of the rpoB gene was amplified from several new planctomycetal isolates and their intra and inter-species relationships to other members of this group were determined based on a 96.3 % species border and 98.2 % for intraspecies resolution.

The prevalence of chronic fatigue syndrome/ myalgic encephalomyelitis: a meta-analysis

Purpose

To perform a meta-analysis to examine variability among prevalence estimates for CFS/ME, according to the method of assessment used.

Methods

Databases were systematically searched for studies on CFS/ME prevalence in adults that applied the 1994 Centers for Disease Control (CDC) case definition.1 Estimates were categorized into two methods of assessment: self-reporting of symptoms versus clinical assessment of symptoms. Meta-analysis was performed to pool prevalences by assessment using random effects modeling. This was stratified by sample setting (community or primary care) and heterogeneity was examined using the I² statistic.

Results

Of 216 records found, 14 studies were considered suitable for inclusion. The pooled prevalence for self-reporting assessment was 3.28% (95% CI: 2.24–4.33) and 0.76% (95% CI: 0.23–1.29) for clinical assessment. High variability was observed among self-reported estimates, while clinically assessed estimates showed greater consistency.

Conclusion

The observed heterogeneity in CFS/ME prevalence may be due to differences in method of assessment. Stakeholders should be cautious of prevalence determined by the self-reporting of symptoms alone. The 1994 CDC case definition appeared to be the most reliable clinical assessment tool available at the time of these studies. Improving clinical case definitions and their adoption internationally will enable better comparisons of findings and inform health systems about the true burden of CFS/ME.

Gemi: PCR primers prediction from multiple alignments

Designing primers and probes for polymerase chain reaction (PCR) is a preliminary and critical step that requires the identification of highly conserved regions in a given set of sequences. This task can be challenging if the targeted sequences display a high level of diversity, as frequently encountered in microbiologic studies. We developed Gemi, an automated, fast, and easy-to-use bioinformatics tool with a user-friendly interface to design primers and probes based on multiple aligned sequences. This tool can be used for the purpose of real-time and conventional PCR and can deal efficiently with large sets of sequences of a large size.

Evolutionary pattern of 5'-UTR of enteroviruses and primer update for the detection of enteroviral RNA in environmental samples

OBJECTIVE

To study the recombination events among enterovirus strains and the development of specific primers for the detection of enteroviruses in environmental samples.

METHODS

Nucleotide sequence analysis of enteroviruses deposited in the international database GenBank (www.ncbi.nlm.nih.gov/Genbank) was conducted to develop specific primers for the detection of these viruses. The specificity and sensitivity of the method were tested using coxackievirus B3 strain Nancy, environmental isolate of human hepatitis A virus and human rotavirus strain WA. Seventy sewage samples were analyzed.

RESULTS

Enterovirus genome was detected in all positive samples. The genome of enterovirus was not detected in negative samples. The level of detection of these viruses was 10(2) TCID(50)/mL.

CONCLUSIONS

The development of new primers is an important issue for the detection of enteroviruses in the environment and the assessment of risk factors to human health.

Core gene set as the basis of multilocus sequence analysis of the subclass Actinobacteridae

Comparative genomic sequencing is shedding new light on bacterial identification, taxonomy and phylogeny. An in silico assessment of a core gene set necessary for cellular functioning was made to determine a consensus set of genes that would be useful for the identification, taxonomy and phylogeny of the species belonging to the subclass Actinobacteridae which contained two orders Actinomycetales and Bifidobacteriales. The subclass Actinobacteridae comprised about 85% of the actinobacteria families. The following recommended criteria were used to establish a comprehensive gene set; the gene should (i) be long enough to contain phylogenetically useful information, (ii) not be subject to horizontal gene transfer, (iii) be a single copy (iv) have at least two regions sufficiently conserved that allow the design of amplification and sequencing primers and (v) predict whole-genome relationships. We applied these constraints to 50 different Actinobacteridae genomes and made 1,224 pairwise comparisons of the genome conserved regions and gene fragments obtained by using Sequence VARiability Analysis Program (SVARAP), which allow designing the primers. Following a comparative statistical modeling phase, 3 gene fragments were selected, ychF, rpoB, and secY with R²>0.85. Selected sets of broad range primers were tested from the 3 gene fragments and were demonstrated to be useful for amplification and sequencing of 25 species belonging to 9 genera of Actinobacteridae. The intraspecies similarities were 96.3–100% for ychF, 97.8–100% for rpoB and 96.9–100% for secY among 73 strains belonging to 15 species of the subclass Actinobacteridae compare to 99.4–100% for 16S rRNA. The phylogenetic topology obtained from the combined datasets ychF+rpoB+secY was globally similar to that inferred from the 16S rRNA but with higher confidence. It was concluded that multi-locus sequence analysis using core gene set might represent the first consensus and valid approach for investigating the bacterial identification, phylogeny and taxonomy.

Pepper mild mottle virus, a plant virus associated with specific immune responses, Fever, abdominal pains, and pruritus in humans

BACKGROUND

Recently, metagenomic studies have identified viable Pepper mild mottle virus (PMMoV), a plant virus, in the stool of healthy subjects. However, its source and role as pathogen have not been determined.

METHODS AND FINDINGS

21 commercialized food products containing peppers, 357 stool samples from 304 adults and 208 stool samples from 137 children were tested for PMMoV using real-time PCR, sequencing, and electron microscopy. Anti-PMMoV IgM antibody testing was concurrently performed. A case-control study tested the association of biological and clinical symptoms with the presence of PMMoV in the stool. Twelve (57%) food products were positive for PMMoV RNA sequencing. Stool samples from twenty-two (7.2%) adults and one child (0.7%) were positive for PMMoV by real-time PCR. Positive cases were significantly more likely to have been sampled in Dermatology Units (p<10(-6)), to be seropositive for anti-PMMoV IgM antibodies (p = 0.026) and to be patients who exhibited fever, abdominal pains, and pruritus (p = 0.045, 0.038 and 0.046, respectively).

CONCLUSIONS

Our study identified a local source of PMMoV and linked the presence of PMMoV RNA in stool with a specific immune response and clinical symptoms. Although clinical symptoms may be imputable to another cofactor, including spicy food, our data suggest the possibility of a direct or indirect pathogenic role of plant viruses in humans.

rpoB sequence-based identification of Mycobacterium avium complex species

The Mycobacterium avium complex (MAC) comprises slowly growing mycobacteria responsible for opportunistic infections and zoonoses. The ability to speciate MAC isolates in the clinical microbiology laboratory is critical for determining the organism implicated in clinical disease and for epidemiological investigation of the source of infection. Investigation of a 711 bp variable fragment of rpoB flanked by the Myco-F/Myco-R primers found a 0.7-5.1 % divergence among MAC reference strains, with Mycobacterium chimaera and Mycobacterium intracellulare being the most closely related. Using a 0.7 % divergence cut-off, 83 % of 100 clinical isolates, which had been previously identified by phenotypic characteristics and 16S-23S rDNA intergenic spacer (ITS) probing, were identified as M. avium, 8 % as M. intracellulare and 2 % as M. chimaera. The uniqueness of seven isolates, exhibiting < 99.3 % rpoB sequence similarity with MAC reference strains, was confirmed by 16S rDNA, ITS and hsp65 sequencing and phylogenetic analyses. Partial rpoB gene sequencing using the Myco-F/Myco-R primers permits one-step identification of MAC isolates at the species level and the detection of potentially novel MAC species.

The rpoB gene as a tool for clinical microbiologists

The rpoB gene, encoding the beta-subunit of RNA polymerase, has emerged as a core gene candidate for phylogenetic analyses and identification of bacteria, especially when studying closely related isolates. Together with the 16S rRNA gene, rpoB has helped to delineate new bacterial species and refine bacterial community analysis, as well as enabling the monitoring of rifampicin resistance-conferring mutations. Sequencing of rpoB enables efficient estimation of bacterial G+C% content, DNA-DNA hybridization value and average nucleotide identity (percentage of the total genomic sequence shared between two strains) when taxonomic relationships have been firmly established. New identification tools targeting a rpoB gene fragment located between positions 2300 and 3300 have been developed recently. Therefore, inclusion of the rpoB gene sequence would be useful when describing new bacterial species.

Amino acid conservation in the gp41 transmembrane protein and natural polymorphisms associated with enfuvirtide resistance across HIV-1 variants

Information about gp41 variability across distinct HIV-1 subtypes is scarce, and yet such knowledge would be desirable for designing new drugs targeting this viral protein. Conserved gp41 residues in viruses derived from 79 individuals infected with distinct HIV-1 subtypes (29 A, 25 B, 8 C, 3 D, 4 F, 4 G, 2 H, 1 J, 1 U, and 2 CRF06_cpx) and naive for entry inhibitors were examined. Conservation of gp41 was also examined in 908, 56, and 3 HIV-1 group M, O, and N sequences, respectively, available at the Los Alamos HIV Sequence Database. Among the 345 residues in the full gp41 protein, 36% showed up to 90% conservation in all 987 group M sequences, as did 40% of 56 group O sequences and 49% of 3 group N sequences. The HR1 region (residues 29-82) showed a higher proportion of highly conserved residues than did the HR2 region (residues 116-161) in all groups (65 vs. 34% in group M, 57 vs. 46% in group O, and 80 vs. 52% in group N). Some secondary resistance mutations to enfuvirtide were found as natural polymorphisms (A30V and Q56K/R in group M, Q56R and S138A in group O, and S138A in group N). In fact, A30V was a signature change in clade G and CRF06_cpx, whereas Q56K/R was a signature change for clades A and J, as well as for CRF04_cpx, CRF09_cpx, CRF11_cpx, and CRF13_cpx. The relative conservation of amino acids in gp41 across HIV-1 variants indirectly highlights the critical role of this protein for HIV infectivity and makes it feasible to design new entry inhibitors with activity against diverse HIV-1 variants.

Polyphenol oxidase (PPO) in wheat and wild relatives: molecular evidence for a multigene family

Wheat polyphenol oxidase (PPO) is the major cause of browning reactions that discolor Asian noodles and other wheat products. It has been hypothesized that genes encoding wheat PPOs may have evolved by gene duplication into a multigene family. Here we characterized PPO genomic sequences from diploid (Triticum monococcum, T. urartu, Aegilops tauschii, and Ae. speltoides), tetraploid (T. turgidum, subspecies dicoccoides and durum) and hexaploid (T. aestivum cultivars Klasic and ID377s) wheat species to gain a better understanding of the structure and organization of PPO genes. DNA fragments were amplified from a highly polymorphic and phylogenetic informative region of the gene. As a result, we obtained highly discriminative sequences. Three distinct PPOs, obtained from the A genome of T. monococcum, provided evidence for gene duplication events (paralogous loci). Furthermore, the number of sequences obtained for bread and durum wheat was higher than the expected number of orthologous loci. Sequence comparison revealed nucleotide and structural diversity, and detected five sequence intron types, all with a common insertion position. This was hypothesized to be homologous to that of intron 2 of previously reported wheat PPOs. A MITE of the Stowaway family accounted for the major difference between the five intervening sequences, and was unique to T. aestivum cv. Klasic. Nucleotide and structural diversity, together with well-resolved phylogenetic trees, provided molecular evidence to support the hypothesis of a PPO multigene family structure and organization.

Polymorphism and drug selected mutations of reverse transcriptase gene in 102 HIV-1 infected patients living in China

Few data are available for genotypic patterns within human immunodeficiency virus-1 (HIV-1) reverse transcriptase (RT) in drug-naive patients and RT inhibitor (RTI) treated patients in China. This study aimed at characterizing the polymorphism of RT HIV-1 in the absence of drug treatment and to identify known and unknown mutations emerging under RTI selective pressure. The HIV-1 RT gene from 21 drug-naive patients and 81 RTI treated patients from three provinces in China was analyzed. Most patients (>80%) received a triple regimen including stavudine (d4T) plus didanosine (ddI) and nevirapine (NVP), or d4T plus lamivudine (3TC) and efavirenz (EFV), or zidovudine (AZT) +ddI + NVP. In untreated patients, four highly polymorphic positions were found (122, 200, 207, and 211). In treated patients, two patterns of resistance associated mutations (RAMs) were observed: (1) K65R (9.8%), L74V (7.4%), M184V (7.4%), Q151M (5%), and thymidine analogue mutations (TAMs) (9.3%) including T215Y (5.5%), in patients who underwent ddI + d4T + NVP. (2) T215Y (23%), M184V (20%), and TAMs (15.4%) in patients receiving d4T + 3TC + EFV. In all cases, a high prevalence of non-nucleoside RTIs (NNRTI) RAMs (41.9%) was found. Four RTI suspected new RAMs were described at position 142, 221, 224, and 228. An association between H221Y and L228H/R with Y181C was noted. These data highlight the predominant spread of NNRTI RAM in China, depict the specific genotypic pattern of RTI selected mutations in China, and suggest the association of newly described mutations with RTI therapy.

PartialrpoBgene sequencing for identification ofLeptospiraspecies

The usual target for sequence-based identification of Leptospira species is the 16S rRNA gene. However, because the 16S rRNA gene is not polymorphic enough, it is necessary to sequence a 1500 bp segment of this gene for accurate identification. Based on the alignment of previously determined rpoB of three Leptospira strains, we designed and tested a primer pair that enabled us to amplify and sequence a 600 bp segment of Leptospira rpoB. This segment was species-specific for the 16 species tested, but was unable to separate Leptospira interrogans serovars accurately. For the 11 L. interrogans serovars tested, only seven genotypes could be determined. We thus think that analysis of partial rpoB may be useful as an initial screening test for the identification of a new isolate of Leptospira and detection or identification of Leptospira in clinical or environmental samples, but not for serovar determination.

Sca1, a previously undescribed paralog from autotransporter protein-encoding genes in Rickettsia species

BACKGROUND

Among the 17 genes encoding autotransporter proteins of the "surface cell antigen" (sca) family in the currently sequenced Rickettsia genomes, ompA, sca5 (ompB) and sca4 (gene D), have been extensively used for identification and phylogenetic purposes for Rickettsia species. However, none of these genes is present in all 20 currently validated Rickettsia species. Of the remaining 14 sca genes, sca1 is the only gene to be present in all nine sequenced Rickettsia genomes. To estimate whether the sca1 gene is present in all Rickettsia species and its usefulness as an identification and phylogenetic tool, we searched for sca1genes in the four published Rickettsia genomes and amplified and sequenced this gene in the remaining 16 validated Rickettsia species.

RESULTS

Sca1 is the only one of the 17 rickettsial sca genes present in all 20 Rickettsia species. R. prowazekii and R. canadensis exhibit a split sca1 gene whereas the remaining species have a complete gene. Within the sca1 gene, we identified a 488-bp variable sequence fragment that can be amplified using a pair of conserved primers. Sequences of this fragment are specific for each Rickettsia species. The phylogenetic organization of Rickettsia species inferred from the comparison of sca1 sequences strengthens the classification based on the housekeeping gene gltA and is similar to those obtained from the analyses of ompA, sca5 and sca4, thus suggesting similar evolutionary constraints. We also observed that Sca1 protein sequences have evolved under a dual selection pressure: with the exception of typhus group rickettsiae, the amino-terminal part of the protein that encompasses the predicted passenger domain, has evolved under positive selection in rickettsiae. This suggests that the Sca1 protein interacts with the host. In contrast, the C-terminal portion containing the autotransporter domain has evolved under purifying selection. In addition, sca1 is transcribed in R. conorii, and might therefore be functional in this species.

CONCLUSION

The sca1 gene, encoding an autotransporter protein that evolves under dual evolution pressure, is the only sca-family gene to be conserved by all Rickettsia species. As such, it is a valuable identification target for these bacteria, especially because rickettsial isolates can be identified by amplification and sequencing of a discriminatory gene fragment using a single primer pair. It may also be used as a phylogenetic tool. However, its current functional status remains to be determined although it was found expressed in R. conorii.