Dissecting the regulatory landscape of the Abd-B gene of the bithorax complex

Boundary bypass activity in the abdominal-B region of the Drosophila bithorax complex is position dependent and regulated

Expression of Abdominal-B (Abd-B) in abdominal segments A5-A8 is controlled by four regulatory domains, iab-5-iab-8. Each domain has an initiator element (which sets the activity state), elements that maintain this state and tissue-specific enhancers. To ensure their functional autonomy, each domain is bracketed by boundary elements (Mcp, Fab-7, Fab-7 and Fab-8). In addition to blocking crosstalk between adjacent regulatory domains, the Fab boundaries must also have bypass activity so the relevant regulatory domains can 'jump over' intervening boundaries and activate the Abd-B promoter. In the studies reported here we have investigated the parameters governing bypass activity. We find that the bypass elements in the Fab-7 and Fab-8 boundaries must be located in the regulatory domain that is responsible for driving Abd-B expression. We suggest that bypass activity may also be subject to regulation.

Boundary Bypass Activity in the Abdominal-B Region of the Drosophila Bithorax Complex is Position Dependent and Regulated

Expression of Abdominal-B ( Abd-B ) in abdominal segments A5 - A8 is controlled by four regulatory domains, iab-5 - iab-8 . Each domain has an initiator element (which sets the activity state), elements that maintain this state and tissue-specific enhancers. To ensure their functional autonomy, each domain is bracketed by boundary elements ( Mcp , Fab-7 , Fab-7 and Fab-8 ). In addition to blocking crosstalk between adjacent regulatory domains, the Fab boundaries must also have bypass activity so the relevant regulatory domains can "jump over" intervening boundaries and activate the Abd-B promoter. In the studies reported here we have investigated the parameters governing bypass activity. We find that the bypass elements in the Fab-7 and Fab-8 boundaries must be located in the regulatory domain that is responsible for driving Abd-B expression. We suggest that bypass activity may also be subject to regulation.

Summary Statement

Boundaries separating Abd-B regulatory domains block crosstalk between domains and mediate their interactions with Abd-B . The latter function is location but not orientation dependent.

Mechanisms of Interaction between Enhancers and Promoters in Three Drosophila Model Systems

In higher eukaryotes, the regulation of developmental gene expression is determined by enhancers, which are often located at a large distance from the promoters they regulate. Therefore, the architecture of chromosomes and the mechanisms that determine the functional interaction between enhancers and promoters are of decisive importance in the development of organisms. Mammals and the model animal Drosophila have homologous key architectural proteins and similar mechanisms in the organization of chromosome architecture. This review describes the current progress in understanding the mechanisms of the formation and regulation of long-range interactions between enhancers and promoters at three well-studied key regulatory loci in Drosophila.

CRISPR/Cas9 and FLP-FRT mediated regulatory dissection of the BX-C of Drosophila melanogaster

The homeotic genes or Hox define the anterior-posterior (AP) body axis formation in bilaterians and are often present on the chromosome in an order collinear to their function across the AP axis. However, there are many cases wherein the Hox are not collinear, but their expression pattern is conserved across the AP axis. The expression pattern of Hox is attributed to the cis-regulatory modules (CRMs) consisting of enhancers, initiators, or repressor elements that regulate the genes in a segment-specific manner. In the Drosophila melanogaster Hox complex, the bithorax complex (BX-C) and even the CRMs are organized in an order that is collinear to their function in the thoracic and abdominal segments. In the present study, the regulatorily inert regions were targeted using CRISPR/Cas9 to generate a series of transgenic lines with the insertion of FRT sequences. These FRT lines are repurposed to shuffle the CRMs associated with Abd-B to generate modular deletion, duplication, or inversion of multiple CRMs. The rearrangements yielded entirely novel phenotypes in the fly suggesting the requirement of such complex manipulations to address the significance of higher order arrangement of the CRMs. The functional map and the transgenic flies generated in this study are important resources to decipher the collective ability of multiple regulatory elements in the eukaryotic genome to function as complex modules.

The Drosophila Fab-7 boundary modulates Abd-B gene activity by guiding an inversion of collinear chromatin organization and alternate promoter use

Hox genes encode transcription factors that specify segmental identities along the anteroposterior body axis. These genes are organized in clusters, where their order corresponds to their activity along the body axis, a feature known as collinearity. In Drosophila, the BX-C cluster contains the three most posterior Hox genes, where their collinear activation incorporates progressive changes in histone modifications, chromatin architecture, and use of boundary elements and cis-regulatory regions. To dissect functional hierarchies, we compare chromatin organization in cell lines and larvae, with a focus on the Abd-B gene. Our work establishes the importance of the Fab-7 boundary for insulation between 3D domains carrying different histone modifications. Interestingly, we detect a non-canonical inversion of collinear chromatin dynamics at Abd-B, with the domain of active histone modifications progressively decreasing in size. This dynamic chromatin organization differentially activates the alternative promoters of the Abd-B gene, thereby expanding the possibilities for fine-tuning of transcriptional output.

Essential role of Cp190 in physical and regulatory boundary formation

Boundaries in animal genomes delimit contact domains with enhanced internal contact frequencies and have debated functions in limiting regulatory cross-talk between domains and guiding enhancers to target promoters. Most mammalian boundaries form by stalling of chromosomal loop-extruding cohesin by CTCF, but most Drosophila boundaries form CTCF independently. However, how CTCF-independent boundaries form and function remains largely unexplored. Here, we assess genome folding and developmental gene expression in fly embryos lacking the ubiquitous boundary-associated factor Cp190. We find that sequence-specific DNA binding proteins such as CTCF and Su(Hw) directly interact with and recruit Cp190 to form most promoter-distal boundaries. Cp190 is essential for early development and prevents regulatory cross-talk between specific gene loci that pattern the embryo. Cp190 was, in contrast, dispensable for long-range enhancer-promoter communication at tested loci. Cp190 is thus currently the major player in fly boundary formation and function, revealing that diverse mechanisms evolved to partition genomes into independent regulatory domains.

Redundant enhancers in the iab-5 domain cooperatively activate Abd-B in the A5 and A6 abdominal segments of Drosophila

The Abdominal-B (Abd-B) gene belongs to the bithorax complex and its expression is controlled by four regulatory domains, iab-5, iab-6, iab-7 and iab-8, each of which is thought to be responsible for directing the expression of Abd-B in one of the abdominal segments from A5 to A8. A variety of experiments have supported the idea that BX-C regulatory domains are functionally autonomous and that each domain is both necessary and sufficient to orchestrate the development of the segment they specify. Unexpectedly, we discovered that this model does not always hold. Instead, we find that tissue-specific enhancers located in the iab-5 domain are required for the proper activation of Abd-B not only in A5 but also in A6. Our findings indicate that the functioning of the iab-5 and iab-6 domains in development of the adult cuticle A5 and A6 in males fit better with an additive model, much like that first envisioned by Ed Lewis.

Genomic organization of the autonomous regulatory domain of eyeless locus in Drosophila melanogaster

In Drosophila, expression of eyeless (ey) gene is restricted to the developing eyes and central nervous system. However, the flanking genes, myoglianin (myo), and bent (bt) have different temporal and spatial expression patterns as compared to the ey. How distinct regulation of ey is maintained is mostly unknown. Earlier, we have identified a boundary element intervening myo and ey genes (ME boundary) that prevents the crosstalk between the cis-regulatory elements of myo and ey genes. In the present study, we further searched for the cis-elements that define the domain of ey and maintain its expression pattern. We identify another boundary element between ey and bt, the EB boundary. The EB boundary separates the regulatory landscapes of ey and bt genes. The two boundaries, ME and EB, show a long-range interaction as well as interact with the nuclear architecture. This suggests functional autonomy of the ey locus and its insulation from differentially regulated flanking regions. We also identify a new Polycomb Response Element, the ey-PRE, within the ey domain. The expression state of the ey gene, once established during early development is likely to be maintained with the help of ey-PRE. Our study proposes a general regulatory mechanism by which a gene can be maintained in a functionally independent chromatin domain in gene-rich euchromatin.

Homeotic Genes: Clustering, Modularity, and Diversity

Hox genes code for transcription factors and are evolutionarily conserved. They regulate a plethora of downstream targets to define the anterior-posterior (AP) body axis of a developing bilaterian embryo. Early work suggested a possible role of clustering and ordering of Hox to regulate their expression in a spatially restricted manner along the AP axis. However, the recent availability of many genome assemblies for different organisms uncovered several examples that defy this constraint. With recent advancements in genomics, the current review discusses the arrangement of Hox in various organisms. Further, we revisit their discovery and regulation in Drosophila melanogaster. We also review their regulation in different arthropods and vertebrates, with a significant focus on Hox expression in the crustacean Parahyale hawaiensis. It is noteworthy that subtle changes in the levels of Hox gene expression can contribute to the development of novel features in an organism. We, therefore, delve into the distinct regulation of these genes during primary axis formation, segment identity, and extra-embryonic roles such as in the formation of hair follicles or misregulation leading to cancer. Toward the end of each section, we emphasize the possibilities of several experiments involving various organisms, owing to the advancements in the field of genomics and CRISPR-based genome engineering. Overall, we present a holistic view of the functioning of Hox in the animal world.

Redundant enhancers in the iab-5 domain cooperatively activate Abd-B in the A5 and A6 abdominal segments of Drosophila.. [DOI: 10.1101/2021.05.22.445252] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The homeotic Abdominal-B (Abd-B) gene belongs to Bithorax complex and is regulated by four regulatory domains named iab-5, iab-6, iab-7 and iab-8, each of which is thought to be responsible for directing the expression of Abd-B in one of the abdominal segments from A5 to A8. It is assumed that male specific features of the adult cuticle in A5 is solely dependent on regulatory elements located in iab-5, while the regulatory elements in the iab-6 are both necessary and sufficient for the proper differentiation of the A6 cuticle. Unexpectedly, we found that this long held assumption is not correct. Instead, redundant tissue-specific enhancers located in the iab-5 domain are required for the proper activation of Abd-B not only in A5 but also in A6. Our study of deletions shows that the iab-5 initiator is essential for the functioning of the iab-5 enhancers in A5, as well as for the correct differentiation of A6. This requirement is circumvented by deletions that remove the initiator and most of the iab-5 regulatory domain sequences. While the remaining iab-5 enhancers are inactive in A5, they are activated in A6 and contribute to the differentiation of this segment. In this case, Abd-B stimulation by the iab-5 enhancers in A6 depends on the initiators in the iab-4 and iab-6 domains.

Summary Statement

In Drosophila, the segmental-specific expression of the homeotic gene Abdominal-B in the abdominal segments is regulated by autonomous regulatory domains. We demonstrated cooperation between these domains in activation of Abdominal-B.

Mechanism and functional role of the interaction between CP190 and the architectural protein Pita in Drosophila melanogaster

Background

Pita is required for Drosophila development and binds specifically to a long motif in active promoters and insulators. Pita belongs to the Drosophila family of zinc-finger architectural proteins, which also includes Su(Hw) and the conserved among higher eukaryotes CTCF. The architectural proteins maintain the active state of regulatory elements and the long-distance interactions between them. In particular, Pita is involved in the formation of several boundaries between regulatory domains that controlled the expression of three hox genes in the Bithorax complex (BX-C). The CP190 protein is recruited to chromatin through interaction with the architectural proteins.

Results

Using in vitro pull-down analysis, we precisely mapped two unstructured regions of Pita that interact with the BTB domain of CP190. Then we constructed transgenic lines expressing the Pita protein of the wild-type and mutant variants lacking CP190-interacting regions. We have demonstrated that CP190-interacting region of the Pita can maintain nucleosome-free open chromatin and is critical for Pita-mediated enhancer blocking activity in BX-C. At the same time, interaction with CP190 is not required for the in vivo function of the mutant Pita protein, which binds to the same regions of the genome as the wild-type protein. Unexpectedly, we found that CP190 was still associated with the most of genome regions bound by the mutant Pita protein, which suggested that other architectural proteins were continuing to recruit CP190 to these regions.

Conclusions

The results directly demonstrate role of CP190 in insulation and support a model in which the regulatory elements are composed of combinations of binding sites that interact with several architectural proteins with similar functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00391-x.

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.

GAGA factor: a multifunctional pioneering chromatin protein

The Drosophila GAGA factor (GAF) is a multifunctional protein implicated in nucleosome organization and remodeling, activation and repression of gene expression, long distance enhancer-promoter communication, higher order chromosome structure, and mitosis. This broad range of activities poses questions about how a single protein can perform so many seemingly different and unrelated functions. Current studies argue that GAF acts as a "pioneer" factor, generating nucleosome-free regions of chromatin for different classes of regulatory elements. The removal of nucleosomes from regulatory elements in turn enables other factors to bind to these elements and carry out their specialized functions. Consistent with this view, GAF associates with a collection of chromatin remodelers and also interacts with proteins implicated in different regulatory functions. In this review, we summarize the known activities of GAF and the functions of its protein partners.

Molecular Lessons from the Drosophila Bithorax Complex

The Genetics Society of America's (GSA's) Edward Novitski Prize recognizes a single experimental accomplishment or a body of work in which an exceptional level of creativity, and intellectual ingenuity, has been used to design and execute scientific experiments to solve a difficult problem in genetics. The 2020 recipient is Welcome W. Bender of Harvard Medical School, recognizing his creativity and ingenuity in revealing the molecular nature and regulation of the bithorax gene complex.

Distinct Elements Confer the Blocking and Bypass Functions of the Bithorax Fab-8 Boundary

Boundaries in the Drosophila bithorax complex (BX-C) enable the regulatory domains that drive parasegment-specific expression of the three Hox genes to function autonomously. The four regulatory domains (iab-5, iab-6, iab-7, and iab-8) that control the expression of the Abdominal-B (Abd-B) gene are located downstream of the transcription unit, and are delimited by the Mcp, Fab-6, Fab-7, and Fab-8 boundaries. These boundaries function to block cross talk between neighboring regulatory domains. In addition, three of the boundaries (Fab-6, Fab-7, and Fab-8) must also have bypass activity so that regulatory domains distal to the boundaries can contact the Abd-B promoter. In the studies reported here, we have undertaken a functional dissection of the Fab-8 boundary using a boundary-replacement strategy. Our studies indicate that the Fab-8 boundary has two separable subelements. The distal subelement blocks cross talk, but cannot support bypass. The proximal subelement has only minimal blocking activity but is able to mediate bypass. A large multiprotein complex, the LBC (large boundary complex), binds to sequences in the proximal subelement and contributes to its bypass activity. The same LBC complex has been implicated in the bypass activity of the Fab-7 boundary.

Separate Polycomb Response Elements control chromatin state and activation of the vestigial gene

Patterned expression of many developmental genes is specified by transcription factor gene expression, but is thought to be refined by chromatin-mediated repression. Regulatory DNA sequences called Polycomb Response Elements (PREs) are required to repress some developmental target genes, and are widespread in genomes, suggesting that they broadly affect developmental programs. While PREs in transgenes can nucleate trimethylation on lysine 27 of the histone H3 tail (H3K27me3), none have been demonstrated to be necessary at endogenous chromatin domains. This failure is thought to be due to the fact that most endogenous H3K27me3 domains contain many PREs, and individual PREs may be redundant. In contrast to these ideas, we show here that PREs near the wing selector gene vestigial have distinctive roles at their endogenous locus, even though both PREs are repressors in transgenes. First, a PRE near the promoter is required for vestigial activation and not for repression. Second, only the distal PRE contributes to H3K27me3, but even removal of both PREs does not eliminate H3K27me3 across the vestigial domain. Thus, endogenous chromatin domains appear to be intrinsically marked by H3K27me3, and PREs appear required to enhance this chromatin modification to high levels at inactive genes.

Changes throughout a Genetic Network Mask the Contribution of Hox Gene Evolution

Hox genes pattern the anterior-posterior axis of animals and are posited to drive animal body plan evolution, yet their precise role in evolution has been difficult to determine. Here, we identified evolutionary modifications in the Hox gene Abd-B that dramatically altered its expression along the body plan of Drosophila santomea. Abd-B is required for pigmentation in Drosophila yakuba, the sister species of D. santomea, and changes to Abd-B expression would be predicted to make large contributions to the loss of body pigmentation in D. santomea. However, manipulating Abd-B expression in current-day D. santomea does not affect pigmentation. We attribute this epistatic interaction to four other genes within the D. santomea pigmentation network, three of which have evolved expression patterns that do not respond to Abd-B. Our results demonstrate how body plans may evolve through small evolutionary steps distributed throughout Hox-regulated networks. Polygenicity and epistasis may hinder efforts to identify genes and mechanisms underlying macroevolutionary traits.

Functional dissection of the developmentally restricted BEN domain chromatin boundary factor Insensitive

Background

Boundaries in the Drosophila bithorax complex delimit autonomous regulatory domains that activate the parasegment (PS)-specific expression of homeotic genes. The Fab-7 boundary separates the iab-6 and iab-7 regulatory domains that control Abd-B expression in PS11 and PS12. This boundary is composed of multiple functionally redundant elements and has two key activities: it blocks crosstalk between iab-6 and iab-7 and facilitates boundary bypass.

Results

Here, we have used a structure–function approach to elucidate the biochemical properties and the in vivo activities of a conserved BEN domain protein, Insensitive, that is associated with Fab-7. Our biochemical studies indicate that in addition to the C-terminal BEN DNA-binding domain, Insv has two domains that mediate multimerization: one is a coiled-coil domain in the N-terminus, and the other is next to the BEN domain. These multimerization domains enable Insv to bind simultaneously to two canonical 8-bp recognition motifs, as well as to a ~ 100-bp non-canonical recognition sequence. They also mediate the assembly of higher-order multimers in the presence of DNA. Transgenic proteins lacking the N-terminal coiled-coil domain are compromised for boundary function in vivo. We also show that Insv interacts directly with CP190, a protein previously implicated in the boundary functions of several DNA-binding proteins, including Su(Hw) and dCTCF. While CP190 interaction is required for Insv binding to a subset of sites on polytene chromosomes, it has only a minor role in the boundary activity of Insv in the context of Fab-7.

Conclusions

The subdivision of eukaryotic chromosomes into discrete topological domains depends upon the pairing of boundary elements. In flies, pairing interactions are specific and typically orientation dependent. They occur in cis between neighboring heterologous boundaries, and in trans between homologous boundaries. One potential mechanism for ensuring pairing-interaction specificity is the use of sequence-specific DNA-binding proteins that can bind simultaneously with two or more recognition sequences. Our studies indicate that Insv can assemble into a multivalent DNA-binding complex and that the N-terminal Insv multimerization domain is critical for boundary function.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0249-2) contains supplementary material, which is available to authorized users.

Boundaries mediate long-distance interactions between enhancers and promoters in the Drosophila Bithorax complex

Drosophila bithorax complex (BX-C) is one of the best model systems for studying the role of boundaries (insulators) in gene regulation. Expression of three homeotic genes, Ubx, abd-A, and Abd-B, is orchestrated by nine parasegment-specific regulatory domains. These domains are flanked by boundary elements, which function to block crosstalk between adjacent domains, ensuring that they can act autonomously. Paradoxically, seven of the BX-C regulatory domains are separated from their gene target by at least one boundary, and must “jump over” the intervening boundaries. To understand the jumping mechanism, the Mcp boundary was replaced with Fab-7 and Fab-8. Mcp is located between the iab-4 and iab-5 domains, and defines the border between the set of regulatory domains controlling abd-A and Abd-B. When Mcp is replaced by Fab-7 or Fab-8, they direct the iab-4 domain (which regulates abd-A) to inappropriately activate Abd-B in abdominal segment A4. For the Fab-8 replacement, ectopic induction was only observed when it was inserted in the same orientation as the endogenous Fab-8 boundary. A similar orientation dependence for bypass activity was observed when Fab-7 was replaced by Fab-8. Thus, boundaries perform two opposite functions in the context of BX-C–they block crosstalk between neighboring regulatory domains, but at the same time actively facilitate long distance communication between the regulatory domains and their respective target genes.

Drosophila bithorax complex (BX-C) is one of a few examples demonstrating in vivo role of boundary/insulator elements in organization of independent chromatin domains. BX-C contains three HOX genes, whose parasegment-specific pattern is controlled by cis-regulatory domains flanked by boundary/insulator elements. Since the boundaries ensure autonomy of adjacent domains, the presence of these elements poses a paradox: how do the domains bypass the intervening boundaries and contact their proper regulatory targets? According to the textbook model, BX-C regulatory domains are able to bypass boundaries because they harbor special promoter targeting sequences. However, contrary to this model, we show here that the boundaries themselves play an active role in directing regulatory domains to their appropriate HOX gene promoter.

The bithorax complex iab-7 Polycomb response element has a novel role in the functioning of the Fab-7 chromatin boundary

Expression of the three bithorax complex homeotic genes is orchestrated by nine parasegment-specific regulatory domains. Autonomy of each domain is conferred by boundary elements (insulators). Here, we have used an in situ replacement strategy to reanalyze the sequences required for the functioning of one of the best-characterized fly boundaries, Fab-7. It was initially identified by a deletion, Fab-7¹, that transformed parasegment (PS) 11 into a duplicate copy of PS12. Fab-7¹ deleted four nuclease hypersensitive sites, HS*, HS1, HS2, and HS3, located between the iab-6 and iab-7 regulatory domains. Transgenic and P-element excision experiments mapped the boundary to HS*+HS1+HS2, while HS3 was shown to be the iab-7 Polycomb response element (PRE). Recent replacement experiments showed that HS1 is both necessary and sufficient for boundary activity when HS3 is also present in the replacement construct. Surprisingly, while HS1+HS3 combination has full boundary activity, we discovered that HS1 alone has only minimal function. Moreover, when combined with HS3, only the distal half of HS1, dHS1, is needed. A ~1,000 kD multiprotein complex containing the GAF protein, called the LBC, binds to the dHS1 sequence and we show that mutations in dHS1, that disrupt LBC binding in nuclear extracts, eliminate boundary activity and GAF binding in vivo. HS3 has binding sites for GAF and Pho proteins that are required for PRE silencing. In contrast, HS3 boundary activity only requires the GAF binding sites. LBC binding with HS3 in nuclear extracts, and GAF association in vivo, depend upon the HS3 GAF sites, but not the Pho sites. Consistent with a role for the LBC in HS3 boundary activity, the boundary function of the dHS1+HS3^mPho combination is lost when the flies are heterozygous for a mutation in the GAF gene. Taken together, these results reveal a novel function for the iab-7 PREs in chromosome architecture.

Polycomb group proteins (PcG) are important epigenetic regulators of developmental genes in all higher eukaryotes. In Drosophila, these proteins are bound to specific regulatory DNA elements called Polycomb group Response Elements (PREs). Drosophila PREs are made up of binding sites for a complex array of DNA binding proteins, including GAF and Pho. In the regulatory region of the bithorax complex (BX-C), the boundary/insulator elements organize the autonomous regulatory domains, and their active or repressed states are regulated by PREs. Here, we studied functional properties of sequences that constitute the Fab-7 boundary and the adjacent iab-7 PRE. It was previously thought that the sole function of the iab-7 PRE is to recruit PcG proteins in parasegments anterior to PS12 and silence the iab-7 domain. However, we found that the iab-7 PRE also functions as a component of the Fab-7 boundary. The boundary activity of the iab-7 PRE sequence depends upon a large complex called the LBC. We show that it is possible to reconstitute a fully functional boundary by combining the LBC binding sequences in HS1 with the iab-7 PRE. Moreover, its boundary function is independent of its PcG silencing activity.

The lncRNA male-specific abdominal plays a critical role in Drosophila accessory gland development and male fertility

Although thousands of long non-coding RNAs (lncRNA) have been identified in the genomes of higher eukaryotes, the precise function of most of them is still unclear. Here, we show that a >65 kb, male-specific, lncRNA, called male-specific abdominal (msa) is required for the development of the secondary cells of the Drosophila male accessory gland (AG). msa is transcribed from within the Drosophila bithorax complex and shares much of its sequence with another lncRNA, the iab-8 lncRNA, which is involved in the development of the central nervous system (CNS). Both lncRNAs perform much of their functions via a shared miRNA embedded within their sequences. Loss of msa, or of the miRNA it contains, causes defects in secondary cell morphology and reduces male fertility. Although both lncRNAs express the same miRNA, the phenotype in the secondary cells and the CNS seem to reflect misregulation of different targets in the two tissues.

In many animals, the male seminal fluid induces physiology changes in the mated female that increase a male’s reproductive success. These changes are often referred to as the post-mating response (PMR). In Drosophila, the seminal fluid proteins responsible for generating the PMR are made in a specialized gland, analogous to the mammalian seminal vesicle and prostate, called the accessory gland (AG). In this work, we show that a male-specific, long, non-coding RNA (lncRNA), called msa, plays a critical role in the development and function of this gland, primarily through a microRNA (miRNA) encoded within its sequence. This same miRNA had previously been shown to be expressed in the central nervous system (CNS) via an alternative promoter, where its ability to repress homeotic genes is required for both male and female fertility. Here, we present evidence that the targets of this miRNA in the AG are likely different from those found in the CNS. Thus, the same miRNA seems to have been selected to affect Drosophila fertility through two different mechanisms. Although many non-coding RNAs have now been identified, very few can be shown to have function. Our work highlights a lncRNA that has multiple biological functions, affecting cellular morphology and fertility.

Architectural protein Pita cooperates with dCTCF in organization of functional boundaries in Bithorax complex

Boundaries in the Bithorax complex (BX-C) of Drosophila delimit autonomous regulatory domains that drive parasegment-specific expression of homeotic genes. BX-C boundaries have two crucial functions: they must block crosstalk between adjacent regulatory domains and at the same time facilitate boundary bypass. The C2H2 zinc-finger protein Pita binds to several BX-C boundaries, including Fab-7 and Mcp To study Pita functions, we have used a boundary replacement strategy by substituting modified DNAs for the Fab-7 boundary, which is located between the iab-6 and iab-7 regulatory domains. Multimerized Pita sites block iab-6↔iab-7 crosstalk but fail to support iab-6 regulation of Abd-B (bypass). In the case of Fab-7, we used a novel sensitized background to show that the two Pita-binding sites contribute to its boundary function. Although Mcp is from BX-C, it does not function appropriately when substituted for Fab-7: it blocks crosstalk but does not support bypass. Mutation of the Mcp Pita site disrupts blocking activity and also eliminates dCTCF binding. In contrast, mutation of the Mcp dCTCF site does not affect Pita binding, and this mutant boundary retains partial function.

Different Evolutionary Strategies To Conserve Chromatin Boundary Function in the Bithorax Complex

Chromatin boundary elements subdivide chromosomes in multicellular organisms into physically independent domains. In addition to this architectural function, these elements also play a critical role in gene regulation. Here we investigated the evolution of a Drosophila Bithorax complex boundary element called Fab-7, which is required for the proper parasegment specific expression of the homeotic Abd-B gene. Using a “gene” replacement strategy, we show that Fab-7 boundaries from two closely related species, D. erecta and D. yakuba, and a more distant species, D. pseudoobscura, are able to substitute for the melanogaster boundary. Consistent with this functional conservation, the two known Fab-7 boundary factors, Elba and LBC, have recognition sequences in the boundaries from all species. However, the strategies used for maintaining binding and function in the face of sequence divergence is different. The first is conventional, and depends upon conservation of the 8 bp Elba recognition sequence. The second is unconventional, and takes advantage of the unusually large and flexible sequence recognition properties of the LBC boundary factor, and the deployment of multiple LBC recognition elements in each boundary. In the former case, binding is lost when the recognition sequence is altered. In the latter case, sequence divergence is accompanied by changes in the number, relative affinity, and location of the LBC recognition elements.

Ever-Changing Landscapes: Transcriptional Enhancers in Development and Evolution

A class of cis-regulatory elements, called enhancers, play a central role in orchestrating spatiotemporally precise gene-expression programs during development. Consequently, divergence in enhancer sequence and activity is thought to be an important mediator of inter- and intra-species phenotypic variation. Here, we give an overview of emerging principles of enhancer function, current models of enhancer architecture, genomic substrates from which enhancers emerge during evolution, and the influence of three-dimensional genome organization on long-range gene regulation. We discuss intricate relationships between distinct elements within complex regulatory landscapes and consider their potential impact on specificity and robustness of transcriptional regulation.

Formation of a Polycomb-Domain in the Absence of Strong Polycomb Response Elements

Polycomb group response elements (PREs) in Drosophila are DNA-elements that recruit Polycomb proteins (PcG) to chromatin and regulate gene expression. PREs are easily recognizable in the Drosophila genome as strong peaks of PcG-protein binding over discrete DNA fragments; many small but statistically significant PcG peaks are also observed in PcG domains. Surprisingly, in vivo deletion of the four characterized strong PREs from the PcG regulated invected-engrailed (inv-en) gene complex did not disrupt the formation of the H3K27me3 domain and did not affect inv-en expression in embryos or larvae suggesting the presence of redundant PcG recruitment mechanism. Further, the 3D-structure of the inv-en domain was only minimally altered by the deletion of the strong PREs. A reporter construct containing a 7.5kb en fragment that contains three weak peaks but no large PcG peaks forms an H3K27me3 domain and is PcG-regulated. Our data suggests a model for the recruitment of PcG-complexes to Drosophila genes via interactions with multiple, weak PREs spread throughout an H3K27me3 domain.

Functional Dissection of the Blocking and Bypass Activities of the Fab-8 Boundary in the Drosophila Bithorax Complex

Functionally autonomous regulatory domains direct the parasegment-specific expression of the Drosophila Bithorax complex (BX-C) homeotic genes. Autonomy is conferred by boundary/insulator elements that separate each regulatory domain from its neighbors. For six of the nine parasegment (PS) regulatory domains in the complex, at least one boundary is located between the domain and its target homeotic gene. Consequently, BX-C boundaries must not only block adventitious interactions between neighboring regulatory domains, but also be permissive (bypass) for regulatory interactions between the domains and their gene targets. To elucidate how the BX-C boundaries combine these two contradictory activities, we have used a boundary replacement strategy. We show that a 337 bp fragment spanning the Fab-8 boundary nuclease hypersensitive site and lacking all but 83 bp of the 625 bp Fab-8 PTS (promoter targeting sequence) fully rescues a Fab-7 deletion. It blocks crosstalk between the iab-6 and iab-7 regulatory domains, and has bypass activity that enables the two downstream domains, iab-5 and iab-6, to regulate Abdominal-B (Abd-B) transcription in spite of two intervening boundary elements. Fab-8 has two dCTCF sites and we show that they are necessary both for blocking and bypass activity. However, CTCF sites on their own are not sufficient for bypass. While multimerized dCTCF (or Su(Hw)) sites have blocking activity, they fail to support bypass. Moreover, this bypass defect is not rescued by the full length PTS. Finally, we show that orientation is critical for the proper functioning the Fab-8 replacement. Though the inverted Fab-8 boundary still blocks crosstalk, it disrupts the topology of the Abd-B regulatory domains and does not support bypass. Importantly, altering the orientation of the Fab-8 dCTCF sites is not sufficient to disrupt bypass, indicating that orientation dependence is conferred by other factors.

Boundary elements in the Bithorax complex have two seemingly contradictory activities. They must block crosstalk between neighboring regulatory domains, but at the same time be permissive (insulator bypass) for regulatory interactions between the domains and the BX-C homeotic genes. We have used a replacement strategy to investigate how they carry out these two functions. We show that a 337 bp fragment spanning the Fab-8 boundary nuclease hypersensitive site is sufficient to fully rescue a Fab-7 boundary deletion. It blocks crosstalk and supports bypass. As has been observed in transgene assays, blocking activity requires the Fab-8 dCTCF sites, while full bypass activity requires the dCTCF sites plus a small part of PTS. In transgene assays, bypass activity typically depends on the orientation of the two insulators relative to each other. A similar orientation dependence is observed for the Fab-8 replacement in BX-C. When the orientation of the Fab-8 boundary is reversed, bypass activity is lost, while blocking is unaffected. Interestingly, unlike what has been observed in mammals, reversing the orientation of only the Fab-8 dCTCF sites does not affect boundary function. This finding indicates that other Fab-8 factors must play a critical role in determining orientation. Taken together, our findings argue that carrying out the paradoxical functions of the BX-C boundaries does not require any unusual or special properties; rather BX-C boundaries utilize generic blocking and insulator bypass activities that are appropriately adapted to their regulatory context. Thus making them a good model for studying the functional properties of boundaries/insulators in their native setting.

The boundary paradox in the Bithorax complex

The parasegment-specific expression of the three Drosophila Bithorax complex homeotic genes is orchestrated by nine functionally autonomous regulatory domains. Functional autonomy depends upon special elements called boundaries or insulators that are located between each domain. The boundaries ensure the independent activity of each domain by blocking adventitious interactions with initiators, enhancers and silencers in the neighboring domains. However, this blocking activity poses a regulatory paradox--the Bithorax boundaries are also able to insulate promoters from regulatory interactions with enhancers and silencers and six of the nine Bithorax regulatory domains are separated from their target genes by at least one boundary element. Here we consider several mechanisms that have been suggested for how the Bithorax regulatory domains are able to bypass intervening boundary elements and direct the appropriate parasegment-specific temporal and spatial expression of their target gene.

The open for business model of the bithorax complex in Drosophila

After nearly 30 years of effort, Ed Lewis published his 1978 landmark paper in which he described the analysis of a series of mutations that affect the identity of the segments that form along the anterior-posterior (AP) axis of the fly (Lewis 1978). The mutations behaved in a non-canonical fashion in complementation tests, forming what Ed Lewis called a "pseudo-allelic" series. Because of this, he never thought that the mutations represented segment-specific genes. As all of these mutations were grouped to a particular area of the Drosophila third chromosome, the locus became known of as the bithorax complex (BX-C). One of the key findings of Lewis' article was that it revealed for the first time, to a wide scientific audience, that there was a remarkable correlation between the order of the segment-specific mutations along the chromosome and the order of the segments they affected along the AP axis. In Ed Lewis' eyes, the mutants he discovered affected "segment-specific functions" that were sequentially activated along the chromosome as one moves from anterior to posterior along the body axis (the colinearity concept now cited in elementary biology textbooks). The nature of the "segment-specific functions" started to become clear when the BX-C was cloned through the pioneering chromosomal walk initiated in the mid 1980s by the Hogness and Bender laboratories (Bender et al. 1983a; Karch et al. 1985). Through this molecular biology effort, and along with genetic characterizations performed by Gines Morata's group in Madrid (Sanchez-Herrero et al. 1985) and Robert Whittle's in Sussex (Tiong et al. 1985), it soon became clear that the whole BX-C encoded only three protein-coding genes (Ubx, abd-A, and Abd-B). Later, immunostaining against the Ubx protein hinted that the segment-specific functions could, in fact, be cis-regulatory elements regulating the expression of the three protein-coding genes. In 1987, Peifer, Karch, and Bender proposed a comprehensive model of the functioning of the BX-C, in which the "segment-specific functions" appear as segment-specific enhancers regulating, Ubx, abd-A, or Abd-B (Peifer et al. 1987). Key to their model was that the segmental address of these enhancers was not an inherent ability of the enhancers themselves, but was determined by the chromosomal location in which they lay. In their view, the sequential activation of the segment-specific functions resulted from the sequential opening of chromatin domains along the chromosome as one moves from anterior to posterior. This model soon became known of as the open for business model. While the open for business model is quite easy to visualize at a conceptual level, molecular evidence to validate this model has been missing for almost 30 years. The recent publication describing the outstanding, joint effort from the Bender and Kingston laboratories now provides the missing proof to support this model (Bowman et al. 2014). The purpose of this article is to review the open for business model and take the reader through the genetic arguments that led to its elaboration.

Functional Requirements for Fab-7 Boundary Activity in the Bithorax Complex

Chromatin boundaries are architectural elements that determine the three-dimensional folding of the chromatin fiber and organize the chromosome into independent units of genetic activity. The Fab-7 boundary from the Drosophila bithorax complex (BX-C) is required for the parasegment-specific expression of the Abd-B gene. We have used a replacement strategy to identify sequences that are necessary and sufficient for Fab-7 boundary function in the BX-C. Fab-7 boundary activity is known to depend on factors that are stage specific, and we describe a novel ∼700-kDa complex, the late boundary complex (LBC), that binds to Fab-7 sequences that have insulator functions in late embryos and adults. We show that the LBC is enriched in nuclear extracts from late, but not early, embryos and that it contains three insulator proteins, GAF, Mod(mdg4), and E(y)2. Its DNA binding properties are unusual in that it requires a minimal sequence of >65 bp; however, other than a GAGA motif, the three Fab-7 LBC recognition elements display few sequence similarities. Finally, we show that mutations which abrogate LBC binding in vitro inactivate the Fab-7 boundary in the BX-C.

Functional role of dimerization and CP190 interacting domains of CTCF protein in Drosophila melanogaster

Background

Insulators play a central role in gene regulation, chromosomal architecture and genome function in higher eukaryotes. To learn more about how insulators carry out their diverse functions, we have begun an analysis of the Drosophila CTCF (dCTCF). CTCF is one of the few insulator proteins known to be conserved from flies to man.

Results

In the studies reported here we have focused on the identification and characterization of two dCTCF protein interaction modules. The first mediates dCTCF multimerization, while the second mediates dCTCF–CP190 interactions. The multimerization domain maps in the N-terminus of the dCTCF protein and likely mediates the formation of tetrameric complexes. The CP190 interaction module encompasses a sequence ~200 amino acids long that spans the C-terminal and mediates interactions with the N-terminal BTB domain of the CP190 protein. Transgene rescue experiments showed that a dCTCF protein lacking sequences critical for CP190 interactions was almost as effective as wild type in rescuing the phenotypic effects of a dCTCF null allele. The mutation did, however, affect CP190 recruitment to specific Drosophila insulator elements and had a modest effect on dCTCF chromatin association. A protein lacking the N-terminal dCTCF multimerization domain incompletely rescued the zygotic and maternal effect lethality of the null and did not rescue the defects in Abd-B regulation evident in surviving adult dCTCF mutant flies. Finally, we show that elimination of maternally contributed dCTCF at the onset of embryogenesis has quite different effects on development and Abd-B regulation than is observed when the homozygous mutant animals develop in the presence of maternally derived dCTCF activity.

Conclusions

Our results indicate that dCTCF–CP190 interactions are less critical for the in vivo functions of the dCTCF protein than the N-terminal dCTCF–dCTCF interaction domain. We also show that the phenotypic consequences of dCTCF mutations differ depending upon when and how dCTCF activity is lost.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0168-7) contains supplementary material, which is available to authorized users.

Pleiotropy of the Drosophila JAK pathway cytokine Unpaired 3 in development and aging

The Janus kinase (JAK) pathway is an essential, highly re-utilized developmental signaling cascade found in most metazoans. In vertebrates, the JAK intracellular cascade mediates signaling by dozens of cytokines and growth factors. In Drosophila, the Unpaired (Upd) family, encoded by three tandemly duplicated genes, is the only class of ligands associated with JAK stimulation. Unpaired has a central role in activation of JAK for most pathway functions, while Unpaired 2 regulates body size through insulin signaling. We show here that the third member of the family, unpaired 3 (upd3), overlaps upd in expression in some tissues and is essential for a subset of JAK-mediated developmental functions. First, consistent with the known requirements of JAK signaling in gametogenesis, we find that mutants of upd3 show an age-dependent impairment of fertility in both sexes. In oogenesis, graded JAK activity stimulated by Upd specifies the fates of the somatic follicle cells. As upd3 mutant females age, defects arise that can be attributed to perturbations of the terminal follicle cells, which require the highest levels of JAK activation. Therefore, in oogenesis, the activities of Upd and Upd3 both appear to quantitatively contribute to specification of those follicle cell fates. Furthermore, the sensitization of upd3 mutants to age-related decline in fertility can be used to investigate reproductive senescence. Second, loss of Upd3 during imaginal development results in defects of adult structures, including reduced eye size and abnormal wing and haltere posture. The outstretched wing and small eye phenotypes resemble classical alleles referred to as outstretched (os) mutations that have been previously ascribed to upd. However, we show that os alleles affect expression of both upd and upd3 and map to untranscribed regions, suggesting that they disrupt regulatory elements shared by both genes. Thus the upd region serves as a genetically tractable model for coordinate regulation of tandemly duplicated gene families that are commonly found in higher eukaryotes.

Deciphering the combinatorial architecture of a Drosophila homeotic gene enhancer

In Drosophila, the 330 kb bithorax complex regulates cellular differentiation along the anterior–posterior axis during development in the thorax and abdomen and is comprised of three homeotic genes: Ultrabithorax, abdominal-A, and Abdominal-B. The expression of each of these genes is in turn controlled through interactions between transcription factors and a number of cis-regulatory modules in the neighboring intergenic regions. In this study, we examine how the sequence architecture of transcription factor binding sites mediates the functional activity of one of these cis-regulatory modules. Using computational, mathematical modeling and experimental molecular genetic approaches we investigate the IAB7b enhancer, which regulates Abdominal-B expression specifically in the presumptive seventh and ninth abdominal segments of the early embryo. A cross-species comparison of the IAB7b enhancer reveals an evolutionarily conserved signature motif containing two FUSHI-TARAZU activator transcription factor binding sites. We find that the transcriptional repressors KNIRPS, KRUPPEL and GIANT are able to restrict reporter gene expression to the posterior abdominal segments, using different molecular mechanisms including short-range repression and competitive binding. Additionally, we show the functional importance of the spacing between the two FUSHI-TARAZU binding sites and discuss the potential importance of cooperativity for transcriptional activation. Our results demonstrate that the transcriptional output of the IAB7b cis-regulatory module relies on a complex set of combinatorial inputs mediated by specific transcription factor binding and that the sequence architecture at this enhancer is critical to maintain robust regulatory function.

Highly conserved ENY2/Sus1 protein binds to Drosophila CTCF and is required for barrier activity

Chromatin insulators affect interactions between promoters and enhancers/silencers and function as barriers for the spreading of repressive chromatin. Drosophila insulator protein dCTCF marks active promoters and boundaries of many histone H3K27 trimethylation domains associated with repressed chromatin. In particular, dCTCF binds to such boundaries between the parasegment-specific regulatory domains of the Bithorax complex. Here we demonstrate that the evolutionarily conserved protein ENY2 is recruited to the zinc-finger domain of dCTCF and is required for the barrier activity of dCTCF-dependent insulators in transgenic lines. Inactivation of ENY2 by RNAi in BG3 cells leads to the spreading of H3K27 trimethylation and Pc protein at several dCTCF boundaries. The results suggest that evolutionarily conserved ENY2 is responsible for barrier activity mediated by the dCTCF protein.

Hox gene regulation in the central nervous system of Drosophila

Hox genes specify the structures that form along the anteroposterior (AP) axis of bilateria. Within the genome, they often form clusters where, remarkably enough, their position within the clusters reflects the relative positions of the structures they specify along the AP axis. This correspondence between genomic organization and gene expression pattern has been conserved through evolution and provides a unique opportunity to study how chromosomal context affects gene regulation. In Drosophila, a general rule, often called “posterior dominance,” states that Hox genes specifying more posterior structures repress the expression of more anterior Hox genes. This rule explains the apparent spatial complementarity of Hox gene expression patterns in Drosophila. Here we review a noticeable exception to this rule where the more-posteriorly expressed Abd-B Hox gene fails to repress the more-anterior abd-A gene in cells of the central nervous system (CNS). While Abd-B is required to repress ectopic expression of abd-A in the posterior epidermis, abd-A repression in the posterior CNS is accomplished by a different mechanism that involves a large 92 kb long non-coding RNA (lncRNA) encoded by the intergenic region separating abd-A and Abd-B (the iab8ncRNA). Dissection of this lncRNA revealed that abd-A is repressed by the lncRNA using two redundant mechanisms. The first mechanism is mediated by a microRNA (mir-iab-8) encoded by intronic sequence within the large iab8-ncRNA. Meanwhile, the second mechanism seems to involve transcriptional interference by the long iab-8 ncRNA on the abd-A promoter. Recent work demonstrating CNS-specific regulation of genes by ncRNAs in Drosophila, seem to highlight a potential role for the iab-8-ncRNA in the evolution of the Drosophila Hox complexes.

Flanking sequence context-dependent transcription factor binding in early Drosophila development

BACKGROUND

Gene expression in the Drosophila embryo is controlled by functional interactions between a large network of protein transcription factors (TFs) and specific sequences in DNA cis-regulatory modules (CRMs). The binding site sequences for any TF can be experimentally determined and represented in a position weight matrix (PWM). PWMs can then be used to predict the location of TF binding sites in other regions of the genome, although there are limitations to this approach as currently implemented.

RESULTS

In this proof-of-principle study, we analyze 127 CRMs and focus on four TFs that control transcription of target genes along the anterio-posterior axis of the embryo early in development. For all four of these TFs, there is some degree of conserved flanking sequence that extends beyond the predicted binding regions. A potential role for these conserved flanking sequences may be to enhance the specificity of TF binding, as the abundance of these sequences is greatly diminished when we examine only predicted high-affinity binding sites.

CONCLUSIONS

Expanding PWMs to include sequence context-dependence will increase the information content in PWMs and facilitate a more efficient functional identification and dissection of CRMs.

A novel function for the Hox gene Abd-B in the male accessory gland regulates the long-term female post-mating response in Drosophila

In insects, products of the male reproductive tract are essential for initiating and maintaining the female post-mating response (PMR). The PMR includes changes in egg laying, receptivity to courting males, and sperm storage. In Drosophila, previous studies have determined that the main cells of the male accessory gland produce some of the products required for these processes. However, nothing was known about the contribution of the gland's other secretory cell type, the secondary cells. In the course of investigating the late functions of the homeotic gene, Abdominal-B (Abd-B), we discovered that Abd-B is specifically expressed in the secondary cells of the Drosophila male accessory gland. Using an Abd-B BAC reporter coupled with a collection of genetic deletions, we discovered an enhancer from the iab-6 regulatory domain that is responsible for Abd-B expression in these cells and that apparently works independently from the segmentally regulated chromatin domains of the bithorax complex. Removal of this enhancer results in visible morphological defects in the secondary cells. We determined that mates of iab-6 mutant males show defects in long-term egg laying and suppression of receptivity, and that products of the secondary cells are influential during sperm competition. Many of these phenotypes seem to be caused by a defect in the storage and gradual release of sex peptide in female mates of iab-6 mutant males. We also found that Abd-B expression in the secondary cells contributes to glycosylation of at least three accessory gland proteins: ovulin (Acp26Aa), CG1656, and CG1652. Our results demonstrate that long-term post-mating changes observed in mated females are not solely induced by main cell secretions, as previously believed, but that secondary cells also play an important role in male fertility by extending the female PMR. Overall, these discoveries provide new insights into how these two cell types cooperate to produce and maintain a robust female PMR.

Adaptive evolution and the birth of CTCF binding sites in the Drosophila genome

Comparative ChIP-seq data reveal adaptive evolution of insulator protein CTCF binding in multiple Drosophila species.

Changes in the physical interaction between cis-regulatory DNA sequences and proteins drive the evolution of gene expression. However, it has proven difficult to accurately quantify evolutionary rates of such binding change or to estimate the relative effects of selection and drift in shaping the binding evolution. Here we examine the genome-wide binding of CTCF in four species of Drosophila separated by between ∼2.5 and 25 million years. CTCF is a highly conserved protein known to be associated with insulator sequences in the genomes of human and Drosophila. Although the binding preference for CTCF is highly conserved, we find that CTCF binding itself is highly evolutionarily dynamic and has adaptively evolved. Between species, binding divergence increased linearly with evolutionary distance, and CTCF binding profiles are diverging rapidly at the rate of 2.22% per million years (Myr). At least 89 new CTCF binding sites have originated in the Drosophila melanogaster genome since the most recent common ancestor with Drosophila simulans. Comparing these data to genome sequence data from 37 different strains of Drosophila melanogaster, we detected signatures of selection in both newly gained and evolutionarily conserved binding sites. Newly evolved CTCF binding sites show a significantly stronger signature for positive selection than older sites. Comparative gene expression profiling revealed that expression divergence of genes adjacent to CTCF binding site is significantly associated with the gain and loss of CTCF binding. Further, the birth of new genes is associated with the birth of new CTCF binding sites. Our data indicate that binding of Drosophila CTCF protein has evolved under natural selection, and CTCF binding evolution has shaped both the evolution of gene expression and genome evolution during the birth of new genes.

A large proportion of the diversity of living organisms results from differential regulation of gene transcription. Transcriptional regulation is thought to differ between species because of evolutionary changes in the physical interactions between regulatory DNA elements and DNA-binding proteins; these can generate variation in the spatial and temporal patterns of gene expression. The mechanisms by which these protein–DNA interactions evolve is therefore an important question in evolutionary biology. Does adaptive evolution play a role, or is the process dominated by neutral genetic drift? Insulator proteins are a special group of DNA-binding proteins—instead of directly serving to activate or repress genes, they can function to coordinate the interactions between other regulatory elements (such as enhancers and promoters). Additionally, insulator proteins can limit the spreading of chromatin condensation and help to demarcate the boundaries of regulatory domains in the genome. In spite of their critical role in genome regulation, little is known about the evolution of interactions between insulator proteins and DNA. Here, we use ChIP-seq to examine the distribution of binding sites for CTCF, a highly conserved insulator protein, in four closely related Drosophila species. We find that genome-wide binding profiles of CTCF are highly dynamic across evolutionary time, with frequent births of new CTCF-DNA interactions, and we demonstrate that this evolutionary process is driven by natural selection. By comparing these with RNA-seq data, we find that gain or loss of CTCF binding impacts the expression levels of nearby genes and correlates with structural evolution of the genome. Together these results suggest a potential mechanism of regulatory re-wiring through adaptive evolution of CTCF binding.

abd-A regulation by the iab-8 noncoding RNA

The homeotic genes in Drosophila melanogaster are aligned on the chromosome in the order of the body segments that they affect. The genes affecting the more posterior segments repress the more anterior genes. This posterior dominance rule must be qualified in the case of abdominal-A (abd-A) repression by Abdominal-B (Abd-B). Animals lacking Abd-B show ectopic expression of abd-A in the epidermis of the eighth abdominal segment, but not in the central nervous system. Repression in these neuronal cells is accomplished by a 92 kb noncoding RNA. This “iab-8 RNA” produces a micro RNA to repress abd-A, but also has a second, redundant repression mechanism that acts only “in cis.” Transcriptional interference with the abd-A promoter is the most likely mechanism.

Although long, noncoding RNAs have been found in many organisms, it has been difficult to assign to them any molecular function. The homeotic gene clusters in the fruit fly, Drosophila melanogaster, contain many such noncoding RNAs. We have characterized one such noncoding RNA, a 92 kb transcription unit from within the bithorax complex. This transcript, called the iab-8 ncRNA, is made in the cells of the central nervous system in the eighth abdominal segment, along with the homeotic transcription factor Abdominal-B. Another homeotic transcription factor, abdominal-A, is repressed in these cells. It has generally been assumed that abdominal-A repression in these cells is mediated by the Abdominal-B protein. However, here we show that it is not Abdominal-B that represses abdominal-A, but the iab-8 ncRNA. This repression is accomplished by two redundant mechanisms; the iab-8 precursor produces a micro RNA, which targets the abdominal-A mRNA, and iab-8 transcription interferes with the abdominal-A promoter, which lies just downstream of the iab-8 ncRNA poly(A) site.

Developmental transcriptional networks are required to maintain neuronal subtype identity in the mature nervous system

During neurogenesis, transcription factors combinatorially specify neuronal fates and then differentiate subtype identities by inducing subtype-specific gene expression profiles. But how is neuronal subtype identity maintained in mature neurons? Modeling this question in two Drosophila neuronal subtypes (Tv1 and Tv4), we test whether the subtype transcription factor networks that direct differentiation during development are required persistently for long-term maintenance of subtype identity. By conditional transcription factor knockdown in adult Tv neurons after normal development, we find that most transcription factors within the Tv1/Tv4 subtype transcription networks are indeed required to maintain Tv1/Tv4 subtype-specific gene expression in adults. Thus, gene expression profiles are not simply “locked-in,” but must be actively maintained by persistent developmental transcription factor networks. We also examined the cross-regulatory relationships between all transcription factors that persisted in adult Tv1/Tv4 neurons. We show that certain critical cross-regulatory relationships that had existed between these transcription factors during development were no longer present in the mature adult neuron. This points to key differences between developmental and maintenance transcriptional regulatory networks in individual neurons. Together, our results provide novel insight showing that the maintenance of subtype identity is an active process underpinned by persistently active, combinatorially-acting, developmental transcription factors. These findings have implications for understanding the maintenance of all long-lived cell types and the functional degeneration of neurons in the aging brain.

For neurons to function properly, they must establish and then maintain their unique, subtype-specific gene expression profiles. These unique gene expression profiles are established during development by networks of DNA–binding proteins, termed transcription factors (TFs). However, how neurons maintain their unique gene expression profiles in the mature and aging brain is largely unknown. Recent advances in inducible genetic technologies now allow us to manipulate gene expression in adult neurons, after normal development. Applying such techniques, we examined the effect of knocking down TF expression in two adult neuronal subtypes. We show that the TF networks that establish unique gene expression profiles during development are then required to maintain them thereafter. Thus, gene expression profiles are not simply “locked-in,” but must be actively maintained by persistent developmental TF networks. However, we found that critical cross-regulatory relationships that had existed between TFs during development were not present in the adult, even between persisting TFs. This highlights important differences between developmental and maintenance transcriptional networks in individual neurons. The dependence of subtype gene expression on active mechanisms represents a potential Achilles heel for long-lived cells, as deterioration of those active mechanisms could lead to functional degeneration of neurons with advancing age.

Molecular dissection of cis-regulatory modules at the Drosophila bithorax complex reveals critical transcription factor signature motifs

At the Drosophila melanogaster bithorax complex (BX-C) over 330kb of intergenic DNA is responsible for directing the transcription of just three homeotic (Hox) genes during embryonic development. A number of distinct enhancer cis-regulatory modules (CRMs) are responsible for controlling the specific expression patterns of the Hox genes in the BX-C. While it has proven possible to identify orthologs of known BX-C CRMs in different Drosophila species using overall sequence conservation, this approach has not proven sufficiently effective for identifying novel CRMs or defining the key functional sequences within enhancer CRMs. Here we demonstrate that the specific spatial clustering of transcription factor (TF) binding sites is important for BX-C enhancer activity. A bioinformatic search for combinations of putative TF binding sites in the BX-C suggests that simple clustering of binding sites is frequently not indicative of enhancer activity. However, through molecular dissection and evolutionary comparison across the Drosophila genus we discovered that specific TF binding site clustering patterns are an important feature of three known BX-C enhancers. Sub-regions of the defined IAB5 and IAB7b enhancers were both found to contain an evolutionarily conserved signature motif of clustered TF binding sites which is critical for the functional activity of the enhancers. Together, these results indicate that the spatial organization of specific activator and repressor binding sites within BX-C enhancers is of greater importance than overall sequence conservation and is indicative of enhancer functional activity.

Genome-wide analysis of the binding of the Hox protein Ultrabithorax and the Hox cofactor Homothorax in Drosophila

Hox genes encode a family of transcription factors that are key developmental regulators with a highly conserved role in specifying segmental diversity along the metazoan body axis. Although they have been shown to regulate a wide variety of downstream processes, direct transcriptional targets have been difficult to identify and this has been a major obstacle to our understanding of Hox gene function. We report the identification of genome-wide binding sites for the Hox protein Ultrabithorax (Ubx) using a YFP-tagged Drosophila protein-trap line together with chromatin immunoprecipitation and microarray analysis. We identify 1,147 genes bound by Ubx at high confidence in chromatin from the haltere imaginal disc, a prominent site of Ubx function where it specifies haltere versus wing development. The functional relevance of these genes is supported by their overlap with genes differentially expressed between wing and haltere imaginal discs. The Ubx-bound gene set is highly enriched in genes involved in developmental processes and contains both high-level regulators as well as genes involved in more basic cellular functions. Several signalling pathways are highly enriched in the Ubx target gene set and our analysis supports the view that Hox genes regulate many levels of developmental pathways and have targets distributed throughout the gene network. We also performed genome-wide analysis of the binding sites for the Hox cofactor Homothorax (Hth), revealing a striking similarity with the Ubx binding profile. We suggest that these binding profiles may be strongly influenced by chromatin accessibility and provide evidence of a link between Ubx/Hth binding and chromatin state at genes regulated by Polycomb silencing. Overall, we define a set of direct Ubx targets in the haltere imaginal disc and suggest that chromatin accessibility has important implications for Hox target selection and for transcription factor binding in general.

Gene expression in time and space: additive vs hierarchical organization of cis-regulatory regions

In higher eukaryotes, individual genes are often intermingled with other genes and spread out across tens to hundreds of kilobases, even though only small portions of their sequence are devoted to protein coding. Yet, in this seemingly extended and tangled mess, the cell is able to precisely regulate gene expression in both time and space. Over the past few decades, numerous elements, like enhancers, silencers and insulators have been found that shed some light on how the precise control of gene expression is achieved. Through these discoveries, an additive model of gene expression was envisioned, where the addition of the patterning details imparted by regulatory elements would create the final pattern of gene expression. Although many genes can be described using this model, recent work in the Drosophila bithorax complex suggests that this model may be somewhat simplistic and, in fact, regulatory elements sometimes seem to communicate with each other to form a functional hierarchy that is far from additive.

Disruption of the abdominal-B promoter tethering element results in a loss of long-range enhancer-directed Hox gene expression in Drosophila

There are many examples within gene complexes of transcriptional enhancers interacting with only a subset of target promoters. A number of molecular mechanisms including promoter competition, insulators and chromatin looping are thought to play a role in regulating these interactions. At the Drosophila bithorax complex (BX-C), the IAB5 enhancer specifically drives gene expression only from the Abdominal-B (Abd-B) promoter, even though the enhancer and promoter are 55 kb apart and are separated by at least three insulators. In previous studies, we discovered that a 255 bp cis-regulatory module, the promoter tethering element (PTE), located 5′ of the Abd-B transcriptional start site is able to tether IAB5 to the Abd-B promoter in transgenic embryo assays. In this study we examine the functional role of the PTE at the endogenous BX-C using transposon-mediated mutagenesis. Disruption of the PTE by P element insertion results in a loss of enhancer-directed Abd-B expression during embryonic development and a homeotic transformation of abdominal segments. A partial deletion of the PTE and neighboring upstream genomic sequences by imprecise excision of the P element also results in a similar loss of Abd-B expression in embryos. These results demonstrate that the PTE is an essential component of the regulatory network at the BX-C and is required in vivo to mediate specific long-range enhancer-promoter interactions.

Initiator elements function to determine the activity state of BX-C enhancers

A >300 kb cis-regulatory region is required for the proper expression of the three bithorax complex (BX-C) homeotic genes. Based on genetic and transgenic analysis, a model has been proposed in which the numerous BX-C cis-regulatory elements are spatially restricted through the activation or repression of parasegment-specific chromatin domains. Particular early embryonic enhancers, called initiators, have been proposed to control this complex process. Here, in order to better understand the process of domain activation, we have undertaken a systematic in situ dissection of the iab-6 cis-regulatory domain using a new method, called InSIRT. Using this method, we create and genetically characterize mutations affecting iab-6 function, including mutations specifically modifying the iab-6 initiator. Through our mutagenesis of the iab-6 initiator, we provide strong evidence that initiators function not to directly control homeotic gene expression but rather as domain control centers to determine the activity state of the enhancers and silencers within a cis-regulatory domain.

Understanding how genes become activated is one of the primary areas of research in modern biology. In order to decipher the DNA components required for this process, scientists have traditionally turned to transgenic reporter assays, where DNA elements are removed from their native environment and placed next to a simplified reporter gene to monitor transcriptional activation. Although this approach is powerful, it can result in artifacts stemming from the channelization of regulatory element activities into predetermined classes. In this manuscript, we investigate the biological role of elements from the Drosophila bithorax complex, called initiators. In transgenic assays, these elements have been categorized as enhancers. However, genetic analysis suggests that, in situ, these elements perform a far more complex function. Here, using a new method to repeatedly target a genetic locus for mutagenesis, we show that initiators function as control elements that coordinate the activity of nearby enhancers and silencers. Overall, our study highlights how gene expression can be controlled through a hierarchical arrangement of cis-regulatory elements.

Selective interactions of boundaries with upstream region of Abd-B promoter in Drosophila bithorax complex and role of dCTCF in this process

Expression of the genes Ubx, abd-A, and Abd-B of the bithorax complex depends on its cis-regulatory region, which is divided into discrete functional domains (iab). Boundary/insulator elements, named Mcp, Fab-6, Fab-7 and Fab-8 (PTS/F8), have been identified at the borders of the iab domains. Recently, binding sites for a Drosophila homolog of the vertebrate insulator protein CTCF have been identified in Mcp, Fab-6 and Fab-8 and also in several regions that correspond to predicted boundaries, Fab-3 and Fab-4 in particular. Taking into account the inability of the yeast GAL4 activator to stimulate the white promoter when the activator and the promoter are separated by a 5-kb yellow gene, we have tested functional interactions between the boundaries. The results show that all dCTCF-containing boundaries interact with each other. However, inactivation of dCTCF binding sites in Mcp, Fab-6 and PTS/F8 only partially reduces their ability to interact, suggesting the presence of additional protein(s) supporting distant interactions between the boundaries. Interestingly, only Fab-6, Fab-7 (which contains no dCTCF binding sites) and PTS/F8 interact with the upstream region of the Abd-B promoter. Thus, the boundaries might be involved in supporting the specific interactions between iab enhancers and promoters of the bithorax complex.

Chromosomal organization at the level of gene complexes

Metazoan genomes primarily consist of non-coding DNA in comparison to coding regions. Non-coding fraction of the genome contains cis-regulatory elements, which ensure that the genetic code is read properly at the right time and space during development. Regulatory elements and their target genes define functional landscapes within the genome, and some developmentally important genes evolve by keeping the genes involved in specification of common organs/tissues in clusters and are termed gene complex. The clustering of genes involved in a common function may help in robust spatio-temporal gene expression. Gene complexes are often found to be evolutionarily conserved, and the classic example is the hox complex. The evolutionary constraints seen among gene complexes provide an ideal model system to understand cis and trans-regulation of gene function. This review will discuss the various characteristics of gene regulatory modules found within gene complexes and how they can be characterized.