ABSTRACT.

The Spodoptera frugiperda, is one of the most deleterious pests of popcorn and the identification of resistant genotypes is determinant in breeding programs. The objective of this study was to select popcorn genotypes resistant to S. frugiperda and the key traits related to the identification of resistance. The popcorn varieties UEM J1, Composto Márcia, Arachida, Composto Gaúcho, and Zapalote Chico (resistant check) were evaluated in a completely randomized design with 100 replications. The experimental unit consisted of one Petri dish, containing plant material and a caterpillar The following traits were evaluated: larval stage duration (LSt), food intake weight(IW), final larva weight (FW), mean larva weight (MW), feces (F), assimilated (A) and metabolized food weight (M), relative consumption rate (RCR), relative metabolic rate (RMR), relative growth rate (RGR), conversion efficiency of ingested food (CEI), apparent digestibility (AD), conversion efficiency of digested food (CED), and leaf area consumed (LAC). The diagnosis of multicollinearity, analysis of canonical variables, genetic divergence, hierarchical clustering, factor analysis and canonical correspondence analysis were carried out to perform multivariate analysis. After the multicollinearity test, the traits FW, IW, RCR, AD, and LAC were maintained for further analysis. The traits IW, FW and AD were determinant in the resistance by antixenosis expressed by the varieties Zapalote Chico and Arachida, for the varieties Composto Gaúcho and Composto Márcia the determining characteristics were RCR and LAC and for the variety UEM J1 the variable LAC showed greater importance. Variety Arachida was considered resistant to S. frugiperda by antixenosis and can be used in the future as a source of favorable alleles to breed resistant popcorn hybrids. The traits relative consumption rate, apparent digestibility and leaf area consumed were considered key traits in the identification of resistance against S. frugiperda in popcorn genotypes.

Keywords:
fall armyworm; genetic resistance; Zea mays

Introduction

The fall armyworm, Spodoptera frugiperda, is considered one of the most deleterious pests of maize in all maize-producing regions (Juárez et al., 2012Juárez, M. L., Murúa, M. G., García, M. G., Ontivero, M., Vera, M. T., Vilardi, J. C., ... Willink, E. (2012). Host association of Spodoptera frugiperda (Lepidoptera: Noctuidae) corn and rice strains in Argentina, Brazil, and Paraguay. Journal of Economic Entomology, 105(2), 573-582. DOI: https://doi.org/10.1603/EC11184
https://doi.org/https://doi.org/10.1603/... , He et al., 2021He, L. M., Wang, T. L., Chen, Y. G., Ge, S. S., Wyckhuys, K. A. G., & Wu, K. M. (2021). Larval diet affects development and reproduction of East Asian strain of the fall armyworm, Spodoptera frugiperda. Journal of Integrative Agriculture, 20(3), 736-744. DOI: https://doi.org/10.1016/S2095-3119(19)62879-0
https://doi.org/https://doi.org/10.1016/... ). In the tropics, it causes significant economic damage because it feeds on maize from the seedling to the reproductive stage (Overton et al., 2021Overton, K., Maino, J. L., Day, R., Umina, P. A., Bett, B., Carnovale, D., … Reynolds, O. L. (2021). Global crop impacts, yield losses and action thresholds for fall armyworm (Spodoptera frugiperda): A review. Crop Protection, 145, 105641. DOI: https://doi.org/10.1016/j.cropro.2021.105641
https://doi.org/https://doi.org/10.1016/... ). In popcorn, which has been less the focus of breeding efforts than field corn, the damage tends to be more severe, since no resistant or tolerant genotypes are available on the market. Considering the increasing nationwide and worldwide demand for and importance of popcorn (Hebei Cici Co., Ltd., 2018Hebei Cici Co., Ltd. (2018). Popcorn market by type (microwave popcorn and ready-to-eat popcorn) and end user (household and commercial) - Global opportunity analysis and industry forecast, 2017-2023. Retrieved on Mar. 10, 2022 from 10, 2022 from https://www.indiampopcorn.com/news/popcorn-market-by-type-microwave-popcorn-and-ready-to-eat-popcorn-and-end-user-household-and-commercial-global-opportunity-analysis-and-industry-forecast-2017-2023/
https://www.indiampopcorn.com/news/popco... ; Caranhato et al., 2022Caranhato, A. L. H., Trindade, R. W. R., Uhdre, R. S., Pinto, R. J. B., Scapim, C. A., & Paterniani, M. E. A. G. Z. (2022). Genetic improvement in popcorn. Revista Brasileira de Milho e Sorgo, 21, 1-25. DOI: https://doi.org/10.18512/rbms2022v21e1272
https://doi.org/https://doi.org/10.18512... ; Stipp et al., 2024Stipp, O. J., Junior, O. P., Rossi, E., Rosa, J. C., Uhdre, R. S., Rizzardi, D. A., ... Pinto, R. J. B. (2024). Agricultural traits and popping expansion of the popcorn hybrid IAC 125 under different plant densities and irrigation water depth levels. Acta Scientiarum. Agronomy , 46(1), e62929. DOI: https://doi.org/10.4025/actasciagron.v46i1.62929
https://doi.org/https://doi.org/10.4025/... ), the possibility of control by the development of resistant genotypes must be investigated and improved, since the number of insecticide applications has increased considerably over the years (Omoto et al., 2016Omoto, C., Bernardi, O., Salmeron, E., Sorgatto, R. J., Dourado, P. M., Crivellari, A., ... Head, G. P. (2016). Field‐evolved resistance to Cry1Ab maize by Spodoptera frugiperda in Brazil. Pest Management Science, 72(9), 1727-1736. DOI: https://doi.org/10.1002/ps.4201
https://doi.org/https://doi.org/10.1002/... ), resulting in excessive costs and increased environmental risks.

In view of the concern to reduce or even eliminate insecticide applications and only a few of studies involving S. frugiperda resistance in popcorn, breeding programs ought to intensify studies to the selection of resistant and tolerant popcorn genotypes, as well as to elucidate key traits related to the identification of resistance (Oliveira et al., 2018Oliveira, N. C. D., Suzukawa, A. K., Pereira, C. B., Santos, H. V., Hanel, A., Albuquerque, F. A. D., & Scapim, C. A. (2018). Popcorn genotypes resistance to fallarmy worm. Ciência Rural, 48(2), 1-6. DOI: https://doi.org/10.1590/0103-8478cr20170378
https://doi.org/https://doi.org/10.1590/... ; Paiva et al., 2016Paiva, L. A., Corrêa, F., Silva, C. L. T., de Lima Moura, T., Silva, F. C., Silva Araújo, M., & Jesus, F. G. (2016). Resistance of corn genotypes to fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae). African Journal of Biotechnology, 15(35), 1877-1882. DOI: https://doi.org/10.5897/AJB2016.15399
https://doi.org/https://doi.org/10.5897/... ).

Although some traits related to identification of S. frugiperda resistance, e.g., larval stage duration, final and mean larva weight, consumption rates and leaf area consumed have been studied in maize breeding programs (Oliveira et al., 2018Oliveira, N. C. D., Suzukawa, A. K., Pereira, C. B., Santos, H. V., Hanel, A., Albuquerque, F. A. D., & Scapim, C. A. (2018). Popcorn genotypes resistance to fallarmy worm. Ciência Rural, 48(2), 1-6. DOI: https://doi.org/10.1590/0103-8478cr20170378
https://doi.org/https://doi.org/10.1590/... ), they have not yet been conclusively proven and used as key traits. A likely reason may be that the trials were analyzed based on a univariate and inconclusive approach. The use of multivariate techniques is required, which can be an adequate and more efficient tool in the analysis of resistance-related data in popcorn (Aaliya et al., 2016Aaliya, K., Qamar, Z., Ahmad, N. I., Ali, Q., Munim, F. A., & Husnain, T. (2016). Transformation, evaluation of gtgene and multivariate genetic analysis for morpho-physiological and yield attributing traits in Zea mays. Genetika, 48(1), 423-433. DOI: http://doi.org/10.2298/GENSR1601423A
https://doi.org/http://doi.org/10.2298/G... ; Veturi et al., 2012Veturi, Y., Kump, K., Walsh, E., Ott, O., Poland, J., Kolkman, J. M., ... Wisser, R. J. (2012). Multivariate mixed linear model analysis of longitudinal data: An information-rich statistical technique for analyzing plant disease resistance. Phytopathology, 102(11), 1016-1025. DOI: https://doi.org/10.1094/PHYTO-10-11-0268
https://doi.org/https://doi.org/10.1094/... ).

Resistance is the set of physical, chemical and morphological traits that will negatively affect the insect’s oviposition and feeding behavior, while tolerance is the set of traits that will cause plants to withstand the attack of insects without substantial reductions in productivity in comparison with other susceptible genotypes (Aaliya et al., 2016Aaliya, K., Qamar, Z., Ahmad, N. I., Ali, Q., Munim, F. A., & Husnain, T. (2016). Transformation, evaluation of gtgene and multivariate genetic analysis for morpho-physiological and yield attributing traits in Zea mays. Genetika, 48(1), 423-433. DOI: http://doi.org/10.2298/GENSR1601423A
https://doi.org/http://doi.org/10.2298/G... ; Veturi et al.,2012Veturi, Y., Kump, K., Walsh, E., Ott, O., Poland, J., Kolkman, J. M., ... Wisser, R. J. (2012). Multivariate mixed linear model analysis of longitudinal data: An information-rich statistical technique for analyzing plant disease resistance. Phytopathology, 102(11), 1016-1025. DOI: https://doi.org/10.1094/PHYTO-10-11-0268
https://doi.org/https://doi.org/10.1094/... ).

In this study, the hypothesis was proposed that resistance to S. frugiperda can be found in popcorn germplasm from tropical regions and that multivariate analysis can discriminate the main traits related to the identification of resistance. The objective of this study was to select popcorn genotypes resistant to S. frugiperda and the key traits related to the identification of resistance.

Material and methods

The trial was carried out in Maringá, Paraná State, Brazil, to evaluate five popcorn varieties (Table 1) for resistance to fall armyworm (Spodoptera frugiperda), as well as the key traits related to the identification of this resistance.

The popcorn genotypes used in the trial consisted of plants of the varieties UEM J1, Composto Márcia, Arachida, Composto Gaúcho, and Zapalote Chico. The first four evaluated varieties were developed by the Specialty Corn Breeding program of the State University of Maringá. The variety used as check, Zapalote Chico, was introduced from Central América and has resistance against fall armyworm on the basis of antixenosis (alimentary avoidance) and antibiosis (lower insect survival after feeding on host tissue) (Cabrera-Toledo, Carballo-Carballo, Mejía-Contreras, García-De los Santos, & Vaquera-Huerta, 2019Cabrera-Toledo, J. M., Carballo-Carballo, A., Mejía-Contreras, J. A., García-de los Santos, G., & Vaquera-Huerta, H. (2019). Caracterización de poblaciones sobresalientes de maíz de la raza zapalote chico. Revista Fitotecnia Mexicana, 42(3), 269-279. DOI: https://doi.org/10.35196/rfm.2019.3.269
https://doi.org/https://doi.org/10.35196... ).

The plants were grown in a greenhouse with automatic irrigation. Crop management was carried out in accordance with the recommendations for corn culture (Oliveira et al., 2018Oliveira, N. C. D., Suzukawa, A. K., Pereira, C. B., Santos, H. V., Hanel, A., Albuquerque, F. A. D., & Scapim, C. A. (2018). Popcorn genotypes resistance to fallarmy worm. Ciência Rural, 48(2), 1-6. DOI: https://doi.org/10.1590/0103-8478cr20170378
https://doi.org/https://doi.org/10.1590/... ). The seeds were separated and three sown in each pot, which contained soil and substrate (3:1). After sowing, the pots were irrigated daily and side-dressed with urea (45% N). The other cultural treatments were applied as required for full crop development, without using any other chemical product, so as not to affect larva growth.

The corn leaves used to feed the caterpillar were collected when the plants were in the eight-leaf (V8) stage, so that all varieties were evaluated when the plants were in the same developmental stage (Oliveira et al., 2018Oliveira, N. C. D., Suzukawa, A. K., Pereira, C. B., Santos, H. V., Hanel, A., Albuquerque, F. A. D., & Scapim, C. A. (2018). Popcorn genotypes resistance to fallarmy worm. Ciência Rural, 48(2), 1-6. DOI: https://doi.org/10.1590/0103-8478cr20170378
https://doi.org/https://doi.org/10.1590/... ; Paiva et al., 2016Paiva, L. A., Corrêa, F., Silva, C. L. T., de Lima Moura, T., Silva, F. C., Silva Araújo, M., & Jesus, F. G. (2016). Resistance of corn genotypes to fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae). African Journal of Biotechnology, 15(35), 1877-1882. DOI: https://doi.org/10.5897/AJB2016.15399
https://doi.org/https://doi.org/10.5897/... ; Sanches et al., 2019Sanches, R. E., Suzukawa, A. K., Contreras-Soto, R. I., Rizzardi, D. A., Kuki, M. C., Zeffa, D. M., ... Scapim, C. A. (2019). Multivariate analysis reveals key traits of fall armyworm resistance in tropical popcorn genotypes. Bragantia, 78(2), 175-182. DOI: https://doi.org/10.1590/1678-4499.20180230
https://doi.org/https://doi.org/10.1590/... ).

The insects required to initiate the trial were hatched from S. frugiperda eggs donated by Embrapa Soybean, in Londrina, Paraná State, Brazil, and from eggs collected in corn fields on the Experimental Farm of Iguatemi - Maringá, Paraná State, Brazil. The larvae hatched from these eggs were fed an artificial diet and three generations were reared to be used in our trial.

The laboratory trial was carried out in a completely randomized design, with five treatments and 100 replications in an air-conditioned chamber at 25 ± 1°C, air humidity of 70 ± 10%, and a 12h photoperiod. Each experimental unit consisted of a sterile acrylic Petri dish (diameter 9.0 cm, height 1.5 cm), lined with filter paper moistened with distilled water to maintain the leaf turgor, containing only one larva per dish, to avoid insect cannibalism. Each treatment consisted of three Petri dishes with moist filter paper and plant material, to calculate the water loss.

After hatching of the fourth larva generation raised on artificial diet, they were distributed in separate Petri dishes. The filter paper was changed, plant material was supplied to feed the larvae and the biological parameters were evaluated daily. Plant material of the same variety was continuously supplied until the end of the larval stage.

The traits were evaluated as proposed by Waldbauer (1968Waldbauer, G. P. (1968). The consumption and utilization of food by insects. Advances in Insect Physiology, 5, 229-288. DOI: https://doi.org/10.1016/S0065-2806(08)60230-1
https://doi.org/https://doi.org/10.1016/... ) with changes made by Scriber and Slansky Jr. (1981Scriber, J. M., & Slansky Jr., F. (1981). The nutritional ecology of immature insects. Annual Review of Entomology, 26(1), 183-211. DOI: https://doi.org/10.1146/annurev.en.26.010181.001151
https://doi.org/https://doi.org/10.1146/... ) as follows:

Larval stage duration: LSt

Food intake weight: IW

Final larva weight: FW

Mean larva weight: MW

Time period (days): T

Feces: F

Assimilated food: A = IW - F

Metabolized food: M = A - FW

Relative consumption rate: RCR = IW/(MW * T)

Relative metabolic rate: RMR = M (MW * T)

Relative growth rate: RGR = FW/ (MW * T)

Conversion efficiency of ingested food: CEI = (FW/IW) * 100

Apparent digestibility: AD = ((A-F)/A) * 100

Conversion efficiency of digested food: CED = FW/A * 100

Leaf area consumed: LAC = IW/SDm, where SDm: Mean surface density.

Food consumption and use were evaluated daily throughout the larval stage, by weighing the fresh food weight, leftover food weight, feces weight and larva weight. The values were measured for each larva and the mean for each repetition was calculated according to the number of live larvae on the day of the evaluation. The final data were obtained from the mean of the 100 replications per treatment.

The larva weight was measured directly by individual weighing during the whole larval stage. The daily weighing of the larvae was only initiated the 5^d after egg hatch. Food intake weight (IW) was calculated indirectly, by subtracting of the statistically corrected leftover food (Lc) weight from the weight of the supplied food (SF) on the day before.

The excreta were collected and weighed individually for each larva during the entire larval stage to obtain the total weight of feces produced (F). The leaf area consumed (LAC) was calculated indirectly from the relationship between food intake (IW) and mean surface density (SDm).

The assumptions of normality of residuals and homogeneity of variances were evaluated by the Shapiro-Wilk and Levene tests, respectively. The data were processed statistically by multivariate analysis of variance (MANOVA), using the statistical software Genes (Cruz, 2016Cruz, C. D. (2016). Genes Software-extended and integrated with the R, Matlab and Selegen. Acta Scientiarum. Agronomy, 38(4), 547-552. DOI: https://doi.org/10.4025/actasciagron.v38i4.32629
https://doi.org/https://doi.org/10.4025/... ) integrated with R software (R Core Team, 2018R Core Team (2013). R: A language and environment for statistical computing. Vienna, AT: R Foundation for Statistical Computing.).

The Variance Inflation Factor (VIF) and the Condition Index (CI) were used as criteria to assess the degree of multicollinearity between the predictive traits. Variance inflation factors of >10 are generally considered evidence of substantial multicollinearity and normally the reason for the removal of certain predictors. In addition, multicollinearity is considered weak when the CI is less than 100 (Montgomery, Peck, & Vining, 2001Montgomery, D. C., Peck, E. A., & Vining, G. G. (2021). Introduction to linear regression analysis. New Jersey, NY: John Wiley & Sons.).

Genetic divergence between the genotypes was evaluated by UPGMA (Unweighted Pair-Group Method using Arithmetic Averages), based on Mahalanobis distance. The groups were established according to the methodology proposed by Mojena (1977Mojena, R. (1977). Hierarchical grouping methods and stopping rules: an evaluation. The Computer Journal, 20(4), 359-363. DOI: https://doi.org/10.1093/comjnl/20.4.359
https://doi.org/https://doi.org/10.1093/... ). Then, canonical variable analysis was carried out with clustering by the Tocher method, factor analysis and later canonical correspondence analysis (Friendly & Sigal, 2017Friendly, M., & Sigal, M. (2017). Graphical methods for multivariate linear models in psychological research: An R tutorial. The Quantitative Methods for Psychology, 13(1), 20-45. DOI: https://doi.org/10.20982/tqmp.13.1.p020
https://doi.org/https://doi.org/10.20982... ).

Factor analysis was performed considering all evaluated traits (Johnson & Wichern, 2002Johnson, R. A., & Wichern, D. W. (2002). Applied multivariate statistical analysis (v. 5, n. 8). Upper Saddle River, NJ: Prentice hall.). The factorial loads were extracted by the principal component method, and the factors established by varimax rotation. In this study, factor loads above 0.90 as well as the highest values for the community were considered (Keith & Reynolds, 2018Keith, T. Z., & Reynolds, M. R. (2018). Using confirmatory factor analysis to aid in understanding the constructs measured by intelligence tests. In D. P. Flanagan, & E. M. McDonough (Eds.), Contemporary intellectual assessment: Theories, tests, and issues (p. 853-900). New York, NY: Guilford Press). Canonical correspondence analysis was performed as described by Greenacre (1984Greenacre, M. J. (1984). Theory and applications of correspondence analysis. London, UK: Academic Press.).

The analyses were performed using the statistical software Genes (Cruz, 2016Cruz, C. D. (2016). Genes Software-extended and integrated with the R, Matlab and Selegen. Acta Scientiarum. Agronomy, 38(4), 547-552. DOI: https://doi.org/10.4025/actasciagron.v38i4.32629
https://doi.org/https://doi.org/10.4025/... ) and SAS (SAS Institute, 2015SAS Institute. (2015). Base SAS 9.2 Procedures Guide. Cary, NC: SAS Institute Inc.) at 1% probability.

Results and discussion

Normality (p > 0.01) and homogeneity of variances (p > 0.01) were reported for all evaluated traits. The MANOVA test showed significant differences between the mean vectors of the genotypes for all evaluated traits (p < 0.01), indicating the existence of genetic variability.

Multicollinearity among traits was assessed by the criteria VIF and CI. Traits with VIF and CI values greater than 10 and 100, respectively, are generally considered evidence of substantial multicollinearity among variables and makes the removal of these predictors necessary (Prunier, Colyn, Legendre, Nimon, & Flamand, 2015Prunier, J. G., Colyn, M., Legendre, X., Nimon, K. F., & Flamand, M. C. (2015). Multicollinearity in spatial genetics: separating the wheat from the chaff using commonality analyses. Molecular Ecology, 24(2), 263-283. DOI: https://doi.org/10.1111/mec.13029
https://doi.org/https://doi.org/10.1111/... ). The traits LSt, MW, F, A, M, RMR, RGR, CEI, and CED had high VIF and CI values and were therefore eliminated. On the other hand, the traits FW, IW, RCR, AD, and LAC had VIF and CI had values below 10 and 100, respectively, and were therefore maintained for subsequent analyses (Table 2).

Multivariate techniques to select traits related to S. frugiperda resistance in popcorn genotypes were also used by Oliveira et al. (2018Oliveira, N. C. D., Suzukawa, A. K., Pereira, C. B., Santos, H. V., Hanel, A., Albuquerque, F. A. D., & Scapim, C. A. (2018). Popcorn genotypes resistance to fallarmy worm. Ciência Rural, 48(2), 1-6. DOI: https://doi.org/10.1590/0103-8478cr20170378
https://doi.org/https://doi.org/10.1590/... ) and Sanches et al. (2019Sanches, R. E., Suzukawa, A. K., Contreras-Soto, R. I., Rizzardi, D. A., Kuki, M. C., Zeffa, D. M., ... Scapim, C. A. (2019). Multivariate analysis reveals key traits of fall armyworm resistance in tropical popcorn genotypes. Bragantia, 78(2), 175-182. DOI: https://doi.org/10.1590/1678-4499.20180230
https://doi.org/https://doi.org/10.1590/... ) and as in this study, the application of multivariate analysis allowed a significant reduction in the number of traits. This raises the chances of a more effective selection, since the smaller number of traits prevents the effect of interrelationships among them, avoiding redundancy and mistakes in the process of selecting promising genotypes (Aaliya et al., 2016Aaliya, K., Qamar, Z., Ahmad, N. I., Ali, Q., Munim, F. A., & Husnain, T. (2016). Transformation, evaluation of gtgene and multivariate genetic analysis for morpho-physiological and yield attributing traits in Zea mays. Genetika, 48(1), 423-433. DOI: http://doi.org/10.2298/GENSR1601423A
https://doi.org/http://doi.org/10.2298/G... ).

The selected traits can be described as directly related to the damage caused by S. frugiperda larvae and can be considered key traits for the resistance of popcorn genotypes, in that resistant genotypes will have a smaller leaf area consumed (LAC) and a lower food intake weight (IW) and relative consumption rate (RCR), as well as less apparent digestibility (AD), which will result in a lower final larva weight (FW) (Souza Crubelati-Mulat et al., 2020Souza Crubelati-Mulat, N. C., Baleroni, A. G., Contreras-Soto, R. I., Ferreira, C. J. B., Castro, C. R., Albuquerque, F. A., & Scapim, C. A. (2020). Evaluation of resistance to Spodoptera frugiperda in sweet and field corn genotypes. Maydica, 64(3), 1-7.).

In this study, the use of multivariate techniques for selection and grouping was used as a way to validate both the selection of resistant genotypes and the influence of the chosen traits on the identification of resistance.

In multivariate procedures, one of the most widely used techniques is hierarchical clustering (Oliveira et al., 2018Oliveira, N. C. D., Suzukawa, A. K., Pereira, C. B., Santos, H. V., Hanel, A., Albuquerque, F. A. D., & Scapim, C. A. (2018). Popcorn genotypes resistance to fallarmy worm. Ciência Rural, 48(2), 1-6. DOI: https://doi.org/10.1590/0103-8478cr20170378
https://doi.org/https://doi.org/10.1590/... ; Paiva et al., 2016Paiva, L. A., Corrêa, F., Silva, C. L. T., de Lima Moura, T., Silva, F. C., Silva Araújo, M., & Jesus, F. G. (2016). Resistance of corn genotypes to fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae). African Journal of Biotechnology, 15(35), 1877-1882. DOI: https://doi.org/10.5897/AJB2016.15399
https://doi.org/https://doi.org/10.5897/... ; Sanches et al., 2019Sanches, R. E., Suzukawa, A. K., Contreras-Soto, R. I., Rizzardi, D. A., Kuki, M. C., Zeffa, D. M., ... Scapim, C. A. (2019). Multivariate analysis reveals key traits of fall armyworm resistance in tropical popcorn genotypes. Bragantia, 78(2), 175-182. DOI: https://doi.org/10.1590/1678-4499.20180230
https://doi.org/https://doi.org/10.1590/... ). By these methods, genotypes are grouped by a process that is repeated at several levels, establishing a dendrogram, with no predetermined optimal number of groups. For this case Cruz, Carneiro, and Regazzi (2014Cruz, C. D., Carneiro, P. C. S., & Regazzi, A. J. (2014). Modelos biométricos aplicados ao melhoramento genético. Viçosa, MG: Editora UFV.) described different forms of representing the clustering structure based on the distance between the genotype pairs, for which UPGMA is the most commonly used method (Sanches et al., 2019; Sathua, Shahi, Mahato, Gayatonde, & Kumar, 2018Sathua, S. K., Shahi, J. P., Mahato, A., Gayatonde, V., & Kumar, P. (2018). Molecular diversity analysis of maize (Zea mays L.) inbreds using SSR markers. Electronic Journal of Plant Breeding, 9(3), 1122-1129. DOI: https://doi.org/10.5958/0975-928X.2018.00140.0
https://doi.org/https://doi.org/10.5958/... ; Vivodík, Petrovičová, Balážová, & Gálová, 2017Vivodík, M., Petrovičová, L., Balážová, Ž., & Gálová, Z. (2017). Genetic diversity of maize accessions (Zea mays L.) Cultivated from Europe using microsatellites markers. In Agrobiodiversity for improving nutrition, health and life quality (p. 524-528). Retrieved on May 5, 2022 from Retrieved on May 5, 2022 from https://agrobiodiversity.uniag.sk/scientificpapers/article/view/134
https://agrobiodiversity.uniag.sk/scient... ).

The relationships between the varieties UEM J1, Zapalote Chico, Composto Márcia, Arachida, and Composto Gaúcho can be observed in a graph of the results of the dendrogram based on Mahalanobis’ generalized distance, grouped by the UPGMA method (Figure 1). The high cophenetic correlation coefficient (CCC = 0.93) indicated an optimal fit between the original data and the dissimilarity matrix and high accuracy of the experimental results. The dendrogram shows two main groups, making the discrimination of the resistance and/or susceptibility of each genotype group possible (Friendly & Sigal, 2017Friendly, M., & Sigal, M. (2017). Graphical methods for multivariate linear models in psychological research: An R tutorial. The Quantitative Methods for Psychology, 13(1), 20-45. DOI: https://doi.org/10.20982/tqmp.13.1.p020
https://doi.org/https://doi.org/10.20982... ; Juárez et al., 2012Juárez, M. L., Murúa, M. G., García, M. G., Ontivero, M., Vera, M. T., Vilardi, J. C., ... Willink, E. (2012). Host association of Spodoptera frugiperda (Lepidoptera: Noctuidae) corn and rice strains in Argentina, Brazil, and Paraguay. Journal of Economic Entomology, 105(2), 573-582. DOI: https://doi.org/10.1603/EC11184
https://doi.org/https://doi.org/10.1603/... ; Omoto et al., 2016Omoto, C., Bernardi, O., Salmeron, E., Sorgatto, R. J., Dourado, P. M., Crivellari, A., ... Head, G. P. (2016). Field‐evolved resistance to Cry1Ab maize by Spodoptera frugiperda in Brazil. Pest Management Science, 72(9), 1727-1736. DOI: https://doi.org/10.1002/ps.4201
https://doi.org/https://doi.org/10.1002/... ).

The first group contained the variety Zapalote Chico, considered by several authors as resistant to S. frugiperda (Oliveira et al., 2018Oliveira, N. C. D., Suzukawa, A. K., Pereira, C. B., Santos, H. V., Hanel, A., Albuquerque, F. A. D., & Scapim, C. A. (2018). Popcorn genotypes resistance to fallarmy worm. Ciência Rural, 48(2), 1-6. DOI: https://doi.org/10.1590/0103-8478cr20170378
https://doi.org/https://doi.org/10.1590/... ; Sanches et al., 2019Sanches, R. E., Suzukawa, A. K., Contreras-Soto, R. I., Rizzardi, D. A., Kuki, M. C., Zeffa, D. M., ... Scapim, C. A. (2019). Multivariate analysis reveals key traits of fall armyworm resistance in tropical popcorn genotypes. Bragantia, 78(2), 175-182. DOI: https://doi.org/10.1590/1678-4499.20180230
https://doi.org/https://doi.org/10.1590/... ), and variety Arachida. The varieties Composto Gaúcho, UEM J1, and Composto Márcia were grouped in the second (Figure 1).

The proximity of the varieties Arachida and Zapalote Chico may suggest resistance of variety Arachida to S. frugiperda, since the grouping was based on the same traits for both varieties. The results show (Table 3) that the varieties Arachida and Zapalote Chico had the lowest values for index LAC, IW, and RCR, aside from the low AD and FW, which indicate antixenosis resistance of these two varieties. The cultivar Zapalote Chico is so far the only one treated as resistant to S. frugiperda, so the identification of new resistant cultivars can be considered a great gain to the breeding since it allows to explore differents sources of resistance.

Moreover, the performance of variety Composto Gaúcho was the worst between the evaluated cultivars, and it was characterized as the most susceptible to S. frugiperda, based on the evaluated traits, since the values of FW, IW, RCR, AD, and LAC were high (Table 3).

The dendrogram was confirmed and analyzed in more detail by canonical variable analysis, grouped by Tocher’s clustering (Figure 2), which showed the presence of three distinct groups. Tocher's optimization is a clustering method based on the formation of groups whose distances within are shorter than the distances between groups. At the end of the process, the number of groups and accessions contained in each group are computed. This method was applied as suggested by Cruz et al. (2014Cruz, C. D., Carneiro, P. C. S., & Regazzi, A. J. (2014). Modelos biométricos aplicados ao melhoramento genético. Viçosa, MG: Editora UFV.) and is an important way of determining different groups based on different traits, together with the techniques of dissimilarity analysis and analysis of canonical variables (Paiva et al., 2016Paiva, L. A., Corrêa, F., Silva, C. L. T., de Lima Moura, T., Silva, F. C., Silva Araújo, M., & Jesus, F. G. (2016). Resistance of corn genotypes to fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae). African Journal of Biotechnology, 15(35), 1877-1882. DOI: https://doi.org/10.5897/AJB2016.15399
https://doi.org/https://doi.org/10.5897/... ).

The analysis of canonical variables explained 95.12% of the total variation between the five traits. When analyzing the dispersion of the scores of the first two canonical variables, there was an agreement with the previous groupings, confirming the results, as well as the choice of the traits used in the study of S. frugiperda resistance. The first canonical variable (VC 1) explained 87.85% of the total variation and the second (VC 2) 25.6% (Figure 2).

Tocher's grouping method, based on the analysis of canonical variables, grouped the varieties Zapalote Chico and Arachida again. Similarly to the variety Zapalote Chico, considered resistant by several authors because of its known resistance classified as antixenosis, i.e., feeding avoidance (Oliveira et al., 2018Oliveira, N. C. D., Suzukawa, A. K., Pereira, C. B., Santos, H. V., Hanel, A., Albuquerque, F. A. D., & Scapim, C. A. (2018). Popcorn genotypes resistance to fallarmy worm. Ciência Rural, 48(2), 1-6. DOI: https://doi.org/10.1590/0103-8478cr20170378
https://doi.org/https://doi.org/10.1590/... ; Sanches et al., 2019Sanches, R. E., Suzukawa, A. K., Contreras-Soto, R. I., Rizzardi, D. A., Kuki, M. C., Zeffa, D. M., ... Scapim, C. A. (2019). Multivariate analysis reveals key traits of fall armyworm resistance in tropical popcorn genotypes. Bragantia, 78(2), 175-182. DOI: https://doi.org/10.1590/1678-4499.20180230
https://doi.org/https://doi.org/10.1590/... ) also in this study, some characteristics of variety Arachida like final larva weight, food intake weight, relative consumption rate, apparent digestibility and leaf area consumed also indicate antixenosis resistance. The varieties UEM J1 and Composto Márcia, which, according to the means of the analyzed variables, were moderately resitant, were contained in the second group. Variety Composto Gaúcho remained isolated from the other evaluated varieties, as the most susceptible to S. frugiperda, based on its performance with regard to the analyzed traits (Table 3).

Factor analysis is a multivariate statistical method that has been applied in agronomic studies for a relatively short time (Keith & Reynolds, 2018Keith, T. Z., & Reynolds, M. R. (2018). Using confirmatory factor analysis to aid in understanding the constructs measured by intelligence tests. In D. P. Flanagan, & E. M. McDonough (Eds.), Contemporary intellectual assessment: Theories, tests, and issues (p. 853-900). New York, NY: Guilford Press). This analysis explains the relationships observed between traits and removes possible redundancies or duplications from a set of correlated phenotypic data (Cattell, 1965Cattell, R. B. (1965). A biometrics invited paper. Factor analysis: An introduction to essentials I. The purpose and underlying models. Biometrics, 21(1), 190-215. DOI: https://doi.org/10.2307/2528364
https://doi.org/https://doi.org/10.2307/... ). The method allows the selection of relevant traits, exploring the relationships and their variations, apart from generating important information about factors and genotypes (Vile et al., 2012Vile, D., Pervent, M., Belluau, M., Vasseur, F., Bresson, J., Muller, B., ... Simonneau, T. (2012). Arabidopsis growth under prolonged high temperature and water deficit: independent or interactive effects?. Plant, Cell & Environment, 35(4), 702-718. DOI: https://doi.org/10.1111/j.1365-3040.2011.02445.x
https://doi.org/https://doi.org/10.1111/... ). In this study, the factor analysis was directed to the identification of key traits related to resistance to S. frugiperda.

By factor analysis, the traits are replaced by a smaller number of latent traits, called factors. These factors group the traits, so that there is little or no variance within groups but maximum variation between groups (Cruz et al., 2014Cruz, C. D., Carneiro, P. C. S., & Regazzi, A. J. (2014). Modelos biométricos aplicados ao melhoramento genético. Viçosa, MG: Editora UFV.). By this technique, associated with the analysis of canonical variables and canonical correspondence, the traits that best discriminate genotypes for a given objective, here S. frugiperda resistance, can be efficiently selected.

In this study, the high commonality (from 0.7012 - A to 0.9912 - RCR), indicated traits with a high relation to the determination of antixenosis resistance, confirming the thesis that the selected traits can be considered key traits (Table 4). In the factor analysis, the first factor was determinant for traits IW, FW, and RCR and the second for AD and LAC (Table 4).

In this study, estimates above 0.90 were considered for one factor and low estimates for the other factor (Figure 3), which shows a high representativeness of the factor for the respective traits (Keith & Reynolds, 2018Keith, T. Z., & Reynolds, M. R. (2018). Using confirmatory factor analysis to aid in understanding the constructs measured by intelligence tests. In D. P. Flanagan, & E. M. McDonough (Eds.), Contemporary intellectual assessment: Theories, tests, and issues (p. 853-900). New York, NY: Guilford Press). The estimates of the other traits were intermediate between both determined factors, which despite having a certain degree of contribution to resistance, were not as expressive as the described traits. In the study of Sanches et al. (2019Sanches, R. E., Suzukawa, A. K., Contreras-Soto, R. I., Rizzardi, D. A., Kuki, M. C., Zeffa, D. M., ... Scapim, C. A. (2019). Multivariate analysis reveals key traits of fall armyworm resistance in tropical popcorn genotypes. Bragantia, 78(2), 175-182. DOI: https://doi.org/10.1590/1678-4499.20180230
https://doi.org/https://doi.org/10.1590/... ), a large number of traits was also reduced, and subsequently divided into two factors by factor analysis.

According to the factor analysis, the chosen traits are directly correlated with the selection of S. frugiperda resistant genotypes. Thus, we proceeded to the analysis of canonical correspondence.

Canonical correspondence analysis is an exploratory technique to simplify the structure of multivariate data variability, in which the traits are arranged in contingency tables, taking correspondence measures between rows and columns of the data matrix into account. According to Ter Braak (1987Ter Braak, C. J. F. (1987). The analysis of vegetation-environment relationships by canonical correspondence analysis. Vegetatio, 69(1-3), 69-77. DOI: https://doi.org/10.1007/978-94-009-4061-1_7
https://doi.org/https://doi.org/10.1007/... ), correspondence analysis is a method to determine an association system between the elements of two or more sets, to explain the association structure of the factors in question. Thus, graphs were constructed with the principal components of the rows and columns, allowing the visualization of the relationship between the sets, where the proximity of the points referring to the row and the column indicates an association and distance indicates repulsion. According to Nyfjäll (2002Nyfjäll, M. (2002). Aspects on correspondence analysis plots under complex survey sampling designs. Uppsala University, Department of Information Science. Retrieved on Aug. 17, 2022 from 17, 2022 from https://urn:nbn:se:uu:diva-42867
https://urn:nbn:se:uu:diva-42867... ), one of the great advantages of canonical correspondence analysis is that relationships can be detected by this technique that would not have been perceived if the analysis were based on trait pairs. In addition, it is highly flexible in the data traits, since no theoretical model of probability distribution must be adopted. A rectangular matrix containing non-negative data is sufficient, which in the field of breeding, makes it possible to masterfully relate the effects of different traits on specific genotypes.

Canonical correspondence analysis explained 99.67% of the total variation between the genotypes and the respective traits evaluated. The first canonical correspondence axis (CCA 1) explained 99.20% of the total variation and the second axis (CCA 2) accounted for 0.47% (Figure 4). Most of the total variation was already explained in the first CCA, which is desirable, for increasing the accuracy between the cluster and the estimated scores (Ter Braak, 1987Ter Braak, C. J. F. (1987). The analysis of vegetation-environment relationships by canonical correspondence analysis. Vegetatio, 69(1-3), 69-77. DOI: https://doi.org/10.1007/978-94-009-4061-1_7
https://doi.org/https://doi.org/10.1007/... ; Nyfjäll, 2002Nyfjäll, M. (2002). Aspects on correspondence analysis plots under complex survey sampling designs. Uppsala University, Department of Information Science. Retrieved on Aug. 17, 2022 from 17, 2022 from https://urn:nbn:se:uu:diva-42867
https://urn:nbn:se:uu:diva-42867... ).

The biplot of canonical correspondence analysis shows that the traits IW, FW, and AD were determinant for the resistance expressed by the varieties Zapalote Chico and Arachida, in which variable AD had a greater influence, due to its spatial angular proximity in the graph with the said varieties. For the varieties Composto Gaúcho and Composto Márcia, the traits that most contributed to the determination of resistance or susceptibility were RCR and LAC. For variety UEM J1, the influence of variable LAC was the highest.

Our results showed that the selected traits were efficient in discriminating the varieties regarding resistance and susceptibility to S. frugiperda by the applied analyses. The identification of key traits in the description of resistant genotypes can guide future studies, allow greater emphasis on specific traits and consequently a more effective selection regarding resistance in popcorn genotypes.

Conclusion

Variety Arachida was identified resistant to S. frugiperda by antixenosis and can be used as a source of favorable alleles for the future development of resistant popcorn hybrids. The traits relative consumption rate (RCR), apparent digestibility (AD) and leaf area consumed (LAC) were efficient and considered key traits for the identification of S. frugiperda resistant genotypes by antixenosis.

Publication Dates

History