MvMRL: a multi-view molecular representation learning method for molecular property prediction.

Effective molecular representation learning is very important for Artificial Intelligence-driven Drug Design because it affects the accuracy and efficiency of molecular property prediction and other molecular modeling relevant tasks. However, previous molecular representation learning studies often suffer from limitations, such as over-reliance on a single molecular representation, failure to fully capture both local and global information in molecular structure, and ineffective integration of multiscale features from different molecular representations. These limitations restrict the complete and accurate representation of molecular structure and properties, ultimately impacting the accuracy of predicting molecular properties. To this end, we propose a novel multi-view molecular representation learning method called MvMRL, which can incorporate feature information from multiple molecular representations and capture both local and global information from different views well, thus improving molecular property prediction. Specifically, MvMRL consists of four parts: a multiscale CNN-SE Simplified Molecular Input Line Entry System (SMILES) learning component and a multiscale Graph Neural Network encoder to extract local feature information and global feature information from the SMILES view and the molecular graph view, respectively; a Multi-Layer Perceptron network to capture complex non-linear relationship features from the molecular fingerprint view; and a dual cross-attention component to fuse feature information on the multi-views deeply for predicting molecular properties. We evaluate the performance of MvMRL on 11 benchmark datasets, and experimental results show that MvMRL outperforms state-of-the-art methods, indicating its rationality and effectiveness in molecular property prediction. The source code of MvMRL was released in https://github.com/jedison-github/MvMRL.

Rapid conversion of porcine pluripotent stem cells into macrophages with chemically defined conditions.

A renewable source of porcine macrophages derived from pluripotent stem cells (PSCs) would be a valuable alternative to primary porcine alveolar macrophages (PAMs) in the research of host-pathogen interaction mechanisms. We developed an efficient and rapid protocol, within 11 days, to derive macrophages from porcine PSCs (pPSCs). The pPSC-derived macrophages (pPSCdMs) exhibited molecular and functional characteristics of primary macrophages. The pPSCdMs showed macrophage-specific surface protein expression and macrophage-specific transcription factors, similar to PAMs. The pPSCdMs also exhibited the functional characteristics of macrophages, such as endocytosis, phagocytosis, porcine respiratory and reproductive syndrome virus infection and the response to lipopolysaccharide stimulation. Furthermore, we performed transcriptome sequencing of the whole differentiation process to track the fate transitions of porcine PSCs involved in the signaling pathway. The activation of transforming growth factor beta signaling was required for the formation of mesoderm and the inhibition of the transforming growth factor beta signaling pathway at the hematopoietic endothelium stage could enhance the fate transformation of hematopoiesis. In summary, we developed an efficient and rapid protocol to generate pPSCdMs that showed aspects of functional maturity comparable with PAMs. pPSCdMs could provide a broad prospect for the platforms of host-pathogen interaction mechanisms.

Super-enhancers conserved within placental mammals maintain stem cell pluripotency.

Despite pluripotent stem cells sharing key transcription factors, their maintenance involves distinct genetic inputs. Emerging evidence suggests that super-enhancers (SEs) can function as master regulatory hubs to control cell identity and pluripotency in humans and mice. However, whether pluripotency-associated SEs share an evolutionary origin in mammals remains elusive. Here, we performed comprehensive comparative epigenomic and transcription factor binding analyses among pigs, humans, and mice to identify pluripotency-associated SEs. Like typical enhancers, SEs displayed rapid evolution in mammals. We showed that BRD4 is an essential and conserved activator for mammalian pluripotency-associated SEs. Comparative motif enrichment analysis revealed 30 shared transcription factor binding motifs among the three species. The majority of transcriptional factors that bind to identified motifs are known regulators associated with pluripotency. Further, we discovered three pluripotency-associated SEs (SE-SOX2, SE-PIM1, and SE-FGFR1) that displayed remarkable conservation in placental mammals and were sufficient to drive reporter gene expression in a pluripotency-dependent manner. Disruption of these conserved SEs through the CRISPR-Cas9 approach severely impaired stem cell pluripotency. Our study provides insights into the understanding of conserved regulatory mechanisms underlying the maintenance of pluripotency as well as species-specific modulation of the pluripotency-associated regulatory networks in mammals.

Identification of different myofiber types in pigs muscles and construction of regulatory networks.

BACKGROUND: Skeletal muscle is composed of muscle fibers with different physiological characteristics, which plays an important role in regulating skeletal muscle metabolism, movement and body homeostasis. The type of skeletal muscle fiber directly affects meat quality. However, the transcriptome and gene interactions between different types of muscle fibers are not well understood. RESULTS: In this paper, we selected 180-days-old Large White pigs and found that longissimus dorsi (LD) muscle was dominated by fast-fermenting myofibrils and soleus (SOL) muscle was dominated by slow-oxidizing myofibrils by frozen sections and related mRNA and protein assays. Here, we selected LD muscle and SOL muscle for transcriptomic sequencing, and identified 312 differentially expressed mRNA (DEmRs), 30 differentially expressed miRNA (DEmiRs), 183 differentially expressed lncRNA (DElRs), and 3417 differentially expressed circRNA (DEcRs). The ceRNA network included ssc-miR-378, ssc-miR-378b-3p, ssc-miR-24-3p, XR_308817, XR_308823, SMIM8, MAVS and FOS as multiple core nodes that play important roles in muscle development. Moreover, we found that different members of the miR-10 family expressed differently in oxidized and glycolytic muscle fibers, among which miR-10a-5p was highly expressed in glycolytic muscle fibers (LD) and could target MYBPH gene mRNA. Therefore, we speculate that miR-10a-5p may be involved in the transformation of muscle fiber types by targeting the MYHBP gene. In addition, PPI analysis of differentially expressed mRNA genes showed that ACTC1, ACTG2 and ACTN2 gene had the highest node degree, suggesting that this gene may play a key role in the regulatory network of muscle fiber type determination. CONCLUSIONS: We can conclude that these genes play a key role in regulating muscle fiber type transformation. Our study provides transcriptomic profiles and ceRNA interaction networks for different muscle fiber types in pigs, providing reference for the transformation of pig muscle fiber types and the improvement of meat quality.

Integrated analysis of the ubiquitination mechanism reveals the specific signatures of tissue and cancer.

BACKGROUND: Ubiquitination controls almost all cellular processes. The dysregulation of ubiquitination signals is closely associated with the initiation and progression of multiple diseases. However, there is little comprehensive research on the interaction and potential function of ubiquitination regulators (UBRs) in spermatogenesis and cancer. METHODS: We systematically characterized the mRNA and protein expression of UBRs across tissues and further evaluated their roles in testicular development and spermatogenesis. Subsequently, we explored the genetic alterations, expression perturbations, cancer hallmark-related pathways, and clinical relevance of UBRs in pan-cancer. RESULTS: This work reveals heterogeneity in the expression patterns of UBRs across tissues, and the expression pattern in testis is the most distinct. UBRs are dynamically expressed during testis development, which are critical for normal spermatogenesis. Furthermore, UBRs have widespread genetic alterations and expression perturbations in pan-cancer. The expression of 79 UBRs was identified to be closely correlated with the activity of 32 cancer hallmark-related pathways, and ten hub genes were screened for further clinical relevance analysis by a network-based method. More than 90% of UBRs can affect the survival of cancer patients, and hub genes have an excellent prognostic classification for specific cancer types. CONCLUSIONS: Our study provides a comprehensive analysis of UBRs in spermatogenesis and pan-cancer, which can build a foundation for understanding male infertility and developing cancer drugs in the aspect of ubiquitination.

FifBase: a comprehensive fertility-associated indicators factor database for domestic animals.

Fertility refers to the ability of animals to maintain reproductive function and give birth to offspring, which is an important indicator to measure the productivity of animals. Fertility is affected by many factors, among which environmental factors may also play key roles. During the past years, substantial research studies have been conducted to detect the factors related to fecundity, including genetic factors and environmental factors. However, the identified genes associated with fertility from countless previous studies are randomly dispersed in the literature, whereas some other novel fertility-related genes are needed to detect from omics-based datasets. Here, we constructed a fertility index factor database FifBase based on manually curated published literature and RNA-Seq datasets. During the construction of the literature group, we obtained 3301 articles related to fecundity for 13 species from PubMed, involving 2823 genes, which are related to 75 fecundity indicators or 47 environmental factors. Eventually, 1558 genes associated with fertility were filtered in 10 species, of which 1088 and 470 were from RNA-Seq datasets and text mining data, respectively, involving 2910 fertility-gene pairs and 58 fertility-environmental factors. All these data were cataloged into FifBase (http://www.nwsuaflmz.com/FifBase/), where the fertility-related factor information, including gene annotation and environmental factors, can be browsed, retrieved and downloaded with the user-friendly interface.

Comprehensive comparison of two types of algorithm for circRNA detection from short-read RNA-Seq.

MOTIVATION: Circular RNA is generally formed by the 'back-splicing' process between the upstream splice acceptor and the downstream donor in/not in the regulation of the corresponding RNA-binding proteins or cis-elements. Therefore, more and more software packages have been developed and they are mostly based on the identification of the back-spliced junction reads. However, recent studies developed two software tools that can detect circRNA candidates by constructing k-mer table or/and de Bruijn graph rather than reads mapping. RESULTS: Here, we compared the precision, sensitivity and detection efficiency between software tools based on different algorithms. Eleven representative detection tools with two types of algorithm were selected for the overall pipeline analysis of RNA-seq datasets with/without RNase R treatment in two cell lines. Precision, sensitivity, AUC, F1 score and detection efficiency metrics were assessed to compare the prediction tools. Meanwhile, the sensitivity and distribution of highly expressed circRNAs before and after RNase R treatment were also revealed by their enrichment, unaffected and depleted candidate frequencies. Eventually, we found that compared to the k-mer based tools, CIRI2 and KNIFE based on reads mapping had relatively superior and more balanced detection performance regardless of the cell line or RNase R (-/+) datasets. AVAILABILITY AND IMPLEMENTATION: All predicted results and source codes can be retrieved from https://github.com/luffy563/circRNA_tools_comparison. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Histone demethylase complexes KDM3A and KDM3B cooperate with OCT4/SOX2 to define a pluripotency gene regulatory network.

The pluripotency gene regulatory network of porcine induced pluripotent stem cells(piPSCs), especially in epigenetics, remains elusive. To determine the biological function of epigenetics, we cultured piPSCs in different culture conditions. We found that activation of pluripotent gene- and pluripotency-related pathways requires the erasure of H3K9 methylation modification which was further influenced by mouse embryonic fibroblast (MEF) served feeder. By dissecting the dynamic change of H3K9 methylation during loss of pluripotency, we demonstrated that the H3K9 demethylases KDM3A and KDM3B regulated global H3K9me2/me3 level and that their co-depletion led to the collapse of the pluripotency gene regulatory network. Immunoprecipitation-mass spectrometry (IP-MS) provided evidence that KDM3A and KDM3B formed a complex to perform H3K9 demethylation. The genome-wide regulation analysis revealed that OCT4 (O) and SOX2 (S), the core pluripotency transcriptional activators, maintained the pluripotent state of piPSCs depending on the H3K9 hypomethylation. Further investigation revealed that O/S cooperating with histone demethylase complex containing KDM3A and KDM3B promoted pluripotency genes expression to maintain the pluripotent state of piPSCs. Together, these data offer a unique insight into the epigenetic pluripotency network of piPSCs.

Early cleavage of preimplantation embryos is regulated by tRNAGln-TTG-derived small RNAs present in mature spermatozoa.

tRNA-derived small RNAs (tsRNAs) from spermatozoa could act as acquired epigenetic factors and contribute to offspring phenotypes. However, the roles of specific tsRNAs in early embryo development remain to be elucidated. Here, using pigs as a research model, we probed the tsRNA dynamics during spermatogenesis and sperm maturation and demonstrated the delivery of tsRNAs from semen-derived exosomes to spermatozoa. By microinjection of antisense sequences into in vitro fertilized oocytes and subsequent single-cell RNA-seq of embryos, we identified a specific functional tsRNA group (termed here Gln-TTGs) that participate in the early cleavage of porcine preimplantation embryos, probably by regulating cell cycle-associated genes and retrotransposons. We conclude that specific tsRNAs present in mature spermatozoa play significant roles in preimplantation embryo development.

PAX7 promotes CD49f-positive dairy goat spermatogonial stem cells' self-renewal.

Spermatogenesis is a complex process that originates from and depends on the spermatogonial stem cells (SSCs). The number of SSCs is rare, which makes the separation and enrichment of SSCs difficult and inefficient. The transcription factor PAX7 maintains fertility in normal spermatogenesis in mice. However, for large animals, much less is known about the SSCs' self-renewal regulation, especially in dairy goats. We isolated and enriched the CD49f-positive and negative dairy goat testicular cells by magnetic-activated cell sorting strategies. The RNA- sequencing and experimental data revealed that cells with a high CD49f and PAX7 expression are undifferentiated spermatogonia in goat testis. Our findings indicated that ZBTB16 (PLZF), PAX7, LIN28A, BMPR1B, FGFR1, and FOXO1 were expressed higher in CD49f-positive cells as compared to negative cells and goat fibroblasts cells. The expression and distribution of PAX7 in dairy goat also have been detected, which gradually decreased in testis tissue along with the increasing age. When the PAX7 gene was overexpressed in dairy goat immortal mGSCs-I-SB germ cell lines, the expression of PLZF, GFRα1, ID4, and OCT4 was upregulated. Together, our data demonstrated that there is a subset of spermatogonial stem cells with a high expression of PAX7 among the CD49f+ spermatogonia, and PAX7 can maintain the self-renewal of CD49f-positive SSCs.

The ACE2 expression in Sertoli cells and germ cells may cause male reproductive disorder after SARS-CoV-2 infection.

The serious coronavirus disease-2019 (COVID-19) was first reported in December 2019 in Wuhan, China. COVID-19 is an infectious disease caused by severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2). Angiotensin converting enzyme 2(ACE2) is the cellular receptor for SARS-CoV-2. Considering the critical roles of testicular cells for the transmission of genetic information between generations, we analyzed single-cell RNA-sequencing (scRNA-seq) data of adult human testis. The mRNA expression of ACE2 was expressed in both germ cells and somatic cells. Moreover, the positive rate of ACE2 in testes of infertile men was higher than normal, which indicates that SARS-CoV-2 may cause reproductive disorders through pathway activated by ACE2 and the men with reproductive disorder may easily to be infected by SARS-CoV-2. The expression level of ACE2 was related to the age, and the mid-aged with higher positive rate than young men testicular cells. Taken together, this research provides a biological background of the potential route for infection of SARS-CoV-2 and may enable rapid deciphering male-related reproductive disorders induced by COVID-19.

CellSim: a novel software to calculate cell similarity and identify their co-regulation networks.

BACKGROUND: Cell direct reprogramming technology has been rapidly developed with its low risk of tumor risk and avoidance of ethical issues caused by stem cells, but it is still limited to specific cell types. Direct reprogramming from an original cell to target cell type needs the cell similarity and cell specific regulatory network. The position and function of cells in vivo, can provide some hints about the cell similarity. However, it still needs further clarification based on molecular level studies. RESULT: CellSim is therefore developed to offer a solution for cell similarity calculation and a tool of bioinformatics for researchers. CellSim is a novel tool for the similarity calculation of different cells based on cell ontology and molecular networks in over 2000 different human cell types and presents sharing regulation networks of part cells. CellSim can also calculate cell types by entering a list of genes, including more than 250 human normal tissue specific cell types and 130 cancer cell types. The results are shown in both tables and spider charts which can be preserved easily and freely. CONCLUSION: CellSim aims to provide a computational strategy for cell similarity and the identification of distinct cell types. Stable CellSim releases (Windows, Linux, and Mac OS/X) are available at: www.cellsim.nwsuaflmz.com , and source code is available at: https://github.com/lileijie1992/CellSim/ .

GED: a manually curated comprehensive resource for epigenetic modification of gametogenesis.

Reproductive infertility affects seventh of couples, which is most attributed to the obstacle of gametogenesis. Characterizing the epigenetic modification factors involved in gametogenesis is fundamental to understand the molecular mechanisms and to develop treatments for human infertility. Although the genetic factors have been implicated in gametogenesis, no dedicated bioinformatics resource for gametogenesis is available. To elucidate the relationship of epigenetic modification and mammalian gametogenesis, we developed a new database, gametogenesis epigenetic modification database (GED), a manually curated database, which aims at providing a comprehensive resource of epigenetic modification of gametogenesis. The database integrates three kinds information of epigenetic modifications during gametogenesis (DNA methylation, histone modification and RNA regulation), and the gametogenesis has been detailed as 16 stages in seven mammal species (Homo sapiens, Mus musculus, Rattus norvegicus, Sus scrofa, Bos taurus, Capra hircus and Ovis aries). Besides, we have predicted the linear pathways of epigenetic modification which were composed of 211 genes/proteins and microRNAs that were involved in gametogenesis. GED is a user-friendly Web site, through which users can obtain the comprehensive epigenetic factor information and molecular pathways by visiting our database freely. GED is free available at http://gametsepi.nwsuaflmz.com.

Dynamic and modular gene regulatory networks drive the development of gametogenesis.

Gametogenesis is a complex process, which includes mitosis and meiosis and results in the production of ovum and sperm. The development of gametogenesis is dynamic and needs many different genes to work synergistically, but it is lack of global perspective research about this process. In this study, we detected the dynamic process of gametogenesis from the perspective of systems biology based on protein-protein interaction networks (PPINs) and functional analysis. Results showed that gametogenesis genes have strong synergistic effects in PPINs within and between different phases during the development. Addition to the synergistic effects on molecular networks, gametogenesis genes showed functional consistency within and between different phases, which provides the further evidence about the dynamic process during the development of gametogenesis. At last, we detected and provided the core molecular modules of different phases about gametogenesis. The gametogenesis genes and related modules can be obtained from our Web site Gametogenesis Molecule Online (GMO, http://gametsonline.nwsuaflmz.com/index.php), which is freely accessible. GMO may be helpful for the reference and application of these genes and modules in the future identification of key genes about gametogenesis. Summary, this work provided a computational perspective and frame to the analysis of the gametogenesis dynamics and modularity in both human and mouse.

EWAS: epigenome-wide association study software 2.0.

Motivation: With the development of biotechnology, DNA methylation data showed exponential growth. Epigenome-wide association study (EWAS) provide a systematic approach to uncovering epigenetic variants underlying common diseases/phenotypes. But the EWAS software has lagged behind compared with genome-wide association study (GWAS). To meet the requirements of users, we developed a convenient and useful software, EWAS2.0. Results: EWAS2.0 can analyze EWAS data and identify the association between epigenetic variations and disease/phenotype. On the basis of EWAS1.0, we have added more distinctive features. EWAS2.0 software was developed based on our 'population epigenetic framework' and can perform: (i) epigenome-wide single marker association study; (ii) epigenome-wide methylation haplotype (meplotype) association study and (iii) epigenome-wide association meta-analysis. Users can use EWAS2.0 to execute chi-square test, t-test, linear regression analysis, logistic regression analysis, identify the association between epi-alleles, identify the methylation disequilibrium (MD) blocks, calculate the MD coefficient, the frequency of meplotype and Pearson's correlation coefficients and carry out meta-analysis and so on. Finally, we expect EWAS2.0 to become a popular software and be widely used in epigenome-wide associated studies in the future. Availability and implementation: The EWAS software is freely available at http://www.ewas.org.cn or http://www.bioapp.org/ewas.

miRNA editing landscape reveals miR-34c regulated spermatogenesis through structure and target change in pig and mouse.

Spermatogenesis has a close relationship with male infertility. MicroRNAs (miRNAs) play crucial roles in their regulation of target genes during spermatogenesis. A huge dataset of high-throughput sequencing all over the world provides the basis to dig the cryptic molecular mechanism. But how to take advantage of the big data and unearth the miRNA regulation is still a challenging problem. Here we integrated transcriptome of spermatogenesis and found miRNA regulate spermatogenesis through miRNA editing. We then compared different species and found that the distributions of miRNA editing site number and editing types among different cell types during spermatogenesis are conservative. Interesting, we further found that nearly half of the editing events occurred in the seed region in both mouse and pig. Finally, we foundmiR-34c, which is edited frequently at all stages during spermatogenesis, regulates its target genes through the RNA structure changing and shows dysfunction when it is edited. Summary, we depicted the overall profile of miRNA editing during spermatogenesis in mouse and pig and reveal miR-34c may play its roles through miRNA editing.

miR-544 Regulates Dairy Goat Male Germline Stem Cell Self-Renewal via Targeting PLZF.

The balance between the self-renewal and differentiation of male germline stem cells (mGSCs) is critical for the initiation and maintenance of mammalian spermatogenesis. The promyelocytic leukemia zinc finger (PLZF), a zinc finger protein, is a critical factor for maintaining the self-renewal of mGSCs, so, evaluation of the PLZF pathway in mGSCs may provide a deeper insight into mammalian spermatogenesis. miRNA was also an important regulating factor for the self-renewal and differentiation of mGSCs; however, there is currently no data indicating that which miRNA regulate the self-renewal and differentiation of mGSCs via PLZF. Here, we predicted the prospective miRNA targeting to PLZF using the online Bioinformatics database-Targetscan, and performed an analysis of the dual-luciferase recombinant vector, psiCHCEKTM-2-PLZF-3'UTR. miR-544 mimics (miR-544m), miR-544 inhibitors (miR-544i), Control (NC, scrambled oligonucleotides transfection), pPLZF-IRES2-EGFP or PLZF siRNA were transfected into mGSCs; the cells proliferation was evaluated by BRDU incorporation assay and flow cytometry, and the mGSC marker, GFRa1, PLZF, KIT, DAZL, and VASA expression were analyzed by RT-qPCR, immunofluorescence and Western blot. The results showed that miR-544 regulates dairy goat male germline stem cell self-renewal via targeting PLZF. Our study identifies a new regulatory pathway for PLZF and expands upon the PLZF regulatory network in mGSCs.

Pathway analysis of genome-wide association study and transcriptome data highlights new biological pathways in colorectal cancer.

Colorectal cancer (CRC) is a common malignancy that meets the definition of a complex disease. Genome-wide association study (GWAS) has identified several loci of weak predictive value in CRC, however, these do not fully explain the occurrence risk. Recently, gene set analysis has allowed enhanced interpretation of GWAS data in CRC, identifying a number of metabolic pathways as important for disease pathogenesis. Whether there are other important pathways involved in CRC, however, remains unclear. We present a systems analysis of KEGG pathways in CRC using (1) a human CRC GWAS dataset and (2) a human whole transcriptome CRC case-control expression dataset. Analysis of the GWAS dataset revealed significantly enriched KEGG pathways related to metabolism, immune system and diseases, cellular processes, environmental information processing, genetic information processing, and neurodegenerative diseases. Altered gene expression was confirmed in these pathways using the transcriptome dataset. Taken together, these findings not only confirm previous work in this area, but also highlight new biological pathways whose deregulation is critical for CRC. These results contribute to our understanding of disease-causing mechanisms and will prove useful for future genetic and functional studies in CRC.

CLU rs2279590 polymorphism contributes to Alzheimer's disease susceptibility in Caucasian and Asian populations.

It is reported that CLU rs2279590 polymorphism is significantly associated with Alzheimer's disease (AD) in European ancestry. Recent studies investigated rs2279590 polymorphism in Asian population (Chinese, Japanese and Korean). Four studies showed negative association and two studies showed weak association between rs2279590 and AD. We believe that the weak association or no association may be caused by the relatively small sample size in Asian population. Here, we reinvestigated the association in Asian population. Meanwhile, to investigate the genetic heterogeneity of the rs2279590 polymorphism in Asian and Caucasian populations, we searched the PubMed and AlzGene databases and selected 11 independent studies (6 studies in Asian population and 5 studies in Caucasian population) including 20,655 individuals (8,605 cases and 12,050 controls) for meta-analysis. Our results showed significant association between rs2279590 polymorphism and AD in Asian population with P = 2.00E-04 and P = 2.00E-04 using additive and recessive models, respectively. We observed no significant heterogeneity between Asian and Caucasian populations. We believe that our results may be helpful to understand the mechanisms of CLU in AD pathogenesis and will be useful for future genetic studies in AD.

CD49f-positive testicular cells in Saanen dairy goat were identified as spermatogonia-like cells by miRNA profiling analysis.

miRNAs, a type of small RNA, play critical roles in mammalian spermatogenesis. Spermatogonia are the foundation of spermatogenesis and are valuable for the study of spermatogenesis. However, the expression profiling of the miRNAs in spermatogonia of dairy goats remains unclear. CD49f has been one of the surface markers used for spermatogonia enrichment by magnetic activated cell sorting (MACS). Therefore, we used a CD49f microbead antibody to purify CD49f-positive and -negative cells of dairy goat testicular cells by MACS and then analysed the miRNA expression in these cells in depth using Illumina sequencing technology. The results of miRNA expression profiling in purified CD49f-positive and -negative testicular cells showed that 933 miRNAs were upregulated in CD49f-positive cells and 916 miRNAs were upregulated in CD49f-negative cells with a twofold increase, respectively; several miRNAs and marker genes specific for spermatogonial stem cells (SSCs) in testis had a higher expression level in CD49f-positive testicular cells, including miR-221, miR-23a, miR-29b, miR-24, miR-29a, miR-199b, miR-199a, miR-27a, and miR-21 and CD90, Gfra1, and Plzf. The bioinformatics analysis of differently expressed miRNAs indicated that the target genes of these miRNAs in CD49f-positive cells were involved in cell-cycle biological processes and the cell-cycle KEGG pathway. In conclusion, our comparative miRNAome data provide useful miRNA profiling data of dairy goat spermatogonia cells and suggest that CD49f could be used to enrich dairy goat spermatogonia-like cells, including SSCs.