A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data

We describe a general likelihood-based 'mixture model' for inferring phylogenetic trees from gene-sequence or other character-state data. The model accommodates cases in which different sites in the alignment evolve in qualitatively distinct ways, but does not require prior knowledge of these patterns or partitioning of the data. We call this qualitative variability in the pattern of evolution across sites "pattern-heterogeneity" to distinguish it from both a homogenous process of evolution and from one characterized principally by differences in rates of evolution. We present studies to show that the model correctly retrieves the signals of pattern-heterogeneity from simulated gene-sequence data, and we apply the method to protein-coding genes and to a ribosomal 12S data set. The mixture model outperforms conventional partitioning in both these data sets. We implement the mixture model such that it can simultaneously detect rate- and pattern-heterogeneity. The model simplifies to a homogeneous model or a rate- variability model as special cases, and therefore always performs at least as well as these two approaches, and often considerably improves upon them. We make the model available within a Bayesian Markov-chain Monte Carlo framework for phylogenetic inference, as an easy-to-use computer program.

Profiling dehydrin gene sequence and physiological parameters in drought tolerant and susceptible spring wheat cultivars.

Physiological and yield traits such as stomatal conductance (mmol m-2s-1), Leaf relative water content (RWC %) and grain yield per plant were studied in a separate experiment. Results revealed that five out of sixteen cultivars viz. Anmol, Moomal, Sarsabz, Bhitai and Pavan, appeared to be relatively more drought tolerant. Based on morphophysiological results, studies were continued to look at these cultivars for drought tolerance at molecular level. Initially, four well recognized primers for dehydrin genes (DHNs) responsible for drought induction in T. durum L., T. aestivum L. and O. sativa L. were used for profiling gene sequence of sixteen wheat cultivars. The primers amplified the DHN genes variably like Primer WDHN13 (T. aestivum L.) amplified the DHN gene in only seven cultivars whereas primer TdDHN15 (T. durum L.) amplified all the sixteen cultivars with even different DNA banding patterns some showing second weaker DNA bands. Third primer TdDHN16 (T. durum L.) has shown entirely different PCR amplification prototype, specially showing two strong DNA bands while fourth primer RAB16C (O. sativa L.) failed to amplify DHN gene in any of the cultivars. Examination of DNA sequences revealed several interesting features. First, it identified the two exon/one intron structure of this gene (complete sequences were not shown), a feature not previously described in the two database cDNA sequences available from T. aestivum L. (gi|21850). Secondly, the analysis identified several single nucleotide polymorphisms (SNPs), positions in gene sequence. Although complete gene sequence was not obtained for all the cultivars, yet there were a total of 38 variable positions in exonic (coding region) sequence, from a total gene length of 453 nucleotides. Matrix of SNP shows these 37 positions with individual sequence at positions given for each of the 14 cultivars (sequence of two cultivars was not obtained) included in this analysis. It demonstrated a considerable diversity for this gene with only three cultivars i.e. TJ-83, Marvi and TD-1 being similar to the consensus sequence. All other cultivars showed a unique combination of SNPs. In order to prove a functional link between these polymorphisms and drought tolerance in wheat, it would be necessary to conduct a more detailed study involving directed mutation of this gene and DHN gene expression.

Distribution, gene sequence and expression in vivo of the plasmid encoded fimbrial antigen of Salmonella serotype Enteritidis

The pefA gene which encoded the serotype associated plasmid (SAP) mediated fimbrial major subunit antigen of Salmonella enterica serotype Typhimurium shared genetic identity with 128 of 706 salmonella isolates as demonstrated by dot (colony) hybridization. Seventy-seven of 113 isolates of Typhimurium and individual isolates of serotypes Bovis-morbificans, Cholerae-suis and Enteritidis phage type 9b hybridized pefA strongly, whereas 48 isolates of Enteritidis hybridized pefA weakly and one Enteritidis isolate of phage type 14b failed to hybridize. Individual isolates of 294 serotypes and 247 individual isolates of serotype Dublin did not hybridize pefA. Southern hybridization of plasmids extracted from Enteritidis demonstrated that the pefA gene probe hybridized strongly an atypical SAP of 80 kb in size harboured by one Enteritidis isolate of phage-type 9b, whereas the typical SAP of 58 kb in size harboured by 48 Enteritidis isolates hybridized weakly. One Enteritidis isolate of phage type 14b which failed to hybridize pefA in dot (colony) hybridization experiments was demonstrated to be plasmid free. A cosmid library of Enteritidis phage type 4 expressed in Escherichia coli K12 was screened by hybridization for the presence of pef sequences. Recombinant clones which were deduced to harbour the entire pef operon elaborated a PEF-like fimbrial structure at the cell surface. The PEF-like fimbrial antigen was purified from one cosmid clone and used in western blot experiments with sera from chickens infected with Enteritidis phage-type 4. Seroconversion to the fimbrial antigen was observed which indicated that the Enteritidis PEF-like fimbrial structure was expressed at some stage during infection. Nucleotide sequence analysis demonstrated that the pefA alleles of Typhimurium and Enteritidis phage-type 4 shared 76% DNA nucleotide and 82% deduced amino acid sequence identity.

Speciation and bursts of evolution

A longstanding debate in evolutionary biology concerns whether species diverge gradually through time or by rapid punctuational bursts at the time of speciation. The theory of punctuated equilibrium states that evolutionary change is characterised by short periods of rapid evolution followed by longer periods of stasis in which no change occurs. Despite years of work seeking evidence for punctuational change in the fossil record, the theory remains contentious. Further there is little consensus as to the size of the contribution of punctuational changes to overall evolutionary divergence. Here we review recent developments which show that punctuational evolution is common and widespread in gene sequence data.

Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon-known as heterotachy-can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.

Phylogenetic mixture models can reduce node-density artifacts

We investigate the performance of phylogenetic mixture models in reducing a well-known and pervasive artifact of phylogenetic inference known as the node-density effect, comparing them to partitioned analyses of the same data. The node-density effect refers to the tendency for the amount of evolutionary change in longer branches of phylogenies to be underestimated compared to that in regions of the tree where there are more nodes and thus branches are typically shorter. Mixture models allow more than one model of sequence evolution to describe the sites in an alignment without prior knowledge of the evolutionary processes that characterize the data or how they correspond to different sites. If multiple evolutionary patterns are common in sequence evolution, mixture models may be capable of reducing node-density effects by characterizing the evolutionary processes more accurately. In gene-sequence alignments simulated to have heterogeneous patterns of evolution, we find that mixture models can reduce node-density effects to negligible levels or remove them altogether, performing as well as partitioned analyses based on the known simulated patterns. The mixture models achieve this without knowledge of the patterns that generated the data and even in some cases without specifying the full or true model of sequence evolution known to underlie the data. The latter result is especially important in real applications, as the true model of evolution is seldom known. We find the same patterns of results for two real data sets with evidence of complex patterns of sequence evolution: mixture models substantially reduced node-density effects and returned better likelihoods compared to partitioning models specifically fitted to these data. We suggest that the presence of more than one pattern of evolution in the data is a common source of error in phylogenetic inference and that mixture models can often detect these patterns even without prior knowledge of their presence in the data. Routine use of mixture models alongside other approaches to phylogenetic inference may often reveal hidden or unexpected patterns of sequence evolution and can improve phylogenetic inference.

The automation of nested clade phylogeographic analysis

ANeCA is a fully automated implementation of Nested Clade Phylogeographic Analysis. This was originally developed by Templeton and colleagues, and has been used to infer, from the pattern of gene sequence polymorphisms in a geographically structured population, the historical demographic processes that have shaped its evolution. Until now it has been necessary to perform large parts of the procedure manually. We provide a program that will take data in Nexus sequential format, and directly output a set of inferences. The software also includes TCS v1.18 and GeoDis v2.2 as part of automation.

Hespellia stercorisuis gen. nov., sp nov and Hespellia porcina sp nov., isolated from swine manure storage pits

Four Gram-positive-staining, strictly anaerobic, non-spore-forming, rod-shaped organisms were isolated from a pig manure storage pit. Comparative 16S rRNA gene sequence analysis revealed that the isolates belonged to two related but distinct groups. Sequence analysis showed that the two groups of isolates were highly related to each other (approx. 97% 16S rRNA gene sequence similarity), forming a distinct cluster within the Clostridium coccoides suprageneric rDNA grouping. Biochemical and physiological studies confirmed the division of the isolates into two related, albeit distinct, groups. Based on both phenotypic and phylogenetic evidence, it is proposed that the unidentified rod-shaped isolates from pig manure should be classified in a novel genus, Hespellia gen. nov., as Hespellia stercorisuis sp. nov. and Hespellia porcina sp. nov. The type species of the novel genus is H. stercorisuis (type strain, PC18(T) = NRRL B-23456(T) = CCUG 46279(T) = ATCC BAA-677(T)) and the type strain of H. porcina is PC80(T) (= NRRL B-23458(T) = ATCC BAA-674(T)).

Characterization of the 23S and 5S rRNA genes and 23S-5S intergenic spacer region (ITS-2) of Photobacterium damselae

The 23S ribosomal RNA (rRNA) gene has been sequenced in strains of the fish pathogens Photobacterium damselae subsp. damselae (ATCC 33539) and subsp. piscicida (ATCC 29690), showing that 3 nucleotide positions are clearly different between subspecies. In addition, the 5S rRNA gene plus the intergenic spacer region between the 23S and 5S rRNA genes (ITS-2) were amplified, cloned and sequenced for the 2 reference strains as well as the field isolates RG91 (subsp. damselae) and DI21 (subsp. piscicida). A 100% similarity was found for the consensus 5S rRNA gene sequence in the 2 subspecies, although some microheterogeneity was detected as inter-cistronic variability within the same chromosome. Sequence analysis of the spacer region between the 23S and 5S rRNA genes revealed 2 conserved and 3 variable nucleotide sequence blocks, and 4 different modular organizations were found. The ITS-2 spacer region exhibited both inter-subspecies and inter-cistronic polymorphism, with a mosaic-like structure. The EMBL accession numbers for the 23S, 5S and ITS-2 sequences are: P. damselae subsp. piscicida 5S gene (AJ274379), P. damselae subsp. damselae 23S gene (Y18520), subsp. piscicida 23S gene (Y17901), R damselae subsp. piscicida ITS-2 (AJ250695, AJ250696), P. damselae subsp. damselae ITS-2 (AJ250697, AJ250698).

Saccharospirillum impatiens gen. nov., sp nov., a novel gamma-Proteobacterium isolated from hypersaline Ekho Lake (East Antarctica)

Five Gram-negative, motile, aerobic to microaerophilic spirilla were isolated from various depths of the hypersaline, heliothermal and meromictic Ekho Lake (East Antarctica). The strains are oxidase- and catalase-positive, metabolize a variety of sugars and carboxylic acids and have an absolute requirement for sodium ions. The predominant fatty acids of the organisms are C-16: (1)omega7c, C-16:0 and C(18:1)omega7c, with C-10:1 3-OH, C-10:0 3-OH, C-12:0 3-OH, C-14:1 3-OH, C-14:0 3-OH and C-19:1 present in smaller amounts. The main polar lipids are diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol and phosphatidylmonomethylamine. The DNA base composition of the strains is 54-55 mol% G + C. 16S rRNA gene sequence comparisons show that the isolates are related to the genera Oceanospirillum, Pseudospirillum, Marinospirillum, Halomonas and Chromohalobacter in the gamma-Proteobacteria. Morphological, physiological and genotypic differences from these previously described genera support the description of a novel genus and species, Saccharospirillum impatiens gen. nov., sp. nov. The type strain is EL-105(T) (= DSM 12546(T) = CECT 5721(T)).

Roseisalinus antarcticus gen. nov., sp nov., a novel aerobic bacteriochlorophyll a-producing alpha-proteobacterium isolated from hypersaline Ekho Lake, Antarctica

A Gram-negative, aerobic to microaerophilic rod was isolated from 10 m depths of the hypersaline, heliothermal and meromictic Ekho Lake (East Antarctica). The strain was oxidase- and catalase-positive, metabolized a variety of carboxylic acids and sugars and produced lipase. Cells had an absolute requirement for artificial sea water, which could not be replaced by NaCl. A large in vivo absorption band at 870 nm indicated production of bacteriochlorophyll a. The predominant fatty acids of this organism were 16:0 and 18:1omega7c, with 3-OH 10:0, 16:1omega7c and 18:0 in lower amounts. The main polar lipids were diphosphatidylglycerol, phosphatidylglycerol and phosphatidylcholine. Ubiquinone 10 was produced. The DNA G + C content was 67 mol%. 16S rRNA gene sequence comparisons indicated that the isolate represents a member of the Roseobacter clade within the alpha-Proteobacteria. The organism showed no particular relationship to any members of this clade but clustered on the periphery of the genera Jannaschia, Octadecabacter and 'Marinosulfonomonas' and the species Ruegeria gelatinovorans. Distinct morphological, physiological and genotypic differences to these previously described taxa supported the description of a new genus and a novel species, for which the name Roseisalinus antarcticus gen. nov., sp. nov. is proposed. The type strain is EL-88(T) (= DSM 11466(T) = CECT 7023(T)).

Actinobaculum urinale sp. nov., from human urine

A hitherto undescribed Actinomyces-like bacterium was isolated from human urine. Based on its biochemical characteristics, the unidentified bacterium did not correspond to any currently described Actinomyces species or related taxa. Comparative 16S rRNA gene sequencing showed that the unknown bacterium exhibits a specific phylogenetic association with the genus Actinobaculum, but a sequence divergence of > 5% from the two currently recognized members of this genus, Actinobaculum schaalii and Actinobaculum suis, demonstrates that it represents a distinct species. Based on both phenotypic and 16S rRNA gene sequence considerations, it is proposed that the unknown bacterium from urine should be classified as a novel species, Actinobaculum urinale sp. nov. The type strain of Actinobaculum urinale is CCUG 46093(T) (= CIP 107424(T)).

Clostridium colicanis sp. nov., from canine faeces

Morphological, biochemical and molecular genetic studies were performed on an unknown, anaerobic, rod-shaped organism isolated from faeces of a canine. The organism was tentatively identified as a member of the genus Clostridium based on its cellular morphology and ability to form endospores but, biochemically, it did not appear to correspond to any recognized species of this genus. Comparative 16S rRNA gene sequence analysis showed that the bacterium represents a previously unrecognized subline within Clostridium rRNA group I (Clostridium sensu stricto), which includes Clostridium butyricum, the type species of the genus. The nearest phylogenetic relatives of the unknown bacterium corresponded to Clostridium absonum, Clostridium baratii, Eubacterium budayi, Eubacterium moniliforme, Eubacterium multiforme and Eubacterium nitritogenes, but 16S rRNA sequence divergence values of > 3% demonstrated that it represents a novel species. Based on the findings presented, a novel species, Clostridium colicanis sp. nov., is described, with the type strain 3WC2(T) (=CCUG 44556(T) =DSM 13634(T)).

Corynebacterium sphenisci sp nov., isolated from wild penguins

Six unidentified Gram-positive, rod-shaped organisms recovered from the cloacae of apparently healthy wild penguins were characterized by phenotypic and molecular taxonomic methods. Chemotaxonomic investigations revealed the presence of a cell wall based on meso-diaminopimelic acid and long-chain cellular fatty acids of the straight-chain saturated and monounsaturated types, consistent with the genus Corynebacterium. Corynomycolic acids, which are characteristic of the genus, were also detected, albeit in small amounts. Comparative 16S rRNA gene sequencing studies showed that the unidentified organisms were phylogenetically related to corynebacteria and represent a novel subline associated with a small subcluster of species that includes Corynebacterium xerosis, Corynebacterium amycolatum and Corynebacterium freneyi. The unknown isolates were readily distinguished from their closest phylogenetic relatives and all other Corynebacterium species with validly published names by using a combination of biochemical and chemotaxonomic criteria. Based on both phenotypic and 16S rRNA gene sequence considerations, it is proposed that the unknown isolates recovered from penguins be classified as a novel species in the genus Corynebacterium, Corynebacterium sphenisci sp. nov. The type strain is CECT 5990(T) (= CCUG 46398(T)).