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

以68种蕨类植物和2种石松类植物的rps12基因为对象,在系统发育背景下,结合最大似然法,使用Hy Phy和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基因处于强烈的负选择压力这一事实表明该基因的功能与结构已经趋于稳定。     

蕨类植物rpoC1内含子缺失及其分子进化速率

rpoC1基因编码RNA聚合酶β￠亚基蛋白,在转录过程中与DNA模板结合,与β亚基形成的β-β￠亚基复合体构成RNA合成的催化中心。以rpoC1基因为研究对象,在贝叶斯因子大于20的条件下,用HyPhy软件位点模型检测到3个正选择位点和541个负选择位点;用PAML软件位点模型检测到10个正选择位点,其中3个位点的后验概率超过99%。此外,基于最大似然法构建64种蕨类植物的系统发育树,结合HyPhy软件分析rpoC1基因的转换率、颠换率、转换率/颠换率、同义替换率、非同义替换率以及同义替换率/非同义替换率,探讨rpoC1基因内含子丢失与分子进化速率的关系。结果表明, rpoC1基因内含子缺失对转换率、颠换率以及非同义替换率有一定影响。     

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

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

蕨类植物rpoC1内含子缺失及其分子进化速率

rpoC1基因编码RNA聚合酶β°亚基蛋白, 在转录过程中与DNA模板结合, 与β亚基形成的β-β°亚基复合体构成RNA合成的催化中心。以rpoC1基因为研究对象, 在贝叶斯因子大于20的条件下, 用HyPhy软件位点模型检测到3个正选择位点和541个负选择位点; 用PAML软件位点模型检测到10个正选择位点, 其中3个位点的后验概率超过99%。此外, 基于最大似然法构建64种蕨类植物的系统发育树, 结合HyPhy软件分析rpoC1基因的转换率、颠换率、转换率/颠换率、同义替换率、非同义替换率以及同义替换率/非同义替换率, 探讨rpoC1基因内含子丢失与分子进化速率的关系。结果表明, rpoC1基因内含子缺失对转换率、颠换率以及非同义替换率有一定影响。     

哺乳动物cd59基因的进化

采用PCR以及BLAST方法从5种哺乳动物中获得了cd59基因的编码区序列, 结合 GenBank中已有的序列, 计算cd59基因在哺乳动物中的核苷酸替换速率. 对非同义替换速率和同义替换速率进行比较的结果显示, cd59在哺乳动物中总体上受到负选择作用; 用PAML软件“位点-特异”模型检测到4个受到正选择作用的位点, 4个位点分布于分子表面, 其中2个位于功能重要的区域; 此外, 用“支-位点-特异”模型在小鼠通过基因复制后形成的cd59a和cd59b上检测正选择引起的加速进化, 并检测到该支系特异的正选择位点1个.     

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

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

蕨类植物psaA基因的分子进化研究

采用“放松分子钟”模型、氨基酸位点正选择模型和分子内共进化网络估算方法,对蕨类植物光合系统Ⅰ核心蛋白PSAA编码基因psaA的进化趋势进行了研究。结果显示,叶绿体基因psaA编码区全序列具备成为蕨类植物系统发育关系重建位点的潜力,与rbcL基因联合后能构建高后验概率的系统发育树;蕨类植物的PSAA蛋白中存在一些曾经历正选择的氨基酸位点,其中29个位点聚合成为16个共进化组,通过共进化网络的方式协同影响光合系统Ⅰ的内部调整,提升其在被子植物兴起后光合环境下的适应能力。本文对蕨类植物进化潜能与分子机理的研究结果为揭示蕨类植物适应新生境提供了科学依据,也为植物系统分类学研究提供了分子依据。     

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

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

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

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

Comparative analysis of inverted repeats of polypod fern (Polypodiales) plastomes reveals two hypervariable regions

Background

Ferns are large and underexplored group of vascular plants (~ 11 thousands species). The genomic data available by now include low coverage nuclear genomes sequences and partial sequences of mitochondrial genomes for six species and several plastid genomes.

Results

We characterized plastid genomes of three species of Dryopteris, which is one of the largest fern genera, using sequencing of chloroplast DNA enriched samples and performed comparative analysis with available plastomes of Polypodiales, the most species-rich group of ferns. We also sequenced the plastome of Adianthum hispidulum (Pteridaceae). Unexpectedly, we found high variability in the IR region, including duplication of rrn16 in D. blanfordii, complete loss of trnI-GAU in D. filix-mas, its pseudogenization due to the loss of an exon in D. blanfordii. Analysis of previously reported plastomes of Polypodiales demonstrated that Woodwardia unigemmata and Lepisorus clathratus have unusual insertions in the IR region. The sequence of these inserted regions has high similarity to several LSC fragments of ferns outside of Polypodiales and to spacer between tRNA-CGA and tRNA-TTT genes of mitochondrial genome of Asplenium nidus. We suggest that this reflects the ancient DNA transfer from mitochondrial to plastid genome occurred in a common ancestor of ferns. We determined the marked conservation of gene content and relative evolution rate of genes and intergenic spacers in the IRs of Polypodiales. Faster evolution of the four intergenic regions had been demonstrated (trnA- orf42, rrn16-rps12, rps7-psbA and ycf2-trnN).

Conclusions

IRs of Polypodiales plastomes are dynamic, driven by such events as gene loss, duplication and putative lateral transfer from mitochondria.

     

The enigma of the gene coding for ribosomal protein S12 in the chloroplasts of Nicotiana.

A 2.9 kbp region from within the inverted repeat of Nicotiana chloroplast DNA hybridized with a chloroplast DNA fragment from Euglena containing the complete rps12 gene coding for ribosomal protein S12. Nucleotide sequencing within this region revealed the existance of two rps12 coding stretches interrupted by 540 bp having class II intron structure. Joining and decoding the exon regions produced a sequence of 85 amino acids colinear and 81% homologous to the S12 protein of Euglena chloroplasts and E. coli, starting from amino acid residue 38 to the stop codon. Immediately upstream of codon 38, conserved intron sequences were located. However, the 5' 37 codon of Nicotiana chloroplast rps12 could not be identified by electron microscopy of RNA-DNA hybrids within a DNA region extending 4000 bp upstream of codon 38, nor by computer search of a completely sequenced region extending for more than 9000 bp upstream of this codon. In E. coli, alteration in rps12 codons 42 or 87 causes streptomycin resistance. However, the nucleotide sequence of the identified rps12 exons in two Nicotiana chloroplast mutants resistant to streptomycin were found to be identical to that of wild type.     

The genes for the ribosomal proteins S12 and S7 are clustered with the gene for the EF-Tu protein on the chloroplast genome of Euglena gracilis.

We characterize a DNA segment of the Euglena gracilis chloroplast DNA fragment Eco . N by nucleotide sequencing and S1 nuclease analysis. We show that this region, which is upstream of the previously sequenced tuf A gene, contains the genes for the ribosomal proteins S12 and S7. The gene arrangement is 5'-rps 12-80 bp spacer-rps 7-174 bp spacer-tuf A, somewhat similar to the str operon of E. coli. The chloroplast S12 and S7 proteins contain 124 and 155 aminoacids, respectively, and are to 68% and 38% homologous with the corresponding E. coli proteins. The region is transcribed into a distronic mRNA of about 1.1 to 1.2 kb. The rps 12 and rps 7 genes, contrary to the tuf A gene, are not split.     

Relative rates of evolution in the coding and control regions of African mtDNAs

Reduced median networks of African haplogroup L mitochondrial DNA (mtDNA) sequences were analyzed to determine the pattern of substitutions in both the noncoding control and coding regions. In particular, we attempted to determine the causes of the previously reported (Howell et al. 2004) violation of the molecular clock during the evolution of these sequences. In the coding region, there was a significantly higher rate of substitution at synonymous sites than at nonsynonymous sites as well as in the tRNA and rRNA genes. This is further evidence for the operation of purifying selection during human mtDNA evolution. For most sites in the control region, the relative rate of substitution was similar to the rate of neutral evolution (assumed to be most closely approximated by the substitution rate at 4-fold degenerate sites). However, there are a number of mutational hot spots in the control region, approximately 3% of the total sites, that have a rate of substitution greater than the neutral rate, at some sites by more than an order of magnitude. It is possible either that these sites are evolving under conditions of positive selection or that the substitution rate at some sites in the control region is strongly dependent upon sequence context. Finally, we obtained preliminary evidence for "nonideal" evolution in the control region, including haplogroup-specific substitution patterns and a decoupling between relative rates of substitution in the control and coding regions.     

Ribosomal protein S14 transcripts are edited in Oenothera mitochondria

The gene encoding ribosomal protein S14 (rps14) in Oenothera mitochondria is located upstream of the cytochrome b gene (cob). Sequence analysis of independently derived cDNA clones covering the entire rps14 coding region shows two nucleotides edited from the genomic DNA to the mRNA derived sequences by C to U modifications. A third editing event occurs four nucleotides upstream of the AUG initiation codon and improves a potential ribosome binding site. A CGG codon specifying arginine in a position conserved in evolution between chloroplasts and E. coli as a UGG tryptophan codon is not edited in any of the cDNAs analysed. An inverted repeat 3' of an unidentified open reading frame is located upstream of the rps14 gene. The inverted repeat sequence is highly conserved at analogous regions in other Oenothera mitochondrial loci.     

Comparative genomics of four Liliales families inferred from the complete chloroplast genome sequence of Veratrum patulum O. Loes. (Melanthiaceae)

The sequence of the chloroplast genome, which is inherited maternally, contains useful information for many scientific fields such as plant systematics, biogeography and biotechnology because its characteristics are highly conserved among species. There is an increase in chloroplast genomes of angiosperms that have been sequenced in recent years. In this study, the nucleotide sequence of the chloroplast genome (cpDNA) of Veratrum patulum Loes. (Melanthiaceae, Liliales) was analyzed completely. The circular double-stranded DNA of 153,699 bp consists of two inverted repeat (IR) regions of 26,360 bp each, a large single copy of 83,372 bp, and a small single copy of 17,607 bp. This plastome contains 81 protein-coding genes, 30 distinct tRNA and four genes of rRNA. In addition, there are six hypothetical coding regions (ycf1, ycf2, ycf3, ycf4, ycf15 and ycf68) and two open reading frames (ORF42 and ORF56), which are also found in the chloroplast genomes of the other species. The gene orders and gene contents of the V. patulum plastid genome are similar to that of Smilax china, Lilium longiflorum and Alstroemeria aurea, members of the Smilacaceae, Liliaceae and Alstroemeriaceae (Liliales), respectively. However, the loss rps16 exon 2 in V. patulum results in the difference in the large single copy regions in comparison with other species. The base substitution rate is quite similar among genes of these species. Additionally, the base substitution rate of inverted repeat region was smaller than that of single copy regions in all observed species of Liliales. The IR regions were expanded to trnH_GUG in V. patulum, a part of rps19 in L. longiflorum and A. aurea, and whole sequence of rps19 in S. china. Furthermore, the IGS lengths of rbcL-accD-psaI region were variable among Liliales species, suggesting that this region might be a hotspot of indel events and the informative site for phylogenetic studies in Liliales. In general, the whole chloroplast genome of V. patulum, a potential medicinal plant, will contribute to research on the genetic applications of this genus.