Comprehensive transcriptome analysis of hypothalamus reveals genes associated with disorders of sex development in pigs

A variant in long noncoding RNA NORSF affects granulosa cells response to transcription factor RFX7

NORSF is a nuclear long noncoding RNA (lncRNA) that contributes to the follicular atresia and restrains 17β-estradiol (E2) release by granulosa cells (GCs). Importantly, it is also a potential candidate gene in the quantitative trait locus (QTLs) for sow fertility traits. We identified NORSF as a candidate (causal) gene affecting sow fertility traits. A novel G-A variant was discovered at -478 nt of the NORSF promoter and termed as g.-478G>A. Association analysis revealed that this variant was associated with sow fertility traits (e.g., the total number of piglets born, the total number of piglets born alive, and the number of healthy piglets). Mechanistically, the g.-478G>A variant reduced the binding activity of the NORSF promoter to its transcription activator regulatory factor X7 (RFX7), leading to decreased NORSF promoter activity and transcription levels in sow GCs (sGCs), and weakened inhibitory effects on the transcription of CYP19A1, which encodes a rate-limiting enzyme for E2 synthesis and E2 release by sGCs. In addition, RFX7 is transcriptionally activated by P53, which restrains E2 release from sGCs via the RFX7/NORSF/CYP19A1 pathway. These findings indicate that the lncRNA NORSF is a causal gene in QTLs for sow fertility traits and define the P53/NORSF/CYP19A1 pathway as a new signaling pathway affecting sow reproduction, which provides a new target for improving female fertility.

Analysis of 206 whole-genome resequencing reveals selection signatures associated with breed-specific traits in Hu sheep

As an invaluable Chinese sheep germplasm resource, Hu sheep are renowned for their high fertility and beautiful wavy lambskins. Their distinctive characteristics have evolved over time through a combination of artificial and natural selection. Identifying selection signatures in Hu sheep can provide a straightforward insight into the mechanism of selection and further uncover the candidate genes associated with breed-specific traits subject to selection. Here, we conducted whole-genome resequencing on 206 Hu sheep individuals, each with an approximate 6-fold depth of coverage. And then we employed three complementary approaches, including composite likelihood ratio, integrated haplotype homozygosity score and the detection of runs of homozygosity, to detect selection signatures. In total, 10 candidate genomic regions displaying selection signatures were simultaneously identified by multiple methods, spanning 88.54 Mb. After annotating, these genomic regions harbored collectively 92 unique genes. Interestingly, 32 candidate genes associated with reproduction were distributed in nine genomic regions detected. Out of them, two stood out as star candidates: BMPR1B and GNRH2, both of which have documented associations with fertility, and a HOXA gene cluster (HOXA1-5, HOXA9, HOXA10, HOXA11 and HOXA13) had also been linked to fertility. Additionally, we identified other genes that are related to hair follicle development (LAMTOR3, EEF1A2), ear size (HOXA1, KCNQ2), fat tail formation (HOXA10, HOXA11), growth and development (FAF1, CCNDBP1, GJB2, GJA3), fat deposition (ACOXL, JAZF1, HOXA3, HOXA4, HOXA5, EBF4), immune (UBR1, FASTKD5) and feed intake (DAPP1, RNF17, NPBWR2). Our results offer novel insights into the genetic mechanisms underlying the selection of breed-specific traits in Hu sheep and provide a reference for sheep genetic improvement programs.

Epigenetics of inflammation in hypothalamus pituitary gonadal and neuroendocrine disorders

The hormone producing hypothalamus, pituitary and gonadal are arranged in hierarchy to form the hypothalamic-pituitary-gonadal axis (HPG axis). The axis is neuroendocrine in nature and releases hormones in response to the inputs from nervous systems. The axis maintains homeostasis and ensures smooth body functions, particularly those related to growth and reproduction. A deregulated HPG axis, such as observed under inflammation and other conditions, is therefore associated with several disorders such as polycystic ovary syndrome, functional hypothalamic amenorrhea etc. Several factors, both genetic as well as environmental, in addition to aging, obesity etc. affect HPG axis with resulting effects on puberty, sexual maturation and reproductive health. More research is now indicative of a role of epigenetics in mediating these HPG-affecting factors. Hypothalamus-secreted gonadotropin-releasing hormone is important for eventual release of sex hormones and it is subjected to several neuronal and epigenetic regulations. Gene promoter methylation as well as histone methylations and acetylations form the backbone of epigenetic regulation of HPG-axis, as the incoming reports suggest. Epigenetic events also mediate several feedback mechanisms within HPG axis and between HPG axis and the central nervous system. In addition, data is emerging for a role of non-coding RNAs, particularly the miRNAs, in regulation and normal functioning of HPG axis. Thus, the epigenetic interactions need better understanding to understand the functioning and regulation of HPG axis.

Hypothalamic transcriptome analysis reveals male-specific differences in molecular pathways related to oxidative phosphorylation between Iberian pig genotypes

The hypothalamus is implicated in controlling feeding and adiposity, besides many other physiological functions, and thus can be of great importance in explaining productive differences between lean and fatty pig breeds. The present study aimed to evaluate the hypothalamic transcriptome of pure Iberian (IBxIB) and Large White x Iberian crossbreds (IBxLW) at 60 days-old, produced in a single maternal environment. Results showed the implication of gender and genotype in the hypothalamic transcriptome, with 51 differentially expressed genes (DEGs) between genotypes and 10 DEGs between genders. Fourteen genotype by sex interactions were found, due to a higher genotype effect on transcriptome found in males. In fact, just 31 DEGs were identified when using only females but 158 using only males. A higher expression of genes related to mitochondrial activity in IBxIB male animals (ND3, ND4, ND5, UQCRC2 and ATP6) was found, which was related to a higher oxidative phosphorylation and greater reactive oxygen species and nitric oxide production. IBxLW male animals showed higher expression of SIRT3 regulator, also related to mitochondrial function. When females were analysed, such differences were not found, since only some differences in genes related to the tricarboxylic acid cycle. Thus, the results indicate a significant effect and interaction of the breed and the sex on the hypothalamic transcriptome at this early age.

Comprehensive Transcriptome Analysis of Follicles from Two Stages of the Estrus Cycle of Two Breeds Reveals the Roles of Long Intergenic Non-Coding RNAs in Gilts

Simple Summary

This study provides new perspectives about the roles of lincRNAs in the estrus expression of gilts, which is correlated with ovarian steroid hormone and follicular development. Follicular tissues from two stages of the estrus cycle of Large White and Mi gilts were used for RNA-seq. Some genes and lincRNAs related to estrus expression in pigs were discovered. PPI and ceRNA networks related to the estrus expression were constructed. These results suggest that the estrus expression may be affected by lincRNAs and their target genes.

Abstract

Visible and long-lasting estrus expression of gilts and sows effectively sends a mating signal. To reveal the roles of Long Intergenic Non-coding RNAs (lincRNAs) in estrus expression, RNA-seq was used to investigate the lincRNAs expression of follicular tissues from Large White gilts at diestrus (LD) and estrus (LE), and Chinese Mi gilts at diestrus (MD) and estrus (ME). Seventy-three differentially expressed lincRNAs (DELs) were found in all comparisons (LE vs. ME, LD vs. LE, and MD vs. ME comparisons). Eleven lincRNAs were differentially expressed in both LD vs. LE and MD vs. ME comparisons. Fifteen DELs were mapped onto the pig corpus luteum number Quantitative Trait Loci (QTL) fragments. A protein–protein interaction (PPI) network that involved estrus expression using 20 DEGs was then constructed. Interestingly, three predicted target DEGs (PTGs) (CYP19A1 of MSTRG.10910, CDK1 of MSTRG.10910 and MSTRG.23984, SCARB1 of MSTRG.1559) were observed in the PPI network. A competitive endogenous RNA (ceRNA) network including three lincRNAs, five miRNAs, and five genes was constructed. Our study provides new insight into the lincRNAs associated with estrus expression and follicular development in gilts.

Transcriptomic Analysis Reveals Functional Interaction of mRNA-lncRNA-miRNA in Steroidogenesis and Spermatogenesis of Gynogenetic Japanese Flounder (Paralichthys olivaceus)

Teleost fishes exhibit extraordinary diversity, plasticity and adaptability with their sex determination and sexual development, and there is growing evidence that non-coding RNAs (ncRNAs) are emerging as critical regulators of reproduction. Japanese flounder (Paralichthys olivaceus) is an important marine cultured fish that presents significant sexual dimorphism with bigger females, in which gynogenesis has been applied for aquaculture industry. In order to reveal the regulatory mechanisms of sexual development in gynogenetic female and sex-reversed neo-male P. olivaceus, the lncRNA-miRNA-mRNA interactions were investigated using high-throughput sequencing. A total of 6772 differentially expressed mRNAs (DEmRNAs), 2284 DElncRNAs, and 244 DEmiRNAs were obtained between gynogenetic female ovaries and sex-reversed neo-male testes. Genes in the steroid hormone biosynthesis and secretion pathway were enriched and mostly significantly upregulated in neo-male testes. Subsequently, network analysis uncovered high functional specificity for gynogenetic P. olivaceus sperm motility, as co-expressed DEmRNAs were significantly enriched in microtubule and cytoskeleton-related biological processes. Clustered miRNAs were characterized in the P. olivaceus genome with examples of the largest conserved let-7 clusters. The 20 let-7 members are distributed in 11 clusters and may not transcribe together with their neighboring miR-125b, with let-7 repressing cyp11a and miR-125b repressing esr2b, both as key steroidogenesis pathway genes. In summary, this study provides comprehensive insights into the mRNA-miRNA-lncRNA functional crosstalk in teleost sexual development and gametogenesis and will expand our understanding of ncRNA biology in teleost gynogenesis.

Analyses of Molecular Characteristics and Enzymatic Activities of Ovine HSD17B3

17β-hydroxysteroid dehydrogenase type 3 (HSD17B3) converts androstenedione (A4) into testosterone (T), which regulates sex steroid production. Because various mutations of the HSD17B3 gene cause disorder of sex differentiation (DSD) in multiple mammalian species, it is very important to reveal the molecular characteristics of this gene in various species. Here, we revealed the open reading frame of the ovine HSD17B3 gene. Enzymatic activities of ovine HSD17B3 and HSD17B1 for converting A4 to T were detected using ovine androgen receptor-mediated transactivation in reporter assays. Although HSD17B3 also converted estrone to estradiol, this activity was much weaker than those of HSD17B1. Although ovine HSD17B3 has an amino acid sequence that is conserved compared with other mammalian species, it possesses two amino acid substitutions that are consistent with the reported variants of human HSD17B3. Substitutions of these amino acids in ovine HSD17B3 for those in human did not affect the enzymatic activities. However, enzymatic activities declined upon missense mutations of the HSD17B3 gene associated with 46,XY DSD, affecting amino acids that are conserved between these two species. The present study provides basic information and tools to investigate the molecular mechanisms behind DSD not only in ovine, but also in various mammalian species.

Whole genome sequencing identifies a cryptic SOX9 regulatory element duplication underlying a case of 46,XX ovotesticular difference of sexual development

Ovotesticular differences of sexual development (OT‐DSD) are rare genetic variances defined by the coexistence of both testicular and ovarian tissues. Various molecular etiologies including SRY translocation or SOX9 pathogenic variants with different modes of inheritance have been associated with 46,XX OT‐DSD. Here we describe a child diagnosed with SRY‐negative 46,XX OT‐DSD after completing a series of complex clinical genetic analyses, including chromosomal microarray, DSD gene panel (sequencing and deletion/duplication analysis), whole exome sequencing, and whole genome sequencing. Of these, only whole genome sequencing reported a pathogenic duplication in a non‐coding region that contains the RevSex regulatory element, which modifies SOX9 expression and is associated with 46,XX OT‐DSD and complete sex reversal. This is the first clinical RevSex duplication detected by clinical whole genome sequencing. We highlight the utility of whole genome sequencing in shortening the diagnostic odyssey and the importance of optimal counseling through a team‐based multi‐specialty approach for patients with DSDs.