Comparative Morphometry of the Sitophilus zeamais Motschulsky Populations (Coleoptera: Curculionidae) Infesting Yellow Maize in Four West African Countries

Sitophilus zeamais is a primary pest of stored maize, which can cause quality and quantity losses. Maize is cultivated in several West African countries where different agro-climatic and agro-ecological conditions exist. These conditions could influence the Sitophilus zeamais morphology. Thus, the purpose of this study is to identify the Sitophilus zeamais morphological groups in the four countries. This was carried out through the morphometric study of 30 male Sitophilus zeamais individuals taken from the yellow maize in each country, by measuring with a binocular magnifying glass, 16 variables of the different tagmas belonging to the insect. Subsequently, statistical analyses were firstly carried out with the raw measurements, such as the principal component analysis, in order to extract the contribution of the variables on each factorial axis and to visualise the correlations between the variables. Taking the case of the transformed measurements, in addition to the principal component analysis, a discriminant factorial analysis was carried out to find out the best factor for gathering individuals according to their morphology, a confusion matrix to summarise the individual’s reclassification in order to deduce the good and bad classification rates, and an ascending hierarchical classification which consists of gathering the individuals according to their similarity and separating them according to their dissimilarity using the truncation method. Therefore, five morphological groups were identified. However, the existence of these morphological groups would not be linked to agro-ecological, climatic conditions and to the geographical remoteness between countries. How to cite this paper: Sarr, A.G.R.J., Dia, C.A.K.M., Diome, T. and Sembene, P.M. (2022) Comparative Morphometry of the Sitophilus zeamais Motschulsky Populations (Coleoptera: Curculionidae) Infesting Yellow Maize in Four West African Countries. Advances in Entomology, 10, 85-98. https://doi.org/10.4236/ae.2022.101007 Received: September 23, 2021 Accepted: January 3, 2022 Published: January 6, 2022 Copyright © 2022 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access A. G. R. J. Sarr et al. DOI: 10.4236/ae.2022.101007 86 Advances in Entomology


Introduction
Maize (Zea mays, Poaceae family) is a cereal cultivated alone or with several crops in various agro-ecological zones [1]. Its production enables people to survive the hunger periods in the rural areas. However, maize is infested in the field and during its storage by insect pests such as the Sitophilus zeamais (Motschulsky). This primary pest of stored maize can cause qualitative and quantitative losses [2] and thus encourage the colonisation of stocks by secondary pests and toxin-producing fungi such as the Aspergillus flavus [3]. These various types of damage reduce the quality of the grain and make it unfit for consumption [4].
However, estimating the losses caused by the S. zeamais is most of the time difficult because of its high variability, which can be explained by the agro-ecological conditions, the harvest period [2] [5] and even by the maize variety [2] [6]. Sitophilus zeamais is a cosmopolitan insect that can adapt itself to high temperatures and humidity [5]. It is widely distributed in tropical regions where it attacks standing crops [7]. In Africa, it is especially widespread on maize [5].
However, we noted that some African countries have different agro-climatic and agro-ecological conditions, while biodemographic profiles affect the biological strength of the insect [8]. Furthermore, the environmental factors, including agronomic practices (types of fertilisers, cultivation practices), influence the grain composition [9], the morphology and even the genetics of the insect population that infests the maize, as noted in the case of Caryedon serratus [10]. This context led us to ask whether the S. zeamais has the same morphology in the different West African agro-ecological areas where maize is cultivated in important quantities. In order to answer this scientific issue, the objective of this study will be to determine the existence of the S. zeamais morphometric groups in West Africa.

Sampling
The study areas were chosen according to the maize-producing countries, the infestation and taking into account the geographical location. Consequently, the sampling was carried out in four West African countries (Table 1). In our sampling strategy, maize was collected from the storage warehouses in those countries and kept in plastic bags. They were taken to the insectarium of the entomology and acarology laboratory of the Animal Biology Department of the in test tubes, which were specific to each population. Each individual of a sample was coded, using the first letter of the genus name in uppercase followed by the first two letters of the country of origin (the first letter in uppercase and the second in lowercase), the colour of the maize type and finally a serial number.
For example, the SSeJ1 code represents a Sitophilus zeamais individual from Senegal, on yellow (Jaune) maize, with order number 1 ( Table 2).

Morphometric Study
The biological material used for the morphometry was composed of male Sitophilus zeamais taken from the maize mass rearing. The choice of males for the morphometric study was explained by the fact that females can change their morphologies and as a result bias the measurements. The male identification compared to the female was based on criteria established by [5] study dealing with the insects found in foodstuffs. The male's rostrum was shorter and thicker, but less smooth than the female's rostrum; the ventral face of the abdomen was concave to the apex of the male, approximately flat in the female. The insect was removed from the alcoholic solution using pliers, placed in distilled water and then transferred on paper towel to be dried (the paper sucked up the water in which the insect was soaked). It was then placed on a graph paper in a binocular magnifying glass to proceed to the measurements. The insect was first dissected using pliers by separating its three tagmas: head, thorax and abdomen. Each tagma was then separated from its different articles. A total of sixteen (16) variables were measured with a reasonable degree of accuracy (±1 mm) ( Table 3) ( Figure 1).

Transformed Measurements
Size effect and data transformation The size effect is explained by a correlation circle that gathers all the variables in a single plan for a given axis. This is an unwanted effect, which metric studies are trying to escape [11]. The principle is therefore to eliminate this effect and to bring back all individuals to the same size, so that only shape differences can be observed on the PCA. The removal of the size effect in this study was carried out -For each individual, we computed the average of all the log-transformed variables. We consider that this average score gives a good idea regarding the "size" for this individual; -Finally, for each individual, we subtracted from each of these raw measurements the average size obtained with the log-transformed data.
Therefore, the size effect was eliminated and only the difference of shape will be observed on the PCA. The elimination of the factor (Size) leads to a decrease in the total discrimination. This transformation was carried out using Excel 2010 version.

Discriminant factor analysis
A discriminant factor analysis (DFA) of the populations with the variables transformed data was carried out using the software R version 3.2.3. It was performed with the variables whose contribution value is the same or bigger than (1/P) × 100. Its purpose is to find out the best factor in order to gather individuals according to their morphology.

Determination of variables discriminating the populations
In order to determine the variables that significantly differentiate the populations morphology, a Shapiro-Wilk normality data test was firstly performed on each transformed log-variable. When the p-value of the test is less than 0.05 (the case of all transformed log-variables in this study), the data of the corresponding transformed-log do not follow the normal law and the Kruskal-Wallis test (non-parametric) was performed for a multiple average rank's comparison of the corresponding variable related to the populations. The Wilcoxon test (non-parametric) coupled with the Bonferroni adjustment was performed to test the average variations significance rank regarding the transformed log-variables per pair of populations.

Confusion matrix for cross-validation results
The confusion matrix summarizes the reclassifications of the individuals to deduce the good and bad classifications rates. This will allow to determine the "correct percentage (%)", which is the ratio between the number of well-ranked individuals and the total number of individuals. The cross-validation is applied on the log-transformed data using the XLSTAT program version 2016.03.30882 (Addinsoft, 2016).

Ascending hierarchical classification (AHC)
The ascending hierarchical classification aims to identify a posteriori classes that are not obvious a priori. It consists in gathering individuals, using the truncation method, with respect to their similarity and separating them according to their dissimilarity. Its purpose is to identify morphologically homogeneous groups from the variables point of view. It was performed with XLSTAT program version 2016. 03

Choice of the Number of Axes to Be Retained
In order to have the number of axes, we accordingly used the elbow criterion.
The purpose of this method is to obtain the maximum kept inertia with the minimum number of factors. In our research, two axes were selected, because the cut-off (the elbow) is found at K = 2 ( Figure 2).

Principal Component Analysis
The

Choice of the Number of Axes to Be Retained after Data Transformation
After the data transformation, the chosen number of axes obtained using the elbow criterion is the same as before the transformation (Figure 4).

Principal Component Analysis of Transformed Variables
The Principal Component Analysis (PCA) revealed after data transformation that the first factorial axis lost 8.08% of its percentage of inertia. On the other hand, the percentage of inertia of the second factorial axis increased up to 3.46%. As a result, the percentage of inertia of the first two axis was equivalent to 39.17%. The best representation quality is always obtained with these first two axis. On the first factorial axis, with a percentage of inertia equivalent to 22.80%, the contribution of some variables is lower than the values obtained before transformation. This is noted with the length of the wing ( The removal of the size effect was materialised by the dispersion of the variables on both sides of the two axes ( Figure 5).

Discrimination after Transformation
The best representation quality is always obtained with the diagram composed of the first two axes (Dim 1 and Dim 2) with 45.76% of inertia. Following the first axis for a percentage of inertia equivalent to 26.37%, two morphometric groups were identified with more or less important areas of introgression. On the one hand, they are individuals from the population of Guinea Conakry, and on the other hand individuals from Senegal, Mali and Ghana. According to the second factorial axis, for a percentage of inertia equivalent to 19.39%, individuals from the population of Ghana are different from individuals belonging to the population of Senegal and Mali, with a zone of introgression between them ( Figure 6).

Variables Allowing Discrimination after Transformation
According to the first axis (Dim), the morphological difference between the individuals belonging to the population of Senegal compared to those belonging to the population of Guinea

Confusion Matrix for the Cross-Validation of Populations
The confusion matrix summarizes the reclassifications of individuals to deduce the good and bad classifications rates. This will allow the determination of the  (Table 5).

Discussion
According to the elbow criterion, two factorial axes were selected for raw measurements and transformed data. These two factorial axes reveal a greater variability in the morphometric properties for the raw measurements compared to the transformed measurements. This difference noted in the percentage of inertia reflects the elimination of the size effect between them. After the elimination of the size effect, individuals from Senegal, Mali and Ghana reveal similarities between them and dissimilarity between the individuals from Guinea Conakry.
However, similarities are noted between individuals from geographically remote countries with different agro-ecological conditions. Geographical remoteness between countries does not seem to be a determining factor in the similarities, as noted in the morphometric identification study of the Sudano-Sahelian population of Caryedon serratus [16]. Therefore, the similarities and dissimilarities that Ghana is a maize exporting country [17]. There are few or no individuals of S. zeamais from other countries. As for the population of Guinea Conakry, although they are maize importing country [17], their homogeneity could be explained by the increase of local maize production marked by the predominance of family farming [18]. Individuals of S. zeamais are originated from this local maize. The heterogeneity of the individuals belonging to Mali and Senegal would come from the importation of one part of the maize. These two countries respectively import 17% and 12% of their maize [19]. An ascending hierarchical classification study of individuals revealed a similarity between certain populations. Five morphometric groups were identified, with the individuals of Guinea having the most varied morphology. However, the similarities noted between certain individuals confirm only the morphological similarities that were identified between certain populations in the discriminant factor analysis of the populations. This is the case of the individuals from Senegal, Ghana and Mali, which are found in morphometric group 4, but morphometric groups 1 and 2 are respectively specific to the populations of Guinea Conakry and Senegal. The existence of these specific morphometric groups could be explained either by the agro-ecological conditions or by the storage conditions in those countries, knowing that according to [20] in Taking the case of Senegal, infestation on stored maize begins during the drying of the ears in the field [21]. In addition, the storage means are mostly rudimentary and aeration system is not respected [21], promoting the multiplication of insects. Moreover, the same warehouses are used upon each harvest operation. The new individuals are generated by the individuals that firstly infested the stocks, which would make the reproduction system closed to other individuals. This reproduction system could be at the origin of this morphometric group that is specific to Senegal.

Conclusion
The determination of S. zeamais morphometric groups in West Africa is of the utmost importance. It enables to find out whether the morphology of this insect depends on its geographical origin and to develop traceability of this insect pest for a better sustainable fighting strategy. Discriminant factorial analysis revealed the existence of morphometric groups that would not be linked to agro-ecological and climatic conditions and to the geographical remoteness of countries. As a matter of fact, a genetic study is needed to find out if the morphometric differences detected are observable at the molecular level using mitochondrial genes.