ABSTRACT
Babassu palm (Attalea speciosa) is a native species from the Cerrado-Amazon transition zone with socio-economic importance for agro-extractivist communities in Northeastern Brazil that rely on products derived from its fruits. Babassu productivity is declining due to unsustainable practices and habitat loss. Hence, this research aimed to estimate genetic parameters related to fruit traits of babassu and indicate superior populations for coconut mesocarp production. Babassu trees and fruits were sampled from natural populations in five cities of Maranhão State - Brazil, totaling 150 individuals and 1,500 fruits. The number of bunches and fruits per plant, plant and stipe height, and various fruit and seed morphological traits were measured and used for estimating genetic parameters by the REML/BLUP method with Selegen Software. Fruit width, fruit weight, and mesocarp weight showed high broad-sense heritability and accuracy, suggesting a strong genetic influence and making them excellent candidates for genetic improvement. The high coefficients of relative variation are suitable for selecting promising babassu genotypes. Babassu fruits sampled on Coroatá and Buritirana had the best characteristics for fruit and mesocarp production. We found significant genotypic correlations between various fruit traits (size and weight) that can simplify the selection of high-yield genotypes. Our results are the first report on genetic parameters of important babassu fruit traits, providing valuable data for future breeding programs aiming to increase productivity.
Index terms:
Fruit biometric traits; plant breeding; plant genetic resources
RESUMO
A palmeira babaçu (A. speciosa), nativa da zona de transição Cerrado-Amazônia, é de vital importância para comunidades agroextrativistas do Nordeste brasileiro, que utilizam subprodutos de seus frutos. Contudo, sua produtividade está em declínio devido a práticas de manejo insustentáveis e perda de habitat. Assim, esta pesquisa teve objetivou estimar parâmetros genéticos relacionados a características de interesse de frutos do babaçu, bem como indicar populações superiores para produção de mesocarpo de coco, acessando populações naturais do Estado do Maranhão, Brasil. Foram amostradas árvores em três localidades de cinco cidades do Maranhão, totalizando 150 indivíduos e 1500 frutos. Foram mensurados o número de cachos e frutos por planta, altura da planta e do estipe, e características morfológicas dos frutos e sementes. Os parâmetros genéticos foram estimados pela metodologia REML/BLUP, pelo Software Selegen. Os caracteres largura do fruto, peso do fruto e peso do mesocarpo apresentaram elevada herdabilidade e acurácia, sugerindo forte influência genética e, que esses caracteres têm potencial para o melhoramento genético da espécie. Os elevados coeficientes de variação relativa observados indicam variação genética suficiente para seleção de genótipos promissores. Frutos amostrados em Coroatá e Buritirana apresentaram melhor desempenho para características para produção. Foram encontradas correlações genotípicas significativas entre vários caracteres de tamanho e peso dos frutos, o que pode simplificar a seleção de genótipos para aumento da produtividade. Este é o primeiro estudo reportando parâmetros genéticos para caracteres de frutos em babaçu, fornecendo informações para seleção de genótipos superiores para programas de melhoramento da espécie.
Termos para indexação:
Características biométricas de frutos; melhoramento de plantas; recursos genéticos vegetais
Introduction
The Cerrado biome, primarily located in Brazil, covers over 2 million km² and is a global biodiversity hotspot (Colli, Vieira, & Dianese, 2020). It hosts more than 12,000 plant species, many are rare and endemic, and has been threatened by agricultural expansion and habitat loss (Mittermeier et al., 2004; Mendonça et al., 2008). Ballesteros-Mejia, Lima and Collevatti (2020) evidenced the highest plant genetic diversity and richness of Cerrado in the northern and central regions.
Babassu is a common name for several palm tree species of the Orbignya and Attalea genera. The species Attalea speciosa (Mart. ex Spreng.) is native to the transition zone between the Cerrado and the Amazon open forests (Fakhouri, Da Silva, & Velasco, 2021). This palm can reach up to 30 m in height and 8 m in length and starts producing fruits around 8 to 10 years (Santos et al., 2022). Babassu fruits consist of four main parts: the epicarp in the outer layer, the mesocarp, the rigid endocarp, and the almonds (3 to 4 per fruit) (Fakhouri, Da Silva, & Velasco, 2021).
Besides the classical utilization in charcoal production, handicrafts, almonds-oil extraction, and mesocarp flour, the scientific literature has evidenced enormous medicinal potential of active compounds of babassu nuts and mesocarp (Silva et al., 2018; Sales et al., 2020; Barroqueiro et al., 2024; Pinheiro, Corrêa, & Araujo 2024), as well as a biofuel production from nuts residuals (Protásio et al., 2014). However, the babassu production chain has a vulnerable socioeconomic background, since many families from Maranhão, Piauí, Tocantins, and Pará states depend on coconut extraction for their subsistence, living in poverty and social exclusion (Lemos, 2014). Women (known as ‘coconut crackers’) play a key role in these communities and family subsistence (Porro, Veiga, & Mota, 2011).
Babassu annual production has gradually reduced in the last years, and the State of Maranhão is responsible for 89% of Brazilian production (Instituto Brasileiro de Geografia e Estatística - IBGE, 2024). The species is facing habitat and genetic losses, caused by cattle grazing advances, rudimentary agricultural practices, and unsustainable land management (Santos et al., 2015; Lemos, 2020; Santos et al., 2022). All these factors highlight the importance of accessing babassu natural genetic variability, not only to reinforce conservation and breeding programs, but also to support agro-extractivist communities obtaining palm trees with better productive traits. Therefore, this research aimed to estimate genetic parameters related to interesting fruit traits of babassu (A. speciosa), as well as indicate superior populations of babassu for coconut mesocarp production, by accessing natural populations of the Maranhão State.
Material and Methods
Sampling area and plant material
The study was conducted in five municipalities of Maranhão State, Brazil (São Luís, Cantanhede, Coroatá, Viana, and Buritirana) (Figure 1).
Location of the five municipalities where Attalea speciosa populations were sampled in Maranhão State, Brazil. Red dots: Buritirana, located on the Western Mesoregion of the state (sampled areas were three private properties: Raimundo Mesquita, Dedeu, and Lucinha). Orange dots: Cantanhede, located on the Northern Mesoregion of the state (sampled areas were Jacaré and Candiba communities). Yellow dots: Coroatá, located on the Eastern Mesoregion (sampled areas were Bem Fica Village and Umaitá Village). Green dots: São Luis, located on the North (sampled areas were Tajipuru and Coquilho). Blue dots: Viana, located on the North Region (samples areas were Quilombo Cocal, Nova Vila Village, and Baias Nós Vamos).
In each municipality, three distinct areas of A. speciosa populations were analyzed, with a minimum distance of 500 m between areas. Ten individuals at the reproductive stage were randomly selected in each area, totaling 150 georeferenced individuals, and ten fruits per plant (a total of 1500 fruits) were sampled in October and November 2021, a transitional period from the dry to the rainy season. Fruits were packed in nylon bags and transported to the Associação das Mulheres Quebradeiras de Coco of Povoado Candiba (Cantanhede, Maranhão, Brazil) for extraction of the mesocarp and seeds. Mesocarp was weighed, packed in plastic bags after extraction, and stored at 10 °C at the Post-Harvest Laboratory (LAPOC), Maranhão State University (São Luís, MA, Brazil).
Morphometric descriptors of trees and fruits
The parameters analyzed in trees were the number of bunches per plant (NBP, unit), number of fruits per bunch (NFB, unit), stipe circumference at ground level (SCG, cm), stipe circumference at breast level (SCB, cm), and plant height (PH, m). A 20-m measuring tape was used to determine SCB. PH was determined using an Opti-Logic 1000LH rangefinder, which also allowed to measure the average distance between plants of 5 m.
For the morphological characterization of the fruits, a digital caliper (Zaas Precision®, Amatools, Piracicaba-SP, Brazil) was used. Analyzed traits were: fruit length (FL, mm), fruit width (FW, mm), average fruit weight (AFW, g), nut length (NL, mm), average nut weight (ANW, g), and average mesocarp weight (AMW, g).
Estimation of genetic parameters and covariance components
Statistical analyses were performed using the software SELEGEN-REML/BLUP (Statistical System and Computerized Genetic Selection by Mixed Linear Models) version 2016, model 83 (Resende et al., 2016). The analysis of deviance (ANADEV) or likelihood ratio test (LRT) were carried out using the model with and without the effect, according to Resende (2007). Results were compared with the chi-square (χ²) value (at 10, 5, and 1% significance level), considering one degree of freedom.
The genetic parameters were estimated using the mixed model Restricted Maximum Likelihood procedure and Best Linear Unbiased Prediction (REML/BLUP) methodology, following the model below (Equation 1):
where: y is the observed data vector, u is the fixed effects vector (overall mean), g is the random effect vector of total genotypic effects, and e corresponds to the vector of errors/residuals (random). X and Z represent the incidence matrices for vectors u and g, effects, respectively.
The estimates of the variance components by REML via the EM algorithm were given by the following mixed model equations (Equations 2, 3, and 4):
C22 is from the form
Where: r(X) is the rank or number of linearly independent columns of X, N g is the number of random elements (individuals), A is the additive genetic relatedness matrix, tr is the matrix of dash operator, given by the sum of the diagonal elements of the matrix; N is the total number of data, and h 2 is the individual broad sense heritability, also calculated by .
Estimated parameters were the broad-sense heritability , selective accuracy (raâ), coefficient of genetic variation (CVg%), experimental variation (CVexp%), and coefficient of relative variation (CVr= CVg%/CVexp%), where the magnitude of and raâ were classified according to Resende et al. (1998) (Table 1).
Classification of the magnitudes of heritability and selective accuracy (raâ) used to estimate the genetic parameters of A. speciosa.
Estimation of selection index
The predicted genotypic values were used to calculate the selection index, based on the sum of ranks to classify the genotypes in terms of each trait (Mulamba & Mock, 1978). Subsequently, the values of each trait for the accessions were summed, resulting in an overall value considered as the selection index (Cruz et al., 2014).
Results and Discussion
CVr contributes to the detection of genetic variability in a population. The CVr for FL, AFW, and AMW were higher than unity (1.47, 1.34, and 1.15, respectively), suggesting higher genetic variation for the selection of promising babassu genotypes based on these traits (Table 2).
Estimates of genetic parameters from morphological traits of A. speciosa. Legend of traits: number of bunches per plant (NBP), number of fruits per bunch (NFB), stipe circumference at breast level (SCB), plant height (PH), fruit length (FL), fruit width (FW), average fruit weight (AFW), nut length (NL), and average mesocarp weight (AMW).
We found high values of broad-sense heritability (>0.5) and high selective accuracy (80%) for FW, AFW, and AMW , confirming their precision and confidence in the estimated genetic values (Table 2). FW, ANW, and AMW are important for productivity, and their elevated heritability and accuracy indicate a substantial genetic influence, making these traits excellent candidates for genetic improvement.
NBP, NBP, PH, FL, and NL presented moderate significant broad-sense heritabilities (among 0.15 and 0.5) and high accuracies (>80%) (Table 2), indicating a good potential for selection. In a study with peach palm (Bactris gasipaes) heritability was reported as >0.2 for the total bunch weight (Farias-Neto et al., 2013).
Selective accuracy represents a correlation between actual and predicted genetic values, and the higher its value, the higher the confidence in the evaluation of individuals (Pivetta et al., 2020). Neves et al. (2023) studied populations of licuri palm (Syagrus coronata), also native to the Cerrado biome, and found similar CVr (obtained from the evaluations of three fruits per bunch) with an accuracy rate of 90%. In a study with 51 progenies of Acrocomia aculeata (coco-de-espinho) a similar CVr was found, with 90% accuracy (Manfio et al., 2011). In our data, accuracy was also high for all traits evaluated, reflecting the great reliability of the estimated genetic values.
SCB was the only trait that showed considered low, but significant, with a high accuracy. Some babassu stipes were affected by fires in the period of data collection, which may have influenced the SCB results. In Maranhão State, the post-burned areas have a massive presence of babassu and other secondary species that are resistant to fire. Babassu palm trees face pressure from large cattle ranchers, but also from family farmers who often remove babassu trees to establish their cultivation gardens, using felling and fire as instruments to clear their cultivation areas (Freitas, 2000; Lemos, 2020).
The analysis of deviance evidenced a significant effect of genotype (p<0.05) in all evaluated traits (Table 3). The high genetic variation reveals the possibility of obtaining genetic gains through the selection criteria.
The highest positive values of genotypic correlations were found between the factors FW and AFW (0.99), and AFW and AMW (0.95). FW and AMW, FL and AMW, FL and AFW, and FL and FW had the same value of 0.91 (Table 4).
Genotypic correlations of A. speciosa traits indicate that fruit size directly influences FW and, consequently, AMW, facilitating the selection process of individuals with interesting agro-extractivist production characteristics. This fact only confirms what the extractivists already observed in the field, as they prefer larger fruits that generate higher yields for seeds and mesocarp. In a study with açaí fruits (Euterpe oleracea) a positive genetic correlation was found between the total weight of the fruits and the number of bunches (Farias Neto et al., 2008). Marçal et al. (2015) also found a positive correlation between the weight of the juçara fruit and the longitudinal diameter of fruits (Euterpe edulis).
The visual variability in fruit dimensions sampled from the different municipalities is shown in Figure 2. The average rank index for the best combinations of the studied traits is separated by the municipalities (Table 5). According to Mulamba and Mock (1978), the lowest values on the average rank indicate a more favorable position of the evaluated traits and the highest values represent an unfavorable combination. Fruits sampled in Coroatá and Buritirana ranked highest, demonstrating that these areas have plants with more favorable combinations of morphological traits (such as fruit weight, fruit size, etc) that are interesting for better production.
Visual variability of A. speciosa fruits collected in different municipalities of the State of Maranhão. (A) Buritirana - medium to large sized fruits, with an oval-ellipsoid shape. (B) Cantanhede - smaller fruits, with a more elongated shape. (C) Coroatá - longer fruits, with a wider base. (D) São Luís - ellipsoid-shaped and uniform fruits. (E) Viana - medium size, rounded-ellipsoid shape. Scale bar: 1 cm.
Fruits from Viana and Cantanhede had the lowest positions in the ranking, and these sampled areas are mostly located in natural reserves and/or protected areas (Figure 1). Studies with palm trees using the evaluation by the average rank index showed considerable genetic variability, favoring breeding programs (Farias Neto et al., 2013; Simiqueli et al., 2018).
The sampled areas were spatially distributed in different microregions of the State (Figure 1) and did not exhibit uniform vegetation cover. Our objective was to capture the greatest possible variability within the state. Some of the sampled areas were greatly transformed by human activity (cassava cultivation, extensive pastures, degraded areas, and secondary forests), but babassu trees also predominate. Gehring et al. (2011) found that babassu becomes dominant in areas affected by constant fires and degraded land, as it is quite tolerant to slash-and-burn practices in traditional agriculture.
Our results are innovative because we identified traits with high and significant heritabilities, important to establishing selection criteria for future babassu breeding programs. This study is the first report on the genetic control of productive traits of babassu and offers the opportunity to select individuals with superior fruit and mesocarp production based not only on their location but also on their morphological traits. However, further studies are necessary to accelerate the process of domestication of this palm.
Conclusions
The babassu populations analyzed exhibited genetic variability and the highest heritability and accuracy in fruit width, fruit weight, and mesocarp weight, indicating strong genetic control on these traits. We identified fruit weight and mesocarp weight as suitable traits for conducting genotypic correlation analyses. Populations of babassu from Coroatá and Buritirana municipalities produced fruits with superior characteristics. These data can assist efforts in selecting superior genotypes and developing technologies for the rational and sustainable exploitation of babassu palm fruits in the future.
Acknowledgments
We acknowledge the infrastructure and technical support of the Post-Harvest Laboratory of Maranhão State University - LAPOC.
References
- Ballesteros-Mejia, L., Lima, J. S., & Collevatti, R.G. (2020). Spatially-explicit analyses reveal the distribution of genetic diversity and plant conservation status in Cerrado biome. Biodiversity and Conservation, 29(5):1537-1554.
- Barroqueiro, Â. T. L. S. et al. (2024). The anti-infective and immunologic effect of babassu (Attalea speciosa, Mart. ex Spreng) reduces mortality induced by MRSA-Staphylococcus aureus Journal of Ethnopharmacology, 320:117363.
- Colli, G. R., Vieira, C. R., & Dianese, J. C. (2020). Biodiversity and conservation of the Cerrado: Recent advances and old challenges. Biodiversity and Conservation, 29(5):1465-1475.
- Cruz, C. D., Carneiro, P. C. S., & Regazzi, A. J. (2014). Modelos biométricos aplicados ao melhoramento genético Viçosa (MG): UFV. 668p.
- Fakhouri, F. M., Da Silva, L. R., & Velasco, J. I. (2021). Attalea speciosa (Orbignya phalerata). In F. Freitas de Lima., L. H. Lescano., & I. Pires de Oliveira. Fruits of the Brazilian Cerrado Cham: Springer. (pp.125-139).
- Farias Neto, J. T. D., Clement, C. R., & RESENDE, M. D. V. D. (2013). Estimativas de parâmetros genéticos e ganho de seleção para produção de frutos em progênies de polinização aberta de pupunheira no Estado do Pará, Brasil. Bragantia, 72(2):122-126.
- Farias Neto, J. T. D. et al. (2008). Estimativas de parâmetros genéticos e ganhos de seleção em progênies de polinização aberta de açaizeiro. Revista Brasileira de Fruticultura, 30(4):1051-1056.
- Freitas, A. C. R. (2000). A crise ecológica dos sistemas de cultivo “corte e queima” na micro-região da pré-Amazônia maranhense. Papers do NAEA, 1(1):1-22.
- Gehring, C. et al. (2011). Allometry of the babassu palm growing on a slash-and-burn agroecosystem of the eastern periphery of Amazonia. Acta Amazonica, 41(1):127-134.
-
Instituto Brasileiro de Geografia e Estatística - IBGE. (2024). Produção da extração vegetal e da silvicultura Available in: https://sidra.ibge.gov.br/tabela/289
» https://sidra.ibge.gov.br/tabela/289 - Lemos, J. J. S. (2014). Maranhão estado (ainda) rico em recursos naturais com população empobrecida. In J. L. A. Silveira. Ensaios sobre a economia maranhense São Luís (MA): EDUFMA, (pp. 171-230).
- Lemos, J. J. S. (2020). Nordeste e suas contradições. In J. J. S. Lemos. Vulnerabilidades induzidas no semiárido Fortaleza (CE): Imprensa Universitária (pp. 09-17).
- Manfio, C. E. et al. (2011). Repetibilidade em características biométricas do fruto de macaúba. Ciência Rural, 41(1):70-76.
- Marçal, T. D. S. et al. (2015). Correlações genéticas e análise de trilha para caracteres de fruto da palmeira juçara. Revista Brasileira de Fruticultura, 37(3):692-698.
- Mendonça, R. D. et al. (2008). Flora Vascular do Bioma Cerrado - Checklist com 12.356 espécies. In S. M. Sano., S. P. Almeida., & J. F. Ribeiro. Cerrado: Ecologia e flora Brasília, Brasil: Embrapa Cerrados (pp. 423-1279).
- Mittermeier, R. A. et al. (2004). Hotspots: Earth’s richest and most endangered terrestrial ecoregions Mexico, Mexico City: CEMEX, 391p.
- Mulamba, N. N., & Mock, J. J. (1978). Improvement of yield potential of the ETO blanco maize (Zea mays L.) population by breeding for plant traits. Egyptian Journal of Genetics and Cytology, 7(1):40-51.
- Neves, G. F. et al. (2023). Genetic diversity of Syagrus coronata and principal component analysis of phenotypic traits: A palm from the brazilian semiarid biome. Biodiversity and Conservation, 32(13):4275-4293.
- Pinheiro, E. M., Corrêa, T. R., & Araujo, J. R. G. (2024). A palmeira de babaçu e a importância do mesocarpo para alimentação humana: uma revisão. Contribuciones a las Ciencias Sociales, 17(7):e8292.
- Pivetta, R. C. et al. (2020). Tracking the occurrence of psychotropic pharmaceuticals in Brazilian wastewater treatment plants and surface water, with an assessment of environmental risks. Science of the Total Environment, 727:138661.
- Porro, N., Veiga, I., & Mota, D. (2011). Traditional communities in the Brazilian Amazon and the emergence of new political identities: the struggle of the quebradeiras de coco babaçu - babassu breaker women. Journal of Cultural Geography, 28(1):123-146.
- Protásio, T. P. et al. (2014). Babassu nut residues: Potential for bioenergy use in the North and Northeast of Brazil. Springer Plus, 3:124.
- Resende, M. D. (2007). Selegen-Reml/Blup: Sistema estatístico e seleção genética computadorizada via modelos lineares mistos Colombo, PR: Embrapa Florestas.
- Resende, M. D. V. et al. (1998). Efeito do tamanho amostral na estimativa da herdabilidade em espécies perenes. Floresta, 28:51-63.
- Resende, M. D. V. (2016). Software Selegen-REML/BLUP: A useful tool for plant breeding. Crop Breeding and Applied Biotechnology, 16(4):330-339.
- Sales, A. R. R. et al. (2020). Caracterização físico-química do oléo de coco babaçu industrial e artesanal e suas aplicações tecnológicas. Brazilian Journal of Development, 6(5):25734-25748.
- Santos, A. M. D. et al. (2022). What is the influence of anthropogenic impact on the population structure of Attalea speciosa Mart. ex Spreng. in the Brazilian Amazonian region? Acta Botanica Basilica, 36:e2020abb0543.
- Santos, M. F. et al. (2015). Genetic variation detected by RAPD markers in natural populations of babassu palm (Attalea speciosa Mart.). Genetics and Molecular Research, 14(2):6124-35.
- Silva, V. C. et al. (2018). dentification of phenolic compounds by LC/MS-MS and antioxidant and anti tyrosinase activities of the Attalea speciosa Mart. ex Spreng. Mesocarp. Journal of Chemical and Pharmaceutical Research, 9(1):268-276.
- Simiqueli, G. F. et al. (2018). Inbreeding depression as a cause of fruit abortion in structured populations of macaw palm (Acrocomia aculeata): Implications for breeding programs. Industrial Crops and Products, 112:652-659.
-
Editor de seção:
Renato Paiva




