Antibiotic resistance determination using Enterococcus faecium whole-genome sequences: a diagnostic accuracy study using genotypic and phenotypic data.

BACKGROUND: DNA sequencing could become an alternative to in vitro antibiotic susceptibility testing (AST) methods for determining antibiotic resistance by detecting genetic determinants associated with decreased antibiotic susceptibility. Here, we aimed to assess and improve the accuracy of antibiotic resistance determination from Enterococcus faecium genomes for diagnosis and surveillance purposes. METHODS: In this retrospective diagnostic accuracy study, we first conducted a literature search in PubMed on Jan 14, 2021, to compile a catalogue of genes and mutations predictive of antibiotic resistance in E faecium. We then evaluated the diagnostic accuracy of this database to determine susceptibility to 12 different, clinically relevant antibiotics using a diverse population of 4382 E faecium isolates with available whole-genome sequences and in vitro culture-based AST phenotypes. Isolates were obtained from various sources in 11 countries worldwide between 2000 and 2018. We included isolates tested with broth microdilution, Vitek 2, and disc diffusion, and antibiotics with at least 50 susceptible and 50 resistant isolates. Phenotypic resistance was derived from raw minimum inhibitory concentrations and measured inhibition diameters, and harmonised primarily using the breakpoints set by the European Committee on Antimicrobial Susceptibility Testing. A bioinformatics pipeline was developed to process raw sequencing reads, identify antibiotic resistance genetic determinants, and report genotypic resistance. We used our curated database, as well as ResFinder, AMRFinderPlus, and LRE-Finder, to assess the accuracy of genotypic predictions against phenotypic resistance. FINDINGS: We curated a catalogue of 228 genetic markers involved in resistance to 12 antibiotics in E faecium. Very accurate genotypic predictions were obtained for ampicillin (sensitivity 99·7% [95% CI 99·5-99·9] and specificity 97·9% [95·8-99·0]), ciprofloxacin (98·0% [96·4-98·9] and 98·8% [95·9-99·7]), vancomycin (98·8% [98·3-99·2] and 98·8% [98·0-99·3]), and linezolid resistance (after re-testing false negatives: 100·0% [90·8-100·0] and 98·3% [97·8-98·7]). High sensitivity was obtained for tetracycline (99·5% [99·1-99·7]), teicoplanin (98·9% [98·4-99·3]), and high-level resistance to aminoglycosides (97·7% [96·6-98·4] for streptomycin and 96·8% [95·8-97·5] for gentamicin), although at lower specificity (60-90%). Sensitivity was expectedly low for daptomycin (73·6% [65·1-80·6]) and tigecycline (38·3% [27·1-51·0]), for which the genetic basis of resistance is not fully characterised. Compared with other antibiotic resistance databases and bioinformatic tools, our curated database was similarly accurate at detecting resistance to ciprofloxacin and linezolid and high-level resistance to streptomycin and gentamicin, but had better sensitivity for detecting resistance to ampicillin, tigecycline, daptomycin, and quinupristin-dalfopristin, and better specificity for ampicillin, vancomycin, teicoplanin, and tetracycline resistance. In a validation dataset of 382 isolates, similar or improved diagnostic accuracies were also achieved. INTERPRETATION: To our knowledge, this work represents the largest published evaluation to date of the accuracy of antibiotic susceptibility predictions from E faecium genomes. The results and resources will facilitate the adoption of whole-genome sequencing as a tool for the diagnosis and surveillance of antimicrobial resistance in E faecium. A complete characterisation of the genetic basis of resistance to last-line antibiotics, and the mechanisms mediating antibiotic resistance silencing, are needed to close the remaining sensitivity and specificity gaps in genotypic predictions. FUNDING: Wellcome Trust, UK Department of Health, British Society for Antimicrobial Chemotherapy, Academy of Medical Sciences and the Health Foundation, Medical Research Council Newton Fund, Vietnamese Ministry of Science and Technology, and European Society of Clinical Microbiology and Infectious Disease.

Methodological Approach to Identify and Expand the Volume of Antimicrobial Resistance (AMR) Data in the Human Health Sector in Low- and Middle-Income Countries in Asia: Implications for Local and Regional AMR Surveillance Systems Strengthening.

Antimicrobial resistance (AMR) is a multifaceted global health problem disproportionately affecting low- and middle-income countries (LMICs). The Capturing data on Antimicrobial resistance Patterns and Trends in Use in Regions of Asia (CAPTURA) project was tasked to expand the volume of AMR and antimicrobial use data in Asia. The CAPTURA project used 2 data-collection streams: facility data and project metadata. Project metadata constituted information collected to map out data sources and assess data quality, while facility data referred to the retrospective data collected from healthcare facilities. A down-selection process, labelled "the funnel approach" by the project, was adopted to use the project metadata in prioritizing and selecting laboratories for retrospective AMR data collection. Moreover, the metadata served as a guide for understanding the AMR data once they were collected. The findings from CAPTURA's metadata add to the current discourse on the limitation of AMR data in LMICs. There is generally a low volume of AMR data generated as there is a lack of microbiology laboratories with sufficient antimicrobial susceptibility testing capacity. Many laboratories in Asia are still capturing data on paper, resulting in scattered or unused data not readily accessible or shareable for analyses. There is also a lack of clinical and epidemiological data captured, impeding interpretation and in-depth understanding of the AMR data. CAPTURA's experience in Asia suggests that there is a wide spectrum of capacity and capability of microbiology laboratories within a country and region. As local AMR surveillance is a crucial instrument to inform context-specific measures to combat AMR, it is important to understand and assess current capacity-building needs while implementing activities to enhance surveillance systems.

Capturing Data on Antimicrobial Resistance Patterns and Trends in Use in Regions of Asia (CAPTURA).

BACKGROUND: In 2015, the UK government established the Fleming Fund with the aim to address critical gaps in surveillance of antimicrobial resistance (AMR) in low- and middle-income countries in Asia and Africa. Among a large portfolio of grants, the Capturing Data on Antimicrobial Resistance Patterns and Trends in Use in Regions of Asia (CAPTURA) project was awarded with the specific objective of expanding the volume of historical data on AMR, consumption (AMC), and use (AMU) in the human healthcare sector across 12 countries in South and Southeast Asia. METHODS: Starting in early 2019, the CAPTURA consortium began working with local governments and >100 relevant data-holding facilities across the region to identify, assess for quality, prioritize, and subsequently retrieve data on AMR, AMC, and AMU. Relevant and shared data were collated and analyzed to provide local overviews for national stakeholders as well as regional context, wherever possible. RESULTS: From the vast information resource generated on current surveillance capacity and data availability, the project has highlighted gaps and areas for quality improvement and supported comprehensive capacity-building activities to optimize local data-collection and -management practices. CONCLUSIONS: The project has paved the way for expansion of surveillance networks to include both the academic and private sector in several countries and has actively engaged in discussions to promote data sharing at the local, national, and regional levels. This paper describes the overarching approach to, and emerging lessons from, the CAPTURA project, and how it contributes to other ongoing efforts to strengthen national AMR surveillance in the region and globally.

Unlocking the Potential of Genomic Data to Inform Typhoid Fever Control Policy: Supportive Resources for Genomic Data Generation, Analysis, and Visualization.

The global response to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic demonstrated the value of timely and open sharing of genomic data with standardized metadata to facilitate monitoring of the emergence and spread of new variants. Here, we make the case for the value of Salmonella Typhi (S. Typhi) genomic data and demonstrate the utility of freely available platforms and services that support the generation, analysis, and visualization of S. Typhi genomic data on the African continent and more broadly by introducing the Africa Centres for Disease Control and Prevention's Pathogen Genomics Initiative, SEQAFRICA, Typhi Pathogenwatch, TyphiNET, and the Global Typhoid Genomics Consortium.

Comparing genomic variant identification protocols for Candida auris.

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies.

Europe-wide expansion and eradication of multidrug-resistant Neisseria gonorrhoeae lineages: a genomic surveillance study.

BACKGROUND: Genomic surveillance using quality-assured whole-genome sequencing (WGS) together with epidemiological and antimicrobial resistance (AMR) data is essential to characterise the circulating Neisseria gonorrhoeae lineages and their association to patient groups (defined by demographic and epidemiological factors). In 2013, the European gonococcal population was characterised genomically for the first time. We describe the European gonococcal population in 2018 and identify emerging or vanishing lineages associated with AMR and epidemiological characteristics of patients, to elucidate recent changes in AMR and gonorrhoea epidemiology in Europe. METHODS: We did WGS on 2375 gonococcal isolates from 2018 (mainly Sept 1-Nov 30) in 26 EU and EEA countries. Molecular typing and AMR determinants were extracted from quality-checked genomic data. Association analyses identified links between genomic lineages, AMR, and epidemiological data. FINDINGS: Azithromycin-resistant N gonorrhoeae (8·0% [191/2375] in 2018) is rising in Europe due to the introduction or emergence and subsequent expansion of a novel N gonorrhoeae multi-antigen sequence typing (NG-MAST) genogroup, G12302 (132 [5·6%] of 2375; N gonorrhoeae sequence typing for antimicrobial resistance [NG-STAR] clonal complex [CC]168/63), carrying a mosaic mtrR promoter and mtrD sequence and found in 24 countries in 2018. CC63 was associated with pharyngeal infections in men who have sex with men. Susceptibility to ceftriaxone and cefixime is increasing, as the resistance-associated lineage, NG-MAST G1407 (51 [2·1%] of 2375), is progressively vanishing since 2009-10. INTERPRETATION: Enhanced gonococcal AMR surveillance is imperative worldwide. WGS, linked to epidemiological and AMR data, is essential to elucidate the dynamics in gonorrhoea epidemiology and gonococcal populations as well as to predict AMR. When feasible, WGS should supplement the national and international AMR surveillance programmes to elucidate AMR changes over time. In the EU and EEA, increasing low-level azithromycin resistance could threaten the recommended ceftriaxone-azithromycin dual therapy, and an evidence-based clinical azithromycin resistance breakpoint is needed. Nevertheless, increasing ceftriaxone susceptibility, declining cefixime resistance, and absence of known resistance mutations for new treatments (zoliflodacin, gepotidacin) are promising. FUNDING: European Centre for Disease Prevention and Control, Centre for Genomic Pathogen Surveillance, Örebro University Hospital, Wellcome.

Rapid Genomic Characterization and Global Surveillance of Klebsiella Using Pathogenwatch.

BACKGROUND: Klebsiella species, including the notable pathogen K. pneumoniae, are increasingly associated with antimicrobial resistance (AMR). Genome-based surveillance can inform interventions aimed at controlling AMR. However, its widespread implementation requires tools to streamline bioinformatic analyses and public health reporting. METHODS: We developed the web application Pathogenwatch, which implements analytics tailored to Klebsiella species for integration and visualization of genomic and epidemiological data. We populated Pathogenwatch with 16 537 public Klebsiella genomes to enable contextualization of user genomes. We demonstrated its features with 1636 genomes from 4 low- and middle-income countries (LMICs) participating in the NIHR Global Health Research Unit (GHRU) on AMR. RESULTS: Using Pathogenwatch, we found that GHRU genomes were dominated by a small number of epidemic drug-resistant clones of K. pneumoniae. However, differences in their distribution were observed (eg, ST258/512 dominated in Colombia, ST231 in India, ST307 in Nigeria, ST147 in the Philippines). Phylogenetic analyses including public genomes for contextualization enabled retrospective monitoring of their spread. In particular, we identified hospital outbreaks, detected introductions from abroad, and uncovered clonal expansions associated with resistance and virulence genes. Assessment of loci encoding O-antigens and capsule in K. pneumoniae, which represent possible vaccine candidates, showed that 3 O-types (O1-O3) represented 88.9% of all genomes, whereas capsule types were much more diverse. CONCLUSIONS: Pathogenwatch provides a free, accessible platform for real-time analysis of Klebsiella genomes to aid surveillance at local, national, and global levels. We have improved representation of genomes from GHRU participant countries, further facilitating ongoing surveillance.

Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool.

The response of the global virus genomics community to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has been unprecedented, with significant advances made towards the 'real-time' generation and sharing of SARS-CoV-2 genomic data. The rapid growth in virus genome data production has necessitated the development of new analytical methods that can deal with orders of magnitude of more genomes than previously available. Here, we present and describe Phylogenetic Assignment of Named Global Outbreak Lineages (pangolin), a computational tool that has been developed to assign the most likely lineage to a given SARS-CoV-2 genome sequence according to the Pango dynamic lineage nomenclature scheme. To date, nearly two million virus genomes have been submitted to the web-application implementation of pangolin, which has facilitated the SARS-CoV-2 genomic epidemiology and provided researchers with access to actionable information about the pandemic's transmission lineages.

A global resource for genomic predictions of antimicrobial resistance and surveillance of Salmonella Typhi at pathogenwatch.

As whole-genome sequencing capacity becomes increasingly decentralized, there is a growing opportunity for collaboration and the sharing of surveillance data within and between countries to inform typhoid control policies. This vision requires free, community-driven tools that facilitate access to genomic data for public health on a global scale. Here we present the Pathogenwatch scheme for Salmonella enterica serovar Typhi (S. Typhi), a web application enabling the rapid identification of genomic markers of antimicrobial resistance (AMR) and contextualization with public genomic data. We show that the clustering of S. Typhi genomes in Pathogenwatch is comparable to established bioinformatics methods, and that genomic predictions of AMR are highly concordant with phenotypic susceptibility data. We demonstrate the public health utility of Pathogenwatch with examples selected from >4,300 public genomes available in the application. Pathogenwatch provides an intuitive entry point to monitor of the emergence and spread of S. Typhi high risk clones.

A community-driven resource for genomic epidemiology and antimicrobial resistance prediction of Neisseria gonorrhoeae at Pathogenwatch.

BACKGROUND: Antimicrobial-resistant (AMR) Neisseria gonorrhoeae is an urgent threat to public health, as strains resistant to at least one of the two last-line antibiotics used in empiric therapy of gonorrhoea, ceftriaxone and azithromycin, have spread internationally. Whole genome sequencing (WGS) data can be used to identify new AMR clones and transmission networks and inform the development of point-of-care tests for antimicrobial susceptibility, novel antimicrobials and vaccines. Community-driven tools that provide an easy access to and analysis of genomic and epidemiological data is the way forward for public health surveillance. METHODS: Here we present a public health-focussed scheme for genomic epidemiology of N. gonorrhoeae at Pathogenwatch ( https://pathogen.watch/ngonorrhoeae ). An international advisory group of experts in epidemiology, public health, genetics and genomics of N. gonorrhoeae was convened to inform on the utility of current and future analytics in the platform. We implement backwards compatibility with MLST, NG-MAST and NG-STAR typing schemes as well as an exhaustive library of genetic AMR determinants linked to a genotypic prediction of resistance to eight antibiotics. A collection of over 12,000 N. gonorrhoeae genome sequences from public archives has been quality-checked, assembled and made public together with available metadata for contextualization. RESULTS: AMR prediction from genome data revealed specificity values over 99% for azithromycin, ciprofloxacin and ceftriaxone and sensitivity values around 99% for benzylpenicillin and tetracycline. A case study using the Pathogenwatch collection of N. gonorrhoeae public genomes showed the global expansion of an azithromycin-resistant lineage carrying a mosaic mtr over at least the last 10 years, emphasising the power of Pathogenwatch to explore and evaluate genomic epidemiology questions of public health concern. CONCLUSIONS: The N. gonorrhoeae scheme in Pathogenwatch provides customised bioinformatic pipelines guided by expert opinion that can be adapted to public health agencies and departments with little expertise in bioinformatics and lower-resourced settings with internet connection but limited computational infrastructure. The advisory group will assess and identify ongoing public health needs in the field of gonorrhoea, particularly regarding gonococcal AMR, in order to further enhance utility with modified or new analytic methods.

Globetrotting strangles: the unbridled national and international transmission of Streptococcus equi between horses.

The equine disease strangles, which is characterized by the formation of abscesses in the lymph nodes of the head and neck, is one of the most frequently diagnosed infectious diseases of horses around the world. The causal agent, Streptococcus equi subspecies equi, establishes a persistent infection in approximately 10 % of animals that recover from the acute disease. Such 'carrier' animals appear healthy and are rarely identified during routine veterinary examinations pre-purchase or transit, but can transmit S. equi to naïve animals initiating new episodes of disease. Here, we report the analysis and visualization of phylogenomic and epidemiological data for 670 isolates of S. equi recovered from 19 different countries using a new core-genome multilocus sequence typing (cgMLST) web bioresource. Genetic relationships among all 670 S. equi isolates were determined at high resolution, revealing national and international transmission events that drive this endemic disease in horse populations throughout the world. Our data argue for the recognition of the international importance of strangles by the Office International des Épizooties to highlight the health, welfare and economic cost of this disease. The Pathogenwatch cgMLST web bioresource described herein is available for tailored genomic analysis of populations of S. equi and its close relative S. equi subspecies zooepidemicus that are recovered from horses and other animals, including humans, throughout the world. This article contains data hosted by Microreact.

Visualizing variation within Global Pneumococcal Sequence Clusters (GPSCs) and country population snapshots to contextualize pneumococcal isolates.

Knowledge of pneumococcal lineages, their geographic distribution and antibiotic resistance patterns, can give insights into global pneumococcal disease. We provide interactive bioinformatic outputs to explore such topics, aiming to increase dissemination of genomic insights to the wider community, without the need for specialist training. We prepared 12 country-specific phylogenetic snapshots, and international phylogenetic snapshots of 73 common Global Pneumococcal Sequence Clusters (GPSCs) previously defined using PopPUNK, and present them in Microreact. Gene presence and absence defined using Roary, and recombination profiles derived from Gubbins are presented in Phandango for each GPSC. Temporal phylogenetic signal was assessed for each GPSC using BactDating. We provide examples of how such resources can be used. In our example use of a country-specific phylogenetic snapshot we determined that serotype 14 was observed in nine unrelated genetic backgrounds in South Africa. The international phylogenetic snapshot of GPSC9, in which most serotype 14 isolates from South Africa were observed, highlights that there were three independent sub-clusters represented by South African serotype 14 isolates. We estimated from the GPSC9-dated tree that the sub-clusters were each established in South Africa during the 1980s. We show how recombination plots allowed the identification of a 20 kb recombination spanning the capsular polysaccharide locus within GPSC97. This was consistent with a switch from serotype 6A to 19A estimated to have occured in the 1990s from the GPSC97-dated tree. Plots of gene presence/absence of resistance genes (tet, erm, cat) across the GPSC23 phylogeny were consistent with acquisition of a composite transposon. We estimated from the GPSC23-dated tree that the acquisition occurred between 1953 and 1975. Finally, we demonstrate the assignment of GPSC31 to 17 externally generated pneumococcal serotype 1 assemblies from Utah via Pathogenwatch. Most of the Utah isolates clustered within GPSC31 in a USA-specific clade with the most recent common ancestor estimated between 1958 and 1981. The resources we have provided can be used to explore to data, test hypothesis and generate new hypotheses. The accessible assignment of GPSCs allows others to contextualize their own collections beyond the data presented here.

Public health surveillance of multidrug-resistant clones of Neisseria gonorrhoeae in Europe: a genomic survey.

BACKGROUND: Traditional methods for molecular epidemiology of Neisseria gonorrhoeae are suboptimal. Whole-genome sequencing (WGS) offers ideal resolution to describe population dynamics and to predict and infer transmission of antimicrobial resistance, and can enhance infection control through linkage with epidemiological data. We used WGS, in conjunction with linked epidemiological and phenotypic data, to describe the gonococcal population in 20 European countries. We aimed to detail changes in phenotypic antimicrobial resistance levels (and the reasons for these changes) and strain distribution (with a focus on antimicrobial resistance strains in risk groups), and to predict antimicrobial resistance from WGS data. METHODS: We carried out an observational study, in which we sequenced isolates taken from patients with gonorrhoea from the European Gonococcal Antimicrobial Surveillance Programme in 20 countries from September to November, 2013. We also developed a web platform that we used for automated antimicrobial resistance prediction, molecular typing (N gonorrhoeae multi-antigen sequence typing [NG-MAST] and multilocus sequence typing), and phylogenetic clustering in conjunction with epidemiological and phenotypic data. FINDINGS: The multidrug-resistant NG-MAST genogroup G1407 was predominant and accounted for the most cephalosporin resistance, but the prevalence of this genogroup decreased from 248 (23%) of 1066 isolates in a previous study from 2009-10 to 174 (17%) of 1054 isolates in this survey in 2013. This genogroup previously showed an association with men who have sex with men, but changed to an association with heterosexual people (odds ratio=4·29). WGS provided substantially improved resolution and accuracy over NG-MAST and multilocus sequence typing, predicted antimicrobial resistance relatively well, and identified discrepant isolates, mixed infections or contaminants, and multidrug-resistant clades linked to risk groups. INTERPRETATION: To our knowledge, we provide the first use of joint analysis of WGS and epidemiological data in an international programme for regional surveillance of sexually transmitted infections. WGS provided enhanced understanding of the distribution of antimicrobial resistance clones, including replacement with clones that were more susceptible to antimicrobials, in several risk groups nationally and regionally. We provide a framework for genomic surveillance of gonococci through standardised sampling, use of WGS, and a shared information architecture for interpretation and dissemination by use of open access software. FUNDING: The European Centre for Disease Prevention and Control, The Centre for Genomic Pathogen Surveillance, Örebro University Hospital, and Wellcome.

Whole-Genome Sequencing for Routine Pathogen Surveillance in Public Health: a Population Snapshot of Invasive Staphylococcus aureus in Europe.

UNLABELLED: The implementation of routine whole-genome sequencing (WGS) promises to transform our ability to monitor the emergence and spread of bacterial pathogens. Here we combined WGS data from 308 invasive Staphylococcus aureus isolates corresponding to a pan-European population snapshot, with epidemiological and resistance data. Geospatial visualization of the data is made possible by a generic software tool designed for public health purposes that is available at the project URL (http://www.microreact.org/project/EkUvg9uY?tt=rc). Our analysis demonstrates that high-risk clones can be identified on the basis of population level properties such as clonal relatedness, abundance, and spatial structuring and by inferring virulence and resistance properties on the basis of gene content. We also show that in silico predictions of antibiotic resistance profiles are at least as reliable as phenotypic testing. We argue that this work provides a comprehensive road map illustrating the three vital components for future molecular epidemiological surveillance: (i) large-scale structured surveys, (ii) WGS, and (iii) community-oriented database infrastructure and analysis tools. IMPORTANCE: The spread of antibiotic-resistant bacteria is a public health emergency of global concern, threatening medical intervention at every level of health care delivery. Several recent studies have demonstrated the promise of routine whole-genome sequencing (WGS) of bacterial pathogens for epidemiological surveillance, outbreak detection, and infection control. However, as this technology becomes more widely adopted, the key challenges of generating representative national and international data sets and the development of bioinformatic tools to manage and interpret the data become increasingly pertinent. This study provides a road map for the integration of WGS data into routine pathogen surveillance. We emphasize the importance of large-scale routine surveys to provide the population context for more targeted or localized investigation and the development of open-access bioinformatic tools to provide the means to combine and compare independently generated data with publicly available data sets.

Microreact: visualizing and sharing data for genomic epidemiology and phylogeography.

Visualization is frequently used to aid our interpretation of complex datasets. Within microbial genomics, visualizing the relationships between multiple genomes as a tree provides a framework onto which associated data (geographical, temporal, phenotypic and epidemiological) are added to generate hypotheses and to explore the dynamics of the system under investigation. Selected static images are then used within publications to highlight the key findings to a wider audience. However, these images are a very inadequate way of exploring and interpreting the richness of the data. There is, therefore, a need for flexible, interactive software that presents the population genomic outputs and associated data in a user-friendly manner for a wide range of end users, from trained bioinformaticians to front-line epidemiologists and health workers. Here, we present Microreact, a web application for the easy visualization of datasets consisting of any combination of trees, geographical, temporal and associated metadata. Data files can be uploaded to Microreact directly via the web browser or by linking to their location (e.g. from Google Drive/Dropbox or via API), and an integrated visualization via trees, maps, timelines and tables provides interactive querying of the data. The visualization can be shared as a permanent web link among collaborators, or embedded within publications to enable readers to explore and download the data. Microreact can act as an end point for any tool or bioinformatic pipeline that ultimately generates a tree, and provides a simple, yet powerful, visualization method that will aid research and discovery and the open sharing of datasets.

Target selection for structural genomics of infectious diseases.

This chapter describes the protocols used to identify, filter, and annotate potential protein targets from an organism associated with infectious diseases. Protocols often combine computational approaches for mining information in public databases or for checking whether the protein has already been targeted for structure determination, with manual strategies that examine the literature for information on the biological role of the protein or the experimental strategies that explore the effects of knocking out the protein. Publicly available computational tools have been cited as much as possible. Where these do not exist, the concepts underlying in-house tools developed for the Center for Structural Genomics of Infectious Diseases have been described.

Crystal structures of the human Dysferlin inner DysF domain.

BACKGROUND: Mutations in dysferlin, the first protein linked with the cell membrane repair mechanism, causes a group of muscular dystrophies called dysferlinopathies. Dysferlin is a type two-anchored membrane protein, with a single C terminal trans-membrane helix, and most of the protein lying in cytoplasm. Dysferlin contains several C2 domains and two DysF domains which are nested one inside the other. Many pathogenic point mutations fall in the DysF domain region. RESULTS: We describe the crystal structure of the human dysferlin inner DysF domain with a resolution of 1.9 Ångstroms. Most of the pathogenic mutations are part of aromatic/arginine stacks that hold the domain in a folded conformation. The high resolution of the structure show that these interactions are a mixture of parallel ring/guanadinium stacking, perpendicular H bond stacking and aliphatic chain packing. CONCLUSIONS: The high resolution structure of the Dysferlin DysF domain gives a template on which to interpret in detail the pathogenic mutations that lead to disease.

Gene3D: Multi-domain annotations for protein sequence and comparative genome analysis.

Gene3D (http://gene3d.biochem.ucl.ac.uk) is a database of protein domain structure annotations for protein sequences. Domains are predicted using a library of profile HMMs from 2738 CATH superfamilies. Gene3D assigns domain annotations to Ensembl and UniProt sequence sets including >6000 cellular genomes and >20 million unique protein sequences. This represents an increase of 45% in the number of protein sequences since our last publication. Thanks to improvements in the underlying data and pipeline, we see large increases in the domain coverage of sequences. We have expanded this coverage by integrating Pfam and SUPERFAMILY domain annotations, and we now resolve domain overlaps to provide highly comprehensive composite multi-domain architectures. To make these data more accessible for comparative genome analyses, we have developed novel search algorithms for searching genomes to identify related multi-domain architectures. In addition to providing domain family annotations, we have now developed a pipeline for 3D homology modelling of domains in Gene3D. This has been applied to the human genome and will be rolled out to other major organisms over the next year.

Evolutionary history of the TBP-domain superfamily.

The TATA binding protein (TBP) is an essential transcription initiation factor in Archaea and Eucarya. Bacteria lack TBP, and instead use sigma factors for transcription initiation. TBP has a symmetric structure comprising two repeated TBP domains. Using sequence, structural and phylogenetic analyses, we examine the distribution and evolutionary history of the TBP domain, a member of the helix-grip fold family. Our analyses reveal a broader distribution than for TBP, with TBP-domains being present across all three domains of life. In contrast to TBP, all other characterized examples of the TBP domain are present as single copies, primarily within multidomain proteins. The presence of the TBP domain in the ubiquitous DNA glycosylases suggests that this fold traces back to the ancestor of all three domains of life. The TBP domain is also found in RNase HIII, and phylogenetic analyses show that RNase HIII has evolved from bacterial RNase HII via TBP-domain fusion. Finally, our comparative genomic screens confirm and extend earlier reports of proteins consisting of a single TBP domain among some Archaea. These monopartite TBP-domain proteins suggest that this domain is functional in its own right, and that the TBP domain could have first evolved as an independent protein, which was later recruited in different contexts.

New functional families (FunFams) in CATH to improve the mapping of conserved functional sites to 3D structures.

CATH version 3.5 (Class, Architecture, Topology, Homology, available at http://www.cathdb.info/) contains 173 536 domains, 2626 homologous superfamilies and 1313 fold groups. When focusing on structural genomics (SG) structures, we observe that the number of new folds for CATH v3.5 is slightly less than for previous releases, and this observation suggests that we may now know the majority of folds that are easily accessible to structure determination. We have improved the accuracy of our functional family (FunFams) sub-classification method and the CATH sequence domain search facility has been extended to provide FunFam annotations for each domain. The CATH website has been redesigned. We have improved the display of functional data and of conserved sequence features associated with FunFams within each CATH superfamily.