Computational Analysis of HLA-presentation of Non-synonymous Recipient Mismatches Indicates Effect on the Risk of Chronic Graft-vs.-Host Disease After Allogeneic HSCT.

Genetic mismatches in protein coding genes between allogeneic hematopoietic stem cell transplantation (allo-HSCT) recipient and donor can elicit an alloimmunity response via peptides presented by the recipient HLA receptors as minor histocompatibility antigens (mHAs). While the impact of individual mHAs on allo-HSCT outcome such as graft-vs.-host and graft-vs.-leukemia effects has been demonstrated, it is likely that established mHAs constitute only a small fraction of all immunogenic non-synonymous variants. In the present study, we have analyzed the genetic mismatching in 157 exome-sequenced sibling allo-HSCT pairs to evaluate the significance of polymorphic HLA class I associated peptides on clinical outcome. We applied computational mismatch estimation approaches based on experimentally verified HLA ligands available in public repositories, published mHAs, and predicted HLA-peptide affinites, and analyzed their associations with chronic graft-vs.-host disease (cGvHD) grades. We found that higher estimated recipient mismatching consistently increased the risk of severe cGvHD, suggesting that HLA-presented mismatching influences the likelihood of long-term complications in the patient. Furthermore, computational approaches focusing on estimation of HLA-presentation instead of all non-synonymous mismatches indiscriminately may be beneficial for analysis sensitivity and could help identify novel mHAs.

Blockade of HMGB1 Attenuates Diabetic Nephropathy in Mice.

Activation of TLR2 or TLR4 by endogenous ligands such as high mobility group box 1 (HMGB1) may mediate inflammation causing diabetic kidney injury. We determined whether blockade of HMGB1 signaling by: (1) supra-physiological production of endogenous secretory Receptor for Advanced Glycation End-products (esRAGE), a receptor for HMGB1; (2) administration of HMGB1 A Box, a specific competitive antagonist, would inhibit development of streptozotocin induced diabetic nephropathy (DN). Wild-type diabetic mice developed albuminuria, glomerular injuries, interstitial fibrosis and renal inflammation. Using an adeno-associated virus vector, systemic over-expression of esRAGE afforded significant protection from all parameters. No protection was achieved by a control vector which expressed human serum albumin. Administration of A Box was similarly protective against development of DN. To determine the mechanism(s) of protection, we found that whilst deficiency of TLR2, TLR4 or RAGE afforded partial protection from development of DN, over-expression of esRAGE provided additional protection in TLR2-/-, modest protection against podocyte damage only in TLR4-/- and no protection in RAGE-/- diabetic mice, suggesting the protection provided by esRAGE was primarily through interruption of RAGE and TLR4 pathways. We conclude that strategies to block the interaction between HMGB1 and its receptors may be effective in preventing the development of DN.

Hidden genomic MHC disparity between HLA-matched sibling pairs in hematopoietic stem cell transplantation.

Matching classical HLA alleles between donor and recipient is an important factor in avoiding adverse immunological effects in HSCT. Siblings with no differences in HLA alleles, either due to identical-by-state or identical-by-descent status, are considered to be optimal donors. We carried out a retrospective genomic sequence and SNP analysis of 336 fully HLA-A, -B, -DRB1 matched and 14 partially HLA-matched sibling HSCT pairs to determine the level of undetected mismatching within the MHC segment as well as to map their recombination sites. The genomic sequence of 34 genes locating in the MHC region revealed allelic mismatching at 1 to 8 additional genes in partially HLA-matched pairs. Also, fully matched pairs were found to have mismatching either at HLA-DPB1 or at non-HLA region within the MHC segment. Altogether, 3.9% of fully HLA-matched HSCT pairs had large genomic mismatching in the MHC segment. Recombination sites mapped to certain restricted locations. The number of mismatched nucleotides correlated with the risk of GvHD supporting the central role of full HLA matching in HSCT. High-density genome analysis revealed that fully HLA-matched siblings may not have identical MHC segments and even single allelic mismatching at any classical HLA gene often implies larger genomic differences along MHC.

Accuracy of Programs for the Determination of Human Leukocyte Antigen Alleles from Next-Generation Sequencing Data.

The human leukocyte antigen (HLA) genes code for proteins that play a central role in the function of the immune system by presenting peptide antigens to T cells. As HLA genes show extremely high genetic polymorphism, HLA typing at the allele level is demanding and is based on DNA sequencing. Determination of HLA alleles is warranted as HLA alleles are major genetic risk factors in autoimmune diseases and are matched in transplantation. Here, we compared the accuracy of several published HLA-typing algorithms that are based on next-generation sequencing (NGS) data. As genome sequencing is becoming increasingly common in research, we wanted to test how well HLA alleles can be deduced from genome data produced in studies with objectives other than HLA typing and in platforms not especially designed for HLA typing. The accuracies were assessed using datasets consisting of NGS data produced using an in-house sequencing platform, including the full 4 Mbp HLA segment, from 94 stem cell transplantation patients and exome sequences from 63 samples of the 1000 Genomes collection. In the patient dataset, none of the software gave perfect results for all the samples and genes when programs were used with the default settings. However, we found that ensemble prediction of the results or modifications of the settings could be used to improve accuracy. For the exome-only data, most of the algorithms did not perform very well. The results indicate that the use of these algorithms for accurate HLA allele determination is not straightforward when based on NGS data not especially targeted to the HLA typing and their accurate use requires HLA expertise.

Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells.

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

Toll-Like Receptor 4 Deficiency Improves Short-term Renal Function but not Long-term Graft Survival in a Fully MHC-Mismatched Murine Model of Renal Allograft Transplantation.

BACKGROUND: We have previously demonstrated that absence of myeloid differentiation primary response gene 88 (MyD88) induced donor-specific kidney allograft tolerance. The upstream pathways of MyD88 that mediate this process, however, remain unclear. Toll-like receptor 4 (TLR4) is an innate immune receptor that is dependent upon MyD88 for activity of its dominant signaling pathway. Here, we investigated the role of TLR4 in kidney allograft rejection using a fully major histocompatibility complex-mismatched, life-sustaining, murine model of renal allograft rejection. METHODS: Donor (BALB/c) and recipient (C57BL/6) mice either both deficient or sufficient for TLR4 underwent heterotopic renal allograft transplantation, with an additional group of mice receiving renal isografts as controls. Survival was assessed up to 100 days posttransplantation. Animals were also sacrificed 14 days posttransplantation for assessment of the acute allograft rejection response. RESULTS: Both wild-type (WT) and TLR4 allografts showed inferior survival compared to isografts, with no difference in survival between the allograft groups. Serum creatinine was lower in TLR4 allografts at day 14 posttransplantation compared with WT allografts, but this was not sustained by day 100. At day 14 posttransplant, increased CD11c dendritic cell accumulation, expression of IL-2 and indoleamine 2,3-dioxygenase were evident in TLR4 compared with WT allografts, whereas expression of inducible nitric oxide synthase was decreased. CONCLUSIONS: Acute kidney allograft rejection was modestly attenuated in TLR4 mice; however, long-term allograft survival and function were not affected in our model. Protection against acute rejection may involve increased accumulation of CD11c cells and indoleamine 2,3-dioxygenase expression.

Population whole-genome bisulfite sequencing across two tissues highlights the environment as the principal source of human methylome variation.

BACKGROUND: CpG methylation variation is involved in human trait formation and disease susceptibility. Analyses within populations have been biased towards CpG-dense regions through the application of targeted arrays. We generate whole-genome bisulfite sequencing data for approximately 30 adipose and blood samples from monozygotic and dizygotic twins for the characterization of non-genetic and genetic effects at single-site resolution. RESULTS: Purely invariable CpGs display a bimodal distribution with enrichment of unmethylated CpGs and depletion of fully methylated CpGs in promoter and enhancer regions. Population-variable CpGs account for approximately 15-20 % of total CpGs per tissue, are enriched in enhancer-associated regions and depleted in promoters, and single nucleotide polymorphisms at CpGs are a frequent confounder of extreme methylation variation. Differential methylation is primarily non-genetic in origin, with non-shared environment accounting for most of the variance. These non-genetic effects are mainly tissue-specific. Tobacco smoking is associated with differential methylation in blood with no evidence of this exposure impacting cell counts. Opposite to non-genetic effects, genetic effects of CpG methylation are shared across tissues and thus limit inter-tissue epigenetic drift. CpH methylation is rare, and shows similar characteristics of variation patterns as CpGs. CONCLUSIONS: Our study highlights the utility of low pass whole-genome bisulfite sequencing in identifying methylome variation beyond promoter regions, and suggests that targeting the population dynamic methylome of tissues requires assessment of understudied intergenic CpGs distal to gene promoters to reveal the full extent of inter-individual variation.

Non-random aneuploidy specifies subgroups of pilocytic astrocytoma and correlates with older age.

Pilocytic astrocytoma (PA) is the most common brain tumor in children but is rare in adults, and hence poorly studied in this age group. We investigated 222 PA and report increased aneuploidy in older patients. Aneuploid genomes were identified in 45% of adult compared with 17% of pediatric PA. Gains were non-random, favoring chromosomes 5, 7, 6 and 11 in order of frequency, and preferentially affecting non-cerebellar PA and tumors with BRAF V600E mutations and not with KIAA1549-BRAF fusions or FGFR1 mutations. Aneuploid PA differentially expressed genes involved in CNS development, the unfolded protein response, and regulators of genomic stability and the cell cycle (MDM2, PLK2),whose correlated programs were overexpressed specifically in aneuploid PA compared to other glial tumors. Thus, convergence of pathways affecting the cell cycle and genomic stability may favor aneuploidy in PA, possibly representing an additional molecular driver in older patients with this brain tumor.

Interrogation of allelic chromatin states in human cells by high-density ChIP-genotyping.

Allele-specific (AS) assessment of chromatin has the potential to elucidate specific cis-regulatory mechanisms, which are predicted to underlie the majority of the known genetic associations to complex disease. However, development of chromatin landscapes at allelic resolution has been challenging since sites of variable signal strength require substantial read depths not commonly applied in sequencing based approaches. In this study, we addressed this by performing parallel analyses of input DNA and chromatin immunoprecipitates (ChIP) on high-density Illumina genotyping arrays. Allele-specificity for the histone modifications H3K4me1, H3K4me3, H3K27ac, H3K27me3, and H3K36me3 was assessed using ChIP samples generated from 14 lymphoblast and 6 fibroblast cell lines. AS-ChIP SNPs were combined into domains and validated using high-confidence ChIP-seq sites. We observed characteristic patterns of allelic-imbalance for each histone-modification around allele-specifically expressed transcripts. Notably, we found H3K4me1 to be significantly anti-correlated with allelic expression (AE) at transcription start sites, indicating H3K4me1 allelic imbalance as a marker of AE. We also found that allelic chromatin domains exhibit population and cell-type specificity as well as heritability within trios. Finally, we observed that a subset of allelic chromatin domains is regulated by DNase I-sensitive quantitative trait loci and that these domains are significantly enriched for genome-wide association studies hits, with autoimmune disease associated SNPs specifically enriched in lymphoblasts. This study provides the first genome-wide maps of allelic-imbalance for five histone marks. Our results provide new insights into the role of chromatin in cis-regulation and highlight the need for high-depth sequencing in ChIP-seq studies along with the need to improve allele-specificity of ChIP-enrichment.

A new LC-MS/MS method for the clinical determination of reduced and oxidized glutathione from whole blood.

Reduced levels of glutathione (γ-glutamylcysteinylglycine, GSH) and the ratio of GSH to glutathione disulfide (GSSG) can serve as important indicators of oxidative stress and disease risk. Measured concentrations of GSH and GSSG vary widely between laboratories, largely due to the instability of GSH during sample handling and variables arising from different analytical methods. We have developed a simple and sensitive liquid chromatography-tandem mass spectrometry (LC-MS/MS) method for measuring whole blood GSH and GSSG that minimizes preanalytic and analytic variability, reliably eliminates interference from ion suppression, and can easily be implemented in clinical laboratories. Samples were deproteinized with sulfosalicylic acid (SSA) and derivatized with N-ethylmaleimide (NEM) in a single preparative step, and the resulting supernatants combined with stable-isotope internal standards (GSH-(13)C, (15)N-NEM and GSSG-(13)C,(15)N), subjected to chromatographic separation using a Hypercarb column, and analyzed by MS/MS in the positive-ion mode. Results showed excellent linearity for both GSH and GSSG over the ranges of physiologic normal, with inter- and intra-assay CV's of 3.1-4.3% and accuracy between 95% and 101%. The lower limits of detection (LLOD) were 0.4µM for GSH and 0.1µM for GSSG and the lower limits of quantitation (LLOQ) were 1.5µM for GSH and 0.1µM for GSSG. Derivatized samples are stable for at least 3 years when stored at -80°C, and underivatized samples for at least 24h at either 4°C or room temperature. Reference intervals were determined for 59 control samples, and were (mean±SD): GSH 900±140µM; GSSG 1.17±0.43µM; GSH/GSSG 880±370.

Interaction between genetic and epigenetic variation defines gene expression patterns at the asthma-associated locus 17q12-q21 in lymphoblastoid cell lines.

Phenotypic variation results from variation in gene expression, which is modulated by genetic and/or epigenetic factors. To understand the molecular basis of human disease, interaction between genetic and epigenetic factors needs to be taken into account. The asthma-associated region 17q12-q21 harbors three genes, the zona pellucida binding protein 2 (ZPBP2), gasdermin B (GSDMB) and ORM1-like 3 (ORMDL3), that show allele-specific differences in expression levels in lymphoblastoid cell lines (LCLs) and CD4+ T cells. Here, we report a molecular dissection of allele-specific transcriptional regulation of the genes within the chromosomal region 17q12-q21 combining in vitro transfection, formaldehyde-assisted isolation of regulatory elements, chromatin immunoprecipitation and DNA methylation assays in LCLs. We found that a single nucleotide polymorphism rs4795397 influences the activity of ZPBP2 promoter in vitro in an allele-dependent fashion, and also leads to nucleosome repositioning on the asthma-associated allele. However, variable methylation of exon 1 of ZPBP2 masks the strong genetic effect on ZPBP2 promoter activity in LCLs. In contrast, the ORMDL3 promoter is fully unmethylated, which allows detection of genetic effects on its transcription. We conclude that the cis-regulatory effects on 17q12-q21 gene expression result from interaction between several regulatory polymorphisms and epigenetic factors within the cis-regulatory haplotype region.

Genetic control of susceptibility to carcinogen-induced colorectal cancer in mice: the Ccs3 and Ccs5 loci regulate different aspects of tumorigenesis.

Colorectal cancer (CRC) is a multistep disease that involves a two-way interaction between a complex genetic pre-disposition component, and a set of poorly understood extrinsic environmental factors. In mice, CRC can be induced by treatment with azoxymethane (AOM). Using a set of AcB/BcA recombinant congenic strains derived from CRC-susceptible A/J and CRC-resistant C57Bl/6J (B6) progenitors, we previously detected the Ccs3 locus (colon cancer susceptibility locus 3) as a major regulator of CRC susceptibility. Phenotyping of additional AcB/BcA strains for susceptibility to AOM-induced CRC has refined the Ccs3 interval to a 6.7 Mb segment on chromosome 3. In addition, the presence of intermediate susceptibility phenotypes in individual AcB/BcA strains suggested additional gene effects regulating CRC susceptibility in A/J and B6 strains. Those were investigated by linkage analysis and whole genome scanning in a set of 208 informative (B6 x A/J)F2 progeny, using tumor multiplicity as a quantitative measure of susceptibility. This analysis validated the important role of Ccs3 in regulating this trait, and additionally detected contribution from a second locus on the distal portion of chromosome 9 (LOD = 3.76), that was given the temporary designation of Ccs5. Ccs5 modulates tumor multiplicity in F2 animals bearing at least one A/J-derived susceptibility allele at Ccs3, with A/J-derived Ccs5 susceptibility alleles being inherited in a recessive manner. There is a strong additive effect of Ccs3 and Ccs5 on tumor multiplicity in F2 mice: mice doubly homozygotes for A/J or B6 alleles at Ccs3 and Ccs5 show tumor numbers similar to those of parental A/J and B6, respectively. Interestingly, the Ccs5 region overlaps several quantitative trait loci previously reported to regulate intestinal homeostasis and susceptibility to intestinal colitis in mice and humans. Our findings identify a novel two-locus system regulating CRC susceptibility in mice, of which the relevance to human CRC can now be tested experimentally.

RNA sequencing reveals the role of splicing polymorphisms in regulating human gene expression.

Expression levels of many human genes are under the genetic control of expression quantitative trait loci (eQTLs). Despite technological advances, the precise molecular mechanisms underlying most eQTLs remain elusive. Here, we use deep mRNA sequencing of two CEU individuals to investigate those mechanisms, with particular focus on the role of splicing control loci (sQTLs). We identify a large number of genes that are differentially spliced between the two samples and associate many of those differences with nearby single nucleotide polymorphisms (SNPs). Subsequently, we investigate the potential effect of splicing SNPs on eQTL control in general. We find a significant enrichment of alternative splicing (AS) events within a set of highly confident eQTL targets discovered in previous studies, suggesting a role of AS in regulating overall gene expression levels. Next, we demonstrate high correlation between the levels of mature (exonic) and unprocessed (intronic) RNA, implying that ∼75% of eQTL target variance can be explained by control at the level of transcription, but that the remaining 25% may be regulated co- or post-transcriptionally. We focus on eQTL targets with discordant mRNA and pre-mRNA expression patterns and use four examples: USMG5, MMAB, MRPL43, and OAS1, to dissect the exact downstream effects of the associated genetic variants.

Allele-specific chromatin remodeling in the ZPBP2/GSDMB/ORMDL3 locus associated with the risk of asthma and autoimmune disease.

Common SNPs in the chromosome 17q12-q21 region alter the risk for asthma, type 1 diabetes, primary biliary cirrhosis, and Crohn disease. Previous reports by us and others have linked the disease-associated genetic variants with changes in expression of GSDMB and ORMDL3 transcripts in human lymphoblastoid cell lines (LCLs). The variants also alter regulation of other transcripts, and this domain-wide cis-regulatory effect suggests a mechanism involving long-range chromatin interactions. Here, we further dissect the disease-linked haplotype and identify putative causal DNA variants via a combination of genetic and functional analyses. First, high-throughput resequencing of the region and genotyping of potential candidate variants were performed. Next, additional mapping of allelic expression differences in Yoruba HapMap LCLs allowed us to fine-map the basis of the cis-regulatory differences to a handful of candidate functional variants. Functional assays identified allele-specific differences in nucleosome distribution, an allele-specific association with the insulator protein CTCF, as well as a weak promoter activity for rs12936231. Overall, this study shows a common disease allele linked to changes in CTCF binding and nucleosome occupancy leading to altered domain-wide cis-regulation. Finally, a strong association between asthma and cis-regulatory haplotypes was observed in three independent family-based cohorts (p = 1.78 x 10(-8)). This study demonstrates the requirement of multiple parallel allele-specific tools for the investigation of noncoding disease variants and functional fine-mapping of human disease-associated haplotypes.

Population genomics in a disease targeted primary cell model.

The common genetic variants associated with complex traits typically lie in noncoding DNA and may alter gene regulation in a cell type-specific manner. Consequently, the choice of tissue or cell model in the dissection of disease associations is important. We carried out an expression quantitative trait loci (eQTL) study of primary human osteoblasts (HOb) derived from 95 unrelated donors of Swedish origin, each represented by two independently derived primary lines to provide biological replication. We combined our data with publicly available information from a genome-wide association study (GWAS) of bone mineral density (BMD). The top 2000 BMD-associated SNPs (P < approximately 10(-3)) were tested for cis-association of gene expression in HObs and in lymphoblastoid cell lines (LCLs) using publicly available data and showed that HObs have a significantly greater enrichment (threefold) of converging cis-eQTLs as compared to LCLs. The top 10 BMD loci with SNPs showing strong cis-effects on gene expression in HObs (P = 6 x 10(-10) - 7 x 10(-16)) were selected for further validation using a staged design in two cohorts of Caucasian male subjects. All 10 variants were tested in the Swedish MrOS Cohort (n = 3014), providing evidence for two novel BMD loci (SRR and MSH3). These variants were then tested in the Rotterdam Study (n = 2090), yielding converging evidence for BMD association at the 17p13.3 SRR locus (P(combined) = 5.6 x 10(-5)). The cis-regulatory effect was further fine-mapped to the proximal promoter of the SRR gene (rs3744270, r(2) = 0.5, P = 2.6 x 10(-15)). Our results suggest that primary cells relevant to disease phenotypes complement traditional approaches for prioritization and validation of GWAS hits for follow-up studies.

Inhibition of transcription in Staphylococcus aureus by a primary sigma factor-binding polypeptide from phage G1.

The primary sigma factor of Staphylococcus aureus, sigma(SA), regulates the transcription of many genes, including several essential genes, in this bacterium via specific recognition of exponential growth phase promoters. In this study, we report the existence of a novel staphylococcal phage G1-derived growth inhibitory polypeptide, referred to as G1ORF67, that interacts with sigma(SA) both in vivo and in vitro and regulates its activity. Delineation of the minimal domain of sigma(SA) that is required for its interaction with G1ORF67 as amino acids 294 to 360 near the carboxy terminus suggests that the G1 phage-encoded anti-sigma factor may occlude the -35 element recognition domain of sigma(SA). As would be predicted by this hypothesis, the G1ORF67 polypeptide abolished both RNA polymerase core-dependent binding of sigma(SA) to DNA and sigma(SA)-dependent transcription in vitro. While G1ORF67 profoundly inhibits transcription when expressed in S. aureus cells in mode of action studies, our finding that G1ORF67 was unable to inhibit transcription when expressed in Escherichia coli concurs with its inability to inhibit transcription by the E. coli holoenzyme in vitro. These features demonstrate the selectivity of G1ORF67 for S. aureus RNA polymerase. We predict that G1ORF67 is one of the central polypeptides in the phage G1 strategy to appropriate host RNA polymerase and redirect it to phage reproduction.

Inherited disorders affecting mitochondrial function are associated with glutathione deficiency and hypocitrullinemia.

Disorders affecting mitochondria, including those that directly affect the respiratory chain function or result from abnormalities in branched amino acid metabolism (organic acidemias), have been shown to be associated with impaired redox balance. Almost all of the evidence underlying this conclusion has been obtained from studies on patient biopsies or animal models. Since the glutathione (iGSH) system provides the main protection against oxidative damage, we hypothesized that untreated oxidative stress in individuals with mitochondrial dysfunction would result in chronic iGSH deficiency. We confirm this hypothesis here in studies using high-dimensional flow cytometry (Hi-D FACS) and biochemical analysis of freshly obtained blood samples from patients with mitochondrial disorders or organic acidemias. T lymphocyte subsets, monocytes and neutrophils from organic acidemia and mitochondrial patients who were not on antioxidant supplements showed low iGSH levels, whereas similar subjects on antioxidant supplements showed normal iGSH. Measures of iROS levels in blood were insufficient to reveal the chronic oxidative stress in untreated patients. Patients with organic acidemias showed elevated plasma protein carbonyls, while plasma samples from all patients tested showed hypocitrullinemia. These findings indicate that measurements of iGSH in leukocytes may be a particularly useful biomarker to detect redox imbalance in mitochondrial disorders and organic acidemias, thus providing a relatively non-invasive means to monitor disease status and response to therapies. Furthermore, studies here suggest that antioxidant therapy may be useful for relieving the chronic oxidative stress that otherwise occurs in patients with mitochondrial dysfunction.

Antimicrobial drug discovery through bacteriophage genomics.

Over evolutionary time bacteriophages have developed unique proteins that arrest critical cellular processes to commit bacterial host metabolism to phage reproduction. Here, we apply this concept of phage-mediated bacterial growth inhibition to antibiotic discovery. We sequenced 26 Staphylococcus aureus phages and identified 31 novel polypeptide families that inhibited growth upon expression in S. aureus. The cellular targets for some of these polypeptides were identified and several were shown to be essential components of the host DNA replication and transcription machineries. The interaction between a prototypic pair, ORF104 of phage 77 and DnaI, the putative helicase loader of S. aureus, was then used to screen for small molecule inhibitors. Several compounds were subsequently found to inhibit both bacterial growth and DNA synthesis. Our results suggest that mimicking the growth-inhibitory effect of phage polypeptides by a chemical compound, coupled with the plethora of phages on earth, will yield new antibiotics to combat infectious diseases.