Tag: GRCm38

A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain

Mouse breeding and husbandry All experimental procedures related to the use of mice were approved by the Institutional Animal Care and Use Committee of the AIBS, in accordance with NIH guidelines. Mice were housed in a room with temperature (21–22 °C) and humidity (40–51%) control within the vivarium of the AIBS…

Continue Reading A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain

RNAseq from P4 murine whole kidneys with Polycystic Kidney Disease

Snap frozen Postnatal day4 murine kidney tissues were crushed using micropestles in RLT Buffer and RNA was extracted from these tissues using QiaShredders followed by Qiagen RNAeasy Mini Kits according to manufacturer’s protocol. RNA concentration and integrity were determined by Aligent bioanalyzer. cDNA library preparation was carried out using NEB…

Continue Reading RNAseq from P4 murine whole kidneys with Polycystic Kidney Disease

biomaRt Ensembl version 83 (mouse)

biomaRt Ensembl version 83 (mouse) 0 Hello, I’m looking for an old Ensemb release (83) for GRCm38 build: library(biomaRt) mouse <- useEnsembl(biomart = “genes”, dataset=”mmusculus_gene_ensembl”, version=83) Error: Specified Ensembl version is not available. Use listEnsemblArchives() to view available versions. Version 83 is not displayed in the list of Ensembl archives…

Continue Reading biomaRt Ensembl version 83 (mouse)

couldn’t find matching transcriptome, returning non-ranged SummarizedExperiment AND unable to find an inherited method for function ‘seqinfo’ for signature ‘”SummarizedExperiment”‘

Dear Michael, I have not been able to run tximeta properly. I have read #38 but could not get any clue. The quant.sf files were generated by the latest nf-core RNA-seq pipeline (3.12.0), as the pipeline did not save the Salmon index, I generated it myself. Salmon used by nf-core…

Continue Reading couldn’t find matching transcriptome, returning non-ranged SummarizedExperiment AND unable to find an inherited method for function ‘seqinfo’ for signature ‘”SummarizedExperiment”‘

RNA sequencing and gene expression analysis in a Mouse model

Introduction Chronic obstructive pulmonary disease (COPD) is a condition that is characterized by persistent respiratory symptoms and airflow limitations that are not fully reversible. The severe complications of the disease may adversely affect its morbidity and mortality.1 According to World Health Organization (WHO) statistics, over 3 million people per year…

Continue Reading RNA sequencing and gene expression analysis in a Mouse model

Transcriptional linkage analysis with in vivo AAV-Perturb-seq

Experimental procedures Plasmid design and cloning AAV genome plasmids (Fig. 1a and Extended Data Figs. 1a,g,h and 5a) were based on Addgene plasmid 60231 (ref. 12). To achieve widespread transgene expression, the hSyn promoter was replaced by the ubiquitous CBh promoter (pAS088). For the triple-colour experiments (Extended Data Fig. 1a),…

Continue Reading Transcriptional linkage analysis with in vivo AAV-Perturb-seq

High-quality single-cell transcriptomics from ovarian histological sections during folliculogenesis

Introduction Single-cell RNA sequencing (RNA-seq) was first achieved by using a quantitative cDNA amplification method and applied to mouse oocytes (Kurimoto et al, 2006; Tang et al, 2009). It has since provided unprecedented opportunities for the study of cellular differentiations, states, and diseases in various biological fields, including developmental biology,…

Continue Reading High-quality single-cell transcriptomics from ovarian histological sections during folliculogenesis

TAPIS installation and usage

download gmap and install cd ~/software wget research-pub.gene.com/gmap/src/gmap-gsnap-2023-07-20.tar.gz tar xvzf gmap-gsnap-2023-07-20.tar.gz cd ~/software/gmap-gsnap-2023-07-20 ./configure –prefix=$HOME/jcbu/software/gmap make -j 20 make check make install build gmap index cd /home/jcbu/refGenome/gencode/mouse/GRCm38.p6_releaseM20/Sequence/gmapIndex nohup /home/jcbu/jcbu/software/gmap/bin/gmap_build -d GRCm38.gmap -D /home/jcbu/refGenome/gencode/mouse/GRCm38.p6_releaseM20/Sequence/gmapIndex \ /home/jcbu/refGenome/gencode/mouse/GRCm38.p6_releaseM20/Sequence/WholeGenomeFasta/GRCm38.p6.genome.fa & 3.alignPacBio.py cd /home/jcbu/YJ/PAIsoseq/ES/TAPIS/gmap/ nohup /usr/bin/python2 /home/jcbu/YJ/PAIsoseq/comp_bio-tapis-44cc05ebc78c/scripts/alignPacBio.m.py \ -p 20 -o Gm_Ctrl_rep1.alignPacBio \ /home/jcbu/refGenome/gencode/mouse/GRCm38.p6_releaseM20/Sequence/gmapIndex/GRCm38.gmap \…

Continue Reading TAPIS installation and usage

Identifying transcription factors for a list of genes

Identifying transcription factors for a list of genes 1 Is there a way to find transcription factors for mouse genes if Ensemble IDs are known? In an article I’m trying to recreate, transcription factor features were retrieved from GTRD database, which is no longer in service. Here is the study…

Continue Reading Identifying transcription factors for a list of genes

scrnaseq: Introduction

Introduction nf-core/scrnaseq is a bioinformatics best-practice analysis pipeline for processing 10x Genomics single-cell RNA-seq data. This is a community effort in building a pipeline capable to support: Alevin-Fry + AlevinQC STARSolo Kallisto + BUStools Cellranger UniverSC Documentation The nf-core/scrnaseq pipeline comes with documentation about the pipeline usage, parameters and output….

Continue Reading scrnaseq: Introduction

Non-cell-autonomous cancer progression from chromosomal instability

Cell culture IMR90, 4T1, CT26, RAW264.7 and B16F10 cell lines were purchased from the American Type Culture Collection and cultured in MEM (IMR90), DMEM (B16F10, RAW264.7) or RPMI (4T1, IMR90, CT26) supplemented with 10% FBS in the presence of penicillin (50 U ml−1) and streptomycin (50 μg ml−1). All cells were found to be…

Continue Reading Non-cell-autonomous cancer progression from chromosomal instability

High discordant alignment with HISAT2 Mus musculus

High discordant alignment with HISAT2 Mus musculus 0 I ran HISAT2 on my reads post adaptor removal, and received very high discordant alignments (>80%) with low concordant alignments. 81987410 reads; of these: 81987410 (100.00%) were paired; of these: 81584433 (99.51%) aligned concordantly 0 times 92014 (0.11%) aligned concordantly exactly 1…

Continue Reading High discordant alignment with HISAT2 Mus musculus

bedtools getfasta chromosome not found

bedtools getfasta chromosome not found 1 I am having the same problem as this post: Trouble with bedtools getfasta When I run the following command, I get this warning for every chromosome and the output fasta is empty: bedtools getfasta -fi GRCm38.p6.genome.fasta -bed e2f5_R1.seacr.peaks.stringent.bed -fo e2f5_R1.seacr.peaks.fasta WARNING. chromosome (1) was…

Continue Reading bedtools getfasta chromosome not found

reference annotation for the human and mouse genomes in 2023

D942–D949 Nucleic Acids Research, 2023, Vol. 51, Database issue Published online 24 November 2022 doi.org/10.1093/nar/gkac1071 GENCODE: reference annotation for the human and mouse genomes in 2023 Adam Frankish 1,* , Sı́lvia Carbonell-Sala2 , Mark Diekhans 3 , Irwin Jungreis 4,5 , Jane E. Loveland 1 , Jonathan M. Mudge1 ,…

Continue Reading reference annotation for the human and mouse genomes in 2023

bbsplit running slow or out of memory?

bbsplit running slow or out of memory? 1 Hello, I have Illumina fastq files from some RNA-seq, ATAC-seq and WES that originated as PDX samples. I am looking to filter out contaminating mouse reads from the human reads in these datasets. I have used Xenome before but wanted to try…

Continue Reading bbsplit running slow or out of memory?

BMPER is a marker of adipose progenitors and adipocytes and a positive modulator of adipogenesis

Human subject characteristics Visceral (greater omental) adipose tissue samples were obtained in consented patients for our University of Utah Institutional Review Board (IRB)-approved obesity biorepository at the time of laparoscopic bariatric surgery (obesity group) or laparoscopic abdominal surgery for non-cancerous, non-infected general surgery indications (lean group). Patients were fasted for…

Continue Reading BMPER is a marker of adipose progenitors and adipocytes and a positive modulator of adipogenesis

Whole genome analysis for 163 gRNAs in Cas9-edited mice reveals minimal off-target activity

Allele design and guide selection For multi-exon genes, a critical region (one or more exons) was identified as shared among all annotated full-length transcripts whose removal was predicted to result in a frame-shift mutation and introduction of premature stop codon greater than 50-nt from the final splice junction increasing the…

Continue Reading Whole genome analysis for 163 gRNAs in Cas9-edited mice reveals minimal off-target activity

Assembly Table.

Assembly Table. A. mellifera (Apr 2011 Amel_4.5/amel5) A. carolinensis (May 2010 AnoCar2.0/anoCar2) A. thaliana (Feb 2011 TAIR10/araTha1) B. taurus (Aug 2006 Btau_3.1/bosTau3) B. taurus (Nov 2014 Bos_taurus_UMD_3.1.1/bosTau8) C. familiaris (May 2005 CanFam2.0/canFam2) C. familiaris (Sep 2011 CanFam3.1/canFam3) C. porcellus (Feb 2008 Cavpor3.0/cavPor3) C. elegans (Oct 2010 WBcel215/ce10) C. elegans (Feb…

Continue Reading Assembly Table.

Prior metabolite extraction fully preserves RNAseq quality and enables integrative multi-‘omics analysis of the liver metabolic response to viral infection

Introduction The metabolome is an incredibly diverse collection of small molecules (<1,500 Da) in biological systems involved in virtually every cellular process, including cellular energy production, macromolecule synthesis, epigenetic modifications, cell signalling and more (for recent reviews see [Citation1–6]). It responds rapidly (in seconds) to both internal (signalling, allostery) and external…

Continue Reading Prior metabolite extraction fully preserves RNAseq quality and enables integrative multi-‘omics analysis of the liver metabolic response to viral infection

Errors after running jobscript for cellranger vdj pipeline

Errors after running jobscript for cellranger vdj pipeline 0 Hello, Below is pasted my jobscript to be submitted to run cellranger vdj pipeline on my fasta sample files. 1 #!/bin/bash 2 #PBS -N RED_aregvdj 3 #PBS -q serial 4 #PBS -l nodes=1:ppn=24 5 #PBS -l walltime=03:00:00 6 #PBS -l pmem=3300mb…

Continue Reading Errors after running jobscript for cellranger vdj pipeline

RHOJ controls EMT-associated resistance to chemotherapy

Mice Rosa26-YFP48, Rosa-tDTomato49, K14creER, Lgr5creER (ref. 50), KrasLSLG12D (ref. 51) and Trp53fl/fl (ref. 52) mice were obtained from the NCI mouse repository and Jackson Laboratories. Rhojfl/fl mice53 were a gift from A. Uemura. All mice used in this study were composed of males and females with mixed genetic background. Mouse…

Continue Reading RHOJ controls EMT-associated resistance to chemotherapy

Bipotent transitional liver progenitor cells contribute to liver regeneration

Mice All mice experiments were performed in accordance with the guidelines of the Institutional Animal Care and Use Committee (IACUC) at the Center for Excellence in Molecular Cell Science, Shanghai Institutes of Biological Sciences, Chinese Academy of Science. The approved animal protocol number is SIBCB-S374-1702-001-C1. The CK19-CreER, HNF4α-DreER, R26-tdTomato (R26-tdT),…

Continue Reading Bipotent transitional liver progenitor cells contribute to liver regeneration

Transcript to gene level count for DEseq(2) use- Salmon/Sailfish/Kallisto etc.

I am working on analysis for an RNASeq experiment with the end goal of using DESeq2 to generate a list of differentially expressed genes. We have 4 biological replicates for 4 conditions (differing genotypes in mice). We have 50 bp single ended Illumina read sets, on total RNA (they had…

Continue Reading Transcript to gene level count for DEseq(2) use- Salmon/Sailfish/Kallisto etc.

Fumarate induces vesicular release of mtDNA to drive innate immunity

Mice Mice were of mixed genetic background C57BL/6 and 129/SvJ. Mice were bred and maintained under specific-pathogen-free conditions at the Breeding Unit (BRU) at the CRUK Cambridge Institute. Fh1fl/fl (ref. 4) and R26-Cre-ERT2 (ref. 5) mice were gifts from E. Gottlieb and D. Winton, respectively. Experimental mice were homozygous for…

Continue Reading Fumarate induces vesicular release of mtDNA to drive innate immunity

Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut

Animal studies Animal care and study protocols were approved by the AAALAC-accredited Institutional Animal Care and Use Committee of the College of Agricultural Life Sciences at the University of Wisconsin-Madison (UW-Madison). All experiments with mice were performed under protocols approved by the UW-Madison Animal Care and Use Committee (Protocol number…

Continue Reading Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut

Nociceptor neurons affect cancer immunosurveillance

Secondary use of biopsies as research specimens The ten melanoma samples used in this study were collected by Sanford Health and classified by a board-certified pathologist. Their secondary use as research specimens (fully de-identified formalin-fixed, paraffin-embedded (FFPE) blocks) was approved under Sanford Health IRB protocol 640 (titled ‘Understanding and improving…

Continue Reading Nociceptor neurons affect cancer immunosurveillance

In vivo transomic analyses of glucose-responsive metabolism in skeletal muscle reveal core differences between the healthy and obese states

Animals and sample preparation Animal experiments were performed as previously described12. C57BL/6J WT mice or ob/ob mice at ten weeks of age were purchased from Japan SLC Inc. (Shizuoka, Japan). The phenotypic data of the mice are summarized in Table S1. Animal experiments were approved by the animal ethics committee…

Continue Reading In vivo transomic analyses of glucose-responsive metabolism in skeletal muscle reveal core differences between the healthy and obese states

Hisat2 – stringtie – deseq2 pipeline for bulk RNA seq

Software official website : Hisat2: Manual | HISAT2 StringTie:StringTie article :Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown | Nature Protocols It is recommended to watch the nanny level tutorial : 1. RNA-seq : Hisat2+Stringtie+DESeq2 – Hengnuo Xinzhi 2. RNA-seq use hisat2、stringtie、DESeq2 analysis – Simple books Basic usage…

Continue Reading Hisat2 – stringtie – deseq2 pipeline for bulk RNA seq

how to build index for cdna?

Hello, I can build index for Mus_musculus.GRCm38.dna_sm.toplevel.fa, but when build for Mus_musculus.GRCm38.cdna.all.fa, there is a bug: “rsem-extract-reference-transcripts Mus_musculus.GRCm38.cdna.all.fa 0 Mus_musculus.GRCm38.cdna.all.fa.gtf None 0 Mus_musculus.GRCm38.cdna.all” failed! Plase check if you provide correct parameters/options for the pipeline! Traceback (most recent call last): File “../indrops.py”, line 1770, in project.build_transcriptome(args.genome_fasta_gz, args.ensembl_gtf_gz, mode=args.mode) File “../indrops.py”, line…

Continue Reading how to build index for cdna?

deeptools plotHeatMap – Convert bed files to gene list?

I might suggest limiting your search to genes: gtf2bed < Mus_musculus.GRCm38.102.gtf | grep -w “gene” > sorted-mm10.genes.bed But that’s up to you. Otherwise, I think you may also get transcripts/exons, which may be more than you want. Again, up to you. If hnf4a-ko-downreg-clusters.bed is the file containing peaks, as described…

Continue Reading deeptools plotHeatMap – Convert bed files to gene list?

Htseq is giving me 0 counts using the GFF3 of miRBase

Hello! I am trying to annotate a miRNA-seq so that it gives me mature miRNAs where I already have 5p and 3p. For this, I have used the index mm10.fa and the miRBase mmu.gff3. I have aligned with HISAT2 and am trying to count with HTSeq, however I get 0…

Continue Reading Htseq is giving me 0 counts using the GFF3 of miRBase

Transcriptional kinetics and molecular functions of long noncoding RNAs

Ethical compliance The research carried out in this study has been approved by the Swedish Board of Agriculture, Jordbruksverket: N343/12. Cell culture Mouse primary fibroblasts were derived from adult (>10 weeks old) CAST/EiJ × C57BL/6J or C57BL/6J × CAST/EiJ mice by skinning, mincing and culturing tail explants (for at least 10 d) in DMEM high…

Continue Reading Transcriptional kinetics and molecular functions of long noncoding RNAs

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

python – How can I fix the dash bio error: devtools cannot load source map dashbio@1.0.1 bundle.js.map?

I am implementing a website in Python with Django framework and using django-plotly-dash to display data. I am trying to use dash_bio’s IGV feature to display some chromosome data, but when I attempt to call the functionality, I receive the following errors and the callback that returns ‘dashbio.igv’ is unable…

Continue Reading python – How can I fix the dash bio error: devtools cannot load source map dashbio@1.0.1 bundle.js.map?

Why are there more lincRNA genes than transcripts in Mus_musculus.GRCm38.100.gtf?

Why are there more lincRNA genes than transcripts in Mus_musculus.GRCm38.100.gtf? 1 Hi All, There are 52446 annotated genes (ENSMUSG IDs) and 142,699 transcripts (ENSMUST IDs) in Mus_musculus.GRCm38.100.gtf. It makes sense that there are WAY more transcripts than genes. My question, however, is – why are there more genes than transcripts…

Continue Reading Why are there more lincRNA genes than transcripts in Mus_musculus.GRCm38.100.gtf?

Characterization and expression of DNA sequences encoding the growth hormone gene in African Pygmy Mouse (Mus minutoides)

Abstract We determined the nucleotide sequence of the growth hormone (Gh) gene in Mus minutoides, one of the smallest mammals, where was predicted to be distinct in the functional regions between M. minutoides and Mus musculus. To investigate the evolutionary characteristics of Gh in M. minutoides, we constructed a phylogenetic…

Continue Reading Characterization and expression of DNA sequences encoding the growth hormone gene in African Pygmy Mouse (Mus minutoides)

Good or bad idea to annotate with newer genome build?

Good or bad idea to annotate with newer genome build? 0 Hello Biostars Community, If the probe sequences on a mouse methylation beadchip array is based on the (mm10) GRCm38 genome build, would it be a bad idea to annotate the probes with the latest genome build (mm39) GRCm39? I…

Continue Reading Good or bad idea to annotate with newer genome build?

STAR align multiple files

STAR align multiple files 1 Hi everybody, I am doing alignment to 36 PE samples using star. to make it little bit easy task I wrote a bash loop to align them all with the same command. here is my loop: for i in $(ls raw_data); do STAR –genomeDir index.150…

Continue Reading STAR align multiple files

IOError [errno 2] No such file or directory: ‘-o’

TopHat error: IOError [errno 2] No such file or directory: ‘-o’ 2 Hello everyone. I’m now running tophat in our server. First, I just simply tried “tophat -p 8 -G <gtf_file> <ref_genome> <fa1><fa2>” and it worked. Then I wrote a for loop scripts but It reported error: [2019-04-03 10:49:16] Beginning…

Continue Reading IOError [errno 2] No such file or directory: ‘-o’