Whole-proteome phylogeny of large dsDNA viruses and parvoviruses through a composition vector method related to dynamical language model

Background  

The vast sequence divergence among different virus groups has presented a great challenge to alignment-based analysis of virus phylogeny. Due to the problems caused by the uncertainty in alignment, existing tools for phylogenetic analysis based on multiple alignment could not be directly applied to the whole-genome comparison and phylogenomic studies of viruses. There has been a growing interest in alignment-free methods for phylogenetic analysis using complete genome data. Among the alignment-free methods, a dynamical language (DL) method proposed by our group has successfully been applied to the phylogenetic analysis of bacteria and chloroplast genomes.     

Genome-based phylogeny of dsDNA viruses by a novel alignment-free method

Sequence alignment is not directly applicable to whole genome phylogeny since several events such as rearrangements make full length alignments impossible. Here, a novel alignment-free method derived from the standpoint of information theory is proposed and used to construct the whole-genome phylogeny for a population of viruses from 13 viral families comprising 218 dsDNA viruses. The method is based on information correlation (IC) and partial information correlation (PIC). We observe that (i) the IC-PIC tree segregates the population into clades, the membership of each is remarkably consistent with biologist's systematics only with little exceptions; (ii) the IC-PIC tree reveals potential evolutionary relationships among some viral families; and (iii) the IC-PIC tree predicts the taxonomic positions of certain “unclassified” viruses. Our approach provides a new way for recovering the phylogeny of viruses, and has practical applications in developing alignment-free methods for sequence classification.     

A fungal phylogeny based on 82 complete genomes using the composition vector method

Background  

Molecular phylogenetics and phylogenomics have greatly revised and enriched the fungal systematics in the last two decades. Most of the analyses have been performed by comparing single or multiple orthologous gene regions. Sequence alignment has always been an essential element in tree construction. These alignment-based methods (to be called the standard methods hereafter) need independent verification in order to put the fungal Tree of Life (TOL) on a secure footing. The ever-increasing number of sequenced fungal genomes and the recent success of our newly proposed alignment-free composition vector tree (CVTree, see Methods) approach have made the verification feasible.     

Capsid assembly of dsDNA viruses

Capsid assembly has been investigated for many dsDNA viruses. Much is known about assembly pathways for these viruses and information is now emerging about the protein/protein and protein/DNA interactions that regulate these pathways. Common themes include progressive transitions in the properties of the capsid subunits as assembly proceeds and the use of accessory factors to aid assembly. Other assembly strategies, such as packaging DNA into a preformed capsid or assembling capsids from subunit oligomers are shared among some but not all viruses of this group. Despite the many similarities, there are sufficient fundamental differences in assembly strategies that it is not possible to describe them all as variations on a single assembly theme.     

Optimal choice of k-mer in composition vector method for genome sequence comparison

Several proteins and genes are members of families that share a public evolutionary. In order to outline the evolutionary relationships and to recognize conserved patterns, sequence comparison becomes an emerging process. The current work investigates critically the k-mer role in composition vector method for comparing genome sequences. Generally, composition vector methods using k-mer are applied under choice of different value of k to compare genome sequences. For some values of k, results are satisfactory, but for other values of k, results are unsatisfactory. Standard composition vector method is carried out in the proposed work using 3-mer string length. In addition, special type of information based similarity index is used as a distance measure. It establishes that use of 3-mer and information based similarity index provide satisfactory results especially for comparison of whole genome sequences in all cases. These selections provide a sort of unified approach towards comparison of genome sequences.     

Whole genome amplification using single-primer PCR

Comprehensive genomic molecular analyses require relatively large DNA amounts that are often not available from forensic, clinical and other crucial biological samples. Numerous methods to amplify the whole genome have been proposed for cancer, forensic and taxonomic research. Unfortunately, when using truly random primers for the initial priming step, all of these procedures suffer from high background problems for sub-nanogram quantities of input DNA. Here we report an approach to eliminate this problem for PCR-based methods even at levels of DNA approaching that of a single cell.     

A comprehensive whole genome bacterial phylogeny using correlated peptide motifs defined in a high dimensional vector space

As whole genome sequences continue to expand in number and complexity, effective methods for comparing and categorizing both genes and species represented within extremely large datasets are required. Methods introduced to date have generally utilized incomplete and likely insufficient subsets of the available data. We have developed an accurate and efficient method for producing robust gene and species phylogenies using very large whole genome protein datasets. This method relies on multidimensional protein vector definitions supplied by the singular value decomposition (SVD) of a large sparse data matrix in which each protein is uniquely represented as a vector of overlapping tetrapeptide frequencies. Quantitative pairwise estimates of species similarity were obtained by summing the protein vectors to form species vectors, then determining the cosines of the angles between species vectors. Evolutionary trees produced using this method confirmed many accepted prokaryotic relationships. However, several unconventional relationships were also noted. In addition, we demonstrate that many of the SVD-derived right basis vectors represent particular conserved protein families, while many of the corresponding left basis vectors describe conserved motifs within these families as sets of correlated peptides (copeps). This analysis represents the most detailed simultaneous comparison of prokaryotic genes and species available to date.     

Bottom-up genome assembly using the Bacillus subtilis genome vector

We established a protocol to construct complete recombinant genomes from their small contiguous DNA pieces and obtained the genomes of mouse mitochondrion and rice chloroplast using a B. subtilis genome (BGM) vector. This method allows the design of any recombinant genomes, valuable not only for fundamental research in systems biology and synthetic biology but also for various applications in the life sciences.     

Multidimensional vector space representation for convergent evolution and molecular phylogeny

With growing amounts of genome data and constant improvement of models of molecular evolution, phylogenetic reconstruction became more reliable. However, our knowledge of the real process of molecular evolution is still limited. When enough large-sized data sets are analyzed, any subtle biases in statistical models can support incorrect topologies significantly because of the high signal-to-noise ratio. We propose a procedure to locate sequences in a multidimensional vector space (MVS), in which the geometry of the space is uniquely determined in such a way that the vectors of sequence evolution are orthogonal among different branches. In this paper, the MVS approach is developed to detect and remove biases in models of molecular evolution caused by unrecognized convergent evolution among lineages or unexpected patterns of substitutions. Biases in the estimated pairwise distances are identified as deviations (outliers) of sequence spatial vectors from the expected orthogonality. Modifications to the estimated distances are made by minimizing an index to quantify the deviations. In this way, it becomes possible to reconstruct the phylogenetic tree, taking account of possible biases in the model of molecular evolution. The efficacy of the modification procedure was verified by simulating evolution on various topologies with rate heterogeneity and convergent change. The phylogeny of placental mammals in previous analyses of large data sets has varied according to the genes being analyzed. Systematic deviations caused by convergent evolution were detected by our procedure in all representative data sets and were found to strongly affect the tree structure. However, the bias correction yielded a consistent topology among data sets. The existence of strong biases was validated by examining the sites of convergent evolution between the hedgehog and other species in mitochondrial data set. This convergent evolution explains why it has been difficult to determine the phylogenetic placement of the hedgehog in previous studies.     

Whole genome sequencing of Theileria parva using target capture

Protozoan parasite isolation and purification are laborious and time-consuming processes required for high quality genomic DNA used in whole genome sequencing. The objective of this study was to capture whole Theileria parva genomes directly from cell cultures and blood samples using RNA baits. Cell culture material was bait captured or sequenced directly, while blood samples were all captured. Baits had variable success in capturing T. parva genomes from blood samples but were successful in cell cultures. Genome mapping uncovered extensive host contamination in blood samples compared to cell cultures. Captured cell cultures had over 81 fold coverage for the reference genome compared to 0–33 fold for blood samples. Results indicate that baits are specific to T. parva, are a good alternative to conventional methods and thus ideal for genomic studies. This study also reports the first whole genome sequencing of South African T. parva.     

A molecular phylogeny for the large African orchid genus Disa

Phylogenetic relationships were inferred for the African subtribe Disinae (Orchidoideae, Orchidaceae), which include the large genus Disa and the small genus Schizodium. One nuclear (ITS) gene region and two plastid (trnLF and matK) gene regions were sequenced for 136 ingroup, representing 70% of all known Disinae species, as well as for 7 outgroup taxa. The combined data matrix contained 4094 characters and was analysed using parsimony and Bayesian inference. Our results show that the generic status of Schizodium can no longer be supported, as it is deeply embedded within the genus Disa. Furthermore, the currently recognised subgenera do not reflect the phylogenetic relationships and should be rejected. Several of the currently recognised sections are monophyletic, others contain misplaced elements, while some are polyphyletic. Morphological divergence, rather than convergence, has hampered previous attempts at a phylogenetic classification of the Disinae. On the basis of our molecular phylogenetic hypothesis, we propose a monotypic subtribe Disinae and a subdivision of the genus Disa into 18 sections.     

The mitochondrial genome sequence and molecular phylogeny of the turkey, Meleagris gallopavo

The mitochondrial genome (mtGenome) has been little studied in the turkey ( Meleagris gallopavo ), a species for which there is no publicly available mtGenome sequence. Here, we used PCR-based methods with 19 pairs of primers designed from the chicken and other species to develop a complete turkey mtGenome sequence. The entire sequence (16 717 bp) of the turkey mtGenome was obtained, and it exhibited 85% similarity to the chicken mtGenome sequence. Thirteen genes and 24 RNAs (22 tRNAs and 2 rRNAs) were annotated. An mtGenome-based phylogenetic analysis indicated that the turkey is most closely related to the chicken, Gallus gallus , and quail, Corturnix japonica . Given the importance of the mtGenome, the present work adds to the growing genomic resources needed to define the genetic mechanisms that underlie some economically significant traits in the turkey.     

Whole genome amplification of buccal cytobrush DNA collected for molecular epidemiology studies

Abstract

When cytobrush buccal cell samples have been collected as a genomic DNA (gDNA) source for an epidemiological study, whole genome amplification (WGA) can be critical to maintain sufficient DNA for genotyping. We evaluated REPLI-g? WGA using gDNA from two paired cytobrushes (cytobush ‘A’ kept in a cell lysis buffer, and ‘B’ dried and kept at room temperature for 3 days, and frozen until DNA extraction) in a pilot study (n=21), and from 144 samples collected by mail in a breast cancer study. WGA success was assessed as the per cent completion/concordance of STR/SNP genotypes. Locus amplification bias was assessed using quantitative PCR of 23 human loci. The pilot study showed > 98% completion but low genotype concordance between cytobrush wgaDNA and paired blood gDNA (82% and 84% for cytobrushes A and B, respectively). Substantial amplification bias was observed with significantly lower human gDNA amplification from cytobrush B than A. Using cytobrush gDNA samples from the breast cancer study (n =20), an independent laboratory demonstrated that increasing template gDNA to the REPLI-g reaction improved genotype performance for 49 SNPs; however, average completion and concordance remained below 90%. To reduce genotype misclassification when cytobrush wgaDNA is used, inclusion of paired gDNA/wgaDNA and/or duplicate wgaDNA samples is critical to monitor data quality.     

Cirripede phylogeny using a novel approach: molecular morphometrics

We present a new method using nucleic acid secondary structure to assess phylogenetic relationships among species. In this method, which we term "molecular morphometrics," the measurable structural parameters of the molecules (geometrical features, bond energies, base composition, etc.) are used as specific characters to construct a phylogenetic tree. This method relies both on traditional morphological comparison and on molecular sequence comparison. Applied to the phylogenetic analysis of Cirripedia, molecular morphometrics supports the most recent morphological analyses arguing for the monophyly of Cirripedia sensu stricto (Thoracica + Rhizocephala + Acrothoracica). As a proof, a classical multiple alignment was also performed, either using or not using the structural information to realign the sequence segments considered in the molecular morphometrics analysis. These methods yielded the same tree topology as the direct use of structural characters as a phylogenetic signal. By taking into account the secondary structure of nucleic acids, the new method allows investigators to use the regions in which multiple alignments are barely reliable because of a large number of insertions and deletions. It thus appears to be complementary to classical primary sequence analysis in phylogenetic studies.     

Studying plant genome variation using molecular markers

The authors studies on the organization and variation of plant genome with the use of molecular markers are briefly reviewed with special emphasis on random amplified polymorphic DNA (RAPD), inter simple sequence repeat (ISSR), sequence characterized amplified region (SCAR), and cleaved amplified polymorphic sequence (CAPS) markers detected with the use of polymerase chain reaction (PCR). These markers have been demonstrated to be promising for identifying cultivars and determining the purity of genetic strains of pea. Genetic relationships between strains, cultivars, and mutants of pea have been studied. The role of molecular markers in molecular genetic mapping and localizing the genes of commercially important characters of pea has been shown. The possibility of the use of molecular markers for studying somaclonal variation and detecting mutagenic factors in plants during long-term spaceflights is considered. The prospects of using DNA markers for understanding the organization and variability of higher plant genomes are discussed.__________Translated from Genetika, Vol. 41, No. 4, 2005, pp. 480–492.Original Russian Text Copyright © 2005 by Gostimsky, Kokaeva, Konovalov.     

Studying plant genome variation using molecular markers

The authors' studies on the organization and variation of plant genome with the use of molecular markers are briefly reviewed with special emphasis on random amplified polymorphic DNA (RAPD), inter simple sequence repeat (ISSR), sequence characterized amplified region (SCAR), and cleaved amplified polymorphic sequence (CAPS) markers detected with the use of polymerase chain reaction (PCR). These markers have been demonstrated to be promising for identifying cultivars and determining the purity of genetic strains of pea. Genetic relationships between strains, cultivars, and mutants of pea have been studied. The role of molecular markers in molecular genetic mapping and localizing the genes of commercially important characters of pea has been shown. The possibility of the use of molecular markers for studying somaclonal variation and detecting mutagenic factors in plants during long-term spaceflights is considered. The prospects of using DNA markers for understanding the organization and variability of higher plant genomes are discussed.     

Whole genome amplification of buccal cytobrush DNA collected for molecular epidemiology studies.

When cytobrush buccal cell samples have been collected as a genomic DNA (gDNA) source for an epidemiological study, whole genome amplification (WGA) can be critical to maintain sufficient DNA for genotyping. We evaluated REPLI-g WGA using gDNA from two paired cytobrushes (cytobush 'A' kept in a cell lysis buffer, and 'B' dried and kept at room temperature for 3 days, and frozen until DNA extraction) in a pilot study (n=21), and from 144 samples collected by mail in a breast cancer study. WGA success was assessed as the per cent completion/concordance of STR/SNP genotypes. Locus amplification bias was assessed using quantitative PCR of 23 human loci. The pilot study showed > 98% completion but low genotype concordance between cytobrush wgaDNA and paired blood gDNA (82% and 84% for cytobrushes A and B, respectively). Substantial amplification bias was observed with significantly lower human gDNA amplification from cytobrush B than A. Using cytobrush gDNA samples from the breast cancer study (n =20), an independent laboratory demonstrated that increasing template gDNA to the REPLI-g reaction improved genotype performance for 49 SNPs; however, average completion and concordance remained below 90%. To reduce genotype misclassification when cytobrush wgaDNA is used, inclusion of paired gDNA/wgaDNA and/or duplicate wgaDNA samples is critical to monitor data quality.