CN101928775A - 一种细胞染色体分析方法 - Google Patents

一种细胞染色体分析方法 Download PDF

Info

Publication number
CN101928775A
CN101928775A CN201010252849XA CN201010252849A CN101928775A CN 101928775 A CN101928775 A CN 101928775A CN 201010252849X A CN201010252849X A CN 201010252849XA CN 201010252849 A CN201010252849 A CN 201010252849A CN 101928775 A CN101928775 A CN 101928775A
Authority
CN
China
Prior art keywords
chrn
karyomit
cell lines
sample
dna
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201010252849XA
Other languages
English (en)
Inventor
张秀清
玄兆伶
陈芳
蒋馥蔓
林静蓉
李培培
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BGI Shenzhen Research Institute
BGI Shenzhen Co Ltd
Original Assignee
BGI Shenzhen Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by BGI Shenzhen Co Ltd filed Critical BGI Shenzhen Co Ltd
Priority to CN201010252849XA priority Critical patent/CN101928775A/zh
Publication of CN101928775A publication Critical patent/CN101928775A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明涉及一种细胞染色体的分析方法,特别是涉及通过测序方法分析羊水细胞染色体数量与标准细胞是否存在差异的方法。

Description

一种细胞染色体分析方法 
技术领域
本发明涉及一种细胞染色体的分析方法,特别是涉及通过测序方法分析羊水细胞染色体的方法。 
背景技术
随着测序技术的不断发展,测序技术在细胞学研究,特别是在细胞染色体分析中得到了越来越广泛的应用(Y.M.Dennis Lo,et al.Quantitative Analysis of Fetal DNA in Maternal Plasma and Serum:Implications for Noninvasive Prenatal Diagnosis.Am.J.Hum.Genet.62:768-775,1998)。本发明人采用测序技术直接对羊水细胞,包括未经培养的羊水细胞进行染色体分析,成功实现了对21、18、13号染色体及其他染色体的分析,进一步拓宽了测序技术在细胞学研究领域中应用。 
发明内容
本发明的一个方面,涉及一种用测序方法分析受试者细胞染色体信息的方法,其包括以下步骤: 
a.将获自细胞的基因组DNA分子随机打断,得到一定大小的DNA片段并测序; 
b.将步骤a中测定的DNA序列与人类基因组参考序列进行严格比对,得到所测DNA序列定位于具体染色体的信息; 
c.对应于特定染色体N,确定所测DNA序列中定位于该染色体上的唯一位置的序列的数量总和,进而得出对应于该染色体N的ChrN%,即所测DNA序列中定位于该染色体N上的唯一位置的序列的数量总和(S1)占所测DNA序列中定位于所有染色体上的唯一位置 的序列的数量总和(S2)的比例:ChrN%=S1/S2; 
d.将该染色体的ChrN%与来自标准细胞中对应染色体的ChrN%进行比较,确定所述细胞的该条染色体与标准细胞是否存在差异。 
本发明中,所述的细胞可以是例如羊水细胞,其中所述的羊水细胞可以为未培养的羊水细胞或培养过的羊水细胞。在本发明的一个实施方案中,为避免培养羊水细胞,羊水细胞为未培养的羊水细胞。 
本发明中,细胞基因组DNA的获取可以采用盐析法、柱层析法、SDS法等常规DNA提取方法,优选采用柱层析法。所谓柱层析法,是指羊水细胞或组织经过细胞裂解液和蛋白酶K的作用后露出裸露的DNA分子,其经过能与带负电的DNA分子结合的硅胶膜柱时,体系中的基因组DNA被可逆吸附,经漂洗液清洗除去蛋白质、脂质等杂质后,用纯化液洗脱获得羊水细胞DNA(具体原理和方法参见Qiagen公司货号为56304的产品说明书和Tiangen公司货号为DP316的产品说明书)。 
本发明中,所述DNA分子的随机打断处理可以采用酶切、雾化、超声、或者HydroShear法。优选地,采用HydroShear法(当含有DNA的溶液通过较小面积的通道时,流体加速,产生的力使DNA突然断裂,流速和通道大小决定DNA片段的大小,具体原理和方法请参见Life Sciences Wiki公司的HydroShear说明书),将DNA分子打断为比较集中的一定大小的片段,通常经过随机打断的主带分布在200~300bp范围。 
本发明中,所采用的测序方法可以为第二代测序方法Illumina/Solexa或ABI/SOLiD。在本发明的一个实施方案中,所述的测序方法为Illumina/Solexa,测序得到的结果为35bp大小的片段。 
当待测的DNA分子来自多个受试样本时,每个样本可以被加上不同的标签序列index,以用于在测序过程中进行区分(Micah Hamady,Jeffrey J Walker,J Kirk Harris et al.Error-correcting barcoded primers forpyrosequencing hundreds of samples in multiplex.Nature Methods,2008,March,Vol.5No.3)。 
本发明中,所述的人类基因组参考序列是人类基因组序列经过屏蔽掉重 复序列后所得到的参考序列,例如NCBI数据库中最新版本的人类基因组参考序列。在本发明的一个实施方案中,所述人类基因组序列是NCBI数据库中版本36(NCBI Build 36)的人类基因组参考序列。 
在本发明中,其中所述与人类基因组参考序列进行严格比对,是指采用的比对方法为不容错的、定位于人类基因组参考序列上的唯一位置的比对。在本发明的一个实施方案中,采用的是Eland比对软件(Illumina提供的软件包),采用的方法为完全不容错比对。 
本发明中,当所述的DNA序列是能够定位于人类基因组参考序列的唯一位置上的DNA序列时,将其定义为唯一序列,以Unique reads表示。在本发明中,为了避免重复序列的干扰,需要去除那些定位于人类基因组参考序列中的串联重复及转座重复位置的DNA序列,只统计那些可以定位到基因组唯一位置的DNA序列,即唯一序列。通常,在所有测序的DNA序列中,大约有四分之一到三分之一的DNA序列可以定位于基因组唯一位置,即唯一序列。而这个统计所得的唯一序列的数量代表分布在基因组染色体上的DNA序列的情况。 
因此,唯一序列能够将来自羊水细胞的DNA分子经打断并测序后的各DNA序列定位于特定染色体。而ChrN%则是将不同染色体上得到的唯一序列进行标准化后得到的数值,该数值只与特定染色体大小有关,而与测序数据量无关,从而可以用ChrN%对46条染色体各自的情况进行分析。因而,所述ChrN%是本发明中进行染色体分析的基础数值。 
在本发明中,所述细胞样本中特定染色体数量与标准细胞是否存在差异可以通过绘制箱式图方法确定,其中将细胞样本中ChrN%位于超出四分位数差1.5倍-3倍或者3倍以上的四分位数差之间的异常值对应的样本,认定为其染色体数量与标准细胞存在差异。 
在本发明中,所述细胞的特定染色体与标准细胞是否存在差异还可以用反映受检细胞样本ChrN%与标准细胞样本ChrN%的偏离程度z score_ChrN来确定的。 
具体的,所述z score_ChrN=(受检样本特定染色体的ChrN%-特定染色体ChrN%平均值(mean_ChrN%))/ChrN%标准差 (S.D._ChrN%)。 
z score_ChrN极大或极小,说明该待检细胞样品中染色体数量偏离正常值非常显著,当达到一定显著性水平时,可认为其染色体数量与标准细胞存在明显差异。 
本发明中,所述来自标准细胞样本的特定染色体ChrN%平均值(mean_ChrN%)可根据来自例如至少10个,例如至少20个,例如至少30个,例如至少50个,例如至少100个标准细胞样本的该染色体ChrN%确定。 
在本发明中,标准细胞样本是指细胞中染色体数量为单倍体二倍的人类细胞样本。正常男性细胞中,有44条常染色体和两条不同的性染色体,记为(46,XY),而正常女性细胞中,有44条常染色体和两条一样的性染色体,记为(46,XX)。 
本发明中,所述来自标准细胞样本的特定染色体标准差(S.D._ChrN%)可根据来自例如至少10个,例如至少20个,例如至少30个,例如至少50个,例如至少100个标准细胞样本的该染色体ChrN%平均值确定。 
在本发明的一个实施方案中,标准细胞样本的个数为20个正常男性和10个正常女性,按顺序编号为1、2...30,其中1~20为正常男性受试样本,21~30为正常女性受试样本,则标准细胞样本的ChrN%平均值(mean_ChrN%)计算如下: 
Figure BSA00000228720400041
(其中N为1号~22号常染色体,M为1~30号正常样本)) 
Figure BSA00000228720400042
(M为1~20号男性正常样本) 
Figure BSA00000228720400043
(M为1~20号男性正常样本) 
Figure BSA00000228720400044
(M为21~30号女性正常样本) 
Figure BSA00000228720400045
(M为21~30号女性正常样本) 
(注:由于测序的波动性以及参考序列的Y染色体存在大量的缺口, 导致即使是正常的女性样本也会有少数DNA序列比对到Y染色体上,但与男性相比,女性的ChY%比男性要小很多,如实施例中女性的ChY%为0.004左右,而男性的ChY%为0.114左右) 
对根据上述方法获得的每一个标准细胞样本ChrN%平均值(mean_ChrN%)计算ChrN%标准差(S.D._ChrN%),计算公式为: 
(其中N为1号~22号常染色体) 
Figure BSA00000228720400052
(M为1~20号男性正常样本) 
Figure BSA00000228720400053
(M为1~20号男性正常样本) 
Figure BSA00000228720400054
(M为21~30号女性正常样本) 
Figure BSA00000228720400055
(M为21~30号女性正常样本) 
由于男性比女性的染色体中缺少一条染色体X,取而代之的是染色体Y,而X染色体全长约155M,而染色体Y仅为59M,因此在对性染色体进行检测时必须建立一套不同性别的ChrX%或ChrY%正态分布曲线,才能根据不同的性别正态分布曲线对性染色体的情况得出最准确的分析结果。 
在本发明的一个实施例中,取30例标准细胞样本进行染色体分析,根据模拟的染色体数量与标准细胞存在差异情况下达到正态分布的显著性要求(例如0.1%)建立正态分布曲线图,从而将z score_ChrN的绝对值确定在小于3的情形,定义为染色体数量相对于标准细胞而言相同。根据上述结果对待检样本进行染色体分析具体如下: 
当z score的绝对值=3时,则该样本99.9%的可能性被拒绝在一个 正态分布的群体中,即为离群样本,即:所检测的细胞其染色体数量相对于标准细胞而言存在差异; 
若z score的绝对值<3,则该样本99.9%的可能性是正常样本,即:所检测的细胞其染色体数量与标准细胞相同; 
若z score的绝对值>3,则该样本99.9%的可能性是异常样本,即:所检测的细胞其染色体数量相对于标准细胞而言存在差异; 
进一步的,在本发明中,当z score的绝对值>3时,对于所述细胞中染色体数量异常的具体情况,可以使用参考值Z(cutoff值)进行判断,参考值Z计算公式如下: 
Z=(mean_ChrN%×0.5×X%)/S.D._ChrN% 
其中,当N为常染色体时,mean_ChrN%和S.D._ChrN%分别为所有标准细胞受试样本的平均值,当N为性染色体时,mean_ChrN%和S.D._ChrN%分别为各自性别的标准细胞受试样本的平均值; 
其中X可以为负整数100(即-100)至正整数100之间的任意整数,例如-100、-90、-80、-70、-60、-50、-40、-30、-20、-10、0、10、20、30、40、50、60、70、80、90、100。 
在本发明的一个实施方案中,当X为100时,表示待测细胞样本比标准细胞的染色体多一条染色体;在本发明的一个实施方案中,当X为-100时,表示待测细胞样本比标准细胞的染色体少一条染色体;当X为50时,表示待测细胞样本比标准细胞的染色体多半条染色体;在本发明的一个实施方案中,当X为-50时,表示待测细胞样本比标准细胞的染色体少半条染色体; 
在本发明中,当计算性染色体的参考值Z时,其中mean_ChrN%和S.D._ChrN%的值为男性和女性样本分开计算的平均值,即: 
对于男性样本可以达到的参考值Z为:(mean_ChrN%(男性)×0.5×X%)/S.D._ChrN%(男性); 
对于女性样本可以达到的参考值Z为:(mean_ChrN%(女性)×0.5×X%)/S.D._ChrN%(女性)。 
当所述的z score_ChrN的绝对值大于等于3且达到参考值Z的绝 对值时,所述细胞的特定染色体数量与标准细胞相比存在X%的差异,例如X为50时,所述细胞的特定染色体比标准细胞多半条,X为100时,所述细胞的特定染色体比标准细胞多一条。 
在本发明的一个实施方案中,针对羊水细胞进行染色体分析的具体方法包括以下步骤: 
1).DNA提取及测序: 
根据Tiangen Micro Kit操作手册提取羊水细胞DNA后,按照修改过的Illumina/Solexa标准建库流程进行建库,在随机打断的DNA分子两端加上测序用的接头。在这个过程中,每个样本被加上不同的标签序列(index),这样多个样本可以在一次测序得到的数据中可以区分开。 
2).比对和统计: 
利用第二代测序方法Illumina/Solexa测序(用其它测序方法如ABI/SOLiD能达到相同或相近的效果)后每个样本得到一定大小片段的DNA序列,将其与人类基因组参考序列进行严格比对,得到所测序列定位于基因组相应位置的信息。 
这样严格的比对条件是因为容错比对或能比对到多个位置的DNA序列不能确定该DNA来源于哪条染色体,不利于后续的数据分析。 
以每条染色体为单位,计算出定位到每条染色体上的唯一序列的数量总和,由此得出每条染色体的ChrN%,即所测DNA序列中定位于该染色体N上的唯一位置的序列的数量总和(S1)占所测DNA序列中定位于所有染色体上的唯一位置的序列的数量总和(S2)的比例:ChrN%=S1/S2。 
这是一个将不同样本测序总量不同时的一个均一化的方法,由于羊水细胞包含46条染色体全部的信息,因此理论上定位到某条染色体上的序列的总数与该染色体长度成正比。 
举例来说,1号染色体是人类基因组中最大的染色体(约247M),而21号染色体是人类基因组中最小的染色体(约47M),因此对于正常的二倍型羊水细胞得到的测序结果在一定的测序总量下近似是一个定值。尽管在一定的测序条件和实验条件下,ChrN%并不完全与染色体 大小成正比,但基本上是定值。 
3).数据分析: 
通过箱式图分析,可能直接判定所述待测细胞样本是否有可能是染色体数量与标准细胞存在差异的样本,将出现异常值的样本直接作为可疑样本,其余的样本则作为标准细胞样本。具体过程如下: 
在本发明中,为了方便进行数据分析,对于根据上述方法计算获得的数据ChrN%,采用绘制箱式图(数理统计中用于区分异常值)的方法,确定可疑样本,具体绘制过程如下: 
1)计算上四分位数(75%),中位数(50%),下四分位数(25%) 
2)计算上四分位数和下四分位数之间的差值,即四分位数差(IQR,interquartile range) 
3)绘制箱式图的上下范围,上限为上四分位数,下限为下四分位数。在箱子内部中位数的位置绘制横线。 
4)大于上四分位数1.5倍四分位数差的值,或者小于下四分位数1.5倍四分位数差的值,划为异常值(outliers)。 
5)异常值之外,最靠近上边缘和下边缘的两个值处,画横线,作为箱式图的触须。 
6)极端异常值,即超出四分位数差3倍距离的异常值,用星点表示;较为温和的异常值,即处于1.5倍-3倍四分位数差之间的异常值,用空心点表示。 
箱式图箱子中间的粗线代表中值,上下边界代表上下四分位数。而异常点是指与上下边界相比,偏离上下分位数之差1.5倍距的点。例如,当检测的样本为标准细胞样本,其染色体对应的ChrN%为一定值,例如为1时,则当某个检测样本其相应染色体对应的ChrN%为1.5,则可以直接认定差异非常显著,为可疑样本,即有可能是染色体数量与标准细胞存在差异的样本。 
如需要,可以根据标准细胞样本对应的特定染色体ChrN%分别确定ChrN%平均值(mean_ChrN%)和标准差(S.D._ChrN%),计算可疑样本对应的染色体的z score_ChrN,公式如下: 
Z score_ChrN=(可疑样本特定染色体的ChrN%-mean_ChrN%)/S.D._ChrN%。 
如果z score_ChrN绝对值大于等于3,所述细胞样本对应的特定染色体的数量与标准细胞相比存在差异。 
进一步的,在本发明中,对于所述细胞中染色体数量异常的具体情况,可以使用参考值Z(cutoff值)进行判断,具体计算公式如下: 
Z=(mean_ChrN%×0.5×X%)/S.D._ChrN% 
其中,当N为常染色体时,mean_ChrN%和S.D._ChrN%分别为所有标准细胞受试样本的平均值,当N为性染色体时,mean_ChrN%和S.D._ChrN%分别为各自性别的标准细胞受试样本的平均值; 
其中X选定为50或100。相应的,当X为100时,表示待测细胞样本比标准细胞的染色体多一条染色体;当X为50时,表示待测细胞样本比标准细胞的染色体多半条染色体; 
在本发明中,当计算性染色体的参考值Z时,其中mean_ChrN%和S.D._ChrN%的值为男性和女性样本分开计算的平均值,即: 
对于男性样本可以达到的参考值Z为:(mean_ChrN%(男性)×0.5×X%)/S.D._ChrN%(男性); 
对于女性样本可以达到的参考值Z为:(mean_ChrN%(女性)×0.5×X%)/S.D._ChrN%(女性)。 
当z score_ChrN的绝对值大于等于3且达到参考值Z的绝对值时,所述细胞的特定染色体数量与标准细胞相比存在X%的差异。 
发明的有益效果 
本发明可以用于细胞,例如羊水细胞的染色体分析。在本发明中,由于待测的羊水细胞可不进行传代培养,直接提取DNA经测序后分析,大大降低了因羊水细胞培养造成的例如羊水细胞难贴壁、羊水细胞数量不够或培养失败等困难。 
利用羊水细胞含有胎儿完整基因组信息的特点,本发明所述方法能对细胞中所有染色体的非整倍性进行分析,而不是只针对21、18、13、 号染色体和性染色体X、Y进行检测。 
此外,尽管本发明的判定方法同样依赖于标准细胞样本建立的近似正态分布,但相对于血浆样本来说,对标准细胞样本的依赖程度大大降低,且当数据量足够时,可直接从数据异常中判断出异常样本。 
采用本发明的方法,可以对大量待测细胞样本进行批量分析,能一次对成百上千个待测细胞样本进行检测,大大节省人力物力。 
附图说明
图1根据53个待测细胞样本的ChrN%绘制的箱式图,其中: 
横坐标为染色体号,纵坐标为ChrN%值; 
图1A为1-6号染色体,图1B为7-12号染色体,图1C为13-17号染色体,图1D为18-22号染色体。 
具体实施方式
下面将结合实施例对本发明的实施方案进行详细描述,但是本领域技术人员将会理解,下列实施例仅用于说明本发明,而不应视为限定本发明的范围。实施例中未注明具体条件者,按照常规条件或制造商建议的条件进行。所用试剂或仪器未注明生产厂商者,均为可以通过市购获得的常规产品。所使用的测序用的接头和标签序列index来源于Illumina公司的Multiplexing Sample Preparation Oligonutide Kit。 
以下括号内为各个试剂或试剂盒的厂家货号。 
实施例1对未培养的羊水细胞进行染色体分析 
1.DNA提取及测序: 
按照Tiangen Micro Kit(DP316)微量基因组操作流程提取羊水细胞DNA,用Qubit(Invitrogen,the Quant-iTTM dsDNA HS Assay Kit)定量,所提取的DNA总量为100~500ng不等。 
提取的羊水细胞DNA是完整的基因组DNA或有部分降解呈smear 状DNA,按照修改后的Illumina/Solexa标准建库流程进行建库,在随机打断的DNA分子两端加上测序用的接头,并被加上不同的标签序列index,然后与flowcell表面互补接头杂交,在一定条件下使核酸分子成簇生长,然后在Illumina Genome AnalyzerⅡ上通过36轮测序循环,得到长度为35-bp的DNA片段序列。 
具体的,将获自羊水细胞的DNA约100~500ng使用Diagenode BioruptorTM随机打断至300bp片段,以起始打断DNA量为100~500ng,进行Illumina/Solexa标准建库,具体流程参照现有技术(参见 http://www.illumina.com/提供的Illumina/Solexa标准建库说明书)。经2100Bioanalyzer(Agilent)确定DNA文库大小及插入片段为300bp,QPCR精确定量后可上机测序。 
本实施例中,对于获自羊水细胞的53例DNA样本分批按照Illumina/Solexa官方公布的Cluster Station和GA Ⅱx(SE sequencing)说明书进行上机测序操作。 
2.比对和统计: 
参照现有技术(参见http://www.illumina.com/提供的Pipeline方法说明书),将步骤1中测得的序列信息经过一个Pipeline过程,去掉测序质量低的序列,最终可以获得针对NCBI版本36的人类基因组参考序列的ELAND比对结果,对所有定位到染色体上的唯一序列数量进行数据统计。 
计算53个样本在22条染色体和X/Y染色体上的ChrN%,做箱式图(见图1),其中计算针对特定样本M中特定染色体N的ChrN%公式如下: 
测试样本M中特定染色体的百分比ChrN%=该样本M中比对到参考序列相应染色体上的唯一序列的总数(S1)/该样本M中比对到参考序列所有染色体上的唯一序列的总数(S2)。 
3.数据分析: 
根据步骤2中获得的箱式图,首先确定是否存在异常点,即是指与上下边界相比,偏离上下分位数之差1.5倍距的点较大偏离的可疑样本,有可能其染色体数量与标准细胞样本存在差异。 
具体的,观察箱式图分布,检出可疑样本8例(样本编号为P1-P8)。将剔除可疑样本后剩余的45例标准细胞样本中随机挑出20个正常男性和10个正常女性的数据作为标准样本建立正态分布,计算每条染色体的ChrN%均值mean(mean_ChrN%)和标准差S.D.值(S.D._ChrN%)(见表1)。 
表1、标准细胞中每条染色体的ChrN%,以及均值(mean)和标准差(S.D.)。 
Figure BSA00000228720400121
Figure BSA00000228720400141
进一步的,为考察所得可疑样本中是否存在多半条或多一条染色体的情形,选定X为50或100,计算出对应的染色体参考值Z(cutoff值)(参见表2): 
Z=(mean_ChrN%×0.5×X%)/S.D._ChrN%,其中N为常染色体1~22号,X为50或100。 
表2、判定待测细胞中染色体数量为三体的参考值Z(cutoff值)。 
Figure BSA00000228720400151
根据计算公式,计算所测可疑样本中每条染色体的z score_ChrN。计算公式如下: 
z score_ChrN=(受检样本某染色体的ChrN%-mean_ChrN%)/S.D._ChrN%。 
表3可疑样本中对应于每条染色体的z score_ChrN值 
Figure BSA00000228720400161
由上述分析可知,53例受检羊水细胞样本中可疑样本共8个,其中对于每个可疑样本的染色体,共捡出z score_ChrN值绝对值>3的染色体数量异常合计8个(参见表3), 
具体为 
1)P1的Chr21,P2的Chr21,P3的Chr21和P4的Chr21; 
2)P5的Chr18,P6的Chr18和P7的Chr18; 
3)P8的Chr13。 
通过查对表2中根据多一条染色体(X=100)计算获得的Z值,可以确定其中P1-P4样本中染色体21的数量和P5-P7中染色体18的数量比标准细胞中的相应染色体数量多一条,P8中染色体13的数量比标准细胞中的相应染色体数量多半条,与传统染色体核型分析结果完全一致。 
实施例2 
另有6例羊水细胞样本(Q1-Q6)在经过同样的处理方法和测序过程后获得上机数据,以实施例1中30例标准细胞样本计算出来的mean_ChrN%和S.D._ChrN%予以计算z score_ChrN。6例样本中检出阳性样本3例(参见表4和表5)。 
表4六例待测样本(Q1-Q6)的ChrN% 
Figure BSA00000228720400171
Figure BSA00000228720400181
表5以实施例1中30例标准细胞样本的mean和S.D.计算得出的6例样本的z score_ChrN 
Figure BSA00000228720400182
Figure BSA00000228720400191
从结果可以看出,Q5的染色体21的数量比标准细胞相应染色体数量多一条,Q3和Q4的性染色体X的数量比标准细胞相应染色体数量少一条,与临床核型分析结果完全一致。 
尽管本发明的具体实施方式已经得到详细的描述,本领域技术人员将会理解。根据已经公开的所有教导,可以对那些细节进行各种修改和替换,这些改变均在本发明的保护范围之内。本发明的全部范围由所附权利要求及其任何等同物给出。 

Claims (7)

1.一种通过测序分析细胞染色体的方法,其包括以下步骤:
a.将获自细胞的DNA分子随机打断,得到一定大小的DNA片段并测序;
b.将步骤a中测定的DNA序列与人类基因组参考序列进行严格比对,得到所测DNA序列定位于具体染色体的信息;
c.对应于特定染色体N,确定所测DNA序列中定位于该染色体上的唯一位置的序列的数量总和,进而得出对应于该染色体N的ChrN%,即所测DNA序列中定位于该染色体N上的唯一位置的序列的数量总和占所测DNA序列中定位于所有染色体上的唯一位置的序列的数量总和的比例:
ChrN%=定位于特定染色体N上的唯一位置的序列的数量总和/定位于所有染色体上的唯一位置的序列的数量总和;
d.将该染色体的ChrN%与来自标准细胞中对应染色体的ChrN%进行比较,确定所述细胞的该条染色体数量与标准细胞是否存在差异。
2.权利要求1所述的方法,其中步骤b中所述严格比对是指不容错的、定位于人类基因组参考序列上的唯一位置的比对;其中所述的人类基因组序列参考序列是人类基因组中经过屏蔽掉重复序列后所得到的序列。
3.权利要求1所述的方法,其中步骤d中所述细胞样本中特定染色体数量与标准细胞是否存在差异是通过绘制箱式图方法确定的,其中将位于超出四分位数差1.5倍-3倍或者3倍以上的四分位数差之间的异常值对应的样本,认定为其染色体数量与标准细胞存在差异。
4.权利要求1所述的方法,其中步骤d中所述细胞的特定染色体数量与标准细胞是否存在差异是用反映受检细胞样本ChrN%与标准细胞样本ChrN%的偏离程度z score_ChrN来确定的,如果zscore_ChrN绝对值大于等于3,所述细胞样本对应的特定染色体的数量与标准细胞相比存在差异。
5.权利要求4所述的方法,其中
所述z score_ChrN=(受检样本特定染色体的ChrN%-特定染色体ChrN%平均值(mean_ChrN%))/ChrN%平均值标准差(S.D._ChrN%);
其中所述特定染色体ChrN%平均值可根据来自至少10个,例如至少20个标准细胞样本的该染色体ChrN%确定;
其中所述特定染色体ChrN%平均值标准差可根据来自至少10个,例如至少20个标准细胞样本的该染色体ChrN%平均值确定。
6.权利要求4所述的方法,其中步骤d中所述细胞的特定染色体数量与标准细胞是否存在差异是用z score_ChrN与参考值Z进行比较加以确定,其中参考值Z根据如下方法确定:
Z=(mean_ChrN%)×0.5×X%/(S.D._ChrN%)
其中X可以为负整数100(即-100)至正整数100之间的任意整数,例如-100、-90、-80、-70、-60、-50、-40、-30、-20、-10、0、10、20、30、40、50、60、70、80、90、100;
当z score_ChrN的绝对值大于等于3且达到参考值Z的绝对值时,所述细胞的特定染色体数量与标准细胞相比存在X%的差异。
7.权利要求1所述的方法,其中所述细胞为人羊水细胞,例如未培养的羊水细胞或培养过的羊水细胞,优选为未培养的羊水细胞。
CN201010252849XA 2010-08-13 2010-08-13 一种细胞染色体分析方法 Pending CN101928775A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010252849XA CN101928775A (zh) 2010-08-13 2010-08-13 一种细胞染色体分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010252849XA CN101928775A (zh) 2010-08-13 2010-08-13 一种细胞染色体分析方法

Publications (1)

Publication Number Publication Date
CN101928775A true CN101928775A (zh) 2010-12-29

Family

ID=43368204

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010252849XA Pending CN101928775A (zh) 2010-08-13 2010-08-13 一种细胞染色体分析方法

Country Status (1)

Country Link
CN (1) CN101928775A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014075228A1 (zh) * 2012-11-13 2014-05-22 深圳华大基因医学有限公司 确定生物样本中染色体数目异常的方法、系统和计算机可读介质
CN103824001A (zh) * 2014-02-27 2014-05-28 北京诺禾致源生物信息科技有限公司 染色体的检测方法和装置
CN109243536A (zh) * 2012-12-19 2019-01-18 阿瑞奥萨诊断公司 卵子供体妊娠中胎儿非整倍性的无创检测

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014075228A1 (zh) * 2012-11-13 2014-05-22 深圳华大基因医学有限公司 确定生物样本中染色体数目异常的方法、系统和计算机可读介质
CN104093858A (zh) * 2012-11-13 2014-10-08 深圳华大基因医学有限公司 确定生物样本中染色体数目异常的方法、系统和计算机可读介质
CN104093858B (zh) * 2012-11-13 2016-08-24 深圳华大基因医学有限公司 确定生物样本中染色体数目异常的方法、系统和计算机可读介质
CN109243536A (zh) * 2012-12-19 2019-01-18 阿瑞奥萨诊断公司 卵子供体妊娠中胎儿非整倍性的无创检测
CN109243536B (zh) * 2012-12-19 2021-12-10 豪夫迈·罗氏有限公司 卵子供体妊娠中胎儿非整倍性的无创检测
CN103824001A (zh) * 2014-02-27 2014-05-28 北京诺禾致源生物信息科技有限公司 染色体的检测方法和装置

Similar Documents

Publication Publication Date Title
Ulirsch et al. The genetic landscape of Diamond-Blackfan anemia
Choin et al. Genomic insights into population history and biological adaptation in Oceania
CN103403183B (zh) 胎儿遗传异常的无创性检测
CN104232777B (zh) 同时确定胎儿核酸含量和染色体非整倍性的方法及装置
CN104830986A (zh) 一种检测胎儿基因信息的方法、装置和系统
Keller et al. Trans-ethnic meta-analysis of white blood cell phenotypes
CN105844116B (zh) 测序数据的处理方法和处理装置
TR201816062T4 (tr) Polimorfizm sayımlarını kullanarak genom fraksiyonlarının çözülmesi.
CN102753703B (zh) 胎儿染色体非整倍性的检测方法
CN105825076B (zh) 消除常染色体内和染色体间gc偏好的方法及检测系统
KR101801871B1 (ko) 산모의 혈청 dna를 이용한 태아의 단일유전자 유전변이의 예측방법
CN104169929A (zh) 用于确定胎儿是否存在性染色体数目异常的方法、系统和计算机可读介质
CN101928775A (zh) 一种细胞染色体分析方法
CN104520437B (zh) 一种染色体非整倍性检测方法及装置
Sun et al. Comparative genomic analysis reveals evolutionary characteristics and patterns of microRNA clusters in vertebrates
CN104093858B (zh) 确定生物样本中染色体数目异常的方法、系统和计算机可读介质
CN105821117B (zh) 多倍染色体检测用的标准品产品及其制备方法
CN105765076B (zh) 一种染色体非整倍性检测方法及装置
CN116287204A (zh) 检测特征基因的突变情况在制备静脉血栓栓塞症风险检测产品中的应用
CN103069004B (zh) 一种细胞染色体分析方法
McCaffrey et al. RNAseq profiling of blood from patients with coronary artery disease: signature of a T cell imbalance
CN104178572B (zh) 测序数据的处理方法及装置
CN109321641B (zh) 一种基于dna片段富集及测序技术的产前无创胎儿染色体检测系统
Teo et al. A novel method for haplotype clustering and visualization
CN104978497B (zh) 一种大样本高通量生物数据关联分析方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
REG Reference to a national code

Ref country code: HK

Ref legal event code: DE

Ref document number: 1151560

Country of ref document: HK

C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20101229

REG Reference to a national code

Ref country code: HK

Ref legal event code: WD

Ref document number: 1151560

Country of ref document: HK