Anita C. Schürch

Combined Species Identification, Genotyping, and Drug Resistance Detection of Mycobacterium tuberculosis Cultures by MLPA on a Bead-Based Array

The population structure of Mycobacterium tuberculosis is typically clonal therefore genotypic lineages can be unequivocally identified by characteristic markers such as mutations or genomic deletions. In addition, drug resistance is mainly mediated by mutations. These issues make multiplexed detection of selected mutations potentially a very powerful tool to characterise Mycobacterium tuberculosis. We used Multiplex Ligation-dependent Probe Amplification (MLPA) to screen for dispersed mutations, which can be successfully applied to Mycobacterium tuberculosis as was previously shown. Here we selected 47 discriminative and informative markers and designed MLPA probes accordingly to allow analysis with a liquid bead array and robust reader (Luminex MAGPIX technology). To validate the bead-based MLPA, we screened a panel of 88 selected strains, previously characterised by other methods with the developed multiplex assay using automated positive and negative calling. In total 3059 characteristics were screened and 3034 (99.2%) were consistent with previous molecular characterizations, of which 2056 (67.2%) were directly supported by other molecular methods, and 978 (32.0%) were consistent with but not directly supported by previous molecular characterizations. Results directly conflicting or inconsistent with previous methods, were obtained for 25 (0.8%) of the characteristics tested. Here we report the validation of the bead-based MLPA and demonstrate its potential to simultaneously identify a range of drug resistance markers, discriminate the species within the Mycobacterium tuberculosis complex, determine the genetic lineage and detect and identify the clinically most relevant non-tuberculous mycobacterial species. The detection of multiple genetic markers in clinically derived Mycobacterium tuberculosis strains with a multiplex assay could reduce the number of TB-dedicated screening methods needed for full characterization. Additionally, as a proportion of the markers screened are specific to certain Mycobacterium tuberculosis lineages each profile can be checked for internal consistency. Strain characterization can allow selection of appropriate treatment and thereby improve treatment outcome and patient management.

CRISPR/Cas9-Mediated Genome Editing of Herpesviruses Limits Productive and Latent Infections

Herpesviruses infect the majority of the human population and can cause significant morbidity and mortality. Herpes simplex virus (HSV) type 1 causes cold sores and herpes simplex keratitis, whereas HSV-2 is responsible for genital herpes. Human cytomegalovirus (HCMV) is the most common viral cause of congenital defects and is responsible for serious disease in immuno-compromised individuals. Epstein-Barr virus (EBV) is associated with infectious mononucleosis and a broad range of malignancies, including Burkitt’s lymphoma, nasopharyngeal carcinoma, Hodgkin’s disease, and post-transplant lymphomas. Herpesviruses persist in their host for life by establishing a latent infection that is interrupted by periodic reactivation events during which replication occurs. Current antiviral drug treatments target the clinical manifestations of this productive stage, but they are ineffective at eliminating these viruses from the infected host. Here, we set out to combat both productive and latent herpesvirus infections by exploiting the CRISPR/Cas9 system to target viral genetic elements important for virus fitness. We show effective abrogation of HCMV and HSV-1 replication by targeting gRNAs to essential viral genes. Simultaneous targeting of HSV-1 with multiple gRNAs completely abolished the production of infectious particles from human cells. Using the same approach, EBV can be almost completely cleared from latently infected EBV-transformed human tumor cells. Our studies indicate that the CRISPR/Cas9 system can be effectively targeted to herpesvirus genomes as a potent prophylactic and therapeutic anti-viral strategy that may be used to impair viral replication and clear latent virus infection.

Exploring the Potential of Next-Generation Sequencing in Detection of Respiratory Viruses

ABSTRACT Efficient detection of human respiratory viral pathogens is crucial in the management of patients with acute respiratory tract infection. Sequence-independent amplification of nucleic acids combined with next-generation sequencing technology and bioinformatics analyses is a promising strategy for identifying pathogens in clinical and public health settings. It allows the characterization of hundreds of different known pathogens simultaneously and of novel pathogens that elude conventional testing. However, major hurdles for its routine use exist, including cost, turnaround time, and especially sensitivity of the assay, as the detection limit is dependent on viral load, host genetic material, and sequencing depth. To obtain insights into these aspects, we analyzed nasopharyngeal aspirates from a cohort of 81 Thai children with respiratory disease for the presence of respiratory viruses using a sequence-independent next-generation sequencing approach and routinely used diagnostic real-time reverse transcriptase PCR (real-time RT-PCR) assays. With respect to the detection of rhinovirus and human metapneumovirus, the next-generation sequencing approach was at least as sensitive as diagnostic real-time RT-PCR in this small cohort, whereas for bocavirus and enterovirus, next-generation sequencing was less sensitive than real-time RT-PCR. The advantage of the sequencing approach over real-time RT-PCR was the immediate availability of virus-typing information. Considering the development of platforms capable of generating more output data at declining costs, next-generation sequencing remains of interest for future virus diagnosis in clinical and public health settings and certainly as an additional tool when screening results from real-time RT-PCR are negative.

Comparative Analysis of Mycobacterium tuberculosis pe and ppe Genes Reveals High Sequence Variation and an Apparent Absence of Selective Constraints

Mycobacterium tuberculosis complex (MTBC) genomes contain 2 large gene families termed pe and ppe. The function of pe/ppe proteins remains enigmatic but studies suggest that they are secreted or cell surface associated and are involved in bacterial virulence. Previous studies have also shown that some pe/ppe genes are polymorphic, a finding that suggests involvement in antigenic variation. Using comparative sequence analysis of 18 publicly available MTBC whole genome sequences, we have performed alignments of 33 pe (excluding pe_pgrs) and 66 ppe genes in order to detect the frequency and nature of genetic variation. This work has been supplemented by whole gene sequencing of 14 pe/ppe (including 5 pe_pgrs) genes in a cohort of 40 diverse and well defined clinical isolates covering all the main lineages of the M. tuberculosis phylogenetic tree. We show that nsSNPs in pe (excluding pgrs) and ppe genes are 3.0 and 3.3 times higher than in non-pe/ppe genes respectively and that numerous other mutation types are also present at a high frequency. It has previously been shown that non-pe/ppe M. tuberculosis genes display a remarkably low level of purifying selection. Here, we also show that compared to these genes those of the pe/ppe families show a further reduction of selection pressure that suggests neutral evolution. This is inconsistent with the positive selection pressure of “classical” antigenic variation. Finally, by analyzing such a large number of genes we were able to detect large differences in mutation type and frequency between both individual genes and gene sub-families. The high variation rates and absence of selective constraints provides valuable insights into potential pe/ppe function. Since pe/ppe proteins are highly antigenic and have been studied as potential vaccine components these results should also prove informative for aspects of M. tuberculosis vaccine design.

Novel cyclovirus in human cerebrospinal fluid, Malawi, 2010-2011.

To identify unknown human viruses, we analyzed serum and cerebrospinal fluid samples from patients with unexplained paraplegia from Malawi by using viral metagenomics. A novel cyclovirus species was identified and subsequently found in 15% and 10% of serum and cerebrospinal fluid samples, respectively. These data expand our knowledge of cyclovirus diversity and tropism.

New viruses in idiopathic human diarrhea cases, the Netherlands

Emerging viral infections can be identified by using a viral metagenomics approach for clinical human material. Diarrhea samples of patients with unexplained gastroenteritis from the Netherlands were analyzed by using viral metagenomics. Novel circular DNA viruses, bufaviruses, and genogroup III picobirnaviruses were identified. These data expand our knowledge of the human virome.

Proteogenomic Analysis of Polymorphisms and Gene Annotation Divergences in Prokaryotes using a Clustered Mass Spectrometry-Friendly Database

Precise annotation of genes or open reading frames is still a difficult task that results in divergence even for data generated from the same genomic sequence. This has an impact in further proteomic studies, and also compromises the characterization of clinical isolates with many specific genetic variations that may not be represented in the selected database. We recently developed software called multistrain mass spectrometry prokaryotic database builder (MSMSpdbb) that can merge protein databases from several sources and be applied on any prokaryotic organism, in a proteomic-friendly approach. We generated a database for the Mycobacterium tuberculosis complex (using three strains of Mycobacterium bovis and five of M. tuberculosis), and analyzed data collected from two laboratory strains and two clinical isolates of M. tuberculosis. We identified 2561 proteins, of which 24 were present in M. tuberculosis H37Rv samples, but not annotated in the M. tuberculosis H37Rv genome. We were also able to identify 280 nonsynonymous single amino acid polymorphisms and confirm 367 translational start sites. As a proof of concept we applied the database to whole-genome DNA sequencing data of one of the clinical isolates, which allowed the validation of 116 predicted single amino acid polymorphisms and the annotation of 131 N-terminal start sites. Moreover we identified regions not present in the original M. tuberculosis H37Rv sequence, indicating strain divergence or errors in the reference sequence. In conclusion, we demonstrated the potential of using a merged database to better characterize laboratory or clinical bacterial strains.

Optimization of Standard In-House 24-Locus Variable-Number Tandem-Repeat Typing for Mycobacterium tuberculosis and Its Direct Application to Clinical Material

ABSTRACT Variable-number tandem-repeat (VNTR) typing with a panel of 24 loci is the current gold standard in the molecular typing of Mycobacterium tuberculosis complex isolates. However, because of technical problems, a part of the loci often cannot be amplified by multiplex PCRs. Therefore, a considerable number of single-locus PCRs have to be performed for the loci with missing results, which impairs the laboratory work flow. Therefore, the original in-house method described by Supply et al. in 2006 was reevaluated. We modified seven primers and the PCR master mixture and obtained a strongly optimized in-house 24-locus VNTR typing method. The percentage of instantly complete 24-locus VNTR patterns detected in the routine flow of typing activities increased to 84.7% from the 72.3% obtained with the typing conducted with the commercially available Genoscreen MIRU-VNTR typing kit. The analytical sensitivity of the optimized in-house method was assessed by serial dilutions of M. tuberculosis in bronchoalveolar lavage fluid. A 1:10 dilution of the different strains tested was the lowest dilution for the detection of a complete 24-locus VNTR pattern. The optimized in-house 24-locus VNTR typing method will reduce the turnaround time of typing significantly and also the financial burden of these activities.

On the (im)possibility of reconstructing plasmids from whole-genome short-read sequencing data

To benchmark algorithms for automated plasmid sequence reconstruction from short-read sequencing data, we selected 42 publicly available complete bacterial genome sequences spanning 12 genera, containing 148 plasmids. We predicted plasmids from short-read data with four programs (PlasmidSPAdes, Recycler, cBar and PlasmidFinder) and compared the outcome to the reference sequences. PlasmidSPAdes reconstructs plasmids based on coverage differences in the assembly graph. It reconstructed most of the reference plasmids (recall=0.82), but approximately a quarter of the predicted plasmid contigs were false positives (precision=0.75). PlasmidSPAdes merged 84 % of the predictions from genomes with multiple plasmids into a single bin. Recycler searches the assembly graph for sub-graphs corresponding to circular sequences and correctly predicted small plasmids, but failed with long plasmids (recall=0.12, precision=0.30). cBar, which applies pentamer frequency analysis to detect plasmid-derived contigs, showed a recall and precision of 0.76 and 0.62, respectively. However, cBar categorizes contigs as plasmid-derived and does not bin the different plasmids. PlasmidFinder, which searches for replicons, had the highest precision (1.0), but was restricted by the contents of its database and the contig length obtained from de novo assembly (recall=0.36). PlasmidSPAdes and Recycler detected putative small plasmids (<10 kbp), which were also predicted as plasmids by cBar, but were absent in the original assembly. This study shows that it is possible to automatically predict small plasmids. Prediction of large plasmids (>50 kbp) containing repeated sequences remains challenging and limits the high-throughput analysis of plasmids from short-read whole-genome sequencing data.

Recovering full-length viral genomes from metagenomes

Infectious disease metagenomics is driven by the question: “what is causing the disease?” in contrast to classical metagenome studies which are guided by “what is out there?” In case of a novel virus, a first step to eventually establishing etiology can be to recover a full-length viral genome from a metagenomic sample. However, retrieval of a full-length genome of a divergent virus is technically challenging and can be time-consuming and costly. Here we discuss different assembly and fragment linkage strategies such as iterative assembly, motif searches, k-mer frequency profiling, coverage profile binning, and other strategies used to recover genomes of potential viral pathogens in a timely and cost-effective manner.