Exogenously Applied Proline Enhances Morph-Physiological Responses and Yield of Drought-Stressed Maize Plants Grown Under Different Irrigation Systems

The exogenous application of osmoprotectants [e.g., proline (Pro)] is an important approach for alleviating the adverse effects of abiotic stresses on plants. Field trials were conducted during the summers of 2017 and 2018 to determine the effects of deficit irrigation and exogenous application of Pro on the productivity, morph-physiological responses, and yield of maize grown under two irrigation systems [surface irrigation (SI) and drip irrigation (DI)]. Three deficit irrigation levels (I100, I85, and I70, representing 100, 85, and 70% of crop evapotranspiration, respectively) and two concentrations of Pro (Pro1 = 2 mM and Pro2 = 4 mM) were used in this study. The plants exposed to drought stress showed a significant reduction in plant height, dry matter, leaf area, chlorophyll content [soil plant analysis development (SPAD)], quantum efficiency of photosystem II [Fv/Fm, Fv/F0, and performance index (PI)], water status [membrane stability index (MSI) and relative water content (RWC)], and grain yield. The DI system increased crop growth and yield and reduced the irrigation water input by 30% compared with the SI system. The growth, water status, and yield of plants significantly decreased with an increase in the water stress levels under the SI system. Under the irrigation systems tested in this study, Pro1 and Pro2 increased plant height by 16 and 18%, RWC by 7 and 10%, MSI by 6 and 12%, PI by 6 and 19%, chlorophyll fluorescence by 7 and 11%, relative chlorophyll content by 9 and 14%, and grain yield by 10 and 14%, respectively, compared with Pro0 control treatment (no Pro). The interaction of Pro2 at I100 irrigation level in DI resulted in the highest grain yield (8.42 t ha–1). However, under the DI or SI system, exogenously applied Pro2 at I85 irrigation level may be effective in achieving higher water productivity and yield without exerting any harmful effects on the growth or yield of maize under limited water conditions. Our results demonstrated the importance of the application of Pro as a tolerance inducer of drought stress in maize.


INTRODUCTION
Maize is one of the main cereal crops worldwide, followed by wheat, which is a staple food in Egypt (OECD/FAO, 2020). It is either consumed freshly or indirectly as corn oil, starch, fructose, glucose, and livestock feed (Ranum et al., 2014). Owing to its low inherent production and high demand of maize, Egypt has become the fourth largest importer of maize in the world (OECD/FAO, 2020). The total maize production in Egypt is approximately 6.4 million tons, and 85,000 hectares of cultivated land is estimated for its production. Its consumption is approximately 16.1 million tons with a self-sufficiency ratio of 42% (Ibrahim et al., 2007). Thus, to minimize the gap between maize production and consumption, it is important to manage irrigation water more efficiently and enhance domestic maize production by following non-traditional procedures, such as growing in cultivated areas or planting high-yield varieties (Ouda et al., 2017).
In arid and semi-arid regions, water shortage particularly affects food security (Misra, 2014;Mancosu et al., 2015). Approximately 70-80% of the available freshwater is required for agriculture (Food and Agriculture Organization [FAO], 2017). Egypt depends on limited water resources from the Nile River (55 billion m 3 of water/year) (MWRI, 2014). Over 95% of irrigated lands in Egypt are managed using surface irrigation (SI) systems, with a low irrigation efficiency of 45-50% (Osman et al., 2016;Abdou et al., 2021;Ali et al., 2021). Therefore, rational water governance seeks to reduce water losses and enhance crop productivity to withstand high evaporative demand. Assuming severe water supply shortages, a deficit and highly efficient irrigation strategy of drip irrigation (DI) system is highly recommended over the SI system to increase the benefits per unit of water.
The application of irrigation water below the evapotranspiration (ET) demand is known as deficit irrigation, which optimizes economic output when there is limited water supply (Pereira et al., 2012;Chai et al., 2016). Plants under deficit irrigation systems receive less irrigation water than the actual amount of water required at plant growth stages and/or during the total crop cultivation period (Badal et al., 2013). Hence, plants are exposed to water stress to some extent under deficit irrigation systems (Wakchaure et al., 2018;Parkash and Singh, 2020). A decrease in the ET rate of plants exposed to water stress results in severe water stress symptoms, such as leaf rolling, diminishing leaf area, and stunted growth and yield of plants (Wang et al., 2014). Drought stress induces several physiological, biochemical, and photosynthetic changes by regulating stomatal closure or controlling CO 2 flow into the mesophyll tissues (Yuan et al., 2016).
Approximately 20-25% of the maize cultivation areas are affected by drought globally (Shahrokhi et al., 2020). Sensitivity to drought stress can lead to dramatic fluctuations in maize yield, which is a common condition in Egypt (Gomaa et al., 2017). Maize exhibits distinct responses to water deficit depending on the developmental stage, crop tolerance level, and severity of the applied water stress treatments (Song et al., 2019;Cheng et al., 2021). The yield reduction due to water stress is primarily attributed to the disruption of physiobiochemical processes, inhibition of photosynthesis, and stunted plant growth and development (Song et al., 2019;Sharma et al., 2020;Cheng et al., 2021).
However, several studies have revealed that plants tend to accumulate various compatible solutes, such as soluble sugars and amino acids [e.g., proline (Pro)], as an adaptive tolerance strategy to increase salinity and/or drought tolerance (Giri, 2011;Semida et al., 2015;Arteaga et al., 2020). Pro is an essential amino acid that accumulates in various plant tissues, particularly in the leaves (Hayat et al., 2012). Proline accumulation plays an indispensable role in the regulation of cell osmosis (El Moukhtari et al., 2020). Additionally, under unfavorable conditions of water stress, Pro stabilizes membranes (Xia et al., 2020) and protects enzymes (Qamar et al., 2019;El-Nashaar et al., 2020), proteins, and macromolecules from denaturation (Agami et al., 2019). The effectiveness of exogenously applied Pro depends on the plant developmental stage and variety, Pro application rate, and the time of application (Semida et al., 2020). Excessive concentrations of free Pro can disrupt cell growth and protein functions. Hayat et al. (2007) observed that the foliar application of Pro (30 mM) negatively affects maize plants.
Using Pro as an antioxidant is an effective, unconventional, and inexpensive approach for mitigating the harmful effects of drought on maize production under different irrigation systems. Further studies are highly recommended to reduce the gap between water consumption and production of maize in arid and semi-arid regions. Therefore, this study aimed to elucidate the ameliorative effects of Pro at different nontoxic application rates on the morph-physiological and yield attributes of maize plants exposed to water stress. Furthermore, irrigation water requirements for maize plants grown under different irrigation systems and water scarcity conditions were determined.

Experimental Site
Field experiments were performed in the summers of 2017 and 2018 at the Agricultural Research Station in the Fayoum University, Fayoum Governorate, Egypt (latitude, 29 • 02 and 29 • 35 N and longitude, 30 • 23 and 31 • 05 E). The experimental treatments were arranged in a spilt-split plot in a completely randomized block design with three replications. The two irrigation systems, SI and DI, were set in the main plots. Three deficit irrigation levels (I 100 , I 85 , and I 70 of ET) were employed in the subplots, and three concentrations of Pro (Pro 0 = 0/control; Pro 1 = 2 mM; and Pro 2 = 4 mM) were used in the sub-sub plots.

Initial Soil Characteristics of the Experimental Site
The soil in the field experiment was sandy loam; the water content at −0.33 and −15 bar pressure was retained by 20.52 and 8.87%, respectively. At a soil depth of 0.0-60 cm, the dry bulk density of the soil was 1.47 mg m −3 , saturated hydraulic conductivity was 2.34 cm h −1 , mean soil pest extract (ECe) was 3.75 dS m −1 , and pH (soil:water suspension, 1:2.5) was 7.77. At the same soil depth, the mean organic matter and CaCO 3 concentration were 13.8 and 87.8 g kg −1 , respectively ( Table 1).
The initial physical and chemical characteristics of the soil at the experimental site were determined according to the methods proposed by Page et al. (1982) and Klute and Dirksen (1986).

Irrigation Systems
The SI and DI systems were implemented and installed at the experimental site as described in the following section.

SI System
Under the SI system, the quantity of irrigation water applied (IWA) for each plot was mainly controlled using a plastic pipe (spiles) with a diameter of 2 inches. For each experimental plot, one spile was constructed to direct the irrigation water. The quantity of IWA was estimated using the following equation (Israelsen and Hansen, 1962): where Q is the irrigation water discharge (L s −1 ), C is the discharge coefficient, A is the cross-sectional area of the spile (cm 2 ), g is gravity acceleration (cm s −2 ), and h is the mean of the influential head of water (cm) above the pipe.

DI System
Under the DI system, irrigation water was provided through polyvinyl chloride pipes. A distance of 70 cm was maintained between the lateral lines, and emitters were spaced at 30 cm. An emitter discharge rate of 2 L h −1 was achieved with pumping at a pressure of 2 bar.

Water Stress Treatments
Three different water stress levels [I 100 , I 85 , and I 70 (%) of ET c ] were implemented. Maize plants were irrigated with the corresponding amounts of irrigation water, which were determined according to the daily reference crop ET as follows (Allen et al., 1998): where E pan is the evaporation rate from Class A pan (mm day −1 ) and K pan is the pan evaporation coefficient (0.85). ET c is the sum of the water evaporation from the soil surface and transpiration (water loss) primarily from the plant leaves. ET o is the reference ET. K c is the crop coefficient; the K c values at the initial, mid, and end stages were 0.70, 1.20, and 0.35, respectively. The irrigation water requirements (IWR) for each plot were determined using the following equation (Abd El-Wahed and Ali, 2013): where IWR is the water requirement for irrigation (m 3 ), A is the irrigated plot area (m 2 ), ET c is the crop ET (mm day −1 ), Ii is the interval between irrigations (days), Kr is the coverage coefficient, and E a is the irrigation efficiency (%).

Antioxidant Applications
The antioxidant Pro was applied at three different concentrations : Pro 0 (control) = 0, Pro 1 = 2 mM, and Pro 2 = 4 mM. It was applied three times as a foliar spray [30, 45, and 60 days after sowing (DAS)].

Meteorological Parameters
The monthly weather data for the summers of 2017 and 2018 were collected from the meteorological station of the Fayoum Governorate, Egypt (Figure 1). The climate in the experimental field was arid and characterized by low or no precipitation. During the period from May to September of each year, the maximum and minimum air temperatures were 40.4 • C and 20.15 • C, respectively, and the relative humidity fluctuated from 35 to 42%. The mean E pan values during the cultivation period (May-September) were 6.48 and 6.50 mm day −1 for the first and second growing seasons, respectively.

Agricultural Management Practices
Maize (Zea mays, hybrid 321) was cultivated on June 5 during the summers of both years (2017 and 2018) in hills that were 30 cm apart, with a distance of 70 cm between the lateral lines. The experimental plots were fertilized using N:P:K (200:100:75; kg ha −1 ). Total superphosphate (15.5% P 2 O 5 ) along with potassium sulfate (48% K 2 O) was applied before cultivation. Nitrogen fertilizer was applied using two equivalent doses, with the first and second doses applied at the first and second irrigation events, respectively. All other agricultural operations required for the growth and development of maize plants were performed similarly in all plots

Physiological Maize Parameters
Ten plant samples (70 DAS) were randomly selected from each plot to measure their physiological responses to the treatments applied.

Relative Water Content
The leaf samples for relative water content (RWC) measurement were randomly collected in the morning (8:00 a.m.-9:00 a.m.). RWC (%) was determined according to the following equation (Hayat et al., 2007): where FW is the fresh weight measured within 2 h after the excision of leaves; TW is the turgid weight computed by soaking the leaves in distilled water and leaving them at room temperature for 16-18 h, followed by rapid and careful dry-blotting on tissue paper. The small leaf pieces were later oven-dried at 70 • C for 48 h to assess the dry weight (DW).

Leaf Membrane Stability Index
MSI (%) was measured using the method described by Premachandra et al. (1990). Small leaf strips (0.2 g) of equal size were prepared and placed in two sets of test tubes, each containing 10 mL of distilled water. The test tubes of the first set were incubated in a water bath at 40 • C for 30 min, and ECe was subsequently estimated (C 1 ), whereas those of the second set were incubated in a boiling water bath at 100 • C for 15 min, followed by ECe measurement (C 2 ). MSI = [1− C1 C2 ] × 100

Relative Chlorophyll Content Values
The relative chlorophyll content SPAD was determined using SPAD 502 (Konicaminolta. Inc., Tokyo), as described by Maxwell and Johnson (2000).

Performance Index
The performance index (PI) of photosynthesis was determined according to the method proposed by Clark et al. (2000). During harvesting (120 DAS), ten plants from each plot were collected to determine their height (cm), stem diameter (cm), leaf number plant −1 , root weight, cob weight (g), 100-grain weight (g), grain yield (t ha −1 ), and biomass yield (t ha −1 ).

Water Productivity
Water productivity (WP) is expressed as the grain yield (kg) per IWA (m 3 ). The values were used to evaluate the variation in different treatments for producing the maximum yield from the water unit consumed by the maize plants. The WP values were calculated (Jensen, 1983) as per the following equation: Irrigation water applied (m 3 ha −1 )

Statistical Analysis
The study was designed in a completely randomized block design (spilt-split plot) with three replications, and the obtained data were statistically analyzed according to the procedures outlined by Gomez and Gomez (1984) using the GenStat statistical package, 12th edition (VSN International Ltd., Oxford, United Kingdom).

Plant Water Status (RWC and MSI)
As shown in Table 3, the two indicators of plant water status, RWC (%) and MSI, showed a significant response to water stress and Pro treatments under the two irrigation systems during the summers of both years, whereas the combined effect of these treatments (IS × I × P) on the aforementioned indicators was not significant. The results showed that maize plants under the DI system exhibited better physiological responses than those under the SI system. Additionally, water stress considerably reduced RWC by 10.51 and 17.38% in maize plants under moderate (I 85 ) and severe (I 70 ) stress levels, respectively, compared with those under no stress (controls; I 100 ). The exogenous application of Pro mitigated the drought-induced inhibitory effects on RWC. Compared with the control (Pro 0 ), the mean values of RWC increased by 15.51 and 21.25% with increased concentrations of Pro 1 and Pro 2 , respectively. Furthermore, the highest mean MSI (67.98%) was observed in adequately irrigated plants (100% ET c ), and the lowest mean MSI (59.69%) was observed in plants exposed to high water deficit levels (I 70 ). Cell membrane integrity is highly susceptible to drought stress. Water stress leads to the disruption of cell membrane stability, thus increasing electrolyte leakage and decreasing membrane integrity. The mean MSI increased by 6.88 and 13.25% for Pro 1 -and Pro 2 -treated plants, respectively, compared with that for non-Pro-treated plants during the summers of both years.

Fv/Fm, PI, and Soil Plant Analysis Development
The data presented in Table 3 revealed that maize plants showed significant differences in Fv/Fm, PI, and SPAD in response to water stress and Pro foliar treatments under the SI and DI systems during the summers of both years; however, the effect of their interaction was not significant. The values for these traits increased by 7.64, 17.06, and 21.26% for Fv/Fm, PI, and SPAD, respectively, under the DI system compared with the SI system. Water stress adversely affects photosynthesis. The Fv/Fm, PI, and SPAD values showed a consistent decrease with increasing water deficit. Compared with sufficiently irrigated plants (I 100 ), the Fv/Fm, PI, and SPAD values decreased by 3.82, 17.25, and 9.59%, respectively, in plants grown under the moderate-deficit irrigation level (I 85 ) and by 9.55, 28.01, and 21.36%, respectively, in those under the high-deficit irrigation level (I 70 ). Furthermore, this study indicated that the harmful effects of drought stress on photosynthesis can be alleviated in maize through foliar Pro application. The exogenous application of Pro (Pro 2 concentration) on maize plants exposed to water stress increased the Fv/Fm, PI, and SPAD values by 12.06, 20.71, and 15.27%, respectively, compared with the control (Pro 0 ). Table 4 showed that the yield and yield components of maize, such as cob weight (g), 100-grain weight (g), grain yield (t ha −1 ), and biomass yield (t ha −1 ), differed significantly in response to water stress and Pro applications under the SI and DI systems during the summers of both years. The effects of interaction were not significant between (IS × I) and (IS × P) for the 100grain weight parameter. Except of grain yield, the interaction effect (IS × IR × P) on all abovementioned traits was not significant (P > 0.05). Maize plants exposed to water stress had a lower yield than fully irrigated plants. Compared with the full irrigation level (I 100 ), the cob weight (g), 100-grain weight (g), grain yield (t ha −1 ), and biomass yield (t ha −1 ) of maize plants decreased by 4.33, 4.45, 6.36, and 9.24%, respectively, at the moderate-deficit irrigation level (I 85 ) and by 7. 49, 14.90, 25.44, and 17.97% at the high-deficit irrigation level (I 70 ). The exogenously applied Pro increased the yield of maize plants under deficit irrigation; the cob weight, 100-grain weight, grain yield, and biomass yield increased (as mean values of both seasons) by 11.70, 21.07, 10.93, and 15.44%, respectively, at moderate Pro (Pro 1 ) application rates and by 15.61, 27.89, 15.99, and 23.05%, respectively, at high Pro (Pro 2 ) application rates, compared with the control (Pro 0 ). Table 5 showed that IWA varied according to the designed irrigation regime and system. The highest values of irrigation water inputs were observed under the SI system at the full irrigation level. The need for irrigation gradually decreased as the water deficit level increased in the DI system because its high-efficiency design saved 29.41% of IWA compared with the SI system.

Water Productivity
As shown in Figure 2, the highest WP (1.80 kg m −3 ) was observed at I 85 with a high Pro application rate (Pro 2 ) under the DI system. In contrast, the lowest (0.79 kg m −3 ) WP was observed at I 100 with a low Pro application rate (Pro 1 ) under the SI system. For each irrigation level and Pro application treatment, DI considerably enhanced WP compared with SI. DI delivers water near each plant, leading to high water-use efficiency. Among the water stress treatments, the irrigation regime of I 85 , in which the reduction in grain yield was lower than the amount of irrigation water saved, showed the highest WP-value compared with other irrigation regimes.  movement or migration of soluble salts, thus reducing the concentrations of salt accumulated in the lower soil layers by 10.14% compared with those accumulated on the soil surface (0-20 cm). In contrast, DI with low irrigation input (I 70 ) promoted the upward movement of salts, thus increasing soil salinity in the same layer (0-20 cm) by 14.38% compared with the initial soil salinity.

DISCUSSION
Water stress adversely affects the growth and development of plants, particularly when the stress conditions remain constant. When plants are exposed to stress, they accumulate an array of metabolites, such as amino acids. Pro plays a highly beneficial role in plants exposed to various stress conditions. Besides acting as an excellent osmolyte, Pro has three fundamental biological roles in stress response, i.e., acting as a metal chelator, antioxidative defense molecule, and signaling molecule. Some reports; however, have indicated the toxic effects of Pro when applied exogenously at high concentrations. In this study, we aimed to (i) investigate the causes underlying the reduction in the growth and yield of maize plants induced by water stress and (ii) determine the crucial role of Pro in mitigating these negative effects at different nontoxic application rates of Pro to improve WP under water scarce conditions. Our results demonstrated significant differences in the vegetative growth of maize plants in response to the two irrigation systems ( Table 2). The DI system was relatively superior in the maize growth parameters, except for root weight, compared with the SI system. Similar to previous reports (Ghamarnia et al., 2013;Mahgoub et al., 2017), the observed increase in the root weight of plants under the SI system could be a result of the downward movement of irrigation water due to gravity, which promotes the penetration of roots through the soil profile to greater depths to extract water from the deeper layers. Maize plants exposed to water stress showed growth retardation. The inhibition of maize growth due to water stress was evident because of the decrease in the length, volume, and water potential of the root system, which disrupted the water extraction by the roots. This results in decreased cell division and elongation, leading to stunted plant growth. These results were consistent with those of previous studies (El-Samnoudi et al., 2019;Song et al., 2019;Sah et al., 2020;Cheng et al., 2021). However, the exogenous application of Pro not only provided osmoprotection but also enhanced the growth of the plants. The beneficial effect of Pro in enhancing maize growth characteristics under low irrigation levels might be attributed to its role in osmoregulation, and maintenance of the tertiary structure of proteins and enzymes, to help growing plants tolerate drought stress. These results were similar to those of the previous reports (Abdelaal et al., 2018Qamar et al., 2019;El-Nashaar et al., 2020).  Additionally, drought adversely affected the plant water retention (RWC and MSI) and photosynthesis ( Table 3). The decrease in water uptake under deficit irrigation regimes was associated with a decrease in leaf water potential. Drought alters the plant water status and stomatal functioning and inhibits chlorophyll biosynthesis, thus reducing the photosynthetic rates. Dehydration in the plant cell protoplasm could be considered as an effect of drought on the RWC of maize leaves. Similar observations were reported in a previous study (Xia et al., 2020). The foliar application of Pro enhanced the plant water status. This was evident by the reduction in water efflux under drought stress, which helped the cells maintain their cell turgor or osmotic balance. These results are consistent with those of another study conducted by Al-Khazrji et al. (2020). It is suggested that the exogenous application of Pro enhances membrane stability, thus preventing electrolyte leakage, as previously documented by Hayat et al. (2012).
Water stress further decreased the values of the physiological parameters in maize seedlings. This was evident by the disruption of plant photosynthetic potential, possibly due to the stomatal closure and/or metabolic destruction, such as the impairment of photosystem 1 (PSI) and photosystem II (PSII), which are chlorophyll-binding protein complexes Zheng et al., 2022). Furthermore, water stress increases the generation of reactive oxygen species (ROS), resulting in oxidative damage to plants and degradation of chlorophyll pigments (Hasanuzzaman et al., 2020;Mansoor et al., 2022). Alternatively, Pro may be responsible for scavenging ROS and other free radicals (Kaul et al., 2008). It may also protect plants from stress injuries by stabilizing membranes and proteins, allowing the transport of mitochondrial electrons, enhancing antioxidant enzyme activity, thus increasing stomatal conductance and facilitating higher CO 2 diffusion through leaves. This promotes higher photosynthetic capacity (Tùmová et al., 2018;Agami et al., 2019;Altuntaş et al., 2020;Semida et al., 2020;Tariq et al., 2021;Abdou et al., 2022). Hence, the foliar application of Pro can enhance the photosynthetic rate of maize under drought stress.
The decreased water uptake from the soil might be responsible for the detrimental effects of water stress on grain and biomass yields of maize (Table 4), resulting in the abovementioned inhibition in plant growth and development, disruption of photosynthetic pigments, and deficits in plant water content. Accordingly, a significant decrease in maize yield under water stress conditions was observed. The results are consistent with those of previous studies (Vazirimehr et al., 2014;Sah et al., 2020;Shah et al., 2020).
The present study demonstrated that the increase in the yield and yield components achieved by the application of Pro was higher than the reduction in the yield caused by drought stress, confirming that Pro treatment successfully compensated for the adverse effects of drought stress on the growth and yield of maize. Furthermore, the yield and yield components of maize under the DI system were higher than those under the SI system (Qamar et al., 2019;Abdelaal et al., 2020;El-Nashaar et al., 2020;Li et al., 2021), and WP increased consistently with increasing concentrations of exogenously applied Pro (Table 5). These results are consistent with those of the previous studies by Semida et al. (2015) and Semida et al. (2020) who found that pressurized irrigation systems via subsurface drip and/or surface drip increase the water-use efficiency and grain yield of maize compared with the SI system.
Notably, the monitoring of the salt distribution pattern in the soil profile (depth, 0-60 cm) based on different applied irrigation systems and treatments showed that the salt distribution pattern was different after plant harvest in the second season compared with the determined initial ECe values before maize planting. The results were consistent with those of the studies by Semida et al. (2015) and Ali et al. (2016) who reported that the highest salinity was observed in the deeper soil layers depending on the wetting front under the SI system. The lowest salinity was observed on the surface; however, under DI, salt accumulation increased on the soil surface.

CONCLUSION
With increasing water scarcity, well-designed deficit DI regimes and the exogenously applied antioxidant Pro can optimize the maize production and WP when the available water is insufficient to provide full irrigation. The DI system efficiently saved irrigation water input by 30% compared with the SI system. Among the tested irrigation levels, the full level (I 100 ) resulted in the highest grain yield. The moderate-deficit irrigation level (I 85 ) showed the maximum WP. The growth, physiological aspects, and grain yield of maize increased significantly after the foliar application of Pro to maize plants exposed to water stress. Therefore, using exogenous Pro at 4 mM for maize plants under the DI system and irrigating with 85% of their ET demand may be a promising agro-management strategy for improving the yield and WP of maize crops grown in arid and semi-aired regions, especially when water scarcity constraints the sustainability of maize production.

DATA AVAILABILITY STATEMENT
The original contributions presented in this study are included in the article, further inquiries can be directed to the corresponding author/s.

ACKNOWLEDGMENTS
KE-T would like to thank the library at Murdoch University, Australia, for providing valuable online resources and comprehensive databases.