Sex reversal has been studied extensively in vertebrate species, particularly in domestic goats, because polled intersex syndrome (PIS) has seriously affected their production efficiency. In the present study, we used histopathologically diagnosed cases of PIS to identify correlated genomic regions and variants using representative selection signatures and performed GWAS using Restriction-Site Associated Resequencing DNA. We identified 171 single-nucleotide polymorphisms (SNPs) that may have contributed to this phenotype, and 53 SNPs were determined to be located in coding regions using a general linear model. The transcriptome data sets of differentially expressed genes (DEGs) in the pituitary tissues of intersexual and nonintersexual goats were examined using high-throughput technology. A total of 10,063 DEGs and 337 long noncoding RNAs were identified. The DEGs were clustered into 56 GO categories and determined to be significantly enriched in 53 signaling pathways by KEGG analysis. In addition, according to qPCR results, PSPO2 and FSH were significantly more highly expressed in sexually mature pituitary tissues of intersexual goats compared to healthy controls (nonintersexual). These results demonstrate that certain novel potential genomic regions may be responsible for intersexual goats, and the transcriptome data indicate that the regulation of various physiological systems is involved in intersexual goat development. Therefore, these results provide helpful data for understanding the molecular mechanisms of intersex syndrome in goats.
Keywords: GWAS; Intersexual goat; RAD; RNAseq.