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

    Field identif ication of morphological and physiological traits in two special mutants with strong tolerance and high sensitivity to drought stress in upland rice (Oryza sativa L.)

    2019-05-10 06:13:34HUANGMinXUYuhuiWANGHuaqi
    Journal of Integrative Agriculture 2019年5期

    HUANG Min, XU Yu-hui, WANG Hua-qi

    College of Agronomy and Biotechnology, China Agricultural University, Beijing 100193, P.R.China

    Abstract The two mutants idr1-1 and 297-28, which were obtained from the radiation mutation of HD297 and IAPAR9, were used as experimental materials in this study for a 2-year (2012 and 2013) experiment about f ield drought resistance identif ication in Beijing, China. Key agronomic traits and water-related physiological indexes were observed and measured, including the leaf anti-dead level (LADL), days to heading, plant height, setting percentage, aboveground biomass, leaf water potential (LWP), net photosynthetic rate (P n) and transpiration rate. The results showed that the mutant idr1-1 that was under drought stress (DS) conditions for 2 years had the highest LADL grades (1.3 and 2.0) among all the materials, and they were 2-3 grades stronger than the wild-type IAPAR9 with an average that was 21.4% higher for the setting percentage than the wild type. Compared with the IAPAR9 for the 2-year average delay in the days to heading and the reduction rates in the plant height, setting percentage, and aboveground biomass under DS compared with the well-watered (WW) treatment, idr1-1 showed 3.2% less delay and 19.1, 16.4, and 6.1% less reduction, respectively. The idr1-1 in the LWP always exhibited the highest performance among all the materials. The P n of idr1-1 under severe and mild DS comparing with that under WW was slightly decreased and even slightly increased, respectively, leading to an average reduction rate of only 0.92%, which was 26.93% less than that of IAPAR9. Under the severe DS, idr1-1 still showed the highest value of 16.88 μmol CO2 m-2 s-1 among all the materials and was signif icantly higher than that of IAPAR9 (11.66 μmol CO2 m-2 s-1). Furthermore, only idr1-1 had the increased and the highest transpiration rate values (7.6 and 6.04 mmol H2O m-2 s-1) under both mild and severe DS compared with the values under WW, when the transpiration rate of all the other materials signif icantly decreased. By contrast, the 297-28 in terms of the LADL grade under DS was the lowest (7.0), and it was four grades weaker than its wildtype HD297 and even one grade weaker than the drought-sensitive paddy rice SN265. For the 2-year average reduction rates in aboveground biomass and plant heights under DS compared with those under the WW, 297-28 was 31.6 and 31.8% higher than HD297, respectively. Meanwhile, 297-28 showed the worst performance for the LWP, P n, and transpiration rate. These results suggest that idr1-1 might be a superior drought tolerant mutant of upland rice found in China. It has a strong ability to maintain and even enhance leaf transpiration while maintaining a high plant water potential under DS, thus supporting a high P n and alleviating the delay in agronomic trait development and yield loss effectively. 297-28 is a much more highly drought-sensitive mutant that is even more sensitive than paddy rice varieties.The two mutants could be used as drought tolerance controls for rice germplasm identif ication and the drought resistant mechanism studies in the future. idr1-1 is also suitable for breeding drought-tolerant and lodging-resistant high-yield rice varieties.

    Keywords: rice mutants, morphological and physiological traits, drought resistance, field identif ication

    1. lntroduction

    Rice is one of the most important grain crops in Asia, and it provides approximately one-third of the caloric intake in the region. According to the FAO, 140.4 million hectares of harvested area produced 667.9 million tons of rice in 2016 (FAO 2017). However, a water shortage has been an important restriction for rice production in many rainfed areas of the world (Pandey et al. 2007). In China, the agriculture industry consumes approximately 70.4% of the water resources, and the proportion of that water used by rice production is nearly 70% (Zhang 2007). Most rice production in China depends on irrigation; however, a less-than-abundant water supply still presents a long-term threat. Water shortages will probably remain a challenge to the sustainable production of rice based on recent climate predictions (Wassmann et al. 2009), which suggest that the intensity and frequency of drought will be aggravated.

    The breeding of drought-tolerant cultivars is an effective and economical way to address this problem, and some of these cultivars have been released, including Morobérékan, IRAT109, B96-3, and IAPAR9. However, the conventional breeding progress towards improved drought tolerance in rice cultivars has been relatively slow (Cooper et al. 1999; Wade et al. 1999; Bernier et al. 2008). One of the constraints is the narrow genetic basis because of the germplasm shortage. The use of radiation-induced mutation is an effective way to enhance the genetic variability in many species (Ahloowalia and Maluszynski 2001; Mostafa et al. 2015; Mondal et al. 2017). In rice, many novel germplasm resources were generated by gamma irradiation. In a functional genomics study, Kong et al. (2014) reported the morphological and physicochemical properties of two starch rice mutants, thereby supporting new genotypes for food processing, and they identif ied defective genes. Aromatic rice with high iron was also obtained by gamma radiation, and this germplasm resource was used to develop a new cultivar (Tran and Ho 2017). In addition, the salt resistance of gamma-induced rice was enhanced using a proper radiation dosage (Shereen et al. 2009). Droughtresistance properties have been demonstrated in sorghum (Sihono 2010; Human et al. 2012), soybeans (Moussa 2011), and sugar beets (Sen and Alikamanoglu 2012), but few studies have focused on rice. In China, many donor parents responsible for the drought tolerance of the released cultivars were derived from IRAT109, IAPAR9, or some landraces. Additionally, material def iciencies have limited gene cloning for drought resistance in rice. For example, the donor parents used to construct populations by positional cloning were mostly derived from CT9993, Swarna, Azucena, Apo, and some other landraces (Lanceras et al. 2004; Venuprasad et al. 2009).

    In the present study, we selected two pure lines of mutants derived from gamma-induced IAPAR9 and HD297 to identify their drought responses. Morphological traits, such as the leaf anti-dead level (LADL), days to heading, plant height, setting percentage and aboveground biomass, and physiological traits, such as the leaf water potential (LWP), net photosynthetic rate (Pn), and transpiration rate, were used to evaluate the drought tolerance under wellwatered (WW) and drought stress (DS) conditions. One mutant called idr1-1 conferred higher drought tolerance characteristics than IAPAR9, the wild type, whereas 297-28 showed greater susceptibility to drought than the HD297 wild type. These two mutants may provide superior germplasm for use in developing drought-resistant varieties and for studying their molecular mechanisms.

    2. Materials and methods

    2.1. Materials for the study

    Six genotypes were used in this study: HD297 (Oryza sativa L. ssp. japonica), 297-28, IAPAR9 (O. sativa L. ssp. indicajaponica intermediate type), idr1-1, SN265 (O. sativa L. ssp. japonica), and IRAT109 (O. sativa L. ssp. japonica). HD297 is an upland rice variety that is planted widely in North China. As a standard control for the national upland rice regional trials (NURRT), HD297 exhibits a good drought tolerance phenotype. The 297-28 is a mutant induced from HD297 by gamma-ray radiation. This mutant shows thinner stems, leaves, and roots and more tillers than HD297 under WW conditions, whereas under drought stress, the leaves wilt easily, leading to sharp decreases in the yield. IAPAR9 is a drought-tolerant upland rice introduced by the Brazilian Agricultural Research Corporation (EMBRAPA) in 1992. The idr1-1 is a semi-dwarf mutant derived from IAPAR9 by γ-ray radiation. It shows stronger drought resistance than the wild type. SN265 is a super high-yield rice variety with low drought resistance in North China, and IRAT109 is a famous international upland rice variety with strong drought resistance from the Cote d'Ivoire; they were the other two controls used in this study.

    2.2. Methods

    Planting methodsThe f ield experiments for the above six genotypes were conducted in the summer and autumn of 2012 and 2013 at Shangzhuang Experimental Station of China Agricultural University, Beijing, China (40°08′13.4′′N, 116°11′06.6′′E). The study f ield is sandy soil. Two treatments were applied, namely, DS and WW, with three replicates for each treatment. A randomized complete block design (RCBD) was applied to all the experimental plots. The commonly used local sowing method, direct sowing, was used to sow the rice; the applied base fertilizer consisted of 48 kg N ha-1, 120 kg P2O5ha-1, 100 kg K2O ha-1, 22.5 kg ZnSO4ha-1, and 30 kg FeSO4ha-1, and a top dressing of 45 kg N ha-1urea was applied at the tillering stage. In 2012, the plot areas of the WW and DS treatments were 4.5 m2(six rows per plot, plot size 2.5 m×1.8 m) and 3.0 m2(four rows per plot, plot size 2.5 m×1.2 m), respectively. In 2013, the plot areas of the WW and DS treatments were 10.92 m2(seven rows per plot, plot size 5.2 m×2.1 m) and 15.75 m2(seven rows per plot, plot size 7.5 m×2.1 m), respectively. Chemical and manual methods were used for weed control, and local and conventional management practices for rice cultivation were followed.

    Drought treatmentFor the WW treatment, irrigation was applied for the entire growth period to maintain a shallow layer of water on the surface and a wet state in the experimental f ield. For the DS treatment, natural drought was applied, in which time full irrigation was applied after sowing to ensure that all the tested materials had normal germinations, and then, from seedling emergence until the end of the experiment, the plant growth was completely dependent on natural precipitation without any irrigation. At the study site in Beijing, the periods from May to June (from the seedling to tillering stages) and from mid-August to early October (from the heading to maturing stages) occur mostly in the dry season when they receive scarce precipitation, and drought may occur between July and early August (from jointing to booting) during some years. During this time without precipitation, the soil moisture sensor TRIME-PICO TDR (IMKO, Ettlingen, Germany) was used to measure the soil water content (%, v/v) at the 0-20 cm layers of the soil surface every 3-5 days. According to our observations over many years, when the topsoil (depth 20 cm) is saturated, the soil volumetric water content is generally above 40 or 35%. When it is lower than 20%, the ground begins to crack, and some of the top rice leaves begin to die. When the soil volumetric water content is lower than 15%, the ground cracks become obvious, and most leaves show complete or partial death. When this content is less than 10%, the whole plant may die. Therefore, in general, volumetric water content ranges of 30-25%, 25-20%, 20-15%, 15-10%, and <10% in the topsoil (depth 20 cm) indicate soil drought stress as follows: “beginning (slight)”, “mild”, “moderate”, “severe”, and “extreme severe”, respectively. When the soil volumetric water content drops below approximately 25%, “mild” drought has begun in the soil. After that, the rice plants were subjected to 2 weeks or more of drought treatment, and their related traits were investigated.

    Trait investigationThe LWP and leaf photosynthesis parameters (LPPs) under diurnal change and the key time points were investigated. The LPPs, which include two major indexes, the Pnand the transpiration rate, were measured after 2 weeks of continuous drought during the heading and grain-f illing stages (~80-100 days after sowing, or DAS). The LPPs were measured using an LI-6400XT Portable Photosynthesis System (Li-Cor Environmental, Lincoln, NE, USA) with the following parameters: a reference cell CO2(Ca) concentration of 380-400 μmol CO2mol-1, a f low of 400 μmol s-1and an instant photosynthetic photon f lux density (PPFD) under natural light. The LWP was determined with a pressure chamber instrument, the Pressure Chamber PMS600 (PMS instruments, Albany, OR, USA). The sampling method was as follows: A plot showing representative growth of rice materials was selected for measurement, and plants with robust growth and the same morphological conditions were selected as candidate plants before the measurement. Candidate plants were marked for further measurements. For each rice material, a total of 3-5 plants were measured, and the measurements were conducted on the uppermost f irst fully expanded leaf for each plant. The daily change monitoring for the LPPs and LWP was performed at four time slots on the same sunny day, from 07:00-08:00 h, 09:00-10:00 h, 14:00-15:00 h, and 17:00-18:00 h. The key time-point measurement was only conducted on sunny days at 09:00-11:00 h (for LPPs) and at 14:00-16:00 h (for LWP); measurements for the WW and DS treatments were conducted on two adjacent days (e.g., measurement for WW was only on the 1st day and for DS, it was only on the 2nd day), and they were repeated three times (taking 6 days).

    The rice leaves become wilted or even died under drought stress. The leaf death performance among different genotypes ref lects the strength of the drought resistance. The LADL was evaluated by taking the ratio of the green leaf length to the total leaf length of each tested leaf by measurement (Zu et al. 2017) or visualization (Xu et al. 2018). The grading standards for the LADL in our lab was on a 1-9 resistance rating scale (1=highly resistant, 3=moderately resistant, 5=moderately sensitive, 7=sensitive, and 9=highly sensitive) (Xu et al. 2018). In this study, when the LPPs and LWP (~100 DAS) were measured, the LADL score was investigated by visual estimation for the DS treatment.

    The days to heading parameter was also recorded. Before each harvest, sampling and seed investigations were performed. In each plot, one row of plants with uniform growth was selected, and a 1-m-long sample segment was accurately measured for sampling. The plants in this sample segment were sampled by removal, along with the roots, and 10 plants were randomly selected to measure the plant height, setting percentage, aboveground biomass, and other indexes. The measurement standards are shown in Table 1.

    Statistical analysis of the study dataEach trait index was presented as the mean±SD, and they were statistically analyzed using SAS 9.1.3 (SAS Ins., Cary, NC, USA). The SAS general linear model (GLM) procedure was used to analyze the signif icance of the variations between the genotypes, between the treatments, and between the replicates (*, P<0.05;**, P<0.01). The DUNCAN method was used for multiple comparisons between varieties. The SAS TTEST procedure and Student's t-test (*, P<0.05;**, P<0.01) were used for comparisons between the WW and DS treatments.

    3. Results

    3.1. Droughtstress degrees and effects

    The drought occurrence degree in the tillage layer of the soil at the heading and grain-f illing stages is shown in Fig. 1. In 2012, the soil had a high water content due to abundant rainfall, and yet a rain-free period occurred from 100- 115 days after sowing (DAS). During this period, the volumetric water content in the soil surface layer ranged between 15 and 25%, and a low value of 15% was maintained for approximately 10 days, thereby causing a mild to moderate drought on the whole. In 2013, no effective precipitation occurred in 80-120 DAS, and the volumetric water content in the soil surface layer was always below 15% (indicating moderate drought) and mostly f luctuated around 10% (indicating severe drought), thereby indicating a generally severe drought. Those 2 years provided an ideal natural drought environment for this study. Herein, 2012 is also referred to as a mild drought and 2013 as a severe drought.

    3.2. Analyzing the agronomic traits of the mutants

    The results for the joint analysis of variance for the effects under the DS and WW treatments are shown in Table 2. The signif icant difference tests between the DS and WW treatments showed that in 2013, the differences in all the traits between the DS treatment and the WW treatment reached the extremely signif icant level or the signif icant level. Specif ically, under the DS treatment, the aboveground biomass, plant height, setting percentage, and LADL decreased extremely signif icant, and the days to heading were signif icantly delayed. In 2012, except for the setting percentage, the other traits under DS treatment were also signif icantly much weaker or signif icantly weaker than those under the WW treatment, with the same trend being observed in 2013. Only the setting percentage did not decrease, but instead it showed an increase, which occurred primarily because during the heading stage under the WW treatment in 2012, 297-28 had lodging problems, whereas the IRAT109 experienced selective damage from rice planthoppers. These problems led to abnormal (decreased) setting percentages in both genotypes under the WW treatment, with setting percentage values that were signif icantly even lower than those under the DS treatment. Compared with the observations in 2012, the degree of difference in each of the traits between the DS and WW treatments was signif icantly higher in 2013. The results for the 2 years showed that compared with the traits under the WW treatment, of all the traits, the aboveground biomass and LADL showed the most sensitivity in responding to drought under the DS treatment, with the reduction rates for both traits being the greatest. Specif ically, the aboveground biomass showed reduction rates of 33.2 and 44.0% in 2012 and 2013, respectively; the LADL showed reduction rates of 34.4% (3.1/9) and 38.9% (3.5/9) in 2012 and 2013, respectively. The days to heading were relatively stable, and they were delayed by 4.3 and 6.4% in 2012 and 2013, respectively. This result can be explained by the fact that during the 2-year experiment (especially in 2013), the f ield drought primarily occurred during the middle and late stages of plant growth and development, which led to damaged functional leaves in the middle and upper parts and directly reduced the biomass and grain yield, whereas the drought before the heading stage was relatively mild and thus had a small impact on the days to heading.

    Table 1 The list of measurement methods for the agronomic traits

    Fig. 1 Soil water moisture content in this study. DS, drought stress. LPPs, leaf photosynthesis parameters; LWP, leaf water potential.

    Table 2 shows that a signif icant difference is present in each of the various traits between the genotypes under the WW and DS treatments during the 2 years. Hence, a comparative analysis of the pros and cons among the genotypes for each trait is possible (Table 3).

    Leaf anti-d ead level (LADL, Table 3)Under the DS treatment, the idr1-1 had the lowest 2-year average LADL value (1.7), and 297-28 had the highest value (7.0). The LADL value of each genotype remained basically the same over the 2 years, and the differences between the genotypes reached a signif icant level. The LADL value of idr1-1 (1.3 and 2.0) was signif icantly lower than that of IAPAR9 (4.0 and 5.0), with a decrease by two to three levels. By contrast, 297-28 exhibited an LADL value of 7.0 in both years, which was signif icantly higher than that of HD297 (2.3 and 3.0) with an increase of four levels; the LADL value of 297-28 was even one level higher than that of the rice variety SN265, exhibiting the characteristics of higher drought sensitivity than common rice varieties.

    Days to heading (Table 3)Although drought primarily occurred during the grain-f illing and maturing stages, the days to heading for each of the genotype was nevertheless delayed. Compared with that under WW treatment, the days to heading under the DS treatment showed a delay. The degree of delay in days to heading in 2013 was signif icantly greater than it was in 2012 for all the genotypes (except 297-28). For the 2-year average delay degree, the result for the idr1-1 (4.0%) was signif icantly less than that of IAPAR9 (7.2%), with a decrease of 3.2%; it was also less than the delay for IRAT109 (6.2%). The HD297 genotype showed the lowest delay in the average days to heading, suggesting a high drought resistance ability, but the days to heading delay of 297-28 was quite similar to that of the wild type because it has a longer vegetative period and luxuriant vegetative growth under WW treatment. SN265 (8.8%) had a signif icantly higher number of delayed average days to heading than the other genotypes. The results showed that after drought stress treatment, the delay degree in the growth and development for drought-sensitive genotypes was generally greater than that for drought-tolerant genotypes.

    Plant height (Table 3)Under the WW environment, the plant height of idr1-1 (57.7 cm) was signif icantly lower than that of IAPAR9 (104.0 cm), indicating that it was a dwarf mutant. Compared with that under the WW treatment, in 2013, the plant height of idr1-1 under the DS treatment decreased signif icantly, but the reduction rate was signif icantly lower than that of IAPAR9. By contrast, the plant height reduction rate for 297-28 was much higher than that of HD297 and it was close to that of SN265, which had the highest reduction rate. In 2012, except for IRAT109, idr1-1 and HD297, which did not have signif icantly decreased plant heights, all the genotypes had decreased plant heights at the extremely signif icant level or at the signif icant level. For the 2-year average of plant height decreases, the reduction rate in the plant height for idr1-1 was 19.1% less than that of IAPAR9, and for 297-28, it was 31.8% more than that of HD297.

    Table 2 Comparison of agronomic traits between the well-watered (WW) and drought stress (DS) treatments using a joint analysis of variance

    Table 3 Analysis of the agronomic traits of the genotypes tested in this study1)

    Setting percentage (Table 3)In 2013, except for 297-28, which had a lower setting percentage under the WW treatment than under the DS treatment, all the genotypes had a signif icantly lower setting percentage under the DS treatment than under the WW treatment. The reduction rate of the setting percentage for idr1-1 reached an extremely signif icant level but was less than half of the reduction rate for IAPAR9. SN265 showed the greatest reduction in the setting percentage; its reduction rate was 60.0% more than that of IRAT109. In 2012, except for 297-28 and IRAT109, all the genotypes experienced a reduction in the setting percentage to varying degrees under DS compared to WW, but the reduction was not signif icant; idr1-1 even had a slightly increased setting percentage, which indicated that the mild drought stress did not signif icantly affect the setting percentage. Based on the average from 2-year data, idr1-1 showed a reduced setting percentage, which was 16.4% less than that of IAPAR9. During both years, as the drought stress intensity increased, the setting percentages for the genotypes decreased by various degrees. Under the DS treatment, the setting percentage of SN265 decreased from 88.5% in 2012 to 33.1% in 2013; it was the largest reduction observed (up to 62.6%) among all the genotypes. IAPAR9 also showed a signif icant reduction (21.0%). By contrast, idr1-1 showed only a marginal decrease (by 2.5%), which was close to that observed for IRAT109 (slight decrease of 1.3%). The setting percentage of 297-28 decreased slightly, but that of HD297 increased instead. These results indicate that there are signif icant differences between the genotypes in the setting percentage traits in terms of the ability to adapt to drought.

    Aboveground biomass (Table 3)During both years of this experiment, all the genotypes showed decreased aboveground biomasses to varying degrees under the DS treatment compared with those under the WW treatments, and these decreases were extremely signif icant in 2013. When the rate of reduction was averaged for the 2 years, HD297 had the smallest average reduction rate (18.0%), and SN265 showed the largest (52.2%). For idr1-1, its 2-year average reduction rate was 19.9%, which was 6.1% lower than that of IAPAR9 (26.0%) and 25.0% lower than that of IRAT109. For 297-28, the 2-year average reduction rate was 49.6%, which was 31.6% higher than that of HD297 and close to that of SN265. These results suggest that the drought resistance capacity of idr1-1 was slightly stronger than that of IAPAR9 and signif icantly stronger than that of IRAT109, whereas the capacity of 297-28 was signif icantly weaker than that of HD297 and comparable to that of SN265.

    3.3. Leaf water potential and photosynthetic physiology of mutants

    Daily changes in the leaf water potential and P n(1) Daily change in the leaf water potential (LWP). The level of LWP (MPa) directly ref lects the water content status of the plant body. The lower the value, the less free water there is in the plant body, which implies more serious drought stress. As shown in Figs. 2 and 3, under the DS and WW treatments, the highest LWP for each genotype occurred before sunrise, during which the differences among the genotypes were the smallest. After sunrise, with the increased light intensity and temperature, the LWPs rapidly decreased from 07:00-08:00 h to 09:00-10:00 h, during which the differences among the genotypes began to occur. The LWPs reached the lowest for all the genotypes at the 14:00-15:00 h, during which the difference among the genotypes was the greatest. Later, the LWPs gradually rose back from the 17:00-18:00 h to sunset, during which the differences among the genotypes experienced a decrease. At each time point, the LWP under the DS treatment was signif icantly lower than that under the WW treatment for the same genotype. The difference in LWPs among different genotypes under the DS treatment was signif icantly higher than that under the WW treatment; the difference reached the largest at the 14:00-15:00 h, and there were still signif icant differences among the genotypes at the 17:00-18:00 h under the DS treatment.

    Fig. 2 Daily change in light intensity and temperature. PPFD, photosynthetic photon f lux density.

    Fig. 3 Daily changes in the leaf water potential (LWP). The measurement day for the well-watered (WW) treatment was 08-22-2013, and for the drought stress (DS) treatment, it was 08-23-2013. The thick solid line represents a genotype under the DS treatment, whereas the thin dashed line represents a genotype under the WW treatment.

    For idr1-1, the LWPs for the entire day under both DS and WW treatments were higher than those of IAPAR9. Unlike the results under the WW treatment, the differences between the two genotypes under the DS treatment were obviously widened starting from early in the morning, reaching a maximum at the 09:00-10:00 h, and it was still present at the 17:00-18:00 h. By contrast, the difference between the two genotypes under the WW treatment only begin to appear after the 09:00-10:00 h, reaching its maximum until 14:00-15:00 h, and it converged at 17:00-18:00 h. Hence, the difference between the two genotypes under the DS treatment appeared signif icantly earlier and lasted longer than under the WW treatment. For 297-28, its LWPs over the entire day under both DS and WW treatments were lower than those of HD297. Signif icantly different from the WW treatment, the gap between 297-28 and HD297 under the DS treatment was signif icantly widened during most of the day from the 09:00-10:00 h to the 17:00-18:00 h; it reached its maximum at the 14:00-15:00 h, and it still existed at 17:00-18:00 h. These f indings indicated that the LWP of the 297-28 mutant was signif icantly weaker than that of HD297 and weaker than that of the reference variety SN265. SN265 is a typical rice variety that had lower LWP values over the entire day relative to the drought-tolerant reference variety IRAT109.

    (2) Daily change in the Pn. Photosynthesis is undoubtedly one of the most sensitive physiological and ecological traits in the plants' transient response to f ield drought stress, and it is the most important physiological basis for the above agronomic traits, especially in terms of the aboveground biomass and grain yield formation. The changes in Pn(Fig. 4) indicate that a signif icantly different trend occurred in Pnunder the DS treatment compared to that under the WW treatment. Under the DS treatment, the Pnvalues of 297-28 and SN265 reached their maximum at the 07:00-08:00 h, and then they rapidly decreased at 09:00-10:00 h with the increase in the light intensity and temperature; the values showed a slowed-down, decreasing rate from the 09:00-10:00 h to the 14:00-15:00 h , and then during the time from the 14:00-15:00 h to the 17:00-18:00 h, photosynthesis again rapidly became weak and stopped. The Pnvalues of other genotypes were at a very high start point at the 07:00-08:00 h, they were slowly increased (IAPAR9, idr1-1) or slowly decreased (HD297, IRAT109) at the 09:00-10:00 h, and then they fell into a decline or continued to decline during the time from the 09:00-10:00 h to 14:00-15:00 h and were rapidly decreased to 0 during the time from the 14:00-15:00 h to 17:00-18:00 h. Under the WW treatment, the Pnvalues were enhanced with the increase in light intensity, and they reached their peak at the 09:00-10:00 h; then after 10:00 h, with the increased temperature and the continuous increase in the light intensity, the Pnvalues decreased signif icantly; at the 17:00-18:00 h, photosynthesis stopped as the light intensity became weak.

    For idr1-1, under the DS treatment, the Pns over the entire day were signif icantly higher than those of IAPAR9 and reached a maximum at the 09:00-10:00 h, the optimal time for dry matter accumulation. Compared with the results under the DS treatment, differences in the Pns between idr1-1 and IAPAR9 under the WW treatment were signif icantly smaller, with crossed overlapping between the two genotypes at different time slots during the day. For 297-28, under the DS treatment, the Pnwas signif icantly lower than that of HD297 even as early as 07:00-08:00 h; with the increased light intensity and temperature, the difference between the two genotypes reached their maximum at the 09:00-10:00 h; later, with the decline in photosynthesis in both genotypes, the difference became small. The Pns over the whole day for 297-28 were lower than those of HD297. Under the WW treatment, the Pncurves for 297-28 and HD297 were close and parallel in the morning, and they exhibited a slight cross in the afternoon, indicating no signif icant difference between them over the entire day.

    Analysis of leaf water potential and leaf photosynthesis parameters at key time points(1) Analysis of leaf water potential (LWP, 14:00-16:00 h). From the daily changes in the LWP (MPa, Fig. 3), the LWPs of the rice plants dropped to the lowest point at approximately 15:00 h (i.e., 14:00-16:00 h) and remained relatively stable. At this time slot, the difference in the LWP among the different varieties reached the maximum, and the LWP values measured at this time slot were suitable for comparisons among the species.

    Fig. 4 Daily change in the net photosynthesis rate (P n). The measurement day for the well-watered (WW) treatment was 08-22-2013, and for the drought stress (DS) treatment, it was 08-23-2013. The thick solid line represents a genotype under the DS treatment, whereas the thin dashed line represents a genotype under the WW treatment.

    As shown in Table 4, except for HD297 and IRAT109 in 2013, all the genotypes had lower LWPs under the DS treatment than under the WW treatment. Specif ically, the reduction rates for 297-28 and SN265 reached an extremely signif icant level, with the 2-year average reduction rates being 28.32 and 25.63%, respectively; by contrast, the idr1-1, IAPAR9, and IRAT109 did not show a signif icant decrease. Additionally, under the DS treatment for the 2-year study period, the LWP values for 297-28 (-1.51 and -1.78) and for SN265 (-1.42 and -1.61) both showed the lowest values, which were signif icantly different from those for the other genotypes. By contrast, idr1-1 showed the highest LWP values (-1.15 and -1.39); its LWP approached that of IAPAR9, IRAT109, and HD297 in 2013 and was signif icantly higher than the results of IAPAR9, IRAT109, and HD297 in 2012. Under the WW treatment for the two study years, the LWP values for 297-28 (-1.21 and -1.35) and SN265 (-1.15 and -1.26) were slightly higher or closer to most upland genotypes (except idr1-1). These results indicate that, under the DS treatment, the responses of the paddy rice varieties were sensitive, showing low LWPs, and the reduction rates were large compared to the WW treatment. By contrast, the LWPs of the upland rice varieties were signif icantly higher than those of the paddy rice varieties; their reduction rates were small and relatively stable, with sometimes even higher levels under the DS treatment than under the WW treatment.

    Table 4 also shows that because of the higher degree of drought stress in 2013 than in 2012, the LWPs for all the genotypes were decreased in 2013 compared to 2012, regardless of the treatment. Specif ically, the reduction rates between the 2 years for 297-28 and SN265 under the DS treatment (17.9 and 13.4%) were signif icantly higher than those under the WW treatment (11.6 and 9.6%). HD297 and IRAT109 showed the opposite trend, and their reduction rates (6.4 and 6.9%) were signif icantly lower than those under the WW treatment (25.4 and 20.2%). The reduction rates for idr1-1 and IAPAR9 under the DS treatment (20.9 and 11.0%) were comparable to those under the WW treatment (23.1 and 8.9%). These f indings indicated that the reduction rates in the LWPs for the drought-sensitive genotypes increased as the degree of drought stress increased over the 2 years, whereas those for the drought-tolerant genotypes showed the opposite trend.

    In summary, 297-28 was the most sensitive to drought, and its sensitivity was higher than that of HD297 and even that of SN265, which was consistent with how its agronomic traits responded to drought. The drought tolerance of idr1-1 was the strongest and was better than that of IAPAR9, and it could adapt to both drought-stressed and non-stressed environments.

    Table 4 Performances of the leaf photosynthesis parameters (09:00-11:00 h) and leaf water potential (14:00-16:00 h)1)

    (2) Analysis of the Pnand transpiration rate (09:00-11:00 h). The parameters Pn(μmol CO2m-2s-1) and the transpiration rate (mmol H2O m-2s-1) are the two most important parameters for use in assessing the strength of photosynthesis. Based on the results of the daily change in Pn(Fig. 4), the strongest photosynthesis of rice leaves under the non-DS treatments occurred at the 09:00-11:00 h, a time during which the differences among the different varieties reached the maximum. The photosynthetic parameters measured during this time slot were suitable for comparisons among varieties.

    The Pns were generally decreased for all the genotypes under the DS treatment (Table 4), but the responses were signif icantly different among these genotypes. The droughtsensitive types 297-28 and SN265 had the lowest Pnvalues among all the genotypes under both severe and mild drought (6.27 and 6.60, respectively, under severe drought; 13.03 and 14.17, respectively, under mild drought), and they showed the highest 2-year average reduction rates (33.18 and 41.13%, respectively). Additionally, HD297, IRAT109, and idr1-1 had average reduction rates of 0.15, 4.42, and 0.92%, respectively; during the mild drought year (2012), HD297 and idr1-1 had slightly higher Pns under DS than under WW. The Pns of all the genotypes under the DS and WW treatments decreased in 2013 compared with 2012, exhibiting the same trend as that of the LWP. Specif ically, idr1-1 and IRAT109 decreased slightly, and HD297, IAPAR9, and 297-28 decreased more obviously, whereas SN265 showed almost no decrease under WW, yet it experienced a sharp decrease under DS.

    During both years, 297-28 had very similar Pnvalues with no signif icant difference relative to HD297 under the WW treatment, but these values were 5.1 lower than those of HD297 on average under the DS treatment, showing a signif icant difference. The 2-year average of the reduction rate in the Pnvalues (from under WW to under DS) for 297-28 was 221 times that for HD297, indicating that 297-28 was highly sensitive to drought stress. idr1-1 showed similar Pnvalues under the DS treatment relative to those under the WW treatment; that is, under the mild DS treatment, it had a slightly higher Pnvalue than it did under the WW treatment, and under the severe DS treatment, its Pnvalue was only slightly decreased (by 3.2% compared with that under the WW treatment and showing no signif icant difference), and it was the highest (16.88) among all the genotypes. Compared with the IAPAR9, idr1-1 had slightly decreased Pnvalues under the WW treatment but showed a signif icant increase by 5.22 under the severe DS treatment, leading to the 2-year average reduction rate of idr1-1, which was 26.93% less than that of IAPAR9. The results showed that idr1-1 was very insensitive to drought stress and was not only superior to IAPAR9 but even better than IRAT109, which was the opposite of the 297-28 outcome.

    The transpiration rate(shown in Table 4) showed basically the same trend as that of Pnbetween the years and between the genotypes. Notably, the transpiration rate values for idr1-1 during the 2 DS treatment years were not lower, but instead, they were slightly higher than under those under the WW treatment, and the values for idr1-1 under the DS treatment during the 2 years were both signif icantly higher than those for IAPAR9 and were the highest (6.04 and 7.60) among all the genotypes.

    From the above results, the Pnand transpiration rate (09:00-11:00 h) had the same performance as the LWP (14:00-16:00 h) for each genotype.

    4. Discussion

    4.1. Drought resistance ability of the two mutants

    Two years of combined experimental data showed that under drought stress, the idr1-1 mutant had a smaller reduction rate for most agronomic traits (Table 3, e.g., plant height, setting percentage, and aboveground biomass) and physiological traits (Table 4, e.g., leaf water potential, net photosynthetic rate, and transpiration rate) and had a smaller degree of delay in days to heading (Table 3) than its wildtype IAPAR9, suggesting that it has stable performance in terms of morphological development under drought stress, and idr1-1 has a stronger drought tolerance ability than the wild type. In addition, of all the tested materials, idr1-1 has the lowest reduction in aboveground biomass and the highest LADL value in terms of morphology (Table 3), a very important and effective drought-related parameter identif ied by Zu et al. (2017). In terms of physiology, idr1-1 showed the highest leaf water potential, net photosynthetic rate, and transpiration rate values (Table 4) under drought stress conditions. These results strongly supported the interpretation that idr1-1 had the strongest drought tolerance of the six tested genotypes, which was consistent with the former study (Zu et al. 2017). Considering the study of Zu et al. (2017) on 13 representative upland rice cultivars, idr1-1 might be a superior and strongly drought-tolerant mutant in upland rice, even for the entire rice cultivar germplasm of China at present. By contrast, under drought stress, the 297-28 mutant showed a greater reduction rate in each of the agronomic traits than HD297, and its physiological capacities, such as photosynthesis and water potential, were signif icantly lower than that of HD297. The drought sensitivity of 297-28 was even higher than that of SN265, a common rice variety, indicating that 297-28 is a highly drought-sensitive mutant.

    4.2. Signif icance in theoretical studies

    Both idr1-1 and 297-28 were viewed as two extreme mutants that exhibit strong drought resistance and high drought sensitivity, respectively, and thus they are ideal materials for studying the mechanism of drought resistance. The idr1-1 showed decreased plant height and reduced aboveground vegetative tissues (Table 3); additionally, our unpublished data indicated that idr1-1 had developed root systems with no signif icant difference in the root length (depth), root quantity, root thickness, and root weight and had a shorter leaf length compared with its wild-type IAPAR9. By contrast, the 297-28 showed exactly the opposite changes: an increased plant height, increased vegetative tissues (Table 3), and an evolution from an upland rice root system into a paddy rice root system (e.g., the root length, the number of thick roots and the root thickness were decreased), and the leaf length was increased, and the leaf width was decreased (data not shown). A study by Lilley and Fukai (1994) showed a correlation between the plant height and the comprehensive score for drought resistance. Is the decreased plant height and developed root system of idr1-1 the cause of its increase in drought resistance, and are the increased plant height and undeveloped roots of 297-28 the cause of its decrease in drought resistance? Wider rice leaves would confer better drought resistance (Farooq et al. 2010), which was fully evident in the wide and thick leaves of idr1-1. The mesophyll cells in idr1-1 have higher turgor pressure to keep the stomata open and to maintain a high transpiration rate level and photosynthesis intensity; thus, the leaves would not curl, even under drought stress. By contrast, the 297-28 leaves are thin, narrow, and long and they lose water more readily to form curling shapes, illustrating the opposite effect. These changes in the two mutants suggest that the root system, plant height, and leaf morphology are involved in the drought resistance mechanism.

    4.3. Signif icance in breeding applications

    Mutants serving as the standard control varieties in drought-resistance identif ication experimentsThe two mutants in this study can be used as the standard control varieties in drought-resistance identif ication experiments. At present, the two mutants have been used as control indicators for the two extremes in drought in a standard gradient drought system in the project of “Precise identif ication of the core germplasm resources in rice: precise identif ication in drought” supported by the National Key Research and Development Program of China (2016YFD0100101). The drought-sensitive mutants 297-28 can accurately indicate the onset of a drought period, whereas the strong drought-resistant mutant idr1-1 can indicate whether the drought intensity reaches the desired level. Using these mutants as representatives of each extreme, the accuracy of the drought-resistance identif ication experiments can be controlled under different test environments.

    Mutants serving as the donor parents carrying droughtresistant and dwarf ing genesA previous study indicated that plant dwarf ing and lodging resistance could help increase the yield (Sasaki et al. 2002). Kumar et al. (2008) noted that when using the hybridization of drought-resistant parents and high-yield yet drought-tolerant parents, there is a better chance for obtaining high-yielding and droughttolerant recombinants from the hybridization offspring. In the past, the hybridization of paddy rice and upland rice has been mostly used in conventional breeding for droughttolerant upland rice. Therefore, as a semi-dwarf, idr1-1 is strongly drought-resistant, and as a high-combining-ability material, it would have inherent advantages in conventional breeding. After hybridization with rice, the chance to screen for high-yield and drought-resistant conventional droughttolerant varieties would be high. Our lab has used the idr1-1 mutant in the heterosis utilization of rice by taking it as a bridge for the transfer of the drought-resistant cytoplasmic male sterility (CMS) line. Unpublished data show that idr1-1 can quickly combine its drought-tolerance performance with heterosis, and it is expected to breed high-yield, droughttolerant, and strongly adaptable hybrid rice combinations that can solve food security issues.

    5. Conclusion

    In this study, comprehensive and comparative observations of physiological traits were performed, including morphological development, agronomic traits, water use, and photosynthesis. The results showed that idr1-1 is a dwarf mutant of the high-stalk upland rice variety IAPAR9, and it confers signif icantly enhanced drought resistance; 297-28 is a mutant of the superior drought-resistant variety HD297 and is highly sensitive to drought, with some drought resistance traits that are even weaker than those of the typical rice SN265. These two mutants could be used as important materials in drought resistance identif ication studies on rice cultivar germplasm resources and in drought resistance mechanism studies in the future.

    Acknowledgements

    This work was supported by the National Key Research and Development Program of China (2016YFD0100101), the 948 Program of China (2006-G51), and the European Commission within the 6th Framework Program (ECFP6) INCO-2003-B1.2 (CEDROME-015468). We thank American Journal Experts (AJE) for their English language editing.

    少妇被粗大猛烈的视频| 性插视频无遮挡在线免费观看| 久久午夜亚洲精品久久| 久久精品国产亚洲av天美| 国产激情偷乱视频一区二区| 成人国产综合亚洲| 久久精品国产亚洲av涩爱 | 日本五十路高清| 国产免费av片在线观看野外av| 又黄又爽又刺激的免费视频.| 老司机午夜福利在线观看视频| 1024手机看黄色片| 精品免费久久久久久久清纯| 亚洲av免费高清在线观看| 级片在线观看| 亚洲中文字幕日韩| 亚洲国产欧洲综合997久久,| 日韩高清综合在线| 国产又黄又爽又无遮挡在线| 在线观看66精品国产| 2021天堂中文幕一二区在线观| 国产真实伦视频高清在线观看 | 国产精品av视频在线免费观看| 深爱激情五月婷婷| 久久久国产成人精品二区| 亚洲午夜理论影院| 色在线成人网| 亚洲av熟女| 精品熟女少妇八av免费久了| 欧美3d第一页| 亚洲精品久久国产高清桃花| 国产主播在线观看一区二区| 国产一区二区在线av高清观看| 国产亚洲av嫩草精品影院| 高清在线国产一区| 国产午夜福利久久久久久| 国产精品一及| 国产成人影院久久av| 婷婷丁香在线五月| av在线老鸭窝| 精品久久久久久久人妻蜜臀av| 午夜免费男女啪啪视频观看 | 少妇熟女aⅴ在线视频| 2021天堂中文幕一二区在线观| 我要看日韩黄色一级片| 精品日产1卡2卡| 欧美绝顶高潮抽搐喷水| 国产免费av片在线观看野外av| 乱人视频在线观看| 免费观看人在逋| 国产精品女同一区二区软件 | 国产激情偷乱视频一区二区| 亚洲精品一区av在线观看| 亚洲av电影在线进入| 热99re8久久精品国产| 午夜日韩欧美国产| 搡老岳熟女国产| 国产一区二区亚洲精品在线观看| 综合色av麻豆| 天堂√8在线中文| 俺也久久电影网| 国产综合懂色| 欧美一区二区精品小视频在线| 深夜精品福利| 国产成人av教育| a级一级毛片免费在线观看| 老熟妇仑乱视频hdxx| 精品久久久久久久久av| 国内毛片毛片毛片毛片毛片| 国产精品爽爽va在线观看网站| 99久久无色码亚洲精品果冻| 久久精品国产亚洲av天美| av天堂中文字幕网| 久久99热这里只有精品18| 欧美xxxx黑人xx丫x性爽| 直男gayav资源| 欧美区成人在线视频| 高清日韩中文字幕在线| 51午夜福利影视在线观看| 亚洲精华国产精华精| 亚洲成av人片免费观看| 成人精品一区二区免费| 精华霜和精华液先用哪个| 制服丝袜大香蕉在线| 性欧美人与动物交配| 国产探花极品一区二区| 久久精品夜夜夜夜夜久久蜜豆| 毛片一级片免费看久久久久 | 免费在线观看成人毛片| 国产视频内射| 美女被艹到高潮喷水动态| 日韩免费av在线播放| 成年女人永久免费观看视频| 亚洲精品456在线播放app | 亚洲片人在线观看| 亚洲欧美日韩高清专用| 亚洲一区二区三区不卡视频| 色噜噜av男人的天堂激情| 宅男免费午夜| 久久久成人免费电影| 国产精品日韩av在线免费观看| 狂野欧美白嫩少妇大欣赏| 18+在线观看网站| 2021天堂中文幕一二区在线观| 国产黄色小视频在线观看| 婷婷丁香在线五月| av在线蜜桃| 女生性感内裤真人,穿戴方法视频| 免费观看的影片在线观看| 国产精品嫩草影院av在线观看 | 天天躁日日操中文字幕| 啦啦啦观看免费观看视频高清| 国产高清视频在线播放一区| 综合色av麻豆| 亚洲精品乱码久久久v下载方式| 一级黄色大片毛片| 90打野战视频偷拍视频| 丰满人妻熟妇乱又伦精品不卡| 特级一级黄色大片| 日本成人三级电影网站| 麻豆久久精品国产亚洲av| 成人av一区二区三区在线看| 免费大片18禁| 观看免费一级毛片| 91久久精品电影网| 亚洲七黄色美女视频| 蜜桃久久精品国产亚洲av| 在线a可以看的网站| 欧美一级a爱片免费观看看| 精品久久久久久成人av| 久久久久精品国产欧美久久久| 久久精品国产亚洲av涩爱 | 看免费av毛片| 波多野结衣高清无吗| 亚洲精品色激情综合| 午夜两性在线视频| 亚洲中文日韩欧美视频| 亚洲人成网站高清观看| 91麻豆精品激情在线观看国产| 97人妻精品一区二区三区麻豆| 麻豆国产97在线/欧美| 97热精品久久久久久| 国产精品亚洲av一区麻豆| 国产色爽女视频免费观看| 国产精品不卡视频一区二区 | 国产一区二区在线观看日韩| 男女床上黄色一级片免费看| 极品教师在线免费播放| 午夜精品久久久久久毛片777| 老司机午夜福利在线观看视频| 国产精品99久久久久久久久| 欧美bdsm另类| 天天躁日日操中文字幕| 亚洲美女搞黄在线观看 | 国产美女午夜福利| 又黄又爽又刺激的免费视频.| 国产午夜福利久久久久久| 男人舔女人下体高潮全视频| а√天堂www在线а√下载| 午夜激情福利司机影院| 国产在线男女| 国产一级毛片七仙女欲春2| 国产精品av视频在线免费观看| 老鸭窝网址在线观看| 首页视频小说图片口味搜索| 别揉我奶头 嗯啊视频| 91狼人影院| 三级国产精品欧美在线观看| 国产精品一区二区性色av| 精华霜和精华液先用哪个| 久久久精品欧美日韩精品| aaaaa片日本免费| 亚洲国产精品合色在线| 国产亚洲欧美在线一区二区| 高清在线国产一区| 亚洲美女黄片视频| 午夜精品久久久久久毛片777| 久久精品国产亚洲av涩爱 | 可以在线观看毛片的网站| 日韩欧美精品v在线| 午夜老司机福利剧场| 麻豆一二三区av精品| 免费看a级黄色片| 91麻豆精品激情在线观看国产| 亚洲中文字幕一区二区三区有码在线看| 91av网一区二区| 久久久久久国产a免费观看| 久久久久国内视频| x7x7x7水蜜桃| 欧美最新免费一区二区三区 | 老司机午夜十八禁免费视频| 欧美最新免费一区二区三区 | 国产精品一区二区免费欧美| 男人舔女人下体高潮全视频| 在线看三级毛片| 亚洲av二区三区四区| 国产精品日韩av在线免费观看| 午夜福利成人在线免费观看| 欧美一区二区精品小视频在线| 岛国在线免费视频观看| 午夜免费男女啪啪视频观看 | 一级a爱片免费观看的视频| 国产精品久久视频播放| 亚洲第一电影网av| 日韩精品中文字幕看吧| 亚洲av成人不卡在线观看播放网| 欧美午夜高清在线| 麻豆成人午夜福利视频| 一个人看的www免费观看视频| 国产黄a三级三级三级人| 91麻豆av在线| 国产精品影院久久| 男人的好看免费观看在线视频| 国产亚洲精品av在线| 我的女老师完整版在线观看| 亚洲欧美精品综合久久99| av在线蜜桃| 一区二区三区四区激情视频 | 91在线精品国自产拍蜜月| 久久久精品大字幕| 国产精品久久久久久精品电影| x7x7x7水蜜桃| 每晚都被弄得嗷嗷叫到高潮| 欧美成人一区二区免费高清观看| 国产精品亚洲美女久久久| 亚州av有码| 国产单亲对白刺激| 久久精品人妻少妇| 一本久久中文字幕| 国产精品98久久久久久宅男小说| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩高清在线视频| 久久久成人免费电影| 欧美色欧美亚洲另类二区| 熟女电影av网| 男女床上黄色一级片免费看| av在线观看视频网站免费| 久久欧美精品欧美久久欧美| 亚洲内射少妇av| 一个人观看的视频www高清免费观看| 99热6这里只有精品| 99热精品在线国产| 国产精品综合久久久久久久免费| 精品午夜福利在线看| 最近在线观看免费完整版| 真人做人爱边吃奶动态| 久久欧美精品欧美久久欧美| 成人午夜高清在线视频| 久久人妻av系列| 大型黄色视频在线免费观看| 免费大片18禁| 亚洲av电影不卡..在线观看| 午夜两性在线视频| 午夜福利在线在线| 噜噜噜噜噜久久久久久91| x7x7x7水蜜桃| 91狼人影院| 亚洲精品日韩av片在线观看| 一夜夜www| 极品教师在线视频| 中文亚洲av片在线观看爽| 亚洲av.av天堂| 久9热在线精品视频| 91狼人影院| 国产精品,欧美在线| 热99re8久久精品国产| 婷婷精品国产亚洲av| 国产高清视频在线观看网站| 国产亚洲欧美98| 日本成人三级电影网站| 国产精品日韩av在线免费观看| 婷婷丁香在线五月| 欧美三级亚洲精品| 亚洲狠狠婷婷综合久久图片| 亚洲男人的天堂狠狠| 99久久久亚洲精品蜜臀av| 一个人免费在线观看电影| 精品人妻一区二区三区麻豆 | 国产av不卡久久| 99热6这里只有精品| 中文亚洲av片在线观看爽| 久久久精品欧美日韩精品| 欧洲精品卡2卡3卡4卡5卡区| 熟女电影av网| 99国产综合亚洲精品| 久久久久精品国产欧美久久久| 国产极品精品免费视频能看的| 亚洲无线在线观看| 日本黄色片子视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲,欧美精品.| 很黄的视频免费| 精品福利观看| 久久久久久久久久黄片| 色视频www国产| 色5月婷婷丁香| 亚洲精品成人久久久久久| 国产av不卡久久| 欧美日韩福利视频一区二区| 少妇的逼水好多| 无遮挡黄片免费观看| 黄色视频,在线免费观看| 如何舔出高潮| 国产精品亚洲美女久久久| 午夜福利在线在线| 精品免费久久久久久久清纯| 天堂影院成人在线观看| 黄色女人牲交| 老司机福利观看| 午夜影院日韩av| 国产精品美女特级片免费视频播放器| 99久久久亚洲精品蜜臀av| 男人的好看免费观看在线视频| h日本视频在线播放| 亚洲av第一区精品v没综合| 日韩欧美在线乱码| 我要看日韩黄色一级片| 久久久久免费精品人妻一区二区| 深爱激情五月婷婷| 欧美+日韩+精品| 我要看日韩黄色一级片| 亚洲无线观看免费| 少妇裸体淫交视频免费看高清| 少妇的逼水好多| 精品久久久久久久久av| 天堂av国产一区二区熟女人妻| 欧美黄色片欧美黄色片| 天堂网av新在线| 亚洲av免费在线观看| 又爽又黄a免费视频| 久久精品久久久久久噜噜老黄 | 91九色精品人成在线观看| 亚洲无线在线观看| netflix在线观看网站| 88av欧美| 免费在线观看日本一区| 成年女人永久免费观看视频| 亚洲中文日韩欧美视频| 国产精品自产拍在线观看55亚洲| 色5月婷婷丁香| 久久久精品欧美日韩精品| 中文字幕高清在线视频| 久久久久久大精品| 午夜免费男女啪啪视频观看 | 人人妻人人看人人澡| 久久婷婷人人爽人人干人人爱| bbb黄色大片| 国产亚洲精品综合一区在线观看| 又黄又爽又刺激的免费视频.| 99久久精品一区二区三区| 搡老岳熟女国产| 日韩有码中文字幕| 又黄又爽又免费观看的视频| 他把我摸到了高潮在线观看| 亚洲欧美日韩高清专用| 亚洲av五月六月丁香网| 在线观看美女被高潮喷水网站 | 精品不卡国产一区二区三区| 我要搜黄色片| 中文字幕高清在线视频| 99热只有精品国产| 免费黄网站久久成人精品 | 亚洲精品456在线播放app | 久久99热这里只有精品18| 首页视频小说图片口味搜索| 国产精品一区二区性色av| 欧美午夜高清在线| 欧美3d第一页| 亚洲第一区二区三区不卡| 女生性感内裤真人,穿戴方法视频| 婷婷丁香在线五月| 久久欧美精品欧美久久欧美| 精品欧美国产一区二区三| 三级男女做爰猛烈吃奶摸视频| 亚洲av一区综合| a级一级毛片免费在线观看| 高清日韩中文字幕在线| 如何舔出高潮| 亚洲专区中文字幕在线| 免费观看的影片在线观看| 国产毛片a区久久久久| 亚洲av电影在线进入| 国产成人a区在线观看| 亚洲国产精品sss在线观看| 非洲黑人性xxxx精品又粗又长| 一级黄色大片毛片| 免费看美女性在线毛片视频| 久久国产精品人妻蜜桃| 村上凉子中文字幕在线| 午夜免费激情av| 特级一级黄色大片| 成年女人永久免费观看视频| 99久久成人亚洲精品观看| 欧美一区二区精品小视频在线| 两个人视频免费观看高清| 国产激情偷乱视频一区二区| 久久精品国产亚洲av香蕉五月| 一级黄片播放器| 午夜福利在线观看免费完整高清在 | 一区二区三区高清视频在线| 在线观看av片永久免费下载| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品合色在线| 嫩草影院精品99| 国内精品久久久久精免费| 国产精品一及| 别揉我奶头 嗯啊视频| 午夜a级毛片| 精品日产1卡2卡| 国产真实伦视频高清在线观看 | 色av中文字幕| 婷婷色综合大香蕉| 国产精品久久久久久久久免 | 观看美女的网站| 欧美日本视频| 听说在线观看完整版免费高清| 亚洲av中文字字幕乱码综合| 精品乱码久久久久久99久播| 色尼玛亚洲综合影院| 久久这里只有精品中国| 久久草成人影院| 最好的美女福利视频网| 在线看三级毛片| 欧美+亚洲+日韩+国产| 亚洲国产色片| 精品国产亚洲在线| 最后的刺客免费高清国语| 一级黄色大片毛片| 欧美黑人巨大hd| 成年版毛片免费区| 一级毛片久久久久久久久女| 色综合婷婷激情| 精品久久国产蜜桃| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久大av| 又爽又黄a免费视频| 草草在线视频免费看| 精品99又大又爽又粗少妇毛片 | 香蕉av资源在线| 极品教师在线视频| 日韩欧美国产一区二区入口| 在现免费观看毛片| av福利片在线观看| 丝袜美腿在线中文| 91九色精品人成在线观看| 午夜福利免费观看在线| 日韩大尺度精品在线看网址| 亚洲,欧美精品.| 久久久久精品国产欧美久久久| 国产伦精品一区二区三区视频9| 桃色一区二区三区在线观看| 波多野结衣高清无吗| 无人区码免费观看不卡| 国产精华一区二区三区| 91在线观看av| 长腿黑丝高跟| 天堂影院成人在线观看| 国产蜜桃级精品一区二区三区| 成熟少妇高潮喷水视频| 18禁黄网站禁片免费观看直播| 美女高潮的动态| av在线观看视频网站免费| 国产中年淑女户外野战色| 好男人电影高清在线观看| 99在线视频只有这里精品首页| 悠悠久久av| 日韩亚洲欧美综合| 我的老师免费观看完整版| av欧美777| 中文字幕精品亚洲无线码一区| 嫩草影视91久久| 国产欧美日韩一区二区精品| 国产精品一区二区性色av| 亚洲成人久久爱视频| 91狼人影院| 国产伦人伦偷精品视频| 久久久久久久久中文| 噜噜噜噜噜久久久久久91| 欧美+日韩+精品| 最近中文字幕高清免费大全6 | 热99re8久久精品国产| 成人av一区二区三区在线看| 好男人在线观看高清免费视频| 男人舔奶头视频| 婷婷丁香在线五月| 亚洲国产欧洲综合997久久,| 欧美+日韩+精品| 日韩高清综合在线| 亚洲精品在线观看二区| www.熟女人妻精品国产| 久久久久久久午夜电影| 精品福利观看| 国产美女午夜福利| 精品国产三级普通话版| 国产伦精品一区二区三区四那| 午夜亚洲福利在线播放| 亚洲精品在线美女| 色视频www国产| 欧美午夜高清在线| 亚洲午夜理论影院| 一区二区三区免费毛片| 97碰自拍视频| 欧美高清性xxxxhd video| 五月玫瑰六月丁香| 国产探花极品一区二区| 男人和女人高潮做爰伦理| 国产av在哪里看| 欧美国产日韩亚洲一区| 欧美性猛交黑人性爽| 国产黄片美女视频| 三级毛片av免费| 一区福利在线观看| 天堂av国产一区二区熟女人妻| 欧美激情在线99| 国产精品一区二区性色av| 免费看日本二区| 亚洲自偷自拍三级| 久久久久久久久久成人| 久久九九热精品免费| 久久国产乱子伦精品免费另类| 国产 一区 欧美 日韩| 欧美性猛交╳xxx乱大交人| 欧美区成人在线视频| 国产中年淑女户外野战色| 国产免费av片在线观看野外av| 午夜视频国产福利| 久久天躁狠狠躁夜夜2o2o| 日本免费一区二区三区高清不卡| 亚洲男人的天堂狠狠| 国产一区二区激情短视频| 亚洲色图av天堂| 欧美xxxx性猛交bbbb| 草草在线视频免费看| 成人精品一区二区免费| 中国美女看黄片| av天堂在线播放| 天堂影院成人在线观看| 国产精品女同一区二区软件 | 十八禁网站免费在线| 天天一区二区日本电影三级| 真人一进一出gif抽搐免费| 亚洲国产精品999在线| 国内精品久久久久精免费| av视频在线观看入口| 欧美一区二区精品小视频在线| 精品99又大又爽又粗少妇毛片 | aaaaa片日本免费| 国产精品一及| 一个人看视频在线观看www免费| 高清毛片免费观看视频网站| 精品久久久久久,| 亚洲欧美日韩高清在线视频| 女人十人毛片免费观看3o分钟| 12—13女人毛片做爰片一| 99热这里只有精品一区| 国产精品爽爽va在线观看网站| 国产精品亚洲一级av第二区| 国产亚洲欧美98| 欧美黑人欧美精品刺激| 国产精品乱码一区二三区的特点| 成人av在线播放网站| 亚洲美女视频黄频| 国产精品人妻久久久久久| 亚洲欧美清纯卡通| 丁香欧美五月| 99久久精品一区二区三区| 国产精品久久视频播放| 欧美日本亚洲视频在线播放| 高潮久久久久久久久久久不卡| 69av精品久久久久久| 日韩精品青青久久久久久| 乱码一卡2卡4卡精品| 午夜亚洲福利在线播放| 麻豆成人av在线观看| 非洲黑人性xxxx精品又粗又长| 国产精品,欧美在线| 日本熟妇午夜| 狂野欧美白嫩少妇大欣赏| 国产精品一及| 在线十欧美十亚洲十日本专区| 午夜精品久久久久久毛片777| 久久婷婷人人爽人人干人人爱| 亚洲国产日韩欧美精品在线观看| 午夜福利视频1000在线观看| 美女免费视频网站| 精品久久久久久久人妻蜜臀av| 性色avwww在线观看| 97人妻精品一区二区三区麻豆| 日日摸夜夜添夜夜添小说| 午夜福利在线观看吧| 日韩欧美精品v在线| av专区在线播放| 99riav亚洲国产免费| 久久国产乱子免费精品| 看片在线看免费视频| 高清毛片免费观看视频网站| 欧美在线黄色| 99久久99久久久精品蜜桃| 久久午夜亚洲精品久久| 国产主播在线观看一区二区| 亚洲无线观看免费| 午夜精品久久久久久毛片777| www.色视频.com| 国产精品亚洲美女久久久| 在线观看午夜福利视频| 村上凉子中文字幕在线| 国产欧美日韩精品一区二区| 99riav亚洲国产免费| 3wmmmm亚洲av在线观看| 亚洲午夜理论影院| 国产欧美日韩一区二区精品| www.色视频.com| 少妇熟女aⅴ在线视频| 很黄的视频免费| 久久精品影院6|