Atypical BCR-ABL1 transcript in mixed phenotype acute leukemia with bone marrow necrosis.

Mixed phenotype acute leukemia (MPAL) is a type of acute leukemia in which encompasses mixed features of myeloid, T-lymphoid, and/or B-lymphoid differentiation. Philadelphia chromosome-positive (Ph+) MPAL is a rare subgroup with a poor prognosis and accounts for ＜1% of adult acute leukemia. Until now, there is still no consensus on how to best treat Ph+ MPAL. Here, we report a 62-year-old male with Ph+ (atypical e13a2 BCR-ABL1 fusion protein) MPAL. This patient presented with recurrent and intense bone pain due to bone marrow necrosis (BMN). Besides, he did not achieve a complete remission for the first two chemotherapies, until he received flumatinib combined with hyper-CVAD (B) (a dose-intensive regimen include methotrexate and cytarabine). To our knowledge, this is the first report to describe the coexistence of BMN and atypical e13a2 BCR-ABL1 transcripts in patients with MPAL. This finding will bring new understandings in the diagnosis and treatment of Ph+ MPAL.

Prognostic Factors for Recurrent Glioma: A Population-Based Analysis.

Background: The overall survival (OS) for patients with recurrent glioma is meager. Also, the effect of radionecrosis and prognostic factors for recurrent glioma remains controversial. In this regard, developing effective predictive models and guiding clinical care is crucial for these patients. Methods: We screened patients with recurrent glioma after radiotherapy and those who received surgery between August 1, 2013, and December 31, 2020. Univariate and multivariate Cox regression analyses determined the independent prognostic factors affecting the prognosis of recurrent glioma. Moreover, nomograms were constructed to predict recurrent glioma risk and prognosis. Statistical methods were used to determine the prediction accuracy and discriminability of the nomogram prediction model based on the area under the curve (AUC), the C-index, the decision curve analysis (DCA), and the calibration curve. In order to distinguish high-risk and low-risk groups for OS, the X-Tile and Kaplan-Meier (K-M) survival curves were employed, and the nomogram prediction model was further validated by the X-Tile and K-M survival curves. Results: According to a Cox regression analysis, independent prognostic factors of recurrent glioma after radiotherapy with radionecrosis were World Health Organization (WHO) grade and gliosis percentage. We utilized a nomogram prediction model to analyze results visually. The C-index was 0.682 (95% CI: 0.616-0.748). According to receiver operating characteristic (ROC) analysis, calibration plots, and DCA, the nomogram prediction model was found to have a high-performance ability, and all patients were divided into low-risk and high-risk groups based on OS (P < .001). Conclusion: WHO grade and gliosis percentage are prognostic factors for recurrent glioma with radionecrosis, and a nomogram prediction model was established based on these two variables. Patients could be divided into high- and low-risk groups with different OS by this model, and it will provide individualized clinical decisions for future treatment.

From Detection to Cure - Emerging Roles for Urinary Tumor DNA (utDNA) in Bladder Cancer.

PURPOSE OF REVIEW: This review sought to define the emerging roles of urinary tumor DNA (utDNA) for diagnosis, monitoring, and treatment of bladder cancer. Building from early landmark studies the focus is on recent studies, highlighting how utDNA could aid personalized care. RECENT FINDINGS: Recent research underscores the potential for utDNA to be the premiere biomarker in bladder cancer due to the constant interface between urine and tumor. Many studies find utDNA to be more informative than other biomarkers in bladder cancer, especially in early stages of disease. Points of emphasis include superior sensitivity over traditional urine cytology, broad genomic and epigenetic insights, and the potential for non-invasive, real-time analysis of tumor biology. utDNA shows promise for improving all phases of bladder cancer care, paving the way for personalized treatment strategies. Building from current research, future comprehensive clinical trials will validate utDNA's clinical utility, potentially revolutionizing bladder cancer management.

Diagnostic accuracy and clinical impact of renal biopsy cytology.

BACKGROUND: The role of fine needle biopsy cytology in the workup of renal mass lesions remains controversial. With advances in imaging technology and clinical management for renal masses, a critical reevaluation of the role of renal biopsy is needed. This study was designed to provide a comprehensive evaluation of the performance and clinical impact of fine needle biopsy in patients with renal masses. METHODS: A 5-year retrospective study of ultrasound or computer tomography (CT)-guided fine needle biopsies of renal masses diagnosed via cytopathology was conducted. Overall diagnostic rate, sensitivity, and diagnostic accuracy were calculated. Further analysis of the impact of fine needle biopsy cytology on clinical management was performed. RESULTS: A total of 227 cases of fine-needle aspiration and/or biopsy (FNA/B) of renal masses were identified, including 76 with subsequent nephrectomies. Complications were rare (<1%). The diagnostic rate and sensitivity of FNA/B were 83.3% and 89.5%, respectively. Diagnostic accuracy was 98.7% at the major categorical level and 94.7% at the tumor subtype level. Subsequent clinical actions were associated with a definitive cytologic diagnosis of malignancy/neoplasia (p < .05) and were affected by tumor subtype (p < .05). CONCLUSION: This study demonstrates that FNA/B of renal masses is a safe and reliable minimally invasive diagnostic tool with excellent accuracy in confirmation of malignancy and subclassification of tumors. Diagnoses made on FNA/B play a key role in guiding a personalized clinical treatment plan.

A rapid and sensitive one-pot platform integrating fluorogenic RNA aptamers and CRISPR-Cas13a for visual detection of monkeypox virus.

The recent global upsurge in Monkeypox virus (MPXV) outbreaks underscores the critical need for rapid and precise diagnostic solutions, particularly in resource-constrained settings. The gold standard diagnostic method, qRT-PCR, is hindered by its time-consuming nature, requirement for nucleic acid purification, expensive equipment, and the need for highly trained personnel. Traditional CRISPR/Cas fluorescence assays, relying on trans-cleavage of ssDNA/RNA reporters labeled with costly fluorophores and quenchers, pose challenges that limit their widespread application, especially for point-of-care testing (POCT). In this study, we utilized a cost-effective and stable fluorogenic RNA aptamer (Mango III), specifically binding and illuminating the fluorophore TO3-3 PEG-Biotin Fluorophore (TO3), as a reporter for Cas13a trans-cleavage activity. We propose a comprehensive strategy integrating RNA aptamer, recombinase-aided amplification (RAA), and CRISPR-Cas13a systems for the molecular detection of MPXV target. Leveraging the inherent collateral cleavage properties of the Cas13a system, we established high-sensitivity and specificity assays to distinguish MPXV from other Orthopoxviruses (OPVs). A streamlined one-pot protocol was developed to mitigate aerosol contamination risks. Our aptamer-coupled RAA-Cas13a one-pot detection method achieved a Limit of Detection (LoD) of 4 copies of target MPXV DNA in just 40 min. Validation using clinical MPX specimens confirmed the rapid and reliable application of our RAA-Cas13a-Apt assays without nucleic acid purification procedure, highlighting its potential as a point-of-care testing solution. These results underscore the user-friendliness and effectiveness of our one-pot RAA-Cas13a-Apt diagnostic platform, poised to revolutionize disease detection and management.

LOX-1 acts as an N6-methyladenosine-regulated receptor for Helicobacter pylori by binding to the bacterial catalase.

The role of N6-methyladenosine (m6A) modification of host mRNA during bacterial infection is unclear. Here, we show that Helicobacter pylori infection upregulates host m6A methylases and increases m6A levels in gastric epithelial cells. Reducing m6A methylase activity via hemizygotic deletion of methylase-encoding gene Mettl3 in mice, or via small interfering RNAs targeting m6A methylases, enhances H. pylori colonization. We identify LOX-1 mRNA as a key m6A-regulated target during H. pylori infection. m6A modification destabilizes LOX-1 mRNA and reduces LOX-1 protein levels. LOX-1 acts as a membrane receptor for H. pylori catalase and contributes to bacterial adhesion. Pharmacological inhibition of LOX-1, or genetic ablation of Lox-1, reduces H. pylori colonization. Moreover, deletion of the bacterial catalase gene decreases adhesion of H. pylori to human gastric sections. Our results indicate that m6A modification of host LOX-1 mRNA contributes to protection against H. pylori infection by downregulating LOX-1 and thus reducing H. pylori adhesion.

Retraction Notice to: Upregulation of lncRNA LINC00460 Facilitates GC Progression through Epigenetically Silencing CCNG2 by EZH2/LSD1 and Indicates Poor Outcomes.

[This retracts the article DOI: 10.1016/j.omtn.2019.12.041.].

Characteristics of Molecular Genetic Mutations and Their Correlation with Prognosis in Adolescent and Adult Patients with Acute Lymphoblastic Leukemia.

INTRODUCTION: The prognosis of acute lymphoblastic leukemia (ALL) in adolescents and adults is poor, and recurrence is an important cause of their death. Changes of genetic information play a vital role in the pathogenesis and recurrence of ALL; however, the impact of molecular genetic mutations on disease diagnosis and prognosis remains unexplored. This study aimed to explore the frequency spectrum of gene mutations and their prognostic significance, along with the minimal residual disease (MRD) level and hematopoietic stem cell transplantation (HSCT), in adolescent and adult patients aged ≥15 years with ALL. METHODS: The basic characteristics, cytogenetics, molecular genetics, MRD level, treatment regimen, and survival outcome of patients with untreated ALL (≥15 years) were collected, and the correlation and survival analysis were performed using the SPSS 25.0 and R software. RESULTS: This study included 404 patients, of which 147 were selected for next-generation sequencing (NGS). NGS results revealed that 91.2% of the patients had at least one mutation, and 67.35% had multiple (≥2) mutations. NOTCH1, PHF6, RUNX1, PTEN, JAK3, TET2, and JAK1 were the most common mutations in T-ALL, whereas FAT1, TET2, NARS, KMT2D, FLT3, and RELN were the most common mutations in B-ALL. Correlation analysis revealed the mutation patterns, which were significantly different between T-ALL and B-ALL. In the prognostic analysis of 107 patients with B-ALL, multivariate analysis showed that the number of mutations ≥5 was an independent risk factor for overall survival and the RELN mutation was an independent poor prognostic factor for event-free survival. DISCUSSION: The distribution of gene mutations and the co-occurrence and repulsion of mutant genes in patients with ALL were closely related to the immunophenotype of the patients. The number of mutations ≥5 and the RELN mutation were significantly associated with poor prognosis in adolescent and adult patients with ALL.

An experimental study of the correlation between P-wave velocity and the physical properties of weakly cemented formations.

Deep water and shallow layers mostly feature weakly cemented formations, with complex geological structures, geological looseness, susceptibility to collapse. In order to obtain information on weakly cemented formation materials, weakly cemented argillaceous siltstone is simulated as the research object and the focus is on analysing the influence of ultrasonic frequency, density, particle size (porosity), and compressive strength on P-wave velocity and establishing the correlation relationship between longitudinal wave velocity and each parameter through indoor simulation experiments. The results showed that there is a linear relationship between P-wave velocity and ultrasonic frequency in terms of positive correlation as well as compressive strength. The nonlinear relationship between P-wave velocity and particle size (porosity) is a negative correlation, while the nonlinear relationship between sound velocity and density is a positive correlation. In addition, the influence of core height on P-wave velocity is analysed; it is found that as the core height increases, the velocity slightly decreases, and each ultrasonic frequency has an ultimate height for sound wave penetration. Through the response relationship between ultrasound and the physical properties of weakly cemented formations, indirect acquisition can be achieved, which is of great significance for the development of oil and gas in weakly cemented formations.

Elevation of cancer antigen 15-3 owing to oncocytic renal neoplasm in a patient without evidence of breast cancer recurrence on follow-up: a case report.

BACKGROUND: Cancer antigen 15-3 is a protein that clinicians commonly measure to monitor outcomes and response to treatment in patients with breast cancer. However, cancer antigen 15-3 can also be elevated in other, benign and malignant conditions. CASE PRESENTATION: A 73-year-old White woman with history of breast cancer presented to her primary care physician with right hip pain, and laboratory testing revealed elevated cancer antigen 15-3. Further workup with radiographic imaging revealed a large mass in her right kidney. The renal mass was subsequently removed, and the cancer antigen 15-3 level returned to normal. CONCLUSIONS: Elevation of cancer antigen 15-3 owing to causes other than breast cancer recurrence can be a potential diagnostic pitfall during a patient's follow-up. It is important for clinicians to be aware of the limitations of cancer markers and to utilize a combination of diagnostic tests for patient evaluation.

Nonlinear error reduction for phase-shifting profilometry considering periodicity and symmetry of a phase histogram.

Phase-shifting profilometry is extensively utilized for three-dimensional (3D) measurement. However, because of gamma nonlinearity, the image intensities of the captured fringe patterns are regrettably distorted. An effective nonlinear error reduction method without requiring parameter estimation is presented in this paper. Differing from the traditional whole-period phase histogram equalization (PHE) method, our method takes into account not only the periodicity but also the symmetry of the phase histogram. Taking a three-step phase-shifting algorithm as an example, the phase error frequency triples the fringe frequency; thus, we first propose a 1/3-period PHE method. Moreover, since the phase error distribution is sinusoidal with symmetry, we further propose a 1/6-period PHE method. Simulations and experiments both indicate that the 1/6-period PHE method, compared with the whole-period PHE and 1/3-period PHE methods, can further reduce the nonlinear error.

Comparison of the efficacy and safety of an oral sulfate solution and 3-L polyethylene glycol on bowel preparation before colonoscopy: a phase III multicenter randomized controlled trial.

BACKGROUND AND AIMS: Adequate bowel preparation is crucial for clear mucosal visualization during colonoscopy. We aimed to comprehensively compare oral sulfate solution (OSS) and 3-L split-dose polyethylene glycol (PEG) for bowel preparation before colonoscopy. METHODS: This randomized, active-controlled, noninferiority study was performed in 10 medical centers. Eligible subjects were enrolled to receive OSS or 3-L PEG in a split-dose regimen. The quality of bowel preparation, adverse reactions, and acceptability were evaluated. The quality of bowel preparation was evaluated using the Boston Bowel Preparation Scale. Safety was evaluated by adverse reactions. The study population was divided into the full analysis set (FAS), the safety set, the modified FAS (mFAS), and the per-protocol set (PPS). RESULTS: Three hundred forty-eight potentially eligible subjects were enrolled. Three hundred forty-four subjects were included in the FAS and safety set, 340 subjects were included in the mFAS, and 328 subjects were included in the PPS. Adequate bowel preparation of the OSS was not inferior to 3-L PEG in the mFAS (98.22% vs 97.66%) and the PPS (98.17% vs 98.78%). There was no significant difference in acceptability between the 2 groups (94.74% vs 94.80%, P = .9798). Overall adverse reactions were similar (50.88% vs 44.51%, P = .2370) between the 2 groups. CONCLUSIONS: The split-dose OSS regimen was not inferior to the split-dose 3-L PEG regimen for the quality of bowel preparation in a Chinese adult population. The safety and acceptability of the 2 groups were similar. (Clinical trial registration number: NCT05465889.).

Targeting N6-methyladenosine reader YTHDF1 with siRNA boosts antitumor immunity in NASH-HCC by inhibiting EZH2-IL-6 axis.

BACKGROUND & AIMS: RNA N6-methyladenosine (m6A) reader protein YTHDF1 has been implicated in cancer; however, its role in hepatocellular carcinoma (HCC), especially in non-alcoholic steatohepatitis-associated HCC (NASH-HCC), remains unknown. Here, we investigated the functional role of YTHDF1 in NASH-HCC and its interplay with the tumor immune microenvironment. METHODS: Hepatocyte-specific Ythdf1-overexpressing mice were subjected to a NASH-HCC-inducing diet. Tumor-infiltrating immune cells were profiled with single-cell RNA-sequencing, flow cytometry, and immunostaining. The molecular target of YTHDF1 was elucidated with RNA-sequencing, m6A-sequencing, YTHDF1 RNA immunoprecipitation-sequencing, proteomics, and ribosome-profiling. Ythdf1 in NASH-HCC models was targeted by lipid nanoparticle (LNP)-encapsulated small-interfering Ythdf1. RESULTS: YTHDF1 is overexpressed in tumor tissues compared to adjacent peri-tumor tissues from patients with NASH-HCC. Liver-specific Ythdf1 overexpression drives tumorigenesis in dietary models of spontaneous NASH-HCC. Single-cell RNA-sequencing and flow cytometry revealed that Ythdf1 induced accumulation of myeloid-derived suppressor cells (MDSCs) and suppressed cytotoxic CD8+ T-cell function. Mechanistically, Ythdf1 expression in NASH-HCC cells induced the secretion of IL-6, which mediated MDSC recruitment and activation, leading to CD8+ T-cell dysfunction. EZH2 mRNA was identified as a key YTHDF1 target. YTHDF1 binds to m6A-modified EZH2 mRNA and promotes EZH2 translation. EZH2 in turn increased expression and secretion of IL-6. Ythdf1 knockout synergized with anti-PD-1 treatment to suppress tumor growth in NASH-HCC allografts. Furthermore, therapeutic targeting of Ythdf1 using LNP-encapsulated small-interfering RNA significantly increased the efficacy of anti-PD-1 blockade in NASH-HCC allografts. CONCLUSIONS: We identified that YTHDF1 promotes NASH-HCC tumorigenesis via EZH2-IL-6 signaling, which recruits and activates MDSCs to cause cytotoxic CD8+ T-cell dysfunction. YTHDF1 may be a novel therapeutic target to improve responses to anti-PD-1 immunotherapy in NASH-HCC. IMPACT AND IMPLICATIONS: YTHDF1, a N6-methyladenosine reader, is upregulated in patients with non-alcoholic steatohepatitis (NASH)-associated hepatocellular carcinoma (HCC); however, its role in modulating the tumor immune microenvironment in NASH-HCC remains unclear. Here, we show that Ythdf1 mediates immunosuppression in NASH-HCC and that targeting YTHDF1 in combination with immune checkpoint blockade elicits robust antitumor immune responses. Our findings suggest novel therapeutic targets for potentiating the efficacy of immune checkpoint blockade in NASH-HCC and provide the rationale for developing YTHDF1 inhibitors for the treatment of NASH-HCC.

CD4+ T-Cell Subsets in Aplastic Anemia, Myelodysplastic Syndrome, and Acute Myelogenous Leukemia Patients: a Comparative Analysis.

BACKGROUND: Systematic and comparative studies on CD4+ T-lymphocytes in aplastic anemia (AA), myelodysplastic syndrome (MDS), and acute myelogenous leukemia (AML) are scarce. This study aimed to investigate the importance of CD4+ T-cells in bone marrow (BM) failure. METHODS: The proportions of Th1, Th2, Th17, and Treg cells in peripheral blood mononuclear cells (PBMCs) were examined by flow cytometry (FCM). The mRNA expression levels of transcription factors were measured using real-time PCR. RESULTS: The proportions of Th1, Th17 cells, and Th1/Th2 in the AA group were higher, whereas Th2 and Tregs were lower compared to controls. The proportions of Th17 and Treg cells accompanied by RORγt, and Foxp3 expression were significantly higher in the MDS group. The proportions of Th1, Th17, and Th1/Th2 were higher, whereas Th2 cells and GATA3 expression were significantly lower in MDS-multilineage dysplasia group, than in control group. The proportions of Th1, Th17, and Th1/Th2 were lower in MDS-excess blasts, and AML groups, than in controls, whereas that of Th2 and Treg cells accompanied by GATA3, and Foxp3 expression were significantly higher. CONCLUSIONS: Imbalance in CD4+ T-cell subsets may play a critical role in the pathogenesis and BM failure in the investigated diseases.

The epidemiological patterns of non-Hodgkin lymphoma: global estimates of disease burden, risk factors, and temporal trends.

Background: The incidence of non-Hodgkin's lymphoma (NHL) has increased steadily over the past few decades. Elucidating its global burden will facilitate more effective disease management and improve patient outcomes. We explored the disease burden, risk factors, and trends in incidence and mortality in NHL globally. Methods: The up-to-date data on age-standardized incidence and mortality rates of NHL were retrieved from the GLOBOCAN 2020, CI5 volumes I-XI, WHO mortality database, and Global Burden of Disease (GBD) 2019, focusing on geographic disparities worldwide. We reported incidence and mortality by sex and age, along with corresponding age-standardized rates (ASRs), the average annual percentage change (AAPC), and future burden estimates to 2040. Results: In 2020, there were an estimated 545,000 new cases and 260,000 deaths of NHL globally. In addition, NHL resulted in 8,650,352 age-standardized DALYs in 2019 worldwide. The age-specific incidence rates varied drastically across world areas, at least 10-fold in both sexes, with the most pronounced increase trend found in Australia and New Zealand. By contrast, North African countries faced a more significant mortality burden (ASR, 3.7 per 100,000) than highly developed countries. In the past decades, the pace of increase in incidence and mortality accelerated, with the highest AAPC of 4.9 (95%CI: 3.6-6.2) and 6.8 (95%CI: 4.3-9.2) in the elderly population, respectively. Considering risk factors, obesity was positively correlated with age-standardized incidence rates (P< 0.001). And North America was the high-risk region for DALYs due to the high body mass index in 2019. Regarding demographic change, NHL incident cases are projected to rise to approximately 778,000 by 2040. Conclusion: In this pooled analysis, we provided evidence for the growing incidence trends in NHL, particularly among women, older adults, obese populations, and HIV-infected people. And the marked increase in the older population is still a public health issue that requires more attention. Future efforts should be directed at cultivating health awareness and formulating effective and locally tailored cancer prevention strategies, especially in most developing countries.

Natural Products for the Immunotherapy of Glioma.

Glioma immunotherapy has attracted increasing attention since the immune system plays a vital role in suppressing tumor growth. Immunotherapy strategies are already being tested in clinical trials, such as immune checkpoint inhibitors (ICIs), vaccines, chimeric antigen receptor T-cell (CAR-T cell) therapy, and virus therapy. However, the clinical application of these immunotherapies is limited due to their tremendous side effects and slight efficacy caused by glioma heterogeneity, antigen escape, and the presence of glioma immunosuppressive microenvironment (GIME). Natural products have emerged as a promising and safe strategy for glioma therapy since most of them possess excellent antitumor effects and immunoregulatory properties by reversing GIME. This review summarizes the status of current immunotherapy strategies for glioma, including their obstacles. Then we discuss the recent advancement of natural products for glioma immunotherapy. Additionally, perspectives on the challenges and opportunities of natural compounds for modulating the glioma microenvironment are also illustrated.

Identification of ORM1, vWF, SPARC, and PPBP as immune-related proteins involved in immune thrombocytopenia by quantitative LC-MS/MS.

BACKGROUND: Immune thrombocytopenia (ITP) is a common autoimmune disease characterized by loss of immune tolerance to platelet autoantigens leading to excessive destruction and insufficient production of platelets. METHOD: Quantitative liquid chromatography tandem mass spectrometry (LC-MS/MS) was performed to detect the differentially expressed proteins in bone marrow samples from active ITP patients and normal controls. RESULT: Our bioinformatic analysis identified two upregulated proteins (ORM1 and vWF) and two downregulated proteins (PPBP and SPARC) related to immune function. The four proteins were all found to be related to the tumor necrosis factor (TNF) -α signalling pathway and involved in the pathogenesis of ITP in KEGG pathway analysis. CONCLUSION: Bioinformatics analysis identified differentially expressed proteins in bone marrow that are involved in the TNF-α signalling pathway and are related to the activation of immune function in ITP patients. These findings could provide new ideas for research on the loss of immune tolerance in ITP patients.

The mushroom body development and learning ability of adult honeybees are influenced by cold exposure during their early pupal stage.

The honeybees are the most important pollinator in the production of crops and fresh produce. Temperature affects the survival of honeybees, and determines the quality of their development, which is of great significance for beekeeping production. Yet, little was known about how does low temperature stress during development stage cause bee death and any sub-lethal effect on subsequent. Early pupal stage is the most sensitive stage to low temperature in pupal stage. In this study, early pupal broods were exposed to 20°C for 12, 16, 24, and 48 h, followed by incubation at 35°C until emergence. We found that 48 h of low temperature duration cause 70% of individual bees to die. Although the mortality at 12 and 16 h seems not very high, the association learning ability of the surviving individuals was greatly affected. The brain slices of honeybees showed that low temperature treatment could cause the brain development of honeybees to almost stop. Gene expression profiles between low temperature treatment groups (T24, T48) and the control revealed that 1,267 and 1,174 genes were differentially expressed respectively. Functional enrichment analysis of differentially expressed genes showed that the differential expression of Map3k9, Dhrs4, and Sod-2 genes on MAPK and peroxisome signaling pathway caused oxidative damage to the honeybee head. On the FoxO signal pathway, InsR and FoxO were upregulated, and JNK, Akt, and Bsk were downregulated; and on the insect hormone synthesis signal pathway, Phm and Spo genes were downregulated. Therefore, we speculate that low temperature stress affects hormone regulation. It was detected that the pathways related to the nervous system were Cholinergic synapse, Dopaminergic synapse, GABAergic synapse, Glutamatergic synapse, Serotonergic synapse, Neurotrophin signaling pathway, and Synaptic vesicle cycle. This implies that the synaptic development of honeybees is quite possibly greatly affected by low temperature stress. Understanding how low temperature stress affects the physiology of bee brain development and how it affects bee behavior provide a theoretical foundation for a deeper comprehension of the temperature adaptation mechanism that underlies the "stenothermic" development of social insects, and help to improve honeybee management strategies to ensure the healthy of colony.

WGCNA based identification of hub genes associated with cold response and development in Apis mellifera metamorphic pupae.

Honeybee is a crucial pollinator in nature, and plays an indispensable role in both agricultural production and scientific research. In recent decades, honeybee was challenged with health problems by biotic and abiotic stresses. As a key ecological factor, temperature has been proved to have an impact on the survival and production efficiency of honeybees. Previous studies have demonstrated that low temperature stress can affect honeybee pupation and shorten adult longevity. However, the molecular mechanism underlying the effects of low temperatures on honeybee growth and development during their developmental period remain poorly understood. In this paper, the weighted gene co-expression analysis (WGCNA) was employed to explore the molecular mechanisms underpinnings of honeybees' respond to low temperatures (20°C) during four distinct developmental stages: large-larvae, prepupae, early-pupae and mid-pupae. Through an extensive transcriptome analysis, thirteen gene co-expression modules were identified and analyzed in relation to honeybee development and stress responses. The darkorange module was found to be associated with low temperature stress, with its genes primarily involved in autophagy-animal, endocytosis and MAPK signaling pathways. Four hub genes were identified within this module, namely, loc726497, loc409791, loc410923, and loc550857, which may contribute to honeybee resistance to low temperature and provide insight into the underlying mechanism. The gene expression patterns of grey60 and black modules were found to correspond to the developmental stages of prepupae and early-pupae, respectively, with the hub genes loc409494, loc725756, loc552457, loc726158, Ip3k and Lcch3 in grey60 module likely involved in brain development, and the hub genes loc410555 in black module potentially related to exoskeleton development. The brown module genes exhibited a distinct pattern of overexpression in mid-pupae specimens, with genes primarily enriched in oxidative phosphorylation, citrate cycle and other pathways, which may be related to the formation of bee flying muscle. No related gene expression module was found for mature larvae stage. These findings provide valuable insights into the developmental process of honeybees at molecular level during the capped brood stage.