Mouse lemurs are the smallest and fastest reproducing primates, as well as one of the most abundant, and they are emerging as a model organism for primate biology, behaviour, health and conservation. Although much has been learnt about their ecology and phylogeny in Madagascar and their physiology, little is known about their cellular and molecular biology. Here we used droplet-based and plate-based single-cell RNA sequencing to create Tabula Microcebus, a transcriptomic atlas of 226,000 cells from 27 mouse lemur organs opportunistically obtained from four donors clinically and histologically characterized. Using computational cell clustering, integration and expert cell annotation, we define and biologically organize more than 750 lemur molecular cell types and their full gene expression profiles. This includes cognates of most classical human cell types, including stem and progenitor cells, and differentiating cells along the developmental trajectories of spermatogenesis, haematopoiesis and other adult tissues. We also describe dozens of previously unidentified or sparsely characterized cell types. We globally compare expression profiles to define the molecular relationships of cell types across the body, and explore primate cell and gene expression evolution by comparing lemur transcriptomes to those of human, mouse and macaque. This reveals cell-type-specific patterns of primate specialization and many cell types and genes for which the mouse lemur provides a better human model than mouse 1 . The atlas provides a cellular and molecular foundation for studying this model primate and establishes a general approach for characterizing other emerging model organisms.

ABSTRACT Mouse lemurs are the smallest, fastest reproducing, and among the most abundant primates, and an emerging model organism for primate biology, behavior, health and conservation. Although much has been learned about their physiology and their Madagascar ecology and phylogeny, little is known about their cellular and molecular biology. Here we used droplet- and plate-based single cell RNA-sequencing to profile 226,000 cells from 27 mouse lemur organs and tissues opportunistically procured from four donors clinically and histologically characterized. Using computational cell clustering, integration, and expert cell annotation, we defined and biologically organized over 750 mouse lemur molecular cell types and their full gene expression profiles. These include cognates of most classical human cell types, including stem and progenitor cells, and the developmental programs for spermatogenesis, hematopoiesis, and other adult tissues. We also described dozens of previously unidentified or sparsely characterized cell types and subtypes. We globally compared cell type expression profiles to define the molecular relationships of cell types across the body, and explored primate cell and gene expression evolution by comparing mouse lemur cell transcriptomes to those of human, mouse, and macaque. This revealed cell type specific patterns of primate specialization, as well as many cell types and genes for which lemur provides a better human model than mouse. The atlas provides a cellular and molecular foundation for studying this primate model organism, and establishes a general approach for other emerging model organisms.

Cell type-specific transcriptional programs that drive differentiation of specialized cell types are key players in development and tissue regeneration. One of the most dramatic changes in the transcription program in Drosophila occurs with the transition from proliferating spermatogonia to differentiating spermatocytes, with >3000 genes either newly expressed or expressed from new alternative promoters in spermatocytes. Here we show that opening of these promoters from their closed state in precursor cells requires function of the spermatocyte-specific tMAC complex, localized at the promoters. The spermatocyte-specific promoters lack the previously identified canonical core promoter elements except for the Inr. Instead, these promoters are enriched for the binding site for the TALE-class homeodomain transcription factors Achi/Vis and for a motif originally identified under tMAC ChIP-seq peaks. The tMAC motif resembles part of the previously identified 14-bp β2UE1 element critical for spermatocyte-specific expression. Analysis of downstream sequences relative to transcription start site usage suggested that ACA and CNAAATT motifs at specific positions can help promote efficient transcription initiation. Our results reveal how promoter-proximal sequence elements that recruit and are acted upon by cell type-specific chromatin binding complexes help establish a robust, cell type-specific transcription program for terminal differentiation.

Repressive H3K27me3 and active H3K4me2/3 together form bivalent chromatin domains, molecular hallmarks of developmental potential. In the male germline, these domains are thought to persist into sperm to establish totipotency in the next generation. However, it remains unknown how H3K27me3 is established on specific targets in the male germline. Here, we demonstrate that a germline-specific Polycomb protein, SCML2, binds to H3K4me2/3-rich hypomethylated promoters in undifferentiated spermatogonia to facilitate H3K27me3. Thus, SCML2 establishes bivalent domains in the male germline of mice. SCML2 regulates two major classes of bivalent domains: Class I domains are established on developmental regulator genes that are silent throughout spermatogenesis, while class II domains are established on somatic genes silenced during late spermatogenesis. We propose that SCML2-dependent H3K27me3 in the male germline prepares the expression of developmental regulator and somatic genes in embryonic development.

Precise epigenetic regulation of the sex chromosomes is vital for the male germline. Here, we analyze meiosis in eight mouse models deficient for various DNA damage response (DDR) factors, including Fanconi anemia (FA) proteins. We reveal a network of FA and DDR proteins in which FA core factors FANCA, FANCB, and FANCC are essential for FANCD2 foci formation, whereas BRCA1 (FANCS), MDC1, and RNF8 are required for BRCA2 (FANCD1) and SLX4 (FANCP) accumulation on the sex chromosomes during meiosis. In addition, FA proteins modulate distinct histone marks on the sex chromosomes: FA core proteins and FANCD2 regulate H3K9 methylation, while FANCD2 and RNF8 function together to regulate H3K4 methylation independently of FA core proteins. Our data suggest that RNF8 integrates the FA-BRCA pathway. Taken together, our study reveals distinct functions for FA proteins and illuminates the male sex chromosomes as a model to dissect the function of the FA-BRCA pathway.

Gametogenesis is dependent on the expression of germline-specific genes. However, it remains unknown how the germline epigenome is distinctly established from that of somatic lineages. Here we show that genes commonly expressed in somatic lineages and spermatogenesis-progenitor cells undergo repression in a genome-wide manner in late stages of the male germline and identify underlying mechanisms. SCML2, a germline-specific subunit of a Polycomb repressive complex 1 (PRC1), establishes the unique epigenome of the male germline through two distinct antithetical mechanisms. SCML2 works with PRC1 and promotes RNF2-dependent ubiquitination of H2A, thereby marking somatic/progenitor genes on autosomes for repression. Paradoxically, SCML2 also prevents RNF2-dependent ubiquitination of H2A on sex chromosomes during meiosis, thereby enabling unique epigenetic programming of sex chromosomes for male reproduction. Our results reveal divergent mechanisms involving a shared regulator by which the male germline epigenome is distinguished from that of the soma and progenitor cells.

Fanconi anemia (FA) is a recessive X-linked and autosomal genetic disease associated with bone marrow failure and increased cancer, as well as severe germline defects such as hypogonadism and germ cell depletion. Although deficiencies in FA factors are commonly associated with germ cell defects, it remains unknown whether the FA pathway is involved in unique epigenetic events in germ cells. In this study, we generated Fancb mutant mice, the first mouse model of X-linked FA, and identified a novel function of the FA pathway in epigenetic regulation during mammalian gametogenesis. Fancb mutant mice were infertile and exhibited primordial germ cell (PGC) defects during embryogenesis. Further, Fancb mutation resulted in the reduction of undifferentiated spermatogonia in spermatogenesis, suggesting that FANCB regulates the maintenance of undifferentiated spermatogonia. Additionally, based on functional studies, we dissected the pathway in which FANCB functions during meiosis. The localization of FANCB on sex chromosomes is dependent on MDC1, a binding partner of H2AX phosphorylated at serine 139 (γH2AX), which initiates chromosome-wide silencing. Also, FANCB is required for FANCD2 localization during meiosis, suggesting that the role of FANCB in the activation of the FA pathway is common to both meiosis and somatic DNA damage responses. H3K9me2, a silent epigenetic mark, was decreased on sex chromosomes, whereas H3K9me3 was increased on sex chromosomes in Fancb mutant spermatocytes. Taken together, these results indicate that FANCB functions at critical stages of germ cell development and reveal a novel function of the FA pathway in the regulation of H3K9 methylation in the germline.

BACKGROUND:The male germline transcriptome changes dramatically during the mitosis-to-meiosis transition to activate late spermatogenesis genes and to transiently suppress genes commonly expressed in somatic lineages and spermatogenesis progenitor cells, termed somatic/progenitor genes.RESULTS:These changes reflect epigenetic regulation. Induction of late spermatogenesis genes during spermatogenesis is facilitated by poised chromatin established in the stem cell phases of spermatogonia, whereas silencing of somatic/progenitor genes during meiosis and postmeiosis is associated with formation of bivalent domains which also allows the recovery of the somatic/progenitor program after fertilization. Importantly, during spermatogenesis mechanisms of epigenetic regulation on sex chromosomes are different from autosomes: X-linked somatic/progenitor genes are suppressed by meiotic sex chromosome inactivation without deposition of H3K27me3.CONCLUSIONS:Our results suggest that bivalent H3K27me3 and H3K4me2/3 domains are not limited to developmental promoters (which maintain bivalent domains that are silent throughout the reproductive cycle), but also underlie reversible silencing of somatic/progenitor genes during the mitosis-to-meiosis transition in late spermatogenesis.

Epigenetic mechanisms precisely regulate sex chromosome inactivation as well as genes that escape the silencing process. In male germ cells, DNA damage response factor RNF8 establishes active epigenetic modifications on the silent sex chromosomes during meiosis, and activates escape genes during a state of sex chromosome-wide silencing in postmeiotic spermatids. During the course of evolution, the gene content of escape genes in postmeiotic spermatids recently diverged on the sex chromosomes. This evolutionary feature mirrors the epigenetic processes of sex chromosomes in germ cells. In this article, we describe how epigenetic processes have helped to shape the evolution of sex chromosome-linked genes. Furthermore, we compare features of escape genes on sex chromosomes in male germ cells to escape genes located on the single X chromosome silenced during X-inactivation in females, clarifying the distinct evolutionary implications between male and female escape genes.

Sex chromosomes are uniquely subject to chromosome-wide silencing during male meiosis, and silencing persists into post-meiotic spermatids. Against this background, a select set of sex chromosome-linked genes escapes silencing and is activated in post-meiotic spermatids. Here, we identify a novel mechanism that regulates escape gene activation in an environment of chromosome-wide silencing in murine germ cells. We show that RNF8-dependent ubiquitination of histone H2A during meiosis establishes active epigenetic modifications, including dimethylation of H3K4 on the sex chromosomes. RNF8-dependent active epigenetic memory, defined by dimethylation of H3K4, persists throughout meiotic division. Various active epigenetic modifications are subsequently established on the sex chromosomes in post-meiotic spermatids. These RNF8-dependent modifications include trimethylation of H3K4, histone lysine crotonylation (Kcr), and incorporation of the histone variant H2AFZ. RNF8-dependent epigenetic programming regulates escape gene activation from inactive sex chromosomes in post-meiotic spermatids. Kcr accumulates at transcriptional start sites of sex-linked genes activated in an RNF8-dependent manner, and a chromatin conformational change is associated with RNF8-dependent epigenetic programming. Furthermore, we demonstrate that this RNF8-dependent pathway is distinct from that which recognizes DNA double-strand breaks. Our results establish a novel connection between a DNA damage response factor (RNF8) and epigenetic programming, specifically in establishing active epigenetic modifications and gene activation.

Sex chromosome inactivation is essential epigenetic programming in male germ cells. However, it remains largely unclear how epigenetic silencing of sex chromosomes impacts the evolution of the mammalian genome. Here we demonstrate that male sex chromosome inactivation is highly conserved between humans and mice and has an impact on the genetic evolution of human sex chromosomes. We show that, in humans, sex chromosome inactivation established during meiosis is maintained into spermatids with the silent compartment postmeiotic sex chromatin (PMSC). Human PMSC is illuminated with epigenetic modifications such as trimethylated lysine 9 of histone H3 and heterochromatin proteins CBX1 and CBX3, which implicate a conserved mechanism underlying the maintenance of sex chromosome inactivation in mammals. Furthermore, our analyses suggest that male sex chromosome inactivation has impacted multiple aspects of the evolutionary history of mammalian sex chromosomes: amplification of copy number, retrotranspositions, acquisition of de novo genes, and acquisition of different expression profiles. Most strikingly, profiles of escape genes from postmeiotic silencing diverge significantly between humans and mice. Escape genes exhibit higher rates of amino acid changes compared with non-escape genes, suggesting that they are beneficial for reproductive fitness and may allow mammals to cope with conserved postmeiotic silencing during the evolutionary past. Taken together, we propose that the epigenetic silencing mechanism impacts the genetic evolution of sex chromosomes and contributed to speciation and reproductive diversity in mammals.

Sex chromosome inactivation in male germ cells is a paradigm of epigenetic programming during sexual reproduction. Recent progress has revealed the underlying mechanisms of sex chromosome inactivation in male meiosis. The trigger of chromosome-wide silencing is activation of the DNA damage response (DDR) pathway, which is centered on the mediator of DNA damage checkpoint 1 (MDC1), a binding partner of phosphorylated histone H2AX (γH2AX). This DDR pathway shares features with the somatic DDR pathway recognizing DNA replication stress in the S phase. Additionally, it is likely to be distinct from the DDR pathway that recognizes meiosis-specific double-strand breaks. This review article extensively discusses the underlying mechanism of sex chromosome inactivation.

PURPOSE:We identified the endogenous retroviruses associated with TTYs (testis specific transcripts linked to the Y) in the AZFb region. We evaluated the relationship between endogenous retroviruses, and TTY expression patterns and function in spermatogenesis.MATERIALS AND METHODS:We identified family members of TTYs in the AZFb region using computational screening. After investigating the relationship between the endogenous retrovirus genome and TTY expression patterns we screened genomic polymerase chain reaction products from TTY13 amplified from 790 Japanese men, including 275 with azoospermia, 285 with oligozoospermia and 230 who were fertile.RESULTS:Computational screening revealed that 3 members of the TTY family, TTY9, 10 and 13, were regulated by endogenous retroviruses in the AZFb region. Homologous recombination between long terminal repeat of the TTY13 associated human endogenous retrovirus-K14C resulted in TTY13 deletion events. These deletions were more common in patients with azoospermia and oligozoospermia than in fertile males. Specifically 15.63% of the azoospermia group, 10.88% of the oligozoospermia group and 0% of fertile controls had only the deletion variant, indicating an association between the homologous recombination rate and the severity of spermatogenesis failure that was statistically significant (p <0.05).CONCLUSIONS:Because of the finding of what are to our knowledge novel microdeletions due to endogenous retrovirus in the AZFb region, our study raises the possibility that specific variations in genomic structure may contribute to some forms of human idiopathic male infertility.

The male-specific region of Y chromosome (MSY) has accumulated a higher density of human endogenous retroviruses (HERVs) and related sequences when compared with other regions of the human genome. Here, we focused on one HERV family, HERV-K14C that seemed to integrate preferentially into the Y chromosome in humans. To identify every copies of HERV-K14C in the human genome, we applied computational screening to map precisely the locus of individual HERV-K14C copies. Interestingly, 29 of all 146 copies were located in Y chromosome, and these 29 copies were mostly dispersed in the palindromic region. Three distinct HERV-K14C-related transcripts were found and were exclusively expressed in human testis tissue. Based on our phylogenetic analysis of the solitary LTRs derived from HERV-K14C on the Y chromosome we suggested that these sequences were generated as pairs of identical sequences. Specifically, analysis of HERV-K14C-related sequences in the palindromic region demonstrated that the Y chromosomal amplicons existed in our common ancestors and the duplicated pairs arose after divergence of great apes approximately 8–10 million years ago. Taken together, our observation suggested that HERV-K14C-related sequences contributed to genomic diversification of Y chromosome during speciation of great ape lineage.

Recently, work has shown that azoospermia factor (AZF) microdeletions result from homologous recombination between almost identical blocks in this gene region. These microdeletions in the Y chromosome are a common molecular genetic cause of spermatogenetic failure leading to male infertility. After completion of the sequencing of the Y chromosome, the classical definition of AZFa, AZFb, and AZFc was modified to five regions, namely AZFa, P5/proximal-P1, P5/distal-P1, P4/distal-P1, and AZFc, as a result of the determination of Y chromosomal structure. Moreover, partial AZFc deletions have also been reported, resulting from recombination in their sub-ampliconic identical pair sequences. These deletions are also implicated in a possible association with Y chromosome haplogroups. In this review, we address Y chromosomal complexity and the modified categories of the AZF deletions. Recognition of the association of Y deletions with male infertility has implications for the diagnosis, treatment, and genetic counseling of infertile men, in particular candidates for intracytoplasmic sperm injection.