Population genetics of group B Streptococcus from maternal carriage in an ethnically diverse community in London

Introduction

Maternal immunization against Group B Streptococcus (GBS) has the potential to significantly reduce the burden of neonatal GBS infections. Population genetics of GBS from maternal carriage can offer key insights into vaccine target distribution.

Methods

In this study we characterized the population structure of GBS isolates from maternal carriage (n = 535) in an ethnically diverse community in London, using whole genome sequencing.

Results

The isolates clustered into nine clonal complexes (CCs) but the majority (95%) belonged to five lineages: CC1 (26%), CC19 (26%), CC23 (20%), CC17 (13%) and CC8/10 (10%). Nine serotypes were identified, the most common were serotypes III (26%), V (21%), II (19%) and Ia (19%). Other serotypes (Ib, IV, VI, VII, IX) represented less than 10% of all isolates each. Intra-lineage serotype diversity was observed in all major CCs but was highest in CC1, which revealed nine serotypes. Nearly all isolates (99%) carried at least one of the four alpha family protein genes (alpha, alp1, alp23, and rib). All isolates were susceptible to penicillin. We found 21% and 13% of isolates to be resistant to clarithromycin and clindamycin, respectively. Prevalence of macrolide-lincosamide-streptogramin B (MLS_B) resistance genes was 22% and they were most common in CC19 (37%) and CC1 (28%), and isolates with serotypes V (38%) and IV (32%). We identified some associations between maternal ethnicity and GBS population structure. Serotype Ib was significantly less common among the South Asian compared to Black women (S. Asian: 3/142, Black: 15/135, p = 0.03). There was also a significantly lower proportion of CC1 isolates among the White other (24/142) in comparison to Black (43/135) and S. Asian (44/142) women (p = 0.04). We found a significantly higher proportion of CC17 isolates among the White other compared to S. Asian women (White other: 32/142, S. Asian: 10/142, p = 0.004).

Conclusion

Our study showed high prevalence of GBS vaccine targets among isolates from pregnant women in London. However, the observed serotype diversity in CC1 and high prevalence of MLS_B resistance genes in CC19 demonstrates presence of high risk lineages, which might act as a reservoir of non-vaccine strains and antimicrobial resistance determinants.

Evidence of widespread endemic populations of highly multidrug resistant Klebsiella pneumoniae in hospital settings in Hanoi, Vietnam: a prospective cohort study

BACKGROUND

Patients with prolonged hospitalisation have a significant risk of carriage of and subsequent infection with extended spectrum β-lactamase (ESBL)-producing and carbapenemase-producing Klebsiella pneumoniae. However, the distinctive roles of the community and hospital environments in the transmission of ESBL-producing or carbapenemase-producing K pneumoniae remain elusive. We aimed to investigate the prevalence and transmission of K pneumoniae within and between the two tertiary hospitals in Hanoi, Viet Nam, using whole-genome sequencing.

METHODS

We did a prospective cohort study of 69 patients in intensive care units (ICUs) from two hospitals in Hanoi, Viet Nam. Patients were included if they were aged 18 years or older, admitted for longer than the mean length of stay in their ICU, and cultured K pneumoniae from their clinical samples. Longitudinally collected samples from patients (collected weekly) and the ICU environment (collected monthly) were cultured on selective media, and whole-genome sequences from K pneumoniae colonies analysed. We did phylogenetic analyses and correlated phenotypic antimicrobial susceptibility testing with genotypic features of K pneumoniae isolates. We constructed transmission networks of patient samples, relating ICU admission times and locations with genetic similarity of infecting K pneumoniae.

FINDINGS

Between June 1, 2017, and Jan 31, 2018, 69 patients were in the ICUs and eligible for inclusion, and a total of 357 K pneumoniae isolates were cultured and successfully sequenced. 228 (64%) of K pneumoniae isolates carried between two and four different ESBL-encoding and carbapenemase-encoding genes, with 164 (46%) isolates carrying genes encoding both, with high minimum inhibitory concentrations. We found a novel co-occurrence of blaKPC-2 and blaNDM-1 in 46·6% of samples from the globally successful ST15 lineage. Despite being physically and clinically separated, the two hospitals shared closely related strains carrying the same array of antimicrobial resistance genes.

INTERPRETATION

These results highlight the high prevalence of ESBL-positive carbapenem-resistant K pneumoniae in ICUs in Viet Nam. Through studying K pneumoniae ST15 in detail, we showed how important resistance genes are contained within these strains that are carried broadly by patients entering the two hospitals directly or through referral.

FUNDING

Medical Research Council Newton Fund, Ministry of Science and Technology, Wellcome Trust, Academy of Medical Sciences, Health Foundation, and National Institute for Health and Care Research Cambridge Biomedical Research Centre.

A2B-COVID: A Tool for Rapidly Evaluating Potential SARS-CoV-2 Transmission Events

Identifying linked cases of infection is a critical component of the public health response to viral infectious diseases. In a clinical context, there is a need to make rapid assessments of whether cases of infection have arrived independently onto a ward, or are potentially linked via direct transmission. Viral genome sequence data are of great value in making these assessments, but are often not the only form of data available. Here, we describe A2B-COVID, a method for the rapid identification of potentially linked cases of COVID-19 infection designed for clinical settings. Our method combines knowledge about infection dynamics, data describing the movements of individuals, and evolutionary analysis of genome sequences to assess whether data collected from cases of infection are consistent or inconsistent with linkage via direct transmission. A retrospective analysis of data from two wards at Cambridge University Hospitals NHS Foundation Trust during the first wave of the pandemic showed qualitatively different patterns of linkage between cases on designated COVID-19 and non-COVID-19 wards. The subsequent real-time application of our method to data from the second epidemic wave highlights its value for monitoring cases of infection in a clinical context.

Superspreaders drive the largest outbreaks of hospital onset COVID-19 infections

SARS-CoV-2 is notable both for its rapid spread, and for the heterogeneity of its patterns of transmission, with multiple published incidences of superspreading behaviour. Here, we applied a novel network reconstruction algorithm to infer patterns of viral transmission occurring between patients and health care workers (HCWs) in the largest clusters of COVID-19 infection identified during the first wave of the epidemic at Cambridge University Hospitals NHS Foundation Trust, UK. Based upon dates of individuals reporting symptoms, recorded individual locations, and viral genome sequence data, we show an uneven pattern of transmission between individuals, with patients being much more likely to be infected by other patients than by HCWs. Further, the data were consistent with a pattern of superspreading, whereby 21% of individuals caused 80% of transmission events. Our study provides a detailed retrospective analysis of nosocomial SARS-CoV-2 transmission, and sheds light on the need for intensive and pervasive infection control procedures.

Patterns of within-host genetic diversity in SARS-CoV-2

Monitoring the spread of SARS-CoV-2 and reconstructing transmission chains has become a major public health focus for many governments around the world. The modest mutation rate and rapid transmission of SARS-CoV-2 prevents the reconstruction of transmission chains from consensus genome sequences, but within-host genetic diversity could theoretically help identify close contacts. Here we describe the patterns of within-host diversity in 1181 SARS-CoV-2 samples sequenced to high depth in duplicate. 95.1% of samples show within-host mutations at detectable allele frequencies. Analyses of the mutational spectra revealed strong strand asymmetries suggestive of damage or RNA editing of the plus strand, rather than replication errors, dominating the accumulation of mutations during the SARS-CoV-2 pandemic. Within- and between-host diversity show strong purifying selection, particularly against nonsense mutations. Recurrent within-host mutations, many of which coincide with known phylogenetic homoplasies, display a spectrum and patterns of purifying selection more suggestive of mutational hotspots than recombination or convergent evolution. While allele frequencies suggest that most samples result from infection by a single lineage, we identify multiple putative examples of co-infection. Integrating these results into an epidemiological inference framework, we find that while sharing of within-host variants between samples could help the reconstruction of transmission chains, mutational hotspots and rare cases of superinfection can confound these analyses.

Genomic epidemiology of COVID-19 in care homes in the east of England

COVID-19 poses a major challenge to care homes, as SARS-CoV-2 is readily transmitted and causes disproportionately severe disease in older people. Here, 1167 residents from 337 care homes were identified from a dataset of 6600 COVID-19 cases from the East of England. Older age and being a care home resident were associated with increased mortality. SARS-CoV-2 genomes were available for 700 residents from 292 care homes. By integrating genomic and temporal data, 409 viral clusters within the 292 homes were identified, indicating two different patterns - outbreaks among care home residents and independent introductions with limited onward transmission. Approximately 70% of residents in the genomic analysis were admitted to hospital during the study, providing extensive opportunities for transmission between care homes and hospitals. Limiting viral transmission within care homes should be a key target for infection control to reduce COVID-19 mortality in this population.

Identifying virulence determinants of multidrug-resistant Klebsiella pneumoniae in Galleria mellonella

Infections caused by Klebsiella pneumoniae are a major public health threat. Extensively drug-resistant and even pan-resistant strains have been reported. Understanding K. pneumoniae pathogenesis is hampered by the fact that murine models of infection offer limited resolution for non-hypervirulent strains which cause the majority of infections. The insect Galleria mellonella larva is a widely used alternative model organism for bacterial pathogens. We have performed genome-scale fitness profiling of a multidrug-resistant K. pneumoniae ST258 strain during infection of G. mellonella, to determine if this model is suitable for large-scale virulence factor discovery in this pathogen. Our results demonstrated a dominant role for surface polysaccharides in infection, with contributions from siderophores, cell envelope proteins, purine biosynthesis genes and additional genes of unknown function. Comparison with a hypervirulent strain, ATCC 43816, revealed substantial overlap in important infection-related genes, as well as additional putative virulence factors specific to ST258, reflecting strain-dependent fitness effects. Our analysis also identified a role for the metalloregulatory protein NfeR (YqjI) in virulence. Overall, this study offers new insight into the infection fitness landscape of K. pneumoniae, and provides a framework for using the highly flexible and easily scalable G. mellonella infection model to dissect molecular virulence mechanisms of bacterial pathogens.

Rapid implementation of SARS-CoV-2 sequencing to investigate cases of health-care associated COVID-19: a prospective genomic surveillance study

BACKGROUND

The burden and influence of health-care associated severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is unknown. We aimed to examine the use of rapid SARS-CoV-2 sequencing combined with detailed epidemiological analysis to investigate health-care associated SARS-CoV-2 infections and inform infection control measures.

METHODS

In this prospective surveillance study, we set up rapid SARS-CoV-2 nanopore sequencing from PCR-positive diagnostic samples collected from our hospital (Cambridge, UK) and a random selection from hospitals in the East of England, enabling sample-to-sequence in less than 24 h. We established a weekly review and reporting system with integration of genomic and epidemiological data to investigate suspected health-care associated COVID-19 cases.

FINDINGS

Between March 13 and April 24, 2020, we collected clinical data and samples from 5613 patients with COVID-19 from across the East of England. We sequenced 1000 samples producing 747 high-quality genomes. We combined epidemiological and genomic analysis of the 299 patients from our hospital and identified 35 clusters of identical viruses involving 159 patients. 92 (58%) of 159 patients had strong epidemiological links and 32 (20%) patients had plausible epidemiological links. These results were fed back to clinical, infection control, and hospital management teams, leading to infection-control interventions and informing patient safety reporting.

INTERPRETATION

We established real-time genomic surveillance of SARS-CoV-2 in a UK hospital and showed the benefit of combined genomic and epidemiological analysis for the investigation of health-care associated COVID-19. This approach enabled us to detect cryptic transmission events and identify opportunities to target infection-control interventions to further reduce health-care associated infections. Our findings have important implications for national public health policy as they enable rapid tracking and investigation of infections in hospital and community settings.

FUNDING

COVID-19 Genomics UK funded by the Department of Health and Social Care, UK Research and Innovation, and the Wellcome Sanger Institute.

Epidemic of carbapenem-resistant Klebsiella pneumoniae in Europe is driven by nosocomial spread

Public health interventions to control the current epidemic of carbapenem-resistant Klebsiella pneumoniae rely on a comprehensive understanding of its emergence and spread over a wide range of geographical scales. We analysed the genome sequences and epidemiological data of >1,700 K. pneumoniae samples isolated from patients in 244 hospitals in 32 countries during the European Survey of Carbapenemase-Producing Enterobacteriaceae. We demonstrate that carbapenemase acquisition is the main cause of carbapenem resistance and that it occurred across diverse phylogenetic backgrounds. However, 477 of 682 (69.9%) carbapenemase-positive isolates are concentrated in four clonal lineages, sequence types 11, 15, 101, 258/512 and their derivatives. Combined analysis of the genetic and geographic distances between isolates with different β-lactam resistance determinants suggests that the propensity of K. pneumoniae to spread in hospital environments correlates with the degree of resistance and that carbapenemase-positive isolates have the highest transmissibility. Indeed, we found that over half of the hospitals that contributed carbapenemase-positive isolates probably experienced within-hospital transmission, and interhospital spread is far more frequent within, rather than between, countries. Finally, we propose a value of 21 for the number of single nucleotide polymorphisms that optimizes the discrimination of hospital clusters and detail the international spread of the successful epidemic lineage, ST258/512.

Complete Assembly of Escherichia coli Sequence Type 131 Genomes Using Long Reads Demonstrates Antibiotic Resistance Gene Variation within Diverse Plasmid and Chromosomal Contexts

The incidence of infections caused by extraintestinal Escherichia coli (ExPEC) is rising globally, which is a major public health concern. ExPEC strains that are resistant to antimicrobials have been associated with excess mortality, prolonged hospital stays, and higher health care costs. E. coli sequence type 131 (ST131) is a major ExPEC clonal group worldwide, with variable plasmid composition, and has an array of genes enabling antimicrobial resistance (AMR). ST131 isolates frequently encode the AMR genes blaCTX-M-14, blaCTX-M-15, and blaCTX-M-27, which are often rearranged, amplified, and translocated by mobile genetic elements (MGEs). Short DNA reads do not fully resolve the architecture of repetitive elements on plasmids to allow MGE structures encoding blaCTX-M genes to be fully determined. Here, we performed long-read sequencing to decipher the genome structures of six E. coli ST131 isolates from six patients. Most long-read assemblies generated entire chromosomes and plasmids as single contigs, in contrast to more fragmented assemblies created with short reads alone. The long-read assemblies highlighted diverse accessory genomes with blaCTX-M-15, blaCTX-M-14, and blaCTX-M-27 genes identified in three, one, and one isolates, respectively. One sample had no blaCTX-M gene. Two samples had chromosomal blaCTX-M-14 and blaCTX-M-15 genes, and the latter was at three distinct locations, likely transposed by the adjacent MGEs: ISEcp1, IS903B, and Tn2 This study showed that AMR genes exist in multiple different chromosomal and plasmid contexts, even between closely related isolates within a clonal group such as E. coli ST131.IMPORTANCE Drug-resistant bacteria are a major cause of illness worldwide, and a specific subtype called Escherichia coli ST131 causes a significant number of these infections. ST131 bacteria become resistant to treatments by modifying their DNA and by transferring genes among one another via large packages of genes called plasmids, like a game of pass-the-parcel. Tackling infections more effectively requires a better understanding of what plasmids are being exchanged and their exact contents. To achieve this, we applied new high-resolution DNA sequencing technology to six ST131 samples from infected patients and compared the output to that of an existing approach. A combination of methods shows that drug resistance genes on plasmids are highly mobile because they can jump into ST131's chromosomes. We found that the plasmids are very elastic and undergo extensive rearrangements even in closely related samples. This application of DNA sequencing technologies illustrates at a new level the highly dynamic nature of ST131 genomes.

Separating Bacteria by Capsule Amount Using a Discontinuous Density Gradient

Capsule is a key virulence factor in many bacterial species, mediating immune evasion and resistance to various physical stresses. While many methods are available to quantify and compare capsule production between different strains or mutants, there is no widely used method for sorting bacteria based on how much capsule they produce. We have developed a method to separate bacteria by capsule amount, using a discontinuous density gradient. This method is used to compare capsule amounts semi-quantitatively between cultures, to isolate mutants with altered capsule production, and to purify capsulated bacteria from complex samples. This method can also be coupled with transposon-insertion sequencing to identify genes involved in capsule regulation. Here, the method is demonstrated in detail, including how to optimize the gradient conditions for a new bacterial species or strain, and how to construct and run the density gradient.

The Capsule Regulatory Network of Klebsiella pneumoniae Defined by density-TraDISort

Klebsiella pneumoniae infections affect infants and the immunocompromised, and the recent emergence of hypervirulent and multidrug-resistant K. pneumoniae lineages is a critical health care concern. Hypervirulence in K. pneumoniae is mediated by several factors, including the overproduction of extracellular capsule. However, the full details of how K. pneumoniae capsule biosynthesis is achieved or regulated are not known. We have developed a robust and sensitive procedure to identify genes influencing capsule production, density-TraDISort, which combines density gradient centrifugation with transposon insertion sequencing. We have used this method to explore capsule regulation in two clinically relevant Klebsiella strains, K. pneumoniae NTUH-K2044 (capsule type K1) and K. pneumoniae ATCC 43816 (capsule type K2). We identified multiple genes required for full capsule production in K. pneumoniae, as well as putative suppressors of capsule in NTUH-K2044, and have validated the results of our screen with targeted knockout mutants. Further investigation of several of the K. pneumoniae capsule regulators identified-ArgR, MprA/KvrB, SlyA/KvrA, and the Sap ABC transporter-revealed effects on capsule amount and architecture, serum resistance, and virulence. We show that capsule production in K. pneumoniae is at the center of a complex regulatory network involving multiple global regulators and environmental cues and that the majority of capsule regulatory genes are located in the core genome. Overall, our findings expand our understanding of how capsule is regulated in this medically important pathogen and provide a technology that can be easily implemented to study capsule regulation in other bacterial species.IMPORTANCE Capsule production is essential for K. pneumoniae to cause infections, but its regulation and mechanism of synthesis are not fully understood in this organism. We have developed and applied a new method for genome-wide identification of capsule regulators. Using this method, many genes that positively or negatively affect capsule production in K. pneumoniae were identified, and we use these data to propose an integrated model for capsule regulation in this species. Several of the genes and biological processes identified have not previously been linked to capsule synthesis. We also show that the methods presented here can be applied to other species of capsulated bacteria, providing the opportunity to explore and compare capsule regulatory networks in other bacterial strains and species.

Dominant role of nucleotide substitution in the diversification of serotype 3 pneumococci over decades and during a single infection

Streptococcus pneumoniae of serotype 3 possess a mucoid capsule and cause disease associated with high mortality rates relative to other pneumococci. Phylogenetic analysis of a complete reference genome and 81 draft sequences from clonal complex 180, the predominant serotype 3 clone in much of the world, found most sampled isolates belonged to a clade affected by few diversifying recombinations. However, other isolates indicate significant genetic variation has accumulated over the clonal complex's entire history. Two closely related genomes, one from the blood and another from the cerebrospinal fluid, were obtained from a patient with meningitis. The pair differed in their behaviour in a mouse model of disease and in their susceptibility to antimicrobials, with at least some of these changes attributable to a mutation that up-regulated the patAB efflux pump. This indicates clinically important phenotypic variation can accumulate rapidly through small alterations to the genotype.

Harnessing the genome: development of a hierarchical typing scheme for meticillin-resistant Staphylococcus aureus

A major barrier to using genome sequencing in medical microbiology is the ability to interpret the data. New schemes that provide information about the importance of sequence variation in both clinical and public health settings are required. Meticillin-resistant Staphylococcus aureus (MRSA) is an important nosocomial pathogen that is being observed with increasing frequency in community settings. Better tools are needed to improve our understanding of its transmissibility and micro-epidemiology in order to develop effective interventions. Using DNA microarray technology we identified a set of 20 binary targets whose presence or absence could be determined by PCR, producing a PCR binary typing scheme (PCR-BT). This was combined with multi-locus sequence type-based, sequence nucleotide polymorphism typing to form a hierarchical typing scheme. When applied to a set of epidemiologically unrelated isolates, a high degree of concordance was observed with PFGE (98.8 %). The scheme was able to detect the presence or absence of an outbreak strain in eight out of nine outbreak investigations, demonstrating epidemiological concordance. PCR-BT was better than PFGE at distinguishing between outbreak strains, particularly where epidemic MRSA-15 was involved. The method developed here is a rapid, digital typing scheme for S. aureus for use in both micro- and macro-epidemiological investigations that has the advantage of being suitable for use in routine diagnostic laboratories. The targets are defined and therefore the types can be defined by any platform capable of detecting the sequences used, including whole genome sequencing.

Genome wide adaptations of Plasmodium falciparum in response to lumefantrine selective drug pressure

The combination therapy of the Artemisinin-derivative Artemether (ART) with Lumefantrine (LM) (Coartem®) is an important malaria treatment regimen in many endemic countries. Resistance to Artemisinin has already been reported, and it is feared that LM resistance (LMR) could also evolve quickly. Therefore molecular markers which can be used to track Coartem® efficacy are urgently needed. Often, stable resistance arises from initial, unstable phenotypes that can be identified in vitro. Here we have used the Plasmodium falciparum multidrug resistant reference strain V1S to induce LMR in vitro by culturing the parasite under continuous drug pressure for 16 months. The initial IC₅₀ (inhibitory concentration that kills 50% of the parasite population) was 24 nM. The resulting resistant strain V1S_LM, obtained after culture for an estimated 166 cycles under LM pressure, grew steadily in 378 nM of LM, corresponding to 15 times the IC₅₀ of the parental strain. However, after two weeks of culturing V1S_LM in drug-free medium, the IC₅₀ returned to that of the initial, parental strain V1S. This transient drug tolerance was associated with major changes in gene expression profiles: using the PFSANGER Affymetrix custom array, we identified 184 differentially expressed genes in V1S_LM. Among those are 18 known and putative transporters including the multidrug resistance gene 1 (pfmdr1), the multidrug resistance associated protein and the V-type H+ pumping pyrophosphatase 2 (pfvp2) as well as genes associated with fatty acid metabolism. In addition we detected a clear selective advantage provided by two genomic loci in parasites grown under LM drug pressure, suggesting that all, or some of those genes contribute to development of LM tolerance – they may prove useful as molecular markers to monitor P. falciparum LM susceptibility.

Transition of Plasmodium sporozoites into liver stage-like forms is regulated by the RNA binding protein Pumilio

Many eukaryotic developmental and cell fate decisions that are effected post-transcriptionally involve RNA binding proteins as regulators of translation of key mRNAs. In malaria parasites (Plasmodium spp.), the development of round, non-motile and replicating exo-erythrocytic liver stage forms from slender, motile and cell-cycle arrested sporozoites is believed to depend on environmental changes experienced during the transmission of the parasite from the mosquito vector to the vertebrate host. Here we identify a Plasmodium member of the RNA binding protein family PUF as a key regulator of this transformation. In the absence of Pumilio-2 (Puf2) sporozoites initiate EEF development inside mosquito salivary glands independently of the normal transmission-associated environmental cues. Puf2- sporozoites exhibit genome-wide transcriptional changes that result in loss of gliding motility, cell traversal ability and reduction in infectivity, and, moreover, trigger metamorphosis typical of early Plasmodium intra-hepatic development. These data demonstrate that Puf2 is a key player in regulating sporozoite developmental control, and imply that transformation of salivary gland-resident sporozoites into liver stage-like parasites is regulated by a post-transcriptional mechanism.

Citrobacter rodentium is an unstable pathogen showing evidence of significant genomic flux

Citrobacter rodentium is a natural mouse pathogen that causes attaching and effacing (A/E) lesions. It shares a common virulence strategy with the clinically significant human A/E pathogens enteropathogenic E. coli (EPEC) and enterohaemorrhagic E. coli (EHEC) and is widely used to model this route of pathogenesis. We previously reported the complete genome sequence of C. rodentium ICC168, where we found that the genome displayed many characteristics of a newly evolved pathogen. In this study, through PFGE, sequencing of isolates showing variation, whole genome transcriptome analysis and examination of the mobile genetic elements, we found that, consistent with our previous hypothesis, the genome of C. rodentium is unstable as a result of repeat-mediated, large-scale genome recombination and because of active transposition of mobile genetic elements such as the prophages. We sequenced an additional C. rodentium strain, EX-33, to reveal that the reference strain ICC168 is representative of the species and that most of the inactivating mutations were common to both isolates and likely to have occurred early on in the evolution of this pathogen. We draw parallels with the evolution of other bacterial pathogens and conclude that C. rodentium is a recently evolved pathogen that may have emerged alongside the development of inbred mice as a model for human disease.

Altered patterns of gene expression underlying the enhanced immunogenicity of radiation-attenuated schistosomes

BACKGROUND

Schistosome cercariae only elicit high levels of protective immunity against a challenge infection if they are optimally attenuated by exposure to ionising radiation that truncates their migration in the lungs. However, the underlying molecular mechanisms responsible for the altered phenotype of the irradiated parasite that primes for protection have yet to be identified.

METHODOLOGY/PRINCIPAL FINDINGS

We have used a custom microarray comprising probes derived from lung-stage parasites to compare patterns of gene expression in schistosomula derived from normal and irradiated cercariae. These were transformed in vitro and cultured for four, seven, and ten days to correspond in development to the priming parasites, before RNA extraction. At these late times after the radiation insult, transcript suppression was the principal feature of the irradiated larvae. Individual gene analysis revealed that only seven were significantly down-regulated in the irradiated versus normal larvae at the three time-points; notably, four of the protein products are present in the tegument or associated with its membranes, perhaps indicating a perturbed function. Grouping of transcripts using Gene Ontology (GO) and subsequent Gene Set Enrichment Analysis (GSEA) proved more informative in teasing out subtle differences. Deficiencies in signalling pathways involving G-protein-coupled receptors suggest the parasite is less able to sense its environment. Reduction of cytoskeleton transcripts could indicate compromised structure which, coupled with a paucity of neuroreceptor transcripts, may mean the parasite is also unable to respond correctly to external stimuli.

CONCLUSIONS/SIGNIFICANCE

The transcriptional differences observed are concordant with the known extended transit of attenuated parasites through skin-draining lymph nodes and the lungs: prolonged priming of the immune system by the parasite, rather than over-expression of novel antigens, could thus explain the efficacy of the irradiated vaccine.

Fasciola hepatica expresses multiple alpha- and beta-tubulin isotypes

We have identified five α-tubulin and six β-tubulin isotypes that are expressed in adult Fasciola hepatica. Amino acid sequence identities ranged between 72 and 95% for fluke α-tubulin and between 65 and 97% for β-tubulin isotypes. Nucleotide sequence identity ranged between 68–77% and 62–80%, respectively, for their coding sequences. Phylogenetic analysis indicated that two of the α-tubulins and two of the β-tubulins were distinctly divergent from the other trematode and nematode tubulin sequences described in this study, whereas the other isotypes segregated within the trematode clades. With regard to the proposed benzimidazole binding site on β-tubulin, three of the fluke isotypes had tyrosine at position 200 of β-tubulin, two had phenylalanine and one had leucine. All had phenylalanine at position 167 and glutamic acid at position 198. When isotype RT-PCR fragment sequences were compared between six individual flukes from the susceptible Cullompton isolate and from seven individual flukes from the two resistant isolates, Sligo and Oberon, these residues were conserved.

The complete genome sequence and comparative genome analysis of the high pathogenicity Yersinia enterocolitica strain 8081

The human enteropathogen, Yersinia enterocolitica, is a significant link in the range of Yersinia pathologies extending from mild gastroenteritis to bubonic plague. Comparison at the genomic level is a key step in our understanding of the genetic basis for this pathogenicity spectrum. Here we report the genome of Y. enterocolitica strain 8081 (serotype 0:8; biotype 1B) and extensive microarray data relating to the genetic diversity of the Y. enterocolitica species. Our analysis reveals that the genome of Y. enterocolitica strain 8081 is a patchwork of horizontally acquired genetic loci, including a plasticity zone of 199 kb containing an extraordinarily high density of virulence genes. Microarray analysis has provided insights into species-specific Y. enterocolitica gene functions and the intraspecies differences between the high, low, and nonpathogenic Y. enterocolitica biotypes. Through comparative genome sequence analysis we provide new information on the evolution of the Yersinia. We identify numerous loci that represent ancestral clusters of genes potentially important in enteric survival and pathogenesis, which have been lost or are in the process of being lost, in the other sequenced Yersinia lineages. Our analysis also highlights large metabolic operons in Y. enterocolitica that are absent in the related enteropathogen, Yersinia pseudotuberculosis, indicating major differences in niche and nutrients used within the mammalian gut. These include clusters directing, the production of hydrogenases, tetrathionate respiration, cobalamin synthesis, and propanediol utilisation. Along with ancestral gene clusters, the genome of Y. enterocolitica has revealed species-specific and enteropathogen-specific loci. This has provided important insights into the pathology of this bacterium and, more broadly, into the evolution of the genus. Moreover, wider investigations looking at the patterns of gene loss and gain in the Yersinia have highlighted common themes in the genome evolution of other human enteropathogens.

The goal of this study was to catalogue all the genes encoded within the Y. enterocolitica genome to help us better understand how this bacterium and related bacteria cause different diseases. There are currently genome sequences (complete gene catalogues) available for two other members of this bacterial lineage, which cause dramatically different diseases: Y. pseudotuberculosis, like Y. enterocolitica, is a gut pathogen (enteropathogen) causing gastroenteritis in humans and animals. Yersinia pestis mostly resides within blood (circulating or in fleas following blood meals) and lymph tissue. It causes bubonic plague in humans and animals, and is historically known as “The Black Death.” A three-way comparison of these genomes revealed a patchwork of genes we have defined as being species- or disease-specific and genes that are common to all three Yersinia species. This has provided us with important information on shared gene functions that define the two enteropathogenic yersinias and those that differentiate them. This will help us to connect what we know about the Y. enterocolitica lifestyle within the gut to the disease it causes and its genetic makeup. We have also provided further evidence of gene-loss by Y. pestis as it has evolved from Y. pseudotuberculosis into a more acute systemic pathogen. Similar patterns of gene loss are seen in other important pathogens such as Salmonella enterica serovar Typhi.

Genome sequence of a proteolytic (Group I) Clostridium botulinum strain Hall A and comparative analysis of the clostridial genomes

Clostridium botulinum is a heterogeneous Gram-positive species that comprises four genetically and physiologically distinct groups of bacteria that share the ability to produce botulinum neurotoxin, the most poisonous toxin known to man, and the causative agent of botulism, a severe disease of humans and animals. We report here the complete genome sequence of a representative of Group I (proteolytic) C. botulinum (strain Hall A, ATCC 3502). The genome consists of a chromosome (3,886,916 bp) and a plasmid (16,344 bp), which carry 3650 and 19 predicted genes, respectively. Consistent with the proteolytic phenotype of this strain, the genome harbors a large number of genes encoding secreted proteases and enzymes involved in uptake and metabolism of amino acids. The genome also reveals a hitherto unknown ability of C. botulinum to degrade chitin. There is a significant lack of recently acquired DNA, indicating a stable genomic content, in strong contrast to the fluid genome of Clostridium difficile, which can form longer-term relationships with its host. Overall, the genome indicates that C. botulinum is adapted to a saprophytic lifestyle both in soil and aquatic environments. This pathogen relies on its toxin to rapidly kill a wide range of prey species, and to gain access to nutrient sources, it releases a large number of extracellular enzymes to soften and destroy rotting or decayed tissues.

Complete genome of acute rheumatic fever-associated serotype M5 Streptococcus pyogenes strain manfredo

Comparisons of the 1.84-Mb genome of serotype M5 Streptococcus pyogenes strain Manfredo with previously sequenced genomes emphasized the role of prophages in diversification of S. pyogenes and the close relationship between strain Manfredo and MGAS8232, another acute rheumatic fever-associated strain.

The multidrug-resistant human pathogen Clostridium difficile has a highly mobile, mosaic genome

We determined the complete genome sequence of Clostridium difficile strain 630, a virulent and multidrug-resistant strain. Our analysis indicates that a large proportion (11%) of the genome consists of mobile genetic elements, mainly in the form of conjugative transposons. These mobile elements are putatively responsible for the acquisition by C. difficile of an extensive array of genes involved in antimicrobial resistance, virulence, host interaction and the production of surface structures. The metabolic capabilities encoded in the genome show multiple adaptations for survival and growth within the gut environment. The extreme genome variability was confirmed by whole-genome microarray analysis; it may reflect the organism's niche in the gut and should provide information on the evolution of virulence in this organism.

Microarray analysis identifies genes preferentially expressed in the lung schistosomulum of Schistosoma mansoni

The lung schistosomulum of Schistosoma mansoni is a validated target of protective immunity elicited in vaccinated mice. To identify genes expressed at this stage we constructed a microarray, representing 3088 contigs and singlets, with cDNA derived from in vitro cultured larvae and used it to screen RNA from seven life-cycle stages. Clustering of genes by expression profile across the life cycle revealed a number of membrane, membrane-associated and secreted proteins up-regulated at the lung stage, that may represent potential immune targets. Two promising secreted molecules have homology to antigens with vaccine and/or immunomodulatory potential in other helminths.

The Chlamydophila abortus genome sequence reveals an array of variable proteins that contribute to interspecies variation

The obligate intracellular bacterial pathogen Chlamydophila abortus strain S26/3 (formerly the abortion subtype of Chlamydia psittaci) is an important cause of late gestation abortions in ruminants and pigs. Furthermore, although relatively rare, zoonotic infection can result in acute illness and miscarriage in pregnant women. The complete genome sequence was determined and shows a high level of conservation in both sequence and overall gene content in comparison to other Chlamydiaceae. The 1,144,377-bp genome contains 961 predicted coding sequences, 842 of which are conserved with those of Chlamydophila caviae and Chlamydophila pneumoniae. Within this conserved Cp. abortus core genome we have identified the major regions of variation and have focused our analysis on these loci, several of which were found to encode highly variable protein families, such as TMH/Inc and Pmp families, which are strong candidates for the source of diversity in host tropism and disease causation in this group of organisms. Significantly, Cp. abortus lacks any toxin genes, and also lacks genes involved in tryptophan metabolism and nucleotide salvaging (guaB is present as a pseudogene), suggesting that the genetic basis of niche adaptation of this species is distinct from those previously proposed for other chlamydial species.

Genomic plasticity of the causative agent of melioidosis, Burkholderia pseudomallei

Burkholderia pseudomallei is a recognized biothreat agent and the causative agent of melioidosis. This Gram-negative bacterium exists as a soil saprophyte in melioidosis-endemic areas of the world and accounts for 20% of community-acquired septicaemias in northeastern Thailand where half of those affected die. Here we report the complete genome of B. pseudomallei, which is composed of two chromosomes of 4.07 megabase pairs and 3.17 megabase pairs, showing significant functional partitioning of genes between them. The large chromosome encodes many of the core functions associated with central metabolism and cell growth, whereas the small chromosome carries more accessory functions associated with adaptation and survival in different niches. Genomic comparisons with closely and more distantly related bacteria revealed a greater level of gene order conservation and a greater number of orthologous genes on the large chromosome, suggesting that the two replicons have distinct evolutionary origins. A striking feature of the genome was the presence of 16 genomic islands (GIs) that together made up 6.1% of the genome. Further analysis revealed these islands to be variably present in a collection of invasive and soil isolates but entirely absent from the clonally related organism B. mallei. We propose that variable horizontal gene acquisition by B. pseudomallei is an important feature of recent genetic evolution and that this has resulted in a genetically diverse pathogenic species.

Complete genomes of two clinical Staphylococcus aureus strains: evidence for the rapid evolution of virulence and drug resistance

Staphylococcus aureus is an important nosocomial and community-acquired pathogen. Its genetic plasticity has facilitated the evolution of many virulent and drug-resistant strains, presenting a major and constantly changing clinical challenge. We sequenced the approximately 2.8-Mbp genomes of two disease-causing S. aureus strains isolated from distinct clinical settings: a recent hospital-acquired representative of the epidemic methicillin-resistant S. aureus EMRSA-16 clone (MRSA252), a clinically important and globally prevalent lineage; and a representative of an invasive community-acquired methicillin-susceptible S. aureus clone (MSSA476). A comparative-genomics approach was used to explore the mechanisms of evolution of clinically important S. aureus genomes and to identify regions affecting virulence and drug resistance. The genome sequences of MRSA252 and MSSA476 have a well conserved core region but differ markedly in their accessory genetic elements. MRSA252 is the most genetically diverse S. aureus strain sequenced to date: approximately 6% of the genome is novel compared with other published genomes, and it contains several unique genetic elements. MSSA476 is methicillin-susceptible, but it contains a novel Staphylococcal chromosomal cassette (SCC) mec-like element (designated SCC(476)), which is integrated at the same site on the chromosome as SCCmec elements in MRSA strains but encodes a putative fusidic acid resistance protein. The crucial role that accessory elements play in the rapid evolution of S. aureus is clearly illustrated by comparing the MSSA476 genome with that of an extremely closely related MRSA community-acquired strain; the differential distribution of large mobile elements carrying virulence and drug-resistance determinants may be responsible for the clinically important phenotypic differences in these strains.

The complete genome sequence and analysis of Corynebacterium diphtheriae NCTC13129

Corynebacterium diphtheriae is a Gram-positive, non-spore forming, non-motile, pleomorphic rod belonging to the genus Corynebacterium and the actinomycete group of organisms. The organism produces a potent bacteriophage-encoded protein exotoxin, diphtheria toxin (DT), which causes the symptoms of diphtheria. This potentially fatal infectious disease is controlled in many developed countries by an effective immunisation programme. However, the disease has made a dramatic return in recent years, in particular within the Eastern European region. The largest, and still on-going, outbreak since the advent of mass immunisation started within Russia and the newly independent states of the former Soviet Union in the 1990s. We have sequenced the genome of a UK clinical isolate (biotype gravis strain NCTC13129), representative of the clone responsible for this outbreak. The genome consists of a single circular chromosome of 2 488 635 bp, with no plasmids. It provides evidence that recent acquisition of pathogenicity factors goes beyond the toxin itself, and includes iron-uptake systems, adhesins and fimbrial proteins. This is in contrast to Corynebacterium's nearest sequenced pathogenic relative, Mycobacterium tuberculosis, where there is little evidence of recent horizontal DNA acquisition. The genome itself shows an unusually extreme large-scale compositional bias, being noticeably higher in G+C near the origin than at the terminus.

Comparative analysis of the genome sequences of Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica

Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica are closely related Gram-negative beta-proteobacteria that colonize the respiratory tracts of mammals. B. pertussis is a strict human pathogen of recent evolutionary origin and is the primary etiologic agent of whooping cough. B. parapertussis can also cause whooping cough, and B. bronchiseptica causes chronic respiratory infections in a wide range of animals. We sequenced the genomes of B. bronchiseptica RB50 (5,338,400 bp; 5,007 predicted genes), B. parapertussis 12822 (4,773,551 bp; 4,404 genes) and B. pertussis Tohama I (4,086,186 bp; 3,816 genes). Our analysis indicates that B. parapertussis and B. pertussis are independent derivatives of B. bronchiseptica-like ancestors. During the evolution of these two host-restricted species there was large-scale gene loss and inactivation; host adaptation seems to be a consequence of loss, not gain, of function, and differences in virulence may be related to loss of regulatory or control functions.

Sequence of Plasmodium falciparum chromosomes 1, 3-9 and 13

Since the sequencing of the first two chromosomes of the malaria parasite, Plasmodium falciparum, there has been a concerted effort to sequence and assemble the entire genome of this organism. Here we report the sequence of chromosomes 1, 3-9 and 13 of P. falciparum clone 3D7--these chromosomes account for approximately 55% of the total genome. We describe the methods used to map, sequence and annotate these chromosomes. By comparing our assemblies with the optical map, we indicate the completeness of the resulting sequence. During annotation, we assign Gene Ontology terms to the predicted gene products, and observe clustering of some malaria-specific terms to specific chromosomes. We identify a highly conserved sequence element found in the intergenic region of internal var genes that is not associated with their telomeric counterparts.

The genome sequence of Schizosaccharomyces pombe

We have sequenced and annotated the genome of fission yeast (Schizosaccharomyces pombe), which contains the smallest number of protein-coding genes yet recorded for a eukaryote: 4,824. The centromeres are between 35 and 110 kilobases (kb) and contain related repeats including a highly conserved 1.8-kb element. Regions upstream of genes are longer than in budding yeast (Saccharomyces cerevisiae), possibly reflecting more-extended control regions. Some 43% of the genes contain introns, of which there are 4,730. Fifty genes have significant similarity with human disease genes; half of these are cancer related. We identify highly conserved genes important for eukaryotic cell organization including those required for the cytoskeleton, compartmentation, cell-cycle control, proteolysis, protein phosphorylation and RNA splicing. These genes may have originated with the appearance of eukaryotic life. Few similarly conserved genes that are important for multicellular organization were identified, suggesting that the transition from prokaryotes to eukaryotes required more new genes than did the transition from unicellular to multicellular organization.

Complete genome sequence of a multiple drug resistant Salmonella enterica serovar Typhi CT18

Salmonella enterica serovar Typhi (S. typhi) is the aetiological agent of typhoid fever, a serious invasive bacterial disease of humans with an annual global burden of approximately 16 million cases, leading to 600,000 fatalities. Many S. enterica serovars actively invade the mucosal surface of the intestine but are normally contained in healthy individuals by the local immune defence mechanisms. However, S. typhi has evolved the ability to spread to the deeper tissues of humans, including liver, spleen and bone marrow. Here we have sequenced the 4,809,037-base pair (bp) genome of a S. typhi (CT18) that is resistant to multiple drugs, revealing the presence of hundreds of insertions and deletions compared with the Escherichia coli genome, ranging in size from single genes to large islands. Notably, the genome sequence identifies over two hundred pseudogenes, several corresponding to genes that are known to contribute to virulence in Salmonella typhimurium. This genetic degradation may contribute to the human-restricted host range for S. typhi. CT18 harbours a 218,150-bp multiple-drug-resistance incH1 plasmid (pHCM1), and a 106,516-bp cryptic plasmid (pHCM2), which shows recent common ancestry with a virulence plasmid of Yersinia pestis.

Genome sequence of Yersinia pestis, the causative agent of plague

The Gram-negative bacterium Yersinia pestis is the causative agent of the systemic invasive infectious disease classically referred to as plague, and has been responsible for three human pandemics: the Justinian plague (sixth to eighth centuries), the Black Death (fourteenth to nineteenth centuries) and modern plague (nineteenth century to the present day). The recent identification of strains resistant to multiple drugs and the potential use of Y. pestis as an agent of biological warfare mean that plague still poses a threat to human health. Here we report the complete genome sequence of Y. pestis strain CO92, consisting of a 4.65-megabase (Mb) chromosome and three plasmids of 96.2 kilobases (kb), 70.3 kb and 9.6 kb. The genome is unusually rich in insertion sequences and displays anomalies in GC base-composition bias, indicating frequent intragenomic recombination. Many genes seem to have been acquired from other bacteria and viruses (including adhesins, secretion systems and insecticidal toxins). The genome contains around 150 pseudogenes, many of which are remnants of a redundant enteropathogenic lifestyle. The evidence of ongoing genome fluidity, expansion and decay suggests Y. pestis is a pathogen that has undergone large-scale genetic flux and provides a unique insight into the ways in which new and highly virulent pathogens evolve.

Massive gene decay in the leprosy bacillus

Leprosy, a chronic human neurological disease, results from infection with the obligate intracellular pathogen Mycobacterium leprae, a close relative of the tubercle bacillus. Mycobacterium leprae has the longest doubling time of all known bacteria and has thwarted every effort at culture in the laboratory. Comparing the 3.27-megabase (Mb) genome sequence of an armadillo-derived Indian isolate of the leprosy bacillus with that of Mycobacterium tuberculosis (4.41 Mb) provides clear explanations for these properties and reveals an extreme case of reductive evolution. Less than half of the genome contains functional genes but pseudogenes, with intact counterparts in M. tuberculosis, abound. Genome downsizing and the current mosaic arrangement appear to have resulted from extensive recombination events between dispersed repetitive sequences. Gene deletion and decay have eliminated many important metabolic activities including siderophore production, part of the oxidative and most of the microaerophilic and anaerobic respiratory chains, and numerous catabolic systems and their regulatory circuits.

Complete DNA sequence of a serogroup A strain of Neisseria meningitidis Z2491

Neisseria meningitidis causes bacterial meningitis and is therefore responsible for considerable morbidity and mortality in both the developed and the developing world. Meningococci are opportunistic pathogens that colonize the nasopharynges and oropharynges of asymptomatic carriers. For reasons that are still mostly unknown, they occasionally gain access to the blood, and subsequently to the cerebrospinal fluid, to cause septicaemia and meningitis. N. meningitidis strains are divided into a number of serogroups on the basis of the immunochemistry of their capsular polysaccharides; serogroup A strains are responsible for major epidemics and pandemics of meningococcal disease, and therefore most of the morbidity and mortality associated with this disease. Here we have determined the complete genome sequence of a serogroup A strain of Neisseria meningitidis, Z2491. The sequence is 2,184,406 base pairs in length, with an overall G+C content of 51.8%, and contains 2,121 predicted coding sequences. The most notable feature of the genome is the presence of many hundreds of repetitive elements, ranging from short repeats, positioned either singly or in large multiple arrays, to insertion sequences and gene duplications of one kilobase or more. Many of these repeats appear to be involved in genome fluidity and antigenic variation in this important human pathogen.

The genome sequence of the food-borne pathogen Campylobacter jejuni reveals hypervariable sequences

Campylobacter jejuni, from the delta-epsilon group of proteobacteria, is a microaerophilic, Gram-negative, flagellate, spiral bacterium-properties it shares with the related gastric pathogen Helicobacter pylori. It is the leading cause of bacterial food-borne diarrhoeal disease throughout the world. In addition, infection with C. jejuni is the most frequent antecedent to a form of neuromuscular paralysis known as Guillain-Barré syndrome. Here we report the genome sequence of C. jejuni NCTC11168. C. jejuni has a circular chromosome of 1,641,481 base pairs (30.6% G+C) which is predicted to encode 1,654 proteins and 54 stable RNA species. The genome is unusual in that there are virtually no insertion sequences or phage-associated sequences and very few repeat sequences. One of the most striking findings in the genome was the presence of hypervariable sequences. These short homopolymeric runs of nucleotides were commonly found in genes encoding the biosynthesis or modification of surface structures, or in closely linked genes of unknown function. The apparently high rate of variation of these homopolymeric tracts may be important in the survival strategy of C. jejuni.

The complete nucleotide sequence of chromosome 3 of Plasmodium falciparum

Analysis of Plasmodium falciparum chromosome 3, and comparison with chromosome 2, highlights novel features of chromosome organization and gene structure. The sub-telomeric regions of chromosome 3 show a conserved order of features, including repetitive DNA sequences, members of multigene families involved in pathogenesis and antigenic variation, a number of conserved pseudogenes, and several genes of unknown function. A putative centromere has been identified that has a core region of about 2 kilobases with an extremely high (adenine + thymidine) composition and arrays of tandem repeats. We have predicted 215 protein-coding genes and two transfer RNA genes in the 1,060,106-base-pair chromosome sequence. The predicted protein-coding genes can be divided into three main classes: 52.6% are not spliced, 45.1% have a large exon with short additional 5' or 3' exons, and 2.3% have a multiple exon structure more typical of higher eukaryotes.

Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence

Countless millions of people have died from tuberculosis, a chronic infectious disease caused by the tubercle bacillus. The complete genome sequence of the best-characterized strain of Mycobacterium tuberculosis, H37Rv, has been determined and analysed in order to improve our understanding of the biology of this slow-growing pathogen and to help the conception of new prophylactic and therapeutic interventions. The genome comprises 4,411,529 base pairs, contains around 4,000 genes, and has a very high guanine + cytosine content that is reflected in the biased amino-acid content of the proteins. M. tuberculosis differs radically from other bacteria in that a very large portion of its coding capacity is devoted to the production of enzymes involved in lipogenesis and lipolysis, and to two new families of glycine-rich proteins with a repetitive structure that may represent a source of antigenic variation.