Combined Proteomic and Metabolomic Analysis Reveals Comprehensive Regulation of Somatostatin DNA Vaccine in Goats.

Somatostatin (SS) plays crucial regulatory roles in animal growth and reproduction by affecting the synthesis and secretion of growth hormone (GH). However, the mechanism by which SS regulates growth and development in goats is still unclear. In order to investigate the regulatory networks of the hypothalamus and pituitary in goats affected by SS DNA vaccines, in this study, we used a previously established oral attenuated Salmonella typhimurium SS DNA vaccine, X9241 (ptCS/2SS-asd), to treat wethers. We analyzed the protein changes in hypothalamic and pituitary tissues using a TMT-based proteomics approach. Additionally, we examined the metabolic profiles of the serum of control and immunized wethers through untargeted metabolomics using liquid chromatography-mass spectrometry (LC-MS). Key signaling pathways were identified based on differentially expressed metabolites (DEMs) and differentially expressed proteins (DEPs). Furthermore, the effect of critical DEPs on signaling pathways was confirmed through Western blotting (WB) experiments, which elucidated the mechanism of active SS immunization in wethers. A proteomics analysis revealed that the expression of 58 proteins in the hypothalamus and 124 in the pituitary gland was significantly altered following SS vaccine treatment (fold change > 1.2 or < 0.83, p < 0.05). In the hypothalamus, many DEPs were associated with gene ontology (GO) terms related to neuronal signaling. In contrast, most DEPs were associated with metabolic pathways. In the pituitary gland, the DEPs were largely related to immune and nutrient metabolism functions, with significant enrichment in KEGG pathways, particularly those involving the metabolic pathway, sphingolipid signaling, and the cGMP-PKG signaling pathway. A metabolomic analysis further showed that active SS immunization in wethers led to significant alterations in seven serum metabolites. Notably, the sphingolipid signaling pathway, secondary bile acid synthesis, sphingolipid metabolism, and lysine synthesis were significantly disrupted. SS vaccines induced marked changes in hypothalamic-pituitary proteins in wethers, facilitating alterations in their growth processes. This study not only provides insights into the mechanism of the SS gene in regulating GH secretion in wethers but also establishes a basis for hormone immunoregulation technology to enhance livestock production performance.

Immunisation of the somatostatin gene alters hypothalamic-pituitary-liver gene expressions and enhances growth in Dazu black goats.

Objective: Somatostatin (SS) plays important regulatory roles in animal growth and reproduction by affecting the synthesis and secretion of growth hormone (GH). However, the mechanism by which SS regulates growth and development in goats is still unclear. Methods: In this study, we randomly selected eight 7-month-old Dazu black goats (DBGs) of similar body weight and equally assigned four bucks as the immunised and negative control groups. The immunised group received the Salmonella typhi attenuated vaccine CSO22 (ptCS/2SS-asd) orally, whilst the negative control group received the empty vector vaccine CSO22 (pVAX-asd) orally. Results: The SS concentration in the serum of goats in the immunised group was significantly lower than that in the negative control group, and the daily gain was significantly higher (p < 0.05). SS-14 DNA vaccine immunisation resulted in significantly higher concentrations of growth-related hormones such as GH-releasing hormone and IGF-1 in the serum of goats (p < 0.05). RNA-seq analysis of hypothalamus of oral SS-14 DNA vaccine and negative control DBGs identified 31 differentially expressed genes (DEGs). Pituitary gland identified 164 DEGs. A total of 246 DEGs were detected in the liver by RNA-seq. Gene ontology (GO) of DEGs was enriched in mitochondrial envelope, extracellular region, receptor binding and cell proliferation. The biological metabolic pathways associated with DEGs were explored by Kyoto Encyclopedia of Genes and Genomes analysis. DEGs were associated with metabolic pathways, oxidative phosphorylation, vitamin digestion and absorption and galactose metabolism. These candidate genes (e.g. DGKK, CYTB, DUSP1 and LRAT) may provide references for exploring the molecular mechanisms by which SS promotes growth and development. Conclusion: Overall, these results demonstrated that the SS DNA vaccine enhanced the growth of DBGs by altering growth-related hormone concentrations and regulating the expression of growth-related genes in the hypothalamic-pituitary-liver axis.

Novel Heredity Basis of the Four-Horn Phenotype in Sheep Using Genome-Wide Sequence Data.

Horns are an important breeding trait for sheep. However, no widely recognized viewpoint on the regulatory genes and mechanisms of horns is available, and the genetic basis of the four-horn phenotype (FHP) is unclear. This work conducted a genome-wide association study with 100 sheep genomes from multiple breeds to investigate the genetic basis of the FHP. The results revealed three significant associations (corrected as p < 1.64 × 10-8) of the InDels (CHR2: g.133,742,709delA, g.133,743,215insC, and g.133,743,940delT) for FHP in the intergenic sequence (IGS) between the MTX2 and the LOC105609047 of CHR2. Moreover, 14 significant associations (corrected as p < 1.42 × 10-9) of SNPs with the FHP phenotype were identified in CHR2 and CHR16, including five (e.g., CHR16: g.40,351,378G > A and g.40,352,577G > A) located in the intron of the ADAMTS12 gene, eight (e.g., CHR2: g.133,727,513C > T and g.133,732,145T > G) in the IGS between MTX2 and LOC105609047, and only one (CHR2: g.133,930,761A > G) in the IGS between HOXD1 and MTX2. Obvious divergence was also observed in genotype patterns between the FHP and others (two horns and hornless) in the HOXD1 and ADAMTS12 gene regions. An extremely significant linkage also occurred between Loci I and Loci II within 100 individuals (LD = -156.02186, p < 0.00001). In summary, our study indicated that the genomic sequences from CHR2 and CHR16 contributed to the FHP in sheep, specifically the key candidate genes HOXD1 and ADAMTS12. These results improved our understanding of the Mendelian genetic basis of the FHP in sheep.

A 1.1 Mb duplication CNV on chromosome 17 contributes to skeletal muscle development in Boer goats.

The Boer goat is one of the top meat breeds in modern animal husbandry and has attracted widespread attention for its unique growth performance. However, the genetic basis of muscle development in the Boer goat remains obscure. In this study, we identified specific structural variants in the Boer goat based on genome-wide selection signals and analyzed the basis of the molecular heredity of related candidate genes in muscle development. A total of 9 959 autosomal copy number variations (CNVs) were identified through selection signal analysis in 127 goat genomes. Specifically, we confirmed that the highest signal CNV (HSV) was a chromosomal arrangement containing an approximately 1.11 Mb (CHIR17: 60062304-61171840 bp) duplicated fragment inserted in reverse orientation and a 5 362 bp deleted region (CHIR17:60145940-60151302 bp) with overlapping genes (e.g., ARHGAP10, NR3C2, EDNRA, PRMT9, and TMEM184C). The homozygous duplicated HSV genotype (+/+) was found in 96% of Boer goats but was not detected in Eurasian goats and was only detected in 4% of indigenous African goats. The expression network of three candidate genes ( ARHGAP10, NR3C2, and EDNRA) regulating dose transcription was constructed by RNA sequencing. Results indicated that these genes were involved in the proliferation and differentiation of skeletal muscle satellite cells (SMSCs) and their overexpression significantly increased the expression of SAA3. The HSV of the Boer goat contributed to superior skeletal muscle growth via the dose effects of overlapping genes.

Comparative analysis of the genetic diversity of the neutral microsatellite loci and second exon of the goat MHC-DQB1 gene.

This study compared and analyzed the genetic diversity and population structure of exon 2 of the DQB1 gene and 13 autosomal neutral microsatellite markers from 14 Chinese goat breeds to explore the potential evolutionary mechanism of the major histocompatibility complex (MHC). A total of 287 haplotypes were constructed from MHC-DQB1 exon 2 from 14 populations, and 82 nucleotide polymorphic sites (SNPs, 31.78%) and 172 heterozygous individuals (79.12%) were identified. The FST values of the microsatellites and MHC-DQB ranged between 0.01831-0.26907 and 0.00892-0.38871, respectively. Furthermore, 14 goat populations showed rich genetic diversity in the microsatellite loci and MHC-DQB1 exon 2. However, the population structure and phylogenetic relationship represented by the two markers were different. Positive selection and Tajima's D test results showed the occurrence of a diversified selection mechanism, which was primarily based on a positive and balancing selection in goat DQB. This study also found that the DQB sequences of bovines exhibited trans-species polymorphism (TSP) among species and families. In brief, this study indicated that positive and balancing selection played a major role in maintaining the genetic diversity of DQB, and TSP of MHC in bovines was common, which enhanced the understanding of the MHC evolution.

Identification of differentially expressed long non-coding RNAs and messenger RNAs involved with muscle development in Dazu black goats through RNA sequencing.

This study aimed to explore the genetic basis of muscle development in goats. The transcriptome dataset for differentially expressed lncRNAs (DELs) and differentially expressed genes (DEGs) of goat muscle at different developmental stages were obtained using RNA-Seq. A total of 447,806,481 and 587,559,465 clean reads in the longissimus dorsi muscle of Dazu black goats between 75d embryonic stage and 1d after birth were generated through Illumina paired-end sequencing, and their mapping rates were 89.82 and 90.99%, respectively. Moreover, 4517 DEGs and 648 DELs were identified, and 4784 lncRNA-mRNA targeting relationships were predicted. Gene function annotation results showed that 4101 DEGs were significantly enriched to 1098 GO terms, and 2014 DEGs were significantly enriched to 40 KEGG pathways, including many GO terms and pathways related to muscle development, such as cell differentiation and Wnt signaling pathway. Then, 10 DELs and 20 DEGs were randomly selected for RT-qPCR verification, and the agreement rate between the verification and RNA-Seq results was 90%, indicating the high reliability of the RNA-Seq data analysis. In conclusion, this study obtained several mRNAs and lncRNAs related to the muscle development of Dazu black goats and identified several targeted regulatory pairs of lncRNA-mRNA. This study may serve as a reference to understand the genetic basis and molecular mechanism of muscle development in goats.

Genome-Wide Selective Analysis of Boer Goat to Investigate the Dynamic Heredity Evolution under Different Stages.

Boer goats, as kemp in meat-type goats, are selected and bred from African indigenous goats under a long period of artificial selection. Their advantages in multiple economic traits, particularly their plump growth, have attracted worldwide attention. The current study displayed the genome-wide selection signature analyses of South African indigenous goat (AF), African Boer (BH), and Australian Boer (AS) to investigate the hereditary basis of artificial selection in different stages. Four methods (principal component analysis, nucleotide diversity, linkage disequilibrium decay, and neighbor-joining tree) implied the genomic diversity changes with different artificial selection intensities in Boer goats. In addition, the θπ, FST, and XP-CLR methods were used to search for the candidate signatures of positive selection in Boer goats. Consequently, 339 (BH vs. AF) and 295 (AS vs. BH) candidate genes were obtained from SNP data. Especially, 10 genes (e.g., BMPR1B, DNER, ITGAL, and KIT) under selection in both groups were identified. Functional annotation analysis revealed that these genes are potentially responsible for reproduction, metabolism, growth, and development. This study used genome-wide sequencing data to identify inheritance by artificial selection. The results of the current study are valuable for future molecular-assisted breeding and genetic improvement of goats.

Genome-Wide Selection Sweep between Wild and Local Pigs from Europe for the Investigation of the Hereditary Characteristics of Domestication in Sus Scrofa.

The phenotypic characteristics of existing domestic pigs (DPs) greatly differ from those of wild boar (WB) populations thousands of years ago. After thousands of years of human domestication, WB and DP have diverged greatly in terms of genetics. Theoretically, worldwide local pigs have independent contributions from their local WBs at the beginning of Sus scrofa domestication. The investigation of the vicissitude of the heredity material between domestic populations and their wild ancestors will help in further understanding the domestication history of domestic animals. In the present study, we performed a genome-wide association scan (GWSA) and phylogeny estimation with a total of 1098 public European Illumina 60K single nucleotide polymorphism data, which included 650 local DPs and 448 WBs. The results revealed that the phylogenetic relationship of WBs corresponds to their geographical distribution and carries large divergence with DPs, and all WB breeds (e.g., HRWB, SBWB, and TIWB) presents a closely linkage with the middle WB (e.g., HRWB, and PLWB). In addition, 64 selected candidate genes (e.g., IDH2, PIP5K1B, SMARCA2, KIF5C, and TJP2) were identified from GWSA. A total of 63 known multiple biological functional pathways were annotated by 22 genes, and ubiquinone and other terpenoid-quinone biosynthesis pathways that belong to the metabolism of cofactors and vitamins were significantly enriched (p < 0.05). The most frequent (28.57%) pathways were classified under metabolism. We confirmed that the middle European WB has made an important genetic contribution to the entire European WB populations. A series of selected genes discovered from this study provides the scientific community with a deeper understanding of the heredity performance of metabolism and emotion and the real purpose behind domestication.

Identification of differentially expressed microRNAs in ovulatory and subordinate follicles in Dazu black goats.

How multiple ovulations happen in prolific goats is still unknown even though studies on ovarian physiology and folliculogenesis have made extensive progress. MicroRNAs (miRNAs) are a class of endogenous regulatory factors and regulate structural gene expression mainly at the post-transcriptional level. In this study, in the follicular phase, six large follicles were collected from prolific Dazu black goat and used to generate RNA libraries for RNA sequencing. Based on the litter size and average number of ovulatory follicles in Dazu black goats, the largest three follicles were sorted as ovulatory follicles, and the remaining as subordinate ones. In total, 418 known miRNAs and 110 novel miRNAs were found, and the expression of six randomly selected miRNAs was validated by quantitative PCR analysis. Nine miRNAs were differently expressed between the ovulatory and subordinate follicles (p < 0.01). Chi-miR-582-5p, novel-130, chi-miR-214-3p, and chi-miR-500-5p were upregulated in the ovulatory group, and chi-miR-383, chi-miR-130b-5p, chi-miR-92a-3p, chi-miR-125b-5p, and novel-9 were downregulated in the same group. Chi-miR-130b-5p and chi-miR-214-3p were predicted to target at LHR (XM_013967581.1), GDF9 (NM_001285708.1), BMP15 (NM_001285588.1), and CYP19A1 (XM_013967046.1). In conclusion, nine miRNAs were differently expressed between ovulatory and subordinate follicles, and chi-miR-130b-5p and chi-miR-214-3p were predicted to regulate the expression of genes involved in gonadotropin hormone signaling and oocyte-derived growth factors. Additional studies are needed to confirm these findings.

A nonsynonymous SNP within the AMH gene is associated with litter size in Dazu black goats.

AMH, KISS1R and GDF9 genes play a vital role in human and animal reproduction and might be used as the genetic markers for the reproduction traits selection. The aim of this study was to screen the single nucleotide polymorphisms (SNPs) within the AMH, KISS1R and GDF9 genes and to determine the correlations between these SNPs and the litter size in goats. Nine single SNPs within these genes were used for genotyping of the 190 Dazu black goat populations by SNaPshot technique. The polymorphisms of nine SNPs within these genes were detected in Dazu black goats. The significant correlation was observed between one SNP (g.89172108A > C) within the AMH gene and the litter size of second born in Dazu black goats (p < 0.05). The SNP was located in exon 4 (XM_018050765.1) of the AMH gene and was one nonsynonymous substitution, which resulted in a change of an amino acid from Glutamine to Proline (Gln38Pro). These results suggested that the nonsynonymous SNP g.89172108A > C of AMH gene could be used as a potential genetic marker for Marker-assisted selection (MAS) in goats breeding programs.

Genetic Basis of Follicle Development in Dazu Black Goat by Whole-Transcriptome Sequencing.

The follicle development (FD) is an important factor determining litter size in animals. Recent studies have found that noncoding RNAs (ncRNAs) play an important role in FD. In particular, the role of the regulatory mechanism of competing endogenous RNAs (ceRNAs) that drive FD has attracted increasing attention. Therefore, this study explored the genetic basis of goat FD by obtaining the complete follicular transcriptome of Dazu black goats at different developmental stages. Results revealed that 128 messenger RNAs (mRNAs), 4 long noncoding RNAs (lncRNAs), 49 microRNAs (miRNAs), and 290 circular RNAs (circRNAs) were significantly differentially expressed (DE) between large and small follicles. Moreover, DEmRNAs were enriched in many signaling pathways related to FD, as well as GO terms related to molecular binding and enzyme activity. Based on the analysis of the ceRNA network (CRN), 34 nodes (1 DElncRNAs, 10 DEcircRNAs, 14 DEmiRNAs, and 9 DEmRNAs) and 35 interactions (17 DEcircRNAs-DEmRNAs, 2 DElncRNAs-DEmiRNAs, and 16 DEmRNA-DEmiRNAs) implied that the CRN could be involved in the FD of goats. In conclusion, we described gene regulation by DERNAs and lncRNA/circRNA-miRNA-mRNA CRNs in the FD of goats. This study provided insights into the genetic basis of FD in precise transcriptional regulation.

Three SNPs within exons of INHA and ACVR2B genes are significantly associated with litter size in Dazu black goats.

The aim of this study was to analyse the association between single-nucleotide polymorphisms within INHA and ACVR2B and litter size in Dazu black goats. In total, twenty-two SNPs were genotyped in 190 individuals by SNaPshot and resequencing. The results showed that three SNPs (SNP_1, SNP_12 and SNP_13 in this study) were detected to have significant additive genetic effect on the recorded goat litter size (p < .05). The SNP_1 (NC_030809.1), a non-synonymous substitution of G for T at chr2-g. 28314990 in the exon 2 of INHA gene (NM_001285606.1), resulted in homozygote 2 (HOM2) contributed 0.25 and heterozygote (HET) contributed 0.12 larger litter than homozygote 1 (HOM1). Meanwhile, SNP_12 (Chr22-g. 11721225 A > T) and SNP_13 (Chr22-g. 11721227 A > C) (NC_030829.1) simultaneously mutated at the first and third position of a triplet AAA (lysine, K) in the exon 4 of ACVR2B gene (XM_018066623.1) had estimated genetic effects of HOM1 (0.00) and HOM2 (0.03) larger than HET (-0.12). In conclusion, one SNPs (chr2-g. 28314990 T > G) within the exon 2 of INHA and two SNPs (Chr22-g. 11721225 A > T and Chr22-g. 11721227 A > C) in the exon 4 of ACVR2B gene were highly recommended as candidate markers of litter size in Dazu black goats. A large-scale association study to assess the impact of these variants on litter size is still necessary.

Immunization with oral KISS1 DNA vaccine inhibits testicular Leydig cell proliferation mainly via the hypothalamic-pituitary-testicular axis and apoptosis-related genes in goats.

This study aimed to analyze the effect and mechanism of immunization of oral KISS1 DNA vaccine on the proliferation of goat testicular Leydig cells. Ten 8-week-old male goats were randomly divided into KISS1 DNA vaccine and control groups for immunization (five goats each group). These goats were sacrificed at 8 weeks after primary immunization, and the tissue samples of hypothalamus, pituitary, and testis and Leydig cell samples were collected for RT-PCR and CCK8 assay. Immunization with the oral KISS1 DNA vaccine effectively inhibited the proliferation of Leydig cells, the expression of hypothalamus KISS1, GPR54, and GnRH mRNA, pituitary GnRHR and LH mRNA, testicular LHR mRNA, and apoptosis-inhibitory gene Bcl-2 mRNA in Leydig cells. By contrast, the immunization enhanced the mRNA expression of apoptosis-promoting gene Bax and Clusterin in Leydig cells. These findings indicate that immunization with the oral KISS1 DNA vaccine can inhibit the proliferation of goat testicular Leydig cells mainly via the hypothalamic-pituitary-testicular axis and apoptosis-related genes.

Comparative and selection sweep analysis of CNV was associated to litter size in Dazu black goats.

This study aims to identify the relative Copy number variation (CNV) associated with the litter size of Dazu black goats based on the unpublished CNV analytical results of our previously published sequencing data, in which the litter-size groups were classified into extreme low- and high-yield groups. Firstly, to compare the existence of valuable CNV in Dazu black goats with different fertility levels with mixed pools. We obtained 4992 and 4888 CNVs from the HY and LY, which overlapping 1461 genes, and classified on the original CNV type. Three genes [LOC108633278, PPP1R12A, and YIPF4] were observed in the intersection between the HY deletion and the LY duplication groups. Secondly, on individuals level, we identified a novel candidate CNV (Chr1_50215501, FST = 0.148, VST = 0.347) from 214 autosomal credible CNVs to be significant with litter size in the Dazu black goat, which located in the CBLB gene. This finding indicates the CBLB gene may affect the litter size of the Dazu black goats through structural variations, and Chr1_50215501 can be an effective genetic marker for marker-assisted selection breeding, and this study was also helps understand the molecular mechanism related to the goat litter size.

SNP screening of the MSTN gene and correlation analysis between genetic polymorphisms and growth traits in Dazu black goat.

The aim of this study was to detect SNPs in myostatin (MSTN) gene of four goat breeds, and analyze the correlation of these markers on body measurement traits in the Dazu black goat breed. In total, twenty polymorphic sites were found in one hundred forty-eight individuals, and all SNPs were distributed in introns 1 and 2, except g. 425 C > T, which was found in the regulatory region. Three SNPs (g. 2732 C > T, g. 2752 G > A and g. 4552 A > C) were polymorphic in all four breeds. None of the tag SNPs (g. 425 C > T, g. 1583 A > G, 2732 C > T, g. 4552 A > C and g. 5167 C > T) were significantly correlated with body measurement traits (p > 0.05) in the Dazu black goat. However, individuals with genotype combination 3 (GtC 3) of tag SNPs had higher birth weight and weaning weight than individuals with the other genotype combinations (p < 0.05). Moreover, the genotype combination 4 (GtC 4) was significantly associated with body length and height at the age of 2 months (p < 0.05), and genotype combination 13 (GtC 13) was significantly correlated with body height at 6 months (p < 0.05). Briefly, the combined tag SNP markers might be useful for goat marker-associated selective breeding.

SNP of AHSA2 gene in three cattle breeds using snapshot technology.

Droughtmaster is a tropical breed of beef cattle that can survive in hot climates and easily adapt to torrid environments. These traits are important in livestock breeding. In this study, we genotyped five single-nucleotide polymorphisms (SNPs) of the AHSA2 gene from 190 cattle belonging to three different breeds (Droughtmaster, Angus and Simmental) by using snapshot technology. This work aimed to identify the valuable molecular marker of heat resistance in cattle. Results showed that Droughtmaster exhibited higher expected heterozygosity and polymorphic information content compared with the two other breeds. The AHSA2-1 locus deviated from the Hardy-Weinberg equilibrium in the Droughtmaster breed (P < 0.05). Two SNPs in Droughtmaster diverged significantly from Angus and Simmental. The SNPs were identified as AHSA2-3 and AHSA2-4, which were closely linked to the three breeds based on pair-wise FST. AHSA2-4 involved a missense mutation. In summary, the GG genotypes in AHSA2-3 and AHSA2-4 may be candidate genotypes associated with heat resistance traits and may serve as valuable genetic markers for breeding of heat-tolerant beef cattle in the future.

Evolutionary relationship and population structure of domestic Bovidae animals based on MHC-linked and neutral autosomal microsatellite markers.

Major histocompatibility complex (MHC) genes are critical for disease resistance or susceptibility responsible for host-pathogen interactions determined mainly by extensive polymorphisms in the MHC genes. Here, we examined the diversity and phylogenetic pattern of MHC haplotypes reconstructed using three MHC-linked microsatellite markers in 55 populations of five Bovidae species and compared them with those based on neutral autosomal microsatellite markers (NAMs). Three-hundred-and-forty MHC haplotypes were identified in 1453 Bovidae individuals, suggesting significantly higher polymorphism and heterozygosity compared with those based on NAMs. The ambitious boundaries in population differentiation (phylogenetic network, pairwise FST and STRUCTURE analyses) within and between species assessed using the MHC haplotypes were different from those revealed by NAMs associated closely with speciation, geographical distribution, domestication and management histories. In addition, the mean FST was significantly correlated negatively with the number of observed alleles (NA), observed (HO) and expected (HE) heterozygosity and polymorphism information content (PIC) (P < 0.05) in the MHC haplotype dataset while there was no correction of the mean FST estimates (P> 0.05) between the MHC haplotype and NAMs datasets. Analysis of molecular variance (AMOVA) revealed a lower percentage of total variance (PTV) between species/groups based on the MHC-linked microsatellites than NAMs. Therefore, it was inferred that individuals within populations accumulated as many MHC variants as possible to increase their heterozygosity and thus the survival rate of their affiliated populations and species, which eventually reduced population differentiation and thereby complicated their classification and phylogenetic relationship inference. In summary, host-pathogen coevolution and heterozygote advantage, rather than demographic history, act as key driving forces shaping the MHC diversity within the populations and determining the interspecific MHC diversity.

Temporal expression of the KISS1/GPR54 system in goats' testes and epididymides and its spatial expression in pubertal goats.

Kisspeptin, encoded by the KISS1 gene, and its receptor GPR54 are essential in puberty onset and male fertility due to their central regulatory roles. However, the roles of KISS1/GPR54 in peripheral tissues remain unclear. This study aimed to investigate the temporal expression patterns of KISS1/GPR54 in goat testes and epididymides and its spatial expression patterns in pubertal goats. Immunohistochemical analysis revealed that kisspeptin/GPR54 were localized in Leydig, Sertoli, and germ cells of pubertal goats' testis, as well as in principal and basal cells of the epididymis. RT-PCR revealed a marked variation in the KISS1/GPR54 expressions in the testes and epididymides from the age of first week to adulthood. KISS1 and GPR54 mRNA levels in testes decreased from the age of first week to two months and then increased from two months to puberty and adulthood. The KISS1 and GPR54 mRNA levels in Leydig cells decreased from the age of one week to two months and increased from two months to puberty, and then decreased from puberty to adulthood. Only GPR54 mRNA levels in the epididymides increased from the age of one week to two months and puberty, and then decreased from puberty to adulthood. RT-PCR analysis showed the different spatial expression patterns of KISS1/GPR54 in pubertal goat tissues. The KISS1 mRNA level was high in the hypothalamus, moderate in pancreas, liver, epididymis and testis; and low in the other tissues. The GPR54 expression was high in the pancreas and testis; moderate in pituitary, hypothalamus and mesenteric lymph node; and low in the other tissues. In conclusion, the KISS1/GPR54 system possessed distinct temporal expression profiles in goats' testes and epididymides, as well as different spatial expression patterns in pubertal goat tissues, which implied the possible local role of this system in goats' testes, epididymides, and other peripheral tissues.

Local expressions and function of Kiss1/GPR54 in goats' testes.

Kisspeptin, encoded by the Kiss1 gene, and its receptor GPR54 have a central regulatory role in the male reproduction. However, their functions in peripheral tissues, such as testes, remain unclear. This study investigated the local expressions and function of Kiss1/GPR54 in goats' testes. The mRNA expression of Kiss1/GPR54 in pubertal goat Leydig cells was detected through reverse transcriptase PCR (RT-PCR), and its protein expression in Leydig cells or the testis was examined by immunohistochemistry and Western blot analysis. Isolated and cultured Leydig cells were treated with different concentration of kisspeptin (0, 1, 10 and 100 µM) and kisspeptin antagonist for 4, 24 or 72 h. The direct effect of kisspeptin on testosterone secretion and Kiss1/GPR54 mRNA expression was evaluated by ELISA and RT-PCR. Kiss1/GPR54 mRNA and protein were expressed in Leydig cells and spermatids, and GPR54 were expressed in Sertoli cells. Kisspeptin treatment significantly stimulated testosterone secretion in Leydig cells, with the highest levels found under 24 h of treatment with 10 µM kisspeptin. Treatment with kisspeptin + kisspeptin antagonist significantly reduced the kisspeptin-stimulated testosterone secretion in Leydig cells. Kisspeptin treatment significantly enhanced the expression of Kiss1/GPR54 mRNA in Leydig cells. These data suggest the local expressions of Kiss1/GPR54 in goats' testes and its autocrine role in Leydig cells, which is helpful in understanding the regulation role of kisspeptin/GPR54 system in other peripheral tissues.

Identification of a Goat Intersexuality-Associated Novel Variant Through Genome-Wide Resequencing and Hi-C.

Background: Polled intersex syndrome (PIS) leads to reproductive disorders in goats and exerts a heavy influence on goat breeding. Since 2001, the core variant of an 11.7 kb deletion at ~129 Mb on chromosome 1 (CHI1) has been widely used as a genetic diagnostic criterion. In 2020, a ~0.48 Mb insertion within the PIS deletion was identified by sequencing in XX intersex goats. However, the suitability of this variation for the diagnosis of intersex goats worldwide and its further molecular genetic mechanism need to be clarified. Results: The whole-genome selective sweep of intersex goats from China was performed with whole-genome next-generation sequencing technology for large sample populations and a case-control study on interbreeds. A series of candidate genes related to the goat intersexuality phenotype were found. We further confirmed that a ~0.48 Mb duplicated fragment (including ERG and KCNJ15) downstream of the ~20 Mb PIS region was reversely inserted into the PIS locus in intersex Chinese goats and was consistent with that in European Saanen and Valais black-necked goats. High-throughput chromosome conformation capture (Hi-C) technology was then used to compare the 3D structures of the PIS variant neighborhood in CHI1 between intersex and non-intersex goats. A newly found structure was validated as an intrachromosomal rearrangement. This inserted duplication changed the original spatial structure of goat CHI1 and caused the appearance of several specific loop structures in the adjacent ~20 kb downstream region of FOXL2. Conclusions: Results suggested that the novel complex PIS variant genome was sufficient as a broad-spectrum clinical diagnostic marker of XX intersexuality in goats from Europe and China. A series of private dense loop structures caused by segment insertion into the PIS deletion might affect the expression of FOXL2 or other neighboring novel candidate genes. However, these structures require further in-depth molecular biological experimental verification. In general, this study provided new insights for future research on the molecular genetic mechanism underlying female-to-male sex reversal in goats.