ClusTrack: Feature Extraction and Similarity Measures for Clustering of Genome-Wide Data Sets

Clustering is a popular technique for explorative analysis of data, as it can reveal subgroupings and similarities between data in an unsupervised manner. While clustering is routinely applied to gene expression data, there is a lack of appropriate general methodology for clustering of sequence-level genomic and epigenomic data, e.g. ChIP-based data. We here introduce a general methodology for clustering data sets of coordinates relative to a genome assembly, i.e. genomic tracks. By defining appropriate feature extraction approaches and similarity measures, we allow biologically meaningful clustering to be performed for genomic tracks using standard clustering algorithms. An implementation of the methodology is provided through a tool, ClusTrack, which allows fine-tuned clustering analyses to be specified through a web-based interface. We apply our methods to the clustering of occupancy of the H3K4me1 histone modification in samples from a range of different cell types. The majority of samples form meaningful subclusters, confirming that the definitions of features and similarity capture biological, rather than technical, variation between the genomic tracks. Input data and results are available, and can be reproduced, through a Galaxy Pages document at http://hyperbrowser.uio.no/hb/u/hb-superuser/p/clustrack. The clustering functionality is available as a Galaxy tool, under the menu option "Specialized analyzis of tracks", and the submenu option "Cluster tracks based on genome level similarity", at the Genomic HyperBrowser server: http://hyperbrowser.uio.no/hb/.     

A statistical model of ChIA-PET data for accurate detection of chromatin 3D interactions

Identification of three-dimensional (3D) interactions between regulatory elements across the genome is crucial to unravel the complex regulatory machinery that orchestrates proliferation and differentiation of cells. ChIA-PET is a novel method to identify such interactions, where physical contacts between regions bound by a specific protein are quantified using next-generation sequencing. However, determining the significance of the observed interaction frequencies in such datasets is challenging, and few methods have been proposed. Despite the fact that regions that are close in linear genomic distance have a much higher tendency to interact by chance, no methods to date are capable of taking such dependency into account. Here, we propose a statistical model taking into account the genomic distance relationship, as well as the general propensity of anchors to be involved in contacts overall. Using both real and simulated data, we show that the previously proposed statistical test, based on Fisher''s exact test, leads to invalid results when data are dependent on genomic distance. We also evaluate our method on previously validated cell-line specific and constitutive 3D interactions, and show that relevant interactions are significant, while avoiding over-estimating the significance of short nearby interactions.     

The Detection of Hamster Connexins: A Comparison of Expression Profiles with Wild-Type Mouse and the Cancer-Prone Min Mouse

The open reading frames of 17 connexins from Syrian hamster (using tissues) and 16 connexins from the Chinese hamster cell line V79, were fully (Cx30, Cx31, Cx37, Cx43 and Cx45) or partially sequenced. We have also detected, and partially sequenced, seven rat connexins that previously were unavailable. The expression of connexin genes was examined in some hamster organs and cultured hamster cells, and compared with wild-type mouse and the cancer-prone Min mouse. Although the expression patterns were similar for most organs and connexins in hamster and mouse, there were also some prominent differences (Cx29 and 30.3 in testis; Cx31.1 and 32 in eye; Cx46 in brain, kidney and testis; Cx47 in kidney). This suggests that some connexins have species-specific expression profiles. In contrast, there were minimal differences in expression profiles between wild type and Min mice. Species-specific expression profiles should be considered in attempts to make animal models of human connexin-associated diseases.     

The influence of fixation on the relative amount of cytoplasmic ribosomes in mouse epidermal basal keratinocytes

The nature and significance of so-called dark keratinocytes in the epidermis during chemical carcinogenesis is still a matter of concern and debate. Based on ultrastructural observations it has been suggested that dark cells most often are shrunken cells. Reports on skin carcinogenesis, however, claim that dark cells are a sign of ongoing tumor promotion and represent those stem cells in the epidermis from which the tumors originate. It is therefore important to find out whether these cells are simply injured and shrunken cells, or vital cells of great importance for carcinogenesis. Dark cells are assumed to be rich in ribosomes. There is evidence, however, that the observed number of dark cells is highly dependent on tissue fixation. In the present ultrastructural study, morphometric methods were used to compare the effects of two different fixation procedures on the amount of cytoplasmic ribosomes in dark cells from both untreated and carcinogen-treated hairless mouse epidermis. The results show that the ultrastructural features of both dark and clear cells vary considerably with different fixation procedures. In acetone-treated controls typical dark cells are only observed when the fixative has a lower osmotic activity than the plasma. With iso-osmolal fixation typical dark cells are not observed. After an abortive two-stage carcinogenesis treatment, in which a single application of 9,10-dimethyl-l,2-benzanthracene (DMBA) in acetone was followed by a single application of 12-O-tetradecanoyl-13-acetate (TPA) in acetone, signs of cell injury could be found after both fixation procedures. With DMBA/TPA and hypo-osmolal fixation the number of dark cells seemed to increase, whereas only signs of cell injury with occurrence of some heavily altered “clear cells” dominated the picture with iso-osmolal fixation. Morphometry showed that both the numerical and the volumetric densities of cytoplasmic ribosomes in basal keratinocytes varied most significantly with the fixation procedure used. The cytoplasmic volumes did not vary in a way that could explain these differences. One might therefore assume that the number of ribosomes depends on the fixative. Large swelling artifacts occurred when a fixative with low osmotic activity was used, leading to compression of neighboring cells. Hence, an increased ribosomal density reported previously in dark cells is probably related to such cell volume artifacts and does not reflect an actually increased quantity of ribosomes. With both fixation procedures, a single application of DMBA followed by one of TPA appeared to produce an increased number of ribosomes in basal keratinocytes. When hypo-osmolal fixation was used, however, treatment with DMBA/TPA did not influence the cytoplasmic volume or the numerical density of ribosomes, in dark cells. This might indicate that so-called dark keratinocytes following DMBA/TPA treatment are functionally inactive cells that appear more vulnerable than active cells to compression during hypo-osmolal fixation.