哺乳动物FoxL2基因的分子进化分析

本研究以FoxL2基因为研究对象,选取Gen Bank公布的26个物种的FoxL2基因氨基酸序列,代表脊椎动物中5纲15目21科24属,重点分析哺乳动物中该基因的结构与功能相关性以及其系统发育情况。结果表明:26条FoxL2基因序列存在非常保守的forkhead结构域,哺乳类比非哺乳动物多出丙氨酸、甘氨酸及脯氨酸重复区,其中多聚丙氨酸束的缩短可能与Ⅰ型BPES相关,而多聚丙氨酸束的延伸使核内蛋白错位和细胞质聚集导致Ⅱ型BPES。FoxL2含有较多磷酸化位点,推测其在体内的活性可能是通过磷酸化调控的。FoxL2基因氨基酸序列构建的系统发育树中亲缘关系较近的物种能够较好的聚在一起。     

冬虫夏草菌3个细胞核蛋白编码基因的分子进化

【目的】分析3个细胞核蛋白编码基因(csp1、MAT1-1-1和MAT1-2-1)在不同冬虫夏草菌株间的分子进化。【方法】从125个冬虫夏草样品中分别扩增csp1、MAT1-1-1和MAT1-2-1基因序列,比较外显子和内含子间以及2个交配型基因间的序列变异程度,比较基于不同基因或基因区域所构建的系统发育树拓扑结构的差异,分析3个基因承受的选择压力和DNA重组情况。【结果】3个蛋白编码基因外显子区的长度在不同菌株间高度保守,具有4.5%?5.7%的变异位点;内含子区的长度在不同菌株间相同或不同,具有1.8%?22%的变异位点。对于2个交配型基因,MAT1-1-1的碱基变异率低于MAT1-2-1。基于外显子与内含子构建的系统发育树的拓扑结构,以及基于2个交配型基因外显子构建的系统发育树的拓扑结构都存在明显差异。3个蛋白编码基因都经历着净化选择作用。基因内部的不同DNA位点间有重组,但3个基因片段之间没有明显的重组发生。【结论】由于冬虫夏草菌不同基因以及基因的不同区域表现出进化上的差异,所以在开展冬虫夏草菌进化相关的研究时,应该联合使用多个不同的基因片段。     

Unc5基因家族多样性进化的研究

Uncoordinated-5(Unc5)基因家族属于经典的轴突导向基因家族。为了探讨Unc5的进化和分化规律,首先通过序列比对和蛋白质结构预测鉴定了Unc5基因家族成员的起源和分布,再利用PAML和DIVERGE软件分析了各基因亚型的进化选择压力和功能歧化。结果表明,Unc5基因家族在脊椎动物中受到了不同程度的选择压力,其中Unc5C基因型受到了纯化选择作用,而Unc5A、Unc5B和Unc5D基因型受到了正选择作用,并且在这3个基因型中分别检测到了8个、1个和6个正选择位点。此外,DIVERGE软件检测出38个I型功能歧化位点和97个Ⅱ型功能歧化位点。这些正选择位点和功能分歧位点为进一步研究Unc5蛋白的结构和功能提供了参考和依据。     

念珠藻属植物hetR基因的适应性进化分析

以念珠藻属(Nostoc)及其近缘类群hetR基因的51条序列为研究对象,对hetR基因的编码蛋白进行生物信息学分析和系统发育分析,并使用分支模型、位点模型和分支-位点模型进行该基因位点的适应性进化研究。系统发育分析结果显示,51条hetR基因蛋白序列可分为4个大分支。适应性进化分析结果表明,在3种进化模型中,大多数分支及藻株都没有检测到统计学上具有显著性的正选择位点,说明检测的位点大多处于负选择压力下。但在普通念珠藻(Nostoc commune,CHAB2802)中检测到正选择位点(126T),提示念珠藻属植物hetR基因发生了适应性改变。     

拟黑多刺蚁脱皮激素受体编码基因的克隆及系统发育分析

为了研究脱皮激素受体在拟黑多刺蚁发育中的功能,本项目采用逆转录PCR方法从拟黑多刺蚁(Polyrhachis vicina Roger)中克隆到脱皮激素受体编码基因PvEcR全长序列,对其编码的氨基酸序列进行生物信息学分析。从拟黑多刺蚁中克隆到脱皮激素受体蛋白编码基因1 737bp的全长序列,该序列编码578个氨基酸,预测的蛋白分子量大小为63.098×103,理论等电点为7.41。NCBI蛋白质数据库中进行Blast搜索表明,拟黑多刺蚁脱皮激素受体蛋白PvEcR存在至少6类保守结构区域,主要为DNA结合位点、配体结合位点和激活因子识别位点。PvEcR蛋白序列磷酸化位点预测共揭示48个可能的磷酸化位点。PvEcR与其他昆虫脱皮激素受体蛋白在氨基酸序列上同源性高达70%以上。基于脱皮激素受体蛋白氨基酸序列构建的系统发育树表明,几乎所有蚁类脱皮激素受体形成一个大分支,其中拟黑多刺蚁与佛罗里达弓背蚁表现出最近的亲缘关系。     

绵羊TFAM基因编码蛋白生物信息学分析

为进一步了解和研究绵羊线粒体转录因子A(TFAM)基因的结构和功能,采用生物信息学方法构建了不同物种的TFAM基因系统发育树,并对绵羊TFAM基因编码蛋白质理化性质、二级结构、亚细胞定位、蛋白磷酸化、三级结构及功能结构域等进行分析。结果表明:绵羊和山羊亲缘关系最为接近,其次是家牛,与原鸡、斑马鱼亲缘关系较远。绵羊TFAM基因共编码246个氨基酸,以赖氨酸含量最高;二级结构以α-螺旋和无规则卷曲占为主;综合分析,此蛋白为碱性、不稳定性亲水蛋白,不存在信号肽和跨膜结构域;其主要在细胞核、线粒体、细胞骨架内发挥其生物学作用。该蛋白的丝氨酸磷酸化位点比较多,而酪氨酸和苏氨酸的磷酸化位点相对较少;且含有2个HMG结构域,在两个HMG结构域之间存在一个由12个氨基酸残基组成的短连接区域。     

桑黄Zn(II)2Cys6锌簇蛋白转录因子及其在不同碳氮源条件下的表达

Zn(II)2Cys6锌簇蛋白转录因子（C6 zinc）在真菌次生代谢产物合成中发挥重要作用。本研究首先利用本地BLAST,从桑黄Sanghuangporus sanghuang全基因组中获得Zn(II)2Cys6锌簇蛋白转录因子编码基因,并利用生物信息学手段分析编码蛋白保守结构域、一级结构及二级结构,构建蛋白系统发育树,最后利用半定量PCR对它们在不同碳源和氮源培养条件下的表达情况进行检测。结果显示:从桑黄基因组中分析获得的11个Zn(II)2Cys6锌簇蛋白均具有Cy6型锌指基序,属于GAL4型锌簇蛋白转录因子;它们均为不稳定亲水性蛋白,具磷酸化修饰位点,糖基化修饰位点较少或没有,定位于细胞核中;其二级结构主要以无规卷曲和α螺旋为主;系统发育树分析结果显示,桑黄Zn(II)2Cys6锌簇蛋白分为2个大分支,其中Ⅰ类分支转录因子的保守结构域分布于蛋白N-端,Ⅱ类分支转录因子的保守结构域则分布于蛋白C-端;11个桑黄Zn(II)2Cys6锌簇蛋白转录因子在不同培养基培养的菌丝体中呈现差异化表达,其中,肌醇培养基和乳糖培养基能够有效促进大部分锌簇蛋白转录因子的表达;另外,基因簇分析显示桑黄锌簇蛋白SHCZ4可能是NRPS-PKS杂合基因簇体系的途经特异性转录因子。该结果将为桑黄次生代谢产物合成调控相关转录因子的研究以及潜在次生代谢相关基因簇的挖掘提供参考依据。     

40种蕨类植物matK基因的系统分类及分子进化研究

采用“放松分子钟”模型、氨基酸位点正选择模型和分子内共进化网络估算方法,对蕨类植物Ⅱ型内含子成熟酶蛋白K(Maturase K,MATK)编码基因matK的进化趋势进行研究。结果显示:matK基因在蕨类植物系统学研究中具有一定的应用价值,与rbcL基因和psaA基因联合后能显著提升系统发育树的可信度;蕨类植物MATK蛋白中存在少数曾经历正选择的位点;MATK蛋白内部有多对氨基酸位点共同构成共进化网络。在被子植物兴起环境改变后,MATK蛋白部分位点发生适应性进化,通过位点间共进化网络协同作用方式提升蕨类植物对新光合环境的适应能力。     

格特隐球菌核基因、交配型位点和线粒体基因的多位点序列分型及重组分析

研究格特隐球菌VGI和VGII基因型菌株间线粒体基因重组与核基因、交配型位点重组间的关系。采用分子鉴定区分受试格特隐球菌的血清型、基因型和交配型;选取包含核基因和线粒体基因10个位点的多位点序列分型（ MLST）进行基因型和系统发育分析;采用同质性检验评价基因系谱间的一致性。多位点的序列和系统发育分析结合同质性检验表明,受试位点中,仅ATP6位点发现VGI和VGII菌株间基因的杂交和重组,该基因系谱与其余位点的基因系谱呈现不一致性。结果显示,受试的部分VGI菌株中的ATP6位点含有VGII菌株的基因序列,表明VGI和VGII菌株间线粒体基因的重组;没有发现核基因及交配型位点中两者间的重组,提示VGI和VGII菌株间线粒体基因的重组现象与交配行为无关。     

弯枝藻属rbcL基因的适应性进化分析

为探讨淡水红藻的叶绿体基因及其适应性进化特征,选取弯枝藻属(Compsopogon)及相近外类群的rbc L基因共17条,利用PAML 4.9软件,对弯枝藻属rbc L基因编码蛋白进行生物信息学分析,并分别采用分支模型、位点模型以及分支-位点模型对基因的选择位点进行检测。结果表明,弯枝藻属rbc L基因编码蛋白的二级结构主要由α螺旋和β折叠构成,结构稳定。采用最大似然法构建的系统发育树表明,内类群为单一物种,分为3个小分支,具有一定地理分布规律。在3种进化模型中均未检测到统计上显著的正选择位点,表明绝大多数位点处于负选择压力下。因此,弯枝藻属rbc L基因未发生适应性进化。     

Significantly different patterns of amino acid replacement after gene duplication as compared to after speciation

We have performed a large-scale analysis of amino acid sequence evolution after gene duplication by comparing evolution after gene duplication with evolution after speciation in over 1,800 phylogenetic trees constructed from manually curated alignments of protein domains downloaded from the PFAM database. The site-specific rate of evolution is significantly altered by gene duplication. A significant increase in the proportion of amino acid substitutions at constrained (slowly evolving) sites after duplication was observed. An increase in the proportion of replacements at normally constrained amino acid sites could result from relaxation of purifying selective pressure. However, the proportion of amino acid replacements involving radical changes in amino acid properties after duplication does not appear to be significantly increased by relaxed selective pressure. The increased proportion of replacements at constrained sites was observed over a relatively large range of protein change (up to 25% amino acid replacements per site). These findings have implications for our understanding of the nature of evolution after duplication and may help to shed light on the evolution of novel protein functions through gene duplication.     

新城疫病毒分离株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基因的变异与分离时间有一定的联系.     

The evolution of lysozyme and alpha-lactalbumin

From the analysis of phylogenetic trees constructed from the amino acid sequences and metal-binding properties of various lysozymes c and alpha-lactalbumins, it was found that before the divergence of the lineages of birds and mammals, calcium-binding lysozyme diverged from non-calcium-binding lysozyme. alpha-Lactalbumin evolved from the calcium-binding lysozyme along the mammalian lineage after the divergence of birds and mammals. Rapid evolution took place, not in the process of acquisition of the activity of alpha-lactalbumin, but after the loss of lysozyme activity, due to the change in the distribution of selective pressure on each amino acid site. A general process for the change in function of a protein during evolution is suggested to be as follows: after duplication of the gene, one of their protein products acquires a new function, besides that already present; the old function is eventually lost.     

Evidence of positive selection at codon sites localized in the C-terminal peptide of ORC6

Origin recognition complex 6 (Orc6) plays a central role in the initiation of DNA replication in all eukaryotic systems. The exact contribution of Orc6 to replication initiation has yet to be elucidated. Here, we analyzed the evolutionary dynamics of Orc6 in 15 vertebrates. Positive selection was detected in the region of exon 6 of the Orc6 gene. Site tests revealed a proportion of codon sites that displayed evidence of positive selection (ω > 1) within the coding sequences of the vertebrate Orc6 gene. Seven positively selected amino acid sites were identified and three were located in exon6. These results suggest that amino acid residues present in the middle region of the protein are more selectively constrained, whereas amino acid residues in the C-terminal peptide of the protein evolve at a faster rate, possibly because of heightened selective pressure during the course of evolution.     

Structure of a cDNA for the pro alpha 2 chain of human type I procollagen. Comparison with chick cDNA for pro alpha 2(I) identifies structurally conserved features of the protein and the gene

Nucleotide sequences were determined for cloned cDNAs encoding for more than half of the pro alpha 2 chain of type I procollagen from man. Comparisons with previously published data on homologous cDNAs from chick embryos made it possible to examine evolution of the gene in two species which have diverged for 250-300 million years. The amino acid sequence of the alpha-chain domain supported previous indications that there is a strong selective pressure to maintain glycine as every third amino acid and to maintain a prescribed distribution of charged amino acids. However, there is little apparent selective pressure on other amino acids. The amino acid sequence of the C-propeptide domain showed less divergence than the alpha-chain domain. The 5' end or N terminus of the human C-propeptide, however, contained an insert of 12 bases coding for 4 amino acids not found in the chick C-propeptide. About 100 amino acid residues from the N terminus, two residues found in the chick sequence were missing from the human. In the second half of the C-propeptide, there was complete conservation of a 37 amino acid sequence and conservation of 50 out of 51 amino acids in the same region, an observation which suggested that the region serves some special purpose such as directing the association of one pro alpha 2(I) C-propeptide with two pro alpha 1(I) C-propeptides so as to produce the heteropolymeric structure of type I procollagen. In addition, comparison of human and chick DNAs for pro alpha 2(I) revealed three different classes of conservation of nucleotide sequence which have no apparent effect on the structure of the protein: a preference for U on the third base position of codons for glycine, proline, and alanine; a high degree of nucleotide conservation in the 51 amino acid highly conserved region of the C-propeptide; a high degree of nucleotide conservation in the 3'-noncoding region. These three classes of nucleotide conservation may reflect unusual features of collagen genes, such as their high GC content or their highly repetitive coding sequences.     

BRCA1蛋白的序列及空间结构的分析

查询人的BRCA1蛋白的氨基酸序列,利用生物信息学的方法进行相似性搜索,获得一系列BRCA1蛋白的氨基酸序列。选择了其中的11条序列,对BRCA1蛋白进行了多重序列分析和进化分析,对BRCA1蛋白的BRCT结构域进行三维同源模型的构建与比较分析。分析结果表明:BRCA1中某些特定部位的氨基酸序列高度保守;确定氨基酸的保守位点并联合进化分析可对基因错义突变的致病性做初步地猜测;相近物种来源的BRCA1具有较近的亲缘关系,而且具有极其相似的三维空间结构。这些为研究BRCA1蛋白的结构与功能关系提供指导意义。     

Molecular evolution of the dotA gene in Legionella pneumophila

The molecular evolution of dotA, which is related to the virulence of Legionella pneumophila, was investigated by comparing the sequences of 15 reference strains (serogroups 1 to 15). It was found that dotA has a complex mosaic structure. The whole dotA gene of Legionella pneumophila subsp. pneumophila serogroups 2, 6, and 12 has been transferred from Legionella pneumophila subsp. fraseri. A discrepancy was found between the trees inferred from the nucleotide and deduced amino acid sequences of dotA, which suggests that multiple hits, resulting in synonymous substitutions, have occurred. Gene phylogenies inferred from three different segments (the 5'-end region, the central, large periplasmic domain, and the 3'-end region) showed impressively dissimilar topologies. This was concordant with the sequence polymorphisms, indicating that each region has experienced an independent evolutionary history, and was evident even within the same domain of each strain. For example, the PP2 domain was found to have a heterogeneous structure, which led us hypothesize that the dotA gene of L. pneumophila may have originated from two or more different sources. Comparisons of synonymous and nonsynonymous substitutions demonstrated that the PP2 domain has been under strong selective pressure with respect to amino acid change. Split decomposition analysis also supported the intragenic recombination of dotA. Multiple recombinational exchange within the dotA gene, encoding an integral cytoplasmic membrane protein that is secreted, probably provided increased fitness in certain environmental niches, such as within a particular biofilm community or species of amoebae.     

Genealogical evidence for positive selection in the nef gene of HIV-1.

The pattern and process of evolution in the nef gene of HIV-1 was analyzed within and among patients. Using a maximum likelihood method that allows for variable intensity of selection pressure among codons, strong positive selection was detected in a hemophiliac patient over 30 mo of infection. By reconstructing the process of allele substitution in this patient using parsimony, the synapomorphic amino acid changes separating each time point were found to have high probabilities of being under positive selection, with selective coefficients of at least 3.6%. Positive selection was also detected among 39 nef sequences from HIV-1 subtype B. In contrast, multiple pairwise comparisons of nonsynonymous and synonymous substitution rates provided no good evidence for positive selection and sliding window analyses failed to detect most positively selected sites. These findings demonstrate that positive selection is an important determinant of nef gene evolution and that genealogy-based methods outperform pairwise methods in the detection of adaptive evolution. Mapping the locations of positively selected sites may also be of use in identifying targets of the immune response and hence aid vaccine design.