Antigenic variation in the Lyme spirochete: detailed functional assessment of recombinational switching at vlsE in the JD1 strain of Borrelia burgdorferi

Borrelia burgdorferi is a causative agent of Lyme disease and establishes long‐term infection in mammalian hosts. Persistence is promoted by the VlsE antigenic variation system, which generates combinatorial diversity of VlsE through unidirectional, segmental gene conversion from an array of silent cassettes. Here we explore the variants generated by the vls system of strain JD1, which has divergent sequence and structural elements from the type strain B31, the only B. burgdorferi strain in which recombinational switching at vlsE has been studied in detail. We first completed the sequencing of the vls region in JD1, uncovering a previously unreported 114 bp inverted repeat sequence upstream of vlsE. A five‐week infection of WT and SCID mice was used for PacBio long read sequencing along with our recently developed VAST pipeline to analyze recombinational switching at vlsE from 40,000 sequences comprising 226,000 inferred recombination events. We show that antigenic variation in B31 and JD1 is highly similar, despite the lack of 17 bp direct repeats in JD1, a somewhat different arrangement of the silent cassettes, divergent inverted repeat sequences and general divergence in the vls sequences. We also present data that strongly suggest that dimerization is required for in vivo functionality of VlsE.     

The role of VlsE antigenic variation in the Lyme disease spirochete: persistence through a mechanism that differs from other pathogens

The linear plasmid, lp28-1, is required for persistent infection by the Lyme disease spirochete, Borrelia burgdorferi. This plasmid contains the vls antigenic variation locus, which has long been thought to be important for immune evasion. However, the role of the vls locus as a virulence factor during mammalian infection has not been clearly defined. We report the successful removal of the vls locus through telomere resolvase-mediated targeted deletion, and demonstrate the absolute requirement of this lp28-1 component for persistence in the mouse host. Moreover, successful infection of C3H/HeN mice with an lp28-1 plasmid in which the left portion was deleted excludes participation of other lp28-1 non-vls genes in spirochete virulence, persistence and the process of recombinational switching at vlsE. Data are also presented that cast doubt on an immune evasion mechanism whereby VlsE directly masks other surface antigens similar to what has been observed for several other pathogens that undergo recombinational antigenic variation.     

Genetic variation at the vlsE locus of Borrelia burgdorferi within ticks and mice over the course of a single transmission cycle

The Lyme disease spirochete, Borrelia burgdorferi, causes a persistent infection in the vertebrate host even though infected animals mount an active immune response against the spirochete. One strategy used by the spirochete to evade vertebrate host immunity is to vary the structure and expression of outer membrane antigens. The vlsE locus represents the best-studied example of antigenic variation in B. burgdorferi. During vertebrate host infection, recombination between the active vlsE locus and silent, partial vlsE copies leads to gene conversion events and the generation of novel alleles at the expression site. In the present study, we followed a population of B. burgdorferi organisms moving through vertebrate host and tick stages to complete one transmission cycle. The major goal of the study was to determine if the vlsE locus was subject to different selective pressure and/or recombination frequency at different stages of the spirochete's life cycle. We report here that the vlsE genetic diversity generated within the rodent host was maintained through the larval and nymphal tick stages. Therefore, naturally infected ticks are likely to transmit spirochete populations with multiple vlsE alleles into naive vertebrate hosts. Although vlsE genetic diversity in mice was maintained through tick stages, the dominant vlsE alleles were different between tick stages as well as between individual ticks. We propose that population-level bottlenecks experienced by spirochetes, especially during the larval-to-nymphal molt, are responsible for individual infected ticks harboring different dominant vlsE alleles. Although vlsE genetic diversity is maintained through tick stages, the VlsE protein is unlikely to be of functional importance in the vector, because the protein was expressed by very few (<1%) bacteria in the vector.     

Detailed Analysis of Sequence Changes Occurring during vlsE Antigenic Variation in the Mouse Model of Borrelia burgdorferi Infection

Lyme disease Borrelia can infect humans and animals for months to years, despite the presence of an active host immune response. The vls antigenic variation system, which expresses the surface-exposed lipoprotein VlsE, plays a major role in B. burgdorferi immune evasion. Gene conversion between vls silent cassettes and the vlsE expression site occurs at high frequency during mammalian infection, resulting in sequence variation in the VlsE product. In this study, we examined vlsE sequence variation in B. burgdorferi B31 during mouse infection by analyzing 1,399 clones isolated from bladder, heart, joint, ear, and skin tissues of mice infected for 4 to 365 days. The median number of codon changes increased progressively in C3H/HeN mice from 4 to 28 days post infection, and no clones retained the parental vlsE sequence at 28 days. In contrast, the decrease in the number of clones with the parental vlsE sequence and the increase in the number of sequence changes occurred more gradually in severe combined immunodeficiency (SCID) mice. Clones containing a stop codon were isolated, indicating that continuous expression of full-length VlsE is not required for survival in vivo; also, these clones continued to undergo vlsE recombination. Analysis of clones with apparent single recombination events indicated that recombinations into vlsE are nonselective with regard to the silent cassette utilized, as well as the length and location of the recombination event. Sequence changes as small as one base pair were common. Fifteen percent of recovered vlsE variants contained “template-independent” sequence changes, which clustered in the variable regions of vlsE. We hypothesize that the increased frequency and complexity of vlsE sequence changes observed in clones recovered from immunocompetent mice (as compared with SCID mice) is due to rapid clearance of relatively invariant clones by variable region-specific anti-VlsE antibody responses.     

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.     

The repertoire of silent pilus genes in Neisseria gonorrhoeae: evidence for gene conversion

To investigate the significance of silent gene loci for pilus antigenic variation in N. gonorrhoeae, we determined the nucleotide sequence of the major silent locus, pilS1. The pilS1 locus contains six tandem pilus gene copies linked by a 39 bp repeat sequence also present in the expression loci. All silent copies lack the common N-terminal coding sequence of pilin, containing instead variant sequence information that constitutes a semivariable (SV) and a hypervariable (HV) domain. The SV and HV domains are interspersed with short, strictly conserved (C) regions flanking small cassettes of variable sequence information. It appears that such minicassettes from silent copies can be duplicated and transferred to other silent or expression genes by means of gene conversion.     

Drosophila valois encodes a divergent WD protein that is required for Vasa localization and Oskar protein accumulation

valois (vls) was identified as a posterior group gene in the initial screens for Drosophila maternal-effect lethal mutations. Despite its early genetic identification, it has not been characterized at the molecular level until now. We show that vls encodes a divergent WD domain protein and that the three available EMS-induced point mutations cause premature stop codons in the vls ORF. We have generated a null allele that has a stronger phenotype than the EMS mutants. The vlsnull mutant shows that vls+ is required for high levels of Oskar protein to accumulate during oogenesis, for normal posterior localization of Oskar in later stages of oogenesis and for posterior localization of the Vasa protein during the entire process of pole plasm assembly. There is no evidence for vls being dependent on an upstream factor of the posterior pathway, suggesting that Valois protein (Vls) instead acts as a co-factor in the process. Based on the structure of Vls, the function of similar proteins in different systems and our phenotypic analysis, it seems likely that vls may promote posterior patterning by facilitating interactions between different molecules.     

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.     

The Mutator Gene Swi8 Effects Specific Mutations in the Mating-Type Region of Schizosaccharomyces Pombe

The swi8(+) gene of Schizosaccharomyces pombe appears to be involved in the termination step of copy synthesis during mating-type (MT) switching. Mutations in swi8 confer a general mutator phenotype and, in particular, generate specific mutations in the MT region. Sequencing of the MT cassettes of the h(90) swi8-137 mutant revealed three altered sites. One is situated at the switching (smt) signal adjacent to the H1 homology box of the expression locus mat1:1. It reduces the rate of MT switching. The alteration at the smt signal arose frequently in other h(90) swi8 strains and is probably caused by gene conversion in which the sequence adjacent to the H1 box of mat2:2 is used as template. This change might be generated during the process of MT switching when hybrid DNA formation is anomalously extended into the more heterologous region flanking the H1 homology box. In addition to the gene conversion at mat1:1, two mutations were found in the H3 homology boxes of the silent cassettes mat2:2 and mat3:3.     

The Mycoplasma hominis vaa gene displays a mosaic gene structure

Mycoplasma hominis contains a variable adherence-associated ( vaa ) gene. To classify variants of the vaa genes, we examined 42 M. hominis isolates by PCR, DNA sequencing and immunoblotting. This uncovered the existence of five gene categories. Comparison of the gene types revealed a modular composition of the Vaa proteins. The proteins constituted a conserved N-terminal part followed by a varying number of interchangeable cassettes encoding approximately 110 amino acids with conserved sequence boxes flanking the cassettes. The interchangeable cassettes showed a high mutual homology and a conserved leucine zipper motif. The smallest product contained only one cassette and the largest five. Additionally, two types of stop mutations caused by substitutions resulting in the expression of truncated Vaa proteins were observed. Our results expand the known potential of the Vaa system in generating antigen variation.     

Phylogenetic and evolutionary implications of nuclear ribosomal DNA variation in dwarf dandelions (Krigia,Lactuceae, Asteraceae)

Restriction site and length variations of nrDNA were examined for 51 populations of seven species ofKrigia. The nrDNA repeat ranged in size from 8.7 to 9.6 kilobase (kb). The transcribed region, including the two ITSs, was 5.35 kb long in all examinedKrigia populations. In contrast, the size of the nontranscribed IGS varied from 3.35 to 4.25 kb. Eight different types of length-variations were identified among the 51 populations, including distinct nrDNA lengths in the tetraploid and diploid populations of bothK. biflora andK. virginica. However, a few variations were detected among populations of the same species or within a cytotype. All populations ofKrigia sect.Cymbia share a 600 bp insertion in IGS near the 18 S gene, and this feature suggests monophyly of the section. AllKrigia spp. had a conjugated type of subrepeat composed of approximately 75 basepairs (bp) and 125 bp. Base modifications in the gene coding regions were highly conserved among species. Forty-five restriction sites from 15 enzymes were mapped, 24 of which were variable among populations. Only four of the variable sites occurred in the rRNA coding region while 20 variable sites were detected in the noncoding regions. Collectively, 25 enzymes generated about 66 restriction sites in each nrDNA; this amounts to about 4.3% of the nrDNA repeat. A total of 50 restriction sites was variable, 28 of which were phylogenetically informative. Phylogenetic analyses of site mutations indicated that two sections ofKrigia, sect.Cymbia and sect.Krigia, are monophyletic. In addition, relationships among several species were congruent with other sources of data, such as cpDNA restriction site variation and morphology. Both length and restriction site variation supported an allopolyploid origin of the hexaploidK. montana. The average sequence divergence value inKrigia nrDNA was 40 times greater than that of the chloroplast DNA. The rapid evolution of nrDNA sequences was primarily due to changes of the IGS sequences.     

Efficient silencing of gene expression with modular trimeric Pol II expression cassettes comprising microRNA shuttles

Expressed polycistronic microRNA (miR) cassettes have useful properties that can be utilized for RNA interference (RNAi)-based gene silencing. To advance their application we generated modular trimeric anti-hepatitis B virus (HBV) Pol II cassettes encoding primary (pri)-miR-31-derived shuttles that target three different viral genome sites. A panel of six expression cassettes, comprising each of the possible ordering combinations of the pri-miR-31 shuttles, was initially tested. Effective silencing of individual target sequences was achieved in transfected cells and transcribed pri-miR trimers generated intended guide strands. There was, however, variation in processing and silencing by each of the shuttles. In some cases the monomers’ position within the trimers influenced processing and this correlated with target silencing. Compromised efficacy could be compensated by substituting the pri-miR-31 backbone with a pri-miR-30a scaffold. Inhibition of HBV replication was achieved in vivo, and in cell culture without disruption of endogenous miR function or induction of the interferon response. A mutant HBV target sequence, with changes in one of the guide cognates, was also silenced by the trimeric cassettes. The modular nature of the cassettes together with compatibility with expression from Pol II promoters should be advantageous for gene silencing applications requiring simultaneous targeting of different sites.     

A simple method for genome-wide screening for advantageous insertions of mobile DNAs in Escherichia coli

Laboratory evolution in Escherichia coli has revealed that fitness typically increases in experimental populations. These changes are sometimes associated with changes in insertion sequence positions, some of which may themselves cause advantageous phenotypes. We have a novel and general method for identifying genes in Escherichia coli, whose knockout by mobile DNA insertions is beneficial in experimental evolution. Insertion sites in favored clones can be identified by reference to genomic information. We have implemented the method using modified Tn10 transposons bearing kanamycin and chloramphenicol resistance cassettes. Results are consistent across replicated experiments, demonstrating that the insertions are themselves creating selective advantages, rather than hitch-hiking with favorable base substitutions. The successful clones have subsequently been confirmed to have a fitness advantage relative to the progenitor strain. In experiments in shaking culture, we find that advantageous insertions usually fall in operons required in the pathways creating flagella. The method allows a rapid genome-wide screening for advantageous insertions in arbitrary environmental conditions. It allows investigation of the extent to which transient mutations generating environment-dependent selective advantages may help to explain the persistence of mobile DNAs in primarily clonal organisms, such as E. coli.     

Alterations of the p53 gene in nasopharyngeal carcinoma.

Nasopharyngeal carcinoma (NPC) is a malignancy which is consistently associated with the Epstein-Barr virus (EBV). The structure of the EBV genome in NPC suggests that NPC is a clonal proliferation of epithelial cells which emerges after EBV infection. The disease develops with high incidence in specific populations in discrete geographic locations, implicating possible genetic or environmental cofactors. Mutations of the p53 gene are among the most frequent genetic changes found in a large variety of human tumors. Mutations in p53 have been shown to abrogate the suppressor function of wild-type p53 and thus contribute to the transformed phenotype. To determine if mutation in p53 participates in the development of the malignant clone in NPC, the structure and sequence of p53 in 42 primary, metastatic, and nude mouse-passaged NPC specimens was analyzed. A high frequency (6 of 9) of mutations was detected in the nude mouse-passaged tumors, while only 2 of 15 metastatic and 0 of the 18 primary tumors harbored mutant p53. The p53 mutations included single-point mutations and more extensive changes such as frame shifts, deletion, duplication, or complete loss of coding sequences. These data indicate that alterations of the p53 gene are unlikely to be involved in the initial genetic events leading to the clonal outgrowth in NPC. However, although it is a rare NPC which can be established in nude mice, this growth advantage appears to be conferred on tumors bearing a mutant p53.     

Central Role of the Holliday Junction Helicase RuvAB in vlsE Recombination and Infectivity of Borrelia burgdorferi

Antigenic variation plays a vital role in the pathogenesis of many infectious bacteria and protozoa including Borrelia burgdorferi, the causative agent of Lyme disease. VlsE, a 35 kDa surface-exposed lipoprotein, undergoes antigenic variation during B. burgdorferi infection of mammalian hosts, and is believed to be a critical mechanism by which the spirochetes evade immune clearance. Random, segmental recombination between the expressed vlsE gene and adjacent vls silent cassettes generates a large number of different VlsE variants within the infected host. Although the occurrence and importance of vlsE sequence variation is well established, little is known about the biological mechanism of vlsE recombination. To identify factors important in antigenic variation and vlsE recombination, we screened transposon mutants of genes known to be involved in DNA recombination and repair for their effects on infectivity and vlsE recombination. Several mutants, including those in BB0023 (ruvA), BB0022 (ruvB), BB0797 (mutS), and BB0098 (mutS-II), showed reduced infectivity in immunocompetent C3H/HeN mice. Mutants in ruvA and ruvB exhibited greatly reduced rates of vlsE recombination in C3H/HeN mice, as determined by restriction fragment polymorphism (RFLP) screening and DNA sequence analysis. In severe combined immunodeficiency (C3H/scid) mice, the ruvA mutant retained full infectivity; however, all recovered clones retained the ‘parental’ vlsE sequence, consistent with low rates of vlsE recombination. These results suggest that the reduced infectivity of ruvA and ruvB mutants is the result of ineffective vlsE recombination and underscores the important role that vlsE recombination plays in immune evasion. Based on functional studies in other organisms, the RuvAB complex of B. burgdorferi may promote branch migration of Holliday junctions during vlsE recombination. Our findings are consistent with those in the accompanying article by Dresser et al., and together these studies provide the first examples of trans-acting factors involved in vlsE recombination.     

Signatures of selection and gene conversion associated with human color vision variation

Trichromatic color vision in humans results from the combination of red, green, and blue photopigment opsins. Although color vision genes have been the targets of active molecular and psychophysical research on color vision abnormalities, little is known about patterns of normal genetic variation in these genes among global human populations. The current study presents nucleotide sequence analyses and tests of neutrality for a 5.5-kb region of the X-linked long-wave "red" opsin gene (OPN1LW) in 236 individuals from ethnically diverse human populations. Our analysis of the recombination landscape across OPN1LW reveals an unusual haplotype structure associated with amino acid replacement variation in exon 3 that is consistent with gene conversion. Compared with the absence of OPN1LW amino acid replacement fixation since divergence from chimpanzee, the human population exhibits a significant excess of high-frequency OPN1LW replacements. Our results suggest that subtle changes in L-cone opsin wavelength absorption may have been adaptive during human evolution.