膜蕨科植物rbcL基因的适应性进化和共进化分析

膜蕨科植物是薄囊蕨类中种类最多的科,主要分布在潮湿的热带地区,拥有陆生、附生、半附生和攀生等多种生态型。为进一步了解膜蕨科植物辐射式物种分化的分子适应机制,该研究在时间框架下采用位点模型对膜蕨科植物rbc L基因的进化式样进行分析。结果表明:共鉴定出6个氨基酸正选择位点(125I、227L、231A、258F、304S和351L),其中位点304S位于环六上,对维持Rubisco功能有重要作用。此外,还计算了Rubisco大亚基内部氨基酸位点之间的共进化关系,共检测出39组(35个氨基酸)共进化位点,其中位点在α螺旋上的占46%,在β折叠上的占14%。膜蕨科植物rbc L基因这种复杂的进化式样可能与其起源较早有关。鉴于此,基于UCLD分子钟模型对膜蕨科植物的分化时间进行了估计,结果显示膜蕨科植物首次发生分歧的时间在三叠纪早期,瓶蕨属和膜蕨属的分歧时间分别发生在侏罗纪早期和白垩纪晚期,并且得出陆生生态型是其它生态型进化的基础,推测最近几次最热事件可能对物种分化的形成产生一定的作用。该研究结果对认识膜蕨科植物如何应对被子植物兴起所导致的陆地生态系统改变具重要意义。     

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

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

凤尾蕨科旱生蕨类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基因的适应性进化有关。     

大黄属(蓼科) 植物ｎｄｈＦ 基因的适应性进化

大黄属(Rheum L.)是蓼科(Polygonaceae)中一个高度分化的大属,广泛分布在亚洲和欧洲的高山和沙漠地区,全世界约60种,其中在青藏高原及其邻近地区发现了约40种。该属种的高度分化曾被推测是第三纪末青藏高原的快速隆升以及第四纪气候的反复变化所引发的适应性辐射导致。为进一步了解大黄属植物辐射式物种分化的分子适应机制,该研究选取34个形态上多样化的大黄属物种,利用系统发育分析软件,在时间框架下采用位点模型和分支模型对大黄属的叶绿体ndhF基因进行了适应性进化分析。结果表明：大黄属植物的分子进化系统树呈现短而平行的辐射式分支式样,显示出典型的物种快速辐射多样化特征;用位点模型检验ndhF基因是否存在经受正向选择(ω>1)时,在氨基酸水平上共鉴定出3个NDHF亚基的正选择位点(188H,465H,551L),对NDHF亚基的二级结构进行分析后发现编码的188H氨基酸位于α螺旋上。大黄属植物可能通过这些结构域的适应性进化,适应青藏高原的快速隆升以及第四纪气候的反复变化而引发的陆地生态系统改变。该研究结果可为今后对该属植物的实验分析提供首选位点。     

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

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

大黄属(蓼科)植物ndhF基因的适应性进化(英文)

大黄属(Rheum L.)是蓼科(Polygonaceae)中一个高度分化的大属,广泛分布在亚洲和欧洲的高山和沙漠地区,全世界约60种,其中在青藏高原及其邻近地区发现了约40种。该属种的高度分化曾被推测是第三纪末青藏高原的快速隆升以及第四纪气候的反复变化所引发的适应性辐射导致。为进一步了解大黄属植物辐射式物种分化的分子适应机制,该研究选取34个形态上多样化的大黄属物种,利用系统发育分析软件,在时间框架下采用位点模型和分支模型对大黄属的叶绿体ndhF基因进行了适应性进化分析。结果表明:大黄属植物的分子进化系统树呈现短而平行的辐射式分支式样,显示出典型的物种快速辐射多样化特征;用位点模型检验ndhF基因是否存在经受正向选择(ω1)时,在氨基酸水平上共鉴定出3个NDHF亚基的正选择位点(188H,465H,551L),对NDHF亚基的二级结构进行分析后发现编码的188H氨基酸位于α螺旋上。大黄属植物可能通过这些结构域的适应性进化,适应青藏高原的快速隆升以及第四纪气候的反复变化而引发的陆地生态系统改变。该研究结果可为今后对该属植物的实验分析提供首选位点。     

细鳞苔科psbA基因的适应性进化分析

在时间框架下,采用机理式模型(Mechanistic model)和MEC模型(Mechanistic-empirical combination model)以及Datamonkey对细鳞苔科psbA基因的进化式样进行了分析.结果均未检测到统计上显著的正选择位点,显示负选择对细鳞苔科psbA基因起主导作用.另外,基于UCLD分子钟估算出的细鳞苔科各分支分歧时间表明,该科物种丰富度的辐射式增长发生在新生代渐新世.     

中坡国家森林公园蕨类植物区系分析

对中坡国家森林公园的蕨类植物分类、区系组成及其特点进行了研究．结果表明：中坡共有蕨类植物28科53属85种,主要科为鳞毛蕨科、水龙骨科、凤尾蕨科和金星蕨科;主要属为凤尾蕨属、鳞毛蕨属、耳蕨属和卷柏属;科的分布类型以世界分布和泛热带分布为主,属以泛热带属和热带亚洲属最多,种以东亚分布类型为主并表现出明显的亚热带特点;生态类型主要为土生,其次是石生类型．蕨类区系与贵州的关系最为密切．蕨类物种密度为每km^26．2种,在纬度相近的9个保护区中仅次于云南西山．中坡拥有桫椤和金毛狗等2种国家重点保护蕨类和7种中国特有的蕨类．桫椤在中坡的发现,对于研究其演变具有一定的意义．     

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

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

湖南壶瓶山国家级保护区药用蕨类植物的调查研究

通过野外调查、标本采集与鉴定和查阅文献,对壶瓶山国家级保护区药用蕨类植物多样性进行统计,并对其资源分布,药用部分、功能分类及民族文化等方面进行了探讨。结果如下：（1）壶瓶山国家级保护区有药用蕨类165种（不包括变种和变型）,隶属于33科57属,分别占武陵山地区和全国药用蕨类总科属种的80.48%7,1.25%,77.10%和67.34%,47.90%3,6.91%;含种数较多的科属有水龙骨科Polypodiaceae,鳞毛蕨科Dryopteridaceae,凤尾蕨科Pteridaceae,卷柏属Selaginella,凤尾蕨属Pteris,鳞盖蕨属Microlepia等。（2）该区药用蕨类植物原始和进化的类群共存,区系以热带成分为主,具有热带亲缘性,体现了壶瓶山药用蕨类植物起源古老,但物种分化强烈的特点。（3）该区的药用蕨类在四个植被垂直带谱都有分布,其中在低于1 000 m的常绿阔叶林带分布最多。（4）该区药用蕨类根据药用部位分为四类,全草类最多,其次是根茎类;按药用功能将分为五类,主要为清热类和解毒类。（5）讨论了壶瓶山国家级保护区独特的地理位置和优厚的自然环境是该区药用蕨类植物多样性高的重要原因,并提出了可持续开发利用与保护的建议。     

A random effects branch-site model for detecting episodic diversifying selection

Adaptive evolution frequently occurs in episodic bursts, localized to a few sites in a gene, and to a small number of lineages in a phylogenetic tree. A popular class of "branch-site" evolutionary models provides a statistical framework to search for evidence of such episodic selection. For computational tractability, current branch-site models unrealistically assume that all branches in the tree can be partitioned a priori into two rigid classes--"foreground" branches that are allowed to undergo diversifying selective bursts and "background" branches that are negatively selected or neutral. We demonstrate that this assumption leads to unacceptably high rates of false positives or false negatives when the evolutionary process along background branches strongly deviates from modeling assumptions. To address this problem, we extend Felsenstein's pruning algorithm to allow efficient likelihood computations for models in which variation over branches (and not just sites) is described in the random effects likelihood framework. This enables us to model the process at every branch-site combination as a mixture of three Markov substitution models--our model treats the selective class of every branch at a particular site as an unobserved state that is chosen independently of that at any other branch. When benchmarked on a previously published set of simulated sequences, our method consistently matched or outperformed existing branch-site tests in terms of power and error rates. Using three empirical data sets, previously analyzed for episodic selection, we discuss how modeling assumptions can influence inference in practical situations.     

Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level

Detecting positive Darwinian selection at the DNA sequence level has been a subject of considerable interest. However, positive selection is difficult to detect because it often operates episodically on a few amino acid sites, and the signal may be masked by negative selection. Several methods have been developed to test positive selection that acts on given branches (branch methods) or on a subset of sites (site methods). Recently, Yang, Z., and R. Nielsen (2002. Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. Mol. Biol. Evol. 19:908-917) developed likelihood ratio tests (LRTs) based on branch-site models to detect positive selection that affects a small number of sites along prespecified lineages. However, computer simulations suggested that the tests were sensitive to the model assumptions and were unable to distinguish between relaxation of selective constraint and positive selection (Zhang, J. 2004. Frequent false detection of positive selection by the likelihood method with branch-site models. Mol. Biol. Evol. 21:1332-1339). Here, we describe a modified branch-site model and use it to construct two LRTs, called branch-site tests 1 and 2. We applied the new tests to reanalyze several real data sets and used computer simulation to examine the performance of the two tests by examining their false-positive rate, power, and robustness. We found that test 1 was unable to distinguish relaxed constraint from positive selection affecting the lineages of interest, while test 2 had acceptable false-positive rates and appeared robust against violations of model assumptions. As test 2 is a direct test of positive selection on the lineages of interest, it is referred to as the branch-site test of positive selection and is recommended for use in real data analysis. The test appeared conservative overall, but exhibited better power in detecting positive selection than the branch-based test. Bayes empirical Bayes identification of amino acid sites under positive selection along the foreground branches was found to be reliable, but lacked power.     

Tryptophan scanning mutagenesis in the TM3 domain of the Torpedo californica acetylcholine receptor beta subunit reveals an alpha-helical structure

We used tryptophan substitutions to characterize the beta M3 transmembrane domain (betaTM3) of the acetylcholine receptor (AChR). We generated 15 mutants with tryptophan substitutions within the betaTM3 domain, between residues R282W and I296W. The various mutants were injected into Xenopus oocytes, and expression levels were measured by [125I]-alpha-bungarotoxin binding. Expression levels of the M288W, I289W, L290W, and F293W mutants were similar to that of wild type, whereas the other mutants (R282W, Y283W, L284W, F286W, I287W, V291W, A292W, S294W, V295W, and I296W) were expressed at much lower levels than that of wild type. None of these tryptophan mutants produced peak currents larger than that of wild type. Five of the mutants, L284W, F286W, I287W, V295W, and I296W, were expressed at levels <15% of the wild type. I296W had the lowest expression levels and did not display any significant ACh-induced current, suggesting that this position is important for the function and assembly of the AChR. Tryptophan substitution at three positions, L284, V291, and A292, dramatically inhibited AChR assembly and function. A periodicity analysis of the alterations in AChR expression at positions 282-296 of the betaTM3 domain was consistent with an alpha-helical structure. Residues known to be exposed to the membrane lipids, including R282, M285, I289, and F293, were all found in all the upper phases of the oscillatory pattern. Mutants that were expressed at lower levels are clustered on one side of a proposed alpha-helical structure. These results were incorporated into a structural model for the spatial orientation of the TM3 of the Torpedo californica beta subunit.     

Adenosine is inherently favored as the branch-site RNA nucleotide in a structural context that resembles natural RNA splicing

We previously used in vitro selection to identify the 7S11 deoxyribozyme, which catalyzes formation of 2',5'-branched RNA using a branch-site adenosine nucleophile and a 5'-triphosphate electrophile. An unanswered question is whether the use of branch-site adenosine is inherently preferred or a chance event during the particular selection experiment. Here we have found that deoxyribozymes newly selected to use uridine as the branch-site RNA nucleotide in a structural context that resembles natural RNA splicing instead prefer a branch-site adenosine, although adenosine was never available during the selection itself. Our results support a chemical basis for nature's choice of the branch-site nucleotide, which is almost always adenosine in group II introns and the spliceosome.     

Dynamic Contacts of U2, RES,Cwc25, Prp8 and Prp45 Proteins with the Pre-mRNA Branch-Site and 3' Splice Site during Catalytic Activation and Step 1 Catalysis in Yeast Spliceosomes

Little is known about contacts in the spliceosome between proteins and intron nucleotides surrounding the pre-mRNA branch-site and their dynamics during splicing. We investigated protein-pre-mRNA interactions by UV-induced crosslinking of purified yeast B^act spliceosomes formed on site-specifically labeled pre-mRNA, and analyzed their changes after conversion to catalytically-activated B* and step 1 C complexes, using a purified splicing system. Contacts between nucleotides upstream and downstream of the branch-site and the U2 SF3a/b proteins Prp9, Prp11, Hsh49, Cus1 and Hsh155 were detected, demonstrating that these interactions are evolutionarily conserved. The RES proteins Pml1 and Bud13 were shown to contact the intron downstream of the branch-site. A comparison of the B^act crosslinking pattern versus that of B* and C complexes revealed that U2 and RES protein interactions with the intron are dynamic. Upon step 1 catalysis, Cwc25 contacts with the branch-site region, and enhanced crosslinks of Prp8 and Prp45 with nucleotides surrounding the branch-site were observed. Cwc25’s step 1 promoting activity was not dependent on its interaction with pre-mRNA, indicating it acts via protein-protein interactions. These studies provide important insights into the spliceosome''s protein-pre-mRNA network and reveal novel RNP remodeling events during the catalytic activation of the spliceosome and step 1 of splicing.     

Analysis of transmembrane segment 7 of the dipeptide transporter hPepT1 by cysteine-scanning mutagenesis

To investigate the involvement of transmembrane segment 7 (TMS7) of hPepT1 in forming the putative central aqueous channel through which the substrate traverses, we individually mutated each of the 21 amino acids in TMS7 to a cysteine and analyzed the mutated transporters using the scanning cysteine accessibility method. Y287C- and M292C-hPepT1 did not express at the plasma membrane. Out of the remaining 19 transporters, three (F293C-, L296C-, and F297C-hPepT1) showed negligible glycyl-sarcosine (gly-sar) uptake activity and may play an important role in defining the overall hPepT1 structure. K278C-hPepT1 showed approximately 40% activity and the 15 other transporters exhibited more than 50% gly-sar uptake when compared with wild type (WT)-hPepT1. Gly-sar uptake for the 16 active transporters containing cysteine mutations was then measured in the presence of 2.5 mM 2-aminoethyl methanethiosulfonate hydrobromide (MTSEA) or 1 mM [2-(trimethylammonium) ethyl] methanethiosulfonate bromide (MTSET). Gly-sar uptake was significantly inhibited for each of the 16 single cysteine mutants in the presence of 2.5 mM MTSEA. In contrast, significant inhibition of uptake was only observed for K278C-, M279C-, V280C-, T281C-, M284C-, L286C-, P291C-, and D298C-hPepT1 in the presence of 1 mM MTSET. MTSET modification of R282C-hPepT1 resulted in a significant increase in gly-sar uptake. To investigate this further, we mutated WT-hPepT1 to R282A-, R282E-, and R282K-hPepT1. R282E-hPepT1 showed a 43% reduction in uptake activity, whereas R282A- and R282K-hPepT1 had activities comparable with WT-hPepT1, suggesting a role for the Arg-282 positive charge in substrate translocation. Most of the amino acids that were MTSET-sensitive upon cysteine mutation, including R282C, are located toward the intracellular end of TMS7. Hence, our results suggest that TMS7 of hPepT1 is relatively solvent-accessible along most of its length but that the intracellular end of the transmembrane domain is particularly so. From a structure-function perspective, we speculate that the extracellular end of TMS7 may shift following substrate binding, providing the basis for channel opening and substrate translocation.     

Episodic positive selection during the evolution of naphthalene dioxygenase to nitroarene dioxygenase

Using different maximum-likelihood models of adaptive evolution, signatures of natural selective pressure, operating across the naphthalene family of dioxygenases, were examined. A lineage- and branch-site specific combined analysis revealed that purifying selection pressure dominated the evolutionary history of the enzyme family. Specifically, episodic positive Darwinian selection pressure, affecting only a few sites in a subset of lineages, was found to be responsible for the evolution of nitroarene dioxygenases (NArDO) from naphthalene dioxygenase (NDO). Site-specific analysis confirmed the absence of diversifying selection pressure at any particular site. Different sets of positively selected residues, obtained from branch-site specific analysis, were detected for the evolution of each NArDO. They were mainly located around the active site, the catalytic pocket and their adjacent regions, when mapped onto the 3D structure of the α-subunit of NDO. The present analysis enriches the current understanding of adaptive evolution and also broadens the scope for rational alteration of substrate specificity of enzyme by directed evolution.     

More than one way to splice an RNA: branching without a bulge and splicing without branching in group II introns.

Domain 6 (D6) of group II introns contains a bulged adenosine that serves as the branch-site during self-splicing. In addition to this adenosine, other structural features in D6 are likely to contribute to the efficiency of branching. To understand their role in promoting self-splicing, the branch-site and surrounding nucleotides were mutagenized. Detailed kinetic analysis on the self-splicing efficiency of the mutants revealed several interesting features. First, elimination of the branch-site does not preclude efficient splicing, which takes place instead through a hydrolytic first step. Second, pairing of the branch-site does not eliminate branching, particularly if the adenosine is involved in a mispair. Third, the G-U pairs that often surround group II intron branch-points contribute to the efficiency of branching. These results suggest that there is a strong driving force for promoting self-splicing by group II introns, which employ a versatile set of different mechanisms for ensuring that splicing is successful. In addition, the behavior of these mutants indicates that a bulged adenosine per se is not the important determinant for branch-site recognition in group II introns. Rather, the data suggest that the branch-site adenosine is recognized as a flipped base, a conformation that can be promoted by a variety of different substructures in RNA and DNA.