Single-cyst transcriptome analysis of Drosophila male germline stem cell lineage

Chromosome-specific maturation of the epigenome in the Drosophila male germline

Spermatogenesis in the Drosophila male germline proceeds through a unique transcriptional program controlled both by germline-specific transcription factors and by testis-specific versions of core transcriptional machinery. This program includes the activation of genes on the heterochromatic Y chromosome, and reduced transcription from the X chromosome, but how expression from these sex chromosomes is regulated has not been defined. To resolve this, we profiled active chromatin features in the testes from wildtype and meiotic arrest mutants and integrate this with single-cell gene expression data from the Fly Cell Atlas. These data assign the timing of promoter activation for genes with germline-enriched expression throughout spermatogenesis, and general alterations of promoter regulation in germline cells. By profiling both active RNA polymerase II and histone modifications in isolated spermatocytes, we detail widespread patterns associated with regulation of the sex chromosomes. Our results demonstrate that the X chromosome is not enriched for silencing histone modifications, implying that sex chromosome inactivation does not occur in the Drosophila male germline. Instead, a lack of dosage compensation in spermatocytes accounts for the reduced expression from this chromosome. Finally, profiling uncovers dramatic ubiquitinylation of histone H2A and lysine-16 acetylation of histone H4 across the Y chromosome in spermatocytes that may contribute to the activation of this heterochromatic chromosome.

Chromosome-specific maturation of the epigenome in the Drosophila male germline

Spermatogenesis in the Drosophila male germline proceeds through a unique transcriptional program controlled both by germline-specific transcription factors and by testis-specific versions of core transcriptional machinery. This program includes the activation of genes on the heterochromatic Y chromosome, and reduced transcription from the X chromosome, but how expression from these sex chromosomes is regulated has not been defined. To resolve this, we profiled active chromatin features in the testes from wildtype and meiotic arrest mutants and integrate this with single-cell gene expression data from the Fly Cell Atlas. These data assign the timing of promoter activation for genes with germline-enriched expression throughout spermatogenesis, and general alterations of promoter regulation in germline cells. By profiling both active RNA polymerase II and histone modifications in isolated spermatocytes, we detail widespread patterns associated with regulation of the sex chromosomes. Our results demonstrate that the X chromosome is not enriched for silencing histone modifications, implying that sex chromosome inactivation does not occur in the Drosophila male germline. Instead, a lack of dosage compensation in spermatocytes accounts for the reduced expression from this chromosome. Finally, profiling uncovers dramatic ubiquitinylation of histone H2A and lysine-16 acetylation of histone H4 across the Y chromosome in spermatocytes that may contribute to the activation of this heterochromatic chromosome.

Emergent dynamics of adult stem cell lineages from single nucleus and single cell RNA-Seq of Drosophila testes

Proper differentiation of sperm from germline stem cells, essential for production of the next generation, requires dramatic changes in gene expression that drive remodeling of almost all cellular components, from chromatin to organelles to cell shape itself. Here, we provide a single nucleus and single cell RNA-seq resource covering all of spermatogenesis in Drosophila starting from in-depth analysis of adult testis single nucleus RNA-seq (snRNA-seq) data from the Fly Cell Atlas (FCA) study. With over 44,000 nuclei and 6000 cells analyzed, the data provide identification of rare cell types, mapping of intermediate steps in differentiation, and the potential to identify new factors impacting fertility or controlling differentiation of germline and supporting somatic cells. We justify assignment of key germline and somatic cell types using combinations of known markers, in situ hybridization, and analysis of extant protein traps. Comparison of single cell and single nucleus datasets proved particularly revealing of dynamic developmental transitions in germline differentiation. To complement the web-based portals for data analysis hosted by the FCA, we provide datasets compatible with commonly used software such as Seurat and Monocle. The foundation provided here will enable communities studying spermatogenesis to interrogate the datasets to identify candidate genes to test for function in vivo.

Enrichment of Undifferentiated Germline and Somatic Cells from Drosophila Testes

The Drosophila male germline provides a strong model system to understand numerous developmental and cell-biological processes, owing to a well-defined anatomy and cell type markers in combination with various genetic tools available for the Drosophila system. A major weakness of this system has been the difficulty of approaches for obtaining material for biochemical assays, proteomics, and genomic or transcriptomic profiling due to small-size and complex tissues. However, the recent development of techniques has started allowing us the usage of a low amount of material for these analyses and now we can strategize many new experiments. The method for enrichment or isolation of rare populations of cells is still challenging and should meaningfully influence the reliability of the results. Here, we provide our semi-optimized protocol of enrichment of undifferentiated germ cells and somatic cells from non-tumorous Drosophila testis, which we have successfully improved after multiple trials.

Age-related ceRNA networks in adult Drosophila ageing

As Drosophila is an extensively used genetic model system, understanding of its regulatory networks has great significance in revealing the genetic mechanisms of ageing and human diseases. Competing endogenous RNA (ceRNA)-mediated regulation is an important mechanism by which circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs) regulate ageing and age-related diseases. However, extensive analyses of the multiomics (circRNA/miRNA/mRNA and lncRNA/miRNA/mRNA) characteristics of adult Drosophila during ageing have not been reported. Here, differentially expressed circRNAs and microRNAs (miRNAs) between 7 and 42-day-old flies were screened and identified. Then, the differentially expressed mRNAs, circRNAs, miRNAs, and lncRNAs between the 7- and 42-day old flies were analysed to identify age-related circRNA/miRNA/mRNA and lncRNA/miRNA/mRNA networks in ageing Drosophila. Several key ceRNA networks were identified, such as the dme_circ_0009500/dme_miR-289-5p/CG31064, dme_circ_0009500/dme_miR-289-5p/frizzled, dme_circ_0009500/dme_miR-985-3p/Abl, and XLOC_027736/dme_miR-985-3p/Abl XLOC_189909/dme_miR-985-3p/Abl networks. Furthermore, real-time quantitative PCR (qPCR) was used to verify the expression level of those genes. Those results suggest that the discovery of these ceRNA networks in ageing adult Drosophila provide new information for research on human ageing and age-related diseases.

Premeiotic pairing of homologous chromosomes during Drosophila male meiosis

In the early stages of meiosis, maternal and paternal chromosomes pair with their homologous partner and recombine to ensure exchange of genetic information and proper segregation. These events can vary drastically between species and between males and females of the same species. In Drosophila, in contrast to females, males do not form synaptonemal complexes (SCs), do not recombine, and have no crossing over; yet, males are able to segregate their chromosomes properly. Here, we investigated the early steps of homolog pairing in Drosophila males. We found that homolog centromeres are not paired in germline stem cells (GSCs) and become paired in the mitotic region before meiotic entry, similarly to females. Surprisingly, male germline cells express SC proteins, which localize to centromeres and promote pairing. We further found that the SUN/KASH (LINC) complex and microtubules are required for homolog pairing as in females. Chromosome movements in males, however, are much slower than in females and we demonstrate that this slow dynamic is compensated in males by having longer cell cycles. In agreement, slowing down cell cycles was sufficient to rescue pairing-defective mutants in female meiosis. Our results demonstrate that although meiosis differs significantly between males and females, sex-specific cell cycle kinetics integrate similar molecular mechanisms to achieve proper centromere pairing.

Dosage Compensation in Drosophila: Its Canonical and Non-Canonical Mechanisms

Dosage compensation equalizes gene expression in a single male X chromosome with that in the pairs of autosomes and female X chromosomes. In the fruit fly Drosophila, canonical dosage compensation is implemented by the male-specific lethal (MSL) complex functioning in all male somatic cells. This complex contains acetyl transferase males absent on the first (MOF), which performs H4K16 hyperacetylation specifically in the male X chromosome, thus facilitating transcription of the X-linked genes. However, accumulating evidence points to an existence of additional, non-canonical dosage compensation mechanisms operating in somatic and germline cells. In this review, we discuss current advances in the understanding of both canonical and non-canonical mechanisms of dosage compensation in Drosophila.

Diverse cell-specific patterns of alternative polyadenylation in Drosophila

Most genes in higher eukaryotes express isoforms with distinct 3' untranslated regions (3' UTRs), generated by alternative polyadenylation (APA). Since 3' UTRs are predominant locations of post-transcriptional regulation, APA can render such programs conditional, and can also alter protein sequences via alternative last exon (ALE) isoforms. We previously used 3'-sequencing from diverse Drosophila samples to define multiple tissue-specific APA landscapes. Here, we exploit comprehensive single nucleus RNA-sequencing data (Fly Cell Atlas) to elucidate cell-type expression of 3' UTRs across >250 adult Drosophila cell types. We reveal the cellular bases of multiple tissue-specific APA/ALE programs, such as 3' UTR lengthening in differentiated neurons and 3' UTR shortening in spermatocytes and spermatids. We trace dynamic 3' UTR patterns across cell lineages, including in the male germline, and discover new APA patterns in the intestinal stem cell lineage. Finally, we correlate expression of RNA binding proteins (RBPs), miRNAs and global levels of cleavage and polyadenylation (CPA) factors in several cell types that exhibit characteristic APA landscapes, yielding candidate regulators of transcriptome complexity. These analyses provide a comprehensive foundation for future investigations of mechanisms and biological impacts of alternative 3' isoforms across the major cell types of this widely-studied model organism.

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

Simple Summary

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

Abstract

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

The continuum of Drosophila embryonic development at single-cell resolution

Drosophila melanogaster is a powerful, long-standing model for metazoan development and gene regulation. We profiled chromatin accessibility in almost 1 million and gene expression in half a million nuclei from overlapping windows spanning the entirety of embryogenesis. Leveraging developmental asynchronicity within embryo collections, we applied deep neural networks to infer the age of each nucleus, resulting in continuous, multimodal views of molecular and cellular transitions in absolute time. We identify cell lineages; infer their developmental relationships; and link dynamic changes in enhancer usage, transcription factor (TF) expression, and the accessibility of TFs' cognate motifs. With these data, the dynamics of enhancer usage and gene expression can be explored within and across lineages at the scale of minutes, including for precise transitions like zygotic genome activation.

Small RNAs and their protein partners in animal meiosis

Meiosis is characterized by highly regulated transitions in gene expression that require diverse mechanisms of gene regulation. For example, in male mammals, transcription undergoes a global shut-down in early prophase I of meiosis, followed by increasing transcriptional activity into pachynema. Later, as spermiogenesis proceeds, the histones bound to DNA are replaced with transition proteins, which are themselves replaced with protamines, resulting in a highly condensed nucleus with repressed transcriptional activity. In addition, two specialized gene silencing events take place during prophase I: meiotic silencing of unsynapsed chromatin (MSUC), and the sex chromatin specific mechanism, meiotic sex chromosome inactivation (MSCI). Notably, conserved roles for the RNA binding protein (RBP) machinery that functions with small non-coding RNAs have been described as participating in these meiosis-specific mechanisms, suggesting that RNA-mediated gene regulation is critical for fertility in many species. Here, we review roles of small RNAs and their associated RBPs in meiosis-related processes such as centromere function, silencing of unpaired chromatin and meiotic recombination. We will discuss the emerging evidence of non-canonical functions of these components in meiosis.

Differential condensation of sister chromatids acts with Cdc6 to ensure asynchronous S-phase entry in Drosophila male germline stem cell lineage

During Drosophila melanogaster male germline stem cell (GSC) asymmetric division, preexisting old versus newly synthesized histones H3 and H4 are asymmetrically inherited. However, the biological outcomes of this phenomenon have remained unclear. Here, we tracked old and new histones throughout the GSC cell cycle through the use of high spatial and temporal resolution microscopy. We found unique features that differ between old and new histone-enriched sister chromatids, including differences in nucleosome density, chromosomal condensation, and H3 Ser10 phosphorylation. These distinct chromosomal features lead to their differential association with Cdc6, a pre-replication complex component, and subsequent asynchronous DNA replication initiation in the resulting daughter cells. Disruption of asymmetric histone inheritance abolishes differential Cdc6 association and asynchronous S-phase entry, demonstrating that histone asymmetry acts upstream of these critical cell-cycle progression events. Furthermore, disruption of these GSC-specific chromatin features leads to GSC defects, indicating a connection between histone inheritance, cell-cycle progression, and cell fate determination.

Mitotic drive in asymmetric epigenetic inheritance

Asymmetric cell division (ACD) produces two daughter cells with distinct cell fates. This division mode is widely used during development and by adult stem cells during tissue homeostasis and regeneration, which can be regulated by both extrinsic cues such as signaling molecules and intrinsic factors such as epigenetic information. While the DNA replication process ensures that the sequences of sister chromatids are identical, how epigenetic information is re-distributed during ACD has remained largely unclear in multicellular organisms. Studies of Drosophila male germline stem cells (GSCs) have revealed that sister chromatids incorporate pre-existing and newly synthesized histones differentially and segregate asymmetrically during ACD. To understand the underlying molecular mechanisms of this phenomenon, two key questions must be answered: first, how and when asymmetric histone information is established; and second, how epigenetically distinct sister chromatids are distinguished and segregated. Here, we discuss recent advances which help our understanding of this interesting and important cell division mode.

Articulating the "stem cell niche" paradigm through the lens of non-model aquatic invertebrates

Stem cells (SCs) in vertebrates typically reside in "stem cell niches" (SCNs), morphologically restricted tissue microenvironments that are important for SC survival and proliferation. SCNs are broadly defined by properties including physical location, but in contrast to vertebrates and other "model" organisms, aquatic invertebrate SCs do not have clearly documented niche outlines or properties. Life strategies such as regeneration or asexual reproduction may have conditioned the niche architectural variability in aquatic or marine animal groups. By both establishing the invertebrates SCNs as independent types, yet allowing inclusiveness among them, the comparative analysis will allow the future functional characterization of SCNs.

OUP accepted manuscript

Eukaryotic chromosomes are spatially segregated into topologically associating domains (TADs). Some TADs are attached to the nuclear lamina (NL) through lamina-associated domains (LADs). Here, we identified LADs and TADs at two stages of Drosophila spermatogenesis – in bam^Δ86 mutant testes which is the commonly used model of spermatogonia (SpG) and in larval testes mainly filled with spermatocytes (SpCs). We found that initiation of SpC-specific transcription correlates with promoters’ detachment from the NL and with local spatial insulation of adjacent regions. However, this insulation does not result in the partitioning of inactive TADs into sub-TADs. We also revealed an increased contact frequency between SpC-specific genes in SpCs implying their de novo gathering into transcription factories. In addition, we uncovered the specific X chromosome organization in the male germline. In SpG and SpCs, a single X chromosome is stronger associated with the NL than autosomes. Nevertheless, active chromatin regions in the X chromosome interact with each other more frequently than in autosomes. Moreover, despite the absence of dosage compensation complex in the male germline, randomly inserted SpG-specific reporter is expressed higher in the X chromosome than in autosomes, thus evidencing that non-canonical dosage compensation operates in SpG.

The rates of stem cell division determine the cell cycle lengths of its lineage

Adult stem cells and their transit-amplifying progeny alter their proliferation rates to maintain tissue homeostasis. To test how the division rates of stem cells and transit-amplifying progeny affect tissue growth and differentiation, we developed a computation strategy that estimates the average cell-cycle lengths (lifespans) of germline stem cells and their progeny from fixed-tissue demography in the Drosophila testis. Analysis of the wild-type data using this method indicated that during the germline transit-amplification, the cellular lifespans extend by nearly 1.3-fold after the first division and shrink by about 2-folds after the second division. Cell-autonomous perturbations of the stem cell lifespan accordingly altered the lifespans of successive transit-amplifying stages. Remarkably, almost 2-fold alterations in the lifespans of stem cells and their immediate daughters did not affect the subsequent differentiation. The results indicate that the early germline division rates can adjust the following division rates and the onset of differentiation.

•

Prediction of cellular lifespan from the demography of transit-amplifying cells

•

Lifespans of spermatogonial cells change anomalously during transit-amplification

•

Anomalous lifespan extension during transit-amplification precedes the onset of Bam

•

Lifespan changes of early TA stages readjust that of the subsequent stages

Trajectory mapping of the early Drosophila germline reveals controls of zygotic activation and sex differentiation

Germ cells in Drosophila melanogaster are specified maternally shortly after fertilization and are transcriptionally quiescent until their zygotic genome is activated to sustain further development. To understand the molecular basis of this process, we analyzed the progressing transcriptomes of early male and female germ cells at the single-cell level between germline specification and coalescence with somatic gonadal cells. Our data comprehensively cover zygotic activation in the germline genome, and analyses on genes that exhibit germline-restricted expression reveal that polymerase pausing and differential RNA stability are important mechanisms that establish gene expression differences between the germline and soma. In addition, we observe an immediate bifurcation between the male and female germ cells as zygotic transcription begins. The main difference between the two sexes is an elevation in X Chromosome expression in females relative to males, signifying incomplete dosage compensation, with a few select genes exhibiting even higher expression increases. These indicate that the male program is the default mode in the germline that is driven to female development with a second X Chromosome.

Absence of X-chromosome dosage compensation in the primordial germ cells of Drosophila embryos

Dosage compensation is a mechanism that equalizes sex chromosome gene expression between the sexes. In Drosophila, individuals with two X chromosomes (XX) become female, whereas males have one X chromosome (XY). In males, dosage compensation of the X chromosome in the soma is achieved by five proteins and two non-coding RNAs, which assemble into the male-specific lethal (MSL) complex to upregulate X-linked genes twofold. By contrast, it remains unclear whether dosage compensation occurs in the germline. To address this issue, we performed transcriptome analysis of male and female primordial germ cells (PGCs). We found that the expression levels of X-linked genes were approximately twofold higher in female PGCs than in male PGCs. Acetylation of lysine residue 16 on histone H4 (H4K16ac), which is catalyzed by the MSL complex, was undetectable in these cells. In male PGCs, hyperactivation of X-linked genes and H4K16ac were induced by overexpression of the essential components of the MSL complex, which were expressed at very low levels in PGCs. Together, these findings indicate that failure of MSL complex formation results in the absence of X-chromosome dosage compensation in male PGCs.

Differential Histone Distribution Patterns in Induced Asymmetrically Dividing Mouse Embryonic Stem Cells

Wnt3a-coated beads can induce asymmetric divisions of mouse embryonic stem cells (mESCs), resulting in one self-renewed mESC and one differentiating epiblast stem cell. This provides an opportunity for studying histone inheritance pattern at a single-cell resolution in cell culture. Here, we report that mESCs with Wnt3a-bead induction display nonoverlapping preexisting (old) versus newly synthesized (new) histone H3 patterns, but mESCs without Wnt3a beads have largely overlapping patterns. Furthermore, H4K20me2/3, an old histone-enriched modification, displays a higher instance of asymmetric distribution on chromatin fibers from Wnt3a-induced mESCs than those from non-induced mESCs. These locally distinct distributions between old and new histones have both cellular specificity in Wnt3a-induced mESCs and molecular specificity for histones H3 and H4. Given that post-translational modifications at H3 and H4 carry the major histone modifications, our findings provide a mammalian cell culture system to study histone inheritance for maintaining stem cell fate and for resetting it during differentiation.