Boudjema Imarazene

Kang Du

Séverine Beille

Elodie Jouanno

Romain Feron

Qiaowei Pan

Jorge Torres-Paz

Céline Lopez-Roques

Adrien Castinel

Lisa Gil

Elodie Jouano

Claire Kuchly

Cécile Donnadieu

Hugues Parrinello

Laurent Journot

Cédric Cabau

Margot Zahm

Christophe Klopp

Tomáš Pavlica

Ahmed Al-Rikabi

Thomas Liehr

Sergey Simanovsky

Joerg Bohlen

Alexandr Sember

Julie Perez

Frédéric Veyrunes

Thomas D Mueller

John H Postlethwait

Manfred Schartl

Amaury Herpin

Sylvie Rétaux

INRAE Yann Guiguen
email: yann.guiguen@inrae.fr

A supernumerary "B-sex" chromosome drives male sex determination in the Pachón cavefish, Astyanax mexicanus

Keywords: Sex determination, cavefish, B chromosome, sex chromosomes, gdf6

niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

SUMMARY

Sex chromosomes are generally derived from a pair of classical type-A chromosomes, and relatively few alternative models have been proposed up to now [START_REF] Wright | How to make a sex chromosome[END_REF][START_REF] Furman | Sex Chromosome Evolution: So Many Exceptions to the Rules[END_REF] . B chromosomes (Bs) are supernumerary and dispensable chromosomes with non-Mendelian inheritance found in many plant and animal species [START_REF] D'ambrosio | B-chrom: a database on B-chromosomes of plants, animals and fungi[END_REF][START_REF] Jones | New species with B chromosomes discovered since 1980[END_REF] , that have often been considered as selfish genetic elements that behave as genome parasites [START_REF] Camacho | CHAPTER 4 -B Chromosomes[END_REF][START_REF] Camacho | B-chromosome evolution[END_REF] . The observation that in some species Bs can be either restricted or predominant in one sex [START_REF] Beladjal | B-chromosomes and male-biased sex ratio with paternal inheritance in the fairy 30 shrimp Branchipus schaefferi (Crustacea, Anostraca)[END_REF][8][START_REF] Néo | Altitudinal variation for B chromosome frequency in the characid fish Astyanax scabripinnis[END_REF][START_REF] Vicente | Sex-ratio distortion associated with the presence of a B chromosome in Astyanax scabripinnis (Teleostei, Characidae)[END_REF][START_REF] Stange | Supernumerary b chromosomes restricted to males in astyanax scabripinnis ( pisces , characidae )[END_REF][START_REF] Clark | Dynamic Sequence Evolution of a Sex-Associated B Chromosome in Lake Malawi Cichlid Fish[END_REF][START_REF] Clark | Changing sex for selfish gain: B chromosomes of Lake Malawi cichlid fish[END_REF][START_REF] Yoshida | B Chromosomes Have a Functional Effect on Female Sex Determination in Lake Victoria Cichlid Fishes[END_REF] raised the interesting hypothesis that Bs could play a role in sex determination [START_REF] Camacho | B Chromosomes and Sex in Animals[END_REF] . The characterization of putative B master sex-determining (MSD) genes, however, has not yet been provided to support this hypothesis. Here, in Astyanax mexicanus cavefish originating from Pachón cave, we show that Bs are strongly male-predominant. Based on a high-quality genome assembly of a B-carrying male, we characterized the Pachón cavefish B sequence and found that it contains two duplicated loci of the putative MSD gene growth differentiation factor 6b (gdf6b). Supporting its role as an MSD gene, we found that the Pachón cavefish gdf6b gene is expressed specifically in differentiating male gonads, and that its knockout induces male-to-female sex reversal in B-carrying males. This demonstrates that gdf6b is necessary for triggering male sex determination in Pachón cavefish. Altogether these results bring multiple and independent lines of evidence supporting the conclusion that the Pachón cavefish B is a "B-sex" chromosome that contains duplicated copies of the gdf6b gene which can promote male sex determination in this species.

RESULTS AND DISCUSSION

Pachón cavefish B chromosomes are male-predominant B chromosomes Supernumerary B chromosomes (Bs) are generally thought to arise from the duplication and assembly of A chromosome sequences [START_REF] Hanlon | B Chromosomes in the Drosophila Genus[END_REF][START_REF] Martis | Selfish supernumerary chromosome reveals its origin as a mosaic of host genome and organellar sequences[END_REF][START_REF] Valente | Origin and evolution of B chromosomes in the cichlid fish Astatotilapia latifasciata based on integrated genomic analyses[END_REF][START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF] and their relationship to sex chromosomes has often been suspected and discussed [START_REF] Camacho | B Chromosomes and Sex in Animals[END_REF] . Some hypotheses state that Bs are derived from sex chromosomes or, alternatively, evolved to become sex chromosomes [START_REF] Camacho | B Chromosomes and Sex in Animals[END_REF][START_REF] Pansonato-Alves | Single origin of sex chromosomes and multiple origins of B chromosomes in fish genus Characidium[END_REF][START_REF] Serrano-Freitas | Satellite DNA content of B chromosomes in the characid fish Characidium gomesi supports their origin from sex chromosomes[END_REF][START_REF] Conte | Origin of a Giant Sex Chromosome[END_REF][START_REF] Zhou | Deciphering neo-sex and B chromosome evolution by the draft genome of Drosophila albomicans[END_REF] . Because Bs have been described in A. mexicanus [START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF][START_REF] Piscor | Microsatellite Organization in the B Chromosome and A Chromosome Complement in Astyanax (Characiformes, Characidae) Species[END_REF][START_REF] Kavalco | Molecular Cytogenetics of Blind Mexican Tetra and Comments on the Karyotypic Characteristics of Genus Astyanax (Teleostei, Characidae)[END_REF] , we performed cytogenetic analyses in 17 males and 11 females of a laboratory population of Pachón cavefish (Pachón) to investigate whether Pachón Bs could be sex-restricted. We found that Pachón Bs are euchromatic mitotically unstable microchromosomes (Figure 1A-B, Figure S1), that are present in one to three copies in most male metaphases (mean number ± SD of Bs per male metaphase = 1.08 ± 0.41), contrasting with a barely detectable B occurrence in female metaphases (mean number ± SD of Bs per female metaphase = 0.05 ± 0.08, Figure 1C and Data S1A). Chromosomal mapping by fluorescence in situ hybridization (FISH) using probes generated from microdissected male Bs painted Bs in males and even the rare B in females, supporting that Pachón male-predominant Bs and the low occurrence female Bs share a similar gDNA content (Figure 1D). In addition, weaker FISH signals were also detected on different terminal parts of some A chromosomes (see small white arrows in Figure 1D), suggesting that Pachón Bs are made up of many duplicated fragments of A chromosomes [START_REF] Hanlon | B Chromosomes in the Drosophila Genus[END_REF][START_REF] Martis | Selfish supernumerary chromosome reveals its origin as a mosaic of host genome and organellar sequences[END_REF][START_REF] Valente | Origin and evolution of B chromosomes in the cichlid fish Astatotilapia latifasciata based on integrated genomic analyses[END_REF][START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF] , and / or that they share repetitive DNAs with the A chromosomes [START_REF] Serrano-Freitas | Satellite DNA content of B chromosomes in the characid fish Characidium gomesi supports their origin from sex chromosomes[END_REF][START_REF] Piscor | Microsatellite Organization in the B Chromosome and A Chromosome Complement in Astyanax (Characiformes, Characidae) Species[END_REF][START_REF] Ebrahimzadegan | Repetitive DNA landscape in essential A and supernumerary B chromosomes of Festuca pratensis Huds[END_REF] . Female or male sex-restricted or sex-predominant Bs have been described in many fishes, like for instance in some cichlids [START_REF] Clark | Dynamic Sequence Evolution of a Sex-Associated B Chromosome in Lake Malawi Cichlid Fish[END_REF][START_REF] Clark | Changing sex for selfish gain: B chromosomes of Lake Malawi cichlid fish[END_REF][START_REF] Yoshida | B Chromosomes Have a Functional Effect on Female Sex Determination in Lake Victoria Cichlid Fishes[END_REF] and characiforms 8,[START_REF] Néo | Altitudinal variation for B chromosome frequency in the characid fish Astyanax scabripinnis[END_REF][START_REF] Mizoguchi | Macro-and Microchromosomes B in Females of Astyanax scabripinnis (Pisces, Characidae)[END_REF][START_REF] Portela-Castro | New occurrence of microchromosomes B in Moenkhausia sanctaefilomenae (Pisces, Characidae) from the Paraná River of Brazil: analysis of the synaptonemal complex[END_REF] . In A. mexicanus, Bs were recently described as being restricted only in some, but not all, males [START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF] , suggesting population differences in the frequency and sex-linkage of Bs that have also been reported in many species [START_REF] Jones | New species with B chromosomes discovered since 1980[END_REF][START_REF] Camacho | CHAPTER 4 -B Chromosomes[END_REF][START_REF] Camacho | B-chromosome evolution[END_REF][START_REF] Houben | B Chromosomes -A Matter of Chromosome Drive[END_REF] . In addition to metaphase spreads (Figure 1E-F), we also detected Bs in pachytene chromosome spreads from testes of two Pachón males (N=52 and N=25).

Single Bs were found in most cells (44/52 and 23/25) and always as unpaired chromosomes (Figure 1G-H). Hence, Pachón Bs are present in meiotic germ cells and they do not pair with A chromosomes in line with what has been found in other species with Bs [START_REF] Jones | Chromosomes with a life of their own[END_REF] . The question of whether Pachón Bs can pair to each other remains open because we detected no case of multiple Bs in these pachytene chromosome spreads.

Characterization of the Pachón cavefish B chromosome sequence

Because the publicly available Pachón genome assembly was obtained from a female [START_REF] Mcgaugh | The cavefish genome reveals candidate genes for eye loss[END_REF] , we sequenced the genome of a B+ Pachón male to assemble its B sequence. Bs are notoriously difficult to assemble [START_REF] Ahmad | The Modern View of B Chromosomes Under the Impact of High Scale Omics Analyses[END_REF][START_REF] Blavet | Sequence of the supernumerary B chromosome of maize provides insight into its drive mechanism and evolution[END_REF] , due to their complex mosaic composition of A chromosome fragments [START_REF] Hanlon | B Chromosomes in the Drosophila Genus[END_REF][START_REF] Martis | Selfish supernumerary chromosome reveals its origin as a mosaic of host genome and organellar sequences[END_REF][START_REF] Valente | Origin and evolution of B chromosomes in the cichlid fish Astatotilapia latifasciata based on integrated genomic analyses[END_REF][START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF] and their high-repeat content [START_REF] Martis | Selfish supernumerary chromosome reveals its origin as a mosaic of host genome and organellar sequences[END_REF][START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF] , and most of the B sequence information is from short-read sequencing of purified Bs [START_REF] Martis | Selfish supernumerary chromosome reveals its origin as a mosaic of host genome and organellar sequences[END_REF][START_REF] Valente | Origin and evolution of B chromosomes in the cichlid fish Astatotilapia latifasciata based on integrated genomic analyses[END_REF] or B+ versus B-devoid (B-) individuals [START_REF] Martis | Selfish supernumerary chromosome reveals its origin as a mosaic of host genome and organellar sequences[END_REF][START_REF] Valente | Origin and evolution of B chromosomes in the cichlid fish Astatotilapia latifasciata based on integrated genomic analyses[END_REF][START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF][START_REF] Navarro-Domínguez | Protein-coding genes in B chromosomes of the grasshopper Eyprepocnemis plorans[END_REF] .

To accurately assemble a high-quality Pachón B sequence, we used a combination of HiFi PacBio and Oxford nanopore long-reads, 10X genomics Illumina short-linked reads, and a Hi-C chromosome contact map. The resulting whole-genome assembly (see assembly metrics in Data S1B) contains 25 large scaffolds corresponding to the 25 Pachón A-chromosome pairs [START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF][START_REF] Kavalco | Molecular Cytogenetics of Blind Mexican Tetra and Comments on the Karyotypic Characteristics of Genus Astyanax (Teleostei, Characidae)[END_REF] and 170 remaining unplaced scaffolds (2.12% of the total assembly size).

To identify the Pachón B sequence we used its male-predominant feature and a poolsequencing (pool-seq) approach to contrast whole-genome sequences of a genomic DNA (gDNA) pool of 91 phenotypic males versus a gDNA pool of 81 phenotypic females. By remapping these male and female pool-seq reads on our male Pachón genome assembly, we identified a single 2.97 Mb unplaced scaffold (HiC_Scaffold_28) that displays a clear malebiased read coverage profile (Figure 2A-A'). The sequence analysis of HiC_Scaffold_28 revealed that it is made up from a complex mosaic of numerous duplicated fragments of A chromosomes (Figure 2B) including complete but also truncated duplicates of A chromosome genes (Table S3). The B also displays a repeat content that is markedly different from A chromosomes (Figure 2C-D, Figures S2A-B, Data S1D). Both its sex-biased profile and its sequence characteristics indicate that HiC_Scaffold_28 is the Pachón B. The contribution of many A chromosome regions to the structure of Pachón Bs is in line with the recent findings that A. mexicanus Bs contain a large number of transposable elements [START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF] . The high B proportion of satellite DNA (Figures S2 and Data S1D ), was also reported in other species [START_REF] Ruiz-Ruano | Quantitative sequence characterization for repetitive DNA content in the supernumerary chromosome of the migratory locust[END_REF] . Our manually curated B gene annotation (Data S1C) identified 63 genes on the Pachón B. Of these, 20 show high-quality annotation over their full length, 11 are truncated compared to their conserved homologs in other fish and one is a chimeric gene. Five genes show multiple copies and constitute almost one-third of the B gene content (Data S1C). These results contrast with earlier studies which reported a much higher number of B genes in A. mexicanus [START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF] . These annotation differences are likely due to indirect assessment of the B gene content based on short-read sequencing of very few B+ versus B-individuals [START_REF] Ahmad | B chromosomes of multiple species have intense evolutionary dynamics and accumulated genes related to important biological processes[END_REF] . This comparison clearly illustrates the need for better, complete and high-quality B assemblies like we provided here for Pachón cavefish, to better understand the structure and gene content of Bs generally.

The Pachón cavefish B contains two copies of a putative master sex determining gene

Among the B genes with well-supported annotation evidence (Data S1C), we identified two duplicated loci of the A-chromosome-3 growth differentiation factor 6b (gdf6b) gene (located at Chr03:863,919-866,170), that is the teleost ohnolog (teleost whole genome duplication [START_REF] Christoffels | Fugu Genome Analysis Provides Evidence for a Whole-Genome Duplication Early During the Evolution of Ray-Finned Fishes[END_REF][START_REF] Vandepoele | Major events in the genome evolution of vertebrates: Paranome age and size differ considerably between ray-finned fishes and land vertebrates[END_REF] paralogous copy) of the gdf6a gene (Figure 3A). Gdf6 genes belong to the TGF-β superfamily within which many master sex-determining (MSD) genes have been found, including TGF-β receptors [START_REF] Feron | Characterization of a Y-specific duplication/insertion of the anti-Mullerian hormone type II receptor gene based on a chromosome-scale genome assembly of yellow perch, Perca flavescens[END_REF][START_REF] Kamiya | A trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, Takifugu rubripes (fugu)[END_REF][START_REF] Rafati | Reconstruction of the birth of a male sex chromosome present in Atlantic herring[END_REF] and ligands [START_REF] Li | A Tandem Duplicate of Anti-Müllerian Hormone with a Missense SNP on the Y Chromosome Is Essential for Male Sex Determination in Nile Tilapia, Oreochromis niloticus[END_REF][START_REF] Hattori | A Y-linked anti-Müllerian hormone duplication takes over a critical role in sex determination[END_REF][START_REF] Pan | Identification of the master sex determining gene in Northern pike (Esox lucius) reveals restricted sex chromosome differentiation[END_REF][START_REF] Reichwald | Insights into Sex Chromosome Evolution and Aging from the Genome of a Short-Lived Fish[END_REF] . Of note, gdf6a on the Y chromosome (gdf6aY) has been characterized as the master sex-determining gene in the turquoise killifish, Nothobranchius furzeri [START_REF] Reichwald | Insights into Sex Chromosome Evolution and Aging from the Genome of a Short-Lived Fish[END_REF] . The Pachón B gdf6b genes (B-gdf6b = B-gdf6b-1 and B-gdf6b-2) were thus retained as potential candidate B MSD genes. Chromosome FISH hybridization with a gdf6b locus probe revealed a gdf6b hybridization signal on the Pachón Bs, along with a single gdf6b labelling on a single A chromosome pair, both in the male-predominant Bs and the low occurrence female Bs (Figure 1E-F, and inset in Figure 1E).

The two B-gdf6b loci are 99.6 % identical in the 21.6 kb region shared by the two genes, and 100% identical in their coding sequences (CDS), with B-gdf6b-2 being derived from an internal B duplication of the B-gdf6b-1 locus (Figure S2C). Such internal B duplications / insertions indicate that the origin of the B structure can be more complex than initially thought.

Comparison of these B-gdf6b loci with the overlapping sequence of their A chromosome counterpart (A-gdf6b) revealed numerous differences in their proximal promoters and also their intron that contains two B-gdf6b specific insertions (Figure S2D). However, differences within the gdf6b CDS were limited to a T-to-C (A-gdf6b-to-B-gdf6b) synonymous substitution (c.591T>C) and two nonsynonymous substitutions, i.e., a T-to-G (A-gdf6b-to-B-gdf6b) transversion (c.180T>G) in exon 1 that switches a A-Gdf6b lysine into a B-Gdf6b asparagine (p.Lys60Asn) and a G-to-A (A-gdf6b-to-B-gdf6b) transition (c.679G>A) in exon 2 that switches the A-Gdf6b serine into a B-Gdf6b glycine (p.Ser227Gly) (Figure S2E). The Lys60Asn does not impact a conserved amino-acid position of Gdf6b proteins, in contrast to the Ser227Gly that impacts a glycine of the TGF-β/BMP propeptide domain that is conserved in Pachón B-Gdf6b and in all vertebrate Gdf6 proteins, but not in Pachón A-Gdf6b (Figure 3B). It is interesting to note that this Ser227Gly modification suggests that the A-gdf6b acquired this mutation after the B-gdf6b copy was duplicated on the B chromosome. This non-conserved Ser227 engages in a central hydrogen bond network at a looptip region in the TGF-β/BMP propeptide domain (Figure 3C) that could affect the stability of this pro-domain. The Gly227Ser exchange in Pachón A-Gdf6b leads to the gain of several hydrogen bonds that are due to the side chain hydroxyl group of the Ser227 residue (Figure 3D-E). The gain of hydrogen bonds can confer a higher folding stability of the proprotein complex and thereby might affect activation of Gdf6b as this requires release of the mature C-terminal domain from the proprotein complex. The mature C-terminal growth factor domain, however, is likely to be unaffected by this mutation. Whether these conformation differences between the A-Gdf6b and B-Gdf6b proproteins could provide a potential functional explanation for a sex-determining role of the male-predominant Pachón B-Gdf6b remains to be explored, but point mutations in other MSD genes like amhr2Y in Takifugu rubripes and amhY in Oreochromis niloticus, have been described to be directly responsible for male sex determination [START_REF] Kamiya | A trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, Takifugu rubripes (fugu)[END_REF][START_REF] Li | A Tandem Duplicate of Anti-Müllerian Hormone with a Missense SNP on the Y Chromosome Is Essential for Male Sex Determination in Nile Tilapia, Oreochromis niloticus[END_REF] .

Based on these B-gdf6b and A-gdf6b loci differences, we developed several B-specific PCR genotyping tests on fin clips (Figure S3A). In our Pachón laboratory population, we found a complete (100%) association between B-specific amplifications and the male phenotype in 723 males, with all the 787 tested females being negative (p-value of association with sex < 2.2e-16). We also found the same complete association in wild-caught Pachón individuals (20 males and 20 females, recognized by external sex-specific traits without sacrifice; p-value of association with sex = 1.87e-09) (Data S1F) showing that this male-predominant B is not the result of a domestication effect. These results strengthen our cytological observations of a malepredominant B. The absence of B-specific amplifications in females despite the cytogenetic detection of rare Bs in females is probably the result of a PCR sensitivity issue as increasing the number of PCR cycles allows the detection of a faint PCR fragment in Pachón females (Figure S3B).

Evidence supporting gdf6b as a potential master sex determining gene in Pachón cavefish

Sex-specific expression patterns during the sex differentiation period and alteration of gonadal development upon knockout are key arguments for the evaluation of a candidate MSD gene.

Due to the high sequence identity of the B-gdf6b and A-gdf6b cDNAs, we were not able to specifically quantify the B-gdf6b expression. Quantification of the expression of gdf6b (B-gdf6b and/or A-gdf6b) showed that it has both a predominant expression in developing and adult male gonads (as well as in male brain, intestine, kidney and swim-bladder; Figure S4A) and a strong sexually dimorphic expression during early differentiation of B+ individuals (Figure 4A). Using in situ hybridization, the expression of gdf6b during the early differentiating period was restricted to gonads of B+ individuals (at 15, 21, 30, and 60 days post-fertilization), with no strict colocalization with the gonadal soma-derived factor gene (gsdf) (Figure 4B-C and Figure S4B). Gsdf is a well-known gonad-restricted, somatic supporting cell lineage marker [START_REF] Gautier | The Proximal Promoter Region of the Zebrafish gsdf Gene Is Sufficient to Mimic the Spatio-Temporal Expression Pattern of the Endogenous Gene in Sertoli and Granulosa Cells1[END_REF][START_REF] Nakamura | Identification and lineage tracing of two populations of somatic gonadal precursors in medaka embryos[END_REF][START_REF] Sawatari | A novel transforming growth factor-β superfamily member expressed in gonadal somatic cells enhances primordial germ cell and spermatogonial proliferation in rainbow trout (Oncorhynchus mykiss)[END_REF] that has also been described as one of the earliest Pachón gonadal sex differentiation marker genes [START_REF] Imarazene | Primordial Germ Cell Migration and Histological and Molecular Characterization of Gonadal Differentiation in Pachón Cavefish Astyanax mexicanus[END_REF] . This result demonstrates that gdf6b mRNA (B-gdf6b and/or A-gdf6b) has an expression profile compatible with a male (B+) MSD function, being expressed in the right place, i.e., only in the differentiating testis, and at the right time, i.e., during early testicular differentiation.

To bring additional and functional evidence that gdf6b could act as an MSD gene, we generated gdf6b knockouts in Pachón cavefish using the genome-editing CRISPR-Cas9 system with two guide RNAs in order to remove a large part of the gdf6b CDS (Figure 4D). This large deletion (~470 bp) includes most of the TGF-β propeptide region and the beginning of the TGF-β like domain resulting in a truncated, likely non-functional, Gdf6b protein (Figure S2F). Among 200 first-generation microinjected individuals (which are mosaic for the genome-edited loci), eighteen B+ individuals (i.e., genotypic males) had a ~470 bp deletion in their A-and/or B-gdf6b exon 2 (Figure S3C), and they were all sex-reversed into phenotypic females (Figure 4G). In contrast, all B+ males and B-females without the gdf6b deletion developed normal testes (B+) or ovaries (B-) (Figure 4E-F). This result shows that gdf6b is necessary to trigger Pachón testicular development in B+ individuals and brings further functional evidence that gdf6b could be used as a male MSD gene in Pachón cavefish. However, because of the high similarity of the B-gdf6b loci with the A-gdf6b locus, we have not been able to specifically knockout the Pachón B-gdf6b. Further studies will bring more functional proof supporting the role of the B-gdf6b genes in sex determination, including the specific B-gdf6b knockout in B+ individuals and the overexpression of gdf6b by transgenesis in B-fish.

Conclusions

Altogether our results bring new pieces of evidence to support a role of some B chromosomes in sex determination. The potential implication of Bs in sex determination had been suspected in many species including some fishes [START_REF] Yoshida | B Chromosomes Have a Functional Effect on Female Sex Determination in Lake Victoria Cichlid Fishes[END_REF][START_REF] Camacho | B Chromosomes and Sex in Animals[END_REF] . Up to now, however, the characterization of potential B MSD genes along with strong functional evidence has only been provided for the bacterial-derived haploidizer gene in the jewel wasp, Nasonia vitripennis B chromosome (named PSR for paternal sex ratio chromosome) [START_REF] Benetta | Genome elimination mediated by gene expression from a selfish chromosome[END_REF] , although maleness in this haplo-diploid organism is determined through elimination of the paternal A chromosome set [START_REF] Nur | A "selfish" B chromosome that enhances its transmission by eliminating the paternal genome[END_REF] . Here, combining a variety of approaches we discovered that Pachón cavefish of the species A. mexicanus carry male-predominant Bs that contain two copies of the gdf6b gene, which itself behaves as an excellent MSD candidate gene. This indeed brings the interesting hypothesis that these Pachón Bs could be considered as "B-sex" chromosomes. However, the question remains open whether Pachón Bs are predominant in males because they are eliminated from female tissues or whether they are by themselves necessary and sufficient to trigger maleness. Despite being male-predominant, Pachón Bs are also found in some females albeit only in very few metaphases (21.6 times less abundant than in males on average) and most often as a single B copy. B frequencies have been described as being highly variable between species, sexes, individuals of the same population, and even in different cells of a single individual [START_REF] Jones | New species with B chromosomes discovered since 1980[END_REF][START_REF] Camacho | CHAPTER 4 -B Chromosomes[END_REF][START_REF] Camacho | B-chromosome evolution[END_REF][START_REF] Houben | B Chromosomes -A Matter of Chromosome Drive[END_REF] . This variation is assumed to result from meiotic and/or mitotic instability of Bs that can be present only in some organs and absent from others [START_REF] Jones | New species with B chromosomes discovered since 1980[END_REF][START_REF] Houben | B Chromosomes -A Matter of Chromosome Drive[END_REF][START_REF] Jones | Transmission and Drive Involving Parasitic B Chromosomes[END_REF][START_REF] Bernardino | B Chromosome Variants of the Grasshopper Xyleus discoideus angulatus Are Potentially Derived from Pericentromeric DNA[END_REF][START_REF] Stevens | Unstable B-chromosomes in a European population of Allium schoenoprasum L. (Liliaceae)[END_REF] . In the plant Aegilops speltoides, a mechanism of programmed elimination of Bs occurs specifically in the roots [START_REF] Ruban | Supernumerary B chromosomes of Aegilops speltoides undergo precise elimination in roots early in embryo development[END_REF] . It results from a B chromatid nondisjunction during mitosis, leading to the micronucleation of Bs and their subsequent degradation at early stages of the proto-root embryonic tissue differentiation [START_REF] Ruban | Supernumerary B chromosomes of Aegilops speltoides undergo precise elimination in roots early in embryo development[END_REF] . Such a mechanism would potentially explain a specific B elimination in Pachón female organs, as it has been hypothesized in another Astyanax species with male-restricted Bs [START_REF] Stange | Supernumerary b chromosomes restricted to males in astyanax scabripinnis ( pisces , characidae )[END_REF] . Further studies are now needed to better understand this sex-specific B drive mechanism, and if it reflects a cause or a consequence of sex determination in Pachón cavefish. Our results also lay a high-quality genome-based foundation in an important emerging fish model for studying the genomic evolution of Bs, including the micro-and macro-evolution of this B chromosome in line with the evolution of sex chromosomes. E) that were often visible on male metaphases, cannot be interpreted as the two different B-gdf6b loci due to their genomic proximity (see Table S3) and the FISH resolution. G-H:

STAR METHODS

RESOURCE AVAILABILITY

Lead contact

Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Yann Guiguen (yann.guiguen@inrae.fr).

Materials availability

To request Pachón fish lines or constructs created in this study, please contact the lead contact.

Data and code availability

Raw sequences and the whole genome assembly of Pachón cavefish have been deposited in the National Center for Biotechnology Information DDBJ/ENA/GenBank databases under the BioProject PRJNA734455. This accession number is listed in the key resources table. This study did not generate new unique code. Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

EXPERIMENTAL MODEL AND SUBJECT DETAILS

Cavefish breeding and sampling

Laboratory stocks of A. mexicanus Pachón cavefish were obtained in 2004 from the Jeffery laboratory at the University of Maryland, College Park, MD. Fish were raised as previously described. Fertilized eggs were provided by CNRS cavefish experimental facilities (Gif sur Yvette, France) and maintained at 24°C until the hatching stage occurring around 24 ± 2 hours post-fertilization (hpf) [START_REF] Hinaux | A developmental staging table for Astyanax mexicanus surface fish and Pachón cavefish[END_REF] . Subsequently, larvae were transferred and raised in the Fish Physiology and Genomics laboratory experimental facilities (LPGP, INRAE, Rennes, France) under standard photoperiod (12 h light / 12 h dark) and at two different temperatures: 21 ± 1 °C and 28 ± 1°C. Animals were fed twice a day, firstly with live artemia (Ocean Nutrition) until 15 days post-fertilization (dpf), then with a commercial diet (BioMar) until adult stage.

For animal dissections and organ sampling, fish were euthanized with a lethal dose of tricaine methanesulfonate (MS 222, 400 mg/l), supplemented by 150 mg/l of sodium bicarbonate.

Phenotypic sex of individuals was determined at 4 months and more, either by macroscopical examination of the gonads when they were enough differentiated, or by histology when gonads were not totally differentiated [START_REF] Imarazene | Primordial Germ Cell Migration and Histological and Molecular Characterization of Gonadal Differentiation in Pachón Cavefish Astyanax mexicanus[END_REF] . Caudal fin clips were collected from all individuals and stored in ethanol 90% at 4°C before genomic DNA (gDNA) extraction. For the chromosome contact map (Hi-C), 80 µl of blood was sampled from three males using a syringe rinsed with EDTA 2%. The fresh blood was slowly frozen in a Freezing Container (Mr. Frosty, Nalgene®) after addition of 15% of dimethyl sulfoxide (DMSO). Karyotypic analyses were carried out in 17 In the field, the phenotypic sex of animals was determined by checking the presence or absence of denticles on the anal fins as described previously [START_REF] Elipot | Astyanax transgenesis and husbandry: how cavefish enters the laboratory[END_REF] , a small fin clip was gently taken and fish were rapidly returned to their natural pond. In addition, pictures of each individual sampled were taken to confirm the phenotypic sex, back in the laboratory, based on the morphological criteria described previously [START_REF] Elipot | Astyanax transgenesis and husbandry: how cavefish enters the laboratory[END_REF] . The permits for field sampling (02241/13, 02438/16, 05389/17 and 1893/19) were delivered by the Mexican authorities (Mexican Secretaría del Medio 16 Ambiente y Recursos Naturales) to Sylvie Rétaux and Patricia Ornelas-Garcia (UNAM, Mexico).

METHOD DETAILS

DNA extraction

For fish genotyping, gDNA was extracted from fin clips stored in 90% ethanol, after lysis in 5% chelex [START_REF] Gharbi | A linkage map for brown trout (Salmo trutta): chromosome homeologies and comparative genome organization with other salmonid fish[END_REF] and 10 mg Proteinase K at 55 °C for 2 h, followed by 10 min at 99 °C. Following extraction, samples were centrifuged and the supernatant containing the gDNA was transferred in clean tubes and stored at -20 °C. For pool-sequencing and TaqMan assay, gDNA was extracted using NucleoSpin Kits for Tissue (Macherey-Nagel, Duren, Germany) according to the supplier's recommendations. For long-read male genome sequencing, high molecular weight (HMW) gDNA was extracted from a mature testis grounded in liquid nitrogen and lysed in TNES-Urea buffer (TNES-Urea: 4M urea; 10 mM Tris-HCl, pH 7.5; 120 mM NaCl; 10 mM EDTA; 5% SDS) for two weeks at room temperature. For HMW gDNA extraction the TNES-Urea solution was supplemented with Proteinase K at a final concentration of 150 µg/ml and incubated at 37 °C overnight. HMW gDNA was extracted with a modified phenol-chloroform protocol as previously described [START_REF] Pan | Identification of the master sex determining gene in Northern pike (Esox lucius) reveals restricted sex chromosome differentiation[END_REF] . The gDNA concentrations for both pool-seq and genome sequencing were quantified with Qubit3 fluorometer (Invitrogen, Carlsbad, CA) and HMW gDNA quality and purity were assessed using spectrophotometry, fluorometry, and capillary electrophoresis.

Primers and probe design

All primers used in this study including PCR genotyping, qPCR gene expression and cDNA cloning were designed using Primer3web software [START_REF] Kõressaar | Primer3_masker: integrating masking of template sequence with primer design software[END_REF]

B polymerase chain reaction (PCR) genotyping

Genetic sex of Pachón cave individuals was determined by PCR tests using the fact that males have a B-predominant chromosome with three different sets of primers (see Table S6 for primer sequences) based on differences between the A-gdf6b and B-gdf6b loci. Three sets of primers (P) were designed (see Figure S3A) with two primer sets designed to amplify specifically the two B-gdf6b copies based either on a single base variation between the A/B gdf6b CDS at position 679 bp (P1-P2), or based on primers located on both sides of the A/B breakpoints downstream of the B-gdf6b gene (P3-P4). The third set of primers (P5-P6) was designed on gaps/indels variations between A-gdf6b and the two B-gdf6b genes in the proximal promoter of gdf6b genes. Another set of primers (P9-P10) was designed as a PCR positive control with primers located on both sides of the A/B breakpoints downstream of the A-gdf6b gene. Primer

10× Genomics sequencing

10X Chromium Library was prepared according to 10X Genomics protocols using the Genome Reagent Kits v2. Optimal performance has been characterized on input gDNA with a mean length greater than 50 kb (~144Kb). GEM reactions were performed on 0.625 ng of genomic DNA, and DNA molecules were partitioned and amplified into droplets to introduce 16-bp partition barcodes. GEM reactions were thermally cycled (30 °C for 3 h and 65 °C for 10 min; held at 4 °C) and after amplification, the droplets were fractured. P5 and P7 primers, read 2, and sample index were added during library construction. The library was amplified using 10 cycles of PCR and the DNA was subsequently size selected to 450 bp by performing a double purification on AMPure Xp beads. Library quality was assessed using a Fragment Analyzer and quantified by qPCR using the Kapa Library Quantification Kit. Sequencing has been performed on an Illumina HiSeq3000 using a paired-end read length of 2x150 bp with the Illumina HiSeq3000 sequencing kits.

Oxford nanopore genome sequencing

High molecular weight gDNA purification steps were performed using AMPure XP beads (Beckman Coulter). Library preparation and sequencing were performed using Oxford Nanopore (Oxford Nanopore Technologies) Ligation Sequencing Kit SQK-LSK109 according to manufacturer's instructions "1D gDNA selecting for long reads (SQK-LSK109)''. Five µg of DNA was purified then sheared at 20 kb using the megaruptor1 system (diagenode). A onestep DNA damage repair + END-repair + dA tail of double-stranded DNA fragments was performed on 2 µg of sample. Then adapters were ligated to the library. Library was loaded onto 1 R9.4.1 flowcell and sequenced on a PromethION instrument at 0.02 pM within 72 h.

PacBio HiFi genome sequencing

Library preparation and sequencing were performed according to the manufacturer's instructions "Procedure & Checklist Preparing HiFi SMRTbell Libraries using SMRTbell Express Template Prep Kit 2.0". Fifteen µg of DNA were purified and then sheared at 15 kb using the Megaruptor3 system (Diagenode). Using SMRTbell Express Template prep kit 2.0, a Single strand overhangs removal and then a DNA and END damage repair steps were performed on 10 µg of the sample. Then blunt hairpin adapters were ligated to the library. The library was treated with an exonuclease cocktail to digest unligated DNA fragments. A size selection step using a 12 kb cutoff was performed on the BluePippin Size Selection system (Sage Science) with "0.75% DF Marker S1 3-10 kb Improved Recovery" protocol. Using Binding kit 2.0 kit and sequencing kit 2.0, the primer V2 annealed and polymerase 2.0 bounded library was sequenced by diffusion loading onto 2 SMRTcells on Sequel2 instrument at 50 pM with a 2 h pre-extension and a 30 h movie.

Hi-C sequencing

Hi-C data was generated using the Arima-HiC kit (Ref. Image analyses and base calling were performed using the Illumina NovaSeq Control Software and Real-Time Analysis component (v3.4.4). Demultiplexing was performed using Illumina's conversion software (bcl2fastq v2.20). The quality of the raw data and potential contaminants was assessed using FastQC (v0.14.0) [START_REF] Wingett | FastQ Screen: A tool for multi-genome mapping and quality control[END_REF] from the Babraham Institute and the Illumina software SAV (Sequencing Analysis Viewer).

Genome assembly

Pacbio HiFi reads were assembled with hifiasm 60 version 0.9 using standard parameters. The genome assembly fasta file was extracted from the principal gfa assembly graph file using an awk command line. This assembly was then scaffolded using Hi-C and 10X as a source of linking information. 10X reads were aligned using Long Ranger v2.1.1 (10x Genomics). Hi-C reads were aligned to the draft genome using Juicer 61 with default parameters. A candidate assembly was then generated with 3D de novo assembly (3D-DNA) pipeline [START_REF] Dudchenko | De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds[END_REF] with the -r 0 and --polisher-input-size 100000 parameters. Finally, the candidate assembly was manually reviewed using the Juicebox assembly tools [START_REF] Durand | Juicebox Provides a Visualization System for Hi-C Contact Maps with Unlimited Zoom[END_REF] . Due to the specific structure of the Pachón cave B chromosome, both Hi-C and 10X signals show some uncertainties in the order and orientation of the contigs. To improve the quality of the B chromosome assembly, ONT reads were then aligned to the final version of the genome using minimap2 64 v2.11 with -x map-ont parameter. Both reads spanning contig junctions and reads showing supplementary alignments linking contigs belonging to the B chromosome were analyzed to resolve these ambiguities.

Genome annotation

The cavefish whole genome assembly was annotated using a pipeline adapted from previous studies [START_REF] Powell | Natural hybridization reveals incompatible alleles that cause melanoma in swordtail fish[END_REF][START_REF] Du | The sterlet sturgeon genome sequence and the mechanisms of segmental rediploidization[END_REF] . In brief, RepeatModeler, RepeatProteinMask, and RepeatMasker (open-4.0.7, http://www.repeatmasker.org/) were first used to scan the genome and mask out repeats. Then protein-coding genes were annotated by collecting gene evidence from homology alignment, RNA-seq mapping, and ab initio prediction. For homology alignment, 464,144 protein sequences collected from NCBI were aligned to the assembly using Genewise 67 and Exonerate respectively. RNA-seq data were mapped to the assembly in two independent parallel steps.

First Hisat 68 was used to align RNA-seq reads and then StringTie 69 was used to predict the gene models; in the other step reads were first assembled into transcript sequences using Trinity [START_REF] Grabherr | Full-length transcriptome assembly from RNA-Seq data without a reference genome[END_REF] and then PASA [START_REF] Haas | Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies[END_REF] was used to map the transcripts to the assembly and model the gene structures.

For ab initio prediction and final integrating, Augutus 72 was first trained using the high-quality gene models and then ran in a hint-guide model.

B chromosome annotation

First, repeats were identified and masked from the B chromosome using RepeatModeler, RepeatProteinMask, and RepeatMasker (open-4.0.7, http://www.repeatmasker.org/). To annotate protein-coding genes, we collected all protein sequences of A. mexicanus annotated by NCBI (Genome ID: 13073) and Ensembl (release-104), and aligned them onto the repeatmasked B sequence using GeneWise [START_REF] Birney | GeneWise and Genomewise[END_REF][START_REF] She | GenBlastA: enabling BLAST to identify homologous gene sequences[END_REF] and Exonerate respectively. For each query, the best hit was kept. To determine the best gene model when multiple ones compete for a splice site, we introduced RNA-seq data to evaluate the quality of these homology gene models. Hence RNA-seq data of A. mexicanus from the previous study [START_REF] Pasquier | Gene evolution and gene expression after whole genome duplication in fish: the PhyloFish database[END_REF] were aligned to B using HISAT [START_REF] Kim | HISAT: a fast spliced aligner with low memory requirements[END_REF] and parsed using StringTie [START_REF] Pertea | StringTie enables improved reconstruction of a transcriptome from RNA-seq reads[END_REF] . RNA-score of each homology gene model was then calculated as the match-extend of splice sites to that of StringTie prediction. When multiple homology gene models compete for a splice site, those with lower RNA-score were discarded. In cases when some genes failed to be identified using homology alignment, we also implemented an ab initio gene prediction using Augustus [START_REF] Stanke | AUGUSTUS: ab initio prediction of alternative transcripts[END_REF] , where all the homology and transcriptome evidence were used as hints. The predicted results were included into the final gene set if 1) the splice sites are not occupied and 2) the splice sites match 100% to that of StringTie predictions (RNAscore =100). To further evaluate the quality of the final gene set, we blasted their protein sequences to SWISSPROT (https://www.uniprot.org/) and NR (https://www.ncbi.nlm.nih.gov/), and took the alignment to the best hit to check how much of the query and subject was aligned, respectively (query coverage & subject coverage). Genes with query and subject coverage both >90% were considered as being of good quality.

To characterize the A chromosome content of the B chromosome, sequences of the B chromosome assembly were aligned to the sequences of the 25 A chromosomes with minimap2 64 (v2.11) and the best match of each contig fragment was retained. Overlapping matches were manually filtered considering match lengths and similarities (cigarline and edit distance) in order to build the best non-overlapping matching list. Karyoplots were then plotted using the R package karyoploteR [START_REF] Gel | karyoploteR: an R/Bioconductor package to plot customizable genomes displaying arbitrary data[END_REF] (https://bernatgel.github.io/karyoploter_tutorial/). The median, minimum, and maximum sizes of the 628 B best matches on A chromosomes were respectively 1,087 bp, 44 bp and 41,908 bp.

Male and female Pool-sequencing

DNA was collected from 91 phenotypic Pachón cave males and 81 phenotypic Pachón cave females and was pooled as male and female pools separately. Before pooling, the DNA concentration was normalized in order to obtain an equal amount of each individual genome in

Chromosome conventional cytogenetics

Mitotic or meiotic chromosome spreads were obtained either from regenerating caudal fin tissue as previously described [START_REF] Völker | Direct Chromosome Preparations from Embryos and Larvae[END_REF] , with slight modifications [START_REF] Sember | Karyotype differentiation in 19 species of river loach fishes (Nemacheilidae, Teleostei): extensive variability associated with rDNA and heterochromatin distribution and its phylogenetic and ecological interpretation[END_REF] and altered time of fin regeneration (one week), or by direct preparation from the cephalic kidney and gonads [START_REF] Ràb | Cold-blooded vertebrates[END_REF] . In the latter, the quality of chromosomal spreading was enhanced by a previously described dropping method [START_REF] Bertollo | Direct Chromosome Preparation from Freshwater Teleost Fishes[END_REF] .

Chromosomes were stained with 5% Giemsa solution (pH 6.8) (Merck, Darmstadt, Germany)

for conventional cytogenetic analyses, or left unstained for other methods. For FISH, slides were dehydrated in an ethanol series (70%, 80% and 96%, 3 min each) and stored at -20 °C before analysis. Constitutive heterochromatin was visualized by C-banding [START_REF] Haaf | An early stage of ZW/ZZ sex chromosome differentiation in Poecilia sphenops var. melanistica (Poeciliidae, Cyprinodontiformes)[END_REF]

Chromosome microdissection and FISH mapping

Twelve copies of B chromosome from A. mexicanus male individual (male 10) and twelve copies encompassing two B chromosomes (per cell) from a male individual (male 9) were manually microdissected as previously described [START_REF] Ahmed | Glass needle-based chromosome microdissection-How to set up probes for molecular cytogenetics? Video[END_REF] under an inverted microscope (Zeiss Axiovert 135) using a sterile glass needle attached to a mechanical micromanipulator (Zeiss).

The chromosomes were subsequently amplified by degenerate oligonucleotide primed-PCR (DOP-PCR) following previously described protocols [START_REF] Yang | Generation of Paint Probes by Flow-Sorted and Microdissected Chromosomes[END_REF] . One µl of the resulting amplification product was used as a template DNA for a labeling DOP-PCR reaction, with Spectrum Orange-dUTP and Spectrum Green-dUTP, for Male 9 (2B) and Male 10 (1B), respectively (both Vysis, Downers Grove, USA). The amplification was done in 30 cycles, following previously described protocols [START_REF] Yang | Animal Probes and ZOO-FISH[END_REF] . Depending on the experimental scheme, the final probe mixture contained i) both painting probes (200 ng each) or ii) a single painting probe (200 ng) and a labeled 4,368 bp long fragment containing gdf6b gene and its promoter (300 ng; see below).

To block the shared repetitive sequences, the probe also contained 4-5 µg of unlabelled competitive DNA prepared from female gDNA (on male preparations) or male gDNA (on female preparations). Male and female gDNAs were isolated from liver and spleen using MagAttract HMW DNA kit (Qiagen) and C0t-1 DNA (i.e., fraction of gDNA enriched with highly and moderately repetitive sequences) was then generated from them according to previously described protocols [START_REF] Zwick | A rapid procedure for the isolation of C0t-1 DNA from plants[END_REF] . The complete B probe mixture was dissolved in the final volume 20 µl (in case of two painting probes) or 14 µl (in case of one painting probe and a gdf6b gene probe) of hybridization mixture (50% formamide and 10% dextran sulfate in 2× SSC).

FISH and whole-chromosome painting

The FISH (Fluorescence in situ hybridization on chromosomes) experiments were done using a combination of two previously published protocols [START_REF] Sember | Karyotype differentiation in 19 species of river loach fishes (Nemacheilidae, Teleostei): extensive variability associated with rDNA and heterochromatin distribution and its phylogenetic and ecological interpretation[END_REF][START_REF] Yano | Highly conserved Z and molecularly diverged W chromosomes in the fish genus Triportheus (Characiformes, Triportheidae)[END_REF] , with slight modifications. Briefly, the aging of slides took place overnight at 37 °C and then 60 min at 60 °C, followed by treatments with RNase A (200 µg/ml in 2× SSC, 60-90 min, 37 °C) (Sigma-Aldrich) and then pepsin (50 µg/ml in 10 mM HCl, 3 min, 37 °C). Subsequently, the slides were incubated in 1% formaldehyde in PBS (10 min) to stabilize the chromatin structure. Denaturation of chromosomes was done in 75% formamide in 2× SSC (pH 7.0) (Sigma-Aldrich) at 72 °C, for 3 min. The hybridization mixture was denatured for 8 min (86 °C) and then pre-hybridized at 37 °C for 45 min to outcompete the repetitive fraction. After application of the probe cocktail on the slide, the hybridization took place in a moist chamber at 37 °C for 72 h. Subsequently, non-specific hybridization was removed by post-hybridization washes: two times in 1× SSC (pH 7.0) (65 °C, 5 min each) and once in 4× SSC in 0.01% Tween 20 (42°C, 5 min). Slides were then washed in PBS (1 min), passed through an ethanol series, and mounted in antifade containing 1.5 µg/ml DAPI (Cambio, Cambridge, United Kingdom).

Synaptonemal complex immunostaining

Pachytene chromosome spreads from two Pachón cave A. mexicanus males (male 2 and 4)

were prepared from testes following the protocol for Danio rerio [START_REF] Kochakpour | Immunofluorescent microscopic study of meiosis in zebrafish[END_REF][START_REF] Kochakpour | Sex-specific crossover patterns in Zebrafish ( Danio rerio )[END_REF] with some modifications.

Briefly, dissected testes were suspended in 200-600 µl (based on the testes size and cell density) of cold PBS. Cell suspensions were applied onto poly-l-lysine slides (ThermoFisher), with 1:30 (v/v) dilution in hypotonic solution (PBS: H2O, 1:2 v/v). After 20 min at room temperature (RT), slides were fixed with freshly prepared cold 2% formaldehyde (pH 8.0 -8.5) for 3 min at RT. Slides were then washed three times in 0.1% Tween-20 (pH 8.0-8.5), 1 min each, and left to dry (1 h). Afterwards, immunofluorescence analysis of synaptonemal complexes took place, using antibodies against the proteins SYCP3 (lateral elements of synaptonemal complexes) and MLH1 (mismatch repair protein; marker for visualization of recombination sites). The primary antibodies -rabbit anti-SYCP3 (1:300; Abcam, Cambridge, UK) and mouse anti-MLH1 (1:50, Abcam) -were diluted (v/v) in 3% BSA (bovine serum albumin) in 0.05% Triton X-100/ PBS. After application onto the slides, the incubation was carried out overnight in a humid chamber at 37°C. Next day, slides were washed three times in 0.1% Tween-20 in PBS, 10 min each and secondary antibodies, diluted (v/v) in 3% BSA (bovine serum albumin) in 0.05% Triton X-100/ PBS, were applied. Specifically, we used goat antirabbit Alexa 488 (1:300; Abcam) and goat anti-mouse Alexa555 (1:100; Abcam), and the slides were incubated for 3 h at 37 °C. Then, after washing in 0,1% Tween-20 in PBS (10 min) and brief washing in 0.01% Tween-20 in distilled H2O, slides were mounted in antifade containing DAPI, as described above.

Microscopy and image analysis

At least 50 metaphase spreads per individual were analyzed to confirm the diploid chromosome number (2n), karyotype structure, and FISH results. Giemsa-stained preparations were analyzed under Axio Imager Z2 microscope (Zeiss, Oberkochen, Germany), equipped with an automatic Metafer-MSearch scanning platform. Photographs of the chromosomes were captured under 100× objective using CoolCube 1 b/w digital camera (MetaSystems, Altlussheim, Germany). The karyotypes were arranged using Ikaros software (MetaSystems, Altlussheim, Germany). Chromosomes were classified according to their centromere positions [START_REF] Levan | Nomenclature for Centromeric Position on Chromosomes[END_REF] , modified as metacentric (m), submetacentric (sm), subtelocentric (st), or acrocentric (a). FISH preparations were inspected using an Olympus BX53 epifluorescence microscope (Olympus, Tokyo, Japan), equipped with an appropriate fluorescence filter set.

Black-and-white images were captured under 100× objective for each fluorescent dye with a cooled DP30BW CCD camera (Olympus) using Olympus Acquisition Software. The digital images were then pseudo-colored (blue for DAPI, red for Cy3, green for FITC) and merged in DP Manager (Olympus). Composed images were then optimized and arranged using Adobe Photoshop CS6.

Phylogeny and synteny of Gdf6 proteins

Phylogeny and synteny relationships of Gdf6 genes were inferred from a Genomicus instance [START_REF] Louis | Genomicus: five genome browsers for comparative genomics in eukaryota[END_REF] in which synteny and phylogeny have been reconciled with the Scorpios pipeline 93 on 4 tetrapod species, one sarcopterygii species, one holostei species and 13 teleosts (see Figure 3 for species names).

Three-dimensional protein modelling

Three-dimensional models for the proprotein of A. mexicanus A-Gdf6b and B-Gdf6b were obtained by homology modeling. The amino acid sequences of full-length A-Gdf6b and B-Gdf6b were submitted to automated homology modeling using the hm_build macro of the software package YASARA [START_REF] Krieger | YASARA View-molecular graphics for all devices-from smartphones to workstations[END_REF] . Modeling includes alignment of the two target sequences against sequences from Uniprot via PSI-Blast to build a position-specific scoring function matrix/profile, which is then used to search the PDB structure data bank for suitable modeling targets. Five templates were identified with sufficiently high scores, i.e. PDB entries 4YCG

(proprotein complex of GDF2 (alternative naming BMP9)) [START_REF] Mi | Structure of bone morphogenetic protein 9 procomplex[END_REF] , 5HLY (proprotein complex of Activin A) [START_REF] Wang | Structure and activation of pro-activin A[END_REF] , 6Z3J (the mature C-terminal growth factor domain of GDF5 in complex with repulsive guidance molecule B) [START_REF] Malinauskas | Repulsive guidance molecules lock growth differentiation factor 5 in an inhibitory complex[END_REF] as well as 6Z3M (which is like 6Z3J but includes neogenin in complex with GDF5 and the repulsive guidance molecule B) [START_REF] Malinauskas | Repulsive guidance molecules lock growth differentiation factor 5 in an inhibitory complex[END_REF] , and 3QB4 (which is only the C-terminal growth factor domain of GDF5) [START_REF] Klammert | GDF-5 can act as a context-dependent BMP-2 antagonist[END_REF] . Several initial models were built on the basis of these template structures, missing sequence elements were modeled in an automated procedure through YASARA on the basis of an indexed PDB structure database. By this scheme 13 models were built, three only covered the C-terminal growth factor domain comprising residues 293 to 398, while 10 models covered the proprotein complex consisting of residues 3 to 398 (models on the basis of the 5HLY entry) and of residues 66 to 398 (models on the basis of template 4YCG). These models were individually refined by a short molecular dynamics simulation in explicit water to optimize hydrogen bonding and protein packing. Due to the overall low Z-score of the individual models, YASARA used the various models to form a hybrid homology model combining the elements with the highest-scoring factors into a single 3D model. This hybrid model was then used for further analysis.

Expression analysis by Real-Time PCR

For gene expression studies, mRNA transcripts levels were quantified during gonadal development from 10 dpf to 90 dpf, male and female gametogenesis stages, and finally in 10 adult organs including gonads as described previously [START_REF] Imarazene | Primordial Germ Cell Migration and Histological and Molecular Characterization of Gonadal Differentiation in Pachón Cavefish Astyanax mexicanus[END_REF] . All samples were frozen in liquid nitrogen and stored at -80 °C until RNA extraction. Total RNA extraction from gonads, trunks, and adult organs, followed by cDNA synthesis, and expression analysis by RT-PCR were carried out as previously described [START_REF] Imarazene | Primordial Germ Cell Migration and Histological and Molecular Characterization of Gonadal Differentiation in Pachón Cavefish Astyanax mexicanus[END_REF] . Specific primers were designed for gdf6b in the most divergent sequence regions between the two paralogous gdf6a and gdf6b genes.

RNAScope in situ hybridization of gdf6b

RNA in situ hybridization (ISH) assays have been carried out using the RNAscope®

Knockout of Pachón A. mexicanus gdf6b

Pachón cave A. mexicanus inactivated for gdf6b were generated using the CRISPR/Cas9 method. Guide RNAs (sgRNAs) targeting two sites located in exon 2 of gdf6b were designed using ZiFiT software [START_REF] Sander | ZiFiT (Zinc Finger Targeter): an updated zinc finger engineering tool[END_REF] S3A, Data S1E) flanking the sgRNAs target sites leading to a ~470 bp deletion on the exon 2 of the gdf6b genes (Figure S10). The genetic sex of the mutants was determined by specific primers (P5-P6) on the gdf6b promoter with gap/indel variation between A-gdf6b and B-gdf6b (see STAR methods above and Figure S3A).

Histology

Gonads were fixed in Bouin's fixative solution for 48 h and then dehydrated serially in aqueous 70% and 95% ethanol, ethanol/butanol (5∶95), and butanol. Tissues were embedded in paraffin blocks that were cutted serially into 5 µm sections, and were stained with hematoxylin-eosinsafran (HES) (Microm Microtech, Brignais, France).

QUANTIFICATION AND STATISTICAL ANALYSIS

Statistical analyses

For the sex genotyping marker based on the heterozygous and specific site of the B chromosome on the exon 2 (position 679 bp of the B-gdf6b CDS), the significance of the correlation between this polymorphism and the male phenotypic sex was tested with the

Figure 1 :

 1 Figure 1: Karyological characterization of male-predominant supernumerary B chromosomes (Bs) in Pachón cave Astyanax mexicanus. A-B: Representative C-banding patterns of single B+ male (A) and B-female (B) from Pachón cave. The Bs (black arrow) lack C-bands, suggesting that Pachón cavefish Bs are largely euchromatic. See also Figure S1. C: Boxplots of the average number of Bs per metaphase in Pachón cavefish males and females. Horizontal lines indicate the median, the box indicates the interquartile range (IQR), and the whiskers the range of values that are within 1.5 x IQR. Statistical significance between males and females were tested with the Wilcoxon Rank Test (*** = P < 0.001). See also Data S1A. D: Fluorescence in situ hybridization (FISH) of male Pachón cave mitotic metaphase labeled with combined microdissected B probes. Yellow arrows point to the strong labelling of Bs and the small white arrows to the lighter labelling of some different parts of A chromosomes. E-F: FISH co-labelling of a 1B male (E) and a 1B female (F) metaphase with microdissected B (green) and gdf6b-specific (red) probes. Bs are indicated by yellow arrows and the white arrowheads point to pairs of A chromosome sister chromatids labelled by the gdf6b probe. Only one A-chromosome gdf6b signal was detected in panel E. The two gdf6b signals (see inset in

  Synaptonemal complex (SC) analysis showing that Pachón cave Bs (yellow arrow) do not pair with the 25 fully synapsed standard bivalents of A chromosomes. SCs were visualized by anti-SYCP3 antibody (green), the recombination sites were identified by anti-MLH1 antibody (red) and chromosomes were counterstained by DAPI (blue). G: Merged image. H: SYCP3 visualization only. Scale bars = 5 µm.

Figure 2 :

 2 Figure 2: Genomic characterization of Pachón cavefish B chromosome (B). A: Read depth ratio of male and female Pachón genomic pools showing a strong coverage bias in a single scaffold, Hi_scaffold_28 (enlarged in A' inset showing male and female read coverage). B: Karyoplots of the A chromosome regions duplicated on the Pachón cavefish B (Chr B) showing that the Pachón B is made from a complex mosaic of duplicated A chromosomal fragments. C-D: Comparison of the repeat landscapes of the Pachón B (C) and whole-genome including the B (D), showing that the Pachón B has a very different repeat element (color code provided in inset of panel C) content compared to A chromosomes. Short interspersed repeated sequences (SINEs), long interspersed nuclear elements (LINEs), long terminal repeats (LTRs), DNA repeat elements (DNAs), and terminal inverted repeat sequences (TIRs). See also Figure S2A-B and Data S1D for additional details.

Figure 3 :

 3 Figure 3: Gdf6b protein evolution and structure. A: Simplified phylogeny (left panel) and synteny (middle panel) relationships of the gdf6a and gdf6b genes, with the additional B-gdf6b Pachón cavefish duplication, showing that they are duplicated paralogues stemming from the teleost whole genome duplication (Ts3R). Species names are given on the right side of panel B. B: Corresponding multiple alignments of Gdf6 protein-coding sequences around the A-Gdf6b lysine to B-Gdf6b asparagine switch in Exon 1 (p.Lys60Asn) and the A-Gdf6b serine to B-Gdf6b glycine switch in Exon 2 (p.Ser227Gly). C: Ribbon plot homology model of A-Gdf6b proprotein dimer (bottom panel view is rotated by 90° around the x-axis). The prodomains are shown in light and dark grey, the furin processing site (R274-R-K-R-R278) is indicated in magenta, the activity-containing mature C-terminal domain is shown in green and cyan. The two residues differing between A-Gdf6b and B-Gdf6b are presented as spheres with their carbon atoms colored in green and cyan. D-E: Comparative magnifications of the structure of A-Gdf6b (D) and B-Gdf6b (E) around the p.Ser227Gly switch. Amino acid residues interacting with Ser227 or Gly227 are shown as sticks, and hydrogen bonds as yellow stippled lines. As shown, the side chain hydroxyl group of Ser227 (shown with carbon atoms colored in green and transparent spheres highlighting the van der Waals spheres of the atoms)

Figure 4 :

 4 Figure 4: Gene expression and functional evidence supporting a role of gdf6b as a potential master sex-determining gene in Pachón cavefish. A: Expression profiles of gdf6b in male and female trunks during early development from 10 to 90 days post-fertilization (dpf, males: solid line; females: dashed line) showing a significant over-expression in males compared to females starting from 16 dpf. Results are presented as log10 mean ± standard errors; black and white dots represent the individual values of relative expression in males and females, respectively. Statistical significance between males and females was tested with the Wilcoxon Rank Sum Test (Wilcoxon-Mann-Whitney Test) and only significant differences are shown (*** = P < 0.01; ** = P < 0.01; * = P < 0.05). B-C: Gonadal expression of gdf6b (in red) and the Sertoli and granulosa supporting cell marker gsdf (in white) in male (B) and female (C) differentiating gonads at 15 dpf showing that gdf6b is specifically expressed in male gonads with no strict colocalization with gsdf in male. See Figure S6B for additional stages of

  males and 11 females of Pachón cave Astyanax mexicanus. Fin samples (a narrow strip of the tail fin) were taken from the live specimens anesthetized by MS-222 (Merck KGaA, Darmstadt, Germany), while for direct preparations (chromosomes from kidneys and gonads), fishes were euthanized first using 2-phenoxyethanol (Sigma-Aldrich, St. Louis, MO, USA). All animal protocols were carried out in strict accordance with the French and European legislations (French decree 2013-118 and directive 2010-63-UE) applied for ethical use and care of laboratory animals used for scientific purposes. Sylvie Retaux and CNRS institutional authorizations for maintaining and handling A. mexicanus in experimental procedures were 91-116 and 91272105, respectively. For karyotypic analysis, all handling of fish individuals followed European standards in agreement with §17 of the Act No. 246/1992 to prevent fish suffering. The procedures involving fish were supervised by the Institutional Animal Care and Use Committee of the Institute of Animal Physiology and Genetics CAS, v.v.i., the supervisor´s permit number CZ 02361 certified and issued by the Ministry of Agriculture of the Czech Republic. Several sampling campaigns were carried out in the field in Mexico between 2013 and 2019, resulting in a collection of tail fin clips from wild-caught individuals from the Pachón cave.

  version 4.1.0 and are listed in Data S1E.

  sets P1-P2 and P3-P4 produce a single PCR fragment only in males (B+), primer set P5-P6 amplifies two bands in males (all B+) and only a single band in females (all B-), and primer set P9-P10 amplifies a single band in all individuals. For PCR reactions with the P1-P2 primer set, HiDi Taq DNA polymerase (myPOLS Biotec) was used for detecting a single nucleotide variation. PCRs were performed in a total volume of 12.5 µl containing 0.2 µM of each primer, a final concentration of 20 ng/µl gDNA, 200 µM of dNTPs mix, 1X of HiDi buffer (10X), and 2.5 U per reaction of HiDi DNA polymerase. Cycling conditions were as follows: 95 °C for 2 min, then 35 cycles of (95 °C for 15 seconds (secs) + 60 °C for 10 sec + 72 °C for 30 sec), and 72 °C for 5 min. For PCR reactions with the P3-P4 and P9-P10 primer sets, PCR reactions were performed in a final volume of 50 µl containing 0.5 µM of each primer, a final concentration of 20 ng/µl gDNA, 10 µM dNTPs mix, 1X of 10X AccuPrime™ buffer, and 0.5 µl per reaction of AccuPrime™ Taq DNA polymerase. Cycling conditions were as follows: 94 °C for 2 min, then 35 cycles of (94 °C for 30 sec + 58 °C for 30 sec + 68 °C for 1 min and 30 sec). For PCR reactions with the P5-P6 primer set, PCR reactions were performed in a final volume of 25 µl containing 0.5 µM of each primer, a final concentration of 20 ng/µl gDNA, 10 µM dNTPs mixture, 1X of Jumpstart™ buffer (10X), and 0.5 µl per reaction of Jumpstart™ Taq DNA polymerase. Cycling conditions were as follows: 95 °C for 2 min, then 35 cycles of (95 °C for 1 min + 60 °C for 30 sec + 72 °C for 1 min), and 72 °C for 5 min.

  510008), according to the manufacturer's protocols using 10 µl of blood as starting material, the Truseq DNA PCR-Free kit and Truseq DNA UD Indexes (Illumina, ref. 20015962, ref. 20020590), and the KAPA library Amplification kit (Roche, ref. KK2620). Hi-C library was sequenced in paired-end 2x150 bp mode on Novaseq6000 (Illumina), using half a lane of a SP flow cell (ref. 20027464).

  , with chromosomes being counterstained by 4′,6-diamidino-2-phenolindole (DAPI), 1.5 µg/ml in antifade (Cambio, Cambridge, United Kingdom).gdf6b probe synthesis for FISH mappinggDNA was extracted using NucleoSpin Kits for Tissue (Macherey-Nagel, Duren, Germany) as described above. A gdf6b fragment comprising the two exons, the intron and 2,260 bp of the proximal promoter (with a total size of 4,368 bp) was amplified by PCR in a total volume of 50 µl. The mixture contained 0.5 µM of each primer, a final concentration of 20 ng/µl gDNA, 1X of 10X AccuPrime™ PCR Buffer II, 1U/reaction of AccuPrime™ Taq DNA Polymerase, High Fidelity (Thermofisher), was adjusted to 50 µl with autoclaved and distilled water.Cycling conditions were as follows: 94 °C for 45 sec, then 35 cycles of (94 °C for 15 sec + 64 °C for 30 sec + 68 °C for 5 min and 30 sec), and 68 °C for 5 min. The resulting PCR product was cloned into TOPO TA cloning Kit XL (Thermofisher) and after sequence verification it was purified using NucleoSpin plasmid DNA purification kit (Machery-Nagel, Düren, Germany) according to the supplier's indications. This Pachón cave gdf6b cloned DNA fragment was labeled by nick translation with Cy3-dUTP using Cy3 NT Labeling Kit (Jena Bioscience, Jena, Germany). The optimal fragment size of the probe (approx. 200-500 bp) was achieved after 30 min of incubation at 15 °C.

Pearson' s

 s Chi-squared test with Yates' continuity correction. For gene expression, normality of data residuals, homogeneity of variances and homoscedasticity were verified before performing parametric or non-parametric tests. Consequently, statistical analyses were carried out only with non-parametric tests using RStudio (Open Source version) considering the level of significance at p<0.05. For comparisons between two groups, we used Wilcoxon signed rank test. All data are shown as Mean ± Standard Error of the Mean (SEM). patterns (panels A'-D'). B numbers were found to be variable among individuals (from 0 to 3 Bs) with all males having a single (panels A-A') or multiple Bs (panels B-B') in most of their metaphases, most females having no B (panels C-C'), and only a few females having rare B positive metaphases (panels D-D') Notice also the lack of C-bands on Bs suggesting that these Pachón cavefish male-predominant Bs are largely euchromatic. Scale bar = 10 µm. Male and female numbers referred to Data S1A.

Figure S2 :

 S2 Figure S2: Genomic repeated elements, gdf6b evolution and gdf6b sequence alignments. Related to Figure 2. A. Comparison of the GC content and repeated elements in Pachón cavefish B chromosome (B) versus Pachón cavefish A chromosomes. Percentage (in log scale) of GC content (GC level), short interspersed repeated sequences (SINEs), long interspersed nuclear elements (LINEs), long terminal repeats (LTRs), DNA repeat elements (DNAs),

Figure S3 :

 S3 Figure S3: PCR genotyping of the Pachón cavefish male-predominant B chromosome and the gdf6b knockout mutant individuals. Related to Figure 4 and STAR Methods. A.Absence / presence of the Pachón cavefish B chromosome was detected with different primer pairs based on differences between the A-gdf6b and B-gdf6b loci. Three sets of primers (P) were designed with two primer sets designed to amplify specifically the B-gdf6b copy based either on a single base variation between the A/B gdf6b CDS at position 679 bp (P1-P2 primer set), or based on primers located on both sides of the A/B breakpoints downstream of the B-gdf6b gene (P3-P4 primer set). The third set of primers (P5-P6 primer set) was designed on gaps/indels variations between A-gdf6b and B-gdf6b in the proximal promoter of gdf6b genes.Primer sets P1-P2 and P3-P4 produce a single PCR fragment only in B+ individuals, and primer set P5-P6 amplifies two bands in B+ individuals and only a single band in B-individuals. Gene knockout (KO) was performed by genome editing using the CRISPR-cas9 method with two guide RNAs target sites (#TS1 and #TS2) designed in order to target gdf6b exon 2 (E2). KO individuals were genotyped using primers (P7-P8 primer set) flanking the 2 target sites that induced a 470 bp deletion in gdf6b mutant individuals (Mut) compared to the wild type (WT) gdf6b sequence. B. Increasing PCR cycle numbers allows the detection of a faint PCR fragment specific to the B chromosome in Pachón cavefish females. PCRs were carried out in one male and two female Pachón cavefish with the protocol described for the P3-P4 primer set and with increasing PCR cycle numbers. The B chromosome was detected using primers specifically designed to amplify the B-gdf6b loci (primer set P3-P4) and a control was incorporated with primers specifically designed to amplify the A-gdf6b locus (primer set P9-P10). NC: PCR negative control. C. Numbers of gdf6b knockout generated by CRISPR-Cas9 method including the number of fish injected, and the number and percentage of sex-reversed males obtained (in bold red type). Out of the 200 micro-injected eggs (at 1 cell stage), 122 adult fishes were obtained including 80 genetic males and 42 genetic females. Among the 80 genetic males, we found 62 phenotypic males (77.5%), and 18 phenotypic females (22.5%) displaying a 470 bp deletion on the exon 2 of their A-gdf6b and/or B-gdf6b gene.

  technology (ACD Biotechne™) on 7 µm cross-sections of 15, 21, 30 and 60 dpf Pachon cavefish fixed in paraformaldehyde 4% overnight at 4 °C and embedded in paraffin after serial dehydration in increasing methanol solutions. Specific probes for Pachon cavefish gdf6b and gsdf were synthesized by ACD Biotechne™. Sections were collected on Super frost+ slides, heated at 60 °C for 1 h and dewaxed 2 x 5 min in xylene followed by 2x 2 min in Ethanol 100% at RT. Fluorescent ISH was carried out with the Multiplex Fluorescent Reagent Kit v2 (ACD

	Biotechne™, ref: 323100) according to manufacturer's protocol. Following hybridization,
	nuclei were labelled with DAPI (4′,6-diamidino-2-phenylindole) staining and slides were
	mounted with ProLong™ Gold Antifade Mountant (Invitrogen™) and observed with a Leica
	TCS SP8 laser scanning confocal microscope.

ACKNOWLEDGEMENTS

We thank Victor Simon, Stéphane Père, Krystel Saroul, Pierre-Lo Sudan and Amélie Patinote for taking care and handling cavefish, and Manon Thomas and the LPGP TEFOR infrastructure platform for acquisition of the RNAScope confocal images. This project was supported by funds from the "Agence Nationale de la Recherche" (ANR/DFG, PhyloSex project, 2014-2016) to Y.G and M.S. S.R was supported by grants from an Equipe FRM (Fondation pour la Recherche Médicale, DEQ20150331745) and MITI CNRS (Mission pour les Initiatives Transverses et Interdisciplinaires). J.H.P was supported by a NIH grant (R35 GM139635). Sequencing was supported by France Génomique as part of an "Investissement d'avenir" program managed by ANR (contract ANR-10-INBS-09) and by the GET-PACBIO program (« Programme opérationnel FEDER-FSE MIDI-PYRENEES ET GARONNE 2014-2020 »). The CytoEvol platform at ISEM was supported by the Labex CeMEB. J.B, T.P and A.S were supported by RVO: 67985904 of IAPG CAS, Liběchov. T.P was supported by the projects of the Czech Ministry of Education (SVV 260571/2021). S.S was supported by the Russian Foundation for Basic Research (RFBR) (18-34-00638). B.I PhD fellowship was supported by the Doctoral School of Ecology, Geosciences, Agronomy, Nutrition of the University of Rennes 1 and INRAE. We are grateful to the genotoul bioinformatics platform Toulouse Occitanie (Bioinfo Genotoul, https://doi.org/10.15454/1.5572369328961167E12) for providing help, computing and storage resources. Funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

AUTHOR CONTRIBUTIONS

Conceptualization, S.R., A.H. and Y.G.; Methodology, B.I., S.R., A.H., A.S., and Y.G.; Formal Analysis, T.D.M., K.D., R.F., Q.P., C.C., M.Z., and C.K.; Investigation, B.I., S.B., E.J., J.T-P., C.L-R., A.C., L.G., C.K., C.D., H.P., L.J., T.P., A.A-R., T.L., S.S., J.B., A.S., J.P., F.V.;

Writing -Original Draft, B.I., S.R., and Y.G.; Writing -Review & Editing, A.S., J.H.P., T.D.M., and M.S.; Visualization, B.I., K.D., E.J., R.F., Q.P., C.C., A.S., T.D.M., and Y.G.; Funding Acquisition, S.R., M.S., J.H.P., and Y.G.; Supervision, S.R., A.H., and Y.G.

DECLARATION OF INTERESTS

The authors declare no competing interests FIGURE the final pool. Pool-sequencing libraries were prepared using the Illumina TruSeq Nano DNA HT Library Prep Kit (Illumina, San Diego, CA) according to the manufacturer's protocol. After the fragmentation of each gDNA pool (200 ng/pool) by sonication using an M220 Focusedultrasonicator (COVARIS), the size selection was performed using SPB beads retaining fragments of 550 bp. Following the 3′ ends of blunt fragments mono-adenylation and the ligation to specific paired-end adaptors, the amplification of the construction was performed using Illumina-specific primers. Library quality was verified with a Fragment Analyzer (Advanced Analytical Technologies) and then quantified by qPCR using the Kapa Library Quantification Kit (Roche Diagnostics Corp, Indianapolis, IN). The enriched male and female pool libraries were then sequenced using a paired-end multiplexed sequencing mode on a NovaSeq S4 lane (Illumina, San Diego, CA), combining the two pools on the same lane and producing 2×150 bp with Illumina NovaSeq Reagent Kits according to the manufacturer's instructions. Sequencing produced 288 million paired reads and 267 million paired reads for the male and female pool libraries, respectively.

Pool-sequencing analysis

Characterization of genomic regions enriched for sex-biased signals between Pachón cave males and females, consisting of coverage and Single Nucleotide Variations (SNVs) was performed as described previously [START_REF] Feron | Characterization of a Y-specific duplication/insertion of the anti-Mullerian hormone type II receptor gene based on a chromosome-scale genome assembly of yellow perch, Perca flavescens[END_REF][START_REF] Pan | Identification of the master sex determining gene in Northern pike (Esox lucius) reveals restricted sex chromosome differentiation[END_REF] . Pachón cave A. mexicanus paired-end reads from male and female pool-seq pools were mapped onto our own Pachón cave genome assembly using BWA mem version 0.7.17 [START_REF] Li | Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM[END_REF] . The resulting BAM files were sorted and the duplicate reads due to PCR amplification during library preparation were removed using Picard tools version 2.18.2 (http://broadinstitute.github.io/picard) with default parameters. Then, for each pool and each genomic position, a file containing the nucleotide composition was generated using samtools mpileup 77 version 1.8, and popoolation2 78 mpileup2sync version 1201. This file was then analyzed with custom software (PSASS version 2.0.0; doi: 10.5281/zenodo.2615936) to compute: (a) the position and density of sex-specific SNVs, defined as SNVs heterozygous in one sex but homozygous in the other, and (b) the average read depths for male and female pools along the genome to look for regions present in one sex but absent in the other (i.e., sex-specific insertions). All PSASS analyses were run with default parameters except for the range of frequency for a sex-linked SNV in the homogametic sex, --range-hom, that was set to 0.01 instead of 0.05, and the size of the sliding window, --window-size, that was set at 50,000 instead of 100,000.

LEGEND OF SUPPLEMENTARY DATAFILE Data S1: Supplementary information on metaphase numbers (A), genome assembly metrics and annotation (B, C, D), primer names (E) and sex linkage of B-gdf6b (F) in

Pachón cavefish. Related to STAR Methods and Figures 1 and2