The carbohydrate response element (ChoRE) is a cis-acting sequence found in the promoters of genes induced transcriptionally by glucose. The ChoRE is composed of two E box-like motifs that are separated by 5 bp and is recognized by two basic helix-loop-helix/leucine zipper (bHLH/LZ) proteins, ChREBP and Mlx, which heterodimerize to bind DNA. In this study, we demonstrate that two ChREBP/Mlx heterodimers interact to stabilize binding to the tandem E box-like motifs in the ChoRE. Based on a model structure that we generated of ChREBP/Mlx bound to the ChoRE, we hypothesized that intermolecular interactions between residues within the Mlx loop regions of adjacent heterodimers are responsible for stabilizing the complex. We tested this hypothesis by preparing Mlx variants in which the loop region was replaced with that of another family member or mutated at several key residues. These Mlx variants retained their ability to bind to a single perfect E-box motif as a heterodimer with ChREBP, but no longer bound to the ChoRE nor supported glucose responsive activity. In summary, our results support a model in which the loop regions of Mlx play an important functional role in mediating the coordinate binding of ChREBP/Mlx heterodimers to the ChoRE.

We describe PReMod, a new database of genome-wide cis-regulatory module (CRM) predictions for both the human and the mouse genomes. The prediction algorithm, described previously in Blanchette et al. (2006) Genome Res., 16, 656-668, exploits the fact that many known CRMs are made of clusters of phylogenetically conserved and repeated transcription factors (TF) binding sites. Contrary to other existing databases, PReMod is not restricted to modules located proximal to genes, but in fact mostly contains distal predicted CRMs (pCRMs). Through its web interface, PReMod allows users to (i) identify pCRMs around a gene of interest; (ii) identify pCRMs that have binding sites for a given TF (or a set of TFs) or (iii) download the entire dataset for local analyses. Queries can also be refined by filtering for specific chromosomal regions, for specific regions relative to genes or for the presence of CpG islands. The output includes information about the binding sites predicted within the selected pCRMs, and a graphical display of their distribution within the pCRMs. It also provides a visual depiction of the chromosomal context of the selected pCRMs in terms of neighboring pCRMs and genes, all of which are linked to the UCSC Genome Browser and the NCBI. PReMod: http://genomequebec.mcgill.ca/PReMod.

The NAR online Molecular Biology Database Collection is a public resource that contains links to the databases described in this issue of Nucleic Acids Research, previous NAR database issues, as well as a selection of other molecular biology databases that are freely available on the web and might be useful to the molecular biologist. The 2007 update includes 968 databases, 110 more than the previous one. Many databases that have been described in earlier issues of NAR come with updated summaries, which reflect recent progress and, in some instances, an expanded scope of these databases. The complete database list and summaries are available online on the Nucleic Acids Research web site http://nar.oxfordjournals.org/.

The variation resources within the University of California Santa Cruz Genome Browser include polymorphism data drawn from public collections and analyses of these data, along with their display in the context of other genomic annotations. Primary data from dbSNP is included for many organisms, with added information including genomic alleles and orthologous alleles for closely related organisms. Display filtering and coloring is available by variant type, functional class or other annotations. Annotation of potential errors is highlighted and a genomic alignment of the variant's flanking sequence is displayed. HapMap allele frequencies and linkage disequilibrium (LD) are available for each HapMap population, along with non-human primate alleles. The browsing and analysis tools, downloadable data files and links to documentation and other information can be found at http://genome.ucsc.edu/.

The generation of a large collection of defined transposon insertion mutants is of general interest to the Caenorhabditis elegans research community and has been supported by the European Union. We describe here a semi-automated high-throughput method for mutant production and screening, using the heterologous transposon Mos1. The procedure allows routine culture of several thousand independent nematode strains in parallel for multiple generations before stereotyped molecular analyses. Using this method, we have already generated >17 500 individual strains carrying Mos1 insertions. It could be easily adapted to forward and reverse genetic screens and may influence researchers faced with making a choice of model organism.

The Definitive Haplotype Database (D-HaploDB) is a web-accessible resource of genome-wide definitive haplotypes determined from a collection of Japanese complete hydatidiform moles (CHMs), each of which carries a genome derived from a single sperm. Currently, the database contains genotypes for 281 439 common SNPs from 74 CHMs which were determined by a high-throughput array-based oligonucleotide hybridization technique. The database also presents maps of haplotype blocks and linkage disequilibrium bins together with tagSNPs that might prove useful for association studies of disease genes. Cryptic relatedness among the samples in this study is unlikely, because the formation of a CHM is a maternal event of rare sporadic occurrence, and its genotype is that of the incoming sperm. This is demonstrated by the absence of long extended shared haplotypes (ESHs). The D-HaploDB is freely accessible via the Internet at http://orca.gen.kyushu-u.ac.jp.

5-fluorouracil (5FU) is an effective anti-cancer drug, yet its mechanism of action remains unclear. Here, we examine the effect of 5FU on pre-mRNA splicing in vivo. Using RT-PCR, we show that the splicing of a number of pre-mRNAs is inhibited in HeLa cells that have been exposed to a low dose of 5FU. It appears that this inhibitory effect is not due to its incorporation into pre-mRNA, because partially or fully 5FU-substituted pre-mRNA, when injected into Xenopus oocytes, is spliced just as well as is the unsubstituted pre-mRNA. Detailed analyses of 5FU-treated cells indicate that 5FU is incorporated into U2 snRNA at important naturally occurring pseudouridylation sites. Remarkably, 5FU incorporation effectively blocks the formation of important pseudouridines in U2 snRNA, as only a trace of pseudouridine is detected when cells are exposed to a low dose of 5FU for 5 days. Injection of the hypopseudouridylated HeLa U2 snRNA into U2-depleted Xenopus oocytes fails to reconstitute pre-mRNA splicing, whereas control U2 isolated from untreated or uracil-treated HeLa cells completely reconstitutes the splicing. Our results demonstrate for the first time that 5FU incorporates into a spliceosomal snRNA at natural pseudouridylation sites in vivo, thereby inhibiting snRNA pseudouridylation and splicing. This mechanism may contribute substantially to 5FU-mediated cell death.

The homotrimeric fusion protein nucleocapsid (NC)-dUTPase combines domains that participate in RNA/DNA folding, reverse transcription, and DNA repair in Mason-Pfizer monkey betaretrovirus infected cells. The structural organization of the fusion protein remained obscured by the N- and C-terminal flexible segments of dUTPase and the linker region connecting the two domains that are invisible in electron density maps. Small-angle X-ray scattering reveals that upon oligonucleotide binding the NC domains adopt the trimeric symmetry of dUTPase. High-resolution X-ray structures together with molecular modeling indicate that fusion with NC domains dramatically alters the conformation of the flexible C-terminus by perturbing the orientation of a critical beta-strand. Consequently, the C-terminal segment is capable of double backing upon the active site of its own monomer and stabilized by non-covalent interactions formed with the N-terminal segment. This co-folding of the dUTPase terminal segments, not observable in other homologous enzymes, is due to the presence of the fused NC domain. Structural and genomic advantages of fusing the NC domain to a shortened dUTPase in betaretroviruses and the possible physiological consequences are envisaged.

Highly accurate knockdown functional analyses based on RNA interference (RNAi) require the possible most complete hydrolysis of the targeted mRNA while avoiding the degradation of untargeted genes (off-target effects). This in turn requires significant improvements to target selection for two reasons. First, the average silencing activity of randomly selected siRNAs is as low as 62%. Second, applying more than five different siRNAs may lead to saturation of the RNA-induced silencing complex (RISC) and to the degradation of untargeted genes. Therefore, selecting a small number of highly active siRNAs is critical for maximizing knockdown and minimizing off-target effects. To satisfy these needs, a publicly available and transparent machine learning tool is presented that ranks all possible siRNAs for each targeted gene. Support vector machines (SVMs) with polynomial kernels and constrained optimization models select and utilize the most predictive effective combinations from 572 sequence, thermodynamic, accessibility and self-hairpin features over 2200 published siRNAs. This tool reaches an accuracy of 92.3% in cross-validation experiments. We fully present the underlying biophysical signature that involves free energy, accessibility and dinucleotide characteristics. We show that while complete silencing is possible at certain structured target sites, accessibility information improves the prediction of the 90% active siRNA target sites. Fast siRNA activity predictions can be performed on our web server at http://optirna.unl.edu/.

Fragile X syndrome, the most frequent form of inherited mental retardation, is due to the absence of expression of the Fragile X Mental Retardation Protein (FMRP), an RNA binding protein with high specificity for G-quartet RNA structure. FMRP is involved in several steps of mRNA metabolism: nucleocytoplasmic trafficking, translational control and transport along dendrites in neurons. Fragile X Related Protein 1 (FXR1P), a homologue and interactor of FMRP, has been postulated to have a function similar to FMRP, leading to the hypothesis that it can compensate for the absence of FMRP in Fragile X patients. Here we analyze the ability of three isoforms of FXR1P, expressed in different tissues, to bind G-quartet RNA structure specifically. Only the longest FXR1P isoform was found to be able to bind specifically the G-quartet RNA, albeit with a lower affinity as compared to FMRP, whereas the other two isoforms negatively regulate the affinity of FMRP for G-quartet RNA. This result is important to decipher the molecular basis of fragile X syndrome, through the understanding of FMRP action in the context of its multimolecular complex in different tissues. In addition, we show that the action of FXR1P is synergistic rather than compensatory for FMRP function.

We have carried out a systematic analysis of the contribution of a set of selected features that include three new features to the accuracy of operon prediction. Our analyses have led to a number of new insights about operon prediction, including that (i) different features have different levels of discerning power when used on adjacent gene pairs with different ranges of intergenic distance, (ii) certain features are universally useful for operon prediction while others are more genome-specific and (iii) the prediction reliability of operons is dependent on intergenic distances. Based on these new insights, our newly developed operon-prediction program achieves more accurate operon prediction than the previous ones, and it uses features that are most readily available from genomic sequences. Our prediction results indicate that our (non-linear) decision tree-based classifier can predict operons in a prokaryotic genome very accurately when a substantial number of operons in the genome are already known. For example, the prediction accuracy of our program can reach 90.2 and 93.7% on Bacillus subtilis and Escherichia coli genomes, respectively. When no such information is available, our (linear) logistic function-based classifier can reach the prediction accuracy at 84.6 and 83.3% for E.coli and B.subtilis, respectively.

The family of box H/ACA snoRNA is an abundant class of non-protein-coding RNAs, which play important roles in the post-transcriptional modification of rRNAs and snRNAs. Here we report the characterization in the human genome of 202 sequences derived from box H/ACA snoRNAs. Most of them were retrogenes formed using the L1 integration machinery. About 96% of the box H/ACA RNA-related sequences are found in corresponding locations on the chimpanzee and human chromosomes, while the mouse shares approximately 50% of these human sequences, suggesting that some of the H/ACA RNA-related sequences in primate occurred after the rodent/primate divergence. Of the H/ACA RNA-related sequences, 49% are found in intronic regions of protein-coding genes and 64 H/ACA-related sequences can be folded to the typical secondary structure of the box H/ACA snoRNA family, while 30 of them were recognized as functional homologs of their corresponding box H/ACA snoRNAs previously reported. Of the 64 sequences with the typical secondary structure of the box H/ACA RNA family, 11 were found in EST databases and 5 among which were shown to be expressed in more than one human tissue. Notably, U107f is nested in an intron of a protein gene coding for nudix-type motif 13, but expressed from the opposite strand, and the searching of EST databases revealed it can be expressed in liver and spleen, even in melanotic melanoma.

While methylcytosines serve as the fifth base encoding epigenetic information, they are also a dangerous endogenous mutagen due to their intrinsic instability. Methylcytosine undergoes spontaneous deamination, at a rate much higher than cytosine, to generate thymine. In mammals, two repair enzymes, thymine DNA glycosylase (TDG) and methyl-CpG binding domain 4 (MBD4), have evolved to counteract the mutagenic effect of methylcytosines. Both recognize G/T mismatches arising from methylcytosine deamination and initiate base-excision repair that corrects them to G/C pairs. However, the mechanism by which the methylation status of the repaired cytosines is restored has remained unknown. We show here that the DNA methyltransferase Dnmt3a interacts with TDG. Both the PWWP domain and the catalytic domain of Dnmt3a are able to mediate the interaction with TDG at its N-terminus. The interaction affects the enzymatic activity of both proteins: Dnmt3a positively regulates the glycosylase activity of TDG, while TDG inhibits the methylation activity of Dnmt3a in vitro. These data suggest a mechanistic link between DNA repair and remethylation at sites affected by methylcytosine deamination.

Regulation of endogenous genes by DNA-binding polyamides requires effective nuclear localization. Previous work employing confocal microscopy to study uptake of fluorophore-labeled polyamides has demonstrated the difficulty of predicting a priori the nuclear uptake of a given polyamide. The data suggest that dye identity influences uptake sufficiently such that a dye-conjugate cannot be used as a proxy for unlabeled analogs. Polyamides capable of nuclear localization unaided by fluorescent dyes are desirable due to size and other limitations of fluorophores. Recently, a polyamide-fluorescein conjugate targeted to the hypoxia response element (HRE) was found to inhibit vascular endothelial growth factor (VEGF) expression in cultured HeLa cells. The current study uses inhibition of VEGF expression as a biological read-out for effective nuclear localization of HRE-targeted polyamides. We synthesized a focused library of non-fluorescent, HRE-targeted polyamides in which the C-terminus 'tail' has been systematically varied. Members of this library bind the HRE with affinities comparable or superior to that of the fluorescein-labeled analog. Although most library members demonstrate modest or no biological activity, two non-fluorescent polyamides are reported with activity rivaling that of the previously reported fluorescein-labeled polyamide. We also show evidence that promoter occupancy by HIF-1, the transcription factor that binds the HRE, is inhibited by HRE-targeted polyamides.

The small GTPase RAB4 regulates endocytic recycling, a process that contributes to Major Histocompatibility Complex (MHC)-mediated antigen presentation by specialized antigen presenting cells (APC) of the immune system. The gene encoding the RAB4B isoform of RAB4 was singled out by two complementary genome-wide screens. One of these consisted of a computer scan to identify genes containing characteristic MHC class II-related regulatory sequences. The second was the use of chromatin immunoprecipitation coupled to microarrays (ChIP-on-chip) to identify novel targets of a transcriptional co-activator called the MHC class II transactivator (CIITA). We show that the RAB4B gene is regulated by a typical MHC class II-like enhancer that is controlled directly by both CIITA and the multiprotein transcription factor complex known as the MHC class II enhanceosome. RAB4B expression is thus activated by the same regulatory machinery that is known to be essential for the expression of MHC class II genes. This molecular link between the transcriptional activation of RAB4B and MHC class II genes implies that APC boost their antigen presentation capacity by increasing RAB4-mediated endocytic recycling.

Histone chaperones are thought to be important for maintaining the physiological activity of histones; however, their exact roles are not fully understood. The physiological function of template activating factor (TAF)-I, one of the histone chaperones, also remains unclear; however, its biochemical properties have been well studied. By performing microarray analyses, we found that TAF-I stimulates the transcription of a sub-set of genes. The transcription of endogenous genes that was up-regulated by TAF-I was found to be additively stimulated by histone acetylation. On performing an experiment with a cell line containing a model gene integrated into the chromosome, TAF-I was found to stimulate the model gene transcription in a histone chaperone activity-dependent manner additively with histone acetylation. TAF-I bound to the core histones and remodeled the chromatin structure independent of the N-terminal histone tail and its acetylation level in vitro. These results suggest that TAF-I remodel the chromatin structure through its interaction with the core domain of the histones, including the histone fold, and this mechanism is independent of the histone acetylation status.

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.