Intersexuality is a congenital reproductive disorder that usually occurs in hornless goats, hindering breeding of goats with hornless traits and the development of the goat industry. In this study, we aimed to identify differentially expressed genes in intersex and normal goat gonads by comparing gene transcription profiles of intersex and normal goat gonads. As intersex goats are genetically based on females, we chose female goats as controls. The goats in the control group and the experimental group were both over one‐year old. We evaluated the anatomical characteristics of the reproductive organs of five intersex goats using histopathological methods. The gonads were found to be ovarian and testicular types. RNA‐Seq technology was used to identify differentially expressed genes in gonads and normal goat ovary tissues. Transcription analysis results were verified by qPCR. The results showed that 2,748 DEGs were upregulated and 3,327 DEGs were downregulated in intersex ovaries unlike in controls, whereas 2006 DEGs were upregulated and 2032 DEGs were downregulated in the interstitial testes. Many of these genes play important roles in mammalian sex determination and sex differentiation, such as SOX9, WT1, GATA4, DMRT1, DHH, AMH, CYP19A1 and FST. We found that many DEGs are involved in biological developmental regulation by GO and KEGG enrichment analyses, and that most genes associated with the steroid synthesis pathway were downregulated. The DEGs identified in this study may be involved in the regulation of intersex goat sex determination and differentiation, and may increase our understanding of the molecular mechanisms of mammalian sex differentiation.

The conditions for sex reversal in vertebrate species have been extensively studied, and the results highlighted numerous key factors involved in sex differentiation. However, the transcriptomes in hypothalamic and pituitary tissues from intersex goats have rarely been studied. The aim of this study was to screen candidate genes and signalling pathways related to sex reversal in Huai goats by analyzing gene expression in hypothalamic and pituitary tissues via transcriptome sequencing and bioinformatics analyses.In total, 612 and 139 differentially expressed genes (DEGs) were identified between the intersex and non-intersex groups in the hypothalamus and pituitary, respectively.The DEGs in the hypothalamus and pituitary were significantly enriched in 41 and 16 signalling pathways, respectively, including the calcium signalling pathway, neuroactive ligand-receptor interaction signalling pathway, and oestrogen signalling pathway, which might be related to intersex sex development disorders. A candidate gene from the tachykinin family (TACR1) was significantly enriched in the calcium signalling pathway. Thirty-one DEGs were shared between these two comparisons and were enriched in several acetyl-CoA-related processes and the oestrogen signalling pathway. The results of the real-time PCR analysis show that the transcriptome sequencing results were reliable. The transcriptome data indicate that the regulation of various physiological systems is involved in intersex goat development. Therefore, these results provide helpful data enhancing our understanding of the molecular mechanisms underlying intersex syndrome in goats.