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

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

Info

Publication number
CN103069004A
CN103069004A CN2010800685719A CN201080068571A CN103069004A CN 103069004 A CN103069004 A CN 103069004A CN 2010800685719 A CN2010800685719 A CN 2010800685719A CN 201080068571 A CN201080068571 A CN 201080068571A CN 103069004 A CN103069004 A CN 103069004A
Authority
CN
China
Prior art keywords
chrn
chromosome
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.)
Granted
Application number
CN2010800685719A
Other languages
English (en)
Other versions
CN103069004B (zh
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.)
Shenzhen Huada Medical Laboratory
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
Publication of CN103069004A publication Critical patent/CN103069004A/zh
Application granted granted Critical
Publication of CN103069004B publication Critical patent/CN103069004B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869—Methods for sequencing
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6813—Hybridisation assays
    • C12Q1/6827—Hybridisation assays for detection of mutation or polymorphism
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/156—Polymorphic or mutational markers

Abstract

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

Description

一种细胞染色体分析方法 技术领域
本发明涉及一种细胞染色体的分析方法， 特别是涉及通过测序方法 分析羊水细胞染色体的方法。 背景技术
胎儿染色体非整倍性
胎儿染色体非整倍性是指胎儿基因组出现染色体数量不是二倍型 的情况。 正常情况下， 人的基因组是 44条常染色体， 2条性染色体， 男 性核型为 （46， XY) ， 女性核型为 （46, XX) 。 胎儿染色体非整倍性， 可以是较正常的二倍体胎儿多了一条染色体， 即胎儿基因组为 47条染色 体， 以胎儿 21号三体为例， 患儿较正常的二倍体胎儿多了一条 21号染色 体， 核型为 47, XX (或 XY) , +21; 可以较正常的二倍体胎儿少了一条 染色体， 即胎儿基因組为 45条染色体， 以特纳综合征患者为例， 患儿较 正常的二倍体胎儿缺了一条 X染色体， 核型为 45， X0; 可以较正常的二 倍体胎儿少了一部分染色体，如易位型 21三体，核型为 45, XX, der(14; 21) (qlO; qlO) , 又如猫叫综合征，核型为 46, XX (XY) ,del (5) (pl3)0 据不完全统计， 全世界染色体非整倍性的胎儿出生率为 1/160， 其 中 21三体（T21, 唐氏综合征）胎儿的出生率是 1/δΟΟ, I8三体（Τΐδ, Edwards综合征）胎儿的出生率是 1/6000， 13三体（T13， Patau综合 征）胎儿的出生率是 1/10000， 其他类型的染色体非整倍性患儿由于难 以完成胚胎发育的某些阶段而停滞， 临床表现为早孕期的不明流产 ( Deborah A. Driscoll, M. D. , and Susan Gross, M. D. Prenatal Screening for Aneuploidy [J] . N Engl J Med 2009;360:2556-62 ) 。
羊水细胞培养的现状 羊水细胞是羊水中漂浮着的来源于胎儿皮肤、 消化道、 呼吸道等脱 落的上皮细胞。 利用羊水细胞培养， 富集胎儿有核细胞的数量， 制备染 色体标本， 分析胎儿染色体核型， 是传统临床检测染色体核型异常胎儿 的金标准， 这种技术可靠、 准确， 能观察到染色体数量和结构的异常， 但其缺点是耗时长， 需要 10天至 3 周才能得到结果， 且培养失败率约 为 1. 00°/。（THEIN AT, ABDEL- FATTAH SA， KYLE PM , et al. an assessment of the use of interphase FISH wi th chromosome specific probes as an al ternat ive to cytogenet ics in prenatal diagnos i s [J] . PrenatDiagn, 2000, 20 (4) : 275-280· ) 。
羊水细胞培养失败率高是由于羊水细胞大部分是衰老和固缩的脱 落细胞， 使得羊水细胞培养较其他组织培养更困难（马长俊， 陈园茶， 霍沛丹。 羊水培养与染色体制片方法 [ J 3. 生殖与避孕， 1985，5 (1) : 53 - 4. ) 。 而羊水细胞培养的成功在染色体非整倍性的检测过程 中具有关键作用。 由于羊水细胞培养的要求较高，部分孕妇的羊水中有 活力的细胞相对较少，收获到的细胞分裂相较少，染色体形态不佳，难以 进行足够的细胞计数分析，给检测带来囷难。 加上羊膜腔穿刺有一定的 风险， 约 2 % ~ 3 %孕妇可出现子宫收缩、 腹部肿胀、 压痛、 阴道出血、 感染、 漏水或胎儿损伤等并发症。 若羊水细胞培养失败或收获后计数分 析不够而要求孕妇行笫 2 次穿刺，一般难以接受； 且一般穿刺时间为孕 16-20周， 若培养失败一次， 很多孕妇孕周已超过检测期， 不得不进行 风险更高的脐带血穿刺。 并且羊水培养后核型分析需要大量的人力、 物 力， 很多医院根本无力承担， 对羊水穿刺技术在临床上的推广使用也造 成很大的囷难。
随着测序技术的不断发展，测序技术在染色体数目的检测分析中得 到了越来越广泛的应用。 Dennis Lo 等人以孕妇外周血为实验材料， 通 过大规模测序， 基于数理统计方法进行染色体数目异常的检测 ( Y. M. Denni s Lo, et al. Quant itat ive Analys is of Fetal DNA in Maternal Plasma and Serum: Impl icat ions for Noninvas ive Prenatal Diagnos i s. Am. J. Hum. Genet. 62: 768 775, 1998 ) 。 但是这种方法 除开技术本身存在一定的缺陷外，应用上也不能完全取代羊水穿刺的作 用。 由于血浆中游离 DNA是片段化的 DM, 本身并不能拼凑出完整的基 因組， 因此对 21、 18、 13号染色体以外的其他染色体非整倍性或者染 色体易位、 镶嵌体等无法进行检测或检测精度不高。 发明内容
为了克服外周取血测序的漏检现象， 同时解决羊水细胞培养失败率 高的现状，本发明结合羊水穿刺进行核型分析和血浆游离 DNA测序方法 的优点， 采用基于羊水细胞的大规模测序进行染色体非整倍性的检测方 法， 即抽取羊水细胞， 提取 DM, 进行高通量测序， 对获得的数据进行 分析， 得出检测结果， 由此完成了本发明。
本发明的一个方面， 涉及一种用测序方法分析受试者细胞染 色体信息的方法， 其包括以下步骤：
a.将获自细胞的基因组 DNA分子随机打断， 得到一定大小的 DNA片段并测序；
b.将步驟 a中测定的 DNA序列与人类基因组参考序列进行严 格比对， 得到所测 DNA序列定位于具体染色体的信息；
c.对应于特定染色体 N， 确定所测 DNA序列中定位于该染色 体上的唯一位置的序列的数量总和， 进而得出对应于该染色体 N 的 ChrN%， 即所测 DNA序列中定位于该染色体 N上的唯一位置的 序列的数量总和（S 1)占所测 DNA序列中定位于所有染色体上的唯 一位置的序列的数量总和（S2)的比例： Chr N%-S 1 /S2;
d .将该染色体的 C h r N。/。与来自标准细胞中对应染色体的 ChrN/。进行比较，确定所述细胞的该条染色体与标准细胞是否存在 差异。
本发明中， 所述的细胞可以是例如羊水细胞， 其中所述的羊水 细胞可以为未培养的羊水细胞或培养过的羊水细胞。在本发明的一个实 施方案中， 为避免培养羊水细胞， 羊水细胞为未培养的羊水细胞。
本发明中，细胞基因组 DM的获取可以采用盐析法、柱层析法、 SDS 法等常规 DNA提取方法， 优选采用柱层析法。 所谓柱层析法， 是指羊水 细胞或组织经过细胞裂解液和蛋白酶 K的作用后露出棵露的 DNA分子， 其经过能与带负电的 DNA分子结合的硅胶膜柱时， 体系中的基因組 DNA 被可逆吸附， 经漂洗液清洗除去蛋白质、 脂质等杂质后， 用純化液洗脱 获得羊水细胞 DNA (具体原理和方法参见 Qiagen公司货号为 56304的产 品说明书和 Tiangen公司货号为 DP316的产品说明书）。
本发明中， 所述 DNA分子的随机打断处理可以采用酶切、 雾 化、超声、或者 HydroShear法。优选地， 采用 HydroShear法（当 舍有 DNA的溶液通过较小面积的通道时， 流体加速， 产生的力使 DNA突然断裂， 流速和通道大小决定 DNA片段的大小， 具体原理 和方法请参见 Life Sciences Wiki公司的 HydroShear说明书） , 将 DM分子打断为比较集中的一定大小的片段，通常经过随机打断的主 带分布在 200 ~ 300bp范围。
本发明中， 所采用的测序方法可以为第二代测序方法 Il lumina/Solexa或 ABI/S0LiD。 在本发明的一个实施方案中， 所述的 测序方法为 I Uumina/Solexa, 测序得到的结果为 35bp大小的片段。
当待测的 DM分子来自多个受试样本时，每个样本可以被加上不同 的标签序列 index, 以用于在测序过程中进行（Micah Hamady, Jeffrey J Walker, J Kirk Harris et al. Error-correct ing barcoded primers forpyrosequencing hundreds of samples in mul t iplex. Nature Methods, 2008, March, Vol. 5 No. 3)。
本发明中，所述的人类基因組参考序列是人类基因组序列经过屏蔽 掉重复序列后所得到的参考序列， 例如 NCBI数据库中最新版本的人类 基因组参考序列。 在本发明的一个实施方案中， 所述人类基因组序列是 NCBI数据库中版本 36 ( NCBI Bui ld 36 ) 的人类基因组参考序列。
在本发明中， 其中所述与人类基因组参考序列进行严格比对， 是指 采用的比对方法为不容错的、 定位于人类基因组参考序列上的唯一 位置的比对。 在本发明的一个实施方案中， 采用的是 Eland 比对软件 ( I l lumina提供的软件包） ， 采用的方法为完全不容错比对。
本发明中，当所述的 DM序列是能够定位于人类基因組参考序列的 唯一位置上的 DNA序列时， 将其定义为唯一序列， 以 Unique reads表 示。 在本发明中， 为了避免重复序列的干扰， 需要去除那些定位于人类 基因组参考序列中的串联重复及转座重复位置的 DNA序列，只统计那 些可以定位到基因组唯一位置的 DNA序列， 即唯一序列。 通常， 在 所有测序的 DNA序列中， 大约有四分之一到三分之一的 DNA序列可 以定位于基因組唯一位置， 即唯一序列。 而这个统计所得的唯一序列的 数量代表分布在基因组染色体上的 DNA序列的情况。
因此，唯一序列能够将来自羊水细胞的 DNA分子经打断并测序后的 各 DNA序列定位于特定染色体。 而 ChrN%则是将不同染色体上得到的唯 一序列进行标准化后得到的数值， 该数值只与特定染色体大小有关， 而 与测序数据量无关，从而可以用（¾1： 对 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为正常女性受试样本， 则标准细胞 样本的 Chr 平均值 （mean_ChrN%) 计算如下：
(其中 N为 1号~ 22号常染色体，
M为 1 ~ 30号正常样本） ）
1 20
Mean_ChrX% (男性） =丄 |； Chrx M为 1 ~ 20号男性正常样本）
1 20
Mean— ChrY%(男性） = 丄 (Μ为 1 ~ 20号男性正常样
20 r{
本）
1 30
Mean_ChrX°/o (女性） = 丄 c½r/ ( M为 21 ~ 30号女性正常样 本） "ί 30
Mean-ChrY% (女性） = 丄 |； ChrY一 M ( M为 21 ~ 30号女性正常样 本）
(注： 由于测序的波动性以及参考序列的 Y染色体存在大量 的缺口， 导致即使是正常的女性样本也会有少数 DNA序列比对到 Y染色体上， 但与男性相比， 女性的 ChY%比男性要小很多， 如实 施例中女性的 0^%为 0.004左右， 而男性的（¾¥%为 0.114左右） 对根据上述方法获得的每一个标准细胞样本 ChrN°/。平均值 ( mean_ChrN% )计算 Chr 标准差 ( S. D. _C rN%) ， 计算公式为：
1 30
S. D. _ChrN%=~ ^ChrN—M - mean C r )2 (其中 N 为 1 号 ~ 22 号常染色体）
S.D.— ChrX°/。（男性） (M 为 1 ~ 20
号男性正常样本）
1 20
S. D. _ChrY% (男性） - 丄 X (C Y - mean CArY)2 ( M为 1 ~ 20号 男性正常样本）
1 30
S. D. _ChrX% (女性） - 丄；£ (C½rX— - mean C rl)2 ( M为 21 ~ 30
10 2l
号女性正常样本）
30
S. D. -ChrY°/o (女性） = 丄 " (C½rY - wean ChrY)2 ( M为 21 ~ 30
10 i 一 一
号女性正常样本） 由于男性比女性的染色体中缺少一条染色体 X， 取而代之的 是染色体 Y, 而 X染色体全长约 155M, 而染色体 Y仅为 59M, 因 此在对性染色体进行检测时必须建立一套不同性别的 C h r X %或 ChrY %正态分布曲线， 才能根据不同的性别正态分布曲线对性染 色体的情况得出最准确的分析结果。
在本发明的一个实施例中， 取 30例标准细胞样本进行染色体 分析，根据模拟的染色体数量与标准细胞存在差异情况下达到正态 分布的显著性要求（例如 0.1% )建立正态分布曲线图， 从而将 ζ score-ChrN的绝对值确定在小于 3的情形， 定义为染色体数量相 对于标准细胞而言相同。根据上述结果对待检样本进行染色体分析 具体如下：
当 z score的绝对值 =3时， 则该样本 99.9%的可能性被拒绝在 一个正态分布的群体中， 即为离群样本， 即： 所检测的细胞其染色 体数量相对于标准细胞而言存在差异， 即为染色体非整倍性；
若 z ore的绝对值 <3，则该样本 99.9%的可能性是正常样本， 即： 所检测的细胞其染色体数量与标准细胞相同；
若 z score的绝对值>3,则该样本 99.9%的可能性是异常样本， 即： 所检测的细胞其染色体数量相对于标准细胞而言存在差异， 即 为染色体非整倍性；
进一步的， 在本发明中， 当 z score的绝对值 >3时， 对于所 述细胞中染色体非整倍性的具体情况，可以使用参考值 Z( cutoff 值） 进行判断， 参考值 Z计算公式如下：
Z= (mean-ChrN %x0.5 X%) / S. D. -C rN %
其中， 当 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% (男性） xO. 5xX%) /S. D. -ChrN°/o (男性） ；
对于女性样本可以达到的参考值 Z为： ( mean_ChrN% (女性） xO. 5xX%) /S. D. -ChrN% (女性 ) 。
当所述的 z score-ChrN的绝对值大于等于 3且达到参考值 Z 的绝对值时， 所述细胞的特定染色体数量与标准细胞相比存在 X% 的差异，例如 X为 50时， 所述细胞的特定染色体比标准细胞多半 条， X为 100时， 所述细胞的特定染色体比标准细胞多一条。
在本发明的一个实施方案中， 针对羊水细胞进行染色体分析 的具体方法包括以下步骤：
1 ) . DNA提取及测序：
根据 Tiangen Micro Ki t操作手册提取羊水细胞 DNA后， 按照 修改过的 I l lumina/Solexa标准建库流程进行建库，在随机打断的 DNA分子两端加上测序用的接头。 在这个过程中， 每个样本被加上 不同的标签序列（index ), 这样多个样本可以在一次测序得到的数 据中可以区分开。
2) . 比对和统计: 利用第二代测序方法 I l lumina/Solexa测序（用其它测序方法 如 ABI/SOLiD能达到相同或相近的效果）后每个样本得到一定大小 片段的 DM序列， 将其与人类基因组参考序列进行严格比对， 得到 所测序列定位于基因组相应位置的信息。 这样严格的比对条件是因为容错比对或能比对到多个位置的 DNA序列不能确定该 DNA来源于哪条染色体， 不利于后续的数据分 析。 以每条染色体为单位，计算出定位到每条染色体上的唯一序列 的数量总和， 由此得出每奈染色体的 ChrN%， 即所测 DNA序列中定 位于该染色体 N上的唯一位置的序列的数量总和（S1)占所测 DNA 序列中定位于所有染色体上的唯一位置的序列的数量总和（S 2)的 比例： ChrN%=Sl/S2。 这是一个将不同样本测序总量不同时的一个均一化的方法， 由 于羊水细胞包含 46条染色体全部的信息， 因此理论上定位到某条 染色体上的序列的总数与该染色体长度成正比。
举例来说， 1号染色体是人类基因组中最大的染色体（约 247M ), 而 21号染色体是人类基因组中最小的染色体（约 47M ), 因此对于 正常的二倍型羊水细胞得到的测序结果在一定的测序总量下近似 是一个定值。 尽管在一定的测序条件和实验条件下， ChrN%并不完 全与染色体大小成正比， 但基本上是定值。
3 ). 数据分析：
通过箱式图分析，可能直接判定所述待测细胞样本是否有可能 是染色体数量与标准细胞存在差异的样本，将出现异常值的样本直 接作为可疑样本，其余的样本则作为标准细胞样本。具体过程如下： 在本发明中， 为了方便进行数据分析， 对于根据上述方法计算 获得的数据 ChrN %， 采用绘制箱式图 （数理统计中用于区分异常 值） 的方法， 确定可疑样本， 具体绘制过程如下：
1 ) 计算上四分位数 Π5% ) , 中位数（50% ) ， 下四分位数 ( 25% )
2 ) 计算上四分位数和下四分位数之间的差值， 即四分位数 差 ( IQR, interquart i le range )
3 ) 绘制箱式图的上下范围， 上限为上四分位数， 下限为下 四分位数。 在箱子内部中位数的位置绘制横线。
4) 大于上四分位数 1.5倍四分位数差的值，或者小于下四分 位数 1.5倍四分位数差的值， 划为异常值 （outliers) 。
5) 异常值之外， 最靠近上边缘和下边缘的两个值处， 画横 线， 作为箱式图的触须。
6) 极端异常值， 即超出四分位数差 3倍距离的异常值， 用星 点表示； 较为温和的异常值， 即处于 1.5倍 -3倍四分位数差之间的 异常值， 用空心点表示。
箱式图箱子中间的粗线代表中值， 上下边界代表上下四分位 数。 而异常点是指与上下边界相比， 偏离上下分位数之差 1.5倍 距的点。 例如， 当检测的样本为标准细胞样本， 其染色体对应的 ChrN%为一定值， 例如为 1 时， 则当某个检测样本其相应染色体对 应的 （¾^%为 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 x0.5xX ) / S. D. -ChrN %
其中， 当 N为常染色体时， mean— 01^%和 S.D. _ChrN%分别为 所有标准细胞受试样本的平均值，当 N为性染色体时， mean_ChrN% 和 S. D. -ChrN%分别为各自性别的标准细胞受试样本的平均值； 其中 X选定为 50或 100。 相应的， 当 X为 100时， 表示待测 细胞样本比标准细胞的染色体多一条染色体； 当 X为 50时，表示 待测细胞样本比标准细胞的染色体多半条染色体；
在本发明中，当计算性染色体的参考值 Z时，其中 mean_ChrN°/。 和 S. D. -01^%的值为男性和女性样本分开计算的平均值， 即： 对于男性样本可以达到的参考值 Z为： （mean-Chr （男性） x0. 5xX°/o) /S. D. _ChrN% (男性） ；
对于女性样本可以达到的参考值 Z为： （ mean— ChrN% (女性） x0. 5xX%) /S. D. _ChrN% (女性 ) 。
当 z score-ChrN的绝对值大于等于 3且达到参考值 Z的绝对 值时， 所述细胞的特定染色体数量与标准细胞相比存在 X%的差 异。 发明的有益效果
本发明可以用于细胞， 例如羊水细胞的染色体分析。 在本发 明中， 由于待测的羊水细胞可不进行传代培养，直接提取 DNA经测 序后分析， 大大降低了因羊水细胞培养造成的例如羊水细胞难贴 壁、 羊水细胞数量不够或培养失败等困难。
利用羊水细胞含有胎儿完整基因组信息的特点，本发明所述方 法能对细胞中所有染色体的非整倍性进行分析， 而不是只针对 21、 18、 13、 号染色体和性染色体 X、 Y进行检测。
此外 ,尽管本发明的判定方法同样依赖于标准细胞样本建立的 近似正态分布， 但相对于血浆样本来说， 对标准细胞样本的依赖程 度大大降低， 且当数据量足够时， 可直接从数据异常中判断出异常 样本。
采用本发明的方法， 可以对大量待测细胞样本进行批量分析， 能一次对成百上千个待测细胞样本进行检测， 大大节省人力物力。 附图说明
图 1 根据 53个待测细胞样本的 011^%绘制的箱式图， 其中： 横坐标为染色体号， 纵坐标为 ChrN%值；
图 1A为 1-6号染色体，图 1B为 7-12号染色体，图 1C为 13-17 号染色体， 图 1D为 18-22号染色体。 具体实施方式
下面将结合实施例对本发明的实施方案进行详细描述， 但是 本领域技术人员将会理解， 下列实施例仅用于说明本发明， 而不 应枧为限定本发明的范围。 实施例中未注明具体条件者， 按照常 规条件或制造商建议的条件进行。 所用试剂或仪器未注明生产厂 商者， 均为可以通过市购获得的常规产品。 所使用的测序用的接 头和标签序列 index 来源于 I l lumina 公司的 Mul t iplexing Sample Preparat ion 01 igonut ide Ki t；。
以下括号内为各个试剂或试剂盒的厂家货号。 实施例 1 对未培养的羊水细胞进行染色体分析
1. DNA提取及测序：
按照 Tiangen Micro Ki t (DP 316)微量基因组操作流程提取羊 7j 细包 DNA， 用 Qubi t ( Invi trogen , the Quant— iT dsDNA HS As say Ki t )定量， 所提取的 DNA总量为 100 ~ 500ng不等。
提取的羊水细胞 DNA 是完整的基因组 DNA 或有部分降解呈 smear状 DNA, 按照修改后的 11 lumina/Solexa标准建库流程进行 建库， 在随机打断的 DNA分子两端加上测序用的接头， 并被加上不 同的标签序列 index,然后与 f lowcel 1表面互补接头杂交，在一定 条件下使核酸分子成簇生长， 然后在 I l lumina Genome Analyzer II上通过 36轮测序循环， 得到长度为 35-bp的 DNA片段序列。
具体的，将获自羊水细胞的 DNA约 100 ~ 500ng使用 Diagenode Bioruptor 随机打断至 300bp片段， 以起始打断 DNA量为 100 ~ 500ng， 进行 Illumina/Solexa标准建库， 具体流程参照现有技术 (参见 http: //www, illumina. com/提供的 Illumina/Solexa标准 建库说明书）。 经 2100 Bioanalyzer (Agilent)确定 DNA文库大小 及插入片段为 300bp， QPCR精确定量后可上机测序。
本实施例中， 对于获自羊水细胞的 53 例 DNA样本分批按照 Illuraina/Sol ex 官方公布的 Cluster Stat ion 和 GA II x ( SE sequencing)说明书进行上机测序操作。
2. 比对和统计：
参照现有技术 （参见 http: //www, i 1 lumina. com/提供的 Pipeline 方法说明书）， 将步骤 1 中测得的序列信息经过一个 Pipeline 过程， 去掉测序质量低的序列， 最终可以获得针对 NCBI 版本 36的人类基因组参考序列的 ELAND比对结果，对所有定位到染 色体上的唯一序列数量进行数据统计。
计算 53个样本在 22条染色体和 X/Y染色体上的 ChrN%, 做箱 式图 （见图 2), 其中计算针对特定样本 M中特定染色体 N的 ChrN% 公式: ^下：
测试样本 M中特定染色体的百分比 ChrN%-该样本 M中比对到参 考序列相应染色体上的唯一序列的总数（S1) /该样本 M中比对到 参考序列所有染色体上的唯一序列的总数（S2)。
3. 数据分析：
根据步驟 2中获得的箱式图， 首先确定是否存在异常点， 即是 指与上下边界相比， 偏离上下分位数之差 1.5倍距的点较大偏离 的可疑样本， 有可能其染色体数量与标准细胞样本存在差异。
具体的， 观察箱式图分布， 检出可疑样本 8 例 （样本编号为 Pl-P8)。 将剔除可疑样本后剩余的 45例标准细胞样本中随机挑出 20 个正常男性和 10 个正常女性的数据作为标准样本建立正态分 布， 计算每条染色体的 ChrN°/。均值 mean ( mean.ChrN% )和标准差 S. D.值 （ S. D. _ChrN% ) (见表 1 )。
表 1、标准细胞中每条染色体的 ChrN%，以及均值（ mean；)和标准差（ S 序号 C rl½ Chr2% Chr3½ Chr5½ Chr6% Chr7¾ ChtU Chr9¾ Chrl0¾ ChrllB/o Chrl2%
1 7.850 9.295 7.542 7.381 6.876 6.660 5.499 5.468 3.942 4.747 - 4.603 4.744
2 7.965 9, 088 7,483 6.907 6.686 6.482 5.494 5, 331 4, 040 4t834 4.645 4.699
3 7.935 9.121 7.414 6.989 6.649 6.491 5.500 5.324 4.035 4.834 4.704 4,695
4 7.866 9.237 7.618 7.353 6.875 6.585 5, 530 5.424 3.976 4.722 4.546 4.743
5 7.847 9.179 7.371 7,100 6,784 6.509 5.535 5, 374 3.988 4.843 4.587 4.721
6 7.752 9.247 7.617 7.337 6.871 6.600 5.573 5,401 3.960 4.776 4, 564 4.742
7 7.920 9.149 7.501 7.17S 6.826 6.607 5.515 5, 360 4.003 4.792 4, 598 4.748
8 8.089 8.953 7， 289 6, 614 6.532 6.317 5.462 5.237 4.059 4.922 4.711 4.698
9 8.155 9.005 7.190 6.459 6.565 6.267 5.529 5.238 4.154 4.950 4.819 4.573
10 7.961 9.133 7.362 6.768 6.627 6.418 5.520 5.280 3.987 4.822 4.674 4.734
11 8, 079 8,980 7, 217 6.602 6.493 6, 294 5.504 5.219 4.078 4.921 4.770 A, 712
12 7.953 9.205 7.499 7.173 6.786 6.533 5.535 5, 343 3.977 4.789 4.654 4.664
13 7.986 9, 051 7.360 6.848 6.701 6, 446 5.541 5, 376 4.065 4.822 4.701 4.689
14 8. Ill 9.040 7.210 6.905 6.548 6, 364 5.472 5.341 4.012 4.884 4.725 4.677
15 8.032 9.002 7.325 6.818 6.571 6.369 5.537 5.345 4, 064 4.913 4.672 4.682
16 8.075 8.977 7.199 6,664 6.515 6.349 5.480 5.311 4.068 4.856 4.741 4.703
17 7, 878 9.184 7.502 7.221 6, 793 6.592 5.523 5,405 3.990 4.763 A.588 4, 727
18 7.873 9.165 7.502 7.194 6.775 6.553 5,496 5.384 4.025 4.786 4.629 4.755
19 7.911 9.119 7.574 7.286 6.830 6.507 5.539 5.365 3.949 4.7S1 4.577 4.695
20 8.013 9.186 7.394 6.822 6.634 6.384 5.475 5.297 4.033 4.849 4.657 4.678
21 7.739 8.991 7.232 6.822 6.503 6.260 5.374 5.214 3.941 4.797 4.629 4.591
22 7.887 8, 962 7.178 6.730 6.471 6.259 5, 353 5,215 3, 978 4.754 4, 608 4.566
23 7.921 8.903 7, 258 6.803 6.438 6.285 5.372 5.234 3.973 A.739 4.624 4, 626
24 7.807 8, 900 7.271 6.954 6.526 6.320 5.357 5.283 3.968 4.772 4.585 4.609
25 7.681 9.020 7.359 7.172 6.685 6.440 5.421 5.255 3.922 4.667 4.502 4.665
26 7.892 8.827 7.205 6.728 6.529 6.305 5.363 5.123 3.914 4.847 4.616 4.557
27 8.071 8.730 6.993 6.186 6.159 6.034 5.267 5.150 4, 044 4.919 4.790 4.556
28 7.878 8,771 7, 059 6.452 6.418 6.210 5.318 5, 208 3.962 4.912 4.593 4.613
29 7.803 8.985 7.335 6.934 6.613 6.382 5.<m 5.268 3.892 4.744 4, 510 4.658
30 7.992 8.818 7.126 6,540 6.363 6.167 5.329 5.089 4.028 4.820 4.689 4.548 mean 7.931 9.041 7.339 6.898 6.621 6.400 5, 461 5.295 4.001 4.819 4.644 4.669
S.D. 0.116 0.146 0.163 0.300 0.172 0.148 0, 082 0.091 0.057 0.069 0.078 0, 065
Chrl3 Chrl5 Chrl6 Chrl7 Chrl8 Chrl9 Chr20 Chr21
序号 Chrl4½ Chr22% C rX% ChrYJi
1 3.917 3.294 2.762 2.401 2.380 3.051 1, 120 1.966 1,295 0.920 2, 163 0.126
2 3.743 3.279 2.885 2.595 2,650 2.980 1.431 2.148 1, 296 1.088 2.128 0.123
3 3.770 3.268 2.914 2.565 2.602 3.027 1.434 2.112 1.270 1.073 2.154 0.119
4 3,944 3.254 2.801 2.416 2.402 3.080 1.140 2.000 1,279 0.914 2.170 0.126
5 3.880 3.234 2.866 2.545 2.576 3.049 1.354 2.064 1.275 1, 050 2.154 0.112 6 3, 909 3.322 2.792 2.393 2.389 3.089 1.139 1.966 1.309 0.910 2.217 0.127
7 3.868 3.305 2.833 2.451 2.481 3.002 1.266 2.033 1.327 0, 942 2.165 0.129
8 3.585 3.258 2.863 2.996 1.685 2.262 1.296 1.223 2.091 0.116
9 3.590 3.205 2.926 2.860 2.960 2.879 1.669 2.271 1.339 1, 240 2.046 0. Ill
10 3.714 3.305 2.873 2.677 2.767 3.001 1.539 2.168 1.307 1.154 2.085 0.122
11 3.638 3.303 2.787 2.872 2.947 1, 687 2.244 1.294 1.275 2.041 0.115
12 3, 823 3.285 2.8寸29 2, 493 2.504 3.012 1.323 2.068 1, 293 1.025 2.116 0.120 t—
13 3.734 3.282 2.819 2.670 3.002 1.479 2.100 1.299 1.097 2.135 0.121
14 3.647 3.274 2.894 2, 746 2.758 2, 939 1.501 2.256 1.308 1.147 2.120 0.121
15 3.720 3.262 2.891 2.663 2.730 2.976 1.551 2.176 1, 287 1.164 2.122 0.130
16 3.584 3, 269 2, 970 2.767 2.868 2.965 1.667 2.186 1.270 1.277 2.125 0.114
17 3.864 3.246 2.827 2, 475 2.476 3, 036 1, 228 1, 024 1.297 0, 963 2.277 0, 120
18 3.853 3, 315 2.506 2.489 3.040 1, 229 2.071 1.299 0.963 2.129 o.
19 3.950 3.312 2, 845 2.457 2.484 3.065 1.234 2.031 1.265 0.962 2.136 0.125
20 3.731 3.285 2. m 2.645 2.720 2.906 1.485 2.155 1.299 1.151 2.200 0.114
21 3.690 3.224 2, 829 2.522 2.618 2.919 1.384 2, 119 1.269 1.107 4.224
22 3.642 3.223 2.877 2.578 2.673 2.910 1.422 2， 137 1, 259 1.087 4.226 0.004
23 3.680 3.234 2.824 2.536 2.585 2.968 1.380 2.098 1, 265 1.062 4.187 0.004
24 3.707 3.223 2.809 2.501 2.580 2, 933 1, 372 2, 039 1, 262 1.042 4.177 0.004
25 3.829 3.214 2.740 2.388 2.390 2.978 1, 211 1.975 1.259 0.929 , 293 0.004
26 3.679 3.232 2.859 2.579 2.707 2.918 1, 485 2.081 1, 262 1, 129 4.159 0.003
27 3.404 3.207 2.912 2, 899 2.968 2.893 1.802 2.344 1.263 1.330 4.077 0.003
28 3.525 3.190 2.869 2.771 2.852 2.903 1.647 2.215 1, 278 1.227 4.127 0, 003
29 3.786 3.177 2.778 2.456 2.498 2.980 1.258 2.037 1.271 0.986 4.216 0, 003
30 3.569 3.224 2.913 2.699 2.816 2.898 1, 603 2.168 1.284 1.179 4.133 0.003 mean 3.733 3.257 2.858 2.595 2.644 2.978 1.424 2.117 L 286 1.087 - 一
S. D. 0.135 0.040 0.055 0.147 0.175 0.060 0.187 0.099 0.021 0.120 ― - mean- ChrX/
y½-M 2.139 0.121
S. D.- ChrX/
Y - 0.055 0.006 mean- ChrX/
Y¾-F 4.182 0, 003
S. D. -
ChrX/
Y%-F 0.062 0.001 进一步的， 为考察所得可疑样本中是否存在多半条或多一条染 色体的情形， 选定 X为 50或 100， 计算出对应的染色体参考值 Z (cutoff值） (参见表 2):
Z= ( mean_ChrN%xO.5xX%)/S. D. _ChrN%，其中 N 为常染色体 1~22号, X为 50或 100。
表 2、 判定待测细胞中染色体数量为三体的参考值 Z (cutoff值）
根据计算公式， 计算所测可疑样本中每条染色体的 z score-ChrN。 计算公式如下：
z score.ChrN = (受检样本某染色体的 CtirN%—mean-OirN% ) -61-
0CZ100/010ZN3/X3d £Z£6雕 ΪΟΖ OAV 表 3 可疑样本中对应于每条染色体的 z score-ChrN值
由上述分析可知， 53例受检羊水细胞样本中可疑样本共 8个， 其中对于每个可疑样本的染色体， 共捡出 z score.ChrN 值绝对 值>3的染色体数量异常合计 8个 （参见表 3 ) ，
具体为
1 ) P1的 Chr21， P2的 Chr21， P3的 Chr21和 P4的 Chr21;
2 ) P5的 Chrl8， P6的 Chrl8和 P7的 Chrl 8; 3 ) P8的 Chrl 3。 通过查对表 2 中根据多一条染色体（X-100 ) 计算获得的 Z 值，可以确定其中 P1-P4样本中染色体 21的数量和 P5-P7中染色 体 18的数量比标准细胞中的相应染色体数量多一条， P8中染色 体 13 的数量比标准细胞中的相应染色体数量多半条， 即 P1-P4 为 T21 (唐氏综合症） ， P5- P7为 T18 ( Edwards综合征） ， P8为 镶嵌体的 T13 (镶嵌体的 Patau综合征）。 与传统染色体核型分析 结果完全一致。 实施例 1
另有 6例羊水细胞样本（Q1-Q6 )在经过同样的处理方法和测 序过程后获得上机数据，以实施例 1中 30例标准细胞样本计算出 来的 1^&11_( 11^%和 S. D. -Chr鄉予以计算 z score—ChrN。 6例样 本中检出阳性样本 3例 （参见表 4和表 5 ) 。
表 4 六例待测样本（Q1-Q6 ) 的 ChrN %
Chrl7% 2. 519897 2. 481007 2. 684055 2. 561488 2. 725292 2. 495129
Chrl8% 3. 026389 2. 939323 3. 027751 2. 994205 2. 893438 2. 85936
Chrl9% 1. 291162 1. 240504 1. 79644 1. 317938 1. 482326 1. 235699
Chr20% 2. 096249 2. 063225 2. 174512 2. 076472 2. 173705 2. 049467
Cht21% 1. 290966 1. 267192 1. 297270 1. 291611 1. 896099 1. 297050
Chr22% 0. 992960 0. 989674 1. 125718 1. 017386 1. 164497 0. 995698
ChrX% 2. 173990 4. 134076 2. 117655 2. 177186 4. 104752 4. 197133
ChrY% 0. 116611 0. 003010 0. 003148 0. 001590 0. 003956 0. 002508 表 5 以实施例 1 中 30 例阴性样本的 mean-ChrN°/。和 S. D. _ChrN%计算得出的 6例样本的 z score-ChrN
从结果可以看出， Q5与标准细胞相比多一条 21号染色体， 为 T21; Q3、 Q4与标准细胞相比少一条 X染色体， 为 45X0 (特纳 综合征）。 与传统染色体核型分析结果完全一致。 尽管本发明的具体实施方式已经得到详细的描述， 本领域技 术人员将会理解。 根据已经公开的所有教导， 可以对那些细节进 行各种修改和替换， 这些改变均在本发明的保护范围之内。 本发 明的全部范围由所附权利要求及其任何等同物给出。

Claims (7)

  1. 权 利 要 求
    1. 一种通过测序分析细胞染色体的方法， 其包括以下步驟： a.将获自细胞的 DNA分子随机打断， 得到一定大小的 DNA片 段并测序；
    b.将步驟 a中测定的 DNA序列与人类基因组参考序列进行严 格比对， 得到所测 DNA序列定位于具体染色体的信息；
    c.对应于特定染色体 N, 确定所测 DNA序列中定位于该染色 体上的唯一位置的序列的数量总和， 进而得出对应于该染色体 N 的 ChrN% , 即所测 DNA序列中定位于该染色体 N上的唯一位置的 序列的数量总和占所测 DNA序列中定位于所有染色体上的唯一位 置的序列的数量总和的比例：
    ChrN%=定位于特定染色体 N上的唯一位置的序列的数量总和 / 定位于所有染色体上的唯一位置的序列的数量总和；
    d.将该染色体的 ChrN%与来自标准细胞中对应染色体的 ChrN%进行比较，确定所述细胞的该奈染色体数量与标准细胞是否 存在差异。
  2. 2. 权利要求 1所述的方法， 其中步骤 b中所述严格比对是指 不容错的、 定位于人类基因组参考序列上的唯一位置的比对； 其 中所述的人类基因組序列参考序列是人类基因组中经过屏蔽掉重 复序列后所得到的序列。
  3. 3. 权利要求 1所述的方法， 其中步骤 d中所述细胞样本中特 定染色体数量与标准细胞是否存在差异是通过绘制箱式图方法确 定的，其中将细胞样本中 ChrN%位于超出四分位数差 1. 5倍 - 3倍或 者 3倍以上的四分位数差之间的异常值对应的样本，认定为其染色 体数量与标准细胞存在差异。
  4. 4. 权利要求 1所述的方法， 其中步驟 d中所述细胞的特定染 色体数量与标准细胞是否存在差异是用反映受检细胞样本 ChrN% 与标准细胞样本 ChrN°/。的偏离程度 z scoreXhrN来确定的， 如果 z score— ChrN绝对值大于等于 3， 所述细胞样本对应的特定染色体的 数量与标准细胞相比存在差异。
  5. 5. 权利要求 4所述的方法， 其中
    所述 z score_ChrN= (受检样本特定染色体的 ChrN% ― 特 定染色体 ChrN%平均值 （mean_ChrN% ) ) / ChrN%平均值标准差 ( S. D. _ChrN% ) ；
    其中所述特定染色体 ChrN 均值可根据来自至少 10个， 例 如至少 20个标准细胞样本的该染色体 ChrN%确定；
    其中所述特定染色体 ChrN 均值标准差可根据来自至少 10 个， 例如至少 20个标准细胞样本的该染色体 ChrN 均值确定。
  6. 6. 权利要求 4所述的方法， 其中步裸 d中所述细胞的特定染 色体数量与标准细胞是否存在差异是用 z scoreXhrN与参考值 Z 进行比较加以确定， 其中参考值 Z根据如下方法确定：
    Z = ( mean_ChrN%) χθ. 5xX % / ( 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-Chr 的绝对值大于等于 3且达到参考值 Z的绝对 值时， 所述细胞的特定染色体数量与标准细胞相比存在 X%的差 异。
  7. 7. 权利要求 1所述的方法， 其中所述细胞为羊水细胞， 例如 未培养的羊水细胞或培养过的羊水细胞， 优选为未培养的羊水细 胞。
CN201080068571.9A 2010-08-13 2010-08-13 一种细胞染色体分析方法 Active CN103069004B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2010/001230 WO2012019323A1 (zh) 2010-08-13 2010-08-13 一种细胞染色体分析方法

Publications (2)

Publication Number Publication Date
CN103069004A true CN103069004A (zh) 2013-04-24
CN103069004B CN103069004B (zh) 2015-03-11

Family

ID=45567260

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201080068571.9A Active CN103069004B (zh) 2010-08-13 2010-08-13 一种细胞染色体分析方法

Country Status (7)

Country Link
US (1) US20130210002A1 (zh)
EP (1) EP2604700B1 (zh)
CN (1) CN103069004B (zh)
DK (1) DK2604700T3 (zh)
ES (1) ES2552343T3 (zh)
HK (1) HK1185112A1 (zh)
WO (1) WO2012019323A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104205106A (zh) * 2013-03-28 2014-12-10 深圳华大基因研究院 确定胎儿染色体非整倍性的方法、系统和计算机可读介质
CN108256296B (zh) * 2017-12-29 2021-05-25 北京科迅生物技术有限公司 数据处理装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101675169A (zh) * 2006-06-14 2010-03-17 阿耳特弥斯保健公司 使用样品拆分和dna标签进行稀有细胞分析

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7659054B1 (en) * 2000-05-23 2010-02-09 Nuvelo, Inc. Methods for genetic analysis of DNA to detect sequence variances
JP2004531271A (ja) * 2001-06-22 2004-10-14 ユニバーシティ オブ ジュネーブ 染色体不均衡により引き起こされる疾患を検出する方法
WO2008070249A2 (en) * 2006-10-30 2008-06-12 Cytotrend (Beijing) Biotech Engineering Co., Ltd A method of detecting genomic aberrations for prenatal diagnosis
US20090029377A1 (en) * 2007-07-23 2009-01-29 The Chinese University Of Hong Kong Diagnosing fetal chromosomal aneuploidy using massively parallel genomic sequencing
CA3069081C (en) * 2008-09-20 2023-05-23 The Board Of Trustees Of The Leland Stanford Junior University Noninvasive diagnosis of fetal aneuploidy by sequencing
US8574842B2 (en) * 2009-12-22 2013-11-05 The Board Of Trustees Of The Leland Stanford Junior University Direct molecular diagnosis of fetal aneuploidy

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101675169A (zh) * 2006-06-14 2010-03-17 阿耳特弥斯保健公司 使用样品拆分和dna标签进行稀有细胞分析

Also Published As

Publication number Publication date
US20130210002A1 (en) 2013-08-15
DK2604700T3 (en) 2015-11-02
EP2604700B1 (en) 2015-07-22
EP2604700A4 (en) 2013-10-30
HK1185112A1 (zh) 2014-02-07
WO2012019323A1 (zh) 2012-02-16
CN103069004B (zh) 2015-03-11
ES2552343T3 (es) 2015-11-27
EP2604700A1 (en) 2013-06-19

Similar Documents

Publication Publication Date Title
CN103403183B (zh) 胎儿遗传异常的无创性检测
Yu et al. Genome-wide, single-cell DNA methylomics reveals increased non-CpG methylation during human oocyte maturation
US10741270B2 (en) Size-based analysis of cell-free tumor DNA for classifying level of cancer
Decker et al. Telomere length in Hutchinson-Gilford progeria syndrome
CN104232777B (zh) 同时确定胎儿核酸含量和染色体非整倍性的方法及装置
Martín-Burriel et al. Genetic diversity, structure, and breed relationships in Iberian cattle
CN105349637A (zh) 一种与鸡蛋蛋重相关的snp分子标记及其应用
KR20140140122A (ko) 복제 수 변이를 검측하기 위한 방법 및 시스템
Rohlfs et al. Detecting coevolution through allelic association between physically unlinked loci
CN104830986A (zh) 一种检测胎儿基因信息的方法、装置和系统
CN104169929A (zh) 用于确定胎儿是否存在性染色体数目异常的方法、系统和计算机可读介质
CN104120181A (zh) 对染色体测序结果进行gc校正的方法及装置
AU2012380221A1 (en) Method, system and computer readable medium for determining base information in predetermined area of fetus genome
CN105039571B (zh) 染色体非整倍体(t21、t18、t13)检测试剂盒的质控品及其应用
CN105825076A (zh) 消除常染色体内和染色体间gc偏好的方法及检测系统
CN105821117B (zh) 多倍染色体检测用的标准品产品及其制备方法
Herzog et al. The tissue-specific aspect of genome-wide DNA methylation in newborn and placental tissues: implications for epigenetic epidemiologic studies
CN103069004A (zh) 一种细胞染色体分析方法
CN105765076B (zh) 一种染色体非整倍性检测方法及装置
CN104520437A (zh) 一种染色体非整倍性检测方法及装置
CN101928775A (zh) 一种细胞染色体分析方法
Li et al. Mitogenome evidence shows two radiation events and dispersals of matrilineal ancestry from northern coastal China to the Americas and Japan
CN108546755A (zh) 用于脆性x综合征致病基因检测的校准品及其应用
Zhu et al. The predominance of a specific genotype of Cryptococcus neoformans var. Grubii in China and Japan
Hawkins et al. Genetic reference materials and their application to haematology

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
ASS Succession or assignment of patent right

Owner name: SHENZHEN BEIJING GENOMICS INSTITUTE HEALTH TECHNOL

Free format text: FORMER OWNER: BGI-SHENZHEN CO., LTD.

Effective date: 20130923

Free format text: FORMER OWNER: BGI-SHENZHEN

Effective date: 20130923

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20130923

Address after: 518083 science and Technology Pioneer Park, comprehensive building, Beishan Industrial Zone, Yantian District, Guangdong, 9F-7, 8, Shenzhen

Applicant after: Shenzhen Huada Gene Health Technology Co., Ltd.

Address before: North Road No. 146, building 11F-3 Industrial Zone in Yantian District of Shenzhen city of Guangdong Province in 518083

Applicant before: BGI-Shenzhen Co., Ltd.

Applicant before: BGI-Shenzhen

C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB02 Change of applicant information

Address after: 518083 science and Technology Pioneer Park, comprehensive building, Beishan Industrial Zone, Yantian District, Guangdong, 9F-7, 8, Shenzhen

Applicant after: Shenzhen BGI Medicine Co., Ltd.

Address before: 518083 science and Technology Pioneer Park, comprehensive building, Beishan Industrial Zone, Yantian District, Guangdong, 9F-7, 8, Shenzhen

Applicant before: Shenzhen Huada Gene Health Technology Co., Ltd.

COR Change of bibliographic data

Free format text: CORRECT: APPLICANT; FROM: SHENZHEN BEIJING GENOMICS INSTITUTE HEALTH TECHNOLOGY CO., LTD. TO: SHENZHEN BGI MEDICINE CO., LTD.

C14 Grant of patent or utility model
GR01 Patent grant
C56 Change in the name or address of the patentee

Owner name: SHENZHEN BGI CORPORATION

Free format text: FORMER NAME: SHENZHEN BGI MEDICINE CO., LTD.

CP03 Change of name, title or address

Address after: Yantian District of Shenzhen City, Guangdong province 518083 Hongan street No. 21 China Comprehensive Park 7 Building 7 layer -14 layer

Patentee after: BGI SHENZHEN CO LTD

Address before: 518083 science and Technology Pioneer Park, comprehensive building, Beishan Industrial Zone, Yantian District, Guangdong, 9F-7, 8, Shenzhen

Patentee before: Shenzhen BGI Medicine Co., Ltd.

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160714

Address after: Beishan Industrial Zone Building in Yantian District of Shenzhen city in Guangdong province 518083 1, 9, 3 and 11 floor of Building 1, 5 floor

Patentee after: Shenzhen Huada clinical inspection center Co., Ltd.

Address before: Yantian District of Shenzhen City, Guangdong province 518083 Hongan street No. 21 China Comprehensive Park 7 Building 7 layer -14 layer

Patentee before: BGI SHENZHEN CO LTD

CP03 Change of name, title or address

Address after: 518083 1st, 3rd, 4th and 5th floors of 11 buildings in Beishan Industrial Zone, Yantian Street, Yantian District, Shenzhen City, Guangdong Province

Patentee after: Shenzhen Huada clinical laboratory center

Address before: 518083 Shenzhen Yantian Beishan Industrial Zone Complex Building 1, 9 Floors and 11 1, 3, 5 Floors

Patentee before: Shenzhen Huada clinical inspection center Co., Ltd.

CP03 Change of name, title or address
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 518083 1st, 3rd, 4th and 5th floors of 11 buildings in Beishan Industrial Zone, Yantian Street, Yantian District, Shenzhen City, Guangdong Province

Patentee after: Shenzhen Huada Medical Laboratory

Address before: 518083 1st, 3rd, 4th and 5th floors of 11 buildings in Beishan Industrial Zone, Yantian Street, Yantian District, Shenzhen City, Guangdong Province

Patentee before: SHENZHEN HUADA CLINIC EXAMINATION CENTER