Current approaches to whole genome phylogenetic analysis

It has long been known that evolutionary trees (phylogenies) can be estimated by comparing the DNA or protein sequences of homologous genes across different organisms. More recently, attempts have been made to estimate phylogenies by comparing entire genomes. These attempts have focused largely on comparisons of gene content and gene order. Many different methods have been proposed for making these comparisons. These include primarily maximum parsimony and distance methods, although more recently maximum likelihood and Bayesian methods are being developed. This paper discusses each of these approaches in turn, including their merits and limitations, and any software which is available to make use of them.     

Conservation of ribosomal protein gene ordering in 16 complete genomes

The organization of ribosomal proteins in 16 prokaryotic genomes was studied as an example of comparative genome analyses of gene systems. Hypothetical ribosomal protein-containing operons were constructed. These operons also contained putative genes and other non-ribosomal genes. The correspondences among these genes across different organisms were clarified by sequence homology computations. In this way a cross tabulation of 70 ribosomal proteins genes was constructed. On average, these were organized into 9-14 operons in each genome. There were also 25 non-ribosomal or putative genes in these mainly ribosomal protein operons. Hence the table contains 95 genes in total. It was found that: (i) the conservation of the block of about 20 r-proteins in the L3 and L4 operons across almost the entire eubacteria and archaebacteria is remarkable; (ii) some operons only belong to eubacteria or archaebacteria; (iii) although the ribosomal protein operons are highly conserved within domain, there are fine variations in some operons across different organisms within each domain, and these variations are informative on the evolutionary relations among the organisms. This method provides a new potential for studying the origin and evolution of old species.     

Conservation of ribosomal protein gene ordering in 16 complete genomes

The organization of ribosomal proteins in 16 prokaryotic genomes was studied as an example of comparative genome analyses of gene systems. Hypothetical ribosomal protein-containing operons were constructed. These operons also contained putative genes and other non-ribosomal genes. The correspondences among these genes across different organisms were clarified by sequence homology computations. In this way a cross tabulation of 70 ribosomal proteins genes was constructed. On average, these were organized into 9-14 operons in each genome. There were also 25 non-ribosomal or putative genes in these mainly ribosomal protein operons. Hence the table contains 95 genes in total. It was found that: (i) the conservation of the block of about 20 r-proteins in the L3 and L4 operons across almost the entire eubacteria and ar-chaebacteria is remarkable; (ii) some operons only belong to eubacteria or archaebacte-ria; (iii) although the ribosomal protein operons are highly conserved within domain, there are fine variat     

Computational modeling of gene structure in Arabidopsis thaliana

Computational gene identification by sequence inspection remains a challenging problem. For a typical Arabidopsis thaliana gene with five exons, at least one of the exons is expected to have at least one of its borders predicted incorrectly by ab initio gene finding programs. More detailed analysis for individual genomic loci can often resolve the uncertainty on the basis of EST evidence or similarity to potential protein homologues. Such methods are part of the routine annotation process. However, because the EST and protein databases are constantly growing, in many cases original annotation must be re-evaluated, extended, and corrected on the basis of the latest evidence. The Arabidopsis Genome Initiative is undertaking this task on the whole-genome scale via its participating genome centers. The current Arabidopsis genome annotation provides an excellent starting point for assessing the protein repertoire of a flowering plant. More accurate whole-genome annotation will require the combination of high-throughput and individual gene experimental approaches and computational methods. The purpose of this article is to discuss tools available to an individual researcher to evaluate gene structure prediction for a particular locus.     

Quantitative determination of gene strand bias in prokaryotic genomes

Comparative genometrics of microorganisms is a relatively new area, in which genome properties are translated into numerical indexes. Such indexes can be used for a comprehensive and comparative analysis of microbial genomes, contributing to the understanding of their evolution. This work presents a new method for quantitative determination of gene strand bias in prokaryotic chromosomes, in which data transformation of gene position skew leads to a numerical index that can be applied to quantitative comparisons of genome organization. It was applied in the comparative analysis of 49 completely sequenced Firmicutes genomes, allowing the distinction of groups defined according to their patterns of gene strand preference. The resulting groups revealed that, regarding gene strand bias, reduced genomes are, in general, the more disordered among Firmicutes, while genomes of extremophile organisms comprehend those with the highest degree of genome organization in this phylum.     

基于转录组解析淡玫红鹅膏环肽基因家族特征及其系统发育

【背景】淡玫红鹅膏(Amanita pallidorosea)是鹅膏属檐托鹅膏组的一种剧毒鹅膏菌,其子实体内含有丰富的鹅膏环肽毒素,但其编码毒环肽和相关肽的基因家族"MSDIN"尚待深入研究。【目的】探究淡玫红鹅膏中编码毒环肽及相关肽的基因家族成员的多样性、保守性及其系统发育。【方法】采用Illumina HiSeq2000平台对淡玫红鹅膏转录组进行测序,使用TBLASTn软件对MSDIN基因家族进行检索,并设计特异性引物进行PCR验证,通过生物信息学比对分析MSDIN基因家族成员的种类和序列构成,通过重建分子系统发育树了解其演化历程。【结果】通过对MSDIN基因家族的查找,从转录组数据获得了60条环肽编码基因,经PCR验证其可编码32条环肽,包括α-鹅膏毒肽、β-鹅膏毒肽和羧基二羟鬼笔毒肽。本研究报道了8条新环肽。分子系统发育树分析显示,鹅膏环肽分为鹅膏毒肽、鬼笔毒肽和未知功能环肽3大支。系统发育结合保守序列推测出7条潜在的新毒肽。【结论】淡玫红鹅膏具有丰富的鹅膏环肽资源,利用转录组测序能够系统挖掘鹅膏环肽新资源,为其整体结构解析奠定基础。     

Identification and categorization of horizontally transferred genes in prokaryotic genomes

Horizontal gene transfer （HGT）, a process through which genomes acquire genetic materials from distantly related organisms, is believed to be one of the major forces in prokaryotic genome evolution.However, systematic investigation is still scarce to clarify two basic issues about HGT： （1） what types of genes are transferred; and （2） what influence HGT events over the organization and evolution of biological pathways. Genome-scale investigations of these two issues will advance the systematical understanding of HGT in the context of prokaryotic genome evolution. Having investigated 82 genomes, we constructed an HGT database across broad evolutionary timescales. We identified four function categories containing a high proportion of horizontally transferred genes： cell envelope, energy metabolism, regulatory functions, and transport/binding proteins. Such biased function distribution indicates that HGT is not completely random;instead, it is under high selective pressure, required by function restraints in organisms. Furthermore, we mapped the transferred genes onto the connectivity structure map of organism-specific pathways listed in Kyoto Encyclopedia of Genes and Genomes （KEGG）. Our results suggest that recruitment of transferred genes into pathways is also selectively constrained because of the tuned interaction between original pathway members. Pathway organization structures still conserve well through evolution even with the recruitment of horizontally transferred genes. Interestingly, in pathways whose organization were significantly affected by HGT events, the operon-like arrangement of transferred genes was found to be prevalent. Such results suggest that operon plays an essential and directional role in the integration of alien genes into pathways.     

Identification of horizontal gene transfer and recombination of PsaA gene in streptococcus mitis group

Pneumococcal surface adhesin A (psaA) gene is universally confirmed as one of the Streptococcus pneumoniae adhesion genes, but it is disputed whether the psaA gene is a Streptococcus pneumoniae species‐specific gene. In the present study, the presence of the psaA gene in 34 streptococcus mitis group isolates was identified by the PCR approach and a comparison of sequencing PCR products (Streptococcus pneumoniae R6 as the control strain). Also, the evolutionary scenarios of these psaA genes in these streptococcus mitis group isolates were analyzed by a phylogenetic tree based on the housekeeping genes (sodA and rnpB) and psaA genes. As a result, a high degree of conservation of open reading frame sequences in all six Streptococcus pneumoniae strains (100% similarity) and in the other species of the streptococcus mitis group (92.6–100% similarity) was revealed. Further genetics research based on housekeeping genes and psaA gene phylogenies showed that the psaA gene was of vertical inheritance only in Streptococcus pneumoniae; however, high‐frequency horizontal psaA gene transfer and recombination occurred in the other species of the streptococcus mitis group. These findings confirmed that the psaA gene was not a Streptococcus pneumoniae species‐specific gene, and high‐frequency HGT and recombination events may explain the presence of the psaA gene in the other species of the streptococcus mitis group.     

用matK序列分析探讨木兰属植物的系统发育关系

用木兰科Magnoliaceae 57种植物的matK基因序列构建了该科的系统发育分支图。结果表明: (1)木兰属Magnolia L.是一个因为性状的趋同演化而建立的多系类群; (2)木兰亚属subgen. Magnolia和玉兰亚属subgen. Yulania (Spach) Reichenb.亲缘关系较远, 支持将后者从该属中分出建立玉兰属Yulania Spach, 木兰亚属作为木兰属保留; (3)木兰亚属的sect. Splendentes Dandy ex Vazquez组与皱种组sect. Rytidospermum Spach的两个美洲种M. macrophylla Michaux和M. dealbata Zucc.亲缘关系较近, 荷花玉兰组sect. Theorhodon Spach与常绿组sect. Gwillimia DC.的亲缘关系较近; (4)盖裂木属Talauma Juss.可以成立, 而其分布于亚洲的Blumiana Blume组可归入木兰属; (5)拟单性木兰属Parakmeria Hu &; Cheng、华盖木属Manglietiastrum Law以及单性木兰属Kmeria (Pierre) Dandy形成一个单系群, 与玉兰亚属和含笑属Michelia L.的亲缘关系较近。花的着生位置不足以作为木兰科的分族依据, 含笑族Michelieae和木兰族Magnolieae的特征及其界定应做修改。将玉兰亚属从木兰属分出后, 木兰属与含笑属无性状交叉,成为两个区别明显的属。     

Periodicity in prokaryotic and eukaryotic genomes identified by power spectrum analysis

We used a power spectrum method to identify periodic patterns in nucleotide sequence, and characterized nucleotide sequences that confer periodicities to prokaryotic and eukaryotic genomes and genomes. A 10-bp periodicity was prevalent in hyperthermophilic bacteria and archaebacteria, and an 11-bp periodicity was prevalent in eubacteria. The 10-bp periodicity was also prevalent in the eukaryotes such as the worm Caenorhabditis elegans. Additionally, in the worm genome, a 68-bp periodicity in chromosome I, a 59-bp periodicity in chromosome II, and a 94-bp periodicity in chromosome III were found. In human chromosomes 21 and 22, approximately 167- or 84-bp periodicity was detected along the entire length of these chromosomes. Because the 167-bp is identical to the length of DNA that forms two complete helical turns in nucleosome organization, we speculated that the respective sequences may correspond to arrays of a special compact form of nucleosomes clustered in specific regions of the human chromosomes. This periodic element contained a high frequency of TGG. TGG-rich sequences are known to form a specific subset of folded DNA structures, and therefore, the sequences might have potential to form specific higher order structures related to the clustered occurrence of a specific form of the speculated nucleosomes.     

Phylogenetic analysis of coat protein gene of BYDV-MAV strain from wheat

Barley yellow dwarf disease is globally the most important viral disease of wheat. The full-length nucleotide sequence of coat protein (CP) gene of 12 isolates revealed the presence of three distinct clusters. Pakistani isolate of MAV (MAV-PK) has maximum similarity of 99.23% with MAV isolate of Morocco and PAV-Australia following 99.22 and 99.22% with PAV-France. Similar degree of similarity was found in comparison of amino acid sequence. The finding of this study is that MAV-PK has similarity with both MAV-France and PAV-Australia, which is due to the reason that both MAV and PAV belong to the same group and both share maximum nucleotide homology. Low genetic diversity was found not only between MAV isolates but also between MAV and PAV isolates because phylogenetic analysis was done on the CP gene which is highly conserved region in genome of Barley yellow dwarf viruses (BYDVs). Divergence in MAV-PK was due to this recombination which is now most prevalent in Pakistan. MAV-PK has maximum similarity with MAV-Morocco followed by MAV-Sweden and MAV-Cz, which seems to indicate that Pakistani isolate of MAV evolved as the result of recombination between MAV isolates of the USA and PAV isolates of Australia and France. At the same time, recombination of MAV-CZ and MAV-Sweden also occur. This work can be successfully utilised in epidemiological studies of MAV isolate in Pakistan. Further analysis of variation level in these isolates will help scientists to formulate appropriate management strategies like incorporation of BdV 2 gene in wheat against BYDVs.     

金针菇CA基因家族分析及其对CO2的响应

金针菇栽培通过提高栽培室CO₂浓度来抑制菌盖生长,提高商品质量。碳酸酐酶能调节CO₂/HCO₃ ^-在细胞内的平衡,它可能是响应CO₂胁迫的关键酶。本研究利用我们已完成的金针菇Flammulina filiformis的基因组测序数据,以双色蜡蘑碳酸酐酶家族基因为参照,通过本地BLAST,得到7个金针菇碳酸酐酶家族的基因,分别命名为CA-1、CA-2、CA-3、CA-4、CA-5、CA-6、CA-7。基因结构、保守结构域和系统进化分析结果均支持这7个序列是碳酸酐酶家族的编码基因。高浓度CO₂胁迫处理金针菇子实体12h,CA-1、CA-2表达量与CO₂浓度正相关,CA-5负相关,CA-3、CA-4和CA-6没有响应,CA-7无明显规律。据此推测,CA-1、CA-2和CA-5可能是金针菇响应CO₂胁迫的关键基因之一,参与调控菌盖生长发育。     

Detection of operons

Operons are clusters of genes that are transcribed as a single message, and regulated by the same gene expression machinery. They are found primarily in prokaryotic genomes. Because genes in the same operon are likely to have related functions, identification of the operon structure is potentially useful for assigning gene function. We report the development and benchmarking of two different methods for detecting operons, based on an analysis of 42 fully sequenced prokaryotic organisms. The Gene Neighbor method (GNM) utilizes the relatively high conservation of gene order in operons, compared with genes in general. The Gene Gap Method (GGM) makes use of the relatively short gap between genes in operons compared with that otherwise found between adjacent genes. The methods have been benchmarked using KEGG pathway data and RegulonDB Escherichia coli operon data. With optimum parameters, the specificity of the GNM is 93% and the sensitivity is 70%. For the GGM, the specificity is 95% and the sensitivity is 68%. Together, the two methods have a sensitivity of 87.2%, while joint predictions have a sensitivity of 50% and a specificity of 98%. The methods are used to infer possible functions for some hypothetical genes in prokaryotic genomes. The methods have proven a useful addition to structure information in deriving protein function in a structural genomics project.     

The evolutionary origin of Xanthomonadales genomes and the nature of the horizontal gene transfer process

Determining the influence of horizontal gene transfer (HGT) on phylogenomic analyses and the retrieval of a tree of life is relevant for our understanding of microbial genome evolution. It is particularly difficult to differentiate between phylogenetic incongruence due to noise and that resulting from HGT. We have performed a large-scale, detailed evolutionary analysis of the different phylogenetic signals present in the genomes of Xanthomonadales, a group of Proteobacteria. We show that the presence of phylogenetic noise is not an obstacle to infer past and present HGTs during their evolution. The scenario derived from this analysis and other recently published reports reflect the confounding effects on bacterial phylogenomics of past and present HGT. Although transfers between closely related species are difficult to detect in genome-scale phylogenetic analyses, past transfers to the ancestor of extant groups appear as conflicting signals that occasionally might make impossible to determine the evolutionary origin of the whole genome.     

副猪嗜血杆菌aroA基因鉴定及遗传进化分析

[目的]细菌aroA基因参与芳香族氨基酸的生物合成,被成功应用于细菌分类和基因失活致弱突变菌株的构建.副猪嗜血杆菌(Hps)是感染猪出现多发性浆膜炎和关节炎的一种病原细菌,鉴定该菌aroA全基因序列将有助于鉴定遗传进化关系和突变分析.[方法]利用PCR和细菌基因组步移技术鉴定Hps的aroA基因序列,进而对不同血清型菌株该基因序列进行鉴定,并与其它革兰氏阴性细菌进行比对和遗传进化分析.[结果]自Hps血清5型基因组DNA中获得包含完整aroA基因的3.7 kb基因片段,其中aroA基因全长1314 bp,编码产物长度437 aa,分子量大小47.9 kDa,该基因上游为磷酸烯醇式丙酮酸羧化酶基因.自本试验选择的Hps不同血清型菌株中均可扩增出包含完整aroA基因的1476 bp片段,且这些不同血清型菌株间核酸序列同源性在97.7%以上.Hps血清5型aroA基因序列与巴氏杆菌科其它成员核酸序列同源性为70.6%-78.9%,与E.coli和S.typhi-murium的同源性分别为66.4%和67.2%.[结论]本试验首次对Hps的15个血清型国际参考菌株及地方分离株aroA全基因序列进行了鉴定,序列比较结果显示aroA基因在革兰氏阴性细菌中具有较高的同源性.aroA基因鉴定对构建基因失活突变菌株以研究Hps生物学特性奠定了基础.     

Comparative analysis of two-component signal transduction system in two streptomycete genomes

Species of the genus Streptomyces are major bacteria responsible for producing most natural antibiotics. Streptomyces coelicolor A3(2) and Streptomyces avermitilis were sequenced in 2002 and 2003, respectively. Two-component signal transduction systems (TCSs), consisting of a histidine sensor kinase (SK) and a cognate response regulator (RR), form the most common mechanism of transmembrane signal transduction in prokaryotes. TCSs in S. coelicolor A3(2) have been analyzed in detail. Here, we identify and classify the SK and RR of S. avermitilis and compare the TCSs with those of S. coelicolor A3(2) by computational approaches. Phylogenetic analysis of the cognate SK-RR pairs of the two species indicated that the cognate SK-RR pairs fall into four classes according to the distribution of their orthologs in other organisms. In addition to the cognate SK-RR pairs, some potential partners of non-cognate SK-RR were found, including those of unpaired SK and orphan RR and the cross-talk between different components in either strain. Our study provides new clues for further exploration of the molecular regulation mechanism of streptomycetes with industrial importance.     

人类基因组中巢式基因对的系统分析

基因组上一个基因的所有或大部分外显子位于另一基因内含子和UTR中被称为巢式基因(Nested gene)。巢式基因对(Nested gene pair)是由主基因(Host gene)和巢式基因组成的基因对。文章针对人类基因组上的顺式巢式基因对(简称为巢式基因对)进行全基因组鉴定,并通过比较巢式基因对在人和小鼠基因组上的保守性分布,研究巢式基因对的进化方式。保守性分析表明基因转座、序列原位突变和主基因转录起始/终止位点突变是巢式基因对进化的主要原因,其中主基因转录起始/终止位点突变是巢式基因对独特的一种进化方式。Gene On-tology分析揭示大部分巢式基因与主基因的产物在功能上没有相关性。     

真核生物转座子鉴定和分类计算方法

重复序列是真核生物基因组的重要组成成分,根据其序列特征及在基因组中的存在形式,可以进一步分为串联重复、片段重复和散在重复。其中,散在重复大多起源于转座子。根据转座介质的不同,转座子又可分为DNA和逆转录转座子。转座子的转座和扩增对基因的进化和基因组的稳定具有显著的影响;同时与其他类型的重复序列相比,转座子的结构和分类更为复杂多样,使得对转座子的鉴定和分类更为复杂和困难。鉴于此,文章简要概括了转座子的功能及分类,总结了真核生物转座子鉴定、分类和注释的3个步骤:(1)重复序列库的构建;(2)重复序列的校正和分类;(3)基因组注释。着重介绍了每一步骤所采用的不同计算方法,比较了不同方法的优缺点。只有把多种方法结合起来使用才能实现全基因组转座子的精确鉴定、分类和注释,这将为转座子的全基因组鉴定和分类提供借鉴意义。