Michael Moorhouse

Natural population dynamics and expansion of pathogenic clones of Staphylococcus aureus

The population structure of Staphylococcus aureus carried by healthy humans was determined using a large strain collection of nonclinical origin (n = 829). High-throughput amplified fragment length polymorphism (AFLP) analysis revealed 3 major and 2 minor genetic clusters of S. aureus, which were corroborated by multilocus sequence typing. Major AFLP cluster I comprised 44.4% of the carriage isolates and showed additional heterogeneity whereas major AFLP groups II and III presented 2 homogeneous clusters, including 47.3% of all carriage isolates. Coanalysis of invasive S. aureus strains and epidemic methicillin-resistant S. aureus (MRSA) revealed that all major clusters contained invasive and multiresistant isolates. However, clusters and subclusters with overrepresentation of invasive isolates were also identified. Bacteremia in elderly adults, for instance, was caused by a IVa cluster-derived strain significantly more often than by strains from other AFLP clusters. Furthermore, expansion of multiresistant clones or clones associated with skin disease (impetigo) was detected, which suggests that epidemic potential is present in pathogenic strains of S. aureus. In addition, the virulence gene encoding Panton-Valentine leukocidin was significantly enriched in S. aureus strains causing abscesses and arthritis in comparison with the carriage group. We provide evidence that essentially any S. aureus genotype carried by humans can transform into a life-threatening human pathogen but that certain clones are more virulent than others.

Delayed and accelerated aging share common longevity assurance mechanisms

Mutant dwarf and calorie-restricted mice benefit from healthy aging and unusually long lifespan. In contrast, mouse models for DNA repair-deficient progeroid syndromes age and die prematurely. To identify mechanisms that regulate mammalian longevity, we quantified the parallels between the genome-wide liver expression profiles of mice with those two extremes of lifespan. Contrary to expectation, we find significant, genome-wide expression associations between the progeroid and long-lived mice. Subsequent analysis of significantly over-represented biological processes revealed suppression of the endocrine and energy pathways with increased stress responses in both delayed and premature aging. To test the relevance of these processes in natural aging, we compared the transcriptomes of liver, lung, kidney, and spleen over the entire murine adult lifespan and subsequently confirmed these findings on an independent aging cohort. The majority of genes showed similar expression changes in all four organs, indicating a systemic transcriptional response with aging. This systemic response included the same biological processes that are triggered in progeroid and long-lived mice. However, on a genome-wide scale, transcriptomes of naturally aged mice showed a strong association to progeroid but not to long-lived mice. Thus, endocrine and metabolic changes are indicative of “survival” responses to genotoxic stress or starvation, whereas genome-wide associations in gene expression with natural aging are indicative of biological age, which may thus delineate pro- and anti-aging effects of treatments aimed at health-span extension.

Transcriptome analysis reveals cyclobutane pyrimidine dimers as a major source of UV-induced DNA breaks

Photolyase transgenic mice have opened new avenues to improve our understanding of the cytotoxic effects of ultraviolet (UV) light on skin by providing a means to selectively remove either cyclobutane pyrimidine dimers (CPDs) or pyrimidine (6‐4) pyrimidone photoproducts. Here, we have taken a genomics approach to delineate pathways through which CPDs might contribute to the harmful effects of UV exposure. We show that CPDs, rather than other DNA lesions or damaged macromolecules, comprise the principal mediator of the cellular transcriptional response to UV. The most prominent pathway induced by CPDs is that associated with DNA double‐strand break (DSB) signalling and repair. Moreover, we show that CPDs provoke accumulation of γ‐H2AX, P53bp1 and Rad51 foci as well as an increase in the amount of DSBs, which coincides with accumulation of cells in S phase. Thus, conversion of unrepaired CPD lesions into DNA breaks during DNA replication may comprise one of the principal instigators of UV‐mediated cytotoxicity.

Differential Role of Basal Keratinocytes in UV-Induced Immunosuppression and Skin Cancer

ABSTRACT Cyclobutane pyrimidine dimers (CPDs) and 6-4 photoproducts (6-4PPs) comprise major UV-induced photolesions. If left unrepaired, these lesions can induce mutations and skin cancer, which is facilitated by UV-induced immunosuppression. Yet the contribution of lesion and cell type specificity to the harmful biological effects of UV exposure remains currently unclear. Using a series of photolyase-transgenic mice to ubiquitously remove either CPDs or 6-4PPs from all cells in the mouse skin or selectively from basal keratinocytes, we show that the majority of UV-induced acute effects to require the presence of CPDs in basal keratinocytes in the mouse skin. At the fundamental level of gene expression, CPDs induce the expression of genes associated with repair and recombinational processing of DNA damage, as well as apoptosis and a response to stress. At the organismal level, photolyase-mediated removal of CPDs, but not 6-4PPs, from the genome of only basal keratinocytes substantially diminishes the incidence of skin tumors; however, it does not affect the UVB-mediated immunosuppression. Taken together, these findings reveal a differential role of basal keratinocytes in these processes, providing novel insights into the skins acute and chronic responses to UV in a lesion- and cell-type-specific manner.

DNA Microarray Analysis for Human Congenital Heart Disease

Right ventricular hypertrophy and failure are prominent features in cyanotic congenital heart disease, tetralogy of Fallot (TF). Patients with TF require primary cardiac surgery at a very young age. To gain insight into the underlying molecular mechanisms of right ventricular hypertrophy and to identify gene(s) involved in TF, differential gene expression profile was assessed using expression-based microarray technology on right ventricular biopsies from young TF patients who underwent primary correction. By using quantitative immunohistochemistry, expression of vascular endothelial growth factor (VEGF), flk-1, and extracellular matrix (ECM) proteins (collagens and fibronectin) as well as vessel counts and myocyte cell size was evaluated in TF patients in relation to age-matched controls. Among 236 genes showing altered expression pattern in TF patients, VEGF (1.8-fold) and ECM markers were clearly upregulated (fibronectin, 2.4-fold; collagen Iα, 7.5-fold; and collagen III, 4.4-fold); flk-1 and most matrix metalloproteinases (MMPs) remained unchanged, except the levels of MMP-13 and-17 declined. Tissue inhibitors of metalloproteinases showed a downregulated pattern. Staining of VEGF in cardiomyocytes and of ECM proteins (fibronectin, collagen I and III) in interstitial as well as in perivascular area was increased (p<0.01) in TF patients. Morphometric analysis revealed enhanced vascular density (p<0.05) with unchanged wall thickness and enlarged myocyte cross-sectional areas (p<0.01) with linear correlation (p<0.01) with the age in TF-1 patients. We conclude that the upregulation of genes encoding VEGF and ECM proteins are the key events contributing to right ventricular hypertrophy and stunted angiogenesis in patients with TF.

HeatMapper: powerful combined visualization of gene expression profile correlations, genotypes, phenotypes and sample characteristics

BackgroundAccurate interpretation of data obtained by unsupervised analysis of large scale expression profiling studies is currently frequently performed by visually combining sample-gene heatmaps and sample characteristics. This method is not optimal for comparing individual samples or groups of samples. Here, we describe an approach to visually integrate the results of unsupervised and supervised cluster analysis using a correlation plot and additional sample metadata.ResultsWe have developed a tool called the HeatMapper that provides such visualizations in a dynamic and flexible manner and is available from http://www.erasmusmc.nl/hematologie/heatmapper/.ConclusionThe HeatMapper allows an accessible and comprehensive visualization of the results of gene expression profiling and cluster analysis.

ImmunoGlobulin galaxy (IGGalaxy) for simple determination and quantitation of immunoglobulin heavy chain rearrangements from NGS.

BackgroundSequence analysis of immunoglobulin heavy chain (IGH) gene rearrangements and frequency analysis is a powerful tool for studying the immune repertoire, immune responses and immune dysregulation in health and disease. The challenge is to provide user friendly, secure and reproducible analytical services that are available for both small and large laboratories which are determining VDJ repertoire using NGS technology.ResultsIn this study we describe ImmunoGlobulin Galaxy (IGGalaxy)- a convenient web based application for analyzing next-generation sequencing results and reporting IGH gene rearrangements for both repertoire and clonality studies. IGGalaxy has two analysis options one using the built in igBLAST algorithm and the second using output from IMGT; in either case repertoire summaries for the B-cell populations tested are available. IGGalaxy supports multi-sample and multi-replicate input analysis for both igBLAST and IMGT/HIGHV-QUEST. We demonstrate the technical validity of this platform using a standard dataset, S22, used for benchmarking the performance of antibody alignment utilities with a 99.9 % concordance with previous results. Re-analysis of NGS data from our samples of RAG-deficient patients demonstrated the validity and user friendliness of this tool.ConclusionsIGGalaxy provides clinical researchers with detailed insight into the repertoire of the B-cell population per individual sequenced and between control and pathogenic genomes. IGGalaxy was developed for 454 NGS results but is capable of analyzing alternative NGS data (e.g. Illumina, Ion Torrent). We demonstrate the use of a Galaxy virtual machine to determine the VDJ repertoire for reference data and from B-cells taken from immune deficient patients. IGGalaxy is available as a VM for download and use on a desktop PC or on a server.

Absence of common somatic alterations in genes on 1p and 19q in oligodendrogliomas

A common and histologically well defined subtype of glioma are the oligodendroglial brain tumors. Approximately 70% of all oligodendrogliomas have a combined loss of the entire 1p and 19q chromosomal arms. This remarkably high frequency suggests that the remaining arms harbor yet to be identified tumor suppressor genes. Identification of these causal genetic changes in oligodendrogliomas is important because they form direct targets for treatment. In this study we therefore performed targeted resequencing of all exons, microRNAs, splice sites and promoter regions residing on 1p and 19q on 7 oligodendrogliomas and 4 matched controls. Only one missense mutation was identified in a single sample in the ARHGEF16 gene. This mutation lies within- and disrupts the conserved PDZ binding domain. No similar ARHGEF16 mutations or deletions were found in a larger set of oligodendrogliomas. The absence of common somatic changes within genes located on 1p and 19q in three out of four samples indicates that no additional “second hit” is required to drive oncogenic transformation on either chromosomal arm.

Extensive Variation in Gene Copy Number at the Killer Immunoglobulin-Like Receptor Locus in Humans

Killer immunoglobulin-like receptors (KIRs) are involved in the regulation of natural killer cell cytotoxicity. Within the human genome seventeen KIR genes are present, which all contain a large number of allelic variants. The high level of homology among KIR genes has hampered KIR genotyping in larger cohorts, and determination of gene copy number variation (CNV) has been difficult. We have designed a multiplex ligation-dependent probe amplification (MLPA) technique for genotyping and CNV determination in one single assay and validated the results by next-generation sequencing and with a KIR gene-specific short tandem repeat assay. In this way, we demonstrate in a cohort of 120 individuals a high level of CNV for all KIR genes except for the framework genes KIR3DL3 and KIR3DL2. Application of our MLPA assay in segregation analyses of families from the Centre d’Etude du Polymorphisme Humaine, previously KIR-genotyped by classical techniques, confirmed an earlier reported duplication and resulted in the identification of a novel duplication event in one of these families. In summary, our KIR MLPA assay allows rapid and accurate KIR genotyping and CNV detection, thus rendering improved transplantation programs and oncology treatment feasible, and enables more detailed studies on the role of KIRs in human (auto)immunity and infectious disease.

TF Target Mapper: a BLAST search tool for the identification of Transcription Factor target genes.

BackgroundIn the current era of high throughput genomics a major challenge is the genome-wide identification of target genes for specific transcription factors. Chromatin immunoprecipitation (ChIP) allows the isolation of in vivo binding sites of transcription factors and provides a powerful tool for examining gene regulation. Crosslinked chromatin is immunoprecipitated with antibodies against specific transcription factors, thus enriching for sequences bound in vivo by these factors in the immunoprecipitated DNA. Cloning and sequencing the immunoprecipitated sequences allows identification of transcription factor target genes. Routinely, thousands of such sequenced clones are used in BLAST searches to map their exact location in the genome and the genes located in the vicinity. These genes represent potential targets of the transcription factor of interest. Such bioinformatics analysis is very laborious if performed manually and for this reason there is a need for developing bioinformatic tools to automate and facilitate it.ResultsIn order to facilitate this analysis we generated TF Target Mapper (T ranscription F actor Target Mapper). TF Target Mapper is a BLAST search tool allowing rapid extraction of annotated information on genes around each hit. It combines sequence cleaning/filtering, pattern searching and BLAST searches with extraction of information on genes located around each BLAST hit and comparisons of the output list of genes or gene ontology IDs with user-implemented lists. We successfully applied and tested TF Target Mapper to analyse sequences bound in vivo by the transcription factor GATA-1. We show that TF Target Mapper efficiently extracted information on genes around ChIPed sequences, thus identifying known (e.g. α-globin and ζ-globin) and potentially novel GATA-1 gene targets.ConclusionTF Target Mapper is a very efficient BLAST search tool that allows the rapid extraction of annotated information on the genes around each hit. It can contribute to the comprehensive bioinformatic transcriptome/regulome analysis, by providing insight into the mechanisms of action of specific transcription factors, thus helping to elucidate the pathways these factors regulate.