Non-canonical Drosophila X chromosome dosage compensation and repressive topologically associated domains

Cell-Type Specific Variation in X-Chromosome Dosage Compensation in Drosophila

Male Drosophila melanogaster require dosage compensation to equalize X-linked gene expression with autosomal expression. Leveraging the single-nucleus Fly Cell Atlas (FCA) dataset, which includes 388,918 nuclei across diverse tissues, we investigated cell-type-specific patterns of X-chromosome dosage compensation. Our analysis identified a continuum of cell groups based on their X-to-autosome (X/A) expression ratios ranging from anti-compensated to effectively compensated and overcompensated. Anti-compensation was predominantly observed in male reproductive tissues, while overcompensation was prevalent in neural cells. The expression levels of the dosage compensation machinery's non-coding RNAs, RoX1 and RoX2 , correlated with compensation levels, but were insufficient to fully explain the observed patterns of compensation. These findings reveal the complexity of dosage compensation and suggest that its regulation by the RoX RNAs is nonlinear, implicating potential alternative mechanisms in certain cell types.

Set2 and H3K36 regulate the Drosophila male X chromosome in a context-specific manner, independent from MSL complex spreading

Dosage compensation in Drosophila involves upregulating male X-genes two-fold. This process is carried out by the MSL (male-specific lethal) complex, which binds high-affinity sites and spreads to surrounding genes. Current models of MSL spreading focus on interactions betwen MSL3 (male-specific lethal 3) and Set2-dependent histone marks like trimethylated H3 lysine-36 (H3K36me3). However, Set2 could affect DC via another target, or there could be redundancy between canonical H3.2 and variant H3.3 histones. Furthermore, it is important to parse male-specific effects from those that are X-specific. To discriminate among these possibilities, we employed genomic approaches in H3K36 'residue' and Set2 'writer' mutants. The results confirm a role for Set2 in X-gene regulation, but show that expression trends in males are often mirrored in females. Instead of global, male-specific reduction of X-genes in Set2 or H3K36 mutants, we observe heterogeneous effects. Interestingly, we identified groups of differentially expressed genes (DEGs) whose changes were in opposite directions following loss of H3K36 or Set2, suggesting that H3K36me states have reciprocal functions. In contrast to H4K16R controls, differential expression analysis of combined H3.2K36R/H3.3K36R mutants showed neither consistent reduction in X-gene expression, nor correlation with MSL3 binding. Motif analysis of the DEGs implicated BEAF-32 and other insulator proteins in Set2/H3K36-dependent regulation. Overall, the data are inconsistent with the prevailing model wherein H3K36me3 is essential for spreading the MSL complex to genes along the male X. Rather, we propose that Set2 and H3K36 support DC indirectly, via processes that are utilized by MSL but common to both sexes.

Set2 and H3K36 regulate the Drosophila male X chromosome in a context-specific manner, independent from MSL complex spreading

Dosage compensation in Drosophila involves upregulating male X-genes two-fold. This process is carried out by the MSL (male-specific lethal) complex, which binds high-affinity sites and spreads to surrounding genes. Current models of MSL spreading focus on interactions of MSL3 (male-specific lethal 3) with histone marks; in particular, Set2-dependent H3 lysine-36 trimethylation (H3K36me3). However, Set2 might affect DC via another target, or there could be redundancy between canonical H3.2 and variant H3.3 histones. Further, it is difficult to parse male-specific effects from those that are simply X-specific. To discriminate among these possibilities, we employed genomic approaches in H3K36 (residue) and Set2 (writer) mutants. The results confirm a role for Set2 in X-gene regulation, but show that expression trends in males are often mirrored in females. Instead of global male-specific reduction of X-genes in Set2/H3K36 mutants, the effects were heterogeneous. We identified cohorts of genes whose expression was significantly altered following loss of H3K36 or Set2, but the changes were in opposite directions, suggesting that H3K36me states have reciprocal functions. In contrast to H4K16R controls, analysis of combined H3.2K36R/H3.3K36R mutants neither showed consistent reduction in X-gene expression, nor any correlation with MSL3 binding. Examination of other developmental stages/tissues revealed additional layers of context-dependence. Our studies implicate BEAF-32 and other insulator proteins in Set2/H3K36-dependent regulation. Overall, the data are inconsistent with the prevailing model wherein H3K36me3 directly recruits the MSL complex. We propose that Set2 and H3K36 support DC indirectly, via processes that are utilized by MSL but common to both sexes.

Integrating massive RNA-seq data to elucidate transcriptome dynamics in Drosophila melanogaster

The volume of ribonucleic acid (RNA)-seq data has increased exponentially, providing numerous new insights into various biological processes. However, due to significant practical challenges, such as data heterogeneity, it is still difficult to ensure the quality of these data when integrated. Although some quality control methods have been developed, sample consistency is rarely considered and these methods are susceptible to artificial factors. Here, we developed MassiveQC, an unsupervised machine learning-based approach, to automatically download and filter large-scale high-throughput data. In addition to the read quality used in other tools, MassiveQC also uses the alignment and expression quality as model features. Meanwhile, it is user-friendly since the cutoff is generated from self-reporting and is applicable to multimodal data. To explore its value, we applied MassiveQC to Drosophila RNA-seq data and generated a comprehensive transcriptome atlas across 28 tissues from embryogenesis to adulthood. We systematically characterized fly gene expression dynamics and found that genes with high expression dynamics were likely to be evolutionarily young and expressed at late developmental stages, exhibiting high nonsynonymous substitution rates and low phenotypic severity, and they were involved in simple regulatory programs. We also discovered that human and Drosophila had strong positive correlations in gene expression in orthologous organs, revealing the great potential of the Drosophila system for studying human development and disease.

Increased gene dosage and mRNA expression from chromosomal duplications in Caenorhabditis elegans

Isolation of copy number variations and chromosomal duplications at high frequency in the laboratory suggested that Caenorhabditis elegans tolerates increased gene dosage. Here, we addressed if a general dosage compensation mechanism acts at the level of mRNA expression in C. elegans. We characterized gene dosage and mRNA expression in 3 chromosomal duplications and a fosmid integration strain using DNA-seq and mRNA-seq. Our results show that on average, increased gene dosage leads to increased mRNA expression, pointing to a lack of genome-wide dosage compensation. Different genes within the same chromosomal duplication show variable levels of mRNA increase, suggesting feedback regulation of individual genes. Somatic dosage compensation and germline repression reduce the level of mRNA increase from X chromosomal duplications. Together, our results show a lack of genome-wide dosage compensation mechanism acting at the mRNA level in C. elegans and highlight the role of epigenetic and individual gene regulation contributing to the varied consequences of increased gene dosage.

Sex-specific variation in R-loop formation in Drosophila melanogaster

R-loops are three-stranded nucleotide structures consisting of a DNA:RNA hybrid and a displaced ssDNA non-template strand. Previous work suggests that R-loop formation is primarily determined by the thermodynamics of DNA:RNA binding, which are governed by base composition (e.g., GC skew) and transcription-induced DNA superhelicity. However, R-loops have been described at genomic locations that lack these properties, suggesting that they may serve other context-specific roles. To better understand the genetic determinants of R-loop formation, we have characterized the Drosophila melanogaster R-loop landscape across strains and between sexes using DNA:RNA immunoprecipitation followed by high-throughput sequencing (DRIP-seq). We find that R-loops are associated with sequence motifs that are G-rich or exhibit G/C skew, as well as highly expressed genes, tRNAs, and small nuclear RNAs, consistent with a role for DNA sequence and torsion in R-loop specification. However, we also find motifs associated with R-loops that are A/T-rich and lack G/C skew as well as a subset of R-loops that are enriched in polycomb-repressed chromatin. Differential enrichment analysis reveals a small number of sex-biased R-loops: while non-differentially enriched and male-enriched R-loops form at similar genetic features and chromatin states and contain similar sequence motifs, female-enriched R-loops form at unique genetic features, chromatin states, and sequence motifs and are associated with genes that show ovary-biased expression. Male-enriched R-loops are most abundant on the dosage-compensated X chromosome, where R-loops appear stronger compared to autosomal R-loops. R-loop-containing genes on the X chromosome are dosage-compensated yet show lower MOF binding and reduced H4K16ac compared to R-loop-absent genes, suggesting that H4K16ac or MOF may attenuate R-loop formation. Collectively, these results suggest that R-loop formation in vivo is not fully explained by DNA sequence and topology and raise the possibility that a distinct subset of these hybrid structures plays an important role in the establishment and maintenance of epigenetic differences between sexes.

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.

When Down Is Up: Heterochromatin, Nuclear Organization and X Upregulation

Organisms with highly differentiated sex chromosomes face an imbalance in X-linked gene dosage. Male Drosophila solve this problem by increasing expression from virtually every gene on their single X chromosome, a process known as dosage compensation. This involves a ribonucleoprotein complex that is recruited to active, X-linked genes to remodel chromatin and increase expression. Interestingly, the male X chromosome is also enriched for several proteins associated with heterochromatin. Furthermore, the polytenized male X is selectively disrupted by the loss of factors involved in repression, silencing, heterochromatin formation or chromatin remodeling. Mutations in many of these factors preferentially reduce male survival or enhance the lethality of mutations that prevent normal recognition of the X chromosome. The involvement of primarily repressive factors in a process that elevates expression has long been puzzling. Interestingly, recent work suggests that the siRNA pathway, often associated with heterochromatin formation and repression, also helps the dosage compensation machinery identify the X chromosome. In light of this finding, we revisit the evidence that links nuclear organization and heterochromatin to regulation of the male X chromosome.

Rapid Evolution of Autosomal Binding Sites of the Dosage Compensation Complex in Drosophila melanogaster and Its Association With Transcription Divergence

How pleiotropy influences evolution of protein sequence remains unclear. The male-specific lethal (MSL) complex in Drosophila mediates dosage compensation by 2-fold upregulation of the X chromosome in males. Nevertheless, several MSL proteins also bind autosomes and likely perform functions not related to dosage compensation. Here, we study the evolution of MOF, MSL1, and MSL2 biding sites in Drosophila melanogaster and its close relative Drosophila simulans. We found pervasive expansion of the MSL binding sites in D. melanogaster, particularly on autosomes. The majority of these newly-bound regions are unlikely to function in dosage compensation and associated with an increase in expression divergence between D. melanogaster and D. simulans. While dosage-compensation related sites show clear signatures of adaptive evolution, these signatures are even more marked among autosomal regions. Our study points to an intriguing avenue of investigation of pleiotropy as a mechanism promoting rapid protein sequence evolution.