Exploring differential exon usage via short- and long-read RNA sequencing strategies.

Alternative splicing produces various mRNAs, and thereby various protein products, from one gene, impacting a wide range of cellular activities. However, accurate reconstruction and quantification of full-length transcripts using short-reads is limited, due to their length. Long-reads sequencing technologies may provide a solution by sequencing full-length transcripts. We explored the use of both Illumina short-reads and two long Oxford Nanopore Technology (cDNA and Direct RNA) RNA-Seq reads for detecting global differential splicing during mouse embryonic stem cell differentiation, applying several bioinformatics strategies: gene-based, isoform-based and exon-based. We detected the strongest similarity among the sequencing platforms at the gene level compared to exon-based and isoform-based. Furthermore, the exon-based strategy discovered many differential exon usage (DEU) events, mostly in a platform-dependent manner and in non-differentially expressed genes. Thus, the platforms complemented each other in the ability to detect DEUs (i.e. long-reads exhibited an advantage in detecting DEUs at the UTRs, and short-reads detected more DEUs). Exons within 20 genes, detected in one or more platforms, were here validated by PCR, including key differentiation genes, such as Mdb3 and Aplp1. We provide an important analysis resource for discovering transcriptome changes during stem cell differentiation and insights for analysing such data.

Whole-genome sequencing reveals that variants in the Interleukin 18 Receptor Accessory Protein 3'UTR protect against ALS.

The noncoding genome is substantially larger than the protein-coding genome but has been largely unexplored by genetic association studies. Here, we performed region-based rare variant association analysis of >25,000 variants in untranslated regions of 6,139 amyotrophic lateral sclerosis (ALS) whole genomes and the whole genomes of 70,403 non-ALS controls. We identified interleukin-18 receptor accessory protein (IL18RAP) 3' untranslated region (3'UTR) variants as significantly enriched in non-ALS genomes and associated with a fivefold reduced risk of developing ALS, and this was replicated in an independent cohort. These variants in the IL18RAP 3'UTR reduce mRNA stability and the binding of double-stranded RNA (dsRNA)-binding proteins. Finally, the variants of the IL18RAP 3'UTR confer a survival advantage for motor neurons because they dampen neurotoxicity of human induced pluripotent stem cell (iPSC)-derived microglia bearing an ALS-associated expansion in C9orf72, and this depends on NF-κB signaling. This study reveals genetic variants that protect against ALS by reducing neuroinflammation and emphasizes the importance of noncoding genetic association studies.

Comparison of seven single cell whole genome amplification commercial kits using targeted sequencing.

Advances in whole genome amplification (WGA) techniques enable understanding of the genomic sequence at a single cell level. Demand for single cell dedicated WGA kits (scWGA) has led to the development of several commercial kit. To this point, no robust comparison of all available kits was performed. Here, we benchmark an economical assay, comparing all commercially available scWGA kits. Our comparison is based on targeted sequencing of thousands of genomic loci, including highly mutable regions, from a large cohort of human single cells. Using this approach we have demonstrated the superiority of Ampli1 in genome coverage and of RepliG in reduced error rate. In summary, we show that no single kit is optimal across all categories, highlighting the need for a dedicated kit selection in accordance with experimental requirements.

Highly conserved and cis-acting lncRNAs produced from paralogous regions in the center of HOXA and HOXB clusters in the endoderm lineage.

Long noncoding RNAs (lncRNAs) have been shown to play important roles in gene regulatory networks acting in early development. There has been rapid turnover of lncRNA loci during vertebrate evolution, with few human lncRNAs conserved beyond mammals. The sequences of these rare deeply conserved lncRNAs are typically not similar to each other. Here, we characterize HOXA-AS3 and HOXB-AS3, lncRNAs produced from the central regions of the HOXA and HOXB clusters. Sequence-similar orthologs of both lncRNAs are found in multiple vertebrate species and there is evident sequence similarity between their promoters, suggesting that the production of these lncRNAs predates the duplication of the HOX clusters at the root of the vertebrate lineage. This conservation extends to similar expression patterns of the two lncRNAs, in particular in cells transiently arising during early development or in the adult colon. Functionally, the RNA products of HOXA-AS3 and HOXB-AS3 regulate the expression of their overlapping HOX5-7 genes both in HT-29 cells and during differentiation of human embryonic stem cells. Beyond production of paralogous protein-coding and microRNA genes, the regulatory program in the HOX clusters therefore also relies on paralogous lncRNAs acting in restricted spatial and temporal windows of embryonic development and cell differentiation.

Community-Wide Experimental Evaluation of the PROSS Stability-Design Method.

Recent years have seen a dramatic improvement in protein-design methodology. Nevertheless, most methods demand expert intervention, limiting their widespread adoption. By contrast, the PROSS algorithm for improving protein stability and heterologous expression levels has been successfully applied to a range of challenging enzymes and binding proteins. Here, we benchmark the application of PROSS as a stand-alone tool for protein scientists with no or limited experience in modeling. Twelve laboratories from the Protein Production and Purification Partnership in Europe (P4EU) challenged the PROSS algorithm with 14 unrelated protein targets without support from the PROSS developers. For each target, up to six designs were evaluated for expression levels and in some cases, for thermal stability and activity. In nine targets, designs exhibited increased heterologous expression levels either in prokaryotic and/or eukaryotic expression systems under experimental conditions that were tailored for each target protein. Furthermore, we observed increased thermal stability in nine of ten tested targets. In two prime examples, the human Stem Cell Factor (hSCF) and human Cadherin-Like Domain (CLD12) from the RET receptor, the wild type proteins were not expressible as soluble proteins in E. coli, yet the PROSS designs exhibited high expression levels in E. coli and HEK293 cells, respectively, and improved thermal stability. We conclude that PROSS may improve stability and expressibility in diverse cases, and that improvement typically requires target-specific expression conditions. This study demonstrates the strengths of community-wide efforts to probe the generality of new methods and recommends areas for future research to advance practically useful algorithms for protein science.

Human genetics and neuropathology suggest a link between miR-218 and amyotrophic lateral sclerosis pathophysiology.

Motor neuron-specific microRNA-218 (miR-218) has recently received attention because of its roles in mouse development. However, miR-218 relevance to human motor neuron disease was not yet explored. Here, we demonstrate by neuropathology that miR-218 is abundant in healthy human motor neurons. However, in amyotrophic lateral sclerosis (ALS) motor neurons, miR-218 is down-regulated and its mRNA targets are reciprocally up-regulated (derepressed). We further identify the potassium channel Kv10.1 as a new miR-218 direct target that controls neuronal activity. In addition, we screened thousands of ALS genomes and identified six rare variants in the human miR-218-2 sequence. miR-218 gene variants fail to regulate neuron activity, suggesting the importance of this small endogenous RNA for neuronal robustness. The underlying mechanisms involve inhibition of miR-218 biogenesis and reduced processing by DICER. Therefore, miR-218 activity in motor neurons may be susceptible to failure in human ALS, suggesting that miR-218 may be a potential therapeutic target in motor neuron disease.

LATS1 and LATS2 suppress breast cancer progression by maintaining cell identity and metabolic state.

Deregulated activity of LArge Tumor Suppressor (LATS) tumor suppressors has broad implications on cellular and tissue homeostasis. We examined the consequences of down-regulation of either LATS1 or LATS2 in breast cancer. Consistent with their proposed tumor suppressive roles, expression of both paralogs was significantly down-regulated in human breast cancer, and loss of either paralog accelerated mammary tumorigenesis in mice. However, each paralog had a distinct impact on breast cancer. Thus, LATS2 depletion in luminal B tumors resulted in metabolic rewiring, with increased glycolysis and reduced peroxisome proliferator-activated receptor γ (PPARγ) signaling. Furthermore, pharmacological activation of PPARγ elicited LATS2-dependent death in luminal B-derived cells. In contrast, LATS1 depletion augmented cancer cell plasticity, skewing luminal B tumors towards increased expression of basal-like features, in association with increased resistance to hormone therapy. Hence, these two closely related paralogs play distinct roles in protection against breast cancer; tumors with reduced expression of either LATS1 or LATS2 may rewire signaling networks differently and thus respond differently to anticancer treatments.

p53 is essential for DNA methylation homeostasis in naïve embryonic stem cells, and its loss promotes clonal heterogeneity.

DNA methylation is a key regulator of embryonic stem cell (ESC) biology, dynamically changing between naïve, primed, and differentiated states. The p53 tumor suppressor is a pivotal guardian of genomic stability, but its contributions to epigenetic regulation and stem cell biology are less explored. We report that, in naïve mouse ESCs (mESCs), p53 restricts the expression of the de novo DNA methyltransferases Dnmt3a and Dnmt3b while up-regulating Tet1 and Tet2, which promote DNA demethylation. The DNA methylation imbalance in p53-deficient (p53-/-) mESCs is the result of augmented overall DNA methylation as well as increased methylation landscape heterogeneity. In differentiating p53-/- mESCs, elevated methylation persists, albeit more mildly. Importantly, concomitant with DNA methylation heterogeneity, p53-/- mESCs display increased cellular heterogeneity both in the "naïve" state and upon induced differentiation. This impact of p53 loss on 5-methylcytosine (5mC) heterogeneity was also evident in human ESCs and mouse embryos in vivo. Hence, p53 helps maintain DNA methylation homeostasis and clonal homogeneity, a function that may contribute to its tumor suppressor activity.

Sensing Cellular Metabolic Activity via a Molecular-Controlled Semiconductor Resistor.

Over the last decade, we have developed a molecular-controlled semiconductor resistor (MOCSER) device that is highly sensitive to variations in its surface potentials. This device was applied as a molecular sensor both in the gas phase and in solutions. The device is based on an AlGaAs/GaAs structure. In the current work, we developed an electronic biosensor for real-time, label-free monitoring of cellular metabolic activity by culturing HeLa cells directly on top of the device's conductive channel. Several properties of GaAs make it attractive for developing biosensors, among others its high electron mobility and ability to control the device's properties by proper epitaxial growing. However, GaAs is very reactive and sensitive to oxidation in aqueous solutions, and its arsenic residues are highly toxic. Nevertheless, we have managed to overcome this inherent chemical instability by developing a surface-protecting layer using polymerized (3-mercaptopropyl)-trimethoxysilane (MPTMS). To improve cell adhesion and biocompatibility, the MPTMS-coated devices were further modified with an additional layer of (3-aminopropyl)-trimethoxysilane (APTMS). HeLa cells were found to grow successfully on these devices, and MOCSER devices cultured with these cells were stable and sensitive to cellular metabolic activity. The sensitivity of the MOCSER device results from the sensing of extracellular acidification in the microenvironment of the cell-MOCSER interspace. We have found that this sensitivity is maintained only when the device is partially covered with the cellular layer, whereas at full coverage the sensitivity is lost. This phenomenon is related to the negatively charged cellular membrane potentials that lead to a reduction in the channel's conductivity. We propose that the coated MOCSER device can be applied for real-time and continuous monitoring of cellular viability and activity.

Cap-independent translation by DAP5 controls cell fate decisions in human embryonic stem cells.

Multiple transcriptional and epigenetic changes drive differentiation of embryonic stem cells (ESCs). This study unveils an additional level of gene expression regulation involving noncanonical, cap-independent translation of a select group of mRNAs. This is driven by death-associated protein 5 (DAP5/eIF4G2/NAT1), a translation initiation factor mediating IRES-dependent translation. We found that the DAP5 knockdown from human ESCs (hESCs) resulted in persistence of pluripotent gene expression, delayed induction of differentiation-associated genes in different cell lineages, and defective embryoid body formation. The latter involved improper cellular organization, lack of cavitation, and enhanced mislocalized apoptosis. RNA sequencing of polysome-associated mRNAs identified candidates with reduced translation efficiency in DAP5-depleted hESCs. These were enriched in mitochondrial proteins involved in oxidative respiration, a pathway essential for differentiation, the significance of which was confirmed by the aberrant mitochondrial morphology and decreased oxidative respiratory activity in DAP5 knockdown cells. Further analysis identified the chromatin modifier HMGN3 as a cap-independent DAP5 translation target whose knockdown resulted in defective differentiation. Thus, DAP5-mediated translation of a specific set of proteins is critical for the transition from pluripotency to differentiation, highlighting the importance of cap-independent translation in stem cell fate decisions.

A generic, cost-effective, and scalable cell lineage analysis platform.

Advances in single-cell genomics enable commensurate improvements in methods for uncovering lineage relations among individual cells. Current sequencing-based methods for cell lineage analysis depend on low-resolution bulk analysis or rely on extensive single-cell sequencing, which is not scalable and could be biased by functional dependencies. Here we show an integrated biochemical-computational platform for generic single-cell lineage analysis that is retrospective, cost-effective, and scalable. It consists of a biochemical-computational pipeline that inputs individual cells, produces targeted single-cell sequencing data, and uses it to generate a lineage tree of the input cells. We validated the platform by applying it to cells sampled from an ex vivo grown tree and analyzed its feasibility landscape by computer simulations. We conclude that the platform may serve as a generic tool for lineage analysis and thus pave the way toward large-scale human cell lineage discovery.

Using Synthetic Mouse Spike-In Transcripts to Evaluate RNA-Seq Analysis Tools.

One of the key applications of next-generation sequencing (NGS) technologies is RNA-Seq for transcriptome genome-wide analysis. Although multiple studies have evaluated and benchmarked RNA-Seq tools dedicated to gene level analysis, few studies have assessed their effectiveness on the transcript-isoform level. Alternative splicing is a naturally occurring phenomenon in eukaryotes, significantly increasing the biodiversity of proteins that can be encoded by the genome. The aim of this study was to assess and compare the ability of the bioinformatics approaches and tools to assemble, quantify and detect differentially expressed transcripts using RNA-Seq data, in a controlled experiment. To this end, in vitro synthesized mouse spike-in control transcripts were added to the total RNA of differentiating mouse embryonic bodies, and their expression patterns were measured. This novel approach was used to assess the accuracy of the tools, as established by comparing the observed results versus the results expected of the mouse controlled spiked-in transcripts. We found that detection of differential expression at the gene level is adequate, yet on the transcript-isoform level, all tools tested lacked accuracy and precision.

DNA lesion identity drives choice of damage tolerance pathway in murine cell chromosomes.

DNA-damage tolerance (DDT) via translesion DNA synthesis (TLS) or homology-dependent repair (HDR) functions to bypass DNA lesions encountered during replication, and is critical for maintaining genome stability. Here, we present piggyBlock, a new chromosomal assay that, using piggyBac transposition of DNA containing a known lesion, measures the division of labor between the two DDT pathways. We show that in the absence of DNA damage response, tolerance of the most common sunlight-induced DNA lesion, TT-CPD, is achieved by TLS in mouse embryo fibroblasts. Meanwhile, BP-G, a major smoke-induced DNA lesion, is bypassed primarily by HDR, providing the first evidence for this mechanism being the main tolerance pathway for a biologically important lesion in a mammalian genome. We also show that, far from being a last-resort strategy as it is sometimes portrayed, TLS operates alongside nucleotide excision repair, handling 40% of TT-CPDs in repair-proficient cells. Finally, DDT acts in mouse embryonic stem cells, exhibiting the same pattern­mutagenic TLS included­despite the risk of propagating mutations along all cell lineages. The new method highlights the importance of HDR, and provides an effective tool for studying DDT in mammalian cells.

Dynamic and static maintenance of epigenetic memory in pluripotent and somatic cells.

Stable maintenance of gene regulatory programs is essential for normal function in multicellular organisms. Epigenetic mechanisms, and DNA methylation in particular, are hypothesized to facilitate such maintenance by creating cellular memory that can be written during embryonic development and then guide cell-type-specific gene expression. Here we develop new methods for quantitative inference of DNA methylation turnover rates, and show that human embryonic stem cells preserve their epigenetic state by balancing antagonistic processes that add and remove methylation marks rather than by copying epigenetic information from mother to daughter cells. In contrast, somatic cells transmit considerable epigenetic information to progenies. Paradoxically, the persistence of the somatic epigenome makes it more vulnerable to noise, since random epimutations can accumulate to massively perturb the epigenomic ground state. The rate of epigenetic perturbation depends on the genomic context, and, in particular, DNA methylation loss is coupled to late DNA replication dynamics. Epigenetic perturbation is not observed in the pluripotent state, because the rapid turnover-based equilibrium continuously reinforces the canonical state. This dynamic epigenetic equilibrium also explains how the epigenome can be reprogrammed quickly and to near perfection after induced pluripotency.

Hellebrin and its aglycone form hellebrigenin display similar in vitro growth inhibitory effects in cancer cells and binding profiles to the alpha subunits of the Na+/K+-ATPase.

BACKGROUND: Surface-expressed Na+/K+-ATPase (NaK) has been suggested to function as a non-canonical cardiotonic steroid-binding receptor that activates multiple signaling cascades, especially in cancer cells. By contrast, the current study establishes a clear correlation between the IC50in vitro growth inhibitory concentration in human cancer cells and the Ki for the inhibition of activity of purified human α1ß1 NaK. METHODS: The in vitro growth inhibitory effects of seven cardiac glycosides including five cardenolides (ouabain, digoxin, digitoxin, gitoxin, uzarigenin-rhamnoside, and their respective aglycone forms) and two bufadienolides (gamabufotalin-rhamnoside and hellebrin, and their respective aglycone forms) were determined by means of the MTT colorimetric assay and hellebrigenin-induced cytotoxic effects were visualized by means of quantitative videomicroscopy. The binding affinity of ten of the 14 compounds under study was determined with respect to human α1ß1, α2ß1 and α3ß1 NaK complexes. Lactate releases and oxygen consumption rates were also determined in cancer cells treated with these various cardiac glycosides. RESULTS: Although cardiotonic steroid aglycones usually display weaker binding affinity and in vitro anticancer activity than the corresponding glycoside, the current study demonstrates that the hellebrin / hellebrigenin pair is at odds with respect to this rule. In addition, while some cardiac steroid glycosides (e.g., digoxin), but not the aglycones, display a higher binding affinity for the α2ß1 and α3ß1 than for the α1ß1 complex, both hellebrin and its aglycone hellebrigenin display ~2-fold higher binding affinity for α1ß1 than for the α2ß1 and α3ß1 complexes. Finally, the current study highlights a common feature for all cardiotonic steroids analyzed here, namely a dramatic reduction in the oxygen consumption rate in cardenolide- and bufadienolide-treated cells, reflecting a direct impact on mitochondrial oxidative phosphorylation. CONCLUSIONS: Altogether, these data show that the binding affinity of the bufadienolides and cardenolides under study is usually higher for the α2ß1 and α3ß1 than for the α1ß1 NaK complex, excepted for hellebrin and its aglycone form, hellebrigenin, with hellebrigenin being as potent as hellebrin in inhibiting in vitro cancer cell growth.

TAF4b and TAF4 differentially regulate mouse embryonic stem cells maintenance and proliferation.

TAF4b is a cell type-specific subunit of the general transcription factor TFIID. Here, we show that TAF4b is highly expressed in embryonic stem cells (ESC) and is down-regulated upon differentiation. To examine the role of TAF4b in ESC, we applied a knockdown (KD) approach. TAF4b depletion is associated with morphological changes and reduced expression of the self-renewal marker alkaline phosphatase. In contrast, KD of TAF4, a ubiquitously expressed TAF4b paralog, retained and even stabilized ESC stemness. Retinoic acid-induced differentiation was facilitated in the absence of TAF4b but was significantly delayed by TAF4 KD. Furthermore, TAF4b supports, whereas TAF4 inhibits, ESC proliferation and cell cycle progression. We identified a subset of TAF4b target genes preferentially expressed in ESC and controlling the cell cycle. Among them are the germ cell-specific transcription factor Sohlh2 and the protein kinase Yes1, which was recently shown to regulate ESC self-renewal. Interestingly, Sohlh2 and Yes1 are also targets of the pluripotency factor Oct4, and their regulation by Oct4 is TAF4b-dependent. Consistent with that, TAF4b but not TAF4 interacts with Oct4. Our findings suggest that TAF4b cooperates with Oct4 to regulate a subset of genes in ESC, whereas TAF4 is required for later embryonic developmental stages.

Differential regulation of NF-kappaB by elongation factors is determined by core promoter type.

NF-kappaB transcription factors activate genes important for immune response, inflammation, and cell survival. P-TEFb and DSIF, which are positive and negative transcription elongation factors, respectively, both regulate NF-kappaB-induced transcription, but the mechanism underlying their recruitment to NF-kappaB target genes is unknown. We show here that upon induction of NF-kappaB, a subset of target genes is regulated differentially by either P-TEFb or DSIF. The regulation of these genes and their occupancy by these elongation factors are dependent on the NF-kappaB enhancer and the core promoter type. Converting a TATA-less promoter to a TATA promoter switches the regulation of NF-kappaB from DSIF to P-TEFb. Accumulation or displacement of DSIF and P-TEFb is dictated by the formation of distinct initiation complexes (TFIID dependent or independent) on the two types of core promoter. The underlying mechanism for the dissociation of DSIF from TATA promoters upon NF-kappaB activation involves the phosphorylation of RNA polymerase II by P-TEFb. The results highlight a regulatory link between the initiation and the elongation phases of the transcription reaction and broaden our comprehension of the NF-kappaB pathway.

Increased myelinating capacity of embryonic stem cell derived oligodendrocyte precursors after treatment by interleukin-6/soluble interleukin-6 receptor fusion protein.

Neurosphere cells (NSc) derived from embryonic stem cells have characteristics of neural stem cells and can differentiate into oligodendrocyte precursors. Culture of NSc with IL6RIL6 chimera (soluble interleukin-6 receptor fused to interleukin-6) enhances their differentiation into oligodendrocytes with longer and more numerous branches and with peripheral accumulation of myelin basic protein (MBP) in myelin membranes indicating maturation. Gene expression profiling reveals that one of the proteins strongly induced by IL6RIL6 is a regulator of microtubule dynamics, stathmin-like 2 (SCG10/Stmn2), and gene silencing shows that Stmn2 plays an important role in the development of the mature oligodendrocyte morphology. IL6RIL6 acts as an effective stimulator of the myelinating function of ES cell-derived oligodendrocyte precursors, as observed upon transplantation of the IL6RIL6- pretreated cells into brain slices of MBP-deficient shiverer mice.

Elongation inhibition by DRB sensitivity-inducing factor is regulated by the A20 promoter via a novel negative element and NF-kappaB.

A20 is an immediate-early NF-kappaB target gene. Prior to NF-kappaB stimulation, the A20 promoter is bound by the polymerase II machinery to allow rapid transcription activation. Here we show that the basal A20 transcription is repressed at the level of elongation in a promoter-specific fashion. Immunodepletion in vitro and RNA interference in cultured cells suggest that the basal elongation inhibition is conferred by DRB sensitivity-inducing factor (DSIF). We have identified a negative upstream promoter element called ELIE that controls DSIF activity. Remarkably, following NF-kappaB stimulation, inhibition of the A20 promoter by DSIF persists, but it is now regulated by NF-kappaB rather than ELIE. Similar regulation by DSIF is shown for another NF-kappaB-responsive gene, the IkappaBalpha gene. These findings reveal an intimate and dynamic relationship between DSIF inhibition of elongation and promoter-bound transcription factors. The potential significance of the differential regulation of DSIF activity by cis-acting elements is discussed.