进化距离与相对替代率的比较研究

系统树是描述物种、人种甚至基因间亲缘关系和演化的重要工具,必须以进化距离或(相对)替代率作为重要的参数。但以哪一个参数构建的树更能反映真实的系统树呢?事实上迄今并无人对此作过认真的研究。本文以模拟数据并用方差分析法检验两个参数的异同并讨论其包含的生物学意义。研究结果表明,当氨基酸的替代率和核苷酸的替代率分别为0.18和0.13时,它们的进化距离分别为0.199和0.143。经方差分析证实,若检验的氨基酸和核苷酸最大数目均为75只时,不论以替代率或进化距离中那一个作为构树参数,拓扑树事实上几乎只有一个。这就是说,该时拓扑树可靠性很大,而且随着它们替代率的减少,则检验的氨基酸和核苷酸的数目会随之增加。但是一旦氨基酸的替代率和核苷酸的替代率超过上述数字,则两个参数构建的树在拓扑长度上是不等价的。经分析,若进化距离大致上与进化时间成线性关系的话,则应选用进化距离。用进化距离重建的系统树事实上支持中性学说;若进化距离与进化时间显著地不存在线性关系的话,则可选用替代率,该情况表明中性学说不适用,似更倾向于新达尔文主义。     

DNA序列进化过程中核苷酸替代的非独立性研究

本文评述了DNA序列间核苷酸替代数的估计方法,并通过对七个物种中组蛋白基因的比较对DNA进化的模型进行了考察。发现H2A基因第三位点上的碱基组成在物种间变异很大,并且跟H2A基因第一位点、H4基因第一、三位点及H2A上游,下游序列中的碱基组成有强正相关,提示DNA序列进化过程中存在着物种特异的区域性约束力。可能的原因是高等真核生物中GC含量升高,或者是染色体重组使这些同源序列位于不同的等质区段,从而受到不同的选择突变压。密码内各位点上核苷酸替代的相关性分析表明不同位点的替代是非独立的,其原因可能是一次替代事件引起多个位点的变化。文中讨论了这些结果对进化树推断的意义。     

DNA序列基于k-字的数值刻画及其应用

借助DNA序列中k-字的频数,将序列转化成一个340维向量,进而计算物种间的进化距离。作为应用:分别以15个物种的β球蛋白基因、13种汉坦病毒的S片段以及26个闭壳龟线粒体基因为例,构建系统发生树,所得结果与前人的结论一致,说明了该方法的有效性。     

真核生物DNA非编码区的组分分析

在全基因组水平上,用直方图、混沌表示灰度图、距离差异度和信息熵差异度四种方法,研究了拟南芥、线虫、果蝇的DNA内含子、基因间隔区DNA、外显子三种区域的核苷酸短序列组分及组分复杂度.结果表明：a.不同基因组之间,不管基因数目多少,用4种方法得到的外显子部分其组分复杂度都比较接近,而非编码区部分的组分复杂度却很大.这一点定量地说明了物种之间的复杂程度,主要不体现在编码区部分,而体现在非编码区部分.b.同一基因组中,内含子的核苷酸短序列组分复杂度都是相似的,外显子和intergenic DNA部分的组分复杂度也是相似的.c.内含子和intergenic DNA在转录、剪切、二级结构等方面有很大的不同,但它们在核苷酸短序列组分上的差异却很小,说明内含子和intergenic DNA在转录、剪切、二级结构上的不同并不通过核苷酸短序列组分来进行限制.     

高效液相色谱法在脱氧核糖核酸酶解动力学

为了更好地对脱氧核糖核酸酶解动力学过程进行研究,建立脱氧核糖核酸(DNA)酶解液中4种脱氧核苷酸(腺嘌呤脱氧核苷酸(dAMP)、鸟嘌呤脱氧核苷酸(dGMP)、胞嘧啶脱氧核苷酸(dCMP)、胸腺嘧啶脱氧核苷酸(dTMP))的高效液相测定方法,能将酶解液中4种脱氧核苷酸完全分离并准确定量.在此基础上,对DNA酶解的动力学进行初步研究,其反应机理为不存在底物和产物抑制的双底物顺序反应,动力学方程为x=1/bln(1+abt)(其中a=0.372 3ρ0-0.974;b=-0.049 3ρ20+1.115 3ρ0 - 1.110 3),该方程可以很好地描述DNA酶解过程,误差仅为3.31%.     

牛科动物HSL基因序列分析及其分子进化研究

在对牛科中4种动物即牦牛、瘤牛、普通牛和水牛HSL基因外显子Ⅰ部分核苷酸序列进行测定的基础上,与Gen-Bank中其他物种相应基因核苷酸序列、氨基酸序列进行了比对分析,并构建了牦牛与其他物种间分子系统进化树。结果表明:牦牛与普通牛、瘤牛、水牛、猪、人、小鼠、大鼠7个物种HSL基因外显子Ⅰ部分核苷酸序列间保守性较高,同源性大小依次为99.8%、99.6%、97.4%、90.6%、88.4%、83.5%、82.3%。相应氨基酸序列间保守性更高,同源性分别为100%、100%、98.2%、94.0%、92.2%、89.8%、89.8%。牦牛与各物种该基因部分核苷酸序列间碱基变异类型主要表现为碱基转换和颠换,无碱基插入和缺失发生,碱基转换的频率高于颠换的频率;在核苷酸水平上的多数碱基替换都是同义替换;序列间单碱基变异位点大多出现在同一位点,多发生在密码子第3位,其次是第1位,最少发生在第2位,符合分子进化的中性学说。HSL基因外显子Ⅰ部分核苷酸序列进行多序列对位排列构建的各物种间分子系统进化树结果表明,普通牛和瘤牛首先聚为一类,再分别与牦牛、水牛、猪、人聚类,最后与大鼠、小鼠聚为一类。该聚类结果与动物学上的分类结果一致,表明HSL基因外显子Ⅰ部分核苷酸序列适合于构建物种间分子系统进化树。研究表明,牦牛、普通牛和瘤牛3个物种间的遗传距离大小相近,牦牛和水牛间的遗传距离与普通牛、瘤牛和水牛间的遗传距离大小相当。牦牛、普通牛和瘤牛3个物种间的遗传距离远小于它们各自与水牛这一物种的遗传距离,它们三者之间的亲缘关系也相对于它们各自与水牛间的亲缘关系都较近,故将牦牛、普通牛和瘤牛划分在同一个属——牛属(Bos)更为合理。     

青海省6种裂腹鱼COI基因作为DNA条形码的适用性初探

青海省裂腹鱼鱼类至少有20种,占土著鱼类的40%以上,具有重要的生态价值。由于生态环境的恶化和人为因素的干扰,很多物种群已濒临灭绝。快速和准确的物种鉴定对于这些物种的保护至关重要,而基于形态学的传统分类法很难满足这一需求。因此,本研究通过DNA条形码技术初步探讨了在青海省裂腹鱼物种鉴定中的适用性。本研究中,测序获得了青海26尾裂腹鱼(6个物种)的细胞色素c氧化酶亚基I (COI)基因,并经GenBank数据库进行比对;基于K2P模型分析COI序列变异;运用贝叶斯法(BI)和最大似然法(ML)构建系统发育树。结果显示,6种裂腹鱼COI序列检测得到15个单倍型,且各物种间无共享的单倍型。基于K2P模型,最大种内遗传距离和最小种间遗传距离分别为(0.621±0.297)%和(2.792±0.644)%。种间平均遗传距离为(12.205±1.307)%,约为种内平均遗传距离((0.327±0.162)%)的37倍,表明各物种的COI序列间已经形成明显"条形码间隙"。AMOVA分析显示,遗传变异主要来自种间,约占97.05%;FST=0.970 54,p<0.01,说明各物种间分化程度极高。此外,基于BI和ML方法构建的系统树具有一致的拓扑结构,分辨率较高;各物种均单独聚为一个发育枝,拓扑结构合理。以上结果表明,COI基因作为DNA条形码在青海裂腹鱼物种鉴定中具有较高的适用性。     

青海省6种裂腹鱼COI基因作为DNA条形码的适用性初探

青海省裂腹鱼鱼类至少有20种,占土著鱼类的40%以上,具有重要的生态价值。由于生态环境的恶化和人为因素的干扰,很多物种群已濒临灭绝。快速和准确的物种鉴定对于这些物种的保护至关重要,而基于形态学的传统分类法很难满足这一需求。因此,本研究通过DNA条形码技术初步探讨了在青海省裂腹鱼物种鉴定中的适用性。本研究中,测序获得了青海26尾裂腹鱼(6个物种)的细胞色素c氧化酶亚基I (COI)基因,并经GenBank数据库进行比对;基于K2P模型分析COI序列变异;运用贝叶斯法(BI)和最大似然法(ML)构建系统发育树。结果显示,6种裂腹鱼COI序列检测得到15个单倍型,且各物种间无共享的单倍型。基于K2P模型,最大种内遗传距离和最小种间遗传距离分别为(0.621±0.297)%和(2.792±0.644)%。种间平均遗传距离为(12.205±1.307)%,约为种内平均遗传距离((0.327±0.162)%)的37倍,表明各物种的COI序列间已经形成明显"条形码间隙"。AMOVA分析显示,遗传变异主要来自种间,约占97.05%;FST=0.970 54,p0.01,说明各物种间分化程度极高。此外,基于BI和ML方法构建的系统树具有一致的拓扑结构,分辨率较高;各物种均单独聚为一个发育枝,拓扑结构合理。以上结果表明,COI基因作为DNA条形码在青海裂腹鱼物种鉴定中具有较高的适用性。     

嗡蜣螂属Onthophagus十二种蜣螂线粒体COI基因的DNA条形码研究

对嗡蜣螂属Onthophagus 12种蜣螂的线粒体COI基因3’端部分序列(731 bp)进行了比较,结果显示,COI序列的变异位点213个,简约信息位点167个.碱基替代主要发生在第3位点(64次),占替代总数的83.12％.除掘嗡蜣螂O.fodiens与婪嗡蜣螂O.lenzi小于2％外,其余种间遗传距离在8.1％ ～15.8％之间,种内遗传距离为0 ～0.2％.单倍型多样性(Hd)和核苷酸序列多样性(Pi)分别为0.944±0.030和0.10518±0.0045.滑动窗口分析表明,可变位点频率在240～290 bp、675 bp附近较高.NJ树聚类结果与传统形态学分类相吻合:外群代表种分化最早,种间聚成一分支,种内个体优先聚集种下.本文认为COI基因适合作为嗡蜣螂属物种鉴定的DNA条形码.     

DNA动力学柔性的统计力学模型

考虑碱基对之间的非紧邻相互作用、涨落的序列依赖效应和非对称涨落,提出了DNA构象的统计力学模型,给出了DNA柔性的新定义。作为模型的应用,对12种三核苷酸重复序列的动力学柔性作了预测。理论预测与其它方法得到的结论比较,有很好的一致性。对模型和结论的理论意义作了讨论。     

Manganese-histidine cluster as the functional center of the water oxidation complex in photosynthesis

The recent model of Kambara and Govindjee for water oxidation [Kambara T. and Govindjee (1985) Proc. Natl. Acad. Sci. U.S.A., 82:6119–6123] has been extended in this paper by examining all the data in order to identify the most likely candidate for the redox-active ligand (RAL), suggested to operate between the water oxidizing complex (WOC) and Z, the electron donor to the reaction center P680. We have concluded that a very suitable candidate for RAL is the imidazole moiety of a histidine residue. The electrochemical data available on imidazole derivatives play heavily in this identification of RAL. Thus, we suggest that histidine might play the role of an electron mediator between the WOC and Z. A model of S-states in terms of their plausible chemical identity is presented here.Abbreviations J electronic spin of ion - P680 reaction center chlorophyll - RAL Redox active ligand - S_n state of the oxygen-evolving system - WOC water oxidation complex - Z electron donor to P680 Dedicated to Prof. L.N.M. Duysens on the occasion of his retirement     

Genetic biasing through cultural transmission: Do simple Bayesian models of language evolution generalise?

The recent Bayesian approaches to language evolution and change seem to suggest that genetic biases can impact on the characteristics of language, but, at the same time, that its cultural transmission can partially free it from these same genetic constraints. One of the current debates centres on the striking differences between sampling and a posteriori maximising Bayesian learners, with the first converging on the prior bias while the latter allows a certain freedom to language evolution. The present paper shows that this difference disappears if populations more complex than a single teacher and a single learner are considered, with the resulting behaviours more similar to the sampler. This suggests that generalisations based on the language produced by Bayesian agents in such homogeneous single agent chains are not warranted. It is not clear which of the assumptions in such models are responsible, but these findings seem to support the rising concerns on the validity of the “acquisitionist” assumption, whereby the locus of language change and evolution is taken to be the first language acquirers (children) as opposed to the competent language users (the adults).     

A codon-based model of host-specific selection in parasites,with an application to the influenza A virus

Parasites sometimes expand their host range by acquiring a new host species. After a host change event, the selective regime acting on a given parasite gene may change as a result of host-specific adaptive alterations of protein functionality or host-specific immune-mediated selection. We present a codon-based model that attempts to include these effects by allowing the position-specific substitution process to change in conjunction with a host change event. Following maximum-likelihood parameter estimation, we employ an empirical Bayesian procedure to identify candidate sites potentially involved in host-specific adaptation. We discuss the applicability of the model to the more general problem of ascertaining whether the selective regime differs in two groups of related organisms. The utility of the model is illustrated on a data set of nucleoprotein sequences from the influenza A virus obtained from avian and human hosts.     

Evolve III: A discrete events model of an evolutionary ecosystem

Evolve III is a discrete events model of an evolutionary ecosystem. The model includes three levels of organization: population, organism and genetic structure. Each of these components was modeled independently, so that selective replacement of subsystems can be used to create families of models capable of testing alternative hypotheses about the real system. To demonstrate the use of the model we describe an experiment on the relationship between adaptability of populations and the variability of the environment. Populations cultured in a constant environment usually dominated those cultured in a variable environment when both were placed in a variable environment at an early stage of development, whereas the opposite is the case at later stages of development. This agrees with experiments on laboratory microcosms and lends credence to the potential predictive value of the model.     

Protein evolution: The deciphering of latent facets — correlation of synthesis profiles of ribosomally directed proteins and enzyme directed peptides

Key facets pertaining to the evolution of proteins have been probed, using as springboard, the relevant data bases constructed from (i) 60 ribosomally directed proteins, whose 3D structures are known and having 10,000 residues and (ii) from 73 enzyme directed peptides, comprising of 524 residues. The preference profiles, both in terms of the choice of neighbours and the placement of the peptide bonds, have been delineated with respect to each of the 20 coded amino acids. By and large, the preference profile from both the sets are similar, thus giving importance to the nature of the side chains of the coded amino acids. The predictive power of the preference profile has been tested with good results, thus demonstrating the evidence of common preference pathways for peptide formation during evolution. The ribosomally directed protein synthesis, controlled by the genome, proceeds by the addition of single residues at a time. On the other hand, the enzyme directed peptide synthesis largely operates in a more energy conscious block mode, where each constituent of a large enzyme ensemble is engaged in precisely assembling the modules and transfering them to the adjacent one, thus realizing a sequence specific peptide synthesis. Of significance is the fact that, in spite of such divergence in assembly, predictions for neighbour preferences in ribosomally directed protein synthesis work well when applied to enzyme directed peptide synthesis. The findings here are significant since they provide (i) a clear picture of directed peptide synthesis in the absence of direct genomic control, (ii) evidence for the preferred formation of peptide bonds using protein templates, (iii) they also provide evidence for the presence of protein like structures, with catalytic activity, prior to the freezing of the genetic code arising from dominance of the information system and (iv) a logical approach to the evolution of a hierarchical pattern. Presented at the National Symposium on Evolution of Lite.     

Compensatory evolution in diploid populations

Compensatory mutations are individually deleterious but harmless in appropriate combinations either at more than two sites within a gene or on separate genes. Considering that dominance effects of selection and heterodimer formation of gene products may affect the rate of compensatory evolution, we investigate compensatory neutral mutation models for diploid populations. Our theoretical analysis on the average time until fixation of compensatory mutations shows that these factors play an important role in reducing the fixation time of compensatory mutations if mutation rates are not low. Compensatory evolution of heterodimers is shown to occur more easily if the deleterious effects of single mutants are recessive.     

Equilibrium properties of a multi-locus, haploid-selection, symmetric-viability model

Under haploid selection, a multi-locus, diallelic, two-niche Levene (1953) model is studied. Viability coefficients with symmetrically opposing directional selection in each niche are assumed, and with a further simplification that the most and least favored haplotype in each niche shares no alleles in common, and that the selection coefficients monotonically increase or decrease with the number of alleles shared. This model always admits a fully polymorphic symmetric equilibrium, which may or may not be stable.We show that a stable symmetric equilibrium can become unstable via either a supercritical or subcritical pitchfork bifurcation. In the supercritical bifurcation, the symmetric equilibrium bifurcates to a pair of stable fully polymorphic asymmetric equilibria; in the subcritical bifurcation, the symmetric equilibrium bifurcates to a pair of unstable fully polymorphic asymmetric equilibria, which then connect to either another pair of stable fully polymorphic asymmetric equilibria through saddle-node bifurcations, or to a pair of monomorphic equilibria through transcritical bifurcations. As many as three fully polymorphic stable equilibria can coexist, and jump bifurcations can occur between these equilibria when model parameters are varied.In our Levene model, increasing recombination can act to either increase or decrease the genetic diversity of a population. By generating more hybrid offspring from the mating of purebreds, recombination can act to increase genetic diversity provided the symmetric equilibrium remains stable. But by destabilizing the symmetric equilibrium, recombination can ultimately act to decrease genetic diversity.     

Heterogeneous Nature and Distribution of Interruptions in Dinucleotides May Indicate the Existence of Biased Substitutions Underlying Microsatellite Evolution

Some aspects of microsatellite evolution, such as the role of base substitutions, are far from being fully understood. To examine the significance of base substitutions underlying the evolution of microsatellites we explored the nature and the distribution of interruptions in dinucleotide repeats from the human genome. The frequencies that we inferred in the repetitive sequences were statistically different from the frequencies observed in other noncoding sequences. Additionally, we detected that the interruptions tended to be towards the ends of the microsatellites and 5'-3' asymmetry. In all the estimates nucleotides forming the same repetitive motif seem to be affected by different base substitution rates in AC and AG. This tendency itself could generate patterning and similarity in flanking sequences and reconcile these phenomena with the high mutation rate found in flanking sequences without invoking convergent evolution. Nevertheless, our data suggest that there is a regional bias in the substitution pattern of microsatellites. The accumulation of random substitutions alone cannot explain the heterogeneity and the asymmetry of interruptions found in this study or the relative frequency of different compound microsatellites in the human genome. Therefore, we cannot rule out the possibility of a mutational bias leading to convergent or parallel evolution in flanking sequences.     

MADS-BOX GENES: DEVELOPMENT AND EVOLUTION OF PLANT BODY PLANS

We review functional data on MADS-box genes, recent phylogenetic analyses of these coding regions, and their roles in the development and evolution of key morphological innovations in plants. We map the origin of important morphological structures in particular diverse stages of the life cycle in different plant clades onto organismal phylogenies, and present relevant molecular genetic aspects of development related to the MADS-box genes. We focus on reproductive structures of the sporophyte because most functional characterizations have been done of MADS-box genes involved in flower development. We discuss MADS-box evolution in flowering plants, but we also review studies in the other nonflowering vascular plants, gymnosperms (conifers and gnetales), and ferns and preliminary data from the algae. We suggest that floral (e.g. flowering time, inflorescence, and flower meristem identity) MADS-box and nonfloral plant MADS-box genes should be the focus of future comparative research. Cloning and functional analyses of MADS-box genes in bryophytes, particularly in the experimental system Physcomitrella patens (Hedw.) B.S.G., are needed. The ABC model of floral organ specification is an excellent general representation of an important network of genes; however, formal analytical tools are required to integrate data on complex gene interaction in comparative analyses. This and other analytical approaches to constructing gene network models will help to frame homology hypotheses in an evolutionary and developmental framework.     

The ‘division of labour’ model of eye evolution

The ‘division of labour’ model of eye evolution is elaborated here. We propose that the evolution of complex, multicellular animal eyes started from a single, multi-functional cell type that existed in metazoan ancestors. This ancient cell type had at least three functions: light detection via a photoreceptive organelle, light shading by means of pigment granules and steering through locomotor cilia. Located around the circumference of swimming ciliated zooplankton larvae, these ancient cells were able to mediate phototaxis in the absence of a nervous system. This precursor then diversified, by cell-type functional segregation, into sister cell types that specialized in different subfunctions, evolving into separate photoreceptor cells, shading pigment cells (SPCs) or ciliated locomotor cells. Photoreceptor sensory cells and ciliated locomotor cells remained interconnected by newly evolving axons, giving rise to an early axonal circuit. In some evolutionary lines, residual functions prevailed in the specialized cell types that mirror the ancient multi-functionality, for instance, SPCs expressing an opsin as well as possessing rhabdomer-like microvilli, vestigial cilia and an axon. Functional segregation of cell types in eye evolution also explains the emergence of more elaborate photosensory–motor axonal circuits, with interneurons relaying the visual information.