Sagar Chhangawala

Multi-platform assessment of transcriptome profiling using RNA-seq in the ABRF next-generation sequencing study.

High-throughput RNA sequencing (RNA-seq) greatly expands the potential for genomics discoveries, but the wide variety of platforms, protocols and performance capabilitites has created the need for comprehensive reference data. Here we describe the Association of Biomolecular Resource Facilities next-generation sequencing (ABRF-NGS) study on RNA-seq. We carried out replicate experiments across 15 laboratory sites using reference RNA standards to test four protocols (poly-A–selected, ribo-depleted, size-selected and degraded) on five sequencing platforms (Illumina HiSeq, Life Technologies PGM and Proton, Pacific Biosciences RS and Roche 454). The results show high intraplatform (Spearman rank R > 0.86) and inter-platform (R > 0.83) concordance for expression measures across the deep-count platforms, but highly variable efficiency and cost for splice junction and variant detection between all platforms. For intact RNA, gene expression profiles from rRNA-depletion and poly-A enrichment are similar. In addition, rRNA depletion enables effective analysis of degraded RNA samples. This study provides a broad foundation for cross-platform standardization, evaluation and improvement of RNA-seq.

Aerobic glycolysis promotes T helper 1 cell differentiation through an epigenetic mechanism

Metabolic support for T cell functions For immunological T cells, responding to infections is energetically demanding. T cells rewire their metabolism so that they rely more heavily on aerobic glycolysis. This helps them to support important effector functions such as secreting the cytokine interferon γ (IFNγ). Peng et al. now provide insight into how aerobic glycolysis promotes T cell effector function. Activated T cells express the aerobic glycolysis–supporting enzyme lactate dehydrogenase A (LDHA), allowing these cells to maintain high amounts of acetyl–coenzyme A, which in turn promotes histone acetylation and transcription of cytokines such as IFNγ. Engineered mice whose T cells lacked LDHA were protected from IFNγ-dependent pathologies that often characterize autoinflammatory diseases. Science, this issue p. 481 Immunological T cell effector functions are promoted by metabolic activity–modifying levels of histone acetylation. Aerobic glycolysis (the Warburg effect) is a metabolic hallmark of activated T cells and has been implicated in augmenting effector T cell responses, including expression of the proinflammatory cytokine interferon-γ (IFN-γ), via 3′ untranslated region (3′UTR)–mediated mechanisms. Here, we show that lactate dehydrogenase A (LDHA) is induced in activated T cells to support aerobic glycolysis but promotes IFN-γ expression independently of its 3′UTR. Instead, LDHA maintains high concentrations of acetyl–coenzyme A to enhance histone acetylation and transcription of Ifng. Ablation of LDHA in T cells protects mice from immunopathology triggered by excessive IFN-γ expression or deficiency of regulatory T cells. These findings reveal an epigenetic mechanism by which aerobic glycolysis promotes effector T cell differentiation and suggest that LDHA may be targeted therapeutically in autoinflammatory diseases.

Cancer Immunosurveillance by Tissue-Resident Innate Lymphoid Cells and Innate-like T Cells

Malignancy can be suppressed by the immune system in a process termed immunosurveillance. However, to what extent immunosurveillance occurs in spontaneous cancers and the composition of participating cell types remains obscure. Here, we show that cell transformation triggers a tissue-resident lymphocyte response in oncogene-induced murine cancer models. Non-circulating cytotoxic lymphocytes, derived from innate, T cell receptor (TCR)αβ, and TCRγδ lineages, expand in early tumors. Characterized by high expression of NK1.1, CD49a, and CD103, these cells share a gene-expression signature distinct from those of conventional NK cells, T cells, and invariant NKT cells. Generation of these lymphocytes is dependent on the cytokine IL-15, but not the transcription factor Nfil3 that is required for the differentiation of tumor-infiltrating NK cells, and IL-15 deficiency, but not Nfil3 deficiency, results in accelerated tumor growth. These findings reveal a tumor-elicited immunosurveillance mechanism that engages unconventional type-1-like innate lymphoid cells and type 1 innate-like T cells.

Geospatial Resolution of Human and Bacterial Diversity with City-Scale Metagenomics

SUMMARY The panoply of microorganisms and other species present in our environment influence human health and disease, especially in cities, but have not been profiled with metagenomics at a city-wide scale. We sequenced DNA from surfaces across the entire New York City (NYC) subway system, the Gowanus Canal, and public parks. Nearly half of the DNA (48%) does not match any known organism; identified organisms spanned 1,688 bacterial, viral, archaeal, and eukaryotic taxa, which were enriched for harmless genera associated with skin (e.g., Acinetobacter). Predicted ancestry of human DNA left on subway surfaces can recapitulate U.S. Census demographic data, and bacterial signatures can reveal a station’s history, such as marine-associated bacteria in a hurricane-flooded station. Some evidence of pathogens was found (Bacillus anthracis), but a lack of reported cases in NYC suggests that the pathogens represent a normal, urban microbiome. This baseline metagenomic map of NYC could help long-term disease surveillance, bioterrorism threat mitigation, and health management in the built environment of cities.

SPOP mutation leads to genomic instability in prostate cancer

Genomic instability is a fundamental feature of human cancer often resulting from impaired genome maintenance. In prostate cancer, structural genomic rearrangements are a common mechanism driving tumorigenesis. However, somatic alterations predisposing to chromosomal rearrangements in prostate cancer remain largely undefined. Here, we show that SPOP, the most commonly mutated gene in primary prostate cancer modulates DNA double strand break (DSB) repair, and that SPOP mutation is associated with genomic instability. In vivo, SPOP mutation results in a transcriptional response consistent with BRCA1 inactivation resulting in impaired homology-directed repair (HDR) of DSB. Furthermore, we found that SPOP mutation sensitizes to DNA damaging therapeutic agents such as PARP inhibitors. These results implicate SPOP as a novel participant in DSB repair, suggest that SPOP mutation drives prostate tumorigenesis in part through genomic instability, and indicate that mutant SPOP may increase response to DNA-damaging therapeutics. DOI: http://dx.doi.org/10.7554/eLife.09207.001

GuideScan software for improved single and paired CRISPR guide RNA design

We present GuideScan software for the design of CRISPR guide RNA libraries that can be used to edit coding and noncoding genomic regions. GuideScan produces high-density sets of guide RNAs (gRNAs) for single- and paired-gRNA genome-wide screens. We also show that the trie data structure of GuideScan enables the design of gRNAs that are more specific than those designed by existing tools.

The impact of read length on quantification of differentially expressed genes and splice junction detection

BackgroundThe initial next-generation sequencing technologies produced reads of 25 or 36 bp, and only from a single-end of the library sequence. Currently, it is possible to reliably produce 300 bp paired-end sequences for RNA expression analysis. While read lengths have consistently increased, people have assumed that longer reads are more informative and that paired-end reads produce better results than single-end reads. We used paired-end 101 bp reads and trimmed them to simulate different read lengths, and also separated the pairs to produce single-end reads. For each read length and paired status, we evaluated differential expression levels between two standard samples and compared the results to those obtained by qPCR.ResultsWe found that, with the exception of 25 bp reads, there is little difference for the detection of differential expression regardless of the read length. Once single-end reads are at a length of 50 bp, the results do not change substantially for any level up to, and including, 100 bp paired-end. However, splice junction detection significantly improves as the read length increases with 100 bp paired-end showing the best performance. We performed the same analysis on two ENCODE samples and found consistent results confirming that our conclusions have broad application.ConclusionsA researcher could save substantial resources by using 50 bp single-end reads for differential expression analysis instead of using longer reads. However, splicing detection is unquestionably improved by paired-end and longer reads. Therefore, an appropriate read length should be used based on the final goal of the study.

Oncogenic BRAF disrupts thyroid morphogenesis and function via twist expression

Thyroid cancer is common, yet the sequence of alterations that promote tumor formation are incompletely understood. Here, we describe a novel model of thyroid carcinoma in zebrafish that reveals temporal changes due to BRAFV600E. Through the use of real-time in vivo imaging, we observe disruption in thyroid follicle structure that occurs early in thyroid development. Combinatorial treatment using BRAF and MEK inhibitors reversed the developmental effects induced by BRAFV600E. Adult zebrafish expressing BRAFV600E in thyrocytes developed invasive carcinoma. We identified a gene expression signature from zebrafish thyroid cancer that is predictive of disease-free survival in patients with papillary thyroid cancer. Gene expression studies nominated TWIST2 as a key effector downstream of BRAF. Using CRISPR/Cas9 to genetically inactivate a TWIST2 orthologue, we suppressed the effects of BRAFV600E and restored thyroid morphology and hormone synthesis. These data suggest that expression of TWIST2 plays a role in an early step of BRAFV600E-mediated transformation. DOI: http://dx.doi.org/10.7554/eLife.20728.001

RSK1 activation promotes invasion in nodular melanoma

Materials and methods Expression of 141 signaling proteins was evaluated by protein pathway array in 3 Radial Growth Phase (RGP)/ SSM and 3 Vertical Growth Phase (VGP)/NM cell lines. The impact of p90-ribosomal-S6-kinase (RSK1) and its inhibition on proliferation, migration and invasion was assessed in SSM and NM cell lines, and confirmed using NM cells treated with a RSK inhibitor (BI-D1870) in microarray profiling studies. The effect of constitutive RSK1 activation in vivo was further studied using a zebrafish model.

Selenoprotein H is an essential regulator of redox homeostasis that cooperates with p53 in development and tumorigenesis

Significance Dietary selenium and selenoproteins play important roles in regulating redox processes that impact human health. The human genome includes 25 genes for selenoproteins, which have diverse roles in redox homeostasis, thyroid hormone metabolism, endoplasmic reticulum quality control, selenium transport, and other functions. Selenoprotein H (seph) is a recently identified nucleolar oxidoreductase with DNA-binding properties whose function is not well understood. In this work, we used a unique combination of unbiased metabolomic and transcriptomic approaches in zebrafish to discover that seph is an essential regulator of redox homeostasis that regulates p53. In addition, we demonstrate the seph-deficient adults are prone to chemically induced carcinogenesis. Our results suggest that seph suppresses oxidative stress and DNA damage in the nucleolus. Selenium, an essential micronutrient known for its cancer prevention properties, is incorporated into a class of selenocysteine-containing proteins (selenoproteins). Selenoprotein H (SepH) is a recently identified nucleolar oxidoreductase whose function is not well understood. Here we report that seph is an essential gene regulating organ development in zebrafish. Metabolite profiling by targeted LC-MS/MS demonstrated that SepH deficiency impairs redox balance by reducing the levels of ascorbate and methionine, while increasing methionine sulfoxide. Transcriptome analysis revealed that SepH deficiency induces an inflammatory response and activates the p53 pathway. Consequently, loss of seph renders larvae susceptible to oxidative stress and DNA damage. Finally, we demonstrate that seph interacts with p53 deficiency in adulthood to accelerate gastrointestinal tumor development. Overall, our findings establish that seph regulates redox homeostasis and suppresses DNA damage. We hypothesize that SepH deficiency may contribute to the increased cancer risk observed in cohorts with low selenium levels.