Frank Wright

TOPALi v2

Summary: TOPALi v2 simplifies and automates the use of several methods for the evolutionary analysis of multiple sequence alignments. Jobs are submitted from a Java graphical user interface as TOPALi web services to either run remotely on high-performance computing clusters or locally (with multiple cores supported). Methods available include model selection and phylogenetic tree estimation using the Bayesian inference and maximum likelihood (ML) approaches, in addition to recombination detection methods. The optimal substitution model can be selected for protein or nucleic acid (standard, or protein-coding using a codon position model) data using accurate statistical criteria derived from ML co-estimation of the tree and the substitution model. Phylogenetic software available includes PhyML, RAxML and MrBayes. Availability: Freely downloadable from http://www.topali.org for Windows, Mac OS X, Linux and Solaris. Contact: iain.milne@scri.ac.uk

Codon usage in the G+C-rich Streptomyces genome.

The codon usage (CU) patterns of 64 genes from the Gram+ prokaryotic genus Streptomyces were analysed. Despite the extremely high overall G+C content of the Streptomyces genome (estimated at 0.74), individual genes varied in G+C content from 0.610 to 0.797, and had third codon position G+C contents (GC3s) that varied from 0.764 to 0.983. The variation in GC3s explains a significant proportion of the variation in CU patterns. This is consistent with an evolutionary model of the Streptomyces genome where biased mutation pressure has led to a high average G+C content with random variation about the mean, although the variation observed is greater than that expected from a simple binomial model. The only gene in the sample that can be confidently predicted to be highly expressed, EF-Tu of Streptomyces coelicolor A3(2) (GC3s = 0.927), shows a preference for a third position C in several of the four codon families, and for CGY and GGY for Arg and Gly codons, respectively (Y = pyrimidine); similar CU patterns are found in highly expressed genes of the G+C-rich Micrococcus luteus genome. It thus appears that codon usage in Streptomyces is determined predominantly by mutation bias, with weak translational selection operating only in highly expressed genes. We discuss the possible consequences of the extreme codon bias of Streptomyces and consider how it may have evolved. A set of CU tables is provided for use with computer programs that locate protein-coding regions.

TOPALi: software for automatic identification of recombinant sequences within DNA multiple alignments

SUMMARY TOPALi is a new Java graphical analysis application that allows the user to identify recombinant sequences within a DNA multiple alignment (either automatically or via manual investigation). TOPALi allows a choice of three statistical methods to predict the positions of breakpoints due to past recombination. The breakpoint predictions are then used to identify putative recombinant sequences and their relationships to other sequences. In addition to its sophisticated interface, TOPALi can import many sequence formats, estimate and display phylogenetic trees and allow interactive analysis and/or automatic HTML report generation. AVAILABILITY TOPALi is freely available from http://www.bioss.ac.uk/software.html

Identification and localisation of the NB-LRR gene family within the potato genome.

BackgroundThe potato genome sequence derived from the Solanum tuberosum Group Phureja clone DM1-3 516 R44 provides unparalleled insight into the genome composition and organisation of this important crop. A key class of genes that comprises the vast majority of plant resistance (R) genes contains a nucleotide-binding and leucine-rich repeat domain, and is collectively known as NB-LRRs.ResultsAs part of an effort to accelerate the process of functional R gene isolation, we performed an amino acid motif based search of the annotated potato genome and identified 438 NB-LRR type genes among the ~39,000 potato gene models. Of the predicted genes, 77 contain an N-terminal toll/interleukin 1 receptor (TIR)-like domain, and 107 of the remaining 361 non-TIR genes contain an N-terminal coiled-coil (CC) domain. Physical map positions were established for 370 predicted NB-LRR genes across all 12 potato chromosomes. The majority of NB-LRRs are physically organised within 63 identified clusters, of which 50 are homogeneous in that they contain NB-LRRs derived from a recent common ancestor.ConclusionsBy establishing the phylogenetic and positional relationship of potato NB-LRRs, our analysis offers significant insight into the evolution of potato R genes. Furthermore, the data provide a blueprint for future efforts to identify and more rapidly clone functional NB-LRR genes from Solanum species.

High diversity of arbuscular mycorrhizal fungi in a boreal herb-rich coniferous forest

* Here, the diversity of arbuscular mycorrhizal (AM) fungi was determined in a boreal herb-rich coniferous forest in relation to environmental variables. * Root samples of five plant species (Fragaria vesca, Galeobdolon luteum, Hepatica nobilis, Oxalis acetosella and Trifolium pratense) were analysed from stands differing in age and forest management intensity. * Thirty-four Glomeromycota taxa (small-subunit ribosomal RNA gene (SSU rDNA) sequence groups) were detected from 90 root samples (911 clones), including eight new taxa. Sequence groups related to Glomus intraradices were most common (MO-G3 and MO-G13). Samples of H. nobilis were colonized by more AM fungal taxa (3.68 +/- 0.31) than those of O. acetosella (2.69 +/- 0.34), but did not differ significantly in this respect from those of F. vesca (3.15 +/- 0.38). Effects of forest management, host plant species (except above) or season on the number or composition of fungal taxa in root samples were not detected, and neither were they explained by environmental variables (vegetation, soil and light conditions). * This is the most taxon-rich habitat described to date in terms of root-colonizing Glomeromycota. The data demonstrate the importance of temperate coniferous forests as habitats for AM fungi and plants. Lack of obvious fungal community patterns suggests more complex effects of biotic and abiotic factors, and possibly no adverse effect of common forest management practices on AM fungal diversity.

Molecular Characterization of Pigmented and Nonpigmented Isolates of Mycobacterium avium subsp. paratuberculosis

ABSTRACT Five pigmented isolates of Mycobacterium avium subsp. paratuberculosis were examined by pulsed-field gel electrophoresis (PFGE), IS900 restriction fragment length polymorphism (IS900-RFLP), and IS1311 polymorphism analysis using PCR. All of the pigmented isolates exhibited one of three distinct PFGE profiles with SnaBI, designated 9, 10, and 11, and with SpeI, designated 7, 8, and 9, which generated three multiplex profiles designated [9-7], [10-8], and [11-9]. All of the pigmented isolates had the same IS900-RFLP BstEII and PvuII profiles. The IS900-RFLP BstEII profile was new, but the IS900-RFLP PvuII profile corresponded to PvuII type 6 of a sheep strain described by Cousins and colleagues (D. V. Cousins, S. N. Williams, A. Hope, and G. J. Eamens, Aust. Vet. J. 78:184-190, 2000). IS1311-PCR analysis typed all of the pigmented isolates as sheep (S) strains. The genetic relationship between pigmented and nonpigmented isolates was investigated by using multiplex PFGE data from the analysis of both the 5 pigmented isolates and 88 nonpigmented isolates of M. avium subsp. paratuberculosis from a variety of host species and geographic locations. It was possible to classify the isolates into two distinct types designated type I, comprising the pigmented isolates, and type II, comprising the nonpigmented isolates, which exhibit a very broad host range.

TOPAL 2.0: improved detection of mosaic sequences within multiple alignments

MOTIVATION The Dss statistic was proposed by McGuire et al. (Mol. Biol. Evol., 14, 1125-1131, 1997) for scanning data sets for the presence of recombination, an important step in some phylogenetic analyses. The statistic, however, could not distinguish well between among-site rate variation and recombination, and had no statistical test for significant values. This paper addresses these shortfalls. RESULTS A modification to the Dss statistic is proposed which accounts for rate variation to a large extent. A statistical test, based on parametric bootstrapping, is also suggested. AVAILABILITY The TOPAL package (version 2) may be accessed from http:/ /www.bioss.sari.ac.uk/frank/Genetics and by anonymous ftp from typ://ftp.bioss.sari.ac.uk in the directory pub/phylogeny/topal. CONTACT frank@bioss.sari.ac.uk

Genetic Diversity and Phylogenetic Relationships among Strains of Prevotella (Bacteroides) ruminicola from the Rumen

A high degree of genetic diversity among 29 strains of Prevotella (Bacteroides) ruminicola from the rumen was revealed by comparing restriction fragment length polymorphisms in 16S rRNA genes, sodium dodecyl sulfate-polyacrylamide gel profiles of total-cell proteins, and G + C contents of chromosomal DNAs. In order to obtain information on phylogenetic relationships, the sequences of a 389-bp region of the 16S rRNA gene, including variable regions 4 and 5, were compared for 10 strains. These 10 strains formed a single group when their sequences were compared with 16S ribosomal DNA sequences from other species, including Bacteroides spp. from the human colon. On the other hand, the great genetic distances between many P. ruminicola strains, including P. ruminicola subsp. brevis B(1)4 and GA33 and P. ruminicola 23T (T = type strain), support the hypothesis that these organisms should be reclassified into new species. We identified signature oligonucleotides based on 16S ribosomal DNA sequences that distinguished strains related to strains 23T, B(1)4, GA33, and M384, as well as an oligonucleotide that specifically recognized all but one of the Bacteroides and Prevotella strains tested. On the basis of the priming activities of these signature oligonucleotides in PCR reactions and on other criteria, we concluded that 12 of the original 29 strains were related to strain 23T, 4 were related to strain B(1)4, and 4 were related to strain GA33. While there are clear grounds for subdividing the species P. ruminicola on the basis of genotypic differences, it is appropriate to delay formal reclassification until further work on the phenotypic differentiation of the new groups is completed.

Expression profiling of potato germplasm differentiated in quality traits leads to the identification of candidate flavour and texture genes

Quality traits such as flavour and texture are assuming a greater importance in crop breeding programmes. This study takes advantage of potato germplasm differentiated in tuber flavour and texture traits. A recently developed 44 000-element potato microarray was used to identify tuber gene expression profiles that correspond to differences in tuber flavour and texture as well as carotenoid content and dormancy characteristics. Gene expression was compared in two Solanum tuberosum group Phureja cultivars and two S. tuberosum group Tuberosum cultivars; 309 genes were significantly and consistently up-regulated in Phureja, whereas 555 genes were down-regulated. Approximately 46% of the genes in these lists can be identified from their annotation and amongst these are candidates that may underpin the Phureja/Tuberosum trait differences. For example, a clear difference in the cooked tuber volatile profile is the higher level of the sesquiterpene α-copaene in Phureja compared with Tuberosum. A sesquiterpene synthase gene was identified as being more highly expressed in Phureja tubers and its corresponding full-length cDNA was demonstrated to encode α-copaene synthase. Other potential ‘flavour genes’, identified from their differential expression profiles, include those encoding branched-chain amino acid aminotransferase and a ribonuclease suggesting a mechanism for 5′-ribonucleotide formation in potato tubers on cooking. Major differences in the expression levels of genes involved in cell wall biosynthesis (and potentially texture) were also identified, including genes encoding pectin acetylesterase, xyloglucan endotransglycosylase and pectin methylesterase. Other gene expression differences that may impact tuber carotenoid content and tuber life-cycle phenotypes are discussed.

Detection of Recombination in DNA Multiple Alignments with Hidden Markov Models

Conventional phylogenetic tree estimation methods assume that all sites in a DNA multiple alignment have the same evolutionary history. This assumption is violated in data sets from certain bacteria and viruses due to recombination, a process that leads to the creation of mosaic sequences from different strains and, if undetected, causes systematic errors in phylogenetic tree estimation. In the current work, a hidden Markov model (HMM) is employed to detect recombination events in multiple alignments of DNA sequences. The emission probabilities in a given state are determined by the branching order (topology) and the branch lengths of the respective phylogenetic tree, while the transition probabilities depend on the global recombination probability. The present study improves on an earlier heuristic parameter optimization scheme and shows how the branch lengths and the recombination probability can be optimized in a maximum likelihood sense by applying the expectation maximization (EM) algorithm. The novel algorithm is tested on a synthetic benchmark problem and is found to clearly outperform the earlier heuristic approach. The paper concludes with an application of this scheme to a DNA sequence alignment of the argF gene from four Neisseria strains, where a likely recombination event is clearly detected.