Nucleic Acids Research

A subgroup of SGS3-like proteins act redundantly in RNA-directed DNA methylation.

Publication Date: 2012 Feb 1 PMID: 22302148
Authors: Xie, M. - Ren, G. - Costa-Nunes, P. - Pontes, O. - Yu, B.
Journal: Nucleic Acids Res

Plant specific SGS3-like proteins are composed of various combinations of an RNA-binding XS domain, a zinc-finger zf-XS domain, a coil-coil domain and a domain of unknown function called XH. In addition to being involved in de novo 2 (IDN2) and SGS3, the Arabidopsis genome encodes 12 uncharacterized SGS3-like proteins. Here, we show that a group of SGS3-like proteins act redundantly in RNA-directed DNA methylation (RdDM) pathway in Arabidopsis. Transcriptome co-expression analyses reveal significantly correlated expression of two SGS3-like proteins, factor of DNA methylation 1 (FDM1) and FDM2 with known genes required for RdDM. The fdm1 and fdm2 double mutations but not the fdm1 or fdm2 single mutations significantly impair DNA methylation at RdDM loci, release transcriptional gene silencing and dramatically reduce the abundance of siRNAs originated from high copy number repeats or transposons. Like IDN2 and SGS3, FDM1 binds dsRNAs with 5' overhangs. Double mutant analyses also reveal that IDN2 and three uncharacterized SGS3-like proteins FDM3, FDM4 and FDM5 have overlapping function with FDM1 in RdDM. Five FDM proteins and IDN2 define a group of SGS3-like proteins that possess all four-signature motifs in Arabidopsis. Thus, our results demonstrate that this group of SGS3-like proteins is an important component of RdDM. This study further enhances our understanding of the SGS3 gene family and the RdDM pathway.

post to: CiteULike

cn.MOPS: mixture of Poissons for discovering copy number variations in next-generation sequencing data with a low false discovery rate.

Publication Date: 2012 Feb 1 PMID: 22302147
Authors: Klambauer, G. - Schwarzbauer, K. - Mayr, A. - Clevert, D. A. - Mitterecker, A. - Bodenhofer, U. - Hochreiter, S.
Journal: Nucleic Acids Res

Quantitative analyses of next-generation sequencing (NGS) data, such as the detection of copy number variations (CNVs), remain challenging. Current methods detect CNVs as changes in the depth of coverage along chromosomes. Technological or genomic variations in the depth of coverage thus lead to a high false discovery rate (FDR), even upon correction for GC content. In the context of association studies between CNVs and disease, a high FDR means many false CNVs, thereby decreasing the discovery power of the study after correction for multiple testing. We propose 'Copy Number estimation by a Mixture Of PoissonS' (cn.MOPS), a data processing pipeline for CNV detection in NGS data. In contrast to previous approaches, cn.MOPS incorporates modeling of depths of coverage across samples at each genomic position. Therefore, cn.MOPS is not affected by read count variations along chromosomes. Using a Bayesian approach, cn.MOPS decomposes variations in the depth of coverage across samples into integer copy numbers and noise by means of its mixture components and Poisson distributions, respectively. The noise estimate allows for reducing the FDR by filtering out detections having high noise that are likely to be false detections. We compared cn.MOPS with the five most popular methods for CNV detection in NGS data using four benchmark datasets: (i) simulated data, (ii) NGS data from a male HapMap individual with implanted CNVs from the X chromosome, (iii) data from HapMap individuals with known CNVs, (iv) high coverage data from the 1000 Genomes Project. cn.MOPS outperformed its five competitors in terms of precision (1-FDR) and recall for both gains and losses in all benchmark data sets. The software cn.MOPS is publicly available as an R package at http://www.bioinf.jku.at/software/cnmops/ and at Bioconductor.

post to: CiteULike

Fidelity of tRNA 5'-maturation: a possible basis for the functional dependence of archaeal and eukaryal RNase P on multiple protein cofactors.

Publication Date: 2012 Jan 31 PMID: 22298511
Authors: Chen, W. Y. - Singh, D. - Lai, L. B. - Stiffler, M. A. - Lai, H. D. - Foster, M. P. - Gopalan, V.
Journal: Nucleic Acids Res

RNase P, which catalyzes tRNA 5'-maturation, typically comprises a catalytic RNase P RNA (RPR) and a varying number of RNase P proteins (RPPs): 1 in bacteria, at least 4 in archaea and 9 in eukarya. The four archaeal RPPs have eukaryotic homologs and function as heterodimers (POP5*RPP30 and RPP21*RPP29). By studying the archaeal Methanocaldococcus jannaschii RPR's cis cleavage of precursor tRNA(Gln) (pre-tRNA(Gln)), which lacks certain consensus structures/sequences needed for substrate recognition, we demonstrate that RPP21*RPP29 and POP5*RPP30 can rescue the RPR's mis-cleavage tendency independently by 4-fold and together by 25-fold, suggesting that they operate by distinct mechanisms. This synergistic and preferential shift toward correct cleavage results from the ability of archaeal RPPs to selectively increase the RPR's apparent rate of correct cleavage by 11 140-fold, compared to only 480-fold for mis-cleavage. Moreover, POP5*RPP30, like the bacterial RPP, helps normalize the RPR's rates of cleavage of non-consensus and consensus pre-tRNAs. We also show that archaeal and eukaryal RNase P, compared to their bacterial relatives, exhibit higher fidelity of 5'-maturation of pre-tRNA(Gln) and some of its mutant derivatives. Our results suggest that protein-rich RNase P variants might have evolved to support flexibility in substrate recognition while catalyzing efficient, high-fidelity 5'-processing.

post to: CiteULike

An on-bead tailing/ligation approach for sequencing resin-bound RNA libraries.

Publication Date: 2012 Jan 31 PMID: 22298510
Authors: Wiesmayr, A. - Fournier, P. - Jaschke, A.
Journal: Nucleic Acids Res

Nucleic acids possess the unique property of being enzymatically amplifiable, and have therefore been a popular choice for the combinatorial selection of functional sequences, such as aptamers or ribozymes. However, amplification typically requires known sequence segments that serve as primer binding sites, which can be limiting for certain applications, like the screening of on-bead libraries. Here, we report a method to amplify and sequence on-bead RNA libraries that requires not more than five known nucleotides. A key element is the attachment of the starting nucleoside to the synthesis resin via the nucleobase, which leaves the 3'-OH group accessible to subsequent enzymatic manipulations. After split-and-mix synthesis of the oligonucleotide library and deprotection, a poly(A)-tail can be efficiently added to this free 3'-hydroxyl terminus by Escherichia coli poly(A) polymerase that serves as an anchored primer binding site for reverse transcription. The cDNA is joined to a DNA adapter by T4 DNA ligase. PCR amplification yielded single-band products that could be cloned and sequenced starting from individual polystyrene beads. The method described here makes the selection of functional RNAs from on-bead RNA libraries more attractive due to increased flexibility in library design, higher yields of full-length sequence on bead and robust sequence determination.

post to: CiteULike

Disparity in the DNA translocase domains of SWI/SNF and ISW2.

Publication Date: 2012 Jan 31 PMID: 22298509
Authors: Dechassa, M. L. - Hota, S. K. - Sen, P. - Chatterjee, N. - Prasad, P. - Bartholomew, B.
Journal: Nucleic Acids Res

An ATP-dependent DNA translocase domain consisting of seven conserved motifs is a general feature of all ATP-dependent chromatin remodelers. While motifs on the ATPase domains of the yeast SWI/SNF and ISWI families of remodelers are highly conserved, the ATPase domains of these complexes appear not to be functionally interchangeable. We found one reason that may account for this is the ATPase domains interact differently with nucleosomes even though both associate with nucleosomal DNA 17-18 bp from the dyad axis. The cleft formed between the two lobes of the ISW2 ATPase domain is bound to nucleosomal DNA and Isw2 associates with the side of nucleosomal DNA away from the histone octamer. The ATPase domain of SWI/SNF binds to the same region of nucleosomal DNA, but is bound outside of the cleft region. The catalytic subunit of SWI/SNF also appears to intercalate between the DNA gyre and histone octamer. The altered interactions of SWI/SNF with DNA are specific to nucleosomes and do not occur with free DNA. These differences are likely mediated through interactions with the histone surface. The placement of SWI/SNF between the octamer and DNA could make it easier to disrupt histone-DNA interactions.

post to: CiteULike

Polyadenylation helps regulate functional tRNA levels in Escherichia coli.

Publication Date: 2012 Jan 28 PMID: 22287637
Authors: Mohanty, B. K. - Maples, V. F. - Kushner, S. R.
Journal: Nucleic Acids Res

Here we demonstrate a new regulatory mechanism for tRNA processing in Escherichia coli whereby RNase T and RNase PH, the two primary 3' --> 5' exonucleases involved in the final step of 3'-end maturation, compete with poly(A) polymerase I (PAP I) for tRNA precursors in wild-type cells. In the absence of both RNase T and RNase PH, there is a >30-fold increase of PAP I-dependent poly(A) tails that are 2-fold decrease in growth rate. Only 7 out of 86 tRNAs are not regulated by this mechanism and are also not substrates for RNase T, RNase PH or PAP I. Surprisingly, neither PNPase nor RNase II has any effect on tRNA poly(A) tail length. Our data suggest that the polyadenylation of tRNAs by PAP I likely proceeds in a distributive fashion unlike what is observed with mRNAs.

post to: CiteULike

Chromatin domain boundary element search tool for Drosophila.

Publication Date: 2012 Jan 28 PMID: 22287636
Authors: Srinivasan, A. - Mishra, R. K.
Journal: Nucleic Acids Res

Chromatin domain boundary elements prevent inappropriate interaction between distant or closely spaced regulatory elements and restrict enhancers and silencers to correct target promoters. In spite of having such a general role and expected frequent occurrence genome wide, there is no DNA sequence analysis based tool to identify boundary elements. Here, we report chromatin domain Boundary Element Search Tool (cdBEST), to identify boundary elements. cdBEST uses known recognition sequences of boundary interacting proteins and looks for 'motif clusters'. Using cdBEST, we identified boundary sequences across 12 Drosophila species. Of the 4576 boundary sequences identified in Drosophila melanogaster genome, >170 sequences are repetitive in nature and have sequence homology to transposable elements. Analysis of such sequences across 12 Drosophila genomes showed that the occurrence of repetitive sequences in the context of boundaries is a common feature of drosophilids. We use a variety of genome organization criteria and also experimental test on a subset of the cdBEST boundaries in an enhancer-blocking assay and show that 80% of them indeed function as boundaries in vivo. These observations highlight the role of cdBEST in better understanding of chromatin domain boundaries in Drosophila and setting the stage for comparative analysis of boundaries across closely related species.

post to: CiteULike

Transcriptional regulation of gene expression during osmotic stress responses by the mammalian target of rapamycin.

Publication Date: 2012 Jan 28 PMID: 22287635
Authors: Ortells, M. C. - Morancho, B. - Drews-Elger, K. - Viollet, B. - Laderoute, K. R. - Lopez-Rodriguez, C. - Aramburu, J.
Journal: Nucleic Acids Res

Although stress can suppress growth and proliferation, cells can induce adaptive responses that allow them to maintain these functions under stress. While numerous studies have focused on the inhibitory effects of stress on cell growth, less is known on how growth-promoting pathways influence stress responses. We have approached this question by analyzing the effect of mammalian target of rapamycin (mTOR), a central growth controller, on the osmotic stress response. Our results showed that mammalian cells exposed to moderate hypertonicity maintained active mTOR, which was required to sustain their cell size and proliferative capacity. Moreover, mTOR regulated the induction of diverse osmostress response genes, including targets of the tonicity-responsive transcription factor NFAT5 as well as NFAT5-independent genes. Genes sensitive to mTOR-included regulators of stress responses, growth and proliferation. Among them, we identified REDD1 and REDD2, which had been previously characterized as mTOR inhibitors in other stress contexts. We observed that mTOR facilitated transcription-permissive conditions for several osmoresponsive genes by enhancing histone H4 acetylation and the recruitment of RNA polymerase II. Altogether, these results reveal a previously unappreciated role of mTOR in regulating transcriptional mechanisms that control gene expression during cellular stress responses.

post to: CiteULike

Performance comparison and evaluation of software tools for microRNA deep-sequencing data analysis.

Publication Date: 2012 Jan 28 PMID: 22287634
Authors: Li, Y. - Zhang, Z. - Liu, F. - Vongsangnak, W. - Jing, Q. - Shen, B.
Journal: Nucleic Acids Res

With the development of next-generation sequencing (NGS) techniques, many software tools have emerged for the discovery of novel microRNAs (miRNAs) and for analyzing the miRNAs expression profiles. An overall evaluation of these diverse software tools is lacking. In this study, we evaluated eight software tools based on their common feature and key algorithms. Three deep-sequencing data sets were collected from different species and used to assess the computational time, sensitivity and accuracy of detecting known miRNAs as well as their capacity for predicting novel miRNAs. Our results provide useful information for researchers to facilitate their selection of the optimal software tools for miRNA analysis depending on their specific requirements, i.e. novel miRNAs discovery or miRNA expression profile analysis of sequencing data sets.

post to: CiteULike

DNA robustly stimulates FANCD2 monoubiquitylation in the complex with FANCI.

Publication Date: 2012 Jan 28 PMID: 22287633
Authors: Sato, K. - Toda, K. - Ishiai, M. - Takata, M. - Kurumizaka, H.
Journal: Nucleic Acids Res

FANCI and FANCD2 form a complex, and play essential roles in the repair of interstrand DNA crosslinks (ICLs) by the Fanconi anemia (FA) pathway. FANCD2 is monoubiquitylated by the FA core complex, composed of 10 FA proteins including FANCL as the catalytic E3 subunit. FANCD2 monoubiquitylation can be reconstituted with purified minimal components, such as FANCI, E1, UBE2T (E2) and FANCL (E3) in vitro; however, its efficiency is quite low as compared to the in vivo monoubiquitylation of FANCD2. In this study, we found that various forms of DNA, such as single-stranded, double-stranded and branched DNA, robustly stimulated the FANCD2 monoubiquitylation in vitro up to a level comparable to its in vivo monoubiquitylation. This stimulation of the FANCD2 monoubiquitylation strictly required FANCI, suggesting that FANCD2 monoubiquitylation may occur in the FANCI-FANCD2 complex. A FANCI mutant that was defective in DNA binding was also significantly defective in FANCD2 monoubiquitylation in vitro. In the presence of 5' flapped DNA, a DNA substrate mimicking the arrested replication fork, about 70% of the input FANCD2 was monoubiquitylated, while less than 1% FANCD2 monoubiquitylation was observed in the absence of the DNA. Therefore, DNA may be the unidentified factor required for proper FANCD2 monoubiquitylation.

post to: CiteULike

Association of the interferon-beta gene with pericentromeric heterochromatin is dynamically regulated during virus infection through a YY1-dependent mechanism.

Publication Date: 2012 Jan 28 PMID: 22287632
Authors: Josse, T. - Mokrani-Benhelli, H. - Benferhat, R. - Shestakova, E. - Mansuroglu, Z. - Kakanakou, H. - Billecocq, A. - Bouloy, M. - Bonnefoy, E.
Journal: Nucleic Acids Res

Nuclear architecture as well as gene nuclear positioning can modulate gene expression. In this work, we have analyzed the nuclear position of the interferon-beta (IFN-beta) locus, responsible for the establishment of the innate antiviral response, with respect to pericentromeric heterochromatin (PCH) in correlation with virus-induced IFN-beta gene expression. Experiments were carried out in two different cell types either non-infected (NI) or during the time course of three different viral infections. In NI cells, we showed a monoallelic IFN-beta promoter association with PCH that strongly decreased after viral infection. Dissociation of the IFN-beta locus away from these repressive regions preceded strong promoter transcriptional activation and was reversible within 12 h after infection. No dissociation was observed after infection with a virus that abnormally maintained the IFN-beta gene in a repressed state. Dissociation induced after virus infection specifically targeted the IFN-beta locus without affecting the general structure and nuclear distribution of PCH clusters. Using cell lines stably transfected with wild-type or mutated IFN-beta promoters, we identified the proximal region of the IFN-beta promoter containing YY1 DNA-binding sites as the region regulating IFN-beta promoter association with PCH before as well as during virus infection.

post to: CiteULike

R-SAP: a multi-threading computational pipeline for the characterization of high-throughput RNA-sequencing data.

Publication Date: 2012 Jan 28 PMID: 22287631
Authors: Mittal, V. K. - McDonald, J. F.
Journal: Nucleic Acids Res

The rapid expansion in the quantity and quality of RNA-Seq data requires the development of sophisticated high-performance bioinformatics tools capable of rapidly transforming this data into meaningful information that is easily interpretable by biologists. Currently available analysis tools are often not easily installed by the general biologist and most of them lack inherent parallel processing capabilities widely recognized as an essential feature of next-generation bioinformatics tools. We present here a user-friendly and fully automated RNA-Seq analysis pipeline (R-SAP) with built-in multi-threading capability to analyze and quantitate high-throughput RNA-Seq datasets. R-SAP follows a hierarchical decision making procedure to accurately characterize various classes of transcripts and achieves a near linear decrease in data processing time as a result of increased multi-threading. In addition, RNA expression level estimates obtained using R-SAP display high concordance with levels measured by microarrays.

post to: CiteULike

In vivo screening of modified siRNAs for non-specific antiviral effect in a small fish model: number and localization in the strands are important.

Publication Date: 2012 Jan 28 PMID: 22287630
Authors: Schyth, B. D. - Bramsen, J. B. - Pakula, M. M. - Larashati, S. - Kjems, J. - Wengel, J. - Lorenzen, N.
Journal: Nucleic Acids Res

Small interfering RNAs (siRNAs) are promising new active compounds in gene medicine but the induction of non-specific immune responses following their delivery continues to be a serious problem. With the purpose of avoiding such effects chemically modified siRNAs are tested in screening assay but often only examining the expression of specific immunologically relevant genes in selected cell populations typically blood cells from treated animals or humans. Assays using a relevant physiological state in biological models as read-out are not common. Here we use a fish model where the innate antiviral effect of siRNAs is functionally monitored as reduced mortality in challenge studies involving an interferon sensitive virus. Modifications with locked nucleic acid (LNA), altritol nucleic acid (ANA) and hexitol nucleic acid (HNA) reduced the antiviral protection in this model indicative of altered immunogenicity. For LNA modified siRNAs, the number and localization of modifications in the single strands was found to be important and a correlation between antiviral protection and the thermal stability of siRNAs was found. The previously published sisiRNA will in some sequences, but not all, increase the antiviral effect of siRNAs. The applied fish model represents a potent tool for conducting fast but statistically and scientifically relevant evaluations of chemically optimized siRNAs with respect to non-specific antiviral effects in vivo.

post to: CiteULike

DNA helicase and helicase-nuclease enzymes with a conserved iron-sulfur cluster.

Publication Date: 2012 Jan 28 PMID: 22287629
Authors: Wu, Y. - Brosh, R. M. Jr
Journal: Nucleic Acids Res

Conserved Iron-Sulfur (Fe-S) clusters are found in a growing family of metalloproteins that are implicated in prokaryotic and eukaryotic DNA replication and repair. Among these are DNA helicase and helicase-nuclease enzymes that preserve chromosomal stability and are genetically linked to diseases characterized by DNA repair defects and/or a poor response to replication stress. Insight to the structural and functional importance of the conserved Fe-S domain in DNA helicases has been gleaned from structural studies of the purified proteins and characterization of Fe-S cluster site-directed mutants. In this review, we will provide a current perspective of what is known about the Fe-S cluster helicases, with an emphasis on how the conserved redox active domain may facilitate mechanistic aspects of helicase function. We will discuss testable models for how the conserved Fe-S cluster might operate in helicase and helicase-nuclease enzymes to conduct their specialized functions that help to preserve the integrity of the genome.

post to: CiteULike

Structure-function analysis and genetic interactions of the yeast branchpoint binding protein Msl5.

Publication Date: 2012 Jan 28 PMID: 22287628
Authors: Chang, J. - Schwer, B. - Shuman, S.
Journal: Nucleic Acids Res

Saccharomyces cerevisiae Msl5 (branchpoint binding protein) orchestrates spliceosome assembly by binding the branchpoint sequence 5'-UACUAAC and establishing cross intron-bridging interactions with other components of the splicing machinery. Reciprocal tandem affinity purifications verify that Msl5 exists in vivo as a heterodimer with Mud2 and that the Msl5-Mud2 complex is associated with the U1 snRNP. By gauging the ability of mutants of Msl5 to complement msl5Delta, we find that the Mud2-binding (amino acids 35-54) and putative Prp40-binding (PPxY(100)) elements of the Msl5 N-terminal domain are inessential, as are the C-terminal proline-rich domain (amino acids 382-476) and two zinc-binding CxxCxxxxHxxxxC motifs (amino acids 273-286 and 299-312). A subset of conserved branchpoint RNA-binding amino acids in the central KH-QUA2 domain (amino acids 146-269) are essential pairwise (Ile198-Arg190; Leu256-Leu259) or in trios (Leu169-Arg172-Leu176), whereas other pairs of RNA-binding residues are dispensable. We used our collection of viable Msl5 mutants to interrogate synthetic genetic interactions, in cis between the inessential structural elements of the Msl5 polypeptide and in trans between Msl5 and yeast splicing factors (Mud2, Nam8 and Tgs1) that are optional for vegetative growth. The results suggest a network of important but functionally buffered protein-protein and protein-RNA interactions between the Mud2-Msl5 complex at the branchpoint and the U1 snRNP at the 5' splice site.

post to: CiteULike

Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation.

Publication Date: 2012 Jan 28 PMID: 22287627
Authors: McCarthy, D. J. - Chen, Y. - Smyth, G. K.
Journal: Nucleic Acids Res

A flexible statistical framework is developed for the analysis of read counts from RNA-Seq gene expression studies. It provides the ability to analyse complex experiments involving multiple treatment conditions and blocking variables while still taking full account of biological variation. Biological variation between RNA samples is estimated separately from the technical variation associated with sequencing technologies. Novel empirical Bayes methods allow each gene to have its own specific variability, even when there are relatively few biological replicates from which to estimate such variability. The pipeline is implemented in the edgeR package of the Bioconductor project. A case study analysis of carcinoma data demonstrates the ability of generalized linear model methods (GLMs) to detect differential expression in a paired design, and even to detect tumour-specific expression changes. The case study demonstrates the need to allow for gene-specific variability, rather than assuming a common dispersion across genes or a fixed relationship between abundance and variability. Genewise dispersions de-prioritize genes with inconsistent results and allow the main analysis to focus on changes that are consistent between biological replicates. Parallel computational approaches are developed to make non-linear model fitting faster and more reliable, making the application of GLMs to genomic data more convenient and practical. Simulations demonstrate the ability of adjusted profile likelihood estimators to return accurate estimators of biological variability in complex situations. When variation is gene-specific, empirical Bayes estimators provide an advantageous compromise between the extremes of assuming common dispersion or separate genewise dispersion. The methods developed here can also be applied to count data arising from DNA-Seq applications, including ChIP-Seq for epigenetic marks and DNA methylation analyses.

post to: CiteULike

Evolution and function of CAG/polyglutamine repeats in protein-protein interaction networks.

Publication Date: 2012 Jan 28 PMID: 22287626
Authors: Schaefer, M. H. - Wanker, E. E. - Andrade-Navarro, M. A.
Journal: Nucleic Acids Res

Expanded runs of consecutive trinucleotide CAG repeats encoding polyglutamine (polyQ) stretches are observed in the genes of a large number of patients with different genetic diseases such as Huntington's and several Ataxias. Protein aggregation, which is a key feature of most of these diseases, is thought to be triggered by these expanded polyQ sequences in disease-related proteins. However, polyQ tracts are a normal feature of many human proteins, suggesting that they have an important cellular function. To clarify the potential function of polyQ repeats in biological systems, we systematically analyzed available information stored in sequence and protein interaction databases. By integrating genomic, phylogenetic, protein interaction network and functional information, we obtained evidence that polyQ tracts in proteins stabilize protein interactions. This happens most likely through structural changes whereby the polyQ sequence extends a neighboring coiled-coil region to facilitate its interaction with a coiled-coil region in another protein. Alteration of this important biological function due to polyQ expansion results in gain of abnormal interactions, leading to pathological effects like protein aggregation. Our analyses suggest that research on polyQ proteins should shift focus from expanded polyQ proteins into the characterization of the influence of the wild-type polyQ on protein interactions.

post to: CiteULike

A novel non-homologous recombination-mediated mechanism for Escherichia coli unilateral flagellar phase variation.

Publication Date: 2012 Jan 28 PMID: 22287625
Authors: Liu, B. - Hu, B. - Zhou, Z. - Guo, D. - Guo, X. - Ding, P. - Feng, L. - Wang, L.
Journal: Nucleic Acids Res

Flagella contribute to the virulence of bacteria through chemotaxis, adhesion to and invasion of host surfaces. Flagellar phase variation is believed to facilitate bacterial evasion of the host immune response. In this study, the flnA gene that encodes Escherichia coli H17 flagellin was examined by whole genome sequencing and genetic deletion analysis. Unilateral flagellar phase variation has been reported in E. coli H3, H47 and H17 strains, although the mechanism for phase variation in the H17 strain has not been previously understood. Analysis of phase variants indicated that the flagellar phase variation in the H17 strain was caused by the deletion of an approximately 35 kb DNA region containing the flnA gene from diverse excision sites. The presence of covalently closed extrachromosomal circular forms of this excised 35 kb region was confirmed by the two-step polymerase chain reaction. The deletion and complementation test revealed that the Int1157 integrase, a tyrosine recombinase, mediates the excision of this region. Unlike most tyrosine recombinases, Int1157 is suggested to recognize diverse sites and mediate recombination between non-homologous DNA sequences. This is the first report of non-homologous recombination mediating flagellar phase variation.

post to: CiteULike

Crystal structure of a c-kit promoter quadruplex reveals the structural role of metal ions and water molecules in maintaining loop conformation.

Publication Date: 2012 Jan 28 PMID: 22287624
Authors: Wei, D. - Parkinson, G. N. - Reszka, A. P. - Neidle, S.
Journal: Nucleic Acids Res

We report here the 1.62 A crystal structure of an intramolecular quadruplex DNA formed from a sequence in the promoter region of the c-kit gene. This is the first reported crystal structure of a promoter quadruplex and the first observation of localized magnesium ions in a quadruplex structure. The structure reveals that potassium and magnesium ions have an unexpected yet significant structural role in stabilizing particular quadruplex loops and grooves that is distinct from but in addition to the role of potassium ions in the ion channel at the centre of all quadruplex structures. The analysis also shows how ions cluster together with structured water molecules to stabilize the quadruplex arrangement. This particular quadruplex has been previously studied by NMR methods, and the present X-ray structure is in accord with the earlier topology assignment. However, as well as the observations of potassium and magnesium ions, the crystal structure has revealed a highly significant difference in the dimensions of the large cleft in the structure, which is a plausible target for small molecules. This difference can be understood by the stabilizing role of structured water networks.

post to: CiteULike

RNA folding with soft constraints: reconciliation of probing data and thermodynamic secondary structure prediction.

Publication Date: 2012 Jan 28 PMID: 22287623
Authors: Washietl, S. - Hofacker, I. L. - Stadler, P. F. - Kellis, M.
Journal: Nucleic Acids Res

Thermodynamic folding algorithms and structure probing experiments are commonly used to determine the secondary structure of RNAs. Here we propose a formal framework to reconcile information from both prediction algorithms and probing experiments. The thermodynamic energy parameters are adjusted using 'pseudo-energies' to minimize the discrepancy between prediction and experiment. Our framework differs from related approaches that used pseudo-energies in several key aspects. (i) The energy model is only changed when necessary and no adjustments are made if prediction and experiment are consistent. (ii) Pseudo-energies remain biophysically interpretable and hold positional information where experiment and model disagree. (iii) The whole thermodynamic ensemble of structures is considered thus allowing to reconstruct mixtures of suboptimal structures from seemingly contradicting data. (iv) The noise of the energy model and the experimental data is explicitly modeled leading to an intuitive weighting factor through which the problem can be seen as folding with 'soft' constraints of different strength. We present an efficient algorithm to iteratively calculate pseudo-energies within this framework and demonstrate how this approach can be used in combination with SHAPE chemical probing data to improve secondary structure prediction. We further demonstrate that the pseudo-energies correlate with biophysical effects that are known to affect RNA folding such as chemical nucleotide modifications and protein binding.

post to: CiteULike

Mycobacterium smegmatis RqlH defines a novel clade of bacterial RecQ-like DNA helicases with ATP-dependent 3'-5' translocase and duplex unwinding activities.

Publication Date: 2012 Jan 28 PMID: 22287622
Authors: Ordonez, H. - Unciuleac, M. - Shuman, S.
Journal: Nucleic Acids Res

The Escherichia coli RecQ DNA helicase participates in a pathway of DNA repair that operates in parallel to the recombination pathway driven by the multisubunit helicase-nuclease machine RecBCD. The model mycobacterium Mycobacterium smegmatis executes homologous recombination in the absence of its helicase-nuclease machine AdnAB, though it lacks a homolog of E. coli RecQ. Here, we identify and characterize M. smegmatis RqlH, a RecQ-like helicase with a distinctive domain structure. The 691-amino acid RqlH polypeptide consists of a RecQ-like ATPase domain (amino acids 1-346) and tetracysteine zinc-binding domain (amino acids 435-499), separated by an RqlH-specific linker. RqlH lacks the C-terminal HRDC domain found in E. coli RecQ. Rather, the RqlH C-domain resembles bacterial ComF proteins and includes a phosphoribosyltransferase-like module. We show that RqlH is a DNA-dependent ATPase/dATPase that translocates 3'-5' on single-stranded DNA and has 3'-5' helicase activity. These functions inhere to RqlH-(1-505), a monomeric motor unit comprising the ATPase, linker and zinc-binding domains. RqlH homologs are distributed widely among bacterial taxa. The mycobacteria that encode RqlH lack a classical RecQ, though many other Actinobacteria have both RqlH and RecQ. Whereas E. coli K12 encodes RecQ but lacks a homolog of RqlH, other strains of E. coli have both RqlH and RecQ.

post to: CiteULike

Silencing of toxic gene expression by Fis.

Publication Date: 2012 Jan 28 PMID: 22287621
Authors: Karambelkar, S. - Swapna, G. - Nagaraja, V.
Journal: Nucleic Acids Res

Bacteria and bacteriophages have evolved DNA modification as a strategy to protect their genomes. Mom protein of bacteriophage Mu modifies the phage DNA, rendering it refractile to numerous restriction enzymes and in turn enabling the phage to successfully invade a variety of hosts. A strong fortification, a combined activity of the phage and host factors, prevents untimely expression of mom and associated toxic effects. Here, we identify the bacterial chromatin architectural protein Fis as an additional player in this crowded regulatory cascade. Both in vivo and in vitro studies described here indicate that Fis acts as a transcriptional repressor of mom promoter. Further, our data shows that Fis mediates its repressive effect by denying access to RNA polymerase at mom promoter. We propose that a combined repressive effect of Fis and previously characterized negative regulatory factors could be responsible to keep the gene silenced most of the time. We thus present a new facet of Fis function in Mu biology. In addition to bringing about overall downregulation of Mu genome, it also ensures silencing of the advantageous but potentially lethal mom gene.

post to: CiteULike

Regulation of Enteric vapBC Transcription: Induction by VapC Toxin Dimer-Breaking.

Publication Date: 2012 Jan 27 PMID: 22287572
Authors: Winther, K. S. - Gerdes, K.
Journal: Nucleic Acids Res

Toxin-antitoxin (TA) loci encode inhibitors of translation, replication or cell wall synthesis and are common elements of prokaryotic plasmids and chromosomes. Ten TA loci of Escherichia coli K-12 encode mRNases that cumulatively contribute to persistence (multidrug tolerance) of the bacterial cells. The mechanisms underlying induction and reversion of the persistent state are not yet understood. The vapBC operon of Salmonalla enterica serovar Typhimurium LT2 encodes VapC, a tRNase that reversibly inhibits translation by site-specific cleavage of tRNA(fMet). VapB is an antitoxin that interacts with and neutralizes VapC via its C-terminal tail and regulate TA operon transcription via its N-terminal DNA binding domain that recognize operators in the vapBC promoter region. We show here that transcription of the vapBC operon of S. enterica is controlled by a recently discovered regulatory theme referred to as 'conditional cooperativity': at low T/A ratios, the TA complex binds cooperatively to the promoter region and represses TA operon transcription whereas at high T/A ratios, the excess toxin leads to destabilization of the TA-operator complex and therefore, induction of transcription. We present evidence that an excess of VapC toxin leads to operator complex destabilization by breaking of toxin dimers.

post to: CiteULike

A human XRCC4-XLF complex bridges DNA.

Publication Date: 2012 Jan 27 PMID: 22287571
Authors: Andres, S. N. - Vergnes, A. - Ristic, D. - Wyman, C. - Modesti, M. - Junop, M.
Journal: Nucleic Acids Res

DNA double-strand breaks pose a significant threat to cell survival and must be repaired. In higher eukaryotes, such damage is repaired efficiently by non-homologous end joining (NHEJ). Within this pathway, XRCC4 and XLF fulfill key roles required for end joining. Using DNA-binding and -bridging assays, combined with direct visualization, we present evidence for how XRCC4-XLF complexes robustly bridge DNA molecules. This unanticipated, DNA Ligase IV-independent bridging activity by XRCC4-XLF suggests an early role for this complex during end joining, in addition to its more well-established later functions. Mutational analysis of the XRCC4-XLF C-terminal tail regions further identifies specialized functions in complex formation and interaction with DNA and DNA Ligase IV. Based on these data and the crystal structure of an extended protein filament of XRCC4-XLF at 3.94 A, a model for XRCC4-XLF complex function in NHEJ is presented.

post to: CiteULike

The conserved 5' apical hairpin stem loops of bamboo mosaic virus and its satellite RNA contribute to replication competence.

Publication Date: 2012 Jan 25 PMID: 22278884
Authors: Chen, H. C. - Kong, L. R. - Yeh, T. Y. - Cheng, C. P. - Hsu, Y. H. - Lin, N. S.
Journal: Nucleic Acids Res

Satellite RNAs associated with Bamboo mosaic virus (satBaMVs) depend on BaMV for replication and encapsidation. Certain satBaMVs isolated from natural fields significantly interfere with BaMV replication. The 5' apical hairpin stem loop (AHSL) of satBaMV is the major determinant in interference with BaMV replication. In this study, by in vivo competition assay, we revealed that the sequence and structure of AHSL, along with specific nucleotides (C(60) and C(83)) required for interference with BaMV replication, are also involved in replication competition among satBaMV variants. Moreover, all of the 5' ends of natural BaMV isolates contain the similar AHSLs having conserved nucleotides (C(64) and C(86)) with those of interfering satBaMVs, suggesting their co-evolution. Mutational analyses revealed that C(86) was essential for BaMV replication, and that replacement of C(64) with U reduced replication efficiency. The non-interfering satBaMV interfered with BaMV replication with the BaMV-C64U mutant as helper. These findings suggest that two cytosines at the equivalent positions in the AHSLs of BaMV and satBaMV play a crucial role in replication competence. The downregulation level, which is dependent upon the molar ratio of interfering satBaMV to BaMV, implies that there is competition for limited replication machinery.

post to: CiteULike

Sample richness and genetic diversity as drivers of chimera formation in nSSU metagenetic analyses.

Publication Date: 2012 Jan 25 PMID: 22278883
Authors: Fonseca, V. G. - Nichols, B. - Lallias, D. - Quince, C. - Carvalho, G. R. - Power, D. M. - Creer, S.
Journal: Nucleic Acids Res

Eukaryotic diversity in environmental samples is often assessed via PCR-based amplification of nSSU genes. However, estimates of diversity derived from pyrosequencing environmental data sets are often inflated, mainly because of the formation of chimeric sequences during PCR amplification. Chimeras are hybrid products composed of distinct parental sequences that can lead to the misinterpretation of diversity estimates. We have analyzed the effect of sample richness, evenness and phylogenetic diversity on the formation of chimeras using a nSSU data set derived from 454 Roche pyrosequencing of replicated, large control pools of closely and distantly related nematode mock communities, of known intragenomic identity and richness. To further investigate how chimeric molecules are formed, the nSSU gene secondary structure was analyzed in several individuals. For the first time in eukaryotes, chimera formation proved to be higher in both richer and more genetically diverse samples, thus providing a novel perspective of chimera formation in pyrosequenced environmental data sets. Findings contribute to a better understanding of the nature and mechanisms involved in chimera formation during PCR amplification of environmentally derived DNA. Moreover, given the similarities between biodiversity analyses using amplicon sequencing and those used to assess genomic variation, our findings have potential broad application for identifying genetic variation in homologous loci or multigene families in general.

post to: CiteULike

DNMT1 modulates gene expression without its catalytic activity partially through its interactions with histone-modifying enzymes.

Publication Date: 2012 Jan 25 PMID: 22278882
Authors: Clements, E. G. - Mohammad, H. P. - Leadem, B. R. - Easwaran, H. - Cai, Y. - Van Neste, L. - Baylin, S. B.
Journal: Nucleic Acids Res

While DNA methyltransferase1 (DNMT1) is classically known for its functions as a maintenance methyltransferase enzyme, additional roles for DNMT1 in gene expression are not as clearly understood. Several groups have shown that deletion of the catalytic domain from DNMT1 does not abolish repressive activity of the protein against a reporter gene. In our studies, we examine the repressor function of catalytically inactive DNMT1 at endogenous genes. First, potential DNMT1 target genes were identified by searching for genes up-regulated in HCT116 colon cancer cells genetically disrupted for DNMT1 (DNMT1(-/-) hypomorph cells). Next, the requirement for DNMT1 activity for repression of these genes was assessed by stably restoring expression of wild-type or catalytically inactive DNMT1. Both wild-type and mutant proteins are able to occupy the promoters and repress the expression of a set of target genes, and induce, at these promoters, both the depletion of active histone marks and the recruitment of a H3K4 demethylase, KDM1A/LSD1. Together, our findings show that there are genes for which DNMT1 acts as a transcriptional repressor independent from its methyltransferase function and that this repressive function may invoke a role for a scaffolding function of the protein at target genes.

post to: CiteULike

Amplified microRNA detection by templated chemistry.

Publication Date: 2012 Jan 25 PMID: 22278881
Authors: Harcourt, E. M. - Kool, E. T.
Journal: Nucleic Acids Res

MicroRNAs (miRNAs) are a class of RNAs that play important regulatory roles in the cell. The detection of microRNA has attracted significant interest recently, as abnormal miRNA expression has been linked to cancer and other diseases. Here, we present a straightforward method for isothermal amplified detection of miRNA that involves two separate nucleic acid-templated chemistry steps. The miRNA first templates the cyclization of an oligodeoxynucleotide from a linear precursor containing a 5'-iodide and a 3'-phosphorothioate. The sequence is amplified through rolling circle amplification with 29 DNA polymerase and then detected via a second amplification using fluorogenic templated probes. Tests showed that the cyclization proceeds in approximately 50% yield over 24 h and is compatible with the conditions required for rolling circle polymerization, unlike enzymatic ligations which required non-compatible buffer conditions. The polymerization yielded 188-fold amplification, and separate experiments showed approximately 15-fold signal amplification from the templated fluorogenic probes. When all components are combined, results show miRNA detection down to 200 pM in solution, and correlation of the detected signal with the initial concentration of miRNA. The doubly templated double-amplification method demonstrates a new approach to detection of rolling circle products and significant advantages in ease of operation for miRNA detection.

post to: CiteULike

The nucleoid-associated proteins H-NS and FIS modulate the DNA supercoiling response of the pel genes, the major virulence factors in the plant pathogen bacterium Dickeya dadantii.

Publication Date: 2012 Jan 24 PMID: 22275524
Authors: Ouafa, Z. A. - Reverchon, S. - Lautier, T. - Muskhelishvili, G. - Nasser, W.
Journal: Nucleic Acids Res

Dickeya dadantii is a pathogen infecting a wide range of plant species. Soft rot, the visible symptom, is mainly due to the production of pectate lyases (Pels) that can destroy the plant cell walls. Previously we found that the pel gene expression is modulated by H-NS and FIS, two nucleoid-associated proteins (NAPs) modulating the DNA topology. Here, we show that relaxation of the DNA in growing D. dadantii cells decreases the expression of pel genes. Deletion of fis aggravates, whereas that of hns alleviates the negative impact of DNA relaxation on pel expression. We further show that H-NS and FIS directly bind the pelE promoter and that the response of D. dadantii pel genes to stresses that induce DNA relaxation is modulated, although to different extents, by H-NS and FIS. We infer that FIS acts as a repressor buffering the negative impact of DNA relaxation on pel gene transcription, whereas H-NS fine-tunes the response of virulence genes precluding their expression under suboptimal conditions of supercoiling. This novel dependence of H-NS effect on DNA topology expands our understanding of the role of NAPs in regulating the global bacterial gene expression and bacterial pathogenicity.

post to: CiteULike

Global regulation of gene expression by OxyR in an important human opportunistic pathogen.

Publication Date: 2012 Jan 24 PMID: 22275523
Authors: Wei, Q. - Le Minh, P. N. - Dotsch, A. - Hildebrand, F. - Panmanee, W. - Elfarash, A. - Schulz, S. - Plaisance, S. - Charlier, D. - Hassett, D. - Haussler, S. - Cornelis, P.
Journal: Nucleic Acids Res

Most bacteria control oxidative stress through the H(2)O(2)-responsive transactivator OxyR, a member of the LTTR family (LysR Type Transcriptional Regulators), which activates the expression of defensive genes such as those encoding catalases, alkyl hydroperoxide reductases and superoxide dismutases. In the human opportunistic pathogen Pseudomonas aeruginosa, OxyR positively regulates expression of the oxidative stress response genes katA, katB, ahpB and ahpCF. To identify additional targets of OxyR in P. aeruginosa PAO1, we performed chromatin immunoprecipitation in combination with whole genome tiling array analyses (ChIP-chip). We detected 56 genes including all the previously identified defensive genes and a battery of novel direct targets of OxyR. Electrophoretic mobility shift assays (EMSAs) for selected newly identified targets indicated that approximately 70% of those were bound by purified oxidized OxyR and their regulation was confirmed by quantitative real-time polymerase chain reaction. Furthermore, a thioredoxin system was identified to enzymatically reduce OxyR under oxidative stress. Functional classification analysis showed that OxyR controls a core regulon of oxidative stress defensive genes, and other genes involved in regulation of iron homeostasis (pvdS), quorum-sensing (rsaL), protein synthesis (rpsL) and oxidative phosphorylation (cyoA and snr1). Collectively, our results indicate that OxyR is involved in oxidative stress defense and regulates other aspects of cellular metabolism as well.

post to: CiteULike

P1 and P2 protein heterodimer binding to the P0 protein of Saccharomyces cerevisiae is relatively non-specific and a source of ribosomal heterogeneity.

Publication Date: 2012 Jan 24 PMID: 22275522
Authors: Cardenas, D. - Revuelta-Cervantes, J. - Jimenez-Diaz, A. - Camargo, H. - Remacha, M. - Ballesta, J. P.
Journal: Nucleic Acids Res

The ribosomal stalk is formed by four acidic phosphoproteins in Saccharomyces cerevisiae, P1alpha, P1beta, P2alpha and P2beta, which form two heterodimers, P1alpha/P2beta and P1beta/P2alpha, that preferentially bind to sites A and B of the P0 protein, respectively. Using mutant strains carrying only one of the four possible P1/P2 combinations, we found a specific phenotype associated to each P1/P2 pair, indicating that not all acidic P proteins play the same role. The absence of one P1/P2 heterodimer reduced the rate of cell growth by varying degrees, depending on the proteins missing. Synthesis of the 60S ribosomal subunit also decreased, particularly in strains carrying the unusual P1alpha-P2alpha or P1beta-P2beta heterodimers, although the distinct P1/P2 dimers are bound with similar affinity to the mutant ribosome. While in wild-type strains the B site bound P1beta/P2alpha in a highly specific manner and the A site bound the four P proteins similarly, both the A and B binding sites efficiently bound practically any P1/P2 pair in mutant strains expressing truncated P0 proteins. The reported results support that while most ribosomes contain a P1alpha/P2beta-P0-P1beta/P2alpha structure in normal conditions, the stalk assembly mechanism can generate alternative compositions, which have been previously detected in the cell.

post to: CiteULike

Direct and specific chemical control of eukaryotic translation with a synthetic RNA-protein interaction.

Publication Date: 2012 Jan 24 PMID: 22275521
Authors: Goldfless, S. J. - Belmont, B. J. - de Paz, A. M. - Liu, J. F. - Niles, J. C.
Journal: Nucleic Acids Res

Sequence-specific RNA-protein interactions, though commonly used in biological systems to regulate translation, are challenging to selectively modulate. Here, we demonstrate the use of a chemically-inducible RNA-protein interaction to regulate eukaryotic translation. By genetically encoding Tet Repressor protein (TetR)-binding RNA elements into the 5'-untranslated region (5'-UTR) of an mRNA, translation of a downstream coding sequence is directly controlled by TetR and tetracycline analogs. In endogenous and synthetic 5'-UTR contexts, this system efficiently regulates the expression of multiple target genes, and is sufficiently stringent to distinguish functional from non-functional RNA-TetR interactions. Using a reverse TetR variant, we illustrate the potential for expanding the regulatory properties of the system through protein engineering strategies.

post to: CiteULike

Flp and Cre expressed from Flp-2A-Cre and Flp-IRES-Cre transcription units mediate the highest level of dual recombinase-mediated cassette exchange.

Publication Date: 2012 Jan 23 PMID: 22270085
Authors: Anderson, R. P. - Voziyanova, E. - Voziyanov, Y.
Journal: Nucleic Acids Res

Recombinase-mediated cassette exchange (RMCE) is a powerful tool for unidirectional integration of DNA fragments of interest into a pre-determined genome locale. In this report, we examined how the efficiency of dual RMCE catalyzed by Flp and Cre depends on the nature of transcription units that express the recombinases. The following recombinase transcription units were analyzed: (i) Flp and Cre genes expressed as individual transcription units located on different vectors, (ii) Flp and Cre genes expressed as individual transcription units located on the same vector, (iii) Flp and Cre genes expressed from a single promoter and separated by internal ribosome entry sequence and (iv) Flp and Cre coding sequences separated by the 2A peptide and expressed as a single gene. We found that the highest level of dual RMCE (35-45% of the transfected cells) can be achieved when Flp and Cre recombinases are expressed as Flp-2A-Cre and Flp-IRES-Cre transcription units. In contrast, the lowest level of dual RMCE ( approximately 1% of the transfected cells) is achieved when Flp and Cre are expressed as individual transcription units. The analysis shows that it is the relative Flp-to-Cre ratio that critically affects the efficiency of dual RMCE. Our results will be helpful for maximizing the efficiency of dual RMCE aimed to engineer and re-engineer genomes.

post to: CiteULike

Biogenesis of mammalian microRNAs by a non-canonical processing pathway.

Publication Date: 2012 Jan 23 PMID: 22270084
Authors: Havens, M. A. - Reich, A. A. - Duelli, D. M. - Hastings, M. L.
Journal: Nucleic Acids Res

Canonical microRNA biogenesis requires the Microprocessor components, Drosha and DGCR8, to generate precursor-miRNA, and Dicer to form mature miRNA. The Microprocessor is not required for processing of some miRNAs, including mirtrons, in which spliceosome-excised introns are direct Dicer substrates. In this study, we examine the processing of putative human mirtrons and demonstrate that although some are splicing-dependent, as expected, the predicted mirtrons, miR-1225 and miR-1228, are produced in the absence of splicing. Remarkably, knockout cell lines and knockdown experiments demonstrated that biogenesis of these splicing-independent mirtron-like miRNAs, termed 'simtrons', does not require the canonical miRNA biogenesis components, DGCR8, Dicer, Exportin-5 or Argonaute 2. However, simtron biogenesis was reduced by expression of a dominant negative form of Drosha. Simtrons are bound by Drosha and processed in vitro in a Drosha-dependent manner. Both simtrons and mirtrons function in silencing of target transcripts and are found in the RISC complex as demonstrated by their interaction with Argonaute proteins. These findings reveal a non-canonical miRNA biogenesis pathway that can produce functional regulatory RNAs.

post to: CiteULike

On the assessment of statistical significance of three-dimensional colocalization of sets of genomic elements.

Publication Date: 2012 Jan 20 PMID: 22266657
Authors: Witten, D. M. - Noble, W. S.
Journal: Nucleic Acids Res

A growing body of experimental evidence supports the hypothesis that the 3D structure of chromatin in the nucleus is closely linked to important functional processes, including DNA replication and gene regulation. In support of this hypothesis, several research groups have examined sets of functionally associated genomic loci, with the aim of determining whether those loci are statistically significantly colocalized. This work presents a critical assessment of two previously reported analyses, both of which used genome-wide DNA-DNA interaction data from the yeast Saccharomyces cerevisiae, and both of which rely upon a simple notion of the statistical significance of colocalization. We show that these previous analyses rely upon a faulty assumption, and we propose a correct non-parametric resampling approach to the same problem. Applying this approach to the same data set does not support the hypothesis that transcriptionally coregulated genes tend to colocalize, but strongly supports the colocalization of centromeres, and provides some evidence of colocalization of origins of early DNA replication, chromosomal breakpoints and transfer RNAs.

post to: CiteULike

MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data.

Publication Date: 2012 Feb 1 PMID: 22266656
Authors: Shen, S. - Won Park, J. - Huang, J. - Dittmar, K. A. - Lu, Z. X. - Zhou, Q. - Carstens, R. P. - Xing, Y.
Journal: Nucleic Acids Res

Ultra-deep RNA sequencing has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We develop MATS (multivariate analysis of transcript splicing), a Bayesian statistical framework for flexible hypothesis testing of differential alternative splicing patterns on RNA-Seq data. MATS uses a multivariate uniform prior to model the between-sample correlation in exon splicing patterns, and a Markov chain Monte Carlo (MCMC) method coupled with a simulation-based adaptive sampling procedure to calculate the P-value and false discovery rate (FDR) of differential alternative splicing. Importantly, the MATS approach is applicable to almost any type of null hypotheses of interest, providing the flexibility to identify differential alternative splicing events that match a given user-defined pattern. We evaluated the performance of MATS using simulated and real RNA-Seq data sets. In the RNA-Seq analysis of alternative splicing events regulated by the epithelial-specific splicing factor ESRP1, we obtained a high RT-PCR validation rate of 86% for differential exon skipping events with a MATS FDR of <10%. Additionally, over the full list of RT-PCR tested exons, the MATS FDR estimates matched well with the experimental validation rate. Our results demonstrate that MATS is an effective and flexible approach for detecting differential alternative splicing from RNA-Seq data.

post to: CiteULike

Revealing stable processing products from ribosome-associated small RNAs by deep-sequencing data analysis.

Publication Date: 2012 Jan 20 PMID: 22266655
Authors: Zywicki, M. - Bakowska-Zywicka, K. - Polacek, N.
Journal: Nucleic Acids Res

The exploration of the non-protein-coding RNA (ncRNA) transcriptome is currently focused on profiling of microRNA expression and detection of novel ncRNA transcription units. However, recent studies suggest that RNA processing can be a multi-layer process leading to the generation of ncRNAs of diverse functions from a single primary transcript. Up to date no methodology has been presented to distinguish stable functional RNA species from rapidly degraded side products of nucleases. Thus the correct assessment of widespread RNA processing events is one of the major obstacles in transcriptome research. Here, we present a novel automated computational pipeline, named APART, providing a complete workflow for the reliable detection of RNA processing products from next-generation-sequencing data. The major features include efficient handling of non-unique reads, detection of novel stable ncRNA transcripts and processing products and annotation of known transcripts based on multiple sources of information. To disclose the potential of APART, we have analyzed a cDNA library derived from small ribosome-associated RNAs in Saccharomyces cerevisiae. By employing the APART pipeline, we were able to detect and confirm by independent experimental methods multiple novel stable RNA molecules differentially processed from well known ncRNAs, like rRNAs, tRNAs or snoRNAs, in a stress-dependent manner.

post to: CiteULike

ATM regulates proteasome-dependent subnuclear localization of TRF1, which is important for telomere maintenance.

Publication Date: 2012 Jan 20 PMID: 22266654
Authors: McKerlie, M. A. - Lin, S. - Zhu, X. D.
Journal: Nucleic Acids Res

Ataxia telangiectasia mutated (ATM), a PI-3 kinase essential for maintaining genomic stability, has been shown to regulate TRF1, a negative mediator of telomerase-dependent telomere extension. However, little is known about ATM-mediated TRF1 phosphorylation site(s) in vivo. Here, we report that ATM phosphorylates S367 of TRF1 and that this phosphorylation renders TRF1 free of chromatin. We show that phosphorylated (pS367)TRF1 forms distinct non-telomeric subnuclear foci and that these foci occur predominantly in S and G2 phases, implying that their formation is cell cycle regulated. We show that phosphorylated (pS367)TRF1-containing foci are sensitive to proteasome inhibition. We find that a phosphomimic mutation of S367D abrogates TRF1 binding to telomeric DNA and renders TRF1 susceptible to protein degradation. In addition, we demonstrate that overexpressed TRF1-S367D accumulates in the subnuclear domains containing phosphorylated (pS367)TRF1 and that these subnuclear domains overlap with nuclear proteasome centers. Taken together, these results suggest that phosphorylated (pS367)TRF1-containing foci may represent nuclear sites for TRF1 proteolysis. Furthermore, we show that TRF1 carrying the S367D mutation is unable to inhibit telomerase-dependent telomere lengthening or to suppress the formation of telomere doublets and telomere loss in TRF1-depleted cells, suggesting that S367 phosphorylation by ATM is important for the regulation of telomere length and stability.

post to: CiteULike

The plasticity of WDR5 peptide-binding cleft enables the binding of the SET1 family of histone methyltransferases.

Publication Date: 2012 Jan 20 PMID: 22266653
Authors: Zhang, P. - Lee, H. - Brunzelle, J. S. - Couture, J. F.
Journal: Nucleic Acids Res

In mammals, the SET1 family of lysine methyltransferases (KMTs), which includes MLL1-5, SET1A and SET1B, catalyzes the methylation of lysine-4 (Lys-4) on histone H3. Recent reports have demonstrated that a three-subunit complex composed of WD-repeat protein-5 (WDR5), retinoblastoma-binding protein-5 (RbBP5) and absent, small, homeotic disks-2-like (ASH2L) stimulates the methyltransferase activity of MLL1. On the basis of studies showing that this stimulation is in part controlled by an interaction between WDR5 and a small region located in close proximity of the MLL1 catalytic domain [referred to as the WDR5-interacting motif (Win)], it has been suggested that WDR5 might play an analogous role in scaffolding the other SET1 complexes. We herein provide biochemical and structural evidence showing that WDR5 binds the Win motifs of MLL2-4, SET1A and SET1B. Comparative analysis of WDR5-Win complexes reveals that binding of the Win motifs is achieved by the plasticity of WDR5 peptidyl-arginine-binding cleft allowing the C-terminal ends of the Win motifs to be maintained in structurally divergent conformations. Consistently, enzymatic assays reveal that WDR5 plays an important role in the optimal stimulation of MLL2-4, SET1A and SET1B methyltransferase activity by the RbBP5-ASH2L heterodimer. Overall, our findings illustrate the function of WDR5 in scaffolding the SET1 family of KMTs and further emphasize on the important role of WDR5 in regulating global histone H3 Lys-4 methylation.

post to: CiteULike

A computational pipeline to discover highly phylogenetically informative genes in sequenced genomes: application to Saccharomyces cerevisiae natural strains.

Publication Date: 2012 Jan 20 PMID: 22266652
Authors: Ramazzotti, M. - Berna, L. - Stefanini, I. - Cavalieri, D.
Journal: Nucleic Acids Res

The quest for genes representing genetic relationships of strains or individuals within populations and their evolutionary history is acquiring a novel dimension of complexity with the advancement of next-generation sequencing (NGS) technologies. In fact, sequencing an entire genome uncovers genetic variation in coding and non-coding regions and offers the possibility of studying Saccharomyces cerevisiae populations at the strain level. Nevertheless, the disadvantageous cost-benefit ratio (the amount of details disclosed by NGS against the time-expensive and expertise-demanding data assembly process) still precludes the application of these techniques to the routinely assignment of yeast strains, making the selection of the most reliable molecular markers greatly desirable. In this work we propose an original computational approach to discover genes that can be used as a descriptor of the population structure. We found 13 genes whose variability can be used to recapitulate the phylogeny obtained from genome-wide sequences. The same approach that we prove to be successful in yeasts can be generalized to any other population of individuals given the availability of high-quality genomic sequences and of a clear population structure to be targeted.

post to: CiteULike

The porphyrin TmPyP4 unfolds the extremely stable G-quadruplex in MT3-MMP mRNA and alleviates its repressive effect to enhance translation in eukaryotic cells.

Publication Date: 2012 Jan 20 PMID: 22266651
Authors: Morris, M. J. - Wingate, K. L. - Silwal, J. - Leeper, T. C. - Basu, S.
Journal: Nucleic Acids Res

We report that the cationic porphyrin TmPyP4, which is known mainly as a DNA G-quadruplex stabilizer, unfolds an unusually stable all purine RNA G-quadruplex (M3Q) that is located in the 5'-UTR of MT3-MMP mRNA. When the interaction between TmPyP4 and M3Q was monitored by UV spectroscopy a 22-nm bathochromic shift and 75% hypochromicity of the porphin major Soret band was observed indicating direct binding of the two molecules. TmPyP4 disrupts folded M3Q in a concentration-dependent fashion as was observed by circular dichroism (CD), 1D (1)H NMR and native gel electrophoresis. Additionally, when TmPyP4 is present during the folding process it inhibits the M3Q RNA from adopting a G-quadruplex structure. Using a dual reporter gene construct that contained the M3Q sequence alone or the entire 5'-UTR of MT3-MMP mRNA, we report here that TmPyP4 can relieve the inhibitory effect of the M3Q G-quadruplex. However, the same concentrations of TmPyP4 failed to affect translation of a mutated construct. Thus, TmPyP4 has the ability to unfold an RNA G-quadruplex of extreme stability and modulate activity of a reporter gene presumably via the disruption of the G-quadruplex.

post to: CiteULike

Effects of Friedreich's ataxia GAA repeats on DNA replication in mammalian cells.

Publication Date: 2012 Jan 19 PMID: 22262734
Authors: Chandok, G. S. - Patel, M. P. - Mirkin, S. M. - Krasilnikova, M. M.
Journal: Nucleic Acids Res

Friedreich's ataxia (FRDA) is a common hereditary degenerative neuro-muscular disorder caused by expansions of the (GAA)n repeat in the first intron of the frataxin gene. The expanded repeats from parents frequently undergo further significant length changes as they are passed on to progeny. Expanded repeats also show an age-dependent instability in somatic cells, albeit on a smaller scale than during intergenerational transmissions. Here we studied the effects of (GAA)n repeats of varying lengths and orientations on the episomal DNA replication in mammalian cells. We have recently shown that the very first round of the transfected DNA replication occurs in the lack of the mature chromatin, does not depend on the episomal replication origin and initiates at multiple single-stranded regions of plasmid DNA. We now found that expanded GAA repeats severely block this first replication round post plasmid transfection, while the subsequent replication cycles are only mildly affected. The fact that GAA repeats affect various replication modes in a different way might shed light on their differential expansions characteristic for FRDA.

post to: CiteULike

Comprehensive literature review and statistical considerations for microarray meta-analysis.

Publication Date: 2012 Jan 19 PMID: 22262733
Authors: Tseng, G. C. - Ghosh, D. - Feingold, E.
Journal: Nucleic Acids Res

With the rapid advances of various high-throughput technologies, generation of '-omics' data is commonplace in almost every biomedical field. Effective data management and analytical approaches are essential to fully decipher the biological knowledge contained in the tremendous amount of experimental data. Meta-analysis, a set of statistical tools for combining multiple studies of a related hypothesis, has become popular in genomic research. Here, we perform a systematic search from PubMed and manual collection to obtain 620 genomic meta-analysis papers, of which 333 microarray meta-analysis papers are summarized as the basis of this paper and the other 249 GWAS meta-analysis papers are discussed in the next companion paper. The review in the present paper focuses on various biological purposes of microarray meta-analysis, databases and software and related statistical procedures. Statistical considerations of such an analysis are further scrutinized and illustrated by a case study. Finally, several open questions are listed and discussed.

post to: CiteULike

A powerful test for multiple rare variants association studies that incorporates sequencing qualities.

Publication Date: 2012 Jan 19 PMID: 22262732
Authors: Daye, Z. J. - Li, H. - Wei, Z.
Journal: Nucleic Acids Res

Next-generation sequencing data will soon become routinely available for association studies between complex traits and rare variants. Sequencing data, however, are characterized by the presence of sequencing errors at each individual genotype. This makes it especially challenging to perform association studies of rare variants, which, due to their low minor allele frequencies, can be easily perturbed by genotype errors. In this article, we develop the quality-weighted multivariate score association test (qMSAT), a new procedure that allows powerful association tests between complex traits and multiple rare variants under the presence of sequencing errors. Simulation results based on quality scores from real data show that the qMSAT often dominates over current methods, that do not utilize quality information. In particular, the qMSAT can dramatically increase power over existing methods under moderate sample sizes and relatively low coverage. Moreover, in an obesity data study, we identified using the qMSAT two functional regions (MGLL promoter and MGLL 3'-untranslated region) where rare variants are associated with extreme obesity. Due to the high cost of sequencing data, the qMSAT is especially valuable for large-scale studies involving rare variants, as it can potentially increase power without additional experimental cost. qMSAT is freely available at http://qmsat.sourceforge.net/.

post to: CiteULike

Integration of DNA into bacterial chromosomes from plasmids without a counter-selection marker.

Publication Date: 2012 Jan 18 PMID: 22259038
Authors: Heap, J. T. - Ehsaan, M. - Cooksley, C. M. - Ng, Y. K. - Cartman, S. T. - Winzer, K. - Minton, N. P.
Journal: Nucleic Acids Res

Most bacteria can only be transformed with circular plasmids, so robust DNA integration methods for these rely upon selection of single-crossover clones followed by counter-selection of double-crossover clones. To overcome the limited availability of heterologous counter-selection markers, here we explore novel DNA integration strategies that do not employ them, and instead exploit (i) activation or inactivation of genes leading to a selectable phenotype, and (ii) asymmetrical regions of homology to control the order of recombination events. We focus here on the industrial biofuel-producing bacterium Clostridium acetobutylicum, which previously lacked robust integration tools, but the approach we have developed is broadly applicable. Large sequences can be delivered in a series of steps, as we demonstrate by inserting the chromosome of phage lambda (minus a region apparently unstable in Escherichia coli in our cloning context) into the chromosome of C. acetobutylicum in three steps. This work should open the way to reliable integration of DNA including large synthetic constructs in diverse microorganisms.

post to: CiteULike

A potent 2'-O-methylated RNA-based microRNA inhibitor with unique secondary structures.

Publication Date: 2012 Jan 17 PMID: 22259037
Authors: Haraguchi, T. - Nakano, H. - Tagawa, T. - Ohki, T. - Ueno, Y. - Yoshida, T. - Iba, H.
Journal: Nucleic Acids Res

MicroRNAs (miRNAs) are involved in various biological processes and human diseases. The development of strong low-molecular weight inhibitors of specific miRNAs is thus expected to be useful in providing tools for basic research or in generating promising new therapeutic drugs. We have previously described the development of 'Tough Decoy (TuD) RNA' molecules, which achieve the long-term suppression of specific miRNA activity in mammalian cells when expressed from a lentivirus vector. In our current study, we describe new synthetic miRNA inhibitors, designated as S-TuD (Synthetic TuD), which are composed of two fully 2'-O-methylated RNA strands. Each of these strands includes a miRNA-binding site. Following the hybridization of paired strands, the resultant S-TuD forms a secondary structure with two stems, which resembles the corresponding TuD RNA molecule. By analyzing the effects of S-TuD against miR-21, miR-200c, miR-16 and miR-106b, we have elucidated the critical design features of S-TuD molecules that will provide optimum inhibitory effects following transfection into human cell lines. We further show that the inhibitory effects of a single transfection of S-TuD-miR200c are quite long-lasting (>7 days) and induce partial EMT, the full establishment of which requires 11 days when using a lentivirus vector that expresses TuD-miR200c continuously.

post to: CiteULike

Analysis of alternative splicing of cassette exons at single-cell level using two fluorescent proteins.

Publication Date: 2012 Jan 17 PMID: 22259036
Authors: Gurskaya, N. G. - Staroverov, D. B. - Zhang, L. - Fradkov, A. F. - Markina, N. M. - Pereverzev, A. P. - Lukyanov, K. A.
Journal: Nucleic Acids Res

Alternative splicing plays a major role in increasing proteome complexity and regulating gene expression. Here, we developed a new fluorescent protein-based approach to quantitatively analyze the alternative splicing of a target cassette exon (skipping or inclusion), which results in an open-reading frame shift. A fragment of a gene of interest is cloned between red and green fluorescent protein (RFP and GFP)-encoding sequences in such a way that translation of the normally spliced full-length transcript results in expression of both RFP and GFP. In contrast, alternative exon skipping results in the synthesis of RFP only. Green and red fluorescence intensities can be used to estimate the proportions of normal and alternative transcripts in each cell. The new method was successfully tested for human PIG3 (p53-inducible gene 3) cassette exon 4. Expected pattern of alternative splicing of PIG3 minigene was observed, including previously characterized effects of UV light irradiation and specific mutations. Interestingly, we observed a broad distribution of normal to alternative transcript ratio in individual cells with at least two distinct populations with approximately 45% and >95% alternative transcript. We believe that this method is useful for fluorescence-based quantitative analysis of alternative splicing of target genes in a variety of biological models.

post to: CiteULike

mRNA knockdown by single strand RNA is improved by chemical modifications.

Publication Date: 2012 Jan 16 PMID: 22253019
Authors: Haringsma, H. J. - Li, J. J. - Soriano, F. - Kenski, D. M. - Flanagan, W. M. - Willingham, A. T.
Journal: Nucleic Acids Res

While RNAi has traditionally relied on RNA duplexes, early evaluation of siRNAs demonstrated activity of the guide strand in the absence of the passenger strand. However, these single strands lacked the activity of duplex RNAs. Here, we report the systematic use of chemical modifications to optimize single-strand RNA (ssRNA)-mediated mRNA knockdown. We identify that 2'F ribose modifications coupled with 5'-end phosphorylation vastly improves ssRNA activity both in vitro and in vivo. The impact of specific chemical modifications on ssRNA activity implies an Ago-mediated mechanism but the hallmark mRNA cleavage sites were not observed which suggests ssRNA may operate through a mechanism beyond conventional Ago2 slicer activity. While currently less potent than duplex siRNAs, with additional chemical optimization and alternative routes of delivery, chemically modified ssRNAs could represent a powerful RNAi platform.

post to: CiteULike

Complex activities of the human Bloom's syndrome helicase are encoded in a core region comprising the RecA and Zn-binding domains.

Publication Date: 2012 Jan 16 PMID: 22253018
Authors: Gyimesi, M. - Harami, G. M. - Sarlos, K. - Hazai, E. - Bikadi, Z. - Kovacs, M.
Journal: Nucleic Acids Res

Bloom's syndrome DNA helicase (BLM), a member of the RecQ family, is a key player in homologous recombination (HR)-based error-free DNA repair processes. During HR, BLM exerts various biochemical activities including single-stranded (ss) DNA translocation, separation and annealing of complementary DNA strands, disruption of complex DNA structures (e.g. displacement loops) and contributes to quality control of HR via clearance of Rad51 nucleoprotein filaments. We performed a quantitative mechanistic analysis of truncated BLM constructs that are shorter than the previously identified minimal functional module. Surprisingly, we found that a BLM construct comprising only the two conserved RecA domains and the Zn(2+)-binding domain (residues 642-1077) can efficiently perform all mentioned HR-related activities. The results demonstrate that the Zn(2+)-binding domain is necessary for functional interaction with DNA. We show that the extensions of this core, including the winged-helix domain and the strand separation hairpin identified therein in other RecQ-family helicases, are not required for mechanochemical activity per se and may instead play modulatory roles and mediate protein-protein interactions.

post to: CiteULike

The rotation-coupled sliding of EcoRV.

Publication Date: 2012 Jan 12 PMID: 22241781
Authors: Dikic, J. - Menges, C. - Clarke, S. - Kokkinidis, M. - Pingoud, A. - Wende, W. - Desbiolles, P.
Journal: Nucleic Acids Res

It has been proposed that certain type II restriction enzymes (REs), such as EcoRV, track the helical pitch of DNA as they diffuse along DNA, a so-called rotation-coupled sliding. As of yet, there is no direct experimental observation of this phenomenon, but mounting indirect evidence gained from single-molecule imaging of RE-DNA complexes support the hypothesis. We address this issue by conjugating fluorescent labels of varying size (organic dyes, proteins and quantum dots) to EcoRV, and by fusing it to the engineered Rop protein scRM6. Single-molecule imaging of these modified EcoRVs sliding along DNA provides us with their linear diffusion constant (D(1)), revealing a significant size dependency. To account for the dependence of D(1) on the size of the EcoRV label, we have developed four theoretical models describing different types of motion along DNA and find that our experimental results are best described by rotation-coupled sliding of the protein. The similarity of EcoRV to other type II REs and DNA binding proteins suggests that this type of motion could be widely preserved in other biological contexts.

post to: CiteULike

A comprehensive framework for prioritizing variants in exome sequencing studies of Mendelian diseases.

Publication Date: 2012 Jan 12 PMID: 22241780
Authors: Li, M. X. - Gui, H. S. - Kwan, J. S. - Bao, S. Y. - Sham, P. C.
Journal: Nucleic Acids Res

Exome sequencing strategy is promising for finding novel mutations of human monogenic disorders. However, pinpointing the casual mutation in a small number of samples is still a big challenge. Here, we propose a three-level filtration and prioritization framework to identify the casual mutation(s) in exome sequencing studies. This efficient and comprehensive framework successfully narrowed down whole exome variants to very small numbers of candidate variants in the proof-of-concept examples. The proposed framework, implemented in a user-friendly software package, named KGGSeq (http://statgenpro.psychiatry.hku.hk/kggseq), will play a very useful role in exome sequencing-based discovery of human Mendelian disease genes.

post to: CiteULike

End-to-end attraction of duplex DNA.

Publication Date: 2012 Jan 12 PMID: 22241779
Authors: Maffeo, C. - Luan, B. - Aksimentiev, A.
Journal: Nucleic Acids Res

Recent experiments [Nakata, M. et al., End-to-end stacking and liquid crystal condensation of 6 to 20 basepair DNA duplexes. Science 2007; 318:1276-1279] have demonstrated spontaneous end-to-end association of short duplex DNA fragments into long rod-like structures. By means of extensive all-atom molecular dynamic simulations, we characterized end-to-end interactions of duplex DNA, quantitatively describing the forces, free energy and kinetics of the end-to-end association process. We found short DNA duplexes to spontaneously aggregate end-to-end when axially aligned in a small volume of monovalent electrolyte. It was observed that electrostatic repulsion of 5'-phosphoryl groups promoted the formation of aggregates in a conformation similar to the B-form DNA double helix. Application of an external force revealed that rupture of the end-to-end assembly occurs by the shearing of the terminal base pairs. The standard binding free energy and the kinetic rates of end-to-end association and dissociation processes were estimated using two complementary methods: umbrella sampling simulations of two DNA fragments and direct observation of the aggregation process in a system containing 458 DNA fragments. We found the end-to-end force to be short range, attractive, hydrophobic and only weakly dependent on the ion concentration. The relation between the stacking free energy and end-to-end attraction is discussed as well as possible roles of the end-to-end interaction in biological and nanotechnological systems.

post to: CiteULike

The structural basis for the oligomerization of the N-terminal domain of SATB1.

Publication Date: 2012 Jan 12 PMID: 22241778
Authors: Wang, Z. - Yang, X. - Chu, X. - Zhang, J. - Zhou, H. - Shen, Y. - Long, J.
Journal: Nucleic Acids Res

Special AT-rich sequence-binding protein 1 (SATB1) is a global chromatin organizer and gene expression regulator essential for T-cell development and breast cancer tumor growth and metastasis. The oligomerization of the N-terminal domain of SATB1 is critical for its biological function. We determined the crystal structure of the N-terminal domain of SATB1. Surprisingly, this domain resembles a ubiquitin domain instead of the previously proposed PDZ domain. Our results also reveal that SATB1 can form a tetramer through its N-terminal domain. The tetramerization of SATB1 plays an essential role in its binding to highly specialized DNA sequences. Furthermore, isothermal titration calorimetry results indicate that the SATB1 tetramer can bind simultaneously to two DNA targets. Based on these results, we propose a molecular model whereby SATB1 regulates the expression of multiple genes both locally and at a distance.

post to: CiteULike

Stoichiometry of MutS and MutL at unrepaired mismatches in vivo suggests a mechanism of repair.

Publication Date: 2012 Jan 12 PMID: 22241777
Authors: Elez, M. - Radman, M. - Matic, I.
Journal: Nucleic Acids Res

Mismatch repair (MMR) is an evolutionarily conserved DNA repair system, which corrects mismatched bases arising during DNA replication. MutS recognizes and binds base pair mismatches, while the MutL protein interacts with MutS-mismatch complex and triggers MutH endonuclease activity at a distal-strand discrimination site on the DNA. The mechanism of communication between these two distal sites on the DNA is not known. We used functional fluorescent MMR proteins, MutS and MutL, in order to investigate the formation of the fluorescent MMR protein complexes on mismatches in real-time in growing Escherichia coli cells. We found that MutS and MutL proteins co-localize on unrepaired mismatches to form fluorescent foci. MutL foci were, on average, 2.7 times more intense than the MutS foci co-localized on individual mismatches. A steric block on the DNA provided by the MutHE56A mutant protein, which binds to but does not cut the DNA at the strand discrimination site, decreased MutL foci fluorescence 3-fold. This indicates that MutL accumulates from the mismatch site toward strand discrimination site along the DNA. Our results corroborate the hypothesis postulating that MutL accumulation assures the coordination of the MMR activities between the mismatch and the strand discrimination site.

post to: CiteULike

Comprehensive literature review and statistical considerations for GWAS meta-analysis.

Publication Date: 2012 Jan 12 PMID: 22241776
Authors: Begum, F. - Ghosh, D. - Tseng, G. C. - Feingold, E.
Journal: Nucleic Acids Res

Over the last decade, genome-wide association studies (GWAS) have become the standard tool for gene discovery in human disease research. While debate continues about how to get the most out of these studies and on occasion about how much value these studies really provide, it is clear that many of the strongest results have come from large-scale mega-consortia and/or meta-analyses that combine data from up to dozens of studies and tens of thousands of subjects. While such analyses are becoming more and more common, statistical methods have lagged somewhat behind. There are good meta-analysis methods available, but even when they are carefully and optimally applied there remain some unresolved statistical issues. This article systematically reviews the GWAS meta-analysis literature, highlighting methodology and software options and reviewing methods that have been used in real studies. We illustrate differences among methods using a case study. We also discuss some of the unresolved issues and potential future directions.

post to: CiteULike

Structural bias in T4 RNA ligase-mediated 3'-adapter ligation.

Publication Date: 2012 Jan 12 PMID: 22241775
Authors: Zhuang, F. - Fuchs, R. T. - Sun, Z. - Zheng, Y. - Robb, G. B.
Journal: Nucleic Acids Res

T4 RNA ligases are commonly used to attach adapters to RNAs, but large differences in ligation efficiency make detection and quantitation problematic. We developed a ligation selection strategy using random RNAs in combination with high-throughput sequencing to gain insight into the differences in efficiency of ligating pre-adenylated DNA adapters to RNA 3'-ends. After analyzing biases in RNA sequence, secondary structure and RNA-adapter cofold structure, we conclude that T4 RNA ligases do not show significant primary sequence preference in RNA substrates, but are biased against structural features within RNAs and adapters. Specifically, RNAs with less than three unstructured nucleotides at the 3'-end and RNAs that are predicted to cofold with an adapter in unfavorable structures are likely to be poorly ligated. The effect of RNA-adapter cofold structures on ligation is supported by experiments where the ligation efficiency of specific miRNAs was changed by designing adapters to alter cofold structure. In addition, we show that using adapters with randomized regions results in higher ligation efficiency and reduced ligation bias. We propose that using randomized adapters may improve RNA representation in experiments that include a 3'-adapter ligation step.

post to: CiteULike

Kinetics and mechanism of G-quadruplex formation and conformational switch in a G-quadruplex of PS2.M induced by Pb2+

Publication Date: 2012 Jan 12 PMID: 22241774
Authors: Liu, W. - Zhu, H. - Zheng, B. - Cheng, S. - Fu, Y. - Li, W. - Lau, T. C. - Liang, H.
Journal: Nucleic Acids Res

DNA sequences with guanine repeats can form G-quartets that adopt G-quadruplex structures in the presence of specific metal ions. Using circular dichroism (CD) and ultraviolet-visible (UV-Vis) spectroscopy, we determined the spectral characteristics and the overall conformation of a G-quadruplex of PS2.M with an oligonucleotide sequence, d(GTG(3)TAG(3)CG(3)TTG(2)). UV-melting curves demonstrate that the Pb(2+)-induced G-quadruplex formed unimolecularly and the highest melting temperature (T(m)) is 72 degrees C. The analysis of the UV titration results reveals that the binding stoichiometry of Pb(2+) ions to PS2.M is two, suggesting that the Pb(2+) ions coordinate between adjacent G-quartets. Binding of ions to G-rich DNA is a complex multiple-pathway process, which is strongly affected by the type of the cations. Kinetic studies suggest that the Pb(2+)-induced folding of PS2.M to G-quadruplex probably proceeds through a three-step pathway involving two intermediates. Structural transition occurs after adding Pb(NO(3))(2) to the Na(+)- or K(+)-induced G-quadruplexes, which may be attributed to the replacement of Na(+) or K(+) by Pb(2+) ions and the generation of a more compact Pb(2+)-PS2.M structure. Comparison of the relaxation times shows that the Na(+)-->Pb(2+) exchange is more facile than the K(+)-->Pb(2+) exchange process, and the mechanisms for these processes are proposed.

post to: CiteULike

Conformational and thermodynamic properties modulate the nucleotide excision repair of 2-aminofluorene and 2-acetylaminofluorene dG adducts in the NarI sequence.

Publication Date: 2012 Jan 12 PMID: 22241773
Authors: Jain, V. - Hilton, B. - Patnaik, S. - Zou, Y. - Chiarelli, M. P. - Cho, B. P.
Journal: Nucleic Acids Res

Nucleotide excision repair (NER) is a major repair pathway that recognizes and corrects various lesions in cellular DNA. We hypothesize that damage recognition is an initial step in NER that senses conformational anomalies in the DNA caused by lesions. We prepared three DNA duplexes containing the carcinogen adduct N-(2'-deoxyguanosin-8-yl)-7-fluoro-2-acetylaminofluorene (FAAF) at G(1), G(2) or G(3) of NarI sequence (5'-CCG(1)G(2)CG(3)CC-3'). Our (19)F-NMR/ICD results showed that FAAF at G(1) and G(3) prefer syn S- and W-conformers, whereas anti B-conformer was predominant for G(2). We found that the repair of FAAF occurs in a conformation-specific manner, i.e. the highly S/W-conformeric G(3) and -G(1) duplexes incised more efficiently than the B-type G(2) duplex (G(3) approximately G(1 )> G(2)). The melting and thermodynamic data indicate that the S- and W-conformers produce greater DNA distortion and thermodynamic destabilization. The N-deacetylated N-(2'-deoxyguanosin-8-yl)-7-fluoro-2-aminofluorene (FAF) adducts in the same NarI sequence are repaired 2- to 3-fold less than FAAF: however, the incision efficiency was in order of G(2) approximately G(1 )> G(3), a reverse trend of the FAAF case. We have envisioned the so-called N-acetyl factor as it could raise conformational barriers of FAAF versus FAF. The present results provide valuable conformational insight into the sequence-dependent UvrABC incisions of the bulky aminofluorene DNA adducts.

post to: CiteULike

SLiCE: a novel bacterial cell extract-based DNA cloning method.

Publication Date: 2012 Jan 12 PMID: 22241772
Authors: Zhang, Y. - Werling, U. - Edelmann, W.
Journal: Nucleic Acids Res

We describe a novel cloning method termed SLiCE (Seamless Ligation Cloning Extract) that utilizes easy to generate bacterial cell extracts to assemble multiple DNA fragments into recombinant DNA molecules in a single in vitro recombination reaction. SLiCE overcomes the sequence limitations of traditional cloning methods, facilitates seamless cloning by recombining short end homologies (>/=15 bp) with or without flanking heterologous sequences and provides an effective strategy for directional subcloning of DNA fragments from Bacteria Artificial Chromosomes (BACs) or other sources. SLiCE is highly cost effective as a number of standard laboratory bacterial strains can serve as sources for SLiCE extract. In addition, the cloning efficiencies and capabilities of these strains can be greatly improved by simple genetic modifications. As an example, we modified the DH10B Escherichia coli strain to express an optimized lambda prophage Red recombination system. This strain, termed PPY, facilitates SLiCE with very high efficiencies and demonstrates the versatility of the method.

post to: CiteULike

New methods for finding common insertion sites and co-occurring common insertion sites in transposon- and virus-based genetic screens.

Publication Date: 2012 Jan 12 PMID: 22241771
Authors: Bergemann, T. L. - Starr, T. K. - Yu, H. - Steinbach, M. - Erdmann, J. - Chen, Y. - Cormier, R. T. - Largaespada, D. A. - Silverstein, K. A.
Journal: Nucleic Acids Res

Insertional mutagenesis screens in mice are used to identify individual genes that drive tumor formation. In these screens, candidate cancer genes are identified if their genomic location is proximal to a common insertion site (CIS) defined by high rates of transposon or retroviral insertions in a given genomic window. In this article, we describe a new method for defining CISs based on a Poisson distribution, the Poisson Regression Insertion Model, and show that this new method is an improvement over previously described methods. We also describe a modification of the method that can identify pairs and higher orders of co-occurring common insertion sites. We apply these methods to two data sets, one generated in a transposon-based screen for gastrointestinal tract cancer genes and another based on the set of retroviral insertions in the Retroviral Tagged Cancer Gene Database. We show that the new methods identify more relevant candidate genes and candidate gene pairs than found using previous methods. Identification of the biologically relevant set of mutations that occur in a single cell and cause tumor progression will aid in the rational design of single and combinatorial therapies in the upcoming age of personalized cancer therapy.

post to: CiteULike

Structural and biochemical characterization of HP0315 from Helicobacter pylori as a VapD protein with an endoribonuclease activity.

Publication Date: 2012 Jan 12 PMID: 22241770
Authors: Kwon, A. R. - Kim, J. H. - Park, S. J. - Lee, K. Y. - Min, Y. H. - Im, H. - Lee, I. - Lee, K. Y. - Lee, B. J.
Journal: Nucleic Acids Res

VapD-like virulence-associated proteins have been found in many organisms, but little is known about this protein family including the 3D structure of these proteins. Recently, a relationship between the Cas2 family of ribonucleases associated with the CRISPR system of microbial immunity and VapD was suggested. Here, we show for the first time the structure of a member of the VapD family and present a relationship of VapD with Cas2 family and toxin-antitoxin (TA) systems. The crystal structure of HP0315 from Helicobacter pylori was solved at a resolution of 2.8 A. The structure of HP0315, which has a modified ferredoxin-like fold, is very similar to that of the Cas2 family. Like Cas2 proteins, HP0315 shows endoribonuclease activity. HP0315-cleaved mRNA, mainly before A and G nucleotides preferentially, which means that HP0315 has purine-specific endoribonuclease activity. Mutagenesis studies of HP0315 revealed that D7, L13, S43 and D76 residues are important for RNase activity, in contrast, to the Cas2 family. HP0315 is arranged as an operon with HP0316, which was found to be an antitoxin-related protein. However, HP0315 is not a component of the TA system. Thus, HP0315 may be an evolutionary intermediate which does not belong to either the Cas2 family or TA system.

post to: CiteULike

Distal chromatin structure influences local nucleosome positions and gene expression.

Publication Date: 2012 Jan 12 PMID: 22241769
Authors: Jansen, A. - van der Zande, E. - Meert, W. - Fink, G. R. - Verstrepen, K. J.
Journal: Nucleic Acids Res

The positions of nucleosomes across the genome influence several cellular processes, including gene transcription. However, our understanding of the factors dictating where nucleosomes are located and how this affects gene regulation is still limited. Here, we perform an extensive in vivo study to investigate the influence of the neighboring chromatin structure on local nucleosome positioning and gene expression. Using truncated versions of the Saccharomyces cerevisiae URA3 gene, we show that nucleosome positions in the URA3 promoter are at least partly determined by the local DNA sequence, with so-called 'antinucleosomal elements' like poly(dA:dT) tracts being key determinants of nucleosome positions. In addition, we show that changes in the nucleosome positions in the URA3 promoter strongly affect the promoter activity. Most interestingly, in addition to demonstrating the effect of the local DNA sequence, our study provides novel in vivo evidence that nucleosome positions are also affected by the position of neighboring nucleosomes. Nucleosome structure may therefore be an important selective force for conservation of gene order on a chromosome, because relocating a gene to another genomic position (where the positions of neighboring nucleosomes are different from the original locus) can have dramatic consequences for the gene's nucleosome structure and thus its expression.

post to: CiteULike

Sensitive and label-free biosensing of RNA with predicted secondary structures by a triplex affinity capture method.

Publication Date: 2012 Jan 12 PMID: 22241768
Authors: Carrascosa, L. G. - Gomez-Montes, S. - Avino, A. - Nadal, A. - Pla, M. - Eritja, R. - Lechuga, L. M.
Journal: Nucleic Acids Res

A novel biosensing approach for the label-free detection of nucleic acid sequences of short and large lengths has been implemented, with special emphasis on targeting RNA sequences with secondary structures. The approach is based on selecting 8-aminoadenine-modified parallel-stranded DNA tail-clamps as affinity bioreceptors. These receptors have the ability of creating a stable triplex-stranded helix at neutral pH upon hybridization with the nucleic acid target. A surface plasmon resonance biosensor has been used for the detection. With this strategy, we have detected short DNA sequences (32-mer) and purified RNA (103-mer) at the femtomol level in a few minutes in an easy and level-free way. This approach is particularly suitable for the detection of RNA molecules with predicted secondary structures, reaching a limit of detection of 50 fmol without any label or amplification steps. Our methodology has shown a marked enhancement for the detection (18% for short DNA and 54% for RNA), when compared with the conventional duplex approach, highlighting the large difficulty of the duplex approach to detect nucleic acid sequences, especially those exhibiting stable secondary structures. We believe that our strategy could be of great interest to the RNA field.

post to: CiteULike

G-quadruplex structure and stability illuminated by 2-aminopurine phasor plots.

Publication Date: 2012 Jan 12 PMID: 22241767
Authors: Buscaglia, R. - Jameson, D. M. - Chaires, J. B.
Journal: Nucleic Acids Res

The use of time-resolved fluorescence measurements in studies of telomeric G-quadruplex folding and stability has been hampered by the complexity of fluorescence lifetime distributions in solution. The application of phasor diagrams to the analysis of time-resolved fluorescence measurements, collected from either frequency-domain or time-domain instrumentation, allows for rapid characterization of complex lifetime distributions. Phasor diagrams are model-free graphical representations of transformed time-resolved fluorescence results. Simplification of complex fluorescent decays by phasor diagrams is demonstrated here using a 2-aminopurine substituted telomeric G-quadruplex sequence. The application of phasor diagrams to complex systems is discussed with comparisons to traditional non-linear regression model fitting. Phasor diagrams allow for the folding and stability of the telomeric G-quadruplex to be monitored in the presence of either sodium or potassium. Fluorescence lifetime measurements revealed multiple transitions upon folding of the telomeric G-quadruplex through the addition of potassium. Enzymatic digestion of the telomeric G-quadruplex structure, fluorescence quenching and Forster resonance energy transfer were also monitored through phasor diagrams. This work demonstrates the sensitivity of time-resolved methods for monitoring changes to the telomeric G-quadruplex and outlines the phasor diagram approach for analysis of complex time-resolved results that can be extended to other G-quadruplex and nucleic acid systems.

post to: CiteULike

Quadruplex-single nucleotide polymorphisms (Quad-SNP) influence gene expression difference among individuals.

Publication Date: 2012 Jan 11 PMID: 22238381
Authors: Baral, A. - Kumar, P. - Halder, R. - Mani, P. - Yadav, V. K. - Singh, A. - Das, S. K. - Chowdhury, S.
Journal: Nucleic Acids Res

Non-canonical guanine quadruplex structures are not only predominant but also conserved among bacterial and mammalian promoters. Moreover recent findings directly implicate quadruplex structures in transcription. These argue for an intrinsic role of the structural motif and thereby posit that single nucleotide polymorphisms (SNP) that compromise the quadruplex architecture could influence function. To test this, we analysed SNPs within quadruplex motifs (Quad-SNP) and gene expression in 270 individuals across four populations (HapMap) representing more than 14 500 genotypes. Findings reveal significant association between quadruplex-SNPs and expression of the corresponding gene in individuals (P < 0.0001). Furthermore, analysis of Quad-SNPs obtained from population-scale sequencing of 1000 human genomes showed relative selection bias against alteration of the structural motif. To directly test the quadruplex-SNP-transcription connection, we constructed a reporter system using the RPS3 promoter-remarkable difference in promoter activity in the 'quadruplex-destabilized' versus 'quadruplex-intact' promoter was noticed. As a further test, we incorporated a quadruplex motif or its disrupted counterpart within a synthetic promoter reporter construct. The quadruplex motif, and not the disrupted-motif, enhanced transcription in human cell lines of different origin. Together, these findings build direct support for quadruplex-mediated transcription and suggest quadruplex-SNPs may play significant role in mechanistically understanding variations in gene expression among individuals.

post to: CiteULike

The RNA helicase RHAU (DHX36) unwinds a G4-quadruplex in human telomerase RNA and promotes the formation of the P1 helix template boundary.

Publication Date: 2012 Jan 11 PMID: 22238380
Authors: Booy, E. P. - Meier, M. - Okun, N. - Novakowski, S. K. - Xiong, S. - Stetefeld, J. - McKenna, S. A.
Journal: Nucleic Acids Res

Human telomerase RNA (hTR) contains several guanine tracts at its 5'-end that can form a G4-quadruplex structure. Previous evidence suggests that a G4-quadruplex within this region disrupts the formation of an important structure within hTR known as the P1 helix, a critical element in defining the template boundary for reverse transcription. RNA associated with AU-rich element (RHAU) is an RNA helicase that has specificity for DNA and RNA G4-quadruplexes. Two recent studies identify a specific interaction between hTR and RHAU. Herein, we confirm this interaction and identify the minimally interacting RNA fragments. We demonstrate the existence of multiple quadruplex structures within the 5' region of hTR and find that these regions parallel the minimal sequences capable of RHAU interaction. We confirm the importance of the RHAU-specific motif in the interaction with hTR and demonstrate that the helicase activity of RHAU is sufficient to unwind the quadruplex and promote an interaction with 25 internal nucleotides to form a stable P1 helix. Furthermore, we have found that a 5'-terminal quadruplex persists following P1 helix formation that retains affinity for RHAU. Finally, we have investigated the functional implications of this interaction and demonstrated a reduction in average telomere length following RHAU knockdown by small interfering RNA (siRNA).

post to: CiteULike

Young intragenic miRNAs are less coexpressed with host genes than old ones: implications of miRNA-host gene coevolution.

Publication Date: 2012 Jan 11 PMID: 22238379
Authors: He, C. - Li, Z. - Chen, P. - Huang, H. - Hurst, L. D. - Chen, J.
Journal: Nucleic Acids Res

MicroRNAs (miRNAs) have emerged as key regulators of gene expression. Intragenic miRNAs account for approximately 50% of mammalian miRNAs. Classic studies reported that they are usually coexpressed with host genes. Here, using genome-wide miRNA and gene expression profiles from five sample sets, we show that evolutionarily conserved ('old') intragenic miRNAs tend to be coexpressed with host genes, but non-conserved ('young') ones rarely do so. This result is robust: in all sample sets, the coexpression rate of young miRNAs is significantly lower than that of conserved ones even after controlling for abundance. As a result, although young miRNAs dominate in human genome, the majority of intragenic miRNAs that show coexpression with host genes are phylogenetically old ones. For younger miRNAs, extrapolation of their expression profiles from those of their host genes should be treated with caution. We propose a model to explain this phenomenon in which the majority of young miRNAs are unlikely to be coexpressed with host genes; however, for some fraction of young miRNAs coexpression with their host genes, initially imbued by chromatin level effects, is advantageous and these are the ones likely to embed into the system and evolve ever higher levels of coexpression, possibly by evolving piggybacking mechanisms.

post to: CiteULike

A novel phage-encoded transcription antiterminator acts by suppressing bacterial RNA polymerase pausing.

Publication Date: 2012 Jan 11 PMID: 22238378
Authors: Berdygulova, Z. - Esyunina, D. - Miropolskaya, N. - Mukhamedyarov, D. - Kuznedelov, K. - Nickels, B. E. - Severinov, K. - Kulbachinskiy, A. - Minakhin, L.
Journal: Nucleic Acids Res

Gp39, a small protein encoded by Thermus thermophilus phage P23-45, specifically binds the host RNA polymerase (RNAP) and inhibits transcription initiation. Here, we demonstrate that gp39 also acts as an antiterminator during transcription through intrinsic terminators. The antitermination activity of gp39 relies on its ability to suppress transcription pausing at poly(U) tracks. Gp39 also accelerates transcription elongation by decreasing RNAP pausing and backtracking but does not significantly affect the rates of catalysis of individual reactions in the RNAP active center. We mapped the RNAP-gp39 interaction site to the beta flap, a domain that forms a part of the RNA exit channel and is also a likely target for lambda phage antiterminator proteins Q and N, and for bacterial elongation factor NusA. However, in contrast to Q and N, gp39 does not depend on NusA or other auxiliary factors for its activity. To our knowledge, gp39 is the first characterized phage-encoded transcription factor that affects every step of the transcription cycle and suppresses transcription termination through its antipausing activity.

post to: CiteULike

Structural insights into the redox-switch mechanism of the MarR/DUF24-type regulator HypR.

Publication Date: 2012 Jan 11 PMID: 22238377
Authors: Palm, G. J. - Khanh Chi, B. - Waack, P. - Gronau, K. - Becher, D. - Albrecht, D. - Hinrichs, W. - Read, R. J. - Antelmann, H.
Journal: Nucleic Acids Res

Bacillus subtilis encodes redox-sensing MarR-type regulators of the OhrR and DUF24-families that sense organic hydroperoxides, diamide, quinones or aldehydes via thiol-based redox-switches. In this article, we characterize the novel redox-sensing MarR/DUF24-family regulator HypR (YybR) that is activated by disulphide stress caused by diamide and NaOCl in B. subtilis. HypR controls positively a flavin oxidoreductase HypO that confers protection against NaOCl stress. The conserved N-terminal Cys14 residue of HypR has a lower pK(a) of 6.36 and is essential for activation of hypO transcription by disulphide stress. HypR resembles a 2-Cys-type regulator that is activated by Cys14-Cys49' intersubunit disulphide formation. The crystal structures of reduced and oxidized HypR proteins were resolved revealing structural changes of HypR upon oxidation. In reduced HypR a hydrogen-bonding network stabilizes the reactive Cys14 thiolate that is 8-9 A apart from Cys49'. HypR oxidation breaks these H-bonds, reorients the monomers and moves the major groove recognition alpha4 and alpha4' helices approximately 4 A towards each other. This is the first crystal structure of a redox-sensing MarR/DUF24 family protein in bacteria that is activated by NaOCl stress. Since hypochloric acid is released by activated macrophages, related HypR-like regulators could function to protect pathogens against the host immune defense.

post to: CiteULike

C7orf30 is necessary for biogenesis of the large subunit of the mitochondrial ribosome.

Publication Date: 2012 Jan 11 PMID: 22238376
Authors: Rorbach, J. - Gammage, P. A. - Minczuk, M.
Journal: Nucleic Acids Res

Defects of the translation apparatus in human mitochondria are known to cause disease, yet details of how protein synthesis is regulated in this organelle remain to be unveiled. Here, we characterize a novel human protein, C7orf30 that contributes critically to mitochondrial translation and specifically associates with the large subunit of the mitochondrial ribosome (mt-LSU). Inactivation of C7orf30 in human cells by RNA interference results in respiratory incompetence owing to reduced mitochondrial translation rates without any appreciable effects on the steady-state levels of mitochondrial mRNAs and rRNAs. Ineffective translation in C7orf30-depleted cells or cells overexpressing a dominant-negative mutant of the protein results from aberrant assembly of mt-LSU and consequently reduced formation of the monosome. These findings lead us to propose that C7orf30 is a human assembly and/or stability factor involved in the biogenesis of the large subunit of the mitochondrial ribosome.

post to: CiteULike

C7orf30 specifically associates with the large subunit of the mitochondrial ribosome and is involved in translation.

Publication Date: 2012 Jan 11 PMID: 22238375
Authors: Wanschers, B. F. - Szklarczyk, R. - Pajak, A. - van den Brand, M. A. - Gloerich, J. - Rodenburg, R. J. - Lightowlers, R. N. - Nijtmans, L. G. - Huynen, M. A.
Journal: Nucleic Acids Res

In a comparative genomics study for mitochondrial ribosome-associated proteins, we identified C7orf30, the human homolog of the plant protein iojap. Gene order conservation among bacteria and the observation that iojap orthologs cannot be transferred between bacterial species predict this protein to be associated with the mitochondrial ribosome. Here, we show colocalization of C7orf30 with the large subunit of the mitochondrial ribosome using isokinetic sucrose gradient and 2D Blue Native polyacrylamide gel electrophoresis (BN-PAGE) analysis. We co-purified C7orf30 with proteins of the large subunit, and not with proteins of the small subunit, supporting interaction that is specific to the large mitoribosomal complex. Consistent with this physical association, a mitochondrial translation assay reveals negative effects of C7orf30 siRNA knock-down on mitochondrial gene expression. Based on our data we propose that C7orf30 is involved in ribosomal large subunit function. Sequencing the gene in 35 patients with impaired mitochondrial translation did not reveal disease-causing mutations in C7orf30.

post to: CiteULike

BC1-FMRP interaction is modulated by 2'-O-methylation: RNA-binding activity of the tudor domain and translational regulation at synapses.

Publication Date: 2012 Jan 11 PMID: 22238374
Authors: Lacoux, C. - Di Marino, D. - Pilo Boyl, P. - Zalfa, F. - Yan, B. - Ciotti, M. T. - Falconi, M. - Urlaub, H. - Achsel, T. - Mougin, A. - Caizergues-Ferrer, M. - Bagni, C.
Journal: Nucleic Acids Res

The brain cytoplasmic RNA, BC1, is a small non-coding RNA that is found in different RNP particles, some of which are involved in translational control. One component of BC1-containing RNP complexes is the fragile X mental retardation protein (FMRP) that is implicated in translational repression. Peptide mapping and computational simulations show that the tudor domain of FMRP makes specific contacts to BC1 RNA. Endogenous BC1 RNA is 2'-O-methylated in nucleotides that contact the FMRP interface, and methylation can affect this interaction. In the cell body BC1 2'-O-methylations are present in both the nucleus and the cytoplasm, but they are virtually absent at synapses where the FMRP-BC1-mRNA complex exerts its function. These results strongly suggest that subcellular region-specific modifications of BC1 affect the binding to FMRP and the interaction with its mRNA targets. We finally show that BC1 RNA has an important role in translation of certain mRNAs associated to FMRP. All together these findings provide further insights into the translational regulation by the FMRP-BC1 complex at synapses.

post to: CiteULike

The molecular basis of ATM-dependent dimerization of the Mdc1 DNA damage checkpoint mediator.

Publication Date: 2012 Jan 10 PMID: 22234878
Authors: Jungmichel, S. - Clapperton, J. A. - Lloyd, J. - Hari, F. J. - Spycher, C. - Pavic, L. - Li, J. - Haire, L. F. - Bonalli, M. - Larsen, D. H. - Lukas, C. - Lukas, J. - Macmillan, D. - Nielsen, M. L. - Stucki, M. - Smerdon, S. J.
Journal: Nucleic Acids Res

Mdc1 is a large modular phosphoprotein scaffold that maintains signaling and repair complexes at double-stranded DNA break sites. Mdc1 is anchored to damaged chromatin through interaction of its C-terminal BRCT-repeat domain with the tail of gammaH2AX following DNA damage, but the role of the N-terminal forkhead-associated (FHA) domain remains unclear. We show that a major binding target of the Mdc1 FHA domain is a previously unidentified DNA damage and ATM-dependent phosphorylation site near the N-terminus of Mdc1 itself. Binding to this motif stabilizes a weak self-association of the FHA domain to form a tight dimer. X-ray structures of free and complexed Mdc1 FHA domain reveal a 'head-to-tail' dimerization mechanism that is closely related to that seen in pre-activated forms of the Chk2 DNA damage kinase, and which both positively and negatively influences Mdc1 FHA domain-mediated interactions in human cells prior to and following DNA damage.

post to: CiteULike

Structural mechanism of the phosphorylation-dependent dimerization of the MDC1 forkhead-associated domain.

Publication Date: 2012 Jan 10 PMID: 22234877
Authors: Liu, J. - Luo, S. - Zhao, H. - Liao, J. - Li, J. - Yang, C. - Xu, B. - Stern, D. F. - Xu, X. - Ye, K.
Journal: Nucleic Acids Res

MDC1 is a key mediator of the DNA-damage response in mammals with several phosphorylation-dependent protein interaction domains. The function of its N-terminal forkhead-associated (FHA) domain remains elusive. Here, we show with structural, biochemical and cellular data that the FHA domain mediates phosphorylation-dependent dimerization of MDC1 in response to DNA damage. Crystal structures of the FHA domain reveal a face-to-face dimer with pseudo-dyad symmetry. We found that the FHA domain recognizes phosphothreonine 4 (pT4) at the N-terminus of MDC1 and determined its crystal structure in complex with a pT4 peptide. Biochemical analysis further revealed that in the dimer, the FHA domain binds in trans to pT4 from the other subunit, which greatly stabilizes the otherwise unstable dimer. We show that T4 is phosphorylated primarily by ATM upon DNA damage. MDC1 mutants with the FHA domain deleted or impaired in its ability to dimerize formed fewer foci at DNA-damage sites, but the localization defect was largely rescued by an artificial dimerization module, suggesting that dimerization is the primary function of the MDC1 FHA domain. Our results suggest a novel mechanism for the regulation of MDC1 function through T4 phosphorylation and FHA-mediated dimerization.

post to: CiteULike

Use of ChIP-Seq data for the design of a multiple promoter-alignment method.

Publication Date: 2012 Jan 9 PMID: 22230796
Authors: Erb, I. - Gonzalez-Vallinas, J. R. - Bussotti, G. - Blanco, E. - Eyras, E. - Notredame, C.
Journal: Nucleic Acids Res

We address the challenge of regulatory sequence alignment with a new method, Pro-Coffee, a multiple aligner specifically designed for homologous promoter regions. Pro-Coffee uses a dinucleotide substitution matrix estimated on alignments of functional binding sites from TRANSFAC. We designed a validation framework using several thousand families of orthologous promoters. This dataset was used to evaluate the accuracy for predicting true human orthologs among their paralogs. We found that whereas other methods achieve on average 73.5% accuracy, and 77.6% when trained on that same dataset, the figure goes up to 80.4% for Pro-Coffee. We then applied a novel validation procedure based on multi-species ChIP-seq data. Trained and untrained methods were tested for their capacity to correctly align experimentally detected binding sites. Whereas the average number of correctly aligned sites for two transcription factors is 284 for default methods and 316 for trained methods, Pro-Coffee achieves 331, 16.5% above the default average. We find a high correlation between a method's performance when classifying orthologs and its ability to correctly align proven binding sites. Not only has this interesting biological consequences, it also allows us to conclude that any method that is trained on the ortholog data set will result in functionally more informative alignments.

post to: CiteULike

Correction of RT-qPCR data for genomic DNA-derived signals with ValidPrime.

Publication Date: 2012 Jan 6 PMID: 22228834
Authors: Laurell, H. - Iacovoni, J. S. - Abot, A. - Svec, D. - Maoret, J. J. - Arnal, J. F. - Kubista, M.
Journal: Nucleic Acids Res

Genomic DNA (gDNA) contamination is an inherent problem during RNA purification that can lead to non-specific amplification and aberrant results in reverse transcription quantitative PCR (RT-qPCR). Currently, there is no alternative to RT(-) controls to evaluate the impact of the gDNA background on RT-PCR data. We propose a novel method (ValidPrime) that is more accurate than traditional RT(-) controls to test qPCR assays with respect to their sensitivity toward gDNA. ValidPrime measures the gDNA contribution using an optimized gDNA-specific ValidPrime assay (VPA) and gDNA reference sample(s). The VPA, targeting a non-transcribed locus, is used to measure the gDNA contents in RT(+) samples and the gDNA reference is used to normalize for GOI-specific differences in gDNA sensitivity. We demonstrate that the RNA-derived component of the signal can be accurately estimated and deduced from the total signal. ValidPrime corrects with high precision for both exogenous (spiked) and endogenous gDNA, contributing approximately 60% of the total signal, whereas substantially reducing the number of required qPCR control reactions. In conclusion, ValidPrime offers a cost-efficient alternative to RT(-) controls and accurately corrects for signals derived from gDNA in RT-qPCR.

post to: CiteULike

A highly efficient and effective motif discovery method for ChIP-seq/ChIP-chip data using positional information.

Publication Date: 2012 Jan 6 PMID: 22228832
Authors: Ma, X. - Kulkarni, A. - Zhang, Z. - Xuan, Z. - Serfling, R. - Zhang, M. Q.
Journal: Nucleic Acids Res

Identification of DNA motifs from ChIP-seq/ChIP-chip [chromatin immunoprecipitation (ChIP)] data is a powerful method for understanding the transcriptional regulatory network. However, most established methods are designed for small sample sizes and are inefficient for ChIP data. Here we propose a new k-mer occurrence model to reflect the fact that functional DNA k-mers often cluster around ChIP peak summits. With this model, we introduced a new measure to discover functional k-mers. Using simulation, we demonstrated that our method is more robust against noises in ChIP data than available methods. A novel word clustering method is also implemented to group similar k-mers into position weight matrices (PWMs). Our method was applied to a diverse set of ChIP experiments to demonstrate its high sensitivity and specificity. Importantly, our method is much faster than several other methods for large sample sizes. Thus, we have developed an efficient and effective motif discovery method for ChIP experiments.

post to: CiteULike

XRCC4's interaction with XLF is required for coding (but not signal) end joining.

Publication Date: 2012 Jan 6 PMID: 22228831
Authors: Roy, S. - Andres, S. N. - Vergnes, A. - Neal, J. A. - Xu, Y. - Yu, Y. - Lees-Miller, S. P. - Junop, M. - Modesti, M. - Meek, K.
Journal: Nucleic Acids Res

XRCC4 and XLF are structurally related proteins important for DNA Ligase IV function. XRCC4 forms a tight complex with DNA Ligase IV while XLF interacts directly with XRCC4. Both XRCC4 and XLF form homodimers that can polymerize as heterotypic filaments independently of DNA Ligase IV. Emerging structural and in vitro biochemical data suggest that XRCC4 and XLF together generate a filamentous structure that promotes bridging between DNA molecules. Here, we show that ablating XRCC4's affinity for XLF results in DNA repair deficits including a surprising deficit in VDJ coding, but not signal end joining. These data are consistent with a model whereby XRCC4/XLF complexes hold DNA ends together-stringently required for coding end joining, but dispensable for signal end joining. Finally, DNA-PK phosphorylation of XRCC4/XLF complexes disrupt DNA bridging in vitro, suggesting a regulatory role for DNA-PK's phosphorylation of XRCC4/XLF complexes.

post to: CiteULike

Genome-wide occupancy links Hoxa2 to Wnt-beta-catenin signaling in mouse embryonic development.

Publication Date: 2012 Jan 5 PMID: 22223247
Authors: Donaldson, I. J. - Amin, S. - Hensman, J. J. - Kutejova, E. - Rattray, M. - Lawrence, N. - Hayes, A. - Ward, C. M. - Bobola, N.
Journal: Nucleic Acids Res

The regulation of gene expression is central to developmental programs and largely depends on the binding of sequence-specific transcription factors with cis-regulatory elements in the genome. Hox transcription factors specify the spatial coordinates of the body axis in all animals with bilateral symmetry, but a detailed knowledge of their molecular function in instructing cell fates is lacking. Here, we used chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq) to identify Hoxa2 genomic locations in a time and space when it is actively instructing embryonic development in mouse. Our data reveals that Hoxa2 has large genome coverage and potentially regulates thousands of genes. Sequence analysis of Hoxa2-bound regions identifies high occurrence of two main classes of motifs, corresponding to Hox and Pbx-Hox recognition sequences. Examination of the binding targets of Hoxa2 faithfully captures the processes regulated by Hoxa2 during embryonic development; in addition, it uncovers a large cluster of potential targets involved in the Wnt-signaling pathway. In vivo examination of canonical Wnt-beta-catenin signaling reveals activity specifically in Hoxa2 domain of expression, and this is undetectable in Hoxa2 mutant embryos. The comprehensive mapping of Hoxa2-binding sites provides a framework to study Hox regulatory networks in vertebrate developmental processes.

post to: CiteULike

Visualization of a DNA-PK/PARP1 complex.

Publication Date: 2012 Jan 5 PMID: 22223246
Authors: Spagnolo, L. - Barbeau, J. - Curtin, N. J. - Morris, E. P. - Pearl, L. H.
Journal: Nucleic Acids Res

The DNA-dependent protein kinase (DNA-PK) and Poly(ADP-ribose) polymerase-1 (PARP1) are critical enzymes that reduce genomic damage caused by DNA lesions. They are both activated by DNA strand breaks generated by physiological and environmental factors, and they have been shown to interact. Here, we report in vivo evidence that DNA-PK and PARP1 are equally necessary for rapid repair. We purified a DNA-PK/PARP1 complex loaded on DNA and performed electron microscopy and single particle analysis on its tetrameric and dimer-of-tetramers forms. By comparison with the DNA-PK holoenzyme and fitting crystallographic structures, we see that the PARP1 density is in close contact with the Ku subunit. Crucially, PARP1 binding elicits substantial conformational changes in the DNA-PK synaptic dimer assembly. Taken together, our data support a functional, in-pathway role for DNA-PK and PARP1 in double-strand break (DSB) repair. We also propose a NHEJ model where protein-protein interactions alter substantially the architecture of DNA-PK dimers at DSBs, to trigger subsequent interactions or enzymatic reactions.

post to: CiteULike

RNA self-cleavage activated by ultraviolet light-induced oxidation.

Publication Date: 2012 Jan 19 PMID: 21989404
Authors: Ariza-Mateos, A. - Prieto-Vega, S. - Diaz-Toledano, R. - Birk, A. - Szeto, H. - Mena, I. - Berzal-Herranz, A. - Gomez, J.
Journal: Nucleic Acids Res

A novel UV-C-light-induced ribozyme activity was discovered within the highly structured 5'-genomic regions of both Hepatitis C Virus (HCV) and the related Classic Swine Fever Virus (CSFV). Cleavage is mediated by exposure to UV-C light but not by exogenous oxygen radicals. It is also very selective, occurring at base positions HCV C(79) and CSFV A(45) in some molecules and at the immediately adjacent 5'-positions HCV U(78) and CSFV U(44) in others. Among other reaction products, the majority of biochemically active products detected contained 3'-phosphate and 5'-phosphate-end groups at the newly generated termini, along with a much lower amount of 3'-hydroxyl end group. While preservation of an E-loop RNA structure in the vicinity of the cleavage site was a requisite for HCV RNA self-cleavage, this was not the case for CSFV RNA. The short size of the reactive domains ( approximately 33 nt), which are compatible with primitive RNA motifs, and the lack of sequence homology, indicate that as-yet unidentified UV-activated ribozymes are likely to be found throughout structured RNAs, thereby providing clues to whether early RNA self-cleavage events were mediated by photosensitive RNA structures.

post to: CiteULike