# Insight into the Influence of Cultivar Type, Cultivation Year, and Site on the Lignans and Related Phenolic Profiles, and the Health-Promoting Antioxidant Potential of Flax (Linum usitatissimum L.) Seeds

Laurine Garros, Samantha Drouet, Cyrielle Corbin, Cédric Decourtil, Thibaud Fidel, Julie Lebas de Lacour, Emilie Leclerc, Sullivan Renouard, Duangjai Tungmunnithum, Joël Doussot, et al.

► To cite this version:

Laurine Garros, Samantha Drouet, Cyrielle Corbin, Cédric Decourtil, Thibaud Fidel, et al. Insight into the Influence of Cultivar Type, Cultivation Year, and Site on the Lignans and Related Phenolic Profiles, and the Health-Promoting Antioxidant Potential of Flax (Linum usitatissimum L.) Seeds. Molecules, 2018, 23 (10), 10.3390/molecules23102636. [https://doi.org/10.3390/molecules23102636](https://doi.org/10.3390/molecules23102636)

HAL Id: hal-01939131

https://hal.science/hal-01939131

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Insight into the Influence of Cultivar Type, Cultivation Year, and Site on the Lignans and Related Phenolic Profiles, and the Health-Promoting Antioxidant Potential of Flax (Linum usitatissimum L.) Seeds

Laurine Garros<sup>1,2,3,†</sup>, Samantha Drouet<sup>1,2,†</sup>, Cyrielle Corbin<sup>1,2</sup>, Cédric Decourtil<sup>1,2</sup>, Thibaud Fidel<sup>1,2</sup>, Julie Lebas de Lacour<sup>1,2</sup>, Emilie A. Leclerc<sup>1,2</sup>, Sullivan Renouard<sup>1,2</sup>, Duangjai Tungmunnithum<sup>1,2,4</sup>, Joël Doussot<sup>1,2,5</sup>, Bilal Haider Abassi<sup>1,2,6</sup>, Benoit Maunit<sup>2,3</sup>, Éric Lainé<sup>1,2</sup>, Ophélie Fliniaux<sup>7</sup>, François Mesnard<sup>7</sup> and Christophe Hano<sup>1,2,*</sup>

Received: 23 September 2018; Accepted: 11 October 2018; Published: 14 October 2018

check for
updates

**Abstract:** Flaxseeds are a functional food representing, by far, the richest natural grain source of lignans, and accumulate substantial amounts of other health beneficial phenolic compounds (i.e., flavonols, hydroxycinnamic acids). This specific accumulation pattern is related to their numerous beneficial effects on human health. However, to date, little data is available concerning the relative impact of genetic and geographic parameters on the phytochemical yield and composition. Here, the major influence of the cultivar over geographic parameters on the flaxseed phytochemical accumulation yield and composition is evidenced. The importance of genetic parameters on the lignan accumulation was further confirmed by gene expression analysis monitored by RT-qPCR. The corresponding antioxidant activity of these flaxseed extracts was evaluated, both in vitro, using ferric reducing antioxidant power (FRAP), oxygen radical absorbance capacity (ORAC), and iron chelating assays, as well as in vivo, by monitoring the impact of UV-induced oxidative stress on the lipid membrane peroxidation of yeast cells. Our results, both the in vitro and in vivo studies, confirm that flaxseed extracts are an effective protector against oxidative stress. The results point out that secoisolariciresinol diglucoside, caffeic acid glucoside, and *p*-coumaric acid glucoside are the main contributors to the antioxidant capacity. Considering the health benefits of these compounds, the present study demonstrates that the flaxseed cultivar type could greatly influence

<footnote n="1">¹ Laboratoire de Biologie des Ligneux et des Grandes Cultures (LBLGC) EA1207 INRA USC1328, Plant LIGNANS Team, Université d'Orléans, 28000 Chartres, France; laurine.garros@univ-orleans.fr (L.G.); samantha.drouet@univ-orleans.fr (S.D.); cyrielle.corbin@univ-orleans.fr (C.C.); cedric.decourtil@univ-orleans.fr (C.D.); thibaud.fidel@univ-orleans.fr (T.F.); julie.lebas-de-lacour@univ-orleans.fr (J.L.d.L.); emilie.leclerc@univ-orleans.fr (E.A.L.); sullivan.renouard@univ-orleans.fr (S.R.); duangjai.tun@mahidol.ac.th (D.T.); joel.doussot@lecnam.net (J.D.); bhabbasi@qau.edu.pk (B.H.A.); eric.laine@univ-orleans.fr (É.L.)</footnote>

<footnote n="2">² COSM'ACTIFS, Bioactifs et Cosmétiques, CNRS GDR3711, 45067 Orléans Cedex 2, France; benoit.maunit@univ-orleans.fr</footnote>

<footnote n="3">³ Institut de Chimie Organique et Analytique (ICOA) UMR7311, Université d'Orléans-CNRS, 45067 Orléans CEDEX 2, France</footnote>

<footnote n="4">⁴ Department of Pharmaceutical Botany, Faculty of Pharmacy, Mahidol University, 447 Sri-Ayuthaya Road, Rajathevi, Bangkok 10400, Thailand</footnote>

<footnote n="5">5 Le CNAM, Ecole Sciences Industrielles et Technologies de l'Information (SITI), Chimie Alimentation Santé Environnement Risque (CASER), 75141 Paris Cedex 3, France</footnote>

<footnote n="6">⁶ Department of Biotechnology, Quaid-i-Azam University, 45320 Islamabad, Pakistan</footnote>

<footnote n="7">⁷ Biologie des Plantes et Innovation (BOPI) EA 3900, Université de Picardie Jules Verne, 80000 Amiens, France; ophelie.fliniaux@u-picardie.fr (O.F.); francois.mesnard@u-picardie.fr (F.M.)</footnote>

<footnote>* Correspondence: hano@univ-orleans.fr; Tel.: +33-237-309-753</footnote>

<footnote>† These two authors contributed equally to this work and should be considered both as first authors.</footnote>
the phytochemical intakes and, therefore, the associated biological activities. We recommend that this
crucial parameter be considered in epidemiological studies dealing with flaxseeds.

**Keywords:** cultivar; environment; flax; flavonol; genetic; hydroxycinnamic acid; lignan; seed

**1. Introduction**

The consumption of fruit, vegetables, and grains has been associated with lower risks of chronic and degeneration diseases [1]. Considering their numerous beneficial effects on human health, during the last decades, there has been an increasing interest in their uses, and flaxseeds are, therefore, considered as functional food [2]. Flaxseeds are the richest natural grain source of lignan and accumulate a substantial amount of other phenolic compounds (e.g., flavonols, hydroxycinnamic acids). In flaxseed, the foremost part of these phytochemicals is accumulated under the form of a macromolecular complex (also known as lignan macromolecule) composed of the lignan secoisolariciresinol diglucoside (SDG, Figure 1A) as the main component, and of flavonol herbacetin diglucoside (HDG, Figure 1A), as well as hydroxycinnamic acid derivatives: *p*-coumaric acid glucoside (CouG, Figure 1A), caffeic acid glucoside (CafG, Figure 1A), and ferulic acid glucoside (FerG, Figure 1A), ester-linked together to hydroxymethylglutaryl spacers (Figure 1B) [3,4].

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Structure of phenolic compounds involved in the lignan macromolecular complex. (A) Structure of the complex components: (1) secoisolariciresinol (R = H), or secoisolariciresinol diglucoside (SDG, R = *β*-d-glucose), (2) herbacetin (R = H) or herbacetin diglucoside (HDG, R = *β*-d-glucose), (3) *p*-coumaric acid (R = H), or *p*-coumaric acid glucoside (CouG, R = *β*-d-glucose), (4) caffeic acid (R = H) or caffeic acid glucoside (CafG, R = *β*-d-glucose), (5) ferulic acid (R = H) or ferulic acid glucoside (FerG, R = *β*-d-glucose), (6) hydroxymethylglutaric acid (HMGA). (B) Schematic representation of lignan macromolecule, where a unit of SDG or HDG is ester-linked to another unit, thanks to HMGA, which can be replaced by one hydroxycinnamic acid glucoside (HCAG) unit (CouG, CafG, or FerG) in terminal position of the chain.</figcaption></figure>
The beneficial effects of lignans on human health are well recognized [5,6]. Particularly, the chemopreventive actions of SDG toward cancer, diabetes mellitus, and cardiovascular diseases have been largely described [5,7,8]. The pharmacological activity of this compound is thought to be due to its high antioxidant capacity [9–11] and to its phytoestrogenic activity [12]. Flavonols and hydroxycinnamic acids, the other constituents of the flaxseed lignan macromolecule, also display a wide range of health-promoting effects. The favorable actions on cardiovascular health of vegetable-rich diets have been ascribed to flavonols, and hydroxycinnamic acids have revealed powerful antioxidant properties and might be of particular interest for dermatologic applications [13,14].

Although both in vivo and in vitro data are globally in favor of a chemopreventive effect of lignans, epidemiological studies are much less conclusive, and the mechanism by which phytoestrogenic lignans prevent cancers still remains unclear [7] and requires further elucidation. This could be explained by the fact that our current knowledge concerning the genetic and environmental factors affecting productivity and yield stability of these phenolic compounds, in flaxseeds, remains partial, and little is known about the variation in antioxidant capacities of different flaxseed cultivars. Moreover, no study has put efforts toward linking the lignan content of different cultivars and the expression of genes involved in their biosynthetic pathway.

Herein we present a complete dataset concerning the relative impact of cultivar, edaphic, and climatic parameters on productivity of the main constituents of the lignan macromolecule of flaxseeds, in relation to their antioxidant capacities determined using both in vitro and in vivo systems. Such data could be useful to predict, more precisely, the accumulation and, therefore, the nutritional intakes of these compounds, with health benefits for pharmaceutical, nutraceutical, and/or cosmetic applications.

## 2. Materials and Methods

### 2.1. *Chemicals*

All chemicals were of analytical grade quality and purchased from Thermo (Illkirch, France). The deionized water was produced using a milli-Q water purification system (Merck Millipore, Molsheim, France). SDG and HDG standard were purchased from LGC Standards (Molsheim, France). The hydroxycinnamic acid glucosides—*p*-coumaric acid glucoside, caffeic acid glucoside, ferulic acid glucoside—were synthesized according to Beejmohun et al. (2004) [15] and Beejmohun et al. (2006) [16]. Prior to their use for HPLC or LC-MS analysis, all solutions were filtered through 0.45 µm nylon syringe membranes (Merck Millipore, Molsheim, France).

### 2.2. *Plant Materials and Cultivation*

Flax cultivars Astral, Baïkal, Baladin, Barbara, and Oliver were provided by Laboulet Semences (Airaines, France), Coopérative Linière Terre de Lin (Saint-Pierre-le-Viger, France) and Arvalis-Institut Technique du Lin (Boigneville, France). Flax was grown up to seed at the following locations in France: Eure (Gamaches-en-Vexin, GAM, 49°16'14''N/1°37'02''E/89 m), Somme (Airaines, AIR, 49°57'57''N/1°56'39''/70 m), and Eure-et-Loir (Chartres, CHA, 48°27'21.05''N/1°29'3.06''E/141 m). Sowings were performed on March 30th of each year, with 450 seeds per m². Fields were fertilized, immediately after sowing, with 80 units of nitrogen, 60 units of potassium, and 60 units of phosphorus per hectare (Figure S1). The soils of these sites were of clay loam type balanced, well-structured with a granulometry of ca. 25% 2000–63 µm, 50% 63–2 µm, and 25% <2 µm particles, and a pH around 7.8. The final harvest took place on August 15th of each year at the same ripening stage for each cultivar. Throughout the experiments, no visible disease or insect attack occurred at either location. During the growing period, the experimental stations received 182.8 mm (year 2003), 305.0 mm (year 2004), 269.8 mm (year 2005) for GAM, 387.4 mm for AIR (year 2005), and 308.6 mm for CHA (year 2005) of rainfall over the growing period. The day temperatures at an elevation of 2 m averaged 15.72 °C (year 2003), 13.67 °C (year 2004), 13.96 °C (year 2005) for GAM, 13.48 °C for AIR (year 2005),
and 14.45 °C for CHA (year 2005) over the growing period. All these meteorological characteristics are displayed in Table S1 and Figure S2.

## 2.3. *Gene Expression Analysis by RT-qPCR*

Total RNA was extracted from 100 mg of frozen plant material in liquid nitrogen as described by Hano et al. (2006) [17]. Expression patterns of *LuDIR5*, *LuPLR1*, and *LuUGT74S1* were analyzed using RT-qPCR, using specific primers described by Dalisay et al. (2015) [18]. For reverse transcription, 50 ng of total RNA was incubated for 60 min at 50 °C with 1× RT buffer, 0.5 mM of each dNTP, 1 µM of oligo-dT primers, 1 unit of RiboLock, and 4 units of Omniscript Reverse Transcriptase in a total volume of 20 µL (Qiagen, Hilden, Germany). qPCR was performed with a PikoReal™ Real-Time PCR System (Thermo Fisher Scientific, Villebon-sur-Yvette, France) using DyNAmo ColorFlash SYBR Green qPCR (ThermoScientific) and specific primers. Two reference genes (CYC and *ETIF5A*) were used for data normalization [19]. The qPCR parameters were as follows: an initial denaturation at 95 °C for 5 min, then 40 three-step cycles of 94 °C for 10 s, primer annealing at 65 °C for 10 s, and extension at 72 °C for 30 s. After 40 cycles, an additional extension step was performed at 72 °C for 90 s. The presence of a single amplicon was confirmed by the observation of a single peak in the melting curve obtained after amplification. Expression levels were calculated and normalized using 2^−ΔΔCt^ method [20]. Reactions were performed in three biological and two technical replicates.

## 2.4. *Extraction, HPLC, and LC-ESI-MS Analysis*

Extractions (4 biological and 2 technical replicates), quantification of compounds was carried out on a Varian liquid chromatographic system (Agilent Technology, Les Ulis, France), as well as LC-ESI-MS analyses using a Waters 2695 Alliance coupled with a single quadrupole mass spectrometer ZQ (Waters-Micromass, Manchester, UK), equipped with an electrospray ion source (ESI-MS), were performed as described in Corbin et al. (2015) [21].

## 2.5. *Determination of the Ferric-Reducing Antioxidant Power (FRAP)*

Ferric-reducing antioxidant power (FRAP) was measured as described by Benzie & Strain, (1996) [22] with little modification. Briefly, 10 µL of the extracted sample was mixed with 190 µL of FRAP (10 mM TPTZ; 20 mM FeCl₃·6H₂O, and 300 mM acetate buffer pH 3.6; ratio 1:1:10 (v/v/v)). Incubation lasted 15 min at room temperature. Absorbance of the reaction mixture was measured at 630 nm with a BioTek ELX800 Absorbance Microplate Reader (Thermo Fisher Scientific, Villebon-sur-Yvette, France). Assays were made in triplicate and antioxidant capacity was expressed as Trolox C equivalent antioxidant capacity (TAEC).

## 2.6. *Determination of Oxygen Radical Absorbance Capacity (ORAC)*

Oxygen radical absorbance capacity (ORAC) assay was performed as described by Prior et al. (2003) [23]. Briefly, 10 µL of the extracted sample was mixed with 190 µL of fluorescein (0.96 µM) in 75 mM phosphate buffer pH 7.4, and incubated for at least 20 minutes at 37 °C with intermittent shaking. Then, 20 µL of 119.4 mM 2,2'-azobis-amidinopropane (ABAP, Sigma Aldrich, Saint-Quentin Fallavier, France) was added and the fluorescence intensity was measured every 5 min for 2.5 h at 37 °C using a fluorescence spectrophotometer (Bio-Rad, Marnes-la-Coquette, France) set with an excitation at 485 nm and emission at 535 nm. Assays were made in triplicate, and antioxidant capacity was expressed as Trolox C equivalent antioxidant capacity (TAEC).

## 2.7. *Determination of the Iron-Chelating Capacity*

The iron-chelating capacity was determined as described by Mladenka et al. (2011) [24]. Briefly, 10 µL of extract sample were mixed with ferrous iron at a final concentration of 50 µM in HEPES (pH 6.8) buffer and 50 µL ferrozine (5 mM aqueous solution). All experiments were performed in
96-well microplates. Each sample was measured with and without (blank) the addition of ferrozine. Absorbance was measured at 550 nm immediately after addition of ferrozine, and 5 min later with a BioTek ELX800 Absorbance Microplate Reader (Thermo Fisher Scientific, Villebon-sur-Yvette, France). Chelating activity values were expressed in µM of fixed iron.

**2.8. Yeast Cells Cultivation and Treatments**

Yeast (Saccharomyces cerevisiae) strain MAV203 (Invitrogen, Thermo Fisher Scientific Villebon-sur-Yvette, France) were used. Cells were grown aerobically at 30 °C in an orbital shaker (150 rpm) in complete 2.0% (w/v) glucose YPD medium (Sigma Aldrich, Saint-Quentin Fallavier, France). All extracts evaporated under nitrogen flow, dissolved in DMSO at 50 µg/mL, and added to the cells 6 h before oxidative stress induction at a final concentration of 1 mg/mL. The final concentration of DMSO applied on the cell was 1 % (v/v). For the control sample, DMSO to 0.1% of the final volume, was added. Cells were irradiated with 106.5 J/m² UV-C (254 nm) under a Vilber VL-6.C filtered lamp (Thermo Fisher Scientific, Villebon-sur-Yvette, France), as described by Bisquert et al. (2018) [25], and then incubated overnight at 30 °C before membrane lipid peroxidation determination.

**2.9. Determination of Membrane Lipid Peroxidation Using Thiobarbituric Acid-Reactive Substances (TBARS) Assay**

Measurement of membrane lipid peroxide was carried out with the thiobarbituric acid (TBA; Sigma Aldrich, Saint-Quentin Fallavier, France) method described by Hano et al. (2008) [26]. Briefly, ca. 10⁷ cells were ground using a mortar and pestle in distilled water, and centrifuged at 10,000×g for 10 min. Supernatant fractions (75 µL) were mixed with 25 µL of 3% (w/v) SDS, 50 µL of 3% TBA (w/v) in 50 mM NaOH, and 50 µL of 23% (v/v) of HCl throughout mixing between each addition. The mixture was heated at 80 °C for 20 min. After cooling on ice, the absorbance at 532 nm (A532) was measured, and non-specific absorbance at 600 nm (A600) was subtracted.

**2.10. Statistical Treatment of Data**

All data presented in this study are the means and the standard deviations of at least three independent replicates. ANOVAs and Pearson correlations were performed using R software version 3.0.2. PCA was performed with XL-STAT2017 software (Addinsoft, Paris, France), with each parameter considered as a discrete variable; the initial dataset was then converted into principal components (PCs), and it was possible to graphically display the relationships among the considered parameters. Gene expression and SDG content were represented using MeV4 software. All statistical tests were considered significant at *p* < 0.05.

**3. Results and Discussion**

**3.1. Influence of Genetic Variations on the Accumulation of the Main Constituents of the Lignan Macromolecule**

The flax cultivars, herein studied, showed an SDG content ranging from 8.23 to 21.85 mg/g of dry weight (DW) (Table 1). Barbara and Oliver are high SDG-producing cultivars, Baladin presents an intermediate content, whereas Astral and Baïkal are poor in SDG, as compared to the other cultivars. A similar range of variation in SDG content has been reported in a flax germplasm collection by Diederichsen and Fu (2008) [27]. Lower SDG content was reported by Zimmermann et al. (2007, 2006) [28,29] for cultivars grown in Spain and Germany. Nevertheless, it should be noted that these authors employed an extraction method based on acid hydrolysis, which is known to be potentially destructive for SDG [30], leading to a possible underestimation in the actual contents. SDG is the main component of the lignan macromolecule accumulated in flaxseed, but other compounds, such as hydroxycinnamic acid glucosides (caffeic acid glucoside (CafG), *p*-coumaric acid glucoside (CouG), and ferulic acid glucoside (FerG), as well as the flavonol herbacetin diglucoside (HDG), are also incorporated in substantial amounts in this macromolecule [4,31,32]. Here, the whole set of these
compounds was assayed. In our hands, CouG contents ranged from 4.78 to 10.48 mg/g DW, and FerG content from 1.03 to 2.28 mg/g DW (Table 1). These results sound consistent with those described by Westcott and Muir (1996), Jonhson et al. (2000), and Eliasson et al. (2003) [33–35] for cultivars grown respectively in Canada, Denmark, and Sweden. To date, only semi-quantitative evaluation of the HDG variations in flax cultivars have been studied through NMR [36], therefore, to the best of our knowledge, the present work is the first study focusing on the quantitative variations in HDG contents in linseed cultivars. Concerning the quantitative variations in CafG, only Wang et al. (2017) [37] reported very low contents ranging from 2.40 to 8.70 µg/g DW for Chinese cultivars. Here, HDG content ranged from 0.75 to 1.18 mg/g DW, and CafG contents from 0.80 to 1.90 mg/g DW (Table 1).

<table><caption>**Table 1.** Influence of the cultivar (C), cultivation site (L), and year (Y) on the accumulation of the main constituents of the lignan macromolecule in flaxseeds.</caption><thead><tr><th>Cultivar</th><th>Location</th><th>Year</th><th>SDG <sup>a</sup></th><th>HDG <sup>a</sup></th><th>FerG <sup>a</sup></th><th>CouG <sup>a</sup></th><th>CafG <sup>a</sup></th></tr></thead><tbody><tr><td rowspan="4">Astral</td><td>AIR_05</td><td>12.85 ± 0.14</td><td>0.98 ± 0.06</td><td>1.65 ± 0.06</td><td>5.88 ± 0.09</td><td>0.85 ± 0.04</td></tr><tr><td>CHA_05</td><td>12.53 ± 0.11</td><td>1.05 ± 0.03</td><td>1.90 ± 0.08</td><td>6.05 ± 0.08</td><td>0.98 ± 0.06</td></tr><tr><td>GAM_03</td><td>11.73 ± 0.06</td><td>0.88 ± 0.04</td><td>1.58 ± 0.07</td><td>4.80 ± 0.08</td><td>1.23 ± 0.06</td></tr><tr><td>GAM_04</td><td>11.68 ± 0.09</td><td>1.10 ± 0.05</td><td>1.58 ± 0.04</td><td>4.78 ± 0.05</td><td>1.25 ± 0.05</td></tr><tr><td rowspan="4">Barbara</td><td>GAM_05</td><td>13.48 ± 0.13</td><td>0.85 ± 0.01</td><td>1.57 ± 0.02</td><td>6.07 ± 0.22</td><td>0.87 ± 0.02</td></tr><tr><td>AIR_05</td><td>21.68 ± 0.17</td><td>0.93 ± 0.06</td><td>1.95 ± 0.03</td><td>10.48 ± 0.12</td><td>1.83 ± 0.04</td></tr><tr><td>CHA_05</td><td>20.88 ± 0.07</td><td>1.05 ± 0.03</td><td>2.18 ± 0.07</td><td>8.63 ± 0.65</td><td>1.63 ± 0.04</td></tr><tr><td>GAM_03</td><td>20.63 ± 0.10</td><td>0.95 ± 0.03</td><td>1.78 ± 0.03</td><td>9.95 ± 0.11</td><td>1.33 ± 0.06</td></tr><tr><td rowspan="4">Baladin</td><td>GAM_04</td><td>21.85 ± 0.34</td><td>0.85 ± 0.03</td><td>1.78 ± 0.02</td><td>9.85 ± 0.26</td><td>1.53 ± 0.06</td></tr><tr><td>GAM_05</td><td>21.85 ± 0.81</td><td>0.88 ± 0.01</td><td>1.66 ± 0.01</td><td>9.84 ± 0.07</td><td>1.53 ± 0.01</td></tr><tr><td>AIR_05</td><td>16.03 ± 0.11</td><td>1.15 ± 0.03</td><td>2.10 ± 0.05</td><td>7.88 ± 0.07</td><td>1.48 ± 0.07</td></tr><tr><td>CHA_05</td><td>15.68 ± 0.19</td><td>1.18 ± 0.03</td><td>2.28 ± 0.07</td><td>7.58 ± 0.10</td><td>1.43 ± 0.09</td></tr><tr><td rowspan="4">Baikal</td><td>GAM_03</td><td>18.20 ± 0.38</td><td>1.03 ± 0.03</td><td>2.03 ± 0.04</td><td>7.88 ± 0.12</td><td>1.33 ± 0.02</td></tr><tr><td>GAM_04</td><td>16.45 ± 0.13</td><td>1.08 ± 0.07</td><td>2.08 ± 0.02</td><td>7.33 ± 0.21</td><td>1.40 ± 0.08</td></tr><tr><td>GAM_05</td><td>16.20 ± 0.38</td><td>0.95 ± 0.01</td><td>1.92 ± 0.01</td><td>6.91 ± 0.04</td><td>1.19 ± 0.01</td></tr><tr><td>AIR_05</td><td>8.33 ± 0.13</td><td>1.03 ± 0.04</td><td>1.75 ± 0.07</td><td>4.85 ± 0.07</td><td>0.85 ± 0.01</td></tr><tr><td rowspan="4">Oliver</td><td>CHA_05</td><td>8.23 ± 0.07</td><td>1.18 ± 0.02</td><td>1.75 ± 0.05</td><td>4.85 ± 0.06</td><td>1.05 ± 0.03</td></tr><tr><td>GAM_03</td><td>8.23 ± 0.17</td><td>1.15 ± 0.05</td><td>1.88 ± 0.07</td><td>4.90 ± 0.07</td><td>0.80 ± 0.05</td></tr><tr><td>GAM_04</td><td>8.45 ± 0.10</td><td>1.05 ± 0.08</td><td>1.98 ± 0.03</td><td>4.98 ± 0.08</td><td>0.98 ± 0.04</td></tr><tr><td>GAM_05</td><td>8.40 ± 0.04</td><td>0.98 ± 0.01</td><td>1.79 ± 0.02</td><td>4.88 ± 0.05</td><td>0.83 ± 0.01</td></tr><tr><td rowspan="4">F values</td><td>AIR_05</td><td>21.00 ± 0.25</td><td>0.75 ± 0.03</td><td>1.55 ± 0.03</td><td>10.15 ± 0.09</td><td>1.90 ± 0.05</td></tr><tr><td>CHA_05</td><td>19.50 ± 0.19</td><td>0.85 ± 0.03</td><td>1.33 ± 0.06</td><td>9.03 ± 0.26</td><td>1.75 ± 0.06</td></tr><tr><td>GAM_03</td><td>19.95 ± 0.11</td><td>0.98 ± 0.04</td><td>1.18 ± 0.06</td><td>9.38 ± 0.07</td><td>1.33 ± 0.06</td></tr><tr><td>GAM_04</td><td>19.93 ± 0.11</td><td>1.03 ± 0.07</td><td>1.03 ± 0.03</td><td>9.35 ± 0.09</td><td>1.45 ± 0.05</td></tr><tr><td></td><td>GAM_05</td><td>21.55 ± 0.37</td><td>0.83 ± 0.01</td><td>1.06 ± 0.01</td><td>9.12 ± 0.08</td><td>1.37 ± 0.01</td></tr><tr><td colspan="2">Cultivar (C)</td><td>284.62 ***</td><td>4.06 *</td><td>19.18 ***</td><td>95.33 ***</td><td>14.76 ***</td></tr><tr><td colspan="2">Location (L)</td><td>0.02</td><td>1.21</td><td>1.04</td><td>0.13</td><td>0.61</td></tr><tr><td colspan="2">Year (Y)</td><td>0.03</td><td>2.09</td><td>0.68</td><td>0.02</td><td>0.48</td></tr><tr><td colspan="2">C × L</td><td>194.69 ***</td><td>3.57 *</td><td>21.90 ***</td><td>75.96 ***</td><td>11.94 ***</td></tr><tr><td colspan="2">C × Y</td><td>198.99 ***</td><td>4.58 *</td><td>17.25 ***</td><td>59.81 ***</td><td>11.19 ***</td></tr><tr><td colspan="2">Y × L</td><td>0.02</td><td>1.484</td><td>0.536</td><td>0.062</td><td>0.441</td></tr><tr><td colspan="2">C × L × Y</td><td>148.70 ***</td><td>4.23 *</td><td>16.18 ***</td><td>51.26 ***</td><td>9.62 ***</td></tr></tbody></table>

In flaxseed, the lignan biosynthesis involves the dirigent protein (LuDIR5; Figure 2A)-mediated stereoselective coupling of two E-coniferyl alcohol moieties, resulting in the formation of (−)-pinoresinol [18,38]. The two following reaction steps leading to the conversion of (−)-pinoresinol to (−)-lariciresinol, and (−)-lariciresinol to (+)-secoisolariciresinol, are catalyzed by the same bifunctional enzyme pinoresinol-lariciresinol reductase (LuPLR1, Figure 2A) [17,39,40]. Secoisolariciresinol is then glycosylated into SDG under the control of UDP-glycosyltransferase (LuUGT74S1, Figure 2A) glycosylating the C-9 and C-9' hydroxyl positions [41,42]. SDG is stored as a

<footnote><sup>a</sup> All contents are given in mg/g DW. Values are mean ± SD of 4 independent replicates. ANOVA, F represents the effect. Significance level: * p < 0.05; ** p < 0.01; *** p < 0.001.</footnote>
3-hydroxy-3-methylglutaryl ester-linked complex (HMG-SDG), as shown in Figure 1. Formation of the HMG-SDG ester-linked oligomers, occurs by linking hydroxylmethylglutaryl (HMG) to C-6a and C-6a' position, via action of HMG CoA-transferase [43].

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Expression profile of flax lignan biosynthetic gene and SDG accumulation in five flaxseed cultivars. (A) Biosynthetic pathway leading to the formation of (+)-SDG in flaxseed. (B) Expression of *LuDIR5*, *LuPLR1*, and *LuUGT74S1* determined by RT-qPCR (normalized with *CYC* and *ETIF5A* reference genes) visualized using MeV4 (*n* = 3) and SDG content measured by HPLC and visualized using MeV4 (*n* = 3). (C) Pearson correlation matrix between (+)-SDG accumulation and the corresponding biosynthetic gene expression. Significance level: \* *p* < 0.05; \*\* *p* < 0.01; \*\*\* *p* < 0.001.</figcaption></figure>

Correlation analysis between the different constituents of the flax lignan macromolecule revealed significant positive correlations between the CafG, CouG, and SDG contents, on the one hand, and between HDG and FerG, on the other hand (Table 2). This correlation was in agreement with our previous results [36]. On the contrary, significant negative correlations were noted between the SDG vs HDG and FerG yields (Table 2), which confirmed our previous observations [36]. From a metabolic point of view, *p*-coumaric acid (Figure 1A) is a branch point leading to the biosynthesis of either flavonoids or lignans [44]. Therefore, caffeic acid and *p*-coumaric acid (Figure 1A) could be considered as more direct precursors for the HDG biosynthesis, whereas ferulic acid (Figure 1A) constitutes a precursor for SDG biosynthesis. These biosynthetic links could explain, in part, the observed correlations. Studying the possible metabolic channel regulation of the carbon allocation between these two branches, during flaxseed development, could be of particular interest.

As a step forward, lignan biosynthetic gene expression analysis performed on immature flaxseed (developmental stage 2; [17]) by RT-qPCR using the 3 specific genes involved in SDG biosynthesis (*LuDIR5*, *LuPLR1*, and *LuUGT74S1*; Figure 2A,B) appeared in good agreement with the HPLC quantification (Figure 2C). High expression of *LuPLR1* was detected in high SDG-producing cultivars, Barbara and Oliver, whereas Astral and Baikal cultivars, accumulating lower SDG content, showed a lower expression of these biosynthetic genes (Table 1). The steady state levels of the key *LuPLR1* transcripts [40], and the two other biosynthetic genes (*LuDIR5* and *LuUGT74S1*) are correlated with the SDG content measured in the corresponding mature seeds (Figure 2C), confirming the great
influence of genetic parameters (i.e., the cultivar), and indicated that most of the regulation occurred at transcriptional level.

<table><caption>**Table 2.** Correlation analysis using Pearson correlation coefficient (PCC).</caption><thead><tr><th>Variables</th><th>SDG</th><th>HDG</th><th>FerG</th><th>CouG</th><th>CafG</th><th>FRAP</th><th>ORAC</th><th>Iron<br>Chelation<br>inhibition</th></tr></thead><tbody><tr><td>SDG</td><td>-0.515 **</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HDG</td><td>-0.231 ns</td><td>0.510 **</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>FerG</td><td>0.966 ***</td><td>-0.476 *</td><td>-0.215 ns</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CouG</td><td>0.835 ***</td><td>-0.360 ns</td><td>-0.063 ns</td><td>0.832 ***</td><td></td><td></td><td></td><td></td></tr><tr><td>CafG</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>FRAP</td><td>0.676 ***</td><td>-0.624 ***</td><td>-0.545 **</td><td>0.639 ***</td><td>0.671 ***</td><td></td><td></td><td></td></tr><tr><td>ORAC</td><td>0.669 ***</td><td>-0.325 ns</td><td>-0.170 ns</td><td>0.676 ***</td><td>0.717 ***</td><td>0.573 **</td><td></td><td></td></tr><tr><td>Iron<br>Chelation</td><td>0.758 ***</td><td>-0.627 ***</td><td>-0.692 ***</td><td>0.768 ***</td><td>0.661 ***</td><td>0.817 ***</td><td>0.665 ***</td><td></td></tr><tr><td>MDA<br>inhibition</td><td>0.867 ***</td><td>-0.666 ***</td><td>-0.482 *</td><td>0.806 ***</td><td>0.721 ***</td><td>0.774 ***</td><td>0.617 ***</td><td>0.875 ***</td></tr></tbody></table>

Significance level: * p < 0.05; ** p < 0.01; *** p < 0.001; ns: not significant.

## 3.2. Influence of Geographic Parameters on the Accumulation of the Main Constituents of Lignan Macromolecule

It is well accepted that environmental conditions, such as the climate of the culture year and the location (soil conditions), could also greatly affect the accumulation of phenolic compounds, as previously observed by Oomah et al. (1996) [45] for the accumulation of total flavonoids in flaxseeds. Here, three different locations have been selected to provide access to the potential influence of edaphic condition on lignan accumulation in flaxseed. Bordered by four different seas, three mountain ranges, and the edge of the central European lowlands, France is known to be a country with very diverse climatic conditions, resulting in very different weather patterns. Here, the three selected experimental sites are representative of the major flax-growing areas in France, i.e., the western part, and the present contrasting weather patterns. The CHA site is characterized by the highest temperatures and the lowest rainfall during the seed maturation phase. On the contrary, AIR location presents the lowest temperatures and the highest rainfall observed during the same period. The last location, GAM, is considered as an intermediate in terms of climate. The impact of these different conditions, on the composition and amount of the main constituents of the lignan macromolecule accumulated in the seed of the five selected cultivars, are presented in Table 1. Analysis of the variance revealed that cultivar was the main contributor for the observed variability (cultivars, C, Table 1). Edaphic factor (location L, Table 1) has no significant effect on the accumulation of these phytochemicals, whereas significant interactions with genetic factor were noted, but evidenced the prominent effect of genetic background at a particular location according to *F* values (Table 1).

Nonetheless, the location constitutes a complex variable, differing by both climatic and edaphic parameters, thus, to evaluate the sole contribution of climate, we decided to compare flaxseeds grown at the same site, GAM (i.e., the same edaphic parameters) but in different cultivation years (i.e., different climatic parameters). Here, we chose to consider three consecutive years with very contrasting weather patterns and, for this reason, the 2003–2005 period was selected. Indeed, it must be noted that the summer of 2003 was the hottest and driest in recent decades, and must be regarded as extremely unusual. The 2003–2005 period was also the warmest period recorded in France since 1950, whereas the low rainfall observed from June 2004 to December 2005 led to a dramatic soil water deficit for 2005, with a soil humidity index close to 0.25 for GAM region (considering that a soil humidity index of 1 is for water-saturated soil whereas 0 is for water-depleted soil; see Table S1, Figure S2 for complete meteorological condition descriptions). As flax is known to be a water-demanding crop during its flowering period (i.e., June), we therefore decided to evaluate how these climate changes, leading to water deficiency during this period, have affected the flaxseed metabolism. The results are reported in
Table 1, and the analysis of the variance evidenced the genetic background (cultivars, C, Table 1) as the sole significant factor influencing the SDG, FerG, CouG, and CafG content (Table 1). The climatic parameters considered here (cultivation year, Y, Table 1) did not influence the accumulation of any molecules in the analyzed cultivars, whereas significant interaction between genetic and climatic parameters (C × Y, Table 1) was noted, but with lower *F* values as compared to genetic parameter alone (C, Table 1), indicating that the main contributions have to be attributed to this latter parameter. Our results are in good agreement with the results of Saastamoinen et al. (2013) [46], who also reported a lower impact of the cultivation year compared to the cultivar parameter on SDG accumulation. On the contrary, Wescott et al. (2002) [47] reported that the cultivation year could also influence SDG yield. This apparently contradictory result can be due to the complexity of the climatic variable, that could also be influenced by the nature of the soil considered (edaphic parameters). The nature of the soil could greatly affect the influence of the drought period as its ability to retain water greatly relied on its composition and granulometry. Indeed, a high soil ability to retain water could alleviate the effect of temporary drought, and differences in this feature could explain such apparent discrepancies, moreover, the rainfall regime differs between Scandinavian and Canadian summers, making it more probable that drought occurs during the latter.

All these phytochemical profiles were subjected to principal component analysis. The resulting biplot representation accounts for 99.43% (F1 + F2) of the initial variability of the data (Figure 3). Discrimination occurs mainly in the first dimension, and SDG content was the main contributor for this F1 axis that explains 89.43% of the initial variability. The concentrations of hydroxycinnamic acid glucosides (particularly CouG) were the main contributors for the second dimension (F2 axis), accounting for only 10% of the initial variability (Figure 3). PCA showed a significant grouping of samples as a function of their SDG content. Using this analysis, the different cultivars could also be easily discriminated. This PCA confirmed the prominence of the genetic background over the environmental (edaphic and climatic) factors studied here.

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Correlation circle for principal component analysis. The SDG, HDG, CafG, CouG, and FerG contents for 5 cultivars (Astral, Baïkal, Baladin, Barbara, and Oliver) growing at 3 different locations (GAM, AIR, and CHA) and over 3 different years (03 (2003), 04 (2004), or 05 (2005)) were submitted for analysis by the PCA algorithm in Excel-XLSTAT software, using the Pearson correlation matrix (at a significance level of *p* < 0.05).</figcaption></figure>
### 3.3. Evaluation and Comparison of In Vitro and In Vivo Antioxidant Capacities

To evaluate the influence of genetic and epigenetic variables on the health benefit potential of these flaxseeds, the antioxidant capacity of the corresponding extracts was then evaluated using both in vitro and in vivo assays. On the basis of the chemical reaction involved, the major antioxidant capacity assays can be roughly divided into two categories: i) hydrogen atom transfer (HAT) reaction-based assay, such as ORAC assay, or ii) electron transfer (ET) reaction-based assay, such as FRAP assay (Table 3).

In our hands, the antioxidant capacity of our flaxseed extracts revealed by these two different assays ranged from 217.35 (Baladin, AIR_05) to 355.75 (Oliver, GAM_05) µM of Trolox C equivalent antioxidant capacity (TEAC) using FRAP assay, and from 269.97 (Baikal, GAM_04) to 375.76 (Oliver, GAM_05) µM TEAC using ORAC assay (Table 3). The antioxidant capacity of polyphenolic compounds, such as lignans, has been previously attributed to their capacity for HAT, from their OH groups to the free radicals [48]. However, the radical scavenging capacity of these extracts occurring through an ET-based mechanism cannot be excluded, according to the high antioxidant values calculated from the FRAP assay (Table 3). Here these two in vitro antioxidant assays were significantly correlated with the presence of SDG, CouG and CafG (Table 2).

Besides these two mechanisms involved in the scavenging of reactive oxygen species, transient metal ion chelation is also considered as an antioxidant mechanism, since the Fenton reaction, responsible for the hydroxyl radical formation and, subsequently, radical chain reaction propagation, could be inhibited through this chelating mechanism [49,50]. Here, we evidenced that flaxseed extracts displayed an efficient iron (Fe²⁺)-chelating activity, ranging from 7.18 µM (Baikal, GAM_03) to 14.72 µM (Oliver, GAM_04) of fixed iron (Table 3), that could also contribute to their antioxidant activity, largely described in the literature [9,10]. In good agreement with the recent rationalization of the iron-chelating capacity of SDG and its aglycone form secoisolariciresinol, high SDG quantities associated with elevated contents of CouG and CafG, appeared to significantly contribute to the development of a high iron-chelating capacity of the corresponding flaxseed extracts (Table 2).

It is necessary to emphasize that the assays described herein are strictly predictive results based on the chemical reaction in vitro, however, they not necessary bear a great similarity to biological systems. The validity of these data has to be, therefore, considered as limited to a strict chemical sense with context interpretation. For this reason, in order to better reflect the in vivo situation, the antioxidant activity of these extracts was further investigated for their capacity to inhibit membrane lipid peroxidation induced by UV-C in yeast cells. Yeast cells have been proven to be an excellent model to evaluate in vivo antioxidant capacity in a cellular oxidative stress context [51]. Indeed, baker's yeast (Saccharomyces cerevisiae) is an attractive and reliable model. This organism is a true eukaryote, and the mechanisms of defense and adaptation to oxidative stress are well understood [25,52]. The in vivo anti-lipoperoxidation activity (inhibition of malondialdehyde (MDA) formation), determined using the TBARS assay, ranged from 21.68% (Baikal, GAM_04) to 47.02% (Oliver, AIR_05) (Table 3). Interestingly, a strong and significant correlation was observed between this cellular antioxidant capacity and the SDG (PCC = 0.867), CouG (PCC = 0.806), and CafG (PCC = 0.721) contents (Table 2). However, we can note that, since the contents of SDG, CouG, and CafG are highly correlated, these parameters are not independent, and it is, therefore, difficult to definitely judge their respective contribution to this biological activity (cellular antioxidant capacity) by single correlation analysis. In yeast models, a similar protective effect against oxidative stress was previously observed on yeast cells treated with thiamine [52] and melatonin [25]. To the best of our knowledge, this the first time that this system is applied to characterize a flax extract. Our results are in agreement with those obtained using in vitro assays, and highlighted the great in vivo antioxidant potential of flaxseed extracts as already proposed by Wang et al. (2017) [37], using another cellular antioxidant assay in HepG2 cells.
<table><caption>**Table 3.** Influence of the cultivar (C), cultivation site (L), and year (Y) on the in vitro and in vivo antioxidant activities of flaxseed extracts.</caption><thead><tr><th>Cultivar</th><th>Location_Year</th><th>FRAP <sup>a</sup></th><th>ORAC <sup>a</sup></th><th>Iron Chelation <sup>b</sup></th><th>MDA<br>Inhibition <sup>c</sup></th></tr></thead><tbody><tr><td rowspan="4">Astral</td><td>AIR_05</td><td>252.55 ± 2.45</td><td>281.55 ± 11.16</td><td>9.57 ± 0.25</td><td>37.10 ± 0.28</td></tr><tr><td>CHA_05</td><td>264.41 ± 10.28</td><td>317.87 ± 8.00</td><td>10.11 ± 0.41</td><td>36.95 ± 0.43</td></tr><tr><td>GAM_03</td><td>322.81 ± 1.41</td><td>263.13 ± 11.91</td><td>9.66 ± 0.12</td><td>38.47 ± 0.75</td></tr><tr><td>GAM_04</td><td>276.68 ± 4.33</td><td>312.34 ± 1.12</td><td>9.31 ± 0.53</td><td>37.71 ± 1.18</td></tr><tr><td rowspan="4">Barbara</td><td>GAM_05</td><td>252.55 ± 5.23</td><td>263.92 ± 4.56</td><td>9.57 ± 0.22</td><td>38.17 ± 1.07</td></tr><tr><td>AIR_05</td><td>332.55 ± 4.57</td><td>339.45 ± 1.95</td><td>11.97 ± 0.19</td><td>45.95 ± 2.26</td></tr><tr><td>CHA_05</td><td>317.61 ± 9.33</td><td>341.29 ± 16.09</td><td>11.35 ± 0.72</td><td>44.58 ± 1.08</td></tr><tr><td>GAM_03</td><td>292.81 ± 5.84</td><td>329.45 ± 8.65</td><td>11.79 ± 0.06</td><td>43.05 ± 3.23</td></tr><tr><td rowspan="4">Baladin</td><td>GAM_04</td><td>296.41 ± 8.20</td><td>336.55 ± 2.14</td><td>12.32 ± 0.53</td><td>46.41 ± 0.43</td></tr><tr><td>GAM_05</td><td>331.48 ± 5.23</td><td>309.45 ± 3.72</td><td>10.99 ± 0.28</td><td>45.65 ± 1.94</td></tr><tr><td>AIR_05</td><td>217.35 ± 4.67</td><td>286.82 ± 2.79</td><td>8.69 ± 0.31</td><td>33.44 ± 2.16</td></tr><tr><td>CHA_05</td><td>240.55 ± 6.31</td><td>334.18 ± 5.95</td><td>8.24 ± 0.44</td><td>33.28 ± 0.97</td></tr><tr><td rowspan="4">Baïkal</td><td>GAM_03</td><td>254.01 ± 2.30</td><td>307.61 ± 4.18</td><td>8.16 ± 0.43</td><td>32.06 ± 0.64</td></tr><tr><td>GAM_04</td><td>288.41 ± 11.27</td><td>321.82 ± 2.60</td><td>8.87 ± 0.06</td><td>28.85 ± 0.43</td></tr><tr><td>GAM_05</td><td>255.08 ± 14.24</td><td>281.03 ± 2.70</td><td>7.80 ± 0.59</td><td>34.50 ± 2.16</td></tr><tr><td>AIR_05</td><td>262.01 ± 3.21</td><td>278.66 ± 5.76</td><td>8.42 ± 0.37</td><td>25.34 ± 1.18</td></tr><tr><td rowspan="4">Oliver</td><td>CHA_05</td><td>274.95 ± 3.49</td><td>281.03 ± 2.79</td><td>8.16 ± 0.40</td><td>24.89 ± 1.19</td></tr><tr><td>GAM_03</td><td>227.61 ± 1.41</td><td>305.24 ± 2.32</td><td>7.18 ± 0.56</td><td>23.66 ± 1.62</td></tr><tr><td>GAM_04</td><td>243.75 ± 5.42</td><td>269.97 ± 4.09</td><td>7.62 ± 0.55</td><td>21.68 ± 3.13</td></tr><tr><td>GAM_05</td><td>243.78 ± 1.27</td><td>286.29 ± 2.51</td><td>8.07 ± 0.12</td><td>26.11 ± 0.86</td></tr><tr><td rowspan="5">F values</td><td>AIR_05</td><td>350.68 ± 0.81</td><td>368.66 ± 21.39</td><td>14.45 ± 0.25</td><td>47.02 ± 2.16</td></tr><tr><td>CHA_05</td><td>306.28 ± 5.27</td><td>374.45 ± 6.42</td><td>14.54 ± 0.44</td><td>46.87 ± 1.19</td></tr><tr><td>GAM_03</td><td>302.41 ± 6.17</td><td>278.92 ± 5.39</td><td>13.74 ± 0.38</td><td>43.21 ± 0.75</td></tr><tr><td>GAM_04</td><td>355.75 ± 8.20</td><td>326.55 ± 2.70</td><td>14.72 ± 0.47</td><td>45.34 ± 2.37</td></tr><tr><td>GAM_05</td><td>343.08 ± 5.28</td><td>375.76 ± 7.91</td><td>14.10 ± 0.38</td><td>46.87 ± 1.40</td></tr><tr><td>Genetic (C)</td><td></td><td>11.91 ***</td><td>5.37 ***</td><td>188.00 ***</td><td>161.33 ***</td></tr><tr><td>Location (L)</td><td></td><td>0.02</td><td>1.05</td><td>0.04</td><td>0.03</td></tr><tr><td>Year (Y)</td><td></td><td>0.06</td><td>0.15</td><td>0.04</td><td>0.08</td></tr><tr><td>C × L</td><td></td><td>7.20 **</td><td>4.59 **</td><td>133.03 ***</td><td>105.31 ***</td></tr><tr><td>C × Y</td><td></td><td>7.34 **</td><td>3.40 *</td><td>133.16 ***</td><td>130.60 ***</td></tr><tr><td>Y × L</td><td></td><td>0.03</td><td>0.64</td><td>0.04</td><td>0.06</td></tr><tr><td>C × L × Y</td><td></td><td>4.91 *</td><td>3.41 *</td><td>114.80 ***</td><td>109.27 ***</td></tr></tbody></table>

4. Conclusions

During the last decades, flaxseeds have emerged as one of the key sources of antioxidant phytochemicals. Knowledge about the variation in the accumulation of these valuable constituents is, hence, of particular interest. This study constitutes the first work devoted to the influence of genetic, edaphic, and climatic parameters on the main compounds constituting the so-called lignan macromolecule of flaxseeds, and the antioxidant activities of the obtained extracts. Our results evidenced the predominant influence of genetic factors (cultivar) on the accumulation of the constituents of the lignan macromolecule in flaxseeds. The results of gene expression suggest a transcriptional regulation of this accumulation, knowledge of which would help to manipulate the phenolic contents of flax. Elucidating the complete transcription regulation of lignan biosynthesis in flax would, therefore, help to better control their accumulation. In our hands, other environmental parameters, such as geographic and climatic variables, did not result in significant changes in the lignan macromolecule accumulation. Both in vitro and in vivo antioxidant activity relied on SDG, CafG, and CouG accumulations. Future works using purified compounds will be conducted to further elucidate their respective contribution to the cellular antioxidant capacity observed with

<footnote><sup>a</sup> expressed in mM of Trolox C equivalent antioxidant capacity (TEAC); <sup>b</sup> expressed in µM of fixed Fe<sup>2+</sup>; <sup>c</sup> expressed in % inhibition of MDA formation relative to control cells; values are mean ± SD of 4 independent replicates. ANOVA, F represents the effect. Significance level: * p < 0.05; ** p < 0.01; *** p < 0.001.</footnote>
flaxseed extracts. Considering the health benefits of these compounds, the present study evidenced the importance of a better knowledge of the flax cultivar type that could greatly influence the phytochemical intakes and the associated biological activities. Therefore, we recommend that this crucial parameter be considered in epidemiological studies dealing with flaxseeds.

**Supplementary Materials:** The following are available online, Table S1: Meteorological characteristics of the cultivation site; Figure S1: Scheme describing cultivation conditions; Figure S2: Climatic data for the trial sites Airaines (AIR), Gamaches-en-Vexin (GAM) and Chartres (CHA) for the years 2003 (03), 2004 (04), and 2005 (05). Precipitations are expressed as cumulative monthly rainfall in mm, and temperatures are an average of daily temperature in °C.

**Author Contributions:** Conceptualization, C.H., F.M. and E.L.; Methodology, L.G., S.D., C.C., C.D., J.L.L., E.A.L.; Software, T.F.; Validation, C.H., E.A.L., S.R., E.L. and F.M.; Formal Analysis, C.H., D.T., J.D., B.H.A., B.M. and O.F.; Investigation L.G., S.D., and C.C.; Resources, C.H. and E.L.; Data Curation, C.H., E.L. and S.R.; Writing—Original Draft Preparation, C.H.; Writing—Review & Editing, C.H., O.F., D.T., S.R., E.L., B.H.A and F.M.; Visualization, L.G., S.D. and C.H.; Supervision, C.H., B.M. and E.L.; Project Administration, C.H.; Funding Acquisition, C.H., E.L. and B.M.

**Funding:** This research was supported by Cosmetosciences, a global training and research program dedicated to the cosmetic industry. Located in the heart of the Cosmetic Valley, this program led by University of Orleans is funded by the Region Centre-Val de Loire. This research was also supported by the Conseil Departemental d'Eure et Loir. BHA acknowledges research fellowship of Le Studium-Institute for Advanced Studies, Loire Valley, Orléans, France. DT gratefully acknowledges the support of French government via the French Embassy in Thailand in the form of Junior Research Fellowship Program 2018.

**Acknowledgments:** The authors wish to thank Laboulet Semences, Coopérative Linière Terre de Lin, Arvalis-Institut Technique du Lin, and Graines de Lin 28 for the donation of flaxseed used in this study and agronomical data on cultivars and access to experimental fields.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Nayak, B.; Liu, R.H.; Tang, J. Effect of Processing on Phenolic Antioxidants of Fruits, Vegetables, and Grains—A Review. *Crit. Rev. Food Sci. Nutr.* **2015**, *55*, 887–919. [CrossRef] [PubMed]

2. Oomah, B.D. Flaxseed as a functional food source. *J. Sci. Food Agric.* **2001**, *81*, 889–894. [CrossRef]

3. Kamal-Eldin, A.; Peerlkamp, N.; Johnsson, P.; Andersson, R.; Andersson, R.E.; Lundgren, L.N.; Åman, P. An oligomer from flaxseed composed of secoisolariciresinoldiglucoside and 3-hydroxy-3-methyl glutaric acid residues. *Phytochemistry* **2001**, *58*, 587–590. [CrossRef]

4. Struijs, K.; Vincken, J.-P.; Doeswijk, T.G.; Voragen, A.G.J.; Gruppen, H. The chain length of lignan macromolecule from flaxseed hulls is determined by the incorporation of coumaric acid glucosides and ferulic acid glucosides. *Phytochemistry* **2009**, *70*, 262–269. [CrossRef] [PubMed]

5. Westcott, N.; Muir, A. Flax seed lignan in disease prevention and health promotion. *Phytochem. Rev.* **2003**, *2*, 401–417. [CrossRef]

6. McCann, M.J.; Gill, C.I.R.; McGlynn, H.; Rowland, I.R. Role of Mammalian Lignans in the Prevention and Treatment of Prostate Cancer Mark. *Nutr. Cancer* **2005**, *52*, 1–14. [CrossRef] [PubMed]

7. Lainé, E.; Hano, C.; Lamblin, F.F. *Phytoestrogens: Lignans*; Knasmüller, S., DeMarini, D.M., Johnson, I., Gerhäuser, C., Eds.; WILEY-VCH: Weinheim, Germany, 2009; ISBN 9783527320585.

8. Hano, C.; Renouard, S.; Molinié, R.; Corbin, C.; Barakzoy, E.; Doussot, J.; Lamblin, F.; Lainé, E. Flaxseed (*Linum usitatissimum* L.) extract as well as (+)-secoisolariciresinol diglucoside and its mammalian derivatives are potent inhibitors of α-amylase activity. *Bioorg. Med. Chem. Lett.* **2013**, *23*, 3007–3012. [CrossRef] [PubMed]

9. Prasad, K. Hydroxyl radical-scavenging property of secoisolariciresinol diglucoside (SDG) isolated from flax-seed. *Mol. Cell. Biochem.* **1997**, *168*, 117–123. [CrossRef] [PubMed]

10. Kitts, D.D.; Yuan, Y.V.; Wijewickreme, A.N.; Thompson, L.U. Antioxidant activity of the flaxseed lignan secoisolariciresinol diglycoside and its mammalian lignan metabolites enterodiol and enterolactone. *Mol. Cell. Biochem.* **1999**, *202*, 91–100. [CrossRef] [PubMed]

11. Hano, C.; Corbin, C.; Drouet, S.; Quéro, A.; Rombaut, N.; Savoire, R.; Molinié, R.; Thomasset, B.; Mesnard, F.; Lainé, E. The lignan (+)-secoisolariciresinol extracted from flax hulls is an effective protectant of linseed oil and its emulsion against oxidative damage. *Eur. J. Lipid Sci. Technol.* **2017**, *119*. [CrossRef]
12. Adlercreutz, H.; Mousavi, Y.; Clark, J.; Höckerstedt, K.; Hämäläinen, E.; Wähälä, K.; Mäkela, T.; Hase, T. Dietary phytoestrogen and cancer: In vitro and In vivo studies. *J. Steroid Biochem. Mol. Biol.* **1992**, *41*, 8012–8020. [CrossRef]

13. Schoenrock, U.; Untiedt, S.; Kux, U.; Inoue, K. Application of Ferulic Acid Glucosides as Anti-Irritants in Cosmetic and Topical Dermatological Preparations. German Patent AN 1997:708582, 1997.

14. Kosinska, A.; Penkacik, K.; Wiczkowski, W.; Amarowicz, R. Presence of caffeic acid in flaxseed lignan macromolecule. *Plant Foods Hum. Nutr.* **2011**, *66*, 270–274. [CrossRef] [PubMed]

15. Beejmohun, V.; Grand, E.; Mesnard, F.; Fliniaux, M.A.; Kovensky, J. First synthesis of (1,2-13C2)-monolignol glucosides. *Tetrahedron Lett.* **2004**, *45*, 8745–8747. [CrossRef]

16. Beejmohun, V.; Grand, E.; Lesur, D.; Mesnard, F.; Fliniaux, M.A.; Kovensky, J. Synthesis and purification of [1,2-13C2]coniferin. *J. Label. Compd. Radiopharm.* **2006**, *49*, 463–470. [CrossRef]

17. Hano, C.; Martin, I.; Fliniaux, O.; Legrand, B.; Gutierrez, L.; Arroo, R.R.J.; Mesnard, F.; Lamblin, F.; Lainé, E. Pinoresinol-lariciresinol reductase gene expression and secoisolariciresinol diglucoside accumulation in developing flax (*Linum usitatissimum*) seeds. *Planta* **2006**, *224*, 1291–1301. [CrossRef] [PubMed]

18. Dalisay, D.S.; Kim, K.W.; Lee, C.; Yang, H.; Rübel, O.; Bowen, B.P.; Davin, L.B.; Lewis, N.G. Dirigent Protein-Mediated Lignan and Cyanogenic Glucoside Formation in Flax Seed: Integrated Omics and MALDI Mass Spectrometry Imaging. *J. Nat. Prod.* **2015**, *78*, 1231–1242. [CrossRef] [PubMed]

19. Huis, R.; Hawkins, S.; Neutelings, G. Selection of reference genes for quantitative gene expression normalization in flax (*Linum usitatissimum* L.). *BMC Plant Biol.* **2010**, *10*, 71. [CrossRef] [PubMed]

20. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. *Methods* **2001**, *25*, 402–408. [CrossRef] [PubMed]

21. Corbin, C.; Fidel, T.; Leclerc, E.A.; Barakzoy, E.; Sagot, N.; Falguières, A.; Renouard, S.; Blondeau, J.; Ferroud, C.; Doussot, J.; et al. Development and validation of an efficient ultrasound assisted extraction of phenolic compounds from flax (*Linum usitatissimum* L.) seeds. *Ultrason. Sonochem.* **2015**, *26*, 176–185. [CrossRef] [PubMed]

22. Benzie, I.; Strain, J. The ferric reducing ability of plasma (FRAP) as a measure of "antioxidant power": The FRAP assay. *Anal. Biochem.* **1996**, *239*, 70–76. [CrossRef] [PubMed]

23. Prior, R.L.; Hoang, H.; Gu, L.; Wu, X.; Bacchiocca, M.; Howard, L.; Hampsch-Woodill, M.; Huang, D.; Ou, B.; Jacob, R. Assays for hydrophilic and lipophilic antioxidant capacity (oxygen radical absorbance capacity (ORACFL)) of plasma and other biological and food samples. *J. Agric. Food Chem.* **2003**, *51*, 3273–3279. [CrossRef] [PubMed]

24. Mladěnka, P.; MacÁková, K.; Filipský, T.; Zatloukalová, L.; Jahodář, L.; Bovicelli, P.; Silvestri, I.P.; Hrdina, R.; Saso, L. In vitro analysis of iron chelating activity of flavonoids. *J. Inorg. Biochem.* **2011**, *105*, 693–701. [CrossRef] [PubMed]

25. Bisquert, R.; Muñiz-Calvo, S.; Guillamón, J.M. Protective role of intracellular Melatonin against oxidative stress and UV radiation in Saccharomyces cerevisiae. *Front. Microbiol.* **2018**, *9*, 1–11. [CrossRef] [PubMed]

26. Hano, C.; Addi, M.; Fliniaux, O.; Bensaddek, L.; Duverger, E.; Mesnard, F.; Lamblin, F.; Lainé, E. Molecular characterization of cell death induced by a compatible interaction between Fusarium oxysporum f. sp. linii and flax (*Linum usitatissimum*) cells. *Plant Physiol. Biochem.* **2008**, *46*, 590–600. [CrossRef] [PubMed]

27. Diederichsen, A.; Fu, Y.B. Flax Genetic Diversity as the Raw Material for Future Success. In Proceedings of the International Conference on Flax and Other Bast Plants, Saskatoon, SK, Canada, 21–23 July 2008; pp. 270–279.

28. Zimmermann, R.; Bauermann, U.; Morales, F. Effects of growing site and nitrogen fertilization on biomass production and lignan content of linseed (*Linum usitatissimum* L.). *J. Sci. Food Agric.* **2006**, *86*, 415–419. [CrossRef]

29. Zimmermann, R.; Bauermann, U.; Spedding, C. Effects of nitrogen fertilisation and two growing sites on biomass production and lignan content of linseed (*Linum usitatissimum* L.): Second year. *Acta Agron. Hungarica* **2007**, *55*, 173–181. [CrossRef]

30. Li, H.B.; Wong, C.C.; Cheng, K.W.; Chen, F. Antioxidant properties in vitro and total phenolic contents in methanol extracts from medicinal plants. *LWT-Food Sci. Technol.* **2008**, *41*, 385–390. [CrossRef]

31. Struijs, K.; Vincken, J.P.; Verhoef, R.; van Oostveen-van Casteren, W.H.M.; Voragen, A.G.J.; Gruppen, H. The flavonoid herbacetin diglucoside as a constituent of the lignan macromolecule from flaxseed hulls. *Phytochemistry* **2007**, *68*, 1227–1235. [CrossRef] [PubMed]
32. Struijs, K.; Vincken, J.P.; Verhoef, R.; Voragen, A.G.J.; Gruppen, H. Hydroxycinnamic acids are ester-linked directly to glucosyl moieties within the lignan macromolecule from flaxseed hulls. *Phytochemistry* **2008**, *69*, 1250–1260. [CrossRef] [PubMed]

33. Westcott, N.D.; Muir, A.D. Variation in the concentration of the flax seed lignan concentration with variety, location and year. In Proceedings of the 56th Flax Institute of the United States Conference, Fargo, ND, USA, 20–22 March 1996; pp. 77–80.

34. Johnsson, P.; Kamal-Eldin, A.; Lundgren, L.N.; Aman, P. HPLC method for analysis of secoisolariciresinol diglucoside in flaxseeds. *J. Agric. Food Chem.* **2000**, *48*, 5216–5219. [CrossRef] [PubMed]

35. Eliasson, C.; Kamal-Eldin, A.; Andersson, R.; Aman, P. High-performance liquid chromatographic analysis of secoisolariciresinol diglucoside and hydroxycinnamic acid glucosides in flaxseed by alkaline extraction. *J. Chromatogr. A* **2003**, *1012*, 151–159. [CrossRef]

36. Ramsay, A.; Fliniaux, O.; Fang, J.; Molinie, R.; Roscher, A.; Grand, E.; Guillot, X.; Kovensky, J.; Fliniaux, M.A.; Schneider, B.; et al. Development of an NMR metabolomics-based tool for selection of flaxseed varieties. *Metabolomics* **2014**, *10*, 1258–1267. [CrossRef]

37. Wang, H.; Wang, J.; Qiu, C.; Ye, Y.; Guo, X.; Chen, G.; Li, T.; Wang, Y.; Fu, X.; Liu, R.H. Comparison of phytochemical profiles and health benefits in fiber and oil flaxseeds (*Linum usitatissimum* L.). *Food Chem.* **2017**, *214*, 227–233. [CrossRef] [PubMed]

38. Corbin, C.; Drouet, S.; Markulin, L.; Auguin, D.; Lainé, É.; Davin, L.B.; Cort, J.R.; Lewis, N.G.; Hano, C. A genome-wide analysis of the flax (*Linum usitatissimum* L.) dirigent protein family: From gene identification and evolution to differential regulation. *Plant Mol. Biol.* **2018**, *97*, 73–101. [CrossRef] [PubMed]

39. von Heimendahl, C.B.I.; Schäfer, K.M.; Eklund, P.; Sjöholm, R.; Schmidt, T.J.; Fuss, E. Pinoresinol-lariciresinol reductases with different stereospecificity from Linum album and Linum usitatissimum. *Phytochemistry* **2005**, *66*, 1254–1263. [CrossRef] [PubMed]

40. Renouard, S.; Tribalatc, M.; Lamblin, F.; Mongelard, G.; Fliniaux, O.; Corbin, C.; Marosevic, D.; Pilard, S.; Demailly, H.; Gutierrez, L.; et al. RNAi-mediated pinoresinol lariciresinol reductase gene silencing in flax (*Linum usitatissimum* L.) seed coat: Consequences on lignans and neolignans accumulation. *J. Plant Physiol.* **2014**, *171*, 1372–1377. [CrossRef] [PubMed]

41. Ghose, K.; Selvaraj, K.; McCallum, J.; Kirby, C.W.; Sweeney-Nixon, M.; Cloutier, S.J.; Deyholos, M.; Datla, R.; Fofana, B. Identification and functional characterization of a flax UDP-glycosyltransferase glucosylating secoisolariciresinol (SECO) into secoisolariciresinol monoglucoside (SMG) and diglucoside (SDG). *BMC Plant Biol.* **2014**, *14*, 82. [CrossRef] [PubMed]

42. Fofana, B.; Ghose, K.; McCallum, J.; You, F.M.; Cloutier, S. UGT74S1 is the key player in controlling secoisolariciresinol diglucoside (SDG) formation in flax. *BMC Plant Biol.* **2017**, *17*, 1–13. [CrossRef] [PubMed]

43. Ford, J.D.; Huang, K.; Wang, H.; Davin, L.B.; Lewis, N.G. Biosynthetic Pathway to the Cancer Chemopreventive Secoisolariciresinol Diglucoside-Hydroxymethyl Glutaryl Ester-Linked Lignan Oligomers in Flax (*Linum usitatissimum*) Seed. *J. Nat. Prod.* **2001**, *2*, 1388–1397. [CrossRef]

44. Žuk, M.; Kulma, A.; Dymińska, L.; Szoltysek, K.; Prescha, A.; Hanuza, J.; Szopa, J. Flavonoid engineering of flax potentiate its biotechnological application. *BMC Biotechnol.* **2011**, *11*, 10. [CrossRef] [PubMed]

45. Dave Oomah, B.; Mazza, G.; Kenaschuk, E.O. Flavonoid content of flaxseed. Influence of cultivar and environment. *Euphytica* **1996**, *90*, 163–167. [CrossRef]

46. Saastamoinen, M.; Pihlava, J.M.; Eurola, M.; Klemola, A.; Jauhiainen, L.; Hietaniemi, V. Yield, SDG lignan, cadmium, lead, oil and protein contents of linseed (*Linum usitatissimum* L.) cultivated in trials and at different farm conditions in the south-western part of Finland. *Agric. Food Sci.* **2013**, *22*, 296–306. [CrossRef]

47. Westcott, N.D.; Muir, A.D.; Lafond, G.; McAndrew, D.W.; May, W.; Irvine, B.; Grant, C.; Shirliffe, S.; Bruulsema, T.W. Factors Affecting the Concentration of a Nutraceutical Lignan in Flaxseed. In Proceedings of the Symposium on Fertilizing Crops for Functional Food, Indianapolis, IN, USA, 11 November 2002; pp. 1–3.

48. Podloucká, P.; Berka, K.; Fabre, G.; Paloncýová, M.; Duroux, J.L.; Otyepka, M.; Trouillas, P. Lipid bilayer membrane affinity rationalizes inhibition of lipid peroxidation by a natural lignan antioxidant. *J. Phys. Chem. B* **2013**, *117*, 5043–5049. [CrossRef] [PubMed]

49. Donoso-Fierro, C.; Becerra, J.; Bustos-Concha, E.; Silva, M. Chelating and antioxidant activity of lignans from Chilean woods (Cupressaceae). *Holzforschung* **2009**, *63*, 559–563. [CrossRef]
50. Fucassi, F.; Heikal, A.; Mikhalovska, L.I.; Standen, G.; Allan, I.U.; Mikhalovsky, S.V.; Cragg, P.J. Metal chelation by a plant lignan, secoisolariciresinol diglucoside. *J. Incl. Phenom. Macrocycl. Chem.* **2014**, *80*, 345–351. [CrossRef]

51. Steels, E.L.; Learmonth, R.P.; Watson, K. Stress tolerance and membrane lipid unsaturation in *Saccharomyces cerevisiae* grown aerobically or anaerobically. *Microbiology* **1994**, *140*, 569–576. [CrossRef] [PubMed]

52. Wolak, N.; Kowalska, E.; Kozik, A.; Rapala-Kozik, M. Thiamine increases the resistance of baker's yeast *Saccharomyces cerevisiae* against oxidative, osmotic and thermal stress, through mechanisms partly independent of thiamine diphosphate-bound enzymes. *FEMS Yeast Res.* **2014**, *14*, 1249–1262. [CrossRef] [PubMed]

<figure><img src="image_5.png" /></figure>

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).