Stool metatranscriptomics: A technical guideline for mRNA stabilisation and isolation.

BACKGROUND: The complex microbiome of the gut has an enormous impact on human health. Analysis of the transcriptional activity of microorganisms through mRNA sequencing (metatranscriptomics) opens a completely new window into their activity in vivo, but it is highly challenging due to numerous technical and bioinformatical obstacles. Here we present an optimized pipeline for extraction of high quality mRNA from stool samples. RESULTS: Comparison of three commercially available RNA extraction kits with the method of Zoetendal revealed that the Powermicrobiome Kit (MoBio) performed best with respect to RNA yield and purity. Next, the influence of the stabilization reagent during sample storage for up to 15 days was studied. RIN analysis and qRT-PCR of spiked-in and indigenous genes revealed that RNA Later preserved mRNA integrity most efficiently, while samples conserved in RNA Protect showed substantial mRNA decay. Using the optimized pipeline developed here, recovery rates for spiked-in E.coli cells expressing fluorescing proteins were 8.7-9.7% for SuperfolderGFP and 14.7-17.8% for mCherry. The mRNA of stabilized stool samples as well as of snap-frozen controls was sequenced with Illumina Hiseq, yielding on average 74 million reads per sample. PCoA analysis, taxonomic classification using Kraken and functional classification using bwa showed that the transcriptomes of samples conserved in RNA Later were unchanged for up to 6 days even at room temperature, while RNA Protect was inefficient for storage durations exceeding 24 h. However, our data indicate that RNA Later introduces a bias which is then maintained throughout storage, while RNA Protect conserved samples are initially more similar to the snap frozen controls. RNA Later conserved samples had a reduced abundance of e.g. Prevotellaceae transcripts and were depleted for e.g. COG category "Carbohydrate transport and metabolism". CONCLUSION: Since the overall similarity between all stool transcriptional profiles studied here was >0.92, these differences are unlikely to affect global comparisons, but should be taken into account when rare but critically important members of the stool microbiome are being studied.

r2cat: synteny plots and comparative assembly.

SUMMARY: Recent parallel pyrosequencing methods and the increasing number of finished genomes encourage the sequencing and investigation of closely related strains. Although the sequencing itself becomes easier and cheaper with each machine generation, the finishing of the genomes remains difficult. Instead of the desired whole genomic sequence, a set of contigs is the result of the assembly. In this applications note, we present the tool r2cat (related reference contig arrangement tool) that helps in the task of comparative assembly and also provides an interactive visualization for synteny inspection.

The complete genome sequence of Corynebacterium pseudotuberculosis FRC41 isolated from a 12-year-old girl with necrotizing lymphadenitis reveals insights into gene-regulatory networks contributing to virulence.

BACKGROUND: Corynebacterium pseudotuberculosis is generally regarded as an important animal pathogen that rarely infects humans. Clinical strains are occasionally recovered from human cases of lymphadenitis, such as C. pseudotuberculosis FRC41 that was isolated from the inguinal lymph node of a 12-year-old girl with necrotizing lymphadenitis. To detect potential virulence factors and corresponding gene-regulatory networks in this human isolate, the genome sequence of C. pseudotuberculosis FCR41 was determined by pyrosequencing and functionally annotated. RESULTS: Sequencing and assembly of the C. pseudotuberculosis FRC41 genome yielded a circular chromosome with a size of 2,337,913 bp and a mean G+C content of 52.2%. Specific gene sets associated with iron and zinc homeostasis were detected among the 2,110 predicted protein-coding regions and integrated into a gene-regulatory network that is linked with both the central metabolism and the oxidative stress response of FRC41. Two gene clusters encode proteins involved in the sortase-mediated polymerization of adhesive pili that can probably mediate the adherence to host tissue to facilitate additional ligand-receptor interactions and the delivery of virulence factors. The prominent virulence factors phospholipase D (Pld) and corynebacterial protease CP40 are encoded in the genome of this human isolate. The genome annotation revealed additional serine proteases, neuraminidase H, nitric oxide reductase, an invasion-associated protein, and acyl-CoA carboxylase subunits involved in mycolic acid biosynthesis as potential virulence factors. The cAMP-sensing transcription regulator GlxR plays a key role in controlling the expression of several genes contributing to virulence. CONCLUSION: The functional data deduced from the genome sequencing and the extended knowledge of virulence factors indicate that the human isolate C. pseudotuberculosis FRC41 is equipped with a distinct gene set promoting its survival under unfavorable environmental conditions encountered in the mammalian host.

Next-generation-sequencing of recurrent childhood high hyperdiploid acute lymphoblastic leukemia reveals mutations typically associated with high risk patients.

20% of children suffering from high hyperdiploid acute lymphoblastic leukemia develop recurrent disease. The molecular mechanisms are largely unknown. Here, we analyzed the genetic landscape of five patients at relapse, who developed recurrent disease without prior high-risk indication using whole-exome- and whole-genome-sequencing. Oncogenic mutations of RAS pathway genes (NRAS, KRAS, FLT3, n=4) and deactivating mutations of major epigenetic regulators (CREBBP, EP300, each n=2 and ARID4B, EZH2, MACROD2, MLL2, each n=1) were prominent in these cases and virtually absent in non-recurrent cases (n=6) or other pediatric acute lymphoblastic leukemia cases (n=18). In relapse nucleotide variations were detected in cell fate determining transcription factors (GLIS1, AKNA). Structural genomic alterations affected genes regulating B-cell development (IKZF1, PBX1, RUNX1). Eleven novel translocations involved the genes ART4, C12orf60, MACROD2, TBL1XR1, LRRN4, KIAA1467, and ELMO1/MIR1200. Typically, patients harbored only single structural variations, except for one patient who displayed massive rearrangements in the context of a germline tumor suppressor TP53 mutation and a Li-Fraumeni syndrome-like family history. Another patient harbored a germline mutation in the DNA repair factor ATM. In summary, the relapse patients of our cohort were characterized by somatic mutations affecting the RAS pathway, epigenetic and developmental programs and germline mutations in DNA repair pathways.

Genomics and drug profiling of fatal TCF3-HLF-positive acute lymphoblastic leukemia identifies recurrent mutation patterns and therapeutic options.

TCF3-HLF-positive acute lymphoblastic leukemia (ALL) is currently incurable. Using an integrated approach, we uncovered distinct mutation, gene expression and drug response profiles in TCF3-HLF-positive and treatment-responsive TCF3-PBX1-positive ALL. We identified recurrent intragenic deletions of PAX5 or VPREB1 in constellation with the fusion of TCF3 and HLF. Moreover somatic mutations in the non-translocated allele of TCF3 and a reduction of PAX5 gene dosage in TCF3-HLF ALL suggest cooperation within a restricted genetic context. The enrichment for stem cell and myeloid features in the TCF3-HLF signature may reflect reprogramming by TCF3-HLF of a lymphoid-committed cell of origin toward a hybrid, drug-resistant hematopoietic state. Drug response profiling of matched patient-derived xenografts revealed a distinct profile for TCF3-HLF ALL with resistance to conventional chemotherapeutics but sensitivity to glucocorticoids, anthracyclines and agents in clinical development. Striking on-target sensitivity was achieved with the BCL2-specific inhibitor venetoclax (ABT-199). This integrated approach thus provides alternative treatment options for this deadly disease.

Phylogenetic comparative assembly.

BACKGROUND: Recent high throughput sequencing technologies are capable of generating a huge amount of data for bacterial genome sequencing projects. Although current sequence assemblers successfully merge the overlapping reads, often several contigs remain which cannot be assembled any further. It is still costly and time consuming to close all the gaps in order to acquire the whole genomic sequence. RESULTS: Here we propose an algorithm that takes several related genomes and their phylogenetic relationships into account to create a graph that contains the likelihood for each pair of contigs to be adjacent. Subsequently, this graph can be used to compute a layout graph that shows the most promising contig adjacencies in order to aid biologists in finishing the complete genomic sequence. The layout graph shows unique contig orderings where possible, and the best alternatives where necessary. CONCLUSIONS: Our new algorithm for contig ordering uses sequence similarity as well as phylogenetic information to estimate adjacencies of contigs. An evaluation of our implementation shows that it performs better than recent approaches while being much faster at the same time.

ProgeniDB: a novel cell lineage database for generation associated phenotypic behavior in cell-based assays.

ProgeniDB is a web-accessible cell lineage database, which comprises data derived from time-lapse microscopy image sequences, using our novel encoding software tool called ProgeniTRAK. This database provides information on qualitative and quantitative changes in cellular dynamics in response to anti-cancer drugs variation, which is critical for pharmacodynamic modeling and validation. The lineage data for each progenitor cell centers around the construction of a bifurcation map, locating in both space and time critical cellular events such as cell division and cell death. This database provides the opportunity to select a set of lineages based on different experimental constraints and observe the associated phenotypic behavior of progenitor cells. The output generated from this database depicts the cellular dynamic of selected population which can be downloaded to explore and exploit in relation to hypothesis-driven drug discovery. A case study is presented where we interrogate drug-induced perturbations on ensuing cell lineages. The database is available at http://biodiversity.cs.cf.ac.uk/cymart/progenidb.html.