学习资料

基因组不同版本介绍

下载人类基因组信息有三大去处，ensembl、RefSeq（NCBI）和UCSC。ensembl、RefSeq（NCBI）和UCSC下载的同一版本的DNA序列是一样的，目前比较常用的有hg19、hg38、GRCh37、GRCh38，其中hg系列是UCSC的叫法，GRCh系列是NCBI和ENSEMBL的叫法。这些都属于genome builds，也就是基因序列，同一版本的序列是一样的（其中hg19对应GRCh37，hg38对应GRCh38），比如hg38和GRCh38的序列就是一样的。同一系列，序号大的是序号小的升级版，比如GRCh38是GRCh37的更新版。

ensembl、RefSeq（NCBI）和UCSC下载的同一版本的DNA序列虽然是一样的，但是它们的注释是不同的，更新频率也不一样。NCBI 的注释是“ refseq”数据集(UCSC 的”refGene” track)， UCSC 和 ensemble 注释都将其作为其中一个子集。而UCSC 的注释比较混乱。同样ID的基因会出现在不同链或不同染色体位置上。Ensembl的注释通常比UCSC更多(所以噪音更多一点) ，但是它们的ID处理的比较好，不会像UCSC的那么混乱，因此，它们的 ID 通常更容易转换成别的(例如，基因名称、GO和pathway)。Ensembl还经常更新它的注释，更新一次作为一个版本发布，因此可以方便地告诉别人你使用的版本，增强分析的可重复性。

Ensembl的注释噪音多怎么理解呢？有哪些潜在的影响呢？
以ensembl release-74版本为例，尽管 其包含63,677个基因注释条目，但只有22,810个条目(大约三分之一)是蛋白编码基因。有17,057个条目是不同类型的 rna，包括 rRNA (566)、snoRNA (1549) 、snRNA (2067)、miRNA (3361) 、misc_rna (2174)和 lincRNA (7340)。ensemble还包含了15,583个假基因。对于大多数RNA-Seq 测序，一般只有富集和测序mRNAs，一般不会比对到诸如 miRNAs 或 lincRNAs 之类的RNA上。R74基因组包含819个经过逆转录产生的转录本，这些转录本是通过逆转录产生的，随后又重新整合到基因组中，这些转录本通常是表达不活跃的。在这种情况下，mRNA产生的reads，尤其是junction reads，可能比对到甚至唯一比对到这种“处理过的转录本”上。因此，该基因的mRNA的真实表达量可能被低估了。使用更全面的注释数据库的另一个缺点是p值，因为原始p值的需要多重校验。而多重校验后的p值的大小主要取决于模型中基因的数量。如果感兴趣的基因在不同的注释中定义不一致，建议使用不同的基因注释模型分析 RNA-Seq 的数据。

总的来说，可以遵循这么一个原则：当进行强调可重复性和稳健的基因表达估计的研究时（比如使用RNAseq来作为芯片分析的替代的时候），优先选较为简单的基因组注释，如 RefGene。但较为简单的基因组注释可能不利于识别或描述新的转录或调控机制。因此，当进行更具探索性的研究时，更全面的注释更优，比如选择Ensembl。而UCSC，emmmm不太建议使用。

Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna.toplevel.fa.gz 615M
Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna_rm.toplevel.fa.gz 123M
Zea_mays.Zm-B73-REFERENCE-NAM-5.0.dna_sm.toplevel.fa.gz 641M

   <species>.<assembly>.<sequence type>.<id type>.<id>.fa.gz

When you look at all the genome files available from Ensembl. You are presented with a bunch of options. Which one is the best to use/download?

You have a combination of choices.

注意：

unmasked、masked、soft-masked的基因组序列解释

从ensembl下载基因组fasta文件时，其中sequence type包括：unmasked genomic DNA sequences、masked genomic DNA和soft-masked genomic DNA。

Haplotypes and patches

Ensembl中的所有基因组组件都是单倍体，对于大多数物种而言，只有一条通过基因组的路径。

小鼠基因组

参考