Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities

PRDM16 determines specification of ventricular cardiomyocytes by suppressing alternative cell fates

PRDM16 is a transcription factor with histone methyltransferase activity expressed at the earliest stages of cardiac development. Pathogenic mutations in humans lead to cardiomyopathy, conduction abnormalities, and heart failure. PRDM16 is specifically expressed in ventricular but not atrial cardiomyocytes, and its expression declines postnatally. Because in other tissues PRDM16 is best known for its role in binary cell fate decisions, we hypothesized a similar decision-making function in cardiomyocytes. Here, we demonstrated that cardiomyocyte-specific deletion of Prdm16 during cardiac development results in contractile dysfunction and abnormal electrophysiology of the postnatal heart, resulting in premature death. By combined RNA+ATAC single-cell sequencing, we found that PRDM16 favors ventricular working cardiomyocyte identity, by opposing the activity of master regulators of ventricular conduction and atrial fate. Myocardial loss of PRDM16 during development resulted in hyperplasia of the (distal) ventricular conduction system. Hence, PRDM16 plays an indispensable role during cardiac development by driving ventricular working cardiomyocyte identity.

Epigenetic changes in pyloric caeca of Atlantic salmon fed diets containing increasing levels of lipids and choline

An earlier study of ours investigating the effect of dietary lipid levels on the choline requirement of Atlantic salmon showed increasing severity of intestinal steatosis with increasing lipid levels. As choline is involved in epigenetic regulation by being the key methyl donor, pyloric caeca samples from the study were analysed for epigenetic effects of dietary lipid and choline levels. The diets varied in lipid levels between 16% and 28%, and choline levels between 1.9 and 2.3 g/kg. The diets were fed for 8 weeks to Atlantic salmon of 25 g of initial weight. Using reduced representation bisulfite sequencing (RRBS), this study revealed that increasing dietary lipid levels induced methylation differences in genes involved in membrane transport and signalling pathways, and in microRNAs important for the regulation of lipid homoeostasis. Increasing choline levels also affected genes involved in fatty acid biosynthesis and transport, lipolysis, and lipogenesis, as well as important immune genes. Our observations confirmed that choline is involved in epigenetic regulation in Atlantic salmon, as has been reported for higher vertebrates. This study showed the need for the inclusion of biomarkers of epigenetic processes in studies that must be conducted to define optimal choline levels in diets for Atlantic salmon.

N6-methyladenosine methylation analysis of long noncoding RNAs and mRNAs in 5-FU-resistant colon cancer cells

N6 methyladenosine (m6A), methylation at the sixth N atom of adenosine, is the most common and abundant modification in mammalian mRNAs and non-coding RNAs. Increasing evidence shows that the alteration of m6A modification level could regulate tumour proliferation, metastasis, self-renewal, and immune infiltration by regulating the related expression of tumour genes. However, the role of m6A modification in colorectal cancer (CRC) drug resistance is unclear. Here, MeRIP-seq and RNA-seq techniques were utilized to obtain mRNA, lncRNA expression, and their methylation profiles in 5-Fluorouracil (5-FU)-resistant colon cancer HCT-15 cells and control cells. In addition, we performed detailed bioinformatics analysis as well as in vitro experiments of lncRNA to explore the function of lncRNA with differential m6A in CRC progression and drug resistance. In this study, we obtained the m6A methylomic landscape of CRC cells and resistance group cells by MeRIP-seq and RNA-seq. We identified 3698 differential m6A peaks, of which 2224 were hypermethylated, and 1474 were hypomethylated. Among the lncRNAs, 60 were hypermethylated, and 38 were hypomethylated. GO and KEGG analysis annotations showed significant enrichment of endocytosis and MAPK signalling pathways. Moreover, knockdown of lncRNA ADIRF-AS1 and AL139035.1 promoted CRC proliferation and invasive metastasis in vitro. lncRNA- mRNA network showed that ADIRF-AS1 and AL139035.1 May play a key role in regulating drug resistance formation. We provide the first m6A methylation profile in 5-FU resistance CRC cells and analyse the functions of differential m6A-modified mRNAs and lncRNAs. Our results indicated that differential m6A RNAs were significantly associated with MAPK signalling and endocytosis after induction of 5-FU resistance. Knockdown of LncRNA ADIRF-AS1 and AL139035.1 promotes CRC progression and might be critical in regulating drug resistance formation.

Benchmarking tools for transcription factor prioritization

Spatiotemporal regulation of gene expression is controlled by transcription factor (TF) binding to regulatory elements, resulting in a plethora of cell types and cell states from the same genetic information. Due to the importance of regulatory elements, various sequencing methods have been developed to localise them in genomes, for example using ChIP-seq profiling of the histone mark H3K27ac that marks active regulatory regions. Moreover, multiple tools have been developed to predict TF binding to these regulatory elements based on DNA sequence. As altered gene expression is a hallmark of disease phenotypes, identifying TFs driving such gene expression programs is critical for the identification of novel drug targets. In this study, we curated 84 chromatin profiling experiments (H3K27ac ChIP-seq) where TFs were perturbed through e.g., genetic knockout or overexpression. We ran nine published tools to prioritize TFs using these real-world datasets and evaluated the performance of the methods in identifying the perturbed TFs. This allowed the nomination of three frontrunner tools, namely RcisTarget, MEIRLOP and monaLisa. Our analyses revealed opportunities and commonalities of tools that will help to guide further improvements and developments in the field.

Mechanistic insights into zinc oxide nanoparticles induced embryotoxicity via H3K9me3 modulation

The widespread application of nanoparticles (NPs) in various fields has raised health concerns, especially in reproductive health. Our research has shown zinc oxide nanoparticles (ZnONPs) exhibit the most significant toxicity to pre-implantation embryos in mice compared to other common NPs. In patients undergoing assisted reproduction technology (ART), a significant negative correlation was observed between Zn concentration and clinical outcomes. Therefore, this study explores the impact of ZnONPs exposure on pre-implantation embryonic development and its underlying mechanisms. We revealed that both in vivo and in vitro exposure to ZnONPs impairs pre-implantation embryonic development. Moreover, ZnONPs were found to reduce the pluripotency of mouse embryonic stem cells (mESCs), as evidenced by teratoma and diploid chimera assays. Employing multi-omics approaches, including RNA-Seq, CUT&Tag, and ATAC-seq, the embryotoxicity mechanisms of ZnONPs were elucidated. The findings indicate that ZnONPs elevate H3K9me3 levels, leading to increased heterochromatin and consequent inhibition of gene expression related to development and pluripotency. Notably, Chaetocin, a H3K9me3 inhibitor, sucessfully reversed the embryotoxicity effects induced by ZnONPs. Additionally, the direct interaction between ZnONPs and H3K9me3 was verified through pull-down and immunoprecipitation assays. Collectively, these findings offer new insights into the epigenetic mechanisms of ZnONPs toxicity, enhancing our understanding of their impact on human reproductive health.

Folate-mediated transgenerational inheritance of sperm DNA methylation patterns correlate with spinal axon regeneration

In mammals, the molecular mechanisms underlying transgenerational inheritance of phenotypic traits in serial generations of progeny after ancestral environmental exposures, without variation in DNA sequence, remain elusive. We've recently described transmission of a beneficial trait in rats and mice, in which F0 supplementation of methyl donors, including folic acid, generates enhanced axon regeneration after sharp spinal cord injury in untreated F1 to F3 progeny linked to differential DNA methylation levels in spinal cord tissue. To test whether the transgenerational effect of folic acid is transmitted via the germline, we performed whole-genome methylation sequencing on sperm DNA from F0 mice treated with either folic acid or vehicle control, and their F1, F2, and F3 untreated progeny. Transgenerational differentially methylated regions (DMRs) are observed in each consecutive generation and distinguish folic acid from untreated lineages, predominate outside of CpG islands and in regions of the genome that regulate gene expression, including promoters, and overlap at both the differentially methylated position (DMP) and gene levels. These findings indicate that molecular changes between generations are caused by ancestral folate supplementation. In addition, 29,719 DMPs exhibit serial increases or decreases in DNA methylation levels in successive generations of untreated offspring, correlating with a serial increase in the phenotype across generations, consistent with a 'wash-in' effect. Sibship-specific DMPs annotate to genes that participate in axon- and synapse-related pathways.

SMARCA5-mediated chromatin remodeling is required for germinal center formation

The establishment of long-lasting immunity against pathogens is facilitated by the germinal center (GC) reaction, during which B cells increase their antibody affinity and differentiate into antibody-secreting cells (ASC) and memory cells. These events involve modifications in chromatin packaging that orchestrate the profound restructuring of gene expression networks that determine cell fate. While several chromatin remodelers were implicated in lymphocyte functions, less is known about SMARCA5. Here, using ribosomal pull-down for analyzing translated genes in GC B cells, coupled with functional experiments in mice, we identified SMARCA5 as a key chromatin remodeler in B cells. While the naive B cell compartment remained unaffected following conditional depletion of Smarca5, effective proliferation during B cell activation, immunoglobulin class switching, and as a result GC formation and ASC differentiation were impaired. Single-cell multiomic sequencing analyses revealed that SMARCA5 is crucial for facilitating the transcriptional modifications and genomic accessibility of genes that support B cell activation and differentiation. These findings offer novel insights into the functions of SMARCA5, which can be targeted in various human pathologies.

Inhibition of PRMT5 moderately suppresses prostate cancer growth in vivo but enhances its response to immunotherapy

Protein arginine methylation is a common post-translational modification (PTM) catalyzed by nine protein arginine methyltransferases (PRMTs). As the major symmetric arginine methyltransferase that methylates both histone and non-histone substrates, PRMT5 plays key roles in a number of biological processes critical for development and tumorigenesis. PRMT5 overexpression has been reported in multiple cancer types including prostate cancer (PCa), but the exact biological and mechanistic understanding of PRMT5 in aggressive PCa remains ill-defined. Here, we show that PRMT5 is upregulated in PCa, correlates with worse patient survival, promotes corrupted RNA splicing, and functionally cooperates with an array of pro-tumorigenic pathways to enhance oncogenesis. PRMT5 inhibition via either genetic knockdown or pharmacological inhibition reduces stemness with paralleled differentiation and arrests cell cycle progression without causing appreciable apoptosis. Strikingly, the severity of antitumor effect of PRMT5 inhibition correlates with disease aggressiveness, with AR+ PCa being less affected. Molecular characterization pinpoints MYC, but not (or at least to a lesser degree) AR, as the main partner of PRMT5 to form a positive feedback loop to exacerbate malignancy in both AR+ and AR- PCa cells. Inspired by the surprising finding that PRMT5 negatively correlates with tumor immune infiltration and transcriptionally suppresses an immune-gene program, we further show that although PRMT5 inhibitor (PRMT5i) EPZ015666 or anti-PD-1 immunotherapy alone exhibits limited antitumor effects, combination of PRMT5i with anti-PD-1 displays superior efficacy in inhibiting castration-resistant PCa (CRPC) in vivo. Finally, to expand the potential use of PRMT5i through a synthetic lethality concept, we also perform a global CRISPR/Cas9 knockout screen to unravel that many clinical-grade drugs of known oncogenic pathways can be repurposed to target CRPC when used in combination with PRMT5i at low doses. Collectively, our findings establish a rationale to exploit PRMT5i in combination with immunotherapy or other targeted therapies to treat aggressive PCa.

Enhancer landscape of lung neuroendocrine tumors reveals regulatory and developmental signatures with potential theranostic implications

Well-differentiated low-grade lung neuroendocrine tumors (lung carcinoids or LNETs) are histopathologically classified as typical and atypical LNETs, but each subtype is still heterogeneous at both the molecular level and its clinical manifestation. Here, we report genome-wide profiles of primary LNETs' cis-regulatory elements by H3K27ac ChIP-seq with matching RNA-seq profiles. Analysis of these regulatory landscapes revealed three regulatory subtypes, independent of the typical/atypical classification. We identified unique differentiation signals that delineate each subtype. The "proneuronal" subtype emerges under the influence of ASCL1, SOX4, and TCF4 transcription factors, embodying a pronounced proneuronal signature. The "luminal-like" subtype is characterized by gain of acetylation at markers of luminal cells and GATA2 activation and loss of LRP5 and OTP. The "HNF+" subtype is characterized by a robust enhancer landscape driven by HNF1A, HNF4A, and FOXA3, with notable acetylation and expression of FGF signaling genes, especially FGFR3 and FGFR4, pivotal components of the FGF pathway. Our findings not only deepen the understanding of LNETs' regulatory and developmental diversity but also spotlight the HNF+ subtype's reliance on FGFR signaling. We demonstrate that targeting this pathway with FGF inhibitors curtails tumor growth both in vitro and in xenograft models, unveiling a potential vulnerability and paving the way for targeted therapies. Overall, our work provides an important resource for studying LNETs to reveal regulatory networks, differentiation signals, and therapeutically relevant dependencies.

Transcriptional repression and enhancer decommissioning silence cell cycle genes in postmitotic tissues

The mechanisms that maintain a non-cycling status in postmitotic tissues are not well understood. Many cell cycle genes have promoters and enhancers that remain accessible even when cells are terminally differentiated and in a non-cycling state, suggesting their repression must be maintained long term. In contrast, enhancer decommissioning has been observed for rate-limiting cell cycle genes in the Drosophila wing, a tissue where the cells die soon after eclosion, but it has been unclear if this also occurs in other contexts of terminal differentiation. In this study, we show that enhancer decommissioning also occurs at specific, rate-limiting cell cycle genes in the long-lived tissues of the Drosophila eye and brain, and we propose this loss of chromatin accessibility may help maintain a robust postmitotic state. We examined the decommissioned enhancers at specific rate-limiting cell cycle genes and showed that they encode for dynamic temporal and spatial expression patterns that include shared, as well as tissue-specific elements, resulting in broad gene expression with developmentally controlled temporal regulation. We extend our analysis to cell cycle gene expression and chromatin accessibility in the mammalian retina using a published dataset and find that the principles of cell cycle gene regulation identified in terminally differentiating Drosophila tissues are conserved in the differentiating mammalian retina. We propose a robust, non-cycling status is maintained in long-lived postmitotic tissues through a combination of stable repression at most cell cycle genes, alongside enhancer decommissioning at specific rate-limiting cell cycle genes.

Loss of ARID3A perturbs intestinal epithelial proliferation-differentiation ratio and regeneration

Intestinal stem cells at the crypt divide and give rise to progenitor cells that proliferate and differentiate into various mature cell types in the transit-amplifying (TA) zone. Here, we showed that the transcription factor ARID3A regulates intestinal epithelial cell proliferation and differentiation at the TA progenitors. ARID3A forms an expression gradient from the villus tip to the upper crypt mediated by TGF-β and WNT. Intestinal-specific deletion of Arid3a reduces crypt proliferation, predominantly in TA cells. Bulk and single-cell transcriptomic analysis shows increased enterocyte and reduced secretory differentiation in the Arid3a cKO intestine, accompanied by enriched upper-villus gene signatures of both cell lineages. We find that the enhanced epithelial differentiation in the Arid3a-deficient intestine is caused by increased binding and transcription of HNF1 and HNF4. Finally, we show that loss of Arid3a impairs irradiation-induced regeneration with sustained cell death and reprogramming. Our findings imply that Arid3a functions to fine-tune the proliferation-differentiation dynamics at the TA progenitors, which are essential for injury-induced regeneration.

SETD1B mutations confer apoptosis resistance and BCL2 independence in B cell lymphoma

The translocation t(14;18) activates BCL2 and is considered the initiating genetic lesion in most follicular lymphomas (FL). Surprisingly, FL patients fail to respond to the BCL2 inhibitor, Venetoclax. We show that mutations and deletions affecting the histone lysine methyltransferase SETD1B (KMT2G) occur in 7% of FLs and 16% of diffuse large B cell lymphomas (DLBCL). Deficiency in SETD1B confers striking resistance to Venetoclax and an experimental MCL-1 inhibitor. SETD1B also acts as a tumor suppressor and cooperates with the loss of KMT2D in lymphoma development in vivo. Consistently, loss of SETD1B in human lymphomas typically coincides with loss of KMT2D. Mechanistically, SETD1B is required for the expression of several proapoptotic BCL2 family proteins. Conversely, inhibitors of the KDM5 histone H3K4 demethylases restore BIM and BIK expression and synergize with Venetoclax in SETD1B-deficient lymphomas. These results establish SETD1B as an epigenetic regulator of cell death and reveal a pharmacological strategy to augment Venetoclax sensitivity in lymphoma.

Stabilization of RRBP1 mRNA via an m6A-dependent manner in prostate cancer constitutes a therapeutic vulnerability amenable to small-peptide inhibition of METTL3

Mounting evidence has implicated the RNA m6A methylation catalyzed by METTL3 in a wide range of physiological and pathological processes, including tumorigenesis. The detailed m6A landscape and molecular mechanism of METTL3 in prostate cancer (PCa) remains ill-defined. We find that METTL3 is overexpressed in PCa and correlates with worse patient survival. Functional studies establish METTL3 as an oncoprotein dependent on its m6A enzymatic activity in both AR+ and AR- PCa cells. To dissect the regulatory network of m6A pathway in PCa, we map the m6A landscape in clinical tumor samples using m6A-seq and identify genome-wide METTL3-binding transcripts via RIP-seq. Mechanistically, we discover RRBP1 as a direct METTL3 target in which METTL3 stabilizes RRBP1 mRNA in an m6A-dependent manner. RRBP1 positively correlates with METTL3 expression in PCa cohorts and exerts an oncogenic role in aggressive PCa cells. Leveraging the 3D structural protein-protein interaction between METTL3 and METTL14, we successfully develop two potential METTL3 peptide inhibitors (RM3 and RSM3) that effectively suppress cancer cell proliferation in vitro and tumor growth in vivo. Collectively, our study reveals a novel METTL3/m6A/RRBP1 axis in enhancing aggressive traits of PCa, which can be therapeutically targeted by small-peptide METTL3 antagonists.

Melanoma-specific mutation hotspots in distal, non-coding, promoter-interacting regions implicate novel candidate driver genes

BACKGROUND

To develop targeted treatments, it is crucial to identify the full spectrum of genetic drivers in melanoma, including those in non-coding regions. However, recent efforts to explore non-coding regions have primarily focused on gene-adjacent elements such as promoters and non-coding RNAs, leaving intergenic distal regulatory elements largely unexplored.

METHODS

We used Hi-C chromatin contact data from melanoma cells to map distal, non-coding, promoter-interacting regulatory elements genome-wide in melanoma. Using this "promoter-interaction network", alongside whole-genome sequence and gene expression data from the Pan Cancer Analysis of Whole Genomes, we developed multivariate linear regression models to identify distal somatic mutation hotspots that affect promoter activity.

RESULTS

We identified eight recurrently mutated hotspots that are novel, melanoma-specific, located in promoter-interacting distal regulatory elements, alter transcription factor binding motifs, and affect the expression of genes (e.g., HSPB7, CLDN1, ADCY9 and FDXR) previously implicated as tumour suppressors/oncogenes in various cancers.

CONCLUSIONS

Our study suggests additional non-coding drivers beyond the well-characterised TERT promoter in melanoma, offering new insights into the disruption of complex regulatory networks by non-coding mutations that may contribute to melanoma development. Furthermore, our study provides a framework for integrating multiple levels of biological data to uncover cancer-specific non-coding drivers.

KDM6A regulates immune response genes in multiple myeloma

ABSTRACT

The histone H3 at lysine 27 (H3K27) demethylase lysine demethylase 6A (KDM6A) is a tumor suppressor in multiple cancers, including multiple myeloma (MM). We created isogenic MM cells disrupted for KDM6A and tagged the endogenous protein to facilitate genome-wide studies. KDM6A binds genes associated with immune recognition and cytokine signaling. Most importantly, KDM6A binds and activates NLRC5 and CIITA, which encode regulators of major histocompatibility complex genes. Patient data indicate that NLRC5 and CIITA are downregulated in MM with low KDM6A expression. Chromatin analysis shows that KDM6A binds poised and active enhancers and KDM6A loss led to decreased H3K27ac at enhancers, increased H3K27me3 levels in body of genes bound by KDM6A, and decreased gene expression. Reestablishing histone acetylation with an HDAC3 inhibitor leads to upregulation of major histocompatibility complex expression, offering a strategy to restore immunogenicity of KDM6A-deficient tumors. Loss of Kdm6a in Kirsten rat sarcoma virus (K-RAS)-transformed murine fibroblasts led to increased growth in vivo associated with decreased T-cell infiltration.

ARID1A-BAF coordinates ZIC2 genomic occupancy for epithelial-to-mesenchymal transition in cranial neural crest specification

The BAF chromatin remodeler regulates lineage commitment including cranial neural crest cell (CNCC) specification. Variants in BAF subunits cause Coffin-Siris syndrome (CSS), a congenital disorder characterized by coarse craniofacial features and intellectual disability. Approximately 50% of individuals with CSS harbor variants in one of the mutually exclusive BAF subunits, ARID1A/ARID1B. While Arid1a deletion in mouse neural crest causes severe craniofacial phenotypes, little is known about the role of ARID1A in CNCC specification. Using CSS-patient-derived ARID1A+/- induced pluripotent stem cells to model CNCC specification, we discovered that ARID1A-haploinsufficiency impairs epithelial-to-mesenchymal transition (EMT), a process necessary for CNCC delamination and migration from the neural tube. Furthermore, wild-type ARID1A-BAF regulates enhancers associated with EMT genes. ARID1A-BAF binding at these enhancers is impaired in heterozygotes while binding at promoters is unaffected. At the sequence level, these EMT enhancers contain binding motifs for ZIC2, and ZIC2 binding at these sites is ARID1A-dependent. When excluded from EMT enhancers, ZIC2 relocates to neuronal enhancers, triggering aberrant neuronal gene activation. In mice, deletion of Zic2 impairs NCC delamination, while ZIC2 overexpression in chick embryos at post-migratory neural crest stages elicits ectopic delamination from the neural tube. These findings reveal an essential ARID1A-ZIC2 axis essential for EMT and CNCC delamination.

Congenital microcoria deletion in mouse links Sox21 dysregulation to disease and suggests a role for TGFB2 in glaucoma and myopia

Congenital microcoria (MCOR) is a rare hereditary developmental defect of the iris dilator muscle frequently associated with high axial myopia and high intraocular pressure (IOP) glaucoma. The condition is caused by submicroscopic rearrangements of chromosome 13q32.1. However, the mechanisms underlying the failure of iris development and the origin of associated features remain elusive. Here, we present a 3D architecture model of the 13q32.1 region, demonstrating that MCOR-related deletions consistently disrupt the boundary between two topologically associating domains (TADs). Deleting the critical MCOR-causing region in mice reveals ectopic Sox21 expression precisely aligning with Dct, each located in one of the two neighbor TADs. This observation is consistent with the TADs' boundary alteration and adoption of Dct regulatory elements by the Sox21 promoter. Additionally, we identify Tgfb2 as a target gene of SOX21 and show TGFΒ2 accumulation in the aqueous humor of an MCOR-affected subject. Accumulation of TGFB2 is recognized for its role in glaucoma and potential impact on axial myopia. Our results highlight the importance of SOX21-TGFB2 signaling in iris development and control of eye growth and IOP. Insights from MCOR studies may provide therapeutic avenues for this condition but also for glaucoma and high myopia conditions, affecting millions of people.

Integrated multi-omics analysis of zinc-finger proteins uncovers roles in RNA regulation

RNA interactome studies have revealed that hundreds of zinc-finger proteins (ZFPs) are candidate RNA-binding proteins (RBPs), yet their RNA substrates and functional significance remain largely uncharacterized. Here, we present a systematic multi-omics analysis of the DNA- and RNA-binding targets and regulatory roles of more than 100 ZFPs representing 37 zinc-finger families. We show that multiple ZFPs are previously unknown regulators of RNA splicing, alternative polyadenylation, stability, or translation. The examined ZFPs show widespread sequence-specific RNA binding and preferentially bind proximal to transcription start sites. Additionally, several ZFPs associate with their targets at both the DNA and RNA levels. We highlight ZNF277, a C2H2 ZFP that binds thousands of RNA targets and acts as a multi-functional RBP. We also show that ZNF473 is a DNA/RNA-associated protein that regulates the expression and splicing of cell cycle genes. Our results reveal diverse roles for ZFPs in transcriptional and post-transcriptional gene regulation.

N6-Methyladenosine RNA Modification Regulates Maize Resistance to Maize Chlorotic Mottle Virus Infection

Maize chlorotic mottle virus (MCMV) is one of the main viruses causing significant losses in maize. N6-methyladenosine (m6A) RNA modification has been proven to play important regulatory roles in plant development and stress response. In this study, we found that MCMV infection significantly up-regulated the m6A level in maize, and methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were performed to investigate the distribution of m6A modified peaks and gene expression patterns in MCMV-infected maize plants. The results showed that 1325 differentially methylated genes (DMGs) and 47 differentially methylated and expressed genes (DMEGs) were identified and analyzed. Moreover, the results of virus-induced gene silencing (VIGS) assays showed that ZmECT18 and ZmGST31 were required for MCMV infection, while silencing of ZmMTC, ZmSCI1 or ZmTIP1 significantly promoted MCMV infection in maize. Our findings provided novel insights into the regulatory roles of m6A modification in maize response to MCMV infection.

Transcription factor PATZ1 promotes adipogenesis by controlling promoter regulatory loci of adipogenic factors

White adipose tissue (WAT) is essential for lipid storage and systemic energy homeostasis. Understanding adipocyte formation and stability is key to developing therapies for obesity and metabolic disorders. Through a high-throughput cDNA screen, we identified PATZ1, a POZ/BTB and AT-Hook Containing Zinc Finger 1 protein, as an important adipogenic transcription factor. PATZ1 is expressed in human and mouse adipocyte precursor cells (APCs) and adipocytes. In cellular models, PATZ1 promotes adipogenesis via protein-protein interactions and DNA binding. PATZ1 ablation in mouse adipocytes and APCs leads to a reduced APC pool, decreased fat mass, and hypertrophied adipocytes. ChIP-Seq and RNA-seq analyses show that PATZ1 supports adipogenesis by interacting with transcriptional machinery at the promoter regions of key early adipogenic factors. Mass-spec results show that PATZ1 associates with GTF2I, with GTF2I modulating PATZ1's function during differentiation. These findings underscore PATZ1's regulatory role in adipocyte differentiation and adiposity, offering insights into adipose tissue development.

Integrated multi-omics analysis reveals molecular changes associated with chronic lipid accumulation following contusive spinal cord injury

Functional and pathological recovery after spinal cord injury (SCI) is often incomplete due to the limited regenerative capacity of the central nervous system (CNS), which is further impaired by several mechanisms that sustain tissue damage. Among these, the chronic activation of immune cells can cause a persistent state of local CNS inflammation and damage. However, the mechanisms that sustain this persistent maladaptive immune response in SCI have not been fully clarified yet. In this study, we integrated histological analyses with proteomic, lipidomic, transcriptomic, and epitranscriptomic approaches to study the pathological and molecular alterations that develop in a mouse model of cervical spinal cord hemicontusion. We found significant pathological alterations of the lesion rim with myelin damage and axonal loss that persisted throughout the late chronic phase of SCI. This was coupled by a progressive lipid accumulation in myeloid cells, including resident microglia and infiltrating monocyte-derived macrophages. At tissue level, we found significant changes of proteins indicative of glycolytic, tricarboxylic acid cycle (TCA), and fatty acid metabolic pathways with an accumulation of triacylglycerides with C16:0 fatty acyl chains in chronic SCI. Following transcriptomic, proteomic, and epitranscriptomic studies identified an increase of cholesterol and m6A methylation in lipid-droplet-accumulating myeloid cells as a core feature of chronic SCI. By characterizing the multiple metabolic pathways altered in SCI, our work highlights a key role of lipid metabolism in the chronic response of the immune and central nervous system to damage.

Differential gene expression and chromatin accessibility in the medial prefrontal cortex associated with individual differences in rat behavioral models of opioid use disorder

Opioid use disorder (OUD) is a neuropsychological disease that has a devastating impact on public health. Substantial individual differences in vulnerability exist, the neurobiological substrates of which remain unclear. To address this question, we investigated genome-wide gene transcription (RNA-seq) and chromatin accessibility (ATAC-seq) in the medial prefrontal cortex (mPFC) of male and female rats exhibiting differential vulnerability in behavioral paradigms modeling different phases of OUD: Withdrawal-Induced Anhedonia (WIA), Demand, and Reinstatement. Ingenuity Pathway Analysis (IPA) of RNA-seq revealed greater changes in canonical pathways in Resilient (vs. Saline) rats in comparison to Vulnerable (vs. Saline) rats across 3 paradigms, suggesting brain adaptations that might contribute to resilience to OUD across its trajectory. Analyses of gene networks and upstream regulators implicated processes involved in oligodendrocyte maturation and myelination in WIA, neuroinflammation in Demand, and metabolism in Reinstatement. Motif analysis of ATAC-seq showed changes in chromatin accessibility to a small set of transcription factor (TF) binding sites as a function either of opioid exposure (i.e., morphine versus saline) generally or of individual vulnerability specifically. Some of these were shared across the 3 paradigms and others were unique to each. In conclusion, we have identified changes in biological pathways, TFs, and their binding motifs that vary with paradigm and OUD vulnerability. These findings point to the involvement of distinct transcriptional and epigenetic mechanisms in response to opioid exposure, vulnerability to OUD, and different stages of the disorder.

Novel Long Noncoding RNA HEAT4 Affects Monocyte Subtypes, Reducing Inflammation and Promoting Vascular Healing

BACKGROUND

Activation of the immune system contributes to cardiovascular diseases. The role of human-specific long noncoding RNAs in cardioimmunology is poorly understood.

METHODS

Single-cell sequencing in peripheral blood mononuclear cells revealed a novel human-specific long noncoding RNA called HEAT4 (heart failure-associated transcript 4). HEAT4 expression was assessed in several in vitro and ex vivo models of immune cell activation, as well as in the blood of patients with heart failure (HF), acute myocardial infarction, or cardiogenic shock. The transcriptional regulation of HEAT4 was verified through cytokine treatment and single-cell sequencing. Loss-of-function and gain-of-function studies and multiple RNA-protein interaction assays uncovered a mechanistic role of HEAT4 in the monocyte anti-inflammatory gene program. HEAT4 expression and function was characterized in a vascular injury model in NOD.CB17-Prkdc scid/Rj mice.

RESULTS

HEAT4 expression was increased in the blood of patients with HF, acute myocardial infarction, or cardiogenic shock. HEAT4 levels distinguished patients with HF from people without HF and predicted all-cause mortality in a cohort of patients with HF over 7 years of follow-up. Monocytes, particularly anti-inflammatory CD16+ monocytes, which are increased in patients with HF, are the primary source of HEAT4 expression in the blood. HEAT4 is transcriptionally activated by treatment with anti-inflammatory interleukin-10. HEAT4 activates anti-inflammatory and inhibits proinflammatory gene expression. Increased HEAT4 levels result in a shift toward more CD16+ monocytes. HEAT4 binds to S100A9, causing a monocyte subtype switch, thereby reducing inflammation. As a result, HEAT4 improves endothelial barrier integrity during inflammation and promotes vascular healing after injury in mice.

CONCLUSIONS

These results characterize a novel endogenous anti-inflammatory pathway that involves the conversion of monocyte subtypes into anti-inflammatory CD16+ monocytes. The data identify a novel function for the class of long noncoding RNAs by preventing protein secretion and suggest long noncoding RNAs as potential targets for interventions in the field of cardioimmunology.

Endogenous glucocorticoids are required for normal macrophage activation and gastric Helicobacter pylori immunity

Glucocorticoids are steroid hormones well known for their potent anti-inflammatory effects. However, their immunomodulatory properties are multifaceted. Increasing evidence suggests that glucocorticoid signaling promotes effective immunity and that disruption of glucocorticoid signaling impairs immune function. In this study, we conditionally deleted the glucocorticoid receptor (GR) in the myeloid lineage using the LysM-Cre driver (myGRKO). We examined the impact on macrophage activation and gastric immune responses to Helicobacter pylori, the best-known risk factor of gastric cancer. Our results indicate that, compared with wild type (WT), glucocorticoid receptor knockout (GRKO) macrophages exhibited higher expression of proinflammatory genes in steroid-free conditions. However, when challenged in vivo, GRKO macrophages exhibited aberrant chromatin landscapes and impaired proinflammatory gene expression profiles. Moreover, gastric colonization with H. pylori revealed impaired gastric immune responses and reduced T cell recruitment in myGRKO mice. As a result, myGRKO mice were protected from atrophic gastritis and pyloric metaplasia development. These results demonstrate a dual role for glucocorticoid signaling in preparing macrophages to respond to bacterial infection but limiting their pathogenic activation. In addition, our results support that macrophages are critical for gastric H. pylori immunity.NEW & NOTEWORTHY Signaling by endogenous glucocorticoids primes macrophages toward more robust responses to pathogens. Disruption of glucocorticoid signaling caused dysregulation of the chromatin landscape, blunted proinflammatory gene activation upon bacterial challenge, and impaired the gastric inflammatory response to Helicobacter pylori infection.

Effector MoSDT1 enhances Magnaporthe oryzae virulence and plays a dual role in regulating rice defense

C2H2 zinc effectors are a class of pathogen proteins that play a dual role in plant-pathogen interactions, promoting pathogenicity and enhancing plant defense. In our previous research, we identified Magnaporthe oryzae Systemic Defense Trigger 1 (MoSDT1) as a C2H2 zinc effector that activates rice (Oryza sativa) defense when overexpressed in rice. However, its regulatory roles in pathogenicity and defense require further investigation. In this study, we generated an MoSDT1 overexpressing strain and 2 knockout strains of M. oryzae to assess the impact of MoSDT1 on pathogenicity, rice defense, and phenotypic characteristics. Our analyses revealed that MoSDT1 substantially influenced vegetative growth, conidia size, and conidiation, and was crucial for the virulence of M. oryzae while suppressing rice defense. MoSDT1 localized to the nucleus and cytoplasm of rice, either dependent or independent of M. oryzae delivery. Through RNA-seq, scRNA-seq, and ChIP-seq, we identified that MoSDT1 modulates rice defense by regulating the phosphorylation and ubiquitination of various rice signaling proteins, including transcription factors, transcription repressors, kinases, phosphatases, and the ubiquitin system. These findings provide valuable insights into the regulatory mechanisms of C2H2 zinc finger effector proteins and offer important foundational information for utilizing their target genes in disease resistance breeding and the design of targets for disease management.

Suppression of GATA3 promotes epithelial-mesenchymal transition and simultaneous cellular senescence in human extravillous trophoblasts

The regulatory mechanism of the transcription factor GATA3 in the differentiation and maturation process of extravillous trophoblasts (EVT) in early pregnancy placenta, as well as its relevance to the occurrence of pregnancy disorders, remains poorly understood. This study leveraged single-cell RNA sequencing data from placental organoid models and placental tissue to explore the dynamic changes in GATA3 expression during EVT maturation. The expression pattern exhibited an initial upregulation followed by subsequent downregulation, with aberrant GATA3 localization observed in cases of recurrent miscarriage (RM). By identifying global targets regulated by GATA3 in primary placental EVT cells, JEG3, and HTR8/SVneo cell lines, this study offered insights into its regulatory mechanisms across different EVT cell models. Shared regulatory targets among these cell types and activation of trophoblast cell marker genes emphasized the importance of GATA3 in EVT differentiation and maturation. Knockdown of GATA3 in JEG3 cells led to repression of GATA3-induced epithelial-mesenchymal transition (EMT), as evidenced by changes in marker gene expression levels and enhanced migration ability. Additionally, interference with GATA3 accelerated cellular senescence, as indicated by reduced proliferation rates and increased activity levels for senescence-associated β-galactosidase enzyme, along with elevated expression levels for senescence-associated genes. This study provides comprehensive insights into the dual role of GATA3 in regulating EMT and cellular senescence during EVT differentiation, shedding light on the dynamic changes in GATA3 expression in normal and pathological placental conditions.

Integrative metabolomics-genomics analysis identifies key networks in a stem cell-based model of schizophrenia

Schizophrenia (SCZ) is a neuropsychiatric disorder, caused by a combination of genetic and environmental factors. The etiology behind the disorder remains elusive although it is hypothesized to be associated with the aberrant response to neurotransmitters, such as dopamine and glutamate. Therefore, investigating the link between dysregulated metabolites and distorted neurodevelopment holds promise to offer valuable insights into the underlying mechanism of this complex disorder. In this study, we aimed to explore a presumed correlation between the transcriptome and the metabolome in a SCZ model based on patient-derived induced pluripotent stem cells (iPSCs). For this, iPSCs were differentiated towards cortical neurons and samples were collected longitudinally at various developmental stages, reflecting neuroepithelial-like cells, radial glia, young and mature neurons. The samples were analyzed by both RNA-sequencing and targeted metabolomics and the two modalities were used to construct integrative networks in silico. This multi-omics analysis revealed significant perturbations in the polyamine and gamma-aminobutyric acid (GABA) biosynthetic pathways during rosette maturation in SCZ lines. We particularly observed the downregulation of the glutamate decarboxylase encoding genes GAD1 and GAD2, as well as their protein product GAD65/67 and their biochemical product GABA in SCZ samples. Inhibition of ornithine decarboxylase resulted in further decrease of GABA levels suggesting a compensatory activation of the ornithine/putrescine pathway as an alternative route for GABA production. These findings indicate an imbalance of cortical excitatory/inhibitory dynamics occurring during early neurodevelopmental stages in SCZ. Our study supports the hypothesis of disruption of inhibitory circuits to be causative for SCZ and establishes a novel in silico approach that enables for integrative correlation of metabolic and transcriptomic data of psychiatric disease models.

CDK9 inhibition inhibits multiple oncogenic transcriptional and epigenetic pathways in prostate cancer

BACKGROUND

Cyclin-dependent kinase 9 (CDK9) stimulates oncogenic transcriptional pathways in cancer and CDK9 inhibitors have emerged as promising therapeutic candidates.

METHODS

The activity of an orally bioavailable CDK9 inhibitor, CDKI-73, was evaluated in prostate cancer cell lines, a xenograft mouse model, and patient-derived tumor explants and organoids. Expression of CDK9 was evaluated in clinical specimens by mining public datasets and immunohistochemistry. Effects of CDKI-73 on prostate cancer cells were determined by cell-based assays, molecular profiling and transcriptomic/epigenomic approaches.

RESULTS

CDKI-73 inhibited proliferation and enhanced cell death in diverse in vitro and in vivo models of androgen receptor (AR)-driven and AR-independent models. Mechanistically, CDKI-73-mediated inhibition of RNA polymerase II serine 2 phosphorylation resulted in reduced expression of BCL-2 anti-apoptotic factors and transcriptional defects. Transcriptomic and epigenomic approaches revealed that CDKI-73 suppressed signaling pathways regulated by AR, MYC, and BRD4, key drivers of dysregulated transcription in prostate cancer, and reprogrammed cancer-associated super-enhancers. These latter findings prompted the evaluation of CDKI-73 with the BRD4 inhibitor AZD5153, a combination that was synergistic in patient-derived organoids and in vivo.

CONCLUSION

Our work demonstrates that CDK9 inhibition disrupts multiple oncogenic pathways and positions CDKI-73 as a promising therapeutic agent for prostate cancer, particularly aggressive, therapy-resistant subtypes.

IRF8 defines the epigenetic landscape in postnatal microglia, thereby directing their transcriptome programs

Microglia are innate immune cells in the brain. Transcription factor IRF8 (interferon regulatory factor 8) is highly expressed in microglia. However, its role in postnatal microglia development is unknown. We demonstrate that IRF8 binds stepwise to enhancer regions of postnatal microglia along with Sall1 and PU.1, reaching a maximum after day 14. IRF8 binding correlated with a stepwise increase in chromatin accessibility, which preceded the initiation of microglia-specific transcriptome. Constitutive and postnatal Irf8 deletion led to a loss of microglia identity and gain of disease-associated microglia (DAM)-like genes. Combined analysis of single-cell (sc)RNA sequencing and single-cell transposase-accessible chromatin with sequencing (scATAC-seq) revealed a correlation between chromatin accessibility and transcriptome at a single-cell level. IRF8 was also required for microglia-specific DNA methylation patterns. Last, in the 5xFAD model, constitutive and postnatal Irf8 deletion reduced the interaction of microglia with amyloidβ plaques and the size of plaques, lessening neuronal loss. Together, IRF8 sets the epigenetic landscape, which is required for postnatal microglia gene expression.

Gene expression and chromatin conformation of microglia in virally suppressed people with HIV

The presence of HIV in sequestered reservoirs is a central impediment to a functional cure, allowing HIV to persist despite life-long antiretroviral therapy (ART), and driving a variety of comorbid conditions. Our understanding of the latent HIV reservoir in the central nervous system is incomplete, because of difficulties in accessing human central nervous system tissues. Microglia contribute to HIV reservoirs, but the molecular phenotype of HIV-infected microglia is poorly understood. We leveraged the unique "Last Gift" rapid autopsy program, in which people with HIV are closely followed until days or even hours before death. Microglial populations were heterogeneous regarding their gene expression profiles but showed similar chromatin accessibility landscapes. Despite ART, we detected occasional microglia containing cell-associated HIV RNA and HIV DNA integrated into open regions of the host's genome (∼0.005%). Microglia with detectable HIV RNA showed an inflammatory phenotype. These results demonstrate a distinct myeloid cell reservoir in the brains of people with HIV despite suppressive ART. Strategies for curing HIV and neurocognitive impairment will need to consider the myeloid compartment to be successful.

VIP1 and its close homologs confer mechanical stress tolerance in Arabidopsis leaves

VIP1, an Arabidopsis thaliana basic leucine zipper transcription factor, and its close homologs are imported from the cytoplasm to the nucleus when cells are exposed to mechanical stress. They bind to AGCTG (G/T) and regulate mechanical stress responses in roots. However, their role in leaves is unclear. To clarify this, mutant lines (QM1 and QM2) that lack the functions of VIP1 and its close homologs (bZIP29, bZIP30 and PosF21) were generated. Brushing more severely damaged QM1 and QM2 leaves than wild-type leaves. Genes regulating stress responses and cell wall properties were downregulated in brushed QM2 leaves and upregulated in brushed VIP1-GFP-overexpressing (VIP1-GFPox) leaves compared to wild-type leaves in a transcriptome analysis. The VIP1-binding sequence AGCTG (G/T) was enriched in the promoters of genes downregulated in brushed QM2 leaves compared to wild-type leaves and in those upregulated in brushed VIP1-GFPox leaves. Calmodulin-binding transcription activators (CAMTAs) are known regulators of mechanical stress responses, and the CAMTA-binding sequence CGCGT was enriched in the promoters of genes upregulated in the brushed QM2 leaves and in those downregulated in the brushed VIP1-GFPox leaves. These findings suggest that VIP1 and its homologs upregulate genes via AGCTG (G/T) and influence CAMTA-dependent gene expression to enhance mechanical stress tolerance in leaves.

Overexpression of a pearl millet WRKY transcription factor gene, PgWRKY74, in Arabidopsis retards shoot growth under dehydration and salinity-stressed conditions

Pearl millet (Cenchrus americanus) is a cereal crop that can tolerate high temperatures, drought, and low-fertility conditions where other crops lose productivity. However, genes regulating this ability are largely unknown. Transcription factors (TFs) regulate transcription of their target genes, regulate downstream biological processes, and thus are candidates for regulators of such tolerance of pearl millet. PgWRKY74 encodes a group IIc WRKY TF in pearl millet and is downregulated by drought. PgWRKY74 may have a role in drought tolerance. The objective of this study was to gain insights into the physiological and biochemical functions of PgWRKY74. Yeast one-hybrid and gel shift assays were performed to examine transcriptional activation potential and deoxyribonucleic acid (DNA)-binding ability, respectively. Transgenic Arabidopsis thaliana plants overexpressing PgWRKY74-green fluorescent protein (GFP) fusion gene were generated and tested for growth and stress-responsive gene expression under mannitol and NaCl-stressed conditions. A construct with PgWRKY74 enabled yeast reporter cells to survive on test media in the yeast one-hybrid assays. The electrophoretic mobility of DNA with putative WRKY TF-binding motifs was lower in the presence of a recombinant PgWRKY74 protein than its absence. The PgWRKY74-GFP-overexpressing Arabidopsis plants exhibited smaller rosette areas than did wild-type plants under mannitol-stressed and NaCl-stressed conditions, and exhibited weaker expression of RD29B, which is induced by the stress-related phytohormone abscisic acid (ABA), under the mannitol-stressed condition. PgWRKY74 have transcriptional activation potential and DNA-binding ability, and can negatively regulate plant responses to mannitol and NaCl stresses, possibly by decreasing ABA levels or ABA sensitivity.

BootCellNet, a resampling-based procedure, promotes unsupervised identification of cell populations via robust inference of gene regulatory networks

Recent advances in measurement technologies, particularly single-cell RNA sequencing (scRNA-seq), have revolutionized our ability to acquire large amounts of omics-level data on cellular states. As measurement techniques evolve, there has been an increasing need for data analysis methodologies, especially those focused on cell-type identification and inference of gene regulatory networks (GRNs). We have developed a new method named BootCellNet, which employs smoothing and resampling to infer GRNs. Using the inferred GRNs, BootCellNet further infers the minimum dominating set (MDS), a set of genes that determines the dynamics of the entire network. We have demonstrated that BootCellNet robustly infers GRNs and their MDSs from scRNA-seq data and facilitates unsupervised identification of cell clusters using scRNA-seq datasets of peripheral blood mononuclear cells and hematopoiesis. It has also identified COVID-19 patient-specific cells and their potential regulatory transcription factors. BootCellNet not only identifies cell types in an unsupervised and explainable way but also provides insights into the characteristics of identified cell types through the inference of GRNs and MDS.

CEP112 coordinates translational regulation of essential fertility genes during spermiogenesis through phase separation in humans and mice

Spermiogenesis, the complex transformation of haploid spermatids into mature spermatozoa, relies on precise spatiotemporal regulation of gene expression at the post-transcriptional level. The mechanisms underlying this critical process remain incompletely understood. Here, we identify centrosomal protein 112 (CEP112) as an essential regulator of mRNA translation during this critical developmental process. Mutations in CEP112 are discovered in oligoasthenoteratospermic patients, and Cep112-deficient male mice recapitulate key phenotypes of human asthenoteratozoospermia. CEP112 localizes to the neck and atypical centrioles of mature sperm and forms RNA granules during spermiogenesis, enriching target mRNAs such as Fsip2, Cfap61, and Cfap74. Through multi-omics analyses and the TRICK reporter assay, we demonstrate that CEP112 orchestrates the translation of target mRNAs. Co-immunoprecipitation and mass spectrometry identify CEP112's interactions with translation-related proteins, including hnRNPA2B1, EEF1A1, and EIF4A1. In vitro, CEP112 undergoes liquid-liquid phase separation, forming condensates that recruit essential proteins and mRNAs. Moreover, variants in patient-derived CEP112 disrupt phase separation and impair translation efficiency. Our results suggest that CEP112 mediates the assembly of RNA granules through liquid-liquid phase separation to control the post-transcriptional expression of fertility-related genes. This study not only clarifies CEP112's role in spermatogenesis but also highlights the role of phase separation in translational regulation, providing insights into male infertility and suggesting potential therapeutic targets.

Cpt1a Drives primed-to-naïve pluripotency transition through lipid remodeling

Metabolism has been implicated in cell fate determination, particularly through epigenetic modifications. Similarly, lipid remodeling also plays a role in regulating cell fate. Here, we present comprehensive lipidomics analysis during BMP4-driven primed to naive pluripotency transition or BiPNT and demonstrate that lipid remodeling plays an essential role. We further identify Cpt1a as a rate-limiting factor in BiPNT, driving lipid remodeling and metabolic reprogramming while simultaneously increasing intracellular acetyl-CoA levels and enhancing H3K27ac at chromatin open sites. Perturbation of BiPNT by histone acetylation inhibitors suppresses lipid remodeling and pluripotency transition. Together, our study suggests that lipid remodeling promotes pluripotency transitions and further regulates cell fate decisions, implicating Cpt1a as a critical regulator between primed-naive cell fate control.

RNA polymerase II transcription initiation in holo-TFIID-depleted mouse embryonic stem cells

The recognition of core promoter sequences by TFIID is the first step in RNA polymerase II (Pol II) transcription initiation. Metazoan holo-TFIID is a trilobular complex, composed of the TATA binding protein (TBP) and 13 TBP-associated factors (TAFs). Why and how TAFs are necessary for the formation of TFIID domains and how they contribute to transcription initiation remain unclear. Inducible TAF7 or TAF10 depletion, followed by comprehensive analysis of TFIID subcomplex formation, chromatin binding, and nascent transcription in mouse embryonic stem cells, result in the formation of a TAF7-lacking TFIID or a minimal core-TFIID complex, respectively. These partial complexes support TBP recruitment at promoters and nascent Pol II transcription at most genes early after depletion, but importantly, TAF10 is necessary for efficient Pol II pausing. We show that partially assembled TFIID complexes can sustain Pol II transcription initiation but cannot replace holo-TFIID over several cell divisions and/or development.

Hypoxic stabilization of RIPOR3 mRNA via METTL3-mediated m6A methylation drives breast cancer progression and metastasis

Dysregulated N6-methyladenosine (m6A) modification has been associated with breast cancer pathogenesis. Hypoxia which characterizes solid tumors is known to reprogram the m6A epitranscriptome, but the underlying mechanisms of how this process contributes to breast cancer progression remain poorly understood. Through integrative analyses of m6A-RIP sequencing and RNA sequencing databases, we reveal a cluster of mRNAs with upregulated m6A methylation and expression under hypoxia, that are enriched by many oncogenic pathways, including PI3K-Akt signaling. Furthermore, we identify the mRNA, RIPOR3, as a target of METTL3-mediated m6A methylation in response to hypoxia. We find that m6A methylation stabilizes RIPOR3, increasing its protein expression in a METTL3 catalytic activity-dependent manner, and consequently driving breast tumor growth and metastasis. RIPOR3 is found to be overexpressed in breast cancer cell lines and tumor tissues from breast cancer patients, in whom elevated RIPOR3 is associated with a worse prognosis. Mechanistically, we show that RIPOR3 interacts with EGFR and is essential for the PI3K-Akt pathway activation. In conclusion, we identify RIPOR3 as a hypoxia-stabilized oncogenic driver via METTL3-mediated m6A methylation, thus provide a potential therapeutic target for breast cancer.

Genome structural dynamics: insights from Gaussian network analysis of Hi-C data

Characterization of the spatiotemporal properties of the chromatin is essential to gaining insights into the physical bases of gene co-expression, transcriptional regulation and epigenetic modifications. The Gaussian network model (GNM) has proven in recent work to serve as a useful tool for modeling chromatin structural dynamics, using as input high-throughput chromosome conformation capture data. We focus here on the exploration of the collective dynamics of chromosomal structures at hierarchical levels of resolution, from single gene loci to topologically associating domains or entire chromosomes. The GNM permits us to identify long-range interactions between gene loci, shedding light on the role of cross-correlations between distal regions of the chromosomes in regulating gene expression. Notably, GNM analysis performed across diverse cell lines highlights the conservation of the global/cooperative movements of the chromatin across different types of cells. Variations driven by localized couplings between genomic loci, on the other hand, underlie cell differentiation, underscoring the significance of the four-dimensional properties of the genome in defining cellular identity. Finally, we demonstrate the close relation between the cell type-dependent mobility profiles of gene loci and their gene expression patterns, providing a clear demonstration of the role of chromosomal 4D features in defining cell-specific differential expression of genes.

A systematic analyses of different bioinformatics pipelines for genomic data and its impact on deep learning models for chromatin loop prediction

Genomic data analysis has witnessed a surge in complexity and volume, primarily driven by the advent of high-throughput technologies. In particular, studying chromatin loops and structures has become pivotal in understanding gene regulation and genome organization. This systematic investigation explores the realm of specialized bioinformatics pipelines designed specifically for the analysis of chromatin loops and structures. Our investigation incorporates two protein (CTCF and Cohesin) factor-specific loop interaction datasets from six distinct pipelines, amassing a comprehensive collection of 36 diverse datasets. Through a meticulous review of existing literature, we offer a holistic perspective on the methodologies, tools and algorithms underpinning the analysis of this multifaceted genomic feature. We illuminate the vast array of approaches deployed, encompassing pivotal aspects such as data preparation pipeline, preprocessing, statistical features and modelling techniques. Beyond this, we rigorously assess the strengths and limitations inherent in these bioinformatics pipelines, shedding light on the interplay between data quality and the performance of deep learning models, ultimately advancing our comprehension of genomic intricacies.

Understanding m6A changes in chromophobe renal cell carcinoma and predicting patient outcomes survival

N6-methyladenosine (m6A) is a prevalent mRNA modification known for its implications in various cancer types, yet its role in chromophobe renal cell carcinoma (chRCC) remains largely unexplored. In this study, we performed m6A-SEAL-seq and RNA-seq analyses on tissues from three chRCC subjects, aiming to uncover m6A alterations in chRCC. Our findings revealed reduced expression levels of four m6A regulators in chRCC tissues and highlighted differences in m6A levels compared to normal tissues. Furthermore, we identified specific genes and cancer-related pathways affected by these differences, including notable candidates like NOTCH1 and FGFR1, implicated in chRCC development. Additionally, we developed a predictive model based on the expression level of m6A associated genes, demonstrating promising prognostic capabilities for patient survival prediction. Overall, our study provides valuable insights into the role of m6A in chRCC and its potential as a prognostic indicator.

Relevance of the regulation of the brain-placental axis to the nocturnal bottleneck of mammals

INTRODUCTION

Evolutionary theory suggests that the ancestors of all placental animals were nocturnal. Visual perceptive function of mammalian brain has evolved extensively, but nearly 70 % of today's mammals are still nocturnal. While placental influence on brain development is known, if placenta plays a role in the visual perceptive function of mammalian brain remains untested. The present study aims to test this hypothesis.

METHODS

In this study, single-nuclei RNA sequencing was performed to identify genes expressed in the pig placenta and fetal brain, and then compared with the orthologous genes expressed in the placenta and fetal brain cells of mouse. Differential gene expression analysis was performed to identify placental genes regulated differentially between nocturnal and diurnal animals. Phylogenetic modeling was performed to test correlated evolution between placenta type, and the nocturnal or diurnal activity among different mammals.

RESULTS

The results showed that genes differentially regulated in the fetal brain were related to visual perception whereas the placental genes were related to the nocturnal or diurnal activity in placental animals. Phylogenetic modeling of these genes in thirty-four diverse mammalian species showed evidence for evolutionary link between placenta and the nocturnal/diurnal activity in animals.

DISCUSSION

The findings of this study suggest that the placenta plays a role in the evolution of visual perceptive function of brain to shape the nocturnal or diurnal activity of placental animals.

Erasure of DNA methylation in rat fetal germ cells is sex-specific and sensitive to maternal high-fat diet

In mammals, DNA methylation (DNAme) erasure and reinstatement during embryo development and germline establishment are sensitive to the intrauterine environment. Maternal intake of a high-fat diet (HFD), associated with excessive gestational weight gain, has transgenerational effects on offspring health, which may be mediated by changes in DNAme in the germline. Here, we tested the impact of a maternal HFD on embryonic germline DNAme erasure using a rat strain that expresses green fluorescent protein specifically in germ cells. DNAme was analysed by methyl-seq capture in germ cells collected from male and female F1 gonads at gestational day 16. Our data show that although HFD induced global hypomethylation in both sexes, DNAme erasure in female germ cells was more advanced compared to male germ cells. The delay in DNAme erasure in males and the greater impact of HFD suggest that male germ cells are more vulnerable to alterations by exogenous factors.

Zinc-Finger Protein ZFP488 Regulates the Timing of Oligodendrocyte Myelination and Remyelination

Oligodendrocyte myelination and remyelination after injury are intricately regulated by various intrinsic and extrinsic factors, including transcriptional regulators. Among these, the zinc-finger protein ZFP488 is an oligodendrocyte-enriched transcriptional regulator that promotes oligodendrocyte differentiation in the developing neural tube and in oligodendroglial cell lines. However, the specific in vivo genetic requirements for ZFP488 during oligodendrocyte development and remyelination have not been defined. To address this gap, we generated a lineage-traceable ZFP488 knock-out mouse line, wherein an H2b-GFP reporter replaces the ZFP488-coding region. Using these mice of either sex, we examined the dynamics of ZFP488 expression from the endogenous promoter in the developing central nervous system (CNS). We observed a unique expression pattern in the oligodendrocyte lineage, with ZFP488 expression particularly enriched in differentiated oligodendrocytes. ZFP488 loss resulted in delayed myelination in the developing CNS and impaired remyelination after demyelinating injury in the brain. Integrated transcriptomic and genomic profiling further revealed that ZFP488 loss decreased the expression of myelination-associated genes but not oligodendrocyte progenitor-associated genes, suggesting that ZFP488 serves as a positive regulator of myelination by regulating maturation programs. Thus, our genetic loss-of-function study revealed that ZFP488 regulates a stage-dependent differentiation program that controls the timing of CNS myelination and remyelination.

Butyrate and propionate are microbial danger signals that activate the NLRP3 inflammasome in human macrophages upon TLR stimulation

Short-chain fatty acids (SCFAs) are immunomodulatory compounds produced by the microbiome through dietary fiber fermentation. Although generally considered beneficial for gut health, patients suffering from inflammatory bowel disease (IBD) display poor tolerance to fiber-rich diets, suggesting that SCFAs may have contrary effects under inflammatory conditions. To investigate this, we examined the effect of SCFAs on human macrophages in the presence of Toll-like receptor (TLR) agonists. In contrast to anti-inflammatory effects under steady-state conditions, we found that butyrate and propionate activated the NOD-, LRR-, and pyrin domain-containing protein 3 (NLRP3) inflammasome in the presence of TLR agonists. Mechanistically, these SCFAs prevented transcription of FLICE-like inhibitory protein (cFLIP) and interleukin-10 (IL-10) through histone deacetylase (HDAC) inhibition, triggering caspase-8-dependent NLRP3 inflammasome activation. SCFA-driven NLRP3 activation was potassium efflux independent and did not result in cell death but rather triggered hyperactivation and IL-1β release. Our findings demonstrate that butyrate and propionate are bacterially derived danger signals that regulate NLRP3 inflammasome activation through epigenetic modulation of the inflammatory response.

A Comparative Study of Algorithms Detecting Differential Rhythmicity in Transcriptomic Data

Rhythmic transcripts play pivotal roles in driving the daily oscillations of various biological processes. Genetic or environmental disruptions can lead to alterations in the rhythmicity of transcripts, ultimately impacting downstream circadian outputs, including metabolic processes and even behavior. To statistically compare the differences in transcript rhythms between 2 or more conditions, several algorithms have been developed to analyze circadian transcriptomic data, each with distinct features. In this study, we compared the performance of 7 algorithms that were specifically designed to detect differential rhythmicity (DODR, LimoRhyde, CircaCompare, compareRhythms, diffCircadian, dryR, and RepeatedCircadian). We found that even when applying the same statistical threshold, these algorithms yielded varying numbers of differentially rhythmic transcripts, most likely because each algorithm defines rhythmic and differentially rhythmic transcripts differently. Nevertheless, the output for the differential phase and amplitude were identical between dryR and compareRhyhms, and diffCircadian and CircaCompare, while the output from LimoRhyde2 was highly correlated with that from diffCircadian and CircaCompare. Because each algorithm has unique requirements for input data and reports different information as an output, it is crucial to ensure the compatibility of input data with the chosen algorithm and assess whether the algorithm's output fits the user's needs when selecting an algorithm for analysis.

Large-scale multi-omic analysis identifies noncoding somatic driver mutations and nominates ZFP36L2 as a driver gene for pancreatic ductal adenocarcinoma

Identification of somatic driver mutations in the noncoding genome remains challenging. To comprehensively characterize noncoding driver mutations for pancreatic ductal adenocarcinoma (PDAC), we first created genome-scale maps of accessible chromatin regions (ACRs) and histone modification marks (HMMs) in pancreatic cell lines and purified pancreatic acinar and duct cells. Integration with whole-genome mutation calls from 506 PDACs revealed 314 ACRs/HMMs significantly enriched with 3,614 noncoding somatic mutations (NCSMs). Functional assessment using massively parallel reporter assays (MPRA) identified 178 NCSMs impacting reporter activity (19.45% of those tested). Focused luciferase validation confirmed negative effects on gene regulatory activity for NCSMs near CDKN2A and ZFP36L2 . For the latter, CRISPR interference (CRISPRi) further identified ZFP36L2 as a target gene (16.0 - 24.0% reduced expression, P = 0.023-0.0047) with disrupted KLF9 binding likely mediating the effect. Our integrative approach provides a catalog of potentially functional noncoding driver mutations and nominates ZFP36L2 as a PDAC driver gene.

The androgen receptor in mesenchymal progenitors regulates skeletal muscle mass via Igf1 expression in male mice

Androgens exert their effects primarily by binding to the androgen receptor (AR), a ligand-dependent nuclear receptor. While androgens have anabolic effects on skeletal muscle, previous studies reported that AR functions in myofibers to regulate skeletal muscle quality, rather than skeletal muscle mass. Therefore, the anabolic effects of androgens are exerted via nonmyofiber cells. In this context, the cellular and molecular mechanisms of AR in mesenchymal progenitors, which play a crucial role in maintaining skeletal muscle homeostasis, remain largely unknown. In this study, we demonstrated expression of AR in mesenchymal progenitors and found that targeted AR ablation in mesenchymal progenitors reduced limb muscle mass in mature adult, but not young or aged, male mice, although fatty infiltration of muscle was not affected. The absence of AR in mesenchymal progenitors led to remarkable perineal muscle hypotrophy, regardless of age, due to abnormal regulation of transcripts associated with cell death and extracellular matrix organization. Additionally, we revealed that AR in mesenchymal progenitors regulates the expression of insulin-like growth factor 1 (Igf1) and that IGF1 administration prevents perineal muscle atrophy in a paracrine manner. These findings indicate that the anabolic effects of androgens regulate skeletal muscle mass via, at least in part, AR signaling in mesenchymal progenitors.

Systematic discovery of DNA-binding tandem repeat proteins

Tandem repeat proteins (TRPs) are widely distributed and bind to a wide variety of ligands. DNA-binding TRPs such as zinc finger (ZNF) and transcription activator-like effector (TALE) play important roles in biology and biotechnology. In this study, we first conducted an extensive analysis of TRPs in public databases, and found that the enormous diversity of TRPs is largely unexplored. We then focused our efforts on identifying novel TRPs possessing DNA-binding capabilities. We established a protein language model for DNA-binding protein prediction (PLM-DBPPred), and predicted a large number of DNA-binding TRPs. A subset was then selected for experimental screening, leading to the identification of 11 novel DNA-binding TRPs, with six showing sequence specificity. Notably, members of the STAR (Short TALE-like Repeat proteins) family can be programmed to target specific 9 bp DNA sequences with high affinity. Leveraging this property, we generated artificial transcription factors using reprogrammed STAR proteins and achieved targeted activation of endogenous gene sets. Furthermore, the members of novel families such as MOON (Marine Organism-Originated DNA binding protein) and pTERF (prokaryotic mTERF-like protein) exhibit unique features and distinct DNA-binding characteristics, revealing interesting biological clues. Our study expands the diversity of DNA-binding TRPs, and demonstrates that a systematic approach greatly enhances the discovery of new biological insights and tools.

Meiotic chromatin-associated HSF5 is indispensable for pachynema progression and male fertility

Pachynema progression contributes to the completion of prophase I. Nevertheless, the regulation of this significant meiotic process remains poorly understood. In this study, we identified a novel testis-specific protein HSF5, which regulates pachynema progression during male meiosis in a manner dependent on chromatin-binding. Deficiency of HSF5 results in meiotic arrest and male infertility, characterized as unconventional pachynema arrested at the mid-to-late stage, with extensive spermatocyte apoptosis. Our scRNA-seq data confirmed consistent expressional alterations of certain driver genes (Sycp1, Msh4, Meiob, etc.) crucial for pachynema progression in Hsf5-/- individuals. HSF5 was revealed to primarily bind to promoter regions of such key divers by CUT&Tag analysis. Also, our results demonstrated that HSF5 biologically interacted with SMARCA5, SMARCA4 and SMARCE1, and it could function as a transcription factor for pachynema progression during meiosis. Therefore, our study underscores the importance of the chromatin-associated HSF5 for the differentiation of spermatocytes, improving the protein regulatory network of the pachynema progression.

Global RNA Interaction and Transcriptome Profiles Demonstrate the Potential Anti-Oncogenic Targets and Pathways of RBM6 in HeLa Cells

BACKGROUND

The fate and functions of RNAs are coordinately regulated by RNA-binding proteins (RBPs), which are often dysregulated in various cancers. Known as a splicing regulator, RNA-binding motif protein 6 (RBM6) harbors tumor-suppressor activity in many cancers; however, there is a lack of research on the molecular targets and regulatory mechanisms of RBM6.

METHODS

In this study, we constructed an RBM6 knock-down (shRBM6) model in the HeLa cell line to investigate its functions and molecular targets. Then we applied improved RNA immunoprecipitation coupled with sequencing (iRIP-seq) and whole transcriptome sequencing approaches to investigate the potential role and RNA targets of RBM6.

RESULTS

Using The Cancer Genome Atlas dataset, we found that higher expression of RBM6 is associated with a better prognosis in many cancer types. In addition, we found that RBM6 knockdown promoted cell proliferation and inhibited apoptosis, demonstrating that RBM6 may act as an anti-oncogenic protein in cancer cells. RBM6 can regulate the alternative splicing (AS) of genes involved in DNA damage response, proliferation, and apoptosis-associated pathways. Meanwhile, RBM6 knockdown activated type I interferon signaling pathways and inhibited the expression of genes involved in the cell cycle, cellular responses to DNA damage, and DNA repair pathways. The differentially expressed genes (DEGs) by shRBM6 and their involved pathways were likely regulated by the transcription factors undergoing aberrant AS by RBM6 knockdown. For iRIP-seq analysis, we found that RBM6 could interact with a large number of mRNAs, with a tendency for binding motifs GGCGAUG and CUCU. RBM6 bound to the mRNA of cell proliferation- and apoptosis-associated genes with dysregulated AS after RBM6 knockdown.

CONCLUSIONS

In summary, our study highlights the important role of RBM6, as well as the downstream targets and regulated pathways, suggesting the potential regulatory mechanisms of RBM6 in the development of cancer.