Phenotypic and genetic characterization of Pasteurella multocida and related isolates from rabbits in Switzerland

Several bacteria belonging to the family Pasteurellaceae are potential pathogens in rabbits. In particular, Pasteurella multocida is considered to be important, and outbreaks caused by this species result in considerable economic losses in rabbitries. However, Pasteurellaceae spp. isolated from rabbits are poorly characterized, and thus, proper identification of P. multocida isolates from these animals is problematic and often unsatisfactory, thereby hampering epidemiological investigations. Therefore, 228 isolates from rabbit populations originating from a breeding and fattening organization with group management and postmortem cases with pasteurellosis from individual owners were phenotypically and genotypically analyzed using biochemical tests and repetitive extragenic palindromic polymerase chain reaction (REP-PCR). Furthermore, 41 samples representing observed phenotypes were selected for phylogenetic analysis using 16S ribosomal RNA and rpoB genes. The REP-PCR typing and phylogenetic analyses correlated well and appeared to be distinct molecular methods for characterization of rabbit isolates. Phenotyping, however, diverged from molecular recognition, reflecting the problematic conventional diagnosis of these strains. The fermentation of sorbitol appeared to be an imprecise indicator for P. multocida subspecies classification. According to REP-PCR and sequencing results, 82% of the isolates were characterized as P. multocida subsp. multocida, 3% as P. multocida subsp. septica, and 5% as P. multocida. Further, 5% were identified as Pasteurella canis. The other 5% represented a homogeneous group of unknown species belonging to the Pasteurellaceae. Samples obtained from individual postmortem cases demonstrated a higher phenotypic and genetic heterogeneity than samples from group management rabbits.

Routine phenotypic identification of bacterial species of the family Pasteurellaceae isolated from animals

Pasteurellaceae are bacteria with an important role as primary or opportunistic, mainly respiratory, pathogens in domestic and wild animals. Some species of Pasteurellaceae cause severe diseases with high economic losses in commercial animal husbandry and are of great diagnostic concern. Because of new data on the phylogeny of Pasteurellaceae, their taxonomy has recently been revised profoundly, thus requiring an improved phenotypic differentiation procedure to identify the individual species of this family. A new and simplified procedure to identify species of Actinobacillus, Avibacterium, Gallibacterium, Haemophilus, Mannheimia, Nicoletella, and Pasteurella, which are most commonly isolated from clinical samples of diseased animals in veterinary diagnostic laboratories, is presented in the current study. The identification procedure was evaluated with 40 type and reference strains and with 267 strains from routine diagnostic analysis of various animal species, including 28 different bacterial species. Type, reference, and field strains were analyzed by 16S ribosomal RNA (rrs) and rpoB gene sequencing for unambiguous species determination as a basis to evaluate the phenotypic differentiation schema. Primary phenotypic differentiation is based on beta-nicotinamide adenine dinucleotide (beta-NAD) dependence and hemolysis, which are readily determined on the isolation medium. The procedure divides the 28 species into 4 groups for which particular biochemical reactions were chosen to identify the bacterial species. The phenotypic identification procedure allowed researchers to determine the species of 240 out of 267 field strains. The procedure is an easy and cost-effective system for the rapid identification of species of the Pasteurellaceae family isolated from clinical specimens of animals.

Laribacter hongkongensis: a potential cause of infectious diarrhea

In this study, we describe the isolation of Laribacter hongkongensis, a recently described genus and species of bacterium, in pure culture on charcoal cefoperazone deoxycholate agar from the stool of six patients with diarrhea. Three patients were residents of Hong Kong, and three of Switzerland. In none of the stool samples obtained from these six patients was Salmonella, Shigella, enterohemorrhagic Escherichia coli, Vibrio, Aeromonas, Plesiomonas, or Campylobacter recovered. Rotavirus antigen detection, electron microscopic examination for viruses, and microscopic examinations for ova and cysts were all negative for the stool samples obtained from the three patients in Hong Kong. Enterotoxigenic E. coli was recovered from one of the patients in Hong Kong. Unlike L. hongkongensis type strain HKU1, all the six strains were motile with bipolar flagellae. Sequencing of the 16S ribosomal RNA genes of the six strains showed that they all had sequences with only 0-2 base differences to that of the type strain. Pulsed field gel electrophoresis of the SpeI digested genomic DNA of the six isolates and that of the type strain revealed that the seven isolates were genotypically unrelated strains. More extensive epidemiologic studies should be carried out to ascertain the causative association between L. hongkongensis and diarrhea and to define the reservoir and modes of transmission of L. hongkongensis.

Tularemia in a common marmoset (Callithrix jacchus) diagnosed by 16S rRNA sequencing

We report a case of tularemia in a common marmoset (Callithrix jacchus) diagnosed by determination of the isolate's 16S ribosomal RNA (rRNA) gene sequence. Pathological examination of the animal revealed a multifocal acute necrotizing hepatitis, interstitial nephritis, splenitis, and lymphangitis of the mandibular, retropharyngeal, and cervical and mesenteric lymph nodes. Moreover, multiple foci of acute necrosis were found in the epithelium of the jejunum and the interstitium of the lung. Bacteriological investigations revealed a septicemia. The isolated infectious agent was uncommon, not routinely diagnosed in our laboratory and therefore difficult to identify by conventional tools in a reasonable time and effort. thus, we decided to perform a genetic analysis based on the 16S rRNA gene sequence. Thereby, an infection with Francisella tularensis, the causative agent of tularemia, was unambiguously diagnosed. This shows the great advantage 16S rRNA gene sequencing has as a general identification approach for unusual or rare isolates.

Control of chromosomal rearrangements in {\it Escherichia coli\/}

A complete physical map of Escherichia coli K-12 strain MG1655 was constructed by digesting chromosomal DNA with the infrequently cutting restriction enzymes NotI, SfiI and XbaI and separating the fragments by pulsed field gel electrophoresis. The map was used to compare six K-12 strains of E. coli. Although several differences were noted and localized, the map of MG1655 was representative of all the K-12 strains tested. The maps were also used to analyze chromosomal rearrangements in the E. coli strain MG1655. The spontaneous and UV induced frequencies of tandem duplication formation were measured at several loci distributed around the chromosome. The spontaneous duplication frequency varied from 10$\sp{-5}$ to 10$\sp{-3}$ and increased at least ten-fold following mild UV irradiation treatment. Duplications of several regions of the chromosome, including the serA region and the metE region, were mapped using pulsed field gel electrophoresis. Duplications of serA were found to be large, ranging in size from 600 kb to 2100 kb. Several of the duplications isolated at serA were caused by ectopic recombination between IS5 elements and between IS186 elements. Duplications of the metE region, however, were almost exclusively the result of ectopic recombination between ribosomal RNA cistrons. Duplication frequencies were determined at both serA and metE in wild type and mismatch repair mutant strains (mutL, mutS, uvrD and recF). Even though all of the mismatch repair mutations increased duplication frequency of metE, the largest increases were observed in the mutL and mutS strains. Duplication frequency of serA was increased less dramatically by mutations in mismatch repair. Several duplications of metE isolated in a wild type and a mismatch repair mutant were mapped. The results showed that the same repeated sequences were used for duplication formation in the mismatch repair mutant as were used in the wild type strain. Several isolates showed evidence of multiple rearrangements indicating that mismatch repair may play a role in stabilizing the genome by controlling chromosomal rearrangement. ^

Subgingival microbiota of Sri Lankan tea labourers naïve to oral hygiene measures.

AIM To characterize the subgingival microbiota within a cohort of adult males (n = 32) naïve to oral hygiene practices, and to compare the composition of bacterial taxa present in periodontal sites with various probing depths. MATERIAL AND METHODS Subgingival plaque samples were collected from single shallow pocket [pocket probing depth (PPD)≤3 mm] and deep pocket (PPD≥6 mm) sites from each subject. A polymerase chain reaction based strategy was used to construct a clone library of 16S ribosomal RNA (rRNA) genes for each site. The sequences of ca. 30-60 plasmid clones were determined for each site to identify resident taxa. Microbial composition was compared using a variety of statistical and bioinformatics approaches. RESULTS A total of 1887 cloned 16S rRNA gene sequences were analysed, which were assigned to 318 operational taxonomic units (98% identity cut-off). The subgingival microbiota was dominated by Firmicutes (69.8%), Proteobacteria (16.3%), and Fusobacteria (8.0%). The overall composition of microbial communities in shallow sites was significantly different from those within deep sites (∫-Libshuff, p < 0.001). CONCLUSIONS A taxonomically diverse subgingival microbiota was present within this cohort; however, the structures of the microbial communities present in the respective subjects exhibited limited variation. Deep and shallow sites contained notably different microbial compositions, but this was not correlated with the rate of periodontal progression.

Identification of animal Pasteurellaceae by MALDI-TOF mass spectrometry

Species of the family Pasteurellaceae play an important role as primary or opportunistic, predominantly respiratory, pathogens in domestic and wild animals. Some of them cause severe disease with high economic losses in commercial animal husbandry. Hence, rapid and accurate differentiation of Pasteurellaceae is important and signifies a particular challenge to diagnostic laboratories. Identification and differentiation of Pasteurellaceae is mostly done using phenotypic tests or genetic identification based on sequence similarity of housekeeping genes, such as the rrs gene encoding the 16S ribosomal RNA (16S rRNA). Both approaches are time consuming, laborious, and costly, therefore often delaying the final diagnosis of disease or epidemics. Matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry represents an alternative rapid and reliable method for the differentiation of most members of the family Pasteurellaceae. It is able to differentiate within a few minutes the currently known 18 genera and most of the over 60 species and subspecies of Pasteurellaceae including many members encountered in veterinary diagnostic laboratories. A few closely related species and subspecies that cannot be discriminated by MALDI-TOF are easily identified further by complementary simple tests, such as hemolysis done simultaneously or routinely during pathogen isolation.

Mechanistic insights into tRNA translocation on the ribosome

The ribosome is a highly conserved cellular complex and constitutes the center of protein biosynthesis. As the ribosome consists to about 2/3 of ribosomal RNA (rRNA), the rRNA is involved in most steps of translation. In order to investigate the role of some defined rRNA residues in different aspects of translation we use the atomic mutagenesis approach. This method allows the site-specific incorporation of unnatural nucleosides into the rRNA in the context of the complete 70S from Thermus aquaticus, and thereby exceeds the possibilities of conventional mutagenesis. We first studied ribosome-stimulated EF-G GTP hydrolysis. Here, we could show that the non-bridging phosphate oxygen of A2662, which is part of the Sarcin-Ricin-Loop, is required for EF-G GTPase activation by the ribosome. EF-G GTPase is a crucial step for tRNA translocation from the A- to the P-site, and from the P- to the E-site, respectively. We furthermore used the atomic mutagenesis approach to more precisely characterize the 23S rRNA functional groups involved in E-site tRNA binding. While the ribosomal A- and P-sites have been functionally well characterized in the past, the contribution of the E-site to protein biosynthesis is still poorly understood in molecular terms. Our data disclose the importance of the highly conserved E-site base pair G2421-C2395 for effective translation. Ribosomes with a disrupted G2421-C2395 base pair are defective in tRNA binding to the E-site. This results in an impaired translation of genuine mRNAs, while homo-polymeric templates are not affected. Cumulatively our data emphasize the importance of E-site tRNA occupancy and in particular the intactness of the 23S rRNA base pair G2421-C2395 for productive protein biosynthesis.

DNA probes for the detection of Fasciola hepatica in snails

Fasciola hepatica, also called the large liver fluke, is a trematode which can infect most mammals. Monitoring the infection rate of snails, which function as intermediate hosts and harbour larval stages of F. hepatica, is an important component of epidemiological studies on fascioliasis. For this purpose, DNA probes were generated which can be used for the detection of F. hepatica larvae in snails. Four highly repetitive DNA fragments were cloned in a plasmid vector and tested by Southern blot hybridization to the DNA of various trematodes for specificity and sensitivity. The probes Fhr-I, Fhr-II and Fhr-III hybridized only to F. hepatica DNA. Fhr-IV contained ribosomal RNA gene sequences and cross-hybridize with the DNA from various other trematode species. Squash blot analysis showed that the different probes were able to detect the parasite larvae in trematode-infected snails even as isolated single larvae. No signals were obtained in squash blots of uninfected snails. Probes Fhr-I, Fhr-II and Fhr-III are thus useful specific tools for studying the epidemiology of fascioliasis. The probe Fhr-IV, because of its broader spectrum, can be used to detect the larvae of a wide range of trematode species of waterbirds, which are the causative agents of swimmer's itch.

The sequence comparison index: A novel method for comparing proteins and proteomes

Historically morphological features were used as the primary means to classify organisms. However, the age of molecular genetics has allowed us to approach this field from the perspective of the organism's genetic code. Early work used highly conserved sequences, such as ribosomal RNA. The increasing number of complete genomes in the public data repositories provides the opportunity to look not only at a single gene, but at organisms' entire parts list. ^ Here the Sequence Comparison Index (SCI) and the Organism Comparison Index (OCI), algorithms and methods to compare proteins and proteomes, are presented. The complete proteomes of 104 sequenced organisms were compared. Over 280 million full Smith-Waterman alignments were performed on sequence pairs which had a reasonable expectation of being related. From these alignments a whole proteome phylogenetic tree was constructed. This method was also used to compare the small subunit (SSU) rRNA from each organism and a tree constructed from these results. The SSU rRNA tree by the SCI/OCI method looks very much like accepted SSU rRNA trees from sources such as the Ribosomal Database Project, thus validating the method. The SCI/OCI proteome tree showed a number of small but significant differences when compared to the SSU rRNA tree and proteome trees constructed by other methods. Horizontal gene transfer does not appear to affect the SCI/OCI trees until the transferred genes make up a large portion of the proteome. ^ As part of this work, the Database of Related Local Alignments (DaRLA) was created and contains over 81 million rows of sequence alignment information. DaRLA, while primarily used to build the whole proteome trees, can also be applied shared gene content analysis, gene order analysis, and creating individual protein trees. ^ Finally, the standard BLAST method for analyzing shared gene content was compared to the SCI method using 4 spirochetes. The SCI system performed flawlessly, finding all proteins from one organism against itself and finding all the ribosomal proteins between organisms. The BLAST system missed some proteins from its respective organism and failed to detect small ribosomal proteins between organisms. ^

Diversity of zooplankton assemblages at Station L4, Western English Channel, measured using next generation DNA sequencing

Background: Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel Next Generation Sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify the richness and diversity of a mixed zooplankton assemblage from a productive monitoring site in the Western English Channel. Methodology/Principle Findings: Plankton WP2 replicate net hauls (200 µm) were taken at the Western Channel Observatory long-term monitoring station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,042 sequences were obtained for all samples. The sequences clustered in to 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 138 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 75 taxonomic groups. Conclusions: The percentage of OTUs assigned to major eukaryotic taxonomic groups broadly aligns between the metagenetic and morphological analysis and are dominated by Copepoda. However, the metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for estimating diversity and species richness of zooplankton communities.

Effect of increased pCO2 on bacterial assemblage shifts in response to glucose addition in Fram Strait seawater mesocosms

Ocean acidification may stimulate primary production through increased availability of inorganic carbon in the photic zone, which may in turn change the biogenic flux of dissolved organic carbon (DOC) and the growth potential of heterotrophic bacteria. To investigate the effects of ocean acidification on marine bacterial assemblages, a two-by-three factorial mescosom experiment was conducted using surface sea water from the East Greenland Current in Fram Strait. Pyrosequencing of the V1-V2 region of bacterial 16S ribosomal RNA genes was used to investigate differences in the endpoint (Day 9) composition of bacterial assemblages in mineral nutrient-replete mesocosms amended with glucose (0 µm, 5.3 µm and 15.9 µm) under ambient (250 µatm) or acidified (400 µatm) partial pressures of CO2 (pCO2). All mesocosms showed low richness and diversity by Chao1 estimator and Shannon index, respectively, with general dominance by Gammaproteobacteria and Flavobacteria. Nonmetric multidimensional scaling analysis and two-way analysis of variance of the Jaccard dissimilarity matrix (97% similarity cut-off) demonstrated that the significant community shift between 0 µm and 15.9 µm glucose addition at 250 µatm pCO2 was eliminated at 400 µatm pCO2. These results suggest that the response potential of marine bacteria to DOC input may be altered under acidified conditions.

Environmental conditions and habitats of GOS sites and bottle data of HOT-ALOHA station

The Global Ocean Sampling (GOS) expedition is currently the largest and geographically most comprehensive metagenomic dataset, including samples from the Atlantic, Pacific, and Indian Oceans. This study makes use of the wide range of environmental conditions and habitats encompassed within the GOS sites in order to investigate the ecological structuring of bacterial and archaeal taxon ranks. Community structures based on taxonomically classified 16S ribosomal RNA (rRNA) gene fragments at phylum, class, order, family, and genus rank levels were examined using multivariate statistical analysis, and the results were inspected in the context of oceanographic environmental variables and structured habitat classifications. At all taxon rank levels, community structures of neritic, oceanic, estuarine biomes, as well as other exotic biomes (salt marsh, lake, mangrove), were readily distinguishable from each other. A strong structuring of the communities with chlorophyll a concentration and a weaker yet significant structuring with temperature and salinity were observed. Furthermore, there were significant correlations between community structures and habitat classification. These results were used for further investigation of one-to-one relationships between taxa and environment and provided indications for ecological preferences shaped by primary production for both cultured and uncultured bacterial and archaeal clades.

Metadata und statistic analysis of archael and bacterial sequences originating from sediments of the Håkon Mosby mud volcano (Z1-Z3, all habitats)

DNA extraction was carried out as described on the MICROBIS project pages (http://icomm.mbl.edu/microbis ) using a commercially available extraction kit. We amplified the hypervariable regions V4-V6 of archaeal and bacterial 16S rRNA genes using PCR and several sets of forward and reverse primers (http://vamps.mbl.edu/resources/primers.php). Massively parallel tag sequencing of the PCR products was carried out on a 454 Life Sciences GS FLX sequencer at Marine Biological Laboratory, Woods Hole, MA, following the same experimental conditions for all samples. Sequence reads were submitted to a rigorous quality control procedure based on mothur v30 (doi:10.1128/AEM.01541-09) including denoising of the flow grams using an algorithm based on PyroNoise (doi:10.1038/nmeth.1361), removal of PCR errors and a chimera check using uchime (doi:10.1093/bioinformatics/btr381). The reads were taxonomically assigned according to the SILVA taxonomy (SSURef v119, 07-2014; doi:10.1093/nar/gks1219) implemented in mothur and clustered at 98% ribosomal RNA gene V4-V6 sequence identity. V4-V6 amplicon sequence abundance tables were standardized to account for unequal sampling effort using 1000 (Archaea) and 2300 (Bacteria) randomly chosen sequences without replacement using mothur and then used to calculate inverse Simpson diversity indices and Chao1 richness (doi:10.2307/4615964). Bray-Curtis dissimilarities (doi:10.2307/1942268) between all samples were calculated and used for 2-dimensional non metric multidimensional scaling (NMDS) ordinations with 20 random starts (doi:10.1007/BF02289694). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. NMDS ordinations were compared and tested using Procrustes correlation analysis (doi:10.1007/BF02291478). All analyses were carried out with the R statistical environment and the packages vegan (available at: http://cran.r-project.org/package=vegan), labdsv (available at: http://cran.r-project.org/package=labdsv), as well as with custom R scripts. Operational taxonomic units at 98% sequence identity (OTU0.03) that occurred only once in the whole dataset were termed absolute single sequence OTUs (SSOabs; doi:10.1038/ismej.2011.132). OTU0.03 sequences that occurred only once in at least one sample, but may occur more often in other samples were termed relative single sequence OTUs (SSOrel). SSOrel are particularly interesting for community ecology, since they comprise rare organisms that might become abundant when conditions change.16S rRNA amplicons and metagenomic reads have been stored in the sequence read archive under SRA project accession number SRP042162.

Phylogenetic relationships between the Acantharea and the Polycystinea: A molecular perspective on Haeckel’s Radiolaria

Polycystine radiolaria are among few protistan groups that possess a comprehensive fossil record available for study by micropaleontologists. The Polycystinea and the Acantharea, whose skeletons do not become fossilized, were once members of the class “Radiolaria” (“Radiolaria” sensu lato: Polycystinea, Phaeodarea, and Acantharea) originally proposed by Haeckel but are now included in the superclass Actinopoda. Phylogenetic relationships within this superclass remain largely enigmatic. We investigated the evolutionary relationship of the Acantharea and the Polycystinea to other protists using phylogenetic analyses of 16S-like ribosomal RNA (rRNA) coding regions. We circumvented the need to culture these organisms by collecting and maintaining reproductive stages that contain many copies of their genomic DNA. This strategy facilitated extraction of genomic DNA and its purification from symbiont and prey DNA. Phylogenetic trees inferred from comparisons of 16S-like coding regions do not support a shared history between the Acantharea and the Polycystinea. However, the monophyly of the Acantharea and the separate monophyly of the Polycystinea (Spumellarida) are well supported by our molecular-based trees. The acantharian lineage branches among crown organisms whereas the polycystine lineage diverges before the radiation of the crown groups. We conclude that the Actinopoda does not represent a monophyletic evolutionary assemblage and recommend that this taxonomic designation be discarded.