An atlas of combinatorial transcriptional regulation in mouse and man

Single-cell analysis of transcription factor regulatory networks reveals molecular basis for subtype-specific dysregulation in acute myeloid leukemia

Highly heterogeneous acute myeloid leukemia (AML) exhibits dysregulated transcriptional programs. Transcription factor (TF) regulatory networks underlying AML subtypes have not been elucidated at single-cell resolution. Here, we comprehensively mapped malignancy-related TFs activated in different AML subtypes by analyzing single-cell RNA sequencing data from AMLs and healthy donors. We first identified six modules of regulatory networks which were prevalently dysregulated in all AML patients. AML subtypes featured with different malignant cellular composition possessed subtype-specific regulatory TFs associated with differentiation suppression or immune modulation. At last, we validated that ERF was crucial for the development of hematopoietic stem/progenitor cells by performing loss- and gain-of-function experiments in zebrafish embryos. Collectively, our work thoroughly documents an abnormal spectrum of transcriptional regulatory networks in AML and reveals subtype-specific dysregulation basis, which provides a prospective view to AML pathogenesis and potential targets for both diagnosis and therapy.

The BTB transcription factors ZBTB11 and ZFP131 maintain pluripotency by repressing pro-differentiation genes

In pluripotent cells, a delicate activation-repression balance maintains pro-differentiation genes ready for rapid activation. The identity of transcription factors (TFs) that specifically repress pro-differentiation genes remains obscure. By targeting ∼1,700 TFs with CRISPR loss-of-function screen, we found that ZBTB11 and ZFP131 are required for embryonic stem cell (ESC) pluripotency. ESCs without ZBTB11 or ZFP131 lose colony morphology, reduce proliferation rate, and upregulate transcription of genes associated with three germ layers. ZBTB11 and ZFP131 bind proximally to pro-differentiation genes. ZBTB11 or ZFP131 loss leads to an increase in H3K4me3, negative elongation factor (NELF) complex release, and concomitant transcription at associated genes. Together, our results suggest that ZBTB11 and ZFP131 maintain pluripotency by preventing premature expression of pro-differentiation genes and present a generalizable framework to maintain cellular potency.

The Fra-1/AP-1 Oncoprotein: From the "Undruggable" Transcription Factor to Therapeutic Targeting

The genetic and epigenetic changes affecting transcription factors, coactivators, and chromatin modifiers are key determinants of the hallmarks of cancer. The acquired dependence on oncogenic transcriptional regulators, representing a major determinant of cancer cell vulnerability, points to transcription factors as ideal therapeutic targets. However, given the unavailability of catalytic activities or binding pockets for small-molecule inhibitors, transcription factors are generally regarded as undruggable proteins. Among components of the AP-1 complex, the FOS-family transcription factor Fra-1, encoded by FOSL1, has emerged as a prominent therapeutic target. Fra-1 is overexpressed in most solid tumors, in response to the BRAF-MAPK, Wnt-beta-catenin, Hippo-YAP, IL-6-Stat3, and other major oncogenic pathways. In vitro functional analyses, validated in onco-mouse models and corroborated by prognostic correlations, show that Fra-1-containing dimers control tumor growth and disease progression. Fra-1 participates in key mechanisms of cancer cell invasion, Epithelial-to-Mesenchymal Transition, and metastatic spreading, by driving the expression of EMT-inducing transcription factors, cytokines, and microRNAs. Here we survey various strategies aimed at inhibiting tumor growth, metastatic dissemination, and drug resistance by interfering with Fra-1 expression, stability, and transcriptional activity. We summarize several tools aimed at the design and tumor-specific delivery of Fra-1/AP-1-specific drugs. Along with RNA-based therapeutics targeting the FOSL1 gene, its mRNA, or cognate regulatory circRNAs, we will examine the exploitation of blocking peptides, small molecule inhibitors, and innovative Fra-1 protein degraders. We also consider the possible caveats concerning Fra-1 inhibition in specific therapeutic contexts. Finally, we discuss a recent suicide gene therapy-based approach, aimed at selectively killing the Fra-1-overexpressing neoplastic cells.

Are Transcription Factors Plausible Oncotargets for Triple Negative Breast Cancers?

Simple Summary

Triple negative breast cancer is a type of breast cancer that does not have a selective and effective therapy. It is known that this cancer possesses high abundance of certain proteins called transcription factors, which are essential for their growth. However, inhibiting transcription factors is very difficult with common therapeutics due to their inaccessibility inside the cell and their molecular structure. In this work, we identified the most important transcription factors for the growth of triple negative breast cancers, and that can predict worse clinical outcome. Moreover, we described different strategies that have been utilised to inhibit them. A successful inhibition of these transcription factors could reduce the mortality and convalescence associated with triple negative breast cancers.

Abstract

Breast cancer (BC) is the most diagnosed cancer worldwide and one of the main causes of cancer deaths. BC is a heterogeneous disease composed of different BC intrinsic subtypes such as triple-negative BC (TNBC), which is one of the most aggressive subtypes and which lacks a targeted therapy. Recent comprehensive analyses across cell types and cancer types have outlined a vast network of protein–protein associations between transcription factors (TFs). Not surprisingly, protein–protein networks central to oncogenesis and disease progression are highly altered during TNBC pathogenesis and are responsible for the activation of oncogenic programs, such as uncontrollable proliferation, epithelial-to-mesenchymal transition (EMT) and stemness. From the therapeutic viewpoint, inhibiting the interactions between TFs represents a very significant challenge, as the contact surfaces of TFs are relatively large and featureless. However, promising tools have emerged to offer a solution to the targeting problem. At the clinical level, some TF possess diagnostic and prognostic value in TNBC. In this review, we outline the recent advances in TFs relevant to TNBC growth and progression. Moreover, we highlight different targeting approaches to inhibit these TFs. Furthermore, the validity of such TFs as clinical biomarkers has been explored. Finally, we discuss how research is likely to evolve in the field.

From shallow to deep: some lessons learned from application of machine learning for recognition of functional genomic elements in human genome

Identification of genomic signals as indicators for functional genomic elements is one of the areas that received early and widespread application of machine learning methods. With time, the methods applied grew in variety and generally exhibited a tendency to improve their ability to identify some major genomic and transcriptomics signals. The evolution of machine learning in genomics followed a similar path to applications of machine learning in other fields. These were impacted in a major way by three dominant developments, namely an enormous increase in availability and quality of data, a significant increase in computational power available to machine learning applications, and finally, new machine learning paradigms, of which deep learning is the most well-known example. It is not easy in general to distinguish factors leading to improvements in results of applications of machine learning. This is even more so in the field of genomics, where the advent of next-generation sequencing and the increased ability to perform functional analysis of raw data have had a major effect on the applicability of machine learning in OMICS fields. In this paper, we survey the results from a subset of published work in application of machine learning in the recognition of genomic signals and regions in human genome and summarize some lessons learnt from this endeavor. There is no doubt that a significant progress has been made both in terms of accuracy and reliability of models. Questions remain however whether the progress has been sufficient and what these developments bring to the field of genomics in general and human genomics in particular. Improving usability, interpretability and accuracy of models remains an important open challenge for current and future research in application of machine learning and more generally of artificial intelligence methods in genomics.

Identification of ZBTB26 as a Novel Risk Factor for Congenital Hypothyroidism

Congenital primary hypothyroidism (CH; OMIM 218700) is characterized by an impaired thyroid development, or dyshormonogenesis, and can lead to intellectual disability and growth retardation if untreated. Most of the children with congenital hypothyroidism present thyroid dysgenesis, a developmental anomaly of the thyroid. Various genes have been associated with thyroid dysgenesis, but all known genes together can only explain a small number of cases. To identify novel genetic causes for congenital hypothyroidism, we performed trio whole-exome sequencing in an affected newborn and his unaffected parents. A predicted damaging de novo missense mutation was identified in the ZBTB26 gene (Zinc Finger A and BTB Domain containing 26). An additional cohort screening of 156 individuals with congenital thyroid dysgenesis identified two additional ZBTB26 gene variants of unknown significance. To study the underlying disease mechanism, morpholino knock-down of zbtb26 in Xenopus laevis was carried out, which demonstrated significantly smaller thyroid anlagen in knock-down animals at tadpole stage. Marker genes expressed in thyroid tissue precursors also indicated a specific reduction in the Xenopus ortholog of human Paired-Box-Protein PAX8, a transcription factor required for thyroid development, which could be rescued by adding zbtb26. Pathway and network analysis indicated network links of ZBTB26 to PAX8 and other genes involved in thyroid genesis and function. GWAS associations of ZBTB26 were found with height. Together, our study added a novel genetic risk factor to the list of genes underlying congenital primary hypothyroidism and provides additional support that de novo mutations, together with inherited variants, might contribute to the genetic susceptibility to CH.

Correlation of Expression Changes between Genes Controlling 5-HT Synthesis and Genes Crh and Trh in the Midbrain Raphe Nuclei of Chronically Aggressive and Defeated Male Mice

Midbrain raphe nuclei (MRNs) contain a large number of serotonergic neurons associated with the regulation of numerous types of psychoemotional states and physiological processes. The aim of this work was to study alterations of the MRN transcriptome in mice with prolonged positive or negative fighting experience and to identify key gene networks associated with the regulation of serotonergic system functioning. Numerous genes underwent alterations of transcription in the MRNs of male mice that either manifested aggression or experienced social defeat in daily agonistic interactions. The expression of the Tph2 gene encoding the rate-limiting enzyme of the serotonin synthesis pathway correlated with the expression of many genes, 31 of which were common between aggressive and defeated mice and were downregulated in the MRNs of mice of both experimental groups. Among these common differentially expressed genes (DEGs), there were genes associated with behavior, learning, memory, and synaptic signaling. These results suggested that, in the MRNs of the mice, the transcriptome changes associated with serotonergic regulation of various processes are similar between the two groups (aggressive and defeated). In the MRNs, more DEGs correlating with Tph2 expression were found in defeated mice than in the winners, which is probably a consequence of deeper Tph2 downregulation in the losers. It was shown for the first time that, in both groups of experimental mice, the changes in the transcription of genes controlling the synthesis and transport of serotonin directly correlate with the expression of genes Crh and Trh, which control the synthesis of corticotrophin- and thyrotropin-releasing hormones. Our findings indicate that CRH and TRH locally produced in MRNs are related to serotonergic regulation of brain processes during a chronic social conflict.

How cancer cells remodel lipid metabolism: strategies targeting transcription factors

Reprogramming of lipid metabolism has received increasing recognition as a hallmark of cancer cells because lipid dysregulation and the alteration of related enzyme profiles are closely correlated with oncogenic signals and malignant phenotypes, such as metastasis and therapeutic resistance. In this review, we describe recent findings that support the importance of lipids, as well as the transcription factors involved in cancer lipid metabolism. With recent advances in transcription factor analysis, including computer-modeling techniques, transcription factors are emerging as central players in cancer biology. Considering the limited number and the crucial role of transcription factors associated with lipid rewiring in cancers, transcription factor targeting is a promising potential strategy for cancer therapy.

ImmReg: the regulon atlas of immune-related pathways across cancer types

Immune system gene regulation perturbation has been found to be a major cause of the development of various types of cancer. Numbers of mechanisms contribute to gene expression regulation, thus, systematically identification of potential regulons of immune-related pathways is critical to cancer immunotherapy. Here, we comprehensively chart the landscape of transcription factors, microRNAs, RNA binding proteins and long noncoding RNAs regulation in 17 immune-related pathways across 33 cancers. The potential immunology regulons are likely to exhibit higher expressions in immune cells, show expression perturbations in cancer, and are significantly correlated with immune cell infiltrations. We also identify a panel of clinically relevant immunology regulons across cancers. Moreover, the regulon atlas of immune-related pathways helps prioritizing cancer-related genes (i.e. ETV7, miR-146a-5p, ZFP36 and HCP5). We further identified two molecular subtypes of glioma (cold and hot tumour phenotypes), which were characterized by differences in immune cell infiltrations, expression of checkpoints, and prognosis. Finally, we developed a user-friendly resource, ImmReg (http://bio-bigdata.hrbmu.edu.cn/ImmReg/), with multiple modules to visualize, browse, and download immunology regulation. Our study provides a comprehensive landscape of immunology regulons, which will shed light on future development of RNA-based cancer immunotherapies.

Proteome-wide profiling of transcriptional machinery on accessible chromatin with biotinylated transposons

To directly and quantitatively identify the transcriptional protein complexes assembled on accessible chromatin, we develop an assay for transposase-accessible chromatin using mass spectrum (ATAC-MS) based on direct transposition of biotinylated adaptors into open chromatin. Coupling with activated gene sequence information by ATAC-seq, ATAC-MS can profile the accessible chromatin-protein machinery. ATAC-MS, combined with fractionation strategies (fATAC-MS), can provide a high-resolution chromatin-transcriptional machinery atlas. ATAC-MS with a novel Tn5-dCas9 fusion protein [dCas9-targeted ATAC-MS (ctATAC-MS)] further facilitates systematic pinpointing of the transcriptional machinery at specific open chromatin regions. We used ATAC-MS and ATAC-seq to investigate transcriptional regulation during C2C12 cell differentiation and demonstrated the role of RFX1 in regulating the proliferation and differentiation of C2C12 cells. Our strategy provides a universal toolbox including ATAC-MS, fATAC-MS, and ctATAC-MS, which enables us to portray the transcriptional regulation machinery atlas in genome scale and investigate the protein-DNA complex at a specific genomic locus.

Proteome-Wide Profiling of Readers for DNA Modification

DNA modifications, represented by 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC), play important roles in epigenetic regulation of biological processes. The specific recognition of DNA modifications by the transcriptional protein machinery is thought to be a potential mechanism for epigenetic-driven gene regulation, and many modified DNA-specific binding proteins have been uncovered. However, the panoramic view of the roles of DNA modification readers at the proteome level remains largely unclear. Here, a recently developed concatenated tandem array of consensus transcription factor (TF) response elements (catTFREs) approach is employed to profile the binding activity of TFs at DNA modifications. Modified DNA-binding activity is quantified for 1039 TFs, representing 70% of the TFs in the human genome. Additionally, the modified DNA-binding activity of 600 TFs is monitored during the mouse brain development from the embryo to the adult stages. Readers of these DNA modifications are predicted, and the hierarchical networks between the transcriptional protein machinery and modified DNA are described. It is further demonstrated that ZNF24 and ZSCAN21 are potential readers of 5fC-modified DNA. This study provides a landscape of TF-DNA modification interactions that can be used to elucidate the epigenetic-related transcriptional regulation mechanisms under physiological conditions.

Functional and regulatory aspects of oxidative stress response in X monosomy

The partial/complete loss of one X chromosome in a human female leads to Turner syndrome (TS). TS individuals display a range of phenotypes including short stature, osteoporosis, ovarian malfunction, diabetes, and thyroid dysfunction. Epigenetic factors and regulatory networks are distinctly different in X monosomy (45, X). In a lifetime, an individual is exposed to a variety of stress conditions. To study whether X monosomy cells display a differential response upon exposure to mild stress as compared to normal 46, XX cells and whether this may contribute to various co-morbidities in aneuploid individuals, we have carried out a transcriptomic analysis of human fibroblasts 45, X and 46, XX after exposure to mild oxidative stress. Under these conditions, over 350 transcripts were seen to be differentially expressed in 45, X and 46, XX cells. Pathways associated with oxidative stress were differentially expressed highlighting the differential regulation of genes and associated phenotypes. It could be seen that X monosomy cells are more susceptible to oxidative stress as compared to normal cells and have altered molecular pathways both in normal conditions and also upon exposure to mild oxidative stress. To explore this aspect in detail, we have mapped the expressions of transcription factors (TFs) in 45, X and 46, XX cells. The network of transcription activating factors is differentially regulated in 45, X and 46, XX cells under stress exposure. It is tempting to speculate that the altered ability of 45, X (Turner) cells to respond to stress may play a significant role in the physiological function and altered phenotypes in Turner syndrome.

ReFeaFi: Genome-wide prediction of regulatory elements driving transcription initiation

Regulatory elements control gene expression through transcription initiation (promoters) and by enhancing transcription at distant regions (enhancers). Accurate identification of regulatory elements is fundamental for annotating genomes and understanding gene expression patterns. While there are many attempts to develop computational promoter and enhancer identification methods, reliable tools to analyze long genomic sequences are still lacking. Prediction methods often perform poorly on the genome-wide scale because the number of negatives is much higher than that in the training sets. To address this issue, we propose a dynamic negative set updating scheme with a two-model approach, using one model for scanning the genome and the other one for testing candidate positions. The developed method achieves good genome-level performance and maintains robust performance when applied to other vertebrate species, without re-training. Moreover, the unannotated predicted regulatory regions made on the human genome are enriched for disease-associated variants, suggesting them to be potentially true regulatory elements rather than false positives. We validated high scoring "false positive" predictions using reporter assay and all tested candidates were successfully validated, demonstrating the ability of our method to discover novel human regulatory regions.

Emerging Complexity in CD4+T Lineage Programming and Its Implications in Colorectal Cancer

The intestinal immune system has the difficult task of protecting a large environmentally exposed single layer of epithelium from pathogens without allowing inappropriate inflammatory responses. Unmitigated inflammation drives multiple pathologies, including the development of colorectal cancer. CD4⁺T cells mediate both the suppression and promotion of intestinal inflammation. They comprise an array of phenotypically and functionally distinct subsets tailored to a specific inflammatory context. This diversity of form and function is relevant to a broad array of pathologic and physiologic processes. The heterogeneity underlying both effector and regulatory T helper cell responses to colorectal cancer, and its impact on disease progression, is reviewed herein. Importantly, T cell responses are dynamic; they exhibit both quantitative and qualitative changes as the inflammatory context shifts. Recent evidence outlines the role of CD4⁺T cells in colorectal cancer responses and suggests possible mechanisms driving qualitative alterations in anti-cancer immune responses. The heterogeneity of T cells in colorectal cancer, as well as the manner and mechanism by which they change, offer an abundance of opportunities for more specific, and likely effective, interventional strategies.

Genetically engineered mouse models of esophageal cancer

BACKGROUND

Esophageal cancer is the most common cause of cancer-related death worldwide with a diverse geographical distribution, poor prognosis, and diagnosis in advanced stages of the disease. Identification of the mechanisms involved in esophageal cancer development is evaluative to improve outcomes for patients. Genetically engineered mouse models (GEMMs) of cancer provide the physiologic, molecular, and histologic features of the human tumors to determine the pathogenesis and treatments for cancer, hence exhibiting a source of tremendous potential for oncology research. The advancement of cancer modeling in mice has improved to the extent that researchers can observe and manipulate the disease process in a specific manner. Despite the significant differences between mice and humans, mice can be great models for human oncology researches due to similarities between them at the molecular and physiological levels. Due to most of the existing esophageal cancer GEMMs do not propose an ideal system for pathogenesis of the disease, genetic risks, and microenvironment exposure, so identification of challenges in GEM modeling and well-developed technologies are required to obtain the most value for patients. In this review, we describe the biology of human and mouse, followed by the exciting esophageal cancer mouse models with a discussion of applicability and challenges of these models for generating new GEMMs in future studies.

NGN2 induces diverse neuron types from human pluripotency

Human neurons engineered from induced pluripotent stem cells (iPSCs) through neurogenin 2 (NGN2) overexpression are widely used to study neuronal differentiation mechanisms and to model neurological diseases. However, the differentiation paths and heterogeneity of emerged neurons have not been fully explored. Here, we used single-cell transcriptomics to dissect the cell states that emerge during NGN2 overexpression across a time course from pluripotency to neuron functional maturation. We find a substantial molecular heterogeneity in the neuron types generated, with at least two populations that express genes associated with neurons of the peripheral nervous system. Neuron heterogeneity is observed across multiple iPSC clones and lines from different individuals. We find that neuron fate acquisition is sensitive to NGN2 expression level and the duration of NGN2-forced expression. Our data reveal that NGN2 dosage can regulate neuron fate acquisition, and that NGN2-iN heterogeneity can confound results that are sensitive to neuron type.

•

NGN2-iNs are molecularly heterogeneous

•

NGN2-iNs subtypes have signatures of central and peripheral nervous system

•

Neural fate acquisition is sensitive to the level and duration of NGN2 expression

Unveiling the Pathogenesis of Psychiatric Disorders Using Network Models

Psychiatric disorders are complex brain disorders with a high degree of genetic heterogeneity, affecting millions of people worldwide. Despite advances in psychiatric genetics, the underlying pathogenic mechanisms of psychiatric disorders are still largely elusive, which impedes the development of novel rational therapies. There has been accumulating evidence suggesting that the genetics of complex disorders can be viewed through an omnigenic lens, which involves contextualizing genes in highly interconnected networks. Thus, applying network-based multi-omics integration methods could cast new light on the pathophysiology of psychiatric disorders. In this review, we first provide an overview of the recent advances in psychiatric genetics and highlight gaps in translating molecular associations into mechanistic insights. We then present an overview of network methodologies and review previous applications of network methods in the study of psychiatric disorders. Lastly, we describe the potential of such methodologies within a multi-tissue, multi-omics approach, and summarize the future directions in adopting diverse network approaches.

Time-sequenced transcriptomes of developing distal mouse limb buds: A comparative tissue layer analysis

BACKGROUND

The development of the amniote limb has been an important model system to study patterning mechanisms and morphogenesis. For proper growth and patterning, it requires the interaction between the distal sub-apical mesenchyme and the apical ectodermal ridge (AER) that involve the separate implementation of coordinated and tissue-specific genetic programs.

RESULTS

Here, we produce and analyze the transcriptomes of both distal limb mesenchymal progenitors and the overlying ectodermal cells, following time-coursed dissections that cover from limb bud initiation to fully patterned limbs. The comparison of transcriptomes within each layer as well as between layers over time, allowed the identification of specific transcriptional signatures for each of the developmental stages. Special attention was given to the identification of genes whose transcription dynamics suggest a previously unnoticed role in the context of limb development and also to signaling pathways enriched between layers.

CONCLUSION

We interpret the transcriptomic data in light of the known development pattern and we conclude that a major transcriptional transition occurs in distal limb buds between E9.5 and E10.5, coincident with the switch from an early phase continuation of the signature of trunk progenitors, related to the initial proximo distal specification, to a late intrinsic phase of development.

Overexpression of Replication-Dependent Histone Signifies a Subset of Dedifferentiated Liposarcoma with Increased Aggressiveness

Liposarcoma (LPS) is an adult soft tissue malignancy that arises from fat tissue, where well-differentiated (WD) and dedifferentiated (DD) forms are the most common. DDLPS represents the progression of WDLPS into a more aggressive high-grade and metastatic form. Although a few DNA copy-number amplifications are known to be specifically found in WD- or DDLPS, systematic genetic differences that signify subtype determination between WDLPS and DDLPS remain unclear. Here, we profiled the genome and transcriptome of 38 LPS tumors to uncover the genetic signatures of subtype differences. Replication-dependent histone (RD-HIST) mRNAs were highly elevated and their regulation was disrupted in a subset of DDLPS, increasing cellular histone molecule levels, as measured using RNA-seq (the averaged fold change of 53 RD-HIST genes between the DD and WD samples was 10.9) and immunohistochemistry. The change was not observed in normal tissues. Integrated whole-exome sequencing, RNA-seq, and methylation analyses revealed that the overexpressed HMGA2 (the fold change between DD and WD samples was 7.3) was responsible for the increased RD-HIST level, leading to aberrant cell proliferation. Therefore, HMGA2-mediated elevation of RD-HISTs were crucial events in determining the aggressiveness of DDLPS, which may serve as a biomarker for prognosis prediction for liposarcoma patients.

Machine learning for profile prediction in genomics

A recent deluge of publicly available multi-omics data has fueled the development of machine learning methods aimed at investigating important questions in genomics. Although the motivations for these methods vary, a task that is commonly adopted is that of profile prediction, where predictions are made for one or more forms of biochemical activity along the genome, for example, histone modification, chromatin accessibility, or protein binding. In this review, we give an overview of the research works performing profile prediction, define two broad categories of profile prediction tasks, and discuss the types of scientific questions that can be answered in each.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Widespread Exaptation of L1 Transposons for Transcription Factor Binding in Breast Cancer

L1 transposons occupy 17% of the human genome and are widely exapted for the regulation of human genes, particularly in breast cancer, where we have previously shown abundant cancer-specific transcription factor (TF) binding sites within the L1PA2 subfamily. In the current study, we performed a comprehensive analysis of TF binding activities in primate-specific L1 subfamilies and identified pervasive exaptation events amongst these evolutionarily related L1 transposons. By motif scanning, we predicted diverse and abundant TF binding potentials within the L1 transposons. We confirmed substantial TF binding activities in the L1 subfamilies using TF binding sites consolidated from an extensive collection of publicly available ChIP-seq datasets. Young L1 subfamilies (L1HS, L1PA2 and L1PA3) contributed abundant TF binding sites in MCF7 cells, primarily via their 5' UTR. This is expected as the L1 5' UTR hosts cis-regulatory elements that are crucial for L1 replication and mobilisation. Interestingly, the ancient L1 subfamilies, where 5' truncation was common, displayed comparable TF binding capacity through their 3' ends, suggesting an alternative exaptation mechanism in L1 transposons that was previously unnoticed. Overall, primate-specific L1 transposons were extensively exapted for TF binding in MCF7 breast cancer cells and are likely prominent genetic players modulating breast cancer transcriptional regulation.

Control of mouse limb initiation and antero-posterior patterning by Meis transcription factors

Meis1 and Meis2 are homeodomain transcription factors that regulate organogenesis through cooperation with Hox proteins. Elimination of Meis genes after limb induction has shown their role in limb proximo-distal patterning; however, limb development in the complete absence of Meis function has not been studied. Here, we report that Meis1/2 inactivation in the lateral plate mesoderm of mouse embryos leads to limb agenesis. Meis and Tbx factors converge in this function, extensively co-binding with Tbx to genomic sites and co-regulating enhancers of Fgf10, a critical factor in limb initiation. Limbs with three deleted Meis alleles show proximal-specific skeletal hypoplasia and agenesis of posterior skeletal elements. This failure in posterior specification results from an early role of Meis factors in establishing the limb antero-posterior prepattern required for Shh activation. Our results demonstrate roles for Meis transcription factors in early limb development and identify their involvement in previously undescribed interaction networks that regulate organogenesis.

Construction of Condition-Specific Gene Regulatory Network Using Kernel Canonical Correlation Analysis

Gene expression profile or transcriptome can represent cellular states, thus understanding gene regulation mechanisms can help understand how cells respond to external stress. Interaction between transcription factor (TF) and target gene (TG) is one of the representative regulatory mechanisms in cells. In this paper, we present a novel computational method to construct condition-specific transcriptional networks from transcriptome data. Regulatory interaction between TFs and TGs is very complex, specifically multiple-to-multiple relations. Experimental data from TF Chromatin Immunoprecipitation sequencing is useful but produces one-to-multiple relations between TF and TGs. On the other hand, co-expression networks of genes can be useful for constructing condition transcriptional networks, but there are many false positive relations in co-expression networks. In this paper, we propose a novel method to construct a condition-specific and combinatorial transcriptional network, applying kernel canonical correlation analysis (kernel CCA) to identify multiple-to-multiple TF-TG relations in certain biological condition. Kernel CCA is a well-established statistical method for computing the correlation of a group of features vs. another group of features. We, therefore, employed kernel CCA to embed TFs and TGs into a new space where the correlation of TFs and TGs are reflected. To demonstrate the usefulness of our network construction method, we used the blood transcriptome data for the investigation on the response to high fat diet in a human and an arabidopsis data set for the investigation on the response to cold/heat stress. Our method detected not only important regulatory interactions reported in previous studies but also novel TF-TG relations where a module of TF is regulating a module of TGs upon specific stress.

Positive Selection in Gene Regulatory Factors Suggests Adaptive Pleiotropic Changes During Human Evolution

Gene regulatory factors (GRFs), such as transcription factors, co-factors and histone-modifying enzymes, play many important roles in modifying gene expression in biological processes. They have also been proposed to underlie speciation and adaptation. To investigate potential contributions of GRFs to primate evolution, we analyzed GRF genes in 27 publicly available primate genomes. Genes coding for zinc finger (ZNF) proteins, especially ZNFs with a Krüppel-associated box (KRAB) domain were the most abundant TFs in all genomes. Gene numbers per TF family differed between all species. To detect signs of positive selection in GRF genes we investigated more than 3,000 human GRFs with their more than 70,000 orthologs in 26 non-human primates. We implemented two independent tests for positive selection, the branch-site-model of the PAML suite and aBSREL of the HyPhy suite, focusing on the human and great ape branch. Our workflow included rigorous procedures to reduce the number of false positives: excluding distantly similar orthologs, manual corrections of alignments, and considering only genes and sites detected by both tests for positive selection. Furthermore, we verified the candidate sites for selection by investigating their variation within human and non-human great ape population data. In order to approximately assign a date to positively selected sites in the human lineage, we analyzed archaic human genomes. Our work revealed with high confidence five GRFs that have been positively selected on the human lineage and one GRF that has been positively selected on the great ape lineage. These GRFs are scattered on different chromosomes and have been previously linked to diverse functions. For some of them a role in speciation and/or adaptation can be proposed based on the expression pattern or association with human diseases, but it seems that they all contributed independently to human evolution. Four of the positively selected GRFs are KRAB-ZNF proteins, that induce changes in target genes co-expression and/or through arms race with transposable elements. Since each positively selected GRF contains several sites with evidence for positive selection, we suggest that these GRFs participated pleiotropically to phenotypic adaptations in humans.

Single-Nucleotide Polymorphisms Associated with the Senescence-Accelerated Phenotype of OXYS Rats: A Focus on Alzheimer's Disease-Like and Age-Related-Macular-Degeneration-Like Pathologies

Alzheimer's disease (AD) and age-related macular degeneration (AMD) are two complex incurable neurodegenerative disorders the common pathogenesis of which is actively discussed. There are overlapping risk factors and molecular mechanisms of the two diseases; at the same time, there are arguments in favor of the notion that susceptibility to each of these diseases is associated with a distinct genetic background. Here we identified single-nucleotide polymorphisms (SNPs) that are specific for senescence-accelerated OXYS rats, which simulate key characteristics of both sporadic AD and AMD. Transcriptomes of the hippocampus, prefrontal cortex, and retina (data of RNA-Seq) were analyzed. We detected SNPs in genes Rims2, AABR07072639.2, Lemd2, and AABR07045405.1, which thus can express significantly truncated proteins lacking functionally important domains. Additionally, 33 mutations in genes-which are related to various metabolic and signaling pathways-cause nonsynonymous amino acid substitutions presumably leading to disturbances in protein structure or functions. Some of the genes carrying these SNPs are associated with aging, neurodegenerative, and mental diseases. Thus, we revealed the SNPs can lead to abnormalities in protein structure or functions and affect the development of the senescence-accelerated phenotype of OXYS rats. Our data are consistent with the latest results of genome-wide association studies that highlight the importance of multiple pathways for the pathogenesis of AD and AMD. Identified SNPs can serve as promising research objects for further studies on the molecular mechanisms underlying this particular rat model as well as for the prediction of potential biomarkers of AD and AMD.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Physiologic Responses to Dietary Sulfur Amino Acid Restriction in Mice Are Influenced by Atf4 Status and Biological Sex

BACKGROUND

Dietary sulfur amino acid restriction (SAAR) improves body composition and metabolic health across several model organisms in part through induction of the integrated stress response (ISR).

OBJECTIVE

We investigate the hypothesis that activating transcription factor 4 (ATF4) acts as a converging point in the ISR during SAAR.

METHODS

Using liver-specific or global gene ablation strategies, in both female and male mice, we address the role of ATF4 during dietary SAAR.

RESULTS

We show that ATF4 is dispensable in the chronic induction of the hepatokine fibroblast growth factor 21 while being essential for the sustained production of endogenous hydrogen sulfide. We also affirm that biological sex, independent of ATF4 status, is a determinant of the response to dietary SAAR.

CONCLUSIONS

Our results suggest that auxiliary components of the ISR, which are independent of ATF4, are critical for SAAR-mediated improvements in metabolic health in mice.

Genome-wide screening of upstream transcription factors using an expression library

The identification of upstream transcription factors regulating the expression of a gene is generally not an easy process.  To facilitate this task, we constructed an expression cDNA library named Transcription Factor Expression Library (TFEL), which is composed of nearly all the transcription factors in the mouse genome. Genome-wide screening using this library (TFEL scan method) enables us to easily identify transcription factors controlling any given promoter or enhancer of interest in a chromosomal context-dependent manner. Thus, TFEL scan method is a powerful approach to explore transcriptional regulatory networks.

Transcriptional circuitry atlas of genetic diverse unstimulated murine and human macrophages define disparity in population-wide innate immunity

Macrophages are ubiquitous custodians of tissues, which play decisive role in maintaining cellular homeostasis through regulatory immune responses. Within tissues, macrophage exhibit extremely heterogeneous population with varying functions orchestrated through regulatory response, which can be further exacerbated in diverse genetic backgrounds. Gene regulatory networks (GRNs) offer comprehensive understanding of cellular regulatory behavior by unfolding the transcription factors (TFs) and regulated target genes. RNA-Seq coupled with ATAC-Seq has revolutionized the regulome landscape influenced by gene expression modeling. Here, we employ an integrative multi-omics systems biology-based analysis and generated GRNs derived from the unstimulated bone marrow-derived macrophages of five inbred genetically defined murine strains, which are reported to be linked with most of the population-wide human genetic variants. Our probabilistic modeling of a basal hemostasis pan regulatory repertoire in diverse macrophages discovered 96 TFs targeting 6279 genes representing 468,291 interactions across five inbred murine strains. Subsequently, we identify core and distinctive GRN sub-networks in unstimulated macrophages to describe the system-wide conservation and dissimilarities, respectively across five murine strains. Our study concludes that discrepancies in unstimulated macrophage-specific regulatory networks not only drives the basal functional plasticity within genetic backgrounds, additionally aid in understanding the complexity of racial disparity among the human population during stress.

Investigating the molecular guts of endoderm formation using zebrafish

The vertebrate endoderm makes major contributions to the respiratory and gastrointestinal tracts and all associated organs. Zebrafish and humans share a high degree of genetic homology and strikingly similar endodermal organ systems. Combined with a multitude of experimental advantages, zebrafish are an attractive model organism to study endoderm development and disease. Recent functional genomics studies have shed considerable light on the gene regulatory programs governing early zebrafish endoderm development, while advances in biological and technological approaches stand to further revolutionize our ability to investigate endoderm formation, function and disease. Here, we discuss the present understanding of endoderm specification in zebrafish compared to other vertebrates, how current and emerging methods will allow refined and enhanced analysis of endoderm formation, and how integration with human data will allow modeling of the link between non-coding sequence variants and human disease.

Role of the bHLH transcription factor TCF21 in development and tumorigenesis

Transcription factors control, coordinate, and separate the functions of distinct network modules spatially and temporally. In this review, we focus on the transcription factor 21 (TCF21) network, a highly conserved basic-helix-loop-helix (bHLH) protein that functions to integrate signals and modulate gene expression. We summarize the molecular and biological properties of TCF21 control with an emphasis on molecular and functional TCF21 interactions. We suggest that these interactions serve to modulate the development of different organs at the transcriptional level to maintain growth homeostasis and to influence cell fate. Importantly, TCF21 expression is epigenetically inactivated in different types of human cancers. The epigenetic modification or activation and/or loss of TCF21 expression results in an imbalance in TCF21 signaling, which may lead to tumor initiation and, most likely, to progression and tumor metastasis. This review focuses on research on the roles of TCF21 in development and tumorigenesis systematically considering the physiological and pathological function of TCF21. In addition, we focus on the main molecular bases of its different roles whose importance should be clarified in future research. For this review, PubMed databases and keywords such as TCF21, POD-1, capsulin, tumors, carcinomas, tumorigenesis, development, and mechanism of action were utilized. Articles were selected within a historical context as were a number of citations from journals with relevant impact.

A computer-guided design tool to increase the efficiency of cellular conversions

Human cell conversion technology has become an important tool for devising new cell transplantation therapies, generating disease models and testing gene therapies. However, while transcription factor over-expression-based methods have shown great promise in generating cell types in vitro, they often endure low conversion efficiency. In this context, great effort has been devoted to increasing the efficiency of current protocols and the development of computational approaches can be of great help in this endeavor. Here we introduce a computer-guided design tool that combines a computational framework for prioritizing more efficient combinations of instructive factors (IFs) of cellular conversions, called IRENE, with a transposon-based genomic integration system for efficient delivery. Particularly, IRENE relies on a stochastic gene regulatory network model that systematically prioritizes more efficient IFs by maximizing the agreement of the transcriptional and epigenetic landscapes between the converted and target cells. Our predictions substantially increased the efficiency of two established iPSC-differentiation protocols (natural killer cells and melanocytes) and established the first protocol for iPSC-derived mammary epithelial cells with high efficiency.

Cellular transcriptomics reveals evolutionary identities of songbird vocal circuits

Birds display advanced behaviors, including vocal learning and problem-solving, yet lack a layered neocortex, a structure associated with complex behavior in mammals. To determine whether these behavioral similarities result from shared or distinct neural circuits, we used single-cell RNA sequencing to characterize the neuronal repertoire of the songbird song motor pathway. Glutamatergic vocal neurons had considerable transcriptional similarity to neocortical projection neurons; however, they displayed regulatory gene expression patterns more closely related to neurons in the ventral pallium. Moreover, while γ-aminobutyric acid-releasing neurons in this pathway appeared homologous to those in mammals and other amniotes, the most abundant avian class is largely absent in the neocortex. These data suggest that songbird vocal circuits and the mammalian neocortex have distinct developmental origins yet contain transcriptionally similar neurons.

Dynamic sex chromosome expression in Drosophila male germ cells

Given their copy number differences and unique modes of inheritance, the evolved gene content and expression of sex chromosomes is unusual. In many organisms the X and Y chromosomes are inactivated in spermatocytes, possibly as a defense mechanism against insertions into unpaired chromatin. In addition to current sex chromosomes, Drosophila has a small gene-poor X-chromosome relic (4th) that re-acquired autosomal status. Here we use single cell RNA-Seq on fly larvae to demonstrate that the single X and pair of 4th chromosomes are specifically inactivated in primary spermatocytes, based on measuring all genes or a set of broadly expressed genes in testis we identified. In contrast, genes on the single Y chromosome become maximally active in primary spermatocytes. Reduced X transcript levels are due to failed activation of RNA-Polymerase-II by phosphorylation of Serine 2 and 5.

A comprehensive survey for human transcription factors on expression, regulation, interaction, phenotype and cancer survival

Transcription factors (TFs) act as key regulators in biological processes through controlling gene expression. Here, we conducted a systematic study for all human TFs on the expression, regulation, interaction, mutation, phenotype and cancer survival. We revealed that the average expression levels of TFs in normal tissues were lower than 50% expression of non-TFs, whereas TF expression was increased in cancers. TFs that are specifically expressed in an individual tissue or cancer may be potential marker genes. For instance, TGIF2LX/Y were preferentially expressed in testis and NEUROG1, PRDM14, SRY, ZNF705A and ZNF716 were specifically highly expressed in germ cell tumors. We found different distributions of target genes and TF co-regulations in different TF families. Some small TF families have huge protein interaction pairs, suggesting their central roles in transcriptional regulation. The bZIP family is a small family involving many signaling pathways. Survival analysis indicated that most TFs significantly affect survival of one or more cancers. Some survival-related TFs were also specifically highly expressed in the corresponding cancer types, which may be potential targets for cancer therapy. Finally, we identified 43 TFs whose mutations were closely correlated to survival, suggesting their cancer-driven roles. The systematic analysis of TFs provides useful clues for further investigation of TF regulatory mechanisms and the role of TFs in diseases.

Structure-function relationships explain CTCF zinc finger mutation phenotypes in cancer

CCCTC-binding factor (CTCF) plays fundamental roles in transcriptional regulation and chromatin architecture maintenance. CTCF is also a tumour suppressor frequently mutated in cancer, however, the structural and functional impact of mutations have not been examined. We performed molecular and structural characterisation of five cancer-specific CTCF missense zinc finger (ZF) mutations occurring within key intra- and inter-ZF residues. Functional characterisation of CTCF ZF mutations revealed a complete (L309P, R339W, R377H) or intermediate (R339Q) abrogation as well as an enhancement (G420D) of the anti-proliferative effects of CTCF. DNA binding at select sites was disrupted and transcriptional regulatory activities abrogated. Molecular docking and molecular dynamics confirmed that mutations in residues specifically contacting DNA bases or backbone exhibited loss of DNA binding. However, R339Q and G420D were stabilised by the formation of new primary DNA bonds, contributing to gain-of-function. Our data confirm that a spectrum of loss-, change- and gain-of-function impacts on CTCF zinc fingers are observed in cell growth regulation and gene regulatory activities. Hence, diverse cellular phenotypes of mutant CTCF are clearly explained by examining structure-function relationships.

Transcriptomic Changes of Murine Visceral Fat Exposed to Intermittent Hypoxia at Single Cell Resolution

Intermittent hypoxia (IH) is a hallmark of obstructive sleep apnea (OSA) and induces metabolic dysfunction manifesting as inflammation, increased lipolysis and insulin resistance in visceral white adipose tissues (vWAT). However, the cell types and their corresponding transcriptional pathways underlying these functional perturbations are unknown. Here, we applied single nucleus RNA sequencing (snRNA-seq) coupled with aggregate RNA-seq methods to evaluate the cellular heterogeneity in vWAT following IH exposures mimicking OSA. C57BL/6 male mice were exposed to IH and room air (RA) for 6 weeks, and nuclei from vWAT were isolated and processed for snRNA-seq followed by differential expressed gene (DEGs) analyses by cell type, along with gene ontology and canonical pathways enrichment tests of significance. IH induced significant transcriptional changes compared to RA across 14 different cell types identified in vWAT. We identified cell-specific signature markers, transcriptional networks, metabolic signaling pathways, and cellular subpopulation enrichment in vWAT. Globally, we also identify 298 common regulated genes across multiple cellular types that are associated with metabolic pathways. Deconvolution of cell types in vWAT using global RNA-seq revealed that distinct adipocytes appear to be differentially implicated in key aspects of metabolic dysfunction. Thus, the heterogeneity of vWAT and its response to IH at the cellular level provides important insights into the metabolic morbidity of OSA and may possibly translate into therapeutic targets.

ZBTB20 regulates WNT/CTNNB1 signalling pathway by suppressing PPARG during hepatocellular carcinoma tumourigenesis

Background & Aims

Zinc finger and BTB domain containing 20 (ZBTB20) has been implicated as a potential oncogene in liver cancer. However, knockout studies have shown it to be a transcriptional repressor of the alpha-foetoprotein (Afp) gene in adult liver, and reduced levels of ZBTB20 allow for upregulation of AFP with increased tumour severity in certain cases of hepatocellular carcinoma (HCC). As there are many discrepancies in the literature regarding its role in liver tumourigenesis, the aim of this study was to elucidate the role of ZBTB20 in HCC tumourigenesis.

Methods

A reverse genetic study using the Sleeping Beauty (SB) transposon system in mice was performed to elucidate the role of ZBTB20 in HCC tumourigenesis. In vitro ZBTB20 gain- and loss-of-function experiments were used to assess the relationship amongst ZBTB20, peroxisome proliferator activated receptor gamma (PPARG) and catenin beta 1 (CTNNB1).

Results

Transgenic overexpression of ZBTB20 in hepatocytes and in the context of transformation related protein (T r p53) inactivation induced hepatic hypertrophy, activation of WNT/CTNNB1 signalling, and development of liver tumours. In vitro overexpression and knockout experiments using CRISPR/Cas9 demonstrated the important role for ZBTB20 in downregulating PPARG, resulting in activation of the WNT/CTNNB1 signalling pathway and its downstream effectors in HCC tumourigenesis.

Conclusions

These findings demonstrate a novel interaction between ZBTB20 and PPARG, which leads to activation of the WNT/CTNNB1 signalling pathway in HCC tumourigenesis.

Lay summary

ZBTB20 has been implicated as a potential oncogene in liver cancer. Herein, we uncover its important role in liver cancer development. We show that it interacts with PPARG to upregulate the WNT/CTNNB1 signalling pathway, leading to tumourigenesis.

Molecular Co-occupancy Identifies Transcription Factor Binding Cooperativity In Vivo

Gene activation requires the cooperative activity of multiple transcription factors at cis-regulatory elements (CREs). Yet, most transcription factors have short residence time, questioning the requirement of their physical co-occupancy on DNA to achieve cooperativity. Here, we present a DNA footprinting method that detects individual molecular interactions of transcription factors and nucleosomes with DNA in vivo. We apply this strategy to quantify the simultaneous binding of multiple transcription factors on single DNA molecules at mouse CREs. Analysis of the binary occupancy patterns at thousands of motif combinations reveals that high DNA co-occupancy occurs for most types of transcription factors, in the absence of direct physical interaction, at sites of competition with nucleosomes. Perturbation of pairwise interactions demonstrates the function of molecular co-occupancy in binding cooperativity. Our results reveal the interactions regulating CREs at molecular resolution and identify DNA co-occupancy as a widespread cooperativity mechanism used by transcription factors to remodel chromatin.

Mapping the epigenomic and transcriptomic interplay during memory formation and recall in the hippocampal engram ensemble

The epigenome and three-dimensional (3D) genomic architecture are emerging as key factors in the dynamic regulation of different transcriptional programs required for neuronal functions. In this study, we used an activity-dependent tagging system in mice to determine the epigenetic state, 3D genome architecture and transcriptional landscape of engram cells over the lifespan of memory formation and recall. Our findings reveal that memory encoding leads to an epigenetic priming event, marked by increased accessibility of enhancers without the corresponding transcriptional changes. Memory consolidation subsequently results in spatial reorganization of large chromatin segments and promoter-enhancer interactions. Finally, with reactivation, engram neurons use a subset of de novo long-range interactions, where primed enhancers are brought in contact with their respective promoters to upregulate genes involved in local protein translation in synaptic compartments. Collectively, our work elucidates the comprehensive transcriptional and epigenomic landscape across the lifespan of memory formation and recall in the hippocampal engram ensemble.

The Multifaceted Output of c-Jun Biological Activity: Focus at the Junction of CD8 T Cell Activation and Exhaustion

c-Jun is a major component of the dimeric transcription factor activator protein-1 (AP-1), a paradigm for transcriptional response to extracellular signaling, whose components are basic-Leucine Zipper (bZIP) transcription factors of the Jun, Fos, activating transcription factor (ATF), ATF-like (BATF) and Jun dimerization protein 2 (JDP2) gene families. Extracellular signals regulate c-Jun/AP-1 activity at multiple levels, including transcriptional and posttranscriptional regulation of c-Jun expression and transactivity, in turn, establishing the magnitude and the duration of c-Jun/AP-1 activation. Another important level of c-Jun/AP-1 regulation is due to the capability of Jun family members to bind DNA as a heterodimer with every other member of the AP-1 family, and to interact with other classes of transcription factors, thereby acquiring the potential to integrate diverse extrinsic and intrinsic signals into combinatorial regulation of gene expression. Here, we review how these features of c-Jun/AP-1 regulation underlie the multifaceted output of c-Jun biological activity, eliciting quite distinct cellular responses, such as neoplastic transformation, differentiation and apoptosis, in different cell types. In particular, we focus on the current understanding of the role of c-Jun/AP-1 in the response of CD8 T cells to acute infection and cancer. We highlight the transcriptional and epigenetic regulatory mechanisms through which c-Jun/AP-1 participates in the productive immune response of CD8 T cells, and how its downregulation may contribute to the dysfunctional state of tumor infiltrating CD8 T cells. Additionally, we discuss recent insights pointing at c-Jun as a suitable target for immunotherapy-based combination approaches to reinvigorate anti-tumor immune functions.

Shifts of Immune Cell Populations Differ in Response to Different Effectors of Beige Remodeling of Adipose Tissue

White adipose tissue (WAT) is a dynamic tissue, which responds to environmental stimuli and dietary cues by changing its morphology and metabolic capacity. The ability of WAT to undergo a beige remodeling has become an appealing strategy to combat obesity and its comorbidities. Here, by using single-cell RNA sequencing, we provide a comprehensive atlas of the cellular dynamics during beige remodeling. We reveal drastic changes both in the overall cellular composition and transcriptional states of individual cell subtypes between Adrb3- and cold-induced beiging. Moreover, we demonstrate that cold induces a myeloid to lymphoid shift of the immune compartment compared to Adrb3 activation. Further analysis showed that, Adrb3 stimulation leads to activation of the interferon/Stat1 pathways favoring infiltration of myeloid immune cells, while repression of this pathway by cold promotes lymphoid immune cell recruitment. These findings highlight that pharmacological mimetics may not provide the same beneficial effects as physiological stimuli.

Gene regulation could be attributed to TCF3 and other key transcription factors in the muscle of pubertal heifers

Puberty is a whole-body event, driven by the hypothalamic integration of peripheral signals such as leptin or IGF-1. In the process of puberty, reproductive development is simultaneous to growth, including muscle growth. To enhance our understanding of muscle function related to puberty, we performed transcriptome analyses of muscle samples from six pre- and six post-pubertal Brahman heifers (Bos indicus). Our aims were to perform differential expression analyses and co-expression analyses to derive a regulatory gene network associate with puberty. As a result, we identified 431 differentially expressed (DEx) transcripts (genes and non-coding RNAs) when comparing pre- to post-pubertal average gene expression. The DEx transcripts were compared with all expressed transcripts in our samples (over 14,000 transcripts) for functional enrichment analyses. The DEx transcripts were associated with "extracellular region," "inflammatory response" and "hormone activity" (adjusted p < .05). Inflammatory response for muscle regeneration is a necessary aspect of muscle growth, which is accelerated during puberty. The term "hormone activity" may signal genes that respond to progesterone signalling in the muscle, as the presence of this hormone is an important difference between pre- and post-pubertal heifers in our experimental design. The DEx transcript with the highest average expression difference was a mitochondrial gene, ENSBTAG00000043574 that might be another important link between energy metabolism and puberty. In the derived co-expression gene network, we identified six hub genes: CDC5L, MYC, TCF3, RUNX2, ATF2 and CREB1. In the same network, 48 key regulators of DEx transcripts were identified, using a regulatory impact factor metric. The hub gene TCF3 was also a key regulator. The majority of the key regulators (22 genes) are members of the zinc finger family, which has been implicated in bovine puberty in other tissues. In conclusion, we described how puberty may affect muscle gene expression in cattle.

Identification of Prognostic Markers in Cholangiocarcinoma Using Altered DNA Methylation and Gene Expression Profiles

Background

Cholangiocarcinoma (CCA) is a rare disease, but it is amongst the most lethal cancers with a median survival under 1 year. Variations in DNA methylation and gene expression have been extensively studied in other cancers for their role in pathogenesis and disease prognosis, but these studies are very limited in CCA. This study focusses on the identification of DNA methylation and gene expression prognostic biomarkers using multi-omics data of CCA tumors from The Cancer Genome Atlas (TCGA).

Method

We have conducted a genome-wide analysis of differential DNA methylation and gene/miRNA expression using data from 36 CCA tumor and 9 normal samples from TCGA. The impact of DNA methylation in promoters and long-range distal enhancers on the regulation and expression of CCA-associated genes was examined using linear regression. Next, we conducted network analyses on genes which are regulated by DNA methylation as well as by miRNA. Finally, we performed Kaplan–Meier and Cox proportional hazards regression analyses in order to identify the role of selected methylation sites and specific genes and miRNAs in patient survival. We also performed real-time quantitative PCR (qPCR) to confirm the change in gene expression in CCA patients’ tumor and adjacent normal samples.

Results

Altered DNA methylation was observed on 12,259 CpGs across all chromosomes, of which 78% were hypermethylated. We observed a strong negative relationship between promoter hypermethylation and corresponding gene expression in 92% of the CpGs. Differential expression analyses revealed altered expression patterns in 3,305 genes and 101 miRNAs. Finally, we identified 17 differentially methylated promoter CpGs, 72 differentially expressed genes, and two miRNAs that are likely associated with patient survival. Pathway analysis suggested that cell division, bile secretion, amino acid metabolism, PPAR signaling, hippo signaling were highly affected by gene expression and DNA methylation alterations. The qPCR analysis further confirmed that MDK, HNF1B, PACS1, and GLUD1 are differentially expressed in CCA.

Conclusion

Based on the survival analysis, we conclude that DEPDC1, FUT4, MDK, PACS1, PIWIL4 genes, miR-22, miR-551b microRNAs, and cg27362525 and cg26597242 CpGs can strongly support their use as prognostic markers of CCA.

Association of the Sp1 binding site and -1997 promoter variations in COL1A1 with osteoporosis risk: The application of meta-analysis and bioinformatics approaches offers a new perspective for future research

As a complex disease, osteoporosis is influenced by several genetic markers. Many studies have examined the link between the Sp1 binding site +1245 G > T (rs1800012) and -1997 G > T (rs1107946) variations in the COL1A1 gene with osteoporosis risk. However, the findings of these studies have been contradictory; therefore, we performed a meta-analysis to aggregate additional information and obtain increased statistical power to more efficiently estimate this correlation. A meta-analysis was conducted with studies published between 1991-2020 that were identified by a systematic electronic search of the Scopus and Clarivate Analytics databases. Studies with bone mineral density (BMD) data and complete genotypes of the single-nucleotide variations (SNVs) for the overall and postmenopausal female population were included in this meta-analysis and analyzed using the R metaphor package. A relationship between rs1800012 and significantly decreased BMD values at the lumbar spine and femoral neck was found in individuals carrying the "ss" versus the "SS" genotype in the overall population according to a random effects model (p < 0.0001). Similar results were also found in the postmenopausal female population (p = 0.003 and 0.0002, respectively). Such findings might be an indication of increased osteoporosis risk in both studied groups in individuals with the "ss" genotype. Although no association was identified between the -1997 G > T and low BMD in the overall population, those individuals with the "GT" genotype showed a higher level of BMD than those with "GG" in the subgroup analysis (p = 0.007). To determine which transcription factor (TF) might bind to the -1997 G > T in COL1A1, 45 TFs were identified based on bioinformatics predictions. According to the GSE35958 microarray dataset, 16 of 45 TFs showed differential expression profiles in osteoporotic human mesenchymal stem cells relative to normal samples from elderly donors. By identifying candidate TFs for the -1997 G > T site, our study offers a new perspective for future research.

Gene Expression Changes in the Ventral Tegmental Area of Male Mice with Alternative Social Behavior Experience in Chronic Agonistic Interactions

Daily agonistic interactions of mice are an effective experimental approach to elucidate the molecular mechanisms underlying the excitation of the brain neurons and the formation of alternative social behavior patterns. An RNA-Seq analysis was used to compare the ventral tegmental area (VTA) transcriptome profiles for three groups of male C57BL/6J mice: winners, a group of chronically winning mice, losers, a group of chronically defeated mice, and controls. The data obtained show that both winners and defeated mice experience stress, which however, has a more drastic effect on defeated animals causing more significant changes in the levels of gene transcription. Four genes (Nrgn, Ercc2, Otx2, and Six3) changed their VTA expression profiles in opposite directions in winners and defeated mice. It was first shown that Nrgn (neurogranin) expression was highly correlated with the expression of the genes involved in dopamine synthesis and transport (Th, Ddc, Slc6a3, and Drd2) in the VTA of defeated mice but not in winners. The obtained network of 31 coregulated genes, encoding proteins associated with nervous system development (including 24 genes associated with the generation of neurons), may be potentially useful for studying their role in the VTA dopaminergic neurons maturation under the influence of social stress.

Multiomics data integration unveils core transcriptional regulatory networks governing cell-type identity

A plethora of computational approaches have been proposed for reconstructing gene regulatory networks (GRNs) from gene expression data. However, gene regulatory processes are often too complex to predict from the transcriptome alone. Here, we present a computational method, Moni, that systematically integrates epigenetics, transcriptomics, and protein–protein interactions to reconstruct GRNs among core transcription factors and their co-factors governing cell identity. We applied Moni to 57 datasets of human cell types and lines and demonstrate that it can accurately infer GRNs, thereby outperforming state-of-the-art methods.

Time- and Dose-Resolved Proteome of PM2.5-Exposure-Induced Lung Injury and Repair in Rats

In recent years, airborne fine particulate matter (PM_2.5) is drawing more public attention due to its various physicochemical features and causing pathological harm, as proven by epidemiological and clinical studies. However, the mechanism of PM_2.5-exposure-induced lung injury has not been fully characterized. Here, we established a PM_2.5-induced rat injury model for both short-term and long-term exposures at different concentrations. We employed the Fast-seq technique to profile 6316 proteins and the catTFRE approach to profile 387 transcription factors (TFs) in the lung tissue. In short-term exposure, we elucidated gradually upregulated proteins enriched in response to oxidative stress, phagosome, and the extracellular matrix (ECM)-receptor interaction pathway. Long-term exposure mainly showed the immune response pathway to be consisting of increased lymphocytes and cytokines. Intriguingly, we found that immune-related proteins were recoverable during short-term exposure. During the process of PM_2.5 exposure, upregulated proteins presented dose-dependence in the lung, including stress response at low dose, minor immune response at middle dose, and severe inflammatory response at high dose. This data set provides a rich resource to facilitate the understanding of PM_2.5-induced lung damage and repair mechanism.