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

    Responses of Lowland, Upland and Aerobic Rice Genotypes to Water Limitation During Different Phases

    2020-07-06 10:28:48PreethiVijayaraghavareddyYinXinyouPaulStruikUdayakumarMakarlaSheshshayeeSreeman
    Rice Science 2020年4期

    Preethi Vijayaraghavareddy, Yin Xinyou, Paul C. Struik, Udayakumar Makarla Sheshshayee Sreeman

    Research Paper

    Responses of Lowland, Upland and Aerobic Rice Genotypes to Water Limitation During Different Phases

    Preethi Vijayaraghavareddy1,2, Yin Xinyou2, Paul C. Struik2, Udayakumar Makarla1, Sheshshayee Sreeman1

    (; )

    Rice yield reduction due to water limitation depends on its severity and duration and on the phenological stage of its occurrence. We exposed three contrasting rice genotypes, IR64, UPLRi7 and Apo (adapted to lowland, upland and aerobic conditions, respectively), to three water regimes (puddle, 100% and 60% field capacity) in pots during the vegetative (GSI), flowering (GSII) and grain filling (GSIII) stages. Stress at all the three stages significantly reduced yield especially in lowland genotype IR64. Effect of water limitation was more severe at GSII than at the other two stages. Stress at GSI stage reduced both source activity (leaf area and photosynthetic rate) and sink capacity (tiller number or panicle number per pot). When stress was imposed at GSII, spikelet fertility was most affected in all the three genotypes. In both GSII and GSIII, although leaf area was constant in all the three water regimes, estimated relative whole-plant photosynthesis was strongly associated with yield reduction. Reduced photosynthesis due to stress at any given stage was found to have direct impact on yield. Compared to the other genotypes, Apo had deeper roots and maintained a better water relation, thus, higher carbon gain and spikelet viability, and ultimately, higher biomass and productivity under water-limited conditions. Therefore, screening for these stage-dependent adaptive mechanisms is crucial in breeding for sustained rice production under water limitation.

    water limitation; phenology; upland; aerobic;; yield

    Rice is a major staple food crop, providing 35%–80% of total calories required for populations of Asia (IRRI, 1997; Yang et al, 2009). It is mainly grown in semi-aquatic environments in paddies, and rice production in such a system requires 3000–4000 L of water to produce 1 kg of grains (Bouman, 2009). Most rice production is in irrigated lowland ecosystems in Asia. However, the industrial and household requirements are rapidly withdrawing fresh water away from agriculture in this region. This urgently necessitates efforts towards devising water saving agronomies as well as improving (or at least sustaining) rice productivity under water-limited conditions (Bouman and Tuong, 2001).

    Many efforts have been made to modify cultivation practices, like SRI (system of rice intensification) (Uphoff et al, 2010), upland (Nishimura et al, 2008), rainfed lowland (Zeigler and Puckridge, 1995) and aerobic (Kadiyala et al, 2012) systems, to reduce water consumption in rice production. The emphasis in these systems has generally been to decrease the ‘E’ component of ‘ET’ (Evapo-Transpiration). However, these systems have one thing in common: yields are reduced compared to the irrigated lowland (paddy) ecosystem. Efforts have been made to identify genotypes with good yields in upland and aerobic conditions. Although these genotypes were selected for cultivation in water-limited conditions, yield penalties are still severe. Under water-limited conditions, reduced field capacity (FC) has a major impact on yield depending on crop growth stage that coincides with the stress, its intensity and duration. Understanding the physiological mechanisms that get affected under stress and assessing the variability in these mechanisms among various genotypes is crucial for further crop improvement.

    Three growth stages, i.e. vegetative (GSI), flowering (GSII) and grain filling (GSIII) stages, are characterized as distinctive growth periods that are differentially influenced by biotic or abiotic stresses (Boonjung and Fukai, 1996; Sarvestani et al, 2008). At all phenological stages, water limitation potentially affects size and capacity of both source and sink for rice growth. Stress at GSI can impede leaf expansion (thus reducing source size) and inhibit tillering (thus reducing both sink and source). The carbon assimilatory capacity of the canopy remains small due to a reduced photosynthetic rate combined with the reduction in canopy leaf area. These factors can lead to lower final biomass production (Praba et al, 2009). GSII is the most critical stage of yield formation, and water limitation at this stage can adversely affect both carbon source and reproductive development, resulting in reduced yields (Jin et al, 2013). GSIII also represents a critical stage, during which reduction in the grain yield can be caused by reduced leaf photosynthesis as well as lower green leaf area by advanced leaf senescence.

    Rice varieties that adapt to different ecosystems with diverse adaptive strategies to grow under water-limited conditions have been developed, with better extraction of water by roots and increased cellular-level tolerance (CLT) under stress (Raju et al, 2014), increased water use efficiency (WUE) (Sheshshayee et al, 2003) andreduced leaf transpiration through wax layers (Zhang et al, 2007). Plants with good root systems can maintain better tissue water potential with cooler canopy (Ramu et al, 2016) and that plants with high CLT have improved spikelet fertility in stress conditions (Raju et al, 2014). Significant progress has also been made in developing suitable cultivars for greater adaptation to diverse rice growing agronomic practices, which leads to the identification of UPLRi7 for upland conditions and Apo for aerobic cultivation (Ouyang et al,2017). IR64, a mega variety in Southeast Asian (Mackill et al, 2018), is predominantly grown in puddled conditions and is known to be sensitive to water limitation (Dharmappa et al, 2019).

    The objective of this study was to understand the stress response of these three rice genotypes known to adapt to upland, aerobic and lowland rice ecosystems, respectively, and to assess their yielding ability under stress imposed at GSI, GSII and GSIII stages.

    Materials and methods

    Plant growth and stress conditions

    Three rice genotypes, which were bred for lowland (IR64), upland (UPLRi7) and aerobic (Apo) conditions, were selected because they had similar phenology. For instance, all the three genotypes took 85–90 d to reach 50% flowering (Ouyang et al,2017). Pot experiments were conducted in the University of Agricultural Sciences, Bengaluru, India. Experiments were done in an open area under a mobile rain-out shelter, which protected plants against rainfall. Whenever required and during nights the shelter was drawn over the experimental area, thus pots were maintained at the specified water regime. Pots (24 L) were filled with equal amounts (20kg) of red soil with recommended fertilizers added. Direct sowing was followed, and thinning was done at 21 days after sowing (DAS) to maintain one plant per pot. Plants were well watered manually until the stress imposition. During treatments, three water regimes were maintained viz., puddle, 100% FC and 60% FC. Puddle treatment was created with standing water of 5 cm above the soil in the pot. 100% FC and 60% FC treatments were maintained following the gravimetric approach (Udayakumar et al, 1998): the pots were weighed daily using a hanging load cell balance (Essae-Teraoka, Japan), the reduction in weight of the pots as a result of evapo-transpiration was recorded, and on the early morning the pots were replenished with the exact same amount of water evapo-transpired in the preceding day to maintain a given FC. There were five replicates for each treatment and arranged in randomized-block design. Three separate experiments were conducted for each of the three phenological stages.

    GSI experiment was conducted from January to April of 2016, in which the three water regimes were maintained at the time of maximum tillering stage (30 to 50 DAS).

    GSII experiment was conducted from July to October 2016, in which the three water regimes were maintained particularly at the flowering (85 DAS). Water was withheld from flag-leaf stage, such that the required stress level was reached at the time of flowering, and maintained for 10 d.

    GSIII experiment was conducted from July to October of 2017, in which the three water regimes were maintained after flowering (100 DAS) until grain maturity.

    Relative water content

    To assess the effect of different water regimes on tissue water status, relative water content (RWC) was measured. Ten leaf discs of 1cm2each were punched out of the second fully expanded leaf counted from the apex. Fresh weight was recorded immediately. The leaf discs were immersed in distilled water in a beaker for 6 h,and were gently blotted on a filter paper to remove the water adhering to the surface. Fresh weight of this tissue was taken to represent turgid weight. Then, the samples were dried in an oven (65 oC–70oC) until constant weight, and then dry weight was recorded. RWC (%) was estimated according to Babitha et al (2015).

    Gas exchange

    Leaf net-photosynthetic rate () and stomatal conductance (s) were measured with a portable gas-exchange system Li-6400 (LiCOR-Inc, Lincoln, Nebraska, USA). The ambient concentration of CO2in the leaf chamber of the Li-6400 was kept at 701.8 mg/m3, and the photosynthetically active radiation (PPFD) was 1500μmol/(m2·s). All measurements were made at 5 d after stress imposition on the second fully expanded leaf from the top for the GSI experiment, and on the flag leaf for the GSII and GSIII experiments.

    Soil and plant analyzer develotrnent (SPAD) value

    SPAD value was calculated by averaging nine readings per leaf using a portable, non-destructive chlorophyll meter (SPAD-502, Minolta, Japan). Mean of the measurements was calculated to arrive at a single value per pot. Measurements were recorded at 5 d after reaching 60% FC.

    Morpho-physiological parameters

    Plant height and number of tillers were recorded before harvest. Plant height was measured from surface of soil to tip of the longest leaf. Roots were separated from shoot to measure root length before drying. All the plant parts were oven-dried at 70oC to record stem weight, leaf weight and root weight, and total biomass per pot (TBM) were calculated. Specific leaf area was calculated by dividing fresh leaf area by its dry biomass. Total leaf area (TLA) per pot was calculated by multiplying final dry leaf biomass with specific leaf area, which was determined from a subsample of leaves.

    Yield parameters

    Panicles were harvested from individual plants after maturity. Grains were manually separated from panicles. Spikelet fertility was calculated as the number of filled grains divided by the total number of grains.Weight of filled grains was recorded to obtain grain yield per pot.

    Statistical analysis

    Two-way ANOVA was conducted by using GenStat (15thedn) (http://www.genstat.co.uk/). The generated least significant difference for each parameter was used to check the significance level in all the three experiments. Standard correlation/regression procedures were performed using Microsoft Excel.

    Results

    Effect of stress at GSI (vegetative stage)

    RWC, a simple proxy for water relations, decreased significantly as stress level increased in all the genotypes (Fig. 1-A). IR64 showed 10.1% and 17.1% reductions in RWC under 100% FC and 60% FC, respectively, relative to the puddle treatment, while UPLRi7(10.0% and 13.9%) and Apo (0.5% and 2.8%) showed significantly lower reductions in RWC (Fig. 1-A). SPAD value was not affected in any of the genotypes (Supplemental Fig. 1-A). Gas exchange parameters like stomatal conductance (s) and net assimilation rate () decreased with increasing stress in all the genotypes, albeit to different extents. Whiledecreased by 33.0% and 54.8% for IR64 under 100% FC and 60% FC compared to puddle, this reduction was much less for UPLRi7 (11.7% and 34.0%) and Apo (19.4% and 27.5%) (Fig.1-B). UPLRi7 maintained higher stomatal conductance compared to the other two genotypes and it was highest with 0.7 mol/(m2·s) under puddle condition.sdecreased significantly in all three genotypes with highest reduction for IR64 and UPLRi7 under 60% FC (Fig. 1-C). Apo recorded lowersunder the two stress regimes and a lower decrease insunder 60% FC compared to the other two genotypes (Fig. 1-C).

    Fig. 1. Effects of water limitation on various physiological characteristics of three contrasting genotypes.

    A,Relative water content (RWC) at the vegetative stage (GSI).B, Photosynthetic rate () at GSI.C,Stomatal conductance (g) at GSI. D, RWC at the flowering stage (GSII).E,at GSII.F,gat GSII. G, RWC at the grain filling stage (GSIII).H,at GSIII.I,Stomatal conductance (g) at GSIII.

    FC, Field capacity. Data indicate Mean ± SE (= 5).* and ** indicate significant differences from puddle at the 0.05 and 0.01 levels within genotypes, respectively.

    The reduction in soil moisture status resulted in a significant reduction in TLA in all the three genotypes (Table 1). UPLRi7 and IR64 showed a stronger reduction in TLA under 100% FC than Apo. Under 60% FC, all the genotypes showed a very strong reduction with IR64 showing the strongest one (57.7%) (Table 1). Both plant height and tiller number were sensitive to water limitation and IR64 showed a higher reduction than the other two genotypes. Although reduction in tiller number was low in Apo, plant height was significantly reduced compared to puddle conditions (Supplemental Table 1).Whereas, both IR64 and UPLRi7 showed significant reduction in tiller number and plant height under 100% FC and 60% FC. Stem weight decreased more in IR64 and Apo than in UPLRi7 with increased water limitation (Supplemental Table 1).

    Reduction in plant height and tiller number closely matched the reduction in total biomass (TBM). Even under 100% FC, an overall reduction in TBM ranging from 16.8% to 25.0% was noticed (Table 1). Under 60% FC, IR64 recorded the strongest reduction in TBM. This effect of water limitation on physiological responses led to an average yield decrease of 13.2% and 27.6% under 100% FC and 60% FC, respectively (Table 2). IR64 recorded the highest reduction by 35.1% under 60% FC. Similarly, UPLRi7 was also sensitive to vegetative stress with a reduction of 17.7% under 100% FC and 31.7% under 60% FC. Apo showed the least reduction in yield under both 100% FC and 60% FC (6.0% and 16.2%, respectively).

    Stress response at GSII (flowering stage)

    There was a significant reduction in RWC under both 100% FC and 60% FC over puddle condition in all the genotypes. This reduction was strong in IR64, especially under 60% FC (Fig. 1-D). Reduction in chlorophyll content was highly significant under 60% FC for IR64 and UPLRi7 (Supplemental Fig. 1). Water limitation at GSII affected both assimilation rate and stomatal conductance (Fig. 2-E and -F). IR64 showed more than 45% reduction inandswhile both UPLRi7 and Apo being drought adaptive also showed >24% of reductions inandsunder 60% FC.

    Table 1.Effects of water limitation on total leaf area (TLA, cm2/plant) and total biomass (TBM, g/plant).

    GSI, Vegetative stage; GSII, Flowering stage; GSIII, Grain filling stage.

    * and ** indicate significant differences from puddle at the 0.05 and 0.01 levels within genotypes, respectively.

    As expected, water limitation during the GSII experiment had no considerable influence on plant height and tiller number (Supplemental Table 1). Similarly, TLA showed no significant changes in response to stress (Table 1). Reduction of TBM was marginal in all the genotypes under 60% FC. A mild water limitation (100% FC) under GSII significantly reduced yield. This reduction in yield was stronger when severity of stress increased to 60% FC. An overall yield loss of 22.3% and 44.4% was recorded under 100% FC and 60% FC, respectively. IR64 had 50.7% yield loss and UPLRi7 had 46.1%. Apo showed a significant reduction of 36.3% under 60% FC over the puddle condition (Table 2).

    Stress effect at GSIII

    RWC was significantly lower for IR64 than for UPLRi7 and Apo when imposed to stress. The latter two genotypes showed a marginal reduction (upto 6.5%) in RWC under 60% FC compared with the puddle treatment (Fig. 1-G). Under 60% FC, RWC of IR64 decreased by more than 19%. Gas exchange parameters recorded a similar trend. IR64 showed marked reduction in bothandsin both the water regimes compared to puddle treatment (Fig. 1-H and -I).

    SPAD values were generally lower in IR64 than in UPLRi7 and Apo at all the phenological stages. Chlorophyll content of IR64 decreased significantly under both 100% FC and 60% FC (Supplemental Fig. 1). On the other hand, there was not much change in leaf area, plant height and tiller number in all the three genotypes (Table 1 and Supplemental Table 1). Shoot weight and TBM were significantly reduced in IR64 and UPLRi7 especially under 60% FC. An overall yield loss of 12.2% to 27.4% was recorded under 100% FC and 60% FC, respectively (Table 2). Under 60% FC, IR64 and UPLRi7 had the highest reduction (41.6% and 32.4%, respectively), whereas Apo had a marginal reduction in yield (9.8%).

    Responses of spikelet fertility to stress at three stages

    When stress was imposed during the GSI stage, a slight reduction in spikelet fertility was noticed (Fig. 2-A). A decrease in spikelet fertility was the most dominant yield component accounting for the observed yield reduction at GSII in all the three genotypes. IR64 had the greatest reduction in spikelet fertility followed by UPLRi7 (Fig. 2-B). Although Apo showed smaller reduction in spikelet fertility than the other genotypes, a reduction of 36.3% under 60% FC was observed irrespective of its adaptability to water-limited conditions. Stress at GSIII reduced the number of filled grains especially under 60% FC. IR64 had the highest reduction (18.0%) and Apo had the lowest reduction (9.5%) under 60% FC (Fig. 2-C).

    Table 2.Yield for the three contrasting genotypes at different growth stages and stress treatments. g/plant

    GSI, Vegetative stage; GSII, Flowering stage; GSIII, Grain filling stage; FC, Field capacity.

    * and ** indicate significant differences from puddle at the 0.05 and 0.01 levels within genotypes, respectively. Values in the parentheses are the reduced percent compared to the puddle condition.

    Aerobic cultivar Apo had higher water extraction ability under stress

    Significant effects of stress treatments at different growth stages on root traits were observed for all the three genotypes. In all the three growth stages, Apo and UPLRi7 maintained higher root growth compared to IR64 under stress. At GSI, root length was increased significantly in all the genotypes under 60% FC. Apo showed the highest increase of 40.8% and 42% under 100% FC and 60% FC, respectively (Fig. 2-D). UPLRi7 also showed a significant increase of 18.2% and 31.3% under 100% FC and 60% FC, respectively, compared to the puddle condition. Interestingly, cultivar with the highest root length showed the lowest root weight (Supplemental Fig. 2). UPLRi7 and Apo had the highest root length under 60% FC but root weight was significantly lower for these two genotypes than for IR64.

    The increasing trend of root length in all the genotypes under stress conditions during GSI was not observed in GSII and GSIII. Root weight was marginally decreased by stress in GSII (Supplemental Fig. 2). At GSIII, IR64 and UPLRi7 showed marginal reduction in root length under 60% FC (Fig. 2-F).

    Source to sink imbalance was the major contributor for yield loss at any given stage

    At all the three growth stages, significant reduction in photosynthesis under 100% FC and 60% FC over puddle was noticed. Reductions inand yield were significantly related in all the three stages (Fig. 3-A). However, a strong linear relationship was shown between reduction in total leaf area and yield only at GSI stage (Fig. 3-B).

    We computed an index as an indication of relative whole-plant photosynthesis by multiplying total leaf area with single leaf photosynthetic rate (the term ‘relative’ is used because the estimations are not absolute values which would require considering the vertical decrease ofwith lowering leaf positions). Reduction in this index was strongly related with reduction in yield in all the three phenological stages (Fig. 3-C).

    Reproductive growth (sink) was examined along with photosynthetic rate (source activity) to explain yield loss at all the three stages. Reduction inand spikelet fertility were significantly correlated at GSII and GSIII (Fig. 3-D).

    Fig.2. Effects of stress on spikelet fertility and root length among the genotypes.

    A,Spikelet fertility at vegetative stage.B, Spikelet fertility at flowering stage.C,Spikelet fertilityat grain filling stage. D,Root length at vegetative stage.E, Root length at at flowering stage.F,Root length at at grain filling stage.

    FC, Field capacity. Data indicate Mean ± SE (= 5). * and ** indicate significant differences from puddle at the 0.05 and 0.01 levels within genotypes, respectively.

    Fig. 3. Linear regression among reductions in several traits.

    A, Linear regression of reductions in assimilation rate () and yield. B, Linear regression of reductions in total leaf area (TLA) and yield. C, Linear regression of reductions in computedand yield. D, Linear regression of reductions in assimilation rate () and spikelet fertility (SF).

    GSI, Vegetative stage; GSII, Flowering stage; GSIII, Grain filling stage.

    Reduction in yield over puddle was calculated for the means under 100% and 60% field capacity. ** indicates significance at the 0.01 level.

    Discussion

    Crop improvement strategies are aiming at identification of specific traits with precise phenotyping methods (Vijayaraghavareddy et al, 2017; Sheshshayee et al, 2018). But the major limitation is that response of plants is widely different at specific phenological stage. Further, these physiological responses could potentially vary between genotypes with known history of adaptation to specific agro-ecosystem.

    Stress at any phenological stage reduces plant productivity

    During all the phenological stages, reduced water availability significantly affected plant performance and productivity. Stress at any phenological stage reduced yields by 15.9% (100% FC) to 33.1% (60% FC), compared to puddle conditions (Table 2). For rice being semi-aquatic, even a mild water limitation reduces yield significantly (Yang et al, 2019). IR64 had the greater reduction in yield under 100% FC in all the three phenological stage because of its adaptability to lowland condition,whereas Apo, a genotype adapted to aerobic condition, showed lower reduction under 100% FC over puddle condition in all the three phenological stages. Though this trend was consistent under 60% FC, reduction in yield at GSII was significantly higher even for aerobic genotype Apo. Therefore, among the growth stages, GSII was most sensitive even for mild water limitation (Yang et al, 2019) with an overall yield loss of 22.3% and 44.4% under 100% FC and 60% FC, respectively.

    Increasing root length assists to maintain tissue turgor

    Maintenance of water relations under low soil moisture condition is an important drought adaptive strategy. In aerobic and water limited conditions, root traits are important for water uptake from deeper layers of soils (Raju et al, 2014). Increase in root length was observed at GSI (Fig. 2-D). But the influence of water limitation has no effect on root length at GSII and GSIII stages (Fig. 2-E and -F). The maximum root growth was observed at GSII and GSIII. It is most unlikely that a period of 15 d stress may have a significant influence on physical root length. However, GSI,which represents a period of active growth, showed a maximum increase in root length under water-limited conditions (Fig. 2-D). Root growth is known to be stimulated by soil moisture deficit and significant genetic variability has also been demonstrated for this response (Kadam et al, 2015), and this will be more pronounced if carbon assimilatory capacity and available carbon resources for remobilisation are not significant. Genotypes with better root system can maintain higher RWC due to efficient water mining ability (Nguyen et al, 1997; Raju et al, 2014). This could be an adaptive mechanism in upland genotype UPLRi7 and aerobic genotype Apo, hence showing significant increase in root length at GSI stage (Fig. 2-D).

    Reduced spikelet fertility is the determining factor for reduced yield in all phenological stages

    In GSI, the most prominent influence of stress was reduced leaf expansion and gas exchange. Decreases in canopy cover when plants experience stress during vegetative stage are well documented and have been attributed to reduced cell expansion and cell division (Boonjung and Fukai, 1996; Sarvestani et al, 2008; Praba et al, 2009). Reduction in turgor has a direct effect on carbon assimilation rates (Turk et al, 1980; Praba et al, 2009). The combined effect of reduced leaf area, tiller number and photosynthetic rate significantly impedes total carbon assimilatory capacity (source), which may lead to decreased spikelet fertility. Spikelet fertility was significantly lower at GSII stage than at the other stages (Fig. 2-B). Spikelet fertility at GSII stage is governed by a number of processes ranging from anther development, anther dehiscence, pollen viability, stigmatic receptivity, pollen efficiency, fertilization and ovule abortion, and all of these processes are extremely sensitive to stress (Shi et al, 2018). Besides, the spikelet water relations could further exacerbate the stress effect (O’Toole et al, 1984; Selote and Khanna-Chopra, 2004). Although a significant reduction in spikelet fertility was also noticed at GSIII (Fig. 2-C), factors like lack of availability of photosynthates for grain filling and negative source-sink interactions were the major contributors to reduce spikelet fertility at GSIII. Hence, yield losses arise from reduced source size resulting from a decrease in amount of photosynthates for grain filling, reduced number of spikelets and/or from reduction in efficiency of grain filling (Roitsch, 1999; Madani et al,2010).

    Reduction in source availability will influence yield

    Influence of source on growth and productivity can be viewed from two perspectives: reduction in canopy cover (Cakir, 2004) and reduction in photosynthetic rate per unit leaf area (Serraj et al, 1999). In other words, source influences growth and productivity because of a combination of these two factors collectively responsible for whole-plant photosynthesis. Yield reduction when stress coincides with grain filling (GSIII) could be resulted from lack of carbohydrates for grain filling and/or due to reduced synthesis of current photosynthates and their translocation (Madani et al, 2010). The reduction in leaf area during GSII and GSIII showed no any association with reduction in yield (Fig. 3-B). However, reduction in photosynthetic rate was strongly associated with reduction in yield (Fig. 3-A). Further, the reduction in the relative whole-plant photosynthesis showed a much stronger reduction in yield in all the stages. A change in source-sink interaction has been suggested as one of the major factors responsible for yield reduction under stress (Roitsch, 1999). Stress during the later phenological stages is well known to induce leaf senescence which decreases green leaf area and thus reduces canopy photosynthesis (Jagadish et al, 2015). However, we noticedno considerable reduction in green leaf area due to stress during GSII and GSIII stages (Table 1), which is against the normally noticed trend. There could be two possible reasons. Firstly, stress was provided for a period of 15 d, which may not be severe enough to cause changes in leaf area. Secondly, total leaf area was computed by multiplying leaf dry weight with specific leaf area. As we had only taken the green and functional leaves for recording leaf weight, this method can cause little error in determining functional leaf area.

    Among the growth stages, moisture stress at GSII significantly affects the yield irrespective of adaptability of a genotype. Reductions in biomass, spikelet fertility and carbon gain seem to be major determining factors of yield under stress at GSI, GSII and GSIII stages, respectively. Reduction in yield was significantly high for lowland cultivar IR64 even at a 100% FC, which is a mild water limitation compared with puddle condition. Although the aerobic genotype Apo maintained much higher yield in water limited conditions (100% FC and 60% FC), there is a need of further improvement to reduce yield loss especially when stress induced at GSII stage. Hence, a genotype which maintained source to sink balance like Apo can be used for further crop improvement programmes. This information would therefore be most useful in improving the performance of rice crops which differ in their adaptability.

    Acknowledgements

    This work was supported by an anonymous private donor, via Wageningen University Fund, to the first author’s PhD fellowship. Dr C. G. van der Linden and Dr P. S. Bindraban are acknowledged for their valuable support.

    Supplemental data

    The following materials are available in the online version of this article at http://www.sciencedirect.com/science/ journal/16726308; http://www.ricescience.org.

    Supplemental Table 1. Effect of water limitation on total plant height, number of tillers and stem weight.

    Supplemental Fig. 1. Effect water limitation on SPAD value at GSI, GSII and GSIII in rice genotypes.

    Supplemental Fig. 2. Effect water limitation on root weight at GSI, GSII and GSIII in rice genotypes.

    Babitha KC, Vemanna R S, Nataraja K N, Udayakumar M. 2015. Overexpression of EcbHLH57 transcription factor fromL. in tobacco confers tolerance to salt, oxidative and drought stress.,10(9): e0137098.

    Boonjung H, Fukai S. 1996. Effects of soil water deficit at different growth stages on rice growth and yield under upland conditions: 2. Phenology, biomass production and yield.,48(1): 47–55.

    Bouman B. 2009. How much water does rice use?, 69: 115–133.

    Bouman B A M, Tuong T P. 2001. Field water management to save water and increase its productivity in irrigated lowland rice.,49(1): 11–30.

    ?akir R. 2004. Effect of water stress at different development stages on vegetative and reproductive growth of corn., 89(1): 1–16.

    Dharmappa P M, Doddaraju P, Malagondanahalli M V, Rangappa R B, Mallikarjuna N M, Rajendrareddy S H, Ramanjinappa R, Mavinahalli R P, Prasad T G, Udayakumar M, Sheshshayee S M. 2019. Introgression of root and water use efficiency traits enhances water productivity: An evidence for physiological breeding in rice (L.).,12(1): 14.

    IRRI. 1997. Rice Almanac. 2nd edn. Los Banos, the Philippines: International Rice Research Institute.

    Jagadish K S V, Kavi Kishor P B, Bahuguna R N, von Wiren N, Sreenivasulu N. 2015. Staying alive or going to die during terminal senescence: An enigma surrounding yield stability.,6: 1070.

    Jin Y, Yang H X, Wei Z, Ma H, Ge X C. 2013. Rice male development under drought stress: Phenotypic changes and stage-dependent transcriptomic reprogramming.,6(5): 1630–1645.

    Kadam N N, Yin X, Bindraban P S, Struik P C, Jagadish K S. 2015. Does morphological and anatomical plasticity during the vegetative stage make wheat more tolerant of water deficit stress than rice?,167(4): 1389–1401.

    Kadiyala M D M, Mylavarapu R S, Li Y C, Reddy G B, Reddy M D. 2012. Impact of aerobic rice cultivation on growth, yield, and water productivity of rice-maize rotation in semiarid tropics.,104(6): 1757–1765.

    Madani A, Shirani-Rad A, Pazoki A, Nourmohammadi G, Zarghami R, Mokhtassi-Bidgoli A. 2010. The impact of source or sink limitations on yield formation of winter wheat (L.) due to post-anthesis water and nitrogen deficiencies.,56(5): 218–227.

    Mackill D J, Khush G S. 2018. IR64: A high-quality and high-yielding mega variety.,11(1): 18.

    Nguyen H T, Babu R C, Blum A. 1997. Breeding for drought resistance in rice: Physiology and molecular genetics consideration., 37(5): 1426–1434.

    Nishimura S, Yonemura S, Sawamoto T, Shirato Y, Akiyama H, Sudo S, Yagi K. 2008. Effect of land use change from paddy rice cultivation to upland crop cultivation on soil carbon budget of a cropland in Japan.,125: 9–20.

    O’Toole J C, Hsiao T C, Namuco O S. 1984. Panicle water relations during water stress.,33(2): 137–143.

    Ouyang W J, Struik P C, Yin X Y, Yang J C. 2017. Stomatal conductance, mesophyll conductance, and transpiration efficiency in relation to leaf anatomy in rice and wheat genotypes under drought.,68(18): 5191–5205.

    Praba M L, Cairns J E, Babu R C, Lafitte H R. 2009. Identification of physiological traits underlying cultivar differences in drought tolerance in rice and wheat.,195(1): 30–46.

    Raju B R, Narayanaswamy B R, Mohankumar M V, Sumanth K K, Rajanna M P, Mohanraju B, Udayakumar M, Sheshshayee M S. 2014. Root traits and cellular level tolerance hold the key in maintaining higher spikelet fertility of rice under water limited conditions.,41(9): 930–939.

    Ramu V S, Swetha T N, Sheela S H, Babitha C K, Rohini S, Reddy M K, Tuteja N, Reddy C P, Prasad T G, Udayakumar M. 2016. Simultaneous expression of regulatory genes associated with specific drought-adaptive traits improves drought adaptation in peanut.,14(3): 1008–1020.

    Roitsch T. 1999. Source-sink regulation by sugar and stress.,2(3): 198–206.

    Sarvestani Z T, Pirdashti H, Sanavy S A M M, Balouchi H. 2008. Study of water stress effects in different growth stages on yield and yield components of different rice (L.) genotypes., 11(10): 1303–1309.

    Selote D S, Khanna-Chopra R. 2004. Drought-induced spikelet sterility is associated with an inefficient antioxidant defence in rice panicles, 121(3): 462–471.

    Serraj R, Sinclair T R, Purcell L C. 1999. Symbiotic N2?xation response to drought., 50(331): 143–155.

    Sheshshayee M S, Bindumadhava H, Shankar A G, Prasad T G, Udayakumar M. 2003. Breeding strategies to exploit water use efficiency for crop improvement., 30(2): 253–268.

    Sheshshayee M S, Vijayaraghavareddy P, Sreevathsa R, Rajendrareddy S, Arakesh S, Bharti P, Dharmappa P, Soolanayakanahally R. 2018. Introgression of physiological traits for a comprehensive improvement of drought adaptation in crop plants.,10(6): 92.

    Shi W J, Li X, Schmidt R C, Struik P C, Yin X Y, Jagadish S V K. 2018. Pollen germination andfertilization in response to high-temperature during flowering in hybrid and inbred rice.,41(6): 1287–1297.

    Turk K J, Hall A E, Asbell C W. 1980. Drought adaptation of cowpea: I. Influence of drought on seed yield.,72(3): 413–420.

    Udayakumar M, Rao R C N, Wright G C, Ramaswamy G C, Ashok RS, Gangadhar G C, Hussain I S A. 1998. Measurement of transpiration efficiency in field conditions., 1: 69–75.

    Uphoff N, KassamA, Harwood R. 2010. SRI as a methodology for raising crop and water productivity: Productive adaptations in rice agronomy and irrigation water management,9(1): 3–11.

    Vijayaraghavareddy P, Vanitha P, Vemanna R, Sheshshayee M S, Makarla U. 2017. Quantification of membrane damage/cell death using Evan’s blue staining technique.,7(16): e2519.

    Yang J C, Huang D F, Duan H, Tan G L, Zhang J H. 2009. Alternate wetting and moderate soil drying increases grain yield and reduces cadmium accumulation in rice grains.,89(10): 1728–1736.

    Yang X L, Wang B F, Chen L, Li P, Cao C G. 2019. The different influences of drought stress at the flowering stage on rice physiological traits, grain yield, and quality., 9(1): 3742.

    Zeigler R S, Puckridge D W. 1995. Improving sustainable productivity in rice-based rainfed lowland systems of South and Southeast Asia.,35(3): 307–324.

    Zhang J Y, Broeckling C D, Sumner L W, Wang Z Y. 2007. Heterologous expression of two Medicago truncatula putative ERF transcription factor genes, WXP1 and WXP2, inled to increased leaf wax accumulation and improved drought tolerance, but differential response in freezing tolerance.,64(3): 265–278.

    12 June 2019;

    24 October 2019

    Sheshshayee Sreeman(msshesh1@uasbangalore.edu.in)

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

    Peer review under responsibility of China National Rice Research Institute

    http://dx.doi.org/10.1016/j.rsci.2020.05.009

    (Managing Editor: Fang Hongmin)

    亚洲,欧美精品.| 亚洲自拍偷在线| 免费在线观看成人毛片| ponron亚洲| 国产精品美女特级片免费视频播放器 | 欧美日韩精品网址| 在线观看舔阴道视频| 香蕉久久夜色| 淫秽高清视频在线观看| 午夜日韩欧美国产| 亚洲熟妇中文字幕五十中出| 久久国产精品人妻蜜桃| 午夜免费观看网址| 观看免费一级毛片| 波多野结衣高清作品| 高潮久久久久久久久久久不卡| 天天躁狠狠躁夜夜躁狠狠躁| 老司机午夜十八禁免费视频| 女人高潮潮喷娇喘18禁视频| 99久久国产精品久久久| 国产一区二区三区视频了| 最近在线观看免费完整版| 国产高清视频在线播放一区| 亚洲成人久久爱视频| 精品人妻1区二区| 国产区一区二久久| 真人一进一出gif抽搐免费| 亚洲人成网站在线播放欧美日韩| cao死你这个sao货| 在线观看免费视频日本深夜| 男插女下体视频免费在线播放| 又黄又爽又免费观看的视频| 色精品久久人妻99蜜桃| 女同久久另类99精品国产91| 亚洲五月婷婷丁香| 老司机靠b影院| 亚洲av中文字字幕乱码综合| 亚洲人成77777在线视频| 特级一级黄色大片| 最近最新免费中文字幕在线| 又紧又爽又黄一区二区| 99在线视频只有这里精品首页| 久久久久久免费高清国产稀缺| 99riav亚洲国产免费| 久久婷婷成人综合色麻豆| 他把我摸到了高潮在线观看| 手机成人av网站| 久久久国产成人免费| 无遮挡黄片免费观看| 国产片内射在线| 国产精品日韩av在线免费观看| 国产精品电影一区二区三区| 黄频高清免费视频| 国产精品亚洲一级av第二区| 99在线视频只有这里精品首页| 两个人的视频大全免费| 我的老师免费观看完整版| 亚洲人成网站在线播放欧美日韩| 国产99久久九九免费精品| 老司机福利观看| 一级作爱视频免费观看| 国产精品久久久久久人妻精品电影| 一区二区三区激情视频| 桃色一区二区三区在线观看| 亚洲美女视频黄频| 成人欧美大片| 99国产精品一区二区蜜桃av| 999精品在线视频| 欧美3d第一页| 国产欧美日韩一区二区三| 亚洲中文日韩欧美视频| 欧美黑人巨大hd| 18禁观看日本| www.精华液| 久久久久久久午夜电影| 757午夜福利合集在线观看| 真人一进一出gif抽搐免费| 亚洲成人精品中文字幕电影| 亚洲第一电影网av| 欧美日韩中文字幕国产精品一区二区三区| 男男h啪啪无遮挡| 人成视频在线观看免费观看| 黄色 视频免费看| 国产v大片淫在线免费观看| 国产精品一区二区三区四区久久| 一区二区三区国产精品乱码| 国产精品久久久久久人妻精品电影| 国产精品亚洲一级av第二区| 在线看三级毛片| 久久婷婷人人爽人人干人人爱| 人妻丰满熟妇av一区二区三区| 亚洲男人天堂网一区| 国产av一区在线观看免费| 国产av不卡久久| 午夜激情av网站| 精品无人区乱码1区二区| 免费av毛片视频| 天堂av国产一区二区熟女人妻 | 国产精品精品国产色婷婷| 小说图片视频综合网站| 夜夜夜夜夜久久久久| 欧美黑人精品巨大| 在线国产一区二区在线| 嫁个100分男人电影在线观看| 日韩欧美国产一区二区入口| 91大片在线观看| 国产成人精品无人区| 久久久久亚洲av毛片大全| 最近最新中文字幕大全电影3| 午夜福利免费观看在线| 在线a可以看的网站| 97超级碰碰碰精品色视频在线观看| 欧美黑人精品巨大| 一本综合久久免费| 香蕉国产在线看| 婷婷六月久久综合丁香| 国产精品亚洲美女久久久| 亚洲欧美日韩高清专用| 黄色女人牲交| 日本熟妇午夜| 亚洲欧美激情综合另类| 美女高潮喷水抽搐中文字幕| 国产精品国产高清国产av| 国产日本99.免费观看| 欧美黑人欧美精品刺激| 白带黄色成豆腐渣| 久久精品国产99精品国产亚洲性色| 国产精品香港三级国产av潘金莲| 国产99久久九九免费精品| 一边摸一边抽搐一进一小说| 国产一区在线观看成人免费| 亚洲成av人片在线播放无| 最近在线观看免费完整版| 91麻豆精品激情在线观看国产| 特级一级黄色大片| 日本撒尿小便嘘嘘汇集6| 哪里可以看免费的av片| 在线观看免费视频日本深夜| 亚洲男人的天堂狠狠| 51午夜福利影视在线观看| 女人高潮潮喷娇喘18禁视频| 日韩国内少妇激情av| 少妇熟女aⅴ在线视频| 久久精品亚洲精品国产色婷小说| 欧美最黄视频在线播放免费| 国产成人av激情在线播放| 91大片在线观看| 日本一二三区视频观看| 亚洲男人天堂网一区| 久久久久久免费高清国产稀缺| √禁漫天堂资源中文www| 亚洲黑人精品在线| 男女下面进入的视频免费午夜| 亚洲国产欧美网| 无人区码免费观看不卡| 久久精品影院6| 一区二区三区高清视频在线| 黑人操中国人逼视频| 亚洲一区中文字幕在线| 激情在线观看视频在线高清| 亚洲国产精品合色在线| 757午夜福利合集在线观看| 天天一区二区日本电影三级| 国产伦一二天堂av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久久久国内视频| 久久中文看片网| 老司机靠b影院| 波多野结衣巨乳人妻| 国产区一区二久久| 亚洲人成77777在线视频| 久久午夜亚洲精品久久| 中文资源天堂在线| 在线播放国产精品三级| 美女扒开内裤让男人捅视频| 国产精品一区二区三区四区久久| 又黄又粗又硬又大视频| 国产av又大| 村上凉子中文字幕在线| 亚洲精品国产一区二区精华液| 亚洲男人天堂网一区| 亚洲精品在线美女| 国产乱人伦免费视频| 久久久久久久久中文| 国产激情久久老熟女| 精品无人区乱码1区二区| 制服诱惑二区| 少妇裸体淫交视频免费看高清 | 韩国av一区二区三区四区| 天天躁狠狠躁夜夜躁狠狠躁| 一边摸一边抽搐一进一小说| 真人做人爱边吃奶动态| 久久精品夜夜夜夜夜久久蜜豆 | 一级毛片精品| 日韩精品免费视频一区二区三区| 国产精品久久视频播放| 亚洲午夜理论影院| 1024香蕉在线观看| 亚洲一区高清亚洲精品| 亚洲成人免费电影在线观看| 最近视频中文字幕2019在线8| 日本一区二区免费在线视频| 欧美精品亚洲一区二区| 90打野战视频偷拍视频| 伦理电影免费视频| 2021天堂中文幕一二区在线观| 国产黄色小视频在线观看| 婷婷丁香在线五月| 成人av一区二区三区在线看| 欧美在线一区亚洲| 99国产精品99久久久久| 在线永久观看黄色视频| 在线看三级毛片| 成人高潮视频无遮挡免费网站| 在线观看一区二区三区| 国产麻豆成人av免费视频| 757午夜福利合集在线观看| 91麻豆精品激情在线观看国产| 禁无遮挡网站| 不卡av一区二区三区| 草草在线视频免费看| 好男人在线观看高清免费视频| 看片在线看免费视频| 久久国产乱子伦精品免费另类| 大型av网站在线播放| 黄色a级毛片大全视频| av中文乱码字幕在线| 亚洲国产欧美一区二区综合| 级片在线观看| 无遮挡黄片免费观看| 日韩高清综合在线| 国产亚洲精品一区二区www| 国产免费av片在线观看野外av| 99久久无色码亚洲精品果冻| 99久久99久久久精品蜜桃| 久久这里只有精品中国| x7x7x7水蜜桃| 亚洲国产高清在线一区二区三| 麻豆国产97在线/欧美 | 嫩草影视91久久| 久久久久久久精品吃奶| 色综合婷婷激情| 日本 av在线| 一级作爱视频免费观看| 久久国产精品影院| 免费观看人在逋| 欧美久久黑人一区二区| 欧美在线黄色| 亚洲人成伊人成综合网2020| 在线观看www视频免费| 久久中文字幕一级| 啦啦啦观看免费观看视频高清| 人人妻,人人澡人人爽秒播| 亚洲精华国产精华精| 人人妻人人看人人澡| 久久这里只有精品中国| 又爽又黄无遮挡网站| 欧美三级亚洲精品| 免费搜索国产男女视频| 又黄又粗又硬又大视频| 亚洲无线在线观看| 啦啦啦观看免费观看视频高清| 美女大奶头视频| 午夜精品在线福利| 一级黄色大片毛片| 亚洲七黄色美女视频| 日日夜夜操网爽| 99久久精品热视频| 午夜福利成人在线免费观看| 欧美日本亚洲视频在线播放| 免费在线观看黄色视频的| 亚洲精品一卡2卡三卡4卡5卡| 最新在线观看一区二区三区| 国产精品久久久久久久电影 | 欧美+亚洲+日韩+国产| 国产午夜精品论理片| 国产成年人精品一区二区| 99久久综合精品五月天人人| 精品人妻1区二区| 成人手机av| 超碰成人久久| 国产一区二区三区在线臀色熟女| 日本一本二区三区精品| 国产精品一区二区精品视频观看| 亚洲国产精品合色在线| 久久久久久九九精品二区国产 | 可以免费在线观看a视频的电影网站| 人妻久久中文字幕网| 在线国产一区二区在线| 在线a可以看的网站| 亚洲av成人不卡在线观看播放网| 国产免费av片在线观看野外av| 国产三级黄色录像| 国产亚洲欧美在线一区二区| 久久人妻福利社区极品人妻图片| 国产成人精品久久二区二区91| 99久久无色码亚洲精品果冻| 99热6这里只有精品| 国产人伦9x9x在线观看| 亚洲自偷自拍图片 自拍| 欧美日本亚洲视频在线播放| 亚洲精品国产精品久久久不卡| 国产亚洲av嫩草精品影院| 国产久久久一区二区三区| 岛国在线观看网站| 丁香欧美五月| 香蕉丝袜av| av超薄肉色丝袜交足视频| 三级毛片av免费| 特大巨黑吊av在线直播| 日日摸夜夜添夜夜添小说| 久久精品91蜜桃| 亚洲欧美日韩无卡精品| 99国产综合亚洲精品| 人妻丰满熟妇av一区二区三区| 巨乳人妻的诱惑在线观看| 亚洲人成伊人成综合网2020| 午夜福利在线在线| 亚洲人成电影免费在线| 久久久久久人人人人人| 99在线视频只有这里精品首页| 无限看片的www在线观看| 在线永久观看黄色视频| 国产精品久久视频播放| 久久久久精品国产欧美久久久| 不卡av一区二区三区| 男女下面进入的视频免费午夜| 亚洲欧美一区二区三区黑人| 日韩欧美在线二视频| 久久久久久免费高清国产稀缺| www.999成人在线观看| 在线观看日韩欧美| 国产亚洲精品一区二区www| 两个人免费观看高清视频| 成人手机av| 99久久99久久久精品蜜桃| av在线播放免费不卡| 久久九九热精品免费| 欧美日韩国产亚洲二区| 日韩三级视频一区二区三区| 亚洲 欧美一区二区三区| 亚洲精品色激情综合| 天堂影院成人在线观看| 久久精品国产清高在天天线| 亚洲国产精品久久男人天堂| 久久99热这里只有精品18| 亚洲av电影不卡..在线观看| 久久精品aⅴ一区二区三区四区| 精品无人区乱码1区二区| 亚洲专区中文字幕在线| 精品一区二区三区四区五区乱码| 久久婷婷成人综合色麻豆| 91国产中文字幕| 国产亚洲av高清不卡| 亚洲五月天丁香| 亚洲一区中文字幕在线| а√天堂www在线а√下载| 成人永久免费在线观看视频| 一夜夜www| 国产伦在线观看视频一区| 最近视频中文字幕2019在线8| 亚洲天堂国产精品一区在线| 人人妻人人澡欧美一区二区| 免费人成视频x8x8入口观看| 伦理电影免费视频| 波多野结衣高清无吗| 在线观看午夜福利视频| 久久久久国产精品人妻aⅴ院| 免费高清视频大片| 日本在线视频免费播放| 日韩成人在线观看一区二区三区| 亚洲成av人片免费观看| 国产精品免费视频内射| 制服丝袜大香蕉在线| 久久久久久久久中文| 久久亚洲精品不卡| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品综合一区在线观看 | 成人手机av| 老司机在亚洲福利影院| 精品国产乱码久久久久久男人| 亚洲男人的天堂狠狠| 大型黄色视频在线免费观看| 岛国在线免费视频观看| 成人一区二区视频在线观看| 美女午夜性视频免费| 啪啪无遮挡十八禁网站| 国产激情欧美一区二区| 国产不卡一卡二| 亚洲人成电影免费在线| 女同久久另类99精品国产91| 国内毛片毛片毛片毛片毛片| 欧美中文综合在线视频| av天堂在线播放| 精品国产美女av久久久久小说| 波多野结衣巨乳人妻| 久久99热这里只有精品18| 国产97色在线日韩免费| 亚洲人成网站在线播放欧美日韩| 最近视频中文字幕2019在线8| 久久午夜综合久久蜜桃| 婷婷六月久久综合丁香| 久久精品国产99精品国产亚洲性色| 精品国产超薄肉色丝袜足j| 亚洲黑人精品在线| 午夜日韩欧美国产| 亚洲18禁久久av| 国产69精品久久久久777片 | 熟妇人妻久久中文字幕3abv| 国产高清视频在线观看网站| 国产单亲对白刺激| 身体一侧抽搐| 老鸭窝网址在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲黑人精品在线| 国产成人精品久久二区二区免费| 五月伊人婷婷丁香| 午夜免费激情av| 成人欧美大片| 国产主播在线观看一区二区| 夜夜夜夜夜久久久久| 亚洲 国产 在线| 禁无遮挡网站| 国产精品99久久99久久久不卡| 99精品欧美一区二区三区四区| 国产熟女午夜一区二区三区| 黑人操中国人逼视频| 亚洲成人国产一区在线观看| 久久久久国产一级毛片高清牌| 亚洲一区二区三区色噜噜| 99国产精品一区二区蜜桃av| 国产亚洲精品av在线| 哪里可以看免费的av片| 日韩欧美在线二视频| 国产男靠女视频免费网站| tocl精华| 精品不卡国产一区二区三区| 日本免费a在线| 天堂av国产一区二区熟女人妻 | 午夜激情av网站| 黄色女人牲交| 最新在线观看一区二区三区| 91字幕亚洲| 禁无遮挡网站| 久热爱精品视频在线9| 国内毛片毛片毛片毛片毛片| 久久久精品大字幕| 久久精品aⅴ一区二区三区四区| 国产成人一区二区三区免费视频网站| 久久久久国产精品人妻aⅴ院| АⅤ资源中文在线天堂| 91大片在线观看| 在线观看一区二区三区| 国产成人aa在线观看| 妹子高潮喷水视频| 亚洲一码二码三码区别大吗| 亚洲午夜理论影院| a级毛片a级免费在线| 啦啦啦观看免费观看视频高清| 国产不卡一卡二| 免费看日本二区| 久久久久亚洲av毛片大全| av福利片在线| 精品国产亚洲在线| 丝袜人妻中文字幕| 天天一区二区日本电影三级| 午夜亚洲福利在线播放| 一本综合久久免费| 天天添夜夜摸| 中文资源天堂在线| 久久久久久久午夜电影| 午夜精品一区二区三区免费看| 五月伊人婷婷丁香| 夜夜夜夜夜久久久久| 日韩精品青青久久久久久| 久久午夜亚洲精品久久| 精品福利观看| а√天堂www在线а√下载| 久久久久免费精品人妻一区二区| 久久天堂一区二区三区四区| 欧美黄色淫秽网站| 又粗又爽又猛毛片免费看| 欧美乱妇无乱码| 国产精品99久久99久久久不卡| 国产精品精品国产色婷婷| 亚洲午夜理论影院| 久久中文字幕人妻熟女| 一级毛片高清免费大全| 亚洲精品美女久久av网站| 国产高清videossex| 99re在线观看精品视频| 性欧美人与动物交配| 国产精品影院久久| 在线观看免费午夜福利视频| 国产午夜福利久久久久久| 亚洲精品久久国产高清桃花| 国产精品免费视频内射| 麻豆av在线久日| 午夜激情av网站| 中亚洲国语对白在线视频| 欧美日韩精品网址| 国产精品久久久久久人妻精品电影| 亚洲国产欧美网| 一个人免费在线观看的高清视频| 50天的宝宝边吃奶边哭怎么回事| 久久这里只有精品中国| 亚洲第一电影网av| 国产免费av片在线观看野外av| 国产成人欧美在线观看| 亚洲精品一区av在线观看| 国产69精品久久久久777片 | 身体一侧抽搐| 午夜精品在线福利| 黄色 视频免费看| 桃红色精品国产亚洲av| 欧美三级亚洲精品| 亚洲成av人片免费观看| 久久久久久免费高清国产稀缺| 免费看十八禁软件| 精品不卡国产一区二区三区| 18美女黄网站色大片免费观看| 国产精品国产高清国产av| 日韩av在线大香蕉| 国产精品国产高清国产av| 国产精品久久久久久人妻精品电影| 中文资源天堂在线| 淫秽高清视频在线观看| 香蕉久久夜色| 母亲3免费完整高清在线观看| 亚洲av五月六月丁香网| 丝袜人妻中文字幕| 99久久久亚洲精品蜜臀av| 97超级碰碰碰精品色视频在线观看| 久久久国产欧美日韩av| 久久伊人香网站| 国产精品一及| 最近最新中文字幕大全电影3| 国产欧美日韩一区二区三| 中文在线观看免费www的网站 | 色综合欧美亚洲国产小说| 国产av在哪里看| 又爽又黄无遮挡网站| 亚洲国产看品久久| www.精华液| 久久精品国产亚洲av香蕉五月| 特大巨黑吊av在线直播| 国产一区二区在线观看日韩 | 中文亚洲av片在线观看爽| 免费看十八禁软件| 熟女电影av网| 禁无遮挡网站| 小说图片视频综合网站| 男插女下体视频免费在线播放| 国产乱人伦免费视频| 日韩av在线大香蕉| 757午夜福利合集在线观看| 亚洲一码二码三码区别大吗| 免费无遮挡裸体视频| 免费搜索国产男女视频| 变态另类成人亚洲欧美熟女| 久久久久九九精品影院| 亚洲一区二区三区色噜噜| 一进一出抽搐动态| 2021天堂中文幕一二区在线观| 成年免费大片在线观看| 美女午夜性视频免费| 国产精品永久免费网站| 深夜精品福利| 黄色毛片三级朝国网站| 黄色视频,在线免费观看| 757午夜福利合集在线观看| 亚洲精品久久国产高清桃花| 久久精品人妻少妇| 99久久精品国产亚洲精品| 中文字幕人成人乱码亚洲影| 黄色丝袜av网址大全| 十八禁网站免费在线| 亚洲av中文字字幕乱码综合| 国内毛片毛片毛片毛片毛片| 亚洲最大成人中文| 久久久国产成人免费| 草草在线视频免费看| 欧美乱码精品一区二区三区| 51午夜福利影视在线观看| 国产黄片美女视频| 国产精华一区二区三区| 我的老师免费观看完整版| 欧美人与性动交α欧美精品济南到| 国产黄片美女视频| 日韩欧美免费精品| 亚洲av五月六月丁香网| 日本一二三区视频观看| 两个人的视频大全免费| 国产麻豆成人av免费视频| 欧美又色又爽又黄视频| 久久这里只有精品19| 亚洲成人精品中文字幕电影| 亚洲av电影不卡..在线观看| 久久中文看片网| 波多野结衣巨乳人妻| 一夜夜www| 一边摸一边做爽爽视频免费| 亚洲人成电影免费在线| 美女大奶头视频| 最新美女视频免费是黄的| 亚洲精品粉嫩美女一区| 国产成人精品久久二区二区免费| а√天堂www在线а√下载| 久久精品人妻少妇| 欧美久久黑人一区二区| 99国产精品一区二区蜜桃av| 午夜福利在线在线| 日韩大尺度精品在线看网址| 午夜福利在线观看吧| 搡老岳熟女国产|