Sequence analysis of the cis-regulatory regions of the bithorax complex of Drosophila

The bithorax complex (BX-C) of Drosophila, one of two complexes that act as master regulators of the body plan of the fly, has now been entirely sequenced and comprises approximate to 315,000 bp, only 1.4% of which codes for protein. Analysis of this sequence reveals significantly overrepresented DNA motifs of unknown, as well as known, functions in the nonprotein-coding portion of the sequence. The following types of motifs in that portion are analyzed: (i) concatamers of mono-, di-, and trinucleotides; (ii) tightly clustered hexanucleotides (spaced less than or equal to 5 bases apart); (iii) direct and reverse repeats longer than 20 bp; and (iv) a number of motifs known from biochemical studies to play a role in the regulation of the BX-C. The hexanucleotide AGATAC is remarkably overrepresented and is surmised to play a role in chromosome pairing. The positions of sites of highly overrepresented motifs are plotted for those that occur at more than five sites in the sequence, when <0.5 case is expected. Expected values are based on a third-order Markov chain, which is the optimal order for representing the BXCALL sequence

Distributed BLAST in a grid computing context

The Basic Local Alignment Search Tool (BLAST) is one of the best known sequence comparison programs available in bioinformatics. It is used to compare query sequences to a set of target sequences, with the intention of finding similar sequences in the target set. Here, we present a distributed BLAST service which operates over a set of heterogeneous Grid resources and is made available through a Globus toolkit v.3 Grid service. This work has been carried out in the context of the BRIDGES project, a UK e-Science project aimed at providing a Grid based environment for biomedical research. Input consisting of multiple query sequences is partitioned into sub-jobs on the basis of the number of idle compute nodes available and then processed on these in batches. To achieve this, we have implemented our own Java-based scheduler which distributes sub-jobs across an array of resources utilizing a variety of local job scheduling systems

Distribution of chromosome 18 and X centric heterochromatin in the interphase nucleus of cultured human cells

In situ hybridization of human chromosome 18 and X-specific alphoid DNA-probes was performed in combination with three dimensional (3D) and two dimensional (2D) image analysis to study the interphase distribution of the centric heterochromatin (18c and Xc) of these chromosomes in cultured human cells. 3D analyses of 18c targets using confocal laser scanning microscopy indicated a nonrandom disposition in 73 amniotic fluid cell nuclei. The shape of these nuclei resembled rather flat cylinders or ellipsoids targets were preferentially arranged in a domain around the nuclear center, but close to or associated with the nuclear envelope. Within this domain, however, positionings of the two targets occurred independently from each other, i.e., the two targets were observed with similar frequencies at the same (upper or lower) side of the nuclear envelope as those on opposite sides. This result strongly argues against any permanent homologous association of 18c. A 2D analytical approach was used for the rapid evaluation of 18c positions in over 4000 interphase nuclei from normal male and female individuals, as well as individuals with trisomy 18 and Bloom's syndrome. In addition to epithelially derived amniotic fluid cells, investigated cell types included in vitro cultivated fibroblastoid cells established from fetal lung tissue and skin-derived fibroblasts. In agreement with the above 3D observations 18c targets were found significantly closer (P < 0.01) to the center of the 2D nuclear image (CNI) and to each other in all these cultures compared to a random distribution derived from corresponding ellipsoid or cylinder model nuclei. For comparison, a chromosome X-specific alphoid DNA probe was used to investigate the 2D distribution of chromosome X centric heterochromatin in the same cell types. Two dimensional Xc-Xc and Xc-CNI distances fit a random distribution in diploid normal and Bloom's syndrome nuclei, as well as in nuclei with trisomy X. The different distributions of 18c and Xc targets were confirmed by the simultaneous staining of these targets in different colors within individual nuclei using a double in situ hybridization approach

Figure 5: Number of WUs under processing at a particular time.

Development of high-throughput technologies, such as Next-generation sequencing, allows thousands of experiments to be performed simultaneously while reducing resource requirement. Consequently, a massive amount of experiment data is now rapidly generated. Nevertheless, the data are not readily usable or meaningful until they are further analysed and interpreted. Due to the size of the data, a high performance computer (HPC) is required for the analysis and interpretation. However, the HPC is expensive and difficult to access. Other means were developed to allow researchers to acquire the power of HPC without a need to purchase and maintain one such as cloud computing services and grid computing system. In this study, we implemented grid computing in a computer training center environment using Berkeley Open Infrastructure for Network Computing (BOINC) as a job distributor and data manager combining all desktop computers to virtualize the HPC. Fifty desktop computers were used for setting up a grid system during the off-hours. In order to test the performance of the grid system, we adapted the Basic Local Alignment Search Tools (BLAST) to the BOINC system. Sequencing results from Illumina platform were aligned to the human genome database by BLAST on the grid system. The result and processing time were compared to those from a single desktop computer and HPC. The estimated durations of BLAST analysis for 4 million sequence reads on a desktop PC, HPC and the grid system were 568, 24 and 5 days, respectively. Thus, the grid implementation of BLAST by BOINC is an efficient alternative to the HPC for sequence alignment. The grid implementation by BOINC also helped tap unused computing resources during the off-hours and could be easily modified for other available bioinformatics software

Chromosomal organization at the level of gene complexes

Metazoan genomes primarily consist of non-coding DNA in comparison to coding regions. Non-coding fraction of the genome contains cis-regulatory elements, which ensure that the genetic code is read properly at the right time and space during development. Regulatory elements and their target genes define functional landscapes within the genome, and some developmentally important genes evolve by keeping the genes involved in specification of common organs/tissues in clusters and are termed gene complex. The clustering of genes involved in a common function may help in robust spatio-temporal gene expression. Gene complexes are often found to be evolutionarily conserved, and the classic example is the hox complex. The evolutionary constraints seen among gene complexes provide an ideal model system to understand cis and trans-regulation of gene function. This review will discuss the various characteristics of gene regulatory modules found within gene complexes and how they can be characterized