马铃薯Y病毒福建分离物P1基因的分子变异和结构特征

为揭示马铃薯Y病毒(Potato virus Y, PVY)P1基因的分子变异及结构特征, 并查明福建分离物P1基因的变异来源。文章设计一对简并引物从感染PVY的马铃薯病叶扩增、克隆获得福建分离物P1基因的cDNA全长序列, 分析其核苷酸序列和氨基酸序列的特征, 并采用贝叶斯法重建P1基因的系统发育树。结果显示, 12个福建分离物均成功扩增出预期大小(约915 bp)的特异性片段, P1基因的核苷酸序列与已报道的PVY不同株系的核苷酸序列一致性为73%~99%; QK44、XT02、XT08、LH05分离物的309 nt位置均检测到一个强烈的重组信号。12个PVY分离物中, P1蛋白有85个变异的氨基酸位点, 表明其高度变异; P1蛋白的第41~275位是一个高度保守的结构功能域, 含有3个保守的活性位点(H₁₉₂、D₂₀₁、V₂₃₅); 系统发育分析显示, 福建分离物形聚为3种不同的簇, 其对应的卷曲结构(Coiled-coil domain)和空间结构也不相同, 表明不同株系型的P1基因在系统发育关系上分化较大。该研究结果表明PVY P1基因在核苷酸序列、氨基酸序列以及空间结构具有高度变异性, 但行使P1蛋白功能的活性位点所在的氨基酸(H₁₉₂、D₂₀₁和V₂₃₅)均高度保守; 福建分离物P1基因的变异源主要来自碱基突变和基因重组。     

不同寄主来源的马铃薯Y病毒群体遗传结构的比较分析

文章以马铃薯Y病毒(Potato virus Y, PVY) P3和pipo基因为分子标记,比较分析烟草和马铃薯两个寄主的PVY遗传多样性和群体分化,并评估突变、选择、基因流等遗传力所起的作用。通过P3和pipo基因计算获得的烟草和马铃薯群体分化指数F_ST分别为0.116和0.120,且统计检验差异显著,表明烟草和马铃薯寄主的PVY之间中度分化。变异分析结果显示,烟草分离物P3和pipo基因的核苷酸序列一致性分别为85.2%~100%和76.5%~100%,而马铃薯分离物的P3和pipo基因的核苷酸序列一致性分别为95.7%~100%和93.0%~100%,表明烟草PVY变异程度高于马铃薯。同时,P3基因内检测到大量的净化选择位点,表明该基因大部分位点的变异为有害突变,在进化过程中被剔除。此外,P3基因内还检测到两个显著正向选择位点,表明这两个位点的变异为有益突变,有利于病毒的生存竞争。在pipo基因中未检测到显著的选择位点,表明该基因上的变异基本不受自然选择影响。通过P3和pipo基因计算烟草和马铃薯群体间的基因流Nm值分别为1.91和1.83,表明这两个群体间存在较强的基因交流。以上结果表明,来源于烟草和马铃薯寄主的PVY遗传差异显著,突变、自然选择以及基因流都影响两者的遗传多样性及遗传分化程度。     

全基因组序列的系统发育与重组分析鉴定我国马铃薯Y病毒株系

为明确我国马铃薯Y病毒(Potato virus Y,PVY)株系组成及其在不同寄主上的分布特征,文章对31条PVY中国分离物的基因组编码区进行了系统发育及重组分析。系统发育分析结果显示,最大似然法(maximum likelihood,ML)和邻接法(neighbor-joining,NJ)重建的系统发育拓扑结构基本一致,GF_YL20、FZ10、DF、SD1、1116和9703-3等6个分离物均未与已知株系相聚,其他株系来源相同的分离物以较高的自展值相聚成簇,显示出较强的株系特异性。重组分析显示,31个PVY分离物中,共有28个PVY分离物均检测到具有显著的重组信号,这些分离物主要包括PVY~(N-Wi)、PVY~(NTN)和PVY~(NTN-NW)三种株系。进一步分析显示,分离物1104、DF、SD1、1116和9703-3为未报道过的PVY新重组分离物。株系统计结果显示PVY重组株系已上升为优势株系,生产上需要密切注意这些株系,尤其新重组分离物的发展动态。以上结果表明,综合应用PVY全基因组序列的系统发育和重组分析,可以实现PVY株系的准确鉴定。     

马铃薯Y病毒N株HC-pro基因的序列分析与表达

马铃薯Ｙ病毒(Potato virus Y,PVY)是马铃薯Ｙ病毒属(Potyvirus)的典型成员,是烟草、马铃薯和辣椒等茄科作物病毒病的主要病原之一,分布于世界各地.近年烟草上的 PVY在全国范围内特别是陕西、黄淮和东北烟区流行呈上升趋势,且以PVY坏死株系（PVY-N)为主,重病田发生率达40%以上,并常与烟草花叶病毒、黄瓜花叶病毒混合侵染,严重影响烟草品质造成巨大的经济损失.     

辣椒脉斑驳病毒文昌分离物基因组测序及分析

辣椒脉斑驳病毒(Chilli veinal mottle virus,ChiVMV)近年来严重危害海南黄灯笼辣椒的产量,开展针对该病毒的研究以求达到对其有效防控迫在眉睫.本研究分7段克隆了ChiVMV文昌分离物(ChiVMV Wenchang isolate,ChiVMV-WC)的基因组.分段克隆的测序结果通过拼接,共获得ChiVMV-WC基因组9 717个核苷酸(nt)序列(GenBank登录号:GQ981316).在核苷酸序列的第164位~第9 270位存在一个大的开放阅读框(ORF),编码一个多聚蛋白(分子量为350.44 kD).近年来被证实的马铃薯Y病毒科的一个新的ORF (pretty interesting potyviridae ORF,PIPO)存在于ChiVMV-WC基因组核苷酸序列的第2 892位~第3 110位,编码一个由72个氨基酸聚合而成的多肽(分子量为8.26 kD).ChiVMV-WC与ChiVMV其它分离物的基因组核苷酸序列比较发现该病毒存在着较高的变异率.本研究通过与同属内病毒的同源性分析以及系统进化树分析,确定了ChiVMV-WC在生物进化史上的地位.     

表达PVY和PLRV双价外壳蛋白基因马铃薯的抗病性研究

表达马铃薯Y病毒(PVY)和马铃薯卷叶病毒(PLRV)双价外壳蛋白基因的马铃薯(Solanum tubero-sum L.)栽培品种“Favorita”和“虎头”,经摩擦接种PVY和用桃蚜接种PLRV后,观察症状并用ELISA测定病毒滴度。结果表明,两个品种转双价CP基因的各株系,接种病毒后表现无症状或症状轻微,其中PVY和PLRV平均滴度均较不转基因对照植株低。不同品种对PVY和PLRV的抗性比较表明,转双价CP基因的“Favorita”对PVY抗性较明显,而转双价CP基因的“虎头”则对PLRV抗性较对PVY抗性明显。不同转基因株系抗病毒水平不同。“Favorita”9个转双价CP基因株系中有6个株系PVY滴度较未转基因对照降低52.5％～90.0％,而“虎头”7个转双价CP基因株系中有4个株系PLRV含量较对照降低53.0％～98.0％。在抗性株系中还出现一些抗1种病毒或抗2种病毒的抗性较强的单株。     

马铃薯Y病毒P3N-PIPO酵母双杂交诱饵载体的构建及鉴定

P3N-PIPO是马铃薯Y病毒属(Potyvirus)病毒基因组中近年来新鉴定的编码蛋白。为研究P3N-PIPO在病毒与宿主互作过程中的功能,本研究以马铃薯Y病毒(potato virus Y,PVY)和马铃薯作为研究系统,用融合PCR技术扩增得到PVY编码蛋白P3N-PIPO的DNA序列,用于构建酵母双杂交诱饵质粒。为钓取马铃薯组织细胞内不同位置的互作蛋白,本研究构建了p BT3-N-P3N-PIPO、p BT3-C-P3N-PIPO、p BT3-STE-P3N-PIPO、p BT3-SUC-P3N-PIPO和p DHB1-P3N-PIPO 5个诱饵质粒,并进行自激活或毒性检测,最终获得p BT3-STE-P3N-PIPO、p BT3-SUC-P3N-PIPO和p DHB1-P3N-PIPO 3个可用于筛选马铃薯c DNA文库的诱饵质粒。     

马铃薯Y病毒多基因系统发育分析及其在株系鉴定中的应用

为实现马铃薯Y病毒(Potato virus Y, PVY)常见株系的快速鉴定,本文以PVY的P1、HC-pro、VPg和CP 4个基因为研究对象建立了快速准确的多基因联合体系。根据基因的不同组合建立5个不同数据集,分别进行系统发育分析,并通过贝叶斯标签关联显著性(Bayesian tip-association significance, BaTS)分析各数据集中代表分离物与株系的关联性,以确定实现PVY快速鉴定的最佳组合。不同数据集的系统发育及BaTS分析结果显示,除了联合P1、VPg和CP 3个基因数据集外,其他4个数据集均无法实现PVY常见株系的准确鉴定。采用不同建树方法对联合P1、VPg和CP 3个基因数据集比较分析显示,基于ML法和NJ法的系统发育树在拓扑结构上基本一致,均优于基于贝叶斯算法的最大分支置信(maximum clade credibility, MCC)树。同时,以HLJ26分离物为研究对象,对建立的多基因联合体系进行实际应用,结果显示该分离物与PVY^NTN-NW株系的3个分离物SYR-Ⅱ-2-8、SYR-Ⅱ-Be1和SYR-Ⅱ-DrH以高置信值聚为一亚簇,表明该分离物可能属于PVY^NTN-NW株系(SYR-Ⅱ型)。重组分析显示,HLJ26基因组存在4个潜在的重组信号,分别位于P1、HC-pro/P3、VPg和CP的5°-末端,与PVY^NTN-NW株系(SYR-Ⅱ型)的重组位点相一致,表明其属于PVY^NTN-NW株系(SYR-Ⅱ型)。同时,应用多重RT-PCR成功扩增出约为1000 bp和400 bp的2个特异性片段,与PVY^NTN-NW株系(SYR-Ⅱ型)的特异条带大小相一致。这些结果进一步支持了多基因联合体系的鉴定结果。联合P1、VPg和CP 3个基因数据集系统发育分析,可以实现PVY常见株系的准确鉴定。     

呼吸道合胞病毒B亚型分离株的G蛋白基因分析

对一株长春地区B亚型分离株(CC169)的G蛋白基因进行了 序列分析,结果表明：我国呼吸道合胞病毒(RSV)分离株CC169同RSV B亚型原型株CH18537的 核苷酸同源性为94%,核苷酸的有义突变率达65%。由核苷酸推导出氨基酸序列的同源性为894%,氨基酸的变异全部发生在胞外区,并主要集中在一个高度保守区的两端,胞内区和跨 膜区保守不变。氨基酸的变异导致了分离株既有糖基化位点的改变,又有蛋白长度的变异。 此外还初步探讨了我国RSV B 亚型分离株CC169的G蛋白基因同原型株之间的变异与疫苗研制 中的意义。     

一株引起严重急性呼吸道感染病例的A亚型人呼吸道合胞病毒全基因组基因特征分析

人呼吸道合胞病毒(Human respiratory syncytial virus,HRSV)是导致儿童急性呼吸道感染的最重要的呼吸道病毒之一。根据对单克隆抗体的反应,HRSV分为A、B两个亚型。为探讨严重急性呼吸道感染(Severe acute respi-ratory infection,SARI)病例中HRSV全基因组基因特征,本研究对2017年河南省漯河市住院SARI病例中检测到的1株HRSV A亚型病毒通过Sanger测序方法对其全基因组序列进行了测定和分析。通过Sequencher 5.4.5、MEGA 5.05、BioEdit 7.0.5等生物信息学软件进行序列拼接和比对,进行了基因亲缘性关系分析、氨基酸变异和糖基化位点分析。基于HRSV全基因组序列和11个单个蛋白基因序列构建的亲缘性关系分析结果提示本研究中检测到的这株HRSVA病毒(RSVAs/Luohe.Henan/CHN/42.17)属于ON1基因型,该型是我国近年流行的优势基因型。该病毒全基因组序列与35条全球代表株的核苷酸和氨基酸同源性分别为92.69%～99.82%和93.63%～99.67%;G蛋白编码区氨基酸变异最高,而F蛋白相对保守。糖基化位点分析发现,该病毒的F蛋白有6个N-糖基化位点,未发现O-糖基化位点,此结果与原型株long株相同;G蛋白N-糖基化位点有6个,O-糖基化位点为82个,而原型株long株有11个N-糖基化位点,15个O-糖基化位点。本研究对2017年河南省漯河市SARI病例中一株HRSVA病毒全基因组序列进行了测定,与世界其他地区报道的HRSVA亚型病毒全基因组序列进行了对比分析,揭示了SARI病例中我国HRSV优势流行ON1基因型病毒全基因组的核苷酸和氨基酸变异特征,以及G蛋白和F蛋白编码区糖基化情况,丰富了我国HRSV基因数据库,也为HRSV的核酸检测方法的建立、疫苗研发和预防性单克隆抗体的评价提供了核苷酸和氨基酸的基础数据。     

Evaluation of the internal transcribed spacer 2 (ITS2) as a molecular marker for phylogenetic inference using sequence and secondary structure information in blow flies (Diptera: Calliphoridae)

The internal transcribed spacer 2 (ITS2) is a small non-coding region located inside the nuclear ribosomal DNA cluster. ITS2 sequence variability is thought to be appropriate to differentiate species and for phylogenetic reconstructions analyses, which can be further improved if structural information is considered. We evaluated the potential of ITS2 as a molecular marker for phylogenetic inference in Calliphoridae (Diptera: Brachycera) using a broad range of inference methods and different substitution models, accounting or not for structural information. Sequence analyses revealed a hierarchically organized pattern of sequence variation and a small level of nucleotide substitution saturation. Intragenomic variation due to small sequence repeats was found mainly in the most variable domain (IV), but it has no significant impact on the phylogenetic signal at the species level. Inferred secondary structures revealed that GC pairs are more frequently found flanking bulges and loops regions in more conserved domains, thus ensuring structure stability. In the phylogenetic analyses, the use of substitution models accounting for structural information significantly improves phylogenetic inference in both neighbour-joining and Bayesian analyses, although the former provides limited resolution for dealing with highly divergent sequences. For Bayesian analyses, a significant improvement in likelihood was observed when considering structure information, although with small changes in topology and overall support, probably reflecting better evolutionary rates estimates. Based on these findings, ITS2 is a suitable molecular marker for phylogenetic analyses in Calliphoridae, at both species and generic level.     

Granule-bound starch synthase: structure, function, and phylogenetic utility

Interest in the use of low-copy nuclear genes for phylogenetic analyses of plants has grown rapidly, because highly repetitive genes such as those commonly used are limited in number. Furthermore, because low- copy genes are subject to different evolutionary processes than are plastid genes or highly repetitive nuclear markers, they provide a valuable source of independent phylogenetic evidence. The gene for granule-bound starch synthase (GBSSI or waxy) exists in a single copy in nearly all plants examined so far. Our study of GBSSI had three parts: (1) Amino acid sequences were compared across a broad taxonomic range, including grasses, four dicotyledons, and the microbial homologs of GBSSI. Inferred structural information was used to aid in the alignment of these very divergent sequences. The informed alignments highlight amino acids that are conserved across all sequences, and demonstrate that structural motifs can be highly conserved in spite of marked divergence in amino acid sequence. (2) Maximum-likelihood (ML) analyses were used to examine exon sequence evolution throughout grasses. Differences in probabilities among substitution types and marked among-site rate variation contributed to the observed pattern of variation. Of the parameters examined in our set of likelihood models, the inclusion of among-site rate variation following a gamma distribution caused the greatest improvement in likelihood score. (3) We performed cladistic parsimony analyses of GBSSI sequences throughout grasses, within tribes, and within genera to examine the phylogenetic utility of the gene. Introns provide useful information among very closely related species, but quickly become difficult to align among more divergent taxa. Exons are variable enough to provide extensive resolution within the family, but with low bootstrap support. The combined results of amino acid sequence comparisons, maximum-likelihood analyses, and phylogenetic studies underscore factors that might affect phylogenetic reconstruction. In this case, accommodation of the variable rate of evolution among sites might be the first step in maximizing the phylogenetic utility of GBSSI.      

The Genetic Structure of Populations of Potato virus Y in Japan; Based on the Analysis of 20 Full Genomic Sequences

The genetic structure of Potato virus Y (PVY) populations in Japan was analysed using 20 isolates; five were retrieved from the public DNA sequence databases, and an additional 15 complete genomic sequences were determined using field samples collected in Japan. Recombination and phylogenetic analyses of a total of 149 isolates from Japan and other countries showed that PVY has three major lineages (C, N and O); at least one, two and six sublineages in C, N and O lineages, respectively. One recombination pattern was newly found among Japanese PVY^NTN strain isolates, which was most closely related to the PVY^NTN strain isolates previously found in Europe and North America. On the other hand, PVY^O was a complex of several divergent lineages, and there were at least three non‐recombinant subpopulations in Japan. Studies on nucleotide diversities of populations and phylogenetic relationships of the isolates in the PVY sequences showed that Japanese PVY populations were in part distinct from the European and North American populations.     

Secondary Structure Analyses of the Nuclear rRNA Internal Transcribed Spacers and Assessment of Its Phylogenetic Utility across the Brassicaceae (Mustards)

The internal transcribed spacers of the nuclear ribosomal RNA gene cluster, termed ITS1 and ITS2, are the most frequently used nuclear markers for phylogenetic analyses across many eukaryotic groups including most plant families. The reasons for the popularity of these markers include: 1.) Ease of amplification due to high copy number of the gene clusters, 2.) Available cost-effective methods and highly conserved primers, 3.) Rapidly evolving markers (i.e. variable between closely related species), and 4.) The assumption (and/or treatment) that these sequences are non-functional, neutrally evolving phylogenetic markers. Here, our analyses of ITS1 and ITS2 for 50 species suggest that both sequences are instead under selective constraints to preserve proper secondary structure, likely to maintain complete self-splicing functions, and thus are not neutrally-evolving phylogenetic markers. Our results indicate the majority of sequence sites are co-evolving with other positions to form proper secondary structure, which has implications for phylogenetic inference. We also found that the lowest energy state and total number of possible alternate secondary structures are highly significantly different between ITS regions and random sequences with an identical overall length and Guanine-Cytosine (GC) content. Lastly, we review recent evidence highlighting some additional problematic issues with using these regions as the sole markers for phylogenetic studies, and thus strongly recommend additional markers and cost-effective approaches for future studies to estimate phylogenetic relationships.     

RNA-dependent RNA polymerase gene sequence from foot-and-mouth disease virus in Hong Kong

A foot-and-mouth disease virus (FMDV, HKN/2002) was isolated in Hong Kong in 2002. The nucleotide sequence of the 3D(pol) gene encoding the viral RNA-dependent RNA polymerase was determined and compared with that of the same gene from other FMDVs. The 3D(pol) gene was 1410 nucleotides in length encoding a protein of 470 amino acid residues. Sequence comparisons indicated that HKN/2002 belonged to serotype O. An evolutionary tree based on the 3D(pol) sequences of 20 FMDV isolates revealed that the nucleotide sequence of the HKN/2002 3D(pol) gene was most similar to those of isolates found in Taiwan in 1997, suggesting that they share a common ancestor. The amino acid sequence of the HKN/2002 3D(pol) gene was determined and aligned with those of representative isolates from seven other Picornaviridae genera. Eight highly conserved regions were detected, indicating a conserved functional relevance for these motifs. Alignment of 20 FMDV 3D(pol) amino acid sequences revealed a hypermutation region near the N-terminus that may help the virus evade host immune systems.     

Variation patterns of the mitochondrial 16S rRNA gene with secondary structure constraints and their application to phylogeny of cyprinine fishes (Teleostei: Cypriniformes)

The mitochondrial 16S ribosomal RNA (rRNA) gene sequences from 93 cyprinid fishes were examined to reconstruct the phylogenetic relationships within the diverse and economically important subfamily Cyprininae. Within the subfamily a biased nucleotide composition (A>T, C>G) was observed in the loop regions of the gene, and in stem regions apparent selective pressures of base pairing showed a bias in favor of G over C and T over A. The bias may be associated with transition-transversion bias. Rates of nucleotide substitution were lower in stems than in loops. Analysis of compensatory substitutions across these taxa demonstrates 68% covariation in the gene and a logical weighting factor to account for dependence in mutations for phylogenetic inference should be 0.66. Comparisons of varied stem-loop weighting schemes indicate that the down-weightings for stem regions could improve the phylogenetic analysis and the degree of non-independence of stem substitutions was not as important as expected. Bayesian inference under four models of nucleotide substitution indicated that likelihood-based phylogenetic analyses were more effective in improving the phylogenetic performance than was weighted parsimony analysis. In Bayesian analyses, the resolution of phylogenies under the 16-state models for paired regions, incorporating GTR + G + I models for unpaired regions was better than those under other models. The subfamily Cyprininae was resolved as a monophyletic group, as well as tribe Labein and several genera. However, the monophyly of the currently recognized tribes, such as Schizothoracin, Barbin, Cyprinion + Onychostoma lineages, and some genera was rejected. Furthermore, comparisons of the parsimony and Bayesian analyses and results of variable length bootstrap analysis indicates that the mitochondrial 16S rRNA gene should contain important character variation to recover well-supported phylogeny of cyprinid taxa whose divergences occurred within the recent 8 MY, but could not provide resolution power for deep phylogenies spanning 10-19 MYA.     

Angiosperm phylogeny based on matK sequence information

Plastid matK gene sequences for 374 genera representing all angiosperm orders and 12 genera of gymnosperms were analyzed using parsimony (MP) and Bayesian inference (BI) approaches. Traditionally, slowly evolving genomic regions have been preferred for deep-level phylogenetic inference in angiosperms. The matK gene evolves approximately three times faster than the widely used plastid genes rbcL and atpB. The MP and BI trees are highly congruent. The robustness of the strict consensus tree supercedes all individual gene analyses and is comparable only to multigene-based phylogenies. Of the 385 nodes resolved, 79% are supported by high jackknife values, averaging 88%. Amborella is sister to the remaining angiosperms, followed by a grade of Nymphaeaceae and Austrobaileyales. Bayesian inference resolves Amborella + Nymphaeaceae as sister to the rest, but with weak (0.42) posterior probability. The MP analysis shows a trichotomy sister to the Austrobaileyales representing eumagnoliids, monocots + Chloranthales, and Ceratophyllum + eudicots. The matK gene produces the highest internal support yet for basal eudicots and, within core eudicots, resolves a crown group comprising Berberidopsidaceae/Aextoxicaceae, Santalales, and Caryophyllales + asterids. Moreover, matK sequences provide good resolution within many angiosperm orders. Combined analyses of matK and other rapidly evolving DNA regions with available multigene data sets have strong potential to enhance resolution and internal support in deep level angiosperm phylogenetics and provide additional insights into angiosperm evolution.