A weighted measure for the similarity analysis of DNA sequences

Here we propose a weighted measure for the similarity analysis of DNA sequences. It is based on LZ complexity and (0,1) characteristic sequences of DNA sequences. This weighted measure enables biologists to extract similarity information from biological sequences according to their requirements. For example, by this weighted measure, one can obtain either the full similarity information or a similarity analysis from a given biological aspect. Moreover, the length of DNA sequence is not problematic. The application of the weighted measure to the similarity analysis of β-globin genes from nine species shows its flexibility.     

一种用于蛋白质相似性分析的新的相对距离

本文论述了一种新的相对距离,用于分析不同蛋白质序列的相似性分析和构造进化树．此种距离基于Lempel-Zip复杂度,不需要进行序列比对和复杂性算法．为了说明这种距离的合理性,本文对8个物种进行了相似性分析并构造了其进化树．     

Phylogenetic analysis of DNA sequences based on the generalized pseudo-amino acid composition

The main work of this paper is to propose a new theory and method, which is based on the idea of the pseudo-amino acid composition, for phylogenetic analysis of DNA primary sequences. In our method, we revise the part of the occurrence frequency of 20 amino acids in the method of the pseudo-amino acid composition by replacing the frequency of 16 dinucleotides. And we select eight LZ complexity factors of eight (0,1) sequences of a DNA primary sequence as PseAA components. Finally, we characterize a DNA sequence with a 24-dimensional vector. We reconstruct the phylogenetic trees of two datasets. The results show that our method is efficient and significant.     

Housekeeping genes for phylogenetic analysis of eutherian relationships

The molecular relationship of placental mammals has attracted great interest in recent years. However, 2 crucial and conflicting hypotheses remain, one with respect to the position of the root of the eutherian tree and the other the relationship between the orders Rodentia, Lagomorpha (rabbits, hares), and Primates. Although most mitochondrial (mt) analyses have suggested that rodents have a basal position in the eutherian tree, some nuclear data in combination with mt-rRNA genes have placed the root on the so-called African clade or on a branch that includes this clade and the Xenarthra (e.g., anteater and armadillo). In order to generate a new and independent set of molecular data for phylogenetic analysis, we have established cDNA sequences from different tissues of various mammalian species. With this in mind, we have identified and sequenced 8 housekeeping genes with moderately fast rate of evolution from 22 placental mammals, representing 11 orders. In order to determine the root of the eutherian tree, the same genes were also sequenced for 3 marsupial species, which were used as outgroup. Inconsistent with the analyses of nuclear + mt-rRNA gene data, the current data set did not favor a basal position of the African clade or Xenarthra in the eutherian tree. Similarly, by joining rodents and lagomorphs on the same basal branch (Glires hypothesis), the data set is also inconsistent with the tree commonly favored in mtDNA analyses. The analyses of the currently established sequences have helped examination of problematic parts in the eutherian tree at the same time as they caution against suggestions that have claimed that basal eutherian relationships have been conclusively settled.     

A Novel Construction of Genome Space with Biological Geometry

A genome space is a moduli space of genomes. In this space, each point corresponds to a genome. The natural distance between two points in the genome space reflects the biological distance between these two genomes. Currently, there is no method to represent genomes by a point in a space without losing biological information. Here, we propose a new graphical representation for DNA sequences. The breakthrough of the subject is that we can construct the moment vectors from DNA sequences using this new graphical method and prove that the correspondence between moment vectors and DNA sequences is one-to-one. Using these moment vectors, we have constructed a novel genome space as a subspace in R^N. It allows us to show that the SARS-CoV is most closely related to a coronavirus from the palm civet not from a bird as initially suspected, and the newly discovered human coronavirus HCoV-HKU1 is more closely related to SARS than to any other known member of group 2 coronavirus. Furthermore, we reconstructed the phylogenetic tree for 34 lentiviruses (including human immunodeficiency virus) based on their whole genome sequences. Our genome space will provide a new powerful tool for analyzing the classification of genomes and their phylogenetic relationships.     

ITS1 Sequences of Nuclear Ribosomal DNA in Wild Rices and Cultivated Rices of China and Their Phylogenetic Implications

he first internal transcribed spacer (ITS1) of nuclear ribosomal DNA of three wild rice species and two subspecies of cultivated rice, which are distributed in China, was amplified using PCR technique and sequenced with automated fluorescent sequencing. The sequences of ITS1 ranged from 193 bp to 218 bp in size and G/C content varied from 69.3%to 72.7%. In pairwise comparison among the five taxa, sequence site divergence ranged from 1.5 % to 10.6%. Phylogenetic analysis of ITS1 sequences using Wagner parsimony generated a single well-resolved tree, which revealed that Oryza rufipogon was much more closely related to cultivated rice species than to the other two wild species. Oryza granulata was less closely related to either cultivated rice species or the other two wild species, and might be a unique and isolated taxon in the genus Oryza. The phylogenetic relationships of the three wild rice species and two cultivated rice subspecies inferred from ITS1 sequences is highly concordant with those based on the molecular evidence from isozyme, chloroplast DNA (cpDNA), mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) of the genus Oryza.     

中国野生稻及栽培稻核糖体DNA第一转录间隔区序列分析及其系统学意义

用 PCR技术从产于我国的 3种野生稻和亚洲栽培稻的 2个亚种中特异地扩增和测序了 r DNA的第一转录间隔区。普通野生稻 (Oryza rufipogon)、药用野生稻 (O.officinalis)、疣粒野生稻 (O.granu-lata)和栽培稻的两个亚种 (O.sativa ssp.indica,O.sativa ssp.japonica)的 ITS1序列为 1 93bp、1 94bp、2 1 8bp、1 94bp和 1 94bp,它们的 G/ C含量为 69.3%～ 72 .7% ,序列中位点趋异率为 1 .5%～ 1 0 .6%。序列的相似性比较和简约性分支分析的结果表明 ,普通野生稻与栽培稻的两个亚种之间的亲缘关系最为密切 ;药用野生稻与普通野生稻和与栽培稻的两个亚种的相似性都为 82 % ,说明它与 AA基因组有一定的亲缘关系 ;疣粒野生稻与普通野生稻、药用野生稻和栽培稻两个亚种的亲缘关系相对较远 ,它在稻属中可能是一个系统地位较独特的类群。以 ITS1序列构建的 3种野生稻和 2个栽培稻亚种的系统发育关系与前人用同工酶、叶绿体 DNA、线粒体 DNA和核 DNA资料重建的稻属的系统发育关系基本一致     

Mapping the Space of Genomic Signatures

We propose a computational method to measure and visualize interrelationships among any number of DNA sequences allowing, for example, the examination of hundreds or thousands of complete mitochondrial genomes. An "image distance" is computed for each pair of graphical representations of DNA sequences, and the distances are visualized as a Molecular Distance Map: Each point on the map represents a DNA sequence, and the spatial proximity between any two points reflects the degree of structural similarity between the corresponding sequences. The graphical representation of DNA sequences utilized, Chaos Game Representation (CGR), is genome- and species-specific and can thus act as a genomic signature. Consequently, Molecular Distance Maps could inform species identification, taxonomic classifications and, to a certain extent, evolutionary history. The image distance employed, Structural Dissimilarity Index (DSSIM), implicitly compares the occurrences of oligomers of length up to k (herein k = 9) in DNA sequences. We computed DSSIM distances for more than 5 million pairs of complete mitochondrial genomes, and used Multi-Dimensional Scaling (MDS) to obtain Molecular Distance Maps that visually display the sequence relatedness in various subsets, at different taxonomic levels. This general-purpose method does not require DNA sequence alignment and can thus be used to compare similar or vastly different DNA sequences, genomic or computer-generated, of the same or different lengths. We illustrate potential uses of this approach by applying it to several taxonomic subsets: phylum Vertebrata, (super)kingdom Protista, classes Amphibia-Insecta-Mammalia, class Amphibia, and order Primates. This analysis of an extensive dataset confirms that the oligomer composition of full mtDNA sequences can be a source of taxonomic information. This method also correctly finds the mtDNA sequences most closely related to that of the anatomically modern human (the Neanderthal, the Denisovan, and the chimp), and that the sequence most different from it in this dataset belongs to a cucumber.     

Template-directed arrest of mammalian mitochondrial DNA synthesis.

Mammalian mitochondrial DNA often contains a short DNA displacement loop at the heavy-strand origin of replication. This short nascent DNA molecule has been used to study site-specific termination of mitochondrial DNA synthesis in human and mouse cells. We examined D-loop strand termination in two distantly related artiodactyls, the pig and the cow. Porcine mitochondrial DNA was unique among mammals in that it contained only a single species of D-loop single-stranded DNA. Its 3' end mapped to a site 187 nucleotides from the 5' end of the proline tRNA gene. This site was 21 and 47 nucleotides 5' to two very similar sequences (5' ACATATPyATTAT 3') which are closely related to the human and mouse termination-associated sequences noted by Doda et al. (J. N. Doda, D. T. Wright, and D. A. Clayton, Proc. Nat. Acad. Sci. USA 78:616-6120, 1981). Bovine mitochondrial DNA contained three major D-loop DNA species whose 3' ends mapped to three different sites. These sites were not found in the porcine sequence. However, the bovine termination sites were located 60 to 64 base pairs 5' from sequences which were also very similar to the termination-associated sequences present in pigs and other mammals. These results firmly establish the concept that arrest of heavy-strand DNA synthesis is an event determined, at least in part, by template sequence. They also suggest that arrest is determined by sequences which are a considerable physical distance away from the actual termination site.     

Phylogenetic relationships between toxic and non-toxic strains of the genus Microcystis based on 16S to 23S internal transcribed spacer sequence

16S to 23S ribosomal DNA internal transcribed spacer sequences of 47 strains of the genus Microcystis were determined. Derived maximum likelihood and DNA distance trees indicated that Microcystis can be divided into three clusters. The first cluster included toxic and non-toxic strains, the second only toxic ones, and the third only non-toxic ones. The tree topologies were not necessarily correlated with morphospecies distinction or phycobilin pigment composition, and one genotype may have more than one morphotype. Phylogenetic analysis based on intergenic spacer sequences was thought to be effective for understanding relationships among closely related species and strains.     

Derived vocalizations of geladas (Theropithecus gelada) and the evolution of vocal complexity in primates

Primates are intensely social and exhibit extreme variation in social structure, making them particularly well suited for uncovering evolutionary connections between sociality and vocal complexity. Although comparative studies find a correlation between social and vocal complexity, the function of large vocal repertoires in more complex societies remains unclear. We compared the vocal complexity found in primates to both mammals in general and human language in particular and found that non-human primates are not unusual in the complexity of their vocal repertoires. To better understand the function of vocal complexity within primates, we compared two closely related primates (chacma baboons and geladas) that differ in their ecology and social structures. A key difference is that gelada males form long-term bonds with the 2-12 females in their harem-like reproductive unit, while chacma males primarily form temporary consortships with females. We identified homologous and non-homologous calls and related the use of the derived non-homologous calls to specific social situations. We found that the socially complex (but ecologically simple) geladas have larger vocal repertoires. Derived vocalizations of geladas were primarily used by leader males in affiliative interactions with 'their' females. The derived calls were frequently used following fights within the unit suggesting that maintaining cross-sex bonds within a reproductive unit contributed to this instance of evolved vocal complexity. Thus, our comparison highlights the utility of using closely related species to better understand the function of vocal complexity.     

拓扑树间的通经拓扑距离

给出了一种新的系统树间的拓扑距离,使用NJ,MP,UPGMA等3种方法对13种动物的线粒体中14个基因（含组合的）DNA序列数据进行系统树的构建,利用分割拓扑距离和本文给出的通经拓扑距离对这14种系统树这间及其与真树进行比较。结果显示,NJ法对获得已知树的有效率最高,MP法次之,UPGMA法最低。这14种DNA序列所构建的系统树与已知树的拓扑距离基本上是随其DNA序列长度增加而减小,但两者的相关系数并未达到显著水平,分割拓扑距离在总体上可反映树间的拓扑结构差异,但其测度精确度比通经拓扑距离要低。     

A new phylogenetic marker,apolipoprotein B,provides compelling evidence for eutherian relationships

Higher-level relationships within, and the root of Placentalia, remain contentious issues. Resolution of the placental tree is important to the choice of mammalian genome projects and model organisms, as well as for understanding the biogeography of the eutherian radiation. We present phylogenetic analyses of 63 species representing all extant eutherian mammal orders for a new molecular phylogenetic marker, a 1.3kb portion of exon 26 of the apolipoprotein B (APOB) gene. In addition, we analyzed a multigene concatenation that included APOB sequences and a previously published data set (Murphy et al., 2001b) of three mitochondrial and 19 nuclear genes, resulting in an alignment of over 17kb for 42 placentals and two marsupials. Due to computational difficulties, previous maximum likelihood analyses of large, multigene concatenations for placental mammals have used quartet puzzling, less complex models of sequence evolution, or phylogenetic constraints to approximate a full maximum likelihood bootstrap. Here, we utilize a Unix load sharing facility to perform maximum likelihood bootstrap analyses for both the APOB and concatenated data sets with a GTR+Gamma+I model of sequence evolution, tree-bisection and reconnection branch-swapping, and no phylogenetic constraints. Maximum likelihood and Bayesian analyses of both data sets provide support for the superordinal clades Boreoeutheria, Euarchontoglires, Laurasiatheria, Xenarthra, Afrotheria, and Ostentoria (pangolins+carnivores), as well as for the monophyly of the orders Eulipotyphla, Primates, and Rodentia, all of which have recently been questioned. Both data sets recovered an association of Hippopotamidae and Cetacea within Cetartiodactyla, as well as hedgehog and shrew within Eulipotyphla. APOB showed strong support for an association of tarsier and Anthropoidea within Primates. Parsimony, maximum likelihood and Bayesian analyses with both data sets placed Afrotheria at the base of the placental radiation. Statistical tests that employed APOB to examine a priori hypotheses for the root of the placental tree rejected rooting on myomorphs and hedgehog, but did not discriminate between rooting at the base of Afrotheria, at the base of Xenarthra, or between Atlantogenata (Xenarthra+Afrotheria) and Boreoeutheria. An orthologous deletion of 363bp in the aligned APOB sequences proved phylogenetically informative for the grouping of the order Carnivora with the order Pholidota into the superordinal clade Ostentoria. A smaller deletion of 237-246bp was diagnostic of the superordinal clade Afrotheria.     

Evaluating the Phylogenetic Position of Chinese Tree Shrew （Tupaia belangeri chinensis） Based on Complete Mitochondrial Genome： Implication for Using Tree Shrew as an Alternative Experimental Animal to Primates in Biomedical Research

Tree shrew(Tupaia belangeri) is currently placed in Order Scandentia and has a wide distribution in Southeast Asia and Southwest China.Due to its unique characteristics,such as small body size,high brain-to-body mass ratio,short reproductive cycle and life span,and low-cost of maintenance,tree shrew has been proposed to be an alternative experimental animal to primates in biomedical research.However,there are some debates regarding the exact phylogenetic affinity of tree shrew to primates.In this study,we determined the mtDNA entire genomes of three Chinese tree shrews(T.belangeri chinensis) and one Malayan flying lemur(Galeopterus variegatus).Combined with the published data for species in Euarchonta,we intended to discern the phylogenetic relationship among representative species of Dermoptera,Scandentia and Primates.The mtDNA genomes of Chinese tree shrews and Malayan flying lemur shared similar gene organization and structure with those of other mammals.Phylogenetic analysis based on 12 concatenated mitochondrial proteinencoding genes revealed a closer relationship between species of Scandentia and Glires,whereas species of Dermoptera were clustered with Primates.This pattern was consistent with previously reported phylogeny based on mtDNA data,but differed from the one reconstructed on the basis of nuclear genes.Our result suggested that the matrilineal affinity of tree shrew to primates may not be as close as we had thought.The ongoing project for sequencing the entire genome of Chinese tree shrew will provide more information to clarify this important issue.     

Phylogenetic relationships of Empetraceae inferred from sequences of chloroplast gene matK and nuclear ribosomal DNA ITS region

We used sequences of nrDNA ITS and chloroplast gene matK to evaluate the monophyly of Empetrum and Corema and to examine phylogenetic relationships of the Empetraceae. Sequences of these two DNA markers were obtained for 11 plant samples, representing species of Empetrum from both the Southern and Northern Hemispheres, species and subspecies of Corema, and the monotypic Ceratiola. Sequences of four species of Rhododendron were used for rooting purposes. Our results show that species of Empetrum form a clade sister to the clade containing both Corema and Ceratiola. These two clades are strongly supported in both the matK and ITS trees, suggesting that Ceratiola is more closely related to Corema than to Empetrum, and is not of a hybrid origin between the ancestors of the latter two genera. In the matK tree, Corema conradii is more closely related to Ceratiola than to Corema album and C. album subsp. azoricum, whereas in the ITS tree, Ceratiola is allied with Corema album and C. album subsp. azoricum. This suggests that C. conradii might be a hybrid between ancestral populations of Ceratiola and C. album. The monophyly of Empetrum rejects the hypothesis of its independent origin in the two Hemispheres. Our trees also suggest the fact that the modern amphitropical distribution of Empetrum is the result of long distance dispersal, not of the vicarious events.     

The structural complexity of DNA templates—Implications on cellular complexity

A previously formulated procedure for the quantitative evaluation of the complexities of molecules and biostructures is applied to assess the complexities of selected genomic DNA sequences. These include: (1) Several E. coli genes, including lacI, as examples of DNA sequences which are nearly as complex as possible (relative complexity=∼1). This is verified by the Lempel-Ziv (LZ) complexity analysis. (2) The telomere of a yeast chromosome, which has a considerable number of regular features that reduce complexity; the telomere shows indeed a lower structural complexity value. (3) A segment of human DNA, gene p53, which has a certain number of regular features such as 29 interspersed alu elements; these features cause a certain reduction in the complexity of the p53 gene, but do not invalidate the (previous) overall conclusion that template complexity is very high. The close to maximal complexity of the transcribed regions of p53 is validated by the LZ compression analysis. The general conclusion is that DNA base sequence composition is the dominant factor determining cellular complexity. The high complexity of DNA arrived at is a direct consequence of the template character of DNA and reflects the role of genomic DNA as a principal regulating element of a cell. It will be a challenge to find systems of lower complexity with the ability to respond to challenges from the environment to the extent that DNA templated systems do. Cellular complexity and template directed activity are thus highly intertwined properties, at the heart of many developmental, behavioral and evolutionary processes.     

Genome analysis with the conditional multinomial distribution profile

The focus of the research is on the analysis of genome sequences. Based on the inter-nucleotide distance sequence, we propose the conditional multinomial distribution profile for the complete genomic sequence. These profiles can be used to define a very simple, computationally efficient, alignment-free, distance measure that reflects the evolutionary relationships between genomic sequences. We use this distance measure to classify chromosomes according to species of origin, to build the phylogenetic tree of 24 complete genome sequences of coronaviruses. Our results demonstrate the new method is powerful and efficient.     

Asynchronous distance between homologous DNA sequences

The distance between homologous DNA sequences of two species is proposed to be -1/4 ln[det(P)], where P is the conditional probability matrix specifying the proportions of the various nucleotides in the second sequence, corresponding to each of the four nucleotides in the first sequence. A probability model is described which supports this choice of distance. Distance measures based on a constant evolutionary rate assumption are described and compared with the proposed measure. Sampling properties of both types of distance are examined and we conclude by applying the distance measures to mitochondrial DNA sequences of the hominoids.