1.Spatiotemporal Electrical Impedance Tomography for Speech Respiratory Assessment in Cleft Palate: an Interpretable Machine Learning Study
Yang WU ; Xiao-Jing ZHANG ; Hao YU ; Cheng-Hui JIANG ; Bo SUN ; Jia-Feng YAO
Progress in Biochemistry and Biophysics 2026;53(2):485-500
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.
2.Spatiotemporal Electrical Impedance Tomography for Speech Respiratory Assessment in Cleft Palate: an Interpretable Machine Learning Study
Yang WU ; Xiao-Jing ZHANG ; Hao YU ; Cheng-Hui JIANG ; Bo SUN ; Jia-Feng YAO
Progress in Biochemistry and Biophysics 2026;53(2):485-500
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.
3.Construction and Application of a Real-World Cohort of Community-Acquired Pneumonia Based on a Multimodal Large-Scale Traditional Chinese Medicine Big Data Platform
Zhichao WANG ; Xianmei ZHOU ; Fanchao FENG ; Mengqi WANG ; Xin WANG ; Bin KANG ; Xiaofan YU ; Xiaoxiao WANG ; Lei XIAO ; Juan LI ; Zhichao ZHANG ; Ye MA ; Yeqing JI ; Xin TONG ; Zhuoyue WU ; Jia LIU
Journal of Traditional Chinese Medicine 2026;67(9):961-965
This paper introduces a real-world cohort research model for community-acquired pneumonia (CAP) based on the Jiangsu Traditional Chinese Medicine (TCM) Dominant Diseases Diagnosis and Treatment Data Platform. Firstly, data cleaning is performed by standardizing diagnosis, symptoms, treatment and imaging, intelligently extracting unstructured information, and cleaning and constructing a standardized database. Secondly, for cohort establishment, CAP patients across the province are screened in accordance with CAP diagnostic criteria to build a high-quality disease-specific cohort. Lastly, in terms of protocol design, the characteristics of TCM research and the CAP disease profile are considered to determine appropriate inclusion and exclusion criteria, estimate sample size, define interventions, outcomes and economic evaluations, providing a reference for real-world TCM research on CAP.
4.Construction and Application of a Real-World Cohort of Community-Acquired Pneumonia Based on a Multimodal Large-Scale Traditional Chinese Medicine Big Data Platform
Zhichao WANG ; Xianmei ZHOU ; Fanchao FENG ; Mengqi WANG ; Xin WANG ; Bin KANG ; Xiaofan YU ; Xiaoxiao WANG ; Lei XIAO ; Juan LI ; Zhichao ZHANG ; Ye MA ; Yeqing JI ; Xin TONG ; Zhuoyue WU ; Jia LIU
Journal of Traditional Chinese Medicine 2026;67(9):961-965
This paper introduces a real-world cohort research model for community-acquired pneumonia (CAP) based on the Jiangsu Traditional Chinese Medicine (TCM) Dominant Diseases Diagnosis and Treatment Data Platform. Firstly, data cleaning is performed by standardizing diagnosis, symptoms, treatment and imaging, intelligently extracting unstructured information, and cleaning and constructing a standardized database. Secondly, for cohort establishment, CAP patients across the province are screened in accordance with CAP diagnostic criteria to build a high-quality disease-specific cohort. Lastly, in terms of protocol design, the characteristics of TCM research and the CAP disease profile are considered to determine appropriate inclusion and exclusion criteria, estimate sample size, define interventions, outcomes and economic evaluations, providing a reference for real-world TCM research on CAP.
5.Preliminary study on the biological role of EF-hand domain-containing protein 2 in hepatocellular carcinoma
Yanmei ZHANG ; Xiao LI ; Xueqiang JIA ; Juanzi LIU ; Wanqing LI ; Junfeng XUAN ; Shiyu FENG ; Zhaohui SUN ; Weiyun ZHANG
Chinese Journal of Preventive Medicine 2025;59(8):1224-1231
This study investigates the expression pattern and functional significance of EF-hand domain-containing protein 2 (EFHD2) in hepatocellular carcinoma (HCC), with particular focus on its regulatory effects on tumor proliferation, migration, and invasion. Cellular experimental study was completed from June 2024 to January 2025 in the Basic Laboratory of the General Hospital of Southern Theater Command. TCGA database to determine EFHD2 expression and its clinicopathological correlations. GSCA database to assess methylation patterns and immune infiltration. Model of transient overexpression and knockdown of EFHD2 was constructed in hepatocellular carcinoma cells Hep3B, then RT-qPCR and Western blot were applied to verify the transfection efficiency. CCK-8 and colony formation assays for proliferation assessment, Transwell chambers for migration/invasion quantification. Protein-protein interaction networks were constructed via STRING, followed by GO/KEGG enrichment analysis. Statistical analysis was performed using the two independent samples t-test. The results showed that EFHD2 demonstrated significant upregulation in HCC tissues versus normal controls ( P<0.05). Elevated EFHD2 expression correlated with advanced clinical stage ( P<0.05) and poor differentiation ( P<0.05). In the CCK-8 assay, the EFHD2 overexpression group demonstrated significantly higher cell viability than the control group, as evidenced by 450 nm relative absorbance values on Day 1 (0.529±0.019 vs. 0.515±0.016, F=0.041, P=0.320), Day 2 (1.356±0.019 vs. 1.094±0.042, F=3.833, P<0.001), Day 3 (2.817±0.049 vs. 2.143±0.124, F=3.833, P<0.001), and Day 4 (3.848±0.015 vs. 3.430±0.021, F=0.469, P<0.001). The EFHD2 knockdown group showed reduced cell viability compared to controls: Day 1 (0.541±0.020 vs. 0.552±0.015, F=0.098, P=0.423), Day 2 (1.154±0.009 vs. 1.326±0.029, F=2.485, P<0.001), Day 3 (2.453±0.041 vs. 2.653±0.031, F=0.479, P<0.001), and Day 4 (3.685±0.038 vs. 3.836±0.021, F=6.804, P<0.001). In colony formation assays, the overexpression group displayed a significant increase in colony numbers (254.667±23.861 vs. 186.000±16.703, F=0.865, P=0.015), whereas the knockdown group exhibited decreased colony formation (229.000±24.637 vs. 306.667±36.501, F=0.988, P=0.038). In Transwell assays, the EFHD2 overexpression group revealed enhanced migratory capacity [ (605.000±72.670) cells vs. (472.667±28.095) cells, F=2.462, P=0.042] and invasive potential [(767.333±21.221) cells vs. (414.333±16.623) cells, F=0.331, P<0.001]. The knockdown group showed attenuated migration [(311.000±71.084) cells vs. (479.667±50.846) cells, F=0.718, P=0.029] and invasion [(247.667±48.263) cells vs. (345.667±32.130) cells, F=0.727, P=0.043] compared to controls. The network of EFHD2-interacting proteins was further constructed by the STRING database, and the GO and KEGG analysis were used to perform bioinformatics analysis reveal that EFHD2 is mainly involved in actin cytoskeleton regulation. In conclusion, EFHD2 is highly expressed in HCC and is involved in the process of proliferation, migration and invasion of HCC.
6.Association of cadmium internal exposure levels with blood lipid in adults aged 18 to 79 years in China
Haocan SONG ; Saisai JI ; Zheng LI ; Yawei LI ; Feng ZHAO ; Yingli QU ; Yifu LU ; Yingying HAN ; Junxin LIU ; Jiayi CAI ; Tian QIU ; Wenli ZHANG ; Xiao LIN ; Junfang CAI ; Yuebin LYU ; Xiaoming SHI
Chinese Journal of Preventive Medicine 2025;59(8):1254-1263
Objective:To explore the association of blood and urinary cadmium levels with lipid profile levels and dyslipidemia in Chinese adults aged 18 to 79 years.Methods:Based on the China National Human Biomonitoring (CNHBM) program, a cross-sectional survey was conducted from 2017 to 2018 using a multi-stage stratified random sampling method, including a total of 10 713 adults aged 18 to 79 years. Data was obtained through questionnaires, physical examinations, biological sample collection, and laboratory testing. Multiple linear mixed effect model (MLMM) and generalized linear mixed effect model (GLMM) were used to analyze the association of blood and creatinine-corrected urinary cadmium levels with lipid profile levels as well as dyslipidemia among adults.Results:The age of 10 713 participants was (47.23±0.24) years, with 5 372 males accounting for 61.3% of the national population. The weighted mean±standard error (SE) of total cholesterol (TC), triglycerides (TG), low-density lipoprotein cholesterol (LDL-C), and high-density lipoprotein cholesterol (HDL-C) was (5.21±0.03), (1.86±0.03), (2.96±0.03), and (1.43±0.01) mmol/L, respectively. The prevalence rate of hypercholesterolemia, hypertriglyceridemia, mixed hyperlipidemia, low HDL-C, and high LDL-C was 16.0%, 21.6%, 6.6%, 13.5%, and 10.0%, respectively. MLMM showed that, after adjusting for relevant confounders, log-transformed blood cadmium levels were positively associated with increased levels of TC, TG and LDL-C ( P<0.05). When blood cadmium levels were categorized into quartiles, compared to the lowest exposure group ( Q1), participants in the highest blood cadmium exposure group ( Q4) had increases of 0.19 (95% CI: 0.06, 0.32) mmol/L in TC and 0.25 (95% CI: 0.08, 0.43) mmol/L in TG. GLMM indicated that, after adjusting for confounders, higher blood cadmium exposure levels were associated with increased risks of hypercholesterolemia, hypertriglyceridemia, mixed hyperlipidemia, and high LDL-C ( P<0.05). Further analysis by quartiles showed that, compared to the blood cadmium Q1 exposure group, the OR value (95% CI) for the Q4 group was 1.53 (1.12, 2.08) for hypercholesterolemia, 1.54 (1.09, 2.17) for hypertriglyceridemia, 2.24 (1.47, 3.40) for mixed hyperlipidemia, and 1.49 (1.07, 2.09) for high LDL-C. Conclusion:The cadmium internal exposure levels are associated with blood lipid profile levels as well as the incidence of dyslipidemia in Chinese adults aged 18 to 79.
7.Construction and application of the evaluation system for the risk of measles transmission in Henan Province
Binghui DU ; Xiaoxiao ZHANG ; Daxing FENG ; Wenhui WANG ; Mingyu ZHANG ; Zhanpei XIAO ; Yanyang ZHANG
Chinese Journal of Preventive Medicine 2025;59(8):1279-1286
Objective:To construct a risk assessment system for measles transmission in Henan Province and scientifically evaluate the risk levels of measles transmission in each city in Henan Province.Methods:The modified Delphi method was used to conduct two rounds of expert consultations to construct a risk assessment system for measles transmission. Data from all 191 cities, counties and districts in Henan Province were collected. The internal consistency (Cronbach′s α), content validity (content validity index, CVI) and structural validity (factor analysis) of the indicator system were evaluated to optimize the assessment framework. The indicator assignment method was adopted, and the comprehensive risk scores were obtained by adding the scores according to different weights. Results:Both rounds of consultation witnessed a 100% participation rate among all experts. The authority coefficients of experts in the two rounds were 0.920 and 0.925, and concordance coefficients were 0.201 ( χ 2=161.11, P<0.001) and 0.210 ( χ 2=163.80, P<0.001). The constructed assessment system comprised four dimensions—population immunity levels, surveillance quality, importation risk, and technical reserve of emergency response capacities—with a total of 30 indicators. Reliability analysis of the assessment system showed an overall Cronbach′s α of 0.741. Validity analysis revealed that all content validity indices reached 1.000, with principal factors cumulatively accounting for 67.625% of the variance, and all factor loadings exceeded 0.400. The measles transmission risk assessment in Henan Province using this assessment system identified Zhengzhou (92), Xinxiang (91), Xinyang (89), and Pingdingshan (73) as high-risk regions. Conclusion:The risk assessment system developed in this study demonstrates good reliability and validity, effectively reflecting measles transmission risks across Henan Province. The findings highlight the need to strengthen surveillance and control measures in high-risk areas, particularly in Zhengzhou.
8.Explore the feasibility of subsegmentectomy in clinical application
Fei QI ; Hongxiang FENG ; Yu HAN ; Fei XIAO ; Yuhui SHI ; Chaoyang LIANG ; Deruo LIU ; Fanjia KONG ; Zhenrong ZHANG
Chinese Journal of Thoracic and Cardiovascular Surgery 2025;41(1):10-15
Objective:Analyze the imaging and clinical data of cases undergoing thoracoscopic subsegmental resection for non-small cell lung cancer(NSCLC) with a diameter≤2 cm, and explore the clinical outcomes of subsegmental resection.Methods:A retrospective analysis was conducted on the clinical data of 58 patients who underwent thoracoscopic subsegmentectomy in China-Japan Friendship Hospital from January 2020 to July 2024. Three-dimensional reconstruction technology was used for surgical planning before the operation, and thoracoscopic subsegmentectomy was performed, including single lung subsegmentectomy(Group 1), multiple lung subsegmentectomy(Group 2), and combined segmentectomy and subsegmentectomy(Group 3).Results:All patients successfully completed the surgery, with 23 cases of single lung subsegmentectomy, 6 cases of multiple lung subsegmentectomy, and 29 cases of combined segmentectomy and subsegmentectomy. The median intraoperative blood loss was 30.0(20.0, 30.0)ml, the average operation time was(2.03±0.68) h, the average pathological size of the nodules was(10.53±4.45) mm, and the average postoperative tube retention was(2.55±0.92) days. There were 6 cases of postoperative complications, including pulmonary air leakage in 2 cases, cerebral embolism in 1 case, pulmonary embolism in 1 case, pulmonary infection in 1 case, and atrial fibrillation in 1 case. All patients had negative surgical margins in the postoperative pathology. Group 1 had less average intraoperative blood loss than Group 2, with statistically significant differences( P=0.027). Surgical procedures for the upper lobe of the lung mainly involve the resection of combined segments and subsegments, while those for the lower lobe primarily consist of single segmentectomy. Conclusion:Subsegmentectomy is an effective surgical approach when the nodule is small and a clear margin can be ensured, allowing for better preservation of remaining lung tissue. Bleeding during multiple subsegmentectomies is greater than that in single subsegmentectomy and combined segmentectomy with subsegmentectomy, which may be related to the more complex vascular variations in multiple subsegmentectomies.
9.Establishment of preparation process and quality standard for Zhenggu Pills
Wen-ming ZHANG ; Zi-fang FENG ; Li-hong GU ; Ping QIN ; Zhen-hua BIAN ; Min-min HU ; Xiao-wei CHEN
Chinese Traditional Patent Medicine 2025;47(9):2863-2869
AIM To establish the preparation process and quality standard for Zhenggu Pills.METHODS With decoction time,decoction frequency and water addition as influencing factors,comprehensive score for extract yield and transfer rates of epicatechin and naringin as an evaluation index,the decoction process was optimized by orthogonal test.With sugarless paste relative density,medicinal powder fineness,sugarless paste-corn starch ratio,drying temperature and drying time as influencing factors,soft material traits,pill formability,moisture and disintegration time limit as evaluation indices,the formability process was optimized by single factor test.TLC was adopted in the qualitative identification of Dipsaci Radix,salt-processed Psoraleae Fructus,cooked Rhei Radix et Rhizoma and Notoginseng Radix et Rhizoma.HPLC was used for the content determination of paeoniflorin and naringin.RESULTS The optimal decoction process was determined to be 0.5 h for decoction time,two times for decoction frequency,and 10 times for water addition,the comprehensive score was 0.93.The optimal formability process was determined to be 1.21-1.22 for sugarless paste relative density,80 mesh for medicinal powder fineness,1∶0.17-1∶0.18 for sugarless paste-corn starch ratio,70 ℃ for drying temperature,and 24 h for drying time,good soft material traits and pill formability were observable,and moisture and disintegration time limit accored with 2020 edition of Chinese Pharmacopoeia requirements.The TLC spots were clear without negative interference.Two constituents showed good linear relationships within 61.30-490.41 μg/mL(r=0.999 8)and 3.27-26.18 μg/mL(r=0.999 8),whose average recoveries were 100.15%and 98.15%with the RSDs of 0.55%and 2.30%,respectively.CONCLUSION This stable,reliable and specific method can be used for the production and quality evaluation of Zhenggu Pills.
10.Microbial community mediated by microbial agents improves the quality of Epimedium pubescens Maxim.
Lai KUNYANG ; Wan XIUFU ; Xiao JIANCAI ; Wang HONGYANG ; Shi SHANGXUAN ; Yan BINBIN ; Lyu CHAOGENG ; Zhang CHENGCAI ; Zhang YUFEI ; Yuan FENG ; Zhao ZHE ; Zhu SHOUDONG ; Kang CHUANZHI ; Zhang YAN
Science of Traditional Chinese Medicine 2025;3(3):270-281
Background:Optimizing cultivation techniques for traditional Chinese medicine has become a crucial means to improve the quality of medicinal materials.Microbial agents,as environmentally friendly and efficient plant growth promoters and soil conditioners,have increasingly attracted attention in eco-agriculture research.Objective:Our understanding remains limited regarding how the application of microbial agents,alone or in combination,affects changes in the rhizosphere microbiome and its association with the bioactive components of medicinal materials.Methods:In this study,Epimedium pubescens Maxim.was employed as a model plant to examine the effects of 2 microbial agents(Paenibacillus mucilaginosus and Bacillus subtilis)applied individually and in combination on plant growth and the accumulation of bioactive components.Additionally,this study explored the relationship between the rhizosphere microbiome and plant development.Results:The application of microbial agents increased the yield of E.pubescens leaves by 20.30%to 33.66%and enhanced the total flavonol glycosides content by 11.40%to 29.94%.Meanwhile,microbial treatments reshaped the rhizosphere microbiome,promoted the enrichment of beneficial microorganisms(e.g.,Frankia and Paenibacillus),suppressed phytopathogenic fungi such as Didymella and Scytalidium,and enhanced the stability of the soil microbial co-occurrence network.The partial least squares path model suggested that microbial agents not only directly impact the quality of medicinal herbs but also indirectly alter the accumula-tion of bioactive components by modulating the soil microbiome.Conclusion:These findings deepen our understanding of the relationship between medicinal plant quality and rhizosphere micro-biomes as mediated by microbial agents.They also provide a basis for designing and manipulating synthetic microbial communities to promote sustainable development in eco-agriculture.

Result Analysis
Print
Save
E-mail