象甲科线粒体基因组特征及系统发育分析

昆虫线粒体基因组广泛应用于系统发育关系的重新建立、分子进化、谱系地理学及物种诊断等领域。为揭示象甲科昆虫线粒体全基因组序列的主要结构特征,探究其系统发育相关信息,为进化遗传学研究和分子标记选取等提供参考依据,本研究利用比较基因组学和生物信息学方法,对NCBI上已公布的35种象甲科物种线粒体全基因组序列进行了分析。结果显示:(1)象甲科tRNA基因存在排序及数目异常情况,不同物种中蛋白质编码基因和2种rRNAs排列相同,线粒体全基因组具有明显AT偏向;(2)COX1、ATP6、ND5、ND4、ND4L和ND1基因除标准三联密码子外,还存在特殊的起始密码子AAT、TTG和GTG;(3)13种蛋白质编码基因的进化速率顺序为COX3ATP8ND2ND5ND1ND4ND6ND4LND3ATP6CytbCOX1COX2;(4)13个蛋白编码基因和rRNAs基因中,ND5、rrnL、ND4和ND2基因变异位点数较高,可作为备选的分子标记;(5)各亚科的系统发育关系可能为(((小蠹亚科Scolytinae+长小蠹亚科Platypodinae)+(隐喙象亚科Cryptorhynchinae+魔喙象亚科Molytinae+象虫亚科Curculioninae)+((孢喙象亚科Cyclominae+粗喙象亚科Entiminae)+(隐颏象亚科Dryophthorinae+长小蠹亚科))),为象甲科的系统发育分析有提供参考。     

使用伪氨基酸组成和模糊支持向量机预测蛋白质结构类

蛋白质结构类预测是生物信息和蛋白质科学中重要的研究领域.基于Chou提出的伪氨基酸离散模型框架,从蛋白质序列出发,设计一种新的伪氨基酸组成方法表示蛋白质序列样本.抽取氨基酸组合(10-D)在序列中出现的频率和疏水氨基酸模式(6-D)表示蛋白质序列的附加特征,用和传统的氨基酸组成(20-D)一起构成的36维的伪氨基酸组成向量来表示蛋白质序列的特征.使用遗传算法来优化附加特征的权重系数.伪氨基酸组成向量作为输入数据,模糊支持向量机作为预测工具.使用三个常用的标准数据集来验证算法的性能.Jack-knife检验结果说明本方法具有较高的准确率,有望成为潜在的预测蛋白质功能的工具.     

蛋白质序列的一种新的三维图形表示及其应用

基于氨基酸的五字母模型,给出蛋白质序列的一种新的三维图形表示,然后构造一个12维向量来刻画蛋白质序列,这个向量的分量是与12个图形相对应的D／D矩阵的正规化的ALE-指标。最后基于s结构蛋白对冠状病毒进行系统发生分析来阐明该方法的有用性。     

不同特征描述下H1N1病毒血凝素蛋白序列的比较分析

在传统表征蛋白质序列的40维特征向量的基础上,依据氨基酸的种类与理化性质,将蛋白质序列40维向量分解为20维、4维和16维3种子特征向量描述。结合33条H1N1流感病毒血凝素(hemagglutinin,HA)蛋白质序列和统计学相关性分析理论,进行了蛋白质序列两两之间及每条病毒蛋白质序列对应的不同子特征向量之间的相关性分析,发现病毒蛋白质序列之间存在高度相关性,且每条病毒蛋白质序列对应的20维子特征向量与其他两种子特征向量之间均不显著相关,而4维与16维子特征向量之间显著相关。进一步依据不同的特征向量对33条HA蛋白质序列进行分类,研究发现依据40维特征向量与16维特征向量进行的分类结果高度一致。因此,在不影响表征病毒序列特性的前提下,对于已有的表征蛋白质序列的40维特征向量,可以用16维的特征向量进行代替,以减少计算复杂度。     

基于序列特征网络的蛋白质结构类型研究

利用复杂网络的方法来探索序列特征因素对蛋白质结构的影响。由于蛋白质的序列对结构具有重要且复杂的影响,因此将蛋白质的结构以及序列特征之间的关系模拟成一个复杂系统,通过利用互相关系数、标准化互信息和传递熵等方法来建立以序列特征为节点的加权网络,进而利用网络中心性的方法来分析不同蛋白质结构类型对应加权网络的中心性分布的差异,探索不同结构类型蛋白质的序列特征差异。发现不同的蛋白质结构类型对应的序列特征网络既有共性又有差异,文章将针对每一种结构类型的网络中心性分布,以及不同结构类型之间的共性与差异进行详细地讨论。研究结果对蛋白质序列与结构之间关系的研究,特别是结构分类研究具有重要的意义。     

蛋白质二级结构预测: 基于词条的最大熵马尔科夫方法

提出了一种新的蛋白质二级结构预测方法. 该方法从氨基酸序列中提取出和自然语言中的“词”类似的与物种相关的蛋白质二级结构词条, 这些词条形成了蛋白质二级结构词典, 该词典描述了氨基酸序列和蛋白质二级结构之间的关系. 预测蛋白质二级结构的过程和自然语言中的分词和词性标注一体化的过程类似. 该方法把词条序列看成是马尔科夫链, 通过Viterbi算法搜索每个词条被标注为某种二级结构类型的最大概率, 其中使用词网格描述分词的结果, 使用最大熵马尔科夫模型计算词条的二级结构概率. 蛋白质二级结构预测的结果是最优的分词所对应的二级结构类型. 在4个物种的蛋白质序列上对这种方法进行测试, 并和PHD方法进行比较. 试验结果显示, 这种方法的Q3准确率比PHD方法高3.9%, SOV准确率比PHD方法高4.6%. 结合BLAST搜索的局部相似的序列可以进一步提高预测的准确率. 在50个CASP5目标蛋白质序列上进行测试的结果是: Q3准确率为78.9%, SOV准确率为77.1%. 基于这种方法建立了一个蛋白质二级结构预测的服务器, 可以通过http://www.insun.hit.edu.cn:81/demos/biology/index.html来访问.     

日本条螽线粒体基因组序列的测定与分析

日本条螽完整的线粒体基因组序列长16 281 bp,包括13个蛋白质编码基因、22个tRNA基因、2个r RNA基因和1个D-loop区,其基因次序和方向与祖先序列相同。该线粒体基因组排列紧凑,但在ND2和tRNA~(Trp)之间有一段长为650 bp的基因间隔区。为研究螽斯科的系统发育关系,本研究选取日本条螽及其它17个螽斯科物种线粒体基因组的蛋白质编码基因和r RNA基因序列构建贝叶斯系统发生树。     

3种棱型大麦中脱水素6(DHN6)的性质、突变位点与蛋白质进化分析

脱水素(dehydrins,DHNs)是高等植物胚胎发育晚期产生的一类特异多肽,能够逆转细胞脱水所造成的伤害.为明确脱水素性质与功能的关系,本研究从3种棱型大麦分离到Dhn6基因,生物信息学分析其编码的蛋白质序列长度分别为523个(六棱)、502个(四棱)和486个(二棱)氨基酸残基;氨基酸突变位点分析发现该基因具有整体保守性和突变位点偏倚性.蛋白质性质与二级结构分析表明,DHN6是高度亲水的碱性蛋白质,线性结构和自由卷曲为二级结构的主要组分,K-片段参与α-螺旋结构的形成,由此推测DHN6的兼性α-螺旋结构域在稳定膜结构的水合保护体系过程中发挥着重要的生理功能.构建的21个物种DHN6蛋白质系统发生树,发现禾本科植物遗传距离较近,结合核苷酸序列的物种特异性,本研究认为Dhn6基因可作为物种鉴定依据.     

基于物化性质对嗜热蛋白的预测

嗜热蛋白在高温下能保持稳定性和活性,是研究蛋白质热稳定性的理想模型,开发一个蛋白质热稳定性识别的方法将对蛋白质工程和蛋白质的设计很有帮助。目前的研究中,氨基酸的组成及其物化性质一直被认为和蛋白质的热稳定性相关。本研究筛选出可靠的数据集,包括915个嗜热蛋白和793个非嗜热蛋白。利用蛋白质氨基酸的物化性质和氨基酸的组成表征嗜热蛋白,将二肽氨基酸组成整合到9组氨基酸物化性质中使蛋白序列公式化。支持向量机5折叠交叉验证表明:当gap=0时,290个特征产生的精度最高,为92.74%。因此说明对于分析蛋白质的热稳定性,所建立的预测模型将是一个很有效的工具。     

牦牛Eppin基因CDS的克隆及其在附睾不同区段中的表达研究

为了研究牦牛附睾组织中精子成熟的相关机理,并为探讨高原动物的生殖机制提供基本数据。本研究运用基因克隆技术对牦牛附睾Eppin基因CDS全长序列进行克隆,采用生物信息学方法进行分析,Eppin基因和编码序列特征进行了预测和分析。结果表明,牦牛Eppin基因的CDS含有一个405 bp长度的片段,由134个氨基酸编码;牦牛Eppin基因对应的蛋白分子量和理论等电点分别为15.09 ku和8.67 ku,其对应的氨基酸没有跨膜结构,归于近水性蛋白;25个α-螺旋、27个延伸链、2个β-折叠及80个无规则卷曲构成其蛋白质二级结构;牦牛Eppin基因编码氨基酸序列与黄牛、藏羚羊、绵羊等物种间同源性较高,系统进化情况与其亲缘关系远近一致。本研究应用实时荧光定量PCR技术分析Eppin基因在附睾组织3个不同区段(头部,颈部和尾部)中的表达情况,荧光定量PCR结果显示,Eppin基因在牦牛附睾组织3个不同区段中均有不同程度的表达,在附睾头部中表达最高,颈部和尾部表达较低。本研究将为牦牛附睾精子成熟的机制和Eppin基因在牦牛附睾上皮细胞中的功能提供一定的基础数据。     

Cloning of the mitochondrial genome of Rana catesbeiana and the nucleotide sequences of the ND2 and five tRNA genes

The entire mitochondrial genome of Rana catesbeiana was cloned into a plasmid vector pBR322 at the unique BamHI site and the nucleotide sequences of the ND2 gene and of its flanking genes were determined. The ND2 gene was encoded by 1,033 base pairs and, as deduced from the nucleotide sequence, the ND2 product consisted of 344 amino acids with a molecular weight of 37,561. This gene was flanked on the 5' side by the tRNA genes for isoleucine, glutamine, and methionine and on the 3' side by those for tryptophan and alanine. These genes were the same in their organization as those found in the mammalian and Xenopus laevis mitochondrial genomes. A comparison of the putative amino acid sequences of the ND2 proteins of different animal species revealed that six regions in the sequence were well conserved during evolution, suggesting that some of these conserved sequences are crucial for biological activity of the ND2 protein. The nucleotide sequence homologies between the five tRNA genes of R. catesbeiana and their counterparts of mammals and X. laevis were in the range of 55 to 85%, depending on the tRNA and animal species.     

Neutral and Non-Neutral Evolution of Drosophila Mitochondrial DNA

To test hypotheses of neutral evolution of mitochondrial DNA (mtDNA), nucleotide sequences were determined for 1515 base pairs of the NADH dehydrogenase subunit 5 (ND5) gene in the mitochondrial DNA of 29 lines of Drosophila melanogaster and 9 lines of its sibling species Drosophila simulans. In contrast to the patterns for nuclear genes, where D. melanogaster generally exhibits much less nucleotide polymorphism, the number of segregating sites was slightly higher in a global sample of nine ND5 sequences in D. melanogaster (s = 8) than in the nine lines of D. simulans (s = 6). When compared to variation at nuclear loci, the mtDNA variation in D. melanogaster does not depart from neutral expectations. The ND5 sequences in D. simulans, however, show fewer than half the number of variable sites expected under neutrality when compared to sequences from the period locus. While this reduction in variation is not significant at the 5% level, HKA tests with published restriction data for mtDNA in D. simulans do show a significant reduction of variation suggesting a selective sweep of variation in the mtDNA in this species. Tests of neutral evolution based on the ratios of synonymous and replacement polymorphism and divergence are generally consistent with neutral expectations, although a significant excess of amino acid polymorphism within both species is localized in one region of the protein. The rate of mtDNA evolution has been faster in D. melanogaster than in D. simulans and the population structure of mtDNA is distinct in these species. The data reveal how different rates of mtDNA evolution between species and different histories of neutral and adaptive evolution within species can compromise historical inferences in population and evolutionary biology.     

Characterization of protein primary sequences based on partial ordering

In this paper, we present a new approach to characterize protein sequences. Based on orderings of the 20 natural amino acids which reflect some of their physico-chemical properties, we construct an augmented Hasse matrix for each protein sequence. Furthermore, the normalized leading eigenvalues of these matrices are computed and considered as invariants for the protein sequences. Finally, we make a comparison for the similarity/diversity of nine different protein sequences.     

一种快速非比对的蛋白质序列相似性与进化分析方法

本文提出了一种新的快速非比对的蛋白质序列相似性与进化分析方法。在刻画蛋白质序列特征时,首先将氨基酸的10种理化性质通过主成分分析浓缩为6个主成分,并且将每条蛋白质序列里的氨基酸数目作为权重对主成分得分值进行加权平均,然后再融合氨基酸的位置信息构成一个26维的蛋白质序列特征向量,最后利用欧式距离度量蛋白质序列间的相似性及进化关系。通过对3个蛋白质序列数据集的测试表明,本文提出的方法能将每条蛋白质序列准确聚类,并且简便快捷,说明了该方法的有效性。     

Molecular differentiation and phylogeny of entomopathogenic nematodes (rhabditida: heterorhabditidae) based on ND4 gene sequences of Mitochondrial DNA.

We determined partial ND4 gene sequences of mitochondrial DNA from 15 heterorhabditid nematode isolates, representing 5 species collected from different regions of the world, by using polymerase chain reaction (PCR) and direct-sequencing of PCR products. Aligned nucleotide as well as amino acid sequences were used to differentiate nematode species by comparing sequence divergence and to infer phylogeny of the nematodes by using maximum parsimony and likelihood methods. Robustness of our phylogenetic trees was checked by bootstrap tests. The 15 nematode isolates can be divided into 7 haplotypes based on DNA sequences. On a larger scale, the sequence divergence revealed 4 distinct groups corresponding to 4 described species. No sequence divergence was detected from 5 isolates of Heterorhabditis bacteriophora or between Heterorhabditis marelatus to Heterorhabditis hepialius. Our sequence data yielded phylogenetic trees with identical topologies when different tree-building methods were used. Most relationships were also confirmed by using amino acid sequences in maximum parsimony analysis. Our molecular phylogeny of Heterorhabditis species support an existing taxonomy that is based largely on morphology and the sequence divergence of the ND4 gene permits species identification.     

cDNA sequence analysis of ribosomal protein S13 gene in Plutella xylostella （Lepidoptera： Plutellidae）

Ribosomal protein S 13 gene has been cloned and analyzed in many organisms,but there are few documents relating to insects. In this communication, the full-length cDNA sequence of ribosomal protein S 13 gene in the diamondback moth, Plutella xylostella(Lepidoptera: Plutellidae), was determined by using PCR amplification technique. The features of the ribosomal protein S 13 gene sequence were analyzed and the deduced amino acids sequence was compared with those from other insects. The results of multi-alignment of the amino acid sequences between the diamondback moth and other insect species revealed that this gene sequence is highly conserved in insects. Based on maximum likelihood method, a phylogenetic tree was constructed from 10 different species using PHYLIP software. It showed that nematode is one separate lineage and the five insect speciesbe long to another lineage, whereas those species higher than insects form the third one. The pattern of this phylogenetic tree evidently represented the evolution of different species.     

Transforming activity of a 16-amino-acid segment of the bovine papillomavirus E5 protein linked to random sequences of hydrophobic amino acids.

The 44-amino-acid E5 protein of bovine papillomavirus type 1 is the smallest transforming protein yet described. Previous results from our laboratory indicate that a hydrophobic core and specific carboxyl-terminal amino acids are required for the E5 protein to exert its transforming function. In this study, additional substitution mutations were generated in the E5 gene to determine the minimal amino acid sequence requirements for focus formation in mouse C127 cells. In most cases examined, substitution of the hydrophobic middle third of the E5 protein with unrelated hydrophobic sequences severely inhibited transforming activity. However, we have identified one hydrophobic amino acid sequence apparently unrelated to the wild-type one that can replace the middle third of the wild-type E5 protein without affecting the ability of the protein to stably transform cells or interact with cell membranes. Furthermore, a mutant E5 protein in which only the carboxyl-terminal 16 amino acids of the protein have been derived from E5 sequences retains transforming activity. Since several residues in the carboxyl-terminal portion of the E5 protein can be freely substituted with different amino acids (B. H. Horwitz, A. L. Burkhardt, R. Schlegel, and D. DiMaio, Mol. Cell. Biol. 8:4071-4078, 1988), the results reported here imply that much of the specific information necessary for cell transformation can be supplied by a subset of the carboxyl-terminal 16 amino acids of this protein.     

细胞色素c蛋白序列分析与结构比较

首先以马心细胞色素c(Horse Cytc)蛋白的氨基酸序列为查询序列,利用生物信息学方法进行相似性搜索,获得了一系列细胞色素c(Cytc)蛋白的氨基酸序列,然后对Cytc蛋白进行了多重对齐分析、进化分析和三维结构比较分析。分析结果表明:Cytc中某些特定部位的氨基酸残基高度保守;相近物种来源的Cytc具有较近的亲缘关系,而来源于同一物种不同部位的Cytc却具有较远的亲缘关系;来源于不同物种的Cytc,即使具有较远的亲缘关系,却具有极其相似的三维空间结构。这些研究结果将为基于Cytc进行蛋白分子设计与构建提供指导意义。