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

    Grain yield,nitrogen use efficiency and physiological performance of indica/japonica hybrid rice in response to various nitrogen rates

    2023-01-06 10:57:28ZHOUQunYUANRuiZHANGWeiyangGUJunfeiLlULijunZHANGHaoWANGZhiqinYANGJianchang
    Journal of Integrative Agriculture 2023年1期

    ZHOU Qun,YUAN Rui,ZHANG Wei-yang,GU Jun-fei,LlU Li-jun,ZHANG Hao,WANG Zhi-qin,YANG Jian-chang

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

    Abstract Utilizing the heterosis of indica/japonica hybrid rice (IJHR) is an effective way to further increase rice grain yield.Rational application of nitrogen (N) fertilizer plays a very important role in using the heterosis of IJHR to achieve its great yield potential.However,the responses of the grain yield and N utilization of IJHR to N application rates and the underlying physiological mechanism remain elusive.The purpose of this study was to clarify these issues.Three rice cultivars currently used in rice production,an IJHR cultivar Yongyou 2640 (YY2640),a japonica cultivar Lianjing 7 (LJ-7) and an indica cultivar Yangdao 6 (YD-6),were grown in the field with six N rates (0,100,200,300,400,and 500 kg ha–1)in 2018 and 2019.The results showed that with the increase in N application rates,the grain yield of each test cultivar increased at first and then decreased,and the highest grain yield was at the N rate of 400 kg ha–1 for YY2640,with a grain yield of 13.4 t ha–1,and at 300 kg ha–1 for LJ-7 and YD-6,with grain yields of 9.4–10.6 t ha–1.The grain yield and N use efficiency (NUE) of YY2640 were higher than those of LJ-7 or YD-6 at the same N rate,especially at the higher N rates.When compared with LJ-7 or YD-6,YY2640 exhibited better physiological traits,including greater root oxidation activity and leaf photosynthetic rate,higher cytokinin content in the roots and leaves,and more remobilization of assimilates from the stem to the grain during grain filling.The results suggest that IJHR could attain both higher grain yield and higher NUE than inbred rice at either low or high N application rates.Improved shoot and root traits of the IJHR contribute to its higher grain yield and NUE,and a higher content of cytokinins in the IJHR plants plays a vital role in their responses to N application rates and also benefits other physiological processes.

    Keywords: indica/japonica hybrid rice,grain yield,nitrogen use efficiency,cytokinins,non-structural carbohydrate,root activity

    1.lntroduction

    Rice (Oryza sativaL.) is one of the most important food crops in the world,providing nearly 60% of the dietary calories consumed by approximately 3 billion people(FAOSTAT 2019).Rice in China is the staple food for more than 65% of the population (Penget al.2009,2010).The rapid population growth and economic development need greater increases in food production,especially rice production (Zhang 2007;Penget al.2009;Ma and Yuan 2015).However,reports indicate a downward trend in the annual increase rate of rice yield in the world,and the rate of increase has even stagnated in some countries,such as in China (Penget al.2010;Grassiniet al.2013;Xuet al.2016).Due to the limitations of cultivation areas and water resources,the potential for increasing the rice planting area to increase the total yield is very limited,and therefore,it is necessary to increase the rice yield by increasing yield per unit area (Weiet al.2016).

    To further increase yield potential,China has successfully bred “super” rice or “super” hybrid rice cultivars in the last 20 years by the use of intersubspecific heterosis combined with improvements in plant types(Chenget al.2007;Penget al.2008,2009;Ma and Yuan 2015).To date,146 “super” rice cultivars,including 92“super” hybrids,have been released with a growing area of more than 40 million hectares (Chenget al.2015;Cheng 2016;Zhang Het al.2021).Compared with non-super rice cultivars,super rice cultivars can increase grain yield by 15–20% in field experiments and by approximately 10%in large-area production (Zhanget al.2009;Ma and Yuan 2015;Huanget al.2016).To obtain a higher grain yield for a hybrid rice cultivar,a high nitrogen (N) rate is usually applied,which results in a low N use efficiency (NUE) (Xueet al.2013;Guet al.2017).Recently,anindica/japonicahybrid rice (IJHR) cultivar,Yongyou 2640 (YY2640),has been released and it has demonstrated a very high grain yield in the lower reaches of the Yangtze River in China(Menget al.2018;Zhuet al.2020).Several agronomic and physiological traits of the IJHR cultivar YY2640 have been documented,such as a larger sink size characterized by a greater number of spikelets per panicle (Huet al.2018;Zhanget al.2018),greater biomass accumulation and leaf area (Weiet al.2016,2017;Huet al.2018),and better root morphological architecture and physiological characteristics (Menget al.2018).Some reports show that a higher grain yield of the IJHR cultivar is usually obtained at a higher N application rate (Weiet al.2017;Huet al.2018;Zhanget al.2018).However,Huanget al.(2018) observed that hybrid rice,especially IJHR cultivars,exhibited higher yields over the inbred rice at low and moderate N rates.Obviously,it remains controversial whether the IJHR cultivars respond to N application rates.

    Nitrogen is an essential nutrient for plant growth.The application of N fertilizer is a primary practice to increase rice grain yield (Liuet al.2015;Zhanget al.2015).However,over application of N fertilizer in rice production is a serious problem in China (Penget al.2010;Yanet al.2018;Zhang W Yet al.2021).Today the national average amount of N used in China’s rice production is about 180 kg ha–1,75% higher than the global average (Penget al.2010;Zhang Het al.2021).In some districts,such as the Taihu Lake area,the N application rate of rice is more than 270 kg ha–1(Grassiniet al.2013;Zhanget al.2018).Excessive application of N fertilizer can not only increase the loss of fertilizer-N,leading to environmental pollution (Hofmeieret al.2015;Lianget al.2017;Valliereet al.2017;Yanet al.2018),but also increase production cost and plant lodging risk and reduce NUE and grain quantity and quality (Dominiciet al.2014;Chenet al.2015;Roweet al.2015;Keskinenet al.2016).Although there are reports showing that many practices,such as using N-efficient rice cultivars (Zhuet al.2020),adoption of site-specific N management (Penget al.2010) or a higher proportion of N application at the late growth stage(Zhanget al.2018),and applying controlled/slow release N-fertilizer (Sunet al.2019;Zhang W Yet al.2021),can decrease N loss and increase NUE,little is known about whether and how a proper N application rate can synergistically increase grain yield and NUE in IJHR.

    The objectives of this study were to: (1) compare the differences in the responses of yield formation to N application rates between IJHR and inbred rice cultivars,(2) understand the agronomic and physiological performances of IJHR under different N application rates,and (3) explore a rational N application rate to achieve both high grain yield and high NUE for IJHR.Some agronomic and physiological traits,including N uptake,leaf area duration (LAD),crop growth rate (CGR),leaf photosynthetic rate and net assimilation rate (NAR),non-structural carbohydrate (NSC) accumulation and remobilization in the stem,root biomass,root oxidation activity (ROA),and cytokinins (zeatin+zeatin riboside,Z+ZR) in roots were determined.This study will provide new insight into the mechanism underlying the response of IJHR to N application rates and synergistically increasing grain yield and NUE in rice.

    2.Materials and methods

    2.1.Plant materials and growing conditions

    The field experiment was conducted at a research farm of Yangzhou University,Jiangsu Province,China (32°30′N,119°25′E) during the rice growing seasons (May–October)of 2018 and 2019.The soil in the field was a sandy loam(Typic Fluvaquent,Etisols) with 24.2 g kg–1organic matter,102 mg kg–1alkali hydrolysable N,34.1 mg kg–1Olsen-P,and 67.4 mg kg–1exchangeable K.The soil moisture content was 0.191 g g–1at field capacity.The soil bulk density was 1.30 g cm–3.The precipitation,solar radiation and mean temperature during the two rice growing seasons are shown in Appendix A.

    Three rice (Oryza sativaL.) cultivars qualified by the Ministry of Agriculture and Rural Affair of China,Yongyou 2640(YY2640,anindica/japonicahybrid rice,IJHR),Lianjing 7(LJ-7,ajaponicainbred rice) and Yangdao 6 (YD-6,anindicainbred rice),were grown in the paddy field.YY2640 is a typical IJHR cultivar with greater heterosis in grain yield,stronger lodging-resistance,and wider adaptability (Menget al.2018;Zhuet al.2020;Liuet al.2021).Both inbred rice cultivars of LJ-7 and YD-6 are currently used in local production with large areas for their high yield and good quality (Wanget al.2016;Zhanget al.2020).Across the two study years,seeds were raised in the seedbed with a sowing date of 12 May and seedlings were transplanted on 11 June with a hill spacing of 15 cm×25 cm and two seedlings per hill.The three cultivars had similar growth durations from sowing to maturity,ranging from 155 to 157 d.The heading dates(50% of plants headed) of the three cultivars were on 13 August for YY2640 and 16 August for LJ-7 and YD-6 in 2018,and on 15 August for YY2640 and 18 August for LJ-7 and YD-6 in 2019.All plants were harvested from 15 to 17 October in both years.

    2.2.Experimental design and treatments

    The experiment was a 3×6 (three rice cultivars and six N application rates) factorial design with 18 treatment combinations.Each treatment had three plots as repetitions in a complete randomized block design.Plot dimensions were 4.5 m×5 m and plots were separated by a 1-m wide alley with plastic film inserted into the soil to a depth of 0.50 m to form a barrier.The six N rates included zero-N control (CK),100,200,300,400 and 500 kg ha–1.Urea as the N fertilizer was applied at the stages of pre-transplanting (1 day before transplanting),early tillering (7 days after transplanting),panicle initiation (the first appearance of differentiated apex),and spikelet initial differentiation (the appearance of glumous flower primordia at the tips of elongating primary rachis–branches) at a ratio of 4:2:2:2.Phosphorus(240 kg ha–1as single superphosphate) was applied before transplanting as basal fertilizer.Potassium(180 kg ha–1as KCl) was applied at the stages of pretransplanting and jointing at a ratio of 6:4.Weeds,insects,and diseases were controlled by either chemical or manual methods to avoid yield losses.Irrigation canals and ditches were made between the plots.Except for mid-season drainage,the field was kept at 2–3 cm water depth throughout the growth season until 1 week before the final harvest.

    2.3.Sampling and measurements

    The plants of five hills were sampled from each plot at the middle tillering (20 days after transplanting),panicle initiation (40–43 days after transplanting),the heading time (63–66 days after transplanting),and at maturity(125–128 days after transplanting) for measuring the shoot biomass.The sampled plants were separated into three components of leaves,stems,and panicles (the panicles were only determined at heading and maturity),and these separated samples were dried at 80°C for 72 h and then weighed.The contents of non-structural carbohydrate (NSC) in the stems and sheaths at heading and maturity were measured according to the method of Yoshidaet al.(1976).The sample was dried in an oven and ground into a fine powder.In a 15-mL centrifuge tube,100 mg of ground sample was added with 10 mL of ethanol (density 630 g L–1) and kept in a water bath at 80°C for 30 min.The tube was then centrifuged at 2 000 r min–1for 20 min after cooling.The supernatant was collected and the extraction was repeated three times.The alcohol in the supernatant was evaporated on a water bath at 80°C until most of the alcohol was removed and the volume was reduced to about 3 mL.The sugar extract was then diluted to 25 mL with distilled water.The concentration of sugars in the extract was then analyzed as described by Yoshidaet al.(1976).

    The residue left in the centrifuge tube was dried at 80°C for starch extraction.A total of 2 mL of distilled water was added to the tube containing the dried residue.The tube was then shaken in a boiling water bath for 15 min.A total of 2 mL of 9.36 mol L–1HClO4was added to the tube after cooling.The solution was shaken for 15 min.The extract was then made up to about 10 mL and centrifuged at 2 000 r min–1for 20 min.The supernatant was collected and a further 2 mL of 4.68 mol L–1HClO4was added to the residue.The extraction was repeated as above.The supernatants were combined and made up to 50 mL with distilled water.The starch was analyzed by the method of Yoshidaet al.(1976).Leaf area was measured by using a leaf area meter (LI-3000C,Li-Cor;Lincoln,NE,USA).The crop growth rate (CGR),leaf area duration (LAD),and net assimilation rate (NAR) were calculated using the following formulas:

    whereW1andW2andLAI1andLAI2indicate the shoot biomass (dry weight) and leaf area measured at the first (t1)and the second measurement times (t2),respectively.

    The net photosynthetic rate of leaves and cytokinins(zeatin+zeatin riboside,Z+ZR) in leaves were measured at panicle initiation,heading time,and at the stages of early grain filling (11–12 days after heading (DAH)),mid grain filling (23–24 DAH),and late grain filling (35–36 DAH).The measurement of leaf photosynthetic rate was made by using a gas exchange analyzer (Li-Cor 6400 portable photosynthesis measurement system,Li-Cor,Lincoln,NE,USA) at 9:00–11:00 h in the morning,when the photosynthetically active radiation above the canopy was 1 300–1 500 μmol m–2s–1,with six leaves for each treatment.The extraction and determination of Z+ZR in leaves were made according to the protocol of Panet al.(2010).Briefly,the detection of Z+ZR in plants was made by multiple reaction monitoring (MRM) liquid chromatography-tandem mass spectrometry (LC-MS/MS)(TSQ Vantage,Thermo Fisher Scientific,Waltham,MA,USA),and was quantified using a calibration curve with known amounts of standard and based on the ratio of the summed areas of the MRM transitions for Z+ZR to those for their internal standards.Data acquisition and analysis were performed using Xcalibur Data System (Thermo Fisher Scientific,Waltham,MA,USA) (Zhanget al.2018).

    At the aforementioned growth stages,a sampling core was used for root sampling,which removed a 25 cm×15 cm×20 cm (length×width×depth) cube of soil around each individual hill,from two hills in each plot.Roots from three hills in each treatment were used for measuring the root length and root weight,and the remainder of the roots were used for determining root oxidation activity (ROA) and the Z+ZR content in roots.The content of Z+ZR in roots was determined following the method of Panet al.(2010) as described above.The ROA was determined by measuring the oxidation of α-naphthylamine (α-NA).A total of 2 g of fresh roots were transferred into a 150-mL flask containing 50 mL of 20 mg L–1α-NA.The flask was incubated for 2 h at 25°C in an end-over-end shaker.After incubation,the aliquots were filtered through a quantitative filter paper and 2 mL was mixed with 1 mL NaNO3(1.18 mmol L–1) and 1 mL sulphanilic acid,and the resulting color was measured by a spectrophotometer (Lambda 35 UV/VIS Spectrometer,PerkinElmer Life and Analytical Sciences,Inc.,Waltham,MA).ROA was expressed as μg α-NA per gram dry weight(DW) per hour (μg α-NA g–1DW h–1) (Zhanget al.2009).

    All plants were hand harvested on 15–17 October 2018 and 2019 for measuring grain yield and yield components.Grain yield was determined from a harvest area of 5.0 m2in each plot and was adjusted to 14%moisture.Aboveground biomass and yield components(i.e.,the number of panicles per square meter,the number of spikelets per panicle,the percentage of filled grains,and the grain weight) were determined from the plants of 10 hills (excluding the border hills) and sampled randomly from each plot.The percentage of filled grains was defined as the filled grains (specific gravity≥1.06 g cm-3) as a percentage of the total number of spikelets.The N content in plants was determined by micro Kjeldahl digestion,distillation,and titration to calculate the aboveground N uptake.The methods for calculating NUE were based on those of Xueet al.(2013).

    2.4.Statistical analysis

    Analysis of variance was conducted using the Statistical Analysis Software (SPSS,23.0 Software,Inc.,Chicago,IL,USA).Year and replication were included as random effects.N rates,cultivars,and their interactions were included as fixed effects.Data from each sampling date were analyzed separately.Means were tested by least significant difference atP=0.05 (LSD0.05).R Software(Corrplot,version 3.5.1,https://cran.r-project.org) was used to calculate the Pearson correlation coefficients and to perform the graphing.

    3.Results

    3.1.Differences in experimental factors

    ComputedFvalues for the differences in grain yield,N uptake,internal N use efficiency,partial factor productivity of applied N,agronomic N use efficiency and the main agronomic and physiological traits of shoots and roots in rice between or among years,cultivars and N rates are shown in Table 1.All the measurements showed significant differences (P<0.05) among the cultivars,N rates and the interactions between the cultivar and N rate (Table 1).Variations due to years and in interactions between year and cultivar as well as year and N rate were not significant (P≥0.05) (Table 1).Since year was not a significant factor in the experiment,the data were averaged across the two study years.The data for 2018 and 2019 are separately presented in Appendices B–O.

    3.2.Grain yield and its components

    With the increase in N application rates,the grain yield of each test cultivar first increased and then decreased,and the highest grain yields obtained were at the N rate of 400 kg ha–1for YY2640 with a grain yield of 13.4 t ha–1and 300 kg ha–1for both check cultivars (LJ-7 or YD-6) with grain yields of 9.4–10.6 t ha–1(Table 2).At the same N rate,especially at a higher N rate,the grain yield of YY2640 was higher than that of either LJ-7 or YD-6.The number of spikelets per panicle and the total spikelets per unit area increased with the increase in Napplication rates of YY2640,but they first increased and then decreased with the increase in N application rates for LJ-7 and YD-6 (Table 2).At the same N application rate,the number of spikelets per panicle and the total spikelets per unit area were significantly higher for YY2640 than for LJ-7 or YD-6.With the increase in N application rates,the percentage of filled grains and 1 000-grain weight of each test cultivar decreased (Table 2).In addition,no lodging was observed for any of the test cultivars at any of the N application rates.

    Table 1 Analysis of variance of grain yield,nitrogen use efficiency (NUE),and the main agronomic and physiological traits of rice between or among years,cultivars,and N rates

    Table 2 Grain yield and its components under different nitrogen rates

    3.3.Plant N uptake,N accumulation,and NUE

    Changes in plant N uptake at different growth stages of the three test cultivars showed a similar trend,and the N uptake increased with the increase in N application rates and with the plant growth process (Fig.1-A–D).Under the same N rate and at the same growth stage,the plant N uptake of YY2640 was higher than that of either LJ-7 or YD-6 (Fig.1-A–D).The plant N accumulation of the three cultivars also increased with the increase in N application rates,but varied with growth stages (Fig.1-E–G).It first increased and then decreased with the growth process,and reached the maximum from panicle initiation to heading.Generally,the IJHR cultivar YY2640 accumulated more N in the plants than thejaponicainbred LJ-7 orindicainbred YD-6 at the same N rate and during the same growth period,especially from heading to maturity (Fig.1-E–G).

    Fig.1 Differences in nitrogen (N) uptake (A–D) and N accumulation (E–G) during the main growth stages of rice under various N application rates.MT,PI,HD and MA represent mid-tillering,panicle initiation,heading,and maturity,respectively.Vertical bars represent standard deviations of the mean where these exceed the size of the symbol (n=6).Data are presented as the averages of the two study years.Different letters above the columns indicate statistical significance at the P=0.05 level within the same measurement date.

    With the increase in N application rates,internal N use efficiency (IEN),partial factor productivity of applied N (PFPN),physiological N use efficiency (PEN),recovery efficiency of N fertilizer (REN),agronomic N use efficiency(AEN) and N harvest index (HIN) all decreased for all of the test cultivars (Table 3).There were no significant differences inAENorPENbetween YY2640 and LJ-7 at the N rates of 100 and 200 kg ha–1,while YY2640 showed greaterAENandPENthan LJ-7 at the N rates of 300,400 and 500 kg ha–1.At the same N rate,YY2640 always showed higherIEN,PFPN,PEN,REN,AEN,andHINwhen compared with YD-6,and higherIEN,PFPN,REN,andHINin comparison with LJ-7 (Table 3).

    3.4.Shoot physiological traits

    The total leaf area index (LAI,the LAI of total tillers and main stems) at heading,the LAI of productive tillers and main stems at heading,and LAI at maturity increased with the increase in N application rates (Table 4).Values of LAI were higher for YY2640 than for LJ-7 or YD-6 at the same N rate.Similar observations were obtained for the leaf area duration (LAD) from mid-tillering to panicle initiation,from panicle initiation to heading and from heading to maturity (Table 4).

    Table 3 Nitrogen uptake and N use efficiency under different N rates1)

    Table 4 Leaf area index (LAI) and leaf area duration (LAD) of rice under different nitrogen rates1)

    Table 5 Grain-leaf ratios of rice under different nitrogen rates1)

    The grain-leaf ratios,in terms of spikelets/leaf area,filled grains/leaf area and grain yield/leaf area of YY2640 were significantly higher than those of LJ-7 or YD-6 under the same N rate (Table 5).The spikelets/leaf area and filled grains/leaf area of YY2640 and the spikelets/leaf area of LJ-7and YD-6 first increased and then decreased with the increase in N application rates,and reached their maxima at the N rate of 200 kg ha–1for YY2640 and at 100 kg ha–1for both check cultivars (LJ-7 and YD-6).The grain yield/leaf area of YY2640 and the filled grains/leaf area and grain yield/leaf area of LJ-7 and YD-6 decreased with the increase in N application rates (Table 5).

    At the same N rate,shoot biomass was significantly higher for YY2640 than for LJ-7 or YD-6 (Fig.2-A–D).With the increase in N application rates,the shoot biomasses of the three cultivars also increased.Similar changing patterns were observed on both crop growth rate (CGR) and net assimilation rate (NAR) (Fig.2-E–J),implying that higher CGR and NAR contribute to a greater biomass for the IJHR cultivar.

    Fig.2 Shoot biomass (A–D),crop growth rate (E–G),and net assimilation rate (H–J) of rice under various nitrogen rates.MT,PI,HD and MA represent mid-tillering,panicle initiation,heading,and maturity,respectively.Vertical bars represent standard deviations of the mean where these exceed the size of the symbol.Data are presented as the averages of the two study years.Different letters above the columns indicate statistical significance at the P=0.05 level within the same measurement date.

    The photosynthetic rates of the three cultivars significantly increased with the increase in N application rates from panicle initiation to late grain filling,but they showed no significant difference under the rates of 300,400 and 500 kg ha–1(Table 6).YY2640 showed a higher photosynthetic rate than LJ-7 or YD-6 at the same N rate,suggesting that the IJHR cultivar has a higher photosynthetic ability than the inbred rice cultivars at the middle and late growth stages,especially under high N rates.

    Table 6 Leaf photosynthetic rate (μmol CO2 m–2 s–1) of rice under different nitrogen rates

    When the N application rate was lower than 400 kg ha–1,the specific leaf N contents of the three cultivars increased with the increase in N application rates.When the N application rate exceeded 400 kg ha–1,increasing N application rates did not significantly increase specific leaf N contents of any of the test cultivars (Table 7).At the same N rate,the specific leaf N content of YY2640 was significantly higher than that of LJ-7,and showed no significant difference from that of YD-6 at panicle initiation,but was lower than that of YD-6 at the heading time (Table 7).The photosynthetic NUE of the three cultivars decreased with the increase in N application rates at panicle initiation and heading time (Table 7).At the same N rate,photosynthetic NUE of YY2640 was significantly higher than that of YD-6 at panicle initiation,and exhibited no significant difference from that of LJ-7 at panicle initiation,but was higher than those of both check cultivars at the heading time (Table 7).

    When the N application rate was less than 400 kg ha–1for YY2640 and less than 300 kg ha–1for LJ-7 and YD-6,the contents of zeatin (Z)+zeatin riboside (ZR) in leaves increased with the increase in N application rates(Fig.3-A–E).When the N application rate exceeded 400 kg ha–1for YY2640 and was more than 300 kg ha–1for LJ-7 and YD-6,an increase in N application did not significantly decrease the Z+ZR content in leaves of YY2640,while it significantly decreased the Z+ZR content in leaves of both check cultivars.Under the same N application rate,YY2640 produced higher Z+ZR content in leaves than either LJ-7 or YD-6 (Fig.3-A–E).

    Fig.3 Zeatin+zeatin riboside (Z+ZR) contents in the leaves of rice under various nitrogen rates.PI,HD,EF,MF and LF represent panicle initiation,heading,early grain filling,mid grain filling,and late grain filling,respectively.Vertical bars represent ±standard deviations of the mean where these exceed the size of the symbol (n=6).Data are presented as the averages of the two study years.Different letters above the columns indicate statistical significance at the P=0.05 level within the same measurement date.

    In general,NSC accumulation in the stems (culms and sheaths) at the heading time and maturity,and the amount of NSC translocation during grain filling,increased at first and then were decreased with the increase in N application rates,and reached the maxima at the N rate of 400 kg ha–1for YY2640 and at 300 kg ha–1for LJ-7 and YD-6 (Table 8).However,the NSC remobilized from the stems and the NSC contributed to grain yield during the grain-filling period both decreased with increases in N application rates for all the test cultivars.At the same N application rate,YY2640 exhibited higher NSC accumulation in the stems,more remobilization of NSC from stems to grains,and a greater contribution of NSC to grain yield than both check cultivars (Table 8).

    Table 7 Photosynthetic nitrogen use efficiency (NUE) of rice under different N rates1)

    Table 8 Non-structural carbohydrate (NSC) in stems and sheaths,and NSC remobilization of rice under different nitrogen rates1)

    3.5.Root morpho-physiological traits

    Overall,the root length and root biomass increased with the increase in N application rates (Fig.4-A–H).At the same N rate,YY2640 had greater root length and root biomass than LJ-7 or YD-6 at each growth stage.However,the ratio of root to shoot decreased with the increase in N application rates for all the test cultivars at each growth stage,and showed no significant differences among the three cultivars (Fig.4-I–L),suggesting that agreater shoot biomass is closely associated with a greater root biomass for the IJHR.

    Fig.4 Root length (A–D),root biomass (E–H),and the ratio of root to shoot (I–L) of rice under various nitrogen rates.MT,PI,HD and MA represent mid-tillering,panicle initiation,heading,and maturity,respectively.Vertical bars represent standard deviations of the mean where these exceed the size of the symbol (n=6).Data are presented as the averages of the two study years.Different letters above the columns indicate statistical significance at the P=0.05 level within the same measurement date.

    Consistent with the root biomass,the root oxidation activity (ROA) and Z+ZR content in the roots increased with increases in N application rates for the three test cultivars at each growth stage (Fig.5-A–H).At the same N rate,YY2640 had higher ROA and Z+ZR content in roots than LJ-7 or YD-6,especially at the heading time and during grain filling,indicating the stronger root activity of the IJHR cultivar.

    Fig.5 Root oxidation activity (A–D) and zeatin+zeatin riboside (Z+ZR) contents in roots (E–H) of rice under various nitrogen rates.MT,PI,HD,EF and MF represent mid-tillering,panicle initiation,heading,early grain filling and mid grain filling,respectively.Vertical bars represent standard deviations of the mean where these exceed the size of the symbol (n=6).Data are presented as the averages of the two study years.Different letters above the columns indicate statistical significance at the P=0.05 level within the same measurement date.

    3.6.Correlations of shoot and root traits with grain yield and NUE

    Pearson correlation analysis showed that high-effective LAI,effective LAI and total LAI at heading time,root and shoot biomass from heading to maturity,CGR and NAR before heading time,leaf photosynthetic rate during grain filling,root length,and ROA and Z+ZR content in roots and leaves during panicle initiation and grain filling were all significantly and positively correlated with the total spikelet number and grain yield (Fig.6).The ROA during panicle initiation and grain filling was significantly and positively correlated with the contents of Z+ZR in both roots and leaves.Contents of Z+ZR in leaves and roots,ROA,root length,LAD,CGR,NAR and leaf photosynthetic rate during grain filling were significantly and positively correlated with root and shoot biomass,grain yield,recovery efficiency of N fertilizer (REN) and agronomic N use efficiency (AEN).The photosynthetic NUE values were significantly and positively correlated with filled grain percentage,filled grains/leaf area,grain yield/leaf area,remobilized NSC reserve,and internal nitrogen use efficiency (IEN) (Fig.6).

    4.Discussion

    It is generally believed that IJHR has a greater capacity for biomass production and greater yield potential when compared with eitherindicaorjaponicainbred rice,and it usually absorbs more N and has a lower NUE than inbred rice (Weiet al.2016,2017;Guet al.2017).The results herein showed that the IJHR cultivar YY2640 exhibited both higher grain yield and higher NUE than either thejaponicainbred or theindicainbred cultivar when the N application rate was the same,especially at a higher N rate (Tables 2 and 3).These results demonstrate that the IJHR cultivar could not only exhibit greater yield potential,but also efficiently use absorbed N to produce more grain yield,leading to a higher NUE,such asIEN.

    Prior to this study,the information about the synergistic improvement in grain yield and NUE in IJHR was rather limited.Our observations in the present study provide several possible explanations for the synergistic improvement.Firstly,the IJHR cultivar YY2640 exhibited a larger sink size (total number of spikelets per unit area)and a higher ratio of sink size to leaf area (spikelets/leaf area) (Tables 2 and 5).A larger sink size is necessary to achieve a higher grain yield (Katoet al.2007;Penget al.2009),while a higher ratio of sink size to leaf area could make plants expend more N to produce productive organs(spikelets) and less N to produce vegetative organs(leaves),leading to higher NUE (Xueet al.2013;Zhang Het al.2021).We observed that the IJHR cultivar absorbed more N during panicle differentiation and development,which may contribute to its larger sink size and higher ratio of sink size to leaf area (Fig.1).

    Secondly,the IJHR cultivar YY2640 showed stronger shoot activity,which was reflected by its larger LAI and LAD,and higher CGR,NAR,and leaf photosynthetic rate than the inbred cultivars at the same N application rate(Fig.2;Tables 4 and 6).Larger LAI and LAD are very important for producing more biomass and increasing sink size (Liuet al.2014;Huanget al.2019;Zhanget al.2019;Zhang Het al.2021).Higher CGR,NAR,and leaf photosynthetic rate can invest less N in producing more biomass and grain yield (Horieet al.2005;Fageria 2007;Xueet al.2013;Zhang Het al.2021),which is also evidenced by our observations that the IJHR cultivar had higher photosynthetic NUE (Table 7).It is proposed that a higher photosynthetic NUE can invest less N in nonphotosynthetic components and more N in dry mater production,which can contribute to a higher NUE (Juet al.2015;Wanget al.2016),and photosynthetic NUE was observed in the present study to be significantly and positively correlated withIEN(Fig.6).

    Thirdly,the IJHR cultivar accumulated more NSC in stems including sheaths before heading and remobilized more NSC from the stem to the grain during grain filling than the inbred rice cultivars (Table 8).Some reports show that more accumulation of NSC in rice stems before heading could enhance sink activity,such as increases in the activities of key enzymes involved in sucrose-to-starch conversion in the grain at the early grain-filling stage,and more NSC remobilization from the stem to the grain during ripening could provide more carbohydrates for grain filling and increase the harvest index and N harvest index,all of which could result in the improvement of grain filling efficiency and NUE (Fuet al.2011;Zhanget al.2013;Okamuraet al.2018;Zhang Het al.2021).

    Fourthly,the IJHR had better root morpho-physiological traits,which were evidenced by its greater root length and biomass and higher ROA in comparison with inbred rice(Figs.4 and 5).A greater root length and root biomass could allow the plants to absorb more water and nutrients from the soil and thereby ensure shoot growth (Zhanget al.2009;Yanget al.2012;Menget al.2018).A higher ROA has been recognized as the most important physiological trait and is vitally important to maintaining root growth and achieving higher shoot activity (Osakiet al.1997;Yanget al.2012;Chuet al.2014;Menget al.2018).Root length,root biomass,and ROA were herein found to be closely associated with grain yield and NUE(Fig.6).

    Fig.6 Pearson correlations of agronomic and physiological traits with grain yield,yield components,and nitrogen use efficiency of rice.r1,grain yield;r2,total spikelets;r3,filled grain;r4,internal N use efficiency;r5,recovery efficiency of N fertilizer;r6,agronomic N use efficiency;r7,partial factor productivity of applied N;r8,total leaf area index (LAI);r9,high-effective LAI;r10,effective LAI;r11,no.of total spikelets/leaf;r12,no.of filled grains/leaf;r13,grain weight/leaf;r14,remobilized non-structural carbohydrate (NSC)reserve;r15,NSC contribution to grain;r16,shoot dry weight at heading time (HT);r17,root dry weight at HT;r18,root:shoot ratio at HT;r19,root length at HT;r20,specific leaf N content (SLNC) at HT;r21,leaf photosynthetic rate at HT;r22,photosynthetic N use efficiency (PNUE) at HT;r23,shoot dry weight at maturity (MA);r24,shoot biomass from HT to MA;r25,root dry weight at MA;r26,root:shoot ratio at MA;r27,LAD;r28,CGR;r29,net assimilation rare (NAR);r30,leaf photosynthetic rate during grain filling;r31,zeatin+zeatin riboside (Z+ZR) contents in leaves at panicle initiation (PI);r32,Z+ZR contents in roots at PI;r33,Z+ZR contents in leaves during grain filling;r34,Z+ZR contents in roots during grain filling;r35,root oxidation activity (ROA) at PI;r36,ROA at HT;r37,ROA during grain filling.The red and blue circles indicate negative or positive correlations between parameters,respectively.The darker the color,the higher the correlation.*,** and *** indicate the significant differences at the P<0.05,P<0.01 and P<0.001 levels,respectively.For interpretation of the references to colour in this figure legend,the reader is referred to the web version of this article.

    Furthermore,the IJHR cultivar exhibited higher cytokinin (Z+ZR) contents in both roots and leaves(Figs.3 and 5).A higher Z+ZR content in plants during panicle differentiation and development could promote spikelet differentiation and increase spikelet number per panicle,and the Z+ZR content in rice caryopses during filling is significantly and positively correlated with the endosperm cell proliferation rate and grain-filling rate,and consequently,results in the increases in sink size and activity and grain yield (Yanget al.2000,2002;Kyozuka 2007;Xuet al.2015).It is proposed that there is a tight linkage between the N signaling network and cytokinin metabolism and signaling (Guet al.2018).Greater N uptake or assimilation can enhance cytokinin biosynthesis,and in turn,the increase in cytokinin biosynthesis can enhance N uptake and assimilation in plants and make the absorbed N produce more biomass,and accordingly,lead to a higher NUE (Yanget al.2012;Guet al.2018).

    It is noteworthy that although a maximum grain yield was obtained at the N rate of 400 kg ha–1for the IJHR cultivar,the grain yield showed no significant difference among the N rates of 300,400,and 500 kg ha–1(Table 2).Furthermore,each index of NUE for the IJHR cultivar was significantly lower at the N rate (kg ha–1) of 400 or 500 that at the N rate at 300 (Table 3).Therefore,we suggest that the N application rate of 300 kg ha–1when the soil fertility is medium,such as the soil of the present experiment,or the specific leaf N content (g m–2) is approximately 3.4 and 2.9,respectively,at the panicle initiation and the heading time (Table 7),could be considered as a rational N application index to achieve both high grain yield and high NUE for the IJHR cultivar.

    It should be also mentioned that the total growth durations of the three tested cultivars were around 155–157 d.Similar observations were reported by Zhanget al.(2020) and Liuet al.(2021).Generally,the growth duration of mid-season rice cultivars,such as Yangdao 6,in the lower reaches of the Yangtze River in China is approximately 145–150 d.The slightly longer growth durations of the three tested cultivars in the present study may mainly be attributed to the increase in temperature during the late grain-filling period in recent years (Wanget al.2021),which results in slower plant senescence and a longer grain-filling period.

    To reduce the influence of non-treatment factors on the experimental results,the present study design used two seedlings per hill for all the tested cultivars at transplanting.Usually,more seedlings per hill are recommended for inbred rice than for hybrid rice because the tillering ability of the former is lower than that of the latter.Our results demonstrated,however,that panicle number per unit area was much greater for the inbred rice than for the IJHR at the same N rate (Table 2),indicating that both inbred rice cultivars have stronger tillering ability than the hybrid rice cultivar.Further investigation is merited to determine whether more than two seedlings per hill at transplanting could increase grain yield for inbred rice at N application rates of 200 or 300 kg ha–1.

    5.Conclusion

    The IJHR cultivar YY2640 produced much higher grain yield than inbred rice cultivars at the same N application rate,due mainly to its greater total number of spikelets per unit area which resulted from more spikelets per panicle.At the same N application rate and compared with inbred rice cultivars,the IJHR cultivar could not only absorb more N,but also efficiently utilize the absorbed N to produce more grain yield,leading to higherIEN,RENandAEN.Improved shoot and root traits of the IJHR cultivar,including LAI,LAD,CGR,NAR,grain–leaf ratio,leaf photosynthetic rate,photosynthetic NUE,the remobilization of NSC from stems,root length,root biomass,and ROA,contributed to its higher grain yield and NUE.The higher contents of Z+ZR in both roots and leaves of the IJHR played a vital role in its response to N application rates and benefited other physiological processes.An N application rate at 300 kg ha–1in the soil with medium-fertility,or specific leaf N content of approximately 3.4 and 2.9 g m–2,respectively,at the panicle initiation and the heading time,could be used as a rational N application index to achieve both high grain yield and high NUE for the IJHR cultivar.

    Acknowledgements

    We are grateful for grants from the National Natural Science Foundation of China (32071943),the Priority Academic Program Development of Jiangsu Higher Education Institutions,China (PAPD-2020-01),the Postgraduate Research and Innovation Program of Jiangsu Province,China (XKYCX17_052),and the Top Talent Supporting Program of Yangzhou University,China(2015-01).

    Declaration of competing interest

    The authors declare that they have no conflict of interest.

    Appendicesassociated with this paper are available on http://www.ChinaAgriSci.com/V2/En/appendix.htm

    啦啦啦观看免费观看视频高清| 成人国产麻豆网| 欧美日韩国产亚洲二区| 日韩欧美在线乱码| 白带黄色成豆腐渣| 永久免费av网站大全| 午夜福利在线观看免费完整高清在| 九草在线视频观看| 中文亚洲av片在线观看爽| 国产乱人视频| 免费观看a级毛片全部| av在线观看视频网站免费| 99久久九九国产精品国产免费| 久久久久久久久久久丰满| 热99re8久久精品国产| 能在线免费看毛片的网站| 卡戴珊不雅视频在线播放| 一级毛片久久久久久久久女| 国内精品一区二区在线观看| 男女边吃奶边做爰视频| 99在线人妻在线中文字幕| 九九久久精品国产亚洲av麻豆| 国产精品无大码| 少妇人妻一区二区三区视频| 啦啦啦观看免费观看视频高清| 色综合亚洲欧美另类图片| 久久久久久久久久黄片| 国产精品一区二区三区四区久久| 免费人成在线观看视频色| 国产精品综合久久久久久久免费| 亚洲精品,欧美精品| 性色avwww在线观看| 日韩人妻高清精品专区| 成人一区二区视频在线观看| 国产91av在线免费观看| 久久99热这里只频精品6学生 | 波多野结衣巨乳人妻| 日韩亚洲欧美综合| 久久精品影院6| 免费观看a级毛片全部| 国产午夜精品一二区理论片| 在线免费十八禁| 国产精品福利在线免费观看| 国产亚洲av片在线观看秒播厂 | 免费黄网站久久成人精品| 啦啦啦韩国在线观看视频| 久久久久久久久大av| 国产精品综合久久久久久久免费| 亚洲精华国产精华液的使用体验| 亚洲欧美一区二区三区国产| 国产亚洲最大av| 欧美三级亚洲精品| 丝袜喷水一区| 三级毛片av免费| 日日摸夜夜添夜夜爱| 国内揄拍国产精品人妻在线| 国产精品国产三级专区第一集| 成年版毛片免费区| 97超碰精品成人国产| 精品无人区乱码1区二区| 岛国在线免费视频观看| 国产男人的电影天堂91| 欧美一区二区精品小视频在线| 久久久久久久久中文| 啦啦啦观看免费观看视频高清| 国产午夜精品一二区理论片| 丝袜喷水一区| 看黄色毛片网站| 国产在视频线精品| 欧美变态另类bdsm刘玥| 搡女人真爽免费视频火全软件| 欧美三级亚洲精品| 看黄色毛片网站| 成人一区二区视频在线观看| 国产av一区在线观看免费| 亚洲精品影视一区二区三区av| 免费不卡的大黄色大毛片视频在线观看 | av在线老鸭窝| av播播在线观看一区| 午夜福利网站1000一区二区三区| 级片在线观看| 三级男女做爰猛烈吃奶摸视频| 国产一区有黄有色的免费视频 | 国产精品福利在线免费观看| 欧美又色又爽又黄视频| 久久精品国产自在天天线| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲5aaaaa淫片| 亚洲欧美成人综合另类久久久 | 国产精品久久久久久久久免| 国产日韩欧美在线精品| 内射极品少妇av片p| 日本黄色视频三级网站网址| 两个人的视频大全免费| 国内少妇人妻偷人精品xxx网站| 日本免费一区二区三区高清不卡| 亚洲欧美成人综合另类久久久 | 日韩在线高清观看一区二区三区| 国产伦一二天堂av在线观看| 永久网站在线| 夜夜看夜夜爽夜夜摸| 91午夜精品亚洲一区二区三区| 国产伦精品一区二区三区视频9| 国产成人精品婷婷| 一级二级三级毛片免费看| 卡戴珊不雅视频在线播放| 久久精品久久久久久久性| 青春草视频在线免费观看| 真实男女啪啪啪动态图| 国产一区有黄有色的免费视频 | 狂野欧美白嫩少妇大欣赏| 免费看a级黄色片| 亚洲最大成人中文| 99国产精品一区二区蜜桃av| 亚洲最大成人av| 狂野欧美激情性xxxx在线观看| 日韩高清综合在线| 久99久视频精品免费| 麻豆乱淫一区二区| 日本黄色片子视频| 久久久欧美国产精品| 永久免费av网站大全| av免费在线看不卡| 岛国在线免费视频观看| 亚洲欧美清纯卡通| 99热这里只有精品一区| 人人妻人人看人人澡| 亚洲经典国产精华液单| 2022亚洲国产成人精品| 成人国产麻豆网| 在线观看美女被高潮喷水网站| 噜噜噜噜噜久久久久久91| 久久久久网色| 国产久久久一区二区三区| 级片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美成人精品一区二区| 日韩国内少妇激情av| 亚洲av中文字字幕乱码综合| 亚洲国产高清在线一区二区三| 亚洲精品日韩av片在线观看| 男女边吃奶边做爰视频| 天天躁日日操中文字幕| 美女高潮的动态| 午夜久久久久精精品| 国产毛片a区久久久久| 1000部很黄的大片| 国产免费福利视频在线观看| av在线天堂中文字幕| 免费黄网站久久成人精品| 精品熟女少妇av免费看| 免费看美女性在线毛片视频| 乱码一卡2卡4卡精品| 永久免费av网站大全| 欧美性猛交╳xxx乱大交人| 亚洲欧美精品综合久久99| 看免费成人av毛片| 国产男人的电影天堂91| 中文字幕亚洲精品专区| 身体一侧抽搐| 最近视频中文字幕2019在线8| 嘟嘟电影网在线观看| 欧美最新免费一区二区三区| 国产亚洲最大av| 欧美另类亚洲清纯唯美| 国产高潮美女av| 舔av片在线| 中文字幕亚洲精品专区| 久久热精品热| 美女内射精品一级片tv| 国产欧美日韩精品一区二区| 婷婷六月久久综合丁香| 国产真实乱freesex| 色5月婷婷丁香| 久久久成人免费电影| 男人的好看免费观看在线视频| 看黄色毛片网站| 午夜a级毛片| 国内精品宾馆在线| 国产精品人妻久久久久久| 看非洲黑人一级黄片| 午夜精品在线福利| 伊人久久精品亚洲午夜| 在线免费观看的www视频| 亚洲欧美一区二区三区国产| 欧美日韩综合久久久久久| a级毛片免费高清观看在线播放| 一级av片app| 欧美区成人在线视频| 国产视频首页在线观看| 最近最新中文字幕大全电影3| 亚洲精品日韩在线中文字幕| 日本av手机在线免费观看| 如何舔出高潮| 免费看光身美女| 黄色配什么色好看| 免费看光身美女| 人人妻人人澡欧美一区二区| 国产亚洲最大av| 国产v大片淫在线免费观看| 精品不卡国产一区二区三区| 一边亲一边摸免费视频| 日本一二三区视频观看| 久久久精品大字幕| av在线蜜桃| 国产精品爽爽va在线观看网站| 欧美丝袜亚洲另类| 欧美日韩国产亚洲二区| 26uuu在线亚洲综合色| 男女视频在线观看网站免费| 国产精品麻豆人妻色哟哟久久 | .国产精品久久| 五月玫瑰六月丁香| 午夜精品在线福利| 久久久久久久久中文| 亚洲真实伦在线观看| 国产精品爽爽va在线观看网站| 亚洲成人久久爱视频| 国产一区有黄有色的免费视频 | 不卡视频在线观看欧美| 久久久久久久久久久免费av| 少妇丰满av| 久久韩国三级中文字幕| 18禁动态无遮挡网站| 欧美激情久久久久久爽电影| 99久久精品热视频| 午夜老司机福利剧场| 国产精品久久久久久精品电影| 欧美性猛交╳xxx乱大交人| 欧美成人精品欧美一级黄| 精品不卡国产一区二区三区| 国产一区二区亚洲精品在线观看| 好男人在线观看高清免费视频| 99热这里只有精品一区| 成人亚洲欧美一区二区av| 亚洲欧洲日产国产| 国产精品美女特级片免费视频播放器| 亚洲欧美成人精品一区二区| 久久久久精品久久久久真实原创| 99久久无色码亚洲精品果冻| 91av网一区二区| 国产精品国产三级国产专区5o | 大话2 男鬼变身卡| 中文精品一卡2卡3卡4更新| 久久午夜福利片| 男人狂女人下面高潮的视频| 一个人看的www免费观看视频| 久久这里只有精品中国| 欧美成人一区二区免费高清观看| 别揉我奶头 嗯啊视频| 网址你懂的国产日韩在线| 亚洲怡红院男人天堂| 99久久精品国产国产毛片| 国产探花在线观看一区二区| 国产三级在线视频| 男女视频在线观看网站免费| 亚洲人成网站高清观看| 亚洲精品aⅴ在线观看| 国产高清国产精品国产三级 | 亚洲国产精品合色在线| 国产精品,欧美在线| 能在线免费观看的黄片| 国产伦理片在线播放av一区| 长腿黑丝高跟| 久久久久久久亚洲中文字幕| 婷婷色综合大香蕉| 国产中年淑女户外野战色| 免费观看在线日韩| 六月丁香七月| 成年版毛片免费区| 男女边吃奶边做爰视频| 联通29元200g的流量卡| 亚洲18禁久久av| 亚洲精品乱码久久久v下载方式| 午夜福利在线在线| 久久精品夜色国产| 国产极品精品免费视频能看的| av在线播放精品| 国产欧美日韩精品一区二区| 午夜a级毛片| 成年免费大片在线观看| 日韩一本色道免费dvd| 最近手机中文字幕大全| 国产伦一二天堂av在线观看| 国产精品精品国产色婷婷| 成人一区二区视频在线观看| 亚洲国产欧洲综合997久久,| 亚洲人成网站高清观看| 久久这里只有精品中国| 99热这里只有精品一区| 26uuu在线亚洲综合色| 韩国高清视频一区二区三区| 蜜臀久久99精品久久宅男| 久久国内精品自在自线图片| 少妇的逼水好多| 亚洲欧洲国产日韩| 亚洲欧美精品自产自拍| 91午夜精品亚洲一区二区三区| 日韩强制内射视频| 亚洲精品456在线播放app| 欧美激情久久久久久爽电影| 我的女老师完整版在线观看| 狂野欧美白嫩少妇大欣赏| 国产大屁股一区二区在线视频| 十八禁国产超污无遮挡网站| 内射极品少妇av片p| 亚洲欧美日韩无卡精品| 亚洲av一区综合| 国产精品永久免费网站| 精品少妇黑人巨大在线播放 | av在线播放精品| 久久人妻av系列| 最近中文字幕2019免费版| 免费电影在线观看免费观看| 一区二区三区乱码不卡18| 国产精品永久免费网站| 国产免费福利视频在线观看| 国国产精品蜜臀av免费| 久久精品国产亚洲网站| 黄色配什么色好看| 黑人高潮一二区| videossex国产| 日韩一区二区视频免费看| 少妇人妻一区二区三区视频| 午夜老司机福利剧场| kizo精华| 国产成人91sexporn| 天天躁日日操中文字幕| 欧美又色又爽又黄视频| 在线观看美女被高潮喷水网站| 丝袜喷水一区| 久久精品夜夜夜夜夜久久蜜豆| 成人性生交大片免费视频hd| 一个人看的www免费观看视频| 亚洲内射少妇av| 18禁动态无遮挡网站| 国产成人精品久久久久久| 国产精品爽爽va在线观看网站| 级片在线观看| 日韩中字成人| 亚洲精品国产av成人精品| 欧美高清成人免费视频www| 国产精品久久久久久精品电影| 青春草国产在线视频| 日本黄大片高清| 激情 狠狠 欧美| 精品无人区乱码1区二区| 国产成人午夜福利电影在线观看| 精品国产三级普通话版| 午夜爱爱视频在线播放| 精品少妇黑人巨大在线播放 | 非洲黑人性xxxx精品又粗又长| 色综合色国产| 国产亚洲91精品色在线| 成人亚洲精品av一区二区| 人妻夜夜爽99麻豆av| 狠狠狠狠99中文字幕| 国产日韩欧美在线精品| 免费观看a级毛片全部| 在现免费观看毛片| 麻豆成人av视频| 人人妻人人看人人澡| 欧美bdsm另类| 亚洲国产色片| 国内精品一区二区在线观看| 日韩制服骚丝袜av| 欧美xxxx黑人xx丫x性爽| 中文字幕精品亚洲无线码一区| 亚洲人成网站在线观看播放| 欧美区成人在线视频| 国产在视频线在精品| 永久网站在线| 舔av片在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av二区三区四区| 亚洲一级一片aⅴ在线观看| 超碰97精品在线观看| 99视频精品全部免费 在线| 免费看美女性在线毛片视频| 日本一本二区三区精品| 国产精品无大码| 国产淫片久久久久久久久| 国产探花极品一区二区| АⅤ资源中文在线天堂| 国产一级毛片在线| 久久热精品热| 午夜免费激情av| 国产成人freesex在线| 国产国拍精品亚洲av在线观看| 亚洲av不卡在线观看| 成人毛片a级毛片在线播放| 丝袜美腿在线中文| 看黄色毛片网站| av又黄又爽大尺度在线免费看 | 成人一区二区视频在线观看| 成人国产麻豆网| 少妇人妻精品综合一区二区| 色网站视频免费| 国产黄色小视频在线观看| 久久久a久久爽久久v久久| 午夜福利网站1000一区二区三区| 久久久国产成人精品二区| 少妇的逼水好多| 国产精品一及| 超碰97精品在线观看| 成人亚洲欧美一区二区av| 成人特级av手机在线观看| av线在线观看网站| 亚洲精品乱码久久久久久按摩| 18禁在线无遮挡免费观看视频| a级一级毛片免费在线观看| 欧美日韩精品成人综合77777| 久久久国产成人免费| 国产伦在线观看视频一区| 国国产精品蜜臀av免费| 99久久人妻综合| 国产女主播在线喷水免费视频网站 | 麻豆久久精品国产亚洲av| 亚洲国产精品成人综合色| 91精品一卡2卡3卡4卡| 国产精品一区二区三区四区久久| 精品99又大又爽又粗少妇毛片| 日韩av在线免费看完整版不卡| 18禁在线无遮挡免费观看视频| 精品少妇黑人巨大在线播放 | 国产亚洲精品av在线| 国产男人的电影天堂91| 中文字幕人妻熟人妻熟丝袜美| 国产精品国产三级专区第一集| 干丝袜人妻中文字幕| 少妇丰满av| 日韩欧美国产在线观看| 国产av不卡久久| 欧美变态另类bdsm刘玥| 青青草视频在线视频观看| 伦理电影大哥的女人| 午夜福利在线观看吧| 国产免费福利视频在线观看| 99热全是精品| 少妇被粗大猛烈的视频| 国产精品麻豆人妻色哟哟久久 | 久99久视频精品免费| 国产精品久久久久久久电影| 国产精品久久久久久精品电影小说 | 久久韩国三级中文字幕| 成年女人看的毛片在线观看| 韩国av在线不卡| 精品久久久久久久久亚洲| av线在线观看网站| 九色成人免费人妻av| 久久久久免费精品人妻一区二区| 91精品伊人久久大香线蕉| 成年女人看的毛片在线观看| 久久久久性生活片| 黄片wwwwww| av在线老鸭窝| 嫩草影院精品99| 国产精品一区二区三区四区免费观看| 日韩高清综合在线| 日韩精品青青久久久久久| 亚洲欧美精品综合久久99| 七月丁香在线播放| 人人妻人人看人人澡| 亚洲欧美日韩卡通动漫| 熟女电影av网| 精品无人区乱码1区二区| 日韩中字成人| 久久久精品大字幕| 超碰av人人做人人爽久久| 白带黄色成豆腐渣| 午夜精品在线福利| 黄片无遮挡物在线观看| 亚洲天堂国产精品一区在线| 男女国产视频网站| www.av在线官网国产| 禁无遮挡网站| 婷婷色麻豆天堂久久 | 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美成人精品一区二区| 一个人观看的视频www高清免费观看| 在线天堂最新版资源| 国产精品一区www在线观看| 91av网一区二区| 国产综合懂色| 听说在线观看完整版免费高清| 亚洲av二区三区四区| 亚洲国产色片| www日本黄色视频网| 乱码一卡2卡4卡精品| 欧美性猛交黑人性爽| 大香蕉久久网| 在线免费十八禁| 美女大奶头视频| av在线老鸭窝| 深爱激情五月婷婷| av天堂中文字幕网| 中文字幕精品亚洲无线码一区| 在线观看一区二区三区| 国产精品电影一区二区三区| av视频在线观看入口| 亚洲中文字幕一区二区三区有码在线看| a级毛色黄片| 亚洲欧美成人综合另类久久久 | 日日干狠狠操夜夜爽| 天堂中文最新版在线下载 | 亚洲av一区综合| 欧美最新免费一区二区三区| 51国产日韩欧美| 免费播放大片免费观看视频在线观看 | 久久久精品94久久精品| 99久久九九国产精品国产免费| 久久久久网色| 黄片wwwwww| 免费观看精品视频网站| 国产高清有码在线观看视频| 亚洲精品aⅴ在线观看| 亚洲欧美精品综合久久99| 中文精品一卡2卡3卡4更新| 亚洲在久久综合| 亚洲最大成人手机在线| 精品人妻视频免费看| 国产大屁股一区二区在线视频| 最近视频中文字幕2019在线8| 夫妻性生交免费视频一级片| 一级黄片播放器| www.色视频.com| 国产精品人妻久久久久久| 欧美xxxx黑人xx丫x性爽| 波多野结衣高清无吗| 亚洲精品乱码久久久v下载方式| 亚洲国产欧美人成| 黄色日韩在线| 免费在线观看成人毛片| av免费在线看不卡| 亚洲,欧美,日韩| 亚洲四区av| 久久国产乱子免费精品| 亚洲av成人精品一区久久| 国产黄色小视频在线观看| 大话2 男鬼变身卡| 中文字幕亚洲精品专区| 国产精品三级大全| 成年版毛片免费区| 午夜激情福利司机影院| 在线观看一区二区三区| 波多野结衣巨乳人妻| 欧美xxxx黑人xx丫x性爽| 国产女主播在线喷水免费视频网站 | 观看免费一级毛片| 日本免费a在线| 欧美日韩一区二区视频在线观看视频在线 | 搡女人真爽免费视频火全软件| 麻豆一二三区av精品| 国产成人a∨麻豆精品| 成人毛片a级毛片在线播放| 九九爱精品视频在线观看| 国产视频首页在线观看| 中文亚洲av片在线观看爽| 中文字幕精品亚洲无线码一区| 久久久久久伊人网av| 亚洲成av人片在线播放无| 午夜老司机福利剧场| 亚洲自拍偷在线| 国产探花在线观看一区二区| 一个人观看的视频www高清免费观看| 男女啪啪激烈高潮av片| 亚洲天堂国产精品一区在线| 视频中文字幕在线观看| 国产高潮美女av| 99久久人妻综合| 18禁在线播放成人免费| 人妻夜夜爽99麻豆av| 亚洲美女搞黄在线观看| 乱人视频在线观看| 日本黄色视频三级网站网址| 国产精品蜜桃在线观看| 欧美性猛交黑人性爽| 国产精品久久久久久久电影| 晚上一个人看的免费电影| 又粗又硬又长又爽又黄的视频| 成年版毛片免费区| 99热这里只有是精品在线观看| 最近最新中文字幕免费大全7| 欧美日韩国产亚洲二区| 精品一区二区三区视频在线| 日本wwww免费看| 伦精品一区二区三区| 在线免费十八禁| 国产成人91sexporn| 熟妇人妻久久中文字幕3abv| av女优亚洲男人天堂| 精品一区二区免费观看| 国产免费男女视频| 亚洲成人久久爱视频| 国产私拍福利视频在线观看| 久久鲁丝午夜福利片| 国产淫片久久久久久久久| 亚洲伊人久久精品综合 | 91久久精品国产一区二区三区| 国产老妇女一区| 如何舔出高潮| 国产精品爽爽va在线观看网站| 国产午夜福利久久久久久| 18禁动态无遮挡网站| 久久久精品欧美日韩精品| 国产精品一区二区性色av| 国产精品久久电影中文字幕| 一级二级三级毛片免费看| 亚洲精品aⅴ在线观看| av黄色大香蕉| 一区二区三区四区激情视频| 婷婷色麻豆天堂久久 | 一夜夜www| 亚洲av中文av极速乱|