Phylogeny and potential virulence of cryptic clade Escherichia coli species complex isolates derived from an arable field trial.

Analysis of Escherichia coli taxonomy has expanded into a species-complex with the identification of divergent cryptic clades. A key question is the evolutionary trajectory of these clades and their relationship to isolates of clinical or veterinary importance. Since they have some environmental association, we screened a collection of E. coli isolated from a long-term spring barley field trial for their presence. While most isolates clustered into the enteric-clade, four of them clustered into Clade-V, and one in Clade-IV. The Clade -V isolates shared >96% intra-clade average nucleotide sequence identity but <91% with other clades. Although pan-genomics analysis confirmed their taxonomy as Clade -V (E. marmotae), retrospective phylogroup PCR did not discriminate them correctly. Differences in metabolic and adherence gene alleles occurred in the Clade -V isolates compared to E. coli sensu scricto. They also encoded the bacteriophage phage-associated cyto-lethal distending toxin (CDT) and antimicrobial resistance (AMR) genes, including an ESBL, blaOXA-453. Thus, the isolate collection encompassed a genetic diversity, and included cryptic clade isolates that encode potential virulence factors. The analysis has determined the phylogenetic relationship of cryptic clade isolates with E. coli sensu scricto and indicates a potential for horizontal transfer of virulence factors.

Pectobacterium atrosepticum KDPG aldolase, Eda, participates in the Entner-Doudoroff pathway and independently inhibits expression of virulence determinants.

Pectobacterium carotovorum has an incomplete Entner-Doudoroff (ED) pathway, including enzyme 2-keto-3-deoxy-6-phosphogluconate aldolase (Eda) but lacking phosphogluconate dehydratase (Edd), while P. atrosepticum (Pba) has a complete pathway. To understand the role of the ED pathway in Pectobacterium infection, mutants of these two key enzymes, Δeda and Δedd, were constructed in Pba SCRI1039. Δeda exhibited significant decreased virulence on potato tubers and colonization in planta and was greatly attenuated in pectinase activity and the ability to use pectin breakdown products, including polygalacturonic acid (PGA) and galacturonic acid. These reduced phenotypes were restored following complementation with an external vector expressing eda. Quantitative reverse transcription PCR analysis revealed that expression of the pectinase genes pelA, pelC, pehN, pelW, and pmeB in Δeda cultured in pyruvate, with or without PGA, was significantly reduced compared to the wild type, while genes for virulence regulators (kdgR, hexR, hexA, and rsmA) remained unchanged. However, Δedd showed similar phenotypes to the wild type. To our knowledge, this is the first demonstration that disruption of eda has a feedback effect on inhibiting pectin degradation and that Eda is involved in building the arsenal of pectinases needed during infection by Pectobacterium.

The southwestern surgical congress multi-center trial on suspected common duct stones.

BACKGROUND: Choledocholithiasis is present in up to 15% of cholecystectomy patients. Treatment can be surgical, endoscopic, or via interventional radiology. We hypothesized significant heterogeneity between hospitals exists in the approach to suspected common duct stones. METHODS: A retrospective review of patients that had a preoperative MRCP, endoscopic ultrasound, endoscopic retrograde cholangiopancreatogram (ERCP), or intra-operative cholangiogram was performed. Comparisons were by Wilcoxon-Mann-Whitney tests with significance of p < 0.05 for paired variables and p < 0.017 for multiple comparisons. RESULTS: Twelve participating institutions identified 1263 patients (409 men and 854 women) with a median age of 49 years (IQR: 31-94). Liver function tests (LFT's) were elevated in 939 patients (75%), median bilirubin level 1.75 mg/dl (IQ: 0.8-3.7 mg/dl) and median common duct size 7 mm (IQR 5-10 mm). The most common initial procedure was cholecystectomy with IOC at seven institutions, endoscopy at four and MRCP at one. CONCLUSION: Significant variation exists within the surgical community regarding suspected common duct stones. These results underscore the need for a protocol for common duct stones to minimize multiple, redundant interventions.

Tangled history of a multigene family: The evolution of ISOPENTENYLTRANSFERASE genes.

ISOPENTENYLTRANSFERASE (IPT) genes play important roles in the initial steps of cytokinin synthesis, exist in plant and pathogenic bacteria, and form a multigene family in plants. Protein domain searches revealed that bacteria and plant IPT proteins were to assigned to different protein domains families in the Pfam database, namely Pfam IPT (IPTPfam) and Pfam IPPT (IPPTPfam) families, both are closely related in the P-loop NTPase clan. To understand the origin and evolution of the genes, a species matrix was assembled across the tree of life and intensively in plant lineages. The IPTPfam domain was only found in few bacteria lineages, whereas IPPTPfam is common except in Archaea and Mycoplasma bacteria. The bacterial IPPTPfam domain miaA genes were shown as ancestral of eukaryotic IPPTPfam domain genes. Plant IPTs diversified into class I, class II tRNA-IPTs, and Adenosine-phosphate IPTs; the class I tRNA-IPTs appeared to represent direct successors of miaA genes were found in all plant genomes, whereas class II tRNA-IPTs originated from eukaryotic genes, and were found in prasinophyte algae and in euphyllophytes. Adenosine-phosphate IPTs were only found in angiosperms. Gene duplications resulted in gene redundancies with ubiquitous expression or diversification in expression. In conclusion, it is shown that IPT genes have a complex history prior to the protein family split, and might have experienced losses or HGTs, and gene duplications that are to be likely correlated with the rise in morphological complexity involved in fine tuning cytokinin production.

Rapid and robust analytical protocol for E. coli STEC bacteria subspecies differentiation using whole cell MALDI mass spectrometry.

Whole cell MALDI is regularly used for the identification of bacteria to species level in clinical Microbiology laboratories. However, there remains a need to rapidly characterize and differentiate isolates below the species level to support outbreak management. We describe the implementation of a modified preparative approach for MALDI-MS combined with a custom analytical computational pipeline as a rapid procedure for subtyping Shigatoxigenic E. coli (STEC) and accurately identifying strain-specifying biomarkers. The technique was able to differentiate E. coli O157:H7 from other STEC. Within O157 serotype O157:H7 isolates were readily distinguishable from Sorbitol Fermenting O157 isolates. Overall, nine homogeneous groups of isolates were distinguished, each exhibiting distinct profiles of defining mass spectra features. This offers a robust analytical tool useable in reference/diagnostic public health scenarios.

MALDIrppa: quality control and robust analysis for mass spectrometry data.

Summary: This R package helps to implement a robust approach to deal with mass spectrometry (MS) data. It is aimed at alleviating reproducibility issues and pernicious effects of deviating signals on both data pre-processing and downstream data analysis. Based on robust statistical methods, it facilitates the identification and filtering of low-quality mass spectra and atypical peak profiles as well as monitoring and data handling through pre-processing, which extends existing computational tools for high-throughput data. Availability and implementation: MALDIrppa is implemented as a package for the R environment for data analysis and it is freely available to download from the CRAN repository at https://CRAN.R-project.org/package=MALDIrppa. Contact: javier.palarea@bioss.ac.uk.

Evolutionary Dynamics of the Cellulose Synthase Gene Superfamily in Grasses.

Phylogenetic analyses of cellulose synthase (CesA) and cellulose synthase-like (Csl) families from the cellulose synthase gene superfamily were used to reconstruct their evolutionary origins and selection histories. Counterintuitively, genes encoding primary cell wall CesAs have undergone extensive expansion and diversification following an ancestral duplication from a secondary cell wall-associated CesA. Selection pressure across entire CesA and Csl clades appears to be low, but this conceals considerable variation within individual clades. Genes in the CslF clade are of particular interest because some mediate the synthesis of (1,3;1,4)-ß-glucan, a polysaccharide characteristic of the evolutionarily successful grasses that is not widely distributed elsewhere in the plant kingdom. The phylogeny suggests that duplication of either CslF6 and/or CslF7 produced the ancestor of a highly conserved cluster of CslF genes that remain located in syntenic regions of all the grass genomes examined. A CslF6-specific insert encoding approximately 55 amino acid residues has subsequently been incorporated into the gene, or possibly lost from other CslFs, and the CslF7 clade has undergone a significant long-term shift in selection pressure. Homology modeling and molecular dynamics of the CslF6 protein were used to define the three-dimensional dispositions of individual amino acids that are subject to strong ongoing selection, together with the position of the conserved 55-amino acid insert that is known to influence the amounts and fine structures of (1,3;1,4)-ß-glucans synthesized. These wall polysaccharides are attracting renewed interest because of their central roles as sources of dietary fiber in human health and for the generation of renewable liquid biofuels.

The barley genome sequence assembly reveals three additional members of the CslF (1,3;1,4)-ß-glucan synthase gene family.

An important component of barley cell walls, particularly in the endosperm, is (1,3;1,4)-ß-glucan, a polymer that has proven health benefits in humans and that influences processability in the brewing industry. Genes of the cellulose synthase-like (Csl) F gene family have been shown to be involved in (1,3;1,4)-ß-glucan synthesis but many aspects of the biosynthesis are still unclear. Examination of the sequence assembly of the barley genome has revealed the presence of an additional three HvCslF genes (HvCslF11, HvCslF12 and HvCslF13) which may be involved in (1,3;1,4)-ß-glucan synthesis. Transcripts of HvCslF11 and HvCslF12 mRNA were found in roots and young leaves, respectively. Transient expression of these genes in Nicotiana benthamiana resulted in phenotypic changes in the infiltrated leaves, although no authentic (1,3;1,4)-ß-glucan was detected. Comparisons of the CslF gene families in cereals revealed evidence of intergenic recombination, gene duplications and translocation events. This significant divergence within the gene family might be related to multiple functions of (1,3;1,4)-ß-glucans in the Poaceae. Emerging genomic and global expression data for barley and other cereals is a powerful resource for characterising the evolution and dynamics of complete gene families. In the case of the CslF gene family, the results will contribute to a more thorough understanding of carbohydrate metabolism in grass cell walls.

Quantitative high resolution mapping of HvMLH3 foci in barley pachytene nuclei reveals a strong distal bias and weak interference.

In barley (Hordeum vulgare L.), chiasmata (the physical sites of genetic crossovers) are skewed towards the distal ends of chromosomes, effectively consigning a large proportion of genes to recombination coldspots. This has the effect of limiting potential genetic variability, and of reducing the efficiency of map-based cloning and breeding approaches for this crop. Shifting the sites of recombination to more proximal chromosome regions by forward and reverse genetic means may be profitable in terms of realizing the genetic potential of the species, but is predicated upon a better understanding of the mechanisms governing the sites of these events, and upon the ability to recognize real changes in recombination patterns. The barley MutL Homologue (HvMLH3), a marker for class I interfering crossovers, has been isolated and a specific antibody has been raised. Immunolocalization of HvMLH3 along with the synaptonemal complex transverse filament protein ZYP1, used in conjunction with fluorescence in situ hybridization (FISH) tagging of specific barley chromosomes, has enabled access to the physical recombination landscape of the barley cultivars Morex and Bowman. Consistent distal localization of HvMLH3 foci throughout the genome, and similar patterns of HvMLH3 foci within bivalents 2H and 3H have been observed. A difference in total numbers of HvMLH3 foci between these two cultivars has been quantified, which is interpreted as representing genotypic variation in class I crossover frequency. Discrepancies between the frequencies of HvMLH3 foci and crossover frequencies derived from linkage analysis point to the existence of at least two crossover pathways in barley. It is also shown that interference of HvMLH3 foci is relatively weak compared with other plant species.

Processing of Chlamydia abortus polymorphic membrane protein 18D during the chlamydial developmental cycle.

BACKGROUND: Chlamydia possess a unique family of autotransporter proteins known as the Polymorphic membrane proteins (Pmps). While the total number of pmp genes varies between Chlamydia species, all encode a single pmpD gene. In both Chlamydia trachomatis (C. trachomatis) and C. pneumoniae, the PmpD protein is proteolytically cleaved on the cell surface. The current study was carried out to determine the cleavage patterns of the PmpD protein in the animal pathogen C. abortus (termed Pmp18D). METHODOLOGY/PRINCIPAL FINDINGS: Using antibodies directed against different regions of Pmp18D, proteomic techniques revealed that the mature protein was cleaved on the cell surface, resulting in a100 kDa N-terminal product and a 60 kDa carboxy-terminal protein. The N-terminal protein was further processed into 84, 76 and 73 kDa products. Clustering analysis resolved PmpD proteins into three distinct clades with C. abortus Pmp18D, being most similar to those originating from C. psittaci, C. felis and C. caviae. CONCLUSIONS/SIGNIFICANCE: This study indicates that C. abortus Pmp18D is proteolytically processed at the cell surface similar to the proteins of C. trachomatis and C. pneumoniae. However, patterns of cleavage are species-specific, with low sequence conservation of PmpD across the genus. The absence of conserved domains indicates that the function of the PmpD molecule in chlamydia remains to be elucidated.

Nebulized hypertonic saline attenuates acute lung injury following trauma and hemorrhagic shock via inhibition of matrix metalloproteinase-13.

OBJECTIVE: We hypothesized that aerosolized inhaled hypertonic saline given at the onset of resuscitation will decrease acute lung injury following hemorrhagic shock, by inhibiting the release of epithelial derived proinflammatory mediators. DESIGN: Animal study. SETTING: Animal-care facility procedure room in a medical center. SUBJECTS: Adult male Sprague-Dawley rats. INTERVENTIONS: Rats underwent hemorrhagic shock followed by 2 hrs of resuscitation and 1 hr of observation. In the study group, nebulized hypertonic saline was delivered at the end of the shock period and after 1 hr and 2 hrs of resuscitation. MEASUREMENTS AND MAIN RESULTS: Shock provoked acute lung injury, which was attenuated with inhaled hypertonic saline (1.56 ± 0.2 mg protein/mL vs. 0.95 ± 0.3 mg protein/mL bronchoalveolar lavage fluid, shock vs. shock + hypertonic saline, p < .01). Nebulized hypertonic saline reduced inflammation (cytokine-induced neutrophil chemoattractant-1 accumulation in bronchoalveolar lavage fluid 5999 ± 1267 pg/mL vs. 3342 ± 859 pg/mL, shock vs. shock + hypertonic saline, p = .006). Additionally, nebulized hypertonic saline inhibited matrix -metalloproteinase-13 accumulation in the bronchoalveolar lavage fluid (1513 ± 337 pg/mL bronchoalveolar lavage fluid vs. 230 ± 19 pg/mL, shock vs. shock + hypertonic saline, p = .009) and pretreatment with a matrix metalloproteinase-13 inhibitor was sufficient to attenuate postshock acute lung injury (1.42 ± 0.09 mg/mL vs. 0.77 ± 0.23 mg/mL bronchoalveolar lavage protein, shock vs. shock + matrix metalloproteinase-13 inhibitor CL-82198, p = .002). CONCLUSION: Inhaled hypertonic saline attenuates postshock acute lung injury by exerting an anti-inflammatory effect on the pulmonary epithelium, suggesting a new clinical strategy to treat acute lung injury/acute respiratory distress syndrome.

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

BACKGROUND: The 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. RESULTS: As 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. CONCLUSIONS: By 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.

Genetic variability of Chlamydophila abortus strains assessed by PCR-RFLP analysis of polymorphic membrane protein-encoding genes.

This study used PCR-RFLP to investigate the genetic variability of pmp-encoding genes from fifty-two Chlamydophila abortus (C. abortus) strains originating from abortion cases from various geographical regions and host species. Six primer pairs were used to PCR-amplify DNA fragments encoding eighteen pmps. PCR products were digested using four restriction endonucleases and Bayesian methodologies were used to compare RFLP profiles and assign strains to a RFLP genotype. Strains could be assigned to 2 genotypes in the region encoding pmp18D, 3 genotypes in the regions encoding pmp1A-pmp2B, pmp3E-pmp6H and pmp11G-pmp15G, 4 genotypes in the region encoding pmp7G-pmp10G and 5 genotypes in the region encoding pmp16G-pmp17G. In all regions, the majority of strains (88.4-96.1%) had the same genotype as the reference strain S26/3. No correlation could be made between genotype, host species or geographical origin except for the two variant Greek strains, LLG and POS, which formed a discrete genotype in all pmp-encoding regions except pmp18D. Relative rates of evolution calculated for each pmp-encoding gene locus suggest that differing selective pressures and functional constraints may exist on C. abortus polymorphic membrane proteins. These findings suggest that although intraspecies heterogeneity of pmp-encoding genes in C. abortus is low, the sequence heterogeneity should be an important consideration when using pmps as the basis for novel diagnostics or vaccine development.

Comparative visualization of genetic and physical maps with Strudel.

UNLABELLED: Data visualization can play a key role in comparative genomics, for example, underpinning the investigation of conserved synteny patterns. Strudel is a desktop application that allows users to easily compare both genetic and physical maps interactively and efficiently. It can handle large datasets from several genomes simultaneously, and allows all-by-all comparisons between these. AVAILABILITY AND IMPLEMENTATION: Installers for Strudel are available for Windows, Linux, Solaris and Mac OS X at http://bioinf.scri.ac.uk/strudel/.

Trans-species polymorphism and selection in the MHC class II DRA genes of domestic sheep.

Highly polymorphic genes with central roles in lymphocyte mediated immune surveillance are grouped together in the major histocompatibility complex (MHC) in higher vertebrates. Generally, across vertebrate species the class II MHC DRA gene is highly conserved with only limited allelic variation. Here however, we provide evidence of trans-species polymorphism at the DRA locus in domestic sheep (Ovis aries). We describe variation at the Ovar-DRA locus that is far in excess of anything described in other vertebrate species. The divergent DRA allele (Ovar-DRA*0201) differs from the sheep reference sequences by 20 nucleotides, 12 of which appear non-synonymous. Furthermore, DRA*0201 is paired with an equally divergent DRB1 allele (Ovar-DRB1*0901), which is consistent with an independent evolutionary history for the DR sub-region within this MHC haplotype. No recombination was observed between the divergent DRA and B genes in a range of breeds and typical levels of MHC class II DR protein expression were detected at the surface of leukocyte populations obtained from animals homozygous for the DRA*0201, DRB1*0901 haplotype. Bayesian phylogenetic analysis groups Ovar-DRA*0201 with DRA sequences derived from species within the Oryx and Alcelaphus genera rather than clustering with other ovine and caprine DRA alleles. Tests for Darwinian selection identified 10 positively selected sites on the branch leading to Ovar-DRA*0201, three of which are predicted to be associated with the binding of peptide antigen. As the Ovis, Oryx and Alcelaphus genera have not shared a common ancestor for over 30 million years, the DRA*0201 and DRB1*0901 allelic pair is likely to be of ancient origin and present in the founding population from which all contemporary domestic sheep breeds are derived. The conservation of the integrity of this unusual DR allelic pair suggests some selective advantage which is likely to be associated with the presentation of pathogen antigen to T-cells and the induction of protective immunity.

Tablet--next generation sequence assembly visualization.

SUMMARY: Tablet is a lightweight, high-performance graphical viewer for next-generation sequence assemblies and alignments. Supporting a range of input assembly formats, Tablet provides high-quality visualizations showing data in packed or stacked views, allowing instant access and navigation to any region of interest, and whole contig overviews and data summaries. Tablet is both multi-core aware and memory efficient, allowing it to handle assemblies containing millions of reads, even on a 32-bit desktop machine. AVAILABILITY: Tablet is freely available for Microsoft Windows, Apple Mac OS X, Linux and Solaris. Fully bundled installers can be downloaded from http://bioinf.scri.ac.uk/tablet in 32- and 64-bit versions.

TOPALi v2: a rich graphical interface for evolutionary analyses of multiple alignments on HPC clusters and multi-core desktops.

UNLABELLED: 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.

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 alpha-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 alpha-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.

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.