Characterization of sheep spermatogenesis through single-cell RNA sequencing

Adult Human, but Not Rodent, Spermatogonial Stem Cells Retain States with a Foetal-like Signature

Spermatogenesis involves a complex process of cellular differentiation maintained by spermatogonial stem cells (SSCs). Being critical to male reproduction, it is generally assumed that spermatogenesis starts and ends in equivalent transcriptional states in related species. Based on single-cell gene expression profiling, it has been proposed that undifferentiated human spermatogonia can be subclassified into four heterogenous subtypes, termed states 0, 0A, 0B, and 1. To increase the resolution of the undifferentiated compartment and trace the origin of the spermatogenic trajectory, we re-analysed the single-cell (sc) RNA-sequencing libraries of 34 post-pubescent human testes to generate an integrated atlas of germ cell differentiation. We then used this atlas to perform comparative analyses of the putative SSC transcriptome both across human development (using 28 foetal and pre-pubertal scRNA-seq libraries) and across species (including data from sheep, pig, buffalo, rhesus and cynomolgus macaque, rat, and mouse). Alongside its detailed characterisation, we show that the transcriptional heterogeneity of the undifferentiated spermatogonial cell compartment varies not only between species but across development. Our findings associate 'state 0B' with a suppressive transcriptomic programme that, in adult humans, acts to functionally oppose proliferation and maintain cells in a ready-to-react state. Consistent with this conclusion, we show that human foetal germ cells-which are mitotically arrested-can be characterised solely as state 0B. While germ cells with a state 0B signature are also present in foetal mice (and are likely conserved at this stage throughout mammals), they are not maintained into adulthood. We conjecture that in rodents, the foetal-like state 0B differentiates at birth into the renewing SSC population, whereas in humans it is maintained as a reserve population, supporting testicular homeostasis over a longer reproductive lifespan while reducing mutagenic load. Together, these results suggest that SSCs adopt differing evolutionary strategies across species to ensure fertility and genome integrity over vastly differing life histories and reproductive timeframes.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

Single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics offer new insights into the etiological basis for male cattle-yak sterility

The variety of species can be efficiently increased by interspecific hybridization. However, because the males in the hybrid progeny are usually sterile, this heterosis cannot be employed when other cattle and yaks are hybridized. While some system-level studies have sought to explore the etiological basis for male cattle-yak sterility, no systematic cellular analyses of this phenomenon have yet been performed. Here, single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics methods were used to study the differences in testicular tissue between 4-year-old male yak and 4-year-old male cattle-yak, providing new and comprehensive insights into the causes of male cattle-yak sterility. Cattle-yak testes samples detected 6 somatic cell types and one mixed germ cell type. Comparisons of these cell types revealed the more significant differences in Sertoli cells (SCs) and [Leydig cells and myoid cells (LCs_MCs)] between yak and cattle-yak samples compared to other somatic cell clusters. Even though the LCs and MCs from yaks and cattle-yaks were derived from the differentiation of the same progenitor cells, a high degree of overlap between LCs and MCs was observed in yak samples. Still, only a small overlap between LCs and MCs was observed in cattle-yak samples. Functional enrichment analyses revealed that genes down-regulated in cattle-yak SCs were primarily enriched in biological activity, whereas up-regulated genes in these cells were enriched for apoptotic activity. Furthermore, the genes of up-regulated in LCs_MCs of cattle-yak were significantly enriched in enzyme inhibitor and molecular function inhibitor activity. On the other hand, the genes of down-regulated in these cells were enriched for signal receptor binding, molecular function regulation, positive regulation of biological processes, and regulation of cell communication activity. The most significant annotated differences between yak and cattle-yak LCs_MCs were associated with cell-to-cell communication. While yak LCs_MCs regulated spermatogenic cells at spermatogonia, spermatocyte, and spermatid levels, no such relationships were found between cattle-yak LCs_MCs and germ cells. This may suggest that the somatic niche in male cattle-yak testes is a microenvironment that is ultimately not favorable for spermatogenesis.

Integrated single cell transcriptome sequencing analysis reveals species-specific genes and molecular pathways for pig spermiogenesis

Mammalian spermatogenesis is a highly complicated and intricately organized process involving spermatogonia propagation (mitosis) and meiotic differentiation into mature sperm cells (spermiogenesis). In pigs, spermatogonia development and the role of somatic cells in spermatogenesis were previously investigated in detail. However, the characterization of key molecules fundamental to pig spermiogenesis remains less explored. Here we compared spermatogenesis between humans and pigs, focusing on spermiogenesis, by integrative testicular single-cell RNA sequencing (scRNA-seq) analysis. Human and pig testicular cells were clustered into 26 different groups, with cell-type-specific markers and signalling pathways. For spermiogenesis, pseudo-time analysis classified the lineage differentiation routes for round, elongated spermatids and spermatozoa. Moreover, markers and molecular pathways specific to each type of spermatids were examined for humans and pigs, respectively. Furthermore, high-dimensional weighted gene co-expression network analysis (hdWGCNA) identified gene modules specific for each type of human and pig spermatids. Hub genes (pig: SNRPD2.1 related to alternative splicing; human: CATSPERZ, Ca[2+] ion channel) potentially involved in spermiogenesis were also revealed. Taken together, our integrative analysis found that human and pig spermiogeneses involve specific genes and molecular pathways and provided resources and insights for further functional investigation on spermatid maturation and male reproductive ability.

Molecular characterization, expression patterns and cellular localization of BCAS2 gene in male Hezuo pig

Background

Breast carcinoma amplified sequence 2 (BCAS2) participates in pre-mRNA splicing and DNA damage response, which is implicated in spermatogenesis and meiosis initiation in mouse. Nevertheless, the physiological roles of BCAS2 in the testes of large mammals especially boars remain largely unknown.

Methods

In this study, testes were collected from Hezuo pig at three development stages including 30 days old (30 d), 120 days old (120 d), and 240 days old (240 d). BCAS2 CDS region was firstly cloned using RT-PCR method, and its molecular characteristics were identified using relevant bioinformatics software. Additionally, the expression patterns and cellular localization of BCAS2 were analyzed by quantitative real-time PCR (qRT-PCR), Western blot, immunohistochemistry and immunofluorescence.

Results

The cloning and sequence analysis indicated that the Hezuo pig BCAS2 CDS fragment encompassed 678 bp open reading frame (ORF) capable of encoding 225 amino acid residues, and possessed high identities with some other mammals. The results of qRT-PCR and Western blot displayed that BCAS2 levels both mRNA and protein were age-dependent increased (p < 0.01). Additionally, immunohistochemistry and immunofluorescence results revealed that BCAS2 protein was mainly observed in nucleus of gonocytes at 30 d testes as well as nucleus of spermatogonia and Sertoli cells at 120 and 240 d testes. Accordingly, we conclude that BCAS2 is critical for testicular development and spermatogenesis of Hezuo pig, perhaps by regulating proliferation or differentiation of gonocytes, pre-mRNA splicing of spermatogonia and functional maintenance of Sertoli cells, but specific mechanism still requires be further investigated.

Whole-genome resequencing reveals genetic diversity and selection characteristics of dairy goat

The dairy goat is one of the earliest dairy livestock species, which plays an important role in the economic development, especially for developing countries. With the development of agricultural civilization, dairy goats have been widely distributed across the world. However, few studies have been conducted on the specific characteristics of dairy goat. In this study, we collected the whole-genome data of 89 goat individuals by sequencing 48 goats and employing 41 publicly available goats, including five dairy goat breeds (Saanen, Nubian, Alpine, Toggenburg, and Guanzhong dairy goat; n = 24, 15, 11, 6, 6), and three goat breeds (Guishan goat, Longlin goat, Yunshang Black goat; n = 6, 15, 6). Through compared the genomes of dairy goat and non-dairy goat to analyze genetic diversity and selection characteristics of dairy goat. The results show that the eight goats could be divided into three subgroups of European, African, and Chinese indigenous goat populations, and we also found that Australian Nubian, Toggenburg, and Australian Alpine had the highest linkage disequilibrium, the lowest level of nucleotide diversity, and a higher inbreeding coefficient, indicating that they were strongly artificially selected. In addition, we identified several candidate genes related to the specificity of dairy goat, particularly genes associated with milk production traits (GHR, DGAT2, ELF5, GLYCAM1, ACSBG2, ACSS2), reproduction traits (TSHR, TSHB, PTGS2, ESR2), immunity traits (JAK1, POU2F2, LRRC66). Our results provide not only insights into the evolutionary history and breed characteristics of dairy goat, but also valuable information for the implementation and improvement of dairy goat cross breeding program.

Single-cell RNA sequencing uncovers dynamic roadmap and cell-cell communication during buffalo spermatogenesis

Spermatogenesis carries the task of precise intergenerational transmission of genetic information from the paternal genome and involves complex developmental processes regulated by the testicular microenvironment. Studies performed mainly in mouse models have established the theoretical basis for spermatogenesis, yet the wide interspecies differences preclude direct translation of the findings, and farm animal studies are progressing slowly. More than 32,000 cells from prepubertal (3-month-old) and pubertal (24-month-old) buffalo testes were analyzed by using single-cell RNA sequencing (scRNA-seq), and dynamic gene expression roadmaps of germ and somatic cell development were generated. In addition to identifying the dynamic processes of sequential cell fate transitions, the global cell-cell communication essential to maintain regular spermatogenesis in the buffalo testicular microenvironment was uncovered. The findings provide the theoretical basis for establishing buffalo germline stem cells in vitro or culturing organoids and facilitating the expansion of superior livestock breeding.

Sperm Transcriptome Analysis Accurately Reveals Male Fertility Potential in Livestock

Nowadays, selection of superior male candidates in livestock as a source of frozen semen based on sperm quality at the cellular level is not considered accurate enough for predicting the potential of male fertility. Sperm transcriptome analysis approaches, such as messenger RNA levels, have been shown to correlate with fertility rates. Using this technology in livestock growth has become the principal method, which can be widely applied to predict male fertility potential in the livestock industry through the analysis of the sperm transcriptome. It provides the gene expression to validate the function of sperm in spermatogenesis, fertilization, and embryo development, as the parameters of male fertility. This review proposes a transcriptomic analysis approach as a high-throughput method to predict the fertility potential of livestock more accurately in the future.

Ubiquitin C-terminal hydrolase L1 (UCHL1), a double-edged sword in mammalian oocyte maturation and spermatogenesis

BACKGROUND

Recent studies have shown that ubiquitin-mediated cell apoptosis can modulate protein interaction and involve in the progress of oocyte maturation and spermatogenesis. As one of the key regulators involved in ubiquitin signal, ubiquitin C-terminal hydrolase L1 (UCHL1) is considered a molecular marker associated with spermatogonia stem cells. However, the function of UCHL1 was wildly reported to regulate various bioecological processes, such as Parkinson's disease, lung cancer, breast cancer and colon cancer, how UCHL1 affects the mammalian reproductive system remains an open question.

METHODS

We identified papers through electronic searches of PubMed database from inception to July 2022.

RESULTS

Here, we summarize the important function of UCHL1 in controlling mammalian oocyte development, regulating spermatogenesis and inhibiting polyspermy, and we posit the balance of UCHL1 was essential to maintaining reproductive cellular and tissue homeostasis.

CONCLUSION

This study considers the 'double-edged sword' role of UCHL1 during gametogenesis and presents new insights into UCHL1 in germ cells.

Dosage Compensation of the X Chromosome during Sheep Testis Development Revealed by Single-Cell RNA Sequencing

Simple Summary

Male and female mammals carry the same complement of autosomes but differ with respect to their sex chromosomes: females carry XX chromosomes and males carry XY chromosomes. The evolutionary loss of genes from the Y chromosome led to a disparity in the dosage of X chromosomes versus autosomal genes, with males becoming monosomic for X-linked gene products. An imbalance in gene expression may have detrimental consequences. In males, X-linked genes need to be upregulated to levels equal to those of females, which is called dosage compensation. The existence of dosage compensation in germ cells is controversial. In testis, dosage compensation is thought to cease during meiosis. Some studies showed that the X chromosome is inactivated during meiosis and premature transcriptional inactivation of the X and Y chromosome during mid-spermatogenesis is essential for fertility. However, some studies failed to find support for male germline X inactivation. Using single-cell RNA seq data, in this study, we presented a comprehensive transcriptional map of sheep testes at different developmental stages and found that germ cell types in sheep testes show X-chromosome expression similar to that in the autosomes. The dosage compensation of germ cells at different stages was analyzed. MSL complex members are expressed in female flies and orthologs exist in many species, where dosage compensation mechanisms are absent or fundamentally different. This suggests that the MSL complex members also function outside of the dosage compensation machinery. Studies have shown that MSL complex can regulate mammalian X inactivation and activation.

Abstract

Dosage compensation is a mechanism first proposed by Susumu Ohno, whereby X inactivation balances X gene output between males (XY) and females (XX), while X upregulation balances X genes with autosomal gene output. These mechanisms have been actively studied in Drosophila and mice, but research regarding them lags behind in domestic species. It is unclear how the X chromosome is regulated in the sheep male germline. To address this, using single-cell RNA sequencing, we analyzed testes in three important developmental stages of sheep. We observed that the total RNA per cell from X and autosomes peaked in SSCs and spermatogonia and was then reduced in early spermatocytes. Furthermore, we counted the detected reads per gene in each cell type for X and autosomes. In cells experiencing dose compensation, close proximity to MSL (male-specific lethal), which is regulated the active X chromosome and was observed. Our results suggest that there is no dose compensation in the pre-meiotic germ cells of sheep testes and, in addition, MSL1 and MSL2 are expressed in early germ cells and involved in regulating mammalian X-chromosome inactivation and activation.

Single-cell transcriptomics reveals male germ cells and Sertoli cells developmental patterns in dairy goats

Spermatogenesis holds considerable promise for human-assisted reproduction and livestock breeding based on stem cells. It occurs in seminiferous tubules within the testis, which mainly comprise male germ cells and Sertoli cells. While the developmental progression of male germ cells and Sertoli cells has been widely reported in mice, much less is known in other large animal species, including dairy goats. In this study, we present the data of single cell RNA sequencing (scRNA-seq) for 25,373 cells from 45 (pre-puberty), 90 (puberty), and 180-day-old (post-puberty) dairy goat testes. We aimed to identify genes that are associated with key developmental events in male germ cells and Sertoli cells. We examined the development of spermatogenic cells and seminiferous tubules from 15, 30, 45, 60, 75, 90, 180, and 240-day-old buck goat testes. scRNA-seq clustering analysis of testicular cells from pre-puberty, puberty, and post-puberty goat testes revealed several cell types, including cell populations with characteristics of spermatogonia, early spermatocytes, spermatocytes, spermatids, Sertoli cells, Leydig cells, macrophages, and endothelial cells. We mapped the timeline for male germ cells development from spermatogonia to spermatids and identified gene signatures that define spermatogenic cell populations, such as AMH, SOHLH1, INHA, and ACTA2. Importantly, using immunofluorescence staining for different marker proteins (UCHL1, C-KIT, VASA, SOX9, AMH, and PCNA), we explored the proliferative activity and development of male germ cells and Sertoli cells. Moreover, we identified the expression patterns of potential key genes associated with the niche-related key pathways in male germ cells of dairy goats, including testosterone, retinoic acid, PDGF, FGF, and WNT pathways. In summary, our study systematically investigated the elaborate male germ cells and Sertoli cells developmental patterns in dairy goats that have so far remained largely unknown. This information represents a valuable resource for the establishment of goat male reproductive stem cells lines, induction of germ cell differentiation in vitro, and the exploration of sequential cell fate transition for spermatogenesis and testicular development at single-cell resolution.

Proteome Informatics in Tibetan Sheep (Ovis aries) Testes Suggest the Crucial Proteins Related to Development and Functionality

Testis has an indispensable function in male reproduction of domestic animals. Tibetan sheep (Ovis aries) is a locally adapted breed of sheep raised in the Qinghai-Tibet Plateau, with outsized roles in providing the livelihood for millions of residents. Nevertheless, less is known on how protein expression and their functional roles in developmental testes of such breed limit their use in breeding efforts. In this study, we obtained comprehensive protein profiles from testes of Tibetan sheep at three developmental stages (including pre-puberty, post-puberty, and adulthood) using data-independent acquisition-based proteomic strategy to quantitatively identify the differentially abundant proteins (DAPs) associated with testicular development and function and to unravel the molecular basis of spermatogenesis. A total of 6,221 proteins were differentially expressed in an age-dependent manner. The reliability of the gene expression abundance was corroborated by quantitative PCR and targeted parallel reaction monitoring. These DAPs were significantly enriched to biological processes concerning spermatid development and sperm deformation, mitosis, glycolytic process, cell-cell/extracellular matrix (ECM) junctions, cell proliferation, apoptosis, and migration and to the pathways including, developmental process and sexual reproduction-related (such as VEGF, estrogen, insulin, GnRH, Hippo, PI3K-Akt, mTOR, MAPK, and AMPK), and testicular cell events-related pathways (such as tight/gap/adherens junctions, ECM-receptor interaction, regulation of actin cytoskeleton, glycolysis, cell cycle, and meiosis). Based on these bioinformatics analysis, we constructed four protein–protein interaction network, among which the proteins are involved in mitosis, meiosis, spermiogenesis, and testicular microenvironment, respectively. Altogether, these bioinformatics-based sequencing results suggest that many protein-coding genes were expressed in a development-dependent manner in Tibetan sheep testes to contribute to the testicular cell development and their surrounding microenvironment remodeling at various stages of spermatogenesis. These findings have important implications for further understanding of the mechanisms underlying spermatogenesis in sheep and even other plateau-adapted animals.

Transcriptome Analysis Reveals Key Gene Expression Changes in Blue Catfish Sperm in Response to Cryopreservation

The hybrids of female channel catfish (Ictalurus punctatus) and male blue catfish (I. furcatus) account for >50% of US catfish production due to superior growth, feed conversion, and disease resistance compared to both parental species. However, these hybrids can rarely be naturally spawned. Sperm collection is a lethal procedure, and sperm samples are now cryopreserved for fertilization needs. Previous studies showed that variation in sperm quality causes variable embryo hatch rates, which is the limiting factor in hybrid catfish breeding. Biomarkers as indicators for sperm quality and reproductive success are currently lacking. To address this, we investigated expression changes caused by cryopreservation using transcriptome profiles of fresh and cryopreserved sperm. Sperm quality measurements revealed that cryopreservation significantly increased oxidative stress levels and DNA fragmentation, and reduced sperm kinematic parameters. The present RNA-seq study identified 849 upregulated genes after cryopreservation, including members of all five complexes in the mitochondrial electron transport chain, suggesting a boost in oxidative phosphorylation activities, which often lead to excessive production of reactive oxygen species (ROS) associated with cell death. Interestingly, functional enrichment analyses revealed compensatory changes in gene expression after cryopreservation to offset detrimental effects of ultra-cold storage: MnSOD was induced to control ROS production; chaperones and ubiquitin ligases were upregulated to correct misfolded proteins or direct them to degradation; negative regulators of apoptosis, amide biosynthesis, and cilium-related functions were also enriched. Our study provides insight into underlying molecular mechanisms of sperm cryoinjury and lays a foundation to further explore molecular biomarkers on cryo-survival and gamete quality.

Single-cell RNA sequencing of the Mongolia sheep testis reveals a conserved and divergent transcriptome landscape of mammalian spermatogenesis

Spermatogenesis is a highly coordinated and complex process, and is pivotal for transmitting genetic information between mammalian generations. In this study, we investigated the conservation, differences, and biological functions of homologous genes during spermatogenesis in Mongolia sheep, humans, cynomolgus monkey, and mice using single-cell RNA sequencing technology. We compared X chromosome meiotic inactivation events in Mongolia sheep, humans, cynomolgus monkey, and mice to uncover the concerted activity of X chromosome genes. Subsequently, we focused on the dynamics of gene expression, key biological functions, and signaling pathways at various stages of spermatogenesis in Mongolia sheep and humans. Additionally, the ligand-receptor networks of Mongolia sheep and humans in testicular somatic and germ cells at different developmental stages were mapped to reveal conserved germ cell-soma communication using single-cell resolution. These datasets provided novel information and insights to unravel the molecular regulatory mechanisms of Mongolia sheep spermatogenesis and highlight conservation in gene expression during spermatogenesis between Mongolia sheep and humans, providing a foundation for the establishment of a large mammalian disease model of male infertility.

Comparative Toxicotranscriptomics of Single Cell RNA-Seq and Conventional RNA-Seq in TCDD-Exposed Testicular Tissue

In this report, we compare the outcomes and limitations of two methods of transcriptomic inquiry on adult zebrafish testes exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) during sexual differentiation: conventional or bulk RNA-seq (bulk-seq) and single cell RNA sequencing (scRNA-seq) data. scRNA-seq has emerged as a valuable tool for uncovering cell type-specific transcriptome dynamics which exist in heterogeneous tissue. Our lab previously showed the toxicological value of the scRNA-seq pipeline to characterize the sequelae of TCDD exposure in testes, demonstrating that loss of spermatids and spermatozoa, but not other cell types, contributed to the pathology of infertility in adult male zebrafish exposed during sexual differentiation. To investigate the potential for technical artifacts in scRNA-seq such as cell dissociation effects and reduced transcriptome coverage, we compared bulk-sequenced and scRNA-seq-paired samples from control and TCDD-exposed samples to understand what is gained and lost in scRNA-seq vs bulk-seq, both transcriptomically and toxicologically. We hypothesized that the testes may be sensitive to tissue disruption as they contain multiple cell types under constant division and/or maturation, and that TCDD exposure may mediate the extent of sensitivity. Thus, we sought to understand the extent to which this dissociation impacts the toxicological value of data returned from scRNA-seq. We confirm that the required dissociation of individual cells from intact tissue has a significant impact on gene expression, affecting gene pathways with the potential to confound toxicogenomics studies on exposures if findings are not well-controlled and well-situated in context. Additionally, a common scRNA-seq method using cDNA amplified from the 3’ end of mRNA under-detects low-expressing transcripts including transcription factors. We confirm this, and show TCDD-related genes may be overlooked by scRNA-seq, however, this under-detection effect is not mediated by TCDD exposure. Even so, scRNA-seq generally extracted toxicologically relevant information better than the bulk-seq method in the present study. This report aims to inform future experimental design for transcriptomic investigation in the growing field of toxicogenomics by demonstrating the differential information extracted from sequencing cells—despite being from the same tissue and exposure scheme—is influenced by the specific protocol used, with implications for the interpretation of exposure-induced risk.

Identification of Genes Related to Squab Muscle Growth and Lipid Metabolism from Transcriptome Profiles of Breast Muscle and Liver in Domestic Pigeon (Columba livia)

Simple Summary

Domestic pigeon is an important small poultry species raised for high-quality meat production. However, the relevant gene associated with meat growth and lipid metabolism during the period from dehulling to marketing are not known. Therefore, we aim to identify genes related to squab muscle growth and lip metabolism from transcriptome profiles of breast muscle and liver in domestic pigeon. In this study, we totally found that 4465 differentially expressed genes (DEGs) identified in the breast muscle and liver libraries, which include 2585 genes that were up-regulated and 2122 genes that were down-regulated. Most genes are involved in cell proliferation and differentiation, lipid metabolism and energy metabolism according to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs. We also detected 16 DEGs to verify data from RNA-Seq and qPCR, which were consistent in quantitatively estimating the transcription levels of the tested transcripts by qPCR analysis. The results of this study will lay the foundation for understanding the mechanisms of muscle growth and lipid metabolism in domestic pigeons.

Abstract

The improvements in muscle growth rate and meat quality are the major breeding aims in pigeon industry. Liver and muscle are recognized as important sites for fatty acid metabolism; understanding the role of specific transcripts in the breast muscle and liver might lead to the elucidation of interrelated biological processes. In this study, RNA-Sequencing (RNA-Seq) was applied to compare the transcriptomes of breast muscle and liver tissues among pigeons at five developmental periods (0, 1, 2, 3, 4 weeks post-hatching) to identify candidate genes related to muscle growth and lipid metabolism. There were 3142 differentially expressed genes (DEGs) identified in the breast muscle libraries; 1794 genes were up-regulated while 1531 genes were down-regulated. A total of 1323 DEGs were acquired from the liver libraries, with 791 up-regulated genes and 591 down-regulated genes. By pathway enrichment analysis, a set of significantly enriched pathways were identified for the DEGs, which are potentially involved in cell proliferation and differentiation, lipid metabolism and energy metabolism in pigeon breast muscle and liver. Our results are consistent with previous partial reports from domestic animals and poultry and provide some unidentified genes involved in muscle growth and lipid metabolism. The reliability of the sequencing data was verified through qPCR analysis of 16 genes from eight comparison groups (two genes per group). The findings from this study could contribute to future investigations of muscle growth and lipid metabolism mechanisms and establish molecular approaches to improve muscle growth rate and meat quality in domestic pigeon breeding.

Transcription Analysis for Core Networks of lncRNAs–mRNAs: Implication for Potential Role in Sterility of Crassostrea gigas

Simple Summary

This study reveals the expression profiles of lncRNA in the gonads of the Pacific oyster Crassostrea gigas. The potential function of lncRNAs was predicted in the case of antisense and cis-regulatory mechanisms based on their physical positions and their coexpression relationships in the case of trans regulation. Sterility-related DEGs and DELs were chosen for subsequent analysis, demonstrating that trans-regulatory lncRNAs might play a vital role in the gametogenesis of C. gigas. We constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites based on pathway results, in which 28 lncRNAs and their 54 trans-regulatory genes were detected. Among 28 sterility-specific lncRNAs, MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites play the most significant role.

Abstract

Long noncoding RNA (lncRNA), a type of non-protein-coding transcript, is emerging as a crucial regulator of gene expression. However, few roles of lncRNA in the reproductive process of the Pacific oyster (Crassostrea gigas) have been defined, especially in the regulatory mechanism of sterile triploids gametogenesis. To uncover the potential role of lncRNA, the gonads of diploids, sterile triploids, and partially sterile triploids underwent RNA sequencing. A total of 9618 reliable lncRNAs were identified. The target relationship between lncRNA and mRNA was predicted based on cis, trans, and antisense regulation with bioinformatic software. We chose differentially expressed lncRNAs and mRNAs when sterile triploids were compared to partially sterile triploids and diploids for subsequent functional enrichment analysis. Findings revealed that trans-regulatory lncRNAs might play a significant role in the gametogenesis of C. gigas. Combining pathway results, we constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites. Fifty-four genes related to cell division, germline-cell maintenance, and glycogen metabolism were found to be associated with sterility. A total of 28 candidate lncRNAs were predicted to trans-regulate these genes. We speculated that MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites were highly important as they were predicted to regulate more sterility-specific genes than others. Our work collectively identified sterility-related lncRNAs and implicated the potential mechanism of lncRNA-mediated regulation in the gametogenesis of sterile triploid oysters.

Overexpression of bmp4, dazl, nanos3 and sycp2 in Hu Sheep Leydig Cells Using CRISPR/dcas9 System Promoted Male Germ Cell Related Gene Expression

Simple Summary

Male germ cell development plays a crucial role in male reproduction, and gene expression also presents an essential regulatory role in its development. Many studies have been devoted to the induction and differentiation of pluripotent stem cells into germ cells in vitro. However, the culture system for pluripotent stem cells from domestic animals is not stable, especially in sheep. Our study attempted to transdifferentiate sheep somatic cells into germ cells in vitro by the overexpression of key germ cell related genes, with the aim of perfecting the construction of germ cell research models in vitro. Therefore, we explored the expression pattern of four crucial genes, bmp4, dazl, nanos3 and sycp2, in Hu sheep testicular development, and investigated the potential efficiency of overexpression of the four candidate genes using the CRISPR/dcas9 system in Leydig cells. We revealed that the overexpression of bmp4, dazl, nanos3 and sycp2 can promote the expression of male germ cell related genes. To the best of our knowledge, this is the first study to construct an overexpression induction system using CRISPR/dcas9 technology, and to induce sheep somatic cells into germ cells in vitro.

Abstract

Male germ cells directly affect the reproduction of males; however, their accurate isolation and culture in vitro is extremely challenging, hindering the study of germ cell development and function. CRISPR/dcas9, as an efficient gene reprogramming system, has been verified to promote the transdifferentiation of pluripotent stem cells into male germ cells by editing target genes. In our research, we explored the expression pattern of the germ cell related genes bmp4, dazl,nanos3 and sycp2 in Hu sheep testicular development and constructed the overexpression model using the CRISPR/dcas9 system. The results indicated that four genes showed more expression in testis tissue than in other tissues, and that bmp4, dazl and sycp2 present higher expression levels in nine-month-old sheep testes than in three-month-olds, while nanos3 expressed the opposite trend (p < 0.05). In addition, the expression of four potential genes in spermatogenic cells was slightly different, but they were all expressed in sheep Leydig cells. To verify the potential roles of the four genes in the process of inducing differentiation of male germ cells, we performed cell transfection in vitro. We found that the expression of the germ cell related genes Prdm1, Prdm14, Mvh and Sox17 were significantly increased after the overexpression of the four genes in Leydig cells, and the co-transfection effect was the most significant (p < 0.05). Our results illustrate the crucial functions of bmp4, dazl, nanos3 and sycp2 in Hu sheep testis development and verified the effectiveness of the overexpression model that was constructed using the CRISPR/dcas9 system, which provided a basis for further male germ cell differentiation in vitro.

Identification of Key Genes Related to Postnatal Testicular Development Based on Transcriptomic Data of Testis in Hu Sheep

The selection of testis size can improve the reproductive capacity of livestock used for artificial insemination and has been considered as an important strategy for accelerating the breeding process. Although much work has been done to investigate the mechanisms of testis development in various species, there is little information available in regard to the differences in transcriptomic profiling of sheep testes at different developmental stages. In this work, we aimed to identify differentially expressed genes (DEGs) by RNA-Seq in sheep during different growth stages, including 0 month old (infant, M0), 3 months old (puberty, M3), 6 months old (sexual maturity, M6) and 12 months old (body maturity, M12). A total of 4,606 (2,381 up and 2,225 down), 7,500 (4,368 up and 3,132 down), 15 (8 up and seven down) DEGs were identified in M3_vs_M0, M6_vs_M3, and M12_vs_M6 comparison, respectively. Of which, a number of genes were continuously up-regulated and down-regulated with testicular development, including ODF3, ZPBP1, PKDREJ, MYBL1, PDGFA, IGF1, LH, INSL3, VIM, AMH, INHBA, COL1A1, COL1A2, and INHA. GO analysis illustrated that DEGs were mainly involved in testis development and spermatogenesis. KEGG analysis identified several important pathways and verified several reproduction-associated DEGs such as COL1A1, COL1A2, PDGFA, and IGF1. In addition, two gene modules highly associated with testis development and core genes with testis size were identified using weighted gene co-expression network analysis (WGCNA), including hub genes positively associated with testis size such as RANBP9, DNAH17, SPATA4, CIB4 and SPEM1, and those negatively associated with testis size such as CD81, CSK, PDGFA, VIM, and INHBA. This study comprehensively identified key genes related to sheep testicular development, which may provide potential insights for understanding male fertility and better guide in animal breeding.

The Transcriptional Cell Atlas of Testis Development in Sheep at Pre-Sexual Maturity

Sheep testes undergo a dramatic rate of development with structural changes during pre-sexual maturity, including the proliferation and maturation of somatic niche cells and the initiation of spermatogenesis. To explore this complex process, 12,843 testicular cells from three males at pre-sexual maturity (three-month-old) were sequenced using the 10× Genomics Chromium^TM single-cell RNA-seq (scRNA-seq) technology. Nine testicular somatic cell types (Sertoli cells, myoid cells, monocytes, macrophages, Leydig cells, dendritic cells, endothelial cells, smooth muscle cells, and leukocytes) and an unknown cell cluster were observed. In particular, five male germ cell types (including two types of undifferentiated spermatogonia (A_pale and A_dark), primary spermatocytes, secondary spermatocytes, and sperm cells) were identified. Interestingly, A_pale and A_dark were found to be two distinct states of undifferentiated spermatogonia. Further analysis identified specific marker genes, including UCHL1, DDX4, SOHLH1, KITLG, and PCNA, in the germ cells at different states of differentiation. The study revealed significant changes in germline stem cells at pre-sexual maturation, paving the way to explore the candidate factors and pathways for the regulation of germ and somatic cells, and to provide us with opportunities for the establishment of livestock stem cell breeding programs.

Single-cell RNA-sequencing reveals the dynamic process and novel markers in porcine spermatogenesis

Background

Spermatogenesis is the process by which male gametes are formed from spermatogonial stem cells and it is essential for the reliable transmission of genetic information between generations. To date, the dynamic transcriptional changes of defined populations of male germ cells in pigs have not been reported.

Results

To characterize the atlas of porcine spermatogenesis, we profiled the transcriptomes of ~ 16,966 testicular cells from a 150-day-old pig testis through single-cell RNA-sequencing (scRNA-seq). The scRNA-seq analysis identified spermatogonia, spermatocytes, spermatids and three somatic cell types in porcine testes. The functional enrichment analysis demonstrated that these cell types played diverse roles in porcine spermatogenesis. The accuracy of the defined porcine germ cell types was further validated by comparing the data from scRNA-seq with those from bulk RNA-seq. Since we delineated four distinct spermatogonial subsets, we further identified CD99 and PODXL2 as novel cell surface markers for undifferentiated and differentiating spermatogonia, respectively.

Conclusions

The present study has for the first time analyzed the transcriptome of male germ cells and somatic cells in porcine testes through scRNA-seq. Four subsets of spermatogonia were identified and two novel cell surface markers were discovered, which would be helpful for studies on spermatogonial differentiation in pigs. The datasets offer valuable information on porcine spermatogenesis, and pave the way for identification of key molecular markers involved in development of male germ cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00638-3.

Single-Cell Atlas of Adult Testis in Protogynous Hermaphroditic Orange-Spotted Grouper, Epinephelus coioides

Spermatogenesis is a process of self-renewal and differentiation in spermatogonial stem cells. During this process, germ cells and somatic cells interact intricately to ensure long-term fertility and accurate genome propagation. Spermatogenesis has been intensely investigated in mammals but remains poorly understood with regard to teleosts. Here, we performed single-cell RNA sequencing of ~9500 testicular cells from the male, orange-spotted grouper. In the adult testis, we divided the cells into nine clusters and defined ten cell types, as compared with human testis data, including cell populations with characteristics of male germ cells and somatic cells, each of which expressed specific marker genes. We also identified and profiled the expression patterns of four marker genes (calr, eef1a, s100a1, vasa) in both the ovary and adult testis. Our data provide a blueprint of male germ cells and supporting somatic cells. Moreover, the cell markers are candidates that could be used for further cell identification.