Nucleic Acids Research - Advance Access

Single nucleotide polymorphism arrays: a decade of biological, computational and technological advances

LaFramboise, T. — 2009-07-01

Array manufacturers originally designed single nucleotide polymorphism (SNP) arrays to genotype human DNA at thousands of SNPs across the genome simultaneously. In the decade since their initial development, the platform's applications have expanded to include the detection and characterization of copy number variation—whether somatic, inherited, or de novo—as well as loss-of-heterozygosity in cancer cells. The technology's impressive contributions to insights in population and molecular genetics have been fueled by advances in computational methodology, and indeed these insights and methodologies have spurred developments in the arrays themselves. This review describes the most commonly used SNP array platforms, surveys the computational methodologies used to convert the raw data into inferences at the DNA level, and details the broad range of applications. Although the long-term future of SNP arrays is unclear, cost considerations ensure their relevance for at least the next several years. Even as emerging technologies seem poised to take over for at least some applications, researchers working with these new sources of data are adopting the computational approaches originally developed for SNP arrays.

Differential requirement of a distal regulatory region for pre-initiation complex formation at globin gene promoters

2009-06-30

Although distal regulatory regions are frequent throughout the genome, the molecular mechanisms by which they act in a promoter-specific manner remain to be elucidated. The human β-globin locus constitutes an extremely well-established multigenic model to investigate this issue. In erythroid cells, the β-globin locus control region (LCR) exerts distal regulatory function by influencing local chromatin organization and inducing high-level expression of individual β-like globin genes. Moreover, in transgenic mice expressing the entire human β-globin locus, deletion of LCR-hypersensitive site 2 (HS2) can alter β-like globin gene expression. Here, we show that abnormal expression of human β-like globin genes in the absence of HS2 is associated with decreased efficacy of pre-initiation complex formation at the human - and -promoters, but not at the β-promoter. This promoter-specific phenomenon is associated with reduced long-range interactions between the HS2-deleted LCR and human -promoters. We also find that HS2 is dispensable for high-level human β-gene transcription, whereas deletion of this hypersensitive site can alter locus chromatin organization; therefore the functions exerted by HS2 in transcriptional enhancement and locus chromatin organization are distinct. Overall, our data delineate one mechanism whereby a distal regulatory region provides promoter-specific transcriptional enhancement.

Rational engineering of type II restriction endonuclease DNA binding and cleavage specificity

Morgan, R. D., Luyten, Y. A. — 2009-06-30

The type II restriction endonucleases are indispensible tools for molecular biology. Although enzymes recognizing nearly 300 unique sequences are known, the ability to engineer enzymes to recognize any sequence of choice would be valuable. However, previous attempts to engineer new recognition specificity have met limited success. Here we report the rational engineering of multiple new type II specificities. We recently identified a family of MmeI-like type II endonucleases that have highly similar protein sequences but different recognition specificity. We identified the amino-acid positions within these enzymes that determine position specific DNA base recognition at three positions within their recognition sequences through correlations between their aligned amino-acid residues and aligned recognition sequences. We then altered the amino acids at the identified positions to those correlated with recognition of a desired new base to create enzymes that recognize and cut at predictable new DNA sequences. The enzymes so altered have similar levels of endonuclease activity compared to the wild-type enzymes. Using simple and predictable mutagenesis in this family it is now possible to create hundreds of unique new type II restriction endonuclease specificities. The findings suggest a simple mechanism for the evolution of new DNA specificity in Nature.

Analysis of individual remodeled nucleosomes reveals decreased histone-DNA contacts created by hSWI/SNF

Bouazoune, K., Miranda, T. B., Jones, P. A., Kingston, R. E. — 2009-06-30

Chromatin remodeling enzymes use the energy of ATP hydrolysis to alter histone–DNA contacts and regulate DNA-based processes in eukaryotes. Whether different subfamilies of remodeling complexes generate distinct products remains uncertain. We have developed a protocol to analyze nucleosome remodeling on individual products formed in vitro. We used a DNA methyltransferase to examine DNA accessibility throughout nucleosomes that had been remodeled by the ISWI and SWI/SNF families of enzymes. We confirmed that ISWI-family enzymes mainly created patterns of accessibility consistent with canonical nucleosomes. In contrast, SWI/SNF-family enzymes generated widespread DNA accessibility. The protection patterns created by these enzymes were usually located at the extreme ends of the DNA and showed no evidence for stable loop formation on individual molecules. Instead, SWI/SNF family proteins created extensive accessibility by generating heterogeneous products that had fewer histone–DNA contacts than a canonical nucleosome, consistent with models in which a canonical histone octamer has been ‘pushed’ off of the end of the DNA.

A FRET-based assay for characterization of alternative splicing events using peptide nucleic acid fluorescence in situ hybridization

Blanco, A. M., Rausell, L., Aguado, B., Perez-Alonso, M., Artero, R. — 2009-06-26

We describe a quantitative method for detecting RNA alternative splicing variants that combines in situ hybridization of fluorescently labeled peptide nucleic acid (PNA) probes with confocal microscopy Förster resonance energy transfer (FRET). The use of PNA probes complementary to sequences flanking a given splice junction allows to specifically quantify, within the cell, the RNA isoform generating such splice junction by FRET measure. As a proof of concept we analyzed two alternative splicing events originating from lymphocyte antigen 6 (LY6) complex, locus G5B (LY6G5B) pre-mRNA. These are characterized by the removal of the first intron (Fully Spliced Isoform, FSI) or by retention of such intron (Intron-Retained Isoform, IRI). The use of PNA probe pairs labeled with donor (Cy3) and acceptor (Cy5) fluorophores, suitable to FRET, flanking FSI and IRI specific splice junctions specifically detected both mRNA isoforms in HeLa cells. We have observed that the method works efficiently with probes 5–11 nt apart. The data supports that this FRET-based PNA fluorescence in situ hybridization (FP–FISH) method offers a conceptually new approach for characterizing at the subcellular level not only splice variant isoform structure, location and dynamics but also potentially a wide variety of close range RNA–RNA interactions.

Fibroblast growth factor 1 induced during myogenesis by a transcription-translation coupling mechanism

2009-06-26

Fibroblast growth factor 1 (FGF1) is involved in muscle development and regeneration. The FGF1 gene contains four tissue-specific promoters allowing synthesis of four transcripts with distinct leader regions. Two of these transcripts contain internal ribosome entry sites (IRESs), which are RNA elements allowing mRNA translation to occur in conditions of blockade of the classical cap-dependent mechanism. Here, we investigated the function and the regulation of FGF1 during muscle differentiation and regeneration. Our data show that FGF1 protein expression is induced in differentiating myoblasts and regenerating mouse muscle, whereas siRNA knock-down demonstrated FGF1 requirement for myoblast differentiation. FGF1 induction occurred at both transcriptional and translational levels, involving specific activation of both promoter A and IRES A, whereas global cap-dependent translation was inhibited. Furthermore, we identified, in the FGF1 promoter A distal region, a cis-acting element able to activate the IRES A-driven translation. These data revealed a mechanism of molecular coupling of mRNA transcription and translation, involving a unique process of IRES activation by a promoter element. The crucial role of FGF1 in myoblast differentiation provides physiological relevance to this novel mechanism. This finding also provides a new insight into the molecular mechanisms linking different levels of gene expression regulation.

Discovering ligands for a microRNA precursor with peptoid microarrays

Chirayil, S., Chirayil, R., Luebke, K. J. — 2009-06-26

We have screened peptoid microarrays to identify specific ligands for the RNA hairpin precursor of miR-21, a microRNA involved in cancer and heart disease. Microarrays were printed by spotting a library of 7680 N-substituted oligoglycines (peptoids) onto glass slides. Two compounds on the array specifically bind RNA having the sequence and predicted secondary structure of the miR-21 precursor hairpin and have specific affinity for the target in solution. Their binding induces a conformational change around the hairpin loop, and the most specific compound recognizes the loop sequence and a bulged uridine in the proximal duplex. Functional groups contributing affinity and specificity were identified, and by varying a critical methylpyridine group, a compound with a dissociation constant of 1.9 µM for the miR-21 precursor hairpin and a 20-fold discrimination against a closely-related hairpin was created. This work describes a systematic approach to discovery of ligands for specific pre-defined novel RNA structures. It demonstrates discovery of new ligands for an RNA for which no specific lead compounds were previously known by screening a microarray of small molecules.

JunB mediates enhancer/promoter activity of COL1A2 following TGF-{beta} induction

Ponticos, M., Harvey, C., Ikeda, T., Abraham, D., Bou-Gharios, G. — 2009-06-26

Transcriptional control of the genes coding for collagen type I is regulated by a complex interaction between a distal enhancer and a proximal promoter. In this study, we have dissected the molecular mechanism of this interaction by defining a specific sequence within the enhancer that respond in fibroblasts to transforming growth factor-β (TGF-β). We show that TGF-β activates COL1A2 gene via a non-canonical (Smad-independent) signalling pathway, which requires enhancer/promoter co-operation. This interaction involves exchange of cJun/Jun B transcription factor occupancy of a critical enhancer site resulting in the stabilization of enhancer/promoter coalescence. Moreover, using transgenesis, we show that interference in this mechanism results in the abolition of COL1A2 fibroblast expression in vivo. These data are therefore relevant to the control of collagen type I in vivo both in embryonic development, in adult connective tissue homeostasis, and in tissue repair and scarring pathologies.

Structural characterization of cationic lipid-tRNA complexes

2009-06-26

Despite considerable interest and investigations on cationic lipid–DNA complexes, reports on lipid–RNA interaction are very limited. In contrast to lipid–DNA complexes where lipid binding induces partial B to A and B to C conformational changes, lipid–tRNA complexation preserves tRNA folded state. This study is the first attempt to investigate the binding of cationic lipid with transfer RNA and the effect of lipid complexation on tRNA aggregation and condensation. We examine the interaction of tRNA with cholesterol (Chol), 1,2-dioleoyl-3-trimethylammonium-propane (DOTAP), dioctadecyldimethylammoniumbromide (DDAB) and dioleoylphosphatidylethanolamine (DOPE), at physiological condition, using constant tRNA concentration and various lipid contents. FTIR, UV-visible, CD spectroscopic methods and atomic force microscopy (AFM) were used to analyze lipid binding site, the binding constant and the effects of lipid interaction on tRNA stability, conformation and condensation. Structural analysis showed lipid–tRNA interactions with G–C and A–U base pairs as well as the backbone phosphate group with overall binding constants of K_Chol = 5.94 (± 0.8) x 10⁴ M^–1, K_DDAB = 8.33 (± 0.90) x 10⁵ M^–1, K_DOTAP = 1.05 (± 0.30) x 10⁵ M^–1 and K_DOPE = 2.75 (± 0.50) x 10⁴ M^–1. The order of stability of lipid–tRNA complexation is DDAB > DOTAP > Chol > DOPE. Hydrophobic interactions between lipid aliphatic tails and tRNA were observed. RNA remains in A-family structure, while biopolymer aggregation and condensation occurred at high lipid concentrations.

Conserved principles of mammalian transcriptional regulation revealed by RNA half-life

Friedel, C. C., Dolken, L., Ruzsics, Z., H. Koszinowski, U., Zimmer, R. — 2009-06-26

RNA levels in a cell are regulated by the relative rates of RNA synthesis and decay. We recently developed a new approach for measuring both RNA synthesis and decay in a single experimental setting by biosynthetic labeling of newly transcribed RNA. Here, we show that this provides measurements of RNA half-lives from microarray data with a so far unreached accuracy. Based on such measurements of RNA half-lives for human B-cells and mouse fibroblasts, we identified conserved regulatory principles for a large number of biological processes. We show that different regulatory patterns between functionally similar proteins are characterized by differences in the half-life of the corresponding transcripts and can be identified by measuring RNA half-life. We identify more than 100 protein families which show such differential regulatory patterns in both species. Additionally, we provide strong evidence that the activity of protein complexes consisting of subunits with overall long transcript half-lives can be regulated by transcriptional regulation of individual key subunits with short-lived transcripts. Based on this observation, we predict more than 100 key regulatory subunits for human complexes of which 28% could be confirmed in mice (P < 10^–9). Therefore, this atlas of transcript half-lives provides new fundamental insights into many cellular processes.

Senescence delay and repression of p16INK4a by Lsh via recruitment of histone deacetylases in human diploid fibroblasts

Zhou, R., Han, L., Li, G., Tong, T. — 2009-06-26

Lymphoid specific helicase (Lsh) belongs to the family of SNF2/helicases. Disruption of Lsh leads to developmental growth retardation and premature aging in mice. However, the specific effect of Lsh on human cellular senescence remains unknown. Herein, we report that Lsh overexpression delays cell senescence by silencing p16^INK4a in human fibroblasts. The patterns of p16^INK4a and Lsh expression during cell senescence present the inverse correlation. We also find that Lsh requires histone deacetylase (HDAC) activity to repress p16^INK4a and treatment with trichostatin A (TSA) is sufficient to block the repressor effect of Lsh. Moreover, overexpression of Lsh is correlated with deacetylation of histone H3 at the p16 promoter, and TSA treatment in Lsh-expressing cells reverses the acetylation status of histones. Additionally, we demonstrate an interaction between Lsh, histone deacetylase 1 (HDAC1) and HDAC2 in vivo. Furthermore, we demonstrate that Lsh interacts in vivo with the p16 promoter and recruits HDAC1. Our data suggest that Lsh represses endogenous p16^INK4a expression by recruiting HDAC to establish a repressive chromatin structure at the p16^INK4a promoter, which in turn delays cell senescence.

Position of eukaryotic translation initiation factor eIF1A on the 40S ribosomal subunit mapped by directed hydroxyl radical probing

2009-06-26

The universally conserved eukaryotic initiation factor (eIF), eIF1A, plays multiple roles throughout initiation: it stimulates eIF2/GTP/Met-tRNA_i^Met attachment to 40S ribosomal subunits, scanning, start codon selection and subunit joining. Its bacterial ortholog IF1 consists of an oligonucleotide/oligosaccharide-binding (OB) domain, whereas eIF1A additionally contains a helical subdomain, N-terminal tail (NTT) and C-terminal tail (CTT). The NTT and CTT both enhance ribosomal recruitment of eIF2/GTP/Met-tRNA_i^Met, but have opposite effects on the stringency of start codon selection: the CTT increases, whereas the NTT decreases it. Here, we determined the position of eIF1A on the 40S subunit by directed hydroxyl radical cleavage. eIF1A's OB domain binds in the A site, similar to IF1, whereas the helical subdomain contacts the head, forming a bridge over the mRNA channel. The NTT and CTT both thread under Met-tRNA_i^Met reaching into the P-site. The NTT threads closer to the mRNA channel. In the proposed model, the NTT does not clash with either mRNA or Met-tRNA_i^Met, consistent with its suggested role in promoting the ‘closed’ conformation of ribosomal complexes upon start codon recognition. In contrast, eIF1A-CTT appears to interfere with the P-site tRNA-head interaction in the ‘closed’ complex and is likely ejected from the P-site upon start codon recognition.

Structural basis for the sequence-specific RNA-recognition mechanism of human CUG-BP1 RRM3

2009-06-24

The CUG-binding protein 1 (CUG-BP1) is a member of the CUG-BP1 and ETR-like factors (CELF) family or the Bruno-like family and is involved in the control of splicing, translation and mRNA degradation. Several target RNA sequences of CUG-BP1 have been predicted, such as the CUG triplet repeat, the GU-rich sequences and the AU-rich element of nuclear pre-mRNAs and/or cytoplasmic mRNA. CUG-BP1 has three RNA-recognition motifs (RRMs), among which the third RRM (RRM3) can bind to the target RNAs on its own. In this study, we solved the solution structure of the CUG-BP1 RRM3 by hetero-nuclear NMR spectroscopy. The CUG-BP1 RRM3 exhibited a noncanonical RRM fold, with the four-stranded β-sheet surface tightly associated with the N-terminal extension. Furthermore, we determined the solution structure of the CUG-BP1 RRM3 in the complex with (UG)₃ RNA, and discovered that the UGU trinucleotide is specifically recognized through extensive stacking interactions and hydrogen bonds within the pocket formed by the β-sheet surface and the N-terminal extension. This study revealed the unique mechanism that enables the CUG-BP1 RRM3 to discriminate the short RNA segment from other sequences, thus providing the molecular basis for the comprehension of the role of the RRM3s in the CELF/Bruno-like family.

Comparative analysis of activator-E{sigma}54 complexes formed with nucleotide-metal fluoride analogues

Burrows, P. C., Joly, N., Nixon, B. T., Buck, M. — 2009-06-24

Bacterial RNA polymerase (RNAP) containing the major variant ⁵⁴ factor forms open promoter complexes in a reaction in which specialized activator proteins hydrolyse ATP. Here we probe binding interactions between ⁵⁴-RNAP (E⁵⁴) and the ATPases associated with various cellular activities (AAA+) domain of the Escherichia coli activator protein, PspF, using nucleotide-metal fluoride (BeF and AlF) analogues representing ground and transition states of ATP, which allow complexes (that are otherwise too transient with ATP) to be captured. We show that the organization and functionality of the ADP–BeF- and ADP–AlF-dependent complexes greatly overlap. Our data support an activation pathway in which the initial ATP-dependent binding of the activator to the E⁵⁴ closed complex results in the re-organization of E⁵⁴ with respect to the transcription start-site. However, the nucleotide-dependent binding interactions between the activator and the E⁵⁴ closed complex are in themselves insufficient for forming open promoter complexes when linear double-stranded DNA is present in the initial closed complex.

Preparation of selective and segmentally labeled single-stranded DNA for NMR by self-primed PCR and asymmetrical endonuclease double digestion

Nelissen, F. H. T., Girard, F. C., Tessari, M., Heus, H. A., Wijmenga, S. S. — 2009-06-24

We demonstrate a new, efficient and easy-to-use method for enzymatic synthesis of (stereo-)specific and segmental ¹³C/¹⁵N/²H isotope-labeled single-stranded DNA in amounts sufficient for NMR, based on the highly efficient self-primed PCR. To achieve this, new approaches are introduced and combined. (i) Asymmetric endonuclease double digestion of tandem-repeated PCR product. (ii) T4 DNA ligase mediated ligation of two ssDNA segments. (iii) In vitro dNTP synthesis, consisting of in vitro rNTP synthesis followed by enzymatic stereo-selective reduction of the C2' of the rNTP, and a one-pot add-up synthesis of dTTP from dUTP. The method is demonstrated on two ssDNAs: (i) a 36-nt three-way junction, selectively ¹³C₉/¹⁵N₃/²H_{(1',2'',3',4',5',5'')}-dC labeled and (ii) a 39-nt triple-repeat three-way junction, selectively ¹³C₉/¹⁵N₃/²H_{(1',2'',3',4',5',5'')}-dC and ¹³C₉/¹⁵N₂/²H_{(1',2'',3',4',5',5'')}-dT labeled in segment C20-C39. Their NMR spectra show the spectral simplification, while the stereo-selective ²H-labeling in the deoxyribose of the dC-residues, straightforwardly provided assignment of their C1'–H2' and C2'–H2' resonances. The labeling protocols can be extended to larger ssDNA molecules and to more than two segments.

Extracting transcription factor targets from ChIP-Seq data

Tuteja, G., White, P., Schug, J., Kaestner, K. H. — 2009-06-24

ChIP-Seq technology, which combines chromatin immunoprecipitation (ChIP) with massively parallel sequencing, is rapidly replacing ChIP-on-chip for the genome-wide identification of transcription factor binding events. Identifying bound regions from the large number of sequence tags produced by ChIP-Seq is a challenging task. Here, we present GLITR (GLobal Identifier of Target Regions), which accurately identifies enriched regions in target data by calculating a fold-change based on random samples of control (input chromatin) data. GLITR uses a classification method to identify regions in ChIP data that have a peak height and fold-change which do not resemble regions in an input sample. We compare GLITR to several recent methods and show that GLITR has improved sensitivity for identifying bound regions closely matching the consensus sequence of a given transcription factor, and can detect bona fide transcription factor targets missed by other programs. We also use GLITR to address the issue of sequencing depth, and show that sequencing biological replicates identifies far more binding regions than re-sequencing the same sample.

Interplay of DNA supercoiling and catenation during the segregation of sister duplexes

2009-06-24

The discrete regulation of supercoiling, catenation and knotting by DNA topoisomerases is well documented both in vivo and in vitro, but the interplay between them is still poorly understood. Here we studied DNA catenanes of bacterial plasmids arising as a result of DNA replication in Escherichia coli cells whose topoisomerase IV activity was inhibited. We combined high-resolution two-dimensional agarose gel electrophoresis with numerical simulations in order to better understand the relationship between the negative supercoiling of DNA generated by DNA gyrase and the DNA interlinking resulting from replication of circular DNA molecules. We showed that in those replication intermediates formed in vivo, catenation and negative supercoiling compete with each other. In interlinked molecules with high catenation numbers negative supercoiling is greatly limited. However, when interlinking decreases, as required for the segregation of newly replicated sister duplexes, their negative supercoiling increases. This observation indicates that negative supercoiling plays an active role during progressive decatenation of newly replicated DNA molecules in vivo.

Degradation of nanoRNA is performed by multiple redundant RNases in Bacillus subtilis

Fang, M., Zeisberg, W.-M., Condon, C., Ogryzko, V., Danchin, A., Mechold, U. — 2009-06-24

Escherichia coli possesses only one essential oligoribonuclease (Orn), an enzyme that can degrade oligoribonucleotides of five residues and shorter in length (nanoRNA). Firmicutes including Bacillus subtilis do not have an Orn homolog. We had previously identified YtqI (NrnA) as functional analog of Orn in B. subtilis. Screening a genomic library from B. subtilis for genes that can complement a conditional orn mutant, we identify here YngD (NrnB) as a second nanoRNase in B. subtilis. Like NrnA, NrnB is a member of the DHH/DHHA1 protein family of phosphoesterases. NrnB degrades nanoRNA 5-mers in vitro similarily to Orn. Low expression levels of NrnB are sufficient for orn complementation. YhaM, a known RNase present in B. subtilis, degrades nanoRNA efficiently in vitro but requires high levels of expression for only partial complementation of the orn^– strain. A triple mutant (nrnA^–, nrnB^–, yhaM^–) in B. subtilis is viable and shows almost no impairment in growth. Lastly, RNase J1 seems also to have some 5'-to-3' exoribonuclease activity on nanoRNA and thus can potentially finish degradation of RNA. We conclude that, unlike in E. coli, degradation of nanoRNA is performed in a redundant fashion in B. subtilis.

Separation of 1-23-kb complementary DNA strands by urea-agarose gel electrophoresis

Hegedus, E., Kokai, E., Kotlyar, A., Dombradi, V., Szabo, G. — 2009-06-23

Double-stranded (ds), as well as denatured, single-stranded (ss) DNA samples can be analyzed on urea–agarose gels. Here we report that after denaturation by heat in the presence of 8 M urea, the two strands of the same ds DNA fragment of ~1–20-kb size migrate differently in 1 M urea containing agarose gels. The two strands are readily distinguished on Southern blots by ss-specific probes. The different migration of the two strands could be attributed to their different, base composition-dependent conformation impinging on the electrophoretic mobility of the ss molecules. This phenomenon can be exploited for the efficient preparation of strand-specific probes and for the separation of the complementary DNA strands for subsequent analysis, offering a new tool for various cell biological research areas.

BRCA2-dependent homologous recombination is required for repair of Arsenite-induced replication lesions in mammalian cells

Ying, S., Myers, K., Bottomley, S., Helleday, T., Bryant, H. E. — 2009-06-23

Arsenic exposure constitutes one of the most widespread environmental carcinogens, and is associated with increased risk of many different types of cancers. Here we report that arsenite (As[III]) can induce both replication-dependent DNA double-strand breaks (DSB) and homologous recombination (HR) at doses as low as 5 µM (0.65 mg/l), which are within the typical doses often found in drinking water in contaminated areas. We show that the production of DSBs is dependent on active replication and is likely to be the result of conversion of a DNA single-strand break (SSB) into a toxic DSB when encountered by a replication fork. We demonstrate that HR is required for the repair of these breaks and show that a functional HR pathway protects against As[III]-induced cytotoxicity. In addition, BRCA2-deficient cells are sensitive to As[III] and we suggest that As[III] could be exploited as a therapy for HR-deficient tumours such as BRCA1 and BRCA2 mutated breast and ovarian cancers.

A trans-splicing group I intron and tRNA-hyperediting in the mitochondrial genome of the lycophyte Isoetes engelmannii

Grewe, F., Viehoever, P., Weisshaar, B., Knoop, V. — 2009-06-23

Plant mitochondrial genomes show much more evolutionary plasticity than those of animals. We analysed the first mitochondrial DNA (mtDNA) of a lycophyte, the quillwort Isoetes engelmannii, which is separated from seed plants by more than 350 million years of evolution. The Isoetes mtDNA is particularly rich in recombination events, and chloroplast as well as nuclear DNA inserts document the incorporation of foreign sequences already in this most ancestral vascular plant lineage. On the other hand, particularly small group II introns and short intergenic regions reveal a tendency of evolution towards a compact mitochondrial genome. RNA editing reaches extreme levels exceeding 100 pyrimidine exchanges in individual mRNAs and, hitherto unobserved in such frequency, also in tRNAs with 18 C-to-U conversions in the tRNA for proline. In total, some 1500 sites of RNA editing can be expected for the Isoetes mitochondrial transcriptome. As a unique molecular novelty, the Isoetes cox1 gene requires trans-splicing via a discontinuous group I intron demonstrating disrupted, but functional, RNAs for yet another class of natural ribozymes.

Analysis of repair mechanism choice during homologous recombination

Agmon, N., Pur, S., Liefshitz, B., Kupiec, M. — 2009-06-23

Double-strand breaks (DSBs) occur frequently during cell growth. Due to the presence of repeated sequences in the genome, repair of a single DSB can result in gene conversion, translocation, deletion or tandem duplication depending on the mechanism and the sequence chosen as partner for the recombinational repair. Here, we study how yeast cells repair a single, inducible DSB when there are several potential donors to choose from, in the same chromosome and elsewhere in the genome. We systematically investigate the parameters that affect the choice of mechanism, as well as its genetic regulation. Our results indicate that intrachromosomal homologous sequences are always preferred as donors for repair. We demonstrate the occurrence of a novel tri-partite repair product that combines ectopic gene conversion and deletion. In addition, we show that increasing the distance between two repeated sequences enhances the dependence on Rad51 for colony formation after DSB repair. This is due to a role of Rad51 in the recovery from the checkpoint signal induced by the DSB. We suggest a model for the competition between the different homologous recombination pathways. Our model explains how different repair mechanisms are able to compensate for each other during DSB repair.

A thermodynamic approach to PCR primer design

Mann, T., Humbert, R., Dorschner, M., Stamatoyannopoulos, J., Noble, W. S. — 2009-06-23

We developed a primer design method, Pythia, in which state of the art DNA binding affinity computations are directly integrated into the primer design process. We use chemical reaction equilibrium analysis to integrate multiple binding energy calculations into a conservative measure of polymerase chain reaction (PCR) efficiency, and a precomputed index on genomic sequences to evaluate primer specificity. We show that Pythia can design primers with success rates comparable with those of current methods, but yields much higher coverage in difficult genomic regions. For example, in RepeatMasked sequences in the human genome, Pythia achieved a median coverage of 89% as compared with a median coverage of 51% for Primer3. For parameter settings yielding sensitivities of 81%, our method has a recall of 97%, compared with the Primer3 recall of 48%. Because our primer design approach is based on the chemistry of DNA interactions, it has fewer and more physically meaningful parameters than current methods, and is therefore easier to adjust to specific experimental requirements. Our software is freely available at http://pythia.sourceforge.net.

Stochastic noise in splicing machinery

Melamud, E., Moult, J. — 2009-06-22

The number of known alternative human isoforms has been increasing steadily with the amount of available transcription data. To date, over 100 000 isoforms have been detected in EST libraries, and at least 75% of human genes have at least one alternative isoform. In this paper, we propose that most alternative splicing events are the result of noise in the splicing process. We show that the number of isoforms and their abundance can be predicted by a simple stochastic noise model that takes into account two factors: the number of introns in a gene and the expression level of a gene. The results strongly support the hypothesis that most alternative splicing is a consequence of stochastic noise in the splicing machinery, and has no functional significance. The results are also consistent with error rates tuned to ensure that an adequate level of functional product is produced and to reduce the toxic effect of accumulation of misfolding proteins. Based on simulation of sampling of virtual cDNA libraries, we estimate that error rates range from 1 to 10% depending on the number of introns and the expression level of a gene.

Non-coding murine centromeric transcripts associate with and potentiate Aurora B kinase

Ferri, F., Bouzinba-Segard, H., Velasco, G., Hube, F., Francastel, C. — 2009-06-19

Non-coding RNAs are emerging as key players in many fundamental biological processes, including specification of higher-order chromatin structure. We examined the implication of RNA transcribed from mouse centromeric minor satellite repeats in the formation and function of centromere-associated complexes. Here we show that the levels of minor satellite RNA vary during cell-cycle progression, peaking in G2/M phase, concomitant with accumulation of proteins of the chromosomal passenger complex near the centromere. Consistent with this, we describe that murine minor satellite RNA are components of CENP-A-associated centromeric fractions and associate with proteins of the chromosomal passenger complex Aurora B and Survivin at the onset of mitosis. Interactions of endogenous Aurora B with CENP-A and Survivin are sensitive to RNaseA. Likewise, the kinase activity of Aurora B requires an RNA component. More importantly, Aurora B kinase activity can be potentiated by minor satellite RNA. In addition, decreased Aurora B activity after RNA depletion can be specifically rescued by restitution of these transcripts. Together, our data provide new functional evidence for minor satellite transcripts as key partners and regulators of the mitotic kinase Aurora B.

A genome-wide library of CB4856/N2 introgression lines of Caenorhabditis elegans

Doroszuk, A., Snoek, L. B., Fradin, E., Riksen, J., Kammenga, J. — 2009-06-19

Recombinant inbred lines (RILs) derived from Caenorhabditis elegans wild-type N2 and CB4856 are increasingly being used for mapping genes underlying complex traits. To speed up mapping and gene discovery, introgression lines (ILs) offer a powerful tool for more efficient QTL identification. We constructed a library of 90 ILs, each carrying a single homozygous CB4856 genomic segment introgressed into the genetic background of N2. The ILs were genotyped by 123 single-nucleotide polymorphism (SNP) markers. The proportion of the CB4856 segments in most lines does not exceed 3%, and together the introgressions cover 96% of the CB4856 genome. The value of the IL library was demonstrated by identifying novel loci underlying natural variation in two ageing-related traits, i.e. lifespan and pharyngeal pumping rate. Bin mapping of lifespan resulted in six QTLs, which all have a lifespan-shortening effect on the CB4856 allele. We found five QTLs for the decrease in pumping rate, of which four colocated with QTLs found for average lifespan. This suggests pleiotropic or closely linked QTL associated with lifespan and pumping rate. Overall, the presented IL library provides a versatile resource toward easier and efficient fine mapping and functional analyses of loci and genes underlying complex traits in C. elegans.

A computational procedure to identify significant overlap of differentially expressed and genomic imbalanced regions in cancer datasets

2009-06-19

The integration of high-throughput genomic data represents an opportunity for deciphering the interplay between structural and functional organization of genomes and for discovering novel biomarkers. However, the development of integrative approaches to complement gene expression (GE) data with other types of gene information, such as copy number (CN) and chromosomal localization, still represents a computational challenge in the genomic arena. This work presents a computational procedure that directly integrates CN and GE profiles at genome-wide level. When applied to DNA/RNA paired data, this approach leads to the identification of Significant Overlaps of Differentially Expressed and Genomic Imbalanced Regions (SODEGIR). This goal is accomplished in three steps. The first step extends to CN a method for detecting regional imbalances in GE. The second part provides the integration of CN and GE data and identifies chromosomal regions with concordantly altered genomic and transcriptional status in a tumor sample. The last step elevates the single-sample analysis to an entire dataset of tumor specimens. When applied to study chromosomal aberrations in a collection of astrocytoma and renal carcinoma samples, the procedure proved to be effective in identifying discrete chromosomal regions of coordinated CN alterations and changes in transcriptional levels.

Identification of small gains and losses in single cells after whole genome amplification on tiling oligo arrays

2009-06-18

Clinical DNA is often available in limited quantities requiring whole-genome amplification for subsequent genome-wide assessment of copy-number variation (CNV) by array-CGH. In pre-implantation diagnosis and analysis of micrometastases, even merely single cells are available for analysis. However, procedures allowing high-resolution analyses of CNVs from single cells well below resolution limits of conventional cytogenetics are lacking. Here, we applied amplification products of single cells and of cell pools (5 or 10 cells) from patients with developmental delay, cancer cell lines and polar bodies to various oligo tiling array platforms with a median probe spacing as high as 65 bp. Our high-resolution analyses reveal that the low amounts of template DNA do not result in a completely unbiased whole genome amplification but that stochastic amplification artifacts, which become more obvious on array platforms with tiling path resolution, cause significant noise. We implemented a new evaluation algorithm specifically for the identification of small gains and losses in such very noisy ratio profiles. Our data suggest that when assessed with sufficiently sensitive methods high-resolution oligo-arrays allow a reliable identification of CNVs as small as 500 kb in cell pools (5 or 10 cells), and of 2.6–3.0 Mb in single cells.

Activation of interferon regulatory factor-3 via toll-like receptor 3 and immunomodulatory functions detected in A549 lung epithelial cells exposed to misplaced U1-snRNA

Sadik, C. D., Bachmann, M., Pfeilschifter, J., Muhl, H. — 2009-06-18

U1-snRNA is an integral part of the U1 ribonucleoprotein pivotal for pre-mRNA splicing. Toll-like receptor (TLR) signaling has recently been associated with immunoregulatory capacities of U1-snRNA. Using lung A549 epithelial/carcinoma cells, we report for the first time on interferon regulatory factor (IRF)-3 activation initiated by endosomally delivered U1-snRNA. This was associated with expression of the IRF3-inducible genes interferon-β (IFN-β), CXCL10/IP-10 and indoleamine 2,3-dioxygenase. Mutational analysis of the U1-snRNA-activated IFN-β promoter confirmed the crucial role of the PRDIII element, previously proven pivotal for promoter activation by IRF3. Notably, expression of these parameters was suppressed by bafilomycin A₁, an inhibitor of endosomal acidification, implicating endosomal TLR activation. Since resiquimod, an agonist of TLR7/8, failed to stimulate A549 cells, data suggest TLR3 to be of prime relevance for cellular activation. To assess the overall regulatory potential of U1-snRNA-activated epithelial cells on cytokine production, co-cultivation with peripheral blood mononuclear cells (PBMC) was performed. Interestingly, A549 cells activated by U1-snRNA reinforced phytohemagglutinin-induced interleukin-10 release by PBMC but suppressed that of tumor necrosis factor-, indicating an anti-inflammatory potential of U1-snRNA. Since U1-snRNA is enriched in apoptotic bodies and epithelial cells are capable of performing efferocytosis, the present data in particular connect to immunobiological aspects of apoptosis at host/environment interfaces.

Establishment of Histone Modifications after Chromatin Assembly

Scharf, A. N. D., Barth, T. K., Imhof, A. — 2009-06-18

Every cell has to duplicate its entire genome during S-phase of the cell cycle. After replication, the newly synthesized DNA is rapidly assembled into chromatin. The newly assembled chromatin ‘matures’ and adopts a variety of different conformations. This differential packaging of DNA plays an important role for the maintenance of gene expression patterns and has to be reliably copied in each cell division. Posttranslational histone modifications are prime candidates for the regulation of the chromatin structure. In order to understand the maintenance of chromatin structures, it is crucial to understand the replication of histone modification patterns. To study the kinetics of histone modifications in vivo, we have pulse-labeled synchronized cells with an isotopically labeled arginine (¹⁵N₄) that is 4 Da heavier than the naturally occurring ¹⁴N₄ isoform. As most of the histone synthesis is coupled with replication, the cells were arrested at the G1/S boundary, released into S-phase and simultaneously incubated in the medium containing heavy arginine, thus labeling all newly synthesized proteins. This method allows a comparison of modification patterns on parental versus newly deposited histones. Experiments using various pulse/chase times show that particular modifications have considerably different kinetics until they have acquired a modification pattern indistinguishable from the parental histones.

The Myb/SANT domain of the telomere-binding protein TRF2 alters chromatin structure

Baker, A. M., Fu, Q., Hayward, W., Lindsay, S. M., Fletcher, T. M. — 2009-06-16

Eukaryotic DNA is packaged into chromatin, which regulates genome activities such as telomere maintenance. This study focuses on the interactions of a myb/SANT DNA-binding domain from the telomere-binding protein, TRF2, with reconstituted telomeric nucleosomal array fibers. Biophysical characteristics of the factor-bound nucleosomal arrays were determined by analytical agarose gel electrophoresis (AAGE) and single molecules were visualized by atomic force microscopy (AFM). The TRF2 DNA-binding domain (TRF2 DBD) neutralized more negative charge on the surface of nucleosomal arrays than histone-free DNA. Binding of TRF2 DBD at lower concentrations increased the radius and conformational flexibility, suggesting a distortion of the fiber structure. Additional loading of TRF2 DBD onto the nucleosomal arrays reduced the flexibility and strongly blocked access of micrococcal nuclease as contour lengths shortened, consistent with formation of a unique, more compact higher-order structure. Mirroring the structural results, TRF2 DBD stimulated a strand invasion-like reaction, associated with telomeric t-loops, at lower concentrations while inhibiting the reaction at higher concentrations. Full-length TRF2 was even more effective at stimulating this reaction. The TRF2 DBD had less effect on histone-free DNA structure and did not stimulate the t-loop reaction with this substrate, highlighting the influence of chromatin structure on the activities of DNA-binding proteins.

Bidirectional silencing of RNA polymerase I transcription by a strand switch region in Trypanosoma brucei

Haenni, S., Studer, E., Burkard, G. S., Roditi, I. — 2009-06-16

The procyclin genes in Trypanosoma brucei are transcribed by RNA polymerase I as part of 5–10 kb long polycistronic transcription units on chromosomes VI and X. Each procyclin locus begins with two procyclin genes followed by at least one procyclin-associated gene (PAG). In procyclic (insect midgut) form trypanosomes, PAG mRNA levels are about 100-fold lower than those of procyclins. We show that deletion of PAG1, PAG2 or PAG3 results in increased mRNA levels from downstream genes in the same transcription unit. Nascent RNA analysis revealed that most of the effects are due to increased transcription elongation in the knockouts. Furthermore, transient and stable transfections showed that sequence elements on both strands of PAG1 can inhibit Pol I transcription. Finally, by database mining we identified 30 additional PAG-related sequences that are located almost exclusively at strand switch regions and/or at sites where a change of RNA polymerase type is likely to occur.

Homologs of the small RNA SgrS are broadly distributed in enteric bacteria but have diverged in size and sequence

Horler, R. S. P., Vanderpool, C. K. — 2009-06-16

Sugar phosphate stress in Escherichia coli is sensed and managed by the transcriptional regulator SgrR and the small RNA (sRNA) SgrS. SgrS is a dual function RNA that performs base pairing-dependent regulation of mRNA targets and encodes a small protein, SgrT. Homologs of SgrR were analyzed for gene synteny and inter-homolog identity to identify those that are likely to be functionally analogous. These 22 SgrR homologs were used to manually locate adjacent sRNAs functionally analogous to SgrS. SgrS homologs shared little sequence identity with E. coli SgrS, but most shared several structural features. The most conserved feature of SgrS homologs was the base pairing region while the most variable feature was the sgrT-coding sequence. Analyses of predicted interactions between SgrS:ptsG mRNA pairs in different organisms revealed interesting differences in the patterns of base pairing interactions. RNA pairs with more interrupted regions of complementarity had a higher proportion of G:C base pairs than those with longer contiguous stretches of complementarity. The identification of this set of homologous sRNAs and their targets sets the stage for future studies to further elucidate the molecular requirements for regulation by SgrS.

ORMA: a tool for identification of species-specific variations in 16S rRNA gene and oligonucleotides design

2009-06-16

16S rRNA gene is one of the preferred targets for resolving species phylogenesis issues in microbiological-related contexts. However, the identification of single-nucleotide variations capable of distinguishing a sequence among a set of homologous ones can be problematic. Here we present ORMA (Oligonucleotide Retrieving for Molecular Applications), a set of scripts for discriminating positions search and for performing the selection of high-quality oligonucleotide probes to be used in molecular applications. Two assays based on Ligase Detection Reaction (LDR) are presented. First, a new set of probe pairs on cyanobacteria 16S rRNA sequences of 18 different species was compared to that of a previous study. Then, a set of LDR probe pairs for the discrimination of 13 pathogens contaminating bovine milk was evaluated. The software determined more than 100 candidate probe pairs per dataset, from more than 300 16S rRNA sequences, in less than 5 min. Results demonstrated how ORMA improved the performance of the LDR assay on cyanobacteria, correctly identifying 12 out of 14 samples, and allowed the perfect discrimination among the 13 milk pathogenic-related species. ORMA represents a significant improvement from other contexts where enzyme-based techniques have been employed on already known mutations of a single base or on entire subsequences.

Using reads to annotate the genome: influence of length, background distribution, and sequence errors on prediction capacity

Philippe, N., Boureux, A., Brehelin, L., Tarhio, J., Commes, T., Rivals, E. — 2009-06-16

Ultra high-throughput sequencing is used to analyse the transcriptome or interactome at unprecedented depth on a genome-wide scale. These techniques yield short sequence reads that are then mapped on a genome sequence to predict putatively transcribed or protein-interacting regions. We argue that factors such as background distribution, sequence errors, and read length impact on the prediction capacity of sequence census experiments. Here we suggest a computational approach to measure these factors and analyse their influence on both transcriptomic and epigenomic assays. This investigation provides new clues on both methodological and biological issues. For instance, by analysing chromatin immunoprecipitation read sets, we estimate that 4.6% of reads are affected by SNPs. We show that, although the nucleotide error probability is low, it significantly increases with the position in the sequence. Choosing a read length above 19 bp practically eliminates the risk of finding irrelevant positions, while above 20 bp the number of uniquely mapped reads decreases. With our procedure, we obtain 0.6% false positives among genomic locations. Hence, even rare signatures should identify biologically relevant regions, if they are mapped on the genome. This indicates that digital transcriptomics may help to characterize the wealth of yet undiscovered, low-abundance transcripts.

Rapid prototyping of 3D DNA-origami shapes with caDNAno

2009-06-16

DNA nanotechnology exploits the programmable specificity afforded by base-pairing to produce self-assembling macromolecular objects of custom shape. For building megadalton-scale DNA nanostructures, a long ‘scaffold’ strand can be employed to template the assembly of hundreds of oligonucleotide ‘staple’ strands into a planar antiparallel array of cross-linked helices. We recently adapted this ‘scaffolded DNA origami’ method to producing 3D shapes formed as pleated layers of double helices constrained to a honeycomb lattice. However, completing the required design steps can be cumbersome and time-consuming. Here we present caDNAno, an open-source software package with a graphical user interface that aids in the design of DNA sequences for folding 3D honeycomb-pleated shapes A series of rectangular-block motifs were designed, assembled, and analyzed to identify a well-behaved motif that could serve as a building block for future studies. The use of caDNAno significantly reduces the effort required to design 3D DNA-origami structures. The software is available at http://cadnano.org/, along with example designs and video tutorials demonstrating their construction. The source code is released under the MIT license.

Transcripts of unknown function in multiple-signaling pathways involved in human stem cell differentiation

2009-06-16

Mammalian transcriptome analysis has uncovered tens of thousands of novel transcripts of unknown function (TUFs). Classical and recent examples suggest that the majority of TUFs may underlie vital intracellular functions as non-coding RNAs because of their low coding potentials. However, only a portion of TUFs have been studied to date, and the functional significance of TUFs remains mostly uncharacterized. To increase the repertoire of functional TUFs, we screened for TUFs whose expression is controlled during differentiation of pluripotent human mesenchymal stem cells (hMSCs). The resulting six TUFs, named transcripts related to hMSC differentiation (TMDs), displayed distinct transcriptional kinetics during hMSC adipogenesis and/or osteogenesis. Structural and comparative genomic characterization suggested a wide variety of biologically active structures of these TMDs, including a long nuclear non-coding RNA, a microRNA host gene and a novel small protein gene. Moreover, the transcriptional response to established pathway activators indicated that most of these TMDs were transcriptionally regulated by each of the two key pathways for hMSC differentiation: the Wnt and protein kinase A (PKA) signaling pathways. The present study suggests that not only TMDs but also other human TUFs may in general participate in vital cellular functions with different molecular mechanisms.

Translation inhibition reveals interaction of 2'-deoxy and 2'-O-methyl molecular beacons with mRNA targets in living cells

Nitin, N., Rhee, W. J., Bao, G. — 2009-06-15

Understanding the interaction between oligonucleotide probes and RNA targets in living cells is important for biological and clinical studies of gene expression in vivo. Here, we demonstrate that starvation of cells and translation inhibition by blocking the mTOR or PI-3 kinase pathway could significantly reduce the fluorescence signal from 2'-deoxy molecular beacons (MBs) targeting K-ras and GAPDH mRNAs in living cells. However, the intensity and localization of fluorescence signal from MBs targeting nontranslated 28S rRNA remained the same in normal and translation-inhibited cells. We also found that, in targeting K-ras and GAPDH mRNAs, the signal level from MBs with 2'-O-methyl backbone did not change when translation was repressed. Taken together, our findings suggest that MBs with DNA backbone hybridize preferentially with mRNAs in their translational state in living cells, whereas those with 2'-O-methyl chemistry tend to hybridize to mRNA targets in both translational and nontranslational states. This work may thus provide a significant insight into probe design for detection of RNA molecules in living cells and RNA biology.

PRMT5 is required for cell-cycle progression and p53 tumor suppressor function

Scoumanne, A., Zhang, J., Chen, X. — 2009-06-15

Protein arginine methyltransferases (PRMTs) mediate the transfer of methyl groups to arginines in proteins involved in signal transduction, transcriptional regulation and RNA processing. Tumor suppressor p53 coordinates crucial cellular processes, including cell-cycle arrest and DNA repair, in response to stress signals. Post-translational modifications and interactions with co-factors are important to regulate p53 transcriptional activity. To explore whether PRMTs modulate p53 function, we generated multiple cell lines in which PRMT1, CARM1 and PRMT5 are inducibly knocked down. Here, we showed that PRMT5, but not PRMT1 or CARM1, is essential for cell proliferation and PRMT5 deficiency triggers cell-cycle arrest in G1. In addition, PRMT5 is required for p53 expression and induction of p53 targets MDM2 and p21 upon DNA damage. Importantly, we established that PRMT5 knockdown prevents p53 protein synthesis. Furthermore, we found that PRMT5 regulates the expression of translation initiation factor eIF4E and growth suppression mediated upon PRMT5 knockdown is independent of p53 but is dependent on eIF4E. Taken together, we uncovered that arginine methyltransferase PRMT5 is a major pro-survival factor regulating eIF4E expression and p53 translation.

cDNA display: a novel screening method for functional disulfide-rich peptides by solid-phase synthesis and stabilization of mRNA-protein fusions

2009-06-15

We report a robust display technology for the screening of disulfide-rich peptides, based on cDNA–protein fusions, by developing a novel and versatile puromycin-linker DNA. This linker comprises four major portions: a ‘ligation site’ for T4 RNA ligase, a ‘biotin site’ for solid-phase handling, a ‘reverse transcription primer site’ for the efficient and rapid conversion from an unstable mRNA–protein fusion (mRNA display) to a stable mRNA/cDNA–protein fusion (cDNA display) whose cDNA is covalently linked to its encoded protein and a ‘restriction enzyme site’ for the release of a complex from the solid support. This enables not only stabilizing mRNA–protein fusions but also promoting both protein folding and disulfide shuffling reactions. We evaluated the performance of cDNA display in different model systems and demonstrated an enrichment efficiency of 20-fold per selection round. Selection of a 32-residue random library against interleukin-6 receptor generated novel peptides containing multiple disulfide bonds with a unique linkage for its function. The peptides were found to bind with the target in the low nanomolar range. These results show the suitability of our method for in vitro selections of disulfide-rich proteins and other potential applications.

Studying human telomerase gene transcription by a chromatinized reporter generated by recombinase-mediated targeting of a bacterial artificial chromosome

Wang, S., Zhao, Y., Leiby, M. A., Zhu, J. — 2009-06-15

The endogenous human telomerase reverse transcriptase (hTERT) gene is repressed in somatic cells. To study the mechanisms of its repression, we developed a strategy of retrovirus-directed Cre recombinase-mediated BAC targeting, or RMBT, to generate single-copy integrations of BAC at pre-engineered chromosomal sites. This technique involved retroviral transduction of acceptor loci, containing an HSV thymidine kinase marker, and subsequent integration of BAC constructs into the acceptor sites, utilizing the loxP and lox511 sites present in the vector backbones. The BAC reporter, with a Renilla luciferase cassette inserted downstream of the hTERT promoter, was retrofitted with a puromycin marker. Through puromycin selection and ganciclovir counter-selection, a targeting efficiency of over 50% was achieved. We demonstrated that the activity and chromatin structures of the hTERT promoter in chromosomally integrated BAC reporter recapitulated its endogenous counterpart of the host cells. Therefore, we have established a genetically amendable platform to study chromatin and epigenetic regulation of the hTERT gene. The highly efficient and versatile RMBT technique has general applicability for studying largely unexplored chromatin-dependent mechanisms of promoter regulation of various genes.

Overestimation of alternative splicing caused by variable probe characteristics in exon arrays

Gaidatzis, D., Jacobeit, K., Oakeley, E. J., Stadler, M. B. — 2009-06-15

In higher eukaryotes, alternative splicing is a common mechanism for increasing transcriptome diversity. Affymetrix exon arrays were designed as a tool for monitoring the relative expression levels of hundreds of thousands of known and predicted exons with a view to detecting alternative splicing events. In this article, we have analyzed exon array data from many different human and mouse tissues and have uncovered a systematic relationship between transcript-fold change and alternative splicing as reported by the splicing index. Evidence from dilution experiments and deep sequencing suggest that this effect is of technical rather than biological origin and that it is driven by sequence features of the probes. This effect is substantial and results in a 12-fold overestimation of alternative splicing events in genes that are differentially expressed. By cross-species exon array comparison, we could further show that the systematic bias persists even across species boundaries. Failure to consider this effect in data analysis would result in the reproducible false detection of apparently conserved alternative splicing events. Finally, we have developed a software in R called COSIE (Corrected Splicing Indices for Exon arrays) that for any given set of new exon array experiments corrects for the observed bias and improves the detection of alternative splicing (available at www.fmi.ch/groups/gbioinfo).

Detection of single nucleotide variations in expressed exons of the human genome using RNA-Seq

Chepelev, I., Wei, G., Tang, Q., Zhao, K. — 2009-06-15

Whole-genome resequencing is still a costly method to detect genetic mutations that lead to altered forms of proteins and may be associated with disease development. Since the majority of disease-related single nucleotide variations (SNVs) are found in protein-coding regions, we propose to identify SNVs in expressed exons of the human genome using the recently developed RNA-Seq technique. We identify 12 176 and 10 621 SNVs, respectively, in Jurkat T cells and CD4⁺ T cells from a healthy donor. Interestingly, our data show that one copy of the TAL-1 proto-oncogene has a point mutation in 3' UTR and only the mutant allele is expressed in Jurkat cells. We provide a comprehensive dataset for further understanding the cancer biology of Jurkat cells. Our results indicate that this is a cost-effective and efficient strategy to systematically identify SNVs in the expressed regions of the human genome.

Differential Spo0A-mediated effects on transcription and replication of the related Bacillus subtilis phages Nf and {varphi}29 explain their different behaviours in vivo

Castilla-Llorente, V., Meijer, W. J. J., Salas, M. — 2009-06-15

Members of groups 1 (e.g. 29) and 2 (e.g. Nf) of the 29 family of phages infect the spore forming bacterium Bacillus subtilis. Although classified as lytic phages, the lytic cycle of 29 can be suppressed and its genome can become entrapped into the B. subtilis spore. This constitutes an alternative infection strategy that depends on the presence of binding sites for the host-encoded protein Spo0A in the 29 genome. Binding of Spo0A to these sites represses 29 transcription and prevents initiation of DNA replication. Although the Nf genome can also become trapped into B. subtilis spores, in vivo studies showed that its lytic cycle is less susceptible to spo0A-mediated suppression than that of 29. Here we have analysed the molecular mechanism underlying this difference showing that Spo0A differently affects transcription and replication initiation of the genomes of these phages. Thus, whereas Spo0A represses all three main early promoters of 29, it only represses one out of the three equivalent early promoters of Nf. In addition, contrary to 29, Spo0A does not prevent the in vitro initiation of Nf DNA replication. Altogether, the differences in Spo0A-mediated regulation of transcription and replication between 29 and Nf explain their different behaviours in vivo.

A novel nucleoid-associated protein of Mycobacterium tuberculosis is a sequence homolog of GroEL

Basu, D., Khare, G., Singh, S., Tyagi, A., Khosla, S., Mande, S. C. — 2009-06-15

The Mycobacterium tuberculosis genome sequence reveals remarkable absence of many nucleoid-associated proteins (NAPs), such as HNS, Hfq or DPS. In order to characterize the nucleoids of M. tuberculosis, we have attempted to identify NAPs, and report an interesting finding that a chaperonin-homolog, GroEL1, is nucleoid associated. We report that M. tuberculosis GroEL1 binds DNA with low specificity but high affinity, suggesting that it might have naturally evolved to bind DNA. We are able to demonstrate that GroEL1 can effectively function as a DNA-protecting agent against DNase I or hydroxyl-radicals. Moreover, Atomic Force Microscopic studies reveal that GroEL1 can condense a large DNA into a compact structure. We also provide in vivo evidences that include presence of GroEL1 in purified nucleoids, in vivo crosslinking followed by Southern hybridizations and immunofluorescence imaging in M. tuberculosis confirming that GroEL1: DNA interactions occur in natural biological settings. These findings therefore reveal that M. tuberculosis GroEL1 has evolved to be associated with nucleoids.

Translation of intronless RNAs is strongly stimulated by the Epstein-Barr virus mRNA export factor EB2

2009-06-15

The Epstein–Barr virus protein (EB2) allows the nuclear export of a particular subset of early and late viral RNAs derived from intronless genes. EB2 is conserved among most herpesvirus members and its presence is essential for the production of infectious particles. Here we show that, besides its role as a nuclear export factor, EB2 strongly stimulates translation of unspliced mRNAs without affecting overall cellular translation. Interestingly, this effect can be reversed by the addition of an intron within the gene. The spliced mRNA is then efficiently exported and translated even in the absence of EB2. Moreover, we show that EB2 associates with translating ribosomes and increases the proportion of its target RNA in the polyribosomal fraction. Finally, testing of EB2 homolog proteins derived from EBV-related herpesviruses, shows that, even if they play similar roles within the replication cycle of their respective virus, their mechanisms of action are different.

Gains and unexpected lessons from genome-scale promoter mapping

Shavkunov, K. S., Masulis, I. S., Tutukina, M. N., Deev, A. A., Ozoline, O. N. — 2009-06-15

Potential promoters in the genome of Escherichia coli were searched by pattern recognition software PlatProm and classified on the basis of positions relative to gene borders. Beside the expected promoters located in front of the coding sequences we found a considerable amount of intragenic promoter-like signals with a putative ability to drive either antisense or alternative transcription and revealed unusual genomic regions with extremely high density of predicted transcription start points (promoter ‘islands’), some of which are located in coding sequences. PlatProm scores converted into probability of RNA polymerase binding demonstrated certain correlation with the enzyme retention registered by ChIP-on-chip technique; however, in ‘dense’ regions the value of correlation coefficient is lower than throughout the entire genome. Experimental verification confirmed the ability of RNA polymerase to interact and form multiple open complexes within promoter ‘island’ associated with appY, yet transcription efficiency was lower than might be expected. Analysis of expression data revealed the same tendency for other promoter ‘islands’, thus assuming functional relevance of non-productive RNA polymerase binding. Our data indicate that genomic DNA of E. coli is enriched by numerous unusual promoter-like sites with biological role yet to be understood.

Activation of alternative Jdp2 promoters and functional protein isoforms in T-cell lymphomas by retroviral insertion mutagenesis

Rasmussen, M. H., Wang, B., Wabl, M., Nielsen, A. L., Pedersen, F. S. — 2009-06-15

Retroviral insertional mutagenesis has been instrumental for the identification of genes important in cancer development. The molecular mechanisms involved in retroviral-mediated activation of proto-oncogenes influence the distribution of insertions within specific regions during tumorigenesis and hence may point to novel gene structures. From a retroviral tagging screen on tumors of 1767 SL3-3 MLV-infected BALB/c mice, intron 2 of the AP-1 repressor Jdp2 locus was found frequently targeted by proviruses resulting in upregulation of non-canonical RNA subspecies. We identified several promoter regions within 1000 bp upstream of exon 3 that allowed for the production of Jdp2 protein isoforms lacking the histone acetylase inhibitory domain INHAT present in canonical Jdp2. The novel Jdp2 isoforms localized to the nucleus and over-expression in murine fibroblast cells induced cell death similar to canonic Jdp2. When expressed in the context of oncogenic NRAS both full length Jdp2 and the shorter isoforms increased anchorage-independent growth. Our results demonstrate a biological function of Jdp2 lacking the INHAT domain and suggest a post-genomic application for the use of retroviral tagging data in identifying new gene products with a potential role in tumorigenesis.

Structural implication of splicing stochastics

Melamud, E., Moult, J. — 2009-06-15

Even though nearly every human gene has at least one alternative splice form, very little is so far known about the structure and function of resulting protein products. It is becoming increasingly clear that a significant fraction of all isoforms are products of noisy selection of splice sites and thus contribute little to actual functional diversity, and may potentially be deleterious. In this study, we examine the impact of alternative splicing on protein sequence and structure in three datasets: alternative splicing events conserved across multiple species, alternative splicing events in genes that are strongly linked to disease and all observed alternative splicing events. We find that the vast majority of all alternative isoforms result in unstable protein conformations. In contrast to that, the small subset of isoforms conserved across species tends to maintain protein structural integrity to a greater extent. Alternative splicing in disease-associated genes produces unstable structures just as frequently as all other genes, indicating that selection to reduce the effects of alternative splicing on this set is not especially pronounced. Overall, the properties of alternative spliced proteins are consistent with the outcome of noisy selection of splice sites by splicing machinery.

Experimentally based contact energies decode interactions responsible for protein-DNA affinity and the role of molecular waters at the binding interface

Temiz, N. A., Camacho, C. J. — 2009-06-15

A major obstacle towards understanding the molecular basis of transcriptional regulation is the lack of a recognition code for protein–DNA interactions. Using high-quality crystal structures and binding data on the promiscuous family of C₂H₂ zinc fingers (ZF), we decode 10 fundamental specific interactions responsible for protein–DNA recognition. The interactions include five hydrogen bond types, three atomic desolvation penalties, a favorable non-polar energy, and a novel water accessibility factor. We apply this code to three large datasets containing a total of 89 C₂H₂ transcription factor (TF) mutants on the three ZFs of EGR. Guided by molecular dynamics simulations of individual ZFs, we map the interactions into homology models that embody all feasible intra- and intermolecular bonds, selecting for each sequence the structure with the lowest free energy. These interactions reproduce the change in affinity of 35 mutants of finger I (R² = 0.998), 23 mutants of finger II (R² = 0.96) and 31 finger III human domains (R² = 0.94). Our findings reveal recognition rules that depend on DNA sequence/structure, molecular water at the interface and induced fit of the C₂H₂ TFs. Collectively, our method provides the first robust framework to decode the molecular basis of TFs binding to DNA.

Estradiol-regulated microRNAs control estradiol response in breast cancer cells

2009-06-14

Estradiol (E2) regulates gene expression at the transcriptional level by functioning as a ligand for estrogen receptor alpha (ER) and estrogen receptor beta (ERβ). E2-inducible proteins c-Myc and E2Fs are required for optimal ER activity and secondary estrogen responses, respectively. We show that E2 induces 21 microRNAs and represses seven microRNAs in MCF-7 breast cancer cells; these microRNAs have the potential to control 420 E2-regulated and 757 non-E2-regulated mRNAs at the post-transcriptional level. The serine/threonine kinase, AKT, alters E2-regulated expression of microRNAs. E2 induced the expression of eight Let-7 family members, miR-98 and miR-21 microRNAs; these microRNAs reduced the levels of c-Myc and E2F2 proteins. Dicer, a ribonuclease III enzyme required for microRNA processing, is also an E2-inducible gene. Several E2-regulated microRNA genes are associated with ER-binding sites or located in the intragenic region of estrogen-regulated genes. We propose that the clinical course of ER-positive breast cancers is dependent on the balance between E2-regulated tumor-suppressor microRNAs and oncogenic microRNAs. Additionally, our studies reveal a negative-regulatory loop controlling E2 response through microRNAs as well as differences in E2-induced transcriptome and proteome.

Classification and energetics of the base-phosphate interactions in RNA

Zirbel, C. L., Sponer, J. E., Sponer, J., Stombaugh, J., Leontis, N. B. — 2009-06-14

Structured RNA molecules form complex 3D architectures stabilized by multiple interactions involving the nucleotide base, sugar and phosphate moieties. A significant percentage of the bases in structured RNA molecules in the Protein Data Bank (PDB) hydrogen-bond with phosphates of other nucleotides. By extracting and superimposing base-phosphate (BPh) interactions from a reduced-redundancy subset of 3D structures from the PDB, we identified recurrent phosphate-binding sites on the RNA bases. Quantum chemical calculations were carried out on model systems representing each BPh interaction. The calculations show that the centers of each cluster obtained from the structure superpositions correspond to energy minima on the potential energy hypersurface. The calculations also show that the most stable phosphate-binding sites occur on the Watson–Crick edge of guanine and the Hoogsteen edge of cytosine. We modified the ‘Find RNA 3D' (FR3D) software suite to automatically find and classify BPh interactions. Comparison of the 3D structures of the 16S and 23S rRNAs of Escherichia coli and Thermus thermophilus revealed that most BPh interactions are phylogenetically conserved and they occur primarily in hairpin, internal or junction loops or as part of tertiary interactions. Bases that form BPh interactions, which are conserved in the rRNA 3D structures are also conserved in homologous rRNA sequence alignments.

Developmentally regulated MAPK pathways modulate heterochromatin in Saccharomyces cerevisiae

Mazor, Y., Kupiec, M. — 2009-06-11

Variegated expression of genes contributes to phenotypic variation within populations of genetically identical cells. Such variation plays a role in development and host pathogen interaction and can be important in adaptation to harsh environments. The expression state of genes placed near telomeres shows a variegated pattern of inheritance due to heterochromatin formation, a phenomenon that is called telomere position effect (TPE). We show that in budding yeast, TPE is controlled by the a1/2 developmental repressor, which dictates developmental decisions in response to environmental changes. Two a1/2 repressed genes, STE5, a MAPK scaffold and HOG1, a stress-activated MAPK, are the targets of this heterochromatin regulation pathway. We provide new evidence that link MAPK signaling and heterochromatin formation in yeast. Our results show that the same components that regulate gene expression states in euchromatic regions regulate heterochromatic expression states and that stress can play a part in turning on or off genes placed in heterochromatic regions.

A novel androgen receptor-binding element modulates Cdc6 transcription in prostate cancer cells during cell-cycle progression

Jin, F., Fondell, J. D. — 2009-06-11

The androgen receptor (AR) plays a pivotal role in the onset and progression of prostate cancer by promoting cellular proliferation. Recent studies suggest AR is a master regulator of G1-S progression and possibly a licensing factor for DNA replication yet the mechanisms remain poorly defined. Here we report that AR targets the human Cdc6 gene for transcriptional regulation. Cdc6 is an essential regulator of DNA replication in eukaryotic cells and its mRNA expression is inversely modulated by androgen or antiandrogen treatment in androgen-sensitive prostate cancer cells. AR binds at a distinct androgen-response element (ARE) in the Cdc6 promoter that is functionally required for androgen-dependent Cdc6 transcription. We found that peak AR occupancy at the novel ARE occurs during the G1/S phase concomitant with peak Cdc6 mRNA expression. We also identified several of the coactivators and corepressors involved in AR-dependent Cdc6 transcriptional regulation in vivo and further characterized ligand-induced alterations in histone acetylation and methylation at the Cdc6 promoter. Significantly, AR silencing in prostate cancer cells markedly decreases Cdc6 expression and androgen-dependent cellular proliferation. Collectively, our results suggest that Cdc6 is a key regulatory target for AR and provide new insights into the mechanisms of prostate cancer cell proliferation.

Fur controls iron homeostasis and oxidative stress defense in the oligotrophic alpha-proteobacterium Caulobacter crescentus

da Silva Neto, J. F., Braz, V. S., Italiani, V. C. S., Marques, M. V. — 2009-06-11

In most bacteria, the ferric uptake regulator (Fur) is a global regulator that controls iron homeostasis and other cellular processes, such as oxidative stress defense. In this work, we apply a combination of bioinformatics, in vitro and in vivo assays to identify the Caulobacter crescentus Fur regulon. A C. crescentus fur deletion mutant showed a slow growth phenotype, and was hypersensitive to H₂O₂ and organic peroxide. Using a position weight matrix approach, several predicted Fur-binding sites were detected in the genome of C. crescentus, located in regulatory regions of genes not only involved in iron uptake and usage but also in other functions. Selected Fur-binding sites were validated using electrophoretic mobility shift assay and DNAse I footprinting analysis. Gene expression assays revealed that genes involved in iron uptake were repressed by iron-Fur and induced under conditions of iron limitation, whereas genes encoding iron-using proteins were activated by Fur under conditions of iron sufficiency. Furthermore, several genes that are regulated via small RNAs in other bacteria were found to be directly regulated by Fur in C. crescentus. In conclusion, Fur functions as an activator and as a repressor, integrating iron metabolism and oxidative stress response in C. crescentus.

Repressor CopG prevents access of RNA polymerase to promoter and actively dissociates open complexes

Hernandez-Arriaga, A. M., Rubio-Lepe, T. S., Espinosa, M., del Solar, G. — 2009-06-11

Replication of the promiscuous plasmid pMV158 requires expression of the initiator repB gene, which is controlled by the repressor CopG. Genes repB and copG are co-transcribed from promoter P_cr. We have studied the interactions between RNA polymerase, CopG and the promoter to elucidate the mechanism of repression by CopG. Complexes formed at 0°C and at 37°C between RNA polymerase and P_cr differed from each other in stability and in the extent of the DNA contacted. The 37°C complex was very stable (half-life of about 3 h), and shared features with typical open complexes generated at a variety of promoters. CopG protein repressed transcription from P_cr at two different stages in the process leading to the initiation complex. First, CopG hindered binding of RNA polymerase to the promoter. Second, CopG was able to displace RNA polymerase once the enzyme has formed a stable complex with P_cr. A model for the CopG-mediated disassembly of the stable RNA polymerase–P_cr promoter complex is presented.

Functional importance of Crenarchaea-specific extra-loop revealed by an X-ray structure of a heterotetrameric crenarchaeal splicing endonuclease

2009-06-10

Archaeal splicing endonucleases (EndAs) are currently classified into three groups. Two groups require a single subunit protein to form a homodimer or homotetramer. The third group requires two nonidentical protein components for the activity. To elucidate the molecular architecture of the two-subunit EndA system, we studied a crenarchaeal splicing endonuclease from Pyrobaculum aerophilum. In the present study, we solved a crystal structure of the enzyme at 1.7-Å resolution. The enzyme adopts a heterotetrameric form composed of two catalytic and two structural subunits. By connecting the structural and the catalytic subunits of the heterotetrameric EndA, we could convert the enzyme to a homodimer that maintains the broad substrate specificity that is one of the characteristics of heterotetrameric EndA. Meanwhile, a deletion of six amino acids in a Crenarchaea-specific loop abolished the endonuclease activity even on a substrate with canonical BHB motif. These results indicate that the subunit architecture is not a major factor responsible for the difference of substrate specificity between single- and two-subunit EndA systems. Rather, the structural basis for the broad substrate specificity is built into the crenarchaeal splicing endonuclease itself.

Conditional gene silencing of multiple genes with antisense RNAs and generation of a mutator strain of Escherichia coli

Nakashima, N., Tamura, T. — 2009-06-10

In this study, we describe a method of simultaneous conditional gene silencing of up to four genes in Escherichia coli by using antisense RNAs. We used antisense RNAs with paired termini, which carried flanking inverted repeats to create paired double-stranded RNA termini; these RNAs have been proven to have high silencing efficacy. To express antisense RNAs, we constructed four IPTG-inducible vectors carrying different but compatible replication origins. When the lacZ antisense RNA was expressed using these vectors, lacZ expression was successfully silenced by all the vectors, but the expression level of the antisense RNA and silencing efficacy differed depending on the used vectors. All the vectors were co-transformable; the antisense RNAs against lacZ, ackA, pta and pepN were co-expressed, and silencing of all the target genes was confirmed. Furthermore, when antisense RNAs were targeted to the mutator genes mutS, mutD (dnaQ) and ndk, which are involved in DNA replication or DNA mismatch repair, spontaneous mutation frequencies increased over 2000-fold. The resulting mutator strain is useful for random mutagenesis of plasmids. The method provides a robust tool for investigating functional relationships between multiple genes or altering cell phenotypes for biotechnological and industrial applications.

Plasmid-based lacZ{alpha} assay for DNA polymerase fidelity: application to archaeal family-B DNA polymerase

Jozwiakowski, S. K., Connolly, B. A. — 2009-06-10

The preparation of a gapped pUC18 derivative, containing the lacZ reporter gene in the single-stranded region, is described. Gapping is achieved by flanking the lacZ gene with sites for two related nicking endonucleases, enabling the excision of either the coding or non-coding strand. However, the excised strand remains annealed to the plasmid through non-covalent Watson–Crick base-pairing; its removal, therefore, requires a heat–cool cycle in the presence of an exactly complementary competitor DNA. The gapped plasmids can be used to assess DNA polymerase fidelity using in vitro replication, followed by transformation into Escherichia coli and scoring the blue/white colony ratio. Results found with plasmids are similar to the well established method based on gapped M13, in terms of background (~0.08% in both cases) and the mutation frequencies observed with a number of DNA polymerases, providing validation for this straightforward and technically uncomplicated approach. Several error prone variants of the archaeal family-B DNA polymerase from Pyrococcus furiosus have been investigated, illuminating the potential of the method.

Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch

Villa, A., Wohnert, J., Stock, G. — 2009-06-10

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA's secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (1 µs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.

DNA binding and synapsis by the large C-terminal domain of {phi}C31 integrase

McEwan, A. R., Rowley, P. A., Smith, M. C. M. — 2009-06-10

The integrase (Int) from phage C31 acts on the phage and host-attachment sites, attP and attB, to form an integrated prophage flanked by attL and attR. Excision (attL x attR recombination) is prevented, in the absence of accessory factors, by a putative coiled-coil motif in the C-terminal domain (CTD). Int has a serine recombinase N-terminal domain, required for synapsis of recombination substrates and catalysis. We show here that the coiled-coil motif mediates protein–protein interactions between CTDs, but only when bound to DNA. Although the histidine-tagged CTD (hCTD) was monomeric in solution, hCTD bound cooperatively to three of the recombination substrates (attB, attL and attR). Furthermore, when provided with attP and attB, hCTD brought these substrates together in a synaptic complex. Substitutions in the coiled-coil motif that greatly reduce Int integration activity, L460P and Y475H, prevented CTD–CTD interactions and led to defective DNA binding and no detectable DNA synapsis. A substitution, E449K, in full length Int confers the ability to perform excision in addition to integration as it has gained the ability to synapse attL x attR. hCTD^E449K was similar to hCTD in DNA binding but unable to form the CTD synapse suggesting that the CTD synapse is not essential but could be part of the mechanism that controls directionality.

Imaging individual microRNAs in single mammalian cells in situ

Lu, J., Tsourkas, A. — 2009-06-10

MicroRNAs (miRNAs) are potent negative regulators of gene expression that have been implicated in most major cellular processes. Despite rapid advances in our understanding of miRNA biogenesis and mechanism, many fundamental questions still remain regarding miRNA function and their influence on cell cycle control. Considering recent reports on the impact of cell-to-cell fluctuations in gene expression on phenotypic diversity, it is likely that looking at the average miRNA expression of cell populations could result in the loss of important information connecting miRNA expression and cell function. Currently, however, there are no efficient techniques to quantify miRNA expression at the single-cell level. Here, a method is described for the detection of individual miRNA molecules in cancer cells using fluorescence in situ hybridization. The method combines the unique recognition properties of locked nucleic acid probes with enzyme-labeled fluorescence. Using this approach, individual miRNAs are identified as bright, photostable fluorescent spots. In this study, miR-15a was quantified in MDA-MB-231 and HeLa cells, while miR-155 was quantified in MCF-7 cells. The dynamic range was found to span over three orders of magnitude and the average miRNA copy number per cell was within 17.5% of measurements acquired by quantitative RT-PCR.

The effects of unnatural base pairs and mispairs on DNA duplex stability and solvation

Hwang, G. T., Hari, Y., Romesberg, F. E. — 2009-06-10

In an effort to develop unnatural DNA base pairs we examined six pyridine-based nucleotides, d3MPy, d4MPy, d5MPy, d34DMPy, d35DMPy and d45DMPy. Each bears a pyridyl nucleobase scaffold but they are differentiated by methyl substitution, and were designed to vary both inter- and intra-strand packing within duplex DNA. The effects of the unnatural base pairs on duplex stability demonstrate that the pyridine scaffold may be optimized for stable and selective pairing, and identify one self pair, the pair formed between two d34DMPy nucleotides, which is virtually as stable as a dA:dT base pair in the same sequence context. In addition, we found that the incorporation of either the d34DMPy self pair or a single d34DMPy paired opposite a natural dA significantly increases oligonucleotide hybridization fidelity at other positions within the duplex. Hypersensitization of the duplex to mispairing appears to result from global and interdependent solvation effects mediated by the unnatural nucleotide(s) and the mispair. The results have important implications for our efforts to develop unnatural base pairs and suggest that the unnatural nucleotides might be developed as novel biotechnological tools, diagnostics, or therapeutics for applications where hybridization stringency is important.

The Hin recombinase assembles a tetrameric protein swivel that exchanges DNA strands

Dhar, G., McLean, M. M., Heiss, J. K., Johnson, R. C. — 2009-06-10

Most site-specific recombinases can be grouped into two structurally and mechanistically different classes. Whereas recombination by tyrosine recombinases proceeds with little movements by the proteins, serine recombinases exchange DNA strands by a mechanism requiring large quaternary rearrangements. Here we use site-directed crosslinking to investigate the conformational changes that accompany the formation of the synaptic complex and the exchange of DNA strands by the Hin serine recombinase. Efficient crosslinking between residues corresponding to the ‘D-helix’ region provides the first experimental evidence for interactions between synapsed subunits within this region and distinguishes between different tetrameric conformers that have been observed in crystal structures of related serine recombinases. Crosslinking profiles between cysteines introduced over the 35 residue E-helix region that constitutes most of the proposed rotating interface both support the long helical structure of the region and provide strong experimental support for a subunit rotation mechanism that mediates DNA exchange.

RNA-binding specificity of E. coli NusA

Prasch, S., Jurk, M., Washburn, R. S., Gottesman, M. E., Wohrl, B. M., Rosch, P. — 2009-06-10

The RNA sequences boxA, boxB and boxC constitute the nut regions of phage . They nucleate the formation of a termination-resistant RNA polymerase complex on the chromosome. The complex includes E. coli proteins NusA, NusB, NusG and NusE, and the N protein. A complex that includes the Nus proteins and other factors forms at the rrn leader. Whereas RNA-binding by NusB and NusE has been described in quantitative terms, the interaction of NusA with these RNA sequences is less defined. Isotropic as well as anisotropic fluorescence equilibrium titrations show that NusA binds only the nut spacer sequence between boxA and boxB. Thus, nutR boxA5-spacer, nutR boxA16-spacer and nutR boxA69-spacer retain NusA binding, whereas a spacer mutation eliminates complex formation. The affinity of NusA for nutL is 50% higher than for nutR. In contrast, rrn boxA, which includes an additional U residue, binds NusA in the absence of spacer. The K_d values obtained for rrn boxA and rrn boxA-spacer are 19-fold and 8-fold lower, respectively, than those for nutR boxA-spacer. These differences may explain why requires an additional protein, N, to suppress termination. Knowledge of the different affinities now describes the assembly of the anti-termination complex in quantitative terms.

ST1710-DNA complex crystal structure reveals the DNA binding mechanism of the MarR family of regulators

Kumarevel, T., Tanaka, T., Umehara, T., Yokoyama, S. — 2009-06-09

ST1710, a member of the multiple antibiotic resistance regulator (MarR) family of regulatory proteins in bacteria and archaea, plays important roles in development of antibiotic resistance, a global health problem. Here, we present the crystal structure of ST1710 from Sulfolobus tokodaii strain 7 complexed with salicylate, a well-known inhibitor of MarR proteins and the ST1710 complex with its promoter DNA, refined to 1.8 and 2.10 Å resolutions, respectively. The ST1710–DNA complex shares the topology of apo-ST1710 and MarR proteins, with each subunit containing a winged helix-turn-helix (wHtH) DNA binding motif. Significantly large conformational changes occurred upon DNA binding and in each of the dimeric monomers in the asymmetric unit of the ST1710–DNA complex. Conserved wHtH loop residues interacting with the bound DNA and mutagenic analysis indicated that R89, R90 and K91 were important for DNA recognition. Significantly, the bound DNA exhibited a new binding mechanism.

QUBIC: a qualitative biclustering algorithm for analyses of gene expression data

Li, G., Ma, Q., Tang, H., Paterson, A. H., Xu, Y. — 2009-06-09

Biclustering extends the traditional clustering techniques by attempting to find (all) subgroups of genes with similar expression patterns under to-be-identified subsets of experimental conditions when applied to gene expression data. Still the real power of this clustering strategy is yet to be fully realized due to the lack of effective and efficient algorithms for reliably solving the general biclustering problem. We report a QUalitative BIClustering algorithm (QUBIC) that can solve the biclustering problem in a more general form, compared to existing algorithms, through employing a combination of qualitative (or semi-quantitative) measures of gene expression data and a combinatorial optimization technique. One key unique feature of the QUBIC algorithm is that it can identify all statistically significant biclusters including biclusters with the so-called ‘scaling patterns’, a problem considered to be rather challenging; another key unique feature is that the algorithm solves such general biclustering problems very efficiently, capable of solving biclustering problems with tens of thousands of genes under up to thousands of conditions in a few minutes of the CPU time on a desktop computer. We have demonstrated a considerably improved biclustering performance by our algorithm compared to the existing algorithms on various benchmark sets and data sets of our own. QUBIC was written in ANSI C and tested using GCC (version 4.1.2) on Linux. Its source code is available at: http://csbl.bmb.uga.edu/~maqin/bicluster. A server version of QUBIC is also available upon request.

Kinetic and thermodynamic characterization of single-mismatch discrimination using single-molecule imaging

Gunnarsson, A., Jonsson, P., Zhdanov, V. P., Hook, F. — 2009-06-09

A single-molecule detection setup based on total internal reflection fluorescence (TIRF) microscopy has been used to investigate association and dissociation kinetics of unlabeled 30mer DNA strands. Single-molecule sensitivity was accomplished by letting unlabeled DNA target strands mediate the binding of DNA-modified and fluorescently labeled liposomes to a DNA-modified surface. The liposomes, acting as signal enhancer elements, enabled the number of binding events as well as the residence time for high affinity binders (K_d < 1 nM, k_off < 0.01 s^–1) to be collected under equilibrium conditions at low pM concentrations. The mismatch discrimination obtained from the residence time data was shown to be concentration and temperature independent in intervals of 1–100 pM and 23–46°C, respectively. This suggests the method as a robust means for detection of point mutations at low target concentrations in, for example, single nucleotide polymorphism (SNP) analysis.

Using DNA mechanics to predict in vitro nucleosome positions and formation energies

2009-06-09

In eukaryotic genomes, nucleosomes function to compact DNA and to regulate access to it both by simple physical occlusion and by providing the substrate for numerous covalent epigenetic tags. While competition with other DNA-binding factors and action of chromatin remodeling enzymes significantly affect nucleosome formation in vivo, nucleosome positions in vitro are determined by steric exclusion and sequence alone. We have developed a biophysical model, DNABEND, for the sequence dependence of DNA bending energies, and validated it against a collection of in vitro free energies of nucleosome formation and a set of in vitro nucleosome positions mapped at high resolution. We have also made a first ab initio prediction of nucleosomal DNA geometries, and checked its accuracy against the nucleosome crystal structure. We have used DNABEND to design both strong and weak histone- binding sequences, and measured the corresponding free energies of nucleosome formation. We find that DNABEND can successfully predict in vitro nucleosome positions and free energies, providing a physical explanation for the intrinsic sequence dependence of histone–DNA interactions.

Thermodynamic characterization of tandem mismatches found in naturally occurring RNA

Christiansen, M. E., Znosko, B. M. — 2009-06-09

Although all sequence symmetric tandem mismatches and some sequence asymmetric tandem mismatches have been thermodynamically characterized and a model has been proposed to predict the stability of previously unmeasured sequence asymmetric tandem mismatches [Christiansen,M.E. and Znosko,B.M. (2008) Biochemistry, 47, 4329–4336], experimental thermodynamic data for frequently occurring tandem mismatches is lacking. Since experimental data is preferred over a predictive model, the thermodynamic parameters for 25 frequently occurring tandem mismatches were determined. These new experimental values, on average, are 1.0 kcal/mol different from the values predicted for these mismatches using the previous model. The data for the sequence asymmetric tandem mismatches reported here were then combined with the data for 72 sequence asymmetric tandem mismatches that were published previously, and the parameters used to predict the thermodynamics of previously unmeasured sequence asymmetric tandem mismatches were updated. The average absolute difference between the measured values and the values predicted using these updated parameters is 0.5 kcal/mol. This updated model improves the prediction for tandem mismatches that were predicted rather poorly by the previous model. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA duplexes containing tandem mismatches, and, furthermore, should allow for improved prediction of secondary structure from sequence.

Extensive DNA mimicry by the ArdA anti-restriction protein and its role in the spread of antibiotic resistance

2009-06-08

The ardA gene, found in many prokaryotes including important pathogenic species, allows associated mobile genetic elements to evade the ubiquitous Type I DNA restriction systems and thereby assist the spread of resistance genes in bacterial populations. As such, ardA contributes to a major healthcare problem. We have solved the structure of the ArdA protein from the conjugative transposon Tn916 and find that it has a novel extremely elongated curved cylindrical structure with defined helical grooves. The high density of aspartate and glutamate residues on the surface follow a helical pattern and the whole protein mimics a 42-base pair stretch of B-form DNA making ArdA by far the largest DNA mimic known. Each monomer of this dimeric structure comprises three alpha–beta domains, each with a different fold. These domains have the same fold as previously determined proteins possessing entirely different functions. This DNA mimicry explains how ArdA can bind and inhibit the Type I restriction enzymes and we demonstrate that 6 different ardA from pathogenic bacteria can function in Escherichia coli hosting a range of different Type I restriction systems.

The incorporation of the novel histone variant H2AL2 confers unusual structural and functional properties of the nucleosome

2009-06-08

In this work we have studied the properties of the novel mouse histone variant H2AL2. H2AL2 was used to reconstitute nucleosomes and the structural and functional properties of these particles were studied by a combination of biochemical approaches, atomic force microscopy (AFM) and electron cryo-microscopy. DNase I and hydroxyl radical footprinting as well as micrococcal and exonuclease III digestion demonstrated an altered structure of the H2AL2 nucleosomes all over the nucleosomal DNA length. Restriction nuclease accessibility experiments revealed that the interactions of the H2AL2 histone octamer with the ends of the nucleosomal DNA are highly perturbed. AFM imaging showed that the H2AL2 histone octamer was complexed with only ~130 bp of DNA. H2AL2 reconstituted trinucleosomes exhibited a type of a ‘beads on a string’ structure, which was quite different from the equilateral triangle 3D organization of conventional H2A trinucleosomes. The presence of H2AL2 affected both the RSC and SWI/SNF remodeling and mobilization of the variant particles. These unusual properties of the H2AL2 nucleosomes suggest a specific role of H2AL2 during mouse spermiogenesis.

A physical and functional link between splicing factors promotes pre-mRNA 3' end processing

2009-06-08

Polypyrimidine tract-binding protein (PTB) is a splicing regulator that also plays a positive role in pre-mRNA 3' end processing when bound upstream of the polyadenylation signal (pA signal). Here, we address the mechanism of PTB stimulatory function in mRNA 3' end formation. We identify PTB as the protein factor whose binding to the human β-globin (HBB) 3' UTR is abrogated by a 3' end processing-inactivating mutation. We show that PTB promotes both in vitro 3' end cleavage and polyadenylation and recruits directly the splicing factor hnRNP H to G-rich sequences associated with several pA signals. Increased binding of hnRNP H results in stimulation of polyadenylation through a direct interaction with poly(A) polymerase. Therefore, our results provide evidence of a concerted regulation of pA signal recognition by splicing factors bound to auxiliary polyadenylation sequence elements.

Conferring a template-dependent polymerase activity to terminal deoxynucleotidyltransferase by mutations in the Loop1 region

Romain, F., Barbosa, I., Gouge, J., Rougeon, F., Delarue, M. — 2009-06-05

Terminal deoxynucleotidyltransferase (Tdt) and DNA polymerase µ (pol µ) are two eukaryotic highly similar proteins involved in DNA processing and repair. Despite their high sequence identity, they differ widely in their activity: pol µ has a templated polymerase activity, whereas Tdt has a non-templated one. Loop1, first described when the Tdt structure was solved, has been invoked as the major structural determinant of this difference. Here we describe attempts to transform Tdt into pol µ with the minimal number of mutations in and around Loop1. First we describe the effect of mutations on six different positions chosen to destabilize Tdt Loop1 structure, either by alanine substitution or by deletion; they result at most in a reduction of Tdt activity, but adding Co⁺⁺ restores most of this Tdt activity. However, a deletion of the entire Loop1 as in pol does confer a limited template-dependent polymerase behavior to Tdt while a chimera bearing an extended pol µ Loop1 reproduces pol µ behavior. Finally, 16 additional substitutions are reported, targeted at the two so-called ‘sequence determinant’ regions located just after Loop1 or underneath. Among them, the single-point mutant F401A displays a sequence-specific replicative polymerase phenotype that is stable upon Co⁺⁺ addition. These results are discussed in light of the available crystal structures.

Using genome-wide measurements for computational prediction of SH2-peptide interactions

Wunderlich, Z., Mirny, L. A. — 2009-06-05

Peptide-recognition modules (PRMs) are used throughout biology to mediate protein–protein interactions, and many PRMs are members of large protein domain families. Recent genome-wide measurements describe networks of peptide–PRM interactions. In these networks, very similar PRMs recognize distinct sets of peptides, raising the question of how peptide-recognition specificity is achieved using similar protein domains. The analysis of individual protein complex structures often gives answers that are not easily applicable to other members of the same PRM family. Bioinformatics-based approaches, one the other hand, may be difficult to interpret physically. Here we integrate structural information with a large, quantitative data set of SH2 domain–peptide interactions to study the physical origin of domain–peptide specificity. We develop an energy model, inspired by protein folding, based on interactions between the amino-acid positions in the domain and peptide. We use this model to successfully predict which SH2 domains and peptides interact and uncover the positions in each that are important for specificity. The energy model is general enough that it can be applied to other members of the SH2 family or to new peptides, and the cross-validation results suggest that these energy calculations will be useful for predicting binding interactions. It can also be adapted to study other PRM families, predict optimal peptides for a given SH2 domain, or study other biological interactions, e.g. protein–DNA interactions.

Combining use of a panel of ssDNA aptamers in the detection of Staphylococcus aureus

2009-06-04

In this article, a panel of ssDNA aptamers specific to Staphylococcus aureus was obtained by a whole bacterium-based SELEX procedure and applied to probing S. aureus. After several rounds of selection with S. aureus as the target and Streptococcus and S. epidermidis as counter targets, the highly enriched oligonucleic acid pool was sequenced and then grouped under different families on the basis of the homology of the primary sequence and the similarity of the secondary structure. Eleven sequences from different families were selected for further characterization by confocal imaging and flow cytometry analysis. Results showed that five aptamers demonstrated high specificity and affinity to S. aureus individually. The five aptamers recognize different molecular targets by competitive experiment. Combining these five aptamers had a much better effect than the individual aptamer in the recognition of different S. aureus strains. In addition, the combined aptamers can probe single S. aureus in pyogenic fluids. Our work demonstrates that a set of aptamers specific to one bacterium can be used in combination for the identification of the bacterium instead of a single aptamer.

Cavities in protein-DNA and protein-RNA interfaces

Sonavane, S., Chakrabarti, P. — 2009-06-03

An analysis of cavities present in protein–DNA and protein–RNA complexes is presented. In terms of the number of cavities and their total volume, the interfaces formed in these complexes are akin to those in transient protein–protein heterocomplexes. With homodimeric proteins protein–DNA interfaces may contain cavities involving both the protein subunits and DNA, and these are more than twice as large as cavities involving a single protein subunit and DNA. A parameter, cavity index, measuring the degree of surface complementarity, indicates that the packing of atoms in protein–protein/DNA/RNA is very similar, but it is about two times less efficient in the permanent interfaces formed between subunits in homodimers. As within the tertiary structure and protein–protein interfaces, protein–DNA interfaces have a higher inclination to be lined by β-sheet residues; from the DNA side, base atoms, in particular those in minor grooves, have a higher tendency to be located in cavities. The larger cavities tend to be less spherical and solvated. A small fraction of water molecules are found to mediate hydrogen-bond interactions with both the components, suggesting their primary role is to fill in the void left due to the local non-complementary nature of the surface patches.

Studying vertebrate topoisomerase 2 function using a conditional knockdown system in DT40 cells

Johnson, M., Phua, H. H., Bennett, S. C., Spence, J. M., Farr, C. J. — 2009-06-03

DT40 is a B-cell lymphoma-derived avian cell line widely used to study cell autonomous gene function because of the high rates with which DNA constructs are homologously recombined into its genome. Here, we demonstrate that the power of the DT40 system can be extended yet further through the use of RNA interference as an alternative to gene targeting. We have generated and characterized stable DT40 transfectants in which both topo 2 genes have been in situ tagged using gene targeting, and from which the mRNA of both topoisomerase 2 isoforms can be conditionally depleted through the tetracycline-induced expression of short hairpin RNAs. The cell cycle phenotype of topo 2-depleted DT40 cells has been compared with that previously reported for other vertebrate cells depleted either of topo 2 through gene targeting, or depleted of both isoforms simultaneously by transient RNAi. In addition, the DT40 knockdown system has been used to explore whether excess catenation arising through topo 2 depletion is sufficient to trigger the G2 catenation (or decatenation) checkpoint, proposed to exist in differentiated vertebrate cells.

Intra-tumor heterogeneity of MLH1 promoter methylation revealed by deep single molecule bisulfite sequencing

Varley, K. E., Mutch, D. G., Edmonston, T. B., Goodfellow, P. J., Mitra, R. D. — 2009-06-03

A single tumor may contain cells with different somatic mutations. By characterizing this genetic heterogeneity within tumors, advances have been made in the prognosis, treatment and understanding of tumorigenesis. In contrast, the extent of epigenetic intra-tumor heterogeneity and how it influences tumor biology is under-explored. We have characterized epigenetic heterogeneity within individual tumors using next-generation sequencing. We used deep single molecule bisulfite sequencing and sample-specific DNA barcodes to determine the spectrum of MLH1 promoter methylation across an average of 1000 molecules in each of 33 individual samples in parallel, including endometrial cancer, matched blood and normal endometrium. This first glimpse, deep into each tumor, revealed unexpectedly heterogeneous patterns of methylation at the MLH1 promoter within a subset of endometrial tumors. This high-resolution analysis allowed us to measure the clonality of methylation in individual tumors and gain insight into the accumulation of aberrant promoter methylation on both alleles during tumorigenesis.

From Corynebacterium glutamicum to Mycobacterium tuberculosis--towards transfers of gene regulatory networks and integrated data analyses with MycoRegNet

Krawczyk, J., Kohl, T. A., Goesmann, A., Kalinowski, J., Baumbach, J. — 2009-06-03

Year by year, approximately two million people die from tuberculosis, a disease caused by the bacterium Mycobacterium tuberculosis. There is a tremendous need for new anti-tuberculosis therapies (antituberculotica) and drugs to cope with the spread of tuberculosis. Despite many efforts to obtain a better understanding of M. tuberculosis' pathogenicity and its survival strategy in humans, many questions are still unresolved. Among other cellular processes in bacteria, pathogenicity is controlled by transcriptional regulation. Thus, various studies on M. tuberculosis concentrate on the analysis of transcriptional regulation in order to gain new insights on pathogenicity and other essential processes ensuring mycobacterial survival. We designed a bioinformatics pipeline for the reliable transfer of gene regulations between taxonomically closely related organisms that incorporates (i) a prediction of orthologous genes and (ii) the prediction of transcription factor binding sites. In total, 460 regulatory interactions were identified for M. tuberculosis using our comparative approach. Based on that, we designed a publicly available platform that aims to data integration, analysis, visualization and finally the reconstruction of mycobacterial transcriptional gene regulatory networks: MycoRegNet. It is a comprehensive database system and analysis platform that offers several methods for data exploration and the generation of novel hypotheses. MycoRegNet is publicly available at http://mycoregnet.cebitec.uni-bielefeld.de.

An integrative genomics approach identifies Hypoxia Inducible Factor-1 (HIF-1)-target genes that form the core response to hypoxia

Benita, Y., Kikuchi, H., Smith, A. D., Zhang, M. Q., Chung, D. C., Xavier, R. J. — 2009-06-02

The transcription factor Hypoxia-inducible factor 1 (HIF-1) plays a central role in the transcriptional response to oxygen flux. To gain insight into the molecular pathways regulated by HIF-1, it is essential to identify the downstream-target genes. We report here a strategy to identify HIF-1-target genes based on an integrative genomic approach combining computational strategies and experimental validation. To identify HIF-1-target genes microarrays data sets were used to rank genes based on their differential response to hypoxia. The proximal promoters of these genes were then analyzed for the presence of conserved HIF-1-binding sites. Genes were scored and ranked based on their response to hypoxia and their HIF-binding site score. Using this strategy we recovered 41% of the previously confirmed HIF-1-target genes that responded to hypoxia in the microarrays and provide a catalogue of predicted HIF-1 targets. We present experimental validation for ANKRD37 as a novel HIF-1-target gene. Together these analyses demonstrate the potential to recover novel HIF-1-target genes and the discovery of mammalian-regulatory elements operative in the context of microarray data sets.

Measuring spatial preferences at fine-scale resolution identifies known and novel cis-regulatory element candidates and functional motif-pair relationships

Yokoyama, K. D., Ohler, U., Wray, G. A. — 2009-06-02

Transcriptional regulation is mediated by the collective binding of proteins called transcription factors to cis-regulatory elements. A handful of factors are known to function at particular distances from the transcription start site, although the extent to which this occurs is not well understood. Spatial dependencies can also exist between pairs of binding motifs, facilitating factor-pair interactions. We sought to determine to what extent spatial preferences measured at high-scale resolution could be utilized to predict cis-regulatory elements as well as motif-pairs binding interacting proteins. We introduce the ‘motif positional function’ model which predicts spatial biases using regression analysis, differentiating noise from true position-specific overrepresentation at single-nucleotide resolution. Our method predicts 48 consensus motifs exhibiting positional enrichment within human promoters, including fourteen motifs without known binding partners. We then extend the model to analyze distance preferences between pairs of motifs. We find that motif-pairs binding interacting factors often co-occur preferentially at multiple distances, with intervals between preferred distances often corresponding to the turn of the DNA double-helix. This offers a novel means by which to predict sequence elements with a collective role in gene regulation.

Gene silencing in the marine diatom Phaeodactylum tricornutum

De Riso, V., Raniello, R., Maumus, F., Rogato, A., Bowler, C., Falciatore, A. — 2009-05-31

Diatoms are a major but poorly understood phytoplankton group. The recent completion of two whole genome sequences has revealed that they contain unique combinations of genes, likely recruited during their history as secondary endosymbionts, as well as by horizontal gene transfer from bacteria. A major limitation for the study of diatom biology and gene function is the lack of tools to generate targeted gene knockout or knockdown mutants. In this work, we have assessed the possibility of triggering gene silencing in Phaeodactylum tricornutum using constructs containing either anti-sense or inverted repeat sequences of selected target genes. We report the successful silencing of a GUS reporter gene expressed in transgenic lines, as well as the knockdown of endogenous phytochrome (DPH1) and cryptochrome (CPF1) genes. To highlight the utility of the approach we also report the first phenotypic characterization of a diatom mutant (cpf1). Our data open the way for reverse genetics in diatoms and represent a major advance for understanding their biology and ecology. Initial molecular analyses reveal that targeted downregulation likely occurs through transcriptional and post-transcriptional gene silencing mechanisms. Interestingly, molecular players involved in RNA silencing in other eukaryotes are only poorly conserved in diatoms.

The flexibility of locally melted DNA

Forties, R. A., Bundschuh, R., Poirier, M. G. — 2009-05-31

Protein-bound duplex DNA is often bent or kinked. Yet, quantification of intrinsic DNA bending that might lead to such protein interactions remains enigmatic. DNA cyclization experiments have indicated that DNA may form sharp bends more easily than predicted by the established worm-like chain (WLC) model. One proposed explanation suggests that local melting of a few base pairs introduces flexible hinges. We have expanded this model to incorporate sequence and temperature dependence of the local melting, and tested it for three sequences at temperatures from 23°C to 42°C. We find that small melted bubbles are significantly more flexible than double-stranded DNA and can alter DNA flexibility at physiological temperatures. However, these bubbles are not flexible enough to explain the recently observed very sharp bends in DNA.

Tracking transcription factor complexes on DNA using total internal reflectance fluorescence protein binding microarrays

Bonham, A. J., Neumann, T., Tirrell, M., Reich, N. O. — 2009-05-31

We have developed a high-throughput protein binding microarray (PBM) assay to systematically investigate transcription regulatory protein complexes binding to DNA with varied specificity and affinity. Our approach is based on the novel coupling of total internal reflectance fluorescence (TIRF) spectroscopy, swellable hydrogel double-stranded DNA microarrays and dye-labeled regulatory proteins, making it possible to determine both equilibrium binding specificities and kinetic rates for multiple protein:DNA interactions in a single experiment. DNA specificities and affinities for the general transcription factors TBP, TFIIA and IIB determined by TIRF–PBM are similar to those determined by traditional methods, while simultaneous measurement of the factors in binary and ternary protein complexes reveals preferred binding combinations. TIRF–PBM provides a novel and extendible platform for multi-protein transcription factor investigation.

Sexy gene conversions: locating gene conversions on the X-chromosome

Lawson, M. J., Zhang, L. — 2009-05-31

Gene conversion can have a profound impact on both the short- and long-term evolution of genes and genomes. Here, we examined the gene families that are located on the X-chromosomes of human (Homo sapiens), chimpanzee (Pan troglodytes), mouse (Mus musculus) and rat (Rattus norvegicus) for evidence of gene conversion. We identified seven gene families (WD repeat protein family, Ferritin Heavy Chain family, RAS-related Protein RAB-40 family, Diphosphoinositol polyphosphate phosphohydrolase family, Transcription Elongation Factor A family, LDOC1-related family, Zinc Finger Protein ZIC, and GLI family) that show evidence of gene conversion. Through phylogenetic analyses and synteny evidence, we show that gene conversion has played an important role in the evolution of these gene families and that gene conversion has occurred independently in both primates and rodents. Comparing the results with those of two gene conversion prediction programs (GENECONV and Partimatrix), we found that both GENECONV and Partimatrix have very high false negative rates (i.e. failed to predict gene conversions), which leads to many undetected gene conversions. The combination of phylogenetic analyses with physical synteny evidence exhibits high resolution in the detection of gene conversions.

Targeting cyclin B1 through peptide-based delivery of siRNA prevents tumour growth

2009-05-29

The development of short interfering RNA (siRNA), has provided great hope for therapeutic targeting of specific genes responsible for patholological disorders. However, the poor cellular uptake and bioavailability of siRNA remain a major obstacle to their clinical development and most strategies that propose to improve siRNA delivery remain limited for in vivo applications. In this study, we report a novel peptide-based approach, MPG-8 an improved variant of the amphipathic peptide carrier MPG, that forms nanoparticles with siRNA and promotes their efficient delivery into primary cell lines and in vivo upon intra-tumoral injection. Moreover, we show that functionalization of this carrier with cholesterol significantly improves tissue distribution and stability of siRNA in vivo, thereby enhancing the efficiency of this technology for systemic administration following intravenous injection without triggering any non-specific inflammatory response. We have validated the therapeutic potential of this strategy for cancer treatment by targeting cyclin B1 in mouse tumour models, and demonstrate that tumour growth is compromised. The robustness of the biological response achieved through this approach, infers that MPG 8-based technology holds a strong promise for therapeutic administration of siRNA.

A novel method for tissue-specific RNAi rescue in Drosophila

Schulz, J. G., David, G., Hassan, B. A. — 2009-05-29

Targeted gene silencing by RNA interference allows the study of gene function in plants and animals. In cell culture and small animal models, genetic screens can be performed—even tissue-specifically in Drosophila—with genome-wide RNAi libraries. However, a major problem with the use of RNAi approaches is the unavoidable false-positive error caused by off-target effects. Until now, this is minimized by computational RNAi design, comparing RNAi to the mutant phenotype if known, and rescue with a presumed ortholog. The ultimate proof of specificity would be to restore expression of the same gene product in vivo. Here, we present a simple and efficient method to rescue the RNAi-mediated knockdown of two independent genes in Drosophila. By exploiting the degenerate genetic code, we generated Drosophila RNAi Escape Strategy Construct (RESC) rescue proteins containing frequent silent mismatches in the complete RNAi target sequence. RESC products were no longer efficiently silenced by RNAi in cell culture and in vivo. As a proof of principle, we rescue the RNAi-induced loss of function phenotype of the eye color gene white and tracheal defects caused by the knockdown of the heparan sulfate proteoglycan syndecan. Our data suggest that RESC is widely applicable to rescue and validate ubiquitous or tissue-specific RNAi and to perform protein structure–function analysis.

Increasing the relative expression of endogenous non-coding Steroid Receptor RNA Activator (SRA) in human breast cancer cells using modified oligonucleotides

2009-05-29

Products of the Steroid Receptor RNA Activator gene (SRA1) have the unusual property to modulate the activity of steroid receptors and other transcription factors both at the RNA (SRA) and the protein (SRAP) level. Balance between these two genetically linked entities is controlled by alternative splicing of intron-1, whose retention alters SRAP reading frame. We have previously found that both fully-spliced SRAP-coding and intron-1-containing non-coding SRA RNAs co-exist in breast cancer cell lines. Herein, we report a significant (Student's t-test, P < 0.003) higher SRA–intron-1 relative expression in breast tumors with higher progesterone receptor contents. Using an antisense oligoribonucleotide, we have successfully reprogrammed endogenous SRA splicing and increased SRA RNA–intron-1 relative level in T5 breast cancer cells. This increase is paralleled by significant changes in the expression of genes such as plasminogen urokinase activator and estrogen receptor beta. Estrogen regulation of other genes, including the anti-metastatic NME1 gene, is also altered. Overall, our results suggest that the balance coding/non-coding SRA transcripts not only characterizes particular tumor phenotypes but might also, through regulating the expression of specific genes, be involved in breast tumorigenesis and tumor progression.

Light-activated RNA interference using double-stranded siRNA precursors modified using a remarkable regiospecificity of diazo-based photolabile groups

Shah, S., Jain, P. K., Kala, A., Karunakaran, D., Friedman, S. H. — 2009-05-28

Diazo-based precursors of photolabile groups have been used extensively for modifying nucleic acids, with the intention of toggling biological processes with light. These processes include transcription, translation and RNA interference. In these cases, the photolabile groups have been typically depicted as modifying the phosphate backbone of RNA and DNA. In this work we find that these diazo-based reagents in fact react very poorly with backbone phosphates. Instead, they show a remarkable specificity for terminal phosphates and very modest modification of the nucleobases. Furthermore, the photo deprotection of these terminal modifications is shown to be much more facile than nucleobase modified sites. In this study we have characterized this regiospecificity using RNA duplexes and model nucleotides, analyzed using LC/MS/MS. We have also applied this understanding of the regio-specificity to our technique of light activated RNA interference (LARI). We examined 27-mer double-stranded precursors of siRNA (‘dsRNA’), and have modified them using the photo-cleavable di-methoxy nitro phenyl ethyl group (DMNPE) group. By incorporating terminal phosphates in the dsRNA, we are able to guide DMNPE to react at these terminal locations. These modified dsRNA duplexes show superior performance to our previously described DMNPE-modified siRNA, with the range of expression that can be toggled by light increasing by a factor of two.

Specific gene silencing by artificial trans-encoded small noncoding RNAs in bacteria

Cheng, R., Miao, C., Gong, Q., Gu, Y., Lu, X., Han, F., Yu, W. — 2009-05-27

Recently, numerous small noncoding RNAs (sRNAs) with important regulatory roles have been identified in bacteria. As their eukaryotic counterparts, a major class of bacterial trans-encoded sRNAs, acts by basepairing with target mRNAs, resulting in changes in translation and stability of the mRNA. RNA interference (RNAi) has become an extraordinarily powerful RNA silencing tool for elucidating and manipulating gene functions in eukaryotes. However, such an effective RNA silencing tool remains to be developed for prokaryotes. In this study, we described firstly the use of artificial trans-encoded sRNAs (atsRNAs) for specific gene silencing in bacteria. Based on the common structural characteristics of natural bacterial trans-encoded sRNAs, we developed the designing principle of atsRNA. Most of the atsRNAs effectively suppressed the expression of exogenous EGFP gene and endogenous uidA gene in Escherichia coli. Further studies demonstrated that the mRNA base-pairing region and AU rich Hfq binding site were crucial for the activity of atsRNA. The atsRNA-mediated gene silencing was Hfq dependent. atsRNA led to translational repression and RNase-E-dependent degradation of target mRNA, and the translation inhibition was the primary event for gene silencing. Our findings demonstrated that atsRNA was an effective RNA tool for specific gene silencing in bacteria.

Efficient silencing of gene expression with modular trimeric Pol II expression cassettes comprising microRNA shuttles

Ely, A., Naidoo, T., Arbuthnot, P. — 2009-05-27

Expressed polycistronic microRNA (miR) cassettes have useful properties that can be utilized for RNA interference (RNAi)-based gene silencing. To advance their application we generated modular trimeric anti-hepatitis B virus (HBV) Pol II cassettes encoding primary (pri)-miR-31-derived shuttles that target three different viral genome sites. A panel of six expression cassettes, comprising each of the possible ordering combinations of the pri-miR-31 shuttles, was initially tested. Effective silencing of individual target sequences was achieved in transfected cells and transcribed pri-miR trimers generated intended guide strands. There was, however, variation in processing and silencing by each of the shuttles. In some cases the monomers’ position within the trimers influenced processing and this correlated with target silencing. Compromised efficacy could be compensated by substituting the pri-miR-31 backbone with a pri-miR-30a scaffold. Inhibition of HBV replication was achieved in vivo, and in cell culture without disruption of endogenous miR function or induction of the interferon response. A mutant HBV target sequence, with changes in one of the guide cognates, was also silenced by the trimeric cassettes. The modular nature of the cassettes together with compatibility with expression from Pol II promoters should be advantageous for gene silencing applications requiring simultaneous targeting of different sites.

The structure of CrgA from Neisseria meningitidis reveals a new octameric assembly state for LysR transcriptional regulators

2009-05-27

LysR-type transcriptional regulators (LTTRs) form the largest family of bacterial regulators acting as both auto-repressors and activators of target promoters, controlling operons involved in a wide variety of cellular processes. The LTTR, CrgA, from the human pathogen Neisseria meningitidis, is upregulated during bacterial–host cell contact. Here, we report the crystal structures of both regulatory domain and full-length CrgA, the first of a novel subclass of LTTRs that form octameric rings. Non-denaturing mass spectrometry analysis and analytical ultracentrifugation established that the octameric form of CrgA is the predominant species in solution in both the presence and absence of an oligonucleotide encompassing the CrgA-binding sequence. Furthermore, analysis of the isolated CrgA–DNA complex by mass spectrometry showed stabilization of a double octamer species upon DNA binding. Based on the observed structure and the mass spectrometry findings, a model is proposed in which a hexadecameric array of two CrgA oligomers binds to its DNA target site.

High-affinity triplex targeting of double stranded DNA using chemically modified peptide nucleic acid oligomers

Hansen, M. E., Bentin, T., Nielsen, P. E. — 2009-05-27

While sequence-selective dsDNA targeting by triplex forming oligonucleotides has been studied extensively, only very little is known about the properties of PNA–dsDNA triplexes—mainly due to the competing invasion process. Here we show that when appropriately modified using pseudoisocytosine substitution, in combination with (oligo)lysine or 9-aminoacridine conjugation, homopyrimidine PNA oligomers bind complementary dsDNA targets via triplex formation with (sub)nanomolar affinities (at pH 7.2, 150 mM Na⁺). Binding affinity can be modulated more than 1000-fold by changes in pH, PNA oligomer length, PNA net charge and/or by substitution of pseudoisocytosine for cytosine, and conjugation of the DNA intercalator 9-aminoacridine. Furthermore, 9-aminoacridine conjugation also strongly enhanced triplex invasion. Specificity for the fully matched target versus one containing single centrally located mismatches was more than 150-fold. Together the data support the use of homopyrimidine PNAs as efficient and sequence selective tools in triplex targeting strategies under physiological relevant conditions.

Using high-density exon arrays to profile gene expression in closely related species

Lin, L., Liu, S., Brockway, H., Seok, J., Jiang, P., Wong, W. H., Xing, Y. — 2009-05-27

Global comparisons of gene expression profiles between species provide significant insight into gene regulation, evolutionary processes and disease mechanisms. In this work, we describe a flexible and intuitive approach for global expression profiling of closely related species, using high-density exon arrays designed for a single reference genome. The high-density probe coverage of exon arrays allows us to select identical sets of perfect-match probes to measure expression levels of orthologous genes. This eliminates a serious confounding factor in probe affinity effects of species-specific microarray probes, and enables direct comparisons of estimated expression indexes across species. Using a newly designed Affymetrix exon array, with eight probes per exon for approximately 315 000 exons in the human genome, we conducted expression profiling in corresponding tissues from humans, chimpanzees and rhesus macaques. Quantitative real-time PCR analysis of differentially expressed candidate genes is highly concordant with microarray data, yielding a validation rate of 21/22 for human versus chimpanzee differences, and 11/11 for human versus rhesus differences. This method has the potential to greatly facilitate biomedical and evolutionary studies of gene expression in nonhuman primates and can be easily extended to expression array design and comparative analysis of other animals and plants.

Construction and functional analyses of a comprehensive {sigma}54 site-directed mutant library using alanine-cysteine mutagenesis

Xiao, Y., Wigneshweraraj, S. R., Weinzierl, R., Wang, Y.-P., Buck, M. — 2009-05-27

The ⁵⁴ factor associates with core RNA polymerase (RNAP) to form a holoenzyme that is unable to initiate transcription unless acted on by an activator protein. ⁵⁴ is closely involved in many steps of activator-dependent transcription, such as core RNAP binding, promoter recognition, activator interaction and open complex formation. To systematically define ⁵⁴ residues that contribute to each of these functions and to generate a resource for site specific protein labeling, a complete mutant library of ⁵⁴ was constructed by alanine–cysteine scanning mutagenesis. Amino acid residues from 3 to 476 of Cys(-)⁵⁴ were systematically mutated to alanine and cysteine in groups of two adjacent residues at a time. The influences of each substitution pair upon the functions of ⁵⁴ were analyzed in vivo and in vitro and the functions of many residues were revealed for the first time. Increased ⁵⁴ isomerization activity seldom corresponded with an increased transcription activity of the holoenzyme, suggesting the steps after ⁵⁴ isomerization, likely to be changes in core RNAP structure, are also strictly regulated or rate limiting to open complex formation. A linkage between core RNAP-binding activity and activator responsiveness indicates that the ⁵⁴-core RNAP interface changes upon activation.

Methylation detection oligonucleotide microarray analysis: a high-resolution method for detection of CpG island methylation

2009-05-27

Methylation of CpG islands associated with genes can affect the expression of the proximal gene, and methylation of non-associated CpG islands correlates to genomic instability. This epigenetic modification has been shown to be important in many pathologies, from development and disease to cancer. We report the development of a novel high-resolution microarray that detects the methylation status of over 25 000 CpG islands in the human genome. Experiments were performed to demonstrate low system noise in the methodology and that the array probes have a high signal to noise ratio. Methylation measurements between different cell lines were validated demonstrating the accuracy of measurement. We then identified alterations in CpG islands, both those associated with gene promoters, as well as non-promoter-associated islands in a set of breast and ovarian tumors. We demonstrate that this methodology accurately identifies methylation profiles in cancer and in principle it can differentiate any CpG methylation alterations and can be adapted to analyze other species.

Directed evolution of an orthogonal nucleoside analog kinase via fluorescence-activated cell sorting

Liu, L., Li, Y., Liotta, D., Lutz, S. — 2009-05-27

Nucleoside analogs (NAs) represent an important category of prodrugs for the treatment of viral infections and cancer, yet the biological potency of many analogs is compromised by their inefficient activation through cellular 2'-deoxyribonucleoside kinases (dNKs). We herein report the directed evolution and characterization of an orthogonal NA kinase for 3'-deoxythymidine (ddT), using a new FACS-based screening protocol in combination with a fluorescent analog of ddT. Four rounds of random mutagenesis and DNA shuffling of Drosophila melanogaster 2'-deoxynucleoside kinase, followed by FACS analysis, yielded an orthogonal ddT kinase with a 6-fold higher activity for the NA and a 20-fold k_cat/K_M preference for ddT over thymidine, an overall 10 000-fold change in substrate specificity. The contributions of individual amino acid substitutions in the ddT kinase were evaluated by reverse engineering, enabling a detailed structure–function analysis to rationalize the observed changes in performance. Based on our results, kinase engineering with fluorescent NAs and FACS should prove a highly versatile method for evolving selective kinase:NA pairs and for studying fundamental aspects of the structure–function relationship in dNKs.

A conformational change in the helicase core is necessary but not sufficient for RNA unwinding by the DEAD box helicase YxiN

Karow, A. R., Klostermeier, D. — 2009-05-27

Cooperative binding of ATP and RNA to DEAD-box helicases induces the closed conformation of their helicase core, with extensive interactions across the domain interface. The bound RNA is bent, and its distortion may constitute the first step towards RNA unwinding. To dissect the role of the conformational change in the helicase core for RNA unwinding, we characterized the RNA-stimulated ATPase activity, RNA unwinding and the propensity to form the closed conformer for mutants of the DEAD box helicase YxiN. The ATPase-deficient K52Q mutant forms a closed conformer upon binding of ATP and RNA, but is deficient in RNA unwinding. A mutation in motif III slows down the catalytic cycle, but neither affects the propensity for the closed conformer nor its global conformation. Hence, the closure of the cleft in the helicase core is necessary but not sufficient for RNA unwinding. In contrast, the G303A mutation in motif V prevents a complete closure of the inter-domain cleft, affecting ATP binding and hydrolysis and is detrimental to unwinding. Possibly, the K52Q and motif III mutants still introduce a kink into the backbone of bound RNA, whereas G303A fails to kink the RNA substrate.

Physical and functional interactions between Escherichia coli MutL and the Vsr repair endonuclease

2009-05-27

DNA mismatch repair (MMR) and very-short patch (VSP) repair are two pathways involved in the repair of T:G mismatches. To learn about competition and cooperation between these two repair pathways, we analyzed the physical and functional interaction between MutL and Vsr using biophysical and biochemical methods. Analytical ultracentrifugation reveals a nucleotide-dependent interaction between Vsr and the N-terminal domain of MutL. Using chemical crosslinking, we mapped the interaction site of MutL for Vsr to a region between the N-terminal domains similar to that described before for the interaction between MutL and the strand discrimination endonuclease MutH of the MMR system. Competition between MutH and Vsr for binding to MutL resulted in inhibition of the mismatch-provoked MutS- and MutL-dependent activation of MutH, which explains the mutagenic effect of Vsr overexpression. Cooperation between MMR and VSP repair was demonstrated by the stimulation of the Vsr endonuclease in a MutS-, MutL- and ATP-hydrolysis-dependent manner, in agreement with the enhancement of VSP repair by MutS and MutL in vivo. These data suggest a mobile MutS–MutL complex in MMR signalling, that leaves the DNA mismatch prior to, or at the time of, activation of downstream effector molecules such as Vsr or MutH.

A biochemically active MCM-like helicase in Bacillus cereus

2009-05-27

The mini-chromosome maintenance (MCM) proteins serve as the replicative helicases in archaea and eukaryotes. Interestingly, an MCM homolog was identified, by BLAST analysis, within a phage integrated in the bacterium Bacillus cereus (Bc). BcMCM is only related to the AAA region of MCM-helicases; the typical amino-terminus is missing and is replaced by a segment with weak homology to primases. We show that BcMCM displays 3'->5' helicase and ssDNA-stimulated ATPase activity, properties that arise from its conserved AAA domain. Isolated BcMCM is a monomer in solution but likely forms the functional oligomer in vivo. We found that the BcMCM amino-terminus can bind ssDNA and harbors a zinc atom, both hallmarks of the typical MCM amino-terminus. No BcMCM-catalyzed primase activity could be detected. We propose that the divergent amino-terminus of BcMCM is a paralog of the corresponding region of MCM-helicases. A divergent amino terminus makes BcMCM a useful model for typical MCM-helicases since it accomplishes the same function using an apparently unrelated structure.

Processing of thymine glycol in a clustered DNA damage site: mutagenic or cytotoxic

Bellon, S., Shikazono, N., Cunniffe, S., Lomax, M., O'Neill, P. — 2009-05-25

Localized clustering of damage is a hallmark of certain DNA-damaging agents, particularly ionizing radiation. The potential for genetic change arising from the effects of clustered damage sites containing combinations of AP sites, 8-oxo-7,8-dihydroguanine (8-oxoG) or 5,6-dihydrothymine is high. To date clusters containing a DNA base lesion that is a strong block to replicative polymerases, have not been explored. Since thymine glycol (Tg) is non-mutagenic but a strong block to replicative polymerases, we have investigated whether clusters containing Tg are highly mutagenic or lead to potentially cytotoxic lesions, when closely opposed to either 8-oxoG or an AP site. Using a bacterial plasmid-based assay and repair assays using cell extracts or purified proteins, we have shown that DNA double-strand breaks (DSBs) arise when Tg is opposite to an AP site, either through attempted base excision repair or at replication. In contrast, 8-oxoG opposite to Tg in a cluster ‘protects’ against DSB formation but does enhance the mutation frequency at the site of 8-oxoG relative to that at a single 8-oxoG, due to the decisive role of endonucleases in the initial stages of processing Tg/8-oxoG clusters, removing Tg to give an intermediate with an abasic site or single-strand break.

Mismatch repair and nucleotide excision repair proteins cooperate in the recognition of DNA interstrand crosslinks

Zhao, J., Jain, A., Iyer, R. R., Modrich, P. L., Vasquez, K. M. — 2009-05-25

DNA interstrand crosslinks (ICLs) are among the most cytotoxic types of DNA damage, thus ICL-inducing agents such as psoralen, are clinically useful chemotherapeutics. Psoralen-modified triplex-forming oligonucleotides (TFOs) have been used to target ICLs to specific genomic sites to increase the selectivity of these agents. However, how TFO-directed psoralen ICLs (Tdp-ICLs) are recognized and processed in human cells is unclear. Previously, we reported that two essential nucleotide excision repair (NER) protein complexes, XPA–RPA and XPC–RAD23B, recognized ICLs in vitro, and that cells deficient in the DNA mismatch repair (MMR) complex MutSβ were sensitive to psoralen ICLs. To further investigate the role of MutSβ in ICL repair and the potential interaction between proteins from the MMR and NER pathways on these lesions, we performed electrophoretic mobility-shift assays and chromatin immunoprecipitation analysis of MutSβ and NER proteins with Tdp-ICLs. We found that MutSβ bound to Tdp-ICLs with high affinity and specificity in vitro and in vivo, and that MutSβ interacted with XPA–RPA or XPC–RAD23B in recognizing Tdp-ICLs. These data suggest that proteins from the MMR and NER pathways interact in the recognition of ICLs, and provide a mechanistic link by which proteins from multiple repair pathways contribute to ICL repair.

The DNA-recognition mode shared by archaeal feast/famine-regulatory proteins revealed by the DNA-binding specificities of TvFL3, FL10, FL11 and Ss-LrpB

2009-05-25

The DNA-binding mode of archaeal feast/famine-regulatory proteins (FFRPs), i.e. paralogs of the Esherichia coli leucine-responsive regulatory protein (Lrp), was studied. Using the method of systematic evolution of ligands by exponential enrichment (SELEX), optimal DNA duplexes for interacting with TvFL3, FL10, FL11 and Ss-LrpB were identified as TACGA[AAT/ATT]TCGTA, GTTCGA[AAT/ATT]TCGAAC, CCGAAA[AAT/ATT]TTTCGG and TTGCAA[AAT/ATT]TTGCAA, respectively, all fitting into the form abcdeWWWedcba. Here W is A or T, and e.g. a and a are bases complementary to each other. Apparent equilibrium binding constants of the FFRPs and various DNA duplexes were determined, thereby confirming the DNA-binding specificities of the FFRPs. It is likely that these FFRPs recognize DNA in essentially the same way, since their DNA-binding specificities were all explained by the same pattern of relationship between amino-acid positions and base positions to form chemical interactions. As predicted from this relationship, when Gly36 of TvFL3 was replaced by Thr, the b base in the optimal DNA duplex changed from A to T, and, when Thr36 of FL10 was replaced by Ser, the b base changed from T to G/A. DNA-binding characteristics of other archaeal FFRPs, Ptr1, Ptr2, Ss-Lrp and LysM, are also consistent with the relationship.

A complex-centric view of protein network evolution

Yosef, N., Kupiec, M., Ruppin, E., Sharan, R. — 2009-05-22

The recent availability of protein–protein interaction networks for several species makes it possible to study protein complexes in an evolutionary context. In this article, we present a novel network-based framework for reconstructing the evolutionary history of protein complexes. Our analysis is based on generalizing evolutionary measures for single proteins to the level of whole subnetworks, comprehensively considering a broad set of computationally derived complexes and accounting for both sequence and interaction changes. Specifically, we compute sets of orthologous complexes across species, and use these to derive evolutionary rate and age measures for protein complexes. We observe significant correlations between the evolutionary properties of a complex and those of its member proteins, suggesting that protein complexes form early in evolution and evolve as coherent units. Additionally, our approach enables us to directly quantify the extent to which gene duplication has played a role in the evolution of complexes. We find that about one quarter of the sets of orthologous complexes have originated from evolutionary cores of homodimers that underwent duplication and divergence, testifying to the important role of gene duplication in protein complex evolution.

Modulation of alternative splicing by long-range RNA structures in Drosophila

Raker, V. A., Mironov, A. A., Gelfand, M. S., Pervouchine, D. D. — 2009-05-22

Accurate and efficient recognition of splice sites during pre-mRNA splicing is essential for proper transcriptome expression. Splice site usage can be modulated by secondary structures, but it is unclear if this type of modulation is commonly used or occurs to a significant degree with secondary structures forming over long distances. Using phlyogenetic comparisons of intronic sequences among 12 Drosophila genomes, we elucidated a group of 202 highly conserved pairs of sequences, each at least nine nucleotides long, capable of forming stable stem structures. This set was highly enriched in alternatively spliced introns and introns with weak acceptor sites and long introns, and most occurred over long distances (>150 nucleotides). Experimentally, we analyzed the splicing of several of these introns using mini-genes in Drosophila S2 cells. Wild-type splicing patterns were changed by mutations that opened the stem structure, and restored by compensatory mutations that re-established the base-pairing potential, demonstrating that these secondary structures were indeed implicated in the splice site choice. Mechanistically, the RNA structures masked splice sites, brought together distant splice sites and/or looped out introns. Thus, base-pairing interactions within introns, even those occurring over long distances, are more frequent modulators of alternative splicing than is currently assumed.

RFXB and its splice variant RFXBSV mediate the antagonism between IFN{gamma} and TGF{beta} on COL1A2 transcription in vascular smooth muscle cells

Fang, M., Kong, X., Li, P., Fang, F., Wu, X., Bai, H., Qi, X., Chen, Q., Xu, Y. — 2009-05-22

Cytokines secreted by infiltrating immune cells during atherogenesis modulate vascular remodeling. One exemplary event is the antagonism between transformed growth factor (TGF-β) and interferon gamma (IFN-) on the transcriptional control of type I collagen gene (COL1A2). Previously we have reported that IFN- up-regulates regulatory factor for X-box B (RFXB) to repress collagen transcription while down-regulates the expression of RFXBSV, a splice variant of RFXB that blocks collagen repression in fibroblasts. Here we demonstrate that TGF-β abrogated COL1A2 repression by IFN- through altering the relative expression of RFXB and RFXBSV. Unlike RFXB, RFXBSV did not bind to the collagen promoter and competed with RFXB for the co-repressor histone deacetylase 2 (HDAC2), limiting HDAC2 recruitment to the collagen transcription start site as evidenced by chromatin immunoprecipitation assays. Over-expression of RFXB by lentiviral infection in HASMCs enhanced HDAC2 enlistment, promoted histone deacetylation surrounding the collagen site by IFN-, and blocked the TGF-β antagonism, a pattern reversed by RFXBSV infection. On the contrary, silencing of RFXB, but not both RFXB and RFXBSV, expression promoted the TGF-β antagonism. Thus, we have identified a novel mechanism whereby TGF-β antagonizes the IFN- repression of collagen transcription in HASMCs and as such provided new insights into antiatherogenic strategies.

In silico selection of RNA aptamers

Chushak, Y., Stone, M. O. — 2009-05-21

In vitro selection of RNA aptamers that bind to a specific ligand usually begins with a random pool of RNA sequences. We propose a computational approach for designing a starting pool of RNA sequences for the selection of RNA aptamers for specific analyte binding. Our approach consists of three steps: (i) selection of RNA sequences based on their secondary structure, (ii) generating a library of three-dimensional (3D) structures of RNA molecules and (iii) high-throughput virtual screening of this library to select aptamers with binding affinity to a desired small molecule. We developed a set of criteria that allows one to select a sequence with potential binding affinity from a pool of random sequences and developed a protocol for RNA 3D structure prediction. As verification, we tested the performance of in silico selection on a set of six known aptamer–ligand complexes. The structures of the native sequences for the ligands in the testing set were among the top 5% of the selected structures. The proposed approach reduces the RNA sequences search space by four to five orders of magnitude—significantly accelerating the experimental screening and selection of high-affinity aptamers.

The role of DNA damage response pathways in chromosome fragility in Fragile X syndrome

Kumari, D., Somma, V., Nakamura, A. J., Bonner, W. M., D'Ambrosio, E., Usdin, K. — 2009-05-21

FRAXA is one of a number of fragile sites in human chromosomes that are induced by agents like fluorodeoxyuridine (FdU) that affect intracellular thymidylate levels. FRAXA coincides with a >200 CGG•CCG repeat tract in the 5' UTR of the FMR1 gene, and alleles prone to fragility are associated with Fragile X (FX) syndrome, one of the leading genetic causes of intellectual disability. Using siRNA depletion, we show that ATR is involved in protecting the genome against FdU-induced chromosome fragility. We also show that FdU increases the number of -H2AX foci seen in both normal and patient cells and increases the frequency with which the FMR1 gene colocalizes with these foci in patient cells. In the presence of FdU and KU55933, an ATM inhibitor, the incidence of chromosome fragility is reduced, suggesting that ATM contributes to FdU-induced chromosome fragility. Since both ATR and ATM are involved in preventing aphidicolin-sensitive fragile sites, our data suggest that the lesions responsible for aphidicolin-induced and FdU-induced fragile sites differ. FRAXA also displays a second form of chromosome fragility in absence of FdU, which our data suggest is normally prevented by an ATM-dependent process.

HMGA2 exhibits dRP/AP site cleavage activity and protects cancer cells from DNA-damage-induced cytotoxicity during chemotherapy

2009-05-21

HMGA proteins are not translated in normal human somatic cells, but are present in high copy numbers in pluripotent embryonic stem cells and most neoplasias. Correlations between the degree of malignancy, patient prognostic index and HMGA levels have been firmly established. Intriguingly, HMGA2 is also found in rare tumor-inducing cells which are resistant to chemotherapy. Here, we demonstrate that HMGA1a/b and HMGA2 possess intrinsic dRP and AP site cleavage activities, and that lysines and arginines in the AT-hook DNA-binding domains function as nucleophiles. We also show that HMGA2 can be covalently trapped at genomic abasic sites in cancer cells. By employing a variety of cell-based assays, we provide evidence that the associated lyase activities promote cellular resistance against DNA damage that is targeted by base excision repair (BER) pathways, and that this protection directly correlates with the level of HMGA2 expression. In addition, we demonstrate an interaction between human AP endonuclease 1 and HMGA2 in cancer cells, which supports our conclusion that HMGA2 can be incorporated into the cellular BER machinery. Our study thus identifies an unexpected role for HMGA2 in DNA repair in cancer cells which has important clinical implications for disease diagnosis and therapy.

The Walker B motif in avian FANCM is required to limit sister chromatid exchanges but is dispensable for DNA crosslink repair

Rosado, I. V., Niedzwiedz, W., Alpi, A. F., Patel, K. J. — 2009-05-21

FANCM, the most highly conserved component of the Fanconi Anaemia (FA) pathway can resolve recombination intermediates and remodel synthetic replication forks. However, it is not known if these activities are relevant to how this conserved protein activates the FA pathway and promotes DNA crosslink repair. Here we use chicken DT40 cells to systematically dissect the function of the helicase and nuclease domains of FANCM. Our studies reveal that these domains contribute distinct roles in the tolerance of crosslinker, UV light and camptothecin-induced DNA damage. Although the complete helicase domain is critical for crosslink repair, a predicted inactivating mutation of the Walker B box domain has no impact on FA pathway associated functions. However, this mutation does result in elevated sister chromatid exchanges (SCE). Furthermore, our genetic dissection indicates that FANCM functions with the Blm helicase to suppress spontaneous SCE events. Overall our results lead us to reappraise the role of helicase domain associated activities of FANCM with respect to the activation of the FA pathway, crosslink repair and in the resolution of recombination intermediates.

Base extrusion is found at helical junctions between right- and left-handed forms of DNA and RNA

Kim, D., Reddy, S., Kim, D. Y., Rich, A., Lee, S., Kim, K. K., Kim, Y.-G. — 2009-05-21

Base extrusion is a major structural feature at the junction between B- and Z-DNA (the B–Z junction) where a base pair is broken, and the two bases are extruded from the double helix. Despite the demonstration of base extrusion at the B–Z junction, it is not clear whether a similar base extrusion occurs at other types of junctions involving the left-handed Z conformation. Here, we investigate structural changes of bases at three Z-form junctions: DNA B–Z and Z–Z and RNA A–Z junctions. By monitoring fluorescently labeled duplex nucleic acids using 2-aminopurines at various positions relative to the junction point, we show that base extrusion occurs not only at the DNA B–Z junction, but also at the RNA A–Z and DNA Z–Z junctions. Our data suggest that base extrusion is a general feature of Z-form nucleic-acid junctions.

Ets-1 p51 and p42 isoforms differentially modulate Stromelysin-1 promoter according to induced DNA bend orientation

Leprivier, G., Baillat, D., Begue, A., Hartmann, B., Aumercier, M. — 2009-05-21

The Stromelysin-1 gene promoter contains a palindrome of two Ets-binding sites (EBS) that bind the p51 and p42 isoforms of the human Ets-1-transcription factor. A previous study established that full gene transactivation is associated with a ternary complex consisting of two p51 bound to the two EBS on the promoter. p42, only able to bind one of the two EBS, induces only very weak activity. Here, we investigate the mechanism by which the Stromelysin-1 promoter discriminates between p51 and p42. The differential stoichiometry of the two Ets-1 isoforms arises from the Stromelysin-1 EBS palindrome. The ternary complex requires the presence of two inhibitory domains flanking the DNA-binding domain and the ability to form an intramolecular autoinhibition module. Most importantly, the p51-ternary and the p42-binary complexes induce DNA curvatures with opposite orientations. These results establish that differential DNA bending, via p51 and p42 differential binding, is correlated with the Stromelysin-1 promoter activation process.

High-throughput sequence-based epigenomic analysis of Alu repeats in human cerebellum

2009-05-20

DNA methylation, the only known covalent modification of mammalian DNA, occurs primarily in CpG dinucleotides. 51% of CpGs in the human genome reside within repeats, and 25% within Alu elements. Despite that, no method has been reported for large-scale ascertainment of CpG methylation in repeats. Here we describe a sequencing-based strategy for parallel determination of the CpG-methylation status of thousands of Alu repeats, and a computation algorithm to design primers that enable their specific amplification from bisulfite converted genomic DNA. Using a single primer pair, we generated amplicons of high sequence complexity, and derived CpG-methylation data from 31 178 Alu elements and their 5' flanking sequences, altogether representing over 4 Mb of a human cerebellum epigenome. The analysis of the Alu methylome revealed that the methylation level of Alu elements is high in the intronic and intergenic regions, but low in the regions close to transcription start sites. Several hypomethylated Alu elements were identified and their hypomethylated status verified by pyrosequencing. Interestingly, some Alu elements exhibited a strikingly tissue-specific pattern of methylation. We anticipate the amplicons herein described to prove invaluable as epigenome representations, to monitor epigenomic alterations during normal development, in aging and in diseases such as cancer.

One RNA aptamer sequence, two structures: a collaborating pair that inhibits AMPA receptors

Huang, Z., Pei, W., Han, Y., Jayaseelan, S., Shekhtman, A., Shi, H., Niu, L. — 2009-05-20

RNA is ideally suited for in vitro evolution experiments, because a single RNA molecule possesses both genotypic (replicable sequence) and phenotypic (selectable shape) properties. Using systematic evolution of ligands by exponential enrichment (SELEX), we found a single 58-nt aptamer sequence that assumes two structures with different functions, both of which are required to inhibit the GluR2 AMPA receptor channel. Yet, the two structures, once formed during transcription, appear to be incapable of interconverting through unfolding and refolding, presumably due to their extraordinary structural stability. Thus, our results suggest more broadly that natural RNA molecules can evolve to acquire alternative structures and associated functions. Such divergence of RNA phenotype may precede gene duplication at the genome level.

Quantitative evaluation of DNA methylation by optimization of a differential-high resolution melt analysis protocol

Malentacchi, F., Forni, G., Vinci, S., Orlando, C. — 2009-05-19

DNA methylation is a key regulator of gene transcription. Alterations in DNA methylation patterns are common in most cancers, occur early in carcinogenesis and can be detected in body fluids. Reliable and sensitive quantitative assays are required to improve the diagnostic role of methylation in the management of cancer patients. Here we present an optimized procedure, based on differential-high resolution melting analysis (D-HRMA), for the rapid and accurate quantification of methylated DNA. Two sets of primers are used in a single tube for the simultaneous amplification of the methylated (M) and unmethylated (Um) DNA sequences in D-HRMA. After HRM, differential fluorescence was calculated at the specific melting temperature after automatic subtraction of UM-DNA fluorescence. Quantification was calculated by interpolation on an external standard curve generated by serial dilutions of M-DNA. To optimize the protocol, nine primer sets were accurately selected on the basis of the number of CpG on promoters of hTERT and Bcl2 genes. The use of optimized D-HRMA allowed us to detect up to 0.025% M-DNA. D-HRMA results of DNA from 85 bladder cancers were comparable to those obtained with real time quantitative methylation specific PCR. In addition, D-HRMA appears suitable for rapid and efficient measurements in ‘in vitro’ experiments on methylation patterns after treatment with demethylating drugs.

Specific sequences within arginine-glycine-rich domains affect mRNA-binding protein function

McBride, A. E., Conboy, A. K., Brown, S. P., Ariyachet, C., Rutledge, K. L. — 2009-05-19

The discovery of roles for arginine methylation in intracellular transport and mRNA splicing has focused attention on the methylated arginine–glycine (RG)-rich domains found in many eukaryotic RNA-binding proteins. Sequence similarity among these highly repetitive RG domains, combined with interactions between RG-rich proteins, raises the question of whether these regions are general interaction motifs or whether there is specificity within these domains. Using the essential Saccharomyces cerevisiae mRNA-binding protein Npl3 (ScNpl3) as a model system, we first tested the importance of the RG domain for protein function. While Npl3 lacking the RG domain could not support growth of cells lacking Npl3, surprisingly, expression of the RG domain alone supported partial growth of these cells. To address the specificity of this domain, we created chimeric forms of ScNpl3 with RG-rich domains of S. cerevisiae nucleolar proteins, Gar1 and Nop1 (ScGar1, ScNop1), or of the Candida albicans Npl3 ortholog (CaNpl3). Whereas the CaNpl3 RG chimeric protein retained nearly wild-type function in S. cerevisiae, the ScGar1 and ScNop1 RG domains significantly reduced Npl3 function and self-association, indicating RG domain specificity. Nuclear localization of Npl3 also requires specific RG sequences, yet heterologous RG domains allow similar modulation of Npl3 transport by arginine methylation.

R-loops do not accumulate in transcription-defective hpr1-101 mutants: implications for the functional role of THO/TREX

Gomez-Gonzalez, B., Aguilera, A. — 2009-05-18

To get further insight into the effect that THO/TREX and R-loops have in transcription-associated recombination and transcription, we analyzed the ability to form R-loops of hpr1-101, a THO mutation that impairs transcription and mRNP biogenesis without triggering hyper-recombination. Human AID, a cytidine deaminase that acts on ssDNA displaced by RNA-DNA hybrids, strongly induced both hyper-recombination and hyper-mutation in hpr1-101, similar to hpr1 mutants. However, in contrast to hpr1, AID-induced mutations in hpr1-101 occur at similar frequencies in both the transcribed and non-transcribed strands, implying that the enhanced AID action in these mutants is not caused by co-transcriptional R-loops. These results indicate for the first time that THO has a transcriptional function that is not mediated by R-loops, providing a new perspective for the understanding of the coupling of transcription with mRNP biogenesis and export.

Identification of candidate regulatory SNPs by combination of transcription-factor-binding site prediction, SNP genotyping and haploChIP

Ameur, A., Rada-Iglesias, A., Komorowski, J., Wadelius, C. — 2009-05-18

Disease-associated SNPs detected in large-scale association studies are frequently located in non-coding genomic regions, suggesting that they may be involved in transcriptional regulation. Here we describe a new strategy for detecting regulatory SNPs (rSNPs), by combining computational and experimental approaches. Whole genome ChIP-chip data for USF1 was analyzed using a novel motif finding algorithm called BCRANK. 1754 binding sites were identified and 140 candidate rSNPs were found in the predicted sites. For validating their regulatory function, seven SNPs found to be heterozygous in at least one of four human cell samples were investigated by ChIP and sequence analysis (haploChIP). In four of five cases where the SNP was predicted to affect binding, USF1 was preferentially bound to the allele containing the consensus motif. Allelic differences in binding for other proteins and histone marks further reinforced the SNPs regulatory potential. Moreover, for one of these SNPs, H3K36me3 and POLR2A levels at neighboring heterozygous SNPs indicated effects on transcription. Our strategy, which is entirely based on in vivo data for both the prediction and validation steps, can identify individual binding sites at base pair resolution and predict rSNPs. Overall, this approach can help to pinpoint the causative SNPs in complex disorders where the associated haplotypes are located in regulatory regions. Availability: BCRANK is available from Bioconductor (http://www.bioconductor.org/).

Detection of intergenic non-coding RNAs expressed in the main developmental stages in Drosophila melanogaster

2009-05-18

How many intergenically encoded non-coding RNAs (ncRNAs) are expressed during various developmental stages in Drosophila? Previous analyses in one or a few developmental stages indicated abundant expression of intergenic ncRNAs. However, some reported that ncRNAs have been recently falsified, and, in general, the false positive rate for ncRNA detection is unknown. In this report, we used reverse transcription-PCR (RT-PCR), a more robust method, to detect ncRNAs from the intergenic regions that are expressed in four major developmental stages (6–8 h embryo, 20–22 h embryo, larvae and adult). We tested 1027 regions, ~10% of all intergenic regions, and detected transcription by RT–PCR. The results from 18 342 RT–PCR experiments revealed evidence for transcription in 72.7% of intergenic regions in the developmental process. The early developmental stage appears to be associated with more abundant ncRNAs than later developmental stages. In the early stage, we detected 43.6% of intergenic regions that encode transcripts in the triplicate RT–PCR experiments, yielding an estimate of 5006 intergenic regions in the entire genome likely encoding ncRNAs. We compared the RT–PCR-related approach with previous tiling array-based approach and observed that the latter method is insensitive to short ncRNAs, especially the molecules less than 120 bp. We measured false positive rates for the analyzed genomic approaches including the RT–PCR and tiling array method.

Human PSF binds to RAD51 and modulates its homologous-pairing and strand-exchange activities

Morozumi, Y., Takizawa, Y., Takaku, M., Kurumizaka, H. — 2009-05-15

RAD51, a eukaryotic recombinase, catalyzes homologous-pairing and strand-exchange reactions, which are essential steps in homologous recombination and recombinational repair of double strand breaks. On the other hand, human PSF was originally identified as a component of spliceosomes, and its multiple functions in RNA processing, transcription and DNA recombination were subsequently revealed. In the present study, we found that PSF directly interacted with RAD51. PSF significantly enhanced RAD51-mediated homologous pairing and strand exchange at low RAD51 concentrations; however, in contrast, it inhibited these RAD51-mediated recombination reactions at the optimal RAD51 concentration. Deletion analyses revealed that the N-terminal region of PSF possessed the RAD51- and DNA-binding activities, but the central region containing the RNA-recognition motifs bound neither RAD51 nor DNA. These results suggest that PSF may have dual functions in homologous recombination and RNA processing through its N-terminal and central regions, respectively.

Inhibition of translesion DNA polymerase by archaeal reverse gyrase

Valenti, A., Perugino, G., Nohmi, T., Rossi, M., Ciaramella, M. — 2009-05-14

Reverse gyrase is a unique DNA topoisomerase endowed with ATP-dependent positive supercoiling activity. It is typical of microorganisms living at high temperature and might play a role in maintenance of genome stability and repair. We have identified the translesion DNA polymerase SsoPolY/Dpo4 as one partner of reverse gyrase in the hyperthermophilic archaeon Sulfolobus solfataricus. We show here that in cell extracts, PolY and reverse gyrase co-immunoprecipitate with each other and with the single strand binding protein, SSB. The interaction is confirmed in vitro by far-western and Surface Plasmon Resonance. In functional assays, reverse gyrase inhibits PolY, but not the S. solfataricus B-family DNA polymerase PolB1. Mutational analysis shows that inhibition of PolY activity depends on both ATPase and topoisomerase activities of reverse gyrase, suggesting that the intact positive supercoiling activity is required for PolY inhibition. In vivo, reverse gyrase and PolY are degraded after induction of DNA damage. Inhibition by reverse gyrase and degradation might act as a double mechanism to control PolY and prevent its potentially mutagenic activity when undesired. Inhibition of a translesion polymerase by topoisomerase-induced modification of DNA structure may represent a previously unconsidered mechanism of regulation of these two-faced enzymes.

p53-dependent antiviral RNA-interference facilitates tumor-selective viral replication

2009-05-14

RNA-interference (RNAi) is a potent tool for specific gene silencing. In this study, we developed an adenovirus for conditional replication in p53-dysfunctional tumor cells that uses p53-selective expression of a microRNA-network directed against essential adenoviral genes. Compared to a control virus that expressed a scrambled microRNA-network, antiviral RNAi selectively attenuated viral replication in cells with transcriptionally active p53, but not in p53-dysfunctional tumor cells where both viruses replicated equivalently. Since these results were confirmed by an in vivo comparison of both viruses after infection of p53-knockout and normal mice, we could demonstrate that attenuated replication was indeed a result of p53-selective exhibition of antiviral RNAi. Addressing the therapeutic applicability, we could show that the application of RNAi-controlled virus efficiently lysed p53-dysfunctional tumors in vitro and in vivo but resulted in drastically reduced load of virus-DNA in the liver of treated mice. We have generated a broadly applicable adenovirus for selective destruction of p53-dysfunctional tumors and thereby demonstrate that virus-encoded RNAi-networks represent an efficient and versatile tool to modify viral functions. RNAi-networks can be applied to all transcriptionally regulated DNA-viruses to remodulate viral tropism and thus provide means to generate specifically replicating vectors for clinical applications.

CUGBP2 directly interacts with U2 17S snRNP components and promotes U2 snRNA binding to cardiac troponin T pre-mRNA

Goo, Y.-H., Cooper, T. A. — 2009-05-14

CUGBP2 (ETR-3/NAPOR/BRUNOL3) promotes inclusion of cardiac troponin T (cTNT) exon 5 via binding between positions 21 and 74 of the downstream intron. The molecular mechanism by which CUGBP2 activates cTNT exon 5 inclusion is unknown. Our results suggest that CUGBP2 promotes exon inclusion by a novel mechanism in which CUGBP2 directly interacts with components of the activated U2 snRNP and enhances binding of U2 snRNP to the branch site located upstream of the exon. Using an in vitro splicing assay, we show that recombinant CUGBP2 enhances complex A formation of a cTNT pre-mRNA. Enhanced complex A assembly requires both the upstream and downstream introns consistent with dual requirements for the downstream CUGBP2-binding site and an upstream branch site for U2 snRNP binding. We also show that CUGBP2 enhances binding of U2 snRNA to the cTNT pre-mRNA consistent with enhanced complex A assembly. Purification of CUGBP2-interacting proteins using tandem affinity purification leads to the demonstration that the core 17S U2 snRNP components, SF3b145 and SF3b49 bind directly to CUGBP2. We conclude that CUGBP2 activates exon inclusion by forming direct interactions with components of the 17S snRNP complex and recruits and/or stabilizes binding of U2 snRNP.

ROS1 5-methylcytosine DNA glycosylase is a slow-turnover catalyst that initiates DNA demethylation in a distributive fashion

Ponferrada-Marin, M. I., Roldan-Arjona, T., Ariza, R. R. — 2009-05-13

Arabidopsis ROS1 belongs to a family of plant 5-methycytosine DNA glycosylases that initiate DNA demethylation through base excision. ROS1 displays the remarkable capacity to excise 5-meC, and to a lesser extent T, while retaining the ability to discriminate effectively against C and U. We found that replacement of the C5-methyl group by halogen substituents greatly decreased excision of the target base. Furthermore, 5-meC was excised more efficiently from mismatches, whereas excision of T only occurred when mispaired with G. These results suggest that ROS1 specificity arises by a combination of selective recognition at the active site and thermodynamic stability of the target base. We also found that ROS1 is a low-turnover catalyst because it binds tightly to the abasic site left after 5-meC removal. This binding leads to a highly distributive behaviour of the enzyme on DNA substrates containing multiple 5-meC residues, and may help to avoid generation of double-strand breaks during processing of bimethylated CG dinucleotides. We conclude that the biochemical properties of ROS1 are consistent with its proposed role in protecting the plant genome from excess methylation.

Histone H2a mRNA interacts with Lin28 and contains a Lin28-dependent posttranscriptional regulatory element

Xu, B., Huang, Y. — 2009-05-13

Lin28 has been shown to block the processing of let-7 microRNAs implicated in the regulation of cell growth and differentiation. Here, we show that Lin28 also specifically associates with ribonucleoprotein particles containing the replication-dependent histone H2a mRNA in mouse embryonic stem cells. We further show that the coding region of H2a mRNA harbors high affinity binding sequences for Lin28 and that these sequences stimulate the expression of reporter genes in a Lin28-dependent manner. We suggest that a key function of Lin28 in the maintenance of pluripotency is to promote the expression of the H2a gene (and perhaps also other replication-dependent histone genes) at the posttranscriptional level in order to coordinate histone production with the unique proliferative properties of embryonic stem cells.

Role of PCNA-dependent stimulation of 3'-phosphodiesterase and 3'-5' exonuclease activities of human Ape2 in repair of oxidative DNA damage

Burkovics, P., Hajdu, I., Szukacsov, V., Unk, I., Haracska, L. — 2009-05-13

Human Ape2 protein has 3' phosphodiesterase activity for processing 3'-damaged DNA termini, 3'–5' exonuclease activity that supports removal of mismatched nucleotides from the 3'-end of DNA, and a somewhat weak AP-endonuclease activity. However, very little is known about the role of Ape2 in DNA repair processes. Here, we examine the effect of interaction of Ape2 with proliferating cell nuclear antigen (PCNA) on its enzymatic activities and on targeting Ape2 to oxidative DNA lesions. We show that PCNA strongly stimulates the 3'–5' exonuclease and 3' phosphodiesterase activities of Ape2, but has no effect on its AP-endonuclease activity. Moreover, we find that upon hydrogen-peroxide treatment Ape2 redistributes to nuclear foci where it colocalizes with PCNA. In concert with these results, we provide biochemical evidence that Ape2 can reduce the mutagenic consequences of attack by reactive oxygen species not only by repairing 3'-damaged termini but also by removing 3'-end adenine opposite from 8-oxoG. Based on these findings we suggest the involvement of Ape2 in repair of oxidative DNA damage and PCNA-dependent repair synthesis.

Characterization of sINR, a strict version of the Initiator core promoter element

Yarden, G., Elfakess, R., Gazit, K., Dikstein, R. — 2009-05-13

The proximal promoter consists of binding sites for transcription regulators and a core promoter. We identified an overrepresented motif in the proximal promoter of human genes with an Initiator (INR) positional bias. The core of the motif fits the INR consensus but its sequence is more strict and flanked by additional conserved sequences. This strict INR (sINR) is enriched in TATA-less genes that belong to specific functional categories. Analysis of the sINR-containing DHX9 and ATP5F1 genes showed that the entire sINR sequence, including the strict core and the conserved flanking sequences, is important for transcription. A conventional INR sequence could not substitute for DHX9 sINR whereas, sINR could replace a conventional INR. The minimal region required to create the major TSS of the DHX9 promoter includes the sINR and an upstream Sp1 site. In a heterologous context, sINR substituted for the TATA box when positioned downstream to several Sp1 sites. Consistent with that the majority of sINR promoters contain at least one Sp1 site. Thus, sINR is a TATA-less-specific INR that functions in cooperation with Sp1. These findings support the idea that the INR is a family of related core promoter motifs.

Diverse transcription influences can be insulated by the Drosophila SF1 chromatin boundary

Majumder, P., Roy, S., Belozerov, V. E., Bosu, D., Puppali, M., Cai, H. N. — 2009-05-12

Chromatin boundaries regulate gene expression by modulating enhancer–promoter interactions and insulating transcriptional influences from organized chromatin. However, mechanistic distinctions between these two aspects of boundary function are not well understood. Here we show that SF1, a chromatin boundary located in the Drosophila Antennapedia complex (ANT-C), can insulate the transgenic miniwhite reporter from both enhancing and silencing effects of surrounding genome, a phenomenon known as chromosomal position effect or CPE. We found that the CPE-blocking activity associates with different SF1 sub-regions from a previously characterized insulator that blocks enhancers in transgenic embryos, and is independent of GAF-binding sites essential for the embryonic insulator activity. We further provide evidence that the CPE-blocking activity cannot be attributed to an enhancer-blocking activity in the developing eye. Our results suggest that SF1 contains multiple non-overlapping activities that block diverse transcriptional influences from embryonic or adult enhancers, and from positive and negative chromatin structure. Such diverse insulating capabilities are consistent with the proposed roles of SF1 to functionally separate fushi tarazu (ftz), a non-Hox gene, from the enhancers and the organized chromatin of the neighboring Hox genes.

Intra- and inter-molecular recombination of mitochondrial DNA after in vivo induction of multiple double-strand breaks

Bacman, S. R., Williams, S. L., Moraes, C. T. — 2009-05-12

To investigate mtDNA recombination induced by multiple double strand breaks (DSBs) we used a mitochondria-targeted form of the ScaI restriction endonuclease to introduce DSBs in heteroplasmic mice and cells in which we were able to utilize haplotype differences to trace the origin of recombined molecules. ScaI cleaves multiple sites in each haplotype of the heteroplasmic mice (five in NZB and three in BALB mtDNA) and prolonged expression causes severe mtDNA depletion. After a short pulse of restriction enzyme expression followed by a long period of recovery, mitochondrial genomes with large deletions were detected by PCR. Curiously, we found that some ScaI sites were more commonly involved in recombined molecules than others. In intra-molecular recombination events, deletion breakpoints were close to or upstream of ScaI cleavage sites, confirming the recombinogenic character of DSBs in mtDNA. A region adjacent to the D-loop was preferentially involved in recombination of all molecules. Sequencing through NZB and BALB haplotype markers in recombined molecules enabled us to show that in addition to intra-molecular mtDNA recombination, rare inter-molecular mtDNA recombination events can also occur. This study underscores the role of DSBs in the generation of mtDNA rearrangements and supports the existence of recombination hotspots.

The RNA binding protein Cwc2 interacts directly with the U6 snRNA to link the nineteen complex to the spliceosome during pre-mRNA splicing

McGrail, J. C., Krause, A., O'Keefe, R. T. — 2009-05-12

Intron removal during pre-messenger RNA (pre-mRNA) splicing involves arrangement of snRNAs into conformations that promote the two catalytic steps. The Prp19 complex [nineteen complex (NTC)] can specify U5 and U6 snRNA interactions with pre-mRNA during spliceosome activation. A candidate for linking the NTC to the snRNAs is the NTC protein Cwc2, which contains motifs known to bind RNA, a zinc finger and RNA recognition motif (RRM). In yeast cells mutation of either the zinc finger or RRM destabilize Cwc2 and are lethal. Yeast cells depleted of Cwc2 accumulate pre-mRNA and display reduced levels of U1, U4, U5 and U6 snRNAs. Cwc2 depletion also reduces U4/U6 snRNA complex levels, as found with depletion of other NTC proteins, but without increase in free U4. Purified Cwc2 displays general RNA binding properties and can bind both snRNAs and pre-mRNA in vitro. A Cwc2 RRM fragment alone can bind RNA but with reduced efficiency. Under splicing conditions Cwc2 can associate with U2, U5 and U6 snRNAs, but can only be crosslinked directly to the U6 snRNA. Cwc2 associates with U6 both before and after the first step of splicing. We propose that Cwc2 links the NTC to the spliceosome during pre-mRNA splicing through the U6 snRNA.

Long-range activation of FKBP51 transcription by the androgen receptor via distal intronic enhancers

Makkonen, H., Kauhanen, M., Paakinaho, V., Jaaskelainen, T., Palvimo, J. J. — 2009-05-11

Androgen receptor (AR) is a ligand-controlled transcription factor frequently deregulated in prostate carcinomas. Since there is scarce information on the action of AR on the chromatin level, we have elucidated the molecular mechanisms underlying the androgen-dependent regulation of immunophilin FKBP51 in prostate cancer cells. In comparison to the canonical AR target PSA, FKBP51 is more rapidly and strongly induced by androgen, with the regulation occurring merely at the transcriptional level. FKBP51 locus harbors 13 in silico-predicted androgen response elements (AREs), with most of them located downstream from transcription start site (TSS) and capable of binding AR in vitro. Chromatin immunoprecipitation assays in VCaP and LNCaP prostate cancer cells indicate that activation of the locus by the AR relies on four major intronic sites, with the compound ARE-containing sites ≥90 kb downstream from the TSS playing critical roles. Binding of agonist-loaded AR onto these sites in vivo was accompanied with significant recruitment of RNA polymerase II and BRM-containing chromatin remodeling complexes to the FKBP51 locus, which resulted in changes in the histone density of the locus. Our results indicate that very distal AREs act as genuine and robust enhancers, highlighting the importance of long-range regulation of transcription by the AR.

Structural insights into CUG repeats containing the 'stretched U-U wobble': implications for myotonic dystrophy

Kiliszek, A., Kierzek, R., Krzyzosiak, W. J., Rypniewski, W. — 2009-05-11

Tracks containing CUG repeats are abundant in human gene transcripts. Their biological role includes modulation of pre-mRNA splicing, mRNA transport and regulation of translation. Expanded forms of CUG runs are associated with pathogenesis of several neurodegenerative diseases, including myotonic dystrophy type 1. We have analysed two crystal structures of RNA duplexes containing the CUG repeats: G(CUG)₂C and (CUG)₆. The first of the structures, analysed at 1.23 Å resolution, is of an oligomer designed by us. The second model was obtained after ‘detwinning’ the 1.58 Å X-ray data previously deposited in the PDB. The RNA duplexes are in the A-form in which all the C–G pairs form Watson–Crick interactions while all the uridine pairs can be described as U•U cis wobble having only one hydrogen bond between the bases. The residue, which accepts the H-bond, is inclined towards the minor groove. This previously unreported base pairing can be described as ‘stretched U–U wobble’. The regular hydrogen-bonding pattern of interactions with the solvent, the electrostatic charge distribution and surface features indicate the ligand binding potential of the CUG tracks.

The formation pathway of i-motif tetramers

Leroy, J.-L. — 2009-05-11

The i-motif is a four-stranded structure formed by two intercalated parallel duplexes containing hemiprotonated C•C⁺ pairs. In order to describe the sequence of reactions by which four C-rich strands associate, we measured the formation and dissociation rates of three [TC_n]₄ tetramers (n = 3, 4 and 5), their dissociation constant and the reaction order for tetramer formation by NMR. We find that TC_n association results in the formation of several tetramers differing by the number of intercalated C•C⁺ pairs. The formation rates of the fully and partially intercalated species are comparable but their lifetimes increase strongly with the number of intercalated C•C⁺ pairs, and for this reason the single tetramer detected at equilibrium is that with optimal intercalation. The tetramer half formation times vary as the power –2 of the oligonucleotide concentration indicating that the reaction order for i-motif formation is 3. This observation is inconsistent with a model supposing association of two preformed duplex and suggests that quadruplex formation proceeds via sequential strand association into duplex and triplex intermediate species and that triplex formation is rate limiting.

Analyzing the forces binding a restriction endonuclease to DNA using a synthetic nanopore

2009-05-11

Restriction endonucleases are used prevalently in recombinant DNA technology because they bind so stably to a specific target sequence and, in the presence of cofactors, cleave double-helical DNA specifically at a target sequence at a high rate. Using synthetic nanopores along with molecular dynamics (MD), we have analyzed with atomic resolution how a prototypical restriction endonuclease, EcoRI, binds to the DNA target sequence—GAATTC—in the absence of a Mg²⁺ ion cofactor. We have previously shown that there is a voltage threshold for permeation of DNA bound to restriction enzymes through a nanopore that is associated with a nanonewton force required to rupture the complex. By introducing mutations in the DNA, we now show that this threshold depends on the recognition sequence and scales linearly with the dissociation energy, independent of the pore geometry. To predict the effect of mutation in a base pair on the free energy of dissociation, MD is used to qualitatively rank the stability of bonds in the EcoRI–DNA complex. We find that the second base in the target sequence exhibits the strongest binding to the protein, followed by the third and first bases, with even the flanking sequence affecting the binding, corroborating our experiments.

Bacillus subtilis polynucleotide phosphorylase 3'-to-5' DNase activity is involved in DNA repair

2009-05-11

In the presence of Mn²⁺, an activity in a preparation of purified Bacillus subtilis RecN degrades single-stranded (ss) DNA with a 3' -> 5' polarity. This activity is not associated with RecN itself, because RecN purified from cells lacking polynucleotide phosphorylase (PNPase) does not show the exonuclease activity. We show here that, in the presence of Mn²⁺ and low-level inorganic phosphate (P_i), PNPase degrades ssDNA. The limited end-processing of DNA is regulated by ATP and is inactive in the presence of Mg²⁺ or high-level P_i. In contrast, the RNase activity of PNPase requires Mg²⁺ and P_i, suggesting that PNPase degradation of RNA and ssDNA occur by mutually exclusive mechanisms. A null pnpA mutation (pnpA) is not epistatic with recA, but is epistatic with recN and ku, which by themselves are non-epistatic. The addA5, recO, recQ (recJ), recU and recG mutations (representative of different epistatic groups), in the context of pnpA, demonstrate gain- or loss-of-function by inactivation of repair-by-recombination, depending on acute or chronic exposure to the damaging agent and the nature of the DNA lesion. Our data suggest that PNPase is involved in various nucleic acid metabolic pathways, and its limited ssDNA exonuclease activity plays an important role in RecA-dependent and RecA-independent repair pathways.

One RNA plays three roles to provide catalytic activity to a group I intron lacking an endogenous internal guide sequence

Vaidya, N., Lehman, N. — 2009-05-11

Catalytic RNA molecules possess simultaneously a genotype and a phenotype. However, a single RNA genotype has the potential to adopt two or perhaps more distinct phenotypes as a result of differential folding and/or catalytic activity. Such multifunctionality would be particularly significant if the phenotypes were functionally inter-related in a common biochemical pathway. Here, this phenomenon is demonstrated by the ability of the Azoarcus group I ribozyme to function when its canonical internal guide sequence (GUG) has been removed from the 5' end of the molecule, and added back exogenously in trans. The presence of GUG triplets in non-covalent fragments of the ribozyme allow trans-splicing to occur in both a reverse splicing assay and a covalent self-assembly assay in which the internal guide sequence (IGS)-less ribozyme can put itself together from two of its component pieces. Analysis of these reactions indicates that a single RNA fragment can perform up to three distinct roles in a reaction: behaving as a portion of a catalyst, behaving as a substrate, and providing an exogenous IGS. This property of RNA to be multifunctional in a single reaction pathway bolsters the probability that a system of self-replicating molecules could have existed in an RNA world during the origins of life on the Earth.

Functional mapping of the interaction between TDP-43 and hnRNP A2 in vivo

2009-05-08

Nuclear factor TDP-43 has been reported to play multiple roles in transcription, pre-mRNA splicing, mRNA stability and mRNA transport. From a structural point of view, TDP-43 is a member of the hnRNP protein family whose structure includes two RRM domains flanked by the N-terminus and C-terminal regions. Like many members of this family, the C-terminal region can interact with cellular factors and thus serve to modulate its function. Previously, we have described that TDP-43 binds to several members of the hnRNP A/B family through this region. In this work, we set up a coupled minigene/siRNA cellular system that allows us to obtain in vivo data to address the functional significance of TDP-43-recruited hnRNP complex formation. Using this method, we have finely mapped the interaction between TDP-43 and the hnRNP A2 protein to the region comprised between amino acid residues 321 and 366. Our results provide novel details of protein–protein interactions in splicing regulation. In addition, we provide further insight on TDP-43 functional properties, particularly the lack of effects, as seen with our assays, of the disease-associated mutations that fall within the TDP-43 321-366 region: Q331K, M337V and G348C.

Molecular characterization of SMILE as a novel corepressor of nuclear receptors

Xie, Y.-B., Nedumaran, B., Choi, H.-S. — 2009-05-08

SMILE (small heterodimer partner interacting leucine zipper protein) has been identified as a coregulator in ER signaling. In this study, we have examined the effects of SMILE on other NRs (nuclear receptors). SMILE inhibits GR, CAR and HNF4-mediated transactivation. Knockdown of SMILE gene expression increases the transactivation of the NRs. SMILE interacts with GR, CAR and HNF4 in vitro and in vivo. SMILE and these NRs colocalize in the nucleus. SMILE binds to the ligand-binding domain or AF2 domain of the NRs. Competitions between SMILE and the coactivators GRIP1 or PGC-1 have been demonstrated in vitro and in vivo. Furthermore, an intrinsic repressive activity of SMILE is observed in Gal4-fusion system, and the intrinsic repressive domain is mapped to the C-terminus of SMILE, spanning residues 203–354. Moreover, SMILE interacts with specific HDACs (histone deacetylases) and SMILE-mediated repression is released by HDAC inhibitor trichostatin A, in a NR-specific manner. Finally, ChIP (chromatin immunoprecipitation) assays reveal that SMILE associates with the NRs on the target gene promoters. Adenoviral overexpression of SMILE represses GR-, CAR- and HNF4-mediated target gene expression. Overall, these results suggest that SMILE functions as a novel corepressor of NRs via competition with coactivators and the recruitment of HDACs.

Dynamics of RecA filaments on single-stranded DNA

van Loenhout, M. T. J., van der Heijden, T., Kanaar, R., Wyman, C., Dekker, C. — 2009-05-08

RecA, the key protein in homologous recombination, performs its actions as a helical filament on single-stranded DNA (ssDNA). ATP hydrolysis makes the RecA–ssDNA filament dynamic and is essential for successful recombination. RecA has been studied extensively by single-molecule techniques on double-stranded DNA (dsDNA). Here we directly probe the structure and kinetics of RecA interaction with its biologically most relevant substrate, long ssDNA molecules. We find that RecA ATPase activity is required for the formation of long continuous filaments on ssDNA. These filaments both nucleate and extend with a multimeric unit as indicated by the Hill coefficient of 5.4 for filament nucleation. Disassembly rates of RecA from ssDNA decrease with applied stretching force, corresponding to a mechanism where protein-induced stretching of the ssDNA aids in the disassembly. Finally, we show that RecA–ssDNA filaments can reversibly interconvert between an extended, ATP-bound, and a compressed, ADP-bound state. Taken together, our results demonstrate that ATP hydrolysis has a major influence on the structure and state of RecA filaments on ssDNA.

Stochastic sampling of the RNA structural alignment space

Harmanci, A. O., Sharma, G., Mathews, D. H. — 2009-05-08

A novel method is presented for predicting the common secondary structures and alignment of two homologous RNA sequences by sampling the ‘structural alignment’ space, i.e. the joint space of their alignments and common secondary structures. The structural alignment space is sampled according to a pseudo-Boltzmann distribution based on a pseudo-free energy change that combines base pairing probabilities from a thermodynamic model and alignment probabilities from a hidden Markov model. By virtue of the implicit comparative analysis between the two sequences, the method offers an improvement over single sequence sampling of the Boltzmann ensemble. A cluster analysis shows that the samples obtained from joint sampling of the structural alignment space cluster more closely than samples generated by the single sequence method. On average, the representative (centroid) structure and alignment of the most populated cluster in the sample of structures and alignments generated by joint sampling are more accurate than single sequence sampling and alignment based on sequence alone, respectively. The ‘best’ centroid structure that is closest to the known structure among all the centroids is, on average, more accurate than structure predictions of other methods. Additionally, cluster analysis identifies, on average, a few clusters, whose centroids can be presented as alternative candidates. The source code for the proposed method can be downloaded at http://rna.urmc.rochester.edu.

Requirement for XLF/Cernunnos in alignment-based gap filling by DNA polymerases {lambda} and {micro} for nonhomologous end joining in human whole-cell extracts

2009-05-06

XLF/Cernunnos is a core protein of the nonhomologous end-joining pathway of DNA double-strand break repair. To better define the role of Cernunnos in end joining, whole-cell extracts were prepared from Cernunnos-deficient human cells. These extracts effected little joining of DNA ends with cohesive 5' or 3' overhangs, and no joining at all of partially complementary 3' overhangs that required gap filling prior to ligation. Assays in which gap-filled but unligated intermediates were trapped using dideoxynucleotides revealed that there was no gap filling on aligned DSB ends in the Cernunnos-deficient extracts. Recombinant Cernunnos protein restored gap filling and end joining of partially complementary overhangs, and stimulated joining of cohesive ends more than twentyfold. XLF-dependent gap filling was nearly eliminated by immunodepletion of DNA polymerase , but was restored by addition of either polymerase or polymerase µ. Thus, Cernunnos is essential for gap filling by either polymerase during nonhomologous end joining, suggesting that it plays a major role in aligning the two DNA ends in the repair complex.

FHL2 interacts with and acts as a functional repressor of Id2 in human neuroblastoma cells

Han, W., Wu, Z., Zhao, Y., Meng, Y., Si, Y., Yang, J., Fu, X., Yu, L. — 2009-05-05

Inhibitor of differentiation 2 (Id2) is a natural inhibitor of the basic helix–loop–helix transcription factors. Although Id2 is well known to prevent differentiation and promote cell-cycle progression and tumorigenesis, the molecular events that regulate Id2 activity remain to be investigated. Here, we identified that Four-and-a-half LIM-only protein 2 (FHL2) is a novel functional repressor of Id2. Moreover, we demonstrated that FHL2 can directly interact with all members of the Id family (Id1–4) via an N-terminal loop–helix structure conserved in Id proteins. FHL2 antagonizes the inhibitory effect of Id proteins on basic helix–loop–helix protein E47-mediated transcription, which was abrogated by the deletion mutation of Ids that disrupted their interaction with FHL2. We also showed a competitive nature between FHL2 and E47 for binding Id2, whereby FHL2 prevents the formation of the Id2–E47 heterodimer, thus releasing E47 to DNA and restoring its transcriptional activity. FHL2 expression was remarkably up-regulated during retinoic acid-induced differentiation of neuroblastoma cells, during which the expression of Id2 was opposite to that. Ectopic FHL2 expression in neuroblastoma cells markedly reduces the transcriptional and cell-cycle promoting functions of Id2. Altogether, these results indicate that FHL2 is an important repressor of the oncogenic activity of Id2 in neuroblastoma cells.

Structural and dynamic characterization of the upper part of the HIV-1 cTAR DNA hairpin

2009-05-05

First strand transfer is essential for HIV-1 reverse transcription. During this step, the TAR RNA hairpin anneals to the cTAR DNA hairpin; this annealing reaction is promoted by the nucleocapsid protein and involves an initial loop–loop interaction between the apical loops of TAR and cTAR. Using NMR and probing methods, we investigated the structural and dynamic properties of the top half of the cTAR DNA (mini-cTAR). We show that the upper stem located between the apical and the internal loops is stable, but that the lower stem of mini-cTAR is unstable. The residues of the internal loop undergo slow motions at the NMR time-scale that are consistent with conformational exchange phenomena. In contrast, residues of the apical loop undergo fast motions. The lower stem is destabilized by the slow interconversion processes in the internal loop, and thus the internal loop is responsible for asymmetric destabilization of mini-cTAR. These findings are consistent with the functions of cTAR in first strand transfer: its apical loop is suitably exposed to interact with the apical loop of TAR RNA and its lower stem is significantly destabilized to facilitate the subsequent action of the nucleocapsid protein which promotes the annealing reaction.

Contribution of ribosomal residues to P-site tRNA binding

Shoji, S., Abdi, N. M., Bundschuh, R., Fredrick, K. — 2009-05-05

Structural studies have revealed multiple contacts between the ribosomal P site and tRNA, but how these contacts contribute to P-tRNA binding remains unclear. In this study, the effects of ribosomal mutations on the dissociation rate (k_off) of various tRNAs from the P site were measured. Mutation of the 30S P site destabilized tRNAs to various degrees, depending on the mutation and the species of tRNA. These data support the idea that ribosome-tRNA interactions are idiosyncratically tuned to ensure stable binding of all tRNA species. Unlike deacylated elongator tRNAs, N-acetyl-aminoacyl-tRNAs and tRNA^fMet dissociated from the P site at a similar low rate, even in the presence of various P-site mutations. These data provide evidence for a stability threshold for P-tRNA binding and suggest that ribosome-tRNA^fMet interactions are uniquely tuned for tight binding. The effects of 16S rRNA mutation G1338U were suppressed by 50S E-site mutation C2394A, suggesting that G1338 is particularly important for stabilizing tRNA in the P/E site. Finally, mutation C2394A or the presence of an N-acetyl-aminoacyl group slowed the association rate (k_on) of tRNA dramatically, suggesting that deacylated tRNA binds the P site of the ribosome via the E site.

Comprehensive prediction of novel microRNA targets in Arabidopsis thaliana

Alves-Junior, L., Niemeier, S., Hauenschild, A., Rehmsmeier, M., Merkle, T. — 2009-05-05

MicroRNAs (miRNAs) are 20–24 nt long endogenous non-coding RNAs that act as post-transcriptional regulators in metazoa and plants. Plant miRNA targets typically contain a single sequence motif with near-perfect complementarity to the miRNA. Here, we extended and applied the program RNAhybrid to identify novel miRNA targets in the complete annotated Arabidopsis thaliana transcriptome. RNAhybrid predicts the energetically most favorable miRNA:mRNA hybrids that are consistent with user-defined structural constraints. These were: (i) perfect base pairing of the duplex from nucleotide 8 to 12 counting from the 5'-end of the miRNA; (ii) loops with a maximum length of one nucleotide in either strand; (iii) bulges with no more than one nucleotide in size; and (iv) unpaired end overhangs not longer than two nucleotides. G:U base pairs are not treated as mismatches, but contribute less favorable to the overall free energy. The resulting hybrids were filtered according to their minimum free energy, resulting in an overall prediction of more than 600 novel miRNA targets. The specificity and signal-to-noise ratio of the prediction was assessed with either randomized miRNAs or randomized target sequences as negative controls. Our results are in line with recent observations that the majority of miRNA targets are not transcription factors.

Folding of the cocaine aptamer studied by EPR and fluorescence spectroscopies using the bifunctional spectroscopic probe C

Cekan, P., Jonsson, E. O., Sigurdsson, S. Th. — 2009-04-30

The cocaine aptamer is a DNA molecule that binds cocaine at the junction of three helices. The bifunctional spectroscopic probe Ç was incorporated independently into three different positions of the aptamer and changes in structure and dynamics upon addition of the cocaine ligand were studied. Nucleoside Ç contains a rigid nitroxide spin label and can be studied directly by electron paramagnetic resonance (EPR) spectroscopy and fluorescence spectroscopy after reduction of the nitroxide to yield the fluoroside Ç^f. Both the EPR and the fluorescence data for aptamer 2 indicate that helix III is formed before cocaine binding. Upon addition of cocaine, increased fluorescence of a fully base-paired Ç^f, placed at the three-way junction in helix III, was observed and is consistent with a helical tilt from a coaxial stack of helices II and III. EPR and fluorescence data clearly show that helix I is formed upon addition of cocaine, concomitant with the formation of the Y-shaped three-way helical junction. The EPR data indicate that nucleotides in helix I are more mobile than nucleotides in regular duplex regions and may reflect increased dynamics due to the short length of helix I.

A glycolytic burst drives glucose induction of global histone acetylation by picNuA4 and SAGA

2009-04-30

Little is known about what enzyme complexes or mechanisms control global lysine acetylation in the amino-terminal tails of the histones. Here, we show that glucose induces overall acetylation of H3 K9, 18, 27 and H4 K5, 8, 12 in quiescent yeast cells mainly by stimulating two KATs, Gcn5 and Esa1. Genetic and pharmacological perturbation of carbon metabolism, combined with ¹H-NMR metabolic profiling, revealed that glucose induction of KAT activity directly depends on increased glucose catabolism. Glucose-inducible Esa1 and Gcn5 activities predominantly reside in the picNuA4 and SAGA complexes, respectively, and act on chromatin by an untargeted mechanism. We conclude that direct metabolic regulation of globally acting KATs can be a potent driving force for reconfiguration of overall histone acetylation in response to a physiological cue.

The ATR-Chk1 pathway plays a role in the generation of centrosome aberrations induced by Rad51C dysfunction

2009-04-29

Rad51C is a central component of two complexes formed by five Rad51 paralogs in vertebrates. These complexes are involved in repairing DNA double-strand breaks through homologous recombination. Despite accumulating evidence suggesting that the paralogs may prevent aneuploidy by controlling centrosome integrity, Rad51C's role in maintaining chromosome stability remains unclear. Here we demonstrate that Rad51C deficiency leads to both centrosome aberrations in an ATR-Chk1-dependent manner and increased aneuploidy in human cells. While it was reported that Rad51C deficiency did not cause centrosome aberrations in interphase in hamster cells, such aberrations were observed in interphase in HCT116 cells with Rad51C dysfunction. Caffeine treatment and down-regulation of ATR, but not that of ATM, reduced the frequency of centrosome aberrations in the mutant cells. Silencing of Rad51C by RNA interference in HT1080 cells resulted in similar aberrations. Treatment with a Chk1 inhibitor and silencing of Chk1 also reduced the frequency in HCT116 mutants. Accumulation of Chk1 at the centrosome and nuclear foci of H2AX were increased in the mutants. Moreover, the mutant cells had a higher frequency of aneuploidy. These findings indicate that the ATR-Chk1 pathway plays a role in increased centrosome aberrations induced by Rad51C dysfunction.

Identification of Apurinic/apyrimidinic endonuclease 1 (APE1) as the endoribonuclease that cleaves c-myc mRNA

2009-04-28

Endonucleolytic cleavage of the coding region determinant (CRD) of c-myc mRNA appears to play a critical role in regulating c-myc mRNA turnover. Using ³²P-labeled c-myc CRD RNA as substrate, we have purified and identified two endoribonucleases from rat liver polysomes that are capable of cleaving the transcript in vitro. A 17-kDa enzyme was identified as RNase1. Apurinic/apyrimidinic (AP) DNA endonuclease 1 (APE1) was identified as the 35-kDa endoribonuclease that preferentially cleaves in between UA and CA dinucleotides of c-myc CRD RNA. APE1 was further confirmed to be the 35-kDa endoribonuclease because: (i) the endoribonuclease activity of the purified 35-kDa native enzyme was specifically immuno-depleted with APE1 monoclonal antibody, and (ii) recombinant human APE1 generated identical RNA cleavage patterns as the native liver enzyme. Studies using E96A and H309N mutants of APE1 suggest that the endoribonuclease activity for c-myc CRD RNA shares the same active center with the AP-DNA endonuclease activity. Transient knockdown of APE1 in HeLa cells led to increased steady-state level of c-myc mRNA and its half-life. We conclude that the ability to cleave RNA dinucleotides is a previously unidentified function of APE1 and it can regulate c-myc mRNA level possibly via its endoribonuclease activity.

S-Adenosyl homocysteine and DNA ends stimulate promiscuous nuclease activities in the Type III restriction endonuclease EcoPI

Peakman, L. J., Szczelkun, M. D. — 2009-04-28

In the absence of the methyl donor S-adenosyl methionine and under certain permissive reaction conditions, EcoPI shows non-specific endonuclease activity. We show here that the cofactor analogue S-adenosyl homocysteine promotes this promiscuous DNA cleavage. Additionally, an extensive exonuclease-like processing of the DNA is also observed that can even result in digestion of non-specific DNA in trans. We suggest a model for how DNA communication events initiating from non-specific sites, and in particular free DNA ends, could produce the observed cleavage patterns.

Highly efficient incorporation of the fluorescent nucleotide analogs tC and tCO by Klenow fragment

2009-04-28

Studies of the mechanisms by which DNA polymerases select the correct nucleotide frequently employ fluorescently labeled DNA to monitor conformational rearrangements of the polymerase–DNA complex in response to incoming nucleotides. For this purpose, fluorescent base analogs play an increasingly important role because they interfere less with the DNA–protein interaction than do tethered fluorophores. Here we report the incorporation of the 5'-triphosphates of two exceptionally bright cytosine analogs, 1,3-diaza-2-oxo-phenothiazine (tC) and its oxo-homolog, 1,3-diaza-2-oxo-phenoxazine (tC^O), into DNA by the Klenow fragment. Both nucleotide analogs are polymerized with slightly higher efficiency opposite guanine than cytosine triphosphate and are shown to bind with nanomolar affinity to the DNA polymerase active site, according to fluorescence anisotropy measurements. Using this method, we perform competitive binding experiments and show that they can be used to determine the dissociation constant of any given natural or unnatural nucleotide. The results demonstrate that the active site of the Klenow fragment is flexible enough to tolerate base pairs that are size-expanded in the major groove. In addition, the possibility to enzymatically polymerize a fluorescent nucleotide with high efficiency complements the tool box of biophysical probes available to study DNA replication.

Remodeling of chromatin structure within the promoter is important for bmp-2-induced fgfr3 expression

Sun, F., Chen, Q., Yang, S., Pan, Q., Ma, J., Wan, Y., Chang, C.-H., Hong, A. — 2009-04-28

Fibroblast growth factor receptor 3 (FGFR3) plays an important role in cartilage development. Although upregulation of FGFR3 expression in response to bone morphogenetic protein-2 (BMP-2) has been reported, the molecular mechanisms remain unknown. In this study, we used in vivo approaches to characterize BMP-2-induced alterations in the chromatin organization of the FGFR3 core promoter. Chromatin immunoprecipitation analysis demonstrated that the binding of Brg1, a component of the SWI/SNF remodeling complex, may selectively remodel a chromatin region (encompassing nucleotide –90 to +35), uncovering the transcription start site and three Sp1-binding sites, as revealed by nuclease digestion hypersensitivity assays. We then showed an increase in the association of Sp1 with the proximal promoter, followed by the recruitment of p300, resulting in a change of the histone ‘code’, such as in phosphorylation and methylation. Collectively, our study results suggest a model for BMP-2-induced FGFR3 expression in which the core promoter architecture is specifically regulated.

Cellular responses to DNA double-strand breaks after low-dose {gamma}-irradiation

Asaithamby, A., Chen, D. J. — 2009-04-28

DNA double-strand breaks (DSBs) are a serious threat to genome stability and cell viability. Although biological effects of low levels of radiation are not clear, the risks of low-dose radiation are of societal importance. Here, we directly monitored induction and repair of single DSBs and quantitatively analyzed the dynamics of interaction of DNA repair proteins at individual DSB sites in living cells using 53BP1 fused to yellow fluorescent protein (YFP-53BP1) as a surrogate marker. The number of DSBs formed was linear with dose from 5 mGy to 1 Gy. The DSBs induced by very low radiation doses (5 mGy) were repaired with efficiency similar to repair of DSBs induced at higher doses. The YFP-53BP1 foci are dynamic structures: 53BP1 rapidly and reversibly interacted at these DSB sites. The time frame of recruitment and affinity of 53BP1 for DSB sites were indistinguishable between low and high doses, providing mechanistic evidence for the similar DSB repair after low- and high-dose radiation. These findings have important implications for estimating the risk associated with low-dose radiation exposure on human health.

IHF-binding sites inhibit DNA loop formation and transcription initiation

Huo, Y.-X., Zhang, Y.-T., Xiao, Y., Zhang, X., Buck, M., Kolb, A., Wang, Y.-P. — 2009-04-24

Transcriptional activation of enhancer and ⁵⁴-dependent promoters requires efficient interactions between enhancer-binding proteins (EBP) and promoter bound ⁵⁴-RNA polymerase (E⁵⁴) achieved by DNA looping, which is usually facilitated by the integration host factor (IHF). Since the lengths of the intervening region supporting DNA-loop formation are similar among IHF-dependent and IHF-independent promoters, the precise reason(s) why IHF is selectively important for the frequency of transcription initiation remain unclear. Here, using kinetic cyclization and in vitro transcription assays we show that, in the absence of IHF protein, the DNA fragments containing an IHF-binding site have much less looping-formation ability than those that lack an IHF-binding site. Furthermore, when an IHF consensus-binding site was introduced into the intervening region between promoter and enhancer of the target DNA fragments, loop formation and DNA-loop-dependent transcriptional activation are significantly reduced in a position-independent manner. DNA-looping-independent transcriptional activation was unaffected. The binding of IHF to its consensus site in the target promoters clearly restored efficient DNA looping formation and looping-dependent transcriptional activation. Our data provide evidence that one function for the IHF protein is to release a communication block set by intrinsic properties of the IHF DNA-binding site.

Both helix topology and counterion distribution contribute to the more effective charge screening in dsRNA compared with dsDNA

Pabit, S. A., Qiu, X., Lamb, J. S., Li, L., Meisburger, S. P., Pollack, L. — 2009-04-24

The recent discovery of the RNA interference mechanism emphasizes the biological importance of short, isolated, double-stranded (ds) RNA helices and calls for a complete understanding of the biophysical properties of dsRNA. However, most previous studies of the electrostatics of nucleic acid duplexes have focused on DNA. Here, we present a comparative investigation of electrostatic effects in RNA and DNA. Using resonant (anomalous) and non-resonant small-angle X-ray scattering, we characterized the charge screening efficiency and counterion distribution around short (25 bp) dsDNA and RNA molecules of comparable sequence. Consistent with theoretical predictions, we find counterion mediated screening to be more efficient for dsRNA than dsDNA. Furthermore, the topology of the RNA A-form helix alters the spatial distribution of counterions relative to B-form DNA. The experimental results reported here agree well with ion-size-corrected non-linear Poisson–Boltzmann calculations. We propose that differences in electrostatic properties aid in selective recognition of different types of short nucleic acid helices by target binding partners.

High-resolution genome-wide cytosine methylation profiling with simultaneous copy number analysis and optimization for limited cell numbers

2009-04-22

Many genome-wide assays involve the generation of a subset (or representation) of the genome following restriction enzyme digestion. The use of enzymes sensitive to cytosine methylation allows high-throughput analysis of this epigenetic regulatory process. We show that the use of a dual-adapter approach allows us to generate genomic representations that includes fragments of <200 bp in size, previously not possible when using the standard approach of using a single adapter. By expanding the representation to smaller fragments using HpaII or MspI, we increase the representation by these isoschizomers to more than 1.32 million loci in the human genome, representing 98.5% of CpG islands and 91.1% of refSeq promoters. This advance allows the development of a new, high-resolution version of our HpaII-tiny fragment Enrichment by Ligation-mediated PCR (HELP) assay to study cytosine methylation. We also show that the MspI representation generates information about copy-number variation, that the assay can be used on as little as 10 ng of DNA and that massively parallel sequencing can be used as an alternative to microarrays to read the output of the assay, making this a powerful discovery platform for studies of genomic and epigenomic abnormalities.

A genome-wide screen for essential yeast genes that affect telomere length maintenance

Ungar, L., Yosef, N., Sela, Y., Sharan, R., Ruppin, E., Kupiec, M. — 2009-04-22

Telomeres are structures composed of repetitive DNA and proteins that protect the chromosomal ends in eukaryotic cells from fusion or degradation, thus contributing to genomic stability. Although telomere length varies between species, in all organisms studied telomere length appears to be controlled by a dynamic equilibrium between elongating mechanisms (mainly addition of repeats by the enzyme telomerase) and nucleases that shorten the telomeric sequences. Two previous studies have analyzed a collection of yeast deletion strains (deleted for nonessential genes) and found over 270 genes that affect telomere length (Telomere Length Maintenance or TLM genes). Here we complete the list of TLM by analyzing a collection of strains carrying hypomorphic alleles of most essential genes (DAmP collection). We identify 87 essential genes that affect telomere length in yeast. These genes interact with the nonessential TLM genes in a significant manner, and provide new insights on the mechanisms involved in telomere length maintenance. The newly identified genes span a variety of cellular processes, including protein degradation, pre-mRNA splicing and DNA replication.

Oligonucleotides suppress PKB/Akt and act as superinductors of apoptosis in human keratinocytes

2009-04-22

DNA oligonucleotides (ODN) applied to an organism are known to modulate the innate and adaptive immune system. Previous studies showed that a CpG-containing ODN (CpG-1-PTO) and interestingly, also a non-CpG-containing ODN (nCpG-5-PTO) suppress inflammatory markers in skin. In the present study it was investigated whether these molecules also influence cell apoptosis. Here we show that CpG-1-PTO, nCpG-5-PTO, and also natural DNA suppress the phosphorylation of PKB/Akt in a cell-type-specific manner. Interestingly, only epithelial cells of the skin (normal human keratinocytes, HaCaT and A-431) show a suppression of PKB/Akt. This suppressive effect depends from ODN lengths, sequence and backbone. Moreover, it was found that TGF-induced levels of PKB/Akt and EGFR were suppressed by the ODN tested. We hypothesize that this suppression might facilitate programmed cell death. By testing this hypothesis we found an increase of apoptosis markers (caspase 3/7, 8, 9, cytosolic cytochrome c, histone associated DNA fragments, apoptotic bodies) when cells were treated with ODN in combination with low doses of staurosporin, a well-known pro-apoptotic stimulus. In summary the present data demonstrate DNA as a modulator of apoptosis which specifically targets skin epithelial cells.

Structural basis for m7G-cap hypermethylation of small nuclear, small nucleolar and telomerase RNA by the dimethyltransferase TGS1

Monecke, T., Dickmanns, A., Ficner, R. — 2009-04-22

The 5'-cap of spliceosomal small nuclear RNAs, some small nucleolar RNAs and of telomerase RNA was found to be hypermethylated in vivo. The Trimethylguanosine Synthase 1 (TGS1) mediates this conversion of the 7-methylguanosine-cap to the 2,2,7-trimethylguanosine (m₃G)-cap during maturation of the RNPs. For mammalian UsnRNAs the generated m^2,2,7G-cap is one part of a bipartite import signal mediating the transport of the UsnRNP-core complex into the nucleus. In order to understand the structural organization of human TGS1 as well as substrate binding and recognition we solved the crystal structure of the active TGS1 methyltransferase domain containing both, the minimal substrate m⁷GTP and the reaction product S-adenosyl-l-homocysteine (AdoHcy). The methyltransferase of human TGS1 harbors the canonical class 1 methyltransferase fold as well as an unique N-terminal, -helical domain of 40 amino acids, which is essential for m⁷G-cap binding and catalysis. The crystal structure of the substrate bound methyltransferase domain as well as mutagenesis studies provide insight into the catalytic mechanism of TGS1.

Evidence of genome-wide G4 DNA-mediated gene expression in human cancer cells

Verma, A., Yadav, V. K., Basundra, R., Kumar, A., Chowdhury, S. — 2009-02-11

Guanine-rich DNA of a particular sequence adopts four-stranded structural forms known as G-quadruplex or G4 DNA. Though in vitro formation of G4 DNA is known for several years, in vivo presence of G4 DNA was only recently noted in eukaryote telomeres. Recent bioinformatics analyses showing prevalence of G4 DNA within promoters of human and related species seems to implicate G4 DNA in a genome-wide cis-regulatory role. Herein we demonstrate that G4 DNA may present regulatory sites on a genome-wide scale by showing widespread effect on gene expression in response to the established intracellular G4 DNA-binding ligands. This is particularly relevant to genes that harbor conserved potential G4 DNA (PG4 DNA) forming sequence across human, mouse and rat promoters of orthologous genes. Genes with conserved PG4 DNA in promoters show co-regulated expression in 79 human and 61 mouse normal tissues (z-score > 3.5; P < 0.0001). Conservation of G4 DNA across related species also emphasizes the biological importance of G4 DNA and its role in transcriptional regulation of genes; shedding light on a relatively novel mechanism of regulation of gene expression in eukaryotes.

This paper has been withdrawn

2007-06-12