Dataset Information

DCAF1 regulates Treg senescence via the ROS axis during immunological aging

ABSTRACT: Total RNA of regulatory T (Treg) cells and conventional CD4 T (Tconv) cells from young WT (2 months), aged WT (>18 months) and young Dcaf1-deficient mice was extracted and RNA-seq libraries were generated. Reads (32-45 Million reads per sample) were analyzed with Salmon software to align and quantify the transcript expression. R packages in Bioconductor, tximport and tximportData were used to aggregate transcript-level quantifications to the gene level, with the R package biomaRt for gene and transcripts mapping. The option "lengthScaledTPM" for countsFromAbundance in tximport was used to obtain the estimated counts at the gene level using abundance estimates scaled based on the average transcript length over samples and the library size. The option "lengthScaledTPM" for countsFromAbundance in tximport was used to obtain the estimated counts at the gene level using abundance estimates scaled based on the average transcript length over samples and the library size. For the differential expression (DE) analysis of RNA-seq data, gene-level count matrix was passed into by DESeq2 package as input directly from the tximport package. The normalized gene expression data was retrieved from DESeq2 analysis after regulated log (rlog) transformation (‘rlog’ in DESeq2). The z-score at gene-level average of normalized expression matrix was used to generate heatmap in Gene-E from Broad Institute (www.broadinstitute.org/ /GENE-E/). Gene Set Enrichment Analysis was performed using the Java application available from Broad Institute (www.broadinstitute.org/gsea/). Gene set databases including Hallmarks (h.all.v6.1.symbols.gmt) and KEGG (c2.cp.kegg.v6.1.symbols.gmt) from the Molecular Signatures Database (MSigDB) were used in the analysis. The aging-program gene set was from DEMAGALHAES_AGING_UP in MSigDB. One thousand gene set permutations were performed. FDR<0.05 was used for enriched terms, as is recommended when performing permutations by gene set.

ORGANISM(S): Mus musculus

PROVIDER: GSE130419 | GEO | 2020/10/10

REPOSITORIES: GEO

ACCESS DATA

Similar Datasets

Project description:We investigated the ability of HDAC inhibitors (HDACi) to target CML stem cells. Treatment with HDACi combined with IM effectively induced apoptosis in quiescent CML progenitors resistant to elimination by IM alone, and eliminated CML stem cells capable of engrafting immunodeficient mice. In vivo administration of HDACi with IM markedly diminished LSC in a transgenic mouse model of CML. The interaction of IM and HDACi inhibited genes regulating hematopoietic stem cell maintenance and survival. HDACi treatment represents a novel and effective strategy to target LSC in CML patients receiving tyrosine kinase inhibitors. CML CD34+CD38- cells were selected using flow cytometry sorting and treated with IM, LBH and the combination of IM and LBH or cultured without exposure to drugs (controls) for 24 hours (n=3 each). Total RNA from 5000 cells was extracted using the RNeasy kit (Qiagen), amplified and labeled using GeneChip Two-Cycle Target Labeling and Control Reagents (Affymetrix, Santa Clara, CA). 15 Âµg cRNA from each sample was hybridized to Affymetrix GeneChip Human Genome U133 Plus 2.0 Arrays. Microarray data analyses were performed using R (version 2.9) with genomic analysis packages from Bioconductor (version 2.4). Expression data were normalized using the robust multiarray average (RMA) algorithm, with background adjustment, quantile normalization and median polish summarization. Probesets with low expression levels or low variability across samples were filtered. For genes with multiple probesets, the gene level expression was set to be the median of the probesets. Linear regression was used to model the gene expression with the consideration of 2x2 factorial design and matched samples. Differentially expressed genes were identified by calculating empirical Bayes moderated t-statistic, and p-values were adjusted by FDR using the âLIMMAâ package. Gene Set Enrichment Analyses (GSEA) was performed using GSEA software version 2.04 [http://www.broadinstitute.org/gsea/] to detect enrichment of predetermined gene sets using t-scores and gene sets in C2 (curated gene sets) category from the Molecular Signature Database (MsigDB). Gene sets representing common functional categories were categorized and grouped. We also analyzed enrichment of gene sets with common transcription factor binding sites (586 sets) from MsigDB.

Project description:Purpose: muscle transcriptomics of subjects supplemented with Urolithin A at different doses or placebo for 4 month. Method: RNA-seq was performed using Illumina HiSeq 4000 sequencing; single read 1 x 50 bp . The quantification of mRNA from the RNA-seq FASTQ files was performed using Salmon. Sample-wise quant.sf files containing raw transcript-level read estimates were read into R, v. 4.0.3 and were combined into a data matrix. Transcripts with very low total counts (< 10) across all samples were filtered out. The data was transformed using the variance stabilizing transformation (VST) method of R package DESeq2, v. 1.30.0. Top 10,000 transcripts with the highest variance across all samples were used for principal component analysis (PCA) using DESeq2. Data transformation and PCA was also done separately for each treatment group. Based on the PCAs, probable outlier samples were excluded and new PCAs were plotted without these samples. The raw transcript-level read count estimates were read in R and summarized to gene-level counts based on the provided transcript and gene ID annotations using summarizeToGene function of R package tximport, v. 1.18.0. DESeqDataSetFromTximport function of DESeq2 was then used for constructing a DESeqDataSet object for DE analysis. Pre-filtering was applied before the DE analysis by excluding genes with < 10 total counts across samples. Subset DE analysis was performed, contrasting Visit time D120 with baseline (BL) and by adjusting for the subject effect. The normalization and DE analysis was done separately for the three different treatment groups. Independent filtering option of DESeq2 was enabled (default), filtering out genes with very low counts and thus unlikely to show significant evidence. R package biomaRt. v. 2.46.0 was used for annotating the results with HGNC gene symbols, gene descriptions and gene biotypes. DESeq2-normalised expression values of all the samples in the given comparison were added to the result tables. Non-adjusted p-value 0.05 was used to filter the results by statistical significance. Results were also generated using DESeq2 function lfcShrink that allows for the shrinkage of the log2 fold change (LFC) estimates toward zero when the information for a gene is low (such as in those cases with low counts or high dispersion values) but has little effect on genes with high counts. The shrinked log2FC values were subsequently used for visualisation and ranking the genes.

Project description:Tight regulation of hematopoietic stem cell (HSC) homeostasis is essential for life-long hematopoiesis, for preventing blood cancers and for averting bone marrow failure. The underlying mechanisms are incompletely understood. Here, we identify production of inositol-tetrakisphosphate (IP4) by inositoltrisphosphate 3-kinase B (ItpkB) as essential for HSC quiescence and function. Young ItpkB-/- mice accumulated phenotypic HSC and showed extramedullary hematopoiesis. ItpkB-/- HSC were less quiescent and proliferated more than wildtype controls. They downregulated quiescence and stemness associated mRNAs, but upregulated activation, oxidative metabolism, protein synthesis and lineage associated transcripts. Although they showed no significant homing defects, ItpkB-/- HSC had a severely reduced competitive long-term repopulating potential. Aging ItpkB-/- mice lost hematopoietic stem and progenitor cells and died with severe anemia. Wildtype HSC normally repopulated ItpkB-/- hosts, incidating a HSC-intrinsic ItpkB requirement. ItpkB-/- HSC had reduced cobblestone-area forming cell activity in vitro and showed increased stem-cell-factor activation of the phosphoinositide 3-kinase (PI3K) effector Akt, reversed by exogenous provision of the known PI3K/Akt antagonist IP4. They also showed transcriptome changes consistent with hyperactive Akt/mTOR signaling. Thus, we propose that ItpkB ensures HSC quiescence by limiting cytokine-induced PI3K signaling in HSC. For each of 3 replicate ItpkB-/- or wt samples, we enriched Lin- cells from BM of 4 pooled age-matched mice with Rapidspheres (Stemcell Technologies), FACS-sorted ?10,000 LSK CD34-CD150+CD48-Flk2- LT-HSC into lysis buffer and prepared RNA with RNeasy Micro kits (Quiagen). RNA sequencing was done using an Illumina HISeq Analyzer 2000, Casava v1.8.2 genome analyzer pipeline, TopHat v1.4.1/Bowtie2 genome alignment and Partek v6.6 mRNA annotation software. Statistical analyses were done with edgeR (Bioconductor package), excluding genes with false discovery rates >0.15, fold-change magnitudes ?1.4 and log2(counts per million) ?4 to avoid undefined values and the poorly defined log fold-changes for low counts close to 0. Unsupervised clustering of 441 significantly changed genes was done with dChip using rank correlation and a centroid linkage method. Scatter plots were generated in Spotfire. GSEA was performed with gene set permutation, using gene sets from MSigDB (www.broadinstitute.org/gsea/msigdb/index.jsp) or manually curated from, excluding genes without HUGO approved symbols

Project description:Purpose: Cachexia is a metabolic phenotype affecting many cancer patients. The immune environment in adipose tissue during the development of cancer associated cachexia has not been characterized. To help address this gap, we isolated visceral adipose tissue from cancer-bearing and sham-operated mice and performed bulk RNA-sequencing. Methods: mRNA expression was evaluated by RNA-sequencing in visceral adipose tissue (VAT) of 10-week old mice which had been subjected to sham surgery (sham) or had received orthotopic pancreatic tumors four weeks prior to harvest (KPC). Paired-end 150 bp single-index sequencing was performed using an Illumina HiSeq system. Transcript quantification was performed with Salmon, using quasi-mapping mode with selective alignment (validateMappings). We set the rangeFactorizationBins to 4 and passed both the seqBias and gcBias flags to correct for sequence-specific bias and fragment-level GC bias, respectively. Transcript quantifications were collapsed to gene-level quantifications using the tximport package and the Ensembl based annotation package EnsDb.Mmusculus.v79. Differential expression analysis was then performed using DESeq2 with fold change shrinkage by ashr. For downstream analysis and visualization, only genes with a mean TPM (transcripts per million) > 2 in VAT from both sham-operated and tumor-bearing mice were considered. Results: We identified over 11,000 genes which were expressed at an average of at least 2 TPM (transcripts per million) in both the sham and KPC adipose tissue samples. Differential expression analysis with DESeq revealed 161 upregulated genes (p < 0.05) and 27 downregulated genes (p < 0.05) in the adipose tissue of cancer bearing mice compared to that of sham operated mice. Gene set enrichment analysis suggested an increase in inflammatory pathways, including type 1 and type 2 interferon pathways, in the adipose tissue of cancer bearing mice. Conclusions: This study represents a comprehensive characterization of adipose tissue in the healthy vs. cachexia setting. The RNA-seq data suggests increased immune infiltration into adipose tissue in the setting of a distant pancreatic tumor. Together with flow cytometry and cytokine abundance analyses, we conclude that reguatory Th2 immunity is maintained in the adipose tissue during this model of cancer associated cachexia.

Dataset Information

DCAF1 regulates Treg senescence via the ROS axis during immunological aging

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure