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

Extracellular cancer‑associated fibroblasts: A novel subgroup in the cervical cancer microenvironment that exhibits tumor‑promoting roles and prognosis biomarker functions

Tumor invasion and metastasis are the processes that primarily cause adverse outcomes in patients with cervical cancer. Cancer-associated fibroblasts (CAFs), which participate in cancer progression and metastasis, are novel targets for the treatment of tumors. The present study aimed to assess the heterogeneity of CAFs in the cervical cancer microenvironment through single-cell RNA sequencing. After collecting five cervical cancer samples and obtaining the CAF-associated gene sets, the CAFs in the cervical cancer microenvironment were divided into myofibroblastic CAFs and extracellular (ec)CAFs. The ecCAFs appeared with more robust pro-tumorigenic effects than myCAFs according to enrichment analysis. Subsequently, through combining the ecCAF hub genes and bulk gene expression data for cervical cancer obtained from The Cancer Genome Atlas and Gene Ontology databases, univariate Cox regression and least absolute shrinkage and selection operator analyses were performed to establish a CAF-associated risk signature for patients with cancer. The established risk signature demonstrated a stable and strong prognostic capability in both the training and validation cohorts. Subsequently, the association between the risk signature and clinical data was evaluated, and a nomogram to facilitate clinical application was established. The risk score was demonstrated to be associated with both the tumor immune microenvironment and the therapeutic responses. Moreover, the signature also has predictive value for the prognosis of head and neck squamous cell carcinoma, and bladder urothelial carcinoma, which were also associated with human papillomavirus infection. In conclusion, the present study assessed the heterogeneity of CAFs in the cervical cancer microenvironment, and a subgroup of CAFs that may be closely associated with tumor progression was defined. Moreover, a signature based on the hub genes of ecCAFs was shown to have biomarker functionality in terms of predicting survival rates, and therefore this CAF subgroup may become a therapeutic target for cervical cancer in the future.

Curation, nomenclature, and topological classification of receptor-like kinases from 528 plant species for novel domain discovery and functional inference

Receptor-like kinases (RLKs) are the most numerous signal transduction components in plants and play important roles in determining how different plants adapt to their ecological environments. Research on RLKs has focused mainly on a small number of typical RLK members in a few model plants. There is an urgent need to study the composition, distribution, and evolution of RLKs at the holistic level to increase our understanding of how RLKs assist in the ecological adaptations of different plant species. In this study, we collected the genome assemblies of 528 plant species and constructed an RLK dataset. Using this dataset, we identified and characterized 524 948 RLK family members. Each member underwent systematic topological classification and was assigned a gene ID based on a unified nomenclature system. Furthermore, we identified two novel extracellular domains in some RLKs, designated Xiao and Xiang. Evolutionary analysis of the RLK family revealed that the RLCK-XVII and RLCK-XII-2 classes were present exclusively in dicots, suggesting that diversification of RLKs between monocots and dicots may have led to differences in downstream cytoplasmic responses. We also used an interaction proteome to help empower data mining for inference of new RLK functions from a global perspective, with the ultimate goal of understanding how RLKs shape the adaptation of different plants to the environments/ecosystems. The assembled RLK dataset, together with annotations and analytical tools, forms an integrated foundation of multiomics data that is publicly accessible via the metaRLK web portal (http://metaRLK.biocloud.top).

FSP1 inhibition enhances olaparib sensitivity in BRCA-proficient ovarian cancer patients via a nonferroptosis mechanism

Poly ADP-ribose polymerase inhibitors (PARPis) exhibit promising efficacy in patients with BRCA mutations or homologous repair deficiency (HRD) in ovarian cancer (OC). However, less than 40% of patients have HRD, it is vital to expand the indications for PARPis in BRCA-proficient patients. Ferroptosis suppressor protein 1 (FSP1) is a key protein in a newly identified ferroptosis-protective mechanism that occurs in parallel with the GPX4-mediated pathway and is associated with chemoresistance in several cancers. Herein, FSP1 is reported to be negatively correlated with the prognosis in OC patients. Combination therapy comprising olaparib and iFSP1 (a FSP1 inhibitor) strongly inhibited tumour proliferation in BRCA-proficient OC cell lines, patient-derived organoids (PDOs) and xenograft mouse models. Surprisingly, the synergistic killing effect could not be reversed by ferroptosis inhibitors, indicating that mechanisms other than ferroptosis were responsible for the synergistic lethality. In addition, cotreatment was shown to induce increased γH2A.X foci and to impair nonhomologous end joining (NHEJ) activity to a greater extent than did any single drug. Mass spectrometry and immunoprecipitation analyses revealed that FSP1 interacted with Ku70, a classical component recruited to and occupying the end of double-strand breaks (DSBs) in the NHEJ process. FSP1 inhibition decreased Ku70 PARylation, impaired subsequent DNA-PKcs recruitment to the Ku complex at DSB sites and was rescued by restoring PARylation. These findings unprecedentedly reveal a novel role of FSP1 in DNA damage repair and provide new insights into how to sensitize OC patients to PARPi treatment.

Haplotype-resolved genome assembly of Phanera championii reveals molecular mechanisms of flavonoid synthesis and adaptive evolution

Phanera championii is a medicinal liana plant that has successfully adapted to hostile karst habitats. Despite extensive research on its medicinal components and pharmacological effects, the molecular mechanisms underlying the biosynthesis of critical flavonoids and its adaptation to karst habitats remain elusive. In this study, we performed high-coverage PacBio and Hi-C sequencing of P. championii, which revealed its high heterozygosity and phased the genome into two haplotypes: Hap1 (384.60 Mb) and Hap2 (383.70 Mb), encompassing a total of 58 612 annotated genes. Comparative genomes analysis revealed that P. championii experienced two whole-genome duplications (WGDs), with approximately 59.59% of genes originating from WGD events, thereby providing a valuable genetic resource for P. championii. Moreover, we identified a total of 112 genes that were strongly positively selected. Additionally, about 81.60 Mb of structural variations between the two haplotypes. The allele-specific expression patterns suggested that the dominant effect of P. championii was the elimination of deleterious mutations and the promotion of beneficial mutations to enhance fitness. Moreover, our transcriptome and metabolome analysis revealed alleles in different tissues or different haplotypes collectively regulate the synthesis of flavonoid metabolites. In summary, our comprehensive study highlights the significance of genomic and morphological adaptation in the successful adaptation of P. championii to karst habitats. The high-quality phased genomes obtained in this study serve as invaluable genomic resources for various applications, including germplasm conservation, breeding, evolutionary studies, and elucidation of pathways governing key biological traits of P. championii.

Effects of different temperatures on Leiocassis longirostris gill structure and intestinal microbial composition

Fish are poikilothermic vertebrates and their physiological activities are affected by water temperature. In recent years, extreme weather has occurred frequently, and temperature changes have adversely affected the growth of farmed fish. To explore the changes in gill tissue structure caused by changing the water temperature and the relationship between the intestinal microbiota and the Leiocassis longirostris host adaptation mechanism, gill tissue sections and intestinal microbial 16S rRNA amplicon sequencing were conducted under different temperature stress (low temperature 4 °C, normal temperature 26 °C and high temperature 32 °C). The results showed that heat stress and cold stress caused injury and swelling, terminal congestion, cell vacuolation, and necrosis of the gill tissue of L. longirostris. For intestinal microbiota, the abundances of Pseudomonadota and Bacillota increased at the cold stress, while the abundances of Fusobacteriota and Bacteroidota increased at the heat stress. The number of opportunistic bacteria, mainly Aeromonas and Acinetobacter, was the highest under cold stress. In addition, the richness of the intestinal microbiota decreased significantly at heat and cold stresses, while evenness increased. Prediction of intestinal microbiota function showed that most common functions, such as metabolism of cofactors and vitamins, energy metabolism and replication and repair, were decreased significantly at heat stress and cold stress, and phylogenetic relationship analysis revealed significant differences among the groups. In conclusion, the change of temperature altered the gill tissue structure, and affected the structure and homeostasis of the intestinal microbiota, thus affecting the survival time of L. longirostris, and cold stress had a greater effect than heat stress.

The immunomodulatory impact of naturally derived neem leaf glycoprotein on the initiation progression model of 4NQO induced murine oral carcinogenesis: a preclinical study

Introduction

Murine tumor growth restriction by neem leaf glycoprotein (NLGP) was established in various transplanted models of murine sarcoma, melanoma and carcinoma. However, the role of NLGP in the sequential carcinogenic steps has not been explored. Thus, tongue carcinogenesis in Swiss mice was induced by 4-nitroquinoline-1-oxide (4NQO), which has close resemblance to human carcinogenesis process. Interventional role of NLGP in initiation-promotion protocol established during 4NQO mediated tongue carcinogenesis in relation to systemic immune alteration and epithelial-mesenchymal transition (EMT) is investigated.

Methods

4NQO was painted on tongue of Swiss mice every third day at a dose of 25µl of 5mg/ml stock solution. After five consecutive treatment with 4NQO (starting Day7), one group of mice was treated with NLGP (s.c., 25µg/mice/week), keeping a group as PBS control. Mice were sacrificed in different time-intervals to harvest tongues and studied using histology, immunohistochemistry, flow-cytometry and RT-PCR on different immune cells and EMT markers (e-cadherin, vimentin) to elucidate their phenotypic and secretory status.

Results

Local administration of 4NQO for consecutive 300 days promotes significant alteration in tongue mucosa including erosion in papillae and migration of malignant epithelial cells to the underlying connective tissue stroma with the formation of cell nests (exophytic-hyperkeratosis with mild dysplasia). Therapeutic NLGP treatment delayed pre-neoplastic changes promoting normalization of mucosa by maintaining normal structure. Flow-cytometric evidences suggest that NLGP treatment upregulated CD8+, IFNγ+, granzyme B+, CD11c+ cells in comparison to 4NQO treated mice with a decrease in Ki67+ and CD4+FoxP3+ cells in NLGP treated cohort. RT-PCR demonstrated a marked reduction of MMP9, IL-6, IL-2, CD31 and an upregulation in CCR5 in tongues from 4NQO+NLGP treated mice in comparison to 4NQO treated group. Moreover, 4NQO mediated changes were associated with reduction of e-cadherin and simultaneous up-regulation of vimentin expression in epithelium that was partially reversed by NLGP.

Discussion

Efficacy of NLGP was tested first time in sequential carcinogenesis model and proved effective in delaying the initial progression. NLGP normalizes type 1 immunity including activation of the CD8+T effector functions, reduction of regulatory T cell functions, along with changes in EMT to make the host systemically alert to combat the carcinogenic threat.

The aldehyde dehydrogenase 2 rs671 variant enhances amyloid β pathology

In the ALDH2 rs671 variant, a guanine changes to an adenine, resulting in a dramatic decrease in the catalytic activity of the enzyme. Population-based data are contradictory about whether this variant increases the risk of Alzheimer's disease. In East Asian populations, the prevalence of the ALDH2 rs671 variant is 30-50%, making the National Human Brain Bank for Development and Function (the largest brain bank in East Asia) an important resource to explore the link between the ALDH2 rs671 polymorphism and Alzheimer's disease pathology. Here, using 469 postmortem brains, we find that while the ALDH2 rs671 variant is associated with increased plaque deposits and a higher Aβ40/42 ratio, it is not an independent risk factor for Alzheimer's disease. Mechanistically, we show that lower ALDH2 activity leads to 4-HNE accumulation in the brain. The (R)-4-HNE enantiomer adducts to residue Lys53 of C99, favoring Aβ40 generation in the Golgi apparatus. Decreased ALDH2 activity also lowers inflammatory factor secretion, as well as amyloid β phagocytosis and spread in brains of patients with Alzheimer's disease. We thus define the relationship between the ALDH2 rs671 polymorphism and amyloid β pathology, and find that ALDH2 rs671 is a key regulator of Aβ40 or Aβ42 generation.

Comparative transcriptome analysis reveals nicotine metabolism is a critical component for enhancing stress response intensity of innate immunity system in tobacco

The pyridine alkaloid nicotine acts as one of best-studied plant resistant traits in tobacco. Previous research has shown that NtERF199 and NtERF189, acting as master regulators within the NIC1 and NIC2 locus, quantitatively contribute to nicotine accumulation levels in N. tabacum. Genome editing-created Nic1(Nterf199) and Nic2 (Nterf189) double mutant provides an ideal platform for precisely dissecting the defensive role of nicotine and the connection between the nicotine biosynthetic pathway with other putative metabolic networks. Taking this advantage, we performed a comparative transcriptomic analysis to reevaluate the potential physiological and metabolic changes in response to nicotine synthesis defect by comparing the nic1nic2 and NIC1NIC2 plants. Our findings revealed that nicotine reduction could systematically diminishes the expression intensities of genes associated with stimulus perception, signal transduction and regulation, as well as secondary metabolic flux. Consequently, this global expression reduction might compromise tobacco adaptions to environmental fitness, herbivore resistances, and plant growth and development. The up-regulation of a novel set of stress-responsive and metabolic pathway genes might signify a newly established metabolic reprogramming to tradeoff the detrimental effect of nicotine loss. These results offer additional compelling evidence regarding nicotine's critical defensive role in nature and highlights the tight link between nicotine biosynthesis and gene expression levels of quantitative resistance-related genes for better environmental adaptation.

PCDHA9 as a candidate gene for amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a devastating neurodegenerative disease. To identify additional genetic factors, we analyzed exome sequences in a large cohort of Chinese ALS patients and found a homozygous variant (p.L700P) in PCDHA9 in three unrelated patients. We generated Pcdhα9 mutant mice harboring either orthologous point mutation or deletion mutation. These mice develop progressive spinal motor loss, muscle atrophy, and structural/functional abnormalities of the neuromuscular junction, leading to paralysis and early lethality. TDP-43 pathology is detected in the spinal motor neurons of aged mutant mice. Mechanistically, we demonstrate that Pcdha9 mutation causes aberrant activation of FAK and PYK2 in aging spinal cord, and dramatically reduced NKA-α1 expression in motor neurons. Our single nucleus multi-omics analysis reveals disturbed signaling involved in cell adhesion, ion transport, synapse organization, and neuronal survival in aged mutant mice. Together, our results present PCDHA9 as a potential ALS gene and provide insights into its pathogenesis.

Forsythiaside A suppresses renal fibrosis and partial epithelial-mesenchymal transition by targeting THBS1 through the PI3K/AKT signaling pathway

Renal fibrosis is a key feature of chronic kidney disease (CKD) progression, whereas no proven effective anti-fibrotic treatments. Forsythiaside A (FTA), derived from Forsythia suspense, has been found to possess nephroprotective properties. However, there is limited research on its anti-fibrotic effects, and its mechanism of action remains unknown. This study aimed to investigate the suppressive effects of FTA on renal fibrosis and explore the underlying mechanisms. In vitro, we established a HK2 cell model induced by transforming growth factor β1 (TGF-β1), and in vivo, we used a mice model induced by unilateral ureteral obstruction (UUO). CCK-8 assay, qRT-PCR, Western blotting, immunofluorescence, flow cytometry, histological staining, immunohistochemistry, TUNEL assay, RNA transcriptome sequencing, and molecular docking were performed. The results showed that FTA (40 μM or 80 μM) treatment improved cell viability and suppressed TGF-β1-induced fibrotic changes and partial epithelial-mesenchymal transition (EMT). Furthermore, FTA treatment reversed the activation of the PI3K/AKT signaling pathway, and THBS1 was identified as the target gene. We found that THBS1 knockdown suppressed the activation of the PI3K/AKT signaling pathway and reduced the fibrosis and partial EMT-related protein level. Conversely, THBS1 overexpression activated the PI3K/AKT signaling pathway and exacerbated renal fibrosis and partial EMT. In vivo, mice were administered FTA (30 or 60 mg/kg) for 2 weeks, and the results demonstrated that FTA administration significantly mitigated tubular injury, tubulointerstitial fibrosis, partial EMT, and apoptosis. In conclusion, FTA inhibited renal fibrosis and partial EMT by targeting THBS1 and inhibiting activation of the PI3K/AKT pathway.

A genome-wide association study based on the China Kadoorie Biobank identifies genetic associations between snoring and cardiometabolic traits

Despite the high prevalence of snoring in Asia, little is known about the genetic etiology of snoring and its causal relationships with cardiometabolic traits. Based on 100,626 Chinese individuals, a genome-wide association study on snoring was conducted. Four novel loci were identified for snoring traits mapped on SLC25A21, the intergenic region of WDR11 and FGFR, NAA25, ALDH2, and VTI1A, respectively. The novel loci highlighted the roles of structural abnormality of the upper airway and craniofacial region and dysfunction of metabolic and transport systems in the development of snoring. In the two-sample bi-directional Mendelian randomization analysis, higher body mass index, weight, and elevated blood pressure were causal for snoring, and a reverse causal effect was observed between snoring and diastolic blood pressure. Altogether, our results revealed the possible etiology of snoring in China and indicated that managing cardiometabolic health was essential to snoring prevention, and hypertension should be considered among snorers.

Graphene promotes the growth of Vigna angularis by regulating the nitrogen metabolism and photosynthesis

Graphene has promising applications in agriculture and forestry. In the current study, six different concentrations of graphene (0mg/L, 0.01mg/L, 0.10mg/L, 1.00mg/L, 10.00mg/L, and 100.00mg/L) were used to investigate its effect on the growth and development of V. angularis plants in soil culture. The results showed that the group treated with 1.00mg/L graphene (G-1) had significantly increased plant height (19.86%), stem diameter (24.33%), and leaf area (13.69%), compared to the control group (CK). Moreover, all concentrations of graphene had positive effects on the total root length, total root surface area, and the number of root tips of V. angularis. Compared to the CK group, the G-1 group had significantly increased leaf water potential (37.89%), leaf conductivity (2.25%), and SOD, POD, and CAT activities (47.67%, 35.22%, and 199.3%, respectively). The G-1 group also showed improved leaf net photosynthetic rate, chlorophyll content, and soluble sugar content (51.28%, 24.25%, and 38.35%, respectively), compared to the CK group. Additionally, 1.00mg/L graphene led to a 23.88% increase in the podding rate and a 17.04% increase in the yield of V. angularis plants. The rhizosphere soil of V. angularis treated with 1.00mg/L graphene had a 25.14% increase in hydrolyzable nitrogen content and a 66.67% increase in available phosphorus content. RNA-seq data indicated that 1.00mg/L graphene induced the expression of photosynthesis and nitrogen transmembrane transport genes, including ATP synthase subunit b, photosystem I reaction center subunit XI, photosystem I reaction center subunit IV A, ferredoxin, and psbP-like protein 1, as well as genes for photosynthesis antenna proteins, glutamine synthetase, glutamate dehydrogenase 1, cyanate hydratase, protein fluG-like, and NRT1/PTR family, suggesting that graphene promoted the growth and development of V. angularis by enhancing the photosynthesis and nitrogen metabolism processes in V. angularis plants. Our results indicated that a suitable concentration of graphene could significantly promote the growth of V. angularis plants in soil.

Unraveling the signaling roles of indole in an opportunistic pathogen Pseudomonas aeruginosa strain Jade-X

Pseudomonas aeruginosa, which can produce several toxins and form biofilm, is listed among the priority pathogens. Indole is a ubiquitous aromatic pollutant and signaling molecule produced by tryptophanase in bacteria. Herein, the impacts of indole on a newly isolated P. aeruginosa strain Jade-X were systematically investigated. Indole (0.5-2.0 mM) enhanced the biofilm production by 1.33-2.31-fold after 24 h incubation at 30 °C. However, the effects indole on biofilm formation were intricate and closely intertwined with factors such as incubation temperature, bacterial growth stage, and indole concentration. The twitching motility was enhanced by 1.15-1.99-fold by indole, potentially facilitating surface exploration and biofilm development. Indole reduced the production of virulence factors (pyocyanin and pyoverdine) as well as altered the surface properties (zeta potential and hydrophobicity). Transcriptional analysis revealed that indole (1.0 mM) significantly downregulated mexGHI-opmD efflux genes (4.73-6.91-fold) and virulence-related genes (pqs, pch, and pvd clusters, and flagella-related genes), while upregulating pili-related genes in strain Jade-X. The quorum sensing related signal regulators, including RhlR, LasR, and MvfR (PqsR), were not altered by indole, while other six transcriptional regulators (AmrZ, BfmR, PchR, QscR, SoxR, and SphR) were significantly affected, implying that indole effects might be regulated in a complex and delicate manner. This study should provide new insights into our understanding of indole signaling roles.

Plant genomic resources at National Genomics Data Center: assisting in data-driven breeding applications

Genomic data serve as an invaluable resource for unraveling the intricacies of the higher plant systems, including the constituent elements within and among species. Through various efforts in genomic data archiving, integrative analysis and value-added curation, the National Genomics Data Center (NGDC), which is a part of the China National Center for Bioinformation (CNCB), has successfully established and currently maintains a vast amount of database resources. This dedicated initiative of the NGDC facilitates a data-rich ecosystem that greatly strengthens and supports genomic research efforts. Here, we present a comprehensive overview of central repositories dedicated to archiving, presenting, and sharing plant omics data, introduce knowledgebases focused on variants or gene-based functional insights, highlight species-specific multiple omics database resources, and briefly review the online application tools. We intend that this review can be used as a guide map for plant researchers wishing to select effective data resources from the NGDC for their specific areas of study.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-023-00134-4.

Carbapenem-resistant Klebsiella pneumoniae capsular types, antibiotic resistance and virulence factors in China: a longitudinal, multi-centre study

Epidemiological knowledge of circulating carbapenem-resistant Klebsiella pneumoniae (CRKP) is needed to develop effective strategies against this public health threat. Here we present a longitudinal analysis of 1,017 CRKP isolates recovered from patients from 40 hospitals across China between 2016 and 2020. Virulence gene and capsule typing revealed expansion of CRKP capsule type KL64 (59.5%) alongside decreases in KL47 prevalence. Hypervirulent CRKP increased in prevalence from 28.2% in 2016 to 45.7% in 2020. Phylogenetic and spatiotemporal analysis revealed Beijing and Shanghai as transmission hubs accounting for differential geographical prevalence of KL47 and KL64 strains across China. Moderate frequency capsule or O-antigen loss was also detected among isolates. Non-capsular CRKP were more susceptible to phagocytosis, attenuated during mouse infections, but showed increased serum resistance and biofilm formation. These findings give insight into CRKP serotype prevalence and dynamics, revealing the importance of monitoring serotype shifts for the future development of immunological strategies against CRKP infections.

Transcriptome Study of Bursaphelenchus xylophilus Treated with Fomepizole Reveals a Serine/Threonine-Protein Phosphatase Gene that Is Substantially Linked with Vitality and Pathogenicity

Bursaphelenchus xylophilus, the pine wood nematode (PWN), is the causal agent of pine wilt disease (PWD), which causes enormous economic loss annually. According to our previous research, fomepizole, as a selective inhibitor of PWN alcohol dehydrogenase (ADH), has the potential to be a preferable lead compound for developing novel nematicides. However, the underlying molecular mechanism is still unclear. The result of molecular docking showed that the stronger interactions between fomepizole and PWN ADH at the active site of ADH were attributed to hydrogen bonds. Low-dose fomepizole had a substantial negative impact on the egg hatchability, development, oviposition, and lifespan of PWN. Transcriptome analysis indicated that 2,124 upregulated genes and 490 downregulated genes in fomepizole-treated PWN were obtained. Kyoto Encyclopedia of Genes and Genomes enrichment analysis of differentially expressed genes indicated that fomepizole could be involved in controlling PWN vitality mainly by regulating key signaling pathways, such as the ribosome, hippo signaling pathway, and lysosome. Remarkably, the results of RNA interference indicated that the downregulated serine/threonine-protein phosphatase gene (stpp) could reduce the egg hatchability, development, oviposition, and lifespan of PWN, which was closely similar to the consequences of nematodes with low-dose fomepizole treatment. In addition, the silencing of stpp resulted in weakness of PWN pathogenicity, which indicated that stpp could be a potential drug target to control PWN.

Multi-omic profiling of clear cell renal cell carcinoma identifies metabolic reprogramming associated with disease progression

Clear cell renal cell carcinoma (ccRCC) is a complex disease with remarkable immune and metabolic heterogeneity. Here we perform genomic, transcriptomic, proteomic, metabolomic and spatial transcriptomic and metabolomic analyses on 100 patients with ccRCC from the Tongji Hospital RCC (TJ-RCC) cohort. Our analysis identifies four ccRCC subtypes including De-clear cell differentiated (DCCD)-ccRCC, a subtype with distinctive metabolic features. DCCD cancer cells are characterized by fewer lipid droplets, reduced metabolic activity, enhanced nutrient uptake capability and a high proliferation rate, leading to poor prognosis. Using single-cell and spatial trajectory analysis, we demonstrate that DCCD is a common mode of ccRCC progression. Even among stage I patients, DCCD is associated with worse outcomes and higher recurrence rate, suggesting that it cannot be cured by nephrectomy alone. Our study also suggests a treatment strategy based on subtype-specific immune cell infiltration that could guide the clinical management of ccRCC.

Whole-Exome Sequencing Reveals the Genomic Profile and IL6ST Variants as a Prognostic Biomarker of Paraneoplastic Pemphigus-Associated Unicentric Castleman Disease

Unicentric Castleman disease (UCD) is a rare lymphoproliferative disorder. Paraneoplastic pemphigus (PNP) is a major complication associated with poor UCD prognosis. However, the genomic profiles and prognostic biomarkers of PNP-associated UCD remain unclear. In this study, we performed whole-exome sequencing analysis for 28 matched tumor-normal pairs and 9 tumor-only samples to define the genomic landscape of Chinese patients with PNP-associated UCD. An integrative analysis was performed to identify somatic variants, the mutational signatures, and key pathways in tumors. Besides, we analyzed the relationship among mutated genes, clinical characteristics, and prognosis. Sixty-one somatic mutant genes were identified in >1 patient with PNP-associated UCD. Specifically, IL6ST and PDGFRB were the most frequently mutated genes (32%), followed by DPP6 (18%) and MUC4 (18%). Signaling molecules and interactions, cellular processes, and signal transduction pathways were enriched. Furthermore, we found that poor overall survival was related to IL6ST variants (P = .02). Finally, we classified PNP-associated UCD into 4 genomic subgroups: IL6ST, PDGFRB, IL6ST-PDGFRB, and an unknown subgroup. In summary, we defined the molecular profile of PNP-associated UCD and identified a potential molecular biomarker for predicting prognosis, which may provide therapeutic targets for treating this severe disorder.

AEBP1 promotes papillary thyroid cancer progression by activating BMP4 signaling

Papillary thyroid cancer (PTC) is the most prevalent endocrine cancer worldwide. Approximately 30 % of PTC patients will progress into the advanced or metastatic stage and have a relatively poor prognosis. It is well known that epithelial-mesenchymal transition (EMT) plays a pivotal role in thyroid cancer metastasis, resistance to therapy, and recurrence. Clarifying the molecular mechanisms of EMT in PTC progression will help develop the targeted therapy of PTC. The aberrant expression of some transcription factors (TFs) participated in many pathological processes of cancers including EMT. In this study, by performing bioinformatics analysis, adipocyte enhancer-binding protein 1 (AEBP1) was screened as a pivotal TF that promoted EMT and tumor progression in PTC. In vitro experiments indicated that knockout of AEBP1 can inhibit the growth and invasion of PTC cells and reduce the expression of EMT markers including N-cadherin, TWIST1, and ZEB2. In the xenograft model, knockout of AEBP1 inhibited the growth and lung metastasis of PTC cells. By performing RNA-sequencing, dual-luciferase reporter assay, and chromatin immunoprecipitation assay, Bone morphogenetic protein 4 (BMP4) was identified as a downstream target of AEBP1. Over-expression of BMP4 can rescue the inhibitory effects of AEBP1 knockout on the growth, invasion, and EMT phenotype of PTC cells. In conclusion, these findings demonstrated that AEBP1 plays a critical role in PTC progression by regulating BMP4 expression and the AEBP1-BMP4 axis may present novel therapeutic targets for PTC treatment.

The Hypothalamic Epigenetic Landscape in Dietary Obesity

The hypothalamus in the brain plays a pivotal role in controlling energy balance in vertebrates. Nutritional excess through high-fat diet (HFD) feeding can dysregulate hypothalamic signaling at multiple levels. Yet, it remains largely unknown in what magnitude HFD feeding may impact epigenetics in this brain region. Here, it is shown that HFD feeding can significantly alter hypothalamic epigenetic events, including posttranslational histone modifications, DNA methylation, and chromatin accessibility. The authors comprehensively analyze the chromatin immunoprecipitation-sequencing (ChIP-seq), methylated DNA immunoprecipitation-sequencing (MeDIP-seq), single nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq), and RNA-seq data of the hypothalamus of C57 BL/6 mice fed with a chow or HFD for 1 to 6 months. The chromatins are categorized into 6 states using the obtained ChIP-seq data for H3K4me3, H3K27ac, H3K9me3, H3K27me3, and H3K36me3. A 1-month HFD feeding dysregulates histone modifications and DNA methylation more pronouncedly than that of 3- or 6-month. Besides, HFD feeding differentially impacts chromatin accessibility in hypothalamic cells. Thus, the epigenetic landscape is dysregulated in the hypothalamus of dietary obesity mice.

Genomic insight into the origin, domestication, dispersal, diversification and human selection of Tartary buckwheat

BACKGROUND

Tartary buckwheat, Fagopyrum tataricum, is a pseudocereal crop with worldwide distribution and high nutritional value. However, the origin and domestication history of this crop remain to be elucidated.

RESULTS

Here, by analyzing the population genomics of 567 accessions collected worldwide and reviewing historical documents, we find that Tartary buckwheat originated in the Himalayan region and then spread southwest possibly along with the migration of the Yi people, a minority in Southwestern China that has a long history of planting Tartary buckwheat. Along with the expansion of the Mongol Empire, Tartary buckwheat dispersed to Europe and ultimately to the rest of the world. The different natural growth environments resulted in adaptation, especially significant differences in salt tolerance between northern and southern Chinese Tartary buckwheat populations. By scanning for selective sweeps and using a genome-wide association study, we identify genes responsible for Tartary buckwheat domestication and differentiation, which we then experimentally validate. Comparative genomics and QTL analysis further shed light on the genetic foundation of the easily dehulled trait in a particular variety that was artificially selected by the Wa people, a minority group in Southwestern China known for cultivating Tartary buckwheat specifically for steaming as a staple food to prevent lysine deficiency.

CONCLUSIONS

This study provides both comprehensive insights into the origin and domestication of, and a foundation for molecular breeding for, Tartary buckwheat.

Sequencing of cerebrospinal fluid cell-free DNA facilitated early differential diagnosis of intramedullary spinal cord tumors

Pre-surgery differential diagnosis is valuable for personalized treatment planning in intramedullary spinal cord tumors. This study assessed the performance of sequencing cell-free DNA (cfDNA) in cerebrospinal fluid (CSF) for differential diagnosis of these tumors. Prospectively enrolling 45 patients with intramedullary spinal cord lesions, including diffuse midline glioma (DMG), H3K27-altered (14/45), glioblastoma (1/45), H3-wildtype-astrocytoma (10/45), ependymoma (11/45), and other lesions (9/45), CSF samples were collected via lumbar puncture (41/45), intraoperative extraction (3/45), and Ommaya reservoir (1/45). Then, these samples underwent targeted sequencing along with paired tissue DNA. DMG, H3K27-altered patients exhibited a higher ctDNA positivity (85.7%, 12/14) compared to patients with H3-wildtype-astrocytoma (0/8, P = 0.0003), ependymoma (2/10, P = 0.003), and glioneuronal tumor (0/3, P = 0.009). The histological-grade-IV (P = 0.0027), Ki-67 index ≥10% (P = 0.014), and tumor reaching spinal cord surface (P = 0.012) are also associated with higher ctDNA positivity. Interestingly, for patients with TERT promoter mutant tumors, TERT mutation was detectable in the CSF cfDNA of one DMG case, but not other five cases with histological-grade-II tumors. Shared copy number variants were exclusively observed in DMG, H3K27-altered, and showed a strong correlation (Correlation = 0.95) between CSF and tissue. Finally, H3K27M mutations in CSF exhibited high diagnostic efficiency for DMG, H3K27-altered (Sensitivity = 85.7%, Specificity = 100.0%, AUC = 0.929). Notably, H3K27M was detectable in CSF from patients with recurrent tumors, making it easily applicable for postoperative monitoring. In conclusion, the molecular profile from ctDNA released into CSF of malignant tumors was more frequently detected compared to relatively benign ones. Sequencing of ctDNA in CSF exhibited high efficiency for the differential diagnosis of DMG, H3K27-altered.

Metabolomics and metagenomics reveal the impact of γδ T inhibition on gut microbiota and metabolism in periodontitis-promoting OSCC

During the process of periodontitis-promoting oral squamous cell carcinoma (OSCC), the periodontitis microbiota can facilitate OSCC development by activating γδ T cells. Inhibiting γδ T cells through immunotherapy has been shown to significantly alleviate various types of cancer. However, the underlying mechanism by which inhibiting γδ T cells influenced cancer treatment has not been fully elucidated. In this study, a mouse model of OSCC with periodontitis was established, and γδ T cells were inhibited by antibodies. Gut samples from the mice were collected and analyzed by metabolomics, metagenomics, and 16S rRNA. Integrative analysis of the gut metabolome and microbiome revealed that targeting γδ T resulted in changes in the levels of metabolites associated with cancer in the gut. Although there was no difference in the α diversity of microbiota, β diversity was significantly different, with a more heterogeneous community structure in the mice receiving targeted γδ T immunotherapy. Statistical analysis of the gut microbiota at the species level revealed a significant enrichment of Lactobacillus murinus, which was significantly correlated with γδ T abundance. The functional analysis revealed that inhibiting γδ T could impact the functional gene. A comprehensive analysis revealed that L. murinus is especially associated with changes in adenine, which also had connection with the concentration of IL-17 and the abundance of γδ T.IMPORTANCEThis study revealed the effect of γδ T cells on gut microbial dysbiosis and identify potential links between gut microbiota and metabolism, providing new insights into the role of γδ T during the process of periodontitis-induced OSCC, and identifying relevant biomarkers for future research and clinical monitoring protocol development.

Analysis of Antibiotic Resistance Genes in Water Reservoirs and Related Wastewater from Animal Farms in Central China

This study aimed to explore the phenotype and relationship of drug resistance genes in livestock and poultry farm wastewater and drinking water reservoirs to provide evidence for the transmission mechanisms of drug resistance genes, in order to reveal the spread of drug resistance genes in wastewater from intensive farms in Central China to urban reservoirs that serve as drinking water sources and provide preliminary data for the treatment of wastewater from animal farms to reduce the threat to human beings. DNA extraction and metagenomic sequencing were performed on eight groups of samples collected from four water reservoirs and four related wastewaters from animal farms in Central China. Metagenomic sequencing showed that the top 20 AROs with the highest abundance were vanT_gene, vanY_gene, adeF, qacG, Mtub_rpsL_STR, vanY_gene_, vanW_gene, Mtub_murA_FOF, vanY_gene, vanH_gene, FosG, rsmA, qacJ, RbpA, vanW_gene, aadA6, vanY_gene, sul4, sul1, and InuF. The resistance genes mentioned above belong to the following categories of drug resistance mechanisms: antibiotic target replacement, antibiotic target protection, antibiotic inactivation, and antibiotic efflux. The resistomes that match the top 20 genes are Streptococcus agalactiae and Streptococcus anginosus; Enterococcus faecalis; Enterococcus faecium; Actinomyces viscosus and Bacillus cereus. Enterococcus faecium; Clostridium tetani; Streptococcus agalactiae and Streptococcus anginosus; Streptococcus agalactiae and Streptococcus anginosus; Acinetobacter baumannii, Bifidobacterium bifidum, Bifidobacterium breve, Bifidobacterium longum, Corynebacterium jeikeium, Corynebacterium urealyticum, Mycobacterium kansasii, Mycobacterium tuberculosis, Schaalia odontolytica, and Trueperella pyogenes; Mycobacterium avium and Mycobacterium tuberculosis; Aeromonas caviae, Enterobacter hormaechei, Vibrio cholerae, Vibrio metoecus, Vibrio parahaemolyticus, and Vibrio vulnificus; Pseudomonas aeruginosa and Pseudomonas fluorescens; Staphylococcus aureus and Staphylococcus equorum; M. avium, Achromobacter xylosoxidans, and Acinetobacter baumannii; Sphingobium yanoikuyae, Acinetobacter indicus, Morganella morganii, Proteus mirabilis, Proteus vulgaris, Providencia rettgeri, and Providencia stuartii. Unreported drug resistance genes and drug-resistant bacteria in Central China were identified in 2023. In the transmission path of drug resistance genes, the transmission path from aquaculture wastewater to human drinking water sources cannot be ignored. For the sake of human health and ecological balance, the treatment of aquaculture wastewater needs to be further strengthened, and the effective blocking of drug resistance gene transmission needs to be considered.

A landscape of gene expression regulation for synovium in arthritis

The synovium is an important component of any synovial joint and is the major target tissue of inflammatory arthritis. However, the multi-omics landscape of synovium required for functional inference is absent from large-scale resources. Here we integrate genomics with transcriptomics and chromatin accessibility features of human synovium in up to 245 arthritic patients, to characterize the landscape of genetic regulation on gene expression and the regulatory mechanisms mediating arthritic diseases predisposition. We identify 4765 independent primary and 616 secondary cis-expression quantitative trait loci (cis-eQTLs) in the synovium and find that the eQTLs with multiple independent signals have stronger effects and heritability than single independent eQTLs. Integration of genome-wide association studies (GWASs) and eQTLs identifies 84 arthritis related genes, revealing 38 novel genes which have not been reported by previous studies using eQTL data from the GTEx project or immune cells. We further develop a method called eQTac to identify variants that could affect gene expression by affecting chromatin accessibility and identify 1517 regions with potential regulatory function of chromatin accessibility. Altogether, our study provides a comprehensive synovium multi-omics resource for arthritic diseases and gains new insights into the regulation of gene expression.

Transcriptome analysis reveals cadmium exposure enhanced the isoquinoline alkaloid biosynthesis and disease resistance in Coptis chinensis

Coptis chinensis Franch is a perennial herb from the Ranunculaceae family with a long history of medicinal use. As the medicinal part, the rhizome of coptis often accumulates excessive cadmium (Cd) even at low concentrations in the soil, which not only compromises its medicinal safety but also raises concerns about adverse effects on human health. Therefore, effective strategies are needed to mitigate this accumulation and ensure its safe use in traditional medicine. This study utilized transcriptome profiling and physiological analysis to explore molecular mechanisms associated with ecological significance and the active accumulation of Cd in C. chinensis. The response to Cd in C. chinensis was assessed through RNA sequencing, Cd determination and isoquinoline alkaloid measurement using its roots, stems, and leaves. The transcriptome revealed, a total of 2667, 2998, or 2815 up-regulated deferentially expressed genes in roots, stems or leaves in response to Cd exposure. Furthermore, we identified phenylpropanoid and isoquinoline alkaloid biosynthesis as the key pathways response to Cd exposure, which suggests that C. chinensis may improve its tolerance to Cd through regulating the phenylpropanoid biosynthesis pathway. Under Cd exposure, plant-pathogen interaction in leaves was identified as the key pathway, which indicates that upregulation of genes involved in plant-pathogen interaction could enhance disease resistance in C. chinensis. WGCNA analysis identified WRKY8 (Cluster-55763.31419) and WRKY47 (Cluster-55763.221590) as potential regulators of secondary metabolic synthesis and plant-pathogen interaction pathway in C. chinensis triggered by Cd. The measurement of berberine, coptisine, palmatine, and epiberberine also demonstrated that Cd simulated the four isoquinoline alkaloids in roots. Therefore, our study not only presented a transcriptome expression profiles that revealed significant upregulation of genes involved in metal transport and detoxification pathways but also suggested a possible mechanism to cope with Cd accumulation. This knowledge provides a new insight into gene manipulation for controlling Cd accumulation, enhancing resistance and promoting synthesis of secondary metabolites with potential medicinal properties in other medicinal plant species.

Single-cell sequencing identifies inflammation-promoting fibroblast-neutrophil interaction in peri-implantitis

AIM

To reveal the cellular composition and molecular environment of the periodontal and peri-implant inflammatory infiltrates through a single-cell sequencing technique, which may explain the pathological difference between these two diseases. A special focus was placed on the phenotypes and potential roles of neutrophils and fibroblasts in peri-implant/periodontal tissue immunity.

MATERIALS AND METHODS

High-throughput single-cell transcriptomic profiling of peri-implant tissues from patients with peri-implantitis as well as periodontal tissues from patients with periodontitis and healthy donors was performed. Immunofluorescence analysis was carried out to further validate the identified cell subtypes and their involvement in peri-implantitis and periodontitis.

RESULTS

Based on our single-cell resolution analysis, a quantified proportional increase of neutrophil (Neu) subtypes was shown in peri-implantitis. Among these, a predominance of Neutro_CXCR2 was revealed. We also found the involvement of inflammation-promoting fibroblasts as well as a predominance of CXCL8+ fibroblast-CXCR2+ neutrophil interaction in peri-implantitis.

CONCLUSIONS

Our study indicated that the predominance of CXCL8+ fibroblast-CXCR2+ neutrophil interaction might underline the enhanced host response in peri-implantitis compared with periodontitis. This information offers a molecular basis by which fibroblast and neutrophil subtypes might be diagnostically and therapeutically targeted in peri-implantitis.

Genomic epidemiology of hypervirulent carbapenem-resistant Klebsiella pneumoniae at Jinshan local hospital, Shanghai, during 2014-2018

BACKGROUND

Hypervirulent carbapenem-resistant Klebsiella pneumoniae (Hv-CRKP) triggered a significant public health challenge. This study explored the prevalence trends and key genetic characteristics of Hv-CRKP in one Shanghai suburbs hospital during 2014-2018.

METHODS

During five years, Hv-CRKP strains identified from 2579 CRKP by specific PCR, were subjected to performed short- and long-read sequencing technology; epidemiological characteristics, antimicrobial-resistance genes (ARGs), virulence determinants, detailed plasmid profiles and conjugation efficiency were comprehensively investigated.

RESULTS

155 Hv-CRKP and 31 non-Hv-CRKP strains were sequenced. Hv-CRKP strains exhibited significant resistance to six common antibiotic classes (>92%). ST11 steadily increased and became the most prevalent ST (85.2%), followed by ST15 (8.5%), ST65 (2.6%), ST23 (1.9%), and ST86 (0.6%). ST11-KL64 (65.2%) rapidly increased from 0 in 2014 to 93.9% in 2018. blaKPC-2 was the primary carbapenemase gene (97.4%). Other ARGs switched from aac(3)-IId to aadA2 in aminoglycoside and from sul1 to sul2 in sulfanilamide. The time-dated phylogenetic tree was divided into four independent evolutionary clades. Clade 1 and 3 strains were mostly limited in the ICU, whereas Clade 2 strains were distributed among multiple departments. Compared to ybt14 in ICEKp12 in Clade 1, Clade 3 strains harbored ybt9 in ICEKp3 and blaCTX-M-65. Hv-CRKP infected more wards than non-Hv-CRKP and showed greater transmission capacity. Three plasmids containing crucial carbapenemase genes demonstrated their early transmission across China.

CONCLUSION

The Hv-CRKP ST11-KL64 has rapidly replaced ST11-KL47 and emerged as the predominant epidemic subtype in various hospital wards, highlighting the importance of conducting comprehensive early surveillance for Hv-CRKP, especially in respiratory infections.

Spatial transcriptomics reveals heterogeneity of histological subtypes between lepidic and acinar lung adenocarcinoma

BACKGROUND

Patients who possess various histological subtypes of early-stage lung adenocarcinoma (LUAD) have considerably diverse prognoses. The simultaneous existence of several histological subtypes reduces the clinical accuracy of the diagnosis and prognosis of early-stage LUAD due to intratumour intricacy.

METHODS

We included 11 postoperative LUAD patients pathologically confirmed to be stage IA. Single-cell RNA sequencing (scRNA-seq) was carried out on matched tumour and normal tissue. Three formalin-fixed and paraffin-embedded cases were randomly selected for 10× Genomics Visium analysis, one of which was analysed by digital spatial profiler (DSP).

RESULTS

Using DSP and 10× Genomics Visium analysis, signature gene profiles for lepidic and acinar histological subtypes were acquired. The percentage of histological subtypes predicted for the patients from samples of 11 LUAD fresh tissues by scRNA-seq showed a degree of concordance with the clinicopathologic findings assessed by visual examination. DSP proteomics and 10× Genomics Visium transcriptomics analyses revealed that a negative correlation (Spearman correlation analysis: r = -.886; p = .033) between the expression levels of CD8 and the expression trend of programmed cell death 1(PD-L1) on tumour endothelial cells. The percentage of CD8+ T cells in the acinar region was lower than in the lepidic region.

CONCLUSIONS

These findings illustrate that assessing patient histological subtypes at the single-cell level is feasible. Additionally, tumour endothelial cells that express PD-L1 in stage IA LUAD suppress immune-responsive CD8+ T cells.

Pluripotency state transition of embryonic stem cells requires the turnover of histone chaperone FACT on chromatin

The differentiation of embryonic stem cells (ESCs) begins with the transition from the naive to the primed state. The formative state was recently established as a critical intermediate between the two states. Here, we demonstrate the role of the histone chaperone FACT in regulating the naive-to-formative transition. We found that the Q265K mutation in the FACT subunit SSRP1 increased the binding of FACT to histone H3-H4, impaired nucleosome disassembly in vitro, and reduced the turnover of FACT on chromatin in vivo. Strikingly, mouse ESCs harboring this mutation showed elevated naive-to-formative transition. Mechanistically, the SSRP1-Q265K mutation enriched FACT at the enhancers of formative-specific genes to increase targeted gene expression. Together, these findings suggest that the turnover of FACT on chromatin is crucial for regulating the enhancers of formative-specific genes, thereby mediating the naive-to-formative transition. This study highlights the significance of FACT in fine-tuning cell fate transition during early development.

An Inflammatory Bowel Diseases Integrated Resources Portal (IBDIRP)

IBD, including ulcerative colitis and Crohn's disease, is a chronic and debilitating gastrointestinal disorder that affects millions of people worldwide. Research on IBD has generated massive amounts of data, including literature, metagenomics, metabolomics, bioresources and databases. We aim to create an IBD Integrated Resources Portal (IBDIRP) that provides the most comprehensive resources for IBD. An integrated platform was developed that provides information on different aspects of IBD research resources, such as single-nucleotide polymorphisms (SNPs), genes, transcriptome, microbiota, metabolomics, single cells and other resources. Valuable and comprehensive IBD-related data were collected from PubMed, Google, GMrepo, gutMega, gutMDisorder, Single Cell Portal and other sources. Then, the data were systematically sorted, and these resources were manually curated. We systematically sorted and cataloged more than 320 unique risk SNPs associated with IBD in the SNP section. We presented over 289 IBD-related genes based on the database collection in the gene section. We also obtained 153 manually curated IBD transcriptomics data, including 12 388 samples, on the Gene Expression Omnibus database. The sorted IBD-related microbiota data from three primary microbiome databases (GMrepo, gutMega and gutMDisorder) were available for download. We selected 23 149 IBD-related taxonomic records from these databases. Additionally, we collected 24 IBD metabolomics studies with 2896 participants in the metabolomics section. We introduced two interactive single-cell data plug-in units that provided data visualization based on cells and genes. Finally, we listed 18 significant IBD web resources, such as the official European Crohn's and Colitis Organisation and International Organization for the Study of IBD websites, IBD scoring tools, IBD genetic and multi-omics resources, IBD biobanks and other useful research resources. The IBDIRP website is the first integrated resource for global IBD researchers. This portal will help researchers by providing comprehensive knowledge and enabling them to reinforce the multidimensional impression of IBD. The IBDIRP website is accessible via www.ibdirp.com Database URL: www.ibdirp.com.

Genetic associations with longevity are on average stronger in females than in males

It is long observed that females tend to live longer than males in nearly every country. However, the underlying mechanism remains elusive. In this study, we discovered that genetic associations with longevity are on average stronger in females than in males through bio-demographic analyses of genome-wide association studies (GWAS) dataset of 2178 centenarians and 2299 middle-age controls of Chinese Longitudinal Healthy Longevity Study (CLHLS). This discovery is replicated across North and South regions of China, and is further confirmed by North-South discovery/replication analyses of different and independent datasets of Chinese healthy aging candidate genes with CLHLS participants who are not in CLHLS GWAS, including 2972 centenarians and 1992 middle-age controls. Our polygenic risk score analyses of eight exclusive groups of sex-specific genes, analyses of sex-specific and not-sex-specific individual genes, and Genome-wide Complex Trait Analysis using all SNPs all reconfirm that genetic associations with longevity are on average stronger in females than in males. Our discovery/replication analyses are based on genetic datasets of in total 5150 centenarians and compatible middle-age controls, which comprises the worldwide largest sample of centenarians. The present study's findings may partially explain the well-known male-female health-survival paradox and suggest that genetic variants may be associated with different reactions between males and females to the same vaccine, drug treatment and/or nutritional intervention. Thus, our findings provide evidence to steer away from traditional view that "one-size-fits-all" for clinical interventions, and to consider sex differences for improving healthcare efficiency. We suggest future investigations focusing on effects of interactions between sex-specific genetic variants and environment on longevity as well as biological function.

Cerebrospinal fluid ctDNA testing shows an advantage over plasma ctDNA testing in advanced non-small cell lung cancer patients with brain metastases

Background

Brain metastases (BM), including brain parenchyma metastases (BPM) and leptomeningeal metastases (LM), are devastating metastatic complications in advanced cancer patients. Next-generation sequencing (NGS) is emerging as a new promising tool for profiling cancer mutation, which could facilitate the diagnosis of cancer. This retrospective study aimed to investigate the molecular genetic characteristics of non-small cell lung cancer (NSCLC) patients with BPM and LM using NGS.

Methods

Cerebrospinal fluid (CSF) samples and paired plasma samples were collected from 37 patients of NSCLC-BM. We profiled genetic mutation characteristics using NGS from NSCLC-BM by comparing CSF circulating tumour DNA (ctDNA) with plasma ctDNA and primary tumour tissues.

Results

Among the 37 patients with NSCLC-BM, 28 patients had LM with or without BPM, while 9 patients only had BPM. Driver and drug-resistant mutations in primary tumours with LM included: EGFR L858R (10, 35.7%), EGFR 19del (6, 21.4%), EGFR L858R+MET (1, 3.6%), EGFR L858R+S768I (1, 3.6%), ALK (2, 7.1%), ROS1 (1, 3.6%), negative (5, 17.9%), and unknown (2, 7.1%). In patients with NSCLC-LM, the detection rate and abundance of ctDNA in the CSF were significantly higher than those in paired plasma. The main driver mutations of NSCLC-LM remained highly consistent with those of the primary tumours, along with other unique mutations. Circulating tumour DNA was negative in the CSF samples of BPM patients. Patients with BMP had a higher ratio of EGFR 19del than L858R mutation (55.6% vs 11.1.%), whereas NSCLC patients with LM had a higher ratio of EGFR L858R than 19del mutation (50.0% vs 25.0%). Most patients with positive plasma ctDNA results were male (p = 0.058) and in an unstable state (p = 0.003).

Conclusion

Our study indicated that the CSF ctDNA detected by NGS may reflect the molecular characteristics and heterogeneity of NSCLC-LM. Timely screening of patients with NSCLC for CSF ctDNA, especially for patients with positive plasma ctDNA, may facilitate the early detection of LM. Furthermore, patients with the EGFR 19del may have a higher risk of developing BPM.

Multiomics profiling reveals the benefits of gamma-delta (γδ) T lymphocytes for improving the tumor microenvironment, immunotherapy efficacy and prognosis in cervical cancer

BACKGROUND

As an unconventional subpopulation of T lymphocytes, γδ T cells can recognize antigens independently of major histocompatibility complex restrictions. Recent studies have indicated that γδ T cells play contrasting roles in tumor microenvironments-promoting tumor progression in some cancers (eg, gallbladder and leukemia) while suppressing it in others (eg, lung and gastric). γδ T cells are mainly enriched in peripheral mucosal tissues. As the cervix is a mucosa-rich tissue, the role of γδ T cells in cervical cancer warrants further investigation.

METHODS

We employed a multiomics strategy that integrated abundant data from single-cell and bulk transcriptome sequencing, whole exome sequencing, genotyping array, immunohistochemistry, and MRI.

RESULTS

Heterogeneity was observed in the level of γδ T-cell infiltration in cervical cancer tissues, mainly associated with the tumor somatic mutational landscape. Definitely, γδ T cells play a beneficial role in the prognosis of patients with cervical cancer. First, γδ T cells exert direct cytotoxic effects in the tumor microenvironment of cervical cancer through the dynamic evolution of cellular states at both poles. Second, higher levels of γδ T-cell infiltration also shape the microenvironment of immune activation with cancer-suppressive properties. We found that these intricate features can be observed by MRI-based radiomics models to non-invasively assess γδ T-cell proportions in tumor tissues in patients. Importantly, patients with high infiltration levels of γδ T cells may be more amenable to immunotherapies including immune checkpoint inhibitors and autologous tumor-infiltrating lymphocyte therapies, than to chemoradiotherapy.

CONCLUSIONS

γδ T cells play a beneficial role in antitumor immunity in cervical cancer. The abundance of γδ T cells in cervical cancerous tissue is associated with higher response rates to immunotherapy.

Integration of RRBS and RNA-seq unravels the regulatory role of DNMT3A in porcine Sertoli cell proliferation

DNMT3A participates in de novo methylation, yet its impact on the proliferation of testicular Sertoli cells remains unclear. Development-specific methylation has been proven to be associated with cellular development. Therefore, in this study, we simulated DNMT3A expression pattern during testicular development by DNMT3A interference. Then, RRBS and RNA-seq were used to decipher DNMT3A regulatory mechanisms on Sertoli cell proliferation. Immunofluorescence staining revealed the expression of DNMT3A in the Sertoli cells of the prepubertal testis. DNMT3A was demonstrated to inhibit the cell cycle and proliferation of Sertoli cells, while promoting cell apoptosis. After transfected with DNMT3A interference, a total of 560 DEGs and 2,091 DMGs produced by DNMT3A interference were identified between two treated groups, respectively. Integrating the results from RRBS and RNA-seq, the overlapping genes between DMGs and DEGs were found to be enriched in the Gene Ontology (GO) terms related to cellular development and the Apelin signaling pathway. The present study demonstrated the impact of DNMT3A on the proliferation of porcine testicular Sertoli cells, suggesting that DNMT3A primarily acts through the Apelin signaling pathway. These findings provide valuable insights into how DNMT3A influences testicular development and health, offering new perspectives.

Genetic dissection of major QTL for grain number per spike on chromosomes 5A and 6A in bread wheat (Triticum aestivum L.)

Grain number per spike (GNS) is a crucial component of grain yield and plays a significant role in improving wheat yield. To identify quantitative trait loci (QTL) associated with GNS, a recombinant inbred line (RIL) population derived from the cross of Zhongkemai 13F10 and Chuanmai 42 was employed to conduct QTL mapping across eight environments. Based on the bulked segregant exome sequencing (BSE-Seq), genomic regions associated with GNS were detected on chromosomes 5A and 6A. According to the constructed genetic maps, two major QTL QGns.cib-5A (LOD = 4.35-8.16, PVE = 8.46-14.43%) and QGns.cib-6A (LOD = 3.82-30.80, PVE = 5.44-12.38%) were detected in five and four environments, respectively. QGns.cib-6A is a QTL cluster for other seven yield-related traits. QGns.cib-5A and QGns.cib-6A were further validated using linked Kompetitive Allele Specific PCR (KASP) markers in different genetic backgrounds. QGns.cib-5A exhibited pleiotropic effects on productive tiller number (PTN), spike length (SL), fertile spikelet number per spike (FSN), and ratio of grain length to grain width (GL/GW) but did not significantly affect thousand grain weight (TGW). Haplotype analysis revealed that QGns.cib-5A and QGns.cib-6A were the targets of artificial selection during wheat improvement. Candidate genes for QGns.cib-5A and QGns.cib-6A were predicted by analyzing gene annotation, spatiotemporal expression patterns, and orthologous and sequence differences. These findings will be valuable for fine mapping and map-based cloning of genes underlying QGns.cib-5A and QGns.cib-6A.

Co-expression network analysis reveals PbTGA4 and PbAPRR2 as core transcription factors of drought response in an important timber species Phoebe bournei

Phoebe bournei is one of the main afforestation tree species in subtropical regions of China and is famous for its timber. Its distribution and growth are significantly impaired by water conditions. Thus, it is essential to understand the mechanism of the stress response in P. bournei. Here, we analyzed the phenotypic changes and transcriptomic rearrangement in the leaves and roots of P. bournei seedlings grown for 0 h, 1 h, 24 h, and 72 h under simulated drought conditions (10% PEG 6000). The results showed that drought stress inhibited plant photosynthesis and increased oxidoreductase activity and abscisic acid (ABA) accumulation. Spatio-temporal transcriptomic analysis identified 2836 and 3704 differentially expressed genes (DEGs) in leaves and roots, respectively. The responsive genes in different organs presented various expression profiles at different times. Gene co-expression network analysis identified two core transcription factors, TGA4 and APRR2, from two modules that showed a strong positive correlation with ABA accumulation. Our study investigated the different responses of aboveground and belowground organs of P. bournei to drought stress and provides critical information for improving the drought resistance of this timber species.

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

The National Genomics Data Center (NGDC), which is a part of the China National Center for Bioinformation (CNCB), provides a family of database resources to support the global academic and industrial communities. With the rapid accumulation of multi-omics data at an unprecedented pace, CNCB-NGDC continuously expands and updates core database resources through big data archiving, integrative analysis and value-added curation. Importantly, NGDC collaborates closely with major international databases and initiatives to ensure seamless data exchange and interoperability. Over the past year, significant efforts have been dedicated to integrating diverse omics data, synthesizing expanding knowledge, developing new resources, and upgrading major existing resources. Particularly, several database resources are newly developed for the biodiversity of protists (P10K), bacteria (NTM-DB, MPA) as well as plant (PPGR, SoyOmics, PlantPan) and disease/trait association (CROST, HervD Atlas, HALL, MACdb, BioKA, BioKA, RePoS, PGG.SV, NAFLDkb). All the resources and services are publicly accessible at https://ngdc.cncb.ac.cn.

Phylogenomic analyses and chromosome ploidy identification reveal multiple cryptic species in Allium sikkimense complex (Amaryllidaceae)

Polyploidization is a process that typically leads to instantaneous reproductive isolation and has, therefore, been considered as one of the major evolutionary forces in the species-rich Hengduan Mountains (HM), yet this topic remains poorly studied in the region. Allium sikkimense and its relatives (about eight species) compose a natural diploid-polyploid complex with the highest diversity in the HM and adjacent areas. A combination of nuclear ribosomal DNA (nrDNA), plastome, transcriptome, and ploidy identification through chromosome counting and flow cytometry is employed to reconstruct the phylogenetic relationships in this complex and to investigate the frequency and the evolutionary significance of polyploidy in the complex. The plastome failed to resolve the phylogenetic relationships of the different species in the A. sikkimense complex, and the phylogenetic tree based on nrDNA also has limited resolution. However, our study reveals a well-resolved phylogenetic framework for species in the A. sikkimense complex using more than 1,000 orthologous genes from the transcriptome data. Previously recognized morphospecies A. sikkimense are non-monophyletic and comprise at least two independently evolved lineages (i.e., cryptic species), each forming a clade with different diploid species in this complex. The embedded pattern of octoploid A. jichouense and tetraploid A. sp. nov. within different polyploid samples of A. sikkimense supports a possible scenario of budding speciation (via niche divergence). Furthermore, our results reveal that co-occurring species in the A. sikkimense complex usually have different ploidy levels, suggesting that polyploidy is an important process for reproductive isolation of sympatric Allium species. Phylogenetic network analyses suggested that the phylogenetic relationships of the A. sikkimense complex, allowing for reticulation events, always fit the dataset better than a simple bifurcating tree. In addition, the included or exserted filaments, which have long been used to delimit species, are highly unreliable taxonomically due to their extensive parallel and convergent evolution.

Hepatoid adenocarcinoma in ureter with next-generation sequencing: A case report and literature review

BACKGROUND

Hepatoid adenocarcinoma (HAC) is rare in the urinary system, with only 7 reported cases in upper urinary tract. This report aimed to explore the genetic characteristics of ureteral HAC for first time, and to describe the treatment prognosis of ureteral HAC.

CASE PRESENTATION

We present a rare case of ureteral HAC in a 53-year-old female, showing elevated serum levels of AFP and CEA, prolonged chronic irritation may be an important cause of her ureteral HAC. Radical nephroureterectomy was performed, the serum levels of AFP and CEA decreased significantly, and metastasis in lymph nodes was found at 9 months after surgery, she had no related symptoms after 18 months postoperatively without adjuvant chemotherapy. Three driver somatic mutations in cancer were identified by NGS testing, including: TP53D281H, KMT2DL1211Ifs*2, KMT2DT1843Nfs*5, demonstrating that ureteral HAC has the similar mutational features to upper tract urothelial carcinoma. Homologous-recombination deficiency (HRD) was positive in this tumor with no mutations in HRD-related genes, which was possibly induced by the copy number deletion of SETD2 gene.

CONCLUSIONS

We report a rare case of ureteral HAC with elevated serum levels of AFP and CEA. NGS testing demonstrated that ureteral HAC has the similar mutational features to upper tract urothelial carcinoma, which is an important guide for the diagnosis and treatment of ureteral HAC.

Transcriptome analysis of axillary buds in low phosphorus stress and functional analysis of TaWRKY74s in wheat

BACKGROUND

Wheat is one of the main grain crops in the world, and the tiller number is a key factor affecting the yield of wheat. Phosphorus is an essential element for tiller development in wheat. However, due to decreasing phosphorus content in soil, there has been increasing use of phosphorus fertilizer, while imposing risk of soil and water pollution. Hence, it is important to identify low phosphorus tolerance genes and utilize them for stress resistance breeding in wheat.

RESULTS

We subjected the wheat variety Kenong 199 (KN199) to low phosphorus stress and observed a reduced tiller number. Using transcriptome analysis, we identified 1651 upregulated genes and 827 downregulated of genes after low phosphorus stress. The differentially expressed genes were found to be enriched in the enzyme activity regulation related to phosphorus, hormone signal transduction, and ion transmembrane transport. Furthermore, the transcription factor analysis revealed that TaWRKY74s were important for low phosphorus tolerance. TaWRKY74s have three alleles: TaWRKY74-A, TaWRKY74-B, and TaWRKY74-D, and they all belong to the WRKY family with conserved WRKYGQK motifs. These proteins were found to be located in the nucleus, and they were expressed in axillary meristem, shoot apical meristem(SAM), young leaves, leaf primordium, and spikelet primordium. The evolutionary tree showed that TaWRKY74s were closely related to OsWRKY74s in rice. Moreover, TaWRKY74s-RNAi transgenic plants displayed significantly fewer tillers compared to wild-type plants under normal conditions. Additionally, the tiller numebr of the RNAi transgenic plants was also significantly lower than that of the wild-type plants under low-phosphorus stress, and increased the decrease amplitude. This suggestd that TaWRKY74s are related to phosphorus response and can affect the tiller number of wheat.

CONCLUSIONS

The results of this research showed that TaWRKY74s were key genes in wheat response to low phosphorus stress, which might regulate wheat tiller number through abscisic acid (ABA) and auxin signal transduction pathways. This research lays the foundation for further investigating the mechanism of TaWRKY74s in the low phosphorus environments and is significant for wheat stress resistance breeding.

Genotype and Phenotype of Renal Hypouricemia: A Single-Center Study from China

BACKGROUND

Renal hypouricemia (RHUC), a rare inherited disorder characterized by impaired uric acid reabsorption and subsequent profound hypouricemia, occurs mainly due to variants in SLC22A12 or SLC2A9. Only anecdotal cases and one small-scale RHUC screening study have been reported in the Chinese population.

METHODS

A total of 19 patients with RHUC from 17 unrelated families were recruited from our center. The medical history, clinical manifestations, biochemical exam, and clinical outcomes were collected. Next-generation sequencing-based targeted gene sequencing or whole exon sequencing was performed.

RESULTS

A total of 22 variants in SLC22A12 or SLC2A9 were found in 19 patients. The variant c.944G>A (p.W315X) in SLC2A9 was identified in three patients. Three variants c.165C>A (p.D55E), c.1549_1555delGAGACCC (p.E517Rfs*17), and c.1483T>C (p.W495R) in SLC22A12 and three variants c.1215+1G>A (splicing variant), c.643A>C (p.T215P), and c.227C>A (p.S76X) in SLC2A9 were novel. A proportion of 10 out of 19 patients presented with exercise-induced acute kidney injury (EIAKI). The renal outcome was favorable. Five patients had nephrolithiasis, in whom three had hypercalciuria.

CONCLUSION

The current study reported six novel variants in SLC22A12 and SLC2A9 genes of Chinese patients with RHUC. The variant c.944G>A (p.W315X) in SLC2A9 may be common in Chinese patients. EIAKI is the main clinical phenotype associated with RHUC in our cohort, with a favorable outcome. Hypercalciuria presented in some RHUC patients is a new finding.

Seasonal expressions of nerve growth factor (NGF), and its receptor TrkA and p75 in the scent glands of muskrats (Ondatra zibethicus)

NGF, also known as nerve growth factor, is crucial for the survival and differentiation of the nervous system, in addition to being involved in a number of non-neuronal systems. The aim of this work was to investigate the immunolocalization and expression patterns of NGF, its receptor, tyrosine kinase receptor A (TrkA), and p75 in the scent glands of muskrats (Ondatra zibethicus) throughout the breeding and non-breeding seasons. The scent gland mass showed considerable seasonal variations, with higher values during the breeding season and comparatively lower levels during the non-breeding season. While no immunostaining was observed in the interstitial cells, NGF, TrkA, and p75 were immunolocalized in the scent glandular cells and epithelial cells during both breeding and non-breeding seasons. NGF, TrkA, and p75 protein and mRNA expression levels were higher in the scent glands during breeding season compared to the non-breeding season. Circulating levels of follicle-stimulating hormone (FSH), luteinizing hormone (LH), testosterone (T), and T in the scent gland were all significantly higher throughout the breeding season. The relative levels of the hormones in the plasma and the scent glands as well as NGF, TrkA, and p75 were positively associated with each other. Additionally, transcriptome analysis of the scent glands revealed that differentially expressed genes may be linked to steroid biosynthesis, the estrogen signaling pathway, and neurotransmitter transmembrane transporter function. These results suggest a potential role for NGF, TrkA, and p75 in controlling seasonal variations in the muskrats' scent gland functioning.

METTL4-mediated N6-methyladenine DNA modification regulates thermotolerance in Arabidopsis thaliana

DNA N6-methyladenine (6 mA) is an evolutionarily conserved DNA modification in procaryotes and eukaryotes. The DNA 6 mA methylation is tightly controlled by 6 mA regulatory proteins. DNA N6-adenine methyltransferase 1 (DAMT-1) has been identified as a DNA 6 mA methyltransferase in animals. In plants, DNA 6 mA methylation has been found, however, the DNA 6 mA methyltransferases and their function in plants are largely unknown. In our study, we find METTL4 is a DNA 6 mA methyltransferase in Arabidopsis thaliana. Both in vitro and in vivo evidences support the DNA 6 mA methyltransferase activity of METTL4. mettl4 mutant is hypersensitive to heat stress, suggesting DNA 6 mA methylation plays important role in heat stress adaption. RNA-seq and 6 mA IP-qPCR analysis show that METTL4 participates in heat stress tolerance by regulating expression of heat responsive genes. Our study find METTL4 is a plant DNA 6 mA methyltransferase and illustrates its function in regulating heat stress response.

Chromatin remodeling analysis reveals the RdDM pathway responds to low-phosphorus stress in maize

Chromatin in eukaryotes folds into a complex three-dimensional (3D) structure that is essential for controlling gene expression and cellular function and is dynamically regulated in biological processes. Studies on plant phosphorus signaling have concentrated on single genes and gene interactions. It is critical to expand the existing signaling pathway in terms of its 3D structure. In this study, low-Pi treatment led to greater chromatin volume. Furthermore, low-Pi stress increased the insulation score and the number of TAD-like domains, but the effects on the A/B compartment were not obvious. The methylation levels of target sites (hereafter as RdDM levels) peaked at specific TAD-like boundaries, whereas RdDM peak levels at conserved TAD-like boundaries shifted and decreased sharply. The distribution pattern of RdDM sites originating from the Helitron transposons matched that of genome-wide RdDM sites near TAD-like boundaries. RdDM pathway genes were upregulated in the middle or early stages and downregulated in the later stages under low-Pi conditions. The RdDM pathway mutant ddm1a showed increased tolerance to low-Pi stress, with shortened and thickened roots contributing to higher Pi uptake from the shallow soil layer. ChIP-seq results revealed that ZmDDM1A could bind to Pi- and root development-related genes. Strong associations were found between interacting genes in significantly different chromatin-interaction regions and root traits. These findings not only expand the mechanisms by which plants respond to low-Pi stress through the RdDM pathway but also offer a crucial framework for the analysis of biological issues using 3D genomics.

A cytosine analogue 5-azacitidine improves the accumulation of licochalcone A in licorice Glycyrrhiza inflata

Licochalcone A (LCA) is a characteristic compound of Glycyrrhiza inflata with anti-inflammatory, antioxidant and antitumor activities. However, G. inflata produces LCA in low quantities that does not meet the market demand. In this study, we found that DNA methylation inhibitor 5-azacitidine (5-azaC) successfully improved the LCA contents in G. inflata seedlings. Transcriptome analysis revealed a series of differentially expressed genes (DEGs), including transcription factors such as MYB, ERF, WRKY, and some structural genes related to flavonoid biosynthesis. However, whole genome bisulfite sequencing (BS-seq) results showed little effect of the 5-azaC treatment on the alteration of DNA methylation on these genes, indicating the possibility that 5-azaC acts as a stimulus, but not an epigenetic modulation factor to improve the LCA content in G. inflata. Additionally, we applied the 5-azaC treatment to field plants and hairy roots and successfully increased the LCA contents in both cases. This research demonstrates the feasibility of 5-azaC treatments in future applications to improve plant production of LCA.

Discovery of Taurocholic Acid Sodium Hydrate as a Novel Repurposing Drug for Intervertebral Disc Degeneration by Targeting MAPK3

OBJECTIVE

Nowadays, more than 90% of people over 50 years suffer from intervertebral disc degeneration (IDD), but there are exist no ideal drugs. The aim of this study is to identify a new drug for IDD.

METHODS

An approved small molecular drug library including 2040 small molecular compounds was used here. We found that taurocholic acid sodium hydrate (NAT) could induce chondrogenesis and osteogenesis in mesenchymal stem cells (MSCs). Then, an in vivo mouse model of IDD was established and the coccygeal discs transcriptome analysis and surface plasmon resonance analysis (SPR) integrated with liquid chromatography-tandem mass spectrometry assay (LC-MS) were performed in this study to study the therapy effect and target proteins of NAT for IDD. Micro-CT was used to evaluate the cancellous bone. The expression of osteogenic (OCN, RNX2), chondrogenic (COL2A1, SOX9), and the target related (ERK1/2, p-ERK1/2) proteins were detected. The alkaline phosphatase staining was performed to estimate osteogenic differentiation. Blood routine and blood biochemistry indexes were analyzed for the safety of NAT.

RESULTS

The results showed that NAT could induce chondrogenesis and osteogenesis in MSCs. Further experiments confirmed NAT could ameliorate the secondary osteoporosis and delay the development of IDD in mice. Transcriptome analysis identified 128 common genes and eight Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways for NAT. SPR-LC-MS assay detected 57 target proteins for NAT, including MAPK3 (mitogen-activated protein kinase 3), also known as ERK1 (extracellular regulated protein kinase 1). Further verification experiment confirmed that NAT significantly reduced the expression of ERK1/2 phosphorylation.

CONCLUSION

NAT would induce chondrogenesis and osteogenesis of MSCs, ameliorate the secondary osteoporosis and delay the progression of IDD in mice by targeting MAPK3.Furthermore, MAPK3, especially the phosphorylation of MAPK3, would be a potential therapeutic target for IDD treatment.

Leveraging Temporal Wnt Signal for Efficient Differentiation of Intestinal Stem Cells in an Organoid Model

The homeostasis of the intestinal epithelium heavily relies on the self-renewal and differentiation of intestinal stem cells (ISCs). Although the orchestration of these processes by signaling pathways such as the Wnt, BMP, Notch, and MAPK signals has been extensively studied, the dynamics of their regulation remains unclear. Our study explores how the Wnt signaling pathway temporally regulates the differentiation of ISCs into various cell types in an intestinal organoid system. We report that the duration of Wnt exposure following Notch pathway inactivation significantly influences the differentiation direction of intestinal epithelial cells toward multiple secretory cell types, including goblet cells, enteroendocrine cells (EECs), and Paneth cells. This temporal regulation of Wnt signaling adds another layer of complexity to the combination of niche signals that govern cell fate. By manipulating this temporal signal, we have developed optimized protocols for the efficient in vitro differentiation of ISCs into EECs and goblet cells. These findings provide critical insights into the dynamic regulation of ISC differentiation and offer a robust platform for future investigations into intestinal biology and potential therapeutic applications.

Single-cell RNA sequencing uncovers the cell type-dependent transcriptomic changes in the retrosplenial cortex after peripheral nerve injury

The retrosplenial cortex (RSC) is a vital area for storing remote memory and has recently been found to undergo broad changes after peripheral nerve injury. However, little is known about the role of RSC in pain regulation. Here, we examine the involvement of RSC in the pain of mice with nerve injury. Notably, reducing the activities of calcium-/calmodulin-dependent protein kinase type II-positive splenial neurons chemogenetically increases paw withdrawal threshold and extends thermal withdrawal latency in mice with nerve injury. The single-cell or single-nucleus RNA-sequencing results predict enhanced excitatory synaptic transmissions in RSC induced by nerve injury. Local infusion of 1-naphthyl acetyl spermine into RSC to decrease the excitatory synaptic transmissions relieves pain and induces conditioned place preference. Our data indicate that RSC is critical for regulating physiological and neuropathic pain. The cell type-dependent transcriptomic information would help understand the molecular basis of neuropathic pain.

Molecule characterization of chemosensory and metabolism-related genes in the proboscis of Athetis lepigone

Introduction: The moth species Athetis lepigone (Möschler) (Lepidoptera: Noctuidae), which has recently been identified as a pest of summer maize (Zea mays L.) in China, has demonstrated a rapid proliferation with in the Huang-Huai-Hai Plain region since its initial discovery in Hebei Province in 2005. It has become a prevalent pest of corn crops, and its ability to adapt quickly to its surroundings is currently being investigated. One of the key characteristics of its siphoning mouthparts is not only the feeding apparatus itself but also the chemosensory organs that enable the detection of chemical signals from the surrounding environment. However, there is a lack of comprehensive research on the genes responsible for chemosensory and metabolic mechanisms in the proboscises of male and female A. lepigone adults. Methods: In this study, we utilized transcriptome analysis to identify a total of fifty chemosensory genes from six distinct families, including 19 odorant-binding proteins (OBPs), 22 chemosensory proteins (CSPs), one co-receptor (Orco), six odorant receptors (ORs), four ionotropic receptors (IRs), and two sensory neuron membrane proteins (SNMPs) in the proboscis. Notably, seven OBPs, two CSPs, and one OR were discovered for the first time. Additionally, fourteen genes related to metabolism, including cytochrome P450 (CYPs) and carboxylesterases (CXEs), were also identified. Furthermore, a qualitative analysis was conducted on the relative transcript levels of eight related genes. The expression of 21 annotated chemosensory and metabolic genes was compared between A. lepigone adults and larvae using qRT-PCR, revealing tissue specificity. The majority of genes exhibited predominant expression in the antennae and proboscis during the adult stage, while showing slight expression in the combination of sixth-instar larval head oral appendages (maxilla, labium, and antenna) and pheromone gland-ovipositors of female adults. Results/discussion: Our study points to a new pest control strategies that these newly discovered genes have the potential to serve as targets for enhancing future pest control, including mating disruption and the use of food attractants. And it would be advantageous to ascertain the distribution of chemosensory gene expression and gain insights into the functionalities of these genes, thereby establishing a novel theoretical framework for the advancement of eco-friendly pesticides and efficient pest management strategies in the future.