TOPCONS (http://topcons.net/) is a widely used web server for consensus prediction of membrane protein topology. We hereby present a major update to the server, with some substantial improvements, including the following: (i) TOPCONS can now efficiently separate signal peptides from transmembrane regions. (ii) The server can now differentiate more successfully between globular and membrane proteins. (iii) The server now is even slightly faster, although a much larger database is used to generate the multiple sequence alignments. For most proteins, the final prediction is produced in a matter of seconds. (iv) The user-friendly interface is retained, with the additional feature of submitting batch files and accessing the server programmatically using standard interfaces, making it thus ideal for proteome-wide analyses. Indicatively, the user can now scan the entire human proteome in a few days. (v) For proteins with homology to a known 3D structure, the homology-inferred topology is also displayed. (vi) Finally, the combination of methods currently implemented achieves an overall increase in performance by 4% as compared to the currently available best-scoring methods and TOPCONS is the only method that can identify signal peptides and still maintain a state-of-the-art performance in topology predictions.

Global network modeling of distal regulatory interactions is essential in understanding the overall architecture of gene expression programs. Here, we developed a Bayesian probabilistic model and computational method for global causal network construction with breast cancer as a model. Whereas physical regulator binding was well supported by gene expression causality in general, distal elements in intragenic regions or loci distant from the target gene exhibited particularly strong functional effects. Modeling the action of long-range enhancers was critical in recovering true biological interactions with increased coverage and specificity overall and unraveling regulatory complexity underlying tumor subclasses and drug responses in particular. Transcriptional cancer drivers and risk genes were discovered based on the network analysis of somatic and genetic cancer-related DNA variants. Notably, we observed that the risk genes were functionally downstream of the cancer drivers and were selectively susceptible to network perturbation by tumorigenic changes in their upstream drivers. Furthermore, cancer risk alleles tended to increase the susceptibility of the transcription of their associated genes. These findings suggest that transcriptional cancer drivers selectively induce a combinatorial misregulation of downstream risk genes, and that genetic risk factors, mostly residing in distal regulatory regions, increase transcriptional susceptibility to upstream cancer-driving somatic changes.

Small silencing RNAs, including microRNAs, endogenous small interfering RNAs (endo-siRNAs) and Piwi-interacting RNAs (piRNAs), have been shown to play important roles in fine-tuning gene expression, defending virus and controlling transposons. Loss of small silencing RNAs or components in their pathways often leads to severe developmental defects, including lethality and sterility. Recently, non-templated addition of nucleotides to the 3' end, namely tailing, was found to associate with the processing and stability of small silencing RNAs. Next Generation Sequencing has made it possible to detect such modifications at nucleotide resolution in an unprecedented throughput. Unfortunately, detecting such events from millions of short reads confounded by sequencing errors and RNA editing is still a tricky problem. Here, we developed a computational framework, Tailor, driven by an efficient and accurate aligner specifically designed for capturing the tailing events directly from the alignments without extensive post-processing. The performance of Tailor was fully tested and compared favorably with other general-purpose aligners using both simulated and real datasets for tailing analysis. Moreover, to show the broad utility of Tailor, we used Tailor to reanalyze published datasets and revealed novel findings worth further experimental validation. The source code and the executable binaries are freely available at https://github.com/jhhung/Tailor.

A crucial step in the cellular adaptation to oxygen deficiency is the binding of hypoxia-inducible factors (HIFs) to hypoxia response elements (HREs) of oxygen-regulated genes. Genome-wide HIF-1α/2α/β DNA-binding studies revealed that the majority of HREs reside distant to the promoter regions, but the function of these distal HREs has only been marginally studied in the genomic context. We used chromatin immunoprecipitation (ChIP), gene editing (TALEN) and chromosome conformation capture (3C) to localize and functionally characterize a 82 kb upstream HRE that solely drives oxygen-regulated expression of the newly identified HIF target gene PAG1. PAG1, a transmembrane adaptor protein involved in Src signalling, was hypoxically induced in various cell lines and mouse tissues. ChIP and reporter gene assays demonstrated that the -82 kb HRE regulates PAG1, but not an equally distant gene further upstream, by direct interaction with HIF. Ablation of the consensus HRE motif abolished the hypoxic induction of PAG1 but not general oxygen signalling. 3C assays revealed that the -82 kb HRE physically associates with the PAG1 promoter region, independent of HIF-DNA interaction. These results demonstrate a constitutive interaction between the -82 kb HRE and the PAG1 promoter, suggesting a physiologically important rapid response to hypoxia.

DNA3'pp5'G caps synthesized by the 3'-PO4/5'-OH ligase RtcB have a strong impact on enzymatic reactions at DNA 3'-OH ends. Aprataxin, an enzyme that repairs A5'pp5'DNA ends formed during abortive ligation by classic 3'-OH/5'-PO4 ligases, is also a DNA 3' de-capping enzyme, converting DNAppG to DNA3'p and GMP. By taking advantage of RtcB's ability to utilize certain GTP analogs to synthesize DNAppN caps, we show that aprataxin hydrolyzes inosine and 6-O-methylguanosine caps, but is not adept at removing a deoxyguanosine cap. We report a 1.5 Å crystal structure of aprataxin in a complex with GMP, which reveals that: (i) GMP binds at the same position and in the same anti nucleoside conformation as AMP; and (ii) aprataxin makes more extensive nucleobase contacts with guanine than with adenine, via a hydrogen bonding network to the guanine O6, N1, N2 base edge. Alanine mutations of catalytic residues His147 and His149 abolish DNAppG de-capping activity, suggesting that the 3' de-guanylylation and 5' de-adenylylation reactions follow the same pathway of nucleotidyl transfer through a covalent aprataxin-(His147)-NMP intermediate. Alanine mutation of Asp63, which coordinates the guanosine ribose hydroxyls, impairs DNAppG de-capping.

Proteins that function in the same pathways, protein complexes or the same environmental conditions can show similar patterns of sequence conservation across phylogenetic clades. In species that no longer require a specific protein complex or pathway, these proteins, as a group, tend to be lost or diverge. Analysis of the similarity in patterns of sequence conservation across a large set of eukaryotes can predict functional associations between different proteins, identify new pathway members and reveal the function of previously uncharacterized proteins. We used normalized phylogenetic profiling to predict protein function and identify new pathway members and disease genes. The phylogenetic profiles of tens of thousands conserved proteins in the human, mouse, Caenorhabditis elegans and Drosophila genomes can be queried on our new web server, PhyloGene. PhyloGene provides intuitive and user-friendly platform to query the patterns of conservation across 86 animal, fungal, plant and protist genomes. A protein query can be submitted either by selecting the name from whole-genome protein sets of the intensively studied species or by entering a protein sequence. The graphic output shows the profile of sequence conservation for the query and the most similar phylogenetic profiles for the proteins in the genome of choice. The user can also download this output in numerical form.

Telomerase is the enzyme that maintains the length of telomeres. It is minimally constituted of two components: a core reverse transcriptase protein (hTERT) and an RNA (hTR). Despite its significance as an almost universal cancer target, the understanding of the structure of telomerase and the optimization of specific inhibitors have been hampered by the limited amount of enzyme available. Here, we present a breakthrough method to produce unprecedented amounts of recombinant hTERT and to reconstitute human telomerase with purified components. This system provides a decisive tool to identify regulators of the assembly of this ribonucleoprotein complex. It also enables the large-scale screening of small-molecules capable to interfere with telomerase assembly. Indeed, it has allowed us to identify a compound that inhibits telomerase activity when added prior to the assembly of the enzyme, while it has no effect on an already assembled telomerase. Therefore, the novel system presented here may accelerate the understanding of human telomerase assembly and facilitate the discovery of potent and mechanistically unique inhibitors.

Early steps of gene expression are a composite of promoter recognition, promoter activation, RNA synthesis and RNA processing, and it is known that SUMOylation, a post-translational modification, is involved in transcription regulation. We previously found that SUMO-1 marks chromatin at the proximal promoter regions of some of the most active housekeeping genes during interphase in human cells, but the SUMOylated targets on the chromatin remained unclear. In this study, we found that SUMO-1 marks the promoters of ribosomal protein genes via modification of the Scaffold Associated Factor B (SAFB) protein, and the SUMOylated SAFB stimulated both the binding of RNA polymerase to promoters and pre-mRNA splicing. Depletion of SAFB decreased RNA polymerase II binding to promoters and nuclear processing of the mRNA, though mRNA stability was not affected. This study reveals an unexpected role of SUMO-1 and SAFB in the stimulatory coupling of promoter binding, transcription initiation and RNA processing.

Enhancer of zeste homolog 2 (EZH2) is a key epigenetic regulator that catalyzes the trimethylation of H3K27 and is modulated by post-translational modifications (PTMs). However, the precise regulation of EZH2 PTMs remains elusive. We, herein, report that EZH2 is acetylated by acetyltransferase P300/CBP-associated factor (PCAF) and is deacetylated by deacetylase SIRT1. We identified that PCAF interacts with and acetylates EZH2 mainly at lysine 348 (K348). Mechanistically, K348 acetylation decreases EZH2 phosphorylation at T345 and T487 and increases EZH2 stability without disrupting the formation of polycomb repressive complex 2 (PRC2). Functionally, EZH2 K348 acetylation enhances its capacity in suppression of the target genes and promotes lung cancer cell migration and invasion. Further, elevated EZH2 K348 acetylation in lung adenocarcinoma patients predicts a poor prognosis. Our findings define a new mechanism underlying EZH2 modulation by linking EZH2 acetylation to its phosphorylation that stabilizes EZH2 and promotes lung adenocarcinoma progression.

Genetic instability, provoked by exogenous mutagens, is well linked to initiation of cancer. However, even in unstressed cells, DNA undergoes a plethora of spontaneous alterations provoked by its inherent chemical instability and the intracellular milieu. Base excision repair (BER) is the major cellular pathway responsible for repair of these lesions, and as deficiency in BER activity results in DNA damage it has been proposed that it may trigger the development of sporadic cancers. Nevertheless, experimental evidence for this model remains inconsistent and elusive. Here, we performed a proteomic analysis of BER deficient human cells using stable isotope labelling with amino acids in cell culture (SILAC), and demonstrate that BER deficiency, which induces genetic instability, results in dramatic changes in gene expression, resembling changes found in many cancers. We observed profound alterations in tissue homeostasis, serine biosynthesis, and one-carbon- and amino acid metabolism, all of which have been identified as cancer cell 'hallmarks'. For the first time, this study describes gene expression changes characteristic for cells deficient in repair of endogenous DNA lesions by BER. These expression changes resemble those observed in cancer cells, suggesting that genetically unstable BER deficient cells may be a source of pre-cancerous cells.

Flowering is a hot topic in Plant Biology and important progress has been made in Arabidopsis thaliana toward unraveling the genetic networks involved. The increasing complexity and the explosion of literature however require development of new tools for information management and update. We therefore created an evolutive and interactive database of flowering time genes, named FLOR-ID (Flowering-Interactive Database), which is freely accessible at http://www.flor-id.org. The hand-curated database contains information on 306 genes and links to 1595 publications gathering the work of >4500 authors. Gene/protein functions and interactions within the flowering pathways were inferred from the analysis of related publications, included in the database and translated into interactive manually drawn snapshots.

PhytoPath (www.phytopathdb.org) is a resource for genomic and phenotypic data from plant pathogen species, that integrates phenotypic data for genes from PHI-base, an expertly curated catalog of genes with experimentally verified pathogenicity, with the Ensembl tools for data visualization and analysis. The resource is focused on fungi, protists (oomycetes) and bacterial plant pathogens that have genomes that have been sequenced and annotated. Genes with associated PHI-base data can be easily identified across all plant pathogen species using a BioMart-based query tool and visualized in their genomic context on the Ensembl genome browser. The PhytoPath resource contains data for 135 genomic sequences from 87 plant pathogen species, and 1364 genes curated for their role in pathogenicity and as targets for chemical intervention. Support for community annotation of gene models is provided using the WebApollo online gene editor, and we are working with interested communities to improve reference annotation for selected species.

Dendritic cells (DC) are professional antigen presenting cells that develop from hematopoietic stem cells through successive steps of lineage commitment and differentiation. Multipotent progenitors (MPP) are committed to DC restricted common DC progenitors (CDP), which differentiate into specific DC subsets, classical DC (cDC) and plasmacytoid DC (pDC). To determine epigenetic states and regulatory circuitries during DC differentiation, we measured consecutive changes of genome-wide gene expression, histone modification and transcription factor occupancy during the sequel MPP-CDP-cDC/pDC. Specific histone marks in CDP reveal a DC-primed epigenetic signature, which is maintained and reinforced during DC differentiation. Epigenetic marks and transcription factor PU.1 occupancy increasingly coincide upon DC differentiation. By integrating PU.1 occupancy and gene expression we devised a transcription factor regulatory circuitry for DC commitment and subset specification. The circuitry provides the transcription factor hierarchy that drives the sequel MPP-CDP-cDC/pDC, including Irf4, Irf8, Tcf4, Spib and Stat factors. The circuitry also includes feedback loops inferred for individual or multiple factors, which stabilize distinct stages of DC development and DC subsets. In summary, here we describe the basic regulatory circuitry of transcription factors that drives DC development.

Telomerase is a specialized reverse transcriptase that is responsible for telomere length maintenance. As in other organisms, the minimal components required for an active human telomerase are the template-providing telomerase RNA (hTR) and the enzymatic entity telomerase reverse transcriptase (hTERT). Here, we explored the structure of hTR and the hTERT-induced conformational changes within hTR in living cells. By employing an in vivo DMS chemical probing technique, we showed that the pseudoknot and associated triple helical scaffold form stably in vivo independently of hTERT. In fact, the dimethyl-sulfate (DMS) modification pattern suggests that hTR alone is capable of adopting a conformation that is suited to interact with hTERT. However, in the absence of hTERT the template region of hTR is only weakly accessible to DMS-modifications. The predominant change after binding of hTERT to hTR is the exposure of the template region.

Long non coding RNAs (lncRNAs) have emerged as important regulators of various biological processes. LncRNAs also behave as response elements or targets of signaling pathway(s) mediating cellular function. Wnt signaling is important in regulating mammalian spermatogenesis. Mrhl RNA negatively regulates canonical Wnt pathway and gets down regulated upon Wnt signaling activation in mouse spermatogonial cells. Also, mrhl RNA regulates expression of genes pertaining to Wnt pathway and spermatogenesis by binding to chromatin. In the present study, we delineate the detailed molecular mechanism of Wnt signaling induced mrhl RNA down regulation in mouse spermatogonial cells. Mrhl RNA has an independent transcription unit and our various experiments like Chromatin Immunoprecipitation (in cell line as well as mouse testis) and shRNA mediated down regulation convincingly show that β-catenin and TCF4, which are the key effector proteins of the Wnt signaling pathway are required for down regulation of mrhl RNA. We have identified Ctbp1 as the co-repressor and its occupancy on mrhl RNA promoter depends on both β-catenin and TCF4. Upon Wnt signaling activation, Ctbp1 mediated histone repression marks increase at the mrhl RNA promoter. We also demonstrate that Wnt signaling induced mrhl RNA down regulation results in an up regulation of various meiotic differentiation marker genes.

SMC5/6 is a highly conserved protein complex related to cohesin and condensin, which are the key components of higher-order chromatin structures. The SMC5/6 complex is essential for proliferation in yeast and is involved in replication fork stability and processing. However, the precise mechanism of action of SMC5/6 is not known. Here we present evidence that the NSE1/NSE3/NSE4 sub-complex of SMC5/6 binds to double-stranded DNA without any preference for DNA-replication/recombination intermediates. Mutations of key basic residues within the NSE1/NSE3/NSE4 DNA-binding surface reduce binding to DNA in vitro. Their introduction into the Schizosaccharomyces pombe genome results in cell death or hypersensitivity to DNA damaging agents. Chromatin immunoprecipitation analysis of the hypomorphic nse3 DNA-binding mutant shows a reduced association of fission yeast SMC5/6 with chromatin. Based on our results, we propose a model for loading of the SMC5/6 complex onto the chromatin.

The transcription start site (TSS) determines the length and composition of the 5' UTR and therefore can have a profound effect on translation. Yet, little is known about the mechanism underlying start site selection, particularly from promoters lacking conventional core elements such as TATA-box and Initiator. Here we report a novel mechanism of start site selection in the TATA- and Initiator-less promoter of miR-22, through a strictly localized downstream element termed DTIE and an upstream distal element. Changing the distance between them reduced promoter strength, altered TSS selection and diminished Pol II recruitment. Biochemical assays suggest that DTIE does not serve as a docking site for TFIID, the major core promoter-binding factor. TFIID is recruited to the promoter through DTIE but is dispensable for TSS selection. We determined DTIE consensus and found it to be remarkably prevalent, present at the same TSS downstream location in ≈20.8% of human promoters, the vast majority of which are TATA-less. Analysis of DTIE in the tumor suppressor p53 confirmed a similar function. Our findings reveal a novel mechanism of transcription initiation from TATA-less promoters.

Eukaryotic translation initiation factor eIF4AI, the founding member of DEAD-box helicases, undergoes ATP hydrolysis-coupled conformational changes to unwind mRNA secondary structures during translation initiation. However, the mechanism of its coupled enzymatic activities remains unclear. Here we report that a gating mechanism for Pi release controlled by the inter-domain linker of eIF4AI regulates the coupling between ATP hydrolysis and RNA unwinding. Molecular dynamic simulations and experimental results revealed that, through forming a hydrophobic core with the conserved SAT motif of the N-terminal domain and I357 from the C-terminal domain, the linker gated the release of Pi from the hydrolysis site, which avoided futile hydrolysis cycles of eIF4AI. Further mutagenesis studies suggested this linker also plays an auto-inhibitory role in the enzymatic activity of eIF4AI, which may be essential for its function during translation initiation. Overall, our results reveal a novel regulatory mechanism that controls eIF4AI-mediated mRNA unwinding and can guide further mechanistic studies on other DEAD-box helicases.

Repetitive nucleic acid sequences are often prone to form secondary structures distinct from B-DNA. Prominent examples of such structures are DNA triplexes. We observed that certain intrastrand triplex motifs are highly conserved and abundant in prokaryotic genomes. A systematic search of 5246 different prokaryotic plasmids and genomes for intrastrand triplex motifs was conducted and the results summarized in the ITxF database available online at http://bioinformatics.uni-konstanz.de/utils/ITxF/. Next we investigated biophysical and biochemical properties of a particular G/C-rich triplex motif (TM) that occurs in many copies in more than 260 bacterial genomes by CD and nuclear magnetic resonance spectroscopy as well as in vivo footprinting techniques. A characterization of putative properties and functions of these unusually frequent nucleic acid motifs demonstrated that the occurrence of the TM is associated with a high degree of genomic instability. TM-containing genomic loci are significantly more rearranged among closely related Escherichia coli strains compared to control sites. In addition, we found very high frequencies of TM motifs in certain Enterobacteria and Cyanobacteria that were previously described as genetically highly diverse. In conclusion we link intrastrand triplex motifs with the induction of genomic instability. We speculate that the observed instability might be an adaptive feature of these genomes that creates variation for natural selection to act upon.

The infection of Pseudomonas aeruginosa by the giant bacteriophage phiKZ is resistant to host RNA polymerase (RNAP) inhibitor rifampicin. phiKZ encodes two sets of polypeptides that are distantly related to fragments of the two largest subunits of cellular multisubunit RNAPs. Polypeptides of one set are encoded by middle phage genes and are found in the phiKZ virions. Polypeptides of the second set are encoded by early phage genes and are absent from virions. Here, we report isolation of a five-subunit RNAP from phiKZ-infected cells. Four subunits of this enzyme are cellular RNAP subunits homologs of the non-virion set; the fifth subunit is a protein of unknown function. In vitro, this complex initiates transcription from late phiKZ promoters in rifampicin-resistant manner. Thus, this enzyme is a non-virion phiKZ RNAP responsible for transcription of late phage genes. The phiKZ RNAP lacks identifiable assembly and promoter specificity subunits/factors characteristic for eukaryal, archaeal and bacterial RNAPs and thus provides a unique model for comparative analysis of the mechanism, regulation and evolution of this important class of enzymes.