【Objective】 Mining the key drought-resistant genes of maize, revealing its drought-resistant molecular mechanism, and providing genetic resources and theoretical guidance for the cultivation of new drought-resistant maize varieties.【Method】Transcriptome data combined with weighted gene co-expression network (WGCNA) and screening methods for physiological and biochemical indicators of drought resistance were used to identify ZmPAL genes associated with drought resistance and rewatering. Genome-wide analysis of the genes encoding PAL was performed using bioinformatics methods. Quantitative real-time fluorescence PCR (qRT-PCR) was used to detect the expression of ZmPAL genes under drought treatment conditions, as well as the expression characteristics of ZmPAL5 among different inbred lines and the expression patterns in different tissues. Finally, genetic transformation was used to analyze the drought resistance function of ZmPAL5 in maize, and the deletion-type Arabidopsis mutant was analyzed for drought resistance with the help of CRISPR/Cas9 technology for the PAL5 homologous gene.【Result】Nineteen maize ZmPAL genes were identified, six of which were clustered on chromosome 5 and encoded proteins that were mostly hydrophilic acidic proteins and relatively evolutionarily conserved in the PAL family of genes. The promoter region of ZmPAL genes contained a large number of cis-acting elements associated with hormonal and abiotic stress responses. Six core genes were identified, four of which were significantly up-regulated for expression after drought treatment. In particular, ZmPAL5 showed an 8.57-fold increase in expression after drought stress. The expression level of ZmPAL5 was found to be significantly higher in the drought-resistant inbred line Zheng 8713 than in the drought-sensitive inbred line B73 under both drought stress and rewatering treatments. Meanwhile, ZmPAL5, a constitutively expressed gene, showed a high level of expression in young stems. Overexpressed ZmPAL5 maize grew well under drought stress, and its relative water content, lignin, chlorophyll, soluble protein, proline content, and activities of superoxide dismutase, peroxidase, catalase, and ascorbate peroxidase were 1.52, 1.49, 1.47, 1.43, 1.44, 1.41, 1.53, 1.41, and 1.35 times, but the malondialdehyde content was 0.65 times that of the wild type. The PAL5-deficient Arabidopsis mutant was sensitive to drought. Under drought stress, its physiological and biochemical indexes showed the opposite trend to those of overexpression of ZmPAL5 maize. 【Conclusion】 Six core genes (ZmPAL3, ZmPAL5, ZmPAL6, ZmPAL8, ZmPAL11, and ZmPAL13) were screened in response to drought stress, in which the expression of ZmPAL5 was positively correlated with drought resistance. ZmPAL5 positively regulated the drought resistance and resilience of the plant by influencing the content of osmotically regulated substances and antioxidant enzyme activities.

【Objective】 In order to provide theoretical basis for molecular breeding of soybean drought resistance, the different evaluation indexes of drought resistance were comprehensively used to screen soybean germplasm with drought-resistance, and the candidate drought-tolerant genes were identified. 【Method】 In 2018, 2019, 2020 and 2021, a total of 188 soybean germplasm were used to determine pod number per plant, biomass per plant and yield per plant under well-watered and drought stressed conditions. Drought resistance index (DI), improved drought resistance index (IDI), weighted drought resistance coefficient (WDC) and weighted drought resistance index (WDI) were used to identify drought resistance of soybean germplasm. The single nucleotide polymorphisms (SNP) loci significantly associated with these parameters were detected by genome-wide association study (GWAS), and the candidate genes for drought resistance were screened by RNA-seq and qRT-PCR analysis of soybean seedling leaves under drought stress. 【Result】 The DI, IDI, WDC and WDI of 188 soybean germplasm varied widely, and five classification criteria for each drought resistance parameter were determined by hierarchical classification method. Among them, Liaodou 15, Liaodou 69, Liaodou 14, Jinzhangzi Huangdou, Zhonghuang 606, Kexin 3 and Koreane 4 were identified as first-grade drought resistant by all evaluation methods. By using GWAS for DI, IDI, WDC and WDI, a total of 15 significantly SNP loci were detected under multiple environments, and the contribution rate of these loci to phenotypic variation ranged from 12.46% to 25.60%. There are 226 annotated genes within 200 kb intervals of upstream and downstream for the significant SNP loci. According to RNA-seq and qRT-PCR analysis of drought-resistant cultivar Liaodou 14 and drought-sensitive cultivar Liaodou 21 under drought stress, a total of 32 annotated genes were significantly differentially expressed by drought stress. Among them, eight genes including Glyma.02G182900, Glyma.04G012400, Glyma.06G258900, Glyma.15G100900, Glyma.01G172600, Glyma.04G012300, Glyma.01G172200 and Glyma.04G010300, encodes calcium-dependent protein kinase, universal stress protein A-like protein, G-type lectin S-receptor-like serine/threonine-protein kinase, protein phosphatase 2C, isoflavone reductase, isoflavone reductase homolog, auxin-like protein, and bZIP transcription factor, respectively. 【Conclusion】 Seven germplasm were identified from 188 soybean germplasm by comprehensive application of different drought tolerance parameters. A total of 15 SNP loci significantly associated with drought tolerance parameters were identified by GWAS, and eight candidate genes were identified.

【Objective】Based on the current “dual carbon” strategic goal, this study aimed to clarify the current characteristics, spatio-temporal pattern and influencing factors of agricultural net carbon sink, so as to provide the important support for accelerating agricultural sink increase and emission reduction.【Method】Based on the scientific reconstruction of the index system, the carbon sink/carbon emission factor method was used to measure and analyze the current situation of China’s agricultural net carbon sink. Then the spatial autocorrelation model was used to discuss the spatial dependence and spatial heterogeneity. Finally, the least-squares method was used to analyze the main factors affecting the change of its intensity. 【Result】From 2005 to 2022, the total amount of agricultural net carbon sink in China was in an obvious upward trend, although there were some interannual fluctuations, and its evolutionary characteristics could be roughly divided into four stages, namely, “continuous rise”, “fluctuating decline”, “rapid rise”, and “slow rise”; the intensity of agricultural net carbon sink was also in an obvious upward trend, with only a slight difference in the trajectory of the evolution, and the difference in its growth rate could be roughly categorized into four stages: “continuous rapid growth”, “slow growth”, “fluctuating ups and downs”, and “slow growth”. 2022, the amount of agricultural net carbon sink had a large interprovincial difference, with Inner Mongolia being the first and Shanghai being the last, and compared with the year of 2005, all the provinces had a significant increase. In 2022, the net carbon sink intensity of agriculture would be the highest in Henan and the lowest in Qinghai, with all provinces showing different degrees of increase compared with 2005. China’s provincial agricultural net carbon sink intensity as a whole showed obvious spatial dependence, but there was also a local spatial clustering phenomenon, more than 70% of the provinces showed obvious spatial clustering characteristics, and the number of provinces located in the high-high clustering and the low-low clustering was approaching. The structure of arable land use, urbanization level, rural residents' income level and the internal industrial structure of agriculture all had a significant impact on the intensity of agricultural net carbon sink; specifically, the higher the ratio of sown area of grain crops, or the higher the urbanization rate, or the higher the income level of rural residents, or the larger the ratio of plantation industry to animal husbandry, the higher the intensity of net carbon sink in agriculture.【Conclusion】The total amount and intensity of China’s agricultural net carbon sink were in a fluctuating upward trend and there were obvious inter-provincial differences. The intensity of China’s agricultural net carbon sink showed obvious spatial dependence and spatial heterogeneity. The intensity of the agricultural net carbon sink was affected by the structure of arable land use, the level of urbanization, the level of rural residents' income, and the structure of the internal industries of agriculture. The measures should be taken to promote the enhancement of sink and emission reductions and to promote the enhancement of agricultural net carbon sink in agriculture, such as establishing a sound policy support system for the development of low-carbon agriculture, strengthening inter-provincial exchanges and cooperation, and increasing financial support for agriculture.

【Objective】Reasonably increasing planting density combined with appropriate nitrogen (N) application rate is an important technical approach for increasing maize yield and resource use efficiency. Understanding the interactive effects of planting density and N rate on maize growth, evapotranspiration (ET) and water use efficiency (WUE) during the growing season, could provide a basis for improving its use efficiency when increasing planting density and controlling N input in maize production. 【Method】Field experiments were conducted during 2022 to 2023 in Jilin Province. Two maize cultivars, Liangyu 99 (LY99) and Demeiya 3 (DMY3), were used in this study. Three planting densities of 50 000, 70 000 and 90 000 plants/hm² and four N application rates of 0, 100, 200 and 300 kg N·hm^-2 were designed to investigate the effects of planting density and N application rate on grain yield and water productivity of different maize cultivars, as well as the dry matter (DM), soil water content, ET and WUE at various growth stages. 【Result】Planting density significantly affected DM and grain yield of maize, but the response trends varied between cultivars. Grain yields of LY99 with 70 000 plants/hm² was 11.1% and 18.3% higher than that with 50 000 and 90 000 plants/hm², respectively. The average yield of DMY3 planted with 70 000 plants/hm² and 90 000 plants/hm² was 10.5% and 9.3% higher than that of 50 000 plants/hm², respectively. Nitrogen fertilization significantly increased DM and grain yield of maize, and also showed significant interactive effects with cultivar or planting density. Compared with N0, grain yields of LY99 were increased by 38.0% to 60.7% under N1, and the yield increases for DMY3 were 24.4% to 38.2%. Notably, the yield responses to N rates were more pronounced for LY99 compared with DMY3. For both cultivars, the yield differences between low N rate and high N rate enlarged with increasing planting density, with LY99 showing a more distinct performance. The water consumption and utilization of maize plants were also significantly affected by planting density, N rate and their interaction. During the growing season, the total ET of DMY3 continually increased with increasing density, while that of LY99 showed the highest values with 70 000 plants/hm² among different densities. In each density condition, the ET of both cultivars increased with increasing N application rates. The WUE of maize plants showed complex responses to planting density and N rate at different growth stages, due to the varied annual precipitation and distribution patterns. The average increase of water productivity of LY99 under planting 50 000 and 70 000 plants/hm² was 8.6% and 10.4% compared with 90 000 plants/hm² respectively_. DMY3 had the highest water productivity when planting 70 000 plants/hm², which increased by 5.8% and 5.3% compared with 50 000 and 90 000 plants/hm², respectively. The water productivity showed different responses to N rate among the three densities. In general, the difference of nitrogen application under low density was small, but it increased significantly under medium and high density. Compared wtih DMY3, LY99 showed higher increases for water productivity when N fertilizer was applied under medium and high density conditions. The correlation analysis showed that interactive effects of planting density and N rate significantly affected maize yield and water productivity by influencing the water utilization at various growth stages. 【Conclusion】Planting density and N rate had significant interactive effects on maize yield and water utilization in the rain-fed region of Northeast China. The two maize cultivars used in this study could obtain high grain yield and water productivity under a moderately higher density of 70 000 plants/hm² combined with 200 kg N·hm^-2 rate.

【Objective】Exploring the genetic loci and related genes that control cottonseed size traits to lay a foundation for subsequent study on the molecular mechanism cottonseed size formation. 【Method】The upland cotton recombinant inbred line (RIL) population composed of 300 lines was used as the research material. Seven phenotypic traits including cottonseed index (SI), seed length-cutting acreage (SLA), seed length-cutting perimeter (SLP), seed length (SL), seed width (SW), length-width ratio (LWR) and seed roundness (SR) were evaluated in four environments. The RIL population was genotyped by liquid phase chip strategy. The high-quality single nucleotide polymorphism (SNP) markers and phenotypic data were subjected to perform genome-wide association study (GWAS), and quantitative trait nucleotides (QTNs) associated with cottonseed size-related traits were mined. The genetic effects of QTNs were analyzed to identify candidate genes. 【Result】Seven cottonseed size-related traits showed a continuous normal distribution in four environments, which expressed a sizable phenotypic variation. The coefficient of variation ranged from 1.82% to 10.70%. The influencing effect on trait formation were basically as genotype>environment>genotype × environment, indicating suitability for GWAS analysis of these results. Correlation analysis showed that the seed index was significantly correlated with SLA, SLP, SL and SW, and LWR was significantly correlated with SR, indicating the possible existence of pleiotropic loci. GWAS was performed using the 3VmrMLM model, and a total of 47 QTNs were associated with these seven traits. A total of 11 QTNs were associated on chromosome A07, of which three physical loci in the region of 71.99-72.87 Mb, A07:71993462, A07:72067994 and A07:72198802 were very close and simultaneously associated with SI, SLA, SLP, SL and SW in four environments. The average value of R2 between markers was>0.8 (P<0.001), showing a large linkage disequilibrium. Genetic effect analysis showed that there were two haplotypes in this region. Among these cottonseed size relating traits, haplotype Ⅱ and haplotype I were significantly different, indicating that these loci directly affected cottonseed size traits and could be used for molecular marker-assisted selection. The expression patterns of the genes in the interval were analyzed using TM-1 transcriptome data. The results revealed that Gh_A07G1767 was preferentially expressed and Gh_A07G1766 specifically expressed at the stage of cottonseed development. These results speculated that these genes may play an important role in the growth and development of cottonseed.【Conclusion】47 QTNs were identified, and two candidate genes related to cottonseed development were screened.

【Objective】An accurate and rapid indoor evaluation system was established by using soybeans with different resistance levels to Phomopsis seed decay as test materials. And then 170 soybean germplasm accessions were employed to screened out disease-resistant varieties, so as to provide methods and material basis for high-throughput assessment of Phomopsis seed decay in soybean and cultivation of resistant varieties.【Method】In terms of establishing a reliable evaluation method for Phomopsis seed decay, Qihuang 34, Williams, Zhongzuo 09-560, z13-631-2, ZDD26268, Chenxiqingpidou 1 and Tongxianhuangdou were selected as experimental materials. For each soybean accession, the seeds with uniform size and undamaged seed coat were germinated in the dark after disinfection. At different germination stages, the pathogen of Phomopsis seed decay was inoculated for 24 h, 48 h, 72 h and 96 h. The mycelium coverage rate and seed decay rate of seed surface under different infection time were counted to determine the optimal identification period for evaluating Phomopsis seed decay in soybean. Then, the resistance of 170 different soybean germplasms in natural population was identified by using the coverage rate of mycelium on the surface of seeds and the decay rate of seeds as evaluation indexes. The high disease resistance varieties were screened based on 5 disease resistance levels.【Result】The soybean accessions showed the most significant differences in disease resistance levels after 96 h of germination when mycelium coverage rate and seed decay rate of soybean surface were used as evaluation indexes. Further comparison of the incidence of 24 h, 48 h, 72 h and 96 h after infection showed that the difference in disease resistance between different varieties after infection for 72 h was the most obvious. Therefore, it was the most suitable period, 72 h of infection at the bud stage after 96 h of germination, for evaluating the resistance level of different soybean varieties to Phomopsis seed decay. The resistance of 170 soybean varieties to Phomopsis seed decay was identified and classified into five disease resistance grades, namely, high resistance, medium resistance, medium susceptibility, susceptibility and high susceptibility. Among them, there were 30 varieties of grade I (high resistance to disease), 51 varieties of grade Ⅱ (medium resistance to disease), 71 varieties of grade Ⅲ (medium disease susceptibility), 4 varieties of grade Ⅳ (disease susceptibility) and 14 varieties of grade V (high disease susceptibility), idicating that there are extensive variations in the resistance to Phomopsis seed decay of soybean germplasm resources in China.【Conclusion】In this study, the most optimum stage of disease identification was considered as soybean seeds after 96 h germination to infect the Phomopsis longicolla for 72 h. After that, the mycelium coverage rate and seed decay rate of soybean surface were counted as evaluation parameters. The evaluation system has high accuracy and reliability, which can provide an effective method for high-throughput identification of different varieties in the laboratory. And 30 highly resistant varieties were further screened to provide a material basis for the breeding of resistant varieties.

【Objective】 Analyzing the evolutionary patterns of genetic diversity of wheat germplasm resources, providing richer and more diverse original parental materials for parental selection and variety selection in wheat breeding in Shanxi Province.【Method】 Using 323 landraces and 105 cultivated varieties as natural populations, a 55K SNP chip was used to perform whole genome scanning on 428 natural populations, analyzing genetic diversity, genetic structure, principal components, genetic clustering, and phylogenetic relationships among varieties. 【Result】 The distribution of SNP loci on 21 chromosomes ranged from 329 to 1 639, with an average of 1 152. The distribution range of 7 partially homologous groups is 2 154-3 852, with an average of approximately 3 456. The distribution pattern of the genome is: B genome>A genome>D genome. Genomic annotation polymorphism markers have the highest distribution among gene regions, accounting for about 50%. Analysis shows that SNP loci cover 21 chromosomes, 7 homologous groups and 3 genomes, but their distribution varies, with a polymorphism rate of 45.60%. The average observed heterozygosity of the entire population (0.0185) was lower than the expected heterozygosity (0.4992). The changes in the average shannon wiener index and polymorphism information content of the entire natural population were not significant. Comparing the diversity parameters of natural populations, it was found that the genetic diversity of the population is not high, the genetic diversity of cultivated varieties is slightly higher than that landraces. The population structure analysis of natural populations divides the population into two major groups. Group I has 307 materials, mainly landraces. Group Ⅱ has 121 materials, mainly cultivated varieties. The natural groups were divided into five groups by both principal component and cluster analysis. The average genetic distance between the varieties in group I is 0.21831, with a range of 0.00127-0.72461. The average genetic distance between varieties in group Ⅱ is 0.14619, with a range of 0.00038-0.76489. The varieties in group Ⅲ the average genetic distance between the varieties of group Ⅳ is 0.16521, with a range of 0.00049-0.43033. The average genetic distance between varieties of group Ⅳ is 0.17643, with a range of 0.00118-0.60496. The average genetic distance between varieties of group V is 0.12039, with a range of 0.00042-0.37032. It can be seen that the variation of genetic distance between wheat varieties is large in Shanxi Province. However, the average genetic distance value is low, the clustering classification differentiation is obvious. The genetic relationship between varieties in the middle of the group is relatively close. Comparison shows that the average genetic distance of group I and group Ⅳ is higher than that of group Ⅱ, group Ⅲ and group V. The genetic distance variation of group I and group Ⅳ is higher than that of group Ⅲ and group V. It can be seen that the genetic distance of cultivated varieties is generally greater than that of landraces.【Conclusion】 The 55K SNP chip was used to analyze the genetic diversity of Shanxi winter wheat germplasm resources, clarifying the distribution characteristics of genetic diversity at the genomic level between Shanxi wheat cultivated varieties and landraces. The introduction of exogenous genes into cultivated varieties is beneficial for improving genetic diversity, while the genetic diversity of landraces is relatively low. At the same time, the genetic relationships of very few varieties are polarized, so it should be rationally used differently in subsequent utilization.

【Objective】Analyzing the yield and yield related traits of Sichuan wheat varieties from 2000 to 2020, providing reference for genetic improvement of yield in Sichuan wheat varieties. 【Method】From 2019 to 2022, a community trial design was used to measure the yield and related traits of 145 wheat varieties in Sichuan Province since 2001 to 2016, as well as 60 high-yield wheat varieties (Varieties with top yields in regional trials in Sichuan Province over the years) since 2000 to 2020. This data was used to analyze the trend of yield and yield related trait changes in Sichuan wheat cultivars cultivated from 2000 to 2020. 【Result】145 Sichuan wheat varieties from 2001 to 2016 have an average annual genetic gain of 37.20 kg·hm^-2 or 0.66% in yield. Grain number per spike and effective spike number per unit area showed an increasing trend, while thousand grain weight and plant height showed a decreasing trend. Correlation analysis showed that effective spike number per unit area was positively correlated with yield. Path analysis showed that the continuous increase of effective spike number per unit area (annual increase 0.42×10⁴/hm² or 0.13%) was the main factor for the increase of yield potential of high-yielding varieties. The average annual yield genetic gain of 60 high-yield wheat varieties from 2000 to 2020 was 61.10 kg·hm^-2 or 0.89%, the effective spike number per unit area showed an increasing trend, the plant height showed a decreasing trend, and the grain number per spike and thousand grain weight had almost no change. Correlation analysis shows that there was a significant positive correlation between yield and the number of effective ears per unit area. Path analysis showed that the continuous increase in effective spike number per unit area (with an average annual increase of 1.80×10⁴/hm² or 0.51%) was also a major factor in improving the yield potential of 60 high-yield wheat varieties in Sichuan from 2000 to 2020. 【Conclusion】The improvement and breeding of wheat yield heritage in Sichuan Province has made some progress, especially the improvement effect of high yield breeding is remarkable, and the yield level of wheat varieties in Sichuan Province is gradually increasing. The continuous increase in effective ears per unit area was the main factor for improving the yield potential of Sichuan wheat varieties. High grain number per spike and thousand grain weight are important foundations for high yield in Sichuan wheat, but their genetic improvement is in a bottleneck period. Increasing the effective spike number per unit area is the key to furtherly improve the yield of wheat in Sichuan.

【Objective】Soluble sugar content is one of the important quality traits of vegetable soybean. The genetic variation and genetic mechanism of soluble sugar content in fresh soybean seeds were studied to provide a basis for germplasm innovation and quality breeding of vegetable soybean.【Method】Using 133 soybean landraces from the Northeast region, the North region, the Huanghuaihai region and the South region, the soluble sugar content of fresh seeds was determined in the three environments of Lianjiang in spring, Fuqing in spring and autumn, in 2021. In combination with 82 187 high-quality SNP markers, the whole genome association analysis of soluble sugar content was conducted based on the mixed linear model MLM (Q+K), and the SNP loci with significantly related to soluble sugar content were identified. The candidate intervals were selected by the significant SNP loci and the extension of 119.07 kb linkage disequilibrium decay distance at both ends. The candidate genes were predicted according to the annotation and tissue expression information of the genes in the candidate intervals.【Result】The variation range of soluble sugar content in fresh seeds under three environments was 3.37-33.84 mg·g^-1, the genetic variation coefficient was 24.59%-32.69%, and the heritability of soluble sugar content was 68.14%. 6, 8 and 22 SNPs were significantly associated with the soluble sugar content of fresh seeds were detected in Lianjiang in spring, Fuqing in spring and autumn, respectively, and phenotypic variation was 12.43%-29.27%. A total of 86 genes were obtained in the candidate regions of 9 significant SNP loci with higher interpretation rate of phenotypic variation, and 9 candidate genes were further screened by gene annotation and tissue expression information. These candidate genes are mainly involved in biological processes such as transcription factors, glycoprotein families and carbohydrate synthesis and transport. Among them, Glyma.01g016500, Glyma.13g042100, Glyma.16g131800 and Glyma.16g155300 were more highly expressed in soybean seeds and pods, which can be used as the most potential candidate genes for soluble sugar in fresh soybean seeds.【Conclusion】Through genome-wide association analysis, 36 SNPs significantly associated with soluble sugar content in fresh seeds were detected, and 9 candidate genes were further screened out, which may be involved in the regulation of soluble sugar content in fresh soybean seeds. Among them, Glyma.01g016500, Glyma.13g042100, Glyma.16g131800 and Glyma.16g155300 can be the key candidate genes for regulating soluble sugar content in fresh soybean seeds.

Seed shattering is a major factor limiting rice production, and breeding new rice varieties with moderate seed shattering is a key challenge faced by rice breeders worldwide. Rice is the most important cereal crop in China, plays a vital role for national food security. Seed shattering is one of the most important traits during rice domestication, and the abscission zone is the important region to control seed shattering. Compared with wild rice, cultivar has eliminated the seed shattering with partially developed abscission layer. Seed shattering not only has a direct impact on the yield, but also affects the way of its mechanical harvest. In order to breed rice varieties with moderate seed shattering in agricultural production, it is necessary to mine and utilize important seed shattering genes and introduce them into excellent rice varieties for genetic improvement, so as to breed new rice varieties suitable for mechanical harvesting with moderate seed shattering. Several seed shattering genes had been identified by map-based cloning, such as SH4/SHA1, qSH1, OsSh1/ObSH3, and their functional mechanisms had been analyzed. At the same time, new rice materials with moderate seed shattering have been successfully developed through CRISPR/Cas9 gene editing technology, gamma ray mutagenesis technology and gene introduction methods. Seed shattering has an important effect on grain yield and rice harvesting methods, in this paper, we reviewed the methods, physiologic basis, the identification of seed shattering genes and genetic mechanism of seed shattering in rice. At the same time, it is proposed that by using the important genes in excellent rice germplasm resources, could provide reference for exploring the mechanism of rice seed shattering, and breed new rice varieties suitable for mechanical harvesting with moderate seed shattering.

【Objective】Taikemai33, derived from a cross between Zhengmai366 and Huaiyin9908, is a new released wheat cultivar with high quality, high yield, and excellent disease resistance, which has a broad genetic base, and a high potential for application in wheat production. The objective of this study is to dissect the genetic composition of Taikemai33 to provide information for parental selection to use this cultivar to develop more new wheat cultivars. 【Method】Taikemai33 and its pedigree parents including Zhengmai366, Huaiyin9908, Yumai47, PH82-2-2, Yumai13, Yumai 2 hao, Bainong3217, Yanda24, Xiannong39, Fengchan 3 hao and Funo were screened using the 55K wheat SNP chip to dissect the genomic composition of Taikemai33 to evaluate the genetic contributions of each parental line to Taikemai33. 【Result】The similarity coefficient between Taikemai33 and its pedigree parents ranged from 0.72 to 0.93, and the genetic composition of Taikemai33 was highly similar to Zhengmai366, the pedigree mother parent, with a genetic similarity coefficient of 0.93. SNP marker analysis showed that the pedigree parents contributed different proportion to the genome of Taikemai33, with the pedigree mother contributed 66.57%, whereas the pedigree father contributed 33.43%, indicating that Taikemai 33 inherits more genetic materials from the maternal lineage. Furthermore, the pedigree mother contributed 71.0%, 85.0% and 49.4% to subgenome A, B and D of Taikeimai33, whereas those were 29.0%, 15.0% and 50.6% contributed by the pedigree father. For each chromosome, the pedigree mother contributed more on chromosome 1A, 2A, 3A, 4A, 7A, 1B to 7B, 1D and 2D, whereas the pedigree father contributed more on chromosome 5A, 4D, 6D and 7D. The contributions of the pedigree parents on 6A, 3D and 5D were equal. Taikemai33 genotype map showed that the contribution loci of the pedigree mother were distributed in clusters on chromosome 1A, 5A, 7A, 2B, 7B, 2D, with those from the pedigree father were on chromosome 4A, 5A, 6D, 7D. Interestingly, among the polymorphic SNP loci, between Zhenmai366 and Huaiyin9908, Taikemai33 showed 109 loci that were absent in both parents, distributing on 19 chromosomes except 1A and 6A. Chromosome 4A, 2B, 6B and 7D of Taikemai33 confer most of the polymorphic SNPs in clusters with cluster number of 10, 9, 11, and 9. 【Conclusion】We constructed the genotype map and dissected the genetic composition of Taikemai33, determined the loci contributed by the pedigree parents and identified that Taikemai33 inherited more genetic materials from the pedigree mother and conferring some specific loci different with the pedigree parents.

【Objective】FWL (Fruit Weight2.2-Like) gene is a negative regulator of cell proliferation, which not only regulates plant organogenesis and organ size, but also participates in the regulation of metal ion transport accumulation and signal transduction. Analyzing of the function of OsFWL3 gene is helpful to reveal the transport mechanism of trace metal elements in crops. It provides theoretical support for reducing heavy metal accumulation and improving crop quality. 【Method】The gene information, genome structure and phylogenetic tree of OsFWLs family were analyzed by bioinformatics method, and the expression profile of OsFWL3 gene was predicted. Two OsFWL3 knockout lines were obtained using CRISPR/Cas9 gene editing technology. Then wild type and Osfwl3 mutants were treated with ZnSO₄ at seedling stage and filling stage, respectively. The phenotypes of plants and grains after treatment were analyzed, and the content variation of metal elements such as Zn was determined to explore the effects of OsFWL3 on the transport and accumulation of metal ions and seed quality. 【Result】The gene function of OsFWLs family is similar to some extent. OsFWL3 gene is highly expressed in anther and panicle, indicating that it is closely related to reproductive development of rice. The number of primary branches, grian length, grain thickness and 100-grain weight of Osfwl3 mutants are significantly larger than WT. OsFWL3 affects the content and distribution of Zn and other metal ions in rice seedlings and grains. The deletion of OsFWL3 gene affects the competitive transport of Zn, Cd and Mn from underground to above-ground, lower grain to central grain and husk to brown rice. 【Conclusion】OsFWL3 gene affects the distribution of Zn and other metal ions in rice grains and plants, and it plays an important role in regulating the growth and development of rice plants and grain size.

【Objective】 Maize stalk rot is one of the common and most devastating diseases in major maize production areas in China. Under natural conditions, maize stalk rot is mostly caused by the synergistic infection of various pathogens. The purpose of this study was to screen elite maize germplasms resistant to the synergistic infection with F. graminearum and F. verticillioides. Meanwhile, QTN associated with resistance to the infection by synergistic infection of two Fusarium species and the candidate resistance genes were identified, which will provide gene resources and theoretical reference for molecular breeding of maize variety resistant to synergistic infection of different pathogens.【Method】 Using a maize natural population as experimental materials, the maize stalk rot phenotypes were investigated upon simultaneous infection with F. graminearum and F. verticillioides. Furthermore, GWAS analysis was conducted to identify significant resistance SNP and to predict candidate genes for stalk rot resistance.【Result】 Through the phenotypic analysis of stalk rot caused by synergistic infection under both field and laboratory conditions, it was found that inbred lines from different sources and subgroups showed significant phenotypic variation. More specifically, the field assay results showed that the inbred lines collected from China were more resistant and that from USA were more susceptible to synergistic infection. Moreover, the inbred lines of tropical and subtropical subgroup were more resistant, whereas the inbred lines of Mixed subgroup more susceptible. The seedling assay results under laboratory condition showed that the lines collected from USA were more resistant, whereas that from CIMMYT were more susceptible. The lines of SS subgroup were more resistant, and that of Mixed subgroup were more susceptible. By integrating field and laboratory phenotypic data, 29 and 16 lines with higher levels of resistance to synergistic infection were screened out respectively, and 6 resistant lines were identified under both conditions. Moreover, based on the field phenotype GWAS, 18 QTNs associated with the resistance were identified, and 93 candidate genes associated with stalk rot resistance to synergistic infection were mined. Among these candidate genes, four genes showed haplotype variation, whose expression levels were up-regulated in disease-resistant lines.【Conclusion】 Using the natural population of maize with great diversity in genetic background, 6 resistant lines to the synergistic infection with Fusarium spp. were identified under two conditions, which can be used as potential germplasm resources for maize stalk rot resistance in the future. Four candidate genes that might be involved in the resistance to synergistic infection were identified by GWAS, which will provide genetic resources for the breeding maize varieties with enhanced resistance to stalk rot caused by F. graminearum and F. verticillioides.

【Objective】Genetic improvement for efficient utilization of maize nutrients represents a crucial method to ensure national food security. Exploring quantitative trait locus (QTL) and related candidate genes of nitrogen use efficiency can provide a theoretical basis for improving the efficiency of nitrogen fertilizer in maize and cultivating high-yield and high-efficiency maize varieties. 【Method】In this study, QTL mapping analysis in one recombinant inbred line (RIL) population constructed by KA105 and KB024 was performed for grain yield under two different nitrogen treatments, including the derived traits partial factor productivity from applied nitrogen (PFPN), low nitrogen tolerance coefficient (LNTC) and nitrogen agronomic efficiency (NAE). Concurrently, integrating the seedling transcriptome data of the parent KA105 under nitrogen treatment, differentially expressed genes were identified, and candidate genes associated with maize nitrogen use efficiency were mined through co-expression analysis. Subsequently, the selected candidate genes were validated using qRT-PCR. 【Result】Through mapping analysis, a total of 36 QTLs distributed across different chromosomes were detected, explaining 1.63% to 17.26% of the phenotypic variation. Among these, eight major QTLs with a phenotypic variation explanation rate exceeding 10% were identified, along with seven genetically stable QTLs commonly identified across different traits or environments. Notably, qNNGYP1 located on chromosome 1 has been repeatedly detected in previous studies, with a phenotypic explanation rate of up to 11.73%. Additionally, other QTLs (qNNGYP1, qPFPN1) co-located in this interval across different environments, suggesting it as a focal region for further investigation. Combining transcriptome data of seedlings under low nitrogen stress, 39 differentially expressed genes within these QTL intervals were identified, and 6 key genes were identified through co-expression network prediction. The result of qRT-PCR indicated that the expression trends of the candidate genes under both nitrogen treatments were consistent with the transcriptome data. Specifically, GRMZM2G366873 was involved in the regulation of auxin homeostasis and may participate in maize responses to low nitrogen stress, drought stress, and boron stress through auxin signal transduction, also regulating ear length. GRMZM2G414192 was involved in the response of the photosynthetic system to low nitrogen stress and was regulated by brassinosteroids. GRMZM2G414043 was associated with maize grain length and biomass, while GRMZM2G040642 may be involved in the long-distance signal transduction of nitrogen. 【Conclusion】In summary, a total of 36 QTLs were identified, distributed across chromosomes 1, 4, 5, 7, 8, and 9, including eight major QTLs (PVE>10%). The candidate genes GRMZM2G366873, GRMZM2G414192, GRMZM2G414043, and GRMZM2G040642 were identified as potential genes for maize nitrogen efficiency.

With the continuous improvement in rice cultivation techniques, China has maintained a high rice production level of about 210 million tons over the past decade. Direct-seeding rice cultivation technology, recognized for its efficiency and simplicity, has been favored by Chinese farmers. However, controversies persist regarding direct-seeding rice compared to transplanted rice in national-scale production. Thus, this study employed meta-analysis techniques to quantify disparities in grain yield, economic benefit, rice quality, lodging characteristic, and greenhouse gas emissions between direct-seeding rice and transplanted rice. Our results indicated that direct-seeding rice significantly reduced grain yield by an average of 6.3% relative to transplanted rice, which was main due to the reduced total spikelet (-3.8%) and filled-grain percentage (-1.8%). In different planting systems in China, the yield of direct-seeding rice had significantly decreased compared to transplanted rice, and the direct-seeding rice-induced reductions in yield of single rice (-10.9%) and late rice (-13.1%) were higher than those of middle rice (-4.8%) and early rice (-4.4%). The grain yield reductions for direct-seeding rice were from 10% to 20% in Jilin, Liaoning, Xinjiang, Ningxia, Shandong, Jiangsu, and Zhejiang provinces, meanwhile Heilongjiang and Jiangxi provinces saw reductions of 5% to 10%, but it had no significant effect in other provinces. Direct-seeding rice resulted in comparable net economic return relative to transplanted rice (p> 0.05). Direct-seeding rice reduced milled rice rate (-3.1%) and gel consistency (-3.5%), improved appearance quality (chalkiness percentage and chalkiness degree, which decreased by 25.3% and 22.5%, respectively), whereas no significant effects were observed on nutrition quality and taste value. Direct-seeding rice increased lodging index at base of the first (+12.4%) and third (+10.3%) internodes, but not at the second internode, indicating an increase in risk of lodging relative to transplanted rice. In terms of greenhouse gas emissions, direct-seeding rice fields showed reductions in methane emissions (-42.8%), global warming potential (-36.2%), and greenhouse gas intensity (-41.1%) compared to transplanted rice fields, while promoting nitrous oxide emissions (+29.1%). In addition, a review was recounted on nitrogen utilization and its loss, water and energy use efficiency, and weed incidence. Finally, the recommendations for the future advancement of direct-seeding rice were proposed, main focusing on rice variety breeding, rice cultivation technique optimization, rice planting area layout, as well as policies and services with the goal of technological innovation and regionalized application of direct-seeding rice cultivation technology in China.

【Objective】 Starch is the main component of wheat kernel and plays an important role in processing. The gelatinization characteristic of starch is an important index to evaluate its quality. The genetic variation of starch gelatinization was studied to provide basis for improving wheat quality. 【Method】 Seven starch gelatinization traits, including gelatinization temperature, peak time, peak viscosity, trough viscosity, final viscosity, decay value and recovery value, were phenotypically determined in 205 winter wheat varieties. Genome-wide association analysis was performed using 90K chip, and haplotype analysis was performed on the stable and significant sites found. 【Result】 The seven characteristics, such as pasting temperature, showed abundant variation in different environments, and the coefficient of variation of attenuation value was the largest (29.31%-31.14%). There were significant differences among genotype, environment and genotype × environment, and the generalized heritability was 0.69-0.86. Through genome-wide association analysis, we found 198 loci that showed significant associations with seven traits. It was distributed in 20 other linked groups except 6D chromosome. There were 58 sites that were stable in 2 or more environments, involving all 7 traits, such as pasting temperature (10), peak time (5), peak viscosity (12), trough viscosity (10), final viscosity (7), break down (4) and set back (10), which could explain 5.54%-22.21% of genetic variation, twenty-one new sites were identified. By haplotype analysis of multiple effector sites that exist in multiple environments and have high phenotypic contribution, Four haplotypes, Hap1 (66.84%), Hap2 (16.84%), Hap3 (9.70%) and Hap4 (6.63%), were found at Kukri_c17417_407 on chromosome 4A, which were significantly related to peak viscosity and break down. Where Hap2 is the peak viscosity and high break down. (P<0.0001). The distribution frequency of varieties (lines) containing haplotype Hap2 in different ecological regions was from high to low as Huanghuai winter wheat region>foreign varieties>Southwest winter wheat region>Middle and lower reaches of Yangtze River winter wheat region>Northern winter wheat region. There were 11 single cause multieffect sites, among which there were 3 multiple effect sites associated with final viscosity, set back, peak time and trough viscosity. Jagger_c4026_328 and other 11 stable genetic loci located on 1B, 2A, 3A, 3B, 4A, 4B, 5B and 6B were mined, and 11 candidate genes that might be related to wheat starch gelatinization traits were screened. 【Conclusion】 In this study, RVA parameters had high heritability, and the RVA parameters of wheat starch were different in different environments. In this study, RVA parameters had high heritability, and the RVA parameters of wheat starch were different in different environments. 58 stable loci were detected that were significantly associated with starch gelatinization traits, and 4 different haplotypes were identified on chromosome 4A that were significantly associated with peak viscosity and break down, and 11 candidate genes related to starch gelatinization were screened, which could provide help for marker-assisted high-quality wheat breeding.

【Objective】Based on the field staging experiments, this study revealed the different responses of winter wheat and summer maize sowing date adjustments in growth and development, photosynthetic physiological characteristics, grain filling, yield formation and quality to climate warming in northern area of North China, providing scientific basis for agricultural production measures to cope with climate change in North China Plain.【Method】Different sowing date experiments of winter wheat and summer maize have been conducted at Hebei Gucheng Agricultural Meteorology National Observation and Research Station in northern area of North China from 2017 to 2023, which were set up in four sowing dates, including 10 d early sowing, 10 d late sowing, 20 d late sowing, and control. The growth process, above-ground dry matter accumulation and distribution, photosynthetic characteristics of leaves, grain filling rate, yield agronomic traits, and grain nutrients of winter wheat and summer maize were observed. 【Result】 The whole growth period of winter wheat was shortened with the delay of sowing date, mainly because of the shortening of seedling stage before winter. There was a parabolic relationship between the whole growth period and sowing date of summer maize. The seedling stage was shortened by 1.3 d, and the flowering stage and grain formation-filling stage were extended by 1.5 d and 1.6 d for every 10 d delay of sowing date. The grain filling characteristics of winter wheat and summer maize were not sensitive to sowing date adjustments, and the grain filling rate of summer maize was little different during different sowing dates, but the grain formation period, the filling end date and the peak date were successively delayed due to the delay of sowing date, and the duration of filling days was shortened by 4.0 d for every 10 d delay of sowing date. Under the background of warm autumn and winter in northern wheat region, the sowing duration of winter wheat was extended, while its influence on the yield was obviously weakened. The delay of sowing date with the increase of sowing seed amount would result in yield increase slightly. The yield of summer maize decreased significantly with the delay of sowing date, and the decline rate of theoretical yield was 1 381.50 kg·hm^-2 for every 10 d delay of sowing date, but the yield of winter wheat and summer maize showed a jumping decease for 20 d late sowing. The grain distribution rate increased by 1.67% for winter wheat, decreased by 1.57% for summer maize with every 10 d delay of sowing date. As a result, the harvest index increased by 0.017 for winter wheat, and decreased by 0.016 for summer maize with every 10 d delay of sowing date. The leaf photosynthetic rates (Pn) of winter wheat and summer maize were also different in response to sowing date, they were similar for winter wheat during different sowing dates, while decreased by 1.21 μmol·m^-2·s^-1 for summer maize for every 10 d delay after sowing date. Sowing date adjustments had no significant effects on grain quality of winter wheat and summer maize in northern area of North China. 【Conclusion】Extending suitable sowing date range and sowing date delay of winter wheat in North China Plain were positive and effective measures to adapt to climate warming. The early sowing of summer maize in North China Plain might avoid the negative effects of high temperature and heat damage, and would promote the increase of yield.

【Objective】 The heading date plays a crucial role in influencing the regional adaptation and yield of rice (Oryza sativa L.). The identification of early heading genes contributes significantly to enhancing and fine-tuning the regulatory network that controls rice heading, which provides valuable genetic resources for molecular breeding with the goal of achieving early maturity and high yield in rice. 【Method】 CL33, a chromosome segment substitution line with early heading, and 93-11, its recipient parent with late heading, were used as research materials to investigate and analyze their major agronomic traits. Two DNA pools were constructed, comprising plants exhibiting extremely early and late heading. Whole-genome resequencing and BSA-Seq analyses were then conducted to locate the genomic region associated with the heading date. In the subsequent steps, InDel markers within this identified region were developed for fine mapping. The gene LOC-Os08g07740 emerged as the primary candidate gene within localization intervals, determined through gene prediction, candidate gene analyses, as well as references to relevant literatures. This candidate gene was subsequently cloned and analyzed for allelic variations. Moreover, we explored the genetic and phylogenetic relationships of the LOC_Os08g07740 gene within the three rice subgroups, Indica, Japonica and O. rufipogon. This analysis involved studying genomic data within approximately 40 kb upstream and downstream of the gene utilizing bioinformatics software.【Result】 Under both natural long-day and short-day conditions in Nanning, Guangxi, CL33 exhibited a 20-25 days shorter than its recipient parent 93-11. Moreover, under natural long-day conditions, the agronomic traits of CL33 were largely similar to those of 93-11, with the exception of a shortened spike length and a reduced number of grains per spike. Genetic analysis revealed that the early heading trait in CL33 was controlled by a recessive gene. This gene was finely localized within a 100 kb region between the markers PSM8-6 and PSM8-8 on the short arm of rice chromosome 8, encompassing 15 predicted candidate genes. Significantly, the candidate gene ORF13 (LOC_Os08g07740), which shared alleles with known heading date genes like DTH8/Ghd8 emerged as a key candidate. Sequencing and sequence alignment of ORF13 demonstrated an 888 bp coding sequence without introns, encoding a protein of 295 amino acids. Compared to the recipient parent 93-11, LOC_Os08g07740 in CL33 featured a 6 bp insertion and a 9 bp deletion between the 535-536th and 820-821st base pairs, respectively, resulting in consequential amino acid sequence alterations. Hence, it was tentatively named OsEHD8 as the target candidate gene. Genetic evolutionary analyses indicated a significant decrease in genetic diversity within the LOC_Os08g07740 gene in Indica and Japonica compared to O. rufipogon, with a 62.53% decrease in Indica and a 53.76% decrease in Japonica. Nevertheless, the differences in genetic diversity between Indica and Japonica were not statistically significant. The LOC_Os08g07740 gene featured a total of 13 haplotypes, with the Hap_2 possibly representing the common ancestor of the three subgroups. Geographical isolation or environmental differences might have led to the fixation of different haplotypes in the Indica and Japonica subgroups. These findings suggested that the LOC_Os08g07740 gene underwent directional selection in the three subgroups.【Conclusion】 OsEHD8, identified as a functional allele of DTH8/Ghd8, played a key role in promoting early heading in rice under both natural long-day and short-day conditions. Moreover, the chromosomal segment substitution line CL33, which carried the OsEHD8 allele, exhibited no significant differences in other agronomic traits compared to the recipient parent 93-11 under natural long-day conditions, except for a shorter spike length and a reduction in grains per spike.

【Objective】This study aimed to explore the differences in genome-wide DNA methylation patterns and their relationships with associated gene expression in different heat-tolerant sesame varieties under high temperature stress, in order to gain a deeper understanding of the regulatory mechanisms of DNA methylation in sesame's response to high temperature stress, and to provide a theoretical basis for heat tolerance breeding in sesame. 【Method】Two sesame varieties, Zhengtaizhi 3 (heat-tolerant) and Shandong White Sesame (heat-sensitive), were selected as experimental materials and cultivated under high temperature (41 ℃) and control (30 ℃) conditions for 10 days. Nanopore sequencing technology was used to conduct methylation sequencing of the genomic DNA of these two sesame varieties, and transcriptome sequencing was performed to analyze changes in the expression of associated genes. Minimap 2 software was utilized for reference genome sequence alignment, and Tombo software was employed to detect 5mC, CpG, and 6mA methylation sites. Differentially methylated regions (DMRs) were identified based on a genome segmentation approach. Finally, functional annotation and pathway analysis of DMR-associated differentially expressed genes (DMR-DEGs) were conducted using GO, COG, and KEGG databases. 【Result】Under high temperature stress, significant changes were observed in the genome-wide DNA methylation patterns of both Zhengtaizhi 3 and Shandong White Sesame. Specifically, the m6A and cytosine methylation (mC) contents of Zhengtaizhi 3 increased, while those of Shandong White Sesame decreased. A total of 621 DMRs (Zhengtaizhi 3) and 374 DMRs (Shandong White Sesame) were identified across the entire genome, mainly distributed in promoter and intergenic regions. Further analysis revealed that these DMRs were significantly associated with 113 DMR-DEGs (Zhengtaizhi 3) and 56 DMR-DEGs (Shandong White Sesame), respectively, and that demethylated DMRs were closely related to upregulated gene expression. Functional annotation results indicated that these DMR-DEGs were primarily involved in biological processes such as carbohydrate transport and metabolism, posttranslational modification, protein turnover, signal transduction, and secondary metabolite biosynthesis. 【Conclusion】This study revealed the differences in genome-wide DNA methylation patterns and their relationships with associated gene expression in different heat-tolerant sesame varieties under high temperature stress. Zhengtaizhi 3, a heat-tolerant sesame variety, regulated the expression of related genes by increasing DNA methylation levels under high temperature stress, while Shandong White Sesame, a heat-sensitive variety, exhibited a decreasing trend in methylation levels. In particular, the dynamic changes in CpG site methylation played a crucial role in regulating sesame's response to high temperature stress. These findings provide new insights and theoretical support for understanding the mechanisms of sesame heat tolerance and for heat tolerance breeding.

【Objective】 Improving the nitrogen transport and accumulation of nutrient organs after wheat flowering, as well as improving the nitrogen accumulation characteristics of wheat grains, is beneficial for increasing the nitrogen accumulation and yield of wheat. This study aimed to clarify the physiological mechanisms underlying the differences in nitrogen utilization efficiency among wheat plants at different yield levels, so as to provide a theoretical basis for reducing yield differences and improving nitrogen utilization efficiency and to achieve high yield and efficient production of wheat.【Method】 The field experiments were conducted in Shiwang Village, Xiaomeng Town, Shandong Province for two consecutive years from 2020 to 2022. Using Tobacco Farmer 1212 as the planting material, three yield levels of wheat fields with a perennial winter wheat yield of around 10 500 kg·hm^-2 (S), 9 000 kg·hm^-2 (H), and 7 500 kg·hm^-2 (M) were selected to compare and analyze the differences in plant nitrogen accumulation and transport, grain nitrogen accumulation characteristics, and grain yield in wheat fields with different yield levels.【Result】 Compared with H and M wheat fields, S wheat field significantly increased the number of spikes per unit area and thousand grain weight, with grain yield 19.64%-27.91% and 51.68%-80.87% higher than H and M wheat fields, respectively, resulting in the highest nitrogen absorption efficiency, nitrogen partial productivity and nitrogen harvest index. The nitrogen accumulation in the nutrient organs of S wheat field during the flowering period increased by 14.22-42.11 kg·hm^-2 and 53.74-103.16 kg·hm^-2, respectively. The nitrogen accumulation in mature nutrient organs was S, H>M. Compared with H and M wheat fields, S wheat field significantly increased the nitrogen transport of pre flowering nutrient organs and nitrogen accumulation after flowering, and increased the contribution rate of nitrogen accumulation in post flowering nutrient organs to grains, thereby achieving the highest nitrogen accumulation in mature grains. The S wheat field significantly increased the content of free amino acids and soluble proteins in flag leaves after flowering, for promoting the inter source sink transport of nitrogen. According to the logistic equation fitting, the S wheat field significantly increased the maximum and average accumulation rates of grain nitrogen, prolonged the duration of grain nitrogen accumulation, and was the main reason for obtaining the highest grain nitrogen accumulation.【Conclusion】 Suitable soil environment in S wheat field could promote the transfer of stored nitrogen from vegetative organs to grain, and increase the rate of grain nitrogen accumulation and prolong the duration of grain nitrogen accumulation, which was the main reason for the highest grain yield and nitrogen use efficiency.

【Objective】To develop high-yielding and FHB-resistant wheat cultivars in the Yellow and Huai River Valley Winter Wheat Zone (YHWZ), simultaneously improving of yield and resistance was conducted in this study.【Method】Using the elite parent dwarf male sterile (DMS) wheat combined with double haploid (DH) technology and molecular marker assisted selection (MAS) of Fhb1 (DMS wheat molecular breeding strategy), DH lines were developed using Sumai 3 as a donor (FHB- resistant parent) and Zhoumai 16’s DMS wheat, Zhoumai 16, Lunxuan 136 and Lunxuan 6 as recipient parents. The agronomic traits (plant height, heading date, yield, etc.) and FHB resistance were evaluated for these DH lines.【Result】A total of 51 Fhb1-DH lines characterized by facultative growth habit, semi-dwarf and white grains were selected using this strategy. The average number of infected spikelets of 51 lines were 5.7 and 7.3 at the 2020Henan and 2020Beijing sites, respectively, and average disease severities were 27.7% and 35.2%, which is not different from moderately susceptible control Huaimai 20. There was no significant difference in grain yield per hm² between the mean performance of the 51 lines and the control Zhoumai 18. DH116 (Lunxuan 20), a promising line from the 51 lines, was further evaluated for FHB resistance and agronomic traits in multiple environments. The resistance of Lunxuan 20 to FHB was significantly improved, and no significant difference was found in the number of infected spikelets or disease severity between Lunxuan 20 and moderately or highly resistant controls at four sites. Lunxuan 20 showed slightly greater grain yield per hm², and significantly higher number of spikelets per spike and thousand grain weight (P<0.05), earlier heading date and shorter plant height (P<0.05) than the control Zhoumai 18 in two environments. The grain yield per hm² of Lunxuan 20 was 4.6% and 1.7% higher than the control cultivar Bainong 207 in the two list trials of Henan Province, and 3.5% higher than Bainong 207 in the demonstration trial. Resistance of Lunxuan 20 to FHB ranged from moderate susceptibility to moderate resistance in two-year list tests using the single-floret injection and spray inoculation methods. Lunxuan 20 carries the semi-dwarfing gene Rht-D1b at the Rht-D1 locus, and the recessive alleles vrn-A1, vrn-B1 and vrn-D1 associated with the winter growth habit at the Vrn-A1, Vrn-B1 and Vrn-D1 loci. Based on the wheat 660K single nucleotide polymorphisms (SNPs), 64.7% of the SNPs were shared by Lunxuan 20 and its parents, and the direct genetic contributions of Zhoumai 16, Lunxuan 136, Lunxuan 6 and Sumai 3 to Lunxuan 20 were 69.8%, 12.6%, 6.1% and 11.5%, respectively.【Conclusion】A high-yielding and FHB-resistant wheat cultivar Lunxuan 20 was bred using the DMS wheat molecular breeding strategy.

【Objective】Pepper mild mottle virus (PMMoV) is one of the major viruses that harm peppers worldwide. The 126 kDa protein is an important pathogenic factor encoded by PMMoV, but its pathogenic mechanism remains unclear. This study aims to screen for pepper host factors that interact with the 126 kDa protein, and to provide a theoretical basis for elucidating the pathogenic mechanism of PMMoV.【Method】Firstly, the bait vector pGBK-126 kDa was constructed using homologous recombination. Using pepper leaves as experimental materials, total RNA was extracted from pepper leaves using the Trizol method, and a pepper yeast cDNA library was prepared. Subsequently, the cDNA library was screened using pGBK-126 kDa, and the screening results were subjected to sequence alignment and bioinformatics functional analysis using NCBI and Uniprot. Based on the alignment and analysis results, host factors that may participate in plant disease resistance pathways were selected, and their full-length CDSs were cloned and constructed into the pGADT7 vector. Yeast two-hybrid (Y2H), BiFC, and LCI were used to further verify the interaction between 126 kDa and host factors. Finally, the role of transiently overexpressed host proteins during PMMoV infection was analyzed.【Result】High-quality pepper RNA was extracted without degradation. A high-quality yeast cDNA library was obtained, and the bait plasmid pGBK-126 kDa was successfully constructed. A total of 18 pepper host factors that interact with the 126 kDa protein were screened. Bioinformatics analysis revealed that these 18 host factors are widely involved in multiple pathways such as plant enzyme systems, regulation of material and energy metabolism, DNA-binding transcription, hormone synthesis, and defense responses. Among them, three host factors (LA2, PDHE1, BXL1) showed interactions with 126 kDa in one-to-one Y2H interaction verification, indicating the reliability of the initial screening results. The interaction between 126 kDa and BXL1 was further verified in vitro and in vivo using BiFC and LCI. Transient overexpression of BXL1 significantly inhibited PMMoV infection.【Conclusion】The pGBK-126 kDa bait plasmid was successfully constructed. Based on this plasmid, 18 interacting host factors were obtained when screening the yeast cDNA library, which are widely involved in multiple pathways of plant life activities. The screening results were verified to be reliable. Among them, BXL1 interacts with 126 kDa both in vitro and in vivo, and can inhibit PMMoV infection. The results can provide a good theoretical and material basis for further exploring of the infection mechanism of PMMoV.

【Objective】Exploring efficient nodulation soybean germplasm adapted to the ecological conditions of the Bashang area, identifying genetic loci and candidate genes regulating soybean-rhizobium symbiotic nodulation, and improving soybean symbiotic nitrogen fixation efficiency.【Method】This study utilized a natural population of 260 soybean germplasms as the research object, rhizobium strain USDA110 was inoculated under outdoor potted conditions in the Bashang of Hebei Province. The single plant nodule number and single plant nodule dry weight data were used as phenotypic values. Combined with genotype data of the 260 germplasms, a genome-wide association analysis was conducted to explore genes related to soybean-rhizobium symbiotic nodulation.【Result】A total of 18 SNPs significantly associated with soybean nodule number were detected, located on chromosomes 2, 7, 8, 13, 18, and 19. Among them, the significant associated locus BARC_2.01_Chr02_43161654_A_G on chromosome 2 was identified as the main locus controlling soybean nodule number (LOD=3.89). Linkage disequilibrium analysis within the 200 kb interval upstream and downstream of this locus containing BARC_2.01_Chr02_43161654_A_G identified 10 candidate genes regulating soybean nodule number. There was a significant difference in the number of nodules among the materials corresponding to different haplotypes of Glyma.02G243200 (P<0.05), the expression pattern of this gene was queried in the SoyBase database, and it was expressed in root hairs, indicating that Glyma.02G243200 may be a key gene influencing soybean nodule number. Additionally, six SNPs significantly associated with soybean nodule dry weight were identified, located on chromosomes 6, 18, and 20. Among them, the significant associated loci BARC_2.01_Chr06_6069381_G_A and BARC_2.01_Chr06_6192925_T_C on chromosome 6 were identified as the main loci controlling soybean nodule dry weight (LOD=3.49 and LOD=3.35, respectively). Linkage disequilibrium analysis within the 100 kb interval upstream of BARC_2.01_ Chr06_6069381_G_A and downstream of BARC_2.01_Chr06_6192925_T_C identified 14 candidate genes regulating soybean nodule dry weight. Haplotype analysis revealed significant differences in nodule dry weight for the genes Glyma.06G079600 and Glyma.06G079900 between different haplotype materials (P<0.01, P<0.001), the expression pattern of this gene was queried in the SoyBase database, and they were expressed in roots, indicating that these two genes may be key genes influencing soybean nodule dry weight.【Conclusion】This study identified a candidate gene significantly associated with nodule number on chromosome 2 and two candidate genes significantly associated with nodule dry weight on chromosome 6, providing new genetic resources and references for genetic improvement of soybean nodulation traits.

【Objective】 The rice-crab co-culture is a predominant three-dimensional ecological cultivation model in northern Chinese rice fields. The water-soil interface is a key area for material cycling within this system. In this paper, studying the diversity and structure of microbial communities at this interface aimed to understand its evolution and support research on the ecological health of water and soil in rice-crab co-culture systems. 【Method】 Eight long-term rice monoculture systems and eight long-term rice-crab co-culture systems (>20 years) were selected in a typical rice-crab co-cultivation area in Panjin, Liaoning province. Based on measurements of physicochemical indicators of rice field water and soil, as well as high-throughput sequencing technology of the 16S rRNA gene, the study compared the effects of two rice cultivation systems on the properties and bacterial community structures of the water-soil interface in paddy fields. 【Result】 (1) The introduction of crabs significantly reduced the unique microbial communities at the water-soil interface. Specifically, the number of unique operational taxonomic units (OTUs) in surface water and interfacial soil decreased by 27.0% and 71.2%, respectively. However, the introduction of crabs had no significant effect on alpha diversity in surface water, but it reduced the richness of bacterial communities in interfacial soil. (2) The introduction of crabs significantly altered the structure and composition of the water-soil interface bacterial community. The introduction of crabs significantly increased the relative abundance of Proteobacteria (30.4%) and decreased the relative abundance of Acidobacteria (39.9%) in surface water. Simultaneously, it increased the relative abundance of Planctomycetes (21.1%) and decreased the relative abundance of Ignavibacteriae (15.1%) and Nitrospirae (21.7%) in interfacial soil. (3) Proteobacteria and Bacteroidetes were not only core species at the water-soil interface of rice field systems, but also key species in co-occurrence networks, playing important roles in stabilizing ecological networks. (4) The introduction of crabs into rice fields increased the complexity and stability of the bacterial co-occurrence network in the interface soil, but decreased it in the paddy field water. (5) Linear regression analysis showed that NO₃^-- N concentration in surface water and interface soil pH were the main driving factors influencing the diversity and stability of their respective bacterial community structures. 【Conclusion】 The introduction of crabs significantly altered the microbial community structure and diversity at the water-soil interface in rice fields. The increase in nutrient salts in the water posed a risk of reducing the stability of the aquatic microbial community. However, the rice-crab co-culture shaped a more stable bacterial community in the interface soil, which facilitated nutrient cycling and enhanced crop nutrient utilization efficiency.

【Objective】Hardness, a structural feature of seed physical dormancy, is an important trait in soybean domestication. Although hardness is beneficial for seeds to survive in unfavorable environments, it will seriously reduce the emergence rate of soybean in the field, and detrimental to yield and processing quality. Analyzing the QTL and candidate genes using bulked segregant analysis sequencing (BSA-Seq), can provide a theoretical reference for understanding the molecular mechanism of hard seededness in soybean.【Method】The hard seed mutant M_zp661 was obtained from the seeds of Zhongpin 661 induced by ethyl methane sulfonate (EMS), and was crossed with cultivated soybean Zhonghuang 13 (male parent) to construct recombinant inbred line (RIL) population. The progeny lines were investigated for seed hardness, water absorption capacity and anatomical structure of seed coats. Two types of extreme lines in the RIL population, with hard seeds or with imbibed seeds, were selected to construct DNA mixed pools respectively, and then BSA-Seq technology was used to detect genotype differences in extreme-mixed pools and parents. Euclidean distance (ED), delta SNP-index, and delta InDel-index methods were applied to associate hard seed genetic loci of soybean. Combining with bioinformatics analysis, transcriptome data of different soybean tissues and gene annotation information, candidate genes within significant association regions were predicted.【Result】In the progenies of M_zp661, all areas of imbibitive seeds had the penetration ability, and the seed volume increased continuously with the soaking time. However, no changes were observed for hard seeds over 36 hours. With the prolonged of soaking time, the seed coat of hard seeds began to shrink locally and gradually spread to other parts, and finally cotyledons recovered their imbibition ability. The hard seed not only has smooth and compact seed coat, but also has regular network structure of cuticle and thicker palisade layer, while numbers of stomata and loose structures, tiny cracks and thinner palisade layer were existed in the imbibed seeds. These results suggest that the seed hardness of M_zp661 may be caused by the impermeability of the seed coat. ED, delta SNP-index and delta InDel-index association analysis methods not only identified the reported seed physical dormancy locus qHS1, but also simultaneously detected the candidate region Chr.06: 45897227-47746047, which contains a total of 189 genes. Further, transcriptome data and gene annotation predicted that Glyma.06G275300, which is specifically and highly expressed in seeds, might be the candidate gene for this associated region to regulate soybean seed hardness.【Conclusion】Seed hardness of soybean mutant M_zp661 was caused by the impermeability of the seed coat, and Glyma.06G275300 was predicted as a candidate gene affecting the structure of seed coat using BSA-Seq.

【Objective】 This study aims to discover new quantitative trait loci (QTLs) for resistance to the brown planthopper (BPH) in rice by utilizing a genetic population and to assess the pyramiding effects of these QTLs, thereby providing crucial genetic materials and resources for breeding insect-resistant rice. 【Method】 A recombinant inbred line (RIL) population derived from a cross between susceptible NPB and resistant Jiafuzhan was used, combined with seedling stage resistance evaluation and extreme mixed pool sequencing (BSA-seq) to map BPH resistance QTLs. Further, fine mapping, candidate gene identification, pyramiding effect analysis, and mixed-pool transcriptome sequencing (BSR-seq) were employed to elucidate the physiological and molecular mechanisms mediated by these QTLs. 【Result】 BSA-seq identified two major resistance QTLs on chromosomes 1 (30-32 Mb) and 4 (5-7 Mb), named QBPH1 and QBPH4, respectively. Interval-linked molecular markers confirmed the authenticity of these QTLs. QBPH4 coincides with previously cloned BPH3 and BPH15, while QBPH1 represents a newly discovered QTL. High-density markers and recombinant analysis further narrowed the QBPH1 region to 30.61-30.65 Mb. This analysis identified Os01g53294 and Os01g53330 as reliable candidate genes, which code for a respiratory burst oxidase protein B and an anthocyanin 5,3-O-glucosyltransferase, respectively. In evaluations of seedling resistance, brown planthopper honeydew area and weight, insect weight gain, lethality, and preference, comparisons between QBPH1 and QBPH4 single genes and their polymerized lines revealed no significant enhancement in resistance. Both QBPH1 and QBPH4 mediated antibiosis and antixenosis mechanisms with varying effects. BSR-seq analysis highlighted significant enrichment of differentially expressed genes (DEGs) involved in transcriptional regulation, protein phosphorylation, and redox processes among different QBPH1 alleles. Additionally, genes associated with jasmonic acid (JA) synthesis and signaling pathways were significantly upregulated in resistant materials, confirmed by RT-qPCR experiments. 【Conclusion】 A novel BPH resistance QTL, QBPH1, was successfully identified on chromosome 1 of rice. QBPH1 mediated both antibiosis and avoidance against BPH, though its aggregation effect with another QTL was not significant. QBPH1 may mediate defense mechanisms against BPH through involvement in the JA pathway. Based on this, Os01g53294 and Os01g53330 have been validated as reliable candidate genes for QBPH1.

【Objective】 The essence of conventional wheat breeding involves the continuous multi-generation selection and maintenance of hybrid F₁ trait heterosis. Exploring the heterosis and combining ability of lodging-related traits among spring wheat varieties (lines) in the Yellow River irrigation area of Ningxia, and clarifying significantly agronomic indexes related to anti-lodging traits in wheat that can provide a certain theoretical basis for utilizing heterosis in lodging traits of spring wheat, and for the screening and maintenance of its descendants in Ningxia.【Method】 This study utilized 14 spring wheat varieties (lines) from Ningxia Yellow River Irrigation Area as parent materials. According to NCⅡ incomplete diallel crossbreeding design, forty-five combinations were produced. Meanwhile, heterosis, combining ability, and correlation analyses were performed on 13 traits related to lodging resistance, such as plant height and culm type index, in their parents and their F₁ generation.【Result】 Significant differences in lodging resistance trait were observed among different spring wheat varieties (lines). Meanwhile, there is a certain level of heterosis in lodging resistance traits was evident in the F₁ generation. The main stem thrust resistance of NZ42, M6445, and M8887 showed high general combining ability (GCA) values of 11.68, 8.00, and 10.67 respectively. Significant combinations for main stem thrust resistance included M6445×M8887, M6445×MJ48, NZ42×N2038, NZ42×M7723, NZ42×NZ39, and H3015×Ningchun50. Correlation analysis revealed that plant height was significantly positively correlated with the first internode length and second internode weight, and highly significantly positively correlated with the second and fourth internode lengths. Besides, the main stem thrust resistance was highly significantly positively correlated with fresh weight of the main stem, stem diameter, bending moment and culm type index, but significantly negatively correlated with the fourth internode length. However, it was not significantly correlated with the second internode weight and fullness.【Conclusion】 The excellent parent varieties identified were NZ42, M6445, and M8887. Meanwhile, the excellent combinations of NZ42 × NZ39, M6445 × MJ48, and M6445 × M8887 were also observed in this study. In addition, Heterosis was present among most lodging resistance traits in different combinations, but the heterobeltiosis was less pronounced. The resistance to lodging was influenced by both additive and non-additive effects of parental genes, with a greater influence from the maternal genetic background. Furthermore, the plant height and main stem thrust resistance are significantly correlated with the fourth internode length, bending moment, and culm type index, serving as important references for selection of lodging-resistant generations.

【Objective】Based on whole genome identification and analysis of soybean LOX gene family members, to understand the taxonomic evolutionary relationships of each member, to study the expression specificity of each gene member in different tissues and their response to abiotic stress, which provided a theoretical basis for further research on the molecular characteristics, evolutionary process, and function of the LOX gene family. 【Method】Based on the LOX protein sequences of rice and Arabidopsis species in Ensembl database, BLASTP alignment of homologous LOX protein sequences in soybean whole genome database was performed, and MEGA X software was used to construct a phylogenetic tree; Using website MEME for protein conserved motif analysis; Using online software GSDS 2.0 to analyze gene structure; Using TBtools for chromosome localization drawing; Analyze soybean LOX family replication genes using McscanX; Using the PlantCARE website to predict the promoter elements of soybean LOX gene family; Draw gene expression heatmaps of soybean under different tissues and abiotic stress using TBtools, and develop molecular markers for the excellent allele variant GmLOX15A1^-G/A significantly correlated with 100-seed weight.【Result】A total of 43 LOX genes were identified in soybean, unevenly distributed on 13 chromosomes. Collinearity analysis indicates that the GmLOX gene has undergone extensive replication during the evolutionary process. Meanwhile, 39 different types of cis regulatory elements were detected in the LOX gene promoter, indicating that they may be involved in different pathways such as growth and development, light response, stress response, and hormone induction. Expression pattern analysis revealed that the LOX gene has different levels of expression in different tissues of soybean, indicating that members of this family have tissue and spatiotemporal expression specificity. Under drought stress conditions, the GmLOX gene was significantly differentially expressed in soybean roots and leaves (P<0.05). Among them, GmLOX3A3, GmLOX7A1, GmLOX20B1, GmLOX13A1, and GmLOX20A2 were significantly upregulated or downregulated in roots and leaves, suggesting that the GmLOX gene may play an important role in response to stress. At the same time, it was found that GmLOX15A1 is highly expressed in grain tissue and there is an excellent G/A allele variation in the seventh exon of the gene coding region. Molecular markers were developed for this variant site, and the correlation between different haplotypes of GmLOX15A1 and 100-seed weight was analyzed using 1 200 soybean germplasm resources from different ecological regions over a period of 2 years. The results showed that compared to the GmLOX15A1^-A genotype, the average 100-seed weight of soybean germplasm carrying the GmLOX15A1^-G allele gene increased by 2.33 g (P<0.001). 【Conclusion】A total of 43 members of the LOX family were identified in soybeans, which can be divided into 3 subfamilies. The promoter region of the GmLOX gene contains a large number of cis acting elements that respond to hormones and stress, playing different roles in drought stress response. Among them, GmLOX15A1 is highly expressed in grain tissue and there is an excellent G/A allele variation in the seventh exon of the coding region of this gene. Compared with the GmLOX15A1^-A genotype, the average 100-seed weight of soybean germplasm carrying the GmLOX15A1^-G allele gene is significantly increased by 2.33 g. This locus can be used as an excellent haplotype for genetic improvement of soybean grain size.

【Objective】Branching is a significant factor influencing alfalfa yield, and the SPL family of transcription factors represents a crucial class of regulatory genes involved in branching (tillering) development in a diverse range of plants. The objective of this reserch is to investigate the biological function of MsSPL17 in alfalfa and to elucidate the role of MsSPL17 in regulating the development of alfalfa meristems. This will provide a valuable reference for alfalfa high-yield biological breeding.【Method】Bioinformatics methods were used to anlyze MsSPL17 sequences and also constructing a phylogenetic tree. The tissue expression specificity of MsSPL17 in alfalfa was analyzed by real-time fluorescence quantitative PCR (qPCR). The subcellular localization of MsSPL17 protein was determined by tobacco transient expression system, and the transcriptional self-activation activity of MsSPL17 was verified. The transgenic alfalfa was obtained by Agrobacterium-mediated transformation and phenotypic analysis was carried out. Transcriptome analysis was utilized to screen for differentially expressed genes in transgenic lines and validate them for further research.【Result】MsSPL17 contained an open reading frame of 1 011 bp, encoding a protein composed of 366 amino acids, belonging to the SBP protein family. Phylogenetic analysis showed that the evolution of MsSPL17 and its homologous genes was highly similar to the differentiation of species, indicating that it is a functional conserved gene. MsSPL17 expressed in all tissues, including stems, nodes, leaves and tops during the critical period of alfalfa growth and development, implied the necessary regulating function of this gene in alfalfa branching. Subcellular localization assay showed that MsSPL17 protein was localized in the nucleus. Transcriptional self-activation assay showed that MsSPL17 did not have self-activation activity and could be used in interacting proteins screening. MsSPL17 transgenic silenced lines exhibited a notable phenotype, including an increase in branch number and stem node number, a reduction in internode length, and an enhancement in nutritional quality.【Conclusion】MsSPL17 was successfully cloned, it expressed in key tissues of alfalfa branching development. The protein encoded by MsSPL17 was localized in the nucleus and demonstrated no transcriptional self-activation activity. Transgenic lines exhibiting multi-branching traits were obtained, and the number of branches increased significantly in yield, while the crude protein content increased in quality.

Gene editing techniques have made gene edited (GE) organisms enter commercial applications from laboratories. In 2022, the Ministry of Agriculture and Rural Affairs specifically issued the “Guidelines for Safety Evaluation of Genetically Edited Plants for Agricultural Use (Trial)” for the safety evaluation of GE plants without introducing exogenous genes. In 2023, China granted the first biosafety certificate for GE soybean AE15-18-1, marking the official start of the commercialization process of GE crops in China. GE organisms are different from traditional genetically modified organisms (GMOs) containing exogenous DNA sequences, making common GM detection strategies inapplicable to the detection of GE organisms. As the industrialization of GE crops progresses positively, how to efficiently and accurately detect whether a product is gene-edited and its editing characteristics is an important basis for the commercial use and intellectual property protection of GE products. There is an urgent need to develop detection technologies suitable for GE products. With the goal of detecting whether the target sequence has been edited, many detection technologies have been developed based on PCR, sequencing, and other technologies, and are widely used in the screening of GE products in the research and development process. After industrialization, safety supervision and intellectual property protection require not only the detection of whether the sample has been edited but also the rapid identification of the nucleotide sequence characteristics of the sample to determine its origin and identity. Subsequently, precise quantification of the GE components is necessary to determine whether quantitative labeling is required. Currently, it is difficult to quickly identify the identity of GE products with only a few base insertions, deletions, and single nucleotide variations (SNV) using conventional PCR or sequencing technologies. It is even more challenging to accurately quantify the content of GE components. Aiming at the rapid identification of the DNA sequence characteristics after editing and precise quantification, based on the molecular characteristics of GE products, this paper reviews the application of the gel electrophoresis-based PCR method, the sequencing-based method, the real-time PCR-based method, the digital PCR-based method, the editing enzyme-based method, and the instrument-based method in detection of GE organisms, and expounds the advantages and disadvantages of each method during detection. This review initially explores the detection and quantification strategies suitable for GE organisms and provides a reference for subsequent development of detection methods for GE organisms.

【Background】To determine the safe planting limit of winter wheat based on agricultural climate indicators is crucial for the scientific and rational utilization of resources, avoiding freezing disasters, and ensuring stable and high yields of winter wheat. However, in the north of China, which is located in the sensitive area of winter wheat planting, the fluctuation of safe winter wheat planting has been intensified due to the increase of extreme weather events caused by global climate change. It is urgent to clarify the agroclimatic factors affecting the safe planting of winter wheat on a large regional scale and to determine their threshold ranges. 【Objective】The research on the agricultural climatic factors and their thresholds for the safe planting of winter wheat was conducted to provide a scientific basis for the sustainable production and planning of winter wheat in response to climate change. 【Method】The northern China was selected as the research area, which was highly sensitive to the safe planting of winter wheat. Based on the spatial distribution of winter wheat with medium and high spatial resolution and ground meteorological observation data, this research utilized methods such as kernel density estimation, geographic detector to reveal the spatial pattern characteristics of the actual northern limit of winter wheat planting, to quantitatively analyze the influence of agricultural climate factors on the formation of the actual northern limit of winter wheat planting, and to explore the threshold of key climate factors. 【Result】(1) The actual northern limit of winter wheat planting, with a total length of about 2 200 km, fluctuated from southwest to northeast. However, agricultural climate factors exhibited more significant fluctuations along the line of Pingning-Xunyi-Tongchuan-Baishui- Heyang-Hancheng-Jishan. (2) The negative accumulative temperature during winter, average temperature of the coldest month, extreme minimum temperature of the year, and accumulative temperature before winter were crucial factors (q >0.45) in shaping the actual northern limit of winter wheat planting. Agricultural precipitation factors had a minor effect (q <0.19) on winter wheat planting, but interacted strongly with temperature factors (q >0.57). (3) Specific meteorological parameters for the northern limit of winter wheat safe planting in northern China were established: negative accumulated temperature in overwintering period≥-620 ℃·d, coldest monthly mean temperature≥-8 ℃, annual extreme minimum temperature≥-22 ℃, and accumulated temperature before overwintering≥529 ℃·d. (4) The potential northern limit for winter wheat planting has moved about 107 km northward compared to the actual limit, with approximately 23.39×10³ km² of expansion area. 【Conclusion】This study identified the key agricultural climate indicators and thresholds influencing safe winter wheat planting in northern China, which provided a basis for determining potential safe planting areas for winter wheat. The research results could provide the theoretical reference and technical support for how winter wheat planting could adapt to climate change and adjust agricultural planting layout reasonably.

【Objective】 This study aimed to clarify the response of maize growth and yield with different nitrogen application rates to intercropped green manure, so as to provide the basis for establishing maize production pattern with nitrogen fertilizer saving via green manure in arid irrigation area.【Method】 A split plot experiment design was adopted, and two planting patterns of maize intercropping with common vetch (MV) and sole maize (SM) were set up in the main plot. Five nitrogen application rates were set up in the split plot: the recommended nitrogen application rates were N360: 360 kg N·hm^-2, application rates reduced with 25% (N270): 270 kg N·hm^-2, and that reduced with 50% (N180): 180 kg N·hm^-2, reduced with 75% (N90): 90 kg N·hm^-2, and no nitrogen application (N0). Leaf area index (LAI), daily leaf volume (LAD), light energy utilization rate, dry matter accumulation and yield of maize were investigated.【Result】 In 75-150 days after maize emergence, the LAI, total LAD and light use efficiency of maize with reduced nitrogen application rate (N270, N180, N90, N0) were increased by intercropped green manure. The maize LAI, total LAD and light use efficiency under M||V with N270, N180, N90, and N0 were significantly increased by 9.8%, 8.2% and 4.6% than that under SM, respectively. The maize LAI, total LAD and light use efficiency under M||V with N270 was not different with that with N360, in 75-150 days after maize emergence, and which was greater than that of other treatments. Somehow, the maize LAI, total LAD and light use efficiency of SM was decreased with reduced nitrogen application rate. The dry matter accumulation of maize under N270, N180, N90, and N0 was significantly increased by intercropped with green manure than that of sole maize, in 90-150 d days after maize emergence, which was increased by 4.6%. In 75-130 days after maize emergence, the growth rate of maize under M||V was significantly increased by 5.9% compared that with SM. With the same reduced nitrogen application rate (N270, N180, N90, and N0), the maize grain yield under M||V was significantly increased by 8.5% than that under SM. The grain yield of intercropped maize under N270 was significantly higher than that of sole maize under N360. Compared with SM, the kernel number per ear and 1000-grain weight of maize under M||V were significantly increased by 6.0% and 6.3%, respectively. The results of the fitting curve showed that the greatest grain yield of maize intercropped with green manure was 14 876.4 kg·hm^-2 with nitrogen application rate 261.4 kg·hm^-2. While the greatest grain yield of sole maize was 14 012.5 kg·hm^-2 with nitrogen application rate of 348.6 kg·hm^-2. The effect of reducing nitrogen application on the fresh yield and nitrogen accumulation of green manure was significant. The fresh yield of green manure under M||VN270 was significantly higher than that under M||VN360, M||VN90 and M||VN0 by 3.8%, 4.6%, and 9.7%, respectively, while the nitrogen accumulation with M||VN 270 was significantly higher than that with M||VN90 and M||VN0 by 5.3% and 11.9 %, respectively. The total dry matter accumulation of maize with reduced nitrogen applied rate was increased by intercropped green manure, and then the number of grains per ear and thousand grains weight of maize was increased. Therefore, the intercropped maize could obtain a greater level of grain yield with the condition of reduced nitrogen application rate.【Conclusion】 The treatment of 25% reduced nitrogen rate with intercropped green manure could make maize obtain higher yield than sole maize with the recommended nitrogen application rate. This pattern could be used as a recommended production method for chemical nitrogen fertilizer reduced application of maize in the research region and even similar ecological region.

【Objective】Polymerizing soybean high oil genotypes aims at breeding varieties with higher oil content to improve economic efficiency and nutritional value. It is of great significance to increase agricultural output, reduce processing costs and meet global demand for vegetable oil growth.【Method】Glyma.18G027100 C2 gene family was identified by bioinformatic analysis method at the whole genome level. A total of 66 soybean C2 gene family members were identified, named GmC2-01.1-GmC2-20.2 according to chromosome position. Tissue pattern analysis revealed that 7 genes were highly expressed in grains among 66 C2 family genes (GmC2-03.6, GmC2-02.7, GmC2-07.2, GmC2-18.1, GmC2-18.4, GmC2-19.1 and GmC2-20.2). In order to analyze the effect sites of these genes on soybean oil content, SNP sites in the coding regions of these genes were obtained from SFGB database. Correlation analysis of oil content in two years showed that GmC2-18.1 has SNP loci that significantly affect oil content. The genetic diversity of GmC2-18.1 coding region was analyzed by 12 extreme materials. There was a G/A mutation at 2 038 273 bp in coding region of Wm82.a2.v1 version, which regulated seed oil content. It was preliminarily speculated that this gene played a role in seed development or nutrient accumulation. Then, SNP/InDel molecular markers were developed for GmC2-18.1^-G/A gene combined with InDel natural allelic variation site 225 bp upstream of the start codon of GmSWEET39, T/C natural allelic variation site at 8 381 058 bp in coding region of GmST1, A/C natural allelic variation site at the third exon of 41 854 422 bp in coding region of GmMFT. 1 200 soybean germplasm resources from three ecological regions in China were identified by markers in 2 years.【Result】Analysis of variance showed that GmC2-18.1^-G, GmSWEET39^-Deletion, GmST1^-T and GmMFT^-A significantly increased oil content by 1.72, 1.95, 1.58 and 2.06 percentage points (P<0.01). The results showed that the average oil content of soybean seeds carrying GmC2-18.1^-G, GmSWEET39^-Deletion, GmST1^-T and GmMFT^-A high-oil allele type (PFAT-1) was 22.89%, which increased by about 4.5% compared with that carrying GmC2-18.1^-A, GmSWEET39^-Insertion, GmST1^-C and GmMFT^-C low-oil allele type (PFAT-14). 5 percentage points, the contribution rate to oil content is about 21.69%. 【Conclusion】Based on the markers developed above, 115 PFAT-1 high oil alleles were screened.

【Objective】 Rice blast is one of the main diseases that threaten rice yield and quality. OsBSK1-2 has been found to be involved in rice blast resistance regulation. Previously, OsCBSX4, a cysteine sulfide-β-Synthase, was identified by screening the proteins that interact with OsBSK1-2. The purpose of this study was to verify the interaction between OsBSK1-2 and OsCBSX4, and clarify the function and molecular mechanism of OsCBSX4 in rice blast resistance, providing a theoretical basis for rice disease resistance breeding.【Method】 Co-immunoprecipitation, bimolecular fluorescence complementation and luciferase complementation assays were used to determine the interaction between OsBSK1-2 and OsCBSX4. Then, quantitative PCR and agrobacterium- mediated transient transformation in N. benthamiana were used to detect the gene expression pattern and protein localization of OsCBSX4. Subsequently, OsCBSX4-knockingout and OsCBSX4-overexpressing plants were generated via CRISPR/Cas9 technology and Agrobacterium-mediated genetic transformation, respectively, and their resistance to rice blast was determined by inoculation with M. oryzae. Moreover, the immune responses induced by chitin and M. oryzae in the oscbsx4 mutant were analyzed using ROS burst and DAB staining assays respectively. In the end, the interaction between OsCBSX4 and OsRbohB was verified via bimolecular fluorescence complementation and split-luciferase complementation imaging assays, and the impact of the metabolite of OsCBSX4 on rice blast resistance was determined using the detached leaf inoculation method. The studies above will reveal the immune function and molecular mechanism of OsCBSX4.【Result】 The interaction between OsCBSX4 and OsBSK1-2 was verified by CoIP, LCI and BiFC assays. Compared to the wild type, the oscbsx4 mutant showed more disease lesions after inoculation with M. oryzae by the spray method, as well as a larger lesion area and greater fungal growth after inoculation with M. oryzae by the punch method, suggesting that knocking out OsCBSX4 decreased rice blast resistance. Moreover, the expression of the pathogenesis-related genes, OsPR5 and OsPR10, and the H₂O₂ accumulation induced by M. oryzae infection and the ROS burst induced by chitin treatment were reduced in the oscbsx4 mutant. Compared with the wild type, OsCBSX4 overexpressing plants showed a smaller lesion area and less fungal growth after inoculation with M. oryzae by the punch method, suggesting that overexpression of OsCBSX4 increased rice blast resistance. In addition, we found that OsCBSX4 can interact with OsRbohB, a key regulator of rice ROS production and treating rice with L-cysteine, a metabolite of OsCBSX4, does not affect rice blast resistance. 【Conclusion】 OsCBSX4 is an important component of OsBSK1-2 signaling and positively regulates rice blast resistance. OsCBSX4 may mediate ROS production by interacting with OsRbohB, thereby regulating rice immunity.

【Objective】 Jasmonic acid (JA) is one of key regulatory hormones during potato tuber development. Study on JA regulating tuber development mechanism would provide an important theoretical basis for tuber yield and quality formation.【Method】 In this study, potato stolons cultured in vitro were exogenously applied with JA (0, 0.5, 5, 50 μmol·L^-1) to analyze the phenotype, tissue microstructure, carbohydrate accumulation and proteome change of tubers.【Result】 With the increase of JA concentration, the tuber number formed per stolon, tuber diameter, fresh and dry weight, cell area of perimedullary zone, starch and soluble sugar content were gradually increased at 0.5 and 5 μmol·L^-1 JA (P<0.05), whereas the tuber diameter, fresh and dry weight, starch content were significantly decreased at 50 μmol·L^-1 JA (P<0.05). Lipoxygenase activity was gradually decreased with JA concentration increasing (P<0.05). Thirty-five differentially abundant proteins (P<0.05 and fold change≥2.5) closely related to tuber development regulating by JA were identified by 2-DE and MALDI-TOF/TOF-MS. They were mainly involved in bioenergy and metabolism (28.6%), cell defense and rescue (28.6%), protein biogenesis and storage (11.4%), signaling (8.6%), transcription (8.6%), unknown (8.6%) and miscellaneous (5.6%). The differential expression patterns of these proteins were clustered into three categories by hierarchical clustering analysis. Cluster 1 included 17 proteins that down-regulated at 0.5 μmol·L^-1 JA and up-regulated at 5 μmol·L^-1 JA, which mainly involved in bioenergy and metabolism, protein biogenesis and storage, signaling and transcription. Cluster 2 included 10 proteins that gradually up-regulated with JA concentration increasing, which were mainly involved in cell defense and rescue, bioenergy and metabolism, transcription. Cluster 3 included 8 proteins that down-regulated after JA treatments, which were mainly involved in bioenergy and metabolism, cell defense and rescue, protein biogenesis and storage. 【Conclusion】 Low JA concentration (0.5, 5 μmol·L^-1) might promote tuber morphogenesis by inducing cell enlargement in tuber perimedullary zone, intracellular sucrose and polysaccharide accumulation, and cell defense ability, whereas high JA concentration (50 μmol·L^-1) showed inhibitory effects.

【Objective】 Grain size and weight are the important factors affecting the yield of maize. The EIN3/EIL gene family is a sort of key transcription factors in the ethylene signaling transduction pathway, and the functions and regulatory mechanisms of the EIN3/EIL gene ZmEIL9 were analyzed in maize kernel development to elucidate its molecular mechanisms.【Method】 The expression patterns of ZmEIL9 in maize kernel at different developmental stages were analyzed by bioinformatics and RT-qPCR. The multiple sequence alignment of ZmEIL9 and its homologs from different species was performed, and the phylogenetic trees was constructed based on the neighbor-joining method. The sequence characteristics of ZmEIL9 protein were analyzed, and subcellular localization of ZmEIL9 was performed. The insertion mutants of Mu transposon and CRISPR/Cas9 knockout mutants of ZmEIL9 were screened, and the agronomic traits including grain filling rate, storage substances such as starch granule and protein content were analyzed. 【Result】 According to the members of EIN3/EIL family in maize, phylogenetic trees showed that ZmEIL9 was closely related to ZmEIL1 and SbEIL1. In the transcriptomic database of maize inbred line B73, the expression levels of ZmEIL9 were higher in the grain at early and late developmental stages. However, the expression levels were higher in inbred line N04 at the middle and late developmental stages. ZmEIL9 encoded 644 amino acids in the inbred lines Dan232 and N04, while its homolog in inbred line B73 has 642 amino acids. Subcellular localization analysis indicated that ZmEIL9 was localized in nucleus. The ZmEIL9 mutants with different Mu transposon insertion sites and CRISPR/Cas9 knockout mutants with amino acid frameshift mutations were obtained, respectively. The plant height, grain length, and 100-grain weight of Mu mutants and knockout mutants were significantly lower than those of its wild counterpart. The grain dry weights at different developmental stages were also analyzed, and the grain filling rates of Zmeil9 mutant were lower than those of the wild type. The starch granules of Zmeil9 mutant were significantly smaller and had an irregular shape based on scanning electron microscopy (SEM) observations. The contents of total starch and the concentration of zein protein in the Zmeil9 mutant were significantly lower than those in the control. 【Conclusion】ZmEIL9 plays an important regulatory role in the kernel development of maize.

【Objective】Based on the long-term experiment in the North China Plain (NCP), the differences in soil nutrient and aggregate nutrient distribution between diversified crops and wheat-maize rotation systems were investigated. Additionally, it provided a comprehensive evaluation of soil quality indices (SQI), offering a scientific basis for enhancing soil quality and productivity in the NCP. 【Method】Four diversified crop rotation systems were evaluated, including spring sweet potato-winter wheat-summer maize (Psw-WM), spring peanut-winter wheat-summer maize (Pns-WM), spring sorghum-winter wheat-summer maize (Ps-WM), with winter wheat-summer maize (WM-WM) serving as the control. The soil samples from the 0-40 cm depth were collected during the second rotation in 2022, at the flowering and harvesting stages of winter wheat. The soil enzymes activities, aggregate stability, organic matter, and concentrations of nitrogen, phosphorus, and potassium in soil and aggregates of different sizes (>2.00 mm, 0.50-2.00 mm, 0.25-0.50 mm, and <0.25 mm) were assessed. The SQI for each crop rotation system was then comprehensively evaluated. 【Result】Compared with WM-WM, the three other crop rotations increased soil inorganic nitrogen content. Psw-WM significantly enhanced organic matter in the 0-20 cm layer, total nitrogen in soil aggregates (>2.00 mm, 0-10 cm), and organic matter in soil aggregates (>2.00 mm and 0.50-2.00 mm, 0-10 cm), which also increased cellulase, catalase, and alkaline protease activities. Pns-WM improved organic matter in the 20-40 cm layer and available potassium in soil aggregates (0.25-0.50 mm and >2.00 mm, 10-20 cm), as well as organic matter in soil aggregates (0-10 cm, >2.00 mm and 10-20 cm, >0.50 mm), which also increased sucrase, urease, and alkaline protease activities. Psw-WM improved the stability of 0-10 cm soil aggregates, while Pns-WM improved the stability of 0-30 cm soil aggregates. Both Pns-WM and Psw-WM significantly improved the SQI, with Pns-WM showing a higher improvement than Psw-WM. The path analysis revealed that the average weight diameter (MWD) of aggregates was a direct and significant affecting SQI. It also had a significant indirect positive effect on SQI by influencing inorganic nitrogen. Additionally, the increased organic matter led to a higher proportion of large aggregates, which significantly affected SQI indirectly. 【Conclusion】Legume (peanut) and root crop (sweet potato) rotations with wheat-maize rotations could significantly improve soil quality and enhance the soil nutrient supply capacity in the NCP.