system: [SYSTEM]

user: 3. 系统和方法 本文建立的系统包括3个核心模块：特征统计模块、特征分析模块和可视化模块。处理过程为把若干段等长的序列分别进行特征统计，再从统计结果中选出具有代表性的特征，最后对选取得特征集合进行可视化，显示各序列中特征片段的分布特征。该系统的结果和处理过程如图1所示。 3.1. 数据处理模块 在诸多对DNA数据分析测量的研究中，使用了诸如数据挖掘，基于统计，基于概率等的各种方法来进行探索，而在对基因回文结构的研究探索中，期待能够发现回文结构并对其位置以及序列形式进行记录和统计分析，以得出此类结构是否存在某种特殊现象及分布规律，故而选择基于统计的方式来探索其 Table 1 Four kinds of corresponding relations of base fragmen 片段本身 A C G T C A 反向片段 A C T G C A 互补片段 T G C A G T 反向互补片段 T G A C G T 表1. 碱基片段的四种对应关系 Table 2 Conjugate symmetry of base pai 片段本身 互补片段 反向互补片段 A A G C T T T T C G A A A A G C T T C C T A G G G G A T C C C C T A G G A C T A G T T G A T C A A C T A G T 表2. 碱基片段的共轭对称关系 图1. 系统模型 特性。通过对对应序列之间的匹配，寻找两条序列之间是否存在由同种类型的碱基片段形成的回文结构，并观察其在同一数据集中的特征。最后用统一的数据结构保存得到的信息，为统计出现在数据集中的位置分布情况以及出现次数以及可视化提供数据。 本文中的匹配模式，即以数据集中的单条或成对DNA序列作为处理对象，在序列中通过滑动窗口的形式依次选择碱基片段，对于序列中选取碱基片段位置之后出现的对应形式碱基片段记录位置，从而实现匹配处理，得到关于数据集中每个出现碱基片段的所有形式的位置数据。匹配过程如图2所示。 模块的输入： 以FASTA格式保存记录的两组若干等数量等长度DNA序列集合，两组序列皆取自某基因左右两端各500 bp长DNA序列。每条编号的序列皆有另一条编号相同的序列与其对应。 模块的匹配处理： 从两组数据集中依次选取左右两条对应编号的DNA序列，使左右两条序列各自对自身序列进行匹配，记录片段位置信息，同时将左侧序列与右侧序列进行匹配，初始匹配片段长度从3开始，直至出现长度为n的碱基片段使得该序列中不存在有任意长度为n + 1的片段在对应序列中存在与其匹配形式的碱基片段。 以如下序列处理过程为例： 如步骤①所示，初始匹配片段长度为3，选取前三个碱基TAC作为匹配片段，由前的匹配方法可知对应四中片段为TAC，ATG，CAT，GTA。 通过往后检索，即可匹配到在右侧序列位置5、位置11、位置20出现了反向片段本身CTT，未出现其它对应形式片段。 ① 左：T A CT C G A C G AC T TCT A C G A C T A A G G C G T A T C C T C G A G 右：G C G TCATAAACATCTT G ACCATTGTGAATTCTAGAA C 记录片段TAC对应4种匹配形式的位置后，将匹配片段位置向下移一位，进行下一个片段的匹配，如步骤②所示，移至ACT片段处。 ② 左：T A C T C G A C G AC T TCT A C G A C T A A G G C G T A T C C T C G A G 右：G C G TCATAAACATCTT G ACCATTGTGAATTCTAGAA C 如步骤③所示，当滑动窗口内选中的片段ACT已经在对应序列进行匹配处理过时，跳过该片段，滑动窗口向下移一位，移至CTT片段处。 ③ 左：T A C T C G A C G AC T TCT A C G A C T A A G G C G T A T C C T C G A G 右：G C G TCATAAACATCTT G ACCATTGTGAATTCTAGAA C 当长度为3的片段匹配到序列最后一位时，长度加1，从序列前段开始再次匹配。如步骤4所示。 图2. 匹配模式模型 ④ 左：T A C T C G A C G AC T TCT A C G A C T A A G G C G T A T C C T C G A G 右：G C G TCATAAACATCTT G ACCATTGTGAATTCTAGAA C 当长度n等于序列长度的二分之一时，停止匹配。进行下一个序列的匹配处理。 直到数据集中的最后一对序列处理完毕后，该数据集的匹配处理过程结束。 模块的输出： 数据集所对应的结果数据集中，保存左侧序列自匹配结果，右侧序列自匹配结果以及左右对应异匹配结果所有有相关记录存在的序列的结果数据，每条序列的结果数据中通过匹配片段长度的不同将结果数据分离组织，不同长度匹配片段结果数据集中分别保存该长度下有匹配位置记录的片段结果数据。 处理匹配结果及分析： 数据集1的数据为read序列在H2H基因上map到两个位置，且位置相距2万bp以上的对称长度500 bp的左右序列结果如表3所示，数据集2的数据为随机在人类基因组上选取的相隔2万bp以上的对称长度500bp的左右序列，结果如表4所示。 在实际的处理过程当中，采用了多组数据进行分析处理，在通过对各组数据最终结果进行分析和比较后发现： (1) 碱基片段“AAGCTT”在进行处理的数组DNA序列中出现频率较高，其反向互补是其本身，存在着共轭对称的情况，同时在生物学中是限制酶的切割位点。 (2) 对于所有数据集中的序列，每条序列中都存在长度至少为8或以上的片段有出现自身，互补，反向，反向互补四种形式之一的片段与其对应。同时可以发现存在对应片段情况的碱基片段平均长度为11到13，可以推测如果出现回文结构发夹结构等现象互补长度不会太长。 Table 3 Processing result of dataset 数据集中序列数量 3000对 最长片段长度 165 bp 最短片段长度 8 bp 平均片段长度 12 bp 出现最多的片段 AAGCTT：1492对序列中出现 CCATGG：1184对序列中出现 ATTTTT：745对序列中出现 表3. 数据集1处理结果表 Table 4 Processing result of dataset 数据集中序列数量 3000对 最长片段长度 76 bp 最短片段长度 8 bp 平均片段长度 12 bp 出现最多的片段 AAGCTT：1400对序列中出现 CCATGG：1036对序列中出现 TTTTTT：711对序列中出现 表4. 数据集2处理结果表 (3) 在数据集1和数据集2中，可以发现得到的统计结果极其相似，而出现次数最多的AAGCTT以及CCATGG都是共轭对称同时它们都是某种限制酶。 3.2. 可视化分析方法 在获得片段位置相关数据并进行统计分析后，得到了每组数据中片段的出现情况，并重点观察分析了出现次数最多的片段，为了进一步对序列中有对应形式的片段位置的确认以及进行观察，需要建立一套能够更为直观地对其中的精细特征进行观察的可视化模型。模型的结构图如图3所示。 输入：匹配模块输出的位置结果及选定需可视化的片段； 处理过程：在输入数据集片段位置结果以及选定片段后，首先需要将所有选定长度的片段位置结果提取出来，再从中选定指定的片段对应形式的位置，将其可视化。 而在匹配过程中，以X轴表示序列长度，Y轴表示不同序列，如数据集左右各N条序列，在第n组左右对应序列中，其长度为L，片段M在自匹配左侧序列中位置a，b (0 < a， b < L)存在对应形式，在异匹配右侧序列中位置c，d (0 < c， d < L)存在有对应形式片段，则在坐标(n， −a)，(n， −b)，(n， c)，(n， d)作标志记录，当数据集中所有序列都进行过处理和显示后，生成图示，处理结束。 输出：针对自匹配所对应的片段位置图示以及自匹配和异匹配位置数据对应的片段对比位置图。 3.3. 可视化结果及分析 在获得片段位置相关数据并进行统计分析后，得到了每组数据中片段的出现情况，并重点观察分析了出现次数最多的片段，如图4至图9所示。 通过上述可视化图示，可以观察出下列现象及特征： (1) 由数据集1中的图示可见，在数据集1中ATTTTT片段位置在序列中分布较为随机和均匀，而出现次数最多的AAGCTT以及其次的CCATGG在分布上具有十分明显的特征，它们在分布上自匹配和异 图3. 可视化模型结构 图4. 数据集1中AAGCTT片段匹配结果图示 图5. 数据集1中CCATGG片段匹配结果图示 图6. 数据集1中ATTTTT片段匹配结果图示 图7. 数据集2中AAGCTT片段异匹配结果图示 图8. 数据集2中CCATGG片段匹配结果图示 图9. 数据集2中TTTTTT片段匹配结果图示 匹配都有很明显的分界线，在序列大概50 bp左右极少出现对应形式片段，大多出现在50~250 bp区间段，之后成彗星尾部状出现次数及频率逐渐减少。 (2) 从数据集2中的图示来看，呈现出了与数据集1中相似的特征与情况，AAGCTT与CCATGG分布十分相似，而TTTTTT则无明显现象。 (3) 通过对数据集1和数据集2的对比中，两组数据呈现了相似的分布现象，出现次数最多的AAGCTT和CCATGG在另一组数据集中分布相似，说明了该特征不仅仅只局限于某一数据集。 (4) 在对选取的其它数据集如部分编码区与非编码区的分析可视化中，并未发现较明显的类似现象与特征，这说明数据集1和数据集2所在的H2H基因可能具有一定的特殊性。