Gene regulation in the postgenomic era: technology takes the wheel

Tumor-derived lactate promotes resistance to bevacizumab treatment by facilitating autophagy enhancer protein RUBCNL expression through histone H3 lysine 18 lactylation (H3K18la) in colorectal cancer

Bevacizumab plays an important role in the first and second line treatment for metastatic colorectal cancer (CRC). And induction of hypoxia and the tumors response to it plays an important role in determining the efficacy of antiangiogenic therapy while the connection between them remains unclear. Here, we found that lactate accumulated in the tumor environment of CRC and acted as substrates for histone lactylation, and this process was further induced by cellular enhanced glycolysis in hypoxia. We determined that CRC patients resistant to bevacizumab treatment presented with elevated levels of histone lactylation and inhibition of histone lactylation efficiently suppressed CRC tumorigenesis, progression and survival in hypoxia. Histone lactylation promoted the transcription of RUBCNL/Pacer, facilitating autophagosome maturation through interacting with BECN1 (beclin 1) and mediating the recruitment and function of the class III phosphatidylinositol 3-kinase complex, which had a crucial role in hypoxic cancer cells proliferation and survival. Moreover, combining inhibition of histone lactylation and macroautophagy/autophagy with bevacizumab treatment demonstrated remarkable treatment efficacy in bevacizumab-resistance patients-derived pre-clinical models. These findings delivered a new exploration and important supplement of metabolic reprogramming-epigenetic regulation, and provided a new strategy for improving clinical efficacy of bevacizumab in CRC by inhibition of histone lactylation.Abbreviations: 2-DG: 2-deoxy-D-glucose; BECN1: beclin 1; CQ: chloroquine; CRC: colorectal cancer; DMOG: dimethyloxalylglycine; H3K18la: histone H3 lysine 18 lactylation; MAP1LC3B/LC3B: microtubule associated protein 1 light chain 3 beta; Nala: sodium lactate; PDO: patient-derived orgnoid; PDX: patient-derived xenograft; RUBCNL/Pacer: rubicon like autophagy enhancer; SQSTM1/p62: sequestosome 1.

Histone H3 proline 16 hydroxylation regulates mammalian gene expression

Histone post-translational modifications (PTMs) are important for regulating various DNA-templated processes. Here, we report the existence of a histone PTM in mammalian cells, namely histone H3 with hydroxylation of proline at residue 16 (H3P16oh), which is catalyzed by the proline hydroxylase EGLN2. We show that H3P16oh enhances direct binding of KDM5A to its substrate, histone H3 with trimethylation at the fourth lysine residue (H3K4me3), resulting in enhanced chromatin recruitment of KDM5A and a corresponding decrease of H3K4me3 at target genes. Genome- and transcriptome-wide analyses show that the EGLN2-KDM5A axis regulates target gene expression in mammalian cells. Specifically, our data demonstrate repression of the WNT pathway negative regulator DKK1 through the EGLN2-H3P16oh-KDM5A pathway to promote WNT/β-catenin signaling in triple-negative breast cancer (TNBC). This study characterizes a regulatory mark in the histone code and reveals a role for H3P16oh in regulating mammalian gene expression.

LGR4/GPR48 inactivation leads to aniridia-genitourinary anomalies-mental retardation syndrome defects

AGR syndrome (the clinical triad of aniridia, genitourinary anomalies, and mental retardation, a subgroup of WAGR syndrome for Wilm's tumor, aniridia, genitourinary anomalies, and mental retardation) is a rare syndrome caused by a contiguous gene deletion in the 11p13-14 region. However, the mechanisms of WAGR syndrome pathogenesis are elusive. In this study we provide evidence that LGR4 (also named GPR48), the only G-protein-coupled receptor gene in the human chromosome 11p12-11p14.4 fragment, is the key gene responsible for the diseases of AGR syndrome. Deletion of Lgr4 in mouse led to aniridia, polycystic kidney disease, genitourinary anomalies, and mental retardation, similar to the pathological defects of AGR syndrome. Furthermore, Lgr4 inactivation significantly increased cell apoptosis and decreased the expression of multiple important genes involved in the development of WAGR syndrome related organs. Specifically, deletion of Lgr4 down-regulated the expression of histone demethylases Jmjd2a and Fbxl10 through cAMP-CREB signaling pathways both in mouse embryonic fibroblast cells and in urinary and reproductive system mouse tissues. Our data suggest that Lgr4, which regulates eye, kidney, testis, ovary, and uterine organ development as well as mental development through genetic and epigenetic surveillance, is a novel candidate gene for the pathogenesis of AGR syndrome.

Solvated protein-DNA docking using HADDOCK

Interfacial water molecules play an important role in many aspects of protein-DNA specificity and recognition. Yet they have been mostly neglected in the computational modeling of these complexes. We present here a solvated docking protocol that allows explicit inclusion of water molecules in the docking of protein-DNA complexes and demonstrate its feasibility on a benchmark of 30 high-resolution protein-DNA complexes containing crystallographically-determined water molecules at their interfaces. Our protocol is capable of reproducing the solvation pattern at the interface and recovers hydrogen-bonded water-mediated contacts in many of the benchmark cases. Solvated docking leads to an overall improvement in the quality of the generated protein-DNA models for cases with limited conformational change of the partners upon complex formation. The applicability of this approach is demonstrated on real cases by docking a representative set of 6 complexes using unbound protein coordinates, model-built DNA and knowledge-based restraints. As HADDOCK supports the inclusion of a variety of NMR restraints, solvated docking is also applicable for NMR-based structure calculations of protein-DNA complexes.

A genome-scale resource for in vivo tag-based protein function exploration in C. elegans

Understanding the in vivo dynamics of protein localization and their physical interactions is important for many problems in biology. To enable systematic protein function interrogation in a multicellular context, we built a genome-scale transgenic platform for in vivo expression of fluorescent- and affinity-tagged proteins in Caenorhabditis elegans under endogenous cis regulatory control. The platform combines computer-assisted transgene design, massively parallel DNA engineering, and next-generation sequencing to generate a resource of 14,637 genomic DNA transgenes, which covers 73% of the proteome. The multipurpose tag used allows any protein of interest to be localized in vivo or affinity purified using standard tag-based assays. We illustrate the utility of the resource by systematic chromatin immunopurification and automated 4D imaging, which produced detailed DNA binding and cell/tissue distribution maps for key transcription factor proteins.

Regulation of antigen receptor gene assembly by genetic-epigenetic crosstalk

Many aspects of gene function are coordinated by changes in the epigenome, which include dynamic revisions of chromatin modifications, genome packaging, subnuclear localization, and chromosome conformation. All of these mechanisms are used by developing lymphocytes to regulate the assembly of functional antigen receptor genes by V(D)J recombination. This somatic rearrangement of the genome must be tightly regulated to ensure proper B and T cell development and to avoid chromosomal translocations that cause lymphoid tumors. V(D)J recombination is controlled by a complex interplay between cis-acting regulatory elements that use transcription factors as liaisons to communicate with epigenetic pathways. Genetic-epigenetic crosstalk is a key strategy employed by precursor lymphocytes to modulate chromatin configurations at Ig and Tcr loci and thereby permit or deny access to a single V(D)J recombinase complex. This article describes our current knowledge of how genetic elements orchestrate crosstalk with epigenetic mechanisms to regulate recombinase accessibility via localized, regional, or long-range changes in chromatin.

Pushing the limits of what is achievable in protein-DNA docking: benchmarking HADDOCK's performance

The intrinsic flexibility of DNA and the difficulty of identifying its interaction surface have long been challenges that prevented the development of efficient protein-DNA docking methods. We have demonstrated the ability our flexible data-driven docking method HADDOCK to deal with these before, by using custom-built DNA structural models. Here we put our method to the test on a set of 47 complexes from the protein-DNA docking benchmark. We show that HADDOCK is able to predict many of the specific DNA conformational changes required to assemble the interface(s). Our DNA analysis and modelling procedure captures the bend and twist motions occurring upon complex formation and uses these to generate custom-built DNA structural models, more closely resembling the bound form, for use in a second docking round. We achieve throughout the benchmark an overall success rate of 94% of one-star solutions or higher (interface root mean square deviation ≤4 A and fraction of native contacts >10%) according to CAPRI criteria. Our improved protocol successfully predicts even the challenging protein-DNA complexes in the benchmark. Finally, our method is the first to readily dock multiple molecules (N > 2) simultaneously, pushing the limits of what is currently achievable in the field of protein-DNA docking.

Universal function-specificity of codon usage

Synonymous codon usage has long been known as a factor that affects average expression level of proteins in fast-growing microorganisms, but neither its role in dynamic changes of expression in response to environmental changes nor selective factors shaping it in the genomes of higher eukaryotes have been fully understood. Here, we propose that codon usage is ubiquitously selected to synchronize the translation efficiency with the dynamic alteration of protein expression in response to environmental and physiological changes. Our analysis reveals that codon usage is universally correlated with gene function, suggesting its potential contribution to synchronized regulation of genes with similar functions. We directly show that coexpressed genes have similar synonymous codon usages within the genomes of human, yeast, Caenorhabditis elegans and Escherichia coli. We also demonstrate that perturbing the codon usage directly affects the level or even direction of changes in protein expression in response to environmental stimuli. Perturbing tRNA composition also has tangible phenotypic effects on the cell. By showing that codon usage is universally function-specific, our results expand, to almost all organisms, the notion that cells may need to dynamically alter their intracellular tRNA composition in order to adapt to their new environment or physiological role.

SetDB1 contributes to repression of genes encoding developmental regulators and maintenance of ES cell state

Transcription factors that play key roles in regulating embryonic stem (ES) cell state have been identified, but the chromatin regulators that help maintain ES cells are less well understood. A high-throughput shRNA screen was used to identify novel chromatin regulators that influence ES cell state. Loss of histone H3 Lys 9 (H3K9) methyltransferases, particularly SetDB1, had the most profound effects on ES cells. Chromatin immunoprecipitation (ChIP) coupled with massively parallel DNA sequencing (ChIP-Seq) and functional analysis revealed that SetDB1 and histone H3K9-methylated nucleosomes occupy and repress genes encoding developmental regulators. These SetDB1-occupied genes are a subset of the "bivalent" genes, which contain nucleosomes with H3K4me3 (H3K4 trimethylation) and H3K27me3 modifications catalyzed by Trithorax and Polycomb group proteins, respectively. These genes are subjected to repression by both Polycomb group proteins and SetDB1, and loss of either regulator can destabilize ES cell state.

Mechanisms of polycomb gene silencing: knowns and unknowns

Polycomb proteins form chromatin-modifying complexes that implement transcriptional silencing in higher eukaryotes. Hundreds of genes are silenced by Polycomb proteins, including dozens of genes that encode crucial developmental regulators in organisms ranging from plants to humans. Two main families of complexes, called Polycomb repressive complex 1 (PRC1) and PRC2, are targeted to repressed regions. Recent studies have advanced our understanding of these complexes, including their potential mechanisms of gene silencing, the roles of chromatin modifications, their means of delivery to target genes and the functional distinctions among variant complexes. Emerging concepts include the existence of a Polycomb barrier to transcription elongation and the involvement of non-coding RNAs in the targeting of Polycomb complexes. These findings have an impact on the epigenetic programming of gene expression in many biological systems.

PDA: an automatic and comprehensive analysis program for protein-DNA complex structures

BACKGROUND

Knowledge of protein-DNA interactions at the structural-level can provide insights into the mechanisms of protein-DNA recognition and gene regulation. Although over 1400 protein-DNA complex structures have been deposited into Protein Data Bank (PDB), the structural details of protein-DNA interactions are generally not available. In addition, current approaches to comparison of protein-DNA complexes are mainly based on protein sequence similarity while the DNA sequences are not taken into account. With the number of experimentally-determined protein-DNA complex structures increasing, there is a need for an automatic program to analyze the protein-DNA complex structures and to provide comprehensive structural information for the benefit of the whole research community.

RESULTS

We developed an automatic and comprehensive protein-DNA complex structure analysis program, PDA (for protein-DNA complex structure analyzer). PDA takes PDB files as inputs and performs structural analysis that includes 1) whole protein-DNA complex structure restoration, especially the reconstruction of double-stranded DNA structures; 2) an efficient new approach for DNA base-pair detection; 3) systematic annotation of protein-DNA interactions; and 4) extraction of DNA subsequences involved in protein-DNA interactions and identification of protein-DNA binding units. Protein-DNA complex structures in current PDB were processed and analyzed with our PDA program and the analysis results were stored in a database. A dataset useful for studying protein-DNA interactions involved in gene regulation was generated using both protein and DNA sequences as well as the contact information of the complexes. WebPDA was developed to provide a web interface for using PDA and for data retrieval.

CONCLUSION

PDA is a computational tool for structural annotations of protein-DNA complexes. It provides a useful resource for investigating protein-DNA interactions. Data from the PDA analysis can also facilitate the classification of protein-DNA complexes and provide insights into rational design of benchmarks. The PDA program is freely available at http://bioinfozen.uncc.edu/webpda.

A novel short hairpin RNA (shRNA) expression system promotes Sox9-dependent gene silencing

Cartilage development and function are dependent on a temporally integrated program of gene expression. With the advent of RNA interference (RNAi), artificial control of these complex programs becomes a possibility, limited only by the ability to regulate and express the RNAi. Using existing methods for production of RNAi's, we have constructed a plasmid-based short hairpin RNA (shRNA) expression system under control of the human pol III H1 promoter and supplemented this promoter with DNA binding sites for the cartilage-specific transcription factor Sox9. The resulting shRNA expression system displays robust, Sox9-dependent gene silencing. Dependence on Sox9 expression was confirmed by electrophoretic mobility shift assays. The ability of the system to regulate heterologously expressed Sox9 was demonstrated by Western blot, as a function of both Sox9 to shRNA ratio, as well as time from transfection. This novel expression system supports auto-regulatory gene silencing, providing a tissue-specific feedback mechanism for temporal control of gene expression. Its applications for both basic mechanistic studies and therapeutic purposes should facilitate the design and implementation of innovative tissue engineering strategies.

H2AZ is enriched at polycomb complex target genes in ES cells and is necessary for lineage commitment

Elucidating how chromatin influences gene expression patterns and ultimately cell fate is fundamental to understanding development and disease. The histone variant H2AZ has emerged as a key regulator of chromatin function and plays an essential but unknown role during mammalian development. Here, genome-wide analysis reveals that H2AZ occupies the promoters of developmentally important genes in a manner that is remarkably similar to that of the Polycomb group (PcG) protein Suz12. By using RNAi, we demonstrate a role for H2AZ in regulating target gene expression, find that H2AZ and PcG protein occupancy is interdependent at promoters, and further show that H2AZ is necessary for ES cell differentiation. Notably, H2AZ occupies a different subset of genes in lineage-committed cells, suggesting that its dynamic redistribution is necessary for cell fate transitions. Thus, H2AZ, together with PcG proteins, may establish specialized chromatin states in ES cells necessary for the proper execution of developmental gene expression programs.