新城疫病毒分离株HN和P基因分子进化及其相关研究

为探讨新城疫病毒(Newcastle disease virus,NDV)血凝素-神经氨酸酶(HN)和磷蛋白(P)基因遗传特性以及相互关系,将1997～2005年间国内分离到12株NDV毒株,分别进行HN和P基因克隆测序,结合15个已发表的国内外不同时期的NDV毒株HN和P基因,计算所有毒株HN和P基因的不同核苷酸和氨基酸片段进化距离,利用统计软件进行了不同片段间进化距离的方差分析,HN或P基因核苷酸进化距离与毒株分离时间、HN或P基因片段与其全长间以及HN和P基因全长间的相关分析.统计分析显示:NDV HN或P基因不同核苷酸和氨基酸序列片段变异程度不一样;不同毒株间HN或P基因片段与其全长间以及HN和P基因全长间无论是核苷酸还是氨基酸遗传变异高度相关.以上说明,NDV HN和P基因虽以不同的方式进化,但是HN和P基因遗传变异的趋势是相同的.HN和P基因的变异与分离时间有一定的联系.     

犬圆环病毒全基因组克隆与序列分析

犬圆环病毒(Dog circovirus,DogCV)是近年来新发现的一种哺乳动物圆环病毒。为研究DogCV中国流行毒株基因组特性和遗传进化,以犬血清样品提取的DNA为模板,采用重叠PCR技术首次获得DogCV中国流行毒株的全基因组序列,命名为JZ98/2014。序列分析表明JZ98/2014全基因组长度为2063nt,编码3个主要开放阅读框:ORF V1(Rep蛋白,303个氨基酸)、ORF C1(Cap蛋白,270个氨基酸)和ORF C2(106个氨基酸)。同源性比较显示JZ98/2014与美国、欧洲流行毒株的全基因组序列同源性为82.1%-89.5%,而Rep和Cap基因同源性分别为82.1%-89.5%和84.6%-89.1%。遗传进化树分析显示,目前世界流行的DogCV存在多个分支,而JZ98/2014与美国、欧洲流行毒株处于不同分支。     

陕西省9株乙脑病毒的分离及E基因序列分析

分析陕西省分离的9株乙脑病毒基因组序列特征。使用乙脑病毒全基因组序列测定引物进行RT-PCR扩增,扩增产物进行测序,拼接后获得基因组序列。利用MEGA 4.1、MegAlin、MEGA7.0等软件进行毒株的系统进化分析,并与P3株、减毒活疫苗SA14-14-2株及覆盖5个基因型别的其他乙脑病毒进行E基因序列比对。9株分离株3株分离自猪舍、6株分离自羊舍,其中4株获得全基因组序列,5株测得E基因序列。基于E基因序列进行毒株核苷酸、氨基酸同源性比较,结果显示分离株均与基因I型GI-b亚型毒株核苷酸和氨基酸同源性最高,核苷酸同源性范围为96.5%～99.7%、氨基酸同源性范围99.2%～100.0%;与SA14-14-2株核苷酸同源性范围为87.5%～88.9%、氨基酸同源性范围96.3%～97.2%;与P3株核苷酸同源性范围为87.6%～88.1%、氨基酸同源性范围96.7%～97.6%。分析09年（陕南地区）分离株与18年（关中地区）分离株的E基因核苷酸差异率为1.8%～2.9%、氨基酸差异率为0%～0.8%。陕西省自然界中循环的乙脑病毒以基因Ⅰ型为主,与P3株在抗原毒力关键位点无差异,...     

利用分子进化树计算蛋白质的特异进化速率

本文提出了一种计算蛋白质绝对进化距离和进化速率的方法,它根据现有同源蛋白质的序列构建分子进化树,并推断进化过程中各结点处的共同祖先序列,根据某成员与某结点处共同祖先序列的氨基酸差异百分率,计算该蛋白质序列的特异进化距离和进化速率。比较我们的算法和Ｄａｙｈｏｆｆ等的模拟统计方法表明,我们的算法在一定范围内是正确的。结合计算哺乳动物红细胞生成素的进化速率,讨论了本算法在分子进化研究中的应用。     

中华蜜蜂囊状幼虫病毒北京分离株VP1蛋白基因序列特征及原核表达

【目的】研究中华蜜蜂囊状幼虫病毒（Chinese sacbrood virus, CSBV）VP1蛋白的分子进化特征及遗传多样性。【方法】利用RT-PCR方法,克隆了8株CSBV北京分离株VP1蛋白的基因编码区。【结果】序列分析表明,VP1蛋白基因编码区开放阅读框长945 bp,编码315个氨基酸,推测编码蛋白的相对分子量和等电点分别为35.42 kDa和9.23,具有亲水性和免疫原性。序列同源性分析表明,不同年份CSBV北京分离株VP1蛋白氨基酸序列间差异较小,仅个别氨基酸存在差异。北京分离株与辽宁分离株及越南分离株VP1核苷酸序列一致性达93%,与印度及韩国分离株VP1核苷酸序列一致性达92%,与英国分离株VP1核苷酸序列一致性最低,为88%。序列分析同时表明,CSBV北京分离株VP1蛋白序列存在特有的序列特征,同其他地区分离株比较,北京分离株VP1蛋白序列中存在着氨基酸的插入突变。序列替换率分析表明,亚洲型分离株间序列替换率低于亚洲分离株与欧洲分离株间的替换率。构建原核表达载体pEASY-E1-VP1,经IPTG诱导,CSBV VP1蛋白在大肠杆菌Escherichia coli BL21（DE3）pLysS菌株中表达。【结论】本研究提示CSBV不同分离株基因序列存在变异,结果为进一步研究CSBV致病性分化的分子机理奠定了基础。     

一株H5N2亚型鹦鹉源禽流感病毒分子进化分析

2005年在广东进行流行病学调查时分离到一株鹦鹉源禽流感病毒,经鉴定为H5N2亚型禽流感病毒(A/Parrot/Guangdong/268/2005)。该毒株的HA裂解位点附近的氨基酸序列为RETRGLF,只含有一个碱性氨基酸,符合低致病性禽流感病毒的HA裂解位点附近氨基酸序列的分子特征;与H5N2亚型禽流感代表毒株相比,该毒株HA和NA基因的糖基化位点、HA基因的受体结合位点编码区、NA基因的耐药性位点均未发生变异。将该毒株全基因组序列与GenBank已公布的19株H5N2亚型禽流感病毒株的相应序列进行比较分析并绘制系统进化树后发现:其与低致病性禽流感毒株A/Pheasant/NJ/1355/1998(H5N2)-like的亲缘关系最近,位于以A/Chicken/Pennsylvania/1/1983(H5N2)为代表的美洲进化分支。     

近20年中国部分地区鸡源H9N2亚型禽流感病毒HA基因遗传演化及其变异频率北大核心CSCD

【目的】通过比较不同时期的H9N2亚型禽流感流行毒株HA基因的分子特征和变异频率,揭示免疫压力下病毒的遗传演化趋势。【方法】选取源于课题组的40株鸡源H9N2毒株,以及从Gen Bank下载的136株中国鸡源H9N2流行毒株和7株经典毒株的序列,利用Lasergen 7.1和MEGA 5.1等软件,对其HA基因进行系统演化、分子特征和变异频率分析。【结果】系统发育分析表明,近20年的鸡源H9N2流行株分属于BJ94、Y280和S2等谱系,优势流行株的分布与年代密切相关。氨基酸序列比较显示,H9N2病毒不同谱系之间具有各自的特征,且存在着明显的氨基酸变异积累。以Ck/BJ/1/1994 HA基因为参照,1994–2014年间,H9N2流行株核苷酸和氨基酸的年均进化率分别为5.73×10^(–3)和4.25×10^(–3)。其中,2011–2014年的核苷酸(氨基酸)年均进化率为6.35×10^(–3)(5.32×10^(–3)),明显高于2006–2010年5.22×10^(–3)(3.70×10^(–3)),更显著高于疫苗推广初期1999–2005年的0.74×10^(–3)(0.50×10^(–3))。【结论】H9N2疫苗株和流行毒株的不匹配是病毒变异频率加快的重要原因。     

嗜肝病毒家族成员间进化关系的初步研究

比较分析了嗜肝病毒家族四成员间的全基因核苷酸序列,得出了它们之间的相对进化距离,DHBV 分野最早,GSHV,WHV 次之,HBV 最晚;其进化方式属趋 异性进化.进一步比较各 ORF 的保守性,发现各 ORF 间及同一 ORF 内的不同区域核苷酸发生替换的频率有明显差别;讨论了各 ORF 的特点及意义.     

传染性支气管炎病毒(IBV)国内分离毒株的分子流行病学分析

选择9个来自我国不同地区的传染性支气管炎病毒(IBV)地方毒株,用逆转录-聚合酶链反应(RT-PCR)分别扩增出它们的完整S1基因,再将各扩增产物与T载体连接转化,筛选出相应毒株的阳性克隆,采用双脱氧链终止法测定基因序列。将这些国内分离毒株、我国常用疫苗毒株和国际上其他血清型的代表参考毒株一起,用DNAstar软件分别进行基因与氨基酸系统发生进化关系分析,结果将这些毒株分成3群：Ⅰ群内的毒株与疫苗毒株H120和M41的同源性很高,Ⅱ群内的毒株与其他毒株的同源性较低,Ⅲ群内的毒株与疫苗毒株有一定的同源性。3个群的毒株以疫苗毒株H120为核心,形成进化距离的梯度,提示弱毒疫苗的使用可能是我国IBV流行毒株的主要来源之一。分析结果还显示,我国IBV的分离毒株没有明显的时间和地理分布规律可循。     

2010年云南省1株ECHO-9病毒全基因序列的遗传特性分析

为了解云南省手足口病患者标本中分离到的一株ECHO-9病毒基因组特征,对2010年云南省ECHO-9病毒分离株MSH-KM812-2010全基因组序列测序,并与GenBank中其它ECHO-9病毒株基因组序列比对和分析。MSH-KM812-2010基因组长为7 424bp,编码2 203个氨基酸,其结构基因区与其它ECHO-9病毒核苷酸和氨基酸的同源性高于其它型别肠道病毒;而在非结构区,与其它肠道病毒血清型的同源性高于ECHO-9病毒株。VP1系统进化分析显示ECHO-9病毒株可形成A、B和C三个进化分支,MSH-KM812-2010株以及其它中国分离株属于C簇。三个分支之间核苷酸序列的差异大于15.0%可将ECHO-9病毒分为三个基因型。通过重组检测软件3(RDP3)和基本局部比对搜索工具(BLAST)比对分析,发现在非结构区可能存在重组。本文首次对我国分离的ECHO-9病毒全基因组序列的测定和分析,对了解ECHO-9病毒遗传特性具有重要意义。     

Mapping of sequences required for mouse neurovirulence of poliovirus type 2 Lansing.

Intracerebral inoculation of mice with poliovirus type 2 Lansing induces a fatal paralysis, while most other poliovirus strains are unable to cause disease in the mouse. To determine the molecular basis for Lansing virus neurovirulence, we determined the complete nucleotide sequence of the Lansing viral genome from cloned cDNA. The deduced amino acid sequence was compared with that of two mouse-avirulent strains. There are 83 amino acid differences between the Lansing and Sabin type 2 strain and 179 differences between the Lansing and Mahoney type 1 strain scattered throughout the genome. To further localize Lansing sequences important for mouse neurovirulence, four intertypic recombinants were isolated by exchanging DNA restriction fragments between the Lansing 2 and Mahoney 1 infectious poliovirus cDNA clones. Plasmids were transfected into HeLa cells, and infectious recombinant viruses were recovered. All four recombinant viruses, which contained the Lansing capsid region and different amounts of the Mahoney genome, were neurovirulent for 18- to 21-day-old Swiss-Webster mice by the intracerebral route. The genome of neurovirulent recombinant PRV5.1 contained only nucleotides 631 to 3413 from Lansing, encoding primarily the viral capsid proteins. Therefore, the ability of Lansing virus to cause paralysis in mice is due to the viral capsid. The Lansing capsid sequence differs from that of the mouse avirulent Sabin 2 strain at 32 of 879 amino acid positions: 1 in VP4, 5 in VP2, 4 in VP3, and 22 in VP1.     

Complete nucleotide sequences of all three poliovirus serotype genomes. Implication for genetic relationship, gene function and antigenic determinants

The complete nucleotide sequences of the genomes of the type 2 ( P712 , Ch, 2ab ) and type 3 (Leon 12a1b ) poliovirus vaccine strains were determined. Comparison of the sequences with the previously established genome sequence of type 1 (LS-c, 2ab ) poliovirus vaccine strain revealed that 71% of the nucleotides in the genome RNAs were common, that the 5' and 3' termini of the genomes were highly homologous, and that more than 80% of the nucleotide differences in the coding region occurred in the third letter position of in-phase codons, resulting in a low frequency of amino acid difference. These results strongly suggested that the serotypes of poliovirus derived from a common prototype. A comparison of the amino acid sequences predicted from the genome sequences showed highest variation in the capsid protein region, whereas non-structural proteins are highly conserved. Initiation of polyprotein synthesis occurs in all three strains more than 740 nucleotides downstream from the 5' end. An analysis of the non-coding region suggests that small peptides that could potentially originate from this region are conserved. The amino acid sequences immediately surrounding the cleavage signals, however, show a higher than average degree of variation. The analysis of the amino acid sequences of the capsid protein VP1 of all serotypes has led to the prediction of potential antigenic sites on the virion involved in neutralization.     

Identification of a consistent pattern of mutations in neurovirulent variants derived from the sabin vaccine strain of poliovirus type 2.

Complete nucleotide sequencing of the RNAs of two unrelated neurovirulent isolates of Sabin-related poliovirus type 2 revealed that two nucleotides and one amino acid (amino acid 143 in the major capsid protein VP1) consistently departed from the sequences of the nonneurovirulent poliovirus type 2 712 and Sabin vaccine strains. This pattern of mutation appeared to be a feature common to all neurovirulent variants of poliovirus type 2.     

Similarities in protein amino acid composition in connection with principles of protein evolution

The amino acid composition of human alcohol dehydrogenase (ADH) was compared with alcohol dehydrogenases from different organisms and with other proteins. Similar amino acid sequences in human ADH (template protein) and in other proteins were determined by means of an original computer program. Analysis of amino acid motifs reveals that the ADHs from evolutionary more close organisms have more common amino acid sequences. The quantity measure of amino acid similarity was the number of similar motifs in analyzed protein per protein length. This value was measured for ADHs and for different proteins. For ADHs, this quotient was higher than for proteins with different functions; for vertebrates it correlated with evolutionary closeness. The similar operation of motif comparison was made with the help of program complex “MEME”. The analysis of ADHs revealed 4 motifs common to 6 of 10 tested organisms and no such motifs for proteins of different function. The conclusion is that general amino composition is more important for protein function than amino acid order and for enzymes of similar function it better correlates with evolutionary distance between organisms.     

Evolution of the nucleoprotein gene of influenza A virus

Nucleotide sequences of 24 nucleoprotein (NP) genes isolated from a wide range of hosts, geographic regions, and influenza A virus serotypes and 18 published NP gene sequences were analyzed to determine evolutionary relationships. The phylogeny of NP genes was determined by a maximum-parsimony analysis of nucleotide sequences. Phylogenetic analysis showed that NP genes have evolved into five host-specific lineages, including (i) Equine/Prague/56 (EQPR56), (ii) recent equine strains, (iii) classic swine (H1N1 swine, e.g., A/Swine/Iowa/15/30) and human strains, (iv) gull H13 viruses, and (v) avian strains (including North American, Australian, and Old World subgroups). These NP lineages match the five RNA hybridization groups identified by W. J. Bean (Virology 133:438-442, 1984). Maximum nucleotide differences among the NPs was 18.5%, but maximum amino acid differences reached only 10.8%, reflecting the conservative nature of the NP protein. Evolutionary rates varied among lineages; the human lineage showed the highest rate (2.54 nucleotide changes per year), followed by the Old World avian lineage (2.17 changes per year) and the recent equine lineage (1.22 changes per year). The per-nucleotide rates of human and avian NP gene evolution (1.62 x 10(-3) to 1.39 x 10(-3) changes per year) are lower than that reported for human NS genes (2.0 x 10(-3) changes per year; D. A. Buonagurio, S. Nakada, J. D. Parvin, M. Krystal, P. Palese, and W. M. Fitch, Science 232:980-982, 1986). Of the five NP lineages, the human lineage showed the greatest evolution at the amino acid level; over a period of 50 years, human NPs have accumulated 39 amino acid changes. In contrast, the avian lineage showed remarkable conservatism; over the same period, avian NP proteins changed by 0 to 10 amino acids. The specificity of the H13 NP in gulls and its distinct evolutionary separation from the classic avian lineage suggests that H13 NPs may have a large degree of adaptation to gulls. The presence of avian and human NPs in some swine isolates demonstrates the susceptibility of swine to different virus strains and supports the hypothesis that swine may serve as intermediates for the introduction of avian influenza virus genes into the human virus gene pool. EQPR56 is relatively distantly related to all other NP lineages, which suggests that this NP is rooted closest to the ancestor of all contemporary NPs. On the basis of estimation of evolutionary rates from nucleotide branch distances, current NP lineages are at least 100 years old, and the EQPR56 NP is much older.(ABSTRACT TRUNCATED AT 400 WORDS)     

Guanidine-selected mutants of poliovirus: mapping of point mutations to polypeptide 2C.

Sequence analysis of the genomic RNA of interstrain guanidine-resistant and antibody-resistant variant recombinants of poliovirus type 1 mapped the resistance of mutants capable of growth in 2.0 mM guanidine hydrochloride to a region located 3' of nucleotide 4444. This region of the viral genome specifies the nonstructural protein 2C. The sequence of genomic RNA encoding 2C from six independently isolated mutants resistant to 2.0 mM guanidine was determined. All six isolates contained a mutation in 2C at the same position in all cases, resulting in two types of amino acid changes. Dependent mutants were examined and found to contain two amino acid changes each within 2C. Mutants resistant to 0.53 mM guanidine were isolated and found to lack the mutations seen in variants resistant to 2.0 mM guanidine. A comparison of the amino acid sequences of the 2C proteins of poliovirus, foot-and-mouth disease virus, rhinovirus types 2 and 14, and encephalomyocarditis virus revealed a strong homology over regions totaling 115 residues. All of the mutations observed in guanidine-selected mutants were contained within this region. The amino acid region containing the mutations observed in poliovirus mutants resistant to 2.0 mM guanidine was compared with the homologous region in the other picornaviruses; a strong correlation was found between the amino acid present at this position and the sensitivity of the virus to 2.0 mM guanidine.     

Natural genetic exchanges between vaccine and wild poliovirus strains in humans

In a previous study of poliovirus vaccine-derived strains isolated from patients with vaccine-associated paralytic poliomyelitis (VAPP) (9, 11), we reported that a high proportion (over 50%) of viruses had a recombinant genome. Most were intertypic vaccine/vaccine recombinants. However, some had restriction fragment length polymorphism (RFLP) profiles different from those of poliovirus vaccine strains. We demonstrate here that five such recombinants, of 88 VAPP strains examined, carried sequences of wild (nonvaccine) origin. To identify the parental wild donor of these sequences, we used RFLP profiles and nucleotide sequencing to look for similarity in the 3D polymerase-coding region of 61 wild, cocirculating poliovirus isolates (43 type 1, 16 type 2, and 2 type 3 isolates). In only one case was the donor identified, and it was a wild type 1 poliovirus. For the other four vaccine/wild recombinants, the wild parent could not be identified. The possibility that the wild sequences were of a non-poliovirus-enterovirus origin could not be excluded. Another vaccine/wild recombinant, isolated in Belarus from a VAPP case, indicated that the poliovirus vaccine/wild recombination is not an isolated phenomenon. We also found wild polioviruses (2 of 15) carrying vaccine-derived sequences in the 3' moiety of their genome. All these results suggest that genetic exchanges with wild poliovirus and perhaps with nonpoliovirus enteroviruses, are also a natural means of evolution for poliovirus vaccine strains.     

Evolution of the Borrelia burgdorferi outer surface protein OspC.

The genes coding for outer surface protein OspC from 22 Borrelia burgdorferi strains isolated from patients with Lyme borreliosis were cloned and sequenced. For reference purposes, the 16S rRNA genes from 17 of these strains were sequenced after being cloned. The deduced OspC amino acid sequences were aligned with 12 published OspC sequences and revealed the presence of 48 conserved amino acids. On the basis of the alignment, OspC could be divided into an amino-terminal relatively conserved region and a relatively variable region in the central portion. The distance tree obtained divided the ospC sequences into three groups. The first group contained ospC alleles from all (n = 13) sensu stricto strains, the second group contained ospC alleles from seven Borrelia afzelii strains, and the third group contained ospC alleles from five B. afzelii and all (n = 9) Borrelia garinii strains. The ratio of the mean number of synonymous (dS) and nonsynonymous (dN) nucleotide substitutions per site calculated for B. burgdorferi sensu stricto, B. garinii, and B. afzelii ospC alleles suggested that the polymorphism of OspC is due to positive selection favoring diversity at the amino acid level in the relatively variable region. On the basis of the comparison of 16S rRNA gene sequences, Borrelia hermsii is more closely related to B. afzelii than to B. burgdorferi sensu stricto and B. garinii. In contrast, the phylogenetic tree obtained for the B. hermsii variable major protein, Vmp33, and 18 OspC amino acid sequences suggested that Vmp33 and OspC from B. burgdorferi sensu stricto strains share a common evolutionary origin.     

Unbiased estimation of evolutionary distance between nucleotide sequences

A new algorithm for estimating the number of nucleotide substitutions per site (i.e., the evolutionary distance) between two nucleotide sequences is presented. This algorithm can be applied to many estimation methods, such as Jukes and Cantor's method, Kimura's transition/transversion method, and Tajima and Nei's method. Unlike ordinary methods, this algorithm is always applicable. Numerical computations and computer simulations indicate that this algorithm gives an almost unbiased estimate of the evolutionary distance, unless the evolutionary distance is very large. This algorithm should be useful especially when we analyze short nucleotide sequences. It can also be applied to amino acid sequences, for estimating the number of amino acid replacements.      

Linkage of Amino Acid Variation and Evolution of Human Immunodeficiency Virus Type 1 gp120 Envelope Glycoprotein (Subtype B) with Usage of the Second Receptor

To clarify the relationship between the amino acid variations of the gp120 of human immunodeficiency virus type 1 (HIV-1) and the chemokine receptors that are used as the second receptor for HIV, we evaluated amino acid site variation of gp120 between the X4 strains (use CXCR4) and the R5 strains (use CCR5) from 21 sequences of subtype B. Our analysis showed that residues 306 and 322 in the V3 loop and residue 440 in the C4 region were associated with usage of the second receptor. The polymorphism at residue 440 is clearly associated with the usage of the second receptor: The amino acid at position 440 was a basic amino acid in the R5 strains, and a nonbasic and smaller amino acid in the X4 strains, while the V3 loop of the X4 strains was more basic than that of the R5 strains. This suggests that residue 440 in the C4 region, which is close to the V3 loop in the three-dimensional structure, is critical in determining which second receptor is used. Analysis of codon frequency suggests that, in almost all cases, the difference at residue 440 between basic amino acids in the R5 strains and nonbasic amino acids in the X4 strains could be due to a single nucleotide change. These findings predict that the evolutionary changes in amino acid residue 440 may be correlated with evolutionary changes in the V3 loop. One possibility is that a change in electric charge at residue 440 compensates for a change in electric charge in the V3 loop. The amino acid polymorphism at position 440 can be useful to predict the cell tropism of a strain of HIV-1 subtype B.