Single-cell transcriptomics is rapidly advancing our understanding of the cellular composition of complex tissues and organisms. A major limitation in most analysis pipelines is the reliance on ...manual annotations to determine cell identities, which are time-consuming and irreproducible. The exponential growth in the number of cells and samples has prompted the adaptation and development of supervised classification methods for automatic cell identification.
Here, we benchmarked 22 classification methods that automatically assign cell identities including single-cell-specific and general-purpose classifiers. The performance of the methods is evaluated using 27 publicly available single-cell RNA sequencing datasets of different sizes, technologies, species, and levels of complexity. We use 2 experimental setups to evaluate the performance of each method for within dataset predictions (intra-dataset) and across datasets (inter-dataset) based on accuracy, percentage of unclassified cells, and computation time. We further evaluate the methods' sensitivity to the input features, number of cells per population, and their performance across different annotation levels and datasets. We find that most classifiers perform well on a variety of datasets with decreased accuracy for complex datasets with overlapping classes or deep annotations. The general-purpose support vector machine classifier has overall the best performance across the different experiments.
We present a comprehensive evaluation of automatic cell identification methods for single-cell RNA sequencing data. All the code used for the evaluation is available on GitHub ( https://github.com/tabdelaal/scRNAseq_Benchmark ). Additionally, we provide a Snakemake workflow to facilitate the benchmarking and to support the extension of new methods and new datasets.
RELATIVE OF EARLY FLOWERING 6 (REF6/JMJ12), a Jumonji C (JmjC)-domain-containing H3K27me3 histone demethylase, finds its target loci in Arabidopsis genome by directly recognizing the CTCTGYTY motif ...via its zinc-finger (ZnF) domains. REF6 tends to bind motifs located in active chromatin states that are depleted for heterochromatic modifications. However, the underlying mechanism remains unknown. Here, we show that REF6 preferentially bind to hypo-methylated CTCTGYTY motifs in vivo, and that CHG methylation decreases REF6 DNA binding affinity in vitro. In addition, crystal structures of ZnF-clusters in complex with DNA oligonucleotides reveal that 5-methylcytosine is unfavorable for REF6 binding. In drm1 drm2 cmt2 cmt3 (ddcc) quadruple mutants, in which non-CG methylation is significantly reduced, REF6 can ectopically bind a small number of new target loci, most of which are located in or neighbored with short TEs in euchromatic regions. Collectively, our findings reveal that DNA methylation, likely acting in combination with other epigenetic modifications, may partially explain why REF6 binding is depleted in heterochromatic loci.
Selection of high-quality embryos is important to achieve successful pregnancy in assisted reproductive technology (ART). Recently, it has been debated whether RNA-sequencing (RNA-Seq) should be ...applied to ART to predict embryo quality. However, information on genes that can serve as markers for pregnant expectancy is limited. Furthermore, there is no information on which transcriptome of trophectoderm (TE) or inner cell mass (ICM) is more highly correlated with pregnant expectancy. Here, we performed RNA-Seq analysis of TE and ICM of human blastocysts, the pregnancy expectation of which was retrospectively determined using the clinical outcomes of 1,890 cases of frozen-thawed blastocyst transfer. We identified genes that were correlated with the expected pregnancy rate in ICM and TE, respectively, with a larger number of genes identified in TE than in ICM. Downregulated genes in the TE of blastocysts that were estimated to have lower expectation of pregnancy included tight junction-related genes such as CXADR and ATP1B1, which have been implicated in peri-implantation development. Moreover, we identified dozens of differentially expressed genes by regrouping the blastocysts based on the maternal age and the Gardner score. Additionally, we showed that aneuploidy estimation using RNA-Seq datasets does not correlate with pregnancy expectation. Thus, our study provides an expanded list of candidate genes for the prediction of pregnancy in human blastocyst embryos.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Primary cutaneous CD8+ aggressive epidermotropic cytotoxic T-cell lymphoma (pcAECyTCL) is a rare variant of cutaneous T-cell lymphoma with an aggressive clinical course and a very poor prognosis. ...Until now, neither a systematic characterization of genetic alterations driving pcAECyTCL has been performed, nor effective therapeutic regimes for patients have been defined. Here, we present the first highresolution genetic characterization of pcAECyTCL by using wholegenome and RNA sequencing. Our study provides a comprehensive description of genetic alterations (i.e., genomic rearrangements, copy number alterations and small-scale mutations) with pathogenic relevance in this lymphoma, including events that recurrently impact genes with important roles in the cell cycle, chromatin regulation and the JAKSTAT pathway. In particular, we show that mutually exclusive structural alterations involving JAK2 and SH2B3 predominantly underlie pcAECyTCL. In line with the genomic data, transcriptome analysis uncovered upregulation of the cell cycle, JAK2 signaling, NF-κB signaling and a high inflammatory response in this cancer. Functional studies confirmed oncogenicity of JAK2 fusions identified in pcAECyTCL and their sensitivity to JAK inhibitor treatment. Our findings strongly suggest that overactive JAK2 signaling is a central driver of pcAECyTCL, and consequently, patients with this neoplasm would likely benefit from therapy with JAK2 inhibitors such as Food and Drug Adminstration-approved ruxolitinib.
Ovo-like transcriptional repressor 1 (OVOL1) is a key mediator of epithelial lineage determination and mesenchymal-epithelial transition (MET). The cytokines transforming growth factor-β (TGF-β) and ...bone morphogenetic proteins (BMP) control the epithelial-mesenchymal plasticity (EMP) of cancer cells, but whether this occurs through interplay with OVOL1 is not known. Here, we show that OVOL1 is inversely correlated with the epithelial-mesenchymal transition (EMT) signature, and is an indicator of a favorable prognosis for breast cancer patients. OVOL1 suppresses EMT, migration, extravasation, and early metastatic events of breast cancer cells. Importantly, BMP strongly promotes the expression of OVOL1, which enhances BMP signaling in turn. This positive feedback loop is established through the inhibition of TGF-β receptor signaling by OVOL1. Mechanistically, OVOL1 interacts with and prevents the ubiquitination and degradation of SMAD family member 7 (SMAD7), which is a negative regulator of TGF-β type I receptor stability. Moreover, a small-molecule compound 6-formylindolo(3,2-b)carbazole (FICZ) was identified to activate OVOL1 expression and thereby antagonizing (at least in part) TGF-β-mediated EMT and migration in breast cancer cells. Our results uncover a novel mechanism by which OVOL1 attenuates TGF-β/SMAD signaling and maintains the epithelial identity of breast cancer cells.
Exacerbations are major contributors to morbidity and mortality in patients with chronic obstructive pulmonary disease (COPD), and respiratory bacterial and viral infections are an important trigger. ...However, using conventional diagnostic techniques, a causative agent is not always found. Metagenomic next-generation sequencing (mNGS) allows analysis of the complete virome, but has not yet been applied in COPD exacerbations.
To study the respiratory virome in nasopharyngeal samples during COPD exacerbations using mNGS.
88 nasopharyngeal swabs from 63 patients from the Bergen COPD Exacerbation Study (2006-2010) were analysed by mNGS and in-house qPCR for respiratory viruses. Both DNA and RNA were sequenced simultaneously using an Illumina library preparation protocol with in-house adaptations.
By mNGS, 24/88 samples tested positive. Sensitivity and specificity, as compared with PCR, were 96% and 98% for diagnostic targets (23/24 and 1093/1120, respectively). Additional viral pathogens detected by mNGS were herpes simplex virus type 1 and coronavirus OC43. A positive correlation was found between Cq value and mNGS viral normalized species reads (log value) (p = 0.002). Patients with viral pathogens had lower percentages of bacteriophages (p<0.001). No correlation was found between viral reads and clinical markers.
The mNGS protocol used was highly sensitive and specific for semi-quantitative detection of respiratory viruses. Excellent negative predictive value implicates the power of mNGS to exclude any pathogenic respiratory viral infectious cause in one test, with consequences for clinical decision making. Reduced abundance of bacteriophages in COPD patients with viral pathogens implicates skewing of the virome during infection, with potential consequences for the bacterial populations, during infection.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
In Huntington disease, cellular toxicity is particularly caused by toxic protein fragments generated from the mutant huntingtin (HTT) protein. By modifying the HTT protein, we aim to reduce ...proteolytic cleavage and ameliorate the consequences of mutant HTT without lowering total HTT levels. To that end, we use an antisense oligonucleotide (AON) that targets HTT pre-mRNA and induces partial skipping of exon 12, which contains the critical caspase-6 cleavage site. Here, we show that AON-treatment can partially restore the phenotype of YAC128 mice, a mouse model expressing the full-length human HTT gene including 128 CAG-repeats. Wild-type and YAC128 mice were treated intracerebroventricularly with AON12.1, scrambled AON or vehicle starting at 6 months of age and followed up to 12 months of age, when MRI was performed and mice were sacrificed. AON12.1 treatment induced around 40% exon skip and protein modification. The phenotype on body weight and activity, but not rotarod, was restored by AON treatment. Genes differentially expressed in YAC128 striatum changed toward wild-type levels and striatal volume was preserved upon AON12.1 treatment. However, scrambled AON also showed a restorative effect on gene expression and appeared to generally increase brain volume.
Abstract
Objective
To identify OA subtypes based on cartilage transcriptomic data in cartilage tissue and characterize their underlying pathophysiological processes and/or clinically relevant ...characteristics.
Methods
This study includes n = 66 primary OA patients (41 knees and 25 hips), who underwent a joint replacement surgery, from which macroscopically unaffected (preserved, n = 56) and lesioned (n = 45) OA articular cartilage were collected Research Arthritis and Articular Cartilage (RAAK) study. Unsupervised hierarchical clustering analysis on preserved cartilage transcriptome followed by clinical data integration was performed. Protein–protein interaction (PPI) followed by pathway enrichment analysis were done for genes significant differentially expressed between subgroups with interactions in the PPI network.
Results
Analysis of preserved samples (n = 56) resulted in two OA subtypes with n = 41 (cluster A) and n = 15 (cluster B) patients. The transcriptomic profile of cluster B cartilage, relative to cluster A (DE-AB genes) showed among others a pronounced upregulation of multiple genes involved in chemokine pathways. Nevertheless, upon investigating the OA pathophysiology in cluster B patients as reflected by differentially expressed genes between preserved and lesioned OA cartilage (DE-OA-B genes), the chemokine genes were significantly downregulated with OA pathophysiology. Upon integrating radiographic OA data, we showed that the OA phenotype among cluster B patients, relative to cluster A, may be characterized by higher joint space narrowing (JSN) scores and low osteophyte (OP) scores.
Conclusion
Based on whole-transcriptome profiling, we identified two robust OA subtypes characterized by unique OA, pathophysiological processes in cartilage as well as a clinical phenotype. We advocate that further characterization, confirmation and clinical data integration is a prerequisite to allow for development of treatments towards personalized care with concurrently more effective treatment response.
In contrast to mouse, human female germ cells develop asynchronously. Germ cells transition to meiosis, erase genomic imprints, and reactivate the X chromosome. It is unknown if these events all ...appear asynchronously, and how they relate to each other. Here we combine exome sequencing of human fetal and maternal tissues with single-cell RNA-sequencing of five donors. We reconstruct full parental haplotypes and quantify changes in parental allele-specific expression, genome-wide. First we distinguish primordial germ cells (PGC), pre-meiotic, and meiotic transcriptional stages. Next we demonstrate that germ cells from various stages monoallelically express imprinted genes and confirm this by methylation patterns. Finally, we show that roughly 30% of the PGCs are still reactivating their inactive X chromosome and that this is related to transcriptional stage rather than fetal age. Altogether, we uncover the complexity and cell-to-cell heterogeneity of transcriptional and epigenetic remodeling in female human germ cells.
Primary cutaneous anaplastic large cell lymphoma (pcALCL), a hematological neoplasm caused by skin-homing CD30+ malignant T cells, is part of the spectrum of primary cutaneous CD30+ ...lymphoproliferative disorders. To date, only a small number of molecular alterations have been described in pcALCL and, so far, no clear unifying theme that could explain the pathogenetic origin of the disease has emerged among patients. In order to clarify the pathogenetic basis of pcALCL, we performed high-resolution genetic profiling (genome/transcriptome) of this lymphoma (n=12) by using whole-genome sequencing, whole-exome sequencing and RNA sequencing. Our study, which uncovered novel genomic rearrangements, copy number alterations and small-scale mutations underlying this malignancy, revealed that the cell cycle, T-cell physiology regulation, transcription and signaling via the PI-3-K, MAPK and G-protein pathways are cellular processes commonly impacted by molecular alterations in patients with pcALCL. Recurrent events affecting cancer-associated genes included deletion of PRDM1 and TNFRSF14, gain of EZH2 and TNFRSF8, small-scale mutations in LRP1B, PDPK1 and PIK3R1 and rearrangements involving GPS2, LINC-PINT and TNK1. Consistent with the genomic data, transcriptome analysis uncovered upregulation of signal transduction routes associated with the PI-3-K, MAPK and G-protein pathways (e.g., ERK, phospholipase C, AKT). Our molecular findings suggest that inhibition of proliferation-promoting pathways altered in pcALCL (particularly PI-3-K/AKT signaling) should be explored as potential alternative therapy for patients with this lymphoma, especially, for cases that do not respond to first-line skin-directed therapies or with extracutaneous disease.