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

    Morphological and Molecular Characterization of Novel Salt-tolerant Rice Germplasms from the Philippines and Bangladesh

    2019-05-23 01:46:00WilsonAalaJrGlennGregorio
    Rice Science 2019年3期

    Wilson F. Aala Jr, Glenn B. Gregorio

    ?

    Morphological and Molecular Characterization of Novel Salt-tolerant Rice Germplasms from the Philippines and Bangladesh

    Wilson F. Aala Jr, Glenn B. Gregorio

    ()

    To screen for new sources of salinity tolerance, 688 traditional rice varieties from the Philippines and Bangladesh were obtained, and their tolerance to hypersaline conditions at the seedling stage was examined. A total of 29 Philippine lines and 15 Bangladeshi lines were scored as salt-tolerant. Morphological assessment (plant height, biomass and Na-K ratio) revealed that among the 44 salt-tolerant accessions, Casibon, Kalagnon and Ikogan had significantly higher relative shoot length difference, relative shoot growth reduction and shoot Na-K ratio than the tolerant check FL478. Additionally, AC and Akundo exhibited significantly higher Na-K ratios than the other genotypes. The genetic diversity of the 44 genotypes was assessed using 34 simple sequence repeat markers. A total of 133 alleles were detected across all loci. Cluster analysis showed that AC, Akundo and Kuplod were clustered along with FL478, indicating a strong genetic relatedness between these genotypes. IR29 (susceptible check) was singly separated. The haplotype analysis revealed that none of the 44 genotypes had a similar allele combination as FL478. These accessions are of interest since each genotype might be different from the classical salinity-tolerant Pokkali.

    rice;simple sequence repeat; Na-K ratio; salinity tolerance; traditional landrace

    Rice feeds about half of the world’s population, and constitutes 30% to 80% of the daily caloric requirements in Asia (Lafitte et al, 2004). Rice is generally grown on coastal areas flanked by sea water, and as thus salinity is one of the major abiotic stresses affecting its growth and yield (Mohamed et al, 2007). In Bangladesh, salinity affects 30% of its total agricultural lands, with over 1.02 × 106hm2of coastal land constantly under the direct effects of salinity (Mazid et al, 2016). The Philippines, which also harbours coastal agricultural areas, have reported a total of 0.4 × 106hm2of land being affected by salinity, with half categorized as severe salinity (PhilRice, 2011).

    Among the cereals, rice ranks as the most saline-sensitive, losing 10% of its yield at 3 dSm-1of salinity, and moderate levels of salinity (7 dSm-1) could potentially reduce its yield by 50% (Umali, 1993). Excess soil salinity induces osmotic stress, ion toxicity and nutrient imbalance (Flowers and Flowers, 2005), leading to a reduction in shoot length, fresh weight and dry weight (Khan et al, 1997; Roy et al, 2002; Janaguiraman et al, 2003). It is most widely accepted that salt injury is due to an imbalance in the shoot Na-K ratio, which reduces tillering capacity and grain yield (Munns and Tester, 2008; Horie et al, 2012).

    Rice is most sensitive to salt stress at the early seedling stage (Gregorio et al, 1997), and hence seedlingstage screening of genotypes for salt tolerance remains valuable due to its rapidity and scalability, which is difficult for the vegetative and reproductive stages (Bhowmik et al, 2009). However, salinity tolerance is heavily affected by the environment, and narrow-sense heritability estimates from such experiments tend to be low (Gregorio and Senadhira, 1993). The efficacy of mass phenetic screening could be increased if aided by suitable genetic techniques (Schmidt et al, 2014). Molecular markers such as simple sequence repeat (SSR), which can tag major genes more commonly known as quantitative trait loci (QTL), holds much value in marker-assisted selection (MAS) (Rines et al, 2006;Ganie et al, 2014). A major QTL on chromosome 1, termed as, has been identified to be tightly associated with salinity tolerance (Mohammadi-Nejad et al, 2009; Islam et al, 2011), and SSR markers linked to this QTL have also been reported in Bonilla et al (2002). Use of such markers to delineate test genotypes and aid in selection have also been reported previously (Moreno et al, 2003; Xu and Crouch, 2008).

    Genetic and phenetic screening relies on the wide array of alleles present in the gene pool representing a multitude of possible traits for salt tolerance (Reddy et al, 2017). Previous studies involving genetic and physiological screening have yielded highly salt-tolerant lines like Pokkali and FL478 (Thomson et al, 2010), although novel sources of salt-tolerant genotypes that could be used as the parents in breeding programs are few (Lisa et al, 2004). In this regard, landraces have been shown to be good sources of new alleles (Hoisington et al, 1999;Guevarra et al, 2001), which remains greatly uncharacterized and unused. The purpose of this study was to identify and characterize new sources of salinity tolerance from traditional landraces in the Philippines and Bangladesh by utilizing morphological and molecular screening techniques (SSR analysis).

    Materials and Methods

    Plant materials

    A total of 688 rice accessions, 505 Philippine traditional landraces [including FL478 (tolerant check) and IR29 (susceptible check)] and 183 Bangladeshilandraces, were collected from the International Rice Research Institute Gene Bank. Seed dormancy was broken by incubation at 50oC for 5 d. The seeds were then sterilized using a 0.2% bleach solution for 5 min, and allowed to pre-germinate in a dark room for 3 d.

    Initial seedling stage screening for salt-tolerant accessions

    Two pre-germinated seeds were sown in each hole of a styrofoam seedling float wherein four holes were allotted per accession. Each float was placed on top of an 8L plastic tray containing tap water and placed inside the IRRI Phytotron (temperatures of 29 oC day /21 oC night, minimum relative humidity of 70%). After 5 d, the water in each tray was replaced with a half-strength Simple Nutrient Addition Program A and B (SNAP) culture solution (0.5% SNAP A, 0.5% SNAP B, 99% H2O, 200 mg FeSO4) salinized to EC=12dSm-1using analytic sodium chloride, and pH was maintained at 5.0. After 14 d, the degree of salt stress injury was evaluated using the modified standard evaluating score (SES) (Table 1). After scoring, the salinity level was increased to EC of 18dSm-1, and scoring was carried out again.

    Table 1. Modified standard evaluation score for assessing the visual symptoms of salt injury at the seedling stage.

    Replicated screening of selected accessions

    The replicated screening process was essentially similar with the initial screening. Three salinized and normal replicates were prepared using the randomized complete block design. The tolerant genotypes from the initial phenotyping, as well as IR29 and FL478, were seeded in seedling floats and placed in trays containing Peter’s culture solution (1 g solid Peter’s fertilizer, 200 mg FeSO4, 1 L H2O, pH=5). The pre-germinated seeds were allowed to grow in the normal solution (culture solution only) for one week with regular checking of pH for acclimatization purposes. Salinization to an EC value of 12dSm-1was done at 7 d after seeding. Scoring was carried out 14 d after salinization, and salinization to an EC of 18dSm-1was done thereafter. Scoring was carried out after two weeks or prior to sample acquisition (Table 2).

    Physio-morphometric screening

    Fresh samples were obtained after scoring of the samples were carried out. The plants were cut at the shoot-root interface, and measured bya standard meter stick with millimeter gradations. For the biomass analysis, shoot fresh and dry weights were obtained. The fresh weight (three plants per accession) was measured using an analytical balance. The samples were then placed inside paper envelopes with proper labels and dried inside an oven for 3 d with an ambient temperature of 65oC. The dry weight of each sample was then measured using an analytical balance.

    Table 2. Physio-morphometric data for the 44 salt-tolerant accessions under salinized conditions (EC = 18 dSm-1).

    CN, Control number; ACCNO, Accession number at IRRI genebank; BD, Bangladesh; PH, the Phillippines; SES, Standard evaluation score; SL, Shoot length; SFW, Shoot fresh weight; SDW, Shoot dry weight; RSLD, Relative shoot length difference; RSWC, Relative shoot water content; RSDR, Relative shoot dry weight reduction.

    Means in the same column with the same lowercase letters are not significantly different at the 0.05 level.

    Na-K ratio wasmeasured after the sensitive check IR29 registered a score of 7. The whole shoot was washed thrice with nanopure water,and dried in an oven for 3 d with an ambient temperature of 65oC. The samples were then finely cut, and 10 mg of each sample was placed in individual vials with proper labels.Sample processing for Na-K ratio determination was based on Yoshida et al (1976) with modifications. For each vial, 10 mL of 100 mmol/L acetic acid was dispensed, and the vials were then placed in a hot water bath set at 90oC for 2 h. The samples were then allowed to cool down for 12 h. Prior to reading, filtering of the shoot particles, if any, was carried out so as not to impede the flow of the sample into the feed tube. Samples were diluted (10× dilution) by adding 9mL of deionized water to 1mL of the sample solution. Once ready, a feed tube was placed inside the vials and a flame photometer was used to read the sodium and potassium content of the samples.

    Analyses of physio-morphometric data

    Mean shoot length (SL), shoot fresh weight (SFW), shoot dry weight (SDW), SES score, relative shoot length difference (RSLD), relative shoot water content (RSWC), relative shoot dry weight reduction (RSDR), and Na-K ratio were analyzed using the analysis of variance (ANOVA) in R software.

    RSLD (%) = (SLN– SLS) / SLN× 100

    RSWC (%) = (SFWS– SDWS) / SFWN× 100

    RSDR (%) =(SDWN– SDWS) / SDWN× 100

    Where N indicates normal treatment and S indicates salinized treatment. The two-tailed-test in R software was utilized to compare means between the normal and salinized conditions. Duncan’s multi-range test in R software was used to separate the means of accessions in variables that showed statistical significance in the ANOVA. Linear correlation analysis (Pearson) was done for SES score and Na-K ratio.

    Molecular screening

    Total DNA was extracted from 21-day-old leaf samples as described in Delaporta et al (1983). DNA quality and quantity were assessed using 1% agarose gel electrophoresis using 50, 100 and 150 ng/μL Lambda DNA as reference.

    Genetic relationship for the 44 selected traditional varieties was determined using 38 previously reported SSR markers which included 16markers (RM315, RM490, RM1287, RM10694, AP3206f, RM3412b, RM10748, RM493, RM10793, RM562, RM7075, RM10745, RM10764, RM10772, RM140 and RM495) on chromosome 1and 2 SSR markers each for chromosomes 2 to 12. Information about all the markers were collected from the Gramene database (http://www.gramene.org).

    Amplification of the targeted loci were done in a 20 μL reaction volume containing 300 ng template DNA, 0.05 mmol/L dNTP mix, 0.25 μmol/L forward and reverse primers, 1× PCR buffer, 0.45 mmol/L MgCl2, 4UDNA polymerase, sterile nanopure water to volume and an overlay of sterile mineral oil. The amplification profile consisted of an initial denaturation at 94oC for 5 min; followed by 35 cycles of 94oC denaturation for 1 min, 55oC to 65oC annealing for 45 s, 72oC extension for 1 min; and a final extension at 72oC for 5 min in a G-Storm Thermal Cycler (Life Science Research Products).

    Amplification products were resolved in 6% or 8% native polyacrylamide gel and documented using the Alpha Imager system (Fisher Scientific) (Supplmental Fig. 1). Band scores were summarized and exported to the PowerMarker v.3.0 software (Liu and Muse, 2005) as a tab-delimited text file. From this, the allele number, expected heterozygosity (genetic diversity), and polymorphism information content (PIC) were obtained. A frequency-based distance algorithm was employed to determine the relative similarity of the accessions to each other, and an unweighted pair-group mean analysis using the similarity index as described in Nei et al (1983) was employed to construct a dendrogram using PYLIP v. 3.965 (Felsentein, 1989).

    Results

    Physio-morphometric screening

    Among the 688 traditional accessions, 44 were found to be moderately to highly salt-tolerantwith EC of 18 dSm-1. ANOVA showed that there were statistical significances within replicates, and between genotypes (<0.05) for SES score (Table 3). ANOVA for RSLD showed statistical significance across genotypesat the 0.01 level whereasRSDR showed statistical significance between genotypes at the 0.05 level. ANOVA for shoot length, RSWC, SDW, SFW and Na-K ratio indicated statistical significances across genotypes and for the combined treatment (<0.01) (Tables 3 and 4). Two-tailed-test between the normal and salinized treatments showed statistical significance for all the test traits (<0.05). Duncan’s multi-range test revealed that none of the test accessions had a mean SES score that was statistically similar to IR29 upon salinization, showing that the initial screening was effective in identifying putative salt tolerant cultivars. Betalga, Mararnao, Ikogan and Anangka had statistically lower mean RSLD than FL478 (Table 2). Thirty accessions and FL478 had mean RSWC which were statistically different from IR29. The mean RSDR of Rao, Ikogan, Casibon and Kalagnonwas statistically lower than that of FL478.The mean Na-K ratio of all the test accessions, excluding IR4493-5-5-3, were statistically lower than that of IR29 under saline conditions. Only AC had a mean Na-K ratio that was significantly lower than FL478.

    Table 3. Mean square values of each physiomorphometric data.

    SES, Standard evaluation score; RSWC, Relative shoot water content; RSDR, Relative shoot dry weight reduction; SDW, Shoot dry weight; SFW, Shoot fresh weight; SL, Shoot length; RSLD, Relative shoot length difference.

    * and ** indicate significant differences at the 0.05 and 0.01 levels, respectively.

    The frequency distribution of the 44 accessions and the two checks showed a bell-shaped curve for RSLD and RSDR which reflected the quantitative nature of these traits (Fig.1-A and B). Using a simple linear correlation analysis, the mean Na-K ratio was found to have a strong positive relationship with the mean SES score (=0.72). Computing for the2value revealed that Na-K ratio accounted for 51.4% of the variation within the 46 samples with respect to the SES scores (Fig.1-C).

    Fig. 1. General distributions of selected accessions for relative shoot length difference (A), relative shoot dry weight reduction (B) and scatter plot of Na-K ratio to standard evaluation score (SES) (C).

    Table 4. Mean square values using the combined treatment as an error term in each of the physio-morphometric data.

    RSWC, Relative shoot water content; SDW, Shoot dry weight; SFW, Shoot fresh weight; SL, Shoot length.

    ** indicates significant difference at the 0.01 level. Combined treatment = Normal vs. Saline.

    Microsattelite marker-assisted screening

    For each marker, the genotype number, major allele frequency, allele number, gene diversity, observed heterozygosity and PIC were obtained for each locus (Table 5). Across the 34 loci considered, a total of 133 alleles were detected ranging from 2 to 10 alleles per locus. Heterozygous genotypes were observed using RM490, RM7075, RM13197, RM152, RM413, RM10764, RM105, RM171 and RM495. The mean gene diversity was observed to be 0.4910, while the mean PIC value was found to be 0.4579. PIC values ranged from 0.0416 to 0.8555 wherein a total of 18 markers showed PIC values higher than 0.5000.

    Utilizing the polymorphisms detected across the 34 loci, a dendogram was constructed and bootstrap trees were generated to form the consensus tree as shown in Fig.2. FL478 was found to form a sub-clade along with AC, Kuplod and Akundo (boxed). IR29 formed a separate clade from the rest of the 45 accessions. Most Bangladeshi lines co-clustered together, and a similar observation could be made for the Philippine lines. A genotype analysis was also carried out (Fig.3) which showed that none of the 44 test genotypes had allele combination similar to FL478 which was used as the reference genotype. AC exhibited the most number of common alleles (21/34), while IR5494 had the lowest number of common alleles (11/34).

    Fig. 2. Consensus tree for 46 rice accessions using 34 simple sequence repeat markers.

    Boxed accessions clustered tightly along with FL478.

    Table 5. Summary statistics for the 34 markers used to screen the selected rice accessions.

    MAF, Major allele frequency; GN, Genotype number; AN, Allele number; GD, Genetic diversity; H, Observed heterozygozity; PIC, Polymorphism information content.

    Fig. 3. Haplotype analysis for the 46 rice accessions using FL478 as the reference genotype.

    Black boxes represent similarity to FL478 allele.

    Discussion

    Mining for genetic and phenotypic variation for salinity tolerance is key to generating better hybrids, and traditional varieties which are mostly uncharacterized holds much potential since each cultivar might possess unique adaptations which are different from Pokkali or FL478, the most commonly used salt-tolerant parents by rice breeders (Mishra et al, 2003; Mohammadi-Nejad et al, 2009; Hoang et al, 2016). In this study, a total of 688 Philippine and Bangladeshi rice accessionswere screened for salt stress response at seedling stage. Based on the SES scores at 18 dSm-1, 44 accessions were identified as salt tolerant. The salinity intensity of 18dSm-1was chosen since rice salt tolerance allows 50% emergence at such levels (Wahhab, 1961), and coastal areas in Bangladesh and the Philippines experience a typical EC of 18dSm-1(Iqbal and Shoajib, 1993). The pairwise mean analysis of the SES scores in the salinized treatment showed that none of the landraces have a mean SES score that was statistically similar tothat of the susceptible check IR29. This shows the discriminatory ability of the standardized evaluation score in evaluating salinity tolerance at seedling stage. Indeed, the efficiency of a standardized score in screening salt-tolerant cultivars had also been shown previously (Lisa et al, 2004;Sexcion et al, 2009).

    Key physiomorphological indicators of salt stress (shoot length, fresh weight, dry weight and Na-K ratio)are correlated with salinity tolerance at seedling stage in rice (Bhowmik et al, 2009). Treatment with NaCl resulted in a decline in the mean shoot length of the genotypes when compared with the control, which is in agreement with previous reports (Javed and Khan, 1975; Saxena and Pandey, 1981; Zeng and Shannon, 2000). However, since plant height is highly cultivar-specific, RSLD was also computed. Pairwise comparison of mean RSLD across genotypes revealed that Betalga, Maranao, Ikogan and Anangka were statistically lower than FL478, wherein Ikogan had a mean RSLD of 0.00. These four landraces might have better Na+sequestering mechanism than FL478 since RSLD could be attributed to leaf senescence. Leaf senescence in the setting of salinization is mainly due to an excess of ions, mostly Na+and Cl-, in the leaf tissue (Munns and Tester, 2008). For most cases, senescence starts at the leaf tips of older leaves and progresses to younger leaves. Closing of stomates due to increased influx of Na+could also increase leaf temperature, leading to reduced shoot length growth (Rajendran et al, 2009; Sirault et al, 2009).

    Shoot fresh weight and dry weight were significantly reduced upon treatment with NaCl. In both parameters, Akundo ranked the highest for SFW (1.708) and SDW (0.413). This was followed by FL478 with SFW of 1.303 and SDW of 0.340 (Table 1). For SFW, AKundo registered a significantly higher value than FL478, whereas their SDW values were not significantly different. IR5494 ranked the lowest (right after IR29) for both SFW and SDW. In wheat, biomass reduction is an initial response to extreme salt stress due to osmotic stress (James et al, 2002). This osmotic imbalance hinders the normal uptake of water and nutrients, leading to reduced tissue development. To detect the difference in performance of the test accessions under normal and salinized conditions for shoot dry weight, RSDR was calculated. The multicomponent test for RSDR showed that FL478, IR29 and 40 of the test accessions had non-significantly different means. This observation appears to contradict the notion that tolerant accessions should have lower biomass reduction than susceptible cultivars. However, cultivars which are sensitive to salt stress (glyophytes) induce early plant vigour in an attempt to avoid salinity-induced damage by rapid growth (Kumar et al, 2013; Reddy et al, 2017). Fast growing rice seedlings were shown to sequester less Na+than slow growers (Yeo et al, 1990). A similar finding is also reported by Asfaw (2011) in sorghum wherein seedling shoot dry weight and fresh weight show no significant difference across genotypes. Lee et al (2003) also reported the same occurrence in Pokkali (tolerant) and IR29, wherein their mean RSWRare non-significantly different under saline conditions at the seedling stage,and IR29 shows significantly lower relative reduction values than Pokkali, whichis attributed to its dwarf stature relative to the tall-type Pokkali. Further scrutiny of the RSDR data showed that Kalagnon, Casibon, Ikogan and Roa have significantly lower mean RSDR than the other accessions. Kalagnon and Casibon showed non-significant difference inshoot lengh comparedto FL478, while Roa and Ikogan hadsignificantly higher values than FL478, which might indicate that RSDR could not be easily correlated with plant height.

    Increased influx of Na+in the roots perturbs the osmotic potential of plant cells and induces physiological drought, which dampens the ability of the plant to acquire water (Verslues et al, 2006). Thus, RSWC was calculated to assess the osmotic stress-induced damage in the test accessions as a function of the relative water content of the plants. Here, Casibon, Ampipit and Binutiti registered significantly superior values when compared to FL478. The ability to tolerate the increased extracellular levels of Na+could be due to: synthesis of osmoprotectants (tertiary amines, dimethylsulfoniopropionate and sulfonium compounds) which stabilizes extracelluar and vacuolar Na+(Turkana and Demiral, 2009); oxidative stability from the upregulation of enzymatic antioxidants (superoxide dismutase, catalase and glutathione reductase) (Demiral and Turkan, 2005); or maintaining ion homeostasis by mobilizing Na+into vacuoles where damage could be contained which is one of the mechanisms observed in Pokkali (Kader and Lindberg, 2005; Ghosh et al, 2011).

    In terms of ion homeostasis, the ratio of Na+and K+in the shoot has been previously established to have a strong inverse correlation with salinity tolerance (Moradi et al, 2007), and exhibit a strong positive linear relationship with SES, which is depicted in Fig.1-c. In this study, the mean Na-K ratio of almost all the test accessions showed non-significant difference with that of FL478. In a previous study, Akundo also accumulates less Na+in its leaves and is thereby protected from osmotic damage (Rahman et al, 2016). Additionally, Chakol, which also registers a shoot Na-K ratio that is significantly lower than FL478, is previously identified to be a popular traditional variety with an expected yield of 3–4 t/hm2by the Bangladesh Rice Research Institute (BRRI) (Hamid and Islam, 1984). Landraces with good yield traits are highly sought due to their ability to be crossed with established lines (Reddy et al, 2017). Notably, AC was the only test accession that obtained a mean Na-K ratio (0.33) that was lower than FL478 (1.75). In a recent study using Na-K ratio to assess salinity tolerance in Bangladeshi landraces, Komol Bhog (0.750), Sona Toly (0.379), Nakraji (0.850), Dud Sail (3.47) and Tal Mugur (1.64) are salt tolerant (Tahjib-Ul-Arif et al, 2018). Among the 44 genotypes screened in this study, three Bangladeshi landraces [Moisdol (0.65), Kuplod (0.76) and Akundo (0.39)] and three Philippine landraces [AC (0.33), Anangka (0.50) and Betalga (0.43)] showed similar values. Dud Sail with SES of 3.0 obtained a relatively high Na-K ratio (3.47), similar to Asha, Sada Rupa, Binutiti, Reppeng, IR5494, IR4493-5-5-3, Anangka, Dukab and Binangahon. Differing tolerance mechanisms aside from ion homeostasis (Na-K ratio) might have contributed to this observation although additional studies on measuring such mechanisms have to be carried out.

    Clearly apparent from the physio-morphometric data is the degree of variation relative to the selected indicators of salt stress. Such varied responses are favourable in developing salt-tolerant cultivars since conventional breeding generally employs lines with extreme characters for the target trait as parentals. Such was the case for the salt-tolerant FL478, which was derived from a cross of Pokkali and IR29 (Walia et al, 2005). One way to complement such efforts is by utilizing molecular markers linked with the salinity tolerance major QTL(Mohammadi-Nejad et al, 2008) for genotypic selection. Here, 16 riceSSR markers on chromosome 1, along with the other markers mentioned above, were utilized to determine the allelic diversity of the test population.

    Based on the constructed consensus tree (Fig. 2), Casibon, Kalagnon, Roa and Ikogan, which showed superior RSDR values against FL478, did not cluster close to each other although they did fall under clade II to which FL478 also belongs, hinting at divergent mechanisms for biomass accumulation under salt stress. Four loci, namely RM10772 (chromosome 1), RM178 (chromosome 1), RM171 (chromosome 7), and RM455 (chromosome 10), were common among these three accessions (Fig. 3). RM10772 has already been mapped to theQTL on rice chromosome 1 (Lapitan et al, 2007;Ganie et al, 2014). Similarly, Akundo, Kuplod and AC cluster together under clade II. RM105, RM10764 and RM413 were found to be common loci among them. Lal Moni, Chakol, Chorua Kartiksail and Aguni Kartiksail clustered close to each other forming clade I. These accessions also obtained significantly lower Na-K ratio values compared to FL478. Clade I being clustered away from FL478 indicates that another mechanism could be complimenting the good Na-K transport system of these accessions. In a recent genotyping study on Bangladeshi landraces using SSR markers, AP3206fwas reported to discriminate cultivars in terms of overall salt tolerance and obtain a high PIC value of 0.7 (Rashid et al, 2018).However, in the present study, AP3206f showed a monomorphic banding pattern excluding IR29 which was heterozygous in this locus. This might be attributed to the replicated screening that was carried out in this experiment, generating a homogeneous population of putatively tolerant group.

    Apparent from all these is the sheer complexity of merging physio-morphometric data with the molecular data which shows just how intricate salt tolerance is, both at the phenotypic and genetic level. Additionally, even with the vast number of SSR markers available to screen for salt tolerant varieties, there is yet to be a definitive locus that could clearly and irrefutably distinguish tolerant from sensitive genotypes. It is therefore a crucial task to characterize as many rice genotypes as possible to fully grasp the full extent of this trait.

    Conclusions

    Due to the complexity of salt stress in plants, novel sources of tolerance must be explored to continuously improve critical crops, including rice. From this study, a total of 44 genotypes of varying degrees of tolerance were screened using physio-morphometric markers of salt stress and microsatellite markers in theQTL. Among these landraces, Casibon, Kalagnon, Ikogan, AC, Akundo, Maranao, Betalga, Anangka, Roa, Ampipit and Binutiti outperformed most accessions and performed better than FL478 (tolerant check at IRRI) in terms of biomass and Na-K ratio. Chakol, which had good performance and yield parameters under salt stress based on a previous selection, could be a candidate for breeding programs. Cluster analysis grouped the accessions sharing the most number of FL478 alleles under cluster II, although it was unable to reflect the performance of accessions in terms of RSLD, RSWC and RSDR. Screening for salt-tolerance in these accessions at the maturity stage could greatly contribute to the gene pool of rice relative to salt stress tolerance.

    Acknowledgements

    We humbly thank the International Rice Research Institute for the germplasm samples. We also thank Ms. Meddy Arceta, Ms. Aiza Vispo, Mr. Junrey Amas for their help during the screening phase, and Ms. Kirsten Gonzalvo for her statistical expertise.

    SUPPlemental DATA

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

    Supplemental Fig. 1.Representative gel documentation and banding pattern of AP3206f, RM3421b and RM10793 used to detect the genetic polymorphisms in the 44 test accessions.

    Asfaw K G. 2011. Effects of salinity on seedling biomass production and relative water content of twenty sorghum (L. Moench) accessions., 3(3): 242–249.

    Bhowmik S K, Titov S, Islam M M, Siddika A, Sultana S, Haquel M D S. 2009. Phenotypic and genotypic screening of rice genotypes at seedling stage for salt tolerance., 8(23): 6490–6494.

    Bonilla P S, Dvorak J, Mackill D, Deal K, Gregorio G. 2002. RFLP and SSLP mapping of salinity tolerance genes in chromosome 1 of rice (L.) using recombinant inbred lines.,65(1): 68–76.

    Delaporta S L, Wood J, Hicks J B. 1983. A plant DNA minipreparation: version II, 1(4): 19–21.

    Demiral T, Türkan I. 2005. Comparative lipid peroxidation, antioxidant defence systems and proline content in roots of two rice cultivars in salt tolerance., 53: 247–257.

    Felsentein J. 1989. PHYLIP: Phylogeny inference package (version 3.2)., 5:164–166.

    Flowers T, Flowers S. 2005. Why does salinity pose such a difficult problem for plant breeders?, 78(1): 15–24.

    Ganie S A, Karmakar J, Roychowdhury R, Mondal T K, Dey N. 2014. Assessment of genetic diversity in salt-tolerant rice and its wild relative for ten SSR loci and one allele mining primer ofgene located on 1st chromosome., 300: 1741–1747.

    Ghosh N, Adak M K, Ghosh P D, Gupta S, Sen Gupta D N, Mandal C. 2011. Differential responses of two rice varieties to salt stress., 5: 89–103.

    Gregorio G B, Senadhira D. 1993. Genetic analysis of salinity tolerance in rice (L.)., 86: 333–338.

    Gregorio G B, Senadhira D, Mendoza R D. 1997. Screening rice for salinity tolerance.: IRRI Discussion Paper Series No.22. Los Ba?os, Laguna, the Philippines: International Rice Research Institute.

    Gregorio GB, Senadhira D, Mendoza RD, Manigbas NL, Roxas JP, Guerta CQ. 2002. Progress in breeding for salinity tolerance and associated abiotic stresses in rice., 76(2/3): 91–101.

    Guevarra E, Loresto G C, Jackson M T. 2001. Use of conserved rice germplasm., 124: 51–56.

    Ha S, Vankova R, Yamaguchi-Shinozaki K, Shinozaki K, Tran L S. 2012. Cytokinins: metabolism and function in plant adaptation to environmental stresses., 17(3): 172–179.

    Hamid M A, Islam M R. 1984. Rice cultivation in the tidal swamps of Bangladesh.: Workshop on Research Priorities in Tidal Swamp Rice. Los Ba?os, the Philippines: International Rice Research Institute: 69–77.

    Hoang T M L, Tran T N, Nguyen T K T, Williams B, Wurm P, Bellairs S, Mundree S. 2016. Improvement of salinity stress tolerance in rice: challenges and opportunities., 6(54): 1–23.

    Hoisington D, Khairallah M, Reeves T, Ribaut JM, Skovmand B, Taba S, Warburton M. 1999. Plant genetic resources: What can they contribute toward increased crop productivity?, 96: 5937–5943.

    Horie T, Karahara I, Katsuhara M. 2012. Salinity tolerance mechanisms in glycophytes: An overview with the central focus on rice plants., 5(1): 1–18.

    Iqbal A, Shoajib J. 1993. Evaluation of soil and water salinity for cropping in the coastal region of Bangladesh., 10: 12–18.

    Islam M R, Salam M A, Hassan L, Collard B C Y, Singh R K, Gregorio G B. 2011. QTL mapping for salinity tolerance at seedling stage in rice (L.)., 23(2): 137–146.

    James R A, Rivelli A R, Munns R, Caemmerer S V. 2002. Factors affecting CO2assimilation, leaf injury and growth in salt-stressed durum wheat., 29: 1393–1403.

    Janaguiraman D M, Ramadass R, Devi D D. 2003. Effect of salt stress on germination and seedling growth in rice genotypes., 90(1/2/3): 50–53.

    Javed A S, Khan M F A. 1975. Effect of sodium chloride and sodium sulphate on IRRI rice., 13: 705–710.

    Kader M A, Lindberg S. 2005. Uptake of sodium in protoplasts of salt-sensitive and salt-tolerant cultivars of rice,L. determined by the fluorescent dye SBFI., 56: 3149–3158.

    Khan M S A, Hamid A, Karim M A. 1997. Effect of sodium chloride on germination and seedling characters of different types of rice (L.)., 179: 163–169.

    Kumar K, Kumar M, Kim S R, Ryu H, Cho Y G. 2013. Insights into genomics of salt stress response in rice., 6(1): 27.

    Lafitte H, Ismail A, Bennett J. 2004. Abiotic stress tolerance in rice for Asia: Progress and the future.: Proceedings of a Symposium on the New Directions for a Diverse Planet at the 4th Meeting of the International Crop Science Congress. September 26 to October 1, 2004. Brisbane, Australia: International Crop Science Congress: 1–17.

    Lapitan V C, Brar D S, Abe T, Redo?a E D. 2007. Assessment of genetic diversity of Philippine rice cultivars carrying good quality traits using SSR markers., 57(4): 263–270.

    Lisa L A, Seraj Z I, Elahi C M F, Das K C, Biswas K, Islam M R, Salam M A, Gomosta A R. 2004. Genetic variation in microsatellite DNA, physiology and morphology of coastal saline rice (L.) landraces of Bangladesh., 263: 213–228.

    Lee KS, Choi WY, Ko JC, Kim TS, Gregorio G B. 2003. Salinity tolerance ofandrice (L.) at the seedling stage., 216: 1043–1046.

    Liu K, Muse S V. 2005. Powermarker: An integrated analysis environment for genetic marker analysis., 21(9): 2128–2129.

    Mazid M A, Malek M A, Hossain M. 2016. The BRAC approach to small farmer innovations.: Technological and Institutional Innovations for Marginalized Smallholders in Agricultural Development. Switzerland: Springer International Publishing: 99–112.

    Mishra B, Singh R, Senadhira D. 2003. Advances in breeding salt-tolerant rice varieties.: Khush G S, Brar D S, Hardy B. Advances in Rice Genetics. Proceedings of a Symposium on Rice Genetics at the 4th meeting of the International Rice Genetics Symposium; October 22–27, 2000. Los Ba?os, the Philippines: International Rice Research Institute: 5–7.

    Mohamed A M, Natarajan S, Vanathi D, Ramasamy S, Sathyamoorthi K. 2007. Lowland rice in costal saline soils: A review., 28(4): 245–253.

    Mohammadi-Nejad G, Arzani A, Rezai A M, Singh R K, Gregorio G B. 2008. Assessment of rice genotypes for salt tolerance using microsatellite markers associated with theQTL., 7(6): 730–736.

    Mohammadi-Nejad G, Singh R K, Arzani A, Rezaie A M, Sabouri H, Gregorio G B. 2009. Evaluation of salinity tolerance in rice (L.) genotypes, 4: 199–208.

    Moradi F, Ismail A M. 2007. Responses of photosynthesis, chlorophyll fluorescence and ros-scavenging systems to salt stress during seedling and reproductive stages in rice., 99: 1161–1173.

    Moreno L, Ordo?ez S, Dela Cruz I, Redo?a E. 2003. Relationship of parental genetic diversity with heterosis in two-line and three-line Philippine rice hybrids.: Khush G S, Brar D S, Hardy B. Advances in Rice Genetics. Proceedings of a Symposium on Rice Genetics at the 4th Meeting of the International Rice Research Institute: 12–14.

    Munns R, Tester M. 2008. Mechanisms of salinity tolerance., 59: 651–681.

    Natarajan S, Ganapthy M, Krishnakumar S, Dhanalakshmi R, Saliha B. 2005. Grouping of rice genotypes for salinity tolerance based upon grain yield and Na:K ratio under coastal environment, 1(2): 162–165.

    PhilRice (Philippine Rice Research Institute). 2011. Management of salt-affected soils for rice production., 77: 1–28.

    Qin F, Kodaira K S, Maruyama K, Mizoi J, Tran L S, Fujita Y, Morimoto K, Shinozaki K, Yamaguchi-Shinozaki K. 2011. SPINDLY, a negative regulator of gibberellic acid signaling, is involved in the plant abiotic stress response.,157: 1900–1913.

    Rahman M A, Thomson M J, Shah-E-Alam M, de Ocampo M, Egdane J, Ismail A M. 2016. Exploring novel genetic sources of salinity tolerance in rice through molecular and physiological characterization., 117(6): 1083–1097.

    Rajendran K, Tester M, Roy S J. 2009. Quantifying the three main components of salinity tolerance in cereals, 32(3): 237–249.

    Rashid MM, Imran S, Islam MA, Hassan L. 2018. Genetic diversity analysis of rice landraces (L.) for salt tolerance using SSR markers in Bangladesh., 3(2): 460–466.

    Reddy I N B L, Kim BK, Yoon IS, Kim KH. 2017. Salt tolerance in rice: Focus on mechanisms and approaches, 24(3): 123–144.

    Rines H W, Molnar S J, Tinker N A, Phillips R L. 2006. Cereals and millets.: Kole C. Genome Mapping and Molecular Breeding in Plants. Berlin: Springer: 211–242.

    Roy S K, Patra S K, Sarkar K K. 2002. Studies on the effect of salinity stress on rice (L.) at seedling stage., 6(3): 254–259.

    Nei M, Chesser RK. 1983. Estimation of fication indices and gene diversities., 47(3): 253–259.

    Saxena M K, Pandey U K. 1981. Physiological studies on salt tolerance of ten rice varieties growth and yield aspect., 24: 61–68.

    Schmidt R, Caldanab C, Mueller-Roeberab B, Schippers J H M. 2014. The contribution of SERF1 to root-to-shoot signaling during salinity stress in rice., 9(1): 1–10.

    Sexcion F S H, Egdane J A, Ismail A M, Dionisio-Sese M L. 2009. Morpho-physiological traits associated with tolerance of salinity during seedling stage in rice (L.)., 34(2): 27–37.

    Sirault X R R, James R A, Furbank R T. 2009. A new screeningmethod for osmotic component of salinity tolerance in cereals using infrared thermography., 36(11): 970–977.

    Tahjib-Ul-Arif MD, Sayed MA, Islam MM, Siddiqui MN, Begum SN, Hossain MA. 2018. Screening of rice landraces (L.) for seedling stage salinity tolerance using morpho-physiological and molecular markers., 40(70): 1–12.

    Thomson M J, de Ocampo M, Egdane J, Rahman M A, Sajise A G, Adorada D L, Tumimbang-Raiz E, Blumwald E, Seraj Z I, Singh R K, Gregorio G B, Ismail A M. 2010. Characterizing thequantitative trait locus for salinity tolerance in rice., 3(2/3): 148–160.

    Türkana I, Demiral T. 2009. Recent developments in understanding salinity tolerance., 67: 2–9.

    Umali D L. 1993. Irrigation-Induced Salinity a Growing Problem for Development and the Environment. Washington DC, USA: the Word Bank.

    Verslues P E, Agarwal M, Katiyar-Agarwal S, Zhu J H, Zhu J K.2006. Methods and concepts in quantifying resistance to drought, salt and freezing, abiotic stresses that affect plant water status., 45(4): 523–539.

    Wahhab A. 1961. Salt tolerance of various varieties of agricultural crops at germination stage.: Salinity Problems in the Arid Zone, Proc. Teheran Symp., 14: 185–192.

    Walia H, Wilson C, Condamine P, Liu X, Ismail A M, Zeng L, Wanamaker S I, Mandal J, Xu J, Cui X, Close T J. 2005. Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage., 139: 822–835.

    Xu Y, Crouch J. 2008. Marker-assisted selection in plant breeding: From publications to practice., 48: 391–407.

    Yeo A R, Yeo M E, Flowers S A, Flowers T J. 1990. Screening of rice () genotypes for physiological characters contributing to salinity resistance, and their relationship to overall performance,79(3):377–384.

    Yoshida S, Forno D A, Cock J H, Gomez K A. 1976. Laboratory Manual for Physiological Studies of Rice. 3rd edn. Los Banos, the Philippines: International Rice Research Institute: 83.

    Zeng L H, Shannon M C. 2000. Salinity effects on seedling growth and yield components rice., 40(4): 996–1003.

    Copyright ? 2019, China National Rice Research Institute. Hosting by Elsevier B V

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

    8 May 2018;

    28 September 2018

    Wilson F. Aala Jr (wfaala@up.edu.ph)

    (Managing Editor: Fang Hongmin)

    国产片特级美女逼逼视频| 亚洲怡红院男人天堂| 亚洲欧美精品自产自拍| 国产在线一区二区三区精| 欧美成人a在线观看| 国产一区二区三区av在线| 自拍偷自拍亚洲精品老妇| 精品一区二区免费观看| 欧美另类一区| 亚洲国产欧美人成| 国产成年人精品一区二区| 亚洲伊人久久精品综合| av免费在线看不卡| 国产欧美另类精品又又久久亚洲欧美| 少妇的逼好多水| 嫩草影院精品99| 免费观看无遮挡的男女| 插阴视频在线观看视频| 热99国产精品久久久久久7| 亚洲人与动物交配视频| 亚洲精品456在线播放app| 熟妇人妻不卡中文字幕| 老师上课跳d突然被开到最大视频| 一级毛片aaaaaa免费看小| 欧美区成人在线视频| 国内揄拍国产精品人妻在线| 王馨瑶露胸无遮挡在线观看| av黄色大香蕉| 亚洲av日韩在线播放| 高清日韩中文字幕在线| 日韩免费高清中文字幕av| 别揉我奶头 嗯啊视频| 黄色配什么色好看| 亚洲美女视频黄频| 最近手机中文字幕大全| 午夜精品一区二区三区免费看| 最近中文字幕2019免费版| 搡女人真爽免费视频火全软件| 亚洲欧美日韩卡通动漫| 日本三级黄在线观看| 亚洲国产日韩一区二区| 国内揄拍国产精品人妻在线| 欧美zozozo另类| 男人和女人高潮做爰伦理| 丰满少妇做爰视频| 在现免费观看毛片| 中文字幕制服av| 另类亚洲欧美激情| 免费观看在线日韩| 国产精品麻豆人妻色哟哟久久| 免费av观看视频| 亚洲成人av在线免费| 人妻夜夜爽99麻豆av| 99热国产这里只有精品6| 国产精品三级大全| 午夜福利视频精品| 天堂俺去俺来也www色官网| 韩国av在线不卡| 麻豆乱淫一区二区| 成人鲁丝片一二三区免费| 亚洲图色成人| 人人妻人人看人人澡| 国产精品一区二区三区四区免费观看| 精品熟女少妇av免费看| 免费高清在线观看视频在线观看| 又大又黄又爽视频免费| 精品午夜福利在线看| 超碰av人人做人人爽久久| 亚洲国产最新在线播放| 一级黄片播放器| 久久精品国产鲁丝片午夜精品| 日韩成人伦理影院| 亚洲第一区二区三区不卡| 各种免费的搞黄视频| 亚洲精品国产av成人精品| 在线观看人妻少妇| 男的添女的下面高潮视频| 国产精品三级大全| 干丝袜人妻中文字幕| 欧美人与善性xxx| 日本与韩国留学比较| 偷拍熟女少妇极品色| 国产精品一二三区在线看| 亚洲国产精品999| 亚洲成人精品中文字幕电影| av卡一久久| 日本三级黄在线观看| 99精国产麻豆久久婷婷| 大片电影免费在线观看免费| 狂野欧美激情性bbbbbb| 亚洲,欧美,日韩| av免费观看日本| 女人十人毛片免费观看3o分钟| 国产91av在线免费观看| 亚洲国产高清在线一区二区三| 各种免费的搞黄视频| 亚洲精品日本国产第一区| 国产伦精品一区二区三区四那| 亚洲一区二区三区欧美精品 | 国产精品久久久久久精品电影小说 | 国产淫片久久久久久久久| 亚洲三级黄色毛片| 精品人妻偷拍中文字幕| 视频中文字幕在线观看| av黄色大香蕉| 国产乱人视频| 久久久久久久精品精品| 伊人久久国产一区二区| av网站免费在线观看视频| a级毛片免费高清观看在线播放| 三级国产精品片| 国产精品99久久久久久久久| 婷婷色综合大香蕉| 国产亚洲av嫩草精品影院| 国产日韩欧美在线精品| 成人毛片60女人毛片免费| 久久久久久国产a免费观看| 久久人人爽人人片av| 麻豆成人av视频| 在线观看一区二区三区| 日本一二三区视频观看| 欧美变态另类bdsm刘玥| 深夜a级毛片| 性插视频无遮挡在线免费观看| 精品人妻一区二区三区麻豆| 22中文网久久字幕| 观看美女的网站| 国产91av在线免费观看| av在线蜜桃| 亚洲精品国产色婷婷电影| 欧美丝袜亚洲另类| 大香蕉久久网| 国产欧美另类精品又又久久亚洲欧美| 少妇人妻精品综合一区二区| 白带黄色成豆腐渣| 精品久久久久久久久av| 成年免费大片在线观看| 三级男女做爰猛烈吃奶摸视频| 在线观看三级黄色| 热re99久久精品国产66热6| 色综合色国产| 国产国拍精品亚洲av在线观看| 亚洲成人久久爱视频| 尾随美女入室| 99热网站在线观看| 国产 精品1| 高清视频免费观看一区二区| 日韩大片免费观看网站| 色视频www国产| 99久久九九国产精品国产免费| av免费在线看不卡| 免费在线观看成人毛片| av福利片在线观看| 精品一区在线观看国产| 亚洲图色成人| 亚洲精品日韩在线中文字幕| 精品久久久久久久末码| 免费人成在线观看视频色| av在线观看视频网站免费| 国产免费视频播放在线视频| 美女被艹到高潮喷水动态| 日产精品乱码卡一卡2卡三| 国产av国产精品国产| 国产综合懂色| 2021天堂中文幕一二区在线观| 一个人观看的视频www高清免费观看| av在线老鸭窝| av播播在线观看一区| 精品久久久久久久末码| 日韩欧美一区视频在线观看 | av在线蜜桃| 深爱激情五月婷婷| 亚洲天堂国产精品一区在线| 只有这里有精品99| 国产亚洲精品久久久com| 日韩在线高清观看一区二区三区| 国产成人91sexporn| 下体分泌物呈黄色| .国产精品久久| 国产精品三级大全| 国产在线一区二区三区精| 午夜福利高清视频| 国产视频首页在线观看| 日本一本二区三区精品| 夜夜看夜夜爽夜夜摸| 久久久久久久午夜电影| 久久精品国产鲁丝片午夜精品| 午夜福利高清视频| 欧美日韩国产mv在线观看视频 | 婷婷色麻豆天堂久久| 国产成人免费无遮挡视频| 国产成人精品一,二区| 国产极品天堂在线| 午夜老司机福利剧场| 成人欧美大片| 亚洲一区二区三区欧美精品 | 国产精品不卡视频一区二区| 久久久午夜欧美精品| a级毛片免费高清观看在线播放| 在线免费观看不下载黄p国产| 80岁老熟妇乱子伦牲交| 五月伊人婷婷丁香| 少妇熟女欧美另类| 日本-黄色视频高清免费观看| 噜噜噜噜噜久久久久久91| 亚洲不卡免费看| 久久久亚洲精品成人影院| 国产黄片视频在线免费观看| a级一级毛片免费在线观看| 五月伊人婷婷丁香| 久久人人爽人人片av| 老司机影院毛片| 丰满乱子伦码专区| 国产亚洲av片在线观看秒播厂| 亚洲精品日韩在线中文字幕| 老司机影院成人| 内射极品少妇av片p| 自拍欧美九色日韩亚洲蝌蚪91 | 黑人高潮一二区| 国产亚洲精品久久久com| 99久久精品一区二区三区| 在线观看一区二区三区| 国产69精品久久久久777片| 日本免费在线观看一区| 女人十人毛片免费观看3o分钟| 国产淫语在线视频| 久久精品久久久久久噜噜老黄| 亚洲成人久久爱视频| 成人亚洲精品一区在线观看 | 国产免费视频播放在线视频| 日韩欧美一区视频在线观看 | 久久国产乱子免费精品| 精品人妻熟女av久视频| 欧美日韩综合久久久久久| 激情 狠狠 欧美| 涩涩av久久男人的天堂| 久久99蜜桃精品久久| 亚洲最大成人手机在线| 亚洲经典国产精华液单| 国产综合精华液| 男女那种视频在线观看| 国产一区二区三区av在线| 一个人观看的视频www高清免费观看| 老师上课跳d突然被开到最大视频| 成年女人看的毛片在线观看| 国产黄片视频在线免费观看| 国产又色又爽无遮挡免| 日本午夜av视频| 久久久久九九精品影院| 一级爰片在线观看| 日韩欧美 国产精品| 国产黄a三级三级三级人| 久久久久久久国产电影| 国产亚洲5aaaaa淫片| 日韩免费高清中文字幕av| 日韩国内少妇激情av| 高清av免费在线| 蜜桃久久精品国产亚洲av| 2022亚洲国产成人精品| 国产精品一区www在线观看| 国产精品国产av在线观看| 亚洲一区二区三区欧美精品 | av.在线天堂| av免费在线看不卡| 高清毛片免费看| 日产精品乱码卡一卡2卡三| 亚洲四区av| 精品酒店卫生间| 搡老乐熟女国产| 天堂网av新在线| 日韩不卡一区二区三区视频在线| 如何舔出高潮| 亚洲自拍偷在线| 熟女人妻精品中文字幕| 国产免费一级a男人的天堂| 卡戴珊不雅视频在线播放| 精品久久久久久久人妻蜜臀av| 七月丁香在线播放| 高清在线视频一区二区三区| 搡女人真爽免费视频火全软件| 韩国高清视频一区二区三区| 午夜福利在线在线| 日韩 亚洲 欧美在线| 欧美区成人在线视频| 美女xxoo啪啪120秒动态图| 成人亚洲精品av一区二区| 欧美+日韩+精品| 精品99又大又爽又粗少妇毛片| 人人妻人人爽人人添夜夜欢视频 | 亚洲最大成人av| 免费看日本二区| 久久久久久久久久久免费av| 在线观看三级黄色| 免费黄色在线免费观看| av在线亚洲专区| 欧美日韩视频精品一区| 啦啦啦中文免费视频观看日本| 蜜桃亚洲精品一区二区三区| 大陆偷拍与自拍| 国产成年人精品一区二区| 亚洲欧美中文字幕日韩二区| 中国美白少妇内射xxxbb| 日本午夜av视频| 嫩草影院入口| 涩涩av久久男人的天堂| 99视频精品全部免费 在线| 一边亲一边摸免费视频| 欧美 日韩 精品 国产| 日韩av在线免费看完整版不卡| 午夜福利视频1000在线观看| 熟女电影av网| 亚洲精品影视一区二区三区av| 特级一级黄色大片| 国产精品精品国产色婷婷| 成人欧美大片| 国产一区二区亚洲精品在线观看| a级一级毛片免费在线观看| 777米奇影视久久| 最近最新中文字幕大全电影3| 亚洲美女视频黄频| 直男gayav资源| 国产一区二区在线观看日韩| 18禁动态无遮挡网站| 国产精品偷伦视频观看了| 少妇 在线观看| 91久久精品电影网| 99久久精品国产国产毛片| 国产老妇女一区| 久久韩国三级中文字幕| 久久久久久久久久久免费av| 亚洲婷婷狠狠爱综合网| 国产精品一区二区在线观看99| 久久99热这里只频精品6学生| 91精品一卡2卡3卡4卡| 欧美成人a在线观看| 国产成人一区二区在线| 久久人人爽av亚洲精品天堂 | 亚洲精品自拍成人| 在线免费观看不下载黄p国产| 男的添女的下面高潮视频| 久久久久网色| 日韩欧美 国产精品| 午夜老司机福利剧场| av福利片在线观看| 久久久久九九精品影院| 欧美xxⅹ黑人| 亚洲国产色片| 日韩一本色道免费dvd| 在线 av 中文字幕| 成人毛片60女人毛片免费| 日本欧美国产在线视频| 91久久精品电影网| 国产精品99久久99久久久不卡 | 久久人人爽人人片av| 男女啪啪激烈高潮av片| 最近的中文字幕免费完整| 天天躁夜夜躁狠狠久久av| 午夜爱爱视频在线播放| 亚洲激情五月婷婷啪啪| 天堂俺去俺来也www色官网| av在线观看视频网站免费| 免费在线观看成人毛片| 国产精品无大码| 亚洲成人一二三区av| 国产精品成人在线| 亚洲精品中文字幕在线视频 | 人妻制服诱惑在线中文字幕| 色播亚洲综合网| www.av在线官网国产| 国产黄片视频在线免费观看| 亚洲色图av天堂| 一个人观看的视频www高清免费观看| 少妇高潮的动态图| 好男人视频免费观看在线| 夫妻性生交免费视频一级片| 插阴视频在线观看视频| 1000部很黄的大片| 少妇人妻精品综合一区二区| 自拍偷自拍亚洲精品老妇| 性色av一级| 国产午夜精品久久久久久一区二区三区| 色视频www国产| 91aial.com中文字幕在线观看| 99热这里只有是精品在线观看| 国产极品天堂在线| 久热久热在线精品观看| 精品国产三级普通话版| videos熟女内射| 最近最新中文字幕免费大全7| 国产高清三级在线| 一级a做视频免费观看| 可以在线观看毛片的网站| 久久久精品94久久精品| 伦精品一区二区三区| 日本爱情动作片www.在线观看| 国产成人精品福利久久| 午夜老司机福利剧场| 亚洲国产高清在线一区二区三| 亚洲最大成人手机在线| 男女啪啪激烈高潮av片| 免费看av在线观看网站| 搡女人真爽免费视频火全软件| 亚洲精品一区蜜桃| 国产一区亚洲一区在线观看| 日韩欧美一区视频在线观看 | 精品午夜福利在线看| 精品视频人人做人人爽| 搡老乐熟女国产| 欧美+日韩+精品| 18禁在线播放成人免费| 免费人成在线观看视频色| 天天躁夜夜躁狠狠久久av| 日韩欧美 国产精品| 精品少妇久久久久久888优播| 日本三级黄在线观看| 国产乱人视频| 久久久久九九精品影院| 欧美激情国产日韩精品一区| 亚州av有码| 亚洲成人一二三区av| 成人毛片60女人毛片免费| 日本av手机在线免费观看| 夜夜爽夜夜爽视频| 日韩伦理黄色片| 免费黄色在线免费观看| 天堂中文最新版在线下载 | 亚洲精品国产av成人精品| 在线精品无人区一区二区三 | 最近最新中文字幕免费大全7| 国产v大片淫在线免费观看| 蜜臀久久99精品久久宅男| 久久精品久久久久久久性| 黄片wwwwww| 伊人久久国产一区二区| 2018国产大陆天天弄谢| 少妇的逼水好多| 国产精品精品国产色婷婷| 久久鲁丝午夜福利片| 毛片一级片免费看久久久久| 欧美日韩一区二区视频在线观看视频在线 | 国产久久久一区二区三区| 亚洲av成人精品一二三区| 男人狂女人下面高潮的视频| 国产极品天堂在线| 在线免费十八禁| 欧美日韩综合久久久久久| 午夜免费鲁丝| 午夜视频国产福利| 午夜激情久久久久久久| 欧美高清性xxxxhd video| 汤姆久久久久久久影院中文字幕| 日本欧美国产在线视频| 日本免费在线观看一区| 亚洲国产色片| 深夜a级毛片| 中文字幕免费在线视频6| 最近最新中文字幕免费大全7| 成人欧美大片| 最近中文字幕高清免费大全6| .国产精品久久| 久久久色成人| 黄片无遮挡物在线观看| 直男gayav资源| 国产成人精品婷婷| 麻豆成人av视频| 国产精品久久久久久av不卡| 我的老师免费观看完整版| 国产乱来视频区| 日韩中字成人| 天堂网av新在线| 听说在线观看完整版免费高清| 国产精品.久久久| 日本欧美国产在线视频| 日本与韩国留学比较| 免费黄网站久久成人精品| 黑人高潮一二区| 亚洲精华国产精华液的使用体验| 91精品伊人久久大香线蕉| 日本黄大片高清| 青春草国产在线视频| 26uuu在线亚洲综合色| av在线老鸭窝| 听说在线观看完整版免费高清| 99re6热这里在线精品视频| 久久6这里有精品| av天堂中文字幕网| 色婷婷久久久亚洲欧美| 日日啪夜夜撸| 精品99又大又爽又粗少妇毛片| 国产成人91sexporn| 在线观看一区二区三区激情| 色哟哟·www| 亚洲国产精品999| av在线播放精品| 亚洲经典国产精华液单| 亚洲av免费高清在线观看| av网站免费在线观看视频| 韩国av在线不卡| 哪个播放器可以免费观看大片| 九九爱精品视频在线观看| 丰满人妻一区二区三区视频av| 成人欧美大片| 青青草视频在线视频观看| 大香蕉久久网| 美女cb高潮喷水在线观看| 成人毛片a级毛片在线播放| 国产乱来视频区| 3wmmmm亚洲av在线观看| 国产综合懂色| 国国产精品蜜臀av免费| 国产精品99久久久久久久久| 精品久久久久久久末码| 中文在线观看免费www的网站| 精品午夜福利在线看| 人人妻人人看人人澡| 嫩草影院新地址| 久久女婷五月综合色啪小说 | 欧美一级a爱片免费观看看| 成人特级av手机在线观看| 精品久久久久久久人妻蜜臀av| 国产综合精华液| 亚洲av二区三区四区| 禁无遮挡网站| 国产黄频视频在线观看| 日本三级黄在线观看| 亚洲国产欧美在线一区| 亚洲欧美成人精品一区二区| 卡戴珊不雅视频在线播放| 国产免费一级a男人的天堂| 亚洲成人一二三区av| 亚洲国产av新网站| 亚洲三级黄色毛片| 久久综合国产亚洲精品| 久久久久久久久久人人人人人人| 最后的刺客免费高清国语| 91久久精品电影网| 嘟嘟电影网在线观看| 亚洲av男天堂| 国产欧美日韩精品一区二区| 校园人妻丝袜中文字幕| 国产有黄有色有爽视频| 在线免费观看不下载黄p国产| 精华霜和精华液先用哪个| 国产日韩欧美亚洲二区| 欧美国产精品一级二级三级 | 国产 精品1| 国产 一区精品| av黄色大香蕉| av一本久久久久| 国产成人精品福利久久| 亚洲av免费高清在线观看| 黄片无遮挡物在线观看| av黄色大香蕉| 丰满乱子伦码专区| 色综合色国产| 国产成人aa在线观看| 九九久久精品国产亚洲av麻豆| 人人妻人人看人人澡| 午夜老司机福利剧场| 日本与韩国留学比较| 国产日韩欧美在线精品| 免费观看无遮挡的男女| 亚洲激情五月婷婷啪啪| videos熟女内射| 视频中文字幕在线观看| 纵有疾风起免费观看全集完整版| 免费黄网站久久成人精品| 亚洲精品国产成人久久av| 美女脱内裤让男人舔精品视频| 久久精品久久久久久久性| 在线观看三级黄色| 伦理电影大哥的女人| 51国产日韩欧美| 免费av毛片视频| 亚洲,欧美,日韩| 国产老妇伦熟女老妇高清| 美女主播在线视频| 波多野结衣巨乳人妻| 18禁裸乳无遮挡免费网站照片| av播播在线观看一区| 中文字幕久久专区| 亚洲精品影视一区二区三区av| 免费高清在线观看视频在线观看| 亚洲人成网站高清观看| 69av精品久久久久久| 神马国产精品三级电影在线观看| 久久久a久久爽久久v久久| 久久精品久久精品一区二区三区| 成人毛片a级毛片在线播放| 伊人久久国产一区二区| 在线观看美女被高潮喷水网站| 永久免费av网站大全| 别揉我奶头 嗯啊视频| 亚洲av电影在线观看一区二区三区 | 观看美女的网站| av网站免费在线观看视频| av在线app专区| 午夜免费鲁丝| 久久久亚洲精品成人影院| 最近最新中文字幕大全电影3| av国产精品久久久久影院| 亚洲精品自拍成人| 亚洲av成人精品一区久久| 久久97久久精品| 91午夜精品亚洲一区二区三区| 最近最新中文字幕大全电影3| 国产av国产精品国产| 男人舔奶头视频| 免费少妇av软件| 亚洲婷婷狠狠爱综合网| 精品人妻视频免费看| 小蜜桃在线观看免费完整版高清| 婷婷色综合www| 免费播放大片免费观看视频在线观看| 亚洲精华国产精华液的使用体验|