Effects of gibberellic acid on Tifton 85 bermudagrass (Cynodon spp.) in constructed wetland systems

This study aimed to evaluate 1) the influence of gibberellic acid (GA3) in the development of Tifton 85 bermudagrass grown in constructed wetland systems (CWs) and 2) the plant's capacity to remove nutrients and sodium from synthetic municipal wastewater (SMW). The experiment was carried out in Viçosa, Minas Gerais, Brazil, and consisted of foliar applications of GA3 set in randomized blocks design, with four replicates and 6 treatments as following: NC (control with plants); 0 μM GA3; N1: 5 μM GA3; N2: 25 μM GA3; N3: 50 and N4: 100 μM GA3 per CWs, NC* (control with no plants): 0 μM GA3. The study was conducted over two crop cycles in the spring 2016. The parameters used to evaluate the performance of the Tifton 85 bermudagrass were its plant height, productivity, chlorophyll measurement, number of internodes, nutrients and Na removals. Chemical analyses of the effluents were conducted. In response to the application of GA3, the increase in height of Tifton 85 bermudagrass in the first crop cycle was higher than the increase in height in the second crop cycle. The decrease in plant growth in response to GA3 in the second crop cycle may be linked to the age of the plant tissue and climatic conditions. The greater growth of the plants cultivated in the CWs allows a more efficient removal of pollutants, using simple management and low cost. The results suggest that applying 50 μM of GA3 to the development of Tifton 85 bermudagrass provides higher dry matter yield and removal of nitrogen, phosphorus, and sodium for the first crop cycle in CWs. However, in the second crop cycle, the application of GA3 had no effect on dry matter production and nutrient removal by Tifton 85 bermudagrass in CWs.


Introduction
The disposal of wastewater into rivers, lakes and other water bodies leads to a reduction in quality of water resources, which are of prime importance to human and aquatic life. The discharge of wastewater contains high levels of organic nitrogen and phosphorus compounds that enable growth of algae and aquatic plants, which are responsible for eutrophication of water bodies. Constructed wetland systems (CWs) are treatment systems planned and designed to use natural processes involving wetland vegetation, soils, and their associated microbial assemblages to improve water quality [1]. In CWs-under appropriate environmental conditionsphysical, chemical and biochemical processes occur to treat sanitary sewage, industrial effluents and many types of wastewater [2]. Many countries have adopted these systems [3] because they require low-cost investments combined with their efficiency and feasibility [4].
The efficiency of sewage treatment using CWs is obtained through the hydraulic interaction of the substrate, plants, microorganisms and solar radiation. Pollutant removal takes place in a support media, which is responsible for the filtration and the formation of the biofilm next to the roots of the plants. The microorganisms adhere to the biofilm, degrade carbonaceous and nitrogenous matter, and the plants absorb the nutrients and heavy metals, thus promoting phytoremediation [5].
CWs have a high removal efficiency and appropriate effluent quality for final disposal, especially regarding organic pollutants (biochemical oxygen demand, chemical oxygen demand (COD), volatile suspended solids, and oils and greases). However, their efficiency of nitrogen (N) and phosphorus (P) removal is relatively low [3].
The selection of wetland plants is crucial to the success of the wastewater treatment because the plants uptake macro and micronutrients from the CWs for their development. The macronutrients N and P are essential for plant growth and development. Nitrogen, besides constituting proteins, is required to produce the molecule chlorophyll. P is important for the energy reserve and structural integrity of tissues.
Reed (Phragmites australis), cattail (Typha spp.), lakeshore bulrush (Schoenoplectus lacustris) and Canna indica are some of the most common wetland species studied in treatment wetlands [6]. Ornamental species [7] or species of zootechnical interest such as elephant grass (Pennisetum purpureum) and Tifton 85 bermudagrass (Cynodon spp.) [8] have been used in previous studies for economic return. Tifton 85 bermudagrass is a cross between bermudagrass (Cynodon dactylon) and a close tropical relative called stargrass (Cynodon nlemfuensis) with high zootechnical and economic potential due to its good palatability, digestibility, and expressive crude protein concentrations. In addition, it presents an attractive commercial value for the hay production, reaching higher values in regions with little forage availability.
Tifton 85 bermudagrass has interesting characteristics for wastewater treatment in CWs, as it is a perennial crop that presents a high growth rate, which in turn allows frequent cuts making it possible to extract large amounts of nutrients from the system [9].
One of the main mechanisms responsible for the removal of N in CWs is the uptake by plants when frequent cuts are performed [10]. Therefore, it is important to create favorable conditions for plants, such as adequate climatic conditions, plant density, cutting height, and the use of cultural practices that stimulate growth and development of plants to potentiate the extraction of nutrients in the CWs, thus generating a better-quality effluent.
In this context, since the nutrient removal is due to the fast plant growth and dry matter production, it is necessary to look for alternatives that potentiate the growth rate of the plants.
The growth rate of plants can be increased with the use of growth-regulating substances [11]. Gibberellic acid (GA 3 ) is one of the classes of growth phytohormones that influence plant growth, increase elongation and cell division, which is evidenced by increased length and number of cells [12]. However, the effect of GA 3 is influenced by several factors, such as environmental conditions, GA 3 concentration, number of applications, and the season of application and species or cultivar [11].
One of the roles played by plants in CWs is nutrient removal. Nutrient uptake affects plant growth rate. However, there are no references in the scientific literature regarding the effect of GA 3 on plants grown in CWs. Therefore, it is necessary to investigate the effects GA 3 has on productive characteristics of Tifton 85 bermudagrass, and pollutant removal efficiency in CWs.
This study aimed to evaluate 1) the influence of GA 3 in the development of Tifton 85 bermudagrass grown in CWs; and 2) the plant's capacity to remove nutrients and sodium (Na) from synthetic municipal wastewater (SMW).

Characterization of the experimental area
The experiment was carried out in a greenhouse located at the experimental area of the water resources reference center belonging to the Department of Agricultural Engineering at the Federal University of Viçosa, (altitude: 650 m; latitude 20˚45'14'' S, longitude 42˚52'53'' W) Viçosa, Minas Gerais state, Brazil. The climate, according to the Köppen climate classification, is Cwa (humid subtropical climate with dry winter and hot summer).
The experimental system consisted of a set of 10 L polyethylene containers (reactors) with 26.50 cm height, 23.58 cm diameter, with a surface area of 0.0437 m 2 , cultivated with Tifton 85 bermudagrass (Cynodon dactylon Pers. x Cynodon nlemfuensis Vanderyst, Viçosa, Brazil), which were then submitted to SMW applications (Fig 1).
The CWs were filled up to 21.2 cm height with gneissic gravel 0 (D60 = 7 mm, coefficient of uniformity D60/D10 = 1.3, and initial void volume of 45.8%) and a free edge of 5.3 cm. Prior to filling the CWs, the rocks were immersed in a bleach solution and rinsed with running water.

Experimental design and treatments
The different treatments and their respective doses of GA 3 are listed in Table 1. A randomized complete block design was used, with four replications and five treatments with Tifton 85 bermudagrass exposed to different doses of GA 3, and one treatment control with no plants to assess the effluent, totaling 24 experimental plots.

Synthetic municipal wastewater and operation of constructed wetland systems
The SMW was used in the experiment to guarantee greater control of chemical and organic compounds, and for sanitary purposes, since the operation of CWs and handling of SMW was done manually.
The SMW used in this study was modified from Nopens [13]; the concentrations of COD, N and P were consistent with the ranges recommended by Von Sperling for municipal wastewater [14], and the concentration of potassium (K) was around the values described by Santos [15]. The ingredients were diluted in drinking water. The compounds used to prepare SMW are shown in Table 2.
Municipal wastewater is characterized by a high concentration of salts, especially K and Na. The concentrations of K and Na in the influent of the present study are similar to those found at the municipal wastewater treatment plant of Janaúba, Minas Gerais [15]. The concentration of Na in raw sewage ranges between 30 to 50 mg L -1 [16].
The CWs worked as a batch system, being fed from the top. Influent SMW was applied until the reactors were filled up to the point where the upper layer of gravel was reached. Then the volume of input of each CWs was recorded. SMW remained in reactors during a cycle time, then the SMW was collected from the bottom to characterize a vertical flow.
SMW was changed every 7 days (cycle time) during the acclimatization period of the experiment, and every 3.5 days after that. The cycle time in batch systems is analogous to the hydraulic residence time (HRT) in continuous systems.  The salts and the other ingredients were diluted into drinking water to prepare SMW; weekly evaluations of the parameters of the influent were conducted, totaling 11 analyses. The averages and standard deviations of the actual SMW composition in the present study are shown in Table 3.
After the acclimatization period, the differentiation of treatments with HRT of 3.5 d began and the average rates of SMW application were calculated (Table 4).

Experiment development and application of gibberellic acid (GA 3 )
Tifton 85 bermudagrass sprigs, formed by a segment of stolon with two buds with leaves and roots removed, were collected and standardized at 0.10 m. Then four sprigs were planted in each reactor.
The stabilization period lasted 48 days. At the end of this period, the plants were cut at a standard height of 0.10 m above the level of the support media. 0.1% (v/v) A spreader was added to the different doses of GA 3 (2 drops of spreader in 25 ml GA 3 solution) to facilitate the Table 2. Compounds and their theoretical concentrations to prepare 1 liter of synthetic municipal wastewater. Source: Adapted from Nopens [13]. contact and absorption of the phytohormone by plants. The solutions were applied on the day of the standardization cut and at the end of the first cycle. The first and second cycles lasted 18 and 23 days respectively and were carried out in the spring 2016 (southern hemisphere). The end of each cycle was pinpointed by when the Tifton 85 bermudagrass started lodging. Both cycles received two applications of the different doses of GA 3 ; the second application was seven days after the plants being cut. The experiment was carried out for 89 days.

Salts
The weather conditions (Fig 2) were measured with two thermohygrometers installed in the center of the greenhouse. Maximum and minimum and average temperatures were recorded to monitor the weather changes in all the stages of the experiment.
The weather conditions inside the greenhouse were not controlled, so in the first and second crop cycles, the average temperatures were 27.9 and 30.2˚C, and the relative humidity was 58% and 56%, respectively.

Experimental evaluations
Yield and vegetative analyses. The parameters analyzed in Tifton 85 bermudagrass and their methodologies are shown in Table 5. The forage yield was evaluated in two cycles when were performed cuts to 0.10 m at the level of the support media. The aerial parts of the material were collected and taken into a drying oven with forced ventilation at 65˚C for 72 hours. At the end of the cycles, two grass stems per reactor were evaluated for each of the following parameters: plant height and number of internodes.
An estimate of chlorophyll content in leaves was determined at the end of each cycle using a chlorophyll meter. First, the chlorophyll meter was calibrated according to the specifications of the equipment, then SPAD readings were measured on one leaf per stem, on two distinct stems per reactor, on the middle third of the stem and the leaf blade.
Nutrients and sodium absorption by the Tifton 85 bermudagrass. The analyses of TN in the aerial part of the plants were quantified by the Kjeldahl semi-micro method with the addition of salicylic acid, adapted from Kiehl [17]. Absortion of K, TP and Na by plants followed the methods suggested by Silva [18].
Based on the analyses results, the nutrients/Na accumulation by Tifton 85 bermudagrass was calculated using Eq 1 [19].
where, A TN, K, TP or Na -nutrient or Na accumulation (g), CNutri-content of nutrients or Na in dry matter of the aerial part of Tifton 85 bermudagrass (g g -1 ); DMY-dry matter of the aerial part of Tifton 85 bermudagrass. Effects of gibberellic acid on Tifton 85 bermudagrass (Cynodon spp.) in constructed wetland systems Effects of gibberellic acid on Tifton 85 bermudagrass (Cynodon spp.) in constructed wetland systems The nutrients/Na uptake by Tifton 85 bermudagrass was calculated using Eq 2 [19].
where, U TN, K, TP or Na -nutrient or Na uptake (g m -3 ), A TN, K, TP or Na −nutrient or Na accumulation (g), Vol-CWs volume (m 3 ). The nutrients/Na uptake rate by Tifton 85 bermudagrass was calculated using Eq 3 [19].

Effluent evaluations.
In the experiments, influents and effluents were evaluated weekly, i.e., every two cycle times (HRT of 3.5 d). Analyses of oxidation-reduction potential, pH, electrical conductivity, evapotranspiration, COD, total nitrogen, total phosphorus, K, Na and total dissolved solids were performed to understand the current reactions and the nutrients and Na removal efficiency in CWs. The parameters analyzed are presented in Table 6.
The Tifton 85 bermudagrass evapotranspiration (ET c ) was determined in each plot and after each cut (0.10 m). ET c was calculated as the difference between the volumetric inflow and outflow divided by the surface area of the plot for data collected during the exchange of SMW, Eq 4 [20].
where, ET c (mm d −1 ), Q in is the plot volumetric inflow (mm 3 d −1 ), Q out , is the measured plot volumetric outflow (mm 3 d −1 ), and SA is the measured plot surface (mm 2 ). The analyses of COD, TP , K, and Na were performed according to the Standard Methods for the Examination Water and Wastewater [21]. TN of effluents was quantified according to the Kjeldahl semi-micro process with the addition of salicylic acid.
To obtain the actual value of the mass reduction in each parameter evaluated, the effluent concentrations of each CWs were corrected by ET c (initial volume applied minus the final collected volume) and were calculated using Eq 5 [19].

Statistical analysis
All numerical data were subjected to variance analysis [23]. When they were significant at the 5% level, it proceeded to regression analysis and it was adopted the model with the highest coefficient of determination (R 2 ) and that expressed the behavior of the phenomenon. The regression coefficients were subjected to t-test at the 5% level of significance. When the regression model presented the level above 5%, the probability found for the model was used. When the regression model did not explain the phenomenon, the Dunnett's test was used to compare the GA 3 treatments with the control.

Yield and vegetative development of Tifton 85 bermudagrass
In the first cycle, the parameters plant height, dry matter yield (DMY) and Tifton 85 bermudagrass SPAD index were significantly influenced by various doses of GA 3 (Table 7). In the second cycle, plant height was the only significant effect for the levels tested. The productivity of Tifton 85 bermudagrass had no significant effect for the second cycle; the same was true for the SPAD index and number of internodes. One of the roles of GA 3 is to promote elongation of the stem [12]. This effect was observed in this study for plant height of Tifton 85 bermudagrass (Figs 3 and 4). The increase in plant height is due to the lengthening of the stem internodes, as there was no difference in the number of internodes (Table 7). Comparing the best treatment for plant height (100 μM GA 3 ) with the control, it is observed an increase of 35.10% and 20.90% in plant height for the first and second cycles, respectively.
In the two cycles observed, the control without the application of GA 3 showed lower growth than the other treatments (Fig 3). The Tifton 85 bermudagrass is well suited for the conditions of CWs and may be able to develop aerenchyma in hypoxic conditions [24]. However, it is not considered an aquatic plant, and, in municipal wastewater treatment systems, this plant is grown in media with high concentrations of salts (K and Na). Thus, one can deduce that the grass is growing in abiotic stress conditions. When plants grow under abiotic stress conditions, they have their levels of phytohormones altered and, therefore respond with reduced growth [25].
When studying the effect of exogenous application of GA 3 on dwarf rice, a significant increase in plant height was observed compared to a treatment without the phytohormone [26].
The increase in height of Tifton 85 bermudagrass in the first cycle was approximately 3.3 x 0.5 , being x the doses of GA 3 , though the increase in the same parameter for the second cycle was 1.3 x 05 . The reduced response to GA 3 may be linked to the age of plants because the GA 3 response in older tissue was less prominent. Gibberellins are produced in young tissues of the shoot system and developing seeds [27]. The effect of GA 3 on plants is influenced by weather conditions [28]; possibly high temperatures in the second cycle negatively affected the induction of Tifton 85 bermudagrass growth.
The DMY of aerial parts of Tifton 85 bermudagrass in the first cycle presented a square root behavior for the various doses of GA 3 (Fig 5). The highest productivity found in the regression model was 4097.70 kg ha -1 depending on the optimal dose of 53 μM GA 3 . It was observed in the control treatment (0 μM GA 3 ) 2682.60 kg ha -1 yield, i.e., 33.80% lower than the yield obtained with the treatment of 50 μM of GA 3 .
The increase of the DMY in aerial parts of Tifton 85 bermudagrass with the exogenous application of increasing doses of GA 3 indicated that the phytohormone plays an important role in the absorption of water and nutrients and processing the accumulation of carbohydrates within plant tissues. The same was observed by Zang during foliar application of GA 3 and evaluation of dry matter in blueberry leaves [29].
The GA 3 stimulates cell division and cell elongation and the exogenous addition of this phytohormone leads to an increase in plant growth, improving the availability of endogenous GA 3 . This view is supported by Kaur [30] and Tuna [31], working with chickpea and corn,  Effects of gibberellic acid on Tifton 85 bermudagrass (Cynodon spp.) in constructed wetland systems respectively. Arabidopsis thaliana plants treated with GA 3 showed a higher shoot development compared to the control treatment [32].
Plants have an important role in the removal of nutrients in CWs, however, the removal capacity can be improved if crop conditions that allow rapid growth and dry matter production are created [33]. In view of this, the application of GA 3 is a viable practice to improve the development of plants, as well as to create conditions to increase the removal of nutrients (pollutants).
As for the SPAD index evaluated in the first cycle (Fig 6), an increasing and linear behavior among the different doses of GA 3 was observed , indicating a lower chlorophyll content in the treatment without exogenous application of GA 3 .
Since the plants were grown in adverse conditions (flooded media and the presence of salts), this lower amount of chlorophyll in the control could be attributed to the formation of proteolytic enzymes such as chlorophyllase, an enzyme responsible for the degradation of chlorophyll [31], which may damage the photosynthetic system. The application of GA 3 improved chlorophyll levels in plants grown in CWs. The absence of effect of GA 3 for DMY and SPAD index in the second cycle may be related to the period of application of the phytohormone. Some studies reported that GA 3 applied during the spring season promotes a high forage growth and that response to the hormone is reduced in summer due to the high temperatures [28]. As can be seen in Fig 2, the maximum and average temperatures in the second crop cycle (near the summer season) were higher than the temperatures in the first cycle. Other factors which may have influenced the absence of effect of GA 3 in the second cycle are the age of the plant, space restriction for root growth, and a lower GA 3 response in older tissues.

Absorption of nutrients and sodium by Tifton 85 bermudagrass
The parameters, total nitrogen and total phosphorus accumulation, and total nitrogen uptake rate by Tifton 85 bermudagrass had a significant effect under the different doses of GA 3 compared to the treatment without the application of GA 3 ( Table 8). The same parameters were evaluated in the second cycle and were not significant at the 5% level. That absence of significance could be correlated with the absence of the GA 3 effect for the productive parameters of the second cycle.
The absence of significant parameters evaluated in the second crop cycle may have occurred due to the room limitation within the reactors for root growth as a function of plant age. Also, the height of the support media was 21.20 cm, while Matos, cultivating Tifton 85 bermudagrass in CWs, reported that they reached a depth of 30 cm [34]. Another factor that may have influenced the absence of significant responses was GA 3 lower response in older tissues since this hormone is produced in young tissues [27].  Where: ns non-significant � and �� , significant at 5 and 1%, respectively, the F test. Total nitrogen accumulation (g) (A TN ), potassium accumulation (A K ), total phosphorus accumulation (A TP ) and sodium accumulation (A Na ). Total nitrogen uptake rate (g m -3 d -1 ) (UR TN ), potassium uptake rate (UR K ), total phosphorus uptake rate (UR TP ) and sodium uptake rate (UR Na ). https://doi.org/10.1371/journal.pone.0206378.t008 Effects of gibberellic acid on Tifton 85 bermudagrass (Cynodon spp.) in constructed wetland systems It is worth mentioning that the experiment aimed to verify the effect of the GA 3 application in the same cutting stage of the control without application, but one of the suitable practices of forage cultivation is to establish the cut by height and not by the time of cultivation, so it would be possible to perform more cuts when applied GA 3 which would result in greater removals of the pollutants from the CWs.
In Fig 7, the effects of the A TN and UR TN parameters were adjusted to the square root regression model with p = 0.1411 and p = 0.1417 respectively, making it possible to understand the relationship between GA 3 doses and TN removal by plants.
For the generated model (Fig 7) the highest accumulation and TN uptake rate was 0.56 g and 3.92 g m -3 d -1 respectively, for the 57 μM dose of GA 3 .
In the first cycle, GA 3 promoted an increase in the DMY of the aerial part of the plants. Associated to this effect, the plants accumulated larger amounts of N resulting in higher rates of removal of this nutrient (Fig 7). Mustard plants accumulated higher amounts of N in the shoots and seeds in treatments with the exogenous application of GA 3 [35].
The application of GA 3 in environmental conditions with a high concentration of atmospheric CO 2 promoted a lower concentration of N in Arabidopsis thaliana. However, for environments with different concentrations of CO 2 , the DMY was higher for the GA 3 treatments compared to the control. As a result, plants treated with GA 3 had a higher accumulation of nitrogen [32].
The accumulation of TP in the aerial part of Tifton 85 bermudagrass is presented in Table 9. Because it did not present a significant regression model to explain the effect of the doses of GA 3 , the treatments with the control were compared by the Dunnett's test, where only level 1 (5 μM GA 3 ) differed from the control level without GA 3 application.
The treatment with exogenous application of GA 3 in a maize crop under salinity conditions increased foliar TP [31]. The higher growth and productivity of Tifton 85 bermudagrass treated with GA 3 compared to the control promoted a greater accumulation of TP for the N1 treatment. From the environmental point of view, it is an advantage since it showed that the removal of TP in CWs can be improved with an application of GA 3 .

The efficiency of constructed wetland systems in the removal of pollutants from synthetic municipal wastewater
The influent had an average oxidation-reduction potential (ORP) of 127.2 mV. Only in three peaks (Fig 8), was it possible to observe ORP above the input value, indicating a small effect of the oxygenation promoted by the root system. This may be an indicator of the formation of aerenchyma by the roots of Tifton 85 bermudagrass. This formation in the plant tissue can be characterized as a mechanism of adaptation in flooded systems. Other studies have described and presented histological sections with the presence of aerenchyma in roots of Tifton 85 bermudagrass in flooded systems [24,36].
ORP evaluated throughout the experimental period remained between 0-204.75 mV. A critical factor for nitrification would be an ORP above 500 mV [37].
In function of the ORP evaluated in the present study, it was possible to indicate that the conditions of the CWs were anoxic/facultative. A classification for the ORP of the medium, ORP above 400 mV has the predominance of O 2 (oxidized medium); between 400-200 mV, the presence of O 2 and nitrate (NO 3 -); and below 200 mV indicates low oxygen presence tending to anoxic/facultative conditions [38].
After collecting the effluents, the pH and EC were determined, and their values are presented in Fig 9. The pH of the inlet was on average 6.7 and the EC of 336.1 μS cm -1 . Except for the NC � , the treatments presented a similar behavior for the parameter pH and EC. In both cycles the pH and EC were higher for the NC � , this allows to understand the importance of the removal of the pollutants (nutrients and Na) and other salts by the plants.
The increase in pH in the different treatments may be related to the use of carbon dioxide by plants and algae [5].
The pH values ranged from 6.5 to 8.2 for treatments and control plots, what indicated that SMW (influent and effluent) presented ideal conditions for the survival of microorganisms and for the degradation of organic matter [39]. The optimum range for nutrient absorption is between 5.5 and 6.5, pH above 7.0 causes a great restriction in the availability of micronutrients and P [40].
The daily average evaporation and ET c was calculated through the determination of the volume of the input (influent) and the output (effluent) (Fig 10). Subsequently, the mass differences were calculated to express the removal of nutrients and Na. The system without Tifton 85 bermudagrass showed an average evaporation of 2.0 mm d -1 and the system with Tifton 85 bermudagrass grass showed an average ET c of the treatments of 5.2 mm d -1 .
It is possible to observe the importance of plant cultivation in CWs since they transform the SMW into water into the atmosphere through transpiration. Thus, in addition to removing the pollutants (nutrients), CWs with plants reduce the volume of sewage for the final disposal.
The behavior of TN and TP concentrations in the effluent is shown in Fig 11. Initially, the removal of TN and TP by plants compared to the control without plants was perceptible. The removal of these elements in CWs reduces the N and P disposal into rivers and lakes, these nutrients are essential for the growth and multiplication of algae that are responsible for the eutrophication of the water bodies. Therefore, plants play an important role in the uptake of N in the forms of NO 3 and ammonium (NH 4 + ); the ammonium ion can also be adsorbed by plants roots. The sum of these N forms is described as TN [41]. The presence of plants in CWs allows a greater removal of N [42,43].
In the first cycle, the COD removal was higher for the CWs without plants (NC � ), and the CWs with the treatment levels N3 and N4 were the ones with the lowest COD removals compared to the two controls (Table 10).
The treatments N3 and N4 were those that obtained higher plant height and DMY in comparison with NC (1st cycle) ( Table 10). As a result, it was possible to deduce that these treatments presented a greater growth of the roots. The application of GA 3 can stimulate root growth [12]. Together with the greater development of the aerial part, these treatments may have released more exudates through the root system, resulting in higher concentrations of COD in the effluents. Exudates are secretions that release ions, free oxygen, water, enzymes, Effects of gibberellic acid on Tifton 85 bermudagrass (Cynodon spp.) in constructed wetland systems mucilage and a diversity of primary and secondary metabolites with carbon in their composition [44].
In the second cycle, the NC � treatment presented a lower COD concentration when compared to N1 and N4. The treatments that received gibberellin (5, 50 and 100 μM GA 3 ) did not differentiate from the NC control. This lower COD efficiency removal in treatments with the presence of Tifton 85 bermudagrass may be related to the release of root exudates and plant senescence.
For the nutrients, Na, and TDS all treatments were superior to NC � , demonstrating the importance of plant species in CWs.
The treatment with 50 μM of GA 3 showed a greater removal of TP and Na from the SMW compared to NC and NC � controls by the Dunnet's test. This fact is interesting because P is one of the main nutrients responsible for eutrophication of rivers and lakes and this nutrient must be removed to comply with the legislation of disposal of effluents in water bodies. Although the efficiency of removal of TP is greater than 90% in N3, it does not comply with the legislation of Minas Gerais [45] that ranges from 0.020 to 0.025 mg L -1 , being necessary to associate with another treatment system or recirculate therein system. The total efficiency of CWs was obtained through the system-substrate-plant-micro-organism interaction and solar radiation. In this study, the reason for the high efficiency of nutrient removal by plants was due to the low load applied, even though this load was sufficiently great to verify the effect of GA 3 on nutrient removal by plants, thus meeting the initial objective of the work.
The efficiency of municipal wastewater treatment of subsurface vertical flow constructed wetland cultivated with vetiver grass, observed a removal of TN of 93.9% and for TP of 90.5% for application rates of 0.8 and 0.3 g m -2 d -1 , respectively [46]. These efficiencies and application rates were similar to the ones in the present study.
The TN uptake rate by the Tifton 85 bermudagrass treated with 50 μM of GA 3 was 0.69 g m -2 d -1 of TN, as a function of an application rate of 1.02 g m -2 d -1 of TN through SMW (Table 11). The TN removal in the present study was higher than the values observed by Matos (2008), Fia (2011) and   [9,34,47]. These authors found removals of 0.46, 0.42 and 0.40 g m -2 d -1 TN from the aerial parts of Tifton 85 bermudagrass, respectively. Only Fia (2011) worked with a high rate of TN application, around 25.70 g m -2 d -1 TN, via swine wastewater [47].
Jesus (2016) observed removals of 0.79 g m -2 d -1 when treating municipal wastewater with TN application rates equal to 9.20 g m -2 d -1 in CWs with Tifton 85 bermudagrass. These values of application and removal rates were higher than the ones in the present study [48].
For different application rates and different influents (Table 12), the TP uptake rate of 0.11 g m -2 d -1 by the Tifton 85 bermudagrass treated with 50 μM of GA 3 was higher than that For each cycle, treatment averages followed by letters (a) and/or (b) did not differ from NC and/or NC � controls, respectively by the Dunnett's test at the 5% level of significance. COD-chemical oxygen demand; TN-total nitrogen; TP-total phosphorus; K-potassium; Na-sodium; TDS-total dissolved solids. NC-0, N1-5, N2-25, N3-50 and N4-100 μM GA 3 ; NC � -control without plants in CWs.

References
Although the use of GA 3 showed potential for nutrient removal via biomass of Tifton 85 bermudagrass cultivated in CWs in the treatment of SMW, it would be interesting to conduct further field testes with higher load applications.
An application of GA 3 in spring and another in autumn resulted in the increased dry mass of azevem and white clover [49]. Thus, responses to GA 3 strongly depend on the time of year in which the plants were treated. And for perennial crops such as Tifton 85 bermudagrass, GA 3 may be interacting with a fixed strategy for seasonal plant growth [50].

Conclusions
The gibberellic acid promoted higher growth in height of Tifton 85 bermudagrass in constructed wetland systems.
The application of 50 μM of GA 3 provided higher dry matter yield and nitrogen removal by Tifton 85 bermudagrass in constructed wetland systems for the first crop cycle.
The application of 50 μM GA 3 was the best dose for the global removal of P and Na from synthetic municipal wastewater in the first crop cycle of Tifton 85 bermudagrass in constructed wetland systems.
The gibberellic acid applied in the second crop cycle had no effect on dry matter yield and nutrient removal by Tifton 85 bermudagrass grass in constructed wetland systems.
Supporting information S1