云南HPV-16型E6、E7基因的突变分析

人乳头瘤病毒(Human papillomavirus,HPV)16型(HPV-16)是引起宫颈癌的一种主要高危型病毒,其2个致癌基因E6和E7的核酸序列变异可能会影响其对宿主细胞的致癌性,已有研究表明其序列突变呈现地域差异性。因此,研究不同地域HPV-16这2个基因的变化情况是宫颈癌流行病学调研的主要内容,也可为研究E6和E7的致癌性积累数据。研究以NCBI登录号为NC_001526.2的HPV-16型病毒的序列为参照,采用Neighbor-joining方法对云南地区74例HPV-16样本的E6、E7的DNA序列构建进化树,结果显示:只有亚洲和欧洲变异亚型,而没有发现非洲1、非洲2、亚-美洲和北美洲这4种变异亚型。DNA序列分析显示:E6的碱基突变以T178G(D25E,59.46%)和T350G(L83V,8.11%)为主,E7的碱基突变主要以A647G(N29S,59.46%)和T846C(同义突变,60.81%)为主。发现E6的新突变有A95G(同义突变,1.35%)和A135G(K11R,1.35%);E7的新突变有C625T(L22F,1.35%)、C627T(同义突变,12.16%)、G689A(G43E,1.35%)、T748G(S63A,1.35%)。此外还发现有一个共突变现象:T178G(D25E,59.46%)-A647G(N29S,59.46%)-T843C(同义突变,21.62%)-T846C(同义突变,60.81%)。     

不同地区HPV16E7基因的克隆及序列差异分析

采用PCR扩增、pGEM-T载体克隆和核苷酸序列分析的方法对一 例武汉地区及两例五峰县高发区宫颈癌患者体内HPV16型的E7基因编码区进行序列分析并与 野生型(德国标准株)及已发表的HPV16湖北株(HPVHB)进行了比较.结果发现武汉地区HPV16 型E7基因仅第54位出现一个同义突变,而高发区HPV16型E7基因存在差异,第77位氨基酸由 精氨酸(Arg)变为半胱氨酸(Cys),第96位由谷氨酰氨酸(Gln)变为精氨酸(Arg),E7蛋白的二级 结构及亲、疏水性也相应改变,与野生型有较大差异.     

新疆维吾尔族妇女宫颈癌组织中乳头瘤病毒16型E6基因的克隆和序列分析

为了分析新疆南部地区维吾尔族妇女宫颈癌组织中HPV16型E6基因结构特点，从中国新疆南部地区维吾尔族妇女宫颈癌活检组织标本中提取DNA，以宫颈癌活检组织标本DNA为模板进行PCR扩增，获得HPV16 E6基因，将其克隆到pUCm-T载体上，并对其进行基因全序列分析.PCR检测结果显示宫颈癌组织中HPV16 E6阳性率为82.35%（14/17）；测序结果显示，新疆株HPV16 E6基因全长456 bp，大小与德国标准株一致.E6基因的第247位碱基发生T→G突变，并由此引起所编码的氨基酸亦发生改变.上述结果表明，中国新疆南部地区维吾尔族妇女宫颈癌患者组织中HPV16 E6的基因结构与德国标准株HPV16 E6基因之间存在差异.     

尖锐湿疣病变的人乳头瘤病毒6型L1序列多态性分析

采用PCR方法从协和医院尖锐湿疣患者皮损中检测人乳头瘤病毒(HPV),并通过限制性片断长度多态性分析分型发现,HPV6型为主要感染型别,其次是HPV11型.根据临床特征选择主要致病型HPV6 8个分离株,扩增L1晚期基因,构建重组测序质粒,双脱氧法测序并分析L1区的核苷酸和氨基酸序列变异状况.结果表明,HPV6L1基因序列发生碱基替换的区域主要有四个区,包括SR1(5911～6104),SR2(6217～6273),SR3(6540～6661)和SR4(7062～7250).少数的碱基替换导致错义突变,推导的蛋白质一级结构中有0～3个氨基酸发生变异,但抗原性强弱与原型HPV6基本一致.     

不同地区HPV16 E7基因的克隆及序列差异分析

采用PCR扩增、pGEM T载体克隆和核苷酸序列分析的方法对一例武汉地区及两例五峰县高发区宫颈癌患者体内HPV16型的E7基因编码区进行序列分析并与野生型 (德国标准株 )及已发表的HPV16湖北株 (HPVHB)进行了比较。结果发现武汉地区HPV16型E7基因仅第 5 4位出现一个同义突变 ,而高发区HPV16型E7基因存在差异 ,第 77位氨基酸由精氨酸 (Arg)变为半胱氨酸 (Cys) ,第 96位由谷氨酰氨酸 (Gln)变为精氨酸 (Arg) ,E7蛋白的二级结构及亲、疏水性也相应改变 ,与野生型有较大差异     

人乳头瘤病毒16 E6蛋白通过Wnt/β-catenin通路促进食管癌Eca109细胞增殖、迁移的实验研究

人乳头瘤病毒16(Human papillomavirus 16,HPV16)感染与食管癌的发生密切相关,HPV16编码的HPV16E6蛋白是主要的致癌蛋白,已经在宫颈癌细胞中被证实能够促进癌细胞的增殖、迁移,同时也可激活Wnt/β-catenin通路。但HPV16 E6蛋白对食管癌细胞增殖、迁移的调控作用及机制尚未阐明。为研究HPV16 E6蛋白对食管癌Eca109细胞增殖、迁移及细胞中Wnt/β-catenin通路的调节作用,本研究培养食管癌Eca109细胞并分组,空白对照组用不含药物及质粒的DMEM处理,空白pcDNA3.1组转染空白的pcDNA3.1质粒,HPV16 E6组转染HPV16 E6基因的pcDNA3.1质粒,HPV16 E6+XAV组转染HPV16 E6基因的pcDNA3.1质粒、同时加入β-catenin抑制剂XAV939;并检测细胞的增殖、迁移活力及细胞中增殖基因、迁移基因、Wnt/β-catenin通路基因的表达。结果显示,转染24h后,HPV16 E6组细胞的增殖活力、迁移活力以及细胞中c-myc、cyclinD1、bcl-2、HMGA-2、N-cadherin、Wnt2、β-catenin的蛋白表达量均明显高于空白对照组、空白pcDNA3.1组;HPV16 E6+XAV组细胞的增殖活力、迁移活力以及细胞中c-myc、cyclinD1、bcl-2、HMGA-2、N-cadherin、β-catenin的蛋白表达量均明显低于HPV16 E6组,Wnt2的蛋白表达量与HPV16 E6组比较无明显差异。本研究揭示,HPV16 E6蛋白能够促进食管癌Eca109细胞的增殖、迁移且该作用与激活Wnt/β-catenin通路有关。     

板角山羊mtDNA Cyt b基因变异特点及系统进化研究

测定了板角山羊品种13个个体的细胞色素b基因全序列(1140 bp),比较分析了群体中细胞色素b基因的碱基组成和序列间碱基的变异情况,结果显示:在该品种(群体)中细胞色素b基因序列中6个变异位点上观察到11次T-C间和2次A-G间的碱基转换,除了有2次T-C间碱基转换发生在密码子第2位点为非同义突变以外,其余的11次碱基转换发生在密码子第3位点,均为同义突变;有1次T-G间碱基颠换发生在密码子第2位点,为非同义突变;观察到5种单倍型,单倍型多样度为0.8077.并以绵羊为外群,与山羊属其他种的同源区序列构建系统发生树.结果显示, 在系统地位上板角山羊与胃石山羊有较近的亲缘关系.     

辽宁省地区人乳头瘤病毒33型基因多态性分析

研究人乳头瘤病毒（Human papillomavirus,HPV）33型E6和E7基因在辽宁省地区的基因多态性分布情况,增加对HPV33基因多态性地理分布情况的了解,为HPV33的诊断,靶向药和疫苗的研发提供理论依据。使用巢式PCR扩增方法对从已确诊HPV33阳性患者宫颈细胞提取的DNA样本进行E6和E7基因的扩增并测序,随后应用MegaX软件将所得测序结果与参考序列比对确定突变位点并进行系统发育树的构建。应用PAML4.9软件中的Codeml程序找出突变中的正向选择位点。分别应用ABC Pred server,ProPred-I server和ProPred server寻找理想的B细胞免疫表位,人类白细胞抗原（Human leukocyte antigen,HLA）I类结合表位和人类白细胞抗原（HLA）II类结合表位。共得到HPV33E6和E7序列各136个。与HPV33参考序列（M12732.1）进行比对后,E6基因中观测到17个单核苷酸突变,同义突变7个,非同义突变10个。E7基因中观测到9个单核苷酸突变,其中同义突变3个,非同义突变6个。系统发育分析显示HPV33E6和E7...     

中国湖北地区HPV16型E6和E7基因变异和进化发生学研究

研究中国湖北地区宫颈癌患者的人乳头瘤病毒16型E6和E7的变异以及HPV16变异体的分布。从宫颈癌患者手术切除标本提取组织DNA,用HPV16 E6和E7特异性引物进行PCR扩增,对扩增的部分E6和E7产物片段进行测序分析。在80例宫颈癌组织DNA中有41例发生E6基因178位核苷酸的突变,突变频率58.75%,相应核苷酸改变为Asp-Glu,E7 647在31例测序样品中有22例发生核苷酸序列A到G改变,使29位氨基酸由Asn变为Ser,突变频率70.97%,结果显示在E6和E7基因的178位和647位核苷酸存在高频率的碱基变异。对E6和E7基因的进化树分析表明,中国湖北地区流行的HPV16病毒株主要为亚洲型变异体(As),其次为欧洲型(E),没有发现非洲-1型(Af-1),非洲-2型(Af-2)和亚洲美洲型(AA)HPV16变异体,中国湖北地区流行的As变异体是否有更高的致宫颈癌的风险还有待于进一步对不同阶段CIN和正常宫颈上皮样品的E6和E7基因进行序列分析和对变异体蛋白进行功能研究。     

蕨类植物psbD基因的适应性进化和共进化分析

D2蛋白是植物光系统Ⅱ复合体（PSⅡ）核心蛋白之一,由叶绿体psbD基因编码。为了深入理解核心薄囊蕨类植物在阴生环境下的“辐射”式演化,我们对12种蕨类植物的psbD基因进行了克隆和测序,然后联合已公布的其他8种蕨类植物的psbD序列,基于ω值（非同义替换率ds和同义替换率ds的比值）探讨了该基因经受的选择压力。发现D2蛋白在大多数分支和位点受到强烈的负选择,但是树蕨类分支的psbD进化速率低且ω值较高。借助多种模型进行的共进化分析显示,树蕨类D2蛋白的168R、245H和272M两两组成具有共进化关系的氨基酸位点对。     

Rates of Conservative and Radical Nonsynonymous Nucleotide Substitutions in Mammalian Nuclear Genes

To understand the process and mechanism of protein evolution, it is important to know what types of amino acid substitutions are more likely to be under selection and what types are mostly neutral. An amino acid substitution can be classified as either conservative or radical, depending on whether it involves a change in a certain physicochemical property of the amino acid. Assuming Kimura's two-parameter model of nucleotide substitution, I present a method for computing the numbers of conservative and radical nonsynonymous (amino acid altering) nucleotide substitutions per site and estimate these rates for 47 nuclear genes from mammals. The results are as follows. (1) The average radical/conservative rate ratio is 0.81 for charge changes, 0.85 for polarity changes, and 0.49 when both polarity and volume changes are considered. (2) The radical/conservative rate ratio is positively correlated with the nonsynonymous/synonymous rate ratio for charge changes or when both polarity and volume changes are considered. (3) Both the conservative/synonymous rate ratio and the radical/synonymous rate ratio are lower in the rodent lineage than in the primate or artiodactyl lineage, suggesting more intense purifying selection in the rodent lineage, for both conservative and radical nonsynonymous substitutions. (4) Neglecting transition/transversion bias would cause an underestimation of both radical and conservative rates and the ratio thereof. (5) Transversions induce more dramatic genetic alternations than transitions in that transversions produce more amino acid altering changes and among which, more radical changes. Received: 6 April 1999 / Accepted: 16 August 1999     

Diversifying selection in human papillomavirus type 16 lineages based on complete genome analyses

Human papillomavirus type 16 (HPV16) is the primary etiological agent of cervical cancer, the second most common cancer in women worldwide. Complete genomes of 12 isolates representing the major lineages of HPV16 were cloned and sequenced from cervicovaginal cells. The sequence variations within the open reading frames (ORFs) and noncoding regions were identified and compared with the HPV16R reference sequence. This whole-genome approach gives us unprecedented precision in detailing sequence-level changes that are under selection on a whole-viral-genome scale. Of 7,908 base pair nucleotide positions, 313 (4.0%) were variable. Within the 2,452 amino acids (aa) comprising 8 ORFs, 243 (9.9%) amino acid positions were variable. In order to investigate the molecular evolution of HPV16 variants, maximum likelihood models of codon substitution were used to identify lineages and amino acid sites under selective pressure. Five codon sites in the E5 (aa 48, 65) and E6 (aa 10, 14, 83) ORFs were demonstrated to be under diversifying selective pressure. The E5 ORF had the overall highest nonsynonymous/synonymous substitution rate (omega) ratio (M3 = 0.7965). The E2 gene had the next-highest omega ratio (M3 = 0.5611); however, no specific codons were under positive selection. These data indicate that the E6 and E5 ORFs are evolving under positive Darwinian selection and have done so in a relatively short time period. Whether response to selective pressure upon the E5 and E6 ORFs contributes to the biological success of HPV16, its specific biological niche, and/or its oncogenic potential remains to be established.     

Analysis of substitution events in HIV-1 vif gene of the Korean clade

Nucleotide and amino acid substitution pattern in vif gene of the Korean clade of HIV-1 isolated from Koreans were analyzed using consensus sequences. At nucleotide level, transition/transversion substitution ratio was 1.88, and nonsynonymous/synonymous substitution ratio was 2.67, suggesting a divergent evolution in the Korean clade. At amino acid level, there were 17 substitutions and G-->E substitution at position 37 may be responsible for change in predicted secondary structure.     

Molecular evolution of hemagglutinin genes of H1N1 swine and human influenza A viruses

Summary The hemagglutinin (HA) genes of influenza type A (H1N1) viruses isolated from swine were cloned into plasmid vectors and their nucleotide sequences were determined. A phylogenetic tree for the HA genes of swine and human influenza viruses was constructed by the neighbor-joining method. It showed that the divergence between swine and human HA genes might have occurred around 1905. The estimated rates of synonymous (silent) substitutions for swine and human influenza viruses were almost the same. For both viruses, the rate of synonymous substitution was much higher than that of nonsynonymous (amino acid altering) substitution. It is the case even for only the antigenic sites of the HA. This feature is consistent with the neutral theory of molecular evolution. The rate of nonsynonymous substitution for human influenza viruses was three times the rate for swine influenza viruses. In particular, nonsynonymous substitutions at antigenic sites occurred less frequently in swine than in humans. The difference in the rate of nonsynonymous substitution between swine and human influenza viruses can be explained by the different degrees of functional constraint operating on the amino acid sequence of the HA in both hosts.     

Molecular evolution of the major capsid protein VP1 of enterovirus 70.

Nucleotide sequences of the genome RNA encoding capsid protein VP1 (918 nucleotides) of 18 enterovirus 70 (EV70) isolates collected from various parts of the world in 1971 to 1981 were determined, and nucleotide substitutions among them were studied. The genetic distances between isolates were calculated by the pairwise comparison of nucleotide difference. Regression analysis of the genetic distances against time of isolation of the strains showed that the synonymous substitution rate was very high at 21.53 x 10(-3) substitution per nucleotide per year, while the nonsynonymous rate was extremely low at 0.32 x 10(-3) substitution per nucleotide per year. The rate estimated by the average value of synonymous and nonsynonymous substitutions (W.-H. Li, C.-C. Wu, and C.-C. Luo, Mol. Biol. Evol. 2:150-174, 1985) was 5.00 x 10(-3) substitution per nucleotide per year. Taking the average value of synonymous and nonsynonymous substitutions as genetic distances between isolates, the phylogenetic tree was inferred by the unweighted pairwise grouping method of arithmetic average and by the neighbor-joining method. The tree indicated that the virus had evolved from one focal place, and the time of emergence was estimated to be August 1967 +/- 15 months, 2 years before first recognition of the pandemic of acute hemorrhagic conjunctivitis. By superimposing every nucleotide substitution on the branches of the phylogenetic tree, we analyzed nucleotide substitution patterns of EV70 genome RNA. In synonymous substitutions, the proportion of transitions, i.e., C<==>U and G<==>A, was found to be extremely frequent in comparison with that reported on other viruses or pseudogenes. In addition, parallel substitutions (independent substitutions at the same nucleotide position on different branches, i.e., different isolates, of the tree) were frequently found in both synonymous and nonsynonymous substitutions. These frequent parallel substitutions and the low nonsynonymous substitution rate despite the very high synonymous substitution rate described above imply a strong restriction on nonsynonymous substitution sites of VP1, probably due to the requirement for maintaining the rigid icosahedral conformation of the virus.     

Circumsporozoite Protein Genes of Malaria Parasites (Plasmodium Spp.): Evidence for Positive Selection on Immunogenic Regions

The circumsporozoite (CS) protein is a cell surface protein of the sporozoite, the stage of the life cycle of malaria parasites (Plasmodium spp.) that infects the vertebrate host. Analysis of DNA sequences supports the hypothesis that in Plasmodium falciparum, positive Darwinian selection favors diversity in the T-cell epitopes (peptides presented to T cells by host MHC molecules) of the CS protein. In gene regions encoding T cell epitopes of P. falciparum, the rate of nonsynonymous nucleotide substitution is significantly higher than that of synonymous substitution, whereas this is not true of other gene regions. Furthermore nonsynonymous nucleotide substitutions in these regions cause a change of amino acid residue charge significantly more frequently than expected by chance. By contrast, in Plasmodium cynomolgi, the same regions show no evidence of positive selection, and residue charge is conserved. The CS protein has a central repeat region, which is the target of host antibodies. In P. falciparum, the amino acid sequence of the repeat region is conserved within and between alleles. In P. cynomolgi, on the other hand, there is evidence that positive selection has favored evolution of two different repeat types within a given allele.     

A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes

A new method is proposed for estimating the number of synonymous and nonsynonymous nucleotide substitutions between homologous genes. In this method, a nucleotide site is classified as nondegenerate, twofold degenerate, or fourfold degenerate, depending on how often nucleotide substitutions will result in amino acid replacement; nucleotide changes are classified as either transitional or transversional, and changes between codons are assumed to occur with different probabilities, which are determined by their relative frequencies among more than 3,000 changes in mammalian genes. The method is applied to a large number of mammalian genes. The rate of nonsynonymous substitution is extremely variable among genes; it ranges from 0.004 X 10(-9) (histone H4) to 2.80 X 10(-9) (interferon gamma), with a mean of 0.88 X 10(-9) substitutions per nonsynonymous site per year. The rate of synonymous substitution is also variable among genes; the highest rate is three to four times higher than the lowest one, with a mean of 4.7 X 10(-9) substitutions per synonymous site per year. The rate of nucleotide substitution is lowest at nondegenerate sites (the average being 0.94 X 10(-9), intermediate at twofold degenerate sites (2.26 X 10(-9)). and highest at fourfold degenerate sites (4.2 X 10(-9)). The implication of our results for the mechanisms of DNA evolution and that of the relative likelihood of codon interchanges in parsimonious phylogenetic reconstruction are discussed.     

Evolution of nucleic acids coding for ribonucleases: the mRNA sequence of mouse pancreatic ribonuclease

The cDNA of mouse pancreatic mRNA has been cloned. After the library was screened with a rat ribonuclease cDNA probe, the positive clones were isolated and sequenced. There were no differences from the previously determined protein sequence. The mRNA codes for a preribonuclease of 149 amino acid residues including a signal peptide of 25 amino acids. The 3' noncoding region has a length of 260 bp, and the total mRNA length is approximately 940 bp. Comparison with the rat pancreatic ribonuclease sequence showed a high rate of nucleotide substitution. Within the coding region, nonsynonymous and synonymous substitution rates are 4.3 X 10(-9) and 15 X 10(-9) nucleotide substitutions/site/year, respectively. The latter value is one of the highest rates observed in the molecular evolution of mammalian nuclear genes. In the signal sequences the synonymous substitution rate is much lower and about the same as the nonsynonymous rate. Signal sequences of other mouse and rat proteins also exhibit little difference between synonymous and nonsynonymous rates. The sequences of rat and mouse pancreatic ribonuclease messengers were compared with those of bovine pancreatic, seminal, and brain ribonuclease. While the 3' noncoding regions of rat and mouse are very similar, as are those of the three bovine messengers, there is no significant similarity between both rodent and the three bovine messengers for the greater part of these regions. There is a duplication of approximately 50 nucleotides in the 3' noncoding region of the bovine messengers, with a region rich in A and C in between. The presence of this structural feature may be correlated with recent gene duplications that have occurred in the bovine genome.     

Protein Evolutionary Rates Correlate with Expression Independently of Synonymous Substitutions in Helicobacter pylori

In free-living microorganisms, such as Escherichia coli and Saccharomyces cerevisiae, both synonymous and nonsynonymous substitution frequencies correlate with expression levels. Here, we have tested the hypothesis that the correlation between amino acid substitution rates and expression is a by-product of selection for codon bias and translational efficiency in highly expressed genes. To this end, we have examined the correlation between protein evolutionary rates and expression in the human gastric pathogen Helicobacter pylori, where the absence of selection on synonymous sites enables the two types of substitutions to be uncoupled. The results revealed a statistically significant negative correlation between expression levels and nonsynonymous substitutions in both H. pylori and E. coli. We also found that neighboring genes located on the same, but not on opposite strands, evolve at significantly more similar rates than random gene pairs, as expected by co-expression of genes located in the same operon. However, the two species differ in that synonymous substitutions show a strand-specific pattern in E. coli, whereas the weak similarity in synonymous substitutions for neighbors in H. pylori is independent of gene orientation. These results suggest a direct influence of expression levels on nonsynonymous substitution frequencies independent of codon bias and selective constraints on synonymous sites. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Nicolas Galtier]     

Nucleic acid composition,codon usage,and the rate of synonymous substitution in protein-coding genes

Summary Based on the rates of synonymous substitution in 42 protein-codin gene pairs from rat and human, a correlation is shown to exist between the frequency of the nucleotides in all positions of the codon and the synonymous substitution rate. The correlation coefficients were positive for A and T and negative for C and G. This means that AT-rich genes accumulate more synonymous substitutions than GC-rich genes. Biased patterns of mutation could not account for this phenomenon. Thus, the variation in synonymous substitution rates and the resulting unequal codon usage must be the consequence of selection against A and T in synonymous positions. Most of the varition in rates of synonymous substitution can be explained by the nucleotide composition in synonymous positions. Codon-anticodon interactions, dinucleotide frequencies, and contextual factors influence neither the rates of synonymous substitution nor codon usage. Interestingly, the nucleotide in the second position of codons (always a nonsynonymous position) was found to affect the rate of synonymous substitution. This finding links the rate of nonsynonymous substitution with the synonymous rate. Consequently, highly conservative proteins are expected to be encoded by genes that evolve slowly in terms of synonymous substitutions, and are consequently highly biased in their codon usage.