The trial was established as a random block design and consisted of six different treatments (Table 1), each of which was performed in triplicate. Pots of 4.5 L were used to ensure effective control over the variables. Each pot was filled with 3 Kg of the following mixture: 60% soil; 30% sand and 10% composting cow dung. The soil was collected from ESAC’s agricultural land and had an approximate pH of 6.4 with an electrical conductivity close to 0.37 ms/cm. Chemically, it had the following constitution: P — 44 mg/Kg; K — 75 mg/Kg; calcium oxide (CaO) — 140 mg/Kg; magnesium oxide (MgO) — 35 mg/Kg. In total, 18 pots (3 per treatment) and 36 plants (2 per pot) were used.

A drip irrigation system was provided. At an early stage of the trial, when the plant’s roots did not occupy a large volume, the plants were watered four times a week using 400 mL per pot. The amount of water was adjusted to the needs of the plants according to temperature in the greenhouse, in order to compensate any losses by evapotranspiration.

The sowing of the bean crop took place on 24th of February 2020. Two seeds were placed on the soil surface per pot, spaced 10 cm apart. The gap between pots was 50 cm, the density of plants was 200 000 per hectare, and was determined using the methodology proposed by Lee & Herbek [29]. After sowing, at the beginning of each stage development, the plants were monitored according to the Gepts & Fernandez scale [30].

The Calmar fertilizer® was applied in the trial 0 DAS at the dosage of 100 kg/ha. To that end, 2 g of product were applied in the soil surface of each pot in T1, T4 and T5 treatments.

Approximately seven weeks later, foliar treatments with Profertil® and Albit® started. Both products were diluted following the supplier’s instructions, i.e., 5 mL of Profertil® to one litre of water and 0.2 mL of Albit® to another litre of water. Afterwards, using 500 mL sprayers, the Profertil® diluted solution was applied in T2 and T5 and the Albit® diluted solution in T3 and T4. In total, 3 foliar treatments were applied between 9 and 10 AM; the first was applied at 47 DAS, when the plants reached their reproductive phase (pre-flowering). The second was applied during the pod development stage, at 61 DAS. The final application was at 76 DAS during the pod-filling stage. Approximately six days later, the plants reached their final stage of maturity with almost no leaves, consequently, the foliar treatments ended.

During the trial, the crop was maintained by removing weeds from the pots by hand, so these did not compete with the beans.

Evapotranspiration depends on the crop coefficient (Kc), in turn, this depends on the vegetative and reproductive phases of the bean life cycle. Usually, Kc hits its maximum value when a crop reaches the flowering stage. During this stage, approximately 49 DAS, the watering plan was changed to twice a day. Every day, approximately 300 mL of water were administered per pot by the irrigation system, at 9 AM and 5 PM, totalling a daily volume of 600 mL. The implementation of this plan significantly reduced the gap between irrigations, and because the pots did not have a saucer, some essential nutrients may have been lost to leaching. In order to restore some of these, at 64 DAS, the NPK ENTEC© 13-10-20 (Eurochem – DEIBA) complex granular fertilizer was applied in the soil surface of every pot.

A CCM-200 (Opti-Sciences) chlorophyll content meter was used to quantify the CCI of the leaves. Measurements were made at 70 DAS for all treatments during the pod-filling, reproductive stage. One plant was randomly selected per treatment, of which, a fully expanded leaf was chosen to sample. Then, three readings (n = 3) were performed in distinct regions of the leaf: central, proximal (leaf base) and distal (leaf apex).

On the 4th of June 2020, 101 DAS, the greenhouse trial ended, and all the pods were harvested and grouped according to their treatment and left to dry at room temperature. Three days later, the pods were counted by plant/treatment. The observations of the pod morphology consisted in measuring pod the length, width, and the ratio between both. A calliper was used to measure the length and width of the pods according to Ningsih et al. [31] methodology. The beans were removed from the pods, counted, and placed according to their treatment in small aluminium packages. Later, the bean mass in each treatment was recorded using an analytical balance (PB3002-S Mettler Toledo).

Before the dry bean laboratory analyses, 20 beans from each treatment were randomly selected and placed in separate vials. Those beans were ground according to the Association of Official Agricultural Chemists (AOAC) 922.02 method [32] by using a coffee mill (KG39 – DeLonghi) and a 0.5 mm sieve, later, the ground bean samples were placed back in their correspondent vials. In order to analyse beans nutritional value and mineral characterization, the ground samples needed to be dried. During this process the moisture content corresponding to the samples of each treatment was determined in triplicate (n = 3), according to the method 930.04 of AOAC [32], by drying 2 g of ground sample at 105 $\pm$ 5 °C over a period of 2 h in an oven (Memmert, Büchenbach, Germany). After getting moisture analyses values and the mass of dry beans for each treatment, their yield was obtained using the method described by Lee & Herbek [29] extrapolated to kg/ha. For the cooked bean analyses, a sample of 15 beans of the two treatments with higher yield (kg/ha) were randomly selected and soaked in water for 12 h, then they were cooked on a pressure cooker during approximately 18 min. Afterwards, using a pestle and mortar, the beans were partially crushed. The method applied for the cooked beans moisture analysis was the same applied for ground bean [32].

The nutritionally analyses of the dry bean, consisted in the determination of protein, ash, fat, and fibre contents of the samples previously prepared corresponding to the different treatments. All these analyses were performed in triplicate (n = 3) according to the methods described by the AOAC [32]. The crude protein was determined by a macro-Kjeldhal method (N $\times$ 6.25) by digesting 0.5 g of sample in a digestor (VELP Scientifica, Usmate Velate MB, Italy) and distilled by a Kjeldhal tube distiller (VELP Scientifica, Usmate Velate MB, Italy). The ash content of the dry bean was obtained by incinerating 2 g of sample at 550 $\pm$ 15 °C in a muffle furnace (INDUZIR, Portugal) for 2 hours. The crude fat was determined by the extraction of 2 g of sample with diethyl ether using a Soxhlet device (Behr Labor Technik, Dusseldorf, Germany). The crude fibre was determined by using the fraction of the dry bean remaining after digestion of 2 g of sample with 12.5 g/L sulfuric acid using a fibre extractor (Labconco raw fibre extractor, United States of America (USA)) followed by vacuum filtration (General Electric, Boston, MA, USA) and digestion of the same sample with sodium 12.5 g/L accompanied by a new vacuum filtration, using the same equipments. The total carbohydrates (g/100 g) were calculated according to Saupi et al. [33]: 100 – (% fat +% protein +% fibre +% ash). The energetic value (Kcal/100 g) of the dried beans was calculated using the Atwater factors: 9 Kcal/g of fat, 4 kcal/g of protein and 4 Kcal/g of carbohydrate, as described by Giuntini et al. [34]. The methods applied for the analysis of ash and protein contents of dry bean were the same applied for cooked bean [32], however, in this case, the analyses were only performed in duplicate (n = 2).

The dry bean mineral characterization reads, and measurements were performed in triplicate (n = 3). In each analysis, 0.5 g of ground sample were taken from the vials corresponding to the different treatments and submitted to the drying process at 105 $\pm$ 5 °C, later, the samples were incinerated at 550 $\pm$ 15 °C on a muffle furnace (INDUZIR, Portugal) for 14 h. After cooling, the samples in the crucibles were digested with hydrochloric acid 20% (m/v). Afterwards, crucibles were placed in a water bath (Memmert, Büchenbach, Germany) for 30 min at 100 °C. The samples were filtered to 50 mL volumetric flasks and the remaining volume was filled with distilled water, in order to obtain a saturated solution.

Quantification of Zn, Fe, Mn, and Cu was made directly from the saturated solution by flame atomic absorption spectrometry (PerkinElmer PinAAcle 900 T, Waltham, MA, USA) adopting the methodology of Lucas & Sequeira [35]. To quantify Ca, Mg and K, the solution was diluted (1:10) and then 2.5 mL of 0.75% strontium chloride added to each sample, this reading was performed using the same methodology [35].

For the determination of P, 5 mL of saturated solution and 10 mL of ammonium molybdate-vanadate were pipetted in to 50 mL volumetric flasks. The remaining volume was filled with distilled water. The measurements of P were made by spectrophotometry (Pye-Unicam SP6-350 visible spectrophotometer) adopting the methodology of Ribas et al. [36].

The Shapiro-Wilk test was performed for the normal distribution of the dry bean analyses obtained data. Later, that data were treated using a one-factor analysis of variance (one-way ANOVA) to evaluate the effect of treatments on the parameters studied. The significance of differences between evaluated mean values was determined using the Tukey test, at a significance level of p 0.05. The software used was the IBM SPSS® Version 25 for Windows (IBM, Chicago, USA). In the analysis of the cooked beans, as the obtained data only resulted from two observations per treatment (n = 2), the mean was the only statistical method performed.

By monitoring the day that the plants reached each phase and stage of development, it was possible to calculate their length in days. Those periods are presented in Table 2.

The total length of the crop life cycle was 101 DAS, obtained by adding the number of days involved in the different development stages (Table 2). This value fits into the interval of 65 to 120 days corresponding to the duration of the Phaseolus vulgaris L. life cycle [37]. Between 2005 and 2006 in the Minas Gerais State - Brazil, Melo et al. [38] made several open field trials using the BRSMG Realce variety and obtained a life cycle of around 83 DAS, 18 days shorter than the present crop that was establish in greenhouse. Amalfitano et al. [39] studied broad bean (Vicia faba L.) on two farming systems (open field, greenhouse), the results revealed a higher harvest precocity for the crops established in greenhouse compared to open field ones. This fact could indicate that a protected environment can provide a quicker life cycle to the plants, however, that did not happen in the present study when comparing with the trials made by Melo et al. [38]. As claimed by Oliveira et al. [40] the bean life cycle length may vary according to the edapho-climatic conditions of the producing region. That said, according to the Köppen-Geiger classification, there are clear differences between the climate where the current trial took place and the tropical savanna climate type predominant in the region where the Melo et al. [38] trials were carried out. Also, the present culture may have been influenced by the number of sunshine hours, as stated by Song & Jin [41], variations in the number of sunshine hours can affect crops growth process.

By analysing the results obtained from CCI readings (Table 3) in the different treatments, it was observed that T2, T4 and T5 differed significantly from the control. Overall, T2 produced the highest mean value of CCI, i.e., about 34% higher than the control, followed by treatments T4 and T5 respectively, where the CCI values increased about 30% as compared to control. There were no significant differences between T3 and T1, compared to control.

The increase of CCI in T4 may indicate a possible positive synergy between the two treatments used, even though when used individually, they did not differ from the control. The increase of CCI found in T5 may be related to the effect of Profertil® and might indicate that this product had a greater impact on the increase and maintenance of chlorophyll in the plant tissue, when compared to Calmar®.

There is a strong relationship between leaf chlorophyll content and nitrogen content, since approximately 70% of the plant total nitrogen is trapped in chloroplast enzymes that participate in the degradation of the molecular structure of chlorophyll [42]. Chlorophyll-measuring devices such as the CCM-200 have been described several times as potentially suitable for assessing total nitrogen in the leaves of many crops [43]. That said, higher values of CCI may be associated with higher levels of N, which could result in greater photosynthetic activity.

In other studies where bio-inputs were used on several crops, similar results were found to the present study. Ali et al. [18] and Goñi et al. [44] using Ascophyllum nodosum formulations in a tomato crop (Solanum lycopersicum L.) obtained a higher CCI in the treated leaves as compared to control. The single use of Phymatolithon calcareum extracts in tomato (S. lycopersicum) [45] and mango crops (Mangifera indica L.) [46] did not promote significant changes of CCI on the treated leaves when compared to the control.

Conversely, Gins et al. [47] concluded contrary to what occurred in the present study, the CCI recorded in the leaves of Passiflora (Passiflora incarnata L.) treated with Albit® were significantly different from their control.

The values corresponding to pod length (PL), width (PW), ratio between length and width (L:W) and number of beans per pod (NBP) are presented in Table 4. In a primary analysis of the PL data, it was verified that the highest values were obtained in treatments T5, T4 and T3, these differed significantly from the control, but not between them. There were no significant differences observed between treatments on the PW parameter. The lowest L:W ratio was obtained on T0 which differed statistically from all of the treatments. Regarding the NBP parameter, once again, the control treatment produced the lowest values, however, there were no significant differences between any treatments.

It is notable that the application of bio-inputs in this crop may have contributed to the increase of PL and consequently L:W parameters, compared to control. However, no significant differences were found between T1, T2, T3, T4 and T5. The correlation between NGP and the L:W ratio could be important for estimating bean yield. Although, the Pearson correlation (r = 0.183) between both variables was low according to the scale of Cohen et al. [48]. The application of bio-inputs may have induced changes in the L:W ratio, but ultimately did not influence the NBP, since no significant differences were found between control treatment and the others.

In studies carried out by other authors, significant differences were also found in fruit dimensions of treated plants with bio-inputs. Eris et al. [49]reported significant increases in fruit length on pepper plants (Capsicum annuum L.) treated with an extract of A. nodosum. Likewise, Ali [50] used an Ascophyllum-based product and obtained significant increases in the length of fava pods (Vicia faba L.). Aguiar [51] recorded a significant increase in the diameter of grape berries (Vitis vinifera L.) using an extract of P. calcareum-based product. Similarly, Amatussi et al. [20] also achieved a significant increase in the diameter of tomato fruits (Solanum lycopersicum L.) with a P. calcareum extract. Morais et al. [52] obtained significant increases in the length of strawberry fruit (Fragaria x ananassa Duch.) in plants inoculated with PGPB (Pedobacter sp. CC1 and Bacillus safensis).

From the results of the bean plants productivity (Table 5), it was observed that T1, T3 and T4 produced the highest values for the NP (number of pods), with an average increase of 34%, as compared to control. As for the NB (number of beans), T1 stands out with a value approximately 47% higher than T0, differing statistically from all other treatments. The results obtained for the NP and NB parameters revealed a strong correlation (r = 0.921) between them, according to Cohen et al. [48]. For BM (bean mass) (g), the highest values were recorded in T4 and T1, showing an average increase of around 42% compared to control. It was also possible to verify the existence of a high correlation (r = 0.787) between NB and BM [48]. Given the existence of these correlations, it was expected that treatments with the highest values of NP, NB and BM would correspond to the greatest yields in kg/ha. That was verified in treatments T4 and T1 that stood out, with a yield 43% higher than control. In general, the treatments subjected to the application of bio-inputs showed significant increases in all parameters.

The application of Profertil® may have inhibited Calmar®, since a higher yield was obtained from the solo application of the latter in T1, when compared to their combination in T5. The clear difference in yields between T1 and T2 (Table 5) was indicative of inhibition or interference between the two products.

In the cultivation of the ground cherry (Physalis peruviana L)., Oliveira [53] achieved a higher fruit yield with P. calcareum treatments. Russo et al. [54] applied Albit® to P. vulgaris leaves and reported significant increases in pod productivity as compared to both control and treatment with proline. Tandon & Dubey [17] obtained a higher number of pods and seed yield per plant in soy crop (Glycine max L.) when applying a A. nodosum commercial extract.

The results obtained from dry bean nutritional analyses are presented in Table 6. Regarding the ash content, amongst all of the treatments, T2 and T5 presented the highest values. However, these did not differ significantly from their control. The highest values of fat content were obtained in T5, but these were not statistically significant. In fibre analyses, T1 samples reported significantly higher values than control. The protein content levels obtained for T4 and T5 were highest. Analyses of carbohydrates analysis for T1, T2 and T3 produced the highest values, which were significantly differently. In general, the calorific data were very uniform. In every treatment, there was a decrease in some parameter as compared to control, as was the case of T1 for ash, fat and protein content, T2 for protein content, T3 for ash, fat and protein content, T4 for ash, fibre and carbohydrate content and T5 for carbohydrate content.

A correlation can be established between the amount of protein and total nitrogen present in dry beans samples, since in most foods non-protein nitrogen is represented by a small fraction that is not significant [55]. The common bean establishes symbiotic relations with nitrogen-fixing bacteria, e.g., Rhizobium spp. and Bradyrhizobium spp. [56]. When successful, this symbiosis allows these bacteria to invade the roots of legumes, forming nodules, stimulating the growth and development of plants, partly due to the absorption of nitrogen and other macro- and micro-nutrients [57]. There are several studies reporting that these metabolic changes in legumes lead to a nutritional increase of their seeds, mainly in protein content [58, 59]. In the present study, by analyzing more thoroughly the protein levels recorded for each treatment, it was demonstrated that for T1, T2 and T3, the protein values were significantly different and lower than the control. This could indicate that the application of bio-inputs impaired, in some way, the biological nitrogen fixation and consequently the symbiotic relation between plants and rhizobium, reducing the number of nodules in the rhizosphere. On the other hand, significantly higher protein content values were obtained in T4 and T5, when compared to T0 which may indicate that the bio-inputs applied favoured the biological nitrogen fixation reaction. Zhai [57] and Khan et al. [60] using A. nodosum soluble extract in alfalfa (Medicago sativa L.) crop reported an increase in the number of nodules in the treated plants compared to the control, suggesting that those extracts can enhance the symbiotic relationship between the plant root and nitrogen-fixing bacteria.

The literature about possible nutritional changes caused by the bio-inputs application in agricultural crops is scarce, however, there are some publications where results relate with the present study. Kocira et al. [61] used a commercial product derived from E. maxima on P. vulgaris, obtaining a higher protein content in treated plants seeds compared to control, as it occurred in the present study in T5 and in opposition to what occurred in T3. In Eryashev et al. [62], the use of Albit® in pea (Pisum sativum L.) promoted significant increases in protein content compared to control, as it occurred in T4 and contrary to T3 of the present study.

The results of the dry bean are presented in Table 7. Regarding the P content, all treatments differed significantly from the control except T1. The highest content was found in T5, i.e., about 16% higher than the control. T3 values were overall the lowest. As for K, a considerably lower value was obtained for T5, which was the only one to differ statistically from T0. The analysis of Ca, Mg and Cu indicated that the levels achieved by the several treatments did not differ significantly from each other, or from the control. For Zn content, the highest value was obtained in the control followed closely by T5. The remaining treatments were considerably lower, although they did not differ significantly from each other, they differed from T0 and T5. Regarding Fe amount, treatments subject to the application of bio-inputs did not differ statistically from each other, however, they were much lower than T0. For Mn, none of the treatments were significantly different from T0.

Taken together, the results indicated that none of the treatments where bio-inputs were applied contributed to mineral enrichment of the dry bean, since only T2, T4 and T5 recorded significantly higher values than control for P content. The most evident difference for these data was associated with the Fe content in the control. Comparing this value with the reference content for the most consumed bean varieties in Portugal, as described in PortFIR website [63], it is noted that this amount of Fe is above the normal composition of the dried bean, since the usual values for the varieties of cowpea, red and white, are given as approximately 52, 64, and 61 mg/kg, respectively.

An excess of this micronutrient may be associated with stress caused by Fe deficiency in control plants. Faced with an Fe deficiency, dicots stimulate responses aimed to improve the Fe accumulation from the soil. Fe deficiency can induce a 5–10-fold increases of Ferric-chelate reductase activity [64]. The function of this enzyme is to reduce Fe (3${}^{+}$) to Fe (2${}^{+}$), so that it can be transported by iron-regulated transporters into the cells [65]. In the study of Cohen et al. [64], it was described that Fe deficiency induced the expression of the iron carrier gene, which could facilitate the transport of divalent cations such as Cd (2${}^{+}$) and Zn (2${}^{+}$), besides Fe (2${}^{+}$). This may also be an explanation for the higher levels of Zn obtained in T0, as compared to other treatments. Another hypothesis for this imbalance could reside in the fact that these treatments prevent the uptake of iron, however, there is a large number of studies in which the use of bio-inputs in agricultural crops promoted a greater iron assimilation and other micronutrients, compared with control [66, 67, 68].

Some studies with similar analyses were found, i.e., Amatussi [45] that used a commercial product based on P. calcareum in a S. lycopersicum crop, found that the micronutrient content did not differ significantly from the relevant control. Di Stasio et al. [69] used two commercial extracts of Ascophyllum nodosum in S. lycopersicum, achieving increased P levels as compared to control, agreeing with the observations in this trial in T2 and T5. However, these authors also obtained increases in levels of K, Ca, Mg, Zn and Fe. Petropoulos et al. [70] used a commercial product containing the Bacillus subtilis bacteria and A. nodosum extracts in a greenhouse study of P. vulgaris reporting increases in Cu and Ca content in dry beans as compared to control and other treatments. Their report did not concur with the present trial, since there was no statistical difference between any of the treatments for both elements.

From the cooked bean analysis (Table 8), it is noticeable that the water content was higher in T0, followed by T4 and T1. Samples belonging to T1 and T4 produced approximately 47 and 37% more ash than the control. Regarding the protein content, T4 and T1 were approximately 17 and 12% higher than control. These results indicate that the application of bio-inputs in this crop can reduce cooked bean water content and increase its ash and protein value.

Nutritional values of the T0, T1 and T4 cooked bean were compared with PortFIR website [63] reference values for some commonly consumed bean varieties in Portugal (Table 9).

It was observed that the T0 water content was lower than the three varieties presented above, its protein content was below the reference value for the cowpea and above the value for the butter and white beans, respectively. The ash content in the T0 was lower than the three bean varieties. The T1 water content was lower than the three bean types in Table 8, however, its protein content was higher. Regarding the ash content, T1 was only superior to the variety of cowpea. Doing the same comparison, the T4 treatment behaved in the same way as T1. By these results, it can be stated that the BRSMG Realce variety has the potential to be included in the Portuguese diet, since after the cooking process, it has similar nutritional values to the usually consumed bean varieties.