# Structure of sweet potato (Ipomoea batatas) diversity in West Africa covaries with a climatic gradient

Kodjo Glato, Atsou Aidam, Ndjido Ardo Kane, Diallo Bassirou, Marie Couderc, Leila Zekraoui, Nora Scarcelli, Adeline Barnaud, Yves Vigouroux

► To cite this version:

Kodjo Glato, Atsou Aidam, Ndjido Ardo Kane, Diallo Bassirou, Marie Couderc, et al. Structure of sweet potato (Ipomoea batatas) diversity in West Africa covaries with a climatic gradient. PLoS ONE, 2017, 12 (5), pp.e0177697. 10.1371/journal.pone.0177697. hal-03344108

HAL Id: hal-03344108

https://hal.umontpellier.fr/hal-03344108

Submitted on 16 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Structure of sweet potato (*Ipomoea batatas*) diversity in West Africa covaries with a climatic gradient

Kodjo Glato<sup>1,2,3,4*</sup>, Atsou Aidam<sup>1</sup>, Ndjido Ardo Kane<sup>2,3</sup>, Diallo Bassirou<sup>2,3</sup>, Marie Couderc<sup>4</sup>, Leila Zekraoui<sup>2,4</sup>, Nora Scarcelli<sup>4</sup>, Adeline Barnaud<sup>2,3</sup>, Yves Vigouroux<sup>4</sup>

1 University of Lomé, Lomé, Togo, 2 Institut Sénégalais de Recherches Agricoles (ISRA), Dakar, Sénégal,

3 Laboratoire Mixte International Adaptation des Plantes et microorganismes aux Stress Environnementaux (LAPSE), Dakar, Sénégal, 4 Institut de Recherche pour le Développement (IRD), Montpellier, France

* glatokodjo@gmail.com

## Abstract

Sub-Saharan agriculture has been identified as vulnerable to ongoing climate change. Adaptation of agriculture has been suggested as a way to maintain productivity. Better knowledge of intra-specific diversity of varieties is prerequisites for the successful management of such adaptation. Among crops, root and tubers play important roles in food security and economic growth for the most vulnerable populations in Africa. Here, we focus on the sweet potato. The Sweet potato (*Ipomoea batatas*) was domesticated in Central and South America and was later introduced into Africa and is now cultivated throughout tropical Africa. We evaluated its diversity in West Africa by sampling a region extending from the coastal area of Togo to the northern Sahelian region of Senegal that represents a range of climatic conditions. Using 12 microsatellite markers, we evaluated 132 varieties along this gradient. Phenotypic data from field trials conducted in three seasons was also obtained. Genetic diversity in West Africa was found to be 18% lower than in America. Genetic diversity in West Africa is structured into five groups, with some groups found in very specific climatic areas, e.g. under a tropical humid climate, or under a Sahelian climate. We also observed genetic groups that occur in a wider range of climates. The genetic groups were also associated with morphological differentiation, mainly the shape of the leaves and the color of the stem or root. This particular structure of diversity along a climatic gradient with association to phenotypic variability can be used for conservation strategies. If such structure is proved to be associated with specific climatic adaptation, it will also allow developing strategies to adapt agriculture to ongoing climate variation in West Africa.

## Introduction

Sub-Saharan agriculture is identified as being particularly vulnerable to ongoing climate change [1]. This is considered unther the climatic context in which African countries have to increase agricultural production to feed their growing populations. Understanding genetic

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Glato K, Aidam A, Kane NA, Bassirou D, Couderc M, Zekraoui L, et al. (2017) Structure of sweet potato (*Ipomoea batatas*) diversity in West Africa covaries with a climatic gradient. PLoS ONE 12(5): e0177697. https://doi.org/10.1371/journal.pone.0177697

**Editor:** Roberto Papa, Università Politecnica delle Marche, ITALY

**Received:** December 26, 2015

**Accepted:** May 2, 2017

**Published:** May 26, 2017

**Copyright:** © 2017 Glato et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This study was partly funded by the Agropolis Foundation's "Agropolis Resource Center for Crop Adaptation and Diversity" (ARCAD) project, the French embassy in Togo (SCAC), Institut Sénégalais de la Recherche Agricole (ISRA), Dakar, Sénégal; the Institut de la Recherche pour le Développement (IRD) in Dakar, Sénégal and Montpellier, France; the West African
Economic and Monetary Union (UEMOA) grant N°
11514/2011/DDS/DESFP/PAES-436 to Kodjo Glato
and the University of Lomé, Togo. We thank C.
Roullier and Professor K. Akpagana for early on
advice on this study, V. Lebot for helpful
comments and the farmers who willingly
participated in the study. The funders had no role in
study design, data collection and analysis, decision
to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared
that no competing interest exist.

diversity, the structure of diversity and the capacity of adaptation of African crops could help
design strategies to face both increased demand and a changing environment. Tuber and root
crops represent a large proportion of the diet in Africa and consequently are very important
for such an endeavor [2]. Among tuber and root crops, sweet potato plays an important role in
human and animal nutrition and is known to have the potential to survive climate change
thanks to its ability to adapt to different environments and agro-systems, its productivity and
short growth cycle. In addition, sweet potato is receiving increased attention as a rich source of
vitamin A [3, 2]. Sweet potato is a staple food, ranking seventh among food crops in annual
production in the world [4] and fifth for its caloric contribution in developing countries after
rice, wheat, maize and cassava [5]. Although Africa occupies the second position, annual pro-
duction (19%) remains far below that in Asia, which alone produces 75% of the world annual
production [6]. Cultivation of the sweet potato is still low in West Africa, and lacks an inven-
tory of diversity and a conservation plan, as well as a study of its sensitivity to pests and dis-
eases [7].

Local landraces of sweet potato need to be collected, characterized and conserved to
improve African production. The sweet potato originated in tropical America [8, 9] and was
dispersed around the world principally by human mediated migration; and is now cultivated
in most tropical regions [10]. It was introduced to Africa by the Portuguese in the 16th century
[11, 12]. Probably first introduced as a tuber in Tanzania [10], the plant then dispersed from
East to West Africa [10]. Sub-Saharan Africa is considered to be a secondary center of diversity
[13]. However, its genetic characterization is still partial.

Among the tools used for assessing the diversity of African sweet potato, agro-morphologi-
cal analysis was conducted in Tanzania [14], Kenya [15] and Burkina-Faso [16]. Unfortu-
nately, morphological characteristics can vary because of trait plasticity [17] and it has been
shown that the environment can affect the phenotypic traits [18]. Consequently, comparisons
between studies of morphological diversity are difficult. Molecular markers make it possible to
conduct such comparisons across countries and continents. Analysis of sweet potato diversity
using SSRs markers was recently initiated [9, 19, 20].

The questions we addressed in this study are: i) how is sweet potato structured along the
environmental gradient from the Tropics to the Sahel? ii) Are some genetic groups associated
with specific phenotypic differences? To understand the genetic structure and diversity of
sweet potato in West Africa, we focused on two countries, Togo and Senegal, which cover a
range of climates. We also performed a phenotypic evaluation to compare genetic structure
and morphological differences between samples. This analysis enabled us to identify specific
genetic groups showing phenotypic differences and occurring in different climatic conditions.

**Materials and methods**

Plant collection

A total of 132 accessions were collected: 112 from Togo in 2012, and 20 from Senegal in 2014 (S1 Table). The collection was done with permission of doctors Bassirou Diallo (ISRA) from Senegal and Nenonene Amen Yawo (ESA) from Togo. Samples were collected with farmers consent. The aim of the collection was to maximize the variability observed in each village [21]. A mean of 3.56 samples was collected per village. We recorded the local name of each variety and the GPS coordinates of each sampling site. For each accession, a stem was sampled and kept isolated for further genetic and phenotypic experiments. In Senegal, we collected accessions in four villages in the region where most sweet potato is cultivated (region of Saint-Louis). We used the same sampling protocol as in Togo, except that young leaves were
collected instead of stems. The GPS coordinates of each sampling site were recorded. A mean of 5 samples were collected per village

## Phenotyping

Field trials were conducted only on Togo accessions in the experimental Agronomic station of University of Lomé (altitude 27 masl. 06°10′27.2″ E 001°12′ 39.6″N). The field studies did not involve endangered or protected species. The soil in the field is classified as ferralitic. The climate is characterized by two rainy seasons, leading to two periods of sweet potato cultivation (April-July and September-December). From one to three nodes of each of the 112 accessions were planted (236 cuttings with an average 2.10 per accession). The trial was performed on a 300 m² plot (30 m × 10 m). The space between plants and between rows was 1 m. The experiment was repeated three times, once in 2012 and twice in 2013. The stems from the 2012 field experiment were harvested for use in the 2013 field experiment. The field was weeded manually with a hoe and watered if necessary. The morphological characteristics of the leaf and stem were recorded 90 days after planting.

A total of 21 morphological descriptors [22], were used as agro-morphological characteristics of these three organs: leaf, stem, and root. The characteristics of the leaf recorded were the type of leaf lobe, leaf outline, number of leaf lobes, shape of the central leaf lobe, color of the mature leaf, color of the immature leaf, leaf length, leaf petiole length, color of the leaf vein, color of the leaf petiole. The characteristics of the stem recorded were vine internode length, diameter between nodes, color, port and pubescence. As roots are the idle part of the plant, we collected root characteristics at maturity. Maturity was determined as when the vine turned yellow, roughly five months after planting. The root traits recorded were root shape, main color of the skin, secondary color of the skin, main color of the flesh, secondary color of the flesh, distribution of the secondary color of the flesh.

## DNA extraction and SSR genotyping

We genotyped all 132 individuals sweet potato plants from Togo and Senegal (S1 Table). Leaves were collected from the Togo samples during the 2012 field experiment. DNA was extracted from dried young leaves (10 to 25 mg) following the method of Doyle and Doyle [23] adding 2% (w/v) of cetyl-tri-methyl-ammonium bromide (CTAB) and 2% (w/v) of sodium sulfite. The pellet was left to dry at room temperature and resuspended in 100 µl of 1x Tris EDTA buffer (TE). The concentrations were assessed using NanoDrop 2000/2000C. DNA was checked on 0.8% agarose gel. The DNA solution was stored at -20°C until genotyping.

Twelve microsatellite (SSRs) previously designed for sweet potato [24] were used (S2 Table). The amplification conditions were slightly modified and the fluorochromes for each marker (new conditions listed in S2 Table). All loci were independently amplified using the multiplex PCR kit (Qiagen) in a final volume of 10 µl; 5 µl master mix (1x); 2 µl Q solution (1x); 2 µl DNA; 0.2 µl R primer; 0.2 µl F primer and 0.6 µl Rnase free solution. The PCR amplification program (S2 Table) was from Roullier et al.[24] with slight modifications. The PCR was performed on a Techné Gradient TC Plus Thermocycler using the following program: 15 min at 95°C, 35 cycles of 30 s at 94°C, 90 s at hybridization temperature varying from 58°C to 60°C (S2 Table), 60 s at 72°C and 30 min at 60°C. Amplifications were performed on 96-well plates with two negative water controls and four positive DNA controls. The amplification product (4 µl) was checked on 2% agarose gel. SSRs were diluted and pooled for migration on an ABI 3130 XL 16 capillary sequencer. GeneMapper software version 4.1 was used to score alleles. Two different people scored the alleles. Since reading of size might be difficult for hexaploid plants [9], we performed a second independent run and an independent score of alleles.
The second dataset was used to assess the reproducibility of the results. Analyses were computed independently on the two datasets, and overall agreement was shown between datasets (see details in the supplementary file). In addition, we used previously genotyped accessions from Tropical America and Oceania [24]. To allow comparison with our sampling number, we randomly sampled 132 accessions both for Tropical America and Oceania. The number of alleles recorded was used to statistically compare diversity in the three areas, West Africa, Tropical America and Oceania.

## Genetic data analyses

**Diversity.** Each SSR allele was binary coded (0 for absence and 1 for presence) following previous approaches used for this hexaploid species [24]. We first compared the diversity observed in Africa with the diversity found in America and Oceania using the same number of samples randomly selected from Roullier *et al.* [25], using both the number of alleles and polymorphic information content (PIC). A Wilcoxon paired test with the number of alleles was used to assess differences in diversity. The PIC of the African accessions was calculated using the binary matrix. For this hexaploid species, the PIC was calculated considering an allele as a locus: $PIC = 1 - (P_1^2 + P_2^2 + 2P_1^2P_2^2)$ where $P_1$ and $P_2$ corresponded to the frequency of the presence (1) and absence (0) of an allele. So, $PIC = 1 - P_1^2$ with $P_1$ the frequency of allele i.

The neighbor joining tree method was used with the Euclidean genetic distance among accessions as implemented in PowerMarker 3.25 [26]. The phylogram was visualized using the international tree of life: Itoel Program v2 [27].

**Genetic structure.** Discriminant analysis of principal components (DAPC) implemented in the ADENEGET package [28] using R software 3.1.1 [29] was used to assess the distribution of genetic diversity. For DAPC, we first performed a principal component analysis (PCA) to capture the large number of allele absence/presence (0/1) matrices on the most significant PCA axes. The number of PCA axes retained was based on (1) the a score optimization procedure proposed in the ADEGENET package [28], (2) the cumulative percentage of variance explained by the first selected axes. We then assessed the optimal number of clusters supported by the genetic data using the k-means procedure [28]. A Bayesian information criterion (BIC) was calculated with each model ($K = 1$ to 20). Membership probabilities for each individual were plotted using the compoplot function of R. We also used STRUCTURE 2.3.4 software [30, 31]. We ran the admixture model with a number of populations from $K = 1$ to 20 with 15 runs using a 100.000 burn in period and 400.000 Markov Chain Monte Carlo (MCMC) steps.

**Analysis of the composition of the genetic groups.** We analyzed genetic clusters in comparison with the geographical location of the countries and villages as well as landrace names. To assess relationships between agro-morphological traits and genetic clustering, we performed Kruskall-Wallis nonparametric tests. We ran this test with the 21 morphological traits recorded. We considered the genetic group to be associated with morphological differences when the test was significant with a Bonferroni corrected p-value of 0.0024. Wilcoxon paired tests were then performed for each significant trait to assess which specific genetic group differed significantly from the others.

We also assessed the relationship between the clustering and climate variables. To do so, for each accession we downloaded the BIOCLIM data [32] composed of 19 bioclimatic variables and monthly climate data (temperature and precipitation). We performed a principal component analysis to obtain synthetic uncorrelated climatic contrast. The PCA axes were then used to check for a relationship between climate variables and inferred cluster groups using Kruskall-Wallis nonparametric tests.
# Results

## Diversity of sweet potato in West Africa

A total of 90 alleles were observed. The number of alleles per locus ranged from 4 to 12 with a mean of 7.5. Locus ib292 had the highest number of alleles and loci J263 and J315 the lowest (S3 Table). The size of the amplified products ranged from 83 bp to 250 bp. Allelic frequency among the accessions ranged from 0.007 to 0.990 and PIC from 0.11 to 0.31 with a mean of 0.19 (Fig 1). The 11 markers used in our study were the same as those used in a previous study [24]. With these 11 SSR markers, we found an average of 7.36 alleles for the West African sample, 9.09 alleles for the Oceanian sample and 10.90 alleles for the American sample (Fig 2). The number of alleles ranged from seven to 22 for Tropical America accessions, four to 17 for Oceanian accessions and four to 12 for West African accessions (S4 Table). Allelic frequency ranged from 0.007 to 0.990 in the three regions. The Wilcoxon paired test revealed significant differences in diversity between America and West Africa (p-value = 0.0083) and Oceania (p-value = 0.048). However, no significant difference was found between West Africa and Oceania (S5 Table).

## Structure of sweet potato populations

With the a-score optimization function applied to the dataset, the number of PCA to be retained ranged from two to four. The cumulative variance explained by the 12 first PCs axes was 74.7%, and consequently we retained these 12 axes (S1 Fig). The Bayesian information criterion (BIC) values assessed from two to 12 clusters decreased steadily with increasing K

<figure><img src="image_4.png" /><figcaption>**Fig 1. Polymorphic Index Content (PIC) of SSR markers.** The values indicate the PIC of the loci used to analyze our 132 African accessions. Each locus value corresponds to the mean PIC of all the alleles of the corresponding markers. The last bar represents the mean PIC of the 12 markers.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0177697.g001
Mean number
of alleles

<figure><img src="image_5.png" /><figcaption>**Fig 2.** Mean number of alleles in the different regions. America has the highest number of alleles, approximately 11, followed by Oceania with 9 and Africa with 7.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0177697.g002

(S2 Fig). We assessed the composition of the group at different K values (3, 4, 5, 6) and K = 5 appeared to be the best number of genetic groups to explain the genetic data on the basis of significant statistical analysis, individual position in the group, and membership probability. Indeed, the DAPC clearly separated individual into groups with K = 3, 4, and 5. But with K = 6, the subdivision of the groups was less clear (S3 Fig). The analysis using STRUCTURE tend also favored five clusters (S4 Fig). In terms of clustering between DAPC and STRUCTURE, the inferences of ancestry between groups were very similar (S6 Table).

When K = 5 (Fig 3), the first group contained 18 accessions, 15 from Togo and three from Senegal. The second group, called “Tropical accessions”, contained 27 accessions all from Togo. The third group was composed of 60 accessions of which only four were from Senegal. All the accessions from the northern part of Togo were in this group, named “Tropical arid accessions”. Group four contained eight accessions all from Senegal called “Sahelian accessions”. Group five contained 19 accessions including five from Senegal and the rest from Togo. Group one and group five were composed of accessions with no clear specific geographical origin, called “Tropical-Sahelian accessions group 1” and “Tropical-Sahelian accessions group 2” respectively. The membership probabilities of each individual in each group were very high for each accession (S5 Fig, S6 and S7 Tables). The result of the phylogram among accessions was also in agreement with the five genetic groups (Fig 4).

Local name and genetic groups

A total of 31 different local names were recorded: 18 in Togo and 13 in Senegal. The names given to the accessions by the farmers did not match the genetic groups. Indeed, 76% of accessions with the same local name were in fact attributed to different groups. For example, four
<figure><img src="image_6.png" /><figcaption>**Fig 3.** **DAPC genetic groups with k = 5.** In blue: Tropical-Saharan accessions group1 with 18 accessions (15 from Togo and 3 from Senegal); in grey: Tropical accessions with 27 accessions (all from Togo); in yellow: Tropical arid accessions with 60 samples (56 from Togo and 4 from Senegal); in pink: Sahelian accessions with 8 accessions (all from Senegal) and in red, Tropical-Saharan accessions group 2 composed of 19 accessions (14 from Togo and 5 from Senegal).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0177697.g003

accessions called “dankadou” were distributed between the Tropical arid accessions group and
Tropical-Saharan accessions group 2. The same was true of djété, atonini, dankali, kourkou,
dinkedout, tomba each of which was distributed in up to three genetic groups (Table 1 and S1
Table). The local name and the genetic group agreed in only 25% of the cases.
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Dendrogram of NJ cluster analysis of 132 sweet potato accessions. A represents Tropical-Sahelian accessions group 1; B Tropical accessions, C Tropical arid accessions, D Sahelian accessions, and E Tropical-Sahelian accessions group 2. The number indicates the number of accessions. The bootstrap value is given on the branch of the tree.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0177697.g004

## Village richness

Different genotypes from different genetic clusters were generally found within the villages. No village had all five genetic groups (Table 1 and S1 Table). Among the 37 villages sampled, only one had four genetic groups, five villages had three genetic groups and 15 villages had two genetic groups. The Tropical arid accessions group was the most widespread across villages and accessions from this group were found in 67% of the villages (Table 1, S1 Table and Fig 5). The least represented group was “Sahelian accessions”, which was found in only 8% of the villages. The other tropical accessions, “Tropical-Sahelian accessions group1” and “Tropical-Sahelian accessions group 2”, were present in 25% of the villages.

## Genetic group and morphological traits

The different genetic groups were also distinguished by the 21 agro-morphological traits (Table 2). For successive values of K from 3 to 5, six agro-morphological traits distinguished a genetic group whatever the K-value: the main color of the skin of the root, leaf outline, the
<table><caption>Table 1. Distribution of the genetic group per village.</caption><thead><tr><th>Villages</th><th>Latitude</th><th>Longitude</th><th>Country</th><th>I</th><th>II</th><th>III</th><th>IV</th><th>V</th></tr></thead><tbody><tr><td>GBAVE</td><td>N 06'20'08.0"</td><td>E 001°01'40.0"</td><td>TOGO</td><td></td><td></td><td>1</td><td></td><td></td></tr><tr><td>AGBODJEKPOE</td><td>N 06'39'57.8"</td><td>E 001° 09'44.7"</td><td>TOGO</td><td>1</td><td>1</td><td>4</td><td></td><td></td></tr><tr><td>GNAMADJI</td><td>N 06' 15'24.4"</td><td>E 001°19'05.3"</td><td>TOGO</td><td>4</td><td>3</td><td>8</td><td></td><td></td></tr><tr><td>GBATOKOPE</td><td>N 06'14'31.6"</td><td>E 001° 32'05.5"</td><td>TOGO</td><td></td><td></td><td>1</td><td></td><td></td></tr><tr><td>BADJA</td><td>N 06'23'81.8"</td><td>E 000'59'67.6"</td><td>TOGO</td><td></td><td></td><td>1</td><td></td><td>1</td></tr><tr><td>ATSANSI DEDZI</td><td>N 06'27'26.9"</td><td>E 001'32'50.9"</td><td>TOGO</td><td>1</td><td>1</td><td>1</td><td></td><td></td></tr><tr><td>ASSAHOUN</td><td>N 06'27'56.0"</td><td>E 000'54'16.1"</td><td>TOGO</td><td></td><td>2</td><td></td><td></td><td></td></tr><tr><td>MOM-HAGOU</td><td>N 06'29'30.8"</td><td>E 001'33'30.9"</td><td>TOGO</td><td></td><td></td><td>1</td><td></td><td></td></tr><tr><td>KAGNIKPEDJI</td><td>N 06'38'48.4"</td><td>E 001'10'22.4"</td><td>TOGO</td><td>1</td><td></td><td></td><td></td><td></td></tr><tr><td>BAKAKOPE</td><td>N 06'39'43.5"</td><td>E 000'54'20.5"</td><td>TOGO</td><td>1</td><td></td><td></td><td></td><td></td></tr><tr><td>ADRALAKOPE</td><td>N 06'41'22.0"</td><td>E 001'08'59.1"</td><td>TOGO</td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>GNIGBE</td><td>N 06'42'06.6"</td><td>E 001'08'14.3"</td><td>TOGO</td><td></td><td>2</td><td>8</td><td></td><td>2</td></tr><tr><td>AGAMAHE</td><td>N 06'42'40.5"</td><td>E 001'10'19.1"</td><td>TOGO</td><td></td><td></td><td></td><td></td><td>2</td></tr><tr><td>ADAKAPE</td><td>N 06'49'04.1"</td><td>E 001'10'52.5"</td><td>TOGO</td><td></td><td></td><td></td><td></td><td>1</td></tr><tr><td>KPALIME TSEVIE</td><td>N 06'55'22.4"</td><td>E 000'38'31.0"</td><td>TOGO</td><td></td><td>1</td><td>4</td><td></td><td></td></tr><tr><td>YOKELE</td><td>N 06'56'42.1"</td><td>E 000'39'45.5"</td><td>TOGO</td><td></td><td>4</td><td></td><td></td><td></td></tr><tr><td>AKATA</td><td>N 07'02'14.6"</td><td>E 000'42'20.5"</td><td>TOGO</td><td>1</td><td>3</td><td>1</td><td></td><td></td></tr><tr><td>KPELE-TSIKO</td><td>N 07'07'37.6"</td><td>E 000'42'13.3"</td><td>TOGO</td><td></td><td>2</td><td></td><td></td><td>1</td></tr><tr><td>DANYI-N'DIGBE</td><td>N 07'08'30.1"</td><td>E 000'40'35.5"</td><td>TOGO</td><td>2</td><td></td><td></td><td></td><td></td></tr><tr><td>DANYI PEYEYEME</td><td>N 07'12'42.1"</td><td>E 000'41'56.2"</td><td>TOGO</td><td></td><td>1</td><td>3</td><td></td><td></td></tr><tr><td>KPETE</td><td>N 07'25'05.6"</td><td>E 000'52'57.6"</td><td>TOGO</td><td></td><td>1</td><td>1</td><td></td><td></td></tr><tr><td>TCHEBEBE</td><td>N 08'26'16.0"</td><td>E 000'59'30.4"</td><td>TOGO</td><td></td><td></td><td>2</td><td></td><td>3</td></tr><tr><td>TCHALO</td><td>N 08'55'43.9"</td><td>E 001'07'03.9"</td><td>TOGO</td><td></td><td>4</td><td>3</td><td></td><td></td></tr><tr><td>BOTCHOLEYO</td><td>N 08'56'02.0"</td><td>E 001'05'48.8"</td><td>TOGO</td><td></td><td></td><td>1</td><td></td><td></td></tr><tr><td>AMAOUDE</td><td>N 09'08'44.2"</td><td>E 001'09'24.0"</td><td>TOGO</td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>ATCHANGBADE</td><td>N 09'28'39.4"</td><td>E 001'08'17.2"</td><td>TOGO</td><td></td><td></td><td>1</td><td></td><td>1</td></tr><tr><td>DJAMDE</td><td>N 09'30'50.6"</td><td>E 001'02'40.0"</td><td>TOGO</td><td>3</td><td></td><td>2</td><td></td><td></td></tr><tr><td>TANDJOUARE</td><td>N 10'40'15.4"</td><td>E 000'12'19.6"</td><td>TOGO</td><td></td><td>1</td><td>2</td><td></td><td></td></tr><tr><td>NAKI-EST</td><td>N 10'43'31.3"</td><td>E 000'22'37.6"</td><td>TOGO</td><td></td><td></td><td>2</td><td></td><td>1</td></tr><tr><td>NABOULPIONGUE</td><td>N 10'56'20.8"</td><td>E 000'08'01.3"</td><td>TOGO</td><td></td><td></td><td>2</td><td></td><td>1</td></tr><tr><td>GABONBONG</td><td>N 11'01'22.1"</td><td>E 000'04'56.7"</td><td>TOGO</td><td></td><td></td><td>3</td><td></td><td>1</td></tr><tr><td>KASSENA</td><td>N 08'52'60"</td><td>E 001'04'60"</td><td>TOGO</td><td>1</td><td></td><td></td><td></td><td></td></tr><tr><td>KPADAPE</td><td>N 06'51'00"</td><td>E 000'36'00"</td><td>TOGO</td><td></td><td>1</td><td></td><td></td><td></td></tr><tr><td>GUIDICK</td><td>N 16'07'16.4"</td><td>W 015'53'52.9"</td><td>SENEGAL</td><td></td><td></td><td>2</td><td></td><td>1</td></tr><tr><td>SANEITE</td><td>N 16'14'06"</td><td>W 015'47'30.2"</td><td>SENEGAL</td><td>2</td><td></td><td>2</td><td>3</td><td>4</td></tr><tr><td>MBANN</td><td>N 16'17'02.2"</td><td>W 015'47'02.9"</td><td>SENEGAL</td><td>1</td><td></td><td></td><td>1</td><td></td></tr><tr><td>DAGANA</td><td>N 16'32'27.1"</td><td>W 015'30'35.6"</td><td>SENEGAL</td><td></td><td></td><td></td><td>4</td><td></td></tr></tbody></table>

Each village is shown together with its passport data and country. The color of the genetic groups (I to V) shows if it is present in each village. The number of samples per genetic group in each village is also shown.

https://doi.org/10.1371/journal.pone.0177697.t001

number of lobes, the type of central leaf lobe, the color of the leaf veins and the color of the immature leaf. When K = 3, seven traits were associated. In addition to the six common traits, we found the color of the leaf petiole. For K = 4 and K = 5, nine traits were associated. In each case, in addition to the six common traits, we found the secondary color of the skin, pubescence, type of leaf lobe for K = 4, and petiole color, pubescence, and type of leaf lobe type for K = 5 (Table 2).
<figure><img src="image_8.png" /><figcaption>**Fig 5. Geographical distribution of the genetic groups.** In the villages sampled, each color represents one genetic group. The number of colors in each circle corresponds to the number of genetic groups found in the village concerned. The size of the circle is proportional to the number of samples. The total surface area of each color is also proportional to number of samples in that genetic group.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0177697.g005

Wilcoxon-paired test tested with *K* = 5 genetic groups, revealed specific agro-morphological traits for each genetic group (S8 Table). The trait was considered as specific to a group when the p-value of the test was < 0.05 for paired tests. The result showed that “Tropical Sahelian accession group 2” was distinguished by the main color of the skin of the root and the color of the immature leaf. “Tropical accession” also had very specific traits: the number of lobes, type of central lobe of the leaf, and pubescence (S8 Table).

## Genetic group and climate

We found a significant relationship between the cluster and climatic variables. The PCA on the climatic variable retrieved the three first axes with of 52.6%, 24.2% and 18.4% variance explained, respectively. The first axis mainly represents the south/north gradient from a tropical humid zone (Danyi-Peyeyeme) to the northern Togo tropical arid zone (Village of Gabonbong), and then to the drier part of Senegal (S6 Fig).

We found an association between the first axis of the PCA and the DAPC group (Kruskal-Wallis, *p* < 10⁻⁶). This relationship held even when we only considered the accession from Togo (Kruskal-Wallis, *p* < 0.002). Only the third axis was significant when the whole sample was tested (Kruskal-Wallis, *p* < 0.001) but not if we only included Togo (Kruskal-Wallis, *p* < 0.26).
<table><caption>Table 2. Kruskall-Wallis test between genetic group and agro-morphological traits.</caption><thead><tr><th>Traits recorded</th><th>Codes</th><th>K = 3</th><th>K = 4</th><th>K = 5</th></tr></thead><tbody><tr><td>Shape of the roots</td><td>FGt</td><td>0.41</td><td>0.53</td><td>0.57</td></tr><tr><td>Main color of the skin</td><td>CPt</td><td>2.5*10<sup>-6</sup></td><td>2.9*10<sup>-6</sup></td><td>1.9*10<sup>-6</sup></td></tr><tr><td>Secondary color of the skin</td><td>CSt</td><td>5.3*10<sup>-3</sup></td><td>3.3*10<sup>-4</sup></td><td>0.01</td></tr><tr><td>Main color of the flesh</td><td>CPc</td><td>0.047</td><td>0.01</td><td>7.2*10<sup>-3</sup></td></tr><tr><td>Secondary color of the flesh</td><td>Csc</td><td>0.15</td><td>0.02</td><td>0.08</td></tr><tr><td>Distribution of secondary color of the flesh</td><td>Dcsc</td><td>0.76</td><td>0.86</td><td>0.92</td></tr><tr><td>Length of vine internodes</td><td>Lnoe</td><td>0.58</td><td>0.62</td><td>0.78</td></tr><tr><td>Diameter between nodes</td><td>Dianoe</td><td>0.03</td><td>0.08</td><td>0.11</td></tr><tr><td>Color of the stem</td><td>Cpliane</td><td>0.01</td><td>0.01</td><td>0.03</td></tr><tr><td>Port</td><td>Port</td><td>0.48</td><td>5.2*10<sup>-3</sup></td><td>0.06</td></tr><tr><td>Pubescence</td><td>Psomet</td><td>0.04</td><td>2.07*10<sup>-7</sup></td><td>2.8*10<sup>-7</sup></td></tr><tr><td>Type of leaf lobe</td><td>Tylobe</td><td>2.9*10<sup>-3</sup></td><td>1.7*10<sup>-3</sup></td><td>6.1*10<sup>-4</sup></td></tr><tr><td>Leaf outline</td><td>FGfeuille</td><td>3.0*10<sup>-4</sup></td><td>1.3*10<sup>-4</sup></td><td>1.3*10<sup>-4</sup></td></tr><tr><td>Number of leaf lobes</td><td>Nlobe</td><td>4.8*10<sup>-4</sup></td><td>2.4*10<sup>-4</sup></td><td>6.3*10<sup>-5</sup></td></tr><tr><td>Shape of central leaf lobe</td><td>Flobec</td><td>1.6*10<sup>-3</sup></td><td>4.3*10<sup>-4</sup></td><td>1.3*10<sup>-4</sup></td></tr><tr><td>Color of mature leaf</td><td>Cfm</td><td>0.74</td><td>0.24</td><td>0.5</td></tr><tr><td>Color of immature leaf</td><td>Cfim</td><td>1.1*10<sup>-4</sup></td><td>1.5*10<sup>-4</sup></td><td>1.0*10<sup>-4</sup></td></tr><tr><td>Length of leaf</td><td>Lfeui</td><td>0.24</td><td>0.25</td><td>0.25</td></tr><tr><td>Length of leaf petiole</td><td>Lpétio</td><td>0.07</td><td>0.03</td><td>0.01</td></tr><tr><td>Color of leaf vein</td><td>CNerv</td><td>5.2*10<sup>-4</sup></td><td>1.1*10<sup>-4</sup></td><td>2.3*10<sup>-4</sup></td></tr><tr><td>Color of leaf petiole</td><td>Cpétiole</td><td>3.5*10<sup>-5</sup></td><td>3.3*10<sup>-3</sup></td><td>1.2*10<sup>-4</sup></td></tr><tr><td>Total</td><td></td><td>7</td><td>9</td><td>9</td></tr></tbody></table>

Twenty-one agro-morphological traits were used. The p-value for each trait is reported. The colored boxes show significant p-values < 0.0024 (Bonferroni corrected p-value) for the traits. The number of group was assumed according to K value: from K = 3 to K = 5.

https://doi.org/10.1371/journal.pone.0177697.t002

# Discussion

## Relatively low genetic diversity in West African sweet potato

The SSR markers revealed a high level of polymorphism. These values are comparable to those reported in similar studies on the sweet potato in East Africa [33], with a mean of 6.1 alleles per locus compared to 7.4 found in the present study. Yada *et al.* [18] reported an average of four alleles per locus in 192 accessions from Uganda with 10 SSRs used. Gwandu *et al.* [34] found a mean of four alleles in 57 sweet potato from Tanzania. Based on the mean allele value per locus, which was 7.4 in the present study and those listed elsewhere, we tend to have a higher level of diversity in West Africa. However, one should be cautious, because different sets of SSR markers are used, and type and number of repetitions for the microsatellite motif impact the number of alleles and consequently the diversity assessment [35] Nevertheless, a higher diversity in our study is not unexpected, because we sampled a greater climatic contrast in West Africa and more than one country. Indeed, the climate in Togo is contrasted and diversified: arid in the north and humid in the south. In Senegal (Saint Louis) the climate is Sahelian. In Africa, the main hypothesis is that sweet potato was introduced by explorers from Spain and Portugal in the 16th century in east Africa in Tanzania [11, 12, 10]. Under this scenario, West Africa sweet potato is derived from East African sweet potato. The lower diversity in East Africa than in West Africa suggests that sweet potato was not simply a sub-sample from East Africa but might have been re-introduced into West Africa later on. However, confirming this hypothesis would require not only a similar set of markers but also a very similar
sampling strategy. A thorough analysis of sweet potato diversity of the whole Africa remains to be done.

## The genetic groups are partly structured by climate

The genetic groups were closely linked to geographic regions with different climatic conditions. Two clusters were specific to a country: the group called “Sahelian accessions” was linked to Senegal and the group called “Tropical accessions” to Togo. The three other groups were roughly associated with climate: “Tropical-arid accessions”, “Tropical-Saharan accessions group 1” and “Tropical-Saharan accessions group 2”. In the tropical-arid accessions, most of the accessions sampled in Togo were from the northern part of the country, where the climate is arid. The two other clusters “Tropical-Sahelian accessions group 1” and “Tropical-Sahelian accessions group 2” contained accessions either from Togo or Senegal and the climatic conditions were less well defined.

The two groups “Tropical accessions” and “Sahelian accessions” were certainly shaped by ecological conditions, and might reflect specific adaptation to climate and soil. However, genetic bottlenecks created during the introduction and later selection may also have played a role in this genetic differentiation. Environmental and climatic variations have also been shown to be associated with the diversity of cassava in Amazonia [36].

The three other genetic groups: “Tropical arid accessions”, “Tropical-Sahelian accessions group 1” and “Tropical-Saharan accessions group 2” contained accessions from both Togo and Senegal. The relationship between a genetic group and specific geographical features has been previously reported [24]. Some studies done in a smaller region did not always report such results, for example in Uganda [18], in Mozambique [37] and in Puerto Rico [38]. However, we covered a wide range of climates, from the tropical humid area to the Sahel, and we consequently have a good power to show such genetic, geographic and climatic congruence.

## Village richness within genetic groups

We observed that each village cultivated accessions from different genetic groups, but did not differentiate the varieties with specific local names. This can be explained by frequent exchanges between villages in Togo because farmers go to other villages to look for work. They return home with stems they found that had attractive traits or with roots of a plant children like because of its sweetness. Such practices both increase and spread genetic diversity between villages and are particularly important in preserving plant diversity and avoiding genetic erosion. The geographical distribution of allelic diversity appears to be a practical and cost-efficient way to preserve genetic diversity [2].

## Do local names distinguish true genetic groups?

Our study has shown that local names used for sweet potato landraces were not correlated with genetic groups (S1 Table). However, the genetic group was well correlated with very specific morphological traits. The way the plant is named and its importance to the farmers probably explains this observation. As sweet potato is propagated vegetatively, farmers use phenotypic traits [9] or organoleptic traits to differentiate and name their varieties [9, 39]. Remarkably, the number of names recorded in our study was relatively small, with only 31 local names for 132 accessions. Etymologically, only six out of 31 appear to be directly related to names of specific varieties. The main trait used to identify and or name varieties was the color of the skin of the tuber and the color of the vine. Indeed, farmers use phenotype to distinguish the plant [16]. In our study, the farmers were sometimes unable to identify the sweet potato varieties from above ground plant parts alone (vine or leaves) without the root tuber.
Similar observations were made by Shigeta [40] for sweet potato landraces in Ethiopia, and by Elias *et al.* [41], for cassava varieties among Makushi Ameridins in Guyana. But, most of the names used for sweet potato varieties clearly identify the species in different languages. Sweet potato is not a major crop in Togo or Senegal, and some farmers may thus not pay much attention to the phenotypic variability of the plant.

## Genetic clusters and morphologic traits

In this study, we demonstrated that the morphological traits of sweet potato were associated with genetic groups. The most significantly associated traits are qualitative like color of the skin root and quantitative like number of leaves. Consequently, a different set of traits found in each genetic group may reflect random variation of phenotype or selection by farmers. At this stage, we could not clearly differentiate these two hypotheses. The importance of these phenotype traits for breeders is mentioned in many studies [16, 42]. Identifying specific morphological traits related to the genetic groups is a very useful way to study diversity in the field, and morphology could help design conservation strategies that maximize genetic diversity.

## Conclusion

The genetic diversity of West African sweet potato is structured by geography, climate, genetics and morphology. This original structure (genetic, morphologic and climate) allows us to suggest the hypothesis that a genetic group associated today with a specific set of climate variable, might performed well in geographic area showing projected future climate with similar characteristics. This hypothesis of specific climatic adaptation could be tested and if validated suggest possible strategy of crop adaptation to future conditions. The high diversity observed in West Africa might suggests that sweet potato was also introduced in the region rather than simply being the result of diffusion after its introduction in East Africa. This suggest another interesting hypothesis to be tested about the diffusion of sweet potato in Africa. More thorough sampling throughout Africa is now needed to assess such hypotheses.

## Supporting information

**S1 Fig. Method used to choose the number of PCA axes.** The number of PCA chosen could be based on (a) a score function which indicated the appropriate number of PCA axes to be retained. In our case, the graph proposed 3 axes (PCs); or (b) the percentage of variance explained; 74.7% of variance was chosen which corresponded to 12 axes (PCs).
(TIFF)

**S2 Fig. The choice of the number of groups.** The graph shows the Bayesian information criterion (BIC) for increasing values of K (1 to 20) which corresponds to the number of clusters to choose. The lowest value of the graph was not clear. To choose K, which describes our dataset well, we ran successive Ks (3, 4, 5, 6) to run DAPC. Finally K = 5 was chosen because it was the best assignment for each sample.
(TIFF)

**S3 Fig. DAPC representation of successive Ks.** The method we used to choose the appropriate value of K, which explained our dataset well. The group subdivision obtained was clear when we moved from K = 3 to K = 4 and K = 5. But when K = 6 the groups were not clearly delimited in the two datasets. The two datasets showed approximately the same thing.
(TIFF)
**S4 Fig.** *Number of group assumed by STRUCTURE and delta K (ΔK) value.* The result shows that the maximum ΔK value is found with K = 5. This result confirms that K = 5 is the best and appropriate number of groups to be returned.
(PDF)

**S5 Fig.** *Accessions linked based on DAPC analysis with SSR data.* The diagram shows probable membership of the five groups (K1, K2, K3, K4 and K5) determined by the DAPC analysis. An individual is represented as a vertical bar. The colors correspond to the five groups (red: group one, yellow: group two, green: group three, blue: group four, and pink: group five).
(TIFF)

**S6 Fig.** *Principal component analysis to distinguish the villages according to climatic conditions based on BIOCLIM data values.* The sampled village are arranged following the first axis according to a climatic gradient from the humid south (DANYI PEYEYEME) to the arid tropical north (GABONBONG) of Togo, followed by the Saharan climat of Senegal (DAGANA).
(TIF)

**S1 Table.** *Passport data of the sweet potato samples used on the study.* Field code, extraction code, genetic group, local name, village, latitude, longitude and country are shown.
(PDF)

**S2 Table.** *PCR conditions, number of alleles per locus on the two data analyze.* The table reports the number of alleles per locus as recorded according to the two data analyze with Gene-Mapper software. Slight difference was observed on the total number of alleles: 90 for the first dataset and 83 for the second. But the similar result was funded for all test running with the two dataset.
(PDF)

**S3 Table.** *Number of alleles in 132 samples by SSRs locus in different countries.* America, Oceania and West Africa with the two datasets. The table shows the SSRs markers used, the numbers of alleles of each markers, the total number of allele found on each continent and their means. (i) Initial dataset (s) second dataset (1) represents the number of alleles on the 132 samples chosen randomly on data base of Roullier *et al* (2011).
(PDF)

**S4 Table.** *Diversity among countries.* The table shows the p-value obtained with Wilcoxon-Paired test between different dataset. The difference is significant when p-value < 0.05.
(PDF)

**S5 Table.** *Comparison of DAPC and STRUCTURE result for k = 5.* The table reports the membership probability according to structure analyse of each individual and genetic group appurtenance as defined by DAPC results.
(PDF)

**S6 Table.** *Kruskall-Wallis test between genetic group and agro-morphological traits.* A total of 21 agro-morphological traits were used. The two datasets (initial dataset and second dataset) were used. The value corresponds to p-value of each traits. In grey significant traits p-value < 0.0024 (Bonferroni corrected p-value). The k-mean cluster, K from 3 to 5 for each dataset was tested. The total of significant traits of each model for the two datasets was mentioned.
(PDF)
**S7 Table.** Individuals found in the same group when the two datasets were run. This table shows that the two datasets produce approximately the same results: 92% of the samples are allocated to same group in the two datasets. The color or the letter for the same column of the same K for the two datasets was not identical, this show that the samples have changed cluster. Using this method, only 12 individuals have changed their place for K = 3, 6 for K = 4 and 16 for K = 5.

(PDF)

**S8 Table.** Wilcoxon-paired test between the genetic groups. A p-value < 5% was considered to be significant. The trait is specific to one group when the p-value is significant in each paired comparison with the same trait.

(PDF)

## Acknowledgments

We thank C. Roullier and Professor K. Akpagana for early on advice on this study, V. Lebot for helpful comments and the farmers who willingly participated in the study.

## Author Contributions

**Conceptualization:** KG AA NAK AB YV.

**Methodology:** KG DB MC LZ NS AB.

**Software:** KG AA NAK NS AB YV.

**Writing – original draft:** KG AA NAK AB YV.

## References

1. Burke MB, Lobell BD, Guarino L. Shifts in African crop climate by 2050, and the implications for crop improvement and genetic resources conservation. Global. Environment. Change 2009; 19: 317–325. https://doi.org/10.1016/j.gloenvcha.2009.04.003

2. Lebot V. Coping with insularity: the need for crop genetic improvement to strengthen adaptation to climatic change and food security in pacific. Environment Development and Sustainability 2013. https://doi.org/10.1007/s10668-013-9445-1

3. Mitra S. Nutritional status of orange-fleshed sweet potato's in alleviating vitamin A malnutrition through a food-based approach. Journal of Nutrition and Food Sciences 2012; 2: 160.

4. Bovel-Benjamin AC. Sweet potato: a review of its past, present and future role in human nutrition. Advances in Food and Nutrition Research. 2007; 57: 1–59.

5. CIP (International Potato Center). Annual Sub-project Progress Report. Biofortification of Sweetpotato 2005. http://www.cipotato.org/icers/3303.doc.

6. FAOSTAT (Food and Agriculture Statistic). Agricultural production crop primary database. Food and agricultural Organization of the United Nation, Rome 2014. http://faostat.fao.org/ accessed in April 2015 at 10 a.m

7. Champagne A, Legendre L, Malapa R, Lebot V. Chemotype profiling to guide breeder and explore traditional selection of root crops in Vanuatu, South Pacific. Journal of food composition and analysis. 2009; 23: 763–771.

8. Austin DF. The taxonomy, evolution and genetic diversity of sweet potato and related wilds species. In: Exploration, maintenance, and utilization of sweet potato genetic resources -Report of the First Sweet Potato Planning Conference, International Potato Center (CIP) 1987; Lima, Peru. pp 27–58.

9. Roullier C, Kambouo R, Paofa J, Mckey DB, Lebot V. On the origin of sweet potato (*Ipomoea batatas* (L.) Lam) genetic diversity in New Guinea, a secondary center of diversity. Heredity 2013; 110 (6): 94–604.

10. Lebot V. Sweet potato. Root and tuber crops: Handbook of plant breeding. J. E. Springer, in Bradshaw, ed., New York. 2010. pp. 97–125
11. Zhang DP, Carbajulca D, Ojeda L, Rossel G, Milla S, Herrera C et al. Microsatellite analysis of genetic diversity in sweet potato varieties from Latin America. CIP program report 1999–2000. International Potato Center, Lima, Peru, 2000. pp 295–301.

12. Zhang D, Rossel G, Kriegner A, Hijmans R. AFLP assessment of diversity in sweet potato from Latin America and the pacific region. Its implication on the dispersal of the crop. Genetic Resources Crop Evolution 2004; 51: 115–120.

13. Huamàn Z, Zhang DP. Sweet potato. In: biodiversity in trust. Cambridge University Press, Cambridge UK. 1997. pp. 29–38.

14. Tairo F, Mneney E, and Kullaya A. Morphological and agronomical characterization of Sweet potato [*Ipomoea batatas* (L.) Lam.] germplasm collection from Tanzania. African Journal of Plant Science. August 2008; Vol. 2 (8), 077–085.

15. Gichuru V, Aritua V, Lubega GW, Edema R, Adipala E, Rubaihayo PR A. Preliminary analyse of diversity among East African sweet potato landraces using morphological and Simple Sequence Repeats (SSRs) markers. Acta Horticulturae 2006; 703: 159–164.

16. Koussao S, Gracen V, Assante I, Danquah EY, Yeremy T, Ouedraogo TJ B et al. Diversity analysis of sweet potato (*Ipomoea batatas* (L.) Lam) gerplasm from Burkina-Faso using morphological and simple sequence repeat markers. African Journal of Biotechnology 2014; 13(6): 729–742.

17. Price TD, Qvarnstrom A and Irwin DE. The role of phenotypic plasticity in driving genetic evolution. Proceeding of the Royal Society of the London. Serie B. Biological Sciences 2003; 270:1433–1440.

18. Yada B, Tukamuhabwa P, Wajala B, Kim DJ, Skiton RA, Alajo A et al. Characterization of Uganda sweet potato germplasm using fluorescence repeat markers. Horticultural Sciences. 2010b; 45(2): 225–230.

19. Li PWY, Sun X, Han J. Using microsatellite (SSR) and morphological markers to assess the genetic diversity of 12 falcata (*Medicagosativaspp. falcata*) populations from Eurasia. Africa. Journal. Biotechnology. 2009; 8(10): 2102–2108.

20. Karuri HW, Ateka EM, Amata R, Nyende AB, Muigai AWT, Mwasame E et al. Evaluating diversity among Kenyan sweet potato genotypes using morphological and SSR markers. International Journal of Agriculture and. Biology. 2010; 12: 33–38.

21. Lebot V, Herail C, Gunua T, Pardales JR, Prana MS, Thongjiem M et al. Isozyme and RAPD variation in *Phytophthora colocasiae* Raciborski isolates from South East Asia and Oceania. Plant Pathology 2003; 52: 303–313.

22. Bioversity International Rome (Italy). International Potato Center (CIP). 2010; 11p

23. Doyle JJ and Doyle JL. Rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bulletin 1987; 19:11–15.

24. Roullier C, Rossel G, Tay D, MCkey D, Lebot V. Combining chloroplast and nuclear microsatellites to investing origin and dispersal of New World Sweet potato landrace. Molecular Ecology 2011; 20: 3963–3977. https://doi.org/10.1111/j.1365-294X.2011.05229.x PMID: 21880085

25. Roullier C, Benoit L, McKey D and Lebot V. Historical collections reveal patterns of diffusion of sweet potato in Oceania obscured by modern plant movements and recombination. PNAS, Proceedings of the National Academy of Sciences USA. 2012. https://doi.org/10.1073/pnas.1211049110.

26. Lui K, Muse SV. Power marker: Integrated analysis environment for genetic marker data. Bioinformatics 2005; 21(9): 2128–2129. https://doi.org/10.1093/bioinformatics/bti282 PMID: 15705655

27. Letunic I, Bork P. Interactive tree of life v2: online annotation and display of philogenetic trees made easy. Nucleic Acids Research 2011; 39 W475. https://doi.org/10.1093/nar/gkr201 PMID: 21470960

28. Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: A new method of principal for the analyse of genetically structured populations. BMC Genetics. 2010; 11:94, 15p. https://doi.org/10.1186/1471-2156-11-94 PMID: 20950446

29. R Development Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2015; ISBN 3-900051-07-0, URL: http://www.R-project.org.

30. Evanno G, Regnaud S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulating study. 2005; Molecular Ecology 14: 2611–2620. https://doi.org/10.1111/j.1365-294X.2005.02553.x PMID: 15969739

31. Faush D, Stephens M, Pritchard JK. Inference of population structure: Extension to linked loci and correlated allele frequencies. Genetics, 164: 1567–1587 PMID: 12930761

32. Hijmans RJ, Cameron SE, Parra JL., Jones PG. and Jarvis A,. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 2005; 25: 1965–1978
33. Tumwegamire S, Rubaihayo PR, LaBoonte DR, Diaz F, Kapinga R, Mwanga ROM, et al. Genetic diversity in white and Orange-fleshed Sweet potato Farmer varieties from East Africa evaluated by Simple Sequence Repeat Markers. Crop Science 2011; 51. 1132–1142.

34. Gwandu C, Tairo F, Mneney E, Kullaya A. Characterization of Tanzania elite sweet potato genotypes for sweet potato (*Ipomoea batatas* L..) virus disease (SPVD) resistance and high dry matter content using simple sequence repeat (SSR). African Journal of Biotechnology 2012; 11 (40): 9582–9590.

35. Vigouroux Y, Jaqueth JS, Matsuoka Y, Smith OS, Beavis WD, Smith JS, et al. Rate and pattern of mutation at microsatellite loci in maize. Molecular Biology of Evolution 2002. 19: 1251–1260

36. Rival L, Mckey D. Domestication and diversity in manioc (ManihotesculentaCrantz ssp. Esculenta, Euphorbiaceae).Current anthropology. 2008; 49 (6): 1119–1128.

37. Maquia I, Muocha I, Naico A, Martins N, Gouvea M, Andrade I et al. Molecular, morphological and agro-nomic characterization of sweet potato (*Ipomoea batatas* L.) germplasm collection from Mozambique: genotype selection for drought prone regions. South African Journal of Botany. 2013;142–151.

38. Rodriguez-Bonilla L, Cueva HE, Montero RM., Bird-Pico F, Luciano RD, Dimuth S. Assessment of genetic diversity of sweet potato in Puerto Rico. Plos ONE 2014. 9(12): 1–14.

39. Kosh-Komba E. Gestion paysanne, Diversité agro-morphologique et génétique de Manihot esculentacrantz cultivée en République Centrafricaine. Université de Lomé, Thèse 2013; 133p

40. Shigeta M. Creating landrace diversity: the case of the Ari people and, *Ensete ventricosum*, in Ethiopia. In: Ellen R., Fukui K. (Eds.), Redefining Nature: Ecology, Culture and Domestication. Oxford University Press, Oxford, 1996; 233–268.

41. Elias M, Rival L, Mckey D. Perception and management of cassava (*Manihot esculenta crantz*) diversity among Makushi Amerindian of Guyana (South America) Journal of Ethnobiology 2000; 20(2):239–265.

42. Pagnotta MA, Mondini I, Codianni P, Fares C. Agronomical quality, and molecular characterization of twenty Italian emmer wheat (Triticumd coccon) accessions. Genetic Resource Crop Evolution. 2009; 56:299–310