Despite the increasing knowledge about DNA methylation, the understanding of human epigenome evolution is in its infancy. Using whole genome bisulfite sequencing we identified hundreds of differentially methylated regions (DMRs) in humans compared to non-human primates and estimated that ∼25% of these regions were detectable throughout several human tissues. Human DMRs were enriched for specific histone modifications and the majority were located distal to transcription start sites, highlighting the importance of regions outside the direct regulatory context. We also found a significant excess of endogenous retrovirus elements in human-specific hypomethylated.We reported for the first time a close interplay between inter-species genetic and epigenetic variation in regions of incomplete lineage sorting, transcription factor binding sites and human differentially hypermethylated regions. Specifically, we observed an excess of human-specific substitutions in transcription factor binding sites located within human DMRs, suggesting that alteration of regulatory motifs underlies some human-specific methylation patterns. We also found that the acquisition of DNA hypermethylation in the human lineage is frequently coupled with a rapid evolution at nucleotide level in the neighborhood of these CpG sites. Taken together, our results reveal new insights into the mechanistic basis of human-specific DNA methylation patterns and the interpretation of inter-species non-coding variation.

Polyamines are ubiquitous cations that are involved in regulating fundamental cellular processes such as cell growth and proliferation; hence, their intracellular concentration is tightly regulated. Antizyme and antizyme inhibitor have a central role in maintaining cellular polyamine levels. Antizyme is unique in that it is expressed via a novel programmed ribosomal frameshifting mechanism. Conventional computational tools are unable to predict a programmed frameshift, resulting in misannotation of antizyme transcripts and proteins on transcript and genomic sequences. Correct annotation of a programmed frameshifting event requires manual evaluation. Our goal was to provide an accurately curated and annotated Reference Sequence (RefSeq) data set of antizyme transcript and protein records across a broad taxonomic scope that would serve as standards for accurate representation of these gene products. As antizyme and antizyme inhibitor proteins are functionally connected, we also curated antizyme inhibitor genes to more fully represent the elegant biology of polyamine regulation. Manual review of genes for three members of the antizyme family and two members of the antizyme inhibitor family in 91 vertebrate organisms resulted in a total of 461 curated RefSeq records.

During ribosomal translation, the two ribosomal subunits remain associated through intersubunit bridges, despite rapid large-scale intersubunit rotation. The absence of large barriers hindering rotation is a prerequisite for rapid rotation. Here, we investigate how such a flat free-energy landscape is achieved, in particular considering the large shifts the bridges undergo at the periphery. The dynamics and energetics of the intersubunit contact network are studied using molecular dynamics simulations of the prokaryotic ribosome in intermediate states of spontaneous translocation. Based on observed occupancies of intersubunit contacts, residues were grouped into clusters. In addition to the central contact clusters, peripheral clusters were found to maintain strong steady interactions by changing contacts in the course of rotation. The peripheral B1 bridges are stabilized by a changing contact pattern of charged residues that adapts to the rotational state. In contrast, steady strong interactions of the B4 bridge are ensured by the flexible helix H34 following the movement of protein S15. The tRNAs which span the subunits contribute to the intersubunit binding enthalpy to an almost constant degree, despite their different positions in the ribosome. These mechanisms keep the intersubunit interaction strong and steady during rotation, thereby preventing dissociation and enabling rapid rotation.

The heterogeneous nuclear ribonucleoprotein A1 (hnRNP-A1) has been implicated in telomere protection and telomerase activation. Recent evidence has further demonstrated that hnRNP-A1 plays a crucial role in maintaining newly replicated telomeric 3' overhangs and facilitating the switch from replication protein A (RPA) to protection of telomeres 1 (POT1). The role of hnRNP-A1 in telomere protection also involves DNA-dependent protein kinase catalytic subunit (DNA-PKcs), although the detailed regulation mechanism has not been clear. Here we report that hnRNP-A1 is phosphorylated by DNA-PKcs during the G2 and M phases and that DNA-PK-dependent hnRNP-A1 phosphorylation promotes the RPA-to-POT1 switch on telomeric single-stranded 3' overhangs. Consequently, in cells lacking hnRNP-A1 or DNA-PKcs-dependent hnRNP-A1 phosphorylation, impairment of the RPA-to-POT1 switch results in DNA damage response at telomeres during mitosis as well as induction of fragile telomeres. Taken together, our results indicate that DNA-PKcs-dependent hnRNP-A1 phosphorylation is critical for capping of the newly replicated telomeres and prevention of telomeric aberrations.

Tyrosine site-specific recombinases, which promote one class of biologically important phosphoryl transfer reactions in DNA, exemplify active site mechanisms for stabilizing the phosphate transition state. A highly conserved arginine duo (Arg-I; Arg-II) of the recombinase active site plays a crucial role in this function. Cre and Flp recombinase mutants lacking either arginine can be rescued by compensatory charge neutralization of the scissile phosphate via methylphosphonate (MeP) modification. The chemical chirality of MeP, in conjunction with mutant recombinases, reveals the stereochemical contributions of Arg-I and Arg-II. The SP preference of the native reaction is specified primarily by Arg-I. MeP reaction supported by Arg-II is nearly bias-free or RP-biased, depending on the Arg-I substituent. Positional conservation of the arginines does not translate into strict functional conservation. Charge reversal by glutamic acid substitution at Arg-I or Arg-II has opposite effects on Cre and Flp in MeP reactions. In Flp, the base immediately 5' to the scissile MeP strongly influences the choice between the catalytic tyrosine and water as the nucleophile for strand scission, thus between productive recombination and futile hydrolysis. The recombinase active site embodies the evolutionary optimization of interactions that not only favor the normal reaction but also proscribe antithetical side reactions.

The bacterial transcription termination factor Rho-a ring-shaped molecular motor displaying directional, ATP-dependent RNA helicase/translocase activity-is an interesting therapeutic target. Recently, Rho from Mycobacterium tuberculosis (MtbRho) has been proposed to operate by a mechanism uncoupled from molecular motor action, suggesting that the manner used by Rho to dissociate transcriptional complexes is not conserved throughout the bacterial kingdom. Here, however, we demonstrate that MtbRho is a bona fide molecular motor and directional helicase which requires a catalytic site competent for ATP hydrolysis to disrupt RNA duplexes or transcription elongation complexes. Moreover, we show that idiosyncratic features of the MtbRho enzyme are conferred by a large, hydrophilic insertion in its N-terminal 'RNA binding' domain and by a non-canonical R-loop residue in its C-terminal 'motor' domain. We also show that the 'motor' domain of MtbRho has a low apparent affinity for the Rho inhibitor bicyclomycin, thereby contributing to explain why M. tuberculosis is resistant to this drug. Overall, our findings support that, in spite of adjustments of the Rho motor to specific traits of its hosting bacterium, the basic principles of Rho action are conserved across species and could thus constitute pertinent screening criteria in high-throughput searches of new Rho inhibitors.

We present a new publicly accessible web-service, RiboSoft, which implements a comprehensive hammerhead ribozyme design procedure. It accepts as input a target sequence (and some design parameters) then generates a set of ranked hammerhead ribozymes, which target the input sequence. This paper describes the implemented procedure, which takes into consideration multiple objectives leading to a multi-objective ranking of the computer-generated ribozymes. Many ribozymes were assayed and validated, including four ribozymes targeting the transcript of a disease-causing gene (a mutant version of PABPN1). These four ribozymes were successfully tested in vitro and in vivo, for their ability to cleave the targeted transcript. The wet-lab positive results of the test are presented here demonstrating the real-world potential of both hammerhead ribozymes and RiboSoft. RiboSoft is freely available at the website http://ribosoft.fungalgenomics.ca/ribosoft/.

Collaborator of alternative reading frame protein (CARF) associates directly with ARF, p53, and/or human double minute 2 protein (HDM2), a ubiquitin-protein ligase, without cofactors and regulates cell proliferation by forming a negative feedback loop. Although ARF, p53, and HDM2 also participate in the regulation of ribosome biogenesis, the involvement of CARF in this process remains unexplored. In this study, we demonstrate that CARF associates with 5'-3' exoribonuclease 2 (XRN2), which plays a major role in both the maturation of rRNA and the degradation of a variety of discarded pre-rRNA species. We show that overexpression of CARF increases the localization of XRN2 in the nucleoplasm and a concomitant suppression of pre-rRNA processing that leads to accumulation of the 5' extended from of 45S/47S pre-rRNA and 5'-01, A0-1 and E-2 fragments of pre-rRNA transcript in the nucleolus. This was also observed upon XRN2 knockdown. Knockdown of CARF increased the amount of XRN2 in the nucleolar fraction as determined by cell fractionation and by immnocytochemical analysis. These observations suggest that CARF regulates early steps of pre-rRNA processing during ribosome biogenesis by controlling spatial distribution of XRN2 between the nucleoplasm and nucleolus.

Life sciences are yielding huge data sets that underpin scientific discoveries fundamental to improvement in human health, agriculture and the environment. In support of these discoveries, a plethora of databases and tools are deployed, in technically complex and diverse implementations, across a spectrum of scientific disciplines. The corpus of documentation of these resources is fragmented across the Web, with much redundancy, and has lacked a common standard of information. The outcome is that scientists must often struggle to find, understand, compare and use the best resources for the task at hand.Here we present a community-driven curation effort, supported by ELIXIR-the European infrastructure for biological information-that aspires to a comprehensive and consistent registry of information about bioinformatics resources. The sustainable upkeep of this Tools and Data Services Registry is assured by a curation effort driven by and tailored to local needs, and shared amongst a network of engaged partners.As of November 2015, the registry includes 1785 resources, with depositions from 126 individual registrations including 52 institutional providers and 74 individuals. With community support, the registry can become a standard for dissemination of information about bioinformatics resources: we welcome everyone to join us in this common endeavour. The registry is freely available at https://bio.tools.

Cancer stem cells (CSCs) have been identified in several solid malignancies and are now emerging as a plausible target for drug discovery. Beside the questionable existence of CSCs specific markers, the expression of CD133 was reported to be responsible for conferring CSC aggressiveness. Here, we identified two G-rich sequences localized within the introns 3 and 7 of the CD133 gene able to form G-quadruplex (G4) structures, bound and stabilized by small molecules. We further showed that treatment of patient-derived colon CSCs with G4-interacting agents triggers alternative splicing that dramatically impairs the expression of CD133. Interestingly, this is strongly associated with a loss of CSC properties, including self-renewing, motility, tumor initiation and metastases dissemination. Notably, the effects of G4 stabilization on some of these CSC properties are uncoupled from DNA damage response and are fully recapitulated by the selective interference of the CD133 expression.In conclusion, we provided the first proof of the existence of G4 structures within the CD133 gene that can be pharmacologically targeted to impair CSC aggressiveness. This discloses a class of potential antitumoral agents capable of targeting the CSC subpopulation within the tumoral bulk.

IID (Integrated Interactions Database) is the first database providing tissue-specific protein-protein interactions (PPIs) for model organisms and human. IID covers six species (S. cerevisiae (yeast), C. elegans (worm), D. melonogaster (fly), R. norvegicus (rat), M. musculus (mouse) and H. sapiens (human)) and up to 30 tissues per species. Users query IID by providing a set of proteins or PPIs from any of these organisms, and specifying species and tissues where IID should search for interactions. If query proteins are not from the selected species, IID enables searches across species and tissues automatically by using their orthologs; for example, retrieving interactions in a given tissue, conserved in human and mouse. Interaction data in IID comprises three types of PPI networks: experimentally detected PPIs from major databases, orthologous PPIs and high-confidence computationally predicted PPIs. Interactions are assigned to tissues where their proteins pairs or encoding genes are expressed. IID is a major replacement of the I2D interaction database, with larger PPI networks (a total of 1,566,043 PPIs among 68,831 proteins), tissue annotations for interactions, and new query, analysis and data visualization capabilities. IID is available at http://ophid.utoronto.ca/iid.

DNA polymerase gamma (POLG) is essential for replication and repair of mitochondrial DNA (mtDNA). Mutations in POLG cause mtDNA instability and a diverse range of poorly understood human diseases. Here, we created a unique Polg animal model, by modifying polg within the critical and highly conserved polymerase domain in zebrafish. polg(+/-) offspring were indistinguishable from WT siblings in multiple phenotypic and biochemical measures. However, polg(-/-) mutants developed severe mtDNA depletion by one week post-fertilization (wpf), developed slowly and had regenerative defects, yet surprisingly survived up to 4 wpf. An in vivo mtDNA polymerase activity assay utilizing ethidium bromide (EtBr) to deplete mtDNA, showed that polg(+/-) and WT zebrafish fully recover mtDNA content two weeks post-EtBr removal. EtBr further reduced already low levels of mtDNA in polg(-/-) animals, but mtDNA content did not recover following release from EtBr. Despite significantly decreased respiration that corresponded with tissue-specific levels of mtDNA, polg(-/-) animals had WT levels of ATP and no increase in lactate. This zebrafish model of mitochondrial disease now provides unique opportunities for studying mtDNA instability from multiple angles, as polg(-/-) mutants can survive to juvenile stage, rather than lose viability in embryogenesis as seen in Polg mutant mice.

Rice is the most important staple food for a large part of the world's human population and also a key model organism for plant research. Here, we present Information Commons for Rice (IC4R; http://ic4r.org), a rice knowledgebase featuring adoption of an extensible and sustainable architecture that integrates multiple omics data through community-contributed modules. Each module is developed and maintained by different committed groups, deals with data collection, processing and visualization, and delivers data on-demand via web services. In the current version, IC4R incorporates a variety of rice data through multiple committed modules, including genome-wide expression profiles derived entirely from RNA-Seq data, resequencing-based genomic variations obtained from re-sequencing data of thousands of rice varieties, plant homologous genes covering multiple diverse plant species, post-translational modifications, rice-related literatures and gene annotations contributed by the rice research community. Unlike extant related databases, IC4R is designed for scalability and sustainability and thus also features collaborative integration of rice data and low costs for database update and maintenance. Future directions of IC4R include incorporation of other omics data and association of multiple omics data with agronomically important traits, dedicating to build IC4R into a valuable knowledgebase for both basic and translational researches in rice.

RNase H1-dependent antisense oligonucleotides (ASOs) are chemically modified to enhance pharmacological properties. Major modifications include phosphorothioate (PS) backbone and different 2'-modifications in 2-5 nucleotides at each end (wing) of an ASO. Chemical modifications can affect protein binding and understanding ASO-protein interactions is important for better drug design. Recently we identified many intracellular ASO-binding proteins and found that protein binding could affect ASO potency. Here, we analyzed the structure-activity-relationships of ASO-protein interactions and found 2'-modifications significantly affected protein binding, including La, P54nrb and NPM. PS-ASOs containing more hydrophobic 2'-modifications exhibit higher affinity for proteins in general, although certain proteins, e.g. Ku70/Ku80 and TCP1, are less affected by 2'-modifications. We found that Hsp90 protein binds PS-ASOs containing locked-nucleic-acid (LNA) or constrained-ethyl-bicyclic-nucleic-acid ((S)-cEt) modifications much more avidly than 2'-O-methoxyethyl (MOE). ASOs bind the mid-domain of Hsp90 protein. Hsp90 interacts with more hydrophobic 2' modifications, e.g. (S)-cEt or LNA, in the 5'-wing of the ASO. Reduction of Hsp90 protein decreased activity of PS-ASOs with 5'-LNA or 5'-cEt wings, but not with 5'-MOE wing. Together, our results indicate Hsp90 protein enhances the activity of PS/LNA or PS/(S)-cEt ASOs, and imply that altering protein binding of ASOs using different chemical modifications can improve therapeutic performance of PS-ASOs.

Initiation of DNA replication depends upon recognition of genomic sites, termed origins, by AAA+ ATPases. In prokaryotes a single factor binds each origin, whereas in eukaryotes this role is played by a six-protein origin recognition complex (ORC). Why eukaryotes evolved a multisubunit initiator, and the roles of each component, remains unclear. In Trypanosoma brucei, an ancient unicellular eukaryote, only one ORC-related initiator, TbORC1/CDC6, has been identified by sequence homology. Here we show that three TbORC1/CDC6-interacting factors also act in T. brucei nuclear DNA replication and demonstrate that TbORC1/CDC6 interacts in a high molecular complex in which a diverged Orc4 homologue and one replicative helicase subunit can also be found. Analysing the subcellular localization of four TbORC1/CDC6-interacting factors during the cell cycle reveals that one factor, TbORC1B, is not a static constituent of ORC but displays S-phase restricted nuclear localization and expression, suggesting it positively regulates replication. This work shows that ORC architecture and regulation are diverged features of DNA replication initiation in T. brucei, providing new insight into this key stage of eukaryotic genome copying.

ISL1 is expressed in cardiac progenitor cells and plays critical roles in cardiac lineage differentiation and heart development. Cardiac progenitor cells hold great potential for clinical and translational applications. However, the mechanisms underlying ISL1 function in cardiac progenitor cells have not been fully elucidated. Here we uncover a hierarchical role of ISL1 in cardiac progenitor cells, showing that ISL1 directly regulates hundreds of potential downstream target genes that are implicated in cardiac differentiation, through an epigenetic mechanism. Specifically, ISL1 promotes the demethylation of tri-methylation of histone H3K27 (H3K27me3) at the enhancers of key downstream target genes, including Myocd and Mef2c, which are core cardiac transcription factors. ISL1 physically interacts with JMJD3, a H3K27me3 demethylase, and conditional depletion of JMJD3 leads to impaired cardiac progenitor cell differentiation, phenocopying that of ISL1 depletion. Interestingly, ISL1 is not only responsible for the recruitment of JMJD3 to specific target loci during cardiac progenitor differentiation, but also modulates its demethylase activity. In conclusion, ISL1 and JMJD3 partner to alter the cardiac epigenome, instructing gene expression changes that drive cardiac differentiation.

Organisms use molecular chaperones to combat the unfolding and aggregation of proteins. While protein chaperones have been widely studied, here we demonstrate that DNA and RNA exhibit potent chaperone activity in vitro Nucleic acids suppress the aggregation of classic chaperone substrates up to 300-fold more effectively than the protein chaperone GroEL. Additionally, RNA cooperates with the DnaK chaperone system to refold purified luciferase. Our findings reveal a possible new role for nucleic acids within the cell: that nucleic acids directly participate in maintaining proteostasis by preventing protein aggregation.

Persistent hepatitis C virus (HCV) infection leads to chronic hepatitis C (CHC), which often progresses to liver cirrhosis (LC) and hepatocellular carcinoma (HCC). The molecular mechanisms that establish CHC and cause its subsequent development into LC and HCC are poorly understood. We have identified a cytoplasmic double-stranded RNA binding protein, Stau1, which is crucial for HCV replication. In this study, Stau1 specifically interacted with the variable-stem-loop region in the 3' NTR and domain IIId of the HCV-IRES in the 5' NTR, and promoted HCV replication and translation. Stau1 coimmunoprecipitates HCV NS5B and a cell factor, protein kinase R (PKR), which is critical for interferon-induced cellular antiviral and antiproliferative responses. Like Stau1, PKR displayed binding specificity to domain IIId of HCV-IRES. Stau1 binds to PKR and strongly inhibits PKR-autophosphorylation. We demonstrated that the transport of HCV RNA on the polysomes is Stau1-dependent, being mainly localized in the monosome fractions when Stau1 is downregulated and exclusively localized in the polysomes when Stau1 is overexpressed. Our findings suggest that HCV may appropriate Stau1 to its advantage to prevent PKR-mediated inhibition of eIF2α, which is required for the synthesis of HCV proteins for translocation of viral RNA genome to the polysomes for efficient translation and replication.

Ligand-based Virtual Screening (VS) methods aim at identifying molecules with a similar activity profile across phenotypic and macromolecular targets to that of a query molecule used as search template. VS using 3D similarity methods have the advantage of biasing this search toward active molecules with innovative chemical scaffolds, which are highly sought after in drug design to provide novel leads with improved properties over the query molecule (e.g. patentable, of lower toxicity or increased potency). Ultrafast Shape Recognition (USR) has demonstrated excellent performance in the discovery of molecules with previously-unknown phenotypic or target activity, with retrospective studies suggesting that its pharmacophoric extension (USRCAT) should obtain even better hit rates once it is used prospectively. Here we present USR-VS (http://usr.marseille.inserm.fr/), the first web server using these two validated ligand-based 3D methods for large-scale prospective VS. In about 2 s, 93.9 million 3D conformers, expanded from 23.1 million purchasable molecules, are screened and the 100 most similar molecules among them in terms of 3D shape and pharmacophoric properties are shown. USR-VS functionality also provides interactive visualization of the similarity of the query molecule against the hit molecules as well as vendor information to purchase selected hits in order to be experimentally tested.

Epistasis plays an essential role in the development of complex diseases. Interaction methods face common challenge of seeking a balance between persistent power, model complexity, computation efficiency, and validity of identified bio-markers. We introduce a novel W-test to identify pairwise epistasis effect, which measures the distributional difference between cases and controls through a combined log odds ratio. The test is model-free, fast, and inherits a Chi-squared distribution with data adaptive degrees of freedom. No permutation is needed to obtain the P-values. Simulation studies demonstrated that the W-test is more powerful in low frequency variants environment than alternative methods, which are the Chi-squared test, logistic regression and multifactor-dimensionality reduction (MDR). In two independent real bipolar disorder genome-wide associations (GWAS) datasets, the W-test identified significant interactions pairs that can be replicated, including SLIT3-CENPN, SLIT3-TMEM132D, CNTNAP2-NDST4 and CNTCAP2-RTN4R The genes in the pairs play central roles in neurotransmission and synapse formation. A majority of the identified loci are undiscoverable by main effect and are low frequency variants. The proposed method offers a powerful alternative tool for mapping the genetic puzzle underlying complex disorders.