Database Resources of the National Genomics Data Center, China National Center for Bioinformation in 2022

The methyltransferase SETD3 regulates mRNA alternative splicing through interacting with hnRNPK

The methyltransferase SETD3 is an enzyme essential for catalyzing histidine-73 methylation on β-Actin, thereby promoting its polymerization and regulating muscle contraction. Although increasing evidence suggests that SETD3 is involved in multiple physiological or pathological events, its biological functions remain incompletely understood. In this study, we utilize in situ proximity labeling combined with mass spectrometry analysis to detect potential interacting partners of SETD3. Unexpectedly, we find that many splicing factors are associated with SETD3. Genome-wide RNA sequencing reveals that SETD3 regulates pre-mRNA splicing events, predominantly influencing exon skipping. Biochemical and bioinformatic analyses suggest that SETD3 interacts with hnRNPK, and they collaboratively regulate exon skipping in a common subset of genes. Functionally, we demonstrate that SETD3 and hnRNPK are required for retention of exon 7 skipping in the FNIP1 gene. This promotes FNIP1-mediated nuclear translocation of the transcription factor TFEB and the subsequent induction of lysosomal and mitochondrial biogenesis. Overall, this study uncovers a novel function of SETD3 in modulating mRNA exon splicing.

Revealing the intricate temporal dynamics and adaptive responses of prokaryotic and eukaryotic microbes in the coastal South China Sea

This comprehensive two-year investigation in the coastal South China Sea has advanced our understanding of marine microbes at both community and genomic levels. By combining metagenomics and metatranscriptomics, we have revealed the intricate temporal dynamics and remarkable adaptability of microbial communities and phytoplankton metagenome-assembled genomes (MAGs) in response to environmental fluctuations. We observed distinct seasonal shifts in microbial community composition and function: cyanobacteria were predominant during warmer months, whereas photosynthetic protists were more abundant during colder seasons. Notably, metabolic marker KOs of photosynthesis were consistently active throughout the year, underscoring the persistent role of these processes irrespective of seasonal changes. Our analysis reveals that environmental parameters such as temperature, salinity, and nitrate concentrations profoundly influence microbial community composition, while temperature and silicate have emerged as crucial factors shaping their functional traits. Through the recovery and analysis of 37 phytoplankton MAGs, encompassing nine prokaryotic cyanobacteria and 28 eukaryotic protists from diverse phyla, we have gained insights into their genetic diversity and metabolic capabilities. Distinct profiles of photosynthesis-related pathways including carbon fixation, carotenoid biosynthesis, photosynthesis-antenna proteins, and photosynthesis among the MAGs indicated their genetic adaptations to changing environmental conditions. This study not only enhances our understanding of microbial dynamics in coastal marine ecosystems but also sheds light on the ecological roles and adaptive responses of different microbial groups to environmental changes.

In vitro monitoring of drug resistance emergence during stepwise induction of bedaquiline and clofazimine, alone and in combination: a phenotypic and genotypic analysis

OBJECTIVES

The co-resistance between bedaquiline and clofazimine raises significant concerns, as they are commonly co-administered as core drugs in drug-resistant TB regimens. The present study aimed to monitor drug resistance-associated gene mutations and the phenotypic change in Mycobacterium tuberculosis (Mtb) under a stepwise drug resistance induction in vitro using bedaquiline, clofazimine or combined drugs.

METHODS

Drug-resistant Mtb strains were gradually induced in vitro on a drug-containing solid medium with a 2-fold increasing concentration of bedaquiline, clofazimine and their combination. The MIC of the induced drug-resistant Mtb strains was determined. The drug resistance-associated genes, including Rv0678, Rv1979c, atpE and pepQ, were sequenced and analysed.

RESULTS

Unlike exposure to bedaquiline alone or the combination of these two drugs, clofazimine alone resulted in drug resistance gene mutations occurring later, specifically in the fourth round of induction as opposed to the second round of induction. Besides, nucleotide deletion or insertion in Rv0678 was the main mutation type for induction under the two-drug combination, while single-nucleotide polymorphisms (SNPs) in Rv0678 were the major mutation types when induced by bedaquiline or clofazimine alone. Rv0678 mutation happened at a relatively lower bedaquiline concentration exposure alone, while atpE mutation occurred at a higher bedaquiline concentration. Regardless of the drug exposure manner, a strong correlation between bedaquiline MICs and clofazimine MICs was observed in all drug resistance strains.

CONCLUSIONS

Combined exposure to bedaquiline and clofazimine developed Rv0678 mutation as early as exposure to bedaquiline alone. However, rather than SNPs, deletion and insertion were the dominant mutation types in dual-drug exposure strain.

Database Resources of the National Genomics Data Center, China National Center for Bioinformation in 2025

The National Genomics Data Center (NGDC), which is a part of the China National Center for Bioinformation (CNCB), offers a comprehensive suite of database resources to support the global scientific community. Amidst the unprecedented accumulation of multi-omics data, CNCB-NGDC is committed to continually evolving and updating its core database resources through big data archiving, integrative analysis and value-added curation. Over the past year, CNCB-NGDC has expanded its collaborations with international databases and established new subcenters focusing on biodiversity, traditional Chinese medicine and tumor genetics. Substantial efforts have been made toward encompassing a broad spectrum of multi-omics data, developing innovative resources and enhancing existing resources. Notably, new resources have been developed for single-cell omics (scTWAS Atlas), genome and variation (VDGE), health and disease (CVD Atlas, CPMKG, Immunosenescence Inventory, HemAtlas, Cyclicpepedia, IDeAS), biodiversity and biosynthesis (RefMetaPlant, MASH-Ocean) and research tools (CCLHunter). All resources and services are publicly accessible at https://ngdc.cncb.ac.cn.

PU.1 eviction at lymphocyte-specific chromatin domains mediates glucocorticoid response in acute lymphoblastic leukemia

The epigenetic landscape plays a critical role in cancer progression, yet its therapeutic potential remains underexplored. Glucocorticoids are essential components of treatments for lymphoid cancers, but resistance, driven in part by epigenetic changes at glucocorticoid-response elements, poses a major challenge to effective therapies. Here we show that glucocorticoid treatment induces distinct patterns of chromosomal organization in glucocorticoid-sensitive and resistant acute lymphoblastic leukemia xenograft models. These glucocorticoid-response elements are primed by the pioneer transcription factor PU.1, which interacts with the glucocorticoid receptor. Eviction of PU.1 promotes receptor binding, increasing the expression of genes involved in apoptosis and facilitating a stronger therapeutic response. Treatment with a PU.1 inhibitor enhances glucocorticoid sensitivity, demonstrating the clinical potential of targeting this pathway. This study uncovers a mechanism involving PU.1 and the glucocorticoid receptor, linking transcription factor activity with drug response, and suggesting potential therapeutic strategies for overcoming resistance.

MARK2 variants cause autism spectrum disorder via the downregulation of WNT/β-catenin signaling pathway

Microtubule affinity-regulating kinase 2 (MARK2) contributes to establishing neuronal polarity and developing dendritic spines. Although large-scale sequencing studies have associated MARK2 variants with autism spectrum disorder (ASD), the clinical features and variant spectrum in affected individuals with MARK2 variants, early developmental phenotypes in mutant human neurons, and the pathogenic mechanism underlying effects on neuronal development have remained unclear. Here, we report 31 individuals with MARK2 variants and presenting with ASD, other neurodevelopmental disorders, and distinctive facial features. Loss-of-function (LoF) variants predominate (81%) in affected individuals, while computational analysis and in vitro expression assay of missense variants supported the effect of MARK2 loss. Using proband-derived and CRISPR-engineered isogenic induced pluripotent stem cells (iPSCs), we show that MARK2 loss leads to early neuronal developmental and functional deficits, including anomalous polarity and dis-organization in neural rosettes, as well as imbalanced proliferation and differentiation in neural progenitor cells (NPCs). Mark2+/- mice showed abnormal cortical formation and partition and ASD-like behavior. Through the use of RNA sequencing (RNA-seq) and lithium treatment, we link MARK2 loss to downregulation of the WNT/β-catenin signaling pathway and identify lithium as a potential drug for treating MARK2-associated ASD.

Zinc finger transcription factors BnaSTOP2s regulate sulfur metabolism and confer Sclerotinia sclerotiorum resistance in Brassica napus

Rapeseed (Brassica napus L.) exhibits high-sulfur requirements to achieve optimal growth, development, and pathogen resistance. Despite the importance of sulfur, the mechanisms regulating its metabolism and disease resistance are not fully understood. In this study, we found that the zinc finger transcription factors BnaSTOP2s play a pivotal role in sulfur metabolism and Sclerotinia sclerotiorum resistance. Our findings indicate that BnaSTOP2s are involved in sulfur metabolism, as evidenced by extensive protein interaction screening. BnaSTOP2s knockout reduced the content of essential sulfur-containing metabolites, including glucosinolate and glutathione, which is consistent with the significantly lowered transcriptional levels of BnaMYB28s and BnaGTR2s, key factors involved in glucosinolate synthesis and transportation, respectively. Comprehensive RNA-seq analysis revealed the substantial effect of BnaSTOP2s on sulfur metabolism from roots to siliques, which serve as pivotal sources and sinks for sulfur metabolism, respectively. Furthermore, we found that leaf lesion size significantly decreased and increased in the BnaSTOP2-OE and Bnastop2 mutants, respectively, compared with the wild-type during S. sclerotiorum infection, suggesting a vital role of BnaSTOP2s in plant defense response. In conclusion, BnaSTOP2s act as global regulators of sulfur metabolism and confer resistance to S. sclerotiorum infection in B. napus. Thus, they have potential implications for improving crop resilience.

Potassium ion channel modulation at cancer-neural interface enhances neuronal excitability in epileptogenic glioblastoma multiforme

The central nervous system (CNS) is increasingly recognized as a critical modulator in the oncogenesis of glioblastoma multiforme (GBM), with interactions between cancer and local neuronal circuits frequently leading to epilepsy; however, the relative contributions of these factors remain unclear. Here, we report a coordinated intratumor shift among distinct cancer subtypes within progenitor-like families of epileptic GBM patients, revealing an accumulation of oligodendrocyte progenitor (OPC)-like subpopulations at the cancer-neuron interface along with heightened electrical signaling activity in the surrounding neuronal networks. The OPC-like cells associated with epilepsy express KCND2, which encodes the voltage-gated K+ channel KV4.2, enhancing neuronal excitability via accumulation of extracellular K+, as demonstrated in patient-derived ex vivo slices, xenografting models, and engineering organoids. Together, we uncovered the essential local circuitry, cellular components, and molecular mechanisms facilitating cancer-neuron interaction at peritumor borders. KCND2 plays a crucial role in mediating nervous system-cancer electrical communication, suggesting potential targets for intervention.

Hypoxia adaptation mechanism in rats' peripheral auditory system in high altitude migration: a time series transcriptome analysis

High altitude is characterized by low oxygen, low pressure, and high radiation. When migrates from low to high altitudes, the body's tissues and organs experience hypoxic stress and will present acoustic adaptation as the protective response. However, the mechanisms of acoustic adaptation at high altitudes remain unclear. In this study, cochlear tissues from Wistar rats were collected at 15, 30, 60, 120, and 180 days after high-altitude migration. Transcriptome sequencing was conducted and DESeq algorithm revealed expression patterns of Differentially Expressed Genes(DEGs) after high altitude migration. Day 60 is a critical stage for cochlear tissue "damage" and "repair" in high-altitude conditions. Transmission Electron Microscopy (TEM) observations of structures also support the findings. A time-series gene co-expression network algorithm was used to investigate gene regulatory patterns and key genes after migration. Immunofluorescence, immunohistochemistry, and qPCR were per-formed for key gene validation and localization. At Day 60, the peak DEG count occurs in rats migrating to high altitude, aligning with the critical phase for cochlear tissue damage and repair at high altitudes. Repair hinges on synaptic plasticity and myelination-linked processes, influencing modules M4 to M6. Module M4's activation gradually diminishes from its peak. However, the 'damage' effect is orchestrated by inflammation-related processes in modules M3 to M5, with module M3's activation also waning. Key gene module M4, pivotal for repair during this pivotal phase, encompasses Sptbn5, Cldn1, Gfra2, and Lims2 as its core genes. Immunohistochemistry reveals Sptbn5's presence in cochlear neurons, hair cells, Schwann cells and stria vascularis tissue. Cldn1 and Gfra2 predominantly localize within the cochlear neuron region. These results may suggest new directions for future research on acoustic acclimatization to high altitude.

Targeting HTR2B suppresses nonfunctioning pituitary adenoma growth and sensitizes cabergoline treatment via inhibiting Gαq/PLC/PKCγ/STAT3 axis

BACKGROUND

Managing nonfunctioning pituitary adenomas (NFPAs) is difficult due to limited drug treatments. Cabergoline's (CAB) effectiveness for NFPAs is debated. This study explores the role of HTR2B in NFPAs and its therapeutic potential.

METHODS

We conducted screening of bulk RNA-sequencing data to analyze HTR2B expression levels in NFPA samples. In vitro and in vivo experiments were performed to evaluate the effects of HTR2B modulation on tumor growth and cell cycle regulation. Mechanistic insights into the HTR2B-mediated signaling pathway were elucidated using pharmacological inhibitors and molecular interaction assays.

RESULTS

Elevated HTR2B expression was detected in NFPA samples, which was associated with increased tumor survival. Inhibition of HTR2B activity resulted in the suppression of tumor growth through modulation of the G2M cell cycle. The inhibition of HTR2B with PRX-08066 was found to block STAT3 phosphorylation and nuclear translocation by interfering with the Gαq/PLC/PKC pathway. A direct interaction between PKC-γ and STAT3 was critical for STAT3 activation. CAB was shown to activate pSTAT3 via HTR2B, reducing its therapeutic potential. However, the combination of an HTR2B antagonist with CAB significantly inhibited tumor cell proliferation in HTR2B-expressing pituitary tumor cell lines, a xenografted pituitary tumor model, and patient-derived samples. Analysis of patient-derived data indicated that a distinct molecular pattern characterized by upregulated HTR2B/PKC-γ and downregulated BTG2/GADD45A may benefit from combination treatment with CAB and PRX-08066.

CONCLUSIONS

HTR2B is a potential therapeutic target for NFPAs, and its inhibition could improve CAB efficacy. A dual therapy approach may be beneficial for NFPA patients with high HTR2B expression.

WTAP weakens oxaliplatin chemosensitivity of colorectal cancer by preventing PANoptosis

As the most abundant post-transcriptional modification in eukaryotes, N6-methyladenosine (m6A) plays a crucial role in cancer cell proliferation, invasion and chemoresistance. However, its specific effects on chemosensitivity to oxaliplatin-based regimens and the impact of these drugs on m6A methylation levels in colorectal cancer (CRC) remain largely unexplored. In this study, we demonstrated that the m6A methyltransferase Wilms tumor 1-associating protein (WTAP) weakens oxaliplatin chemosensitivity in HCT116 and DLD1 cells. Mechanistically, oxaliplatin treatment upregulated WTAP expression, preventing multiple forms of cell death simultaneously, a process known as PANoptosis, by decreasing intracellular oxidative stress through maintaining the expression of nuclear factor erythroid-2-related factor 2 (NRF2), a major antioxidant response element, in an m6A-dependent manner. In addition, high WTAP expression in CRC patients is associated with a poor prognosis and reduced benefit from standard chemotherapy by clinical data analysis of The Cancer Genome Atlas (TCGA) database and patient cohort study. These findings suggest that targeting WTAP-NRF2-PANoptosis axis could enhance the antitumor efficacy of oxaliplatin-based chemotherapy in CRC treatment.

The transcription factor EjNAC5 regulates loquat fruit chilling lignification

Changes in both lignin biosynthesis and DNA methylation have been reported to be associated with chilling stress in plants. When stored at low temperatures, red-fleshed loquat is prone to lignification, with increased lignin content and fruit firmness, which has deleterious effects on taste and eating quality. Here, we found that 5 °C storage mitigated the increasing firmness and lignin content of red-fleshed 'Dahongpao' ('DHP') loquat fruit that occurred during 0 °C storage. EjNAC5 was identified by integrating RNA sequencing with whole-genome bisulfite sequencing analysis of 'DHP' loquat fruit. The transcript levels of EjNAC5 were positively correlated with changes in firmness and negatively correlated with changes in DNA methylation level of a differentially methylated region in the EjNAC5 promoter. In white-fleshed 'Baisha' ('BS') loquat fruit, which do not undergo chilling-induced lignification at 0 °C, the transcripts of EjNAC5 remained low and the methylation level of the differentially methylated region in the EjNAC5 promoter was higher, compared with 'DHP' loquat fruit. Transient overexpression of EjNAC5 in loquat fruit and stable overexpression in Arabidopsis and liverwort led to an increase in lignin content. Furthermore, EjNAC5 interacts with EjERF39 and EjHB1 and activates the transcription of Ej4CL1 and EjPRX12 genes involved in lignin biosynthesis. This regulatory network involves different transcription factors from those involved in the lignification pathway. Our study indicates that EjNAC5 promoter methylation modulates EjNAC5 transcript levels and identifies novel EjNAC5-EjERF39-Ej4CL1 and EjNAC5-EjHB1-EjPRX12 regulatory modules involved in chilling induced-lignification.

Case report: Single-cell RNA sequencing of PBMCs highlights monocyte gene expression alterations in a type A HBV-ACLF patient

Hepatitis B Virus-related acute-on-chronic liver failure (HBV-ACLF) is a severe complication with high fatality rates. However, the underlying mechanisms are still elusive and require further investigation. In this report, we described a case of type A HBV-ACLF in which significant changes were found in monocyte gene expression through single-cell RNA sequencing (scRNA-seq). Furthermore, we observed a shifted M1/M2 polarization as well as dynamic changes in HBV-ACLF markers expression within the circulating monocyte population. The co-expression of HBV-ACLF markers (MERTK, THBS1, PPARγ, and SEMA6B) in the circulating monocyte population suggests that monocytes could play an essential role in the development of HBV-ACLF. By analyzing a public HBV-ACLF cohort with bulk RNA-seq data (64 patients), we showed that the expression level of monocytes marker CD163 gradually increased among normal control individuals (NC, n = 15), patients with liver cirrhosis (LC, n = 10), patients with chronic hepatitis B infection (CHB, n = 10), patients with acute-on-chronic hepatic dysfunction (ACHD, n = 10), and patients with HBV-ACLF (n = 20). Furthermore, the representative HBV-ACLF marker THBS1 was significantly correlated with CD163 in this large clinical cohort. It indicated that the dramatic alteration in monocytes may not be limited to our type A HBV-ACLF patient alone but rather a common phenomenon in HBV-ACLF patients. Collectively, our scRNA-seq analysis showed that the pro-inflammatory status of monocytes had shifted into an anti-inflammatory status in this patient, indicating successful treatment and benign prognosis. Although scRNA-seq is still a time-consuming procedure and difficult to apply in daily clinical practice, this report preliminarily shows the promising potential utility of scRNA-seq in HBV-ACLF patients, by which altered status of monocytes could be unbiasedly detected.

Photosynthetic characteristics and genetic mapping of a yellow-green leaf mutant jym165 in soybean

BACKGROUND

Leaves are important sites for photosynthesis and can convert inorganic substances into organic matter. Photosynthetic performance is an important factor affecting crop yield. Leaf colour is closely related to photosynthesis, and leaf colour mutants are considered an ideal material for studying photosynthesis.

RESULTS

We obtained a yellow-green leaf mutant jym165, using ethyl methane sulfonate (EMS) mutagenesis. Physiological and biochemical analyses indicated that the contents of chlorophyll a, chlorophyll b, carotenoids, and total chlorophyll in the jym165 mutant decreased significantly compared with those in Jiyu47 (JY47). The abnormal chloroplast development of jym165 led to a decrease in net photosynthetic rate and starch content compared with that of JY47. However, quality traits analysis showed that the sum of oil and protein contents in jym165 was higher than that in JY47. In addition, the regional yield (seed spacing: 5 cm) of jym165 increased by 2.42% compared with that of JY47 under high planting density. Comparative transcriptome analysis showed that the yellow-green leaf phenotype was closely related to photosynthesis and starch and sugar metabolism pathways. Genetic analysis suggests that the yellow-green leaf phenotype is controlled by a single recessive nuclear gene. Using Mutmap sequencing, the candidate regions related of leaf colour was narrowed to 3.44 Mb on Chr 10.

CONCLUSIONS

Abnormal chloroplast development in yellow-green mutants leads to a decrease in the photosynthetic pigment content and net photosynthetic rate, which affects the soybean photosynthesis pathway and starch and sugar metabolism pathways. Moreover, it has the potentiality to increase soybean yield under dense planting conditions. This study provides a useful reference for studying the molecular mechanisms underlying photosynthesis in soybean.

A-to-G/C/T and C-to-T/G/A dual-function base editor for creating multi-nucleotide variants

MNVs (multi-nucleotide variants) are critical genetic variants associated with various genetic diseases. However, tools for precisely installing MNVs are limited. In this study, we present the development of a dual-base editor, BDBE, by integrating TadA-dual and engineered human N-methylpurine DNA glycosylase (eMPG) into nCas9 (D10A). Our results demonstrate that BDBE effectively converts A-to-G/C/T (referred to as A-to-B) and C-to-T/G/A (referred to as C-to-D) simultaneously, yielding nine types of dinucleotides from adjacent CA nucleotides while maintaining minimal off-target effects. Notably, BDBE4 exhibits exceptional performance across multiple human cell lines and successfully simulated all nine dinucleotide MNVs from the gnomAD database. These findings indicate that BDBE significantly expands the product range of base editors and offers a valuable resource for advancing MNV research.

LncRNA-mRNA co-expression analysis reveals aquaporin-9-promoted neutrophil extracellular trap formation and inflammatory activation in sepsis

Sepsis is a life-threatening condition caused by an excessive inflammatory response to an infection. However, the precise regulatory mechanism of sepsis remains unclear. Using a strand-specific RNA-sequencing, we identified 115 hub differentially expressed long noncoding RNAs (lncRNAs) and 443 mRNAs in septic patients, primarily participated in crucial pathways including neutrophil extracellular trap (NET) formation and toll-like receptor signaling. Notably, NETs related gene aquaporin-9 (AQP9) and its associated lncRNAs exhibited significant upregulation in septic neutrophils. Functional experiments revealed AQP9 interacts with its lncRNAs to augment the formation of neutrophil NETs. In murine sepsis models, AQP9 inhibition with phloretin reduced proinflammatory cytokine production and lung damage. These findings provide crucial insights into the regulatory role of AQP9 in sepsis, unraveling its interaction with associated lncRNAs in transmitting downstream signals, holding promise in informing the development of novel therapeutic strategies aimed at ameliorating the debilitating effects of sepsis.

Transcriptome analysis reveals regulatory mechanism of postharvest softening in kiwiberry

BACKGROUND

Kiwiberry is an emerging edible fruit with market potential owing to its advantages of small size, thin and hairless skin, and sweet taste. However, kiwiberry is highly susceptible to softening after harvest, which poses a challenge for storage and transport. To reveal the underlying cause of kiwiberry softening, it is essential to investigate the characteristics of postharvest fruit and the molecular mechanisms that affect changes in fruit firmness.

RESULTS

Morphological observations and analysis of physical parameters showed that the skin of kiwiberry did not change markedly from the 1st to the 7th day after harvest, while the colour of the inner pericarp gradually turned yellow. By the 9th day of room temperature storage, the kiwiberries began to rot. The hardness decreased rapidly from the 1st to the 5th day postharvest, reaching the low level on the 5th day. The starch and pectin contents of kiwiberry showed a downward trend with increasing storage time. Transcriptome sequencing and weighted gene co-expression network analysis identified 29 key genes associated with the changes in the hardness of kiwiberry after harvest. Gene Ontology enrichment analysis indicated that these 29 genes are mainly involved in pectin metabolism, starch synthesis, starch decomposition, and starch metabolism. In addition, three transcription factors, AGL31, HAT14, and ALC, were identified to be strongly positively correlated with the 29 genes that affect the hardness changes of kiwiberry after harvest, and 28 of the 29 key genes were predicted to be regulated by HAT14.

CONCLUSIONS

These results reveal the changes in morphological characteristics and physiological indicators during the postharvest ripening and softening of kiwiberry stored under room temperature conditions. Transcriptome analysis identified 29 key genes and three transcription factors that affect the firmness changes of postharvest kiwiberry. The results of this study thus provide insight into the transcriptional regulatory mechanism of kiwiberry softening during storage to improve the postharvest quality.

Clopidogrel ameliorates high-fat diet-induced hepatic steatosis in mice through activation of the AMPK signaling pathway and beyond

Introduction

Metabolic dysfunction-associated steatotic liver disease (MASLD) frequently confers an increased risk of vascular thrombosis; however, the marketed antiplatelet drugs are investigated for the prevention and treatment of MASLD in patients with these coexisting diseases.

Methods

To determine whether clopidogrel could ameliorate high-fat diet (HFD)-induced hepatic steatosis in mice and how it works, mice were fed on normal diet or HFD alone or in combination with or without clopidogrel for 14 weeks, and primary mouse hepatocytes were treated with palmitate/oleate alone or in combination with the compounds examined for 24 h. Body weight, liver weight, insulin resistance, triglyceride and total cholesterol content in serum and liver, histological morphology, transcriptomic analysis of mouse liver, and multiple key MASLD-associated genes and proteins were measured, respectively.

Results and discussion

Clopidogrel mitigated HFD-induced hepatic steatosis (as measured with oil red O staining and triglyceride kit assay) and reduced elevations in serum aminotransferases, liver weight, and the ratio of liver to body weight. Clopidogrel downregulated the expression of multiple critical lipogenic (Acaca/Acacb, Fasn, Scd1, Elovl6, Mogat1, Pparg, Cd36, and Fabp4), profibrotic (Col1a1, Col1a2, Col3a1, Col4a1, Acta2, and Mmp2), and proinflammatory (Ccl2, Cxcl2, Cxcl10, Il1a, Tlr4, and Nlrp3) genes, and enhanced phosphorylation of AMPK and ACC. However, compound C (an AMPK inhibitor) reversed enhanced phosphorylation of AMPK and ACC in clopidogrel-treated primary mouse hepatocytes and alleviated accumulation of intracellular lipids. We concluded that clopidogrel may prevent and/or reverse HFD-induced hepatic steatosis in mice, suggesting that clopidogrel could be repurposed to fight fatty liver in patients.

The BET inhibitor JQ1 suppresses tumor survival by ABCB5-mediated autophagy in uveal melanoma

Uveal melanoma (UM), the most common adult ocular tumor, is aggressive and resistant to treatment, posing threat to patients' lives. The novel, effective therapies and the exploration of chemosensitizer for UM are imperative. The anticancer efficacy was evaluated with and without JQ1 treatment or ABCB5 gene silencing or overexpression. RNA sequencing identified downstream effectors in JQ1-treated cells. Integrated analysis of The Cancer Genome Atlas data (TCGA) and immunohistochemistry (IHC) revealed the oncogenic role of ABCB5. Functional analyses of JQ1 and defective ABCB5 were conducted using flow cytometry, transmission electron microscopy (TEM), IHC and western blot. The effects of JQ1 were validated in a heterotopic tumor model derived from OCM-1 cells. JQ1 inhibited cell proliferation, migration and invasion, induced cell cycle arrest and promoted apoptosis. JQ1 also suppressed the survival of UM in heterotopic tumor model. RNA sequencing indicated that JQ1 down-regulated the expressions of ABCB5 and autophagy-related genes, which was confirmed in vitro and in vivo by western blot. ABCB5, a marker associated with cancer stem cells and chemo-resistance, exhibited heightened expression in UM tissues, linked to immune infiltration. Notably, disrupting ABCB5 expression impeded UM cell proliferation and interfered with autophagy. Moreover, the overexpression of ABCB5 promoted cell proliferation, migration and invasion, and rescued autophagy related gene expression. Of note, JQ1 enhanced the sensitivity of OCM-1 cells to chemotherapy. Thus JQ1 inhibits UM survival via ABCB5-mediated autophagy and enhances chemo-sensitivity, suggesting potential for BET-based approaches in UM clinical management.

IL-1β-induced epithelial cell and fibroblast transdifferentiation promotes neutrophil recruitment in chronic rhinosinusitis with nasal polyps

Neutrophilic inflammation contributes to multiple chronic inflammatory airway diseases, including asthma and chronic rhinosinusitis with nasal polyps (CRSwNP), and is associated with an unfavorable prognosis. Here, using single-cell RNA sequencing (scRNA-seq) to profile human nasal mucosa obtained from the inferior turbinates, middle turbinates, and nasal polyps of CRSwNP patients, we identify two IL-1 signaling-induced cell subsets-LY6D+ club cells and IDO1+ fibroblasts-that promote neutrophil recruitment by respectively releasing S100A8/A9 and CXCL1/2/3/5/6/8 into inflammatory regions. IL-1β, a pro-inflammatory cytokine involved in IL-1 signaling, induces the transdifferentiation of LY6D+ club cells and IDO1+ fibroblasts from primary epithelial cells and fibroblasts, respectively. In an LPS-induced neutrophilic CRSwNP mouse model, blocking IL-1β activity with a receptor antagonist significantly reduces the numbers of LY6D+ club cells and IDO1+ fibroblasts and mitigates nasal inflammation. This study implicates the function of two cell subsets in neutrophil recruitment and demonstrates an IL-1-based intervention for mitigating neutrophilic inflammation in CRSwNP.

Screening of functional maternal-specific chromatin regulators in early embryonic development of zebrafish

The early stages of embryonic development rely on maternal products for proper regulation. However, screening for functional maternal-specific factors is challenging due to the time- and labor-intensive nature of traditional approaches. Here, we combine a computational pipeline and F0 null mutant technology to screen for functional maternal-specific chromatin regulators in zebrafish embryogenesis and identify Mcm3l, Mcm6l, and Npm2a as playing essential roles in DNA replication and cell division. Our results contribute to understanding the molecular mechanisms underlying early embryo development and highlight the importance of maternal-specific chromatin regulators in this critical stage.

Control of grain size and weight by the RNA-binding protein EOG1 in rice and wheat

Grain size is one of the important yield traits in crops. Understanding the molecular and genetic mechanisms of grain-size control is important for yield improvement. Here, we report that the enhancer of GS2AA (EOG1) encodes an RNA-binding protein, which can bind mRNAs of several grain-size genes and influence their abundance. The eog1-1 mutant produces large and heavy grains by promoting cell proliferation in the spikelet hull. OsGSK3 physically interacts with and phosphorylates EOG1, thereby influencing the stability of EOG1. Genetic analyses support that EOG1 and OsGSK3 share overlapped function in grain size and weight control but does so independently of GS2. Notably, genome editing of wheat homologs TaEOG1A/B/D causes large and heavy grains. Thus, our findings identify a genetic and molecular mechanism whereby the OsGSK3-EOG1 module regulates grain size and weight in rice, suggesting that this pathway has the potential for grain-size improvement in key crops.

Patient-derived rhabdomyosarcoma cells recapitulate the genetic and transcriptomic landscapes of primary tumors

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in childhood and adolescence. The availability of appropriate and well-characterized preclinical models for RMS is limited, posing a challenge for investigating the molecular mechanisms and evaluating new targeted compounds in preclinical settings. Here, we collected 51 RMS specimens (referred to as ZJUCH-RMS cohort) and established 9 patient-derived cells (PDCs) and validated the identity of these cells by the expression of RMS-specific markers. Whole-transcriptome analysis identified high-confidence mutations in ZJUCH-RMS cohort including RAS, TP53, ARID1A, MYOD1, and MYCN. Further studies showed that RMS PDCs retained the genetic alterations and the expression of RMS hallmark and dependency genes in matched primary tumors and acted as valuable tools to assess drug responses and pharmacogenomic interactions. Our study provides unique PDCs that are available for preclinical studies of RMS and further advances the feasibility of RMS PDCs as valuable tools for developing personalized treatments for patients.

A clinically effective model based on cell-free DNA methylation and low-dose CT for risk stratification of pulmonary nodules

Accurate, non-invasive, and cost-effective tools are needed to assist pulmonary nodule diagnosis and management due to increasing detection by low-dose computed tomography (LDCT). We perform genome-wide methylation sequencing on malignant and non-malignant lung tissues and designed a panel of 263 differential DNA methylation regions, which is used for targeted methylation sequencing on blood cell-free DNA (cfDNA) in two prospectively collected and retrospectively analyzed multicenter cohorts. We develop and optimize an integrative model for risk stratification of pulmonary nodules based on 40 cfDNA methylation biomarkers, age, and five simple computed tomography (CT) imaging features using machine learning approaches and validate its good performance in two cohorts. Using the two-threshold strategy can effectively reduce unnecessary invasive surgeries, overtreatment costs, and injury for patients with benign nodules while advising immediate treatment for patients with lung cancer, which can potentially improve the overall diagnosis of lung cancer following LDCT/CT screening.

Decoding molecular features of bovine oocyte fate during antral follicle growth via single-cell multi-omics analysis†

Antral follicle size is a useful predictive marker of the competency of enclosed oocytes for yielding an embryo following in vitro maturation and fertilization. However, the molecular mechanisms underpinning oocyte developmental potential during bovine antral follicle growth are still unclear. Here, we used a modified single-cell multi-omics approach to analyze the transcriptome, DNA methylome, and chromatin accessibility in parallel for oocytes and cumulus cells collected from bovine antral follicles of different sizes. Transcriptome profiling identified three types of oocytes (small, medium, and large) that underwent different developmental trajectories, with large oocytes exhibiting the largest average follicle size and characteristics resembling metaphase-II oocytes. Differential expression analysis and real-time polymerase chain reaction assay showed that most replication-dependent histone genes were highly expressed in large oocytes. The joint analysis of multi-omics data revealed that the transcription of 20 differentially expressed genes in large oocytes was associated with both DNA methylation and chromatin accessibility. In addition, oocyte-cumulus interaction analysis showed that inflammation, DNA damage, and p53 signaling pathways were active in small oocytes, which had the smallest average follicle sizes. We further confirmed that p53 pathway inhibition in the in vitro maturation experiments using oocytes obtained from small antral follicles could improve the quality of oocytes and increased the blastocyte rate after in vitro fertilization and culture. Our work provides new insights into the intricate orchestration of bovine oocyte fate determination during antral folliculogenesis, which is instrumental for optimizing in vitro maturation techniques to optimize oocyte quality.

The evolution of ectomycorrhizal symbiosis and host-plant switches are the main drivers for diversification of Amanitaceae (Agaricales, Basidiomycota)

BACKGROUND

Evolutionary radiation is widely recognized as a mode of species diversification, but the drivers of the rapid diversification of fungi remain largely unknown. Here, we used Amanitaceae, one of the most diverse families of macro-fungi, to investigate the mechanism underlying its diversification.

RESULTS

The ancestral state of the nutritional modes was assessed based on phylogenies obtained from fragments of 36 single-copy genes and stable isotope analyses of carbon and nitrogen. Moreover, a number of time-, trait-, and paleotemperature-dependent models were employed to investigate if the acquisition of ectomycorrhizal (ECM) symbiosis and climate changes promoted the diversification of Amanitaceae. The results indicate that the evolution of ECM symbiosis has a single evolutionary origin in Amanitaceae. The earliest increase in diversification coincided with the acquisition of the ECM symbiosis with angiosperms in the middle Cretaceous. The recent explosive diversification was primarily triggered by the host-plant switches from angiosperms to the mixed forests dominated by Fagaceae, Salicaceae, and Pinaceae or to Pinaceae.

CONCLUSIONS

Our study provides a good example of integrating phylogeny, nutritional mode evolution, and ecological analyses for deciphering the mechanisms underlying fungal evolutionary diversification. This study also provides new insights into how the transition to ECM symbiosis has driven the diversification of fungi.

Single-cell analysis reveals alternations between the aged and young mice prostates

BACKGROUND

Aging of the male prostate is an inevitable process in which the prostate undergoes hyperplasia, and this growth may lead to compression of the urethra, resulting in voiding dysfunction and associated symptoms, and an increased risk of prostate cancer. Despite the significance of prostate aging, the molecular mechanisms involved are still not fully understood.

METHODS

Prostate split by lobes from young (2 months) and aged (24 months) mice were collected for single-cell RNA sequencing (scRNA-seq) analysis. Tissues from both anterior prostate (AP) and ventral/dorsal/lateral prostate (VDLP) were included in the study. Data analysis included unsupervised clustering using the uniform manifold approximation and projection (UMAP) algorithm to identify distinct cell types based on marker gene expression. Differential gene expression analysis was performed to identify age-related changes in gene expression across different cell types. Functional enrichment analysis was conducted to elucidate biological pathways associated with differentially expressed genes. Additionally, cellular interactions and developmental trajectories were analyzed to characterize cellular dynamics during prostate aging.

RESULTS

The single-cell transcriptome analysis of the mouse prostate during aging revealed heterogeneity across various cell types and their changes during the aging process. We found a significant increase in the proportion of mesenchymal and immune cells in aged mice. Our study unveiled alterations in genes and pathways associated with cellular senescence, oxidative stress, and regeneration in epithelial cells. Furthermore, we observed that basal cells may undergo epithelial-mesenchymal transition (EMT) to become mesenchymal cells, particularly prominent in aged mice. Additionally, immune cells, notably macrophages and T cells, exhibited a heightened inflammatory response in aged mice.

CONCLUSION

In summary, our study provides a comparative analysis of the single-cell transcriptome of the aged and young mice prostates, elucidating cellular and molecular changes between the aged and young mice prostates.

Interleukin-34-orchestrated tumor-associated macrophage reprogramming is required for tumor immune escape driven by p53 inactivation

As the most frequent genetic alteration in cancer, more than half of human cancers have p53 mutations that cause transcriptional inactivation. However, how p53 modulates the immune landscape to create a niche for immune escape remains elusive. We found that cancer stem cells (CSCs) established an interleukin-34 (IL-34)-orchestrated niche to promote tumorigenesis in p53-inactivated liver cancer. Mechanistically, we discovered that Il34 is a gene transcriptionally repressed by p53, and p53 loss resulted in IL-34 secretion by CSCs. IL-34 induced CD36-mediated elevations in fatty acid oxidative metabolism to drive M2-like polarization of foam-like tumor-associated macrophages (TAMs). These IL-34-orchestrated TAMs suppressed CD8+ T cell-mediated antitumor immunity to promote immune escape. Blockade of the IL-34-CD36 axis elicited antitumor immunity and synergized with anti-PD-1 immunotherapy, leading to a complete response. Our findings reveal the underlying mechanism of p53 modulation of the tumor immune microenvironment and provide a potential target for immunotherapy of cancer with p53 inactivation.

A MSTNDel73C mutation with FGF5 knockout sheep by CRISPR/Cas9 promotes skeletal muscle myofiber hyperplasia

Mutations in the well-known Myostatin (MSTN) produce a 'double-muscle' phenotype, which makes it commercially invaluable for improving livestock meat production and providing high-quality protein for humans. However, mutations at different loci of the MSTN often produce a variety of different phenotypes. In the current study, we increased the delivery ratio of Cas9 mRNA to sgRNA from the traditional 1:2 to 1:10, which improves the efficiency of the homozygous mutation of biallelic gene. Here, a MSTNDel73C mutation with FGF5 knockout sheep, in which the MSTN and FGF5 dual-gene biallelic homozygous mutations were produced via the deletion of 3-base pairs of AGC in the third exon of MSTN, resulting in cysteine-depleted at amino acid position 73, and the FGF5 double allele mutation led to inactivation of FGF5 gene. The MSTNDel73C mutation with FGF5 knockout sheep highlights a dominant 'double-muscle' phenotype, which can be stably inherited. Both F0 and F1 generation mutants highlight the excellent trait of high-yield meat with a smaller cross-sectional area and higher number of muscle fibers per unit area. Mechanistically, the MSTNDel73C mutation with FGF5 knockout mediated the activation of FOSL1 via the MEK-ERK-FOSL1 axis. The activated FOSL1 promotes skeletal muscle satellite cell proliferation and inhibits myogenic differentiation by inhibiting the expression of MyoD1, and resulting in smaller myotubes. In addition, activated ERK1/2 may inhibit the secondary fusion of myotubes by Ca2+-dependent CaMKII activation pathway, leading to myoblasts fusion to form smaller myotubes.

A chromosome-level genome assembly of Cape hare (Lepus capensis)

The Cape hare (Lepus capensis) is among the most widely distributed hare species globally, inhabiting extensive regions across Africa, the Middle East, and Central Asia. However, evolutionary and genetic research on L. capensis was seriously impeded by the absence of a reference genome. Here, we assembled and constructed a chromosome-level genome of L. capensis (with scaffolds anchored to 25 chromosomes and a total assembled length of 2.9 Gb, achieving a contig N50 length of 124.44 Mb) using PacBio HiFi sequencing and Hi-C assembly technology. Evaluation using BUSCO indicated the genome assembly to be 98.2% complete. The de novo prediction revealed that repetitive sequences constitute 46.13% of the entire genome, and long interspersed nuclear elements (LINEs) constituted the largest portion. We annotated a total of 13, 868 protein-coding genes using transcriptomes from two tissues (muscle and skin). This high-quality reference genome serves as a valuable genomic resource for advancing genetic studies in this species.

Species-level understanding of the bacterial community in Daqu based on full-length 16S rRNA gene sequences

Daqu is used as the fermentation starter of Baijiu and contributes diversified functional microbes for saccharifying grains and converting sugars into ethanol and aroma components in Baijiu products. Daqu is mainly classified into three types, namely low (LTD), medium (MTD) and high (HTD) temperature Daqu, according to the highest temperatures reached in their fermentation processes. In this study, we used the PacBio small-molecule real-time (SMRT) sequencing technology to determine the full-length 16 S rRNA gene sequences from the metagenomes of 296 samples of different types of Daqu collected from ten provinces in China, and revealed the bacterial diversity at the species level in the Daqu samples. We totally identified 310 bacteria species, including 78 highly abundant species (with a relative abundance >0.1% each) which accounted for 91.90% of the reads from all the Daqu samples. We also recognized the differentially enriched bacterial species in different types of Daqu, and in the Daqu samples with the same type but from different provinces. Specifically, Lactobacillales, Enterobacterales and Bacillaceae were significantly enriched in the LTD, MTD and HTD groups, respectively. The potential co-existence and exclusion relationships among the bacteria species involved in all the Daqu samples and in the LTD, MTD and HTD samples from a specific region were also identified. These results provide a better understanding of the bacterial diversity in different types of Daqu at the species level.

Directed differentiation of human embryonic stem cells into conjunctival epithelial cells

Severe conjunctival damage can lead to extensive ocular cicatrisation, fornix shortening, and even ocular surface failure, resulting in significant vision impairment. Conjunctival reconstruction is the primary therapeutic strategy for these clinical conjunctival diseases. However, there have been limited studies on induced differentiation of conjunctival epithelial cells derived from stem cells. In this study, we established a chemical defined differentiation protocol from human embryonic stem cells (hESCs) into conjunctival epithelial cells. hES cell line H1 was used for differentiation, and RT-qPCR, immunofluorescence staining, Periodic-acid-Schiff staining (PAS), and transcriptome analysis were employed to identify the differentiated cells. Here, to imitate the development of the vertebrate conjunctiva, hESCs were induced using a three-step process involving first chetomin was used to induce ocular surface ectoderm, then nicotinamide was used to induce ocular surface epithelial progenitor cells, and finally epidermal growth factor, keratinocyte growth factor and other factors were used to differentiate mature conjunctival epithelial cells. hESC-derived conjunctival epithelial cells expressed mature conjunctival epithelial lineage markers (including PAX6, P63, K13). The presence of goblet cells was confirmed by positive PAS. Transcriptome analysis revealed that hESC-derived conjunctival epithelial cells possessed a more naïve phenotype, and exhibited greater proliferation capacity compared to mature human conjunctival epithelial cells, suggesting their potential as alternative seed cells for conjunctival reconstruction.

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

Integrative physiological and transcriptome analysis unravels the mechanism of low nitrogen use efficiency in burley tobacco

Burley tobacco, a chlorophyll-deficient mutant with impaired nitrogen use efficiency (NUE), generally requires three to five times more nitrogen fertilization than flue-cured tobacco to achieve a comparable yield, which generates serious environmental pollution and negatively affects human health. Therefore, exploring the mechanisms underlying NUE is an effective measure to reduce environmental pollution and an essential direction for burley tobacco plant improvement. Physiological and genetic factors affecting tobacco NUE were identified using two tobacco genotypes with contrasting NUE in hydroponic experiments. Nitrogen use inefficient genotype (TN90) had lower nitrogen uptake and transport efficiencies, reduced leaf and root biomass, lower nitrogen assimilation and photosynthesis capacity, and lower nitrogen remobilization ability than the nitrogen use efficient genotype (K326). Transcriptomic analysis revealed that genes associated with photosynthesis, carbon fixation, and nitrogen metabolism are implicated in NUE. Three nitrate transporter genes in the leaves (NPF2.11, NPF2.13, and NPF3.1) and three nitrate transporter genes (NPF6.3, NRT2.1, and NRT2.4) in roots were down-regulated by nitrogen starvation, all of which showed lower expression in TN90 than in K326. In addition, the protein-protein interaction (PPI) network diagram identified eight key genes (TPIP1, GAPB, HEMB, PGK3, PSBO, PSBP2, PSAG, and GLN2) that may affect NUE. Less advantageous changes in nitrogen uptake, nitrogen assimilation in combination with nitrogen remobilization, and maintenance of photosynthesis in response to nitrogen deficiency led to a lower NUE in TN90. The key genes (TPIP1, GAPB, PGK3, PSBO, PSBP2, PSAG, and GLN2) were associated with improving photosynthesis and nitrogen metabolism in tobacco plants grown under N-deficient conditions.

Analytical and clinical validation of a NGS panel in detecting targetable variants from ctDNA of metastatic NSCLC patients

BACKGROUND

Circulating tumor DNA (ctDNA) has emerged as a promising biomarker for noninvasive cancer diagnostics, particularly in the context of metastatic non-small-cell lung cancer (NSCLC). Detecting targetable variants through ctDNA analysis offers the potential to guide treatment decisions, especially in cases where tissue samples are insufficient or unavailable.

METHOD

In this study, we developed and validated a next-generation sequencing panel targeting 101 cancer-related genes (101-test) to detect somatic variants in ctDNA from a large cohort of Chinese patients with metastatic NSCLC. The performance of the 101-test was assessed by evaluating its limit of detection (LOD), accuracy, and precision in identifying molecular variants. Additionally, the concordance between ctDNA and tissue samples for detecting targetable variants was analyzed in 904 patients.

RESULTS

The 101-test demonstrated a LOD of 0.38% for single-nucleotide variants (SNVs), 0.33% for insertions and deletions (InDels), and 0.33% for fusions. Sensitivity was 98.3% for SNVs, 100% for InDels, and 100% for fusions when compared to digital droplet PCR (ddPCR)/breakpoint PCR reference methods. The by-variant sensitivity for somatic variants was 97.5%, with a specificity of 99.9% between tumor-only and tumor-normal analyses. In a real-world cohort, the concordance between ctDNA and tissue samples for identifying targetable variants was 72.2% (457/633). Notably, the EGFR S768I variant, recently recommended by clinical guidelines, achieved an 80% concordance rate. Furthermore, 4.3% of patients (27/633) with targetable variants were identified exclusively through ctDNA testing.

CONCLUSION

The ctDNA-based 101-test is a highly sensitive and specific tool for detecting targetable variants in metastatic NSCLC, particularly in cases with insufficient tissue samples. The findings support the use of ctDNA testing as a reliable and complementary method to traditional tissue-based molecular analysis, enhancing the precision of treatment strategies for NSCLC patients.

Diabetes Mellitus Inhibits Hair Follicle Regeneration by Inducing Macrophage Reprogramming-Mediated Pyroptosis

Background

Diabetes mellitus (DM) is known to inhibit skin self-renewal and hair follicle stem cell (HFSC) activation, which may be key in the formation of chronic diabetic wounds. This study aimed to investigate the reasons behind the suppression of HFSC activation in DM mice.

Methods

Type 1 DM (T1DM) was induced in 6-week-old mice via streptozotocin, and hair follicle growth was subsequently monitored. RNA sequencing, bioinformatics analyses, qRT‒PCR, immunostaining, and cellular experiments were carried out to investigate the underlying mechanisms involved.

Results

T1DM inhibited HFSC activation, which correlated with an increase in caspase-dependent programmed cell death. Additionally, T1DM triggered apoptosis and pyroptosis, predominantly in HFSCs and epidermal regions, with pyroptosis being more pronounced in the inner root sheath of hair follicles. Notably, significant cutaneous immune imbalances were observed, particularly in macrophages. Cellular experiments demonstrated that M1 macrophages inhibited HaCaT cell proliferation and induced cell death, whereas high-glucose environments alone did not have the same effect.

Conclusion

T1DM inhibits HFSC activation via macrophage reprogramming-mediated caspase-dependent pyroptosis, and there is a significant regional characterization of cell death. Moreover, T1DM-induced programmed cell death in the skin may be more closely related to immune homeostasis imbalance than to hyperglycemia itself. These findings shed light on the pathogenesis of diabetic ulcers and provide a theoretical basis for the use of hair follicle grafts in wound repair.

Identification of key lncRNAs and mRNAs in muscle development pathways of Tan sheep

The study aimed to identify the long noncoding RNA (lncRNA) responsible for regulating muscle development in Tan sheep. RNA-seq analysis was conducted on longissimus dorsi samples from 1-day-old and 60-day-old Tan sheep to investigate the molecular processes involved in muscle development. A total of 5517 lncRNAs and 2885 mRNAs were found to be differentially expressed in the 60-day-old Tan sheep. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these differentially expressed lncRNAs and mRNAs were linked to pathways crucial for muscle development, such as MAPK, cAMP, and calcium-mediated signaling pathways. Key genes like CDKN1A, MAPK14, TGFB1, MEF2C, MYOD1, and CD53 were identified as significant players in muscle development. The study validated the RNA-seq results through RT-qPCR, confirming the consistency of expression levels of differentially expressed lncRNAs and mRNAs. These findings indicate that lncRNA-mRNA networks produce a remarked effect on modulating muscle development in Tan sheep, such as lncRNAs (MSTRG.12808.1/MSTRG.22662.3/MSTRG.18310.1) and mRNAs (MSTRG.10027/MSTRG.10029/MSTRG.10258/MSTRG.11011/MSTRG.10354), laying the groundwork for future research in this area.

Neolithic to Bronze Age human maternal genetic history in Yunnan, China

Yunnan in southwest China is a geographically and ethnically complex region at the intersection of southern China and Southeast Asia, and a focal point for human migrations. To clarify its maternal genetic history, we generated 152 complete mitogenomes from 17 Yunnan archaeological sites. Our results reveal distinct genetic histories segregated by geographical regions. Maternal lineages of ancient populations from northwestern and northern Yunnan exhibit closer affinities with past and present-day populations from northern East Asia and Tibet, providing important genetic evidence for the migration and interaction of populations along the Tibetan-Yi corridor since the Neolithic. Between 5500 to 1800 years ago, central Yunnan populations maintained their internal genetic relationships, including a 7000-year-old basal lineage of the rare and widely dispersed haplogroup M61. At the Xingyi site, changes in mitochondrial DNA haplogroups occurred between the Late Neolithic and Bronze Age, with haplogroups shifting from those predominant in the Yellow River region to those predominant in coastal southern China. These results highlight the high diversity of Yunnan populations during the Neolithic to Bronze Age.

De novo genome of thick-tailed pygmy jerboa Salpingotus crassicauda, Cardiocraniinae, Dipodidae

Pygmy jerboas are one of the smallest taxa of rodents. They exhibit distinctly different morphological and biological characteristics from other subfamilies, such as more restricted distribution, species richness, reproductive ability, and population size. Agricultural expansion and the development of new energy projects in recent years lead to sharp decline of their natural populations. Here, we assembled and annotated the first reference genome for the subfamily Cardiocraniinae using Illunima and Nanopore sequencing from the thick-tailed pygmy jerboa, Salpingotus crassicauda. The final genome is 2.44 Gb in size, with a contig N50 length of 13.71 Mb and a BUSCO completeness of 96.35%. A total of 23,344 protein-coding genes were annotated in the final genome. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and 2 rRNA. These genomic assemblies provide resources in studying phylogeny and adaptive evolution of Dipodidae, as well as implementing conservation management of jerboas.

Integrative Analysis of Transcriptome and Metabolome Reveals the Pivotal Role of the NAM Family Genes in Oncidium hybridum Lodd. Pseudobulb Growth

Oncidium hybridum Lodd. is an important ornamental flower that is used as both a cut flower and a potted plant around the world; additionally, its pseudobulbs serve as essential carriers for floral organs and flower development. The NAM gene family is crucial for managing responses to various stresses as well as regulating growth in plants. However, the mechanisms by which NAM genes regulate the development of pseudobulbs remain unclear. In this study, a total of 144 NAM genes harboring complete structural domains were identified in O. hybridum. The 144 NAM genes were systematically classified into 14 distinct subfamilies via phylogenetic analysis. Delving deeper into the conserved motifs revealed that motifs 1-6 exhibited remarkable conservation, while motifs 7-10 presented in a few NAM genes only. Notably, NAM genes sharing identical specific motifs were classified into the same subfamily, indicating functional relatedness. Furthermore, the examination of occurrences of gene duplication indicated that the NAM genes display 16 pairs of tandem duplications along with five pairs of segmental duplications, suggesting their role in genetic diversity and potential adaptive evolution. By conducting a correlation analysis integrating transcriptomics and metabolomics at four stages of pseudobulb development, we found that OhNAM023, OhNAM030, OhNAM007, OhNAM019, OhNAM083, OhNAM047, OhNAM089, and OhNAM025 exhibited significant relationships with the endogenous plant hormones jasmonates (JAs), hinting at their potential involvement in hormonal signaling. Additionally, OhNAM089, OhNAM025, OhNAM119, OhNAM055, and OhNAM136 showed strong links with abscisic acid (ABA) and abscisic acid glucose ester (ABA-GE), suggesting the possible regulatory function of these NAM genes in plant growth and stress responses. The 144 NAM genes identified in this study provide a basis for subsequent research and contribute to elucidating the intricate molecular mechanisms of NAM genes in Oncidium and potentially in other species.

Haplotype-resolved genome assembly of the upas tree (Antiaris toxicaria)

The upas tree (Antiaris toxicaria Lesch.) is a medically important plant that contains various specialized metabolites with significant bioactivity. The lack of a reference genome hinders the in-depth study as well as rational exploitation and conservation of this plant. Here, we present the first holotype-resolved chromosome-scale genome of the upas tree. The assembled genome consisted of 26 chromosomes that contain 1.34 Gb of sequencing data with a contig N50 length of 60 Mb. Genome annotation identified 43,500 protein-coding genes in the upas tree genome, of which 98.75% were functionally annotated. This high-quality reference genome will lay the foundation for further studies on the evolution and functional genomics of the upas tree.

Comparative DNA methylation reveals epigenetic adaptation to high altitude in snub-nosed monkeys

DNA methylation plays a crucial role in environmental adaptations. Here, using whole-genome bisulfite sequencing, we generated comprehensive genome-wide DNA methylation profiles for the high-altitude Yunnan snub-nosed monkey ( Rhinopithecus bieti) and the closely related golden snub-nosed monkey ( R. roxellana). Our findings indicated a slight increase in overall DNA methylation levels in golden snub-nosed monkeys compared to Yunnan snub-nosed monkeys, suggesting a higher prevalence of hypermethylated genomic regions in the former. Comparative genomic methylation analysis demonstrated that genes associated with differentially methylated regions were involved in membrane fusion, vesicular formation and trafficking, hemoglobin function, cell cycle regulation, and neuronal differentiation. These results suggest that the high-altitude-related epigenetic modifications are extensive, involving a complete adaptation process from the inhibition of single Ca 2+ channel proteins to multiple proteins collaboratively enhancing vesicular function or inhibiting cell differentiation and proliferation. Functional assays demonstrated that overexpression or down-regulation of candidate genes, such as SNX10, TIMELESS, and CACYBP, influenced cell viability under stress conditions. Overall, this research suggests that comparing DNA methylation across closely related species can identify novel candidate genomic regions and genes associated with local adaptations, thereby deepening our understanding of the mechanisms underlying environmental adaptations.

Decoding the pathogenesis of spermatogenic failure in cryptorchidism through single-cell transcriptomic profiling

Cryptorchidism, commonly known as undescended testis, affects 1%-9% of male newborns, posing infertility and testis tumor risks. Despite its prevalence, the detailed pathophysiology underlying male infertility within cryptorchidism remains unclear. Here, we profile and analyze 46,644 single-cell transcriptomes from individual testicular cells obtained from adult males diagnosed with cryptorchidism and healthy controls. Spermatogenesis compromise in cryptorchidism links primarily to spermatogonium self-renewal and differentiation dysfunctions. We illuminate the involvement of testicular somatic cells, including immune cells, thereby unveiling the activation and degranulation of mast cells in cryptorchidism. Mast cells are identified as contributors to interstitial fibrosis via transforming growth factor β1 (TGF-β1) and cathepsin G secretion. Furthermore, significantly increased levels of secretory proteins indicate mast cell activation and testicular fibrosis in the seminal plasma of individuals with cryptorchidism compared to controls. These insights serve as valuable translational references, enriching our comprehension of testicular pathogenesis and informing more precise diagnosis and targeted therapeutic strategies for cryptorchidism.

Profiles and interactions of gut microbiome and intestinal microRNAs in pediatric Crohn's disease

Gut dysbiosis is closely related to dysregulated microRNAs (miRNAs) in the intestinal epithelial cells, which plays an important role in the pathogenesis of Crohn's disease (CD). We investigated the relationship between fecal gut microbiome (GM) and intestinal tissue miRNAs in different stages of pediatric CD. Metagenomic analysis and miRNA sequencing were conducted to examine the GM and intestinal miRNA profiles of CD patients before and after clinical induction therapy and the controls. Twenty-seven newly diagnosed, therapy-naïve pediatric patients with active CD and 11 non-inflammatory bowel disease (IBD) controls were recruited in this study. Among CD patients, 11 patients completed induction treatment and reached clinical remission. Both GM and miRNA profiles were significantly changed between CD patients and controls. Seven key bacteria were identified at species level including Defluviitalea raffinosedens, Thermotalea metallivorans, Roseburia intestinalis, Dorea sp. AGR2135, Escherichia coli, Shigella sonnei, and Salmonella enterica, the exact proportions of which were further validated by real-time quantitative PCR analysis. Eight key miRNAs were also identified including hsa-miR-215-5p, hsa-miR-194-5p, hsa-miR-12135, hsa-miR-509-3-5p, hsa-miR-212-5p, hsa-miR-4448, hsa-miR-501-3p, and hsa-miR-503-5p. The functional enrichment analysis of differential miRNAs indicated the significantly altered cyclin protein, cyclin-dependent protein, and cell cycle pathway. The close interactions between seven key bacteria and eight key miRNAs were further investigated by miRNA target prediction. The association between specific miRNA expressions and key gut bacteria at different stages of CD supported their important roles as potential molecular biomarkers. Understanding the relationship between them will help us to explore the molecular mechanisms of CD.

IMPORTANCE

Since previous studies have focused on the change of the fecal gut microbiome and intestinal tissue miRNA in pediatric Crohn's disease (CD), the relationship between them in different stages is still not clear. This is the first study to explore the gut microbiota and miRNA and their correlations with the Pediatric Crohn's Disease Activity Index (PCDAI). Crohn's Disease Endoscopic Index of Severity (CDEIS), and calprotectin, by applying two omics approach in three different groups (active CD, CD in remission with exclusive enteral nutrition or infliximab induction therapy, and the healthy controls). Both gut microbiome structure and the miRNA profiles were significantly changed in the different stage of CD. Seven key gut microbiome at species and eight key miRNAs were found, and their close interactions were further fully investigated by miRNA target prediction.

Biochar regulates the functions of keystone taxa to reduce p-coumaric acid accumulation in soil

Introduction

Applying biochar (BC) to reduce toxic substance accumulation in soil, either through direct adsorption or modulation of the microbial community, has received considerable attention. However, a knowledge gap exists regarding how BC regulates microbial community structure and functions to mitigate toxic substance accumulation.

Methods

We previously identified p-coumaric acid (p-CA) as a representative autotoxin in tobacco rhizosphere soil. On this basis, this study simulated a soil environment with p-CA accumulation to investigate the impacts of BC on p-CA, soil physicochemical properties, and microbial community structure and function.

Results

The results showed that p-CA could be directly adsorbed onto BC, which followed the pseudo-second-order kinetic model (R 2 = 0.996). A pot experiment revealed that BC significantly reduced soil p-CA, altered soil microbial composition, and enhanced bacterial community diversity. A weighted correlation network analysis showed a close association between taxon 1 in the microbial network and p-CA, suggesting a pivotal role for this taxon in reducing p-CA, with Devosia and Nocardioides identified as potential key contributors to this process. The prediction of possible keystone taxa functions showed that BC increased the relative abundances of aromatic compound degraders. Mantel tests indicated that soil organic matter exerted the greatest influence on keystone taxa functions and hub genera.

Discussion

These findings suggest that BC may either directly chemisorb p-CA or indirectly facilitate p-CA degradation by regulating the functioning of keystone taxa. The results of this study provide a novel perspective for further investigation of the mechanisms through which BC reduces the accumulation of toxic substances in soil.

The response patterns of r- and K-strategist bacteria to long-term organic and inorganic fertilization regimes within the microbial food web are closely linked to rice production

Soil microbial food web is crucial for maintaining crop production, while its community structure varies among fertilization regimes. Currently, the mechanistic understanding of the relationships between microbial food web and crop production under various nutrient fertilizations is poor. This knowledge gap limits our capacity to achieve precision agriculture for ensuring yield stability. In this study, we investigated the abiotic (i.e., soil chemical properties) and biotic factors (i.e., microbial food web, including bacteria, fungi, archaea and nematodes) that were closely associated with rice (Oryza sativa L.) production, using soils from seven fertilization regimes in distinct sampling locations (i.e., bulk vs rhizosphere soil) at a long-term experimental site. Organic manure alone fertilization (M) and integrated fertilization (NPKM) combining manure with inorganic fertilizers increased soil pH by 0.21-0.41 units and organic carbon content by 49.1 %-65.2 % relative to the non-fertilization (CK), which was distinct with inorganic fertilization. The principal coordinate analysis (PCoA) revealed that soil microbial and nematode communities were primarily shaped by fertilization rather than sampling locations. Organic fertilization (M, NPKM) increased the relative abundance of both r-strategist bacteria, specific taxa within the fungal (i.e., Pezizales) and nematode communities (i.e., omnivores-predators), whereas inorganic fertilization increased K-strategist bacteria abundances relative to the CK. Correspondingly, network analysis showed that the keystone taxa in the amplicon sequence variants (ASVs) enriched by organic manure and inorganic fertilization were mainly affiliated with r- and K-strategist bacteria, respectively. Structural equation modeling (SEM) analysis found that r- and K-strategist bacteria were positively correlated with rice production under organic and inorganic fertilization, respectively. Our results demonstrate that the response patterns of r/K-strategists to nutrient fertilization largely regulate rice yield, suggesting that the enhanced soil fertility and r-strategists contribute to the highest crop production in NPKM fertilization.

Microglia in the hypothalamic paraventricular nucleus sense hemodynamic disturbance and promote sympathetic excitation in hypertension

Hypertension is usually accompanied by elevated sympathetic tonicity, but how sympathetic hyperactivity is triggered is not clear. Recent advances revealed that microglia-centered neuroinflammation contributes to sympathetic excitation in hypertension. In this study, we performed a temporospatial analysis of microglia at both morphological and transcriptomic levels and found that microglia in the hypothalamic paraventricular nucleus (PVN), a sympathetic center, were early responders to hypertensive challenges. Vasculature analyses revealed that the PVN was characterized by high capillary density, thin vessel diameter, and complex vascular topology relative to other brain regions. As such, the PVN was susceptible to the penetration of ATP released from the vasculature in response to hemodynamic disturbance after blood pressure increase. Mechanistically, ATP ligation to microglial P2Y12 receptor was responsible for microglial inflammatory activation and the eventual sympathetic overflow. Together, these findings identified a distinct vasculature pattern rendering vulnerability of PVN pre-sympathetic neurons to hypertension-associated microglia-mediated inflammatory insults.

A 3-component module maintains sepal flatness in Arabidopsis

As in origami, morphogenesis in living systems heavily relies on tissue curving and folding through the interplay between biochemical and biomechanical cues. By contrast, certain organs maintain their flat posture over several days. Here, we identified a pathway that is required for the maintenance of organ flatness, taking the sepal, the outermost floral organ, in Arabidopsis as a model system. Through genetic, cellular, and mechanical approaches, our results demonstrate that the global gene expression regulator VERNALIZATION INDEPENDENCE 4 (VIP4) fine-tunes the mechanical properties of sepal cell walls and maintains balanced growth on both sides of the sepals, mainly by orchestrating the distribution pattern of AUXIN RESPONSE FACTOR 3 (ARF3). vip4 mutation results in softer cell walls and faster cell growth on the adaxial sepal side, which eventually cause sepals to bend outward. Downstream of VIP4, ARF3 works through modulating auxin to downregulate pectin methylesterase VANGUARD1, resulting in decreased cell wall stiffness. Thus, our work unravels a 3-component module that relates hormonal patterns to organ curvature and actively maintains sepal flatness during its growth.

Integrated landscape of plasma metabolism and proteome of patients with post-traumatic deep vein thrombosis

Deep vein thrombosis (DVT) is a leading cause of morbidity and mortality after trauma. Here, we integrate plasma metabolomics and proteomics to evaluate the metabolic alterations and their function in up to 680 individuals with and without DVT after trauma (pt-DVT). We identify 28 metabolites and 2 clinical parameter clusters associated with pt-DVT. Then, we develop a panel of 9 metabolites (hexadecanedioic acid, pyruvic acid, L-Carnitine, serotonin, PE(P-18:1(11Z)/18:2(9Z,12Z)), 3-Hydroxycapric acid, 5,6-DHET, 3-Methoxybenzenepropanoic acid and pentanenitrile) that can predict pt-DVT with high performance, which can be verified in an independent cohort. Furthermore, the integration analysis of metabolomics and proteomics data indicates that the upregulation of glycolysis/gluconeogenesis-TCA cycle may promote thrombosis by regulating ROS levels in red blood cells, suggesting that interfering with this process might be potential therapeutic strategies for pt-DVT. Together, our study comprehensively delineates the metabolic and hematological dysregulations for pt-DVT, and provides potential biomarkers for early detection.

Unraveling the key role of chromatin structure in cancer development through epigenetic landscape characterization of oral cancer

Epigenetic alterations, such as those in chromatin structure and DNA methylation, have been extensively studied in a number of tumor types. But oral cancer, particularly oral adenocarcinoma, has received far less attention. Here, we combined laser-capture microdissection and muti-omics mini-bulk sequencing to systematically characterize the epigenetic landscape of oral cancer, including chromatin architecture, DNA methylation, H3K27me3 modification, and gene expression. In carcinogenesis, tumor cells exhibit reorganized chromatin spatial structures, including compromised compartment structures and altered gene-gene interaction networks. Notably, some structural alterations are observed in phenotypically non-malignant paracancerous but not in normal cells. We developed transformer models to identify the cancer propensity of individual genome loci, thereby determining the carcinogenic status of each sample. Insights into cancer epigenetic landscapes provide evidence that chromatin reorganization is an important hallmark of oral cancer progression, which is also linked with genomic alterations and DNA methylation reprogramming. In particular, regions of frequent copy number alternations in cancer cells are associated with strong spatial insulation in both cancer and normal samples. Aberrant methylation reprogramming in oral squamous cell carcinomas is closely related to chromatin structure and H3K27me3 signals, which are further influenced by intrinsic sequence properties. Our findings indicate that structural changes are both significant and conserved in two distinct types of oral cancer, closely linked to transcriptomic alterations and cancer development. Notably, the structural changes remain markedly evident in oral adenocarcinoma despite the considerably lower incidence of genomic copy number alterations and lesser extent of methylation alterations compared to squamous cell carcinoma. We expect that the comprehensive analysis of epigenetic reprogramming of different types and subtypes of primary oral tumors can provide additional guidance to the design of novel detection and therapy for oral cancer.