Gene conversion: a mechanism for generation of heterogeneity in the tprK gene of Treponema pallidum during infection

The tprK gene sequence of Treponema pallidum subspecies pallidum (T. pallidum) is heterogeneous within and among isolates. Heterogeneity in the tprK open reading frame is localized in seven discrete variable (V) regions, and variability results from apparent base changes, insertions or deletions. The TprK V regions are the focus of anti-TprK antibodies arising during infection. To test our hypothesis that V region sequences change during infection and passage, we developed a clonal isolate from the Chicago strain of T. pallidum and confirmed V region diversification during passage of this isolate. We describe the sequence anatomy of the seven V regions of tprK and the identification of putative donor sites for new V region sequences, and we propose a model for generation of new V regions by segmental gene conversion. These findings suggest that antigenic variation of TprK occurs in T. pallidum and may be important in immune evasion and persistence.     

Comparative investigation of the genomic regions involved in antigenic variation of the TprK antigen among treponemal species, subspecies, and strains

Although the three Treponema pallidum subspecies (T. pallidum subsp. pallidum, T. pallidum subsp. pertenue, and T. pallidum subsp. endemicum), Treponema paraluiscuniculi, and the unclassified Fribourg-Blanc treponeme cause clinically distinct diseases, these pathogens are genetically and antigenically highly related and are able to cause persistent infection. Recent evidence suggests that the putative surface-exposed variable antigen TprK plays an important role in both treponemal immune evasion and persistence. tprK heterogeneity is generated by nonreciprocal gene conversion between the tprK expression site and donor sites. Although each of the above-mentioned species and subspecies has a functional tprK antigenic variation system, it is still unclear why the level of expression and the rate at which tprK diversifies during infection can differ significantly among isolates. To identify genomic differences that might affect the generation and expression of TprK variants among these pathogens, we performed comparative sequence analysis of the donor sites, as well as the tprK expression sites, among eight T. pallidum subsp. pallidum isolates (Nichols Gen, Nichols Sea, Chicago, Sea81-4, Dal-1, Street14, UW104, and UW126), three T. pallidum subsp. pertenue isolates (Gauthier, CDC2, and Samoa D), one T. pallidum subsp. endemicum isolate (Iraq B), the unclassified Fribourg-Blanc isolate, and the Cuniculi A strain of T. paraluiscuniculi. Synteny and sequence conservation, as well as deletions and insertions, were found in the regions harboring the donor sites. These data suggest that the tprK recombination system is harbored within dynamic genomic regions and that genomic differences might be an important key to explain discrepancies in generation and expression of tprK variants among these Treponema isolates.     

Subfamily I Treponema pallidum repeat protein family: sequence variation and immunity

A 12-membered Treponema pallidum repeat (Tpr) protein family has been identified in T. pallidum subsp. pallidum, the causative agent of syphilis. The subfamily I Tpr proteins (C, D, F, and I) possess conserved sequence at the N- and C-termini and central regions that differentiate the members. These proteins may be important in the immune response during syphilis infection and in protective immunity. Strong antibody responses have been observed toward some of the subfamily I Tpr proteins during infection with different syphilis isolates. Some sequence variation has also been identified in one subfamily I Tpr member, TprD, among T. pallidum subsp. pallidum isolates. In this study, we examined sequences in the remaining subfamily I Tpr proteins among strains. Both TprF and TprI were conserved among T. pallidum subsp. pallidum isolates.While some heterogeneity was identified in TprC. We further examined the immune response and protective capacity of TprF protein in this paper. We demonstrate that the N-terminal conserved region of the subfamily I Tpr proteins elicits strong antibody and T-cell responses during infection, and immunization with this region attenuates syphilitic lesion development upon infectious challenge.     

Segregation of B and T cell epitopes of Treponema pallidum repeat protein K to variable and conserved regions during experimental syphilis infection

Robust immune responses clear millions of treponemes to resolve lesions of primary and secondary syphilis, but cannot clear the treponemes that lead to debilitating and sometimes fatal tertiary syphilis. It is also known that the rabbit model and humans can be reinfected with heterologous isolates. How some treponemes are able to escape the immune system is unknown. In our laboratories rabbits immunized with the Seattle Nichols strain Treponema pallidum repeat protein K (TprK) were previously shown to have attenuated lesion development following challenge. In other isolates, TprK was shown to have seven discrete variable regions, with sequence variation among and within isolates. Using overlapping synthetic 20-aa peptides, we demonstrate that during experimental infection with the Nichols strain, the T cell responses are directed to conserved regions, while the Ab responses are directed primarily to variable regions. Abs from rabbits immunized with recombinant TprK recognized conserved and variable regions, suggesting that the conserved regions are inherently as immunogenic as the variable regions. TprK variability may allow some treponemes to escape recognition from Abs. The variable region heterogeneity may help explain the lack of protection against heterologous isolates.     

Multiple alleles of Treponema pallidum repeat gene D in Treponema pallidum isolates

Two new tprD alleles have been identified in Treponema pallidum: tprD2 is found in 7 of 12 T. pallidum subsp. pallidum isolates and 7 of 8 non-pallidum isolates, and tprD3 is found in one T. pallidum subsp. pertenue isolate. Antibodies against TprD2 are found in persons with syphilis, demonstrating that tprD2 is expressed during infection.     

Characterization of the serum requirement for macrophage-mediated killing of Treponema pallidum ssp. pallidum: Relationship to the development of opsonizing antibodies

Abstract Serum pools were collected from rabbits bled at various times after intra-testicular infection with Treponema pallidum ssp. pallidum . These were tested for their ability to opsonize T. pallidum and promote killing of the organisms by macrophages. Compared to normal sera, significant opsonization was first seen on day 10 of infection as measured by both ingestion ( P < 0.001) and macrophage-mediated killing ( P = 0.006); significant levels of functional antibodies persisted through 300 days of infection. Although opsonic activity peaked early in infection, antibodies that promoted optimal macrophage-mediated killing developed much later, suggesting that these two functions may represent activities of antibodies with differing specificities or affinities. The initial development of antibodies that augment both phagocytosis and killing corresponds with the in vivo clearance of treponemes from the primary site of infection. These observations support the hypothesis that macrophages are the major effector mechanism for elimination of T. pallidum during early syphilis infection.     

Longitudinal TprK profiling of in vivo and in vitro-propagated Treponema pallidum subsp. pallidum reveals accumulation of antigenic variants in absence of immune pressure

Immune evasion by Treponema pallidum subspecies pallidum (T. pallidum) has been attributed to antigenic variation of its putative outer-membrane protein TprK. In TprK, amino acid diversity is confined to seven variable (V) regions, and generation of sequence diversity within the V regions occurs via a non-reciprocal segmental gene conversion mechanism where donor cassettes recombine into the tprK expression site. Although previous studies have shown the significant role of immune selection in driving accumulation of TprK variants, the contribution of baseline gene conversion activity to variant diversity is less clear. Here, combining longitudinal tprK deep sequencing of near clonal Chicago C from immunocompetent and immunosuppressed rabbits along with the newly developed in vitro cultivation system for T. pallidum, we directly characterized TprK alleles in the presence and absence of immune selection. Our data confirm significantly greater sequence diversity over time within the V6 region during syphilis infection in immunocompetent rabbits compared to immunosuppressed rabbits, consistent with previous studies on the role of TprK in evasion of the host immune response. Compared to strains grown in immunocompetent rabbits, strains passaged in vitro displayed low level changes in allele frequencies of TprK variable region sequences similar to that of strains passaged in immunosuppressed rabbits. Notably, we found significantly increased rates of V6 allele generation relative to other variable regions in in vitro cultivated T, pallidum strains, illustrating that the diversity within these hypervariable regions occurs in the complete absence of immune selection. Together, our results demonstrate antigenic variation in T. pallidum can be studied in vitro and occurs even in the complete absence of immune pressure, allowing the T. pallidum population to continuously evade the immune system of the infected host.     

Genotypic variations in field isolates of Theileria species infecting giraffes (Giraffa camelopardalis tippelskirchi and Giraffa camelopardalis reticulata) in Kenya

Recently, mortalities among giraffes, attributed to infection with unique species of piroplasms were reported in South Africa. Although haemoparasites are known to occur in giraffes of Kenya, the prevalence, genetic diversity and pathogenicity of these parasites have not been investigated.In this study, blood samples from 13 giraffes in Kenya were investigated microscopically and genomic DNA extracted. PCR amplicons of the hyper-variable region 4 (V4) of Theileria spp. small subunit ribosomal RNA (18S rRNA) gene were hybridized to a panel of genus- and species-specific oligonucleotide probes by reverse line blot (RLB). Two newly designed oligonucleotide probes specific for previously identified Theileria spp. of giraffes found single infections in eight of the specimens and mixed infections in the remaining five samples.Partial 18S rRNA genes were successfully amplified from 9 samples and the PCR amplicons were cloned. A total of 28 plasmid clones representing the Kenyan isolates were analyzed in the present study and compared with those of closely-related organisms retrieved from GenBank. In agreement with RLB results, the nucleotide sequence alignment indicated the presence of mixed infections in the giraffes. In addition, sequence alignment with the obtained 18S rRNA gene sequences revealed extensive microheterogeneities within and between isolates, characterized by indels in the V4 regions and point mutations outside this region. Phylogeny with 18S rRNA gene sequences from the detected parasites and those of related organisms places Theileria of giraffes into two major groups, within which are numerous clades that include the isolates reported in South Africa. Collectively, these data suggest the existence of at least two distinct Theileria species among giraffes, and extensive genetic diversity within the two parasite groups.     

The Effects of Alignment Quality,Distance Calculation Method,Sequence Filtering,and Region on the Analysis of 16S rRNA Gene-Based Studies

Pyrosequencing of PCR-amplified fragments that target variable regions within the 16S rRNA gene has quickly become a powerful method for analyzing the membership and structure of microbial communities. This approach has revealed and introduced questions that were not fully appreciated by those carrying out traditional Sanger sequencing-based methods. These include the effects of alignment quality, the best method of calculating pairwise genetic distances for 16S rRNA genes, whether it is appropriate to filter variable regions, and how the choice of variable region relates to the genetic diversity observed in full-length sequences. I used a diverse collection of 13,501 high-quality full-length sequences to assess each of these questions. First, alignment quality had a significant impact on distance values and downstream analyses. Specifically, the greengenes alignment, which does a poor job of aligning variable regions, predicted higher genetic diversity, richness, and phylogenetic diversity than the SILVA and RDP-based alignments. Second, the effect of different gap treatments in determining pairwise genetic distances was strongly affected by the variation in sequence length for a region; however, the effect of different calculation methods was subtle when determining the sample''s richness or phylogenetic diversity for a region. Third, applying a sequence mask to remove variable positions had a profound impact on genetic distances by muting the observed richness and phylogenetic diversity. Finally, the genetic distances calculated for each of the variable regions did a poor job of correlating with the full-length gene. Thus, while it is tempting to apply traditional cutoff levels derived for full-length sequences to these shorter sequences, it is not advisable. Analysis of β-diversity metrics showed that each of these factors can have a significant impact on the comparison of community membership and structure. Taken together, these results urge caution in the design and interpretation of analyses using pyrosequencing data.     

Ancient diversity and geographical sub-structuring in African buffalo Theileria parva populations revealed through metagenetic analysis of antigen-encoding loci

An infection and treatment protocol involving infection with a mixture of three parasite isolates and simultaneous treatment with oxytetracycline is currently used to vaccinate cattle against Theileria parva. While vaccination results in high levels of protection in some regions, little or no protection is observed in areas where animals are challenged predominantly by parasites of buffalo origin. A previous study involving sequencing of two antigen-encoding genes from a series of parasite isolates indicated that this is associated with greater antigenic diversity in buffalo-derived T. parva. The current study set out to extend these analyses by applying high-throughput sequencing to ex vivo samples from naturally infected buffalo to determine the extent of diversity in a set of antigen-encoding genes. Samples from two populations of buffalo, one in Kenya and the other in South Africa, were examined to investigate the effect of geographical distance on the nature of sequence diversity. The results revealed a number of significant findings. First, there was a variable degree of nucleotide sequence diversity in all gene segments examined, with the percentage of polymorphic nucleotides ranging from 10% to 69%. Second, large numbers of allelic variants of each gene were found in individual animals, indicating multiple infection events. Third, despite the observed diversity in nucleotide sequences, several of the gene products had highly conserved amino acid sequences, and thus represent potential candidates for vaccine development. Fourth, although compelling evidence for population differentiation between the Kenyan and South African T. parva parasites was identified, analysis of molecular variance for each gene revealed that the majority of the underlying nucleotide sequence polymorphism was common to both areas, indicating that much of this aspect of genetic variation in the parasite population arose prior to geographic separation.     

Assessment of cell-surface exposure and vaccinogenic potentials of Treponema pallidum candidate outer membrane proteins

Syphilis, a sexually transmitted infection caused by the spirochetal bacterium Treponema pallidum, remains a global public health problem. T. pallidum is believed to be an extracellular pathogen and, as such, the identification of T. pallidum outer membrane proteins that could serve as targets for opsonic or bactericidal antibodies has remained a high research priority for vaccine development. However, the identification of T. pallidum outer membrane proteins has remained highly elusive. Recent studies and bioinformatics have implicated four treponemal proteins as potential outer membrane proteins (TP0155, TP0326, TP0483 and TP0956). Indirect immunofluorescence assays performed on treponemes encapsulated within agarose gel microdroplets failed to provide evidence that any of these four molecules were surface-exposed in T. pallidum. Second, recombinant fusion proteins corresponding to all four candidate outer membrane proteins were used separately, or in combination, to vaccinate New Zealand White rabbits. Despite achieving high titers (>1:50,000) of serum antibodies, none of the rabbits displayed chancre immunity after intradermal challenge with viable T. pallidum.     

GAPDH gene diversity in spirochetes: a paradigm for genetic promiscuity.

In this study we have determined gap sequences from nine different spirochetes. Phylogenetic analyses of these sequences in the context of all other available eubacterial and a selection of eukaryotic Gap sequences demonstrated that the eubacterial glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene diversity encompasses at least five highly distinct gene families. Within these gene families, spirochetes show an extreme degree of sequence divergence that is probably the result of several lateral gene transfer events between spirochetes and other eubacterial phyla, and early gene duplications in the eubacterial ancestor. A Gap1 sequence from the syphilis spirochete Treponema pallidum has recently been shown to be closely related to GapC sequences from Euglenozoa. Here we demonstrate that several other spirochetal species are part of this cluster, supporting the conclusion that an interkingdom gene transfer from spirochetes to Euglenozoa must have occurred. Furthermore, we provide evidence that the GAPDH genes present in the protists Parabasalia may also be of spirochetal descent.     

Molecular evolution of the tprC, D, I, K, G, and J genes in the pathogenic genus Treponema

We investigated the evolution of 6 genes from the Treponema pallidum repeat (tpr) gene family, which encode potential virulence factors and are assumed to have evolved through gene duplication and gene conversion events. The 6 loci (tprC, D, G, J, I, and K) were sequenced and analyzed in several members of the genus Treponema, including the 3 subspecies of human T. pallidum (T. pallidum subsp. pallidum, pertenue, and endemicum), Treponema paraluiscuniculi (rabbit syphilis), and the unclassified Fribourg-Blanc (simian) isolate. Phylogenetic methods, recombination analysis, and measures of nucleotide diversity were used to investigate the evolutionary history of the tpr genes. Numerous instances of gene conversion were detected by all 3 methods including both homogenizing gene conversion that involved the entire length of the sequence as well as site-specific conversions that affected smaller regions. We determined the relative age and directionality of the gene conversion events whenever possible. Our data are also relevant to a discussion of the evolution of the treponemes themselves. Higher levels of variation exist between the human subspecies than within them, supporting the classification of the human treponemes into 3 subspecies. In contrast to published theories, the divergence and diversity of T. pallidum subsp. pertenue relative to the other subspecies does not support a much older origin of yaws at the emergence of modern human, nor is the level of divergence seen in T. pallidum subsp. pallidum consistent with a very recent (< 500 years) origin of this subspecies. In general, our results demonstrate that intragenomic recombination has played a significant role in the evolution of the studied tpr genes and emphasize that efforts to infer evolutionary history of the treponemes can be complicated if past recombination events are not recognized.     

Computer-assisted analysis of envelope protein sequences of seven human immunodeficiency virus isolates: prediction of antigenic epitopes in conserved and variable regions.

Independent isolates of human immunodeficiency virus (HIV) exhibit a striking genomic diversity, most of which is located in the viral envelope gene. Since this property of the HIV group of viruses may play an important role in the pathobiology of the virus, we analyzed the predicted amino acid sequences of the envelope proteins of seven different HIV strains, three of which represent sequential isolates from a single patient. By using a computer program that predicts the secondary protein structure and superimposes values for hydrophilicity, surface probability, and flexibility, we identified several potential antigenic epitopes in the envelope proteins of the seven different viruses. Interestingly, the majority of the predicted epitopes in the exterior envelope protein (gp120) were found in regions of high sequence variability which are interspersed with highly conserved regions among the independent viral isolates. A comparison of the sequential viral isolates revealed that changes concerning the secondary structure of the protein occurred only in regions which were predicted to be antigenic, predominantly in highly variable regions. The membrane-associated protein gp41 contains no highly variable regions; about 80% of the amino acids were found to be conserved, and only one hydrophilic area was identified as likely to be accessible to antibody recognition. These findings give insight into the secondary and possible tertiary structure of variant HIV envelope proteins and should facilitate experimental approaches directed toward the identification and fine mapping of HIV envelope proteins.     

Footprint of positive selection in Treponema pallidum subsp. pallidum genome sequences suggests adaptive microevolution of the syphilis pathogen

In the rabbit model of syphilis, infection phenotypes associated with the Nichols and Chicago strains of Treponema pallidum (T. pallidum), though similar, are not identical. Between these strains, significant differences are found in expression of, and antibody responses to some candidate virulence factors, suggesting the existence of functional genetic differences between isolates. The Chicago strain genome was therefore sequenced and compared to the Nichols genome, available since 1998. Initial comparative analysis suggested the presence of 44 single nucleotide polymorphisms (SNPs), 103 small (≤3 nucleotides) indels, and 1 large (1204 bp) insertion in the Chicago genome with respect to the Nichols genome. To confirm the above findings, Sanger sequencing was performed on most loci carrying differences using DNA from Chicago and the Nichols strain used in the original T. pallidum genome project. A majority of the previously identified differences were found to be due to errors in the published Nichols genome, while the accuracy of the Chicago genome was confirmed. However, 20 SNPs were confirmed between the two genomes, and 16 (80.0%) were found in coding regions, with all being of non-synonymous nature, strongly indicating action of positive selection. Sequencing of 16 genomic loci harboring SNPs in 12 additional T. pallidum strains, (SS14, Bal 3, Bal 7, Bal 9, Sea 81-3, Sea 81-8, Sea 86-1, Sea 87-1, Mexico A, UW231B, UW236B, and UW249C), was used to identify "Chicago-" or "Nichols -specific" differences. All but one of the 16 SNPs were "Nichols-specific", with Chicago having identical sequences at these positions to almost all of the additional strains examined. These mutations could reflect differential adaptation of the Nichols strain to the rabbit host or pathoadaptive mutations acquired during human infection. Our findings indicate that SNPs among T. pallidum strains emerge under positive selection and, therefore, are likely to be functional in nature.     

Characterization of the cytoplasmic filament protein gene (cfpA) of Treponema pallidum subsp. pallidum.

Treponema pallidum and other members of the genera Treponema, Spirochaeta, and Leptonema contain multiple cytoplasmic filaments that run the length of the organism just underneath the cytoplasmic membrane. These cytoplasmic filaments have a ribbon-like profile and consist of a major cytoplasmic filament protein subunit (CfpA, formerly called TpN83) with a relative molecular weight of approximately 80,000. Degenerate DNA primers based on N-terminal and CNBr cleavage fragment amino acid sequences of T. pallidum subsp. pallidum (Nichols) CfpA were utilized to amplify a fragment of the encoding gene (cfpA). A 6.8-kb EcoRI fragment containing all but the 5' end of cfpA was identified by hybridization with the resulting PCR product and cloned into Lambda ZAP II. The 5' region was obtained by inverse PCR, and the complete gene sequence was determined. The cfpA sequence contained a 2,034-nucleotide coding region, a putative promoter with consensus sequences (5'-TTTACA-3' for -35 and 5'-TACAAT-3' for -10) similar to the sigma70 recognition sequence of Escherichia coli and other organisms, and a putative ribosome-binding site (5'-AGGAG-3'). The deduced amino acid sequence of CfpA indicated a protein of 678 residues with a calculated molecular mass of 78.5 kDa and an estimated pI of 6.15. No significant homology to known proteins or structural motifs was found among known prokaryotic or eukaryotic sequences. Expression of a LacZ-CfpA fusion protein in E. coli was detrimental to survival and growth of the host strain and resulted in the formation of short, irregular filaments suggestive of partial self-assembly of CfpA. The cytoplasmic filaments of T. pallidum and other spirochetes appear to represent a unique form of prokaryotic intracytoplasmic inclusions.