Quantification of escape from X chromosome inactivation with single-cell omics data reveals heterogeneity across cell types and tissues.

Several X-linked genes escape from X chromosome inactivation (XCI), while differences in escape across cell types and tissues are still poorly characterized. Here, we developed scLinaX for directly quantifying relative gene expression from the inactivated X chromosome with droplet-based single-cell RNA sequencing (scRNA-seq) data. The scLinaX and differentially expressed gene analyses with large-scale blood scRNA-seq datasets consistently identified the stronger escape in lymphocytes than in myeloid cells. An extension of scLinaX to a 10x multiome dataset (scLinaX-multi) suggested a stronger escape in lymphocytes than in myeloid cells at the chromatin-accessibility level. The scLinaX analysis of human multiple-organ scRNA-seq datasets also identified the relatively strong degree of escape from XCI in lymphoid tissues and lymphocytes. Finally, effect size comparisons of genome-wide association studies between sexes suggested the underlying impact of escape on the genotype-phenotype association. Overall, scLinaX and the quantified escape catalog identified the heterogeneity of escape across cell types and tissues.

Transcriptome regulation by PARP13 in basal and antiviral states in human cells.

The RNA-binding protein PARP13 is a primary factor in the innate antiviral response, which suppresses translation and drives decay of bound viral and host RNA. PARP13 interacts with many proteins encoded by interferon-stimulated genes (ISG) to activate antiviral pathways including co-translational addition of ISG15, or ISGylation. We performed enhanced crosslinking immunoprecipitation (eCLIP) and RNA-seq in human cells to investigate PARP13's role in transcriptome regulation for both basal and antiviral states. We find that the antiviral response shifts PARP13 target localization, but not its binding preferences, and that PARP13 supports the expression of ISGylation-related genes, including PARP13's cofactor, TRIM25. PARP13 associates with TRIM25 via RNA-protein interactions, and we elucidate a transcriptome-wide periodicity of PARP13 binding around TRIM25. Taken together, our study implicates PARP13 in creating and maintaining a cellular environment poised for an antiviral response through limiting PARP13 translation, regulating access to distinct mRNA pools, and elevating ISGylation machinery expression.

Microenvironmental modulation in tandem with human stem cell transplantation enhances functional recovery after chronic complete spinal cord injury.

While rapid advancements in regenerative medicine strategies for spinal cord injury (SCI) have been made, most research in this field has focused on the early stages of incomplete injury. However, the majority of patients experience chronic severe injury; therefore, treatments for these situations are fundamentally important. Here, we hypothesized that environmental modulation via a clinically relevant hepatocyte growth factor (HGF)-releasing scaffold and human iPS cell-derived neural stem/progenitor cells (hNS/PCs) transplantation contributes to functional recovery after chronic complete transection SCI. Effective release of HGF from a collagen scaffold induced progressive axonal elongation and increased grafted cell viability by activating microglia/macrophages and meningeal cells, inhibiting inflammation, reducing scar formation, and enhancing vascularization. Furthermore, hNS/PCs transplantation enhanced endogenous neuronal regrowth, the extension of graft axons, and the formation of circuits around the lesion and lumbar enlargement between host and graft neurons, resulting in the restoration of locomotor and urinary function. This study presents an effective therapeutic strategy for severe chronic SCI and provides evidence for the feasibility of regenerative medicine strategies using clinically relevant materials.

Single-cell transcriptomics of human gut T cells identifies cytotoxic CD4+CD8A+ T cells related to mouse CD4 cytotoxic T cells.

Cytotoxic CD4+ T cells (CD4-CTLs) show the presence of cytolytic granules, which include the enzymes granzyme and perforin. The cells have a pathogenic and protective role in various diseases, including cancer, viral infection, and autoimmune disease. In mice, cytotoxic CD4+ T cells express CD8αα+ and reside in the intestine (mouse CD4+CTLs; mCD4-CTLs). The population of cytotoxic CD4+ T cells in the human intestine is currently unknown. Moreover, it is unclear how cytotoxic CD4 T cells change in patients with inflammatory bowel disease (IBD). Here, we aimed to identify cytotoxic CD4+ T cells in the human intestine and analyze the characteristics of the population in patients with IBD using single-cell RNA-seq (scRNA-seq). In CD4+ T cells, granzyme and perforin expression was high in humanMAIT (hMAIT) cells and hCD4+CD8A+ T cell cluster. Both CD4 and CD8A were expressed in hTreg, hMAIT, and hCD4+CD8A+ T cell clusters. Next we performed fast gene set enrichment analysis to identify cell populations that showed homology to mCD4CTLs. The analysis identified the hCD4+CD8A+ T cell cluster (hCTL-like population; hCD4-CTL) similar to mouse CTLs. The percentage of CD4+CD8A+ T cells among the total CD4+ T cells in the inflamed intestine of the patients with Crohn's disease was significantly reduced compared with that in the noninflamed intestine of the patients. In summary, we identified cytotoxic CD4+CD8+ T cells in the small intestine of humans. The integration of the mouse and human sc-RNA-seq data analysis highlight an approach to identify human cell populations related to mouse cell populations, which may help determine the functional properties of several human cell populations in mice.

SCAFE: a software suite for analysis of transcribed cis-regulatory elements in single cells.

MOTIVATION: Cell type-specific activities of cis-regulatory elements (CRE) are central to understanding gene regulation and disease predisposition. Single-cell RNA 5'end sequencing (sc-end5-seq) captures the transcription start sites (TSS) which can be used as a proxy to measure the activity of transcribed CREs (tCREs). However, a substantial fraction of TSS identified from sc-end5-seq data may not be genuine due to various artifacts, hindering the use of sc-end5-seq for de novo discovery of tCREs. RESULTS: We developed SCAFE-Single-Cell Analysis of Five-prime Ends-a software suite that processes sc-end5-seq data to de novo identify TSS clusters based on multiple logistic regression. It annotates tCREs based on the identified TSS clusters and generates a tCRE-by-cell count matrix for downstream analyses. The software suite consists of a set of flexible tools that could either be run independently or as pre-configured workflows. AVAILABILITY AND IMPLEMENTATION: SCAFE is implemented in Perl and R. The source code and documentation are freely available for download under the MIT License from https://github.com/chung-lab/SCAFE. Docker images are available from https://hub.docker.com/r/cchon/scafe. The submitted software version and test data are archived at https://doi.org/10.5281/zenodo.7023163 and https://doi.org/10.5281/zenodo.7024060, respectively. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Deep sequencing of short capped RNAs reveals novel families of noncoding RNAs.

In eukaryotes, capped RNAs include long transcripts such as messenger RNAs and long noncoding RNAs, as well as shorter transcripts such as spliceosomal RNAs, small nucleolar RNAs, and enhancer RNAs. Long capped transcripts can be profiled using cap analysis gene expression (CAGE) sequencing and other methods. Here, we describe a sequencing library preparation protocol for short capped RNAs, apply it to a differentiation time course of the human cell line THP-1, and systematically compare the landscape of short capped RNAs to that of long capped RNAs. Transcription initiation peaks associated with genes in the sense direction have a strong preference to produce either long or short capped RNAs, with one out of six peaks detected in the short capped RNA libraries only. Gene-associated short capped RNAs have highly specific 3' ends, typically overlapping splice sites. Enhancers also preferentially generate either short or long capped RNAs, with 10% of enhancers observed in the short capped RNA libraries only. Enhancers producing either short or long capped RNAs show enrichment for GWAS-associated disease SNPs. We conclude that deep sequencing of short capped RNAs reveals new families of noncoding RNAs and elucidates the diversity of transcripts generated at known and novel promoters and enhancers.

An era of single-cell genomics consortia.

The human body consists of 37 trillion single cells represented by over 50 organs that are stitched together to make us who we are, yet we still have very little understanding about the basic units of our body: what cell types and states make up our organs both compositionally and spatially. Previous efforts to profile a wide range of human cell types have been attempted by the FANTOM and GTEx consortia. Now, with the advancement in genomic technologies, profiling the human body at single-cell resolution is possible and will generate an unprecedented wealth of data that will accelerate basic and clinical research with tangible applications to future medicine. To date, several major organs have been profiled, but the challenges lie in ways to integrate single-cell genomics data in a meaningful way. In recent years, several consortia have begun to introduce harmonization and equity in data collection and analysis. Herein, we introduce existing and nascent single-cell genomics consortia, and present benefits to necessitate single-cell genomic consortia in a regional environment to achieve the universal human cell reference dataset.

CDK1 dependent phosphorylation of hTERT contributes to cancer progression.

The telomerase reverse transcriptase is upregulated in the majority of human cancers and contributes directly to cell transformation. Here we report that hTERT is phosphorylated at threonine 249 during mitosis by the serine/threonine kinase CDK1. Clinicopathological analyses reveal that phosphorylation of hTERT at threonine 249 occurs more frequently in aggressive cancers. Using CRISPR/Cas9 genome editing, we introduce substitution mutations at threonine 249 in the endogenous hTERT locus and find that phosphorylation of threonine 249 is necessary for hTERT-mediated RNA dependent RNA polymerase (RdRP) activity but dispensable for reverse transcriptase and terminal transferase activities. Cap Analysis of Gene Expression (CAGE) demonstrates that hTERT phosphorylation at 249 regulates the expression of specific genes that are necessary for cancer cell proliferation and tumor formation. These observations indicate that phosphorylation at threonine 249 regulates hTERT RdRP and contributes to cancer progression in a telomere independent manner.

ELTA: Enzymatic Labeling of Terminal ADP-Ribose.

ADP-ribosylation refers to the addition of one or more ADP-ribose groups onto proteins. The attached ADP-ribose monomers or polymers, commonly known as poly(ADP-ribose) (PAR), modulate the activities of the modified substrates or their binding affinities to other proteins. However, progress in this area is hindered by a lack of tools to investigate this protein modification. Here, we describe a new method named ELTA (enzymatic labeling of terminal ADP-ribose) for labeling free or protein-conjugated ADP-ribose monomers and polymers at their 2'-OH termini using the enzyme OAS1 and dATP. When coupled with various dATP analogs (e.g., radioactive, fluorescent, affinity tags), ELTA can be used to explore PAR biology with techniques routinely used to investigate DNA or RNA function. We demonstrate that ELTA enables the biophysical measurements of protein binding to PAR of a defined length, detection of PAR length from proteins and cells, and enrichment of sub-femtomole amounts of ADP-ribosylated peptides from cell lysates.

Identifying and characterizing functional 3' nucleotide addition in the miRNA pathway.

Over the past decade, modifications to microRNAs (miRNAs) via 3' end nucleotide addition have gone from a deep-sequencing curiosity to experimentally confirmed drivers of a range of regulatory activities. Here we overview the methods that have been deployed by researchers seeking to untangle these diverse functional roles and include characterizing not only the nucleotidyl transferases catalyzing the additions but also the nucleotides being added, and the timing of their addition during the miRNA pathway. These methods and their further development are key to clarifying the diverse and sometimes contradictory functional findings presently attributed to these nucleotide additions.

Semi-quantitative Detection of RNA-dependent RNA Polymerase Activity of Human Telomerase Reverse Transcriptase Protein.

Human telomerase reverse transcriptase (TERT) is the catalytic subunit of telomerase, and it elongates telomere through RNA-dependent DNA polymerase activity. Although TERT is named as a reverse transcriptase, structural and phylogenetic analyses of TERT demonstrate that TERT is a member of right-handed polymerases, and relates to viral RNA-dependent RNA polymerases (RdRPs) as well as viral reverse transcriptase. We firstly identified RdRP activity of human TERT that generates complementary RNA stand to a template non-coding RNA and contributes to RNA silencing in cancer cells. To analyze this non-canonical enzymatic activity, we developed RdRP assay with recombinant TERT in 2009, thereafter established in vitro RdRP assay for endogenous TERT. In this manuscript, we describe the latter method. Briefly, TERT immune complexes are isolated from cells, and incubated with template RNA and rNTPs including radioactive rNTP for RdRP reaction. To eliminate single-stranded RNA, reaction products are treated with RNase I, and the final products are analyzed with polyacrylamide gel electrophoresis. Radiolabeled RdRP products can be detected by autoradiography after overnight exposure.

Methods for studying microRNA functions during stress.

Constituting 5 % of the human genome, microRNAs represent a sizeable class of gene regulators that is predicted to control the expression of at least 60 % of all protein-coding RNAs. Dysregulation of microRNA function results in developmental defects and pathological diseases such as cancers and neurological disorders. Intriguingly, many phenotypes of microRNA deficiencies are subdued in normal condition but manifested apparently upon stress. Here, we outline experimental methods to monitor the level, targets, and activity of microRNAs as the first few steps to characterize how microRNA functions are altered upon stress.

Telomerase reverse transcriptase regulates microRNAs.

MicroRNAs are small non-coding RNAs that inhibit the translation of target mRNAs. In humans, most microRNAs are transcribed by RNA polymerase II as long primary transcripts and processed by sequential cleavage of the two RNase III enzymes, DROSHA and DICER, into precursor and mature microRNAs, respectively. Although the fundamental functions of microRNAs in RNA silencing have been gradually uncovered, less is known about the regulatory mechanisms of microRNA expression. Here, we report that telomerase reverse transcriptase (TERT) extensively affects the expression levels of mature microRNAs. Deep sequencing-based screens of short RNA populations revealed that the suppression of TERT resulted in the downregulation of microRNAs expressed in THP-1 cells and HeLa cells. Primary and precursor microRNA levels were also reduced under the suppression of TERT. Similar results were obtained with the suppression of either BRG1 (also called SMARCA4) or nucleostemin, which are proteins interacting with TERT and functioning beyond telomeres. These results suggest that TERT regulates microRNAs at the very early phases in their biogenesis, presumably through non-telomerase mechanism(s).

PAPD5-mediated 3' adenylation and subsequent degradation of miR-21 is disrupted in proliferative disease.

Next-generation sequencing experiments have shown that microRNAs (miRNAs) are expressed in many different isoforms (isomiRs), whose biological relevance is often unclear. We found that mature miR-21, the most widely researched miRNA because of its importance in human disease, is produced in two prevalent isomiR forms that differ by 1 nt at their 3' end, and moreover that the 3' end of miR-21 is posttranscriptionally adenylated by the noncanonical poly(A) polymerase PAPD5. PAPD5 knockdown caused an increase in the miR-21 expression level, suggesting that PAPD5-mediated adenylation of miR-21 leads to its degradation. Exoribonuclease knockdown experiments followed by small-RNA sequencing suggested that PARN degrades miR-21 in the 3'-to-5' direction. In accordance with this model, microarray expression profiling demonstrated that PAPD5 knockdown results in a down-regulation of miR-21 target mRNAs. We found that disruption of the miR-21 adenylation and degradation pathway is a general feature in tumors across a wide range of tissues, as evidenced by data from The Cancer Genome Atlas, as well as in the noncancerous proliferative disease psoriasis. We conclude that PAPD5 and PARN mediate degradation of oncogenic miRNA miR-21 through a tailing and trimming process, and that this pathway is disrupted in cancer and other proliferative diseases.

New perspectives on the diversification of the RNA interference system: insights from comparative genomics and small RNA sequencing.

Our understanding of the pervasive involvement of small RNAs in regulating diverse biological processes has been greatly augmented by recent application of deep-sequencing technologies to small RNA across diverse eukaryotes. We review the currently known small RNA classes and place them in context of the reconstructed evolutionary history of the RNA interference (RNAi) protein machinery. This synthesis indicates that the earliest versions of eukaryotic RNAi systems likely utilized small RNA processed from three types of precursors: (1) sense-antisense transcriptional products, (2) genome-encoded, imperfectly complementary hairpin sequences, and (3) larger noncoding RNA precursor sequences. Structural dissection of PIWI proteins along with recent discovery of novel families (including Med13 of the Mediator complex) suggest that emergence of a distinct architecture with the N-terminal domains (also occurring separately fused to endoDNases in prokaryotes) formed via duplication of an ancestral unit was key to their recruitment as primary RNAi effectors and use of small RNAs of certain preferred lengths. Prokaryotic PIWI proteins are typically components of several RNA-directed DNA restriction or CRISPR/Cas systems. However, eukaryotic versions appear to have emerged from a subset that evolved RNA-directed RNAi. They were recruited alongside RNaseIII domains and RNA-dependent RNA polymerase (RdRP) domains, also from prokaryotic systems, to form the core eukaryotic RNAi system. Like certain regulatory systems, RNAi diversified into two distinct but linked arms concomitant with eukaryotic nucleocytoplasmic compartmentalization. Subsequent elaboration of RNAi proceeded via diversification of the core protein machinery through lineage-specific expansions and recruitment of new components from prokaryotes (nucleases and small RNA-modifying enzymes), allowing for diversification of associating small RNAs.

Does an emergency visit to the ER make microRNAs stronger during stress?

In this issue of Molecular Cell, Wu et al. (2013) report the identification of GW182-independent microRNA complexes that confer stronger repression upon serum starvation; interestingly, these complexes are associated with polyribosomes and the endoplasmic reticulum.

CpG site-specific alteration of hydroxymethylcytosine to methylcytosine beyond DNA replication.

Hydroxymethylcytosines (hmC), one of several reported cytosine modifications, was recently found to be enriched in embryonic stem cells and neuronal cells, and thought to play an important role in regulating gene expression and cell specification. However, unlike methylcytosines (mC), the fate of hmC beyond DNA replication is not well understood. Here, to monitor the status of hmC during DNA replication, we prepared a stable episomal vector-based monitoring system called MoCEV in 293T cells. The MoCEV system containing fully hydroxymethylated-cytosine fragments revealed a significant modification towards mC after several rounds of DNA replication. Strikingly this modification was specifically observed at the CpG sites (71.9% of cytosines), whereas only 1.1% of modified cytosines were detected at the non-CpG sites. Since the unmodified MoCEV did not undergo any DNA methylation during cell division, the results strongly suggest that somatic cells undergo hmC to mC specifically at the CpG sites during cell division.

Poly(ADP-ribose) regulates post-transcriptional gene regulation in the cytoplasm.

Since its discovery in 1963, poly(ADP-ribose) (pADPr) has been shown to play important functions in the nucleus of multicellular eukaryotes. Each of these functions centers upon DNA metabolism, including DNA-damage repair, chromatin remodeling, transcription and telomere functions. We recently described two novel functions for pADPr in the cytoplasm, both of which involve RNA metabolism - 1) the assembly of cytoplasmic stress granules, cellular macrostructures that aggregate translationally stalled mRNA/protein complexes, and 2) modulation of microRNA activities. Multiple stress granule-localized, post-transcriptional gene regulators, including microRNA-binding argonaute family members, are substrates for pADPr modification and are increasingly modified by pADPr upon stress. Interestingly, the cytoplasmic RNA regulatory functions for PARPs are likely mediated through activities of catalytically inactive PARP-13/ARTD13/ZC3HAV1/ZAP and mono/poly(ADP-ribose)-synthesizing enzymes, including PARP-5a/ARTD5/TNKS1, PARP-12/ARTD12/ZC3HDC1 and PARP-15/ARTD7/BAL3. These data are consistent with other recent work, which suggests that mono(ADP-ribosyl)ated residues can be poly(ADP-ribosyl)ated by different enzymes.

pre-miRNA profiles obtained through application of locked nucleic acids and deep sequencing reveals complex 5'/3' arm variation including concomitant cleavage and polyuridylation patterns.

Recent research hints at an underappreciated complexity in pre-miRNA processing and regulation. Global profiling of pre-miRNA and its potential to increase understanding of the pre-miRNA landscape is impeded by overlap with highly expressed classes of other non coding (nc) RNA. Here, we present a data set excluding these RNA before sequencing through locked nucleic acids (LNA), greatly increasing pre-miRNA sequence counts with no discernable effect on pre-miRNA or mature miRNA sequencing. Analysis of profiles generated in total, nuclear and cytoplasmic cell fractions reveals that pre-miRNAs are subject to a wide range of regulatory processes involving loci-specific 3'- and 5'-end variation entailing complex cleavage patterns with co-occurring polyuridylation. Additionally, examination of nuclear-enriched flanking sequences of pre-miRNA, particularly those derived from polycistronic miRNA transcripts, provides insight into miRNA and miRNA-offset (moRNA) production, specifically identifying novel classes of RNA potentially functioning as moRNA precursors. Our findings point to particularly intricate regulation of the let-7 family in many ways reminiscent of DICER1-independent, pre-mir-451-like processing, introduce novel and unify known forms of pre-miRNA regulation and processing, and shed new light on overlooked products of miRNA processing pathways.

Nuclear pore complex protein mediated nuclear localization of dicer protein in human cells.

Human DICER1 protein cleaves double-stranded RNA into small sizes, a crucial step in production of single-stranded RNAs which are mediating factors of cytoplasmic RNA interference. Here, we clearly demonstrate that human DICER1 protein localizes not only to the cytoplasm but also to the nucleoplasm. We also find that human DICER1 protein associates with the NUP153 protein, one component of the nuclear pore complex. This association is detected predominantly in the cytoplasm but is also clearly distinguishable at the nuclear periphery. Additional characterization of the NUP153-DICER1 association suggests NUP153 plays a crucial role in the nuclear localization of the DICER1 protein.