# High-density genetic map and identification of QTLs for responses to temperature and salinity stresses in the model brown alga *Ectocarpus*

Komlan Avia, Susana Coelho, Gabriel J Montecinos, Alexandre Cormier, Fiona Lerck, Stéphane Mauger, Sylvain Faugeron, Myriam Valero, J. Mark Cock, Pierre Boudry

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

## ► To cite this version:

Komlan Avia, Susana Coelho, Gabriel J Montecinos, Alexandre Cormier, Fiona Lerck, et al. High-density genetic map and identification of QTLs for responses to temperature and salinity stresses in the model brown alga *Ectocarpus*. Scientific Reports, 2017, 7 (1), pp.43241. 10.1038/srep43241. hal-03145860

HAL Id: hal-03145860

https://hal.inrae.fr/hal-03145860

Submitted on 18 Feb 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

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

OPEN

# High-density genetic map and identification of QTLs for responses to temperature and salinity stresses in the model brown alga *Ectocarpus*

Received: 05 September 2016

Accepted: 20 January 2017

Published: 03 March 2017

Komlan Avia<fnref n="1" />,², Susana M. Coelho<fnref n="1" />, Gabriel J. Montecinos<fnref n="1" />,²,³, Alexandre Cormier<fnref n="1" />, Fiona Lerck<fnref n="1" />, Stéphane Mauger², Sylvain Faugeron²,³, Myriam Valero², J. Mark Cock<fnref n="1" /> & Pierre Boudry⁴

Deciphering the genetic architecture of adaptation of brown algae to environmental stresses such as temperature and salinity is of evolutionary as well as of practical interest. The filamentous brown alga *Ectocarpus* sp. is a model for the brown algae and its genome has been sequenced. As sessile organisms, brown algae need to be capable of resisting the various abiotic stressors that act in the intertidal zone (e.g. osmotic pressure, temperature, salinity, UV radiation) and previous studies have shown that an important proportion of the expressed genes is regulated in response to hyposaline, hypersaline or oxidative stress conditions. Using the double digest RAD sequencing method, we constructed a dense genetic map with 3,588 SNP markers and identified 39 QTLs for growth-related traits and their plasticity under different temperature and salinity conditions (tolerance to high temperature and low salinity). GO enrichment tests within QTL intervals highlighted membrane transport processes such as ion transporters. Our study represents a significant step towards deciphering the genetic basis of adaptation of *Ectocarpus* sp. to stress conditions and provides a substantial resource to the increasing list of tools generated for the species.

Brown algae (Phaeophyceae) represent one of only five eukaryotic lineages that have evolved complex multicellularity independently. They are very distantly related to other major lineages such as animals or green plants<fnref n="1" />,². These seaweeds constitute the dominant vegetation in the intertidal and subtidal zone of coastal ecosystems. They display a large variability of morphologies and some brown algae, such as kelps, may form extensive forests that offer suitable habitats for numerous species and hence support ocean biodiversity³,⁴. Brown algae have a wide range of uses as food, cosmetics or fertilizers and there is increasing attention on new biotechnological applications such as biofuels or cell-wall polysaccharides⁵⁻⁷. Despite their evolutionary, ecological and economic importance, many aspects of the biology of brown algae remain poorly understood. In particular, the range of phenotypic variation of such morphologically simple organisms, as well as the sources of this variation, their impact on fitness and the limits they may impose on adaptive responses to environmental conditions are largely unknown.

A large number of seaweed species are able to colonize the intertidal zone, and some of these species live only in this habitat. The intertidal zone is often characterized by the presence of multiple stresses including high temperature, low or high salinity, extreme irradiation with sunlight, wave mechanical force, desiccation and osmotic shocks due to tidal emersion-immersion cycles. There is evidence that the vertical distribution of organisms is related to their ability to tolerate different types of stress. This is particularly the case in the intertidal zone, where the difference between tolerant and susceptible species lies in the extent to which they can recover photosynthesis

<footnote n="1">¹Algal Genetics Group, UMR 8227, CNRS, Sorbonne Universités, UPMC, Station Biologique Roscoff, CS 90074, 29688 Roscoff, France. ²UMI 3614 Evolutionary Biology and Ecology of Algae, CNRS, Sorbonne Universités, UPMC, Pontificia Universidad Católica de Chile, Universidad Austral de Chile, Station Biologique Roscoff, CS 90074, 29688 Roscoff, France. ³Centro de Conservación Marina and CeBiB, Facultad de Ciencias Biológicas, Pontificia Universidad Católica de Chile, Casilla 114-D, Santiago, Chile. ⁴Ifremer, Laboratoire des Sciences de l'Environnement Marin (UMR 6539 LEMAR, UBO, CNRS, IRD, Ifremer), Centre Bretagne - ZI de la Pointe du Diable, CS 10070, 29280 Plouzané, France. Correspondence and requests for materials should be addressed to K.A. (email: komlan.avia@sb-roscoff.fr)</footnote>
when re-immersed in seawater⁸. Coastal habitats may present other stresses, such as low salinity in estuarine environments, high concentrations of heavy metals and other chemical pollutants near ports, mines or coastal cities, spatial temporal variations in upwelling and other coastal oceanographic processes⁹,¹⁰. Adding to this complexity, variations in environmental conditions and their impact can even be observed at the microhabitat scale¹¹. Remarkably, a number of seaweed species are able to cope with these strong environmental gradients. From an evolutionary point of view, two mechanisms can explain this tolerance to abiotic stressors: either evolution towards a generalist, all-purpose phenotype by increasing phenotypic plasticity, or local adaptation to specific sets of environmental conditions. It is however difficult to determine the relative importance of phenotypic plasticity and local adaptation, respectively, in determining differences among populations¹². Plasticity may mimic local adaptation resulting in differences between populations that are not due to genetic factors¹³. While an extensive literature is available on the range of phenotypic plasticity in responses of seaweeds to both biotic and abiotic factors¹⁴, the genetic basis of their responses to environmental heterogeneity remains poorly explored. Several mechanisms involved in the responses to stressors such as heat, salinity and dehydration have been shown to be conserved among land plant taxa, including the presence of common actors such as reactive oxygen species (ROS), ion fluxes, activation of kinases and a cascade of reactions leading to the expression of transcription factors¹⁵⁻¹⁷. However, the extent to which mechanisms of stress tolerance known for those taxa are used by brown algae is largely unknown. Although some aspects might be conserved¹⁸,¹⁹, other important processes are likely to be found only in brown algae because of their unique evolutionary history. To understand the stress tolerance features characteristic of this group, it is therefore essential to study this phenomenon using brown algal models rather than extrapolating information from distant models in other phylogenetic groups. Such knowledge is of evolutionary as well as of applied interest, as novel biomolecules and metabolic pathways may be discovered. Furthermore, with the current expansion of the algal aquaculture sector, brown algae need to be well characterized for future domestication and selective breeding to provide the industry with optimized strains.

*Ectocarpus* sp., a small filamentous alga, has been established as the model organism for the brown algae, inducing the generation of several genetic and genomic resources such as a fully sequenced genome¹, a genetic map²⁰ and transcriptomic data¹⁹,²¹. Previously named *Ectocarpus siliculosus*, it has recently become clear that this taxa corresponds to a separate species that has not been described previously²² and we therefore refer to the species provisionally as *Ectocarpus* sp. 1c²³, simplified as *Ectocarpus* sp. in the text below. As with many other members of the brown algal group, this species alternates between a haploid (gametophyte) and a diploid (sporophyte) phase, and thus its life cycle involves sequential development of two successive complex organisms²⁴. A recent study in the *Ectocarpus* group has shown that haploid and diploid organisms also exhibit clear niche partitioning. However, mainly asexual populations also exhibit similar habitat usage²⁵. Although brown algae are mainly marine species, rare cases of adaptation to fresh water have been reported²⁶. An *Ectocarpus* species, *Ectocarpus subulatus*, is able to live in both marine and freshwater, suggesting extreme physiological plasticity²³,²⁷⁻²⁹. Freshwater colonization is a rare event that is thought to induce rapid evolutionary radiations³⁰. The colonization event involving this brown alga is still under investigation but it does not appear to have involved extensive radiation, as only about 1% of these species have colonized freshwater³¹. Recent analyses showed that the transition to low salinity has been accompanied by fundamental morphological, transcriptomic and metabolomic changes but might be still reversible³¹, suggesting high phenotypic plasticity rather than local adaptation. This transition is also dependent on the host interactions with its complex bacterial communities²⁹. Although *E. subulatus* is a different species, such host-bacterial interactions may also play a role in the adaptation of *Ectocarpus* sp. to its environment (although this aspect has not been investigated in the study presented here).

Analysis of the *Ectocarpus* sp. genome suggested that this species has evolved effective mechanisms to survive its harsh intertidal environment such as a complex photosynthetic system and a large number of ion channels¹. These characteristics may have contributed to its ability to colonize freshwater environments. However, even with the current improved assembly and annotation of the *Ectocarpus* genome, about 39% of the predicted genes lack functional annotations based on available databases³². Hence, an important proportion of the genes in the genome do not have a predicted function and, therefore, approaches such as detection of quantitative trait loci (QTLs) and analyses of polymorphism of candidate genes³³ are important to explore the genetic basis of adaptation in this species.

QTL analysis provides insights into the genetic determinism of traits and is an important aid for marker assisted selection³⁴,³⁵. QTL mapping has been used successfully to improve our understanding of the genetic basis of traits in land plants³⁶ and in animals, including marine animals³⁷, but remains rarely investigated in algae³⁸. With recent advances in sequencing technology, molecular markers such as single nucleotide polymorphisms (SNPs) can now be identified in abundance at continually decreasing cost, even in complex genomes. Methods such as Restriction site Associated DNA (RAD) sequencing allow the complexity of the genome to be reduced in order to produce high confidence SNPs cost effectively for various applications³⁹,⁴⁰.

In our study, we generated 3,588 RAD-seq-based SNP markers to construct a dense genetic map for *Ectocarpus* sp. and used it to identify QTLs for growth-related traits and their plasticity under different temperature and salinity conditions. Several candidate QTLs for high temperature and low salinity tolerance were detected and GO (Gene Ontology) enrichment tests within these QTLs contributed to our understanding of the ability of *Ectocarpus* sp. to adapt to high temperature or low salinity.

## Results

**Identification of phenotypic traits under temperature and salinity stress.** To study the genetic basis of growth under temperature and salinity stresses, we grew a family of 89 *Ectocarpus* sp. haploid progeny, together with their two parental strains (six clones for each strain) for 12 days in growth cabinets under three temperature (13 °C as control, 26 °C and 28 °C; hereafter T13, T26 and T28) and three salinity conditions (34% as control, 20% and 15%; hereafter S34, S20 and S15). Temperature and salinity stresses were applied
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Distributions of the progeny BLUPs of the studied traits. T13, T26 and T28 correspond to growth rate phenotypes for temperatures 13°C, 26°C and 28°C respectively. S34, S20 and S15 correspond to growth rate phenotypes for 34‰ salinity, 20‰ and 15‰ respectively. Plasticity traits are represented by a letter followed by a number and a slash, then another letter and number. For example, S34/S20 corresponds to plasticity represented by slope of the reaction norm between salinity 34‰ condition and salinity 20‰ condition. Surv_S15 and Surv_S20 correspond to survival trait for conditions salinity 15‰ and salinity 20‰. Frequency distributions of the progeny BLUPs of the analyzed traits; blue arrows indicate mean values for the male parent and red arrows indicate mean values for the female parent; the green dashed line represents the mean value for the progeny population; “Slope” represents the slope values of the reaction norm between two conditions indicated at the top of each graph; “% survival” represents the proportion of survival among replicates of each progeny for salinity conditions 15‰ and 20‰.</figcaption></figure>

under normal, control salinity (34‰) and temperature (13°C) conditions, respectively. An initial series of test experiments aimed at selecting optimal temperature and salinity conditions for QTL detection had shown that 30°C and 10‰ salinity were extreme conditions for the mapping family, leading to high mortality. We therefore selected our stress conditions so that they had significant impacts on the growth phenotype of the offspring without killing a high proportion of individuals (we considered up to 30% of mortality as acceptable, i.e. none of the conditions should induce a complete loss of more than 26 progeny).

Growth, expressed as a ratio of the thallus area on day 12 to the thallus area on day 0, varied significantly across the studied progeny under all three salinity conditions (p-value between $4 \times 10^{-4}$ and $7 \times 10^{-16}$, ANOVA) but only under two of the three temperature conditions (p-value = 0.005 for 26°C, p-value < 1 × $10^{-10}$ for 28°C and p = 0.2 for 13°C, ANOVA). A significant sex-dependent difference was observed under 15‰ salinity (adjusted p-value Tukey comparison test: $3 \times 10^{-7}$) with males showing more growth than females (growth ratio of 0.47 vs 0.30 respectively).

The impact of the applied stresses was clearly visible through their effects on the distribution of the traits, with a reduced growth rate for the stress conditions and a much lower population mean for the growth ratio (Fig. 1). We observed substantial variance in all traits across the progeny, with the fastest growers having growth ratios 2 to 5-fold larger than the slowest growers. Differences between the mean values for each parent were more limited.

Transgressive individuals (i.e. displaying more extreme phenotypes than the parental strains) were observed in the mapping family under almost all conditions (Fig. 1).

By treating different experimental conditions as different environments, we observed significant genotype by environment (G × E) interactions (genotype × salinity p-value = $2 \times 10^{-6}$ and genotype × temperature p-value = $4 \times 10^{-4}$), which is potentially indicative of a genetic basis of variation of phenotypic plasticity$^{41,42}$. Reaction norms of the progeny for temperature and salinity variations were also analyzed. The results
<table><caption>**Table 1. Trait statistics.** T13, T26 and T28 correspond to growth rate phenotypes for temperatures 13 °C, 26 °C and 28 °C respectively. S34, S20 and S15 correspond to growth rate phenotypes for 34‰ salinity, 20‰ and 15‰ respectively. Surv_S15 and Surv_S20 correspond to survival trait for conditions salinity 15‰ and salinity 20‰. Salinity, Temperature and Survival correspond to the joint analysis of the three salinity conditions, the three temperature conditions and the two survival conditions respectively. $V_g$ and $CV_g$ corresponds to genetic variances and their coefficients of variation ($CV_g$ expressed as $\sqrt{V_g/\text{trait mean}}$), $H^2$ (s.e.) corresponds to broad-sense heritabilities and their standard errors estimated via jackknife resampling. Min and Max are based on average values for individuals and not the raw replicate values.</caption><thead><tr><th>Trait</th><th>Mean</th><th>Min</th><th>Max</th><th>V<sub>g</sub></th><th>CV<sub>g</sub></th><th>H<sup>2</sup> (s.e.)</th></tr></thead><tbody><tr><td>S34</td><td>1.49</td><td>0.67</td><td>2.16</td><td>0.039</td><td>0.13</td><td>0.42 (0.11)</td></tr><tr><td>S20</td><td>0.64</td><td>0.024</td><td>1.3</td><td>0.054</td><td>0.36</td><td>0.78 (0.04)</td></tr><tr><td>S15</td><td>0.37</td><td>0.008</td><td>0.9</td><td>0.033</td><td>0.49</td><td>0.78 (0.05)</td></tr><tr><td>T13</td><td>1.56</td><td>1.17</td><td>2.42</td><td>0.014</td><td>0.08</td><td>0.20 (0.14)</td></tr><tr><td>T26</td><td>1.18</td><td>0.4</td><td>2.07</td><td>0.024</td><td>0.13</td><td>0.36 (0.14)</td></tr><tr><td>T28</td><td>0.34</td><td>0.03</td><td>0.86</td><td>0.019</td><td>0.41</td><td>0.62 (0.08)</td></tr><tr><td>Surv_S20</td><td>—</td><td>—</td><td>—</td><td>1.02</td><td>—</td><td>0.65 (0.09)</td></tr><tr><td>Surv_S15</td><td>—</td><td>—</td><td>—</td><td>1.73</td><td>—</td><td>0.76 (0.06)</td></tr><tr><td>Salinity</td><td>1</td><td>0.008</td><td>2.16</td><td>0.009</td><td>0.09</td><td>0.31 (0.19)</td></tr><tr><td>Temperature</td><td>1.03</td><td>0.03</td><td>2.42</td><td>0.002</td><td>0.04</td><td>0.1 (0.18)</td></tr><tr><td>Survival</td><td>—</td><td>—</td><td>—</td><td>1.42</td><td>—</td><td>0.72 (0.05)</td></tr></tbody></table>

showed that several reaction norms had different slopes and therefore crossed, suggesting significant $G \times E$ (Supplementary Fig. S1). We therefore used the slopes of these reaction norms as plasticity traits. There were three temperature plasticity traits (T13/T26, T13/T28 and T26/T28) and three salinity plasticity traits (S34/S20, S34/S15 and S20/S15, Fig. 1).

In total, 15 traits were tested: 3 for temperature, 3 for salinity, 3 for temperature plasticity, 3 for salinity plasticity, 2 for survival, plus sex (Fig. 1).

Broad-sense heritability estimates for growth related-traits and survival varied between 0.20 and 0.78 for individual traits and between 0.10 and 0.72 when combining all conditions of temperature and salinity (Table 1).

**Construction of a ddRAD-seq-based genetic map.** Barcoded ddRAD libraries were constructed for the two parents and their 89 progeny using the enzymes *Pst*1 and *Hha*1 and were pooled for sequencing. Of the initial total of 470 million reads, the Process_radtags module of the Stacks pipeline together with additional cleaning steps retained 339 million reads. The Ref_map process in Stacks generated 183,336 consensus sequences from this sequence data, of which 181,655 were unique, and the Genotypes program identified 11,740 putative markers within these regions. After application of a series of filters (see materials and methods), 6,275 of the 11,740 markers identified by Stacks were retained. During the filtering, nine individuals (out of the 89 progeny) with excessive missing genotypes were removed.

For the construction of the genetic map, distorted markers that did not show the expected 1:1 segregation pattern (at a threshold of 0.05%) were also filtered out. When two or more markers exhibited similar segregation patterns (based on a similarity of ≥0.95), only one of the markers was retained. Markers with suspect linkages were also filtered out. Finally, a total of 3,588 markers and 80 individuals were used for the construction of the genetic map. Maps obtained with R/qlt and JoinMap using this data were compared and found to be very similar.

The 3,588 SNPs were distributed over 28 linkage groups (LGs), with a minimum logarithm of the odd (LOD) score of 4.0. The total length of the map was 2,585.7 centimorgans (cM), which accounted for 98.35% of the estimated genome (Ge) length (Ge = 2629 cM). The average spacing between two adjacent markers was 0.7 cM and the largest gap was 20.3 cM (on LG28). The lengths of the 28 LGs ranged from 41.8 cM with 54 markers to 152.3 with 217 markers (Table 2, Supplementary Fig. S2 and Supplementary Dataset).

Considering the 214 Mbp genome size of *Ectocarpus* sp.¹, the estimated global recombination rate was 12.28 cM/Mb. The average recombination rate per LG varied between 11.11 cM/Mb for LG28 and 17.93 cM/Mb for LG4. At a 1 Mbp scale along the different LGs, the recombination rate varied considerably with broad patterns of decreases and increases (Supplementary Fig. S3). However, there were no substantial differences between the LGs.

**Identification of QTLs for responses to temperature and salinity stresses.** Combination of genotyping and phenotyping data allowed the identification of a total of 39 QTLs, with between 0 and 10 QTLs detected per trait across a total of 15 traits tested. The highest number of QTLs (10) was observed for the plasticity trait S34/S20. No QTLs were detected for the traits T13 and S34. Each QTL is associated among others with the observed phenotypic variance ($R^2$) and an estimated additive effect (a) that it accounts for. The detected QTLs explained 6% to 34% of the observed phenotypic variance, depending on the trait (Table 3). In some cases, QTLs for different traits mapped approximately to the same genomic region (Fig. 2). For all the temperature and salinity traits (i.e. growth ratios), and for the plasticity traits, there were both cases where the allele transmitted by the male parent was favorable (increasing the growth ratio) and cases where it was unfavorable (decreasing the
<table><caption>Table 2. Linkage group statistics for the genetic map.</caption><thead><tr><th>LG</th><th>Number of markers</th><th>LG length (cM)</th><th>Average spacing between markers (cM)</th><th>Maximum spacing between markers (cM)</th></tr></thead><tbody><tr><td>1</td><td>217</td><td>152.3</td><td>0.7</td><td>11.2</td></tr><tr><td>2</td><td>194</td><td>96.5</td><td>0.5</td><td>5.3</td></tr><tr><td>3</td><td>202</td><td>90.4</td><td>0.4</td><td>8.5</td></tr><tr><td>4</td><td>180</td><td>116.7</td><td>0.7</td><td>7.6</td></tr><tr><td>5</td><td>168</td><td>117.9</td><td>0.7</td><td>6.7</td></tr><tr><td>6</td><td>171</td><td>110.4</td><td>0.6</td><td>8.8</td></tr><tr><td>7</td><td>170</td><td>113.4</td><td>0.7</td><td>9.3</td></tr><tr><td>8</td><td>152</td><td>88.7</td><td>0.6</td><td>6.3</td></tr><tr><td>9</td><td>155</td><td>87.8</td><td>0.6</td><td>7.6</td></tr><tr><td>10</td><td>138</td><td>92</td><td>0.7</td><td>8</td></tr><tr><td>11</td><td>120</td><td>92</td><td>0.8</td><td>10</td></tr><tr><td>12</td><td>127</td><td>112.4</td><td>0.9</td><td>8.9</td></tr><tr><td>13</td><td>124</td><td>92.6</td><td>0.8</td><td>7.6</td></tr><tr><td>14</td><td>124</td><td>96.7</td><td>0.8</td><td>6.3</td></tr><tr><td>15</td><td>119</td><td>83.1</td><td>0.7</td><td>7.8</td></tr><tr><td>16</td><td>106</td><td>80.7</td><td>0.8</td><td>12.8</td></tr><tr><td>17</td><td>111</td><td>64.9</td><td>0.6</td><td>7.6</td></tr><tr><td>18</td><td>116</td><td>74.8</td><td>0.7</td><td>11.5</td></tr><tr><td>19</td><td>107</td><td>78.5</td><td>0.7</td><td>8.7</td></tr><tr><td>20</td><td>112</td><td>81.8</td><td>0.7</td><td>10.1</td></tr><tr><td>21</td><td>87</td><td>128.9</td><td>1.5</td><td>12.8</td></tr><tr><td>22</td><td>103</td><td>71.4</td><td>0.7</td><td>18.3</td></tr><tr><td>23</td><td>91</td><td>72.4</td><td>0.8</td><td>8.8</td></tr><tr><td>24</td><td>90</td><td>71.6</td><td>0.8</td><td>8.8</td></tr><tr><td>25</td><td>96</td><td>71.9</td><td>0.8</td><td>7</td></tr><tr><td>26</td><td>79</td><td>112.5</td><td>1.4</td><td>16.9</td></tr><tr><td>27</td><td>75</td><td>91.5</td><td>1.2</td><td>19.4</td></tr><tr><td>28</td><td>54</td><td>41.8</td><td>0.8</td><td>20.3</td></tr><tr><td>overall</td><td>3588</td><td>2585.7</td><td>0.7</td><td>20.3</td></tr></tbody></table>

growth ratio), indicating antagonistic QTLs (see values of additive effect *a* in Table 3). In contrast, for the two survival traits, the male allele was always unfavorable (lowering the survival value).

The sex of each individual was determined using PCR primers based on the coding sequence of the sex-linked gene Esi0068_0003$^{21,43}$. Using this characteristic as a binary trait, a major sex-associated locus was mapped to LG13. This result confirmed the position of the previously identified sex-determining region$^{21}$.

Twenty-one significant epistatic interactions were detected, involving 38 unique linkage map positions. Among these interactions, only six involved previously detected QTL positions (Supplementary Table S1). For example, in the case of trait T28, the result of epistasis between LG1 and LG14 (Supplementary Fig. S4) showed that the locus on LG1 (at 6cM) had an effect only in the presence of the B genotype (A and B genotypes derived from the male and female parents, respectively) at the locus on LG14 (at 64cM). In the presence of the A genotype at the locus on LG14, no significant difference was observed between the growth ratio of individuals carrying the A genotype and those carrying the B genotype at the locus on LG1.

We also detected five significant QTL × sex interactions, where the effect of the locus was different for the two sexes, representing potential sexually antagonistic alleles (Supplementary Fig. S5).

**Gene ontology enrichment test.** A genome-wide list of candidate genes was constructed by retaining all loci in the 1.5-LOD support confidence intervals of the complete set of QTLs. This list of 562 genes included 263 genes that had associated GO terms. In the re-annotated genome of *Ectocarpus*$^{32}$, 17,418 genes were identified among which 10,688 had predicted functions and 7,383 had associated GO terms. We built a reference list, consisting of genes located in 10 kb windows around the SNPs localized on the genetic map. This gave a list of 2,710 genes among which, 1,156 had associated GO terms (Supplementary Dataset). A GO enrichment test was carried out by comparing the candidate list of 263 genes which had associated GO terms with the reference list of 1,156 genes also having associated GO terms. Six terms were significantly enriched in the candidate gene list compared with the reference at *p* = 0.05, with one term having a p-value < 0.01. Four of those 6 terms were linked to transport processes, the most significant being anion transport (Table 4, Supplementary Fig. S6 and Supplementary Dataset).
<table><caption>**Table 3. QTL mapping results.** T13, T26 and T28 correspond to growth rate phenotypes for temperatures 13°C, 26°C and 28°C respectively. S34, S20 and S15 correspond to growth rate phenotypes for 34‰ salinity, 20‰ and 15‰ salinity respectively. Plasticity traits are represented by a letter followed by a number and a slash, then another letter and number. For example, S34/S20 corresponds to plasticity represented by slope of the reaction norm between the two conditions salinity 34‰ condition and salinity 20‰ condition. Surv_S15 and Surv_S20 correspond to survival trait for conditions salinity 15‰ and salinity 20‰. In the method column, MapQTL_MQM means that the QTL was detected by both MapQTL and R/qlt, RQTL_Scanone means that the QTL was only detected by R/qlt. a: estimated additive effect = (mu_A - mu_B)/2. Mu_A: the estimated mean of the distribution of the quantitative trait associated with the “a” genotype which is the genotype of the male parent. R²: the percentage of the variance explained for by the QTL. 1.5LODLeft and 1.5LODRight represent left and right limits of the 1.5-LOD support interval.</caption><thead><tr><th>Trait</th><th>LG</th><th>Peak (cM)</th><th>LOD</th><th>R<sup>2</sup> (%)</th><th>a</th><th>1.5LOD Left</th><th>1.5LOD Right</th><th>Method</th></tr></thead><tbody><tr><td rowspan="15">Salinity</td><td rowspan="3">S15</td><td>2</td><td>60</td><td>3.83</td><td>23</td><td>0.110</td><td>49</td><td>64</td><td>RQTL_Scanone</td></tr><tr><td>6</td><td>52.686</td><td>4.97</td><td>10</td><td>0.077</td><td>50</td><td>52.7</td><td>MapQTL_MQM</td></tr><tr><td>23</td><td>12.585</td><td>3.96</td><td>9.4</td><td>-0.048</td><td>7.6</td><td>12.6</td><td>MapQTL_MQM</td></tr><tr><td rowspan="6">S20</td><td>1</td><td>54.123</td><td>5.81</td><td>11.2</td><td>-0.075</td><td>53</td><td>55</td><td>MapQTL_MQM</td></tr><tr><td>6</td><td>63.964</td><td>4.62</td><td>9.6</td><td>0.068</td><td>62</td><td>64</td><td>MapQTL_MQM</td></tr><tr><td>21</td><td>74.591</td><td>4.4</td><td>10.7</td><td>0.117</td><td>70.7</td><td>76</td><td>MapQTL_MQM</td></tr><tr><td>24</td><td>71.624</td><td>4.89</td><td>10.2</td><td>-0.099</td><td>68.6</td><td>71.6</td><td>MapQTL_MQM</td></tr><tr><td>26</td><td>36.69</td><td>3.38</td><td>6.8</td><td>-0.090</td><td>31.4</td><td>36.7</td><td>MapQTL_MQM</td></tr><tr><td>26</td><td>56.356</td><td>3.2</td><td>6.4</td><td>0.090</td><td>51.4</td><td>56.3</td><td>MapQTL_MQM</td></tr><tr><td rowspan="4">S20/S15</td><td>7</td><td>74.093</td><td>4.18</td><td>27</td><td>-0.033</td><td>68.4</td><td>86</td><td>MapQTL_MQM</td></tr><tr><td>2</td><td>60.456</td><td>9.43</td><td>32.1</td><td>-0.029</td><td>60</td><td>64</td><td>MapQTL_MQM</td></tr><tr><td>7</td><td>109.287</td><td>3.83</td><td>11</td><td>-0.007</td><td>108</td><td>113</td><td>MapQTL_MQM</td></tr><tr><td>22</td><td>47.593</td><td>4.1</td><td>10.6</td><td>0.016</td><td>44</td><td>47.6</td><td>MapQTL_MQM</td></tr><tr><td rowspan="15">Temperature</td><td rowspan="6">S34/S15</td><td>22</td><td>58.867</td><td>3.53</td><td>9.7</td><td>-0.015</td><td>50</td><td>59</td><td>MapQTL_MQM</td></tr><tr><td>1</td><td>13.812</td><td>3.95</td><td>5.8</td><td>-0.008</td><td>11.2</td><td>13.9</td><td>MapQTL_MQM</td></tr><tr><td>1</td><td>54.123</td><td>15.83</td><td>33.6</td><td>0.020</td><td>54</td><td>56</td><td>MapQTL_MQM</td></tr><tr><td>4</td><td>70.291</td><td>7.9</td><td>13.3</td><td>0.018</td><td>65</td><td>70.3</td><td>MapQTL_MQM</td></tr><tr><td>4</td><td>83.254</td><td>8.62</td><td>15.2</td><td>-0.022</td><td>82.3</td><td>83.3</td><td>MapQTL_MQM</td></tr><tr><td>11</td><td>30.17</td><td>7.61</td><td>12.7</td><td>0.015</td><td>28.8</td><td>30.2</td><td>MapQTL_MQM</td></tr><tr><td rowspan="6">S34/S20</td><td>16</td><td>0</td><td>4.93</td><td>8</td><td>0.011</td><td>0</td><td>0</td><td>MapQTL_MQM</td></tr><tr><td>16</td><td>21.773</td><td>6.43</td><td>10.4</td><td>-0.013</td><td>10.8</td><td>21.8</td><td>MapQTL_MQM</td></tr><tr><td>26</td><td>9.778</td><td>5.63</td><td>8.5</td><td>-0.012</td><td>8.6</td><td>11.5</td><td>MapQTL_MQM</td></tr><tr><td>26</td><td>56.356</td><td>6.61</td><td>10.5</td><td>-0.016</td><td>50.5</td><td>56.4</td><td>MapQTL_MQM</td></tr><tr><td>26</td><td>80.621</td><td>6.39</td><td>9.6</td><td>0.012</td><td>74.4</td><td>80.7</td><td>MapQTL_MQM</td></tr><tr><td>23</td><td>17</td><td>3.75</td><td>5.6</td><td>-0.010</td><td>7</td><td>37.7</td><td>RQTL_Scanone</td></tr><tr><td rowspan="4">Survival_S15</td><td>11</td><td>91.953</td><td>6.48</td><td>13</td><td>-10.408</td><td>91.9</td><td>92.3</td><td>MapQTL_MQM</td></tr><tr><td>19</td><td>50.116</td><td>3.52</td><td>6.3</td><td>-7.211</td><td>45.5</td><td>50.4</td><td>MapQTL_MQM</td></tr><tr><td>22</td><td>37.558</td><td>4.6</td><td>6.7</td><td>-13.560</td><td>36.5</td><td>41</td><td>MapQTL_MQM</td></tr><tr><td>23</td><td>17</td><td>3.75</td><td>5.6</td><td>-0.010</td><td>7</td><td>37.7</td><td>RQTL_Scanone</td></tr><tr><td rowspan="4">Survival_S20</td><td>9</td><td>32.541</td><td>4.18</td><td>19.5</td><td>0.188</td><td>31</td><td>34</td><td>MapQTL_MQM</td></tr><tr><td>2</td><td>57.891</td><td>9.54</td><td>21.3</td><td>0.058</td><td>57.5</td><td>57.9</td><td>MapQTL_MQM</td></tr><tr><td>4</td><td>17.573</td><td>3.65</td><td>6.8</td><td>-0.040</td><td>11.3</td><td>18.6</td><td>MapQTL_MQM</td></tr><tr><td>4</td><td>91.605</td><td>13.35</td><td>33.6</td><td>0.081</td><td>91.6</td><td>94.3</td><td>MapQTL_MQM</td></tr><tr><td rowspan="4">T26</td><td rowspan="4">T28</td><td>14</td><td>0</td><td>5.8</td><td>11.4</td><td>-0.045</td><td>0</td><td>0</td><td>MapQTL_MQM</td></tr><tr><td>18</td><td>15.042</td><td>4.83</td><td>9.3</td><td>0.037</td><td>15</td><td>15.1</td><td>MapQTL_MQM</td></tr><tr><td>24</td><td>37.7</td><td>3.72</td><td>7</td><td>0.014</td><td>33</td><td>63</td><td>RQTL_Scanone</td></tr><tr><td>5</td><td>99.023</td><td>3.65</td><td>13.4</td><td>-0.010</td><td>94</td><td>101</td><td>MapQTL_MQM</td></tr><tr><td>T13/T26</td><td>26</td><td>111.259</td><td>3.26</td><td>12.7</td><td>0.012</td><td>90</td><td>112</td><td>MapQTL_MQM</td></tr><tr><td>T13/T28</td><td>6</td><td>63.964</td><td>4.63</td><td>18.4</td><td>0.069</td><td>63.9</td><td>64.1</td><td>MapQTL_MQM</td></tr><tr><td>T26/T28</td><td>6</td><td>63.964</td><td>4.63</td><td>18.4</td><td>0.069</td><td>63.9</td><td>64.1</td><td>MapQTL_MQM</td></tr><tr><td>Sex</td><td>13</td><td>40.1</td><td>52.1</td><td>75</td><td>0.500</td><td>39</td><td>41</td><td>RQTL_Scanone</td></tr></tbody></table>

Discussion

Growth under different temperature and salinity conditions showed a distribution with the presence of transgressive segregants among the offspring in almost all cases, even when there was very little difference between the parents. Transgressive individuals may have characteristics allowing them to occupy new ecological habitats or to perform more effectively in already occupied environments⁴⁴. These transgressive segregants are likely to be the result of complementary action of different alleles of genes from the male and female parent⁴⁵,⁴⁶ and their
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Genetic map showing the localization of detected QTLs for temperature, salinity and their associated plasticity and survival traits. Blue color indicates salinity and its plasticity traits, red indicates temperature and associate plasticity and green indicates survival traits. The brown color on LG13 represents the sex-determining region and the QTL LOD score curve from a Kruskal-Wallis non-parametric test in MapQTL. T13, T26 and T28 correspond to growth rate phenotypes for temperatures 13 °C, 26 °C and 28 °C respectively. S34, S20 and S15 correspond to growth rate phenotypes for 34‰, 20‰ and 15‰ salinity, respectively. Plasticity traits are represented by a letter followed by a number and a slash, then another letter and number. For example, S34/S20 corresponds to plasticity represented by the slope of the reaction norm between 34‰ and 20‰ salinity. Surv_S15 and Surv_S20 correspond to survival traits for 15‰ and 20‰ salinity, respectively.</figcaption></figure>

<table><caption>**Table 4.** Pearson's correlation coefficients between main studied traits. See Table 2 for trait details.
*Significant at 0.05. ****Significance <0.0001.</caption><thead><tr><th></th><th>T13</th><th>T26</th><th>T28</th><th>S34</th><th>S20</th></tr></thead><tbody><tr><td>T26</td><td>0.095</td><td></td><td></td><td></td><td></td></tr><tr><td>T28</td><td>-0.017</td><td>0.003</td><td></td><td></td><td></td></tr><tr><td>S34</td><td>0.038</td><td>0.217*</td><td>-0.05</td><td></td><td></td></tr><tr><td>S20</td><td>-0.2</td><td>-0.26*</td><td>0.163</td><td>0.025</td><td></td></tr><tr><td>S15</td><td>0.149</td><td>0.061</td><td>0.202</td><td>0.024</td><td>0.531****</td></tr></tbody></table>

presence provides the potential for improvement of the traits using selective breeding or as a result of natural selection.

The two parents of the mapping family originated from the west coast of South America, between the south of Peru and the north of Chile, representing about 2,000 km of coastal ecosystem where seasonal variation in seawater temperature is high (between 13 °C in winter and 25 °C in summer on average, Fig. 3). The parental strains are therefore presumably adapted to living in a changing environment. The annual average temperature of the original location of the male parent is lower than that of the female parent (Fig. 3). And it is noteworthy that it was only at the highest temperature that the female parent performed better than the male one. This may be explained by the mean temperatures at the parental locations, but since 28 °C is not a common mean temperature encountered there, this is rather speculative. Nevertheless, their occurrence in rockpools or low tides during summer months can expose these *Ectocarpus* strains to very high temperatures. Furthermore, the average seasonal temperature of seawater in the distribution zone of the species can increase by 4 °C to 10 °C during *El Niño*-Southern Oscillation events⁴⁷.

The distributions of growth ratios varied little at 13 °C and 26 °C but significantly decreased at 28 °C. Even at 28 °C, the majority of the offspring could maintain non-null growth. A temperature of 30 °C may exceed the tolerance of this family as observed during preliminary tests, at least under our laboratory conditions, but it is possible that the response of the population to extreme temperatures would show some differences under field conditions compared to the observations in laboratory conditions, particularly because several stressors act in synergy in the wild.

For salinity stress, the female parent always performed better than the male parent. Typical salinity of seawater in the distribution range of the species is 34–36‰. However, our data showed that several strains were able to maintain active growth even at 15‰. Although most of the species of the *Ectocarpus* genera are predominantly marine, some strains (e.g. *E. subulatus* isolates) are tolerant to low salinity (i.e. occur in estuaries) or are even adapted to a freshwater environment³¹,⁴⁸. The maintenance of a capacity to respond to large variations in water salinity in the extremely arid region of the Atacama Desert, where there are no estuaries, is likely a result of an adaptive evolution in other regions of the coast. Interestingly, this low salinity tolerance in *Ectocarpus* species may explain why there are cases of colonization of freshwater environments, albeit rare. Dittami *et al*.,³¹ suggested that, while high stress tolerance and plasticity may be prerequisites for the colonization of freshwater, genomic alterations are likely to have occurred and these alterations would have produced permanent changes in metabolite profiles to stabilize the transition.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Sea surface temperature (sst) in the distribution range of *Ectocarpus* sp. for year 2014. The blue arrow represents the original location of the male parent (San Juan de Marcona, Peru) and the red arrow represents the original location of the female parent (Arica, Chile). Created with the SeaWiFS Data Analysis System (SeaDAS) visualization and data processing software version 7.2, using seasonal (winter and summer) and annual MODIS sst data from 2014: http://seadas.gsfc.nasa.gov/installers/; http://oceandata.sci.gsfc.nasa.gov/MODIS-Aqua/Mapped. NASA Goddard Space Flight Center, Ocean Biology Processing Group; (2014): Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Ocean Color Data, NASA OB.DAAC, Greenbelt, MD, USA. http://doi.org/10.5067/ORBVIEW-2/SEAWIFS_OC.2014.0. Accessed 2015/05/11. Maintained by NASA Ocean Biology Distributed Active Archive Center (OB.DAAC), Goddard Space Flight Center, Greenbelt, MD, USA.</figcaption></figure>

Broad-sense heritability estimates were low to moderate but increased with the degree of stress for the environmental conditions we tested (from 0.20 to 0.78). The heritability estimates were lower when temperature and salinity stress were combined (0.10 and 0.31). It is well known that estimated heritability of a trait can vary either with environmental conditions or between populations. Hoffmann and Merilä⁴⁹ summarized hypotheses to explain the relationship between environmental effects and heritability. It is possible that phenotypic variations between genotypes would be difficult to detect unless resources become limiting, regardless of the history of selection⁵⁰. In contrast, unfavorable conditions could also decrease the importance of observed genetic variation if increased environmental variance reduces the adaptive value of individual genotypic combinations. This might apply to our study, when we combined stress conditions for temperature and for salinity, hence increasing the overall environmental heterogeneity. Furthermore, although we were able to replicate genotypes in this studied species, traits are nonetheless measured with some error and this might lead to an underestimate of heritability if phenotypic variance is higher under unfavorable conditions⁵¹.

With 3,588 SNP markers, the genetic map described here is considerably denser than the previously published map for *Ectocarpus* sp., which was based on 406 microsatellite markers²⁰. Overall, the distribution of the markers on the 28 LG was homogenous. The 28 LGs detected here are coherent with cytogenetic studies carried out using European strains of *E. siliculosus* that indicated an approximate number of 25 chromosomes⁵².

Assignation of the SNP markers to sequence scaffolds allowed 500 of the 1,561 scaffolds of the genome sequence to be anchored onto linkage groups, significantly improving the large-scale assembly of the genome³². Integration of this new information along with the data from the microsatellite-based linkage map²⁰ has allowed a total of 531 scaffolds to be anchored onto linkage groups, corresponding to 177 Mbp of sequence data and 90.5% of the total genome. The improved genetic map will represent an important resource for future genetic and comparative genomic studies using *Ectocarpus* sp. as a model species.

The quality of a genetic map is influenced by the density of the markers used and the size of the mapping family. As far as marker density is concerned, next generation sequencing technologies represent a means to tremendously improve the generation of genetic markers and the ddRAD sequencing approach used in this study produced enough markers to cover the major part of the genome. However, the mapping family is a "pseudo F2" haploid population where each individual is derived from a separate but unique meiosis. This, plus the relatively small size of the family, constitute factors that may limit the power of QTL detection⁵³.

We estimated a global recombination rate of 12.28 cM/Mb for *Ectocarpus* sp., an estimate that is higher than values obtained for most terrestrial plant species⁵⁴. If we examine data available for *Saccharina japonica*⁵⁵,⁵⁶, the only other brown algal species for which both a genetic map and a complete genome sequence currently exist, the global recombination rate is 4.93 cM/Mb. This value is more similar to those observed for angiosperms. The genome size of *Saccharina japonica* is 537 cM⁵⁶, which is 2.5 times larger than that of *Ectocarpus* sp. These estimates are in accordance with the observation of negative correlation between global recombination rate and genome size, pointing to a putative role of LTR retrotransposons⁵⁴. More estimates are needed for additional brown algae to provide improved insights into global recombination rates in this group. In any case, recombination is a very complex feature because it is known to vary substantially across eukaryote genomes and even
<table><caption>**Table 5.** List of the top 20 GO terms identified by the test for GO term enrichment in genes located within the QTL intervals identified in this study. See materials and methods for details.</caption><thead><tr><td>GO ID</td><td>Term</td><td>Annotated</td><td>Significant</td><td>Expected</td><td>Rank in classic</td><td>p-value classic</td></tr></thead><tbody><tr><td>GO:0015698</td><td>inorganic anion transport</td><td>3</td><td>2</td><td>0.15</td><td>1</td><td>0.007</td></tr><tr><td>GO:0006820</td><td>anion transport</td><td>4</td><td>2</td><td>0.2</td><td>2</td><td>0.014</td></tr><tr><td>GO:0015031</td><td>protein transport</td><td>22</td><td>4</td><td>1.1</td><td>3</td><td>0.02</td></tr><tr><td>GO:0008104</td><td>protein localization</td><td>24</td><td>4</td><td>1.2</td><td>5</td><td>0.027</td></tr><tr><td>GO:0006396</td><td>RNA processing</td><td>6</td><td>2</td><td>0.3</td><td>8</td><td>0.032</td></tr><tr><td>GO:0006886</td><td>intracellular protein transport</td><td>17</td><td>3</td><td>0.85</td><td>9</td><td>0.048</td></tr><tr><td>GO:0006605</td><td>protein targeting</td><td>8</td><td>2</td><td>0.4</td><td>12</td><td>0.056</td></tr><tr><td>GO:0072594</td><td>establishment of protein localization</td><td>8</td><td>2</td><td>0.4</td><td>13</td><td>0.056</td></tr><tr><td>GO:1902580</td><td>single-organism cellular localization</td><td>9</td><td>2</td><td>0.45</td><td>15</td><td>0.07</td></tr><tr><td>GO:0006811</td><td>ion transport</td><td>21</td><td>3</td><td>1.05</td><td>16</td><td>0.082</td></tr><tr><td>GO:0016482</td><td>cytoplasmic transport</td><td>11</td><td>2</td><td>0.55</td><td>17</td><td>0.1</td></tr><tr><td>GO:0006810</td><td>transport</td><td>103</td><td>8</td><td>5.14</td><td>22</td><td>0.128</td></tr><tr><td>GO:0051234</td><td>establishment of localization</td><td>104</td><td>8</td><td>5.19</td><td>23</td><td>0.133</td></tr><tr><td>GO:0009148</td><td>pyrimidine nucleoside triphosphate biosynthesis</td><td>3</td><td>1</td><td>0.15</td><td>25</td><td>0.142</td></tr><tr><td>GO:0043094</td><td>cellular metabolic compound salvage</td><td>3</td><td>1</td><td>0.15</td><td>26</td><td>0.142</td></tr><tr><td>GO:0007005</td><td>mitochondrion organization</td><td>3</td><td>1</td><td>0.15</td><td>24</td><td>0.142</td></tr><tr><td>GO:0044743</td><td>intracellular protein transmembrane importation</td><td>3</td><td>1</td><td>0.15</td><td>27</td><td>0.142</td></tr><tr><td>GO:0009147</td><td>pyrimidine nucleoside triphosphate metabolism</td><td>3</td><td>1</td><td>0.15</td><td>28</td><td>0.142</td></tr><tr><td>GO:0017038</td><td>protein import</td><td>3</td><td>1</td><td>0.15</td><td>29</td><td>0.142</td></tr><tr><td>GO:0051179</td><td>localization</td><td>106</td><td>8</td><td>5.29</td><td>32</td><td>0.144</td></tr></tbody></table>

between individuals within populations, for reasons that are still to be fully understood⁵⁷. In our study, recombination rates varied along and between the different *Ectocarpus* sp. LGs (average rate per LG between 11.11 cM/Mb for LG28 and 17.93 cM/Mb for LG4), at a 1 Mb scale. Although no clear differential patterns between the LGs were observed, this result should be considered with caution for several reasons. First, resolution of cross-based estimates of recombination rates like those we provide here depends on the number of recombination events occurring in a generation, and the generation of the family analyzed here involved only one meiosis per individual. Second, because the scale at which the recombination rates are examined may influence the detected relationships between recombination and genomic features and, consequently, patterns at coarse scales may not accurately reflect relationships at finer scales⁵⁸.

Our experimental setup allowed us to detect 38 QTLs for growth-related traits and their plasticity under different temperature and salinity conditions and 1 QTL for sex. Eighteen of the 28 LGs contained QTLs for these traits. Although several of the QTLs had large confidence intervals, some showed relatively small ones. For example, the QTLs for T28 on LG2 and LG18 had small confidence intervals, which contained respectively 41 and 6 annotated or predicted genes. There was only one case of co-localization of QTLs (for temperature and salinity on LG2 for traits T28 and S15). The rarity of co-localization of QTLs was in agreement with the low to moderate correlation between traits, with most correlations being statistically non-significant (Table 5). More co-localizations were observed when survival and plasticity traits were considered. In three cases (LG1, LG2 and LG26), plasticity QTLs co-localized with main traits. These co-localizations could be due to the presence of genes that control metabolic pathways that influence multiple traits or to the presence of pleiotropic genes.

Phenotypic plasticity is one of the strategies used by organisms to cope with variability in their environment but its genetic basis is often poorly documented⁵⁹. Plasticity often involves altering gene expression and physiology in response to environmental cues⁶⁰. The QTL approach developed here could allow evolutionary questions such as the evolution of the reaction norm as a response to the changing environment to be addressed. The detection of plasticity QTLs in *Ectocarpus* sp. indicates that factors controlling plasticity in this species are likely to be genetically-based⁴¹. We also identified several antagonistic QTLs (i.e. QTLs with alleles that have opposing effects). It is thought that, in such cases, trait differences may have arisen through drift and/or through selection that fluctuates in direction⁴⁶,⁶¹. Antagonistic QTLs are also considered to be able to generate transgressive segregants⁴⁶, which is in accordance with our results.

Epistasis is considered to be an important genetic basis of complex phenotypes⁶²,⁶³. More than one trait was involved in all the 21 cases of epistatic interactions detected in this study, suggesting pleiotropy not only for single loci but also for epistatic effects. Among those 21 epistatic interactions, we observed only six cases where one of the markers had a main QTL effect on the trait. This suggests that epistasis is an important mechanism for the control of the response to stress in *Ectocarpus* sp., at least when growth ratio is used as a proxy for fitness.

Few QTL × sex interactions were identified in our study. These interactions were sex-antagonistic since the effects of the QTL alleles were opposite for the two sexes. Such interactions indicate that the effect of a QTL is different between the two sexes and this can be important in evolutionary as well as in breeding contexts. Evolutionary arguments implicate QTL × sex interactions in the maintenance of variation for complex traits⁶⁴. From a biological point of view, QTL × sex interactions can be explained by invoking interactions between genes with different effects in different sexual environments (considering sexes as different environments and hence
QTL × sex interactions to be QTL × environment interactions)⁶⁵. However, we cannot exclude the possibility that the observed QTL × sex interactions observed in our study may just reflect the different origins of the parental strains and their specific trait inheritance, since we cannot separate the effects of their sex and their origin. The detection of main effects, epistasis and sexually-dimorphic QTLs suggest complex genetic regulation of growth under stress conditions.

A GO enrichment test identified six terms significantly enriched in putative genes lying in the QTL intervals, among which, four were linked to ion and protein transport processes. Of the 562 putative genes in the QTL intervals, 102 were also identified as being significantly differentially expressed between low and high salinity by Dittami *et al.*³¹ (Supplementary Dataset), but more interestingly, the authors also observed an accumulation of intracellular ions at different levels between freshwater and marine strains, suggesting an important role for ion transport in the adaptation to low salinities. The list of putative stress-related genes corresponding to the QTL intervals identified in the current study should be treated with caution because there are uncertainties linked to positions of QTLs and because growth is a complex phenotype integrating numerous metabolic pathways. Taking into account these limitations, the data presented here is however in accordance both with the suggested importance of intracellular ion composition and osmotic balance in adaptation to stresses such as low salinity and with the large number of ion channels observed in the genome of the species¹,³¹. Furthermore, transport processes such as ion transport are well studied in land plants and have been shown to be involved in the response to osmotic stresses¹⁶,⁶⁶. In seagrasses, several genes have been also shown to respond to elevated sea surface temperatures⁶⁷ or to be under selection along the salinity gradient of the intertidal shore⁶⁸.

In our study, we built a high-density genetic map using ddRAD-seq approaches and detected several QTLs for growth rate under high temperature and low salinity stresses as well as related plasticity and survival traits. GO enrichment tests suggested a role for ion transport genes among the main loci controlling the QTLs. This study represents a significant step towards deciphering the genetic architecture of adaptation of *Ectocarpus* sp. to stress conditions and, furthermore, adds a substantial resource to the increasing list of resources generated for the species. It also opens new perspectives in population genomics of adaptation in brown algal species. As mentioned in the introduction, brown algae have also been shown to display important interactions with complex bacterial communities, which influence their ability to respond to environmental stressors and this aspect makes the picture more complex. Several studies have shown that algal-bacterial associations are beneficial for the algal hosts, increasing their fitness⁶⁹,⁷⁰ and some brown algal species such as *Ectocarpus fasciculatus* and *Pylaiella littoralis* are absolutely dependent on their associated bacteria⁷¹,⁷². These complex interactions were not included in this study but should be examined in future studies in order to have a more complete understanding of the adaptation of *Ectocarpus* sp. to its stressful environments.

## Materials and Methods

**Mapping family.** The mapping family was generated by crossing the male *Ectocarpus* sp. strain Ec 32 (whose genome has been fully sequenced) with a compatible female strain, Ec 568. See Heesch *et al.*²⁰ for a detailed description of the production of the family. Briefly, a single F1 hybrid sporophyte isolated from the cross (Ec 569) was used to produce a family of sibling haploid gametophytes, each individual being derived from a separate meiotic event. Gametes from each individual were then allowed to germinate parthenogenetically to produce haploid partheno-sporophytes. Eighty-nine progeny of this haploid family, plus haploid partheno-sporophytes derived from the parental gametophytes, were cultivated for the experiments.

**Culture conditions and phenotyping.** For both the temperature and salinity stress experiments, preliminary tests were carried out to determine conditions that resulted in phenotypic differences without causing extreme mortality. Three temperature conditions were selected: 13 °C (control condition) and two higher stress temperatures, 26 °C and 28 °C. Similarly, three salinity conditions were also selected: 34‰ (normal sea salinity) as the control and two low salinity conditions, 20‰ and 15‰. For each strain, 6 replicates were used and individuals were randomly distributed in 6-well plates containing 6 ml of autoclaved, Provasoli-enriched seawater. Natural seawater was used for temperature experiments (salinity is stable at around 33.5‰ to 34.5‰) while artificial seawater was used for salinity experiments. Note that differences might exist between natural and artificial seawater regarding nutrients. Experiments were carried out in growth chambers and the positions of the plates were permutated daily to reduce microenvironment effects. Light conditions were 12 h/12 h day/night (~2000 lux light intensity) and the temperature in the growth chamber was set to 13 °C for the salinity experiment. The plates were maintained in the growth chambers for 12 days. Each individual was photographed under a binocular microscope before and after the experiment and thallus area was measured using the program ImageJ. Growth was estimated as the ratio of the thallus area on day 12 to the thallus area on day 0. By treating different experimental conditions as different environments, we determined reaction norms of the mapping progeny and used the slopes of these reaction norms as plasticity traits. The relation of these plasticity traits to specific environmental variations was indicated by giving the two environmental conditions separated by a slash mark (e.g.: S34/S20 is the plasticity trait between salinity conditions 34‰ and 20‰).

**DNA extraction and generation of RAD-seq data.** For DNA extraction, the haploid family together with the parental lines were grown in a growth chamber in sterile seawater medium for 6 weeks under a 12-hour light photoperiod and 13 °C. As the cultures were not fully axenic, 5 ml of a mix of antibiotics (9 mg/ml of Penicillin G, 4.5 mg/ml of Streptomycin and 0.9 mg/ml of chloramphenicol) were added to each liter of seawater used for the culture to reduce the bacterial load. Harvested individuals were frozen and lyophilized, and DNA was extracted using the NucleoSpin® 96 Plant II kit (Macherey-Nagel GmbH & Co. KG, Germany), according to the manufacturer's instructions. DNA quality was checked on agarose gels and the quantity measured by PicoGreen® (Fisher Scientific).
Preparation of the double digest RAD library (ddRAD-seq) was carried out following Peterson et al.⁴⁰. The genomic DNA was digested with the restriction enzymes PstI and HhaI (New England Biolabs, https://www.neb.com/), which had been selected based on in silico digestion simulations, and adapters were ligated using the T4 DNA ligase (New England Biolabs, https://www.neb.com/). After pooling the samples which had been individually barcoded with a unique adapter and cleaning with AMPure XP beads (Beckman Coulter Genomics), size selection (between 150 and 600 bp) was carried out with a Pippin-Prep kit (Sage Science, Beverly, MA, USA). A PCR amplification was then carried out using the Q5® hot Start High-Fidelity DNA polymerase kit (New England Biolabs) to increase the concentration of the libraries and to add Illumina flowcell annealing sequences, multiplexing indices and sequencing primer annealing regions, to all fragments. Compared with the original protocol, the barcodes were modified by adding an extra nucleotide, extending to six nucleotides (see Supplementary Dataset). Both an Agilent® 2100 Bioanalyzer® (Agilent Technologies) and qPCR were used to quantify the libraries generated from each individual. Libraries with distinct multiplexing indices were then combined in equimolar ratios to compose a final pool of libraries for sequencing. An Illumina HiSeq 2500 platform (Rapid Run Mode) was used to generate a total of 470 million (48.5 Gb) high quality (Q30 = 82%) 100 bp paired-end sequence reads at a cluster density of 850 K/mm² (HiSeq Rapid SBS Kit v2, 200 cycles). As a test, a subset of the library was also sequenced with Illumina MiSeq, providing 93 million additional reads (MiSeq Reagent Kit v3, 150 cycles). The raw sequence data in fastq format are stored in the Sequence Read Archive (SRA) at National Center for Biotechnology Information (NCBI) and can be accessed at NCBI homepage (https://www.ncbi.nlm.nih.gov/; bioproject accession PRJNA371840 or study accession SRP099384).

**Analysis of ddRAD-seq data.** The ddRAD-seq sequencing data was analyzed using the Stacks pipeline (version 1.40)⁷³. Raw sequence reads were checked for intact barcodes and the restriction enzyme sites. A sliding window of 25% of the length of a read was set to check sequence quality and reads with <90% base call accuracy were discarded. Because Stacks removes whole reads when adapter contamination is detected, we first used the program Cutadapt (version 1.8.3) to check for the presence of adapters and cut only adapter sequences. Paired-end sequencing of shorter fragments generated overlapping reads that were treated with the program PEAR (version 0.9.10)⁷⁴. PEAR identifies paired-end non-overlapping reads and generates a single consensus sequence from overlapping read pairs. These single consensus sequences were added to the singleton rem1 and rem2 sequences produced by Stacks for a unique group of singleton sequences. Paired-end reads as well as singletons were then trimmed to 95 bp with the program TRIMMOMATIC (version 0.32)⁷⁵. The paired-end sequences and the singletons were then aligned to the *Ectocarpus* reference genome¹ using Bowtie 2 (version 0.12.9)⁷⁶. Output aligned .sam files were subsequently imported into the Stacks pipeline and the final Genotypes program exported haplotypes encoded as genotypes at the different loci.

**Construction of the genetic map.** After Stacks analysis, genotypes were exported in a generic format that was imported into Excel for filtering. As we used a haploid population, only aa/bb segregation was expected. All non-expected segregation patterns, as well as loci with log likelihood values below -30, were removed. Loci with over 25% of missing data were also removed. In the finale matrix (see Supplementary Dataset), A genotypes derived from the male parent and B genotypes derived from the female parent. The R package R/qlt (version 1.39-5)⁷⁷ was used for the construction of the genetic map, in conjunction with either R/ASMap (version 0.4-7) or JoinMap® version 4.1. At this stage, individuals with more than 40% of missing genotypes were removed. Further filtering was done within the mapping programs. Markers that showed significant segregation distortion (≤5%) were removed. Recombination fraction estimates were used to correct erroneous genotypes or switched alleles. Markers were also filtered based on their similarity, on whether they were located on the same contig and based on suspect linkage patterns. In JoinMap, the maximum likelihood mapping algorithm was applied.

Estimated genome (Ge) length was obtained by adding 2s (s being the average spacing between markers along the LGs) to the length of each LG to account for chromosome ends beyond the terminal markers of each LG⁷⁸. The map coverage was then estimated as the ratio between the observed map length and Ge.

The R/xoi package (version 0.66-9) was used to obtain a smoothed estimate of the recombination rate along the LGs, in 1 Mbp sliding windows.

**Statistical analyses of phenotypic data.** Statistical treatments were done in R. BLUPs were obtained by fitting several linear mixed models (LMM) with the R package lme4 (version 1.1-12) and selecting the optimal model with AIC:

$$Y_{ij} = \mu + G_i + R_j + \varepsilon_{ij} \qquad (1)$$

where $Y_{ij}$ represents the value of the trait under investigation for the genotype $i$ at replicate $j$, $\mu$ the general mean of the trait, $G_i$ the random effect of the genotype $i$, $R_j$ the random effect of the replicate $j$ and $\varepsilon_{ij}$ the random residual effect.

Broad-sense heritability of individual traits was calculated as:

$$H^2 = \frac{\sigma_G^2}{\sigma_G^2 + \frac{\sigma_j^2}{J}} \qquad (2)$$

where $\sigma_G^2$ is the genetic variance, $\sigma_j^2$ the environmental variance and $J$ the number of replicates.
Some strains did not survive at salinities of 15‰ and 20‰ and hence we introduced a binomial survival trait. The survival traits were analyzed by fitting generalized linear mixed-effects models (GLMM) and using $\frac{\pi^2}{3}$ as an approximate to the environmental variance as suggested by Nakagawa and Schielzeth$^{79}$.

The different temperature, salinity or survival conditions were also combined as a single phenotype (temperature, salinity or survival). Variances and broad-sense heritabilities were obtained by fitting LMMs or GLMMs and selecting the optimal model:

$$Y_{ijk} = \mu + G_i + R_j + E_k + GE_{ik} + \varepsilon_{ijk} \quad (3)$$

where $E_k$ is the fixed effect of environment (different conditions of temperature, salinity or survival) and $G_{ik}$ the random interaction effect between genotype i and environment k. The sex of each strain was included in the data. Broad-sense heritability of individual traits was calculated as:

$$H^2 = \frac{\sigma_G^2}{\sigma_G^2 + \frac{\sigma_{GE}^2}{K} + \frac{\sigma_{IK}^2}{IK}} \quad (4)$$

where $\sigma_{GE}^2$ is the $G \times E$ interaction variance, K the number of environments and J the number of replicates.

**QTL mapping.** QTL mapping was performed using the R package R/qlt$^{77}$ as well as MapQTL$^{®}$ version 5. In R/qlt, the scanone function was used to perform interval mapping with the Haley-Knott regression method, including sex as a cofactor. 5% significance thresholds of the LOD scores were obtained by 5,000 permutations. In MapQTL, composite interval mapping was performed with the MQM method. The method uses markers as cofactors in an approximate multiple-QTL model with additive and dominant gene actions. Cofactors were selected in a first round as loci with LOD $\ge$2.5 after a simple interval mapping. Then this list was extended with markers on a given linkage group selected every 10 or 20 cM and optimized via the automatic cofactor selection method. The method uses a backward elimination procedure to determine which markers show a significant association and which do not, and ultimately produces a list of significant cofactor markers. For survival traits and sex, the “binary” or the “nonparametric” models were applied. BLUP values were used for QTL mapping. For QTL confidence intervals, 1.5-LOD supports were used. Mapchart 2.3$^{80}$ was used to draw the QTL on the linkage map.

QTL × sex interactions were detected by using “sex” as a covariate. QTL × QTL additive epistatic interactions were detected by the scantwo function of R/qlt, performed by the Haley-Knott regression. Four models were used for this: a full model where the two QTLs are allowed to interact, an additive model where the QTLs are assumed to act additively, a single-QTL model, which corresponds to the result of the single-QTL genome scan and a null model assuming no QTL$^{81}$.

**Gene ontology enrichment test.** Taking advantage of the fact that the genome of the species has been sequenced and carefully annotated, we tested GO enrichment among the putative genes within the QTL intervals. For that, we made two lists: a candidate list consisting of all loci in the 1.5-LOD support confidence intervals of the QTLs, based on the sequenced reference genome and a reference list including all the putative genes in the mapped scaffolds, refining the search only on a 10 Kb window around each SNP on the genetic map. Only loci with clear GO identifications were retained from the two lists. The candidate list was then compared to the reference list using the R/topGO package (version 2.26.0), which is designed to facilitate semi-automated enrichment analysis for GO terms. The test was based on gene count and used the “classic” algorithm with the “Fisher” statistic.

## References

1. Cock, J. M. et al. The *Ectocarpus* genome and the independent evolution of multicellularity in brown algae. *Nature* **465**, 617–621, doi: 10.1038/nature09016 (2010).

2. Cock, J. M. & Collén, J. In *Evolutionary Transitions to Multicellular Life, Advances in Marine Genomics* 2 (eds I. Ruiz-Trillo & A. M. Nedelcu) Ch. *Independent Emergence of Complex Multicellularity in the Brown and Red Algae* 335–361 (Springer, 2015).

3. Dayton, P. K. Ecology of kelp communities. *Annual review of ecology and systematics*, 215–245 (1985).

4. Steneck, R. S. et al. Kelp forest ecosystems: biodiversity, stability, resilience and future. *Environmental conservation* **29**, 436–459 (2002).

5. Wei, N., Quarterman, J. & Jin, Y.-S. Marine macroalgae: an untapped resource for producing fuels and chemicals. *Trends in biotechnology* **31**, 70–77 (2013).

6. Smit, A. J. Medicinal and pharmaceutical uses of seaweed natural products: a review. *Journal of applied phycology* **16**, 245–262 (2004).

7. Michel, G., Tonon, T., Scornet, D., Cock, J. M. & Kloareg, B. Central and storage carbon metabolism of the brown alga *Ectocarpus siliculosus*: insights into the origin and evolution of storage carbohydrates in Eukaryotes. *New Phytologist* **188**, 67–81, doi: 10.1111/j.1469-8137.2010.03345.x (2010).

8. Davison, I. R. & Pearson, G. A. Stress tolerance in intertidal seaweeds. *Journal of Phycology* **32**, 197–211 (1996).

9. Bustamante, R. H. et al. Gradients of intertidal primary productivity around the coast of South Africa and their relationships with consumer biomass. *Oecologia* **102**, 189–201 (1995).

10. Menge, B. A. Top-down and bottom-up community regulation in marine rocky intertidal habitats. *Journal of Experimental Marine Biology and Ecology* **250**, 257–289 (2000).

11. Mota, C. F., Engelen, A. H., Serrao, E. A. & Pearson, G. A. Some don't like it hot: microhabitat-dependent thermal and water stresses in a trailing edge population. *Functional Ecology* **29**, 640–649 (2015).

12. Sanford, E. & Kelly, M. W. Local adaptation in marine invertebrates. *Annual Review of Marine Science* **3**, 509–535 (2011).

13. Kawecki, T. J. & Ebert, D. Conceptual issues in local adaptation. *Ecology Letters* **7**, 1225–1241, doi: 10.1111/j.1461-0248.2004.00684.x (2004).
14. Padilla, D. K. & Savedo, M. M. A systematic review of phenotypic plasticity in marine invertebrate and plant systems. *Advances in marine biology* **65**, 67–94 (2013).

15. Apel, K. & Hirt, H. Reactive oxygen species: Metabolism, oxidative stress, and signal transduction. *Annual Review of Plant Biology* **55**, 373–399, doi: 10.1146/annurev.arplant.55.031903.141701 (2004).

16. Wang, W. X., Vinocur, B. & Altman, A. Plant responses to drought, salinity and extreme temperatures: towards genetic engineering for stress tolerance. *Planta* **218**, 1–14, doi: 10.1007/s00425-003-1105-5 (2003).

17. Wahid, A., Gelani, S., Ashraf, M. & Foolad, M. R. Heat tolerance in plants: An overview. *Environmental and Experimental Botany* **61**, 199–223, doi: 10.1016/j.envexbot.2007.05.011 (2007).

18. Dittami, S. M. et al. Integrative analysis of metabolite and transcript abundance during the short-term response to saline and oxidative stress in the brown alga *Ectocarpus siliculosus*. *Plant, cell & environment* **34**, 629–642 (2011).

19. Dittami, S. M. et al. Global expression analysis of the brown alga *Ectocarpus siliculosus* (Phaeophyceae) reveals large-scale reprogramming of the transcriptome in response to abiotic stress. *Genome Biology* **10**, doi: 10.1186/gb-2009-10-6-r66 (2009).

20. Heesch, S. et al. A sequence-tagged genetic map for the brown alga *Ectocarpus siliculosus* provides large-scale assembly of the genome sequence. *New Phytologist* **188**, doi: 10.1111/j.1469-8137.2010.03273.x (2010).

21. Ahmed, S. et al. A haploid system of sex determination in the brown alga *Ectocarpus* sp. *Current Biology* **24**, 1945–1957 (2014).

22. Montecinos, A. E. et al. Species delimitation and phylogeographic analyses in the *Ectocarpus* subgroup *siliculosi* (Ectocarpales, Phaeophyceae). *Journal of Phycology* (2016).

23. Peters, A. F., Couceiro, L., Tsiamis, K., Küpper, F. C. & Valero, M. Barcoding of cryptic stages of marine brown algae isolated from incubated substratum reveals high diversity in Acinetosporaceae (Ectocarpales, Phaeophyceae). *Cryptogamie, Algologie* **36**, 3–29 (2015).

24. Peters, A. F., Scornet, D., Muller, D. G., Kloareg, B. & Cock, J. M. Inheritance of organelles in artificial hybrids of the isogamous multicellular chromist alga *Ectocarpus siliculosus* (Phaeophyceae). *European Journal of Phycology* **39**, 235–242, doi: 10.1080/09670260410001683241 (2004).

25. Couceiro, L. et al. Evolution and maintenance of haploid-diploid life cycles in natural populations: The case of the marine brown alga *Ectocarpus*. *Evolution* **69**, 1808–1822, doi: 10.1111/evo.12702 (2015).

26. Bold, H. C. & Wynnne, M. *Introduction to the algae: structure and reproduction* (Prentice Hall, 1978).

27. McCauley, L. A. & Wehr, J. D. Taxonomic reappraisal of the freshwater brown algae *Bodanella*, *Ectocarpus*, *Heribaudiella*, and *Pleurocladia* (Phaeophyceae) on the basis of rbc L sequences and morphological characters. *Phycologia* **46**, 429–439 (2007).

28. West, J. & Kraft, G. *Ectocarpus siliculosus* (Dillwyn) Lyngb. from Hopkins River Falls, Victoria-the first record of a freshwater brown alga in Australia. *Muelleria* **9**, 29–33 (1996).

29. Dittami, S. M. et al. Host-microbe interactions as a driver of acclimation to salinity gradients in brown algal cultures. *The ISME journal* **10**, 51–63 (2016).

30. Lee, C. E. & Bell, M. A. Causes and consequences of recent freshwater invasions by saltwater animals. *Trends in Ecology & Evolution* **14**, 284–288 (1999).

31. Dittami, S. M. et al. Towards deciphering dynamic changes and evolutionary mechanisms involved in the adaptation to low salinities in *Ectocarpus* (brown algae). *The plant journal* **71**, 366–377 (2012).

32. Cormier, A. et al. Re-annotation, improved large-scale assembly and establishment of a catalogue of noncoding loci for the genome of the model brown alga *Ectocarpus*. *New Phytologist* (2016).

33. Mackay, T. F., Stone, E. A. & Ayroles, J. F. The genetics of quantitative traits: challenges and prospects. *Nature Reviews Genetics* **10**, 565–577 (2009).

34. Price, A. H. Believe it or not, QTLs are accurate! *Trends in plant science* **11**, 213–216 (2006).

35. Collard, B. C. Y., Jahufer, M. Z. Z., Brouwer, J. B. & Pang, E. C. K. An introduction to markers, quantitative trait loci (QTL) mapping and marker-assisted selection for crop improvement: The basic concepts. *Euphytica* **142**, 169–196, doi: 10.1007/s10681-005-1681-5 (2005).

36. Ren, Z. H. et al. A rice quantitative trait locus for salt tolerance encodes a sodium transporter. *Nature Genetics* **37**, 1141–1146, doi: 10.1038/ng1643 (2005).

37. Sauvage, C. et al. QTL for resistance to summer mortality and OsHV-1 load in the Pacific oyster (*Crassostrea gigas*). *Animal genetics* **41**, 390–399 (2010).

38. Shan, T. F., Pang, S. J., Li, J., Li, X. & Su, L. Construction of a high-density genetic map and mapping of a sex-linked locus for the brown alga Undaria pinnatifida (Phaeophyceae) based on large scale marker development by specific length amplified fragment (SLAF) sequencing. *Bmc Genomics* **16**, doi: 10.1186/s12864-015-2184-y (2015).

39. Baird, N. A. et al. Rapid SNP discovery and genetic mapping using sequenced RAD markers. *PloS one* **3**, e3376 (2008).

40. Peterson, B. K., Weber, J. N., Kay, E. H., Fisher, H. S. & Hoekstra, H. E. Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. *PloS one* **7**, e37135 (2012).

41. Lacaze, X., Hayes, P. & Korol, A. Genetics of phenotypic plasticity: QTL analysis in barley, *Hordeum vulgare*. *Heredity* **102**, 163–173 (2009).

42. Assmann, S. M. Natural variation in abiotic stress and climate change responses in *Arabidopsis*: implications for twenty-first-century agriculture. *International Journal of Plant Sciences* **174**, 3–26 (2013).

43. Lipinska, A. P. et al. Development of PCR-based markers to determine the sex of kelp s. *Plos One* **10**, doi: 10.1371/journal. pone.0140535 (2015).

44. Lewontin, R. & Birch, L. Hybridization as a source of variation for adaptation to new environments. *Evolution*, 315–336 (1966).

45. Tanksley, S. QTL analysis of transgressive segregation in an interspecific tomato cross. *Genetics* **134**, 585–596 (1993).

46. Rieseberg, L. H., Widmer, A., Arntz, A. M. & Burke, B. The genetic architecture necessary for transgressive segregation is common in both natural and domesticated populations. *Philosophical Transactions of the Royal Society of London B: Biological Sciences* **358**, 1141–1147 (2003).

47. Peters, A. & Breeman, A. Temperature tolerance and latitudinal range of brown algae from temperate Pacific South America. *Marine Biology* **115**, 143–150 (1993).

48. Russell, G. & Bolton, J. J. Euryhaline ecotypes of *Ectocarpus siliculosus* (Dillw.) Lyngb. *Estuarine and Coastal Marine Science* **3**, 91–94 (1975).

49. Hoffmann, A. A. & Merilä, J. Heritable variation and evolution under favourable and unfavourable conditions. *Trends in Ecology & Evolution* **14**, 96–101 (1999).

50. Ward, P. J. Parent-offspring regression and extreme environments. *Heredity* **72**, 574–581 (1994).

51. Ebert, D., Vampolsky, L. & Stearns, S. C. Genetics of life history in *Daphnia magna*. 1. Heritabilities at two food levels. *Heredity* **70**, 335–343 (1993).

52. Müller, D. G. Generationswechsel, kernphasenwechsel und sexualität der braunalge *Ectocarpus siliculosus* im kulturversuch. *Planta* **75**, 39–54 (1967).

53. Buckler, E. S. et al. The genetic architecture of maize flowering time. *Science* **325**, 714–718, doi: 10.1126/science.1174276 (2009).

54. Tiley, G. P. & Burleigh, G. The relationship of recombination rate, genome structure, and patterns of molecular evolution across angiosperms. *BMC evolutionary biology* **15**, 194 (2015).

55. Zhang, N. et al. Construction of a high density SNP linkage map of kelp (Saccharina japonica) by sequencing Taq I site associated DNA and mapping of a sex determining locus. *BMC genomics* **16**, 1 (2015).
56. Ye, N. et al. Saccharina genomes provide novel insight into kelp biology. *Nature Communications* **6**, doi: 10.1038/ncomms7986 (2015).

57. Kong, A. et al. Fine-scale recombination rate differences between sexes, populations and individuals. *Nature* **467**, 1099–1103 (2010).

58. Paape, T. et al. Fine scale population recombination rates, hotspots and correlates of recombination in the *Medicago truncatula* genome. *Genome biology and evolution* evs046 (2012).

59. Gratani, L. Plant phenotypic plasticity in response to environmental factors. *Advances in botany* **2014** (2014).

60. El-Soda, M., Malosetti, M., Zwaan, B. J., Koornneef, M. & Aarts, M. G. Genotype× environment interaction QTL mapping in plants: lessons from *Arabidopsis*. *Trends in Plant Science* **19**, 390–398 (2014).

61. Orr, H. A. Testing natural selection vs. genetic drift in phenotypic evolution using quantitative trait locus data. *Genetics* **149**, 2099–2104 (1998).

62. Phillips, P. C. Epistasis—the essential role of gene interactions in the structure and evolution of genetic systems. *Nature Reviews Genetics* **9**, 855–867 (2008).

63. Cordell, H. J. Epistasis: what it means, what it doesn't mean, and statistical methods to detect it in humans. *Human molecular genetics* **11**, 2463–2468 (2002).

64. Gillespie, J. H. & Turelli, M. Genotype-environment interactions and the maintenance of polygenic variation. *Genetics* **121**, 129–138 (1989).

65. Abasht, B. et al. Fatness QTL on chicken chromosome 5 and interaction with sex. *Genetics Selection Evolution* **38**, 297–311 (2006).

66. Munns, R. & Tester, M. Mechanisms of salinity tolerance. *Annu. Rev. Plant Biol.* **59**, 651–681 (2008).

67. Bergmann, N. et al. Population-specificity of heat stress gene induction in northern and southern eelgrass Zostera marina populations under simulated global warming. *Molecular Ecology* **19**, 2870–2883 (2010).

68. Coyer, J. et al. Genomic scans detect signatures of selection along a salinity gradient in populations of the intertidal seaweed *Fucus serratus* on a 12 km scale. *Marine genomics* **4**, 41–49 (2011).

69. Goecke, F., Labes, A., Wiese, J. & Imhoff, J. F. Chemical interactions between marine macroalgae and bacteria. *Marine Ecology Progress Series* **409**, 267–299 (2010).

70. Hollants, J., Leliaert, F., De Clerck, O. & Willems, A. What we can learn from sushi: a review on seaweed-bacterial associations. *FEMS microbiology ecology* **83**, 1–16 (2013).

71. Pedersén, M. Identification of a cytokinin, 6-(3 Methyl-2-butenylamino) purine, in sea water and the effect of cytokinins on brown algae. *Physiologia Plantarum* **28**, 101–105 (1973).

72. Pedersén, M. *Ectocarpus fasciculatus*: marine brown alga requiring kinetin. (1968).

73. Catchen, J., Hohenlohe, P. A., Bassham, S., Amores, A. & Cresko, W. A. Stacks: an analysis tool set for population genomics. *Molecular Ecology* **22**, 3124–3140, doi: 10.1111/mec.12354 (2013).

74. Zhang, J., Kobert, K., Flouri, T. & Stamatakis, A. PEAR: a fast and accurate Illumina Paired-End reAd mergeR. *Bioinformatics* **30**, 614–620 (2014).

75. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. *Bioinformatics* btu170 (2014).

76. Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. *Nature Methods* **9**, 357–U354, doi: 10.1038/nmeth.1923 (2012).

77. Broman, K. W., Wu, H., Sen, S. & Churchill, G. A. R/qlt: QTL mapping in experimental crosses. *Bioinformatics* **19**, 889–890, doi: 10.1093/bioinformatics/btg112 (2003).

78. Fishman, L., Kelly, A. J., Morgan, E. & Willis, J. H. A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions. *Genetics* **159**, 1701–1716 (2001).

79. Nakagawa, S. & Schielzeth, H. Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists. *Biological Reviews* **85**, 935–956, doi: 10.1111/j.1469-185X.2010.00141.x (2010).

80. Voorrips, R. E. MapChart: Software for the graphical presentation of linkage maps and QTLs. *The Journal of Heredity* **93**, 77–78 (2002).

81. Broman, K. W. & Sen, S. *A Guide to QTL Mapping with R/qlt*. Vol. 46 (Springer, 2009).

## Acknowledgements

This work was supported by the Centre National de la Recherche Scientifique, the University Pierre and Marie Curie, the Agence Nationale de la Recherche (project Idealg ANR-10-BTBR-04-01 and ANR-12-JSV7-0008), the Brittany Region grant SAD 2015-ECTOQTL 9119 to Komlan Avia, the International Research Network "Diversity, Evolution and Biotechnology of Marine Algae" (GDRI N° 0803) and the European Research Council (ERC, grant agreement 638240). We are very grateful to Philippe Potin and Monique Ras for their management and support. We are grateful to Akira Peters for providing different strains for our experiments and for his advice. We thank Claire Dagin-Thiébaut, Charlotte Roby and Delphine Scornet for precious advice during experiments. We thank Gildas Le Corguillé for bioinformatics support. We thank Ian Probert and Jérôme Coudret for technical support. We thank Thomas Broquet for his comments on the paper. We thank the SeaDAS Development Group at NASA GSFC, Greenbelt, Maryland, the Ocean Biology Processing Group (OBPG) at NASA GSFC, Greenbelt, Maryland and Brockmann Consult GmbH, Hamburg, Germany for distributing the SeaDAS software and the MODIS data. We thank the UMR 8199 LIGAN-PM Genomics platform (Lille, France) which belongs to the 'Federation de Recherche' 3508 Labex EGID (European Genomics Institute for Diabetes; ANR-10-LABX-46) and was supported by the ANR Equipex 2010 session (ANR-10-EQPX-07-01; 'LIGAN-PM'). The LIGAN-PM Genomics platform (Lille, France) is also supported by the FEDER and the Region Nord-Pas-de-Calais-Picardie. We thank the anonymous reviewers for their comments that led to significant improvement of the manuscript.

## Author Contributions

S.M.C., J.M.C., M.V. and P.B. conceived the study, managed its organization and edited the draft versions of the manuscript. K.A. performed research (all the different steps), analyzed data and wrote the first draft of the manuscript, F.L. performed research (preliminary tests and phenotyping of the mapping family), S.M. performed research (the ddRAD-seq experiment), G.M. analyzed data for his master's internship, S.F. edited the manuscript, S.M.C. and A.C. analyzed data. All authors discussed the results, commented on the manuscript, read and approved the final manuscript.

## Additional Information

Supplementary information accompanies this paper at http://www.nature.com/srep

**Competing Interests:** The authors declare no competing financial interests.
**How to cite this article:** Avia, K. *et al*. High-density genetic map and identification of QTLs for responses to temperature and salinity stresses in the model brown alga *Ectocarpus*. *Sci. Rep.* 7, 43241; doi: 10.1038/srep43241 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017