JASPAR 2020: update of the open-access database of transcription factor binding profiles

Krüppel-like factor 4 in transcriptional control of the three unique isoforms of Agouti-related peptide in mice

Agouti-related peptide (AgRP/Agrp) within the hypothalamic arcuate nucleus (ARC) contributes to the control of energy balance, and dysregulated Agrp may contribute to metabolic adaptation during prolonged obesity. In mice, three isoforms of Agrp are encoded via distinct first exons. Agrp-A (ENSMUST00000005849.11) contributed 95% of total Agrp in mouse ARC, whereas Agrp-B (ENSMUST00000194654.2) dominated in placenta (73%). Conditional deletion of Klf4 from Agrp-expressing cells (Klf4Agrp-KO mice) reduced Agrp mRNA and increased energy expenditure but had no effects on food intake or the relative abundance of Agrp isoforms in the ARC. Chronic high-fat diet feeding masked these effects of Klf4 deletion, highlighting the context-dependent contribution of KLF4 to Agrp control. In the GT1-7 mouse hypothalamic cell culture model, which expresses all three isoforms of Agrp (including Agrp-C, ENSMUST00000194091.6), inhibition of extracellular signal-regulated kinase (ERK) simultaneously increased KLF4 binding to the Agrp promoter and stimulated Agrp expression. In addition, siRNA-mediated knockdown of Klf4 reduced expression of Agrp. We conclude that the expression of individual isoforms of Agrp in the mouse is dependent upon cell type and that KLF4 directly promotes the transcription of Agrp via a mechanism that is superseded during obesity.NEW & NOTEWORTHY In mice, three distinct isoforms of Agouti-related peptide are encoded via distinct first exons. In the arcuate nucleus of the hypothalamus, Krüppel-like factor 4 stimulates transcription of the dominant isoform in lean mice, but this mechanism is altered during diet-induced obesity.

Id2 epigenetically controls CD8+ T-cell exhaustion by disrupting the assembly of the Tcf3-LSD1 complex

CD8+ T-cell exhaustion is a state of dysfunction that promotes tumor progression and is marked by the generation of Slamf6+ progenitor exhausted (Texprog) and Tim-3+ terminally exhausted (Texterm) subpopulations. Inhibitor of DNA binding protein 2 (Id2) has been shown to play important roles in T-cell development and CD8+ T-cell immunity. However, the role of Id2 in CD8+ T-cell exhaustion is unclear. Here, we found that Id2 transcriptionally and epigenetically regulates the generation of Texprog cells and their conversion to Texterm cells. Genetic deletion of Id2 dampens CD8+ T-cell-mediated immune responses and the maintenance of stem-like CD8+ T-cell subpopulations, suppresses PD-1 blockade and increases tumor susceptibility. Mechanistically, through its HLH domain, Id2 binds and disrupts the assembly of the Tcf3-Tal1 transcriptional regulatory complex, and thus modulates chromatin accessibility at the Slamf6 promoter by preventing the interaction of Tcf3 with the histone lysine demethylase LSD1. Therefore, Id2 increases the abundance of the permissive H3K4me2 mark on the Tcf3-occupied E-boxes in the Slamf6 promoter, modulates chromatin accessibility at the Slamf6 promoter and epigenetically regulates the generation of Slamf6+ Texprog cells. An LSD1 inhibitor GSK2879552 can rescue the Id2 knockout phenotype in tumor-bearing mice. Inhibition of LSD1 increases the abundance of Slamf6+Tim-3- Texprog cells in tumors and the expression level of Tcf1 in Id2-deleted CD8+ T cells. This study demonstrates that Id2-mediated transcriptional and epigenetic modification drives hierarchical CD8+ T-cell exhaustion, and the mechanistic insights gained may have implications for therapeutic intervention with tumor immune evasion.

Adaptive Selection of Cis-regulatory Elements in the Han Chinese

Cis-regulatory elements have an important role in human adaptation to the living environment. However, the lag in population genomic cohort studies and epigenomic studies, hinders the research in the adaptive analysis of cis-regulatory elements in human populations. In this study, we collected 4,013 unrelated individuals and performed a comprehensive analysis of adaptive selection of genome-wide cis-regulatory elements in the Han Chinese. In total, 12.34% of genomic regions are under the influence of adaptive selection, where 1.00% of enhancers and 2.06% of promoters are under positive selection, and 0.06% of enhancers and 0.02% of promoters are under balancing selection. Gene ontology enrichment analysis of these cis-regulatory elements under adaptive selection reveals that many positive selections in the Han Chinese occur in pathways involved in cell-cell adhesion processes, and many balancing selections are related to immune processes. Two classes of adaptive cis-regulatory elements related to cell adhesion were in-depth analyzed, one is the adaptive enhancers derived from neanderthal introgression, leads to lower hyaluronidase level in skin, and brings better performance on UV-radiation resistance to the Han Chinese. Another one is the cis-regulatory elements regulating wound healing, and the results suggest the positive selection inhibits coagulation and promotes angiogenesis and wound healing in the Han Chinese. Finally, we found that many pathogenic alleles, such as risky alleles of type 2 diabetes or schizophrenia, remain in the population due to the hitchhiking effect of positive selections. Our findings will help deepen our understanding of the adaptive evolution of genome regulation in the Han Chinese.

Sall4 regulates posterior trunk mesoderm development by promoting mesodermal gene expression and repressing neural genes in the mesoderm

The trunk axial skeleton develops from paraxial mesoderm cells. Our recent study demonstrated that conditional knockout of the stem cell factor Sall4 in mice by TCre caused tail truncation and a disorganized axial skeleton posterior to the lumbar level. Based on this phenotype, we hypothesized that, in addition to the previously reported role of Sall4 in neuromesodermal progenitors, Sall4 is involved in the development of the paraxial mesoderm tissue. Analysis of gene expression and SALL4 binding suggests that Sall4 directly or indirectly regulates genes involved in presomitic mesoderm differentiation, somite formation and somite differentiation. Furthermore, ATAC-seq in TCre; Sall4 mutant posterior trunk mesoderm shows that Sall4 knockout reduces chromatin accessibility. We found that Sall4-dependent open chromatin status drives activation and repression of WNT signaling activators and repressors, respectively, to promote WNT signaling. Moreover, footprinting analysis of ATAC-seq data suggests that Sall4-dependent chromatin accessibility facilitates CTCF binding, which contributes to the repression of neural genes within the mesoderm. This study unveils multiple mechanisms by which Sall4 regulates paraxial mesoderm development by directing activation of mesodermal genes and repression of neural genes.

Essential role of PLD2 in hypoxia-induced stemness and therapy resistance in ovarian tumors

BACKGROUND

Hypoxia in solid tumors is an important source of chemoresistance that can determine poor patient prognosis. Such chemoresistance relies on the presence of cancer stem cells (CSCs), and hypoxia promotes their generation through transcriptional activation by HIF transcription factors.

METHODS

We used ovarian cancer (OC) cell lines, xenograft models, OC patient samples, transcriptional databases, induced pluripotent stem cells (iPSCs) and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq).

RESULTS

Here, we show that hypoxia induces CSC formation and chemoresistance in ovarian cancer through transcriptional activation of the PLD2 gene. Mechanistically, HIF-1α activates PLD2 transcription through hypoxia response elements, and both hypoxia and PLD2 overexpression lead to increased accessibility around stemness genes, detected by ATAC-seq, at sites bound by AP-1 transcription factors. This in turn provokes a rewiring of stemness genes, including the overexpression of SOX2, SOX9 or NOTCH1. PLD2 overexpression also leads to decreased patient survival, enhanced tumor growth and CSC formation, and increased iPSCs reprograming, confirming its role in dedifferentiation to a stem-like phenotype. Importantly, hypoxia-induced stemness is dependent on PLD2 expression, demonstrating that PLD2 is a major determinant of de-differentiation of ovarian cancer cells to stem-like cells in hypoxic conditions. Finally, we demonstrate that high PLD2 expression increases chemoresistance to cisplatin and carboplatin treatments, both in vitro and in vivo, while its pharmacological inhibition restores sensitivity.

CONCLUSIONS

Altogether, our work highlights the importance of the HIF-1α-PLD2 axis for CSC generation and chemoresistance in OC and proposes an alternative treatment for patients with high PLD2 expression.

Comprehensive bioinformatics analysis and systems biology approaches to identify the interplay between COVID-19 and pericarditis

Background

Increasing evidence indicating that coronavirus disease 2019 (COVID-19) increased the incidence and related risks of pericarditis and whether COVID-19 vaccine is related to pericarditis has triggered research and discussion. However, mechanisms behind the link between COVID-19 and pericarditis are still unknown. The objective of this study was to further elucidate the molecular mechanisms of COVID-19 with pericarditis at the gene level using bioinformatics analysis.

Methods

Genes associated with COVID-19 and pericarditis were collected from databases using limited screening criteria and intersected to identify the common genes of COVID-19 and pericarditis. Subsequently, gene ontology, pathway enrichment, protein-protein interaction, and immune infiltration analyses were conducted. Finally, TF-gene, gene-miRNA, gene-disease, protein-chemical, and protein-drug interaction networks were constructed based on hub gene identification.

Results

A total of 313 common genes were selected, and enrichment analyses were performed to determine their biological functions and signaling pathways. Eight hub genes (IL-1β, CD8A, IL-10, CD4, IL-6, TLR4, CCL2, and PTPRC) were identified using the protein-protein interaction network, and immune infiltration analysis was then carried out to examine the functional relationship between the eight hub genes and immune cells as well as changes in immune cells in disease. Transcription factors, miRNAs, diseases, chemicals, and drugs with high correlation with hub genes were predicted using bioinformatics analysis.

Conclusions

This study revealed a common gene interaction network between COVID-19 and pericarditis. The screened functional pathways, hub genes, potential compounds, and drugs provided new insights for further research on COVID-19 associated with pericarditis.

CHOP-mediated Gasdermin E expression promotes pyroptosis, inflammation, and mitochondrial damage in renal ischemia-reperfusion injury

In clinical practice, renal ischemia-reperfusion injury (IRI) is a common cause of acute kidney injury (AKI), often leading to acute renal failure or end-stage renal disease (ESRD). The current understanding of renal IRI mechanisms remains unclear, and effective therapeutic strategies and clear targets are lacking. Therefore, the need to find explicit and effective ways to reduce renal IRI remains a scientific challenge. The current study explored pyroptosis, a type of inflammation-regulated programmed cell death, and the role of Gasdermins E (GSDME)-mediated pyroptosis, mitochondrial damage, and inflammation in renal IRI. The analysis of human samples showed that the expression levels of GSDME in normal human renal tissues were higher than those of GSDMD. Moreover, our study demonstrated that GSDME played an important role in mediating pyroptosis, inflammation, and mitochondrial damage in renal IRI. Subsequently, GSDME-N accumulated in the mitochondrial membrane, leading to mitochondrial damage and activation of caspase3, which generated a feed-forward loop of self-amplification injury. However, GSDME knockout resulted in the amelioration of renal IRI. Moreover, the current study found that the transcription factor CHOP was activated much earlier in renal IRI. Inhibition of BCL-2 by CHOP leaded to casapse3 activation, resulting in mitochondrial damage and apoptosis; not only that, but CHOP positively regulated GSDME thereby causing pyroptosis. Therefore, this study explored the transcriptional mechanisms of GSDME during IRI development and the important role of CHOP/Caspase3/GSDME mechanistic axis in regulating pyroptosis in renal IRI. This axis might serve as a potential therapeutic target.

Comprehensive transcriptomic analysis revealing the regulatory dynamics and networks of the pituitary-testis axis in sheep across developmental stages

Spermatogenesis is a complex process intricately regulated by the hypothalamic-pituitary-testis (HPT) axis. However, research on the regulatory factors governing the HPT axis remains limited. This study addresses this gap by conducting a comprehensive analysis of transcriptomes from the pituitary and testis tissues across various developmental stages, encompassing embryonic day (E120), neonatal period (P0), pre-puberty (P90), and post-puberty day (P270). Utilizing edgeR and WGCNA, we identified stage-specific genes in both the pituitary and testis throughout the four developmental stages. Notably, 380, 242, 34, and 479 stage-specific genes were identified in the pituitary, while 886, 297, 201, and 3,678 genes were identified in the testis. Subsequent analyses unveiled associations between these stage-specific genes and crucial pathways such as the cAMP signaling pathway, GnRH secretion, and male gamete generation. Furthermore, leveraging single-cell data from the pituitary and testis, we identified some signaling pathways involving BMP, HGF, IGF, and TGF-β, highlighting mutual regulation between the pituitary and testis at different developmental stages. This study sheds light on the pivotal role of the pituitary-testis axis in the reproductive process of sheep across four distinct developmental stages. Additionally, it delves into the intricate regulatory networks governing reproduction, offering novel insights into the dynamics of the pituitary-testis axis within the reproductive system.

Endothelial downregulation of nuclear m6A reader YTHDC1 promotes pulmonary vascular remodeling in sugen hypoxia model of pulmonary hypertension

Background

Pulmonary hypertension (PH) is characterized with vascular remodeling, which is intiated by vascular endothelial dysfunction. N6-methyladenosine (m6A) modification mediates gene expression in many ways including mediating RNA degradation, splicing, nuclear export et al. m6A modification have been found to be associated with the development of PH. However, the role of m6A regulators in pulmonary artery endothelial cells (PAECs) dysfunction of PH is still under research.

Methods

The expression levels of m6A regulators in PAECs were analyzed with the single-cell sequencing Data(scRNA). Next, the target differentially expressed genes (DEGs) of m6A regulators in PAECs were functionally annotated. The analysis of cellular interactions included the examination of receptor-ligand pairs regulated by m6A regulators. Pseudo-time trajectory analyses and a ceRNA network involving lncRNAs, miRNAs, and mRNAs were conducted in PAECs. Furthermore, microarray data (GSE180169) for Sugen Hypoxia PH (SuHx PH) mouse models was screened for DEGs and m6A regulators in PAECs. Moreover, the expression of YTHDC1 in the lung samples of SuHx PH models was determined using immunofluorescence. In vitro, the mRNA expression of YTHDC1 in HPAECs under hypoxia conditions was detected. The effect of YTHDC1 recombinant protein on HPAEC proliferation was detected by Cell Counting Kit-8 (CCK8).

Results

Dysregulation of m6A regulators was observed in mouse PAECs. The m6A reader of YTHDC1 was decreased in PAECs in scRNA data and RNAseq data of isolated PAECs of SuHx PH models. Downregulation of YTHDC1 was caused by hypoxia in PAECs in vitro and similar results was observed in PAECs of SuHx PH mouse models. Next, YTHDC1 recombinant protein was found to inhibit HPAECs proliferation. The DEGs targeted by YTHDC1 were enriched in angiogenesis, endothelial cell migration, fluid shear stress, and stem cell maintenance. Analysis indicates that interactions among endothelial cells, smooth muscle cells, fibroblasts, and immune cells, mediated by specific YTHDC1 target genes (e.g., PTPRC-MRC1, ITBG2-ICAM1, COL4A1-CD44), contribute to PH development. Also, the YTHDC1 expression were consistent with Thioredoxin interacting protein (TXNIP). What's more, the predicted transcription factors showed that NFKB1, Foxd3 may be involved in the regulation of YTHDC1. Lastly, our data suggest that YTHDC1 may be involved in regulating PAECs dysfunction through lncRNA/miRNA/mRNA network.

Conclusion

For the first time, we analyzed changes in the expression and biological functions of m6A regulators in SuHx PH mouse models. We causatively linked YTHDC1 to PAECs dysfunction, providing novel insight into and opportunities to diagnose and treat PH.

Transcriptome analysis in Aegilops tauschii unravels further insights into genetic control of stripe rust resistance

MAIN CONCLUSION

The Aegilops tauschii resistant accession prevented the pathogen colonization by controlling the sugar flow and triggering the hypersensitive reaction. This study suggested that NBS-LRRs probably induce resistance through bHLH by controlling JA- and SA-dependent pathways. Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst) is one of wheat's most destructive fungal diseases that causes a severe yield reduction worldwide. The most effective and economically-friendly strategy to manage this disease is genetic resistance which can be achieved through deploying new and effective resistance genes. Aegilops tauschii, due to its small genome and co-evolution with Pst, can provide detailed information about underlying resistance mechanisms. Hence, we used RNA-sequencing approach to identify the transcriptome variations of two contrasting resistant and susceptible Ae. tauschii accessions in interaction with Pst and differentially expressed genes (DEGs) for resistance to stripe rust. Gene ontology, pathway analysis, and search for functional domains, transcription regulators, resistance genes, and protein-protein interactions were used to interpret the results. The genes encoding NBS-LRR, CC-NBS-kinase, and phenylalanine ammonia-lyase, basic helix-loop-helix (bHLH)-, basic-leucine zipper (bZIP)-, APETALA2 (AP2)-, auxin response factor (ARF)-, GATA-, and LSD-like transcription factors were up-regulated exclusively in the resistant accession. The key genes involved in response to salicylic acid, amino sugar and nucleotide sugar metabolism, and hypersensitive response contributed to plant defense against stripe rust. The activation of jasmonic acid biosynthesis and starch and sucrose metabolism pathways under Pst infection in the susceptible accession explained the colonization of the host. Overall, this study can fill the gaps in the literature on host-pathogen interaction and enrich the Ae. tauschii transcriptome sequence information. It also suggests candidate genes that could guide future breeding programs attempting to develop rust-resistant cultivars.

A low-input high resolution sequential chromatin immunoprecipitation method captures genome-wide dynamics of bivalent chromatin

BACKGROUND

Bivalent chromatin is an exemplar of epigenetic plasticity. This co-occurrence of active-associated H3K4me3 and inactive-associated H3K27me3 histone modifications on opposite tails of the same nucleosome occurs predominantly at promoters that are poised for future transcriptional upregulation or terminal silencing. We know little of the dynamics, resolution, and regulation of this chromatin state outside of embryonic stem cells where it was first described. This is partly due to the technical challenges distinguishing bone-fide bivalent chromatin, where both marks are on the same nucleosome, from allelic or sample heterogeneity where there is a mix of H3K4me3-only and H3K27me3-only mononucleosomes.

RESULTS

Here, we present a robust and sensitive method to accurately map bivalent chromatin genome-wide, along with controls, from as little as 2 million cells. We optimized and refined the sequential ChIP protocol which uses two sequential overnight immunoprecipitation reactions to robustly purify nucleosomes that are truly bivalent and contain both H3K4me3 and H3K27me3 modifications. Our method generates high quality genome-wide maps with strong peak enrichment and low background, which can be analyzed using standard bioinformatic packages. Using this method, we detect 8,789 bivalent regions in mouse embryonic stem cells corresponding to 3,918 predominantly CpG rich and developmentally regulated gene promoters. Furthermore, profiling Dppa2/4 knockout mouse embryonic stem cells, which lose both H3K4me3 and H3K27me3 at approximately 10% of bivalent promoters, demonstrated the ability of our method to capture bivalent chromatin dynamics.

CONCLUSIONS

Our optimized sequential reChIP method enables high-resolution genome-wide assessment of bivalent chromatin together with all required controls in as little as 2 million cells. We share a detailed protocol and guidelines that will enable bivalent chromatin landscapes to be generated in a range of cellular contexts, greatly enhancing our understanding of bivalent chromatin and epigenetic plasticity beyond embryonic stem cells.

Cholesterol lowering depletes atherosclerotic lesions of smooth muscle cell-derived fibromyocytes and chondromyocytes

Drugs that lower plasma apolipoprotein B (ApoB)-containing lipoproteins are central to treating advanced atherosclerosis and provide partial protection against clinical events. Previous research showed that lowering ApoB-containing lipoproteins stops plaque inflammation, but how these drugs affect the heterogeneous population of plaque cells derived from smooth muscle cells (SMCs) is unknown. SMC-derived cells are the main cellular component of atherosclerotic lesions and the source of structural components that determine the size of plaques and their propensity to rupture and trigger thrombosis, the proximate cause of heart attack and stroke. Using lineage tracing and single-cell techniques to investigate the full SMC-derived cellular compartment in progressing and regressing plaques in mice, here we show that lowering ApoB-containing lipoproteins reduces nuclear factor kappa-light-chain-enhancer of activated B cells signaling in SMC-derived fibromyocytes and chondromyocytes and leads to depletion of these abundant cell types from plaques. These results uncover an important mechanism through which cholesterol-lowering drugs can achieve plaque regression.

Long noncoding RNA lncPostn links TGF-β and p53 signaling pathways to transcriptional regulation of cardiac fibrosis

Cardiac fibroblasts are essential for the homeostasis of the extracellular matrix, whose remodeling in many cardiovascular diseases leads to fibrosis. Long noncoding RNAs (lncRNAs) are associated with cardiac pathologies, but their functions in cardiac fibroblasts and contributions to cardiac fibrosis remain unclear. Here, we aimed to identify fibroblast-enriched lncRNAs essential in myocardial infarction (MI)-induced fibrosis and explore the molecular mechanisms responsible for their functions. Global lncRNA profiling was performed in post-MI mouse heart ventricles and transforming growth factor-β (TGF-β)-treated primary cardiac fibroblasts and confirmed in published data sets. We identified the cardiac fibroblast-enriched lncPostn, whose expression is stimulated in cardiac fibrosis induced by MI and the extracellular growth factor TGF-β. The promoter of lncPostn contains a functional TGF-β response element, and lncPostn knockdown suppresses TGF-β-stimulated cardiac fibroblast activation and improves cardiac functions post-MI. LncPostn stabilizes and recruits EP300 to the profibrotic periostin's promoter, representing a major mechanism for its transcriptional activation. Moreover, both MI and TGF-β enhance lncPostn expression while suppressing the cellular growth gatekeeper p53. TGF-β and p53 knockdown-induced profibrotic gene expression and fibrosis occur mainly through lncPostn and show additive effects. Finally, levels of serum lncPostn are significantly increased in patients' postacute MI and show a strong correlation with fibrosis markers, revealing a potential biomarker of cardiac fibrosis. Our findings identify the fibroblast-enriched lncPostn as a potent profibrotic factor, providing a transcriptional link between TGF-β and p53 signaling pathways to regulate fibrosis in cardiac fibroblasts.NEW & NOTEWORTHY Cardiac fibroblasts are essential for the homeostasis of the extracellular matrix, whose remodeling in many cardiovascular diseases leads to fibrosis. Long noncoding RNAs are functional and contribute to the biological processes of cardiovascular development and disorders. Our findings identify the fibroblast-enriched lncPostn as a potent profibrotic factor and demonstrate that serum lncPostn level may serve as a potential biomarker of human cardiac fibrosis postacute myocardial infarction.

KLF9‑regulated FBXO31 inhibits the progression of endometrial cancer and enhances the sensitivity of endometrial cancer cells to cisplatin

Endometrial cancer (EC) is one of the most common malignancies with an increasing annual incidence. F-box only protein 31 (FBXO31) plays a significant regulatory role in several types of cancer. The transcription factor Krüppel-like factor 9 (KLF9) of FBXO31 is reduced in EC as a tumor suppressor. However, their particular regulatory role and mechanism in EC have not been previously reported. Therefore, the UALCAN database was used to predict the expression levels of FBXO31 in EC. In addition, the regulatory effect of FBXO31 on EC cell proliferation, invasion, migration, apoptosis and cisplatin (DDP) sensitivity was investigated at the cellular level. The association between KLF9 and FBXO31 was predicted using the JASPAR database and verified by chromatin immunoprecipitation and luciferase reporter assays. Finally, the regulatory effects of KLF9 and FBXO31 overexpression or silencing were also explored. The results demonstrated that FBXO31 was poorly expressed in EC. Additionally, FBXO31 overexpression inhibited the malignant progression of EC cells and enhanced their sensitivity to DDP. Furthermore, KLF9 promoted FBXO31 transcription. Overall, the present study suggested that the KLF9-mediated regulation of FBXO31 could inhibit the progression of EC and enhance the sensitivity of EC cells to DDP.

Discovery of common molecular signatures and drug repurposing for COVID-19/Asthma comorbidity: ACE2 and multi-partite networks

Angiotensin-converting enzyme 2 (ACE2) is identified as the functional receptor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the ongoing global coronavirus disease-2019 (COVID-19) pandemic. This study aimed to elucidate potential therapeutic avenues by scrutinizing approved drugs through the identification of the genetic signature associated with SARS-CoV-2 infection in individuals with asthma. This exploration was conducted through an integrated analysis, encompassing interaction networks between the ACE2 receptor and common host (co-host) factors implicated in COVID-19/asthma comorbidity. The comprehensive analysis involved the identification of common differentially expressed genes (cDEGs) and hub-cDEGs, functional annotations, interaction networks, gene set variation analysis (GSVA), gene set enrichment analysis (GSEA), and module construction. Interaction networks were used to identify overlapping disease modules and potential drug targets. Computational biology and molecular docking analyzes were utilized to discern functional drug modules. Subsequently, the impact of the identified drugs on the expression of hub-cDEGs was experimentally validated using a mouse model. A total of 153 cDEGs or co-host factors associated with ACE2 were identified in the COVID-19 and asthma comorbidity. Among these, seven significant cDEGs and proteins - namely, HRAS, IFNG, JUN, CDH1, TLR4, ICAM1, and SCD-were recognized as pivotal host factors linked to ACE2. Regulatory network analysis of hub-cDEGs revealed eight top-ranked transcription factors (TFs) proteins and nine microRNAs as key regulatory factors operating at the transcriptional and post-transcriptional levels, respectively. Molecular docking simulations led to the proposal of 10 top-ranked repurposable drug molecules (Rapamycin, Ivermectin, Everolimus, Quercetin, Estradiol, Entrectinib, Nilotinib, Conivaptan, Radotinib, and Venetoclax) as potential treatment options for COVID-19 in individuals with comorbid asthma. Validation analysis demonstrated that Rapamycin effectively inhibited ICAM1 expression in the HDM-stimulated mice group (p < 0.01). This study unveils the common pathogenesis and genetic signature underlying asthma and SARS-CoV-2 infection, delineated by the interaction networks of ACE2-related host factors. These findings provide valuable insights for the design and discovery of drugs aimed at more effective therapeutics within the context of lung disease comorbidities.

A type 1 immunity-restricted promoter of the IL-33 receptor gene directs antiviral T-cell responses

The pleiotropic alarmin interleukin-33 (IL-33) drives type 1, type 2 and regulatory T-cell responses via its receptor ST2. Subset-specific differences in ST2 expression intensity and dynamics suggest that transcriptional regulation is key in orchestrating the context-dependent activity of IL-33-ST2 signaling in T-cell immunity. Here, we identify a previously unrecognized alternative promoter in mice and humans that is located far upstream of the curated ST2-coding gene and drives ST2 expression in type 1 immunity. Mice lacking this promoter exhibit a selective loss of ST2 expression in type 1- but not type 2-biased T cells, resulting in impaired expansion of cytotoxic T cells (CTLs) and T-helper 1 cells upon viral infection. T-cell-intrinsic IL-33 signaling via type 1 promoter-driven ST2 is critical to generate a clonally diverse population of antiviral short-lived effector CTLs. Thus, lineage-specific alternative promoter usage directs alarmin responsiveness in T-cell subsets and offers opportunities for immune cell-specific targeting of the IL-33-ST2 axis in infections and inflammatory diseases.

Nucleus pulposus cells regulate macrophages in degenerated intervertebral discs via the integrated stress response-mediated CCL2/7-CCR2 signaling pathway

Lower back pain (LBP), which is a primary cause of disability, is largely attributed to intervertebral disc degeneration (IDD). Macrophages (MΦs) in degenerated intervertebral discs (IVDs) form a chronic inflammatory microenvironment, but how MΦs are recruited to degenerative segments and transform into a proinflammatory phenotype remains unclear. We evaluated chemokine expression in degenerated nucleus pulposus cells (NPCs) to clarify the role of NPCs in the establishment of an inflammatory microenvironment in IDD and explored the mechanisms. We found that the production of C-C motif chemokine ligand 2 (CCL2) and C-C motif chemokine ligand 7 (CCL7) was significantly increased in NPCs under inflammatory conditions, and blocking CCL2/7 and their receptor, C-C chemokine receptor type 2(CCR2), inhibited the inductive effects of NPCs on MΦ infiltration and proinflammatory polarization. Moreover, activation of the integrated stress response (ISR) was obvious in IDD, and ISR inhibition reduced the production of CCL2/7 in NPCs. Further investigation revealed that activating Transcription Factor 3 (ATF3) responded to ISR activation, and ChIP-qPCR verified the DNA-binding activity of ATF3 on CCL2/7 promoters. In addition, we found that Toll-like receptor 4 (TLR4) inhibition modulated ISR activation, and TLR4 regulated the accumulation of mitochondrial reactive oxygen species (mtROS) and double-stranded RNA (dsRNA). Downregulating the level of mtROS reduced the amount of dsRNA and ISR activation. Deactivating the ISR or blocking CCL2/7 release alleviated inflammation and the progression of IDD in vivo. Moreover, MΦ infiltration and IDD were inhibited in CCR2-knockout mice. In conclusion, this study highlights the critical role of TLR4/mtROS/dsRNA axis-mediated ISR activation in the production of CCL2/7 and the progression of IDD, which provides promising therapeutic strategies for discogenic LBP.

The nuclear factor ID3 endows macrophages with a potent anti-tumour activity

Macrophage activation is controlled by a balance between activating and inhibitory receptors1-7, which protect normal tissues from excessive damage during infection8,9 but promote tumour growth and metastasis in cancer7,10. Here we report that the Kupffer cell lineage-determining factor ID3 controls this balance and selectively endows Kupffer cells with the ability to phagocytose live tumour cells and orchestrate the recruitment, proliferation and activation of natural killer and CD8 T lymphoid effector cells in the liver to restrict the growth of a variety of tumours. ID3 shifts the macrophage inhibitory/activating receptor balance to promote the phagocytic and lymphoid response, at least in part by buffering the binding of the transcription factors ELK1 and E2A at the SIRPA locus. Furthermore, loss- and gain-of-function experiments demonstrate that ID3 is sufficient to confer this potent anti-tumour activity to mouse bone-marrow-derived macrophages and human induced pluripotent stem-cell-derived macrophages. Expression of ID3 is therefore necessary and sufficient to endow macrophages with the ability to form an efficient anti-tumour niche, which could be harnessed for cell therapy in cancer.

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.

Cell layer-specific expression of the homeotic MADS-box transcription factor PhDEF contributes to modular petal morphogenesis in petunia

Floral homeotic MADS-box transcription factors ensure the correct morphogenesis of floral organs, which are organized in different cell layers deriving from distinct meristematic layers. How cells from these distinct layers acquire their respective identities and coordinate their growth to ensure normal floral organ morphogenesis is unresolved. Here, we studied petunia (Petunia × hybrida) petals that form a limb and tube through congenital fusion. We identified petunia mutants (periclinal chimeras) expressing the B-class MADS-box gene DEFICIENS in the petal epidermis or in the petal mesophyll, called wico and star, respectively. Strikingly, wico flowers form a strongly reduced tube while their limbs are almost normal, while star flowers form a normal tube but greatly reduced and unpigmented limbs, showing that petunia petal morphogenesis is highly modular. These mutants highlight the layer-specific roles of PhDEF during petal development. We explored the link between PhDEF and petal pigmentation, a well-characterized limb epidermal trait. The anthocyanin biosynthesis pathway was strongly downregulated in star petals, including its major regulator ANTHOCYANIN2 (AN2). We established that PhDEF directly binds to the AN2 terminator in vitro and in vivo, suggesting that PhDEF might regulate AN2 expression and therefore petal epidermis pigmentation. Altogether, we show that cell layer-specific homeotic activity in petunia petals differently impacts tube and limb development, revealing the relative importance of the different cell layers in the modular architecture of petunia petals.

Single-cell multiomics decodes regulatory programs for mouse secondary palate development

Perturbations in gene regulation during palatogenesis can lead to cleft palate, which is among the most common congenital birth defects. Here, we perform single-cell multiome sequencing and profile chromatin accessibility and gene expression simultaneously within the same cells (n = 36,154) isolated from mouse secondary palate across embryonic days (E) 12.5, E13.5, E14.0, and E14.5. We construct five trajectories representing continuous differentiation of cranial neural crest-derived multipotent cells into distinct lineages. By linking open chromatin signals to gene expression changes, we characterize the underlying lineage-determining transcription factors. In silico perturbation analysis identifies transcription factors SHOX2 and MEOX2 as important regulators of the development of the anterior and posterior palate, respectively. In conclusion, our study charts epigenetic and transcriptional dynamics in palatogenesis, serving as a valuable resource for further cleft palate research.

Identification of Senescence-Associated Biomarkers in Diabetic Glomerulopathy Using Integrated Bioinformatics Analysis

Background

Cellular senescence is thought to play a significant role in the onset and development of diabetic nephropathy. The goal of this study was to explore potential biomarkers associated with diabetic glomerulopathy from the perspective of senescence.

Methods

Datasets about human glomerular biopsy samples related to diabetic nephropathy were systematically obtained from the Gene Expression Omnibus database. Hub senescence-associated genes were investigated by differential gene analysis and Least Absolute Shrinkage and Selection Operator analysis. Cluster analysis was employed to identify senescence molecular subtypes. A single-cell dataset was used to validate the above findings and further evaluate the senescence environment. The relationship between these genes and the glomerular filtration rate was explored based on the Nephroseq database. These gene expressions have also been explored in various kidney diseases.

Results

Twelve representative senescence-associated genes (VEGFA, IQGAP2, JUN, PLAT, ETS2, ANG, MMP14, VEGFC, SERPINE2, CXCR2, PTGES, and EGF) were finally identified. Biological changes in immune inflammatory response, cell cycle regulation, metabolic regulation, and immune microenvironment have been observed across different molecular subtypes. The above results were also validated based on single-cell analysis. Additionally, we also identified several significantly altered cell communication pathways, including COLLAGEN, PTN, LAMININ, SPP1, and VEGF. Finally, almost all these genes could well predict the occurrence of diabetic glomerulopathy based on receiver operating characteristic analysis and are associated with the glomerular filtration rate. These genes are differently expressed in various kidney diseases.

Conclusion

The present study identified potential senescence-associated biomarkers and further explored the heterogeneity of diabetic glomerulopathy that might provide new insights into the diagnosis, assessment, management, and personalized treatment of DN.

The chromatin landscape of healthy and injured cell types in the human kidney

There is a need to define regions of gene activation or repression that control human kidney cells in states of health, injury, and repair to understand the molecular pathogenesis of kidney disease and design therapeutic strategies. Comprehensive integration of gene expression with epigenetic features that define regulatory elements remains a significant challenge. We measure dual single nucleus RNA expression and chromatin accessibility, DNA methylation, and H3K27ac, H3K4me1, H3K4me3, and H3K27me3 histone modifications to decipher the chromatin landscape and gene regulation of the kidney in reference and adaptive injury states. We establish a spatially-anchored epigenomic atlas to define the kidney's active, silent, and regulatory accessible chromatin regions across the genome. Using this atlas, we note distinct control of adaptive injury in different epithelial cell types. A proximal tubule cell transcription factor network of ELF3, KLF6, and KLF10 regulates the transition between health and injury, while in thick ascending limb cells this transition is regulated by NR2F1. Further, combined perturbation of ELF3, KLF6, and KLF10 distinguishes two adaptive proximal tubular cell subtypes, one of which manifested a repair trajectory after knockout. This atlas will serve as a foundation to facilitate targeted cell-specific therapeutics by reprogramming gene regulatory networks.

Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet

BACKGROUND

Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult.

RESULTS

Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPβ for methylated motifs and the preference of c-Myc for unmethylated E-box motifs.

CONCLUSIONS

Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

The Goldfish Genome and Its Utility for Understanding Gene Regulation and Vertebrate Body Morphology

Goldfish, widely viewed as an ornamental fish, is a member of Cyprinidae family and has a very long history in research for both genetics and physiology studies. Among Cyprinidae, the chromosomal locations of orthologs and the amino acid sequences are usually highly conserved. Adult goldfish are 1000 times larger than adult zebrafish (who are in the same family of fishes), which can make it easier to perform several types of experiments compared to their zebrafish cousins. Comparing mutant phenotypes in orthologous genes between goldfish and zebrafish can often be very informative and provide a deeper insight into the gene function than studying the gene in either species alone. Comparative genomics and phenotypic comparisons between goldfish and zebrafish will provide new opportunities for understanding the development and evolution of body forms in the vertebrate lineage.

Longitudinal modeling of human neuronal aging reveals the contribution of the RCAN1-TFEB pathway to Huntington's disease neurodegeneration

Aging is a common risk factor in neurodegenerative disorders. Investigating neuronal aging in an isogenic background stands to facilitate analysis of the interplay between neuronal aging and neurodegeneration. Here we perform direct neuronal reprogramming of longitudinally collected human fibroblasts to reveal genetic pathways altered at different ages. Comparative transcriptome analysis of longitudinally aged striatal medium spiny neurons (MSNs) in Huntington's disease identified pathways involving RCAN1, a negative regulator of calcineurin. Notably, RCAN1 protein increased with age in reprogrammed MSNs as well as in human postmortem striatum and RCAN1 knockdown rescued patient-derived MSNs of Huntington's disease from degeneration. RCAN1 knockdown enhanced chromatin accessibility of genes involved in longevity and autophagy, mediated through enhanced calcineurin activity, leading to TFEB's nuclear localization by dephosphorylation. Furthermore, G2-115, an analog of glibenclamide with autophagy-enhancing activities, reduced the RCAN1-calcineurin interaction, phenocopying the effect of RCAN1 knockdown. Our results demonstrate that targeting RCAN1 genetically or pharmacologically can increase neuronal resilience in Huntington's disease.

KLF4 facilitates chromatin accessibility remodeling in porcine early embryos

Chromatin accessibility remodeling driven by pioneer factors is critical for the development of early embryos. Current studies have illustrated several pioneer factors as being important for agricultural animals, but what are the pioneer factors and how the pioneer factors remodel the chromatin accessibility in porcine early embryos is not clear. By employing low-input DNase-seq (liDNase-seq), we profiled the landscapes of chromatin accessibility in porcine early embryos and uncovered a unique chromatin accessibility reprogramming pattern during porcine preimplantation development. Our data revealed that KLF4 played critical roles in remodeling chromatin accessibility in porcine early embryos. Knocking down of KLF4 led to the reduction of chromatin accessibility in early embryos, whereas KLF4 overexpression promoted the chromatin openness in porcine blastocysts. Furthermore, KLF4 deficiency resulted in mitochondrial dysfunction and developmental failure of porcine embryos. In addition, we found that overexpression of KLF4 in blastocysts promoted lipid droplet accumulation, whereas knockdown of KLF4 disrupted this process. Taken together, our study revealed the chromatin accessibility dynamics and identified KLF4 as a key regulator in chromatin accessibility and cellular metabolism during porcine preimplantation embryo development.

MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants.

The profiles and clinical significance of extraocular muscle-expressed lncRNAs and mRNAs in oculomotor nerve palsy

Introduction

Oculomotor nerve palsy (ONP) arises from primary abnormalities in the central neural pathways that control the extraocular muscles (EOMs). Long non-coding RNAs (lncRNAs) have been found to be involved in the pathogenesis of various neuroparalytic diseases. However, little is known about the role of lncRNAs in ONP.

Methods

We collected medial rectus muscle tissue from ONP and constant exotropia (CXT) patients during strabismus surgeries for RNA sequencing analysis. Differentially expressed mRNAs and lncRNAs were revealed and included in the functional enrichment analysis. Co-expression analysis was conducted between these differentially expressed mRNAs and lncRNAs, followed by target gene prediction of differentially expressed lncRNAs. In addition, lncRNA-microRNA and lncRNA-transcription factor-mRNA interaction networks were constructed to further elaborate the pathological changes in medial rectus muscle of ONP. Furthermore, RT-qPCR was applied to further validate the expression levels of important lncRNAs and mRNAs, whose clinical significance was examined by receiver operating characteristic (ROC) curve analysis.

Results

A total of 618 differentially expressed lncRNAs and 322 differentially expressed mRNAs were identified. The up-regulated mRNAs were significantly related to cholinergic synaptic transmission (such as CHRM3 and CHRND) and the components and metabolism of extracellular matrix (such as CHI3L1 and COL19A1), while the down-regulated mRNAs were significantly correlated with the composition (such as MYH7 and MYL3) and contraction force (such as MYH7 and TNNT1) of muscle fibers. Co-expression analysis and target gene prediction revealed the strong correlation between MYH7 and NR_126491.1 as well as MYOD1 and ENST00000524479. Moreover, the differential expressions of lncRNAs (XR_001739409.1, NR_024160.1 and XR_001738373.1) and mRNAs (CDKN1A, MYOG, MYOD1, MYBPH, TMEM64, STATH, and MYL3) were validated by RT-qPCR. ROC curve analysis showed that lncRNAs (XR_001739409.1, NR_024160.1, and NR_002766.2) and mRNAs (CDKN1A, MYOG, MYOD1, MYBPH, TMEM64, and STATH) might be promising biomarkers of ONP.

Conclusions

These results may shed light on the molecular biology of EOMs of ONP, as well as the possible correlation of lncRNAs and mRNAs with clinical practice.

Mutant GGGGCC RNA prevents YY1 from binding to Fuzzy promoter which stimulates Wnt/β-catenin pathway in C9ALS/FTD

The GGGGCC hexanucleotide repeat expansion mutation in the chromosome 9 open reading frame 72 (C9orf72) gene is a major genetic cause of amyotrophic lateral sclerosis and frontotemporal dementia (C9ALS/FTD). In this study, we demonstrate that the zinc finger (ZF) transcriptional regulator Yin Yang 1 (YY1) binds to the promoter region of the planar cell polarity gene Fuzzy to regulate its transcription. We show that YY1 interacts with GGGGCC repeat RNA via its ZF and that this interaction compromises the binding of YY1 to the FuzzyYY1 promoter sites, resulting in the downregulation of Fuzzy transcription. The decrease in Fuzzy protein expression in turn activates the canonical Wnt/β-catenin pathway and induces synaptic deficits in C9ALS/FTD neurons. Our findings demonstrate a C9orf72 GGGGCC RNA-initiated perturbation of YY1-Fuzzy transcriptional control that implicates aberrant Wnt/β-catenin signalling in C9ALS/FTD-associated neurodegeneration. This pathogenic cascade provides a potential new target for disease-modifying therapy.

Essential transcription factors for induced neuron differentiation

Neurogenins are proneural transcription factors required to specify neuronal identity. Their overexpression in human pluripotent stem cells rapidly produces cortical-like neurons with spiking activity and, because of this, they have been widely adopted for human neuron disease models. However, we do not fully understand the key downstream regulatory effectors responsible for driving neural differentiation. Here, using inducible expression of NEUROG1 and NEUROG2, we identify transcription factors (TFs) required for directed neuronal differentiation by combining expression and chromatin accessibility analyses with a pooled in vitro CRISPR-Cas9 screen targeting all ~1900 TFs in the human genome. The loss of one of these essential TFs (ZBTB18) yields few MAP2-positive neurons. Differentiated ZBTB18-null cells have radically altered gene expression, leading to cytoskeletal defects and stunted neurites and spines. In addition to identifying key downstream TFs for neuronal differentiation, our work develops an integrative multi-omics and TFome-wide perturbation platform to rapidly characterize essential TFs for the differentiation of any human cell type.

Discovering common pathogenetic processes between COVID-19 and tuberculosis by bioinformatics and system biology approach

Introduction

The coronavirus disease 2019 (COVID-19) pandemic, stemming from the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has persistently threatened the global health system. Meanwhile, tuberculosis (TB) caused by Mycobacterium tuberculosis (M. tuberculosis) still continues to be endemic in various regions of the world. There is a certain degree of similarity between the clinical features of COVID-19 and TB, but the underlying common pathogenetic processes between COVID-19 and TB are not well understood.

Methods

To elucidate the common pathogenetic processes between COVID-19 and TB, we implemented bioinformatics and systematic research to obtain shared pathways and molecular biomarkers. Here, the RNA-seq datasets (GSE196822 and GSE126614) are used to extract shared differentially expressed genes (DEGs) of COVID-19 and TB. The common DEGs were used to identify common pathways, hub genes, transcriptional regulatory networks, and potential drugs.

Results

A total of 96 common DEGs were selected for subsequent analyses. Functional enrichment analyses showed that viral genome replication and immune-related pathways collectively contributed to the development and progression of TB and COVID-19. Based on the protein-protein interaction (PPI) network analysis, we identified 10 hub genes, including IFI44L, ISG15, MX1, IFI44, OASL, RSAD2, GBP1, OAS1, IFI6, and HERC5. Subsequently, the transcription factor (TF)-gene interaction and microRNA (miRNA)-gene coregulatory network identified 61 TFs and 29 miRNAs. Notably, we identified 10 potential drugs to treat TB and COVID-19, namely suloctidil, prenylamine, acetohexamide, terfenadine, prochlorperazine, 3'-azido-3'-deoxythymidine, chlorophyllin, etoposide, clioquinol, and propofol.

Conclusion

This research provides novel strategies and valuable references for the treatment of tuberculosis and COVID-19.

Revealing spatial multimodal heterogeneity in tissues with SpaTrio

Capturing and depicting the multimodal tissue information of tissues at the spatial scale remains a significant challenge owing to technical limitations in single-cell multi-omics and spatial transcriptomics sequencing. Here, we developed a computational method called SpaTrio that can build spatial multi-omics data by integrating these two datasets through probabilistic alignment and enabling further analysis of gene regulation and cellular interactions. We benchmarked SpaTrio using simulation datasets and demonstrated its accuracy and robustness. Next, we evaluated SpaTrio on biological datasets and showed that it could detect topological patterns of cells and modalities. SpaTrio has also been applied to multiple sets of actual data to uncover spatially multimodal heterogeneity, understand the spatiotemporal regulation of gene expression, and resolve multimodal communication among cells. Our data demonstrated that SpaTrio could accurately map single cells and reconstruct the spatial distribution of various biomolecules, providing valuable multimodal insights into spatial biology.

Epigenetic modifications regulate cultivar-specific root development and metabolic adaptation to nitrogen availability in wheat

The breeding of crops with improved nitrogen use efficiency (NUE) is crucial for sustainable agriculture, but the involvement of epigenetic modifications remains unexplored. Here, we analyze the chromatin landscapes of two wheat cultivars (KN9204 and J411) that differ in NUE under varied nitrogen conditions. The expression of nitrogen metabolism genes is closely linked to variation in histone modification instead of differences in DNA sequence. Epigenetic modifications exhibit clear cultivar-specificity, which likely contributes to distinct agronomic traits. Additionally, low nitrogen (LN) induces H3K27ac and H3K27me3 to significantly enhance root growth in KN9204, while remarkably inducing NRT2 in J411. Evidence from histone deacetylase inhibitor treatment and transgenic plants with loss function of H3K27me3 methyltransferase shows that changes in epigenetic modifications could alter the strategy preference for root development or nitrogen uptake in response to LN. Here, we show the importance of epigenetic regulation in mediating cultivar-specific adaptation to LN in wheat.

Epicardioid single-cell genomics uncovers principles of human epicardium biology in heart development and disease

The epicardium, the mesothelial envelope of the vertebrate heart, is the source of multiple cardiac cell lineages during embryonic development and provides signals that are essential to myocardial growth and repair. Here we generate self-organizing human pluripotent stem cell-derived epicardioids that display retinoic acid-dependent morphological, molecular and functional patterning of the epicardium and myocardium typical of the left ventricular wall. By combining lineage tracing, single-cell transcriptomics and chromatin accessibility profiling, we describe the specification and differentiation process of different cell lineages in epicardioids and draw comparisons to human fetal development at the transcriptional and morphological levels. We then use epicardioids to investigate the functional cross-talk between cardiac cell types, gaining new insights into the role of IGF2/IGF1R and NRP2 signaling in human cardiogenesis. Finally, we show that epicardioids mimic the multicellular pathogenesis of congenital or stress-induced hypertrophy and fibrotic remodeling. As such, epicardioids offer a unique testing ground of epicardial activity in heart development, disease and regeneration.

Multimodal spatiotemporal phenotyping of human retinal organoid development

Organoids generated from human pluripotent stem cells provide experimental systems to study development and disease, but quantitative measurements across different spatial scales and molecular modalities are lacking. In this study, we generated multiplexed protein maps over a retinal organoid time course and primary adult human retinal tissue. We developed a toolkit to visualize progenitor and neuron location, the spatial arrangements of extracellular and subcellular components and global patterning in each organoid and primary tissue. In addition, we generated a single-cell transcriptome and chromatin accessibility timecourse dataset and inferred a gene regulatory network underlying organoid development. We integrated genomic data with spatially segmented nuclei into a multimodal atlas to explore organoid patterning and retinal ganglion cell (RGC) spatial neighborhoods, highlighting pathways involved in RGC cell death and showing that mosaic genetic perturbations in retinal organoids provide insight into cell fate regulation.

Antiviral innate immune memory in alveolar macrophages following SARS-CoV-2 infection

Pathogen encounter results in long-lasting epigenetic imprinting that shapes diseases caused by heterologous pathogens. The breadth of this innate immune memory is of particular interest in the context of respiratory pathogens with increased pandemic potential and wide-ranging impact on global health. Here, we investigated epigenetic imprinting across cell lineages in a disease relevant murine model of SARS-CoV-2 recovery. Past SARS-CoV-2 infection resulted in increased chromatin accessibility of type I interferon (IFN-I) related transcription factors in airway-resident macrophages. Mechanistically, establishment of this innate immune memory required viral pattern recognition and canonical IFN-I signaling and augmented secondary antiviral responses. Past SARS-CoV-2 infection ameliorated disease caused by the heterologous respiratory pathogen influenza A virus. Insights into innate immune memory and how it affects subsequent infections with heterologous pathogens to influence disease pathology could facilitate the development of broadly effective therapeutic strategies.

Epigenetic regulation during cancer transitions across 11 tumour types

Chromatin accessibility is essential in regulating gene expression and cellular identity, and alterations in accessibility have been implicated in driving cancer initiation, progression and metastasis1-4. Although the genetic contributions to oncogenic transitions have been investigated, epigenetic drivers remain less understood. Here we constructed a pan-cancer epigenetic and transcriptomic atlas using single-nucleus chromatin accessibility data (using single-nucleus assay for transposase-accessible chromatin) from 225 samples and matched single-cell or single-nucleus RNA-sequencing expression data from 206 samples. With over 1 million cells from each platform analysed through the enrichment of accessible chromatin regions, transcription factor motifs and regulons, we identified epigenetic drivers associated with cancer transitions. Some epigenetic drivers appeared in multiple cancers (for example, regulatory regions of ABCC1 and VEGFA; GATA6 and FOX-family motifs), whereas others were cancer specific (for example, regulatory regions of FGF19, ASAP2 and EN1, and the PBX3 motif). Among epigenetically altered pathways, TP53, hypoxia and TNF signalling were linked to cancer initiation, whereas oestrogen response, epithelial-mesenchymal transition and apical junction were tied to metastatic transition. Furthermore, we revealed a marked correlation between enhancer accessibility and gene expression and uncovered cooperation between epigenetic and genetic drivers. This atlas provides a foundation for further investigation of epigenetic dynamics in cancer transitions.

VDR promotes testosterone synthesis in mouse Leydig cells via regulation of cholesterol side chain cleavage cytochrome P450 (Cyp11a1) expression

BACKGROUND

The vitamin D receptor (VDR) mediates the pleiotropic biological actions that include osteoporosis, immune responses and androgen synthesis.VDR is widely expressed in testis cells such as Leydig cells, Sertoli cells, and sperm. The levels of steroids are critical for sexual development. In the early stage of steroidogenesis, cholesterol is converted to pregnenolone (precursor of most steroid hormones) by cholesterol side-chain lyase (CYP11A1), which eventually synthesizes the male hormone testosterone.

OBJECTIVE

This study aims to reveal how VDR regulates CYP11A1 expression and affects testosterone synthesis in murine Leydig cells.

METHODS

The levels of VDR, CYP11A1 were determined by quantitative real-time polymerase chain reaction (RT-qPCR) or western blot. Targeted relationship between VDR and Cyp11a1 was evaluated by dual-luciferase reporter assay. The levels of testosterone concentrations in cell culture media serum by enzyme-linked immunosorbent assay (ELISA).

RESULTS

Phylogenetic and motif analysis showed that the Cyp11a1 family had sequence loss, which may have special biological functions during evolution. The results of promoter prediction showed that vitamin D response element (VDRE) existed in the upstream promoter region of murine Cyp11a1. Dual-luciferase assay confirmed that VDR could bind candidate VDREs in upstream region of Cyp11a1, and enhance gene expression. Tissue distribution and localizatio analysis showed that Cyp11a1 was mainly expressed in testis, and dominantly existed in murine Leydig cells. Furthermore, over-expression VDR and CYP11A1 significantly increased testosterone synthesis in mice Leydig cells.

CONCLUSIONS

Active vitamin D3 (VD3) and Vdr interference treatment showed that VD3/VDR had a positive regulatory effect on Cyp11a1 expression and testosterone secretion. VDR promotes testosterone synthesis in male mice by up-regulating Cyp11a1 expression, which played an important role for male reproduction.

GATA transcription factors drive initial Xist upregulation after fertilization through direct activation of long-range enhancers

X-chromosome inactivation (XCI) balances gene expression between the sexes in female mammals. Shortly after fertilization, upregulation of Xist RNA from one X chromosome initiates XCI, leading to chromosome-wide gene silencing. XCI is maintained in all cell types, except the germ line and the pluripotent state where XCI is reversed. The mechanisms triggering Xist upregulation have remained elusive. Here we identify GATA transcription factors as potent activators of Xist. Through a pooled CRISPR activation screen in murine embryonic stem cells, we demonstrate that GATA1, as well as other GATA transcription factors can drive ectopic Xist expression. Moreover, we describe GATA-responsive regulatory elements in the Xist locus bound by different GATA factors. Finally, we show that GATA factors are essential for XCI induction in mouse preimplantation embryos. Deletion of GATA1/4/6 or GATA-responsive Xist enhancers in mouse zygotes effectively prevents Xist upregulation. We propose that the activity or complete absence of various GATA family members controls initial Xist upregulation, XCI maintenance in extra-embryonic lineages and XCI reversal in the epiblast.

Peroxiredoxin II regulates exosome secretion from dermal mesenchymal stem cells through the ISGylation signaling pathway

BACKGROUND

Exosomes are small extracellular vesicles that play important roles in intercellular communication and have potential therapeutic applications in regenerative medicine. Dermal mesenchymal stem cells (DMSCs) are a promising source of exosomes due to their regenerative and immunomodulatory properties. However, the molecular mechanisms regulating exosome secretion from DMSCs are not fully understood.

RESULTS

In this study, the role of peroxiredoxin II (Prx II) in regulating exosome secretion from DMSCs and the underlying molecular mechanisms were investigated. It was discovered that depletion of Prx II led to a significant reduction in exosome secretion from DMSCs and an increase in the number of intracellular multivesicular bodies (MVBs), which serve as precursors of exosomes. Mechanistically, Prx II regulates the ISGylation switch that controls MVB degradation and impairs exosome secretion. Specifically, Prx II depletion decreased JNK activity, reduced the expression of the transcription inhibitor Foxo1, and promoted miR-221 expression. Increased miR-221 expression inhibited the STAT signaling pathway, thus downregulating the expression of ISGylation-related genes involved in MVB degradation. Together, these results identify Prx II as a critical regulator of exosome secretion from DMSCs through the ISGylation signaling pathway.

CONCLUSIONS

Our findings provide important insights into the molecular mechanisms regulating exosome secretion from DMSCs and highlight the critical role of Prx II in controlling the ISGylation switch that regulates DMSC-exosome secretion. This study has significant implications for developing new therapeutic strategies in regenerative medicine. Video Abstract.

Single-cell chromatin accessibility and transcriptomic characterization of Behcet's disease

Behect's disease is a chronic vasculitis characterized by complex multi-organ immune aberrations. However, a comprehensive understanding of the gene-regulatory profile of peripheral autoimmunity and the diverse immune responses across distinct cell types in Behcet's disease (BD) is still lacking. Here, we present a multi-omic single-cell study of 424,817 cells in BD patients and non-BD individuals. This study maps chromatin accessibility and gene expression in the same biological samples, unraveling vast cellular heterogeneity. We identify widespread cell-type-specific, disease-associated active and pro-inflammatory immunity in both transcript and epigenomic aspects. Notably, integrative multi-omic analysis reveals putative TF regulators that might contribute to chromatin accessibility and gene expression in BD. Moreover, we predicted gene-regulatory networks within nominated TF activators, including AP-1, NF-kB, and ETS transcript factor families, which may regulate cellular interaction and govern inflammation. Our study illustrates the epigenetic and transcriptional landscape in BD peripheral blood and expands understanding of potential epigenomic immunopathology in this disease.

KDeep: a new memory-efficient data extraction method for accurately predicting DNA/RNA transcription factor binding sites

This paper addresses the crucial task of identifying DNA/RNA binding sites, which has implications in drug/vaccine design, protein engineering, and cancer research. Existing methods utilize complex neural network structures, diverse input types, and machine learning techniques for feature extraction. However, the growing volume of sequences poses processing challenges. This study introduces KDeep, employing a CNN-LSTM architecture with a novel encoding method called 2Lk. 2Lk enhances prediction accuracy, reduces memory consumption by up to 84%, reduces trainable parameters, and improves interpretability by approximately 79% compared to state-of-the-art approaches. KDeep offers a promising solution for accurate and efficient binding site prediction.

Core clock gene BMAL1 and RNA-binding protein MEX3A collaboratively regulate Lgr5 expression in intestinal crypt cells

The intestinal epithelium is highly regenerative. Rapidly proliferating LGR5+ crypt base columnar (CBC) cells are responsible for epithelial turnover needed to maintain intestinal homeostasis. Upon tissue damage, loss of LGR5+ CBCs can be compensated by activation of quiescent +4 intestinal stem cells (ISCs) or early progenitor cells to restore intestinal regeneration. LGR5+ CBC self-renewal and ISC conversion to LGR5+ cells are regulated by external signals originating from the ISC niche. In contrast, little is known about intrinsic regulatory mechanisms critical for maintenance of LGR5+ CBC homeostasis. We found that LGR5 expression in intestinal crypt cells is controlled by the circadian core clock gene BMAL1 and the BMAL1-regulated RNA-binding protein MEX3A. BMAL1 directly activated transcription of Mex3a. MEX3A in turn bound to and stabilized Lgr5 mRNA. Bmal1 depletion reduced Mex3a and Lgr5 expression and led to increased ferroptosis, which consequently decreased LGR5+ CBC numbers and increased the number of crypt cells expressing +4 ISC marker BMI1. Together, these findings reveal a BMAL1-centered intrinsic regulatory pathway that maintains LGR5 expression in the crypt cells and suggest a potential mechanism contributing to ISC homeostasis.

Super enhancer-driven core transcriptional regulatory circuitry crosstalk with cancer plasticity and patient mortality in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is a clinically aggressive subtype of breast cancer. Core transcriptional regulatory circuitry (CRC) consists of autoregulated transcription factors (TFs) and their enhancers, which dominate gene expression programs and control cell fate. However, there is limited knowledge of CRC in TNBC. Herein, we systemically characterized the activated super-enhancers (SEs) and interrogated 14 CRCs in breast cancer. We found that CRCs could be broadly involved in DNA conformation change, metabolism process, and signaling response affecting the gene expression reprogramming. Furthermore, these CRC TFs are capable of coordinating with partner TFs bridging the enhancer-promoter loops. Notably, the CRC TF and partner pairs show remarkable specificity for molecular subtypes of breast cancer, especially in TNBC. USF1, SOX4, and MYBL2 were identified as the TNBC-specific CRC TFs. We further demonstrated that USF1 was a TNBC immunophenotype-related TF. Our findings that the rewiring of enhancer-driven CRCs was related to cancer immune and mortality, will facilitate the development of epigenetic anti-cancer treatment strategies.

Integrative single-nucleus multi-omics analysis prioritizes candidate cis and trans regulatory networks and their target genes in Alzheimer's disease brains

BACKGROUND

The genetic underpinnings of late-onset Alzheimer's disease (LOAD) are yet to be fully elucidated. Although numerous LOAD-associated loci have been discovered, the causal variants and their target genes remain largely unknown. Since the brain is composed of heterogenous cell subtypes, it is imperative to study the brain on a cell subtype specific level to explore the biological processes underlying LOAD.

METHODS

Here, we present the largest parallel single-nucleus (sn) multi-omics study to simultaneously profile gene expression (snRNA-seq) and chromatin accessibility (snATAC-seq) to date, using nuclei from 12 normal and 12 LOAD brains. We identified cell subtype clusters based on gene expression and chromatin accessibility profiles and characterized cell subtype-specific LOAD-associated differentially expressed genes (DEGs), differentially accessible peaks (DAPs) and cis co-accessibility networks (CCANs).

RESULTS

Integrative analysis defined disease-relevant CCANs in multiple cell subtypes and discovered LOAD-associated cell subtype-specific candidate cis regulatory elements (cCREs), their candidate target genes, and trans-interacting transcription factors (TFs), some of which, including ELK1, JUN, and SMAD4 in excitatory neurons, were also LOAD-DEGs. Finally, we focused on a subset of cell subtype-specific CCANs that overlap known LOAD-GWAS regions and catalogued putative functional SNPs changing the affinities of TF motifs within LOAD-cCREs linked to LOAD-DEGs, including APOE and MYO1E in a specific subtype of microglia and BIN1 in a subpopulation of oligodendrocytes.

CONCLUSIONS

To our knowledge, this study represents the most comprehensive systematic interrogation to date of regulatory networks and the impact of genetic variants on gene dysregulation in LOAD at a cell subtype resolution. Our findings reveal crosstalk between epigenetic, genomic, and transcriptomic determinants of LOAD pathogenesis and define catalogues of candidate genes, cCREs, and variants involved in LOAD genetic etiology and the cell subtypes in which they act to exert their pathogenic effects. Overall, these results suggest that cell subtype-specific cis-trans interactions between regulatory elements and TFs, and the genes dysregulated by these networks contribute to the development of LOAD.

PD-1 instructs a tumor-suppressive metabolic program that restricts glycolysis and restrains AP-1 activity in T cell lymphoma

The PDCD1-encoded immune checkpoint receptor PD-1 is a key tumor suppressor in T cells that is recurrently inactivated in T cell non-Hodgkin lymphomas (T-NHLs). The highest frequencies of PDCD1 deletions are detected in advanced disease, predicting inferior prognosis. However, the tumor-suppressive mechanisms of PD-1 signaling remain unknown. Here, using tractable mouse models for T-NHL and primary patient samples, we demonstrate that PD-1 signaling suppresses T cell malignancy by restricting glycolytic energy and acetyl coenzyme A (CoA) production. In addition, PD-1 inactivation enforces ATP citrate lyase (ACLY) activity, which generates extramitochondrial acetyl-CoA for histone acetylation to enable hyperactivity of activating protein 1 (AP-1) transcription factors. Conversely, pharmacological ACLY inhibition impedes aberrant AP-1 signaling in PD-1-deficient T-NHLs and is toxic to these cancers. Our data uncover genotype-specific vulnerabilities in PDCD1-mutated T-NHL and identify PD-1 as regulator of AP-1 activity.

Human-genome single nucleotide polymorphisms affecting transcription factor binding and their role in pathogenesis

Single nucleotide polymorphisms (SNPs) are the most common type of variation in the human genome. The vast majority of SNPs identified in the human genome do not have any effect on the phenotype; however, some can lead to changes in the function of a gene or the level of its expression. Most SNPs associated with certain traits or pathologies are mapped to regulatory regions of the genome and affect gene expression by changing transcription factor binding sites. In recent decades, substantial effort has been invested in searching for such regulatory SNPs (rSNPs) and understanding the mechanisms by which they lead to phenotypic differences, primarily to individual differences in susceptibility to diseases and in sensitivity to drugs. The development of the NGS (next-generation sequencing) technology has contributed not only to the identification of a huge number of SNPs and to the search for their association (genome-wide association studies, GWASs) with certain diseases or phenotypic manifestations, but also to the development of more productive approaches to their functional annotation. It should be noted that the presence of an association does not allow one to identify a functional, truly disease-associated DNA sequence variant among multiple marker SNPs that are detected due to linkage disequilibrium. Moreover, determination of associations of genetic variants with a disease does not provide information about the functionality of these variants, which is necessary to elucidate the molecular mechanisms of the development of pathology and to design effective methods for its treatment and prevention. In this regard, the functional analysis of SNPs annotated in the GWAS catalog, both at the genome-wide level and at the level of individual SNPs, became especially relevant in recent years. A genome-wide search for potential rSNPs is possible without any prior knowledge of their association with a trait. Thus, mapping expression quantitative trait loci (eQTLs) makes it possible to identify an SNP for which - among transcriptomes of homozygotes and heterozygotes for its various alleles - there are differences in the expression level of certain genes, which can be located at various distances from the SNP. To predict rSNPs, approaches based on searches for allele-specific events in RNA-seq, ChIP-seq, DNase-seq, ATAC-seq, MPRA, and other data are also used. Nonetheless, for a more complete functional annotation of such rSNPs, it is necessary to establish their association with a trait, in particular, with a predisposition to a certain pathology or sensitivity to drugs. Thus, approaches to finding SNPs important for the development of a trait can be categorized into two groups: (1) starting from data on an association of SNPs with a certain trait, (2) starting from the determination of allele-specific changes at the molecular level (in a transcriptome or regulome). Only comprehensive use of strategically different approaches can considerably enrich our knowledge about the role of genetic determinants in the molecular mechanisms of trait formation, including predisposition to multifactorial diseases.

Mild sleep restriction increases endothelial oxidative stress in female persons

Sleep restriction is associated with increased cardiovascular risk, which is more pronounced in female than male persons. We reported recently first causal evidence that mild, prolonged sleep restriction mimicking "real-life" conditions impairs endothelial function, a key step in the development and progression of cardiovascular disease, in healthy female persons. However, the underlying mechanisms are unclear. In model organisms, sleep restriction increases oxidative stress and upregulates antioxidant response via induction of the antioxidant regulator nuclear factor (erythroid-derived 2)-like 2 (Nrf2). Here, we assessed directly endothelial cell oxidative stress and antioxidant responses in healthy female persons (n = 35) after 6 weeks of mild sleep restriction (1.5 h less than habitual sleep) using randomized crossover design. Sleep restriction markedly increased endothelial oxidative stress without upregulating antioxidant response. Using RNA-seq and a predicted protein-protein interaction database, we identified reduced expression of endothelial Defective in Cullin Neddylation-1 Domain Containing 3 (DCUN1D3), a protein that licenses Nrf2 antioxidant responses, as a mediator of impaired endothelial antioxidant response in sleep restriction. Thus, sleep restriction impairs clearance of endothelial oxidative stress that over time increases cardiovascular risk.Trial Registration: NCT02835261 .

Seizure enhances SUMOylation and zinc-finger transcriptional repression in neuronal nuclei

A single episode of pilocarpine-induced status epilepticus can trigger the development of spontaneous recurrent seizures in a rodent model for epilepsy. The initial seizure-induced events in neuronal nuclei that lead to long-term changes in gene expression and cellular responses likely contribute toward epileptogenesis. Using a transgenic mouse model to specifically isolate excitatory neuronal nuclei, we profiled the seizure-induced nuclear proteome via tandem mass tag mass spectrometry and observed robust enrichment of nuclear proteins associated with the SUMOylation pathway. In parallel with nuclear proteome, we characterized nuclear gene expression by RNA sequencing which provided insights into seizure-driven transcriptional regulation and dynamics. Strikingly, we saw widespread downregulation of zinc-finger transcription factors, specifically proteins that harbor Krüppel-associated box (KRAB) domains. Our results provide a detailed snapshot of nuclear events induced by seizure activity and demonstrate a robust method for cell-type-specific nuclear profiling that can be applied to other cell types and models.