Extending beyond traditional forage: potential nutritional benefits of native plants in extreme arid insular regions

Abstract

The scarcity, unstable nutritional quality and environmental cost of imported forages in arid insular regions like Fuerteventura in the Canary Islands (Spain) need exploring sustainable local alternatives. This study evaluated the nutritional quality of twelve native and endemic plant species categorized into legumes, grasses, and a mixed group, cultivated under controlled conditions. The bromatological profiles, focusing on fiber, protein, lipids, and minerals, showed significant differences among plants in key parameters of forage quality: neutral detergent fiber (NDF; 24.2–71.3%), acid detergent fiber (ADF; 9.0–40.5%), acid detergent lignin (ADL; 2.0–15.8%), crude protein (CP; 6.1–20.9%), total lipids (TL; 1.5–6.3%), ash content (25.4–88.6%), enzymatic organic matter digestibility (EOMD; 5.9–10.9 MJ/kg), metabolizable energy (ME; 5.9–10.9 MJ/kg), and relative feed value (RFV; 74.8–317.9). Among lipids, all species had a high proportion of polyunsaturated fatty acids (PUFA; 34.7–63.1% of total fatty acids), mainly α-linolenic acid (ALA; 18:3 n-3; 24.8–54.4%) and linoleic acid (LA; 18:2 n-6; 6.4–25.0%). Other beneficial lipid molecules for animal health such as γ-linolenic acid (GLA; 18:3 n-6), stearidonic acid (SDA; 18:4n-3) and phytosterols (PTS) were detected in specific species. Mineral composition analysis revealed that only Ca, Na, Fe and Cu levels were near or above the established maximum tolerable levels (MTLs) in some species. According to literature, most of the species had a similar or slightly lower nutritional value compared to conventional forages such as alfalfa. Thus, the evaluated native species pool could serve as alternative feed for ruminants during forage shortages, suggesting their combined use to improve livestock health and product quality. This research emphasizes the untapped potential of native plant biodiversity to enhance sustainable agro-livestock practices in arid regions, supporting livestock nutrition and conserving unique botanical heritage.

1 Introduction

Forage crops play a crucial role in supporting livestock by significantly contributing to the health, productivity, and yield of animals (Chand et al., 2022). These crops have primarily comprised herbaceous legumes such as alfalfa (Medicago sativa) and various grass species including barley (Hordeum vulgare) and oats (Avena sativa). However, identifying alternative forages from native flora thriving in numerous arid and semi-arid regions where feed resources are scarce is receiving increasing attention in recent decades in response to the growing global food demand (e.g., Díaz et al., 2018; Louhaichi et al., 2021). These species usually exhibit greater adaptability than conventional crops to adverse conditions such as water stress and salinity, making them a valuable tool for adapting to projected climate scenarios of increased aridity while potentially meeting the nutritional needs of livestock over extended periods (Álvarez et al., 2007). Additionally, local forages can provide specific nutritional elements that positively impact the health and well-being of the herd, enhancing the value of derived meat and dairy products (El Shaer, 2010; Elgersma et al., 2013; Capstaff and Miller, 2018; Davis et al., 2022). Obtaining new forage crops requires an exhaustive study of the ecology and physiology of these plants, particularly their nutritional value, which can offer valuable insights into the most suitable candidates for domestication (Norman and Masters, 2023).

Considering that forages, whether fresh or preserved, are an essential part of the animal diet globally, providing high nutritional quality species is critical to ensure the health and productivity of livestock (Capstaff and Miller, 2018). Key parameters like crude protein (CP), digestibility, metabolizable energy (ME), fiber content, and lipid profile largely reflect the nutritional status of a forage (Clapham et al., 2005; Lee, 2018; Chand et al., 2022). These quality parameters vary significantly among different crops depending on their biochemical composition and prevailing agroclimatic conditions (Capstaff and Miller, 2018; Chand et al., 2022). Although there are notable differences in the nutritional value and palatability among forage plants from various botanical families, global data are sparse, and information is even more limited for native species (Lee, 2018; Louhaichi et al., 2021). In this context, a recent review indicates that legumes typically offer higher protein content than grasses, which may be more digestible, while information on other plant groups is scarce (Lee, 2018; Chand et al., 2022).

Forages are also a significant source of fatty acids (FAs) in ruminant diets, and its concentrations in different species are crucial for the quality of dairy and meat products (Elgersma et al., 2013; Davis et al., 2022). For example, a study evaluated the effects of dietary supplementation with fresh Camelina sativa in goats, finding that the milk contained higher levels of polyunsaturated fatty acids (PUFA) and conjugated linolenic acid (CLA) than traditional forages, and this milk also produced cheese with better sensory properties (Colonna et al., 2021). Therefore, increasing the available knowledge about forage compositional characteristics would aid in a more comprehensive understanding of its nutritional value for ruminants (Castro-Montoya and Dickhoefer, 2020; Norman and Masters, 2023).

In the Canary Islands, as in other arid and semi-arid regions of the world, there is a deficit in food production for livestock, becoming the most dependent Spanish region on imports for this sector (PFORCA, 2014; ISTAC, 2022). Particularly, in Fuerteventura Island, one of the more arid territories, nearly 30,000 tons of food are imported annually to feed around 84,000 goats and sheep (ISTAC, 2022), significantly impacting the environment due to the carbon footprint and greenhouse gas emissions from mainland transport. In contrast, this archipelago is an oceanic hotspot harboring an endemic plant diversity of more than 600 taxa, representing over 50% of the total native flora (Caujapé-Castells et al., 2022). Fuerteventura alone accounts for approximately 15% of the Canary Islands’ total endemic plant species (Biota, 2024). This taxonomic and genetic diversity provides numerous options with considerable potential as livestock forage. However, information about the performance and quality of these plants is limited and restricted to a few varieties (Méndez et al., 2003; Chinea et al., 2014; Álvarez et al., 2007; 2018). More specifically, to date, quantitative data on the lipid composition of native and endemic plant species are almost null (Guil-Guerrero et al., 2000a, b). The provision of new functional forage resources would allow increasing the pool of plant raw materials for rationing and opens the possibility of recovering abandoned farmland, contributing to environmental sustainability, soil conservation and the fight against desertification (Chinea et al., 2014; Hunter et al., 2017; Díaz et al., 2018). Given this potential, we pose the following questions: Will this biodiversity reservoir be able to meet the nutritional requirements of livestock? Will these species provide nutritional compounds that are potentially beneficial to animal and human health, thereby adding value to derived products such as milk and cheese?

Considering the above factors, this study aimed to assess the bromatological characteristics of several species from the native and endemic flora of Fuerteventura, focusing on key quality indicators like protein, fiber content, and lipid profiles, to identify forage crops with high nutritional quality traits that could enhance livestock productivity. The working hypothesis posits that evaluating nutritional qualities can serve as a tool to select the best single or combined candidate species from the native flora for potential cultivation in the arid and semi-arid zones of the Canary Archipelago.

2 Material and methods

2.1 Plant tissues sampling

The plants were grown at the facilities of a local government experimental farm located in the southwest of Fuerteventura Island (Canary Islands, Spain). This farm features a 1,200 m² mesh nursery dedicated exclusively to the reproduction of native flora (Supplementary Figure S1). Twelve native and endemic species from the island of Fuerteventura, regarded in oral tradition as goat feed, were selected for the present study. These included four legumes (Lotus lancerottensis, Bituminaria bituminosa, Coronilla viminalis, and Retama rhodorhizoides), three grasses (Cenchrus ciliaris, Phalaris coerulescens, and Tricholaena teneriffae), and five other herbaceous and shrub species from different families denominated as mixed group (Lavatera acerifolia, Periploca laevigata, Campylanthus salsoloides, Echium decaisnei, and Crambe sventenii) (Figure 1).

Figure 1

Biota (2024) (C. viminalis, L. lancerottensis, P. coerulescens, C. salsoloides, C. sventenii, E. decaisnei, L. acerifolia, P. laevigata, by © G. García Casanova; B. bituminosa, R. rhodorhizoides, C. ciliaris, P. laevigata by © E. Ojeda-Land; C. sventenii by © S. Scholz) and Endémicas Canarias (2024) (B. bituminosa, C. viminalis, R. rhodorhizoides, C. ciliaris, P. coerulescens, T. teneriffae by © G. García Casanova).

All plants were germinated under controlled conditions and grown in a substrate made from a mix of Fuerteventura soil (Typic Torrifluvents; Soil Survey Staff, 2014) with a sandy-loam texture (clay 198.1 g/kg, silt 594.0 g/kg, and sand 207.9 g/kg), peat, and gravel in a 1:1:1 ratio, and irrigated with desalinated brackish water (EC_iw ~ 0.4 dS/m). The soil had a very low organic matter and nitrogen content, typical of arid region soils, requiring the use of basal fertilization (~ 6 g of NPK 17:6:14 fertilizer per 3 L volume container) to prevent nutritional deficiencies.

Plant tissue sampling took place in early November 2021 when the plants were between 10-18 months old. Samples (mainly leaves) from five individuals of each of the twelve selected species were collected. Once in the laboratory, the samples were pre-washed with distilled water, and divided into two subsamples. One of the subsamples was dried at 60°C for 72 hours and ground (1 mm) for fiber, protein, and minerals analysis, while the other subsample was stored at -80°C until their lipid analysis.

2.2 Analysis of plant tissues

The nutritive composition of the leaf tissues was determined using standard procedures (AOAC, 1990, 2000). The nutritional potential of the forages was assessed based on the levels of CP, ash, neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), enzymatic digestibility of organic matter (EOMD), relative feed value (RFV), ME, lipid content, lipid classes and fatty acid profiles, and minerals. CP was determined by dried combustion and measurement in the resulting gases with a thermoconductivity cell from CN828 equipment (LECO Corporation, St. Joseph, MI, USA). The ash content after incineration at 550°C, NDF, and ADF were analyzed using the Van Soest fractionation analysis (Van Soest et al., 1991). ADL was extracted using 72% sulfuric acid after the ADF procedure described by Robertson and Van Soest (1981), and EOMD following de la Roza and Argamentería (1992). The RFV was calculated from the ADF and NDF as follows: RFV = [88.9 - (0.779 * ADF)] * [(120/NDF)/1.29] (Putnam et al., 2008). The ME was calculated from the ash content and the in vitro digestibility of organic matter (IVOMD) as follows: ME = K * (100 - ash) * (IVOMD/100), where K = 0.16 for forages and the IVOMD was predicted from EOMD values (ARC, 1984; de la Roza and Argamentería, 1992).

Lipid characterization was based on the content of total lipids (TL), lipid classes and FAs profiles. TL were extracted by homogenization in chloroform/methanol (2:1 v/v) containing 0.01% butylated hydroxytoluene (BHT) as an antioxidant according to the method described by Folch et al. (1957) with slight modifications (Christie and Han, 2010). The organic solvent was evaporated under a stream of nitrogen, and the lipid content was determined gravimetrically. The lipid residue was then resuspended in chloroform/methanol (2:1 v/v) containing 0.05% BHT and stored under an inert nitrogen atmosphere at -20°C until further analysis.

The lipid classes were determined by high-performance thin-layer chromatography (HPTLC) with dual one-dimensional development (Olsen and Henderson, 1989), using 1-propanol/chloroform/methyl acetate/methanol/potassium chloride at 0.25% (5:5:5:2:1.8, v/v) to separate the polar lipids (PL), and hexane/diethyl ether/acetic acid (20:5:0.5, v/v) for the neutral lipids (NL). Lipid classes were quantified by calibrated densitometry by means of a CAMAG TLC Visualizer dual-wavelength flying-spot scanner (Camag, Muttenz, Switzerland), as described by Reis et al. (2019). Lipid class identification was performed by comparison with external lipid standards (cod roe lipid extract and a mix of digalactosyldiacylglycerol (DGDG), monogalactosyldiacylglycerol (MGDG), and sulfoquinovosyldiacylglycerol (SQDG) (Avanti Polar Lipids, Inc., Alabaster, Alabama, USA) placed on the same HPTLC plate (Supplementary Figure S2). The proportion of each lipid class was expressed as a percentage of the TL in the analyzed sample.

Another aliquot of TL was subjected to acid-catalyzed transesterification according to Christie and Han (2010). The resulting fatty acid methyl esters (FAMEs) were purified by thin-layer chromatography (TLC) and separated and quantified using a TRACE-GC Ultra gas chromatograph (Thermo Fisher Scientific Inc., Milan, Italy) under the chromatographic conditions described in Galindo et al. (2022). Individual FAMEs were identified by reference to a well-characterized commercial standard mixture (Mix C4-C24 and PUFA No. 3 from menhaden oil (Supelco Inc., Merck KGaA, Darmstadt, Germany)), and their relative abundance was expressed as a percentage of the total FAs in the analyzed sample (Supplementary Figure S3). Before transesterification, nonadecanoic acid (19:0) was added to the lipid fraction as an internal standard. When necessary, the identity of FAME was confirmed by GC-MS (DSQII, Thermo Scientific), under the same chromatographic conditions.

Nutritional quality of plant FA composition was assessed by calculating the index of atherogenicity (IA), the index of thrombogenicity (IT) following Ulbricht and Southgate (1991), and the hypocholesterolemic (h) to hypercholesterolemic (H) FAs ratio (hH) as described by Santos-Silva et al. (2002):

The elements analyzed for the study of mineral composition were N, P, K, Ca, Mg, Na, Fe, Mn, Cu, Zn, and B. All elements were determined by digestion combined with inductively coupled plasma mass spectrometry (ICP-MS), except N, which was measured using LECO CN828 equipment. The concentrations of these elements were compared with the established maximum tolerable levels (MTLs), defined as the highest level of a particular element in an animal’s diet over a specified period that does not lead to health deterioration or a decrease in animal performance (NRC, 2005).

2.3 Statistical analysis

Statistical methods were implemented using R version 4.3.0 (R Core Team, 2023) and IBM SPSS Statistics V26.0. The significance level for all tests was set at p< 0.05. Prior to each analysis, the assumptions of normality (Shapiro-Wilk test) and homogeneity of variance (Levene’s test) were verified within groups, and, where necessary, appropriate variance stabilizing transformations were performed. For all parameters, differences between the plant species were evaluated using a one-way ANOVA followed by the Tukey HSD post-hoc test. When transformations failed, Welch and Brown-Forsythe test followed by the Dunnett T3 test, or non-parametric Kruskal Wallis and Mann-Whitney U tests with Bonferroni adjustment were performed. Principal component analysis (PCA) for the main nutritional parameters of plant species was carried out. Two hierarchical cluster heatmaps were created using Euclidean distance, one for the lipid classes and the other for the main FAs groups, to identify patterns among species.

3 Results

3.1 Forage fiber, protein and ash contents

Table 1 presents the main nutritional parameters related to fiber, protein and ash contents of the analyzed plant species. Levels of NDF, ADF, and ADL varied strongly between families and species (Table 1). Poaceae species exhibited higher average values of NDF and ADF compared to the other groups (e.g., 64.8, 41.4 and 28.0% for NDF in Poaceae, Fabaceae and mixed groups, respectively). T. teneriffae was the grass species with the highest concentrations of NDF and ADF among all analyzed plants. In contrast, ADL levels were up to two-three times higher in R. rhodorhizoides and E. decaisnei than in the other species. Content of CP varied significantly from ~ 6% (T. teneriffae) to ~ 32% (P. coerulescens), both from the Poaceae family. Generally, legumes surpassed grasses in terms of CP, except for P. coerulescens which presented the highest proportions of all analyzed species. Species within the mixed group showed CP levels close to those of legumes (average values 14.2 vs. 17.5%), with C. salsoloides reaching an average concentration of 18.9% (Table 1). Ash content ranged from ~ 5% (R. rhodorhizoides) to 19% (C. sventenii). The average EOMD for all species in the mixed group was 80.0%, which was significantly higher compared to legumes (64.1%) and grasses (40.1%). C. sventenii, B. bituminosa, and C. ciliaris notably stood out within their respective groups. The energy potential of forage, calculated as ME, ranged from 8.1 to 8.6 MJ/kg DM on average per group. Among all studied species, B. bituminosa recorded the highest significant values (~ 10.9 MJ/kg DM), unlike specimens of R. rhodorhizoides which averaged 5.9 MJ/kg DM. Levels of RFV ranged from 75 in T. teneriffae to 318 in C. salsoloides. The mixed group had the highest average RFV indices (260), followed by Fabaceae (165) and Poaceae (108).

3.2 Forage lipid composition

3.2.1 Total lipid content

The total lipid (TL) content of the analyzed species ranged between 1.5 and 6.3% of DM (Figure 2). Within the legume group, TL tended to be higher in C. viminalis compared to the rest (> 5% vs. 3-4%, respectively) whereas all grasses presented values from 1.5 to 3% of DM, with significantly lower levels in C. ciliaris. C. salsoloides showed the highest TL proportions (> 6%) of all analyzed plants.

Figure 2

The multivariate PCA depicted in Figure 3 combines the above-mentioned nutritional factors and confirms the observations regarding nutritional differences between the analyzed species. Based on this analysis, 68.67% of the variance was explained by the two principal components. The first component (PC1) showed a strong positive influence with RFV and EOMD (score = 1.936 and 1.920), while ADF and NDF were the variables with the highest negative values (-1.956 and -1.954, respectively). In the second component (PC2), the variables with the highest weight were ME (score = 1.758) and ADL (score = -1.776). Poaceae species were clearly differentiated from the others, mainly due to the high levels of NDF and ADF. Species from the mixed group showed a closer relationship with EOMD and RFV, except for P. laevigata and E. decaisnei. This distinction was mainly due to higher levels of ADL and ADF in those two species. Most legumes showed a positive correlation with ME and CP, with the exception of R. rhodorhizoides, which also exhibited high levels of ADL and ADF (Figure 3).

Figure 3

3.2.2 Lipid classes profiles

The total lipids of the analyzed species presented highly variable proportions of total polar lipids (TPL; from 38 to 59% of TL) and total neutral lipids (TNL; from 41 to 62% of TL) (Supplementary Table S1). Concentration of TPL was greater than TNL, except B. bituminosa, in P. coerulescens (Poaceae) and C. salsoloides (Boraginaceae) (Supplementary Table S1). Among the PL, three main fractions clearly stood out: MGDG (8.4-22.5% of TL), DGDG (8.8-17.9% of TL), and SQDG (5.3-12.3% TL). In legumes MGDG levels were approximately 7.5-fold higher than in grasses, whereas SQDG + PE were up to four times lower in legumes. By contrast, DGDG proportions were similar in both families. Regarding TNL, free fatty acids (FFA) were generally more abundant in legumes than in grasses (Supplementary Table S1). Significant quantities of phytosterols (PTS) were present in all analyzed species (> 5% of TL), particularly in R. rhodorhizoides, C. ciliaris, T. teneriffae, and C. seventenii where average amounts up 10% of TL were registered.

A heatmap of lipid classes (Figure 4) identified three clusters based on the abundance patterns in the analyzed plants. Species in cluster 1 (e.g., C. sventenii and R. rhodorhizoides) were dominated by SQDG + phosphatidylethanolamine (PE), DGDG, and PTS compounds. Cluster 2 (i.e., L. acerifolia and C. ciliaris) stood out for their high contents of MGDG. The abundance of phosphatidylcholine (PC), phosphatidylglycerol (PG), and diacylglycerols (DAG) + pigments (P) was characteristic of species in cluster 3 (e.g., P. coerulescens and C. salsoloides).

Figure 4

3.2.3 Fatty acids profiles

Fatty acids profiles and quality indices for the studied plant species are shown in Supplementary Tables S2, S3, respectively. In all species, FA groups were predominantly represented by PUFA (from 46.5 to 63.1% of total FA), followed by saturated fatty acids (SFA; from 24.3 to 33.7% of total FA), and to a lesser extent, monounsaturated fatty acids (MUFA; from 7.1 to 18.4% of total FA). An exception was E. decaisnei, which had higher levels of SFA than PUFA (43.6 vs. 34.7%, respectively). Palmitic acid (16:0) was the most abundant saturated fatty acid, representing more than 50% of the SFA in all cases. Within PUFA, linoleic acid (LA; 18:2 n-6) and alpha-linolenic acid (ALA; 18:3 n-3) dominated, comprising up to 25% and 55% of total FAs, respectively (Supplementary Table S2). It is noteworthy that E. decaisnei was the only species where gamma-linolenic acid (GLA; 18:3 n-6) was detected, with an average value of 0.8 ± 0.2% of total FA and where the amount of stearidonic acid (SDA; 18:4n-3) was also prominent (2.3± 0.3% of total FA), compared to the other species (Supplementary Table S2). N-3/n-6 PUFA ratio was > 1 in all the analyzed species, being especially relevant in B. bituminosa, P. coerulescens, and C. salsoloides (5.4 ± 0.3, 8.1 ± 0.5 and 5.1 ± 0.7, respectively) (Supplementary Table S2).

The heatmap of FA families highlighted three familial clusters (Figure 5). In cluster 1, high concentrations of SFA and MUFA dominated, with species such as T. teneriffae and C. sventenii exhibiting the highest values, respectively. Cluster 2 is distinguished by an abundance of SFA, with E. decaisnei being the species with the greatest representation of this component. Finally, cluster 3 is characterized by a predominance of PUFA, n-3, and n-6 (e.g, L. acerifolia and P. coesulescens) displaying the highest average concentrations.

Figure 5

B. bituminosa, R. rhodorhizoides, P. coerulescens, and C. salsoloides showed the lowest IA and IT indices (~ 0.3 and 0.1, respectively) of all species. Regarding the hH ratio, R. rhodorhizoides, T. teneriffae, C. sventenii, and P. laevigata presented ratios ≥ 1, with the legume species showing the highest hH value at 1.7 (Supplementary Table S3).

3.3 Forage mineral composition

The mean values of macronutrients showed significant variations among species (Table 2). Thus, Ca concentrations were higher in legumes than in grasses (8.8 and 4.7 g/kg, respectively), but the highest concentrations were reached in E. decaisnei and P. laevigata (up to 22.7 and 34.2 g/kg, respectively). Regardless of the plant group, Mg levels remained low with average values around 3 g/kg. Potassium dominated in E. decaisnei and C. viminalis, not exceeding 35 g/kg in any specimen. The lowest K values were found in C. salsoloides, C. sventenii, R. rhodorhizoides, and T. teneriffae (ranging from 3.6 to 8 g/kg). Sodium concentrations were high for E. decaisnei and P. laevigata (~ 26 g/kg), followed by two legumes species C. viminalis and L. lancerottensis. not exceeding 5 g/kg in the remaining species. Content of P was similar in all species, ranging between 0.3 and 2.3 g/kg, except in E. decaisnei which presented the highest average value of 3.4 g/kg. Nitrogen followed the same pattern as CP, peaking in P. coerulescens (Table 2).

The mean micronutrient concentrations in the collected plant species are presented in Table 3. Similarly to macroelements, significant differences were detected between groups and species. Iron levels showed high variability, with differences of up to approximately two orders of magnitude, led by L. lancerottensis and at the opposite end C. salsoloides (803.6 vs. 9.3 mg/kg). In most species, average Mn values fluctuated between ~ 46 and 81 mg/kg, except for R. rhodorhizoides, P. coerulescens, C. salsoloides and C. sventenii where it accounted for 1-2 mg/kg. Copper and Zn followed similar trends, with the highest concentrations in legumes and mixed group (average ~ 37.0 and 45.3 mg/kg for Cu and Zn, respectively). In both cases, E. decaisnei presented the highest concentrations of these elements. Levels of B also exhibited high variability, ranging from less than 10 mg/kg in grasses such as P. coerulescens and T. teneriffae to levels exceeding 100 mg/kg in the legume C. viminalis (Table 3).

4 Discussion

4.1 Implications of the fiber, protein and ash profile on forage quality

Forage quality generally improves with increases in CP, EOMD, RFV, and ME, and reductions in NDF, ADF, and ash content (Suyama et al., 2007a; Díaz et al., 2018; Pirasteh-Anosheh et al., 2023). The levels of NDF, ADF and ash found in our present work coincide largely with those reported by other authors for species such as B. bituminosa, C. viminales, L. lancerottensis and E. decaisnei under natural conditions (Chinea et al., 2014). However, CP values were slightly higher than those previously reported, possibly due to the more favorable nursery conditions. When compared to a protein-rich forage like alfalfa (15-20% CP; Annicchiarico et al., 2015; Yang et al., 2021), four of the studied species (L. lancerottensis, B. bituminosa, C. viminalis and C. salsoloides) exhibited CP levels within this range (Table 1), demonstrating their potential as high-quality forage. The critical threshold of CP required for maintaining an adult non-pregnant sheep is approximately 7-9% (Pirasteh-Anosheh et al., 2023), and all examined native species, except T. teneriffae, surpassed this minimum, reinforcing their suitability for forage based on nitrogen and protein benchmarks.

The fiber percentages in our study (24.2-64.5% NDF and 9.0-40.5% ADF) also reflect the nutritional diversity within the studied forages. Eight out of twelve native species analyzed had an NDF content below 45%, and nine out of twelve had an ADF content below 35%, positioning them well within the ranges considered suitable for high-quality forage (Pirasteh-Anosheh et al., 2023). This has a positive impact on intake, energy density, milk production, cattle health, and feed cost (Christensen et al., 2015; Anderson et al., 2023). All forages had acceptable levels of ME (i.e., > 7 MJ/kg DM in most cases) for beef cattle and goats fed at a maintenance energy level (NRC, 1996). Considering that a 53 kg sheep needs forage with an ME content of 10-12 MJ/kg (Pirasteh-Anosheh et al., 2023)., only B. bituminosa met this requirement. This matches findings from Pirasteh-Anosheh et al. (2023) who reported that none of the selected halophyte species reached the required ME content.

The most valuable species in terms of forage quality within legumes were B. bituminosa and C. viminalis, both P. coerulescens and C. ciliaris within grasses, and C. salsoloides, E. decaisnei and C. sventenii for the mix group. Another criterion based on the quality parameters analyzed indicates that legumes can be classified as “premium” quality, grasses as “good-utility” and the rest between “premium” and “good” (Putnam et al., 2008; USDA Agricultural Marketing Service, 2024), which makes most of them suitable for dairy cattle and for low production or maintenance level beef cattle (Suyama et al., 2007b), reflecting their high nutritional content.

4.2 Evaluation of lipids and their impact on forage quality, animal and human health

The TL content of the studied plants was in line with the generally low lipid levels of forages stated in the literature (< 8% DM; Davis et al., 2022). In recent years, there has been growing interest in under-exploited but promising plant species as alternative sources of vegetable oils. Some of these species can contain significant amounts of oil and, in addition, a high proportion of nutritionally desirable fatty acids, such as omega-3 PUFA (Manso et al., 2009; Timilsena et al., 2017). The PUFA levels recorded in our work were slightly higher than those determined by Cerci et al. (2011) in fresh alfalfa forages (30-50%), except in the case of P. laevigata, E. decaisnei, and C. sventenii (Figure 4). Previous studies have provided evidence of the differences in the FA profiles among various legumes, grasses, and herbaceous forage species, in which grasses tend to have higher levels of the C18 omega-3 PUFA, ALA than legumes, and conversely, lower contents of the C18 omega-6 PUFA, LA (Clapham et al., 2005; Elgersma et al., 2013; Whetsell and Rayburn, 2022). In this context, both C18 PUFA have been described to exert several benefits to animal health and well-being, and must be incorporated through the diet as they cannot be synthesized de novo (Morales-Almaráz et al., 2011; Davis et al., 2022). However, the increasing concern on the excessive intake of omega-6 fatty acids in the occidental diet, have raised the interest to search for plants and forage particularly rich in omega-3 fatty acids, with higher n-3/n-6 ratios, both for human and farmed animals feeding (Almeida et al., 2019; Zárate et al., 2017). It should also be stressed that E. decaisnei was the only species containing substantial amounts of two fatty acids of increasing nutraceutical value: GLA (0.8% of total FA) and SDA (2.3 of total FA). In this regard, species of the Echium genus are known to be rich sources of GLA and SDA (Guil-Guerrero et al., 2000a, b), two physiologically valuable FA for mammals’ health (López-Martínez et al., 2004; Wang et al., 2022; Berti et al., 2007), that have shown promising results in the cultivation of aquatic species (Díaz-López et al., 2009, 2010) and terrestrial farmed animals (Navarro et al., 2008).

It is well established that in the feeding of ruminants, the most important lipids are those containing fatty acids bound to glycerol (i.e., phospholipids, glycolipids and triglycerides). Glycerolipids are the most abundant type of lipid in our forages species, represented by major phospholipids including PC, phosphatidylserine (PS), phosphatidylinositol (PI), PG, and PE, and galactolipids such as SQDG, MGDG, and DGDG (Liu et al., 2017). Among the NLs, DAG and TAG were prominent. Like FAs, the relative proportions of these compounds fluctuate among plant species and are strongly influenced by the nutrient levels in the plants (Liu et al., 2017). SQDG, MGDG, and DGDG have been described as compounds with anti-inflammatory and antithrombotic properties (Ibañez and Cifuentes, 2013). Similarly, the levels of PTS in the forage can be very important as it has been demonstrated to reduce cholesterol levels in the form of low-density lipoprotein (LDL) in humans, thereby reducing the risk of cardiovascular diseases and inflammatory processes (Ibañez and Cifuentes, 2013; Kendel et al., 2015; Simonetti et al., 2021).

4.3 Implications of mineral composition in ruminant diets

The mineral composition of forage plant tissue shows significant variations, which could have direct implications on dietary formulations for ruminants. The levels of macro- and micro-minerals generally fall within the range that meets ruminants’ nutritional requirements. Species such as C. viminalis, L. lancerottensis, E. decaisnei, L. acerifolia, and P. laevigata had mineral contents that exceed the normally adequate values, particularly for Ca, Na, Fe and Cu (NRC, 2005). High Ca levels may cause health risks such as milk fever in dairy cows (Lean et al., 2006). The MTLs for dietary Ca in ruminants is 1.5% DM (NRC, 2005), yet E. decaisnei and P. laevigata recorded Ca concentrations of 22.7 and 34.2 g/kg DM, respectively, which can depress feed intake and milk yield. Sodium levels in some specimens reached 28 g/kg DM, potentially reducing appetite and increasing urine output, as well as the risk of kidney damage (ARC, 1980; NRC, 2007). Excessive dietary intake of Na is usually manageable as long as non-saline water is readily available for animals (Meyer et al., 1955; Nelson et al., 1955). Iron toxicity could lead to diarrhea and metabolic acidosis, in our case only L. lancerottensis exceeded the MTLs for this element (~ 500 mg/kg DM; NRC, 2005). Additionally, E. decaisnei exhibited Cu levels slightly above the safe threshold at 40 mg/kg DM (NRC, 2005), which could result in liver necrosis. Chronic Cu toxicity, particularly challenging to diagnose in cattle, involves symptoms such as hemoglobinuria and jaundice, typically more obvious in sheep, and is often indicated by high liver, kidney, and serum Cu concentrations (López-Alonso and Miranda, 2020).

Our results reveal strong interspecific variations in all parameters, indicating the potential benefits of using a diverse mix of species to optimize livestock nutritional intake. This is especially crucial in arid regions where plant species naturally exhibit unique nutritional profiles (Lee, 2018). Selecting appropriate forage species can optimize mineral intake and prevent deficiencies or toxicities that could impair livestock health and productivity.

Globally, most of the evaluated species have an adequate nutritional value, in addition to presenting a FA profile and certain lipid classes, such as PTS, beneficial for animal health and well-being. In particular, the species B. bituminosa, P. coerulescens, E. decaisnei, and C. sventenii could be promising candidates to be introduced as part of the diet of ruminants. The present work is limited to evaluating plant specimens under controlled conditions, so future research should focus on studying these species in field conditions, under different agronomic managements, that allow evaluating the influence of seasonal and environmental factors on the nutritional composition.

5 Conclusions

This study demonstrates that native and endemic plant species from arid environments like Fuerteventura Island have nutritional profiles capable of meeting the dietary requirements of livestock, making them viable forage options. The bromatological characterization of the twelve selected species shows they offer a sustainable alternative to conventional forages with potential nutritional benefits for livestock health and productivity. However, the strong interspecific variability suggests the combined use of different species as livestock feed.

By selecting native species with the most optimal nutritional profiles, it is possible to cultivate forages that not only withstand extreme conditions but also enhance livestock productivity and sustainability. This approach could be especially beneficial to regions heavily reliant on imported forages, thus reducing environmental impact and improving food security. Overall, this study provides a crucial foundation for the strategic selection and cultivation of native forage species in arid regions. Our findings advocate for an integrated approach that includes these native species in the agricultural matrix of arid regions, which could serve a dual purpose: enhancing livestock nutrition and conserving the unique botanical heritage of these ecosystems. Special attention should be paid to endemic species, some of which are endangered, as the discovery of potential benefits of their consumption in livestock health could contribute to their cultivation and simultaneously to their conservation. This paradigm shift toward the use of local and biodiverse plant species as forage represents a crucial step towards ecological sustainability, highlighting the untapped potential of arid region flora and paving the way for future agricultural innovations that address the nutritional needs of livestock in challenging environments.

References

• 1

  AlmeidaO. C.FerrazM. V.Jr.SusinI.GentilR. S.PolizelD. M.FerreiraE. M.et al. (2019). Plasma and milk fatty acid profiles in goats fed diets supplemented with oils from soybean, linseed or fish. Small Rumin. Res.170, 125–130. doi: 10.1016/j.smallrumres.2018.11.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  ÁlvarezS.FresnoM.MéndezN.CastroN.FernándezJ. R.Sanz SampelayoM. R. (2007). Alternatives for improving physical, chemical, and sensory characteristics of goat cheeses: The use of arid land forage in the diet. J. Dairy Sci.90, 2181–2188. doi: 10.3168/jds.2006-506

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  Álvarez RíosS.MéndezP.FresnoM. (2018). Effect of feeding goats with leguminous shrubs (Chamaecytisus proliferus spp. palmensis and Bituminaria bituminosa) on milk and cheese properties. J. Appl. Anim. Res.46, 1443–1451. doi: 10.1080/09712119.2018.1531762

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  AndersonA. W.GullU.BenesS. E.SinghS.HutmacherR. B.BrummerE. C.et al. (2023). Salinity and cultivar effects on alfalfa forage yield and nutritive value in a Mediterranean climate. Grassl. Res.2, 153–166. doi: 10.1002/glr2.12061

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  AnnicchiaricoP.BarrettB.BrummerE. C.JulierB.MarshallA. H. (2015). Achievements and challenges in improving temperate perennial forage legumes. Crit. Rev. Plant Sci.34, 327–380. doi: 10.1080/07352689.2014.898462

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  AOAC (1990). Official Methods of Analysis of the Association of Official Analytical Chemists. 15th Edn (Arlington, VA, USA: Association of Official Analytical Chemists).

  • Google Scholar
  • View reference in article

• 7

  AOAC (2000). Official Methods of Analysis. 17th ed (Washington, DC, USA: AOAC).

  • Google Scholar
  • View reference in article

• 8

  ARC (1980). Agricultural Research Council. The Nutrient Requirements of Ruminant Livestock.Technical Review (Bureaux, Farnham Royal, UK: Commonwealth Agric.).

  • Google Scholar
  • View reference in article

• 9

  ARC (1984). Agricultural Research Council. The Nutrient Requirements of Ruminant Livestock. Supplement No1 (Bureaux, Slough, UK: Commonwealth Agric.).

  • Google Scholar
  • View reference in article

• 10

  BertiM.JohnsonB. L.DashS.FischerS.WilckensR.HeviaF. (2007). “Echium: A source of stearidonic acid adapted to thenorthern Great Plains in the US,” in Issues in new crops and new uses. Eds. JanickJ.WhipkeyA. (Alexandria, VA, USA: ASHS Press), 120–125.

  • Google Scholar
  • View reference in article

• 11

  Biota (2024). Banco de Datos de Biodiversidad de Canarias. Gobierno de Canarias. Available at: https://www.biodiversidadcanarias.es/biota (Accessed February 22, 2024).

  • Google Scholar
  • View reference in article

• 12

  CapstaffN. M.MillerA. J. (2018). Improving the yield and nutritional quality of forage crops. Front. Plant Sci.9. doi: 10.3389/fpls.2018.00535

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  Castro-MontoyaJ. M.DickhoeferU. (2020). The nutritional value of tropical legume forages fed to ruminants as affected by their growth habit and fed form: A systematic review. Anim. Feed Sci. Technol.269, 114641. doi: 10.1016/j.anifeedsci.2020.114641

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  Caujapé-CastellsJ.García-VerdugoC.SanmartínI.Fuertes-AguilarJ.RomeirasM. M.Zurita-PérezN.et al. (2022). The late Pleistocene endemicity increase hypothesis and the origins of diversity in the Canary Islands Flora. J. Biogeogr.49, 1469–1480. doi: 10.1111/jbi.14394

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  CercıI. H.CıftcıM.BahsıM.KılıncU. (2011). Cholesterol and fatty acid composition of lamb serum and offal as affected by alfalfa and concentrate. Vet. Arh.81, 575–584.

  • Google Scholar
  • View reference in article

• 16

  ChandS.SinghalR. K.GovindasamyP. (2022). Agronomical and breeding approaches to improve the nutritional status of forage crops for better livestock productivity. Grass Forage Sci.7, 11–32. doi: 10.1111/gfs.12557

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  ChineaE.BatistaC.MoraJ. L.García-CiudadA. (2014). Calidad nutritiva de especies autóctonas de Canarias con interés forrajero para zonas áridas. Rev. la Facultad Agronomía31, 100–128. Available at: https://produccioncientificaluz.org/index.php/agronomia/article/view/27151.

  • Google Scholar
  • View reference in article

• 18

  ChristensenR. G.YangS. Y.EunJ. S.YoungA. J.HallJ. O.MacAdamJ. W. (2015). Effects of feeding birdsfoot trefoil hay on neutral detergent fiber digestion, nitrogen utilization efficiency, and lactational performance by dairy cows. J. Dairy Sci.98, 7982–7992. doi: 10.3168/jds.2015-9348

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  ChristieW. W.HanX. (2010). Lipid analysis: isolation, separation, identification and lipidomic analysis. Bridgwater: Oily Press an imprint PJ Barnes Associates. pp.55–56. doi: 10.1533/9780857097866

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  ClaphamW. M.FosterJ. G.NeelJ. P. S.FeddersJ. M. (2005). Fatty acid composition of traditional and novel forages. J. Agric. Food Chem.53, 10068–10073. doi: 10.1021/jf0517039

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  ColonnaM. A.GiannicoF.TufarelliV.LaudadioV.SelvaggiM.De MastroG.et al. (2021). Dietary Supplementation with camelina sativa (L. Crantz) forage in autochthonous ionica goats: effects on milk and caciotta cheese chemical, fatty acid composition and sensory properties. Animals11, 1589. doi: 10.3390/ani11061589

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  DavisH.MagistraliA.ButlerG.StergiadisS. (2022). Nutritional benefits from fatty acids in organic and grass-fed beef. Foods11, 646. doi: 10.3390/foods11050646

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  de la RozaB.ArgamenteríaA. (1992). Organic matter digestibility of forages for ruminants. Prediction by sodium hydroxide or sodium chlorite pretreatment in conjunction with enzymatic treatment. A. Biotech. Forum Europe9, 294–297.

  • Google Scholar
  • View reference in article

• 24

  DíazF. J.GrattanS. R.ReyesJ. A.de la Roza-DelgadoB.BenesS. E.JiménezC.et al. (2018). Using saline soil and marginal quality water to produce alfalfa in arid climates. Agric. Water Manage.199, 11–21. doi: 10.1016/j.agwat.2017.12.003

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  Díaz-LópezM.PérezM. J.AcostaN. G.JerezS.Dorta-GuerraR.TocherD. R.et al. (2010). Effects of dietary fish oil substitution by Echium oil on enterocyte and hepatocyte lipid metabolism of gilthead seabream (Sparus aurata L.). Comp. Biochem. Physiol. B. Biochem. Mol. Biol.155, 371–379. doi: 10.1016/j.cbpb.2009.12.004

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  Díaz-LópezM.PérezM. J.AcostaN. G.TocherD. R.JerezS.LorenzoA.et al. (2009). Effect of dietary substitution of fish oil by Echium oil on growth, plasma parameters and body lipid composition in gilthead seabream (Sparus aurata L.). Aquac. Nutr.15, 500–512. doi: 10.1111/j.1365-2095.2008.00616.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  ElgersmaA.SøegaardK.JensenS. K. (2013). Fatty acids, α-tocopherol, β-carotene, and lutein contents in forage legumes, forbs, and a grass-clover mixture. J. Agric. Food Chem.61, 11913–11920. doi: 10.1021/jf403195v

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  El ShaerH. M. (2010). Halophytes and salt-tolerant plants as potential forage for ruminants in the Near East region. Small Rumin. Res.91, 3–12. doi: 10.1016/j.smallrumres.2010.01.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  Endémicas Canarias (2024). Available at: https://endemicascanarias.com/es/ (Accessed December 10, 2024)

  • Google Scholar
  • View reference in article

• 30

  FolchJ.LeesM.StanleyG. H. S. (1957). A simple method for the insolation and purification of total lipides from animal tissues. J. Biol. Chem.226, 497–509. doi: 10.1016/S0021-9258(18)64849-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  GalindoA.ReisD. B.RodríguezI.PérezJ. A.Abdul-JalbarB.ZárateR.et al. (2022). Lipid characterization of 14 macroalgal species from Madeira Archipelago: Implications for animal and human nutrition. Bot. Mar.65, 51–67. doi: 10.1515/bot-2021-0024

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  Guil-GuerreroJ. L.García-MarotoF.Campra-MadridP.Gómez-MercadoF. (2000b). Occurrence and characterization of oils rich in gamma linolenic acid Part II: fatty acids and squalene from Macaronesian Echium leaves. Phytochem.54, 525–529. doi: 10.1016/S0031-9422(00)00139-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  Guil-GuerreroJ. L.Gómez-MercadoF.García-MarotoF.Campra-MadridP. (2000a). Occurrence and characterization of oils rich in gamma-linolenic acid Part I: Echium seeds from Macaronesia. Phytochem.53, 451–456. doi: 10.1016/S0031-9422(99)00549-X

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  HunterM. C.SmithR. G.SchipanskiM. E.AtwoodL. W.MortensenD. A. (2017). Agriculture in 2050: Recalibrating targets for sustainable intensification. BioSci.67, 386–391. doi: 10.1093/biosci/bix010

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  IbañezE.CifuentesA. (2013). Benefits of using algae as natural sources of functional ingredients. J. Sci. Food Agric.93, 703–709. doi: 10.1002/jsfa.6023

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  ISTAC (2022). Instituto Canario de Estadística. Available online at: http://www.gobiernodecanarias.org/istac/ (Accessed September 1, 2022).

  • Google Scholar
  • View reference in article

• 37

  KendelM.Wielgosz-CollinG.BertrandS.RoussakisC.BourgougnonN. B.BedouxG. (2015). Lipid composition, fatty acids and sterols in the seaweeds Ulva armoricana, and Solieria chordalis from Brittany (France): An analysis from nutritional, chemotaxonomic, and antiproliferative activity perspectives. Mar. Drugs13, 5606–5 28. doi: 10.3390/md13095606

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  LeanI. J.DeGarisP. J.McNeilD. M.BlockE. (2006). Hypocalcemia in dairy cows: Meta-analysis and dietary cation anion difference theory revisited. J. Dairy Sci.89, 669–68. doi: 10.3168/jds.S0022-0302(06)72130-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  LeeM. A. (2018). A global comparison of the nutritive values of forage plants grown in contrasting environments. J. Plant Res.131, 641–654. doi: 10.1007/s10265-018-1024-y

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  LiuM. Y.BurgosA.MaL.ZhangQ.TangD.RuanJ. (2017). Lipidomics analysis unravels the effect of nitrogen fertilization on lipid metabolism in tea plant (Camellia sinensis L.). BMC Plant Biol.17, 165. doi: 10.1186/s12870-017-1111-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  López-AlonsoM.MirandaM. (2020). Copper supplementation, A challenge in cattle. Animals10, 1890. doi: 10.3390/ani10101890

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  López-MartínezJ. C.Campra-MadridP.Guil-GuerreroJ. L. (2004). [amp]]gamma;Linolenic acid enrichment from Borago officinalis and Echium fastuosum seed oils and fatty acids by low temperature crystallization. J. Biosci. Bioeng.97, 294–298. doi: 10.1016/S1389-1723(04)70208-X

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  LouhaichiM.GamounM.HassanS.AbdallahM. A. B. (2021). Characterizing biomass yield and nutritional value of selected indigenous range species from arid Tunisia. Plants10, 2031. doi: 10.3390/plants10102031

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  MansoT.BodasR.CastroT.JimenoV.ManteconA. R. (2009). Animal performance and fatty acid composition of lambs fed with different vegetable oils. Meat Sci.83, 511–516. doi: 10.1016/j.meatsci.2009.06.035

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  MéndezP.ÁlvarezS.DíazA.DíazC.FresnoM. (2003). Nutritive value of vinagrera (Rumex lunaria), an endemic forage shrub of Canary Islands. Grassland Sci. Europe8, 229–231.

  • Google Scholar
  • View reference in article

• 46

  MeyerJ. H.WeirW. C.IttnerN. R.SmithJ. D. (1955). The influence of high sodium chloride intake by fattening sheep and cattle. J. Anim. Sci.14, 412–418. doi: 10.2527/jas1955.142412x

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  Morales-AlmarázE.de la Roza-DelgadoB.GonzálezA.SoldadoA.RodríguezM. L.PeláezM.et al. (2011). Effect of feeding system on unsaturated fatty acid level in milk of dairy cows. Renew. Agr. Food Syst.26, 224–229. doi: 10.1017/S1742170511000019

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  NavarroM. C. G.Ramos MoralesE.de la Torre AdarveG.Fernández NavarroJ. R.Rodríguez OsorioM.Gil ExtremeraF.et al. (2008). Growth of pre-ruminant kid goats and the composition of carcass fat deposits: effects of providing a PUFA-rich fat in the milk replacer and influence of the kidding season. Food Sci.Techn. Int.14, 85–94. doi: 10.1177/1082013208094680

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  NelsonA. B.MacVicarR. W.ArcherW.Jr.MeiskeJ. C. (1955). Effect of a high salt intake on the digestibility of ration constituents and on nitrogen, sodium, and chloride retention by steers and wethers. J. Anim. Sci.14, 825–830. doi: 10.2527/1955.143825x

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  NormanH. C.MastersD. G. (2023). Livestock preference and feeding value as key determinants for forage improvement - why not ask the consumers? Anim. Prod. Sci.63, 1161–1176. doi: 10.1071/AN23140

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  NRC (1996). National Research Council. Nutrient Requirements of Beef Cattle. seventh ed (Washington, DC, USA: National Academy Press).

  • Google Scholar
  • View reference in article

• 52

  NRC (2005). National Research Council. Mineral Tolerance of Animals. 2nd rev. ed (Washington, DC: The National Academic Press).

  • Google Scholar
  • View reference in article

• 53

  NRC (2007). National Research Council. Nutrient requirement of Small Ruminants: sheep, goats, cervids and new world camelids (Washington DC: National Academic Press).

  • Google Scholar
  • View reference in article

• 54

  OlsenR. E.HendersonR. J. (1989). The rapid analysis of neutral and polar marine lipids using double-development HPTLC and scanning densitometry. J. Exp. Mar. Biol. Ecol.129, 189–197. doi: 10.1016/0022-0981(89)90056-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  PFORCA (2014). Plan Forrajero de Canarias. Available online at: https://www.gobiernodecanarias.org/ (Accessed September 1, 2022).

  • Google Scholar
  • View reference in article

• 56

  Pirasteh-AnoshehH.RanjbarG.AkramN. A.GhafarM. A.PanicoA. (2023). Forage potential of several halophytic species grown on saline soil in arid environments. Environ. Res.219, 114954. doi: 10.1016/j.envres.2022.114954

  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  PutnamD. H.RobinsonP.DePetersE. (2008). “Forage quality and testing,” in Irrigated Alfalfa Management for Mediterranean and Desert Zones. Chapter 16. Eds. SummersC. G.PutnamD. H. (Oakland, CA, USA: University of California/Agricultural and Natural Resources).

  • Google Scholar
  • View reference in article

• 58

  R Core Team (2023). R: A language and environment for statistical computing. Available online at: https://www.r-project.org/ (Accessed January 18, 2023).

  • Google Scholar
  • View reference in article

• 59

  ReisD. B.AcostaN. G.AlmansaE.GarridoD.AndradeJ. P.SykesA. V.et al. (2019). Effect of Artemia inherent fatty acid metabolism on the bioavailability of essential fatty acids for Octopus vulgaris paralarvae development. Aquaculture500, 264–271. doi: 10.1016/j.aquaculture.2018.10.021

  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  RobertsonJ. B.Van SoestP. J. (1981). “The detergent system of analysis and its application to human foods,” in The Analysis of Dietary Fiber in Foods. Eds. JamesW. P. T.TheanderO. (Marcel Dekker, New York, NY), 123–158.

  • Google Scholar
  • View reference in article

• 61

  Santos-SilvaJ.BessaR. J. B.Santos-SilvaF. (2002). Effect of genotype, feeding system and slaughter weight on the quality of light lambs: II. Fatty acid composition of meat. Livest. Prod. Sci.77, 187–194. doi: 10.1016/S0301-6226(02)00059-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 62

  SimonettiG.Di FilippoP.PomataD.RiccardiC.BuiarelliF.SonegoE.et al. (2021). Characterization of seven sterols in five different types of cattle feedstuffs. Food Chem.340, 127926. doi: 10.1016/j.foodchem.2020.127926

  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  SuyamaH.BenesS. E.RobinsonP. H.GetachewG.GrattanS. R.GrieveC. M. (2007b). Biomass yield and nutritional quality of forage species under long-term irrigation with saline-sodic drainage water: field evaluation. Anim. Feed Sci. Technol.135, 329–345. doi: 10.1016/j.anifeedsci.2006.08.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 64

  SuyamaH.BenesS. E.RobinsonP. H.GrattanS. R.GrieveC. M.GetachewG. (2007a). Forage yield and quality under irrigation with saline-sodic drainage water: greenhouse evaluation. Agric. Water Manage.88, 159–172. doi: 10.1016/j.agwat.2006.10.011

  • CrossRef
  • Google Scholar
  • View reference in article

• 65

  TimilsenaY. P.WangB.AdhikariR.AdhikariB. (2017). Advances in microencapsulation of polyunsaturated fatty acids (PUFAs)-rich plant oils using complex coacervation: A review. Food Hydrocoll.69, 369–381. doi: 10.1016/j.foodhyd.2017.03.007

  • CrossRef
  • Google Scholar
  • View reference in article

• 66

  UlbrichtT. L. V.SouthgateD. A. T. (1991). Coronary heart disease: seven dietary factors. Lancet338, 985–992. doi: 10.1016/0140-6736(91)91846-M

  • CrossRef
  • Google Scholar
  • View reference in article

• 67

  USDA Agricultural Marketing Service. Available online at: https://www.ams.usda.gov/market-news/livestock-poultry-grain (Accessed February 22, 2024).

  • View reference in article

• 68

  Van SoestP. J.RobertsonJ. B.LewisB. A. (1991). Methods of dietary fiber, neutral detergent fiber, and non-starch polysaccharides in relation to animal nutrition. J. Dairy Sci.74, 3583–3597. doi: 10.3168/jds.S0022-0302(91)78551-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 69

  WangW.JinJ.XuH.ShiY.BoerschM.YinY. (2022). Comparative analysis of the main medicinal substances and applications of Echium vulgare L. and Echium plantagineum L.: A review. J. Ethnopharmacol.285, 114894. doi: 10.1016/j.jep.2021.114894

  • CrossRef
  • Google Scholar
  • View reference in article

• 70

  WhetsellM.RayburnE. (2022). Variation in fatty acids concentration in grasses, legumes, and forbs in the allegheny plateau. Agronomy12, 1693. doi: 10.3390/agronomy12071693

  • CrossRef
  • Google Scholar
  • View reference in article

• 71

  YangC.ZhangF.JiangX.YangX.HeF.WangZ.et al. (2021). Identification of genetic loci associated with crude protein content and fiber composition in alfalfa (Medicago sativa L.) using QTL mapping. Front. Plant Sci.12. doi: 10.3389/fpls.2021.608940

  • CrossRef
  • Google Scholar
  • View reference in article

• 72

  ZárateR.el Jaber-VazdekisN.TejeraN.PérezJ. A.RodríguezC. (2017). Significance of long chain polyunsaturated fatty acids in human health. Clin. Transl. Res.6, 1–19. doi: 10.1186/s40169-017-0153-6

  • CrossRef
  • Google Scholar
  • View reference in article