Single-cell transcriptome analysis of the in vivo response to viral infection in the cave nectar bat Eonycteris spelaea.

Bats are reservoir hosts of many zoonotic viruses with pandemic potential. We utilized single-cell transcriptome sequencing (scRNA-seq) to analyze the immune response in bat lungs upon in vivo infection with a double-stranded RNA virus, Pteropine orthoreovirus PRV3M. Bat neutrophils were distinguished by high basal IDO1 expression. NK cells and T cells were the most abundant immune cells in lung tissue. Three distinct CD8+ effector T cell populations could be delineated by differential expression of KLRB1, GFRA2, and DPP4. Select NK and T clusters increased expression of genes involved in T cell activation and effector function early after viral infection. Alveolar macrophages and classical monocytes drove antiviral interferon signaling. Infection expanded a CSF1R+ population expressing collagen-like genes, which became the predominant myeloid cell type post-infection. This work uncovers features relevant to viral disease tolerance in bats, lays a foundation for future experimental work, and serves as a resource for comparative immunology studies.

Loss of PTEN Is Associated with Resistance to Anti-PD-1 Checkpoint Blockade Therapy in Metastatic Uterine Leiomyosarcoma.

Response to immune checkpoint blockade in mesenchymal tumors is poorly characterized, but immunogenomic dissection of these cancers could inform immunotherapy mediators. We identified a treatment-naive patient who has metastatic uterine leiomyosarcoma and has experienced complete tumor remission for >2 years on anti-PD-1 (pembrolizumab) monotherapy. We analyzed the primary tumor, the sole treatment-resistant metastasis, and germline tissue to explore mechanisms of immunotherapy sensitivity and resistance. Both tumors stained diffusely for PD-L2 and showed sparse PD-L1 staining. PD-1+ cell infiltration significantly decreased in the resistant tumor (p = 0.039). Genomically, the treatment-resistant tumor uniquely harbored biallelic PTEN loss and had reduced expression of two neoantigens that demonstrated strong immunoreactivity with patient T cells in vitro, suggesting long-lasting immunological memory. In this near-complete response to PD-1 blockade in a mesenchymal tumor, we identified PTEN mutations and reduced expression of genes encoding neoantigens as potential mediators of resistance to immune checkpoint therapy.

The correlation between TRIM28 expression and immune checkpoints in CRPC.

This study delves into the unexplored realm of castration-resistant prostate cancer (CRPC) by investigating the role of TRIM28 and its intricate molecular mechanisms using high-throughput single-cell transcriptome sequencing and advanced bioinformatics analysis. Our comprehensive examination unveiled dynamic TRIM28 expression changes, particularly in immune cells such as macrophages and CD8+ T cells within CRPC. Correlation analyses with TCGA data highlighted the connection between TRIM28 and immune checkpoint expression and emphasized its pivotal influence on the quantity and functionality of immune cells. Using TRIM28 knockout mouse models, we identified differentially expressed genes and enriched pathways, unraveling the potential regulatory involvement of TRIM28 in the cGAS-STING pathway. In vitro, experiments further illuminated that TRIM28 knockout in prostate cancer cells induced a notable anti-tumor immune effect by inhibiting M2 macrophage polarization and enhancing CD8+ T cell activity. This impactful discovery was validated in an in situ transplant tumor model, where TRIM28 knockout exhibited a deceleration in tumor growth, reduced proportions of M2 macrophages, and enhanced infiltration of CD8+ T cells. In summary, this study elucidates the hitherto unknown anti-tumor immune role of TRIM28 in CRPC and unravels its potential regulatory mechanism via the cGAS-STING signaling pathway. These findings provide novel insights into the immune landscape of CRPC, offering promising directions for developing innovative therapeutic strategies.

Differential mRNA profiles reveal the potential roles of genes involved in lactate stimulation in mouse macrophages.

Lactate is a glycolysis end product, and its levels are markedly associated with disease severity, morbidity, and mortality in sepsis. It modulates key functions of immune cells, including macrophages. In this investigation, transcriptomic analysis was performed using lactic acid, sodium lactate, and hydrochloric acid-stimulated mouse bone marrow-derived macrophages (iBMDM), respectively, to identify lactate-associated signaling pathways. After 24 h of stimulation, 896 differentially expressed genes (DEG) indicated were up-regulation, whereas 792 were down-regulated in the lactic acid group, in the sodium lactate group, 128 DEG were up-regulated, and 41 were down-regulated, and in the hydrochloric acid group, 499 DEG were up-regulated, and 285 were down-regulated. Subsequently, clinical samples were used to further verify the eight genes with significant differences, among which Tssk6, Ypel4, Elovl3, Trp53inp1, and Cfp were differentially expressed in patients with high lactic acid, indicating their possible involvement in lactic acid-induced inflammation and various physiological diseases caused by sepsis. However, elongation of very long chain fatty acids protein 3 (Elovl3) was negatively correlated with lactic acid content in patients. The results of this study provide a necessary reference for better understanding the transcriptomic changes caused by lactic acid and explain the potential role of high lactic acid in the regulation of macrophages in sepsis.

Whole transcriptome sequencing indicated the Anti-tumor immunity of NLRP3 in breast cancer.

Breast cancer (BC) is a prevalent cancer of the female reproductive system and a major contributor to cancer-related mortality. The activation of NLRP3, a key inflammasome, has been extensively associated with tumor-related molecular and cellular processes; however, the regulatory mechanisms and specific role of NLRP3 in breast cancer remain incompletely elucidated. This study aimed to evaluate the molecular mechanisms of NLRP3-related genes in BC. Utilizing bioinformatics methods, the present research analyzed the TCGA-BRCA dataset, which included four groups of transcriptome sequencing data as follows, normal (WT), NLRP3 knockout (KO), non-knockout-BRCA (BC-WT), and NLRP3-knockout-BRCA (BC-KO). Results indicated that NLRP3 was significantly down-regulated in TCGA-BRCA. Key module genes were mainly enriched in leukocyte cell-cell adhesion and cytokine-cytokine receptor interaction. Moreover, correlation analysis showed that NLRP3 was positively associated with cancer-associated fibroblasts and negatively associated with CD4+ Th1 T-cells. In addition, the DEGs1 and DEGs2 overlapping indicated 505 feature genes, with Chac1 (negative) and Ugt8a (positive) had the strongest correlation with differential immune cells (class-switched memory B cells). Pathway intersection revealed 13 co-KEGG pathways. The BC-KO group indicated markedly reduced levels of four genes (Ccl19, Ccl20, Ccl21a, and H2-Oa) and increased levels of two genes (Il2ra and H2-Ob). This study delved into the role of NLRP3 in BC, exploring its regulatory mechanisms and the impact gene knockout. Bioinformatics approaches identified NLRP3-associated genes, their enriched pathways, and interactions within the tumor microenvironment (TME), providing novel insights into NLRP3 function, TME dynamics, and potential targets for BC prevention and treatment.

Novel MED15::ATF1 fusion in a pediatric melanoma with spitzoid features and aggressive presentation.

Childhood melanoma is a rare and biologically heterogeneous pediatric malignancy. The differential diagnosis of pediatric melanoma is usually broad, including a wide variety of spindle cell or epithelioid neoplasms. Different molecular alterations affecting the MAPK and PI3K/AKT/mTOR pathways, tumor suppressor genes, and telomerase reactivation have been implicated in melanoma tumorigenesis and progression. Here, we report a novel MED15::ATF1 fusion in a pediatric melanoma with spitzoid features and an aggressive clinical course.

Impact of LuxS on virulence and pathogenicity in Klebsiella pneumoniae exhibiting varied mucoid phenotypes.

How the LuxS/AI-2 quorum sensing (QS) system influences the pathogenicity of K. pneumoniae is complicated by the heterogeneity of the bacterial mucoid phenotypes. This study aims to explore the LuxS-mediated regulation of the pathogenicity of K. pneumoniae with diverse mucoid phenotypes, including hypermucoid, regular-mucoid, and nonmucoid. The wild-type, luxS knockout, and complemented strains of three K. pneumoniae clinical isolates with distinct mucoid phenotypes were constructed. The results revealed the downregulation of virulence genes of regular-mucoid, and nonmucoid but not hypermucoid strains. The deletion of luxS reduced the pathogenicity of the regular-mucoid, and nonmucoid strains in mice; while in hypermucoid strain, luxS knockout reduced virulence in late growth but enhanced virulence in the early growth phase. Furthermore, the absence of luxS led the regular-mucoid and nonmucoid strains to be more sensitive to the host cell defense, and less biofilm-productive than the wild-type at both the low and high-density growth state. Nevertheless, luxS knockout enhanced the resistances to adhesion and phagocytosis by macrophage as well as serum-killing, of hypermucoid K. pneumoniae at its early low-density growth state, while it was opposite to those in its late high-density growth phase. Collectively, our results suggested that LuxS plays a crucial role in the pathogenicity of K. pneumoniae, and it is highly relevant to the mucoid phenotypes and growth phases of the strains. LuxS probably depresses the capsule in the early low-density phase and promotes the capsule, biofilm, and pathogenicity during the late high-density phase, but inhibits lipopolysaccharide throughout the growth phase, in K. pneumoniae.IMPORTANCECharacterizing the regulation of physiological functions by the LuxS/AI-2 quorum sensing (QS) system in Klebsiella pneumoniae strains will improve our understanding of this important pathogen. The genetic heterogeneity of K. pneumoniae isolates complicates our understanding of its pathogenicity, and the association of LuxS with bacterial pathogenicity has remained poorly addressed in K. pneumoniae. Our results demonstrated strain and growth phase-dependent variation in the contributions of LuxS to the virulence and pathogenicity of K. pneumoniae. Our findings provide new insights into the important contribution of the LuxS/AI-2 QS system to the networks that regulate the pathogenicity of K. pneumoniae. Our study will facilitate our understanding of the regulatory mechanisms of LuxS/AI-2 QS on the pathogenicity of K. pneumoniae under the background of their genetic heterogeneity and help develop new strategies for diminished bacterial virulence within the clinical K. pneumoniae population.

Downregulation of HbFPS1 affects rubber biosynthesis of Hevea brasiliensis suffering from tapping panel dryness.

Tapping panel dryness (TPD) is a century-old problem that has plagued the natural rubber production of Hevea brasiliensis. TPD may result from self-protective mechanisms of H. brasiliensis in response to stresses such as excessive hormone stimulation and mechanical wounding (bark tapping). It has been hypothesized that TPD impairs rubber biosynthesis; however, the underlying mechanisms remain poorly understood. In the present study, we firstly verified that TPD-affected rubber trees exhibited lower rubber biosynthesis activity and greater rubber molecular weight compared to healthy rubber trees. We then demonstrated that HbFPS1, a key gene of rubber biosynthesis, and its expression products were downregulated in the latex of TPD-affected rubber trees, as revealed by transcriptome sequencing and iTRAQ-based proteome analysis. We further discovered that the farnesyl diphosphate synthase HbFPS1 could be recruited to small rubber particles by HbSRPP1 through protein-protein interactions to catalyze farnesyl diphosphate (FPP) synthesis and facilitate rubber biosynthesis initiation. FPP content in the latex of TPD-affected rubber trees was significantly decreased with the downregulation of HbFPS1, ultimately resulting in abnormal development of rubber particles, decreased rubber biosynthesis activity, and increased rubber molecular weight. Upstream regulator assays indicated that a novel regulator, MYB2-like, may be an important regulator of downregulation of HbFPS1 in the latex of TPD-affected rubber trees. Our findings not only provide new directions for studying the molecular events involved in rubber biosynthesis and TPD syndrome and contribute to rubber management strategies, but also broaden our knowledge of plant isoprenoid metabolism and its regulatory networks.

Multiple transcriptome analyses reveal mouse testis developmental dynamics.

The testes are the organs of gamete production and testosterone synthesis. Up to date, no model system is available for mammalian testicular development, and only few studies have characterized the mouse testis transcriptome from no more than three postnatal ages. To describe the transcriptome landscape of the developing mouse testis and identify the potential molecular mechanisms underlying testis maturation, we examined multiple RNA-seq data of mouse testes from 3-week-old (puberty) to 11-week-old (adult). Sperm cells appeared as expected in 5-week-old mouse testis, suggesting the proper sample collection. The principal components analysis revealed the genes from 3w to 4w clustered away from other timepoints, indicating they may be the important nodes for testicular development. The pairwise comparisons at two adjacent timepoints identified 7,612 differentially expressed genes (DEGs), resulting in 58 unique mRNA expression patterns. Enrichment analysis identified functions in tissue morphogenesis (3-4w), regulation of peptidase activity (4-5w), spermatogenesis (7-8w), and antigen processing (10-11w), suggesting distinct functions in different developmental periods. 50 hub genes and 10 gene cluster modules were identified in the testis maturation process by protein-protein interaction (PPI) network analysis, and the miRNA-lncRNA-mRNA, miRNA-circRNA-mRNA and miRNA-circRNA-lncRNA-mRNA competing endogenous RNA (ceRNA) networks were constructed. The results suggest that testis maturation is a complex developmental process modulated by various molecules, and that some potential RNA-RNA interactions may be involved in specific developmental stages. In summary, this study provides an update on the molecular basis of testis development, which may help to understand the molecular mechanisms of mouse testis development and provide guidance for mouse reproduction.

Integrated bioinformatics analysis of retinal ischemia/reperfusion injury in rats with potential key genes.

The tissue damage caused by transient ischemic injury is an essential component of the pathogenesis of retinal ischemia, which mainly hinges on the degree and duration of interruption of the blood supply and the subsequent damage caused by tissue reperfusion. Some research indicated that the retinal injury induced by ischemia-reperfusion (I/R) was related to reperfusion time.In this study, we screened the differentially expressed circRNAs, lncRNAs, and mRNAs between the control and model group and at different reperfusion time (24h, 72h, and 7d) with the aid of whole transcriptome sequencing technology, and the trend changes in time-varying mRNA, lncRNA, circRNA were obtained by chronological analysis. Then, candidate circRNAs, lncRNAs, and mRNAs were obtained as the intersection of differentially expression genes and trend change genes. Importance scores of the genes selected the key genes whose expression changed with the increase of reperfusion time. Also, the characteristic differentially expressed genes specific to the reperfusion time were analyzed, key genes specific to reperfusion time were selected to show the change in biological process with the increase of reperfusion time.As a result, 316 candidate mRNAs, 137 candidate lncRNAs, and 31 candidate circRNAs were obtained by the intersection of differentially expressed mRNAs, lncRNAs, and circRNAs with trend mRNAs, trend lncRNAs and trend circRNAs, 5 key genes (Cd74, RT1-Da, RT1-CE5, RT1-Bb, RT1-DOa) were selected by importance scores of the genes. The result of GSEA showed that key genes were found to play vital roles in antigen processing and presentation, regulation of the actin cytoskeleton, and the ribosome. A network included 4 key genes (Cd74, RT1-Da, RT1-Bb, RT1-DOa), 34 miRNAs and 48 lncRNAs, and 81 regulatory relationship axes, and a network included 4 key genes (Cd74, RT1-Da, RT1-Bb, RT1-DOa), 9 miRNAs and 3 circRNAs (circRNA_10572, circRNA_03219, circRNA_11359) and 12 regulatory relationship axes were constructed, the subcellular location, transcription factors, signaling network, targeted drugs and relationship to eye diseases of key genes were predicted. 1370 characteristic differentially expressed mRNAs (spec_24h mRNA), 558 characteristic differentially expressed mRNAs (spec_72h mRNA), and 92 characteristic differentially expressed mRNAs (spec_7d mRNA) were found, and their key genes and regulation networks were analyzed.In summary, we screened the differentially expressed circRNAs, lncRNAs, and mRNAs between the control and model groups and at different reperfusion time (24h, 72h, and 7d). 5 key genes, Cd74, RT1-Da, RT1-CE5, RT1-Bb, RT1-DOa, were selected. Key genes specific to reperfusion time were selected to show the change in biological process with the increased reperfusion time. These results provided theoretical support and a reference basis for the clinical treatment.

Comparative transcriptomics reveals a highly polymorphic Xanthomonas HrpG virulence regulon.

BACKGROUND: Bacteria of the genus Xanthomonas cause economically significant diseases in various crops. Their virulence is dependent on the translocation of type III effectors (T3Es) into plant cells by the type III secretion system (T3SS), a process regulated by the master response regulator HrpG. Although HrpG has been studied for over two decades, its regulon across diverse Xanthomonas species, particularly beyond type III secretion, remains understudied. RESULTS: In this study, we conducted transcriptome sequencing to explore the HrpG regulons of 17 Xanthomonas strains, encompassing six species and nine pathovars, each exhibiting distinct host and tissue specificities. We employed constitutive expression of plasmid-borne hrpG*, which encodes a constitutively active form of HrpG, to induce the regulon. Our findings reveal substantial inter- and intra-specific diversity in the HrpG* regulons across the strains. Besides 21 genes directly involved in the biosynthesis of the T3SS, the core HrpG* regulon is limited to only five additional genes encoding the transcriptional activator HrpX, the two T3E proteins XopR and XopL, a major facility superfamily (MFS) transporter, and the phosphatase PhoC. Interestingly, genes involved in chemotaxis and genes encoding enzymes with carbohydrate-active and proteolytic activities are variably regulated by HrpG*. CONCLUSIONS: The diversity in the HrpG* regulon suggests that HrpG-dependent virulence in Xanthomonas might be achieved through several distinct strain-specific strategies, potentially reflecting adaptation to diverse ecological niches. These findings enhance our understanding of the complex role of HrpG in regulating various virulence and adaptive pathways, extending beyond T3Es and the T3SS.

LncRNA and mRNA expression characteristic and bioinformatic analysis in myocardium of diabetic cardiomyopathy mice.

BACKGROUND: Diabetic cardiomyopathy (DCM) is becoming a very well-known clinical entity and leads to increased heart failure in diabetic patients. Long non-coding RNAs (LncRNAs) play an important role in the pathogenesis of DCM. In the present study, the expression profiles of lncRNAs and mRNAs were illuminated in myocardium from DCM mice, with purpose of exploring probable pathological processes of DCM involved by differentially expressed genes in order to provide a new direction for the future researches of DCM. RESULTS: The results showed that a total of 93 differentially expressed lncRNA transcripts and 881 mRNA transcripts were aberrantly expressed in db/db mice compared with the controls. The top 6 differentially expressed lncRNAs like up-regulated Hmga1b, Gm8909, Gm50252 and down-regulated Msantd4, 4933413J09Rik, Gm41414 have not yet been reported in DCM. The lncRNAs-mRNAs co-expression network analysis showed that LncRNA 2610507I01Rik, 2310015A16Rik, Gm10503, A930015D03Rik and Gm48483 were the most relevant to differentially expressed mRNAs. CONCLUSION: Our results showed that db/db DCM mice exist differentially expressed lncRNAs and mRNAs in hearts. These differentially expressed lncRNAs may be involved in the pathological process of cardiomyocyte apoptosis and fibrosis in DCM.

Endothelial pyroptosis-driven microglial activation in choroid plexus mediates neuronal apoptosis in hemorrhagic stroke rats.

BACKGROUND: Spontaneous intracerebral hemorrhage (ICH) is associated with alarmingly high rates of disability and mortality, and current therapeutic options are suboptimal. A critical component of ICH pathology is the initiation of a robust inflammatory response, often termed "cytokine storm," which amplifies the secondary brain injury following the initial hemorrhagic insult. The precise sources and consequences of this cytokine-driven inflammation are not fully elucidated, necessitating further investigation. METHODS: To address this knowledge gap, our study conducted a comprehensive cytokine profiling using Luminex® assays, assessing 23 key cytokines. We then employed single-cell RNA sequencing and spatial transcriptomics at three critical time points post-ICH: the hyperacute, acute, and subacute phases. Integrating these multimodal analyses allowed us to identify the cellular origins of cytokines and elucidate their mechanisms of action. RESULTS: Luminex® cytokine assays revealed a significant upregulation of IL-6 and IL-1ß levels at the 24-h post-ICH time point. Through the integration of scRNA-seq and spatial transcriptomics in the hemorrhagic hemisphere of rats, we observed a pronounced activation of cytokine-related signaling pathways within the choroid plexus. Initially, immune cell presence was sparse, but it surged 24 h post-ICH, particularly in the choroid plexus, indicating a substantial shift in the immune microenvironment. We traced the source of IL-1ß and IL-6 to endothelial cells, establishing a link to pyroptosis. Endothelial pyroptosis post-ICH induced the production of IL-1ß and IL-6, which activated microglial polarization characterized by elevated expression of Msr1, Lcn2, and Spp1 via the NF-κB pathway in the choroid plexus. Furthermore, we identified neuronal populations undergoing apoptosis, mediated by the Lcn2-SLC22A17 pathway in response to IL-1ß and IL-6 signaling. Notably, the inhibition of pyroptosis using VX-765 significantly mitigated neurological impairments. CONCLUSIONS: Our study provides evidence that endothelial pyroptosis, characterized by the release of IL-1ß and IL-6, triggers microglial polarization through NF-κB pathway activation, ultimately leading to microglia-mediated neuronal apoptosis in the choroid plexus post-ICH. These findings suggest that targeted therapeutic strategies aimed at mitigating endothelial cell pyroptosis and neutralizing inflammatory cytokines may offer neuroprotection for both microglia and neurons, presenting a promising avenue for ICH treatment.

Phase II trial of imatinib mesylate in patients with PDGFRA/B-negative hypereosinophilic syndrome.

The role of imatinib in PDGFRA/B-negative hypereosinophilic syndromes (HES) is controversial because of the heterogeneity of HES and the scarcity of prospective studies. We conducted a phase II clinical trial to evaluate the efficacy of imatinib in PDGFRA/B-negative HES. Thirty-two patients were treated with imatinib (100-400 mg daily), and the molecular basis of their response was identified using whole-exome sequencing (WES) and whole-transcriptome sequencing (WTS). The haematological response rate was 46.9%, with a complete haematological response (CHR) rate of 18.8%. The median time to response was 1.5 months. Among the six patients who achieved CHR, five maintained it until the 24th cycle of imatinib and one lost response after 20 months. The median progression-free survival was 4.3 months. WES and WTS were conducted for 11 patients. The number of non-silent mutations did not differ between responders and non-responders. Nine differentially expressed genes, including SNORD15A, were downregulated in responders. STAT5B::RARA, PAK2::PIGX, and FIP1L1::CHIC2 fusions were identified in patients with sustained responses, and RNF130::BRAF and WNK1::KDM5A fusions were identified in non-responders. Imatinib, along with an appropriate biomarker, could be a promising option for PDGFRA/B-negative HES.

Breadth versus depth: whole transcriptome sequencing has reduced sensitivity for detection of clinically relevant fusions compared to RNA comprehensive genomic profiling.

While there is great potential for unbiased next-generation sequencing (NGS) approaches-eg, whole transcriptome sequencing (WTS)-for exploration, discovery, and clinical application in the realm of oncology, there are limitations that should be considered when relying on these methodologies for clinical decision making. When using WTS for the detection of clinically relevant gene fusions in tumor specimens, a key consideration is whether a limited coverage depth (approximately 30-50X) is sufficient for detecting these events, especially in samples with low tumor purity. We demonstrate the reduced sensitivity of both a commercial WTS assay for the detection of clinically relevant fusions in analytical validation control samples and of a research use only (RUO) WTS assay for the detection of clinically relevant fusions in real-world clinical samples compared to RNA comprehensive genomic profiling (CGP). Notably, the RUO WTS assay would not have reported 30% (6/20) of fusions detected using RNA CGP assays in fusion-positive tumor samples, highlighting a potential disadvantage of broader sequencing.

The profile of cytokines against bacterial infection in dental pulp.

BACKGROUND: Dental pulp in a confined environment, with little connection to the outside and only a small distribution of immune cells, provides a good research model for investigating how cells respond to bacterial infections through cytokines. METHODS: The data of single-cell transcriptome sequencing of healthy and inflamed pulp tissue were downloaded from the GEO dataset. The expression character of 79 cytokines was analyzed based on the expression matrix. RESULTS: The cytokine secretion profiles of the two populations of pulp cells in healthy dental pulp were associated with vascularization and nervous system development, as well as immune cell regulation. For the three populations of pulp stem cells with stem cell activity in the dental pulp, the secretion of cytokines related to nervous system development, regulation of endothelial cell proliferation and migration, and regulation of immune cell function comprised the characteristics that we observed. The cytokines secreted by T cells and macrophages were more of an immune reserve against pathogenic microorganisms. In the inflammatory state, the spectrum of cytokines secreted by various types of cells in the dental pulp tended to be identical, such that it mainly resisted pathogenic microorganisms. CONCLUSIONS: The cytokine secretion profiles of various cell types in healthy and inflamed dental pulp at the single-cell level are summarized.

Screening of potential drugs for the treatment of diabetic kidney disease using single-cell transcriptome sequencing and connectivity map data.

OBJECTIVE: To explore the feasibility of screening potential drugs for the treatment of diabetic kidney disease (DKD) using a single-cell transcriptome sequencing dataset and Connectivity Map (CMap) database screening. METHODS: A DKD single-nucleus transcriptome sequencing dataset was analyzed using Seurat 4.0 to obtain specific podocyte subclusters and differentially expressed genes (DEGs) related to DKD. These DEGs were subsequently subjected to a search against the CMap database to screen for drug candidates. Cell and animal experiments were conducted to evaluate the efficacy of the top 3 drug candidates. RESULTS: Initially, we analyzed the DKD single-nucleus transcriptome sequencing dataset to obtain intrinsic renal cells such as podocytes, endothelial cells, mesangial cells, proximal tubular cells, collecting duct cells and immune cells. Podocytes were further divided into four subclusters, among which the proportion of POD_1 podcytes was significantly greater in DKD kidneys than in control kidneys (34.0 % vs. 3.4 %). The CMap database was searched using the identified DEGs in the POD_1 subcluster, and the drugs, including tozasertib, paroxetine, and xylazine, were obtained. Cell-based experiments showed that tozasertib, paroxetine and xylazine had no significant podocyte toxicity in the concentration range of 0.01-50 µM. Tozasertib, paroxetine, and xylazine all reversed the advanced glycation end products (AGEs)-induced decrease in podocyte marker levels, but the effect of paroxetine was more prominent. Animal experiments showed that paroxetine decreased urine ALB/Cr levels in DKD model mice by approximately 51.5 % (115.7 mg/g vs. 238.8 mg/g, P < 0.05). Histopathological assessment revealed that paroxetine attenuated basement membrane thickening, restored the number of foot processes of podocytes, and reduced foot process fusion. In addition, paroxetine also attenuated renal tubular-interstitial fibrosis. Mechanistically, paroxetine inhibited the expression of GRK2 and NLRP3, decreased the phosphorylation level of p65, restored NRF2 expression, and relieved inflammation and oxidative stress. CONCLUSION: This strategy based on single-cell transcriptome sequencing and CMap data can facilitate the identification and aid the rapid development of clinical DKD drugs. Paroxetine, screened by this strategy, has excellent renoprotective effects.

Study on the causes of changes in colour during Hibiscus syriacus flowering based on transcriptome and metabolome analyses.

BACKGROUND: The flower colour of H. syriacus 'Qiansiban' transitions from fuchsia to pink-purple and finally to pale purple, thereby enhancing the ornamental value of the cultivars. However, the molecular mechanism underlying this change in flower colour in H. syriacus has not been elucidated. In this study, the transcriptomic data of H. syriacus 'Qiansiban' at five developmental stages were analysed to investigate the impact of flavonoid components on flower colour variation. Additionally, five cDNA libraries were constructed from H. syriacus 'Qiansiban' during critical blooming stages, and the transcriptomes were sequenced to investigate the molecular mechanisms underlying changes in flower colouration. RESULTS: High-performance liquid chromatography‒mass spectrometry detected five anthocyanins in H. syriacus 'Qiansiban', with malvaccin-3-O-glucoside being the predominant compound in the flowers of H. syriacus at different stages, followed by petunigenin-3-O-glucoside. The levels of these five anthocyanins exhibited gradual declines throughout the flowering process. In terms of the composition and profile of flavonoids and flavonols, a total of seven flavonoids were identified: quercetin-3-glucoside, luteolin-7-O-glucoside, Santianol-7-O-glucoside, kaempferol-O-hexosyl-C-hexarbonoside, apigenin-C-diglucoside, luteolin-3,7-diglucoside, and apigenin-7-O-rutinoside. A total of 2,702 DEGs were identified based on the selected reference genome. Based on the enrichment analysis of differentially expressed genes, we identified 9 structural genes (PAL, CHS, FLS, DRF, ANS, CHI, F3H, F3'5'H, and UFGT) and 7 transcription factors (3 MYB, 4 bHLH) associated with flavonoid biosynthesis. The qRT‒PCR results were in good agreement with the high-throughput sequencing data. CONCLUSION: This study will establish a fundamental basis for elucidating the mechanisms underlying alterations in the flower pigmentation of H. syriacus.

Transcriptome analysis of lipid biosynthesis during kernel development in two walnut (Juglans regia L.) varieties of 'Xilin 3' and 'Xiangling'.

BACKGROUND: Walnut is an oilseed tree species and an ecologically important woody tree species that is rich in oil and nutrients. In light of differences in the lipid content, fatty acid composition and key genes expression patterns in different walnut varieties, the key gene regulatory networks for lipid biosynthesis in different varieties of walnuts were intensively investigated. RESULTS: The kernels of two walnut varieties, 'Xilin 3' (X3) and 'Xiangling' (XL) were sampled at 60, 90, and 120 days post-anthesis (DPA) to construct 18 cDNA libraries, and the candidate genes related to oil synthesis were identified via sequencing and expression analysis. A total of 106 differentially expressed genes associated with fatty acid biosynthesis, fatty acid elongation, unsaturated fatty acid biosynthesis, triglyceride assembly, and oil body storage were selected from the transcriptomes. Weighted gene co-expression network analysis (WGCNA), correlation analysis and quantitative validation confirmed the key role of the FAD3 (109002248) gene in lipid synthesis in different varieties. CONCLUSIONS: These results provide valuable resources for future investigations and new insights into genes related to oil accumulation and lipid metabolism in walnut seed kernels. The findings will also aid future molecular studies and ongoing efforts to genetically improve walnut.

Transcriptomic analysis reveals the regulatory mechanisms of messenger RNA (mRNA) and long non-coding RNA (lncRNA) in response to waterlogging stress in rye (Secale cereale L.).

BACKGROUND: Waterlogging stress (WS) negatively impacts crop growth and productivity, making it important to understand crop resistance processes and discover useful WS resistance genes. In this study, rye cultivars and wild rye species were subjected to 12-day WS treatment, and the cultivar Secale cereale L. Imperil showed higher tolerance. Whole transcriptome sequencing was performed on this cultivar to identify differentially expressed (DE) messenger RNAs (DE-mRNAs) and long non-coding RNAs (DE-lncRNAs) involved in WS response. RESULTS: Among the 6 species, Secale cereale L. Imperil showed higher tolerance than wild rye species against WS. The cultivar effectively mitigated oxidative stress, and regulated hydrogen peroxide and superoxide anion. A total of 728 DE-mRNAs and 60 DE-lncRNAs were discovered. Among these, 318 DE-mRNAs and 32 DE-lncRNAs were upregulated, and 410 DE-mRNAs and 28 DE-lncRNAs were downregulated. GO enrichment analysis discovered metabolic processes, cellular processes, and single-organism processes as enriched biological processes (BP). For cellular components (CC), the enriched terms were membrane, membrane part, cell, and cell part. Enriched molecular functions (MF) terms were catalytic activity, binding, and transporter activity. LncRNA and mRNA regulatory processes were mainly related to MAPK signaling pathway-plant, plant hormone signal transduction, phenylpropanoid biosynthesis, anthocyanin biosynthesis, glutathione metabolism, ubiquitin-mediated proteolysis, ABC transporter, Cytochrome b6/f complex, secondary metabolite biosynthesis, and carotenoid biosynthesis pathways. The signalling of ethylene-related pathways was not mainly dependent on AP2/ERF and WRKY transcription factors (TF), but on other factors. Photosynthetic activity was active, and carotenoid levels increased in rye under WS. Sphingolipids, the cytochrome b6/f complex, and glutamate are involved in rye WS response. Sucrose transportation was not significantly inhibited, and sucrose breakdown occurs in rye under WS. CONCLUSIONS: This study investigated the expression levels and regulatory functions of mRNAs and lncRNAs in 12-day waterlogged rye seedlings. The findings shed light on the genes that play a significant role in rye ability to withstand WS. The findings from this study will serve as a foundation for further investigations into the mRNA and lncRNA WS responses in rye.