Noninvasive stratification of nonalcoholic fatty liver disease by whole transcriptome cell-free mRNA characterization.

Hepatic fibrosis stage is the most important determinant of outcomes in patients with nonalcoholic fatty liver disease (NAFLD). There is an urgent need for noninvasive tests that can accurately stage fibrosis and determine efficacy of interventions. Here, we describe a novel cell-free (cf)-mRNA sequencing approach that can accurately and reproducibly profile low levels of circulating mRNAs and evaluate the feasibility of developing a cf-mRNA-based NAFLD fibrosis classifier. Using separate discovery and validation cohorts with biopsy-confirmed NAFLD (n = 176 and 59, respectively) and healthy subjects (n = 23), we performed serum cf-mRNA RNA-Seq profiling. Differential expression analysis identified 2,498 dysregulated genes between patients with NAFLD and healthy subjects and 134 fibrosis-associated genes in patients with NAFLD. Comparison between cf-mRNA and liver tissue transcripts revealed significant overlap of fibrosis-associated genes and pathways indicating that the circulating cf-mRNA transcriptome reflects molecular changes in the livers of patients with NAFLD. In particular, metabolic and immune pathways reflective of known underlying steatosis and inflammation were highly dysregulated in the cf-mRNA profile of patients with advanced fibrosis. Finally, we used an elastic net ordinal logistic model to develop a classifier that predicts clinically significant fibrosis (F2-F4). In an independent cohort, the cf-mRNA classifier was able to identify 50% of patients with at least 90% probability of clinically significant fibrosis. We demonstrate a novel and robust cf-mRNA-based RNA-Seq platform for noninvasive identification of diverse hepatic molecular disruptions and for fibrosis staging with promising potential for clinical trials and clinical practice.NEW & NOTEWORTHY This work is the first study, to our knowledge, to utilize circulating cell-free mRNA sequencing to develop an NAFLD diagnostic classifier.

Local sequence assembly reveals a high-resolution profile of somatic structural variations in 97 cancer genomes.

Genomic structural variations (SVs) are pervasive in many types of cancers. Characterizing their underlying mechanisms and potential molecular consequences is crucial for understanding the basic biology of tumorigenesis. Here, we engineered a local assembly-based algorithm (laSV) that detects SVs with high accuracy from paired-end high-throughput genomic sequencing data and pinpoints their breakpoints at single base-pair resolution. By applying laSV to 97 tumor-normal paired genomic sequencing datasets across six cancer types produced by The Cancer Genome Atlas Research Network, we discovered that non-allelic homologous recombination is the primary mechanism for generating somatic SVs in acute myeloid leukemia. This finding contrasts with results for the other five types of solid tumors, in which non-homologous end joining and microhomology end joining are the predominant mechanisms. We also found that the genes recursively mutated by single nucleotide alterations differed from the genes recursively mutated by SVs, suggesting that these two types of genetic alterations play different roles during cancer progression. We further characterized how the gene structures of the oncogene JAK1 and the tumor suppressors KDM6A and RB1 are affected by somatic SVs and discussed the potential functional implications of intergenic SVs.

Forward genetic screen of human transposase genomic rearrangements.

BACKGROUND: Numerous human genes encode potentially active DNA transposases or recombinases, but our understanding of their functions remains limited due to shortage of methods to profile their activities on endogenous genomic substrates. RESULTS: To enable functional analysis of human transposase-derived genes, we combined forward chemical genetic hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1) screening with massively parallel paired-end DNA sequencing and structural variant genome assembly and analysis. Here, we report the HPRT1 mutational spectrum induced by the human transposase PGBD5, including PGBD5-specific signal sequences (PSS) that serve as potential genomic rearrangement substrates. CONCLUSIONS: The discovered PSS motifs and high-throughput forward chemical genomic screening approach should prove useful for the elucidation of endogenous genome remodeling activities of PGBD5 and other domesticated human DNA transposases and recombinases.

TEMP: a computational method for analyzing transposable element polymorphism in populations.

Insertions and excisions of transposable elements (TEs) affect both the stability and variability of the genome. Studying the dynamics of transposition at the population level can provide crucial insights into the processes and mechanisms of genome evolution. Pooling genomic materials from multiple individuals followed by high-throughput sequencing is an efficient way of characterizing genomic polymorphisms in a population. Here we describe a novel method named TEMP, specifically designed to detect TE movements present with a wide range of frequencies in a population. By combining the information provided by pair-end reads and split reads, TEMP is able to identify both the presence and absence of TE insertions in genomic DNA sequences derived from heterogeneous samples; accurately estimate the frequencies of transposition events in the population and pinpoint junctions of high frequency transposition events at nucleotide resolution. Simulation data indicate that TEMP outperforms other algorithms such as PoPoolationTE, RetroSeq, VariationHunter and GASVPro. TEMP also performs well on whole-genome human data derived from the 1000 Genomes Project. We applied TEMP to characterize the TE frequencies in a wild Drosophila melanogaster population and study the inheritance patterns of TEs during hybrid dysgenesis. We also identified sequence signatures of TE insertion and possible molecular effects of TE movements, such as altered gene expression and piRNA production. TEMP is freely available at github: https://github.com/JialiUMassWengLab/TEMP.git.

Sequence features and chromatin structure around the genomic regions bound by 119 human transcription factors.

Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) has become the dominant technique for mapping transcription factor (TF) binding regions genome-wide. We performed an integrative analysis centered around 457 ChIP-seq data sets on 119 human TFs generated by the ENCODE Consortium. We identified highly enriched sequence motifs in most data sets, revealing new motifs and validating known ones. The motif sites (TF binding sites) are highly conserved evolutionarily and show distinct footprints upon DNase I digestion. We frequently detected secondary motifs in addition to the canonical motifs of the TFs, indicating tethered binding and cobinding between multiple TFs. We observed significant position and orientation preferences between many cobinding TFs. Genes specifically expressed in a cell line are often associated with a greater occurrence of nearby TF binding in that cell line. We observed cell-line-specific secondary motifs that mediate the binding of the histone deacetylase HDAC2 and the enhancer-binding protein EP300. TF binding sites are located in GC-rich, nucleosome-depleted, and DNase I sensitive regions, flanked by well-positioned nucleosomes, and many of these features show cell type specificity. The GC-richness may be beneficial for regulating TF binding because, when unoccupied by a TF, these regions are occupied by nucleosomes in vivo. We present the results of our analysis in a TF-centric web repository Factorbook (http://factorbook.org) and will continually update this repository as more ENCODE data are generated.

Factorbook.org: a Wiki-based database for transcription factor-binding data generated by the ENCODE consortium.

The Encyclopedia of DNA Elements (ENCODE) consortium aims to identify all functional elements in the human genome including transcripts, transcriptional regulatory regions, along with their chromatin states and DNA methylation patterns. The ENCODE project generates data utilizing a variety of techniques that can enrich for regulatory regions, such as chromatin immunoprecipitation (ChIP), micrococcal nuclease (MNase) digestion and DNase I digestion, followed by deeply sequencing the resulting DNA. As part of the ENCODE project, we have developed a Web-accessible repository accessible at http://factorbook.org. In Wiki format, factorbook is a transcription factor (TF)-centric repository of all ENCODE ChIP-seq datasets on TF-binding regions, as well as the rich analysis results of these data. In the first release, factorbook contains 457 ChIP-seq datasets on 119 TFs in a number of human cell lines, the average profiles of histone modifications and nucleosome positioning around the TF-binding regions, sequence motifs enriched in the regions and the distance and orientation preferences between motif sites.

Circulating cell-free messenger RNA secretome characterization of primary sclerosing cholangitis.

BACKGROUND: Primary sclerosing cholangitis (PSC) is a rare chronic cholestatic liver disease characterized by multifocal bile duct strictures. To date, underlying molecular mechanisms of PSC remain unclear, and therapeutic options are limited. METHODS: We performed cell-free messenger RNA (cf-mRNA) sequencing to characterize the circulating transcriptome of PSC and noninvasively investigate potentially bioactive signals that are associated with PSC. Serum cf-mRNA profiles were compared among 50 individuals with PSC, 20 healthy controls, and 235 individuals with NAFLD. Tissue and cell type-of-origin genes that are dysregulated in subjects with PSC were evaluated. Subsequently, diagnostic classifiers were developed using PSC dysregulated cf-mRNA genes. RESULTS: Differential expression analysis of the cf-mRNA transcriptomes of PSC and healthy controls resulted in identification of 1407 dysregulated genes. Furthermore, differentially expressed genes between PSC and healthy controls or NAFLD shared common genes known to be involved in liver pathophysiology. In particular, genes from liver- and specific cell type-origin, including hepatocyte, HSCs, and KCs, were highly abundant in cf-mRNA of subjects with PSC. Gene cluster analysis revealed that liver-specific genes dysregulated in PSC form a distinct cluster, which corresponded to a subset of the PSC subject population. Finally, we developed a cf-mRNA diagnostic classifier using liver-specific genes that discriminated PSC from healthy control subjects using gene transcripts of liver origin. CONCLUSIONS: Blood-based whole-transcriptome cf-mRNA profiling revealed high abundance of liver-specific genes in sera of subjects with PSC, which may be used to diagnose patients with PSC. We identified several unique cf-mRNA profiles of subjects with PSC. These findings may also have utility for noninvasive molecular stratification of subjects with PSC for pharmacotherapy safety and response studies.

TGF-ß1 enhances Kv2.1 potassium channel protein expression and promotes maturation of cerebellar granule neurons.

Members of the transforming growth factor-ß (TGF-ß) family of cytokines are involved in diverse physiological processes. Although TGF-ß is known to play multiple roles in the mammalian central nervous system (CNS), its role in neuronal development has not been explored. We have studied the effects of TGF-ß1 on the electrophysiological properties and maturation of rat primary cerebellar granule neurons (CGNs). We report that incubation with TGF-ß1 increased delayed rectifier potassium current (I(K) ) amplitudes in a dose- and time-dependent manner, but did not affect the kinetic properties of the channel. Exposure to TGF-ß1 (20 ng/ml) for 36 h led to a 37.2% increase in I(K) amplitudes. There was no significant change in mRNA levels for the key Kv2.1 channel protein, but translation blockade abolished the increase in protein levels and channel activity, arguing that TGF-ß1 increases I(K) amplitudes by upregulating translation of the Kv2.1 channel protein. Although TGF-ß1 treatment did not affect the activity of protein kinase A (PKA), and constitutive activation of PKA with forskolin failed to increase I(K) amplitudes, inhibition of PKA prevented channel upregulation, demonstrating that basal PKA activity is required for TGF-ß1 stimulation of I(K) channel activity. TGF-ß1 also promoted the expression of the γ-aminobutyric acid (GABA(A) ) receptor α6 subunit, a marker of mature CGNs, and calcium influx during depolarizing stimuli was reduced by TGF-ß1. The effects of TGF-ß1 were only observed during a narrow developmental time-window, and were lost as CGNs matured. These findings suggest that TGF-ß1 upregulates K(+) channel expression and I(K) currents and thereby promotes CGN maturation.

The antidepressant citalopram inhibits delayed rectifier outward K⁺ current in mouse cortical neurons.

Citalopram, a selective serotonin (5-HT) reuptake inhibitor (SSRI) as well as an antidepressant, is thought to exert its effects by increasing synaptic 5-HT levels. However, few studies have addressed the possibility that citalopram has other molecular mechanisms of action. We examined the effects of citalopram on delayed rectifier outward K(+) current (I(K) ) in mouse cortical neurons. Extracellular citalopram reversibly inhibited I(K) in a dose-dependent manner and significantly shifted both steady-state activation and inactivation curves toward hyperpolarization. Neither 5-HT itself nor antagonists of 5-HT and dopamine receptors could abolish citalopram-induced inhibition of I(K) . In addition, intracellular application of GTPγ-S similarly failed to prevent the inhibition of I(K) by citalopram. When applied intracellularly, citalopram had no effect on I(K) and did not influence the reduction of I(K) induced by extracellular citalopram. The effect of citalopram was use dependent, but not frequency dependent, and it did not require channel opening. Electrophysiological recordings in acute cortical slice showed that citalopram significantly reduced the action potential (AP) firing frequency of cortical neurons and increased action potential duration (APD). The selective Kv2.1 subunit blocker Jingzhaotoxin-III (JZTX-III) did not abolish citalopram-induced I(K) inhibition. Transfection of HEK293 cells with Kv2.1 or Kv2.2 constructs indicated that citalopram mainly inhibited Kv2.2 current. We suggest that citalopram-induced inhibition of I(K) in mouse cortical neurons is independent of G-protein-coupled receptors and might exert its antidepressant effects by enhancing presynaptic efficiency. Our results may help to explain some of the unknown therapeutic effects of citalopram.

Survey of extracellular communication of systemic and organ-specific inflammatory responses through cell free messenger RNA profiling in mice.

BACKGROUND: Inflammatory and immune responses are essential and dynamic biological processes that protect the body against acute and chronic adverse stimuli. While conventional protein markers have been used to evaluate systemic inflammatory response, the immunological response to stimulation is complex and involves modulation of a large set of genes and interacting signalling pathways of innate and adaptive immune systems. There is a need for a non-invasive tool that can comprehensively evaluate and monitor molecular dysregulations associated with inflammatory and immune responses in circulation and in inaccessible solid organs. METHODS: Here we utilized cell-free messenger RNA (cf-mRNA) RNA-Seq whole transcriptome profiling and computational biology to temporally assess lipopolysaccharide (LPS) induced and JAK inhibitor modulated inflammatory and immune responses in mouse plasma samples. FINDINGS: Cf-mRNA profiling displayed a pattern of systemic immune responses elicited by LPS and dysregulation of associated pathways. Moreover, attenuation of several inflammatory pathways, including STAT and interferon pathways, were observed following the treatment of JAK inhibitor. We further identified the dysregulation of liver-specific transcripts in cf-mRNA which reflected changes in the gene-expression pattern in this generally inaccessible biological compartment. INTERPRETATION: Using a preclinical mouse model, we demonstrated the potential of plasma cf-mRNA profiling for systemic and organ-specific characterization of drug-induced molecular alterations that are associated with inflammatory and immune responses. FUNDING: Molecular Stethoscope.

Brain natriuretic peptide modulates the delayed rectifier outward K(+) current and promotes the proliferation of mouse Schwann cells.

Brain natriuretic peptide (BNP) may act as a neuromodulator via its associated receptors (natriuretic peptide receptors, NPRs) in the central nervous system (CNS), but few studies have reported its activity in the peripheral nervous system (PNS). In this study, we observed that BNP increased the tetraethylammonium chloride (TEA)-sensitive delayed rectifier outward potassium current (I(K)) in mouse Schwann cells (SCs) using whole-cell recording techniques. At concentrations of 1-100 nM, BNP reversibly activated I(K) in a dose-dependent manner, with modulating its steady-state activation and inactivation properties. The effect of BNP on I(K) was abolished by preincubation with the specific antagonist of NPR-A, and could not be mimicked by application of NPR-C agonist. These results were supported by immunocytochemical findings indicating that NPR-A was expressed in SCs. The application of 8-Br-guanosine 3',5'-monophosphate (8-Br-cGMP) mimicked the effect of BNP on I(K), but BNP was unable to further increase I(K) after the application of cyclic guanosine monophosphate (cGMP). Genistein blocked I(K) and also completely eliminated the effects of BNP and cGMP on I(K). The selective K(V)2.1 subunit blocker, Jingzhaotoxin-III (JZTX-III), reduced I(K) amplitude by 30%, but did not abolish the increase effect of BNP on I(K) amplitude. In addition, BNP significantly stimulated SCs proliferation and this effect could be partly inhibited by TEA. Together these results suggest that BNP modulated I(K) probably via cGMP- and tyrosine kinase-dependent pathways by activation of NPR-A. This effect of BNP on I(K) in SCs might partly explain its effect on cell proliferation.

Noninvasive characterization of Alzheimer's disease by circulating, cell-free messenger RNA next-generation sequencing.

The lack of accessible noninvasive tools to examine the molecular alterations occurring in the brain limits our understanding of the causes and progression of Alzheimer's disease (AD), as well as the identification of effective therapeutic strategies. Here, we conducted a comprehensive profiling of circulating, cell-free messenger RNA (cf-mRNA) in plasma of 126 patients with AD and 116 healthy controls of similar age. We identified 2591 dysregulated genes in the cf-mRNA of patients with AD, which are enriched in biological processes well known to be associated with AD. Dysregulated genes included brain-specific genes and resembled those identified to be dysregulated in postmortem AD brain tissue. Furthermore, we identified disease-relevant circulating gene transcripts that correlated with the severity of cognitive impairment. These data highlight the potential of high-throughput cf-mRNA sequencing to evaluate AD-related pathophysiological alterations in the brain, leading to precision healthcare solutions that could improve AD patient management.

Non-invasive characterization of human bone marrow stimulation and reconstitution by cell-free messenger RNA sequencing.

Circulating cell-free mRNA (cf-mRNA) holds great promise as a non-invasive diagnostic biomarker. However, cf-mRNA composition and its potential clinical applications remain largely unexplored. Here we show, using Next Generation Sequencing-based profiling, that cf-mRNA is enriched in transcripts derived from the bone marrow compared to circulating cells. Further, longitudinal studies involving bone marrow ablation followed by hematopoietic stem cell transplantation in multiple myeloma and acute myeloid leukemia patients indicate that cf-mRNA levels reflect the transcriptional activity of bone marrow-resident hematopoietic lineages during bone marrow reconstitution. Mechanistically, stimulation of specific bone marrow cell populations in vivo using growth factor pharmacotherapy show that cf-mRNA reflects dynamic functional changes over time associated with cellular activity. Our results shed light on the biology of the circulating transcriptome and highlight the potential utility of cf-mRNA to non-invasively monitor bone marrow involved pathologies.

Author Correction: PGBD5 promotes site-specific oncogenic mutations in human tumors.

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

PGBD5 promotes site-specific oncogenic mutations in human tumors.

Genomic rearrangements are a hallmark of human cancers. Here, we identify the piggyBac transposable element derived 5 (PGBD5) gene as encoding an active DNA transposase expressed in the majority of childhood solid tumors, including lethal rhabdoid tumors. Using assembly-based whole-genome DNA sequencing, we found previously undefined genomic rearrangements in human rhabdoid tumors. These rearrangements involved PGBD5-specific signal (PSS) sequences at their breakpoints and recurrently inactivated tumor-suppressor genes. PGBD5 was physically associated with genomic PSS sequences that were also sufficient to mediate PGBD5-induced DNA rearrangements in rhabdoid tumor cells. Ectopic expression of PGBD5 in primary immortalized human cells was sufficient to promote cell transformation in vivo. This activity required specific catalytic residues in the PGBD5 transposase domain as well as end-joining DNA repair and induced structural rearrangements with PSS breakpoints. These results define PGBD5 as an oncogenic mutator and provide a plausible mechanism for site-specific DNA rearrangements in childhood and adult solid tumors.