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

HAL
open science

# Genetic Diversity in the UV Sex Chromosomes of the Brown Alga Ectocarpus

Komlan Avia, Agnieszka P Lipinska, Laure Mignerot, Alejandro Montecinos, Mahwash Jamy, Sophia Ahmed, Myriam Valero, Akira F Peters, J. Mark Cock, Denis Roze, et al.

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

## To cite this version:

Komlan Avia, Agnieszka P Lipinska, Laure Mignerot, Alejandro Montecinos, Mahwash Jamy, et al. Genetic Diversity in the UV Sex Chromosomes of the Brown Alga Ectocarpus. Genes, 2018, 9 (6), pp.286. 10.3390/genes9060286. hal-03145778

HAL Id: hal-03145778

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

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_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Genetic Diversity in the UV Sex Chromosomes of the Brown Alga *Ectocarpus*

Komlan Avia<sup>1,2,*†</sup>,<sup>1</sup> ID, Agnieszka P. Lipinska<sup>1,†</sup>, Laure Mignerot<sup>1</sup>, Alejandro E. Montecinos<sup>2,3</sup>, Mahwash Jamy<sup>1</sup>, Sophia Ahmed<sup>1</sup>, Myriam Valero<sup>2 ID</sup>, Akira F. Peters<sup>4 ID</sup>, J. Mark Cock<sup>1 ID</sup>, Denis Roze<sup>2</sup> and Susana M. Coelho<sup>1 ID</sup>

Received: 1 May 2018; Accepted: 5 June 2018; Published: 6 June 2018

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

**Abstract:** Three types of sex chromosome system exist in nature: diploid XY and ZW systems and haploid UV systems. For many years, research has focused exclusively on XY and ZW systems, leaving UV chromosomes and haploid sex determination largely neglected. Here, we perform a detailed analysis of DNA sequence neutral diversity levels across the U and V sex chromosomes of the model brown alga *Ectocarpus* using a large population dataset. We show that the U and V non-recombining regions of the sex chromosomes (SDR) exhibit about half as much neutral diversity as the autosomes. This difference is consistent with the reduced effective population size of these regions compared with the rest of the genome, suggesting that the influence of additional factors such as background selection or selective sweeps is minimal. The pseudoautosomal region (PAR) of this UV system, in contrast, exhibited surprisingly high neutral diversity and there were several indications that genes in this region may be under balancing selection. The PAR of *Ectocarpus* is known to exhibit unusual genomic features and our results lay the foundation for further work aimed at understanding whether, and to what extent, these structural features underlie the high level of genetic diversity. Overall, this study fills a gap between available information on genetic diversity in XY/ZW systems and UV systems and significantly contributes to advancing our knowledge of the evolution of UV sex chromosomes.

**Keywords:** UV sex chromosomes; pseudoautosomal region; brown algae; neutral diversity

## 1. Introduction

Morphologically distinct sex chromosomes have evolved multiple times independently in both plants and animals [1]. Sex chromosome evolution has been mainly studied in male-heterogametic (XX/XY) and female heterogametic (ZZ/ZW) sex determination systems. A typical sex chromosome pair derives from a pair of autosomes through the acquisition of genes involved in sex determination. If more than one locus involved in sex determination is located on the chromosome, recombination between these loci is expected to be suppressed, leading to the establishment of a non-recombining

<footnote n="1">¹ Sorbonne Université, UPMC Univ Paris 06, CNRS, Algal Genetics Group, Integrative Biology of Marine Models, Station Biologique de Roscoff, CS 90074, 29688 Roscoff, France; alipinska@sb-roscoff.fr (A.P.L.); laure.mignerot@sb-roscoff.fr (L.M.); mahwash.jamy@ebc.uu.se (M.J.); S.M.Ahmed@leeds.ac.uk (S.A.); cock@sb-roscoff.fr (J.M.C.); coelho@sb-roscoff.fr (S.M.C.)</footnote>

<footnote n="2">² Evolutionary Biology and Ecology of Algae, CNRS, Sorbonne Universités, UPMC, University of Paris VI, UC, UACH, UMI 3614, 29688 Roscoff, France; jano.montecinos@gmail.com (A.E.M.); myriam.valero@sb-roscoff.fr (M.V.); denis.roze@sb-roscoff.fr (D.R.)</footnote>

<footnote n="3">³ Facultad de Ciencias, Instituto de Ciencias Ambientales y Evolutivas, Universidad Austral de Chile, Casilla 567, Valdivia, Chile</footnote>

<footnote n="4">⁴ Bezhin Rosko, 29250 Santec, France; akirapeters@gmail.com</footnote>

<footnote>* Correspondence: komlan.avia@sb-roscoff.fr; Tel.: +33-256-452-138</footnote>

<footnote>† These authors contributed equally to this work.</footnote>
region on the nascent sex chromosome, the sex-determining region or SDR. The formation of this non-recombining region has important consequences for the evolution of this part of the genome. Repetitive DNA can accumulate, leading to an increase in SDR size (see review by Bachtrog [1]). There is also a tendency for genes within the SDR to degenerate as a consequence either of an accumulation of deleterious mutations or of a lower rate of adaptation [1,2]. At a later stage, deletion of non-functional DNA from within the SDR may lead to a decrease in the size of the SDR. Furthermore, the SDR can progressively expand into the flanking regions of the chromosome (the pseudoautosomal regions, PAR), so that it encompasses an increasingly greater proportion of the sex chromosome.

Evolutionary processes at a given site in a genome are influenced by selection acting on closely linked sites, an effect called Hill-Robertson interference [3,4]. This selection interference reduces the effective population size ($N_e$) experienced by the site in question [5]. As this effect is expected to be maximal in regions experiencing little or no recombination, diversity in SDRs such as Y-linked regions will be reduced compared to that of autosomes or pseudo-autosomal regions. In addition, loci on the Y chromosome are expected to experience a $N_e$ that is one-quarter that of autosomal, and one-third that of X-linked genes. Moreover, since the level of neutral polymorphism maintained at equilibrium is proportional to the product of $N_e$ and the neutral mutation rate, $\mu$ ($\pi = 4Ne\mu$) [6], diversity should be lower for Y-linked genes than for their X-linked counterparts and diversity in both X and Y genes should be lower than for autosomal genes [7]. Note, however, that this does not apply to genes located in the PAR, which should have the same $N_e$ as autosomal genes. Accordingly, in *Silene latifolia* for example, diversity has been shown to be reduced in the Y-linked regions relative to X-linked regions [8]. A similar situation was observed in *Drosophila* [9] and in Saudi-Arabian hamadryas baboons [10].

In contrast to the SDR, the PARs of sex chromosomes maintain similarity between alleles of the same gene because they undergo homologous pairing and recombination. Therefore, genes in the PARs are expected to evolve in a similar manner to autosomal genes [11,12]. However, because of the proximity to the SDR, PARs are expected to display specific evolutionary dynamics [13]. One of these specificities is that linkage with the SDR widens the conditions allowing the maintenance of polymorphism at loci under sexually antagonistic selection and also increases neutral diversity due to longer coalescence times [14–16]. Diversity in PARs is therefore predicted to be high in those regions in close proximity with the SDR [17]. Increased genetic diversity, and overall footprints of balancing selection due to sexually-antagonistic selection, have been observed for several *Silene latifolia* PAR loci, although there was no evidence for an effect of proximity to the SDR [18,19].

While information (theoretical and empirical) is available for XY and ZW sex chromosome systems, we know very little about evolutionary process, and in particular about patterns of genetic diversity, in a third type of sexual system that exists in nature, UV sex chromosomes (see review by Wilson Sayres [20]). In UV systems, which are very common in non-vascular plants and red, green and brown algae, sexes are expressed during the haploid stage of the life cycle, and females carry a U chromosome whereas males carry a V chromosome [21–23]. UV sexual systems (Figure 1) have specific evolutionary and genetic properties, including the absence of homozygous or heterozygous sexes and the absence of masking of deleterious mutations during the haploid phase when sex is expressed. Another significant feature that distinguishes UV from diploid sex chromosomes is that the $N_e$ for U- and V-specific regions is expected to be half that of autosomes [21] and both the U and the V are theoretically subject to the same mutation rate (unlike XY or ZW systems where, for example, the Y can have a higher mutation rate than the X [24,25]). The PAR region is expected to have the same $N_e$ as autosomes.
<figure><img src="image_5.png" /></figure>

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

**Figure 1.** Comparisons of XX/XY and UV sex determination systems. (A) XX/XY sex determination system. The sexual individuals are diploid, and the sex of an offspring is determined after fertilization, depending on the sex chromosome contributed by the sperm. Please note that the haploid phase of the life cycle is limited to the gametic stage. Note also that ZW systems function in a similar manner to XY systems, with diploid phase sex determination, but it is the female that is the heterogametic sex; (B) UV sex determination system. The diploid, asexual generation (sporophyte) carries both the U and the V sex chromosomes which are passed on to the haploid spores after meiosis. Spores that receive the V sex chromosome develop into a male gametophyte whereas spores carrying U sex chromosome will produce a female gametophyte. Egg and sperm produced by the gametophyte fuse to return to the diploid generation. In UV sex determination systems, the sexual individuals are haploid and sex chromosomes function in the haploid state.

The only detailed studies focusing on the structure and evolution of both the SDR and PAR regions of UV systems have been carried out in the brown algal model organism *Ectocarpus* sp. [26–28]. In this organism, the U and V-specific regions are small, and exhibit mild degeneration despite the action of haploid purifying selection [26]. SDR genes were shown to evolve rapidly, mainly due to relaxed purifying selection [29]. Remarkably, the relatively large *Ectocarpus* PARs exhibit unique features. Although they recombine normally, these regions differ from autosomes in terms of their gene density, transposable element content and genetic structure [28]. Moreover, the PAR is significantly enriched, compared to autosomes, in genes expressed specifically or predominantly during the diploid, sporophyte phase of the life cycle (hereinafter called sporophyte-biased or SP-biased genes), and these genes have been shown to evolve faster than unbiased genes [28]. A model was proposed to explain this enrichment phenomenon, giving SP-biased genes an advantage to spread when they were partially linked to the SDR and had a positive effect on fitness in one of the sexes [28]. The model assumes that the evolution of the PAR in haploid systems is under the influence of differential selection pressures in males and females acting on alleles that are advantageous during the sporophyte generation of the life cycle.

Here we used extensive Double Digest Restriction Associated DNA sequencing (ddRAD-seq) data combined with a gene-by-gene approach to perform a comprehensive analysis of the genetic diversity
across the UV sex chromosome of the brown alga *Ectocarpus*. We show that the level of neutral diversity in the U and V SDR haplotypes is about half that of the autosomes. This observation is in line with theoretical predictions based on the reduced effective population size of the SDR (the U and V SDR each have half the $N_e$ of the autosomes), suggesting that the influence of additional factors such as background selection or selective sweeps is minimal. Interestingly, genetic diversity in the PAR region was surprisingly elevated and there were several indications that genes in this region may be under balancing selection.

## 2. Materials and Methods

### 2.1. *Quantitative Trait Loci (QTL) Mapping of the Sex Locus in Ectocarpus Siliculosus*

The populations analyzed in this study belong to the species *Ectocarpus siliculosus*, a sister species to *Ectocarpus* sp. for which a reference genome sequence and genetic map with detailed coordinates of the positions of the SDR and pseudoautosomal regions are available [27,30–32]. The *Ectocarpus* sp. reference genome strain still lacks a formal species name. It was referred to as *Ectocarpus* 7 in a recent phylogenetic analysis paper [33] and we therefore also referred to it as *Ectocarpus* 7 throughout this paper. To confirm that the borders of the SDR were the same in the two species, in the absence of a complete genome sequence for *E. siliculosus*, we generated a genetic map for this species [34] focusing specifically on the sex chromosome, to investigate the location of the SDR, as described below.

A diploid sporophyte (Ec236) was generated by crossing two compatible *E. siliculosus* strains from a Naples population (EA1 and RB1) (Table S1). From this sporophyte, 152 haploid gametophytes were isolated, each arising from a unique meiotic event [35,36]. The sex of each individual was determined using sex-specific PCR markers [37]. Molecular methods, ddRAD sequencing of this population and detailed analysis of the genetic map obtained are described in the Supplementary Materials and Methods.

#### QTL Mapping of the Sex Locus

The genetic map obtained with the 152 haploid gametophyte progeny derived from the diploid sporophyte Ec236 was used to map the *E. siliculosus* SDR. Using the sex of the progeny as a binary trait, the SDR location was determined as a QTL in the R package R/qtl [38] with the scanone function and the "binary" model. To confirm its position, we also used MapQTL [39] with the Kruskal-Wallis non-parametric method. The R/xoi package (version 0.67-4) [40] was used to obtain a smoothed estimate of the recombination rate along the linkage groups (LGs), in 1 Mbp sliding windows.

### 2.2. *Sequencing of the Individuals from Different Natural Populations for Population Genomics Analyses*

#### 2.2.1. Field Sample Unialgal Collections

We selected several populations from different geographical origins in order to test the repeatability of the observed pattern of neutral variation between genomic regions in face of the population history or environment conditions of this cosmopolitan species. The samples used were previously obtained from natural populations collected along the European coast of the Atlantic and Mediterranean Sea and the Pacific coast of Chile by Couceiro, et al. [41] and Montecinos, et al. [42] (Table S2). Samples were maintained in the lab as unialgal cultures as described in Couceiro, et al. [41].

#### 2.2.2. Analyses of Neutral Diversity in Non-Recombining versus Recombining Regions

ddRAD-sequencing is a reduced-representation genome sequencing method that involves digestion of the genomic DNA with two different restriction enzymes followed by size selection, PCR amplification and sequencing of the obtained library on a sequencing platform such as Illumina MiSeq or HiSeq [43]. ddRAD-seq data were generated for 49 diploid *E. siliculosus* individuals (sporophytes), representing three European and one South American population and 6 additional haploid individuals (gametophytes) from two European and one South American population. Sample information and accession numbers are given in Table S2. Sequencing methods followed the protocols described in [30]. Sequences consisted
of paired-end reads obtained using the Illumina HiSeq 2500 system (Illumina Inc., San Diego, CA, USA). Final reads were trimmed to 70 bp. After quality (QUAL >30 and minimum genotype quality =40) and missing data filtering (not more than 40% of missing data per sample and not more than 40% of missing data per locus), 39 samples were retained for further analysis.

Quality-filtered reads were mapped to the *Ectocarpus* 7 reference genome using BWA [44] with the parameters "bwa mem -M -c 2". The genetic map of the reference *Ectocarpus* 7 genome strain provided high-quality annotation and allowed the compartmentalization of the mapped reads to autosomal, pseudoautosomal and non-recombining regions. Genotypes were called using samtools mpileup v.1.6 [45] and filtered using vcftools v.0.1.15 [46]. Only high-quality genotype calls (Phred-scaled mapping quality and genotype quality ≥20) were retained and sites with more than 25% missing data were excluded from the downstream analysis. Additionally, since coding sequences can experience positive or negative evolution that will affect their diversity patterns (see [47]), we excluded regions overlapping with exons to minimize the effect of selection and focused on neutrally evolving sites. We used samtools mpileup to report non-variant sites as well as polymorphic sites in order to concatenate the sequenced portion of the genome and assemble a "reduced" genomic sequence.

The method described above was validated by employing an alternative approach, which is described in the Supplementary Materials. Briefly, this method was based on using the Stacks pipeline [48] to carry out a *de novo* analysis of the ddRAD-seq data. The main difference between the first method described above and this alternative method was that, in the latter, only bi-allelic variant positions were called with the Stacks pipeline. Because both approaches gave consistently similar outputs, we focus here on the results obtained using the first method.

### 2.2.3. Analyses of Non-Recombining and Recombining Coding Regions

In addition to the ddRAD sequencing data, we used cDNA obtained from a total of 20 individuals (11 males and 9 females) from an *E. siliculosus* population from Naples (Italy) to provide further information about genetic diversity in the coding regions of autosomes and the sex chromosome (Table S3). Total RNA was extracted using a chloroform-isoamyl alcohol protocol (adapted from [45]). The SuperScript IV Reverse Transcription System kit (Thermo Fisher Scientific Inc., Villebon sur Yvette, France) was used to synthesize cDNA using random hexaprimers and oligo primers according to manufacturer's instructions. Primers were designed using Primer3 (Table S4). The coding regions of two autosomal, six PAR and five SDR genes (three male-specific and two female-specific) (Table 1) were amplified and sequenced on a 3130xl-3 capillary sequencer (Applied Biosystems Inc., Station Biologique de Roscoff, France). Amplicon sequences were processed with CodonCode Aligner v5.1.5 (http://www.codoncode.com).

<table><caption>**Table 1.** Selected autosomal, pseudoautosomal regions (PAR) and sex-determining region (SDR) genes analyzed in this study. In brackets after the length of the studied region, is the total length of the coding sequence.</caption><thead><tr><th>Gene</th><th>Functional Description</th><th>No. of Sequences</th><th>Length of Region Studied (Total CDS)</th><th>Segregating Sites</th><th>No. of Haplotypes</th></tr></thead><tbody><tr><td colspan="6"><strong>Autosomes</strong></td></tr><tr><td>Ec-18_002220</td><td>Alpha tubulin</td><td>19</td><td>141 (1362)</td><td>1</td><td>2</td></tr><tr><td>Ec-20_003070</td><td>Similar to G-protein coupled receptors</td><td>20</td><td>747 (1650)</td><td>1</td><td>2</td></tr><tr><td colspan="6"><strong>PAR</strong></td></tr><tr><td>Ec-13_003040</td><td>Expressed unknown protein</td><td>16</td><td>258 (450)</td><td>28</td><td>10</td></tr><tr><td>Ec-13_000140</td><td>Expressed unknown protein</td><td>19</td><td>603 (1200)</td><td>8</td><td>2</td></tr><tr><td>Ec-13_001070</td><td>Expressed unknown protein</td><td>12</td><td>390 (1728)</td><td>6</td><td>5</td></tr><tr><td>Ec-13_003030</td><td>Expressed unknown protein</td><td>9</td><td>330 (648)</td><td>5</td><td>3</td></tr><tr><td>Ec-13_004000</td><td>Tetratricopeptide TPR_2 repeat protein</td><td>13</td><td>423 (2997)</td><td>0</td><td>1</td></tr><tr><td>Ec-13_002700</td><td>Expressed unknown protein</td><td>18</td><td>606 (993)</td><td>4</td><td>4</td></tr><tr><td colspan="6"><strong>SDR</strong></td></tr><tr><td>Ec-sdr_f_000010</td><td>STE20-like serine/threonine kinase</td><td>8</td><td>597 (1265)</td><td>0</td><td>1</td></tr><tr><td>Ec-sdr_f_000090</td><td>GTPase activating protein</td><td>9</td><td>543 (2547)</td><td>0</td><td>1</td></tr><tr><td>Ec-13_001710</td><td>GTPase activating protein</td><td>11</td><td>627 (1944)</td><td>2</td><td>3</td></tr><tr><td>Ec-13_001910</td><td>STE20-like serine/threonine kinase</td><td>9</td><td>1062 (1314)</td><td>1</td><td>2</td></tr><tr><td>Ec-13_001980</td><td>Ubiquitin C-terminal hydrolase</td><td>11</td><td>570 (1110)</td><td>0</td><td>1</td></tr></tbody></table>
## 2.3. Data Analysis

We used vcftools [46] and the concatenated genomic sequences described above to calculate nucleotide diversity (π) and Tajima's *D* for non-overlapping 1-kb windows across the different genome compartments (autosomes, PAR and SDR). Since the SDR regions of both the U and the V are rather small (ca. 1Mbp), the windowed π and Tajima's *D* values calculated separately for the U and V SDRs were grouped subsequently and represented as one global (concatenated) SDR region to increase the statistical power of the downstream analyses. We performed the analysis of the *E. siliculosus* species (all populations pooled) as a whole and also separately for each population where indicated.

Nucleotide diversity and its partition between synonymous (*p*<sub>S</sub>) and non-synonymous (*p*<sub>N</sub>) mutations and synonymous (*D*<sub>S</sub>) and non-synonymous (*D*<sub>N</sub>) divergence rates for the coding regions in the population from Naples were calculated in DnaSP v5.10.01 [49] using the reference genome strain (*Ectocarpus* 7) as an outgroup. Coalescent simulations were performed with recombination (where applicable) under the standard neutral model with 10,000 iterations to obtain values and 95% confidence intervals of Tajima's *D* and Fay and Wu's *H*. The direction of selection (DoS) statistic [50] was calculated using the following formula:

$$DoS = D_N/[D_N + D_S] - p_N/[p_N + p_S]$$

We used the Wilcoxon test to compare the neutral diversity (π) and Tajima's *D* values generated per window across different genomic regions. All statistical analyses were performed in RStudio (R version 3.3.2) with graphs produced using the R package ggplot2 [51].

## 2.4. Data Availability

Data availability and accession numbers are described in Table S2. All sequences have been deposited in the Sequence Read Archive (SRA) under the ID SRP149054 (BioProject ID: PRJNA473288).

# 3. Results

## 3.1. Identification of the PAR and SDR in the Sex Chromosome of *E. siliculosus*

*E. siliculosus* is a closely related species to the *Ectocarpus* 7 genome-sequenced strain, whose assembly is of very high quality [27,31]. Knowledge of the exact position of the borders of the SDR and PAR is particularly important in the context of this study because theoretical models predict increased neutral diversity at the borders of the SDR and PAR regions [13,17]. To conduct the population genetic tests on the different genomic compartments, we needed to confirm that sister species *E. siliculosus* and *Ectocarpus* 7 (reference genome strain) shared the same PAR-SDR boundary. Previous studies had shown that the SDR of *E. siliculosus* and *Ectocarpus* 7 contain exactly the same genes [29]; however, the analysis focused on genes and did not determine whether the borders of the SDR were the same in the two species. We therefore used a newly generated genetic map for *E. siliculosus* [34] focusing specifically on the sex chromosome, to investigate the location of the SDR (Figure S1). A sex-specific QTL peak was detected at 50.9 cM on the linkage group 2 (362 markers over 104.3 cM). Based on the mapping of *de novo* assembled tags of *E. siliculosus* onto the reference genome, the SDR boundaries (position 2,775,867 bp to position 3,674,342 bp on the chromosome 13 of the reference genome of *Ectocarpus* 7 [52]) were found to be located overall at the same positions in the two species (Figure S1).

Taken together these analyses confirmed that the position of the SDR of the sex chromosome is similar in *E. siliculosus* and *Ectocarpus* 7 (Figure S1). Therefore, we concluded that the ddRAD-seq reads from *E. siliculosus* populations can be assigned to the different genomic compartments based on the *Ectocarpus* 7 genome annotation.
## 3.2. *ddRAD-seq Data*

We used *ddRAD-seq* data generated from samples collected from four populations of *E. siliculosus* (Table S2), three from Europe (Ribadeo and Gandario, in Spain and Roscoff in France) and one from South America (Pan de Azucar, Chile) to assess the extent to which genetic diversity in the male and female SDR, PAR and autosomes differed in different genomic regions.

Once demultiplexed and cleaned, the *ddRAD* sequencing generated between 14.6 and 28.9 million sequence reads per population (Table S2). Based on the number of uniquely mapped reads to the *Ectocarpus* 7 reference genome, we estimated that the data covered ca. 13% of the reference genome sequence. This proportion of the genome captured is in the range of typical proportions captured by Restriction Associated DNA (RAD) sequencing methods (e.g., [53]) and RAD data has been shown to provide useful information for analyses such as neutral diversity studies [54,55]. After applying stringent filtering (see material and methods), more than 2 million sites were scored across the genome for all individuals including 187,062 SNPs (Table S5).

## 3.3. *DNA Neutral Diversity*

Sequence diversity was estimated for the autosomes, the PAR and the male and female SDRs (Figure 2A). The windowed values were calculated separately for the U and V SDRs (Figure S2) and grouped subsequently to be represented as one global SDR region. Neutral diversity ($\pi$) was relatively similar across all autosomes (mean $\pi_A$ = $3.23 \times 10^{-3} \pm 3.47 \times 10^{-5}$ SE) (Figure S3A). Genetic diversity on the SDR was approximately half that of the autosomes (mean $\pi_{SDR}$ = $0.00221 \pm 3.97 \times 10^{-4}$ SE) and this difference was significant (Wilcoxon test, *p* = 0.0005) (Figure 2A), corresponding approximately to the equilibrium neutral expectations for the population size of these regions (the U and V SDR each have half of the $N_e$ of the autosomes). Remarkably, however, the PAR exhibited significantly higher diversity (mean $\pi_{PAR}$ = $4.39 \times 10^{-3} \pm 3.11 \times 10^{-4}$ SE; Wilcoxon test, *p* = 0.0004) (Figure 2A and Figure S3A).

<figure><img src="image_7.png" /><figcaption>**Figure 2.** Population genetics statistics for the three *Ectocarpus siliculosus* genomic compartments: autosomes, pseudoautosomal regions (PAR) and sex-determining region (SDR). (A) Boxplots of neutral diversity ($\pi$), calculated in 1 kb windows without overlap; (B) Boxplots of Tajima's *D*, calculated in 1 kb windows without overlap. The mean values are represented by the diamond shape. Letters above the boxplots denote significant differences (Wilcoxon test, *p* < 0.0005).</figcaption></figure>
Although $\pi$ values were higher overall in the PAR compared to the autosomes (Figure 2A), sliding window analysis along the sex chromosome (Figure 3A) did not show any clear bias towards elevated $\pi$ regions being located close to the SDR.

Another factor that could influence genetic diversity is gene density. Higher gene density implies more nearby sites potentially under selection and could influence the levels of diversity in the linked sites in a negative manner. Negative correlations between gene density and local levels of neutral diversity have recently been described in *Heliconius* species [47] We found a weak positive correlation between local gene density and neutral diversity along the PAR (Spearman's $rho = 0.27, p = 0.016$), but no correlation (Spearman's $rho = 0.03, p = 0.915$) for an autosomal chromosome of similar size (chr4). Positive correlation with gene density on the PAR could indicate increased polymorphism due to balancing selection. Recent work in *Silene* has shown that some of the loci on the PAR that exhibit high diversity may be under balancing selection [19], which is reflected by positive Tajima's *D* values [56]. To test the possibility that the high diversity values found in the *Ectocarpus* PAR region reflect balancing selection we therefore performed Tajima's *D* tests on the different genomic compartments (autosomes, PAR and SDRs) (Figures 2B and 3B).

Overall Tajima's *D* was negative for all autosomes (mean Tajima's *D* = -0.389 ± 0.014 SE, Figure S3B). Strikingly, however, Tajima's *D* showed positive values for the PAR that were significantly higher than values obtained for autosomes (mean Tajima's *D* = 0.130 ± 0.090 SE, Wilcoxon test, $p = 2.1 \times 10^{-8}$) (Figure 2B). Tajima's *D* was also elevated in the SDR (mean Tajima's *D* = 0.0761 ± 0.196 SE) but was not significantly different from that of the autosomes.

### 3.4. Comparison of DNA Neutral Genetic Diversity Pattern between the Four Study Populations

The neutral diversity and Tajima's *D* pattern along the sex chromosome of *E. siliculosus* was analyzed separately for each of the four study populations (Figure S4A,B respectively). The pattern of diversity in the PAR was similar despite their different geographical origin and diverse ecological environments (Figure S4A) (Wilcoxon pairwise test between populations, $p > 0.2$). The pattern of Tajima's *D* on the other hand was more variable (Figure S4B), and this result may reflect differences among population history such as different demographic processes.

### 3.5. PAR Nucleotide Diversity and Generation-Biased Genes

The PAR of *Ectocarpus* exhibits unusual characteristics compared to the autosomes, including higher transposable element content, lower gene density and higher $D_N/D_S$ rates [28,57]. Of the 455 genes on the PAR, 177 show generation-biased expression patterns (out of 6202 genome-wide [58]) and the PAR regions are significantly enriched in sporophyte-biased genes (82 compared to 2097 genome-wide) [28,57]. Interestingly, when we plotted diversity in sliding windows and compared it with the position of these generation-biased genes, the pattern of neutral diversity on the PAR was significantly correlated with the distribution of the generation-biased genes in the windows of analysis (Figure 3A) (Kendall rank correlation test: tau = 0.18, *p* = 0.02). For Tajima's *D* per window (Figure 3B), the correlation was not significant (Kendall rank correlation test: tau = -0.03, *p* = 0.7).

Based on the linkage map, the average recombination rate for the PAR was 19.2 cM/Mb. Four autosomes of similar size (chr4, chr5, chr21 and chr26) had average recombination levels in a similar range (between 17 and 20.6 cM/Mb). As shown before for *Ectocarpus* 7 [28], the PAR of *E. siliculosus* did not exhibit a significantly higher recombination rate on average than the autosomes. To investigate the possibility of a direct link between recombination rate and nucleotide diversity, we plotted recombination rate together with neutral diversity along the sex chromosome as well as for one representative autosome with a similar number of sliding windows for the $\pi$ estimates (Figure S5). We also analyzed the correlation between neutral diversity and the recombination rate for those two chromosomes (Figure S6). In the sex chromosome, recombination rate was lower around the position of the SDR (as expected) and at the one of the telomeres (Figure S5). Correlation analysis indicated a weak negative correlation ($r = -0.18$) between recombination rate and nucleotide diversity,
which was barely significant ($p = 0.049$). In the autosome, we observed more variation in the recombination rate than in the pattern of nucleotide diversity along the chromosome but no significant correlation between nucleotide diversity and recombination rate was observed (Figures S5 and S6).

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Population genetic statistics for the sex chromosome of *Ectocarpus siliculosus*. (A) Sliding window analysis of neutral diversity ($\pi$) in 1 kb non-overlapping windows along the concatenated ddRAD sequences of the sex chromosome. Values of $\pi$ are indicated by black dots; (B) Sliding window analysis of Tajima's D in 1 kb windows (indicated by black dots) along the concatenated sex chromosome sequence. The number of genes with differential expression in the sporophyte or gametophyte generation (generation-biased genes or GBGs) normalized by the physical distance covered by the concatenated 1 kb ddRAD window are marked in orange. Global gene density along the concatenated sex chromosome sequence is represented by the heatmap. The position of the sex-determining region is shaded in gray.</figcaption></figure>

### 3.6. Genetic Diversity in a Selected Subset of Autosomal, PAR and SDR Genes

ddRAD-seq data used in this study gave a broad overview of neutral diversity (because we removed regions corresponding to exons from the data before analysis) but did not provide information specifically for genes. To investigate diversity patterns at the gene level in the PAR and SDR, and in particular to test if the signal of high values of $\pi$ in the PAR measured using the ddRAD-seq data was due to footprints of balancing selection on generation-biased genes, we determined $\pi$ for a subset of six single-copy PAR genes and five single-copy SDR genes (two female-specific and three male-specific, Figure 4A) to study. For the PAR subset, we chose four generation-biased genes and two genes without significant bias in expression. We also sequenced two single-copy autosomal genes.

All 13 autosomal, PAR and SDR genes were successfully amplified from at least eight *E. siliculosus* individuals and aligned to the reference sequence from *Ectocarpus* 7 [52]. Nucleotide diversity statistics for all genes were studied based on synonymous sites (Tables 1 and 2 and Figure 4).
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 4.** Neutral diversity of PAR and SDR genes in *Ectocarpus siliculosus*. (A) Physical position of the studied PAR genes on the sex chromosome. Gene names are indicated, the sex-determining region (SDR) is marked in red (female) and blue (male); (B) Diversity of the studied autosomal, PAR and SDR genes at synonymous sites. The level of differential expression (fold change) of the PAR genes in the sporophyte (SP) and gametophyte (GA) generations is represented by the colored circles (green denotes strong GA-bias and red stands for strong SP-bias).</figcaption></figure>

<table><caption>**Table 2.** Measurements used to infer the evolutionary forces acting on autosomal, PAR and SDR genes. Significant *p*-values (*p* < 0.05) for the neutrality tests and for the bias in expression between generations are shown in boldface.</caption><thead><tr><th>Gene</th><th>π<sub>Syn</sub> #</th><th>D<sub>N</sub>/D<sub>S</sub></th><th>DoS</th><th>Tajima's D</th><th>p-Value</th><th>Fay and Wu's H</th><th>p-Value</th><th>log<sub>2</sub>FC (SP/GA) *</th></tr></thead><tbody><tr><td colspan="9"><strong>Autosomes</strong></td></tr><tr><td>Ec-18_002220</td><td>0.01138</td><td>0</td><td>No replacements</td><td>0.417</td><td>0.368</td><td>0.257</td><td>0.196</td><td>-0.52</td></tr><tr><td>Ec-20_003070</td><td>0.00106</td><td>0.13</td><td>0.316</td><td>-0.592</td><td>0.501</td><td><strong>-1.516</strong></td><td>0.015</td><td>-0.78</td></tr><tr><td colspan="9"><strong>PAR</strong></td></tr><tr><td>Ec-13_003040</td><td>0.032</td><td>0.512</td><td>-0.223</td><td><strong>-1.740</strong></td><td>0.0002</td><td>-3.467</td><td>0.102</td><td><strong>2.25</strong></td></tr><tr><td>Ec-13_000140</td><td>0.016</td><td>0.253</td><td>0.105</td><td><strong>2.033</strong></td><td>0.0007</td><td>-1.392</td><td>0.167</td><td><strong>6.52</strong></td></tr><tr><td>Ec-13_001070</td><td>0.02</td><td>0</td><td>0</td><td>-0.120</td><td>0.443</td><td>-0.848</td><td>0.171</td><td><strong>1.19</strong></td></tr><tr><td>Ec-13_003030</td><td>0.008</td><td>0.494</td><td>0.158</td><td><strong>-1.678</strong></td><td>0.013</td><td><strong>-4.278</strong></td><td>0.0073</td><td><strong>0.98</strong></td></tr><tr><td>Ec-13_004000</td><td>0</td><td>0.312</td><td>no polymorphism</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td><strong>0.71</strong></td></tr><tr><td>Ec-13_002700</td><td>0.002</td><td>0.109</td><td>-0.265</td><td>-1.381</td><td>0.050</td><td>-0.863</td><td>0.101</td><td>-0.44</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th>Gene</th><th>π<sub>Syn</sub> #</th><th>D<sub>N</sub>/D<sub>S</sub></th><th>DoS</th><th>Tajima's D</th><th>p-Value</th><th>Fay and Wu's H</th><th>p-Value</th><th>log<sub>2</sub>FC (SP/GA) *</th></tr></thead><tbody><tr><td colspan="9">SDR</td></tr><tr><td>Ec-sdr_f_000010</td><td>0</td><td>0</td><td>no polymorphism</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>-0.27</td></tr><tr><td>Ec-sdr_f_000090</td><td>0</td><td>0.6</td><td>no polymorphism</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0.41</td></tr><tr><td>Ec-13_001710</td><td>0.00351</td><td>0.692</td><td>0.190</td><td>0.199</td><td>0.399</td><td>0.255</td><td>0.412</td><td>-0.53</td></tr><tr><td>Ec-13_001910</td><td>0.00229</td><td>0.077</td><td>0.222</td><td>1.401</td><td>0.122</td><td>-0.139</td><td>0.210</td><td>0.16</td></tr><tr><td>Ec-13_001980</td><td>0</td><td>2.077</td><td>no polymorphism</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>-0.54</td></tr></tbody></table>

Table 2 and Figure 4 summarize the diversity statistics for all genes, based on synonymous sites. Consistent with our ddRAD-seq data, PAR genes tended to exhibit high neutral diversity with a mean and standard error of π = 0.013 ± 0.0046, whereas SDR genes exhibited the lowest diversity (0.0012 ± 0.0007), though it should be noted that this difference was not found to be significant by a pairwise Wilcoxon test (p = 0.28). The two autosomal genes showed an average π = 0.0062. Consistent with the ddRAD-seq data, no obvious correlation could be observed between the neutral diversity of PAR genes in relation to their distance from the SDR (Figure 4), i.e., the genes with the highest diversity were not necessarily located closer to the SDR border. We noted, however, that all five sporophyte-biased genes showed higher π values, whereas the two genes that were not SP-biased showed lower π (Figure 4). Note however, that statistical power was low for these analyses because of the limited number of genes studied.

### 3.7. Evolutionary Histories of the Selected Autosomal, PAR and SDR Genes

Most of the seven selected PAR and SDR genes that were polymorphic exhibited negative Tajima's *D* values, except for one gene in the PAR (Ec-13_000140) and two genes in the male SDR (Ec-13_001710 and Ec-13_001910) (Table 2). The positive DoS values for the latter two genes suggest that they are evolving under positive selection. However, a scenario of random differences being fixed due to the smaller effective population size (N<sub>e</sub>) of the SDR (compared to the autosomes) cannot be ruled out [59]. Out of the six PAR genes, only Ec-13_001070 appeared to be evolving neutrally. Ec-13_004000 might also be evolving neutrally but the relatively high D<sub>N</sub>/D<sub>S</sub> value of 0.312 suggests positive selection.

The PAR gene Ec-13_003030 is also most likely evolving under positive selection (or has undergone a selective sweep) given the positive DoS value, and strongly negative Tajima's *D* and Fay and Wu's *H* values. Similarly, a positive DoS value suggests that Ec-13_000140 is also under adaptive selection. However, Ec-13_000140 also has a significantly positive Tajima's *D* (*p* < 0.01) which strongly suggests that the gene may be under balancing selection, a phenomenon that cannot be detected by the DoS statistic. Please note that on the other hand, Ec-13_003040 and Ec-13_002700 seem to be evolving under relaxed purifying selection, as they have much higher polymorphism than divergence (indicated by the negative DoS value). Concerning the autosomal genes, one lacked replacements sites and hence D<sub>N</sub>/D<sub>S</sub> could not be estimated. The second one showed a low D<sub>N</sub>/D<sub>S</sub> value and none of them showed a significant Tajima's *D* value.

## 4. Discussion

### 4.1. Non-Recombining Regions in U and V Chromosomes Exhibit Reduced Neutral Diversity

Reduced neutral diversity is a hallmark of Y chromosomes [60] and other non-recombining chromosomes [61]. Reduction in neutral diversity is expected to be proportional to the effective population size of the non-recombining chromosomes compared to the autosomes, leading to a prediction of 1/4 as much neutral diversity on the Y chromosome compared to the autosomes in XY systems or 1/2 as much neutral diversity on the V or U SDRs compared to the autosomes in UV

<footnote>* SP—sporophyte, GA—gametophyte, FC—fold change in expression. # π<sub>Syn</sub> values for synonymous sites.</footnote>
sex-determination systems, assuming a 1:1 sex ratio and similar levels of reproductive success in males and females. The data presented in this study was in line with this prediction, with the neutral diversity of the combined U and V sex-determining regions being 0.68 of the autosomal chromosomes.

In several species, the observed levels of $\pi$ for the SDR have been lower than would have been predicted based solely on this neutral effect. In *Rumex*, for example, the low level diversity observed in the Y region is consistent with an important role of linked selection, with either effects of purifying selection alone or combined with positive selection driving loss of diversity [62]. Similarly, neutral diversity on the human Y chromosome [60] and on the W chromosome of birds [63] was drastically lower than that of autosomes (5–10 times and 8–13 times, respectively). These patterns of loss of diversity could be driven by selective sweeps due to sexual selection acting on testis-specific Y chromosome genes [64]; however, they can also occur in the absence of sexual selection as shown in the example of the flycatcher's W chromosome [63].

In contrast, the agreement between the theoretical prediction and the observed level of neutral diversity for *Ectocarpus* suggests that evolution of its SDR is not significantly driven by background selection or selective sweeps. In other words, since the observed pattern of diversity could be explained by neutral processes (e.g., stochastic processes caused by genetic drift), it is not necessary to propose additional evolutionary forces such as background selection or selective sweeps. This finding is congruent with our previous studies that failed to detect signatures of positive selection acting on a set of SDR genes that have been conserved across several brown algal species [29]. Please note that the *Ectocarpus* SDR is rather small with only 20 and 22 genes being sex-linked on the V and U chromosomes, respectively [26] and this species displays limited levels of sexual dimorphism [57]. In the absence of recombination, Hill-Robertson interference should decrease the local $N_e$ and thereby diversity due to effects of linked selection, but the magnitude of this effect will depend on the number of linked selected sites [62,65,66]. Therefore, a small SDR may be less affected by Hill-Robertson interference than an SDR with many genes where there is a larger scope for selection.

## 4.2. *Increased Nucleotide Diversity on the PAR Compared with Autosomes*

Our data indicate that the PAR of the sex chromosome had higher median neutral diversity than any of the autosomes. Given that the PAR is also enriched in sporophyte-biased genes, one interesting possibility is that the two phenomena are connected. Such an association is supported by the observation that all of the five genes that exhibited SP-biased expression in our gene-by-gene analysis also presented high $\pi_{syn}$ values and there was evidence that at least one of the genes was evolving under balancing selection.

Two of the five SP-biased genes (Ec-13_002700 and Ec-13_003040) also presented sex-biased gene expression and, overall, the PAR has been shown to be enriched in female-biased genes [28]. These observations open up the possibility that polymorphism may be maintained both by generation-antagonistic selection and by sexually-antagonistic selection, as has been proposed for *Silene* [19]. A recently proposed model to explain the spread of generation-biased alleles in the PAR of *Ectocarpus* [28] assumes that the evolution of the PAR in haploid systems is influenced by differential selection pressures in males and females acting on alleles that are advantageous during the sporophyte generation of the life cycle. One consequence of this model is that loci on the PAR that are subject to sexually-antagonistic selection would tend to be more polymorphic.

We have also considered alternative explanations for the elevated level of neutral diversity on the PAR. For example, we cannot exclude an effect of recombination rate on the variation of the nucleotide diversity [47, 67, 68]. Our analysis showed a weak negative correlation between recombination rate and neutral diversity, but this correlation was barely significant. No significant correlation was observed for a representative autosome. A positive correlation between DNA diversity and recombination rate has been found in many organisms including animals [69–73], plants [74–76] and fungi [77]. Although such patterns have been attributed to the action of natural selection, neutral explanations for an observed positive correlation between recombination and diversity have also
been proposed, when, in addition to intraspecific variation, recombination rate was correlated
with divergence (e.g., [78]). There are also cases where such a correlation has not been found.
For example, in maize, single nucleotide polymorphism and recombination rate were found to be
uncorrelated [79–81]. Additional studies will be required to depict more precisely the relationship
between neutral diversity and recombination rate in *Ectocarpus* genomes.

Finally, the higher diversity of the PAR could also potentially be explained by a higher mutation rate. Elevated local mutation rate (as Ds) was a strong predictor of intraspecific diversity in *Heliconius* species [47] but did not explain higher levels of π in the PAR of *Silene* [19]. In *Ectocarpus*, PAR genes have been shown to have significantly higher Ds than autosomal genes [28], so an analysis of mutation rates across the genome would be valuable to understand the role of this parameter in driving higher genetic diversity across the PAR. A thorough comparison of nucleotide variations between parents and progeny of a mapping population for instance would allow the mutation rate within the different chromosomal regions to be evaluated. Such an analysis has been carried out for the collared flycatcher [82].

4.3. *Lack of a Regional Pattern of PAR Diversity in Relation to the SDR*

Neutral sites very closely linked to any balanced polymorphism are expected to have higher diversity than surrounding genome regions [83]. Neutral diversity is thus expected to be elevated in the regions of the PAR that are closely linked to the SDR [17]. However, recombination during each generation will break down the association and the effect on diversity will become negligible unless both the recombination rate and the effective population size are small. This prediction is consistent with the observation that the peaks of polymorphism near sites known to be under long-term balancing selection are often confined to the gene itself. Therefore, under the neutral null hypothesis, only PAR genes very closely linked to the SDR should be affected [17]. Both ddRAD-seq and gene-by-gene analysis failed to find evidence for a pattern of elevated π close to the SDR border in *Ectocarpus*. The failure to find such a pattern could be due to the low resolution of ddRAD-seq markers or may be explained by a disruption of the PAR-SDR association due to gene movement and chromosomal rearrangements. The latter explanation was proposed to explain neutral diversity patterns in *Silene* [8]. Consistent with this, we showed recently that brown algal UV sex chromosomes are very dynamic with a lot of gene trafficking and rearrangements taking place [29]. Another possibility, again also suggested for *Silene* [8], is that the PAR used to be a non-recombining region, and recombination started recently. Under this hypothesis, the increased number of variants we found on the PAR could reflect a long evolutionary history of complete sex linkage. This hypothesis would also explain the intermediate characteristics of the PAR (between the autosome and SDR) in terms of GC content, TE content, gene density, Ds, gene size etc. [28]. Additional information about the PARs of more brown algal species is required to fully understand the evolutionary history of these interesting chromosomal regions.

**5. Conclusions**

We have used ddRAD-seq and gene-by-gene analyses to investigate, for the first time, patterns
of neutral diversity in a UV sex determination system. Our results confirm theoretical predictions
for UV SDRs, which predict that the level of genetic diversity in such regions should be about
half that of the autosomes. This correlation between the theoretical prediction and experimental
measurements suggests that the evolution of the SDR in *Ectocarpus* is not significantly impacted by
Hill-Robertson effects.

In contrast, the PAR exhibited a higher level of diversity than the autosomes and high levels
of Tajima's D, suggesting balancing selection. This observation is interesting because the PAR of
*Ectocarpus* has been shown to have an unusual structure with low gene density, high transposable
element content and an enrichment in both generation biased and sex-biased genes [28,57]. More work
is needed to understand whether, and to what extent, these structural features underlie the high level
of genetic diversity.
**Supplementary Materials:** The following are available online at http://www.mdpi.com/2073-4425/9/6/286/s1,
Table S1: Samples used for genetic mapping of the sex locus of *E. siliculosus*; Table S2: Natural populations
used for ddRAD sequencing and ddRAD-seq data summary; Table S3: Samples used in the analysis of the
population genetics statistics in coding regions; Table S4: PCR primers used for amplification of the PAR and
SDR genes; Table S5: Summary statistics of the VCF files; Figure S1. Mapping of the sex-determining region
in *Ectocarpus siliculosus*; Figure S2. Population genetics statistics for the three *Ectocarpus siliculosus* genomic
compartments: autosomes, PAR and SDR; Figure S3. Population genetics statistics for the individual autosomes,
PAR and SDR in *Ectocarpus siliculosus*; Figure S4. Comparison of neutral diversity patterns along the sex
chromosome of *Ectocarpus siliculosus* between the four study populations; Figure S5. The landscape of neutral
diversity (π) and recombination rates along the sex chromosome and autosomes in *Ectocarpus siliculosus*; Figure S6.
Correlation between neutral diversity and recombination rate.

**Author Contributions:** A.M., A.F.P., M.J., L.M., S.A. and K.A. prepared the biological material and performed experiments. A.L., M.J. and K.A. performed the computational analysis. A.L., K.A. and S.M.C. wrote the manuscript with valuable input from M.V., D.R. and J.M.C. S.M.C. coordinated the study. All authors read and approved the final manuscript.

**Funding:** This work was supported by the CNRS, Sorbonne Université, the project IDEALG (France: ANR-10-BTBR-04), the Brittany Region grant SAD 2015-ECTOQTL 9119 to K.A., the Becas-Chile doctoral grant to A.E.M. (CONICYT, advanced human resources program) and the ERC (grant agreement 638240). Field collections in Chile and Italy (2012/2013) were supported by two grants from the EU FP7 "capacities" specific program ASSEMBLE to A.F.P. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

**Acknowledgments:** We thank Thomas Broquet and Christophe Destombe for fruitful discussions, Dieter G. Muller for the isolation and maintenance of the *Ectocarpus siliculosus* strains. We thank Gildas Le Corguillé from the Roscoff Bioinformatic platform: ABiMS (http://abims.sb-roscoff.fr) for bioinformatics support. We thank Jérôme Coudret, Chloé Jollivet and Stéphane Mauger for technical support. Finally, we thank the two anonymous reviewers for their constructive comments, which helped us to improve the manuscript.

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

**References**

1. Bachtrog, D. Y-chromosome evolution: Emerging insights into processes of Y-chromosome degeneration. *Nat. Rev. Genet.* **2013**, *14*, 113. [CrossRef] [PubMed]

2. Charlesworth, B.; Charlesworth, D. The degeneration of Y chromosomes. *Philos. Trans. R. Soc. B Biol. Sci.* **2000**, *355*, 1563–1572. [CrossRef] [PubMed]

3. Cutter, A.D.; Payseur, B.A. Genomic signatures of selection at linked sites: Unifying the disparity among species. *Nat. Rev. Genet.* **2013**, *14*, 262. [CrossRef] [PubMed]

4. Hill, W.G.; Robertson, A. The effect of linkage on limits to artificial selection. *Genet. Res.* **1966**, *8*, 269–294. [CrossRef] [PubMed]

5. Charlesworth, B. The effects of deleterious mutations on evolution at linked sites. *Genetics* **2012**, *190*, 5–22. [CrossRef] [PubMed]

6. Kimura, M. The number of heterozygous nucleotide sites maintained in a finite population due to steady flux of mutations. *Genetics* **1969**, *61*, 893–903. [PubMed]

7. Hellborg, L.; Ellegren, H. Low levels of nucleotide diversity in mammalian Y chromosomes. *Mol. Biol. Evolut.* **2004**, *21*, 158–163. [CrossRef] [PubMed]

8. Qiu, S.; Bergero, R.; Forrest, A.; Kaiser, V.B.; Charlesworth, D. Nucleotide diversity in *Silene latifolia* autosomal and sex-linked genes. *Proc. R. Soc. Lond. B Biol. Sci.* **2010**, *277*, 3283–3290. [CrossRef] [PubMed]

9. McAllister, B.F.; Charlesworth, B. Reduced Sequence Variability on the NeoY Chromosome of *Drosophila americana* americana. *Genetics* **1999**, *153*, 221–233. [PubMed]

10. Handley, L.L.; Hammond, R.; Emaresi, G.; Reber, A.; Perrin, N. Low Y chromosome variation in Saudi-Arabian hamadryas baboons (*Papio hamadryas hamadryas*). *Heredity* **2006**, *96*, 298. [CrossRef] [PubMed]

11. Ellis, N.; Taylor, A.; Bengtsson, B.O.; Kidd, J.; Rogers, J. Population structure of the human pseudoautosomal boundary. *Nature* **1990**, *344*, 663. [CrossRef] [PubMed]

12. Toder, R.; Graves, J.A.M. CSF2RA, ANT3, and STS are autosomal in marsupials: Implications for the origin of the pseudoautosomal region of mammalian sex chromosomes. *Mamm. Genome* **1998**, *9*, 373–376. [CrossRef] [PubMed]
13. Otto, S.P.; Pannell, J.R.; Peichel, C.L.; Ashman, T.-L.; Charlesworth, D.; Chippindale, A.K.; Delph, L.F.; Guerrero, R.F.; Scarpino, S.V.; McAllister, B.F. About PAR: The distinct evolutionary dynamics of the pseudoautosomal region. *Trends Genet.* **2011**, *27*, 358–367. [CrossRef] [PubMed]

14. Charlesworth, B.; Jordan, C.Y.; Charlesworth, D. The evolutionary dynamics of sexually antagonistic mutations in pseudoautosomal regions of sex chromosomes. *Evolution* **2014**, *68*, 1339–1350. [CrossRef] [PubMed]

15. Jordan, C.Y.; Charlesworth, D. The potential for sexually antagonistic polymorphism in different genome regions. *Evolution* **2012**, *66*, 505–516. [CrossRef] [PubMed]

16. Kirkpatrick, M.; Guerrero, R.F. Signatures of sex-antagonistic selection on recombining sex chromosomes. *Genetics* **2014**, *197*, 531–541. [CrossRef] [PubMed]

17. Kirkpatrick, M.; Guerrero, R.F.; Scarpino, S.V. Patterns of neutral genetic variation on recombining sex chromosomes. *Genetics* **2010**, *184*, 1141–1152. [CrossRef] [PubMed]

18. Guirao-Rico, S.; Sánchez-Gracia, A.; Charlesworth, D. Sequence diversity patterns suggesting balancing selection in partially sex-linked genes of the plant *Silene latifolia* are not generated by demographic history or gene flow. *Mol. Ecol.* **2017**, *26*, 1357–1370. [CrossRef] [PubMed]

19. Qiu, S.; Bergero, R.; Charlesworth, D. Testing for the footprint of sexually antagonistic polymorphisms in the pseudoautosomal region of a plant sex chromosome pair. *Genetics* **2013**, *194*, 663–672. [CrossRef] [PubMed]

20. Wilson Sayres, M.A. Genetic diversity on the sex chromosomes. *Genome Biol. Evolut.* **2018**, *10*, 1064–1078. [CrossRef] [PubMed]

21. Bachtrog, D.; Kirkpatrick, M.; Mank, J.E.; McDaniel, S.F.; Pires, J.C.; Rice, W.; Valenzuela, N. Are all sex chromosomes created equal? *Trends Genet.* **2011**, *27*, 350–357. [CrossRef] [PubMed]

22. Bull, J.J. Sex chromosomes in haploid dioecy: A unique contrast to Muller's theory for diploid dioecy. *Am. Nat.* **1978**, *112*, 245–250. [CrossRef]

23. Cole, K.M.; Sheath, R.G. *Biology of the Red Algae*; Cambridge University Press: Cambridge, UK, 1990.

24. Lindblad-Toh, K.; Wade, C.M.; Mikkelsen, T.S.; Karlsson, E.K.; Jaffe, D.B.; Kamal, M.; Clamp, M.; Chang, J.L.; Kulbokas Iii, E.J.; Zody, M.C.; et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. *Nature* **2005**, *438*, 803. [CrossRef] [PubMed]

25. Makova, K.D.; Li, W.-H. Strong male-driven evolution of DNA sequences in humans and apes. *Nature* **2002**, *416*, 624. [CrossRef] [PubMed]

26. Ahmed, S.; Cock, J.M.; Pessia, E.; Luthringer, R.; Cormier, A.; Robuchon, M.; Sterck, L.; Peters, A.F.; Dittami, S.M.; Corre, E.; et al. A haploid system of sex determination in the brown alga *Ectocarpus* sp. *Curr. Biol.* **2014**, *24*, 1945–1957. [CrossRef] [PubMed]

27. Cock, J.M.; Sterck, L.; Rouze, P.; Scornet, D.; Allen, A.E.; Amoutzias, G.; Anthouard, V.; Artiguenave, F.; Aury, J.-M.; Badger, J.H.; et al. The *Ectocarpus* genome and the independent evolution of multicellularity in brown algae. *Nature* **2010**, *465*, 617–621. [CrossRef] [PubMed]

28. Luthringer, R.; Lipinska, A.P.; Roze, D.; Cormier, A.; Macaisne, N.; Peters, A.F.; Cock, J.M.; Coelho, S.M. The Pseudoautosomal Regions of the U/V Sex Chromosomes of the Brown Alga *Ectocarpus* Exhibit Unusual Features. *Mol. Biol. Evolut.* **2015**, *32*, 2973–2985. [CrossRef] [PubMed]

29. Lipinska, A.P.; Toda, N.R.; Heesch, S.; Peters, A.F.; Cock, J.M.; Coelho, S.M. Multiple gene movements into and out of haploid sex chromosomes. *Genome Biol.* **2017**, *18*, 104. [CrossRef] [PubMed]

30. Avia, K.; Coelho, S.M.; Montecinos, G.J.; Cormier, A.; Lerck, F.; Mauger, S.; Faugeron, S.; Valero, M.; Cock, J.M.; Boudry, P. High-density genetic map and identification of QTLs for responses to temperature and salinity stresses in the model brown alga *Ectocarpus*. *Sci. Rep.* **2017**, *7*, 43241. [CrossRef] [PubMed]

31. Cormier, A.; Avia, K.; Sterck, L.; Derrien, T.; Wucher, V.; Andres, G.; Monsoor, M.; Godfroy, O.; Lipinska, A.; Perrineau, M.M. Re-annotation, improved large-scale assembly and establishment of a catalogue of noncoding loci for the genome of the model brown alga *Ectocarpus*. *New Phytol.* **2017**, *214*, 219–232. [CrossRef] [PubMed]

32. Heesch, S.; Cho, G.Y.; Peters, A.F.; Le Corguille, G.; Falentin, C.; Boutet, G.; Coedel, S.; Jubin, C.; Samson, G.; Corre, E.; et al. A sequence-tagged genetic map for the brown alga *Ectocarpus siliculosus* provides large-scale assembly of the genome sequence. *New Phytol.* **2010**, *188*, 42–51. [CrossRef] [PubMed]

33. Montecinos, A.E.; Couceiro, L.; Peters, A.F.; Desrut, A.; Valero, M.; Guillemin, M.L. Species delimitation and phylogeographic analyses in the *Ectocarpus* subgroup *siliculosi* (Ectocarpales, Phaeophyceae). *J. Phycol.* **2017**, *53*, 17–31. [CrossRef] [PubMed]
34. Mignerot, L.; Avia, K.; Luthringer, R.; Lipinska, A.P.; Cock, J.M.; Coelho, S.M. Genetic and cellular characterization of parthenogenesis in the brown alga *Ectocarpus*. In preparation. 2018.

35. Coelho, S.M.; Scornet, D.; Rousvoal, S.; Peters, N.; Dartevelle, L.; Peters, A.F.; Cock, J.M. Genetic crosses between *Ectocarpus* strains. *Cold Spring Harbor Protocols* **2012**, *2012*, pdb. prot067942. [CrossRef] [PubMed]

36. Coelho, S.M.; Scornet, D.; Rousvoal, S.; Peters, N.T.; Dartevelle, L.; Peters, A.F.; Cock, J.M. How to cultivate *Ectocarpus*. *Cold Spring Harbor Protocols* **2012**, *2012*, pdb. prot067934. [CrossRef] [PubMed]

37. Lipinska, A.P.; Ahmed, S.; Peters, A.F.; Faugeron, S.; Cock, J.M.; Coelho, S.M. Development of PCR-based markers to determine the sex of kelps. *PLoS ONE* **2015**, *10*, e0140535. [CrossRef] [PubMed]

38. Broman, K.W.; Wu, H.; Sen, Ś.; Churchill, G.A. R/qtl: QTL mapping in experimental crosses. *Bioinformatics* **2003**, *19*, 889–890. [CrossRef] [PubMed]

39. Van Ooijen, J. MapQTL®5. In *Software for the Mapping of Quantitative Trait Loci in Experimental Populations*; Kyazma BV: Wageningen, The Netherlands, 2004.

40. Broman, K.W.; Kwak, I.-Y. xoi: Tools for Analyzing Crossover Interference. R Package Version 0.67-4. Available online: https://github.com/kbroman/xoi (accessed on 5 June 2018).

41. Couceiro, L.; Le Gac, M.; Hunsperger, H.M.; Mauger, S.; Destombe, C.; Cock, J.M.; Ahmed, S.; Coelho, S.M.; Valero, M.; Peters, A.F. Evolution and maintenance of haploid-diploid life cycles in natural populations: The case of the marine brown alga *Ectocarpus*. *Evolution* **2015**, *69*, 1808–1822. [CrossRef] [PubMed]

42. Montecinos, A.E.; Guillemin, M.L.; Couceiro, L.; Peters, A.F.; Stoeckel, S.; Valero, M. Hybridization between two cryptic filamentous brown seaweeds along the shore: Analysing pre-and postzygotic barriers in populations of individuals with varying ploidy levels. *Mol. Ecol.* **2017**, *26*, 3497–3512. [CrossRef] [PubMed]

43. 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* **2012**, *7*, e37135. [CrossRef] [PubMed]

44. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. *arXiv*, 2013.

45. Pearson, G.; Lago-Leston, A.; Valente, M.; Serrão, E. Simple and rapid RNA extraction from freeze-dried tissue of brown algae and seagrasses. *Eur. J. Phycol.* **2006**, *41*, 97–104. [CrossRef]

46. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T. The variant call format and VCFtools. *Bioinformatics* **2011**, *27*, 2156–2158. [CrossRef] [PubMed]

47. Martin, S.H.; Most, M.; Palmer, W.J.; Salazar, C.; McMillan, W.O.; Jiggins, F.M.; Jiggins, C.D. Natural selection and genetic diversity in the butterfly *Heliconius melpomene*. *Genetics* **2016**, *203*, 525–541. [CrossRef] [PubMed]

48. Catchen, J.; Hohenlohe, P.A.; Bassham, S.; Amores, A.; Cresko, W.A. Stacks: An analysis tool set for population genomics. *Mol. Ecol.* **2013**, *22*, 3124–3140. [CrossRef] [PubMed]

49. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. *Bioinformatics* **2009**, *25*, 1451–1452. [CrossRef] [PubMed]

50. Stoletzki, N.; Eyre-Walker, A. Estimation of the neutrality index. *Mol. Biol. Evolut.* **2010**, *28*, 63–70. [CrossRef] [PubMed]

51. Wickham, H. ggplot2. In *Wiley Interdisciplinary Reviews: Computational Statistics*; Wiley: Hoboken, NJ, USA, 2011; Volume 3, pp. 180–185.

52. *Ectocarpus siliculosus* V2. Available online: http://bioinformatics.psb.ugent.be/orcae/overview/EctsiV2 (accessed on 5 June 2018).

53. Willing, E.-M.; Hoffmann, M.; Klein, J.D.; Weigel, D.; Dreyer, C. Paired-end RAD-seq for de novo assembly and marker design without available reference. *Bioinformatics* **2011**, *27*, 2187–2193. [CrossRef] [PubMed]

54. Catchen, J.M.; Hohenlohe, P.A.; Bernatchez, L.; Funk, W.C.; Andrews, K.R.; Allendorf, F.W. Unbroken: RADseq remains a powerful tool for understanding the genetics of adaptation in natural populations. *Mol. Ecol. Resour.* **2017**, *17*, 362–365. [CrossRef] [PubMed]

55. Lowry, D.B.; Hoban, S.; Kelley, J.L.; Lotterhos, K.E.; Reed, L.K.; Antolin, M.F.; Storfer, A. Breaking RAD: An evaluation of the utility of restriction site-associated DNA sequencing for genome scans of adaptation. *Mol. Ecol. Resour.* **2017**, *17*, 142–152. [CrossRef] [PubMed]

56. Tajima, F. Statistical methods for testing the neutral mutation hypothesis by DNA polymorphism. *Genetics* **1989**, *123*, 585–595. [PubMed]

57. Lipinska, A.P.; Luthringer, R.; Roze, D.; Cormier, A.; Peters, A.F.; Cock, J.M.; Coelho, S.M. Unusual features of the pseudoautosomal region of a U/V pair of sex chromosomes. *Eur. J. Phycol.* **2015**, *50*, 68.
58. Lipinska, A.P.; Serrano-Serrano, M.; Peters, A.F.; Kogame, K.; Cock, J.M.; Coelho, S.M. Rapid turnover of life-cycle-related genes in the brown algae. *bioRxiv* **2018**. [CrossRef]

59. Charlesworth, B.; Campos, J.L. The relations between recombination rate and patterns of molecular variation and evolution in *Drosophila*. *Annu. Rev. Genet.* **2014**, *48*, 383–403. [CrossRef] [PubMed]

60. Wilson Sayres, M.A.; Lohmueller, K.E.; Nielsen, R. Natural selection reduced diversity on human Y chromosomes. *PLoS Genet.* **2014**, *10*, e1004064. [CrossRef] [PubMed]

61. Smeds, L.; Kawakami, T.; Burri, R.; Bolivar, P.; Husby, A.; Qvarnström, A.; Uebbing, S.; Ellegren, H. Genomic identification and characterization of the pseudoautosomal region in highly differentiated avian sex chromosomes. *Nat. Commun.* **2014**, *5*, 5448. [CrossRef] [PubMed]

62. Hough, J.; Wang, W.; Barrett, S.C.; Wright, S.I. Hill-Robertson interference reduces genetic diversity on a young plant Y-chromosome. *Genetics* **2017**, *207*, 685-695. [CrossRef] [PubMed]

63. Smeds, L.; Warmuth, V.; Bolivar, P.; Uebbing, S.; Burri, R.; Suh, A.; Nater, A.; Bureš, S.; Garamszegi, L.Z.; Hogner, S. Evolutionary analysis of the female-specific avian W chromosome. *Nat. Commun.* **2015**, *6*, 7330. [CrossRef] [PubMed]

64. Singh, N.D.; Koerich, L.B.; Carvalho, A.B.; Clark, A.G. Positive and purifying selection on the *Drosophila* Y chromosome. *Mol. Biol. Evolut.* **2014**, *31*, 2612–2623. [CrossRef] [PubMed]

65. Comeron, J.M.; Williford, A.; Kliman, R. The Hill-Robertson effect: Evolutionary consequences of weak selection and linkage in finite populations. *Heredity* **2008**, *100*, 19. [CrossRef] [PubMed]

66. McVean, G.A.; Charlesworth, B. The effects of Hill-Robertson interference between weakly selected mutations on patterns of molecular evolution and variation. *Genetics* **2000**, *155*, 929–944. [PubMed]

67. Kawakami, T.; Smeds, L.; Backström, N.; Husby, A.; Qvarnström, A.; Mugal, C.F.; Olason, P.; Ellegren, H. A high-density linkage map enables a second-generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution. *Mol. Ecol.* **2014**, *23*, 4035–4058. [CrossRef] [PubMed]

68. Lander, E.S.; Linton, L.M.; Birren, B.; Nusbaum, C.; Zody, M.C.; Baldwin, J.; Devon, K.; Dewar, K.; Doyle, M.; FitzHugh, W. Initial sequencing and analysis of the human genome. *Nature* **2001**, *409*, 860–921. [CrossRef] [PubMed]

69. Aguade, M.; Miyashita, N.; Langley, C.H. Reduced variation in the yellow-achaete-scute region in natural populations of *Drosophila melanogaster*. *Genetics* **1989**, *122*, 607–615. [PubMed]

70. Begun, D.J.; Aquadro, C.F. Levels of naturally occurring DNA polymorphism correlate with recombination rates in *D. melanogaster*. *Nature* **1992**, *356*, 519. [CrossRef] [PubMed]

71. Nachman, M.W. Patterns of DNA variability at X-linked loci in *Mus domesticus*. *Genetics* **1997**, *147*, 1303–1316. [PubMed]

72. Nachman, M.W.; Bauer, V.L.; Crowell, S.L.; Aquadro, C.F. DNA variability and recombination rates at X-linked loci in humans. *Genetics* **1998**, *150*, 1133–1141. [PubMed]

73. Stephan, W.; Langley, C.H. Molecular genetic variation in the centromeric region of the X chromosome in three *Drosophila ananassae* populations. I. Contrasts between the vermilion and forked loci. *Genetics* **1989**, *121*, 89–99. [PubMed]

74. Dvorák, J.; Luo, M.-C.; Yang, Z.-L. Restriction fragment length polymorphism and divergence in the genomic regions of high and low recombination in self-fertilizing and cross-fertilizing *Aegilops* species. *Genetics* **1998**, *148*, 423–434.

75. Kraft, T.; Säll, T.; Magnusson-Rading, I.; Nilsson, N.-O.; Halldén, C. Positive correlation between recombination rates and levels of genetic variation in natural populations of sea beet (*Beta vulgaris* subsp. maritima). *Genetics* **1998**, *150*, 1239–1244. [PubMed]

76. Stephan, W.; Langley, C.H. DNA polymorphism in *Lycopersicon* and crossing-over per physical length. *Genetics* **1998**, *150*, 1585–1593. [PubMed]

77. Badouin, H.; Gladieux, P.; Gouzy, J.; Siguenza, S.; Aguileta, G.; Snirc, A.; Le Prieur, S.; Jeziorski, C.; Branca, A.; Giraud, T. Widespread selective sweeps throughout the genome of model plant pathogenic fungi and identification of effector candidates. *Mol. Ecol.* **2017**, *26*, 2041–2062. [CrossRef] [PubMed]

78. Hellmann, I.; Ebersberger, I.; Ptak, S.E.; Pääbo, S.; Przeworski, M. A neutral explanation for the correlation of diversity with recombination rates in humans. *Am. J. Hum. Genet.* **2003**, *72*, 1527–1535. [CrossRef] [PubMed]
79. Tenaillon, M.I.; Sawkins, M.C.; Anderson, L.K.; Stack, S.M.; Doebley, J.; Gaut, B.S. Patterns of diversity and recombination along chromosome 1 of maize (Zea mays ssp. mays L.). *Genetics* **2002**, *162*, 1401–1413. [PubMed]

80. Tenaillon, M.I.; Sawkins, M.C.; Long, A.D.; Gaut, R.L.; Doebley, J.F.; Gaut, B.S. Patterns of DNA sequence polymorphism along chromosome 1 of maize (Zea mays ssp. mays L.). *Proc. Natl. Acad. Sci. USA* **2001**, *98*, 9161–9166. [CrossRef] [PubMed]

81. Tenaillon, M.I.; U'ren, J.; Tenaillon, O.; Gaut, B.S. Selection versus demography: A multilocus investigation of the domestication process in maize. *Mol. Biol. Evolut.* **2004**, *21*, 1214–1225. [CrossRef] [PubMed]

82. Smeds, L.; Qvarnström, A.; Ellegren, H. Direct estimate of the rate of germline mutation in a bird. *Genome Res.* **2016**, *26*, 1211–1218. [CrossRef] [PubMed]

83. Hudson, R.R.; Kaplan, N.L. The coalescent process in models with selection and recombination. *Genetics* **1988**, *120*, 831–840. [PubMed]

<figure><img src="image_11.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/).