转录组测序揭示翼盖蕨(Didymochlaena trancatula)的全基因组复制历史

全基因组复制在动植物中普遍存在, 被认为是促进物种进化的重要动力之一。作为蕨类植物的单种科物种, 翼盖蕨(Didymochlaena trancatula)是真水龙骨类I的基部类群, 在蕨类中具有独特的演化地位。本研究基于高通量测序, 通过同义替换率(Ks)分析、相对定年分析揭示翼盖蕨的全基因组复制发生情况。Ks分析表明, 翼盖蕨至少经历了两次全基因组复制事件, 其中一次发生于59-62 million years ago (Mya), 另一次发生于90-94 Mya, 这两次全基因组复制事件分别和白垩纪第三纪的Cretaceous-Tertiary (C-T)大灭绝事件以及翼盖蕨的物种分化时间相吻合。进一步对两次全基因组复制保留的基因进行功能注释和富集分析, 结果显示与转录及代谢调控相关的基因优势被保留。翼盖蕨的全基因组复制事件可能促进了该物种的分化及其对极端环境的适应性。     

全基因组复制事件的绝对定年揭示莲座蕨属植物的迟滞演化

全基因组复制在维管植物的物种形成过程中普遍存在, 被认为是物种适应极端环境的重要机制之一。确定全基因组复制事件的发生时间对理解生物的适应性演化具有重要意义。然而, 在维管植物, 特别是蕨类植物中, 全基因组复制事件的发生时间及其演化意义仍知之甚少。本研究以蕨类植物重要基部类群——福建莲座蕨(Angiopteris fokiensis)为例, 基于不同采样点(广东、广西、上海)的3个转录组学数据, 利用同义替换率(Ks)和绝对定年的方法分析全基因组复制事件的发生时间和物种单位时间内的分子演化速率, 并对事件发生后保留下的基因进行基因功能注释和富集分析。结果表明, 福建莲座蕨在159‒165 Mya发生了一次全基因组复制事件, 该复制事件优先保留的基因主要与营养代谢、信号传导、适应调节和组织结构生长相关。另外, 福建莲座蕨的分子演化速率为1.66 × 10^‒9 (同义替换/位点/年), 是除裸子植物外, 陆生植物中已知演化速率最缓慢的类群。综合以上研究结果, 我们推测福建莲座蕨全基因组复制的发生可能与裸子植物繁盛、核心被子植物集中兴起或托阿尔阶灭绝事件有关。而复制后显著保留基因可能促进了莲座蕨属(Angiopteris)植物的遗传和形态创新, 从而帮助其快速适应环境的剧烈变化。进一步对该类群植物演化速率缓慢的原因进行讨论, 推测莲座蕨属缓慢的演化速率可能与其本身世代周期长、基因组较大及其生长环境稳定有关。本研究通过分析福建莲座蕨的全基因组复制历史和复制基因的保留模式, 推测全基因组复制事件对促进演化速率较慢的植物适应极端环境变化具有重要意义, 可为理解其他陆生植物的适应性演化提供更多启发。     

多叶斑叶兰繁殖体系建立及基于转录组的发育调控途径功能基因研究

多叶斑叶兰,是兰科斑叶兰属濒危野生植物、国家二级保护植物,具有极高的观赏和药用价值。多叶斑叶兰由于分布种群小,传播扩散能力弱,自然繁殖受到极大限制。以野生多叶斑叶兰茎段作为外植体,建立高效直接的植株再生体系。结合高通量转录组测序技术与生物信息学分析技术,深入挖掘参与多叶斑叶兰器官发育过程的功能基因。结果表明, Morel+2.0 mg/L 6-BA +0.5 mg/LKT+1.0 mg/LNAA+1g/L蛋白胨+25g/L 蔗糖+7.0 g/L Agar+1.0 g/L活性炭+30 g/L 香蕉+50 g/L土豆为最佳的芽诱导培养基;Morel+3 mg/L 6-BA+0.5mg/L NAA+0.5 mg/LKT+0.01 mg/L TDZ+2g/L蛋白胨+25g/L 蔗糖+7.0 g/L Agar+1.0 g/L活性炭30g/L 香蕉+50g/L 土豆为最佳芽增殖培养基;1/2 Morel+1.0 mg/L IBA+0.1 mg/L NAA+1 g/L+花宝2 号+25g/L 蔗糖+7.0g/L Agar +1.0 g/L活性炭+1g·L-1蛋白胨为最佳的生根培养基。转录组测序分析组装获得170688个Unigene,平均长度为584 bp, N50为833 bp。17352个Unigene比对注释到NR、Swiss-Prot、KOG、GO、KEGG数据库。野生苗与组培苗差异表达Unigene GO与KEGG功能富集分析显示, Unigene主要参与调控植物激素信号转导、植物形态发育、次生代谢过程以及能量代谢过程等生物学功能。进一步分析获得511个参与植物器官发育调控相关的转录因子。成功建立了多叶斑叶兰种质资源保存与高效离体再生体系,结合高通量转录组学技术,获得全面完整的多叶斑叶兰转录组信息特征,为后期多叶斑叶兰快速扩繁、遗传转化以及功能基因鉴定、遗传发育及其调控机制研究奠定基础。     

Investigating ancient duplication events in the Arabidopsis genome

The complete genomic analysis of Arabidopsis thaliana has shown that a major fraction of the genome consists of paralogous genes that probably originated through one or more ancient large-scale gene or genome duplication events. However, the number and timing of these duplications still remains unclear, and several different hypotheses have been put forward recently. Here, we reanalyzed duplicated blocks found in the Arabidopsis genome described previously and determined their date of divergence based on silent substitution estimations between the paralogous genes and, where possible, by phylogenetic reconstruction. We show that methods based on averaging protein distances of heterogeneous classes of duplicated genes lead to unreliable conclusions and that a large fraction of blocks duplicated much more recently than assumed previously. We found clear evidence for one large-scale gene or even complete genome duplication event somewhere between 70 to 90 million years ago. Traces pointing to a much older (probably more than 200 million years) large-scale gene duplication event could be detected. However, for now it is impossible to conclude whether these old duplicates are the result of one or more large-scale gene duplication events. abbreviations dA, fraction of amino acid substitutions; Kn, number of nonsynonymous substitutions per nonsynonymous site; Ks, number of synonymous substitutions per synonymous site; MYA, million years ago     

Variation in the pattern of synonymous and nonsynonymous difference between two fungal genomes

The proportion of synonymous nucleotide differences per synonymous site (p(S)) and the proportion of nonsynonymous differences per nonsynonymous site (p(N)) were computed at 1,993,217 individual codons in 4,133 protein-coding genes between the two yeast species Saccharomyces cerevisiae and Saccharomyces paradoxus. When the modified Nei-Gojobori method was used, significantly more codons with p(N) > p(S) were observed than expected, based on random pairing of observed p(S) and p(N) values. However, this finding was most likely explained by the presence of a strong negative correlation between the number of synonymous differences and the number of nonsynonymous differences at codons with at least one difference. As a result of this correlation, codons with p(N) > p(S) were characterized not only by unusually high p(N) but also by unusually low p(S). On the other hand, the number of codons with p(N)>p(S) (where p(S) is the mean p(S) for all codons) was very similar to the random expectation, and the observed number of 30-codon windows with p(N) > p(S) was significantly lower than the random expectation. These results imply that the occurrence of a certain number of codons or codon windows with p(N) > p(S) is expected given the nature of nucleotide substitution and need not imply the action of positive Darwinian selection.     

Site-to-site variation of synonymous substitution rates

We develop a new model for studying the molecular evolution of protein-coding DNA sequences. In contrast to existing models, we incorporate the potential for site-to-site heterogeneity of both synonymous and nonsynonymous substitution rates. We demonstrate that within-gene heterogeneity of synonymous substitution rates appears to be common. Using the new family of models, we investigate the utility of a variety of new statistical inference procedures, and we pay particular attention to issues surrounding the detection of sites undergoing positive selection. We discuss how failure to model synonymous rate variation in the model can lead to misidentification of sites as positively selected.     

Evolution of BK Virus Based on Complete Genome Data

The human polyomavirus BK virus (BKV) is ubiquitous in humans, infecting children asymptomatically. BKV is the only primate polyomavirus that has subtypes (I–IV) distinguishable by immunological reactivity. Nucleotide (nt) variations in a major capsid protein (VP1) gene region (designated the epitope region), probably responsible for antigenic diversity, have been used to classify BKV isolates into subtypes. Here, with all the protein-encoding gene sequences, we attempted to elucidate the evolutionary relationships among 28 BKV isolates belonging to subtypes I, III, and IV (no isolate belonging to subtype II, a minor one, was included). First, using the GTR + Γ + I model, maximum likelihood trees were reconstructed for individual viral genes as well as for concatenated viral genes. On the resultant trees, the 28 BKV isolates were consistently divided into three clades corresponding to subtypes I, III, and IV, although bootstrap probabilities are not always high. Then we used more sophisticated likelihood models, one of which takes account of codon structure, to elucidate the phylogenetic relationships among BKV subtypes, but the phylogeny of the deep branchings remained ambiguous. Furthermore, the possibility of positive selection in the evolution of BKV was examined using the nonsynonymous/synonymous rate ratio as a measure of selection. An analysis based on entire genes could not detect any strong evidence for positive selection, but that based on the epitope region identified a few sites potentially under positive selection (these sites were among those showing subtype linked polymorphisms). These author Yuriko Nishimoto and Tomokazu Takasaka contributed equally to this article. [Reviewing Editor: Dr. Rasmus Nielsen]     

On the varied pattern of evolution of 2 fungal genomes: a critique of Hughes and Friedman

A number of statistical tests have been proposed to detect positive Darwinian selection affecting a few amino acid sites in a protein, exemplified by an excess of nonsynonymous nucleotide substitutions. These tests are often more powerful than pairwise sequence comparison, which averages synonymous (d(S)) and nonsynonymous (d(N)) rates over the whole gene. In a recent study, however, Hughes AL and Friedman R (2005. Variation in the pattern of synonymous and nonsynonymous difference between two fungal genomes. Mol Bio Evol. 22: 1320-1324) argue that d(S) and d(N) are expected to fluctuate along the sequence by chance and that an excess of nonsynonymous differences in individual codons is no evidence for positive selection. The authors compared codons in protein-coding genes from the genomes of 2 yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus. They calculated the proportions of synonymous and nonsynonymous differences per site (p(S) and p(N)) in every codon and discovered that p(N) is often greater than p(S) and that among some codons p(S) and p(N) are negatively correlated. The authors argued that these results invalidate previous tests of codons under positive selection. Here I discuss several errors of statistics in the analysis of Hughes and Friedman, including confusion of statistics with parameters, arbitrary data filtering, and derivation of hypotheses from data. I also apply likelihood ratio tests of positive selection to the yeast data and illustrate empirically that Hughes and Friedman's criticisms on such tests are not valid.     

A large variation in the rates of synonymous substitution for RNA viruses and its relationship to a diversity of viral infection and transmission modes

RNA viruses successfully adapt to various environments by repeatedly producing new mutants, often through generating a number of nucleotide substitutions. To estimate the degree of variation in mutation rates of RNA viruses and to understand the source of such variation, we studied the synonymous substitution rate because synonymous substitution is exempt from functional constraints at the protein level, and its rate reflects the mutation rate to a great extent. We estimated the synonymous substitution rates for a total of 49 different species of RNA viruses, and we found that the rates had tremendous variation by 5 orders of magnitude (from 1.3 x 10(-7) to 6.2 x 10(-2) /synonymous site/year). Comparing the synonymous substitution rates with the replication frequencies and replication error rates for the RNA viruses, we found that the main source of the rate variation was differences in the replication frequency because the rates of replication error were roughly constant over different RNA viruses. Moreover, we examined a relationship between viral life strategies and synonymous substitution rates to understand which viral life strategies affect replication frequencies. The results show that the variation of synonymous substitution rates has been influenced most by either the difference in the infection modes or the differences in the transmission modes. In conclusion, the variation of mutation rates for RNA viruses is caused by different replication frequencies, which are affected strongly by the infection and transmission modes.     

氨基酸的置换数与生物大分子进化改变量的测度

以蛋白质分子的氨基酸置换数或核酸分子的核苷酸置换数为衡量尺度 ,说明生物大分子随时间的改变 (即分子进化速率 )保持相对恒定     

用RAPD分析四种无尾类的系统演化关系

利用RAPD技术检测了分属无尾目3个科（雨蛙科、蟾蜍科、蛙科）的黑眶蟾蜍（Bufomelanostictus)、中国雨蛙（Hyla chinensis）、泽陆蛙（Rana limnocharis）、沼水蛙（R.guentheri)的系统发生关系。经19个随机引物对4个物种基因组DNA进行扩增,选择其中扩增谱带清晰的16个引物进行分析,计算不同科间及同一科内不同种间的遗传距离,结果表明：16个引物获得的RAPD谱带均表现出不同程度的多态性;泽陆蛙与沼水蛙间的亲缘关系最近,而黑眶蟾蜍与中国雨蛙之间的亲缘关系较黑眶蟾蜍与蛙科的泽陆蛙、沼水蛙之间以及中国雨蛙与泽陆蛙、沼水蛙之间的亲缘关系近,从基因组DNA水平上也说明雨蛙科与蟾蜍科间的亲缘关系更近,与蛙科的亲缘关系更远,这与形态学、染色体和线粒体DNA多态性研究的分析结果一致,从而进一步从分子水平上为无属目这3科的系统演化提供了新的证据。