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

    Geographic patterns of genetic diversity and fertility restoration ability of Asian and African origin pearl millet populations

    2022-03-30 08:51:34SudarshanPatilShashiKumarGupta
    The Crop Journal 2022年2期

    K. Sudarshan Patil, Shashi Kumar Gupta

    a International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Patancheru, Hyderabad 502324, Telangana, India

    b Professor Jayashankar Telangana State Agricultural University (PJTSAU), Rajendranagar, Hyderabad 500030, Telangana, India

    Keywords:CMS systems Fertility restoration Morphological diversity Pearl millet Populations

    ABSTRACT A set of 45 diverse Asian and African origin/bred populations of pearl millet were investigated multi-locationally for morphological traits including grain yield and assessed for their fertility restoration ability on three cytoplasmic male sterility (CMS) systems in two contrasting seasons. Significant genetic variation was found for all the yield linked traits.Multivariate analysis for grain yield and its component traits grouped these populations into five clusters.Most of the Asian origin populations grouped into single cluster (Cluster III) and separated from the rest of the African origin populations indicating the presence of correspondence between clustering pattern and geographical origin of the populations. Clusters dominated with Asian populations (As-As) had linkages with early flowering, short plant height, more number of tillers, small and thin panicles, small seeded and low grain yield compared to the clusters dominated by African populations (Af-Af) or African origin and Asian bred (Af-As) populations.Genetically related populations having common parentage were found grouped in same clusters.Fertility restoration/maintainer frequency of 45 populations on three diverse CMS systems revealed that overall fertility restoration frequency was highest for A1(86%)followed by A4(37%)and for A5(7%)CMS system.Five populations were identified as potential sources for developing maintainer lines for all three CMS systems and eight populations were identified specifically for A4 and A5 CMS systems. A set of 11 and four populations were identified for restorer line development exclusively for A1 and A5 CMS system,respectively. Six populations were identified for the development of dual restorers for both A1 and A4 CMS system.

    1. Introduction

    Pearl millet(Pennisetum glaucum(L.)R.Br.)is an important food and fodder crop in the arid and semi-arid tropics of Asia and Africa.It serves as a staple cereal and most common source of energy and micronutrients for >90 million people living in the arid and semiarid regions. Pearl millet originated and was domesticated in Sub-Saharan region of West and Central Africa,and later migrated to eastern Africa,semi-arid regions of South Asia,and other parts of the world. ICRISAT (International Crops Research Institute for the Semi-Arid Tropics) conserves the world’s largest collection of 22,888 pearl millet germplasm accessions from 51 countries with 6610 accessions from India [1] and 13,185 accessions from Africa[2]. India and Africa together contribute ~98% of the world pearl millet cultivation. Pearl millet cultivation in Africa (~22 Mha) is mainly based on open pollinated varieties (OPVs), landraces and improved populations. Whereas, in India, hybrids occupy about 70% (~5 Mha) of the area under pearl millet; the rest (~2.5 Mha)is under OPVs or landraces [3]. The pearl millet productivity in India increased at the rate of 3% per annum during 1990–2017[4] due to exploitation of heterosis through hybrid development.In contrast, the productivity has not significantly changed in the African continent in the last three decades from 1988(691 kg ha-1)to 2018(718 kg ha-1),though area has increased from 15.8 to 22.1 Mha[5].Hence,it is necessary to develop new series of high yielding cultivars by involving under-utilized genetic diversity.

    Pearl millet breeding in India during 1930 s started with the development of OPVs from locally adapted material through mass selection and progeny evaluation. But, they provided marginal improvement in the yield level. Later, with the establishment of ICRISAT during 1970s resulted into the development of some high yielding OPVs following different population improvement programs utilizing diverse range of Asian and African germplasms. It also leads to the development of several trait-based composites[e.g.early composite(EC),medium composite(MC),late composite(LC), smut-resistant composite (SRC), high-tillering composite(HTC), bold-seeded composite (BSC), dwarf composite (DC), high head volume composite (HHVC) and many others] to support the breeding programs. Several high yielding OPVs such as WC-C75,ICMV 221, ICTP 8203, HC 4, Raj 171, JBV 2 and JBV 3, HC 20, RCB 2, CZP 9802 and Pusa 266 were developed from these trait-based composites through recurrent selection [6]. However, many of these composites/OPVs, breeding populations and germplasms were not phenotypically well characterized which limits their utilization in current breeding programs in development of the new range of high yielding cultivars. Initially, during the beginning of hybrid breeding era in pearl millet(1960s),trait based composites and populations were observed for traits and fertility/sterility and were involved in seed parent or restorer parent development programs, but this information was never documented.

    Understanding the structure of diversity and identification of distinct materials with desirable traits provides the foundation for effective and sustained pearl millet population breeding, and hybrid development in this cross pollinated crop. Access to wellcharacterized, adapted germplasm and knowledge about the geographic patterns of genetic variation for various traits in a target region are important prerequisites for successful plant breeding[7]. Previous diversity studies involving diverse germplasms reported wide range of valuable diversity for various agronomic traits,stress tolerance and nutritional traits in African populations[8–12],and in Asian germplasms[1,13,14]and also in ICRISAT bred accessions [15,16]. These studies emphasized that phenotypic characterization is the first step for the assessment, description and classification of germplasm collections to enhance their utilization in pearl millet breeding. The present study systematically characterizes the promising Asian and African germplasms bred/originated in Asia and Africa in the last 40 years and documents their morphological traits to use them in current pearl millet breeding programs.

    In addition,these populations were also assessed for their ability to restore fertility or sterility maintenance on A1, A4, and A5CMS systems. A1CMS system is the most commonly used CMS and most of the hybrids cultivated in farm are based on A1CMS till date. The A4and A5CMS systems were found to be a bit more stable than the A1CMS system. However, commercial hybrids couldn’t be developed using A4and A5CMS system due to nonavailability of suitable restores in these two CMS systems [17].Hence, the present study aimed to characterize diverse Asian and African origin/bred pearl millet germplasm for agromorphological traits and to identify promising restorers/maintainer populations for three CMS systems to derive new restorer and maintainer lines to use them as parents in future crop improvement programs.

    2. Materials and methods

    A set of 45 pearl millet populations maintained by pearl millet breeding program of ICRISAT, Patancheru, Hyderabad, India were investigated. They involved genetic materials of African and Asian origins and were developed either at ICRISAT or in collaboration with partners. These populations are expected to be in Hardy-Weinberg equilibrium (HWE) as they were maintained through different cycles of random mating in isolation at ICRISAT. These 45 populations involved 10 Asian origin-Asian bred (abbreviated as As-As), 7 African origin-African bred (Af-Af), 3(Asian × African) origin-Asian bred [(As × Af)-As], and 25 African origin-Asian bred (Af-As) populations. Pedigree details of these populations are given in Table S1. This set also involved some of popular OPVs released in Africa such as, GB 8735, SOSAT-C88,Ugandi and ICMV IS 92222 and some of popular OPVs cultivated in Asia(India)such as,ICTP 8203,CZP 86,Raj 171,ICMV 221,ICMV 155, and WRajPop.

    Fertility restoration of these 45 populations was assessed on A1,A4and A5CMS systems.A-lines of A1,A4and A5cytoplasm in three different genetic backgrounds representing different flowering duration such as ICMA 94222 (early flowering habit), ICMA 89111 (medium flowering) and ICMA 91555 (late flowering) were used in order to synchronize the flowering duration of the populations for facilitating the ease of crossing program[18].A set of 135 top crosses (45 populations × 3 CMS systems) were developed by crossing the male-sterile line with bulked pollen from 20 to 25 random plants of each population during rainy (July–October) season of 2015. Some populations being photo-sensitive, so couldn’t be crossed in rainy season of 2015, hence were planted again and crossed in summer (March–June) season of 2016. Soon after pollination, the crossed panicles were covered with parchment paper bag to avoid foreign pollen contamination and labeled properly.The crossed panicles were left covered with parchment paper bags until harvest, at physiological maturity panicles were harvested,sundried and threshed to collect top cross hybrid seeds.

    2.1. Experimental layout

    2.1.1. Morphological diversity

    The experiment was laid out in alpha-lattice design with three replications at two locations in India during rainy season of 2015.The experimental locations were ICRISAT, Patancheru, Telangana(17°30′N, 78°27′E, and 545 m above sea level) and Regional Agricultural Research Station (RARS), Palem, Telangana (16°53′N,78°23′E,and 545 m above sea level).The main block for each replication had eight sub-blocks with a size of six entries per block.Each population was sown in six rows, 4 m in length with 75 cm inter-row spacing and an intra-row spacing of 15 cm. All the recommended agronomic practices were followed at both the locations to support good crop growth. There were no major incidences of diseases and pests. All the panicles were harvested at physiological maturity and threshed after sun drying for 10–15 days.

    Observations on seven phenotypic traits were recorded. Data was recorded on 20 randomly selected plants from each population to appropriately represent the existing diversity present in the populations since populations are heterogeneous and heterozygous in nature.Phenotypic traits such as plant height(cm),number of productive tillers per plant,panicle length(cm)and panicle girth(cm) were recorded on plant basis. Days to 50% flowering, grain yield (kg plot-1) and 1000-grain weight (g) were recorded on plot basis. Grain yield (kg plot-1) was later converted to kg ha-1for analysis.

    2.1.2. Fertility restoration

    Field evaluation trial of 135 top cross hybrids was conducted at ICRISAT, Patancheru during two contrasting seasons, the rainy(July–October) season of 2016 and summer (March–June) season of 2017.Each hybrid was evaluated in four-row plots of four meter length with about 100–120 plants per plot. One of the tillers of each plant was selfed with parchment paper bag at boot leaf stage.At physiological maturity, bags were opened and observed physically for seed setting to confirm fertile/sterile reaction.Those plant panicles with complete seed set were scored as male-fertile(F)and with no seed set as male sterile (S). Frequency of restoration or maintenance in percentage were calculated using the following formulas:

    2.2. Statistical analysis

    2.2.1. Analysis of variance (ANOVA)

    Combined ANOVA was carried out using PROC MIXED program[19], considering location, genotypes and replication as fixed and block as random. In order to pool the data across two locations and to make the error variance homogeneous, individual location variances were estimated and modeled to error distribution using residual maximum likelihood(REML)procedure.The analysis procedures used the linear mathematical models with the general notations of b = lattice blocks, r = number of replications and g = number of genotypes. The mathematical model followed for analysis of variance for the alpha-lattice design was:

    where,μ is the grand mean,eiis the fixed effect of location i,glis the fixed effect of genotype l,(e/r)ijis the fixed effect of replication j nested with in location i, (e/r/b)ijkis the random effect of block k nested with in replication j and location i and is ~NID(0,σ2b),(eg)ilis the fixed effect of the interaction between genotype l in location i;εijklis the random residual effect and is ~NID(0,σ2ε).

    2.2.2. Multivariate analysis

    Multivariate analysis suggested by Mahalanobis [20] was used to assess phenotypic diversity in the material. A hierarchical cluster analysis was performed using population means of all the variables for all the 45 populations based on dissimilarity matrix according to Ward’s minimum variance method [21] as per the PROC-CLUSTER program in SAS v. 9.4 [19]. This method first computes a matrix based on Euclidean distances among group means and produces a dendrogram depicting successive fusion of individuals, and conclude at the stage in which all the individuals of the same group form a cluster [22].

    3. Results and discussion

    3.1. Differentiation for morphological traits

    Morphological evaluation of the 45 populations revealed wide range of variation for all the traits(Table 1;Fig.1).Details of mean performance of populations across the locations are provided in Table S2. The wide range of performance per se for all the traits under investigation provides an opportunity to select populations for different traits in breeding programs as per requirements/traits performance of different pearl millet growing regions.

    Combined ANOVA revealed significant variation for all the traits in the populations evaluated in this study (Table 2). Significant variation due to locations (environments) was observed for all the traits except for days to 50%flowering,indicating that the populations were evaluated under diverse environments.Location × genotype interaction was found non-significant for most of the traits except for plant height and number of productive tillers, indicating that the ranking of populations based on their performance remains same for other traits,like days to 50%flowering,panicle length,panicle girth,1000-grain weight and grain yield across the locations. Previous reports by Murthy et al. [23], Ouendeba et al. [24], and Ousmane et al. [9] also reported significant variation for location, and genotypic effect for most of the traits,and non-significant variation for interaction effect(location × genotype) while investigating diversity among wide set of germplasms. This non-significant interaction effect might be due to heterogeneous nature of populations which provides them inherent buffering capacity to adapt to different environments[25–27].On the contrary,other phenotypic diversity studies in pearl millet germplasms [8,10,11,14] reported significant variation for location × genotype interaction.

    3.2. Diversity analysis

    The multivariate analysis of 45 populations partitioned them into five different clusters (Fig. 2). Of these five clusters, two clusters(designated as Cluster I and II) had majority of the Af-As populations (24 populations, 96%) and third cluster (designated as Cluster III)had majority of Asian populations(6 As-As populations,60%); fourth cluster (designated as Cluster IV) had only 3 populations of which one was African OPV SOSAT-C88 and other two were(As × Af)-As and Af-As populations. Remaining one cluster (designates as Cluster V) had only one Af-Af population ICMV IS 92222.

    Fig. 1. Box plots showing variation for seven phenotypic traits evaluated for 45 pearl millet populations, combined across locations. (+indicates mean values; — indicates median values; ◆indicates extreme values).

    Among five clusters, the first cluster (Cluster I) was the largest and composed of 20 populations (44%), out of which 14 were Af-As populations and 3 each of As-As and Af-Af populations. Cluster I had three sub-clusters.Sub-cluster Ia had majority of populations developed from SRC bred at Asia using African germplasms. The SRC was formed at ICRISAT in 1978 using 37 parents with lowsusceptible to smut which were derived from an intercross of a wide range of breeding materials selected under disease pressure[28]. Sub-cluster Ib had two Af-As population developed using SRC and New elite Composite (NELC) and two As-As populationsdeveloped using Rajasthan(India)material.Grouping of these two As-As populations with the Af-As populations might be due to their closeness for phenotypic traits as observed for high yielding ability(>2 t ha-1), tillering ability (~2), panicle girth (2.4 cm) and 1000-grain weight(10 g)(C Ib)(Table S2).Sub-cluster Ic had intermixing of the three populations each from Af-As and Af-Af category, one population each from As-As(WRajPop bmr)and(As×Af)-As(ICMS 7704). Among the three Af-As populations two were developed from SRC material and other one was from bold seeded early composite (BSEC). Grouping of these two SRC derived populations in sub-cluster Ib instead of sub-cluster Ia might be due to the involvement of diverse composites (like D1C, NEC, IVC, NELC, togo, MC,and WC-C75)in addition to SRC,and these populations were characterized by high grain yield(>2 t ha-1)compared to the rest of the SRC derived populations of sub-cluster Ia. Among the rest of the three populations, two were Af-Af populations (NWC C2 and WC C3 bulk), derived from Nigerian World composite [29] and other population (SC1 C4 bulk) was developed using Serere composite from Uganda[30].The intermixing of the different category of populations in common cluster seems to be due to the close similarity in their mean values for multiple phenotypic traits such as medium flowering duration (44 days), low tillering ability (2), medium plant height (173 cm), 1000-grain weight (9.8 g), panicle girth(2.3 cm), and medium grain yield (2049 kg ha-1) (Table 1).

    Table 1 Cluster mean values, overall mean and range of yield and its component traits in 45 pearl millet populations.

    Table 2 Analysis of variance of pearl millet populations for grain yield and its component traits across two locations (Patancheru and Palem).

    Cluster II had 14 populations (31%) having three sub-clusters.Sub-cluster IIa had seven populations of which two were developed from West African material and germplasms from Serere composite (Ugandi and ICMP 87237), Sudan (Sudan I and Sudan II), Zimbabwe (ICMP 96601), and one derived from Western Rajasthan (CZP 86). Though, Ugandi and CZP 86 belongs to Af-Af and As-As population category, respectively, and were found grouped with Af-As populations which might be due to their close similarity in the phenotypic traits, as both of them had medium flowering duration (43 days), medium plant height (<200 cm),number of tillers (~2), panicle length, panicle girth and 1000-grain weight as that of the other populations of this cluster.Sub-cluster IIb and IIc had Af-As populations having large seeds developed using Large Grain Population (LaGraP), togo population(known for large seeds), and BSEC of African origin (known for large and bold seeds). Hence, all the populations of the Cluster II were predominantly composed of bold seed with high mean 1000-grain weight (11.6 g), with taller plants (180 cm), high panicle girth(2.7 cm)and medium flowering duration(42 days)compared to the rest of the population clusters (Table 1). The popular OPV variety, GB 8735 (Af-Af) grouped with Af-As populations in Cluster II due to its close trait similarities with this cluster (high yielding, bold seeds with high 1000-seed weight, smaller plant height, medium panicle thickness with medium to small panicle).This indicated that African germplasms were used as sources in Asian pearl millet breeding programs for developing germplasm having early maturity with large and bold grains,hence Af-As populations(ICMV 94135,ICMV 94132,ICMV 88908,and ICMP 96132)grouped into a common sub-cluster IIb in the present study.Manga[31]also reported that landrace from African regions such as Benin,Burkina Faso, Ghana and Togo were utilized to breed for bold seeded and early-maturing germplasms in India which led to the release of promising OPVs such as ICTP 8203 and ICMV 221.

    Fig.2. Clustering pattern of 45 populations based on morphological traits.Full circle,As-As populations;Dotted circle,Af-Af populations;Rectangle,(As×Af)-As and rest are Af-As populations. As-As,Asian origin and Asian bred; Af-Af,African origin and African bred;(As×Af)-As,(Asian ×African)derived and Asian bred populations, and Af-As,African origin and Asian bred populations. a, b, and c indicates sub-clusters of that particular main cluster.

    Cluster III had 7 populations (16%) which were bred and originated in Asia (As-As).This cluster had populations with mean values of days to flowering (40 days), plant height (154 cm), tiller number (3.2), panicle length (2.11 cm), panicle diameter(2.17 cm), 1000-grain weight (8.74 g) and with mean grain yield(1897 kg ha-1)compared to the rest of the clusters.This indicated that Asian populations were characterized by early flowering,shorter plant height, high tillering, thin panicles and small seed size with low grain yield compared to African populations.Though,the cluster had majority of As-As population(86%),one population EC C6 from (Af × As)-As category grouped in this cluster. This might be due to the pedigree of EC C6 population which involves intercrossing of 41 exotic lines with the large number of Indian entries (151 lines) [32] which might have resulted in similar phenotypic traits and adaptation as that of the As-As populations(Table 2).

    Cluster IV had three and Cluster V had only one population bred using African germplasms.Cluster mean of these populations indicated that these were late in flowering, taller, thicker panicles,medium bold seeded and high yielding. Cluster V had only one African population ICMV IS 92222,it had distinct phenotypic characters, like late flowering (49 days), taller plants (231 cm) with large panicle length (47.5 cm) than rest of all the populations of the four clusters.

    Population category-wise trait mean values also indicated that As-As populations were early in flowering, shorter plant height,higher number of tillers, smaller and thin panicles, smaller seed size and low grain yielders compared to the rest of the three population categories (Table 3). Af-Af population category showed higher mean values for all the seven traits whereas, Af-As and(Af × As)-As populations were found to be intermediate for traitvalues between the Af-Af and As-As populations (Table 3). This indicates that African germplasm helped in improving the yielding ability of Asian germplasms through improved panicle and seed traits.Yadav and Rai[6]and Manga[31]stated that Asian cultivars developed using Iniadi landraces (derived from the Togo population cultivated in West African region)are known for their unique combination of several useful traits, such as early maturity, high seed yield, large seed size, resistance to multiple pathotypes of downy mildew, tolerance to terminal drought, high quality stover yield, and higher levels of grain Fe and Zn content [33,34] have been widely adopted by farmers in India.

    Table 3 Population category wise mean values of morphological traits and fertility restoration (%) on three cytoplasmic male sterility systems.

    Present study clearly indicated that Asian populations were early in flowering having short plant height,more tillering,smaller and thinner panicles, smaller seeds and low yielding as compared to the African origin and bred populations(Table 1).Thus,it seems that these characters were primarily responsible for divergencebetween the Asian (Indian) and African germplasm. The continuous cultivation of pearl millet in Asia (India) in arid regions and poor agronomical environments might have resulted in the selection for earliness and resulting into the erosion of alleles for productivity [35]. Previous study on diversity of world collection of pearl millet landraces conserved at ICRISAT by Upadhyaya et al.[1]also reported early flowering,short plant height and high tillering among the Asian germplasms and very late flowering,tall plant height, long panicles and large seeds among African populations.Murthy et al. [23] also found that African populations have longer panicle,with high yielding,less tillering and late maturity as compared to Indian landraces while characterizing world collection of pearl millet genetic stocks.

    Table 4Fertility restoration (%) of 45 populations on three diverse A1, A4 and A5 cytoplasmic male sterility systems for individual as well as across the seasons.

    Overall,the clustering pattern showed Af-As and As-As populations delineated into clear-cut separate clusters. However, Af-Af populations were found intermixed with Af-As populations. It was noted that the populations related by their parentage clustered together. For instance, all the nine populations developed using SRC lines were grouped in Cluster I and those developed using Nigerian world composite (WC C75, NWC C2, and WC C3 bulk)also grouped in Cluster I.Those populations developed using BSEC(2 populations),LaGraP(3 populations),and Togo population(originated from landraces cultivated in the Togo region of West Africa known for early flowering, limited tillering, broad panicles and very large seed size with high Fe and Zn content) having bold seeded trait were found grouped together in cluster II.This implies the presence of certain level of congruence in the clustering of populations sharing common genetic material involved in them and their trait values in the present study.Wilson et al.[36]also found some correspondence between the geographical collection sites of pearl millet landraces from Central Burkina Faso and their inclusion in particular clusters.

    3.3. Fertility restoration/Sterility maintenance reaction

    Across all the 45 pearl millet populations, the fertility restoration ranged from 44%to 100%with a mean of 86%for A1CMS system;from 1%to 85%with a mean of 37%for the A4; and from 0 to 26% with a mean of 7% for A5CMS system was observed over the two seasons(Table 4).Among three CMS systems,all the 45 populations restored fertility on A1CMS system, of which 18 populations showed >90% of restoration for A1CMS across two seasons.However, for A4and A5CMS system, none of the populations showed >90% of restoration. In case of A4CMS system, 8 populations had >60% of restoration and for A5CMS system the highest restoration observed was 26% in one population (ICMP 96201).Two populations ICTP 8203 and ICMV 94135 showed higher restoration percentage on A4than A1cytoplasm. One population GICKV 97652 (JBV 3) showed almost similar level of restoration percentage on both A1(88%)and A4(85%)cytoplasm.On contrary,CZ-IC 9802 showed 100%restoration on A1CMS system,had maintainer reaction on A4and showed poor restoration on A5(10%)CMS system. The frequency of restorers across all the populations was in the order of A1>A4>A5.

    3.4. Restoration differences in Asian and African populations

    Among four different categories of populations, As-As and(As × Af)-As populations showed restoration (>90%), followed by Af-Af (88%) and Af-As (83%) for A1cytoplasm. In case of A4cytoplasm, Af-As populations had higher restoration (40%) compared to rest of the categories of the population. However, none of the populations of different categories showed satisfactory restoration for A5cytoplasm (Table 3; Fig. 3). Rai et al. [18] also reported higher restoration frequency in Asian populations varying form 26%–94% for the A1CMS system, from 4% to 72% for the A4and from 0 to 13%for the A5system as compared to African populations varying form 26% to 84% for A1CMS system, from 12% to 73% for the A4and from 1% to 12% (except one population with 55%) for A5CMS system.

    Fig. 3. Fertility restoration frequency of male sterile plants (restorers) in top cross hybrids of A1, A4 and A5 made with pearl millet populations based on their origin and adaptation.As-As,Asian origin and Asian bred;Af-Af,African origin and African bred;(As×Af)-As,(Asian×African)derived and Asian bred populations;Af-As,African origin and Asian bred population.

    3.5. Seasonal variation in restoration/maintainer frequency

    Thirty-six populations showed ≥10% of variation in fertility restoration frequency between two contrasting seasons. Among them, five populations for A1CMS, six for A4, eight for A5CMS,six for both A1and A4CMS, one for both A1and A5CMS, nine for both A4and A5and one population (Ugandi) for all the three CMS systems showed seasonal variation(≥10%)in restoration frequency. It indicated that the restoration frequency of the populations was affected by season for all the three CMS systems.However,mean restoration frequency was higher during rainy season than in summer, irrespective of cytoplasmic sources and restorer populations (Table 4). The high fertility restoration/seed set percentage during rainy season as compared to post rainy(summer) has been previously reported by Rai and Hash [37], Rai et al. [18], Lakshmana et al. [38] in pearl millet and Biradar [39]in sorghum. Also, Gupta et al. [40] reported that the enhanced expression of fertility restoration during the rainy season might be due to relatively lower temperature and higher humidity than summer season.

    3.6. Restoration/maintenance ability of populations and their utilization in hybrid breeding

    Seed and restorer parents in pearl millet breeding programs have been bred for a specific set of traits separately in B- and Rline hybrid parental lines development breeding programs. For instance, maintainer or B-lines are generally bred for short height(<100 cm)and larger seed size while restorer or R-lines are generally bred for taller height (150–180 cm), more tillers, relatively small seed size, and profuse pollen production [17]. Results of the present study indicated that A1system provides greatest opportunity for development and diversification of the restorer parents or pollen parents (R-lines) followed by A4. While, A5CMS system with high sterility or maintainer provides higher opportunity for diversification of maintainer lines(B-lines)than other two cytoplasmic male sterile sources A1and A4.

    Male sterile A-line or its counterpart maintainers (B-lines)should have high yielding panicle and seed traits, like large and thick panicle,bold seeds to economize the seed production.In case of A1CMS system most of the populations showed lesser maintainer frequency (0 to 66%) as compared to A4(15% to 99%) and A5(74%to 99%)CMS system.However,the A1CMS system is most widely used in the pearl millet hybrid breeding program due to its highly stable maintenance or restoration ability across agroecologies and seasons across the pearl millet growing regions of the world, and >95% of the hybrids released till date are based on A1CMS system. Owing to its highest restoration found in global germplasm, the hybrid breeding programs have preferred A1CMS system over other (A4and A5) CMS systems, though A1has been found bit lesser stable than A4and A5CMS system[18].This necessitates more efforts on the diversification of maintainers for A1CMS system to derive new stable maintainer (A-/B- line) breeding lines.

    In the present study,five populations were identified for developing maintainers for all the three CMS system(A1,A4,and A5)and eight populations were identified as potential sources for developing maintainers for A4and A5CMS systems(Table 5).These 13 populations had desirable traits such as high grain yield(>2500 kg ha-1),medium to early flowering(39–49 days),medium to high panicle girth (2.2–3.0 cm) with medium to long panicle length (20–29.5 cm) and bold grains (>10 g of 1000-grain weight)for maintainer line (B-line) development. Though the five populations (ICMV 155, ICMV 221, ICMV 94132, ICMP 96601, and ICMP 98107) have low maintainer frequency for A1CMS system (15%–19%) and three populations (RCB–IC 948, Raj 171, and Sudan I)have lower maintainer frequency of 42%–45%for A4CMS,yet they can be subjected to recurrent selection for sterility/maintainer reaction to increase the frequency of maintenance to derive new maintainer lines. Previously, Rai et al. [17] demonstrated this increase in the maintainer frequency among Asian and African populations for different CMS systems through recurrent selection.One cycle of selection in favor of sterility/maintainer reaction in Asian population Raj 171 lead to increase in frequency of maintainers from 18%(C0 bulk)to 98%(C1 bulk)for A1CMS system and 49%(C0 bulk)to 99%(C1 bulk)for A4CMS system.However,all of these 13 populations have higher plant height(>150 cm),hence there is a need to transfer the d2gene (dwarfing gene) through the back cross program(using dwarf plant height d2donor lines)to develop diverse maintainer lines with lesser plant height [41,42].

    Table 5 High grain yield and linked traits in populations for developing maintainers (B-lines) in three different CMS systems.

    Table 6 High restoring populations for developing restorers (R-lines) in three different CMS systems.

    The list of potential populations for developing new restorers for three CMS system are given in Table 6.Eleven populations having high restoration frequency(92%–100%)for A1CMS system were identified as source populations for developing restorer parental lines exclusively for A1CMS system. In addition, six populations were identified as dual-restorers (for A1and A4CMS) having a restoration frequency of 88%–98% for A1and 54%–85% for A4CMS system. The development of dual-restorers which can restore fertility of two different CMS systems(A1and A4CMS)will be advantageous for the hybrid breeding programs. ICRISAT designates around nine restorer lines each year,and 48 are dual restorer lines(for both A1and A4CMS)out of the 126 restorers designated since 2006.The dual-restorer populations identified in the present study will be helpful in developing diverse dual-restorer hybrid parental lines. For A5CMS system, four populations showing >15% restoration were identified as potential restorer source populations.These populations can be subjected to recurrent selection to increase the frequency of restoration and develop restorer parents for respective CMS system.Rai et al.[17]increased the frequency of restorers through recurrent selection in Raj 171.One cycle of selection in Raj 171 population increased the restoration frequency from 82% (C0bulk) to 99% (C1bulk) for A1cytoplasm and 51% (C0bulk) to 96%(C1bulk) for A4CMS system.

    At ICRISAT, four populations such as MRC general bulk, AIMP 92901, JBV 3, and GB 8735 for restorer line development program and HHVBC population in case of B-line development have been extensively utilized till date. Whereas, populations such as ICTP 8203,SOSAT-C88,ICMS 7704,and EC C6 were used in both restorer and maintainer breeding programs. In any crop breeding program it is necessary to maintain the distinct gene pools for maintainer and restorer line development. Hence, in the present study, we identified the distinct populations for their utilization separately in development of maintainer and restorer parents to avoid intermixing of gene pools.Based on the present study,it can be inferred that populations having desirable traits for respective B- or R-line development program provides a great opportunity to derive new diverse restorer and maintainer lines to enhance the genetic diversity.

    4. Conclusions

    To enhance introgression of new germplasm in current pearl millet breeding programs, the present study showed existence of significant amount of diversity among the Asian and African origin/bred populations for different traits like days to 50%flowering,plant height, panicle traits, 1000-grain weight and grain yield.Most of the As-As populations found clustered together and were found to share common traits, such as early flowering, less plant height, high tillering, small thin panicles, small-seeded, and low yielding as compared to the African origin/bred populations.These diverse populations can be used as base populations to derive new breeding lines with desired traits to broaden the genetic base of current pearl millet breeding programs. Fertility restoration frequency (%) was found highest for A1(86%) followed by A4(37%)and least for A5(7%) CMS system in 45 populations under investigation. Populations were identified as potential sources for maintainer and restorer line development for A1, A4and A5CMS systems.These populations can serve as potential sources to diversify the genetic background of the maintainers and restorers in the hybrid pearl millet breeding programs.

    CRediT authorship contribution statement

    K. Sudarshan Pati:Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing - original draft,Writing-review&editing.Shashi Kumar Gupta:Conceptualization, Funding acquisition, Project administration, Resources,Software, Supervision, Validation, Visualization, Writing - original draft, Writing - review & editing.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    The research was financially supported by the CGIAR Research Program on Dryland Cereals Scholarship(CRP-DC)and S.M.Sehgal Foundation Fund,USA for carrying out a doctoral study of K.Sudarshan Patil at ICRISAT.We gratefully acknowledge Dr.C.T.Hash(formerly Principal Scientist, ICRISAT, Niger) for helping with pedigrees of populations involved in the study.

    Appendix A. Supplementary data

    Supplementary data for this article can be found online at https://doi.org/10.1016/j.cj.2021.04.013.

    久久鲁丝午夜福利片| 在线观看人妻少妇| 午夜久久久在线观看| 男女免费视频国产| 女性生殖器流出的白浆| 午夜福利,免费看| 久久精品久久久久久久性| 两个人看的免费小视频| 别揉我奶头~嗯~啊~动态视频 | 国产成人一区二区在线| 好男人电影高清在线观看| 一边摸一边抽搐一进一出视频| 亚洲av片天天在线观看| 久久人人爽av亚洲精品天堂| 操出白浆在线播放| 丝袜在线中文字幕| av国产精品久久久久影院| 国产精品人妻久久久影院| 亚洲国产av新网站| 欧美精品一区二区大全| 国产精品久久久av美女十八| 国产黄色视频一区二区在线观看| 夫妻午夜视频| 女人精品久久久久毛片| 国产色视频综合| 丝袜人妻中文字幕| 免费观看人在逋| 久久av网站| 女人高潮潮喷娇喘18禁视频| 女人爽到高潮嗷嗷叫在线视频| 脱女人内裤的视频| 国产精品人妻久久久影院| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品av久久久久免费| 精品久久久久久久毛片微露脸 | 亚洲av国产av综合av卡| 人人妻人人爽人人添夜夜欢视频| 91字幕亚洲| 秋霞在线观看毛片| 国产成人欧美在线观看 | 亚洲九九香蕉| 18禁国产床啪视频网站| 亚洲成人免费av在线播放| 亚洲欧美中文字幕日韩二区| 一区二区三区四区激情视频| 两个人看的免费小视频| 日韩av不卡免费在线播放| 国产主播在线观看一区二区 | 人体艺术视频欧美日本| 少妇被粗大的猛进出69影院| 男人舔女人的私密视频| 超碰成人久久| 欧美日韩视频高清一区二区三区二| 成人国语在线视频| 汤姆久久久久久久影院中文字幕| 免费高清在线观看日韩| 伊人久久大香线蕉亚洲五| 欧美性长视频在线观看| 欧美性长视频在线观看| 国产欧美亚洲国产| 无限看片的www在线观看| 一级毛片我不卡| 天天躁日日躁夜夜躁夜夜| 91九色精品人成在线观看| 亚洲精品国产av成人精品| 久久99热这里只频精品6学生| 久久综合国产亚洲精品| 黄片播放在线免费| www.精华液| 成人三级做爰电影| 国产精品国产三级专区第一集| 成年动漫av网址| 中文字幕制服av| 国产一级毛片在线| 另类精品久久| 黄色 视频免费看| 青青草视频在线视频观看| 一本—道久久a久久精品蜜桃钙片| 人成视频在线观看免费观看| 人人妻人人澡人人爽人人夜夜| 九草在线视频观看| 美女高潮到喷水免费观看| 老汉色av国产亚洲站长工具| 97在线人人人人妻| 搡老岳熟女国产| avwww免费| 婷婷成人精品国产| 午夜精品国产一区二区电影| 亚洲国产欧美日韩在线播放| av线在线观看网站| 欧美日韩亚洲高清精品| 99久久综合免费| 亚洲欧美日韩高清在线视频 | 精品国产一区二区三区四区第35| 美女视频免费永久观看网站| 97精品久久久久久久久久精品| 夜夜骑夜夜射夜夜干| 少妇猛男粗大的猛烈进出视频| 国产高清国产精品国产三级| 成人三级做爰电影| 最新在线观看一区二区三区 | 大码成人一级视频| 国产视频首页在线观看| 狂野欧美激情性xxxx| 别揉我奶头~嗯~啊~动态视频 | 亚洲欧美色中文字幕在线| 狠狠精品人妻久久久久久综合| 精品少妇内射三级| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品自拍成人| 久久精品亚洲av国产电影网| 欧美精品一区二区免费开放| 日本午夜av视频| 亚洲天堂av无毛| 人体艺术视频欧美日本| 国产在线视频一区二区| 成人午夜精彩视频在线观看| 国产不卡av网站在线观看| 大香蕉久久网| 国产精品.久久久| 一区二区日韩欧美中文字幕| 精品一品国产午夜福利视频| 国产精品国产av在线观看| 亚洲精品美女久久av网站| 免费看不卡的av| 精品福利观看| av国产精品久久久久影院| 亚洲av电影在线观看一区二区三区| 日本av手机在线免费观看| 真人做人爱边吃奶动态| 女人精品久久久久毛片| 婷婷色综合大香蕉| 51午夜福利影视在线观看| 你懂的网址亚洲精品在线观看| 99国产综合亚洲精品| 1024视频免费在线观看| 成人黄色视频免费在线看| 欧美黄色淫秽网站| 日韩一卡2卡3卡4卡2021年| 性少妇av在线| 久久热在线av| 亚洲欧美日韩高清在线视频 | 伊人久久大香线蕉亚洲五| 男人添女人高潮全过程视频| 宅男免费午夜| 国产一卡二卡三卡精品| 韩国高清视频一区二区三区| 亚洲精品国产区一区二| 色精品久久人妻99蜜桃| 婷婷色麻豆天堂久久| 只有这里有精品99| 亚洲一码二码三码区别大吗| 亚洲av电影在线进入| 欧美黄色片欧美黄色片| 十八禁网站网址无遮挡| 成人亚洲精品一区在线观看| 国产免费福利视频在线观看| 色网站视频免费| www.av在线官网国产| 19禁男女啪啪无遮挡网站| 亚洲国产日韩一区二区| 捣出白浆h1v1| 两个人免费观看高清视频| 国产日韩欧美亚洲二区| 无遮挡黄片免费观看| 久久久精品免费免费高清| 尾随美女入室| 日韩精品免费视频一区二区三区| 精品国产一区二区三区四区第35| 久久毛片免费看一区二区三区| 性高湖久久久久久久久免费观看| 99国产精品一区二区蜜桃av | 最近中文字幕2019免费版| 色精品久久人妻99蜜桃| 国产极品粉嫩免费观看在线| 18禁裸乳无遮挡动漫免费视频| 一级毛片 在线播放| 亚洲av欧美aⅴ国产| 亚洲精品美女久久av网站| 国产成人精品在线电影| 久久精品成人免费网站| 国产极品粉嫩免费观看在线| 欧美久久黑人一区二区| 国产片内射在线| 亚洲免费av在线视频| 欧美成狂野欧美在线观看| 菩萨蛮人人尽说江南好唐韦庄| 成人三级做爰电影| 一级毛片女人18水好多 | av有码第一页| 一级,二级,三级黄色视频| 亚洲熟女精品中文字幕| 久久久久久亚洲精品国产蜜桃av| 国产精品一区二区免费欧美 | 日韩av免费高清视频| 99久久综合免费| 一边摸一边做爽爽视频免费| 国产亚洲av片在线观看秒播厂| 国产xxxxx性猛交| 嫩草影视91久久| 侵犯人妻中文字幕一二三四区| 亚洲午夜精品一区,二区,三区| 亚洲av在线观看美女高潮| 国产精品久久久人人做人人爽| 女性生殖器流出的白浆| 亚洲人成电影观看| 色婷婷av一区二区三区视频| 成年美女黄网站色视频大全免费| 大香蕉久久网| 国产成人精品久久久久久| 国产精品偷伦视频观看了| 亚洲色图 男人天堂 中文字幕| 欧美日韩成人在线一区二区| av福利片在线| 久久综合国产亚洲精品| 曰老女人黄片| 亚洲欧美精品综合一区二区三区| 曰老女人黄片| 亚洲成人免费电影在线观看 | 99国产综合亚洲精品| 久久久久久久久久久久大奶| 久久九九热精品免费| 亚洲人成电影免费在线| 成人国产av品久久久| 精品久久久精品久久久| av电影中文网址| 丝袜脚勾引网站| 蜜桃国产av成人99| 老熟女久久久| 亚洲精品久久成人aⅴ小说| 18禁国产床啪视频网站| 久久久国产精品麻豆| cao死你这个sao货| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲人成电影免费在线| 操美女的视频在线观看| 亚洲欧洲精品一区二区精品久久久| 国语对白做爰xxxⅹ性视频网站| 中文字幕制服av| 日韩制服丝袜自拍偷拍| 日本五十路高清| 午夜福利乱码中文字幕| 国产三级黄色录像| 日本色播在线视频| 精品一区在线观看国产| 好男人视频免费观看在线| 国产成人一区二区在线| 国产欧美日韩精品亚洲av| 欧美精品人与动牲交sv欧美| 少妇的丰满在线观看| 国产1区2区3区精品| 亚洲自偷自拍图片 自拍| 色播在线永久视频| 在线观看免费高清a一片| 亚洲av片天天在线观看| 啦啦啦啦在线视频资源| 国产精品久久久人人做人人爽| 国产福利在线免费观看视频| 久久久精品国产亚洲av高清涩受| www.精华液| 欧美日韩视频高清一区二区三区二| 一级毛片我不卡| 国产精品久久久久久精品古装| 成人午夜精彩视频在线观看| 亚洲九九香蕉| 黄色片一级片一级黄色片| 国产日韩一区二区三区精品不卡| 国产精品一区二区精品视频观看| 欧美人与性动交α欧美软件| 久久久久久久大尺度免费视频| 精品少妇黑人巨大在线播放| 国产福利在线免费观看视频| 99热网站在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲熟女毛片儿| 91麻豆精品激情在线观看国产 | 国产成人免费观看mmmm| 成人影院久久| 欧美日韩av久久| 免费少妇av软件| 亚洲第一青青草原| videosex国产| 嫁个100分男人电影在线观看 | 欧美+亚洲+日韩+国产| 91麻豆精品激情在线观看国产 | 中文字幕精品免费在线观看视频| 飞空精品影院首页| 91麻豆精品激情在线观看国产 | 日日摸夜夜添夜夜爱| 日韩免费高清中文字幕av| 亚洲精品在线美女| 精品一区二区三卡| 五月开心婷婷网| 亚洲精品自拍成人| 一本—道久久a久久精品蜜桃钙片| 精品第一国产精品| 国产精品一区二区精品视频观看| 乱人伦中国视频| 免费在线观看完整版高清| 亚洲av日韩在线播放| 国产国语露脸激情在线看| 国产精品一区二区在线不卡| 美国免费a级毛片| 老司机靠b影院| 亚洲人成网站在线观看播放| 欧美精品一区二区大全| 亚洲综合色网址| 老司机影院成人| 女人被躁到高潮嗷嗷叫费观| 国产一区有黄有色的免费视频| 亚洲精品久久成人aⅴ小说| 成年美女黄网站色视频大全免费| 亚洲精品乱久久久久久| 国产片特级美女逼逼视频| 国产免费又黄又爽又色| 黄色视频不卡| 男的添女的下面高潮视频| 亚洲成国产人片在线观看| 一本综合久久免费| 亚洲第一av免费看| 美女大奶头黄色视频| 久久ye,这里只有精品| 一边摸一边做爽爽视频免费| 日韩欧美一区视频在线观看| 在线看a的网站| 午夜激情久久久久久久| 欧美日韩精品网址| 国产91精品成人一区二区三区 | 亚洲色图综合在线观看| 精品福利观看| 亚洲第一av免费看| 国产亚洲精品久久久久5区| 欧美中文综合在线视频| 麻豆乱淫一区二区| 亚洲精品美女久久久久99蜜臀 | 午夜福利影视在线免费观看| 亚洲国产成人一精品久久久| 亚洲精品美女久久久久99蜜臀 | 久久久久国产精品人妻一区二区| 日韩视频在线欧美| 久久免费观看电影| 亚洲欧洲精品一区二区精品久久久| 热re99久久精品国产66热6| 成年人免费黄色播放视频| 日本一区二区免费在线视频| 亚洲中文日韩欧美视频| 久久久久久久国产电影| 色网站视频免费| 一本—道久久a久久精品蜜桃钙片| 老司机深夜福利视频在线观看 | 好男人电影高清在线观看| 成在线人永久免费视频| 色综合欧美亚洲国产小说| 亚洲欧美日韩另类电影网站| 国产成人一区二区在线| 欧美激情极品国产一区二区三区| 波多野结衣一区麻豆| 久久国产精品人妻蜜桃| 曰老女人黄片| 久久av网站| 久久久精品国产亚洲av高清涩受| 免费在线观看视频国产中文字幕亚洲 | 午夜免费男女啪啪视频观看| 男女免费视频国产| 大码成人一级视频| 一边摸一边抽搐一进一出视频| 美女午夜性视频免费| 人成视频在线观看免费观看| 国产亚洲午夜精品一区二区久久| 男人爽女人下面视频在线观看| 一区二区日韩欧美中文字幕| 亚洲黑人精品在线| 久久国产精品人妻蜜桃| 午夜视频精品福利| 国产欧美日韩精品亚洲av| av福利片在线| 国产免费视频播放在线视频| 天天躁夜夜躁狠狠躁躁| 久久精品人人爽人人爽视色| 久久精品久久久久久噜噜老黄| 亚洲 欧美一区二区三区| 婷婷成人精品国产| 久久鲁丝午夜福利片| 日韩伦理黄色片| 精品少妇内射三级| 精品人妻一区二区三区麻豆| 国产熟女欧美一区二区| 日韩中文字幕视频在线看片| www日本在线高清视频| 巨乳人妻的诱惑在线观看| 人成视频在线观看免费观看| 夫妻午夜视频| 一边摸一边做爽爽视频免费| 免费在线观看完整版高清| 后天国语完整版免费观看| 国产av国产精品国产| 老熟女久久久| 国产熟女欧美一区二区| 亚洲综合色网址| 不卡av一区二区三区| 另类精品久久| videosex国产| 国产男女内射视频| 亚洲精品在线美女| 欧美激情高清一区二区三区| 国产成人91sexporn| 99久久99久久久精品蜜桃| 涩涩av久久男人的天堂| 国产熟女欧美一区二区| 亚洲精品国产一区二区精华液| 日韩一区二区三区影片| 波野结衣二区三区在线| 国产精品欧美亚洲77777| 一二三四在线观看免费中文在| 亚洲精品美女久久av网站| 中文字幕色久视频| 一级黄片播放器| 亚洲av日韩精品久久久久久密 | 中文字幕精品免费在线观看视频| 精品一品国产午夜福利视频| 欧美精品一区二区免费开放| 国产精品 欧美亚洲| 国产精品久久久久成人av| 好男人电影高清在线观看| 亚洲黑人精品在线| 高清不卡的av网站| 黄色怎么调成土黄色| 99九九在线精品视频| 另类亚洲欧美激情| 国产一区二区在线观看av| 一边摸一边抽搐一进一出视频| 久久天堂一区二区三区四区| 男人操女人黄网站| 2021少妇久久久久久久久久久| 少妇被粗大的猛进出69影院| 亚洲精品久久久久久婷婷小说| 91精品国产国语对白视频| av在线老鸭窝| 国产精品一国产av| 2021少妇久久久久久久久久久| 亚洲成av片中文字幕在线观看| 少妇人妻久久综合中文| 水蜜桃什么品种好| 国产一区二区三区综合在线观看| 国产国语露脸激情在线看| 9热在线视频观看99| 午夜福利在线免费观看网站| av一本久久久久| 亚洲精品一区蜜桃| 亚洲精品久久久久久婷婷小说| 咕卡用的链子| 成年动漫av网址| 欧美日韩成人在线一区二区| 国产免费一区二区三区四区乱码| 一区福利在线观看| 欧美日韩福利视频一区二区| av不卡在线播放| 欧美黑人欧美精品刺激| 国产一区二区在线观看av| 亚洲精品第二区| 久久中文字幕一级| 亚洲欧美日韩高清在线视频 | 国产成人影院久久av| 中文字幕av电影在线播放| 午夜影院在线不卡| 久久九九热精品免费| 中文字幕高清在线视频| 日韩电影二区| 女警被强在线播放| 一二三四在线观看免费中文在| 美女脱内裤让男人舔精品视频| 久久影院123| 极品少妇高潮喷水抽搐| 国产亚洲精品久久久久5区| 欧美另类一区| 亚洲精品国产av蜜桃| 国产成人精品久久二区二区91| e午夜精品久久久久久久| 国产老妇伦熟女老妇高清| 看十八女毛片水多多多| 精品第一国产精品| 日韩视频在线欧美| 麻豆国产av国片精品| 久久久久网色| 女人久久www免费人成看片| 国产精品一区二区免费欧美 | 久久久国产欧美日韩av| 国产精品国产三级专区第一集| 欧美日韩国产mv在线观看视频| 欧美变态另类bdsm刘玥| 日本午夜av视频| 在现免费观看毛片| 91精品三级在线观看| 亚洲国产中文字幕在线视频| 精品亚洲成a人片在线观看| 国产熟女午夜一区二区三区| 国产黄色视频一区二区在线观看| 久久中文字幕一级| 50天的宝宝边吃奶边哭怎么回事| 黄色 视频免费看| 国产熟女欧美一区二区| 一区二区三区四区激情视频| 国产av一区二区精品久久| 首页视频小说图片口味搜索 | 精品国产乱码久久久久久男人| 精品一区二区三区av网在线观看 | 免费高清在线观看视频在线观看| 亚洲国产精品成人久久小说| 大香蕉久久网| 男女高潮啪啪啪动态图| 久久久国产欧美日韩av| 少妇人妻 视频| 老司机靠b影院| 电影成人av| 中文字幕高清在线视频| 18禁国产床啪视频网站| √禁漫天堂资源中文www| 在线观看人妻少妇| 少妇人妻久久综合中文| 久久久久久久大尺度免费视频| 国产精品成人在线| 国产黄色视频一区二区在线观看| 欧美激情高清一区二区三区| 欧美国产精品一级二级三级| 999精品在线视频| 免费看十八禁软件| 水蜜桃什么品种好| cao死你这个sao货| 欧美 亚洲 国产 日韩一| 丝袜美腿诱惑在线| 丝袜美足系列| 亚洲男人天堂网一区| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧美一区二区三区黑人| 国产免费一区二区三区四区乱码| 一二三四社区在线视频社区8| 亚洲精品国产色婷婷电影| 精品一区二区三卡| 国产成人免费观看mmmm| 亚洲国产欧美日韩在线播放| 黄频高清免费视频| 精品久久久精品久久久| 国产男女超爽视频在线观看| 啦啦啦中文免费视频观看日本| 成人午夜精彩视频在线观看| 成人国语在线视频| 老熟女久久久| av福利片在线| 激情视频va一区二区三区| 高清黄色对白视频在线免费看| netflix在线观看网站| 19禁男女啪啪无遮挡网站| 亚洲国产精品一区三区| 亚洲精品久久成人aⅴ小说| 各种免费的搞黄视频| 亚洲国产精品国产精品| 亚洲人成电影免费在线| 亚洲av综合色区一区| 最新的欧美精品一区二区| 国产精品久久久久久精品电影小说| 丝袜美腿诱惑在线| 亚洲中文日韩欧美视频| 九草在线视频观看| 亚洲欧美日韩高清在线视频 | 精品欧美一区二区三区在线| 999久久久国产精品视频| 久久天躁狠狠躁夜夜2o2o | 男女下面插进去视频免费观看| 国产精品熟女久久久久浪| 纵有疾风起免费观看全集完整版| 一本综合久久免费| 中文字幕人妻丝袜制服| 老司机深夜福利视频在线观看 | 亚洲人成77777在线视频| 久久久久久久精品精品| 人妻一区二区av| 久久午夜综合久久蜜桃| 久久久精品免费免费高清| 丰满人妻熟妇乱又伦精品不卡| av欧美777| 国产片内射在线| 宅男免费午夜| 高潮久久久久久久久久久不卡| 亚洲精品美女久久久久99蜜臀 | 国产精品国产av在线观看| 亚洲欧洲日产国产| 精品人妻在线不人妻| 夫妻午夜视频| 欧美激情 高清一区二区三区| 9色porny在线观看| 亚洲精品美女久久久久99蜜臀 | 亚洲第一青青草原| 黄片播放在线免费| 亚洲久久久国产精品| 一级毛片女人18水好多 | 精品一区在线观看国产| 人妻 亚洲 视频| 黑人猛操日本美女一级片| 国产又色又爽无遮挡免| 丰满迷人的少妇在线观看| 久久精品国产亚洲av涩爱| av在线播放精品| 午夜福利免费观看在线| 99久久99久久久精品蜜桃| 丝袜美腿诱惑在线| 亚洲精品日本国产第一区| 天天躁夜夜躁狠狠久久av| 欧美日韩精品网址| 亚洲一区二区三区欧美精品| 国产精品一二三区在线看| 一区在线观看完整版| 日韩制服丝袜自拍偷拍| 成人国语在线视频| 国产真人三级小视频在线观看| 中文字幕最新亚洲高清|