Construction and comparative analysis of two-component system and metabolic network profile based phylogenetic trees

In this study, two-component system (TCS) gene profile and metabolic network gene profile based phylogenetic trees were constructed and compared to each other to evaluate the evolutionary relationship between the bacterial sensing system and metabolism. The gene profiles of the these systems suggested that bacteria employed different evolutionary strategies to optimize the two-component system and metabolic network. In addition, comparative analysis revealed that the TCS based tree showed better family grouping than the metabolic network based tree, which indicated that the TCS and metabolic network have been modified via self-evolution and recruitment methods, respectively.     

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.     

Whole genome-based phylogenetic analysis of free-living microorganisms.

A phylogenetic 'tree of life' has been constructed based on the observed presence and absence of families of protein-encoding genes observed in 11 complete genomes of free-living microorganisms. Past attempts to reconstruct the evolutionary relation-ships of microorganisms have been limited to sets of genes rather than complete genomes. Despite apparent rampant lateral gene transfer among microorganisms, these results indicate a single robust underlying evolutionary history for these organisms. Broadly, the tree produced is very similar to the small subunit rRNA tree although several additional phylogenetic relationships appear to be resolved, including the relationship of Archaeoglobus to the methanogens studied. This result is in contrast to notions that a robust phylogenetic reconstruction of microorganisms is impossible due to their genomes being composed of an incomprehensible amalgam of genes with complicated histories and suggests that this style of genome-wide phylogenetic analysis could become an important method for studying the ancient diversification of life on Earth. Analyses using informational and operational subsets of the genes showed that this 'tree of life' is not dependent on the phylogenetically more consistent informational genes.     

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.     

Chloroplast genome characteristics and phylogenetic analysis of the medicinal plant Blumea balsamifera (L.) DC

Blumea balsamifera (L.) DC., a medicinal plant with high economic value in the Asteraceae family, is widely distributed in China and Southeast Asia. However, studies on the population structure or phylogenetic relationships with other related species are rare owing to the lack of genome information. In this study, through high-throughput sequencing, we found that the chloroplast genome of B. balsamifera was 151,170 bp in length, with a pair of inverted repeat regions (IRa and IRb) comprising 24,982 bp, a large single-copy (LSC) region comprising 82,740 bp, and a small single-copy (SSC) region comprising 18,466 bp. A total of 130 genes were identified in the chloroplast genome of B. balsamifera, including 85 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes; furthermore, sequence analysis identified 53 simple sequence repeats. Whole chloroplast genome comparison indicated that the inverted regions (IR) were more conserved than large single-copy and SSC regions. Phylogenetic analysis showed that B. balsamifera is closely related to Pluchea indica. Conclusively, the chloroplast genome of B. balsamifera was helpful for species identification and analysis of the genetic diversity and evolution in the genus Blumea and family Asteraceae.     

Whole-genome survey and phylogenetic analysis of Gadus macrocephalus

Gadus macrocephalus (Pacific cod) is an economically important species on the northern coast of the Pacific. Although numerous studies on G. macrocephalus exist, there are few reports on its genomic data. Here, we used whole-genome sequencing data to elucidate the genomic characteristics and phylogenetic relationship of G. macrocephalus. From the 19-mer frequency distribution, the genome size was estimated to be 658.22 Mb. The heterozygosity, repetitive sequence content and GC content were approximately 0.62%, 27.50% and 44.73%, respectively. The draft genome sequences were initially assembled, yielding a total of 500,760 scaffolds (N50 = 3565 bp). A total of 789,860 microsatellite motifs were identified from the genomic data, and dinucleotide repeat was the most dominant simple sequence repeat motif. As a byproduct of whole-genome sequencing, the mitochondrial genome was assembled to investigate the evolutionary relationships between G. macrocephalus and its relatives. On the basis of 13 protein-coding gene sequences of the mitochondrial genome of Gadidae species, the maximum likelihood phylogenetic tree showed that complicated relationships and divergence times among Gadidae species. Demographic history analysis revealed changes in the G. macrocephalus population during the Pleistocene by using the pairwise sequentially Markovian coalescent model. These findings supplement the genomic data of G. macrocephalus, and make a valuable contribution to the whole-genome studies on G. macrocephalus.     

鸡IL-18成熟蛋白基因的克隆及分子进化分析

根据已发表的鸡白介素18(IL-18)基因序列设计合成引物,以植物凝集素(PHA)和脂多糖(LPS)激活的AA肉鸡脾细胞mRNA为模板,通过RT-PCR扩增出编码鸡IL-18成熟蛋白的eDNA。将该eDNA克隆于pUCm-T载体,并对其进行测序,结果表明所克隆的核苷酸片段包含了全部成熟蛋白编码基因,成熟蛋白编码区507个核苷酸,编码169个氧基酸。把该基因编码的鸡成熟IL-18蛋白氨基酸序列与已公布的禽及哺乳动物IL-8成熟蛋白基因氧基酸序列进行比较,其同源性分别在96．5％～100％和20．1％～26．6％之间,分子系统进化树分析表明鸡IL-18与哺乳动物IL-18有共同的祖先,亲源关系较近,但在免疫系统选择性压力下,形成独特的种族特异性。鸡IL-18基因的克隆为体外表达鸡IL-18蛋白及作为免疫佐剂应用于预防接种的研究奠定了基础。     

新型水果空心泡叶绿体基因组特征及其系统发育分析北大核心CSCD

为探究空心泡(Rubus rosaefolius)叶绿体基因组特征,本研究以空心泡为试验材料,采用Illumina NovaSeq平台进行高通量测序,获得空心泡完整的叶绿体基因组序列,并进行空心泡叶绿体基因序列特征和系统发育分析。结果表明:空心泡的完整叶绿体基因组总长度为155650 bp,具有典型的四分体结构,包括2个反向重复序列(各25748 bp)、1个大拷贝区(85443 bp)、1个小拷贝区(18711 bp)。空心泡叶绿体全基因组共鉴定出131个基因,包括86个蛋白质编码基因、37个tRNA基因和8个rRNA基因,全基因组的GC含量为36.9%。空心泡叶绿体基因组包含47个散在重复序列、72个简单重复序列(simple sequence repeating,SSR)位点,密码子偏好性为亮氨酸密码子,偏好使用A/U结尾的密码子。系统发育分析表明,空心泡与小叶悬钩子(Rubus taiwanicola)亲缘关系最近,其次是能高悬钩子(Rubus rubroangustifolius)和腺萼悬钩子(Rubus glandulosopunctatus)。空心泡的叶绿体基因组特征及其系统发育分析,为空心泡的遗传多样性研究和叶绿体开发利用提供理论依据。     

副猪嗜血杆菌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生物学特性奠定了基础.     

Role of two-component regulatory systems in intracellular survival of Mycobacterium tuberculosis

The typical two-component regulatory systems (TCSs), consisting of response regulator and histidine kinase, play a central role in survival of pathogenic bacteria under stress conditions such as nutrient starvation, hypoxia, and nitrosative stress. A total of 11 complete paired two-component regulatory systems have been found in Mycobacterium tuberculosis, including a few isolated kinase and regulatory genes. Increasing evidence has shown that TCSs are closely associated with multiple physiological process like intracellular persistence, pathogenicity, and metabolism. This review gives the two-component signal transduction systems in M. tuberculosis and their signal transduction roles in adaption to the environment.     

极小种群植物川柿叶绿体基因组特征与系统发育分析（附录）

川柿(Diospyros sutchuensis)为极小种群和国家重点保护野生植物,分布范围狭窄,种群数量极少。目前,川柿基因组信息缺乏,在柿属(Diospyros)中的系统亲缘关系不明确。该研究通过Illumina平台对川柿叶绿体基因组进行测序,应用Getorganellev1.7.3.4和PGA软件对基因组进行组装和注释,使用DnaSP6.12.03软件进行多序列对比分析,并使用REPuter、Tandem Reapeats Finder和MISA软件进行重复序列分析,使用CodonW1.4和EasyCodemL软件分别进行密码子偏好性和选择压力分析。同时,基于4个不同的叶绿体基因组序列数据集,使用IQtree软件分析川柿与11个柿属物种的系统发育关系。结果表明:(1)川柿叶绿体基因组全长157 917 bp,包含1对26 111 bp的反向重复区、大单拷贝区(87 303 bp)和小单拷贝区(18 392 bp),GC碱基含量为37.4%。(2)川柿叶绿体基因组共注释到113个基因,包括79个蛋白编码基因、30个tRNA基因和4个rRNA基因; 共检测到49个长重复序列、27个串联重复序列和34个简单重复序列; 蛋白编码基因中高频密码子31个,多数密码子末位碱基为A或U,编码亮氨酸的密码子使用最多; 基因组编码区比非编码区更为保守,10个高变热点区域可作为潜在的分子标记; 蛋白编码基因中有8个基因(ndhB、ndhG、ndhI、rbcL、rpoB、petB、petD、rps12)受到正选择压力。(3)系统发育分析显示,川柿与老鸦柿(D. rhombifolia)和乌柿(D. cathayensis)亲缘关系最为密切,它们与海南柿(D. hainanensis)共同形成一个单系分支。该研究结果既为川柿及柿属种质资源鉴定、遗传多样性保护以及种群恢复等提供了叶绿体基因组资源,也为阐明川柿的系统进化提供了重要的分子信息。     

The complete plastid genome of Cotinus coggygria and phylogenetic analysis of the Anacardiaceae

Cotinus coggygria Scop. (Anacardiaceae) is an important ornamental tree with beautiful characteristics that is grown in China. In this study, the complete plastid genome of C. coggygria was sequenced and assembled. This genome was 158,843 bp in size and presented a typical tetrad structure, consisting of a large single-copy region (87,121 bp), a pair of inverted repeat regions (26,829 bp), and a small single-copy region (18,064 bp). A total of 134 genes were annotated, including 88 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. We observed a deletion that caused the loss of the rpl32 gene, and a small expansion of IR regions resulted in the trnH gene accessing IR regions; two copies were obtained. Phylogenetic analysis showed that C. coggygria was most closely related to Pistacia, with 100% bootstrap support within Anacardiaceae. In this study, we report the plastid genome of Cotinus species for the first time, which provides insight into the evolution of the plastid genome in Anacardiaceae and promotes the understanding of Cotinus plants.     

Comparative analysis of the complete mitochondrial genomes of three rockfishes (Scorpaeniformes,Sebastiscus) and insights into the phylogenetic relationships of Sebastidae

Mitochondrial genome is a powerful molecule marker to provide information for phylogenetic relationships and revealing molecular evolution in ichthyological studies. Sebastiscus species, a marine rockfish, are of essential economic value. However, the taxonomic status and phylogenetic relationships of Sebastidae have been controversial so far. Here, the mitochondrial genomes (mitogenomes) of three species, S. tertius, S. albofasciatus, and S. marmoratus, were systemically investigated. The lengths of the mitogenomes’ sequences of S. tertius, S. albofasciatus, and S. marmoratus were 16910, 17056, and 17580 bp, respectively. It contained 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNA (tRNA) genes, and one identical control region (D-loop) among the three species. The genetic distance and K_a/K_s ratio analyses indicated 13 PCGs were suffering purifying selection and the selection pressures were different from certain deep-sea fishes, which were most likely due to the difference in their living environment. The phylogenetic tree was constructed by Bayesian Inference (BI) and Maximum Likelihood (ML). Most interestingly, the results indicated that Sebastidae and Scorpaenidae were grouped into a separate branch, so the taxonomic status of Sebastidae should be classified into subfamily Sebastinae. Our results may lead to a taxonomic revision of Scorpaenoidei.     

Searching for potential drug targets in two-component and phosphorelay signal-transduction systems using three-dimensional cluster analysis

Two-component and phosphorelay signal transduction systems are central components in the virulence and antimicrobial resistance responses of a number of bacterial and fungal pathogens; in some cases, these systems are essential for bacterial growth and viability. Herein, we analyze in detail the conserved surface residue clusters in the phosphotransferase domain of histidine kinases and the regulatory domain of response regulators by using complex structure-based three-dimensional cluster analysis. We also investigate the protein-protein interactions that these residue clusters participate in. The Spo0B-Spo0F complex structure was used as the reference structure, and the multiple aligned sequences of phosphotransferases and response regulators were paired correspondingly. The results show that a contiguous conserved residue cluster is formed around the active site, which crosses the interface of histidine kinases and response regulators. The conserved residue clusters of phosphotransferase and the regulatory domains are directly involved in the functional implementation of two-component signal transduction systems and are good targets for the development of novel antimicrobial agents.     

Extracting sequence motifs and the phylogenetic features of SNARE-dependent membrane traffic

The SNARE proteins are required for membrane fusion during intracellular vesicular transport and for its specificity. Only the unique combination of SNARE proteins (cognates) can be bound and can lead to membrane fusion, although the characteristics of the possible specificity of the binding combinations encoded in the SNARE sequences have not yet been determined. We discovered by whole genome sequence analysis that sequence motifs (conserved sequences) in the SNARE motif domains for each protein group correspond to localization sites or transport pathways. We claim that these motifs reflect the specificity of the binding combinations of SNARE motif domains. Using these motifs, we could classify SNARE proteins from 48 organisms into their localization sites or transport pathways. The classification result shows that more than 10 SNARE subgroups are kingdom specific and that the SNARE paralogs involved in the plasma membrane-related transport pathways have developed greater variations in higher animals and higher plants than those involved in the endoplasmic reticulum-related transport pathways throughout eukaryotic evolution.     

Arabidopsis genome analysis as exemplified by analysis of chromosome 4

During the last decade the small cruciferous plant Arabidopsis thaliana has become a model organism for flowering plants. Sequencing and analysis of the Arabidopsis genome is nearing completion. Beside an overview on methods and strategies for Arabidopsis genome analysis, a summary of the results from the first analysis is presented.This includes an overview on chromosomal organisation and topological features as well as a first comparison with other genomes.     

鸽新城疫病毒BJ-C分离株基因组测序及系统发育分析

【背景】鸽新城疫是由鸽Ⅰ型副黏病毒(pigeon paramyxovirus type Ⅰ,PPMV-1)感染引起的危害最严重的疫病之一,至今尚无有效的防控制剂。【目的】分析鸽新城疫病毒BJ-C株的基因组信息及系统发育关系,为鸽新城疫的防控提供科学依据。【方法】设计首尾重叠的6对特异性引物,利用分段扩增的方法,以鸽新城疫病毒BJ-C株基因组cDNA为模板,分别扩增、测序后进行全基因组序列拼接。以NCBI数据库中发布的新城疫病毒序列为参考,针对鸽新城疫病毒BJ-C株的基因组、F基因建立系统发育树。【结果】鸽新城疫病毒BJ-C株的基因组全长为15192nt。基于全基因组的系统发育分析发现其与PPMV-1/BJ-01/CH株的系统发育关系最近,核苷酸相似性为99.96%,氨基酸相似性高达100%,属于同一个分支,而与LaSota疫苗株等其他新城疫毒株的亲缘关系相对较远。基于F基因序列的系统发育树分析发现BJ-C株F基因与我国的BJP2013株同属一个分支。ClassⅡ类Ⅵ亚型F基因高变区序列(47-420nt)比对结果显示,安徽株Pigeon/Anhui/2369/2012、广东株Pigeon/Guangdong/GZ288/2013、北京株BJP13、浙江株Pigeon/Zhejiang/2036/2012及比利时株PPMV-1/Belgium/11-09620/2011等与BJ-C毒株处于同一个分支,同属Ⅵb亚型。【结论】本研究获得了鸽新城疫病毒BJ-C株全基因组序列,分析了其系统发育关系,确定其属于ClassⅡ类Ⅵb型,为后续防控产品的开发提供了理论依据。     

Genome phylogenetic analysis based on extended gene contents

With the rapid growth of entire genome data, whole-genome approaches such as gene content become popular for genome phylogeny inference, including the tree of life. However, the underlying model for genome evolution is unclear, and the proposed (ad hoc) genome distance measure may violate the additivity. In this article, we formulate a stochastic framework for genome evolution, which provides a basis for defining an additive genome distance. However, we show that it is difficult to utilize the typical gene content data-i.e., the presence or absence of gene families across genomes-to estimate the genome distance. We solve this problem by introducing the concept of extended gene content; that is, the status of a gene family in a given genome could be absence, presence as single copy, or presence as duplicates, any of which can be used to estimate the genome distance and phylogenetic inference. Computer simulation shows that the new tree-making method is efficient, consistent, and fairly robust. The example of 35 microbial complete genomes demonstrates that it is useful not only to study the universal tree of life but also to explore the evolutionary pattern of genomes.