MicroRNAs (miRNAs) are approximately 22 nt-long non-coding RNA molecules, believed to play important roles in gene regulation. We present a comprehensive analysis of the conservation and clustering patterns of known miRNAs in human. We show that human miRNA gene clustering is significantly higher than expected at random. A total of 37% of the known human miRNA genes analyzed in this study appear in clusters of two or more with pairwise chromosomal distances of at most 3000 nt. Comparison of the miRNA sequences with their homologs in four other organisms reveals a typical conservation pattern, persistent throughout the clusters. Furthermore, we show enrichment in the typical conservation patterns and other miRNA-like properties in the vicinity of known miRNA genes, compared with random genomic regions. This may imply that additional, yet unknown, miRNAs reside in these regions, consistent with the current recognition that there are overlooked miRNAs. Indeed, by comparing our predictions with cloning results and with identified miRNA genes in other mammals, we corroborate the predictions of 18 additional human miRNA genes in the vicinity of the previously known ones. Our study raises the proportion of clustered human miRNAs that are <3000 nt apart to 42%. This suggests that the clustering of miRNA genes is higher than currently acknowledged, alluding to its evolutionary and functional implications.

Heterochromatin protein 1 (HP1) was originally described as a non-histone chromosomal protein and is required for transcriptional gene silencing and the formation of heterochromatin. Although it is localized primarily at pericentric heterochromatin, a scattered distribution over a large number of euchromatic loci is also evident. Here, we provide evidence that Drosophila HP1 is essential for the maintenance of active transcription of euchromatic genes functionally involved in cell-cycle progression, including those required for DNA replication and mitosis. Depletion of HP1 in proliferating embryonic cells caused aberrant progression of the cell cycle at S phase and G2/M phase, linked to aberrant chromosome segregation, cytokinesis, and an increase in apoptosis. The chromosomal distribution of Aurora B, and the level of phosphorylation of histone H3 serine 10 were also altered in the absence of HP1. Using chromatin immunoprecipitation analysis, we further demonstrate that the promoters of a number of cell-cycle regulator genes are bound to HP1, supporting a direct role for HP1 in their active transcription. Overall, our data suggest that HP1 is essential for the maintenance of cell-cycle progression and the transcription of cell-cycle regulatory genes. The results also support the view that HP1 is a positive regulator of transcription in euchromatin.

TFIIE and TFIIH are essential for the promoter opening and escape that occurs as RNA polymerase II transits into early elongation. XPB, a subunit of TFIIH, contains an ATP-dependent helicase activity that is used in both of these processes. Here, we show that the smaller beta subunit of TFIIE stimulates the XPB helicase and ATPase activities. The larger alpha subunit can use its known inhibitory activity to moderate the stimulation by the beta subunit. Regions of TFIIE beta required for the helicase stimulation were identified. Mutants were constructed that are defective in stimulating the XPB helicase but still allow intact TFIIE to bind and recruit XPB and TFIIH to form the pre-initiation complex. In a test for the functional significance of the stimulatory effect of TFIIE beta, these mutant forms of TFIIE were shown to be defective in a transcription assay on linear DNA. The data suggest that the beta subunit of TFIIE is an ATPase and helicase co-factor that can assist the XPB subunit of TFIIH during transcription initiation and the transition to early elongation, enhancing the potential diversity of regulatory targets.

Bloom's syndrome (BS) is an autosomal recessive human disorder characterized by genomic instability and a predisposition to a wide variety of cancers. The gene mutated in BS, BLM, encodes a protein containing three domains: an N-terminal domain whose function remains elusive, a helicase domain characterized by seven 'signature' motifs conserved in a wide range of helicases and a C-terminal extension that can be further divided into two sub-domains: RecQ-Ct and HRDC. The RecQ-Ct domain appears essential because two point-mutations altering highly conserved cysteine residues within this domain have been found in BS patients. We report herein that BLM contains a zinc ion. Modelling studies suggest that four conserved cysteine residues within the RecQ-Ct domain coordinate this zinc ion and subsequent mutagenesis studies further confirm this prediction. Biochemical and biophysical studies show that the ATPase, helicase and DNA binding activities of the mutants are severely modified. Structural analysis of both wild-type and mutant proteins reveal that alteration of cysteine residues does not significantly change the overall conformation. The observed defects in ATPase and helicase activities were inferred to result from a compromise of DNA binding. Our results implicate an important role of this zinc binding domain in both DNA binding and protein conformation. They could be pivotal for understanding the molecular basis of BS disease.

Pseudouridine (Psi) are frequently modified residues in RNA. In Eukarya, their formation is catalyzed by enzymes or by ribonucleoprotein complexes (RNPs) containing H/ACA snoRNAs. H/ACA sRNA and putative ORFs for H/ACA sRNP proteins (L7Ae, aCBF5, aNOP10 and aGAR1) were found in Archaea. Here, by using Pyrococcus abyssi recombinant proteins and an in vitro transcribed P.abyssi H/ACA sRNA, we obtained the first complete in vitro reconstitution of an active H/ACA RNP. Both L7Ae and the aCBF5 RNA:Psi synthase bind directly the sRNA; aCBF5 also interacts directly and independently with aNOP10 and aGAR1. Presence of aCBF5, aNOP10 and a U residue at the pseudouridylation site in the target RNA are required for RNA target recruitment. In agreement, we found that the aCBF5-aNOP10 pair is the minimal set of proteins needed for the formation of a particle active for pseudouridylation. However, particles more efficient in targeted pseudouridylation can be formed with the addition of proteins L7Ae and/or aGAR1. Although necessary for optimal activity, the conserved ACA motif in the sRNA was found to be not essential.

T-cells retain cell-type-specific programming for IL-2 inducibility through many rounds of division without being stimulated to transcribe the locus. To understand the layering of controls needed to poise this gene heritably for activation, we have used chromatin immunoprecipitation to map histone modifications across the murine IL2 locus, from -10.2 through +0.25 kb, in induction-competent and incompetent cells. In highly inducible EL4 T-lineage cells, stimulation with PMA/A23187 induced strong acetylation of histone H3 and H4, in parallel with transcriptional activation, from -4.6 through +0.25 kb. However, dimethylation of histone H3/K4 was already fully elevated across the same restricted domain before stimulation, with little change after stimulation. RNA polymerase II binding, in contrast, was only found at the known promoter region after stimulation. Similar patterns of histone modifications were seen also in normal IL-2-inducible T-lineage cells. However, neither acetylated histone H3, H4 nor dimethylated histone H3/K4 marking was detected, with or without stimulation, in expression-incompetent cells (NIH/3T3 or Scid.adh). These results identify a discrete new domain of IL2 regulatory sequence marked by dimethylated histone H3/K4 in expression-permissive T-cells even when they are not transcribing IL2, setting boundaries for histone H3 and H4 acetylation when the IL2 gene is transcriptionally activated.

Methylation of lysine 9 (K9) in the N-terminus tail of histone H3 (H3) in chromatin is associated with transcriptionally silenced genes and is mediated by histone methyltransferases. Murine G9a is a 1263 amino acid H3-K9 methyltransferase that possesses characteristic SET domain and ANK repeats. In this paper, we have used a series of green fluorescent protein-tagged deletion constructs to identify two nuclear localization signals (NLS), the first NLS embedded between amino acids 24 and 109 and the second between amino acids 394 and 401 of murine G9a. Our data show that both long and short G9a isoforms were capable of entering the nucleus to methylate chromatin. Full-length or N-terminus-deleted G9a isoforms were also catalytically active enzymes that methylated recombinant H3 or synthetic peptides representing the N-terminus tail of H3. In vitro methylation reactions using N-terminus tail peptides resulted in tri-methylation of K9 that remained processive, even in G9a enzymes that lacked an N-terminus region by deletion. Co-expression of G9a and H3 resulted in di- and tri-methylation of H3-K9, while siRNA-mediated knockdown of G9a in HeLa cells resulted in reduction of global H3-K9 di- and tri-methylation. A recombinant deletion mutant enzyme fused with maltose-binding protein (MBP-G9aDelta634) was used for steady-state kinetic analysis with various substrates and was compared with full-length G9a (G9aFL). Turnover numbers of MBP-G9aDelta634 for various substrates was approximately 3-fold less compared with G9aFL, while their Michaelis constants (K(m)) for recombinant H3 were similar. The K(AdoMet)m for MBP-G9aDelta634 was approximately 2.3-2.65 microM with various substrates. Catalytic efficiencies (kcat/K(m)) for both MBP-G9aDelta634 and G9aFL were similar, suggesting that the N-terminus is not essential for catalysis. Furthermore, mutation of conserved amino acids R1097A, W1103A, Y1120A, Y1138A and R1162A, or the metal binding C1168A in the catalytic region, resulted in catalytically impaired enzymes, thereby confirming the involvement of the C-terminus of G9a in catalysis. Thus, distinct domains modulate nuclear targeting and catalytic functions of G9a.

Abnormal DNA methylation is observed in certain promoters of neoplastic cells, although the likelihood of methylation for each individual promoter varies. Simultaneous analysis of many promoters in the same sample can allow use of statistical methods for identification of neoplasia. Here we describe an assay for such analysis, based on digestion of genomic DNA with methylation-sensitive restriction enzyme and multiplexed PCR with gene-specific primers (MSRE-PCR). MSRE-PCR includes extensive digestion of genomic DNA (uncut fragments cannot be identified by PCR), can be applied to dilute samples (<1 pg/microl), requires limited amount of starting material (42 pg or genomic equivalent of seven cells) and can identify methylation in a heterogeneous mix containing <2% of cells with methylated fragments. When applied to 53 promoters of breast cancer cell lines MCF-7, MDA-MB-231 and T47D, MSRE-PCR correctly identified the methylation status of genes analyzed by other techniques. For selected genes results of MSRE-PCR were confirmed by methylation-specific PCR and bisulfite sequencing. The assay can be configured for any number of desired targets in any user-defined set of genes.

The functional and dynamical properties of the human topoisomerase I Thr718Ala mutant have been compared to that of the wild-type enzyme using functional assays and molecular dynamics (MD) simulations. At physiological ionic strength, the cleavage and religation rates, evaluated on oligonucleotides containing the preferred topoisomerase I DNA sequence, are almost identical for the wild-type and the mutated enzymes, as is the cleavage/religation equilibrium. On the other hand, the Thr718Ala mutant shows a decreased efficiency in a DNA plasmid relaxation assay. The MD simulation, carried out on the enzyme complexed with its preferred DNA substrate, indicates that the mutant has a different dynamic behavior compared to the wild-type enzyme. Interestingly, no changes are observed in the proximity of the mutation site, whilst a different flexibility is detected in regions contacting the DNA scissile strand, such as the linker and the V-shaped alpha helices. Taken together, the functional and simulation results indicate a direct communication between the mutation site and regions located relatively far away, such as the linker domain, that with their altered flexibility confer a reduced DNA relaxation efficiency. These results provide evidence that the comprehension of the topoisomerase I dynamical properties are an important element in the understanding of its complex catalytic cycle.

A semi-quantitative real-time RT-PCR assay was designed to measure gonococcal pilin antigenic variation (SQ-PCR Av assay). This assay employs 17 hybridization probe sets that quantitate subpopulations of pilin transcripts carrying different silent pilin copy sequences and one set that detects total pilE transcript levels. Mixtures of a DNA standard carrying the silent copy being detected and a clone encoding the starting pilE sequence, which is the majority pilE template, provided amplification curves that closely matched the experimental data and allowed an analysis of the contribution of different silent pilin copies to variation. The SQ-PCR Av assay was verified using DNA sequence analysis to demonstrate that this methodology allowed an accurate analysis of pilin variation. Both assays showed that with a specific starting pilE sequence, only a subset of the silent pilin copies recombine into pilE at a detectable level, and that this limited subset was reproducibly detected in replicate cultures. When an isogenic pilE sequence variant was examined using both assays, a new subset of silent copy sequences were detected recombining into pilE and the overall frequency of variation was increased. Thus, the parental pilE sequence influences the frequency of variation and the repertoire of pilin variants produced.

The X-ray crystal structure of RadB from Thermococcus kodakaraensis KOD1, an archaeal homologue of the RecA/Rad51 family proteins, have been determined in two crystal forms. The structure represents the core ATPase domain of the RecA/Rad51 proteins. Two independent molecules in the type 1 crystal were roughly related by 7-fold screw symmetry whereas non-crystallographic 2-fold symmetry was observed in the type 2 crystal. The dimer structure in the type 1 crystal is extended to construct a helical assembly, which resembles the filamentous structures reported for other RecA/Rad51 proteins. The molecular interface in the type 1 dimer is formed by facing a basic surface patch of one monomer to an acidic one of the other. The empty ATP binding pocket is located at the interface and barely concealed from the outside similarly to that in the active form of the RecA filament. The model assembly has a positively charged belt on one surface bordering the helical groove suitable for facile binding of DNA. Electron microscopy has revealed that, in the absence of ATP and DNA, RadB forms a filament with a similar diameter to that of the hypothetical assembly, although its helical properties were not confirmed.

Identification of tissue-specific gene regulatory networks can yield insights into the molecular basis of a tissue's development, function and pathology. Here, we present a computational approach designed to identify potential regulatory target genes of photoreceptor cell-specific transcription factors (TFs). The approach is based on the hypothesis that genes related to the retina in terms of expression, disease and/or function are more likely to be the targets of retina-specific TFs than other genes. A list of genes that are preferentially expressed in retina was obtained by integrating expressed sequence tag, SAGE and microarray datasets. The regulatory targets of retina-specific TFs are enriched in this set of retina-related genes. A Bayesian approach was employed to integrate information about binding site location relative to a gene's transcription start site. Our method was applied to three retina-specific TFs, CRX, NRL and NR2E3, and a number of potential targets were predicted. To experimentally assess the validity of the bioinformatic predictions, mobility shift, transient transfection and chromatin immunoprecipitation assays were performed with five predicted CRX targets, and the results were suggestive of CRX regulation in 5/5, 3/5 and 4/5 cases, respectively. Together, these experiments strongly suggest that RP1, GUCY2D, ABCA4 are novel targets of CRX.

A CpG-containing DNA oligonucleotide functionalized with the 2-(N-formyl-N-methyl)aminoethyl thiophosphate protecting group (CpG ODN fma1555) was prepared from phosphoramidites 1a-d using solid-phase techniques. The oligonucleotide behaved as a prodrug by virtue of its conversion to the well-studied immunomodulatory CpG ODN 1555 through thermolytic cleavage of the 2-(N-formyl-N-methyl)aminoethyl thiophosphate protecting group. Such a conversion occurred at 37 degrees C with a half-time of 73 h. The immunostimulatory properties of CpG ODN fma1555 were evaluated in two in vivo assays, one of which consisted of mice challenged in the ear with live Leishmania major metacyclic promastigotes. Local intradermal administration of CpG ODN fma1555 was as effective as that of CpG ODN 1555 in reducing the size of Leishmania lesions over time. In a different infectious model, CpG ODN 1555 prevented the death of Tacaribe-infected mice (43% survival) when administered between day 0 and 3 post infection. Administration of CpG ODN fma1555 three days before infection resulted in improved immunoprotection (60-70% survival). Moreover, co-administration of CpG ODN fma1555 and CpG ODN 1555 in this model increased the window for therapeutic treatment against Tacaribe virus infection, and thus supports the use of thermolytic oligonucleotides as prodrugs in the effective treatment of infectious diseases.

GlyProt (http://www.glycosciences.de/glyprot/) is a web-based tool that enables meaningful N-glycan conformations to be attached to all the spatially accessible potential N-glycosylation sites of a known three-dimensional (3D) protein structure. The probabilities of physicochemical properties such as mass, accessible surface and radius of gyration are calculated. The purpose of this service is to provide rapid access to reliable 3D models of glycoproteins, which can subsequently be refined by using more elaborate simulations and validated by comparing the generated models with experimental data.

One-Block CYRCA is an automated procedure for identifying multiple-block alignments from single block queries (http://bioinfo.weizmann.ac.il/blocks/OneCYRCA). It is based on the LAMA and CYRCA block-to-block alignment methods. The procedure identifies whether the query blocks can form new multiple-block alignments (block sets) with blocks from a database or join pre-existing database block sets. Using pre-computed LAMA block alignments and CYRCA sets from the Blocks database reduces the computation time. LAMA and CYRCA are highly sensitive and selective methods that can augment many other sequence analysis approaches.

Our knowledge of metabolism can be represented as a network comprising several thousands of nodes (compounds and reactions). Several groups applied graph theory to analyse the topological properties of this network and to infer metabolic pathways by path finding. This is, however, not straightforward, with a major problem caused by traversing irrelevant shortcuts through highly connected nodes, which correspond to pool metabolites and co-factors (e.g. H2O, NADP and H+). In this study, we present a web server implementing two simple approaches, which circumvent this problem, thereby improving the relevance of the inferred pathways. In the simplest approach, the shortest path is computed, while filtering out the selection of highly connected compounds. In the second approach, the shortest path is computed on the weighted metabolic graph where each compound is assigned a weight equal to its connectivity in the network. This approach significantly increases the accuracy of the inferred pathways, enabling the correct inference of relatively long pathways (e.g. with as many as eight intermediate reactions). Available options include the calculation of the k-shortest paths between two specified seed nodes (either compounds or reactions). Multiple requests can be submitted in a queue. Results are returned by email, in textual as well as graphical formats (available in http://www.scmbb.ulb.ac.be/pathfinding/).

SNPAnalyzer is a software that performs four essential statistical analyses of SNPs in a common computational environment. It is composed of three main modules: (i) data manipulation, (ii) analysis and (iii) visualization. The data manipulation module is responsible for data input and output, and handles genotype, phenotype and genetic distance data. To ensure user convenience, the data format is simple. The analysis module performs statistical calculations and consists of four subcomponents: (i) Hardy-Weinberg equilibrium, (ii) Haplotype Estimation, (iii) linkage disequilibrium (LD) and (iv) quantitative trait locus analysis. The main feature of the analysis module is multiple implementations of different algorithms and indices for haplotype estimation and for LD analysis. This enables users to compare separate results generated by different algorithms, which help to avoid biased results acquired by applying a single statistical algorithm. The performance of all implemented algorithms has been validated using experimentally proven datasets. The visualization module presents most of the analyzed results as figures, rather than as simple text, which aids in the intuitive understanding of complex data. The SNPAnalyzer has been developed using C and C++ and is available at http://www.istech.info/istech/board/login_form.jsp.

The Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) is a relatively simple and inexpensive method for genotyping single nucleotide polymorphisms (SNPs). It requires minimal investment in instrumentation. Here, we describe a web application, 'SNP Cutter,' which designs PCR-RFLP assays on a batch of SNPs from the human genome. NCBI dbSNP rs IDs or formatted SNPs are submitted into the SNP Cutter which then uses restriction enzymes from a pre-selected list to perform enzyme selection. The program is capable of designing primers for either natural PCR-RFLP or mismatch PCR-RFLP, depending on the SNP sequence data. SNP Cutter generates the information needed to evaluate and perform genotyping experiments, including a PCR primers list, sizes of original amplicons and different allelic fragment after enzyme digestion. Some output data is tab-delimited, therefore suitable for database archiving. The SNP Cut-ter is available at http://bioinfo.bsd.uchicago.edu/SNP_cutter.htm.

SNPServer is a real-time flexible tool for the discovery of SNPs (single nucleotide polymorphisms) within DNA sequence data. The program uses BLAST, to identify related sequences, and CAP3, to cluster and align these sequences. The alignments are parsed to the SNP discovery software autoSNP, a program that detects SNPs and insertion/deletion polymorphisms (indels). Alternatively, lists of related sequences or pre-assembled sequences may be entered for SNP discovery. SNPServer and autoSNP use redundancy to differentiate between candidate SNPs and sequence errors. For each candidate SNP, two measures of confidence are calculated, the redundancy of the polymorphism at a SNP locus and the co-segregation of the candidate SNP with other SNPs in the alignment. SNPServer is available at http://hornbill.cspp.latrobe.edu.au/snpdiscovery.html.

We have developed a web tool, PupasView, for the selection of single nucleotide polymorphisms (SNPs) with potential phenotypic effect. PupasView constitutes an interactive environment in which functional information and population frequency data can be used as sequential filters over linkage disequilibrium parameters to obtain a final list of SNPs optimal for genotyping purposes. PupasView is the first resource that integrates phenotypic effects caused by SNPs at both the translational and the transcriptional level. PupasView retrieves SNPs that could affect conserved regions that the cellular machinery uses for the correct processing of genes (intron/exon boundaries or exonic splicing enhancers), predicted transcription factor binding sites and changes in amino acids in the proteins for which a putative pathological effect is calculated. The program uses the mapping of SNPs in the genome provided by Ensembl. PupasView will be of much help in studies of multifactorial disorders, where the use of functional SNPs will increase the sensitivity of the identification of the genes responsible for the disease. The PupasView web interface is accessible through http://pupasview.ochoa.fib.es and through http://www.pupasnp.org.