Cystic fibrosis gene repair: correction of ΔF508 using ZFN and CRISPR/Cas9 guide RNA gene editing tools

Cystic Fibrosis (CF) is an autosomal recessive monogenic disorder caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene with the ΔF508 mutation accounting for approximately 70% of all CF cases worldwide. This thesis investigates whether existing zinc finger nucleases designed in this lab and CRISPR/gRNAs designed in this thesis can mediate efficient homology-directed repair (HDR) with appropriate donor repair plasmids to correct CF-causing mutations in a CF cell line. Firstly, the most common mutation, ΔF508, was corrected using a pair of existing ZFNs, which cleave in intron 9, and the donor repair plasmid pITR-donor-XC, which contains the correct CTT sequence and two unique restriction sites. HDR was initially determined to be <1% but further analysis by next generation sequencing (NGS) revealed HDR occurred at a level of 2%. This relatively low level of repair was determined to be a consequence of distance from the cut site to the mutation and so rather than designing a new pair of ZFNs, the position of the existing intron 9 ZFNs was exploited and attempts made to correct >80% of CF-causing mutations. The ZFN cut site was used as the site for HDR of a mini-gene construct comprising exons 10-24 from CFTR cDNA (with appropriate splice acceptor and poly A sites) to allow production of full length corrected CFTR mRNA. Finally, the ability to cleave closer to the mutation and mediate repair of CFTR using the latest gene editing tool CRISPR/Cas9 was explored. Two CRISPR gRNAs were tested; CRISPR ex10 was shown to cleave at an efficiency of 15% and CRISPR in9 cleaved at 3%. Both CRISPR gRNAs mediated HDR with appropriate donor plasmids at a rate of ~1% as determined by NGS. This is the first evidence of CRISPR induced HDR in CF cell lines.

Dairy-Met: compositional metagenomic analysis of milk and cheeses

Unpasteurised milk and many cheeses contain a diverse microbiological population. These microorganisms play important roles in dairy foods and can, for example, contribute to the development of flavours and aromas, determine safety, cause spoilage or enhance the health of the consumer. It is thus important to understand thoroughly the microorganisms present in these food types. Traditional culture dependent and culture-independent methods have provided much detail regarding the microbial content of dairy foods. However, the development of next-generation DNA sequencing technologies has revolutionised our knowledge of complex microbial environments. Throughout this thesis we observe the benefits of applying these technologies to provide a detailed understanding of the bacterial content of dairy foods, including those present in milk pre- and post-pasteurisation, Irish farmhouse cheeses and commercially produced cheeses which encounter a discolouration defect, as well as to study genomic changes in microbes associated with dairy foods. Through the application of these state-of-the-art technologies we identified the presence of microorganisms not previously associated with dairy foods.

Identification of RNA editing in the human exome and development of DARNED DatabaseSF

RNA editing is a biological phenomena that alters nascent RNA transcripts by insertion, deletion and/or substitution of one or a few nucleotides. It is ubiquitous in all kingdoms of life and in viruses. The predominant editing event in organisms with a developed central nervous system is Adenosine to Inosine deamination. Inosine is recognized as Guanosine by the translational machinery and reverse-transcriptase. In primates, RNA editing occurs frequently in transcripts from repetitive regions of the genome. In humans, more than 500,000 editing instances have been identified, by applying computational pipelines on available ESTs and high-throughput sequencing data, and by using chemical methods. However, the functions of only a small number of cases have been studied thoroughly. RNA editing instances have been found to have roles in peptide variants synthesis by non-synonymous codon substitutions, transcript variants by alterations in splicing sites and gene silencing by miRNAs sequence modifications. We established the Database of RNA EDiting (DARNED) to accommo-date the reference genomic coordinates of substitution editing in human, mouse and fly transcripts from published literatures, with additional information on edited genomic coordinates collected from various databases e.g. UCSC, NCBI. DARNED contains mostly Adenosine to Inosine editing and allows searches based on genomic region, gene ID, and user provided sequence. The Database is accessible at http://darned.ucc.ie RNA editing instances in coding region are likely to result in recoding in protein synthesis. This encouraged me to focus my research on the occurrences of RNA editing specific CDS and non-Alu exonic regions. By applying various filters on discrepancies between available ESTs and their corresponding reference genomic sequences, putative RNA editing candidates were identified. High-throughput sequencing was used to validate these candidates. All predicted coordinates appeared to be either SNPs or unedited.

The impact of a variety of factors on the obesity associated gut microbiota

The obesity pandemic has become perhaps the most prevalent health issue of our time, with more than 10% of the world’s population now being obese. Obesity can be defined as abnormal or excess fat accumulation that may impair health and results from an imbalance between energy intake and energy expenditure. A decrease in physical activity due to an increase in sedentary forms of work, changing modes of transport and increasing urbanization is likely a major contributory factor. Diet is another major factor with the increased availability and intake of calorie dense, high fat foods being of global concern. Notably, with respect to this thesis, over the last decade advances in the field of next generation sequencing (NGS) have facilitated investigations to determine the relationship between the gut microbiota and obesity. This thesis examines the impact of a variety of factors on the obesity associated gut microbiota. Overall the results presented in this thesis highlight that microbial diversity is influenced by diet, exercise, antibiotics and disease state, however it is only through further understanding of the structure and function that we can identify targets that can impact on health.

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.

Food preservation using antifungal lactic acid bacteria

Fungal spoilage of food and feed prevails as a major problem for the food industry. The use antifungal-producing lactic acid bacteria (LAB) may represent a safer, natural alternative to the use of chemical preservatives in foods. A large scale screen was undertaken to identify a variety of LAB with antifungal properties from plant, animal and human sources. A total of 6,720 LAB colonies were isolated and screened for antifungal activity against the indicator Penicillium expansum. 94 broad-spectrum producers were identified through 16S rRNA sequencing with the majority of the population comprising Lactobacillus plantarum isolates. Six broad-spectrum isolates were consequently characterised. Pedicococcus pentosaceous 54 displayed potent anti-mould capabilities in pear, plum and grape models and may represent an ideal candidate for use in the beverage industry. Two antifungal Lb. plantarum isolates were assessed for their technological robustness and potential as biopreservatives in refrigerated foods. Lb. plantarum 16 and 62 displayed high levels of tolerance to freeze-drying, low temperature exposure and high salt concentrations. Both lactobacilli were introduced as supplements into orange juice to retard the growth of the spoilage yeast Rhodotorula mucilaginosa. Furthermore the isolates were applied as adjuncts in yoghurt production to successfully reduce yeast growth. Lb. plantarum 16 proved to be the optimal inhibitor of yeast growth in both food matrices. To date there is limited information available describing the mechanisms behind fungal inhibition by LAB. The effects of concentrated cell-free supernatant (cCFS), derived from Lb. plantarum 16, on the growth of two food-associated moulds was assessed microscopically. cCFS completely inhibited spore, germ tube and hyphal development. A transcriptomic approach was undertaken to determine the impact of antifungal activity on Aspergillus fumigatus Af293. A variety of genes, most notably those involved in cellular metabolism, were found to have their transcription modulated in response to cCFS which is indicative of global cellular shutdown. This study provides the first insights into the molecular targets of antifungal compounds produced by LAB. The genome sequence of the steep water isolate Lb. plantarum 16 was determined. The complete genome of Lb. plantarum16 consists of a single circular chromosome of 3,044,738 base pairs with an average G+C content of 44.74 % in addition to eight plasmids. The genome represents the smallest of this species to date while harbouring the largest plasmid complement. Some features of particular interest include the presence of two prophages, an interrupted plantaricin cluster and a chromosomal and plasmid encoded polysaccharide cluster. The sequence presented here provides a suitable platform for future studies elucidating the mechanisms governing antifungal production.

Aspects of the biology of Mya arenaria and Ensis spp. (Mollusca; Bivalvia) in the Irish Sea and adjacent areas

This study was undertaken to investigate the general biology, including the reproductive cycle and health status, of two clam taxa in Irish waters, with particular reference to the Irish Sea area. Monthly samples of the soft shell clam, Mya arenaria, were collected from Bannow Bay, Co. Wexford, Ireland, for sixteen months, and of the razor clam, Ensis spp. from the Skerries region (Irish Sea) between June 2010 and September 2011. In 2010, M. arenaria in Bannow Bay matured over the summer months, with both sexes either ripe or spawning by August. The gonads of both sexes of E. siliqua developed over autumn and winter 2010, with the first spawning individuals being recorded in January 2011. Two unusually cold winters, followed by a warmer than average spring, appear to have affected M. arenaria and E. siliqua gametogenesis at these sites. It was noted that wet weight of E. siliqua dropped significantly in the summer of both 2010 and 2011, after spawning, which may impact on the economic viability of fishing during this period. Additional samples of M. arenaria were collected at Flaxfort (Ireland), and Ensis spp. at Oxwich (Wales), and the pathology of all clams was examined using both histological and molecular methods. No pathogenic conditions were observed in M. arenaria while Prokaryote inclusions, trematode parasites, Nematopsis spp. and inflammatory pathologies were observed at low incidences in razor clams from Ireland but not from Wales; the first time these conditions have been reported in Ensis spp. in northern European waters. Mya arenaria from sites in Europe and eastern and western North America were investigated for genetic variation using both mitochondrial (cytochrome oxidase I (COI) and 16S ribosomal RNA genes) and nuclear markers (10 microsatellite loci). Both mitochondrial CO1 and all nuclear markers showed reduced levels of variation in certain European samples, with significant differences in haplotype and allelic composition between most samples, particularly those from the two different continents, but with the same common haplotypes or alleles throughout the range. The appearance of certain unique rare haplotypes and microsatellite alleles in the European samples suggest a complicated origin involving North American colonization but also possible southern European Pleistocene refugia. Specimens of Ensis spp. were obtained from five coastal areas around Ireland and Wales and species-specific PCR primers were used to amplify the internal transcribed spacer region 1 (ITS1) and the mitochondrial DNA CO1 gene and all but 15 razor clams were identified as Ensis siliqua. Future investigations should focus on continued monitoring of reproductive biology and pathology of the two clam taxa (in particular, to assess the influence of environmental change), and on genetics of southern European M. arenaria and sequencing the CO1 gene in Ensis individuals to clarify species identity

Proteomic approach to oxidative stress in Daphnia magna

The keystone aquatic organism Daphnia magna is extensively used to assess the toxicity of chemicals. This has recently lead to an increase in the omics literature focusing on daphnids, an increase fuelled by the sequencing of the Daphnia pulex genome. Yet, no omics study has looked directly at oxidative stress (OS) in daphnids, even though OS is of primary importance in the response of aquatic organisms to their changing environment and is often induced by anthropogenic xenobiotics. This thesis thus focuses on the application of redox-proteomics, the study of the oxidative modification of proteins, to D. magna Specifically, daphnids were exposed to copper or paraquat, two well studied prooxidants, and protein carbonyls were labelled with fluorescein-5-thiosemicarbazide prior to twodimensional electrophoresis (2DE). This showed clearly that both compounds affect a different portion of the proteome. The identified proteins indicated that energy metabolism was affected by paraquat, while copper induced a reduction of the heat shock response (heat shock proteins, proteases and chaperones) a counterintuitive result which may be adaptative to metal toxicity in arthropods. The same approach was then applied to the study of the toxicity mechanism of silver nanoparticles (AgNP), an increasingly utilised form of silver with expected environmental toxicity, and its comparison to silver nitrate. The results demonstrate that, although less toxic than silver ions, AgNP toxicity functions through a different mechanism. AgNP toxicity is thus not a product of silver dissolution and increased protein carbonylation indicates that AgNP cause OS. Interestingly three of the four tested compounds altered vitellogenin levels and oxidation. Vitellogenins could thus represent an interesting subproteome for the detection of stress in daphnids. Finally, an experiment with oxidised BSA demonstrates the applicability of solid phase hydrazide in the enrichment of undigested carbonylated proteins.

Helicobacter pylori: comparative genomics and structure-function analysis of the flagellum biogenesis protein HP0958

Helicobacter pylori is a gastric pathogen which infects ~50% of the global population and can lead to the development of gastritis, gastric and duodenal ulcers and carcinoma. Genome sequencing of H. pylori revealed high levels of genetic variability; this pathogen is known for its adaptability due to mechanisms including phase variation, recombination and horizontal gene transfer. Motility is essential for efficient colonisation by H. pylori. The flagellum is a complex nanomachine which has been studied in detail in E. coli and Salmonella. In H. pylori, key differences have been identified in the regulation of flagellum biogenesis, warranting further investigation. In this study, the genomes of two H. pylori strains (CCUG 17874 and P79) were sequenced and published as draft genome sequences. Comparative studies identified the potential role of restriction modification systems and the comB locus in transformation efficiency differences between these strains. Core genome analysis of 43 H. pylori strains including 17874 and P79 defined a more refined core genome for the species than previously published. Comparative analysis of the genome sequences of strains isolated from individuals suffering from H. pylori related diseases resulted in the identification of “disease-specific” genes. Structure-function analysis of the essential motility protein HP0958 was performed to elucidate its role during flagellum assembly in H. pylori. The previously reported HP0958-FliH interaction could not be substantiated in this study and appears to be a false positive. Site-directed mutagenesis confirmed that the coiled-coil domain of HP0958 is involved in the interaction with RpoN (74-284), while the Zn-finger domain is required for direct interaction with the full length flaA mRNA transcript. Complementation of a non-motile hp0958-null derivative strain of P79 with site-directed mutant alleles of hp0958 resulted in cells producing flagellar-type extrusions from non-polar positions. Thus, HP0958 may have a novel function in spatial localisation of flagella in H. pylori

The role of bacterial interspecies communication in polymicrobial infection of the cystic fibrosis lung

There is an increasing appreciation of the polymicrobial nature of bacterial infections associated with Cystic Fibrosis (CF) and of the important role for interactions in influencing bacterial virulence and response to therapy. Patients with CF are co-infected with Pseudomonas aeruginosa, Burkholderia cenocepacia and Stenotrophomonas maltophilia. These latter bacteria produce signal molecules of the diffusible signal factor (DSF) family, which are cis-2-unsaturated fatty acids. Previous studies showed that DSF from S. maltophilia leads to altered biofilm formation and increased tolerance to antibiotics in P. aeruginosa and that these responses require the P. aeruginosa sensor kinase PA1396. The work in this thesis aims of further elucidate the influence and mechanism of DSF signalling on P. aeruginosa and examine the role that such interspecies signalling play in infection of the CF airway. Next generation sequencing technologies targeting the 16S ribosomal RNA gene were applied to DNA and RNA isolated from sputum taken from cohorts of CF and non-CF subjects to characterise the bacterial community. In parallel, metabolomics analysis of sputum provided insight into the environment of the CF airway. This analysis revealed a number of observations including; that differences in metabolites occur in sputum taken from clinically stable CF patients and those with exacerbation and DNA- and RNA-based methods suggested that a strong relationship existed between the abundance of specific strict anaerobes and fluctuations in the level of metabolites during exacerbation. DSF family signals were also detected in the sputum and a correlation with the presence of DSFproducing organisms was observed. To examine the signal transduction mechanisms used by P. aeruginosa, bioinformatics with site directed mutagenesis were employed to identify signalling partners for PA1396. A pathway suggesting a role for a number of proteins in the regulation of several factors following DSF recognition by PA1396 were observed.

Identification and analysis of mechanisms involved in a broad-spectrum disease resistant mutant of the winter wheat cv. Guardian

M66 an X-ray induced mutant of winter wheat (Triticum aestivum) cv. Guardian exhibits broad-spectrum resistance to powdery mildew (Blumeria graminis f. sp. tritici), yellow rust (Puccinia striiformis f. sp. tritici), and leaf rust (Puccinia recondita f. sp. tritici), along with partial resistance to stagnonospora nodorum blotch (caused by the necrotroph Stagonosporum nodorum) and septoria tritici blotch (caused by the hemibiotroph Mycosphaerella graminicola) compared to the parent plant ‘Guardian’. Analysis revealed that M66 exhibited no symptoms of infection following artificial inoculation with Bgt in the glasshouse after adult growth stage (GS 45). Resistance in M66 was associated with widespread leaf flecking which developed during tillering. Flecking also occurred in M66 leaves without Bgt challenge; as a result grain yields were reduced by approximately 17% compared to ‘Guardian’ in the absence of disease. At the seedling stage, M66 exhibited partial resistance. M66, along with Tht mutants (Tht 12, Tht13), also exhibit increased tolerance to environmental stresses (abiotic), such as drought and heat stress at seedling and adult growth stages, However, adult M66 exhibited increased susceptibility to the aphid Schizaphis graminum compared to ‘Guardian’. Resistance to Bgt in M66 was characterized with increased and earlier H2O2 accumulation at the site of infection which resulted in increased papilla formation in epidermal cells, compared to ‘Guardian’. Papilla formation was associated with reduced pathogen ingress and haustorium formation, indicating that the primary cause of resistance in M66 was prevention of pathogen penetration. Heat treatment at 46º C prior to challenge with Bgt also induced partial disease resistance to Blumeria graminis f. sp. tritici in ‘Guardian’ and M66 seedlings. This was characterized by a delay in primary infection, due to increased production of ROS species, such as hydrogen peroxide, ROS-scavenging enzymes and Hsp70, resulting in cross-linking of cell wall components prior to inoculation. This actively prevented the fungus from penetrating the epidermal cell wall. Proteomics analysis using 2-D gel electrophoresis identified primary and secondary disease resistance effects in M66 including detection of ROS scavenging enzymes (4, 24 hai), such as ascorbate peroxidase and a superoxidase dismutase isoform (CuZnSOD) in M66 which were absent from ‘Guardian’. Chitinase (PR protein) was also upregulated (24 hai) in M66 compared to ‘Guardian’.Monosomic and ditelosomic analysis of M66 revealed that the mutation in M66 is located on the long arm of chromosome 2B (2BL). Chromosome 2BL is known to have key genes involved in resistance to pathogens such as those causing stripe rust and powdery mildew. The TaMloB1 gene, an orthologue of the barley Mlo gene, is also located on chromosome 2BL. Sanger sequencing of part of the coding sequence revealed no deletions in the TaMloB1 gene between ‘Guardian’ and M66.

The relevance of bacterial isoprenoid biosynthetic pathways for host-microbe interactions

This thesis was undertaken to investigate the relevance of two bacterial isoprenoid biosynthetic pathways (Mevalonate (MVAL) and 2-C-methyl-D-erythritol 4-phosphate (MEP)) for host-microbe interactions. We determined a significant reduction in microbial diversity in the murine gut microbiota (by next generation sequencing) following oral administration of a common anti-cholesterol drug Rosuvastatin (RSV) that targets mammalian and bacterial HMG-CoA reductase (HMG-R) for inhibition of MVAL formation. In tandem we identified significant hepatic and intestinal off-target alterations to the murine metabolome indicating alterations in inflammation, bile acid profiles and antimicrobial peptide synthesis with implications on community structure of the gastrointestinal microbiota in statin-treated animals. However we found no effect on local Short Chain Fatty Acid biosynthesis (metabolic health marker in our model). We demonstrated direct inhibition of bacterial growth in-vitro by RSV which correlated with reductions in bacterial MVAL formation. However this was only at high doses of RSV. Our observations demonstrate a significant RSV-associated impact on the gut microbiota prompting similar human analysis. Successful deletion of another MVAL pathway enzyme (HMG-CoA synthase (mvaS)) involved in Listeria monocytogenes EGDe isoprenoid biosynthesis determined that the enzyme is non-essential for normal growth and in-vivo pathogenesis of this pathogen. We highlight potential evidence for alternative means of synthesis of the HMG-CoA substrate that could render mvaS activity redundant under our test conditions. Finally, we showed by global gene expression analysis (Massive Analysis of cDNA Ends (MACE RNA-seq) a significant role for the penultimate MEP pathway metabolite (E)-4-hydroxy-3-methyl-2-but-2-enyl pyrophosphate (HMBPP) in significant up regulation of genes of immunity and antigen presentation in THP-1 cells at nanomolar levels. We infected THP-1 cells with wild type or HMBPP under/over-producing L. monoctyogenes EGDe mutants and determined subtle effects of HMBPP upon overall host responses to Listeria infection. Overall our findings provide greater insights regarding bacterial isoprenoid biosynthetic pathways for host-microbe/microbe-host dialogue.

The biodiversity, control and genetic characterisation of the lactococcal 936 group phages

Phages belonging to the 936 group represent one of the most prevalent and frequently isolated phages in dairy fermentation processes using Lactococcus lactis as the primary starter culture. In recent years extensive research has been carried out to characterise this phage group at a genomic level in an effort to understand how the 936 group phages dominate this particular niche and cause regular problems during large scale milk fermentations. This thesis describes a large scale screening of industrial whey samples, leading to the isolation of forty three genetically different lactococcal phages. Using multiplex PCR, all phages were identified as members of the 936 group. The complete genome of thirty eight of these phages was determined using next generation sequencing technologies which identified several regions of divergence. These included the structural region surrounding the major tail protein, the replication region as well as the genes involved in phage DNA packing. For a number of phages the latter genomic region was found to harbour genes encoding putative orphan methyltransferases. Using small molecule real time (SMRT) sequencing and heterologous gene expression, the target motifs for several of these MTases were determined and subsequently shown to actively protect phage DNA from restriction endonuclease activity. Comparative analysis of the thirty eight phages with fifty two previously sequenced members of this group showed that the core genome consists of 28 genes, while the non-core genome was found to fluctuate irrespective of geographical location or time of isolation. This study highlights the continued need to perform large scale characterisation of the bacteriophage populations infecting industrial fermentation facilities in effort to further our understanding dairy phages and ways to control their proliferation.