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.

In situ synthesis of DNA microarray on functionalized cyclic olefin copolymer substrate.

Thermoplastic materials such as cyclic-olefin copolymers (COC) provide a versatile and cost-effective alternative to the traditional glass or silicon substrate for rapid prototyping and industrial scale fabrication of microdevices. To extend the utility of COC as an effective microarray substrate, we developed a new method that enabled for the first time in situ synthesis of DNA oligonucleotide microarrays on the COC substrate. To achieve high-quality DNA synthesis, a SiO(2) thin film array was prepatterned on the inert and hydrophobic COC surface using RF sputtering technique. The subsequent in situ DNA synthesis was confined to the surface of the prepatterned hydrophilic SiO(2) thin film features by precision delivery of the phosphoramidite chemistry using an inkjet DNA synthesizer. The in situ SiO(2)-COC DNA microarray demonstrated superior quality and stability in hybridization assays and thermal cycling reactions. Furthermore, we demonstrate that pools of high-quality mixed-oligos could be cleaved off the SiO(2)-COC microarrays and used directly for construction of DNA origami nanostructures. It is believed that this method will not only enable synthesis of high-quality and low-cost COC DNA microarrays but also provide a basis for further development of integrated microfluidics microarrays for a broad range of bioanalytical and biofabrication applications.

SplicerAV: a tool for mining microarray expression data for changes in RNA processing.

BACKGROUND: Over the past two decades more than fifty thousand unique clinical and biological samples have been assayed using the Affymetrix HG-U133 and HG-U95 GeneChip microarray platforms. This substantial repository has been used extensively to characterize changes in gene expression between biological samples, but has not been previously mined en masse for changes in mRNA processing. We explored the possibility of using HG-U133 microarray data to identify changes in alternative mRNA processing in several available archival datasets. RESULTS: Data from these and other gene expression microarrays can now be mined for changes in transcript isoform abundance using a program described here, SplicerAV. Using in vivo and in vitro breast cancer microarray datasets, SplicerAV was able to perform both gene and isoform specific expression profiling within the same microarray dataset. Our reanalysis of Affymetrix U133 plus 2.0 data generated by in vitro over-expression of HRAS, E2F3, beta-catenin (CTNNB1), SRC, and MYC identified several hundred oncogene-induced mRNA isoform changes, one of which recognized a previously unknown mechanism of EGFR family activation. Using clinical data, SplicerAV predicted 241 isoform changes between low and high grade breast tumors; with changes enriched among genes coding for guanyl-nucleotide exchange factors, metalloprotease inhibitors, and mRNA processing factors. Isoform changes in 15 genes were associated with aggressive cancer across the three breast cancer datasets. CONCLUSIONS: Using SplicerAV, we identified several hundred previously uncharacterized isoform changes induced by in vitro oncogene over-expression and revealed a previously unknown mechanism of EGFR activation in human mammary epithelial cells. We analyzed Affymetrix GeneChip data from over 400 human breast tumors in three independent studies, making this the largest clinical dataset analyzed for en masse changes in alternative mRNA processing. The capacity to detect RNA isoform changes in archival microarray data using SplicerAV allowed us to carry out the first analysis of isoform specific mRNA changes directly associated with cancer survival.

The genomic analysis of erythrocyte microRNA expression in sickle cell diseases.

BACKGROUND: Since mature erythrocytes are terminally differentiated cells without nuclei and organelles, it is commonly thought that they do not contain nucleic acids. In this study, we have re-examined this issue by analyzing the transcriptome of a purified population of human mature erythrocytes from individuals with normal hemoglobin (HbAA) and homozygous sickle cell disease (HbSS). METHODS AND FINDINGS: Using a combination of microarray analysis, real-time RT-PCR and Northern blots, we found that mature erythrocytes, while lacking ribosomal and large-sized RNAs, contain abundant and diverse microRNAs. MicroRNA expression of erythrocytes was different from that of reticulocytes and leukocytes, and contributed the majority of the microRNA expression in whole blood. When we used microRNA microarrays to analyze erythrocytes from HbAA and HbSS individuals, we noted a dramatic difference in their microRNA expression pattern. We found that miR-320 played an important role for the down-regulation of its target gene, CD71 during reticulocyte terminal differentiation. Further investigation revealed that poor expression of miR-320 in HbSS cells was associated with their defective downregulation CD71 during terminal differentiation. CONCLUSIONS: In summary, we have discovered significant microRNA expression in human mature erythrocytes, which is dramatically altered in HbSS erythrocytes and their defect in terminal differentiation. Thus, the global analysis of microRNA expression in circulating erythrocytes can provide mechanistic insights into the disease phenotypes of erythrocyte diseases.

Convergent differential regulation of parvalbumin in the brains of vocal learners.

Spoken language and learned song are complex communication behaviors found in only a few species, including humans and three groups of distantly related birds--songbirds, parrots, and hummingbirds. Despite their large phylogenetic distances, these vocal learners show convergent behaviors and associated brain pathways for vocal communication. However, it is not clear whether this behavioral and anatomical convergence is associated with molecular convergence. Here we used oligo microarrays to screen for genes differentially regulated in brain nuclei necessary for producing learned vocalizations relative to adjacent brain areas that control other behaviors in avian vocal learners versus vocal non-learners. A top candidate gene in our screen was a calcium-binding protein, parvalbumin (PV). In situ hybridization verification revealed that PV was expressed significantly higher throughout the song motor pathway, including brainstem vocal motor neurons relative to the surrounding brain regions of all distantly related avian vocal learners. This differential expression was specific to PV and vocal learners, as it was not found in avian vocal non-learners nor for control genes in learners and non-learners. Similar to the vocal learning birds, higher PV up-regulation was found in the brainstem tongue motor neurons used for speech production in humans relative to a non-human primate, macaques. These results suggest repeated convergent evolution of differential PV up-regulation in the brains of vocal learners separated by more than 65-300 million years from a common ancestor and that the specialized behaviors of learned song and speech may require extra calcium buffering and signaling.

Computational inference of neural information flow networks.

Determining how information flows along anatomical brain pathways is a fundamental requirement for understanding how animals perceive their environments, learn, and behave. Attempts to reveal such neural information flow have been made using linear computational methods, but neural interactions are known to be nonlinear. Here, we demonstrate that a dynamic Bayesian network (DBN) inference algorithm we originally developed to infer nonlinear transcriptional regulatory networks from gene expression data collected with microarrays is also successful at inferring nonlinear neural information flow networks from electrophysiology data collected with microelectrode arrays. The inferred networks we recover from the songbird auditory pathway are correctly restricted to a subset of known anatomical paths, are consistent with timing of the system, and reveal both the importance of reciprocal feedback in auditory processing and greater information flow to higher-order auditory areas when birds hear natural as opposed to synthetic sounds. A linear method applied to the same data incorrectly produces networks with information flow to non-neural tissue and over paths known not to exist. To our knowledge, this study represents the first biologically validated demonstration of an algorithm to successfully infer neural information flow networks.

Predictive Modeling of Loss-of-Heterozygosity Events in Yeast

During mitotic cell cycles, DNA experiences many types of endogenous and exogenous damaging agents that could potentially cause double strand breaks (DSB). In S. cerevisiae, DSBs are primarily repaired by mitotic recombination and as a result, could lead to loss-of-heterozygosity (LOH). Genetic recombination can happen in both meiosis and mitosis. While genome-wide distribution of meiotic recombination events has been intensively studied, mitotic recombination events have not been mapped unbiasedly throughout the genome until recently. Methods for selecting mitotic crossovers and mapping the positions of crossovers have recently been developed in our lab. Our current approach uses a diploid yeast strain that is heterozygous for about 55,000 SNPs, and employs SNP-Microarrays to map LOH events throughout the genome. These methods allow us to examine selected crossovers and unselected mitotic recombination events (crossover, noncrossover and BIR) at about 1 kb resolution across the genome. Using this method, we generated maps of spontaneous and UV-induced LOH events. In this study, we explore machine learning and variable selection techniques to build a predictive model for where the LOH events occur in the genome.

Randomly from the yeast genome, we simulated control tracts resembling the LOH tracts in terms of tract lengths and locations with respect to single-nucleotide-polymorphism positions. We then extracted roughly 1,100 features such as base compositions, histone modifications, presence of tandem repeats etc. and train classifiers to distinguish control tracts and LOH tracts. We found interesting features of good predictive values. We also found that with the current repertoire of features, the prediction is generally better for spontaneous LOH events than UV-induced LOH events.

A complex intronic enhancer regulates expression of the CFTR gene by direct interaction with the promoter.

Genes can maintain spatiotemporal expression patterns by long-range interactions between cis-acting elements. The cystic fibrosis transmembrane conductance regulator gene (CFTR) is expressed primarily in epithelial cells. An element located within a DNase I-hypersensitive site (DHS) 10 kb into the first intron was previously shown to augment CFTR promoter activity in a tissue-specific manner. Here, we reveal the mechanism by which this element influences CFTR transcription. We employed a high-resolution method of mapping DHS using tiled microarrays to accurately locate the intron 1 DHS. Transfection of promoter-reporter constructs demonstrated that the element displays classical tissue-specific enhancer properties and can independently recruit factors necessary for transcription initiation. In vitro DNase I footprinting analysis identified a protected region that corresponds to a conserved, predicted binding site for hepatocyte nuclear factor 1 (HNF1). We demonstrate by electromobility shift assays (EMSA) and chromatin immunoprecipitation (ChIP) that HNF1 binds to this element both in vitro and in vivo. Moreover, using chromosome conformation capture (3C) analysis, we show that this element interacts with the CFTR promoter in CFTR-expressing cells. These data provide the first insight into the three- dimensional (3D) structure of the CFTR locus and confirm the contribution of intronic cis-acting elements to the regulation of CFTR gene expression.

Validation of the detection of Pseudo-nitzschia spp. using specific RNA probes tested in a microarray format: Calibration of signal based on variability of RNA content with environmental conditions.

Harmful algal blooms (HAB) occur worldwide and cause health problems and economic damage to fisheries and tourism. Monitoring for toxic algae is therefore essential but is based primarily on light microscopy, which is time consuming and can be limited by insufficient morphological characters such that more time is needed to examine critical features with electron microscopy. Monitoring with molecular tools is done in only a few places world-wide. EU FP7 MIDTAL (Microarray Detection of Toxic Algae) used SSU and LSU rRNA genes as targets on microarrays to identify toxic species. In order to comply with current monitoring requirements to report cell numbers as the relevant threshold measurement to trigger closure of fisheries, it was necessary to calibrate our microarray to convert the hybridisation signal obtained to cell numbers. Calibration curves for two species of Pseudo-nitzschia for use with the MIDTAL microarray are presented to obtain cell numbers following hybridisation. It complements work presented by Barra et al. (2012b. Environ. Sci. Pollut. Res. doi: 10.1007/s11356-012-1330-1v) for two other Pseudo-nitzschia spp., Dittami and Edvardsen (2012a. J. Phycol. 48, 1050) for Pseudochatonella, Blanco et al. (2013. Harmful Algae 24, 80) for Heterosigma, McCoy et al. (2013. FEMS. doi: 10.1111/1574-6941.12277) for Prymnesium spp., Karlodinium veneficum, and cf. Chatonella spp. and Taylor et al. (2014. Harmful Algae, in press) for Alexandrium.

High-quality RNA extraction from copepods for Next Generation Sequencing: A comparative study

Despite the ecological importance of copepods, few Next Generation Sequencing studies (NGS) have been performed on small crustaceans, and a standard method for RNA extraction is lacking. In this study, we compared three commonly-used methods: TRIzol®, Aurum Total RNA Mini Kit and Qiagen RNeasy Micro Kit, in combination with preservation reagents TRIzol® or RNAlater®, to obtain high-quality and quantity of RNA from copepods for NGS. Total RNA was extracted from the copepods Calanus helgolandicus, Centropages typicus and Temora stylifera and its quantity and quality were evaluated using NanoDrop, agarose gel electrophoresis and Agilent Bioanalyzer. Our results demonstrate that preservation of copepods in RNAlater® and extraction with Qiagen RNeasy Micro Kit were the optimal isolation method for high-quality and quantity of RNA for NGS studies of C. helgolandicus. Intriguingly, C. helgolandicus 28S rRNA is formed by two subunits that separate after heat-denaturation and migrate along with 18S rRNA. This unique property of protostome RNA has never been reported in copepods. Overall, our comparative study on RNA extraction protocols will help increase gene expression studies on copepods using high-throughput applications, such as RNA-Seq and microarrays.

Identification of 5-Fluorouracil-inducible Target Genes Using cDNA Microarray Profiling

The fluoropyrimidine 5-Fluorouracil (5-FU) is widely used in the treatment of cancer. To identify novel downstream mediators of tumor cell response to 5-FU, we used DNA microarray technology to identify genes that are transcriptionally activated by 5-FU treatment in the MCF-7 breast cancer cell line. Of 2400 genes analyzed, 619 were up-regulated by >3-fold. Highly up-regulated genes (>6-fold) with signal intensities of >3000 were analyzed by Northern blot. Genes that were consistently found to be up-regulated were spermine/spermidine acetyl transferase (SSAT), annexin II, thymosin-beta-10, chaperonin-10, and MAT-8. Treatment of MCF-7 cells with the antifolate tomudex and DNA-damaging agent oxaliplatin also resulted in up-regulation of each of these targets. The 5-FU-induced activation of MAT-8, thymosin-beta-10, and chaperonin-10 was abrogated by inactivation of p53 in MCF-7 cells, whereas induction of SSAT and annexin II was significantly reduced in the absence of p53. Moreover, each of these genes contained more than one potential p53-binding site, suggesting that p53 may play an important regulatory role in 5-FU-induced expression of these genes. In addition, we found that basal expression levels of SSAT, annexin II, thymosin beta-10, and chaperonin-10 were increased (by approximately 2-3-fold), and MAT-8 expression dramatically increased (by approximately 10-fold) in a 5-FU-resistant colorectal cancer cell line (H630-R10) compared with the parental H630 cell line, suggesting these genes may be useful biomarkers of resistance. These results demonstrate the potential of DNA microarrays to identify novel genes involved in mediating the response of tumor cells to chemotherapy.

Erythropoietin receptor expression in non-small cell lung carcinoma: a question of antibody specificity

Immunohistochemical studies on formalin-fixed, paraffin-embedded (FFPE) tissue utilizing polyclonal antibodies form the cornerstone of many reports claiming to demonstrate erythropoietin receptor (EPOR) expression in malignant tissue. Recently, Elliott et al. (Blood 2006;107:1892-1895) reported that the antibodies commonly used to detect EPOR expression also detect non-EPOR proteins, and that their binding to EPOR was severely abrogated by two synthetic peptides based on the sequence of heat shock protein (HSP) 70, HSP70-2, and HSP70-5. We have investigated the specificity of the C20 antibody for detecting EPOR expression in non-small cell lung carcinoma (NSCLC) utilizing tissue microarrays. A total of 34 cases were available for study. Antibody absorbed with peptide resulted in marked suppression of cytoplasmic staining compared with nonabsorbed antibody. Four tumors that initially showed a membranous pattern of staining retained this pattern with absorbed antibody. Positive membranous immunoreactivity was also observed in 6 of 30 tumors that originally showed a predominantly cytoplasmic pattern of staining. Using the C20 antibody for Western blots, we detected three main bands, at 100, 66, and 59 kDa. Preincubation with either peptide caused abolition of the 66-kDa band, which contains non-EPOR sequences including heat shock peptides. These results call into question the significance of previous immunohistochemical studies of EPOR expression in malignancy and emphasize the need for more specific anti-EPOR antibodies to define the true extent of EPOR expression in neoplastic tissue

sscMap: An extensible Java application for connecting small-molecule drugs using gene-expression signatures

Background
Connectivity mapping is a process to recognize novel pharmacological and toxicological properties in small molecules by comparing their gene expression signatures with others in a database. A simple and robust method for connectivity mapping with increased specificity and sensitivity was recently developed, and its utility demonstrated using experimentally derived gene signatures.

Results
This paper introduces sscMap (statistically significant connections' map), a Java application designed to undertake connectivity mapping tasks using the recently published method. The software is bundled with a default collection of reference gene-expression profiles based on the publicly available dataset from the Broad Institute Connectivity Map 02, which includes data from over 7000 Affymetrix microarrays, for over 1000 small-molecule compounds, and 6100 treatment instances in 5 human cell lines. In addition, the application allows users to add their custom collections of reference profiles and is applicable to a wide range of other 'omics technologies.

Conclusion
The utility of sscMap is two fold. First, it serves to make statistically significant connections between a user-supplied gene signature and the 6100 core reference profiles based on the Broad Institute expanded dataset. Second, it allows users to apply the same improved method to custom-built reference profiles which can be added to the database for future referencing. The software can be freely downloaded from http://purl.oclc.org/NET/sscMap

Identification of gene networks associated with acute myeloid leukemia by comparative molecular methylation and expression profiling

Around 80% of acute myeloid leukemia (AML) patients achieve a complete remission, however many will relapse and ultimately die of their disease. The association between karyotype and prognosis has been studied extensively and identified patient cohorts as having favourable [e.g. t(8; 21), inv (16)/t(16; 16), t(15; 17)], intermediate [e.g. cytogenetically normal (NK-AML)] or adverse risk [e.g. complex karyotypes]. Previous studies have shown that gene expression profiling signatures can classify the sub-types of AML, although few reports have shown a similar feature by using methylation markers. The global methylation patterns in 19 diagnostic AML samples were investigated using the Methylated CpG Island Amplification Microarray (MCAM) method and CpG island microarrays containing 12,000 CpG sites. The first analysis, comparing favourable and intermediate cytogenetic risk groups, revealed significantly differentially methylated CpG sites (594 CpG islands) between the two subgroups. Mutations in the NPM1 gene occur at a high frequency (40%) within the NK-AML subgroup and are associated with a more favourable prognosis in these patients. A second analysis comparing the NPM1 mutant and wild-type research study subjects again identified distinct methylation profiles between these two subgroups. Network and pathway analysis revealed possible molecular mechanisms associated with the different risk and/or mutation sub-groups. This may result in a better classification of the risk groups, improved monitoring targets, or the identification of novel molecular therapies.

The Delta Np63 Proteins Are Key Allies of BRCA1 in the Prevention of Basal-Like Breast Cancer

Little is known about the origin of basal-like breast cancers, an aggressive disease that is highly similar to BRCA1-mutant breast cancers. p63 family proteins that are structurally related to the p53 suppressor protein are known to function in stem cell regulation and stratified epithelia development in multiple tissues, and p63 expression may be a marker of basal-like breast cancers. Here we report that Delta Np63 isoforms of p63 are transcriptional targets for positive regulation by BRCA1. Our analyses of breast cancer tissue microarrays and BRCA1-modulated breast cancer cell lines do not support earlier reports that p63 is a marker of basal-like or BRCA1 mutant cancers. Nevertheless, we found that BRCA1 interacts with the specific p63 isoform Delta Np63 gamma along with transcription factor isoforms AP-2 alpha and AP-2 gamma. BRCA1 required Delta Np63 gamma and AP-2 gamma to localize to an intronic enhancer region within the p63 gene to upregulate transcription of the Delta Np63 isoforms. In mammary stem/progenitor cells, siRNA- mediated knockdown of Delta Np63 expression resulted in genomic instability, increased cell proliferation, loss of DNA damage checkpoint control, and impaired growth control. Together, our findings establish that transcriptional upregulation of Delta Np63 proteins is critical for BRCA1 suppressor function and that defects in BRCA1-Delta Np63 signaling are key events in the pathogenesis of basal-like breast cancer. Cancer Res; 71( 5); 1933-44. (c) 2011 AACR.