Direct Whole-Genome Sequencing of Sputum Accurately Identifies Drug-Resistant Mycobacterium tuberculosis Faster than MGIT Culture Sequencing.

The current methods available to diagnose antimicrobial-resistant Mycobacterium tuberculosis infections require a positive culture or only test a limited number of resistance-associated mutations. A rapid accurate identification of antimicrobial resistance enables the prompt initiation of effective treatment. Here, we determine the utility of whole-genome sequencing (WGS) of M. tuberculosis directly from routinely obtained diagnostic sputum samples to provide a comprehensive resistance profile compared to that from mycobacterial growth indicator tube (MGIT) WGS. We sequenced M. tuberculosis from 43 sputum samples by targeted DNA enrichment using the Agilent SureSelectXT kit, and 43 MGIT positive samples from each participant. Thirty two (74%) sputum samples and 43 (100%) MGIT samples generated whole genomes. The times to antimicrobial resistance profiles and concordance were compared with Xpert MTB/RIF and phenotypic resistance testing from cultures of the same samples. Antibiotic susceptibility could be predicted from WGS of sputum within 5 days of sample receipt and up to 24 days earlier than WGS from MGIT culture and up to 31 days earlier than phenotypic testing. Direct sputum results could be reduced to 3 days with faster hybridization and if only regions encoding drug resistance are sequenced. We show that direct sputum sequencing has the potential to provide comprehensive resistance detection significantly faster than MGIT whole-genome sequencing or phenotypic testing of resistance from cultures in a clinical setting. This improved turnaround time enables prompt appropriate treatment with associated patient and health service benefits. Improvements in sample preparation are necessary to ensure comparable sensitivities and complete resistance profile predictions in all cases.

Rapid Whole-Genome Sequencing of Mycobacterium tuberculosis Isolates Directly from Clinical Samples.

The rapid identification of antimicrobial resistance is essential for effective treatment of highly resistant Mycobacterium tuberculosis. Whole-genome sequencing provides comprehensive data on resistance mutations and strain typing for monitoring transmission, but unlike for conventional molecular tests, this has previously been achievable only from cultures of M. tuberculosis. Here we describe a method utilizing biotinylated RNA baits designed specifically for M. tuberculosis DNA to capture full M. tuberculosis genomes directly from infected sputum samples, allowing whole-genome sequencing without the requirement of culture. This was carried out on 24 smear-positive sputum samples, collected from the United Kingdom and Lithuania where a matched culture sample was available, and 2 samples that had failed to grow in culture. M. tuberculosis sequencing data were obtained directly from all 24 smear-positive culture-positive sputa, of which 20 were of high quality (>20× depth and >90% of the genome covered). Results were compared with those of conventional molecular and culture-based methods, and high levels of concordance between phenotypical resistance and predicted resistance based on genotype were observed. High-quality sequence data were obtained from one smear-positive culture-negative case. This study demonstrated for the first time the successful and accurate sequencing of M. tuberculosis genomes directly from uncultured sputa. Identification of known resistance mutations within a week of sample receipt offers the prospect for personalized rather than empirical treatment of drug-resistant tuberculosis, including the use of antimicrobial-sparing regimens, leading to improved outcomes.

Whole-genome enrichment and sequencing of Chlamydia trachomatis directly from clinical samples.

BACKGROUND: Chlamydia trachomatis is a pathogen of worldwide importance, causing more than 100 million cases of sexually transmitted infections annually. Whole-genome sequencing is a powerful high resolution tool that can be used to generate accurate data on bacterial population structure, phylogeography and mutations associated with antimicrobial resistance. The objective of this study was to perform whole-genome enrichment and sequencing of C. trachomatis directly from clinical samples. METHODS: C. trachomatis positive samples comprising seven vaginal swabs and three urine samples were sequenced without prior in vitro culture in addition to nine cultured C. trachomatis samples, representing different serovars. A custom capture RNA bait set, that captures all known diversity amongst C. trachomatis genomes, was used in a whole-genome enrichment step during library preparation to enrich for C. trachomatis DNA. All samples were sequenced on the MiSeq platform. RESULTS: Full length C. trachomatis genomes (>95-100% coverage of a reference genome) were successfully generated for eight of ten clinical samples and for all cultured samples. The proportion of reads mapping to C. trachomatis and the mean read depth across each genome were strongly linked to the number of bacterial copies within the original sample. Phylogenetic analysis confirmed the known population structure and the data showed potential for identification of minority variants and mutations associated with antimicrobial resistance. The sensitivity of the method was >10-fold higher than other reported methodologies. CONCLUSIONS: The combination of whole-genome enrichment and deep sequencing has proven to be a non-mutagenic approach, capturing all known variation found within C. trachomatis genomes. The method is a consistent and sensitive tool that enables rapid whole-genome sequencing of C. trachomatis directly from clinical samples and has the potential to be adapted to other pathogens with a similar clonal nature.

Clinical and analytical validation of an 82-gene comprehensive genome-profiling panel for identifying and interpreting variants responsible for inherited retinal dystrophies.

Inherited retinal dystrophies comprise a clinically complex and heterogenous group of diseases characterized by visual impairment due to pathogenic variants of over 300 different genes. Accurately identifying the causative gene and associated variant is crucial for the definitive diagnosis and subsequent selection of precise treatments. Consequently, well-validated genetic tests are required in the clinical practice. Here, we report the analytical and clinical validation of a next-generation sequencing targeted gene panel, the PrismGuide IRD Panel System. This system enables comprehensive genome profiling of 82 genes related to inherited retinal dystrophies. The PrismGuide IRD Panel System demonstrated 100% (n = 43) concordance with Sanger sequencing in detecting single-nucleotide variants, small insertions, and small deletions in the target genes and also in assessing their zygosity. It also identified copy-number loss in four out of five cases. When assessing precision, we evaluated the reproducibility of variant detection with 2,160 variants in 144 replicates and found 100% agreement in terms of single-nucleotide variants (n = 1,584) and small insertions and deletions (n = 576). Furthermore, the PrismGuide IRD Panel System generated sufficient read depth for variant calls across the purine-rich and highly repetitive open-reading frame 15 region of RPGR and detected all five variants tested. These results show that the PrismGuide IRD Panel System can accurately and consistently detect single-nucleotide variants and small insertions and deletions. Thus, the PrismGuide IRD Panel System could serve as useful tool that is applicable in clinical practice for identifying the causative genes based on the detection and interpretation of variants in patients with inherited retinal dystrophies and can contribute to a precise molecular diagnosis and targeted treatments.

The genome of a pathogenic rhodococcus: cooptive virulence underpinned by key gene acquisitions.

We report the genome of the facultative intracellular parasite Rhodococcus equi, the only animal pathogen within the biotechnologically important actinobacterial genus Rhodococcus. The 5.0-Mb R. equi 103S genome is significantly smaller than those of environmental rhodococci. This is due to genome expansion in nonpathogenic species, via a linear gain of paralogous genes and an accelerated genetic flux, rather than reductive evolution in R. equi. The 103S genome lacks the extensive catabolic and secondary metabolic complement of environmental rhodococci, and it displays unique adaptations for host colonization and competition in the short-chain fatty acid-rich intestine and manure of herbivores--two main R. equi reservoirs. Except for a few horizontally acquired (HGT) pathogenicity loci, including a cytoadhesive pilus determinant (rpl) and the virulence plasmid vap pathogenicity island (PAI) required for intramacrophage survival, most of the potential virulence-associated genes identified in R. equi are conserved in environmental rhodococci or have homologs in nonpathogenic Actinobacteria. This suggests a mechanism of virulence evolution based on the cooption of existing core actinobacterial traits, triggered by key host niche-adaptive HGT events. We tested this hypothesis by investigating R. equi virulence plasmid-chromosome crosstalk, by global transcription profiling and expression network analysis. Two chromosomal genes conserved in environmental rhodococci, encoding putative chorismate mutase and anthranilate synthase enzymes involved in aromatic amino acid biosynthesis, were strongly coregulated with vap PAI virulence genes and required for optimal proliferation in macrophages. The regulatory integration of chromosomal metabolic genes under the control of the HGT-acquired plasmid PAI is thus an important element in the cooptive virulence of R. equi.

Islands of linkage in an ocean of pervasive recombination reveals two-speed evolution of human cytomegalovirus genomes.

Human cytomegalovirus (HCMV) infects most of the population worldwide, persisting throughout the host's life in a latent state with periodic episodes of reactivation. While typically asymptomatic, HCMV can cause fatal disease among congenitally infected infants and immunocompromised patients. These clinical issues are compounded by the emergence of antiviral resistance and the absence of an effective vaccine, the development of which is likely complicated by the numerous immune evasins encoded by HCMV to counter the host's adaptive immune responses, a feature that facilitates frequent super-infections. Understanding the evolutionary dynamics of HCMV is essential for the development of effective new drugs and vaccines. By comparing viral genomes from uncultivated or low-passaged clinical samples of diverse origins, we observe evidence of frequent homologous recombination events, both recent and ancient, and no structure of HCMV genetic diversity at the whole-genome scale. Analysis of individual gene-scale loci reveals a striking dichotomy: while most of the genome is highly conserved, recombines essentially freely and has evolved under purifying selection, 21 genes display extreme diversity, structured into distinct genotypes that do not recombine with each other. Most of these hyper-variable genes encode glycoproteins involved in cell entry or escape of host immunity. Evidence that half of them have diverged through episodes of intense positive selection suggests that rapid evolution of hyper-variable loci is likely driven by interactions with host immunity. It appears that this process is enabled by recombination unlinking hyper-variable loci from strongly constrained neighboring sites. It is conceivable that viral mechanisms facilitating super-infection have evolved to promote recombination between diverged genotypes, allowing the virus to continuously diversify at key loci to escape immune detection, while maintaining a genome optimally adapted to its asymptomatic infectious lifecycle.

Development and application of versatile high density microarrays for genome-wide analysis of Streptomyces coelicolor: characterization of the HspR regulon.

BACKGROUND: DNA microarrays are a key resource for global analysis of genome content, gene expression and the distribution of transcription factor binding sites. We describe the development and application of versatile high density ink-jet in situ-synthesized DNA arrays for the G+C rich bacterium Streptomyces coelicolor. High G+C content DNA probes often perform poorly on arrays, yielding either weak hybridization or non-specific signals. Thus, more than one million 60-mer oligonucleotide probes were experimentally tested for sensitivity and specificity to enable selection of optimal probe sets for the genome microarrays. The heat-shock HspR regulatory system of S. coelicolor, a well-characterized repressor with a small number of known targets, was exploited to test and validate the arrays for use in global chromatin immunoprecipitation-on-chip (ChIP-chip) and gene expression analysis. RESULTS: In addition to confirming dnaK, clpB and lon as in vivo targets of HspR, it was revealed, using a novel ChIP-chip data clustering method, that HspR also apparently interacts with ribosomal RNA (rrnD operon) and specific transfer RNA genes (the tRNAGln/tRNAGlu cluster). It is suggested that enhanced synthesis of Glu-tRNAGlu may reflect increased demand for tetrapyrrole biosynthesis following heat-shock. Moreover, it was found that heat-shock-induced genes are significantly enriched for Gln/Glu codons relative to the whole genome, a finding that would be consistent with HspR-mediated control of the tRNA species. CONCLUSIONS: This study suggests that HspR fulfils a broader, unprecedented role in adaptation to stresses than previously recognized -- influencing expression of key components of the translational apparatus in addition to molecular chaperone and protease-encoding genes. It is envisaged that these experimentally optimized arrays will provide a key resource for systems level studies of Streptomyces biology.

Studies of the distribution of Escherichia coli cAMP-receptor protein and RNA polymerase along the E. coli chromosome.

Chromatin immunoprecipitation and high-density microarrays have been used to monitor the distribution of the global transcription regulator Escherichia coli cAMP-receptor protein (CRP) and RNA polymerase along the E. coli chromosome. Our results identify targets occupied by CRP and genes transcribed by RNA polymerase in vivo. Many of the loci of CRP binding are at known CRP regulated promoters. However, our results show that CRP also interacts with thousands of weaker sites across the whole chromosome and that this "background" binding can be used as a probe for organization within the E. coli folded chromosome. In rapidly growing cells, we show that the major sites of RNA polymerase binding are approximately 90 transcription units that include genes needed for protein synthesis. Upon the addition of rifampicin, RNA polymerase is distributed among >500 functional promoters. We show that the chromatin immunoprecipitation and high-density-microarrays methodology can be used to study the redistribution of RNA polymerase induced by environmental stress, revealing previously uncharacterized aspects of RNA polymerase behavior and providing an alternative to the "transcriptomics" approach for studying global transcription patterns.

Transgenes encompassing dual-promoter CpG islands from the human TBP and HNRPA2B1 loci are resistant to heterochromatin-mediated silencing.

The genetic elements that are responsible for establishing a transcriptionally competent, open chromatin structure at a region of the genome that consists only of ubiquitously expressed, housekeeping genes are currently unknown. We demonstrate for the first time through functional analysis in stably transfected tissue culture cells that transgenes containing methylation-free CpG islands spanning the dual divergently transcribed promoters from the human TATA binding protein (TBP)-proteasome component-B1 (PSMB1) and heterogeneous nuclear ribonucleoprotein A2/B1 (HNRPA2B1)-heterochromatin protein 1Hs-gamma (chromobox homolog 3, CBX3) gene loci are sufficient to prevent transcriptional silencing and a variegated expression pattern when integrated within centromeric heterochromatin. In addition, only transgene constructs extending over both the HNRPA2B1 and the CBX3 promoters, and not the HNRPA2B1 promoter alone, were able to confer high and stable long-term EGFP reporter gene expression. These observations suggest that methylation-free CpG islands associated with dual, divergently transcribed promoters possess an independent dominant chromatin opening function and may therefore be major determinants in establishing and maintaining a region of open chromatin at housekeeping gene loci.