- 1Institute of Biotechnology, Addis Ababa University, Addis Ababa, Ethiopia
- 2Department of Plant Breeding, Swedish University of Agricultural Sciences, Lomma, Sweden
- 3Ethiopian Biotechnology Institute, Addis Ababa, Ethiopia
Ethiopia is the center of origin for sorghum [Sorghum bicolor (L.) Moench], where the distinct agro-ecological zones significantly contributed to the genetic diversity of the crops. A large number of sorghum landrace accessions have been conserved ex situ. Molecular characterization of this diverse germplasm can contribute to its efficient conservation and utilization in the breeding programs. This study aimed to investigate the genetic diversity of Ethiopian sorghum using gene-based single nucleotide polymorphism (SNP) markers. In total, 359 individuals representing 24 landrace accessions were genotyped using 3,001 SNP markers. The SNP markers had moderately high polymorphism information content (PIC = 0.24) and gene diversity (H = 0.29), on average. This study revealed 48 SNP loci that were significantly deviated from Hardy–Weinberg equilibrium with excess heterozygosity and 13 loci presumed to be under selection (P < 0.01). The analysis of molecular variance (AMOVA) determined that 35.5% of the total variation occurred within and 64.5% among the accessions. Similarly, significant differentiations were observed among geographic regions and peduncle shape-based groups. In the latter case, accessions with bent peduncles had higher genetic variation than those with erect peduncles. More alleles that are private were found in the eastern region than in the other regions of the country, suggesting a good in situ conservation status in the east. Cluster, principal coordinates (PCoA), and STRUCTURE analyses revealed distinct accession clusters. Hence, crossbreeding genotypes from different clusters and evaluating their progenies for desirable traits is advantageous. The exceptionally high heterozygosity observed in accession SB4 and SB21 from the western geographic region is an intriguing finding of this study, which merits further investigation.
Introduction
Sorghum [Sorghum bicolor (L.) Moench] is the fifth most important cereal crop in the world next to maize, rice, wheat, and barley in terms of both production and harvested area (FAOSTAT, 2019). It is a major food crop for more than 500 million people across Africa, Asia, and Latin America, particularly for those in the semi-arid tropical regions (Ejeta, 2005). It is grown in drought-prone areas where several other crops cannot reliably grow. Recent FAOSTAT data on annual global production of sorghum showed that it covered about 40 million ha of land and produced grains of ca 57.9 million metric tons (MMT) (FAOSTAT, 2019). The United States, Nigeria, and Ethiopia are the leading sorghum-producing countries in the world with a total production of 8.6, 6.7, and 5.2 MMT, respectively (Statista, 2020). In Africa, sorghum is the second most widely cultivated cereal crop, only surpassed by maize (FAOSTAT, 2019).
Ethiopia is considered as one of the centers of origin and diversity of sorghum (De Wet and Harlan, 1971) due to the presence of wild relatives and diversified forms of the crop in the country. The sorghum gene pool in the country has been used as novel sources of germplasm for crop improvements. For example, genotypes harboring genes that confer resistance to ergot and green bug (Wu et al., 2006) as well as high lysine (Singh and Axtell, 1973) and drought-tolerant (Borrell et al., 2000) sorghum genotypes were identified from the Ethiopian accessions.
Studying the genetic diversity of a crop is very important for effective germplasm management, utilization, and genotype selection for crop improvement (Bucheyeki et al., 2009). It is the most important step for conserving and increasing the rate of genetic gain in crop-breeding programs. The level of genetic diversity within a species is commonly used to measure the level of species adaptability and survival in unpredictable environmental conditions (Rao and Hodgkin, 2002; Govindaraj et al., 2015). Similarly, the level of genetic variation within a population is the basis for germplasm selection in plant breeding and is vital for crop improvement (Mohammadi and Prasanna, 2003). Hence, the conservation and utilization of plant genetic variation are crucial to human food security (Rao and Hodgkin, 2002).
Sorghum is a predominantly self-pollinated diploid species (Poehlman and Sleper, 1979) with 2n = 2× = 20 chromosomes. It has a small genome relative to other cereal crops, which is about 730 Mbp (Paterson et al., 2009). Its whole genome was sequenced and made accessible for public use1 (Paterson et al., 2009; McCormick et al., 2018), which facilitated the development of DNA markers, such as single nucleotide polymorphism (SNPs) for various applications, including analyses of population genetics and identification of genomic regions associated with complex traits through quantitative trait loci (QTL) and association mapping (Too et al., 2018; Girma et al., 2019).
The genetic diversity of crop species can be studied through morphological, biochemical, and molecular markers (Rao et al., 1996; Geleta and Labuschagne, 2005; Mehmood et al., 2008; Enyew et al., 2021). Previous studies on the genetic diversity of sorghum have been carried out by using random amplified polymorphism DNA (RAPD) analysis (Ayana et al., 2000; Ruiz-Chután et al., 2019), simple sequence repeat (SSR) markers (Djè et al., 2000; Ghebru et al., 2002; Manzelli et al., 2007; Ali et al., 2008; Wang et al., 2009; Ng’uni et al., 2011, 2012; Burow et al., 2012; Adugna et al., 2013; Adugna, 2014; Mofokeng et al., 2014; Motlhaodi et al., 2014, 2017), and express sequence tags (EST) SSR markers (Ramu et al., 2013), SNP markers (Cuevas et al., 2017; Afolayan et al., 2019; Cuevas and Prom, 2020). More recently, a few studies have been performed on the genetic diversity of Ethiopian sorghum accessions using SNP markers (Girma et al., 2019; Menamo et al., 2021; Wondimu et al., 2021). These studies brought out the contribution of geographic regions and agro-ecological zones for the genetic variation and population structure of sorghum grown in Ethiopia. However, these studies did not consider genetic variation within populations, as the analyses were based on either a single plant per accession or a pool of individual plants treated as a single sample per accession). Ethiopian Biodiversity Institute (EBI) has conserved more than 9,432 sorghum accessions collected from diverse agro-ecologies across the country.2 However, the genetic diversity of most of the accessions in the collection remains molecularly uncharacterized. Therefore, this study analyzes the genetic diversity and population structure of selected Ethiopian sorghum accessions using SNP markers in order to generate highly important information, which together with previous research results, lead to deeper insight on the sorghum gene pool in the country and beyond.
Materials and Methods
Plant Materials
Twenty-four Ethiopian sorghum landrace accessions originally collected by the EBI were obtained from Melkassa Agricultural Research Center (MARC). The accessions were selected to represent three agro-ecological zones according to the classification by Amede et al. (2015) viz. cool/subhumid, cool/semiarid, and warm/semiarid zones (Supplementary Figure 1). Supplementary Table 1 provides details about these accessions, including the sampling locations, as well as major morphological and phenological characteristics. Photographs showing panicle diversity in the Ethiopian sorghum that represents these accessions are provided as Supplementary Figure 2.
Planting, Sampling, and Genomic DNA Extraction
Sorghum seeds representing the 24 accessions were planted using plastic pots filled with soil in a greenhouse at the Department of Plant Breeding, SLU, Sweden. Two weeks after planting, the leaf tissues from individual plants were collected using a sample collection kit provided by LGC-Genomics (Berlin, Germany), as described by Tsehay et al. (2020). Each accession was represented by 15 individual plants, except accession SB10, which was represented by 14 individuals; hence 359 genotypes were sampled in total. The samples were then sent to LGC Genomics (Berlin, Germany) where genomic DNA extraction was conducted for subsequent genotyping. High-quality genomic DNA, suitable for next-generation sequencing (NGS), was extracted using the Sbeadex plant kit.3
SNP Selection, Assay Design, Sequencing, and Genotype Calling
The vast majority of SNPs (97%) used in this study were selected from sorghum genome SNP database SorGSD,4 a web-portal that provides genome-wide SNP markers for diverse sorghum genetic resources (Luo et al., 2016). Among different sorghum lines in the database, Cherekit (an Ethiopian sorghum landrace accession) was targeted for selecting the SNPs. For genotyping, SeqSNP method (an advanced NGS method for genotyping target SNPs) was used. Initially, 12,316 SNPs were targeted for high-specificity (without allowing for off-target hit) assay design, using Sorghum bicolor v3.1.1 genome in Phytozome 12.15 as a reference (Paterson et al., 2009; McCormick et al., 2018). Additionally, 380 SNPs within functionally annotated sorghum genes were identified through the Basic Local Alignment Search Tool (BLAST), searching the genes targeting S. bicolor v3.1.1 genome sequence using Phytozome 12.1 search function were targeted for the assay design. Out of the total 12,696 targets used for the high specificity assay design, 9,495 were totally covered (two oligo probes per target), 1,631 were partially covered (one oligo probe), whereas 1,190 failed.
For the seqSNP genotyping, 5,000 SNPs were selected among the totally covered SNPs, based on their distribution across the sorghum genome. The number of SNPs targeted on chromosome-1 to chromosome-10 included in the order, 532, 521, 572, 506, 497, 515, 437, 446, 465, and 509 (refer to Supplementary Table 2). One hundred fifty-seven of these SNPs belong to 51 functionally annotated genes (Supplementary Table 3). This was followed by the construction of SeqSNP kit LGC, Biosearch Technologies (Berlin, Germany) comprising 10,000 high-specificity oligo probes for the 5,000 target SNPs and construction of a sequencing library. The target sequencing was conducted using Illumina NextSeq 500/550 v2 system with 75 bp single read sequencing mode. In the end, ca 973,000 reads per sample were obtained and the effective target of SNP coverage per sample was 175 times on average. After sequencing, the reads were adapter-clipped and quality-trimmed to get a minimum Phred quality score of 30 over a window of ten bases. After discarding reads shorter than 65 bases. the quality trimmed reads were aligned against the reference genome using Bowtie2 v2.2.3 (Langmead and Salzberg, 2012), and the SNP genotyping pipeline was set to diploid genotyping with a minimum coverage of eight reads per sample per locus. The variant identification and genotype calling were done using Freebayes v1.0.2-16 (Garrison and Marth, 2012).
Data Analysis
The site frequency spectra were analyzed for each accession using DnaSP version 6 (Rozas et al., 2003). The nucleotide diversity (Nei, 1987) and Tajima’s D (Tajima, 1989) were calculated using the PopGenome package (Pfeifer et al., 2014) in R software (R Core Team) to reveal the genome-wide pattern of variation using a sliding window approach (window size = 1 Mb, step size = 200 kb), in line with previous studies in sorghum (Yan et al., 2018), maize, and common bean (Lai et al., 2010; Cortés and Blair, 2018).
The mean effective number of alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (He), Shannon information index (I), and gene flow (Nm) for each SNP marker and accession were estimated using GenAlEx 6.5 (Peakall and Smouse, 2012), and the gene diversity (H) and the polymorphism information content (PIC) were performed using PowerMarker (Liu and Muse, 2005). The Hardy–Weinberg equilibrium (HWE) test was also done using GenAlEx 6.5.
Analysis of Molecular Variance (AMOVA) within and among the accessions as well as at higher hierarchical levels were done using the software, Arlequin ver. 3.5.2.2 (Excoffier and Lischer, 2010). Arlequin was also used for estimating pairwise genetic differentiation between accessions and groups and for detecting outlier SNP markers through a non-hierarchical finite island model. The significance of the differentiation of accessions and groups was tested by 10,000 permutations. The joint distribution of population differentiation (FST) and heterozygosity (heterozygosity within populations)/(1 – FST) were obtained according to Excoffier and Lischer (2010). The loci under selection were identified based on the FST significance level of P < 0.01.
The principal coordinates analysis (PCoA) was done using GenAlEx 6.5. The bootstrap-supported unweighted pair group method with arithmetic mean (UPGMA) clustering based on Nei’s genetic distance (Nei and Takezaki, 1983) was performed using PowerMarker v 3.25 (Liu and Muse, 2005) and the resulting trees were visualized using MEGA-X (Kumar et al., 2018). STRUCTURE v. 2.3.4 software (Pritchard et al., 2000) was used for the Bayesian clustering of the 359 individuals representing the 24 sorghum accessions, at the burn-in period length of 100,000 and a Markov Chain Monte Carlo (MCMC) replications of 100,000. The structure analysis was done for K ranging from two to ten, with ten iterations at each K, to determine the optimum number of clusters (genetic populations). The optimum K value was predicted following the simulation method of Evanno et al. (2005) using STRUCTURE HARVESTER version 0.6.92 (Earl, 2012). A bar plot for the optimum K was determined through Clumpak beta version (Kopelman et al., 2015).
Results
The Quality and Level of Polymorphism of SNP Markers
In the data matrix of 5,000 SNP loci for the 359 sorghum genotypes, missing data accounted for 2.8% and, of those with data, 94.1% were homozygous. Among the 5,000 target SNP loci, 4,301 (86%) were polymorphic, whereas 699 loci (14%) were monomorphic across the 359 genotypes. Among the 4,301 SNP loci, 4,256, 42, and 1 were bi-, tri- and tetra-allelic, respectively, whereas two loci had the combination of SNP and length variants. Tri- and tetra-allelic loci were excluded from further analysis. The filtering of the 4,256 bi-allelic SNP data, based on the percentage of missing data and minimum allele frequency (MAF) resulted in different numbers of loci varying from 2,259 loci with less than1% missing and greater than10% MAF to 3,089 loci with less than 5% missing and greater than 5% MAF. For the population genetics analyses, 3,001 bi-allelic SNP loci with missing data less than 2% and MAF greater than 5% were used (Supplementary Figure 3). About 26% of the markers had MAF between 5 and 10% whereas 31% had MAF between 11 and 20%. About 40% of the markers had MAF greater than 2% (Supplementary Figure 3 and Supplementary Table 4). The SNP markers had a moderately balanced distribution across the chromosomes, ranging from 252 SNPs (8.4%) on chromosome 8–371 SNPs (12.4%) on chromosome 3 (Supplementary Table 4).
The site frequency spectrum revealed high variation in the MAF distribution of the SNPs among the 24 sorghum landrace accessions (Figure 1). All individuals in 67% of the accessions (horizontally, the first 16 accessions in Figure 1) had major alleles in most of their SNP loci although to a different extent. Interestingly, one individual from each of the remaining eight accessions carried minor alleles across most of their SNP loci. The expected site-frequency spectrum determined using a coalescent approach matched the observed frequency distributions fairly well for the accessions, SB2, SB7, and SB19 while it was inversely related to the observed frequency distributionsin accession SB4 (Figure 1). The genome-wide diversity across sorghum landrace accessions was quantified using a sliding window approach (window size = 1 Mbp, step size = 200 kb) to explore the genomic signatures of diversity in sorghum. The analyses resulted in an overall average nucleotide diversity (π) and Tajima’s D of 1.2 and 1.4/Mb, respectively (Figure 2). It is clear from Figure 2 that SNPs representing the centromere regions of each chromosome almost do not exist among the 3,001 SNPs used in this study, and thus the diversity estimates were extremely low or zero. At chromosome level, the highest average nucleotide diversity and Tajima’s D were recorded in chromosome 9 (1.5 and 2.2/Mb) and the lowest in chromosome 7 (1.0/Mb) and chromosome 8 (0.7/Mb), respectively (Figure 2). Previously reported candidate loci for domestication are found at the genomic regions with notably low diversity on chromosomes 2, 4, and 7 (Figure 2). The effective number of alleles found across the 3,001 SNP markers ranged from 1.01 to 1.98 with a mean of 1.16. Observed heterozygosity (Ho) varied from 0.0 to 0.96 with a mean of 0.06 while the mean expected heterozygosity (He) was 0.10 with individual values per locus ranging from 0.01 to 0.49. Similarly, the gene diversity estimates per locus varied from 0.10 to 0.50 with a mean value of 0.29 (Figure 3 and Supplementary Table 4). The average PIC of the loci was 0.24 with individual values ranging from 0.09 to 0.37. In the case of fixation indices, the minimum, maximum, and mean values for FIS were –0.96, 1.00, and 0.45, for FIT they were –0.93, 1.00, and 0.79, and for FST, they were 0.01, 0.95, and 0.63, respectively. The estimates of gene flow (Nm) per locus showed wide variation, ranging from 0.01 to 20.23, with a mean of 0.20 (Figure 3 and Supplementary Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-g001.jpg)
Figure 1. The pattern of site frequency spectrum based on the proportion of the minor allele frequency (MAF) of single nucleotide polymorphism (SNP) in the 24 sorghum accessions.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-g002.jpg)
Figure 2. Genome-wide pattern of diversity in the 359 individual plants representing the 24 sorghum accessions. A sliding window approach (window size = 1 Mb, step size = 200 kb) was used to analyze nucleotide diversity and Tajima’s D. The black vertical lines on chromosomes 2, 4, and 7 show the positions of shrunken2 (sh2), amylose extender1 (ae1), and brittle2 (bt2), respectively, which were previously identified as domestication loci in maize and sorghum that are localized at regions of low diversity. The overall average nucleotide diversity (π) and Tajima’s D were 1.2/Mb and 1.4/Mb, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-g003.jpg)
Figure 3. The mean, minimum (Min), and maximum (Max) values for the number of allele (Na), number of effective allele (Ne), Shannon informative index (I), observed heterozygosity (Ho), unbiased expected heterozygosity (uHe), expected heterozygosity (He), fixation indices (Fis, Fit, Fst), polymorphic information content (PIC), and gene diversity (H) for the 3001 polymorphic SNP loci.
Based on the HWE test, 99.5% of the SNP markers showed significant deviation from HWE (Supplementary Table 4). Among the 3,001 SNP loci, 97.9% were heterozygote-deficient, whereas 1.6% (48 loci) had excess heterozygosity showing significant deviation from HWE (P < 0.05). The candidate genes containing SNP markers showing excess heterozygosity and their annotated functions were retrieved from SorGSD6 and further evaluated. Among the 48 SNP loci that showed excess heterozygosity, nine SNPs lacked one of the three possible genotypes expected in a bi-allelic polymorphic locus under the assumption of HWE. The change in amino acid sequences of the corresponding genes was obtained in all SNPs, except three SNPs (snp_sb001000020838, snp_sb042060612417, and snp_sb042061102446) (Supplementary Table 5).
Genetic Diversity Analysis
The 3,001 polymorphic SNP markers revealed a wide range of variation in the Ethiopian sorghum germplasm, as estimated using different population genetics parameters across the 24 accessions, and are summarized in Table 1. The effective number of alleles of the accessions varied from 1.01 to 1.46 with a mean of 1.21, whereas the mean Shannon’s Information index (I) was 0.25 with individual values ranging from 0.0 (SB14 and SB15) to 0.42 (SB21). The lowest and the highest Ho values varied from 0.01 (SB5, SB14, and SB15 and SB16) to 0.25 (SB21) with a mean of 0.07. Likewise, the He and unbiased expected heterozygosity (uHe) of the accessions ranged from 0.0 to 0.27 and 0.0 to 0.28, respectively with a mean of 0.15 (Table 1). The lowest values were recorded in accessions, SB14, SB15, and SB16, whereas SB21 recorded the highest values for these parameters. The percent polymorphic loci (PPL) of the accessions varied from 0.8 to 91.4% with a mean of 47.7%. The fixation index (F) showed wide variation with values ranging from –0.76 (SB14) to 0.84 (SB3). Overall, accession SB21 showed the highest value for Ne, Ho, He, uHe, I, and the number of locally common alleles (NLCA) and PPL while SB14, SB15, and SB16 showed the lowest values for all genetic diversity parameters analyzed (Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-t001.jpg)
Table 1. Summary of different genetic diversity estimates based on 3,001 SNP markers for each of the 24 sorghum accessions and for a group of accessions grouped according to different agro-ecological zones (cool/semiarid, cool/subhumid, and warm/semiarid), geographical regions (eastern, northern and western), and peduncle shape (bent and erect).
Among the agro-ecological zones, warm/semiarid zones showed the highest values for Ne, Ho, He, uHe, and I whereas cool/subhumid zones showed the lowest in the majority of the genetic diversity estimates. Among the groups of accessions in the three agro-ecological zones, the highest value of PPL, which is equal to 99% and the number of private allele per locus were recorded in warm/semiarid and cool/semiarid zones (Table 1 and Supplementary Table 6). In terms of geographic regions, accessions from the western geographic regions showed the highest values in most of the genetic diversity parameters analyzed (I, Ho, He, and uHe) (Figure 4 and Table 1). For example, the eastern, northern, and western accessions had uHe values of 0.24, 0.21, and 0.37, respectively. Accessions from the eastern geographic region showed the highest value PPL, which is equal to 99.7% and in the number of private alleles. Four private alleles were recorded for the eastern region with MAF ranging from 0.19 to 0.47 whereas two and one private alleles were detected in the accessions originated from the western and northern regions, respectively (Table 1 and Supplementary Table 6). With regard to peduncle shape, accessions with bent peduncles were more diverse than those with erect peduncles as shown by the values of I, Ho He, uHe, and PPL (Figure 4 and Table 1). One hundred thirty-nine alleles were specific to accessions with bent peduncles shape with MAF ranging from 0.09 to 0.32, whereas only one private allele with MAF of 0.14 was specific to accessions with erect peduncles (Supplementary Table 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-g004.jpg)
Figure 4. Graphs displaying mean values of different genetic diversity parameters estimated based on 3,001 SNP markers for a group of sorghum accessions grouped according to their (A) geographic regions, (B) agro-ecological zones, and (C) peduncle shape. Na = No. of different alleles; Ne = effective number of alleles; I = Shannon’s information index; Ho = observed heterozygosity; uHe = unbiased expected heterozygosity; He = expected heterozygosity; F = fixation index; NPA = number of private alleles.
Genetic Differentiation of Accessions and Hierarchical Groups
The results of the AMOVA without grouping the accessions showed that 64.5% of the total variation was observed among accessions and 35.5% within accessions (FST = 0.65; FIS = 0.47, P < 0.001) (Table 2). Additionally, hierarchical AMOVA was conducted by grouping the accessions according to their geographic regions, agro-ecological zones of their collection sites, and their peduncle shape. In this analysis, 19.5% of the total variation was observed among the geographical regions, which is a highly significant differentiation (FCT = 0.20, P < 0.001). Similarly, significant differentiation was found among peduncle shape groups with 4.3% of the total variation between them (FCT = 0.04, P < 0.05) (Table 2). However, only 1.83% of the total variation accounted for the variation among the agro-ecological zones, which is statistically insignificant (FCT = 0.02 and P = 0.17) (Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-t002.jpg)
Table 2. Analysis of molecular variance (AMOVA) for 24 accessions without grouping, and by grouping them based on their geographic regions, agro-ecological zones, and peduncle shapes.
Population Differentiation and Gene Flow
The pairwise population differentiation analysis revealed significant differentiation among all pairs of accessions with FST values ranging from 0.18 to 0.99 (Figure 5 and Supplementary Table 7) except in the case of SB16 vs. SB12, which was not significant (FST = 0.02, P > 0.05). The pairs of accessions with the highest FST value (0.99) were SB14 vs. SB15 and SB15 vs. SB16, corresponding to the lowest estimate of gene flow (Nm = 0; Supplementary Table 7). The mean FST values for the differentiation of each accession from all other accessions varied from 0.47 to 0.81. Accessions SB15, SB14, and SB5 were the most differentiated with FST values of 0.81, 0.80, and 0.78, respectively, whereas SB18 was the least differentiated accession (FST = 0.47) (Figure 5 and Supplementary Table 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-g005.jpg)
Figure 5. Graphical display of pairwise genetic differentiation (FST) among the 24 sorghum accessions. The differentiation between each pair was significant (P < 0.05) except in the case of SB12 vs. SB16.
The analyses of the average number of pairwise differences (πxy) and the pairwise net number of allele differences (Nei’s distance, d) between the accessions revealed a wide variation with πxy ranging from 1.2 (SB12 vs. SB16) to 1,197.2 (SB6 vs. SB15) and d ranging from 0.001 (SB12 v.s SB16) to 0.56 (SB21 vs. SB6 and SB1 vs. SB14) (Figure 6 and Supplementary Table 8). Accessions SB1, SB6, and SB21 also showed a higher pairwise net number of allele differences (Nei’s distance, d) with other accessions (Figure 6 and Supplementary Table 8). In line with the results of the pairwise FST analysis, the average number of pairwise differences and Nei’s distance were the lowest for SB16 vs. SB12 suggesting that these two accessions are genetically very similar. The average number of pairwise differences within accessions also showed wide variation with the values ranging from 10 (SB14) to 840 (SB21). This parameter was very low for SB15 and SB16, as with SB14 (Figure 6 and Supplementary Table 8).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-g006.jpg)
Figure 6. Average number of pairwise differences within and between sorghum accessions: average number of pairwise differences among the accessions (πxy) (above diagonal, in green), average number of pairwise differences within accessions (π) (diagonal, in orange), and pairwise net number of allele differences among the accessions (d) (below diagonal, in blue).
At the geographic region level, the pairwise FST values among each pair of the three groups were significant (P < 0.001). However, accessions from the western region were the most distinct with higher differentiation from those from the northern and eastern regions (FST = 0.40 and 0.35, respectively). Among the three pairs, accessions from the eastern vs. northern regions were the least differentiated (FST = 0.12) (Supplementary Table 7). Similar to that of geographic regions, the FST values among each pair of the accessions from the three agro-ecological zones were also significant (P < 0.001). The accessions belonging to the cool/subhumid group were the most differentiated having FST values of 0.13 and 0.10 against warm/semiarid and cool/semiarid groups, respectively. The warm/semiarid vs. cool/semiarid groups were the least differentiated (FST = 0.08) among the three pairs (Supplementary Table 7). The average number of pairwise differences and Nei’s distance among the geographic regions and agro-ecological zones had a similar pattern with that of pairwise FST-based differentiation, revealing that the western region was the most differentiated group. In the case of pairwise differences within regions, accessions from the western region had the highest variation whereas the lowest was recorded for the northern region. With regard to agro-ecological zones, warm/semiarid and cool/subhumid zones showed the highest and lowest variations, respectively (Supplementary Table 7).
The non-hierarchical finite island model-based analysis involving the examination of the joint distribution of FST and heterozygosity among accessions to detect loci under selection revealed 74 SNP loci that were highly significant (P < 0.01). Among them, 61 loci had low FST value (ranging from –0.02 to 0.35), and hence were considered as candidates for balanced selection. Whereas 13 loci (Table 3) had high FST values (ranging from 0.81 to 0.94), and hence considered as under directional selection. The MAF of these loci ranged from 0.05 to 0.34. The markers were distributed on chromosomes 1, 5, 6, 7, and 9 with over 50% of them located on chromosome 7 (Table 3). The candidate genes containing these SNP markers and their putative functions were identified through BLAST searching the sorghum v3.1 genome at Phytozome 12.1 (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-t003.jpg)
Table 3. The list of 13 SNP loci that were identified as loci under selection and their descriptions.
Cluster Analyses of Individual Genotypes and Accessions
The unweighted pair group method with arithmetic mean-based cluster analysis of the 359 individual genotypes generated a dendrogram of three major clusters, which were denoted by different line colors in Figure 7. The cluster analysis at the accession level resulted in the clustering of the 24 accessions into two groups. In the case of individual genotypes, Cluster I consisted of 316 individuals, whereas Cluster II comprised 43 individuals, respectively. The cluster analysis showed that at least the majority of individuals from the same accessions were clustered together (Figure 7). All individuals of an accession were clustered closely together in several cases. For example, all individuals from accessions, SB21 and SB4 were clustered in Cluster I and Cluster II, respectively. In other cases, a few individuals of an accession were placed under different clusters. For instance, two individuals from accession SB1, one individual from SB6, SB17, and SB22 were separated from the other members of their accession and grouped with other genotypes in different clusters. Except in a few cases, most accessions were clearly clustered based on their geographic regions (Figure 8A). On the other hand, the clustering pattern of the accessions according to their agro-ecological zones or administrative regions was less resolved, as accessions were mostly clustered irrespective of their groups (Figures 8B,C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-g007.jpg)
Figure 7. Unweighted pair group method with arithmetic mean (UPGMA) dendrogram of 359 individuals representing the 24 sorghum accessions generated based on Nei’s genetic distance (Nei and Takezaki, 1983). The individual samples were coded in a way that the first two letters (SB) with either two- or three-digit numbers represent their accessions and the last two-digit numbers represent the codes for the individual plant in that accession. Individuals denoted by the same color and shape belong to the same accession.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-g008.jpg)
Figure 8. The UPGMA dendrogram of the 24 accessions generated based on Nei’s genetic distance (Nei and Takezaki, 1983) calculated using the genotypic data of 3,001 SNP markers. Accessions denoted by the same color labels belong to the same (A) geographic region, (B) agro-ecological zone, and (C) administrative region.
Principal Coordinate Analysis
Principal coordinate analysis was performed to determine the relationship between the sorghum accessions and individuals within the accession, which grouped the accessions into three separate clusters (Figure 9A and Supplementary Figure 4). The first and second coordinates explained 29.5 and 12.4% of the total variation among the accessions, respectively. Similar to the cluster analysis, PCoA revealed that accessions SB1, SB4, SB6, and SB21 are the most differentiated groups being clearly separated from the other accessions along the first principal coordinate (Figure 9A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/799482/fpls-12-799482-HTML/image_m/fpls-12-799482-g009.jpg)
Figure 9. (A) Principal coordinates analysis (PCoA) showing the clustering pattern of the 24 Ethiopian sorghum accessions and accessions denoted by the same color labels and shapes belonging to the same geographic region and (B) a graphical display of the population genetic structure of the 24 sorghum accessions for K = 2. The two colors represent the two clusters (genetic populations) and each color of an accession represents the average proportion of the alleles that placed that accession under the corresponding clusters.
Population Structure
The admixture model-based population structure of the 359 individuals representing the 24 accessions was inferred using STRUCTURE software. The analysis of the STRUCTURE output using STRUCTURE HARVESTER program (Earl, 2012) that implemented ΔK method of Evanno et al. (2005) revealed that the optimal number of genetic clusters is two (Supplementary Figure 5). The results suggest that the 24 sorghum accessions originate from two genetic populations as graphically depicted in Figure 9B. In line with the results of the cluster analysis and PCoA, accessions SB1, SB4, SB6, and SB21 were significantly differentiated groups, as the majority of their alleles belong to a different genetic population (represented by orange in Figure 9B) as compared to the other accessions.
Discussion
SNP Markers and Their Use in Genetic Diversity Analysis of Sorghum Gene Pool
Genetic diversity analysis of crop species is an important step in detecting alleles that could be used for their improvement through breeding. The Ethiopian sorghum gene pool has been used as a novel source of biotic and abiotic stress tolerance, greatly contributing to the improvement of sorghum, globally (Adugna, 2014). The gene pool has been utilized in various studies that aimed at the identification of novel QTLs and genes governing complex traits (Cuevas and Prom, 2013, 2020; Cuevas et al., 2017; Menamo et al., 2021). Since polymorphism within genes or their close vicinity is expected to be the main basis of phenotypic variation, priority was given to SNPs located in genes in the SNP selection process in this study. Because of simplicity and abundance in plant genomes, bi-allelic SNPs are the most commonly used SNPs used in genetic analyses. In the present study, 86% of the genotyped bi-allelic SNP loci were polymorphic, which can be considered high. This is most likely because, the SNP selection was mainly made based on the SNPs recorded for the Ethiopian sorghum genotype, Cherekit at the SorGSD database. The vast majority of the SNP loci (94.1%) were homozygous across the 359 individual samples genotyped, which is not surprising as sorghum is a self-pollinating crop.
The variation in allele frequency distribution among accessions shown by the analysis of site frequency spectrum indicates a high level of genetic diversity in the Ethiopian sorghum. Accessions containing individual genotypes dominated by minor alleles across the loci require further investigations to reveal the phenotypic diversity of desirable traits. The overall nucleotide diversity (π) of 1.2 recorded in this study is in agreement with the result of a previous study on sorghum landraces (Mace et al., 2013). However, it is higher than the values reported in some other studies on sorghum (Morris et al., 2013; Yan et al., 2018). Similarly, the overall Tajima’s D value recorded in this study was 1.4, which is lower and higher than values reported in Morris et al. (2013) and (Mace et al., 2013), respectively. Among the seven starch-related genes, amylose extender1 (ae1), brittle2 (bt2), Opaque2 (O2), shrunken1 (sh1), shrunken2 (sh2), sugary1 (su1), and waxy1 (wx1), previously identified as candidates of domestication loci (Whitt et al., 2002; De Alencar Figueiredo et al., 2010; Morris et al., 2013), three of them (sh2, ae1, and bt2) are found at the genomic regions with notably low diversity on chromosomes 2, 4, and 7, respectively (Figure 2). The bt2 gene on chromosome 7 coding for a starch biosynthesis enzyme has been shown to be a likely domestication locus in sorghum and maize (Whitt et al., 2002; De Alencar Figueiredo et al., 2010; Morris et al., 2013). The low recombination rates in the pericentromeric region or the presence of other loci under selection in this region may be the reason for the low diversity in the present study and previous studies on sorghum (Morris et al., 2013).
The average Ho of 0.06 obtained in the present study was in line with the results of previous studies on sorghum employing SNP markers (Cuevas et al., 2017) and SSR markers (Ng’uni et al., 2011) (Ho = 0.04), (Ramu et al., 2013; Motlhaodi et al., 2014) (Ho = 0.09), (Motlhaodi et al., 2017) (Ho = 0.03). The Ho was expected, as sorghum is a predominantly a self-pollinating crop (Poehlman and Sleper, 1979). Gene diversity (H) and PIC are the most common measures of polymorphism of markers, which shed light on the evolutionary pressure on the alleles and the mutation rate at a locus over time (Shete et al., 2000; Wilkinson et al., 2012). The total genetic diversity in a population can be estimated through the analyses of a large number of informative markers across their genome (Melchiorre et al., 2013). The gene diversity of the SNP markers across all accessions in this study ranged from 0.1 to 0.50 with a mean of 0.29, which is high. Informative markers could be used for genotyping populations for genetic diversity studies, and the informativeness of the markers can be measured by their PIC value (Salem and Sallam, 2016). In the case of bi-allelic SNP markers, the maximum PIC value of 0.375 is attained when both alleles have a frequency of 0.5. In the present bi-allelic SNP-based study, the PIC values ranged from 0.09 to 0.375 with the overall average of 0.24, which is comparable with previous studies on sorghum using SNP markers (Afolayan et al., 2019; Silva et al., 2021; Wondimu et al., 2021). Forty-seven percent of these SNP loci have a PIC value of greater than 0.25 and hence they are highly informative and could be used for various applications including population genetic studies of sorghum.
Selections, both natural and artificial, as well as inbreeding, contribute to the deviation of populations from HWE. In this study, 98% of the loci showed heterozygote deficiency while 1.60% of the loci showed excess heterozygosity. Since sorghum is a predominately self-pollinating species, heterozygote deficiency at the vast majority of the loci can be attributed to inbreeding. However, the small proportion of loci showing excess heterozygosity suggests that they could be under selection or linked to loci under selection. Among the SNP loci that showed excess heterozygosity, nine loci lacked one of the two homozygous genotypes. The data suggest that one of the two alleles at each locus reduces the fitness of homozygous genotypes, or the locus is linked to another locus within its corresponding gene or the nearby gene that has a significant fitness value. Most of these SNP markers are within the coding region of genes. For instance, snp_sb001000687053, snp_sb001000723312, snp_sb042060543510, snp_sb042060515233, and snp_sb042060517985 are within the coding region of senescence-related gene 1, tetratricopeptide repeat (TPR)-like superfamily protein, cysteine proteinases superfamily protein, C-terminal domain phosphatase-like 4, and hydroxyproline-rich glycoprotein family protein, respectively. These genes had major roles in the growth, development, physiology, and biotic and abiotic stress tolerances in plants. For instance, hydroxyproline-rich glycoproteins (HRGPs) play a major role in the growth and development of plants (Showalter et al., 2016) while cysteine proteinases play an important physiological process ranging from seed germination (Becker et al., 1994) to senescence (Valpuesta et al., 1995). Therefore, further study that investigates the effect of these SNPs on the response of sorghum to abiotic and biotic stresses is of high significance.
Genetic Diversity Within Accessions
The average He (0.15), I (0.25), and PPL (47.7%) obtained in the present study suggest low genetic variation within the sorghum accessions. In general, the relatively low genetic variation within landrace accessions in the present and previous studies on sorghum (Ng’uni et al., 2011; Motlhaodi et al., 2017) is likely due to the combination of its inbreeding nature and due to the strict selection criteria of farmers. However, the variation within accessions varied widely. In this regard, accessions from the western region [Benishangul-Gumuz, Gambella, and Southern Nations, Nationalities and Peoples’ Region (SNNPR)] had higher variation than other accessions with SB21 being the most diverse accession followed by SB4. Accessions from this region (SB21, SB1, and SB4) are characterized by bent peduncle and light brown seeds with the exception of the red seed color of SB1.
Mengistu et al. (2020) also reported higher gene diversity and PIC for accessions from the Benshangul-Gumuz, Gambella, and SNNP regions as compared to the other regions in Ethiopia. The higher variation within accessions from these regions may suggest less human selection pressure on the landraces as compared to sorghum grown elsewhere in the country. Since the genetic diversity of populations implies their potential to adapt to environmental changes (Markert et al., 2010), sorghum landraces from this region may serve as a potential source of genes for biotic and abiotic stresses. Another interesting result of this study is a significantly higher Ho in two of the three accessions (SB4 and SB21) from the western regions as compared to all other accessions. Higher Ho suggests a higher outcrossing rate in these accessions, which might have allowed for gene flow through pollen and hence increased the variation within the accessions. The results suggest the western region as an important source of sorghum genotypes with desirable traits, such as tolerance to biotic and abiotic stresses. On the other hand, most accessions from Northern Ethiopia (Tigray and Amhara) had very low variation within accessions. Their average Ho was 0.03, indicating that the vast majority of the loci in the genotypes of these accessions were homozygous. In this group, accessions SB5, SB14, SB15, and SB16 can be regarded as pure lines, as individuals within each accession are almost identical across the whole loci. On the other hand, other accessions in this group (SB6, SB13, and SB22) are more diverse although their heterozygosity is still very low. Since the loss of heterozygosity increases the chance of deleterious recessive alleles being expressed in the progeny (Radosavljević et al., 2015), these accessions may be more susceptible to biotic or abiotic stresses unless they have been selected for tolerance against these stresses over time.
Genetic Differentiation of Accessions and Hierarchical Groups
In this study, most of the total variations (64.5%) were observed among the accessions than within the accessions (35.5%). The lower genetic variation within the accessions is expected in self-pollinating crops like sorghum (Hamrick, 1983). In addition, strict farmers’ selection for crop improvement might have contributed to the lower within-accession variation, which were clearly displayed in accessions, such as SB14 and SB15. Previous genetic diversity studies through SNP and SSR markers also showed a higher genetic variation among sorghum accessions than within the accessions. For instance, SNP-based genetic diversity study on sorghum accessions from Ethiopia showed that the variation among and within the accessions accounted for 59.6 and 40.4% of the total variation, respectively (Mengistu et al., 2020). Similarly, genetic diversity study through SSR markers on sorghum accessions from Zambia revealed 82 and 18% genetic variations among and within the accession variations, respectively (Ng’uni et al., 2011). Motlhaodi et al. (2017) reported a significant genetic variation among 22 accessions of sorghum, which accounted for 66.9% of the total variation while the within accession variation accounted for 23.6%. However, high genetic variation within sorghum accessions were reported on sorghum studied through SNP markers (Afolayan et al., 2019) and SSR markers (Manzelli et al., 2007; Adugna, 2014), suggesting that the accessions are not under selection processes.
Several studies have shown that the diversity of sorghum is associated with geography, agro-ecology, ethnicity, or botanical racial classifications (Barnaud et al., 2007; Ng’uni et al., 2011; Faye et al., 2019; Menamo et al., 2021). Significant genetic variations among the geographic regions and peduncle shape groups were observed in this study as shown by hierarchical AMOVA. Among the geographic regions, the western and eastern regions had higher genetic diversity than the northern region as shown by average He and the percentage of polymorphic loci, which were higher than the overall average (He = 0.24 and PIC = 89%). The western region accessions were the most distinct, with higher differentiation from those of the northern and eastern geographic regions. The major sorghum growing area (northern region) of the country had relatively low genetic variation probably due to intensive farmers’ selection of landraces to cope with the local environmental factors, such as the duration of the rainy season. The diversity of the crop has been reduced over time due to the recurrent drought in this major sorghum-growing region of the country. Overall, farmers in the drought-prone lowland areas tend to use early maturing and high yielding types and or shift their production systems to more vulnerable and low yielding early maturing crop species, such as tef (Eragrostis tef) (Adugna, 2014), which may provide genetic erosion of the sorghum landraces in these regions. High adoption of early maturing improved varieties in drought-prone areas in the northern region was also reported (Tesfaye et al., 2013).
Private alleles represent a unique genetic variability at certain loci of a particular population or hierarchically grouped populations. In this study, private alleles were not detected at the population level, but were recorded in all geographic regions. The Eastern region had a higher number of private alleles as compared to the western and northern regions, and hence it may serve as a rich source of desirable alleles for sorghum improvement. Private alleles generally support the potential to respond to a selection or have evolutionary significance (Petit et al., 1998). Information on private alleles is crucial for selecting highly diverse genotypes that can be used in breeding programs as a source of parental lines for crossbreeding that would eventually lead to new cultivars enriched with desirable alleles (Brondani et al., 2006; De Oliveira Borba et al., 2009; Salem and Sallam, 2016). The presence of more private alleles in the eastern region suggests the good in situ conservation status of sorghum in that location. Hence, further studies that explore the region for highly desirable traits need to be conducted, especially for use in sorghum-breeding programs.
Sorghum genotypes showed a significant genetic differentiation based on their peduncle shape, possibly because the shape of the peduncle influences the mating system, with the architecture of very bent peduncle obstructing pollination with outcrossed pollen. A more interesting finding was that accessions with bent peduncles exhibited higher genetic variation on average than those with erect peduncles. Unlike previous studies on Ethiopian sorghum (Menamo et al., 2021; Wondimu et al., 2021), sorghum accessions were not significantly differentiated according to agro-ecology in this study. However, a high significant genetic difference among the three pairs of agro-ecological zones was observed and the warm/semiarid zones showed the highest genetic diversity among the agro-ecological zones. Private alleles were detected from warm/semiarid and cool semiarid zones. Cool/subhumid zones, however, did not exhibit any private allele.
In the present study, 13 SNP loci were identified as loci under selection through the determination of the joint distribution of FST and heterozygosity. More than 50% of these SNP loci are located on chromosome 7 of the sorghum genome, suggesting that this chromosome carries many genes under natural selection or targeted by farmers directly or indirectly during and after domestication. These SNPs include those located in genes coding for zinc finger CCCH type family protein, DEAD-box ATP-dependent RNA helicase 42 and pentatricopeptide (PPR) repeat-containing proteins, which play a crucial role in plant responses to biotic and abiotic stresses (Peng et al., 2012; Xing et al., 2018; Nidumukkala et al., 2019). Hence, further study on these loci using individual genotypes that carry different alleles may shed more light on their significance in terms of desirable traits.
The Clustering Pattern and Population Structure of the Sorghum Accessions
Unweighted pair group method with arithmetic mean clustering based on Nei’s genetic distance (Nei and Takezaki, 1983) placed the individuals from the 24 accessions into three clusters. In line with the generally low genetic variation within accessions revealed through different analyses, there was a clear clustering pattern of individual genotypes according to their accessions. At the accession level, the cluster analysis generated three distinct clusters that matched the three clusters of the PCoA, which explained 42% of the total variation in its first two principal axes. The STRUCTURE analysis also generally agrees with the observed clustering pattern although it suggested two genetic populations (K = 2) as the best representation of the germplasm studied. Most of the alleles of the most distinct clusters in UPGMA and PCoA analyses (containing SB1, SB4, SB6, and SB21) originate from the first genetic cluster of STRUCTURE analysis (shown orange in Figure 9). Hence, it is interesting to crossbreed individual genotypes in these accessions with genotypes of genetically uniform accessions (e.g., SB14 and SB15), and evaluate the progeny generations for desirable traits.
In this study, the significant differentiation among geographic groups but not among agro-ecological groups revealed through AMOVA was also evident in the cluster analysis at the level of accessions. Based on redundancy analysis in their recent study on sorghum, Menamo et al. (2021) reported that agro-ecology is more important than the administrative region in defining the genetic variation in sorghum, which is not in agreement with the present study. The present study showed that the genetic diversity of Ethiopian sorghum landrace accessions was more structured along the geographical regions than along the administrative regions or agro-ecological zones. The lack of clear genetic differentiation of sorghum along the administrative regions, which was also previously reported (Ayana and Bekele, 2000; Desmae et al., 2016; Wondimu et al., 2021), could be explained by a high gene flow because of extensive exchange of seeds among farmers across adjacent regions where sorghum is a major crop.
Conclusion
In this study, SeqSNP method was used to genotype diverse sorghum accessions using a combination of previously developed and newly identified gene-based SNP markers. Despite the fact that they were gene-based, the SNP markers revealed a comparable genetic variation from the previous studies using SNP markers in sorghum. About half of the SNP markers can be regarded as highly informative and can be prioritized for future population genetics studies. A significant number of loci exhibited excess heterozygosity and/or were presumed to be under selection, some of which are located within genes playing crucial roles in plant responses to biotic and abiotic stresses. Further research on these loci using genotypes carrying different alleles may shed light on their significance in terms of desirable traits. The observed highly significant genetic differentiation among the sorghum accessions will be beneficial to the sorghum breeders in selecting desirable parents for crossbreeding. The sorghum accessions formed three distinct clusters, and it is therefore interesting to crossbreed genotypes from different clusters to evaluate their progeny for desirable traits. In this study, highly significant variations were observed among the geographic regions and peduncle-shaped groups. Compared to the western and northern regions, the eastern region had a higher number of private alleles, and hence it may serve as a rich source of desirable alleles for improving sorghum. Lastly, given that sorghum is generally regarded as a self-pollinating species, an exceptionally high heterozygosity observed in accessions, namely, SB4 and SB21 from the western geographic region, is an interesting result of this study, and should be further investigated.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
MG and ME designed the experiment and analyzed the data. ME conducted the experiment and wrote the draft manuscript. AC, CH, KT, MG, and TF reviewed the manuscript. All authors conceived the study and read and approved the submission of the manuscript for publication.
Funding
This research was financially supported by the Swedish International Development Cooperation Agency (Sida) and the Research and Training Grant awarded to the Addis Ababa University and the Swedish University of Agricultural Sciences (AAU-SLU Biotech; https://sida.aau.edu.et/index.php/biotechnology-phd-program/; accessed on September 25, 2021).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank the Swedish International Development Cooperation Agency (Sida) for financing this research. We would also like to thank the Institute of Biotechnology, Addis Ababa University and Department of Plant Breeding, Swedish University of Agricultural Sciences, for technical support during the course of the study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.799482/full#supplementary-material
Supplementary Figure 1 | Geographical maps of Ethiopia showing (A) the original sampling locations of the sorghum accessions with red, green, and blue colors to highlight the western, northern, and eastern geographic regions, respectively and (B) the agro-ecological zones of Ethiopia as per the Global 16 Class classification system by Amede et al. (2015).
Supplementary Figure 2 | The panicle diversity of sorghum landraces grown in the country.
Supplementary Figure 3 | Minor allelic frequency (MAF) range of 3,001 SNP markers used for genetic diversity and population structure analyses of 359 individual plants representing the 24 sorghum accessions.
Supplementary Figure 4 | Principal coordinates analysis (PCoA) showing the clustering pattern of the 359 individuals of sorghum landraces and individuals denoted by the same color labels and shapes belonging to the same geographic region.
Supplementary Figure 5 | Inferred population structure of 24 sorghum accessions at K = 2. ΔK plot showing its maximum value at K = 2 suggesting two as the optimal number of genetic populations.
Footnotes
- ^ http://genome.jgi-psf.org/Sorbi1/Sorbi1.info.html
- ^ https://www.ebi.gov.et/biodiversity/conservation/genetic-material-holdings/
- ^ https://biosearch-cdn.azureedge.net/assetsv6/sbeadex-plant-data-sheet.pdf
- ^ http://sorgsd.big.ac.cn
- ^ www.phytozome.net
- ^ https://ngdc.cncb.ac.cn/sorgsd/
References
Adugna, A. (2014). Analysis of in situ diversity and population structure in Ethiopian cultivated Sorghum bicolor (L.) landraces using phenotypic traits and SSR markers. SpringerPlus 3, 1–14. doi: 10.1186/2193-1801-3-212
Adugna, A., Snow, A. A., Sweeney, P. M., Bekele, E., and Mutegi, E. (2013). Population genetic structure of in situ wild Sorghum bicolor in its Ethiopian center of origin based on SSR markers. Genet. Resour. Crop Evol. 60, 1313–1328. doi: 10.1007/s10722-012-9921-8
Afolayan, G., Deshpande, S., Aladele, S., Kolawole, A., Angarawai, I., Nwosu, D., et al. (2019). Genetic diversity assessment of sorghum (Sorghum bicolor (L.) Moench) accessions using single nucleotide polymorphism markers. Plant Genet. Resour. 17, 412–420.
Ali, M., Rajewski, J., Baenziger, P., Gill, K., Eskridge, K., and Dweikat, I. (2008). Assessment of genetic diversity and relationship among a collection of US sweet sorghum germplasm by SSR markers. Mol. Breed. 21, 497–509. doi: 10.1007/s11032-007-9149-z
Amede, T., Auricht, C., Boffa, J.-M., Dixon, J. A., Mallawaarachchi, T., Rukuni, M., et al. (2015). The evolving farming and pastoral landscapes in ethiopia: a farming system framework for investment planning and priority setting. Canberra, ACT: ACIAR.
Ayana, A., and Bekele, E. (2000). Geographical patterns of morphological variation in sorghum (Sorghum bicolor (L.) Moench) germplasm from Ethiopia and Eritrea: quantitative characters. Euphytica 115, 91–104.
Ayana, A., Bryngelsson, T., and Bekele, E. (2000). Genetic variation of Ethiopian and Eritrean sorghum (Sorghum bicolor (L.) Moench) germplasm assessed by random amplified polymorphic DNA (RAPD). Genet. Resour. Crop Evol. 47, 471–482. doi: 10.1111/j.1601-5223.2000.t01-1-00249.x
Barnaud, A., Deu, M., Garine, E., Mckey, D., and Joly, H. I. (2007). Local genetic diversity of sorghum in a village in northern Cameroon: structure and dynamics of landraces. Theor. Appl. Genet. 114, 237–248. doi: 10.1007/s00122-006-0426-8
Becker, C., Fischer, J., and Münitz, K. (1994). PCR cloning and expression analysis of cDNAs encoding cysteine proteinases from germinating seeds of Vicia sativa L. Plant Mol. Biol. 26, 1207–1212. doi: 10.1007/BF00040701
Borrell, A. K., Hammer, G. L., and Douglas, A. C. (2000). Does maintaining green leaf area in sorghum improve yield under drought? I. Leaf growth and senescence. Crop Sci. 40, 1026–1037.
Brondani, C., Caldeira, K. D. S., Borba, T. C. O., Pn, R., De Morais, O. P., Castro, E. D. M., et al. (2006). Genetic variability analysis of elite upland rice genotypes with SSR markers. Embrapa Arroz Feijão Artigo periódico indexado 6, 9–17. doi: 10.12702/1984-7033.v06n01a02
Bucheyeki, T. L., Gwanama, C., Mgonja, M., Chisi, M., Folkertsma, R., and Mutegi, R. (2009). Genetic variability characterisation of Tanzania sorghum landraces based on simple sequence repeats (SSRs) molecular and morphological markers. Afr. Crop Sci. J. 17:54201.
Burow, G., Franks, C. D., Xin, Z., and Burke, J. J. (2012). Genetic diversity in a collection of Chinese sorghum landraces assessed by microsatellites. Am. J. Plant Sci. 3, 1722–1729. doi: 10.4236/ajps.2012.312210
Cortés, A. J., and Blair, M. W. (2018). Genotyping by sequencing and genome–environment associations in wild common bean predict widespread divergent adaptation to drought. Front. Plant Sci. 9:128. doi: 10.3389/fpls.2018.00128
Cuevas, H. E., and Prom, L. K. (2013). Assessment of molecular diversity and population structure of the Ethiopian sorghum [Sorghum bicolor (L.) Moench] germplasm collection maintained by the USDA–ARS National Plant Germplasm System using SSR markers. Genet. Resour. Crop Evol. 60, 1817–1830. doi: 10.1007/s10722-013-9956-5
Cuevas, H. E., and Prom, L. K. (2020). Evaluation of genetic diversity, agronomic traits, and anthracnose resistance in the NPGS Sudan Sorghum Core collection. BMC Genomics 21:88. doi: 10.1186/s12864-020-6489-0
Cuevas, H. E., Rosa-Valentin, G., Hayes, C. M., Rooney, W. L., and Hoffmann, L. (2017). Genomic characterization of a core set of the USDA-NPGS Ethiopian sorghum germplasm collection: implications for germplasm conservation, evaluation, and utilization in crop improvement. BMC Genomics 18, 1–17. doi: 10.1186/s12864-016-3475-7
De Alencar Figueiredo, L. F., Sine, B., Chantereau, J., Mestres, C., Fliedel, G., Rami, J.-F., et al. (2010). Variability of grain quality in sorghum: association with polymorphism in Sh2, Bt2, SssI, Ae1, Wx and O2. Theor. Appl. Genet. 121, 1171–1185. doi: 10.1007/s00122-010-1380-z
De Oliveira Borba, T. C., Brondani, R. P. V., Rangel, P. H. N., and Brondani, C. (2009). Microsatellite marker-mediated analysis of the EMBRAPA Rice Core Collection genetic diversity. Genetica 137, 293–304. doi: 10.1007/s10709-009-9380-0
De Wet, J., and Harlan, J. (1971). The origin and domestication ofSorghum bicolor. Econ. Bot. 25, 128–135.
Desmae, H., Jordan, D. R., and Godwin, I. D. (2016). Geographic patterns of phenotypic diversity in sorghum (Sorghum bicolor (L.) Moench) landraces from North Eastern Ethiopia. Afr. J. Agricult. Res. 11, 3111–3122.
Djè, Y., Heuertz, M., Lefebvre, C., and Vekemans, X. (2000). Assessment of genetic diversity within and among germplasm accessions in cultivated sorghum using microsatellite markers. Theor. Appl. Genet. 100, 918–925. doi: 10.1007/s001220051371
Earl, D. A. (2012). STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361. doi: 10.1007/s12686-011-9548-7
Ejeta, G. (2005). “Integrating biotechnology, breeding, and agronomy in the control of the parasitic weed Striga spp in sorghum,” in the wake of the double helix: from the green revolution to the gene revolution, Bologna Bologna, eds R. Tuberosa, R. L. Phillips, and M. Gale (Bologna: Avenue Media), 239–251.
Enyew, M., Feyissa, T., Geleta, M., Tesfaye, K., Hammenhag, C., and Carlsson, A. S. (2021). Genotype by environment interaction, correlation, AMMI, GGE biplot and cluster analysis for grain yield and other agronomic traits in sorghum (Sorghum bicolor L. Moench). PLoS One 16:e0258211. doi: 10.1371/journal.pone.0258211
Evanno, G., Regnaut, S., and Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x
Excoffier, L., and Lischer, H. E. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567. doi: 10.1111/j.1755-0998.2010.02847.x
Faye, J. M., Maina, F., Hu, Z., Fonceka, D., Cisse, N., and Morris, G. P. (2019). Genomic signatures of adaptation to Sahelian and Soudanian climates in sorghum landraces of Senegal. Ecol. Evolut. 9, 6038–6051. doi: 10.1002/ece3.5187
Garrison, E., and Marth, G. (2012). Haplotype-based variant detection from short-read sequencing. arXiv. [Preprint].
Geleta, N., and Labuschagne, M. (2005). Qualitative traits variation in sorghum (Sorghum bicolor (L.) Moench) germplasm from, eastern highlands of Ethiopia. Biodivers. Conserv. 14, 3055–3064. doi: 10.1007/s10531-004-0315-x
Ghebru, B., Schmidt, R., and Bennetzen, J. (2002). Genetic diversity of Eritrean sorghum landraces assessed with simple sequence repeat (SSR) markers. Theor. Appl. Genet. 105, 229–236. doi: 10.1007/s00122-002-0929-x
Girma, G., Nida, H., Seyoum, A., Mekonen, M., Nega, A., Lule, D., et al. (2019). A large-scale genome-wide association analyses of Ethiopian sorghum landrace collection reveal loci associated with important traits. Front. Plant Sci. 10:691. doi: 10.3389/fpls.2019.00691
Govindaraj, M., Vetriventhan, M., and Srinivasan, M. (2015). Importance of genetic diversity assessment in crop plants and its recent advances: an overview of its analytical perspectives. Genet. Res. Int. 2015:431487. doi: 10.1155/2015/431487
Hamrick, J. L. (1983). The distribution of genetic variation within and among natural plant populations. Genet. Conserv. 1983, 335–363.
Kopelman, N. M., Mayzel, J., Jakobsson, M., Rosenberg, N. A., and Mayrose, I. (2015). Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Resour. 15, 1179–1191. doi: 10.1111/1755-0998.12387
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Lai, J., Li, R., Xu, X., Jin, W., Xu, M., Zhao, H., et al. (2010). Genome-wide patterns of genetic variation among elite maize inbred lines. Nat. Genet. 42, 1027–1030. doi: 10.1038/ng.684
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Liu, K., and Muse, S. P. (2005). New Genetic Data Analysis S oftware. Massachusetts, MA: Whitehead Institute.
Luo, H., Zhao, W., Wang, Y., Xia, Y., Wu, X., Zhang, L., et al. (2016). SorGSD: a sorghum genome SNP database. Biotechnol. Biofuels 9, 1–9.
Rozas, J., Sánchez-Delbarrio, J.C., Messeguer, X., and Rozas, R. (2003). DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19, 2496–2497. doi: 10.1093/bioinformatics/btg359
Mace, E. S., Tai, S., Gilding, E. K., Li, Y., Prentis, P. J., Bian, L., et al. (2013). Whole-genome sequencing reveals untapped genetic potential in Africa’s indigenous cereal crop sorghum. Nat. Commun. 4, 1–9. doi: 10.1038/ncomms3320
Manzelli, M., Pileri, L., Lacerenza, N., Benedettelli, S., and Vecchio, V. (2007). Genetic diversity assessment in Somali sorghum (Sorghum bicolor (L.) Moench) accessions using microsatellite markers. Biodivers. Conserv. 16, 1715–1730. doi: 10.1007/s10531-006-9048-3
Markert, J. A., Champlin, D. M., Gutjahr-Gobell, R., Grear, J. S., Kuhn, A., Mcgreevy, T. J., et al. (2010). Population genetic diversity and fitness in multiple environments. BMC Evol. Biol. 10, 1–13. doi: 10.1186/1471-2148-10-205
McCormick, R. F., Truong, S. K., Sreedasyam, A., Jenkins, J., Shu, S., Sims, D., et al. (2018). The Sorghum bicolor reference genome: improved assembly, gene annotations, a transcriptome atlas, and signatures of genome organization. Plant J. 93, 338–354. doi: 10.1111/tpj.13781
Mehmood, S., Bashir, A., Ahmad, A., Akram, Z., Jabeen, N., and Gulfraz, M. (2008). Molecular characterization of regional Sorghum bicolor varieties from Pakistan. Pak. J. Bot 40, 2015–2021.
Melchiorre, M. G., Chiatti, C., Lamura, G., Torres-Gonzales, F., Stankunas, M., Lindert, J., et al. (2013). Social support, socio-economic status, health and abuse among older people in seven European countries. PLoS One 8:e54856. doi: 10.1371/journal.pone.0054856
Menamo, T., Kassahun, B., Borrell, A., Jordan, D., Tao, Y., Hunt, C., et al. (2021). Genetic diversity of Ethiopian sorghum reveals signatures of climatic adaptation. Theor. Appl. Genet. 134, 731–742. doi: 10.1007/s00122-020-03727-5
Mengistu, G., Shimelis, H., Laing, M., Lule, D., Assefa, E., and Mathew, I. (2020). Genetic diversity assessment of sorghum (Sorghum bicolor (L.) Moench) landraces using SNP markers. S. Afr. J. Plant Soil 37, 220–226. doi: 10.1080/02571862.2020.1736346
Mofokeng, A., Shimelis, H., Tongoona, P., and Laing, M. (2014). A genetic diversity analysis of South African sorghum genotypes using SSR markers. S. Afr. J. Plant Soil 31, 145–152. doi: 10.1080/02571862.2014.923051
Mohammadi, S. A., and Prasanna, B. (2003). Analysis of genetic diversity in crop plants—salient statistical tools and considerations. Crop Sci. 43, 1235–1248. doi: 10.1186/s12864-017-3922-0
Morris, G. P., Ramu, P., Deshpande, S. P., Hash, C. T., Shah, T., Upadhyaya, H. D., et al. (2013). Population genomic and genome-wide association studies of agroclimatic traits in sorghum. Proc. Natl. Acad. Sci. 110, 453–458. doi: 10.1073/pnas.1215985110
Motlhaodi, T., Geleta, M., Bryngelsson, T., Fatih, M., Chite, S., and Ortiz, R. (2014). Genetic diversity in’ex-situ’conserved sorghum accessions of Botswana as estimated by microsatellite markers. Austral. J. Crop Sci. 8, 35–43.
Motlhaodi, T., Geleta, M., Chite, S., Fatih, M., Ortiz, R., and Bryngelsson, T. (2017). Genetic diversity in sorghum [Sorghum bicolor (L.) Moench] germplasm from Southern Africa as revealed by microsatellite markers and agro-morphological traits. Genet. Resour. Crop Evol. 64, 599–610. doi: 10.1007/s10722-016-0388-x
Nei, M., and Takezaki, N. (1983). Estimation of genetic distances and phylogenetic trees from DNA analysis. Proc. 5th World Cong. Genet. Appl. Livstock. Prod. 21, 405–412.
Ng’uni, D., Geleta, M., and Bryngelsson, T. (2011). Genetic diversity in sorghum (Sorghum bicolor (L.) Moench) accessions of Zambia as revealed by simple sequence repeats (SSR). Hereditas 148, 52–62. doi: 10.1111/j.1601-5223.2011.02208.x
Ng’uni, D., Geleta, M., Hofvander, P., Fatih, M., and Bryngelsson, T. (2012). Comparative genetic diversity and nutritional quality variation among some important Southern African sorghum accessions [’Sorghum bicolor’(L.) Moench]. Austral. J. Crop Sci. 6, 56–64.
Nidumukkala, S., Tayi, L., Chittela, R. K., Vudem, D. R., and Khareedu, V. R. (2019). DEAD box helicases as promising molecular tools for engineering abiotic stress tolerance in plants. Crit. Rev. Biotechnol. 39, 395–407. doi: 10.1080/07388551.2019.1566204
Paterson, A. H., Bowers, J. E., Bruggmann, R., Dubchak, I., Grimwood, J., Gundlach, H., et al. (2009). The Sorghum bicolor genome and the diversification of grasses. Nature 457, 551–556. doi: 10.1038/nature07723
Peakall, R., and Smouse, P. (2012). GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research–an update. Bioinformatics 28:2537e2539. doi: 10.1093/bioinformatics/bts460
Peng, X., Zhao, Y., Cao, J., Zhang, W., Jiang, H., Li, X., et al. (2012). CCCH-type zinc finger family in maize: genome-wide identification, classification and expression profiling under abscisic acid and drought treatments. PLoS One 7:e40120. doi: 10.1371/journal.pone.0040120
Petit, R. J., El Mousadik, A., and Pons, O. (1998). Identifying populations for conservation on the basis of genetic markers. Conserv. Biol. 12, 844–855.
Pfeifer, B., Wittelsbürger, U., Ramos-Onsins, S. E., and Lercher, M. J. (2014). PopGenome: an efficient Swiss army knife for population genomic analyses in R. Mol. Biol. Evol. 31, 1929–1936. doi: 10.1093/molbev/msu136
Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959. doi: 10.1093/genetics/155.2.945
Radosavljević, I., Satovic, Z., and Liber, Z. (2015). Causes and consequences of contrasting genetic structure in sympatrically growing and closely related species. AoB Plants 7:lv106. doi: 10.1093/aobpla/plv106
Ramu, P., Billot, C., Rami, J.-F., Senthilvel, S., Upadhyaya, H., Reddy, L. A., et al. (2013). Assessment of genetic diversity in the sorghum reference set using EST-SSR markers. Theor. Appl. Genet. 126, 2051–2064. doi: 10.1007/s00122-013-2117-6
Rao, S. A., Rao, K. P., Mengesha, M., and Reddy, V. G. (1996). Morphological diversity in sorghum germplasm from India. Genet. Resour. Crop Evol. 43, 559–567. doi: 10.1007/bf00138832
Rao, V. R., and Hodgkin, T. (2002). Genetic diversity and conservation and utilization of plant genetic resources. Plant Cell Tiss. Org. Cult. 68, 1–19.
Ruiz-Chután, J. A., Salava, J., Janovská, D., Žiarovská, J., Kalousová, M., and Fernández, E. (2019). Assessment of genetic diversity in Sorghum bicolor using RAPD markers. Genetika 51, 789–803.
Salem, K. F., and Sallam, A. (2016). Analysis of population structure and genetic diversity of Egyptian and exotic rice (Oryza sativa L.) genotypes. C. R. Biol. 339, 1–9. doi: 10.1016/j.crvi.2015.11.003
Shete, S., Tiwari, H., and Elston, R. C. (2000). On estimating the heterozygosity and polymorphism information content value. Theor. Popul. Biol. 57, 265–271. doi: 10.1006/tpbi.2000.1452
Showalter, A. M., Keppler, B. D., Liu, X., Lichtenberg, J., and Welch, L. R. (2016). Bioinformatic identification and analysis of hydroxyproline-rich glycoproteins in Populus trichocarpa. BMC Plant Biol. 16, 1–34. doi: 10.1186/s12870-016-0912-3
Silva, K. J. D., Pastina, M. M., Guimarães, C. T., Magalhães, J. V., Pimentel, L. D., Schaffert, R. E., et al. (2021). Genetic diversity and heterotic grouping of sorghum lines using SNP markers. Sci. Agricola 78:39.
Singh, R., and Axtell, J. D. (1973). High Lysine Mutant Gene (hl that Improves Protein Quality and Biological Value of Grain Sorghum 1. Crop Sci. 13, 535–539.
Statista (2020). Sorghum production worldwide in 2019/2020, by leading country (in 1,000 metric tons). Hamburg: Statista.
Tajima, F. (1989). Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123, 585–595.
Tesfaye, E., Sime, M., Tirfesa, A., and Bekele, A. (2013). Socio-economic assessment of moisture stress sorghum growing areas of Miesso and Kobo districts. Addis Ababa: Ethiopian Institute of Agricultural Research.
Too, E. J., Onkware, A. O., Were, B. A. I, Gudu, S., Carlsson, A., and Geleta, M. (2018). Molecular markers associated with aluminium tolerance in Sorghum bicolor. Hereditas 155, 1–13. doi: 10.1186/s41065-018-0059-3
Tsehay, S., Ortiz, R., Johansson, E., Bekele, E., Tesfaye, K., Hammenhag, C., et al. (2020). New transcriptome-based SNP markers for Noug (Guizotia abyssinica) and their conversion to KASP markers for population genetics analyses. Genes 11:1373. doi: 10.3390/genes11111373
Valpuesta, V., Lange, N. E., Guerrero, C., and Reid, M. S. (1995). Up-regulation of a cysteine protease accompanies the ethylene-insensitive senescence of daylily (Hemerocallis) flowers. Plant Mol. Biol. 28, 575–582. doi: 10.1007/BF00020403
Wang, M. L., Zhu, C., Barkley, N. A., Chen, Z., Erpelding, J. E., Murray, S. C., et al. (2009). Genetic diversity and population structure analysis of accessions in the US historic sweet sorghum collection. Theor. Appl. Genet. 120, 13–23. doi: 10.1007/s00122-009-1155-6
Whitt, S. R., Wilson, L. M., Tenaillon, M. I., Gaut, B. S., and Buckler, E. S. (2002). Genetic diversity and selection in the maize starch pathway. Proc. Natl. Acad. Sci. 99, 12959–12962. doi: 10.1073/pnas.202476999
Wilkinson, P. A., Winfield, M. O., Barker, G. L., Allen, A. M., Burridge, A., Coghill, J. A., et al. (2012). CerealsDB 2.0: an integrated resource for plant breeders and scientists. BMC Bioinformatics 13, 1–6. doi: 10.1186/1471-2105-13-219
Wondimu, Z., Dong, H., Paterson, A. H., Worku, W., and Bantte, K. (2021). Genetic diversity, population structure and selection signature in Ethiopian Sorghum (Sorghum bicolor L.[Moench]) germplasm. bioRxiv. [Preprint].
Wu, Y., Huang, Y., Tauer, C., and Porter, D. R. (2006). Genetic diversity of sorghum accessions resistant to greenbugs as assessed with AFLP markers. Genome 49, 143–149. doi: 10.1139/g05-095
Xing, H., Fu, X., Yang, C., Tang, X., Guo, L., Li, C., et al. (2018). Genome-wide investigation of pentatricopeptide repeat gene family in poplar and their expression analysis in response to biotic and abiotic stresses. Sci. Rep. 8, 1–9. doi: 10.1038/s41598-018-21269-1
Keywords: agro-ecological zone, genetic differentiation, geographical region, population structure, sorghum [Sorghum bicolor (L.) Moench]
Citation: Enyew M, Feyissa T, Carlsson AS, Tesfaye K, Hammenhag C and Geleta M (2022) Genetic Diversity and Population Structure of Sorghum [Sorghum Bicolor (L.) Moench] Accessions as Revealed by Single Nucleotide Polymorphism Markers. Front. Plant Sci. 12:799482. doi: 10.3389/fpls.2021.799482
Received: 21 October 2021; Accepted: 03 December 2021;
Published: 05 January 2022.
Edited by:
Andrés J. Cortés, Colombian Corporation for Agricultural Research (AGROSAVIA), ColombiaReviewed by:
Zhenbin Hu, Saint Louis University, United StatesReza Darvishzadeh, Urmia University, Iran
Nemera Geleta Shargie, Agricultural Research Council of South Africa (ARC-SA), South Africa
Copyright © 2022 Enyew, Feyissa, Carlsson, Tesfaye, Hammenhag and Geleta. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Muluken Enyew, muluken.birara@aau.edu.et; mulukenbi@gmail.com; muluken.birara.enyew@slu.se