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

Isoform-resolution single-cell RNA sequencing reveals the transcriptional panorama of adult Baoshan pig testis cells

BACKGROUND

As the primary organ of the male reproductive system, the testis facilitates spermatogenesis and androgen secretion. Due to the complexity of spermatogenesis, elucidating cellular heterogeneity and gene expression dynamics within the porcine testis is critical for advancing reproductive biology. Nevertheless, the cellular composition and regulatory mechanisms of porcine testes remain insufficiently characterized. In this study, we applied integrated long-read (Nanopore) and short-read (Illumina) scRNA-seq to Baoshan pig testes, establishing a comprehensive transcriptional profile to delineate cellular heterogeneity and molecular regulation.

RESULTS

Through systematic analysis of testicular architecture and the temporal progression of spermatogenesis, we characterized 11,520 single cells and 23,402 genes, delineating germ cell developmental stages: proliferative-phase spermatogonia (SPG), early-stage spermatocytes (Early SPC) and late-stage spermatocytes (Late SPC) during meiosis, and spermiogenic-phase round spermatids (RS) followed by elongating/elongated spermatids (ES), culminating in mature spermatozoa (Sperm). We further identified nine distinct testicular cell types, with germ cells spanning all developmental stages and somatic components comprising Sertoli cells, macrophages, and peritubular myoid cells as microenvironmental constituents, revealing the cellular heterogeneity of testicular tissue and dynamic characteristics of spermatogenesis. We obtained the dynamic expression changes of 16 vital marker genes during spermatogenesis and performed immunofluorescence validation on 7 marker genes. Gene ontology analysis revealed that germ cells at various stages were involved in specific biological processes, while cell communication networks highlighted eight pivotal signaling pathways, including MIF, NRG, WNT, VEGF, BMP, CCL, PARs, and ENHO pathways. Long-read sequencing further captured the full integrity and diversity of RNA transcripts, identifying 60% of the novel annotated isoforms and revealing that FSM isoforms exhibited longer transcript lengths, longer coding sequences, longer open reading frames, and a great number of exons, suggesting the complexity of isoforms within the testicular microenvironment.

CONCLUSIONS

Our results provide insight into the cellular heterogeneity, intercellular communication, and gene expression/transcript diversity in porcine testes, and offer a valuable resource for understanding the molecular mechanisms of porcine spermatogenesis.

Single-Cell RNA Sequencing Reveals an Atlas of Meihua Pig Testis Cells

Mammalian spermatogenesis is a complex biological process that is regulated by multiple types of cells. The heterogeneity of these cells poses a challenge for analyzing different cell types at different developmental stages. To characterize the transcriptomic landscape of porcine spermatogenesis and identify potential marker genes for spermatogonia, an unbiased transcriptomic study of spermatogenesis in neonatal and sexually mature six-month-old Meihua pigs was performed using 10× Genomics single-cell RNA sequencing (scRNA-seq). Through the collection of scRNA-seq data from 13,839 cells from Meihua pig testes, three germ cells (spermatogonia, spermatocytes and spermatids) and eight somatic cells (Sertoli cells, Leydig cells, myoid/stromal cells, endothelial cells, T cells/macrophages and erythroblasts) were identified. Pseudo-timing analysis showed that myoid cells and stromal cells originated from common progenitors in Meihua pigs. Functional enrichment analysis revealed that the differentially expressed genes (DEGs) in testicular somatic cells were enriched in the pathways of Ribosome, Oxidative phosphorylation, Protein processing in endoplasmic reticulum, Retrograde endocannabinoid signaling, Cellular senescence and Insulin signaling. Meanwhile, in the three different germ cells, except for pathways which were the same as the first three pathways for somatic cells, DEGs were also enriched in the Spliceosome, Cell cycle, Autophagy and Mitophagy pathways. Furthermore, the candidate marker gene TKTL1 in spermatogonia was identified using immunohistochemistry and immunofluorescence. In conclusion, we collected transcription datasets and constructed single-cell developmental maps of germ cells and somatic cells during the testicular development of Meihua pigs, which provided new insights into the spermatogenesis of Meihua pigs and the development of various types of cells in their testes.

A Review and New Insights Into Spermatogenesis From Single-Cell RNA Sequencing

As the most primitive germ cells in the testis, spermatogonial stem cells (SSCs) not only constantly renew themselves to ensure their quantity, but also differentiate into mature sperm cells to complete spermatogenesis and transmit genetic information to the next generation. Successful spermatogenesis is inseparable from niche regulation, which provides factors that enhance the self-renewal of SSCs to maintain their numbers and directs the appropriate differentiation of spermatogonia. Some progress has been achieved in the definition and isolation of SSCs. However, a high degree of cellular heterogeneity is found in the testis, revealing a combination of various cell types at different developmental stages and a lack of specific molecular markers (especially in domestic animals) for fully screening and purifying SSCs. These factors have considerably hindered further research into the mechanisms of maintenance, self-renewal, and differentiation of SSCs, as well as limited their isolation, purification, and applications. Accumulated studies have recently successfully employed single-cell RNA sequencing (scRNA-seq) as a novel approach to detailing the classification of cell subsets, mining specifically expressed genes in different cell types, and accurate identification of specific cell types. This review summarises the progress of SSCs identification and offers new insights into the SSCs developmental trajectory from single-cell RNA sequencing.

Genetic parameters and genome-wide association studies including the X chromosome for various reproduction and semen quality traits in Nellore cattle

BACKGROUND

The profitability of the beef industry is directly influenced by the fertility rate and reproductive performance of both males and females, which can be improved through selective breeding. When performing genomic analyses, genetic markers located on the X chromosome have been commonly ignored despite the X chromosome being one of the largest chromosomes in the cattle genome. Therefore, the primary objectives of this study were to: (1) estimate variance components and genetic parameters for eighteen male and five female fertility and reproductive traits in Nellore cattle including X chromosome markers in the analyses; and (2) perform genome-wide association studies and functional genomic analyses to better understand the genetic background of male and female fertility and reproductive performance traits in Nellore cattle.

RESULTS

The percentage of the total direct heritability (h2total) explained by the X chromosome markers (h2x) ranged from 3 to 32% (average: 16.4%) and from 9 to 67% (average: 25.61%) for female reproductive performance and male fertility traits, respectively. Among the traits related to breeding soundness evaluation, the overall bull and semen evaluation and semen quality traits accounted for the highest proportion of h2x relative to h2total with an average of 39.5% and 38.75%, respectively. The total number of significant genomic markers per trait ranged from 7 (seminal vesicle width) to 43 (total major defects). The number of significant markers located on the X chromosome ranged from zero to five. A total of 683, 252, 694, 382, 61, and 77 genes overlapped with the genomic regions identified for traits related to female reproductive performance, semen quality, semen morphology, semen defects, overall bulls' fertility evaluation, and overall semen evaluation traits, respectively. The key candidate genes located on the X chromosome are PRR32, STK26, TMSB4X, TLR7, PRPS2, SMS, SMARCA1, UTP14A, and BCORL1. The main gene ontology terms identified are "Oocyte Meiosis", "Progesterone Mediated Oocyte Maturation", "Thermogenesis", "Sperm Flagellum", and "Innate Immune Response".

CONCLUSIONS

Our findings indicate the key role of genes located on the X chromosome on the phenotypic variability of male and female reproduction and fertility traits in Nellore cattle. Breeding programs aiming to improve these traits should consider adding the information from X chromosome markers in their genomic analyses.

Transcriptional Profiling of Testis Development in Pre-Sexually-Mature Hezuo Pig

Spermatogenesis is an advanced biological process, relying on intricate interactions between somatic and germ cells in testes. Investigating various cell types is challenging because of cellular heterogeneity. Single-cell RNA sequencing (scRNA-seq) offers a method to analyze cellular heterogeneity. In this research, we performed 10× Genomics scRNA-seq to conduct an unbiased single-cell transcriptomic analysis in Hezuo pig (HZP) testis at one month of age during prepuberty. We collected 14,276 cells and identified 8 cell types (including 2 germ cells types and 6 somatic cell types). Pseudo-timing analysis demonstrated that Leydig cells (LCs) and myoid cells (MCs) originated from a shared progenitor cell lineage. Moreover, the functional enrichment analyses showed that the genes of differential expression were enriched in spermatogonia (SPG) and were enriched in the cell cycle, reproduction, and spermatogenesis. Expressed genes in spermatocytes (SPCs) were enriched in the cAMP, cell cycle, male gamete generation, reproductive system development, and sexual reproduction, while growth hormone synthesis, gamete generation, reproductive process, and spermine synthase activity were enriched in Sertoli cells (SCs). Additionally, chemokine, B cell receptor, activation of immune response, and enzyme binding were enriched in macrophages. Our study investigated transcriptional alterations across different cell types during spermatogenesis, yielding new understandings of spermatogenic processes and cell development. This research delivers an exploration of spermatogenesis and testicular cell biology in HZP, establishing the groundwork for upcoming breeding initiatives.

A multi-omics database of buffaloes from Yangtze valley reveals diversity of water buffalo (Bubalus bubalis)

Asian water buffalo (Bubalus bubalis) is the fundamental livestock resource for local rural populations and holds a promising prospect of their milk and meat. Xuyi mountain (XYM) and Haizi (HZ) buffaloes from Yangtze valley comprises of species diversity of Asian water buffaloes. Current multi-omics enables identification of causal genes and elucidation of genetic regulatory mechanisms underlying complex traits in buffaloes. Here, we conducted the integrated analysis of metabolome and metagenome of rumen fluid, transcriptome and metabolome of blood, and whole genome sequence data from XYM (n = 7) and HZ (n = 10) male buffaloes. Our results revealed the apparent diversity of multi-layer omics profiles between two buffalo species. The built-up multi-omics database supports the discoveries of diversity in Asian water buffalo and potentially serves valuable resources for studying causal regulatory variants and their mechanisms.

Innovative Insights into Single-Cell Technologies and Multi-Omics Integration in Livestock and Poultry

In recent years, single-cell RNA sequencing (scRNA-seq) has marked significant strides in livestock and poultry research, especially when integrated with multi-omics approaches. These advancements provide a nuanced view into complex regulatory networks and cellular dynamics. This review outlines the application of scRNA-seq in key species, including poultry, swine, and ruminants, with a focus on outcomes related to cellular heterogeneity, developmental biology, and reproductive mechanisms. We emphasize the synergistic power of combining scRNA-seq with epigenomic, proteomic, and spatial transcriptomic data, enhancing molecular breeding precision, optimizing health management strategies, and refining production traits in livestock and poultry. The integration of these technologies offers a multidimensional approach that not only broadens the scope of data analysis but also provides actionable insights for improving animal health and productivity.

Combined Use of Univariate and Multivariate Approaches to Detect Selection Signatures Associated with Milk or Meat Production in Cattle

OBJECTIVES

The aim of this study was to investigate the genomic structure of the cattle breeds selected for meat and milk production and to identify selection signatures between them.

METHODS

A total of 391 animals genotyped at 41,258 SNPs and belonging to nine breeds were considered: Angus (N = 62), Charolais (46), Hereford (31), Limousin (44), and Piedmontese (24), clustered in the Meat group, and Brown Swiss (42), Holstein (63), Jersey (49), and Montbéliarde (30), clustered in the Milk group. The population stratification was analyzed by principal component analysis (PCA), whereas selection signatures were identified by univariate (Wright fixation index, FST) and multivariate (canonical discriminant analysis, CDA) approaches. Markers with FST values larger than three standard deviations from the chromosomal mean were considered interesting. Attention was focused on markers selected by both techniques.

RESULTS

A total of 10 SNPs located on seven different chromosomes (7, 10, 14, 16, 17, 18, and 24) were identified. Close to these SNPs (±250 kb), 165 QTL and 51 genes were found. The QTL were grouped in 45 different terms, of which three were significant (Bonferroni correction < 0.05): milk fat content, tenderness score, and length of productive life. Moreover, genes mainly associated with milk production, immunity and environmental adaptation, and reproduction were mapped close to the common SNPs.

CONCLUSIONS

The results of the present study suggest that the combined use of univariate and multivariate approaches can help to better identify selection signatures due to directional selection.

Single-cell transcriptomes reveal spermatogonial stem cells and the dynamic heterogeneity of spermatogenesis in a seasonal breeding teleost

Seasonal spermatogenesis in fish is driven by spermatogonial stem cells (SSCs), which undergo a complex cellular process to differentiate into mature sperm. In this study, we characterized spermatogenesis in the large yellow croaker (Larimichthys crocea), a marine fish of significant commercial value, based on a high-resolution single-cell RNA-sequencing atlas of testicular cells from three distinct developmental stages: juvenile, adult differentiating and regressed testes. We detailed a continuous developmental trajectory of spermatogenic cells, from spermatogonia to spermatids, elucidating the molecular events involved in spermatogenesis. We uncovered dynamic heterogeneity in cellular compositions throughout the annual reproductive cycle, accompanied by strong molecular signatures within specific testicular cells. Notably, we identified a distinct population of SSCs and observed a critical metabolic transition from glycolysis to oxidative phosphorylation, enhancing our understanding of the biochemical and molecular characteristics of SSCs. Additionally, we elucidated the interactions between somatic cells and spermatogonia, illuminating the mechanisms that regulate SSC development. Overall, this work enhances our understanding of spermatogenesis in seasonal breeding teleosts and provides essential insights for the further conservation and culture of SSCs.

Germ cell progression through zebrafish spermatogenesis declines with age

Vertebrate spermatogonial stem cells maintain sperm production over the lifetime of an animal, but fertility declines with age. Although morphological studies have informed our understanding of typical spermatogenesis, the molecular and cellular mechanisms underlying the maintenance and decline of spermatogenesis are not yet understood. We used single-cell RNA sequencing to generate a developmental atlas of the aging zebrafish testis. All testes contained spermatogonia, but we observed a progressive decline in spermatogenesis that correlated with age. Testes from some older males only contained spermatogonia and a reduced population of spermatocytes. Spermatogonia in older males were transcriptionally distinct from spermatogonia in testes capable of robust spermatogenesis. Immune cells including macrophages and lymphocytes drastically increased in abundance in testes that could not complete spermatogenesis. Our developmental atlas reveals the cellular changes as the testis ages and defines a molecular roadmap for the regulation of spermatogenesis.

Single-Cell RNA Sequencing Reveals an Atlas of Hezuo Pig Testis Cells

Spermatogenesis is a complex biological process crucial for male reproduction and is characterized by intricate interactions between testicular somatic cells and germ cells. Due to the cellular heterogeneity of the testes, investigating different cell types across developmental stages has been challenging. Single-cell RNA sequencing (scRNA-seq) has emerged as a valuable approach for addressing this limitation. Here, we conducted an unbiased transcriptomic study of spermatogenesis in sexually mature 4-month-old Hezuo pigs using 10× Genomics-based scRNA-seq. A total of 16,082 cells were collected from Hezuo pig testes, including germ cells (spermatogonia (SPG), spermatocytes (SPCs), spermatids (SPTs), and sperm (SP)) and somatic cells (Sertoli cells (SCs), Leydig cells (LCs), myoid cells (MCs), endothelial cells (ECs), and natural killer (NK) cells/macrophages). Pseudo-time analysis revealed that LCs and MCs originated from common progenitors in the Hezuo pig. Functional enrichment analysis indicated that the differentially expressed genes (DEGs) in the different types of testicular germ cells were enriched in the PI3K-AKT, Wnt, HIF-1, and adherens junction signaling pathways, while the DEGs in testicular somatic cells were enriched in ECM-receptor interaction and antigen processing and presentation. Moreover, genes related to spermatogenesis, male gamete generation, sperm part, sperm flagellum, and peptide biosynthesis were expressed throughout spermatogenesis. Using immunohistochemistry, we verified several stage-specific marker genes (such as UCHL1, WT1, SOX9, and ACTA2) for SPG, SCs, and MCs. By exploring the changes in the transcription patterns of various cell types during spermatogenesis, our study provided novel insights into spermatogenesis and testicular cells in the Hezuo pig, thereby laying the foundation for the breeding and preservation of this breed.

Advances in single-cell transcriptomics in animal research

Understanding biological mechanisms is fundamental for improving animal production and health to meet the growing demand for high-quality protein. As an emerging biotechnology, single-cell transcriptomics has been gradually applied in diverse aspects of animal research, offering an effective method to study the gene expression of high-throughput single cells of different tissues/organs in animals. In an unprecedented manner, researchers have identified cell types/subtypes and their marker genes, inferred cellular fate trajectories, and revealed cell‒cell interactions in animals using single-cell transcriptomics. In this paper, we introduce the development of single-cell technology and review the processes, advancements, and applications of single-cell transcriptomics in animal research. We summarize recent efforts using single-cell transcriptomics to obtain a more profound understanding of animal nutrition and health, reproductive performance, genetics, and disease models in different livestock species. Moreover, the practical experience accumulated based on a large number of cases is highlighted to provide a reference for determining key factors (e.g., sample size, cell clustering, and cell type annotation) in single-cell transcriptomics analysis. We also discuss the limitations and outlook of single-cell transcriptomics in the current stage. This paper describes the comprehensive progress of single-cell transcriptomics in animal research, offering novel insights and sustainable advancements in agricultural productivity and animal health.

Single-Cell Transcriptomic Analysis Revealed the Cell Population Changes and Cell-Cell Communication in the Liver of a Carnivorous Fish in Response to High-Carbohydrate Diet

BACKGROUND

Carnivorous fish have a low carbohydrate utilization ability, and the physiologic and molecular basis of glucose intolerance has not been fully illustrated.

OBJECTIVES

This study aimed to use largemouth bass as a model to investigate the possible mechanism of glucose intolerance in carnivorous fish with the help of single-nuclei RNA sequencing (snRNA-seq).

METHODS

Two diets were formulated, a low-carbohydrate (LC) diet and a high-carbohydrate (HC) diet. The feeding trial lasted for 6 wk, and then, growth performance, biochemical parameters, liver histology, and snRNA-seq were performed.

RESULTS

Growth performance of fish was not affected by the HC diet, while liver glucolipid metabolism disorder and liver injury were observed. A total of 13,247 and 12,848 cells from the liver derived from 2 groups were isolated and sequenced, and 7 major liver cell types were annotated by the marker genes. Hepatocytes and cholangiocytes were lower and hepatic stellate cells (HSCs) and immune cells were higher in the HC group than those in the LC group. Reclustering analysis identified 7 subtypes of hepatocytes and immune cells, respectively. The HSCs showed more cell communication with other cell types, and periportal hepatocytes showed more cell communication with other hepatocyte subtypes. Cell-cell communication mainly focused on cell junction-related signaling pathways. Uncovered by the pseudotime analysis, midzonal hepatocytes were differentiated into 2 major branches-biliary epithelial hepatocytes and hepatobiliary hybrid progenitor. Cell junction and liver fibrosis-related genes were highly expressed in the HC group. HC diet induced the activation of HSCs and, therefore, led to the liver fibrosis of largemouth bass.

CONCLUSIONS

HC diet induces liver glucolipid metabolism disorder and liver injury of largemouth bass. The increase and activation of HSCs might be the main reason for the liver injury. In adaption to HC diet, midzonal hepatocytes differentiates into 2 major branches-biliary epithelial hepatocytes and hepatobiliary hybrid progenitors.

Applications of single-cell RNA sequencing in spermatogenesis and molecular evolution

Spermatogenic cell heterogeneity is determined by the complex process of spermatogenesis differentiation. However, effectively revealing the regulatory mechanisms underlying mammalian spermatogenic cell development and differentiation via traditional methods is difficult. Advances in technology have led to the emergence of many single-cell transcriptome sequencing protocols, which have partially addressed these challenges. In this review, we detail the principles of 10x Genomics technology and summarize the methods for downstream analysis of single-cell transcriptome sequencing data. Furthermore, we explore the role of single-cell transcriptome sequencing in revealing the heterogeneity of testicular ecological niche cells, delineating the establishment and disruption of testicular immune homeostasis during human spermatogenesis, investigating abnormal spermatogenesis in humans, and, ultimately, elucidating the molecular evolution of mammalian spermatogenesis.

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.

Single-Cell Transcriptome Sequencing Reveals Molecular Expression Differences and Marker Genes in Testes during the Sexual Maturation of Mongolian Horses

This study aimed to investigate differences in testicular tissue morphology, gene expression, and marker genes between sexually immature (1-year-old) and sexually mature (10-year-old) Mongolian horses. The purposes of our research were to provide insights into the reproductive physiology of male Mongolian horses and to identify potential markers for sexual maturity. The methods we applied included the transcriptomic profiling of testicular cells using single-cell sequencing techniques. Our results revealed significant differences in tissue morphology and gene expression patterns between the two age groups. Specifically, 25 cell clusters and 10 cell types were identified, including spermatogonial and somatic cells. Differential gene expression analysis highlighted distinct patterns related to cellular infrastructure in sexually immature horses and spermatogenesis in sexually mature horses. Marker genes specific to each stage were also identified, including APOA1, AMH, TAC3, INHA, SPARC, and SOX9 for the sexually immature stage, and PRM1, PRM2, LOC100051500, PRSS37, HMGB4, and H1-9 for the sexually mature stage. These findings contribute to a deeper understanding of testicular development and spermatogenesis in Mongolian horses and have potential applications in equine reproductive biology and breeding programs. In conclusion, this study provides valuable insights into the molecular mechanisms underlying sexual maturity in Mongolian horses.

Integrative analysis of single-nucleus RNA-seq and bulk RNA-seq reveals germline cells development dynamics and niches in the Pacific oyster gonad

Gametogenesis drives the maturation of germ cell precursors into functional gametes, facilitated by interactions with the niche environment. However, the molecular mechanisms, especially in invertebrates, remain incompletely understood. In this study, the gonadal microenvironment and gametogenic processes in the Pacific oyster, a model for diffuse gonadal organization and periodic gametogenesis, are investigated. We combine single-nucleus RNA-seq and bulk RNA-seq to analyze gonadal microenvironments in oysters. Twenty-three male and nineteen female gonadal cell clusters are identified, revealing four male and three female germ cell types, alongside follicular cells in females and Sertoli/Leydig cells in males. The NOTCH and BMP (bone morphogenetic protein) signaling pathways play a significant role in the male germline niche, suggesting similarities with mammalian germ cell microenvironment. This study offers valuable insights into germ cell developmental transitions and microenvironmental characteristics.

Revelation of genetic diversity and genomic footprints of adaptation in Indian pig breeds

In the present study, the genetic diversity measures among four Indian domestic breeds of pig namely Agonda Goan, Ghurrah, Ghungroo, and Nicobari, of different agro-climatic regions of country were explored and compared with European commercial breeds, European wild boar and Chinese domestic breeds. The double digest restriction site-associated DNA sequencing (ddRADseq) data of Indian pigs (102) and Landrace (10 animals) were generated and whole genome sequencing data of exotic pigs (60 animals) from public data repository were used in the study. The principal component analysis (PCA), admixture analysis and phylogenetic analysis revealed that Indian breeds were closer in ancestry to Chinese breeds than European breeds. European breeds exhibited highest genetic diversity measures among all the considered breeds. Among Indian breeds, Agonda Goan and Ghurrah were found to be more genetically diverse than Nicobari and Ghungroo. The selection signature regions in Indian pigs were explored using iHS and XP-EHH, and during iHS analysis, it was observed that genes related to growth, reproduction, health, meat quality, sensory perception and behavior were found to be under selection pressure in Indian pig breeds. Strong selection signatures were recorded in 24.25-25.25 Mb region of SSC18, 123.25-124 Mb region of SSC15 and 118.75-119.5 Mb region of SSC2 in most of the Indian breeds upon pairwise comparison with European commercial breeds using XP-EHH. These regions were harboring some important genes such as EPHA4 for thermotolerance, TAS2R16, FEZF1, CADPS2 and PTPRZ1 for adaptability to scavenging system of rearing, TRIM36 and PGGT1B for disease resistance and CCDC112, PIAS1, FEM1B and ITGA11 for reproduction.

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.

Germ cells do not progress through spermatogenesis in the infertile zebrafish testis

Vertebrate spermatogonial stem cells maintain sperm production over the lifetime of an animal but fertility declines with age. While morphological studies have greatly informed our understanding of typical spermatogenesis, the molecular and cellular mechanisms underlying spermatogenesis are not yet understood, particularly with respect to the onset of fertility. We used single-cell RNA sequencing to generate a developmental atlas of the zebrafish testis. Using 5 timepoints across the adult life of a zebrafish, we described cellular profiles in the testis during and after fertility. While all germ cell stages of spermatogenesis are detected in testes from fertile adult zebrafish, testes from older infertile males only contained spermatogonia and a reduced population of spermatocytes. These remaining germ cells are transcriptionally distinct from fertile spermatogonia. Immune cells including macrophages and lymphocytes drastically increase in abundance in infertile testes. Our developmental atlas reveals the cellular changes as the testis ages and defines a molecular roadmap for the regulation of male fertility.