中国缘蝽科部分种类细胞色素b基因序列及系统进化研究

以线粒体细胞色素b(Cyt b)基因作为分子标记,首次对缘蝽科4亚科14种昆虫进行序列测定,获得Cyt b基因412 bp的序列片段,该片段中碱基T、C、A、G的平均含量分别为34.2%、11.4%、35.7 %和18.7 %,A T平均含量为69.9 %,明显高于G C含量(30.1 %).密码子第3位点的A T含量高达82.8%,亚科间序列变异大,有193个核苷酸位点发生变异,碱基替换多发生于第3位点.以筛豆龟蝽为外群构建系统发育树,表明在亚科级关系上,姬缘蝽亚科最原始,蛛缘蝽亚科次之,巨缘蝽亚科和缘蝽亚科亲缘关系较近,为较进化种类.     

番石榴实蝇和瓜实蝇的线粒体L-rRNA基因测序及双翅目昆虫L-rRNA基因的结构特征与系统发育关系分析

RNA二级结构比一级结构包含更多的信息, 在物种系统发育分析中更能反应真实情况。因而, 本研究对番石榴实蝇Bactrocera correcta和瓜实蝇B. cucurbitae的L-rRNA全基因进行了测序, 构建了双翅目L-rRNA基因二级结构的模式图, 并分析了其结构特征, 而后基于H45~H47茎环结构参数和结构序列及L-rRNA结构序列对双翅目13科的系统发育关系进行了初步分析。结果表明： 双翅目昆虫具有保守的L-rRNA二级结构; 不同结构区碱基组成和分布不均匀, Ⅳ区和Ⅴ区的完全保守碱基含量最高, GC含量最高, 绝大部分科Ⅵ区的AT斜率均小于0; 科内特有碱基和科间少数保守碱基大部分为G或C。瘿蚊科与双翅目其他科具有较远的系统发育关系, 而丽蝇科、 寄蝇科和食蚜蝇科有相近的系统发育关系; 虻科与网翅虻科同在一个小分支; 蠓科与蚊科同在一个大的系统发育分支。利用单方面参数不易得到理想的进化结果, 准确系统发育分析需要结合多方面的参数。     

赤潮棕囊藻(Phaeocystis globosa)rDNA ITS区序列变异与二级结构分析

测定了南海球形棕囊藻香港株P1、P2和湛江株ZhJ1的rDNAITS区序列(含5.8srDNA),结合Gen Bank的13条同源序列,比对长度为904bp,变异位点271个,简约信息位点221个,平均(A+T)(34.5%)<(G+C)(65.4%).藻株P1、P2和ZhJ1序列存在变异位点20个,序列间相似性为97.9%～98.5%.ITS序列在种间和种内的解析度高于18srDNA和28srDNA基因;构建的NJ树、MP树、贝叶斯推断系统树的结构是一致的,不同种类的棕囊藻单独聚类,不同地理来源的球形棕囊藻混杂分布但相同地理来源的藻株多聚类在一起.RNA二级结构显示,不同藻种间5.8srDNA区结构基本一致,表现出属的特异性;ITS1、2区结构表现较大的种间差异,表明ITS区RNA二级结构可为棕囊藻分类鉴定提供有用的分子结构信息.     

基于COⅠ部分序列对皿蛛科盖蛛属(Linyphiidae,Neriene)几种蜘蛛亲缘关系的研究

将自测的皿蛛科Linyphiidae盖蛛属Netiene 15种蜘蛛与从GenBank中下载的1种盖蛛的COⅠmtDNA序列片段进行了同源性比较,以皿蛛科朱氏盾蛛Frontinella zhui的作为外群,用NJ、MP、ML和贝叶斯法构建了分子系统树。在所比较的818bp的序列中,有371个变异位点,143个简约信息位点;A、T、G、C的平均含量为23．1％、42．9％、17．8％、16．2％,所有种类的A＋T含量比较高,平均为66．0％。分子系统树基本支持已有的形态学研究,但鹤嘴盖蛛Nerien macella却没有被归入哈氏盖蛛组,可能是由于鹤嘴盖蛛与同组其他3种蜘蛛属于不同的动物地理区系所致;分析结果同时支持将中华盖蛛Neriene sinensis和黑斑盖蛛Neriene nigripectoris归入盾形盖蛛组。     

基于线粒体16S rRNA基因对中国狼蛛科(Lycosidae)系统发育研究

将自测的中国狼蛛科Lycosidae4亚科6属26种和从GenBank中检索到的北美2种豹蛛的mtD-NA-16S rRNA序列进行比较;以漏斗蛛科1种蜘蛛作为外群,对碱基序列的组成和遗传距离进行了分析,采用Bayesian方法和最大简约法(MP)构建分子系统树。研究结果表明:16SrRNA基因的部分序列为340bp到360bp,A T含量平均为75%,存在较强的A T含量偏向性;序列共有157个碱基存在变异,其中79个简约信息位点。狼蛛科各属间的遗传距离介于0.026￣0.200之间。2种建树方法均表明:科内的属及属内的种优先聚在一起;水狼蛛属相对马蛛属是狼蛛科中较为原始的类群,分化较早;獾蛛属作为1个单系群与熊蛛属合为1个并系,属于狼蛛亚科。狼蛛科6属间的分子系统关系为(Pirata(Hippasa(Trochsa Arctosa(Pardosa Wadicosa))))。     

2008～2009年珠海市H3N2亚型流感病毒HA1基因特性分析

为了解2008~2009年珠海市H3N2亚型流感病毒HA1基因变异情况,选择珠海市2008~2009年期间不同时间点的经狗肾传代细胞(MDCK)培养分离的H3N2亚型流感毒株20株,提取病毒RNA,通过RT-PCR扩增HA1基因片段,将产物纯化并测序,推导氨基酸序列,进行基因进化特性分析。与同时期的疫苗株比较,2008年珠海市流行的H3N2亚型流感毒株HA1区抗原决定簇的氨基酸位点变异数少于4个;2009年珠海市流行的H3N2亚型流感毒株除09-0056外,HA1区存在5个位于抗原决定簇内的变异氨基酸位点。2008年H3N2亚型流感毒株的HA1区的糖基化位点与疫苗株一致;2009年H3N2亚型流感毒株HA1区丢失第144位糖基化位点。2008~2009年H3N2亚型流感毒株RBS氨基酸序列未见明显变异。与2008年H3N2亚型流感毒株比较,2009年H3N2亚型流感毒株HA1区抗原决定簇内存在多个位点的氨基酸替换。这些说明2008年珠海市流行的H3N2亚型流感病毒不是新变种;2009年流行的H3N2亚型流感病毒为新的变异株,这可能是H3N2亚型流感病毒在2009年6-9月为珠海地区季节性流感流行优势株的原因。     

哺乳动物cd59基因的进化

采用PCR以及BLAST方法从5种哺乳动物中获得了cd59基因的编码区序列, 结合 GenBank中已有的序列, 计算cd59基因在哺乳动物中的核苷酸替换速率. 对非同义替换速率和同义替换速率进行比较的结果显示, cd59在哺乳动物中总体上受到负选择作用; 用PAML软件“位点-特异”模型检测到4个受到正选择作用的位点, 4个位点分布于分子表面, 其中2个位于功能重要的区域; 此外, 用“支-位点-特异”模型在小鼠通过基因复制后形成的cd59a和cd59b上检测正选择引起的加速进化, 并检测到该支系特异的正选择位点1个.     

组蛋白H3Ser10磷酸化在草地贪夜蛾Sf9细胞有丝分裂中的动态分布

【目的】探讨表观遗传标记组蛋白H3Ser10磷酸化在草地贪夜蛾Spodoptera frugiperda Sf9细胞系细胞有丝分裂中的功能。【方法】序列比对分析组蛋白H3保守性。通过固相法合成特定位点磷酸化修饰的一段组蛋白H3的肽段【RK(p S)TGGKAPRKQLC】进行抗体制备;培养Sf9细胞,通过细胞爬片制备有丝分裂的片子,统计不同时期的细胞数目;以免疫荧光标记检测组蛋白H3Ser10磷酸化抗体在不同时期的定位特点。【结果】序列比对分析发现组蛋白H3第1-60位氨基酸在大多数物种中高度保守。在草地贪夜蛾Sf9细胞中,组蛋白H3Ser10磷酸化发生在早前期细胞的核膜附近,成点状分布;随着细胞周期的进行,早中期达到最高水平,在整条染色体上集中分布。有丝分裂后期开始去磷酸化,并于末期完成去磷酸化。【结论】组蛋白H3Ser10磷酸化与Sf9细胞有丝分裂中染色质的凝缩相关。     

三种蜘蛛酯酶同工酶的比较研究

应用聚丙烯酰胺凝胶电泳对狼蛛科的拟水狼蛛、蟹蛛科的三突花蛛和肖蛸科的鳞纹肖蛸3种蜘蛛的酯酶同工酶进行了比较研究。结果表明,不同科的蜘蛛酯酶同工酶种问差异性大并有着明显的种簇特异性,推测它们的酯酶同工酶酶谱中的区带组受不同的基因位点控制,且各自的基因位点数不等;同种蜘蛛的雌蛛和雄蛛之间也有各自的酯酶同工酶谱型,但差异小,其控制基因位点大体相同。这样,我们从分子的水平上讨论了酯酶同工酶的差异性可以用来作为识别物种的附加指标。     

Fast learning optimized prediction methodology (FLOPRED) for protein secondary structure prediction

Computational methods are rapidly gaining importance in the field of structural biology, mostly due to the explosive progress in genome sequencing projects and the large disparity between the number of sequences and the number of structures. There has been an exponential growth in the number of available protein sequences and a slower growth in the number of structures. There is therefore an urgent need to develop computational methods to predict structures and identify their functions from the sequence. Developing methods that will satisfy these needs both efficiently and accurately is of paramount importance for advances in many biomedical fields, including drug development and discovery of biomarkers. A novel method called fast learning optimized prediction methodology (FLOPRED) is proposed for predicting protein secondary structure, using knowledge-based potentials combined with structure information from the CATH database. A neural network-based extreme learning machine (ELM) and advanced particle swarm optimization (PSO) are used with this data that yield better and faster convergence to produce more accurate results. Protein secondary structures are predicted reliably, more efficiently and more accurately using FLOPRED. These techniques yield superior classification of secondary structure elements, with a training accuracy ranging between 83?% and 87?% over a widerange of hidden neurons and a cross-validated testing accuracy ranging between 81?% and 84?% and a segment overlap (SOV) score of 78?% that are obtained with different sets of proteins. These results are comparable to other recently published studies, but are obtained with greater efficiencies, in terms of time and cost.     

Determination of Melting Sequences in DNA and DNA-Protein Complexes by Difference Spectra

A graphical formula is presented for determining the base ratio of melted DNA. By use of this formula, the composition of sequences which melt in different portions of the melting curves of Clostridium DNA, Escherichia coli DNA, and mouse DNA were determined. As the DNA melts, the per cent of adenine and thymine (AT) in the melted sequences decreases linearly with temperature. The average composition of sequences which melt in a given part of the melting curve is proportional to the base ratio of the DNA. The concentration and average composition of sequences were determined for three parts of the melting curves of the DNA samples, and a frequency distribution curve was constructed. The curve is symmetrical and has a maximum at about 56% AT. The distribution of GC-rich sequences on the E. coli chromosome was estimated by shearing, partially melting, and fractionating the DNA on hydroxylapatite. GC-rich sequences appear to occur every thousand base pairs, and have a maximum length of about 180 base pairs. The graphical formula was applied to the determination of the composition of sequences which melt in different parts of the melting curve of chromatin. Throughout the melting curve, the composition of the melting sequences is about 60% AT, which appears to suggest that relatively long sequences are melting simultaneously. Their melting temperature may be a function of the composition of the protein on different parts of the DNA. The problem of light scattering in DNA-protein and DNA was also investigated. A formula is presented which corrects for light scattering by relating the intensity of the scattered light to the rate of change of absorbance of DNA with wavelength.     

基于混沌游走方法的Rh血型系统中RHD基因的分析

利用基于经典HP模型的蛋白质序列混沌游走方法（chaos game representation,CGR）,给出了RHD基因的蛋白质序列CGR图,可视作蛋白质序列二级结构的一个特征图谱描述．对临床上的血型鉴别有一定的参考价值．另外．还根据由Jeffrey在1990年提出的描绘DNA序列的CGR方法,给出了RHD基因的DNA序列的CGR图．并且根据RHD基因DNA序列的CGR图算出了尺日D基因相应的马尔可夫两步转移概率矩阵,从概率矩阵表可以看出RHD基因对编码氨基酸的三联子的第3个碱基的使用偏好性．     

Molecular evolution of a portion of the mitochondrial 16S ribosomal gene region in scleractinian corals

Relationships among families and suborders of scleractinian corals are poorly understood because of difficulties 1) in making inferences about the evolution of the morphological characters used in coral taxonomy and 2) in interpreting their 240-million-year fossil record. Here we describe patterns of molecular evolution in a segment of the mitochondrial (mt) 16S ribosomal gene from taxa of 14 families of corals and the use of this gene segment in a phylogenetic analysis of relationships within the order. We show that sequences obtained from scleractinians are homologous to other metazoan 16S ribosomal sequences and fall into two distinct clades defined by size of the amplified gene product. Comparisons of sequences from the two clades demonstrate that both sets of sequences are evolving under similar evolutionary constraints: they do not differ in nucleotide composition, numbers of transition and transversion substitutions, spatial patterns of substitutions, or in rates of divergence. The characteristics and patterns observed in these sequences as well as the secondary structures, are similar to those observed in mt 16S ribosomal DNA sequences from other taxa. Phylogenetic analysis of these sequences shows that they are useful for evaluating relationships within the order. The hypothesis generated from this analysis differs from traditional hypotheses for evolutionary relationships among the Scleractinia and suggests that a reevaluation of evolutionary affinities in the order is needed. Received: 4 September 1996 / Accepted: 7 April 1997     

Circular dichroism of films of polynucleotides

The CD spectra of films of the lithium salt of E. coli and calf thymus DNA, and alternating d-AT : AT were measured as a function of relative humidity. Films of the ammonium acetate salt of DNA were also measured. The ammonium films yield the previously reported A-form CD spectra. A possible explanation for the small magnitude of the 260-nm band of the A-form film spectra compared to double-stranded RNA spectra is that the film DNA is in a different conformation than RNA within the A family of conformations. At relative humidities of 92% or lower, a negative nonconservative CD spectrum with negative minima near 270 and 210 nm is observed with the lithium films. The magnitude of the minima varies from film to film. In films of DNA the magnitude ranges from a delta epsilon of ?5 to ?35; d-AT : AT films show magnitudes to ?300. CD spectra of this type are designated Ψ spectra. Similar spectra have been reported from reconstituted complexes of DNA and polylysine or f-1 histone. If the origins of the film and protein–DNA complex spectra are similar, the complex spectra are not the result of specific secondary structural changes induced in the DNA by the protein fraction. Theoretical analysis suggests that Ψ spectra are not the result of changes in the secondary or tertiary structure of DNA. Instead, the previously proposed explanation based on liquid crystals is favored. The DNA could form asymmetric structures with long-range periodicity. It is likely that the observed CD spectra of f-1 complexes are artifacts of DNA aggregation. The possibility that some other previously published spectra of protein–DNA complexes also reflect artifacts is suggested.     

Structural conservation and variation in the D-loop-containing region of vertebrate mitochondrial DNA

The nucleotide sequences of the D-loop-containing regions of three rat mitochondrial DNAs (mtDNAs), two from the species Rattus norvegicus and one from R. rattus, were determined. Comparisons made among these sequences and with the mouse sequence showed that, on the basis of both base composition and frequency of nucleotide alterations, three domains could be defined within the D-loop-containing region: a central conserved segment, poor in L-strand adenine, flanked by two divergent, adenine-rich regions. Deletions and insertions were found to occur at an unexpectedly high frequency in these sequences and the conserved sequence block called CSB-1 was found not to be intact in the R. rattus sequence. Although in comparisons of more distantly related mtDNAs the D-loop region is the most divergent on the molecule, it does not diverge more than typical protein genes between R. norvegicus and R. rattus, and its central conserved domain appears to be one of the molecule's most conserved regions. The most variable domain borders the tRNAPhe gene and contains the L and H-strand promoters and the 5' terminus for H-strand DNA synthesis. Within this region we have found sequences in all the mtDNAs we have examined, including those of human, two artiodactyls and Xenopus, that are capable of folding into cloverleaf structures. In the other divergent domain of the same mtDNAs, we find sequences capable of assuming similar secondary structural configurations at or near the sites for the termination of D-loop DNA synthesis. The evolutionary preservation of the potential to form such structures despite the high primary-structural divergence of the regions they occur in, suggests the structures are of principal importance for some processes occurring in the D-loop-containing region.     

Identification of related proteins with weak sequence identity using secondary structure information

Molecular modeling of proteins is confronted with the problem of finding homologous proteins, especially when few identities remain after the process of molecular evolution. Using even the most recent methods based on sequence identity detection, structural relationships are still difficult to establish with high reliability. As protein structures are more conserved than sequences, we investigated the possibility of using protein secondary structure comparison (observed or predicted structures) to discriminate between related and unrelated proteins sequences in the range of 10%-30% sequence identity. Pairwise comparison of secondary structures have been measured using the structural overlap (Sov) parameter. In this article, we show that if the secondary structures likeness is >50%, most of the pairs are structurally related. Taking into account the secondary structures of proteins that have been detected by BLAST, FASTA, or SSEARCH in the noisy region (with high E: value), we show that distantly related protein sequences (even with <20% identity) can be still identified. This strategy can be used to identify three-dimensional templates in homology modeling by finding unexpected related proteins and to select proteins for experimental investigation in a structural genomic approach, as well as for genome annotation.     

猪胸膜肺炎放线杆菌转铁结合蛋白基因（rbp8）分析及建模

为了深入研究猪胸膜肺炎放线杆菌（Actinobacillus pleuropneumonie,App）转铁结合蛋白基因（Transferrin BindingProtein8,脚）的生物学特性,采用生物信息学方法,对GenBank中的5株App的TbpB的核酸及其氨基酸序列进行比对,选取其中的中国湖北分离株（JL03）对其分子结构、理化性质及功能域、蛋白质二级和三级结构等重要参数进行了预测和分析,并在三级结构的基础上进行了同源建模。结果表明,不同APP菌株之间核酸序列相似性较大,而氨基酸序列存在较大差异,二级结构以延伸链和随机卷曲为主要构件,其空间结构与脑膜炎双球菌GNAl870蛋白相似性较高,以此为模板成功构建了三维结构分子模型,为TbpB基因功能的深入研究提供了线索和参考依据。     

Estimating the significance of sequence order in protein secondary structure and prediction

MOTIVATION: How critical is the sequence order information in predicting protein secondary structure segments? We tried to get a rough insight on it from a theoretical approach using both a prediction algorithm and structural fragments from Protein Databank (PDB). RESULTS: Using reverse protein sequences and PDB structural fragments, we theoretically estimated the significance of the order for protein secondary structure and prediction. On average: (1) 79% of protein sequence segments resulted in the same prediction in both normal and reverse directions, which indicated a relatively high conservation of secondary structure propensity in the reverse direction; (2) the reversed sequence prediction alone performed less accurately than the normal forward sequence prediction, but comparably high (2% difference); (3) the commonly predicted regions showed a slightly higher prediction accuracy (4%) than the normal sequences prediction; and (4) structural fragments which have counterparts in reverse direction in the same protein showed a comparable degree of secondary structure conservation (73% identity with reversed structures on average for pentamers). CONTACT: jong@biosophy.org; dietmann@ebi.ac.uk; heger@ebi.ac.uk; holm@ebi.ac.uk     

草鱼(Ctenopharyngodon idellus)线粒体DNA控制区及其两侧tRNA基因的克隆与结构分析

动物线粒体ＤＮＡ控制区是线粒体基因组复制与基因表达的最主要的调控区．采用杂交和测序的方法对草鱼线粒体ＤＮＡ控制区进行定位、克隆并测定了控制区及其旁侧的ｔＲＮＡＰｈｅ、ｒＲＮＡＰｒｏ和ｒＲＮＡＴｈｒ三个基因的序列,与多种脊椎动物的相应序列进行了比较,并进行了结构分析．草鱼线粒体控制区全长９２７ｂｐ,含有与酵母和爪蟾线粒体启动子相似的序列,其ＣＳＢⅠ、ＣＳＢⅡ和ＣＳＢⅢ序列与其他几种动物的ＣＳＢ比较相当保守,ＴＡＳ与其回文基序可形成稳定的茎环结构,成为Ｈ－链复制的终止信号．草鱼线粒体ｔＲＮＡＰｈｅ、ｔＲＮＡＰｒｏ和ｔＲＮＡＴｈｒ可折叠成三叶草形二级结构,其基因具有许多不同于细胞质ｔＲＮＡ基因的结构特点,可能反映了线粒体ｔＲＮＡ与线粒体核糖体具有不寻常的作用方式