Revista Brasileira de Engenharia Agrícola e Ambiental
Revista Brasileira de Engenharia Agrícola e Ambiental
DOI: http://dx.doi.org/10.1590/1807-1929/agriambi.v27n12p958-965
HIGHLIGHTS:
Salinity reduces growth and gas exchange of pomegranate.
Foliar-applied N dose of up to 1.31 g L-1 reduces the harmful effect of salinity on pomegranate growth.
Foliar nitrogen fertilization up to 1.59 g L-1 increases net photosynthesis.
ABSTRACT: Salinity is one of the abiotic stresses that affect gas exchange and growth of pomegranate. In this context,
the application of nitrogen fertilizer through the leaves can minimize these effects. Thus, the objective of present
study was to evaluate the effect of foliar nitrogen fertilization on pomegranate seedlings irrigated with brackish water.
The design used was randomized blocks, in an incomplete factorial scheme (Central Composite Design) with five
electrical conductivities of irrigation water - ECw (0.50, 1.15, 2.75, 4.35, and 5.00 dS m-1) and five doses of foliar
nitrogen fertilization - FNF (0, 0.33, 1.15, 1.97, and 2.30 g L-1), with four replicates and two plants per experimental
plot. Plant height, number of leaves, stem diameter, and gas exchange (stomatal conductance, net photosynthesis,
intercellular CO2 concentration, transpiration rate, instantaneous carboxylation efficiency, instantaneous water use
efficiency, and intrinsic water use efficiency) were evaluated at 60 days after the beginning of irrigation with saline
water. Salinity of irrigation water negatively affects the gas exchange of pomegranate seedlings. Foliar nitrogen
fertilization up to 1.31 g L-1 improves plant height and number of leaves in pomegranate seedlings under salt stress.
Foliar nitrogen fertilization up to 1.59 g L-1 increases the net photosynthesis of pomegranate seedlings.
Key words: Punica granatum L., photosynthesis, nitrogen, salinity
RESUMO: A salinidade é um dos estresses abióticos que afeta as trocas gasosas e o crescimento de romãzeira. Neste
sentido, à aplicação de adubo nitrogenado via foliar poderá minimizar estes efeitos. Com isso, o objetivo do presente
estudo foi avaliar o efeito da aplicação foliar com nitrogênio em mudas de romã irrigadas com águas salobra. O
delineamento utilizado foi de blocos casualizados, em esquema fatorial incompleto (Composto Central de Box),
sendo cinco condutividade elétrica da água - CEa (0,50; 1,15; 2,75; 4,35 e 5,00 dS m-1) e cinco doses de adubação
nitrogenada via foliar – ANF (0; 0,33; 1,15; 1,97 e 2,30 g L-1), com quatro repetições. A salinidade da água de irrigação
afeta negativamente as trocas gasosas de mudas de romãzeira. A adubação nitrogenada foliar até 1,31 g L-1 melhora o
crescimento em altura e no número de folhas em mudas de romãzeira sob estresse salino. Adubação nitrogenada
foliar até a dose de 1,59 g L-1 aumenta a fotossíntese das mudas de romãzeira.
Palavras-chave: Punica granatum L., fotossíntese, nitrogênio, salinidade
Figure 1. Data of maximum and minimum air temperature, precipitation, and average relative air humidity during the
experimental period in greenhouse
Pomegranate seeds cv. Mollar from fruits harvested in the Foliar nitrogen fertilization started at 25 DAE, with
orchard of the Universidade Federal Rural do Semi-Árido previous dissolution of the fertilizer in distilled water, followed
(UFERSA) were used. After acquiring the fruits, the seeds by application with an atomizer. In total, seven applications
were extracted and processed manually with the aid of a #¼ were performed, at interval of 10 days. The total volume applied
mesh sieve to remove the sarcotesta, to overcome the possible was 175 mL per plant, which provided 0, 58.16, 200, 341.84,
dormancy of the seeds, aiming at the standardization of and 400 mg of N per plant according to the increase of the
germination and the establishment of seedlings. Three seeds evaluated doses, based on study of Fátima et al. (2023).
were used per polyethylene bag, which was kept close to the The polyethylene bags had a capacity of 1.15 dm3, filled
field capacity from sowing to germination. Seed germination with substrate formed by 85% soil, 10% fine sand and 5% aged
began four days after sowing, extending to 25 days, when manure (Sakazaki et al., 2019). The physical and chemical
seedling emergence was established. At 25 days after sowing, characteristics of the substrate were evaluated according to
thinning was performed, leaving only one plant per bag, which methodologies proposed by EMBRAPA (2017) and Richards
was considered the most vigorous. (1954) (Table 2).
The electrical conductivities of irrigation water (ECw) The evaluations were carried out at 90 days after the
were obtained by dissolving sodium chloride (NaCl) in beginning of the treatment application. Gas exchange was
water from the supply system (0.5 dS m-1), until the required measured on the fourth leaf from the apex to the base, between
conductivities were obtained, the values being measured with a 9 and 10 a.m. with an infrared gas analyzer – IRGA (LI-6400XT,
microprocessor-based portable conductivity meter (model CD- LI-COR®, Nebraska, USA) with an air flow of 300 µmol s-1
860, Instrutherm®). The application of brackish water started 25 and humidity between 50-60%, 400 µmol mol-1 of CO2 and
days after emergence (DAE), with daily manual irrigation. The 1200 µmol m-2 s-1 coupled light source. Stomatal conductance
crop water demand was determined by the drainage lysimetry (gs – mol H2O m-2 s-1), net photosynthesis (A – μmol CO2 m-2
(Bernardo et al., 2019) method (Eq. 1), with the volume of s-1), intercellular CO2 concentration (Ci – mmol CO2 mol-1),
water lost through evapotranspiration the previous day being transpiration rate (E – mmol H2O m-2 s-1), instantaneous water
replaced daily, thus maintaining the soil at field capacity. use efficiency (WUE – [(μmol CO2 m-2 s-1) (mmol H2O m-2
s-1)-1]), intrinsic water use efficiency (iWUE - [(µmol CO2 m-2
( Va − Vd ) s-1) (mmol H2O m-2 s-1)-1]), and instantaneous carboxylation
VI = (1) efficiency (iCE – [(μmol CO2 m-2 s-1) (mmol CO2 mol-1)-1])
(1 − LF ) were evaluated.
To evaluate the growth of pomegranate seedlings, plant
where: height was analyzed, considering the distance from the collar
VI - volume of water to be used in the next irrigation (mL); to the apex of the plant, evaluated with a ruler graduated in
Va - volume of water applied in the previous irrigation cm; stem diameter was measured with a digital caliper and
event (mL); expressed in mm; and the number of leaves was determined
Vd - volume drained (mL); and, by counting fully formed leaves.
LF - leaching fraction (0.15). Data were subjected to normality (Shapiro-Wilk) and
homogeneity of variances (Bartlett) tests. Subsequently, an
Every 15 days a 15% leaching fraction was applied based analysis of variance (p ≤ 0.05) was performed, and a regression
on the volume applied in this period, in order to reduce the analysis was carried out in cases of significance. The statistical
accumulation of substrate salts. program R (R Core Team, 2021) was used.
Nitrogen doses were based on the need for 300 mg per
plant, proposed by Novais et al. (1991) for a 1 dm3 pot, the
Results and Discussion
highest evaluated dose being 400 mg per plant. The commercial
product Nitrotecnia-20 (Carbotecnia®), based on urea, with As observed in the summary of the analysis of variance
99 g L-1 of N, was used. (Table 3), a significant effect was found for the interaction
Table 2. Physical and chemical composition of the components of substrate used in the experiment
OM - Organic matter; SB - Sum of bases (Na+ + K+ + Ca2+ + Mg2+); CEC - Cation exchange capacity = SB + (H+ + Al3+); ECse - Electrical conductivity of the saturation extract;
SARse - Sodium adsorption ratio of the saturation extract = Na+ × [(Ca2+ + Mg2+)/2]1/2; ESP - Exchangeable sodium percentage (100 × Na+/CEC)
Table 3. Summary of the analysis of variance for stomatal conductance (gs), net photosynthesis (A), intercellular CO2
concentration (Ci), instantaneous carboxylation efficiency (iCE), intrinsic water use efficiency (iWUE), instantaneous water
use efficiency (WUE), plant height (PH), number of leaves (NL), and stem diameter (SD) of Punica granatum L. under foliar
nitrogen fertilization (DFN) and salinity of brackish irrigation water (ECw)
between brackish water and foliar nitrogen rates for water this efficiency is reduced. This behavior is a response to the
use efficiency - WUE, plant height - PH and number of leaves deleterious effects caused by salt stress, which induces this low
- NL. No significant effect was observed for the transpiration water consumption, limiting the absorption of toxic ions, as
variable - E. Individual effect was observed for the other observed by Nóbrega et al. (2022) in Mesosphaerum suaveolens.
analyzed variables (Table 3). The intrinsic water use efficiency (iWUE) increased up to
For stomatal conductance (gs) and net photosynthesis (A), the ECw of 2.20 dS m-1, reaching a value of 105.87 (Figure 2E).
there was a linear reduction as a function of the salinity of The greater iWUE observed in this study is associated with a
brackish irrigation water, with the highest values (0.076 mol plant mechanism of adaptation to stress conditions, leading to
m-2 s-1 and 7.142 μmol CO2 m-2 s-1, respectively) at the ECw of the maintenance of water in the tissues, as reported by Huang et
0.5 dS m-1, with decreases of 25.3 and 41.5%, at the highest al. (2015). The negative effects of salinity on gas exchange have
salinity of 5.00 dS m-1 (Figures 2A and 2B). The occurrence of been observed in plants of Solanum lycopersicum (Talebnejad &
these reductions is common in plants subjected to salt stress, Sepaskhah, 2016), Passiflora edulis (Silva et al., 2019), Psidium
closing their stomata as a defense mechanism of the plant. guajava L. (Bezerra et al., 2019), and Punica granatum L.
Consequently, CO2 absorption decreases, reducing the rate (Soares et al., 2021), results similar to those found in this study.
of net photosynthesis. This fact was observed by Soares et al. The foliar nitrogen fertilization promoted increases in
(2021) in pomegranate seedlings and by Silva et al. (2019) in gs, with the highest value obtained being 0.0747 mol m-2 s-1 at
yellow passion fruit, under saline conditions. the dose of 1.6 g L-1 of N, representing gains of 45.6% when
The intercellular CO2 concentration (Ci) was higher (277.70 compared to plants of the control treatment (dose 0), indicating
μmol CO2 m-2 s-1) in seedlings subjected to ECw of 0.5 dS m-1, the beneficial effect of foliar application of N on the gs of
followed by decreases up to ECw of 4.1 dS m-1 (232 μmol CO2 pomegranate plants (Figure 3A). This result may be related to
m-2 s-1), reaching decreases equivalent to 16.1% when compared the positive effect of nitrogen on plants, which is an important
to the values obtained in seedlings subjected to ECw of 0.5 dS component in the synthesis of photosynthetic pigments and
m-1 (Figure 2C). The reduction of Ci up to the salinity of 4.1 enzymes, as it improves stomatal regulation and the light
dS m-1 shows that the decrease in photosynthesis may be due saturation point, which explains the increases found for gs in
not only to stomatal factors, but also to factors of non-stomatal this study (Wang et al., 2016).
nature, such as biochemical alterations in the reduction of Foliar nitrogen fertilization stimulated A and Ci, with the
energy supply, which supposedly can affect the activity of highest values of 6.3457 and 270.26 μmol CO2 m-2 s-1 being
ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO), observed at doses of 1.6 and 1.2 g L-1, representing gains
as reported by Sá et al. (2019). equivalent to 46.9 and 18.8%, respectively (Figures 3B and
The instantaneous carboxylation efficiency (iCE) was 3C). This positive effect indicates that the foliar application
reduced with increasing salinity, with the highest value (0.026) of nitrogen promoted improvements in the assimilation and
at the ECw of 0.5 dS m-1, followed by decreases that reached internal concentration of CO2 in pomegranate seedlings. This
36.1% at the salinity of 5.00 dS m-1 (Figure 2D). This marked can be explained by the fact that N is a nutrient that acts in
reduction in iCE with the increase in salinity levels is due to the several physiological processes, being involved in the synthesis
deleterious effects on the absorption and assimilation of CO2 of chlorophyll, influencing the photosynthetic capacity of the
by pomegranate seedlings. This directly reflects the reduction plant (Cerqueira et al., 2019).
of photosynthesis, promoting metabolic changes in the Calvin The intrinsic water use efficiency (iWUE) was reduced
cycle, preventing carbon from being fixed (Sousa et al., 2016). up to the dose of 1.47 g L-1 of N, increasing at the highest
The instantaneous water use efficiency (WUE) (Z = 3.46 concentrations (Figure 3D). The occurrence of this effect
– 0.13**x - 0.08**x2 - 0.18**y - 0.26**y2 + 0.29**xy; R2 = 0.54) may be associated with the positive effect of N on stomatal
was higher at the dose of 0.03 g L -1 of N and ECw of opening, as observed in gs, favoring the diffusion of CO2 in
0.52 dS m-1. The occurrence of this effect indicates that N the cells of the leaf mesophyll, consequently facilitating the loss
promotes improvements in WUE when pomegranate plants of water through transpiration, reducing water use efficiency
are subjected to high salinity, so when the ECw increases (Talebnejad & Sepakhah, 2016).
For the interaction between factors, foliar nitrogen nitrogen improves the performance of pomegranate seedlings.
fertilization and brackish water for plant height and number This is due to the fact that this nutrient is part of several organic
of leaves, foliar nitrogen fertilization at doses of 1.28 and compounds important for plant metabolism (amino acids,
1.31 g L-1 reduced the deleterious effects of salinity up to the proteins, proline), favoring osmotic adjustment (Cerqueira et al.,
ECw of 0.52 and 1.73 dS m-1, respectively (Figures 4A and 4B). 2019). Foliar nitrogen fertilization attenuates salt stress in yellow
The attenuating effect of nitrogen on growth in height and passion fruit seedlings – Passiflora edulis (Pereira et al., 2022)
number of leaves is an indication that foliar application of and sugar apple – Annona squamosa L. (Fátima et al., 2023).
*, ** - Significant at p ≤ 0.05 and p ≤ 0.01, respectively; ns - Not significant, by F test *, ** - Significant at p ≤ 0.05 and p ≤ 0.01, respectively, by F test
Figure 3. Stomatal conductance - gs (A), net photosynthesis - Figure 4. Plant height (A), number of leaves (B) as a function of
A (B), intercellular CO2 concentration - Ci (C), and intrinsic interaction between foliar nitrogen fertilization and salinity of
water use efficiency - iWUE (D) of Punica granatum L. under irrigation water and stem diameter (C and D) of Punica granatum
foliar nitrogen fertilization L. under foliar nitrogen fertilization and salinity of irrigation water
Stem diameter was negatively influenced by salinity, with a Coulombier, N.; Nicolau, E.; Le Dean, L.; Barthelemy, V.; Schreiber,
37.9% decrease at the highest salinity (5.00 dS m-1) compared N.; Brun, P.; Lebouvier, N.; Jauffrais. T. Effects of nitrogen
to the control (Figure 4C). The reduction in water potential availability on the antioxidant activity and carotenoid content of
limits water absorption by plants, resulting in less metabolic the microalgae Nephroselmis sp. Marine Drugs, v.18, p.1-22, 2020.
activity for cell division, which explains the decrease in stem https://doi.org/10.3390/md18090453
diameter with increased salinity (Veloso et al., 2018). EMBRAPA - Manual de métodos de análise de solo. Manual de
Foliar nitrogen application promoted an increase in stem análises químicas de solos, plantas e fertilizantes. 3.ed. Brasília:
diameter up to the dose of 1.50 g L-1 of N, with an increase of Embrapa Solos, 2017. 574p.
23% or equivalent to 1.20 mm being observed when compared Fátima, R. T. de; Nóbrega, J. S.; Ribeiro, J. E. da S.; Celedônio, W.
to plants subjected to dose 0 (Figure 4D). This increase in stem F.; Ferreira, J. T. A.; Pereira, W. E.; Souto, A. G. de L.; Lima, G.
diameter may be related to the action of nitrogen in favoring S. de. Morphophysiology and quality of custard apple seedlings
the assimilation of CO2 and, consequently, promoting greater irrigated with saline water and foliar nitrogen. Journal of Plant
growth and development of plants (Basra et al., 2014). Nutrition, v.46, p.1-12, 2023. https://doi.org/10.1080/01904167.
2022.2155553
Conclusions Huang, C. J.; Wei, G.; Jie, Y. C.; Xu, J. J.; Zhão, S. Y.; Wang, L. C.; Anjum,
S. A. Responses of gas exchange, chlorophyll synthesis and ROS-
1. Foliar nitrogen fertilization at doses of 1.28 and 1.31 g L-1 scavenging systems to salinity stress in two ramie (Boehmeria
of N (Nitrotecnia-20) attenuates the effect of salinity on growth nivea L.) cultivars. Photosynthetica, v.53, p.455-463, 2015. https://
in height and number of leaves in pomegranate seedlings; doi.org/10.1007/s11099-015-0127-0
2. Irrigation water salinity above 1.7 dS m-1 drastically Lima, G. S. de; Pinheiro, F. W. A.; Dias, A. S.; Gheyi, H. R.; Silva, S. S.
reduces gas exchange in pomegranate seedlings. da; Soares, L. A. dos A.; Silva, A. A. R. da; Fernandes, P. D.; Dantas,
3. Net photosynthesis of pomegranate seedlings is J. S. Water status, cell damage and gas exchanges in West Indian
stimulated by foliar application of 1.58 g L-1 of N. cherry (Malpighia emarginata) under salt stress and nitrogen
fertilization. Australian Journal of Crop Science, v.14, p.319-324,
Acknowledgements 2020. https://doi.org/10.21475/ajcs.20.14.02.p2320
Novais, R. F.; Neves J. C. L.; Barros N. F. Ensaio em ambiente
The Coordenação de Aperfeiçoamento de Pessoal de
controlado. In: Oliveira, A. J. (ed.) Métodos de pesquisa em
Nível Superior (CAPES) and the Conselho Nacional de
fertilidade do solo. Brasília: Embrapa-SEA, 1991. p.189-253.
Desenvolvimento Cientifico e Tecnológico (CNPq) for the
Nóbrega, J. S.; Figueiredo, F. R. A.; Silva, T. I. da; Fátima, R. T.
granting of scholarships to graduate students of the work.
de; Ferreira, J. T. A.; Ribeiro, J. E. da S.; Bruno, R. de L. A.
Ecophysiology of Mesosphaerum suaveolens (L.) Kuntze
Literature Cited (Lamiaceae) under saline stress and salicylic acid. Ciência Rural,
v.59, p.1-9, 2022. https://doi.org/10.1590/0103-8478cr20210389
Alvares, C. A.; Stape. J. L.; Sentelhas, P. C.; Gonçalves, J. L. de M.;
Nóbrega, J. S.; Silva, T. I. da; Lopes, A. S.; Costa, R. N. M.; Ribeiro,
Leonardo. J.; Sparovek, G. Köppen’s climate classification map for
J. E. da S.; Silva, E. C. da; Bezerra, A. C.; Silva, A. V. da; Dias, T.
Brazil. Meteorologische Zeitschrift, v.22, p.711-728, 2013. https://
J. Foliar nitrogen fertilization attenuating harmful effects of salt
doi.org/10.1127/0941-2948/2013/0507
Basra, S. M. A.; Iqbal S.; Afzal, I. Evaluating the response of nitrogen stress on purple basil. Revista Brasileira de Engenharia Agrícola
application on growth, development and yield of quinoa e Ambiental, v.27, p.472-479, 2023. https://doi.org/10.1590/1807-
genotypes. International Journal of Agriculture & Biology, v.16, 1929/agriambi.v27n6p472-479
p.886-892, 2014. Pereira, M. B.; Nóbrega, J. S.; Fátima, R. T. de; Ferreira, J. T. A.;
Bernardo, S.; Mantovani, E. C.; Silva, D. D. da.; Soares, A. A. Manual Figueiredo, F. R. A.; Lopes, M. F. de Q.; Ribeiro, J. E. da S.; Pereira,
de irrigação. 8.ed. Viçosa: UFV, 2019. 545p. W. E. Growth and photosynthetic pigments of passion fruit
Bezerra, I. B.; Nobre, R. G.; Gheyi, H. R.; Lima, G. S. de; Barbosa, (Passiflora edulis) seedlings under foliar fertilization with nitrogen
J. L. Physiological indices and growth of ‘Paluma’ guava and irrigated with saline water. Dyna, v.89, p.60-67, 2022. https://
under saline water irrigation and nitrogen fertigation. Revista doi.org/10.15446/dyna.v89n224.100919
Caatinga, v.31, p.808-816, 2019. https://doi.org/10.1590/1983- R Core Team. R: A language and environment for statistical computing.
21252018v31n402rc Vienna: R Foundation for Statistical Computing, 2021. Available
Catola, S.; Marino, G.; Emiliani, G.; Huseynova, T.; Musayey, M.; on: <https://www.r-project.org/>. Accessed on: Jul. 2022.
Akparov, Z.; Maserti, B. E. Physiological and metabolomic analysis Richards. L. A. Diagnosis and improvement of saline and alkaline
of Punica granatum (L.) under drought stress. Planta, v.243, p.441- soils. Washington: United States Salinity Laboratory Staff, 1954.
449, 2016. https://doi.org/10.1007/s00425-015-2414-1 160p. Agriculture Handbook, 60
Cerqueira, G.; Santos, M. C.; Marchiori, P. E. R.; Silveira, N. Sá, F. V. da S.; Gheyi, H. R.; Lima, G. S. de; Paiva, E. P. de; Silva, L. de
M.; Machado, E. C.; Ribeiro, R. V. Leaf nitrogen supply A.; Moreira, R. C. L.; Fernandes, P. D.; Dias, A. S. Ecophysiology of
improves sugarcane photosynthesis under low temperature. West Indian cherry irrigated with saline water under phosphorus
Photosynthetica, v.57, p.18-26, 2019. http://doi.org/10.32615/ and nitrogen doses. Bioscience Journal, v.35, p.211-221, 2019.
ps.2019.033 https://doi.org/10.14393/BJ-v35n1a2019-41742
Sakazaki, R. T.; Araújo, W. F.; Monteiro Neto, J. L. L.; Chagas, P. C.; Sousa, J. R. M. de; Gheyi, H. R.; Brito, M. E. B.; Xavier, D. A.; Furtado, G.
Chagas, E. A.; Murga-Orrillo, H.; Bardales-Lozano, R. M.; Abanto- de F. Impact of saline conditions and nitrogen fertilization on citrus
Rodriguez, C. Shade nets and substrates in seedling production of production and gas exchanges. Revista Caatinga, v.29, p.415-424,
Annona squamosa L. in the Roraima Cerrado. Semina: Ciências 2016. https://doi.org/10.1590/1983-21252016v29n218rc
Agrárias, v.40, p.2535-2544, 2019. https://doi.org/10.5433/1679- Talebnejad, R.; Sepaskhah, A. R. Physiological characteristics, gas exchange,
0359.2019v40n6p2535 and plant ion relations of quinoa to different saline groundwater depths
Silva, A. A. R. da; Lima, G. S. de; Azevedo, C. A. V. de; Gheyi, H. R.; and water salinity. Archives of Agronomy and Soil Science, v.62, p.1347-
Souza, L. de P.; Veloso, L. L. de S. A. Gas exchanges and growth of 1367, 2016. https://doi.org/10.1080/03650340.2016.1144925
passion fruit seedlings under salt stress and hydrogen peroxide. Veloso, L. L. de S. A.; Nobre, R. G.; Souza, C. M. A. de; Fátima, R. T. de; Souza,
Pesquisa Agropecuária Tropical, v.49, p.1-10, 2019. https://doi. L. de P.; Elias, J. J.; Azevêdo, F. L. de; Santos, J. B. dos. Morphophysiology
org/10.1590/1983-40632019v4955671 of guava cv. Paluma with water of different salt concentrations and
Soares, L. A. dos A.; Oliveira, S. G. de; Lima, G. S. de; Fernandes, P. proline doses. Semina: Ciências Agrárias, v.39, p.1877-1886, 2018.
D.; Araújo, R. H. C. R.; Fernandes, E. A. Physiological changes of https://doi.org/10.5433/1679-0359.2018v39n5p1877
pomegranate seedlings under salt stress and nitrogen fertilization. Wang, X.; Wang, L.; Shangguan, Z. Leaf gas exchange and fluorescence
Revista Brasileira de Engenharia Agrícola e Ambiental, v.25, of two winter wheat varieties in response to drought stress
p.453-459, 2021. http://dx.doi.org/10.1590/1807-1929/agriambi. and nitrogen supply. Plos One, v.11, p.1-15, 2016. https://doi.
v25n7p453-459 org/10.1371/journal.pone.0165733