The use of high-throughput DNA sequencing in the investigation of antigenic variation: application to Neisseria species

A droplet digital PCR assay to measure pilin antigenic variation frequency in Neisseria gonorrhoeae

The strict human pathogen Neisseria gonorrhoeae (gonococcus [Gc]) infects an estimated 82 million individuals globally and is a World Health Organization-designated bacterial pathogen of public health importance due to escalating antimicrobial resistance. Gc vaccines have been hindered by Gc's ability to evade immune surveillance in part by varying its major surface antigens like the type IV pilus. We developed a quick and precise method for measuring pilin antigenic variation (Av) frequency using droplet digital PCR (ddPCR) technology. Two fluorescent probes were designed to detect either the hypervariable tail region of silent pilin locus pilS3-copy 1 (S3C1) or a sequence conserved in all pilE variants (CYS2). The appropriate frequency of pilin antigenic variation is measured by the proportion of pilE amplicons carrying the recombinant S3C1 copy relative to the total pilE amplicons measured by CYS2. The ddPCR assay is specific for RecA-dependent pilin antigenic variation. The reduced frequency of pilin Av in strains lacking RecA-modulating recombination protein RecX and the DNA helicase RecQ confirms the ability of the assay to measure changes in pilin Av frequency. We used the ddPCR assay to determine that pilin Av frequency is altered by the colony densities on a solid medium. The ddPCR assay is an accurate, efficient way to measure Gc pilin Av frequency.

IMPORTANCE

Gonorrhea is a sexually transmitted infectious disease of the human genital and nasopharyngeal mucosa caused by the host-restricted bacterium Neisseria gonorrhoeae. The rise of antibiotic-resistant gonorrhea is an urgent global threat to public health. Pilus antigenic variation is a gene conversion process that allows N. gonorrhoeae to evade host immune surveillance, and a mechanistic understanding of this process is crucial to understanding N. gonorrhoeae pathogenesis. This report shows that we can adopt a digital PCR methodology to quickly and accurately measure pilin antigenic variation.

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.

PacBio Amplicon Sequencing Method To Measure Pilin Antigenic Variation Frequencies of Neisseria gonorrhoeae

Gene diversification is a common mechanism pathogens use to alter surface structures to aid in immune avoidance. Neisseria gonorrhoeae uses a gene conversion-based diversification system to alter the primary sequence of the gene encoding the major subunit of the pilus, pilE Antigenic variation occurs when one of the nonexpressed 19 silent copies donates part of its DNA sequence to pilE We have developed a method using Pacific Biosciences (PacBio) amplicon sequencing and custom software to determine pilin antigenic variation frequencies. The program analyzes 37 variable regions across the strain FA1090 1-81-S2 pilE gene and can be modified to determine sequence variation from other starting pilE sequences or other diversity generation systems. Using this method, we measured pilin antigenic variation frequencies for various derivatives of strain FA1090 and showed we can also analyze pilin antigenic variation frequencies during macrophage infection.IMPORTANCE Diversity generation systems are used by many unicellular organism to provide subpopulations of cell with different properties that are available when needed. We have developed a method using the PacBio DNA sequencing technology and a custom computer program to analyze the pilin antigenic variation system of the organism that is the sole cause of the sexually transmitted infection, gonorrhea.

Analysis of Pilin Antigenic Variation in Neisseria meningitidis by Next-Generation Sequencing

Many pathogenic microbes evade host immune surveillance by varying the surface antigens, a process termed antigenic variation. While the process of pilin antigenic variation has been extensively studied in the human pathogen Neisseria gonorrhoeae (gonococcus [Gc]), relatively few studies of pilin antigenic variation have been conducted with Neisseria meningitidis (meningococcus [Mc]). Mc is usually a commensal organism that colonizes the human nasopharynx, but when it translocates to the bloodstream or meninges, it results in the severe and often deadly meningococcal disease. The type IV pili of Mc isolates play a critical role in host surface adherence, and its major pilin component (PilE) can undergo antigenic variation. In this study, Roche 454 pyrosequencing was used to examine the pilin antigenic variation of Mc strain 8013, as well as 8013 recA, recX, recQ, rep, and recJ mutants, Gc orthologues which have been shown to play a role in pilin antigenic variation. This study confirms that the Mc recA, rep, and recJ genes are essential for pilin antigenic variation. While the Mc recQ and recX gene products contribute to normal frequencies of antigenic variation, the loss of these factors does not alter the types of pilin variants produced. Overall, this study shows that the mechanisms of pilin antigenic variation are conserved between Gc and Mc.IMPORTANCE Antigenic variation is a strategy used by many pathogens to escape host immune surveillance and establish persistent infections. This study successfully applies next-generation sequencing to study pilin antigenic variation in the human pathogen Neisseria meningitidis This assay provides an affordable and efficient solution for quantifying antigenic variation frequency in mutant strains and for defining the recombination products of the process. We determined that there is a nonuniformity of silent donor copies used during meningococcus antigenic variation, and by the analysis of selected mutants deficient for specific recombination pathways, we show for the first time that the processes are conserved between N. meningitidis and Neisseria gonorrhoeae.

Analysis of recombinational switching at the antigenic variation locus of the Lyme spirochete using a novel PacBio sequencing pipeline

The Lyme disease spirochete evades the host immune system by combinatorial variation of VlsE, a surface antigen. Antigenic variation occurs via segmental gene conversion from contiguous silent cassettes into the vlsE locus. Because of the high degree of similarity between switch variants and the size of vlsE, short-read NGS technologies have been unsuitable for sequencing vlsE populations. Here we use PacBio sequencing technology coupled with the first fully-automated software pipeline (VAST) to accurately process NGS data by minimizing error frequency, eliminating heteroduplex errors and accurately aligning switch variants. We extend earlier studies by showing use of almost all of the vlsE SNP repertoire. In different tissues of the same mouse, 99.6% of the variants were unique, suggesting that dissemination of Borrelia burgdorferi is predominantly unidirectional with little tissue-to-tissue hematogenous dissemination. We also observed a similar number of variants in SCID and wild-type mice, a heatmap of location and frequency of amino acid changes on the 3D structure and note differences observed in SCID versus wild type mice that hint at possible amino acid function. Our observed selection against diversification of residues at the dimer interface in wild-type mice strongly suggests that dimerization is required for in vivo functionality of vlsE.

Analyzing Neisseria gonorrhoeae Pilin Antigenic Variation Using 454 Sequencing Technology

Many pathogens use homologous recombination to vary surface antigens in order to avoid immune surveillance. Neisseria gonorrhoeae, the bacterium responsible for the sexually transmitted infection gonorrhea, achieves this in part by changing the sequence of the major subunit of the type IV pilus in a process termed pilin antigenic variation (Av). The N. gonorrhoeae chromosome contains one expression locus (pilE) and many promoterless, partial-coding silent copies (pilS) that act as reservoirs for variant pilin information. Pilin Av occurs by high-frequency gene conversion reactions, which transfer pilS sequences into the pilE locus. We have developed a 454 sequencing-based assay to analyze the frequency and characteristics of pilin Av that allows a more robust analysis of pilin Av than previous assays. We used this assay to analyze mutations and conditions previously shown to affect pilin Av, confirming many but not all of the previously reported phenotypes. We show that mutations or conditions that cause growth defects can result in Av phenotypes when analyzed by phase variation-based assays. Adapting the 454 sequencing to analyze pilin Av demonstrates the utility of this technology to analyze any diversity generation system that uses recombination to develop biological diversity.

IMPORTANCE

Measuring and analyzing complex recombination-based systems constitute a major barrier to understanding the mechanisms used to generate diversity. We have analyzed the contributions of many gonococcal mutations or conditions to the process of pilin antigenic variation.

Positive Selection Pressure Drives Variation on the Surface-Exposed Variable Proteins of the Pathogenic Neisseria

Pathogenic species of Neisseria utilize variable outer membrane proteins to facilitate infection and proliferation within the human host. However, the mechanisms behind the evolution of these variable alleles remain largely unknown due to analysis of previously limited datasets. In this study, we have expanded upon the previous analyses to substantially increase the number of analyzed sequences by including multiple diverse strains, from various geographic locations, to determine whether positive selective pressure is exerted on the evolution of these variable genes. Although Neisseria are naturally competent, this analysis indicates that only intrastrain horizontal gene transfer among the pathogenic Neisseria principally account for these genes exhibiting linkage equilibrium which drives the polymorphisms evidenced within these alleles. As the majority of polymorphisms occur across species, the divergence of these variable genes is dependent upon the species and is independent of geographical location, disease severity, or serogroup. Tests of neutrality were able to detect strong selection pressures acting upon both the opa and pil gene families, and were able to locate the majority of these sites within the exposed variable regions of the encoded proteins. Evidence of positive selection acting upon the hypervariable domains of Opa contradicts previous beliefs and provides evidence for selection of receptor binding. As the pathogenic Neisseria reside exclusively within the human host, the strong selection pressures acting upon both the opa and pil gene families provide support for host immune system pressure driving sequence polymorphisms within these variable genes.

Neisseria meningitidis Type IV Pili Composed of Sequence Invariable Pilins Are Masked by Multisite Glycosylation

The ability of pathogens to cause disease depends on their aptitude to escape the immune system. Type IV pili are extracellular filamentous virulence factors composed of pilin monomers and frequently expressed by bacterial pathogens. As such they are major targets for the host immune system. In the human pathogen Neisseria meningitidis, strains expressing class I pilins contain a genetic recombination system that promotes variation of the pilin sequence and is thought to aid immune escape. However, numerous hypervirulent clinical isolates express class II pilins that lack this property. This raises the question of how they evade immunity targeting type IV pili. As glycosylation is a possible source of antigenic variation it was investigated using top-down mass spectrometry to provide the highest molecular precision on the modified proteins. Unlike class I pilins that carry a single glycan, we found that class II pilins display up to 5 glycosylation sites per monomer on the pilus surface. Swapping of pilin class and genetic background shows that the pilin primary structure determines multisite glycosylation while the genetic background determines the nature of the glycans. Absence of glycosylation in class II pilins affects pilus biogenesis or enhances pilus-dependent aggregation in a strain specific fashion highlighting the extensive functional impact of multisite glycosylation. Finally, molecular modeling shows that glycans cover the surface of class II pilins and strongly decrease antibody access to the polypeptide chain. This strongly supports a model where strains expressing class II pilins evade the immune system by changing their sugar structure rather than pilin primary structure. Overall these results show that sequence invariable class II pilins are cloaked in glycans with extensive functional and immunological consequences.

During infection pathogens and their host engage in a series of measures and counter-measures to promote their own survival: pathogens express virulence factors, the immune system targets these surface structures and pathogens modify them to evade detection. Like numerous bacterial pathogens, Neisseria meningitidis express type IV pili, long filamentous adhesive structures composed of pilins. Intriguingly the amino acid sequences of pilins from most hypervirulent strains do not vary, raising the question of how they evade the immune system. This study shows that the pilus structure is completely coated with sugars thus limiting access of antibodies to the pilin polypeptide chain. We propose that multisite glycosylation and thus variation in the type of sugar mediates immune evasion in these strains.

pilS loci in Neisseria gonorrhoeae are transcriptionally active

Piliation is an important virulence determinant for Neisseria gonorrhoeae. PilE polypeptide is the major protein subunit in the pilus organelle and engages in extensive antigenic variation due to recombination between pilE and a pilS locus. pilS were so-named as they are believed to be transcriptionally silent, in contrast to the pilE locus. In this study, we demonstrate the presence of a small, pil-specific RNA species. Through using a series of pilE deletion mutants, we show by Northern blotting and quantitative reverse transcriptase PCR analysis (qRT-PCR), that these smaller RNA species are not derived from the primary pilE transcript following some processing events, but rather, arose through transcription of the pilS loci. Small transcriptome analysis, in conjunction with analysis of pilS recombinants, identified both sense and anti-sense RNAs originating from most, but not all, of the pilS gene copies. Focusing on the MS11 pilS6 locus, we identified by site-directed mutagenesis a sense promoter located immediately upstream of pilS6 copy 2, as well as an anti-sense promoter immediately downstream of pilS6 copy 1. Whole transcriptome analysis also revealed the presence of pil-specific sRNA in both gonococci and meningococci. Overall, this study reveals an added layer of complexity to the pilE/pilS recombination scheme by demonstrating pil-specific transcription within genes that were previously thought to be transcriptionally silent.

The genetics of Neisseria species

Neisseria gonorrhoeae and Neisseria meningitidis are closely related organisms that cause the sexually transmitted infection gonorrhea and serious bacterial meningitis and septicemia, respectively. Both species possess multiple mechanisms to alter the expression of surface-exposed proteins through the processes of phase and antigenic variation. This potential for wide variability in surface-exposed structures allows the organisms to always have subpopulations of divergent antigenic types to avoid immune surveillance and to contribute to functional variation. Additionally, the Neisseria are naturally competent for DNA transformation, which is their main means of genetic exchange. Although bacteriophages and plasmids are present in this genus, they are not as effective as DNA transformation for horizontal genetic exchange. There are barriers to genetic transfer, such as restriction-modification systems and CRISPR loci, that limit particular types of exchange. These host-restricted pathogens illustrate the rich complexity of genetics that can help define the similarities and differences of closely related organisms.

Sequence, distribution and chromosomal context of class I and class II pilin genes of Neisseria meningitidis identified in whole genome sequences

Background

Neisseria meningitidis expresses type four pili (Tfp) which are important for colonisation and virulence. Tfp have been considered as one of the most variable structures on the bacterial surface due to high frequency gene conversion, resulting in amino acid sequence variation of the major pilin subunit (PilE). Meningococci express either a class I or a class II pilE gene and recent work has indicated that class II pilins do not undergo antigenic variation, as class II pilE genes encode conserved pilin subunits. The purpose of this work was to use whole genome sequences to further investigate the frequency and variability of the class II pilE genes in meningococcal isolate collections.

Results

We analysed over 600 publically available whole genome sequences of N. meningitidis isolates to determine the sequence and genomic organization of pilE. We confirmed that meningococcal strains belonging to a limited number of clonal complexes (ccs, namely cc1, cc5, cc8, cc11 and cc174) harbour a class II pilE gene which is conserved in terms of sequence and chromosomal context. We also identified pilS cassettes in all isolates with class II pilE, however, our analysis indicates that these do not serve as donor sequences for pilE/pilS recombination. Furthermore, our work reveals that the class II pilE locus lacks the DNA sequence motifs that enable (G4) or enhance (Sma/Cla repeat) pilin antigenic variation. Finally, through analysis of pilin genes in commensal Neisseria species we found that meningococcal class II pilE genes are closely related to pilE from Neisseria lactamica and Neisseria polysaccharea, suggesting horizontal transfer among these species.

Conclusions

Class II pilins can be defined by their amino acid sequence and genomic context and are present in meningococcal isolates which have persisted and spread globally. The absence of G4 and Sma/Cla sequences adjacent to the class II pilE genes is consistent with the lack of pilin subunit variation in these isolates, although horizontal transfer may generate class II pilin diversity. This study supports the suggestion that high frequency antigenic variation of pilin is not universal in pathogenic Neisseria.