ObjectiveTo establish a reliable chronic obstructive pulmonary disease (COPD) mouse model based on a self-developed multichannel automatic control system for long-term continuous cigarette smoke exposure in small animals using a novel continuous cigarette smoke exposure method, and to conduct phenotypic evaluation and analysis, thereby providing an animal experimental basis for investigating COPD pathogenesis and prevention strategies. MethodsTwenty male C57BL/6J mice aged 6 weeks were randomly and equally divided into a control group and a model group. The model group (n=10) underwent 6 h of continuous cigarette smoke exposure daily (6 cigarettes per day for 12 consecutive weeks), while the control group (n=10) received no intervention. Body weight was monitored biweekly. Post-exposure, in vivo micro-CT imaging was performed. After euthanasia, serum and bronchoalveolar lavage fluid (BALF) levels of interleukin-6 (IL-6) and tumor necrosis factor-α (TNF-α) were quantified by ELISA. Lung tissues underwent H&E and Masson's trichrome staining to observe changes in lung morphology and inflammatory cell infiltration, and the mean linear intercept (MLI) was calculated, thereby comprehensively evaluating the clinical features of COPD in the mouse model. ResultsCompared with the control group, the model group showed significantly reduced body weight (P<0.01) from the fourth week. Compared with the control group, IL-6 level in the serum and BALF of the model group increased by 27.2% and 140.0%, respectively (P<0.01). TNF-α level in the serum and bronchoalveolar lavage fluid of the model group increased by 16.7% (P<0.01) and 19.3% (P<0.05), respectively. Histopathological examination revealed alveolar wall thinning, septal rupture, emphysematous bullae formation, reduced alveolar count, bronchial wall thickening with lumen narrowing, and inflammatory cell infiltration. MLI was significantly elevated (P<0.01). Masson's staining confirmed collagen deposition and bronchial remodeling. Micro-CT demonstrated localized high-density shadows exhibiting typical features of chronic bronchitis. Conclusion The self-developed device enables long-term continuous smoke exposure, and the successfully established COPD mouse model exhibits pathological features highly consistent with clinical manifestations, offering an efficient and reliable tool for COPD research.

ObjectiveGastrodia elata has evolved ecological types with shortened rhizome internodes and diversified flower and fruit coloration in response to different altitudes. Studying the genetic mechanisms of different ecotype germplasm is significant for guiding variety breeding in different cultivation areas. MethodsThe bHLH gene family was identified based on the whole-genome datasets of G. elata f. elata and G. elata f. glauca. Subsequently， the gene family members were subject to analysis， including gene structure， chromosomal localization， cis-acting elements， gene synteny， and phylogeny. Combined with transcriptome data and quantitative Real-time PCR， the expression patterns of bHLH genes in the stems of the different G. elata ecotype germplasm were analyzed. Finally， correlation analysis was conducted between gene expression patterns and color to obtain the key bHLH genes regulating the color formation of stem. ResultsA total of 63 bHLH genes were identified in both G elata f. elata and G. elata f. glauca， unevenly distributed across 17 chromosomes and clustered into 16 subfamilies， with significant expansion in some family members. Obvious inversions of bHLH genes on the same chromosome and interchromosomal translocations were detected in the two ecotype germplasm. Among these genes， 12 bHLH genes （such as bHLH62-3 and bHLH74） were associated with the bright yellow color of G elata f. elata stem， while 9 bHLH genes （such as PIL13， UNE12， and bHLH130） were correlated with the red color of G. elata f. glauca stem. Compared to G. elata f. glauca， the bHLH48 expression level was significantly higher in flowers and scale leaves of G elata f. elata， and the bHLH62-3 expression level was significantly higher in all organs of G elata f. elata. ConclusionsFunctional pathway divergence of the bHLH family members has occurred across different chromosomes in G elata f. elata and G. elata f. glauca. Through synergism or antagonism with other genes， 21 bHLH genes participate in the coloration metabolic pathway regulation of stems， flowers， and fruits. Specifically， bHLH62-3 is involved in regulating stem color differentiation in the anthocyanin biosynthesis pathway of G. elata， thus relevant to the color formation of stem. Additionally， GebHLH48 positively regulates flowering-related pathways to promote the early-flowering phenotype of G. elata f. elata. These findings have laid the foundation for analyzing the genetic regulatory mechanisms underlying the color formation of the G. elata stem.

ObjectiveGastrodia elata has evolved ecological types with shortened rhizome internodes and diversified flower and fruit coloration in response to different altitudes. Studying the genetic mechanisms of different ecotype germplasm is significant for guiding variety breeding in different cultivation areas. MethodsThe bHLH gene family was identified based on the whole-genome datasets of G. elata f. elata and G. elata f. glauca. Subsequently， the gene family members were subject to analysis， including gene structure， chromosomal localization， cis-acting elements， gene synteny， and phylogeny. Combined with transcriptome data and quantitative Real-time PCR， the expression patterns of bHLH genes in the stems of the different G. elata ecotype germplasm were analyzed. Finally， correlation analysis was conducted between gene expression patterns and color to obtain the key bHLH genes regulating the color formation of stem. ResultsA total of 63 bHLH genes were identified in both G elata f. elata and G. elata f. glauca， unevenly distributed across 17 chromosomes and clustered into 16 subfamilies， with significant expansion in some family members. Obvious inversions of bHLH genes on the same chromosome and interchromosomal translocations were detected in the two ecotype germplasm. Among these genes， 12 bHLH genes （such as bHLH62-3 and bHLH74） were associated with the bright yellow color of G elata f. elata stem， while 9 bHLH genes （such as PIL13， UNE12， and bHLH130） were correlated with the red color of G. elata f. glauca stem. Compared to G. elata f. glauca， the bHLH48 expression level was significantly higher in flowers and scale leaves of G elata f. elata， and the bHLH62-3 expression level was significantly higher in all organs of G elata f. elata. ConclusionsFunctional pathway divergence of the bHLH family members has occurred across different chromosomes in G elata f. elata and G. elata f. glauca. Through synergism or antagonism with other genes， 21 bHLH genes participate in the coloration metabolic pathway regulation of stems， flowers， and fruits. Specifically， bHLH62-3 is involved in regulating stem color differentiation in the anthocyanin biosynthesis pathway of G. elata， thus relevant to the color formation of stem. Additionally， GebHLH48 positively regulates flowering-related pathways to promote the early-flowering phenotype of G. elata f. elata. These findings have laid the foundation for analyzing the genetic regulatory mechanisms underlying the color formation of the G. elata stem.

ObjectiveThe biosynthesis of heterophyllin B （HB）， a cyclopeptide from Pseudostellaria heterophylla， is regulated by various abiotic stresses. Elucidating the transcriptional regulatory mechanism underlying HB biosynthesis is of great guiding significance for the directional improvement of P. heterophylla varieties and the enhancement of HB content. MethodsBased on transcriptome data from different tissues of P. heterophylla， transcription factors （TFs） specifically upregulated and highly expressed in the phloem of tuberous roots were screened through a combination of Mfuzz time-series clustering， transcription factor family prediction， and correlation analysis. Quantitative real-time polymerase chain reaction （Real-time PCR） was employed to analyze expression patterns of candidate TFs under abscisic acid （ABA） induction， and the dual-luciferase reporter assay was applied to verify their regulatory effects on HB precursor genes. ResultsContent determination showed that HB accumulated at the highest in the phloem of P. heterophylla tuberous roots （34 μg·g^-1 fresh weight）. Transcriptome analysis identified 15 868 up-regulated differentially expressed genes （DEGs）， among which 4 375 were specifically up-regulated in the phloem. From these， 25 TFs highly expressed in the phloem were screened. Correlation analysis revealed that the expression level of WRKY70 was significantly positively correlated with HB content， the expression of precursor genes prePhHB， and PhPOP1. Additionally， WRKY70 expression was significantly upregulated under ABA induction. Dual-luciferase reporter assay confirmed that WRKY70 specifically activated the transcriptional activity of the prePhHB promoter （Pro-prePhHB）. ConclusionHB is mainly accumulated in the phloem of P. heterophylla tuberous roots. WRKY70 positively regulates HB biosynthesis by directly activating the promoter activity of prePhHB. This study clarifies the key role of WRKY70 in the regulation of HB biosynthesis and provides an important theoretical basis for the in-depth analysis of the molecular regulatory mechanism underlying HB biosynthesis.

Schistosomiasis is prevalent in 78 countries and territories worldwide, while the eastern and western parts of sub-Sahara Africa bear the highest disease burden due to schistosomiasis. Recently, climate change, international trade and travel, urbanization and war have increased the risk of cross-boundary importation and transmission of schistosomiasis, and schistosomiasis has increasingly become a public health concern in non-endemic countries and territories. Biomphalaria straminea, the intermediate host of Schistosoma mansoni, has colonized in southern China and its habitats continue to move northward. In addition, cross-boundary imported cases of schistosomiasis have been reported occasionally in China. However, the real number of cases may be underestimated greatly due to insufficient diagnostic capacity and weak awareness of case reporting for overseas imported schistosomiasis in healthcare facilities. It is necessary to establish a multi-party collaborative mechanism, improve corresponding systems and technical specifications, reinforce surveillance and early warning, and border management, enhance technical reserves and capability building, and improve the awareness of schistosomiasis prevention and healthcare-seeking among entry-exit personnel, in order to effectively address the threat of cross-boundary imported schistosomiasis.

Objective To analyze the structural and phylogenetic characteristics of the mitochondrial genome from Bulinus globosus, so as to provide a theoretical basis for classification and identification of species within the Bulinus genus, and to provide insights into understanding of Bulinus-schistosomes interactions and the mechanisms of parasite transmission. Methods B. globosus samples were collected from the Ruya River basin in Zimbabwe. Mitochondrial DNA was extracted from B. globosus samples and the corresponding libraries were constructed for high-throughput sequencing on the Illumina NovaSeq 6000 platform. After raw sequencing data were subjected to quality control using the fastp software, genome assembly was performed using the A5-miseq and SPAdes tools, and genome annotation was conducted using the MITOS online server. Circular maps and sequence plots of the mitochondrial genome were generated using the CGView and OGDRAW software, and the protein conservation motifs and structures were analyzed using the TBtools software. Base composition and codon usage bias were analyzed and visualized using the software MEGA X and the ggplot2 package in the R software. In addition, a phylogenetic tree was created in the software MEGA X after sequence alignment with the software MAFFT 7, and visualized using the software iTOL. Results The mitochondrial genome of B. globosus was a 13 730 bp double-stranded circular molecule, containing 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and 13 protein-coding genes, with a marked AT preference. The mitochondrial genome composition of B. globosus was similar to that of other species within the Bulinus genus. Phylogenetic analysis revealed that the complete mitochondrial genome sequence of B. globosus was clustered with B. truncatus, B. nasutus, and B. ugandae into the same evolutionary clade, and gene superfamily analysis showed that the metabolism-related proteins of B. globosus were highly conserved, notably the cytochrome c oxidase family, which showed a significant consistency. Conclusions This is the first whole mitochondrial genome sequencing to decode the compositional features of the mitochondrial genome of B. globosus from Zimbabwe and its evolutionary relationship within the Bulinus genus, which provides important insights for further understanding of the phylogeny and mitochondrial genome characteristics of the Bulinus genus.

Background/Aims: Early detection and effective prognosis prediction in patients with hepatocellular carcinoma (HCC) provide an avenue for survival improvement, yet more effective approaches are greatly needed. We sought to develop the detection and prognosis models with ultra-sensitivity and low cost based on fragmentomic features of circulating cell free mtDNA (ccf-mtDNA). Methods: Capture-based mtDNA sequencing was carried out in plasma cell-free DNA samples from 1168 participants, including 571 patients with HCC, 301 patients with chronic hepatitis B or liver cirrhosis (CHB/LC) and 296 healthy controls (HC). Results: The systematic analysis revealed significantly aberrant fragmentomic features of ccf-mtDNA in HCC group when compared with CHB/LC and HC groups. Moreover, we constructed a random forest algorithm-based HCC detection model by utilizing ccf-mtDNA fragmentomic features. Both internal and two external validation cohorts demonstrated the excellent capacity of our model in distinguishing early HCC patients from HC and highrisk population with CHB/LC, with AUC exceeding 0.983 and 0.981, sensitivity over 89.6% and 89.61%, and specificity over 98.20% and 95.00%, respectively, greatly surpassing the performance of alpha-fetoprotein (AFP) and mtDNA copy number. We also developed an HCC prognosis prediction model by LASSO-Cox regression to select 20 fragmentomic features, which exhibited exceptional ability in predicting 1-year, 2-year and 3-year survival (AUC=0.8333, 0.8145 and 0.7958 for validation cohort, respectively). Conclusions We have developed and validated a high-performing and low-cost approach in a large clinical cohort based on aberrant ccf-mtDNA fragmentomic features with promising clinical translational application for the early detection and prognosis prediction of HCC patients.

Background/Aims: Early detection and effective prognosis prediction in patients with hepatocellular carcinoma (HCC) provide an avenue for survival improvement, yet more effective approaches are greatly needed. We sought to develop the detection and prognosis models with ultra-sensitivity and low cost based on fragmentomic features of circulating cell free mtDNA (ccf-mtDNA). Methods: Capture-based mtDNA sequencing was carried out in plasma cell-free DNA samples from 1168 participants, including 571 patients with HCC, 301 patients with chronic hepatitis B or liver cirrhosis (CHB/LC) and 296 healthy controls (HC). Results: The systematic analysis revealed significantly aberrant fragmentomic features of ccf-mtDNA in HCC group when compared with CHB/LC and HC groups. Moreover, we constructed a random forest algorithm-based HCC detection model by utilizing ccf-mtDNA fragmentomic features. Both internal and two external validation cohorts demonstrated the excellent capacity of our model in distinguishing early HCC patients from HC and highrisk population with CHB/LC, with AUC exceeding 0.983 and 0.981, sensitivity over 89.6% and 89.61%, and specificity over 98.20% and 95.00%, respectively, greatly surpassing the performance of alpha-fetoprotein (AFP) and mtDNA copy number. We also developed an HCC prognosis prediction model by LASSO-Cox regression to select 20 fragmentomic features, which exhibited exceptional ability in predicting 1-year, 2-year and 3-year survival (AUC=0.8333, 0.8145 and 0.7958 for validation cohort, respectively). Conclusions We have developed and validated a high-performing and low-cost approach in a large clinical cohort based on aberrant ccf-mtDNA fragmentomic features with promising clinical translational application for the early detection and prognosis prediction of HCC patients.

Background/Aims: Early detection and effective prognosis prediction in patients with hepatocellular carcinoma (HCC) provide an avenue for survival improvement, yet more effective approaches are greatly needed. We sought to develop the detection and prognosis models with ultra-sensitivity and low cost based on fragmentomic features of circulating cell free mtDNA (ccf-mtDNA). Methods: Capture-based mtDNA sequencing was carried out in plasma cell-free DNA samples from 1168 participants, including 571 patients with HCC, 301 patients with chronic hepatitis B or liver cirrhosis (CHB/LC) and 296 healthy controls (HC). Results: The systematic analysis revealed significantly aberrant fragmentomic features of ccf-mtDNA in HCC group when compared with CHB/LC and HC groups. Moreover, we constructed a random forest algorithm-based HCC detection model by utilizing ccf-mtDNA fragmentomic features. Both internal and two external validation cohorts demonstrated the excellent capacity of our model in distinguishing early HCC patients from HC and highrisk population with CHB/LC, with AUC exceeding 0.983 and 0.981, sensitivity over 89.6% and 89.61%, and specificity over 98.20% and 95.00%, respectively, greatly surpassing the performance of alpha-fetoprotein (AFP) and mtDNA copy number. We also developed an HCC prognosis prediction model by LASSO-Cox regression to select 20 fragmentomic features, which exhibited exceptional ability in predicting 1-year, 2-year and 3-year survival (AUC=0.8333, 0.8145 and 0.7958 for validation cohort, respectively). Conclusions We have developed and validated a high-performing and low-cost approach in a large clinical cohort based on aberrant ccf-mtDNA fragmentomic features with promising clinical translational application for the early detection and prognosis prediction of HCC patients.

Objective: To investigate the predictive value of EZH2 expression for malignant transformation in oral leukoplakia (OLK) and to provide a reference for clinical practice. Methods: This study was approved by the institutional ethics committee, and informed consent was obtained from all participants. A total of 114 patients diagnosed with OLK by pathological examination and treated at our hospital between November 2020 and July 2022 were initially enrolled. After excluding those with incomplete data or follow-up, 105 participants were included in the final analysis, comprising 14 in the high EZH2 expression group and 91 in the low EZH2 expression group. Histopathological examination of oral mucosa and immunohistochemical detection of EZH2 protein expression were performed. The follow-up period was 30 months; participants were followed until malignant transformation occurred or until the end of follow-up, at which point they were withdrawn from the study. The exposure factor was the level of EZH2 protein expression, and the outcome was the malignant transformation rate of OLK. Differences in EZH2 expression levels and transformation outcomes were analyzed. Results: There were no statistically significant differences between the high and low EZH2 expression groups in terms of age, sex, history of systemic disease, lifestyle habits, psychological status, diet, and sleep conditions (P > 0.05). Lesions in the high EZH2 expression group were mainly located on the ventral tongue, while in the low EZH2 expression group, they were more commonly found on the dorsal tongue and buccal mucosa. The malignant transformation rate was 28.6% (4/14) in the high expression group and 8.8% (8/91) in the low expression group; these differences were not statistically significant (P=0.053). In univariate Cox regression analysis, the risk of malignant transformation in the high EZH2 expression group was 3.647 times that of the low EZH2 expression group (HR = 3.647, 95% CI: 1.097-12.120, P<0.05). Kaplan-Meier survival analysis showed that over the 30-month follow-up period, the cancer-free survival rate in the high EZH2 expression group was 19.8% lower than in the low expression group, and the difference was statistically significant (P<0.05). In multivariate Cox regression analysis, only moderate and severe epithelial dysplasia were identified as independent risk factors for malignant transformation. The risk of malignant transformation in the moderate and severe dysplasia groups was 10.695 and 13.623 times higher, respectively, than in the mild dysplasia group (HR = 10.695, 95% CI: 2.270-50.396, P<0.05; HR=13.623, 95% CI: 1.918-96.774, P<0.05). EZH2 high expression was not an independent risk factor in the multivariate model (HR= 2.528, 95% CI: 0.752-8.500, P = 0.134). Conclusion High EZH2 protein expression is a risk factor for the malignant transformation of OLK but does not have independent predictive value.