1.Identification and Analysis of bHLH Genes Related to Color Formation of Gastrodia elata Stem
Xue JIANG ; Dandan RAN ; Xiuwen WANG ; Xiaobo ZHANG ; Xiaohong OU ; Jie PAN ; Tao ZHOU ; Zhen OUYANG ; Jiao XU
Chinese Journal of Experimental Traditional Medical Formulae 2026;32(8):202-209
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.
2.Identification and Analysis of bHLH Genes Related to Color Formation of Gastrodia elata Stem
Xue JIANG ; Dandan RAN ; Xiuwen WANG ; Xiaobo ZHANG ; Xiaohong OU ; Jie PAN ; Tao ZHOU ; Zhen OUYANG ; Jiao XU
Chinese Journal of Experimental Traditional Medical Formulae 2026;32(8):202-209
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.
3.Clinicopathological analysis of conjunctival lymphangiectasia
Jie RAN ; Duosheng XIA ; Fuli WANG ; Xuandi SU ; Lele ZHENG ; Yujin WANG
International Eye Science 2026;26(7):1276-1280
AIM: To analyze the pathological features, immunophenotype, and imaging findings of conjunctival lymphangiectasia(CL), and to explore the etiological mechanisms and provide a theoretical basis for clinical diagnosis and treatment. METHODS:This single-center descriptive cross-sectional study enrolled postoperative specimens from patients with CL who underwent surgical treatment in the hospital between Feb. 2023 and Sept. 2025. Routine hematoxylin and eosin(HE)staining and immunohistochemical staining(D2-40, CD31, CD34, CK)were performed. Anterior segment optical coherence tomography(AS-OCT)was used to observe the lesion morphology. The pathological results were comprehensively analyzed combined with clinical data. RESULTS: A total of postoperative specimens from 32 eyes of 32 patients with CL were enrolled, including 23 females(72%)and 9 males(28%), with a mean age of 53.03±12.47 y. All patients presented with single or multiple transparent cystic elevations beneath the bulbar conjunctiva. The postoperative pathological manifestations were characterized by dilation of conjunctival lymphatic vessels lined with a single layer of flattened endothelial cells, accompanied by edema and inflammatory infiltration in the surrounding stroma. All cases were positive for D2-40, confirming a lymphatic origin; some cases also expressed CD31 and CD34. AS-OCT revealed the lesions as unilocular or multilocular cystic spaces with low reflectivity. After complete surgical resection, the mean follow-up period was 16.2 mo with no recurrence.CONCLUSION:CL is a benign ocular surface lesion characterized by lymphatic vessel dilation. The endothelium co-expresses lymphatic and some vascular markers, suggesting that CL may belong to the spectrum of vascular malformations. AS-OCT has adjunctive diagnostic value, and surgical resection has definitive therapeutic efficacy.
4.Mosquito monitoring and influence of meteorological factors on mosquito density around Beijing Capital International Airport, 2016-2022
Ran FENG ; Tie-zheng MA ; Si-jie ZHU ; Bo TIAN ; Fei QUAN ; Zhi-lin WU ; Xiao-tao LIU ; Fu-yuan ZHANG ; Song-jian ZHANG
Acta Parasitologica et Medica Entomologica Sinica 2026;33(1):40-48
Objective To understand mosquito density, species composition, and seasonal variations around the Beijing Capital International Airport(BCIA)and provide scientific evidence for mosquito-borne disease prevention and control. Methods Meteorological data were collected from 2016 to 2022, and mosquito density was monitored at seven surveillance sites around the BCIA from May to October each year using CO2-baited mosquito traps to analyze the relationship between mosquito density and meteorological factors. Results In total,68 518 female mosquitoes were captured, with a mosquito density of 10.20 per light · hour. The dominant species was Culex pipiens pallens(73.90%), followed by Aedes albopictus(9.71%)and Ae. vexans(9.29%);the mosquito density was highest in 2016(14.8 per light·hour)and lowest in 2022(4.46 per light·hour). Significant statistical differences were observed in mosquito density among different species(F=18.118, P<0.05). The density of Cx. pipiens pallens and other mosquito species showed a mutually exclusive trend; the peak of mosquito density varied based on years and habitats, with significant statistical differences in mosquito density among different habitats(F=8.504, P<0.05). The highest mosquito density was observed near the Wenyu River(21.50 per light·hour), whereas the lowest was recorded at the BCIA construction site(2.44 per light· hour). The monthly average temperature(r=0.595), monthly average highest temperature(r=0.575), and monthly average lowest temperature(r=0.624)showed moderate positive correlations, whereas the monthly average air pressure(r=-0.484)showed a moderate negative correlation. The monthly average minimum temperature was included in the regression model(F=25.575, P<0.000), and the equation Y=1.029X4-8.181 was used. The monthly average air pressure 2 months prior(b=-2.418, β=-1.619, P<0.05)and the monthly average relative humidity(b=-0.739, β=-1.201, P<0.05)significantly negatively predicted the mosquito density. Notably, the regression equation used was y=2526.170+(-2.418)X9+(-0.739)X13. The exposure-response analysis revealed that the density of mosquitoes was not linearly related to the average monthly air pressure, rainfall, and duration of sunshine. Conclusions The dominant mosquito species around the BCIA was Cx. pipiens pallens, Ae. albopictus, and Ae. vexans. During the control of Cx. pipiens pallens, the impact on other mosquito species must also be considered. Environmental management, breeding sites, and scientific use of pesticides should be prioritized based on the activity periods of mosquitoes. Further studies on mosquitoes and meteorological factors should be conducted to provide references and novel avenues for mosquito control.
5.The Influence of BDNF on Cortical Remodeling After Peripheral Nerve Injury
Bo-Yuan CHEN ; Hao-Ran CHEN ; Tao ZHANG ; Jie ZHANG
Progress in Biochemistry and Biophysics 2026;53(8):2194-2209
Peripheral nerve injury (PNI) severs peripheral connections and induces profound reorganization of primary somatosensory (S1) and motor (M1) cortical maps, a dynamic process that critically shapes the extent and quality of functional recovery. Brain-derived neurotrophic factor (BDNF) has emerged as a central molecular hub linking the initial peripheral insult to subsequent central plastic changes. This review systematically examines BDNF-mediated cortical remodeling following PNI, with a particular emphasis on its spatiotemporally specific regulatory mechanisms and the therapeutic opportunities they present. Immediately after PNI, deafferented cortical territories undergo rapid functional silencing and, over time, become progressively encroached upon by adjacent intact representations. This maladaptive reorganization—characterized by sensory map fusion and motor compensatory invasion—actively impedes successful reinnervation and optimal recovery. Such phenomena have been consistently documented in both non-human primate models and human fMRI studies, spanning sensory and motor modalities. Motor remapping is also prominently observed in facial nerve injury and chronic nerve compression models, underscoring the universality of these plastic changes. Notably, a range of interventions including regional local anesthesia, auditory-tactile substitution, and guided tactile imagery can effectively mitigate this aberrant remodeling, revealing promising avenues for the therapeutic modulation of cortical plasticity during critical post-injury windows. BDNF coordinates post-PNI cortical remodeling through three mechanistically interconnected pathways. First, BDNF drives dendritic and axonal growth via TrkB receptor signaling: retrogradely transported axonal signaling endosomes integrate CREB-dependent transcriptional programs with mTOR-mediated local protein synthesis to bolster dendritic arborization and axonal elongation, while BDNF-induced Limk1 translation precisely fine-tunes actin cytoskeletal dynamics, thereby stabilizing dendritic spine morphology. Second, BDNF governs GABAergic interneuron maturation and perineuronal net formation through the JNK signaling cascade, thereby acting as a “plasticity brake” that preserves the delicate excitatory/inhibitory balance and constrains indiscriminate, maladaptive rewiring. Third, BDNF robustly potentiates NMDAR-dependent long-term potentiation (LTP) by elevating NMDA channel open probability and increasing synaptic receptor density, thus tightly coupling structural remodeling with functional synaptic strengthening. Concurrently, long-term depression pathways are subject to modulation, with NMDAR antagonism exerting time-dependent, bidirectional effects that vary with the post-injury phase. Superimposed on these mechanisms, the functionally antagonistic actions of proBDNF/p75NTR (favoring pruning and growth cone collapse) and mature BDNF/TrkB (promoting survival and stabilization) add further layers of regulatory complexity, highlighting the absolute necessity for precise spatiotemporal control in any therapeutic strategy. Extending our previous experimental findings, we introduce a “cortical reclamation threshold” model. In this framework, early upregulation of BDNF-TrkB signaling enhances inhibitory tone to restrict cortical encroachment by neighboring intact regions, whereas delayed but precisely timed BDNF delivery following surgical nerve repair actively facilitates the reclamation of the original, deafferented cortical maps by regenerating afferent fibers. Importantly, the common BDNF Val66Met polymorphism may significantly modulate individual thresholds, providing a biologically grounded rationale for patient stratification and personalized timing of intervention. Despite this mechanistic promise, clinical translation faces considerable hurdles. These include inherently poor pharmacokinetics of BDNF, severely limited penetration across the blood-brain barrier, dose-dependent neurotoxicity mediated through the p75NTR receptor, and the overarching dual-edged nature of BDNF signaling, which can simultaneously promote both adaptive and maladaptive plasticity. Emerging delivery platforms——including engineered exosomes, polymeric nanoparticles, cryogel microcarriers, and adeno-associated virus (AAV)-based gene therapy——offer potential solutions, but all require rigorous pharmacokinetic, toxicological, and functional validation. Future research priorities should encompass: (1) comprehensive single-cell and spatial transcriptomic mapping to resolve BDNF cellular sources and real-time signaling dynamics across distinct cortical layers and cell types; (2) genotype-stratified dose-response studies to establish safe, effective, and personalized delivery protocols; and (3) large-scale multicenter clinical trials that seamlessly integrate BDNF-targeted interventions with multimodal rehabilitation strategies, such as repetitive transcranial magnetic stimulation and mirror therapy, with longitudinal neuroimaging serving as a biomarker for cortical reclamation. By systematically addressing these priorities, we aim to transform BDNF from a critical endogenous regulator into a precisely controllable and titratable therapeutic target, ultimately enabling a paradigm shift from mere structural repair to true functional neural reconstruction following PNI.
6.Correlation Between Hyperbaric Oxygen Therapy Course and Reduction of Clonazepam in Patients with Paroxysmal Sympathetic Hyperactivity:A Retrospective Cohort Study
Sujuan LIU ; Yong WANG ; Ran LI ; Haidong LI ; Jie LIU
Medical Journal of Peking Union Medical College Hospital 2025;16(2):386-392
Objective To analyze the correlation between the duration of hyperbaric oxygen(HBO)therapy and oral medication clonazepam reduction in patients with paroxysmal sympathetic hyperactivity(PSH).Methods Clinical data of patients with secondary PSH after severe brain injury at Capital Medical University Affiliated Fuxing Hospital from September 2017 to July 2023 were retrospectively included,covering general information,etiology,lesion location,comorbidities,vital signs at admission,PSH attack characteristics,HBO treatment frequency,PSH treatment drugs and dosage.According to the number of HBO treatments,PSHpatients were divided into HBO short course group(10 treatments)and HBO long course group(>10 treatments).Multiple logistic regression and rank correlation analysis were used to investigate the correla-tion between the duration/frequency of HBO treatment and the reduction of clonazepam.Results A total of 75 PSH patients who met the inclusion and exclusion criteria were selected,including 38 cases(50.7%)who re-ceived reduced doses of clonazepam.There were 32 cases(42.7%)in the HBO short course group and 43 ca-ses(52.3%)in the HBO long course group.The reduction rate of clonazepam in the HBO short course group was lower than that in the HBO long course group[31.3%(10/32)vs.65.1%(28/43),crude OR=4.11(95%CI:1.55-10.90),P=0.004].After adjusting for confounding variables,a multivariate Logistic re-gression model showed a significant correlation between HBO long course and clonazepam reduction(OR=3.76,95%CI:1.23-11.55,P=0.021).Rank correlation analysis showed a positive correlation between HBO treatment frequency and clonazepam reduction(rs=0.331,P=0.004).Conclusions Long course HBO treatment was positively correlated with oral reduction of clonazepam in PSH patients,which may help re-duce the side effects caused by clonazepam and become a non-pharmacological treatment option for PSH.
7.n-butanol fraction of ethanol extract of Periploca forrestii Schltr.: its active components, targets and pathways for treating Alcheimer's disease in rats.
Niandong RAN ; Jie LIU ; Jian XU ; Yongping ZHANG ; Jiangtao GUO
Journal of Southern Medical University 2025;45(4):785-798
OBJECTIVES:
To investigate the active components and possible mechanisms of n-butanol fraction of Periploca forrestii Schltr. ethanol extract for treating Alzheimer's disease (AD).
METHODS:
The active components of n-butanol fraction of Periploca forrestii Schltr. ethanol extract were analyzed using UPLC-QE-MS technique. In a SD rat model of AD induced by treatment with AlCl3 and D-gal, the therapeutic effects of low, moderate and high doses of the n-butanol fraction, saline, and donepezil hydrochloride were evaluated using ELISA, HE and Nissl staining, immunohistochemistry and Western blotting. The therapeutic mechanisms of the n-butanol fraction were explored using network pharmacology and molecular docking.
RESULTS:
Seventeen active components were identified from the n-butanol fraction of Periploca forrestii Schltr. ethanol extract, including phenylpropanoids, flavonoids, anthraquinones, triterpenoids, steroids, and volatile oils. In the rat models of AD, treatment with the n-butanol fraction significantly lowed AChE content in the hippocampus, increased the contents of ACh, SOD, CAT, and GSH-Px, enhanced the expressions of neuronal apoptotic factors Bcl-2, PI3K, Akt, p-PI3K, and p-Akt, and reduced the expressions of Bax and caspase-3 proteins. The treatment also dose-dependently up-regulated hippocampal expressions of Nrf-2, HO-1 and BDNF and down-regulated Keap-1, Aβ and Tau expressions. Bioinformatics analysis identified 14 key intersected targets (including TNF, AKT1 and ESR1) between the n-butanol fraction and AD.
CONCLUSIONS
The therapeutic effect of n-butanol fraction of Periploca forrestii Schltr. ethanol extract in AD mice is mediated by its multiple active components that regulate multiple targets and pathways.
Animals
;
Rats, Sprague-Dawley
;
Rats
;
1-Butanol/chemistry*
;
Plant Extracts/pharmacology*
;
Periploca/chemistry*
;
Ethanol/chemistry*
;
Alzheimer Disease/drug therapy*
;
Male
;
Molecular Docking Simulation
;
Apoptosis/drug effects*
8.Association between Organochlorine Exposures and Lung Functions Modified by Thyroid Hormones and Mediated by Inflammatory Factors among Healthy Older Adults.
Xiao Jie GUO ; Hui Min REN ; Ji Ran ZHANG ; Xiao MA ; Shi Lu TONG ; Song TANG ; Chen MAO ; Xiao Ming SHI
Biomedical and Environmental Sciences 2025;38(2):144-153
OBJECTIVE:
To examine the mechanistic of organochlorine-associated changes in lung function.
METHODS:
This study investigated 76 healthy older adults in Jinan, Shandong Province, over a five-month period. Personal exposure to organochlorines was quantified using wearable passive samplers, while inflammatory factors and thyroid hormones were analyzed from blood samples. Participants' lung function was evaluated. After stratifying participants according to their thyroid hormone levels, we analyzed the differential effects of organochlorine exposure on lung function and inflammatory factors across the low and high thyroid hormone groups. Mediation analysis was further conducted to elucidate the relationships among organochlorine exposures, inflammatory factors, and lung function.
RESULTS:
Bis (2-chloro-1-methylethyl) ether (BCIE), was negatively associated with forced vital capacity (FVC, -2.05%, 95% CI: -3.11% to -0.97%), and associated with changes in inflammatory factors such as interleukin (IL)-2, IL-7, IL-8, and IL-13 in the low thyroid hormone group. The mediation analysis indicated a mediating effect of IL-2 (15.63%, 95% CI: 0.91% to 44.64%) and IL-13 (13.94%, 95% CI: 0.52% to 41.07%) in the association between BCIE exposure and FVC.
CONCLUSION
Lung function and inflammatory factors exhibited an increased sensitivity to organochlorine exposure at lower thyroid hormone levels, with inflammatory factors potentially mediating the adverse effects of organochlorines on lung function.
Environmental Exposure
;
Hydrocarbons, Chlorinated/metabolism*
;
China
;
Ethyl Ethers/metabolism*
;
Environmental Monitoring
;
Thyroid Hormones/blood*
;
Lung/physiology*
;
Inhalation Exposure/statistics & numerical data*
;
Air Pollution/statistics & numerical data*
;
Air Pollutants/metabolism*
;
Humans
;
Male
;
Female
;
Middle Aged
;
Aged
9.Comparative Transcriptomic and Metabolomic Analyses Reveal the Mechanism by Which Foam Macrophages Restrict Survival of Intracellular Mycobacterium Tuberculosis.
Xiao PENG ; Yuan Yuan LIU ; Li Yao CHEN ; Hui YANG ; Yan CHANG ; Ye Ran YANG ; Xuan ZHANG ; An Na JIA ; Yong Bo YU ; Yong Li GUO ; Jie LU
Biomedical and Environmental Sciences 2025;38(7):781-791
OBJECTIVES:
This study aimed to investigate the impact of foam macrophages (FMs) on the intracellular survival of Mycobacterium tuberculosis (MTB) and identify the molecular mechanisms influencing MTB survival.
METHODS:
An in vitro FM model was established using oleic acid induction. Transcriptomic and metabolomic analyses were conducted to identify the key molecular pathways involved in FM-mediated MTB survival.
RESULTS:
Induced FMs effectively restricted MTB survival. Transcriptomic and metabolomic profiling revealed distinct changes in gene and metabolite expression in FMs during MTB infection compared with normal macrophages. Integrated analyses identified significant alterations in the cyclic adenosine monophosphate (cAMP) signaling pathway, indicating that its activation contributes to the FM-mediated restriction of MTB survival.
CONCLUSIONS
FMs inhibit MTB survival. The cAMP signaling pathway is a key contributor. These findings enhance the understanding of the role of FMs in tuberculosis progression, suggest potential targets for host-directed therapies, and offer new directions for developing diagnostic and therapeutic strategies against tuberculosis.
Mycobacterium tuberculosis/physiology*
;
Transcriptome
;
Metabolomics
;
Foam Cells/microbiology*
;
Humans
;
Metabolome
;
Tuberculosis/microbiology*
;
Gene Expression Profiling
10.Differences in growth and secondary metabolite accumulation of Panax quinquefolius between understory and field planting in Shandong, China.
Yue WANG ; Xin-Ying MAO ; Yu DING ; Hong-Xia YU ; Zhi-Fang RAN ; Xiao-Li CHEN ; Jie ZHOU
China Journal of Chinese Materia Medica 2025;50(6):1524-1533
In order to compare the differences in growth and secondary metabolite accumulation of Panax quinquefolius between understory and field planting, growth indexes, photosynthetic characteristics, soil enzyme activities, secondary metabolite contents, and antioxidant activities of P. quinquefolius under different planting modes were examined and compared, and One-way analysis of variance(ANOVA) and correlation analyses were carried out by using the software SPSS 25.0 and GraphPad Prism 9.5. The Origin 2021 software was used for plotting. The results showed that compared with those under field planting, the plant height, leaf length, leaf width, photosynthetic rate, and chlorophyll content of P. quinquefolius under understory planting were significantly reduced, and arbuscular mycorrhizal fungi(AMF) infestation rate and infestation intensity, ginsenoside content, and antioxidant activity were significantly increased. The activities of inter-root soil urease, sucrase, and catalase increased, while the activities of non-inter-root soil urease and alkaline phosphatase increased. Correlation analyses showed that the plant height and leaf length of P. quinquefolius plant were significantly positively correlated with net photosynthetic rate, transpiration rate, chlorophyll content, and electron transfer rate(P<0.05), while ginsenoside content was significantly negatively correlated with net photosynthetic rate, chlorophyll content, and electron transfer rate(P<0.05) and significantly positively correlated with AMF infestation rate and infestation intensity(P<0.05). In addition, ginsenoside content was significantly positively correlated with the activities of inter-root soil sucrase, urease, and catalase(P<0.05). This study provides basic data for revealing the mechanism of secondary metabolite accumulation in P. quinquefolius under understory planting and for exploring and practicing the ecological mode of P. quinquefolius under understory planting.
Panax/microbiology*
;
China
;
Secondary Metabolism
;
Soil/chemistry*
;
Photosynthesis
;
Plant Leaves/metabolism*
;
Chlorophyll/metabolism*
;
Mycorrhizae


Result Analysis
Print
Save
E-mail