Objective To identify factors affecting and key environmental factors of the Oncomelania hupensis snail density in the Yangtze River Delta region using machine learning methods. Methods Administrative village-level O. hupensis snail survey data in the Yangtze River Delta (including Shanghai Municipality, Jiangsu Province, Zhejiang Province and Anhui Province) from 2011 to 2021 were retrieved from the Information Management System for Parasitic Disease Control of Chinese Center for Disease Control and Prevention. Environmental factor data were captured from the Google Earth Engine platform, including elevation, slope, terrain, normalized difference vegetation index (NDVI), vegetation type, soil type, total petroleum hydrocarbon (TPH), ammonium nitrogen, inorganic nitrogen, dissolved oxygen, pH of water, chemical oxygen demand (COD) and inorganic phosphorus, and climatic factor data in the study region were retrieved from the Copernicus Climate Data Store, including annual precipitation, aridity index and annual mean temperature (AMT). O. hupensis snail survey data in the Yangtze River Delta region from 2011 to 2021 were randomly divided into a training set (70%) and a test set (30%), and five machine learning models were selected for machine learning model construction and comparative analysis of the O. hupensis snail density using the software R 4.3.0, including random forest (RF), eXtreme gradient boosting (XGBoost), support vector machine (SVM), gradient boosting machine (GBM) and neural network (NN). The XGBoost model was employed to construct a predictive model for the O. hupensis snail density, and the impact of each environmental factor on O. hupensis snail distribution was quantified. The SHapley Additive exPlanations (SHAPs) values were calculated to estimate the average contribution of each variable to the model prediction, and the core environmental factors affecting the O. hupensis snail population density were screened. Results Among the five machine learning models, the XGBoost model exhibited the optimal comprehensive performance, with the coefficient of determination (R²) of 0.855, mean squared error (MSE) of 0.188, root mean squared error (RMSE) of 0.434 and mean absolute error (MAE) of 0.155, respectively. Analysis of factors affecting the O. hupensis snail density with the XGBoost model showed that among the 16 environmental factors, the top four high-impact factors ranked by SHAPs values included annual precipitation, elevation, aridity index and NDVI, with cumulative SHAPs contributions of 75%, which was higher than that of other environmental factors. If NDVI was higher than 0.6, the O. hupensis snail density increased with NDVI and peaked if NDVI was 0.8 (1.60 snails/0.1 m²). The O. hupensis snail density increased with elevation if the elevation ranged from 14 to 40 m, and slowly rose if the annual precipitation ranged from 900 to 1 300 mm, and then increased rapidly to the peak (1.52 snails/0.1 m²) if the annual precipitation ranged from 1 300 to 1 500 mm. In addition, the O. hupensis snail density increased rapidly to the maximum (1.60 snails/0.1 m²) if the aridity index ranged from 0.8 to 1.1, and decreased gradually if the aridity index exceeded 1.1. Conclusions The XGBoost model shows excellent performance in prediction of the O. hupensis snail density and identification of key environmental factors in the Yangtze River Delta region. Annual precipitation, elevation, aridity index and NDVI are key environmental factors affecting the distribution and density of O. hupensis snails in the Yangtze River Delta region.

Objective To examine the impact of extreme temperature and drought stress on the survival of Bulinus globosus, so as to provide the theoretical evidence for the genomic research of Bulinus in absence of reference genes. Methods B. globosus snail samples were collected from Kiwani Shehia in Pemba Island, Zanzibar, Tanzania, and offspring snails were obtained through laboratory breeding and reproduction. A total of 120 10-week-old B. globosus snails from the same generation were selected and randomly assigned into four groups, including the high-temperature drought (HD) group, normal temperature drought (D) group, low-temperature drought (LD) group, and the control (C) group, of 30 snails in each group. Snails in HD, D, and LD groups were placed in beakers containing dry soil at the bottom and subsequently housed in climate chambers at 35, 26 ℃, and 10 ℃, respectively, while snails in Group C were maintained in 500 mL petri dishes containing dechlorinated tap water at 26 ℃. Following 3 days of breeding, living snails in each group were collected, and soft tissues were dissected and isolated. Total RNA was extracted from snail soft tissues for library construction, followed by high-throughput sequencing on the Illumina HiSeq 4000 sequencing system. De novo transcriptome assembly was performed using the Trinity software, and the longest transcripts were selected as unigenes. Gene functional annotations of unigenes were conducted using the Diamond software against Gene Ontology (GO) knowledgebase, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, NCBI non-redundant (NR) protein sequences database, Protein Family (Pfam) database, and UniProtKB/Swiss-Prot (Swiss-Prot) knowledgebase. GO and KEGG enrichment analyses of differentially expressed genes (DEGs) were performed using the topGO and clusterProfiler software, respectively. In addition, four relevant genes were selected for validation using a real-time quantitative PCR (qRT-PCR) assay to verify the reliability of transcriptome sequencing results. Results Following 3 days of breeding, there were 7, 20, 28, and 30 survival B. globosus snails in HD, LD, D, and C groups, with corresponding survival rates of 23.33% (7/30), 66.67% (20/30), 93.33% (28/30), and 100.00% (30/30), respectively (χ² = 52.72, P < 0.001). De novo transcriptome assembly generated 176 942 unigenes, with annotation rates of 0.98%, 13.49%, 26.46%, 12.48%, and 14.39% against GO knowledgebase, KEGG pathway database, NR protein sequences database, Pfam database, and Swiss-Prot knowledgebase, respectively. There were 33 up-regulated and 72 down-regulated genes in Group D, 483 up-regulated and 815 down-regulated genes in Group HD, and 245 up-regulated and 172 down-regulated genes in Group LD relative to in Group C. Following removal of overlapping genes across groups and unmatched genes, 11 candidate genes were identified. GO and KEGG analyses revealed 3 heat shock protein (HSP)-related DEGs in these 11 candidate genes, which were annotated as HSP12.2, HSP70, and HSP20 genes and were all significantly up-regulated in each treatment group. Three immune and nervous system-related DEGs were identified, and were all significantly down-regulated in each treatment group, which were involved in the neural cell adhesion molecule L1-like protein pathway, fibrinogen binding protein pathway, and leukocyte elastase inhibitor-like protein pathway. qRT-PCR assay quantified that the expression trends of four genes related to temperature and drought stress across different treatment groups were highly consistent with transcriptome sequencing data. Conclusion The survival rate of B. globosus significantly reduces under combined stresses of extreme temperature and drought, possibly due to an imbalance in its cellular homeostasis regulatory system.

Objective Objective To analyze the potential spatial factors affecting the genetic differentiation of Oncomelania hupensis robertsoni in Yunnan Province. Methods A total of 13 administrative villages were selected from schistosomiasis-endemic areas of Yunnan Province as O. hupensis snail sampling sites. At least 200 snails were collected in each site, and the spatial variable data of each site were recorded, including longitude, latitude and altitude. Thirty active and Schistosoma japonicum uninfected O. hupensis snails were selected from each sampling site by means of the crawling method and the cercarial shedding method. Genomic DNA was extracted from O. hupensis snails. Following PCR amplification, purification of PCR amplification products and sequencing, the gene sequences of O. hupensis snail samples were spliced and edited using the DNAstar software and the NCBI database to yield the complete mitochondrial sequences of O. hupensis snails at each sampling site, and the mitochondrial genetic distance matrix of O. hupensis robertsoni was calculated at each sampling site. The geographical coordinates of each sampling site were marked using the software ArcGIS 10.2, and the straight-line geographical distance between each sampling site was calculated. The altitude difference, longitude difference and latitude difference between each sampling site were calculated using the Excel software, and the correlation between the mitochondrial genetic distance matrix of O. hupensis robertsoni and each spatial variable matrix was examined by using the Mantel test at 13 sampling sites in Yunnan Province. Results Among the 13 O. hupensis snail sampling sites in Yunnan Province, the largest mitochondrial genetic distance of O. hupensis robertsoni snail populations was seen between Anding Village, Nanjian Yi Autonomous County and Caizhuang Village, Midu County (26.244 2), and the largest geographical distance was seen between Dongyuan Village, Gucheng District and Cangling Village, Chuxiong County (272.64 km). The highest altitude difference was seen between Anding Village, Nanjian Yi Autonomous County and Dongyuan Village, Gucheng District (1 086.10 m), and the largest longitude difference was found between Qiandian Village, Eryuan County and Cangling Village, Chuxiong County (1.86°), while the largest latitude difference was measured between Leqiu Village, Nanjian Yi Autonomous County and Dongyuan Village, Gucheng District (1.81°). In addition, the mitochondrial genetic distance of O. hupensis robertsoni snail populations was positively correlated with altitude at 13 snail sampling sites in Yunnan Province (r = 0.542 8, P < 0.001), and showed no significant correlations with geographical distance (r = 0.093 4, P > 0.05), longitude (r = −0.199 5, P > 0.05) or latitude (r = 0.205 7, P > 0.05). Conclusion Altitude may be a potential spatial factor affecting the genetic differentiation of O. hupensis robertsoni in Yunnan Province.

ObjectiveTo delineate imaging findings using an imaging platform and investigate the correlation between MRI characteristics of spinal cord atrophy and clinical diagnosis in children with spinal cord injury (SCI). MethodsImaging data of 150 children with SCI admitted to Beijing Bo'ai Hospital, China Rehabilitation Research Center, from January, 2002 to March, 2024 were collected and imported into the imaging platform. The anteroposterior and transverse diameters of the middle part of the spinal cord at the cross-section with the most severe atrophy were measured, and the relevant indicators of the previous normal spinal cord segment were measured as controls; the radiomic features were extracted. Clinical data of the children including gender, age, cause of injury, sensory level, motor level, spinal cord injury level, injury severity and disease course were collected. ResultsSpinal cord atrophy was identified in 81 cases (54%), among which 78 cases (96%) were American Spinal Injury Association Impairment Scale (AIS) grade A and 3 cases (4%) were AIS grade C. The upper boundary of the spinal cord atrophy site strongly correlated with the injury level, motor level and sensory level (r > 0.8, P < 0.001). ConclusionMore than half of children with SCI may develop secondary spinal cord atrophy, the vast majority of whom suffer from complete spinal cord injury; the upper boundary of spinal cord atrophy is correlated with the injury level.

Organoid-on-a-chip technology represents a promising interdisciplinary advancement that merges two cutting-edge biomedical platforms: stem cell-derived organoids and microfluidics-based organ-on-a-chip systems. Organoids are self-organizing three-dimensional (3D) cell cultures that mimic the key structural and functional features of in vivo organs. However, traditional organoid culture systems are often static, lacking dynamic environmental cues and suffering from limitations such as batch-to-batch variability, low stability, and low throughput. Organ-on-a-chip platforms, by contrast, utilize microfluidic technologies to simulate the dynamic physiological microenvironment of human tissues and organs, enabling more controlled cell growth and differentiation. By integrating the advantages of organoids and organ-on-a-chip technologies, organoid-on-a-chip systems transcend the limitations of conventional 3D culture models, offering a more physiologically relevant and controllable in vitro platform. In organoid-on-a-chip systems, stem cells or pre-formed organoids are cultured in micro-engineered environments that mimic in vivo conditions, enabling precise control over fluid flow, mechanical forces, and biochemical cues. Specifically, these platforms employ advanced strategies including bio-inspired 3D scaffolds for structural support, precise spatial cell patterning via 3D bioprinting, and integrated biosensors for real-time monitoring of metabolic activities. These synergistic elements recreate complex extracellular matrix signals and ensure high structural fidelity. Based on structural complexity, organoid-on-a-chip systems are classified into single-organoid and multi-organoid types, forming a trajectory from unit biomimicry to systemic simulation. Single-organoid chips focus on highly biomimetic units by integrating vascular, immune, or neural functions. Multi-organoid chips simulate inter-organ crosstalk and systemic homeostasis, advancing complex disease modeling and PK/PD evaluation. This emerging technology has demonstrated broad application potential in multiple fields of biomedicine. Organoid-on-a-chip systems can recapitulate organ developmentin vitro, facilitating research in developmental biology. They mimic organ-specific physiological activities and mechanisms, showing promising applications in regenerative medicine for tissue repair or replacement. In disease modeling, they support the reconstruction of models for neurodegenerative, inflammatory, infectious, metabolic diseases, and cancers. These platforms also enable in vitro drug testing and pharmacokinetic studies (ADME). Patient-derived chips preserve genetic and pathological features, offering potential for precision medicine. Additionally, they reduce species differences in toxicology, providing human-relevant data for environmental, food, cosmetic, and drug safety assessments. Despite progress, organoid-on-a-chip systems face challenges in dynamic simulation, extracellular matrix (ECM) variability, and limited real-time 3D imaging, requiring improved materials and the integration of developmental signals. Current bottlenecks also include the high technical threshold for automation and the lack of standardized validation frameworks for regulatory adoption. Meanwhile, the concept of a “human-on-a-chip” has been proposed to mimic whole-body physiology by integrating multiple organoid modules. This approach enables systemic modeling of drug responses and toxicity, with the potential to reduce animal testing and revolutionize drug development. Future advancements in bio-responsive hydrogels and flexible biosensors will further empower these platforms to bridge the gap between bench-side research and personalized clinical interventions. In conclusion, organoid-on-a-chip technology offers a transformative in vitro model that closely recapitulates the complexity of human tissues and organ systems. It provides an unprecedented platform for advancing biomedical research, clinical translation, and pharmaceutical innovation. Continued development in biomaterials, microengineering, and analytical technologies will be essential to unlocking the full potential of this powerful tool.

ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.

ObjectiveThis Study was conducted to investigate the interaction mechemisms between gutmicrobiota dysregulation and ischemic stroke by establishing a rat model of ischemic stroke and employing fecal microbiota transplantation (FMT). MethodsA preliminary experiment was conducted to establish an antibiotic-induced pseudo-sterile (ABX) rat model through antibiotic treatment, and a cerebral ischemia model was prepared using the middle cerebral artery occlusion (MCAO) method. Fecal microbiota from stroke patients and healthy individuals were transplanted via FMT, followed by behavioral testing. 16S rRNA sequencing was used to analyze the microbial community, hematoxylin and eosin (HE) staining to observe histopathological status, transmission electron microscopy (TEM) to examine the tight junction structure of the small intestine, and enzyme-linked immunosorbent assay (ELISA) to detect levels of inflammatory factors and intestinal barrier-related markers. Results16S rRNA sequencing of fecal samples showed that compared with the normal control group and the metronidazole group, the abundance and diversity of fecal microorganisms in the quadruple antibiotic group were significantly reduced, indicating successful establishment of the ABX model. After transplanting fecal microbiota from stroke patients into ABX rats, significant changes in gut microbiota composition were observed. Behavioral tests revealed that the MCAO model group showed significant decreases in both horizontal movement and vertical exploration abilities. ELISA results indicated that IL-17 concentration in the ABX+mFMT (antibiotic-treated+model fecal microbiota transplantation) group was lower than in the ABX+cFMT (antibiotic-treated+control fecal microbiota transplantation) group, suggesting that IL-17 may serve as a key inflammatory indicator for evaluating the impact of stroke intervention on gut microbiota. Triphenyltetrazolium chloricle staining (TTC) staining suggested that gut microbiota intervention may increase the risk of stroke. HE staining showed that, except for the control group, all groups exhibited ischemic changes and inflammatory infiltration in brain tissues. TEM revealed that microvilli of small intestinal epithelial cells in the ABX+mFMT group were sparser than those in the ABX+cFMT group, indicating that microbial intervention affects intestinal barrier function. ConclusionThe ABX model established using broad-spectrum antibiotics showed no significant differences in physiological characteristics compared to normal rats, and the findings were consistent with those from germ-free rat models. Stroke prognosis appears to be influenced by intestinal dysbiosis, accompanied by significantly elevated levels of the pro-inflammatory cytokine IL-17, which may exacerbate neural injury via the gut-brain axis. Behavioral experiments indicated that transplantation of gut microbiota from stroke rats impaired cognitive function. Furthermore, IL-17 demonstrated sensitivity to alterations in the gut microbiota, suggesting its potential as a key therapeutic target for stroke intervention.

ObjectiveCleft palate (CP) is a common congenital deformity often associated with velopharyngeal insufficiency (VPI), which disrupts the physiological coupling between respiration and speech. Conventional clinical assessments, such as nasometry and spirometry, provide limited static data and fail to visualize the dynamic spatiotemporal distribution of lung ventilation during phonation. This study introduces spatiotemporal electrical impedance tomography (ST-EIT) to evaluate speech-respiratory functional features in CP patients compared to normal controls (NC). The aim is to characterize multi-domain respiratory patterns and to validate an interpretable machine learning framework for providing objective, quantitative evidence for clinical assessment. MethodsSeventy-five participants were enrolled in this study, comprising 37 patients with surgically repaired CP and 38 healthy volunteers matched for age, gender, and body mass index (BMI). All subjects performed standardized sustained phonation tasks while undergoing synchronous monitoring with a 16-electrode EIT system and a pneumotachograph. A comprehensive feature engineering pipeline was developed to extract physiological parameters across 3 complementary domains. (1) Temporal domain: including inspiratory/expiratory phase duration (tPhase), time constants (Tau), and inspiratory-to-expiratory time ratios (TI/TE); (2) airflow domain: comprising mean flow, peak flow, and instantaneous flow at 25%, 50%, and 75% of tidal volume; and (3) spatial domain: quantifying global and regional tidal impedance variation (TIV), global inhomogeneity (GI), and center of ventilation (CoV). Extreme Gradient Boosting (XGBoost) classifiers were trained using 5 distinct data sources (Spirometry, Nasometry, Inspiratory-EIT, Expiratory-EIT, and fused ST-EIT). Model performance was rigorously evaluated via stratified 5-fold cross-validation, and Shapley additive explanations (SHAP) were employed to quantify global and local feature contributions. ResultsThe CP group exhibited a distinct respiratory phenotype compared to controls. In the temporal domain, CP patients showed significantly shorter inspiratory (1.60 s vs.1.85 s, P<0.001) and expiratory phase durations (2.45 s vs. 3.95 s, P<0.001), indicating a rapid, shallow breathing rhythm. In the airflow domain, while inspiratory flows were comparable, the CP group demonstrated significantly elevated mean and peak flows during the expiratory phase (P<0.001), reflecting compensatory respiratory effort. Spatially, CP patients presented significant ventilation redistribution, characterized by higher regional TIV in the right-anterior (ROI1) and left-posterior (ROI4) quadrants, but lower TIV in the left-anterior (ROI2) quadrant. In terms of diagnostic accuracy, the multi-modal ST-EIT model achieved the highest performance (AUC: 0.915±0.012, Accuracy: 0.843±0.019, F1-score: 0.872±0.017), substantially outperforming models based on spirometry (AUC: 0.721) or nasometry (AUC: 0.625) alone. Interpretability analysis revealed that spatial domain features were the most critical, contributing 53.4% to the model’s decision-making, followed by temporal (25.0%) and airflow (21.6%) features. ConclusionST-EIT successfully captures the temporal, airflow, and spatial deviations in CP speech respiration that are undetectable by conventional methods—specifically, rapid phase transitions, hyperdynamic expiratory airflow, and regional ventilation heterogeneity. This study validates ST-EIT as a robust, non-invasive, and radiation-free tool for characterizing speech-respiratory dysfunction, offering high clinical value for bedside screening, rehabilitation planning, and longitudinal monitoring of patients with cleft palate.

Organoid-on-a-chip technology represents a promising interdisciplinary advancement that merges two cutting-edge biomedical platforms: stem cell-derived organoids and microfluidics-based organ-on-a-chip systems. Organoids are self-organizing three-dimensional (3D) cell cultures that mimic the key structural and functional features of in vivo organs. However, traditional organoid culture systems are often static, lacking dynamic environmental cues and suffering from limitations such as batch-to-batch variability, low stability, and low throughput. Organ-on-a-chip platforms, by contrast, utilize microfluidic technologies to simulate the dynamic physiological microenvironment of human tissues and organs, enabling more controlled cell growth and differentiation. By integrating the advantages of organoids and organ-on-a-chip technologies, organoid-on-a-chip systems transcend the limitations of conventional 3D culture models, offering a more physiologically relevant and controllable in vitro platform. In organoid-on-a-chip systems, stem cells or pre-formed organoids are cultured in micro-engineered environments that mimic in vivo conditions, enabling precise control over fluid flow, mechanical forces, and biochemical cues. Specifically, these platforms employ advanced strategies including bio-inspired 3D scaffolds for structural support, precise spatial cell patterning via 3D bioprinting, and integrated biosensors for real-time monitoring of metabolic activities. These synergistic elements recreate complex extracellular matrix signals and ensure high structural fidelity. Based on structural complexity, organoid-on-a-chip systems are classified into single-organoid and multi-organoid types, forming a trajectory from unit biomimicry to systemic simulation. Single-organoid chips focus on highly biomimetic units by integrating vascular, immune, or neural functions. Multi-organoid chips simulate inter-organ crosstalk and systemic homeostasis, advancing complex disease modeling and PK/PD evaluation. This emerging technology has demonstrated broad application potential in multiple fields of biomedicine. Organoid-on-a-chip systems can recapitulate organ developmentin vitro, facilitating research in developmental biology. They mimic organ-specific physiological activities and mechanisms, showing promising applications in regenerative medicine for tissue repair or replacement. In disease modeling, they support the reconstruction of models for neurodegenerative, inflammatory, infectious, metabolic diseases, and cancers. These platforms also enable in vitro drug testing and pharmacokinetic studies (ADME). Patient-derived chips preserve genetic and pathological features, offering potential for precision medicine. Additionally, they reduce species differences in toxicology, providing human-relevant data for environmental, food, cosmetic, and drug safety assessments. Despite progress, organoid-on-a-chip systems face challenges in dynamic simulation, extracellular matrix (ECM) variability, and limited real-time 3D imaging, requiring improved materials and the integration of developmental signals. Current bottlenecks also include the high technical threshold for automation and the lack of standardized validation frameworks for regulatory adoption. Meanwhile, the concept of a “human-on-a-chip” has been proposed to mimic whole-body physiology by integrating multiple organoid modules. This approach enables systemic modeling of drug responses and toxicity, with the potential to reduce animal testing and revolutionize drug development. Future advancements in bio-responsive hydrogels and flexible biosensors will further empower these platforms to bridge the gap between bench-side research and personalized clinical interventions. In conclusion, organoid-on-a-chip technology offers a transformative in vitro model that closely recapitulates the complexity of human tissues and organ systems. It provides an unprecedented platform for advancing biomedical research, clinical translation, and pharmaceutical innovation. Continued development in biomaterials, microengineering, and analytical technologies will be essential to unlocking the full potential of this powerful tool.

ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.