裸子植物中光敏色素PHY-PAS1 结构域的适应性进化

光敏色素是一类红光/远红光受体, 在植物种子萌发到成熟的整个生长发育过程中均起重要的调节作用。光敏色素PHYPAS1 结构域存在于光敏色素基因家族的所有成员中, 对调节发色团的光谱特性和光信号转导非常关键。光敏色素基因家族通过基因重复产生, 而基因重复可能与物种形成有关。PHYP基因是裸子植物光敏色素基因家族发生第1次重复后产生的, 并且以单拷贝形式存在。为了研究不同裸子植物PHYP基因编码蛋白的PHY-PAS1结构域在进化过程中是否受到相同的选择压力以及是否发生了适应性进化, 该研究利用分支模型、位点模型以及分支-位点模型对裸子植物31条PHYP基因序列编码蛋白的PHY-PAS1结构域所受到的选择压力进行了分析。结果表明, 在由PHY-PAS1结构域序列构建的系统树中, 多数分支处于强烈的负选择压力下 (w<1); 有14个分支处于正选择压力下 (w>1), 其中13个分支发生在属内种间; 与之相比, 在较为古老的谱系中相对缺少这种正选择压力。     

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

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

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

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

蕨类植物rbcL基因正选择和负选择位点的鉴定

基于分支模型、位点模型及分支-位点模型对蕨类的rbcL基因所受到的选择压力进行了分析.结果显示:分支模型下检测到大部分分支处于负选择,仅4个分支处于正选择压力下,并且仅2分支具有统计上的显著性;在位点模型下,通过比较模型M1a/M2a和M7/M8,在氨基酸水平上模型M2a和模型M8均鉴定出98L位点被正选择;在模型M8下,鉴定负向选择位点共228个,占总序列的83.82%,从而揭示出负选择对rbcL基因的进化起着非常重要的作用;在分支位点-模型c下鉴定出262A被正选择.98L、262A位点分别位于rbcL羧基末端α/β桶结构域的第3和第8个α螺旋上.蕨类通过该结构域的适应性进化,适应白垩纪被子植物兴起而引发的陆地生态系统改变,研究结果为以后实验分析提供了首选位点.     

裸子植物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蛋白部分位点经历了适应性进化,通过位点间共进化网络协同作用方式辅助裸子植物响应陆生生境。     

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

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

高等植物光敏色素的分子结构、生理功能和进化特征

光敏色素是植物感受外界环境变化的最重要光受体之一,对红光和远红外光非常敏感。本文综述了光敏色素的分子结构、它所包含的结构域和相应功能以及植物各主要类群中光敏色素基因家族的成员组成与进化关系;重点在分子水平上介绍了光敏色素的生理功能与作用机制。最后,基于最新的研究进展提出了将来的研究方向。     

一种在核苷酸水平检测自然选择的新方法

非编码区序列在基因表达调控中起着重要作用,但其在进化过程中是否受到选择作用一直较难检测。最近有一些研究使用平均的核苷酸替换速率与中性序列的核苷酸替换速率的比值(ω)作为检测非编码区总体受选择作用的指标;但是对于非编码区而言,了解具体哪些核苷酸受到选择作用更具有意义。我们借鉴Nielsen&Yang(1998)检测单个氨基酸位点是否受选择作用的思路,在最大似然法的模型下,提出一种在核苷酸位点水平上对自然选择作用检测的方法。本方法能够检测在进化过程中对功能分化有重要贡献的核苷酸位点,包括编码和非编码区。将此方法应用于熟知的受到正选择作用的蛋白编码基因序列(HIV-1包装蛋白基因编码区),均能够检测到那些已知的受到正选择的核苷酸(密码子)位点,说明此方法可以有效地在核苷酸位点水平检测选择作用;又将此方法应用于非编码区(CTGF基因5′UTR),也得到了良好的结果。     

高等植物光敏色素的分子结构、生理功能和进化特征

光敏色素是植物感受外界环境变化的最重要光受体之一, 对红光和远红外光非常敏感。本文综述了光敏色素的分子结构、它所包含的结构域和相应功能以及植物各主要类群中光敏色素基因家族的成员组成与进化关系; 重点在分子水平上介绍了光敏色素的生理功能与作用机制。最后, 基于最新的研究进展提出了将来的研究方向。     

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

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

Adaptive Evolution in the GAF Domain of Phytochromes in Gymnosperms

The GAF domain of phytochrome is essential for photoconversion and signal transduction. In gymnosperms, it exists in all members of the phytochrome family that experience gene duplication. Maximum-likelihood models of codon substitution can provide a framework for constructing likelihood ratio tests of changes in selective pressure and make clear predictions about patterns of genetic change following gene duplication. In this study, 68 gymnosperm GAF sequences were analyzed to identify lineages and sites under positive selection. Our results indicate that (1) positive selection at a few sites (3.6%), rather than relaxation of selective constraints, has played a major role in the evolution of the gymnosperm GAF domain; (2) strong positive selective pressure tends to occur in the recent PHYP lineages of cogeneric species, but is absent in old lineages consisting of distantly related species; and (3) the selective pressure indicated by the ω ratio varies greatly among lineages and sites in the GAF domain.     

Positive selection during the diversification of class I vomeronasal receptor-like (V1RL) genes,putative pheromone receptor genes,in human and primate evolution

Vomeronasal receptors are the major receptors for pheromones in vertebrates, and five putative type 1 vomeronasal receptors (V1RL) have been identified in humans. The evolution of the V1RL1 gene in non-human primates, and patterns of selection on V1RL genes, were investigated. The presumed ortholog of V1RL1 was sequenced from 13 species of nonhuman primate, and in eight of these species V1RL1 was a pseudogene. Phylogenetic reconstructions reveal that V1RL1 pseudogene formation occurred independently in multiple primate lineages. Using maximum likelihood estimates of dN/dS ratios in PAML, we show that V1RL genes have evolved under neutral evolution in lineages in which they became a pseudogene. In contrast, among lineages in which V1RL genes contain an open reading frame, the majority of sites are under purifying selection and a minority are under significant positive selection. These results provide an interesting case where all three categories of selection can be teased apart in the same data set using maximum likelihood methods. The finding of positive selection on V1RL genes during primate evolution provides indirect support for the hypothesis that V1RL genes have a function in species-specific pheromone detection in primates.     

Physicochemical evolution and positive selection of the gymnosperm matK proteins

It is not clear whether matK evolves under Darwinian selection. In this study, the gymnosperm Taxaceae, Cephalotaxaceae and Pinaceae were used to illustrate the physicochemical evolution, molecular adaptation and evolutionary dynamics of gene divergence in matKs. matK sequences were amplified from 27 Taxaceae and 12 Cephalotaxaceae species. matK sequences of 19 Pinaceae species were retrieved from GenBank. The phylogenetic tree was generated using conceptual-translated amino acid sequences. Selective influences were investigated using standard d _N/d _S ratio methods and more sensitive techniques investigating the amino acid property changes resulting from nonsynonymous replacements in a phylogenetic context. Analyses revealed the presence of positive selection in matKs (N-terminal region, RT domain and domain X) of Taxaceae and Pinaceae, and found positive destabilizing selection in N-terminal region and RT domain of Cephalotaxaceae matK. Moreover, various amino acid properties were found to be influenced by destabilizing positive selection. Amino acid sites relating to these properties and to different secondary structures were found and have the potential to affect group II intron maturase function. Despite the evolutionary constraint on the rapidly evolving matK, this protein evolves under positive selection in gymnosperm. Several regions of matK have experienced molecular adaptation which fine-tunes maturase performance.     

Variation in the analysis of positively selected sites using nonsynonymous/synonymous rate ratios: an example using influenza virus

Sites in a gene showing the nonsynonymous/synonymous rate ratio (ω) >1 have been frequently identified to be under positive selection. To examine the performance of such analysis, sites of the ω ratio >1 in the HA1 gene of H3N2 subtype human influenza viruses were identified from seven overlapping sequence data sets in this study. Our results showed that the sites of the ω ratio >1 were of significant variation among the data sets even though they targeted similar clusters, indicating that the analysis is likely to be either of low sensitivity or of low specificity in identifying sites under positive selection. Most (43/45) of the sites showing ω >1 calculated from at least one data set are involved in B-cell epitopes which cover less than a half sites in the protein, suggesting that the analysis is likely to be of low sensitivity rather than of low specificity. It was further found that the analysis sensitivity could not be enhanced by including more sequences or covering longer time intervals. Previously some reports also likely identified only a portion of the sites under positive selection in the viral gene using the ω ratio. Low sensitivity of the analysis may result from that some sites under positive selection in the gene are also under negative (purifying) selection simultaneously for functional constrains, and so their ω ratios could be <1. Theoretically, the sites under the two opposite selection forces at the same time favor only certain nonsynonymous changes, e.g. those changing the antigenicity of the gene and maintaining the gene function. This study also suggested that sometimes we can identify more sites under positive selection using the ω ratio by integrating the positively selected sites estimated from multiple data sets.     

Evolutionary history and complementary selective relaxation of the duplicated PI genes in grasses

Gene duplication plays an important role in the evolution of organisms by allowing functional innovation and the divergence of duplicate genes. Previous studies found two PI-like genes in grass species, suggesting functional divergence between the paralogous copies. Here, we reconstructed the evolutionary history of two PI genes from major lineages of grasses and other monocot species, and demonstrated that two PI genes (PI1 and PI2) arose from a whole genome duplication that occurred in a common ancestor of extant grasses. Molecular evolutionary analyses at the family and tribal levels found strong purifying selection acting on two genes in grasses, consistent with the conserved class B function of the PI genes. Importantly, we detected different patterns of selective relaxation between the duplicated PI genes although no signature of positive selection was found. Likelihood ratio tests revealed that the ω ratio for M domain is significantly higher in PI1 than in PI2 but that for K domain is significantly higher in PI2 than in PI1. These findings imply that complementary selective relaxation occurs in two PI genes after duplication, and provide additional molecular evidence for the subfunctionalization of the duplicated PI genes in grasses.     

Rates and relations of mitochondrial genome evolution across the Echinoidea,with special focus on the superfamily Odontophora

In order to better characterize the placement of genus Tripneustes, as a representative of the Toxopneustidae family within the broader sea urchin mitochondrial (MT) phylogeny, the complete MT genome of Tripneustes gratilla was generated and compared with all published echinoid MT genomes currently available on NCBI GenBank. The MT genome phylogeny supports the existence of the superfamily Odontophora (consisting of the families Strongylocentrotidae, Echinometridae, and Toxopneustidae). A relaxed molecular‐clock time calibration suggests a split between the three key Odontophore MT lineages occurred during the late Eocene/Oligocene. Major global oceanographic changes have been inferred during this time frame, potentially driving species diversification through environmental selection pressures. To test for signatures of selection acting on the mitochondria, the historical rate of gene evolution of individual MT genes was assessed through a branch‐site comparison of nonsynonymous to synonymous substitution ratios (ω). Models of positive selection and neutral evolution, as compared via a likelihood ratio test, show no evidence of strong historical positive selection on mitochondrial genes at the genesis of the Odontophora. However, while pairwise ω comparison revealed signatures of strong negative selection, relatively elevated ω values were observed within the Strongylocentrotus genus.     

Accelerated evolution and loss of a domain of the sperm-egg-binding protein SED1 in ancestral primates

Proteins involved in sperm-egg binding have been shown to evolve rapidly in several groups of invertebrates and vertebrates. Mammalian SED1 (secreted protein containing N-terminal Notch-like type II epidermal growth factor (EGF) repeats and C-terminal discoidin/F5/8 C domains) is a recently identified sperm surface protein that binds the egg zona pellucida and facilitates sperm-egg adhesion. SED1-null male mice are subfertile. Here we examine the SED1 gene from 11 mammalian species and provide evidence that it underwent accelerated evolution in ancestral primates, most likely driven by positive selection. Specifically, the intensity of the positive selection across various protein domains of SED1 was heterogeneous. Although one of the 2 Notch-like EGF domains, which mediate protein-protein binding, was lost in primate SED1, the second EGF domain evolved under strong positive selection favoring polar to nonpolar amino acid replacements. By contrast, the 2 discoidin/F5/8 type C domains, which are involved in protein-cell membrane binding, do not show definite signs of positive selection. The structural modification and occurrence of directional selection in ancestral primates but not any other lineage suggest that the function of SED1 may have changed during primate evolution. These results reveal a different evolutionary pattern of SED1 from that of many other sperm-egg-binding proteins, which often show diversifying selection occurring in multiple lineages.     

A genomics approach to the detection of positive selection in cattle: adaptive evolution of the T-cell and natural killer cell-surface protein CD2

The detection of adaptive evolution at the molecular level is of interest not only as an insight into the process of evolution but also because of its functional implications for genes of interest. Here, we present the first genomics approach to detecting positive selection operating on the Bos taurus lineage, an important domestic species. This analysis led to the identification of the T-cell and natural killer (NK) cell receptor cluster of differentiation 2 (CD2) as having a strong signal of selection. Further detailed investigation of CD2 revealed that this gene was subject to positive selection during the evolution of a number of mammalian lineages. Moreover, we show that selection has operated primarily on the extracellular domain of CD2 and discuss the implications of this for an important regulator of the adaptive immune response.     

Positive selection in MAOA gene is human exclusive: determination of the putative amino acid change selected in the human lineage

Monoamine oxidase A (MAOA) is the X-linked gene responsible for deamination and subsequent degradation of several neurotransmitters and other amines. Among other activities, the gene has been shown to play a role in locomotion, circadian rhythm, and pain sensitivity and to have a critical influence on behavior and cognition. Previous studies have reported a non-neutral evolution of the gene attributable to positive selection in the human lineage. To determine whether this selection was human-exclusive or shared with other species, we performed a population genetic analysis of the pattern of nucleotide variation in non-human species, including bonobo, chimpanzee, gorilla, and orangutan. Footprints of positive selection were absent in all analyzed species, suggesting that positive selection has been recent and unique to humans. To determine which human-unique genetic changes could have been responsible for this differential evolution, the coding region of the gene was compared between human, chimpanzee, and gorilla. Only one human exclusive non-conservative change is present in the gene: Glu151Lys. This human substitution affects protein dimerization according to a three-dimensional structural model that predicts a non-negligible functional shift. This is the only candidate position at present to have been selected to fixation in humans during an episode of positive selection. Divergence analysis among species has shown that, even under positive selection in the human lineage, the MAOA gene did not experience accelerated evolution in any of the analyzed lineages, and that tools such as K_a/K_s would not have detected the selective history of the gene.