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

BACKGROUND: Sequence 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. RESULTS: In 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. CONCLUSIONS: IGGalaxy 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.

Similar recombination-activating gene (RAG) mutations result in similar immunobiological effects but in different clinical phenotypes.

BACKGROUND: V(D)J recombination takes place during lymphocyte development to generate a large repertoire of T- and B-cell receptors. Mutations in recombination-activating gene 1 (RAG1) and RAG2 result in loss or reduction of V(D)J recombination. It is known that different mutations in RAG genes vary in residual recombinase activity and give rise to a broad spectrum of clinical phenotypes. OBJECTIVE: We sought to study the immunologic mechanisms causing the clinical spectrum of RAG deficiency. METHODS: We included 22 patients with similar RAG1 mutations (c.519delT or c.368_369delAA) resulting in N-terminal truncated RAG1 protein with residual recombination activity but presenting with different clinical phenotypes. We studied precursor B-cell development, immunoglobulin and T-cell receptor repertoire formation, receptor editing, and B- and T-cell numbers. RESULTS: Clinically, patients were divided into 3 main categories: T(-)B(-) severe combined immunodeficiency, Omenn syndrome, and combined immunodeficiency. All patients showed a block in the precursor B-cell development, low B- and T-cell numbers, normal immunoglobulin gene use, limited B- and T-cell repertoires, and slightly impaired receptor editing. CONCLUSION: This study demonstrates that similar RAG mutations can result in similar immunobiological effects but different clinical phenotypes, indicating that the level of residual recombinase activity is not the only determinant for clinical outcome. We postulate a model in which the type and moment of antigenic pressure affect the clinical phenotypes of these patients.

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.

Recent advances in i-Gene tools and analysis: microarrays, next generation sequencing and mass spectrometry.

Recent advances in technology and associated methodology have made the current period one of the most exciting in molecular biology and medicine. Underlying these is an appreciation that modern research is driven by increasing large amounts of data being interpreted by interdisciplinary collaborative teams which are often geographically dispersed. The availability of cheap computing power, high speed informatics networks and high quality analysis software has been essential to this as has the application of modern quality assurance methodologies. In this review, we discuss the application of modern 'High-Throughput' molecular biological technologies such as 'Microarrays' and 'Next Generation Sequencing' to scientific and biomedical research as we have observed. Furthermore in this review, we also offer some guidance that enables the reader as to understand certain features of these as well as new strategies and help them to apply these i-Gene tools in their endeavours successfully. Collectively, we term this 'i-Gene Analysis'. We also offer predictions as to the developments that are anticipated in the near and more distant future.

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.

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.

Nasopharyngeal co-colonization with Staphylococcus aureus and Streptococcus pneumoniae in children is bacterial genotype independent.

Bacterial interference between Staphylococcus aureus and Streptococcus pneumoniae in the nasopharynx has been observed during colonization, which might have important clinical implications for the widespread use of pneumococcal conjugate vaccine in young children. This study aimed to determine whether the capacity of Staph. aureus to compete with Strep. pneumoniae is dependent on bacterial genotype. Demographic and microbiological determinants of carriage of specific genotypes of Staph. aureus in children were also studied. Children (n=3198) were sampled in the nasopharynx to detect carriage of Staph. aureus, Strep. pneumoniae and Neisseria meningitidis. Staph. aureus genotypes and pneumococcal sero- and genotypes were determined. Age, gender, zip code, active smoking and co-colonization with N. meningitidis or Strep. pneumoniae, both vaccine- and non-vaccine types, were not associated with colonization by specific Staph. aureus genotypes. Based on the whole-genome typing data obtained, there was no obvious correlation between staphylococcal and pneumococcal genotypes during co-colonization. Passive smoking showed a significant association (P=0.003) with carriage of a specific Staph. aureus cluster. This study suggests that there are no major differences between Staph. aureus clones (with different disease-invoking potential) in their capacity to compete with Strep. pneumoniae subtypes. Further studies should demonstrate whether differences in bacterial interference are due to more subtle genetic changes.

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

BACKGROUND: Accurate 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. RESULTS: We 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/. CONCLUSION: The HeatMapper allows an accessible and comprehensive visualization of the results of gene expression profiling and cluster analysis.

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

BACKGROUND: In 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. RESULTS: In order to facilitate this analysis we generated TF Target Mapper (Transcription Factor 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. alpha-globin and zeta-globin) and potentially novel GATA-1 gene targets. CONCLUSION: TF 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.

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 immuno histochemistry, 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 Ialpha, 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.

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 gamma-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.

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.

Prognostically useful gene-expression profiles in acute myeloid leukemia.

BACKGROUND: In patients with acute myeloid leukemia (AML) a combination of methods must be used to classify the disease, make therapeutic decisions, and determine the prognosis. However, this combined approach provides correct therapeutic and prognostic information in only 50 percent of cases. METHODS: We determined the gene-expression profiles in samples of peripheral blood or bone marrow from 285 patients with AML using Affymetrix U133A GeneChips containing approximately 13,000 unique genes or expression-signature tags. Data analyses were carried out with Omniviz, significance analysis of microarrays, and prediction analysis of microarrays software. Statistical analyses were performed to determine the prognostic significance of cases of AML with specific molecular signatures. RESULTS: Unsupervised cluster analyses identified 16 groups of patients with AML on the basis of molecular signatures. We identified the genes that defined these clusters and determined the minimal numbers of genes needed to identify prognostically important clusters with a high degree of accuracy. The clustering was driven by the presence of chromosomal lesions (e.g., t(8;21), t(15;17), and inv(16)), particular genetic mutations (CEBPA), and abnormal oncogene expression (EVI1). We identified several novel clusters, some consisting of specimens with normal karyotypes. A unique cluster with a distinctive gene-expression signature included cases of AML with a poor treatment outcome. CONCLUSIONS: Gene-expression profiling allows a comprehensive classification of AML that includes previously identified genetically defined subgroups and a novel cluster with an adverse prognosis.