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

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

第三组胚胎晚期丰富蛋白lea3基因研究概述

晚期胚胎富集蛋白(late embryogenesis abundant protein,LEA蛋白)是在高等植物胚胎发育晚期大量积累的一类蛋白,根据其结构特点LEA蛋白一般分为6组,其中第3组LEA蛋白(LEA3)含有11个氨基酸串联重复的基元序列,可以形成α-螺旋结构,能在干旱胁迫的环境中保护生物大分子,减轻水份胁迫对植物造成的伤害,与植物抗逆性密切相关。该文就lea3基因及其蛋白的结构、功能、基因表达和应用等进行简要的综述,并对lea3基因及其蛋白今后的研究方向和应用前景进行了展望。     

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

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

极小种群濒危植物广西火桐、丹霞梧桐的叶绿体基因组特征

广西火桐(Firmiana kwangsiensis)和丹霞梧桐(F. danxiaensis)是我国南方特有物种, 其分布范围狭窄, 种群数量少。为了解其叶绿体基因组结构及系统发生关系, 本文通过高通量测序方法获得广西火桐和丹霞梧桐的浅层基因组数据, 通过生物信息学方法对叶绿体全基因组进行组装, 并对其结构特征进行分析。结果表明: 广西火桐和丹霞梧桐的叶绿体基因组大小分别为160,836 bp和161,253 bp, 具有典型被子植物叶绿体基因组环状四分体结构, 包含长度分别为89,700 bp、90,142 bp的大单拷贝区(large single copy, LSC), 长度分别为19,970 bp、20,067 bp的小单拷贝区(small single copy, SSC)及长度分别为25,583 bp、25,522 bp的2个反向重复序列区(inverted repeat sequence, IR)。两个物种的叶绿体基因组共注释得到131个基因, 包括86个蛋白编码基因、37个tRNA基因和8个rRNA基因。广西火桐的叶绿体基因组中共检测出26个正向重复序列、2个反向重复序列、21个回文重复序列、21个串联重复序列和98个简单重复序列; 丹霞梧桐叶绿体基因组中共检测出23个正向重复序列、5个反向重复序列、21个回文重复序列、30个串联重复序列和107个简单重复序列。系统发生分析结果表明5种梧桐属(Firmiana)植物构成两个强烈支持的分支(支持率100%), 一个分支为广西火桐、美丽火桐(F. pulcherrima)和火桐(F. colorata), 其中广西火桐与美丽火桐构成姐妹群; 另一分支是互为姐妹群的丹霞梧桐和云南梧桐(F. major)。综上所述, 广西火桐和丹霞梧桐的叶绿体基因组结构、基因排列及重复序列具有较高的相似性, 系统进化树将5种梧桐属物种分为两个分支, 其中广西火桐和美丽火桐最近; 而丹霞梧桐与云南梧桐关系最近。本研究鉴定的SSR位点可为梧桐属物种系统发生、进化关系的研究提供遗传信息。     

黄麂Mhc-DRB 基因多态性及其维持机制

利用牛DRB3 特异性引物(LA31 和LA32),通过聚合酶链式反应(PCR)、单链构象多态性(SSCP)以及克隆测序技术,从12 只黄麂个体中共获得20 个DRB 第二外显子等位基因,其中6 个个体具有3 ～ 4 个等位基因,提示利用该引物从黄麂中至少可以扩增出2 个DRB 位点。所有序列均无插入、缺失和终止密码子。基于序列比对(与牛DRB3 和鹿科DRB 基因同源性非常高),以及所检测到的氨基酸变异位点主要位于抗原结合区,推测本文所获得的黄麂序列为表达的、且具有重要功能的DRB 位点。抗原结合区氨基酸位点的非同义替换(d_N )显著大于同义替换(d_S )(P < 0.01),说明历史上黄麂DRB 基因经历过正选择作用。CODMEL 程序中的模型M7 和M8 似然比检测(Likelihood ratio test,LRT)结果同样支持上述推论。进一步利用经验贝叶斯法准确地检测出6 个受正选择作用的氨基酸位点(位点11、37、61、67、71、86),其中的5 个位点位于PBR 区。因此,正选择作用可能是维持黄麂DRB 基因多态性的主要机制之一。基于DRB 外显子2 序列利用邻接法(NJ)
构建了部分偶蹄动物系统发生关系,在NJ 树上,黄麂DRB 基因与其它鹿科动物DRB 基因呈镶嵌式分布,提示跨物种进化是维持黄麂DRB 基因多态性的另一重要机制。此外,黄麂两个等位基因(Mure-DRB1 和Mure-DRB11)和马鹿的两个等位基因(Ceel-DRB34 和Ceel-DRB46)与牛科的等位基因构成一个独立的进化枝,说明黄麂和马鹿的某些DRB 基因具有非常古老的谱系。     

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

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

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

念珠藻(Nostoc)固氮过程关键在于固氮酶的催化,而固氮酶复合物中的铁蛋白(NifH)是由高度保守的nifH基因编码的,该基因是进化史上现存最古老的功能基因之一。该研究选取念珠藻属及近缘类群的nifH基因序列共40条,采用最大似然法构建系统发育树;运行PAML4.9软件,对nifH基因编码蛋白进行生物信息学分析,并使用分支模型、位点模型和分支-位点模型检测该基因的选择位点,探讨nifH基因的适应性进化特征。结果表明:(1)最大似然树显示内类群中该研究物种共分为6个分支(A、B、C、D、E和F),其中D和E是2个大的分支,每个大分支中又各包含2个特殊的小分支A、F和B、C,其中F分支包含新疆古尔班通古特沙漠采集到的9株念珠藻,A分支包含F分支及该研究测定序列的4株葛仙米,B分支包含本研究测定序列的4株地皮菜和3株未定种的念珠藻,C分支包含NCBI数据库中下载的5株念珠藻、鱼腥藻序列和本研究测定序列的1株念珠藻。(2)在所分析的3种进化模型中,仅通过分支-位点模型检测出14个统计学上显著的正选择位点,即1F、2S、3S、4T、5A、6F、7F、8I、9S、10C、17I、27Y、29D和31R位点,表明念珠藻属植物的nifH基因发生了适应性变化,分支-位点模型是研究藻类基因适应性进化较好的模型。     

普通小麦低分子量麦谷蛋白亚基核心编码区与面团强度关系的研究

为分析LMW-GS基因对面团强度的影响, 利用两个重组自交系99G45/京771和Pm97034/京771的F9代,对LMW-GS基因特异位点和与其紧密连锁的Gli-1位点进行分析, 研究对面团强度影响差异显著的Glu-B3位点的等位基因核心编码区的序列差异。结果表明, 3个亲本LMW-GS核心编码区都具有6个半胱氨酸残基, 但PB较GB和JB缺失了一个7氨基酸序列的重复单元, 并且在不同序列中出现了氨基酸代换, 其中有2个代换可能影响氨基酸序列的亲水性, 进而影响面团强度。     

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

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

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

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

Genetic Diversity and Evolution of Satellite RNAs Associated with the Bamboo Mosaic Virus

Satellite RNAs (satRNAs) are subviral agents that depend on cognate helper viruses for genome replication and encapsidation. Their negative impacts on helper viruses have been exploited to control plant viral diseases. SatBaMV is a commonly found satRNA associated with Bamboo mosaic virus (BaMV) that infects diverse bamboo species in the field. To investigate the genetic diversity and evolution of satRNAs, we examined seven satBaMV populations derived from five bamboo species and cultivars from Taiwan, China, and India and one from the greenhouse. We found 3 distinct clades among the seven populations. Clade I is consisted of all satBaMV isolates, except for those from Dendrocalamus latiflorus in Taiwan and Bambusa vulgaris in India, which belong to Clades II and III, respectively. Interestingly, nucleotide diversity was lower for Clade I than II and III. However, the nucleotide diversity did not seem to depend on bamboo species or geographic location. Our population genetic analyses revealed the presence of excessive low-frequency polymorphic sites, which suggests that the satBaMV population was under purifying selection and/or population expansion. Further analysis of P20, the only satBaMV gene that encodes a non-structural protein involved in the long-distance movement of satBaMV, showed evidence of purifying selection. Taken together, our results suggest that purifying selection against defective P20 protein is responsible at least in part for the evolution of the satBaMV genome.     

Molecular evolution of the nontandemly repeated genes of the histone 3 multigene family.

In some species, histone gene clusters consist of tandem arrays of each type of histone gene, whereas in other species the genes may be clustered but not arranged in tandem. In certain species, however, histone genes are found scattered across several different chromosomes. This study examines the evolution of histone 3 (H3) genes that are not arranged in large clusters of tandem repeats. Although H3 amino acid sequences are highly conserved both within and between species, we found that the nucleotide sequence divergence at synonymous sites is high, indicating that purifying selection is the major force for maintaining H3 amino acid sequence homogeneity over long-term evolution. In cases where synonymous-site divergence was low, recent gene duplication appeared to be a better explanation than gene conversion. These results, and other observations on gene inactivation, organization, and phylogeny, indicated that these H3 genes evolve according to a birth-and-death process under strong purifying selection. Thus, we found little evidence to support previous claims that all H3 proteins, regardless of their genome organization, undergo concerted evolution. Further analyses of the structure of H3 proteins revealed that the histones of higher eukaryotes might have evolved from a replication-independent-like H3 gene.     

Identification and characterization of tandem repeats in exon III of dopamine receptor D4 (DRD4) genes from different mammalian species

In this study we have identified and characterized dopamine receptor D4 (DRD4) exon III tandem repeats in 33 public available nucleotide sequences from different mammalian species. We found that the tandem repeat in canids could be described in a novel and simple way, namely, as a structure composed of 15- and 12- bp modules. Tandem repeats composed of 18-bp modules were found in sequences from the horse, zebra, onager, and donkey, Asiatic bear, polar bear, common raccoon, dolphin, harbor porpoise, and domestic cat. Several of these sequences have been analyzed previously without a tandem repeat being found. In the domestic cow and gray seal we identified tandem repeats composed of 36-bp modules, each consisting of two closely related 18-bp basic units. A tandem repeat consisting of 9-bp modules was identified in sequences from mink and ferret. In the European otter we detected an 18-bp tandem repeat, while a tandem repeat consisting of 27-bp modules was identified in a sequence from European badger. Both these tandem repeats were composed of 9-bp basic units, which were closely related with the 9-bp repeat modules identified in the mink and ferret. Tandem repeats could not be identified in sequences from rodents. All tandem repeats possessed a high GC content with a strong bias for C. On phylogenetic analysis of the tandem repeats evolutionary related species were clustered into the same groups. The degree of conservation of the tandem repeats varied significantly between species. The deduced amino acid sequences of most of the tandem repeats exhibited a high propensity for disorder. This was also the case with an amino acid sequence of the human DRD4 exon III tandem repeat, which was included in the study for comparative purposes. We identified proline-containing motifs for SH3 and WW domain binding proteins, potential phosphorylation sites, PDZ domain binding motifs, and FHA domain binding motifs in the amino acid sequences of the tandem repeats. The numbers of potential functional sites varied pronouncedly between species. Our observations provide a platform for future studies of the architecture and evolution of the DRD4 exon III tandem repeat, and they suggest that differences in the structure of this tandem repeat contribute to specialization and generation of diversity in receptor function.     

城市绿化竹子生态适应性评价

针对城市绿化竹种要求,选择生长和抗性指标,用集对分析方法,分析和评价华安竹种园200个丛生竹和散混生竹种（或变种）生态适应性,结果表明：丛生竹种中,黄竹、小叶琴丝竹、糯竹、花竹、乡土竹、妈竹、青竿竹、泰竹等在华安生态适应性好,银丝大眼竹、孟竹、大绿竹、麻竹、白节勒竹、撑篙竹、石角竹、大木竹、小刺竹、勃氏甜龙竹、鱼肚腩竹等竹种生态适应性较差,其它竹种适应性中等.而散混生竹种中,衢县苦竹、实心竹、南平倭竹、矮若竹、乌哺鸡竹、毛竹、高节竹、红哺鸡竹、毛环水竹、斑苦竹、满山爆竹、斑箨茶秆竹等竹种适应性较好;缅甸方竹、四季竹、玉山竹篌竹、花毛竹、金镶玉竹等竹种适应性较差,其它竹种适应性中等.     

Positive diversifying selection on Plasmodium vivax RON2 protein

Plasmodium rhoptry neck protein 2 (RON2), which is released from the neck portion of the merozoite rhoptries and interacts with the microneme protein Apical Membrane Antigen 1 (AMA1), plays a crucial role in erythrocyte invasion. In this study, we sequenced the Plasmodium vivax RON2 gene from 19 P. vivax isolates collected in central China in order to establish whether this protein is under positive diversifying selection, which may occur as a result of protective host immune pressure?. In comparison with the P. vivax Sal-1 reference line, we found 10 amino acid substitutions dispersed throughout the open reading frame as well as indels caused by polymorphism in a repeat unit (21-23 repeats of (Q/E/K/N/H)(G/D)G(H/L/Y/P)G) in the second tandem repeat region located at amino acid positions 541-650. A McDonald-Kreitman test with RON2 sequences from the primate malaria parasite Plasmodium knowlesi, detected significant departure from neutrality in the PvRON2 3' region (nucleotide positions 2668-6609). These results suggest that the PvRON2 gene has evolved under positive diversifying selection.     

Evidence for heterogeneous selective pressures in the evolution of the env gene in different human immunodeficiency virus type 1 subtypes

Recent studies have demonstrated the emergence of human immunodeficiency virus type 1 (HIV-1) subtypes with various levels of fitness. Using heterogeneous maximum-likelihood models of adaptive evolution implemented in the PAML software package, with env sequences representing each HIV-1 group M subtype, we examined the various intersubtype selective pressures operating across the env gene. We found heterogeneity of evolutionary mechanisms between the different subtypes with a category of amino acid sites observed that had undergone positive selection for subtypes C, F1, and G, while these sites had undergone purifying selection in all other subtypes. Also, amino acid sites within subtypes A and K that had undergone purifying selection were observed, while these sites had undergone positive selection in all other subtypes. The presence of such sites indicates heterogeneity of selective pressures within HIV-1 group M subtype evolution that may account for the various levels of fitness of the subtypes.     

Mining of SSR markers from Expressed Sequence Tags of bamboo species

With the ever increasing number of Expressed Sequence Tags (ESTs) from various sequencing projects, ESTs have become valuable and first-hand source of in-silico mining of simple sequence repeats (SSR) markers. We examined a total of 3419 EST sequences from three bamboo species, namely, Phyllostachys edulis, Bambusa oldhamii and Dendrocalamus sinicus for the presence of di- to hexa- microsatellites. The frequency of SSR containing ESTs varied from 5.36% in B. oldhamii to 13.05% in P. edulis. No SSRs were found in D. sinicus. Tri-nucleotide repeats (49.34%) were most frequent in P. edulis, while not much comparable difference in repeats was found in B. oldhamii. Flanking primer pairs were also designed in-silico for the sequences containing SSRs and their position on the genome hypothesized using similarity searching. SSRs located in open reading frame (ORF) were given functional annotation using Gene Ontology. Polymorphic SSRs were also detected using new pipeline- polySSR. Polymorphism level was very low (2.43%) and the position of the polymorphic SSRs was determined. The development of SSRs and the study of polymorphism will help in the further study of intra- and inter- gene flow, genetic structure, variability, linkage mapping and evolutionary relationships in bamboo.     

Natural selection and mammalian BRCA1 sequences: elucidating functionally important sites relevant to breast cancer susceptibility in humans

Comparison of orthologous gene sequences is emerging as a powerful approach to elucidating functionally important positions in human disease genes. Using a diverse array of 132 mammalian BRCA1 (exon 11) sequences, we evaluated the functional significance of specific sites in the context of selection information (purifying, neutral, or diversifying) as well as the ability to extract such information from alignments that index varying degrees of mammalian diversity. Small data sets of either closely related taxa (Primates) or divergent placental taxa were unable to distinguish sites conserved due to purifying selection from sites conserved due to chance (false-positive rate = 65%-99%). Increasing the number of placental taxa to 57 greatly reduced the potential false-positive rate (0%-1.5%). Using the larger data set, we ranked the oncogenic risk of human missense mutations using a novel method that incorporates site-specific selection level and severity of the amino acid change evaluated against the amino acids present in other mammalian taxa. In addition to sites undergoing positive selection in Marsupialia, Laurasiatheria, Euarchontoglires, and Primates, we identified sites most likely to be undergoing divergent selection pressure in different lineages and six pairs of potentially interacting sites. Our results demonstrate the necessity of including large numbers of sequences to elucidate functionally important sites of a protein when using a comparative evolutionary approach.     

酿酒酵母中亚硫酸盐转运基因SSU1的分子进化分析

SSU1基因是涉及亚硫酸外排及SO₂耐受性的重要因素之一。为了研究酿酒酵母(Saccharomyce cerevisiae)中SSU1对SO₂耐受性及其分化机制, 文章探讨了SSU1基因在酿酒酵母中的遗传特征及进化规律。基于SSU1基因序列的聚类分析表明, 酿酒酵母群体可通过该基因分为3个亚群, 且与其分离的地理位置无关; 基于群体数据的McDonald-Kreitman 检验表明, SSU1基因在酿酒酵母中受到适应性选择的作用; Ka/Ks检验表明, 在酿酒酵母中, 不同的亚群间有Ka/Ks显著大于1 的值, 且PAML的支系模型检验到正选择作用在群体中特定的支系上; PAML的支系-位点模型检验获得9个潜在正选择作用位点, 其中有4个发生在受正选择作用的特定支系中; 基于ssu1p蛋白结构的分析表明, 在特定支系存在的正选择作用位点中, 除345(R/K)位点上两氨基酸替换均为碱性氨基酸外, 其他3个位点均是极性氨基酸/疏水性氨基酸之间替换, 考虑不同区域的氨基酸pKa值对其维持正常的功能有着重要的作用, 在该类位点的替换可能影响到ssu1p蛋白对SO₂的转运作用。