Intestinal B cells license metabolic T-cell activation in NASH microbiota/antigen-independently and contribute to fibrosis by IgA-FcR signalling

BACKGROUND & AIMS

The progression of non-alcoholic steatohepatitis (NASH) to fibrosis and hepatocellular carcinoma (HCC) is aggravated by auto-aggressive T cells. The gut-liver axis contributes to NASH, but the mechanisms involved and the consequences for NASH-induced fibrosis and liver cancer remain unknown. We investigated the role of gastrointestinal B cells in the development of NASH, fibrosis and NASH-induced HCC.

METHODS

C57BL/6J wild-type (WT), B cell-deficient and different immunoglobulin-deficient or transgenic mice were fed distinct NASH-inducing diets or standard chow for 6 or 12 months, whereafter NASH, fibrosis, and NASH-induced HCC were assessed and analysed. Specific pathogen-free/germ-free WT and μMT mice (containing B cells only in the gastrointestinal tract) were fed a choline-deficient high-fat diet, and treated with an anti-CD20 antibody, whereafter NASH and fibrosis were assessed. Tissue biopsy samples from patients with simple steatosis, NASH and cirrhosis were analysed to correlate the secretion of immunoglobulins to clinicopathological features. Flow cytometry, immunohistochemistry and single-cell RNA-sequencing analysis were performed in liver and gastrointestinal tissue to characterise immune cells in mice and humans.

RESULTS

Activated intestinal B cells were increased in mouse and human NASH samples and licensed metabolic T-cell activation to induce NASH independently of antigen specificity and gut microbiota. Genetic or therapeutic depletion of systemic or gastrointestinal B cells prevented or reverted NASH and liver fibrosis. IgA secretion was necessary for fibrosis induction by activating CD11b+CCR2+F4/80+CD11c-FCGR1+ hepatic myeloid cells through an IgA-FcR signalling axis. Similarly, patients with NASH had increased numbers of activated intestinal B cells; additionally, we observed a positive correlation between IgA levels and activated FcRg+ hepatic myeloid cells, as well the extent of liver fibrosis.

CONCLUSIONS

Intestinal B cells and the IgA-FcR signalling axis represent potential therapeutic targets for the treatment of NASH.

IMPACT AND IMPLICATIONS

There is currently no effective treatment for non-alcoholic steatohepatitis (NASH), which is associated with a substantial healthcare burden and is a growing risk factor for hepatocellular carcinoma (HCC). We have previously shown that NASH is an auto-aggressive condition aggravated, amongst others, by T cells. Therefore, we hypothesized that B cells might have a role in disease induction and progression. Our present work highlights that B cells have a dual role in NASH pathogenesis, being implicated in the activation of auto-aggressive T cells and the development of fibrosis via activation of monocyte-derived macrophages by secreted immunoglobulins (e.g., IgA). Furthermore, we show that the absence of B cells prevented HCC development. B cell-intrinsic signalling pathways, secreted immunoglobulins, and interactions of B cells with other immune cells are potential targets for combinatorial NASH therapies against inflammation and fibrosis.

A proteo-transcriptomic map of non-alcoholic fatty liver disease signatures

Non-alcoholic fatty liver disease (NAFLD) is a common, progressive liver disease strongly associated with the metabolic syndrome. It is unclear how progression of NAFLD towards cirrhosis translates into systematic changes in circulating proteins. Here, we provide a detailed proteo-transcriptomic map of steatohepatitis and fibrosis during progressive NAFLD. In this multicentre proteomic study, we characterize 4,730 circulating proteins in 306 patients with histologically characterized NAFLD and integrate this with transcriptomic analysis in paired liver tissue. We identify circulating proteomic signatures for active steatohepatitis and advanced fibrosis, and correlate these with hepatic transcriptomics to develop a proteo-transcriptomic signature of 31 markers. Deconvolution of this signature by single-cell RNA sequencing reveals the hepatic cell types likely to contribute to proteomic changes with disease progression. As an exemplar of use as a non-invasive diagnostic, logistic regression establishes a composite model comprising four proteins (ADAMTSL2, AKR1B10, CFHR4 and TREM2), body mass index and type 2 diabetes mellitus status, to identify at-risk steatohepatitis.

Quantitative modeling of human liver reveals dysregulation of glycosphingolipid pathways in nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is an increasingly prevalent disease that is associated with multiple metabolic disturbances, yet the metabolic pathways underlying its progression are poorly understood. Here, we studied metabolic pathways of the human liver across the full histological spectrum of NAFLD. We analyzed whole liver tissue transcriptomics and serum metabolomics data obtained from a large, prospectively enrolled cohort of 206 histologically characterized patients derived from the European NAFLD Registry and developed genome-scale metabolic models (GEMs) of human hepatocytes at different stages of NAFLD. We identified several metabolic signatures in the liver and blood of these patients, specifically highlighting the alteration of vitamins (A, E) and glycosphingolipids, and their link with complex glycosaminoglycans in advanced fibrosis. Furthermore, we derived GEMs and identified metabolic signatures of three common NAFLD-associated gene variants (PNPLA3, TM6SF2, and HSD17B13). The study demonstrates dysregulated liver metabolic pathways which may contribute to the progression of NAFLD.

A de novo paradigm for male infertility

De novo mutations are known to play a prominent role in sporadic disorders with reduced fitness. We hypothesize that de novo mutations play an important role in severe male infertility and explain a portion of the genetic causes of this understudied disorder. To test this hypothesis, we utilize trio-based exome sequencing in a cohort of 185 infertile males and their unaffected parents. Following a systematic analysis, 29 of 145 rare (MAF < 0.1%) protein-altering de novo mutations are classified as possibly causative of the male infertility phenotype. We observed a significant enrichment of loss-of-function de novo mutations in loss-of-function-intolerant genes (p-value = 1.00 × 10^-5) in infertile men compared to controls. Additionally, we detected a significant increase in predicted pathogenic de novo missense mutations affecting missense-intolerant genes (p-value = 5.01 × 10^-4) in contrast to predicted benign de novo mutations. One gene we identify, RBM5, is an essential regulator of male germ cell pre-mRNA splicing and has been previously implicated in male infertility in mice. In a follow-up study, 6 rare pathogenic missense mutations affecting this gene are observed in a cohort of 2,506 infertile patients, whilst we find no such mutations in a cohort of 5,784 fertile men (p-value = 0.03). Our results provide evidence for the role of de novo mutations in severe male infertility and point to new candidate genes affecting fertility.

Transcriptomics Identify Thrombospondin-2 as a Biomarker for NASH and Advanced Liver Fibrosis

BACKGROUND AND AIMS

NAFLD is the most common liver disease worldwide. NASH, the progressive form of NAFLD, and advanced fibrosis are associated with poor outcomes. We searched for their noninvasive biomarkers.

APPROACH AND RESULTS

Global RNA sequencing of liver tissue from 98 patients with biopsy-proven NAFLD was performed. Unsupervised hierarchical clustering well distinguished NASH from nonalcoholic fatty liver (NAFL), and patients with NASH exhibited molecular abnormalities reflecting their pathological features. Transcriptomic analysis identified proteins up-regulated in NASH and/or advanced fibrosis (stage F3-F4), including matricellular glycoprotein thrombospondin-2 (TSP-2), encoded by the thrombospondin 2 (THBS2) gene. The intrahepatic THBS2 expression level showed the highest areas under the receiver operating characteristic curves (AUROCs) of 0.915 and 0.957 for diagnosing NASH and advanced fibrosis, respectively. THBS2 positively correlated with inflammation and ballooning according to NAFLD activity score, serum aspartate aminotransferase and hyaluronic acid (HA) levels, and NAFLD Fibrosis Score (NFS). THBS2 was associated with extracellular matrix and collagen biosynthesis, platelet activation, caspase-mediated cleavage of cytoskeletal proteins, and immune cell infiltration. Serum TSP-2 expression was measured in 213 patients with biopsy-proven NAFLD, was significantly higher in NASH than in NAFL, and increased parallel to fibrosis stage. The AUROCs for predicting NASH and advanced fibrosis were 0.776 and 0.856, respectively, which were comparable to Fibrosis-4 index, serum HA level, and NFS in advanced fibrosis diagnosis. Serum TSP-2 level and platelet count were independent predictors of NASH and advanced fibrosis. Serum TSP-2 levels could stratify patients with NAFLD according to the risk of hepatic complications, including liver cancer and decompensated cirrhotic events.

CONCLUSIONS

TSP-2 may be a useful biomarker for NASH and advanced fibrosis diagnosis in patients with NAFLD.

NASH limits anti-tumour surveillance in immunotherapy-treated HCC

Hepatocellular carcinoma (HCC) can have viral or non-viral causes^1-5. Non-alcoholic steatohepatitis (NASH) is an important driver of HCC. Immunotherapy has been approved for treating HCC, but biomarker-based stratification of patients for optimal response to therapy is an unmet need^6,7. Here we report the progressive accumulation of exhausted, unconventionally activated CD8⁺PD1⁺ T cells in NASH-affected livers. In preclinical models of NASH-induced HCC, therapeutic immunotherapy targeted at programmed death-1 (PD1) expanded activated CD8⁺PD1⁺ T cells within tumours but did not lead to tumour regression, which indicates that tumour immune surveillance was impaired. When given prophylactically, anti-PD1 treatment led to an increase in the incidence of NASH-HCC and in the number and size of tumour nodules, which correlated with increased hepatic CD8⁺PD1⁺CXCR6⁺, TOX⁺, and TNF⁺ T cells. The increase in HCC triggered by anti-PD1 treatment was prevented by depletion of CD8⁺ T cells or TNF neutralization, suggesting that CD8⁺ T cells help to induce NASH-HCC, rather than invigorating or executing immune surveillance. We found similar phenotypic and functional profiles in hepatic CD8⁺PD1⁺ T cells from humans with NAFLD or NASH. A meta-analysis of three randomized phase III clinical trials that tested inhibitors of PDL1 (programmed death-ligand 1) or PD1 in more than 1,600 patients with advanced HCC revealed that immune therapy did not improve survival in patients with non-viral HCC. In two additional cohorts, patients with NASH-driven HCC who received anti-PD1 or anti-PDL1 treatment showed reduced overall survival compared to patients with other aetiologies. Collectively, these data show that non-viral HCC, and particularly NASH-HCC, might be less responsive to immunotherapy, probably owing to NASH-related aberrant T cell activation causing tissue damage that leads to impaired immune surveillance. Our data provide a rationale for stratification of patients with HCC according to underlying aetiology in studies of immunotherapy as a primary or adjuvant treatment.

Transcriptomic profiling across the nonalcoholic fatty liver disease spectrum reveals gene signatures for steatohepatitis and fibrosis

The mechanisms that drive nonalcoholic fatty liver disease (NAFLD) remain incompletely understood. This large multicenter study characterized the transcriptional changes that occur in liver tissue across the NAFLD spectrum as disease progresses to cirrhosis to identify potential circulating markers. We performed high-throughput RNA sequencing on a discovery cohort comprising histologically characterized NAFLD samples from 206 patients. Unsupervised clustering stratified NAFLD on the basis of disease activity and fibrosis stage with differences in age, aspartate aminotransferase (AST), type 2 diabetes mellitus, and carriage of PNPLA3 rs738409, a genetic variant associated with NAFLD. Relative to early disease, we consistently identified 25 differentially expressed genes as fibrosing steatohepatitis progressed through stages F2 to F4. This 25-gene signature was independently validated by logistic modeling in a separate replication cohort (n = 175), and an integrative analysis with publicly available single-cell RNA sequencing data elucidated the likely relative contribution of specific intrahepatic cell populations. Translating these findings to the protein level, SomaScan analysis in more than 300 NAFLD serum samples confirmed that circulating concentrations of proteins AKR1B10 and GDF15 were strongly associated with disease activity and fibrosis stage. Supporting the biological plausibility of these data, in vitro functional studies determined that endoplasmic reticulum stress up-regulated expression of AKR1B10, GDF15, and PDGFA, whereas GDF15 supplementation tempered the inflammatory response in macrophages upon lipid loading and lipopolysaccharide stimulation. This study provides insights into the pathophysiology of progressive fibrosing steatohepatitis, and proof of principle that transcriptomic changes represent potentially tractable and clinically relevant markers of disease progression.

Correlation of Infinium HumanMethylation450K and MethylationEPIC BeadChip arrays in cartilage

DNA methylation of CpG sites is commonly measured using Illumina Infinium BeadChip platforms. The Infinium MethylationEPIC array has replaced the Infinium Methylation450K array. The two arrays use the same technology, with the EPIC array assaying almost double the number of sites than the 450K array. In this study, we compare DNA methylation values of shared CpGs of the same human cartilage samples assayed using both platforms. DNA methylation was measured in 21 human cartilage samples using the both 450K and EPIC arrays. Additional matched 450K and EPIC data in whole tumour and whole blood were downloaded from GEO GSE92580 and GSE86833, respectively. Data were processed using the Bioconductor package Minfi. DNA methylation of six CpG sites was validated for the same 21 cartilage samples by pyrosequencing. In cartilage samples, overall sample correlations of methylation values between arrays were high (Pearson’s r > 0.96). However, 50.5% of CpG sites showed poor correlation (r < 0.2) between arrays. Sites with limited variance and with either very high or very low methylation levels in cartilage exhibited lower correlation values, corroborating prior studies in whole blood. Bisulphite pyrosequencing did not highlight one array as generating more accurate methylation values. For a specific CpG site, the array methylation correlation coefficient differed between cartilage, tumour, and whole blood, reflecting the difference in methylation variance between cell types. Researchers should be cautious when analysing methylation of CpG sites that show low methylation variance within the cell type of interest, regardless of the method used to assay methylation.

Epigenetic modifiers DNMT3A and BCOR are recurrently mutated in CYLD cutaneous syndrome

Patients with CYLD cutaneous syndrome (CCS; syn. Brooke-Spiegler syndrome) carry germline mutations in the tumor suppressor CYLD and develop multiple skin tumors with diverse histophenotypes. Here, we comprehensively profile the genomic landscape of 42 benign and malignant tumors across 13 individuals from four multigenerational families and discover recurrent mutations in epigenetic modifiers DNMT3A and BCOR in 29% of benign tumors. Multi-level and microdissected sampling strikingly reveal that many clones with different DNMT3A mutations exist in these benign tumors, suggesting that intra-tumor heterogeneity is common. Integrated genomic, methylation and transcriptomic profiling in selected tumors suggest that isoform-specific DNMT3A2 mutations are associated with dysregulated methylation. Phylogenetic and mutational signature analyses confirm cylindroma pulmonary metastases from primary skin tumors. These findings contribute to existing paradigms of cutaneous tumorigenesis and metastasis.

High-throughput sequencing identifies aetiology-dependent differences in ductular reaction in human chronic liver disease

Ductular reaction (DR) represents the activation of hepatic progenitor cells (HPCs) and has been associated with features of advanced chronic liver disease; yet it is not clear whether these cells contribute to disease progression and how the composition of their micro-environment differs depending on the aetiology. This study aimed to identify HPC-associated signalling pathways relevant in different chronic liver diseases using a high-throughput sequencing approach. DR/HPCs were isolated using laser microdissection from patient samples diagnosed with HCV or primary sclerosing cholangitis (PSC), as models for hepatocellular or biliary regeneration. Key signals were validated at the protein level for a cohort of 56 patients (20 early and 36 advanced stage). In total, 330 genes were significantly differentially expressed between the HPCs in HCV and PSC. Recruitment and homing of inflammatory cells were distinctly different depending on the aetiology. HPCs in PSC were characterised by a response to oxidative stress (e.g. JUN, VNN1) and neutrophil-attractant chemokines (CXCL5, CXCL6, IL-8), whereas HPCs in HCV were identified by T- and B-lymphocyte infiltration. Moreover, we found that communication between HPCs and macrophages was aetiology driven. In PSC, a high frequency of CCL28-positive macrophages was observed in the portal infiltrate, already in early disease in the absence of advanced fibrosis, while in HCV, HPCs showed a strong expression of the macrophage scavenger receptor MARCO. Interestingly, DR/HPCs in PSC showed more deposition of ECM (e.g. FN1, LAMC2, collagens) compared to HCV, where an increase of pro-invasive genes (e.g. PDGFRA, IGF2) was observed. Additionally, endothelial cells in the vicinity of DR/HPCs showed differential immunopositivity (e.g. IGF2 and INHBA expression). In conclusion, our data shine light on the role of DR/HPCs in immune signalling, fibrogenesis and angiogenesis in chronic liver disease. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

An ancient germ cell-specific RNA-binding protein protects the germline from cryptic splice site poisoning

Male germ cells of all placental mammals express an ancient nuclear RNA binding protein of unknown function called RBMXL2. Here we find that deletion of the retrogene encoding RBMXL2 blocks spermatogenesis. Transcriptome analyses of age-matched deletion mice show that RBMXL2 controls splicing patterns during meiosis. In particular, RBMXL2 represses the selection of aberrant splice sites and the insertion of cryptic and premature terminal exons. Our data suggest a Rbmxl2 retrogene has been conserved across mammals as part of a splicing control mechanism that is fundamentally important to germ cell biology. We propose that this mechanism is essential to meiosis because it buffers the high ambient concentrations of splicing activators, thereby preventing poisoning of key transcripts and disruption to gene expression by aberrant splice site selection.

In humans and other mammals, a sperm from a male fuses with an egg cell from a female to produce an embryo that may ultimately grow into a new individual. Sperm and egg cells are made when certain cells in the body divide in a process called meiosis. Many proteins are required for meiosis to happen and these proteins are made using instructions provided by genes, which are made of a molecule called DNA.

The DNA within a gene is transcribed to make molecules of ribonucleic acid (or RNA for short). The cell then modifies many of these RNAs in a process called splicing before using them as templates to make proteins. During splicing, segments of RNA known as introns are discarded and other segments termed exons are joined together. Some exons may also be removed from RNAs in different combinations to create different proteins from the same gene.

A protein called RBMXL2 is able to bind to RNA molecules and is only made during and after meiosis in humans and most other mammals. RBMXL2 can also bind to other proteins that are known to be involved in controlling splicing of RNAs, but its role in splicing remains unclear.

To address this question, Ehrmann et al. studied the gene that encodes the RBMXL2 protein in mice. Removing this gene prevented male mice from being able to make sperm. Further experiments using a technique called RNA sequencing showed that the RBMXL2 protein helps to ensure that splicing happens correctly by preventing bits of exons and introns in mouse genes from being rearranged. These findings suggest that the gene encoding RBMXL2 is part of a splicing control mechanism that is important for making sperm and egg cells.

The work of Ehrmann et al. could eventually help some couples understand why they have problems conceiving children. Male infertility is poorly understood, and not knowing its causes can harm the mental health of affected men. Furthermore, these findings may help researchers to understand the role of a closely related protein called RBMY that has also been linked to infertility in men, but is much more difficult to study.

Deconstructing Retinal Organoids: Single Cell RNA-Seq Reveals the Cellular Components of Human Pluripotent Stem Cell-Derived Retina

The rapid improvements in single cell sequencing technologies and analyses afford greater scope for dissecting organoid cultures composed of multiple cell types and create an opportunity to interrogate these models to understand tissue biology, cellular behavior and interactions. To this end, retinal organoids generated from human embryonic stem cells (hESCs) were analyzed by single cell RNA-sequencing (scRNA-Seq) at three time points of differentiation. Combinatorial data from all time points revealed the presence of nine clusters, five of which corresponded to key retinal cell types: retinal pigment epithelium (RPE), retinal ganglion cells (RGCs), cone and rod photoreceptors, and Müller glia. The remaining four clusters expressed genes typical of mitotic cells, extracellular matrix components and those involved in homeostasis. The cell clustering analysis revealed the decreasing presence of mitotic cells and RGCs, formation of a distinct RPE cluster, the emergence of cone and rod photoreceptors from photoreceptor precursors, and an increasing number of Müller glia cells over time. Pseudo-time analysis resembled the order of cell birth during retinal development, with the mitotic cluster commencing the trajectory and the large majority of Müller glia completing the time line. Together, these data demonstrate the feasibility and potential of scRNA-Seq to dissect the inherent complexity of retinal organoids and the orderly birth of key retinal cell types. Stem Cells 2019;37:593-598.

Androgen-regulated transcription of ESRP2 drives alternative splicing patterns in prostate cancer

Prostate is the most frequent cancer in men. Prostate cancer progression is driven by androgen steroid hormones, and delayed by androgen deprivation therapy (ADT). Androgens control transcription by stimulating androgen receptor (AR) activity, yet also control pre-mRNA splicing through less clear mechanisms. Here we find androgens regulate splicing through AR-mediated transcriptional control of the epithelial-specific splicing regulator ESRP2. Both ESRP2 and its close paralog ESRP1 are highly expressed in primary prostate cancer. Androgen stimulation induces splicing switches in many endogenous ESRP2-controlled mRNA isoforms, including splicing switches correlating with disease progression. ESRP2 expression in clinical prostate cancer is repressed by ADT, which may thus inadvertently dampen epithelial splice programmes. Supporting this, treatment with the AR antagonist bicalutamide (Casodex) induced mesenchymal splicing patterns of genes including FLNB and CTNND1. Our data reveals a new mechanism of splicing control in prostate cancer with important implications for disease progression.

RRAD, IL4I1, CDKN1A, and SERPINE1 genes are potentially co-regulated by NF-κB and p53 transcription factors in cells exposed to high doses of ionizing radiation

BACKGROUND

The cellular response to ionizing radiation involves activation of p53-dependent pathways and activation of the atypical NF-κB pathway. The crosstalk between these two transcriptional networks include (co)regulation of common gene targets. Here we looked for novel genes potentially (co)regulated by p53 and NF-κB using integrative genomics screening in human osteosarcoma U2-OS cells irradiated with a high dose (4 and 10 Gy). Radiation-induced expression in cells with silenced TP53 or RELA (coding the p65 NF-κB subunit) genes was analyzed by RNA-Seq while radiation-enhanced binding of p53 and RelA in putative regulatory regions was analyzed by ChIP-Seq, then selected candidates were validated by qPCR.

RESULTS

We identified a subset of radiation-modulated genes whose expression was affected by silencing of both TP53 and RELA, and a subset of radiation-upregulated genes where radiation stimulated binding of both p53 and RelA. For three genes, namely IL4I1, SERPINE1, and CDKN1A, an antagonistic effect of the TP53 and RELA silencing was consistent with radiation-enhanced binding of both p53 and RelA. This suggested the possibility of a direct antagonistic (co)regulation by both factors: activation by NF-κB and inhibition by p53 of IL4I1, and activation by p53 and inhibition by NF-κB of CDKN1A and SERPINE1. On the other hand, radiation-enhanced binding of both p53 and RelA was observed in a putative regulatory region of the RRAD gene whose expression was downregulated both by TP53 and RELA silencing, which suggested a possibility of direct (co)activation by both factors.

CONCLUSIONS

Four new candidates for genes directly co-regulated by NF-κB and p53 were revealed.

Pro-inflammatory cytokine and high doses of ionizing radiation have similar effects on the expression of NF-kappaB-dependent genes

The NF-κB transcription factors are activated via diverse molecular mechanisms in response to various types of stimuli. A plethora of functions associated with specific sets of target genes could be regulated differentially by this factor, affecting cellular response to stress including an anticancer treatment. Here we aimed to compare subsets of NF-κB-dependent genes induced in cells stimulated with a pro-inflammatory cytokine and in cells damaged by a high dose of ionizing radiation (4 and 10 Gy). The RelA-containing NF-κB species were activated by the canonical TNFα-induced and the atypical radiation-induced pathways in human osteosarcoma cells. NF-κB-dependent genes were identified using the gene expression profiling (by RNA-Seq) in cells with downregulated RELA combined with the global profiling of RelA binding sites (by ChIP-Seq), with subsequent validation of selected candidates by quantitative PCR. There were 37 NF-κB-dependent protein-coding genes identified: in all cases RelA bound in their regulatory regions upon activation while downregulation of RELA suppressed their stimulus-induced upregulation, which apparently indicated the positive regulation mode. This set of genes included a few "novel" NF-κB-dependent species. Moreover, the evidence for possible negative regulation of ATF3 gene by NF-κB was collected. The kinetics of the NF-κB activation was slower in cells exposed to radiation than in cytokine-stimulated ones. However, subsets of NF-κB-dependent genes upregulated by both types of stimuli were essentially the same. Hence, one should expect that similar cellular processes resulting from activation of the NF-κB pathway could be induced in cells responding to pro-inflammatory cytokines and in cells where so-called "sterile inflammation" response was initiated by radiation-induced damage.

Abundant and equipotent founder cells establish and maintain acute lymphoblastic leukaemia

High frequencies of blasts in primary acute lymphoblastic leukaemia (ALL) samples have the potential to induce leukaemia and to engraft mice. However, it is unclear how individual ALL cells each contribute to drive leukaemic development in a bulk transplant and the extent to which these blasts vary functionally. We used cellular barcoding as a fate mapping tool to track primograft ALL blasts in vivo. Our results show that high numbers of ALL founder cells contribute at similar frequencies to leukaemic propagation over serial transplants, without any clear evidence of clonal succession. These founder cells also exhibit equal capacity to home and engraft to different organs, although stochastic processes may alter the composition in restrictive niches. Our findings enhance the stochastic stem cell model of ALL by demonstrating equal functional abilities of singular ALL blasts and show that successful treatment strategies must eradicate the entire leukaemic cell population.

Transcriptomics and proteomics show that selenium affects inflammation, cytoskeleton, and cancer pathways in human rectal biopsies

Epidemiologic studies highlight the potential role of dietary selenium (Se) in colorectal cancer prevention. Our goal was to elucidate whether expression of factors crucial for colorectal homoeostasis is affected by physiologic differences in Se status. Using transcriptomics and proteomics followed by pathway analysis, we identified pathways affected by Se status in rectal biopsies from 22 healthy adults, including 11 controls with optimal status (mean plasma Se = 1.43 μM) and 11 subjects with suboptimal status (mean plasma Se = 0.86 μM). We observed that 254 genes and 26 proteins implicated in cancer (80%), immune function and inflammatory response (40%), cell growth and proliferation (70%), cellular movement, and cell death (50%) were differentially expressed between the 2 groups. Expression of 69 genes, including selenoproteins W1 and K, which are genes involved in cytoskeleton remodelling and transcription factor NFκB signaling, correlated significantly with Se status. Integrating proteomics and transcriptomics datasets revealed reduced inflammatory and immune responses and cytoskeleton remodelling in the suboptimal Se status group. This is the first study combining omics technologies to describe the impact of differences in Se status on colorectal expression patterns, revealing that suboptimal Se status could alter inflammatory signaling and cytoskeleton in human rectal mucosa and so influence cancer risk.-Méplan, C., Johnson, I. T., Polley, A. C. J., Cockell, S., Bradburn, D. M., Commane, D. M., Arasaradnam, R. P., Mulholland, F., Zupanic, A., Mathers, J. C., Hesketh, J. Transcriptomics and proteomics show that selenium affects inflammation, cytoskeleton, and cancer pathways in human rectal biopsies.

Genome-wide ChIP-seq analysis of human TOP2B occupancy in MCF7 breast cancer epithelial cells

We report the whole genome ChIP seq for human TOP2B from MCF7 cells. Using three different peak calling methods, regions of binding were identified in the presence or absence of the nuclear hormone estradiol, as TOP2B has been reported to play a role in ligand-induced transcription. TOP2B peaks were found across the whole genome, 50% of the peaks fell either within a gene or within 5 kb of a transcription start site. TOP2B peaks coincident with gene promoters were less frequently associated with epigenetic features marking active promoters in estradiol treated than in untreated cells. Significantly enriched transcription factor motifs within the DNA sequences underlying the peaks were identified. These included SP1, KLF4, TFAP2A, MYF, REST, CTCF, ESR1 and ESR2. Gene ontology analysis of genes associated with TOP2B peaks found neuronal development terms including axonogenesis and axon guidance were significantly enriched. In the absence of functional TOP2B there are errors in axon guidance in the zebrafish eye. Specific heparin sulphate structures are involved in retinal axon targeting. The glycosaminoglycan biosynthesis–heparin sulphate/heparin pathway is significantly enriched in the TOP2B gene ontology analysis, suggesting changes in this pathway in the absence of TOP2B may cause the axon guidance faults.

Summary: Gene ontology enrichment analysis of genes associated with human TOP2B peaks, identified by whole genome ChIP seq used to identify regions of binding, highlighted a number of processes in neuronal development including axonogenesis and axon guidance.

A distributed computational search strategy for the identification of diagnostics targets: Application to finding aptamer targets for methicillin-resistant staphylococci

The rapid and cost-effective identification of bacterial species is crucial, especially for clinical diagnosis and treatment. Peptide aptamers have been shown to be valuable for use as a component of novel, direct detection methods. These small peptides have a number of advantages over antibodies, including greater specificity and longer shelf life. These properties facilitate their use as the detector components of biosensor devices. However, the identification of suitable aptamer targets for particular groups of organisms is challenging. We present a semi-automated processing pipeline for the identification of candidate aptamer targets from whole bacterial genome sequences. The pipeline can be configured to search for protein sequence fragments that uniquely identify a set of strains of interest. The system is also capable of identifying additional organisms that may be of interest due to their possession of protein fragments in common with the initial set. Through the use of Cloud computing technology and distributed databases, our system is capable of scaling with the rapidly growing genome repositories, and consequently of keeping the resulting data sets up-to-date. The system described is also more generically applicable to the discovery of specific targets for other diagnostic approaches such as DNA probes, PCR primers and antibodies.

Contribution of natural antisense transcription to an endogenous siRNA signature in human cells

BACKGROUND

Eukaryotic cells express a complex layer of noncoding RNAs. An intriguing family of regulatory RNAs includes transcripts from the opposite strand of protein coding genes, so called natural antisense transcripts (NATs). Here, we test the hypothesis that antisense transcription triggers RNA interference and gives rise to endogenous short RNAs (endo-siRNAs).

RESULTS

We used cloned human embryonic kidney cells (HEK293) followed by short RNAseq to investigate the small genic RNA transcriptome. 378 genes gave rise to short RNA reads that mapped to exons of RefSeq genes. The length profile of short RNAs showed a broad peak of 20-24 nucleotides, indicative of endo-siRNAs. Collapsed reads mapped predominantly to the first and the last exon of genes (74%). RNAs reads were intersected with sequences occupied by RNAPII or bound to Argonaute (AGO1 by crosslinking, ligation, and sequencing of hybrids, CLASH). In the first exon, 94% of the reads correlated with RNAPII occupancy with an average density of 130 (relative units); this decreased to 65%/20 in middle exons and 54%/12 in the last exon. CLASH reads mapping to multi-exon genes showed little distribution bias with an average of about 5 CLASH reads overlapping with 60% of the endo-siRNA reads. However, endo-siRNAs (21-25 nt) intersecting with CLASH reads were enriched at the 5'end and decreased towards the 3'end.We then investigated the 378 genes with particular focus on features indicative for short RNA production; however, found that endo-siRNA numbers did not correlate with gene structures that favor convergent transcription. In contrast, our gene set was found notably over-represented in the NATsDB sense/antisense group as compared to non-overlapping and non-bidirectional groups. Moreover, read counts showed no correlation with the steady-state levels of the related mRNAs and the pattern of endo-siRNAs proved reproducible after an induced mutagenic insult.

CONCLUSIONS

Our results suggest that antisense transcripts contribute to low levels of endo-siRNAs in fully differentiated human cells. A characteristic endo-siRNA footprint is being produced at sites of RNAPII transcription which is also related to AGO1. This endo-siRNA signature represents an intriguing finding and its reproducibility suggests that the production of endo-siRNAs is a regulated process with potential homoeostatic impact.