Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae.

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplications (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Characterization of the sex determining region of channel catfish (Ictalurus punctatus) and development of a sex-genotyping test.

Channel catfish is an important species for aquaculture that exhibits a sexually dimorphic growth in favor of males. Genetic sexing and development of sex markers are crucial for the early identification of sex and of particular genotypes (YY males) for the production of all-male population in channel catfish aquaculture. In this study, we sequenced genomic DNA from pools of males and pools of females to better characterize the sex determining region (SDR) of channel catfish and to develop sex-specific markers for genetic sexing. Performing comparative analyses on male and female pooled genomic reads, we identified a large SDR (∼8.3 Mb) in the middle of channel catfish linkage group 4 (LG04). This non-recombining SDR contains a high-density of male-specific (Y chromosome) fixed single nucleotide polymorphisms (SNPs) along with âˆ¼ 185 kb male-specific insertions or deletions. This SDR contains 95 annotated protein-encoding genes, including the recently reported putative channel catfish master sex determining (MSD) gene, breast cancer anti-estrogen resistance protein 1 (bcar1), located at one edge of the SDR. No sex-specific SNPs and/or indels were found in the coding sequence of bcar1, but one male-specific SNP was identified in its first intron. Based on this genomic information, we developed a PCR-based sex-specific genetic test. Genotyping results confirmed strong linkage between phenotypic sexes and the identified SDR in channel catfish. Our results confirm, using a Pool-Seq approach, that channel catfish is male heterogametic (XX-XY) with a large SDR on the LG04 sex chromosome. Furthermore, our genotyping primers can be used to identify XX, XY, and YY fish that will facilitate future research on sex determination and aquaculture applications in channel catfish.

Evidence for selection and spatially distinct patterns found in a putative zona pellucida gene in Pacific cod, and implications for management.

Genetic differentiation has been observed in marine species even when no obvious barriers to gene flow exist, and understanding such differentiation is essential for effective fisheries management. Highly differentiated outlier loci can provide information on how genetic variation might not only contribute to local adaptation but may also be affected by historical demographic events. A locus which aligned to a predicted zona pellucida sperm-binding protein 3 gene (ZP3) in Atlantic cod (Gadus morhua) was previously identified as the highest outlier based on F ST in a RADseq study of Pacific cod (Gadus macrocephalus) across the West Coast of North America. However, because of the limited length of the RAD sequence and restricted geographic area of sampling, no conclusion on the functional significance of the observed variation was possible. In other marine species, ZP3 is involved in reproductive isolation, local adaptation, and has neofunctionalized as an antifreeze gene, and so it may provide important insights in functional population structure of Pacific cod. Here, we sequenced a 544-bp region of ZP3 in 230 Pacific cod collected from throughout their geographic range. We observed striking patterns of spatial structuring of ZP3 haplotypes, with a sharp break near Kodiak, Alaska, USA where populations within ~200 km of each other are nearly fixed for different haplotypes, contrasting a pattern of isolation by distance at other genetic markers in this region (F ST = 0.003). Phylogenetic analysis of ZP3 haplotypes revealed that the more southern haplotypes appear to be ancestral, with the northern haplotype evolving more recently, potentially in response to a novel selective pressure as Pacific cod recolonized northern latitudes after glaciation. The sharp break in haplotype frequencies suggests strong selective pressures are operating on small spatial scales and illustrates that selection can create high divergence even in marine species with ample opportunities for gene flow.

Prion protein polymorphisms in Michigan white-tailed deer (Odocoileus virginianus).

Chronic Wasting Disease (CWD), a well-described transmissible spongiform encephalopathy of the Cervidae family, is associated with the aggregation of an abnormal isoform (PrPCWD) of the naturally occurring host prion protein (PrPC). Variations in the PrP gene (PRNP) have been associated with CWD rate of infection and disease progression. We analysed 568 free-ranging white-tailed deer (Odocoileus virginianus) from 9 CWD-positive Michigan counties for PRNP polymorphisms. Sampling included 185 CWD-positive, 332 CWD non-detected, and an additional 51 CWD non-detected paired to CWD-positives by sex, age, and harvest location. We found 12 polymorphic sites of which 5 were non-synonymous and resulted in a change in amino acid composition. Thirteen haplotypes were predicted, of which 11 have previously been described. Using logistic regression, consistent with other studies, we found haplotypes C (OR = 0.488, 95% CI = 0.321-0.730, P < 0.001) and F (OR = 0.122, 95% CI = 0.007-0.612, P < 0.05) and diplotype BC (OR = 0.340, 95% CI = 0.154-0.709, P < 0.01) were less likely to be found in deer infected with CWD. As has also been documented in other studies, the presence of a serine at amino acid 96 was less likely to be found in deer infected with CWD (P < 0.001, OR = 0.360 and 95% CI = 0.227-0.556). Identification of PRNP polymorphisms associated with reduced vulnerability to CWD in Michigan deer and their spatial distribution can help managers design surveillance programmesand identify and prioritize areas for CWD management.