Identification of unusual Campylobacter-like isolates from poultry products as Helicobacter pullorum

Twenty-six unclassified Campylobacter -like strains previously isolated from 15 chicken carcasses and caecal contents, together with two more strains isolated from chicken faeces on a different occasion, were identified as Helicobacter pullorum using various phenotypic identification methods. API Campy identification kits and a 16-test identification scheme developed for campylobacters failed to identify these bacteria, or identified them as Campylobacter spp. Eighteen strains (including the two isolated on a different occasion) were chosen for examination using a more comprehensive probabilistic identification scheme. Using this method, 14 of the 18 strains were identified as H. pullorum with ID scores >95% ; two strains were also identified as H. pullorum with lower ID scores. Of the remaining two strains, one was not identified with this scheme and the other was misidentified to the H. acinonyx pylori complex. Whole cell protein profiling by SDS-PAGE confirmed the identity of these isolates as H. pullorum , affirming the value of a polyphasic approach for accurately identifying campylobacteria. The comparatively high prevalence of H. pullorum in poultry determined in this study (60%) suggests that routine isolation and identification methods should be amended to enable a thorough evaluation of its role in human gastroenteritis and avian hepatitis. Some phenotypic characters useful in identifying poultry campylobacteria are presented which could be utilized, along with other technique(s), for improved differentiation of the campylobacteria that are found in poultry.     

Species diversity of campylobacteria in a wild bird community in Sweden

AIMS: To analyse the occurrence and host species distribution of campylobacteria species in shorebirds, geese and cattle on grazed coastal meadows in Sweden. METHODS AND RESULTS: Species identification was performed through a polyphasic approach, incorporating Amplified Fragment Length Polymorphism (AFLP) profiling, 16S RNA gene sequence analysis together with extensive phenotypic characterization. From 247 sampled birds and 71 cattle, we retrieved 113 urease positive thermophilic Campylobacter (UPTC) and 16 Campylobacter jejuni ssp. jejuni isolates. Furthermore, 18 isolates of Helicobacter canadensis, and five isolates that potentially represent a new genus of micro-aerophilic, spiral and Gram-negative bacteria were isolated. The distribution of bacterial species on hosts was uneven: all H. canadensis isolates were retrieved from geese, while all but one of the Campylobacter lari UPTC isolates were found in shorebirds. AFLP type distribution of Camp. lari UPTC isolates among individual, resampled and breeding-paired Redshank birds generally indicated a constant shift in strain populations over time and absence of geographical clustering. CONCLUSIONS: The large number of isolated campylobacteria, including species that are zoonotic enteropathogens, indicates that these wild birds potentially may serve as reservoirs of human infections. However, despite a common environment, the different host species largely carried their own campylobacteria populations, indicating that cross-species transmission is rare. SIGNIFICANCE AND IMPACT OF THE STUDY: Our study is one of few that provide data on the occurrence of campylobacteria in wild animals, adding information on the ecology and epidemiology of micro-organisms that are of public health concern.     

Evolution of the courtship phenotype in the bird of paradise genus Parotia (Aves: Paradisaeidae): homology, phylogeny, and modularity

Birds of paradise (Aves: Paradisaeidae) exhibit extreme differences among taxa in courtship-related form (i.e. courtship phenotype). In the genus Parotia , the courtship phenotype is organizationally modular and this property may play an important role in the evolution of phenotypic disparity among taxa. The present study investigates variational aspects of phenotypic modularity in the Parotia by examining the structure and composition of courtship form in a comparative context. First, a module-based model of male display-phenotypes is compiled for four biological species to facilitate phenotypic comparison. Models are constructed using data from existing phenotype ontologies and associated video-vouchers. Next, a phylogenetic analysis of display-phenotype data is performed using a matrix of 47 etho-phenotypic characters coded for eight Parotia and out-group taxa. Analysis yields one tree, length 60 (CI = 0.83; RI = 0.85). The results demonstrate variation among taxa to be greater at higher-levels of phenotypic integration (i.e. among display-modules) than at intermediate and lower-levels (i.e. among phase- and element-modules). Three display-modules and five of six phase-modules were present in the common ancestor and complexity has increased through time as the display-modules became dissociated into subunits that diverged independently. The history of Parotia evolution involves numerous instances of duplication and divergence of etho-phenotypic modular components and likely reflects the same processes that have contributed to the pronounced phenotypic disparity within the entire bird of paradise radiation.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 491–504.     

Numerical classification and identification of Bacillus sphaericus including some strains pathogenic for mosquito larvae

Ninety-one strains of Bacillus sphaericus, including representatives of all the established DNA homology groups, related round-spored and oval-spored species, and six strains pathogenic for mosquito larvae, were examined for 155 characters. Numerical analyses (Jaccard coefficient/average linkage clustering) based on the 88 variable features revealed 14 clusters at the 79% similarity level that contained more than one strain and 17 single member clusters. All insect pathogenic strains were recovered in a single cluster and the classification was in accord with an established classification based on DNA sequence homology. Two frequency matrices for probabilistic identification were constructed and tested. A comprehensive matrix comprising 14 mesophilic, round-spored taxa and 27 tests gave good results for identification of hypothetical median organisms, cluster overlap and identifications of representative strains (based on data generated in the classification study). Reference strains for the 14 taxa and eight additional insect pathogenic strains were examined for the 27 tests and were correctly identified with high scores using this matrix. A second matrix comprising seven taxa and 13 tests also performed well in the theoretical evaluation and correctly identified the reference strains and insect pathogenic strains.     

Dirichlet multinomial mixtures: generative models for microbial metagenomics

We introduce Dirichlet multinomial mixtures (DMM) for the probabilistic modelling of microbial metagenomics data. This data can be represented as a frequency matrix giving the number of times each taxa is observed in each sample. The samples have different size, and the matrix is sparse, as communities are diverse and skewed to rare taxa. Most methods used previously to classify or cluster samples have ignored these features. We describe each community by a vector of taxa probabilities. These vectors are generated from one of a finite number of Dirichlet mixture components each with different hyperparameters. Observed samples are generated through multinomial sampling. The mixture components cluster communities into distinct 'metacommunities', and, hence, determine envirotypes or enterotypes, groups of communities with a similar composition. The model can also deduce the impact of a treatment and be used for classification. We wrote software for the fitting of DMM models using the 'evidence framework' (http://code.google.com/p/microbedmm/). This includes the Laplace approximation of the model evidence. We applied the DMM model to human gut microbe genera frequencies from Obese and Lean twins. From the model evidence four clusters fit this data best. Two clusters were dominated by Bacteroides and were homogenous; two had a more variable community composition. We could not find a significant impact of body mass on community structure. However, Obese twins were more likely to derive from the high variance clusters. We propose that obesity is not associated with a distinct microbiota but increases the chance that an individual derives from a disturbed enterotype. This is an example of the 'Anna Karenina principle (AKP)' applied to microbial communities: disturbed states having many more configurations than undisturbed. We verify this by showing that in a study of inflammatory bowel disease (IBD) phenotypes, ileal Crohn's disease (ICD) is associated with a more variable community.     

The role of genetic constraints on the diversification of Iberian taxa of the genus Aquilegia (Ranunculaceae)

The microevolutionary process of adaptive phenotypic differentiation of quantitative traits between populations or closely‐related taxa depends on the response of populations to the action of natural selection. However, this response can be constrained by the structure of the matrix of additive genetic variance and covariance between traits in each population ( G matrix). In the present study, we obtained additive genetic variance and narrow sense heritability for 25 floral and vegetative traits of three subspecies of Aquilegia vulgaris, and one subspecies of Aquilegia pyrenaica through a common garden crossing experiment. For two vegetative and one floral trait, we also obtained the G matrix and genetic correlations between traits in each subspecies. The amount of genetic variation available in wild populations is not responsible for the larger differentiation of vegetative than floral traits found in this group of columbines. However, the low heritability of some traits constrained their evolution because phenotypic variability among taxa was larger for traits with larger heritability. We confirmed that the process of diversification of the studied taxa involved shifts in the G matrix, mainly determined by changes in the genetic covariance between floral and vegetative traits, probably caused by linkage disequilibrium in narrow endemic taxa. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 252–261.     

Phylogenetic analysis of phenotypic characters of Tunicata supports basal Appendicularia and monophyletic Ascidiacea

With approximately 3000 marine species, Tunicata represents the most disparate subtaxon of Chordata. Molecular phylogenetic studies support Tunicata as sister taxon to Craniota, rendering it pivotal to understanding craniate evolution. Although successively more molecular data have become available to resolve internal tunicate phylogenetic relationships, phenotypic data have not been utilized consistently. Herein these shortcomings are addressed by cladistically analyzing 117 phenotypic characters for 49 tunicate species comprising all higher tunicate taxa, and five craniate and cephalochordate outgroup species. In addition, a combined analysis of the phenotypic characters with 18S rDNA-sequence data is performed in 32 OTUs. The analysis of the combined data is congruent with published molecular analyses. Successively up-weighting phenotypic characters indicates that phenotypic data contribute disproportionally more to the resulting phylogenetic hypothesis. The strict consensus tree from the analysis of the phenotypic characters as well as the single most parsimonious tree found in the analysis of the combined dataset recover monophyletic Appendicularia as sister taxon to the remaining tunicate taxa. Thus, both datasets support the hypothesis that the last common ancestor of Tunicata was free-living and that ascidian sessility is a derived trait within Tunicata. “Thaliacea” is found to be paraphyletic with Pyrosomatida as sister taxon to monophyletic Ascidiacea and the relationship between Doliolida and Salpida is unresolved in the analysis of morphological characters; however, the analysis of the combined data reconstructs Thaliacea as monophyletic nested within paraphyletic “Ascidiacea”. Therefore, both datasets differ in the interpretation of the evolution of the complex holoplanktonic life history of thaliacean taxa. According to the phenotypic data, this evolution occurred in the plankton, whereas from the combined dataset a secondary transition into the plankton from a sessile ascidian is inferred. Besides these major differences, both analyses are in accord on many phylogenetic groupings, although both phylogenetic reconstructions invoke a high degree of homoplasy. In conclusion, this study represents the first serious attempt to utilize the potential phylogenetic information present in phenotypic characters to elucidate the inter-relationships of this diverse marine taxon in a consistent cladistic framework.     

Development of a nutrient medium to culture microorganisms of the genus Campylobacter

The nutrient media developed in the USSR and CSR for the isolation of campylobacteria are multi-purpose as they enable the isolation of campylobacteria not only from material derived from patients with acute intestinal infection, but also from the internal organs of healthy and diseased animals and poultry as well as from a variety of objects in the environment.     

Probabilistic methods surpass parsimony when assessing clade support in phylogenetic analyses of discrete morphological data

Fossil taxa are critical to inferences of historical diversity and the origins of modern biodiversity, but realizing their evolutionary significance is contingent on restoring fossil species to their correct position within the tree of life. For most fossil species, morphology is the only source of data for phylogenetic inference; this has traditionally been analysed using parsimony, the predominance of which is currently challenged by the development of probabilistic models that achieve greater phylogenetic accuracy. Here, based on simulated and empirical datasets, we explore the relative efficacy of competing phylogenetic methods in terms of clade support. We characterize clade support using bootstrapping for parsimony and Maximum Likelihood, and intrinsic Bayesian posterior probabilities, collapsing branches that exhibit less than 50% support. Ignoring node support, Bayesian inference is the most accurate method in estimating the tree used to simulate the data. After assessing clade support, Bayesian and Maximum Likelihood exhibit comparable levels of accuracy, and parsimony remains the least accurate method. However, Maximum Likelihood is less precise than Bayesian phylogeny estimation, and Bayesian inference recaptures more correct nodes with higher support compared to all other methods, including Maximum Likelihood. We assess the effects of these findings on empirical phylogenies. Our results indicate probabilistic methods should be favoured over parsimony.     

Morphological and molecular phylogenetic context of the angiosperms: contrasting the 'top-down' and 'bottom-up' approaches used to infer the likely characteristics of the first flowers

Recent attempts to address the long-debated 'origin' of the angiosperms depend on a phylogenetic framework derived from a matrix of taxa versus characters; most assume that empirical rigour is proportional to the size of the matrix. Sequence-based genotypic approaches increase the number of characters (nucleotides and indels) in the matrix but are confined to the highly restricted spectrum of extant species, whereas morphology-based approaches increase the number of phylogenetically informative taxa (including fossils) at the expense of accessing only a restricted spectrum of phenotypic characters. The two approaches are currently delivering strongly contrasting hypotheses of relationship. Most molecular studies indicate that all extant gymnosperms form a natural group, suggesting surprisingly early divergence of the lineage that led to angiosperms, whereas morphology-only phylogenies indicate that a succession of (mostly extinct) gymnosperms preceded a later angiosperm origin. Causes of this conflict include: (i) the vast phenotypic and genotypic lacuna, largely reflecting pre-Cenozoic extinctions, that separates early-divergent living angiosperms from their closest relatives among the living gymnosperms; (ii) profound uncertainty regarding which (a) extant and (b) extinct angiosperms are most closely related to gymnosperms; and (iii) profound uncertainty regarding which (a) extant and (b) extinct gymnosperms are most closely related to angiosperms, and thus best serve as 'outgroups' dictating the perceived evolutionary polarity of character transitions among the early-divergent angiosperms. These factors still permit a remarkable range of contrasting, yet credible, hypotheses regarding the order of acquisition of the many phenotypic characters, reproductive and vegetative, that distinguish 'classic' angiospermy from 'classic' gymnospermy. The flower remains ill-defined and its mode (or modes) of origin remains hotly disputed; some definitions and hypotheses of evolutionary relationships preclude a role for the flower in delimiting the angiosperms. We advocate maintenance of parallel, reciprocally illuminating programmes of morphological and molecular phylogeny reconstruction, respectively supported by homology testing through additional taxa (especially fossils) and evolutionary-developmental genetic studies that explore genes potentially responsible for major phenotypic transitions.     

Phylogeny of extant and fossil Juglandaceae inferred from the integration of molecular and morphological data sets

It is widely acknowledged that integrating fossils into data sets of extant taxa is imperative for proper placement of fossils, resolution of relationships, and a better understanding of character evolution. The importance of this process has been further magnified because of the crucial role of fossils in dating divergence times. Outstanding issues remain, including appropriate methods to place fossils in phylogenetic trees, the importance of molecules versus morphology in these analyses, as well as the impact of potentially large amounts of missing data for fossil taxa. In this study we used the angiosperm clade Juglandaceae as a model for investigating methods of integrating fossils into a phylogenetic framework of extant taxa. The clade has a rich fossil record relative to low extant diversity, as well as a robust molecular phylogeny and morphological database for extant taxa. After combining fossil organ genera into composite and terminal taxa, our objectives were to (1) compare multiple methods for the integration of the fossils and extant taxa (including total evidence, molecular scaffolds, and molecular matrix representation with parsimony [MRP]); (2) explore the impact of missing data (incomplete taxa and characters) and the evidence for placing fossils on the topology; (3) simulate the phylogenetic effect of missing data by creating "artificial fossils"; and (4) place fossils and compare the impact of single and multiple fossil constraints in estimating the age of clades. Despite large and variable amounts of missing data, each of the methods provided reasonable placement of both fossils and simulated "artificial fossils" in the phylogeny previously inferred only from extant taxa. Our results clearly show that the amount of missing data in any given taxon is not by itself an operational guideline for excluding fossils from analysis. Three fossil taxa (Cruciptera simsonii, Paleoplatycarya wingii, and Platycarya americana) were placed within crown clades containing living taxa for which relationships previously had been suggested based on morphology, whereas Polyptera manningii, a mosaic taxon with equivocal affinities, was placed firmly as sister to two modern crown clades. The position of Paleooreomunnea stoneana was ambiguous with total evidence but conclusive with DNA scaffolds and MRP. There was less disturbance of relationships among extant taxa using a total evidence approach, and the DNA scaffold approach did not provide improved resolution or internal support for clades compared to total evidence, whereas weighted MRP retained comparable levels of support but lost crown clade resolution. Multiple internal minimum age constraints generally provided reasonable age estimates, but the use of single constraints provided by extinct genera tended to underestimate clade ages.     

Discordant 16S and 23S rRNA gene phylogenies for the genus Helicobacter: implications for phylogenetic inference and systematics

Analysis of 16S rRNA gene sequences has become the primary method for determining prokaryotic phylogeny. Phylogeny is currently the basis for prokaryotic systematics. Therefore, the validity of 16S rRNA gene-based phylogenetic analyses is of fundamental importance for prokaryotic systematics. Discrepancies between 16S rRNA gene analyses and DNA-DNA hybridization and phenotypic analyses have been noted in the genus Helicobacter. To clarify these discrepancies, we sequenced the 23S rRNA genes for 55 helicobacter strains representing 41 taxa (>2,700 bases per sequence). Phylogenetic-tree construction using neighbor-joining, parsimony, and maximum likelihood methods for 23S rRNA gene sequence data yielded stable trees which were consistent with other phenotypic and genotypic methods. The 16S rRNA gene sequence-derived trees were discordant with the 23S rRNA gene trees and other data. Discrepant 16S rRNA gene sequence data for the helicobacters are consistent with the horizontal transfer of 16S rRNA gene fragments and the creation of mosaic molecules with loss of phylogenetic information. These results suggest that taxonomic decisions must be supported by other phylogenetically informative macromolecules, such as the 23S rRNA gene, when 16S rRNA gene-derived phylogeny is discordant with other credible phenotypic and genotypic methods. This study found Wolinella succinogenes to branch with the unsheathed-flagellum cluster of helicobacters by 23S rRNA gene analyses and whole-genome comparisons. This study also found intervening sequences (IVSs) in the 23S rRNA genes of strains of 12 Helicobacter species. IVSs were found in helices 10, 25, and 45, as well as between helices 31' and 27'. Simultaneous insertion of IVSs at three sites was found in H. mesocricetorum.     

The probability of genetic parallelism and convergence in natural populations

Genomic and genetic methods allow investigation of how frequently the same genes are used by different populations during adaptive evolution, yielding insights into the predictability of evolution at the genetic level. We estimated the probability of gene reuse in parallel and convergent phenotypic evolution in nature using data from published studies. The estimates are surprisingly high, with mean probabilities of 0.32 for genetic mapping studies and 0.55 for candidate gene studies. The probability declines with increasing age of the common ancestor of compared taxa, from about 0.8 for young nodes to 0.1–0.4 for the oldest nodes in our study. Probability of gene reuse is higher when populations begin from the same ancestor (genetic parallelism) than when they begin from divergent ancestors (genetic convergence). Our estimates are broadly consistent with genomic estimates of gene reuse during repeated adaptation to similar environments, but most genomic studies lack data on phenotypic traits affected. Frequent reuse of the same genes during repeated phenotypic evolution suggests that strong biases and constraints affect adaptive evolution, resulting in changes at a relatively small subset of available genes. Declines in the probability of gene reuse with increasing age suggest that these biases diverge with time.     

Phenotypic Variation in Infants,Not Adults,Reflects Genotypic Variation among Chimpanzees and Bonobos

Studies comparing phenotypic variation with neutral genetic variation in modern humans have shown that genetic drift is a main factor of evolutionary diversification among populations. The genetic population history of our closest living relatives, the chimpanzees and bonobos, is now equally well documented, but phenotypic variation among these taxa remains relatively unexplored, and phenotype-genotype correlations are not yet documented. Also, while the adult phenotype is typically used as a reference, it remains to be investigated how phenotype-genotye correlations change during development. Here we address these questions by analyzing phenotypic evolutionary and developmental diversification in the species and subspecies of the genus Pan. Our analyses focus on the morphology of the femoral diaphysis, which represents a functionally constrained element of the locomotor system. Results show that during infancy phenotypic distances between taxa are largely congruent with non-coding (neutral) genotypic distances. Later during ontogeny, however, phenotypic distances deviate from genotypic distances, mainly as an effect of heterochronic shifts between taxon-specific developmental programs. Early phenotypic differences between Pan taxa are thus likely brought about by genetic drift while late differences reflect taxon-specific adaptations.     

Positive correlation between diversification rates and phenotypic evolvability can mimic punctuated equilibrium on molecular phylogenies

The hypothesis of punctuated equilibrium proposes that most phenotypic evolution occurs in rapid bursts associated with speciation events. Several methods have been developed that can infer punctuated equilibrium from molecular phylogenies in the absence of paleontological data. These methods essentially test whether the variance in phenotypes among extant species is better explained by evolutionary time since common ancestry or by the number of estimated speciation events separating taxa. However, apparent "punctuational" trait change can be recovered on molecular phylogenies if the rate of phenotypic evolution is correlated with the rate of speciation. Strong support for punctuational models can arise even if the underlying mode of trait evolution is strictly gradual, so long as rates of speciation and trait evolution covary across the branches of phylogenetic trees, and provided that lineages vary in their rate of speciation. Species selection for accelerated rates of ecological or phenotypic divergence can potentially lead to the perception that most trait divergence occurs in association with speciation events.     

Bayesian species delimitation combining multiple genes and traits in a unified framework

Delimitation of species based exclusively on genetic data has been advocated despite a critical knowledge gap: how might such approaches fail because they rely on genetic data alone, and would their accuracy be improved by using multiple data types. We provide here the requisite framework for addressing these key questions. Because both phenotypic and molecular data can be analyzed in a common Bayesian framework with our program iBPP, we can compare the accuracy of delimited taxa based on genetic data alone versus when integrated with phenotypic data. We can also evaluate how the integration of phenotypic data might improve species delimitation when divergence occurs with gene flow and/or is selectively driven. These two realities of the speciation process are ignored by currently available genetic approaches. Our model accommodates phenotypic characters that exhibit different degrees of divergence, allowing for both neutral traits and traits under selection. We found a greater accuracy of estimated species boundaries with the integration of phenotypic and genetic data, with a strong beneficial influence of phenotypic data from traits under selection when the speciation process involves gene flow. Our results highlight the benefits of multiple data types, but also draws into question the rationale of species delimitation based exclusively on genetic data.     

Seriation in paleontological data using markov chain Monte Carlo methods

Given a collection of fossil sites with data about the taxa that occur in each site, the task in biochronology is to find good estimates for the ages or ordering of sites. We describe a full probabilistic model for fossil data. The parameters of the model are natural: the ordering of the sites, the origination and extinction times for each taxon, and the probabilities of different types of errors. We show that the posterior distributions of these parameters can be estimated reliably by using Markov chain Monte Carlo techniques. The posterior distributions of the model parameters can be used to answer many different questions about the data, including seriation (finding the best ordering of the sites) and outlier detection. We demonstrate the usefulness of the model and estimation method on synthetic data and on real data on large late Cenozoic mammals. As an example, for the sites with large number of occurrences of common genera, our methods give orderings, whose correlation with geochronologic ages is 0.95.     

Supraspecific taxa as terminals in cladistic analysis: implicit assumptions of monophyly and a comparison of methods

The use of supraspecific terminal taxa to represent groups of species in phylogenetic analyses can result in changes to inferred relationships as compared to a complete species level analysis. These changes in topology result from interactions among (1) the cladistic status of the supraspecific taxa; (2) the method used to represent the taxa as single terminals, and (3) incongruence in the data set. We examine the effects of using supraspecific terminal taxa using a parallel analysis of hypothetical examples and an actual data matrix for the true seals (Mammalia: Phocidae). Incongruence among characters can produce changes in topology by shifting the ‘balance of power’ among groups of characters when supraspecific taxa are represented as single terminals. In the absence of homoplasy, the correct topology is maintained. Of the three methods for representing supraspecific taxa, the ‘ancestral’ method, which explicidy infers the common ancestor of the group corresponding to the taxon, performed the best, always maintaining the correct topology when monophyletic taxa were represented. This agrees with theoretical predictions. The ‘democratic’ and ‘exemplar’ methods, which represent the higher level taxon through a survey of all or one of its extant constituent species, respectively, were not as effective in maintaining the correct topology. Although both occasionally provided correct answers, their occurrences were largely unpredictable. The success of the exemplar method varies with the species selected. The simultaneous representation of two or more higher level taxa produced interactive effects where the resultant topology included different clades than when the taxa were collapsed individually. Interactive effects occurred with all three methods, albeit to a lesser degree for the ancestral method. Changes in topology were observed regardless of whether the higher group was monophyletic or not, but were more prevalent when it was paraphyletic. Unfortunately, there does not seem to be a reliable way to determine when a paraphyletic group has been included in the analysis (e.g. through bootstrap values or indices measuring homoplasy). The implications of these findings for phylogenetic analyses of molecular data are also discussed.