Large-scale integration of single-cell transcriptomic data captures transitional progenitor states in mouse skeletal muscle regeneration.

Skeletal muscle repair is driven by the coordinated self-renewal and fusion of myogenic stem and progenitor cells. Single-cell gene expression analyses of myogenesis have been hampered by the poor sampling of rare and transient cell states that are critical for muscle repair, and do not inform the spatial context that is important for myogenic differentiation. Here, we demonstrate how large-scale integration of single-cell and spatial transcriptomic data can overcome these limitations. We created a single-cell transcriptomic dataset of mouse skeletal muscle by integration, consensus annotation, and analysis of 23 newly collected scRNAseq datasets and 88 publicly available single-cell (scRNAseq) and single-nucleus (snRNAseq) RNA-sequencing datasets. The resulting dataset includes more than 365,000 cells and spans a wide range of ages, injury, and repair conditions. Together, these data enabled identification of the predominant cell types in skeletal muscle, and resolved cell subtypes, including endothelial subtypes distinguished by vessel-type of origin, fibro-adipogenic progenitors defined by functional roles, and many distinct immune populations. The representation of different experimental conditions and the depth of transcriptome coverage enabled robust profiling of sparsely expressed genes. We built a densely sampled transcriptomic model of myogenesis, from stem cell quiescence to myofiber maturation, and identified rare, transitional states of progenitor commitment and fusion that are poorly represented in individual datasets. We performed spatial RNA sequencing of mouse muscle at three time points after injury and used the integrated dataset as a reference to achieve a high-resolution, local deconvolution of cell subtypes. We also used the integrated dataset to explore ligand-receptor co-expression patterns and identify dynamic cell-cell interactions in muscle injury response. We provide a public web tool to enable interactive exploration and visualization of the data. Our work supports the utility of large-scale integration of single-cell transcriptomic data as a tool for biological discovery.

Direct Determination of Pseudouridine in RNA by Mass Spectrometry Coupled with Stable Isotope Labeling.

Pseudouridine (Ψ) is the only "mass-silent" nucleoside produced by post-transcriptional RNA modification. We developed a mass spectrometry (MS)-based technique coupled with in vivo deuterium (D) labeling of uridines for direct determination of Ψs in cellular RNA and applied it to the comprehensive analysis of post-transcriptional modifications in human ribosomal RNAs. The method utilizes human TK6/mouse FM3A cells deficient in uridine monophosphate synthase using a CRISPR-Cas9 technique to turn off de novo uridine synthesis and fully labels uridines with D at uracil positions 5 and 6 by cultivating the cells in a medium containing uridine-5,6-D2. The pseudouridylation reaction in those cells results in the exchange of the D at the C5 of uracil with hydrogen from solvent, which produces a -1 Da mass shift, thus allowing MS-based determination of RNA Ψs. We present here the experimental details of this method and show that it allows the identification of all Ψs in human major nuclear and nucleolar RNAs, including several previously unknown Ψs. Because the method allows direct determination of Ψs at the femtomole level of RNA, it will serve as a useful tool for structure/function studies of a wide variety of noncoding RNAs.

Evaluating the effect of body fluid mixture on the relative expression ratio of blood-specific RNA markers.

The estimation of the time elapsed since a biological stain was deposited at a crime scene can provide crucial information to a forensic investigation, indicating either when a crime was committed, or whether the biological evidence was deposited at the time of a known crime event. This would enable the investigators to limit the number of suspects and to assess alibis. The relative expression ratios (RERs) of body fluid-specific RNA markers are promising molecular tools for indicating the age of biological stains. However, the nature of some forensic samples found at crime scenes could be challenging, as they frequently occur in a mixture of different body fluid types. The research presented here has utilised reverse transcription quantitative PCR (RT-qPCR) to explore the impact of bloodstains being present in mixtures with other body fluids (saliva or semen) on the resulting RERs of blood-specific markers. The expression level of three blood-specific markers (HBA, HBB and miR16) along with two reference genes (18S and U6) were analysed across multiple ageing time points in pure and mixed bloodstains. For some markers, no significant differences were found when comparing RERs in pure and mixed bloodstains, however some RERs were altered in mixed stains. This indicates that the presence of body fluid mixtures may have a significant effect on the RERs of some blood-specific markers. This should therefore be considered when selecting markers for estimating the age of stains, particularly when multiple body fluids are thought to be present.

Neuronal vulnerability and multilineage diversity in multiple sclerosis.

Multiple sclerosis (MS) is a neuroinflammatory disease with a relapsing-remitting disease course at early stages, distinct lesion characteristics in cortical grey versus subcortical white matter and neurodegeneration at chronic stages. Here we used single-nucleus RNA sequencing to assess changes in expression in multiple cell lineages in MS lesions and validated the results using multiplex in situ hybridization. We found selective vulnerability and loss of excitatory CUX2-expressing projection neurons in upper-cortical layers underlying meningeal inflammation; such MS neuron populations exhibited upregulation of stress pathway genes and long non-coding RNAs. Signatures of stressed oligodendrocytes, reactive astrocytes and activated microglia mapped most strongly to the rim of MS plaques. Notably, single-nucleus RNA sequencing identified phagocytosing microglia and/or macrophages by their ingestion and perinuclear import of myelin transcripts, confirmed by functional mouse and human culture assays. Our findings indicate lineage- and region-specific transcriptomic changes associated with selective cortical neuron damage and glial activation contributing to progression of MS lesions.

Chromatin-Associated RNA Sequencing (ChAR-seq).

RNA is a fundamental component of chromatin. Noncoding RNAs (ncRNAs) can associate with chromatin to influence gene expression and chromatin state; many also act at long distances from their transcriptional origin. Yet we know almost nothing about the functions or sites of action for most ncRNAs. Current methods to identify sites of RNA interaction with the genome are limited to the study of a single RNA at a time. Here we describe a protocol for ChAR-seq, a strategy to identify all chromatin-associated RNAs and map their DNA contacts genome-wide. In ChAR-seq, proximity ligation of RNA and DNA to a linker molecule is used to construct a chimeric RNA-DNA molecule that is converted to DNA for sequencing. In a single assay, ChAR-seq can discover de novo chromatin interactions of distinct RNAs, including nascent transcripts, splicing RNAs, and long noncoding RNAs (lncRNAs). Resulting "maps" of genome-bound RNAs should provide new insights into RNA biology. © 2019 by John Wiley & Sons, Inc.

Prognostic value of small nuclear RNAs (snRNAs) for digestive tract pan- adenocarcinomas identified by RNA sequencing data.

Malignant tumors of the digestive tract include esophageal, gastric, and colorectal carcinomas, which all have high global mortality rates. A clinical role for small nuclear RNA (snRNA), a type of small non-coding RNA, has not yet been documented for digestive tract pan-adenocarcinomas. Therefore, the aim of the study was to identify differentially expressed snRNAs and to explore their prognostic implications in pan-adenocarcinomas from the esophagus, stomach, colon, and rectum. The pan-carcinoma RNA-sequencing data of four types of digestive tract cancers with 1, 102 cases obtained from The Cancer Genome Atlas (TCGA) project were analyzed and the differentially expressed snRNAs were evaluated using the edgeR package. The prognostic value of each of the selected snRNAs was determined by univariate and multivariate Cox regression analyses. All the digestive tract pan-adenocarcinomas showed differential expression of three snRNAs: the up-regulated RNU1-106 P and RNU6-850 P and the down-regulated RNU6-529 P. Interestingly, RNU6-101 P appeared to be a risk factor for esophageal adenocarcinoma (ESAD) and RNVU1-4 was potentially a protective factor for stomach adenocarcinoma (STAD) survival. This consistent finding of differential expression of all three snRNAs in all four types of digestive system cancers suggests potential roles for these snRNAs in the tumorigenesis of digestive system cancers. RNU6-101 P could play a pivotal role in the progression of ESAD and RNVU1-4 could perform a protective role in STAD. However, since the current findings were based on RNA-sequencing data mining, more studies are needed for verification.

U6 can be used as a housekeeping gene for urinary sediment miRNA studies of IgA nephropathy.

Recent studies have indicated that urinary sediment miRNAs not only are able to serve as non-invasive diagnostic biomarkers for IgA nephropathy (IgAN) but may also be closely related to several clinical and pathological indicators. However, the lack of a suitable internal reference miRNA has hampered research into urinary sediment miRNAs. To date, U6 has been used as a reference gene in urinary sediment miRNA studies mostly based on the results from studies using tissue samples and cell lines. In a total of 330 IgAN patients, 164 disease control patients and 130 normal control patients, there was no significant difference in U6 levels. We also compared the U6 levels in different types of primary glomerular disease groups (IgA nephropathy, membranous nephropathy, minimal change nephrosis and focal segmental glomerular sclerosis). The results confirmed that there was no significant difference in the expression of U6 in different primary glomerular disease groups. Moreover, treatment had no significant effect on the expression levels of U6 in IgA nephropathy. Therefore, U6 is an excellent housekeeping gene for urinary sediment miRNA studies of IgA nephropathy.

Label-Free Direct Detection of MiRNAs with Poly-Silicon Nanowire Biosensors.

The diagnostic and prognostic value of microRNAs (miRNAs) in diseases becomes promising. Owing to fast response and high sensitivity, silicon nanowire (SiNW) biosensor has been considered a potential tool for miRNAs detection. Here, we describe a booming method to detect miRNAs with poly-silicon nanowire biosensors. Standard and real miRNA samples are applied in this study. The results show a limitation of 1 fM in the detection of standard miRNA sample with our poly-nanowire devices. Meanwhile, one-base mismatched sequence could be distinguished. Furthermore, these poly-SiNW arrays can detect snRNA U6 in total RNA samples extracted from HepG2 cells with a detection limitation of 0.2 µg/mL.

Differential long-term stability of microRNAs and RNU6B snRNA in 12-20 year old archived formalin-fixed paraffin-embedded specimens.

BACKGROUND: The quantitative analysis of microRNA (miRNA) gene expression in archived formalin-fixed, paraffin embedded (FFPE) tissues has been instrumental to identifying their potential roles in cancer biology, diagnosis, and prognosis. However, it remains unclear whether miRNAs remain stable in FFPE tissues stored for long periods of time. METHODS: Here we report Taqman real-time RT-PCR quantification of miR-21, miR-141, miR-221, and RNU6B small nuclear RNA (snRNA) levels from 92 radical prostatectomy specimens stored for 12-20 years in FFPE blocks. The relative stability of each transcript over time was assessed using general linear models. The correlation between transcript quantities, sample age, and RNA integrity number (RIN) were determined utilizing Spearman rank correlation. RESULTS: All transcript levels linearly decreased with sample age, demonstrating a clear loss of miRNA stability and RNU6B snRNA stability over time. The most rapid rates of degradation were observed for RNU6B and miR-21, while miR-141 and miR-221 were more stable. RNA quality was not correlated with sample age or with miR-21, miR-221, or RNU6B snRNA levels. Conversely, miR-141 levels increased with RNA quality. CONCLUSIONS: MiRNA and snRNA levels gradually decreased over an eight year period in FFPE tissue blocks. Sample age was the most consistent feature associated with miRNA stability. The reference snRNA, RUN6B, was more rapidly degraded when compared to miR-141 and miR-221 miRNAs. Various miRNAs demonstrated differential rates of degradation. Quantitative miRNA studies from long-term archived FFPE tissues may therefore benefit from epidemiologic study design or statistical analysis methods that take into account differential storage-dependent transcript degradation.

[Correlation between RNA Expression Level and Early PMI in Human Brain Tissue].

OBJECTIVES: To explore the correlation between the expression levels of several RNA markers in human brain tissue and early postmortem interval （PMI）. METHODS: Twelve individuals with known PMI （range from 4.3 to 22.5 h） were selected and total RNA was extracted from brain tissue. Eight commonly used RNA markers were chosen including ß-actin, GAPDH, RPS29, 18S rRNA, 5S rRNA, U6 snRNA, miRNA-9 and miRNA-125b, and the expression levels were detected in brain tissue by real-time fluorescent quantitative PCR. The internal reference markers with stable expression in early PMI were screened using geNorm software and the relationship between its expression level and some relevant factors such as age, gender and cause of death were analyzed. RNA markers normalized by internal reference were inserted into the mathematic model established by previous research for PMI estimation using R software. Model quality was judged by the error rate calculated with estimated PMI. RESULTS: 5S rRNA, miRNA-9 and miRNA-125b showed quite stable expression and their expression levels had no relation with age, gender and cause of death. The error rate of estimated PMI using ß-actin was 24.6%, while GAPDH was 41.0%. CONCLUSIONS: 5S rRNA, miRNA-9 and miRNA-125b are suitable as internal reference markers of human brain tissue owing to their stable expression in early PMI. The expression level of ß-actin correlates well with PMI, which can be used as an additional index for early PMI estimation.

Profiling of small RNA cargo of extracellular vesicles shed by Trypanosoma cruzi reveals a specific extracellular signature.

Over the last years, an expanding family of small regulatory RNAs (e.g. microRNAs, siRNAs and piRNAs) was recognized as key players in novel forms of post-transcriptional gene regulation in most eukaryotes. However, the machinery associated with Ago/Dicer-dependent small RNA biogenesis was thought to be either entirely lost or extensively simplified in some unicellular organisms including Trypanosoma cruzi, Saccharomyces cerevisiae, Leishmania major and Plasmodium falciparum. Although the biogenesis of small RNAs from non-coding RNAs represent a minor fraction of the normal small RNA transcriptome in eukaryotic cells, they represent the unique small RNA pathways in Trypanosoma cruzi which produce different populations of small RNAs derived from tRNAs, rRNAs, sn/snoRNAs and mRNAs. These small RNAs are secreted included in extracellular vesicles and transferred to other parasites and susceptible mammalian cells. This process represents a novel form of cross-kingdom transfer of genetic material suggesting that secreted vesicles could represent new relevant pieces in life cycle transitions, infectivity and cell-to-cell communication. Here, we provide for the first time a detailed analysis of the small RNA cargo of extracellular vesicles from T. cruzi epimastigotes under nutritional stress conditions compared to the respective intracellular compartment using deep sequencing. Compared with the intracellular compartment, shed extracellular vesicles showed a specific extracellular signature conformed by distinctive patterns of small RNAs derived from rRNA, tRNA, sno/snRNAs and protein coding sequences which evidenced specific secretory small RNA processing pathways.

Impact of Hfq on the Bacillus subtilis transcriptome.

The RNA chaperone Hfq acts as a central player in post-transcriptional gene regulation in several Gram-negative Bacteria, whereas comparatively little is known about its role in Gram-positive Bacteria. Here, we studied the function of Hfq in Bacillus subtilis, and show that it confers a survival advantage. A comparative transcriptome analysis revealed mRNAs with a differential abundance that are governed by the ResD-ResE system required for aerobic and anaerobic respiration. Expression of resD was found to be up-regulated in the hfq- strain. Furthermore, several genes of the GerE and ComK regulons were de-regulated in the hfq- background. Surprisingly, only six out of >100 known and predicted small RNAs (sRNAs) showed altered abundance in the absence of Hfq. Moreover, Hfq positively affected the transcript abundance of genes encoding type I toxin-antitoxin systems. Taken the moderate effect on sRNA levels and mRNAs together, it seems rather unlikely that Hfq plays a central role in RNA transactions in Bacillus subtilis.

Evaluation of suitable reference genes for normalization of microRNA expression by real-time reverse transcription PCR analysis during longan somatic embryogenesis.

Accurate profiling of microRNAs (miRNAs) is an essential step for understanding both developmental and physiological functions of miRNAs. Real-time quantitative PCR (qPCR) is being widely used in miRNA expression studies, but choosing a suitable reference gene is a crucial factor for correct analysis of results. To date, there has been no systematic evaluation of qPCR reference genes for the study of miRNAs during somatic embryogenesis (SE) in the longan tree (Dimocarpus longan). Here, the most stably expressed miRNAs in synchronized longan tree embryogenic cultures at different developmental stages were determined using the geNorm and NormFinder algorithms. Validation qPCR experiments were performed for 24 miRNAs together with a snRNA (U6 snRNA), a rRNA (5S rRNA), and three housekeeping genes. It was found that small RNAs had better expression stability than protein-coding genes, and dlo-miR24 was identified as the most reliable reference gene, followed by dlo-miR168a*, dlo-miR2089*-1 and 5S rRNA. dlo-miR24 was recommended as a normalizer if only a single reference gene was to be used, while the combination of dlo-miR156c, dlo-2089*-1 and 5S rRNA was preferred to normalize miRNA expression data during longan SE.

New insights into the spliceosome by single molecule fluorescence microscopy.

Splicing is an essential eukaryotic process in which introns are excised from precursors to messenger RNAs and exons ligated together. This reaction is catalyzed by a multi-MegaDalton machine called the spliceosome, composed of 5 small nuclear RNAs (snRNAs) and a core set of ∼100 proteins minimally required for activity. Because of the spliceosome's size, its low abundance in cellular extracts, and its highly dynamic assembly pathway, analysis of the kinetics of splicing and the conformational rearrangements occurring during spliceosome assembly and disassembly has proven extraordinarily challenging. Here, we review recent progress in combining chemical biology methodologies with single molecule fluorescence techniques to provide a window into splicing in real time. These methods complement ensemble measurements of splicing in vivo and in vitro to facilitate kinetic dissection of pre-mRNA splicing.

snoU6 and 5S RNAs are not reliable miRNA reference genes in neuronal differentiation.

Accurate profiling of microRNAs (miRNAs) is an essential step for understanding the functional significance of these small RNAs in both physiological and pathological processes. Quantitative real-time PCR (qPCR) has gained acceptance as a robust and reliable transcriptomic method to profile subtle changes in miRNA levels and requires reference genes for accurate normalization of gene expression. 5S and snoU6 RNAs are commonly used as reference genes in microRNA quantification. It is currently unknown if these small RNAs are stably expressed during neuronal differentiation. Panels of miRNAs have been suggested as alternative reference genes to 5S and snoU6 in various physiological contexts. To test the hypothesis that miRNAs may serve as stable references during neuronal differentiation, the expressions of eight miRNAs, 5S and snoU6 RNAs in five differentiating neuronal cell types were analyzed using qPCR. The stabilities of the expressions were evaluated using two complementary statistical approaches (geNorm and Normfinder). Expressions of 5S and snoU6 RNAs were stable under some but not all conditions of neuronal differentiation and thus are not suitable reference genes. In contrast, a combination of three miRNAs (miR-103, miR-106b and miR-26b) allowed accurate expression normalization across different models of neuronal differentiation.

A homozygous mutation in RNU4ATAC as a cause of microcephalic osteodysplastic primordial dwarfism type I (MOPD I) with associated pigmentary disorder.

The designation microcephalic osteodysplastic primordial dwarfism (MOPD) refers to a group of autosomal recessive disorders, comprising microcephaly, growth retardation, and a skeletal dysplasia. The different types of MOPD have been delineated on the basis of clinical, radiological, and genetic criteria. We describe two brothers, born to healthy, consanguineous parents, with intrauterine and postnatal growth retardation, microcephaly with abnormal gyral pattern and partial agenesis of corpus callosum, and skeletal anomalies reminiscent of those described in MOPD type I. This was confirmed by the identification of the homozygous g.55G > A mutation of RNU4ATAC encoding U4atac snRNA. The sibs had yellowish-gray hair, fair skin, and deficient retinal pigmentation. Skin biopsy showed abnormal melanin function but OCA genes were normal. The older sib had an intracranial hemorrhage at 1 week after birth, the younger developed chilblains-like lesions at the age 2½ years old but analysis of the SAMHD1 and TREX1 genes did not show any mutations. To the best of our knowledge, vasculopathy and pigmentary disorders have not been reported in MOPD I.

snRNA-specific role of SMN in trypanosome snRNP biogenesis in vivo.

Pre-mRNA splicing in trypanosomes requires the SMN-mediated assembly of small nuclear ribonucleoproteins (snRNPs). In contrast to higher eukaryotes, the cellular localization of snRNP biogenesis and the involvement of nuclear-cytoplasmic trafficking in trypanosomes are controversial. By using RNAi knockdown of SMN in T. brucei to investigate its functional role in snRNP assembly, we found dramatic changes in the steady-state levels of snRNAs and snRNPs: The SL RNA accumulates, whereas U1, U4, and U5 snRNA levels decrease, and Sm core assembly in particular of the SL RNA is strongly reduced. In addition, SMN depletion blocks U4/U6 di-snRNP formation; the variant Sm core of the U2 snRNP, however, still forms efficiently after SMN knockdown. Concerning the longstanding question, whether nuclear-cytoplasmic trafficking is involved in trypanosomal snRNP biogenesis, fluorescence in situ hybridization (FISH) and immunofluorescence assays revealed that the SL RNA genes and transcripts colocalize with SMN. Remarkably, SMN silencing leads to a nucleoplasmic accumulation of both SL RNA and the Sm proteins. In sum, our data demonstrate an essential and snRNA-selective role of SMN in snRNP biogenesis in vivo and strongly argue for a nucleoplasmic Sm core assembly of the SL RNP.

PsRNA: a computing engine for the comparative identification of putative small RNA locations within intergenic regions.

Small RNAs (sRNAs) are non-coding transcripts exerting their functions in the cells directly. Identification of sRNAs is a difficult task due to the lack of clear sequence and structural biases. Most sRNAs are identified within genus specific intergenic regions in related genomes. However, several of these regions remain un-annotated due to lack of sequence homology and/or potent statistical identification tools. A computational engine has been built to search within the intergenic regions to identify and roughly annotate new putative sRNA regions in Enterobacteriaceae genomes. It utilizes experimentally known sRNA data and their flanking genes/KEGG Orthology (KO) numbers as templates to identify similar sRNA regions in related query genomes. The search engine not only has the capability to locate putative intergenic regions for specific sRNAs, but also has the potency to locate conserved, shuffled or deleted gene clusters in query genomes. Because it uses the KO terms for locating functionally important regions such as sRNAs, any further KO number assignment to additional genes will increase the sensitivity. The PsRNA server is used for the identification of putative sRNA regions through the information retrieved from the sRNA of interest. The computing engine is available online at http://bioserver1.physics.iisc.ernet.in/psrna/ and http://bicmku.in:8081/psrna/.

In-gel digestion for mass spectrometric characterization of RNA from fluorescently stained polyacrylamide gels.

Although current mass spectrometry-based proteomics technology allows for high-throughput analysis of protein components in functional ribonucleoprotein complexes, this technology has had limited application to studies of RNA itself. Here we present a protocol for RNA analysis using polyacrylamide gel electrophoresis coupled with liquid chromatography-tandem mass spectrometry. Specifically, RNAs of interest are subjected to polyacrylamide gel electrophoresis and stained with a fluorescent dye, and RNAs in gel bands are digested with nuclease and then analyzed directly liquid chromatography-mass spectrometry, resulting in highly accurate mass values and reliable information on post-transcriptional modifications. We demonstrate that the method can be applied to the identification and chemical analysis of small RNAs in mouse embryonic stem cell extracts and of small RNAs in the spliceosomal ribonucleoprotein complex pulled down from yeast cells using a tagged protein cofactor as bait. The protocol is relatively simple and allowed us to identify not only three novel methylated nucleotide residues of RNase P RNA, U6 snRNA, and 7SL RNA prepared from mouse ES cells but also various 3'-end forms of U4, U5S, and U6 snRNAs isolated from the yeast spliceosome at the femtomole level. The method is thus a convenient tool for direct analysis of RNAs in various cellular ribonucleoprotein complexes, particularly for the analysis of post-transcriptional modifications and metabolic processing of RNA.