MolBioLib: A C++11 Framework for Rapid Development and Deployment of Bioinformatics Tasks

Summary: We developed MolBioLib to address the need for adaptable next-generation sequencing analysis tools. The result is a compact, portable and extensively tested C++11 software framework and set of applications tailored to the demands of next-generation sequencing data and applicable to many other applications. MolBioLib is designed to work with common file formats and data types used both in genomic analysis and general data analysis. A central relational-database-like Table class is a flexible and powerful object to intuitively represent and work with a wide variety of tabular datasets, ranging from alignment data to annotations. MolBioLib has been used to identify causative single-nucleotide polymorphisms in whole genome sequencing, detect balanced chromosomal rearrangements and compute enrichment of messenger RNAs (mRNAs) on microtubules, typically requiring applications of under 200 lines of code. MolBioLib includes programs to perform a wide variety of analysis tasks, such as computing read coverage, annotating genomic intervals and novel peak calling with a wavelet algorithm. Although MolBioLib was designed primarily for bioinformatics purposes, much of its functionality is applicable to a wide range of problems. Complete documentation and an extensive automated test suite are provided

Atypical Case Of Wolfram Syndrome Revealed Through Targeted Exome Sequencing In A Patient With Suspected Mitochondrial Disease

Background: Mitochondrial diseases comprise a diverse set of clinical disorders that affect multiple organ systems with varying severity and age of onset. Due to their clinical and genetic heterogeneity, these diseases are difficult to diagnose. We have developed a targeted exome sequencing approach to improve our ability to properly diagnose mitochondrial diseases and apply it here to an individual patient. Our method targets mitochondrial DNA (mtDNA) and the exons of 1,600 nuclear genes involved in mitochondrial biology or Mendelian disorders with multi-system phenotypes, thereby allowing for simultaneous evaluation of multiple disease loci. Case Presentation: Targeted exome sequencing was performed on a patient initially suspected to have a mitochondrial disorder. The patient presented with diabetes mellitus, diffuse brain atrophy, autonomic neuropathy, optic nerve atrophy, and a severe amnestic syndrome. Further work-up revealed multiple heteroplasmic mtDNA deletions as well as profound thiamine deficiency without a clear nutritional cause. Targeted exome sequencing revealed a homozygous c.1672C > T (p.R558C) missense mutation in exon 8 of WFS1 that has previously been reported in a patient with Wolfram syndrome. Conclusion: This case demonstrates how clinical application of next-generation sequencing technology can enhance the diagnosis of patients suspected to have rare genetic disorders. Furthermore, the finding of unexplained thiamine deficiency in a patient with Wolfram syndrome suggests a potential link between WFS1 biology and thiamine metabolism that has implications for the clinical management of Wolfram syndrome patients

RNA signatures allow rapid identification of pathogens and antibiotic susceptibilities

With rising rates of drug-resistant infections, there is a need for diagnostic methods that rapidly can detect the presence of pathogens and reveal their susceptibility to antibiotics. Here we propose an approach to diagnosing the presence and drug-susceptibility of infectious diseases based on direct detection of RNA from clinical samples. We demonstrate that species-specific RNA signatures can be used to identify a broad spectrum of infectious agents, including bacteria, viruses, yeast, and parasites. Moreover, we show that the behavior of a small set of bacterial transcripts after a brief antibiotic pulse can rapidly differentiate drug-susceptible and -resistant organisms and that these measurements can be made directly from clinical materials. Thus, transcriptional signatures could form the basis of a uniform diagnostic platform applicable across a broad range of infectious agents

Complex Reorganization and Predominant Non-Homologous Repair Following Chromosomal Breakage in Karyotypically Balanced Germline Rearrangements and Transgenic Integration

We defined the genetic landscape of balanced chromosomal rearrangements at nucleotide resolution by sequencing 141 breakpoints from cytogenetically-interpreted translocations and inversions. We confirm that the recently described phenomenon of “chromothripsis” (massive chromosomal shattering and reorganization) is not unique to cancer cells but also occurs in the germline where it can resolve to a karyotypically balanced state with frequent inversions. We detected a high incidence of complex rearrangements (19.2%) and substantially less reliance on microhomology (31%) than previously observed in benign CNVs. We compared these results to experimentally-generated DNA breakage-repair by sequencing seven transgenic animals, and revealed extensive rearrangement of the transgene and host genome with similar complexity to human germline alterations. Inversion is the most common rearrangement, suggesting that a combined mechanism involving template switching and non-homologous repair mediates the formation of balanced complex rearrangements that are viable, stably replicated and transmitted unaltered to subsequent generations

HIV Testing of At Risk Patients in a Large Integrated Health Care System

OBJECTIVE: Early identification of HIV infection is critical for patients to receive life-prolonging treatment and risk-reduction counseling. Understanding HIV screening practices and barriers to HIV testing is an important prelude to designing successful HIV screening programs. Our objective was to evaluate current practice patterns for identification of HIV. METHODS: We used a retrospective cohort analysis of 13,991 at-risk patients seen at 4 large Department of Veterans Affairs (VA) health-care systems. We also reviewed 1,100 medical records of tested patients. We assessed HIV testing rates among at-risk patients, the rationale for HIV testing, and predictors of HIV testing and of HIV infection. RESULTS: Of the 13,991 patients at risk for HIV, only 36% had been HIV-tested. The prevalence of HIV ranged from 1% to 20% among tested patients at the 4 sites. Approximately 90% of patients who were tested had a documented reason for testing. CONCLUSION: One-half to two-thirds of patients at risk for HIV had not been tested within our selected VA sites. Among tested patients, the rationale for HIV testing was well documented. Further testing of at-risk patients could clearly benefit patients who have unidentified HIV infection by providing earlier access to life-prolonging therapy

Oncogenic BRAF mutation induces DNA methylation changes in a murine model for human serrated colorectal neoplasia

Colorectal cancer is a major cause of cancer death and approximately 20% arises within serrated polyps, which are under-recognized and poorly understood. Human serrated colorectal polyps frequently exhibit both oncogenic BRAF mutation and widespread DNA methylation changes, which are important in silencing genes restraining neoplastic progression. Here, we investigated whether in vivo induction of mutant Braf is sufficient to result in coordinated promoter methylation changes for multiple cancer-related genes. The BrafV637E mutation was induced in murine intestine on an FVB;C57BL/6J background and assessed for morphological and DNA methylation changes at multiple time points from 10 days to 14 months. Extensive intestinal hyperplasia developed by 10 days post-induction of the mutation. By 8 months, most mice had murine serrated adenomas with dysplasia and invasive cancer developed in 40% of mice by 14 months. From 5 months onwards, Braf mutant mice showed extensive, gene-specific increases in DNA methylation even in hyperplastic mucosa without lesions. This demonstrates that persistent oncogenic Braf signaling is sufficient to induce widespread DNA methylation changes. This occurs over an extended period of time, mimicking the long latency followed by rapid progression of human serrated neoplasia. This study establishes for the first time that DNA methylation arises slowly in direct response to prolonged oncogenic Braf signaling in serrated polyps; this finding has implications both for chemoprevention and for understanding the origin of DNA hypermethylation in cancer generally.Catherine E. Bond, Cheng Liu, Futoshi Kawamata, Diane M. McKeone, Winnie Fernando, Saara Jamieson, Sally-Ann Pearson, Alexandra Kane, Susan L. Woods, Tamsin R.M. Lannagan, Roshini Somashekar, Young Lee, Troy Dumenil, Gunter Hartel, Kevin J. Spring, Jennifer Borowsky, Lochlan Fennell, Mark Bettington, Jason Lee, Daniel L. Worthley, Barbara A. Leggett and Vicki L.J. Whitehal

Comparative Genomic Characterization of Francisella tularensis Strains Belonging to Low and High Virulence Subspecies

Tularemia is a geographically widespread, severely debilitating, and occasionally lethal disease in humans. It is caused by infection by a gram-negative bacterium, Francisella tularensis. In order to better understand its potency as an etiological agent as well as its potential as a biological weapon, we have completed draft assemblies and report the first complete genomic characterization of five strains belonging to the following different Francisella subspecies (subsp.): the F. tularensis subsp. tularensis FSC033, F. tularensis subsp. holarctica FSC257 and FSC022, and F. tularensis subsp. novicida GA99-3548 and GA99-3549 strains. Here, we report the sequencing of these strains and comparative genomic analysis with recently available public Francisella sequences, including the rare F. tularensis subsp. mediasiatica FSC147 strain isolate from the Central Asian Region. We report evidence for the occurrence of large-scale rearrangement events in strains of the holarctica subspecies, supporting previous proposals that further phylogenetic subdivisions of the Type B clade are likely. We also find a significant enrichment of disrupted or absent ORFs proximal to predicted breakpoints in the FSC022 strain, including a genetic component of the Type I restriction-modification defense system. Many of the pseudogenes identified are also disrupted in the closely related rarely human pathogenic F. tularensis subsp. mediasiatica FSC147 strain, including modulator of drug activity B (mdaB) (FTT0961), which encodes a known NADPH quinone reductase involved in oxidative stress resistance. We have also identified genes exhibiting sequence similarity to effectors of the Type III (T3SS) and components of the Type IV secretion systems (T4SS). One of the genes, msrA2 (FTT1797c), is disrupted in F. tularensis subsp. mediasiatica and has recently been shown to mediate bacterial pathogen survival in host organisms. Our findings suggest that in addition to the duplication of the Francisella Pathogenicity Island, and acquisition of individual loci, adaptation by gene loss in the more recently emerged tularensis, holarctica, and mediasiatica subspecies occurred and was distinct from evolutionary events that differentiated these subspecies, and the novicida subspecies, from a common ancestor. Our findings are applicable to future studies focused on variations in Francisella subspecies pathogenesis, and of broader interest to studies of genomic pathoadaptation in bacteria