logo
搜索
我的图书馆

发文章

发文工具

其他工具

留言交流
搜索
分享

人类参考基因组知识点(更新ing~)

 菌心说 2021-04-16

一、人类基因组有多大

chr size size2 1 chr1 248956422 249M 2 chr2 242193529 242M 3 chr3 198295559 198M 4 chr4 190214555 190M 5 chr5 181538259 182M 6 chr6 170805979 171M 7 chr7 159345973 159M 8 chrX 156040895 156M 9 chr8 145138636 145M 10 chr9 138394717 138M 11 chr11 135086622 135M 12 chr10 133797422 134M 13 chr12 133275309 133M 14 chr13 114364328 114M 15 chr14 107043718 107M 16 chr15 101991189 102M 17 chr16 90338345 90M 18 chr17 83257441 83M 19 chr18 80373285 80M 20 chr20 64444167 64M 21 chr19 58617616 59M 22 chrY 57227415 57M 23 chr22 50818468 51M 24 chr21 46709983 47M 25 SUM 3088269832 3088M #未考虑M线粒体，其长度较短，为16569，16Kbp,
  • 如上可看出染色体序号越靠前的，长度越大，范围在50M~250M之间；
  • 由于人为二倍体，所以基因组由60亿个碱基组成；
  • 参考基因组一般保存为纯文本格式，即直接记录“A”、“T”、“C”、“G”这样的 ASCII 码字符。
  • 而1个 ASCII 字符，大小是 1B，所以，如果按纯文本保存 30亿个字母(单链)，就是30亿字母 = 3,000,000,000 B = 3 GB。
from NCBI

二、奇怪的染色体name(chrUn,random,alt)

  • 同样以UCSC里的hg38版本为例
wget http://hgdownload.soe./goldenPath/hg38/bigZips/hg38.fa.gz
gunzip hg38.fa.gz
#提取染色体id
grep '^>' hg38.fa > chr.id
wc -l chr.id
#455 chr.id
head chr.id
####
>chr1
>chr10
>chr11
>chr11_KI270721v1_random
>chr12
>chr13
>chr14
>chr14_GL000009v2_random
>chr14_GL000225v1_random
>chr14_KI270722v1_random
  • 如上发现，序列并不是只有25条(22+X+Y+M)，加起来共有455条。其它特殊的序列可分为三类。

  • 在此之前需要简单了解由最初的测序read数据组装成基因组的染色体序列需要经历contigs与scaffolds两个过程，如下图所示。contigs是依靠read间的重叠拼接的序列(a few kbp long)，特点是不含有N碱基；scaffolds则主要依靠read pairs关系进一步拼接contigs，特点是会产生N碱基(a few hundred kbp)；最终由scaffolds拼接成染色体序列。


    read→contigs→scaffolds

    read→chromosomes

2.1 Unlocalized scaffolds(*****random)

  • a sequence found in an assembly that is associated with a specific chromosome but cannot be ordered or oriented on that chromosome.
  • 简单理解：知道这个scaffolds在哪条染色体上，但不知道其在染色体的具体位置及方向
  • format: chr{chromosome number orname}_{sequence_accession}v{sequence_version}_random
grep 'random' chr.id > chr.random wc -l chr.random #42 chr.random head chr.random ### >chr11_KI270721v1_random >chr14_GL000009v2_random >chr14_GL000225v1_random >chr14_KI270722v1_random >chr14_GL000194v1_random >chr14_KI270723v1_random >chr14_KI270724v1_random >chr14_KI270725v1_random >chr14_KI270726v1_random

2.2 Unplaced scaffolds(chrUn******)

  • a sequence found in an assembly that is not associated with any chromosome.
  • 简单理解：不知道这条scaffolds的所属染色体信息
  • format: chrUn_{sequence_accession}v{sequence_version}
grep 'chrUn' chr.id > chr.chrUn
wc -l chr.chrUn
#127 chr.chrUn
head chr.chrUn
###
>chrUn_KI270302v1
>chrUn_KI270304v1
>chrUn_KI270303v1
>chrUn_KI270305v1
>chrUn_KI270322v1
>chrUn_KI270320v1
>chrUn_KI270310v1
>chrUn_KI270316v1
>chrUn_KI270315v1
>chrUn_KI270312v1

2.3 Alternate loci scaffolds(*****alt)

  • a scaffold that provides an alternate representation of a locus found in the primary assembly. These sequences do not represent a complete chromosome sequence although there is no hard limit on the size of the alternate locus; currently these are less than 1 Mb. These could either be NOVEL patch sequences, added through patch releases, or present in the initial assembly release.
  • 简单理解：参考基因组存在的主要依据是人类99.9%的序列是一致的。但是会存在一些序列在不同人群中不一致。例如49%人群该基因组特定位置为序列A，而49%人群则为序列B，都是正常的。但拿其中一种作为参考基因组都可能不太合适，因此标记出Alternate loci scaffolds。
  • format: chr{chromosome number or name}_{sequence_accession}v{sequence_version}_alt
  • Alternate loci scaffolds为hg38版本基因组新添类型Sequence，此前hg19版本还没有。
grep 'alt' chr.id > chr.alt wc -l chr.alt #261 chr.alt head chr.alt ### >chr1_KI270762v1_alt >chr1_KI270766v1_alt >chr1_KI270760v1_alt >chr1_KI270765v1_alt >chr1_GL383518v1_alt >chr1_GL383519v1_alt >chr1_GL383520v2_alt >chr1_KI270764v1_alt >chr1_KI270763v1_alt >chr1_KI270759v1_alt

注意：以上具体的chromosome name均为ucsc的hg版本，与GRCh38略有差异，但基本也是这几种类型sequence

三、编码基因占比多少

  • 在30亿碱基基因组中，能够编码蛋白质的基因总长度只占总长度的5%，而其中转录本exon单元总长度只占总长度的1.5%；
  • 人类染色体共编码2w~3w个蛋白基因，分布于不同染色体中。平均长度有10Kbp长度左右，而实际上基因的长度分布十分广泛(from a few hundred bases to more than 2 million bases)
wget https://hgdownload.soe./goldenPath/hg38/bigZips/genes/hg38.refGene.gtf.gz
awk '{print$1, $10}' hg38.refGene.gtf |sort -k 2|uniq|grep -v alt | grep -v random | grep -v alt | grep -v fix| sort -k 1 > chr.gene
cut -d' ' -f 1 chr.gene | uniq -c
###
   1113 chr10
   1676 chr11
   1392 chr12
    632 chr13
    946 chr14
   1010 chr15
   1146 chr16
   1574 chr17
    434 chr18
   1791 chr19
   2832 chr1
    780 chr20
    414 chr21
    644 chr22
   1817 chr2
   1563 chr3
   1088 chr4
   1313 chr5
   1453 chr6
   1341 chr7
   1029 chr8
   1114 chr9
      1 chrM
   1157 chrX
    143 chrY

四、下载参考基因组

  • 目前常用的基因组版本为GRCh38/37，hg38/19，前者可通过NCBI/Ensembl下载，后者可通过UCSC网站下载。如下图所示GRCh38可认为等同于hg38，GRCh37可认为等同于hg19。


    human genome version
  • 以下载GRCh38/hg38为例，如下

4.1 NCBI

wget -c ftp://ftp.ncbi.nlm./refseq/H_sapiens/annotation/GRCh38_latest/refseq_identifiers/GRCh38_latest_genomic.fna.gz
NCBI

4.2 ensembl

wget -c http://ftp./pub/release-103/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.toplevel.fa.gz
ensembl

4.3 UCSC

wget -c http://hgdownload.soe./goldenPath/hg38/bigZips/hg38.fa.gz
UCSC

五、更新ing~

  • 如有错误欢迎指正；
  • 以及关于生信研究中人类参考基因组其它常见问题，也可评论区留言，让我们一起弄明白，加油~

    本站是提供个人知识管理的网络存储空间，所有内容均由用户发布，不代表本站观点。请注意甄别内容中的联系方式、诱导购买等信息，谨防诈骗。如发现有害或侵权内容，请点击一键举报。
    转藏 分享 献花（0）

    来自： 菌心说 > 《生物信息学，生信，统计，数据分析》

    举报/认领

    0条评论

    发表

    请遵守用户 评论公约

    类似文章 更多