Probabilistic Mixture Models Improve Calibration of Panel-derived Tumor Mutational Burden in the Context of both Tumor-normal and Tumor-only Sequencing.

Background: Tumor mutational burden (TMB) has been investigated as a biomarker for immune checkpoint blockade (ICB) therapy. Increasingly, TMB is being estimated with gene panel-based assays (as opposed to full exome sequencing) and different gene panels cover overlapping but distinct genomic coordinates, making comparisons across panels difficult. Previous studies have suggested that standardization and calibration to exome-derived TMB be done for each panel to ensure comparability. With TMB cutoffs being developed from panel-based assays, there is a need to understand how to properly estimate exomic TMB values from different panel-based assays. Design: Our approach to calibration of panel-derived TMB to exomic TMB proposes the use of probabilistic mixture models that allow for nonlinear relationships along with heteroscedastic error. We examined various inputs including nonsynonymous, synonymous, and hotspot counts along with genetic ancestry. Using The Cancer Genome Atlas cohort, we generated a tumor-only version of the panel-restricted data by reintroducing private germline variants. Results: We were able to model more accurately the distribution of both tumor-normal and tumor-only data using the proposed probabilistic mixture models as compared with linear regression. Applying a model trained on tumor-normal data to tumor-only input results in biased TMB predictions. Including synonymous mutations resulted in better regression metrics across both data types, but ultimately a model able to dynamically weight the various input mutation types exhibited optimal performance. Including genetic ancestry improved model performance only in the context of tumor-only data, wherein private germline variants are observed. Significance: A probabilistic mixture model better models the nonlinearity and heteroscedasticity of the data as compared with linear regression. Tumor-only panel data are needed to properly calibrate tumor-only panels to exomic TMB. Leveraging the uncertainty of point estimates from these models better informs cohort stratification in terms of TMB.

ALK Mutation Status Before and After Alectinib Treatment in Locally Advanced or Metastatic ALK-Positive NSCLC: Pooled Analysis of Two Prospective Trials.

INTRODUCTION: The effectiveness of ALK receptor tyrosine kinase (ALK) inhibitors can be limited by the development of ALK resistance mutations. This exploratory analysis assessed the efficacy of alectinib in patients with NSCLC and ALK point mutations using pooled data from two single-arm phase II studies. METHODS: Studies NP28673 and NP28761 enrolled adults with locally advanced/metastatic ALK-positive NSCLC who had progressed on crizotinib. ALK mutation analysis was conducted on cell-free DNA from 187 patients post-crizotinib/pre-alectinib, and from 49 of these patients who subsequently progressed on alectinib. RESULTS: Baseline characteristics were generally balanced across patient subgroups. At baseline, 34 distinct ALK mutations were identified in 48 of 187 patients (25.7%). Median investigator-assessed progression-free survival was longer in patients without ALK single-nucleotide variants (n = 138) versus those with (n = 48): 10.2 months (95% confidence interval [CI]: 8.1-14.3) versus 5.6 months (95% CI: 4.5-10.9), respectively. Sixteen of 32 patients (50%) with ALK resistance mutations to crizotinib achieved an investigator-assessed response to alectinib (all partial responses); most of these ALK mutations were known to be sensitive to alectinib. Analysis of plasma samples obtained post-progression on alectinib revealed that 26 of 49 (53%) samples harbored 16 distinct ALK mutations, with known alectinib-resistance mutations, I1171 T/N/S, G1202R, and V1180L, observed in 15 of 49 (31%) tumors. CONCLUSIONS: Alectinib appears clinically active against ALK rearrangements and mutations, as well as several ALK variants that can cause resistance to crizotinib. The use of cell-free DNA in plasma samples may be an alternative noninvasive method for monitoring resistance mutations during therapy.

Blood-based tumor mutational burden as a predictor of clinical benefit in non-small-cell lung cancer patients treated with atezolizumab.

Although programmed death-ligand 1-programmed death 1 (PD-L1-PD-1) inhibitors are broadly efficacious, improved outcomes have been observed in patients with high PD-L1 expression or high tumor mutational burden (TMB). PD-L1 testing is required for checkpoint inhibitor monotherapy in front-line non-small-cell lung cancer (NSCLC). However, obtaining adequate tumor tissue for molecular testing in patients with advanced disease can be challenging. Thus, an unmet medical need exists for diagnostic approaches that do not require tissue to identify patients who may benefit from immunotherapy. Here, we describe a novel, technically robust, blood-based assay to measure TMB in plasma (bTMB) that is distinct from tissue-based approaches. Using a retrospective analysis of two large randomized trials as test and validation studies, we show that bTMB reproducibly identifies patients who derive clinically significant improvements in progression-free survival from atezolizumab (an anti-PD-L1) in second-line and higher NSCLC. Collectively, our data show that high bTMB is a clinically actionable biomarker for atezolizumab in NSCLC.

Using Whole Genome Analysis to Examine Recombination across Diverse Sequence Types of Staphylococcus aureus.

Staphylococcus aureus is an important clinical pathogen worldwide and understanding this organism's phylogeny and, in particular, the role of recombination, is important both to understand the overall spread of virulent lineages and to characterize outbreaks. To further elucidate the phylogeny of S. aureus, 35 diverse strains were sequenced using whole genome sequencing. In addition, 29 publicly available whole genome sequences were included to create a single nucleotide polymorphism (SNP)-based phylogenetic tree encompassing 11 distinct lineages. All strains of a particular sequence type fell into the same clade with clear groupings of the major clonal complexes of CC8, CC5, CC30, CC45 and CC1. Using a novel analysis method, we plotted the homoplasy density and SNP density across the whole genome and found evidence of recombination throughout the entire chromosome, but when we examined individual clonal lineages we found very little recombination. However, when we analyzed three branches of multiple lineages, we saw intermediate and differing levels of recombination between them. These data demonstrate that in S. aureus, recombination occurs across major lineages that subsequently expand in a clonal manner. Estimated mutation rates for the CC8 and CC5 lineages were different from each other. While the CC8 lineage rate was similar to previous studies, the CC5 lineage was 100-fold greater. Fifty known virulence genes were screened in all genomes in silico to determine their distribution across major clades. Thirty-three genes were present variably across clades, most of which were not constrained by ancestry, indicating horizontal gene transfer or gene loss.

Genomic characterization of Salmonella Cerro ST367, an emerging Salmonella subtype in cattle in the United States.

BACKGROUND: Within the last decade, Salmonella enterica subsp. enterica serovar Cerro (S. Cerro) has become one of the most common serovars isolated from cattle and dairy farm environments in the northeastern US. The fact that this serovar is commonly isolated from subclinically infected cattle and is rarely associated with human disease, despite its frequent isolation from cattle, has led to the hypothesis that this emerging serovar may be characterized by reduced virulence. We applied comparative and population genomic approaches to (i) characterize the evolution of this recently emerged serovar and to (ii) gain a better understanding of genomic features that could explain some of the unique epidemiological features associated with this serovar. RESULTS: In addition to generating a de novo draft genome for one Salmonella Cerro strain, we also generated whole genome sequence data for 26 additional S. Cerro isolates, including 16 from cattle operations in New York (NY) state, 2 from human clinical cases from NY in 2008, and 8 from diverse animal sources (7 from Washington state and 1 from Florida). All isolates sequenced in this study represent sequence type ST367. Population genomic analysis showed that isolates from the NY cattle operations form a well-supported clade within S. Cerro ST367 (designated here "NY bovine clade"), distinct from isolates from Washington state, Florida and the human clinical cases. A molecular clock analysis indicates that the most recent common ancestor of the NY bovine clade dates back to 1998, supporting the recent emergence of this clone.Comparative genomic analyses revealed several relevant genomic features of S. Cerro ST367, that may be responsible for reduced virulence of S. Cerro, including an insertion creating a premature stop codon in sopA. In addition, patterns of gene deletion in S. Cerro ST367 further support adaptation of this clone to a unique ecological or host related niche. CONCLUSIONS: Our results indicate that the increase in prevalence of S. Cerro ST367 is caused by a highly clonal subpopulation and that S. Cerro ST367 is characterized by unique genomic deletions that may indicate adaptation to specific ecological niches and possibly reduced virulence in some hosts.

Molecular characterization of the predominant influenza A(H1N1)pdm09 virus in Mexico, December 2011-February 2012.

When the A(H1N1)pdm09 pandemic influenza virus moved into the post-pandemic period, there was a worldwide predominance of the seasonal influenza A(H3N2) and B viruses. However, A(H1N1)pdm09 became the prevailing subtype in the 2011-2012 influenza season in Mexico and most of Central America. During this season, we collected nasopharyngeal swabs of individuals presenting with influenza-like illness at our institution in Mexico City. Samples were tested for seasonal A(H3N2) and B influenza viruses, as well as A(H1N1)pdm09 by real-time reverse transcription-polymerase chain reaction. Of 205 samples tested, 46% were positive to influenza, all of them A(H1N1)pdm09. The clinical characteristics of patients showed a similar pattern to the 2009 pandemic cases. Using next generation sequencing, we obtained whole genome sequences of viruses from 4 different patients, and in 8 additional viruses we performed partial Sanger sequencing of the HA segment. Non-synonymous changes found in the Mexican isolates with respect to the prototype isolate H1N1 (A/California/04/2009) included HA S69T, K163R and N260D unique to 2012 Mexican and North American isolates and located within or adjacent to HA antigenic sites; HA S143G, S185T, A197T and S203T previously reported in viruses from the 2010-2011 season, located within or adjacent to HA antigenic sites; and HA E374K located in a relevant site for membrane fusion. All Mexican isolates had an oseltamivir-sensitive genotype. Phylogenetic analysis with all 8 influenza gene segments showed that 2012 Mexican sequences formed a robust, distinct cluster. In all cases, 2012 Mexican sequences tended to group with 2010-2011 Asian and European sequences, but not with 2009 Mexican sequences, suggesting a possible recent common ancestor between these latter regions and the 2012 Mexican viruses. It remains to be defined if these viral changes represent an important antigenic drift that would enable viral immune evasion and/or affect influenza vaccine effectiveness.

Comparative analysis of genome sequences covering the seven cronobacter species.

BACKGROUND: Species of Cronobacter are widespread in the environment and are occasional food-borne pathogens associated with serious neonatal diseases, including bacteraemia, meningitis, and necrotising enterocolitis. The genus is composed of seven species: C. sakazakii, C. malonaticus, C. turicensis, C. dublinensis, C. muytjensii, C. universalis, and C. condimenti. Clinical cases are associated with three species, C. malonaticus, C. turicensis and, in particular, with C. sakazakii multilocus sequence type 4. Thus, it is plausible that virulence determinants have evolved in certain lineages. METHODOLOGY/PRINCIPAL FINDINGS: We generated high quality sequence drafts for eleven Cronobacter genomes representing the seven Cronobacter species, including an ST4 strain of C. sakazakii. Comparative analysis of these genomes together with the two publicly available genomes revealed Cronobacter has over 6,000 genes in one or more strains and over 2,000 genes shared by all Cronobacter. Considerable variation in the presence of traits such as type six secretion systems, metal resistance (tellurite, copper and silver), and adhesins were found. C. sakazakii is unique in the Cronobacter genus in encoding genes enabling the utilization of exogenous sialic acid which may have clinical significance. The C. sakazakii ST4 strain 701 contained additional genes as compared to other C. sakazakii but none of them were known specific virulence-related genes. CONCLUSIONS/SIGNIFICANCE: Genome comparison revealed that pair-wise DNA sequence identity varies between 89 and 97% in the seven Cronobacter species, and also suggested various degrees of divergence. Sets of universal core genes and accessory genes unique to each strain were identified. These gene sequences can be used for designing genus/species specific detection assays. Genes encoding adhesins, T6SS, and metal resistance genes as well as prophages are found in only subsets of genomes and have contributed considerably to the variation of genomic content. Differences in gene content likely contribute to differences in the clinical and environmental distribution of species and sequence types.

Identification and characterization of novel Salmonella mobile elements involved in the dissemination of genes linked to virulence and transmission.

The genetic diversity represented by >2,500 different Salmonella serovars provides a yet largely uncharacterized reservoir of mobile elements that can contribute to the frequent emergence of new pathogenic strains of this important zoonotic pathogen. Currently, our understanding of Salmonella mobile elements is skewed by the fact that most studies have focused on highly virulent or common serovars. To gain a more global picture of mobile elements in Salmonella, we used prediction algorithms to screen for mobile elements in 16 sequenced Salmonella genomes representing serovars for which no prior genome scale mobile element data were available. From these results, selected mobile elements underwent further analyses in the form of validation studies, comparative analyses, and PCR-based population screens. Through this analysis we identified a novel plasmid that has two cointegrated replicons (IncI1-IncFIB); this plasmid type was found in four genomes representing different Salmonella serovars and contained a virulence gene array that had not been previously identified. A Salmonella Montevideo isolate contained an IncHI and an IncN2 plasmid, which both encoded antimicrobial resistance genes. We also identified two novel genomic islands (SGI2 and SGI3), and 42 prophages with mosaic architecture, seven of them harboring known virulence genes. Finally, we identified a novel integrative conjugative element (ICE) encoding a type IVb pilus operon in three non-typhoidal Salmonella serovars. Our analyses not only identified a considerable number of mobile elements that have not been previously reported in Salmonella, but also found evidence that these elements facilitate transfer of genes that were previously thought to be limited in their distribution among Salmonella serovars. The abundance of mobile elements encoding pathogenic properties may facilitate the emergence of strains with novel combinations of pathogenic traits.

Escherichia coli serotype O55:H7 diversity supports parallel acquisition of bacteriophage at Shiga toxin phage insertion sites during evolution of the O157:H7 lineage.

Enteropathogenic Escherichia coli (EPEC) continues to be a leading cause of mortality and morbidity in children around the world. Two EPEC genomes have been fully sequenced: those of EPEC O127:H6 strain E2348/69 (United Kingdom, 1969) and EPEC O55:H7 strain CB9615 (Germany, 2003). The O55:H7 serotype is a recent precursor to the virulent enterohemorrhagic E. coli O157:H7. To explore the diversity of O55:H7 and better understand the clonal evolution of O157:H7, we fully sequenced EPEC O55:H7 strain RM12579 (California, 1974), which was collected 1 year before the first U.S. isolate of O157:H7 was identified in California. Phage-related sequences accounted for nearly all differences between the two O55:H7 strains. Additionally, O55:H7 and O157:H7 strains were tested for the presence and insertion sites of Shiga toxin gene (stx)-containing bacteriophages. Analysis of non-phage-associated genes supported core elements of previous O157:H7 stepwise evolutionary models, whereas phage composition and insertion analyses suggested a key refinement. Specifically, the placement and presence of lambda-like bacteriophages (including those containing stx) should not be considered stable evolutionary markers or be required in placing O55:H7 and O157:H7 strains within the stepwise evolutionary models. Additionally, we suggest that a 10.9-kb region (block 172) previously believed unique to O55:H7 strains can be used to identify early O157:H7 strains. Finally, we defined two subsets of O55:H7 strains that share an as-yet-unobserved or extinct common ancestor with O157:H7 strains. Exploration of O55:H7 diversity improved our understanding of the evolution of E. coli O157:H7 and suggested a key revision to accommodate existing and future configurations of stx-containing bacteriophages into current models.

A whole-genome single nucleotide polymorphism-based approach to trace and identify outbreaks linked to a common Salmonella enterica subsp. enterica serovar Montevideo pulsed-field gel electrophoresis type.

In this study, we report a whole-genome single nucleotide polymorphism (SNP)-based evolutionary approach to study the epidemiology of a multistate outbreak of Salmonella enterica subsp. enterica serovar Montevideo. This outbreak included 272 cases that occurred in 44 states between July 2009 and April 2010. A case-control study linked the consumption of salami made with contaminated black and red pepper to the outbreak. We sequenced, on the SOLiD System, 47 isolates with XbaI PFGE pattern JIXX01.0011, a common pulsed-field gel electrophoresis (PFGE) pattern associated with isolates from the outbreak. These isolates represented 20 isolates collected from human sources during the period of the outbreak and 27 control isolates collected from human, food, animal, and environmental sources before the outbreak. Based on 253 high-confidence SNPs, we were able to reconstruct a tip-dated molecular clock phylogeny of the isolates and to assign four human isolates to the actual outbreak. We developed an SNP typing assay to rapidly discriminate between outbreak-related cases and non-outbreak-related cases and tested this assay on an extended panel of 112 isolates. These results suggest that only a very small percentage of the human isolates with the outbreak PFGE pattern and obtained during the outbreak period could be attributed to the actual pepper-related outbreak (20%), while the majority (80%) of the putative cases represented background cases. This study demonstrates that next-generation-based SNP typing provides the resolution and accuracy needed for outbreak investigations of food-borne pathogens that cannot be distinguished by currently used subtyping methods.

Genome sequencing reveals diversification of virulence factor content and possible host adaptation in distinct subpopulations of Salmonella enterica.

BACKGROUND: Divergence of bacterial populations into distinct subpopulations is often the result of ecological isolation. While some studies have suggested the existence of Salmonella enterica subsp. enterica subclades, evidence for these subdivisions has been ambiguous. Here we used a comparative genomics approach to define the population structure of Salmonella enterica subsp. enterica, and identify clade-specific genes that may be the result of ecological specialization. RESULTS: Multi-locus sequence analysis (MLSA) and single nucleotide polymorphisms (SNPs) data for 16 newly sequenced and 30 publicly available genomes showed an unambiguous subdivision of S. enterica subsp. enterica into at least two subpopulations, which we refer to as clade A and clade B. Clade B strains contain several clade-specific genes or operons, including a ß-glucuronidase operon, a S-fimbrial operon, and cell surface related genes, which strongly suggests niche specialization of this subpopulation. An additional set of 123 isolates was assigned to clades A and B by using qPCR assays targeting subpopulation-specific SNPs and genes of interest. Among 98 serovars examined, approximately 20% belonged to clade B. All clade B isolates contained two pathogenicity related genomic islands, SPI-18 and a cytolethal distending toxin islet; a combination of these two islands was previously thought to be exclusive to serovars Typhi and Paratyphi A. Presence of ß-glucuronidase in clade B isolates specifically suggests an adaptation of this clade to the vertebrate gastrointestinal environment. CONCLUSIONS: S. enterica subsp. enterica consists of at least two subpopulations that differ specifically in genes involved in host and tissue tropism, utilization of host specific carbon and nitrogen sources and are therefore likely to differ in ecology and transmission characteristics.

Transcriptional profiling of the iron starvation response in Bordetella pertussis provides new insights into siderophore utilization and virulence gene expression.

Serological studies of patients with pertussis and the identification of antigenic Bordetella pertussis proteins support the hypothesis that B. pertussis perceives an iron starvation cue and expresses multiple iron source utilization systems in its natural human host environment. Furthermore, previous studies using a murine respiratory tract infection model showed that several of these B. pertussis iron systems are required for colonization and persistence and are differentially expressed over the course of infection. The present study examined genome-wide changes in B. pertussis gene transcript abundance in response to iron starvation in vitro. In addition to known iron source utilization genes, we identified a previously uncharacterized iron-repressed cytoplasmic membrane transporter system, fbpABC, that is required for the utilization of multiple structurally distinct siderophores including alcaligin, enterobactin, ferrichrome, and desferrioxamine B. Expression of type III secretion system genes was also found to be upregulated during iron starvation in both B. pertussis strain Tohama I and Bordetella bronchiseptica strain RB50. In a survey of type III secretion system protein production by an assortment of B. pertussis laboratory-adapted and low-passage clinical isolate strains, iron limitation increased the production and secretion of the type III secretion system-specific translocation apparatus tip protein Bsp22 in all Bvg-proficient strains. These results indicate that iron starvation in the infected host is an important environmental cue influencing not only Bordetella iron transport gene expression but also the expression of other important virulence-associated genes.

Prospective genomic characterization of the German enterohemorrhagic Escherichia coli O104:H4 outbreak by rapid next generation sequencing technology.

An ongoing outbreak of exceptionally virulent Shiga toxin (Stx)-producing Escherichia coli O104:H4 centered in Germany, has caused over 830 cases of hemolytic uremic syndrome (HUS) and 46 deaths since May 2011. Serotype O104:H4, which has not been detected in animals, has rarely been associated with HUS in the past. To prospectively elucidate the unique characteristics of this strain in the early stages of this outbreak, we applied whole genome sequencing on the Life Technologies Ion Torrent PGM™ sequencer and Optical Mapping to characterize one outbreak isolate (LB226692) and a historic O104:H4 HUS isolate from 2001 (01-09591). Reference guided draft assemblies of both strains were completed with the newly introduced PGM™ within 62 hours. The HUS-associated strains both carried genes typically found in two types of pathogenic E. coli, enteroaggregative E. coli (EAEC) and enterohemorrhagic E. coli (EHEC). Phylogenetic analyses of 1,144 core E. coli genes indicate that the HUS-causing O104:H4 strains and the previously published sequence of the EAEC strain 55989 show a close relationship but are only distantly related to common EHEC serotypes. Though closely related, the outbreak strain differs from the 2001 strain in plasmid content and fimbrial genes. We propose a model in which EAEC 55989 and EHEC O104:H4 strains evolved from a common EHEC O104:H4 progenitor, and suggest that by stepwise gain and loss of chromosomal and plasmid-encoded virulence factors, a highly pathogenic hybrid of EAEC and EHEC emerged as the current outbreak clone. In conclusion, rapid next-generation technologies facilitated prospective whole genome characterization in the early stages of an outbreak.

Comparative genomics of the bacterial genus Listeria: Genome evolution is characterized by limited gene acquisition and limited gene loss.

BACKGROUND: The bacterial genus Listeria contains pathogenic and non-pathogenic species, including the pathogens L. monocytogenes and L. ivanovii, both of which carry homologous virulence gene clusters such as the prfA cluster and clusters of internalin genes. Initial evidence for multiple deletions of the prfA cluster during the evolution of Listeria indicates that this genus provides an interesting model for studying the evolution of virulence and also presents practical challenges with regard to definition of pathogenic strains. RESULTS: To better understand genome evolution and evolution of virulence characteristics in Listeria, we used a next generation sequencing approach to generate draft genomes for seven strains representing Listeria species or clades for which genome sequences were not available. Comparative analyses of these draft genomes and six publicly available genomes, which together represent the main Listeria species, showed evidence for (i) a pangenome with 2,032 core and 2,918 accessory genes identified to date, (ii) a critical role of gene loss events in transition of Listeria species from facultative pathogen to saprotroph, even though a consistent pattern of gene loss seemed to be absent, and a number of isolates representing non-pathogenic species still carried some virulence associated genes, and (iii) divergence of modern pathogenic and non-pathogenic Listeria species and strains, most likely circa 47 million years ago, from a pathogenic common ancestor that contained key virulence genes. CONCLUSIONS: Genome evolution in Listeria involved limited gene loss and acquisition as supported by (i) a relatively high coverage of the predicted pan-genome by the observed pan-genome, (ii) conserved genome size (between 2.8 and 3.2 Mb), and (iii) a highly syntenic genome. Limited gene loss in Listeria did include loss of virulence associated genes, likely associated with multiple transitions to a saprotrophic lifestyle. The genus Listeria thus provides an example of a group of bacteria that appears to evolve through a loss of virulence rather than acquisition of virulence characteristics. While Listeria includes a number of species-like clades, many of these putative species include clades or strains with atypical virulence associated characteristics. This information will allow for the development of genetic and genomic criteria for pathogenic strains, including development of assays that specifically detect pathogenic Listeria strains.

Accurate, rapid and high-throughput detection of strain-specific polymorphisms in Bacillus anthracis and Yersinia pestis by next-generation sequencing.

BACKGROUND: In the event of biocrimes or infectious disease outbreaks, high-resolution genetic characterization for identifying the agent and attributing it to a specific source can be crucial for an effective response. Until recently, in-depth genetic characterization required expensive and time-consuming Sanger sequencing of a few strains, followed by genotyping of a small number of marker loci in a panel of isolates at or by gel-based approaches such as pulsed field gel electrophoresis, which by necessity ignores most of the genome. Next-generation, massively parallel sequencing (MPS) technology (specifically the Applied Biosystems sequencing by oligonucleotide ligation and detection (SOLiD™) system) is a powerful investigative tool for rapid, cost-effective and parallel microbial whole-genome characterization. RESULTS: To demonstrate the utility of MPS for whole-genome typing of monomorphic pathogens, four Bacillus anthracis and four Yersinia pestis strains were sequenced in parallel. Reads were aligned to complete reference genomes, and genomic variations were identified. Resequencing of the B. anthracis Ames ancestor strain detected no false-positive single-nucleotide polymorphisms (SNPs), and mapping of reads to the Sterne strain correctly identified 98% of the 133 SNPs that are not clustered or associated with repeats. Three geographically distinct B. anthracis strains from the A branch lineage were found to have between 352 and 471 SNPs each, relative to the Ames genome, and one strain harbored a genomic amplification. Sequencing of four Y. pestis strains from the Orientalis lineage identified between 20 and 54 SNPs per strain relative to the CO92 genome, with the single Bolivian isolate having approximately twice as many SNPs as the three more closely related North American strains. Coverage plotting also revealed a common deletion in two strains and an amplification in the Bolivian strain that appear to be due to insertion element-mediated recombination events. Most private SNPs (that is, a, variant found in only one strain in this set) selected for validation by Sanger sequencing were confirmed, although rare false-positive SNPs were associated with variable nucleotide tandem repeats. CONCLUSIONS: The high-throughput, multiplexing capability, and accuracy of this system make it suitable for rapid whole-genome typing of microbial pathogens during a forensic or epidemiological investigation. By interrogating nearly every base of the genome, rare polymorphisms can be reliably discovered, thus facilitating high-resolution strain tracking and strengthening forensic attribution.

Phase variation and microevolution at homopolymeric tracts in Bordetella pertussis.

BACKGROUND: Bordetella pertussis, the causative agent of whooping cough, is a highly clonal pathogen of the respiratory tract. Its lack of genetic diversity, relative to many bacterial pathogens, could limit its ability to adapt to a hostile and changing host environment. This limitation might be overcome by phase variation, as observed for other mucosal pathogens. One of the most common mechanisms of phase variation is reversible expansion or contraction of homopolymeric tracts (HPTs). RESULTS: The genomes of B. pertussis and the two closely related species, B. bronchiseptica and B. parapertussis, were screened for homopolymeric tracts longer than expected on the basis of chance, given their nucleotide compositions. Sixty-nine such HPTs were found in total among the three genomes, 74% of which were polymorphic among the three species. Nine HPTs were genotyped in a collection of 90 geographically and temporally diverse B. pertussis strains using the polymerase chain reaction/ligase detection reaction (PCR/LDR) assay. Six HPTs were polymorphic in this collection of B. pertussis strains. Of note, one of these polymorphic HPTs was found in the fimX promoter, where a single base insertion variant was present in seven strains, all of which were isolated prior to introduction of the pertussis vaccine. Transcript abundance of fimX was found to be 3.8-fold lower in strains carrying the longer allele. HPTs in three other genes, tcfA, bapC, and BP3651, varied widely in composition across the strain collection and displayed allelic polymorphism within single cultures. CONCLUSION: Allelic polymorphism at homopolymeric tracts is common within the B. pertussis genome. Phase variability may be an important mechanism in B. pertussis for evasion of the immune system and adaptation to different niches in the human host. High sensitivity and specificity make the PCR/LDR assay a powerful tool for investigating allelic variation at HPTs. Using this method, allelic diversity and phase variation were demonstrated at several B. pertussis loci.

Growth phase- and nutrient limitation-associated transcript abundance regulation in Bordetella pertussis.

To survive in a host environment, microbial pathogens must sense local conditions, including nutrient availability, and adjust their growth state and virulence functions accordingly. No comprehensive investigation of growth phase-related gene regulation in Bordetella pertussis has been reported previously. We characterized changes in genome-wide transcript abundance of B. pertussis as a function of growth phase and availability of glutamate, a key nutrient for this organism. Using a Bordetella DNA microarray, we discovered significant changes in transcript abundance for 861 array elements during the transition from log phase to stationary phase, including declining transcript levels of many virulence factor genes. The responses to glutamate depletion exhibited similarities to the responses induced by exit from log phase, including decreased virulence factor transcript levels. However, only 23% of array elements that showed at least a fourfold growth phase-associated difference in transcript abundance also exhibited glutamate depletion-associated changes, suggesting that nutrient limitation may be one of several interacting factors affecting gene regulation during stationary phase. Transcript abundance patterns of a Bvg+ phase-locked mutant revealed that the BvgAS two-component regulatory system is a key determinant of growth phase- and nutrient limitation-related transcriptional control. Several adhesin genes exhibited lower transcript abundance during stationary phase and under glutamate restriction conditions. The predicted bacterial phenotype was confirmed: adherence to bronchoepithelial cells decreased 3.3- and 4.4-fold at stationary phase and with glutamate deprivation, respectively. Growth phase and nutrient availability may serve as cues by which B. pertussis regulates virulence according to the stage of infection or the location within the human airway.

Significant gene order and expression differences in Bordetella pertussis despite limited gene content variation.

Bordetella pertussis, an obligate human pathogen and the agent of whooping cough, is a clonal species, despite the dynamic selection pressures imposed by host immunity and vaccine usage. Because the generation of variation is critical for species evolution, we employed a variety of approaches to examine features of B. pertussis genetic variation. We found a high level of conservation of gene content among 137 B. pertussis strains with different geographical, temporal, and epidemiological associations, using comparative genomic hybridization. The limited number of regions of difference were frequently located adjacent to copies of the insertion element IS481, which is present in high numbers in the B. pertussis chromosome. This repeated sequence appears to provide targets for homologous recombination, resulting in deletion of intervening sequences. Using subtractive hybridization, we searched for previously undetected genes in diverse clinical isolates but did not detect any new genes, indicating that gene acquisition is rare in B. pertussis. In contrast, we found evidence of altered gene order in the several strains that were examined and again found an association of IS481 with sites of rearrangement. Finally, we compared whole-genome expression profiles of different strains and found significant changes in transcript abundance, even in the same strain after as few as 12 laboratory passages. This combination of approaches provides a detailed picture of a pathogenic species with little gene loss or gain but with the capacity to generate variation by rearranging its chromosome and altering gene expression. These findings have broad implications for host adaptation by microbial pathogens.

A novel rhoptry protein in Toxoplasma gondii bradyzoites and merozoites.

The secretory organelles of Toxoplasma gondii orchestrate invasion of the host cell and establish the parasitophorous vacuole. Although much has been learned about the roles played by these organelles in invasion by the tachyzoite stage, little is known about the contents or functions of these organelles during bradyzoite development or pathogenesis. We identified a novel protein that localizes to the rhoptries of the bradyzoite stage, but is absent from the tachyzoite stage. This protein, BRP1, first appears in the nascent rhoptries during the first division of bradyzoite stage development. We observed secretion of BRP1 and other rhoptry proteins into the parasitophorous vacuole during bradyzoite development in vitro, but there was no evidence that this occurs in vivo. Brp1 knockout parasites did not appear to have any developmental or growth defects in vitro, and were able to establish infections in mice both as tachyzoites (via intraperitoneal injection of in vitro-derived tachyzoites) or bradyzoites (via oral gavage using cysts harvested from mouse brain). Mice infected using brain cysts from the brp1 knockout or the control strain developed similar numbers and sizes of brain cysts. Thus BRP1 does not appear to play an essential role in development of the bradyzoite stage, development of brain cysts, or oral infection of new hosts, at least in the mouse model used here. Since we also observed that BRP1 is expressed in the merozoite stages in the gut of infected cats, the coccidian phase of the life cycle may be where BRP1 plays its most important role.

Bordetella pertussis, the causative agent of whooping cough, evolved from a distinct, human-associated lineage of B. bronchiseptica.

Bordetella pertussis, B. bronchiseptica, B. parapertussis(hu), and B. parapertussis(ov) are closely related respiratory pathogens that infect mammalian species. B. pertussis and B. parapertussis(hu) are exclusively human pathogens and cause whooping cough, or pertussis, a disease that has resurged despite vaccination. Although it most often infects animals, infrequently B. bronchiseptica is isolated from humans, and these infections are thought to be zoonotic. B. pertussis and B. parapertussis(hu) are assumed to have evolved from a B. bronchiseptica-like ancestor independently. To determine the phylogenetic relationships among these species, housekeeping and virulence genes were sequenced, comparative genomic hybridizations were performed using DNA microarrays, and the distribution of insertion sequence elements was determined, using a collection of 132 strains. This multifaceted approach distinguished four complexes, representing B. pertussis, B. parapertussis(hu), and two distinct B. bronchiseptica subpopulations, designated complexes I and IV. Of the two B. bronchiseptica complexes, complex IV was more closely related to B. pertussis. Of interest, while only 32% of the complex I strains were isolated from humans, 80% of the complex IV strains were human isolates. Comparative genomic hybridization analysis identified the absence of the pertussis toxin locus and dermonecrotic toxin gene, as well as a polymorphic lipopolysaccharide biosynthesis locus, as associated with adaptation of complex IV strains to the human host. Lipopolysaccharide structural diversity among these strains was confirmed by gel electrophoresis. Thus, complex IV strains may comprise a human-associated lineage of B. bronchiseptica from which B. pertussis evolved. These findings will facilitate the study of pathogen host-adaptation. Our results shed light on the origins of the disease pertussis and suggest that the association of B. pertussis with humans may be more ancient than previously assumed.