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.

Multiplexed Analysis of Retinal Gene Expression and Chromatin Accessibility using scRNA-Seq and scATAC-Seq.

Powerful next generation sequencing techniques offer robust and comprehensive analysis to investigate how retinal gene regulatory networks function during development and in disease states. Single-cell RNA sequencing allows us to comprehensively profile gene expression changes observed in retinal development and disease at a cellular level, while single-cell ATAC-Seq allows analysis of chromatin accessibility and transcription factor binding to be profiled at similar resolution. Here the use of these techniques in the developing retina is described, and MULTI-Seq is demonstrated, where individual samples are labeled with a modified oligonucleotide-lipid complex, enabling researchers to both increase the scope of individual experiments and substantially reduce costs.

Preprocessing choices affect RNA velocity results for droplet scRNA-seq data.

Experimental single-cell approaches are becoming widely used for many purposes, including investigation of the dynamic behaviour of developing biological systems. Consequently, a large number of computational methods for extracting dynamic information from such data have been developed. One example is RNA velocity analysis, in which spliced and unspliced RNA abundances are jointly modeled in order to infer a 'direction of change' and thereby a future state for each cell in the gene expression space. Naturally, the accuracy and interpretability of the inferred RNA velocities depend crucially on the correctness of the estimated abundances. Here, we systematically compare five widely used quantification tools, in total yielding thirteen different quantification approaches, in terms of their estimates of spliced and unspliced RNA abundances in five experimental droplet scRNA-seq data sets. We show that there are substantial differences between the quantifications obtained from different tools, and identify typical genes for which such discrepancies are observed. We further show that these abundance differences propagate to the downstream analysis, and can have a large effect on estimated velocities as well as the biological interpretation. Our results highlight that abundance quantification is a crucial aspect of the RNA velocity analysis workflow, and that both the definition of the genomic features of interest and the quantification algorithm itself require careful consideration.

Exploration of small RNA biomarkers for testicular injury in the serum exosomes of rats.

Testicular injury is often observed in drug development. Serum hormones are usually used as noninvasive biomarkers for testicular injury; however, their sensitivities are low. Therefore, it is difficult to monitor testicular injury in drug development. In recent years, molecules in body fluid exosomes have attracted attention as biomarkers for diseases. In this study, small RNAs in serum exosomes were analyzed to identify noninvasive biomarkers of testicular injury in rats, which are often used in preclinical drug development. The rat models of testicular injury were prepared by a single oral administration of 2000 mg/kg ethylene glycol monomethyl ether, in which spermatocyte degeneration and Sertoli cell vacuolation were observed, or 400 mg/kg carbendazim, in which Sertoli cell vacuolation and seminiferous tubule dilation were observed. Serum exosomal small RNA-seq analysis of these models was performed. The analysis identified 3 small RNAs that fluctuated in common between the models, and miR-423-5p and miR-128-3p were selected as candidate markers. For evaluating these candidate markers in other testicular injury models, the models were prepared by a single oral administration of 60 mg/kg 1,3-dinitrobenzene or 500 mg/kg nitrofurazone, and spermatocyte degeneration and Sertoli cell vacuolation were observed. In qPCR analysis, these exosomal miRNAs were upregulated in all models except for the 1,3-dinitrobenzene model, in which severe hemolysis was observed. By contrast, these miRNAs in whole serum extracts did not significantly change in any of the models. In conclusion, we identified miR-423-5p and miR-128-3p in serum exosomes as noninvasive biomarkers for testicular injury in rats.

Single-Cell RNA Profiling of Glomerular Cells Shows Dynamic Changes in Experimental Diabetic Kidney Disease.

BACKGROUND: Recent single-cell RNA sequencing (scRNA-seq) analyses have offered much insight into cell-specific gene expression profiles in normal kidneys. However, in diseased kidneys, understanding of changes in specific cells, particularly glomerular cells, remains limited. METHODS: To elucidate the glomerular cell-specific gene expression changes in diabetic kidney disease, we performed scRNA-seq analysis of isolated glomerular cells from streptozotocin-induced diabetic endothelial nitric oxide synthase (eNOS)-deficient (eNOS-/-) mice and control eNOS-/- mice. RESULTS: We identified five distinct cell populations, including glomerular endothelial cells, mesangial cells, podocytes, immune cells, and tubular cells. Using scRNA-seq analysis, we confirmed the expression of glomerular cell-specific markers and also identified several new potential markers of glomerular cells. The number of immune cells was significantly higher in diabetic glomeruli compared with control glomeruli, and further cluster analysis showed that these immune cells were predominantly macrophages. Analysis of differential gene expression in endothelial and mesangial cells of diabetic and control mice showed dynamic changes in the pattern of expressed genes, many of which are known to be involved in diabetic kidney disease. Moreover, gene expression analysis showed variable responses of individual cells to diabetic injury. CONCLUSIONS: Our findings demonstrate the ability of scRNA-seq analysis in isolated glomerular cells from diabetic and control mice to reveal dynamic changes in gene expression in diabetic kidneys, with variable responses of individual cells. Such changes, which might not be apparent in bulk transcriptomic analysis of glomerular cells, may help identify important pathophysiologic factors contributing to the progression of diabetic kidney disease.

A web server for comparative analysis of single-cell RNA-seq data.

Single cell RNA-Seq (scRNA-seq) studies profile thousands of cells in heterogeneous environments. Current methods for characterizing cells perform unsupervised analysis followed by assignment using a small set of known marker genes. Such approaches are limited to a few, well characterized cell types. We developed an automated pipeline to download, process, and annotate publicly available scRNA-seq datasets to enable large scale supervised characterization. We extend supervised neural networks to obtain efficient and accurate representations for scRNA-seq data. We apply our pipeline to analyze data from over 500 different studies with over 300 unique cell types and show that supervised methods outperform unsupervised methods for cell type identification. A case study highlights the usefulness of these methods for comparing cell type distributions in healthy and diseased mice. Finally, we present scQuery, a web server which uses our neural networks and fast matching methods to determine cell types, key genes, and more.

Exploring the single-cell RNA-seq analysis landscape with the scRNA-tools database.

As single-cell RNA-sequencing (scRNA-seq) datasets have become more widespread the number of tools designed to analyse these data has dramatically increased. Navigating the vast sea of tools now available is becoming increasingly challenging for researchers. In order to better facilitate selection of appropriate analysis tools we have created the scRNA-tools database (www.scRNA-tools.org) to catalogue and curate analysis tools as they become available. Our database collects a range of information on each scRNA-seq analysis tool and categorises them according to the analysis tasks they perform. Exploration of this database gives insights into the areas of rapid development of analysis methods for scRNA-seq data. We see that many tools perform tasks specific to scRNA-seq analysis, particularly clustering and ordering of cells. We also find that the scRNA-seq community embraces an open-source and open-science approach, with most tools available under open-source licenses and preprints being extensively used as a means to describe methods. The scRNA-tools database provides a valuable resource for researchers embarking on scRNA-seq analysis and records the growth of the field over time.

Parabiosis and single-cell RNA sequencing reveal a limited contribution of monocytes to myofibroblasts in kidney fibrosis.

Fibrosis is the common final pathway of virtually all chronic injury to the kidney. While it is well accepted that myofibroblasts are the scar-producing cells in the kidney, their cellular origin is still hotly debated. The relative contribution of proximal tubular epithelium and circulating cells, including mesenchymal stem cells, macrophages, and fibrocytes, to the myofibroblast pool remains highly controversial. Using inducible genetic fate tracing of proximal tubular epithelium, we confirm that the proximal tubule does not contribute to the myofibroblast pool. However, in parabiosis models in which one parabiont is genetically labeled and the other is unlabeled and undergoes kidney fibrosis, we demonstrate that a small fraction of genetically labeled renal myofibroblasts derive from the circulation. Single-cell RNA sequencing confirms this finding but indicates that these cells are circulating monocytes, express few extracellular matrix or other myofibroblast genes, and express many proinflammatory cytokines. We conclude that this small circulating myofibroblast progenitor population contributes to renal fibrosis by paracrine rather than direct mechanisms.

Population snapshots predict early haematopoietic and erythroid hierarchies.

The formation of red blood cells begins with the differentiation of multipotent haematopoietic progenitors. Reconstructing the steps of this differentiation represents a general challenge in stem-cell biology. Here we used single-cell transcriptomics, fate assays and a theory that allows the prediction of cell fates from population snapshots to demonstrate that mouse haematopoietic progenitors differentiate through a continuous, hierarchical structure into seven blood lineages. We uncovered coupling between the erythroid and the basophil or mast cell fates, a global haematopoietic response to erythroid stress and novel growth factor receptors that regulate erythropoiesis. We defined a flow cytometry sorting strategy to purify early stages of erythroid differentiation, completely isolating classically defined burst-forming and colony-forming progenitors. We also found that the cell cycle is progressively remodelled during erythroid development and during a sharp transcriptional switch that ends the colony-forming progenitor stage and activates terminal differentiation. Our work showcases the utility of linking transcriptomic data to predictive fate models, and provides insights into lineage development in vivo.

Clinical and parasitological factors in parasite persistence after treatment and clinical cure of cutaneous leishmaniasis.

BACKGROUND: The determinants of parasite persistence or elimination after treatment and clinical resolution of cutaneous leishmaniasis (CL) are unknown. We investigated clinical and parasitological parameters associated with the presence and viability of Leishmania after treatment and resolution of CL caused by L. Viannia. METHODS: Seventy patients who were treated with meglumine antimoniate (n = 38) or miltefosine (n = 32) and cured, were included in this study. Leishmania persistence and viability were determined by detection of kDNA and 7SLRNA transcripts, respectively, before, at the end of treatment (EoT), and 13 weeks after initiation of treatment in lesions and swabs of nasal and tonsillar mucosa. RESULTS: Sixty percent of patients (42/70) had evidence of Leishmania persistence at EoT and 30% (9/30) 13 weeks after treatment initiation. A previous episode of CL was found to be a protective factor for detectable Leishmania persistence (OR: 0.16, 95%CI: 0.03-0.92). kDNA genotyping could not discern differences between parasite populations that persisted and those isolated at diagnosis. CONCLUSIONS: Leishmania persist in skin and mucosal tissues in a high proportion of patients who achieved therapeutic cure of CL. This finding prompts assessment of the contribution of persistent infection in transmission and endemicity of CL, and in disease reactivation and protective immunity.

Nucleoid and cytoplasmic localization of small RNAs in Escherichia coli.

Bacterial small RNAs (sRNAs) regulate protein production by binding to mRNAs and altering their translation and degradation. sRNAs are smaller than most mRNAs but larger than many proteins. Therefore it is uncertain whether sRNAs can enter the nucleoid to target nascent mRNAs. Here, we investigate the intracellular localization of sRNAs transcribed from plasmids in Escherichia coli using RNA fluorescent in-situ hybridization. We found that sRNAs (GlmZ, OxyS, RyhB and SgrS) have equal preference for the nucleoid and cytoplasm, and no preferential localization at the cell membrane. We show using the gfp mRNA (encoding green fluorescent protein) that non-sRNAs can be engineered to have different proportions of nucleoid and cytoplasmic localization by altering their length and/or translation. The same localization as sRNAs was achieved by decreasing gfp mRNA length and translation, which suggests that sRNAs and other RNAs may enter the densely packed DNA of the nucleoid if they are sufficiently small. We also found that the Hfq protein, which binds sRNAs, minimally affects sRNA localization. Important implications of our findings for engineering synthetic circuits are: (i) sRNAs can potentially bind nascent mRNAs in the nucleoid, and (ii) localization patterns and distribution volumes of sRNAs can differ from some larger RNAs.

Parasitological Confirmation and Analysis of Leishmania Diversity in Asymptomatic and Subclinical Infection following Resolution of Cutaneous Leishmaniasis.

BACKGROUND: The contribution of individuals with subclinical infection to the transmission and endemicity of cutaneous leishmaniasis (CL) is unknown. Immunological evidence of exposure to Leishmania in residents of endemic areas has been the basis for defining the human population with asymptomatic infection. However, parasitological confirmation of subclinical infection is lacking. METHODS: We investigated the presence and viability of Leishmania in blood and non-invasive mucosal tissue samples from individuals with immunological evidence of subclinical infection in endemic areas for CL caused by Leishmania (Viannia) in Colombia. Detection of Leishmania kDNA was conducted by PCR-Southern Blot, and parasite viability was confirmed by amplification of parasite 7SLRNA gene transcripts. A molecular tool for genetic diversity analysis of parasite populations causing persistent subclinical infection based on PCR amplification and sequence analysis of an 82bp region between kDNA conserved blocks 1 and 2 was developed. PRINCIPAL FINDINGS: Persistent Leishmania infection was demonstrated in 40% (46 of 114) of leishmanin skin test (LST) positive individuals without active disease; parasite viability was established in 59% of these (27 of 46; 24% of total). Parasite burden quantified from circulating blood monocytes, nasal, conjunctival or tonsil mucosal swab samples was comparable, and ranged between 0.2 to 22 parasites per reaction. kDNA sequences were obtained from samples from 2 individuals with asymptomatic infection and from 26 with history of CL, allowing genetic distance analysis that revealed diversity among sequences and clustering within the L. (Viannia) subgenus. CONCLUSIONS: Our results provide parasitological confirmation of persistent infection among residents of endemic areas of L. (Viannia) transmission who have experienced asymptomatic infection or recovered from CL, revealing a reservoir of infection that potentially contributes to the endemicity and transmission of disease. kDNA genotyping establishes proof-of-principle of the feasibility of genetic diversity analysis in previously inaccessible and unexplored parasite populations in subclinically infected individuals.

Rhizomelic chondrodysplasia punctata with maternal systemic lupus erythromatosus.

We report Rhizomelic Chondrodysplasia Punctata (RDCP), a rare, autosomal recessive disorder with rhizomelic shortening of limbs, congenital cataracts and seizures but without any biochemical abnormality. The mother of the baby developed Systemic Lupus Erythromatosus (SLE) with Ro/SSA antibodies 11 months after delivery. Ro/SSA antibodies may generate calreticulin antibodies causing characteristic skeletal changes.

Natural infection of North African gundi (Ctenodactylus gundi) by Leishmania tropica in the focus of cutaneous leishmaniasis, Southeast Tunisia.

North African gundis (Ctenodactylus gundi) were trapped in the Leishmania (L.) tropica focus of cutaneous leishmaniasis, situated in southeast Tunisia and evaluated for Leishmania infection by real-time kinetoplast DNA polymerase chain reaction (PCR). Species identification was performed by internal transcribed spacer one (ITS1)-PCR-restriction fragment length polymorphism (RFLP) and high-resolution melting (HRM) analysis of the 7SL RNA gene. Real-time PCR on blood was positive in 6 of 13 (46.2%) tested gundis. Leishmania tropica was identified in five infected gundis and Leishmania major in one specimen. Alignments of the ITS-1 DNA sequences and 7S-HRM curves analysis indicated that similar genotypes were present in humans, a sandfly, and gundis from the same region suggesting a potential role of this rodent as reservoir host of L. tropica in southeast Tunisia.

Prevalence of canine leishmaniasis in Beichuan County, Sichuan, China and phylogenetic evidence for an undescribed Leishmania sp. in China based on 7SL RNA.

BACKGROUND: Leishmaniasis is a vector-borne disease, which is still endemic in the west and northwest area of China. Canines are the major reservoirs of Leishmania, the etiological agent of human visceral leishmaniasis. Phlebotomus chinensis is the main transmission vector of zoonotic visceral leishmaniasis (ZVL). METHODS: In this study, rK39 dip-stick, ELISA and PCR methods were used to investigate the prevalence of canine leishmaniasis (CanL) in Beichuan County, Sichuan Province, China. RESULTS: Among the 86 dogs which were included in the study, 13 dogs were positive using the dip-stick test (15.12%), while 8 dogs were positive using ELISA (9.30%) and 19 dogs were positive for PCR (22.03%). In total, 32 dogs were positive for one or more tests (37.21%). Interestingly, phylogenetic analysis based on the partial 7SL RNA fragment provided evidence that an undescribed Leishmania species, which is clearly a causative agent of CanL and human visceral leishmaniasis, does exist in China. This result is consistent with our previous study. CONCLUSIONS: Our work confirmed that canine leishmaniasis is still prevalent in Beichuan County. Further control is urgently needed, as canine leishmaniasis is of great public health importance. The phylogenetic analysis based on 7SL RNA segment provides evidence for the existence of an undescribed Leishmania sp. in China.

Salivary anti-Ro60 and anti-Ro52 antibody profiles to diagnose Sjogren's Syndrome.

Simple and non-invasive saliva-based diagnostics may be useful for the identification, understanding, and monitoring of autoimmune and infectious diseases. Previously, Luciferase Immunoprecipitation Systems (LIPS) were used for sensitive detection of patient serum autoantibodies in Sjögren's Syndrome (SjS), a chronic autoimmune disease affecting the salivary and lacrimal glands. Here we explored the ability of LIPS to diagnose SjS based on IgG autoantibodies in patient saliva. From LIPS testing, anti-Ro60 autoantibodies were detected in the saliva of 70% (19/27) of SjS patients with 96% specificity. Positive anti-Ro60 autoantibodies were also found in 70% of the matched serum samples (96% specificity). LIPS detected Ro52 autoantibodies in the saliva and serum of 67% of SjS patients with 100% specificity. Overall, the autoantibody titers in saliva were approximately 4000-fold lower by volume than serum, but still distinguished seropositive patients from controls. These results suggest that LIPS salivary-based testing for SjS autoantibodies is a practical alternative to serum and compatible with point-of-care testing.

7SL RNA is retained in HIV-1 minimal virus-like particles as an S-domain fragment.

HIV-1 is known to package several small cellular RNAs in addition to its genome. Previous work consistently demonstrated that the host structural RNA 7SL is abundant in HIV-1 virions but has yielded conflicting results regarding whether 7SL is present in minimal, assembly-competent virus-like particles (VLPs). Here, we demonstrate that minimal HIV-1 VLPs retain 7SL RNA primarily as an endoribonucleolytic fragment, referred to as 7SL remnant (7SLrem). Nuclease mapping showed that 7SLrem is a 111-nucleotide internal portion of 7SL, with 5' and 3' ends corresponding to unpaired loops in the 7SL two-dimensional structure. Analysis of VLPs comprised of different subsets of Gag domains revealed that all NC-positive VLPs contained intact 7SL while the presence of 7SLrem correlated with the absence of the NC domain. Because 7SLrem, which maps to the 7SL S domain, was not detectable in infected cells, we propose a model whereby the species recruited to assembling VLPs is intact 7SL RNA, with 7SLrem produced by an endoribonuclease in the absence of NC. Since recruitment of 7SL RNA was a conserved feature of all tested minimal VLPs, our model further suggests that 7SL's recruitment is mediated, either directly or indirectly, through interactions with conserved features of all tested VLPs, such as the C-terminal domain of CA.

Oral manifestations of Sjögren's syndrome.

Sjögren's syndrome is a common autoimmune rheumatic disease. The most common symptoms of Sjögren's syndrome are extreme tiredness, along with dry eyes (keratoconjunctivitis sicca) and dry mouth (xerostomia). Saliva plays an essential role in numerous functions of the mouth. Xerostomia can be caused by medications, chronic diseases like Sjögren's syndrome, and medical treatments, such as radiation therapy and bone marrow transplant. Xerostomia can eventually lead to difficulty in swallowing, severe and progressive tooth decay, or oral infections. Despite having excellent oral hygiene, individuals with Sjögren's syndrome have elevated levels of dental caries, along with the loss of many teeth, early in the disease. Sjögren's syndrome alters the protein profile and brings about a change in the composition of saliva. There is an increase in the levels of lactoferrin, beta(2)-microglobulin, sodium, lysozyme C, and cystatin C, and a decrease in salivary amylase and carbonic anhydrase. Up to 90% of individuals with Sjögren's syndrome have antibodies targeting the Ro 60 and La autoantigens. Natural aging, regardless of Sjögren's syndrome, is also another factor that brings about a significant change in the composition of saliva. The most prevailing cause of xerostomia in elderly persons is the use of anticholinergic medications. Currently, there is no cure for Sjögren's syndrome, and treatment is mainly palliative.

7SL RNA, but not the 54-kd signal recognition particle protein, is an abundant component of both infectious HIV-1 and minimal virus-like particles.

The virion incorporation of 7SL, the RNA component of the host signal recognition particle (SRP), has been shown for several simple retroviruses. Data here demonstrate that 7SL is also packaged by HIV-1, in sevenfold molar excess of genomic RNA. Viral determinants of HIV-1 genome and primer tRNA packaging were not required for 7SL incorporation, as virus-like particles with only minimal assembly components efficiently packaged 7SL. The majority of 7SL within cells resides in ribonucleoprotein complexes bound by SRP proteins, and most SRP protein exists in signal recognition particles. However, Western blot comparison of virion and cell samples revealed that there is at least 25-fold less SRP p54 protein per 7SL RNA in HIV-1 particles than in cells. Comparing 7SL:actin mRNA ratios in virions and cells revealed that 7SL RNA appears selectively enriched in virions.

Nonrandom packaging of host RNAs in moloney murine leukemia virus.

Moloney murine leukemia virus (MLV) particles contain both viral genomic RNA and an assortment of host cell RNAs. Packaging of virus-encoded RNA is selective, with virions virtually devoid of spliced env mRNA and highly enriched for unspliced genome. Except for primer tRNA, it is unclear whether packaged host RNAs are randomly sampled from the cell or specifically encapsidated. To address possible biases in host RNA sampling, the relative abundances of several host RNAs in MLV particles and in producer cells were compared. Using 7SL RNA as a standard, some cellular RNAs, such as those of the Ro RNP, were found to be enriched in MLV particles in that their ratios relative to 7SL differed little, if at all, from their ratios in cells. Some RNAs were underrepresented, with ratios relative to 7SL several orders of magnitude lower in virions than in cells, while others displayed intermediate values. At least some enriched RNAs were encapsidated by genome-defective nucleocapsid mutants. Virion RNAs were not a random sample of the cytosol as a whole, since some cytoplasmic RNAs like tRNA(Met) were vastly underrepresented, while U6 spliceosomal RNA, which functions in the nucleus, was enriched. Real-time PCR demonstrated that env mRNA, although several orders of magnitude less abundant than unspliced viral RNA, was slightly enriched relative to actin mRNA in virions. These data demonstrate that certain host RNAs are nearly as enriched in virions as genomic RNA and suggest that Psi- mRNAs and some other host RNAs may be specifically excluded from assembly sites.