小鼠参考基因组id转换gtf文件chb注释官网下载相应的基因中注释文件gse155802 getmatrixgenecode 参考基因组下载序列比对 hg19 索引文件mapping referece

本文主要是介绍小鼠参考基因组id转换gtf文件chb注释官网下载相应的基因中注释文件gse155802 getmatrixgenecode 参考基因组下载序列比对 hg19 索引文件mapping referece,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

 参考基因组在很多网站都可以下载到

Index of /pub/release-93/gtf/mus_musculus (ensembl.org)http://ftp.ensembl.org/pub/release-93/gtf/mus_musculus/ 我想要研究的基因组

我怎么知道是这个的呢 ,因为这篇文献里面提到

References - 3.0.0 (November 19, 2018)

  • Human reference, GRCh38 (Ensembl 93)
  • Human reference, hg19 (Ensembl 87)https://support.10xgenomics.com/single-cell-gene-expression/software/release-notes/build#hg19_3.0.0
  • Mouse reference, mm10 (Ensembl 93)

所以我就在cellranger官网上找到了这个Build Notes for Reference Packages -Software -Single Cell Gene Expression -Official 10x Genomics Support

点进去之后,出现代码

Mouse reference, mm10 (includes mouse V(D)J genes)
wget ftp://ftp.ensembl.org/pub/release-93/fasta/mus_musculus/dna/Mus_musculus.GRCm38.dna.primary_assembly.fa.gz
gunzip Mus_musculus.GRCm38.dna.primary_assembly.fa.gzwget ftp://ftp.ensembl.org/pub/release-93/gtf/mus_musculus/Mus_musculus.GRCm38.93.gtf.gz
gunzip Mus_musculus.GRCm38.93.gtf.gz

根据代码下载的地址  就可以找到自己要复现的基因组地址了!这里我只需要gtf格式就行啦

Index of /pub/release-93/gtf/mus_musculus (ensembl.org)

手动下载之后

如何处理gtf文件呢?在R语言中读取GTF文件的最好方法 · 大专栏 (dazhuanlan.com)

链接里给了好多方法  哈哈哈

我就用第一个吧 ,开始!

####参考基因组id转换gtf文件 https://www.dazhuanlan.com/insafe/topics/975998
#BiocManager::install("rtracklayer")
library(rtracklayer)#https://support.10xgenomics.com/single-cell-gene-expression/software/release-notes/build
#把Mus_musculus.GRCm38.93.gtf.gz下载下来 为什么下载这个呢?http://ftp.ensembl.org/pub/release-93/gtf/mus_musculus/
gtf_df=rtracklayer::import('Mus_musculus.GRCm38.93.gtf.gz')
getwd()

更改seurat行名 seurat基因名

如何给Seurat对象的基因重命名?【基因名转换】2022-07-27 - 简书 (jianshu.com)

#创建函数 改名字
RenameGenesSeurat <- function(obj , 
                              newnames ) { 
  # Replace gene names in different slots of a Seurat object. Run this before integration. Run this before integration. 
  # It only changes obj@assays$RNA@counts, @data and @scale.data.
  print("Run this before integration. It only changes obj@assays$RNA@counts, @data and @scale.data.")
  RNA <- obj@assays$RNA
  
  if (nrow(RNA) == length(newnames)) {
    if (length(RNA@counts)) RNA@counts@Dimnames[[1]]            <- newnames
    if (length(RNA@data)) RNA@data@Dimnames[[1]]                <- newnames
    if (length(RNA@scale.data)) RNA@scale.data@Dimnames[[1]]    <- newnames
  } else {"Unequal gene sets: nrow(RNA) != nrow(newnames)"}
  obj@assays$RNA <- RNA
  return(obj)
}

#教程地址
#https://cloud.tencent.com/developer/article/1697249#https://bioconductor.org/packages/release/data/experiment/vignettes/scRNAseq/inst/doc/scRNAseq.html
#https://mp.weixin.qq.com/mp/appmsgalbum?__biz=MzI1Njk4ODE0MQ==&action=getalbum&album_id=1326587538303434752&scene=173&from_msgid=2247484689&from_itemidx=1&count=3&nolastread=1#wechat_redirect
rm(list = ls()) 
Sys.setenv(R_MAX_NUM_DLLS=999)
options(stringsAsFactors = F)##########三对三的数据
1#准备原始分析数据  先手动下载  去浏览器下载文件到自己的文件夹下,然后解压
#https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi  
# 3个样本,共9个文件。需要分到9个文件夹里,并且重命名
getwd()
path="G:/silicosis/geo/GSE155802_scRNA-seq_AT2_sorted_mice_lung_chb/" #空间转录组
dir.create(path,recursive = T)
setwd(path)
getwd()??list.files
dirs_and_files=list.files(path = "./",all.files = T,pattern = ".csv",include.dirs = F,recursive = T)fs=grep(x=dirs_and_files,"GSE155802_RAW/GSM.+/",value = T)
fs
getwd()

这篇关于小鼠参考基因组id转换gtf文件chb注释官网下载相应的基因中注释文件gse155802 getmatrixgenecode 参考基因组下载序列比对 hg19 索引文件mapping referece的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/209268

相关文章

mybatis执行insert返回id实现详解

《mybatis执行insert返回id实现详解》MyBatis插入操作默认返回受影响行数,需通过useGeneratedKeys+keyProperty或selectKey获取主键ID,确保主键为自... 目录 两种方式获取自增 ID:1. ​​useGeneratedKeys+keyProperty(推

浅谈mysql的not exists走不走索引

《浅谈mysql的notexists走不走索引》在MySQL中,​NOTEXISTS子句是否使用索引取决于子查询中关联字段是否建立了合适的索引,下面就来介绍一下mysql的notexists走不走索... 在mysql中,​NOT EXISTS子句是否使用索引取决于子查询中关联字段是否建立了合适的索引。以下

关于集合与数组转换实现方法

《关于集合与数组转换实现方法》:本文主要介绍关于集合与数组转换实现方法,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1、Arrays.asList()1.1、方法作用1.2、内部实现1.3、修改元素的影响1.4、注意事项2、list.toArray()2.1、方

使用Python实现可恢复式多线程下载器

《使用Python实现可恢复式多线程下载器》在数字时代,大文件下载已成为日常操作,本文将手把手教你用Python打造专业级下载器,实现断点续传,多线程加速,速度限制等功能,感兴趣的小伙伴可以了解下... 目录一、智能续传:从崩溃边缘抢救进度二、多线程加速:榨干网络带宽三、速度控制:做网络的好邻居四、终端交互

Python中注释使用方法举例详解

《Python中注释使用方法举例详解》在Python编程语言中注释是必不可少的一部分,它有助于提高代码的可读性和维护性,:本文主要介绍Python中注释使用方法的相关资料,需要的朋友可以参考下... 目录一、前言二、什么是注释?示例:三、单行注释语法:以 China编程# 开头,后面的内容为注释内容示例:示例:四

MySQL之InnoDB存储引擎中的索引用法及说明

《MySQL之InnoDB存储引擎中的索引用法及说明》:本文主要介绍MySQL之InnoDB存储引擎中的索引用法及说明,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐... 目录1、背景2、准备3、正篇【1】存储用户记录的数据页【2】存储目录项记录的数据页【3】聚簇索引【4】二

全面解析MySQL索引长度限制问题与解决方案

《全面解析MySQL索引长度限制问题与解决方案》MySQL对索引长度设限是为了保持高效的数据检索性能,这个限制不是MySQL的缺陷,而是数据库设计中的权衡结果,下面我们就来看看如何解决这一问题吧... 目录引言:为什么会有索引键长度问题?一、问题根源深度解析mysql索引长度限制原理实际场景示例二、五大解决

MySQL中的索引结构和分类实战案例详解

《MySQL中的索引结构和分类实战案例详解》本文详解MySQL索引结构与分类,涵盖B树、B+树、哈希及全文索引,分析其原理与优劣势,并结合实战案例探讨创建、管理及优化技巧,助力提升查询性能,感兴趣的朋... 目录一、索引概述1.1 索引的定义与作用1.2 索引的基本原理二、索引结构详解2.1 B树索引2.2

利用Python脚本实现批量将图片转换为WebP格式

《利用Python脚本实现批量将图片转换为WebP格式》Python语言的简洁语法和库支持使其成为图像处理的理想选择,本文将介绍如何利用Python实现批量将图片转换为WebP格式的脚本,WebP作为... 目录简介1. python在图像处理中的应用2. WebP格式的原理和优势2.1 WebP格式与传统

python3如何找到字典的下标index、获取list中指定元素的位置索引

《python3如何找到字典的下标index、获取list中指定元素的位置索引》:本文主要介绍python3如何找到字典的下标index、获取list中指定元素的位置索引问题,具有很好的参考价值,... 目录enumerate()找到字典的下标 index获取list中指定元素的位置索引总结enumerat