Correction: Defining the microbial transcriptional response to colitis through integrated host and microbiome profiling.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Sequencing of human genomes with nanopore technology.

Whole-genome sequencing (WGS) is becoming widely used in clinical medicine in diagnostic contexts and to inform treatment choice. Here we evaluate the potential of the Oxford Nanopore Technologies (ONT) MinION long-read sequencer for routine WGS by sequencing the reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome and severe immune dysregulation. We develop and apply a novel reference panel-free analytical method to infer and then exploit phase information which improves single-nucleotide variant (SNV) calling performance from otherwise modest levels. In the clinical sample, we identify and directly phase two non-synonymous de novo variants in SAMD9L, (OMIM #159550) inferring that they lie on the same paternal haplotype. Whilst consensus SNV-calling error rates from ONT data remain substantially higher than those from short-read methods, we demonstrate the substantial benefits of analytical innovation. Ongoing improvements to base-calling and SNV-calling methodology must continue for nanopore sequencing to establish itself as a primary method for clinical WGS.

The Light Chain IgLV3-21 Defines a New Poor Prognostic Subgroup in Chronic Lymphocytic Leukemia: Results of a Multicenter Study.

Purpose: Unmutated (UM) immunoglobulin heavy chain variable region (IgHV) status or IgHV3-21 gene usage is associated with poor prognosis in chronic lymphocytic leukemia (CLL) patients. Interestingly, IgHV3-21 is often co-expressed with light chain IgLV3-21, which is potentially able to trigger cell-autonomous BCR-mediated signaling. However, this light chain has never been characterized independently of the heavy chain IgHV3-21.Experimental Design: We performed total RNA sequencing in 32 patients and investigated IgLV3-21 prognostic impact in terms of treatment-free survival (TFS) and overall survival (OS) in 3 other independent cohorts for a total of 813 patients. IgLV3-21 presence was tested by real-time PCR and confirmed by Sanger sequencing.Results: Using total RNA sequencing to characterize 32 patients with high-risk CLL, we found a high frequency (28%) of IgLV3-21 rearrangements. Gene set enrichment analysis revealed that these patients express higher levels of genes responsible for ribosome biogenesis and translation initiation (P < 0.0001) as well as MYC target genes (P = 0.0003). Patients with IgLV3-21 rearrangements displayed a significantly shorter TFS and OS (P < 0.05), particularly those with IgHV mutation. In each of the three independent validation cohorts, we showed that IgLV3-21 rearrangements-similar to UM IgHV status-conferred poor prognosis compared with mutated IgHV (P < 0.0001). Importantly, we confirmed by multivariate analysis that this was independent of IgHV mutational status or subset #2 stereotyped receptor (P < 0.0001).Conclusions: We have demonstrated for the first time that a light chain can affect CLL prognosis and that IgLV3-21 light chain usage defines a new subgroup of CLL patients with poor prognosis. Clin Cancer Res; 24(20); 5048-57. ©2018 AACR.

Genomic Response to Vitamin D Supplementation in the Setting of a Randomized, Placebo-Controlled Trial.

BACKGROUND: Vitamin D deficiency has been associated with multiple diseases, but the causal relevance and underlying processes are not fully understood. Elucidating the mechanisms of action of drug treatments in humans is challenging, but application of functional genomic approaches in randomized trials may afford an opportunity to systematically assess molecular responses. METHODS: In the Biochemical Efficacy and Safety Trial of Vitamin D (BEST-D), a double-blind, placebo-controlled, dose-finding, randomized clinical trial, 305 community-dwelling individuals aged over 65 years were randomly allocated to treatment with vitamin D3 4000 IU, 2000 IU or placebo daily for 12 months. Genome-wide genotypes at baseline, and transcriptome and plasma levels of cytokines (IFN-γ, IL-10, IL-8, IL-6 and TNF-α) at baseline and after 12 months, were measured. The trial had >90% power to detect 1.2-fold changes in gene expression. FINDINGS: Allocation to vitamin D for 12-months was associated with 2-fold higher plasma levels of 25-hydroxy-vitamin D (25[OH]D, 4000 IU regimen), but had no significant effect on whole-blood gene expression (FDR < 5%) or on plasma levels of cytokines compared with placebo. In pre-specified analysis, rs7041 (intron variant, GC) had a significant effect on circulating levels of 25(OH)D in the low dose, but not in the placebo or high dose vitamin D regimen. A gene expression quantitative trait locus analysis (eQTL) demonstrated evidence of 31,568 cis-eQTLs (unique SNP-probe pairs) among individuals at baseline and 34,254 after supplementation for 12 months (any dose). No significant associations involving vitamin D supplementation response eQTLs were found. INTERPRETATION: We performed a comprehensive functional genomics and molecular analysis of vitamin D supplementation in a randomized, placebo-controlled trial. Although this study was limited to mostly Caucasian individuals aged over 65 years, the results differ from many previous studies and do not support a strong effect of vitamin D on long-term transcriptomic changes in blood or on plasma cytokine levels. The trial demonstrates the feasibility of applying functional genomic and genetic approaches in randomized trials to assess molecular and individual level responses. KEY RESULT: Supplementation with high-dose vitamin D in older people for 12 months in a randomized, placebo-controlled trial had no significant effect on gene expression or on plasma concentrations of selected cytokines. TRIAL REGISTRATION: SRCTN registry (Number 07034656) and the European Clinical Trials Database (EudraCT Number 2011-005763-24).

Investigation into the role of the germline epigenome in the transmission of glucocorticoid-programmed effects across generations.

BACKGROUND: Early life exposure to adverse environments affects cardiovascular and metabolic systems in the offspring. These programmed effects are transmissible to a second generation through both male and female lines, suggesting germline transmission. We have previously shown that prenatal overexposure to the synthetic glucocorticoid dexamethasone (Dex) in rats reduces birth weight in the first generation (F1), a phenotype which is transmitted to a second generation (F2), particularly through the male line. We hypothesize that Dex exposure affects developing germ cells, resulting in transmissible alterations in DNA methylation, histone marks and/or small RNA in the male germline. RESULTS: We profile epigenetic marks in sperm from F1 Sprague Dawley rats expressing a germ cell-specific GFP transgene following Dex or vehicle treatment of the mothers, using methylated DNA immunoprecipitation sequencing, small RNA sequencing and chromatin immunoprecipitation sequencing for H3K4me3, H3K4me1, H3K27me3 and H3K9me3. Although effects on birth weight are transmitted to the F2 generation through the male line, no differences in DNA methylation, histone modifications or small RNA were detected between germ cells and sperm from Dex-exposed animals and controls. CONCLUSIONS: Although the phenotype is transmitted to a second generation, we are unable to detect specific changes in DNA methylation, common histone modifications or small RNA profiles in sperm. Dex exposure is associated with more variable 5mC levels, particularly at non-promoter loci. Although this could be one mechanism contributing to the observed phenotype, other germline epigenetic modifications or non-epigenetic mechanisms may be responsible for the transmission of programmed effects across generations in this model.

UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy.

Unique Molecular Identifiers (UMIs) are random oligonucleotide barcodes that are increasingly used in high-throughput sequencing experiments. Through a UMI, identical copies arising from distinct molecules can be distinguished from those arising through PCR amplification of the same molecule. However, bioinformatic methods to leverage the information from UMIs have yet to be formalized. In particular, sequencing errors in the UMI sequence are often ignored or else resolved in an ad hoc manner. We show that errors in the UMI sequence are common and introduce network-based methods to account for these errors when identifying PCR duplicates. Using these methods, we demonstrate improved quantification accuracy both under simulated conditions and real iCLIP and single-cell RNA-seq data sets. Reproducibility between iCLIP replicates and single-cell RNA-seq clustering are both improved using our proposed network-based method, demonstrating the value of properly accounting for errors in UMIs. These methods are implemented in the open source UMI-tools software package.

Transcriptomic profiling of purified patient-derived dopamine neurons identifies convergent perturbations and therapeutics for Parkinson's disease.

While induced pluripotent stem cell (iPSC) technologies enable the study of inaccessible patient cell types, cellular heterogeneity can confound the comparison of gene expression profiles between iPSC-derived cell lines. Here, we purified iPSC-derived human dopaminergic neurons (DaNs) using the intracellular marker, tyrosine hydroxylase. Once purified, the transcriptomic profiles of iPSC-derived DaNs appear remarkably similar to profiles obtained from mature post-mortem DaNs. Comparison of the profiles of purified iPSC-derived DaNs derived from Parkinson's disease (PD) patients carrying LRRK2 G2019S variants to controls identified significant functional convergence amongst differentially-expressed (DE) genes. The PD LRRK2-G2019S associated profile was positively matched with expression changes induced by the Parkinsonian neurotoxin rotenone and opposed by those induced by clioquinol, a compound with demonstrated therapeutic efficacy in multiple PD models. No functional convergence amongst DE genes was observed following a similar comparison using non-purified iPSC-derived DaN-containing populations, with cellular heterogeneity appearing a greater confound than genotypic background.

11ß-hydroxysteroid dehydrogenase-1 deficiency alters the gut microbiome response to Western diet.

The enzyme 11ß-hydroxysteroid dehydrogenase (11ß-HSD) interconverts active glucocorticoids and their intrinsically inert 11-keto forms. The type 1 isozyme, 11ß-HSD1, predominantly reactivates glucocorticoids in vivo and can also metabolise bile acids. 11ß-HSD1-deficient mice show altered inflammatory responses and are protected against the adverse metabolic effects of a high-fat diet. However, the impact of 11ß-HSD1 on the composition of the gut microbiome has not previously been investigated. We used high-throughput 16S rDNA amplicon sequencing to characterise the gut microbiome of 11ß-HSD1-deficient and C57Bl/6 control mice, fed either a standard chow diet or a cholesterol- and fat-enriched 'Western' diet. 11ß-HSD1 deficiency significantly altered the composition of the gut microbiome, and did so in a diet-specific manner. On a Western diet, 11ß-HSD1 deficiency increased the relative abundance of the family Bacteroidaceae, and on a chow diet, it altered relative abundance of the family Prevotellaceae Our results demonstrate that (i) genetic effects on host-microbiome interactions can depend upon diet and (ii) that alterations in the composition of the gut microbiome may contribute to the aspects of the metabolic and/or inflammatory phenotype observed with 11ß-HSD1 deficiency.

Long noncoding RNAs in B-cell development and activation.

Long noncoding RNAs (lncRNAs) are potentially important regulators of cell differentiation and development, but little is known about their roles in B lymphocytes. Using RNA-seq and de novo transcript assembly, we identified 4516 lncRNAs expressed in 11 stages of B-cell development and activation. Most of these lncRNAs have not been previously detected, even in the closely related T-cell lineage. Comparison with lncRNAs previously described in human B cells identified 185 mouse lncRNAs that have human orthologs. Using chromatin immunoprecipitation-seq, we classified 20% of the lncRNAs as either enhancer-associated (eRNA) or promoter-associated RNAs. We identified 126 eRNAs whose expression closely correlated with the nearest coding gene, thereby indicating the likely location of numerous enhancers active in the B-cell lineage. Furthermore, using this catalog of newly discovered lncRNAs, we show that PAX5, a transcription factor required to specify the B-cell lineage, bound to and regulated the expression of 109 lncRNAs in pro-B and mature B cells and 184 lncRNAs in acute lymphoblastic leukemia.

Defining the microbial transcriptional response to colitis through integrated host and microbiome profiling.

The gut microbiome is significantly altered in inflammatory bowel diseases, but the basis of these changes is not well understood. We have combined metagenomic and metatranscriptomic profiling of the gut microbiome to assess modifications to both bacterial community structure and transcriptional activity in a mouse model of colitis. By using transcriptomic analysis of colonic tissue and luminal RNA derived from the host, we have also characterised how host transcription relates to the microbial transcriptional response in inflammation. In colitis, increased abundance and transcription of diverse microbial gene families involved in responses to nutrient deprivation, antimicrobial peptide production and oxidative stress support an adaptation of multiple commensal genera to withstand a diverse set of environmental stressors in the inflammatory environment. These data are supported by a transcriptional signature of activated macrophages and granulocytes in the gut lumen during colitis, a signature that includes the transcription of the key antimicrobial genes S100a8 and S100a9 (calprotectin). Genes involved in microbial resistance to oxidative stress, including Dps/ferritin, Fe-dependent peroxidase and glutathione S-transferase were identified as changing to a greater extent at the level of transcription than would be predicted by DNA abundance changes, implicating a role for increased oxygen tension and/or host-derived reactive oxygen species in driving transcriptional changes in commensal microbes.

CGAT: a model for immersive personalized training in computational genomics.

How should the next generation of genomics scientists be trained while simultaneously pursuing high quality and diverse research? CGAT, the Computational Genomics Analysis and Training programme, was set up in 2010 by the UK Medical Research Council to complement its investment in next-generation sequencing capacity. CGAT was conceived around the twin goals of training future leaders in genome biology and medicine, and providing much needed capacity to UK science for analysing genome scale data sets. Here we outline the training programme employed by CGAT and describe how it dovetails with collaborative research projects to launch scientists on the road towards independent research careers in genomics.

Identification of a candidate prognostic gene signature by transcriptome analysis of matched pre- and post-treatment prostatic biopsies from patients with advanced prostate cancer.

BACKGROUND: Although chemotherapy for prostate cancer (PCa) can improve patient survival, some tumours are chemo-resistant. Tumour molecular profiles may help identify the mechanisms of drug action and identify potential prognostic biomarkers. We performed in vivo transcriptome profiling of pre- and post-treatment prostatic biopsies from patients with advanced hormone-naive prostate cancer treated with docetaxel chemotherapy and androgen deprivation therapy (ADT) with an aim to identify the mechanisms of drug action and identify prognostic biomarkers. METHODS: RNA sequencing (RNA-Seq) was performed on biopsies from four patients before and ~22 weeks after docetaxel and ADT initiation. Gene fusion products and differentially-regulated genes between treatment pairs were identified using TopHat and pathway enrichment analyses undertaken. Publically available datasets were interrogated to perform survival analyses on the gene signatures identified using cBioportal. RESULTS: A number of genomic rearrangements were identified including the TMPRSS2/ERG fusion and 3 novel gene fusions involving the ETS family of transcription factors in patients, both pre and post chemotherapy. In total, gene expression analyses showed differential expression of at least 2 fold in 575 genes in post-chemotherapy biopsies. Of these, pathway analyses identified a panel of 7 genes (ADAM7, FAM72B, BUB1B, CCNB1, CCNB2, TTK, CDK1), including a cell cycle-related geneset, that were differentially-regulated following treatment with docetaxel and ADT. Using cBioportal to interrogate the MSKCC-Prostate Oncogenome Project dataset we observed a statistically-significant reduction in disease-free survival of patients with tumours exhibiting alterations in gene expression of the above panel of 7 genes (p = 0.015). CONCLUSIONS: Here we report on the first "real-time" in vivo RNA-Seq-based transcriptome analysis of clinical PCa from pre- and post-treatment TRUSS-guided biopsies of patients treated with docetaxel chemotherapy plus ADT. We identify a chemotherapy-driven PCa transcriptome profile which includes the down-regulation of important positive regulators of cell cycle progression. A 7 gene signature biomarker panel has also been identified in high-risk prostate cancer patients to be of prognostic value. Future prospective study is warranted to evaluate the clinical value of this panel.

Population and single-cell genomics reveal the Aire dependency, relief from Polycomb silencing, and distribution of self-antigen expression in thymic epithelia.

Promiscuous gene expression (PGE) by thymic epithelial cells (TEC) is essential for generating a diverse T cell antigen receptor repertoire tolerant to self-antigens, and thus for avoiding autoimmunity. Nevertheless, the extent and nature of this unusual expression program within TEC populations and single cells are unknown. Using deep transcriptome sequencing of carefully identified mouse TEC subpopulations, we discovered a program of PGE that is common between medullary (m) and cortical TEC, further elaborated in mTEC, and completed in mature mTEC expressing the autoimmune regulator gene (Aire). TEC populations are capable of expressing up to 19,293 protein-coding genes, the highest number of genes known to be expressed in any cell type. Remarkably, in mouse mTEC, Aire expression alone positively regulates 3980 tissue-restricted genes. Notably, the tissue specificities of these genes include known targets of autoimmunity in human AIRE deficiency. Led by the observation that genes induced by Aire expression are generally characterized by a repressive chromatin state in somatic tissues, we found these genes to be strongly associated with H3K27me3 marks in mTEC. Our findings are consistent with AIRE targeting and inducing the promiscuous expression of genes previously epigenetically silenced by Polycomb group proteins. Comparison of the transcriptomes of 174 single mTEC indicates that genes induced by Aire expression are transcribed stochastically at low cell frequency. Furthermore, when present, Aire expression-dependent transcript levels were 16-fold higher, on average, in individual TEC than in the mTEC population.

IRF5:RelA interaction targets inflammatory genes in macrophages.

Interferon Regulatory Factor 5 (IRF5) plays a major role in setting up an inflammatory macrophage phenotype, but the molecular basis of its transcriptional activity is not fully understood. In this study, we conduct a comprehensive genome-wide analysis of IRF5 recruitment in macrophages stimulated with bacterial lipopolysaccharide and discover that IRF5 binds to regulatory elements of highly transcribed genes. Analysis of protein:DNA microarrays demonstrates that IRF5 recognizes the canonical IRF-binding (interferon-stimulated response element [ISRE]) motif in vitro. However, IRF5 binding in vivo appears to rely on its interactions with other proteins. IRF5 binds to a noncanonical composite PU.1:ISRE motif, and its recruitment is aided by RelA. Global gene expression analysis in macrophages deficient in IRF5 and RelA highlights the direct role of the RelA:IRF5 cistrome in regulation of a subset of key inflammatory genes. We map the RelA:IRF5 interaction domain and suggest that interfering with it would offer selective targeting of macrophage inflammatory activities.

Long non-coding RNAs and enhancer RNAs regulate the lipopolysaccharide-induced inflammatory response in human monocytes.

Early reports indicate that long non-coding RNAs (lncRNAs) are novel regulators of biological responses. However, their role in the human innate immune response, which provides the initial defence against infection, is largely unexplored. To address this issue, here we characterize the long non-coding RNA transcriptome in primary human monocytes using RNA sequencing. We identify 76 enhancer RNAs (eRNAs), 40 canonical lncRNAs, 65 antisense lncRNAs and 35 regions of bidirectional transcription (RBT) that are differentially expressed in response to bacterial lipopolysaccharide (LPS). Crucially, we demonstrate that knockdown of nuclear-localized, NF-κB-regulated, eRNAs (IL1ß-eRNA) and RBT (IL1ß-RBT46) surrounding the IL1ß locus, attenuates LPS-induced messenger RNA transcription and release of the proinflammatory mediators, IL1ß and CXCL8. We predict that lncRNAs can be important regulators of the human innate immune response.

Targeting polycomb to pericentric heterochromatin in embryonic stem cells reveals a role for H2AK119u1 in PRC2 recruitment.

The mechanisms by which the major Polycomb group (PcG) complexes PRC1 and PRC2 are recruited to target sites in vertebrate cells are not well understood. Building on recent studies that determined a reciprocal relationship between DNA methylation and Polycomb activity, we demonstrate that, in methylation-deficient embryonic stem cells (ESCs), CpG density combined with antagonistic effects of H3K9me3 and H3K36me3 redirects PcG complexes to pericentric heterochromatin and gene-rich domains. Surprisingly, we find that PRC1-linked H2A monoubiquitylation is sufficient to recruit PRC2 to chromatin in vivo, suggesting a mechanism through which recognition of unmethylated CpG determines the localization of both PRC1 and PRC2 at canonical and atypical target sites. We discuss our data in light of emerging evidence suggesting that PcG recruitment is a default state at licensed chromatin sites, mediated by interplay between CpG hypomethylation and counteracting H3 tail modifications.

CGAT: computational genomics analysis toolkit.

Computational genomics seeks to draw biological inferences from genomic datasets, often by integrating and contextualizing next-generation sequencing data. CGAT provides an extensive suite of tools designed to assist in the analysis of genome scale data from a range of standard file formats. The toolkit enables filtering, comparison, conversion, summarization and annotation of genomic intervals, gene sets and sequences. The tools can both be run from the Unix command line and installed into visual workflow builders, such as Galaxy.