Next Article in Journal
Net Conversion of Human-Edible Vitamins and Minerals in the U.S. Southern Great Plains Beef Production System
Next Article in Special Issue
YPEL3 Negatively Regulates Endometrial Function via the Wnt/β-Catenin Pathways during Early Pregnancy in Goats
Previous Article in Journal
Effects of Dietary Supplementation with Glutamine on the Immunity and Intestinal Barrier Gene Expression in Broiler Chickens Infected with Salmonella Enteritidis
Previous Article in Special Issue
CYP19A1 May Influence Lambing Traits in Goats by Regulating the Biological Function of Granulosa Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Inner Mongolia Key Laboratory of Basic Veterinary Science, Inner Mongolia Agriculture University, Hohhot 010018, China
2
Department of Psychosomatic Medicine, Inner Mongolia Medical University, Hohhot 010030, China
3
Inner Mongolia Academy of Agriculture & Animal Husbandry Sciences, Hohhot 010000, China
4
Inner Mongolia Saikexing Institutes of Breeding and Reproductive Biotechnologies in Domestic Animal, Hohhot 011517, China
5
Research Center for Animal Genetic Resources of Mongolia Plateau, College of Life Science, Inner Mongolia University, Hohhot 010021, China
*
Authors to whom correspondence should be addressed.
Animals 2022, 12(17), 2169; https://doi.org/10.3390/ani12172169
Submission received: 21 June 2022 / Revised: 9 August 2022 / Accepted: 17 August 2022 / Published: 24 August 2022
(This article belongs to the Collection Reproductive Management of Sheep and Goats)

Abstract

:

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.

1. Introduction

In mammals, sex is determined by the heteromorphic sex chromosomes: XY chromosomes in males and XX chromosomes in females. Different genomic doses of sex chromosomes lead to unbalanced gene expression. The somatic expression of X-linked genes must be adjusted so that males and females produce similar levels of most proteins encoded on the variable chromosome [1]. Susumo Ohno suggested that to rectify this imbalance, the expression of X-chromosome genes should be increased two-fold to match the output of the diploid autosomal complement, i.e., to enable an X-to-autosome ratio (X:A) of 1 [2]. This is known as ‘dosage compensation’. Dosage compensation is achieved by increasing the X-chromosome expression levels of dosage-sensitive X-linked genes in both sexes and through the random silencing of one X chromosome in females (X-chromosome inactivation, XCl) [3,4].
Dosage compensation in animal somatic tissues has been extensively studied [4,5,6]. Two previous microarray studies fully supported dosage compensation in both humans and mice [5,7]. However, an RNA sequencing (RNA-seq) study in humans and mice claimed that microarray-based expression was not suitable for the comparison of the expression levels of different genes and there was no evidence of dosage compensation. Subsequently, the same RNA-seq data were reanalyzed in humans and mice and dosage compensation was found after taking into account the skewed gene content of the X chromosome [8,9]. X-chromosome upregulation in sheep has been supported by results from RNA sequencing (RNA seq) [4]. During male meiosis in some animal species, the X and Y become transcriptionally silenced in a process known as meiotic sex chromosome inactivation (MSCI) [10]. Following the initiation of MSCl, chromosomes undergo various genetic modification. However, whether dosage compensation occurs in germ cells is actively debated [11]. For instance, research has revealed the two-fold upregulation of X-linked genes in the germline of adult Drosophila testes [11,12]. Other research has found that the expression of the X chromosome is greatly suppressed in the male germline [13]. Further studies have suggested there is little or no evidence for dosage compensation and X-chromosome-specific inactivation during meiosis [14].
Despite the fact that dosage compensation in mammalian cells is fundamentally different compared to Drosophila [15], at least the core MSL complex consisting of MSL1, 2, 3 and MOF is conserved in mammalian species. The MSL complex in humans consists of hMSL1, hMSL2, hMSL3 and hMOF (hMYST1), whereas in mice, Hampin (Homolog of drosophila MSL1), MSL2, MRG15 and MYST1 were detected [16,17,18,19,20,21]. Remarkably, there is also evidence for a function of MOF and/or the MSL complex in upregulating the active X chromosome [8,22]. The study of the expression of the MSL complex during sheep spermatogenesis will further be essential for clarifying the role of the MSL complex in regulating mammalian X inactivation as well as activation.
The mechanism of germ cell and somatic cell dosage compensation in sheep testis remains unknown. In this study, we used RNA-seq data to study germ cell and somatic cell dosage compensation in sheep testes in different developmental stages. Meanwhile, we perform single-cell RNA-seq analysis to examine the timing of X:A balancing in newborn, juvenile and adult sheep testes.

2. Materials and Methods

2.1. Animal Care

The Institutional Animal Care and Use Committee of Inner Mongolia Agricultural University approved the experimental protocol employed in this study.

2.2. Reagents

All the reagents were purchased from Sigma–Aldrich Co. (St. Louis, MO, USA), unless otherwise indicated.

2.3. Animals and Sample Collection

Three healthy male Hu sheep at the ages of 0 days, 90 days and 720 days were used in this study. Hu sheep fed at the Aowoteke Animal Husbandry Technology Co., Ltd., Chifeng, China, were selected and castrated for testis collection. After the castrations were performed, the testes were removed and washed several times using DPBS, immediately transferred to Tissue Storage Solution (130-100-008, Miltenyi, Miltenyi Biotec B.V & CO, Mönchengladbach, Germany) at 4 °C and digested via two-step enzyme digestion within 8 h [8].

2.4. Preparation of Single Testicular Cell Suspension

The testis tissues were washed four times with DPBS and cut into small pieces. Small pieces of testicular tissue were washed four times with DPBS and placed in a 1.5 mL centrifuge tube. Single testicular cells were isolated using enzymatic digestion. In brief, the tissues were incubated in dissociation buffer containing collagenase IV 10 mg/mL, Dispase II 250 ug/mL, DNase I 1:1000 and 10% FBS in Modified Eagle’s Medium (DMEM, Gibco Laboratories, Grand City, NY, USA) for 30 min at 37 °C. The cell resuspensions were washed and resuspended and then filtered with through 40 μm strainers. The cell viability was determined via trypan blue staining and the cell numbers were determined using a hemocytometer (Countstar, New York, NY, USA). A cell suspension at a concentration of ~1000 cells/µL and cell viability >80% was used for single-cell RNA sequencing.

2.5. Single-Cell RNA-Sequencing Performance, Library Preparation and Sequencing

Single-cell RNA sequencing was performed using the 10× Genomics system. The single cells were loaded on a Chromium Single Cell Controller (10× Genomics, Pleasanton, CA, USA) to generate single-cell gel beads in emulsion (GEMs) by using the Single Cell 30 Library and Gel Bead Kit V2 (10× Genomics, 120237). The lysis and barcoded reverse transcription of polyadenylated mRNA from single cells were performed inside each GEM and they were cleaned. cDNA was amplified, then fragmented and fragment end repaired. An Agilent 2100 bioanalyzer (Agilent Technologies Inc., Santa Clara, CA, USA) and quantitative real-time PCR (qRT-PCR) were used to quantify the library. Finally, the qualified cDNA libraries were sequenced using the Illumina Hiseq 2000 (Illumina, Inc., Hayward, CA, USA).

2.6. Mapping, Cell Identification and Clustering Analysis

After low-quality reads were removed (nFeature_RNA < 200 or nFeature_RNA > 6000), Cell Ranger software (Cell Ranger v6.1.1, 10× Genomics, Pleasanton, CA, USA) was used to carry out gene expression quantification and cell-type identification based on the cell barcode and UMI information of read 1 in every single cell. The corresponding read 2 was mapped to the reference genome (Ovis aries, Oar_rambouillet_v1.0) using STAR software to determine the source genes of reads that completed the gene expression quantification. We used Seurat software and the PCA to classify the cell types. The filtering criteria of effective cells, including the number of genes per cell and the number of reads per cell, were kept in a range of 1–99% and the mitochondrial percentage was less than 5%. For the top 10 principal components with the largest variance explained in the PCA results, UMAP (uniform manifold approximation and projection) was used to visualize the single-cell clustering. The Louvain–Jaccard clustering method with a resolution of 0.4 was used to cluster and classify the cell types. Single-cell RNA-sequencing data were submitted to SRA repository through the access code: PRJNA865464.

2.7. Calculating Relative RNA Content from Each Cell Type

The expression level of each gene was calculated as the mean of the three biological replicates performed for each tissue in units of fragments per kilobase of transcript per million mapped reads (FPKM). Only genes with FPKM > 1 were included in our analysis [23]. After removing testis-specific or testis-biased genes and genes with effective lengths less than 100 bp, in total, 433 X genes and 13,003 autosome genes were present in the sheep genome. The relative RNA content per cell was estimated using the total RNA expression of X- and autosome-related genes in each cell divided by the number of cells. As each mRNA was randomly linked to one Unique Molecular Index (UMI/counts) after reverse transcription, the number of mRNA could be calculated by counting different UMIs. This is a proxy for the relative RNA content per cell, rather than a measurement of the actual number of RNA molecules present.

2.8. Comparing RNA Output by Chromosome and Spermatogenic Stage

In this study, we also detected the RNA content produced by individual X chromosome and autosomal genes. The total RNA counts from every X and autosomal gene in every cell type were calculated using sum of UMIs and by dividing this by nGene (the total number of genes on the X chromosomes/autosomes). We then log transformed these counts with y = Log2(counts + 1) and performed nondirectional Wilcoxon tests, with Holm-corrected p values indicating whether genes from the X chromosome were likely to have equal median counts to genes from the autosomes [12].

2.9. Immunofluorescence

Rabbit anti-MSL1 antibody (1:400 dilution, BS-178556R. Bioss, Beijing, China) was used for primary antibody incubation. Alexa Fluor 555-conjugated goat anti-rabbit IgG H&L (1:400) (BS-0295G-AF555. Bioss, Beijing, China) was used to produce immunofluorescence. Finally, DAPI was used for staining. Images were obtained using a confocal laser microscope (LSM, Zeiss, Oberkochen, Germany).

3. Results

3.1. Mapping, Cell Identification and Clustering Analysis

We obtained sheep testis tissues from newborn (0 days), juvenile (90 days) and adult (720 days) sheep. We sequenced 12,465 newborn, 12,388 juvenile and 8900 adult cells, and 11,000, 9811 and 5785 cells, respectively, were filtered for further analysis based on the filtering criteria. In total, 82.4%, 71.3% and 85.2% reads, respectively, were mapped to the sheep genome. In addition, 16,849, 17,571 and 17,441 genes were detected, and approximately 1314, 1151 and 1124 genes, respectively, were identified, on average, in each cell.
After sequencing, cell partitioning took place through the analysis of single-cell transcriptomes. We selected the marker genes of each testicular cell type to identify the cell type of each cluster [24,25,26,27,28]. Clustered corresponding cell types were obtained (Figure 1A). Newborn-stage testis cells were divided into 10 clusters, including 8 somatic cell clusters (macrophages, T cells, VSMs, endothelial, Sertoli cells, mesenchymal cells, Leydig cells and immature Leydig cells) and 2 germ cell clusters (SSCs and spermatogonia). In juvenile and adult sheep, the testis cells were divided into 15 clusters, including 8 somatic cell clusters and 7 germ cell clusters (SSCs, spermatogonia, early primary S’cytes, late primary S’cytes, spermatid-1, spermatids-2 and sperm) (Figure 1B,C and Figure S1).

3.2. X and Autosomal Total Transcription

To investigate the transcription regulation of X chromosomes and autosomes in sheep testes, we counted the Unique Molecular Indices (UMIs, a proxy for RNA content) from X chromosomes and autosomes in all the cells. We observed that the overall patterns for X and autosomes looked very similar in newborns, juveniles and adults (Figure 2 and Table 1), the ratios of X/autosome were the highest in somatic cells, intermediate in SSCs and spermatogonia cells and the lowest in meiotic (early primary S’cytes and late primary S’cytes) and post-meiotic (spermatids and sperm) cells. The relatively reduced ratio of X-chromosome RNA from spermatocytes to sperm suggested that the X chromosome was downregulated in excess of autosomes in meiotic and post-meiotic cells. It was observed that this X/autosome ratio was lower than that in somatic cells. If a cell type was dosage compensated, the ratio of X/A (X/autosome) would be similar to or higher than that found in somatic cells. If a cell type lacked dosage compensation, the ratio would be lower than that of somatic cells.

3.3. X Chromosome and Autosomal Relative RNA Production

As the total RNA content in cells cannot fully reflect the gene expression trend, we also detected the RNA content produced by individual genes of X chromosomes and autosomals at different stages of sperm development. We counted the detected reads per gene in each cell type for X chromosomes and autosomes. We carried out an analysis at different developmental stages and compared the median reads per gene between the X chromosomes and autosomes as a measure of X-to-autosome dosage compensation. The results showed that the median reads per gene from X genes were lower than those from autosomal genes; from this, we concluded that cell types are unaffected by dose compensation (Figure 3A–C and Table 2, adjusted p values < 0.05).

3.4. MSL Expression in Germ Cells

The core MSL complex consisting of MSL1, 2, 3 and MOF is conserved in mammalian species. There is no related study on sheep MSL complex. In our data, we found that only MSL1 and MSL2 were expressed in early germ cells, including SSCs, spermatogonia and early primary S’cytes (Figure 4A). The immunofluorescence staining results showed that MSL1 was located in the spermatogonia and SSCs of newborn stage; meanwhile, it was expressed in the SSCs, spermatogonia, early primary S’cyte and late primary S’cyte of the juvenile stage and adult stage (Figure 4B).

4. Discussion

In this study, we selected three important periods of sheep testes development. The two most important stages of testicular development are pre-maturity and post-maturity [29]. The newborn and juvenile period belongs to pre-maturity; in the former, seminiferous tubules contain infantile spermatogonia and undifferentiated Sertoli cells, and in the latter, both types of cells increase in number and spermatocytes and sperm appears gradually [30,31]. In adults, after maturity, spermatogonia immediately continue the proliferation and differentiation effect, representing an acceleration in the germ cell developmental process and establishing the functional seminiferous epithelium [32]. Firstly, we established an approach to successfully identify all types of testis cells from sheep at different stages with single-cell RNA seq.
In this work, we directly observed dosage compensation in terms of X:autosome count ratios via single-cell RNA seq. The results showed that the X:A ratio in germ cells was lower than that in somatic cells and the total RNA per cell from the X and autosomes reached its peak in SSCs and spermatogonia. Analyses showed a decreased X expression in early primary S’cytes, indicating that the sperm that undergoes meiotic sex chromosome inactivation (MSCI) brings in inactive X chromosomes. In mammals, the X chromosomes are inactivated during meiosis, possibly as a defense mechanism against insertions into unpaired chromatin [33]. This result is consistent with recent findings regarding Drosophila germ cells [12]. Due to the limitations of single-cell RNA seq, this study did not take into account any individual variation in gene expression at different developmental stages.
Somatic dosage compensation mechanisms in major genetic model organisms clearly do not function in the germline. Rastelli and Kuroda showed that MLE does not localize to X chromosomes and MSL3 is not expressed in the Drosophila germline cells of adult testes [34]. However, Gupta and colleagues assayed X-chromosome transcripts using DNA microarrays and revealed the two-fold upregulation of X-linked genes in the germline of adult Drosophila testes [7]. Recently, some transcriptome studies have shown that dosage compensation occurs in pre-meiotic germ cells, while others have shown that no dosage compensation occurs in the germline of adult Drosophila testes [11,13,14].
By comparing the RNA output, we found evidence that X chromosomes produced lower numbers of reads per gene than autosomes in germ cells. This finding is contrary to Witt‘s [12]. Our results showed that there was no dose compensation in the germ cells of sheep testis at different developmental stages. This is the first study of X-chromosome dosage compensation in sheep germ cells.
The MSL complex in male flies mediates the global acetylation of histone H4 lysine 16 (H4K16ac) on single male X chromosomes, leading to upregulation [35,36]. Remarkably, there is also evidence that clarifies the role of the MSL complex in regulating mammalian X inactivation as well as activation [37]. Studies have shown that MSLs do not function in the germline [34]. Notably, apart from MSL2, other MSL-complex-member orthologs exist in many species, in which dosage compensation mechanisms are absent or fundamentally different [38]. At present, there is no relevant study on MSL in sheep testis. In our study, only MSL1 and MSL2 were found. MSLs were located in premeiotic germ cells, which may be involved in regulating mammalian X-chromosome inactivation and activation.

5. Conclusions

Although many questions remain, our study strongly supports the absence of dose compensation in sheep pre-meiotic germ cells. MSL is not related to dose compensation but may be involved in regulating mammalian X-chromosome inactivation and activation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani12172169/s1, Figure S1: The expression patterns of marker genes projected on the UMAP plot. Red indicates high expression and gray indicates low or no expression.

Author Contributions

G.C. and Y.S. conceived and supervised the study and had substantial input into the analysis and all drafts. J.S. conducted the study and created a first draft of the paper. Y.Z., H.S., C.W., D.W., Y.Y., X.L. (Xiunan Li) participated in the experiment. W.Q., H.L. (Haijun Li) and X.L. (Xihe Li) participated in the data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Inner Mongolia Scientific and Technological Achievement Transformation guide project: Industrialization of sex control propagation technology of dairy goats and sheep (NO. 2020CG0078), Science and Technology Planning Project of Inner Mongolia (NO. 2019GG241, NO. 219GG348, NO. 2020ZD0007, NO. 2020GG0055, NO. 2021GG0062) and National Natural Science Foundation of China (NO. 31860689).

Institutional Review Board Statement

All the sheep experimental procedures and protocols were evaluated and authorized by the Animal Care and Use Committee of the Institute of Animal Husbandry, Inner Mongolia Agricultural university (approval number: IMG-2020-06).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Acknowledgments

We thank Wei An for critical comments and support. We thank Aowoteke Animal Husbandry Technology Co., Ltd., (Chifeng, Inner mongolia, China) for providing the experimental animals. We thank Liheng Bai and He Lv helped to revise the manuscript.

Conflicts of Interest

The authors declare that there is no conflict of interest.

References

  1. Huang, K.; Zeng, Q.; Feng, Y.; Hu, Y.; An, Q.; Li, T.; Qin, L.J.; Liu, J.Y.; Xue, Z.; Fan, G. X-chromosome dosage compensation dynamics in human early embryos. bioRxiv 2020. [Google Scholar] [CrossRef]
  2. Ohno, S. Sex Chromosomes and Sex-Linked Genes; Springer: Berlin/Heidelberg, Germany, 1966. [Google Scholar]
  3. Julien, P.; Brawand, D.; Soumillon, M.; Necsulea, A.; Liechti, A.; Schütz, F.; Daish, T.; Grützner, F.; Kaessmann, H.; Barton, N.H. Mechanisms and Evolutionary Patterns of Mammalian and Avian Dosage Compensation. PLoS Biol. 2012, 10, e1001328. [Google Scholar] [CrossRef] [PubMed]
  4. Duan, J.E.; Flock, K.; Jue, N.; Zhang, M.; Jones, A.; Al Seesi, S.; Mandoiu, I.; Pillai, S.; Hoffman, M.; O’Neill, R. Dosage Compensation and Gene Expression of the X Chromosome in Sheep. G3 Genesgenetics 2018, 9, 305–314. [Google Scholar] [CrossRef] [PubMed]
  5. Di, K.N.; Disteche, C.M. Dosage compensation of the active X chromosome in mammals. Nat. Genet. 2006, 38, 47–53. [Google Scholar]
  6. Duan, J.E.; Jue, N.K.; Jiang, Z.; O’Neill, R.; Tian, X.C. 144 Dosage Compensation and X-Linked Gene Expression in Bovine In Vivo and In Vitro Embryos. Reprod. Fertil. Dev. 2015, 28, 202. [Google Scholar] [CrossRef]
  7. Gupta, V.; Parisi, M.; Sturgill, D.; Nuttall, R.; Doctolero, M.; Dudko, O.K.; Malley, J.D.; Eastman, P.S.; Oliver, B. Global analysis of X-chromosome dosage compensation. J. Biol. 2006, 5, 3. [Google Scholar] [CrossRef]
  8. Deng, X.; Hiatt, J.B.; Nguyen, D.K.; Ercan, S.; Sturgill, D.; Hillier, L.W.; Schlesinger, F.; Davis, C.A.; Reinke, V.J.; Gingeras, T.R.; et al. Evidence for compensatory upregulation of expressed X-linked genes in mammals, Caenorhabditis elegans and Drosophila melanogaster. Nat. Genet. 2011, 43, 1179–1185. [Google Scholar] [CrossRef]
  9. Kharchenko, P.V.; Xi, R.; Park, P.J. Evidence for dosage compensation between the X chromosome and autosomes in mammals. Nat. Genet. 2011, 43, 1167–1169. [Google Scholar] [CrossRef]
  10. Goto, T.; Monk, M. Regulation of X-Chromosome Inactivation in Development in Mice and Humans. Microbiol. Mol. Biol. Rev. 1998, 62, 362. [Google Scholar] [CrossRef]
  11. Shi, Z.; Lim, C.; Tran, V.; Cui, K.; Zhao, K.; Chen, X. Single-cyst transcriptome analysis of Drosophila male germline stem cell lineage. Development 2020, 147, dev.184259. [Google Scholar] [CrossRef]
  12. Witt, E.; Shao, Z.; Hu, C.; Krause, H.M.; Zhao, L. Single-cell RNA-sequencing reveals pre-meiotic X-chromosome dosage compensation in Drosophila testis. PLoS Genet. 2021, 17, e1009728. [Google Scholar] [CrossRef] [PubMed]
  13. Argyridou, E.; Parsch, J. Regulation of the X Chromosome in the Germline and Soma of Drosophila melanogaster Males. Genes 2018, 9, 242. [Google Scholar] [CrossRef] [PubMed]
  14. Meiklejohn, C.D.; Landeen, E.L.; Cook, J.M.; Kingan, S.B.; Presgraves, D.C. Sex Chromosome-Specific Regulation in the Drosophila Male Germline But Little Evidence for Chromosomal Dosage Compensation or Meiotic Inactivation. PLoS Biol. 2011, 9, e1001126. [Google Scholar] [CrossRef] [PubMed]
  15. Chow, J.; Heard, E. X inactivation and the complexities of silencing a sex chromosome. Curr. Opin. Cell Biol. 2009, 21, 359–366. [Google Scholar] [CrossRef]
  16. Dmitriev, R.I.; Pestov, N.B.; Korneenko, T.V.; Gerasimova, A.V.; Zhao, H.; Modyanov, N.N.; Kostina, M.B.; Shakhparonov, M.I. Tissue Specificity of Alternative Splicing of Transcripts Encoding Hampin, a New Mouse Protein Homologous to the Drosophila MSL-1 Protein. Russ. J. Bioorganic Chem. 2005, 31, 325–331. [Google Scholar] [CrossRef]
  17. Taipale, M.; Rea, S.; Richter, K.; Vilar, A.; Lichter, P.; Imhof, A.; Akhtar, A. hMOF Histone Acetyltransferase Is Required for Histone H4 Lysine 16 Acetylation in Mammalian Cells. Mol. Cell. Biol. 2005, 25, 6798–6810. [Google Scholar] [CrossRef]
  18. Peng, Z.; Du, J.; Sun, B.; Dong, X.; Xu, G.; Zhou, J.; Huang, Q.; Liu, Q.; Quan, H.; Ding, J. Structure of human MRG15 chromo domain and its binding to Lys36-methylated histone H3. Nucleic Acids Res. 2006, 34, 6621–6628. [Google Scholar]
  19. Tominaga, K.; Kirtane, B.; Jackson, J.G.; Ikeno, Y.; Ikeda, T.; Hawks, C.; Smith, J.R.; Matzuk, M.M.; Pereira-Smith, O.M. MRG15 Regulates Embryonic Development and Cell Proliferation. Mol. Cell. Biol. 2005, 25, 2924–2937. [Google Scholar] [CrossRef]
  20. Dmitriev, R.I.; Korneenko, T.V.; Bessonov, A.A.; Shakhparonov, M.I.; Modyanov, N.N.; Pestov, N.B. Characterization of hampin/MSL1 as a node in the nuclear interactome. Biochem. Biophys. Res. Commun. 2007, 355, 1051–1057. [Google Scholar] [CrossRef]
  21. Wu, L.; Zee, B.M.; Wang, Y.; Garcia, B.A.; Dou, Y. The RING Finger Protein MSL2 in the MOF Complex Is an E3Ubiquitin Ligase for H2B K34 and Is Involved in Crosstalk with H3 K4 and K79 Methylation. Mol. Cell 2011, 43, 132–144. [Google Scholar] [CrossRef]
  22. Deng, X.; Berletch, J.B.; Ma, W.; Nguyen, D.K.; Hiatt, J.B.; Noble, W.S.; Shendure, J.; Disteche, C.M. Mammalian X upregulation is associated with enhanced transcription initiation, RNA half-life, and MOF-mediated H4K16 acetylation. Dev. Cell 2013, 25, 55–68. [Google Scholar] [CrossRef] [PubMed]
  23. Sangrithi, M.N.; Royo, H.; Mahadevaiah, S.K.; Ojarikre, O.; Bhaw, L.; Sesay, A.; Peters, A.H.; Stadler, M.; Turner, J.M. Non-Canonical and Sexually Dimorphic X Dosage Compensation States in the Mouse and Human Germline. Dev. Cell 2017, 40, 289–301.e283. [Google Scholar] [CrossRef] [PubMed]
  24. Yang, H.; Ma, J.; Wan, Z.; Wang, Q.; Zhang, Y. Characterization of sheep spermatogenesis through single cell RNA sequencing. FASEB J. 2020, 35, e21187. [Google Scholar] [CrossRef] [PubMed]
  25. Guo, J.; Sosa, E.; Chitiashvili, T.; Nie, X.; Rojas, E.J.; Oliver, E.; DonorConnect; Plath, K.; Hotaling, J.M.; Stukenborg, J.B. Single-cell analysis of the developing human testis reveals somatic niche cell specification and fetal germline stem cell establishment. Cell Stem Cell 2021, 28, 764–778. [Google Scholar] [CrossRef] [PubMed]
  26. The adult human testis transcriptional cell atlas. Cell Res. 2018, 28, 1141–1157. [CrossRef] [PubMed]
  27. Daniel, G.C.; Ma, Q.; Manske, G.L.; Niederriter, S.A.; Xianing, Z.; Simone, M.; Lindsay, M.; Caleb, S.; Gurczynski, S.J.; Moore, B.B. A Comprehensive Roadmap of Murine Spermatogenesis Defined by Single-Cell RNA-Seq. Dev. Cell 2018, 46, 651–667. [Google Scholar]
  28. Lukassen, S.; Bosch, E.; Ekici, A.B.; Winterpacht, A. Characterization of germ cell differentiation in the male mouse through single-cell RNA sequencing. Sci. Rep. 2018, 8, 6521. [Google Scholar] [CrossRef] [PubMed]
  29. Turner, J. Meiotic Silencing in Mammals. Annu. Rev. Genet. 2015, 49, 395–412. [Google Scholar] [CrossRef] [PubMed]
  30. Sharpe, R.; Mckinnell, C.; Kivlin, C.; Fisher, J. Proliferation and functional maturation of Sertoli cells, and their relevance to disorders of testis function in adulthood. Reproduction 2003, 125, 769–784. [Google Scholar] [CrossRef]
  31. Meyts, E.R.D.; Almstrup, K.; Nielsen, J.E.; Skakkebæk, N.E. The Testis in Childhood Between Birth and Puberty; Springer: London, UK, 2013; pp. 69–75. [Google Scholar] [CrossRef]
  32. Skandalakis, J.; Gray, S.; Rowe, J.J. Testis and epididymis. J. Reprod. Fertil. 1983, 58, 395–399. [Google Scholar]
  33. Gu, L.; Walters, J.R. Evolution of Sex Chromosome Dosage Compensation in Animals: A Beautiful Theory, Undermined by Facts and Bedeviled by Details. Genome Biol. Evol. 2017, 9, 2461–2476. [Google Scholar] [CrossRef] [PubMed]
  34. Rastelli, L.; Kuroda, M.I. An analysis of maleless and histone H4 acetylation in Drosophila melanogaster spermatogenesis. Mech. Dev. 1998, 71, 107–117. [Google Scholar] [CrossRef]
  35. Bone, J.R.; Lavender, J.; Richman, R.; Palmer, M.J.; Turner, B.M.; Kuroda, M.I. Acetylated histone H4 on the male X chromosome is associated with dosage compensation in Drosophila. Genes Dev. 1994, 8, 96–104. [Google Scholar] [CrossRef] [PubMed]
  36. Hamada, F.N.; Park, P.J.; Gordadze, P.R.; Kuroda, M.I. Global regulation of X chromosomal genes by the MSL complex in Drosophila melanogaster. Genes Dev. 2005, 19, 2289–2294. [Google Scholar] [CrossRef] [PubMed]
  37. Keller, C.I.; Akhtar, A. The MSL complex: Juggling RNA-protein interactions for dosage compensation and beyond. Curr. Opin. Genet. Dev. 2015, 31, 1–11. [Google Scholar] [CrossRef]
  38. Raja, S.J.; Charapitsa, I.; Conrad, T.; Vaquerizas, J.M.; Gebhardt, P.; Holz, H.; Kadlec, J.; Fraterman, S.; Luscombe, N.M.; Akhtar, A. The nonspecific lethal complex is a transcriptional regulator in Drosophila. Mol. Cell 2010, 38, 827–841. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Clustering and cell−type assignment of single cells in Seurat. (A) An illustration of the major cell types in the testis. (B). A UMAP projection of every cell type identified in the data. (C). Expression distribution of marker genes in each cell types.
Figure 1. Clustering and cell−type assignment of single cells in Seurat. (A) An illustration of the major cell types in the testis. (B). A UMAP projection of every cell type identified in the data. (C). Expression distribution of marker genes in each cell types.
Animals 12 02169 g001aAnimals 12 02169 g001b
Figure 2. P0, P90, P720 relative RNA content from each cell type in testis. For every cell, the number of RNA molecules detected from the X chromosome and autosomes approximated by the summed number of Unique Molecular Indices counted per cell. Relative RNA content from each cell type in newborn (P0). Relative RNA content from each cell type in juvenile (P90). Relative RNA content from each cell type in adult. In juvenile and adult (P720).
Figure 2. P0, P90, P720 relative RNA content from each cell type in testis. For every cell, the number of RNA molecules detected from the X chromosome and autosomes approximated by the summed number of Unique Molecular Indices counted per cell. Relative RNA content from each cell type in newborn (P0). Relative RNA content from each cell type in juvenile (P90). Relative RNA content from each cell type in adult. In juvenile and adult (P720).
Animals 12 02169 g002
Figure 3. Dosage compensation equalizes X and autosomal transcription in select cell types. (AC) are the median log-transformed (log2(counts+1) counts for every gene in X and autosomes of P0, P90 and P720.
Figure 3. Dosage compensation equalizes X and autosomal transcription in select cell types. (AC) are the median log-transformed (log2(counts+1) counts for every gene in X and autosomes of P0, P90 and P720.
Animals 12 02169 g003
Figure 4. Expression of MSL1 in sheep testis. (A). Expression patterns for MSL in testes scRNA-seq data. (B). MSL1 immunostaining of sheep testis sections at each age. Bar = 50 μm. DAPI, bule; MSL1, red.
Figure 4. Expression of MSL1 in sheep testis. (A). Expression patterns for MSL in testes scRNA-seq data. (B). MSL1 immunostaining of sheep testis sections at each age. Bar = 50 μm. DAPI, bule; MSL1, red.
Animals 12 02169 g004
Table 1. Mean counts per cell from the X chromosome and autosomes in every cell type.
Table 1. Mean counts per cell from the X chromosome and autosomes in every cell type.
Cell TypeX CountsAutosome CountsX/Autosome Ratio
P0
CD16+ Macrophage128.53534641.21190.0277
Endothelial130.22666034.09770.0216
Immature Leydig cells142.45365328.61180.0267
Leydig cells75.91122764.38170.0275
Mesenchymal70.15642510.04740.0280
T cells75.71432641.50000.0287
VSMs110.24243778.68690.0292
Sertoli cells114.59224240.66680.0270
SSC108.94515488.49450.0198
Spermatogonia34.92311683.30770.0207
P30
CD16+ Macrophage82.33592801.88280.0294
Endothelial132.48905846.19230.0227
Immature Leydig cells101.37843810.38510.0266
Leydig cells57.71912056.52870.0281
Mesenchymal88.03813548.59050.0248
T cells43.54221807.60240.0241
VSMs85.27363195.18570.0267
Sertoli cells124.17194505.43550.0276
SSC157.95638479.79610.0186
Spermatogonia107.16975025.11500.0213
Early primary S’cytes36.08185056.93080.0071
Late primary S’cytes5.5809594.57320.0094
Spermatid-131.20004163.82280.0075
Spermatid-220.10001843.76670.0109
Sperm34.98651622.84750.0216
P90
CD16+ Macrophage62.80872341.62420.0268
Endothelial90.95684067.57600.0224
Immature Leydig cells107.80003975.80000.0271
Leydig cells54.37942121.62860.0256
Mesenchymal cells17.0870791.78260.0216
T cells34.36021462.01580.0235
VSMs58.27912217.84280.0263
Sertoli cells56.29932160.63610.0261
SSC140.05267333.11840.0191
Spermatogonia92.80734542.23740.0204
Early primary S’cytes37.28295132.58540.0073
Late primary S’cytes4.0830531.48570.0077
Spermatid-166.828510,416.33660.0064
Spermatid-2113.713216,945.30000.0067
Sperm14.23882233.58530.0064
Table 2. X and autosomal count ratios detected per gene for every cell type.
Table 2. X and autosomal count ratios detected per gene for every cell type.
Cell TypeX/Autosome Ratio of Raw CountsRAW p ValueAdjusted p Value
P0
CD16+ Macrophage0.587.22 × 10−465.1 × 10−45
Endothelial0.579.69 × 10−445.8 × 10−43
Immature Leydig cells0.552.23 × 10−791.8 × 10-78
Leydig cells0.5400
Mesenchymal cells0.542.27 × 10−361.1 × 10−35
T cells0.530.000122070.00024
VSMs0.435.78 × 10−182.3 × 10−17
Sertoli cells0.5500
SSCs0.391.21 × 10−163.6 × 10−16
Spermatogonia0.410.0002441410.00024
P30
CD16+ Macrophage0.539.69 × 10−447.8 × 10−43
Endothelial0.486.88 × 10−2329.6 × 10−231
Immature Leydig cells0.494.95 × 10−262 × 10−25
Leydig cells0.5500
Mesenchymal cells0.565.92 × 10−191.8 × 10−18
T cells0.532.54 × 10−155.1 × 10−15
VSMs0.454.37 × 10−523.9 × 10−51
Sertoli cells0.581.39 × 10−851.4 × 10−84
SSCs0.371.50 × 10−359 × 10−35
Spermatogonia0.422.62 × 10−1453.4 × 10−144
Early primary S’cytes0.147.70 × 10−283.9 × 10−27
Late primary S’cytes0.196.02 × 10−1087.2 × 10−107
Spermatid-10.154.72 × 10−955.2 × 10−94
Spermatid-20.221.67 × 10−111.7 × 10−11
Sperm0.432.46 × 10−381.7 × 10−37
P720
CD16+ Macrophage0.521.29863 × 10−506.5 × 10−50
Endothelial0.475.74 × 10−1917.5 × 10−190
Immature Leydig cells0.430.06250.062
Leydig cells0.5200
Mesenchymal cells0.512.38419 × 10−070.00000048
T cells0.542.4637 × 10−1052.7 × 10−104
VSMs0.443.1914 × 10−622.9 × 10−61
Sertoli cells0.535.85977 × 10−502.3 × 10−49
SSCs0.383.67866 × 10−141.1 × 10−13
Spermatogonia0.412.00458 × 10−601.6 × 10−59
Early primary S’cytes0.142.77632 × 10−601.9 × 10−59
Late primary S’cytes0.154.5374 × 10−2246.4 × 10−223
Spermatid-10.132.0614 × 10−521.2 × 10−51
Spermatid-20.131.6026 × 10−881.6 × 10−87
Sperm0.132.199 × 10−1262.6 × 10−125
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Su, J.; Zhang, Y.; Su, H.; Wang, C.; Wang, D.; Yang, Y.; Li, X.; Qi, W.; Li, H.; Li, X.; et al. Dosage Compensation of the X Chromosome during Sheep Testis Development Revealed by Single-Cell RNA Sequencing. Animals 2022, 12, 2169. https://doi.org/10.3390/ani12172169

AMA Style

Su J, Zhang Y, Su H, Wang C, Wang D, Yang Y, Li X, Qi W, Li H, Li X, et al. Dosage Compensation of the X Chromosome during Sheep Testis Development Revealed by Single-Cell RNA Sequencing. Animals. 2022; 12(17):2169. https://doi.org/10.3390/ani12172169

Chicago/Turabian Style

Su, Jie, Yue Zhang, Hong Su, Caiyun Wang, Daqing Wang, Yanyan Yang, Xiunan Li, Wangmei Qi, Haijun Li, Xihe Li, and et al. 2022. "Dosage Compensation of the X Chromosome during Sheep Testis Development Revealed by Single-Cell RNA Sequencing" Animals 12, no. 17: 2169. https://doi.org/10.3390/ani12172169

APA Style

Su, J., Zhang, Y., Su, H., Wang, C., Wang, D., Yang, Y., Li, X., Qi, W., Li, H., Li, X., Song, Y., & Cao, G. (2022). Dosage Compensation of the X Chromosome during Sheep Testis Development Revealed by Single-Cell RNA Sequencing. Animals, 12(17), 2169. https://doi.org/10.3390/ani12172169

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop