Dynamic changes in RNA-protein interactions and RNA secondary structure in mammalian erythropoiesis.

Two features of eukaryotic RNA molecules that regulate their post-transcriptional fates are RNA secondary structure and RNA-binding protein (RBP) interaction sites. However, a comprehensive global overview of the dynamic nature of these sequence features during erythropoiesis has never been obtained. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approach to reveal the global landscape of RNA secondary structure and RBP-RNA interaction sites and the dynamics of these features during this important developmental process. We identify dynamic patterns of RNA secondary structure and RBP binding throughout the process and determine a set of corresponding protein-bound sequence motifs along with their dynamic structural and RBP-binding contexts. Finally, using these dynamically bound sequences, we identify a number of RBPs that have known and putative key functions in post-transcriptional regulation during mammalian erythropoiesis. In total, this global analysis reveals new post-transcriptional regulators of mammalian blood cell development.

RNA-Binding Proteins PCBP1 and PCBP2 Are Critical Determinants of Murine Erythropoiesis.

We previously demonstrated that the two paralogous RNA-binding proteins PCBP1 and PCBP2 are individually essential for mouse development: Pcbp1-null embryos are peri-implantation lethal, while Pcbp2-null embryos lose viability at midgestation. Midgestation Pcbp2-/- embryos revealed a complex phenotype that included loss of certain hematopoietic determinants. Whether PCBP2 directly contributes to erythropoietic differentiation and whether PCBP1 has a role in this process remained undetermined. Here, we selectively inactivated the genes encoding these two RNA-binding proteins during differentiation of the erythroid lineage in the developing mouse embryo. Individual inactivation of either locus failed to impact viability or blood formation. However, combined inactivation of the two loci resulted in midgestational repression of erythroid/hematopoietic gene expression, loss of blood formation, and fetal demise. Orthogonal ex vivo analyses of primary erythroid progenitors selectively depleted of these two RNA-binding proteins revealed that they mediate a combination of overlapping and isoform-specific impacts on hematopoietic lineage transcriptome, impacting both mRNA representation and exon splicing. These data lead us to conclude that PCBP1 and PCBP2 mediate functions critical to differentiation of the erythroid lineage.

Multi-omic profiling of pituitary thyrotropic cells and progenitors.

BACKGROUND: The pituitary gland is a neuroendocrine organ containing diverse cell types specialized in secreting hormones that regulate physiology. Pituitary thyrotropes produce thyroid-stimulating hormone (TSH), a critical factor for growth and maintenance of metabolism. The transcription factors POU1F1 and GATA2 have been implicated in thyrotrope fate, but the transcriptomic and epigenomic landscapes of these neuroendocrine cells have not been characterized. The goal of this work was to discover transcriptional regulatory elements that drive thyrotrope fate. RESULTS: We identified the transcription factors and epigenomic changes in chromatin that are associated with differentiation of POU1F1-expressing progenitors into thyrotropes using cell lines that represent an undifferentiated Pou1f1 lineage progenitor (GHF-T1) and a committed thyrotrope line that produces TSH (TαT1). We compared RNA-seq, ATAC-seq, histone modification (H3K27Ac, H3K4Me1, and H3K27Me3), and POU1F1 binding in these cell lines. POU1F1 binding sites are commonly associated with bZIP transcription factor consensus binding sites in GHF-T1 cells and Helix-Turn-Helix (HTH) or basic Helix-Loop-Helix (bHLH) factors in TαT1 cells, suggesting that these classes of transcription factors may recruit or cooperate with POU1F1 binding at unique sites. We validated enhancer function of novel elements we mapped near Cga, Pitx1, Gata2, and Tshb by transfection in TαT1 cells. Finally, we confirmed that an enhancer element near Tshb can drive expression in thyrotropes of transgenic mice, and we demonstrate that GATA2 enhances Tshb expression through this element. CONCLUSION: These results extend the ENCODE multi-omic profiling approach to the pituitary gland, which should be valuable for understanding pituitary development and disease pathogenesis.

Single-cell transcriptomic analysis of adult mouse pituitary reveals sexual dimorphism and physiologic demand-induced cellular plasticity.

The anterior pituitary gland drives highly conserved physiologic processes in mammalian species. These hormonally controlled processes are central to somatic growth, pubertal transformation, fertility, lactation, and metabolism. Current cellular models of mammalian anteiror pituitary, largely built on candidate gene based immuno-histochemical and mRNA analyses, suggest that each of the seven hormones synthesized by the pituitary is produced by a specific and exclusive cell lineage. However, emerging evidence suggests more complex relationship between hormone specificity and cell plasticity. Here we have applied massively parallel single-cell RNA sequencing (scRNA-seq), in conjunction with complementary imaging-based single-cell analyses of mRNAs and proteins, to systematically map both cell-type diversity and functional state heterogeneity in adult male and female mouse pituitaries at single-cell resolution and in the context of major physiologic demands. These quantitative single-cell analyses reveal sex-specific cell-type composition under normal pituitary homeostasis, identify an array of cells associated with complex complements of hormone-enrichment, and undercover non-hormone producing interstitial and supporting cell-types. Interestingly, we also identified a Pou1f1-expressing cell population that is characterized by a unique multi-hormone gene expression profile. In response to two well-defined physiologic stresses, dynamic shifts in cellular diversity and transcriptome profiles were observed for major hormone producing and the putative multi-hormone cells. These studies reveal unanticipated cellular complexity and plasticity in adult pituitary, and provide a rich resource for further validating and expanding our molecular understanding of pituitary gene expression programs and hormone production.

The Transcription Factor NR4A2 Plays an Essential Role in Driving Prolactin Expression in Female Pituitary Lactotropes.

Differentiation of the hormone-producing cells of the pituitary represents an informative model of cell fate determination. The generation and maintenance of 2 pituitary lineages, the growth hormone (GH)- producing somatotropes and the prolactin (PRL)- producing lactotropes, are dependent on the pituitary-specific transcription factor, POU1F1. While POU1F1 is expressed in both cell types, and plays a role in activation of both the Gh and Prl genes, expression of Gh and Prl is restricted to somatotropes and lactotropes, respectively. These observations imply the existence of additional factors that contribute to the somatotrope and lactotrope identities and their hormone expressions. Prior transcriptome analysis of primary somatotropes and lactotropes isolated from the mouse pituitary identified enrichment of a transcription factor, Nr4a2, in the lactotropes. Nr4a2 was shown in a cell culture model to bind the Prl promoter at a position adjacent to Pou1f1 and to synergize with Pou1f1 in driving Prl transcription. Here we demonstrate in vivo the role of Nr4a2 as an enhancer of Prl expression by conditional gene inactivation of the Nr4a2 gene in mouse lactotropes. We demonstrate that nuclear orphan receptor transcription factor (NR4A2) binding at the Prl promoter is dependent on actions of POU1F1; while POU1F1 is essential to loading polymerase (Pol) II on the Prl promoter, Nr4a2 plays a role in enhancing Pol II release into the Prl gene body. These studies establish an in vivo role of Nr4a2 in enhancing Prl expression in mouse lactotropes, explore its mechanism of action, and establish a system for further study of the lactotrope lineage in the pituitary.

A cytosine-rich splice regulatory determinant enforces functional processing of the human α-globin gene transcript.

The establishment of efficient and stable splicing patterns in terminally differentiated cells is critical to maintenance of specific functions throughout the lifespan of an organism. The human α-globin (hα-globin) gene contains 3 exons separated by 2 short introns. Naturally occurring α-thalassemia mutations that trigger aberrant splicing have revealed the presence of cryptic splice sites within the hα-globin gene transcript. How cognate (functional) splice sites are selectively used in lieu of these cryptic sites has remained unexplored. Here we demonstrate that the preferential selection of a cognate splice donor essential to functional splicing of the hα-globin transcript is dependent on the actions of an intronic cytosine (C)-rich splice regulatory determinant and its interacting polyC-binding proteins. Inactivation of this determinant by mutation of the C-rich element or by depletion of polyC-binding proteins triggers a dramatic shift in splice donor activity to an upstream, out-of-frame, cryptic donor. The essential role of the C-rich element in hα-globin gene expression is supported by its coevolution with the cryptic donor site in primate species. These data lead us to conclude that an intronic C-rich determinant enforces functional splicing of the hα-globin transcript, thus acting as an obligate determinant of hα-globin gene expression.

Transcriptome Analyses of Female Somatotropes and Lactotropes Reveal Novel Regulators of Cell Identity in the Pituitary.

The differentiation of the hormone-producing cell lineages of the anterior pituitary represents an informative model of mammalian cell fate determination. The generation and maintenance of two of these lineages, the GH-producing somatotropes and prolactin (PRL)-producing lactotropes, are dependent on the pituitary-specific transcription factor POU1F1. Whereas POU1F1 is expressed in both cell types, and plays a direct role in the activation of both the Gh and Prl genes, GH expression is restricted to somatotropes and PRL expression is restricted to lactotropes. These observations imply the existence of additional, cell type-enriched factors that contribute to the somatotrope and lactotrope cell identities. In this study, we use transgenic mouse models to facilitate sorting of somatotrope and lactotrope populations based on the expression of fluorescent markers expressed under Gh and Prl gene transcriptional controls. The transcriptomic analyses reveal a concordance of gene expression profiles in the two populations. The limited number of divergent mRNAs between the two populations includes a set of transcription factors that may have roles in pituitary lineage divergence and/or in regulating expression of cell type-specific genes after differentiation. Four of these factors were validated for lineage enrichment at the level of protein expression, two somatotrope enriched and two lactotrope enriched. Three of these four factors were shown to have corresponding activities in appropriate enhancement or repression of landmark genes in a cell culture model system. These studies identify novel regulators of the somatotropes and lactotropes, and they establish a useful database for further study of these lineages in the anterior pituitary.

Domain-focused CRISPR screen identifies HRI as a fetal hemoglobin regulator in human erythroid cells.

Increasing fetal hemoglobin (HbF) levels in adult red blood cells provides clinical benefit to patients with sickle cell disease and some forms of ß-thalassemia. To identify potentially druggable HbF regulators in adult human erythroid cells, we employed a protein kinase domain-focused CRISPR-Cas9-based genetic screen with a newly optimized single-guide RNA scaffold. The screen uncovered the heme-regulated inhibitor HRI (also known as EIF2AK1), an erythroid-specific kinase that controls protein translation, as an HbF repressor. HRI depletion markedly increased HbF production in a specific manner and reduced sickling in cultured erythroid cells. Diminished expression of the HbF repressor BCL11A accounted in large part for the effects of HRI depletion. Taken together, these results suggest HRI as a potential therapeutic target for hemoglobinopathies.

Poly(C)-Binding Protein Pcbp2 Enables Differentiation of Definitive Erythropoiesis by Directing Functional Splicing of the Runx1 Transcript.

Formation of the mammalian hematopoietic system is under a complex set of developmental controls. Here, we report that mouse embryos lacking the KH domain poly(C) binding protein, Pcbp2, are selectively deficient in the definitive erythroid lineage. Compared to wild-type controls, transcript splicing analysis of the Pcbp2-/- embryonic liver reveals accentuated exclusion of an exon (exon 6) that encodes a highly conserved transcriptional control segment of the hematopoietic master regulator, Runx1. Embryos rendered homozygous for a Runx1 locus lacking this cassette exon (Runx1ΔE6) effectively phenocopy the loss of the definitive erythroid lineage in Pcbp2-/- embryos. These data support a model in which enhancement of Runx1 cassette exon 6 inclusion by Pcbp2 serves a critical role in development of hematopoietic progenitors and constitutes a critical step in the developmental pathway of the definitive erythropoietic lineage.

PolyC-binding proteins enhance expression of the CDK2 cell cycle regulatory protein via alternative splicing.

The PolyC binding proteins (PCBPs) impact alternative splicing of a subset of mammalian genes that are enriched in basic cellular functions. Here, we focus our analysis on PCBP-controlled cassette exon-splicing within the cell cycle control regulator cyclin-dependent kinase-2 (CDK2) transcript. We demonstrate that PCBP binding to a C-rich polypyrimidine tract (PPT) preceding exon 5 of the CDK2 transcript enhances cassette exon inclusion. This splice enhancement is U2AF65-independent and predominantly reflects actions of the PCBP1 isoform. Remarkably, PCBPs' control of CDK2 ex5 splicing has evolved subsequent to mammalian divergence via conversion of constitutive exon 5 inclusion in the mouse CDK2 transcript to PCBP-responsive exon 5 alternative splicing in humans. Importantly, exclusion of exon 5 from the hCDK2 transcript dramatically represses the expression of CDK2 protein with a corresponding perturbation in cell cycle kinetics. These data highlight a recently evolved post-transcriptional pathway in primate species with the potential to modulate cell cycle control.

Long-range looping of a locus control region drives tissue-specific chromatin packing within a multigene cluster.

The relationships of higher order chromatin organization to mammalian gene expression remain incompletely defined. The human Growth Hormone (hGH) multigene cluster contains five gene paralogs. These genes are selectively activated in either the pituitary or the placenta by distinct components of a remote locus control region (LCR). Prior studies have revealed that appropriate activation of the placental genes is dependent not only on the actions of the LCR, but also on the multigene composition of the cluster itself. Here, we demonstrate that the hGH LCR 'loops' over a distance of 28 kb in primary placental nuclei to make specific contacts with the promoters of the two GH genes in the cluster. This long-range interaction sequesters the GH genes from the three hCS genes which co-assemble into a tightly packed 'hCS chromatin hub'. Elimination of the long-range looping, via specific deletion of the placental LCR components, triggers a dramatic disruption of the hCS chromatin hub. These data reveal a higher-order structural pathway by which long-range looping from an LCR impacts on local chromatin architecture that is linked to tissue-specific gene regulation within a multigene cluster.

αCP binding to a cytosine-rich subset of polypyrimidine tracts drives a novel pathway of cassette exon splicing in the mammalian transcriptome.

Alternative splicing (AS) is a robust generator of mammalian transcriptome complexity. Splice site specification is controlled by interactions of cis-acting determinants on a transcript with specific RNA binding proteins. These interactions are frequently localized to the intronic U-rich polypyrimidine tracts (PPT) located 5' to the majority of splice acceptor junctions. αCPs (also referred to as polyC-binding proteins (PCBPs) and hnRNPEs) comprise a subset of KH-domain proteins with high affinity and specificity for C-rich polypyrimidine motifs. Here, we demonstrate that αCPs promote the splicing of a defined subset of cassette exons via binding to a C-rich subset of polypyrimidine tracts located 5' to the αCP-enhanced exonic segments. This enhancement of splice acceptor activity is linked to interactions of αCPs with the U2 snRNP complex and may be mediated by cooperative interactions with the canonical polypyrimidine tract binding protein, U2AF65. Analysis of αCP-targeted exons predicts a substantial impact on fundamental cell functions. These findings lead us to conclude that the αCPs play a direct and global role in modulating the splicing activity and inclusion of an array of cassette exons, thus driving a novel pathway of splice site regulation within the mammalian transcriptome.

Cytoplasmic poly(A) binding protein-1 binds to genomically encoded sequences within mammalian mRNAs.

The functions of the major mammalian cytoplasmic poly(A) binding protein, PABPC1, have been characterized predominantly in the context of its binding to the 3' poly(A) tails of mRNAs. These interactions play important roles in post-transcriptional gene regulation by enhancing translation and mRNA stability. Here, we performed transcriptome-wide CLIP-seq analysis to identify additional PABPC1 binding sites within genomically encoded mRNA sequences that may impact on gene regulation. From this analysis, we found that PABPC1 binds directly to the canonical polyadenylation signal in thousands of mRNAs in the mouse transcriptome. PABPC1 binding also maps to translation initiation and termination sites bracketing open reading frames, exemplified most dramatically in replication-dependent histone mRNAs. Additionally, a more restricted subset of PABPC1 interaction sites comprised A-rich sequences within the 5' UTRs of mRNAs, including Pabpc1 mRNA itself. Functional analyses revealed that these PABPC1 interactions in the 5' UTR mediate both auto- and trans-regulatory translational control. In total, these findings reveal a repertoire of PABPC1 binding that is substantially broader than previously recognized with a corresponding potential to impact and coordinate post-transcriptional controls critical to a broad array of cellular functions.

Functional characterization of a human POU1F1 mutation associated with isolated growth hormone deficiency: a novel etiology for IGHD.

POU1F1, a pituitary-specific POU-homeo domain transcription factor, plays an essential role in the specification of the somatotroph, lactotroph and thyrotroph lineages and in the activation of GH1, PRL and TSHß transcription. Individuals with mutations in POU1F1 present with combined deficiency of GH, PRL and TSH. Here, we identified a heterozygous missense mutation with evidence of pathogenicity, at the POU1F1 locus, in a large family in which an isolated growth hormone deficiency segregates as an autosomal dominant trait. The corresponding p.Pro76Leu mutation maps to a conserved site within the POU1F1 transactivation domain. Bandshift assays revealed that the mutation alters wild-type POU1F1 binding to cognate sites within the hGH-LCR and hGH1 promoter, but not to sites within the PRL promoter, and it selectively increases binding affinity to sites within the hGH-LCR. Co-immunoprecipitation studies reveal that this substitution enhances interactions of POU1F1 with three of its cofactors, PITX1, LHX3a and ELK1, and that residue 76 plays a critical role in these interactions. The insertion of the mutation at the mouse Pou1f1 locus results in a dramatic loss of protein expression despite normal mRNA concentrations. Mice heterozygous for the p.Pro76Leu mutation were phenotypically normal while homozygotes demonstrated a dwarf phenotype. Overall, this study unveils the involvement of POU1F1 in dominantly inherited isolated GH deficiency and demonstrates a significant impact of the Pro76Leu mutation on DNA-binding activities, alterations in transactivating functions and interactions with cofactors. Our data further highlight difficulties in modeling human genetic disorders in the mouse despite apparent conservation of gene expression pathways and physiologic functions.

The Poly(C) Binding Protein Pcbp2 and Its Retrotransposed Derivative Pcbp1 Are Independently Essential to Mouse Development.

RNA-binding proteins participate in a complex array of posttranscriptional controls essential to cell type specification and somatic development. Despite their detailed biochemical characterizations, the degree to which each RNA-binding protein impacts mammalian embryonic development remains incompletely defined, and the level of functional redundancy among subsets of these proteins remains open to question. The poly(C) binding proteins, PCBPs (αCPs and hnRNP E proteins), are encoded by a highly conserved and broadly expressed gene family. The two major Pcbp isoforms, Pcbp2 and Pcbp1, are robustly expressed in a wide range of tissues and exert both nuclear and cytoplasmic controls over gene expression. Here, we report that Pcbp1-null embryos are rendered nonviable in the peri-implantation stage. In contrast, Pcbp2-null embryos undergo normal development until midgestation (12.5 to 13.5 days postcoitum), at which time they undergo a dramatic loss in viability associated with combined cardiovascular and hematopoietic abnormalities. Mice heterozygous for either Pcbp1 or Pcbp2 null alleles display a mild and nondisruptive defect in initial postpartum weight gain. These data reveal that Pcbp1 and Pcbp2 are individually essential for mouse embryonic development and have distinct impacts on embryonic viability and that Pcpb2 has a nonredundant in vivo role in hematopoiesis. These data further provide direct evidence that Pcbp1, a retrotransposed derivative of Pcpb2, has evolved an essential function(s) in the mammalian genome.

Resistance of mRNAs with AUG-proximal nonsense mutations to nonsense-mediated decay reflects variables of mRNA structure and translational activity.

Nonsense-mediated mRNA decay (NMD) is a surveillance pathway that recognizes and selectively degrades mRNAs carrying premature termination codons (PTCs). The level of sensitivity of a PTC-containing mRNA to NMD is multifactorial. We have previously shown that human ß-globin mRNAs carrying PTCs in close proximity to the translation initiation AUG codon escape NMD. This was called the 'AUG-proximity effect'. The present analysis of nonsense codons in the human α-globin mRNA illustrates that the determinants of the AUG-proximity effect are in fact quite complex, reflecting the ability of the ribosome to re-initiate translation 3' to the PTC and the specific sequence and secondary structure of the translated ORF. These data support a model in which the time taken to translate the short ORF, impacted by distance, sequence, and structure, not only modulates translation re-initiation, but also impacts on the exact boundary of AUG-proximity protection from NMD.

An autoregulatory pathway establishes the definitive chromatin conformation at the pit-1 locus.

The transcription factor Pit-1 (POU1-F1) plays a dominant role in cell lineage expansion and differentiation in the anterior pituitary. Prior studies of the mouse Pit-1 (mPit-1) gene revealed that this master regulatory locus is activated at embryonic day 13.5 (E13.5) by an early enhancer (EE), whereas its subsequent expression throughout adult life is maintained by a more distal definitive enhancer (DE). Here, we demonstrate that the sequential actions of these two enhancers are linked to corresponding shifts in their proximities to the Pit-1 promoter. We further demonstrate that the looping of the definitive enhancer to the mPit-1 promoter is critically dependent on a self-sustaining autoregulatory mechanism mediated by the Pit-1 protein. These Pit-1-dependent actions are accompanied by localized recruitment of CBP and enrichment for H3K27 acetylation within the Pit-1 locus. These data support a model in which the sequential actions of two developmentally activated enhancers are linked to a corresponding shift in higher-order chromatin structures. This shift establishes an autoregulatory circuit that maintains durable expression of Pit-1 throughout adult life.

Autonomous actions of the human growth hormone long-range enhancer.

The human growth hormone (hGH) gene is controlled by a long-range enhancer, HSI, located 14.5 kb 5' to the hGH promoter. HSI establishes a domain of noncoding transcription that is 'looped' to the hGH promoter as an essential step in initiating hGH gene expression. Thus, defining how HSI generates its domain of noncoding transcription is central to understanding its long-range function. Here, we demonstrate that activation of noncoding transcription reflects an HSI-autonomous activity fully independent of interactions with linked gene promoters and occurring in spatial and temporal synchrony with initiation of GH expression in the embryonic pituitary. HSI establishes its noncoding transcription start sites (TSS) over a defined distance from its core determinants and in a manner independent of local primary sequences. The interval between HSI and it TSS co-maps with a domain of disordered and/or highly mobile nucleosomes specific to the pituitary locus. Thus, a localized chromatin reconfiguration by HSI and consequent establishment of an adjacent domain of noncoding transcription constitute initiating events in long-range enhancer function within the hGH locus.

Long-range enhancer activity determines Myc sensitivity to Notch inhibitors in T cell leukemia.

Notch is needed for T-cell development and is a common oncogenic driver in T-cell acute lymphoblastic leukemia. The protooncogene c-Myc (Myc) is a critical target of Notch in normal and malignant pre-T cells, but how Notch regulates Myc is unknown. Here, we identify a distal enhancer located >1 Mb 3' of human and murine Myc that binds Notch transcription complexes and physically interacts with the Myc proximal promoter. The Notch1 binding element in this region activates reporter genes in a Notch-dependent, cell-context-specific fashion that requires a conserved Notch complex binding site. Acute changes in Notch activation produce rapid changes in H3K27 acetylation across the entire enhancer (a region spanning >600 kb) that correlate with Myc expression. This broad Notch-influenced region comprises an enhancer region containing multiple domains, recognizable as discrete H3K27 acetylation peaks. Leukemia cells selected for resistance to Notch inhibitors express Myc despite epigenetic silencing of enhancer domains near the Notch transcription complex binding sites. Notch-independent expression of Myc in resistant cells is highly sensitive to inhibitors of bromodomain containing 4 (Brd4), a change in drug sensitivity that is accompanied by preferential association of the Myc promoter with more 3' enhancer domains that are strongly dependent on Brd4 for function. These findings indicate that altered long-range enhancer activity can mediate resistance to targeted therapies and provide a mechanistic rationale for combined targeting of Notch and Brd4 in leukemia.

Tissue specific CTCF occupancy and boundary function at the human growth hormone locus.

The robust and tissue-specific activation of the human growth hormone (hGH) gene cluster in the pituitary and placenta constitutes an informative model for analysis of gene regulation. The five-gene hGH cluster is regulated by two partially overlapping sets of DNase I hypersensitive sites (HSs) that constitute the pituitary (HSI, II, III and V) and placental (HSIII, IV, and V) locus control regions (LCRs). The single placenta-specific LCR component, HSIV, is located at -30 kb to the cluster. Here we generate a series of hGH/BAC transgenes specifically modified to identify structural features of the hGH locus required for its appropriate placental expression. We find that placental specificity is dependent on the overall multigene configuration of the cluster whereas the distance between the cluster and its LCR impacts the level of placental expression. We further observe that a major function of the placental hGH LCR is to insulate the transgene locus from site-of-integration effects. This insulation activity is linked to placenta-specific occupancy of the chromatin architectural protein, CTCF, at HSIV. These data reveal a remarkable combination of structural configurations and regulatory determinants that must work in concert to insure robust and tightly controlled expression from a complex multigene locus.