GenomeCAT: a versatile tool for the analysis and integrative visualization of DNA copy number variants.

BACKGROUND: The analysis of DNA copy number variants (CNV) has increasing impact in the field of genetic diagnostics and research. However, the interpretation of CNV data derived from high resolution array CGH or NGS platforms is complicated by the considerable variability of the human genome. Therefore, tools for multidimensional data analysis and comparison of patient cohorts are needed to assist in the discrimination of clinically relevant CNVs from others. RESULTS: We developed GenomeCAT, a standalone Java application for the analysis and integrative visualization of CNVs. GenomeCAT is composed of three modules dedicated to the inspection of single cases, comparative analysis of multidimensional data and group comparisons aiming at the identification of recurrent aberrations in patients sharing the same phenotype, respectively. Its flexible import options ease the comparative analysis of own results derived from microarray or NGS platforms with data from literature or public depositories. Multidimensional data obtained from different experiment types can be merged into a common data matrix to enable common visualization and analysis. All results are stored in the integrated MySQL database, but can also be exported as tab delimited files for further statistical calculations in external programs. CONCLUSIONS: GenomeCAT offers a broad spectrum of visualization and analysis tools that assist in the evaluation of CNVs in the context of other experiment data and annotations. The use of GenomeCAT does not require any specialized computer skills. The various R packages implemented for data analysis are fully integrated into GenomeCATs graphical user interface and the installation process is supported by a wizard. The flexibility in terms of data import and export in combination with the ability to create a common data matrix makes the program also well suited as an interface between genomic data from heterogeneous sources and external software tools. Due to the modular architecture the functionality of GenomeCAT can be easily extended by further R packages or customized plug-ins to meet future requirements.

Comprehensive genotyping and clinical characterisation reveal 27 novel NKX2-1 mutations and expand the phenotypic spectrum.

BACKGROUND: NKX2-1 encodes a transcription factor with large impact on the development of brain, lung and thyroid. Germline mutations of NKX2-1 can lead to dysfunction and malformations of these organs. Starting from the largest coherent collection of patients with a suspected phenotype to date, we systematically evaluated frequency, quality and spectrum of phenotypic consequences of NKX2-1 mutations. METHODS: After identifying mutations by Sanger sequencing and array CGH, we comprehensively reanalysed the phenotype of affected patients and their relatives. We employed electrophoretic mobility shift assay (EMSA) to detect alterations of NKX2-1 DNA binding. Gene expression was monitored by means of in situ hybridisation and compared with the expression level of MBIP, a candidate gene presumably involved in the disorders and closely located in close genomic proximity to NKX2-1. RESULTS: Within 101 index patients, we detected 17 point mutations and 10 deletions. Neurological symptoms were the most consistent finding (100%), followed by lung affection (78%) and thyroidal dysfunction (75%). Novel symptoms associated with NKX2-1 mutations comprise abnormal height, bouts of fever and cardiac septum defects. In contrast to previous reports, our data suggest that missense mutations in the homeodomain of NKX2-1 not necessarily modify its DNA binding capacity and that this specific type of mutations may be associated with mild pulmonary phenotypes such as asthma. Two deletions did not include NKX2-1, but MBIP, whose expression spatially and temporarily coincides with NKX2-1 in early murine development. CONCLUSIONS: The high incidence of NKX2-1 mutations strongly recommends the routine screen for mutations in patients with corresponding symptoms. However, this analysis should not be confined to the exonic sequence alone, but should take advantage of affordable NGS technology to expand the target to adjacent regulatory sequences and the NKX2-1 interactome in order to maximise the yield of this diagnostic effort.

Distribution of segmental duplications in the context of higher order chromatin organisation of human chromosome 7.

BACKGROUND: Segmental duplications (SDs) are not evenly distributed along chromosomes. The reasons for this biased susceptibility to SD insertion are poorly understood. Accumulation of SDs is associated with increased genomic instability, which can lead to structural variants and genomic disorders such as the Williams-Beuren syndrome. Despite these adverse effects, SDs have become fixed in the human genome. Focusing on chromosome 7, which is particularly rich in interstitial SDs, we have investigated the distribution of SDs in the context of evolution and the three dimensional organisation of the chromosome in order to gain insights into the mutual relationship of SDs and chromatin topology. RESULTS: Intrachromosomal SDs preferentially accumulate in those segments of chromosome 7 that are homologous to marmoset chromosome 2. Although this formerly compact segment has been re-distributed to three different sites during primate evolution, we can show by means of public data on long distance chromatin interactions that these three intervals, and consequently the paralogous SDs mapping to them, have retained their spatial proximity in the nucleus. Focusing on SD clusters implicated in the aetiology of the Williams-Beuren syndrome locus we demonstrate by cross-species comparison that these SDs have inserted at the borders of a topological domain and that they flank regions with distinct DNA conformation. CONCLUSIONS: Our study suggests a link of nuclear architecture and the propagation of SDs across chromosome 7, either by promoting regional SD insertion or by contributing to the establishment of higher order chromatin organisation themselves. The latter could compensate for the high risk of structural rearrangements and thus may have contributed to their evolutionary fixation in the human genome.

Subtelomeric 6p monosomy and 12q trisomy in a patient with a 46,XX,der(6)t(6;12)(p25.3;q24.31) karyotype: Phenotypic overlap with Mutchinick syndrome.

We report on a patient with partial monosomy 6p and partial trisomy 12q identified by fluorescent in situ hybridization (FISH) and array-based comparative genomic hybridization (aCGH). She had a complex phenotype characterized by mental retardation (MR), psychomotor developmental delay, speech disorder, hypertelorism, eye anomalies, hearing loss, low-set malformed ears, thin upper lip, heart defect, clinodactyly, pes valgus, and skeletal anomalies. There is phenotypic overlap between our case and Mutchinick syndrome. This is the first report of a combined partial monosomy 6p and partial trisomy 12q due to an unbalanced translocation between subtelomeric regions of these chromosomes.

Genome-Wide Analysis of Interchromosomal Interaction Probabilities Reveals Chained Translocations and Overrepresentation of Translocation Breakpoints in Genes in a Cutaneous T-Cell Lymphoma Cell Line.

In classical models of tumorigenesis, the accumulation of tumor promoting chromosomal aberrations is described as a gradual process. Next-generation sequencing-based methods have recently revealed complex patterns of chromosomal aberrations, which are beyond explanation by these classical models of karyotypic evolution of tumor genomes. Thus, the term chromothripsis has been introduced to describe a phenomenon, where temporarily and spatially confined genomic instability results in dramatic chromosomal rearrangements limited to segments of one or a few chromosomes. Simultaneously arising and misrepaired DNA double-strand breaks are also the cause of another phenomenon called chromoplexy, which is characterized by the presence of chained translocations and interlinking deletion bridges involving several chromosomes. In this study, we demonstrate the genome-wide identification of chromosomal translocations based on the analysis of translocation-associated changes in spatial proximities of chromosome territories on the example of the cutaneous T-cell lymphoma cell line Se-Ax. We have used alterations of intra- and interchromosomal interaction probabilities as detected by genome-wide chromosome conformation capture (Hi-C) to infer the presence of translocations and to fine-map their breakpoints. The outcome of this analysis was subsequently compared to datasets on DNA copy number alterations and gene expression. The presence of chained translocations within the Se-Ax genome, partly connected by intervening deletion bridges, indicates a role of chromoplexy in the etiology of this cutaneous T-cell lymphoma. Notably, translocation breakpoints were significantly overrepresented in genes, which highlight gene-associated biological processes like transcription or other gene characteristics as a possible cause of the observed complex rearrangements. Given the relevance of chromosomal aberrations for basic and translational research, genome-wide high-resolution analysis of structural chromosomal aberrations will gain increasing importance.

Inactivation of RUNX3/p46 Promotes Cutaneous T-Cell Lymphoma.

The key role of RUNX3 in physiological T-cell differentiation has been extensively documented. However, information on its relevance for the development of human T-cell lymphomas or leukemias is scarce. Here, we show that alterations of RUNX3 by either heterozygous deletion or methylation of its distal promoter can be observed in the tumor cells of 15 of 21 (71%) patients suffering from Sézary syndrome, an aggressive variant of cutaneous T-cell lymphoma. As a consequence, mRNA levels of RUNX3/p46, the isoform controlled by the distal promoter, are significantly lower in Sézary syndrome tumor cells. Re-expression of RUNX3/p46 reduces cell viability and promotes apoptosis in a RUNX3/p46low cell line of cutaneous T-cell lymphoma. Based on this, we present evidence that RUNX3 can act as a tumor suppressor in a human T-cell malignancy and suggest that this effect is predominantly mediated through transcripts from its distal promoter, in particular RUNX3/p46.

Genomic loss of the putative tumor suppressor gene E2A in human lymphoma.

The transcription factor E2A is essential for lymphocyte development. In this study, we describe a recurrent E2A gene deletion in at least 70% of patients with Sézary syndrome (SS), a subtype of T cell lymphoma. Loss of E2A results in enhanced proliferation and cell cycle progression via derepression of the protooncogene MYC and the cell cycle regulator CDK6. Furthermore, by examining the gene expression profile of SS cells after restoration of E2A expression, we identify several E2A-regulated genes that interfere with oncogenic signaling pathways, including the Ras pathway. Several of these genes are down-regulated or lost in primary SS tumor cells. These data demonstrate a tumor suppressor function of E2A in human lymphoid cells and could help to develop new treatment strategies for human lymphomas with altered E2A activity.

An integrative approach for analyzing the interplay of genetic and epigenetic changes in tumors.

The accumulation of chromosomal aberrations is a characteristic feature of tumor development. However, an understanding of tumorigenesis that assumes that changes in DNA copy number always cause equivalent changes in the corresponding RNA and protein levels is an oversimplification and completely ignores the individual genetic and epigenetic context in which an aberration has to be evaluated. We present a brief introduction to various techniques dedicated to the genome-wide analysis of genetic and epigenetic changes, and illustrate how complementary information derived from these various DNA array-based technologies can lead to a better understanding of the consequences of chromosomal aberrations.