BEN-solo factors partition active chromatin to ensure proper gene activation in Drosophila

A Drosophila model of mucopolysaccharidosis IIIB

Mucopolysaccharidosis type IIIB is a rare lysosomal storage disorder caused by defects in alpha-N-acetylglucosaminidase (NAGLU) and characterized by severe effects in the central nervous system. Mutations in NAGLU cause accumulation of partially degraded heparan sulfate in lysosomes. The consequences of these mutations on whole-genome gene expression and their causal relationships to neural degeneration remain unknown. Here, we used the functional Drosophila melanogaster ortholog of NAGLU, Naglu, to develop a fly model for MPS IIIB induced by gene deletion (NagluKO), missense (NagluY160C), and nonsense (NagluW422X) mutations. We used the Drosophila activity monitoring system to analyze activity and sleep and found sex- and age-dependent hyperactivity and sleep defects in mutant flies. Fluorescence microscopy on mutant fly brains using Lysotracker dye revealed a significant increase in acidic compartments. Differentially expressed genes determined from RNA sequencing of fly brains are involved in biological processes that affect nervous system development. A genetic interaction network constructed using known interacting partners of these genes consists of 2 major subnetworks, one of which is enriched in genes associated with synaptic function and the other with neurodevelopmental processes. Our data indicate that lysosomal dysfunction arising from disruption of heparan sulfate breakdown has widespread effects on the steady state of intracellular vesicle transport, including vesicles associated with synaptic transmission. Evolutionary conservation of fundamental biological processes predicts that the Drosophila model of mucopolysaccharidosis type IIIB can serve as an in vivo system for the future development of therapies for mucopolysaccharidosis type IIIB and related disorders.

The BEN domain protein LIN-14 coordinates neuromuscular positioning during epidermal maturation

Development and function of an organism depend on coordinated inter-tissue interaction. How such interactions are maintained during tissue renewal and reorganization remains poorly understood. Here, we find that Caenorhabditis elegans BEN domain transcription factor LIN-14 is required in epidermis for maintaining the position of motor neurons and muscles during developmental tissue reorganization. lin-14 loss of function (lf) mutants display highly penetrant ventral neuromuscular mispositioning. These defects arise post-embryonically during first larval (L1) stage as the maturing epidermis replaces the embryonic ventral epidermis. Tissue-specific and temporally controlled depletion experiments indicate LIN-14 acts within the epidermis for ventral neuromuscular positioning. lin-14(lf) mutants show defects in formation of epidermis-muscle attachment complex hemidesmosomes in the maturing ventral epidermis, leading to detachment of muscles and motor neurons as well as movement defects. Our findings reveal a cell non-autonomous role for LIN-14 in coordinating inter-tissue interaction and neuromuscular positioning during epidermal maturation.

Chromosome structure in Drosophila is determined by boundary pairing not loop extrusion

Two different models have been proposed to explain how the endpoints of chromatin looped domains ('TADs') in eukaryotic chromosomes are determined. In the first, a cohesin complex extrudes a loop until it encounters a boundary element roadblock, generating a stem-loop. In this model, boundaries are functionally autonomous: they have an intrinsic ability to halt the movement of incoming cohesin complexes that is independent of the properties of neighboring boundaries. In the second, loops are generated by boundary:boundary pairing. In this model, boundaries are functionally non-autonomous, and their ability to form a loop depends upon how well they match with their neighbors. Moreover, unlike the loop-extrusion model, pairing interactions can generate both stem-loops and circle-loops. We have used a combination of MicroC to analyze how TADs are organized, and experimental manipulations of the even skipped TAD boundary, homie, to test the predictions of the 'loop-extrusion' and the 'boundary-pairing' models. Our findings are incompatible with the loop-extrusion model, and instead suggest that the endpoints of TADs in flies are determined by a mechanism in which boundary elements physically pair with their partners, either head-to-head or head-to-tail, with varying degrees of specificity. Although our experiments do not address how partners find each other, the mechanism is unlikely to require loop extrusion.

lncRNA read-through regulates the BX-C insulator Fub-1

Though long non-coding RNAs (lncRNAs) represent a substantial fraction of the Pol II transcripts in multicellular animals, only a few have known functions. Here we report that the blocking activity of the Bithorax complex (BX-C) Fub-1 boundary is segmentally regulated by its own lncRNA. The Fub-1 boundary is located between the Ultrabithorax (Ubx) gene and the bxd/pbx regulatory domain, which is responsible for regulating Ubx expression in parasegment PS6/segment A1. Fub-1 consists of two hypersensitive sites, HS1 and HS2. HS1 is an insulator while HS2 functions primarily as an lncRNA promoter. To activate Ubx expression in PS6/A1, enhancers in the bxd/pbx domain must be able to bypass Fub-1 blocking activity. We show that the expression of the Fub-1 lncRNAs in PS6/A1 from the HS2 promoter inactivates Fub-1 insulating activity. Inactivation is due to read-through as the HS2 promoter must be directed toward HS1 to disrupt blocking.

Unanticipated broad phylogeny of BEN DNA-binding domains revealed by structural homology searches

Organization of protein sequences into domain families is a foundation for cataloging and investigating protein functions. However, long-standing strategies based on primary amino acid sequences are blind to the possibility that proteins with dissimilar sequences could have comparable tertiary structures. Building on our recent findings that in silico structural predictions of BEN family DNA-binding domains closely resemble their experimentally determined crystal structures, we exploited the AlphaFold2 database for comprehensive identification of BEN domains. Indeed, we identified numerous novel BEN domains, including members of new subfamilies. For example, while no BEN domain factors had previously been annotated in C. elegans, this species actually encodes multiple BEN proteins. These include key developmental timing genes of orphan domain status, sel-7 and lin-14, the latter being the central target of the founding miRNA lin-4. We also reveal that the domain of unknown function 4806 (DUF4806), which is widely distributed across metazoans, is structurally similar to BEN and comprises a new subtype. Surprisingly, we find that BEN domains resemble both metazoan and non-metazoan homeodomains in 3D conformation and preserve characteristic residues, indicating that despite their inability to be aligned by conventional methods, these DNA-binding modules are probably evolutionarily related. Finally, we broaden the application of structural homology searches by revealing novel human members of DUF3504, which exists on diverse proteins with presumed or known nuclear functions. Overall, our work strongly expands this recently identified family of transcription factors and illustrates the value of 3D structural predictions to annotate protein domains and interpret their functions.

Stress-sensitive dynamics of miRNAs and Elba1 in Drosophila embryogenesis

Early-life stress can result in life-long effects that impact adult health and disease risk, but little is known about how such programming is established and maintained. Here, we show that such epigenetic memories can be initiated in the Drosophila embryo before the major wave of zygotic transcription, and higher-order chromatin structures are established. An early short heat shock results in elevated levels of maternal miRNA and reduced levels of a subgroup of zygotic genes in stage 5 embryos. Using a Dicer-1 mutant, we show that the stress-induced decrease in one of these genes, the insulator-binding factor Elba1, is dependent on functional miRNA biogenesis. Reduction in Elba1 correlates with the upregulation of early developmental genes and promotes a sustained weakening of heterochromatin in the adult fly as indicated by an increased expression of the PEV w^m4h reporter. We propose that maternal miRNAs, retained in response to an early embryonic heat shock, shape the subsequent de novo heterochromatin establishment that occurs during early development via direct or indirect regulation of some of the earliest expressed genes, including Elba1.

The metazoan landscape of mitochondrial DNA gene order and content is shaped by selection and affects mitochondrial transcription

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation.

Identification and characterization of BEND2 as a key regulator of meiosis during mouse spermatogenesis

The chromatin state, which undergoes global changes during spermatogenesis, is critical to meiotic initiation and progression. However, the key regulators involved and the underlying molecular mechanisms remain to be uncovered. Here, we report that mouse BEND2 is specifically expressed in spermatogenic cells around meiotic initiation and that it plays an essential role in meiotic progression. Bend2 gene knockout in male mice arrested meiosis at the transition from zygonema to pachynema, disrupted synapsis and DNA double-strand break repair, and induced nonhomologous chromosomal pairing. BEND2 interacted with chromatin-associated proteins that are components of certain transcription-repressor complexes. BEND2-binding sites were identified in diverse chromatin states and enriched in simple sequence repeats. BEND2 inhibited the expression of genes involved in meiotic initiation and regulated chromatin accessibility and the modification of H3K4me3. Therefore, our study identified BEND2 as a previously unknown key regulator of meiosis, gene expression, and chromatin state during mouse spermatogenesis.

Distinct structural bases for sequence-specific DNA binding by mammalian BEN domain proteins

The BEN domain is a recently recognized DNA binding module that is present in diverse metazoans and certain viruses. Several BEN domain factors are known as transcriptional repressors, but, overall, relatively little is known of how BEN factors identify their targets in humans. In particular, X-ray structures of BEN domain:DNA complexes are only known for Drosophila factors bearing a single BEN domain, which lack direct vertebrate orthologs. Here, we characterize several mammalian BEN domain (BD) factors, including from two NACC family BTB-BEN proteins and from BEND3, which has four BDs. In vitro selection data revealed sequence-specific binding activities of isolated BEN domains from all of these factors. We conducted detailed functional, genomic, and structural studies of BEND3. We show that BD4 is a major determinant for in vivo association and repression of endogenous BEND3 targets. We obtained a high-resolution structure of BEND3-BD4 bound to its preferred binding site, which reveals how BEND3 identifies cognate DNA targets and shows differences with one of its non-DNA-binding BEN domains (BD1). Finally, comparison with our previous invertebrate BEN structures, along with additional structural predictions using AlphaFold2 and RoseTTAFold, reveal distinct strategies for target DNA recognition by different types of BEN domain proteins. Together, these studies expand the DNA recognition activities of BEN factors and provide structural insights into sequence-specific DNA binding by mammalian BEN proteins.

An alternative D. melanogaster 7SK snRNP

BACKGROUND

The 7SK small nuclear RNA (snRNA) found in most metazoans is a key regulator of P-TEFb which in turn regulates RNA polymerase II elongation. Although its primary sequence varies in protostomes, its secondary structure and function are conserved across evolutionary distant taxa.

RESULTS

Here, we describe a novel ncRNA sharing many features characteristic of 7SK RNAs, in D. melanogaster. We examined the structure of the corresponding gene and determined the expression profiles of the encoded RNA, called snRNA:7SK:94F, during development. It is probably produced from the transcription of a lncRNA which is processed into a mature snRNA. We also addressed its biological function and we show that, like dm7SK, this alternative 7SK interacts in vivo with the different partners of the P-TEFb complex, i.e. HEXIM, LARP7 and Cyclin T. This novel RNA is widely expressed across tissues.

CONCLUSION

We propose that two distinct 7SK genes might contribute to the formation of the 7SK snRNP complex in D. melanogaster.

Drosophila architectural protein CTCF is not essential for fly survival and is able to function independently of CP190

CTCF is the most likely ancestor of proteins that contain large clusters of C2H2 zinc finger domains (C2H2) and is conserved among most bilateral organisms. In mammals, CTCF functions as the main architectural protein involved in the organization of topology-associated domains (TADs). In vertebrates and Drosophila, CTCF is involved in the regulation of homeotic genes. Previously, it was found that null mutations in the dCTCF gene died as pharate adults, which failed to eclose from their pupal case, or shortly after hatching of adults. Here, we obtained several new null dCTCF mutations and found that the complete inactivation of dCTCF appears is limited mainly to phenotypic manifestations of the Abd-B gene and fertility of adult flies. Many modifiers that are not associated with an independent phenotypic manifestation can significantly enhance the expressivity of the null dCTCF mutations, indicating that other architectural proteins are able to functionally compensate for dCTCF inactivation in Drosophila. We also mapped the 715-735 aa region of dCTCF as being essential for the interaction with the BTB (Broad-Complex, Tramtrack, and Bric a brac) and microtubule-targeting (M) domains of the CP190 protein, which binds to many architectural proteins. However, the mutational analysis showed that the interaction with CP190 was not important for the functional activity of dCTCF in vivo.

Mapping of functional elements of the Fab-6 boundary involved in the regulation of the Abd-B hox gene in Drosophila melanogaster

The autonomy of segment-specific regulatory domains in the Bithorax complex is conferred by boundary elements and associated Polycomb response elements (PREs). The Fab-6 boundary is located at the junction of the iab-5 and iab-6 domains. Previous studies mapped it to a nuclease hypersensitive region 1 (HS1), while the iab-6 PRE was mapped to a second hypersensitive region HS2 nearly 3 kb away. To analyze the role of HS1 and HS2 in boundary we generated deletions of HS1 or HS1 + HS2 that have attP site for boundary replacement experiments. The 1389 bp HS1 deletion can be rescued by a 529 bp core Fab-6 sequence that includes two CTCF sites. However, Fab-6 HS1 cannot rescue the HS1 + HS2 deletion or substitute for another BX-C boundary - Fab-7. For this it must be combined with a PRE, either Fab-7 HS3, or Fab-6 HS2. These findings suggest that the boundary function of Fab-6 HS1 must be bolstered by a second element that has PRE activity.