Skip to main content

Multielement analysis coupled with chemometrics modelling for geographical origin classification of teff [Eragrostis tef (Zuccagni) Trotter] grains from Amhara Region, Ethiopia


Teff [Eragrostis tef (Zuccagni) Trotter] is an indigenous crop in Ethiopia, and Amhara region is the predominant teff producing region in the country. This study was aimed at developing an analytical methodology useful to determine the geographical origin of teff produced in the Amhara Region, based on multielement analysis combined with multivariate statistical techniques. For this, a total of 72 teff grain samples were collected from three zones (West Gojjam, East Gojjam, and Awi) and analysed for K, Na, Mg, Ca, Mn, Cu, Fe, Co, Ni, Zn, Cr, and Cd contents using inductively coupled plasma-optical emission spectroscopy (ICP-OES). The digestion and ICP-OES analysis method were accurate, with percentage recovery ranging 85.5 to 109% across the different metals analysed. Principal component analysis (PCA) and linear discriminant analysis (LDA) were applied to discriminate samples based on their production regions. Magnesium, Ca, Fe, Mn, and Zn were the most discriminating elements among the samples. The LDA model provided 96% correct classification of samples into production regions and varietal types, with an average prediction ability of 92%. Hence, the multielement analysis combined with statistical modeling can be used in the authentication of the geographical origin and varietal type of teff from Amhara region.

Peer Review reports


Teff [Eragrostis tef (Zuccagni) Trotter] is a native cereal grain to Africa and is used as a staple crop in Ethiopia, Eritrea, Djibouti, Sudan's south-eastern region, and northern Kenya [1]. Teff is the most intensively produced cereal grain for human consumption in Ethiopia, followed by maize [2]. In comparison to other cereals like wheat, rice, oat, and barley, teff has a higher nutritional and fibre contents [3, 4]. Teff is a gluten-free cereal grain, and hence it is an ideal substitute for wheat and other cereals in food products [5, 6].

Teff is a warm seasonal grain that can be grown in wide climate conditions, such as in temperatures ranging from 10 to 27 °C and altitudes up to 2800 m above sea level [3, 7, 8]. In Ethiopia, teff covers about 2.8 million hectares, accounting for 25–30% of the total cereal acreage in the country [9, 10].

Amhara region is one of the major teff producer regions in the country. Within the region, East Gojjam, West Gojjam, and Awi zones are known for cultivating a significant amount of teff, and contributing the major share of the region’s teff production quantity. Although several varieties of teff are registered and available in the country, only three varieties are cultivated widely in the three zones. These are white teff, which includes Kuncho (DZ-CR-387) and Tseday (DZ-CR-37) varieties, red teff, locally called “Qey teff”, and mixed or ivory-colored teff, locally called “Sergegna teff”. The white Kuncho variety is mostly farmed in East and West Gojjam zones, while the white Tsedey, Qey, and Sergegna teff varieties are widely grown in Awi zone. The Kuncho and Tseday teff types are white in color and indistinguishable by visual observation.

In Ethiopia, teff grains are marketed by the names of their cultivation zones and districts. Consequently, the price depends on the cultivation region and varietal type [11, 12]. For example, white teff from Bichena district of East Gojjam zone is the most appreciated teff by consumers and is more expensive than others. However, inferior teff types are mislabelled and sold with the name of Bichena teff. Thus, a simple method that distinguishes the cultivation region of teff is in high demand to protect against fraudulent activities and mislabeling of inferior teffs in the market. It has to be noted that mislabelling would negatively affect both the consumer, farmer and legitimate merchant. In this respect, it is important to establish methods that facilitate recognition of the geographical origin of teff to protect both consumers and producers.

Multi-element analysis has been effectively used to discriminate the geographical origin of various food products and provide information about the cultivation environment [13, 14]. The discrimination of the geographical origin of various foodstuffs has been achieved using a combination of multivariate data analysis tools and the diverse elemental contents of food products. For instance, cereal grains such as rice, wheat, soybean, millet, tea, and coffee were classified to their respective region of origin based on their elemental composition together with some chemometrics models [15,16,17,18,19,20].

Thus, this study aimed to determine the levels of metals in teff grain from 21 sub- districts of three zones (East Gojjam, West Gojjam, and Awi zones) of Amhara region of Ethiopia by using inductively coupled plasma optical emission spectrometry (ICP-OES). The study also attempted to determine the variabilities of metal contents across the zones and evaluate the adequacy of the muti-element data in predicting the terroir of teff grains.

Chemicals, reagents and materials

Chemicals and reagents

The reagents HClO4 (70%) and HNO3 (69–72%) (Blulux,India) were used for the digestion of teff grains. CPI International, ICP quality control standards of multielements (USA) of K, Ca, Na, Fe, Mg, Mn, Zn, Ni, Cu, Cr, Co, and Cd were used in the study. The standard working solutions of the metals were prepared freshly from the stock standard solution.


Kjeldahl system (Velp Scientifica, Italy) was used for the digestion of teff samples. Inductively coupled plasma optical emission spectrometry (Perkin Elmer, Optima 8000) was used for the determination of multi-elements.

Sample collection and pre-treatment

During the 2020 harvest season, seventy two teff samples were directly collected from farmers in the three administrative zones (East Gojjam, West Gojjam, and Awi zones) of Amhara region, Ethiopia. Based on the information obtained from Amhara Region Bureau of Agriculture, these three zones are the country's most significant teff producing zones, which contribute voluminous teff grains to the market. From each zone, major teff producing districts and sub-districts were selected for sampling (Fig. 1). From East Gojjam zone, three districts (Aneded, Shebelberenta, and Enemay) and nine sub-districts (Jama Didik, Gudalema, Adisge Yegewera, Yeidwuha Town, Weregona Akababiw, Geda Iyesus, Weyira Gurazam, Hulet Amiba Dibisa, and Enekorna Adis Amba) were selected. From West Gojjam zone, two districts (Goinj kolela and Adet/Yilmanadensa) and seven sub-districts (Zanat, Gonj, Kore Tenkere, Menta Deber, Kilelet, Adet Zuria, and Mosbo) were selected. Where, five sub-districts (Gisayita, Gagasta, Ahiti, Zigem Town, and Guder Jawi) were chosen from Zigem district of Awi zone.

Fig. 1
figure 1

Source Reta et al. [21], published in Plos One, where the journal applies the Creative Commons Attribution (CC BY) license to make articles legally available for reuse, without permission or fees, for virtually any purpose

Map of Amhara Region showing teff sampling zones and sub-districts.

From each sampling site, 500 g of teff grains were collected.The samples were washed with tap water,rinsed twice with distilled water, and further rinsed with deionized water to eliminate adsorbed dust and particulate matter. The samples were oven dried for 12 h to remove moisture and obtain constant mass. The dried samples were ground using an electric grinder, sieved with a 0.5 mm pore-sized sieve and stored in bottles under airtight conditions until digestion and analysis.

Optimization of digestion procedure

Wet acid digestion is one of the methods involved in getting free metal ions in dissolved form from complex organic matrix. The digestion procedure was optimized by changing different parameters like volume ratio of reagents, temperature, and time (Table 1). The acid digestion was carried out by using Kjeldahl digestion system. It was assumed that digestion is complete when the solution becomes clear and colourless.

Table 1 Reagent volume, temperature and time attempted during optimization of digestion of 0.5 g of teff sample

Instrument calibration and sample analysis

Ten working standard solutions that contain the elements (Na, K, Mg, Ca, Cr, Mn, Fe, Co, Ni, Cu, Zn, and Cd) were prepared and used to calibrate the ICP-OES before the analysis of samples.The concentrations of working standard solutions were 0.05, 1, 2, 4, 6, 8, 10, 12, 14, and 16 mg/L for Cr, Mn, Co, Ni, Cu, Zn, and Cd while for Na, K, Mg, Ca, and Fe were 0.05, 1, 10, 20, 30, 40, 50, 60, 70, and 80 mg/L. A single measurement was performed for the calibration standards and the recovery experiment, while triple measurements were carried out on the digested teff samples. The accuracy of the method was assessed through recovery experiments, where samples were spiked with a known concentration of the standard solution of each metal and digested and analysed following the optimized procedure.

Method detection and quantification limits

For the determination of the limit of detection (LOD) of the analytical method, the standard deviation (SD) of the blank signals were multiplied by three and devided by the slop of the calibration equations. The limit of quantification (LOQ) was calculated as ten times the standard deviation of the blank devided by the slop of the calibration equations [22]. The wave length, method detection, quantification limit, correlation coefficient, and calibration curve equations are given in Table 2. From the correlation coefficients, it is possible to say that the change in absorbance with the concentration of each metal was in a good positive correlation and linearly fit.

Table 2 The wavelength used with ICP-OES, instrumental detection limit (IDL), limit of detection (LOD), limit of quantification (LOQ), correlation coefficient, and calibration equations

Statistical data analysis

Multivariate analysis methods used in this study included one-way analysis of variance (ANOVA), principal component analysis (PCA), linear discriminant analysis (LDA), which were performed with software SPSS 25.0 (SPSS Inc., Chicago, USA) and SIMCA 13 (Umetrics, Sweden). One-way ANOVA was used to test the presence of significant differences in the mean concentration of metals among samples. Duncan’s multiple comparisons was carried out to compare mean values of elements among samples. Differences were considered as statistically significant when p < 0.05. Principal component analysis (PCA) was used to reduce data dimensionality, to visualize sample trends and identify the most discriminating elements among samples [23, 24]. Linear discriminant analysis (LDA) was used to construct classification models useful to predict the geographical origin of teff samples.

Results and discussion

Analytical characteristics of the method

The efficiency of the analytical method was verified by testing its accuracy, sensitivity, repeatability, and linearity. The accuracy of the optimized digestion procedure was evaluated by analyzing the digests of spiked samples with standard solutions of each metal. The percentage recoveries varied in the range 85.5–109%, which signifies that the method possesses good accuracy.

Concentration of metals in teff samples

The teff samples collected from 21 sub-districts in three administrative zones had variable concentrations of metals, as indicated in Table 3. The concentration of K was the highest followed by Mg and Ca,while Cd, Co, and Cr were the lowest in concentration. The concentration (mg/kg) of metals in the teff samples ranged 1746.9–4850.2 K, 159.8–291.6 Na, 1391.7–3403.1 Mg, 951.6–2985.5 Ca, 0.6–17.7 Cr, 25.9–337.4 Mn, 94.6–900 Fe, 0.1–9.9 Co, 1.2—15.6 Ni, 4.7–16.7 Cu, 25.5–187.3 Zn, and 0.1–4.1 Cd.

Table 3 The concentration of metals (mean ± standard deviation) in mg/kg dry weight in teff samples from 21 sub-districts of three zones

Multi-elements variation among different regions

In general, the pattern of concentration of metals in teff samples from the three administrative zones was in decreasing order of K > Mg > Ca > Fe > Na > Mn > Zn > Cu > Cr > Ni > Co > Cd. The concentrations of Fe, Mn, Zn, Cu, Ca, Co, Mg, Na, and K in this study were in agreement with the previous studies [25,26,27]. The concentration of K was the most abundant element, followed by Mg and Ca in teff samples corroborated with the previous studies [26, 28,29,30]. A similar trend, i.e., a higher concentration of K, has also been reported in teff grains from other areas of Ethiopia [31]. The high concentrations of K, Ca, and Mg in teff grains could be due to the fact that these elements are readily available from the soil to the plant to translocate into the different parts of the plant [31, 32].

Regarding teff grains from the three zones, the mean concentrations of K and Na were highest in samples from East Gojjam, followed by Awi zone. Samples from the Awi zone were rich in the mean concentration of Mg, Ca, Fe, and Zn followed by samples from the West Gojjam zone (Table 4). The mean concentrations of the major elements in the teff samples were in the order of K > Mg > Ca > Na (Table 4). Similarly, the mean concentrations of trace elements in teff grains from Awi and West Gojjam zones were ordered as Cu > Ni > Cr > Co > Cd. Whereas the comparative order of these metals in samples from the East Gojjam zone was recorded as Cu > Cr > Ni > Co > Cd. This indicates the uptake of metals by plants takes place through different and complex biochemical processes. The uptake processes vary based on the nature of the metal, ability of the plants to absorb metals from the soil, the availability of the minerals in soluble and usable forms, the abundance of particular minerals in the particular areas, and the degree of contamination of the soil with heavy metals.

Table 4 Mean, maximum (Max), minimum (Min) concentration (mg/kg) and standard error (SE) the elements determined in teff grains from three zones of Amhara Region, Ethiopia

Statistical analysis with one-way ANOVA (α = 0.05) was used to test the presence of significant differences in the mean elemental contents of teff from the three administrative zones (Table 4). Samples from Awi zone contained significantly higher mean concentrations of Mg, Ca, Fe, and Zn than teff samples from East Gojjam and West Gojjam zones. On the other hand, no statistically significant difference was observed between the samples from East Gojjam and West Gojjam zones in their mean concentrations of Mg, Ca, Fe, and Zn. Concerning the other elements, except Mn and Cr, no significant difference was observed in the levels of Na, K, Ni, Cu, Cd, and Co among samples from the three zones. Significantly higher concentration of Mn was found in teff grains from Awi zone, while higher Cr in teff grains from East Gojjam. The observed significant variation among production zones might be due to the differences in the growing environmental conditions and crop inherent potential to take nutrients from the soil, resulting in differences in their metabolic activities like accumulations of different metals in their body [9, 33, 34].

The variations in the metal contents of teff grains from the different sub-districts was assessed using the standard deviation (SD) values (Table 3). Most of the sub-districts have been found to show less variation in their metal contents, which indicated that the sampling areas within a given sub-districts have similar agroecological conditions. Whereas, a higher variation was noted for some of the sub-districts, which indiced the existence of different agroecological conditions within the sub-districts. Some reports has indicated that variations in mineral contents in plants grown in different areas occurs when factors such as agricultural practices, soil chemistry, and plant genotype are varied [9, 33, 34].

Multi-elements variation among different varieties

Comparing the four varieties of teff collected from the three zones,Tseday teff variety was found to be rich in Mg (3403 mg/kg), Ca (2986 mg/kg), Fe (900 mg/kg), and Zn (187 mg/kg). The concentration of K (4850 mg/ kg) was highest in Kuncho teff variety when compared with the others (Table 3). The differences in the concentration levels of metals in varieties of teff arise mainly due to differences genetic properties, soil type, agricultural practices, and other environmental factors [31]. In the Amhara region, some farmers in some district use crop rotation systems to improve soil fertility, while other districts constantly cultivate the same crop type for many years at the same place [35]. Secondly, the soil types of the three sampling zones are different. Most East Gojjam teff samples are grown in black soil, while others are produced in grey and red soils [36, 37].

The concentrations of the trace elements (Cu, Zn, Co, Ni, Mn, and Cr) are in good agreement with the results reported by Habte et al. [26] and Baye et al. [38], as indicated in Table 5. Similarly, the determined concentrations of Mn and Cu lies within the range of values reported by Gebregewergis [39] and Habte et al. [26]. On the other hand, the determined concentrations of Ca and Zn are higher than that reported by Baye et al.[38], and the amounts of Co, Cu, and Cr were higher than that reported by Dame [28]. The literature generally indicates that the amount of metals in teff grains varies considerably with geographical origin, which can be attributed to variations in growing soil types, environmental conditions, and other factors, as it is also supported by previous findings of other researchers [31, 33, 34].

Table 5 Elemental composition (mg /kg dry weight) of teff samples reported in the literature

Principal component analysis

A preliminary investigation based on principal component analysis (PCA) was applied to visualize sample trends and evaluate the discriminatory characteristics of the determined elements. The data set was subjected to different scaling tools (Pareto, SNV, UV, freeze, and others) using the statistical software tool SIMCA 13 (Umetrics, Sweden). The Pareto scaling was found to provide excellent model parameters (Q2 (cum) = 0.568 and R2x (cum) = 0.922 with two principal components). Together, the first component (PC1) and second component (PC2) explained 82.8% of the variation in the data. The score plots and loading plots are shown in Figs. 2 and 3, respectively. The separation of the teff samples by the first two components is illustrated by the scores plot in Fig. 2, which indicates the presence of groups and some other important patterns in the data. The teff samples from Awi zone are clustered from the other zone teff samples and lie on the positive side of PC2, where two of the samples were found to be outliers. Teff samples from East Gojjam are presented on the negative side of PC1 and PC2, where two of them are scattered on the positive side of PC1 and three samples were found as outliers. At the same time, teff samples from West Gojjam zone lie on the negative side of PC1 and PC2. The elements responsible for separation in the first component were Cu, Mg,and Na, whereas separation along principal component 2 was influenced primarily by Ni, Ca, and Fe (Fig. 3).

Fig. 2
figure 2

Score plot of the first two principal components of the principal component analysis model, indicating the distribution of teff samples from three administrative zones of Amhara Region, Ethiopia

Fig. 3
figure 3

Loadings plot for the first two principal components obtained from principal component analysis metals in teff grains from three administrative zones of Amhara Region, Ethiopia

Teff grains from Awi zone were distinguished by their higher concentrations of Ca, Mg, and Fe. These observations are in agreement with the results obtained with one-way ANOVA. Similarly, teff grains from East Gojjam and West Gojjam were separated mainly due to their higher concentrations of Na and K. The loadings plot (Fig. 3) displays how the individual elements correlate with each other and contribute to the PCA model. Elements far away from the origin have a strong impact on the model, whereas those closer to the center have a weaker influence. From the plot, it is evident that K, Mg, Ca, and Fe play the largest role in discriminating the teff samples among the administrative zones. Hence,these elements are potentially good chemical descriptors for the construction of classification models for teff grains from the three zones.

Linear discriminant analysis at zones and districts

Evaluating the contributions of different factors towards variations in elemental contents of teff grains from different zones and varieties was not within the scope of this study. Instead, the study assumed that the different teff growing zones may have differences in one or more of different factors like local climate (rainfall, temperature, humidity, and solar radiation), soil type and agronomic practices that can result in significant variations in the elemental contents of teff grains cultivated in different areas. Hence, attempt was made to evaluate the adequacy of the overall variations in elemental concentrations for use as descriptors in chemometric models for the geographical origin determination of the grains.

Linear discriminant analysis was applied to construct classification models for teff grains cultivated in the studied three zones. Accordingly, two discriminant functions were computed, where the first accounted for 69% and the second 31% of the total variance present in the data set. The LDA model provided 96% correct classification, where only 3 out of the 72 samples were incorrectly classified. Two samples from Awi zone were incorrectly classified, one as East Gojjam and the other as West Gojjam, as well as one sample from East Gojjam was incorrectly classified as from West Gojjam. The separation between the geographical origins of teff samples in the discriminant space was examined by plotting two function scores in Fig. 4. It was clearly observed that samples from the three administrative zones were well distinguished from each other, demonstrating that elements contained enough information to identify the zones where the teff samples were grown.

Fig. 4
figure 4

Scatter plot showing the distribution of teff samples on the space created by the scores of the first two discriminant functions

The distribution of the teff samples on the plane of the two functions is shown in Fig. 4. The first function effectively distinguished between East Gojjam and Awi zone teff samples. This function is correlated mainly with Ca and Fe, followed by Mg, Zn, and Mn (Table 6). Teff samples from Awi zone occupy the space to the positive side of the first function, hence characterized by higher contents of Ca and Fe than samples from East Gojjam.The second function distinguished West Gojjam from East Gojjam and Awi zone teff samples. This function is correlated mainly and negatively with Cr and Ni, followed by Cu, K, and Na in Table 6. Teff samples from West Gojjam occupy the space to the positive side of the second function, hence characterized by lower contents of Cr and Ni than samples from East Gojjam and Awi zones. Hence, the Ca, Fe, Mg, and Zn were the main discriminating factors among teff samples from Awi zone Table 6.

Table 6 Canonical discriminant function coefficients of the linear discriminant model constructed with the metal contents of teff grain samples from three administrative zones of Amhara Region, Ethiopia

The reliability of using the determined elements as descriptors with LDA modeling, in developing an analytical method for determining the origin of teff from the three zones, was evaluated from the prediction ability of the model. For this, randomly selected 24 samples, out of the total 72 samples, were used as a validation set, while the remaining 48 samples were used as a training set to construct the LDA model. The analysis was carried out five times by changing the 24 random samples, and the average of the five iterations was calculated. Consequently, the model provided a prediction ability of 91.7%, equivalent to two samples out of the 24 samples whose group membership was incorrectly predicted.

Additionally, leave-one-out cross-validation analysis was performed, where each sample was classified by the LDA model constructed with the remaining sample set. In this, 91.7% of the cross-validated samples were correctly classified.

Regarding the growing districts, application of LDA provided 93.1% correct classification of the teff grain samples into the respective growing districts. Only 5 from the total 72 samples were incorrectly classified. These were, 2 samples from Aneded, one incorrectly classified as Shebelberenta and the other as Gonjikolela, one sample from Shebelberenta, which was incorrectly classified as Gonjikolela, and one sample from Zigam, which was incorrectly classified as Shebelberenta.

Three clusters were evident from the LDA scores plot of teff grains from the six districts (Fig. 5). These are teff grains from Zigam, Enamey and a group of samples from Shebelberenta, Gonjikolela and Adet. Enamey is discriminated from the other districts by the first function, which explained 53% of the data variance. Magnesium is more associated with this function (Table 7). Hence, teff from Enamey is differentiated by its higher content of Mg. Teff from Zigam is discriminated from the remaining districts by the second function, which accounted for 21% of the data variance.The elements Ca, Mn, Fe, and Zn are correlated with the second function (Table 7). Hence, Ca, Mn, Fe, and Zn were good descriptors of teff samples from Zigam district.

Fig. 5
figure 5

Scatter plot of teff samples from six districts on the plane of the first two discriminant functions

Table 7 Correlation between the analyzed elements and discriminant functions

Linear discriminant analysis of teff varieties

The potential of using the elemental concentrations in the construction of discriminant models useful for the authentication of the four varieties of teff (Kuncho, Tseday, Key, Sergegna) was also investigated. From the computed three discriminant functions, the first two accounted for 93% of the variation in the data set.The LDA model provided 93% correct classification of the 72 teff grain samples into the four varities. Where, 4 samples of Tseday variety were incorrectly classified, 3 as Key and one as Kuncho, and one sample of Key variety incorrectly classified as Sergegna. Distinct discrimination between Kuncho and Tseday varieties was observed by the first function (Fig. 6), which accounts for 92.6% of the total variance in elemental concentrations among varieties of teff. The elements Ca, and Fe were the highest contributors to the first function, and Tseday variety was distinguished from Kuncho by its higher content of Ca and Fe (Table 8). Hence, Ca and Fe were the main discriminating factors among the teff samples depending on varieties.

Fig. 6
figure 6

Scatter plot of four varieties of teff samples on the space of the first two discriminant functions

Table 8 Correlation between elements in varietal types of teff samples and discriminant functions


An efficient digestion procedure for the determination of metals in the teff samples was optimized and validated through the spiking method, and a good percentage recovery was obtained from 85.5 to 109% for all the metals of interest. The levels of metals in teff samples were in the order K > Mg > Ca > Fe > Na > Mn > Zn > Cu > Cr > Ni > Co > Cd. One-way ANOVA revealed the presence of significant difference (p < 0.05) in the mean concentrations of metals, except for K, Na, Co, Ni, Cu and Cd,among the sampling zones. With the aid of statistical pattern recognition,metals Ca and Fe, followed by Mg, Zn, and Mn were identified as the most discriminating metals for teff samples from three sampling zones. Ca and Fe were selected as suitable discriminant marker for teff samples originating from Awi zone,as well as varietal types of teff. While, K was the marker element for teff from East Gojjam zone. Application of LDA successfully classified the teff samples in to three administrative zones. The LDA model provided a prediction ability of 91.7%, equivalent to two samples out of the 24 samples whose group membership was incorrectly predicted. Additionally, leave-one-out cross-validation analysis of the entire 72 samples resulted in 91.7% correct classification. Hence, multielement data combined with LDA modeling can be used in the authentication of the geographical origin and varietal types of teff from Amhara region.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.


  1. Bultosa G. Physicochemical characteristics of grain and flour in 13 tef [Eragrostis tef (Zucc.) Trotter] grain varieties. J Appl Sci Res. 2007;3:2042–51.

    CAS  Google Scholar 

  2. Belay G, Tefera H, Getachew A, Assefa K, Metaferia G. Highly client-oriented breeding with farmer participation in the Ethiopian cereal tef [Eragrostis tef (Zucc.) Trotter]. African J Agric Res. 2008;3:022–8.

    Google Scholar 

  3. Gebremariam MM, Zarnkow M, Becker T. Teff (Eragrostis tef) as a raw material for malting, brewing and manufacturing of gluten-free foods and beverages: a review. J Food Sci Technol. 2014;51:2881–95.

    Article  CAS  PubMed  Google Scholar 

  4. Kimani W, Zhang L-M, Wu X-Y, Hao H-Q, Jing H-C. Genome-wide association study reveals that different pathways contribute to grain quality variation in sorghum (Sorghum bicolor). BMC Genomics. 2020;21:1–19.

    Article  Google Scholar 

  5. Girma A, Meareg A. Determination of the iron content of soil and cultivated white teff from parzete and zeghi kebele, debatie woreda metekel zone, benshangul gumuz, Ethiopia. Int J Innovat Pharmaceut Sci Res. 2017;5:11–32.

    Google Scholar 

  6. do Nascimento KdO. Paes SdND, de Oliveira IR, Reis IP, Augusta IM: Teff: suitability for different food applications and as a raw material of gluten-free, a literature review. J Food Nutrit Res. 2018;6:74–81.

    Article  Google Scholar 

  7. Ahmed H, Satheesh N, Dibaba K. Functional, physical and sensory properties of cookies prepared from okara, red teff and wheat flours. Croatian J Food Sci Technol. 2018;10:23–32.

    Article  Google Scholar 

  8. Stich B, Haussmann BI, Pasam R, Bhosale S, Hash CT, Melchinger AE, Parzies HK. Patterns of molecular and phenotypic diversity in pearl millet [Pennisetum glaucum (L.) R. Br.] from West and Central Africa and their relation to geographical and environmental parameters. BMC Plant Biol. 2010;10:1–10.

    Article  Google Scholar 

  9. Kibatu G, Chacha R, Kiende R. Determination of major, minor and trace elements in Tef using portable total x-ray fluorescence (TXRF) spectrometer. EC Nutrition. 2017;9:51–9.

    Google Scholar 

  10. Merga M. Progress, achievements and challenges of tef breeding in Ethiopia. J Agric Sci Food Res. 2018;9:1–8.

    Google Scholar 

  11. Fernández-Cáceres PL, Martín MJ, Pablos F, González AG. Differentiation of tea (Camellia sinensis) varieties and their geographical origin according to their metal content. J Agric Food Chem. 2001;49:4775–9.

    Article  PubMed  Google Scholar 

  12. Villarreal D, Laffargue A, Posada H, Bertrand B, Lashermes P, Dussert S. Genotypic and environmental effects on coffee (Coffea arabica L.) bean fatty acid profile: Impact on variety and origin chemometric determination. J Agric Food Chem. 2009;57:11321–7.

    Article  CAS  PubMed  Google Scholar 

  13. Maione C, Batista BL, Campiglia AD, Barbosa F Jr, Barbosa RM. Classification of geographic origin of rice by data mining and inductively coupled plasma mass spectrometry. Comput Electron Agric. 2016;121:101–7.

    Article  Google Scholar 

  14. Kalinova J, Moudry J. Content and quality of protein in proso millet (Panicum miliaceum L.) varieties. Plant Foods Human Nutr. 2006;61:43–7.

    Article  Google Scholar 

  15. Luo D, Dong H, Luo H, Xian Y, Wan J, Guo X, Wu Y. The application of stable isotope ratio analysis to determine the geographical origin of wheat. Food Chem. 2015;174:197–201.

    Article  CAS  PubMed  Google Scholar 

  16. Wu Y, Luo D, Dong H, Wan J, Luo H, Xian Y, Guo X, Qin F, Han W, Wang L. Geographical origin of cereal grains based on element analyser-stable isotope ratio mass spectrometry (EA-SIRMS). Food Chem. 2015;174:553–7.

    Article  CAS  PubMed  Google Scholar 

  17. Gonzálvez A, Armenta S, De La Guardia M. Geographical traceability of “Arròs de Valencia” rice grain based on mineral element composition. Food Chem. 2011;126:1254–60.

    Article  Google Scholar 

  18. Zhao H, Guo B, Wei Y, Zhang B. Multi-element composition of wheat grain and provenance soil and their potentialities as fingerprints of geographical origin. J Cereal Sci. 2013;57:391–7.

    Article  CAS  Google Scholar 

  19. Endaye M, Atlabachew M, Mehari B, Alemayehu M, Mengistu DA, Kerisew B. Combining multi-element analysis with statistical modeling for tracing the origin of green coffee beans from Amhara region. Ethiopia Biol Trace Element Res. 2020;195:669–78.

    Article  CAS  Google Scholar 

  20. Zhao H, Zhang S, Zhang Z. Relationship between multi-element composition in tea leaves and in provenance soils for geographical traceability. Food Control. 2017;76:82–7.

    Article  CAS  Google Scholar 

  21. Reta C, Atlabachew M, Asmellash T, Hilluf W, Yayinie M, Wubieneh TA. Polyphenol and flavonoid content in major Teff [Eragrostis tef (Zuccagni) Trotte r] varieties in Amhara region. Ethiopia Plos One. 2022;17:e0272010.

    Article  CAS  PubMed  Google Scholar 

  22. Mekassa B, Chandravanshi BS. Levels of selected essential and non-essential metals in seeds of korarima (Aframomum corrorima) cultivated in Ethiopia. Braz J Food Technol. 2015;18:102–11.

    Article  CAS  Google Scholar 

  23. Wang F, Zhao H, Yu C, Tang J, Wu W, Yang Q. Determination of the geographical origin of maize (Zea mays L.) using mineral element fingerprints. J Sci Food Agric. 2020;100:1294–300.

    Article  CAS  PubMed  Google Scholar 

  24. Liang K, Liang S, Lu L, Zhu D, Cheng L. Geographical origin traceability of foxtail millet based on the combination of multi-element and chemical composition analysis. Int J Food Prop. 2018;21:1769–77.

    Article  CAS  Google Scholar 

  25. Gebregewergis A. Levels of selected metals in white teff grain samples collected from there different areas of Ethiopia by using microwave plasma atomic emission spectroscopy (MP-AES). Int J Novel Res Physics Chem Math. 2020;7:13–24.

    Google Scholar 

  26. Habte A, Lelago A, Doda MB, Kanido CK. Determination of selected essential and non-essential metals in three types of teff (Eragrostis tef) Cultivated in Southern Ethiopia. J Agric Chem Environ. 2020;9:255.

    CAS  Google Scholar 

  27. Kebede Z. Levels of essential elements in three teff [Eragrostis tef (Zucc.) Trotter] varieties. Addis Ababa University. 2009.

  28. Dame ZT. Analysis of major and trace elements in teff (Eragrostis tef). J King Saud University Sci. 2020;32:145–8.

    Article  Google Scholar 

  29. Lovis L. Alternatives to wheat flour in baked goods. Cereal Foods World. 2003;48:61.

    Google Scholar 

  30. Yun Y-R, Park S-H. Antioxidant activities of brown teff hydrolysates produced by protease treatment. J Nutr Health. 2018;51:599–606.

    Article  CAS  Google Scholar 

  31. Abebe A, Chandravanshi BS. Levels of essential and non-essential metals in the raw seeds and processed food (roasted seeds and bread) of maize/corn (Zea mays L.) cultivated in selected areas of Ethiopia. Bull Chem Soc Ethiopia. 2017;31:185–99.

    Article  CAS  Google Scholar 

  32. Atlabachew M, Chandravan BS, Redi M. Profile of major, minor and toxic metals in soil and Khat (Catha edulis Forsk) cultivars in Ethiopia. Trends Appl Sci Res. 2011;6:640–55.

    Article  CAS  Google Scholar 

  33. Martínez-Ballesta M, Dominguez-Perles R, Moreno D, Muries B, Alcaraz-López C, Bastías E, García-Viguera C, Carvajal M. Minerals in plant food: effect of agricultural practices and role in human health. A Rev Agron Sustain Develop. 2010;30:295–309.

    Article  Google Scholar 

  34. Ncube B, Finnie J, Van Staden J. Quality from the field: the impact of environmental factors as quality determinants in medicinal plants. S Afr J Bot. 2012;82:11–20.

    Article  Google Scholar 

  35. Gizaw B, ZerihunTsegay GT, Aynalem E, Abatneh E, Amsalu G. Traditional knowledge on teff (Eragrostis tef) farming practice and role of crop rotation to enrich plant growth promoting microbes for soil fertility in East Showa: Ethiopia. Agric Res Technol. 2018;16:39–55.

    Google Scholar 

  36. Mihretie FA, Tsunekawa A, Haregeweyn N, Adgo E, Tsubo M, Masunaga T, Meshesha DT, Ebabu K, Nigussie Z, Sato S. Exploring teff yield variability related with farm management and soil property in contrasting agro-ecologies in Ethiopia. Agric Syst. 2022;196:103338.

    Article  Google Scholar 

  37. Abewa A, Adgo E, Yitaferu B, Alemayehu G, Assefa K, Solomon JK, Payne W. Teff grain physical and chemical quality responses to soil physicochemical properties and the environment. Agronomy. 2019;9:283.

    Article  CAS  Google Scholar 

  38. Baye K, Mouquet-Rivier C, Icard-Vernière C, Picq C, Guyot JP. Changes in mineral absorption inhibitors consequent to fermentation of E thiopian injera: implications for predicted iron bioavailability and bioaccessibility. Int J Food Sci Technol. 2014;49:174–80.

    Article  CAS  Google Scholar 

  39. Gebregewergis A, Chandravanshi BS, Redi-Abshiro M. levels of selected metals in teff grain samples collected from three different areas of Ethiopia by microwave plasma-atomic emission spectroscopy. Bull Chem Soc Ethiop. 2020;34:449–62.

    Article  CAS  Google Scholar 

  40. Umeta M, West CE, Fufa H. Content of zinc, iron, calcium and their absorption inhibitors in foods commonly consumed in Ethiopia. J Food Compos Anal. 2005;18:803–17.

    Article  CAS  Google Scholar 

Download references


The authors are thankful to the College of Science, Bahir Dar University, Ethiopia, for giving access to its laboratory facilities. We would also like to thank the Bureau of Agriculture of the Amhara Region, Ethiopia, as well as its district and sub-district offices, for giving us information for this work.


Not applicable.

Author information

Authors and Affiliations



MA and TA performed the conceptualization and supervision of the study. MA and BM carried out the data analysis and interpretation. CR did the laboratory work and drafted the manuscript. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Minaleshewa Atlabachew or Bewketu Mehari.

Ethics declarations

Ethics approval and consent to participate

This research on Eragrostis tef (Zucc.) Trotter complies with the research ethics guidelines of Bahir Dar University, Ethiopia. The grain samples were collected directly from private land (farmers farmland) upon permission granted by the farmland owners. The plants were identified by Dr. Ali Seid and specimens are deposited in the Ethiopian Herbarium with voucher numbers Ali and Minale 001 Eragrostis tef (Zucc.) Trotter Var. Name Tsedey, Ali and Minale 002 Eragrostis tef (Zucc.) Trotter Var. Name Kuncho, Ali and Minale 003 Eragrostis tef (Zucc.) Trotter Var. Name Sergena and Ali and Minale 004 Eragrostis tef (Zucc.) Trotter Var. Name Qey tef.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Reta, C., Asmellash, T., Atlabachew, M. et al. Multielement analysis coupled with chemometrics modelling for geographical origin classification of teff [Eragrostis tef (Zuccagni) Trotter] grains from Amhara Region, Ethiopia. BMC Chemistry 17, 50 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Teff
  • Multielement
  • Geographical origin
  • Multivariate
  • Amhara
  • Ethiopia