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

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

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

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

蕨类植物叶绿体基因ycf94的分子进化研究

ycf94基因是近年来在叶绿体基因组中新发现的一个基因,在蕨类植物中表现高度保守。该研究共选取94种蕨类植物,在系统发育背景下,对ycf94基因的结构特征、密码子偏好性、进化速率和适应性进化进行分析。结果表明, ycf94基因的密码子偏好性较弱,偏好使用以A/U结尾的密码子,且不同物种间的偏好性存在一定差异。密码子偏好性的形成主要受到突变压的影响,同时也存在其他因素的作用;基于凤尾蕨科和其他蕨类中ycf94基因的结构特征存在区别,对两者的分子替换速率进行了比较,表明颠换率、非同义替换率和ω值间存在显著差异;仅检测出1个正选择位点74A,强烈的负选择作用表明ycf94基因的结构和功能基本趋于稳定。这为蕨类系统发育分析提供了新依据,并提供了解析ycf94基因功能的线索。     

厚囊蕨类植物VA菌根的初步研究

对１０种厚囊蕨类植物的ＶＡ菌根状况以及这些植物根际土壤中的ＶＡ菌根真菌进行了研究,结果发现１０种植物的根部全部感染有ＶＡ菌根。从这１０种植物的根际土壤中分离鉴定了分属于３个属的９种ＶＡ菌根真菌,其中,球囊霉属和无梗囊霉真菌是这些植物的根际土壤中的优势类群。ＶＡ菌根对厚囊蕨类植物的进化可能具有重要的意义。     

蕨类植物叶绿体rps12基因的分子进化研究

以68种蕨类植物和2种石松类植物的rps12基因为对象,在系统发育背景下,结合最大似然法,使用HyPhy和PAML软件对该基因进行进化速率和适应性进化研究。结果显示:位于IR区的外显子2~3,其替换率明显降低,rps12基因编码序列的替换率也随之降低,且rps12基因密码子第3位的GC含量明显升高;在蕨类植物的进化过程中,3′-rps12更倾向定位于IR区,以保持较低的替换率;rps12基因编码的123个氨基酸位点中,共检测到4个正选择位点和116个负选择位点。研究结果表明基因序列进入到IR区后,显示出降低的替换率;强烈的负选择压力表明RPS12蛋白的高度保守性以及rps12基因的功能和结构已经趋于稳定。     

蕨类植物叶绿体rps4基因的适应性进化分析

在原核生物和植物叶绿体中,RPS4(ribosomal protein small subunit4)在核糖体30S小亚基形成起始过程中发挥重要作用;该蛋白在植物中由叶绿体rps4基因编码。为验证蕨类植物在白垩纪适应被子植物兴起而发生分化的观点,本文以23种蕨类植物为研究对象,利用分支模型、位点模型和分支位点模型对其叶绿体rps4基因进化适应性进行分析。分支模型检测到4个可能存在正选择的分支;位点模型和分支位点模型虽然没有检测出正选择位点,但是位点模型检测出了85个负选择位点。通过研究我们仅仅得出a、b两个代表水龙骨类的分支处于正选择压力下,这与水龙骨类在白垩纪发生辐射式演化的理论相一致。同时rps4基因处于强烈的负选择压力这一事实表明该基因的功能与结构已经趋于稳定。     

蕨类植物和裸子植物的起源与进化

最近笔者拜读了由中山大学生物系与南京大学生物系合编的《植物学》(系统、分类部分),高等教育出版社,1984年3月第四次印刷本,发现在该书的157—188页,即于论述“蕨类植物的起源与进化”和“裸子植物的起源与进化”两节中,存在一些值得与作者商榷的问题,现评论如下,谨希同行专家广为讨论。     

凤尾蕨科旱生蕨类rbcL基因的适应性进化和共进化分析

核酮糖-1,5-二磷酸羧化酶/加氧酶（Rubisco,EC 4.1.1.39）是植物参与光合作用的关键酶,其大亚基由叶绿体rbcL基因编码。为深入理解凤尾蕨科植物对干旱生境的分子适应机制,本研究以53种凤尾蕨科旱生植物的rbcL基因为对象,展开适应性进化和共进化研究。采用位点间可变ω比值模型以及SLAC、REL和FEL等方法进行的适应性进化分析显示：在氨基酸水平上共有15个正选择位点（66S、84E、139L、235G、245I、252A、273Y、295K、296N、299M、307G、330E、349S、365F、404A）,其中位点245I、252A和273Y对维持Rubisco功能起重要作用。共进化分析共鉴定出2组共进化位点,分别由139L、273Y、295K和273Y、295K、349S组成,这些氨基酸位点间的共进化方式与蛋白质的疏水性和分子量都显著相关。以上结果一方面支持基于ω比值检验DNA编码序列发生适应性进化的有效性,另一方面也提示凤尾蕨科植物对干旱生境的适应可能与rbcL基因的适应性进化有关。     

蕨类植物的VA菌根与蕨类植物系统演化关系的研究

通过对云南热带、亚热带生长的２５６种蕨类植物ＶＡ菌根的调查,发现蕨类植物ＶＡ菌根营养者所占的比例低于被子植物;在真蕨类植物中,植物具有由ＶＡ菌根营养经兼性ＶＡ菌根营养向自养方向进化的趋势。     

裸子植物psbA基因分子进化式样的研究

为阐明裸子植物对陆生生境生态响应的分子机制,以新近的裸子植物分类系统为指导,基于psb A基因编码全序列对4亚纲53种代表植物进行分子进化分析。首先,依据"放松分子钟"模型重建裸子植物在时间尺度下系统发育关系;其次,采用6个模型(MEC/JTT、MEC/cp REV、M5、M7、M8、M8a)估测氨基酸位点ω值,并对各模型结果进行统计检测;随后,利用Bootstrap方法检PSBA蛋白内部氨基酸位点的共进化动态。结果表明,系统树提示的物种分化历程支持前期分类结果;光合系统反应中心核心PSBA蛋白有3个氨基酸位点(13、19和243)曾经受正选择压力;PSBA蛋白内部有多对氨基酸位点间构成了共进化网络。因此,psb A基因编码序列具有作为描绘裸子植物系统发育关系标记的潜力,PSBA蛋白部分位点经历了适应性进化,通过位点间共进化网络协同作用方式辅助裸子植物响应陆生生境。     

Phylogenetic relationships of ferns deduced from rbcL gene sequence

Part of the large subunit of the ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco) gene (rbcL) was sequenced from three fern species: Adiantum capillus-veneris, Botrypus strictus, and Osmunda cinnamomea var. fokiensis. This region included 1,333 base pairs, about 90% of the gene. Maximum likelihood analysis of the deduced amino acid sequences indicated that (1) Botrypus (Ophioglossaceae) clustered monophyletically with other ferns (Adiantum, Angiopteris, Osmunda); the closest relative to Botrypus among the three species was Osmunda, which did not support the hypothesis that the Ophioglossaceae are linked to the progymnosperm-seed plant lineage. (2) Eusporangiate ferns containing Botrypus (Ophioglossaceae) and Angiopteris (Marattiaceae) were a paraphyletic group. (3) Seed plants and the four fern species examined formed a monophyletic group, but ferns and bryophytes (liverwort) did not. Variations in rates of substitution for synonymous and nonsynonymous codons were found in fern lineages.Correspondence to: M. Hasebe     

Adaptive evolution in the rat olfactory receptor gene family

Summary Comparison of DNA sequences of the rat (Rattus norvegicus) olfactory receptor gene family revealed an unusual pattern of nucleotide substitution in the gene region encoding the second extracellular domain (E2) of the protein. In this domain, nonsynonymous nucleotide differences between members of this subfamily that caused a change in amino acid residue polarity were over four times more frequent than nonsynonymous differences that did not cause a polarity change. This nonrandom pattern of nucleotide substitution is evidence of past directional selection favoring diversification of the E2 domain among members of this subfamily. This in turn suggests that E2 may play some important role in the functions unique to each member of the olfactory receptor family, and that it may perhaps be an odorant binding domain.Offprint requests to: A.L. Hughes     

Variable evolutionary rates in the molecular evolution of mammalian growth hormones

In mammals pituitary growth hormone (GH) shows a slow basal rate of evolution (0.22 ± 0.03 × 10^–9 substitutions/amino acid site/year) which appears to have increased by at least 25–50-fold on two occasions, during the evolution of primates (to at least 10.8 ± 1.3 X 10^–9 substitutions/amino acid site/year) and artiodactyl ruminants (to at least 5.6 ± 1.3 X 10^–9 substitutions/amino acid site/year). That these rate increases are real, and not due to inadvertent comparison of nonorthologous genes, was established by showing that features of the GH gene sequences that are not expressed as mature hormone do not show corresponding changes in evolutionary rate. Thus, analysis of nonsynonymous substitutions in the coding sequence for the mature protein confirmed the rate increases seen in the primate and ruminant GHs, but analysis of nonsynonymous substitutions in the signal peptide sequence, synonymous substitutions in the coding sequence for signal peptide or mature protein, and 5 and 3 untranslated sequences showed no statistically significant changes in evolutionary rate. Evidence that the increases in evolutionary rate are probably due to positive selection is provided by the observation that in the cases of both ruminant and primate GHs the periods of rapid evolution were followed by a return to a slow rate similar to the basal rate seen in other mammalian GHs. Differences between the biological properties of GHs have been identified which may relate to these periods of rapid adaptive molecular evolution. On the basis of sequence data currently available (but excluding rodent GHs which show an intermediate rate, the basis of which is not clear) for most (90%) of evolutionary time mammalian GHs have been in the slow phase of evolution, with possibly most of the few amino acid substitutions that have occurred being neutral in nature. But most (80%) of the amino acid substitutions that have been introduced into GH during the course of mammalian evolution have been accepted during the rapid phases and were adaptive in nature.     

A hardware interpretation of the evolution of the genetic code

A quantitative rationale for the evolution of the genetic code is developed considering the principle of minimal hardware. This principle defines an optimal code as one that minimizes for a given amount of information encoded, the product of the number of physical devices used by the average complexity of each device. By identifying the number of different amino acids, number of nucleotide positions per codon and number of base types that can occupy each such position with, respectively, the amount of information, number of devices and the complexity, we show that optimal codes occur for 3, 7 and 20 amino acids with codons having a single, two and three base positions per codon, respectively. The advantage of a code of exactly 4 symbols is deduced, as well as a plausible evolutionary pathway from a code of doublets to triplets. The present day code of 20 amino acids encoded by 64 codons is shown to be the most optimal in an absolute sense. Using a tetraplet code further evolution to a code in which there would be 55 amino acids is in principle possible, but such a code would deviate slightly more than the present day code from the minimal hardware configuration. The change from a triplet code to a tetraplet code would occur at about 32 amino acids. Our conclusions are independent of, but consistent with, the observed physico-chemical properties of the amino acids and codon structures. These correlations could have evolved within the constrains imposed by the minimal hardware principle.     

Parallel evolution of IDH2 gene in cetaceans,primates and bats

Cetaceans and primates both have large brains that require large amounts of aerobic energy metabolism. In bats, the cost of flight makes locomotion energetically demanding. These mammalian groups may represent three independent evolutionary origins of an energy-demanding lifestyle in mammals. IDH2 encodes an enzyme in the tricarboxylic acid cycle in the mitochondrion, which plays a key role in aerobic energy metabolism. In this study, we cloned and sequenced this gene in two cetaceans, and 19 bat species, and compared the data with available primate sequences to test its evolution. We found significant signals of parallel evolution in this gene among these three groups. Parallel evolution of this gene may reflect their parallel evolution towards a higher demand for energy.     

Alleles versus mutations: Understanding the evolution of genetic architecture requires a molecular perspective on allelic origins

Perspectives on the role of large‐effect quantitative trait loci (QTL) in the evolution of complex traits have shifted back and forth over the past few decades. Different sets of studies have produced contradictory insights on the evolution of genetic architecture. I argue that much of the confusion results from a failure to distinguish mutational and allelic effects, a limitation of using the Fisherian model of adaptive evolution as the lens through which the evolution of adaptive variation is examined. A molecular‐based perspective reveals that allelic differences can involve the cumulative effects of many mutations plus intragenic recombination, a model that is supported by extensive empirical evidence. I discuss how different selection regimes could produce very different architectures of allelic effects under a molecular‐based model, which may explain conflicting insights on genetic architecture from studies of variation within populations versus between divergently selected populations. I address shortcomings of genome‐wide association study (GWAS) practices in light of more suitable models of allelic evolution, and suggest alternate GWAS strategies to generate more valid inferences about genetic architecture. Finally, I discuss how adopting more suitable models of allelic evolution could help redirect research on complex trait evolution toward addressing more meaningful questions in evolutionary biology.     

Rates of molecular evolution in tree ferns are associated with body size,environmental temperature,and biological productivity

Variation in rates of molecular evolution (heterotachy) is a common phenomenon among plants. Although multiple theoretical models have been proposed, fundamental questions remain regarding the combined effects of ecological and morphological traits on rate heterogeneity. Here, we used tree ferns to explore the correlation between rates of molecular evolution in chloroplast DNA sequences and several morphological and environmental factors within a Bayesian framework. We revealed direct and indirect effects of body size, biological productivity, and temperature on substitution rates, where smaller tree ferns living in warmer and less productive environments tend to have faster rates of molecular evolution. In addition, we found that variation in the ratio of nonsynonymous to synonymous substitution rates (dN/dS) in the chloroplast rbcL gene was significantly correlated with ecological and morphological variables. Heterotachy in tree ferns may be influenced by effective population size associated with variation in body size and productivity. Macroevolutionary hypotheses should go beyond explaining heterotachy in terms of mutation rates and instead, should integrate population‐level factors to better understand the processes affecting the tempo of evolution at the molecular level.