Transcripts of enriched germ cells responding to heat shock as potential markers for porcine semen quality

Integration of ssGWAS and ROH analyses for uncovering genetic variants associated with reproduction traits in Large White pigs

The low heritability of reproduction traits such as total number born (TNB), number born alive (NBA) and adjusted litter weight until 21 days at weaning (ALW) poses a challenge for genetic improvement. In this study, we aimed to identify genetic variants that influence these traits and evaluate the accuracy of genomic selection (GS) using these variants as genomic features. We performed single-step genome-wide association studies (ssGWAS) on 17 823 Large White (LW) pigs, of which 2770 were genotyped by 50K single nucleotide polymorphism (SNP) chips. Additionally, we analyzed runs of homozygosity (ROH) in the population and tested their effects on the traits. The genomic feature best linear unbiased prediction (GFBLUP) was then carried out in an independent population of 350 LW pigs using identified trait-related SNP subsets as genomic features. As a result, our findings identified five, one and four SNP windows that explaining more than 1% of genetic variance for ALW, TNB, and NBA, respectively and discovered 358 hotspots and nine ROH islands. The ROH SSC1:21814570-27186456 and SSC11:7220366-14276394 were found to be significantly associated with ALW and NBA, respectively. We assessed the genomic estimated breeding value accuracy through 20 replicates of five-fold cross-validation. Our findings demonstrate that GFBLUP, incorporating SNPs located in effective ROH (p-value < 0.05) as genomic features, might enhance GS accuracy for ALW compared with GBLUP. Additionally, using SNPs explaining more than 0.1% of the genetic variance in ssGWAS for NBA as genomic features might improve the GS accuracy, too. However, it is important to note that the incorporation of inappropriate genomic features can significantly reduce GS accuracy. In conclusion, our findings provide valuable insights into the genetic mechanisms of reproductive traits in pigs and suggest that the ssGWAS and ROH have the potential to enhance the accuracy of GS for reproductive traits in LW pigs.

Transcriptome Analysis Reveals Spermatogenesis-Related CircRNAs and LncRNAs in Goat Spermatozoa

Mammalian spermatozoa comprises both coding and non-coding RNAs, which are traditionally believed to be a residual of spermatogenesis. The differential expression level of spermatozoal RNAs is also observed between fertile and infertile human, thereby anticipated as potential molecular marker of male fertility. This study investigated the transcriptome profile of goat (Capra hircus) spermatozoa. The sperm transcriptome was analyzed by three different methods viz. RLM-RACE, long-read RNA sequencing (RNAseq) in Nanopore™ platform, and short-read RNAseq in Illumina™ platform. The Illumina™ sequencing discovered 16,604 transcripts with 357 mRNAs having FPKM (fragments per kilobase per million mapped reads) of more than five. The spermatozoal RNA suite included mRNA (94%), rRNA (3%), miscRNA (1%), circRNA (1%), miRNA (1%), etc. This study also predicted circRNAs (127), lncRNAs (655), and imprinted genes (160) that have potential role in male reproduction. The gene ontology analysis revealed the involvement of spermatozoal RNA in regulating male meiosis (TET3, STAT5B), capacitation (ACRBP, CATSPER4), sperm motility (GAS8, TEKT2), spermatogenesis (ADAMTS2, CREB3L4), etc. The spermatozoal RNA were also associated with different biological pathways viz. Wnt signaling pathway, cAMP signaling pathway, AMPK signaling pathway, and MAPK signaling pathways having potential role in spermatogenesis. Overall, this study enlightened the suite of spRNA transcripts in goat and their relevance in male fertility for diagnostic approach.

Germ cell apoptosis and expression of Bcl-2 and Bax in porcine testis under normal and heat stress conditions

The aim of this study was to examine whether an elevated ambient temperature (37-40°C) had an effect on the apoptosis of germ cells and the expression of Bcl-2 and Bax in porcine testis. Six boars were used. Three boars were subjected to an elevated ambient temperature (37-40°C, 7days, 3h per day) as a heat stress (HS) group. The other 3 boars were kept in a room temperature house (20-27°C) as a control group. All boars were castrated and the testes were harvested. TUNEL assay was used for the detection of apoptotic cells. Immunohistochemistry, Western blotting and quantitative real-time PCR were used to analyze protein and mRNA levels of Bcl-2 and Bax in response to heat treatment. The results showed that apoptotic signals increased under heat stress conditions compared with the control (P<0.01), and the cell types most affected by heat treatment were spermatocytes and spermatids. In both the control and experimental groups, Bcl-2 was expressed in the cytoplasm and nucleus of spermatogonia, spermatocytes and differentiating spermatids and Bcl-2 preferentially localized close to the seminiferous tubule's luminal surface in late spermatocytes and spermatids. Compared with the control group, the expression levels of Bcl-2 protein and mRNA significantly increased in heat treatment group, while the expression levels of Bax protein and mRNA did not show significant changes between the control and experimental group. Low to moderate Bax immunoreactivity staining was observed in all kinds of germ cells in the control group. Strong staining was observed in spermatogonia, and low to moderate Bax staining was observed in spermatocytes and spermatids. A redistribution of Bax from a cytoplasmic to perinuclear or nuclear localization could be observed in the spermatogonia, spermatocytes and spermatids obtained in the heat treated group. These results showed that elevated ambient temperatures induced germ cell apoptosis. In response to heat stress, the expression of Bcl-2 increased and a redistribution of Bax from a cytoplasmic to a perinuclear or nuclear localization. This indicates that Bcl-2 and Bax may be involved in regulation of germ cell apoptosis induced by heat stress in boars.

The human RHOX gene cluster: target genes and functional analysis of gene variants in infertile men

The X-linked reproductive homeobox (RHOX) gene cluster encodes transcription factors preferentially expressed in reproductive tissues. This gene cluster has important roles in male fertility based on phenotypic defects of Rhox-mutant mice and the finding that aberrant RHOX promoter methylation is strongly associated with abnormal human sperm parameters. However, little is known about the molecular mechanism of RHOX function in humans. Using gene expression profiling, we identified genes regulated by members of the human RHOX gene cluster. Some genes were uniquely regulated by RHOXF1 or RHOXF2/2B, while others were regulated by both of these transcription factors. Several of these regulated genes encode proteins involved in processes relevant to spermatogenesis; e.g. stress protection and cell survival. One of the target genes of RHOXF2/2B is RHOXF1, suggesting cross-regulation to enhance transcriptional responses. The potential role of RHOX in human infertility was addressed by sequencing all RHOX exons in a group of 250 patients with severe oligozoospermia. This revealed two mutations in RHOXF1 (c.515G > A and c.522C > T) and four in RHOXF2/2B (-73C > G, c.202G > A, c.411C > T and c.679G > A), of which only one (c.202G > A) was found in a control group of men with normal sperm concentration. Functional analysis demonstrated that c.202G > A and c.679G > A significantly impaired the ability of RHOXF2/2B to regulate downstream genes. Molecular modelling suggested that these mutations alter RHOXF2/F2B protein conformation. By combining clinical data with in vitro functional analysis, we demonstrate how the X-linked RHOX gene cluster may function in normal human spermatogenesis and we provide evidence that it is impaired in human male fertility.

Neural expression of the transcription factor THAP1 during development in rat

Loss of function mutations in THAP1 has been associated with primary generalized and focal dystonia in children and adults. THAP1 encodes a transcription factor (THAP1) that harbors an atypical zinc finger domain and plays a critical role in G(1)-S cell cycle control. Current thinking suggests that dystonia may be a neurodevelopmental circuit disorder. Hence, THAP1 may participate in the development of the nervous system. Herein, we report the neurodevelopmental expression patterns of Thap1 transcript and THAP1 protein from the early postnatal period through adulthood in the rat brain, spinal cord and dorsal root ganglia (DRG). We detected Thap1 transcript and THAP1-immunoreactivity (IR) in the cerebral cortex, cerebellum, striatum, substantia nigra, thalamus, spinal cord and DRG. Thap1 transcript expression was higher in the brain than in spinal cord and DRG at P1 and P7 and declined to similar levels at P14 and later time points in all regions except the cerebellum, where it remained high through adulthood. In the brain, THAP1 expression was highest in early development, particularly in the cerebellum at P7. In addition to Purkinje cells in the cerebellum, THAP1-IR was also localized to pyramidal neurons in the cerebral cortex, relay neurons in the thalamus, medium spiny and cholinergic neurons in the striatum, dopaminergic neurons in the substantia nigra, and pyramidal and interneurons in the hippocampus. In the cerebellar cortex, THAP1-IR was prominently distributed in the perikarya and proximal dendrites of Purkinje cells at early time-points. In contrast, it was more diffusely distributed throughout the dendritic arbor of adult Purkinje cells producing a moderate diffuse staining pattern in the molecular layer. At all time points, nuclear IR was weaker than cytoplasmic IR. The prominent cytoplasmic and developmentally regulated expression of THAP1 suggests that THAP1 may function as part of a cell surface-nucleus signaling cascade involved in terminal neural differentiation.