Combined analysis of variation in core, accessory and regulatory genome regions provides a super-resolution view into the evolution of bacterial populations

The use of whole-genome phylogenetic analysis has revolutionized our understanding of the evolution and spread of many important bacterial pathogens due to the high resolution view it provides. However, the majority of such analyses do not consider the potential role of accessory genes when inferring evolutionary trajectories. Moreover, the recently discovered importance of the switching of gene regulatory elements suggests that an exhaustive analysis, combining information from core and accessory genes with regulatory elements could provide unparalleled detail of the evolution of a bacterial population. Here we demonstrate this principle by applying it to a worldwide multi-host sample of the important pathogenic E. coli lineage ST131. Our approach reveals the existence of multiple circulating subtypes of the major drug–resistant clade of ST131 and provides the first ever population level evidence of core genome substitutions in gene regulatory regions associated with the acquisition and maintenance of different accessory genome elements.

Biodiscovery of natural products from microbes associated with Irish coastal sponges

Marine sponges have been an abundant source of new metabolites in recent years. The symbiotic association between the bacteria and the sponge has enabled scientists to access the bacterial diversity present within the bacterial/sponge ecosystem. This study has focussed on accessing the bacterial diversity in two Irish coastal marine sponges, namely Amphilectus fucorum and Eurypon major. A novel species from the genus Aquimarina has been isolated from the sponge Amphilectus fucorum. The study has also resulted in the identification of an α–Proteobacteria, Pseudovibrio sp. as a potential producer of antibiotics. Thus a targeted based approach to specifically cultivate Pseudovibrio sp. may prove useful for the development of new metabolites from this particular genus. Bacterial isolates from the marine sponge Haliclona simulans were screened for anti–fungal activity and one isolate namely Streptomyces sp. SM8 displayed activity against all five fungal strains tested. The strain was also tested for anti–bacterial activity and it showed activity against both against B. subtilis and P. aeruginosa. Hence a combinatorial approach involving both biochemical and genomic approaches were employed in an attempt to identify the bioactive compounds with these activities which were being produced by this strain. Culture broths from Streptomyces sp. SM8 were extracted and purified by various techniques such as reverse–phase HPLC, MPLC and ash chromatography. Anti–bacterial activity was observed in a fraction which contained a hydroxylated saturated fatty acid and also another compound with a m/z 227 but further structural elucidation of these compounds proved unsuccessful. The anti–fungal fractions from SM8 were shown to contain antimycin–like compounds, with some of these compounds having different retention times from that of an antimycin standard. A high–throughput assay was developed to screen for novel calcineurin inhibitors using yeast as a model system and three putative bacterial extracts were found to be positive using this screen. One of these extracts from SM8 was subsequently analysed using NMR and the calcineurin inhibition activity was con rmed to belong to a butenolide type compound. A H. simulans metagenomic library was also screened using the novel calcineurin inhibitor high–throughput assay system and eight clones displaying putative calcineurin inhibitory activity were detected. The clone which displayed the best inhibitory activity was subsequently sequenced and following the use of other genetic based approaches it became clear that the inhibition was being caused by a hypothetical protein with similarity to a hypothetical Na+/Ca2+ exchanger protein. The Streptomyces sp. SM8 genome was sequenced from a fragment library using Roche 454 pyrosequencing technology to identify potential secondary metabolism clusters. The draft genome was annotated by IMG/ER using the Prodigal pipeline. The Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession AMPN00000000. The genome contains genes which appear to encode for several polyketide synthases (PKS), non–ribosomal peptide synthetases (NRPS), terpene and siderophore biosynthesis and ribosomal peptides. Transcriptional analyses led to the identification of three hybrid clusters of which one is predicted to be involved in the synthesis of antimycin, while the functions of the others are as yet unknown. Two NRPS clusters were also identified, of which one may be involved in gramicidin biosynthesis and the function of the other is unknown. A Streptomyces sp. SM8 NRPS antC gene knockout was constructed and extracts from the strain were shown to possess a mild anti–fungal activity when compared to the SM8 wild–type. Subsequent LCMS analysis of antC mutant extracts confirmed the absence of the antimycin in the extract proving that the observed anti–fungal activity may involve metabolite(s) other than antimycin. Anti–bacterial activity in the antC gene knockout strain against P. aeruginosa was reduced when compared to the SM8 wild–type indicating that antimycin may be contributing to the observed anti–bacterial activity in addition to the metabolite(s) already identified during the chemical analyses. This is the first report of antimycins exhibiting anti–bacterial activity against P. aeruginosa. One of the hybrid clusters potentially involved in secondary metabolism in SM8 that displayed high and consistent levels of gene–expression in RNA studies was analysed in an attempt to identify the metabolite being produced by the pathway. A number of unusual features were observed following bioinformatics analysis of the gene sequence of the cluster, including a formylation domain within the NRPS cluster which may add a formyl group to the growing chain. Another unusual feature is the lack of AT domains on two of the PKS modules. Other unusual features observed in this cluster is the lack of a KR domain in module 3 of the cluster and an aminotransferase domain in module 4 for which no clear role has been hypothesised.

Assessment of the immunomodulatory effects of raw/farm milk and probiotic bacteria

The overall aims of this study were to investigate the differences between raw/farm milk and pasteurised milk with respect to potential immune modifying effects following consumption and investigate the bacterial composition of raw milk compared to pasteurised milk. Furthermore, in this thesis, panels of potential probiotic bacteria from the Bifidobacterium and Lactobacillus genera were investigated. The overall bacterial composition of raw milk was compared with pasteurised milk using samples obtained from commercial milk producers around Ireland using next generation sequencing technology (454 pyrosequencing). Here the presence of previously unrecognised and diverse bacterial populations in unpasteurised cow’s milk was identified. Futhermore the bacterial content of pasteurised milk was found to be more diverse than previously thought. The global response of the adenocarcinoma cell line HT-29 to raw milk and pasteurised milk exposures were also characterised using whole genome microarray technology. Over one thousand differentially expressed genes were identified which were found to be involved in a plethora of cellular functions. Interestingly a reduction in immune related activity (e.g. Major histocompatability complex class II signalling and T and B cell proliferation) was identified in cells exposed to pasteurised milk compared with raw milk exposures. Further studies comparing human cell response to raw versus pasteurised milk was performed using peripheral blood mononuclear cells (PBMC) from healthy donors. A reduction in CD14 was identified following raw milk exposures compared with pasteurised milk and the pattern of cytokine production may indicate that gram positive bacteria in the raw milk were contributing to the differences in the cellular response to raw versus pasteurised milk. Panels of potentially probiotic bacteria (comprising of lactobacilli and bifidobacteria) were further assessed for immunomodulatory capabilities using cell culture based models. Gene expression and cytokine production were used to evaluate stimulated and unstimulated (LPS) cellular responses as well as interaction mechanisms

Functional genomics of commensal Lactobacilli

Catabolic flexibility affords a bacterium the ability to utilise different sugar sources as carbon for energy. This is important for commensal lactobacilli like Lactobacillus ruminis which can be exposed to a variety of carbohydrates in vivo. However, little is known about the fermentation capabilities, metabolic pathways, genetic diversity or potential survival mechanisms used by L. ruminis in vivo. A combination of in vitro and in silico techniques was used to identify the catabolic pathways of L. ruminis. I also compared 16 L. ruminis strains using a panel of biochemical and survival assays, genetically, whole genome sequencing and RNA sequencing. Multi locus sequence typing revealed that strains clustered according to their host sources. Transcriptome analysis by RNAseq of two motile strains under three growth conditions, including swarming, identified the up-regulation of carbohydrate-related genes under swarming conditions. This suggests that carbohydrate flexibility may have an uncharacterised role in L. ruminis swarming. Following on from the assessment of L. ruminis catabolic flexibility, the porcine diet was supplemented with galactooligosaccharides or L. ruminis ATCC 25644 plus galactooligosaccharides. Supplementation of the porcine diet with galactooligosaccharide had no effect on microbiota diversity. In contrast, the L. ruminis plus galactooligosaccharide treatment significantly reduced the microbiota diversity. Diet is a major factor that affects the diversity of the gut microbiota. In order to get a more thorough understanding of diet and gut health in animals such as racehorses and domesticated herbivores, I determined the core microbiota of animals consuming different feeds. Interestingly, the gut microbiota diversity correlated with the host phylogeny of the animal. The genome of Lactobacillus equi (2.19 Mb), isolated from a healthy Irish thoroughbred was also sequenced and annotated, and comprised 2,263 predicted genes. The large repertoire of predicted carbohydrate-related genes may offer L. equi an advantage in the complex and harsh hindgut environment. In summary, this thesis uses functional genomics to assess the effect that carbohydrates have on commensal lactobacilli and the microbiota as a whole.

Characterisation of bacteriophage-host interactions in Lactococcus lactis

Lactococcus lactis is used extensively world-wide for the production of fermented dairy products. Bacteriophages (phages) infecting L. lactis can result in slow or incomplete fermentations, or may even cause total fermentation failure. Therefore, bacteriophages disrupting L. lactis fermentation are of economic concern. This thesis employed a multifaceted approach to investigate various molecular aspects of phage-host interaction in L. lactis. The genome sequence of an Irish dairy starter strain, the prophage-cured L. lactis subsp. cremoris UC509.9, was studied. The 2,250,427 bp circular chromosome represents the smallest among its sequenced lactococcal equivalents. The genome displays clear genetic adaptation to the dairy niche in the form of extensive reductive evolution. Gene prediction identified 2066 protein-encoding genes, including 104 which showed significant homology to transposase-specifying genes. Over 9 % of the identified genes appear to be inactivated through stop codons or frame shift mutations. Many pseudogenes were found in genes that are assigned to carbohydrate and amino acid transport and metabolism orthologous groups, reflecting L. lactis UC509.9’s adaptation to the lactose and casein-rich dairy environment. Sequence analysis of the eight plasmids of L. lactis revealed extensive adaptation to the dairy environment. Key industrial phenotypes were mapped and novel lactococcal plasmid-associated genes highlighted. In addition to chromosomally-encoded bacteriophage resistance systems, six functional such systems were identified, including two abortive infection systems, AbiB and AbiD1, explaining the observed phage resistance of L. lactis UC509.9 Molecular analysis suggests that the constitutive expression of AbiB is not lethal to cells, suggesting the protein is expressed in an un/inactivated form. Analysis of 936 species phage sk1-escape mutants of AbiB revealed that all such mutants harbour mutations in orf6, which encodes the major capsid protein. Results suggest that the major capsid protein is required for activation of the AbiB system, although this requires furrther investigations. Temporal transcriptomes of L. lactis UC509.9 undergoing lytic infection with either one of two distinct bacteriophages, Tuc2009 and c2, was determined and compared to the transcriptome of uninfected UC509.9 cells. Whole genome microarrays performed at various time-points post-infection demonstrated a rather modest impact on host transcription. Alterations in the UC509.9 transcriptome during lytic infection appear phage-specific, with a relatively small number of differentially transcribed genes shared between infection with either Tuc2009 or c2. Transcriptional profiles of both bacteriophages during lytic infection was shown to generally correlate with previous studies and allowed the confirmation of previously predicted promoter sequences. Bioinformatic analysis of genomic regions encoding the presumed cell wall polysaccharide (CW PS) biosynthesis gene cluster of several strains of L. lactis was performed. Results demonstrate the presence of three dominant genetic types of this gene cluster, termed type A, B and C. These regions were used for the development of a multiplex PCR to identify CW PS genotype of various lactococcal strains. Analysis of 936 species phage receptor binding protein phylogeny (RBP) and CW PS genotype revealed an apparent correlation between RBP phylogeny and CW PS type, thereby providing a partial explanation for the observed narrow host range of 936 phages. Further analysis of the genetic locus encompassing the presumed CW PS biosynthesis operon of eight strains identified as belonging to the CW PS C (geno)type, revealed the presence of a variable region among the examined strains. The obtained comparative analysis allowed for the identification of five subgroups of the C type, named C1 to C5. We purified an acidic polysaccharide from the cell wall of L. lactis 3107 (C2 subtype) and confirmed that it is structurally different from the CW PS of the C1 subtype L. lactis MG1363. Combinations of genes from the variable region of C2 subtype were amplified from L. lactis 3107 and introduced into a mutant of the C1 subtype L. lactis NZ9000 (a direct derivative of MG1363) deficient in CW PS biosynthesis. The resulting recombinant mutant synthesized a CW PS with a composition characteristic for that of the C2 subtype L. lactis 3107 and not the wildtype C1 L. lactis NZ9000. The recombinant mutant exhibited a changed phage resistance/sensitivity profile consistent with that of L. lactis 3107, which unambiguously demonstrated that L. lactis 3107 CW PS is the host cell surface receptor of two bacteriophages belonging to the P335 species as well as phages that are member of the 936 species. The research presented in this thesis has significantly advanced our understanding of L. lactis bacteriophage-host interactions in several ways. Firstly, the examination of plasmidencoded bacteriophage resistance systems has allowed inferences to be made regarding the mode of action of AbiB, thereby providing a platform for further elucidation of the molecular trigger of this system. Secondly, the phage infection transcriptome data presented, in addition to previous work, has made L. lactis a model organism in terms of transcriptomic studies of bacteriophage-host interactions. And finally, the research described in this thesis has for the first time explicitly revealed the nature of a carbohydrate bacteriophage receptor in L. lactis, while also providing a logical explanation for the observed narrow host ranges exhibited by 936 and P335 phages. Future research in discerning the structures of other L. lactis CW PS, combined with the determination of the molecular interplay between receptor binding proteins of these phages and CW PS will allow an in depth understanding of the mechanism by which the most prevalent lactococcal phages identify and adsorb to their specific host.

Molecular analysis of evolutionary relationships of Phaseolus dumosus with the gene pool of common bean

Valuable genetic variation for bean breeding programs is held within the common bean secondary gene pool which consists of Phaseolus albescens, P. coccineus, P. costaricensis, and P. dumosus. However, the use of close relatives for bean improvement is limited due to the lack of knowledge about genetic variation and genetic plasticity of many of these species. Characterisation and analysis of the genetic diversity is necessary among beans' wild relatives; in addition, conflicting phylogenies and relationships need to be understood and a hypothesis of a hybrid origin of P. dumosus needs to be tested. This thesis research was orientated to generate information about the patterns of relationships among the common bean secondary gene pool, with particular focus on the species Phaseolus dumosus. This species displays a set of characteristics of agronomic interest, not only for the direct improvement of common bean but also as a source of valuable genes for adaptation to climate change. Here I undertake the first comprehensive study of the genetic diversity of P. dumosus as ascertained from both nuclear and chloroplast genome markers. A germplasm collection of the ancestral forms of P. dumosus together with wild, landrace and cultivar representatives of all other species of the common bean secondary gene pool, were used to analyse genetic diversity, phylogenetic relationships and structure of P. dumosus. Data on molecular variation was generated from sequences of cpDNA loci accD-psaI spacer, trnT-trnL spacer, trnL intron and rps14-psaB spacer and from the nrDNA the ITS region. A whole genome DArT array was developed and used for the genotyping of P. dumosus and its closes relatives. 4208 polymorphic markers were generated in the DArT array and from those, 742 markers presented a call rate >95% and zero discordance. DArT markers revealed a moderate genetic polymorphism among P. dumosus samples (13% of polymorphic loci), while P. coccineus presented the highest level of polymorphism (88% of polymorphic loci). At the cpDNA one ancestral haplotype was detected among all samples of all species in the secondary genepool. The ITS region of P. dumosus revealed high homogeneity and polymorphism bias to P. coccineus genome. Phylogenetic reconstructions made with Maximum likelihood and Bayesian methods confirmed previously reported discrepancies among the nuclear and chloroplast genomes of P. dumosus. The outline of relationships by hybridization networks displayed a considerable number of interactions within and between species. This research provides compelling evidence that P. dumosus arose from hybridisation between P. vulgaris and P. coccineus and confirms that P. costaricensis has likely been involved in the genesis or backcrossing events (or both) in the history of P. dumosus. The classification of the specie P. persistentus was analysed based on cpDNA and ITS sequences, the results found this species to be highly related to P. vulgaris but not too similar to P. leptostachyus as previously proposed. This research demonstrates that wild types of the secondary genepool carry a significant genetic variation which makes this a valuable genetic resource for common bean improvement. The DArT array generated in this research is a valuable resource for breeding programs since it has the potential to be used in several approaches including genotyping, discovery of novel traits, mapping and marker-trait associations. Efforts should be made to search for potential populations of P. persistentus and to increase the collection of new populations of P. dumosus, P. albescens and P. costaricensis that may provide valuable traits for introgression into common bean and other Phaseolus crops.

The application of next generation sequencing to profile microbe related cheese quality defects

High throughput next generation sequencing, together with advanced molecular methods, has considerably enhanced the field of food microbiology. By overcoming biases associated with culture dependant approaches, it has become possible to achieve novel insights into the nature of food-borne microbial communities. In this thesis, several different sequencing-based approaches were applied with a view to better understanding microbe associated quality defects in cheese. Initially, a literature review provides an overview of microbe-associated cheese quality defects as well as molecular methods for profiling complex microbial communities. Following this, 16S rRNA sequencing revealed temporal and spatial differences in microbial composition due to the time during the production day that specific commercial cheeses were manufactured. A novel Ion PGM sequencing approach, focusing on decarboxylase genes rather than 16S rRNA genes, was then successfully employed to profile the biogenic amine producing cohort of a series of artisanal cheeses. Investigations into the phenomenon of cheese pinking formed the basis of a joint 16S rRNA and whole genome shotgun sequencing approach, leading to the identification of Thermus species and, more specifically, the pathway involved in production of lycopene, a red coloured carotenoid. Finally, using a more traditional approach, the effect of addition of a facultatively heterofermentative Lactobacillus (Lactobacillus casei) to a Swiss-type cheese, in which starter activity was compromised, was investigated from the perspective of its ability to promote gas defects and irregular eye formation. X-ray computed tomography was used to visualise, using a non-destructive method, the consequences of the undesirable gas formation that resulted. Ultimately this thesis has demonstrated that the application of molecular techniques, such as next generation sequencing, can provide a detailed insight into defect-causing microbial populations present and thereby may underpin approaches to optimise the quality and consistency of a wide variety of cheeses.

Ruthenibacterium lactatiformans gen. nov., sp.nov., an anaerobic, lactate-producing member of the family Ruminococcaceae isolated from human faeces

Two novel strains of Gram-stain-negative, rod-shaped, obligately anaerobic, non-spore-forming, non-motile bacteria were isolated from the faeces of healthy human subjects. The strains, designated as 585-1T and 668, were characterized by mesophilic fermentative metabolism, production of d-lactic acid, succinic acid and acetic acid as end products of d-glucose fermentation, prevalence of C18 : 1 ω9, C18 : 1 ω9 aldehyde, C16 : 0 and C16 : 1 ω7c fatty acids, presence of glycine, glutamic acid, lysine, alanine and aspartic acid in the petidoglycan peptide moiety and lack of respiratory quinones. Whole genome sequencing revealed the DNA G+C content was 56.4–56.6 mol%. The complete 16S rRNA gene sequences of the two strains shared 91.7/91.6 % similarity with Anaerofilum pentosovorans FaeT, 91.3/91.2 % with Gemmiger formicilis ATCC 27749T and 88.9/88.8 % with Faecalibacterium prausnitzii ATCC 27768T. On the basis of chemotaxonomic and genomic properties it was concluded that the strains represent a novel species in a new genus within the family Ruminococcaceae , for which the name Ruthenibacterium lactatiformans gen. nov., sp. nov. is proposed. The type strain of Ruthenibacterium lactatiformans is 585-1T (=DSM 100348T=VKM B-2901T).