Enhancing Erucic Acid and Wax Ester Production in Brassica carinata through Metabolic Engineering for Industrial Applications.

Metabolic engineering enables oilseed crops to be more competitive by having more attractive properties for oleochemical industrial applications. The aim of this study was to increase the erucic acid level and to produce wax ester (WE) in seed oil by genetic transformation to enhance the industrial applications of B. carinata. Six transgenic lines for high erucic acid and fifteen transgenic lines for wax esters were obtained. The integration of the target genes for high erucic acid (BnFAE1 and LdPLAAT) and for WEs (ScWS and ScFAR) in the genome of B. carinata cv. 'Derash' was confirmed by PCR analysis. The qRT-PCR results showed overexpression of BnFAE1 and LdPLAAT and downregulation of RNAi-BcFAD2 in the seeds of the transgenic lines. The fatty acid profile and WE content and profile in the seed oil of the transgenic lines and wild type grown in biotron were analyzed using gas chromatography and nanoelectrospray coupled with tandem mass spectrometry. A significant increase in erucic acid was observed in some transgenic lines ranging from 19% to 29% in relation to the wild type, with a level of erucic acid reaching up to 52.7%. Likewise, the transgenic lines harboring ScFAR and ScWS genes produced up to 25% WE content, and the most abundant WE species were 22:1/20:1 and 22:1/22:1. This study demonstrated that metabolic engineering is an effective biotechnological approach for developing B. carinata into an industrial crop.

Affordable and robust phenotyping framework to analyse root system architecture of soil-grown plants.

The phenotypic analysis of root system growth is important to inform efforts to enhance plant resource acquisition from soils; however, root phenotyping remains challenging because of the opacity of soil, requiring systems that facilitate root system visibility and image acquisition. Previously reported systems require costly or bespoke materials not available in most countries, where breeders need tools to select varieties best adapted to local soils and field conditions. Here, we report an affordable soil-based growth (rhizobox) and imaging system to phenotype root development in glasshouses or shelters. All components of the system are made from locally available commodity components, facilitating the adoption of this affordable technology in low-income countries. The rhizobox is large enough (approximately 6000 cm2 of visible soil) to avoid restricting vertical root system growth for most if not all of the life cycle, yet light enough (approximately 21 kg when filled with soil) for routine handling. Support structures and an imaging station, with five cameras covering the whole soil surface, complement the rhizoboxes. Images are acquired via the Phenotiki sensor interface, collected, stitched and analysed. Root system architecture (RSA) parameters are quantified without intervention. The RSAs of a dicot species (Cicer arietinum, chickpea) and a monocot species (Hordeum vulgare, barley), exhibiting contrasting root systems, were analysed. Insights into root system dynamics during vegetative and reproductive stages of the chickpea life cycle were obtained. This affordable system is relevant for efforts in Ethiopia and other low- and middle-income countries to enhance crop yields and climate resilience sustainably.

Genetic diversity and population structure of Zymoseptoria tritici in Ethiopia as revealed by microsatellite markers.

Septoria tritici blotch (STB), caused by Zymoseptoria tritici (formerly: Mycosphaerella graminicola or Septoria tritici), is one of the most devastating diseases of wheat globally. Understanding genetic diversity of the pathogen has supreme importance in developing best management strategies. However, there is dearth of information on the genetic structure of Z. tritici populations in Ethiopia. Therefore, the present study was targeted to uncover the genetic diversity and population structure of Z. tritici populations from the major wheat-growing areas of Ethiopia. Totally, 182 Z. tritici isolates representing eight populations were analyzed with 14 microsatellite markers. All the microsatellite loci were polymorphic and highly informative, and hence useful genetic tools to depict the genetic diversity and population structure of the pathogen. A wide range of diversity indices including number of observed alleles, effective number of alleles, Shannon's diversity index, number of private alleles, Nei's gene diversity and percentage of polymorphic loci (PPL) were computed to determine genetic variation within populations. A high within-populations genetic diversity was confirmed with gene diversity index and PPL values ranging from 0.34 - 0.58 and 79-100% with overall mean of 0.45 and 94%, respectively. Analysis of molecular variance (AMOVA) revealed a moderate genetic differentiation where 92% of the total genetic variation resides within populations, leaving only 8% among populations. Cluster (UPGMA), PCoA and STRUCTURE analyses did not group the populations into sharply genetically distinct clusters according to their geographical origins, likely due to high gene flow (Nm = 5.66) and reproductive biology of the pathogen. All individual samples shared alleles from two subgroups (K = 2) evidencing high potential of genetic admixture. In conclusion, the microsatellite markers used in the present study were highly informative and thus, helped to dissect the genetic structures of Z. tritici populations in Ethiopia. Among the studied populations, those of East Shewa, Arsi, South West Shewa and Bale showed a high genetic diversity, and hence these areas can be considered as hot spots for investigations planned on the pathogen and host-pathogen interactions. Therefore, the present study not only enriches missing information in Ethiopia but also provides new insights into the epidemiology and genetic structure of Z. tritici in Africa where the agro-climatic conditions and the wheat cropping systems are different from other parts of the world. Such baseline information is useful for designing and implementing durable and effective management strategies.

Genetic Diversity of Durum Wheat (Triticum turgidum L. ssp. durum, Desf) Germplasm as Revealed by Morphological and SSR Markers.

Ethiopia is considered a center of origin and diversity for durum wheat and is endowed with many diverse landraces. This research aimed to estimate the extent and pattern of genetic diversity in Ethiopian durum wheat germplasm. Thus, 104 durum wheat genotypes representing thirteen populations, three regions, and four altitudinal classes were investigated for their genetic diversity, using 10 grain quality- and grain yield-related phenotypic traits and 14 simple sequence repeat (SSR) makers. The analysis of the phenotypic traits revealed a high mean Shannon diversity index (H' = 0.78) among the genotypes and indicated a high level of phenotypic variation. The principal component analysis (PCA) classified the genotypes into three groups. The SSR markers showed a high mean value of polymorphic information content (PIC = 0.50) and gene diversity (h = 0.56), and a moderate number of alleles per locus (Na = 4). Analysis of molecular variance (AMOVA) revealed a high level of variation within populations, regions, and altitudinal classes, accounting for 88%, 97%, and 97% of the total variation, respectively. Pairwise genetic differentiation and Nei's genetic distance analyses identified that the cultivars are distinct from the landrace populations. The distance-based (Discriminant Analysis of Principal Component (DAPC) and Minimum Spanning Network (MSN)) and model-based population stratification (STRUCTURE) methods of clustering grouped the genotypes into two clusters. Both the phenotypic data-based PCA and the molecular data-based DAPC and MSN analyses defined distinct groupings of cultivars and landraces. The phenotypic and molecular diversity analyses highlighted the high genetic variation in the Ethiopian durum wheat gene pool. The investigated SSRs showed significant associations with one or more target phenotypic traits. The markers identify landraces with high grain yield and quality traits. This study highlights the usefulness of Ethiopian landraces for cultivar development, contributing to food security in the region and beyond.

Genetic architecture of adult-plant resistance to stripe rust in bread wheat (Triticum aestivum L.) association panel.

Stripe rust, caused by Puccinia striiformis f. sp. tritici, is a severe disease in wheat worldwide, including Ethiopia, causing up to 100% wheat yield loss in the worst season. The use of resistant cultivars is considered to be the most effective and durable management technique for controlling the disease. Therefore, the present study targeted the genetic architecture of adult plant resistance to yellow rust in 178 wheat association panels. The panel was phenotyped for yellow rust adult-plant resistance at three locations. Phonological, yield, yield-related, and agro-morphological traits were recorded. The association panel was fingerprinted using the genotyping-by-sequencing (GBS) platform, and a total of 6,788 polymorphic single nucleotide polymorphisms (SNPs) were used for genome-wide association analysis to identify effective yellow rust resistance genes. The marker-trait association analysis was conducted using the Genome Association and Prediction Integrated Tool (GAPIT). The broad-sense heritability for the considered traits ranged from 74.52% to 88.64%, implying the presence of promising yellow rust resistance alleles in the association panel that could be deployed to improve wheat resistance to the disease. The overall linkage disequilibrium (LD) declined within an average physical distance of 31.44 Mbp at r2 = 0.2. Marker-trait association (MTA) analysis identified 148 loci significantly (p = 0.001) associated with yellow rust adult-plant resistance. Most of the detected resistance quantitative trait loci (QTLs) were located on the same chromosomes as previously reported QTLs for yellow rust resistance and mapped on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 3D, 4A, 4B, 4D, 5A, 5B, 6A, 6B, 7A, and 7D. However, 12 of the discovered MTAs were not previously documented in the wheat literature, suggesting that they could represent novel loci for stripe rust resistance. Zooming into the QTL regions in IWGSC RefSeq Annotation v1 identified crucial disease resistance-associated genes that are key in plants' defense mechanisms against pathogen infections. The detected QTLs will be helpful for marker-assisted breeding of wheat to increase resistance to stripe rust. Generally, the present study identified putative QTLs for field resistance to yellow rust and some important agronomic traits. Most of the discovered QTLs have been reported previously, indicating the potential to improve wheat resistance to yellow rust by deploying the QTLs discovered by marker-assisted selection.

In Vitro Regeneration of Korarima (Aframomum corrorima (Braun) P. C. M. Jansen): A Threatened Spice and Medicinal Herb from Ethiopia.

Developing an in vitro regeneration system is very important to increase production and productivity of plants as well as for the conservation of rare and threatened medicinal plants like korarima (Aframomum corrorima (Braun) P. C. M. Jansen). To date, no study dealing with in vitro indirect regeneration system of korarima has been reported. Thus, in this study, we developed an efficient and reproducible protocol for in vitro regeneration of korarima via callus. The procedure involved soaking seeds in 50% H2SO4 for 16 h that resulted in 92.5% germination on plant growth regulators (PGRs)-free half-strength Murashige and Skoog (MS) basal medium after a month. Shoot and rhizome induction rate of 93.75% was obtained on the MS medium containing 1.5 mg/l BAP in combination with 0.1 mg/l IBA after five weeks. Whitish yellow friable callus was obtained from rhizome culture taken from in vitro grown plantlets. The MS medium containing 2.0 mg/l 2, 4D in combination with 0.5 mg/l kinetin, resulted in 77.5% callus induction. The shoot regeneration rate of 45% was obtained from callus on the MS medium containing 2.0 mg/l TDZ in combination with 0.5 mg/l IBA. The mean shoot number of 10.83 per explant was obtained upon multiplication on the MS medium containing 1.5 mg/l BAP with a mean shoot height of 5.37 cm. The best rooting responses were obtained on half MS medium supplemented with 0.5 mg/l IAA resulting in a mean number of root of 18.59, mean root length of 9.71 cm, and mean shoot height of 7.32 cm. The plantlets showed 75% survival efficiency after acclimatization. The present regeneration protocol offers a conceivable system towards effective conservation and genetic improvement of the crop by increasing the efficiency of genetic transformation.

Novel GBS-Based SNP Markers for Finger Millet and Their Use in Genetic Diversity Analyses.

Eleusine coracana (L.) Gaertn., commonly known as finger millet, is a multipurpose crop used for food and feed. Genomic tools are required for the characterization of crop gene pools and their genomics-led breeding. High-throughput sequencing-based characterization of finger millet germplasm representing diverse agro-ecologies was considered an effective method for determining its genetic diversity, thereby suggesting potential candidates for breeding. In this study, the genotyping-by-sequencing (GBS) method was used to simultaneously identify novel single nucleotide polymorphism (SNP) markers and genotype 288 finger millet accessions collected from Ethiopia and Zimbabwe. The accessions were characterized at individual and group levels using 5,226 bi-allelic SNPs, with a minimum allele frequency (MAF) of above 0.05, distributed across 2,500 scaffolds of the finger millet reference genome. The polymorphism information content (PIC) of the SNPs was 0.23 on average, and a quarter of them have PIC values over 0.32, making them highly informative. The grouping of the 288 accessions into seven populations based on geographic proximity and the potential for germplasm exchange revealed a narrow range of observed heterozygosity (Ho; 0.09-0.11) and expected heterozygosity (He) that ranged over twofold, from 0.11 to 0.26. Alleles unique to the different groups were also identified, which merit further investigation for their potential association with desirable traits. The analysis of molecular variance (AMOVA) revealed a highly significant genetic differentiation among groups of accessions classified based on the geographic region, country of origin, days to flowering, panicle type, and Al tolerance (p < 0.01). The high genetic differentiation between Ethiopian and Zimbabwean accessions was evident in the AMOVA, cluster, principal coordinate, and population structure analyses. The level of genetic diversity of finger millet accessions varies moderately among locations within Ethiopia, with accessions from the northern region having the lowest level. In the neighbor-joining cluster analysis, most of the improved cultivars included in this study were closely clustered, probably because they were developed using genetically less diverse germplasm and/or selected for similar traits, such as grain yield. The recombination of alleles via crossbreeding genetically distinct accessions from different regions of the two countries can potentially lead to the development of superior cultivars.

Finger millet RNA-seq reveals differential gene expression associated with tolerance to aluminum toxicity and provides novel genomic resources.

Eleusine coracana, finger millet, is a multipurpose crop cultivated in arid and semi-arid regions of Africa and Asia. RNA sequencing (RNA-seq) was used in this study to obtain valuable genomic resources and identify genes differentially expressed between Al-tolerant and Al-susceptible genotypes. Two groups of finger millet genotypes were used: Al-tolerant (215836, 215845, and 229722) and Al-susceptible (212462, 215804 and 238323). The analysis of the RNA-seq data resulted in 198,546 unigenes, 56.5% of which were annotated with significant hits in one or more of the following six databases: NR (48.8%), GO (29.7%), KEGG (45%), PlantTFDB (19.0%), Uniprot (49.2%), and NT (46.2%). It is noteworthy that only 220 unigenes in the NR database had significant hits against finger millet sequences suggesting that finger millet's genomic resources are scarce. The gene expression analysis revealed that 322 genes were significantly differentially expressed between the Al-tolerant and Al-susceptible genotypes, of which 40.7% were upregulated while 59.3% were downregulated in Al-tolerant genotypes. Among the significant DEGs, 54.7% were annotated in the GO database with the top hits being ATP binding (GO:0005524) and DNA binding (GO:0003677) in the molecular function, DNA integration (GO:0015074) and cell redox homeostasis in the biological process, as well as cellular anatomical entity and intracellular component in the cellular component GO classes. Several of the annotated DEGs were significantly enriched for their corresponding GO terms. The KEGG pathway analysis resulted in 60 DEGs that were annotated with different pathway classes, of which carbohydrate metabolism and signal transduction were the most prominent. The homologs of a number of significant DEGs have been previously reported as being associated with Al or other abiotic stress responses in various crops, including carboxypeptidase SOL1, HMA3, AP2, bZIP, C3H, and WRKY TF genes. A more detailed investigation of these and other DEGs will enable genomic-led breeding for Al tolerance in finger millet. RNA-seq data analysis also yielded 119,073 SNP markers, the majority of which had PIC values above 0.3, indicating that they are highly informative. Additionally, 3,553 single-copy SSR markers were identified, of which trinucleotide SSRs were the most prevalent. These genomic resources contribute substantially to the enrichment of genomic databases for finger millet, and facilitate future research on this crop.

Genetic variability and population structure of Ethiopian chickpea (Cicer arietinum L.) germplasm.

Evaluation of the genetic diversity and an understanding of the genetic structure and relationships of chickpea genotypes are valuable to design efficient germplasm conservation strategies and crop breeding programs. Information is limited, in these regards, for Ethiopian chickpea germplasms. Therefore, the present study was carried out to estimate the genetic diversity, population structure, and relationships of 152 chickpea genotypes using simple sequence repeats (SSR) markers. Twenty three SSR markers exhibited polymorphism producing a total of 133 alleles, with a mean of 5.8 alleles per locus. Analyses utilizing various genetic-based statistics included pairwise population Nei's genetic distance, heterozygosity, Shannon's information index, polymorphic information content, and percent polymorphism. These analyses exemplified the existence of high genetic variation within and among chickpea genotypes. The 152 genotypes were divided into two major clusters based on Nei's genetic distances. The exotic genotypes were grouped in one cluster exclusively showing that these genotypes are distinct to Ethiopian genotypes, while the patterns of clustering of Ethiopian chickpea genotypes based on their geographic region were not consistent because of the seed exchange across regions. Model-based population structure clustering identified two discrete populations. These finding provides useful insight for chickpea collections and ex-situ conservation and national breeding programs for widening the genetic base of chickpea.

Genome-wide association mapping identifies yellow rust resistance loci in Ethiopian durum wheat germplasm.

Durum wheat is an important cereal grown in Ethiopia, a country which is also its center for genetic diversity. Yellow (stripe) rust caused by Puccinia striiformis fsp tritici is one of the most devastating diseases threatening Ethiopian wheat production. To identify sources of genetic resistance and combat this pathogen, we conducted a genome wide association study of yellow rust resistance on 300 durum wheat accessions comprising 261 landraces and 39 cultivars. The accessions were evaluated for their field resistance using a modified Cobb scale at Meraro, Kulumsa and Chefe Donsa in the 2015 and 2016 main growing seasons. Analysis of the 35K Axiom Array genotyping data of the panel resulted in a total of 8,797 polymorphic SNPs of which 7,093 were used in subsequent analyses. Population structure analysis suggested two groups in which the cultivars clearly stood out separately from the landraces. Eleven SNPs significantly associated with yellow rust resistance were identified on four chromosomes (1A, 1B, 2B, and 5A) which defined at least five genomic loci. Six of the SNPs were consistently identified on chromosome 1B singly at each and combined overall environments which explained 62.6-64.0% of the phenotypic variation (R2). Resistant allele frequency ranged from 14.0-71.0%; Zooming in to the identified resistance loci revealed the presence of disease resistance related genes involved in the plant defense system such as the ABC transporter gene family, disease resistance protein RPM1 (NBS-LRR class), Receptor kinases and Protein kinases. This study has provided SNPs for tracking the loci associated with yellow rust resistance and a diversity panel which can be used for association study of other agriculturally important traits in durum wheat.

Novel Expressed Sequence Tag-Derived and Other Genomic Simple Sequence Repeat Markers Revealed Genetic Diversity in Ethiopian Finger Millet Landrace Populations and Cultivars.

Finger millet (Eleusine coracana (L.) Geartn.) is a self-pollinating amphidiploid crop cultivated with minimal input for food and feed, as well as a source of income for small-scale farmers. To efficiently assess its genetic diversity for conservation and use in breeding programs, polymorphic DNA markers that represent its complex tetraploid genome have to be developed and used. In this study, 13 new expressed sequence tag-derived simple sequence repeat (EST-SSR) markers were developed based on publicly available finger millet ESTs. Using 10 polymorphic SSR markers (3 genomic and 7 novel EST-derived), the genetic diversity of 55 landrace accessions and 5 cultivars of finger millet representing its major growing areas in Ethiopia was assessed. In total, 26 alleles were detected across the 10 loci, and the average observed number of alleles per locus was 5.6. The polymorphic information content (PIC) of the loci ranged from 0.045 (Elco-48) to 0.71 (UGEP-66). The level of genetic diversity did not differ much between the accessions with the mean gene diversity estimates ranging only from 0.44 (accession 216054) to 0.68 (accession 237443). Similarly, a narrow range of variation was recorded at the level of regional states ranging from 0.54 (Oromia) to 0.59 (Amhara and Tigray). Interestingly, the average gene diversity of the landrace accessions (0.57) was similar to that of the cultivars (0.58). The analysis of molecular variance (AMOVA) revealed significant genetic variation both within and among accessions. The variation among the accessions accounted for 18.8% of the total variation (F ST = 0.19; P < 0.001). Similarly, significant genetic variation was obtained among the geographic regions, accounting for 6.9% of the total variation (P < 0.001). The results of the cluster, principal coordinate, and population structure analyses suggest a poor correlation between the genetic makeups of finger millet landrace populations and their geographic regions of origin, which in turn suggests strong gene flow between populations within and across geographic regions. This study contributed novel EST-SSR markers for their various applications, and those that were monomorphic should be tested in more diverse finger millet genetic resources.

Genome-Wide Association Study Reveals Novel Genetic Loci for Quantitative Resistance to Septoria Tritici Blotch in Wheat (Triticum aestivum L.).

Septoria tritici blotch, caused by the fungus Zymoseptoria titici, poses serious and persistent challenges to wheat cultivation in Ethiopia and worldwide. Deploying resistant cultivars is a major component of controlling septoria tritici blotch (STB). Thus, the objective of this study was to elucidate the genomic architecture of STB resistance in an association panel of 178 bread wheat genotypes. The association panel was phenotyped for STB resistance, phenology, yield, and yield-related traits in three locations for 2 years. The panel was also genotyped for single nucleotide polymorphism (SNP) markers using the genotyping-by-sequencing (GBS) method, and a total of 7,776 polymorphic SNPs were used in the subsequent analyses. Marker-trait associations were also computed using a genome association and prediction integrated tool (GAPIT). The study then found that the broad-sense heritability for STB resistance ranged from 0.58 to 0.97 and 0.72 to 0.81 at the individual and across-environment levels, respectively, indicating the presence of STB resistance alleles in the association panel. Population structure and principal component analyses detected two sub-groups with greater degrees of admixture. A linkage disequilibrium (LD) analysis in 338,125 marker pairs also detected the existence of significant (p ≤ 0.01) linkage in 27.6% of the marker pairs. Specifically, in all chromosomes, the LD between SNPs declined within 2.26-105.62 Mbp, with an overall mean of 31.44 Mbp. Furthermore, the association analysis identified 53 loci that were significantly (false discovery rate, FDR, <0.05) associated with STB resistance, further pointing to 33 putative quantitative trait loci (QTLs). Most of these shared similar chromosomes with already published Septoria resistance genes, which were distributed across chromosomes 1B, 1D, 2A, 2B, 2D, 3A,3 B, 3D, 4A, 5A, 5B, 6A, 7A, 7B, and 7D. However, five of the putative QTLs identified on chromosomes 1A, 5D, and 6B appeared to be novel. Dissecting the detected loci on IWGSC RefSeq Annotation v2.1 revealed the existence of disease resistance-associated genes in the identified QTL regions that are involved in plant defense responses. These putative QTLs explained 2.7-13.2% of the total phenotypic variation. Seven of the QTLs (R 2 = 2.7-10.8%) for STB resistance also co-localized with marker-trait associations (MTAs) for agronomic traits. Overall, this analysis reported on putative QTLs for adult plant resistance to STB and some important agronomic traits. The reported and novel QTLs have been identified previously, indicating the potential to improve STB resistance by pyramiding QTLs by marker-assisted selection.

Genetic diversity analysis in Plectranthus edulis (Vatke) Agnew populations collected from diverse geographic regions in Ethiopia using inter-simple sequence repeats (ISSRs) DNA marker system.

BACKGROUND: Plectranthus edulis (Vatke) Agnew (Lamiaceae), locally known as Ethiopian potato syno. Ethiopian dinich, is one of the native Ethiopian edible tuber crops that has been significantly contributing to household food security for millions of subsistence farmers. However, its current production is declining to the extent of total extinction from several administrative regions where it used to be widely cultivated. It is one of the less researched crops regardless of being indigenous and its contribution to food security during time of scarcity. Therefore, we intended to assess the level of genetic diversity in 67 accessions, representing nine populations that were collected from diverse agro-ecologies in the country, using ISSR markers and hence, to generate a baseline information that assists marker assisted breeding, conservation and germplasm management efforts. RESULTS: In the present study, ten polymorphic ISSR markers were screened and optimized, that generated an average of 7.4 scorable bands per marker and revealed high overall percent polymorphism (95%), Nei's gene diversity (h = 0.40) and Shannon index (I = 0.62) suggesting ISSR's effectiveness in detecting high levels of genetic diversity. A considerably high overall populations gene diversity (Nei's) (h = 0.32) and Shannon index (I = 0.47) were observed, revealing high potential of the populations for further breeding and conservation efforts particularly for population from Gurage administrative zone, which showed the highest values. Similarly, estimation of pairwise genetic distance revealed the importance of cross breeding population from Awi administrative zone to the rest populations. Analysis of hierarchical molecular variance (AMOVA) showed higher levels of genetic differentiation within populations (92%), and collection regions (94%) suggesting that either clonal mode of propagation in the crop or farmers selection pressure for important agronomic traits or both maintained the original heterozygosity in the crop. UPGMA phylogenetic analysis did not strictly group the populations based on their geographic region of origin, which could be attributed to the widely practiced tuber exchange and hence continuous human mediated exchange of genetic material and sharing of the same genetic base among the geographic regions. CONCLUSIONS: The ISSR markers used in the present study were effective in revealing extent and patterns of genetic diversity in P. edulis populations. However, it is important to couple them with agro-morphological traits or codominant molecular markers to get more reliable information for use in breeding and conservation. Several of the potential administrative zones we covered are useful for P. edulis diversification and conservation. However, the crop is currently highly marginalized and this led to rapid decline in population size and loss of valuable agronomic traits. To address this challenge, there is an urgent need to take counteractive measures.