Tag: DEG

Unsuccessful DE analysis using limma

This might be a bit long, please bare with me. I’m conducting a differential expression analysis using limma – voom. My comparison is regarding response vs non-response to a cancer drug. However, I’m not getting any DE genes, absolute zeros. Someone here once recommended not to use contrast matrix for…

Continue Reading Unsuccessful DE analysis using limma

Problem in getting SNAI1 gene in prostate DEG table

Problem in getting SNAI1 gene in prostate DEG table 0 Hello everyone, would you please help me I have got prostate DEGs by TCGAbiolinks package but I couldn’t get SNAI1 gene in my DEG table. I wonder why it does happen. First I normalized it and then I filtered it,…

Continue Reading Problem in getting SNAI1 gene in prostate DEG table

Mitogenome-wise codon usage pattern from comparative analysis of the first mitogenome of Blepharipa sp. (Muga uzifly) with other Oestroid flies

Outcome of DNA sequencing, assembly, and validation In this study, initially total DNA was isolated from the finely chopped, full-grown pupa of Blepharipa sp. The NanoDrop spectrophotometer (1294 ng/μl) and the Qubit fluorometer (732.8 ng/μl) both found that the concentration of total DNA in the sample at an optimum level for mitochondrial DNA enrichment. The Tape Station profile showed…

Continue Reading Mitogenome-wise codon usage pattern from comparative analysis of the first mitogenome of Blepharipa sp. (Muga uzifly) with other Oestroid flies

Extremely different results for both EdgeR and DESeq2 analysis

Extremely different results for both EdgeR and DESeq2 analysis 1 @373f98d7 Last seen 23 hours ago Singapore Dear all, Upon comparing my results for the analysis between DESeq2 and EdgeR, I have realized that the 2 results obtained after DEG analysis are extremely different from each other. The thresholds I…

Continue Reading Extremely different results for both EdgeR and DESeq2 analysis

Condition shown as ‘Untreated vs. Treated’ in output of results() in DESeq2

Condition shown as ‘Untreated vs. Treated’ in output of results() in DESeq2 1 @bffcbc5f Last seen 16 hours ago United States of America I am trying to find differentially expressed genes using DESeq2 on some RNA-seq data. In the pheno data, there is a column named ‘condition’ with factored values…

Continue Reading Condition shown as ‘Untreated vs. Treated’ in output of results() in DESeq2

RNAseq data DEG analysis – DESeq2 normalized data

RNAseq data DEG analysis – DESeq2 normalized data 1 1) You can’t use because those data are already normalized and log-transformed. 3) RSEM expected_count is best to start off with for differential expression. Login before adding your answer. Traffic: 2089 users visited in the last hour Read more here: Source…

Continue Reading RNAseq data DEG analysis – DESeq2 normalized data

dataframe – uwot is throwing an error running the Monocle3 R package’s “find_gene_module()” function, likely as an issue with how my data is formatted

I am trying to run the Monocle3 function find_gene_modules() on a cell_data_set (cds) but am getting a variety of errors in this. I have not had any other issues before this. I am working with an imported Seurat object. My first error came back stating that the number of rows…

Continue Reading dataframe – uwot is throwing an error running the Monocle3 R package’s “find_gene_module()” function, likely as an issue with how my data is formatted

DataViz question for DEG and scRNA-seq

DataViz question for DEG and scRNA-seq 1 Maybe a dumb question, but I recently came across a single-cell paper that showed “volcano” plots for all clusters in one image by filtering for all genes that were FDR <0.05, the x-axis being represented by different colors characterizing different clusters, the y-axis…

Continue Reading DataViz question for DEG and scRNA-seq

Key Players Identified Through Comparative Transcriptomics.ZIP

With long reproductive timescales, large complex genomes, and a lack of reliable reference genomes, understanding gene function in conifers is extremely challenging. Consequently, our understanding of which genetic factors influence the development of reproductive structures (cones) in monoecious conifers remains limited. Genes with inferred roles in conifer reproduction have mostly…

Continue Reading Key Players Identified Through Comparative Transcriptomics.ZIP

High-Throughput Transcriptome Analysis for Investigating Host-Pathogen Interactions

The protocol presented here describes a complete pipeline to analyze RNA-sequencing transcriptome data from raw reads to functional analysis, including quality control and preprocessing steps to advanced statistical analytical approaches. Welcome to the protocol of high-throughput transcriptome analysis for investigating host-pathogen interactions. This protocol is divided in the following steps….

Continue Reading High-Throughput Transcriptome Analysis for Investigating Host-Pathogen Interactions

In RNA-Seq data modeling process, glmfit function in R used raw RNA-Seq count data?

Hello, I’m newbie in RNA-Seq analysis process. When I processed a RNA-Seq analysis, there are some questions. If you guys have a time, please let me know. In RNA-Seq data analysis (i.e DEG analysis), glmfit function in R used raw RNA-Seq count data for modelling? I TMM normalized for correct…

Continue Reading In RNA-Seq data modeling process, glmfit function in R used raw RNA-Seq count data?

TIP4P/2005f implementation

Hi, I was trying to implement TIP4P/2005f force field in CP2K. I have setup the relevant forcefield and geometry section in the following manner: &FORCE_EVAL  METHOD FIST  &MM    &FORCEFIELD      &SPLINE        EMAX_SPLINE 1.0E+128        RCUT_NB 12.0      &END SPLINE      &BEND   …

Continue Reading TIP4P/2005f implementation

Profiling and functional characterization of maternal mRNA translation during mouse maternal-to-zygotic transition

INTRODUCTION Mammalian life starts with the fusion of two terminally differentiated gametes, sperm and oocyte, resulting in a totipotent zygote. After going through preimplantation development, the zygote reaches blastocyst before implantation. The two most important events taking place during preimplantation development are zygotic genome activation (ZGA) and the first cell…

Continue Reading Profiling and functional characterization of maternal mRNA translation during mouse maternal-to-zygotic transition

Identification of a four-gene signature & PTC.

Introduction Thyroid carcinoma (THCA) is the most common type of endocrine malignancy and its incidence is increasing.1 Based on its histopathological characteristics, thyroid carcinoma can be classified into multiple subtypes, such as papillary thyroid carcinoma (PTC), follicular thyroid carcinoma, and anaplastic thyroid carcinoma.2 PTC is the most common subtype of…

Continue Reading Identification of a four-gene signature & PTC.

DESeq2 input from GDAC firehose

Hi guys, I hope you are fine. I’m not good in English so if you couldn’t understand my question, please feel free to reply. I’m a beginner of bioinformatics. I want to practice differential expressed gene (DEG) analysis in R. The RNA seq data I used was downloaded from broad…

Continue Reading DESeq2 input from GDAC firehose

Identification of differentially expressed genes in AF

Defeng Pan,1,&ast; Yufei Zhou,2,&ast; Shengjue Xiao,1,&ast; Yue Hu,3,&ast; Chunyan Huan,1 Qi Wu,1 Xiaotong Wang,1 Qinyuan Pan,1 Jie Liu,1 Hong Zhu1 1Department of Cardiology, The Affiliated Hospital of Xuzhou Medical University, Xuzhou, Jiangsu, 221004, People’s Republic of China; 2Department of Cardiology, Shanghai Institute of Cardiovascular Diseases, Zhongshan Hospital and Institutes of…

Continue Reading Identification of differentially expressed genes in AF

AccuClear Ultra High Sensitivity dsDNA Quantitation Kit with DNA Standard, trial size (200 assays) 31028-T

Description: trial size (200 assays), AccuClear Ultra High Sensitivity dsDNA Quantitation Kit with DNA Standard Category: Nucleic acid quantitation Shipping temperature: 15&deg, C, C to 30&deg Shipping conditions: n/a Storage: 2&deg, C, C to 8&deg Gene target: AccuClear Sensitivity dsDNA Quantitation Kit with DNA trial size (200 assays) Short name:…

Continue Reading AccuClear Ultra High Sensitivity dsDNA Quantitation Kit with DNA Standard, trial size (200 assays) 31028-T

Monocle3 differential expression failed when active.assay is not “RNA”

after run estimate_size_factors, data with active.assay = ‘integrated’ works too, but no deg in the result. > [email protected] = ‘integrated’ > cds_raw <- as.cell_data_set(seurat_object) Warning: Monocle 3 trajectories require cluster partitions, which Seurat does not calculate. Please run ‘cluster_cells’ on your cell_data_set object > cds <- cluster_cells(cds_raw) > pr_graph_test_res <-…

Continue Reading Monocle3 differential expression failed when active.assay is not “RNA”

Using DESeq2 to analyse multi-variate design resulting in testing the wrong parameter

Enter the body of text here Hi, I am analysing a RNA Seq dataset coming from 3 independent cell isolates (isolate1, isolate2, isolate3), each given 3 different treatments (control, drug1, drug2). We are testing drug 1 against control in the first instance: We also observed that there is some variation…

Continue Reading Using DESeq2 to analyse multi-variate design resulting in testing the wrong parameter

How to generate a Venn diagram from edgeR DEG results?

How to generate a Venn diagram from edgeR DEG results? 0 Dear All, I have done DEG analysis of genes from 4 samples using edgeR. I want to plot them (number of DEGs) using a Venn diagram. I have results from pairwise comparisons as UP and DOWN gene list with…

Continue Reading How to generate a Venn diagram from edgeR DEG results?

Problem for analysing GSEA and ORA

Problem for analysing GSEA and ORA 0 I analysed two different GEO dataset individually and found DEG for each dataset .After that I compare the two dataset DEG and found some common DEG.Now i want to analyse GSEA and ORA for that common DEG.As far i know for analysing GSEA…

Continue Reading Problem for analysing GSEA and ORA

Identification of Hub Genes in Patients with Alzheimer Disease and Obs

Introduction Alzheimer’s disease (AD) ranks first among the common dementia type of the world. According to epidemiological investigation from the International Alzheimer’s disease association, about 45 million people has been suffered from AD, and the number is expected to increase to 131 million in 2050.1 Despite the widespread prevalence of…

Continue Reading Identification of Hub Genes in Patients with Alzheimer Disease and Obs

Re-analysis of multiple scRNAseq data as a combined new one

Re-analysis of multiple scRNAseq data as a combined new one 0 Hi all, There are so many regions in brain.The cortex itself is of many types.I found many single cell RNA seq studies which focus on different regions eg- gse … somatosensory cortex dropseq year 2015 gse….. primary visual cortex…

Continue Reading Re-analysis of multiple scRNAseq data as a combined new one

Sep 15 | GEO Data Mining and Submission – In Person

Event listing from University of Pittsburgh: Wednesday, September 15 from 10:00 AM to 4:00 PM This workshop will focus on Gene Expression Omnibus (GEO) repository. We will learn how to find high throughput gene expression studies – bulk RNA-seq/ scRNA-seq /microarray from GEO, retrieve gene expression count data and generate…

Continue Reading Sep 15 | GEO Data Mining and Submission – In Person

Sep 15 | GEO Data Mining and Submission – Online

Event listing from University of Pittsburgh: Wednesday, September 15 from 10:00 AM to 4:00 PM This workshop will focus on Gene Expression Omnibus (GEO) repository. We will learn how to find high throughput gene expression studies – bulk RNA-seq/ scRNA-seq /microarray from GEO, retrieve gene expression count data and generate…

Continue Reading Sep 15 | GEO Data Mining and Submission – Online

Gene with hundreds of counts very similar, is it possible to use it to get DEGs?

Gene with hundreds of counts very similar, is it possible to use it to get DEGs? 0 Hi, I’m starting in boinfomatics and I’m using Deseq2 for differential expression analysis. I have an RNAseq dataset, where one of the genes I intend to analyze has hundreds of counts ranging from…

Continue Reading Gene with hundreds of counts very similar, is it possible to use it to get DEGs?

GitHub – HaoWuLab-Bioinformatics/wu-group: ECSWalk

GitHub – HaoWuLab-Bioinformatics/wu-group: ECSWalk Files Permalink Failed to load latest commit information. Type Name Latest commit message Commit time ECSWalk: A carcinogenic driver module detection method based on a network model. Data used in this study: In our research, we used data processed by Rafsan Ahmed, et al. and Mark…

Continue Reading GitHub – HaoWuLab-Bioinformatics/wu-group: ECSWalk

create protein database for local Blast

create protein database for local Blast 1 Hello friends, how’re you? I wanted to know how I can do a blast against a database that I downloaded. I’ve downloaded a database of protein sequences from DEG (essential genes) and my idea is to do a blast against this database in…

Continue Reading create protein database for local Blast

List of human protein coding genes with given name (known function?)

List of human protein coding genes with given name (known function?) 2 Hello, To put it simply, I am doing differential expression analysis on human RNA-seq data and I want to focus my analysis of genes that are: 1) Protein coding, so no SNOR or MIR 2) Genes with a…

Continue Reading List of human protein coding genes with given name (known function?)

How to transform the deg gene list from seurat to a gene list input to clusterProfiler compareCluster ?

Sorry for lateness, I wanted to do something similar. This is what I did for reference: Using a Seurat generated gene list for input into ClusterProfiler to see the GO or KEGG terms per cluster. I’ll keep the meat and potatoes of the Seurat vignette in this tutorial: library(dplyr) library(Seurat)…

Continue Reading How to transform the deg gene list from seurat to a gene list input to clusterProfiler compareCluster ?

DEGs with uninsteresed sample files

DEGs with uninsteresed sample files 1 Hi all, I am working on some bulkRNAseq data from GEO and my sample description looks like this. I am interested in A vs B, A vs C and B vs C but I am not interested in sample D. My question is whether…

Continue Reading DEGs with uninsteresed sample files

How filter genes to construct co-expression network?

How filter genes to construct co-expression network? 1 Hi, I am interested to filter data for constructing co-expression network , Which parameter can i use to filter genes? As i know in WGCNA tutorial, it suggests not to use differential expressed genes(DEG) to filter genes. WGCNA Co-expreesion network DEG Filtering…

Continue Reading How filter genes to construct co-expression network?

FindMarkers for ClusterProfiler

FindMarkers for ClusterProfiler 1 Hi, I recently ran FindMarkers to compare DEG between two different clusters in a single-cell RNA-seq analysis This is my code: markers= FindMarkers(obj, ident.1=c(4), ident.2 = c(5)) head(markers) dim(markers) table(markers$avg_log2FC > 0) table(markers4v5$p_val_adj < 0.05 & markers$avg_log2FC > 0) I would like to run ClusterProfiler to…

Continue Reading FindMarkers for ClusterProfiler

Power calculation for microarray data

Power calculation for microarray data 0 I have an initial sample of 228 patients from a microarray study. Recently I have obtained a new set of labels specifying different condition types for only 69 out of the 228 patients. I wanted to run a DEG analysis on this set of…

Continue Reading Power calculation for microarray data

microarray miRNA expression data analysis

I wrote a script on how to analyze the microarray-based miRNA expression data. Here is my code: # general config baseDir <- ‘.’ annotfile <- ‘mirbase_genelist.tsv’ setwd(baseDir) options(scipen = 99) require(limma) # read in the data targets <- read.csv(“/media/mdrcubuntu/46B85615B8560439/microarray_text_files/targets.txt”, sep=””) # retain information about background via gIsWellAboveBG project <- read.maimages(targets,source=”agilent.median”,green.only…

Continue Reading microarray miRNA expression data analysis