• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Wetting alternating with partial drying during grain filling increases lysine biosynthesis in inferior rice grain

    2024-03-07 01:52:04YiJiangWenliTaoWeiyangZhangZhiqinWangJianchangYang
    The Crop Journal 2024年1期

    Yi Jiang, Wenli Tao, Weiyang Zhang, Zhiqin Wang, Jianchang Yang*

    Jiangsu Key Laboratory of Crop Genetics and Physiology/Co-Innovation Center for Modern Production Technology of Grain Crops, Yangzhou University, Yangzhou 225009,Jiangsu, China

    Keywords:Brassinosteroids Inferior grain Lysine biosynthesis Rice Wetting alternating with partial drying

    ABSTRACT Lysine content is a criterion of the nutritional quality of rice.Understanding the process of lysine biosynthesis in early-flowering superior grain (SG) and late-flowering inferior grain (IG) of rice would advance breeding and cultivation to improve nutritional quality.However,little information is available on differences in lysine anabolism between SG and IG and the underlying mechanism,and whether and how irrigation regimes affect lysine anabolism in these grains.A japonica rice cultivar was grown in the field and two irrigation regimes, continuous flooding (CF) and wetting alternating with partial drying (WAPD),were imposed from heading to the mature stage.Lysine content and activities of key enzymes of lysine biosynthesis, and levels of brassinosteroids (BRs) were lower in the IG than in the SG at the early grainfilling stage but higher at middle and late grain-filling stages.WAPD increased activities of these key enzymes, BR levels, and contents of lysine and total amino acids in IG, but not SG relative to CF.Application of 2,4-epibrassinolide to rice panicles in CF during early grain filling reproduced the effects of WAPD, but neither treatment altered the activities of enzymes responsible for lysine catabolism in either SG or IG.WAPD and elevated BR levels during grain filling increased lysine biosynthesis in IG.Improvement in lysine biosynthesis in rice should focus on IG.

    1.Introduction

    In rice,whose grain feeds over half of the world’s population[1],content of the amino acid lysine is a criterion of nutritional quality[2,3].Insufficient intake of lysine not only leads to metabolic,physiological, and functional diseases in humans, but hinders the uptake and use of other amino acids[3–5].Accordingly,increasing grain lysine content is a target of rice breeding and cultivation[3–6].

    Grain on a rice panicle can be divided into two types: superior grain (SG) and inferior grain (IG) [7–9].SG, located in the apical part of the panicle, usually flowers earlier, fills faster, and weighs more than the late-flowering IG, located in the basal part of the panicle [7–9].Although grain filling of SG and IG and its mechanisms have been extensively investigated [7–9], it is not known whether there is a difference in lysine content between the two types of kernel and what mechanism accounts for it.

    Lysine is synthesized in the aspartate pathway in rice kernels,and there are three key enzymes involved in its biosynthesis:aspartate kinase (AK, EC 2.7.2.4), homoserine dehydrogenase(HSDH; EC 1.1.1.3), and dihydropicolinate synthase (DHDPS, EC 4.2.1.52).[10–12].In contrast, lysine α-ketoglutarate reductase(LKR; EC 1.5.1.8) and saccharopine dehydrogenase (SDH, EC 1.5.1.9) are considered bi-functional lysine-degrading enzymes[11–13].Either increasing activities of AK, HSDH, and DHDPS or inhibiting LKR and SDH activities increased lysine content in cereal seeds [10–13].However, little information is available about dymamics of lysine biosynthesis and catabolism in SG and IG during grain filling, and whether such dymamics are associated with lysine content.

    It is believed [5,14,15] that plant hormones, especially brassinosteroids(BRs), function in amino acid biosynthesis in rice grain.A rice cultivar with a higher lysine content showed higher levels of grain 2,4-epibrassinolide (2,4-epiBL) and 2,8-homobrassinolide(2,8-homoBL)than one with a lower lysine content,whereas other phytohormones did not differ between the two cultivars [5,14].Applying 2,4-epiBL, a synthetic BR, increased grain lysine content,whereas applying other plant growth regulators, such as zeatin riboside or gibberellin 3, showed no effects [5,14].It is desirable to know how BRs function in SG and IG during grain filling, and how to regulate their levels to increase lysine biosynthesis there.

    Irrigation influences both the quantity and quality of rice yield[16–18].Alternate wetting and drying(AWD)has been extensively adopted in world rice production[17–19].Although it saves water,its influence on quantity and quality of rice yield and grain lysine content is disputed[19–21].In our earlier work[22–24]a wetting alternating with partial soil drying (WAPD) regime during grain filling,in which leaf photosynthesis was not impaired, accelerated grain filling, thereby increasing grain yield.It remains unclear whether a WAPD regime could increase grain lysine content,especially in IG, by increasing lysine biosynthesis.

    This study aimed to test the hypothesis that a WAPD regime not only promotes rice grain filling but increases lysine biosynthesis by elevating grain BR levels.The experimental plan was to compare SG and IG with respect to grain filling, amino acid contents, activities of lysine biosynthetic enzymes, BR levels, and responses of these characters,leaf photosynthesis, and root activity,to application of exogenous BRs and a BR inhibitor.

    2.Materials and methods

    2.1.Plant materials and cultivation

    A high-yielding rice (Oryza satica L.) cultivar Jinxiangyu 1 (a japonica inbred),was used.Seeds were sown on May 20 and seedlings were transplanted on June 15 in 2020 and 2021.The hill space was 0.25 m × 0.15 m with two seedlings per hill.The full amount of N was applied at 210 kg ha-1as urea on the day of rice transplanting,10 d after transplanting(DAT),38 DAT(young panicle differentiation), and 54 DAT(pistil and stamen differentiation)in the ratio 4:2:2:2.The respective rates of P(P2O5)and K(K2O)were 105 and 90 kg ha-1, applied on the transplanting day.Except for midseason drainage at 32–36 DAT, a water level of 4–5 cm was maintained in the field until heading (July 28 in 2020 and July 29 in 2021), when irrigation treatments were begun.

    The experimental field was located on a farm at Yangzhou University, Yangzhou, China.The soil texture, contents of N, P,and K, and field moisture capacity were described previously[25].The main weather factors during the rice-growing season in both years are presented in Fig.S1.

    Field experiments were also performed in 2018 and 2019 using the indica rice cultivar Yangdao 6.Except for the total N rate of 180 kg ha-1for this cultivar, other growth conditions were the same as Jinxiangyu 1.

    2.2.Irrigation treatment

    Irrigation treatments were applied from heading to the mature stage (74–129 DAT) of rice.They consisted of continuous flooding(CF)and wetting alternating with partial drying(WAPD).Under the CF regime, a 4–5 cm water level was maintained in the field until seven days before the plants were finally harvested,a practice consistent with farmer water management[26].In the WAPD regime,the water potential threshold was-10 kPa at 15–20 cm soil depth.The soil was naturally dried from the 4–5 cm water level to the threshold, then re-watered and naturally dried again (Fig.S2).The soil water potential was monitored with tensiometers as described previously [25].The experiment followed a randomized complete block design with three replications.Plot size was 6.25×5.10 m and plots were separated by a cement ridge(60 cm wide,40 cm tall, and 50 cm deep in the soil) built between plots to prevent water and fertilizer runoff.A removable rain shelter was built over the experimental field.

    2.3.Sampling and measurement of grain contents

    Approximately 600 panicles heading on the same day were tagged in each treatment.The SG and IG on tagged panicles were sampled separately following Yang et al.[8].Eight to ten tagged panicles from each plot were sampled every 3 d from anthesis to the final harvest.Of the sampled kernels (both SG and IG), 80%were used for measurement of amino acids, enzymatic activities,and BR levels, and 0.2 were used for measuring grain filling rate.The contents of lysine,aspartate,and total amino acids were measured following Zhang et al.[27] using a high-speed amino acid analyzer(Biochrom 30,Cambridge,UK).Activities of lysine biosynthetic and catabolic enzymes were determined as described previously:AK and HSDH,Azevedo et al.[28];DHDPS,Varisi et al.[29];LKR and SDH,Gaziola et al.[30].Both 2,4-epiBL and 2,8-homoBL in kernels were extracted and purified following Ding et al.[31], and quantified following Chen et al.[32].The grain-filling processes and increases of lysine and total amino acid contents in grain were fitted with Richards’ equation [33]:

    where W is grain weight,C is content of lysine or total amino acids,A is maximum grain weight,or maximum content of lysine or total amino acids,t is time after anthesis(days),and B,k,and N are coefficients of the exponential expression.The active grain-filling period and active period of lysine or total amino acid increase (D), the mean increase rate (Gmean) and maximum increase rate (Gmax) of lysine or total amino acids, and days to reach the maximum increase rate of lysine or total amino acids(Dmax.g) were calculated using the equation following Zhu et al.[34].

    Leaf water potential and photosynthetic rate of the flag leaf,root biomass, and root oxidation activity (ROA) were measured 10 d after heading (DAH) (D1) and 30 DAH (D2) in 2020 and 10 DAH (D1) and 29 DAH (D2) in 2021, respectively, when soil water potential was approximately-10 kPa in the WAPD.They were also measured 12 DAH (W1) and 32 DAH (W2) in 2020 and 12 DAH(W1) and 31 DAH (W2) in 2021, during the wetting period.Measurement details of these parameters were described previously[25].

    At the mature stage (October 22 in both 2020 and 2021, 129 DAT),plants were harvested for measurement of rice yield components, quantity, and quality.Yield components were determined from the 20 plants and grain yield was measured for plants harvested in an area of 5.4 m2in each plot following Wang et al.[25].Approximately 800 g of grain from each plot was used for quality measurement following Zhang et al.[35].

    2.4.Chemical treatment

    At 4,8,and 12 DAH,10 μmol L-12,4-epiBL,10 μmol L-1brassinazole (BRZ, an inhibitor of BR biosynthesis), or 10 μmol L-12,4-epiBL + 10 μmol L-1BRZ were sprayed on panicles at the rate of 50 mL m-2for each application.Control plants were treated with the same volume of deionized water.Leaves were wrapped with plastic film during spraying to reduce the effect of the chemicals on leaf photosynthesis.Each chemical treatment was applied to three plots as replications with 6.3 m2(3.0 m × 2.1 m) for each plot.Activities of AK, HSDH, DHDPS, SDH, and LKR were determined 6,10,and 14 d after anthesis and contents of lysine and total amino acids in milled rice were measured as described above after final harvest.

    2.5.Statistical analysis

    SPSS 20.0(IBM,Armonk,NY,USA)was used to perform analysis of variance.Treatment means were compared by least significant difference (LSD).Pearson’s correlations among BR levels, enzymatic activities, and grain lysine content were calculated.

    3.Results

    3.1.Grain yield, grain filling rate, and milling, appearance,and eating qualities

    In comparison with the CF, WAPD increased grain yield by 11.6%–12.3%, owing mainly to increases in seed-setting rate and 1000-kernel weight (Table 1), leaf photosynthesis and ROA(Figs.S3,S4).Because irrigation treatments were initiated at heading time, they did not affect numbers of panicles and spikelets(Table 1).Similar influences of irrigation regime on yield were observed in Yangdao 6 (Table S1).

    As shown in Fig.S5,WAPD markedly increased grain filling and grain weight of the IG in comparison with the CF.No differences between the two irrigation regimes in the filling and weight of SG were observed (Fig.S5A–D; Table S2).

    Similarly, WAPD markedly improved milling, appearance, and taste qualities of the IG (Tables S3, S4, S5).The SG did not differ in these quality traits between CF and WAPD.

    3.2.Contents of lysine and total amino acids and their increasing rates

    Contents of lysine and total amino acids in both SG and IG,based on mg g-1dry weight(DW),first increased and reached their peak values at 9–12 d after anthesis(DAA)for the SG and at 18–21 DAA for the IG,and then sharply decreased(Figs.1A,B,2A,B).If the contents were expressed as μg or mg kernel-1, their changes or increase rates (μg kernel-1d-1) showed trends similar to those of grain-filling rates (Figs.1C–F, 2C–F).During the rapid increase period (3–21 DAA for the SG and 9–33 DAA for the IG), the rates of increase of the mean(Gmean)and maximum(Gmax),and the contents of lysine and total amino acids for the IG were greater in the WAPD than in the CF(Figs.1A–F,2A–F;Tables S6,S7).However,in the SG they did not differ between the CF and the WAPD.Similarly,contents of lysine and total amino acids,either mg g-1DW or μg or mg kernel-1, in the milled IG were markedly increased in the WAPD in comparison with those in the CF, but did not differ in the milled SG between the CF and the WAPD (Table S8).

    To determine whether the increase in lysine content was due to the increase in aspartate content in the IG in the WAPD,changes in aspartate content were measured in both SG and IG.As shown in Fig.S6, aspartate content in either SG or IG, expressed as either mg g-1DW or μg kernel-1, did not differ between CF and WAPD(Fig.S6A–D).Thus aspartate, the substrate for lysine synthesis,did not account for the difference in lysine content in the IG between the CF and the WAPD.As shown in Fig.S7, the mean increase rate (Gmean) and maximum increase rate (Gmax) of lysine or total amino acids during the rapid increase period were positively correlated(r=0.92**to 0.99**,Fig.S7A–D)with the contents(μg or mg kernel-1) of lysine and total amino acids in milled rice,suggesting that increased lysine biosynthesis during grain filling contributed to the increases in lysine and total amino acids in the milled IG under the WAPD regime.

    3.3.Activities of key enzymes in lysine anabolism

    Activities of the lysine biosynthetic enzymes AK, HSDH, and DHDPS showed similar changing trends during grain filling, first increasing to maximums at 9 DAA for SG and 18 DAA for IG and then decreasing(Fig.3A–F).Relative to IG,SG showed higher activities of these enzymes at early grain filling (3–12 DAA), but lower activities from 15 DAA.During 3–33 DAA, WAPD increased enzymatic activities in IG but not in SG, relative to the CF (Fig.3A–F).The activities of the three enzymes were closely correlated with lysine and total amino acid contents in the grain during the active grain-filling period (3–33 DAA) (Fig.S8A–F).

    Activities of the lysine catabolic enzymes SDH and LKR and the ratio of SDH to LKR(SDH/LKR),reached their peak value at 18 DAA for the SG and at 33 DAA for the IG (Fig.S9A–F),in contrast to the changing trend of lysine content.Activities of SDH and LKR and SDH/LKR did not differ in either SG or IG between the CF and the WAPD (Fig.S9A–F), indicating that irrigation regime did not influence lysine catabolism.

    3.4.BR levels and their relationship with lysine biosynthesis

    In agreement with the activities of AK, HSDH, and DHDPS, the levels of 2,4-epiBL and 2,8-homoBL reached their maximum at 9 DAA in the SG and at 18 DAA in the IG(Fig.4A–D).WAPD increased 2,4-epiBL and 2,8-homoBL levels in the IG but did not increase them in the SG, relative to the CF (Fig.4A–D).In the active grainfilling period, levels of 2,4-epiBL or 2,8-homoBL were correlated(r = 0.64**to 0.93**, P < 0.01) with the activities of AK, HSDH, and DHDP and the contents of lysine and total amino acids(Fig.S10A–J).Application of 2,4-epiBL to panicles in the CF regime at early grain filling increased the activities of AK, HSDH, and DHDPS, grain weight of IG, and the contents of lysine and total amino acids in IG, whereas applying BRZ produced the opposite effects (Tables 2, 3).The negative effect of applying BRZ on lysine biosynthesis was eliminated by simultaneous application of 2,4-epiBL and BRZ.Although BRZ reduced the activities of AK, HSDH,and DHDPS, grain weight of the SG, and the contents of lysine as well as total amino acids in the SG, application of 2,4-epiBL or 2,4-epiBL + BRZ did not markedly affect lysine biosynthesis and amino acid content in SG (Tables 2, 3).Applying either 2,4-epiBLor BRZ did not affect the activities of SDH and LKR in either SG or IG(Table 2).

    Table 1 Grain yield and its components in rice cv.Jinxiangyu 1.

    Fig.1.Lysine content based on dry weight(A and B)and based on per kernel(C and D)and lysine increase rate(E and F)of rice grain(cv.Jinxiangyu 1).SG,superior grain;IG,inferior grain.CF and WAPD represent respectively continuous flooding and wetting alternating with partial drying,during grain filling.The lysine increase rate is fitted with Richards’ equation [33].Values in (A–D) are mean ± SE.

    In Yangdao 6 the influence of irrigation regimes on lysine content in milled rice, changes in enzymatic activities and BR levels in filling grain,and effects of chemical applications on grain weight and lysine biosynthesis and catabolism in either SG or IG were similar (Tables S9–S13).

    4.Discussion

    4.1.Differences in lysine biosynthesis between SG and IG

    In comparison with SG,IG showed a smaller mean increase rate(Gmean),a smaller maximum increase rate(Gmax),and a lower content of lysine based on either per unit dry weight or per grain at early grain filling and at maturity expressed as μg kernel-1(Figs.1,2; Table S8).Three conditions can result in a lower content of lysine in cereal seeds: insufficient substrate supply, lower ability of biosynthesis,and higher catabolic activity[4–6,12–14].The content (mg g-1DW) of aspartate, the substrate for lysine biosynthesis, did not differ between SG and IG at the early filling stage, and was higher in IG than in SG at later filling stages(Fig.S6).The activities of the two key enzymes in lysine catabolism, SDH and LKR,were lower in IG than in SG before 24 DAA (Fig.S9).Neither the adequate supply of substrate for lysine biosynthesis nor the reduced activities of lysine catabolic enzymes responsible for lysine catabolism explain the lower lysine content in IG during the active grain-filling period.The lower activities of the three lysine biosynthetic enzymes in IG than in SG at the early grainfilling stage and the correlation of these activities with grain lysine content suggests that the lower content of lysine in the IG was due mainly to lower biosynthetic ability at the early grain-filling stage.

    Fig.2.Contents of total amino acids based on dry weight(A and B)and based on per kernel(C and D)and total amino acids increase rate(E and F)of rice grain(cv.Jinxiangyu 1).SG, superior grain; IG, inferior grain.CF and WAPD represent respectively continuous flooding and wetting alternating with partial drying, during grain filling.The total amino acid increase rate was fitted with Richards’ equation [33].Values in (A–D) are presented as mean ± SE.

    The mechanism explaining differences in lysine synthesis between SG and IG during grain filling is not fully understood.Many factors, such as source–sink relationship, carbon and nitrogen anabolism, and hormonal levels, could contribute to differences in lysine biosynthesis [12,27,36,37].Levels of 2,4-epiBL and 2,8-homoBL were much lower in IG than in SG during the active grain-filling period, in agreement with activities of AK,HSDH, and DHDPS and with grain lysine content (Figs.3, 4, S10).Application of 2,4-epiBL increased and that of BRZ reduced activities of these enzymes and lysine content in IG (Tables 2, 3, S12,S13).The lower levels of the BRs in IG than in SG during grain filling and the respective increases and reductions in lysine biosynthetic activity produced by application of the BR and inhibitor suggest that BRs regulate lysine synthesis and that a lower level of BRs leads to lower lysine synthesis in IG.Little is known, however,about how BRs regulate lysine biosynthesis.There are reports[5,14,37,38]that BRs mediate many biological processes including plant growth and development, organogenesis, plant configuration, fertility conversion, and stress resistance.It is possible that BRs could increase signal transduction in the aspartate-to-lysine pathway and expression of genes undertaking lysine synthesis,and accordingly, modulate amino acid biosynthesis [5,38].

    The finding that the per-kernel lysine or total amino acid content was much greater in the milled SG than in the milled IG, but the per-dry-weight content was lower, may be explained by the lower grain dry weight and higher protein content of milled IG than milled SG (Fig.S5; Table S4).Because the content of amino acids including lysine is usually proportional to protein content in a cereal seed and is also influenced by grain weight [3–5], a lower grain dry weight can lead to a lower per-kernel content of amino acids, whereas a higher protein content and lower dry weight can lead to a higher dry-weight content of amino acids in milled IG.

    4.2.The role of a WAPD regime in regulating lysine anabolism

    Fig.3.Activities of aspartate kinase(AK,A and B),homoserine dehydrogenase(HSDH,C and D),and dihydropicolinate synthase(DHDPS,E and F)in rice grain(cv.Jinxiangyu 1).SG,superior grain;IG,inferior grain.CF and WAPD represent respectively continuous flooding and wetting alternating with partial drying,during grain filling.Values are mean ± SE.

    The finding that in comparison with the CF,the WAPD imposed during grain filling increased not only lysine but total amino acid content in IG (Figs.1, 2; Table S8) suggests that promoting lysine synthesis also increases other amino acid biosynthesis in grain.In the WAPD regime, increased grain filling rate, activities of enzymes responsible for lysine synthesis,and BRs levels in IG led to the increases in contents of amino acids including lysine(Figs.3,4,S5).Increases in leaf photosynthesis(Fig.S3C,D)and ROA(Fig.S4C,D)under the WAPD regime could also lead to an increase in amino acid biosynthesis in the grain by supplementing carbohydrates or providing more nutrients [5,21,39].The WAPD regime not only increased seed-setting rate and grain weight,leading to an increase in grain yield,but improved milling,appearance,cooking,and taste qualities of the IG (Tables 1, S1–S5).We conclude that a WAPD regime during grain filling can simultaneously improve rice grain yield and quality, in particular by increasing lysine content in the IG.

    The WAPD regime or application of 2,4-epiBL did not markedly affect lysine biosynthesis and contents of lysine and total amino acids in the SG(Figs.1–3;Tables 2–3,S8-S13).One likely interpretation is that the levels of phytohormones such as BRs, and ability of amino acid biosynthesis are high in SG, and this grain could be insensitive or even unresponsive to the WAPD regime or 2,4-epiBL application.In another study [40], levels of polyamines,numbers of endosperm cells, and cell division rate differed among rice cultivars in IG but not in SG,and an effect of polyamine application on endosperm cell division and grain filling was observed only in IG.These findings suggest that IG has greater potential to increase grain quantity and quality than SG and that improvement in rice grain filling and amino acid biosynthesis should focus on IG.

    Fig.4.Contents of 2,4-epibrassinolide (2,4-epiBL, A and B) and 2,8-homobrassinolide (2,8-homoBL, C and D) in the grain of rice (cv.Jinxiangyu 1).SG, superior grain; IG,inferior grain.CF and WAPD represent respectively continuous flooding and wetting alternating with partial drying, during grain filling.Values are mean ± SE.

    Table 2 Effects of chemical applications on activities of enzymes involved in lysine biosynthesis and catabolism in filling grain of rice (cv.Jinxiangyu 1).

    Usually, higher activitees of AK, HSDH, and DHDPS are accompanied by higher activities of SDH and LKR to counterbalance increased lysine synthesis [11–13].But in the present study, the activities of AK, HSDH, and DHDPS were increased in the IG,whereas the activities of SDH and LKR were not greatly increased in either SG or IG under the WAPD regime (Figs.3, S9).Probably more lysine is incorporated into protein under such conditions[3–5].Thus, an increase in in lysine biosynthesis is not necessarily accompanied by an increase in lysine catabolism,such as under the WAPD regime.

    Applying BRZ reduced grain weight, contents of lysine and total amino acids,and activities of AK,HSDH,and DHDPS,whereas applying 2,4-epiBL or BRZ did not change the activities of SDH and LKR,in both SG and IG(Tables 2,3,S12,S13).A possible explanation is that a higher level of BRs is necessary to maintain a higher grain-filling rate and a higher activity of lysine biosynthesis in the two types of filling kernels.Decreases in BRs levels, such as by application of BRZ,can inhibit grain filling and reduce the activity of lysine biosynthesis,leading to a reduction in grain weight and contents of amino acids not only in IG but in SG.The observation that application of 2,4-epiBL or BRZ did not affect the activities of SDH and LKR may be explained by an assumption that BRs function mainly in mediating lysine biosynthesis and not in lysine catabolism.

    Table 3 Effects of chemical applications on grainweight andcontents of lysine andtotalaminoacidsin milled rice (cv.Jinxiangyu 1).

    5.Conclusions

    A WAPD regime imposed during grain filling elevated BR levels and increased lysine biosynthesis,increasing contents of lysine and total amino acids in IG.BRs function in mediating lysine biosynthesis.Neither WAPD or application of a lysine biosynthesis inhibitor altered lysine catabolism.Improvement of amino acid biosynthesis in rice should focus on IG.

    CRediT authorship contribution statement

    Yi Jiang:Investigation,Formal analysis,Visualization,Writing–original draft.Wenli Tao:Investigation,Formal analysis,Visualization.Weiyang Zhang:Investigation, Data curation, Formal analysis.Zhiqin Wang:Investigation, Formal analysis.Jianchang Yang:Funding acquisition, Project administration, Conceptualization, Supervision, Resources, Methodology, Writing – review &editing.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    This work was supported by the National Natural Science Foundation of China (32071943, 32272198).

    Appendix A.Supplementary data

    Supplementary data for this article can be found online at https://doi.org/10.1016/j.cj.2023.11.008.

    国产激情久久老熟女| 制服诱惑二区| 黑丝袜美女国产一区| 欧美97在线视频| 国产毛片在线视频| 久久久久久久久久人人人人人人| 肉色欧美久久久久久久蜜桃| 两性夫妻黄色片| 亚洲欧美成人综合另类久久久| 日韩中字成人| 国产成人av激情在线播放| 国产精品久久久久久av不卡| 国产精品久久久久久久久免| 国产av一区二区精品久久| 国产一区二区三区综合在线观看| 欧美精品人与动牲交sv欧美| 欧美亚洲 丝袜 人妻 在线| 国产成人91sexporn| 久久综合国产亚洲精品| 九色亚洲精品在线播放| 欧美精品高潮呻吟av久久| 人妻一区二区av| 高清黄色对白视频在线免费看| 香蕉国产在线看| 国产精品亚洲av一区麻豆 | 精品久久久精品久久久| 日韩制服丝袜自拍偷拍| 欧美人与善性xxx| 国产女主播在线喷水免费视频网站| 国产精品一国产av| 精品第一国产精品| 一边摸一边做爽爽视频免费| 国产xxxxx性猛交| 国产在线免费精品| 一区在线观看完整版| 欧美国产精品一级二级三级| 免费在线观看视频国产中文字幕亚洲 | 日韩av免费高清视频| 久久 成人 亚洲| 一边摸一边做爽爽视频免费| 另类亚洲欧美激情| 91午夜精品亚洲一区二区三区| 精品一区二区免费观看| 男男h啪啪无遮挡| 在线观看www视频免费| 一级爰片在线观看| 精品国产露脸久久av麻豆| 69精品国产乱码久久久| 免费观看性生交大片5| 97精品久久久久久久久久精品| 国产一区有黄有色的免费视频| 国产一区二区激情短视频 | 亚洲国产欧美日韩在线播放| 亚洲精品一二三| 亚洲第一区二区三区不卡| 婷婷色综合大香蕉| 我的亚洲天堂| 久久久久久免费高清国产稀缺| 久久国产精品男人的天堂亚洲| 90打野战视频偷拍视频| 91午夜精品亚洲一区二区三区| 人妻系列 视频| 亚洲精品aⅴ在线观看| 日韩中文字幕欧美一区二区 | 人妻系列 视频| 老汉色∧v一级毛片| 亚洲精品自拍成人| 人妻少妇偷人精品九色| 天天操日日干夜夜撸| 欧美97在线视频| 美女午夜性视频免费| 日韩一区二区三区影片| 亚洲男人天堂网一区| 国产成人午夜福利电影在线观看| 亚洲少妇的诱惑av| 国产午夜精品一二区理论片| 国产精品.久久久| 成年女人在线观看亚洲视频| 999久久久国产精品视频| 久久久久久久大尺度免费视频| 三级国产精品片| 国产午夜精品一二区理论片| 欧美97在线视频| 欧美在线黄色| 亚洲欧美一区二区三区国产| 伦精品一区二区三区| 中国国产av一级| 伊人亚洲综合成人网| 成年动漫av网址| 在线看a的网站| 亚洲欧洲日产国产| 男女国产视频网站| 看免费av毛片| 卡戴珊不雅视频在线播放| 少妇人妻久久综合中文| 只有这里有精品99| 一级,二级,三级黄色视频| 亚洲精品美女久久久久99蜜臀 | 日韩大片免费观看网站| 中国三级夫妇交换| av.在线天堂| 久久国产精品大桥未久av| 日产精品乱码卡一卡2卡三| 美女大奶头黄色视频| av有码第一页| 精品国产露脸久久av麻豆| 大话2 男鬼变身卡| 国产探花极品一区二区| 日韩在线高清观看一区二区三区| 久久精品国产亚洲av涩爱| 亚洲国产日韩一区二区| 免费大片黄手机在线观看| 看免费成人av毛片| 免费黄频网站在线观看国产| 建设人人有责人人尽责人人享有的| 久久精品国产综合久久久| 精品久久久精品久久久| av不卡在线播放| 亚洲精品第二区| 97精品久久久久久久久久精品| 一级a爱视频在线免费观看| 老汉色∧v一级毛片| 三上悠亚av全集在线观看| 久久久久久人妻| 黑人巨大精品欧美一区二区蜜桃| 一级爰片在线观看| 欧美黄色片欧美黄色片| 国产一区二区三区综合在线观看| 肉色欧美久久久久久久蜜桃| 亚洲av国产av综合av卡| 桃花免费在线播放| 日本猛色少妇xxxxx猛交久久| 看免费av毛片| 卡戴珊不雅视频在线播放| 毛片一级片免费看久久久久| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美中文字幕日韩二区| 亚洲一区中文字幕在线| 国语对白做爰xxxⅹ性视频网站| 精品午夜福利在线看| 国产精品秋霞免费鲁丝片| 免费高清在线观看日韩| 国产福利在线免费观看视频| 美女高潮到喷水免费观看| 日本黄色日本黄色录像| 黄色视频在线播放观看不卡| 有码 亚洲区| 高清av免费在线| 国产精品人妻久久久影院| 久久99热这里只频精品6学生| 老司机影院成人| 妹子高潮喷水视频| 考比视频在线观看| 精品人妻熟女毛片av久久网站| a级毛片在线看网站| 日韩中字成人| 精品第一国产精品| 侵犯人妻中文字幕一二三四区| 欧美成人午夜免费资源| 国产毛片在线视频| 你懂的网址亚洲精品在线观看| 中文字幕人妻熟女乱码| 黑人欧美特级aaaaaa片| 热re99久久国产66热| 亚洲欧美精品综合一区二区三区 | 老女人水多毛片| 国产探花极品一区二区| 男的添女的下面高潮视频| 国产国语露脸激情在线看| 日本av手机在线免费观看| 一级黄片播放器| 免费黄频网站在线观看国产| 久久久久国产网址| 亚洲精品视频女| www.熟女人妻精品国产| 老司机影院成人| 久久国产精品男人的天堂亚洲| 黄色一级大片看看| 一区二区三区精品91| 精品一区二区免费观看| 久久青草综合色| 女性被躁到高潮视频| 天天影视国产精品| 中文字幕人妻丝袜一区二区 | 在线亚洲精品国产二区图片欧美| 热re99久久精品国产66热6| 人人妻人人添人人爽欧美一区卜| 国产片特级美女逼逼视频| 久久久久久人妻| 女性生殖器流出的白浆| 日韩人妻精品一区2区三区| h视频一区二区三区| 国产成人av激情在线播放| 久久久久网色| 18禁观看日本| 国产精品久久久久久精品古装| 久久这里只有精品19| 亚洲综合色网址| 一边亲一边摸免费视频| 青青草视频在线视频观看| 嫩草影院入口| 欧美国产精品va在线观看不卡| 女人高潮潮喷娇喘18禁视频| 精品一区二区三卡| 久久久国产精品麻豆| 香蕉丝袜av| 香蕉国产在线看| 欧美精品国产亚洲| 男女边吃奶边做爰视频| 一区二区三区乱码不卡18| 亚洲精品日本国产第一区| 成年av动漫网址| 国产一区有黄有色的免费视频| 欧美精品国产亚洲| 日韩av免费高清视频| 午夜免费鲁丝| 午夜福利影视在线免费观看| 韩国av在线不卡| 亚洲国产毛片av蜜桃av| 国产精品.久久久| 日产精品乱码卡一卡2卡三| 十分钟在线观看高清视频www| 国产人伦9x9x在线观看 | 国产毛片在线视频| 美女大奶头黄色视频| 一级毛片 在线播放| 久久韩国三级中文字幕| 韩国高清视频一区二区三区| 国产淫语在线视频| 成年美女黄网站色视频大全免费| 中国三级夫妇交换| 色播在线永久视频| 狠狠精品人妻久久久久久综合| 国产在线免费精品| 国产精品蜜桃在线观看| 亚洲欧美一区二区三区黑人 | 午夜久久久在线观看| 亚洲国产色片| 国产毛片在线视频| 久久精品国产自在天天线| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品电影小说| 999久久久国产精品视频| 热99久久久久精品小说推荐| 亚洲av欧美aⅴ国产| 26uuu在线亚洲综合色| 久久久久久久久免费视频了| 欧美人与性动交α欧美软件| 女人高潮潮喷娇喘18禁视频| 日韩熟女老妇一区二区性免费视频| 日韩av在线免费看完整版不卡| xxx大片免费视频| 亚洲色图综合在线观看| 妹子高潮喷水视频| 精品少妇黑人巨大在线播放| 免费人妻精品一区二区三区视频| 有码 亚洲区| 一区二区三区激情视频| 午夜免费男女啪啪视频观看| 日韩精品有码人妻一区| 啦啦啦啦在线视频资源| 一级,二级,三级黄色视频| xxxhd国产人妻xxx| 久久人人97超碰香蕉20202| 亚洲国产看品久久| 国产成人aa在线观看| 人妻 亚洲 视频| 2022亚洲国产成人精品| 两个人看的免费小视频| 18禁观看日本| 国产男女超爽视频在线观看| 欧美日韩成人在线一区二区| 天堂8中文在线网| 男人舔女人的私密视频| 亚洲av日韩在线播放| 色哟哟·www| 日韩一区二区三区影片| 成人亚洲精品一区在线观看| 波多野结衣一区麻豆| 久久人人爽av亚洲精品天堂| 欧美精品亚洲一区二区| 男人添女人高潮全过程视频| 亚洲人成77777在线视频| 婷婷色综合大香蕉| 搡老乐熟女国产| 9热在线视频观看99| 久久久欧美国产精品| 男女免费视频国产| 大码成人一级视频| 婷婷色综合www| 在线观看一区二区三区激情| 人体艺术视频欧美日本| 美女国产视频在线观看| av电影中文网址| 精品一品国产午夜福利视频| 国产黄色免费在线视频| 久久午夜综合久久蜜桃| 久久精品久久久久久久性| 久久久久精品人妻al黑| 性少妇av在线| 午夜激情av网站| 亚洲欧美成人精品一区二区| 亚洲色图 男人天堂 中文字幕| 久久久国产一区二区| 成人毛片60女人毛片免费| 亚洲精品aⅴ在线观看| 多毛熟女@视频| 久久久久久久国产电影| 国产高清不卡午夜福利| 色婷婷av一区二区三区视频| 毛片一级片免费看久久久久| 国产野战对白在线观看| 精品福利永久在线观看| 免费在线观看黄色视频的| 卡戴珊不雅视频在线播放| 国产在线一区二区三区精| 夫妻性生交免费视频一级片| av免费在线看不卡| 91精品国产国语对白视频| 午夜影院在线不卡| 母亲3免费完整高清在线观看 | 色婷婷av一区二区三区视频| 国产黄色视频一区二区在线观看| 赤兔流量卡办理| 日韩伦理黄色片| 9191精品国产免费久久| 亚洲视频免费观看视频| 五月伊人婷婷丁香| 国产精品嫩草影院av在线观看| 亚洲国产毛片av蜜桃av| 蜜桃在线观看..| 免费女性裸体啪啪无遮挡网站| 亚洲国产av新网站| 9191精品国产免费久久| 久久久久国产精品人妻一区二区| 午夜福利视频在线观看免费| 午夜福利乱码中文字幕| 久久精品国产亚洲av高清一级| 国产 一区精品| 日韩三级伦理在线观看| a级毛片黄视频| 国产又色又爽无遮挡免| 天天躁夜夜躁狠狠久久av| 丰满迷人的少妇在线观看| 超碰97精品在线观看| 日本黄色日本黄色录像| 国产极品粉嫩免费观看在线| 最近手机中文字幕大全| 久久久久久人妻| 日本-黄色视频高清免费观看| av一本久久久久| 九九爱精品视频在线观看| 熟女电影av网| 久久午夜综合久久蜜桃| 十八禁高潮呻吟视频| 日韩av不卡免费在线播放| 大片免费播放器 马上看| 久久这里有精品视频免费| 国产精品一国产av| 春色校园在线视频观看| 18+在线观看网站| 国产黄色免费在线视频| 久久韩国三级中文字幕| 久久久国产精品麻豆| 黑人欧美特级aaaaaa片| 99re6热这里在线精品视频| 丰满乱子伦码专区| 一级爰片在线观看| av网站在线播放免费| 免费少妇av软件| 久久久久视频综合| 美女午夜性视频免费| 精品久久蜜臀av无| 新久久久久国产一级毛片| 亚洲美女视频黄频| 久久人人97超碰香蕉20202| 99香蕉大伊视频| 久久久精品区二区三区| 日本色播在线视频| 精品99又大又爽又粗少妇毛片| 中文字幕人妻熟女乱码| 久久人人97超碰香蕉20202| 少妇被粗大猛烈的视频| 如何舔出高潮| 婷婷色综合大香蕉| 亚洲美女搞黄在线观看| 69精品国产乱码久久久| 九九爱精品视频在线观看| 最新中文字幕久久久久| 成年av动漫网址| 日韩不卡一区二区三区视频在线| 久久精品国产鲁丝片午夜精品| 国产精品秋霞免费鲁丝片| 日韩制服丝袜自拍偷拍| 女人精品久久久久毛片| 国产又爽黄色视频| 夫妻午夜视频| 成人午夜精彩视频在线观看| 久久久精品94久久精品| 国产成人午夜福利电影在线观看| 国产精品女同一区二区软件| 久久鲁丝午夜福利片| 亚洲,一卡二卡三卡| 亚洲,欧美精品.| 久久久久久久久免费视频了| www.av在线官网国产| 天美传媒精品一区二区| 日韩熟女老妇一区二区性免费视频| 少妇猛男粗大的猛烈进出视频| 日韩熟女老妇一区二区性免费视频| 精品少妇内射三级| 欧美日韩一区二区视频在线观看视频在线| 免费日韩欧美在线观看| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久精品电影小说| 男的添女的下面高潮视频| 国产黄色视频一区二区在线观看| 免费人妻精品一区二区三区视频| 香蕉精品网在线| 少妇被粗大猛烈的视频| 十八禁高潮呻吟视频| 日韩大片免费观看网站| 日韩精品免费视频一区二区三区| 搡老乐熟女国产| 国产成人一区二区在线| 日日摸夜夜添夜夜爱| 欧美日本中文国产一区发布| 国产成人a∨麻豆精品| 黑丝袜美女国产一区| 中文字幕精品免费在线观看视频| 亚洲av欧美aⅴ国产| 侵犯人妻中文字幕一二三四区| 久久99精品国语久久久| 国产欧美日韩一区二区三区在线| 国产成人免费观看mmmm| 国产精品香港三级国产av潘金莲 | 91午夜精品亚洲一区二区三区| 中文字幕制服av| 91aial.com中文字幕在线观看| 免费看av在线观看网站| 免费播放大片免费观看视频在线观看| 国产深夜福利视频在线观看| 日韩大片免费观看网站| 精品一品国产午夜福利视频| 两个人看的免费小视频| 女性生殖器流出的白浆| 99久久精品国产国产毛片| 美女脱内裤让男人舔精品视频| 超碰成人久久| 亚洲av国产av综合av卡| 女人被躁到高潮嗷嗷叫费观| 国产一区二区在线观看av| 国产乱人偷精品视频| 精品久久蜜臀av无| 国产又色又爽无遮挡免| 国产av国产精品国产| 国产男女超爽视频在线观看| 中文乱码字字幕精品一区二区三区| 男女边摸边吃奶| 妹子高潮喷水视频| 亚洲国产av影院在线观看| 亚洲国产精品国产精品| 2018国产大陆天天弄谢| 亚洲男人天堂网一区| 国产免费视频播放在线视频| 免费观看在线日韩| 人妻少妇偷人精品九色| 国产一区二区三区av在线| 欧美 日韩 精品 国产| 极品人妻少妇av视频| 欧美97在线视频| 久久久久人妻精品一区果冻| 亚洲中文av在线| 又黄又粗又硬又大视频| 亚洲精品日韩在线中文字幕| 纯流量卡能插随身wifi吗| av免费在线看不卡| av视频免费观看在线观看| 天美传媒精品一区二区| 叶爱在线成人免费视频播放| 国产黄色视频一区二区在线观看| 69精品国产乱码久久久| 久久久a久久爽久久v久久| 国产男人的电影天堂91| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜福利视频在线观看免费| 一级毛片电影观看| 亚洲在久久综合| 久久青草综合色| 叶爱在线成人免费视频播放| 久久久久国产精品人妻一区二区| 日韩中字成人| 亚洲第一av免费看| 国产精品女同一区二区软件| 久久久久久久久久人人人人人人| 女性被躁到高潮视频| 一区二区三区精品91| 伊人久久大香线蕉亚洲五| 七月丁香在线播放| 亚洲视频免费观看视频| 热re99久久国产66热| 男女下面插进去视频免费观看| 天堂俺去俺来也www色官网| 精品国产一区二区久久| 国产又爽黄色视频| 国产成人av激情在线播放| 国产高清不卡午夜福利| 狂野欧美激情性bbbbbb| 日韩一本色道免费dvd| 亚洲国产精品成人久久小说| 国产福利在线免费观看视频| 国产在视频线精品| 亚洲成色77777| 毛片一级片免费看久久久久| 老汉色av国产亚洲站长工具| 久久青草综合色| 亚洲在久久综合| 在线观看美女被高潮喷水网站| 秋霞伦理黄片| 大片电影免费在线观看免费| 少妇精品久久久久久久| 成人黄色视频免费在线看| 亚洲精品成人av观看孕妇| 国产成人精品久久二区二区91 | 国产精品三级大全| 国产日韩欧美视频二区| 波野结衣二区三区在线| 精品少妇久久久久久888优播| 亚洲人成网站在线观看播放| 成人国产麻豆网| 丝袜在线中文字幕| 99国产综合亚洲精品| 综合色丁香网| 久久毛片免费看一区二区三区| 一级片'在线观看视频| 精品亚洲乱码少妇综合久久| 看非洲黑人一级黄片| www日本在线高清视频| 在线天堂中文资源库| 日韩成人av中文字幕在线观看| 亚洲三级黄色毛片| av片东京热男人的天堂| 热re99久久精品国产66热6| 国产福利在线免费观看视频| 久久韩国三级中文字幕| 国产伦理片在线播放av一区| 亚洲伊人久久精品综合| 伊人久久大香线蕉亚洲五| 黄片无遮挡物在线观看| 热re99久久精品国产66热6| 亚洲精品视频女| 人成视频在线观看免费观看| 久久久久久久精品精品| 亚洲国产日韩一区二区| 2018国产大陆天天弄谢| 国产成人aa在线观看| videosex国产| 高清不卡的av网站| 午夜福利乱码中文字幕| 国产精品成人在线| 下体分泌物呈黄色| 成人毛片60女人毛片免费| 9热在线视频观看99| 纵有疾风起免费观看全集完整版| 老司机亚洲免费影院| 久久国产精品大桥未久av| 啦啦啦在线免费观看视频4| 在线天堂最新版资源| 精品一区二区三卡| 亚洲精品国产av成人精品| 精品亚洲成国产av| 在线观看免费日韩欧美大片| 午夜激情久久久久久久| 国产精品成人在线| 久久ye,这里只有精品| av在线app专区| 视频区图区小说| 久久精品久久久久久久性| 欧美精品一区二区大全| 国产免费视频播放在线视频| 国产伦理片在线播放av一区| 精品午夜福利在线看| 精品国产一区二区三区四区第35| 最新中文字幕久久久久| 亚洲三区欧美一区| 久久99精品国语久久久| 水蜜桃什么品种好| 亚洲欧美一区二区三区黑人 | 久久久精品94久久精品| 国产免费又黄又爽又色| 桃花免费在线播放| 国产精品久久久久久久久免| 亚洲精品久久久久久婷婷小说| av女优亚洲男人天堂| 男女免费视频国产| 性色avwww在线观看| 国产高清不卡午夜福利| 美女福利国产在线| 久久久久国产网址| 人人妻人人添人人爽欧美一区卜| 午夜福利,免费看| 两个人看的免费小视频| 在线观看免费视频网站a站| 国产成人精品婷婷| 在线观看一区二区三区激情| 99热网站在线观看| 亚洲图色成人| 精品久久久久久电影网| 久久精品国产亚洲av天美| 五月伊人婷婷丁香| 成年av动漫网址| 亚洲av成人精品一二三区| 熟妇人妻不卡中文字幕| 亚洲伊人久久精品综合| 99久久人妻综合|