Tracking crop varieties using genotyping-by-sequencing markers: a case study using cassava (Manihot esculenta Crantz)
BMC Genetics volume 16, Article number: 115 (2015)
Accurate identification of crop cultivars is crucial in assessing the impact of crop improvement research outputs. Two commonly used identification approaches, elicitation of variety names from farmer interviews and morphological plant descriptors, have inherent uncertainty levels. Genotyping-by-sequencing (GBS) was used in a case study as an alternative method to track released varieties in farmers’ fields, using cassava, a clonally propagated root crop widely grown in the tropics, and often disseminated through extension services and informal seed systems. A total of 917 accessions collected from 495 farming households across Ghana were genotyped at 56,489 SNP loci along with a “reference library” of 64 accessions of released varieties and popular landraces.
Accurate cultivar identification and ancestry estimation was accomplished through two complementary clustering methods: (i) distance-based hierarchical clustering; and (ii) model-based maximum likelihood admixture analysis. Subsequently, 30 % of the identified accessions from farmers’ fields were matched to specific released varieties represented in the reference library. ADMIXTURE analysis revealed that the optimum number of major varieties was 11 and matched the hierarchical clustering results. The majority of the accessions (69 %) belonged purely to one of the 11 groups, while the remaining accessions showed two or more ancestries. Further analysis using subsets of SNP markers reproduced results obtained from the full-set of markers, suggesting that GBS can be done at higher DNA multiplexing, thereby reducing the costs of variety fingerprinting. A large proportion of discrepancy between genetically unique cultivars as identified by markers and variety names as elicited from farmers were observed. Clustering results from ADMIXTURE analysis was validated using the assumption-free Discriminant Analysis of Principal Components (DAPC) method.
We show that genome-wide SNP markers from increasingly affordable GBS methods coupled with complementary cluster analysis is a powerful tool for fine-scale population structure analysis and variety identification. Moreover, the ancestry estimation provides a framework for quantifying the contribution of exotic germplasm or older improved varieties to the genetic background of contemporary improved cultivars.
Agricultural productivity in developing countries is affected by limited access to improved varieties, in addition to biotic, abiotic constraints and sub-optimal agronomic practices [1, 2]. Successful dissemination and adoption of improved varieties from both private and public breeding programs is expected to contribute positively to farm-level productivity and income generation. It is the role of household level impact assessment studies, particularly collection of variety specific adoption data, to determine whether this is happening [3, 4].
Traditionally, estimation of improved variety adoption in socio-economic impact studies relies mostly on: expert opinion of breeders, extension services and other experts; elicited responses from farmers in farmer-level surveys; and morphological descriptors. However, such methods have several inherent uncertainty levels. For example, variety naming systems in the absence of formal seed systems can be quite temporally and spatially variable leading to inconsistencies in the names of a particular variety. Also, environmental conditions and different stages of plant development influence morphological descriptors [5, 6]. Finally, the number of descriptors can be quite limited as varieties are developed to conform to desired ideotypes, thus greatly reducing the power to distinguish consanguineous varieties .
These challenges can be overcome by using molecular markers which are not only unaffected by the environmental factors and crop developmental stages but are also ubiquitous throughout plant genomes. Genome-wide markers, like single nucleotide polymorphisms (SNP), not only facilitate germplasm classification using genetic distance estimates but can also be used to quantify the relative proportion of ancestries derived from various founder genotypes of currently grown cultivars . Such inferences of ancestries are useful in understanding and/or reconstructing the evolution of successful varieties, either landraces or products of formal breeding programs that lack breeding pedigree records or where the varieties are derived from open-pollinated breeding methods . In the context of impact assessment of a specific breeding program, ancestry inferences can be useful in estimating the benefits resulting from the usage of its improved germplasm by other programs . This is because improved germplasm often moves easily throughout the network of plant breeding systems, resulting in research spill-over benefits.
In the past, simple sequence repeats and anonymous markers such as amplified fragment length polymorphisms and randomly amplified DNA polymorphisms have been used in DNA-based fingerprinting applications . However, due to inadequacies of these markers, including limited multiplexing ability, high genotyping costs and low frequency in the genome, they are increasingly being displaced by SNP markers generated from next-generation sequencing using reduced representation library (RRLs) methods. These recent methods rely on restriction enzymes to target a specific and reproducible subset of the genome for sequencing, thus allowing for simultaneous discovery and scoring of large numbers of markers. Genotyping-by-sequencing  is an RRL method that is relatively simple and inexpensive, making it feasible to genotype large populations of individuals. GBS has therefore become very popular, particularly for researchers working on non-model species with limited genomic resources .
Here, we report the use of GBS markers for cultivar identification with the objective of tracking released varieties in farmers’ fields, using cassava (Manihot esculenta Crantz) as a case study. Cassava is a highly heterozygous, clonally propagated species that originates from Latin America . Its starchy storage roots are the main source of calories for over 500 million people in the tropics . Africa is currently the leading producer of the crop accounting for more than 50 % of global production . Its ability to produce reasonable yields in marginal environments, its tolerance to drought and poor soils, and its ability for in-ground storage to allow piece-meal harvesting makes cassava one of the most important food-security crops in the continent . Despite its importance, planting materials are predominantly sourced from the informal seed system often from the farmer’s own harvest or exchange between farmers . Dissemination of new varieties has often been limited to efforts by the extension services connected to national programs and informal diffusion through farmer-to-farmer exchanges. This situation contributes to the challenge of tracking the spread of such varieties.
This study was conducted in three regions of Ghana covering the largest cassava producing area accounting for 61 % of cassava production in the country in 2010 [16, 19]. The three study regions included Brong Ahafo, Ashanti and Eastern (Fig. 1). A total of 495 households were selected using a multi-stage cluster sampling method. These households were distributed across 100 villages from 20 districts in the three study regions. For each of the 495 households surveyed, field sample collection entailed visiting one cassava field for each household with the largest number of varieties. A consent statement was read to the main decision maker of the household to inform him/her about the purpose of the study and to seek his/her permission to visit the cassava field to collect the leaf samples. Data and sample collection proceeded only if the farmer gave the consent. The GPS coordinates of the field were taken and farmers were asked to identify plants representing each of the varieties grown. Apical leaf samples were collected from one plant representing each variety and preserved in silica gel for transportation to a central laboratory at IITA in Ibadan, Nigeria for DNA extraction. Since a major objective of the present study was to identify specific cultivars in farmers’ fields, a “reference library” consisting of 64 clones representing released varieties (n = 16) and key landraces (n = 48) maintained by the Council for Science and Industrial Research-Crops Research Institute (CSIR-CRI) of Ghana were genotyped alongside the accessions from farmers’ fields. It should be noted that many of the released varieties in Ghana are landraces with superior agronomic traits (resistance to cassava mosaic disease, high yield and dry-matter content) and culinary qualities (root friability after boiling). These landraces were officially released as varieties following multi-year and -locational testing (Prof. S. Kantanka and Prof. J. P. Tetteh, personal communication).
DNA extraction and GBS
DNA was isolated from 1045 genotypes, representing 917 accessions collected from farmer’s fields and a library of 64 clones genotyped in duplicate (Additional file 1: Table S1). The Dellaporta method  with modifications described in  was used for high throughput DNA extraction. For genotyping-by-sequencing library preparation, we chose the ApekI restriction enzyme (recognition site: G|CWCG) that produces less variable distributions of read depth and therefore a larger number of scorable SNPs in cassava . Eleven 96-plex GBS libraries were constructed as described in  and sequenced at the Institute of Genomic Diversity at Cornell University using the Illumina HiSeq2500. Raw read sequences were processed through cassava GBS production pipelines developed using TASSEL 5.0 initially generated with about 2500 cassava clones under the NextGen Cassava project (www.nextgencassava.org) ; http://www.maizegenetics.net/#!tassel/c17q9). Resulting hapmap files (SNPs) were filtered with minor allele frequency (MAF) of 0.001 and coverage of 10x. SNPs were further processed by removing those with MAF of less than 0.01 and loci with more than 40 % missing data. The remaining missing SNP data-points were then imputed using GLMNET .
Identification of the cassava varieties was performed using three complementary clustering approaches: (i) pairwise distance-based hierarchical clustering; (ii) model-based maximum likelihood estimation of individual ancestries from multi-locus SNP genotype datasets using ADMIXTURE ; and (iii) Discriminant Analysis of Principal Components (DAPC) .
In the first approach, a pairwise genetic distance (identity-by-state, IBS) matrix was calculated from 56,489 SNP markers in PLINK . A Ward’s minimum variance hierarchical cluster dendrogram was then built from the IBS matrix using the Analyses of Phylogenetics and Evolution (ape) package  implemented in R . The critical distance threshold to declare whether two genotypes are identical was empirically determined from the distribution of pairwise distances between duplicated DNAs from 64 samples. This “calibration principle” approach  was taken because of the possibility of SNP genotype errors resulting from miscalling some heterozygous SNPs with low sequencing read depth as homozygotes .
In the second approach, ADMIXTURE analysis using the same set of 56,489 SNP markers was used to identify ancestries of the sampled cassava accessions. The number of sub-populations, K, was varied from 2 to 18 (K, in this case are considered founders of the currently cultivated varieties in the study regions). The most appropriate K value was selected after considering (i) 10-fold cross-validations whereby the best K exhibits low cross-validation error compared to other K values  and (ii) good correspondence with the clustering pattern obtained by the hierarchical tree.
To develop smaller sets of ancestry informative markers (AIMs) for follow-up studies using lower density genotyping, further ADMIXTURE–based ancestry estimation was carried out using decreasing subsets of SNP markers. These were selected based on Weir and Cockerham  FST, a measure of differences in allele frequencies among the subpopulations detected by ADMIXTURE. For comparative purposes, equivalent numbers of markers were randomly selected, each twenty times with replacement. The objective here was to see how much we can reduce the number of markers while still obtaining cluster assignment results that is close to that obtained from the full set of markers. We used the ‘supervised’ ADMIXTURE method assuming K = 11 . Accuracies of ancestry estimates was determined through correlations between the subsets and the complete set of 56,489 markers.
The model-based clustering approach implemented in ADMIXTURE assumes linkage equilibrium among loci and Hardy-Weinberg equilibrium within ancestral populations . However, such assumptions may be violated in vegetatively propagated species like cassava due to presence of clonal duplicates in germplasm collections. To validate the clustering pattern obtained from ADMIXTURE and the hierarchical clustering algorithms, we carried out Discriminant Analysis of Principal Components (DAPC), an assumption-free multivariate clustering method  using the R package ‘adegenet’  in a two-step process. Firstly, the optimal number of clusters was inferred using k-means analysis  of PCA-transformed genome-wide SNP data. After varying possible number of clusters from 2 to 40, Bayesian Information Criterion (BIC) was used to assess the best supported model i.e. the number and nature of clusters. Secondly, DAPC  was carried out on the clusters identified through k-means using the first 70 principal components. Membership probabilities of each individual for the different groups, akin to the sub-population membership coefficients from ADMIXTURE was obtained from DAPC. The results of DAPC analysis was then compared with those achieved from ADMIXTURE.
Field surveys found that farmers cultivated between one and five different varieties of cassava in their fields, but majority of them (>80 %) grow only one or two varieties (Fig. 2). A large number of unique farmer-elicited variety names (180) were associated with the 917 accessions collected from the three study regions of Ghana. Most of these names occurred five or less instances in the survey. The two most frequent names were “Debor” and “Ankra”, each recorded 90 and 87 times, respectively.
An average genetic distance between repeat genotyping of the 64 accessions in the “reference library” was below 0.05 (Ward’s distance, Additional file 2: Figure S1). We therefore chose 0.05 as the distance threshold below which we can declare that two accessions represent the same clone. The residual distance between same DNA is most likely due to miss-calling of heterozygotes as homozygotes from low sequencing read-depth, as is typical in high-multiplexing, sequence-based genotyping methods .
Genetic relationships among the 1045 genotyped accessions is described using a hierarchical clustering dendrogram (Fig. 3a) while the estimated ancestries (Q) obtained from ADMIXTURE are presented as a barplot (Fig. 3b). Major as well as minor clusters of genetically identical genotypes with genetic distances below the empirically determined distance threshold are clearly discernible. The two most dominant varieties (Cluster I and II) belong to the same branch of the dendrogram, and are therefore likely to share some common farmer-preferred characteristics. According to farmers’ naming system, the first variety which is associated with the most commonly recorded cultivar names (‘Debor’, ‘Ankra’ and ‘Bankye Kokoo’), is quite popular because of its excellent culinary traits, two of which are mealiness after boiling and relatively sweet taste. The remaining clones form a second large cluster that further subdivides into about nine clusters.
After elucidation of these groupings through hierarchical clustering, we turned to the STRUCTURE-like analysis  using the ADMIXTURE program  to assign individuals proportionally to hypothetical founder populations. After varying the number of sub-populations (K) from 2 to 18, the most appropriate number was found to be K = 11, which produced the lowest 10-fold cross-validation error compared to other K values (Fig. 4). These groupings corresponded to the hierarchical clustering dendrogram: each of the major branches of the dendrogram formed a distinct ancestry group.
Groups of clones with predominant ancestry membership to one of each of the identified ADMIXTURE sub-populations (>90 %) were discernible with the exception of group IX that had small admixture from groups II and III (Fig. 3b). The results of ADMIXTURE-based clustering is strongly supported by the topology of the distance-based dendrogram, with most of these genotypes also having very low IBS distance within their respective clusters (Fig. 3a). A large number of individual cassava genotypes (n = 277) that share ancestry between two or more of the identified varieties were also detected. Of these, about 157 had at least 50 % of their ancestry coming from one of the eleven sub-populations, while 120 accessions have multiple ancestries (Table 1). Moreover, the proportions of ancestries in these varieties appear to be consistent with simple crossing or backcrossing to produce F1 hybrids or backcross hybrids that may have occurred either in farmer’s fields or in formal plant breeding programs. Other genotypes show more complex multi-parent ancestries.
Following clustering of accessions into groups of genetically identical clones, actual variety identities were determined by matching each accession to the samples in the CSIR-CRI library. The library contained a total of 64 accessions but based on genetic similarity, these were collapsed into 34 unique cultivars (Additional file 2: Figure S1) of which 16 are released varieties. Using this library, we successfully classified a total of 282 accessions from the farmers’ fields as released varieties, representing about 30 % of the sampled 917 accessions (Table 1). These accessions matched only 8 of the 16 released varieties in CSIR-CRI reference library. Of the identified varieties, the most common was “IFAD”, also known as “UCC” and found in 158 households. The next most common variety was “Nkabom” (n = 65), followed by less common varieties of “Afisiafi”, “Tek Bankye”, “Bankye Broni” and “Doku Duade”, which occurred in 17, 12, 21 and 4 households, respectively. Although “Nkabom” is a released variety in the CSIR-CRI library, it was found to correspond to a superior landrace from Nigeria (TMEB3), one of the first clones discovered to harbor dominant resistance to cassava mosaic disease . It is therefore likely to have been introduced to Ghana through formal germplasm exchange between public breeding programs of the two countries. The least common of the released varieties was “Sika Bankye” found in only two of the surveyed households.
Besides the released varieties, a total of 315 accessions belonging to five different landraces with corresponding clones in the CSIR-CRI library were identified (Table 1). However, we could not match a total of 202 accessions from farmers’ fields to any of the genotypes in the reference library. These belonged to groups VII, X and XI as well as the various hybrid groups (Table 1).
Geographical distribution of the identified varieties
To further place the results from cluster analyses in a geographical context, we projected individual accessions on the map of Ghana (Fig. 1) using the associated GPS co-ordinates. Each accession is represented by a barplot that is colored according to the inferred membership in the K = 11 genetic clusters. The two most common varieties (I and II) are equally well distributed across the three study regions suggesting they are highly preferred by most farmers and have broad adaptation. On the other extreme, varieties VII and X are geographically restricted and found only in one geographic area.
Variation in the geographic distribution of the released varieties was observed in the three study regions. The most common released variety (Variety II in Fig. 1) is well distributed across the three study regions; Variety III was found mostly in the Brong Ahafo region and a few places in the Eastern Region; Variety IV and VIII occurred mostly in the Eastern Region and in small patches of the other regions (Fig. 1). Potential reasons for the geographic clustering of varieties include region specific uses and adaptation as well as being newly evolved or introduced varieties with limited dissemination opportunities. The location and limited number of industrial processing facilities may also restrict distribution of high yielding varieties suitable for processing.
Correspondence between local names and each of the identified cultivars
While many farmer given variety names correspond to specific clones, there are often differences between genetically unique cultivars as identified by 56,489 SNPs and variety names as elicited from farmers (Fig. 5). For example, the most common clone (Variety I) was variously named as “Debor”, “Bankye Kokoo” and “Ankra”, as well as other less common names not shown in Fig. 5. Spatial distribution analysis revealed that these three most commonly used names are geographically structured by regions (Additional file 3: Figure S2) suggesting there are regional differences in the name of the same variety. The naming system was similarly complex for the eight released varieties cultivated by farmers. For instance, the two most common released varieties (IFAD/UCC and Nkabom) were associated with 33 and 25 different names, respectively (Additional file 1: Table S1). Such discrepencies resulting from synonymy and homonymy in clones names is expected to confound tracking of released varieties when relying on use of names alone.
Performance of reduced numbers of SNP markers in cultivar identification
Accuracy of ancestry estimates (mean R2 and SD) obtained with subsets of SNP markers, selected according to increasing FST (0.1, 0.2,…, 0.9, and 0.95), the corresponding equivalent numbers of randomly sampled SNPs, and the full set of 56,489 SNPs are presented in Fig. 6. Markers passing the predetermined FST thresholds were the complete set of SNPs, 43007, 37900, 30962, 24560, 14426, 5359, 2755, 1392, 570 and 324, respectively. We found that the randomly drawn SNP panels slightly but consistently outperformed the SNP panels selected according to FST (average R2 = 0.97 for FST and R2 = 0.99 for the random SNPs). In addition, each of the 20 independently drawn samples for each subset produced very similar results, as indicated by the small standard deviation. A very high correlation with the full SNP data was obtained using SNPs with FST below 0.6 (5359 SNP) and both random and FST-based subsets performed similarly. Using 2755 SNPs (FST > 0.70) resulted in (0.05 units) lower correlations. The accuracy of 324 SNPs with FST > 0.95 was substantially lower (R2 = 0.90), even much less than that obtained from a similar number but randomly drawn SNPs (R2 = 0.96). Our results suggest there is loss of information in predicting admixture when going for markers with larger FST while random samples of SNPs give higher accuracies, though the actual differences are small. New SNPs-based variety identification studies for cassava would therefore require at least 300 informative SNP markers in order to have sufficient power to not only identify varieties but also estimate ancestries of these accessions.
Cluster analysis using DAPC and validation of ADMIXTURE results
We validated maximum likelihood-based clustering results from ADMIXTURE analysis using DAPC method that is considered free of Hardy-Weinberg and linkage disequilibrium assumptions. Model selection using BIC revealed the presence of hierarchical structure in the population, with steep decline from K = 2 up to around K = 10 followed by a more gentle decrease. The lowest BIC value which corresponded to optimal cluster number was obtained at K = 21 (Fig. 4b). Although this number was larger than that found by ADMIXTURE (K = 11), DAPC clustering recapitulated the groupings uncovered by both the distance-based hierarchical clustering topology as well as ancestry estimates achieved by ADMIXTURE. Comparison of the cluster membership results from the DAPC and ADMIXTURE analyses are summarized in Additional file 4: Table S2. A major difference between the results of the two clustering methods was the propensity of the DAPC analysis to assigned entire individuals to a single cluster compared to ADMIXTURE, which was able to assign admixed individuals to multiple clusters. Indeed, whereas a total of 277 genotypes (or 339 when including the CSIR-CRI library accessions) did not meet the 90 % threshold for belonging to single cluster in the ADMIXTURE, only 15 genotypes did not to not exceed the same threshold in the DAPC analysis. Of the accessions belonging purely to clusters I to XI (i.e. ancestry > 90 %) from ADMIXTURE analysis, we found 100 % agreement with their corresponding DAPC clusters, except for cluster VI whose members were assigned to two DAPC groups (9 and 14) in roughly equal proportions (Additional file 4: Table S2). Large number clusters from DAPC mostly corresponded to sets of genetically similar groups of admixed individuals that shared same ancestries (Additional file 4: Table S2).
As a clonally propagated crop, cassava has several special characteristics for consideration, which also makes it an interesting crop for this case study. First, due to its broad tropical distribution and its predominantly outbreeding system, cassava carries considerable heterozygosity . As expected in typical subsistence farming systems, a substantial number of cassava farmers cultivate more than one variety in their fields to take care of diverse needs. This multiplicity of varieties in farmers’ fields, enables cross-breeding, and eventually some of the volunteer seedlings end-up being selected either consciously or unconsciously as new varieties that are subsequently exchanged . This is the most likely explanation for the occurrence of numerous admixtures in the study region.
Second, as a so-called “orphan crop” , cassava improvement has been mainly implemented by public breeding programs and lacks a formal seed system, thereby making varietal dissemination a challenge . Most farmers use their own planting materials (usually stem cuttings from the preceding crop) or they source stem cuttings from neighboring farmers . Even released varieties may be relatively old due to the low rate of variety turn-over . This allows for spontaneous emergence of clonal variants with different phenotypes that may be undistinguishable by molecular markers.
Genetic distance based on use of molecular markers has been proposed as an appropriate tool to identify putative ‘essentially derived varieties’ [11, 43–45]. The concept of essential derivation is often used in relation to protection of breeder’s rights and refers to a variety with slight modifications from an original variety (such as a single gene insertion through transgenic approaches, back-crossing, or induced mutagenesis) . Our study has followed similar principles, but the objective is different: to assist in collecting accurate variety-specific identification data that can be used to study rates of adoption. However, the success of DNA-based sample identification procedure ultimately depends on the availability of a library panel containing representatives of known varieties. Ideally, the library should be as comprehensive as possible and well curated. In our study, we found instances of accessions from farmers that did not have corresponding genotypes in the library (i.e. Cluster groups VII, X and XI). Moreover, several sets of differently named duplicates were found in the library. Most of these duplicates are classified as landraces and were perhaps independently collected from different regions of Ghana and therefore came with different names.
Next generation sequencing-based genotyping methods such as GBS yield thousands to hundreds of thousands of SNP markers, depending on the genome size, choice of restriction enzymes and the level of sample multiplexing. In the present study, we obtained more than half a million SNPs, which were reduced to about 56,000 markers after curation. Cluster analysis using subsets of either randomly selected or FST-selected SNP markers showed that smaller number of markers could produce similar results to those obtained from complete marker data. An ideal set of ancestry informative SNP markers should have one allele that is fixed in one ancestral lineage and not present in the other . Such sets of markers are designed to provide most of the ancestry information using low density cost-effective SNP genotyping arrays and will be valuable for follow-up studies. Numerous ancestry-specific SNP markers have been developed and used in human population studies [47–49]. However, the number of markers required for population assignment will depend on the populations under consideration, their respective level of genetic differentiation and the desired stringency of assignment. Use of array-based genotyping with a fixed set of pre-selected SNPs thus requires an upfront investment and research to determine the genetic structure of the target study population. A more plausible alternative is to use GBS, which at a higher multiplexing level (for instance 384 DNA samples instead of 96) will be cost-effective enough for direct genotype identification. Although higher multiplexing of samples will proportionately reduce the number of scorable SNPs, it is expected that the final number will still be more than sufficient for cultivar identification. In other words, the increasingly affordable sequencing-based DNA fingerprinting methods should be employed as the primary variety identification tool in collection of variety specific adoption data during household level impact assessment studies.
In the present study, the distance-based approach was successfully used to match accessions from farmers’ fields to corresponding varieties in the ‘library’ of released varieties maintained by CSIR-CRI, based on pairwise distance threshold determined empirically from redundant genotyping of a subset of the collection. We then went further by unraveling the underlying population structure of the studied germplasm with the aim of determining the ancestry of individual accessions. In impact assessment studies, the ancestry information is important since it provides a framework for determining the contribution of specific germplasm in development of new varieties and therefore show indirect impact of germplasm originating from a specific breeding program . This was achieved through the analysis of the populations structure from the high-density SNP data using the complementary model-based methods of ADMIXTURE and discriminant analysis of principal components. In the absence of reliable pedigree records or where varieties are selected from open-pollinated seeds, ancestry analysis from DNA markers is the only way uncovering the genetic source of varieties. The DAPC method uncovered more clusters than ADMIXTURE (Additional file 4: Table S2) but whereas the latter method revealed large number of individuals with two or more ancestries, DAPC mostly assigned individuals to single clusters. This is because the DAPC approach relies on discriminant functions that seeks to maximize the diversity between clusters by while minimizes within-cluster diversity . Such method works best with discontinuous population structure such as in island-model but was found to be less efficient in cassava germplasm due to their continuous and complex population structure [51, 52]. In clonal crop species like cassava, varieties are often derived from complex inter-generational crosses, resulting in clusters that tend to dissolve into clinal patterns of genetic differentiation . Still, DAPC cluster assignment generally agreed with the main ADMIXTURE clusters where >90 % ancestries were assigned to specific clusters. In conclusion, this study confirms the reliability and accuracy of high-density SNP markers from sequencing-based genotyping methods for variety identification and tracking adoption of crop varieties.
Availability of supporting data
The SNP data sets supporting the result of this article are available at www.cassavabase.org (ftp://ftp.cassavabase.org/manuscripts/Rabbi_et_al_2015.zip).
Ancestry Informative Markers
Council for Scientific and Industrial Research - Crops Research Institute
Global Positioning System
International Institute of Tropical Agriculture
Minor allele frequency
Reduced Representation Library
Single Nucleotide Polymorphism
Principal Component Analysis
Fermont AM, van Asten PJA, Tittonell P, van Wijk MT, Giller KE. Closing the cassava yield gap: an analysis from smallholder farms in east africa. Field Crop Res. 2009;112(1):24–36.
Asfaw S, Kassie M, Simtowe F, Lipper L. Poverty reduction effects of agricultural technology adoption: a micro-evidende from rural Tanzania. J Develop Studies. 2012;48(9):1288––1305.
Morris ML, Tripp R, Dankyi AA. Adoption and impacts of improved maize production technology: a case study of the Ghana grains development project. Mexico, D.F: CIMMYT; 1999.
Shiferaw B, Kassie M, Jaleta M, Yirga C. Adoption of improved wheat varieties and impacts on household food security in Ethiopia. Food Policy. 2014;44:272–84.
Elhoumaizi M, Saaidi M, Oihabi A, Cilas C. Phenotypic diversity of date-palm cultivars(Phoenix dactylifera L.) from Morocco. Genet Resour Crop Ev. 2002;49(5):483–90.
Racchi ML, Bove A, Turchi A, Bashir G, Battaglia M, Camussi A. Genetic characterization of Libyan date palm resources by microsatellite markers. 3. Biotech. 2014;4(1):21–32.
Duminil J, Di Michele M. Plant species delimitation: a comparison of morphological and molecular markers. Plant Biosystems. 2009;143(3):528–42.
Yoon MS, Song QJ, Choi IY, Specht JE, Hyten DL, Cregan PB. BARCSoySNP23: a panel of 23 selected SNPs for soybean cultivar identification. Theor Appl Genet. 2007;114(5):885–99.
Jones AG, Ardren WR. Methods of parentage analysis in natural populations. Mol Ecol. 2003;12(10):2511–23.
Morris ML, Heisey PW. Estimating the benefits of plant breeding research: methodological issues and practical challenges. Agr Econ. 2003;29(3):241–52.
Ercisli S, Ipek A, Barut E. SSR marker-based DNA fingerprinting and cultivar identification of olives (Olea europaea). Biochem Genet. 2011;49(9–10):555–61.
Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, et al. A robust, simple Genotyping-by-Sequencing (GBS) approach for high diversity species. PLoS One. 2011;6(5), e19379.
Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q, et al. TASSEL-GBS: a high capacity genotyping - by - sequencing analysis pipeline. PLoS One. 2014;9(2), e90346.
Olsen KM, Schaal BA. Evidence on the origin of cassava: phylogeography of Manihot esculenta. Proc Natl Acad Sci. 1999;96(10):5586–91.
FAO. Why cassava? Rome, Italy: FAO; 2008.
FAOSTAT. Statistical database. Rome, Italy: FAO; 2014.
Okogbenin E, Setter TL, Ferguson M, Mutegi R, Ceballos H, Olasanmi B, et al. Phenotypic approaches to drought in cassava: review. Front Physiol. 2013;4:93.
Mtunguja MK, Laswai HS, Muzanila YC, Ndunguru J. Farmer’s Knowledge on Selection and Conservation of Cassava (Manihot esculanta) Genetic Resources in Tanzania. J Biol, Agriculture HealthCare 2014;4(10):74–78
Angelucci F. Analysis of incentives and disincentives for cassava in Ghana. Technical notes series. MAFAP, FAO, Rome; 2013
Dellaporta SL, Wood J, Hicks JB. A plant DNA minipreparation Version II. Plant Mol Biol Rep. 1983;1:19–21.
Rabbi I, Hamblin M, Gedil M, Kulakow P, Ferguson M, Ikpan AS, et al. Genetic mapping using genotyping-by-sequencing in the clonally propagated cassava. Crop Sci. 2014;54(4):1384–96.
Hamblin MT, Rabbi IY. The effects of restriction-enzyme choice on properties of genotyping-by-sequencing libraries: a study in cassava (Manihot esculenta). Crop Sci. 2014;54(6):2603–8.
Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Software. 2010;33(1):1–22.
Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19(9):1655–64.
Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11(1):1–15.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20(2):289–90.
R Core Team. R: A language and environment for statistical computing. In. Vienna, Austria.URL http://www.R-project.org/. R Foundation for Statistical Computing 2013.
Noli E, Teriaca MS, Conti S. Criteria for the definition of similarity thresholds for identifying essentially derived varieties. Plant Breed. 2013;132(6):525–31.
Alexander DH, Lange K. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinformatics. 2011;12:246.
Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38(6):1358–70.
Liu Y, Nyunoya T, Leng S, Belinsky SA, Tesfaigzi Y, Bruse S. Softwares and methods for estimating genetic ancestry in human populations. Hum Genomics 2013, 7:1. doi:10.1186/1479-7364-7-1
Frichot E, Mathieu F, Trouillon T, Bouchard G, Francois O. Fast and efficient estimation of individual ancestry coefficients. Genetics. 2014;196(4):973–83.
Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24(11):1403–5.
Legendre P, Legendre L. Numerical ecology. 2nd ed. Amsterdam: Elsevier; 1998.
Weiss KM, Long JC. Non-Darwinian estimation: my ancestors, my genes’ ancestors. Genome Res. 2009;19(5):703–10.
Akano O, Dixon O, Mba C, Barrera E, Fregene M. Genetic mapping of a dominant gene conferring resistance to cassava mosaic disease. Theor Appl Genet. 2002;105(4):521–5.
Wang W, Feng B, Xiao J, Xia Z, Zhou X, Li P, et al. Cassava genome from a wild ancestor to cultivated varieties. Nat Commun. 2014;5, 10.1038/ncomms6110
Duputié A, DeléTre M, De Granville J-J, McKey D. Population genetics of Manihot esculenta ssp. flabellifolia gives insight into past distribution of xeric vegetation in a postulated forest refugium area in northern Amazonia. Mol Ecol. 2009;18(13):2897–907.
Fauquet CM, Taylor NJ, Tohme J. The global cassava partnership for the 21st century (GCP21). Trop Plant Biol. 2012;5(1):4–8.
Kyamanywa S, Kashaija I, Getu E, Amata R, Senkesha N, Kullaya A. Enhancing food security through improved seed systems of appropriate varieties of cassava, potato and sweetpotato resilient to climate change in Eastern Africa. Nairobi, Kenya: ILRI; 2011. p. 1–28.
Krishna VV, Spielman DJ, Veettil PC, Ghimire S. An empirical examination of the dynamics of varietal turnover in Indian wheat. Washington, DC: IFPRI; 2014. p. 1–32.
Busti A, Caceres M, Calderini O, Arcioni S, Pupilli F. RFLP markers for cultivar identification in tall fescue (Festuca arundinacea Schreb.). Genet Resour Crop Evol. 2004;51(4):443–8.
Wu B, Zhong G-Y, Yue J-Q, Yang R-T, Li C, Li Y-J, et al. Identification of pummelo cultivars by using a panel of 25 selected snps and 12 dna segments. PLoS One. 2014;9(4), e94506.
Rodrigues DH, Neto FdA, Schuster I. Identification of essentially derived soybean cultivars using microsatellite markers. Crop Breed Appl Biotechnol. 2008;8(1):74—78.
Rosenberg NA, Li LM, Ward R, Pritchard JK. Informativeness of genetic markers for inference of ancestry. Am J Hum Genet. 2003;73(6):1402–22.
Huckins LM, Boraska V, Franklin CS, Floyd JAB, Southam L, GCAN, et al. Using ancestry-informative markers to identify fine structure across 15 populations of European origin. Eur J Hum Genet. 2014;22(10):1190–200.
Kosoy R, Nassir R, Tian C, White PA, Butler LM, Silva G, et al. Ancestry informative marker sets for determining continental origin and admixture proportions in common populations in America. Hum Mutat. 2009;30(1):69–78.
Qin P, Li Z, Jin W, Lu D, Lou H, Shen J, et al. A panel of ancestry informative markers to estimate and correct potential effects of population stratification in Han Chinese. Eur J Hum Genet. 2014;22(2):248–53.
Sawler J, Reisch B, Aradhya MK, Prins B, Zhong G-Y, Schwaninger H, et al. Genomics assisted ancestry deconvolution in grape. PLoS One. 2013;8(11), e80791.
de Oliveira EJ, Ferreira CF, da Silva SV, de Jesus ON, Oliveira GA, da Silva MS. Potential of SNP markers for the characterization of Brazilian cassava germplasm. Theor Appl Genet. 2014;127(6):1423–40.
Kawuki RS, Herselman L, Labuschagne MT, Nzuki I, Ralimanana I, Bidiaka M, et al. Genetic diversity of cassava (Manihot esculenta Crantz) landraces and cultivars from southern, eastern and central Africa. Plant Genet Resour. 2013;11(02):170–81.
We sincerely thank Sharon E. Mitchell and her team at the Genomic Diversity Facility – Institute of Biotechnology, Cornell University, for the genotyping service. Field activities for this study were conducted in partnership with the research and technical staff of the Crops Research Institute of Ghana (Bright B. Pepra, Charles Debrah, Benedicta N. Frimpong, Edem Lotsu, Ohene Djan, Habiba Aggrey, Alex A. Appiah); the Agriculture Innovation Consult – Ghana; IITA (Peter Iluebbey, Afolabi Agboona, Ibrahim Sadeeq). Funding for this study was provided by: The Standing Panel on Impact Assessment (SPIA) of the CGIAR through support from the Bill and Melinda Gates Foundation; and The Roots, Tubers and Bananas program of the CGIAR. We thank Tunrayo Alabi of IITA’s GIS unit for generating the variety distribution map and Oluwafemi Alaba for sample processing and DNA extraction. Critical suggestions from two anonymous reviewers greatly improved the manuscript.
The authors declare that they have no competing interests.
Initiated and designed the research study field sampling plan: MKM, TA, PK, BR; Coordinated project: EP, PK, TA; Designed and led DNA fingerprinting: IYR, MAG; Analyzed data: IYR, PK, GG, PR, TA, EP; Involved and supervised field sample and data collection: JYA, JAM; House-hold survey data curation: BR; Wrote manuscript; IYR, GG, MAG, Revised manuscript: GG, PR, PK, MKM; Trained enumerators: MKM, TA, EP, BR, PK; Lead survey team: JAM, AD; Prepared and programmed survey instruments: BR, MKM. All authors have read and approved the final version of the manuscript.
Sample information associated with the 917 accessions from three regions of Ghana and the library of known varieties from CSIR-CRI. ADMIXTURE-based ancestry estimates according to predefined eleven clusters are also provided. We have also attempted to, as far as the library is concerned, classify each of the farmers’ accessions to their matching released varieties in the CSIR-CRI library. (XLSX 151 kb)
Hierarchical clustering dendrogram of the duplicated DNA samples. The dashed red line indicates the threshold for declaring genetic identity (i.e. distance between this threshold is spurious and results from GBS SNP calling error). (PDF 17 kb)
Strong geographical structure associated with the most common names attributed to Variety I. (PDF 114 kb)
An overview of the distribution of accessions present in the ADMIXTURE clusters and the groups identified by the DAPC. (DOCX 83 kb)
About this article
Cite this article
Rabbi, I.Y., Kulakow, P.A., Manu-Aduening, J.A. et al. Tracking crop varieties using genotyping-by-sequencing markers: a case study using cassava (Manihot esculenta Crantz). BMC Genet 16, 115 (2015). https://doi.org/10.1186/s12863-015-0273-1