ALBERT

Description: The advent in high-throughput-sequencing (HTS) technologies has revolutionized conventional biodiversity research by enabling parallel capture of DNA sequences possessing species-level diagnosis. However, polymerase chain reaction (PCR)-based implementation is biased by the efficiency of primer binding across lineages of organisms. A PCR-free HTS approach will alleviate this artefact and significantly improve upon the multi-locus method utilizing full mitogenomes. Here we developed a novel multiplex sequencing and assembly pipeline allowing for simultaneous acquisition of full mitogenomes from pooled animals without DNA enrichment or amplification. By concatenating assemblies from three de novo assemblers, we obtained high-quality mitogenomes for all 49 pooled taxa, with 36 species 〉15 kb and the remaining 〉10 kb, including 20 complete mitogenomes and nearly all protein coding genes (99.6%). The assembly quality was carefully validated with Sanger sequences, reference genomes and conservativeness of protein coding genes across taxa. The new method was effective even for closely related taxa, e.g. three Drosophila spp., demonstrating its broad utility for biodiversity research and mito-phylogenomics. Finally, the in silico simulation showed that by recruiting multiple mito-loci, taxon detection was improved at a fixed sequencing depth. Combined, these results demonstrate the plausibility of a multi-locus mito-metagenomics approach as the next phase of the current single-locus metabarcoding method.

Description: The efficacy and the mutation spectrum of genome editing methods can vary substantially depending on the targeted sequence. A simple, quick assay to accurately characterize and quantify the induced mutations is therefore needed. Here we present TIDE, a method for this purpose that requires only a pair of PCR reactions and two standard capillary sequencing runs. The sequence traces are then analyzed by a specially developed decomposition algorithm that identifies the major induced mutations in the projected editing site and accurately determines their frequency in a cell population. This method is cost-effective and quick, and it provides much more detailed information than current enzyme-based assays. An interactive web tool for automated decomposition of the sequence traces is available. TIDE greatly facilitates the testing and rational design of genome editing strategies.

Description: NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3–0.6 ppm and correlation coefficients ( r 2 ) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.

Description: Bicyclic oxazaphospholidine monomers were used to prepare a series of phosphorothioate (PS)-modified gapmer antisense oligonucleotides (ASOs) with control of the chirality of each of the PS linkages within the 10-base gap. The stereoselectivity was determined to be 98% for each coupling. The objective of this work was to study how PS chirality influences biophysical and biological properties of the ASO including binding affinity ( T m ), nuclease stability, activity in vitro and in vivo , RNase H activation and cleavage patterns (both human and E. coli ) in a gapmer context. Compounds that had nine or more Sp-linkages in the gap were found to be poorly active in vitro , while compounds with uniform Rp-gaps exhibited activity very similar to that of the stereo-random parent ASOs. Conversely, when tested in vivo , the full Rp-gap compound was found to be quickly metabolized resulting in low activity. A total of 31 ASOs were prepared with control of the PS chirally of each linkage within the gap in an attempt to identify favorable Rp/Sp positions. We conclude that a mix of Rp and Sp is required to achieve a balance between good activity and nuclease stability.

Description: In an isogenic cell population, phenotypic heterogeneity among individual cells is common and critical for survival of the population under different environment conditions. DNA modification is an important epigenetic factor that can regulate phenotypic heterogeneity. The single molecule real-time (SMRT) sequencing technology provides a unique platform for detecting a wide range of DNA modifications, including N6-methyladenine (6-mA), N4-methylcytosine (4-mC) and 5-methylcytosine (5-mC). Here we present qDNAmod, a novel bioinformatic tool for genome-wide quantitative profiling of intercellular heterogeneity of DNA modification from SMRT sequencing data. It is capable of estimating proportion of isogenic haploid cells, in which the same loci of the genome are differentially modified. We tested the reliability of qDNAmod with the SMRT sequencing data of Streptococcus pneumoniae strain ST556. qDNAmod detected extensive intercellular heterogeneity of DNA methylation (6-mA) in a clonal population of ST556. Subsequent biochemical analyses revealed that the recognition sequences of two type I restriction–modification (R-M) systems are responsible for the intercellular heterogeneity of DNA methylation initially identified by qDNAmod. qDNAmod thus represents a valuable tool for studying intercellular phenotypic heterogeneity from genome-wide DNA modification.

Description: Combinatorial transcription factor (TF) binding is essential for cell-type-specific gene regulation. However, much remains to be learned about the mechanisms of TF interactions, including to what extent constrained spacing and orientation of interacting TFs are critical for regulatory element activity. To examine the relative prevalence of the ‘enhanceosome’ versus the ‘TF collective’ model of combinatorial TF binding, a comprehensive analysis of TF binding site sequences in large scale datasets is necessary. We developed a motif-pair discovery pipeline to identify motif co-occurrences with preferential distance(s) between motifs in TF-bound regions. Utilizing a compendium of 289 mouse haematopoietic TF ChIP-seq datasets, we demonstrate that haematopoietic-related motif-pairs commonly occur with highly conserved constrained spacing and orientation between motifs. Furthermore, motif clustering revealed specific associations for both heterotypic and homotypic motif-pairs with particular haematopoietic cell types. We also showed that disrupting the spacing between motif-pairs significantly affects transcriptional activity in a well-known motif-pair—E-box and GATA, and in two previously unknown motif-pairs with constrained spacing—Ets and Homeobox as well as Ets and E-box. In this study, we provide evidence for widespread sequence-specific TF pair interaction with DNA that conforms to the ‘enhanceosome’ model, and furthermore identify associations between specific haematopoietic cell-types and motif-pairs.

Description: In silico tools have been developed to predict variants that may have an impact on pre-mRNA splicing. The major limitation of the application of these tools to basic research and clinical practice is the difficulty in interpreting the output. Most tools only predict potential splice sites given a DNA sequence without measuring splicing signal changes caused by a variant. Another limitation is the lack of large-scale evaluation studies of these tools. We compared eight in silico tools on 2959 single nucleotide variants within splicing consensus regions (scSNVs) using receiver operating characteristic analysis. The Position Weight Matrix model and MaxEntScan outperformed other methods. Two ensemble learning methods, adaptive boosting and random forests, were used to construct models that take advantage of individual methods. Both models further improved prediction, with outputs of directly interpretable prediction scores. We applied our ensemble scores to scSNVs from the Catalogue of Somatic Mutations in Cancer database. Analysis showed that predicted splice-altering scSNVs are enriched in recurrent scSNVs and known cancer genes. We pre-computed our ensemble scores for all potential scSNVs across the human genome, providing a whole genome level resource for identifying splice-altering scSNVs discovered from large-scale sequencing studies.

Description: In mammals, RNA interference is primarily a post-transcriptional mechanism. Evidence has accumulated for additional role in transcriptional gene silencing (TGS) but the question for a good paradigm for small interfering antigene RNA (agRNA)-induced chromatin modification remains unanswered. Here, we show that SETDB1, a histone H3-lysine 9 (H3K9)-specific methyltransferase, cooperates with Argonaute-2 (AGO2) and plays an essential role in agRNA-induced TGS. The androgen receptor ( AR ) gene was transcriptionally silenced by agRNA targeted to its promoter, and we show that this repression was mitigated by knockdown of SETDB1 or AGO2 . Chromatin immunoprecipitation demonstrated that agRNA-driven AGO2 was first targeted to the AR promoter, followed by SETDB1. SIN3A and HDAC1/2, the components of the SIN3-HDAC complex, immunoprecipitated with SETDB1, and localized at the agRNA-targeted promoter. Agreeing with the presence of SETDB1, trimethyl-H3K9 was enriched in the AR promoter. Both EZH2 and trimethyl-H3K27 were also present in the targeted locus; accordingly, EZH2 immunoprecipitated with SETDB1. DNA methylation level was not significantly changed, suggesting the absence of de novo methylating activity in agRNA-induced AR promoter. Our results demonstrate that SETDB1, together with AGO2, plays an essential role in TGS through recruiting chromatin remodeler and/or other modifiers, consequently creating a repressive chromatin milieu at the targeted promoter.

Description: LEF/TCFs direct the final step in Wnt/β-catenin signalling by recruiting β-catenin to genes for activation of transcription. Ancient, non-vertebrate TCFs contain two DNA binding domains, a High Mobility Group box for recognition of the Wnt Response Element (WRE; 5'-CTTTGWWS-3') and the C-clamp domain for recognition of the GC-rich Helper motif (5'-RCCGCC-3'). Two vertebrate TCFs (TCF-1/TCF7 and TCF-4/TCF7L2) use the C-clamp as an alternatively spliced domain to regulate cell-cycle progression, but how the C-clamp influences TCF binding and activity genome-wide is not known. Here, we used a doxycycline inducible system with ChIP-seq to assess how the C-clamp influences human TCF1 binding genome-wide. Metabolic pulse-labeling of nascent RNA with 4'Thiouridine was used with RNA-seq to connect binding to the Wnt transcriptome. We find that the C-clamp enables targeting to a greater number of gene loci for stronger occupancy and transcription regulation. The C-clamp uses Helper sites concurrently with WREs for gene targeting, but it also targets TCF1 to sites that do not have readily identifiable canonical WREs. The coupled ChIP-seq/4'Thiouridine-seq analysis identified new Wnt target genes, including additional regulators of cell proliferation. Thus, C-clamp containing isoforms of TCFs are potent transcriptional regulators with an expanded transcriptome directed by C-clamp-Helper site interactions.

Description: Chromatin modifiers and histone modifications are components of a chromatin-signaling network involved in transcription and its regulation. The interactions between chromatin modifiers and histone modifications are often unknown, are based on the analysis of few genes or are studied in vitro . Here, we apply computational methods to recover interactions between chromatin modifiers and histone modifications from genome-wide ChIP-Seq data. These interactions provide a high-confidence backbone of the chromatin-signaling network. Many recovered interactions have literature support; others provide hypotheses about yet unknown interactions. We experimentally verified two of these predicted interactions, leading to a link between H4K20me1 and members of the Polycomb Repressive Complexes 1 and 2. Our results suggest that our computationally derived interactions are likely to lead to novel biological insights required to establish the connectivity of the chromatin-signaling network involved in transcription and its regulation.

Description: As the biomedical impact of small RNAs grows, so does the need to understand competing structural alternatives for regions of functional interest. Suboptimal structure analysis provides significantly more RNA base pairing information than a single minimum free energy prediction. Yet computational enhancements like Boltzmann sampling have not been fully adopted by experimentalists since identifying meaningful patterns in this data can be challenging. Profiling is a novel approach to mining RNA suboptimal structure data which makes the power of ensemble-based analysis accessible in a stable and reliable way. Balancing abstraction and specificity, profiling identifies significant combinations of base pairs which dominate low-energy RNA secondary structures. By design, critical similarities and differences are highlighted, yielding crucial information for molecular biologists. The code is freely available via http://gtfold.sourceforge.net/profiling.html .

Description: The stress-sensitive restriction-modification (RM) system CglI from Corynebacterium glutamicum and the homologous NgoAVII RM system from Neisseria gonorrhoeae FA1090 are composed of three genes: a DNA methyltransferase (M.CglI and M.NgoAVII), a putative restriction endonuclease (R.CglI and R.NgoAVII, or R-proteins) and a predicted DEAD-family helicase/ATPase (N.CglI and N.NgoAVII or N-proteins). Here we report a biochemical characterization of the R- and N-proteins. Size-exclusion chromatography and SAXS experiments reveal that the isolated R.CglI, R.NgoAVII and N.CglI proteins form homodimers, while N.NgoAVII is a monomer in solution. Moreover, the R.CglI and N.CglI proteins assemble in a complex with R 2 N 2 stoichiometry. Next, we show that N-proteins have ATPase activity that is dependent on double-stranded DNA and is stimulated by the R-proteins. Functional ATPase activity and extensive ATP hydrolysis (~170 ATP/s/monomer) are required for site-specific DNA cleavage by R-proteins. We show that ATP-dependent DNA cleavage by R-proteins occurs at fixed positions (6–7 nucleotides) downstream of the asymmetric recognition sequence 5'-GCCGC-3'. Despite similarities to both Type I and II restriction endonucleases, the CglI and NgoAVII enzymes may employ a unique catalytic mechanism for DNA cleavage.

Description: Emerging evidence points to roles for tRNA modifications and tRNA abundance in cellular stress responses. While isolated instances of stress-induced tRNA degradation have been reported, we sought to assess the effects of stress on tRNA levels at a systems level. To this end, we developed a next-generation sequencing method that exploits the paucity of ribonucleoside modifications at the 3'-end of tRNAs to quantify changes in all cellular tRNA molecules. Application of this tRNA-seq method to Saccharomyces cerevisiae identified all 76 expressed unique tRNA species out of 295 coded in the yeast genome, including all isoacceptor variants, with highly precise relative (fold-change) quantification of tRNAs. In studies of stress-induced changes in tRNA levels, we found that oxidation (H 2 O 2 ) and alkylation (methylmethane sulfonate, MMS) stresses induced nearly identical patterns of up- and down-regulation for 58 tRNAs. However, 18 tRNAs showed opposing changes for the stresses, which parallels our observation of signature reprogramming of tRNA modifications caused by H 2 O 2 and MMS. Further, stress-induced degradation was limited to only a small proportion of a few tRNA species. With tRNA-seq applicable to any organism, these results suggest that translational control of stress response involves a contribution from tRNA abundance.

Description: During meiosis programmed DNA double-strand breaks (DSBs) are repaired by homologous recombination using the sister chromatid or the homologous chromosome (homolog) as a template. This repair results in crossover (CO) and non-crossover (NCO) recombinants. Only CO formation between homologs provides the physical linkages guiding correct chromosome segregation, which are essential to produce healthy gametes. The factors that determine the CO/NCO decision are still poorly understood. Using Schizosaccharomyces pombe as a model we show that the Rad51/Dmc1-paralog complexes Rad55-Rad57 and Rdl1-Rlp1-Sws1 together with Swi5-Sfr1 play a major role in antagonizing both the FANCM-family DNA helicase/translocase Fml1 and the RecQ-type DNA helicase Rqh1 to limit hybrid DNA formation and promote Mus81-Eme1-dependent COs. A common attribute of these protein complexes is an ability to stabilize the Rad51/Dmc1 nucleoprotein filament, and we propose that it is this property that imposes constraints on which enzymes gain access to the recombination intermediate, thereby controlling the manner in which it is processed and resolved.

Description: To study target sequence specificity, selectivity, and reaction kinetics of Streptococcus pyogenes Cas9 activity, we challenged libraries of random variant targets with purified Cas9::guide RNA complexes in vitro . Cleavage kinetics were nonlinear, with a burst of initial activity followed by slower sustained cleavage. Consistent with other recent analyses of Cas9 sequence specificity, we observe considerable (albeit incomplete) impairment of cleavage for targets mutated in the PAM sequence or in ‘seed’ sequences matching the proximal 8 bp of the guide. A second target region requiring close homology was located at the other end of the guide::target duplex (positions 13–18 relative to the PAM). Sequences flanking the guide+PAM region had measurable (albeit modest) effects on cleavage. In addition, the first-base Guanine constraint commonly imposed by gRNA expression systems has little effect on overall cleavage efficiency. Taken together, these studies provide an in vitro understanding of the complexities of Cas9–gRNA interaction and cleavage beyond the general paradigm of site determination based on the ‘seed’ sequence and PAM.

Description: Transcripts possessing a 5'-triphosphate are a hallmark of viral transcription and can trigger the host antiviral response. 5'-triphosphates are also found on common host transcripts transcribed by RNA polymerase III (RNAP III), yet how these transcripts remain non-immunostimulatory is incompletely understood. Most microRNAs (miRNAs) are 5'-monophosphorylated as a result of sequential endonucleolytic processing by Drosha and Dicer from longer RNA polymerase II (RNAP II)-transcribed primary transcripts. In contrast, bovine leukemia virus (BLV) expresses subgenomic RNAP III transcripts that give rise to miRNAs independent of Drosha processing. Here, we demonstrate that each BLV pre-miRNA is directly transcribed by RNAP III from individual, compact RNAP III type II genes. Thus, similar to manmade RNAP III-generated short hairpin RNAs (shRNAs), the BLV pre-miRNAs are initially 5'-triphosphorylated. Nonetheless, the derivative 5p miRNAs and shRNA-generated 5p small RNAs (sRNAs) possess a 5'-monophosphate. Our enzymatic characterization and small RNA sequencing data demonstrate that BLV 5p miRNAs are co-terminal with 5'-triphosphorylated miRNA precursors (pre-miRNAs). Thus, these results identify a 5'-tri-phosphatase activity that is involved in the biogenesis of BLV miRNAs and shRNA-generated sRNAs. This work advances our understanding of retroviral miRNA and shRNA biogenesis and may have implications regarding the immunostimulatory capacity of RNAP III transcripts.

Description: Within the field of synthetic biology, a rational design of genetic parts should include a causal understanding of their input-output responses—the so-called transfer function—and how to tune them. However, a commonly adopted strategy is to fit data to Hill-shaped curves without considering the underlying molecular mechanisms. Here we provide a novel mathematical formalization that allows prediction of the global behavior of a synthetic device by considering the actual information from the involved biological parts. This is achieved by adopting an enzymology-like framework, where transfer functions are described in terms of their input affinity constant and maximal response. As a proof of concept, we characterize a set of Lux homoserine-lactone-inducible genetic devices with different levels of Lux receptor and signal molecule. Our model fits the experimental results and predicts the impact of the receptor's ribosome-binding site strength, as a tunable parameter that affects gene expression. The evolutionary implications are outlined.

Description: The restriction endonuclease (REase) NgoAVII is composed of two proteins, R.NgoAVII and N.NgoAVII, and shares features of both Type II restriction enzymes and Type I/III ATP-dependent restriction enzymes (see accompanying paper Zaremba et al. , 2014). Here we present crystal structures of the R.NgoAVII apo-protein and the R.NgoAVII C-terminal domain bound to a specific DNA. R.NgoAVII is composed of two domains: an N-terminal nucleolytic PLD domain; and a C-terminal B3-like DNA-binding domain identified previously in BfiI and EcoRII REases, and in plant transcription factors. Structural comparison of the B3-like domains of R.NgoAVII, EcoRII, BfiI and the plant transcription factors revealed a conserved DNA-binding surface comprised of N- and C-arms that together grip the DNA. The C-arms of R.NgoAVII, EcoRII, BfiI and plant B3 domains are similar in size, but the R.NgoAVII N-arm which makes the majority of the contacts to the target site is much longer. The overall structures of R.NgoAVII and BfiI are similar; however, whilst BfiI has stand-alone catalytic activity, R.NgoAVII requires an auxiliary cognate N.NgoAVII protein and ATP hydrolysis in order to cleave DNA at the target site. The structures we present will help formulate future experiments to explore the molecular mechanisms of intersubunit crosstalk that control DNA cleavage by R.NgoAVII and related endonucleases.

Description: Distinct translational initiation mechanisms between prokaryotes and eukaryotes limit the exploitation of prokaryotic riboswitch repertoire for regulatory RNA circuit construction in mammalian application. Here, we explored programmed ribosomal frameshifting (PRF) as the regulatory gene expression platform for engineered ligand-responsive RNA devices in higher eukaryotes. Regulation was enabled by designed ligand-dependent conformational rearrangements of the two cis-acting RNA motifs of opposite activity in -1 PRF. Particularly, RNA elements responsive to trans-acting ligands can be tailored to modify co-translational RNA refolding dynamics of a hairpin upstream of frameshifting site to achieve reversible and adjustable -1 PRF attenuating activity. Combined with a ligand-responsive stimulator, synthetic RNA devices for synergetic translational-elongation control of gene expression can be constructed. Due to the similarity between co-transcriptional RNA hairpin folding and co-translational RNA hairpin refolding, the RNA-responsive ligand repertoire provided in prokaryotic systems thus becomes accessible to gene-regulatory circuit construction for synthetic biology application in mammalian cells.

Description: Cutaneous photosensitization is a common side effect of drug treatment and can be associated with an increased skin cancer risk. The immunosuppressant azathioprine, the fluoroquinolone antibiotics and vemurafenib—a BRAF inhibitor used to treat metastatic melanoma—are all recognized clinical photosensitizers. We have compared the effects of UVA radiation on cultured human cells treated with 6-thioguanine (6-TG, a DNA-embedded azathioprine surrogate), the fluoroquinolones ciprofloxacin and ofloxacin and vemurafenib. Despite widely different structures and modes of action, each of these drugs potentiated UVA cytotoxicity. UVA photoactivation of 6-TG, ciprofloxacin and ofloxacin was associated with the generation of singlet oxygen that caused extensive protein oxidation. In particular, these treatments were associated with damage to DNA repair proteins that reduced the efficiency of nucleotide excision repair. Although vemurafenib was also highly phototoxic to cultured cells, its effects were less dependent on singlet oxygen. Highly toxic combinations of vemurafenib and UVA caused little protein carbonylation but were nevertheless inhibitory to nucleotide excision repair. Thus, for three different classes of drugs, photosensitization by at least two distinct mechanisms is associated with reduced protection against potentially mutagenic and carcinogenic DNA damage.

Description: The first level of genome packaging in eukaryotic cells involves the formation of dense nucleosome arrays, with DNA coverage near 90% in yeasts. How cells achieve such high coverage within a short time, e.g. after DNA replication, remains poorly understood. It is known that random sequential adsorption of impenetrable particles on a line reaches high density extremely slowly, due to a jamming phenomenon. The nucleosome-shifting action of remodeling enzymes has been proposed as a mechanism to resolve such jams. Here, we suggest two biophysical mechanisms which assist rapid filling of DNA with nucleosomes, and we quantitatively characterize these mechanisms within mathematical models. First, we show that the ‘softness’ of nucleosomes, due to nucleosome breathing and stepwise nucleosome assembly, significantly alters the filling behavior, speeding up the process relative to ‘hard’ particles with fixed, mutually exclusive DNA footprints. Second, we explore model scenarios in which the progression of the replication fork could eliminate nucleosome jamming, either by rapid filling in its wake or via memory of the parental nucleosome positions. Taken together, our results suggest that biophysical effects promote rapid nucleosome filling, making the reassembly of densely packed nucleosomes after DNA replication a simpler task for cells than was previously thought.

Description: Abasic (AP) lesions are the most frequent type of damages occurring in cellular DNA. Here we describe the conformational effects of AP sites substituted for 2'-deoxyadenosine in the first ( ap7 ), second ( ap13 ) or third ( ap19 ) loop of the quadruplex formed in K + by the human telomere DNA 5'-d[AG 3 (TTAG 3 ) 3 ]. CD spectra and electrophoresis reveal that the presence of AP sites does not hinder the formation of intramolecular quadruplexes. NMR spectra show that the structural heterogeneity is substantially reduced in ap7 and ap19 as compared to that in the wild-type. These two ( ap7 and ap19 ) sequences are shown to adopt the hybrid-1 and hybrid-2 quadruplex topology, respectively, with AP site located in a propeller-like loop. All three studied sequences transform easily into parallel quadruplex in dehydrating ethanol solution. Thus, the AP site in any loop region facilitates the formation of the propeller loop. Substitution of all adenines by AP sites stabilizes the parallel quadruplex even in the absence of ethanol. Whereas guanines are the major determinants of quadruplex stability, the presence or absence of loop adenines substantially influences quadruplex folding. The naturally occurring adenine-lacking sites in the human telomere DNA can change the quadruplex topology in vivo with potentially vital biological consequences.

Description: The DNA sequence preferences of nearly all sequence specific DNA binding proteins are influenced by the identities of bases that are not directly contacted by protein. Discrimination between non-contacted base sequences is commonly based on the differential abilities of DNA sequences to allow narrowing of the DNA minor groove. However, the factors that govern the propensity of minor groove narrowing are not completely understood. Here we show that the differential abilities of various DNA sequences to support formation of a highly ordered and stable minor groove solvation network are a key determinant of non-contacted base recognition by a sequence-specific binding protein. In addition, disrupting the solvent network in the non-contacted region of the binding site alters the protein's ability to recognize contacted base sequences at positions 5–6 bases away. This observation suggests that DNA solvent interactions link contacted and non-contacted base recognition by the protein.

Description: In recent years, an increasing number of reports have been focused on the structure and biological role of non-canonical nucleic acid secondary structures. Many of these studies involve the use of oligonucleotides that can often adopt a variety of structures depending on the experimental conditions, and hence change the outcome of an assay. The knowledge of the structure(s) formed by oligonucleotides is thus critical to correctly interpret the results, and gain insight into the biological role of these particular sequences. Herein we demonstrate that size-exclusion HPLC (SE-HPLC) is a simple yet surprisingly powerful tool to quickly and effortlessly assess the secondary structure(s) formed by oligonucleotides. For the first time, an extensive calibration and validation of the use of SE-HPLC to confidently detect the presence of different species displaying various structure and/or molecularity, involving 〉110 oligonucleotides forming a variety of secondary structures (antiparallel, parallel, A-tract bent and mismatched duplexes, triplexes, G-quadruplexes and i-motifs, RNA stem loops), is performed. Moreover, we introduce simple metrics that allow the use of SE-HPLC without the need for a tedious calibration work. We show that the remarkable versatility of the method allows to quickly establish the influence of a number of experimental parameters on nucleic acid structuration and to operate on a wide range of oligonucleotide concentrations. Case studies are provided to clearly illustrate the all-terrain capabilities of SE-HPLC for oligonucleotide secondary structure analysis. Finally, this manuscript features a number of important observations contributing to a better understanding of nucleic acid structural polymorphism.

Description: Single-molecule manipulation (SMM) techniques use applied force, and measured elastic response, to reveal microscopic physical parameters of individual biomolecules and details of biomolecular interactions. A major hurdle in the application of these techniques is the labeling method needed to immobilize biomolecules on solid supports. A simple, minimally-perturbative labeling strategy would significantly broaden the possible applications of SMM experiments, perhaps even allowing the study of native biomolecular structures. To accomplish this, we investigate the use of functionalized locked nucleic acid (LNA) oligomers as biomolecular handles that permit sequence-specific binding and immobilization of DNA. We find these probes form bonds with DNA with high specificity but with varied stability in response to the direction of applied mechanical force: when loaded in a shear orientation, the bound LNA oligomers were measured to be two orders of magnitude more stable than when loaded in a peeling, or unzipping, orientation. Our results show that LNA provides a simple, stable means to functionalize dsDNA for manipulation. We provide design rules that will facilitate their use in future experiments.

Description: Development of single-molecule localization microscopy techniques has allowed nanometre scale localization accuracy inside cells, permitting the resolution of ultra-fine cell structure and the elucidation of crucial molecular mechanisms. Application of these methodologies to understanding processes underlying DNA replication and repair has been limited to defined in vitro biochemical analysis and prokaryotic cells. In order to expand these techniques to eukaryotic systems, we have further developed a photo-activated localization microscopy-based method to directly visualize DNA-associated proteins in unfixed eukaryotic cells. We demonstrate that motion blurring of fluorescence due to protein diffusivity can be used to selectively image the DNA-bound population of proteins. We designed and tested a simple methodology and show that it can be used to detect changes in DNA binding of a replicative helicase subunit, Mcm4, and the replication sliding clamp, PCNA, between different stages of the cell cycle and between distinct genetic backgrounds.

Description: Decitabine (5-aza-2'-deoxycytidine) is a DNA methyltransferase inhibitor and an archetypal epigenetic drug for the therapy of myeloid leukemias. The mode of action of decitabine strictly depends on the incorporation of the drug into DNA. However, DNA incorporation and ensuing genotoxic effects of decitabine have not yet been investigated in human cancer cell lines or in models related to the approved indication of the drug. Here we describe a robust assay for the quantitative determination of decitabine incorporation rates into DNA from human cancer cells. Using a panel of human myeloid leukemia cell lines we show appreciable amounts of decitabine incorporation that closely correlated with cellular drug uptake. Decitabine incorporation was also detectable in primary cells from myeloid leukemia patients, indicating that the assay is suitable for biomarker analyses to predict drug responses in patients. Finally, we also used next-generation sequencing to comprehensively analyze the effects of decitabine incorporation on the DNA sequence level. Interestingly, this approach failed to reveal significant changes in the rates of point mutations and genome rearrangements in myeloid leukemia cell lines. These results indicate that standard rates of decitabine incorporation are not genotoxic in myeloid leukemia cells.

Description: A new functional gene database, FOAM (Functional Ontology Assignments for Metagenomes), was developed to screen environmental metagenomic sequence datasets. FOAM provides a new functional ontology dedicated to classify gene functions relevant to environmental microorganisms based on Hidden Markov Models (HMMs). Sets of aligned protein sequences (i.e. ‘profiles’) were tailored to a large group of target KEGG Orthologs (KOs) from which HMMs were trained. The alignments were checked and curated to make them specific to the targeted KO. Within this process, sequence profiles were enriched with the most abundant sequences available to maximize the yield of accurate classifier models. An associated functional ontology was built to describe the functional groups and hierarchy. FOAM allows the user to select the target search space before HMM-based comparison steps and to easily organize the results into different functional categories and subcategories. FOAM is publicly available at http://portal.nersc.gov/project/m1317/FOAM/ .

Description: The understanding of folding and function of RNA molecules depends on the identification and classification of interactions between ribonucleotide residues. We developed a new method named ClaRNA for computational classification of contacts in RNA 3D structures. Unique features of the program are the ability to identify imperfect contacts and to process coarse-grained models. Each doublet of spatially close ribonucleotide residues in a query structure is compared to clusters of reference doublets obtained by analysis of a large number of experimentally determined RNA structures, and assigned a score that describes its similarity to one or more known types of contacts, including pairing, stacking, base–phosphate and base–ribose interactions. The accuracy of ClaRNA is 0.997 for canonical base pairs, 0.983 for non-canonical pairs and 0.961 for stacking interactions. The generalized squared correlation coefficient (GC 2 ) for ClaRNA is 0.969 for canonical base pairs, 0.638 for non-canonical pairs and 0.824 for stacking interactions. The classifier can be easily extended to include new types of spatial relationships between pairs or larger assemblies of nucleotide residues. ClaRNA is freely available via a web server that includes an extensive set of tools for processing and visualizing structural information about RNA molecules.

Description: The Prp19-associated complex is required for spliceosome activation by stabilizing the binding of U5 and U6 on the spliceosome after the release of U4. The complex comprises at least eight proteins, among which Ntc90 and Ntc77 contain multiple tetratricopeptide repeat (TPR) elements. We have previously shown that Ntc90 is not involved in spliceosome activation, but is required for the recruitment of essential first-step factor Yju2 to the spliceosome. We demonstrate here that Ntc77 has dual functions in both spliceosome activation and the first catalytic step in recruiting Yju2. We have identified an amino-terminal region of Ntc77, which encompasses the N-terminal domain and the first three TPR motifs, dispensable for spliceosome activation but required for stable interaction of Yju2 with the spliceosome. Deletion of this region had no severe effect on the integrity of the NTC, binding of NTC to the spliceosome or spliceosome activation, but impaired splicing and exhibited a dominant-negative growth phenotype. Our data reveal functional roles of Ntc77 in both spliceosome activation and the first catalytic step, and distinct structural domains of Ntc77 required for these two steps.

Description: Here we conducted an integrative multi-omics analysis to understand how cancers harbor various types of aberrations at the genomic, epigenomic and transcriptional levels. In order to elucidate biological relevance of the aberrations and their mutual relations, we performed whole-genome sequencing, RNA-Seq, bisulfite sequencing and ChIP-Seq of 26 lung adenocarcinoma cell lines. The collected multi-omics data allowed us to associate an average of 536 coding mutations and 13,573 mutations in promoter or enhancer regions with aberrant transcriptional regulations. We detected the 385 splice site mutations and 552 chromosomal rearrangements, representative cases of which were validated to cause aberrant transcripts. Averages of 61, 217, 3687 and 3112 mutations are located in the regulatory regions which showed differential DNA methylation, H3K4me3, H3K4me1 and H3K27ac marks, respectively. We detected distinct patterns of aberrations in transcriptional regulations depending on genes. We found that the irregular histone marks were characteristic to EGFR and CDKN1A, while a large genomic deletion and hyper-DNA methylation were most frequent for CDKN2A. We also used the multi-omics data to classify the cell lines regarding their hallmarks of carcinogenesis. Our datasets should provide a valuable foundation for biological interpretations of interlaced genomic and epigenomic aberrations.

Description: Alkylative damage to DNA can be induced by environmental chemicals, endogenous metabolites and some commonly prescribed chemotherapeutic agents. The regioisomeric N 3-, O 2 - and O 4 -ethylthymidine ( N 3-, O 2 - and O 4 -EtdT, respectively) represent an important class of ethylated DNA lesions. Using nonreplicative double-stranded vectors containing an N 3-EtdT, O 2 -EtdT or O 4 -EtdT at a defined site in the template strand, herein we examined the effects of these lesions on DNA transcription mediated by single-subunit T7 RNA polymerase or multisubunit human RNA polymerase II in vitro and in human cells. We found that O 4 -EtdT is highly mutagenic and exclusively induces the misincorporation of guanine opposite the lesion, whereas N 3-EtdT and O 2 -EtdT display promiscuous miscoding properties during transcription. In addition, N 3-EtdT and O 2 -EtdT were found to inhibit strongly DNA transcription in vitro and in certain human cells. Moreover, N 3-EtdT, but not O 2 -EtdT or O 4 -EtdT, is an efficient substrate for transcription-coupled nucleotide excision repair. These findings provide new important insights into how these alkylated DNA lesions compromise the flow of genetic information, which may help to understand the risk of these lesions in living cells.

Description: In the liver Wnt-signaling contributes to the metabolic fate of hepatocytes, but the precise role of the TCF7L2 in this process is unknown. We employed a temporal RNA-Seq approach to examine gene expression 3–96 h following Tcf7l2 silencing in rat hepatoma cells, and combined this with ChIP-Seq to investigate mechanisms of target gene regulation by TCF7L2. Silencing Tcf7l2 led to a time-dependent appearance of 406 differentially expressed genes (DEGs), including key regulators of cellular growth and differentiation, and amino acid, lipid and glucose metabolism. Direct regulation of 149 DEGs was suggested by strong proximal TCF7L2 binding (peak proximity score 〉 10) and early mRNA expression changes (≤18 h). Indirect gene regulation by TCF7L2 likely occurred via alternate transcription factors, including Hnf4a , Foxo1 , Cited2 , Myc and Lef1 , which were differentially expressed following Tcf7l2 knock-down. Tcf7l2 -silencing enhanced the expression and chromatin occupancy of HNF4α, and co-siRNA experiments revealed that HNF4α was required for the regulation of a subset of metabolic genes by TCF7L2, particularly those involved in lipid and amino-acid metabolism. Our findings suggest TCF7L2 is an important regulator of the hepatic phenotype, and highlight novel mechanisms of gene regulation by TCF7L2 that involve interplay between multiple hepatic transcriptional pathways.

Description: 5S Ribosomal RNA (5S rRNA) is a universal component of ribosomes, and the corresponding gene is easily identified in archaeal, bacterial and nuclear genome sequences. However, organelle gene homologs ( rrn5 ) appear to be absent from most mitochondrial and several chloroplast genomes. Here, we re-examine the distribution of organelle rrn5 by building mitochondrion- and plastid-specific covariance models (CMs) with which we screened organelle genome sequences. We not only recover all organelle rrn5 genes annotated in GenBank records, but also identify more than 50 previously unrecognized homologs in mitochondrial genomes of various stramenopiles, red algae, cryptomonads, malawimonads and apusozoans, and surprisingly, in the apicoplast (highly derived plastid) genomes of the coccidian pathogens Toxoplasma gondii and Eimeria tenella . Comparative modeling of RNA secondary structure reveals that mitochondrial 5S rRNAs from brown algae adopt a permuted triskelion shape that has not been seen elsewhere. Expression of the newly predicted rrn5 genes is confirmed experimentally in 10 instances, based on our own and published RNA-Seq data. This study establishes that particularly mitochondrial 5S rRNA has a much broader taxonomic distribution and a much larger structural variability than previously thought. The newly developed CMs will be made available via the Rfam database and the MFannot organelle genome annotator.

Description: Mammalian splicing regulatory protein RNA-binding motif protein 4 (RBM4) has an alanine repeat-containing C-terminal domain (CAD) that confers both nuclear- and splicing speckle-targeting activities. Alanine-repeat expansion has pathological potential. Here we show that the alanine-repeat tracts influence the subnuclear targeting properties of the RBM4 CAD in cultured human cells. Notably, truncation of the alanine tracts redistributed a portion of RBM4 to paraspeckles. The alanine-deficient CAD was sufficient for paraspeckle targeting. On the other hand, alanine-repeat expansion reduced the mobility of RBM4 and impaired its splicing activity. We further took advantage of the putative coactivator activator (CoAA)-RBM4 conjoined splicing factor, CoAZ, to investigate the function of the CAD in subnuclear targeting. Transiently expressed CoAZ formed discrete nuclear foci that emerged and subsequently separated—fully or partially—from paraspeckles. Alanine-repeat expansion appeared to prevent CoAZ separation from paraspeckles, resulting in their complete colocalization. CoAZ foci were dynamic but, unlike paraspeckles, were resistant to RNase treatment. Our results indicate that the alanine-rich CAD, in conjunction with its conjoined RNA-binding domain(s), differentially influences the subnuclear localization and biogenesis of RBM4 and CoAZ.

Description: Efficient transcription of the HIV-1 genome is regulated by Tat, which recruits P-TEFb from the 7SK small nuclear ribonucleoprotein (snRNP) and other nucleoplasmic complexes to phosphorylate RNA polymerase II and other factors associated with the transcription complex. Although Tat activity is dependent on its binding to the viral TAR sequence, little is known about the cellular factors that might also assemble onto this region of the viral transcript. Here, we report that the splicing factor SRSF1 (SF2/ASF) and Tat recognize overlapping sequences within TAR and the 7SK RNA. SRSF1 expression can inhibit Tat transactivation by directly competing for its binding to TAR. Additionally, we provide evidence that SRSF1 can increase the basal level of viral transcription in the absence of Tat. We propose that SRSF1 activates transcription in the early stages of viral infection by recruiting P-TEFb to TAR from the 7SK snRNP. Whereas in the later stages, Tat substitutes for SRSF1 by promoting release of the stalled polymerase and more efficient transcriptional elongation.

Description: We describe the identification and characterization of novel homing endonucleases using genome database mining to identify putative target sites, followed by high throughput activity screening in a bacterial selection system. We characterized the substrate specificity and kinetics of these endonucleases by monitoring DNA cleavage events with deep sequencing. The endonuclease specificities revealed by these experiments can be partially recapitulated using 3D structure-based computational models. Analysis of these models together with genome sequence data provide insights into how alternative endonuclease specificities were generated during natural evolution.

Description: Mtr4 is a conserved Ski2-like RNA helicase and a subunit of the TRAMP complex that activates exosome-mediated 3'-5' turnover in nuclear RNA surveillance and processing pathways. Prominent features of the Mtr4 structure include a four-domain ring-like helicase core and a large arch domain that spans the core. The ‘ratchet helix’ is positioned to interact with RNA substrates as they move through the helicase. However, the contribution of the ratchet helix in Mtr4 activity is poorly understood. Here we show that strict conservation along the ratchet helix is particularly extensive for Ski2-like RNA helicases compared to related helicases. Mutation of residues along the ratchet helix alters in vitro activity in Mtr4 and TRAMP and causes slow growth phenotypes in vivo . We also identify a residue on the ratchet helix that influences Mtr4 affinity for polyadenylated substrates. Previous work indicated that deletion of the arch domain has minimal effect on Mtr4 unwinding activity. We now show that combining the arch deletion with ratchet helix mutations abolishes helicase activity and produces a lethal in vivo phenotype. These studies demonstrate that the ratchet helix modulates helicase activity and suggest that the arch domain plays a previously unrecognized role in unwinding substrates.

Description: Recent evidence points to a role of chromatin in regulation of alternative pre-mRNA splicing (AS). In order to identify novel chromatin regulators of AS, we screened an RNAi library of chromatin proteins using a cell-based high-throughput in vivo assay. We identified a set of chromatin proteins that regulate AS. Using simultaneous genome-wide expression and AS analysis, we demonstrate distinct and non-overlapping functions of these chromatin modifiers on transcription and AS. Detailed mechanistic characterization of one dual function chromatin modifier, the H3K9 methyltransferase EHMT2 (G9a), identified VEGFA as a major chromatin-mediated AS target. Silencing of EHMT2, or its heterodimer partner EHMT1, affects AS by promoting exclusion of VEGFA exon 6a, but does not alter total VEGFA mRNA levels. The epigenetic regulatory mechanism of AS by EHMT2 involves an adaptor system consisting of the chromatin modulator HP1, which binds methylated H3K9 and recruits splicing regulator SRSF1. The epigenetic regulation of VEGFA is physiologically relevant since EHMT2 is transcriptionally induced in response to hypoxia and triggers concomitant changes in AS of VEGFA. These results characterize a novel epigenetic regulatory mechanism of AS and they demonstrate separate roles of epigenetic modifiers in transcription and alternative splicing.

Description: N 6 A methylation is the most abundant RNA modification occurring within messenger RNA. Impairment of methylase or demethylase functions are associated with severe phenotypes and diseases in several organisms. Beside writer and eraser enzymes of this dynamic RNA epigenetic modification, reader proteins that recognize this modification are involved in numerous cellular processes. Although the precise characterization of these reader proteins remains unknown, preliminary data showed that most potential reader proteins contained a conserved YT521-B homology (YTH) domain. Here we define the YTH domain of rat YT521-B as a N 6 -methylated adenosine reader domain and report its solution structure in complex with a N 6 -methylated RNA. The structure reveals a binding preference for NGANNN RNA hexamer and a deep hydrophobic cleft for m 6 A recognition. These findings establish a molecular function for YTH domains as m 6 A reader domains and should guide further studies into the biological functions of YTH-containing proteins in m 6 A recognition.

Description: Homologous non-coding RNAs frequently exhibit domain insertions, where a branch of secondary structure is inserted in a sequence with respect to its homologs. Dynamic programming algorithms for common secondary structure prediction of multiple RNA homologs, however, do not account for these domain insertions. This paper introduces a novel dynamic programming algorithm methodology that explicitly accounts for the possibility of inserted domains when predicting common RNA secondary structures. The algorithm is implemented as Dynalign II, an update to the Dynalign software package for predicting the common secondary structure of two RNA homologs. This update is accomplished with negligible increase in computational cost. Benchmarks on ncRNA families with domain insertions validate the method. Over base pairs occurring in inserted domains, Dynalign II improves accuracy over Dynalign, attaining 80.8% sensitivity (compared with 14.4% for Dynalign) and 91.4% positive predictive value (PPV) for tRNA; 66.5% sensitivity (compared with 38.9% for Dynalign) and 57.0% PPV for RNase P RNA; and 50.1% sensitivity (compared with 24.3% for Dynalign) and 58.5% PPV for SRP RNA. Compared with Dynalign, Dynalign II also exhibits statistically significant improvements in overall sensitivity and PPV. Dynalign II is available as a component of RNAstructure, which can be downloaded from http://rna.urmc.rochester.edu/RNAstructure.html .

Description: Fanconi anaemia (FA) is a genome instability disease caused by defects in the FA DNA repair pathway that senses and repairs damage caused by DNA interstrand crosslinks. At least 8 of the 16 genes found mutated in FA encode proteins that assemble into the FA core complex, a multisubunit monoubiquitin E3 ligase. Here, we show that the RuvBL1 and RuvBL2 AAA+ ATPases co-purify with FA core complex isolated under stringent but native conditions from a vertebrate cell line. Depletion of the RuvBL1-RuvBL2 complex in human cells causes hallmark features of FA including DNA crosslinker sensitivity, chromosomal instability and defective FA pathway activation. Genetic knockout of RuvBL1 in a murine model is embryonic lethal while conditional inactivation in the haematopoietic stem cell pool confers profound aplastic anaemia. Together these findings reveal a function for RuvBL1-RuvBL2 in DNA repair through a physical and functional association with the FA core complex. Surprisingly, depletion of RuvBL1-RuvBL2 leads to co-depletion of the FA core complex in human cells. This suggests that a potential mechanism for the role of RuvBL1-RuvBL2 in maintaining genome integrity is through controlling the cellular abundance of FA core complex.

Description: It is now known that unwanted noise and unmodeled artifacts such as batch effects can dramatically reduce the accuracy of statistical inference in genomic experiments. These sources of noise must be modeled and removed to accurately measure biological variability and to obtain correct statistical inference when performing high-throughput genomic analysis. We introduced surrogate variable analysis (sva) for estimating these artifacts by (i) identifying the part of the genomic data only affected by artifacts and (ii) estimating the artifacts with principal components or singular vectors of the subset of the data matrix. The resulting estimates of artifacts can be used in subsequent analyses as adjustment factors to correct analyses. Here I describe a version of the sva approach specifically created for count data or FPKMs from sequencing experiments based on appropriate data transformation. I also describe the addition of supervised sva (ssva) for using control probes to identify the part of the genomic data only affected by artifacts. I present a comparison between these versions of sva and other methods for batch effect estimation on simulated data, real count-based data and FPKM-based data. These updates are available through the sva Bioconductor package and I have made fully reproducible analysis using these methods available from: https://github.com/jtleek/svaseq .

Description: High-throughput techniques have considerably increased the potential of comparative genomics whilst simultaneously posing many new challenges. One of those challenges involves efficiently mining the large amount of data produced and exploring the landscape of both conserved and idiosyncratic genomic regions across multiple genomes. Domains of application of these analyses are diverse: identification of evolutionary events, inference of gene functions, detection of niche-specific genes or phylogenetic profiling. Insyght is a comparative genomic visualization tool that combines three complementary displays: (i) a table for thoroughly browsing amongst homologues, (ii) a comparator of orthologue functional annotations and (iii) a genomic organization view designed to improve the legibility of rearrangements and distinctive loci. The latter display combines symbolic and proportional graphical paradigms. Synchronized navigation across multiple species and interoperability between the views are core features of Insyght. A gene filter mechanism is provided that helps the user to build a biologically relevant gene set according to multiple criteria such as presence/absence of homologues and/or various annotations. We illustrate the use of Insyght with scenarios. Currently, only Bacteria and Archaea are supported. A public instance is available at http://genome.jouy.inra.fr/Insyght . The tool is freely downloadable for private data set analysis.

Description: BRCA1 is involved in many disparate cellular functions, including DNA damage repair, cell-cycle checkpoint activation, gene transcriptional regulation, DNA replication, centrosome function and others. The majority of evidence strongly favors the maintenance of genomic integrity as a principal tumor suppressor activity of BRCA1. At the same time some functional aspects of BRCA1 are not fully understood. Here, a HAC (human artificial chromosome) module with a regulated centromere was constructed for delivery and expression of the 90 kb genomic copy of the BRCA1 gene into BRCA1-deficient human cells. A battery of functional tests was carried out to demonstrate functionality of the exogenous BRCA1. In separate experiments, we investigated the role of BRCA1 in maintenance of heterochromatin integrity within a human functional kinetochore. We demonstrated that BRCA1 deficiency results in a specific activation of transcription of higher-order alpha-satellite repeats (HORs) assembled into heterochromatin domains flanking the kinetochore. At the same time no detectable elevation of transcription was observed within HORs assembled into centrochromatin domains. Thus, we demonstrated a link between BRCA1 deficiency and kinetochore dysfunction and extended previous observations that BRCA1 is required to silence transcription in heterochromatin in specific genomic loci. This supports the hypothesis that epigenetic alterations of the kinetochore initiated in the absence of BRCA1 may contribute to cellular transformation.

Description: Nucleic acids have become a powerful tool in nanotechnology because of their controllable diverse conformational transitions and adaptable higher-order nanostructure. Using single-stranded DNA probes as the pore-caps for various target recognition, here we present an ultrasensitive universal electrochemical detection system based on graphene and mesoporous silica, and achieve sensitivity with all of the major classes of analytes and simultaneously realize DNA logic gate operations. The concept is based on the locking of the pores and preventing the signal-reporter molecules from escape by target-induced the conformational change of the tailored DNA caps. The coupling of ‘waking up’ gatekeeper with highly specific biochemical recognition is an innovative strategy for the detection of various targets, able to compete with classical methods which need expensive instrumentation and sophisticated experimental operations. The present study has introduced a new electrochemical signal amplification concept and also adds a new dimension to the function of graphene-mesoporous materials hybrids as multifunctional nanoscale logic devices. More importantly, the development of this approach would spur further advances in important areas, such as point-of-care diagnostics or detection of specific biological contaminations, and hold promise for use in field analysis.

Description: Many critical processes in the cell involve direct binding between RNAs and proteins, making it imperative to fully understand the physicochemical principles behind such interactions at the atomistic level. Here, we use molecular dynamics simulations and 15 μs of sampling to study the behavior of amino acids and amino acid sidechain analogs in high-concentration aqueous solutions of standard RNA nucleobases. Structural and energetic analysis of simulated systems allows us to derive interaction propensity scales for different amino acid/nucleobase combinations. The derived scales closely match and greatly extend the available experimental data, providing a comprehensive foundation for studying RNA–protein interactions in different contexts. By using these scales, we demonstrate a statistically significant connection between nucleobase composition of human mRNA coding sequences and nucleobase interaction propensities of their cognate protein sequences. For example, pyrimidine density profiles of mRNAs match uracil-propensity profiles of their cognate proteins with a median Pearson correlation coefficient of R = –0.70. Our results provide support for the recently proposed hypotheses that mRNAs and their cognate proteins may be physicochemically complementary to each other and bind, especially if unstructured, with the complementarity level being negatively influenced by mRNA adenine content. Finally, we utilize the derived scales to refine the complementarity hypothesis and closely examine its physicochemical underpinnings.

Description: The 54 promoters are unique in prokaryotic genome and responsible for transcripting carbon and nitrogen-related genes. With the avalanche of genome sequences generated in the postgenomic age, it is highly desired to develop automated methods for rapidly and effectively identifying the 54 promoters. Here, a predictor called ‘ iPro54-PseKNC ’ was developed. In the predictor, the samples of DNA sequences were formulated by a novel feature vector called ‘pseudo k -tuple nucleotide composition’, which was further optimized by the incremental feature selection procedure. The performance of iPro54-PseKNC was examined by the rigorous jackknife cross-validation tests on a stringent benchmark data set. As a user-friendly web-server, iPro54-PseKNC is freely accessible at http://lin.uestc.edu.cn/server/iPro54-PseKNC . For the convenience of the vast majority of experimental scientists, a step-by-step protocol guide was provided on how to use the web-server to get the desired results without the need to follow the complicated mathematics that were presented in this paper just for its integrity. Meanwhile, we also discovered through an in-depth statistical analysis that the distribution of distances between the transcription start sites and the translation initiation sites were governed by the gamma distribution, which may provide a fundamental physical principle for studying the 54 promoters.

Description: Telomeres at chromosome ends are normally masked from proteins that signal and repair DNA double strand breaks (DSBs). Bulky DNA lesions can cause DSBs if they block DNA replication, unless they are bypassed by translesion (TLS) DNA polymerases. Here, we investigated roles for TLS polymerase , (pol) in preserving telomeres following acute physical UVC exposure and chronic chemical Cr(VI) exposure, which both induce blocking lesions. We report that pol protects against cytotoxicity and replication stress caused by Cr(VI), similar to results with ultraviolet C light (UVC). Both exposures induce ataxia telangiectasia and Rad3-related (ATR) kinase and pol accumulation into nuclear foci and localization to individual telomeres, consistent with replication fork stalling at DNA lesions. Pol-deficient cells exhibited greater numbers of telomeres that co-localized with DSB response proteins after exposures. Furthermore, the genotoxic exposures induced telomere aberrations associated with failures in telomere replication that were suppressed by pol. We propose that pol's ability to bypass bulky DNA lesions at telomeres is critical for proper telomere replication following genotoxic exposures.

Description: The DNA damage response triggers cell-cycle checkpoints, DNA repair and apoptosis using multiple post-translational modifications as molecular switches. However, how ubiquitination regulates ATR signaling in response to replication stress and single-strand break is still unclear. Here, we identified the deubiquitination enzyme (DUB) USP20 as a pivotal regulator of ATR-related DDR pathway. Through screening a panel of DUBs, we identified USP20 as critical for replication stress response. USP20 is phosphorylated by ATR, resulting in disassociation of the E3 ubiquitin ligase HERC2 from USP20 and USP20 stabilization. USP20 in turn deubiquitinates and stabilizes Claspin and enhances the activation of ATR-Chk1 signaling. These findings reveal USP20 to be a novel regulator of ATR-dependent DNA damage signaling.

Description: Sunlight-induced C to T mutation hotspots in skin cancers occur primarily at methylated CpG sites that coincide with sites of UV-induced cyclobutane pyrimidine dimer (CPD) formation. The C or 5-methyl-C in CPDs are not stable and deaminate to U and T, respectively, which leads to the insertion of A by DNA polymerase and defines a probable mechanism for the origin of UV-induced C to T mutations. We have now determined the photoproduct formation and deamination rates for 10 consecutive T= m CG CPDs over a full helical turn at the dyad axis of a nucleosome and find that whereas photoproduct formation and deamination is greatly inhibited for the CPDs closest to the histone surface, it is greatly enhanced for the outermost CPDs. Replacing the G in a T= m CG CPD with A greatly decreased the deamination rate. These results show that rotational position and flanking sequence in a nucleosome can significantly and synergistically modulate CPD formation and deamination that contribute to C to T mutations associated with skin cancer induction and may have influenced the evolution of the human genome.

Description: We present a discriminative learning method for pattern discovery of binding sites in nucleic acid sequences based on hidden Markov models. Sets of positive and negative example sequences are mined for sequence motifs whose occurrence frequency varies between the sets. The method offers several objective functions, but we concentrate on mutual information of condition and motif occurrence. We perform a systematic comparison of our method and numerous published motif-finding tools. Our method achieves the highest motif discovery performance, while being faster than most published methods. We present case studies of data from various technologies, including ChIP-Seq, RIP-Chip and PAR-CLIP, of embryonic stem cell transcription factors and of RNA-binding proteins, demonstrating practicality and utility of the method. For the alternative splicing factor RBM10, our analysis finds motifs known to be splicing-relevant. The motif discovery method is implemented in the free software package Discrover. It is applicable to genome- and transcriptome-scale data, makes use of available repeat experiments and aside from binary contrasts also more complex data configurations can be utilized.

Description: In Escherichia coli , the ATP-bound form of DnaA (ATP–DnaA) promotes replication initiation. During replication, the bound ATP is hydrolyzed to ADP to yield the ADP-bound form (ADP–DnaA), which is inactive for initiation. The chromosomal site DARS2 facilitates the regeneration of ATP–DnaA by catalyzing nucleotide exchange between free ATP and ADP bound to DnaA. However, the regulatory mechanisms governing this exchange reaction are unclear. Here, using in vitro reconstituted experiments, we show that two nucleoid-associated proteins, IHF and Fis, bind site-specifically to DARS2 to activate coordinately the exchange reaction. The regenerated ATP–DnaA was fully active in replication initiation and underwent DnaA–ATP hydrolysis. ADP–DnaA formed heteromultimeric complexes with IHF and Fis on DARS2 , and underwent nucleotide dissociation more efficiently than ATP–DnaA. Consistently, mutant analyses demonstrated that specific binding of IHF and Fis to DARS2 stimulates the formation of ATP–DnaA production, thereby promoting timely initiation. Moreover, we show that IHF– DARS2 binding is temporally regulated during the cell cycle, whereas Fis only binds to DARS2 in exponentially growing cells. These results elucidate the regulation of ATP–DnaA and replication initiation in coordination with the cell cycle and growth phase.

Description: The virF gene of Shigella , responsible for triggering the virulence cascade in this pathogenic bacterium, is transcriptionally repressed by the nucleoid-associated protein H-NS. The primary binding sites of H-NS within the promoter region of virF have been detected here by footprinting experiments in the presence of H-NS or its monomeric DNA-binding domain (H-NS ctd ), which displays the same specificity as intact H-NS. Of the 14 short DNA fragments identified, 10 overlap sequences similar to the H-NS binding motif. The ‘fast’, ‘intermediate’ and ‘slow’ H-NS binding events leading to the formation of the nucleoprotein complex responsible for transcription repression have been determined by time-resolved hydroxyl radical footprinting experiments in the presence of full-length H-NS. We demonstrate that this process is completed in ≤1 s and H-NS protections occur simultaneously on site I and site II of the virF promoter. Furthermore, all ‘fast’ protections have been identified in regions containing predicted H-NS binding motifs, in agreement with the hypothesis that H-NS nucleoprotein complex assembles from a few nucleation sites containing high-affinity binding sequences. Finally, data are presented showing that the 22-bp fragment corresponding to one of the HNS binding sites deviates from canonical B-DNA structure at three TpA steps.

Description: Human DBC1 (Deleted in Breast Cancer 1; KIAA1967; CCAR2) is a protein implicated in the regulation of apoptosis, transcription and histone modifications. Upon DNA damage, DBC1 is phosphorylated by ATM/ATR on Thr454 and this modification increases its inhibitory interaction with SIRT1, leading to p53 acetylation and p53-dependent apoptosis. Here, we report that the inhibition of SIRT1 by DBC1 in the DNA damage response (DDR) also depends on Chk2, the transducer kinase that is activated by ATM upon DNA lesions and contributes to the spreading of DNA damage signal. Indeed we found that inactivation of Chk2 reduces DBC1-SIRT1 binding, thus preventing p53 acetylation and DBC1-induced apoptosis. These events are mediated by Chk2 phosphorylation of the 11S proteasome activator REG on Ser247, which increases REG-DBC1 interaction and SIRT1 inhibition. Overall our results clarify the mechanisms underlying the DBC1-dependent SIRT1 inhibition and link, for the first time, Chk2 and REG to the ATM-DBC1-SIRT1 axis.

Description: The mechanisms of gene amplification in tumour cells are poorly understood and the relationship between extrachromosomal DNA molecules, named double minutes (dmins), and intrachromosomal homogeneously staining regions (hsr) is not documented at nucleotide resolution. Using fluorescent in situ hybridization and whole genome sequencing, we studied a xenografted human oligodendroglioma where the co-amplification of the EGFR and MYC loci was present in the form of dmins at early passages and of an hsr at later passages. The amplified regions underwent multiple rearrangements and deletions during the formation of the dmins and their transformation into hsr. In both forms of amplification, non-homologous end-joining and microhomology-mediated end-joining rather than replication repair mechanisms prevailed in fusions. Small fragments, some of a few tens of base pairs, were associated in contigs. They came from clusters of breakpoints localized hundreds of kilobases apart in the amplified regions. The characteristics of some pairs of junctions suggest that at least some fragments were not fused randomly but could result from the concomitant repair of neighbouring breakpoints during the interaction of remote DNA sequences. This characterization at nucleotide resolution of the transition between extra- and intrachromosome amplifications highlights a hitherto uncharacterized organization of the amplified regions suggesting the involvement of new mechanisms in their formation.

Description: DNA palindromes are hotspots for DNA double strand breaks, inverted duplications and intra-chromosomal translocations in a wide spectrum of organisms from bacteria to humans. These reactions are mediated by DNA secondary structures such as hairpins and cruciforms. In order to further investigate the pathways of formation and cleavage of these structures, we have compared the processing of a 460 base pair (bp) perfect palindrome in the Escherichia coli chromosome with the same construct interrupted by a 20 bp spacer to form a 480 bp interrupted palindrome. We show here that the perfect palindrome can form hairpin DNA structures on the templates of the leading- and lagging-strands in a replication-dependent reaction. In the presence of the hairpin endonuclease SbcCD, both copies of the replicated chromosome containing the perfect palindrome are cleaved, resulting in the formation of an unrepairable DNA double-strand break and cell death. This contrasts with the interrupted palindrome, which forms a hairpin on the lagging-strand template that is processed to form breaks, which can be repaired by homologous recombination.

Description: Nucleic acid-dependent ATPases are involved in nearly all aspects of DNA and RNA metabolism. Previous studies have described a number of mitochondrial helicases. However, double-stranded DNA-dependent ATPases, including translocases or enzymes remodeling DNA-protein complexes, have not been identified in mitochondria of the yeast Saccharomyces cerevisae . Here, we demonstrate that Irc3p is a mitochondrial double-stranded DNA-dependent ATPase of the Superfamily II. In contrast to the other mitochondrial Superfamily II enzymes Mss116p, Suv3p and Mrh4p, which are RNA helicases, Irc3p has a direct role in mitochondrial DNA (mtDNA) maintenance. Specific Irc3p-dependent mtDNA metabolic intermediates can be detected, including high levels of double-stranded DNA breaks that accumulate in irc3 mutants. irc3- related topology changes in rho- mtDNA can be reversed by the deletion of mitochondrial RNA polymerase RPO41, suggesting that Irc3p counterbalances adverse effects of transcription on mitochondrial genome stability.

Description: In Escherichia coli , an increase in the ATP bound form of the DnaA initiator protein results in hyperinitiation and inviability. Here, we show that such replication stress is tolerated during anaerobic growth. In hyperinitiating cells, a shift from anaerobic to aerobic growth resulted in appearance of fragmented chromosomes and a decrease in terminus concentration, leading to a dramatic increase in ori / ter ratio and cessation of cell growth. Aerobic viability was restored by reducing the level of reactive oxygen species (ROS) or by deleting mutM (Fpg glycosylase). The double-strand breaks observed in hyperinitiating cells therefore results from replication forks encountering single-stranded DNA lesions generated while removing oxidized bases, primarily 8-oxoG, from the DNA. We conclude that there is a delicate balance between chromosome replication and ROS inflicted DNA damage so the number of replication forks can only increase when ROS formation is reduced or when the pertinent repair is compromised.

Description: Transcription-coupled DNA repair (TCR) is a subpathway of nucleotide excision repair (NER) dedicated to rapid removal of DNA lesions in the transcribed strand of actively transcribed genes. The precise nature of the TCR signal and how the repair machinery gains access to lesions imbedded in stalled RNA polymerase II (RNAP II) complexes in eukaryotic cells are still enigmatic. RNAP II has an intrinsic capacity for transcription bypass of DNA lesions by incorporation or misincorporation of nucleotides across the lesions. It has been suggested that transcription bypass of lesions, which exposes the lesions, may be required for TCR. Here, we show that E1103G mutation of Rpb1, the largest subunit of RNAP II, which promotes transcription bypass of UV-induced cyclobutane pyrimidine dimers (CPDs), increases survival of UV irradiated yeast cells but attenuates TCR. The increased cell survival is independent of any NER subpathways. In contrast, G730D mutation of Rpb1, which impairs transcription bypass of CPDs, enhances TCR. Our results suggest that transcription bypass of lesions attenuates TCR but enhances cell tolerance to DNA lesions. Efficient stalling of RNAP II is essential for efficient TCR.

Description: The avian bacterial pathogen Mycoplasma gallisepticum is a good model for systems studies due to small genome and simplicity of regulatory pathways. In this study, we used RNA-Seq and MS-based proteomics to accurately map coding sequences, transcription start sites (TSSs) and transcript 3'-ends (T3Es). We used obtained data to investigate roles of TSSs and T3Es in stress-induced transcriptional responses. We identified 1061 TSSs at a false discovery rate of 10% and showed that almost all transcription in M. gallisepticum is initiated from classic TATAAT promoters surrounded by A/T-rich sequences. Our analysis revealed the pronounced operon structure complexity: on average, each coding operon has one internal TSS and T3Es in addition to the primary ones. Our transcriptomic approach based on the intervals between the two nearest transcript ends allowed us to identify two classes of T3Es: strong, unregulated, hairpin-containing T3Es and weak, heat shock-regulated, hairpinless T3Es. Comparing gene expression levels under different conditions revealed widespread and divergent transcription regulation in M. gallisepticum . Modeling suggested that the core promoter structure plays an important role in gene expression regulation. We have shown that the heat stress activation of cryptic promoters combined with the hairpinless T3Es suppression leads to widespread, seemingly non-functional transcription.

Description: Genome-wide assessment of protein–DNA interaction by chromatin immunoprecipitation followed by massive parallel sequencing (ChIP-seq) is a key technology for studying transcription factor (TF) localization and regulation of gene expression. Signal-to-noise-ratio and signal specificity in ChIP-seq studies depend on many variables, including antibody affinity and specificity. Thus far, efforts to improve antibody reagents for ChIP-seq experiments have focused mainly on generating higher quality antibodies. Here we introduce KOIN (knockout implemented normalization) as a novel strategy to increase signal specificity and reduce noise by using TF knockout mice as a critical control for ChIP-seq data experiments. Additionally, KOIN can identify ‘hyper ChIPable regions’ as another source of false-positive signals. As the use of the KOIN algorithm reduces false-positive results and thereby prevents misinterpretation of ChIP-seq data, it should be considered as the gold standard for future ChIP-seq analyses, particularly when developing ChIP-assays with novel antibody reagents.

Description: Intracellular reactive oxygen species as well as the exposure to harsh environmental conditions can cause, in the single chromosome of Bacillus subtilis spores, the formation of apurinic/apyrimidinic (AP) sites and strand breaks whose repair during outgrowth is crucial to guarantee cell viability. Whereas double-stranded breaks are mended by the nonhomologous end joining (NHEJ) system composed of an ATP-dependent DNA Ligase D (LigD) and the DNA-end-binding protein Ku, repair of AP sites would rely on an AP endonuclease or an AP-lyase, a polymerase and a ligase. Here we show that B. subtilis Ku ( Bsu Ku), along with its pivotal role in allowing joining of two broken ends by B. subtilis LigD ( Bsu LigD), is endowed with an AP/deoxyribose 5'-phosphate (5'-dRP)-lyase activity that can act on ssDNA, nicked molecules and DNA molecules without ends, suggesting a potential role in BER during spore outgrowth. Coordination with Bsu LigD makes possible the efficient joining of DNA ends with near terminal abasic sites. The role of this new enzymatic activity of Ku and its potential importance in the NHEJ pathway is discussed. The presence of an AP-lyase activity also in the homolog protein from the distantly related bacterium Pseudomonas aeruginosa allows us to expand our results to other bacterial Ku proteins.

Description: Huntington's disease is a fatal neurodegenerative disease caused by polyglutamine-expansion in huntingtin (HTT). Recent work showed that gene silencing approaches, including RNA interference (RNAi), improve disease readouts in mice. To advance RNAi to the clinic, we designed miHDS1, with robust knockdown of human HTT and minimized silencing of unintended transcripts. In Rhesus macaque , AAV delivery of miHDS1 to the putamen reduced HTT expression with no adverse effects on neurological status including fine and gross motor skills, no immune activation and no induction of neuropathology out to 6 weeks post injection. Others showed safety of a different HTT-targeting RNAi in monkeys for 6 months. Application of miHDS1 to Huntington's patients requires further safety testing in normal rodents, despite the fact that it was optimized for humans. To satisfy this regulatory requirement, we evaluated normal mice after AAV.miHDS1 injection. In contrast to monkeys, neurological deficits occurred acutely in mice brain and was attributed to off-target silencing through interactions of miHDS1 with the 3'UTR of other transcripts. While we resolved miHDS1 toxicity in mouse brain and maintained miHDS1-silencing efficacy, these studies highlight that optimizing nucleic acid-based medicines for safety in humans presents challenges for safety testing in rodents or other distantly related species.

Description: The RNA-binding protein L7Ae, known for its role in translation (as part of ribosomes) and RNA modification (as part of sn/oRNPs), has also been identified as a subunit of archaeal RNase P, a ribonucleoprotein complex that employs an RNA catalyst for the Mg 2+ -dependent 5' maturation of tRNAs. To better understand the assembly and catalysis of archaeal RNase P, we used a site-specific hydroxyl radical-mediated footprinting strategy to pinpoint the binding sites of Pyrococcus furiosus ( Pfu ) L7Ae on its cognate RNase P RNA (RPR). L7Ae derivatives with single-Cys substitutions at residues in the predicted RNA-binding interface (K42C/C71V, R46C/C71V, V95C/C71V) were modified with an iron complex of EDTA-2-aminoethyl 2-pyridyl disulfide. Upon addition of hydrogen peroxide and ascorbate, these L7Ae-tethered nucleases were expected to cleave the RPR at nucleotides proximal to the EDTA-Fe–modified residues. Indeed, footprinting experiments with an enzyme assembled with the Pfu RPR and five protein cofactors (POP5, RPP21, RPP29, RPP30 and L7Ae–EDTA-Fe) revealed specific RNA cleavages, localizing the binding sites of L7Ae to the RPR's catalytic and specificity domains. These results support the presence of two kink-turns, the structural motifs recognized by L7Ae, in distinct functional domains of the RPR and suggest testable mechanisms by which L7Ae contributes to RNase P catalysis.

Description: The intricate network of interactions observed in RNA three-dimensional structures is often described in terms of a multitude of geometrical properties, including helical parameters, base pairing/stacking, hydrogen bonding and backbone conformation. We show that a simple molecular representation consisting in one oriented bead per nucleotide can account for the fundamental structural properties of RNA. In this framework, canonical Watson-Crick, non-Watson-Crick base-pairing and base-stacking interactions can be unambiguously identified within a well-defined interaction shell. We validate this representation by performing two independent, complementary tests. First, we use it to construct a sequence-independent, knowledge-based scoring function for RNA structural prediction, which compares favorably to fully atomistic, state-of-the-art techniques. Second, we define a metric to measure deviation between RNA structures that directly reports on the differences in the base–base interaction network. The effectiveness of this metric is tested with respect to the ability to discriminate between structurally and kinetically distant RNA conformations, performing better compared to standard techniques. Taken together, our results suggest that this minimalist, nucleobase-centric representation captures the main interactions that are relevant for describing RNA structure and dynamics.

Description: To elucidate the molecular mechanism of the integration of long interspersed elements (LINEs), we characterized the 5' ends of more than 200 LINE de novo retrotransposition events into chicken DT40 or human HeLa cells. Human L1 inserts produced 15-bp target-site duplications (TSDs) and zebrafish ZfL2-1 inserts produced 5-bp TSDs in DT40 cells, suggesting that TSD length depends on the LINE species. Further analysis of 5' junctions revealed that the 5'-end-joining pathways of LINEs can be divided into two fundamental types—annealing or direct. We also found that the generation of 5' inversions depends on host and LINE species. These results led us to propose a new model for 5'-end joining, the type of which is determined by the extent of exposure of 3' overhangs generated after the second-strand cleavage and by the involvement of host factors.

Description: The boundaries of R-loops are well-documented at immunoglobulin heavy chain loci in mammalian B cells. Within primary B cells or B cell lines, the upstream boundaries of R-loops typically begin early in the repetitive portion of the switch regions. Most R-loops terminate within the switch repetitive zone, but the remainder can extend a few hundred base pairs further, where G-density on the non-template DNA strand gradually drops to the genome average. Whether the G-density determines how far the R-loops extend is an important question. We previously studied the role of G-clusters in initiating R-loop formation, but we did not examine the role of G-density in permitting the elongation of the R-loop, after it had initiated. Here, we vary the G-density of different portions of the switch region in a murine B cell line. We find that both class switch recombination (CSR) and R-loop formation decrease significantly when the overall G-density is reduced from 46% to 29%. Short 50 bp insertions with low G-density within switch regions do not appear to affect either CSR or R-loop elongation, whereas a longer (150 bp) insertion impairs both. These results demonstrate that G-density is an important determinant of the length over which mammalian genomic R-loops extend.

Description: At equilibrium, empty ribosomes freely transit between the rotated and un-rotated states. In the cell, the binding of two translation elongation factors to the same general region of the ribosome stabilizes one state over the other. These stabilized states are resolved by expenditure of energy in the form of GTP hydrolysis. A prior study employing mutants of a late assembling peripheral ribosomal protein suggested that ribosome rotational status determines its affinity for elongation factors, and hence translational fidelity and gene expression. Here, mutants of the early assembling integral ribosomal protein uL2 are used to test the generality of this hypothesis. rRNA structure probing analyses reveal that mutations in the uL2 B7b bridge region shift the equilibrium toward the rotated state, propagating rRNA structural changes to all of the functional centers of ribosome. Structural disequilibrium unbalances ribosome biochemically: rotated ribosomes favor binding of the eEF2 translocase and disfavor that of the elongation ternary complex. This manifests as specific translational fidelity defects, impacting the expression of genes involved in telomere maintenance. A model is presented describing how cyclic intersubunit rotation ensures the unidirectionality of translational elongation, and how perturbation of rotational equilibrium affects specific aspects of translational fidelity and cellular gene expression.

Description: Although trans -translation mediated by tmRNA-SmpB has long been known as the sole system to relieve bacterial stalled ribosomes, ArfA has recently been identified as an alternative factor for ribosome rescue in Escherichia coli . This process requires hydrolysis of nascent peptidyl-tRNA by RF2, which usually acts as a stop codon-specific peptide release factor. It poses a fascinating question of how ArfA and RF2 recognize and rescue the stalled ribosome. Here, we mapped the location of ArfA in the stalled ribosome by directed hydroxyl radical probing. It revealed an ArfA-binding site around the neck region of the 30S subunit in which the N- and C-terminal regions of ArfA are close to the decoding center and the mRNA entry channel, respectively. ArfA and RF2 sequentially enter the ribosome stalled in either the middle or 3' end of mRNA, whereas RF2 induces a productive conformational change of ArfA only when ribosome is stalled at the 3' end of mRNA. On the basis of these results, we propose that ArfA functions as the sensor to recognize the target ribosome after RF2 binding.

Description: Rad53 is a conserved protein kinase with a central role in DNA damage response and nucleotide metabolism. We observed that the expression of a dominant-lethal form of RAD53 leads to significant expression changes for at least 16 genes, including the RNR3 and the HUG1 genes, both of which are involved in the control of nucleotide metabolism. We established by multiple biophysical and biochemical approaches that Hug1 is an intrinsically disordered protein that directly binds to the small RNR subunit Rnr2. We characterized the surface of interaction involved in Hug1 binding to Rnr2, and we thus defined a new binding region to Rnr2. Moreover, we show that Hug1 is deleterious to cell growth in the context of reduced RNR activity. This inhibitory effect of Hug1 on RNR activity depends on the binding of Hug1 to Rnr2. We propose a model in which Hug1 modulates Rnr2–Rnr1 association by binding Rnr2. We show that Hug1 accumulates under various physiological conditions of high RNR induction. Hence, both the regulation and the mode of action of Hug1 are different from those of the small protein inhibitors Dif1 and Sml1, and Hug1 can be considered as a regulator for fine-tuning of RNR activity.

Description: LepA is a paralog of EF-G found in all bacteria. Deletion of lepA confers no obvious growth defect in Escherichia coli , and the physiological role of LepA remains unknown. Here, we identify nine strains ( dksA , molR1 , rsgA , tatB , tonB , tolR , ubiF , ubiG or ubiH ) in which lepA confers a synthetic growth phenotype. These strains are compromised for gene regulation, ribosome assembly, transport and/or respiration, indicating that LepA contributes to these functions in some way. We also use ribosome profiling to deduce the effects of LepA on translation. We find that loss of LepA alters the average ribosome density (ARD) for hundreds of mRNA coding regions in the cell, substantially reducing ARD in many cases. By contrast, only subtle and codon-specific changes in ribosome distribution along mRNA are seen. These data suggest that LepA contributes mainly to the initiation phase of translation. Consistent with this interpretation, the effect of LepA on ARD is related to the sequence of the Shine–Dalgarno region. Global perturbation of gene expression in the lepA mutant likely explains most of its phenotypes.

Description: Many ribosome-interacting GTPases, with proposed functions in ribosome biogenesis, are also implicated in the cellular regulatory coupling between ribosome assembly process and various growth control pathways. EngA is an essential GTPase in bacteria, and intriguingly, it contains two consecutive GTPase domains (GD), being one-of-a-kind among all known GTPases. EngA is required for the 50S subunit maturation. However, its molecular role remains elusive. Here, we present the structure of EngA bound to the 50S subunit. Our data show that EngA binds to the peptidyl transferase center (PTC) and induces dramatic conformational changes on the 50S subunit, which virtually returns the 50S subunit to a state similar to that of the late-stage 50S assembly intermediates. Very interestingly, our data show that the two GDs exhibit a pseudo-two-fold symmetry in the 50S-bound conformation. Our results indicate that EngA recognizes certain forms of the 50S assembly intermediates, and likely facilitates the conformational maturation of the PTC of the 23S rRNA in a direct manner. Furthermore, in a broad context, our data also suggest that EngA might be a sensor of the cellular GTP/GDP ratio, endowed with multiple conformational states, in response to fluctuations in cellular nucleotide pool, to facilitate and regulate ribosome assembly.

Description: The occurrence of a G-triplex folding intermediate of thrombin binding aptamer (TBA) has been recently predicted by metadynamics calculations, and experimentally supported by Nuclear Magnetic Resonance (NMR), Circular Dichroism (CD) and Differential Scanning Calorimetry (DSC) data collected on a 3' end TBA-truncated 11-mer oligonucleotide (11-mer-3'-t-TBA). Here we present the solution structure of 11-mer-3'-t-TBA in the presence of potassium ions. This structure is the first experimental example of a G-triplex folding, where a network of Hoogsteen-like hydrogen bonds stabilizes six guanines to form two G:G:G triad planes. The G-triplex folding of 11-mer-3'-t-TBA is stabilized by the potassium ion and destabilized by increasing the temperature. The superimposition of the experimental structure with that predicted by metadynamics shows a great similarity, with only significant differences involving two loops. These new structural data show that 11-mer-3'-t-TBA assumes a G-triplex DNA conformation as its stable form, reinforcing the idea that G-triplex folding intermediates may occur in vivo in human guanine-rich sequences. NMR and CD screening of eight different constructs obtained by removing from one to four bases at either the 3' and the 5' ends show that only the 11-mer-3'-t-TBA yields a relatively stable G-triplex.

Description: The accurate construction and interpretation of gene association networks (GANs) is challenging, but crucial, to the understanding of gene function, interaction and cellular behavior at the genome level. Most current state-of-the-art computational methods for genome-wide GAN reconstruction require high-performance computational resources. However, even high-performance computing cannot fully address the complexity involved with constructing GANs from very large-scale expression profile datasets, especially for the organisms with medium to large size of genomes, such as those of most plant species. Here, we present a new approach, GPLEXUS ( http://plantgrn.noble.org/GPLEXUS/ ), which integrates a series of novel algorithms in a parallel-computing environment to construct and analyze genome-wide GANs. GPLEXUS adopts an ultra-fast estimation for pairwise mutual information computing that is similar in accuracy and sensitivity to the Algorithm for the Reconstruction of Accurate Cellular Networks (ARACNE) method and runs ~1000 times faster. GPLEXUS integrates Markov Clustering Algorithm to effectively identify functional subnetworks. Furthermore, GPLEXUS includes a novel ‘condition-removing’ method to identify the major experimental conditions in which each subnetwork operates from very large-scale gene expression datasets across several experimental conditions, which allows users to annotate the various subnetworks with experiment-specific conditions. We demonstrate GPLEXUS’s capabilities by construing global GANs and analyzing subnetworks related to defense against biotic and abiotic stress, cell cycle growth and division in Arabidopsis thaliana .

Description: To reveal the full potential of human pluripotent stem cells, new methods for rapid, site-specific genomic engineering are needed. Here, we describe a system for precise genetic modification of human embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs). We identified a novel human locus, H11 , located in a safe, intergenic, transcriptionally active region of chromosome 22, as the recipient site, to provide robust, ubiquitous expression of inserted genes. Recipient cell lines were established by site-specific placement of a ‘landing pad’ cassette carrying attP sites for phiC31 and Bxb1 integrases at the H11 locus by spontaneous or TALEN-assisted homologous recombination. Dual integrase cassette exchange (DICE) mediated by phiC31 and Bxb1 integrases was used to insert genes of interest flanked by phiC31 and Bxb1 attB sites at the H11 locus, replacing the landing pad. This system provided complete control over content, direction and copy number of inserted genes, with a specificity of 100%. A series of genes, including mCherry and various combinations of the neural transcription factors LMX1a, FOXA2 and OTX2, were inserted in recipient cell lines derived from H9 ESC, as well as iPSC lines derived from a Parkinson’s disease patient and a normal sibling control. The DICE system offers rapid, efficient and precise gene insertion in ESC and iPSC and is particularly well suited for repeated modifications of the same locus.

Description: Distant genomic elements were found to interact within the folded eukaryotic genome. However, the used experimental approach (chromosome conformation capture, 3C) enables neither determination of the percentage of cells in which the interactions occur nor demonstration of simultaneous interaction of 〉2 genomic elements. Each of the above can be done using in-gel replication of interacting DNA segments, the technique reported here. Chromatin fragments released from formaldehyde–cross-linked cells by sodium dodecyl sulfate extraction and sonication are distributed in a polyacrylamide gel layer followed by amplification of selected test regions directly in the gel by multiplex polymerase chain reaction. The fragments that have been cross-linked and separate fragments give rise to multi- and monocomponent molecular colonies, respectively, which can be distinguished and counted. Using in-gel replication of interacting DNA segments, we demonstrate that in the material from mouse erythroid cells, the majority of fragments containing the promoters of active β-globin genes and their remote enhancers do not form complexes stable enough to survive sodium dodecyl sulfate extraction and sonication. This indicates that either these elements do not interact directly in the majority of cells at a given time moment, or the formed DNA–protein complex cannot be stabilized by formaldehyde cross-linking.

Description: The enzyme predominantly used for in vitro run-off RNA synthesis is bacteriophage T7 RNA polymerase. T7 RNA polymerase synthesizes, in addition to run-off products of precise length, transcripts with an additional non-base-paired nucleotide at the 3'-terminus (N + 1 product). This contaminating product is extremely difficult to remove. We recently characterized the single-subunit RNA polymerase from marine cyanophage Syn5 and identified its promoter sequence. This marine enzyme catalyses RNA synthesis over a wider range of temperature and salinity than does T7 RNA polymerase. Its processivity is 〉30 000 nt without significant intermediate products. The requirement for the initiating nucleotide at the promoter is less stringent for Syn5 RNA polymerase as compared to T7 RNA polymerase. A major difference is the precise run-off transcripts with homogeneous 3'-termini synthesized by Syn5 RNA polymerase. Therefore, the enzyme is advantageous for the production of RNAs that require precise 3'-termini, such as tRNAs and RNA fragments that are used for subsequent assembly.

Description: Genetic disorders can be detected by prenatal diagnosis using Chorionic Villus Sampling, but the 1:100 chance to result in miscarriage restricts the use to fetuses that are suspected to have an aberration. Detection of trisomy 21 cases noninvasively is now possible owing to the upswing of next-generation sequencing (NGS) because a small percentage of fetal DNA is present in maternal plasma. However, detecting other trisomies and smaller aberrations can only be realized using high-coverage NGS, making it too expensive for routine practice. We present a method, WISECONDOR (WIthin-SamplE COpy Number aberration DetectOR), which detects small aberrations using low-coverage NGS. The increased detection resolution was achieved by comparing read counts within the tested sample of each genomic region with regions on other chromosomes that behave similarly in control samples. This within-sample comparison avoids the need to re-sequence control samples. WISECONDOR correctly identified all T13, T18 and T21 cases while coverages were as low as 0.15–1.66. No false positives were identified. Moreover, WISECONDOR also identified smaller aberrations, down to 20 Mb, such as del(13)(q12.3q14.3), +i(12)(p10) and i(18)(q10). This shows that prevalent fetal copy number aberrations can be detected accurately and affordably by shallow sequencing maternal plasma. WISECONDOR is available at bioinformatics.tudelft.nl/wisecondor.

Description: In eukaryotic cells, the shortening and removal of the poly(A) tail of cytoplasmic mRNA by deadenylase enzymes is a critical step in post-transcriptional gene regulation. The ribonuclease activity of deadenylase enzymes is attributed to either a DEDD (Asp-Glu-Asp-Asp) or an endonuclease–exonuclease–phosphatase domain. Both domains require the presence of two Mg 2+ ions in the active site. To facilitate the biochemical analysis of deadenylase enzymes, we have developed a fluorescence-based deadenylase assay. The assay is based on end-point measurement, suitable for quantitative analysis and can be adapted for 96- and 384-well microplate formats. We demonstrate the utility of the assay by screening a chemical compound library, resulting in the identification of non-nucleoside inhibitors of the Caf1/CNOT7 enzyme, a catalytic subunit of the Ccr4–Not deadenylase complex. These compounds may be useful tools for the biochemical analysis of the Caf1/CNOT7 deadenylase subunit of the Ccr4–Not complex and indicate the feasibility of developing selective inhibitors of deadenylase enzymes using the fluorescence-based assay.

Description: Deciphering the causal networks of gene interactions is critical for identifying disease pathways and disease-causing genes. We introduce a method to reconstruct causal networks based on exploring phenotype-specific modules in the human interactome and including the expression quantitative trait loci (eQTLs) that underlie the joint expression variation of each module. Closely associated eQTLs help anchor the orientation of the network. To overcome the inherent computational complexity of causal network reconstruction, we first deduce the local causality of individual subnetworks using the selected eQTLs and module transcripts. These subnetworks are then integrated to infer a global causal network using a random-field ranking method, which was motivated by animal sociology. We demonstrate how effectively the inferred causality restores the regulatory structure of the networks that mediate lymph node metastasis in oral cancer. Network rewiring clearly characterizes the dynamic regulatory systems of distinct disease states. This study is the first to associate an RXRB -causal network with increased risks of nodal metastasis, tumor relapse, distant metastases and poor survival for oral cancer. Thus, identifying crucial upstream drivers of a signal cascade can facilitate the discovery of potential biomarkers and effective therapeutic targets.

Description: Recent sequencing technologies that allow massive parallel production of short reads are the method of choice for transcriptome analysis. Particularly, digital gene expression (DGE) technologies produce a large dynamic range of expression data by generating short tag signatures for each cell transcript. These tags can be mapped back to a reference genome to identify new transcribed regions that can be further covered by RNA-sequencing (RNA-Seq) reads. Here, we applied an integrated bioinformatics approach that combines DGE tags, RNA-Seq, tiling array expression data and species-comparison to explore new transcriptional regions and their specific biological features, particularly tissue expression or conservation. We analysed tags from a large DGE data set (designated as ‘TranscriRef’). We then annotated 750 000 tags that were uniquely mapped to the human genome according to Ensembl. We retained transcripts originating from both DNA strands and categorized tags corresponding to protein-coding genes, antisense, intronic- or intergenic-transcribed regions and computed their overlap with annotated non-coding transcripts. Using this bioinformatics approach, we identified ~34 000 novel transcribed regions located outside the boundaries of known protein-coding genes. As demonstrated using sequencing data from human pluripotent stem cells for biological validation, the method could be easily applied for the selection of tissue-specific candidate transcripts. DigitagCT is available at http://cractools.gforge.inria.fr/softwares/digitagct .

Description: Combinatorial interplay among transcription factors (TFs) is an important mechanism by which transcriptional regulatory specificity is achieved. However, despite the increasing number of TFs for which either binding specificities or genome-wide occupancy data are known, knowledge about cooperativity between TFs remains limited. To address this, we developed a computational framework for predicting genome-wide co-binding between TFs (CCAT, C ombinatorial C ode A nalysis T ool), and applied it to Drosophila melanogaster to uncover cooperativity among TFs during embryo development. Using publicly available TF binding specificity data and DNaseI chromatin accessibility data, we first predicted genome-wide binding sites for 324 TFs across five stages of D. melanogaster embryo development. We then applied CCAT in each of these developmental stages, and identified from 19 to 58 pairs of TFs in each stage whose predicted binding sites are significantly co-localized. We found that nearby binding sites for pairs of TFs predicted to cooperate were enriched in regions bound in relevant ChIP experiments, and were more evolutionarily conserved than other pairs. Further, we found that TFs tend to be co-localized with other TFs in a dynamic manner across developmental stages. All generated data as well as source code for our front-to-end pipeline are available at http://cat.princeton.edu .

Description: Recombineering, which is the use of homologous recombination for DNA engineering in Escherichia coli , usually uses antibiotic selection to identify the intended recombinant. When combined in a second step with counterselection using a small molecule toxin, seamless products can be obtained. Here, we report the advantages of a genetic strategy using CcdB as the counterselectable agent. Expression of CcdB is toxic to E. coli in the absence of the CcdA antidote so counterselection is initiated by the removal of CcdA expression. CcdB counterselection is robust and does not require titrations or experiment-to-experiment optimization. Because counterselection strategies necessarily differ according to the copy number of the target, we describe two variations. For multi-copy targets, we use two E. coli hosts so that counterselection is exerted by the transformation step that is needed to separate the recombined and unrecombined plasmids. For single copy targets, we put the ccdA gene onto the temperature-sensitive pSC101 Red expression plasmid so that counterselection is exerted by the standard temperature shift to remove the expression plasmid. To reduce unwanted intramolecular recombination, we also combined CcdB counterselection with Redα omission. These options improve the use of counterselection in recombineering with BACs, plasmids and the E. coli chromosome.

Description: As biotechnology advances rapidly, a tremendous amount of cancer genetic data has become available, providing an unprecedented opportunity for understanding the genetic mechanisms of cancer. To understand the effects of duplications and deletions on cancer progression, two genomes (normal and tumor) were sequenced from each of five stomach cancer patients in different stages (I, II, III and IV). We developed a phylogenetic model for analyzing stomach cancer data. The model assumes that duplication and deletion occur in accordance with a continuous time Markov Chain along the branches of a phylogenetic tree attached with five extended branches leading to the tumor genomes. Moreover, coalescence times of the phylogenetic tree follow a coalescence process. The simulation study suggests that the maximum likelihood approach can accurately estimate parameters in the phylogenetic model. The phylogenetic model was applied to the stomach cancer data. We found that the expected number of changes (duplication and deletion) per gene for the tumor genomes is significantly higher than that for the normal genomes. The goodness-of-fit test suggests that the phylogenetic model with constant duplication and deletion rates can adequately fit the duplication data for the normal genomes. The analysis found nine duplicated genes that are significantly associated with stomach cancer.

Description: Ten-eleven translocation (TET) family enzymes convert 5-methylcytosine to 5-hydroxylmethylcytosine. However, the molecular mechanism that regulates this biological process is not clear. Here, we show the evidence that PGC7 (also known as Dppa3 or Stella) interacts with TET2 and TET3 both in vitro and in vivo to suppress the enzymatic activity of TET2 and TET3. Moreover, lacking PGC7 induces the loss of DNA methylation at imprinting loci. Genome-wide analysis of PGC7 reveals a consensus DNA motif that is recognized by PGC7. The CpG islands surrounding the PGC7-binding motifs are hypermethylated. Taken together, our study demonstrates a molecular mechanism by which PGC7 protects DNA methylation from TET family enzyme-dependent oxidation.

Description: The anti-silencing function protein 1 (Asf1) is a chaperone that forms a complex with histones H3 and H4 facilitating dimer deposition and removal from chromatin. Most eukaryotes possess two different Asf1 chaperones but their specific functions are still unknown. Trypanosomes, a group of early-diverged eukaryotes, also have two, but more divergent Asf1 paralogs than Asf1 of higher eukaryotes. To unravel possible different functions, we characterized the two Asf1 proteins in Trypanosoma brucei . Asf1A is mainly localized in the cytosol but translocates to the nucleus in S phase. In contrast, Asf1B is predominantly localized in the nucleus, as described for other organisms. Cytosolic Asf1 knockdown results in accumulation of cells in early S phase of the cell cycle, whereas nuclear Asf1 knockdown arrests cells in S/G2 phase. Overexpression of cytosolic Asf1 increases the levels of histone H3 and H4 acetylation. In contrast to cytosolic Asf1, overexpression of nuclear Asf1 causes less pronounced growth defects in parasites exposed to genotoxic agents, prompting a function in chromatin remodeling in response to DNA damage. Only the cytosolic Asf1 interacts with recombinant H3/H4 dimers in vitro . These findings denote the early appearance in evolution of distinguishable functions for the two Asf1 chaperons in trypanosomes.

Description: DNA, RNA and proteins are major biological macromolecules that coevolve and adapt to environments as components of one highly interconnected system. We explore here sequence/structure determinants of mechanisms of adaptation of these molecules, links between them, and results of their mutual evolution. We complemented statistical analysis of genomic and proteomic sequences with folding simulations of RNA molecules, unraveling causal relations between compositional and sequence biases reflecting molecular adaptation on DNA, RNA and protein levels. We found many compositional peculiarities related to environmental adaptation and the life style. Specifically, thermal adaptation of protein-coding sequences in Archaea is characterized by a stronger codon bias than in Bacteria. Guanine and cytosine load in the third codon position is important for supporting the aerobic life style, and it is highly pronounced in Bacteria. The third codon position also provides a tradeoff between arginine and lysine, which are favorable for thermal adaptation and aerobicity, respectively. Dinucleotide composition provides stability of nucleic acids via strong base-stacking in ApG dinucleotides. In relation to coevolution of nucleic acids and proteins, thermostability-related demands on the amino acid composition affect the nucleotide content in the second codon position in Archaea.

Description: Alternative splicing (AS), in higher eukaryotes, is one of the mechanisms of post-transcriptional regulation that generate multiple transcripts from the same gene. One particular mode of AS is the skipping event where an exon may be alternatively excluded or constitutively included in the resulting mature mRNA. Both transcript isoforms from this skipping event site, i.e . in which the exon is either included (inclusion isoform) or excluded (skipping isoform), are typically present in one cell, and maintain a subtle balance that is vital to cellular function and dynamics. However, how the prevailing conditions dictate which isoform is expressed and what biological factors might influence the regulation of this process remain areas requiring further exploration. In this study, we have developed a novel computational method, graph-based exon-skipping scanner (GESS), for de novo detection of skipping event sites from raw RNA-seq reads without prior knowledge of gene annotations, as well as for determining the dominant isoform generated from such sites. We have applied our method to publicly available RNA-seq data in GM12878 and K562 cells from the ENCODE consortium and experimentally validated several skipping site predictions by RT-PCR. Furthermore, we integrated other sequencing-based genomic data to investigate the impact of splicing activities, transcription factors (TFs) and epigenetic histone modifications on splicing outcomes. Our computational analysis found that splice sites within the skipping-isoform-dominated group (SIDG) tended to exhibit weaker MaxEntScan-calculated splice site strength around middle, ‘skipping’, exons compared to those in the inclusion-isoform-dominated group (IIDG). We further showed the positional preference pattern of splicing factors, characterized by enrichment in the intronic splice sites immediately bordering middle exons. Finally, our analysis suggested that different epigenetic factors may introduce a variable obstacle in the process of exon–intron boundary establishment leading to skipping events.

Description: Recent advances in technology have led to a dramatic increase in the number of available transcription factor ChIP-seq and ChIP-chip data sets. Understanding the motif content of these data sets is an important step in understanding the underlying mechanisms of regulation. Here we provide a systematic motif analysis for 427 human ChIP-seq data sets using motifs curated from the literature and also discovered de novo using five established motif discovery tools. We use a systematic pipeline for calculating motif enrichment in each data set, providing a principled way for choosing between motif variants found in the literature and for flagging potentially problematic data sets. Our analysis confirms the known specificity of 41 of the 56 analyzed factor groups and reveals motifs of potential cofactors. We also use cell type-specific binding to find factors active in specific conditions. The resource we provide is accessible both for browsing a small number of factors and for performing large-scale systematic analyses. We provide motif matrices, instances and enrichments in each of the ENCODE data sets. The motifs discovered here have been used in parallel studies to validate the specificity of antibodies, understand cooperativity between data sets and measure the variation of motif binding across individuals and species.

Description: The identification of transcription factor binding motifs is important for the study of gene transcriptional regulation. The chromatin immunoprecipitation (ChIP), followed by massive parallel sequencing (ChIP-seq) experiments, provides an unprecedented opportunity to discover binding motifs. Computational methods have been developed to identify motifs from ChIP-seq data, while at the same time encountering several problems. For example, existing methods are often not scalable to the large number of sequences obtained from ChIP-seq peak regions. Some methods heavily rely on well-annotated motifs even though the number of known motifs is limited. To simplify the problem, de novo motif discovery methods often neglect underrepresented motifs in ChIP-seq peak regions. To address these issues, we developed a novel approach called SIOMICS to de novo discover motifs from ChIP-seq data. Tested on 13 ChIP-seq data sets, SIOMICS identified motifs of many known and new cofactors. Tested on 13 simulated random data sets, SIOMICS discovered no motif in any data set. Compared with two recently developed methods for motif discovery, SIOMICS shows advantages in terms of speed, the number of known cofactor motifs predicted in experimental data sets and the number of false motifs predicted in random data sets. The SIOMICS software is freely available at http://eecs.ucf.edu/~xiaoman/SIOMICS/SIOMICS.html .

Description: The regulatory networks of differentiation programs and the molecular mechanisms of lineage-specific gene regulation in mammalian embryos remain only partially defined. We document differential expression and temporal switching of BRG1-associated factor (BAF) subunits, core pluripotency factors and cardiac-specific genes during post-implantation development and subsequent early organogenesis. Using affinity purification of BRG1 ATPase coupled to mass spectrometry, we characterized the cardiac-enriched remodeling complexes present in E8.5 mouse embryos. The relative abundance and combinatorial assembly of the BAF subunits provides functional specificity to Switch/Sucrose NonFermentable (SWI/SNF) complexes resulting in a unique gene expression profile in the developing heart. Remarkably, the specific depletion of the BAF250a subunit demonstrated differential effects on cardiac-specific gene expression and resulted in arrhythmic contracting cardiomyocytes in vitro . Indeed, the BAF250a physically interacts and functionally cooperates with Nucleosome Remodeling and Histone Deacetylase (NURD) complex subunits to repressively regulate chromatin structure of the cardiac genes by switching open and poised chromatin marks associated with active and repressed gene expression. Finally, BAF250a expression modulates BRG1 occupancy at the loci of cardiac genes regulatory regions in P19 cell differentiation. These findings reveal specialized and novel cardiac-enriched SWI/SNF chromatin-remodeling complexes, which are required for heart formation and critical for cardiac gene expression regulation at the early stages of heart development.

Description: Understanding the structure of interphase chromosomes is essential to elucidate regulatory mechanisms of gene expression. During recent years, high-throughput DNA sequencing expanded the power of chromosome conformation capture (3C) methods that provide information about reciprocal spatial proximity of chromosomal loci. Since 2012, it is known that entire chromatin in interphase chromosomes is organized into regions with strongly increased frequency of internal contacts. These regions, with the average size of ~1 Mb, were named topological domains. More recent studies demonstrated presence of unconstrained supercoiling in interphase chromosomes. Using Brownian dynamics simulations, we show here that by including supercoiling into models of topological domains one can reproduce and thus provide possible explanations of several experimentally observed characteristics of interphase chromosomes, such as their complex contact maps.

Description: Tumor angiogenesis is mainly mediated by vascular endothelial growth factor (VEGF), a pro-angiogenic factor produced by cancer cells and active on the endothelium through the VEGF receptor 2 (VEGFR-2). Here we identify a G-rich sequence within the proximal promoter region of vegfr-2 , able to form an antiparallel G-quadruplex (G4) structure. This G4 structure can be efficiently stabilized by small molecules with the consequent inhibition of vegfr-2 expression. Functionally, the G4-mediated reduction of VEGFR-2 protein causes a switching off of signaling components that, converging on actin cytoskeleton, regulate the cellular events leading to endothelial cell proliferation, migration and differentiation. As a result of endothelial cell function impairment, angiogenic process is strongly inhibited by G4 ligands both in vitro and in viv o. Interestingly, the G4-mediated antiangiogenic effect seems to recapitulate that observed by using a specific interference RNA against vegfr-2 , and it is strongly antagonized by overexpressing the vegfr-2 gene. In conclusion, we describe the evidence for the existence of G4 in the promoter of vegfr-2 , whose expression and function can be markedly inhibited by G4 ligands, thereby revealing a new, and so far undescribed, way to block VEGFR-2 as target for anticancer therapy.