# Protein expression reveals a molecular sexual identity of avian primordial germ cells at pre-gonadal stages

Laura Soler-Vasco, Sabine Alves, Aurélien Brionne, Aurore Jacques, Vanessa Guérin, Maeva Cherif-Feildel, Lucie Combes-Soia, Sophie Fouchécourt, Aurore Thélie, Elisabeth Blesbois, et al.

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

► To cite this version:

Laura Soler-Vasco, Sabine Alves, Aurélien Brionne, Aurore Jacques, Vanessa Guérin, et al. Protein expression reveals a molecular sexual identity of avian primordial germ cells at pre-gonadal stages. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-98454-2, hal-03367022

HAL Id: hal-03367022

https://hal.science/hal-03367022

Submitted on 13 Oct 2021

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Protein expression reveals a molecular sexual identity of avian primordial germ cells at pre-gonadal stages

Laura Soler<fnref n="1" />,²,⁶, Sabine Alves<fnref n="1" />,⁶, Aurélien Brionne³, Aurore Jacques<fnref n="1" />, Vanessa Guérin<fnref n="1" />, Maeva Cherif-Feildel<fnref n="1" />, Lucie Combes-Soia<fnref n="1" />,⁴, Sophie Fouchécourt<fnref n="1" />, Aurore Thélie<fnref n="1" />, Elisabeth Blesbois<fnref n="1" />, Michael J. McGrew⁵, Valérie Labas<fnref n="1" />,⁴ & Marina S. Govoroun<fnref n="1" />✉

In poultry, in vitro propagated primordial germ cells (PGCs) represent an important tool for the cryopreservation of avian genetic resources. However, several studies have highlighted sexual differences exhibited by PGCs during in vitro propagation, which may compromise their reproductive capacities. To understand this phenomenon, we compared the proteome of pregonadal migratory male (ZZ) and female (ZW) chicken PGCs propagated in vitro by quantitative proteomic analysis using a GeLC-MS/MS strategy. Many proteins were found to be differentially abundant in chicken male and female PGCs indicating their early sexual identity. Many of the proteins more highly expressed in male PGCs were encoded by genes localised to the Z sex chromosome. This suggests that the known lack of dosage compensation of the transcription of Z-linked genes between sexes persists at the protein level in PGCs, and that this may be a key factor of their autonomous sex differentiation. We also found that globally, protein differences do not closely correlate with transcript differences indicating a selective translational mechanism in PGCs. Male and female PGC expressed protein sets were associated with differential biological processes and contained proteins known to be biologically relevant for male and female germ cell development, respectively. We also discovered that female PGCs have a higher capacity to uptake proteins from the cell culture medium than male PGCs. This study presents the first evidence of an early predetermined sex specific cell fate of chicken PGCs and their sexual molecular specificities which will enable the development of more precise sex-specific in vitro culture conditions for the preservation of avian genetic resources.

The protection of endangered avian breeds requires the preservation of the entire genotype of the animal. Bird oocytes and zygotes are very lipid rich, which hinders their cryopreservation. For this reason, the conservation of avian genetic resources is essentially based on semen cryopreservation. Semen-based breed recovery, however, is not optimal as it requires many successive generations of insemination to restore the breed genotype and does not conserve the female specific W sex chromosome<fnref n="1" />. Hence, there is a need to develop new reproductive technologies that allow the collection, propagation and preservation of avian reproductive germ cells.

The precursors of the germ line are called primordial germ cells (PGCs). PGCs are embryonic unipotent stem cells fated to develop into gametes in the appropriate cellular and molecular environment. Aves is one of the few animal classes from which PGCs can be isolated at early embryonic stages due to their precocious formation in avian species and the particularities of their migration through the avian embryo. Avian PGCs can be isolated at pre-migratory stages from the embryonic germinal crescent, during their migration in the gonads by blood sampling, or from the gonads before the onset of their sex differentiation²,³.

Isolated PGCs from the chicken, *Gallus gallus*, can be propagated indefinitely in vitro, cryopreserved and/or transplanted into host embryos, where they will develop into functional gametes⁴⁻⁶. In poultry, the development

<footnote n="1">¹INRAE, CNRS, Université de Tours, IFCE, PRC, 37380 Nouzilly, France. ²Present address: Toxalim (Research Centre in Food Toxicology), INRAE, Université de Toulouse, ENVT, INP-Purpan, UPS, 31027 Toulouse, France. ³INRAE, Université de Tours, BOA, 37380 Nouzilly, France. ⁴INRAE, CHU de Tours, Université de Tours, PRC, PIXANIM, 37380 Nouzilly, France. ⁵The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush Campus, Midlothian EH25 9RG, UK. ⁶These authors contributed equally: Laura Soler and Sabine Alves. ✉email: marina.govoroun@inrae.fr</footnote>
of these methods has advanced the use of PGCs for the conservation of endangered indigenous breeds and of valuable commercial genetic resources as well as offering an efficient vehicle for genome editing and transgenesis⁶⁻⁸.

In the animal kingdom, PGCs show some degree of sexual plasticity, which is inversely correlated with their level of complexity. In molluscs, fish, amphibians and other animals with poorly differentiated sex chromosomes, male and female PGCs are sexually bipotential and may develop into spermatozoa or ova depending on hormonal and cellular contexts regulated by the environment⁹. In animals with highly differentiated X and Y sex chromosomes, such as mammals, the differentiation of male and female PGCs into functional gametes is highly restricted by chromosomal factors¹⁰. However, several experiments carried out in rodents have shown that PGCs can have some sexual plasticity. In some natural genetic contexts or in experimentally created sexual reversals, XY germ cells may develop into functional oocytes although with a considerable decrease in fertility¹¹. In contrast, XX germ cells never develop into spermatozoa, the second X chromosome being incompatible with spermatogenesis¹⁰. In addition, some Y-linked genes are required for the resumption of second meiotic division in spermatocytes¹².

Birds have a ZZ/ZW system where the female is heterogametic (ZW), and male is homogametic (ZZ). The sex chromosomes Z and W are well differentiated, however, the existing knowledge on the sexual plasticity of avian PGCs is still rudimental and contradictory. Studies on chicken germline intersex chimeras obtained after the direct transplantation of donor PGCs to recipient embryos of the opposite sex suggest that chicken germ lineage development is under the strict control of genetic and epigenetic factors, but might have some degree of sexual plasticity. After the direct transplantation following PGC isolation of a ZZ male donor/ZW female recipient combination, progeny from donor-derived cells is very rare event. No germline transmission is generally found in the case of the reverse ZW female donor/ZZ male recipient combination¹³. Similarly, no donor-derived progeny were observed after intersex transplantation of ZZ or ZW PGCs, when they have been propagated in vitro before transplantation⁶,¹⁴⁻¹⁶. Moreover, *in ovo* treatment of female embryos with aromatase inhibitor induces a female-to-male sex reversal, but these animals produce non-functional spermatozoa¹⁷,¹⁸. However, ZZ donor PGCs propagated in vitro were observed to differentiate into functional oocytes and produce progeny after transplantation to partially sterilised ZW host embryo¹⁹. These data suggest that chicken PGC developmental plasticity may depend on several factors including hormonal environment, the host biological context and in vitro culture conditions, which can reveal sex-specific features of these cells. Indeed, survival and propagation of chicken PGCs in vitro is sex-dependent³,¹⁴ and their reproductive capacity may be compromised in a sex specific manner during in vitro propagation⁴,⁵. However, the molecular nature of this phenomenon is largely unknown.

The objective of the present study is to provide for the first time a description of the intrinsic differences of the molecular phenotypes of male (ZZ) and female (ZW) chicken PGCs isolated and propagated in vitro, using media conditions that support the proliferation of PGCs from both sexes⁵. We investigated the protein differences of ZZ and ZW chicken PGCs by a quantitative bottom-up proteomic approach using a GeLC-MS/MS strategy, followed by a functional analysis of the data. This allowed us to identify and validate differential biologically relevant pathways and cellular functions, pointing to an early sexual molecular identity of ZZ and ZW PGCs.

## Methods

**Ethical statement.** The housing and the breeding procedures to obtain eggs were approved by local official ethics committee (Comité d'Ethique en Expérimentation Animale Val de Loire CEEA - n°19), by the French Ministry of Agriculture and French Veterinary Services regulations under agreement number C37-175-1. The experiments were performed in accordance to European welfare regulations (European directive 2010/63/EU).

**Animals.** Animals were housed at the INRA experimental facility PEAT: https://doi.org/10.15454/1.55723 26250887292E12 (Nouzilly, France). Twenty weeks old male and female chickens of local breed Noire du Berry were obtained from the Club Français de la Poule Noire du Berry. This local breed was used as it is part of a Genetics Resource conservation program developed for French rare endangered breeds. Males were housed in individual battery cages and females in groups of 5 under 14L/10D photoperiod and fed with a standard diet of 12.5 MegaJoules/day supplemented with calcium for females. For fertile egg production, 20-weeks old hens were fertilized by artificial insemination using pooled sperm from 5 males that was collected twice a week by massage²⁰. Eggs were incubated at 37.7° C and 45-51% humidity for 56 h with turning over an angle of 45° either side of the short axis of the egg.

**Culture derivation and preparation of samples for further analyses.** After 56 h of incubation, a small window was made in the egg shell of the pointed end of the egg and 1.0-2.0 µls of blood containing circulating PGCs were aspirated from dorsal aorta using fine glass capillary and placed immediately in culture medium. The culture media and culture conditions were set up as described in⁵ with slight modifications: instead of using Activin A, or Activin A and BMP4, the culture medium was supplemented with 25 ng/ml of BMP4 only. Five male and five female PGCs cultures deriving from individual embryos were in vitro expanded for 1.5-2 months until each culture population reached 30 million cells. For protein extraction, 2.0 × 10⁶ cells were centrifuged at 1000 g for 10 min, washed twice with 50 volumes of Mg²⁺ and Ca²⁺ free PBS buffer and frozen immediately in liquid nitrogen. For RNA extraction 5.0 × 10⁵ cells were centrifuged at 1000 g for 10 min, then resuspended in lysis buffer of NucleoSpin RNA XS kit (MACHEREY-NAGEL, France) and frozen immediately in liquid nitrogen. All samples were stored at -80 °C until proceeding with protein or RNA extraction.

**Sexing of embryos and cells.** Sexing of embryos and PGCs was performed by PCR as described previously²¹⁻²³. For DNA extraction from PGCs, 60,000 cultured PGCs were incubated in 20 µl of proteinase
K 0.25 mg/ml water solution at 55 °C for 1 h. After inactivation of proteinase K at 99° C for 10 min, DNA were stored at 4 °C for further analysis.

**Protein extraction and SDS-PAGE.** Proteins were extracted from samples containing 2.0 × 10⁶ PGCs in 250 µl of 6 M Urea, 4% SDS, 50 mM pH 8.8 Tris-HCl containing protease inhibitors (Roche diagnostic, Paris, France) by sonication on ice. Extracted proteins were allowed to solubilize for 20 min at room temperature under shaking. After centrifugation at 10,000 g for 10 min, protein concentration was quantified from protein extracts (Thermo Scientific Pierce BCA Protein Assay Kit). Male and female PGCs protein extract pools were prepared by combining 10 µg of total protein from each sample. Fifty µg of each pool were reduced using DTT 20 mg/ml and denatured at 95 °C for 5 min before being loaded in a 4–12% SDS-PAGE gel. The gel was run at 30 mA until full migration of the bromophenol blue front.

**In-gel protein digestion and NanoLC-MS/MS for PGC protein identification.** After separation of proteins by SDS-PAGE and Coomassie Blue staining, each lane was cut into 10 slices. After washing of gel slices, proteins were reduced, alkylated and in-gel digested by trypsin, as described in Supplementary Methods. For each protein band, the resultant peptide mixture was analysed by nanoflow liquid chromatography tandem mass spectrometry (nanoLC-MS/MS) in triplicate. MS/MS ion searches were performed using Mascot search engine v 2.3 (Matrix Science, London, UK) and confronted against the “chordata” section of a locally maintained copy of nr NCBI (download Jan. 2018, 9 822 753 sequences). After identification, proteins were validated using Scaffold software (v 4.8, Proteome Software, Portland, USA).

**Label-free protein quantification.** For comparative analysis, we employed Scaffold Q+ software (version 4.8.4, Proteome Software, Portland, USA) to apply two independent quantitative methods: (1) the Spectral Counting (SC) using the “Weighed Spectra” option; (2) the Average Precursor Intensity (API). Thus, values of Normalized Weighed Spectra (NWS) and Normalized Average Precursor Intensity (NAPI) were tabulated using experiment wide protein clustering. Statistical analyses were performed using T-tests, where *p* < 0.05 was considered significant. Limits of an average normalized weighted spectra (NWS) ≥ 3 and male/female fold-change (FC) ≥ 2 and < 0.5 were used to increase validity.

**Bioinformatic analysis.** Biological interpretations were carried out using Gene Ontology (GO) public database with the use of Biological process (BP) and Molecular Function (MF) categories using VISEAGO R package²⁴. Associated Gene terms, from the identified *Gallus gallus* proteins, were retrieved from EnterGene²⁵ and well-annotated species orthologues. Enrichment tests were performed for each studied comparison. All enriched GO terms (*p* < 0.01) were grouped into functional clusters using hierarchical clustering based on Wang’s semantic similarity between GO terms respecting GO graph topology and Ward’s criterion.

Chromosome enrichment tests were performed for each cluster independently using Fisher’s exact test (*p*-value < 0.01) using R version 3.6.1.

For proteins for which automatic analysis assigned their corresponding genes to sex-chromosomes Z or W, blast P and TblastN²⁶ were performed individually in Gallus gallus Ensembl genome browser²⁷ and NCBI²⁸ in order to search for their counterparts on the second sex chromosome.

**Western blot.** After PGC protein extraction, proteins were separated in 10% SDS-PAGE and blotted onto nitrocellulose membranes. To evaluate hydroxysteroid dehydrogenase like 2 (HSDL2), ovalbumin and insulin-like growth factor 2 binding protein 1 (IGF2BP1) abundance rabbit anti-human HSDL2 antibody HPA050453(1:500), rabbit anti-chicken egg albumin antibody C6534 (1:500) (Sigma Aldrich, France) mouse anti-human IGF2BP1antibody PCRP-IGF2BP1-2D4 (DSHB, Iowa, USA) (1:50) and secondary goat anti-rabbit or goat anti-mouse IRDy e800CW secondary antibody (LI-COR Biosciences—GmbH) were used for HSDL2 and Ovalbumin, and IGF2BP1 respectively. Total protein staining of the membrane with SyproRuby blot stain (ThermoFischer) or Revert 700 Total Protein Stain (Li-Cor LI-COR Biosciences—GmbH) prior to blocking served as a loading control. For each lane, band intensity values were normalized onto the respective overall protein staining. Detailed procedure is described in Supplementary Mat&Met. Statistical analysis was performed by Fisher’s F-test to compare two variances, followed by Student’s T-test to reveal significant differences. Differences were considered significant when *p* < 0.05.

**Reverse transcription quantitative real time polymerase chain reaction (RT-qPCR).** Total RNA was isolated from individual PGCs cultures using NucleoSpin RNA XS kit (MCHHEREY-NAGEL, France) and was reverse transcribed (RT) using Thermo Scientific Maxima First Strand cDNA Synthesis Kit for RT-qPCR (ThermoFisher Scientific, France) according manufacture instructions. RT-qPCR was performed using SsoAdvanced Universal SYBR Green Supermix (BIORAD, France) according manufactory instructions. Three technical replicates were performed for each reaction. At least two independent qPCR experiments were carried out for each sample and for each primer pair. The sequences of primers used for cDNAs amplifications are presented in Supplementary Table S1. Efficacy of each primer pairs was between 1.8 and 2.2. The relative amount of transcript expression was calculated by LightCycler 480 software release 1.5.0 using standard curve composed of serial dilutions of the cDNA obtained by RT of total RNA extracted from the pool of PGCs containing equal amounts of five ZZ and five ZW PGCs cultures used for proteomic analysis. RT-qPCR products were verified by sequencing at Genewiz (Germany). For cultured PGCs characterisation the relative transcript expression level of target genes was normalized by the mean of expression ratios of two housekeeping genes: eukaryotic translation
elongation factor 1 alpha 1 (EEF1A) and ribosomal protein L15 (RPL15). For evaluation of the expression of transcripts corresponding to differentially abundant proteins four housekeeping genes were used: EEF1A, actin, beta (ACTB), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), RPL15. Statistical analysis was performed by Fisher's F-test to compare two variances, followed by Student's T-test to reveal significant differences. Differences were considered significant when *p* < 0.05.

**Immunohistochemistry.** After fixation and permeabilization on glass slides, cells were incubated with rabbit anti chicken DDX4 antibodies (1:500) kindly provided by B. Pain²⁹ or with mouse anti-mouse SSEA1 monoclonal antibody (DSHB, Iowa, USA) (1:40). After washing cells were incubated with corresponding cross-Adsorbed Alexa Fluor secondary antibodies. For negative control Rabbit IgG or Mouse IgG (Invitrogen) were used as primary antibodies for DDX4 or SSEA1 analysis, respectively. Detailed procedure is described in supplementary Mat&Met.

**Gonad colonisation test.** The PGCs were labelled with green fluorescent dye PKH67 (SIGMA-Aldrich, France) according manufacture instructions. 1.0 µl of culture medium containing 2000 labelled PGCs was injected into the embryonic dorsal aorta through the small window in the eggshell of the egg incubated for 56 h. The shell window was sealed with Parafilm M (PRAXISDIENST, France) and the eggs were incubated for 90 additional hours. The embryos were removed from the egg, euthanized by decapitation and the gonads-mesonephros complexes were isolated. Green fluorescence was observed using a Zeiss Axiovert 200 microscope (Zeiss, France) microscope and the images were captured by Axiocam type MRm camera and analysed using AxioVision software (Zeiss, France).

## Results and discussion

**Proteomic analysis reveals sexual differences in chicken PGCs.** Here, for the first-time, we conduct a comparative proteomic analysis of functional male and female PGCs propagated in vitro. Prior to the proteomics analysis, the basic molecular phenotype of the cultured male and female PGCs was confirmed using RT-qPCR analysis of selected key genes and immunohistochemistry (Supplementary Fig. S1a, b). As expected, male and female PGCs showed no difference in the expression of genes that are specific for germ cells (Deleted In Azoospermia Like; *DAZL*) and pluripotency factors (SRY-Box Transcription Factor 2; *SOX2*, Nanog Homeobox; *NANOG*) (Supplementary Fig. S1a).

A male molecular phenotype was confirmed in male PGCs by the expression level of the germ cell specific gene, *DDX4*, located in the Z-chromosome. Unlike mammals, birds do not fully compensate for lower gene expression of the Z chromosome genes in females to equalise expression in both sexes. This results in a significantly lower global expression in females (ZW) of Z-linked genes in all somatic cells³⁰. According to the lack of compensation of double dosage of Z-linked chromosomal genes in birds, we expected *DDX4* to be two times more highly expressed in male (ZZ) PGCs than in female (ZW) PGCs. Our gene expression results are in agreement with the double dosage of the Z-linked gene *DEAD-Box Helicase 4 (DDX4)* in males and with the absence of general mechanism of dosage compensation of Z-linked genes in the chicken. Immunohistochemistry analysis confirmed the presence of the germ cell-specific proteins DDX4 and stage-specific embryonic antigen-1 (SSEA-1) in PGCs from both sexes (Supplementary Fig. S1b). Additionally, we confirmed that in vitro propagated PGCs conserved their basic ability to colonize the embryonic gonads when injected in the dorsal aorta of recipient chicken embryos (Supplementary Fig. S1c).

GeLC-MS/MS analysis identified 1332 proteins with unique accession numbers distributed in 777 protein clusters (Supplementary Table S2a, b). Among these proteins, 299 were differentially abundant (*p*-value < 0.05 and male/female (M/F) FC limit 2 < FC or FC < 0.5) in male and female PGCs as determined by at least one of two used quantification methods: WS (columns G to N) and/or API (columns O to V). Among these 299 differentially abundant proteins with unique accession numbers, 154 were preferentially abundant in male PGCs and 155 were more abundant in female PGCs (Supplementary Table S2c, d). Among 154 proteins preferentially abundant in male PGCs, 32 proteins corresponding to 28 gene IDs were quantified by both WS and API, 90 proteins corresponding to 72 gene IDs were quantified by WS only and 32 proteins corresponding to 26 gene IDs by API only. Among 155 proteins preferentially abundant in female PGCs, 41 proteins corresponding to 30 gene IDs were quantified by both WS and API, 97 proteins corresponding to 78 gene IDs were quantified by WS only and 17 proteins corresponding to 16 gene IDs by API only. Both quantification methods were, thus, complementary. Out of the top 20 differentially abundant proteins (Table 1), 14 were identified only in female PGCs and 11 only in male PGCs.

Differential protein abundance found by quantitative proteomic analysis was confirmed by western blot analysis for two proteins, HSDL2 and IGF2BP1, carried out on individual male and female PGCs cultures (Fig. 1, Supplementary Fig. S2), and for ovalbumin, on the pool of male and female individual PGCs cultures (supplementary Fig. S3).

**The male PGC proteome is enriched in proteins encoded in the Z chromosome.** We investigated the chromosomal position of the genes encoding differentially abundant proteins and we observed a weak enrichment of genes encoding female overabundant proteins on chromosomes 7 and 25, whereas a strong enrichment of genes encoding male overabundant proteins on chromosome Z (Supplementary Table S3a). Indeed, 15.5% (17 out of 111 unique chicken ID found for proteins more abundant in male PGCs) of more abundant proteins in male PGCs were encoded by genes located on the Z chromosome, whereas only 0.9% (1 out of 108 unique chicken IDs found for proteins more abundant in female PGCs) of female more abundant proteins was linked to the Z chromosome (Table 2, Supplementary Table S3b). When considering all of the identified proteins, only
<table><thead><tr><th rowspan="2">NCBI<br>Protein accession<br>number</th><th rowspan="2">Gene ID</th><th rowspan="2">Gene Name</th><th rowspan="2">Chr</th><th rowspan="2">Description</th><th colspan="3">Normalized weighed spectra</th><th colspan="3">Normalized precursor intensity</th></tr><tr><th>Male/Female</th><th>Average</th><th>Male/Female</th><th>Average</th><th>Male</th><th>Female</th></tr><tr><th></th><th></th><th></th><th></th><th></th><th>T-Test</th><th>Male</th><th>Female</th><th>T-Test</th><th>Male</th><th>Female</th></tr></thead><tbody><tr><td>NP_001006398.1</td><td>421348</td><td>EPRS</td><td>3</td><td>Multifunctional aminoacyl-tRNA synthetase that catalyzes the aminoacylation of glutamic acid and proline tRNA species</td><td>0.00028</td><td>10.58</td><td>0.00</td><td>0.006</td><td>4.24E+06</td><td>8.55E-01</td></tr><tr><td>NP_989854.1</td><td>395194</td><td>TLN1</td><td>Z</td><td>Cytoskeletal protein concentrated in areas of cell-substratum and cell-cell contacts. Plays a significant role in the assembly of actin filaments and in spreading and migration of various cell types</td><td>0.0094</td><td>9.06</td><td>0.00</td><td>0.0083</td><td>3.34E+06</td><td>0.00</td></tr><tr><td>XP_021238991.1</td><td>110391426</td><td>HUWE1</td><td>unplaced</td><td>Functions as an E3 ubiquitin ligase required for the ubiquitination and subsequent degradation of the anti-apoptotic protein, the p53 tumor suppressor, core histones, and DNA polymerase beta</td><td>0.0059</td><td>10.57</td><td>0.00</td><td>0.0014</td><td>5.07E+06</td><td>5.52E-01</td></tr><tr><td>XP_005533424.1</td><td>107050871</td><td>LOC107050871</td><td>33</td><td>DNA-dependent RNA polymerase that catalyzes the transcription of DNA into RNA using the four ribonucleoside triphosphates as substrates</td><td>0.015</td><td>13.54</td><td>0.00</td><td>0.00015</td><td>5.39E+06</td><td>5.99E-01</td></tr><tr><td>XP_418787.2</td><td>420688</td><td>DNAJC13</td><td>2</td><td>Member of the DnaJ protein family whose members act as co-chaperones of a partner heat-shock protein by binding to the latter and stimulating ATP hydrolysis</td><td>0.00058</td><td>8.42</td><td>0.00</td><td>0.00045</td><td>4.58E+06</td><td>8.55E-01</td></tr><tr><td>XP_004940443.1</td><td>42,811</td><td>MDN1</td><td>3</td><td>Nuclear chaperone required for maturation and nuclear export of pre-60S ribosome subunits</td><td>0.017</td><td>4.39</td><td>0.00</td><td>0.00039</td><td>5.02E+06</td><td>8.55E-01</td></tr><tr><td>XP_004949443.1</td><td>427114</td><td>TTC37</td><td>Z</td><td>Protein rich in tetratricopeptide repeats, which mediate protein-protein interactions and chaperone activity</td><td>0.00096</td><td>2.91</td><td>0.00</td><td>0.00015</td><td>2.53E+06</td><td>5.58E-01</td></tr><tr><td>NP_001026553.1</td><td>426429</td><td>OGDH</td><td>22</td><td>Catalyzes the overall conversion of 2-oxoglutarate (alpha-ketoglutarate) to succinyl-CoA and CO(2) during the Krebs cycle</td><td>0.016</td><td>4.38</td><td>0.00</td><td>0.0046</td><td>1.42E+07</td><td>5.99E-01</td></tr><tr><td>NP_001026640.1</td><td>427694</td><td>SBNO1</td><td>15</td><td>Chromatin DNA and histone binding protein implicated in the regulation of transcription</td><td>0.039</td><td>4.02</td><td>0.00</td><td>0.00054</td><td>4.95E+06</td><td>3.03E-01</td></tr><tr><td>NP_001026029.1</td><td>419257</td><td>TPD52L2</td><td>20</td><td>Member of the tumor protein D52-like family that regulated cell proliferation and metabolic processes</td><td>0.00012</td><td>5.85</td><td>0.00</td><td>0.00029</td><td>6.16E+06</td><td>0.00E+00</td></tr></tbody><tfoot><tr><td colspan="12">Continued</td></tr></tfoot></table>
<table><thead><tr><th rowspan="2">NCBI<br>Protein accession<br>number</th><th rowspan="2">Gene ID</th><th rowspan="2">Gene Name</th><th rowspan="2">Chr</th><th rowspan="2">Description</th><th colspan="3">Normalized weighed spectra</th><th colspan="3">Normalized precursor intensity</th></tr><tr><th>Male/Female</th><th>Average</th><th>Male/Female</th><th>Average</th><th>Male</th><th>Female</th></tr><tr><th></th><th></th><th></th><th></th><th></th><th>T-Test</th><th>Male</th><th>Female</th><th>T-Test</th><th>Male</th><th>Female</th></tr></thead><tbody><tr><td>NP_001004402.1</td><td>422842</td><td>WDR1</td><td>4</td><td>Actin binding protein that induce disassembly of actin filaments in conjunction with ADF/cofilin family proteins. Involved in cytokinesis</td><td>0.00096</td><td>2.91</td><td>0.00</td><td>0.0025</td><td>9.72E+06</td><td>2.96E-01</td></tr><tr><td>BAA24137.1</td><td>107050620</td><td>LOC107050620</td><td>36</td><td>Chicken high molecular mass nuclear antigen, involved in mRNA translation</td><td>0.0005</td><td>29.66</td><td>0.61</td><td>0.0041</td><td>8.15E+06</td><td>9.26E+05</td></tr><tr><td>NP_990427.2</td><td>395984</td><td>HDLBP</td><td>9</td><td>High density lipoprotein binding protein that plays a role in cell sterol metabolism protecting cells from over-accumulation of cholesterol</td><td>0.004</td><td>9.27</td><td>0.51</td><td>0.0074</td><td>4.93E+06</td><td>8.55E-01</td></tr><tr><td>XP_424271.3</td><td>426644</td><td>ERMP1</td><td>Z</td><td>Zinc-binding protease belonging to the peptidase M28 family. Its expression is required in the ovary for the organization of somatic cells and oocytes into discrete follicular structures</td><td>0.0001</td><td>82.91</td><td>5.05</td><td>0.018</td><td>5.28E+07</td><td>7.65E+06</td></tr><tr><td>XP_021273023.1</td><td>110408544</td><td>UBR4</td><td>20</td><td>E3 ubiquitin-protein ligase that is a cytoskeletal component in the cytoplasm and part of the chromatin scaffold in the nucleus</td><td>0.00026</td><td>11.34</td><td>0.84</td><td>0.00044</td><td>5.09E+06</td><td>8.55E-01</td></tr><tr><td>NP_990842.1</td><td>396517</td><td>SLC2A14</td><td>1</td><td>Highly conserved integral membrane protein that transport hexoses such as glucose and fructose into cells</td><td>0.0004</td><td>16.11</td><td>1.19</td><td>0.04</td><td>4.97E+07</td><td>2.16E+06</td></tr><tr><td>XP_015149397.1</td><td>416254</td><td>GEMIN5</td><td>13</td><td>Component of the survival of motor neurons (SMN) complex. The SMN complex plays a critical role in mRNA splicing through the assembly of spliceosomal small nuclear ribonucleoproteins (snRNPs), and may also mediate the assembly and transport of other classes of ribonucleoproteins</td><td>0.00089</td><td>11.12</td><td>1.02</td><td>0.029</td><td>9.70E+06</td><td>1.68E+06</td></tr><tr><td>XP_015128188.1</td><td>103906336;104639 878;1044883201</td><td>STOML2;LO C104014593; LOC101880232</td><td>NA</td><td>Mitochondrial protein that probably regulates the biogenesis and the activity of mitochondria</td><td>0.019</td><td>11.40</td><td>1.50</td><td>0.0076</td><td>2.75E+07</td><td>5.51E+06</td></tr><tr><td>NP_989530.1</td><td>374025</td><td>SMN</td><td>Z</td><td>Plays a catalyst role in the assembly of small nuclear ribonucleoproteins (snRNPs), the building blocks of the spliceosome. Thereby, plays an important role in the splicing of cellular pre-mRNAs</td><td>0.00033</td><td>12.39</td><td>1.80</td><td>0.03</td><td>2.09E+07</td><td>2.03E+06</td></tr></tbody></table>

Continued
<table><thead><tr><th rowspan="2">NCBI<br>Protein accession<br>number</th><th rowspan="2">Gene ID</th><th rowspan="2">Gene Name</th><th rowspan="2">Chr</th><th rowspan="2">Description</th><th colspan="3">Normalized weighed spectra</th><th colspan="3">Normalized precursor intensity</th></tr><tr><th>Male/Female</th><th>Average</th><th>Male/Female</th><th>Average</th><th>Male</th><th>Female</th></tr><tr><th></th><th></th><th></th><th></th><th></th><th>T-Test</th><th>Male</th><th>Female</th><th>T-Test</th><th>Male</th><th>Female</th></tr></thead><tbody><tr><td>XP_418564.2</td><td>420462</td><td>PITRM1</td><td>2</td><td>Metalloendopeptidase of the mitochondrial matrix that functions in peptide cleavage and degradation rather than in protein processing</td><td>0.0049</td><td>2.19</td><td>0.33</td><td>0.0082</td><td>7.93E+06</td><td>5.99E-01</td></tr><tr><td>P02789.2</td><td>396241</td><td>TF</td><td>9</td><td>Iron binding transport proteins responsible for the transport of iron from sites of absorption and heme degradation to those of storage and utilization. Has a role in stimulating cell proliferation</td><td>&lt;0.00010</td><td>0.00</td><td>458.31</td><td>0.00022</td><td>1.20E+00</td><td>5.81E+07</td></tr><tr><td>BAU79380.1</td><td>395882</td><td>BPIFB2</td><td>20</td><td>Member of the lipid transfer/lipopolysaccharide binding protein (LT/LBP) gene family</td><td>&lt;0.00010</td><td>0.00</td><td>80.57</td><td>&lt;0.00010</td><td>0.00E+00</td><td>4.41E+07</td></tr><tr><td>NP_001263315.1</td><td>420898</td><td>OVALX</td><td>2</td><td>Heparin-binding ov-serpin exhibiting antimicrobial activities</td><td>0.00057</td><td>0</td><td>24.69</td><td>0.049</td><td>0.00E+00</td><td>6.75E+07</td></tr><tr><td>P10184.2</td><td>416235</td><td>SPIK5</td><td>13</td><td>Serine protease inhibitor, probably important for the anti-inflammatory and/or antimicrobial protection of mucous epithelia</td><td>&lt;0.00010</td><td>0</td><td>43.50</td><td>&lt;0.00010</td><td>0.00E+00</td><td>1.87E+07</td></tr><tr><td>P67942.1</td><td>416236</td><td>SPIK7</td><td>13</td><td>Serine-type endopeptidase inhibitor</td><td>&lt;0.00010</td><td>0</td><td>27.92</td><td>0.014</td><td>0.00E+00</td><td>3.04E+07</td></tr><tr><td>CAA55385.1</td><td>396151</td><td>OVST</td><td>1</td><td>Serine-type endopeptidase inhibitor</td><td>0.0032</td><td>0</td><td>8.85</td><td>0.00055</td><td>0.00E+00</td><td>1.44E+07</td></tr><tr><td>NP_001026578.1</td><td>426937</td><td>PSMA5</td><td>26</td><td>Component of the 20S core proteasome complex involved in the proteolytic degradation of most intracellular proteins</td><td>0.00049</td><td>0</td><td>2.44</td><td>&lt;0.00010</td><td>1.20E+00</td><td>3.68E+07</td></tr><tr><td>NP_989992.1</td><td>395381</td><td>LOC395381</td><td>5</td><td>Ovomucin, the glycoprotein responsible for the gel properties of egg white</td><td>0.00023</td><td>0</td><td>14.02</td><td>0.0076</td><td>0.00E+00</td><td>1.17E+07</td></tr><tr><td>ABU24464.1</td><td>395220</td><td>ORM1</td><td>17</td><td>Functions as transport protein in the blood stream. Appears to function in modulating the activity of the immune system during the acute-phase reaction</td><td>0.006</td><td>0</td><td>16.35</td><td>0.031</td><td>0.00E+00</td><td>2.75E+07</td></tr><tr><td>XP_021263172.1</td><td>415523</td><td>SNRPA1</td><td>10</td><td>Involved in pre-mRNA splicing as component of the spliceosome</td><td>0.0096</td><td>0</td><td>3.63</td><td>&lt;0.00010</td><td>1.20E+00</td><td>8.84E+06</td></tr><tr><td>XP_009233304.1</td><td>395232</td><td>EIF4A2</td><td>9</td><td>ATP-dependent RNA helicase which is a subunit of the EIF4F complex involved in cap recognition and is required for mRNA binding to ribosome</td><td>0.0093</td><td>0</td><td>5.09</td><td>0.0028</td><td>1.20E+00</td><td>8.17E+06</td></tr></tbody><tfoot><tr><td colspan="12">Continued</td></tr></tfoot></table>
<table><caption>**Table 1.** Top forty differentially abundant proteins in male and female PGCs identified by two quantitative methods: WS and API.</caption><thead><tr><th rowspan="2">NCBI<br>Protein accession<br>number</th><th rowspan="2">Gene ID</th><th rowspan="2">Gene Name</th><th rowspan="2">Chr</th><th rowspan="2">Description</th><th colspan="3">Normalized weighed spectra</th><th colspan="3">Normalized precursor intensity</th></tr><tr><th>Male/Female</th><th>Average</th><th>Male/Female</th><th>Average</th><th>Male</th><th>Female</th></tr><tr><th></th><th></th><th></th><th></th><th></th><th>T-Test</th><th>Male</th><th>Female</th><th>T-Test</th><th>Male</th><th>Female</th></tr></thead><tbody><tr><td>XP_015708132.1</td><td>419058</td><td>PPME1</td><td>1</td><td>Protein phosphatase<br>methylesterase<br>localized to the<br>nucleus. The<br>encoded protein<br>acts on the protein<br>phosphatase-2A<br>catalytic subunit and<br>supports the ERK<br>pathway through<br>dephosphoryla-<br>tion of regulatory<br>proteins</td><td>0.032</td><td>0</td><td>3.95</td><td>0.00029</td><td>0.00E+00</td><td>1.76E+07</td></tr><tr><td>NP_990025.1</td><td>395430</td><td>CDC37</td><td>30</td><td>Molecular chaper-<br>one with specific<br>function in cell<br>signal transduc-<br>tion, binds to<br>numerous kinases<br>and promotes their<br>interaction with the<br>Hsp90 complex</td><td>0.026</td><td>0</td><td>1.25</td><td>0.028</td><td>8.27E-01</td><td>4.63E+06</td></tr><tr><td>NP_001026411.1</td><td>423947</td><td>ACADSB</td><td>6</td><td>Short and branched<br>chain specific acyl-<br>CoA dehydrogenase<br>that catalyzes the<br>removal of one<br>hydrogen from C-2<br>and C-3 of the fatty<br>acyl-CoA thioester,<br>resulting in the<br>formation of trans-<br>2-enoyl-CoA</td><td>0.028</td><td>0</td><td>2.71</td><td>0.0012</td><td>4.31E-01</td><td>6.19E+06</td></tr><tr><td>NP_001007840.1</td><td>416955</td><td>YWHAH</td><td>15</td><td>Adapter protein<br>implicated in the<br>regulation of a<br>large spectrum of<br>both general and<br>specialized signaling<br>pathways</td><td>0.022</td><td>0.02</td><td>2.16</td><td>0.025</td><td>2.02E+06</td><td>1.94E+07</td></tr><tr><td>BAM13279.1</td><td>420897</td><td>OVALY</td><td>2</td><td>Serine-type endo-<br>peptidase inhibitor</td><td>&lt;0.00010</td><td>1.83</td><td>95.43</td><td>0.016</td><td>8.35E+06</td><td>8.44E+07</td></tr><tr><td>XP_014135657.1</td><td>422354</td><td>RAB39B</td><td>4</td><td>Small GTPase<br>Rab involved in<br>autophagy and key<br>regulators of intra-<br>cellular membrane<br>trafficking, from<br>the formation of<br>transport vesicles<br>to their fusion with<br>membranes</td><td>0.041</td><td>0.03</td><td>0.93</td><td>0.026</td><td>1.39E+07</td><td>7.43E+07</td></tr><tr><td>Q7LZQ2.1</td><td>NA</td><td>FG</td><td>NA</td><td>Protein with bac-<br>teriolytic function,<br>also associated with<br>the monocyte-mac-<br>rophage system and<br>immunomodulation</td><td>0.0024</td><td>0.37</td><td>10.39</td><td>0.0022</td><td>4.31E-01</td><td>4.28E+07</td></tr><tr><td>XP_010117667.1</td><td>396228</td><td>SERPINH1</td><td>1</td><td>Serine proteinase<br>inhibitor that func-<br>tions as a chaperone<br>in the biosynthetic<br>pathway of collagen</td><td>0.035</td><td>0.35</td><td>9.93</td><td>0.001</td><td>1.20E+00</td><td>2.38E+07</td></tr><tr><td>NP_035318.1</td><td>417716</td><td>PSMC2</td><td>1</td><td>Component of the<br>26S proteasome,<br>a multiprotein<br>complex involved in<br>the ATP-dependent<br>degradation of ubiq-<br>uitinated proteins</td><td>0.0018</td><td>0.74</td><td>5.25</td><td>0.00015</td><td>1.20E+00</td><td>2.18E+07</td></tr><tr><td>XP_015726668.1</td><td>107317943</td><td>DTYMK</td><td>9</td><td>Protein with thymi-<br>dylate kinase activ-<br>ity that catalyzes the<br>conversion of dTMP<br>to dTDP</td><td>0.008</td><td>1.12</td><td>7.92</td><td>0.00012</td><td>7.66E-01</td><td>4.23E+06</td></tr></tbody></table>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Western blot analysis of HSDL2 (at the top) and IGF2BP1 (at the bottom) in chicken male and female PGCs derived in vitro. Western blot analysis was performed on the protein extracts from 5 individual female (ZW₁⁻₅) and male (ZZ₁⁻₅) PGCs cultures as described in the Sect. *Reverse transcription quantitative real time polymerase chain reaction (RT-QPCR)*. The pictures of nitrocellulose membrane with blotted proteins after total protein staining and after incubation with HSDL2 and IGF2BP1 are represented. For IGF2BP1 two presented parts of the blot correspond to the same blot. The full-length blots for HSDL2 and IGF2BP1 are presented on the Supplementary Fig. S2. The bands of 45 kDa and of 63 kDa corresponding to HSDL2 and IGF2BP1 proteins respectively are indicated with arrows. Mouse liver (ml); chicken liver (cl) present positive controls. The results on the graphs represent the means ± SEM, *p* < 0.001 (***), *p* < 0.02 (*).</figcaption></figure>

0.7% (6 out of 764 unique chicken IDs found for all identified proteins) of Z-linked proteins did not change in abundance according to sex. We did not identify any proteins encoded by W-linked genes among proteins more abundant in female PGCs. This was, however, not surprising, considering that the W chromosome is very small and poor in protein-coding genes.

**Overall transcriptome differences do not explain proteome differences found in male and female PGCs.** The male to female abundance ratios were higher than two-fold for 68% of the differentially expressed proteins encoded by Z-linked genes. For several of these proteins the ratio was much higher than two, suggesting a positive regulation in addition to a double dose effect.

We next compared the transcript and protein expression ratios for 5 selected proteins encoded by Z-linked genes and 13 selected proteins encoded by autosomal genes (Fig. 2). The autosomal genes *DAZL* and *Ribosomal Protein L30* (*RPL30*; a germ cell marker and a housekeeping gene, respectively) were chosen as controls, and, as expected, no changes were observed neither at the protein level nor at the transcript level. We observed that for proteins more abundant in male PGCs, transcripts gene expression levels were also significantly higher (6 genes), the single exception being the nucleoporin 210 (NUP210). In contrast, no differences in transcript expression were observed in genes corresponding to female overabundant proteins (9 genes), except for the core histone macro-H2A (H2AFY). For Z-linked genes analysed (*DDX4*, *Tallin-1*; *TLN1*, *Tetratricopeptide Repeat Domain 37*; *TTC37*, *Endoplasmic Reticulum Metallopeptidase 1*; *ERMP1* and *Nucleoporin 155*; *NUP155*), transcript level differences appeared to be due to the previously mentioned lack of Z chromosome dosage compensation in ZZ PGCs, and not due to gene expression upregulation. For these Z chromosome genes and for the autosomal gene
<table><caption>**Table 2.** List of sex chromosome linked genes encoding proteins identified in cultured in vitro chicken PGCs.</caption><thead><tr><td>NCBI gene name</td><td>NCBI gene ID</td><td>Sex of preferential abundance</td><td>Chromosome</td></tr></thead><tbody><tr><td>TXNL1</td><td>426854</td><td>Male</td><td>Z</td></tr><tr><td>TXN</td><td>396437</td><td>Male</td><td>Z</td></tr><tr><td>TTC37</td><td>427114</td><td>Male</td><td>Z</td></tr><tr><td>TLN1</td><td>395194</td><td>Male</td><td>Z</td></tr><tr><td>TBCA</td><td>416367</td><td>Male</td><td>Z</td></tr><tr><td>TARS</td><td>427427</td><td>Male</td><td>Z</td></tr><tr><td>SMN</td><td>374025</td><td>Male</td><td>Z</td></tr><tr><td>PTGR1</td><td>427337</td><td>Male</td><td>Z</td></tr><tr><td>NUP155</td><td>427443</td><td>Male</td><td>Z</td></tr><tr><td>IQGAP2</td><td>427211</td><td>Male</td><td>Z</td></tr><tr><td>HSDL2</td><td>100858057</td><td>Male</td><td>Z</td></tr><tr><td>HINT2</td><td>395424</td><td>Male</td><td>Z</td></tr><tr><td>ERMP1</td><td>426644</td><td>Male</td><td>Z</td></tr><tr><td>DDX4</td><td>395447</td><td>Male</td><td>Z</td></tr><tr><td>AD012Z</td><td>407092</td><td>Male</td><td>Z</td></tr><tr><td>ACO1</td><td>373916</td><td>Male</td><td>Z</td></tr><tr><td>ACAA2</td><td>426847</td><td>Male</td><td>Z</td></tr><tr><td>RPL17</td><td>426845</td><td>Female</td><td>Z</td></tr><tr><td>RPS6</td><td>396148</td><td>none</td><td>Z</td></tr><tr><td>RPS23</td><td>427323</td><td>none</td><td>Z</td></tr><tr><td>LOC430766/LOC107049323</td><td>430766/107049323</td><td>none</td><td>W/Z</td></tr><tr><td>ISOC1</td><td>415601</td><td>none</td><td>Z</td></tr><tr><td>HNRNPK/HNRNPKL</td><td>427458/426516</td><td>none</td><td>Z/W</td></tr><tr><td>ELAVL2</td><td>770158</td><td>none</td><td>Z</td></tr><tr><td>CLTA</td><td>427284</td><td>none</td><td>Z</td></tr><tr><td>ATP5I</td><td>769146</td><td>none</td><td>Z</td></tr><tr><td>ATP5A1Z/ATP5A1W</td><td>374159/431564</td><td>none</td><td>Z/W</td></tr></tbody></table>

WDR1, the increased transcript levels in male PGCs does not appear to be the only mechanism responsible for increased protein accumulation since their protein male/female ratio was much higher than that of the transcript ratios, with the exception of DDX4 (Fig. 2). Instead, post-transcriptional and post-translational mechanisms are likely to be activated to create the higher abundance of these proteins, and this can also explain the high protein male/female ratio for NUP210, whose transcript was not differentially expressed in male and female PGCs. In agreement with this observation, a previous study that compared gene expression and protein abundance in male and female adult chicken tissues described that 30% of all proteins encoded from Z-linked genes showed a significant change in the male/female ratio compared with the corresponding ratio at the RNA level³¹,³².

The differential expression of both autosomal and Z-linked genes in male PGCs in the absence of surrounding gonadal somatic cells or sexual hormones suggests that their molecular sex differentiation is cell autonomous, as previously suggested in studies performed in murine PGCs and chicken embryonic tissues³²⁻³⁴. The substantial sexual differences in the protein abundances in pregonadal PGCs observed in this study may be the consequence of the cumulative effect of uncompensated Z-linked gene products. This observation may be linked to the proposed role of uncompensated dosage of Z-linked genes in the molecular sexual identity of early chicken embryo³².

**Female PGCs accumulate proteins from the culture medium.** Several of the proteins that were highly overabundant in female PGCs (Table 1, supplementary Table S2c), transferrin, (TF), BPIFB2, SPINK5, SPINK7, ovostatin (OVST), ovomucin, alpha subunit (LOC395381, OVOA), orosomucoid 1 (ovoglycoprotein) ORM1, ovalbumin-related protein Y (OVALY, SERPINB14B), ovalbumin-related protein X (OVALX (SERPINB14C), VMO1, ovalbumin (OVAL, SERPINB14), were not expected to be present in PGCs, as they are oviduct- or serum-specific proteins. Using RT-qPCR analysis, we observed that female PGCs have no or very low mRNA expression for these proteins (Supplementary Fig. S4). We, thus, suspected that female PGCs were able to selectively pinocytosis and/or transport these proteins from the extracellular culture medium. To confirm this hypothesis, we analysed two commercial protein products added to cell culture medium: ovalbumin (albumin from hen egg white; SIGMA) and chicken serum (SIGMA) by nanoLC-MS/MS analysis (Supplementary Methods). All the proteins detected in female PGCs for which we could not find corresponding mRNA expression were identified by this analysis, with the exception of vitelline membrane outer layer 1 homolog (VMO1) and ORM1 (Table 3, Supplementary Table S4). Among these proteins BPIFB2 (also known as protein TENP, G2), OVOA (also known as LOC395381), SPINK5, SPINK7, OVAL (also known as SERPINB14), OVALX, OVALY, OVST were identified in the commercial ovalbumin product while TF was found in both albumin from hen egg white and chicken serum. VMO1, another component of the egg white, was not identified in the ovalbumin
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Transcript and protein male/female expression fold change (FC) in cultured in vitro chicken PGCs. Transcript and protein abundances were evaluated by RTQPCR and two label-free proteomic quantitative methods (Spectral Counting and Average Precursor Intensity) respectively, as described in the Sect. *Label-free protein quantification*. Weighted spectra (WS); Average Precursor Intensity (API); *p* < 0.05 (*); *p* < 0.01 (**); *p* < 0.001 (***); not determined in female PGCs (ndf); not determined in male PGCs (ndm).</figcaption></figure>

fraction used in PGCs culture medium, probably because of its lower abundance in egg white compared with the rest of the proteins which may have prevented its detection³⁵. In view of the fact that these components were also present in the male PGC culture medium, these results suggest that female PGCs are able to differentially import proteins from the medium compared with male PGCs, indicating an intrinsic capacity of female germ cells to transport and accumulate proteins.

**Post-transcriptional RNA and post-translational protein processing are enhanced in female chicken PGCs.** We carried out a functional pathway analysis on proteins found in male and female PGCs. This functional analysis suggested increased transcriptional and translational processing in female PGCs compared to male PGCs, indicated by the significant overrepresentation of the “Biological Process” (BP) GO terms “RNA transcription and its regulation”, “RNA biosynthetic processes”, “regulation of RNA biosynthesis” and “regulation of RNA metabolism” (Fig. 3a; Supplementary Table S5). Based on the functional analysis and literature searches on overabundant proteins in female PGCs, we elaborate on how these proteins are interrelated, as summarized in Fig. 4a and as detailed below.

*Germ cell phenotype of female PGCs is promoted by an enrichment in proteins related to RNA metabolism and protein processing.* In germ cells (and somatic cells in a lesser extent) cytoplasmic untranslated mRNAs accumulate in ribonucleoprotein (RNP) complexes, which are important posttranscriptional regulators of gene expression during development with functions in mRNA degradation, storage, sorting, and transport. A subset of more abundant female PGCs proteins has been found to be components of RNP granules in murine PGCs-like cells including: insulin growth factor 2 binding protein 1 (IGF2BP1, also known as IMP1), Interleukin Enhancer Binding Factor 2 (ILF2), Eukaryotic Translation Initiation Factor 4A3 (EIF4A3) and Heterogeneous Nuclear Ribonucleoprotein K (HNRNPK)³⁶. In these cells, IGF2BP1 is associated with the germ cell marker DAZL in RNP granules, whose role is to carry out translational repression of differentiation, pluripotency and apoptosis related RNAs³⁶,³⁷. Interestingly, female PGCs were enriched in proteins involved in secretion and the ubiquitin-proteasome degradation system compared to male PGCs (RABB1, RAB11B, UBE2V2, UCHL3, PSMA1, PSMA5, PSMA7, PSME3, PSMC2, PSMD12, COPS8). The respective mRNAs of these proteins are known to be components of IGF2BP1-containing RNP granules³⁷. These results suggest that in female PGCs, there is a preponderant cooperation between RNP granules, the protein-processing organelle machinery and the ubiquitin-proteasome degradation system to govern RNA metabolism, protein synthesis, processing, and degradation.

Another protein related to the promotion of 'germ cell-ness' mediated by IGF2BP1 signalling, the post-transcriptional regulator LIN-28 homolog A (LIN28A), was overabundant in female PGCs. The expression
<figure><img src="image_6.png" /><figcaption>b</figcaption></figure>

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

**Figure 3.** Wang GO terms distance clustering heat map plots for significantly enriched GO terms in differentially abundant proteins in male and female cultured PGCs. DA, GO terms enriched in differentially abundant proteins using as background all identified proteins except components of culture medium specified in the Sect. *Female PGCs accumulate proteins from the culture medium* and that are not expressed in PGCs; male, GO terms significantly enriched in proteins preferentially abundant in male PGCs; female, GO terms significantly enriched in proteins preferentially abundant in female PGCs; IC, informational content. Functional analysis and its visualisation were performed using ViSEAGO²⁴. (a), Go terms Biological Process (BP) clustering heatmap plots; (b), Go terms Molecular Function (MF) clustering heatmap plots.

of LIN28A is restricted to PGCs and premeiotic germ cells in foetal mouse and human ovaries³⁸,³⁹, where it is known to promote pluripotency and to be involved in PGC formation. LIN28A has been shown to interact mechanistically and functionally with IGF2BP1 to regulate the insulin-like growth factor 2-mammalian target of rapamycin (IGF2-mTOR) signalling in neural progenitor cells, a pathway that is also essential for chicken PGCs survival and proliferation⁵,⁴⁰. Moreover, LIN28A promotes IGF2BP1 expression, allowing post-transcriptional repression of differentiation genes in foetal mouse neural stem cells⁴¹. Both proteins are also indirectly related through the regulation of pluripotency gene SOX2 in neural precursor cells and in murine in vitro derived PGC-like cells³⁶,⁴². Accordingly, LIN28A and IGF2BP1 could cooperate in chicken PGCs in maintaining the fine balance between germ cell-ness and pluripotency and this mechanism could be enhanced in chicken female PGCs. Sex differences observed in LIN28A protein expression may indicate differences in the pluripotency state of male and female PGCs.

In addition to the possible role of IGF2BP1 in the regulation of pluripotency in chicken PGCs IGF2BP1 may also promote female PGCs migration and gonad colonisation through protein phosphatase methyl esterase 1 (PPME1) overabundant in female PGCs, as PPME1 has been found to be novel effector of IGF2BP1 in promoting migration and invasion of human choriocarcinoma cells⁴³.

*Female PGCs may regulate the activity of the ubiquitin-proteasome system through the activity of COP9 signalosome.* The COP9 signalosome (CSN) is a highly conserved protein complex associated with de-ubiquitina-
<table><caption>**Table 3.** Proteins overabundant in female PGCs and identified by LC-MS/MS Analysis of chicken ovalbumin and chicken serum. **Mean emPAI for 35 proteins from group a1 identified as Ovalbumin.</caption><caption>Albumin from hen egg white (SIGMA)</caption><thead><tr><th>Group ID</th><th>Subgroup ID</th><th>subgroup protein ID</th><th>UniProt Protein ID</th><th>Description</th><th>Gene Name</th><th>emPAI</th></tr></thead><tbody><tr><td>a1</td><td>a1.a1</td><td>a1.a1.a1</td><td>A0A2H4Y814</td><td>CHICK OVA</td><td>OVA</td><td>7.6E+08</td></tr><tr><td></td><td></td><td>35 other proteins for group a1</td><td></td><td>CHICK OVA</td><td>OVA</td><td>5.4E+08**</td></tr><tr><td>a2</td><td>a2.a1</td><td>a2.a1.a1</td><td>P01014</td><td>CHICK Ovalbumin-related protein Y</td><td>OVALY</td><td>153</td></tr><tr><td></td><td>a2.a2</td><td>a2.a2.a1</td><td>A0A1D5P531</td><td>CHICK Uncharacterized protein</td><td>OVALX</td><td>3.6</td></tr><tr><td></td><td>a2.a2</td><td>a2.a2.a2</td><td>A0A1D5PI58</td><td>CHICK Uncharacterized protein</td><td>OVALX</td><td>3.2</td></tr><tr><td>a3</td><td>a3.a1</td><td>a3.a1.a1</td><td>A0A1D5P4L7</td><td>CHICK Ovotransferrin</td><td>N/A</td><td>2</td></tr><tr><td>a3</td><td>a3.a2</td><td>a3.a2.a1</td><td>P02789</td><td>CHICK Ovotransferrin</td><td>TF</td><td>1.5</td></tr><tr><td>a4</td><td>a4.a1</td><td>a4.a1.a1</td><td>P01005</td><td>CHICK Ovomucoid</td><td>N/A</td><td>55</td></tr><tr><td>a4</td><td>a4.a1</td><td>a4.a1.a2</td><td>A0A1D5NYW</td><td>CHICK Ovomucoid</td><td>SPINK7</td><td>35</td></tr><tr><td>a4</td><td>a4.a1</td><td>a4.a1.a3</td><td>B6V1G0</td><td>CHICK Ovomucoid</td><td>N/A</td><td>35</td></tr><tr><td>a5</td><td>a5.a1</td><td>a5.a1.a1</td><td>Q98UI9</td><td>CHICK Mucin-5B</td><td>MUC5B</td><td>0.2</td></tr><tr><td>a5</td><td>a5.a1</td><td>a5.a1.a2</td><td>A0A1D5PPR9</td><td>CHICK Uncharacterized protein</td><td>LOC395381</td><td>0.2</td></tr><tr><td>a6</td><td>a5.a1</td><td>a5.a1.a3</td><td>F1NZY2</td><td>CHICK Uncharacterized protein</td><td>LOC395382</td><td>0.2</td></tr><tr><td>a5</td><td>a6.a1</td><td>a6.a1.a1</td><td>P20740</td><td>CHICK Ovostatin</td><td>N/A</td><td>0.3</td></tr><tr><td>a6</td><td>a6.a1</td><td>a6.a1.a2</td><td>A0A1D5NT83</td><td>CHICK Uncharacterized protein</td><td>OVST</td><td>0.3</td></tr><tr><td>a6</td><td>a6.a1</td><td>a6.a1.a3</td><td>A0A1D5P3R8</td><td>CHICK Uncharacterized protein</td><td>OVST</td><td>0.3</td></tr><tr><td>a6</td><td>a6.a1</td><td>a6.a1.a4</td><td>A0A1L1RUH2</td><td>CHICK Uncharacterized protein</td><td>OVST</td><td>0.3</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.a1</td><td>A0A146J2W9</td><td>CHICK Protein TENP</td><td>TENP</td><td>0.4</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.a2</td><td>A0A146J2X3</td><td>CHICK Protein TENP</td><td>TENP</td><td>0.4</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.a3</td><td>A0A146J2Y8</td><td>CHICK Protein TENP</td><td>TENP</td><td>0.4</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.a4</td><td>A0A146J2Z8</td><td>CHICK Protein TENP</td><td>TENP</td><td>0.4</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.a5</td><td>A0A146J302</td><td>CHICK Protein TENP</td><td>TENP</td><td>0.4</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.a6</td><td>A0A146J377</td><td>CHICK Protein TENP</td><td>TENP</td><td>0.4</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.a7</td><td>A0A146J3I5</td><td>CHICK Protein TENP</td><td>TENP</td><td>0.4</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.a8</td><td>A0A1L1RMA4</td><td>CHICK Protein TENP</td><td>BPIFB2</td><td>0.4</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.a9</td><td>I0J173</td><td>CHICK OvoglobulinG2</td><td>G2</td><td>0.4</td></tr><tr><td>a9</td><td>a9.a1</td><td>a9.a1.b10</td><td>I0J174</td><td>CHICK OvoglobulinG3</td><td>G2</td><td>0.4</td></tr><tr><td colspan="7"><strong>chicken serum (SIGMA)</strong></td></tr><tr><td>a3</td><td>a3.a1</td><td>a3.a1.a1</td><td>A0A1D5P4L7</td><td>CHICK Ovotransferrin</td><td>TF</td><td>182.3</td></tr><tr><td>a3</td><td>a3.a2</td><td>a3.a2.a1</td><td>Q4ADJ6</td><td>CHICK Ovotransferrin</td><td>TFEW</td><td>111.9</td></tr><tr><td>b17</td><td>b17.a1</td><td>b17.a1.a1</td><td>P10184</td><td>CHICK Ovoinhibitor OS</td><td>OIH</td><td>9.0</td></tr><tr><td>b17</td><td>b17.a2</td><td>b17.a2.a1</td><td>F1NMN2</td><td>CHICK Uncharacterized protein</td><td>SPINK5</td><td>9.0</td></tr></tbody></table>

tion activity and protein kinase activities capable of phosphorylating important signalling regulators. NGG1 Interacting Factor 3 Like 1 (NIF3L1) and COP9 Signalosome Subunit 8 (COPS8/CSN8) were overabundant in female PGCs. NIF3L1 is a transcriptional corepressor through its interaction with CSN2/COPS2 (a subunit of the COP9 signalosome complex or CSN), negatively regulating the expression of OKT-3/4, another master regulator of pluripotency and partner of Sox2, and promoting neural differentiation³⁶,⁴⁴. Female PGCs might have an enhanced CSN activity in order to control ubiquitin-proteasome-mediated protein degradation to regulate their state of differentiation.

*Female PGCs control transcription through an enrichment in nucleosome proteins.* Several nucleosome proteins involved in nuclear architecture, chromatin state and epigenetic regulation of transcription were overabundant in female PGCs. Among them, three subtypes of histone H1 (H1.10 Linker Histone; H1FX, Histone H1.03;
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Crosstalk of pathways involving orthologs of chicken proteins preferentially abundant in chicken female (a) and male (b) PGCs based on the literature (for references see Sects. 3.6 and 3.7). Proteins in pink and blue correspond to the orthologs of proteins identified as overabundant in chicken female and male PGCs respectively. Pointed and blunted arrow heads represent activation and inhibition respectively.</figcaption></figure>

HIST1H103, histone cluster 1, H1.11R) and the core histone macroH2A2 (MACROH2A2). Although a higher abundance of H1 linker histone subtypes in female PGCs in comparison to male PGCs after gonad differentiation have been described in the mouse⁴⁵, ours is the first report of sex-related H1 differences in germ cells before gonad differentiation. These results indicate that sexually dimorphic nucleosome regulation of transcription is present in chicken PGCs, before PGC differentiation begins.
**Male chicken PGCs proteome suggests common features with spermatogonia.** Functional analysis of differentially abundant proteins in male PGCs pointed to cytoskeleton organisation differences and increased metabolism, as suggested by the enriched BP GO terms “spindle localisation” and “aerobic respiration”, as well as with the MF GO terms “ATP-ase activity coupled to substances”, “transmembrane movement” and “transmembrane transport” (Fig. 3b). The crosstalk of pathways in which these proteins are involved and which may operate in chicken ZZ PGCs based on the literature is summarized in Fig. 4b.

These GO terms reflect the enrichment of the male chicken PGC proteome in proteins known to play an essential role in the post-natal establishment and maintenance of spermatogonia in several species. PGCs, as opposed to gonocytes (quiescent gonadal embryonic germ cells), share some properties with spermatogonia, in particular the ability to proliferate and migrate towards the somatic niche. Our results suggest that these proteins are necessary in pre-gonadal male PGCs for these processes and that their role may, therefore, be in part similar to that in spermatogonia maintenance.

*The germ cell phenotype of male PGCs is promoted by an enrichment in proteins related to the ubiquitin-proteasome protein degradation system.* Several proteins associated with the Ubiquitin-26S proteasome complex were preferentially abundant in male PGCs, namely Ubiquitin Protein Ligase E3 Component N-Recognin 4 (UBR4), WWE Domain Containing E3 Ubiquitin Protein Ligase 1 (HUWE1), non-ATPase regulatory subunit 2 (PSMD2), 26S proteasome co-factor thioredoxin-like protein 1 (TXNL1)⁴⁶,⁴⁷, and a chicken orthologue of Ubiquitin Specific Peptidase 14 (USP14)⁴⁸. A number of these proteins are known to specifically direct protein degradation in male spermatogonia or sperm cells. HUWE1, known to ubiquitinate core histones, has been recently reported to be essential to the spermatogonia establishment and maintenance by preventing genomic instability via suppressing the DNA Damage Response (DDR) through blocking the accumulation of Histone H2AX in mice⁴⁹. Moreover, UBR4, an E-3 ligase having “N-end rule” specificity has been shown to be involved in the degradation of histone-binding complex HAT1/RBBP4/RBBP7 in Drosophila and mouse myofibers⁵⁰. The proteins in this complex are known to be involved in DNA repair⁵¹⁻⁵³. According to these results, male PGCs could promote the ubiquitin-proteasome system to control core histone half-life. Additionally, UBR4 and USP14 are essential for male but not for female fertility in several species⁵⁴⁻⁵⁶.

*Spinal motor neuron (SMN) is a central functional node in the differential male PGC proteome.* The protein, survival motor neuron (SMN), as well as its partner and translational activator, GEMIN5⁵⁷,⁵⁸, were significantly more abundant in male PGCs. Using literature mining, SMN appears to be a central functional node in the network of the overabundant proteins in male PGCs. SMN protein seems to interrelate different processes essential for male fertility involving intracellular trafficking, protein degradation and translational silencing⁵⁷, and to directly interact with several of its components. SMN is able to interact with the ubiquitin-proteasome-degradation system to control protein degradation including E3-ligases⁵⁹⁻⁶¹. In our dataset, this relationship is suggested by the co-accumulation of SMN with the E3 ligases HUWE1 and UBR4 in male PGCs.

SMN is present in the cell nucleus, localized to subnuclear bodies near coiled bodies containing high concentrations of small ribonucleoproteins (snRNPs), and is involved in the biogenesis of snRNPs, the building blocks of the spliceosome⁶². One of snRNPs components, small nuclear ribonucleoprotein polypeptide F (SNRPF), was also overabundant in male PGCs. SMN is essential for male but not female reproductive function and is involved in the germ cell development and survival and spermatogonia maintenance in the mouse⁶³,⁶⁴. In Drosophila, SMN expression is essential for survival and proliferation of male germline stem cells while it is not expressed in spermatocytes or the somatic niche⁶⁵.

In Drosophila, SMN is implicated in the regulation of germline nuclear organization through the connection of U bodies (structures rich in uridine-rich small nuclear ribonucleoproteins involved in key steps of pre-mRNA processing) and processing bodies (P-bodies; mRNA-degrading RNP granules)⁶⁶. Other two proteins known to be involved in the formation of P-bodies were more abundant in male PGCs, namely E3 SUMO-protein ligase RANBP2 (also known as NUP358)⁶⁷,⁶⁸ and CCR4-NOT transcription complex subunit 1(CNOT1). In mice, RANBP2 is known to be essential for male fertility⁶⁸ and CNOT1, interacting with NANOS2 in both male gonocytes and spermatogonial stem cells (SSCs) maintenance, is involved in suppression of specific RNAs promoting male differentiation fate of germ cells and blocking the female program⁶⁹. This suggests that in chicken male PGCs, SMN, CNOT1 and RANBP2 could be connected through their involvement in the function of P bodies and, consequently, in translational silencing.

SMN is a repressor of the cell-fate Notch signalling pathway⁷⁰. The involvement of Notch signalling in male PGCs homeostasis was indicated with the co-accumulation of Strawberry Notch1 (SBNO1), a helicase-related nuclear factor which is known to promote the proliferation of SCCs in mouse neonatal testis⁷¹. Interestingly, Huwe1 has also been found to regulate Notch signalling in neural progenitors⁷².

Our results indicate that the protein SMN was highly inter-connected with other proteins expressed more abundantly in male PGCs (Fig. 4b), suggesting its central role in male PGCs homeostasis.

*Intracellular trafficking and cytoskeletal remodelling are enhanced in male PGCs.* Pre-gonadal male and female PGCs migrate towards and colonize the gonads during embryonic development. These two fundamental PGCs functions require changes in the actin cytoskeleton as well as enhanced intracellular trafficking. Our results indicate that male PGCs are enriched in proteins related to these processes. Several proteins related to the actin cytoskeleton and cytokinesis were overrepresented in the male PGCs proteome, namely WDR1 (also known as AIP1)⁷³, SMN⁷⁴,⁷⁵, Talin 1 (TLN1) and Talin 2 (TLN2)⁷⁶⁻⁷⁸. WDR1 is essential for the migration of germ cells towards the testicular cords in the post-natal mouse testis⁷³. WDR1 and talins are known to be functionally interconnected during cytoskeletal re-arrangements. For instance, the redistribution of TLN1 and TLN2 is regu-
lated by WDR1 in a process mediated by integrins during platelet activation in the mouse<sup>79</sup>. Integrin β1 receptors are involved in chicken PGCs migration during embryonic development<sup>80</sup>. Integrin-mediated signaling is mediated by the E3 ubiquitin ligase enzyme UBR4, which was also more abundant in male PGCs together with the protein Clathrin. Both proteins form meshwork structures involved in membrane morphogenesis and cytoskeletal organization<sup>81</sup>. Directly linked with this, the biological process “spindle localization” was enriched in male PGCs according to functional analysis, which involved Dynein Cytoplasmic 1 Heavy Chain 1 (DYNC1H1) and Dynactin Subunit 1 (DCTN1). Indeed, these proteins and UBR4 have been related to equivalent processes in neural progenitors<sup>82,83</sup>.

An enhanced membrane and intracellular trafficking, which is a feature compatible with the pro-migrating phenotype of PGCs, was also suggested in male PGCs by proteomic results. For instance, SMN, HUWE1 and UBR4 are involved in intracellular trafficking and promote endocytosis<sup>84,85</sup>. Moreover, male PGCs were enriched in several endosome components including the proteins Clathrin, Vesicle Associated Membrane Protein 3 (VAMP3), DYNC1H1, DCTN1 and DnaJ Heat Shock Protein Family (Hsp40) Member C13 (DNAJC13/RME8). Moreover, Clathrin-mediated trafficking is a major mechanism for internalization of membrane receptors-ligand complexes<sup>86</sup> including those of Insulin<sup>87</sup>, fibroblast growth factor (FGF)<sup>88</sup> and bone morphogenetic proteins (BMPs)<sup>89</sup>, which are components of culture medium essential for survival and proliferation of PGCs<sup>5</sup>. Thus, enhanced trafficking may have a consequence on the sensitivity of PGCs to these molecules and on their signaling as it has been described in other cell models<sup>90-92</sup>.

In conclusion, male and female chicken PGCs propagated in vitro showed early proteome differences in agreement with a cell autonomous sex determination during early embryonic development which is independent of sex hormones. Moreover, our proteomic results show that the absence of a general mechanism of dosage compensation of Z-linked genes in chicken extends to pre-gonadal PGCs. Indeed, this phenomenon had been documented at transcriptional and translational levels in adult chicken tissues, and at a transcriptional level for embryonic tissues, but not in embryonic germ cells. The cumulative effect of the double dosage of proteins encoded by Z-linked genes in male PGCs may be one specific factor governing cell autonomous sex determination of chicken PGCs. In consequence, ZZ and ZW PGCs may develop different molecular strategies to achieve equivalent functions such as the blocking of sex-specific differentiation and repression of opposite sex developmental pathways. Furthermore, specific differences in ZZ and ZW PGCs overabundant protein sets indicate an early programming of male and female germ cell fates. This study provides the first molecular support that PGCs may be subjected to sex specific regulation independent of the in vitro environment. Sex specific adaptation of in vitro culture conditions for chicken PGCs would be useful to improve quality of these cells for the management of avian genetic diversity. Based on our findings, Insulin, FGF and BMP signaling and the pluripotency state of PGCs are potential candidate targets for further investigation.

## Data availability

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE89 partner repository with the dataset identifier PXD022452 and 10.6019/PXD022452.

Received: 26 February 2021; Accepted: 6 September 2021

Published online: 28 September 2021

## References

1. Blesbois, E. et al. Semen cryopreservation for ex situ management of genetic diversity in chicken: Creation of the French avian cryobank. *Poult. Sci.* **86**, 555–564. https://doi.org/10.1093/ps/86.3.555 (2007).

2. Nakamura, Y., Kagami, H. & Tagami, T. Development, differentiation and manipulation of chicken germ cells. *Dev. Growth Differ.* **55**, 20–40. https://doi.org/10.1111/dgd.12026 (2013).

3. Song, Y. et al. Characteristics of long-term cultures of avian primordial germ cells and gonocytes. *Biol. Reprod.* **90**, 15. https://doi.org/10.1095/biolreprod.113.113381 (2014).

4. Nandi, S. et al. Cryopreservation of specialized chicken lines using cultured primordial germ cells. *Poult. Sci.* **95**, 1905–1911. https://doi.org/10.3382/ps/pew133 (2016).

5. Whyte, J. et al. FGF, insulin, and SMAD signaling cooperate for avian primordial germ cell self-renewal. *Stem Cell Rep.* **5**, 1171–1182. https://doi.org/10.1016/j.stemcr.2015.10.008 (2015).

6. van de Lavoir, M. C. et al. Germline transmission of genetically modified primordial germ cells. *Nature* **441**, 766–769. https://doi.org/10.1038/nature04831 (2006).

7. Woodcock, M. E., Idoko-Akoh, A. & McGrew, M. J. Gene editing in birds takes flight. *Mammalian Genome Official J. Int. Mammalian Genome Soc.* **28**, 315–323. https://doi.org/10.1007/s00335-017-9701-z (2017).

8. Whyte, J., Blesbois, E. & McGrew, M. *Increased sustainability in poultry production: New tools and resources for genetic management.*, Vol. 31 214–230 (CABI Publishing, 2016).

9. Ortega-Recalde, O., Goikoetxea, A., Hore, T. A., Todd, E. V. & Gemmell, N. J. The genetics and epigenetics of sex change in fish. *Ann. Rev. Animal Biosci.* https://doi.org/10.1146/annurev-animal-021419-083634 (2019).

10. Taketo, T. The role of sex chromosomes in mammalian germ cell differentiation: Can the germ cells carrying X and Y chromosomes differentiate into fertile oocytes?. *Asian J. Androl.* **17**, 360–366. https://doi.org/10.4103/1008-682x.143306 (2015).

11. Hoekstra, H. E. & Hoekstra, J. M. An unusual sex-determination system in South American field mice (Genus Akodon): The role of mutation, selection, and meiotic drive in maintaining XY females. *Evol. Int. J. Organic Evol.* **55**, 190–197 (2001).

12. Vernet, N. et al. Mouse Y-linked Zfy1 and Zfy2 are expressed during the male-specific interphase between meiosis I and meiosis II and promote the 2nd meiotic division. *PLoS Genet.* **10**, e1004444. https://doi.org/10.1371/journal.pgen.1004444 (2014).

13. Tagami, T. et al. Differentiation of female primordial germ cells in the male testes of chicken (*Gallus gallus* domesticus). *Mol. Reprod. Dev.* **74**, 68–75. https://doi.org/10.1002/mrd.20499 (2007).

14. Macdonald, J., Glover, J. D., Taylor, L., Sang, H. M. & McGrew, M. J. Characterisation and germline transmission of cultured avian primordial germ cells. *PLoS ONE* **5**, e15518. https://doi.org/10.1371/journal.pone.0015518 (2010).

15. Park, T. S. & Han, J. Y. Conservation of migration and differentiation circuits in primordial germ cells between avian species. *J. Reprod. Dev.* **59**, 252–257 (2013).
16. Woodcock, M. E. *et al.* Reviving rare chicken breeds using genetically engineered sterility in surrogate host birds. *Proc. Natl. Acad. Sci. U.S.A.* **116**, 20930–20937. https://doi.org/10.1073/pnas.1906316116 (2019).

17. Elbrecht, A. & Smith, R. G. Aromatase enzyme activity and sex determination in chickens. *Science*. **255**, 467–470 (1992).

18. Vaillant, S. *et al.* Degree of sex reversal as related to plasma steroid levels in genetic female chickens (Gallus domesticus) treated with Fadrozole. *Mol. Reprod. Dev.* **65**, 420–428. https://doi.org/10.1002/mrd.10318 (2003).

19. Liu, C. *et al.* Uniparental chicken offsprings derived from oogenesis of chicken primordial germ cells (ZZ). *Biol. Reprod.* **96**, 686–693. https://doi.org/10.1095/biolreprod.116.144253 (2017).

20. Burrows, W. H. & Quinn, J. P. The collection of spermatozoa from the domestic fowl and turkey. *Poult. Sci.* **16**, 19–24 (1937).

21. Govoroun, M. S. *et al.* Isolation of chicken homolog of the FOXL2 gene and comparison of its expression patterns with those of aromatase during ovarian development. *Develop. Dyn. Official Publ. Am. Assoc. Anatom.* **231**, 859–870. https://doi.org/10.1002/dvdy.20189 (2004).

22. Clinton, M., Haines, L., Belloir, B. & McBride, D. Sexing chick embryos: A rapid and simple protocol. *Br. Poult. Sci.* **42**, 134–138. https://doi.org/10.1080/713655025 (2001).

23. Estoup, A., Perrot, E. & Chourrout, D. Rapid one-tube DNA extraction for reliable detection of fish polymorphic markers and transgenes. *Mol. Mar. Biol. Biotechnol.* **5**, 295–298 (1996).

24. Brionne, A., Juanchich, A. & Hennequet-Antier, C. ViSEAGO: A Bioconductor package for clustering biological functions using Gene Ontology and semantic similarity. *BioData mining* **12**, 16. https://doi.org/10.1186/s13040-019-0204-1 (2019).

25. Gene [Internet. JBethesda (MD): National Library of Medicine (US), National Center for Biotechnology Information, <//www.ncbi.nlm.nih.gov/gene/>]

26. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. *J. Mol. Biol.* **215**, 403–410. https://doi.org/10.1016/s0022-2836(05)80360-2 (1990).

27. Yates, A. D. *et al.* Ensembl 2020. *Nucleic Acids Res.* **48**, D682–D688. https://doi.org/10.1093/nar/gkz966 (2020).

28. Coordinators, N. R. Database resources of the National Center for biotechnology information. *Nucleic Acids Res.* **44**, D7-D19. https://doi.org/10.1093/nar/gkv1290 (2016).

29. Rauci, F., Fuet, A. & Pain, B. In vitro generation and characterization of chicken long-term germ cells from different embryonic origins. *Theriogenology* **84**, 732–742. https://doi.org/10.1016/j.theriogenology.2015.04.032 (2015).

30. Ellegren, H. *et al.* Faced with inequality: Chicken do not have a general dosage compensation of sex-linked genes. *BMC Biol.* **5**, 40. https://doi.org/10.1186/1741-7007-5-40 (2007).

31. Uebbing, S. *et al.* Quantitative mass spectrometry reveals partial translational regulation for dosage compensation in Chicken. *Mol. Biol. Evol.* **32**, 2716–2725. https://doi.org/10.1093/molbev/msv147 (2015).

32. Ayers, K. L. *et al.* RNA sequencing reveals sexually dimorphic gene expression before gonadal differentiation in chicken and allows comprehensive annotation of the W-chromosome. *Genome Biol.* **14**, R26. https://doi.org/10.1186/gb-2013-14-3-r26 (2013).

33. Iwahashi, K. *et al.* Autonomous regulation of sex-specific developmental programming in mouse fetal germ cells. *Biol. Reprod.* **77**, 697–706. https://doi.org/10.1095/biolreprod.107.062851 (2007).

34. Zhao, D. *et al.* Somatic sex identity is cell autonomous in the chicken. *Nature* **464**, 237–242. https://doi.org/10.1038/nature08852 (2010).

35. Guyot, N. *et al.* Proteomic analysis of egg white heparin-binding proteins: towards the identification of natural antibacterial molecules. *Sci. Rep.* **6**, 27974. https://doi.org/10.1038/srep27974 (2016).

36. Chen, H. H. *et al.* DAZL limits pluripotency, differentiation, and apoptosis in developing primordial germ cells. *Stem Cell Rep.* **3**, 892–904. https://doi.org/10.1016/j.stemcr.2014.09.003 (2014).

37. Jonson, L. *et al.* Molecular composition of IMP1 ribonucleoprotein granules. *Mol. Cell. Proteom.* **6**, 798–811. https://doi.org/10.1074/mcp.M600346-MCP200 (2007).

38. West, J. A. *et al.* A role for Lin28 in primordial germ-cell development and germ-cell malignancy. *Nature* **460**, 909–913. https://doi.org/10.1038/nature08210 (2009).

39. Childs, A. J., Kinnell, H. L., He, J. & Anderson, R. A. LIN28 is selectively expressed by primordial and pre-meiotic germ cells in the human fetal ovary. *Stem Cells Develop.* **21**, 2343–2349. https://doi.org/10.1089/scd.2011.0730 (2012).

40. Yang, M. *et al.* Lin28 promotes the proliferative capacity of neural progenitor cells in brain development. *Development (Cambridge, England)* **142**, 1616–1627. https://doi.org/10.1242/dev.120543 (2015).

41. Nishino, J., Kim, S., Zhu, Y., Zhu, H. & Morrison, S. J. A network of heterochronic genes including Imp1 regulates temporal changes in stem cell properties. *Elife* **2**, 00924. https://doi.org/10.7554/eLife.00924 (2013).

42. Cimadamore, F., Amador-Arjona, A., Chen, C., Huang, C. T. & Terskikh, A. V. SOX2-LIN28/let-7 pathway regulates proliferation and neurogenesis in neural precursors. *Proc. Natl. Acad. Sci. U.S.A.* **110**, E3017-E3026. https://doi.org/10.1073/pnas.1220176110 (2013).

43. Hsieh, Y. T., Chou, M. M., Chen, H. C. & Tseng, J. J. IMP1 promotes choriocarcinoma cell migration and invasion through the novel effectors RSK2 and PPME1. *Gynecol. Oncol.* **131**, 182–190. https://doi.org/10.1016/j.ygyno.2013.07.106 (2013).

44. Rizzino, A. Sox2 and Oct-3/4: a versatile pair of master regulators that orchestrate the self-renewal and pluripotency of embryonic stem cells. *Wiley Interdisciplin. Rev. Syst. Biol. Med.* **1**(228–236), 2009. https://doi.org/10.1002/wsbm.12 (2009).

45. Izzo, A. *et al.* Dynamic changes in H1 subtype composition during epigenetic reprogramming. *J. Cell Biol.* https://doi.org/10.1083/jcb.201611012 (2017).

46. Besche, H. C., Haas, W., Gygi, S. P. & Goldberg, A. L. Isolation of mammalian 26S proteasomes and p97/VCP complexes using the ubiquitin-like domain from HHR23B reveals novel proteasome-associated proteins. *Biochemistry* **48**, 2538–2549. https://doi.org/10.1021/bi02198q (2009).

47. Andersen, K. M. *et al.* Thioredoxin Txnl1/TPR32 is a redox-active cofactor of the 26 S proteasome. *J. Biol. Chem.* **284**, 15246–15254. https://doi.org/10.1074/jbc.M900016200 (2009).

48. Peth, A., Besche, H. C. & Goldberg, A. L. Ubiquitinated proteins activate the proteasome by binding to Usp14/Ubp6, which causes 20S gate opening. *Mol. Cell* **36**, 794–804. https://doi.org/10.1016/j.molcel.2009.11.015 (2009).

49. Fok, K. L. *et al.* Huwel regulates the establishment and maintenance of spermatogonia by suppressing DNA damage response. *Endocrinology* **158**, 4000–4016. https://doi.org/10.1210/en.2017-00396 (2017).

50. Hunt, L. C. *et al.* A key role for the ubiquitin ligase UBR4 in myofiber hypertrophy in drosophila and mice. *Cell Rep.* **28**, 1268–1281. e1266. https://doi.org/10.1016/j.celrep.2019.06.094 (2019).

51. Kitange, G. J. *et al.* Retinoblastoma binding protein 4 modulates temozolomide sensitivity in glioblastoma by regulating DNA repair proteins. *Cell Rep.* **14**, 2587–2598. https://doi.org/10.1016/j.celrep.2016.02.045 (2016).

52. Liu, F. *et al.* CRL4B(RBBP7) targets HUWE1 for ubiquitination and proteasomal degradation. *Biochem. Biophys. Res. Commun.* **501**, 440–447. https://doi.org/10.1016/j.bbrc.2018.05.008 (2018).

53. Yang, X. *et al.* Histone acetyltransferase 1 promotes homologous recombination in DNA repair by facilitating histone turnover. *J. Biol. Chem.* **288**, 18271–18282. https://doi.org/10.1074/jbc.M113.473199 (2013).

54. Crimmins, S. *et al.* Transgenic rescue of ataxia mice reveals a male-specific sterility defect. *Dev. Biol.* **325**, 33–42. https://doi.org/10.1016/j.ydbio.2008.09.021 (2009).

55. Kovacs, L., Nagy, A., Pal, M. & Deak, P. Usp14 is required for spermatogenesis and ubiquitin stress responses in *Drosophila melanogaster*. *J. Cell Sci.* **1**, 33. https://doi.org/10.1242/jcs.237511 (2020).
56. Castrillon, D. H. et al. Toward a molecular genetic analysis of spermatogenesis in *Drosophila melanogaster*: Characterization of male-sterile mutants generated by single P element mutagenesis. *Genetics* **135**, 489–505 (1993).

57. Chaytow, H., Huang, Y. T., Gillingwater, T. H. & Faller, K. M. E. The role of survival motor neuron protein (SMN) in protein homeostasis. *Cell. Mol. Life Sci. CMLS* **75**, 3877–3894. https://doi.org/10.1007/s00018-018-2849-1 (2018).

58. Workman, E., Kalda, C., Patel, A. & Battle, D. J. Gemin5 binds to the survival motor neuron mRNA to regulate SMN expression. *J. Biol. Chem.* **290**, 15662–15669. https://doi.org/10.1074/jbc.M115.646257 (2015).

59. Fuller, H. R. et al. Spinal muscular atrophy patient iPSC-derived motor neurons have reduced expression of proteins important in neuronal development. *Front. Cell. Neurosci.* **9**, 506. https://doi.org/10.3389/fncel.2015.00506 (2015).

60. Han, K. J. et al. Monoubiquitination of survival motor neuron regulates its cellular localization and Cajal body integrity. *Hum. Mol. Genet.* **25**, 1392–1405. https://doi.org/10.1093/hmg/ddw021 (2016).

61. Kwon, D. Y. et al. The E3 ubiquitin ligase mind bomb 1 ubiquitinates and promotes the degradation of survival of motor neuron protein. *Mol. Biol. Cell* **24**, 1863–1871. https://doi.org/10.1091/mbc.E13-01-0042 (2013).

62. Li, D. K., Tisdale, S., Lotti, F. & Pellizzoni, L. SMN control of RNP assembly: From post-transcriptional gene regulation to motor neuron disease. *Semin. Cell Dev. Biol.* **32**, 22–29. https://doi.org/10.1016/j.semcdb.2014.04.026 (2014).

63. Chang, W. F. et al. Survival motor neuron protein participates in mouse germ cell development and spermatogonium maintenance. *Int. J. Mol. Sci.* **2**, 1. https://doi.org/10.3390/ijms21030794 (2020).

64. Ottesen, E. W. et al. Severe impairment of male reproductive organ development in a low SMN expressing mouse model of spinal muscular atrophy. *Sci. Rep.* **6**, 20193. https://doi.org/10.1038/srep20193 (2016).

65. Grice, S. J. & Liu, J. L. Survival motor neuron protein regulates stem cell division, proliferation, and differentiation in Drosophila. *PLoS Genet.* **7**, e1002030. https://doi.org/10.1371/journal.pgen.1002030 (2011).

66. Lee, L., Davies, S. E. & Liu, J. L. The spinal muscular atrophy protein SMN affects *Drosophila germline* nuclear organization through the U body-P body pathway. *Dev. Biol.* **332**, 142–155. https://doi.org/10.1016/j.ydbio.2009.05.553 (2009).

67. Ito, K., Takahashi, A., Morita, M., Suzuki, T. & Yamamoto, T. The role of the CNOT1 subunit of the CCR4-NOT complex in mRNA deadenylation and cell viability. *Protein Cell* **2**, 755–763. https://doi.org/10.1007/s13238-011-1092-4 (2011).

68. Nagai, M. et al. Mice lacking Ran binding protein 1 are viable and show male infertility. *FEBS Lett.* **585**, 791–796. https://doi.org/10.1016/j.febslet.2011.02.002 (2011).

69. Suzuki, A., Saba, R., Miyoshi, K., Morita, Y. & Saga, Y. Interaction between NANOS2 and the CCR4-NOT deadenylation complex is essential for male germ cell development in mouse. *PLoS ONE* **7**, e33558. https://doi.org/10.1371/journal.pone.0033558 (2012).

70. Caraballo-Miralles, V. et al. Notch signaling pathway is activated in motoneurons of spinal muscular atrophy. *Int. J. Mol. Sci.* **14**, 11424–11437. https://doi.org/10.3390/ijms140611424 (2013).

71. Shen, C. et al. Strawberry Notch 1 (SBNO1) promotes proliferation of spermatogonial stem cells via the noncanonical Wnt pathway in mice. *Asian J. Androl.* https://doi.org/10.4103/aja.aja_65_18 (2018).

72. Watanabe, Y. et al. Notch and Hippo signaling converge on Strawberry Notch 1 (Sbno1) to synergistically activate Cdx2 during specification of the trophectoderm. *Sci. Rep.* **7**, 46135. https://doi.org/10.1038/srep46135 (2017).

73. Xu, J. et al. AIP1-mediated actin disassembly is required for postnatal germ cell migration and spermatogonial stem cell niche establishment. *Cell Death Dis.* **6**, e1818. https://doi.org/10.1038/cddis.2015.182 (2015).

74. Rossoll, W. et al. Smn, the spinal muscular atrophy-determining gene product, modulates axon growth and localization of beta-actin mRNA in growth cones of motoneurons. *J. Cell Biol.* **163**, 801–812. https://doi.org/10.1083/jcb.200304128 (2003).

75. Hensel, N. & Claus, P. The actin cytoskeleton in SMA and ALS: How does it contribute to motoneuron degeneration?. *Neurosci. Rev. J. Bring. Neurobiol. Neurol. Psychiat.* **24**, 54–72. https://doi.org/10.1177/1073858417705059 (2018).

76. Gough, R. E. & Goult, B. T. The tale of two talins—two isoforms to fine-tune integrin signalling. *FEBS Lett.* **592**, 2108–2125. https://doi.org/10.1002/1873-3468.13081 (2018).

77. Mathew, S. et al. Talin regulates integrin beta1-dependent and-independent cell functions in uretic bud development. *Development* (Cambridge, England) **144**, 4148–4158. https://doi.org/10.1242/dev.149914 (2017).

78. Becam, I. E., Tanentzapf, G., Lepesant, J. A., Brown, N. H. & Huynh, J. R. Integrin-independent repression of cadherin transcription by talin during axis formation in Drosophila. *Nat. Cell Biol.* **7**, 510–516. https://doi.org/10.1038/ncb1253 (2005).

79. Dasgupta, S. K. et al. Wdr1-dependent actin reorganization in platelet activation. *PLoS ONE* **11**, e0162897. https://doi.org/10.1371/journal.pone.0162897 (2016).

80. Huss, D. J. et al. Avian primordial germ cells contribute to and interact with the extracellular matrix during early migration. *Front. Cell Develop. Biol.* **7**, 35 (2019).

81. Nakatani, Y. et al. p600, a unique protein required for membrane morphogenesis and cell survival. *Proc. Natl. Acad. Sci. U.S.A.* **102**, 15093–15098. https://doi.org/10.1073/pnas.0507458102 (2005).

82. Belzil, C. et al. p600 regulates spindle orientation in apical neural progenitors and contributes to neurogenesis in the developing neocortex. *Biol. Open* **3**, 475–485. https://doi.org/10.1242/bio.20147807 (2014).

83. Vallee, R. B., McKenney, R. J. & Ori-McKenney, K. M. Multiple modes of cytoplasmic dynein regulation. *Nat. Cell Biol.* **14**, 224–230. https://doi.org/10.1038/ncb2420 (2012).

84. Tong, J., Taylor, P. & Moran, M. F. Proteomic analysis of the epidermal growth factor receptor (EGFR) interactome and post-translational modifications associated with receptor endocytosis in response to EGF and stress. *Mol. Cell. Proteom. MCP* **13**, 1644–1658. https://doi.org/10.1074/mcp.M114.038596 (2014).

85. Singh, R. N., Howell, M. D., Ottesen, E. W. & Singh, N. N. Diverse role of survival motor neuron protein. *Biochim. Biophys. Acta* **1860**, 299–315. https://doi.org/10.1016/j.bbagrm.2016.12.008 (2017).

86. Le Roy, C. & Wrana, J. L. Clathrin- and non-clathrin-mediated endocytic regulation of cell signalling. *Nat. Rev. Mol. Cell Biol.* **6**, 112–126. https://doi.org/10.1038/nrm1571 (2005).

87. Rieger, L. & O'Connor, R. Controlled signaling-insulin-like growth factor receptor endocytosis and presence at intracellular compartments. *Front. Endocrinol.* **11**, 620013. https://doi.org/10.3389/fendo.2020.620013 (2020).

88. Wiedlocha, A. & Sørensen, V. Signaling, internalization, and intracellular activity of fibroblast growth factor. *Curr. Top. Microbiol. Immunol.* **286**, 45–79. https://doi.org/10.1007/978-3-540-69494-6_3 (2004).

89. Alborzinia, H. et al. Quantitative kinetics analysis of BMP2 uptake into cells and its modulation by BMP antagonists. *J. Cell Sci.* **126**, 117–127. https://doi.org/10.1242/jcs.109777 (2013).

90. Chen, Y., Huang, L., Qi, X. & Chen, C. Insulin receptor trafficking: Consequences for insulin sensitivity and diabetes. *Int. J. Mol. Sci.* https://doi.org/10.3390/ijms20205007 (2019).

91. Hegarty, S. V., Sullivan, A. M. & O'Keeffe, G. W. Endocytosis contributes to BMP2-induced Smad signalling and neuronal growth. *Neurosci. Lett.* **643**, 32–37. https://doi.org/10.1016/j.neulet.2017.02.013 (2017).

92. Tuzon, C. T., Rigueur, D. & Merrill, A. E. Nuclear fibroblast growth factor receptor signaling in skeletal development and disease. *Curr. Osteoporos. Rep.* **17**, 138–146. https://doi.org/10.1007/s11914-019-00512-2 (2019).

## Acknowledgements

We are very grateful to the staff of the experimental unit PEAT for egg production and animal care, Grégoire Harichaux for technical assistance in proteomic experiments, Lydie Oliveira-Correia and Plateforme d'Analyse Protéomique de Paris Sud-Ouest (PAPPSO), Quintin Viboud for technical assistance in RT-qPCR analysis, Dr
Bertrand Pain for DDX4 antibody and DDX4 primer sequences, Dr Joël Gautron and Dr Sophie Réhault-Godbert for primer sequences for VMO1, OVAL, OVALX, BPIF2, SPINK5 genes and for magnum RNA.

## Author contributions

M.S.G. conceived and designed the experiments. M.J.M. brought his expertise in PGCs culture. S.A., V.G., M.C.F. and M.S.G. performed PGCs culture. V.L. designed and supervised proteomic analysis. L.S. performed proteomic analysis. L.C.S. performed protein identification and quantification. A.B. performed bioinformatic analysis of the data. A.T. contributed to the analysis of the data. S.A. performed Western Blot analysis. S.A., A.J., M.C.F. performed RT-QPCR and immunostaining of PGCs. S.A., M.S.G. performed gonad colonization tests. M.S.G. analyzed the data. M.S.G., L.S. wrote the paper. V.L. wrote the Sects. *In-gel Protein Digestion and NanoLC-MS/MS for PGC protein identification and Label-free protein quantification* of the Methods. M.J.M., S.F., E.B. critically reviewed and edited the manuscript. All authors approved the manuscript.

## Funding

This study was supported by grants VALBIODI funded by Région Centre Val-de-Loire, the French National Infrastructure of Research CRB Anim funded by "Investissements d'avenir," ANR-11-INBS-0003, from the European Union's Horizon 2020 Research and Innovation Program under Grant Agreement N°677353. The high resolution mass spectrometer (Project SMHART, 35069) was financed by European Regional Development Fund (ERDF), the Val-de-Loire Region general comity, the French National Research Institute for Agriculture, Food and Environment (INRAE) and the French National Institute of Health and Medical Research (INSERM). Laura Soler has received the support of the EU in the framework of the Marie-Curie FP7 COFUND People Programme, through the award of an AgreenSkills fellowship (under Grant Agreement No. 267196).

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-98454-2.

Correspondence and requests for materials should be addressed to M.S.G.

Reprints and permissions information is available at www.nature.com/reprints.

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

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021