Comparative analysis of genetic diversity in Tunisian apricot germplasm using AFLP and SSR markers

Eighty-one accessions representing apricot germplasm in Tunisia were collected from different areas of cultivation and fingerprinted using amplified fragment length polymorphism (AFLP) and microsatellites (SSR) markers. A total of 339 polymorphic markers were revealed using 5 AFLP primers combinations and 24 SSR loci. AFLP and SSR markers expressed a high level of polymorphism allowing the distinction of the accessions with an efficiency coefficient of discrimination of 100% for AFLP and 97% for SSR markers. Genetic diversity structure was assessed with AFLPs and SSRs markers separately then with combined matrix data by the help of hierarchical clustering elaborated using Wards method based on Nei and Li (1979) distances. Comparison of the obtained dendrograms revealed a phylogeographic structure into two major groups with significant conservation between the observed subgroups in relation with the geographic origin of the accessions. The relative efficiency of the markers in determining the genetic relationships among apricot accessions has been assessed and a combination of AFLPs and SSRs markers was the most effective. In addition, Mantel test based on genetic distances indicated highly significant correlation between AFLP-SSR data and each of the AFLP and SSR ones, with Pearson correlation values of r = 0.873 and r = 0.692, respectively, revealing the higher efficiency of the combination of both molecular techniques (AFLP and SSR) to estimate the levels of genetic variability among apricot germplasm.     

Genetic diversity and population genetic structure of pomegranate (Punica granatum L.) in Iran using AFLP markers

Pomegranate (Punica granatum L.) is one of the oldest known edible fruits. It is native to Iran and spread from Iran to other areas. In this study amplified fragment length polymorphism (AFLP) was used to detect intra- and inter-population genetic diversity of pomegranate. A group of 67 accessions belonged to 4 populations from Iran was studied using eight primer combinations. A total of 221 scorable bands were amplified, of which, 118 (54.13%) were polymorphic. Resolving power (Rp) ranged from 5.70 to 9.21, and the average of polymorphism information content (PIC) per primer pair was 0.40. According to Nei's gene diversity and allelic statistics, Isfahan population had a highest genetic diversity (H = 0.3646, I = 0.5327, N_e = 1.6467). Coefficient of gene differentiation between populations (GST) was 0.124, indicated that mainly proportion of genetic variation (87.6%), was within populations and the remaining (12.4%) of the variation was among populations that, also supported by analysis of molecular variance (AMOVA). The gene flow (Nm) varied from 0.969 to 10.404 between pair-wise populations and was 3.504 among all of the populations. The Jaccard similarity coefficient between individuals ranged from 0.26 to 0.88. The UPGMA dendrogram clustered all 67 accessions into 6 groups. In some cases accessions from same region were grouped together but in most cases, there was gene exchange. To study the genetic relationships among populations, a principal coordinate analysis (PCoA) based on Nei's genetic distances was performed. Results of this study showed that AFLP marker can be a useful tool for investigating the genetic diversity of pomegranate genotypes.     

Genetic relationships in Eriobotrya species as revealed by amplified fragment length polymorphism (AFLP) markers

Although most of the species in Eriobotrya originated in China there are few studies on the genetics and genetic relationship of this genus. The aim of the present study was to clarify genetic relationships among 18 Eriobotrya accessions from China using amplified fragment length polymorphism (AFLP). Two close relatives Raphiolepisindica (L.) Lindl. and Photiniaserrulata Lindl. were used as outgroups. The results showed that 5 selected AFLP primer combinations amplified 297 scorable bands, of which 282 were polymorphic, yielding a polymorphism rate of 95%. Genetic similarity coefficient among Eriobotrya accessions ranged from 0.82 (between E. bengalensis Hook.f. and E. bengalensis forma angustifolia Vidal) to 0.41 (between E. serrate Vidal and E. deflexa Nakai) indicating substantial genetic variations in the Eriobotrya accessions evaluated. The unweighted pair-group method of arithmetic averages (UPGMA) cluster analysis, based on genetic similarity coefficient, produced a dendrogram which grouped 20 accessions into three main clusters. The result of cluster analysis was generally in agreement with the known taxonomic grouping. Furthermore, the clustering of the 18 accessions of Eriobotrya was consistent with systematic classification based on morphological characteristics.     

Variety identification and comparative analysis of genetic diversity in yardlong bean (Vigna unguiculata spp. sesquipedalis) using morphological characters,SSR and ISSR analysis

Genetic diversity and relatedness of 23 yardlong bean (Vigna unguiculata spp. sesquipedalis) accessions and 7 accessions of a hybrid between cowpea (V. unguiculata spp. unguiculata) and yardlong bean (dwarf yardlong bean) in Thailand were estimated using morphological characters, simple sequence repeat (SSR) and inter-simple sequence repeat (ISSR) markers. In addition, two mungbean (Vignaradiata (L.) Wilczek) and two blackgram (Vigna mungo (L.) Hepper) accessions were also used as outgroup species for molecular analysis. Five morphological characters were diverse among most accessions. However, five groups of 2–3 accessions could not be distinguished from one another based on these morphological characters alone. Unweighted pair-group arithmetic average (UPGMA) analysis of these characters separated these 30 accessions into 2 major groups; the yardlong bean group and the dwarf yardlong bean group. Eleven of the sixteen SSR primers yielded clear SSRs, ten of which were polymorphic (90.91% polymorphism), detecting a total of 54 alleles with an average of 4.91 alleles per locus. These 10 polymorphic SSR markers successfully distinguished 28 yardlong bean and dwarf yardlong bean accessions. The polymorphic information content (PIC) among genotypes varied from 0.251 to 0.752 with an average of 0.597. Among the 16 ISSR primers used, a total of 312 ISSR fragments were amplified for these three Vigna species, revealing the polymorphism percentage of 91.03%. The average ISSR PIC value (0.197) with the range of 0.137–0.276 was lower than that of SSR. Nevertheless, the average marker index of this multilocus marker was 3.495, which was higher than that of SSR (0.669), owing to the differences in the effective multiplex ratio. In addition, Mantel test cophenetic correlation coefficient was higher for ISSR (0.566) than that of SSR (0.198). These results indicated higher efficiency of ISSR for estimating the levels of genetic diversity and relationships among yardlong beans and dwarf yardlong beans in this study. Pair-wise coefficients of SSR- and ISSR-based genetic similarity among all yardlong bean and dwarf yardlong bean accessions averaged 0.87 and 0.91, respectively, suggesting a narrow genetic base that emphasizes the need to broaden genetic diversity to ensure continued breeding success. Clustering of genotypes within groups was not similar when SSR and ISSR derived dendrograms from UPGMA analysis were compared. It appeared that ISSR was the most effective marker system in determining the genetic variability and relationships among yardlong bean and dwarf yardlong bean accessions and differentiating three Vigna species. In addition, ISSR was also most useful for variety identification since all 30 yardlong beans and dwarf yardlong bean accessions can be effectively distinguished by only four ISSR primers with the highest PIC values.     

Genetic relatedness (diversity) and cultivar identification by randomly amplified polymorphic DNA (RAPD) markers in teasle gourd (Momordica dioica Roxb.)

Genetic relationship and variation of 29 accessions of teasle gourd (Momordica dioica Roxb.) and 1 accession of Momordica cochinchinensis Spreng. (wild relatives of teasle gourd) were examined by RAPD analysis using 44 dodecamer oligonucleotide primers. A total of 496 fragments were produced by 44 primers of which 95% bands were polymorphic. Using presence or absence of specific RAPD markers or combination of primers, 23 out of 30 accessions were identified. The genetic relatedness or genetic distance based on Nei and Li's genetic similarity varied from 0.86 to 0.65 with an average of 0.74 among 29 M. dioica accessions (when M. cochinchinensis excluded). In the phenetic dendrogram developed from cluster analysis using UPGMA method, M. cochinchinensis was out grouped as single accession, while others showing relatively weak grouping formed four groups. Clustering pattern did not demonstrate any relationship between geographical origin and genetic diversity. A DNA extraction method has been standardized. This is the first report of using RAPD techniques in teasle gourd. It was concluded that RAPD analysis is a useful tool for genotypic identification and estimation of genetic similarity in teasle gourd.     

ISSR polymorphism in Indian wild orange (Citrus indica Tanaka,Rutaceae) and related wild species in North-east India

Inter simple sequence repeats (ISSR) polymorphism in Citrus indica Tanaka (Rutaceae), an endemic and threatened wild species, was examined along with three other closely related wild taxa (C. medica L., C. latipes (Swingle) Tanaka, and C. sp. ‘Memang athur’) by analyzing 53 representative accessions sampled from North-east India. Jaccard's similarity values among 53 accessions of Citrus ranged from 0.46 to 0.97 (average = 0.75). Genetic similarity values among all the 34 accessions of C. indica were found in the range of 0.82 to 0.97 with an average of 0.90. Heterozygosity (Ht = 0.123) and Shannon's information index (I = 0.188) values estimated for C. indica revealed significantly low level of genetic variation within the species. UPGMA dendrogram grouped all 53 accessions of Citrus into four major clusters: Cluster I – C. latipes; Cluster II – C. medica; Cluster III – ‘Memang athur’ and Cluster IV – C. indica. The dendrogram placed all the 34 accessions of C. indica in five sub-clusters under Cluster IV. The placement of C. indica accessions in various sub-clusters and groups in the dendrogram was based on molecular differentiation of individual accessions rather than their geographical origin. Very low genetic diversity and destruction of its natural habitat pose serious threat to C. indica even in the Citrus Gene Sanctuary in Nokrek Biosphere Reserve (NBR) in Meghalaya. Low genetic variability, heterozygosity and Shannon's information index in C. medica, C. latipes and ‘Memang athur’ are also concerns that need to be addressed for developing appropriate strategies to conserve the genetic diversity extant in these valuable genetic resources.     

中国马铃薯部分栽培品种遗传多样性的AFLP分析

 采用AFLP方法分析中国各地79个马铃薯栽培品种, 获得了8对引物组合的扩增谱带, 为马铃薯品种的鉴定提供了分子依据。8对引物组合扩增结果共获得801条带, 平均每对引物组合产生100.13条带。其中493条为多态性条带, 平均多态性检出率为61.2% , 引物组合E11 /M51扩增的多态性比率最高。聚类分析将材料分成3类, 聚类结果与地域无明显相关。     

甘肃中部梨种质资源的AFLP分析

采用AFLP分子标记技术对甘肃中部梨种质资源的遗传多样性进行分析。结果表明,6对AFLP引物在40份甘肃梨种质中共扩增出472条带,其中多态性带为404条,多态率高达86%,显示了甘肃中部梨资源丰富的遗传多样性。任何一对引物可以鉴别所有40份资源,显示了很高的鉴别能力。采用UPGMA聚类分析法构建的系统树在相似系数为0.72时将40份种质分为7个组:西洋梨、新疆梨、白梨(砂梨)、木梨、褐梨组和2个秋子梨组。所有白梨品种与唯一的砂梨品种黄花梨聚为一个大组。形态学上归属不明的品种分别聚类到西洋梨、白梨、木梨和秋子梨组中。研究表明基于AFLP标记的梨资源分类体系可以反映甘肃地方梨品种和类型间的遗传多样性和亲缘关系。     

Genetic diversity and population structure of an endangered Orchid (Dendrobium loddigesii Rolfe) from China revealed by SRAP markers

Dendrobium loddigesii Rolfe is an endangered perennial herb with ornamental and medicinal value. Due to habitat deterioration and human over-exploitation, it has suffered a significant decline in abundance. Determining the level of genetic diversity and pattern of population genetic structure of this species would be helpful for its conservation and management. In this paper, sequence-related amplified polymorphism (SRAP) markers were applied to assess the level and pattern of genetic diversity in seven populations of D. loddigesii. Seventeen SRAP primer combinations generated a total of 231 clear amplification bands encompassing 187 (80.95%) polymorphic bands. A high level of genetic diversity was detected (PPB = 80.52%, H = 0.2743, I = 0.4113) at the species level. There was a moderate genetic differentiation (G_st = 0.304) among populations. Two main clusters were detected by cluster analysis using the unweighted pair-group method with arithmetic average (UPGMA). Mantel test revealed that no significant positive correlation was found between genetic distances and geographic distances (r = 0.2302; P > 0.05). Recommendations for conservation of the endangered species resources are proposed.     

Genetic diversity analysis in the traditional and improved ginger (Zingiber officinale Rosc.) clones cultivated in North-East India

Identification of clonal or genotypic variations is a prerequisite for ginger (Zingiber officinale Rosc.) improvement programmes. Genetic diversity analysis was carried out in a set of forty-nine ginger clones cultivated in North-East India using random amplified polymorphic DNA (RAPD) markers. The set included clones of released varieties and clones collected from various parts of North East India. Jaccard's genetic similarity, cluster analysis and principal component analysis identified five clusters. Cluster V included four clones traditionally cultivated in the Indian state of Meghalaya known for production of high-quality ginger indicating that the clones were a good candidate for ginger improvement. Specific bands for these clones were also identified. Principal component analysis of the molecular data supported grouping of the clones into six hypothetical populations based on their source or location of collection.     

Fingerprinting of three typical macrosperma Italian lentil (Lens culinaris Medik.) landraces using fluorescence-based AFLP markers

Italian lentil landraces are principally cultivated for self or local consumption. Most of them are disappearing, particularly macrosperma types by being less required by the market. A pre-requisite for the conservation and the efficient use of genetic resources is the better understanding of the extent and the distribution of the existing genetic variation, useful for future breeding programmes. Our study was undertaken to analyse and quantify the genetic diversity within and among three macrosperma Italian lentil landraces (Onano, Altamura and Villalba), using fluorescent AFLP markers. AFLP markers generated information to differentiate among closely related genotypes and group within the same cluster individuals belonging to the same landrace. The total genetic diversity (H_T), the genetic diversity within population (H_S) and the extent of differentiation between populations (D_ST) were 0.198, 0.155 and 0.043, respectively. The fixation index (G_ST = 0.219) showed that about 78% of the observed total genetic variation can be attributed to within population differences and around 22% is due to differences among populations. The gene flow estimate (N_m = 1.774) and the mean genetic distance value (0.077) suggested narrow genetic base among the analysed populations, confirming the tendency of Italian lentil landraces to group together. The present study showed that fluorescence-based AFLP technique is a biotechnological tool that can provide significant insights for research in genetic diversity of lentil landraces and their subsequent conservation and utilization in breeding programs.     

Using bulked AFLP analysis to assess genetic diversity in Echinacea species

Echinacea is an allogamous genus, thus its cultivars or populations are genetically heterogeneous. Using amplified fragment length polymorphism (AFLP) to estimate the genetic diversity of Echinacea is generally limited by the large number of individual plants and the higher cost that need to be processed. In the present study, effectiveness of several sizes of DNA bulking (10, 15, 20, 25 and 30 individuals) with 20, 36 and 55 primer pairs was compared using AFLP in determining the genetic diversity of Echinacea species. The results indicated that the use of bulked DNA-based AFLP analysis by using the selected eight primer pairs was capable of detecting genetic diversity between the tested Echinacea species, provided that the potential presence of low frequency variants was ignored and a possible bias in the estimates of genetic similarity was accepted. The assessments showed that a bulk of 15 individuals could detect sufficient AFLP variations at most genomic sites. Additionally, 20 primer pairs could generate sufficient polymorphic fragments to achieve high resolving power of AFLP for the tested Echinacea species.     

Assessment of diversity in Podophyllum hexandrum by genetic and phytochemical markers

For successful conservation and domestication of a species, evaluation of its genetic diversity by different markers is important. Morphological characteristics, phytochemical variation and random amplified polymorphic DNA (RAPD) profiles were generated in different accessions of Podophyllum hexandrum in order to determine the genetic diversity. Random amplified polymorphic DNA (RAPD) analysis revealed a high degree of genetic diversity among the accessions used in the study. There was also high diversity in the concentration of marker compounds in the collected samples as revealed by HPLC analysis. It is shown that the approaches used in the work successfully discriminate between the accessions of this species and thus they constitute interesting tools to analyze molecular, biochemical and phenotypic diversity within this species. Similarity measurement using UPGMA followed by cluster analysis resulted in formation of many groups based on geographical distribution that generally reflected expected trends between the genotypes. There were also some important exceptions like PW-S, an accession from Wastoorwan, Khrew showing close resemblance to PG-S and PG-B collected from Gulmarg but grown at two different gene banks at Srinagar and Bonera. Further an accession PSH-B from Keller was significantly diverse from the rest of the native genotypes phytochemically, morphologically and at molecular level. RAPD data analysis was found to be significant predictor of phytochemical markers in cultivated P. hexandrum germplasm. Twelve accessions grown in gene bank repository were subjected to RAPD analysis and were assessed for content of podophyllotoxin and podophyllotoxin β-d-glycoside by HPLC. Individual regressions of podophyllotoxin and podophyllotoxin β-d-glycoside by RAPD analysis against HPLC has been found to determine linear values. Strong correlation and a strong association of values of the phytochemical variables and the DNA polymorphism data has been recorded.     

Analysis of genetic diversity in Portuguese Ceratonia siliqua L. cultivars using RAPD and AFLP markers

Although carob (Ceratonia siliqua L.) is of great economic importance little is still known about the pattern of genetic variation within this species. Morphological characteristics based on 31 fruit and seeds of continuous characters determinant for agro-industrial uses, were compared with RAPD and AFLP markers for assessing genetic distances in 68 accessions of carob trees, from different cultivars, varieties and eco-geographic regions of Algarve. Eighteen selected RAPD primers applied to the 68 accessions produced a total of 235 fragments ranging from 200 to 2000 bp, of which 93 (40%) were polymorphic. Four AFLP selective primer combinations generated a total of 346 amplification fragments of which 110 were polymorphic. The average level of polymorphism based on four primer combinations was 31.8%. The phenetic trees based on RAPD and AFLP analyses gave high co-phenetic correlation values, and were found to be consistent in general with the analysis of morphological data, carried out on the same accessions. A number of RAPD and AFLP markers were found to be diagnostic for ‘Canela’ cultivar and 13 wild ungrafted trees.     

基于SCoT 标记的芥菜种质遗传多样性与指纹图谱

利用SCoT 标记对46 份芥菜种质资源进行遗传多样性分析。从70 条引物中筛选出21 条用于PCR 扩增,共扩增出200 条条带,其中多态性条带140 条,多态性比率为70%。供试材料相似系数介于0.77～0.95 之间,平均Nei’s 基因多样性为0.198,平均Shannon 信息指数为0.301,平均多态信息含量为0.822 7。利用UPGMA 构建46 份芥菜种质资源的聚类树状图,在相似系数D1=0.780 处,可以将46 份芥菜种质分为两大类;而在D2=0.807 处,第Ⅱ类可进一步细分为5 个亚类。选出清晰的16 个多态性位点,构建46 份芥菜种质的SCoT 指纹图谱,结果由引物SP4、SP19、SP20、SP23、SP30 所构建的DNA 指纹图谱可将供试材料完全区分开。     

Characterization of Tunisian pomegranate (Punica granatum L.) cultivars using amplified fragment length polymorphism analysis

The amplified fragment length polymorphism (AFLP) analysis of DNA was used to characterize 34 pomegranate cultivars. By using a combination of six primers, a total of 327 markers were scored with a mean of 57.5. The high percentage of polymorphic bands (ppb) of 94.7 and the resolving power (R_p) collective rate value of 129.14 were scored. Data proved that the tested primers were informative to discriminate among cultivars and to survey the genetic diversity in this fruit crop. It has been assumed that the local pomegranate germplasm is characterized by a typically continuous genetic diversity. The derived dendrogram proved that cultivars are clustered independently from their geographical origin and their denomination. In addition, AFLP permitted the generation of a nearly unlimited number of molecular markers that are reliable in differentiating the cultivars and/or the polyclonal varieties.     

Molecular diversity of Chinese Cucurbita moschata germplasm collections detected by AFLP markers

Cucurbita moschata is an important vegetable crop. Although a total of 1032 landraces of C. moschata are maintained in China, little is known about their genetic diversity. Molecular characterization is needed to facilitate the use of this C. moschata germplasm collection in breeding. Seventy-four Chinese accessions and 15 accessions from other countries were selected for evaluation based upon variation in fruit traits and geographical origin of molecular diversity with AFLP analysis. Nine pairs of EcoRI/MseI primers produced 500 fragments, of which 75.57% were polymorphic, indicating a high degree of diversity. The accessions from China were classified into two clusters, which were clearly differentiated from the accessions originating from Mexico, Guatemala, Honduras, and Ecuador. Chinese group genetically more closely related to other Asian countries group (India and Japan). In general, the accessions from the Americas had a greater number of unique loci than those from China. The differences are probably due to a limited number of introductions and genetic drift. The Americas are the center of origin of C. moschata and therefore more diverse. With AFLP analysis, the accessions did not clearly group according to fruit shape; however, sub-clusters exist in acorn- and dumbbell-shaped accessions. The assessment of genetic distance, along with some unique traits among the different genotypes, could be useful in further genetic studies and the selection of the most adequate accessions for use in breeding programs.     

An assessment of genetic variability and relationships within Asian pears based on AFLP (amplified fragment length polymorphism) markers

A total of 100 Pyrus L. accessions native mainly to East Asia were subjected to amplified fragment length polymorphism (AFLP) analysis to evaluate genetic variability and relationships among the accessions. Six AFLP primer combinations produced a total of 459 fragments, of which 410 were polymorphic with a polymorphism percentage of 89%. The Dice's similarity coefficient among pear accessions ranged from 0.671 (P. betulaefolia Bge and P. elaeagrifolia Pall.) to 0.947 (‘Umajirou’ and ‘Immuraaki’). Occidental pears generally had low similarities to Asian pears. The dendrogram generated from all the accessions by unweighted pair-group method of arithmetic analysis (UPGMA) cluster analysis clearly distinguished Occidental pears from accessions of East Asia. P. ussuriensis Maxim., P. betulaefolia and P. communis L. clustered separately into independent groups in accordance with their morphological classification. Japanese pear cultivars formed two groups with some Chinese white pears and Chinese sand pears. Chinese white pears and Chinese sand pears independently formed their own groups and also mingled into mixed groups in the dendrogram. Therefore, Chinese white pears were treated as a cultivated group or an ecotype of P. pyrifolia: P. pyrifolia White Pear Group. The information obtained from this study will be of great help for understanding the origin and evolution of Asian pear cultivars.     

黄瓜种质资源遗传多样性及其亲缘关系的AFLP分析

 采用AFLP分子标记技术，对中国黄瓜种质资源遗传多样性及其与外来种质的关系进行了分析， 结果表明8对AFLP引物在70份黄瓜种质中共扩增出425条带，多态性带的比例为66％。供试黄瓜种质的平 均期望杂合度为0．376，中国种质的平均期望杂合度为0．387，明显高于国外种质的n 291。西双版纳黄瓜和印度野生黄瓜具有一些栽培种质没有的特异位点，中国栽培种质的特异位点多于外来栽培种质，后者也有一些中国栽培种质没有的特异位点。聚类分析将70份种质分为三大组群，即西双版纳黄瓜(Cucumis sativus L.vaF．xishuangbannanesis Qi et Yuan)组群，e sativus var．hardwickii野生黄瓜组群和栽培黄瓜组群。西双版纳黄瓜与栽培黄瓜的距离最远，与野生黄瓜次之。按一定的遗传距离可以将中国和外来栽培种质分开。大多数 华南型和华北型种质归属于不同的亚组。这些结果有助于有目的地利用这些变异拓宽育种材料的遗传背景。     

Use of ISSR markers to assess genetic diversity of African edible seeded Citrullus lanatus landraces

We used the twenty primers to evaluate the genetic variability of 80 individuals belonging to four accessions of edible seeded Citrullus lanatus originated from Côte d’Ivoire. Edible seeded C. lanatus, named “egusi” or “pistachio”, had a great importance in nutrition in West Africa. Nevertheless, due to its neglected status no study to our knowledge has been devoted to its genetic variability using DNA markers. The twenty ISSR primers generated 258 bands among which 252 were polymorphic (97.67%). On the whole, the bands generated revealed three types of profile sharply distinct from each other with minor differences within each type. One profile (P1) was most frequent with 65 individuals. Three accessions (NI084, NI127 and NI145) generated the three types of profile and had medium values of genetic diversity (GD = 0.246–0.275, respectively). On the opposite, the accession NI076 only contained individuals of the most represented type of profile (P1) and had the lowest genetic diversity (GD = 0.055 ± 0.017). The pairwise genetic distance between the 80 individuals varied from 0 to 0.61. The Factorial Component Analysis and the dendrogram clearly separated the 80 individuals into three clusters corresponding to the three types of profile. The results showed that clusters were well separated from each other whereas accessions were not. Our results suggest that high number of individuals should be taken into account for sampling missions and conservation strategies because accessions were not well differentiated from each other. Local agricultural practices consisting of frequent seeds exchanges between farmers and the conservation of harvested seeds for next year culture could be one explanation.