Competing endogenous RNAs (ceRNAs) are RNA molecules that sequester shared microRNAs (miRNAs) thereby affecting the expression of other targets of the miRNAs. Whether genetic variants in ceRNA can affect its biological function and disease development is still an open question. Here we identified a large number of genetic variants that are associated with ceRNA's function using Geuvaids RNA-seq data for 462 individuals from the 1000 Genomes Project. We call these loci competing endogenous RNA expression quantitative trait loci or 'cerQTL', and found that a large number of them were unexplored in conventional eQTL mapping. We identified many cerQTLs that have undergone recent positive selection in different human populations, and showed that single nucleotide polymorphisms in gene 3΄UTRs at the miRNA seed binding regions can simultaneously regulate gene expression changes in both cis and trans by the ceRNA mechanism. We also discovered that cerQTLs are significantly enriched in traits/diseases associated variants reported from genome-wide association studies in the miRNA binding sites, suggesting that disease susceptibilities could be attributed to ceRNA regulation. Further in vitro functional experiments demonstrated that a cerQTL rs11540855 can regulate ceRNA function. These results provide a comprehensive catalog of functional non-coding regulatory variants that may be responsible for ceRNA crosstalk at the post-transcriptional level.

Topoisomerase (topo) IIα and IIβ maintain genome stability and are targets for anti-tumor drugs. In this study, we demonstrate that the decatenation checkpoint is regulated, not only by topo IIα, as previously reported, but also by topo IIβ. The decatenation checkpoint is most efficient when both isoforms are present. Regulation of this checkpoint and sensitivity to topo II-targeted drugs is influenced by the C-terminal domain (CTD) of the topo II isoforms and by a conserved non-catalytic tyrosine, Y640 in topo IIα and Y656 in topo IIβ. Deletion of most of the CTD of topo IIα, while preserving the nuclear localization signal (NLS), enhances the decatenation checkpoint and sensitivity to topo II-targeted drugs. In contrast, deletion of most of the CTD of topo IIβ, while preserving the NLS, and mutation of Y640 in topo IIα and Y656 in topo IIβ inhibits these activities. Structural studies suggest that the differential impact of the CTD on topo IIα and topo IIβ function may be due to differences in CTD charge distribution and differential alignment of the CTD with reference to transport DNA. Together these results suggest that topo IIα and topo IIβ cooperate to maintain genome stability, which may be distinctly modulated by their CTDs.

DNA-single strand annealing proteins (SSAPs) are recombinases frequently encoded in the genome of many bacteriophages. As SSAPs can promote homologous recombination among DNA substrates with an important degree of divergence, these enzymes are involved both in DNA repair and in the generation of phage mosaicisms. Here, analysing Sak and Sak4 as representatives of two different families of SSAPs present in phages infecting the clinically relevant bacterium Staphylococcus aureus, we demonstrate for the first time that these enzymes are absolutely required for phage reproduction. Deletion of the genes encoding these enzymes significantly reduced phage replication and the generation of infectious particles. Complementation studies revealed that these enzymes are required both in the donor (after prophage induction) and in the recipient strain (for infection). Moreover, our results indicated that to perform their function SSAPs require the activity of their cognate single strand binding (Ssb) proteins. Mutational studies demonstrated that the Ssb proteins are also required for phage replication, both in the donor and recipient strain. In summary, our results expand the functions attributed to the Sak and Sak4 proteins, and demonstrate that both SSAPs and Ssb proteins are essential for the life cycle of temperate staphylococcal phages.

The transcription factor ICP4 from herpes simplex virus has a central role in regulating the gene expression cascade which controls viral infection. Here we present the crystal structure of the functionally essential ICP4 DNA binding domain in complex with a segment from its own promoter, revealing a novel homo-dimeric fold. We also studied the complex in solution by small angle X-Ray scattering, nuclear magnetic resonance and surface-plasmon resonance which indicated that, in addition to the globular domain, a flanking intrinsically disordered region also recognizes DNA. Together the data provides a rationale for the bi-partite nature of the ICP4 DNA recognition consensus sequence as the globular and disordered regions bind synergistically to adjacent DNA motifs. Therefore in common with its eukaryotic host, the viral transcription factor ICP4 utilizes disordered regions to enhance the affinity and tune the specificity of DNA interactions in tandem with a globular domain.

The Fox-1 RNA recognition motif (RRM) domain is an important member of the RRM protein family. We report a 1.8 Å X-ray structure of the free Fox-1 containing six distinct monomers. We use this and the nuclear magnetic resonance (NMR) structure of the Fox-1 protein/RNA complex for molecular dynamics (MD) analyses of the structured hydration. The individual monomers of the X-ray structure show diverse hydration patterns, however, MD excellently reproduces the most occupied hydration sites. Simulations of the protein/RNA complex show hydration consistent with the isolated protein complemented by hydration sites specific to the protein/RNA interface. MD predicts intricate hydration sites with water-binding times extending up to hundreds of nanoseconds. We characterize two of them using NMR spectroscopy, RNA binding with switchSENSE and free-energy calculations of mutant proteins. Both hydration sites are experimentally confirmed and their abolishment reduces the binding free-energy. A quantitative agreement between theory and experiment is achieved for the S155A substitution but not for the S122A mutant. The S155 hydration site is evolutionarily conserved within the RRM domains. In conclusion, MD is an effective tool for predicting and interpreting the hydration patterns of protein/RNA complexes. Hydration is not easily detectable in NMR experiments but can affect stability of protein/RNA complexes.

Mcm10 is an essential eukaryotic factor required for DNA replication. The replication fork helicase is composed of Cdc45, Mcm2-7 and GINS (CMG). DDK is an S-phase-specific kinase required for replication initiation, and the DNA primase-polymerase in eukaryotes is pol α. Mcm10 forms oligomers in vitro, mediated by the coiled-coil domain at the N-terminal region of the protein. We characterized an Mcm10 mutant at the N-terminal Domain (NTD), Mcm10-4A, defective for self-interaction. We found that the Mcm10-4A mutant was defective for stimulating DDK phosphorylation of Mcm2, binding to eighty-nucleotide ssDNA, and recruiting pol α to Mcm2-7 in vitro. Expression of wild-type levels of mcm10-4A resulted in severe growth and DNA replication defects in budding yeast cells, with diminished DDK phosphorylation of Mcm2. We then expressed the mcm10-4A in mcm5-bob1 mutant cells to bypass the defects mediated by diminished stimulation of DDK phosphorylation of Mcm2. Expression of wild-type levels of mcm10-4A in mcm5-bob1 mutant cells resulted in severe growth and DNA replication defects, along with diminished RPA signal at replication origins. We also detected diminished GINS and pol-α recruitment to the Mcm2-7 complex. We conclude that an intact Mcm10 coiled-coil interaction surface is important for origin melting, helicase assembly, and the recruitment of pol α to Mcm2-7.

Bacterial RNA polymerase is able to initiate transcription with adenosine-containing cofactor NAD+, which was proposed to result in a portion of cellular RNAs being 'capped' at the 5' end with NAD+, reminiscent of eukaryotic cap. Here we show that, apart from NAD+, another adenosine-containing cofactor FAD and highly abundant uridine-containing cell wall precursors, UDP-Glucose and UDP-N-acetylglucosamine are efficiently used to initiate transcription in vitro. We show that the affinity to NAD+ and UDP-containing factors during initiation is much lower than their cellular concentrations, and that initiation with them stimulates promoter escape. Efficiency of initiation with NAD+, but not with UDP-containing factors, is affected by amino acids of the Rifampicin-binding pocket, suggesting altered RNA capping in Rifampicin-resistant strains. However, relative affinity to NAD+ does not depend on the -1 base of the template strand, as was suggested earlier. We show that incorporation of mature cell wall precursor, UDP-MurNAc-pentapeptide, is inhibited by region 3.2 of σ subunit, possibly preventing targeting of RNA to the membrane. Overall, our in vitro results propose a wide repertoire of potential bacterial RNA capping molecules, and provide mechanistic insights into their incorporation.

Methylated cytosines (5mCs) are frequently mutated in the genome. However, no studies have yet comprehensively analysed mutation-methylation associations across cancer types. Here we analyse 916 cancer genomes, together with tissue type-specific methylation and replication timing data. We describe a strong mutation-methylation association across colorectal cancer subtypes, most interestingly in samples with microsatellite instability (MSI) or Polymerase epsilon (POLE) exonuclease domain mutations. By analysing genomic regions with differential mismatch repair (MMR) efficiency, we suggest a possible role for MMR in the correction of 5mC deamination events, potentially accounting for the high rate of 5mC mutation accumulation in MSI tumours. Additionally, we propose that mutant POLE asserts a mutator phenotype specifically at 5mCs, and we find coding mutation hotspots in POLE-mutant cancers at highly-methylated CpGs in the tumour-suppressor genes APC and TP53. Finally, using multivariable regression models, we demonstrate that different cancers exhibit distinct mutation-methylation associations, with DNA repair influencing such associations in certain cancer genomes. Taken together, we find differential associations with methylation that are vital for accurately predicting expected mutation loads across cancer types. Our findings reveal links between methylation and common mutation and repair processes, with these mechanisms defining a key part of the mutational landscape of cancer genomes.

A considerable effort has been devoted to retrieve systematically information for genes and proteins as well as relationships between them. Despite the importance of chemical compounds and drugs as a central bio-entity in pharmacological and biological research, only a limited number of freely available chemical text-mining/search engine technologies are currently accessible. Here we present LimTox (Literature Mining for Toxicology), a web-based online biomedical search tool with special focus on adverse hepatobiliary reactions. It integrates a range of text mining, named entity recognition and information extraction components. LimTox relies on machine-learning, rule-based, pattern-based and term lookup strategies. This system processes scientific abstracts, a set of full text articles and medical agency assessment reports. Although the main focus of LimTox is on adverse liver events, it enables also basic searches for other organ level toxicity associations (nephrotoxicity, cardiotoxicity, thyrotoxicity and phospholipidosis). This tool supports specialized search queries for: chemical compounds/drugs, genes (with additional emphasis on key enzymes in drug metabolism, namely P450 cytochromes-CYPs) and biochemical liver markers. The LimTox website is free and open to all users and there is no login requirement. LimTox can be accessed at: http://limtox.bioinfo.cnio.es.

Many cellular proteins are methylated on lysine residues and this has been most intensively studied for histone proteins. Lysine methylations on non-histone proteins are also frequent, but in most cases the functional significance of the methylation event, as well as the identity of the responsible lysine (K) specific methyltransferase (KMT), remain unknown. Several recently discovered KMTs belong to the so-called seven-β-strand (7BS) class of MTases and we have here investigated an uncharacterized human 7BS MTase currently annotated as part of the endothelin converting enzyme 2, but which should be considered a separate enzyme. Combining in vitro enzymology and analyzes of knockout cells, we demonstrate that this MTase efficiently methylates K36 in eukaryotic translation elongation factor 1 alpha (eEF1A) in vitro and in vivo. We suggest that this novel KMT is named eEF1A-KMT4 (gene name EEF1AKMT4), in agreement with the recently established nomenclature. Furthermore, by ribosome profiling we show that the absence of K36 methylation affects translation dynamics and changes translation speed of distinct codons. Finally, we show that eEF1A-KMT4 is part of a novel family of human KMTs, defined by a shared sequence motif in the active site and we demonstrate the importance of this motif for catalytic activity.

Mechanistic studies in erythroid cells indicate that LDB1, as part of a GATA1/TAL1/LMO2 complex, brings erythroid-expressed genes into proximity with enhancers for transcription activation. The role of co-activators in establishing this long-range interaction is poorly understood. Here we tested the contributions of the RNA Pol II pre-initiation complex (PIC), mediator and cohesin to establishment of locus control region (LCR)/β-globin proximity. CRISPR/Cas9 editing of the β-globin promoter to eliminate the RNA Pol II PIC by deleting the TATA-box resulted in loss of transcription, but enhancer-promoter interaction was unaffected. Additional deletion of the promoter GATA1 site eliminated LDB1 complex and mediator occupancy and resulted in loss of LCR/β-globin proximity. To separate the roles of LDB1 and mediator in LCR looping, we expressed a looping-competent but transcription-activation deficient form of LDB1 in LDB1 knock down cells: LCR/β-globin proximity was restored without mediator core occupancy. Further, Cas9-directed tethering of mutant LDB1 to the β-globin promoter forced LCR loop formation in the absence of mediator or cohesin occupancy. Moreover, ENCODE data and our chromatin immunoprecipitation results indicate that cohesin is almost completely absent from validated and predicted LDB1-regulated erythroid enhancer-gene pairs. Thus, lineage specific factors largely mediate enhancer-promoter looping in erythroid cells independent of mediator and cohesin.

Profiling of proteome dynamics is crucial for understanding cellular behavior in response to intrinsic and extrinsic stimuli and maintenance of homeostasis. Over the last 20 years, mass spectrometry (MS) has emerged as the most powerful tool for large-scale identification and characterization of proteins. Bottom-up proteomics, the most common MS-based proteomics approach, has always been challenging in terms of data management, processing, analysis and visualization, with modern instruments capable of producing several gigabytes of data out of a single experiment. Here, we present ProteoSign, a freely available web application, dedicated in allowing users to perform proteomics differential expression/abundance analysis in a user-friendly and self-explanatory way. Although several non-commercial standalone tools have been developed for post-quantification statistical analysis of proteomics data, most of them are not end-user appealing as they often require very stringent installation of programming environments, third-party software packages and sometimes further scripting or computer programming. To avoid this bottleneck, we have developed a user-friendly software platform accessible via a web interface in order to enable proteomics laboratories and core facilities to statistically analyse quantitative proteomics data sets in a resource-efficient manner. ProteoSign is available at http://bioinformatics.med.uoc.gr/ProteoSign and the source code at https://github.com/yorgodillo/ProteoSign.

Gene expression stochasticity plays a major role in biology, creating non-genetic cellular individuality and influencing multiple processes, including differentiation and stress responses. We have addressed the lack of knowledge about posttranscriptional contributions to noise by determining cell-to-cell variations in the abundance of mRNA and reporter protein in yeast. Two types of structural element, a stem-loop and a poly(G) motif, not only inhibit translation initiation when inserted into an mRNA 5΄ untranslated region, but also generate noise. The noise-enhancing effect of the stem-loop structure also remains operational when combined with an upstream open reading frame. This has broad significance, since these elements are known to modulate the expression of a diversity of eukaryotic genes. Our findings suggest a mechanism for posttranscriptional noise generation that will contribute to understanding of the generally poor correlation between protein-level stochasticity and transcriptional bursting. We propose that posttranscriptional stochasticity can be linked to cycles of folding/unfolding of a stem-loop structure, or to interconversion between higher-order structural conformations of a G-rich motif, and have created a correspondingly configured computational model that generates fits to the experimental data. Stochastic events occurring during the ribosomal scanning process can therefore feature alongside transcriptional bursting as a source of noise.

Horizontal gene transfer is a key process in the evolution of bacteria and also represents a source of genetic variation in eukaryotes. Among elements participating in gene transfer, thousands of small (<10 kb) mobile bacterial plasmids that replicate by the rolling circle mechanism represent a driving force in the spread of antibiotic resistances. In general, these plasmids are built as genetic modules that encode a replicase, an antibiotic-resistance determinant, and a relaxase that participates in their conjugative mobilization. Further, they control their relatively high copy number (∼30 copies per genome equivalent) by antisense RNAs alone or combined with a repressor protein. We report here that the MobM conjugative relaxase encoded by the promiscuous plasmid pMV158 participates in regulation of the plasmid copy number by transcriptional repression of the antisense RNA, thus increasing the number of plasmid molecules ready to be horizontally transferred (mobilization) and/or vertically inherited (replication). This type of crosstalk between genetic modules involved in vertical and horizontal gene flow has not been reported before.

Here, we report a webserver for the improved SDM, used for predicting the effects of mutations on protein stability. As a pioneering knowledge-based approach, SDM has been highlighted as the most appropriate method to use in combination with many other approaches. We have updated the environment-specific amino-acid substitution tables based on the current expanded PDB (a 5-fold increase in information), and introduced new residue-conformation and interaction parameters, including packing density and residue depth. The updated server has been extensively tested using a benchmark containing 2690 point mutations from 132 different protein structures. The revised method correlates well against the hypothetical reverse mutations, better than comparable methods built using machine-learning approaches, highlighting the strength of our knowledge-based approach for identifying stabilising mutations. Given a PDB file (a Protein Data Bank file format containing the 3D coordinates of the protein atoms), and a point mutation, the server calculates the stability difference score between the wildtype and mutant protein. The server is available at http://structure.bioc.cam.ac.uk/sdm2.

Ribosomally synthesized and post-translationally modified peptides (RiPPs) constitute a rapidly growing class of natural products with diverse structures and bioactivities. We have developed RiPPMiner, a novel bioinformatics resource for deciphering chemical structures of RiPPs by genome mining. RiPPMiner derives its predictive power from machine learning based classifiers, trained using a well curated database of more than 500 experimentally characterized RiPPs. RiPPMiner uses Support Vector Machine to distinguish RiPP precursors from other small proteins and classify the precursors into 12 sub-classes of RiPPs. For classes like lanthipeptide, cyanobactin, lasso peptide and thiopeptide, RiPPMiner can predict leader cleavage site and complex cross-links between post-translationally modified residues starting from genome sequences. RiPPMiner can identify correct cross-link pattern in a core peptide from among a very large number of combinatorial possibilities. Benchmarking of prediction accuracy of RiPPMiner on a large lanthipeptide dataset indicated high sensitivity, specificity, accuracy and precision. RiPPMiner also provides interfaces for visualization of the chemical structure, downloading of simplified molecular-input line-entry system and searching for RiPPs having similar sequences or chemical structures. The backend database of RiPPMiner provides information about modification system, precursor sequence, leader and core sequence, modified residues, cross-links and gene cluster for more than 500 experimentally characterized RiPPs. RiPPMiner is available at http://www.nii.ac.in/rippminer.html.

Linker DNA conformational variability has been proposed to direct nucleosome array folding into more or less compact chromatin fibers but direct experimental evidence for such models are lacking. Here, we tested this hypothesis by designing nucleosome arrays with A-tracts at specific locations in the nucleosome linkers to induce inward (AT-IN) and outward (AT-OUT) bending of the linker DNA. Using electron microscopy and analytical centrifugation techniques, we observed spontaneous folding of AT-IN nucleosome arrays into highly compact structures, comparable to those induced by linker histone H1. In contrast, AT-OUT nucleosome arrays formed less compact structures with decreased nucleosome interactions similar to wild-type nucleosome arrays. Adding linker histone H1 further increased compaction of the A-tract arrays while maintaining structural differences between them. Furthermore, restriction nuclease digestion revealed a strongly reduced accessibility of nucleosome linkers in the compact AT-IN arrays. Electron microscopy analysis and 3D computational Monte Carlo simulations are consistent with a profound zigzag linker DNA configuration and closer nucleosome proximity in the AT-IN arrays due to inward linker DNA bending. We propose that the evolutionary preferred positioning of A-tracts in DNA linkers may control chromatin higher-order folding and thus influence cellular processes such as gene expression, transcription and DNA repair.

Inhibition of the Akt kinase activates HPV16 late gene expression by reducing HPV16 early polyadenylation and by activating HPV16 late L1 mRNA splicing. We identified 'hot spots' for RNA binding proteins at the early polyA signal and at splice sites on HPV16 late mRNAs. We observed that hnRNP L was associated with sequences at all HPV16 late splice sites and at the early polyA signal. Akt kinase inhibition resulted in hnRNP L dephosphorylation and reduced association of hnRNP L with HPV16 mRNAs. This was accompanied by an increased binding of U2AF65 and Sam68 to HPV16 mRNAs. Furthermore, siRNA knock-down of hnRNP L or Akt induced HPV16 gene expression. Treatment of HPV16 immortalized keratinocytes with Akt kinase inhibitor reduced hnRNP L binding to HPV16 mRNAs and induced HPV16 L1 mRNA production. Finally, deletion of the hnRNP L binding sites in HPV16 subgenomic expression plasmids resulted in activation of HPV16 late gene expression. In conclusion, the Akt kinase inhibits HPV16 late gene expression at the level of RNA processing by controlling the RNA-binding protein hnRNP L. We speculate that Akt kinase activity upholds an intracellular milieu that favours HPV16 early gene expression and suppresses HPV16 late gene expression.

Poly(ADP-ribose) glycohydrolase (PARG) regulates cellular poly(ADP-ribose) (PAR) levels by rapidly cleaving glycosidic bonds between ADP-ribose units. PARG interacts with proliferating cell nuclear antigen (PCNA) and is strongly recruited to DNA damage sites in a PAR- and PCNA-dependent fashion. Here we identified PARG acetylation site K409 that is essential for its interaction with PCNA, its localization within replication foci and its recruitment to DNA damage sites. We found K409 to be part of a non-canonical PIP-box within the PARG disordered regulatory region. The previously identified putative N-terminal PIP-box does not bind PCNA directly but contributes to PARG localization within replication foci. X-ray structure and MD simulations reveal that the PARG non-canonical PIP-box binds PCNA in a manner similar to other canonical PIP-boxes and may represent a new type of PIP-box. While the binding of previously described PIP-boxes is based on hydrophobic interactions, PARG PIP-box binds PCNA via both stabilizing hydrophobic and fine-tuning electrostatic interactions. Our data explain the mechanism of PARG-PCNA interaction through a new PARG PIP-box that exhibits non-canonical sequence properties but a canonical mode of PCNA binding.

DNA mismatch repair (MMR) is a highly-conserved DNA repair mechanism, whose primary role is to remove DNA replication errors preventing them from manifesting as mutations, thereby increasing the overall genome stability. Defects in MMR are associated with increased cancer risk in humans and other organisms. Here, we characterize the interaction between MMR and a proofreading-deficient allele of the human replicative DNA polymerase delta, PolδD316A;E318A, which has a higher capacity for strand displacement DNA synthesis than wild type Polδ. Human cell lines overexpressing PolδD316A;E318A display a mild mutator phenotype, while nuclear extracts of these cells exhibit reduced MMR activity in vitro, and these defects are complemented by overexpression or addition of exogenous human Exonuclease 1 (EXO1). By contrast, another proofreading-deficient mutant, PolδD515V, which has a weaker strand displacement activity, does not decrease the MMR activity as significantly as PolδD316A;E318A. In addition, PolδD515V does not increase the mutation frequency in MMR-proficient cells. Based on our findings, we propose that the proofreading activity restricts the strand displacement activity of Polδ in MMR. This contributes to maintain the nicks required for EXO1 entry, and in this manner ensures the dominance of the EXO1-dependent MMR pathway.