Loading...

Table of Content

    16 September 2026, Volume 59 Issue 18
    CROP GENETICS & BREEDING·GERMPLASM RESOURCES·MOLECULAR GENETICS
    QTL Mapping and Candidate Gene Analysis of Flag Leaf Size in Common Wheat (Triticum aestivum L.)
    GUAN PanFeng, LI MengXin, WAN ZiDong, GE Shuai, WANG XinYan, CUI DongJie, CHI Qing, SUN Hao, MA RuoNan, LI Bing, XU HangBo, ZHANG ShiJie, GUO ShuQing, LU JiaLing, JIAO Zhen, ZHU YuPan, WEI WenHui
    Scientia Agricultura Sinica. 2026, 59(18):  3949-3972.  doi:10.3864/j.issn.0578-1752.2026.18.001
    Abstract ( 26 )   HTML ( 4 )   PDF (3622KB) ( 13 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Flag leaf size is a critical component of wheat plant architecture, playing a vital role in photosynthesis, yield formation, and stress response. The identification of major quantitative trait loci (QTL) and candidate genes regulating flag leaf size-related traits will provide theoretical support and genetic resources for the molecular design breeding of high-yielding wheat ideotypes.【Method】A recombinant inbred line (RIL) population comprising 184 lines was derived from a cross between the Chinese elite cultivar Bainong 4199 (BN4199, female parent) and the exotic variety Zhengyinmai 2 (ZYM2, male parent). Genotyping of the BN4199/ZYM2-RIL population was performed using the wheat 120K liquid chip (120K-4HWA) to construct a high-density genetic linkage map. Based on the composite interval mapping (CIM) method, QTL mapping was performed to identify stable major-effect loci and to analyze their additive effects using the flag leaf size-related traits, including flag leaf length (FLL), flag leaf width (FLW), flag leaf length/width ratio (FLR), and flag leaf area (FLA), and the best linear unbiased estimate (BLUE) values obtained across four environments (two locations over two years). Furthermore, candidate gene annotation and identification, as well as tissue expression pattern analysis, were conducted using several wheat functional genomics databases. In addition, corresponding kompetitive allele-specific PCR markers (KASP) were developed based on polymorphic SNP sites located within the stable QTL intervals.【Result】A high-density genetic linkage map of the BN4199/ZYM2-RIL population was constructed, comprising 27 linkage groups with 3 666 SNP markers, covering all 21 wheat chromosomes. The map spanned 5 194.09 cM with an average inter-marker distance of 1.42 cM. High collinearity with the Chinese Spring physical map indicated the robustness of this genetic map. A total of 130 QTL were detected, including: 31 for FLL, 34 for FLW, 33 for FLR, and 32 for FLA, distributed across all chromosomes except 3B, 3D, and 4D. Among these, 20 loci on chromosomes 1A, 2A, 2B, 2D, 4B, 5A, 7A, and 7B were environmentally stable QTL, with 10 being novel loci reported for the first time. A total of 3 923 genes were annotated within the 20 stable QTL intervals. Gene Ontology (GO) analysis revealed significant enrichment in transcription elongation factor complex, UDP-glycosyltransferase activity, rhythmic processes, and chromatin-mediated maintenance of transcription. Based on gene function annotation, information from the cloning and characterization of rice orthologs, and analysis via the QTG-miner module of an interactive wheat functional gene discovery platform, 55 candidate genes were prioritized within stable QTL intervals, exhibiting different expression patterns at various developmental stages of wheat root, stem, leaf, spike and grain tissues. In addition, 10 KASP markers closely linked to the novel QTL were developed.【Conclusion】A high-density genetic linkage map was constructed using the wheat 120K liquid chip, and 20 stable major QTL controlling flag leaf size-related traits were detected, 10 of which were novel loci. The KASP molecular markers closely linked to novel QTL were developed, and 55 potential candidate genes were identified in this study.

    The Effect of the Tuber-Specific Protein StSP6A on in vitro Potato Tuber Formation Under High-Temperature Stress
    DING Ning, LI QingQuan, XU Jing, VLADIMIR N. Obraztsov, JIA XiaoXia, HUANG Wei, QI EnFang
    Scientia Agricultura Sinica. 2026, 59(18):  3973-3988.  doi:10.3864/j.issn.0578-1752.2026.18.002
    Abstract ( 13 )   HTML ( 1 )   PDF (4642KB) ( 3 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】High temperatures, as a major abiotic stress factor, affect potato tuber formation. Exploring the mechanisms underlying potato tuber formation under high-temperature stress holds significant theoretical value for elucidating the intrinsic mechanisms by which environmental factors regulate tuber development.【Method】Using in vitro seedlings from 103 potato germplasm accessions as experimental material, the seedlings were subjected to either optimal temperature (17 ℃) or high temperature (28 ℃) in plant growth chamber for 150 days. After collecting and analyzing the phenotypic traits of the seedlings, representative germplasm accessions were selected for transcriptomic analysis, which was validated using real-time quantitative PCR (RT-qPCR). Furthermore, following the construction of an SP6A overexpression vector, genetic transformation was performed using Longshu 10. After subjecting the transgenic material to high-temperature treatment, phenotypic evaluation and molecular-level analysis were conducted.【Result】After treatment at 17 ℃ and 28 ℃, 103 accessions were classified into four temperature response patterns: widespread tuber-forming type (GroupⅠ, 46 accessions; tubers formed at both temperatures), heat-sensitive type (GroupⅡ, 33 accessions; tubers formed only at 17 ℃), high temperatures /high sugar induced sensitive type (GroupⅢ, 11 accessions; tubers formed only at 28 ℃), and non-tuber-forming (GroupⅣ, 13 accessions; no tubers formed at either temperature). Furthermore, a variance analysis indicated that high-temperature stress shortened the number of days to tuber formation, reduced plant height, and significantly decreased the weight and number of test-tube tubers. At the same time, the coefficients of variation for those indicators all increased significantly. Transcriptome analysis revealed that Longshu 11 (which tuberizes only at 17 ℃) and Lily (which tuberizes only at 28 ℃) had 3 398 and 2 197 differentially expressed genes, respectively. with the differentially expressed genes in Longshu 11 focused on molecular functions such as catalytic activity, as well as metabolic processes and pathways, while those in Lily were enriched in cell wall remodeling and carbohydrate metabolism pathways. A combined transcriptomic analysis revealed that the two groups shared 730 differentially expressed genes, which were classified into eight categories based on their expression patterns. Among the 14 key differentially expressed genes, the tuber-specific gene SP6A is the primary factor responsible for the differences in tuber formation in in vitro tubers under high-temperature stress. An increased number of microtubers was observed in the overexpression plants (SP6A-OE) under 17 ℃ treatment. Under 28 ℃ heat treatment, although the formation of in vitro tubers was inhibited in both the SP6A-OE and control, the SP6A-OE line still produced more in vitro tubers. Western blot results indicated that the expression level of the SP6A protein in SP6A-OE plants was lower at 28 ℃ than at 17 ℃. Meanwhile, RT-qPCR results showed that SP6A expression was suppressed at high temperatures, whereas the expression levels of both SP5G and HSFA6b were significantly elevated at 28 ℃.【Conclusion】Under high-temperature stress, different potato germplasm exhibit varying tuber-forming habits. Overexpression of StSP6A alleviates the inhibitory effect caused by heat stress on axillary tuber formation, thereby enhancing the ability of in-vitro seedlings to form tubers under high-temperature conditions.

    TILLAGE & CULTIVATION·PHYSIOLOGY & BIOCHEMISTRY·AGRICULTURE INFORMATION TECHNOLOGY
    Wheat Adaptability Evaluation and Identification of High-Yield and Stable-Yield Cultivars Based on Time-Series Multi-Temporal Canopy Spectroscopy
    YU LiQiang, YAO DongZe, FEI ShuaiPeng, LI HuiLi, GUO XianRui, YAO DanYu, ZHANG HongJun, ZHOU Yang, XIAO YongGui
    Scientia Agricultura Sinica. 2026, 59(18):  3989-4003.  doi:10.3864/j.issn.0578-1752.2026.18.003
    Abstract ( 14 )   HTML ( 1 )   PDF (3416KB) ( 6 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Aiming at the problem that traditional single-temporal spectroscopy is difficult to accurately acquire time-series dynamic phenotyping throughout the whole wheat growth period and has insufficient capability for evaluating high-yield and stable-yield potential, this study integrates the full-growth-period time-series hyperspectral features with the Stacking ensemble learning algorithm to construct a method for evaluating and selecting wheat germplasm with stable yield potential, which takes into account both qualitative classification and quantitative yield estimation.【Method】This study selected 100 winter wheat genotypes as test materials. Canopy hyperspectral images were collected at 16 key growth stages over two consecutive years (2022-2023 and 2023-2024). The temporal mean (cumulative photosynthetic potential), standard deviation (stay-green stability) and maximum value (maximum growth potential) of the time-series spectral vegetation indices were calculated. A 21-dimensional feature space containing temporal statistical features and interaction features was constructed for dimensionality reduction and information enrichment. After applying within-year standardization to eliminate annual sensor and environmental drift, the K-Means algorithm (K=3) was applied to classify 195 valid variety-year samples under an unsupervised framework. On this basis, a stacking ensemble regression model was trained on 92 samples from the first year to predict the yield of 98 unseen genotypes in the second year (Leave-One-Year-Out validation), with Variance Inflation applied to address the regression toward the mean problem.【Result】The principal components analysis (PCA) results showed that the first two PCs successfully interpreted cumulative photosynthetic potential and growth stability, with a cumulative variance contribution rate exceeding 85%. The subsequent K-means clustering classified the varieties into three categories: high-yield and stable-yield, general adaptability and low-yield sensitive. Subgroup-level yield variance analysis showed highly significant differences in measured plot yield among the three groups (F=17.83, P=7.91×10-8), and the average plot yield of the high-yield stable-yield group was 17.61% higher than that of the low-yield sensitive group. Evaluated against within-year yield tertiles as ground truth, the overall classification accuracy (OA) reached 50.77% (Kappa=0.263), and the safety retention rate for high-yielding genotypes reached 93.94%, confirming that this unsupervised framework could effectively identify and eliminate low-performing germplasm. Feature analysis showed that whole-growth-period GNDVI and OSAVI contributed highly to evaluating the high-yielding groups. This physiologically explained that green and red-edge bands had stronger canopy penetration than the red band, effectively overcoming the saturation limitation of traditional NDVI under high canopy closure in the mid-to-late growth stages. In the cross-year independent validation, the predicted yields were significantly correlated with the measured yields (r=0.51, P=1.99×10-7, RMSE=0.38 t·hm-2), and the prediction slope was successfully adjusted to 0.31.【Conclusion】This study verified the effectiveness of time-series dynamic canopy spectroscopy in identifying stable and high-yielding wheat genotypes, and proposed a breeding decision-making framework combining qualitative screening with quantitative ranking, which provided a methodological reference for germplasm selection.

    SPECIAL FOCUS: MECHANISM ANALYSIS, TECHNOLOGY OPTIMIZATION AND DECISION-MAKING SUPPORT FOR PRECISION PREVENTION AND CONTROL OF AGRICULTURAL PESTS
    Ligand Binding Property and Key Amino Acid Analysis of Odorant- Binding Protein OBP11 in Bradysia odoriphaga
    WANG QinLi, TAN ShunYang, XIA Gen, LIU WenQing, YANG YuTing, TIAN LiXia
    Scientia Agricultura Sinica. 2026, 59(18):  4007-4020.  doi:10.3864/j.issn.0578-1752.2026.18.005
    Abstract ( 17 )   HTML ( 4 )   PDF (5427KB) ( 10 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】The objective of this study is to elucidate the molecular characteristics of the odorant-binding protein BodoOBP11, which is highly expressed in the antennae of adult Bradysia odoriphaga, and to investigate its olfactory function in host localization.【Method】Bioinformatics tools were used to analyze the basic characteristics and evolutionary relationships of BodoOBP11. The spatio-temporal expression patterns of BodoOBP11 across different developmental stages and tissues were then assessed using quantitative real-time PCR (qPCR). The BodoOBP11 protein was obtained through prokaryotic expression in E. coli, purification via nickel affinity chromatography, and cleavage using recombinant enterokinase. Using fluorescence competitive binding assays with the fluorescent probe 1-NPN, the binding ability of the protein to 22 volatile compounds from the host plant (Chinese chives) and one putative sex pheromone compound was subsequently evaluated. Finally, protein structure prediction, molecular docking, and site-directed mutagenesis were employed to identify the key binding sites between BodoOBP11 and ligand compounds.【Result】BodoOBP11 has an open reading frame of 426 bp, encoding 141 amino acid residues, including an N-terminal signal peptide of 20 amino acids. It exhibits the typical features of the Classic-C OBP, characterized by six conserved cysteine residues. The predicted molecular weight of the mature protein is 13.53 kDa, with an isoelectric point of 7.71. qPCR analysis revealed distinct expression patterns of BodoOBP11. Across developmental stages, its expression was highest in male adults, followed by female adults and pupae. In different tissues, BodoOBP11 showed significant high expression in the antennae of both sexes, with the highest level detected in male antennae. Fluorescence competitive binding assays demonstrated that BodoOBP11 binds to six volatiles: dipropyl trisulfide, 2-tridecanone, dipropyl disulfide, 2-undecanone, phenethyl isothiocyanate, and α-ionone. The strongest binding affinity was observed for dipropyl trisulfide. Protein structure prediction and molecular docking indicated that BodoOBP11 binds to these six ligands primarily via hydrogen bonding and hydrophobic interactions. Specifically, Thr74 forms hydrogen bonds with 2-tridecanone, 2-undecanone, and α-ionone, while residues Val76, Ile107, Leu127, Met131, Phe138, and Phe140 mainly mediate hydrophobic interactions. Site-directed mutagenesis followed by fluorescence competitive binding assays revealed that mutations of Ile107, Phe138, and Phe140 to alanine significantly reduced the binding affinity of BodoOBP11 to all six ligands, indicating their critical roles in binding.【Conclusion】BodoOBP11 exhibits strong binding affinity to six volatile compounds from Chinese chive, including dipropyl trisulfide, 2-tridecanone, and dipropyl disulfide, etc. Furthermore, residues Ile107, Phe138, and Phe140 were identified as crucial sites for binding these ligands. These findings preliminarily reveal its molecular mechanism in mediating host plant volatile recognition by adult B. odoriphaga using BodoOBP11, offering a theoretical foundation for developing green control technologies based on olfactory interference.

    Effects of Refrigeration and Embryo-Killing Treatment of Corcyra cephalonica Eggs on Propagation Efficiency of Egg Parasitoid Trichogramma dendrolimi
    SUN JiaNing, LIU ChenXi, LIU XiaoXia, LI Zhen
    Scientia Agricultura Sinica. 2026, 59(18):  4021-4032.  doi:10.3864/j.issn.0578-1752.2026.18.006
    Abstract ( 13 )   HTML ( 1 )   PDF (627KB) ( 2 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Background】Efficient artificial propagation of egg parasitoids requires high quality alternative host eggs, so the fresh and nutritional states of host eggs are key factors to the propagation efficiency of egg parasitoids.【Objective】The Trichogramma dendrolimi, an egg parasitoid widely used in agricultural production, and its dominant alternative host Corcyra cephalonica eggs were used as test materials to explore the effects of refrigeration and embryo-killing on the quality of host eggs, so as to provide data support for scientific preservation of host eggs, optimization of artificial propagation technology of egg parasitoids, improvement of propagation efficiency and reduction of control costs.【Method】Four treatments of 1-day-old C. cephalonica eggs were performed: 4 ℃ refrigeration for 3 d followed by UV irradiation for embryo-killing (3 d+K), UV irradiation followed by 4 ℃ storage for 3 d (K+3 d), 4 ℃ storage for 7 d followed by UV irradiation (7 d+K), and UV irradiation followed by 4 ℃ storage for 7 d (K+7 d), with eggs only treated with UV irradiation for embryo-killing (K) as control. The parasitism rate of T. dendrolimi and the eclosion rate, female ratio, and developmental duration of the offspring parasitoids in eggs resulted from different treatments were compared. The contents of triglyceride, total protein, and trehalose in C. cephalonica eggs from each treatment were quantified, as well as the Wolbachia levels in both C. cephalonica eggs and T. dendrolimi emerged from the eggs, and their correlations with the biological indices of T. dendrolimi were analyzed.【Result】Compared with the control (K), the 7 d+K treatment significantly increased about 20% of the parasitism rate, and significantly increased about 72% of the offspring female amount of T. dendrolimi, without significant impact on eclosion rate and developmental duration of the offspring. The K+7 d treatment also could significantly increase 19% parasitism rate of the T. dendrolimi, but decreased the eclosion rate of the offspring wasps. Nutrient analysis showed that 4 ℃ storage (either before or after embryo-killing) significantly increased the triglyceride content in C. cephalonica eggs. Refrigeration treatment for 3 d before UV irradiation (3 d+K) significantly increased the protein content. None of the treatments significantly affected the trehalose content in host eggs. Correlation analysis revealed that triglyceride content in the eggs was significantly positively correlated with the parasitism rate of T. dendrolimi, the female proportion of its offspring and the Wolbachia level in the offspring wasps.【Conclusion】Pretreatment of 1-day-old C. cephalonica eggs with 7 d of refrigeration at 4 ℃ prior to embryo-killing under UV can effectively extend the duration of their freshness. This pretreatment may promote the parasitism efficiency and offspring female ratio of T. dendrolimi by increasing nutrient content, such as triglyceride, in host eggs. This study provides an optimized preservation and processing protocol for efficient and low-cost artificial mass rearing of T. dendrolimi using C. cephalonica eggs as an alternative host.

    Monitoring of Resistance to Omethoate in Field Populations of Aphis gossypii in China and Elucidation of Target Resistance Mechanisms
    LI Ren, CHENG ShenHang, LIANG PingZhuo, ZHANG Lei, WANG ShanNing
    Scientia Agricultura Sinica. 2026, 59(18):  4033-4047.  doi:10.3864/j.issn.0578-1752.2026.18.007
    Abstract ( 16 )   HTML ( 1 )   PDF (6121KB) ( 14 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】The objective of this study is to determine the resistance level of Aphis gossypii populations from the main cotton planting areas of China to omethoate, to elucidate the predominant mutant types in acetylcholinesterase 1 (AChE1) (the target of omethoate) in these populations, and further analyze the impacts on the AChE1 structure by these major mutation types and their relationship with omethoate resistance.【Method】In 2019-2023, A. gossypii field populations were collected from the main cotton planting areas of China, and the resistance levels of these populations to omethoate were determined using the leaf-dipping method. Mutations in AChE1 were detected by single-aphid PCR amplification and sequencing, and the types of mutations carried by an individual aphid were analyzed. Monoclonal strains carrying different AChE1 mutations were established using single-female rearing protocols, and their contributions to omethoate resistance were subsequently assessed through toxicity bioassays and synergist experiments. The mechanisms of A. gossypii resistance to omethoate mediated by these mutations were revealed using homology modeling, molecular docking, and molecular dynamics simulations.【Result】In 2019-2023, 38 A. gossypii field populations in China developed moderate to high levels of resistance to omethoate (17.1-207.8-fold) and remained relatively stable over the years. G221A, A302S, A332V, and S431F mutations in AChE1 were prevalent in these populations, with S431F mutation frequency reaching 100%, predominantly in homozygous form. While the average frequencies of G221A, A302S, and A332V were 37.26%, 51.26%, and 74.64%, respectively. Detection of mutant types from a single aphid revealed that individuals carrying the A332V+S431F combination mutation had the highest frequency (35.24%), followed by G221A+A302S+A332V+S431F (26.80%) and A302S+A332V+S431F (12.12%). The average frequency of the S431F single mutation was 9.56%. The weighted frequencies of these major mutation types showed a certain linear correlation with omethoate resistance (R2=0.46, P<0.001). By constructing monoclonal mutant strains, it was confirmed that A. gossypii strains carrying the combination mutations A332V+S431F, G221A+A302S+S431F, A302S+A332V+S431F, and G221A+A302S+A332V+S431F developed a moderate level of resistance to omethoate, with resistance ratios of 46.1-77.8-fold, significantly higher than that of the S431F single-mutant strain (12.5-fold). Moreover, PBO, TPP, and DEM did not exhibit significant synergistic effects on omethoate toxicity in these mutant strains. Molecular docking and dynamics simulations confirmed that the S431F single mutation and A332V+S431F, G221A+A302S+S431F, A302S+A332V+S431F, and G221A+A302S+A332V+S431F combination mutations significantly increased the geometric center distance between the AChE1 catalytic triad and omethoate by 0.074-0.185 nm compared to the wild-type. These mutations also altered the hydrophobic area and radius of gyration of the AChE1 active center, thereby affecting the active center of A. gossypii AChE1 and its binding conformation with omethoate.【Conclusion】A. gossypii populations in major cotton-planting regions of China have developed moderate to high levels of resistance to omethoate. Four mutations were identified in AChE1 across these populations, and these site mutations altered the active center of AChE1 and its binding affinity for omethoate through different combination types, thereby mediating field resistance to this insecticide. In field control practices, it is advisable to reduce the use of AChE-targeting insecticides to prevent the development of cross-resistance mediated by target mutation and further escalation of resistance levels.

    Resistance Monitoring and Inheritance of Resistance to Spinetoram in Spodoptera litura
    LI Hao, GAO XinJu, LU Peng, ZHAN QianYuan, WANG Ran, SU Qi, QU Cheng
    Scientia Agricultura Sinica. 2026, 59(18):  4048-4055.  doi:10.3864/j.issn.0578-1752.2026.18.008
    Abstract ( 13 )   HTML ( 1 )   PDF (582KB) ( 9 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】This study aims to clarify the current resistance status and underlying mechanisms of field populations of Spodoptera litura to spinetoram, and to provide a scientific basis for field resistance monitoring and the development of integrated resistance management (IRM) strategies.【Method】From 2023 to 2024, 23 field populations of S. litura were collected across China. The resistance levels to spinetoram were determined using the leaf-dipping method. The highly resistant Nanning (NN) population and the susceptible strain (SS) were subjected to reciprocal crosses and mass crosses, and the F1 progeny larvae were bioassayed to analyze the resistance inheritance pattern. The synergistic effects of piperonyl butoxide (PBO), diethyl maleate (DEM), and triphenyl phosphate (TPP) were evaluated, and the activities of carboxylesterase (CarE), glutathione S-transferase (GST), and cytochrome P450 monooxygenase (P450) were measured in both the NN and SS strains.【Result】Field resistance monitoring revealed significant variation in spinetoram resistance among different geographic populations of S. litura. High-level resistance (43.3- to 107.8-fold) was detected in populations from Nanning (Guangxi) and Fuzhou (Fujian) in 2023, and from Sanya (Hainan), Kunming (Yunnan), Weifang (Shandong), and Xuchang (Henan) in 2024. Notably, the resistance ratios for the Nanning (2023) and Sanya (2024) populations reached 106.8- and 107.8-fold, respectively. Moderate resistance levels were observed in the 2023 Ningbo (Zhejiang) and Jingzhou (Hubei) populations, as well as the 2024 Fuzhou (Fujian), Jingzhou (Hubei), Guangzhou (Guangdong), Changsha (Hunan), Nanjing (Jiangsu), Shanghai, and Guilin (Guangxi) populations, with resistance ratios ranging from 10.3- to 31.8-fold. Low resistance levels were found in the 2023 Hangzhou (Zhejiang) population and the 2024 Nanchang (Jiangxi), Hefei (Anhui), Beijing, and Hangzhou (Zhejiang) populations, with resistance ratios between 5.0- and 8.9-fold. Inheritance analysis showed that resistance to spinetoram in the NN population was controlled by an autosomal, incompletely recessive gene. The synergists DEM and TPP exhibited synergistic ratios of 1.8 and 1.5, respectively, suppressing resistance in the NN population. Moreover, the activities of GST and CarE in the NN population were 4.1- and 3.5-fold higher than those in the susceptible strain, respectively.【Conclusion】Resistance to spinetoram has developed in S. litura populations across multiple regions. Enhanced activities of GST and CarE play important roles in the resistance of S. litura to spinetoram.

    The Importance of Different Crop Types to the Hawkmoth-Host Plant Interaction Network in Beijing Area
    LI HaoYuan, YANG CaiQing, CHEN Dan, XU ZhenBang, MEN YaHui, GAO Fei, ZHANG AiBing
    Scientia Agricultura Sinica. 2026, 59(18):  4056-4075.  doi:10.3864/j.issn.0578-1752.2026.18.009
    Abstract ( 9 )   HTML ( 2 )   PDF (4512KB) ( 3 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Hawkmoth larvae feed on the leaves of host plants, causing plant damage or even death, thereby posing a serious threat to regional agricultural and forestry development. This paper aims to reveal the potential importance of different types of economic crops to hawkmoth pests in Beijing area, identify the main hawkmoth pest species associated with each crop type, and propose prevention and control strategies.【Method】A hawkmoth-host plant interaction network in Beijing area was constructed. Network structure was analyzed using metrics such as connectance, nestedness, and modularity. The roles of individual nodes within the network’s modular structure were quantified using inter-module connectivity (c-value) and within-module degree (z-value). The impact of nodes on network robustness was assessed using secondary extinction curves. The Kruskal-Wallis test was employed to determine whether there were significant differences in the effects of different crop types on network robustness and modular structure. For hawkmoth species with high connectivity in the network, their dietary composition was analyzed, and prevention and control recommendations were proposed based on different crop types.【Result】The hawkmoth-host plant interaction network exhibited a connectance C=0.043, nestedness N=3.67, and modularity Q=0.668. The agricultural crop category contained two module hubs, the forestry crop category contained two module hubs, the medicinal crop category contained two module hubs and three connectors, the ornamental crop category contained two module hubs and two connectors, the native tree category contained three module hubs and one connector, the other crop category contained four module hubs and two connectors, and plants in the non-economic crop category were all peripheral nodes. There was no significant difference among crop categories in their impact on network robustness. However, the within-module degree of all economic crop categories was significantly higher than that of non-economic crops (agricultural crops: P=0.0015; forestry crops: P=0.0030; medicinal crops: P=0.0010; ornamental crops: P=0.0069; native trees: P=0.0006; other crops: P=0.0011).【Conclusion】The hawkmoth-host plant interaction network in Beijing area exhibits a sparse and dispersed interaction pattern. Different hawkmoth pests pose high potential risks to various economic crop types. For prevention and control efforts in Beijing area, species with high centrality in the network, such as Psilogramma increta and Agrius convolvuli, should be prioritized for control across all host plants. For all crop categories, in addition to the protection of key agricultural and forestry crops, other economic crops and native trees should also be included in the scope of protection. For individual crop types, specific control measures should target hawkmoth species that primarily feed on them. Furthermore, for hawkmoth species that have distribution records in Beijing area and potential host plants but currently have small population sizes, vigilance should be exercised regarding the risk of population outbreaks.

    Application of Geostatistics in Insect Biological Control and Its Support for Precision Pest Management Decisions
    LI Shu, MA Chi, WANG JiaJia, XIAO Da, MI YingYing, YUE YanLi, WANG Su
    Scientia Agricultura Sinica. 2026, 59(18):  4076-4091.  doi:10.3864/j.issn.0578-1752.2026.18.010
    Abstract ( 16 )   HTML ( 1 )   PDF (2821KB) ( 9 )   Save
    Figures and Tables | References | Related Articles | Metrics

    The spatial distribution of natural enemy populations as well as their spatio-temporal relationships with pests directly determines the performance of agricultural ecosystems’ natural pest suppression function, and acts as a critical foundation for research on green and precision pest management. Traditional research continues to follow the analytical paradigm of “population monitoring - average assessment - uniform intervention”, focusing on the temporal dynamics of pest density. This confines many studies to merely descriptive analysis of spatial patterns, leaving them unable to underpin practices of pesticide reduction & efficiency improvement and green pest control in terms of spatial heterogeneity identification, pest-natural enemy spatial relationship dissection, and management decision support. By integrating spatial location data and population observation values, geostatistics can effectively characterize the spatial structure of pest and natural enemy populations, identify pest high-density aggregation hotspots and natural enemy dominant distribution zones, and further quantitatively diagnose their spatial overlap zones, tracking relationships and spatial mismatch zones. It thus provides an innovative analytical approach to break through the limitations of traditional methods. We provide a systematic review of research progress and applications of geostatistics in biological pest control, focusing on how it has advanced from describing average pest populations to recognizing spatial heterogeneity, from describing the distribution of single populations to diagnosing natural enemy-pest relationships, and further from analyzing spatial patterns to supporting precision pest control decision-making. We summarize differentiated geostatistical application modes marked by “rapid response” and “stable persistence” in annual field crops and perennial horticultural systems. A pest control decision-making pathway of “data capture - spatial diagnosis - management decision - feedback optimization” is proposed, which systematically converts spatial evidence including pest hotspots, natural enemy coverage and spatial mismatches into implementable management measures such as zoned management, targeted natural enemy release and control efficacy evaluation. In the future, multi-source spatial data integration, cross-scale dynamic modeling, ecological mechanism embedding and intelligent decision-making platform construction should be further strengthened to elevate the interpretive capacity, decision-making ability and application potential of geostatistics in biological pest control.

    SOIL & FERTILIZER·WATER-SAVING IRRIGATION·AGROECOLOGY & ENVIRONMENT
    Dynamics of Black Soil Labile Carbon Fractions and Carbon- Hydrolyzing Enzymes Under Long-Term Fertilization Regimes
    HAN XiaoYu, ZHU YuPing, LIU HongFang, SUN ZhiMei, ZHANG XiuZhi, GAO HongJun, PENG Chang, ZHANG ShuXiang
    Scientia Agricultura Sinica. 2026, 59(18):  4092-4104.  doi:10.3864/j.issn.0578-1752.2026.18.011
    Abstract ( 13 )   HTML ( 3 )   PDF (2170KB) ( 3 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】This study investigated the dynamics and interrelationships between soil labile carbon fractions and carbon-hydrolyzing enzymes under long-term fertilization practices, aiming to provide a theoretical basis for enhancing black soil fertility and optimizing fertilization strategies in Northeast China.【Method】A long-term fertilization experiment established in 1989 in Gongzhuling, Jilin Province, was utilized. Soil samples were collected from five treatments: no fertilization (CK), nitrogen and potassium fertilization (NK), nitrogen-phosphorus-potassium fertilization (NPK), NPK plus straw return (NPKS), and NPK plus manure (NPKM). Measurements included soil nutrients, activities of amylase (AMY), sucrase (SUC), cellobiohydrolase (CBH), β-1,4-glucosidase (βG), and content of particulate organic carbon (POC), easily oxidizable organic carbon (EOC), dissolved organic carbon (DOC), and microbial biomass carbon (MBC). Data were analyzed using ANOVA, correlation analysis, Mantel tests, and partial least squares path modeling to elucidate changes and underlying linkages.【Result】Compared with CK, the NPKM treatment demonstrated the most comprehensive improvement in soil organic carbon and its active components. Soil organic carbon (SOC) content increased by 6.7 g·kg-1, while POC, EOC, MBC, and DOC content was significantly increased by 2.5 g·kg-1, 3.0 g·kg-1, 97.5 mg·kg-1, and 14.1 mg·kg-1, respectively, with the geometric mean of active carbon also rising substantially by 50.1%. Under the NPKS treatment, SOC, POC, and MBC content was also increased significantly by 2.3 g·kg-1, 1.5 g·kg-1, and 93.0 mg·kg-1, respectively. In contrast, both the NPK and NK treatments reduced POC and MBC content, with the geometric mean of active carbon decreasing by 24.7% and 26.8%, respectively. The carbon pool management index increased significantly only under the NPKM treatment by 55.9% compared with CK, while the carbon pool activity decreased significantly by 33.6% under the NK treatment. Regarding enzyme activity, all fertilization treatments enhanced AMY activity compared with CK. Additionally, the NK and NPK treatments significantly increased CBH and βG activities, with the NK treatment showing the greatest increase: CBH rose by 167.3 μmol·g-1·h-1 and βG by 457.7 μmol·g-1·h-1. In terms of soil nutrients, the NPKM treatment significantly increased total nitrogen (TN), total phosphorus (TP), available phosphorus (AP), and available potassium (AK) content, with increases of 0.72 g·kg-1, 1.85 g·kg-1, 148.9 mg·kg-1, and 402.1 mg·kg-1, respectively, compared with CK. Conversely, the NK treatment resulted in a significant decrease in pH by 2 units and a reduction in TP content by 0.2 g·kg-1. Furthermore, the NK treatment increased soil C:P and N:P, while other fertilization treatments generally decreased C:N, C:P, and N:P. Correlation and path analysis indicated that soil TN, TP, AP, and AK directly promoted the accumulation of POC, EOC, and MBC. In contrast, higher C:P and N:P stimulated CBH and βG activities, thereby accelerating carbon decomposition and mineralization, which inhibited the accumulation of POC, EOC, and MBC.【Conclusion】Long-term combined application of chemical fertilizers with organic amendments not only improved soil nutrient content but also increased labile organic carbon by optimizing nutrient balance and modulating carbon-hydrolyzing enzyme activities. Thus, integrating organic manure with fertilizers represents the optimal strategy for enhancing black soil fertility in Northeast China.

    Effects of Grass Cover on Slope Runoff Process in Loess Areas of Northern Shaanxi, China
    WANG Bo, CHEN Hao, QUAN BeiPing, CHEN MingChuan, JING XinRui, XU JiaHao, HUANG ChaoYi, LI ZhanBin
    Scientia Agricultura Sinica. 2026, 59(18):  4105-4120.  doi:10.3864/j.issn.0578-1752.2026.18.012
    Abstract ( 17 )   HTML ( 3 )   PDF (3060KB) ( 6 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Grass cover plays a critical ecological role in regulating slope runoff during rainfall. However, quantitative studies on how root systems and soil properties change during grass cover restoration and subsequently affect slope runoff remain limited.【Method】In this study, indoor simulated rainfall experiments were conducted to systematically examine the effects of varying rainfall intensities (40, 80, and 120 mm·h-1) and grass cover levels (0%, 33%, 48%, 56%, and 69%) on surface runoff characteristics, soil physicochemical properties, root parameters, and their combined influence on slope runoff in ryegrass (Lolium perenne) plots (2 m×0.5 m×0.4 m, 15° slope).【Result】(1) Compared with bare slopes, grass cover restoration significantly improved the physicochemical properties of the surface soil (0-10 cm). Specifically, it led to a 2.4% decrease in soil bulk density and increases of 1.4%, 9.4%, 227.1% and 82.5% in total soil porosity, mean weight diameter of aggregates, saturated hydraulic conductivity, and soil organic matter content, respectively, with all improvements enhanced by increasing cover. Meanwhile, root parameters (root length density, root surface area density, and root volume density) of grass plots increased significantly with increasing cover. Regression analysis further indicated that grass root growth directly drove corresponding changes in topsoil properties, and the improvement in soil saturated hydraulic conductivity was the most significant. (2) Under the same rainfall intensity, the initial runoff time increased significantly with greater grass cover, while the average runoff rate decreased, with runoff reduction benefits increasing accordingly (7.7%, 13.3%, 16.9%, and 27.1%). With increasing rainfall intensity, the initial runoff time was markedly shortened, and the greater the grass coverage, the larger the reduction magnitude. Meanwhile, the average runoff rate rose significantly, accompanied by a decline in corresponding runoff reduction benefits (24.7%, 15.9%, and 8.1%); moreover, the weakening trend of the regulatory effect of grass cover on slope runoff with increasing rainfall intensity was particularly obvious at higher coverages (56% and 69%). (3) The changes in plant roots and soil physicochemical properties during grass restoration significantly influenced slope runoff processes. Path model analysis indicated that the direct effect of plant roots on runoff rate (standardized path coefficient: -0.60) was greater than the indirect effect (standardized path coefficient: -0.27). This was attributed to the fact that the increase in grass cover after restoration mainly regulates slope runoff through the physical barrier effect (direct effect) formed by root density, while its indirect effect was primarily achieved by improving topsoil hydro-physical properties and ultimately increasing soil saturated hydraulic conductivity, thereby enhancing water infiltration capacity.【Conclusion】This study revealed the regulatory effect of grass restoration on slope runoff processed by driving changes in topsoil properties through root growth, and this effect was the most significant under high grass coverage. These results provide a scientific basis for erosion control, grassland restoration, and management in the Loess Plateau region.

    HORTICULTURE
    Identification and Ecological Impact Factor Analysis of Nutritional Component Pinoresinol in Apple Fruits
    YUE JingWen, GUO XinXin, ZHOU Jia, SHI TianLe, YANG YuZhang, ZHANG JunKe, LI XingLiang
    Scientia Agricultura Sinica. 2026, 59(18):  4121-4134.  doi:10.3864/j.issn.0578-1752.2026.18.013
    Abstract ( 10 )   HTML ( 2 )   PDF (3227KB) ( 5 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】The formation of apple fruit quality is closely associated with ecological factors. This study aimed to clarify the relationship between key quality attributes and ecological factors, so as to provide a reference for the specialty cultivation of functional fruits.【Method】Apple fruits of the main cultivars Gala and Fuji were collected from four representative regions, namely Ili and Kashgar in Xinjiang, Luochuan in Shaanxi, and Yantai in Shandong. Widely targeted metabolomics was employed to detect metabolites in apple fruits. Metabolites with the highest identification confidence from each sample were selected, and multivariate statistical analysis was used to screen for differential metabolites that were significantly up-regulated in both Ili and Kashgar compared with Luochuan and Yantai. Overlay analysis was then performed to obtain the common high-content nutritional components. A liquid chromatography-tandem mass spectrometry (LC-MS/MS) method was further established for the identified common high-content components to quantitatively validate the metabolomics results. Simulation treatments of ecological factors, including drought, diurnal temperature difference, and UV-B, were conducted on apple fruits to analyze the accumulation of target components and the expression responses of key genes involved in their biosynthesis, thereby identifying the predominant ecological factor promoting the accumulation of target components.【Result】Based on metabolomics analysis, 1 653 and 2 052 metabolites were detected in Gala and Fuji apple fruits from four regions, respectively. A total of 18 and 19 metabolites were significantly more abundant in Gala and Fuji fruits from the Ili and Kashgar regions than in samples from the Luochuan and Yantai regions, respectively. Among these, the lignan metabolites pinoresinol and epipinoresinol were significantly more abundant in both apple cultivars from the Ili and Kashgar regions than that from the other two regions, which were identified as the common high-content nutritional components. An LC-MS/MS method for pinoresinol determination was further established. Quantitative analysis showed that pinoresinol content in Gala apples from Ili and Kashgar reached 110.2 and 80.4 mg·kg-1, respectively, and those in Fuji apples reached 73.4 and 64.5 mg·kg-1, respectively, all of which were more than twice the content in apples from Yantai and Luochuan. Furthermore, pinoresinol content gradually increased during fruit development and reached the highest level at the mature stage. Drought, large diurnal temperature difference, and UV-B promoted pinoresinol accumulation, whereas high salt stress had no significant promoting effect. A total of 35 dirigent protein (DIR) genes involved in pinoresinol biosynthesis were identified in the apple genome. Among them, MdDIR-a1 and MdDIR-a2 were expressed in apple fruits, and both were significantly up-regulated by drought stress.【Conclusion】Drought, large diurnal temperature difference, and strong UV-B constituted an ecological environment favorable for the accumulation of the nutritional component pinoresinol in apple fruits, among which drought stress was the predominant ecological factor promoting its accumulation.

    Comprehensive Evaluation of Jujube Fruit Quality Based on Sensory, Texture, and Physicochemical Analysis
    SU WanLong, REN HaiYan, XUE XiaoFang, SHI MeiJuan, LIU Li, ZHAO AiLing, LI Yi, WANG YongKang
    Scientia Agricultura Sinica. 2026, 59(18):  4135-4148.  doi:10.3864/j.issn.0578-1752.2026.18.014
    Abstract ( 11 )   HTML ( 4 )   PDF (535KB) ( 4 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】This study aimed to establish a standardized and objective evaluation method for the textural quality of jujube fruits through comprehensive analysis of sensory, textural, and physicochemical traits, so as to provide a solid theoretical basis for fruit quality evaluation and selective breeding of jujube germplasm resources.【Method】A total of 123 jujube germplasms with different utilization types were employed as experimental materials. Correlation analysis, principal component analysis (PCA), significance difference analysis and stepwise linear regression were performed on sensory, textural and physicochemical indices of jujube fruits. The core indicators affecting fruit quality were screened, and a digital comprehensive evaluation model for jujube fruit quality was subsequently constructed.【Result】Descriptive statistical analysis revealed extensive phenotypic variations in sensory, textural and physicochemical traits among different jujube germplasms. The coefficient of variation (CV) of sensory indices ranged from 22.13% to 47.80%, that of textural indices from 13.94% to 30.84%, and that of physicochemical indices from 3.95% to 46.68%; pulp moisture content possessed the minimum CV of 3.95%. Correlation analysis demonstrated that the comprehensive taste quality of jujube fruits was strongly correlated with cell wall component contents. PCA identified five core factors affecting fruit quality from 14 textural and physiological indices, including cell wall components, total phenols and flavonoids, total soluble sugars, pulp moisture content, and soluble protein. The respective contribution rates of the first four principal factors were 28.97%, 13.07%, 10.72%, 8.75%, 7.52%, respectively, with a cumulative contribution rate of 69.03% to the total phenotypic variation. The comprehensive score H of tested germplasms varied from -1.48 to 1.52. The top 5% of germplasms ranked by H value had scores varying from 1.18 to 1.52, while the bottom 5% germplasms ranged from -1.48 to -0.97. Germplasms with high H values exhibited dense fruit texture, low juice yield and thick peel in sensory performance; those with low H values showed crisp texture, high juice yield and thin peel. A significant correlation was detected between total sensory scores and comprehensive H values. Significance difference analysis indicated fresh-eating jujube germplasms exhibited crisp texture, thin-peeled and juicy; drying-type germplasms exhibited denser texture, thick-peeled and low-juice; dual-purpose (fresh and dried) cultivars had moderately dense flesh, medium-thick peel and abundant juice. Finally, stepwise linear regression established a predictive model for sensory score: Y=7.16×10-16+0.26CSP+0.285HCE+0.207WSP+0.183ISP+0.229CE+ 0.132TPC+0.108SC+0.083SP, where Y represented the comprehensive sensory score. The model presented high precision with a determination coefficient of R2 = 0.98.【Conclusion】In conclusion, the model constructed from textural and physicochemical parameters achieved favorable evaluation performance on jujube fruit texture quality and eliminates the subjective bias inherent in conventional sensory evaluation.

    Rapid Prediction of Tea Leaf Catechins via Hyperspectral Multi-Task Synergistic Sensing Fusing Plant Physiological Priors
    CHAI ZhiChao, QIAN XiaoNa, YU YongChen, LI XiWang, SUN XiaoLing
    Scientia Agricultura Sinica. 2026, 59(18):  4149-4162.  doi:10.3864/j.issn.0578-1752.2026.18.015
    Abstract ( 12 )   HTML ( 4 )   PDF (5382KB) ( 2 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】To meet the demand for in situ, real-time, and non-destructive monitoring of tea leaf quality in smart tea gardens, this study proposed a multi-component quantitative prediction method based on visible-near infrared hyperspectral imaging (400-1 000 nm) integrated with plant physiological priors, so as to realize rapid assessment of secondary metabolites, such as catechins.【Method】Twelve tea cultivars were used as experiment materials. Leaves were classed as young (2nd-3rd), mature (4th-5th), and old (6th-7th) to represent physiological gradients, giving 429 independent samples in total. For each leaf, hyperspectral reflectance was recorded, and the contents of catechin (C), epicatechin (EC), epigallocatechin (EGC), epicatechin gallate (ECG), epigallocatechin gallate (EGCG), gallocatechin (GC), gallocatechin gallate (GCG), caffeine (CAF), and gallic acid (GA) were quantified by high-performance liquid chromatography. Spectra were preprocessed using Savitzky-Golay (SG) smoothing, standard normal variate (SNV), and first derivative. The generalization performance of partial least squares regression (PLSR), support vector regression (SVR), XGBoost, multilayer perceptron (MLP), and 1D-CNN-SE was compared. Based on the best single-task model, a multi-task XGBoost framework was then developed to simultaneously predict the nine components. Additionally, the successive projections algorithm (SPA) was used to select key wavelengths from the full spectrum (557 bands) to assess the feasibility of simplified modeling.【Result】When leaf maturity was ignored, EGCG and GCG showed a spurious negative correlation due to confounding by maturity, which was a clear instance of Simpson's paradox. This meant that maturity needed to be accounted for, either by stratifying the data or by explicit adjustment, so that models would not learn false correlations. Among the tested models, XGBoost combined with SG+SNV+first-derivative preprocessing achieved the best performance on the small-sample with high-dimensional spectral dataset. The R2 values for EGCG and GCG were approximately 0.60 and 0.59, respectively. In contrast, 1D-CNN-SE achieved calibration set R2 value above 0.90 but dropped to about 0.40-0.50 on the prediction set, showing clear overfitting. The multi-task XGBoost model predicted all nine components in a single run; the R2 value for EGCG remained about 0.57, and the RPD values for EGCG, GCG, and CAF were about 1.48-1.52, with only a slight loss in single-component accuracy but a substantial gain in throughput. After SPA reduced the full spectrum from 557 bands to 18 key wavelengths, the R2 for EGCG remained about 0.54. These key wavelengths fell mainly in the visible region and the near-infrared overtone absorption bands, supporting the design of portable quality sensors based on discrete wavelengths.【Conclusion】The spectral modeling that integrated physiological priors, such as leaf maturity, could enhance the model’s explanatory power and improve its generalization across cultivars and leaf ages. The single- and multi-task XGBoost frameworks, combined with wavelength selection, provided a theoretical basis and practical support for rapid multi-index tea quality assessment in tea plantations.

    ANIMAL SCIENCE·VETERINARY SCIENCE
    Identification and Functional Analysis of Adipose Tissue Related Long Noncoding RNAs in Response to Cold Stress in Pigs
    XU HaoNan, YAO Lei, PENG Cheng, JIN LinYu, CAO GuoQing, LU Chang, LI BuGao
    Scientia Agricultura Sinica. 2026, 59(18):  4163-4182.  doi:10.3864/j.issn.0578-1752.2026.18.016
    Abstract ( 12 )   HTML ( 1 )   PDF (11485KB) ( 6 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Long non-coding RNAs (lncRNAs) are a class of RNA molecules with a length exceeding 200 nucleotides that generally do not encode proteins, and they play important regulatory roles in cold stress-induced adipose thermogenesis. This study aimed to analyze the transcriptomic data of subcutaneous adipose tissue from Tibetan pigs and Bama pigs to identify lncRNAs involved in regulating adipose metabolism under cold stress, thereby laying a foundation for further elucidation of the mechanisms underlying adipose thermogenesis in response to cold stress.【Method】Five-week-old Tibetan and Bama pigs (n = 8 per breed) were assigned to room temperature or 4 ℃ groups (four biological replicates per group). After 4 h, subcutaneous adipose tissue was collected for RNA-seq. Based on the sequencing data, lncRNAs and mRNAs were identified, and differential expression analysis was performed to compare cold stress responses between the two breeds. Target genes of differentially expressed lncRNAs were predicted via cis-regulation, followed by GO/KEGG enrichment analysis using DAVID. WGCNA was used to screen cold stress-related modules and construct co-expression networks of lncRNAs and mRNAs. Finally, RT-qPCR was performed to validate the transcriptome data of differentially expressed lncRNAs.【Result】A total of 2 988 lncRNAs were identified in the subcutaneous adipose tissue of Tibetan pigs and Bama pigs. These lncRNAs generally exhibited short transcript lengths and a small number of exons. Among them, 1 939 lncRNAs were expressed in at least one sample. Following cold stress treatment, 88 and 10 differentially expressed lncRNAs, as well as 360 and 88 differentially expressed mRNAs, were identified in the subcutaneous adipose tissue of Tibetan pigs and Bama pigs, respectively. Enrichment analysis indicated that the differentially expressed mRNAs were mainly enriched in biological processes such as HPV infection, fatty acid metabolism, and antiviral infection, as well as signaling pathways including the phosphatidylinositol 3-kinase-protein kinase B (PI3K-Akt) pathway, the peroxisome proliferator-activated receptor (PPAR) signaling pathway, and the Ras-related protein 1 (Rap1) signaling pathway. Furthermore, pairs of cis-regulated lncRNAs and protein-coding genes were screened. Among these, 108 lncRNAs were located within 10 kb upstream or downstream of 83 protein-coding genes, and 216 differentially expressed lncRNAs were located within 100 kb upstream or downstream of 289 differentially expressed genes. Finally, a co-expression network of lncRNAs and protein-coding genes was constructed using WGCNA, identifying 18 modules. Among these, the yellow module was associated with the cold stress phenotype. Key lncRNAs within this module, such as MSTRG.8042.2 and MSTRG.8043.2, showed significant associations with core lipid metabolism-related genes. It is speculated that these lncRNAs may serve as key regulatory molecules under cold stress, participating coordinately in the host's cold defense and lipid metabolism regulation. The RT-qPCR results confirmed the reliability of the transcriptome data.【Conclusion】This study identified key lncRNAs including MSTRG.8042.2 and MSTRG.8043.2 that affect cold resistance of hosts under cold stress. The results provide a reference for deciphering the molecular mechanism regulating fat metabolism upon cold stress, and offer a novel perspective for elucidating the cold adaptation mechanism in local pig breeds.

    Genome-Wide Association Study Identifies Candidate Genes Associated with Breast Muscle Water-Holding Capacity in Lindian Chickens
    SUN ChangSheng, WANG Xu, MU Fang, BAI Xue, SU ZhiYong, WANG YuXiang, LI YuMao, ZHANG WenPeng, LUAN Peng, WANG Ning, LENG Li
    Scientia Agricultura Sinica. 2026, 59(18):  4183-4194.  doi:10.3864/j.issn.0578-1752.2026.18.017
    Abstract ( 16 )   HTML ( 1 )   PDF (2156KB) ( 4 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Background】Water-holding capacity of chicken meat is an important index to evaluate meat quality, which directly affects meat tenderness, juiciness, processing loss and commodity value. As a high-quality local chicken breed in Northeast China, Lindian chicken is characterized by strong stress resistance and good meat flavor, but the genetic basis of its breast muscle water-holding capacity traits is still unclear.【Objective】The purpose of this study was to analyze the genetic mechanism of breast muscle water-holding capacity in Lindian chicken based on genome-wide association analysis (GWAS), and to explore relevant molecular markers and candidate genes, so as to provide a theoretical basis for molecular breeding.【Method】In this study, a total of 441 male Lindian chickens at 150 days of age were used. Blood samples were collected from the wing vein, after which the chickens were slaughtered and the left pectoral muscles were collected for the measurement of drip loss, cooking loss, and pressing loss. Genotype data were obtained by 10 × whole genome resequencing, and a high-quality SNP set was obtained by FASTP quality control, BWA alignment, GATK variant calling and PLINK screening. Principal component analysis was used to evaluate the population structure, and significant principal components were incorporated into the model as covariates. Based on GEMMA software, linear mixed model and multivariate linear mixed model were used for single-trait and multi-trait GWAS, respectively, and Bonferroni correction was used to determine the genome-wide significance threshold. The genes in the 40 kb region upstream and downstream of the significant SNPs were annotated, and combined with GO annotation and KEGG enrichment analysis, the biological functions and regulatory networks of candidate genes were systematically analyzed.【Result】The average drip loss, cooking loss and pressing loss of breast muscle of Lindian chicken were 2.92 %, 20.14 % and 18.73 %, respectively, and the phenotypes of the three were approximately normally distributed. Correlation analysis showed that there was a significant positive correlation between the three water-holding capacity traits. The correlation coefficient between drip loss and pressing loss was 0.996 (P < 0.001), the correlation coefficient between cooking loss and pressing loss was 0.975 (P < 0.001), and the correlation coefficient between drip loss and cooking loss was 0.964 (P < 0.001). A total of 28 SNPs significantly associated with water-holding capacity were identified by single-trait GWAS, of which 7 SNPs significantly associated with drip loss were annotated to ADTRP and EXT1 genes. Seven SNPs significantly associated with cooking loss were annotated to SLC25A21, CMPK1, CIB2 and DPY19L3 genes. The 14 SNPs significantly associated with pressing loss were annotated to 6 genes, including PCDH9, NUDT13, CDH2 and CALN1. A total of 28 SNPs significantly associated with breast muscle water-holding capacity were identified and annotated to 9 genes such as ILDR1, EXT1, and SLC25A21. Combining the two analysis results, it was found that EXT1, CMPK1, SLC25A21 and CDH2 were identified in both single-trait and multi-trait GWAS. Functional enrichment analysis showed that candidate genes were mainly enriched in functional items and signaling pathways such as nucleotide metabolism, intercellular adhesion regulation, glycosaminoglycan biosynthesis-heparan sulfate / heparin pathway, nucleotide metabolism pathway and cell adhesion molecules.【Conclusion】In this study, combined single-trait and multi-trait GWAS identified EXT1, CMPK1, SLC25A21 and CDH2 as candidate genes associated with breast muscle water-holding capacity in Lindian chickens. These genes may serve as important candidate genes for improving breast muscle water-holding capacity in Lindian chicken breeding.