• <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.

    午夜久久久在线观看| 丝袜脚勾引网站| 99热全是精品| 1024香蕉在线观看| 桃红色精品国产亚洲av| 国产男女内射视频| 又黄又粗又硬又大视频| 午夜免费成人在线视频| 中文字幕人妻丝袜一区二区| 国产熟女午夜一区二区三区| 9色porny在线观看| 欧美精品啪啪一区二区三区 | 99精品欧美一区二区三区四区| 美女大奶头黄色视频| 爱豆传媒免费全集在线观看| 老司机影院成人| 在线观看人妻少妇| 亚洲中文av在线| 麻豆av在线久日| 2018国产大陆天天弄谢| 水蜜桃什么品种好| 亚洲精品一二三| 母亲3免费完整高清在线观看| 亚洲精品乱久久久久久| 超碰成人久久| 美女中出高潮动态图| 免费观看av网站的网址| 亚洲综合色网址| 中文字幕精品免费在线观看视频| 久久国产精品男人的天堂亚洲| 久久国产精品大桥未久av| bbb黄色大片| 啦啦啦中文免费视频观看日本| 在线观看www视频免费| 他把我摸到了高潮在线观看 | 欧美中文综合在线视频| 免费日韩欧美在线观看| 国产成人系列免费观看| 宅男免费午夜| 欧美成狂野欧美在线观看| www.熟女人妻精品国产| 亚洲精品国产区一区二| 久久久水蜜桃国产精品网| 香蕉国产在线看| 国产精品一区二区免费欧美 | 国产国语露脸激情在线看| av国产精品久久久久影院| 999精品在线视频| 国产精品二区激情视频| bbb黄色大片| 国产免费视频播放在线视频| 亚洲少妇的诱惑av| 在线看a的网站| 777米奇影视久久| 欧美国产精品va在线观看不卡| 极品少妇高潮喷水抽搐| 久久中文字幕一级| 日本猛色少妇xxxxx猛交久久| 无限看片的www在线观看| 精品亚洲成a人片在线观看| 精品国产一区二区三区久久久樱花| 国产精品秋霞免费鲁丝片| 91精品国产国语对白视频| 三上悠亚av全集在线观看| 19禁男女啪啪无遮挡网站| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲中文日韩欧美视频| 欧美av亚洲av综合av国产av| 亚洲av日韩在线播放| 777米奇影视久久| 成年人午夜在线观看视频| 美女午夜性视频免费| 久久精品国产亚洲av高清一级| 人妻久久中文字幕网| 一边摸一边抽搐一进一出视频| 久久国产精品男人的天堂亚洲| 国产福利在线免费观看视频| 国产一区二区 视频在线| 激情视频va一区二区三区| 人妻一区二区av| 两性夫妻黄色片| 青青草视频在线视频观看| 另类精品久久| 色94色欧美一区二区| 日本撒尿小便嘘嘘汇集6| 人妻人人澡人人爽人人| 一区二区日韩欧美中文字幕| 国产欧美日韩一区二区精品| 搡老熟女国产l中国老女人| 国产亚洲av片在线观看秒播厂| 国内毛片毛片毛片毛片毛片| 高清av免费在线| 黄色视频在线播放观看不卡| 在线精品无人区一区二区三| 色婷婷av一区二区三区视频| 岛国在线观看网站| 欧美黑人欧美精品刺激| 99热网站在线观看| 国产精品1区2区在线观看. | 纯流量卡能插随身wifi吗| 五月开心婷婷网| 亚洲国产av影院在线观看| 黑人操中国人逼视频| 在线观看人妻少妇| 国产成人影院久久av| 人人澡人人妻人| 黄色 视频免费看| 国产精品.久久久| 国产麻豆69| 一区二区三区乱码不卡18| 亚洲欧美色中文字幕在线| 99国产精品一区二区蜜桃av | 婷婷丁香在线五月| 亚洲中文字幕日韩| 国产免费现黄频在线看| 美女扒开内裤让男人捅视频| 欧美亚洲日本最大视频资源| 欧美激情高清一区二区三区| 久久ye,这里只有精品| 啦啦啦啦在线视频资源| 亚洲性夜色夜夜综合| 日韩中文字幕视频在线看片| 97人妻天天添夜夜摸| a级片在线免费高清观看视频| 亚洲美女黄色视频免费看| 日日爽夜夜爽网站| videos熟女内射| 最近中文字幕2019免费版| 久久久久国产一级毛片高清牌| 欧美日韩精品网址| 19禁男女啪啪无遮挡网站| 色婷婷av一区二区三区视频| 五月天丁香电影| 精品国产一区二区三区四区第35| 超色免费av| 青青草视频在线视频观看| 亚洲全国av大片| 免费观看人在逋| 夫妻午夜视频| 老熟妇仑乱视频hdxx| 叶爱在线成人免费视频播放| 欧美少妇被猛烈插入视频| 国产亚洲一区二区精品| 女人久久www免费人成看片| 最新在线观看一区二区三区| 麻豆乱淫一区二区| 国产精品自产拍在线观看55亚洲 | 日韩制服骚丝袜av| 丝袜脚勾引网站| 亚洲成av片中文字幕在线观看| 免费观看a级毛片全部| 50天的宝宝边吃奶边哭怎么回事| 99热全是精品| 人人妻人人爽人人添夜夜欢视频| 免费看十八禁软件| 777米奇影视久久| 午夜福利视频精品| 97精品久久久久久久久久精品| 18禁观看日本| kizo精华| 亚洲国产精品一区二区三区在线| 别揉我奶头~嗯~啊~动态视频 | 一区二区日韩欧美中文字幕| 国产三级黄色录像| 午夜福利在线观看吧| 777米奇影视久久| 亚洲成人手机| 妹子高潮喷水视频| 美女中出高潮动态图| 亚洲精品中文字幕一二三四区 | 伦理电影免费视频| 亚洲欧美日韩另类电影网站| 深夜精品福利| 一区在线观看完整版| 国产av国产精品国产| 热re99久久精品国产66热6| 免费在线观看完整版高清| 亚洲av欧美aⅴ国产| 久久久久视频综合| 美国免费a级毛片| 亚洲国产毛片av蜜桃av| 女人久久www免费人成看片| 国产欧美日韩综合在线一区二区| 精品卡一卡二卡四卡免费| 自线自在国产av| 亚洲av电影在线进入| 日本一区二区免费在线视频| 久久热在线av| 国产精品99久久99久久久不卡| 一个人免费看片子| 女人久久www免费人成看片| 免费高清在线观看视频在线观看| 法律面前人人平等表现在哪些方面 | 免费在线观看视频国产中文字幕亚洲 | 亚洲精品美女久久久久99蜜臀| 色精品久久人妻99蜜桃| bbb黄色大片| 看免费av毛片| 亚洲欧美日韩另类电影网站| 天天躁狠狠躁夜夜躁狠狠躁| 91av网站免费观看| 99精品久久久久人妻精品| 人人澡人人妻人| 丝袜喷水一区| 国产极品粉嫩免费观看在线| 精品少妇一区二区三区视频日本电影| 亚洲欧美一区二区三区黑人| 91精品国产国语对白视频| 精品一品国产午夜福利视频| 亚洲国产av新网站| 妹子高潮喷水视频| 一级,二级,三级黄色视频| 丝袜美足系列| 捣出白浆h1v1| 成人影院久久| 日韩 欧美 亚洲 中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇 在线观看| 天天添夜夜摸| 国产精品一区二区在线不卡| 国产真人三级小视频在线观看| 国产av又大| 欧美性长视频在线观看| 国产三级黄色录像| 久久人人爽av亚洲精品天堂| 天天影视国产精品| 亚洲伊人久久精品综合| 婷婷成人精品国产| 丝袜人妻中文字幕| 五月天丁香电影| 国产一区二区激情短视频 | 国产一区二区三区综合在线观看| svipshipincom国产片| 好男人电影高清在线观看| 1024视频免费在线观看| 欧美激情久久久久久爽电影 | 国产麻豆69| 秋霞在线观看毛片| 啪啪无遮挡十八禁网站| 国产免费一区二区三区四区乱码| 19禁男女啪啪无遮挡网站| 多毛熟女@视频| 亚洲少妇的诱惑av| 一进一出抽搐动态| 999精品在线视频| 亚洲av欧美aⅴ国产| 午夜日韩欧美国产| 国产欧美亚洲国产| 精品第一国产精品| 国产免费av片在线观看野外av| 女人精品久久久久毛片| 亚洲精品国产区一区二| a在线观看视频网站| 99国产精品一区二区三区| 亚洲精品美女久久av网站| svipshipincom国产片| 精品亚洲成a人片在线观看| 成年人黄色毛片网站| 手机成人av网站| 亚洲国产日韩一区二区| 久久热在线av| 老司机深夜福利视频在线观看 | 亚洲欧美成人综合另类久久久| 91精品国产国语对白视频| 国产精品麻豆人妻色哟哟久久| 交换朋友夫妻互换小说| 超碰成人久久| 精品国产乱码久久久久久男人| 亚洲精品久久午夜乱码| 18禁国产床啪视频网站| 美女扒开内裤让男人捅视频| 国产精品久久久av美女十八| 亚洲av男天堂| 亚洲专区中文字幕在线| 亚洲天堂av无毛| 亚洲欧美精品综合一区二区三区| 国产有黄有色有爽视频| 精品国产国语对白av| 十八禁网站免费在线| √禁漫天堂资源中文www| 嫩草影视91久久| 香蕉国产在线看| 一级a爱视频在线免费观看| 国内毛片毛片毛片毛片毛片| 黑人操中国人逼视频| 国产欧美日韩精品亚洲av| 黑人猛操日本美女一级片| 国产亚洲精品第一综合不卡| 国产一区二区在线观看av| 不卡一级毛片| 国产色视频综合| 狂野欧美激情性bbbbbb| 欧美在线一区亚洲| 在线 av 中文字幕| 久久精品国产亚洲av高清一级| 亚洲av片天天在线观看| 亚洲中文av在线| 天堂8中文在线网| 久热这里只有精品99| 成人国语在线视频| 丝袜喷水一区| 国产成人欧美在线观看 | 国产成人欧美| 一级,二级,三级黄色视频| 我的亚洲天堂| 精品人妻1区二区| 国产精品.久久久| 啦啦啦中文免费视频观看日本| 91九色精品人成在线观看| 黄色 视频免费看| 一本久久精品| 黑人猛操日本美女一级片| 热re99久久精品国产66热6| 国产亚洲精品久久久久5区| 亚洲情色 制服丝袜| 精品亚洲乱码少妇综合久久| 一进一出抽搐动态| 欧美日韩亚洲高清精品| 国产精品偷伦视频观看了| 在线 av 中文字幕| 亚洲七黄色美女视频| 美女福利国产在线| 日韩,欧美,国产一区二区三区| 国产亚洲欧美精品永久| 亚洲第一av免费看| 中文字幕色久视频| av在线老鸭窝| 国产成人欧美在线观看 | 久久久国产一区二区| 亚洲少妇的诱惑av| 色婷婷av一区二区三区视频| 亚洲精品一区蜜桃| 国产麻豆69| 成人av一区二区三区在线看 | 嫩草影视91久久| 成人黄色视频免费在线看| 一级a爱视频在线免费观看| 亚洲精品中文字幕一二三四区 | 国产亚洲精品第一综合不卡| av片东京热男人的天堂| av天堂在线播放| 精品国产乱码久久久久久小说| 国产无遮挡羞羞视频在线观看| 色精品久久人妻99蜜桃| 久久久久国产精品人妻一区二区| 最近最新免费中文字幕在线| 爱豆传媒免费全集在线观看| 欧美精品av麻豆av| 精品少妇内射三级| 搡老熟女国产l中国老女人| 成年人午夜在线观看视频| 人妻人人澡人人爽人人| 激情视频va一区二区三区| 伊人久久大香线蕉亚洲五| 久热这里只有精品99| 国产精品九九99| 亚洲熟女精品中文字幕| 亚洲av电影在线观看一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲成人免费av在线播放| 免费观看人在逋| av片东京热男人的天堂| 免费一级毛片在线播放高清视频 | 国产精品影院久久| 少妇人妻久久综合中文| 国产精品99久久99久久久不卡| 国产xxxxx性猛交| 久久人妻熟女aⅴ| 亚洲av日韩在线播放| 国产精品麻豆人妻色哟哟久久| 免费看十八禁软件| 国产高清videossex| 午夜激情av网站| 国产主播在线观看一区二区| 国产色视频综合| 国产成人精品久久二区二区免费| 人人妻,人人澡人人爽秒播| 性高湖久久久久久久久免费观看| 免费人妻精品一区二区三区视频| 性少妇av在线| 欧美国产精品va在线观看不卡| 一个人免费看片子| 999久久久国产精品视频| 在线观看舔阴道视频| 中文字幕av电影在线播放| 妹子高潮喷水视频| 久久久欧美国产精品| 国产亚洲精品久久久久5区| 国产男人的电影天堂91| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲中文av在线| 首页视频小说图片口味搜索| 高清欧美精品videossex| 亚洲 国产 在线| 精品久久久久久久毛片微露脸 | 亚洲欧美日韩高清在线视频 | 国产亚洲av高清不卡| 国产在线视频一区二区| 亚洲精品av麻豆狂野| 19禁男女啪啪无遮挡网站| 日日摸夜夜添夜夜添小说| 性高湖久久久久久久久免费观看| 黑人操中国人逼视频| 国产av精品麻豆| 人妻人人澡人人爽人人| 满18在线观看网站| 少妇 在线观看| 丝袜人妻中文字幕| 欧美日韩黄片免| 亚洲天堂av无毛| 国产欧美日韩一区二区三区在线| 少妇人妻久久综合中文| 亚洲第一av免费看| 欧美少妇被猛烈插入视频| 人妻 亚洲 视频| 久久ye,这里只有精品| 宅男免费午夜| 日本av手机在线免费观看| 国产欧美日韩精品亚洲av| 亚洲精品av麻豆狂野| 午夜福利免费观看在线| 夜夜骑夜夜射夜夜干| 51午夜福利影视在线观看| 国产精品熟女久久久久浪| 天天操日日干夜夜撸| 水蜜桃什么品种好| 国产色视频综合| 黄色 视频免费看| 国产国语露脸激情在线看| 亚洲精品一区蜜桃| 看免费av毛片| 欧美黄色片欧美黄色片| 国产1区2区3区精品| 90打野战视频偷拍视频| 欧美人与性动交α欧美软件| 美女大奶头黄色视频| 欧美成人午夜精品| 欧美日韩一级在线毛片| 人人妻,人人澡人人爽秒播| 一本—道久久a久久精品蜜桃钙片| 啦啦啦免费观看视频1| 夫妻午夜视频| 久久中文看片网| 色视频在线一区二区三区| 啦啦啦 在线观看视频| 精品卡一卡二卡四卡免费| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久久精品电影小说| 伊人亚洲综合成人网| 一边摸一边做爽爽视频免费| 国产欧美日韩一区二区精品| 日本五十路高清| av又黄又爽大尺度在线免费看| 久久亚洲精品不卡| 青青草视频在线视频观看| 女人久久www免费人成看片| 黑人操中国人逼视频| 亚洲欧美一区二区三区黑人| 亚洲精品成人av观看孕妇| 大香蕉久久成人网| 成人国产一区最新在线观看| 各种免费的搞黄视频| 久久精品aⅴ一区二区三区四区| 欧美老熟妇乱子伦牲交| 精品人妻一区二区三区麻豆| 亚洲五月婷婷丁香| 啦啦啦视频在线资源免费观看| 欧美黑人精品巨大| 老司机午夜十八禁免费视频| 亚洲精品国产区一区二| 人成视频在线观看免费观看| 中文字幕av电影在线播放| 无遮挡黄片免费观看| 日本av免费视频播放| 欧美激情久久久久久爽电影 | 丝袜脚勾引网站| 欧美黑人欧美精品刺激| av福利片在线| 美女国产高潮福利片在线看| 91麻豆精品激情在线观看国产 | 丰满少妇做爰视频| 国产成人精品无人区| 日韩 亚洲 欧美在线| 成年人黄色毛片网站| 亚洲伊人色综图| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区日韩欧美中文字幕| 久久久久久久久免费视频了| 一个人免费在线观看的高清视频 | 黑丝袜美女国产一区| 欧美 亚洲 国产 日韩一| 电影成人av| 亚洲三区欧美一区| 日本一区二区免费在线视频| 亚洲国产毛片av蜜桃av| 国产成人啪精品午夜网站| 欧美老熟妇乱子伦牲交| 国产1区2区3区精品| 国内毛片毛片毛片毛片毛片| 狂野欧美激情性bbbbbb| av福利片在线| 天天影视国产精品| 久久热在线av| 日韩视频一区二区在线观看| 国产精品99久久99久久久不卡| 亚洲熟女精品中文字幕| 国产亚洲av片在线观看秒播厂| 午夜久久久在线观看| 丁香六月天网| 波多野结衣一区麻豆| 超碰97精品在线观看| 久久精品亚洲av国产电影网| 久久久精品区二区三区| 欧美日韩成人在线一区二区| 激情视频va一区二区三区| 在线观看免费日韩欧美大片| 建设人人有责人人尽责人人享有的| 久久天躁狠狠躁夜夜2o2o| 欧美精品亚洲一区二区| 最近最新免费中文字幕在线| 少妇 在线观看| 18在线观看网站| 韩国高清视频一区二区三区| 如日韩欧美国产精品一区二区三区| 超碰97精品在线观看| 一区在线观看完整版| 每晚都被弄得嗷嗷叫到高潮| 久久精品久久久久久噜噜老黄| 男人爽女人下面视频在线观看| 国产av精品麻豆| 在线观看舔阴道视频| 久久中文字幕一级| 超碰成人久久| 高清黄色对白视频在线免费看| 亚洲国产精品成人久久小说| 满18在线观看网站| a级毛片黄视频| 永久免费av网站大全| 水蜜桃什么品种好| 一级,二级,三级黄色视频| 一边摸一边做爽爽视频免费| 免费高清在线观看视频在线观看| 天堂俺去俺来也www色官网| 如日韩欧美国产精品一区二区三区| 亚洲国产欧美日韩在线播放| 免费在线观看黄色视频的| 丝瓜视频免费看黄片| 一进一出抽搐动态| 女人高潮潮喷娇喘18禁视频| 最新在线观看一区二区三区| 国产成人系列免费观看| 久久久国产一区二区| 日韩欧美一区二区三区在线观看 | 黄色 视频免费看| 精品国产一区二区久久| 丁香六月天网| 久久久国产一区二区| 一区二区日韩欧美中文字幕| 亚洲欧洲精品一区二区精品久久久| 午夜福利,免费看| 黄色毛片三级朝国网站| 久久中文字幕一级| 又黄又粗又硬又大视频| 51午夜福利影视在线观看| 国产精品一区二区精品视频观看| 久久久久国内视频| 亚洲一区中文字幕在线| 久久久久久久国产电影| 视频区图区小说| 亚洲精品一区蜜桃| 国产一区二区在线观看av| 新久久久久国产一级毛片| 日韩中文字幕视频在线看片| 国产淫语在线视频| 久久久精品94久久精品| 久久狼人影院| 国产精品欧美亚洲77777| 人妻一区二区av| 青青草视频在线视频观看| 精品少妇黑人巨大在线播放| 少妇猛男粗大的猛烈进出视频| 亚洲人成77777在线视频| 欧美性长视频在线观看| 色精品久久人妻99蜜桃| 超碰97精品在线观看| 国产极品粉嫩免费观看在线| 91成人精品电影| 老司机福利观看| 久久精品国产亚洲av香蕉五月 | 亚洲av美国av| 亚洲成人手机| 国产伦人伦偷精品视频| 一二三四在线观看免费中文在| 久久精品亚洲av国产电影网| 亚洲全国av大片| 日本vs欧美在线观看视频| 搡老岳熟女国产| www.熟女人妻精品国产| www.999成人在线观看| 色婷婷久久久亚洲欧美| 国产极品粉嫩免费观看在线| 国产91精品成人一区二区三区 | 亚洲国产中文字幕在线视频| 亚洲av电影在线进入| 另类精品久久| svipshipincom国产片| 91av网站免费观看| 欧美日韩亚洲高清精品| 夫妻午夜视频| 大陆偷拍与自拍| 亚洲精品国产区一区二| 欧美另类一区| 精品国产乱码久久久久久男人| 午夜视频精品福利| 国产片内射在线|