Differential binding of transcription factors (TFs) at cis-regulatory loci drives the differentiation and function of diverse cellular lineages. Understanding the regulatory interactions that underlie cell fate decisions requires characterizing TF binding sites (TFBS) across multiple cell types and conditions. Techniques, e.g. ChIP-Seq can reveal genome-wide patterns of TF binding, but typically requires laborious and costly experiments for each TF-cell-type (TFCT) condition of interest. Chromosomal accessibility assays can connect accessible chromatin in one cell type to many TFs through sequence motif mapping. Such methods, however, rarely take into account that the genomic context preferred by each factor differs from TF to TF, and from cell type to cell type. To address the differences in TF behaviors, we developed Mocap, a method that integrates chromatin accessibility, motif scores, TF footprints, CpG/GC content, evolutionary conservation and other factors in an ensemble of TFCT-specific classifiers. We show that integration of genomic features, such as CpG islands improves TFBS prediction in some TFCT. Further, we describe a method for mapping new TFCT, for which no ChIP-seq data exists, onto our ensemble of classifiers and show that our cross-sample TFBS prediction method outperforms several previously described methods.

Transposable elements (TEs) comprise a large proportion of long non-coding RNAs (lncRNAs). Here, we employed CRISPR to delete a short interspersed nuclear element (SINE) in Malat1, a cancer-associated lncRNA, to investigate its significance in cellular physiology. We show that Malat1 with a SINE deletion forms diffuse nuclear speckles and is frequently translocated to the cytoplasm. SINE-deleted cells exhibit an activated unfolded protein response and PKR and markedly increased DNA damage and apoptosis caused by dysregulation of TDP-43 localization and formation of cytotoxic inclusions. TDP-43 binds stronger to Malat1 without the SINE and is likely 'hijacked' by cytoplasmic Malat1 to the cytoplasm, resulting in the depletion of nuclear TDP-43 and redistribution of TDP-43 binding to repetitive element transcripts and mRNAs encoding mitotic and nuclear-cytoplasmic regulators. The SINE promotes Malat1 nuclear retention by facilitating Malat1 binding to HNRNPK, a protein that drives RNA nuclear retention, potentially through direct interactions of the SINE with KHDRBS1 and TRA2A, which bind to HNRNPK. Losing these RNA-protein interactions due to the SINE deletion likely creates more available TDP-43 binding sites on Malat1 and subsequent TDP-43 aggregation. These results highlight the significance of lncRNA TEs in TDP-43 proteostasis with potential implications in both cancer and neurodegenerative diseases.

NrS-1 is the first known phage that can infect Epsilonproteobacteria, one of the predominant primary producers in the deep-sea hydrothermal vent ecosystems. NrS-1 polymerase is a multidomain enzyme and is one key component of the phage replisome. The N-terminal Prim/Pol and HBD domains are responsible for DNA polymerization and de novo primer synthesis activities of NrS-1 polymerase. However, the structure and function of the C-terminus (CTR) of NrS-1 polymerase are poorly understood. Here, we report two crystal structures, showing that NrS-1 CTR adopts one unique hexameric ring-shaped conformation. Although the central helicase domain of NrS-1 CTR shares structural similarity with the superfamily III helicases, the folds of the Head and Tail domains are completely novel. Via mutagenesis and in vitro biochemical analysis, we identified many residues important for the helicase and polymerization activities of NrS-1 polymerase. In addition to NrS-1 polymerase, our study may also help us identify and understand the functions of multidomain polymerases expressed by many NrS-1 related phages.

The advent of single-cell sequencing opened a new era in transcriptomic and genomic research. To understand cell composition using single-cell studies, a variety of cell markers have been widely used to label individual cell types. However, the specific database of cell markers for use by the plant research community remains very limited. To overcome this problem, we developed the Plant Cell Marker DataBase (PCMDB, http://www.tobaccodb.org/pcmdb/), which is based on a uniform annotation pipeline. By manually curating over 130 000 research publications, we collected a total of 81 117 cell marker genes of 263 cell types in 22 tissues across six plant species. Tissue- and cell-specific expression patterns can be visualized using multiple tools: eFP Browser, Bar, and UMAP/TSNE graph. The PCMDB also supports several analysis tools, including SCSA and SingleR, which allows for user annotation of cell types. To provide information about plant species currently unsupported in PCMDB, potential marker genes for other plant species can be searched based on homology with the supported species. PCMDB is a user-friendly hierarchical platform that contains five built-in search engines. We believe PCMDB will constitute a useful resource for researchers working on cell type annotation and the prediction of the biological function of individual cells.

Single same cell RNAseq/ATACseq multiome data provide unparalleled potential to develop high resolution maps of the cell-type specific transcriptional regulatory circuitry underlying gene expression. We present CREMA, a framework that recovers the full cis-regulatory circuitry by modeling gene expression and chromatin activity in individual cells without peak-calling or cell type labeling constraints. We demonstrate that CREMA overcomes the limitations of existing methods that fail to identify about half of functional regulatory elements which are outside the called chromatin 'peaks'. These circuit sites outside called peaks are shown to be important cell type specific functional regulatory loci, sufficient to distinguish individual cell types. Analysis of mouse pituitary data identifies a Gata2-circuit for the gonadotrope-enriched disease-associated Pcsk1 gene, which is experimentally validated by reduced gonadotrope expression in a gonadotrope conditional Gata2-knockout model. We present a web accessible human immune cell regulatory circuit resource, and provide CREMA as an R package.

To determine early somatic changes in high-grade serous ovarian cancer (HGSOC), we performed whole genome sequencing on a rare collection of 16 low stage HGSOCs. The majority showed extensive structural alterations (one had an ultramutated profile), exhibited high levels of p53 immunoreactivity, and harboured a TP53 mutation, deletion or inactivation. BRCA1 and BRCA2 mutations were observed in two tumors, with nine showing evidence of a homologous recombination (HR) defect. Combined Analysis with The Cancer Genome Atlas (TCGA) indicated that low and late stage HGSOCs have similar mutation and copy number profiles. We also found evidence that deleterious TP53 mutations are the earliest events, followed by deletions or loss of heterozygosity (LOH) of chromosomes carrying TP53, BRCA1 or BRCA2. Inactivation of HR appears to be an early event, as 62.5% of tumours showed a LOH pattern suggestive of HR defects. Three tumours with the highest ploidy had little genome-wide LOH, yet one of these had a homozygous somatic frame-shift BRCA2 mutation, suggesting that some carcinomas begin as tetraploid then descend into diploidy accompanied by genome-wide LOH. Lastly, we found evidence that structural variants (SV) cluster in HGSOC, but are absent in one ultramutated tumor, providing insights into the pathogenesis of low stage HGSOC.

Transcription factors belonging to the same transcription factor families contain very similar DNA binding domains and hence have the potential to bind to related DNA sequences. However, subtle differences in binding specificities can be detected in vitro with the potential to direct specific responses in vivo. Here, we have examined the binding properties of three Forkhead (FOX) transcription factors, FOXK2, FOXO3 and FOXJ3 in vivo. Extensive overlap in chromatin binding is observed, although underlying differential DNA binding specificity can dictate the recruitment of FOXK2 and FOXJ3 to chromatin. However, functionally, FOXO3-dependent gene regulation is generally mediated not through uniquely bound regions but through regions occupied by both FOXK2 and FOXO3 where both factors play a regulatory role. Our data point to a model whereby FOX transcription factors control gene expression through dynamically binding and generating partial occupancy of the same site rather than mutually exclusive binding derived by stable binding of individual FOX proteins.

Methylation of miRNAs at the 2'-hydroxyl group on the ribose at 3'-end (2'-O-methylation, 2'Ome) is critical for miRNA function in plants and Drosophila. Whether this methylation phenomenon exists for mammalian miRNA remains unknown. Through LC-MS/MS analysis, we discover that majority of miR-21-5p isolated from human non-small cell lung cancer (NSCLC) tissue possesses 3'-terminal 2'Ome. Predominant 3'-terminal 2'Ome of miR-21-5p in cancer tissue is confirmed by qRT-PCR and northern blot after oxidation/β-elimination procedure. Cancerous and the paired non-cancerous lung tissue miRNAs display different pattern of 3'-terminal 2'Ome. We further identify HENMT1 as the methyltransferase responsible for 3'-terminal 2'Ome of mammalian miRNAs. Compared to non-methylated miR-21-5p, methylated miR-21-5p is more resistant to digestion by 3'→5' exoribonuclease polyribonucleotide nucleotidyltransferase 1 (PNPT1) and has higher affinity to Argonaute-2, which may contribute to its higher stability and stronger inhibition on programmed cell death protein 4 (PDCD4) translation, respectively. Our findings reveal HENMT1-mediated 3'-terminal 2'Ome of mammalian miRNAs and highlight its role in enhancing miRNA's stability and function.

RPA is a critical factor for DNA replication and replication stress response. Surprisingly, we found that chromatin RPA stability is tightly regulated. We report that the GDP/GTP exchange factor DOCK7 acts as a critical replication stress regulator to promote RPA stability on chromatin. DOCK7 is phosphorylated by ATR and then recruited by MDC1 to the chromatin and replication fork during replication stress. DOCK7-mediated Rac1/Cdc42 activation leads to the activation of PAK1, which subsequently phosphorylates RPA1 at S135 and T180 to stabilize chromatin-loaded RPA1 and ensure proper replication stress response. Moreover, DOCK7 is overexpressed in ovarian cancer and depleting DOCK7 sensitizes cancer cells to camptothecin. Taken together, our results highlight a novel role for DOCK7 in regulation of the replication stress response and highlight potential therapeutic targets to overcome chemoresistance in cancer.

Mitochondrial nucleases play important roles in accurate maintenance and correct metabolism of mtDNA, the own genetic materials of mitochondria that are passed exclusively from mother to child. MGME1 is a highly conserved DNase that was discovered recently. Mutations in MGME1-coding gene lead to severe mitochondrial syndromes characterized by external ophthalmoplegia, emaciation, and respiratory failure in humans. Unlike many other nucleases that are distributed in multiple cellular organelles, human MGME1 is a mitochondria-specific nuclease; therefore, it can serve as an ideal target for treating related syndromes. Here, we report one HsMGME1-Mn2+ complex and three different HsMGME1-DNA complex structures. In combination with in vitro cleavage assays, our structures reveal the detailed molecular basis for substrate DNA binding and/or unwinding by HsMGME1. Besides the conserved two-cation-assisted catalytic mechanism, structural analysis of HsMGME1 and comparison with homologous proteins also clarified substrate binding and cleavage directionalities of the DNA double-strand break repair complexes RecBCD and AddAB.

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

Chromatin immunoprecipitation with massively parallel DNA sequencing (ChIP-seq) has greatly improved the reliability with which transcription factor binding sites (TFBSs) can be identified from genome-wide profiling studies. Many computational tools are developed to detect binding events or peaks, however the robust detection of weak binding events remains a challenge for current peak calling tools. We have developed a novel Bayesian approach (ChIP-BIT) to reliably detect TFBSs and their target genes by jointly modeling binding signal intensities and binding locations of TFBSs. Specifically, a Gaussian mixture model is used to capture both binding and background signals in sample data. As a unique feature of ChIP-BIT, background signals are modeled by a local Gaussian distribution that is accurately estimated from the input data. Extensive simulation studies showed a significantly improved performance of ChIP-BIT in target gene prediction, particularly for detecting weak binding signals at gene promoter regions. We applied ChIP-BIT to find target genes from NOTCH3 and PBX1 ChIP-seq data acquired from MCF-7 breast cancer cells. TF knockdown experiments have initially validated about 30% of co-regulated target genes identified by ChIP-BIT as being differentially expressed in MCF-7 cells. Functional analysis on these genes further revealed the existence of crosstalk between Notch and Wnt signaling pathways.

Hepatic miR-122 can serve as a pro-apoptotic factor to suppress tumorigenesis. The underlying mechanism, however, remains incompletely understood. Here we present the first evidence that miR-122 promotes hepatocellular carcinoma cell apoptosis through directly silencing the biogenesis of cell survival oncomiR miR-21 at posttranscriptional level. We find that miR-122 is strongly expressed in primary liver cell nucleus but its nuclear localization is markedly decreased in transformed cells particularly in chemoresistant tumor cells. MiRNA profiling and RT-qPCR confirm an inverse correlation between miR-122 and miR-21 in hepatocellular carcinoma tissues/cells, and increasing or decreasing nuclear level of miR-122 respectively reduces or increases miR-21 expression. Mechanistically, nuclear miR-122 suppresses miR-21 maturation via binding to a 19-nt UG-containing recognition element in the basal region of pri-miR-21 and preventing the Drosha-DGCR8 microprocessor's conversion of pri-miR-21 into pre-miR-21. Furthermore, both in vitro and in vivo studies demonstrate that nuclear miR-122 participates in the regulation of HCC cell apoptosis through modulating the miR-21-targeted programmed cell death 4 (PDCD4) signal pathway.

Insects are the largest group of animals on the planet and have a huge impact on human life by providing resources, transmitting diseases, and damaging agricultural crop production. Recently, a large amount of insect genome and gene data has been generated. A comprehensive database is highly desirable for managing, sharing, and mining these resources. Here, we present an updated database, InsectBase 2.0 (http://v2.insect-genome.com/), covering 815 insect genomes, 25 805 transcriptomes and >16 million genes, including 15 045 111 coding sequences, 3 436 022 3'UTRs, 4 345 664 5'UTRs, 112 162 miRNAs and 1 293 430 lncRNAs. In addition, we used an in-house standard pipeline to annotate 1 434 653 genes belonging to 164 gene families; 215 986 potential horizontally transferred genes; and 419 KEGG pathways. Web services such as BLAST, JBrowse2 and Synteny Viewer are provided for searching and visualization. InsectBase 2.0 serves as a valuable platform for entomologists and researchers in the related communities of animal evolution and invertebrate comparative genomics.

Infection with kinetoplastid parasites, including Trypanosoma brucei (T. brucei), Trypanosoma cruzi (T. cruzi) and Leishmania can cause serious disease in humans. Like other kinetoplastid species, mRNAs of these disease-causing parasites must undergo posttranscriptional editing in order to be functional. mRNA editing is directed by gRNAs, a large group of small RNAs. Similar to mRNAs, gRNAs are also precisely regulated. In T. brucei, overexpression of RNase D ribonuclease (TbRND) leads to substantial reduction in the total gRNA population and subsequent inhibition of mRNA editing. However, the mechanisms regulating gRNA binding and cleavage by TbRND are not well defined. Here, we report a thorough structural study of TbRND. Besides Apo- and NMP-bound structures, we also solved one TbRND structure in complexed with single-stranded RNA. In combination with mutagenesis and in vitro cleavage assays, our structures indicated that TbRND follows the conserved two-cation-assisted mechanism in catalysis. TbRND is a unique RND member, as it contains a ZFD domain at its C-terminus. In addition to T. brucei, our studies also advanced our understanding on the potential gRNA degradation pathway in T. cruzi, Leishmania, as well for as other disease-associated parasites expressing ZFD-containing RNDs.

Bromodomain-containing protein 9 (BRD9) is a specific subunit of the non-canonical SWI/SNF (ncBAF) chromatin-remodeling complex, whose function in human embryonic stem cells (hESCs) remains unclear. Here, we demonstrate that impaired BRD9 function reduces the self-renewal capacity of hESCs and alters their differentiation potential. Specifically, BRD9 depletion inhibits meso-endoderm differentiation while promoting neural ectoderm differentiation. Notably, supplementation of NODAL, TGF-β, Activin A or WNT3A rescues the differentiation defects caused by BRD9 loss. Mechanistically, BRD9 forms a complex with BRD4, SMAD2/3, β-CATENIN and P300, which regulates the expression of pluripotency genes and the activity of TGF-β/Nodal/Activin and Wnt signaling pathways. This is achieved by regulating the deposition of H3K27ac on associated genes, thus maintaining and directing hESC differentiation. BRD9-mediated regulation of the TGF-β/Activin/Nodal pathway is also demonstrated in the development of pancreatic and breast cancer cells. In summary, our study highlights the crucial role of BRD9 in the regulation of hESC self-renewal and differentiation, as well as its participation in the progression of pancreatic and breast cancers.

Archaeal ribosomes have many domain-specific features; however, our understanding of these structures is limited. We present 10 cryo-electron microscopy (cryo-EM) structures of the archaeal ribosome from crenarchaeota Sulfolobus acidocaldarius (Sac) at 2.7-5.7 Å resolution. We observed unstable conformations of H68 and h44 of ribosomal RNA (rRNA) in the subunit structures, which may interfere with subunit association. These subunit structures provided models for 12 rRNA expansion segments and 3 novel r-proteins. Furthermore, the 50S-aRF1 complex structure showed the unique domain orientation of aRF1, possibly explaining P-site transfer RNA (tRNA) release after translation termination. Sac 70S complexes were captured in seven distinct steps of the tRNA translocation reaction, confirming conserved structural features during archaeal ribosome translocation. In aEF2-engaged 70S ribosome complexes, 3D classification of cryo-EM data based on 30S head domain identified two new translocation intermediates with 30S head domain tilted 5-6° enabling its disengagement from the translocated tRNA and its release post-translocation. Additionally, we observed conformational changes to aEF2 during ribosome binding and switching from three different states. Our structural and biochemical data provide new insights into archaeal translation and ribosome translocation.

Understanding the role of a given transcription factor (TF) in regulating gene expression requires precise mapping of its binding sites in the genome. Chromatin immunoprecipitation-exo, an emerging technique using λ exonuclease to digest TF unbound DNA after ChIP, is designed to reveal transcription factor binding site (TFBS) boundaries with near-single nucleotide resolution. Although ChIP-exo promises deeper insights into transcription regulation, no dedicated bioinformatics tool exists to leverage its advantages. Most ChIP-seq and ChIP-chip analytic methods are not tailored for ChIP-exo, and thus cannot take full advantage of high-resolution ChIP-exo data. Here we describe a novel analysis framework, termed MACE (model-based analysis of ChIP-exo) dedicated to ChIP-exo data analysis. The MACE workflow consists of four steps: (i) sequencing data normalization and bias correction; (ii) signal consolidation and noise reduction; (iii) single-nucleotide resolution border peak detection using the Chebyshev Inequality and (iv) border matching using the Gale-Shapley stable matching algorithm. When applied to published human CTCF, yeast Reb1 and our own mouse ONECUT1/HNF6 ChIP-exo data, MACE is able to define TFBSs with high sensitivity, specificity and spatial resolution, as evidenced by multiple criteria including motif enrichment, sequence conservation, direct sequence pileup, nucleosome positioning and open chromatin states. In addition, we show that the fundamental advance of MACE is the identification of two boundaries of a TFBS with high resolution, whereas other methods only report a single location of the same event. The two boundaries help elucidate the in vivo binding structure of a given TF, e.g. whether the TF may bind as dimers or in a complex with other co-factors.

Human polynucleotide kinase/phosphatase (PNKP) is a dual specificity 5'-DNA kinase/3'-DNA phosphatase, with roles in base excision repair, DNA single-strand break repair and non-homologous end joining (NHEJ); yet precisely how PNKP functions in the repair of DNA double strand breaks (DSBs) remains unclear. We demonstrate that PNKP is phosphorylated by the DNA-dependent protein kinase (DNA-PK) and ataxia-telangiectasia mutated (ATM) in vitro. The major phosphorylation site for both kinases was serine 114, with serine 126 being a minor site. Ionizing radiation (IR)-induced phosphorylation of cellular PNKP on S114 was ATM dependent, whereas phosphorylation of PNKP on S126 required both ATM and DNA-PK. Inactivation of DNA-PK and/or ATM led to reduced PNKP at DNA damage sites in vivo. Cells expressing PNKP with alanine or aspartic acid at serines 114 and 126 were modestly radiosensitive and IR enhanced the association of PNKP with XRCC4 and DNA ligase IV; however, this interaction was not affected by mutation of PNKP phosphorylation sites. Purified PNKP protein with mutation of serines 114 and 126 had decreased DNA kinase and DNA phosphatase activities and reduced affinity for DNA in vitro. Together, our results reveal that IR-induced phosphorylation of PNKP by ATM and DNA-PK regulates PNKP function at DSBs.

Protein-protein interaction (PPI) network maintains proper function of all organisms. Simple high-throughput technologies are desperately needed to delineate the landscape of PPI networks. While recent state-of-the-art yeast two-hybrid (Y2H) systems improved screening efficiency, either individual colony isolation, library preparation arrays, gene barcoding or massive sequencing are still required. Here, we developed a recombination-based 'library vs library' Y2H system (RLL-Y2H), by which multi-library screening can be accomplished in a single pool without any individual treatment. This system is based on the phiC31 integrase-mediated integration between bait and prey plasmids. The integrated fragments were digested by MmeI and subjected to deep sequencing to decode the interaction matrix. We applied this system to decipher the trans-kingdom interactome between Mycobacterium tuberculosis and host cells and further identified Rv2427c interfering with the phagosome-lysosome fusion. This concept can also be applied to other systems to screen protein-RNA and protein-DNA interactions and delineate signaling landscape in cells.