Loading...

Table of Content

    16 July 2021, Volume 54 Issue 14
    CROP GENETICS & BREEDING·GERMPLASM RESOURCES·MOLECULAR GENETICS
    Production of High Linoleic Acid Rice by Genome Editing
    CHEN XiaoJun,WU KaiShen,HUI Jian,BAI HaiBo,MA SiShuang,LI JingYu
    Scientia Agricultura Sinica. 2021, 54(14):  2931-2940.  doi:10.3864/j.issn.0578-1752.2021.14.001
    Abstract ( 476 )   HTML ( 96 )   PDF (2753KB) ( 386 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】 Linoleic acid is one of essential fatty acid for human, which can accelerate the metabolism and decomposition of fat in the body, reduce the accumulation of cholesterol on the wall of blood vessels, effectively prevent the occurrence of arteriosclerosis, improve human immunity, promote bone development, improve memory and prevent brain function degradation. Rice bran is a by-product of rice processing, and is an edible oil with high nutritional value. In this paper, we precisely edited the gene of ω-3 fatty acid desaturase, which is the key enzyme controlling fatty acid synthesis pathway. It can enrich the upstream linoleic acid, while closing the synthesis pathway of downstream products. 【Method】There are two copies of the key enzyme ω-3 fatty acid desaturase gene ( OsFAD3) in the fatty acid synthesis pathway in rice, which are located on chromosome 11 and chromosome 12, and their cDNA homology is 97.32%. According to the basic principles of gene editing, the specificity of edit depends on the characteristics of the guide RNA (gRNA). In this study, two gRNAs were designed and synthesized in the exon region of two homologous genes, and were constructed into editing vectors respectively. The local receptor material (Fuyuan 4) was successfully transformed by Agrobacterium mediated method. The editing efficiency and genotype of the two sites were analyzed. At the same time, the main agronomic characters of OsFAD3 double mutant were measured, such as grain width and grain length. The contents of 37 fatty acids were figured by GC-MS in the grains of OsFAD3double mutants. 【Result】The results showed two homozygous editing materials were obtained, and other editing materials with different genotypes were obtained. Compared with the control materials, the main agronomic traits such as grain width and grain length of OsFAD3 double mutant materials had no significant changes, however, the relative content of linoleic acid in rice increased by 3.36%. 【Conclusion】One gRNA vector knock out two ω-3 fatty acid desaturase genes at the same time, which improves the relative content of linoleic acid in grains without changing the main agronomic traits of seeds.

    Development and Effect Evaluation of KASP Markers Closely Linked to Major QTLs of Spike Number Per Unit Area and Grain Length in Wheat
    FAN Tao,LI Zhi,JIANG Qing,CHEN ShuLin,OU Xia,CHEN YongYan,REN TianHeng
    Scientia Agricultura Sinica. 2021, 54(14):  2941-2951.  doi:10.3864/j.issn.0578-1752.2021.14.002
    Abstract ( 465 )   HTML ( 48 )   PDF (784KB) ( 285 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Spike numbers per unit area (SN) and kernel length (KL) are both important yield-related traits of wheat. Genetic improvement on SN and KL will help increase the yield of wheat. The corresponding KASP molecular markers were developed for the QSn.sau-2D.2 and the QTL Kl.sau-3D.2, which were identified in previous study for increasing SN and KL, respectively. Both KASP markers were verified and evaluated in the RIL population crossed by Chuannong 18 and T1208, which laid a foundation for better utilization of these two QTLs.【Method】The major QTL QSn.sau-2D.2 for SN and the major QTL QKl.sau-3D.2 for KL were previously identified in the RIL population crossed by Chuannong 18 and T1208. The SNP molecular marker sequence of Wheat 55K SNP array within the two QTLs were used to develop and design the KASP molecular markers. The polymorphic KASP markers were screened with wheat parents and then verified in the RIL population. The selected KASP molecular markers were genotyped and used to identify the corresponding phenotypic traits in Chuannong 18×T1208 RILs population. The effects of the two major QTLs on other agronomic traits were analyzed.【Result】KASP-AX-111151907 and KASP-AX-109962767 have polymorphic between the wheat parents. Subsequently, these two KASP markers were verified in the RIL population linkage with QSn.sau-2D.2 and QKl.sau-3D.2, respectively. KASP-AX-111151907 and KASP-AX-10996276 could divide the genotypes into two groups. According to the phenotype, the average selection rate of KASP-AX-111151907 for multi spike lines was 72.58%, the average selection rate of KASP-AX-111151907 for few spike lines was 71.68%, the average selection rate of KASP-AX -10996276 for long kernel lines was 69.86%, and the average selection rate of KASP-AX-10996276 for short kernel lines was 61.52%, indicating the reliability of the two KASP markers. Moreover, the genotyping results based on KASP molecular markers showed that the two QTLs had significant effects on plant height, 1000 kernel weight, kernel length, kernel width, kernel diameter ratio, spike numbers per unit area and kernel weight per spike, respectively. The validation in the RIL population of Chuannong17×Chuannong11 also indicated that these two KASP markers had a certain effect on the selection of the corresponding traits.【Conclusion】This study developed two KASP markers for two major QTLQSn.sau-2D.2 and QKl.sau-3D.2, respectively. Both KASP markers can be used for the selection of corresponding traits. These two QTLs which are tightly linked to these two KASP markers could significantly improve the SN and KL, respectively. In addition, QSn.sau-2D.2had negative effects on plant height, 1000 kernel weight, kernel length, kernel width, kernel diameter ratio and kernel weight per spike. On the other hand,QKl.sau-3D.2 had positive effect on plant height, 1000 kernel weight, kernel width, kernel diameter ratio and kernel weight per spike, but had negative effects on SN. These two QTLs and the developed KASP markers can be used in the future high yield wheat breeding program.

    Cytological Characters Analysis and Low-Resolution Mapping of Stripe-Leaf MutantA36-S in Foxtail Millet
    ZHANG Shuo,ZHI Hui,TANG ChanJuan,LUO MingZhao,TANG Sha,JIA GuanQing,JIA YanChao,DIAO XianMin
    Scientia Agricultura Sinica. 2021, 54(14):  2952-2964.  doi:10.3864/j.issn.0578-1752.2021.14.003
    Abstract ( 334 )   HTML ( 38 )   PDF (4565KB) ( 294 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】 Foxtail millet is a C4 model plant, and its leaf color mutants are good materials for C4 photosynthetic pathway research. Through the cytological characters analysis and gene initial mapping of the stripe-leaf mutant A36-S in foxtail millet, it laid the foundation for cloning the mutant gene, analyzing the chloroplast biogenesis and development, and further understanding the C4 photosynthetic regulation mechanism in foxtail millet. 【Method】 The stripe-leaf mutant of foxtail millet A36-S was naturally mutated from intermediate material A36 created by breeding. Comparing the phenotypic characteristics of A36-S and its isogenic line A36-N, which showed normal phenotypes, and investigating the agronomic traits, such as plant height, leaf width, leaf length, panicle weight, thousand-grain weight, and seed setting rate. To analyze the photosynthetic characters of A36-S, the chlorophyll content, net photosynthetic rate, intercellular CO2 concentration, stomatal conductance, and transpiration rate of A36-S and A36-Nwere determined. By observing the leaf transverse section and ultrathin section of A36-S and the control variety Yugu1, the leaf anatomical structure characters were analyzed, by counting the numbers and areas of the chloroplasts in mesophyll cells and bundle sheath cells respectively, the chloroplast biogenesis and development were assessed. An F2 segregation population of A36-S×SSR41 were created, and genetic analysis was conducted by counting the number of normal phenotype single plant and stripe-leaf single plant in the population. The DNA mixed pools of normal single plants and stripe-leaf single plants of the F2 segregation population were constructed separately, and the method of Bulked Segregation Analysis (BSA) was used to locate the mutant gene. By screening the stripe-leaf plants in the F2 generation using SSR and In-Del markers, the mutant gene were furtherly located.【Result】The stripe-leaf mutant of foxtail millet A36-S showed the phenotype of irregular white stripe-leaf in the whole growth period. Agronomic traits analysis showed that compared with its isogenic line A36-N, A36-S decreased significantly in plant height, leaf width, panicle weight, thousand-grain weight, and setting percentage. Photosynthetic index measurement showed that the chlorophyll contents of A36-S were also reduced significantly, especially the chlorophyll b content declined more severely, additionally, the net photosynthetic rate was also decreased significantly. Observation of the leaf anatomical structure showed that the chloroplasts number and area were significantly lower than that of the contrast Yugu1, while the changes in Kranz structure were not obvious. Furtherly, the ultrastructure of chloroplast was observed and showed that the chloroplast development situation in different cells was quite different. And the leaf cells ofA36-S could be classified into three types: type Ⅰ cells had normal chloroplasts, type Ⅱ cells had chloroplasts with reduced grana and lamellar structures, while type Ⅲ cells had severely abnormal chloroplasts or even had no chloroplast. Genetic analysis suggested that the stripe-leaf trait of A36-S was controlled by a single recessive gene, and the mutant gene was located to the region from 7.66 Mb to 27.90 Mb of chromosome 4 by F2 population. 【Conclusion】 Stripe-leaf mutant of foxtail millet A36-S represented decreased agronomic traits and photosynthetic capacity, and the number, size, and ultrastructure of leaf cell chloroplast were significantly abnormal. The stripe-leaf trait of A36-S was controlled by a single recessive gene, which was mapped to a region from 7.66 Mb to 27.90 Mb of chromosome 4.

    TILLAGE & CULTIVATION·PHYSIOLOGY & BIOCHEMISTRY·AGRICULTURE INFORMATION TECHNOLOGY
    Assessment of Terrestrial Laser Scanning and Hyperspectral Remote Sensing for the Estimation of Rice Grain Yield
    LI PengLei, ZHANG Xiao, WANG WenHui, ZHENG HengBiao, YAO Xia, ZHU Yan, CAO WeiXing, CHENG Tao
    Scientia Agricultura Sinica. 2021, 54(14):  2965-2976.  doi:10.3864/j.issn.0578-1752.2021.14.004
    Abstract ( 551 )   HTML ( 7 )   PDF (2500KB) ( 234 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Background】Non-destructive and accurate estimation of crop biomass and yield is crucial for the quantitative diagnosis of growth status and planning of national food policies. Hyperspectral and Terrestrial Laser Scanning (TLS) remote sensing provide convenient and effective ways to monitor crop condition.【Objective】The aim of this study was to examine the feasibility of developing models with various independent datasets to build a universal yield monitoring model. The expected results can provide theoretical and technical support for rice growth monitoring and scientific guidance of precision agriculture.【Method】Field plot experiments were conducted in 2016, 2017 and 2018 and involved different study sites, nitrogen (N) rates, planting techniques and rice varieties. Linear regression (LR) and random forest (RF) were evaluated in estimating yield with TLS and spectral data collected since the heading stage, and the feasibility of developing models with various independent datasets was examined to build a universal yield monitoring model.【Result】 The results showed that TLS models exhibited higher estimation accuracies for yield in the heading stage, early-filling stage and late-filling stage (R2: 0.64-0.69) than hyperspectral models (R2: 0.20-0.58). Compared to RF, LR modeling yielded significantly higher validation accuracies for growth stages after heading. While the predictive model was transferred to other datasets, the validation accuracies from the same site (RRMSE: 11.37%-12.41%) were higher than those from a different site (RRMSE: 16.69%-17.85%).【Conclusion】The results suggested that TLS was a promising technique to monitor yield at post-heading stages with high accuracy and to overcome the saturation of canopy reflectance signals encountered in optical remote sensing.

    Effects of Tillage on Soil Moisture and Yield of Wheat-Maize Rotation Field in Weibei Upland Plateau
    GUO XingYu, WANG Hao, YU Qi, WANG Rui, WANG XiaoLi, LI Jun
    Scientia Agricultura Sinica. 2021, 54(14):  2977-2990.  doi:10.3864/j.issn.0578-1752.2021.14.005
    Abstract ( 403 )   HTML ( 6 )   PDF (588KB) ( 151 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】 In this paper, the effects of tillage practices on soil moisture and crop yield in winter wheat-fallow-spring maize rotation fields under different rainfall patterns were studied, so as to provide the theoretical support for the efficient use of rainfall and innovation of tillage systems in dryland region.【Method】 In 2007-2019, a long-term consecutive in-situ conservation tillage experiment was conducted in the Weibei Upland Plateau. Two conservation tillage practices were set, including no-tillage (NT) and subsoiling (ST), using conventional tillage (CT) as a control. Then, the effects of different tillage methods on soil moisture storage, yield and water use efficiency during the fallow and growth periods under different rainfall patterns were analyzed. 【Result】(1) The rainfall pattern and tillage practice as well as their interactions all significantly affected the soil water storage at the end of the fallow period and moisture storage rate, with rainfall pattern as the dominant factor. The soil water storage (mm) in the 0-200 cm soil layer at the end of fallow period values decreased in the following order: humid year (430.6) > dry year (405.9) > normal year (381.5); NT (417.4) > ST (402.3) > CT (398.2); The soil water storage in the 0-200 cm soil layer at the end of fallow period was the highest in the humid year under NT treatment (438.5) and the lowest in the normal year under ST treatment (370.2). The moisture storage rate (%) during the fallow period values decreased in the following order: humid year (27.1) > dry year (26.6) > normal year (25.1); NT (27.6) > ST (26.4) > CT (25.8); NT in the dry year had the highest moisture storage rate (29.1) and CT in the normal year had the lowest moisture storage rate (25.0). (2) The rainfall pattern and tillage practice as well as their interactions all significantly affected yield and water use efficiency (WUE) in winter wheat, of which the influence of tillage practice was the most significant factor of winter wheat yield (kg∙hm-2) and WUE (kg∙hm-2∙mm-1). The winter wheat yield values decreased in the following order: humid year (4985) > dry year (3984); NT (4522) > ST (4468) > CT (4465); The winter wheat yield was the highest in the humid year under NT treatment (5033) and the lowest in the dry year under ST treatment (3957). The winter wheat WUE values decreased in the following order: humid year (15.4) > dry year (14.9); NT (16.2) > ST (15.4) > CT (14.0); The winter wheat WUE was the highest in the humid year under NT treatment (16.5) and the lowest in the dry year under CT treatment (13.9). (3) The rainfall pattern and tillage practice as well as their interactions all significantly affected yield and WUE in spring maize, of which the influence of rainfall pattern was the most significant factor of spring maize yield and WUE. The spring maize yield values (kg∙hm-2) decreased in the following order: humid year (7677) > dry year (6999) > normal year (5887); NT (6900) > ST (6879) > CT (6785); The spring maize yield was the highest in the humid year under ST treatment (8003) and the lowest in the normal year under ST treatment (5723). The spring maize WUE values (kg∙hm-2∙mm-1) decreased in the following order: humid year (18.7) > dry year (17.5) > normal year (14.8); NT (17.8) > ST (17.0) > CT (16.2); The spring maize WUE was the highest in the humid year under NT treatment (19.3) and the lowest in the normal year under CT treatment (13.9). (4) In the winter wheat-fallow-spring maize rotation cycle, the tillage practice significantly affected the crop yield and WUE. The mean yield and WUE under NT and ST were 1.6%, 1.2% and 9.5%, 3.9% higher than that under CT treatment, respectively. 【Conclusion】 Integrating soil moisture and crop yield under each rainfall pattern, conservation tillage could achieve water conservation and yield increase relatively efficiently, among which the effects of water storage and yield of no-tillage were better. Based on the high yield, high efficiency and long-term sustainable development of upland plateau, no-tillage was recommended to realize the aim of water storage, soil water conservation, yield and income increased in Loess Plateau region.

    Effect of the Mechanical Harvesting Methods on the Key Traits of Rapeseed
    BAI ChenYang, HE HanZi, JIA CaiHua, LI XiaoHua, REN YiLin, YE Jun, WANG Bo, KUAI jie, ZHOU GuangSheng
    Scientia Agricultura Sinica. 2021, 54(14):  2991-3003.  doi:10.3864/j.issn.0578-1752.2021.14.006
    Abstract ( 310 )   HTML ( 1 )   PDF (607KB) ( 132 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Harvesting is one of the key links in rapeseed production, which affects the yield and quality. In this experiment, artificially simulated combined harvesting and two-stage harvesting were used and analyzed the impacts on key traits of rapeseeds. This study provides a basis for the selection of mechanical harvesting methods and the determination of supporting parameters.【Method】Two middle-ripening rapeseed varieties with different oil content were selected, and different harvesting periods were used. For combined harvesting the period is 20-44 d in Huanggang and 23-47 d in Xiangyang, and the two-stage harvesting is practiced on different time 20, 26, and 32 d in Huanggang while, 23, 29, 35 d in Xiangyang, after final flowering, whereas 3, 6, 9, and 12 d post-ripening, in Huanggang and Xiangyang, Hubei. Furthermore, key indicators such as1000-grain weight, water content, quality, etc were measured and analyzed. 【Result】The two pilots and two varieties of the combined harvesting reached the physiological maturity about 38 days after final flowering, the required accumulated temperature was about 1 100℃, and the thousand-grain weight and oil content were the largest. Approximately 46 days after the final flowering, the seed moisture and chlorophyll content rapidly decreased to the optimum point, which is a suitable period for mechanized combined harvesting. Moreover, reduction in seed moisture content is highly significantly positively correlated with daily average temperature and atmospheric pressure, while decreased chlorophyll content is highly significantly positively correlated with daily average temperature, but significantly negatively correlated with rainfall. Under D3 (about 33 d after final flowering), the highest average content of 1000-grain weight, oil contents, rate of oil yield, uniformity and oleic acid, and the lowest average content of chlorophyll, linoleic acid and moisture under the two-stage harvest while comparing D3 with D1 and D2, Our study showed that the number of days of post-ripening has no significant effect on the oil content, oil yield, thousand-grain weight, uniformity, oleic acid, and linoleic acid content, but the moisture content rapidly decreased to the optimum point after 6-9 days, which is a suitable time for picking and threshing. The seed moisture content decreased and significantly positively correlated with solar radiation but highly significantly negatively correlated with rainfall. 【Conclusion】During the period of final flowering to the harvesting, the average daily temperature is 19-22℃, the rainfall is 1.0-3.5 mm, and the atmospheric-pressure is 1.45-1.75 kPa. However, under prevailing climatic conditions, the combined harvesting is done on approximately 46th day after final flowering (which is about 8 days after physiological maturity) and the two-staged harvesting is used to cut down on about 33rd day after the final flowering (which is about 5 days before the physiological maturity), then post-ripening practice is carried out (average relative humidity is 64%-80%, the rainfall is 0.7-3.1 mm, and the solar radiation is 10-13 MJ·m-2), to picked the harvested crop after 6-9 days of two-stage harvesting. Under this condition, the moisture content of seed can be reduced to about 15%, and the chlorophyll content can be reduced to below 10 mg·kg-1. There was no significant difference between the two harvesting methods in terms of 1000-grain weight, oil content, rate of oil yield and seed uniformity, however, the oleic acid and linoleic acid contents of seed were higher, seed and oil quality were better in two-stage harvesting as compared with the combined harvesting.

    PLANT PROTECTION
    Resistance Monitoring to Thiamethoxam and Expression Analysis of Cytochrome P450 Genes in Leptinotarsa decemlineata from Xinjiang
    SHI Xin,LI Sha,WANG ZhiMin,FU KaiYun,FU WenJun,JIANG WeiHua
    Scientia Agricultura Sinica. 2021, 54(14):  3004-3016.  doi:10.3864/j.issn.0578-1752.2021.14.007
    Abstract ( 342 )   HTML ( 26 )   PDF (615KB) ( 275 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Colorado potato beetle (CPB), Leptinotarsa decemlineata, is a significant agricultural pest in potato worldwide. The objective of the study is to identify the main cytochrome P450 genes involved in metabolic regulation of neonicotinoid insecticide thiamethoxam in L. decemlineata.【Method】In this study, the resistance level to thiamethoxam was assayed using topical applications in 11 field populations of L. decemlineata from Qapqal (QPQL), Yining (YN), Tacheng (TC), Urumqi (URMQ), Jimusar (JMS) of Xinjiang in 2018 and 2019. The activities of three detoxifying enzymes including cytochrome P450 monooxygenase (P450), glutathioneS-transferase (GST) and esterase (EST), exposed to LD50of thiamethoxam for 72 h were analyzed by in vitro enzyme activity assay. The differentially expressed genes (DEGs) of L. decemlineata adults susceptible and resistant to thiamethoxam were detected by Illumina HiSeqTM 2500 sequencing platform. The expression verification of three P450 genes and the expression change of six P450 genes in different populations and the beetles exposed to LD50 of thiamethoxam for 72 h were performed by real-time quantitative PCR (qRT-PCR).【Result】Resistance monitoring results showed that the adults of YN1, URMQ1, TC and QPQL1 populations in 2018 and YN2 and JMS in 2019 developed low to moderate level resistance to thiamethoxam with resistance ratio (RR) ranging from 5.99 to 20.33 times. The P450 activities of adults from two sensitive populations URMQ2 and URMQ3 and two resistant populations YN2 and JMS after treatment with thiamethoxam LD50for 72 h were significantly increased to 1.76, 2.75, 1.91 and 1.66 times as high as the control, respectively. In addition, the GST (CDNB and DCNB as substrates) of URMQ3 population and the EST activities of URMQ2 and YN2 populations were obviously increased to 1.19, 2.10 and 1.35, 1.91 times of the control, respectively. Based on Illumina RNA sequencing, the raw reads (56 872 051 and 62 249 136), clean reads (55 903 706 and 61 082 076), clean bases (8.39 and 9.16 G) and error of basic group (0.03%) were obtained for the thiamethoxam-susceptible and -resistant samples of adult populations, respectively. Thirteen differentially expressed P450 genes were identified in the thiamethoxam-resistant sample, of which two genes were significantly up-regulated. The expression levels of three P450 genes, namely CYP4C1,CYP9e2and CYP305a1by RNA sequencing were validated by qRT-PCR analysis. Furthermore, CYP9Ya, CYP9Yc, CYP4C1, CYP305a1 and CYP9Z25 were up-regulated in the thiamethoxam-resistant adults, and the expression of CYP9Ya, CYP9Yc increased significantly in the 4th instar larvae and adults under thiamethoxam treatment. It was found by correlation analysis that there was a significant positive correlation between resistance level to thiamethoxam and CYP9Ya expression of adults. 【Conclusion】CYP9Ya may be involved in the production of thiamethoxam resistance in L. decemlineata, and the role of other genes cannot be ruled out.

    Expression, Purification and Functional Analysis of Odorant Binding Protein 11 (OBP11) in Anomala corpulenta
    QIN JianHui,LI JinQiao,ZHAO Xu,LI KeBin,CAO YaZhong,YIN Jiao
    Scientia Agricultura Sinica. 2021, 54(14):  3017-3028.  doi:10.3864/j.issn.0578-1752.2021.14.008
    Abstract ( 300 )   HTML ( 35 )   PDF (4865KB) ( 281 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】The objective of this research is to explore the function of Anomala corpulenta odorant binding protein 11 (AcorOBP11) by studying the binding characteristics of AcorOBP11 with the host plant volatiles, and to lay the foundation for clarifying the olfactory molecular mechanism of AcorOBP11.【Method】Reverse transcription PCR (RT-PCR) was used to clone the full-length ORF ofAcorOBP11 by specific primers, and sequences with high similarity were downloaded for sequence comparison and analysis by BLAST. Prokaryotic expression of AcorOBP11 was conducted by Escherichia coliprotein expression system, the OBP11 was inserted into the expression vector pET28a, and the recombinant plasmid was transferred into E. coli competent cells BL21 (DE3). The target protein was purified by using recombinant enterokinase to remove the His-tag and purified again by using a nickel column. The purified protein was diluted with 50 mmol·L-1 Tris-HCl (pH 7.4) to a dilution of 2 μmol·L-1, and the odorant was diluted with chromatographic grade methanol to a final concentration of 1 mmol·L-1, and 1-NPN was used as a fluorescent probe to determine the binding characteristics of AcorOBP11 to 37 host plant volatiles by fluorescence competitive binding assay. Modeller was used to obtain the three-dimensional structure of AcorOBP11 by using AgamOBP48 (PDB ID: 4ij7) as a template. Autodock semi-flexible docking was used to simulate the binding of AcorOBP11 to host plant volatiles.【Result】The full-length ORF ofAcorOBP11 was amplified, which is 639 bp in total, including 17 amino acids of signal peptide at the N-terminal. With eight conserved cysteine sites, AcorOBP11 belongs to the Plus-C OBP subfamily. The evolutionary tree results indicated that it had the closest evolutionarily relationship with HparOBP9. AcorOBP11 was successfully inserted into pET-28a and expressed at 0.25 mmol·L-1 IPTG and induced expression at 37℃. The target protein was obtained by nickel column purification twice. The results of competitive binding experiments showed that AcorOBP11 has a good binding ability to 13 kinds of host plant volatiles, such as α-ionone, 3-methylbutanal, 2-ethylhexyl acrylate, bornyl acetate, cinene, trans-2-hexenal, 2-ethylhexanal, isovaleric acid. Among them, α-ionone has the best binding ability, and its competitive dissociation constant is 22.78 μmol·L-1. The molecular docking results showed that AcorOBP11 has the lowest binding free energy with α-ionone, which is -24.74 kJ·mol-1, and forms a hydrogen bond at Cys163, indicating that it has the strongest binding ability with α-ionone, which is consistent with the fluorescence competition binding result. 【Conclusion】AcorOBP11 can recognize a variety of host plant volatiles, and it is speculated that it plays an important role in locating host plants for A. corpulenta. The results will provide new insights forA. corpulenta ecological control.

    Knockout of Single Allele of fl(2)d Significantly Decreases the Fecundity and Fertility inPlutella xylostella
    LI FeiFei,WANG BeiBei,LAI YingFang,YANG FeiYing,YOU MinSheng,HE WeiYi
    Scientia Agricultura Sinica. 2021, 54(14):  3029-3042.  doi:10.3864/j.issn.0578-1752.2021.14.009
    Abstract ( 294 )   HTML ( 33 )   PDF (2751KB) ( 235 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】 RNA methylation is the main form of epigenetic modification at post-transcriptional level, which is involved in many important cellular processes. The diamondback moth, Plutella xylostella, is an important oligophagous insect pest, causing serious loss on the production of cruciferous vegetables. However, the function of RNA methylation-related genes in P. xylostella is still unclear. The present study aims to identify and clone the homologous fl(2)d, one of the members of the RNA methylation protein complex (writers), to determine the expression pattern offl(2)d, and to knockout fl(2)d using CRISPR/Cas9 for the investigation of its biological functions in P. xylostella.【Method】The sequence of homologous fl(2)d was identified in the genome database ofP. xylostella, which was used for PCR amplification of the coding sequence (CDS). Quantitative real-time PCR (qRT-PCR) was used to study the relative expression levels of fl(2)d in different developmental stages and adult gonads of P. xylostella. The fl(2)d was edited using CRISPR/Cas9 combined with egg injection. Each of the adults that developed from the injected eggs was used to pair with a wild-type adult for reproduction. Offspring of the same population was forced to inbreed by single-pair mating to establish the mutant strains. The differences of genetic characters, biological parameters and phenotypes between mutants and wild-type individuals were recorded and compared to decipher the function offl(2)d.【Result】The CDS of fl(2)d with length of 912 bp was isolated, the expression of which was high in female pupa, adult and egg, moderate in male adult and pupa, the lowest in larva, and significantly higher in ovary than in testis of adult. The sgRNAs targeting fl(2)d and the Cas9 protein were mixed to inject eggs, and the offsprings carrying mutant alleles were screened for homozygous strains based on single-pair inbreeding for 10 generations. Three types of heterozygous mutant strains both predicted to cause frameshift of the CDS were obtained, with the deletion of 4 (Δfl(2)d213-4), 5 (Δfl(2)d213-5) and 7 (Δfl(2)d213-7) bases. During the screening process, six and two homozygous mutants from Δ fl(2)d213-4 and Δ fl(2)d213-5 strains were identified, respectively. The homozygous mutants of Δ fl(2)d213-4 successfully mated in two pairs, but no eggs were produced. Meanwhile, each two male adults of homozygous mutants of either Δ fl(2)d213-4 or Δ fl(2)d213-5 were mated with the same type of female heterozygous mutant, and also no eggs were produced. The results indicated that individuals with homozygousfl(2)d mutation may have extremely low survival rate and not be able to produce offspring. Through analyzing separation ratio of the genotypes of offspring from the inbreeding of heterozygous mutants and the hybridization between heterozygous mutants and wild-type, it was found that the ratio of heterozygous mutant individuals to wild-type was slightly less than 2 and 1, respectively, indicating that heterozygous mutation of fl(2)d would affect the normal growth and development of P. xylostella, and in some cases would lead to death. The offsprings of mutant individuals, which carry a mutant allele, showed a sex ratio close to 1﹕1 (P<0.05). It was speculated that thefl(2)d might not be involved in sex determination in P. xylostella. For the mating consists of mutant adults, the fecundity and hatchability were significantly lower (P<0.01) than the mating between wild-type adults. Most of the eggs produced from the mutant parents look abnormal, and could not hatch normally due to water loss and shrinkage. Based on the dissection of adult gonads, it was found that the number of attached eggs on the ovary of the mutant female adult and the wild-type female adult that has mated with mutant male adult was less than that of the wild-type virgin female adult, while no obvious abnormality was found for the testis of mutant male adult. Some of the hatched heterozygous mutants showed different degrees of distortion during the whole developmental process, resulting in the failure to complete life cycle. A small part of the heterozygous mutant individuals could develop normally, and thus transmit the mutant allele to their offspring. According to our findings, a model of genetic control ofP. xylostella based on fl(2)d was proposed.【Conclusion】The fl(2)d is involved in the reproductive process and embryonic development ofP. xylostella, mutation of which significantly affects the population size of the offspring, making it an ideal target for the genetic control of P. xylostella.

    SOIL & FERTILIZER·WATER-SAVING IRRIGATION·AGROECOLOGY & ENVIRONMENT
    Advance in Indicator Screening and Methodologies of Soil Quality Evaluation
    LI Xin,ZHANG WenJu,WU Lei,REN Yi,ZHANG JunDa,XU MingGang
    Scientia Agricultura Sinica. 2021, 54(14):  3043-3056.  doi:10.3864/j.issn.0578-1752.2021.14.010
    Abstract ( 2077 )   HTML ( 112 )   PDF (949KB) ( 1091 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】 The objectives of the present study were to synthesize the current information on soil quality assessment method and indicator system, and to present the hot topics and frontiers related to soil quality, so as to, provide references for Chinese scholars and experts in the field of soil quality evaluation research and application. 【Method】 The published articles regarding the selection of soil quality indicators, construction of minimum data set, and selection of soil quality evaluation methods were collected based on Web of Science and CNKI databases using bibliometrics method, and a total of 415 articles and 155 minimum data sets related to soil quality evaluation were screened. Development trend, frontier fields and current problems of global soil quality assessment during the past 30 years were analyzed according to selection frequency of indicators, assessment method and construction of minimum data set. 【Result】 The soil quality evaluation indicator system mainly included 25 physical, 36 chemical, 35 biological and 19 environmental indicators. Soil organic matter, as the core indicator of soil quality, was selected with the highest frequency of 96.6%, followed by pH, total nitrogen, available phosphorus, available potassium, and bulk density, with a frequency more than 50%. The selection frequency of biological indicators such as microbial biomass and soil enzyme activity was less than 25%, while increasing over time. Principal component analysis, minimizing indicator redundancy and reflecting most of the information of original variables, was the most widely used for minimum data set construction method. Soil organic matter, available phosphorus, bulk density, and pH were selected into the minimum data set with a frequency of 67.7%, 43.2%, 34.8%, and 34.2%, respectively, being widely used to characterize soil quality. Nowadays, the most studies on soil quality evaluation focus on the utilization of principal component analysis to select soil quality indicators and establish soil quality index for comprehensive soil quality evaluation, which was suitable for sustainable soil management.【Conclusion】 Soil organic matter, available phosphorus, soil pH, bulk density and soil water content were the main parameters selected for soil quality evaluation. Construction a comprehensive and objective soil quality indicator system and the integration with the information technology would be the focus in future research. The application of evaluation indicators in large-scale soil quality assessment was the trend of future development.

    Contribution of Carbon Sources in Sedimentary Soils Combining Carbon and Nitrogen Isotope with Stable Isotope Model
    LI Na,SUN ZhanXiang,ZHANG YanQing,LIU EnKe,LI FengMing,LI ChunQian,LI Fei
    Scientia Agricultura Sinica. 2021, 54(14):  3057-3064.  doi:10.3864/j.issn.0578-1752.2021.14.011
    Abstract ( 499 )   HTML ( 34 )   PDF (818KB) ( 685 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】To study the sources of deposited soil organic carbon (SOC) under different land use patterns in a typical small watershed in the brown soil hilly area of western Liaoning through eroded sedimentation, and to provide a scientific reference for the reasonable control of soil carbon loss caused by soil erosion in the small watershed. 【Method】Through field sampling of small watersheds in the hilly and gully area of western Liaoning, the sources of deposited soil carbon in the small watersheds were studied and their contribution was quantified. Using GIS and GPS technology to analyze the surface soil of 4 different land use types (cropland, forest, grassland, gully) in the small watershed and 3 locations of the check dam in the small watershed (S1 in front of the dam, S2 in the middle of the dam, S3 behind the dam) 0-100 cm soil profile was sampled to analyze the carbon source of the sedimentary soil based on a mixed carbon and nitrogen isotope model.【Result】Using13C and15N isotopic characteristics and their elemental composition qualitative and quantitative identification of soil organic carbon in eroded sediments in the hilly and gully area of western Liaoning was carried out. The SOC loss was primarily from cropland, accounting for 58.75%, followed by gully (25.49%), forest (9.2%), and grassland (6.49%). 【Conclusion】 The stable isotope SIAR mixing model, as a reliable "fingerprint" tool, could be successfully employed to estimate the contribution of various C sources within a complex ecosystem, The research results can provide theoretical references for soil protection and nutrient loss control in small watersheds affected by water erosion, and for maintaining the sustainability of the ecosystem.

    Response of Water Use and Yield of Dryland Winter Wheat to Nitrogen Application Under Different Rainfall Patterns
    LIU PengZhao,ZHOU Dong,GUO XingYu,YU Qi,ZHANG YuanHong,LI HaoYu,ZHANG Qi,WANG XuMin,WANG XiaoLi,WANG Rui,LI Jun
    Scientia Agricultura Sinica. 2021, 54(14):  3065-3076.  doi:10.3864/j.issn.0578-1752.2021.14.012
    Abstract ( 284 )   HTML ( 28 )   PDF (698KB) ( 226 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Rainfall and nitrogen application are the main factors affecting winter wheat production in Weibei dryland, especially the interannual variability of rainfall is large. Therefore, their yield-increasing effects are significantly affected by the type of rainfall year. In this study, the effects of nitrogen rates on water use, grain yield and protein content in dryland wheat field under different rainfall patterns were analyzed, which provided a theoretical basis for “nitrogen applying according to rain” and ensuring stable yield and high quality of winter wheat in Weibei Dryland.【Method】A 3-year field experiment (2017-2020) of winter wheat (Jinmai 47) was performed with different nitrogen fertilization at five levels (0, 60, 120, 180 and 240 kg·hm-2, represented as N0, N60, N120, N180, and N240, respectively) in Heyang County, located in Weibei dryland of Shaanxi, and the effects of nitrogen application under different rainfall patterns on soil water dynamics, water use efficiency (WUE), wheat yield performance and grain protein content were evaluated. 【Result】Different rainfall patterns had significant impacts on soil water storage before sowing (SWSS), soil water content during growth period, ET, WUE, yield and protein content of winter wheat. (1) There was a linear correlation between rainfall in fallow period (from July to September) and SWSS, with an increment of 0.9 mm SWSS per 1 mm rainfall. In the humid and normal years with adequate rainfall during fallow stage, the SWSs in present winter wheat growth season was not significantly influenced by the increase of nitrogen fertilization in previous growth season. However, in the dry year with less rainfall in fallow stage, the SWSS in present winter wheat growth season decreased significantly by 15.4 mm when nitrogen fertilization in previous growth season was increased by each 100 kg·hm-2. Compared with dry and normal year, the soil water content of 0-200 cm soil layer during growth period of winter wheat could be increased in humid rainfall year, thus evapotranspiration (ET) was increased by 35.7% and 6.6%, respectively. The soil water accumulation of 0-120 cm soil depth during the growth period fluctuated greatly under the influence of rainfall and the growth of winter wheat. However, the soil water content in 160-200 cm deep soil depth showed a stable change trend. Compared with dry and normal year, the WUE in humid pattern was increased by 55.7% and 26.5%, the grain yield was increased by 112.3% and 39.1%, and protein content (PC) was decreased by 8.3% and 5.2%, respectively. (2) Compared with N0 treatment, soil water content in 0-200 cm soil depth was decreased by nitrogen applied during each growth period under humid, dry and normal years. The nitrogen fertilizer application increased ET by 4.6%-14.6%, 6.0%-8.6% and 2.2%-9.5%, increased WUE by 20.7%-39.8%, 4.7%-33.3%, 13.1%-35.4%, increased yield by 7.1%-28.1%, 1.5%-34.1%, 8.5%-28.9%, and increased PC by 5.6%-10.4%, 10.1%-17.7% and 8.5%-15.6%, respectively. (3) The effects of nitrogen rates on grain yield and protein yield followed a quadratic curve relationship, and the fitting equation showed that the optimal nitrogen application rates for stable yield and quality of winter wheat were 189-202, 116-124 and 161-174 kg·hm-2 in humid, dry and normal years, respectively. 【Conclusion】On the whole, the best nitrogen application schemes were 189-202, 116-124 and 161-174 kg·hm-2 in humid, dry and normal years, respectively. And the management model of “nitrogen applying according to rain” was adopted, which was “the amount of basic nitrogen fertilizer was determined by SWSs, while the top dressing of nitrogen fertilizer was determined by rainfall from sowing to jointing stage”. The model could not only meet the requirements of stable yield and high quality of winter wheat, but also ensure water high-efficient use.

    HORTICULTURE
    Effects of Two Microbial Agents on Yield, Quality and Rhizosphere Environment of Autumn Cucumber Cultured in Organic Substrate
    WANG JunZheng,ZHANG Qi,GAO ZiXing,MA XueQiang,QU Feng,HU XiaoHui
    Scientia Agricultura Sinica. 2021, 54(14):  3077-3087.  doi:10.3864/j.issn.0578-1752.2021.14.013
    Abstract ( 355 )   HTML ( 35 )   PDF (528KB) ( 247 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】A bag-cultivated experiment was conducted to test the potential of two microbial agents (Bacillus methylotrophicus and Lactobacillus plantarum) on the yield, quality and rhizosphere environment of autumn cucumber using organic substrate as a growth medium. The results of the experiment were helpful to develop the new functional type of microbial agents to guide the high quality and high efficiency production of cucumber. 【Method】Cucumber cultivar ‘Bonai 526’ was sown in a bag of organic substrate. There were a total of six treatments, including CK1 (Irrigation of water and non-application of microbial agents), NT1 (Irrigation of water and application of Bacillus methylotrophicusagents), NT2 (Irrigation of water and application of Lactobacillus plantarumagents), CK2 (Irrigation of nutrient solution and non-application of microbial agents), FT1 (Irrigation of nutrient solution and application ofBacillus methylotrophicusagents), and FT2 (Irrigation of nutrient solution and application of Lactobacillus plantarumagents). The microbial agents (Bacillus methylotrophicus as ‘VL-10’ and Lactobacillus plantarum as ‘LYS-1’) were added to the cucumber rhizosphere at the rate of 2.5×1010 CFU/plant on 20th, 35th and 50th day. 【Result】Compared with CK1 and CK2 treatment,Bacillus methylotrophicus inoculation in both of the nutrient levels increased the dry matter content by 13.51% and 15.02%, and yield by 20.83% and 15.63%, respectively; Lactobacillus plantarum increased the dry matter content by 11.43% and 8.42%, and the yield by 17.42% and 14.96%, respectively. The fruit quality under FT1 treatment was the best; compared with CK2 treatment, its content of free amino acid, organic acid, soluble sugar, reducing sugar and vitamin C were significantly increased by 10.61%, 28.93%, 22.92% and 39.88%, respectively,, followed by FT2 treatment. The accumulation of phosphorus and potassium were significantly higher under both NT1 (7.43% and 10.60%) and NT2 (13.50% and 8.19%) compared with CK1 treatment (P<0.05). Compared with CK2 treatment, the accumulation of nitrogen, and potassium were significantly higher under FT2 (24.18% and 26.25%), while the phosphorus concentration in plants were significantly higher under FT1 (17.16%). Compared with CK2 treatment, the nitrogen, phosphorus and potassium fertilizer utilization rates under FT1 treatment were significantly increased by 82.85%, 483.90% and 75.60% (P<0.05), respectively, which under FT2 treatment was significantly increased by 102.42%, 367.98% and 120.46%, respectively (P<0.05). Compared with CK1 treatment, the treatments NT1 and NT2 improved soil enzymatic activities of sucrase (100.66% and 116.60%), catalase (3.39% and 4.10%) and alkaline phosphatase (6.99% and 95.08%) at the full fruiting stage, among which catalase and alkaline phosphatase activities were still higher at the end of the experiment. Compared with CK2 treatment, the FT1 treatment significantly improved soil enzymatic activities of urease (3.75% and 13.13%), sucrase (68.62% and 31.68%) and alkaline phosphatase (18.00% and 109.64%) in the full fruiting period and the seedling pulling period, while the urease (4.95% and 6.12 %), sucrase (24.93% and 63.35%) and alkaline phosphatase activity (26.99% and 84.01%) under FT2 treatment (P<0.05). The effect ofBacillus methylotrophicagent on the urease and alkaline phosphatase activities were better than that ofLactobacillus plantarumagent. Moreover, the substrate incubated with Bacillus methylotrophic showed significantly higher available nitrogen content (63.33% and 72.70%) (P<0.05) compared with their respective controls at the full fruiting stage, and the increase was 25.48% and 86.46% at the end of the experiment, respectively. 【Conclusion】In conclusion, adding 7.5×1010CFU/plant of Bacillus methylotrophic to the substrate could improve the rhizosphere environment of cucumber and promote the absorption, assimilation and accumulation of essential elements by the plant. Moreover, it could effectively improve the yield and fruit quality of cucumber.

    Screening of ARF-Aux/IAA Interaction Combinations Involved in Apple Fruit Size
    ZHOU Zhe,BIAN ShuXun,ZHANG HengTao,ZHANG RuiPing,GAO QiMing,LIU ZhenZhen,YAN ZhenLi
    Scientia Agricultura Sinica. 2021, 54(14):  3088-3096.  doi:10.3864/j.issn.0578-1752.2021.14.014
    Abstract ( 490 )   HTML ( 39 )   PDF (3387KB) ( 231 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】The transcriptomics and bioinformatics methods were used to carry out the genome-wide potential interacting MdARFs and MdIAAs pairs screening, so as to build a foundation for clarifying the function of related genes and elucidating the molecular mechanism underlying auxin-regulated apple fruit size.【Method】 Two apple genotypes, Royal Gala (WT) and 35S::miRNA172p transgenic Royal Gala (miR172OX ) were used as test materials in this study. Fruit materials from different developmental stages and tissue types were collected for both genotypes and were subjected to transcriptome sequencing. Clean reads were aligned to the reference genome and the differential expression analysis was performed. Based on the transcriptome data obtained by sequencing the transgenic small fruit and wild-type large fruit, the pairwise expression analysis was performed across MdARFs and MdIAAs families. The amino acid sequences of 23 ARFs and 34 Aux/IAAs were downloaded from Arabidopsis genome, and 21 ARFs and 25 Aux/IAAs were downloaded from tomato genome, which were further compared with candidate MdARFs and MdIAAs to construction phylogenetic trees. The MEME and TBtools were used to carry out the Motif analysis for candidate MdARFs and MdIAAs. Pairs with high interacting possibilities were further confirmed by a protein-protein interacting network constructed in apple to finalize combinations with the highest probability of involvement in fruit development. 【Result】 The whole fruit at 2 weeks post full bloom and the fruit skin, fruit flesh and fruit core at 4 WPFB were collected from WT and miR172OX, respectively. To achieve research objectives, transcriptome sequencing was carried out. A total of 178.19 Gb paired-end reads of 125 bp/150 bp were generated. All indexes indicated that the three biological replicates had highly consistent transcriptome profiles across all tissue types. FPKM values in at least one library was over 2 were used as a standard to eliminate the low expressed genes, so a total of 38MdARFs and 27 MdIAAs were expressed. In our fruit developmental transcriptome data, eight pairs of MdARF-MdIAA were obtained through Pearson correlation analysis, whose Pearson correlation coefficient was over 0.9 or below -0.9. The systematic phylogenetic analysis showed that MdARF6 and MdARF19 belonged to the same branch with AtARFs, which played a role in transcription activation, while MdARF2, MdARF4, and MdARF9 were closely related to transcriptional inhibitory AtARFs. Motif analysis results showed that both the candidate MdARFs and MdIAAs proteins contained Motif 2 and Motif 5, which were corresponded to the conserved domains Motif IV and Motif III in the IAA protein, respectively. After homolog mapping inspection with Arabidopsis, two potential MdARF-MdIAA interacting pairs were selected for future functional identification. 【Conclusion】Among apple MdARF and MdIAA family members, eight pairs of MdARF-MdIAA showed significant correlations in terms of their expression patterns during fruit development. Further homology mapping confirmed two pairs of them, including MdARF4-MdIAA17 and MdARF4-MdIAA19, were most likely to participate in the regulation of apple fruit development through mediating auxin signal transduction.

    Genome-Wide Identification and Expression Pattern Analysis of LRR-RLK Gene Family in Apple
    HUANG JinFeng,LÜ TianXing,WANG Xu,WANG YingDa,WANG DongMei,YAN ZhongYe,LIU Zhi
    Scientia Agricultura Sinica. 2021, 54(14):  3097-3112.  doi:10.3864/j.issn.0578-1752.2021.14.015
    Abstract ( 854 )   HTML ( 68 )   PDF (9601KB) ( 657 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】 The study was carried out to explore the whole genome characteristics and expression patterns ofLRR-RLKs in apple, to reveal the expression specificity of family members in different tissues and their responses to biological and abiotic stress, and further understand its biological function in apple. 【Method】 The members of LRR-RLK gene family in the whole genome of apple were identified based on the local BLAST database and Pfam database. The LRR-RLK amino acid sequence prediction, subcellular localization prediction, domain analysis, phylogenetic tree and chromosome localization were completed by software of ExPASy Proteomics Server, Cell-PLoc, CD-Search Tool, MEGAX and MG2C. In addition, the expression pattern of 12 LRR-RLK genes in different tissues and stress were analyzed by real-time fluorescent quantitative PCR (qRT-PCR).【Result】 A total of 378 LRR-RLKgenes were identified from apple genome. TheseLRR-RLKgenes encoded proteins containing 318-1 827 amino acid, and the theoretical isoelectric point ranged from 6.14 to 9.01. The prediction subcellular localization of apple LRR-RLK proteins was all distributed in the cell membrane. The gene family could be divided into 15 subgroups, containing 1-111 genes. The 378 genes in this family were distributed on all 17 chromosomes of the apple, and the chromosome 7 contained 40 genes. The LRR-RLK gene family had two conserved domains, namely the leucine-rich repeat structure and the protein kinase domain. Irregular curl and α-helix was the main secondary structure in the LRR-RLK gene family, and the rotation of β-turn was very small. It was found that the 12 selected family members were expressed in all tissues (except MD00G1105400), and most genes were expressed at relatively high levels in stem. Seven genes were up-regulated under low temperature conditions, and the expression of MD09G1153800 was the most obvious. The expression of MD09G1153800 was raised to 6.8 times of that under the control. While MD06G1170200 and MD05G1061600 were both down-regulated. Eight genes were up-regulated under drought conditions, and MD00G1105400 was the most obvious one. The expression of MD00G1105400 was raised to 9.6 times of that under the control; under salt conditions, MD04G1150400, MD13G1108000 and MD02G1071800 were always up-regulated. Among them, MD02G1071800 had the highest expression after 24 hours of salt stress treatment, which was 14.9 times of that under the control. After inoculating Botryosphaeria dothidea, the expression of 12 LRR-RLKgenes increased first and then decreased. After inoculating, the expression level of Wangshanhong was higher on the first day, however the expression level of Jiguan was higher on the third day. The expression of MD09G1153800 and MD05G1065800 were up-regulated significantly in Jiguan, in relative to have no change in Wangshanhong, suggesting that these twoLRR-RLKscould serve as candidate genes for further functional characterization. 【Conclusion】 A total of 378 LRR-RLK members were identified from apple whole genome sequences, which could be divided into 15 groups and distributed on 17 chromosomes, and the most of LRR-RLK genes were responsive to stress stimulus and Botryosphaeria dothidea.

    ANIMAL SCIENCE·VETERINARY SCIENCE·RESOURCE INSECT
    Screening, Identification and Functional Analysis of Important LncRNAs for Lactation Traits in Small-Tailed Han Sheep
    WANG JiQing,HAO ZhiYun,SHEN JiYuan,KE Na,HUANG ZhaoChun,LIANG WeiWei,LUO YuZhu,HU Jiang,LIU Xiu,LI ShaoBin
    Scientia Agricultura Sinica. 2021, 54(14):  3113-3123.  doi:10.3864/j.issn.0578-1752.2021.14.016
    Abstract ( 377 )   HTML ( 42 )   PDF (1128KB) ( 297 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Long non-coding RNAs (lncRNAs) are a type of non-coding RNAs with >200 nt in length, which have been shown to regulate mammary gland development and lactation process in dairy cows and dairy goats. However, little is known about the effect of lncRNAs on milk traits in sheep. The aim of the study was to analyze the effect of lncRNAs on milk performance and then provided a theoretical basis for elucidating molecular mechanism of lactation performance in sheep.【Method】Three high-lactating yield and high-milk-fat-content Small-Tailed Han sheep and three low-lactating yield and low-milk-fat-content Small-Tailed Han sheep were selected to profile the expression of lncRNAs in the mammary gland tissues during lactation using RNA-Seq. The enrichment analysis was performed using GO and KEGG databases for the target genes of differentially expressed lncRNAs between the two groups. The expression levels of 16 differentially expressed lncRNAs were verified using reverse transcription-quantitative PCR (RT-qPCR).【Result】A total of 7 239 expressed lncRNAs were identified in the mammary gland tissues of Small-Tailed Han sheep, including 2 262 known lncRNAs and 4 977 novel lncRNAs. The most of lncRNAs were expressed at low levels. The 120 differentially expressed lncRNAs were found between the two groups of Small-Tailed Han sheep, of which 68 lncRNAs were up-regulated in high-lactating performance Small-Tailed Han sheep, while 52 lncRNAs were down-regulated. The target genes of differentially expressed lncRNAs were significantly enriched in sulfur compound metabolic process, thioester biosynthesis process, acyl-CoA biosynthetic process, Rap1 signal pathway and adhesion junction. The lncRNA- miRNA network showed that some target miRNAs of the six most differentially expressed lncRNAs including MSTRG.125242.6 and MSTRG.59580.8, play important roles in mammary gland development and lactation in domestic animals. The RT-qPCR results showed that the expression tendency of 16 lncRNAs was consistent with the RNA-Seq results, which confirmed the accuracy and authenticity of the RNA-Seq data.【Conclusion】The differentially expressed lncRNAs screened were involved in the regulation of mammary gland development and milk performance in sheep and the results will provide a theoretical basis for analyzing the molecular genetic mechanism of lactation performance in sheep.

    Effect and Mechanism of Tea Tree Oil on LPS Induced Mastitis in Dairy Cows
    CHEN Zhi,ZHANG Yi,LU QinYue,GUO JiaHe,LIANG Yan,ZHANG MingYiXing,YANG ZhangPing
    Scientia Agricultura Sinica. 2021, 54(14):  3124-3133.  doi:10.3864/j.issn.0578-1752.2021.14.017
    Abstract ( 385 )   HTML ( 31 )   PDF (1253KB) ( 197 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】Cow mastitis has been one of the biggest challenges in dairy farming and dairy products industry, which restricts the healthy development of dairy industry. Effective prevention and treatment of cow mastitis can provide a good guarantee for the health of cows and the production of high-quality dairy products. This experiment explored the effects of tea tree oil on LPS-induced mastitis in dairy cows, and explored the feasibility of using tea tree oil instead of antibiotics to treat mastitis in dairy cows. This experiment provides a reference for the treatment of dairy cow mastitis with tea tree oil. 【Methods】The cells in good condition, were selected, which were added 50, 100, 200, 500, and 1 000 μg·mL-1 LPS respectively in the culture of these cells to detect the relevant indicators. The tea tree oil and LPS were added to the model for co-culture. The model of dairy cow mastitis induced LPS was established by the CCK-8 method, flow cytometry, real-time fluorescence quantification and ELISA assay. Antagonistic effect of tea tree oil on LPS in dairy cow mastitis cell model: 0.0002%, 0.0004%, 0.0006%, 0.0008%, 0.001%, 0.002%, 0.004%, 0.006%, 0.008% and 0.01% tea tree oil were added to the dairy cow mastitis cell model induced by 200 μg·mL-1 LPS for 12 hours to detect the related indexes. 【Result】CCK-8 method was used to detect the cell proliferation activity. The results showed that under the condition of 100 μg·mL-1 LPS poisoning, the activity of the cells began to decline in varying degrees. There was not a large number of apoptosis in 100 μg·mL-1 LPS after 12 hours of induction, while about 46% of the cells showed early and late apoptosis in 200 μg·mL-1 LPS. 200 μg·mL-1 LPS induced for 12 hours was the best condition for the establishment of mastitis model. The results also showed that when tea tree oil concentration was 0.0004%, 0.0006% and 0.0008%, the apoptosis rate of the cells decreased. Among them, when tea tree oil concentration was 0.0006%, the protection effect was the most obvious. The proportion of living cells was 71.95%, the proportion of early apoptotic cells was 22.15%, and the proportion of late apoptotic cells was 5.11%; compared with the living cells, the proportion of the mastitis model group of tree oil increased by about 22%. After that, the expression of cytokines and apoptotic factors were detected by qPCR in the three groups with protective effect. With the increase of the concentration of tea tree oil, the expression of TNF-α was down regulated more, the expression of IL-6 was down regulated less (P< 0.01), and the expression of STAT1 was up regulated slightly when 0.0004% tea tree oil was added, while down regulated slightly when 0.0006% and 0.0008% tea tree oil were added, and the expression of tea tree oil with 0.0006% concentration was the lowest (P< 0.05). The expression of NF-κB, MAPK and caspase-3 was significantly reduced in the three groups of tea tree oil adding concentration. Among them, the expression of inflammatory response protein in 0.0006% tea tree oil group was the lowest, about 50% of that in the blank control group. The protein expression was almost the same, about 55% of the blank control group (P< 0.05).【Conclusion】Tea tree oil had a certain antagonistic effect on LPS within the appropriate concentration range, which could reduce the proportion of apoptosis, improve the survival proportion of normal cells, and down regulate the expression of inflammatory factors, apoptosis factors and corresponding proteins.

    RESEARCH NOTES
    Investigation of miR-486 Target Genes in Skeletal Muscle of Bashbay Sheep in Different Development Periods
    ZHANG Wei,WANG ShiYin,GAO Li,YANG LiWei,DENG ShuangYi,LIU XiaoNa,SHI GuoQing,GAN ShangQuan
    Scientia Agricultura Sinica. 2021, 54(14):  3134-3148.  doi:10.3864/j.issn.0578-1752.2021.14.018
    Abstract ( 361 )   HTML ( 34 )   PDF (1826KB) ( 231 )   Save
    Figures and Tables | References | Related Articles | Metrics

    【Objective】The aim of this study was to deeply reveal the target genes of miR-486 in skeletal muscle of Bashbay sheep during different development periods, and to get basic data for finally uncover the molecular regulation mechanism of excellent meat traits of this sheep breed, then finally to support to further breeding. 【Method】The skeletal muscle of Bashbay sheep were collected during 40 d, 50 d, 60 d, 80 d, 100 d and 120 d of fetal period, and 1, 2, and 3 months old of new born sheep. The total RNA was extracted, and two cDNA libraries of miR-486 target genes were constructed using mRNA of fetal and postnatal period, respectively, then the library was sequenced applied high-throughput sequencing technology. Based on the function analysis, 10 candidate target genes were selected, and their express patterns in Bashbay sheep skeletal muscle of 10 development periods mentioned above were detected by qRT-PCR. By analyzing the expression pattern of these 10 candidate target genes and miR-486 in skeletal muscle of Bashbay sheep, their target regulation relationship and functions were primarily verified. Finally, 4 candidate target genes were chosen to confirm their regulation relationship using the double luciferase reporter assay and target regulation assay in satellite cell of skeletal muscle of Bashbay sheep. 【Result】Total 123 and 118 target genes of fetal and postnatal period were obtained, respectively. Due to these target genes were not confirmed by further function assay, so they were called candidate target genes temporarily. Among of them, 96 genes were expressed in skeletal muscle of fetal and postnatal period, 27 and 22 genes were specially expressed in these two periods respectively. The result of GO and KEGG analysis showed that these candidate target genes regulated mass pathways related to muscle cell differentiation and development, just like PI3k-Akt, MAPK, Wnt, Adherens junction and Regulation of actin cytoskeleton, etc. All 10 candidate target genes were expressed in skeletal muscle of Bashbay sheep, but their expression patterns were different. Among of them,PTEN, Foxo1, Dock3, PAX7, IGF1R, PIK3I1 and FBN1 were expressed at a relatively high level in the fetal period skeletal muscle of Bashbay sheep, but were significantly down regulated during postnatal period, and OLFM4 showed an opposite expression pattern compared with above 7 genes. The expression of ARHGAP5 and PDCD4 gene did not found significantly change in all 10 development periods of Bashbay sheep. The result of double luciferase reporter assay showed that miR-486 could bind the target sites of PTEN, Foxo1, IGF1R and PIK3I1 efficiently, and significantly suppressed the activity of firefly luciferase. In skeletal muscle satellite cell, mir-486 also could down-regulate mRNA of these four genes significantly, and finally suppressed their biology functions. So it was ultimately confirmed that miR-486 could regulate these 4 target genes indeed in skeletal muscle of Bashbay sheep. 【Conclusion】The investigation deeply revealed the target genes of miR-486 in skeletal muscle of Bashbay sheep during different development periods for the first time, and comprehensively analyzed the biology process and signaling pathway they participated. The further function assay proved that the data of target genes were reliable. These data would help to uncover the molecular regulation mechanisms related to excellent meat traits of Bashbay sheep.