Transcription Factor Binding Motifs Research Articles

ScMultiome analysis identifies a single caudal hindbrain compartment in the developing zebrafish nervous system

BackgroundA key step in nervous system development involves the coordinated control of neural progenitor specification and positioning. A long-standing model for the vertebrate CNS postulates that transient anatomical compartments – known as neuromeres – function to position neural progenitors along the embryonic anteroposterior neuraxis. Such neuromeres are apparent in the embryonic hindbrain – that contains six rhombomeres with morphologically apparent boundaries – but other neuromeres lack clear morphological boundaries and have instead been defined by different criteria, such as differences in gene expression patterns and the outcomes of transplantation experiments. Accordingly, the caudal hindbrain (CHB) posterior to rhombomere (r) 6 has been variably proposed to contain from two to five ‘pseudo-rhombomeres’, but the lack of comprehensive molecular data has precluded a detailed definition of such structures.MethodsWe used single-cell Multiome analysis, which allows simultaneous characterization of gene expression and chromatin state of individual cell nuclei, to identify and characterize CHB progenitors in the developing zebrafish CNS.ResultsWe identified CHB progenitors as a transcriptionally distinct population, that also possesses a unique profile of accessible transcription factor binding motifs, relative to both r6 and the spinal cord. This CHB population can be subdivided along its dorsoventral axis based on molecular characteristics, but we do not find any molecular evidence that it contains multiple pseudo-rhombomeres. We further observe that the CHB is closely related to r6 at the earliest embryonic stages, but becomes more divergent over time, and that it is defined by a unique gene regulatory network.ConclusionsWe conclude that the early CHB represents a single neuromere compartment that cannot be molecularly subdivided into pseudo-rhombomeres and that it may share an embryonic origin with r6.

The C. elegans Myc-family of transcription factors coordinate a dynamic adaptive response to dietary restriction.

Dietary restriction (DR), the process of decreasing overall food consumption over an extended period of time, has been shown to increase longevity across evolutionarily diverse species and delay the onset of age-associated diseases in humans. In Caenorhabditis elegans, the Myc-family transcription factors (TFs) MXL-2 (Mlx) and MML-1 (MondoA/ChREBP), which function as obligate heterodimers, and PHA-4 (orthologous to FOXA) are both necessary for the full physiological benefits of DR. However, the adaptive transcriptional response to DR and the role of MML-1::MXL-2 and PHA-4 remains elusive. We identified the transcriptional signature of C. elegans DR, using the eat-2 genetic model, and demonstrate broad changes in metabolic gene expression in eat-2 DR animals, which requires both mxl-2 and pha-4. While the requirement for these factors in DR gene expression overlaps, we found many of the DR genes exhibit an opposing change in relative gene expression in eat-2;mxl-2 animals compared to wild-type, which was not observed in eat-2 animals with pha-4 loss. Surprisingly, we discovered more than 2000 genes synthetically dysregulated in eat-2;mxl-2, out of which the promoters of down-regulated genes were substantially enriched for PQM-1 and ELT-1/3 GATA TF binding motifs. We further show functional deficiencies of the mxl-2 loss in DR outside of lifespan, as eat-2;mxl-2 animals exhibit substantially smaller brood sizes and lay a proportion of dead eggs, indicating that MML-1::MXL-2 has a role in maintaining the balance between resource allocation to the soma and to reproduction under conditions of chronic food scarcity. While eat-2 animals do not show a significantly different metabolic rate compared to wild-type, we also find that loss of mxl-2 in DR does not affect the rate of oxygen consumption in young animals. The gene expression signature of eat-2 mutant animals is consistent with optimization of energy utilization and resource allocation, rather than induction of canonical gene expression changes associated with acute metabolic stress, such as induction of autophagy after TORC1 inhibition. Consistently, eat-2 animals are not substantially resistant to stress, providing further support to the idea that chronic DR may benefit healthspan and lifespan through efficient use of limited resources rather than broad upregulation of stress responses, and also indicates that MML-1::MXL-2 and PHA-4 may have distinct roles in promotion of benefits in response to different pro-longevity stimuli.

Transfer learning reveals sequence determinants of the quantitative response to transcription factor dosage.

Deep learning approaches have made significant advances in predicting cell type-specific chromatin patterns from the identity and arrangement of transcription factor (TF) binding motifs. However, most models have been applied in unperturbed contexts, precluding a predictive understanding of how chromatin state responds to TF perturbation. Here, we used transfer learning to train and interpret deep learning models that use DNA sequence to predict, with accuracy approaching experimental reproducibility, how the concentration of two dosage-sensitive TFs (TWIST1, SOX9) affects regulatory element (RE) chromatin accessibility in facial progenitor cells. High-affinity motifs that allow for heterotypic TF co-binding and are concentrated at the center of REs buffer against quantitative changes in TF dosage and strongly predict unperturbed accessibility. In contrast, motifs with low-affinity or homotypic binding distributed throughout REs lead to sensitive responses with minimal contributions to unperturbed accessibility. Both buffering and sensitizing features show signatures of purifying selection. We validated these predictive sequence features using reporter assays and showed that a biophysical model of TF-nucleosome competition can explain the sensitizing effect of low-affinity motifs. Our approach of combining transfer learning and quantitative measurements of the chromatin response to TF dosage therefore represents a powerful method to reveal additional layers of the cis-regulatory code.

Molecular regulation of PPARγ/RXRα signaling by the novel cofactor ZFP407.

Cofactors interacting with PPARγ can regulate adipogenesis and adipocyte metabolism by modulating the transcriptional activity and selectivity of PPARγ signaling. ZFP407 was previously demonstrated to regulate PPARγ target genes such as GLUT4, and its overexpression improved glucose homeostasis in mice. Here, using a series of molecular assays, including protein-interaction studies, mutagenesis, and ChIP-seq, ZFP407 was found to interact with the PPARγ/RXRα protein complex in the nucleus of adipocytes. Consistent with this observation, ZFP407 ChIP-seq peaks significantly overlapped with PPARγ ChIP-seq peaks, with more than half of ZFP407 peaks overlapping with PPARγ peaks. Transcription factor binding motifs enriched in these overlapping sites included CTCF, RARα/RXRγ, TP73, and ELK1, which regulate cellular development and function within adipocytes. Site-directed mutagenesis of frequent PPARγ phosphorylation or SUMOylation sites did not prevent its regulation by ZFP407, while mutagenesis of ZFP407 domains potentially necessary for RXR and PPARγ binding abrogated any impact of ZFP407 on PPARγ activity. These data suggest that ZFP407 controls the activity of PPARγ, but does so independently of post-translational modifications, likely by direct binding, establishing ZFP407 as a newly identified PPARγ cofactor. In addition, ZFP407 ChIP-seq analyses identified regions that did not overlap with PPARγ peaks. These non-overlapping peaks were significantly enriched for the transcription factor binding motifs of TBX19, PAX8, HSF4, and ZKSCAN3, which may contribute to the PPARγ-independent functions of ZFP407 in adipocytes and other cell types.

Single-cell analyses reveal transient retinal progenitor cells in the ciliary margin of developing human retina

The emergence of retinal progenitor cells and differentiation to various retinal cell types represent fundamental processes during retinal development. Herein, we provide a comprehensive single cell characterisation of transcriptional and chromatin accessibility changes that underline retinal progenitor cell specification and differentiation over the course of human retinal development up to midgestation. Our lineage trajectory data demonstrate the presence of early retinal progenitors, which transit to late, and further to transient neurogenic progenitors, that give rise to all the retinal neurons. Combining single cell RNA-Seq with spatial transcriptomics of early eye samples, we demonstrate the transient presence of early retinal progenitors in the ciliary margin zone with decreasing occurrence from 8 post-conception week of human development. In retinal progenitor cells, we identified a significant enrichment for transcriptional enhanced associate domain transcription factor binding motifs, which when inhibited led to loss of cycling progenitors and retinal identity in pluripotent stem cell derived organoids.

Transcription of microRNAs is regulated by developmental signaling pathways and transcription factors.

In early embryonic development, the cross-regulation of transcription factors and signaling pathways are critical in mediating developmental and physiological processes. Additionally, many studies have shown the importance of post-transcriptional regulation of signaling and network components mediated by microRNAs (miRNAs); however, how miRNAs are transcriptionally regulated is poorly understood. miRNAs are critical fine-tuners of many biological processes and their dysregulation leads to a variety of diseases and developmental defects. Previously, we have shown that miRNAs are dynamically expressed throughout sea urchin development, suggesting that miRNAs are likely to be under transcriptional regulation. Here, we used pharmacological inhibitors, genetic constructs, and loss-of-function reagents to assess the impact of key signaling pathways (Wnt, Nodal, MAPK, Sonic Hedgehog, Delta/Notch, VEGF, and BMP) and transcription factors (Alx1, Ets1/2, and Tbr) on the transcript levels of the evolutionarily conserved miR-1, miR-31, miR-92 and miR-124; the invertebrate-specific miR-71; and the echinoderm-specific miR-2002, miR-2007, and miR-2012. We also used computational methods to identify potential transcription factor binding sites of these miRNAs. Lists of binding motifs for transcription factors (TFs) were acquired from the MEME-Suite Motif Database and used as inputs for the algorithm FIMO (Find Individual Motif Occurrences), which detects short nucleotide motifs within larger sequences. Based on experimental data on miRNA expression in conjunction with bioinformatic predictions, we propose that the transcription factors Tbr, Alx1, and Ets1 regulate SpmiR-1, SpmiR-31, and SpmiR-71, respectively. We additionally observed significant effects on miRNA levels as a result of perturbations to Wnt, Nodal, MAPK, and Sonic Hedgehog signaling pathways, while no significant change on miRNA levels were observed with perturbations to Delta/Notch, VEGF, or BMP signaling pathways. Overall, this study provides insights into the transcriptional regulation of miRNAs by signaling pathways and transcription factors and contribute to our overall understanding of the genetic regulation of developmental processes.

Type I interferons induce an epigenetically distinct memory B cell subset in chronic viral infection

Memory B cells (MBCs) are key providers of long-lived immunity against infectious disease, yet in chronic viral infection, they do not produce effective protection. How chronic viral infection disrupts MBC development and whether such changes are reversible remain unknown. Through single-cell (sc)ATAC-seq and scRNA-seq during acute versus chronic lymphocytic choriomeningitis viral infection, we identified a memory subset enriched for interferon (IFN)-stimulated genes (ISGs) during chronic infection that was distinct from the T-bet+ subset normally associated with chronic infection. Blockade of IFNAR-1 early in infection transformed the chromatin landscape of chronic MBCs, decreasing accessibility at ISG-inducing transcription factor binding motifs and inducing phenotypic changes in the dominating MBC subset, with a decrease in the ISG subset and an increase in CD11c+CD80+ cells. However, timing was critical, with MBCs resistant to intervention at 4weeks post-infection. Together, our research identifies a key mechanism to instruct MBC identity during viral infection.

Abstract LB293: O-GlcNAc transferase drives the proliferative state in adaptive drug tolerance

Abstract Drivers of the drug tolerant proliferative persister (DTPP) state have not been well investigated. Histone H3 lysine-4 trimethylation (H3K4me3), an active histone mark, might enable slow cycling drug tolerant persisters (DTP) to regain proliferative capacity. This study aimed to determine H3K4me3 transcriptionally active sites identifying a key regulator of DTPPs in melanoma. Deploying a model of adaptive cancer drug tolerance H3K4me3 ChIP-seq data of DTPPs guided identification of top transcription factor binding motifs. These suggested involvement of O-linked N-acetylglucosamine transferase (OGT), which was confirmed by metabolomics analysis and biochemical assays. OGT impact on DTPPs and adaptive resistance was explored in vitro and in vivo. H3K4me3 remodeling was widespread in CPG island regions and DNA binding motifs associated with O-GlcNAc marked chromatin. Accordingly, we observed an upregulation of OGT, O-GlcNAc and its binding partner TET1 in chronically treated cancer cells. Inhibition of OGT led to loss of H3K4me3 and downregulation of genes contributing to drug resistance. Genetic ablation of OGT prevented acquired drug resistance in vivo. Our findings uncover a fundamental mechanism of adaptive drug resistance that governs cancer cell reprogramming towards acquired drug resistance, a process that can be exploited to improve response duration and patient outcomes. Citation Format: Dinoop Ravindran Menon, Heinz Hammerlindl, Gregory Gimenez, Elmar Zuegner, Christoph Magens, Michael Eccles, Mayumi Fujita, Helmut Schaider. O-GlcNAc transferase drives the proliferative state in adaptive drug tolerance [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB293.

The Putative GATA Transcription Factor SbGATA22 as a Novel Regulator of Dhurrin Biosynthesis.

Cyanogenic glucosides are specialized metabolites produced by over 3000 species of higher plants from more than 130 families. The deployment of cyanogenic glucosides is influenced by biotic and abiotic factors in addition to being developmentally regulated, consistent with their roles in plant defense and stress mitigation. Despite their ubiquity, very little is known regarding the molecular mechanisms that regulate their biosynthesis. The biosynthetic pathway of dhurrin, the cyanogenic glucoside found in the important cereal crop sorghum (Sorghum bicolor (L.) Moench), was described over 20 years ago, and yet no direct regulator of the biosynthetic genes has been identified. To isolate regulatory proteins that bind to the promoter region of the key dhurrin biosynthetic gene of sorghum, SbCYP79A1, yeast one-hybrid screens were performed. A bait fragment containing 1204 base pairs of the SbCYP79A1 5' regulatory region was cloned upstream of a reporter gene and introduced into Saccharomyces cerevisiae. Subsequently, the yeast was transformed with library cDNA representing RNA from two different sorghum developmental stages. From these screens, we identified SbGATA22, an LLM domain B-GATA transcription factor that binds to the putative GATA transcription factor binding motifs in the SbCYP79A1 promoter region. Transient assays in Nicotiana benthamiana show that SbGATA22 localizes to the nucleus. The expression of SbGATA22, in comparison with SbCYP79A1 expression and dhurrin concentration, was analyzed over 14 days of sorghum development and in response to nitrogen application, as these conditions are known to affect dhurrin levels. Collectively, these findings suggest that SbGATA22 may act as a negative regulator of SbCYP79A1 expression and provide a preliminary insight into the molecular regulation of dhurrin biosynthesis in sorghum.

Wemics: A Single-Base Resolution Methylation Quantification Method for Enhanced Prediction of Epigenetic Regulation.

DNA methylation, an epigenetic mechanism that alters gene expression without changing DNA sequence, is essential for organism development and key biological processes like genomic imprinting and X-chromosome inactivation. Despite tremendous efforts in DNA methylation research, accurate quantification of cytosine methylation remains a challenge. Here, a single-base methylation quantification approach is introduced by weighting methylation of consecutive CpG sites (Wemics) in genomic regions. Wemics quantification of DNA methylation better predicts its regulatory impact on gene transcription and identifies differentially methylated regions (DMRs) with more biological relevance. Most Wemics-quantified DMRs in lung cancer are epigenetically conserved and recurrently occurred in other primary cancers from The Cancer Genome Atlas (TCGA), and their aberrant alterations can serve as promising pan-cancer diagnostic markers. It is further revealed that these detected DMRs are enriched in transcription factor (TF) binding motifs, and methylation of these TF binding motifs and TF expression synergistically regulate target gene expression. Using Wemics on epigenomic-transcriptomic data from the large lung cancer cohort, a dozen novel genes with oncogenic potential are discovered that are upregulated by hypomethylation but overlooked by other quantification methods. These findings increase the understanding of the epigenetic mechanism by which DNA methylation regulates gene expression.

Cardiac GR Mediates the Diurnal Rhythm in Ventricular Arrhythmia Susceptibility.

Ventricular arrhythmias (VAs) demonstrate a prominent day-night rhythm, commonly presenting in the morning. Transcriptional rhythms in cardiac ion channels accompany this phenomenon, but their role in the morning vulnerability to VAs and the underlying mechanisms are not understood. We investigated the recruitment of transcription factors that underpins transcriptional rhythms in ion channels and assessed whether this mechanism was pertinent to the heart's intrinsic diurnal susceptibility to VA. Assay for transposase-accessible chromatin with sequencing performed in mouse ventricular myocyte nuclei at the beginning of the animals' inactive (ZT0) and active (ZT12) periods revealed differentially accessible chromatin sites annotating to rhythmically transcribed ion channels and distinct transcription factor binding motifs in these regions. Notably, motif enrichment for the glucocorticoid receptor (GR; transcriptional effector of corticosteroid signaling) in open chromatin profiles at ZT12 was observed, in line with the well-recognized ZT12 peak in circulating corticosteroids. Molecular, electrophysiological, and in silico biophysically-detailed modeling approaches demonstrated GR-mediated transcriptional control of ion channels (including Scn5a underlying the cardiac Na+ current, Kcnh2 underlying the rapid delayed rectifier K+ current, and Gja1 responsible for electrical coupling) and their contribution to the day-night rhythm in the vulnerability to VA. Strikingly, both pharmacological block of GR and cardiomyocyte-specific genetic knockout of GR blunted or abolished ion channel expression rhythms and abolished the ZT12 susceptibility to pacing-induced VA in isolated hearts. Our study registers a day-night rhythm in chromatin accessibility that accompanies diurnal cycles in ventricular myocytes. Our approaches directly implicate the cardiac GR in the myocyte excitability rhythm and mechanistically link the ZT12 surge in glucocorticoids to intrinsic VA propensity at this time.

Abstract 2342: BARTsc: A transcription factor analysis suite for single-cell omics data

Abstract In Eukaryotic cells, the expression of a small set of transcription factors (TFs) is usually sufficient to set up a cell-type-specific gene expression program. The underlying complicated epigenetic mechanisms that explain how cells change their states in response to intra- and extra-cellular signals in different scenarios remain to be fully understood. Importantly, epigenetic regulation of gene expression has been considered a central regulatory mechanism of cell fate determination and cellular plasticity. Insights into it will also profoundly impact our understanding of misregulation of gene expression in diseases such as cancer. However, the functional TFs of many cell types remain unclear, which greatly hinders the study of epigenetic regulatory mechanisms. The rapid development of single-cell technology is a promising opportunity that enables us to infer functional TFs for different cell types. Some computational methods have been published in this area. However, existing methods have various pitfalls. For instance, some methods use transcription factor binding motifs (TFBMs) as the reference of TF binding, overlooking widespread non-specific binding. Some other methods hold strict criteria for training data and thus can only predict a small number of TFs. To overcome these issues, we propose to leverage publicly available high-quality ChIP-seq data to predict TFs functional in single-cell omics data. Here, we present BARTsc, a comprehensive suite for a series of TF analysis functionalities, including inference of signature TFs of given cell types/clusters, cross-cell-type/cluster TF activity comparative analysis, and TF clustering based on their cross cell-type distributions. The input of BARTsc can be either scRNA-seq, scATAC-seq, or scMultiome data. We demonstrate that BARTsc outperforms several existing tools in both sensitivity and specificity in identifying key cell type regulators in single-cell datasets of well-studied systems. Applying BARTsc to a scRNA-seq dataset for cancer-associated fibroblasts (CAFs), we identified distinct regulatory activities of some TFs in different subtypes of CAFs, but the TF genes are universally expressed in these subtypes. In addition, by time course BARTsc analysis, we identified a list of TFs that may either promote or suppress CAF differentiation, offering a new approach to find potential candidate therapeutic targets. Citation Format: Hongpan Zhang, Jingyi Wang, Zhenjia Wang, Chongzhi Zang. BARTsc: A transcription factor analysis suite for single-cell omics data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2342.

Abstract 5853: Regulatory network inference of lung adenocarcinoma drug-tolerant persister cells

Abstract Despite the effectiveness of the latest generation of cancer drugs, partial tumor response and subsequent disease relapse is commonly observed in patients. Recent evidence suggests that treatment resistance emerges from drug-tolerant persister cells (DTPC) that survive through non-genetic tumor cell adaptation. Although DTPC are mostly slow-cycling, a rare subset can re-enter the cell cycle under constitutive presence of treatment. Understanding how DTPC dynamically remodel their molecular portraits at the single-cell level may provide high-resolution insights into the molecular circuits associated with adaptive survival to therapy.Here, we use single-cell RNA sequencing (scRNA-seq) to identify potential mechanisms that promote survival of DTPC and decipher their underlying gene regulatory networks (GRNs), in a model of EGFR-mutant lung adenocarcinoma treated with osimertinib (EGFR inhibitor). Gene expression analysis revealed fewer cell states in DTP compared to treatment-naïve cells (TN). We confirmed the acquisition of DTPC hallmarks (e.g., high epithelial-mesenchymal transition, slow-cycling phenotype, deficient DNA repair). We also identified subpopulations of cells that display intermediate transcriptional states: a pre-existing DTP-like state in a bulk of TN cells, and DTPC with enriched cell cycle signatures, evoking a phenotype of cycling DTPC. To characterize the GRN that govern TN cells, DTPC and the intermediate cell states, we applied the network inference tool PYSCENIC, which associates transcription factors (TFs) to their candidate target genes (TGs) on the basis of the scRNA-seq data and enrichment of TF binding motifs in the genomic regions around TGs. The resulting (GRN) consisted of 5 TF modules clustered on the basis of correlated TF activity.This analysis revealed that, compared to TN cells, DTPC and the intermediate cell states display a marked decreased expression of TFs and TGs involved in mRNA splicing, and more specifically in the initiation of spliceosome assembly through E-complex and A-complex formation. Interestingly, the splicing-related TFs all appear to have affinity for binding to a GCCxC consensus motif or close variant in genes (i.e., SNRPB and SNRPE) that code for key effectors of these initiation complexes. This pilot study highlights that scRNA-seq analyses coupled with inference of regulatory networks might enable the identification of molecular mechanisms that potentially trigger cell state transitions during drug-tolerance. The role of altered mRNA splicing in persistence to EGFR-targeted therapies will be explored experimentally in patient-derived cell lines and 3D organoid cultures. Citation Format: Lisa Chabrier, Camille Leonce, Raphael Schneider, Cyril Degletagne, Pierre Saintigny, Anton Crombach, Sandra Ortiz-Cuaran. Regulatory network inference of lung adenocarcinoma drug-tolerant persister cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 5853.

Endothelial gene regulatory elements associated with cardiopharyngeal lineage differentiation

Endothelial cells (EC) differentiate from multiple sources, including the cardiopharyngeal mesoderm, which gives rise also to cardiac and branchiomeric muscles. The enhancers activated during endothelial differentiation within the cardiopharyngeal mesoderm are not completely known. Here, we use a cardiogenic mesoderm differentiation model that activates an endothelial transcription program to identify endothelial regulatory elements activated in early cardiogenic mesoderm. Integrating chromatin remodeling and gene expression data with available single-cell RNA-seq data from mouse embryos, we identify 101 putative regulatory elements of EC genes. We then apply a machine-learning strategy, trained on validated enhancers, to predict enhancers. Using this computational assay, we determine that 50% of these sequences are likely enhancers, some of which are already reported. We also identify a smaller set of regulatory elements of well-known EC genes and validate them using genetic and epigenetic perturbation. Finally, we integrate multiple data sources and computational tools to search for transcriptional factor binding motifs. In conclusion, we show EC regulatory sequences with a high likelihood to be enhancers, and we validate a subset of them using computational and cell culture models. Motif analyses show that the core EC transcription factors GATA/ETS/FOS is a likely driver of EC regulation in cardiopharyngeal mesoderm.

Epigenetic MLH1 silencing concurs with mismatch repair deficiency in sporadic, naturally occurring colorectal cancer in rhesus macaques

BackgroundNaturally occurring colorectal cancers (CRC) in rhesus macaques share many features with their human counterparts and are useful models for cancer immunotherapy; but mechanistic data are lacking regarding the comparative molecular pathogenesis of these cancers.MethodsWe conducted state-of-the-art imaging including CT and PET, clinical assessments, and pathological review of 24 rhesus macaques with naturally occurring CRC. Additionally, we molecularly characterized these tumors utilizing immunohistochemistry (IHC), microsatellite instability assays, DNAseq, transcriptomics, and developed a DNA methylation-specific qPCR assay for MLH1, CACNA1G, CDKN2A, CRABP1, and NEUROG1, human markers for CpG island methylator phenotype (CIMP). We furthermore employed Monte-Carlo simulations to in-silico model alterations in DNA topology in transcription-factor binding site-rich promoter regions upon experimentally demonstrated DNA methylation.ResultsSimilar cancer histology, progression patterns, and co-morbidities could be observed in rhesus as reported for human CRC patients. IHC identified loss of MLH1 and PMS2 in all cases, with functional microsatellite instability. DNA sequencing revealed the close genetic relatedness to human CRCs, including a similar mutational signature, chromosomal instability, and functionally-relevant mutations affecting KRAS (G12D), TP53 (R175H, R273*), APC, AMER1, ALK, and ARID1A. Interestingly, MLH1 mutations were rarely identified on a somatic or germline level. Transcriptomics not only corroborated the similarities of rhesus and human CRCs, but also demonstrated the significant downregulation of MLH1 but not MSH2, MSH6, or PMS2 in rhesus CRCs. Methylation-specific qPCR suggested CIMP-positivity in 9/16 rhesus CRCs, but all 16/16 exhibited significant MLH1 promoter hypermethylation. DNA hypermethylation was modelled to affect DNA topology, particularly propeller twist and roll profiles. Modelling the DNA topology of a transcription factor binding motif (TFAP2A) in the MLH1 promoter that overlapped with a methylation-specific probe, we observed significant differences in DNA topology upon experimentally shown DNA methylation. This suggests a role of transcription factor binding interference in epigenetic silencing of MLH1 in rhesus CRCs.ConclusionsThese data indicate that epigenetic silencing suppresses MLH1 transcription, induces the loss of MLH1 protein, abrogates mismatch repair, and drives genomic instability in naturally occurring CRC in rhesus macaques. We consider this spontaneous, uninduced CRC in immunocompetent, treatment-naïve rhesus macaques to be a uniquely informative model for human CRC.Graphical abstract

Unsupervised AI reveals insect species-specific genome signatures.

Insects are a highly diverse phylogeny and possess a wide variety of traits, including the presence or absence of wings and metamorphosis. These diverse traits are of great interest for studying genome evolution, and numerous comparative genomic studies have examined a wide phylogenetic range of insects. Here, we analyzed 22 insects belonging to a wide phylogenetic range (Endopterygota, Paraneoptera, Polyneoptera, Palaeoptera, and other insects) by using a batch-learning self-organizing map (BLSOM) for oligonucleotide compositions in their genomic fragments (100-kb or 1-Mb sequences), which is an unsupervised machine learning algorithm that can extract species-specific characteristics of the oligonucleotide compositions (genome signatures). The genome signature is of particular interest in terms of the mechanisms and biological significance that have caused the species-specific difference, and can be used as a powerful search needle to explore the various roles of genome sequences other than protein coding, and can be used to unveil mysteries hidden in the genome sequence. Since BLSOM is an unsupervised clustering method, the clustering of sequences was performed based on the oligonucleotide composition alone, without providing information about the species from which each fragment sequence was derived. Therefore, not only the interspecies separation, but also the intraspecies separation can be achieved. Here, we have revealed the specific genomic regions with oligonucleotide compositions distinct from the usual sequences of each insect genome, e.g., Mb-level structures found for a grasshopper Schistocerca americana. One aim of this study was to compare the genome characteristics of insects with those of vertebrates, especially humans, which are phylogenetically distant from insects. Recently, humans seem to be the "model organism" for which a large amount of information has been accumulated using a variety of cutting-edge and high-throughput technologies. Therefore, it is reasonable to use the abundant information from humans to study insect lineages. The specific regions of Mb length with distinct oligonucleotide compositions have also been previously observed in the human genome. These regions were enriched by transcription factor binding motifs (TFBSs) and hypothesized to be involved in the three-dimensional arrangement of chromosomal DNA in interphase nuclei. The present study characterized the species-specific oligonucleotide compositions (i.e., genome signatures) in insect genomes and identified specific genomic regions with distinct oligonucleotide compositions.

Genome evolution and divergence in cis-regulatory architecture is associated with condition-responsive development in horned dung beetles.

Phenotypic plasticity is thought to be an important driver of diversification and adaptation to environmental variation, yet the genomic mechanisms mediating plastic trait development and evolution remain poorly understood. The Scarabaeinae, or true dung beetles, are a species-rich clade of insects recognized for their highly diversified nutrition-responsive development including that of cephalic horns-evolutionarily novel, secondary sexual weapons that exhibit remarkable intra- and interspecific variation. Here, we investigate the evolutionary basis for horns as well as other key dung beetle traits via comparative genomic and developmental assays. We begin by presenting chromosome-level genome assemblies of three dung beetle species in the tribe Onthophagini (> 2500 extant species) including Onthophagus taurus, O. sagittarius, and Digitonthophagus gazella. Comparing these assemblies to those of seven other species across the order Coleoptera identifies evolutionary changes in coding sequence associated with metabolic regulation of plasticity and metamorphosis. We then contrast chromatin accessibility in developing head horn tissues of high- and low-nutrition O. taurus males and females and identify distinct cis-regulatory architectures underlying nutrition- compared to sex-responsive development, including a large proportion of recently evolved regulatory elements sensitive to horn morph determination. Binding motifs of known and new candidate transcription factors are enriched in these nutrition-responsive open chromatin regions. Our work highlights the importance of chromatin state regulation in mediating the development and evolution of plastic traits, demonstrates gene networks are highly evolvable transducers of environmental and genetic signals, and provides new reference-quality genomes for three species that will bolster future developmental, ecological, and evolutionary studies of this insect group.

Deciphering the spatiotemporal transcriptional and chromatin accessibility of human retinal organoid development at the single-cell level

Molecular information on the early stages of human retinal development remains scarce due to limitations in obtaining early human eye samples. Pluripotent stem cell-derived retinal organoids (ROs) provide an unprecedented opportunity for studying early retinogenesis. Using a combination of single cell RNA-seq and spatial transcriptomics we present for the first-time a single cell spatiotemporal transcriptome of RO development. Our data demonstrate that ROs recapitulate key events of retinogenesis including optic vesicle/cup formation, presence of a putative ciliary margin zone, emergence of retinal progenitor cells and their orderly differentiation to retinal neurons. Combining the scRNA- with scATAC-seq data, we were able to reveal cell-type specific transcription factor binding motifs on accessible chromatin at each stage of organoid development, and to show that chromatin accessibility is highly correlated to the developing human retina, but with some differences in the temporal emergence and abundance of some of the retinal neurons.

SpecLoop predicts cell type-specific chromatin loop via transcription factor cooperation

Cell-type-Specific Chromatin Loops (CSCLs) are crucial for gene regulation and cell fate determination. However, the mechanisms governing their establishment remain elusive. Here, we present SpecLoop, a network regularization-based machine learning framework, to investigate the role of transcription factors (TFs) cooperation in CSCL formation. SpecLoop integrates multi-omics data, including gene expression, chromatin accessibility, sequence, protein-protein interaction, and TF binding motif data, to predict CSCLs and identify TF cooperations. Using high resolution Hi-C data as the gold standard, SpecLoop accurately predicts CSCL in GM12878, IMR90, HeLa–S3, K562, HUVEC, HMEC, and NHEK seven cell types, with the AUROC values ranging from 0.8645 to 0.9852 and AUPR values ranging from 0.8654 to 0.9734. Notably SpecLoop demonstrates improved accuracy in predicting long-distance CSCLs and identifies TF complexes with strong predictive ability. Our study systematically explores the TFs and TF pairs associated with CSCL through effective integration of diverse omics data. SpecLoop is freely available at https://github.com/AMSSwanglab/SpecLoop.

Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4.

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here, we determine the roles of CHD4 inenhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility. Its depletion leads to redistribution of transcription factors to previously unoccupied sites. During cellular reprogramming induced by the pioneer factor GATA3, CHD4 activity is necessary to prevent inappropriate chromatin opening. Mechanistically, CHD4 promotes nucleosome positioning over GATA3 binding motifs to compete with transcription factor-DNA interaction. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents unnecessary gene expression by editing chromatin binding activities of transcription factors.