i-cisTarget: an integrative genomics method for the prediction of regulatory features and cis-regulatory modules

Metformin decelerates aging clock in male monkeys

In a rigorous 40-month study, we evaluated the geroprotective effects of metformin on adult male cynomolgus monkeys, addressing a gap in primate aging research. The study encompassed a comprehensive suite of physiological, imaging, histological, and molecular evaluations, substantiating metformin's influence on delaying age-related phenotypes at the organismal level. Specifically, we leveraged pan-tissue transcriptomics, DNA methylomics, plasma proteomics, and metabolomics to develop innovative monkey aging clocks and applied these to gauge metformin's effects on aging. The results highlighted a significant slowing of aging indicators, notably a roughly 6-year regression in brain aging. Metformin exerts a substantial neuroprotective effect, preserving brain structure and enhancing cognitive ability. The geroprotective effects on primate neurons were partially mediated by the activation of Nrf2, a transcription factor with anti-oxidative capabilities. Our research pioneers the systemic reduction of multi-dimensional biological age in primates through metformin, paving the way for advancing pharmaceutical strategies against human aging.

Cross-species single-cell landscapes identify the pathogenic gene characteristics of inherited retinal diseases

Introduction

Inherited retinal diseases (IRDs) affect ∼4.5 million people worldwide. Elusive pathogenic variants in over 280 genes are associated with one or more clinical forms of IRDs. It is necessary to understand the complex interaction among retinal cell types and pathogenic genes by constructing a regulatory network. In this study, we attempt to establish a panoramic expression view of the cooperative work in retinal cells to understand the clinical manifestations and pathogenic bases underlying IRDs.

Methods

Single-cell RNA sequencing (scRNA-seq) data on the retinas from 35 retina samples of 3 species (human, mouse, and zebrafish) including 259,087 cells were adopted to perform a comparative analysis across species. Bioinformatic tools were used to conduct weighted gene co-expression network analysis (WGCNA), single-cell regulatory network analysis, cell-cell communication analysis, and trajectory inference analysis.

Results

The cross-species comparison revealed shared or species-specific gene expression patterns at single-cell resolution, such as the stathmin family genes, which were highly expressed specifically in zebrafish Müller glias (MGs). Thirteen gene modules were identified, of which nine were associated with retinal cell types, and Gene Ontology (GO) enrichment of module genes was consistent with cell-specific highly expressed genes. Many IRD genes were identified as hub genes and cell-specific regulons. Most IRDs, especially the retinitis pigmentosa (RP) genes, were enriched in rod-specific regulons. Integrated expression and transcription regulatory network genes, such as congenital stationary night blindness (CSNB) genes GRK1, PDE6B, and TRPM1, showed cell-specific expression and transcription characteristics in either rods or bipolar cells (BCs). IRD genes showed evolutionary conservation (GNAT2, PDE6G, and SAG) and divergence (GNAT2, MT-ND4, and PDE6A) along the trajectory of photoreceptors (PRs) among species. In particular, the Leber congenital amaurosis (LCA) gene OTX2 showed high expression at the beginning of the trajectory of both PRs and BCs.

Conclusion

We identified molecular pathways and cell types closely connected with IRDs, bridging the gap between gene expression, genetics, and pathogenesis. The IRD genes enriched in cell-specific modules and regulons suggest that these diseases share common etiological bases. Overall, mining of interspecies transcriptome data reveals conserved transcriptomic features of retinas across species and promising applications in both normal retina anatomy and retina pathology.

Single cell multi-omic analysis identifies key genes differentially expressed in innate lymphoid cells from COVID-19 patients

Introduction

Innate lymphoid cells (ILCs) are enriched at mucosal surfaces where they respond rapidly to environmental stimuli and contribute to both tissue inflammation and healing.

Methods

To gain insight into the role of ILCs in the pathology and recovery from COVID-19 infection, we employed a multi-omics approach consisting of Abseq and targeted mRNA sequencing to respectively probe the surface marker expression, transcriptional profile and heterogeneity of ILCs in peripheral blood of patients with COVID-19 compared with healthy controls.

Results

We found that the frequency of ILC1 and ILC2 cells was significantly increased in COVID-19 patients. Moreover, all ILC subsets displayed a significantly higher frequency of CD69-expressing cells, indicating a heightened state of activation. ILC2s from COVID-19 patients had the highest number of significantly differentially expressed (DE) genes. The most notable genes DE in COVID-19 vs healthy participants included a) genes associated with responses to virus infections and b) genes that support ILC self-proliferation, activation and homeostasis. In addition, differential gene regulatory network analysis revealed ILC-specific regulons and their interactions driving the differential gene expression in each ILC.

Discussion

Overall, this study provides mechanistic insights into the characteristics of ILC subsets activated during COVID-19 infection.

DCAF15 control of cohesin dynamics sustains acute myeloid leukemia

The CRL4-DCAF15 E3 ubiquitin ligase complex is targeted by the aryl-sulfonamide molecular glues, leading to neo-substrate recruitment, ubiquitination, and proteasomal degradation. However, the physiological function of DCAF15 remains unknown. Using a domain-focused genetic screening approach, we reveal DCAF15 as an acute myeloid leukemia (AML)-biased dependency. Loss of DCAF15 results in suppression of AML through compromised replication fork integrity and consequent accumulation of DNA damage. Accordingly, DCAF15 loss sensitizes AML to replication stress-inducing therapeutics. Mechanistically, we discover that DCAF15 directly interacts with the SMC1A protein of the cohesin complex and destabilizes the cohesin regulatory factors PDS5A and CDCA5. Loss of PDS5A and CDCA5 removal precludes cohesin acetylation on chromatin, resulting in uncontrolled chromatin loop extrusion, defective DNA replication, and apoptosis. Collectively, our findings uncover an endogenous, cell autonomous function of DCAF15 in sustaining AML proliferation through post-translational control of cohesin dynamics.

SAP30, an oncogenic driver of progression, poor survival, and drug resistance in neuroblastoma

Neuroblastoma is the most devastating extracranial solid malignancy in children. Despite an intense treatment regimen, the prognosis for high-risk neuroblastoma patients remains poor, with less than 40% survival. So far, MYCN amplification status is considered the most prognostic factor but corresponds to only ∼25% of neuroblastoma patients. Therefore, it is essential to identify a better prognosis and therapy response marker in neuroblastoma patients. We applied robust bioinformatic data mining tools, such as weighted gene co-expression network analysis, cisTarget, and single-cell regulatory network inference and clustering on two neuroblastoma patient datasets. We found Sin3A-associated protein 30 (SAP30), a driver transcription factor positively associated with high-risk, progression, stage 4, and poor survival in neuroblastoma patient cohorts. Tumors of high-risk neuroblastoma patients and relapse-specific patient-derived xenografts showed higher SAP30 levels. The advanced pharmacogenomic analysis and CRISPR-Cas9 screens indicated that SAP30 essentiality is associated with cisplatin resistance and further showed higher levels in cisplatin-resistant patient-derived xenograft tumor cell lines. Silencing of SAP30 induced cell death in vitro and led to a reduced tumor burden and size in vivo. Altogether, these results indicate that SAP30 is a better prognostic and cisplatin-resistance marker and thus a potential drug target in high-risk neuroblastoma.

Systematic benchmarking of single-cell ATAC-sequencing protocols

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Integrative Single Cell Atlas Revealed Intratumoral Heterogeneity Generation from an Adaptive Epigenetic Cell State in Human Bladder Urothelial Carcinoma

Intratumor heterogeneity (ITH) of bladder cancer (BLCA) contributes to therapy resistance and immune evasion affecting clinical prognosis. The molecular and cellular mechanisms contributing to BLCA ITH generation remain elusive. It is found that a TM4SF1-positive cancer subpopulation (TPCS) can generate ITH in BLCA, evidenced by integrative single cell atlas analysis. Extensive profiling of the epigenome and transcriptome of all stages of BLCA revealed their evolutionary trajectories. Distinct ancestor cells gave rise to low-grade noninvasive and high-grade invasive BLCA. Epigenome reprograming led to transcriptional heterogeneity in BLCA. During early oncogenesis, epithelial-to-mesenchymal transition generated TPCS. TPCS has stem-cell-like properties and exhibited transcriptional plasticity, priming the development of transcriptionally heterogeneous descendent cell lineages. Moreover, TPCS prevalence in tumor is associated with advanced stage cancer and poor prognosis. The results of this study suggested that bladder cancer interacts with its environment by acquiring a stem cell-like epigenomic landscape, which might generate ITH without additional genetic diversification.

Disruption of maternal vascular remodeling by a fetal endoretrovirus-derived gene in preeclampsia

BACKGROUND

Preeclampsia, one of the most lethal pregnancy-related diseases, is associated with the disruption of uterine spiral artery remodeling during placentation. However, the early molecular events leading to preeclampsia remain unknown.

RESULTS

By analyzing placentas from preeclampsia, non-preeclampsia, and twin pregnancies with selective intrauterine growth restriction, we show that the pathogenesis of preeclampsia is attributed to immature trophoblast and maldeveloped endothelial cells. Delayed epigenetic reprogramming during early extraembryonic tissue development leads to generation of excessive immature trophoblast cells. We find reduction of de novo DNA methylation in these trophoblast cells results in selective overexpression of maternally imprinted genes, including the endoretrovirus-derived gene PEG10 (paternally expressed gene 10). PEG10 forms virus-like particles, which are transferred from the trophoblast to the closely proximate endothelial cells. In normal pregnancy, only a low amount of PEG10 is transferred to maternal cells; however, in preeclampsia, excessive PEG10 disrupts maternal vascular development by inhibiting TGF-beta signaling.

CONCLUSIONS

Our study reveals the intricate epigenetic mechanisms that regulate trans-generational genetic conflict and ultimately ensure proper maternal-fetal interface formation.

Construction of single-cell cross-species chromatin accessibility landscapes with combinatorial-hybridization-based ATAC-seq

Despite recent advances in single-cell genomics, the lack of maps for single-cell candidate cis-regulatory elements (cCREs) in non-mammal species has limited our exploration of conserved regulatory programs across vertebrates and invertebrates. Here, we developed a combinatorial-hybridization-based method for single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) named CH-ATAC-seq, enabling the construction of single-cell accessible chromatin landscapes for zebrafish, Drosophila, and earthworms (Eisenia andrei). By integrating scATAC censuses of humans, monkeys, and mice, we systematically identified 152 distinct main cell types and around 0.8 million cell-type-specific cCREs. Our analysis provided insights into the conservation of neural, muscle, and immune lineages across species, while epithelial cells exhibited a higher organ-origin heterogeneity. Additionally, a large-scale gene regulatory network (GRN) was constructed in four vertebrates by integrating scRNA-seq censuses. Overall, our study provides a valuable resource for comparative epigenomics, identifying the evolutionary conservation and divergence of gene regulation across different species.

An atlas of cells in the human tonsil

Palatine tonsils are secondary lymphoid organs (SLOs) representing the first line of immunological defense against inhaled or ingested pathogens. We generated an atlas of the human tonsil composed of >556,000 cells profiled across five different data modalities, including single-cell transcriptome, epigenome, proteome, and immune repertoire sequencing, as well as spatial transcriptomics. This census identified 121 cell types and states, defined developmental trajectories, and enabled an understanding of the functional units of the tonsil. Exemplarily, we stratified myeloid slan-like subtypes, established a BCL6 enhancer as locally active in follicle-associated T and B cells, and identified SIX5 as putative transcriptional regulator of plasma cell maturation. Analyses of a validation cohort confirmed the presence, annotation, and markers of tonsillar cell types and provided evidence of age-related compositional shifts. We demonstrate the value of this resource by annotating cells from B cell-derived mantle cell lymphomas, linking transcriptional heterogeneity to normal B cell differentiation states of the human tonsil.

Transcriptomic and network analysis identifies shared and unique pathways and immune changes across fibrotic interstitial lung diseases

BACKGROUND

Interstitial lung disease (ILD) encompasses a diverse group of disorders characterized by chronic inflammation and fibrosis of the pulmonary interstitium. Three ILDs, namely idiopathic pulmonary fibrosis (IPF), fibrotic hypersensitivity pneumonitis (fHP), and connective tissue disease-associated ILD (CTD-ILD), exhibit similar progressive fibrosis phenotypes, yet possess distinct etiologies, encouraging us to explore their different underlying mechanisms.

METHODS

Transcriptome data of fibrotic lung tissues from patients with IPF, fHP, and CTD-ILD were subjected to functional annotation, network, and pathway analyses. Additionally, we employed the xCell deconvolution algorithm to predict immune cell infiltration in patients with fibrotic ILDs and healthy controls.

RESULTS

We identified a shared progressive fibrosis-related module in these diseases which was related to extracellular matrix (ECM) degradation and production and potentially regulated by the p53 family transcription factors. In IPF, neuron-related processes emerged as a critical specific mechanism in functional enrichment. In fHP, we observed that B cell signaling and immunoglobulin A (IgA) production may act as predominant processes, which was further verified by B cell infiltration and the central role of CD19 gene. In CTD-ILD, active chemokine processes were enriched, and active dendritic cells (aDCs) were predicted to infiltrate the lung tissues.

CONCLUSIONS

This study revealed shared and specific molecular and cellular pathways among IPF, fHP, and CTD-ILD, providing a basis for understanding their pathogenesis and identifying potential therapeutic targets.

Single-nucleus multi-omic profiling of human placental syncytiotrophoblasts identifies cellular trajectories during pregnancy

The human placenta has a vital role in ensuring a successful pregnancy. Despite the growing body of knowledge about its cellular compositions and functions, there has been limited research on the heterogeneity of the billions of nuclei within the syncytiotrophoblast (STB), a multinucleated entity primarily responsible for placental function. Here we conducted integrated single-nucleus RNA sequencing and single-nucleus ATAC sequencing analyses of human placentas from early and late pregnancy. Our findings demonstrate the dynamic heterogeneity and developmental trajectories of STB nuclei and their correspondence with human trophoblast stem cell (hTSC)-derived STB. Furthermore, we identified transcription factors associated with diverse STB nuclear lineages through their gene regulatory networks and experimentally confirmed their function in hTSC and trophoblast organoid-derived STBs. Together, our data provide insights into the heterogeneity of human STB and represent a valuable resource for interpreting associated pregnancy complications.

Cephalopod-omics: Emerging Fields and Technologies in Cephalopod Biology

Few animal groups can claim the level of wonder that cephalopods instill in the minds of researchers and the general public. Much of cephalopod biology, however, remains unexplored: the largest invertebrate brain, difficult husbandry conditions, and complex (meta-)genomes, among many other things, have hindered progress in addressing key questions. However, recent technological advancements in sequencing, imaging, and genetic manipulation have opened new avenues for exploring the biology of these extraordinary animals. The cephalopod molecular biology community is thus experiencing a large influx of researchers, emerging from different fields, accelerating the pace of research in this clade. In the first post-pandemic event at the Cephalopod International Advisory Council (CIAC) conference in April 2022, over 40 participants from all over the world met and discussed key challenges and perspectives for current cephalopod molecular biology and evolution. Our particular focus was on the fields of comparative and regulatory genomics, gene manipulation, single-cell transcriptomics, metagenomics, and microbial interactions. This article is a result of this joint effort, summarizing the latest insights from these emerging fields, their bottlenecks, and potential solutions. The article highlights the interdisciplinary nature of the cephalopod-omics community and provides an emphasis on continuous consolidation of efforts and collaboration in this rapidly evolving field.

Early-adult methionine restriction reduces methionine sulfoxide and extends lifespan in Drosophila

Methionine restriction (MetR) extends lifespan in various organisms, but its mechanistic understanding remains incomplete. Whether MetR during a specific period of adulthood increases lifespan is not known. In Drosophila, MetR is reported to extend lifespan only when amino acid levels are low. Here, by using an exome-matched holidic medium, we show that decreasing Met levels to 10% extends Drosophila lifespan with or without decreasing total amino acid levels. MetR during the first four weeks of adult life only robustly extends lifespan. MetR in young flies induces the expression of many longevity-related genes, including Methionine sulfoxide reductase A (MsrA), which reduces oxidatively-damaged Met. MsrA induction is foxo-dependent and persists for two weeks after cessation of the MetR diet. Loss of MsrA attenuates lifespan extension by early-adulthood MetR. Our study highlights the age-dependency of the organismal response to specific nutrients and suggests that nutrient restriction during a particular period of life is sufficient for healthspan extension.

A single cell genomics atlas of the Drosophila larval eye reveals distinct photoreceptor developmental timelines

The Drosophila eye is a powerful model system to study the dynamics of cell differentiation, cell state transitions, cell maturation, and pattern formation. However, a high-resolution single cell genomics resource that accurately profiles all major cell types of the larval eye disc and their spatiotemporal relationships is lacking. Here, we report transcriptomic and chromatin accessibility data for all known cell types in the developing eye. Photoreceptors appear as strands of cells that represent their dynamic developmental timelines. As photoreceptor subtypes mature, they appear to assume a common transcriptomic profile that is dominated by genes involved in axon function. We identify cell type maturation genes, enhancers, and potential regulators, as well as genes with distinct R3 or R4 photoreceptor specific expression. Finally, we observe that the chromatin accessibility between cones and photoreceptors is distinct. These single cell genomics atlases will greatly enhance the power of the Drosophila eye as a model system.

Combined inactivation of RB and Hippo converts differentiating Drosophila photoreceptors into eye progenitor cells through derepression of homothorax

The retinoblastoma (RB) and Hippo pathways interact to regulate cell proliferation and differentiation. However, the mechanism of interaction is not fully understood. Drosophila photoreceptors with inactivated RB and Hippo pathways specify normally but fail to maintain their neuronal identity and dedifferentiate. We performed single-cell RNA sequencing to elucidate the cause of dedifferentiation and to determine the fate of these cells. We find that dedifferentiated cells adopt a progenitor-like fate due to inappropriate activation of the retinal differentiation suppressor homothorax (hth) by Yki/Sd. This results in the activation of a distinct Yki/Hth transcriptional program, driving photoreceptor dedifferentiation. We show that Rbf physically interacts with Yki and, together with the GAGA factor, inhibits the hth expression. Thus, RB and Hippo pathways cooperate to maintain photoreceptor differentiation by preventing inappropriate expression of hth in differentiating photoreceptors. Our work highlights the importance of both RB and Hippo pathway activities for maintaining the state of terminal differentiation.

Smurfness-based two-phase model of ageing helps deconvolve the ageing transcriptional signature

Ageing is characterised at the molecular level by six transcriptional 'hallmarks of ageing', that are commonly described as progressively affected as time passes. By contrast, the 'Smurf' assay separates high-and-constant-mortality risk individuals from healthy, zero-mortality risk individuals, based on increased intestinal permeability. Performing whole body total RNA sequencing, we found that Smurfness distinguishes transcriptional changes associated with chronological age from those associated with biological age. We show that transcriptional heterogeneity increases with chronological age in non-Smurf individuals preceding the other five hallmarks of ageing that are specifically associated with the Smurf state. Using this approach, we also devise targeted pro-longevity genetic interventions delaying entry in the Smurf state. We anticipate that increased attention to the evolutionary conserved Smurf phenotype will bring about significant advances in our understanding of the mechanisms of ageing.

Activation of the transcription factor NFAT5 in the tumor microenvironment enforces CD8+ T cell exhaustion

Persistent exposure to antigen during chronic infection or cancer renders T cells dysfunctional. The molecular mechanisms regulating this state of exhaustion are thought to be common in infection and cancer, despite obvious differences in their microenvironments. Here we found that NFAT5, an NFAT family transcription factor that lacks an AP-1 docking site, was highly expressed in exhausted CD8+ T cells in the context of chronic infections and tumors but was selectively required in tumor-induced CD8+ T cell exhaustion. Overexpression of NFAT5 in CD8+ T cells reduced tumor control, while deletion of NFAT5 improved tumor control by promoting the accumulation of tumor-specific CD8+ T cells that had reduced expression of the exhaustion-associated proteins TOX and PD-1 and produced more cytokines, such as IFNɣ and TNF, than cells with wild-type levels of NFAT5, specifically in the precursor exhausted PD-1+TCF1+TIM-3-CD8+ T cell population. NFAT5 did not promote T cell exhaustion during chronic infection with clone 13 of lymphocytic choriomeningitis virus. Expression of NFAT5 was induced by TCR triggering, but its transcriptional activity was specific to the tumor microenvironment and required hyperosmolarity. Thus, NFAT5 promoted the exhaustion of CD8+ T cells in a tumor-selective fashion.

Integrated immunological analysis of single-cell and bulky tissue transcriptomes reveals the role of interactions between M0 macrophages and naïve CD4+ T cells in the immunosuppressive microenvironment of cervical cancer

In recent decades, the incidence and mortality of cervical cancer have declined in developed countries due to the implementation of screening and vaccination programs. However, cervical cancer remains one of the major culprits of cancer-related deaths in young women. Current studies have found that immune cell-related intercellular communication in the tumor microenvironment has a large impact on the construction of the immunosuppressive microenvironment. In this study, we performed a comprehensive immune analysis on bulk RNA-seq and scRNA-seq data obtained from cervical cancer and revealed that two highly plastic cell populations, M0 macrophages and naïve CD4+ T cells, were significantly correlated with prognosis and clinical phenotypes. Notably, signaling between M0 macrophages and naïve CD4+ T cells as well as intracellular transcription factor activity were significantly altered in the tumor state. Furthermore, we identified overlapping genes between the transcription factor target genes of M0 macrophages or naïve CD4+ T cells and the differentially expressed genes in each type of cell, and these overlapping genes were subsequently subjected to an analysis using the LASSO regression model. Finally, we generated a score index that was significantly associated with the clinical prognosis of cervical cancer. In conclusion, interventions to improve the communication between M0 macrophages and naïve CD4+ T cells may help to improve the immunosuppressive microenvironment of cervical cancer and prevent immune evasion. The relevant molecular mechanisms need to be further validated by experimental and cohort studies.

Single-cell dissection of Merkel cell carcinoma heterogeneity unveils transcriptomic plasticity and therapeutic vulnerabilities

Merkel cell carcinoma (MCC), a rare but aggressive skin cancer, remains a challenge in the era of precision medicine. Immune checkpoint inhibitors (ICIs), the only approved therapy for advanced MCC, are impeded by high primary and acquired resistance. Hence, we dissect transcriptomic heterogeneity at single-cell resolution in a panel of patient tumors, revealing phenotypic plasticity in a subset of treatment-naive MCC. The tumor cells in a "mesenchymal-like" state are endowed with an inflamed phenotype that portends a better ICI response. This observation is also validated in the largest whole transcriptomic dataset available from MCC patient tumors. In contrast, ICI-resistant tumors predominantly express neuroepithelial markers in a well-differentiated state with "immune-cold" landscape. Importantly, a subtle shift to "mesenchymal-like" state reverts copanlisib resistance in primary MCC cells, highlighting potential strategies in patient stratification for therapeutics to harness tumor cell plasticity, augment treatment efficacy, and avert resistance.

Damage sensing mediated by serine proteases Hayan and Persephone for Toll pathway activation in apoptosis-deficient flies

The mechanisms by which the innate immune system senses damage have been extensively explored in multicellular organisms. In Drosophila, various types of tissue damage, including epidermal injury, tumor formation, cell competition, and apoptosis deficiency, induce sterile activation of the Toll pathway, a process that requires the use of extracellular serine protease (SP) cascades. Upon infection, the SP Spätzle (Spz)-processing enzyme (SPE) cleaves and activates the Toll ligand Spz downstream of two paralogous SPs, Hayan and Persephone (Psh). However, upon tissue damage, it is not fully understood which SPs establish Spz activation cascades nor what damage-associated molecules can activate SPs. In this study, using newly generated uncleavable spz mutant flies, we revealed that Spz cleavage is required for the sterile activation of the Toll pathway, which is induced by apoptosis-deficient damage of wing epidermal cells in adult Drosophila. Proteomic analysis of hemolymph, followed by experiments with Drosophila Schneider 2 (S2) cells, revealed that among hemolymph SPs, both SPE and Melanization Protease 1 (MP1) have high capacities to cleave Spz. Additionally, in S2 cells, MP1 acts downstream of Hayan and Psh in a similar manner to SPE. Using genetic analysis, we found that the upstream SPs Hayan and Psh contributes to the sterile activation of the Toll pathway. While SPE/MP1 double mutants show more impairment of Toll activation upon infection than SPE single mutants, Toll activation is not eliminated in these apoptosis-deficient flies. This suggests that Hayan and Psh sense necrotic damage, inducing Spz cleavage by SPs other than SPE and MP1. Furthermore, hydrogen peroxide, a representative damage-associated molecule, activates the Psh-Spz cascade in S2 cells overexpressing Psh. Considering that reactive oxygen species (ROS) were detected in apoptosis-deficient wings, our findings highlight the importance of ROS as signaling molecules that induce the activation of SPs such as Psh in response to damage.

LincRNA-Cox2 Regulates Smoke-induced Inflammation in Murine Macrophages

Cigarette smoke (CS) exposure is a risk factor for many chronic diseases, including chronic obstructive pulmonary disease, but the mechanism by which smoke exposure can alter homeostasis and bring about chronic inflammation is poorly understood. Here, we showcase a novel role for smoke in regulating long noncoding RNAs, showing that it activates lincRNA-Cox2, which we previously characterized as functional in inflammatory regulation. Exposing lincRNA-Cox2 murine models to smoke in vivo confirmed lincRNA-Cox2 as a regulator of inflammatory gene expression in response to smoke both systemically and within the lung. We also report that lincRNA-Cox2 negatively regulates genes in smoked bone marrow-derived macrophages exposed to LPS stimulation. In addition to the effects on long noncoding RNAs, we also report dysregulated transcription and splicing of inflammatory protein-coding genes in the bone marrow niche after CS exposure in vivo. Collectively, this work provides insights into how innate immune signaling from gene expression to splicing is altered after in vivo exposure to CS and highlights an important new role for lincRNA-Cox2 in regulating immune genes after smoke exposure.

Orthogonal cytokine engineering enables novel synthetic effector states escaping canonical exhaustion in tumor-rejecting CD8+ T cells

To date, no immunotherapy approaches have managed to fully overcome T-cell exhaustion, which remains a mandatory fate for chronically activated effector cells and a major therapeutic challenge. Understanding how to reprogram CD8+ tumor-infiltrating lymphocytes away from exhausted effector states remains an elusive goal. Our work provides evidence that orthogonal gene engineering of T cells to secrete an interleukin (IL)-2 variant binding the IL-2Rβγ receptor and the alarmin IL-33 reprogrammed adoptively transferred T cells to acquire a novel, synthetic effector state, which deviated from canonical exhaustion and displayed superior effector functions. These cells successfully overcame homeostatic barriers in the host and led-in the absence of lymphodepletion or exogenous cytokine support-to high levels of engraftment and tumor regression. Our work unlocks a new opportunity of rationally engineering synthetic CD8+ T-cell states endowed with the ability to avoid exhaustion and control advanced solid tumors.

Combined inactivation of RB and Hippo pathways converts differentiating photoreceptors into eye progenitor cells through derepression of homothorax

The RB and Hippo pathways interact to regulate cell proliferation and differentiation. However, their mechanism of interaction is not fully understood. Drosophila photoreceptors with inactivated RB and Hippo pathways specify normally but fail to maintain neuronal identity and dedifferentiate. We performed single-cell RNA-sequencing to elucidate the cause of dedifferentiation and the fate of these cells. We find that dedifferentiated cells adopt a progenitor-like fate due to inappropriate activation of the retinal differentiation suppressor homothorax (hth) by Yki/Sd. This results in activation of the Yki/Hth transcriptional program, driving photoreceptor dedifferentiation. We show that Rbf physically interacts with Yki which, together with the GAGA factor, inhibits hth expression. Thus, RB and Hippo pathways cooperate to maintain photoreceptor differentiation by preventing inappropriate expression of hth in differentiating photoreceptors. Our work accentuates the importance of both RB and Hippo pathway activity for maintaining the state of terminal differentiation.

Loss of Kmt2c in vivo leads to EMT, mitochondrial dysfunction and improved response to lapatinib in breast cancer

Deep sequencing of human tumours has uncovered a previously unappreciated role for epigenetic regulators in tumorigenesis. H3K4 methyltransferase KMT2C/MLL3 is mutated in several solid malignancies, including more than 10% of breast tumours. To study the tumour suppressor role of KMT2C in breast cancer, we generated mouse models of Erbb2/Neu, Myc or PIK3CA-driven tumorigenesis, in which the Kmt2c locus is knocked out specifically in the luminal lineage of mouse mammary glands using the Cre recombinase. Kmt2c knock out mice develop tumours earlier, irrespective of the oncogene, assigning a bona fide tumour suppressor role for KMT2C in mammary tumorigenesis. Loss of Kmt2c induces extensive epigenetic and transcriptional changes, which lead to increased ERK1/2 activity, extracellular matrix re-organization, epithelial-to-mesenchymal transition and mitochondrial dysfunction, the latter associated with increased reactive oxygen species production. Loss of Kmt2c renders the Erbb2/Neu-driven tumours more responsive to lapatinib. Publicly available clinical datasets revealed an association of low Kmt2c gene expression and better long-term outcome. Collectively, our findings solidify the role of KMT2C as a tumour suppressor in breast cancer and identify dependencies that could be therapeutically amenable.

Atf4 regulates angiogenic differences between alveolar bone and long bone macrophages by regulating M1 polarization, based on single-cell RNA sequencing, RNA-seq and ATAC-seq analysis

In the repair of maxillofacial bone defects, autogenous craniofacial bone can often provide superior clinical results over long bone grafts. Most current studies have focused on the osteogenic differences between alveolar bone marrow (ABM) and long bone marrow (LBM), however, studies about the angiogenic differences between the two are currently lacking. We downloaded single-cell RNA sequencing (scRNA-seq) of mouse ABM and LBM respectively from the public database, and the data were processed by using Seurat package. CellphoneDB2 results showed that macrophages had the strongest interaction with mesenchymal stem cells (MSCs) and endothelial cells (ECs). ELISA results confirmed that ABM macrophages secreted a higher level of vascular endothelial growth factor A (Vegfa) compared to LBM macrophages, which further promoted angiogenesis of ECs and MSCs. Using SCENIC package, six key transcription factors (TFs) were identified to regulate the difference between ABM and LBM macrophages, and activating transcription factor 4 (Atf4) was confirmed to be more expressed in ABM macrophages by polymerase chain reaction (PCR) and western blot (WB), with predicted target genes including Vegfa. Besides, the result of scRNA-seq implied ABM macrophages more in M1 status than LBM macrophages, which was confirmed by the following experiments. From the results of another assay for transposase accessible chromatin sequencing (ATAC-seq) and RNA-seq about M1 macrophages, Atf4 was also confirmed to regulate the M1 polarization. So, we suspected that Atf4 regulated the different expression of Vegfa between ABM and LBM macrophages by activating M1 polarization. After knocking down Atf4, the expression of M1 polarization markers and Vegfa were downregulated and vasculogenic differences were eliminated, which were subsequently reversed by the addition of LPS/IFN-γ. Our study might provide a new idea to improve the success rate of autologous bone grafting and treatment of oral diseases.

Integrated bioinformatics analysis identifies shared immune changes between ischemic stroke and COVID 19

Although COVID-19 is primarily a respiratory disease, its neurological complications, such as ischemic stroke (IS), have aroused growing concerns and reports. However, the molecular mechanisms that underlie IS and COVID-19 are not well understood. Therefore, we implemented transcriptomic analysis from eight GEO datasets consist of 1191 samples to detect common pathways and molecular biomarkers in IS and COVID-19 that help understand the linkage between them. Differentially expressed genes (DEGs) were detected for IS and COVID-19 separately for finding shared mechanisms and we found that immune-related pathways were outlined with statistical significance. JAK2, which was identified as a hub gene, was supposed to be a potential therapeutic gene targets during the immunological process of COVID-19 and IS. Besides, we found a decrease in the proportion of CD8+ T and T helper 2 cells in the peripheral circulation of both COVID and IS patients, and NCR3 expression was significantly correlated with this change. In conclusion, we demonstrated that transcriptomic analyses reported in this study could make a deeper understanding of the common mechanism and might be promising for effective therapeutic for IS and COVID-19.

Endothelial-to-osteoblast transition in normal mouse bone development

Metastatic prostate cancer (PCa) in bone induces bone-forming lesions. We have previously shown that PCa-induced bone originates from endothelial cells (ECs) that have undergone EC-to-osteoblast (OSB) transition. Here, we investigated whether EC-to-OSB transition also occurs during normal bone formation. We developed an EC and OSB dual-color reporter mouse (DRM) model that marks EC-OSB hybrid cells with red and green fluorescent proteins. We observed EC-to-OSB transition (RFP and GFP co-expression) in both endochondral and intramembranous bone formation during embryonic development and in adults. Co-expression was confirmed in cells isolated from DRM. Bone marrow- and lung-derived ECs underwent transition to OSBs and mineralization in osteogenic medium. RNA-sequencing revealed GATA family transcription factors were upregulated in EC-OSB hybrid cells and knockdown of GATA3 inhibited BMP4-induced mineralization. Our findings support that EC-to-OSB transition occurs during normal bone development and suggest a new paradigm regarding the endothelial origin of OSBs.

Senescent cells perturb intestinal stem cell differentiation through Ptk7 induced noncanonical Wnt and YAP signaling

Cellular senescence and the senescence-associated secretory phenotype (SASP) are implicated in aging and age-related disease, and SASP-related inflammation is thought to contribute to tissue dysfunction in aging and diseased animals. However, whether and how SASP factors influence the regenerative capacity of tissues remains unclear. Here, using intestinal organoids as a model of tissue regeneration, we show that SASP factors released by senescent fibroblasts deregulate stem cell activity and differentiation and ultimately impair crypt formation. We identify the secreted N-terminal domain of Ptk7 as a key component of the SASP that activates non-canonical Wnt / Ca2+ signaling through FZD7 in intestinal stem cells (ISCs). Changes in cytosolic [Ca2+] elicited by Ptk7 promote nuclear translocation of YAP and induce expression of YAP/TEAD target genes, impairing symmetry breaking and stem cell differentiation. Our study discovers secreted Ptk7 as a factor released by senescent cells and provides insight into the mechanism by which cellular senescence contributes to tissue dysfunction in aging and disease.

Immunoediting instructs tumor metabolic reprogramming to support immune evasion

Immunoediting sculpts immunogenicity and thwarts host anti-tumor responses in tumor cells during tumorigenesis; however, it remains unknown whether metabolic programming of tumor cells can be guided by immunosurveillance. Here, we report that T cell-mediated immunosurveillance in early-stage tumorigenesis instructs c-Myc upregulation and metabolic reprogramming in tumor cells. This previously unexplored tumor-immune interaction is controlled by non-canonical interferon gamma (IFNγ)-STAT3 signaling and supports tumor immune evasion. Our findings uncover that immunoediting instructs deregulated bioenergetic programs in tumor cells to empower them to disarm the T cell-mediated immunosurveillance by imposing metabolic tug-of-war between tumor and infiltrating T cells and forming the suppressive tumor microenvironment.

Integrated analysis of mRNA-seq and miRNA-seq reveals the potential roles of Egr1, Rxra and Max in kidney stone disease

Nephrolithiasis is one of the most common and frequent urologic diseases worldwide. The molecular mechanism of kidney stone formation is complex and remains to be illustrated. Transcript factors (TFs) that influenced the expression pattern of multiple genes, as well as microRNAs, important posttranscriptional modulators, play vital roles in this disease progression. Datasets of nephrolithiasis mice and kidney stone patients were acquired from Gene Expression Omnibus repository. TFs were predicted from differentially expressed genes by RcisTarget. The target genes of differential-expressed microRNAs were predicted by miRWalk. MicroRNA-mRNA network and PPI network were constructed. Functional enrichment analysis was performed via Metascape and Cytoscape identified hub genes. The assay of quantitative real-time PCR (q-PCR) and immunochemistry and the datasets of oxalate diet-induced nephrolithiasis mice kidneys and kidney stone patients' samples were utilized to validate the bioinformatic results. We identified three potential key TFs (Egr1, Rxra, Max), which can be modulated by miR-181a-5p, miR-7b-3p and miR-22-3p, respectively. The TFs and their regulated hub genes influenced the progression of nephrolithiasis via altering the expression of genes enriched in the functions of fibrosis, cell proliferation and molecular transportation and metabolism. The expression changes of transcription factors were consistent in q-PCR and immunochemistry results. For regulated hub genes, they showed consistent expression changes in oxalate diet-induced nephrolithiasis mice model and human kidneys with stones. The identified and verified three TFs, which may be modulated by microRNAs in nephrolithiasis disease progression, mainly influence biological processes responding to fibrosis, proliferation and molecular transportation and metabolism. The transcript influence showed consistency in multiple nephrolithiasis mice models and kidney stone patients.

Epigenetic Regulation of Profibrotic Macrophages in Systemic Sclerosis-Associated Interstitial Lung Disease

OBJECTIVE

Systemic sclerosis-associated interstitial lung disease (SSc-ILD) is the leading cause of death in patients with SSc with unclear pathogenesis and limited treatment options. Evidence strongly supports an important role for profibrotic secreted phosphoprotein 1 (SPP1)-expressing macrophages in SSc-ILD. This study was undertaken to define the transcriptome and chromatin structural changes of SPP1 SSc-ILD macrophages in order to better understand their role in promoting fibrosis and to identify transcription factors associated with open chromatin driving their altered phenotype.

METHODS

We performed single-cell RNA sequencing (scRNA-Seq) on 11 explanted SSc-ILD and healthy control lung samples, as well as single-cell assay for transposase-accessible chromatin sequencing on 5 lung samples to define altered chromatin accessibility of SPP1 macrophages. We predicted transcription factors regulating SPP1 macrophages using single-cell regulatory network inference and clustering (SCENIC) and determined transcription factor binding sites associated with global alterations in SPP1 chromatin accessibility using Signac/Seurat.

RESULTS

We identified distinct macrophage subpopulations using scRNA-Seq analysis in healthy and SSc-ILD lungs and assessed gene expression changes during the change of healthy control macrophages into SPP1 macrophages. Analysis of open chromatin validated SCENIC predictions, indicating that microphthalmia-associated transcription factor, transcription factor EB, activating transcription factor 6, sterol regulatory element binding transcription factor 1, basic helix-loop-helix family member E40, Kruppel-like factor 6, ETS variant transcription factor 5, and/or members of the activator protein 1 family of transcription factors regulate SPP1 macrophage differentiation.

CONCLUSION

Our findings shed light on the underlying changes in chromatin structure and transcription factor regulation of profibrotic SPP1 macrophages in SSc-ILD. Similar alterations in SPP1 macrophages may underpin fibrosis in other organs involved in SSc and point to novel targets for the treatment of SSc-ILD, specifically targeting profibrotic macrophages.

Hippo signaling instructs ectopic but not normal organ growth

The Hippo signaling pathway is widely considered a master regulator of organ growth because of the prominent overgrowth phenotypes caused by experimental manipulation of its activity. Contrary to this model, we show here that removing Hippo transcriptional output did not impair the ability of the mouse liver and Drosophila eyes to grow to their normal size. Moreover, the transcriptional activity of the Hippo pathway effectors Yap/Taz/Yki did not correlate with cell proliferation, and hyperactivation of these effectors induced gene expression programs that did not recapitulate normal development. Concordantly, a functional screen in Drosophila identified several Hippo pathway target genes that were required for ectopic overgrowth but not normal growth. Thus, Hippo signaling does not instruct normal growth, and the Hippo-induced overgrowth phenotypes are caused by the activation of abnormal genetic programs.

Genomics analysis of Drosophila sechellia response to Morinda citrifolia fruit diet

Drosophila sechellia is an island endemic host specialist that has evolved to consume the toxic fruit of Morinda citrifolia, also known as noni fruit. Recent studies by our group and others have examined genome-wide gene expression responses of fruit flies to individual highly abundant compounds found in noni responsible for the fruit's unique chemistry and toxicity. In order to relate these reductionist experiments to the gene expression responses to feeding on noni fruit itself, we fed rotten noni fruit to adult female D. sechellia and performed RNA-sequencing. Combining the reductionist and more wholistic approaches, we have identified candidate genes that may contribute to each individual compound and those that play a more general role in response to the fruit as a whole. Using the compound specific and general responses, we used transcription factor prediction analyses to identify the regulatory networks and specific regulators involved in the responses to each compound and the fruit itself. The identified genes and regulators represent the possible genetic mechanisms and biochemical pathways that contribute to toxin resistance and noni specialization in D. sechellia.

Single-cell multimodal analysis identifies common regulatory programs in synovial fibroblasts of rheumatoid arthritis patients and modeled TNF-driven arthritis

BACKGROUND

Synovial fibroblasts (SFs) are specialized cells of the synovium that provide nutrients and lubricants for the proper function of diarthrodial joints. Recent evidence appreciates the contribution of SF heterogeneity in arthritic pathologies. However, the normal SF profiles and the molecular networks that govern the transition from homeostatic to arthritic SF heterogeneity remain poorly defined.

METHODS

We applied a combined analysis of single-cell (sc) transcriptomes and epigenomes (scRNA-seq and scATAC-seq) to SFs derived from naïve and hTNFtg mice (mice that overexpress human TNF, a murine model for rheumatoid arthritis), by employing the Seurat and ArchR packages. To identify the cellular differentiation lineages, we conducted velocity and trajectory analysis by combining state-of-the-art algorithms including scVelo, Slingshot, and PAGA. We integrated the transcriptomic and epigenomic data to infer gene regulatory networks using ArchR and custom-implemented algorithms. We performed a canonical correlation analysis-based integration of murine data with publicly available datasets from SFs of rheumatoid arthritis patients and sought to identify conserved gene regulatory networks by utilizing the SCENIC algorithm in the human arthritic scRNA-seq atlas.

RESULTS

By comparing SFs from healthy and hTNFtg mice, we revealed seven homeostatic and two disease-specific subsets of SFs. In healthy synovium, SFs function towards chondro- and osteogenesis, tissue repair, and immune surveillance. The development of arthritis leads to shrinkage of homeostatic SFs and favors the emergence of SF profiles marked by Dkk3 and Lrrc15 expression, functioning towards enhanced inflammatory responses and matrix catabolic processes. Lineage inference analysis indicated that specific Thy1+ SFs at the root of trajectories lead to the intermediate Thy1+/Dkk3+/Lrrc15+ SF states and culminate in a destructive and inflammatory Thy1- SF identity. We further uncovered epigenetically primed gene programs driving the expansion of these arthritic SFs, regulated by NFkB and new candidates, such as Runx1. Cross-species analysis of human/mouse arthritic SF data determined conserved regulatory and transcriptional networks.

CONCLUSIONS

We revealed a dynamic SF landscape from health to arthritis providing a functional genomic blueprint to understand the joint pathophysiology and highlight the fibroblast-oriented therapeutic targets for combating chronic inflammatory and destructive arthritic disease.

Sensing of the non-essential amino acid tyrosine governs the response to protein restriction in Drosophila

The intake of dietary protein regulates growth, metabolism, fecundity and lifespan across various species, which makes amino acid (AA)-sensing vital for adaptation to the nutritional environment. The general control nonderepressible 2 (GCN2)-activating transcription factor 4 (ATF4) pathway and the mechanistic target of rapamycin complex 1 (mTORC1) pathway are involved in AA-sensing. However, it is not fully understood which AAs regulate these two pathways in living animals and how they coordinate responses to protein restriction. Here we show in Drosophila that the non-essential AA tyrosine (Tyr) is a nutritional cue in the fat body necessary and sufficient for promoting adaptive responses to a low-protein diet, which entails reduction of protein synthesis and mTORC1 activity and increased food intake. This adaptation is regulated by dietary Tyr through GCN2-independent induction of ATF4 target genes in the fat body. This study identifies the Tyr-ATF4 axis as a regulator of the physiological response to a low-protein diet and sheds light on the essential function of a non-essential nutrient.

Heterochronic parabiosis induces stem cell revitalization and systemic rejuvenation across aged tissues

The young circulatory milieu capable of delaying aging in individual tissues is of interest as rejuvenation strategies, but how it achieves cellular- and systemic-level effects has remained unclear. Here, we constructed a single-cell transcriptomic atlas across aged tissues/organs and their rejuvenation in heterochronic parabiosis (HP), a classical model to study systemic aging. In general, HP rejuvenated adult stem cells and their niches across tissues. In particular, we identified hematopoietic stem and progenitor cells (HSPCs) as one of the most responsive cell types to young blood exposure, from which a continuum of cell state changes across the hematopoietic and immune system emanated, through the restoration of a youthful transcriptional regulatory program and cytokine-mediated cell-cell communications in HSPCs. Moreover, the reintroduction of the identified rejuvenating factors alleviated age-associated lymphopoiesis decline. Overall, we provide comprehensive frameworks to explore aging and rejuvenating trajectories at single-cell resolution and revealed cellular and molecular programs that instruct systemic revitalization by blood-borne factors.

ASC proneural factors are necessary for chromatin remodeling during neuroectodermal to neuroblast fate transition to ensure the timely initiation of the neural stem cell program

Background

In both Drosophila and mammals, the achaete-scute (ASC/ASCL) proneural bHLH transcription factors are expressed in the developing central and peripheral nervous systems, where they function during specification and maintenance of the neural stem cells in opposition to Notch signaling. In addition to their role in nervous system development, ASC transcription factors are oncogenic and exhibit chromatin reprogramming activity; however, the impact of ASC on chromatin dynamics during neural stem cell generation remains elusive. Here, we investigate the chromatin changes accompanying neural commitment using an integrative genetics and genomics methodology.

Results

We found that ASC factors bind equally strongly to two distinct classes of cis-regulatory elements: open regions remodeled earlier during maternal to zygotic transition by Zelda and less accessible, Zelda-independent regions. Both classes of cis-elements exhibit enhanced chromatin accessibility during neural specification and correlate with transcriptional regulation of genes involved in a variety of biological processes necessary for neuroblast function/homeostasis. We identified an ASC-Notch regulated TF network that includes likely prime regulators of neuroblast function. Using a cohort of ASC target genes, we report that ASC null neuroblasts are defectively specified, remaining initially stalled, unable to divide, and lacking expression of many proneural targets. When mutant neuroblasts eventually start proliferating, they produce compromised progeny. Reporter lines driven by proneural-bound enhancers display ASC dependency, suggesting that the partial neuroblast identity seen in the absence of ASC genes is likely driven by other, proneural-independent, cis-elements. Neuroblast impairment and the late differentiation defects of ASC mutants are corrected by ectodermal induction of individual ASC genes but not by individual members of the TF network downstream of ASC. However, in wild-type embryos, the induction of individual members of this network induces CNS hyperplasia, suggesting that they synergize with the activating function of ASC to consolidate the chromatin dynamics that promote neural specification.

Conclusions

We demonstrate that ASC proneural transcription factors are indispensable for the timely initiation of the neural stem cell program at the chromatin level by regulating a large number of enhancers in the vicinity of neural genes. This early chromatin remodeling is crucial for both neuroblast homeostasis as well as future progeny fidelity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01300-8.

The Blimp-1 transcription factor acts in non-neuronal cells to regulate terminal differentiation of the Drosophila eye

The formation of a functional organ such as the eye requires specification of the correct cell types and their terminal differentiation into cells with the appropriate morphologies and functions. Here, we show that the zinc-finger transcription factor Blimp-1 acts in secondary and tertiary pigment cells in the Drosophila retina to promote the formation of a bi-convex corneal lens with normal refractive power, and in cone cells to enable complete extension of the photoreceptor rhabdomeres. Blimp-1 expression depends on the hormone ecdysone, and loss of ecdysone signaling causes similar differentiation defects. Timely termination of Blimp-1 expression is also important, as its overexpression in the eye has deleterious effects. Our transcriptomic analysis revealed that Blimp-1 regulates the expression of many structural and secreted proteins in the retina. Blimp-1 may function in part by repressing another transcription factor; Slow border cells is highly upregulated in the absence of Blimp-1, and its overexpression reproduces many of the effects of removing Blimp-1. This work provides insight into the transcriptional networks and cellular interactions that produce the structures necessary for visual function.

Characterization of Accessible Chromatin Regions in Cattle Rumen Epithelial Tissue during Weaning

Weaning in ruminants is characterized by the transition from a milk-based diet to a solid diet, which drives a critical gastrointestinal tract transformation. Understanding the regulatory control of this transformation during weaning can help to identify strategies to improve rumen health. This study aimed to identify regions of accessible chromatin in rumen epithelial tissue in pre- and post-weaning calves and investigate differentially accessible regions (DARs) to uncover regulatory elements in cattle rumen development using the ATAC-seq approach. A total of 126,071 peaks were identified, covering 1.15% of the cattle genome. From these accessible regions, 2766 DARs were discovered. Gene ontology enrichment resulted in GO terms related to the cell adhesion, anchoring junction, growth, cell migration, motility, and morphogenesis. In addition, putative regulatory canonical pathways were identified (TGFβ, integrin-linked kinase, integrin signaling, and regulation of the epithelial–mesenchymal transition). Canonical pathways integrated with co-expression results showed that TGFβ and ILK signaling pathways play essential roles in rumen development through the regulation of cellular adhesions. In this study, DARs during weaning were identified, revealing enhancers, transcription factors, and candidate target genes that represent potential biomarkers for the bovine rumen development, which will serve as a molecular tool for rumen development studies.

ETV6-NCOA2 fusion induces T/myeloid mixed-phenotype leukemia through transformation of nonthymic hematopoietic progenitor cells

Mixed-phenotype acute leukemia is a rare subtype of leukemia in which both myeloid and lymphoid markers are co-expressed on the same malignant cells. The pathogenesis is largely unknown, and the treatment is challenging. We previously reported the specific association of the recurrent t(8;12)(q13;p13) chromosomal translocation that creates the ETV6-NCOA2 fusion with T/myeloid leukemias. Here we report that ETV6-NCOA2 initiates T/myeloid leukemia in preclinical models; ectopic expression of ETV6-NCOA2 in mouse bone marrow hematopoietic progenitors induced T/myeloid lymphoma accompanied by spontaneous Notch1-activating mutations. Similarly, cotransduction of human cord blood CD34+ progenitors with ETV6-NCOA2 and a nontransforming NOTCH1 mutant induced T/myeloid leukemia in immunodeficient mice; the immunophenotype and gene expression pattern were similar to those of patient-derived ETV6-NCOA2 leukemias. Mechanistically, we show that ETV6-NCOA2 forms a transcriptional complex with ETV6 and the histone acetyltransferase p300, leading to derepression of ETV6 target genes. The expression of ETV6-NCOA2 in human and mouse nonthymic hematopoietic progenitor cells induces transcriptional dysregulation, which activates a lymphoid program while failing to repress the expression of myeloid genes such as CSF1 and MEF2C. The ETV6-NCOA2 induced arrest at an early immature T-cell developmental stage. The additional acquisition of activating NOTCH1 mutations transforms the early immature ETV6-NCOA2 cells into T/myeloid leukemias. Here, we describe the first preclinical model to depict the initiation of T/myeloid leukemia by a specific somatic genetic aberration.

JASPAR 2022: the 9th release of the open-access database of transcription factor binding profiles

JASPAR (http://jaspar.genereg.net/) is an open-access database containing manually curated, non-redundant transcription factor (TF) binding profiles for TFs across six taxonomic groups. In this 9th release, we expanded the CORE collection with 341 new profiles (148 for plants, 101 for vertebrates, 85 for urochordates, and 7 for insects), which corresponds to a 19% expansion over the previous release. We added 298 new profiles to the Unvalidated collection when no orthogonal evidence was found in the literature. All the profiles were clustered to provide familial binding profiles for each taxonomic group. Moreover, we revised the structural classification of DNA binding domains to consider plant-specific TFs. This release introduces word clouds to represent the scientific knowledge associated with each TF. We updated the genome tracks of TFBSs predicted with JASPAR profiles in eight organisms; the human and mouse TFBS predictions can be visualized as native tracks in the UCSC Genome Browser. Finally, we provide a new tool to perform JASPAR TFBS enrichment analysis in user-provided genomic regions. All the data is accessible through the JASPAR website, its associated RESTful API, the R/Bioconductor data package, and a new Python package, pyJASPAR, that facilitates serverless access to the data.

Multilayer omics analysis reveals a non-classical retinoic acid signaling axis that regulates hematopoietic stem cell identity

Hematopoietic stem cells (HSCs) rely on complex regulatory networks to preserve stemness. Due to the scarcity of HSCs, technical challenges have limited our insights into the interplay between metabolites, transcription, and the epigenome. In this study, we generated low-input metabolomics, transcriptomics, chromatin accessibility, and chromatin immunoprecipitation data, revealing distinct metabolic hubs that are enriched in HSCs and their downstream multipotent progenitors. Mechanistically, we uncover a non-classical retinoic acid (RA) signaling axis that regulates HSC function. We show that HSCs rely on Cyp26b1, an enzyme conventionally considered to limit RA effects in the cell. In contrast to the traditional view, we demonstrate that Cyp26b1 is indispensable for production of the active metabolite 4-oxo-RA. Further, RA receptor beta (Rarb) is required for complete transmission of 4-oxo-RA-mediated signaling to maintain stem cells. Our findings emphasize that a single metabolite controls stem cell fate by instructing epigenetic and transcriptional attributes.

Decoding gene regulation in the fly brain

The Drosophila brain is a frequently used model in neuroscience. Single-cell transcriptome analysis^1-6, three-dimensional morphological classification⁷ and electron microscopy mapping of the connectome^8,9 have revealed an immense diversity of neuronal and glial cell types that underlie an array of functional and behavioural traits in the fly. The identities of these cell types are controlled by gene regulatory networks (GRNs), involving combinations of transcription factors that bind to genomic enhancers to regulate their target genes. Here, to characterize GRNs at the cell-type level in the fly brain, we profiled the chromatin accessibility of 240,919 single cells spanning 9 developmental timepoints and integrated these data with single-cell transcriptomes. We identify more than 95,000 regulatory regions that are used in different neuronal cell types, of which 70,000 are linked to developmental trajectories involving neurogenesis, reprogramming and maturation. For 40 cell types, uniquely accessible regions were associated with their expressed transcription factors and downstream target genes through a combination of motif discovery, network inference and deep learning, creating enhancer GRNs. The enhancer architectures revealed by DeepFlyBrain lead to a better understanding of neuronal regulatory diversity and can be used to design genetic driver lines for cell types at specific timepoints, facilitating their characterization and manipulation.

OUP accepted manuscript

Mitochondria are subcellular organelles present in almost all eukaryotic cells, which play a central role in cellular metabolism. Different tissues, health and age conditions are characterized by a difference in mitochondrial structure and composition. The visual data mining platform mitoXplorer 1.0 was developed to explore the expression dynamics of genes associated with mitochondrial functions that could help explain these differences. It, however, lacked functions aimed at integrating mitochondria in the cellular context and thus identifying regulators that help mitochondria adapt to cellular needs. To fill this gap, we upgraded the mitoXplorer platform to version 2.0 (mitoXplorer 2.0). In this upgrade, we implemented two novel integrative functions, network analysis and transcription factor enrichment, to specifically help identify signalling or transcriptional regulators of mitochondrial processes. In addition, we implemented several other novel functions to allow the platform to go beyond simple data visualization, such as an enrichment function for mitochondrial processes, a function to explore time-series data, the possibility to compare datasets across species and an IDconverter to help facilitate data upload. We demonstrate the usefulness of these functions in three specific use cases. mitoXplorer 2.0 is freely available without login at http://mitoxplorer2.ibdm.univ-mrs.fr.

Myeloid antigen-presenting cell niches sustain antitumor T cells and license PD-1 blockade via CD28 costimulation

The mechanisms regulating exhaustion of tumor-infiltrating lymphocytes (TIL) and responsiveness to PD-1 blockade remain partly unknown. In human ovarian cancer, we show that tumor-specific CD8+ TIL accumulate in tumor islets, where they engage antigen and upregulate PD-1, which restrains their functions. Intraepithelial PD-1+CD8+ TIL can be, however, polyfunctional. PD-1+ TIL indeed exhibit a continuum of exhaustion states, with variable levels of CD28 costimulation, which is provided by antigen-presenting cells (APC) in intraepithelial tumor myeloid niches. CD28 costimulation is associated with improved effector fitness of exhausted CD8+ TIL and is required for their activation upon PD-1 blockade, which also requires tumor myeloid APC. Exhausted TIL lacking proper CD28 costimulation in situ fail to respond to PD-1 blockade, and their response may be rescued by local CTLA-4 blockade and tumor APC stimulation via CD40L.

Cyclin D2 overexpression drives B1a-derived MCL-like lymphoma in mice

Mantle cell lymphoma (MCL) is an aggressive B cell lymphoma with poor long-term overall survival. Currently, MCL research and development of potential cures is hampered by the lack of good in vivo models. MCL is characterized by recurrent translocations of CCND1 or CCND2, resulting in overexpression of the cell cycle regulators cyclin D1 or D2, respectively. Here, we show, for the first time, that hematopoiesis-specific activation of cyclin D2 is sufficient to drive murine MCL-like lymphoma development. Furthermore, we demonstrate that cyclin D2 overexpression can synergize with loss of p53 to form aggressive and transplantable MCL-like lymphomas. Strikingly, cyclin D2-driven lymphomas display transcriptional, immunophenotypic, and functional similarities with B1a B cells. These MCL-like lymphomas have B1a-specific B cell receptors (BCRs), show elevated BCR and NF-κB pathway activation, and display increased MALT1 protease activity. Finally, we provide preclinical evidence that inhibition of MALT1 protease activity, which is essential for the development of early life-derived B1a cells, can be an effective therapeutic strategy to treat MCL.

Variation in Pleiotropic Hub Gene Expression Is Associated with Interspecific Differences in Head Shape and Eye Size in Drosophila

Revealing the mechanisms underlying the breathtaking morphological diversity observed in nature is a major challenge in Biology. It has been established that recurrent mutations in hotspot genes cause the repeated evolution of morphological traits, such as body pigmentation or the gain and loss of structures. To date, however, it remains elusive whether hotspot genes contribute to natural variation in the size and shape of organs. As natural variation in head morphology is pervasive in Drosophila, we studied the molecular and developmental basis of differences in compound eye size and head shape in two closely related Drosophila species. We show differences in the progression of retinal differentiation between species and we applied comparative transcriptomics and chromatin accessibility data to identify the GATA transcription factor Pannier (Pnr) as central factor associated with these differences. Although the genetic manipulation of Pnr affected multiple aspects of dorsal head development, the effect of natural variation is restricted to a subset of the phenotypic space. We present data suggesting that this developmental constraint is caused by the coevolution of expression of pnr and its cofactor u-shaped (ush). We propose that natural variation in expression or function of highly connected developmental regulators with pleiotropic functions is a major driver for morphological evolution and we discuss implications on gene regulatory network evolution. In comparison to previous findings, our data strongly suggest that evolutionary hotspots are not the only contributors to the repeated evolution of eye size and head shape in Drosophila.

Oncogenic cooperation between TCF7-SPI1 and NRAS(G12D) requires β-catenin activity to drive T-cell acute lymphoblastic leukemia

Spi-1 Proto-Oncogene (SPI1) fusion genes are recurrently found in T-cell acute lymphoblastic leukemia (T-ALL) cases but are insufficient to drive leukemogenesis. Here we show that SPI1 fusions in combination with activating NRAS mutations drive an immature T-ALL in vivo using a conditional bone marrow transplant mouse model. Addition of the oncogenic fusion to the NRAS mutation also results in a higher leukemic stem cell frequency. Mechanistically, genetic deletion of the β-catenin binding domain within Transcription factor 7 (TCF7)-SPI1 or use of a TCF/β-catenin interaction antagonist abolishes the oncogenic activity of the fusion. Targeting the TCF7-SPI1 fusion in vivo with a doxycycline-inducible knockdown results in increased differentiation. Moreover, both pharmacological and genetic inhibition lead to down-regulation of SPI1 targets. Together, our results reveal an example where TCF7-SPI1 leukemia is vulnerable to pharmacological targeting of the TCF/β-catenin interaction.

Fusion transcripts FYN-TRAF3IP2 and KHDRBS1-LCK hijack T cell receptor signaling in peripheral T-cell lymphoma, not otherwise specified

Peripheral T-cell lymphoma (PTCL) is a heterogeneous group of non-Hodgkin lymphomas with poor prognosis. Up to 30% of PTCL lack distinctive features and are classified as PTCL, not otherwise specified (PTCL-NOS). To further improve our understanding of the genetic landscape and biology of PTCL-NOS, we perform RNA-sequencing of 18 cases and validate results in an independent cohort of 37 PTCL cases. We identify FYN-TRAF3IP2, KHDRBS1-LCK and SIN3A-FOXO1 as new in-frame fusion transcripts, with FYN-TRAF3IP2 as a recurrent fusion detected in 8 of 55 cases. Using ex vivo and in vivo experiments, we demonstrate that FYN-TRAF3IP2 and KHDRBS1-LCK activate signaling pathways downstream of the T cell receptor (TCR) complex and confer therapeutic vulnerability to clinically available drugs.