生物信息学辅助定位及延伸辐射诱导未知表达序列标签

研究辐射诱导的基因表达调控对于认识细胞对辐射损伤的应激反应有重要意义.在低剂量辐射诱导新基因RIG1表达序列标签(expression sequence tag,EST)片段的基础上,通过非克隆cDNA文库和RACE(rapidamplification of cDNA end)技术获得了其3′末端.依据实验得到的这两段EST序列所提供的信息,通过生物信息学分析将RIG1基因初步定位在20号染色体.对20号染色体RIG1区基因组序列进行外显子扫描,发现预测的外显子正好与实验得到的EST相吻合.利用预测的外显子设计特异引物,成功地克隆了RIG1基因全长序列.同时,对20号染色体RIG1区的生物信息学分析表明,在RIG1基因的上游存在启动子区,从而确定了RIG1基因的基因组序列.因此,通过生物信息学辅助设计实验,快捷地定位及延伸了未知EST片段RIG1,基本完成了RIG1的全基因、基因组序列及染色体定位研究.     

生物信息学研究新基因LRP15

为了利用生物信息学探索网上克隆基因与预测基因功能的新方法,首先用一段1.8kbDNA片段的人类EST数据库中进行电子杂交并对相互重叠的EST片段组装,通过引物设计进行cDNA末端快速扩增,以高通量基因组序列（HTGS）数据库及SAGE库为基础,进行染色体定位及组织表达分析。通过DAS及RPS－BLAST程序预测其蛋白质结构及功能。结果显示,LRP15 cDNA全长1718bp,含一个780bp的开放读码框架,编码259个氨工酸,该基因定位于染色体3p24,在多种组织中表达。其编码蛋白含有N端富含亮氨酸重复序列及潜在的穿膜序列。研究表明,生物信息学是克隆与预测基因功能的有效方法;LRP15基因可能是一个参与细胞发育调控的新基因。     

小鼠睾丸特异表达基因TSEG-1的克隆及序列分析

从表达序列标签(expressed sequence tags, ESTs)数据库ZooDDD中获得小鼠正常睾丸表达的EST, 通过dbEST数据库检索出与其高度同源的EST序列, 构建EST叠加群(contigs), Biolign软件拼接, GeneScan软件预测contigs对应的基因组序列中的外显子、内含子; 针对开放阅读框设计引物序列, 采用RT-PCR从小鼠睾丸组织中克隆新基因的cDNA, 分析该基因在小鼠各脏器中的mRNA表达, 并对测序结果进行生物信息学分析。结果表明: 在小鼠X染色体的1 668~2 011 kb间克隆出一新基因TSEG-1, 全长为510 bp, 开放阅读框为336 bp, 编码111氨基酸, 分子量12.84258 kDa, 等电点11.4000。RT-PCR证实该基因开放阅读框正确, 在小鼠睾丸组织中特异性表达, 且与小鼠其他cDNA 无同源性, 获得GenBank 登录号EU079024。功能区分析发现TSEG-1蛋白可能为一种跨膜蛋白, 跨膜区位于第41~61氨基酸残基。TSEG-1基因与人类睾丸特异性组蛋白2a变异体基因有较高同源性, 在TSEG-1基因5′-端非编码侧翼预测发现存在1个启动子区域, 范围为680 bp。 TSEG-1蛋白可能有4个抗原性位点, 2个特异性蛋白激酶的磷酸化位点, 其亚细胞定位可能位于线粒体。小鼠睾丸特异性基因TSEG-1的克隆为进一步研究其生物学功能和表达调控奠定了基础。     

EST分析技术在鸡基因组学研究中的应用

表达序列标签（EST）是由大量随机取出的cDNA库克隆经测序得到的组织或细胞基因组的一段cDNA序列,一个EST代表生物体某种组织某一时期的一个表达基因。综述了EST分析技术在鸡基因组研究中的应用。如用于鉴定、发现和预测鸡的新基因,用于基因图谱的绘制,用于筛选基因的单核苷酸多态性（SNP）位点,用于基因表达分析和基因芯片制作等。EST数据库和生物信息学的联合分析技术在推动家鸡后基因组的研究中发挥着重要的作用。     

通过新基因计算机识别与实验确认对NCBI人类基因数据库一些模式参考序列错误的分析与纠正

采用生物信息学分析与实验确认相结合的技术路线,通过所识别的基因在非冗余数据库比对发现了网上公布的计算机注释人类基因组编码序列存在各种类型的多处错误,包括cDNA水平的一个或一段碱基插入、缺失或突变,或是这些错误的不同排列组合,其中以错误插入为多,往往导致编码氨基酸的移码突变。最先举证了NCBIGENOME Annotation Project预测人类新基因的下列错误类型：(1)开放读码框架(0RF)中错误插入一个碱基造成编码氨基酸移码;(2)错误拼接;(3)开放读框中错误插入一个或一段碱基造成该读框提前终止。只编码N-端氨基酸的cDNA序列而不完整;(4)只有编码c一端氨基酸序列的cDNA而不完整;(5)只是正确基因0RF中间的一段编码蛋白cDNA序列而不完整,缺N-端与C-端氨基酸序列,并且将不完整蛋白氨基酸序列的第一个非起始码氨基酸错误地预测为起始码氨基酸,如将L错误地预测为M;(6)开放读框中错误插入一个或一段碱基造成前面出现不该有的终止码,因而编码蛋白缺开头部分氨基酸;(7)可能将污染基因组序列当作完整基因cDNA序列对待而预测出所谓单一外显子基因。即便真是基因,也只是较长单一外显子mRNA中有一小0RF,而0RF起始码上游同一相位确实存在终止码,无其他特点符合基因条件;(8)所预测基因只有0RF,而0RF两端没有任何EST证据,可据此0RF拼接出受EST和人类基因组双重支持的完整基因cDNA(开放读框上游同一相位有终止码),预示所预测0RF参考序列可能不正确;(9)有EST实验证据支持存在基因的人类基因组序列范围内又被预测出一条相似但更小的蛋白编码基因,因而新预测基因有可能是错误的。     

染色体7q32-ter鼻咽癌相关基因的初步研究

为克隆7q32-ter区域等位基因杂合性丢失最小共同缺失区的鼻咽癌相关基因.以STS D7S509为探针,PCR法筛选位于该区的细菌人工染色体(BAC)克隆,应用EST介导的定位-候选克隆策略并结合生物信息学筛选出在鼻咽癌细胞株和活检组织中表达增强的EST AA773454,cDNA克隆测序和生物信息学资源获取全长cDNA,DNA印迹和甲基化分析研究其表达增强的机制.结果表明,克隆的NAG18基因cDNA全长802 bp,编码227个氨基酸,定位于胞核.该基因分别与人、鼠TAXREB107基因及RPL6基因高度同源.其表达增强的机制不是基因拷贝数的丢失和甲基化位点的改变.可以断定NAG18是定位于7q32-ter最小共同缺失区的鼻咽癌相关基因,它是一高度保守的基因,参与了DNA的转录活化.     

多发性骨髓瘤患者13q14.3区域一个候选抑瘤基因MYETS1的克隆与序列分析

为克隆位于多发性骨髓瘤(multiple myeloma,MM)患者染色体13q14.2～13q21.1区域候选抑瘤基因,通过生物信息学分析获取疾病基因定位区域内代表新基因的ESTs,并运用半定量RT-PCR检测它们在正常人与MM患者骨髓组织中的表达水平,发现一条在MM患者骨髓组织中明显表达下调的EST(GenBank收录号:H86826).Northern印迹杂交显示H86826在骨髓组织中转录本大小为1.5kb.通过购买商品化克隆IM-AGE223589测序获得了H86826所代表的基因的1491 bp全长cDNA序列(GenBank收录号:AY368652),人类基因组命名委员会将其命名为MYETS1(myeloma tumor suppressor 1).生物信息学分析其为一个编码分子质量为15.1 kD、等电点为6.13的135个氨基酸的新基因.该基因的功能正在进一步的研究之中.     

PAC579基因组序列外显子计算机分析与预测

为了寻找肝癌相关基因,选择位于人原发性肝细胞肝癌染色体17p13.3杂合性缺失最小共同缺失区内含有D17S926位点的PAC579克隆进行了基因组测序,通过实验分离和克隆了位于PAC579上的4个全长新基因,并将它们在该基因组上进行了外显子定位.使用5种计算机预测外显子的方法和4种剪接位点识别的方法对PAC579前60 kb序列(含上述4个基因)进行扫描,将得到的外显子预测结果和剪切位点结果与实验结果进行了对照与评析,又对未获得基因的区域加以预测.     

绵羊基因组MHC ClassⅢ区段部分基因的分离与鉴定

绵羊基因组MHC段分为ClassⅠ、Class Ⅲ和ClassⅡ（含Ⅱa和Ⅱb两个亚区）3个区段,与另外2个区段相比,Class Ⅲ区的基因信息远少于ClassⅠ和ClassⅡ区.为丰富绵羊基因组 MHC Class Ⅲ 区段基因信息. 本研究用位于中国美利奴羊基因组BAC文库中 MHC ClassⅢ 区段4个BAC克隆的酶切片段制备32P 标记探针,继而采用噬菌斑原位杂交筛选法筛选中国美利奴羊混合组织cDNA 文库,并对分离到的cDNA 阳性克隆进行全序列测定及生物信息学分析.本实验共筛选出 31 个 cDNA 阳性克隆,对其序列进行了测定及分析,确定了序列在ClassⅢ 区段上的位置,并通过在NCBI中的同源检索对其功能做了初步鉴定.由实验可知,利用BAC 文库与 cDNA 文库杂交筛选法对较大区段基因的筛选和分离是有效的.同时,对分离到的表达序列结合生物信息学进行分析,这将有助于对该序列功能的深入研究.     

以绵羊MHC区段BAC克隆酶切片段为探针杂交筛选绵羊混合组织cDNA文库

在已知中国美利奴羊MHC(Major histocompatibility complex)区段BAC(Bacterial artificial chromosome)克隆序列信息和预测的基因注释前提下,用位于中国美利奴羊基因组BAC文库MHC区段的6个BAC克隆酶切片段为探针,以噬菌斑原位杂交筛选法筛选中国美利奴羊混合组织cDNA文库(库库杂交),对分离到的cDNA阳性克隆进行全序列测定,并与相应的已知序列信息和基因注释的BAC克隆比对以及在NCBI Blastn数据库中序列相似性检索,旨在验证基因注释结果的准确性和对基因(序列)功能的初步分析。实验中,经过两轮杂交共筛选出27个cDNA阳性克隆(序列),并发现这些序列均可定位到相应的BAC克隆上,且25条序列处在注释基因的外显子部分;在NCBI数据库中经Blastn序列相似性检索发现,23条序列与牛基因的序列相似性最高,且与免疫功能密切相关。     

Retroposition in a family of carcinoma-associated antigen genes.

The gene encoding the carcinoma-associated antigen defined by the monoclonal antibody GA733 is a member of a family of at least two type I membrane proteins. This study describes the mechanism of evolution of the GA733-1 and GA733-2 genes. A full-length cDNA clone for GA733-1 was obtained by screening a human placental library with a genomic DNA probe. Comparative analysis of the cDNA sequence with the previously determined genomic sequence confirmed that GA733-1 is an intronless gene. The GA733-2 gene encoding the monoclonal antibody-defined antigen was molecularly cloned with a cDNA probe and partially sequenced. Comparison of GA733-2 gene sequences with the previously established cDNA sequence revealed that this gene consists of nine exons. The putative promoter regions of the GA733-1 and GA733-2 genes are unrelated. These findings suggest that the GA733-1 gene was formed by the retroposition of the GA733-2 gene via an mRNA intermediate. Prior to retroposition, the GA733-2 gene had been affected by exon shuffling. Analysis of GA733-2 exons revealed that many delineate structural motifs. The GA733-1 retroposon was localized either to chromosome region 1p32-1p31 or to 1p13-1q12, and the GA733-2 founder gene was localized to chromosome 4q.     

黄龙病菌侵染过程中柑橘CsMYB基因克隆及表达分析

MYB类蛋白是一类与植物防卫反应有关的转录因子家族,本研究利用构建的甜橙健株与感染黄龙病的病株差减（SSH）文库,采用RACE技术克隆了一个MYB类基因的cDNA全长序列,命名为CsMYB（GenBank登录号为HQ841074）。柑橘CsMYB基因的cDNA全长为1 306 bp,生物信息学分析显示,该基因包括一个909 bp的完整开放读码框以及一个典型的26 bp poly-A,编码302个氨基酸,分子量为32.97 kD,等电点为8.5。同时,还有MYB类基因的保守特征区域,即在N端有两个典型的MYB DNA结合域：R2和R3。实时荧光定量PCR分析表明,柑橘CsMYB基因受到黄龙病菌侵染后的不同时期表达量不同,伴随着黄龙病病程的变化呈现不同的表达变化。推测CsMYB基因是一个转录因子,可能参与柑橘对黄龙病菌的防御反应过程。     

Cloning of the cDNA encoding phenylalanyl tRNA synthetase regulatory alpha-subunit-like protein whose expression is down-regulated during differentiation.

Hybrid polar compounds (HPCs), such as suberoylanilide hydroxamic acid (SAHA), induce differentiation of transformed cells. Differential display of RNA was used to identify genes whose expression is changed during SAHA-induced differentiation of murine erythroleukemia (MEL) cells. One such cDNA was identified whose mRNA level decreased by 50% after 8h of SAHA treatment as determined by Northern blot analysis. The full-length cDNA (1944bp in length) was cloned by sequencing of an EST clone and rapid amplification of 5' cDNA ends (5'-RACE). The predicted amino acid sequence is 589 amino acids and shares 45% identity with the yeast cytoplasmic phenylalanyl tRNA synthetase (PheRS) regulatory alpha-subunit. Human EST clones which share over 90% identity of predicted amino acid sequence with this cDNA map to chromosome 2 near the paired box homeotic gene 3 (PAX3) locus, a region syngenic to mouse chromosome 1. This is the first report of the cloning of the full-length cDNA for the murine PheRS regulatory alpha-subunit-like protein. The level of PheRS alpha-subunit-like mRNA is regulated during differentiation but not during cell cycle progression.     

多发性骨髓瘤细胞中一个新的Alu超基因家族成员的克隆与分析

根据Genebank收录的多发性骨髓瘤细胞株（ARH－77）表达上调EST AF497797设计引物,采用半定量RT－PCR证实了多发性骨髓瘤患者及正常人骨髓细胞中该expressed sequence tages(EST)存在表达差异。用EST AF497797作探针筛选人胚肾cDNA文库,获得cDNA克隆经测序并用生物信息学方法对该序列进行了初步分析。AF497797在多发性骨髓瘤患者骨髓中确有较高的表达,而在正常人骨髓细胞中低表达。获得的全长cDNA克隆序列长1248bp(Genebank登录号：AY094612）。生物信息学分析显示该片段全长cDNA编码44个氨基酸的蛋白产物且可能属于Alu家族成员。基因AY094612为一个在多发性骨髓瘤中表达上调的新基因,其改变可能与多发性骨髓瘤的发生与发展有关。     

Structural and functional analysis of a growth-regulated gene, the human calcyclin

Calcyclin was originally defined as a cDNA clone (2A9) whose cognate RNA is growth-regulated and whose sequence shows strong similarities to the sequences of the S-100 protein, a calcium-binding protein, as well as to a subunit of the major cellular substrate for tyrosine kinase. Using the full-length cDNA, we have now isolated from a human genomic library several phages containing calcyclin sequences. One of the phages, ch. 28-10, contains the entire calcyclin gene, plus extensive flanking sequences. The calcyclin gene is a unique copy gene and has 3 exons. The 5' flanking sequence has been characterized, both structurally and functionally. Besides a TATA box, it contains, in the region proximate to the cap site, GC boxes and a sequence with a strong homology to the enhancer core of the SV40 promoter. Other enhancer-like elements are found scattered in both the 5' and 3' flanking regions. The proximate 5' flanking region is very active in driving the transient expression of linked reporters in transfection experiments. Finally, the calcyclin gene has been localized to the long arm of human chromosome 1, near the ski oncogene.     

Cloning of the human cDNA which can complement the defect of the yeast mannosyltransferase I-deficient mutant alg 1

The assembly of the lipid-linked oligosaccharide, Glc(3)Man(9)GlcNAc(2)-P-P-Dol, occurs on the rough ER membrane in an ordered stepwise manner. The process is highly conserved among eukaryotes. In order to isolate the human mannosyltransferase I (MT-I) gene involved in the process, we used the Saccharomyces cerevisiae MT-I gene ( ALG1 ), which has already been cloned. On searching the EST database with the amino acid sequence of the ALG1 gene product, we detected seven related human EST clones. A human fetal brain cDNA library was screened by PCR using gene-specific primers based on the EST nucleotide sequences and a 430 bp cDNA fragment was amplified. The cDNA library was rescreened with this 430 bp cDNA, and two cDNA clones (HR1-3 and HR1-4) were isolated and sequenced. On a homology search of the EST database with the nucleotide sequence of HR1-3, we detected a novel human EST clone, AA675921 (GenBank accession number). Based on the nucleotide sequences of AA675921 and HR1-4, we designed gene-specific PCR primers, which allowed to amplify a 1.8 kb cDNA from human fetal brain cDNA. This cDNA was cloned and shown to contain an ORF encoding a protein of 464 amino acids. We designated this ORF as Hmat-1. The amino acid sequence deduced from the Hmat-1 gene showed several highly conserved regions shared with the yeast and nematode MT-I sequences. Furthermore, this 1.8 kb cDNA successfully complemented the S. cerevisiae alg1-1 mutation, indicating that the Hmat-1 gene encodes the human MT-I and that the function of this enzyme was conserved between yeast and human.