CeO2 nanoparticle dose and exposure modulate soybean development and plant-mediated responses in root-associated bacterial communities.

Agricultural soils are increasingly undergoing inadvertent and purposeful exposures to engineered CeO2 nanoparticles (NPs), which can impact crops and root-associated microbial communities. However, interactions between NP concentration and exposure duration on plant-mediated responses of root-associated bacterial communities are not well understood. Soybeans seedlings were grown in soil with uncoated NPs added at concentrations of 0, 1 or 100 mg kg-1. Total soil exposure durations were either 190 days, starting 106 days before planting or 84 days with NP amendments coinciding with planting. We assessed plant development, bacterial diversity, differential abundance and inferred functional changes across rhizosphere, rhizoplane, and root tissue compartments. Plant non-monotonic dose responses were mirrored in bacterial communities. Most notably, effects were magnified in the rhizoplane under low-dose, short-exposures. Enriched metabolic pathways were primarily related to biosynthesis and degradation/utilization/assimilation, rather than responses to metals or oxidative stress. Our results indicate that plant-mediated bacterial responses were greater than direct NP impacts. Also, we identify needs for modeling non-monotonic legume stress responses that account for coinfection with mutualistic and parasitic bacteroids. Our findings provide new insights regarding effects of applications of soil amendments such as biosolids containing NPs or nano-enabled formulations used in cultivation of legumes and other crops.

A mouse model of wildfire smoke-induced health effects: sex differences in acute and sustained effects of inhalation exposures.

Due to climate change, wildfires have increased in intensity and duration. While wildfires threaten lives directly, the smoke has more far-reaching adverse health impacts. During an extreme 2017 wildfire event, residents of Seeley Lake, Montana were exposed to unusually high levels of wood smoke (WS) causing sustained effects on lung function (decreased FEV1/FVC). Objective: The present study utilized an animal model of WS exposure to research cellular and molecular mechanisms of the resulting health effects. Methods: Mice were exposed to inhaled WS utilizing locally harvested wood to recapitulate community exposures. WS was generated at a rate resulting in a 5 mg/m3 PM2.5 exposure for five days. Results: This exposure resulted in a similar 0.28 mg/m2 particle deposition (lung surface area) in mice that was calculated for human exposure. As with the community observations, there was a significant effect on lung function, increased resistance, and decreased compliance, that was more pronounced in males at an extended (2 months) timepoint and males were more affected than females: ex vivo assays illustrated changes to alveolar macrophage functions (increased TNFα secretion and decreased efferocytosis). Female mice had significantly elevated IL-33 levels in lungs, however, pretreatment of male mice with IL-33 resulted in an abrogation of the observed WS effects, suggesting a dose-dependent role of IL-33. Additionally, there were greater immunotoxic effects in male mice. Discussion: These findings replicated the outcomes in humans and suggest that IL-33 is involved in a mechanism of the adverse effects of WS exposures that inform on potential sex differences.

Limited expression of Nrf2 in neurons across the central nervous system.

Nrf2, encoded by the gene Nfe2l2, is a broadly expressed transcription factor that regulates gene expression in response to reactive oxygen species (ROS) and oxidative stress. It is commonly referred to as a ubiquitous pathway, but this generalization overlooks work indicating that Nrf2 is essentially unexpressed in some neuronal populations. To explore whether this pattern extends throughout the central nervous system (CNS), we quantified Nfe2l2 expression and chromatin accessibility at the Nfe2l2 locus across multiple single cell datasets. In both the mouse and human CNS, Nfe2l2 was repressed in almost all mature neurons, but highly expressed in non-neuronal support cells, and this pattern was robust across multiple human CNS diseases. A subset of key Nrf2 target genes, like Slc7a11, also remained low in neurons. Thus, these data suggest that while most cells express Nfe2l2, with activity determined by ROS levels, neurons actively avoid Nrf2 activity by keeping Nfe2l2 expression low.

An antioxidant response element regulates the HIF1α axis in breast cancer cells.

The redox sensitive transcription factor NRF2 is a central regulator of the transcriptional response to reactive oxygen species (ROS). NRF2 is widely recognized for its ROS-responsive upregulation of antioxidant genes that are essential for mitigating the damaging effects of oxidative stress. However, multiple genome-wide approaches have suggested that NRF2's regulatory reach extends well beyond the canonical antioxidant genes, with the potential to regulate many noncanonical target genes. Recent work from our lab and others suggests HIF1A, which encodes the hypoxia-responsive transcription factor HIF1α, is one such noncanonical NRF2 target. These studies found that NRF2 activity is associated with high HIF1A expression in multiple cellular contexts, HIF1A expression is partially dependent on NRF2, and there is a putative NRF2 binding site (antioxidant response element, or ARE) approximately 30 kilobases upstream of HIF1A. These findings all support a model in which HIF1A is a direct target of NRF2, but did not confirm the functional importance of the upstream ARE in HIF1A expression. Here we use CRISPR/Cas9 genome editing to mutate this ARE in its genomic context and test the impact on HIF1A expression. We find that mutation of this ARE in a breast cancer cell line (MDA-MB-231) eliminates NRF2 binding and decreases HIF1A expression at the transcript and protein levels, and disrupts HIF1α target genes as well as phenotypes driven by these HIF1α targets. Taken together, these results indicate that this NRF2 targeted ARE plays an important role in the expression of HIF1A and activity of the HIF1α axis in MDA-MB-231 cells.

Limited Expression of Nrf2 in Neurons Across the Central Nervous System.

Nrf2 is a broadly expressed transcription factor that regulates gene expression in response to reactive oxygen species (ROS) and oxidative stress. It is commonly referred to as a ubiquitous pathway, but this generalization overlooks work indicating that Nrf2 is essentially unexpressed in some neuronal populations. To explore whether this pattern extends throughout the central nervous system (CNS), we quantified Nrf2 expression and chromatin accessibility at the Nrf2 locus across multiple single cell datasets. In both the mouse and human CNS, Nrf2 was repressed in almost all mature neurons, but highly expressed in non-neuronal support cells, and this pattern was robust across multiple human CNS diseases. A subset of key Nrf2 target genes, like Slc7a11 , also remained low in neurons. Thus, these data suggest that while most cells express Nrf2, with activity determined by ROS levels, neurons actively avoid Nrf2 activity by keeping Nrf2 expression low.

Ribbon boosts ribosomal protein gene expression to coordinate organ form and function.

Cell growth is well defined for late (postembryonic) stages of development, but evidence for early (embryonic) cell growth during postmitotic morphogenesis is limited. Here, we report early cell growth as a key characteristic of tubulogenesis in the Drosophila embryonic salivary gland (SG) and trachea. A BTB/POZ domain nuclear factor, Ribbon (Rib), mediates this early cell growth. Rib binds the transcription start site of nearly every SG-expressed ribosomal protein gene (RPG) and is required for full expression of all RPGs tested. Rib binding to RPG promoters in vitro is weak and not sequence specific, suggesting that specificity is achieved through cofactor interactions. Accordingly, we demonstrate Rib's ability to physically interact with each of the three known regulators of RPG transcription. Surprisingly, Rib-dependent early cell growth in another tubular organ, the embryonic trachea, is not mediated by direct RPG transcription. These findings support a model of early cell growth customized by transcriptional regulatory networks to coordinate organ form and function.

Transcriptional reprogramming by oxidative stress occurs within a predefined chromatin accessibility landscape.

Reactive oxygen species (ROS) are important signaling molecules in many physiological processes, yet excess ROS leads to cell damage and can lead to pathology. Accordingly, cells need to maintain tight regulation of ROS levels, and ROS-responsive transcriptional reprogramming is central to this process. Although it has long been recognized that oxidative stress leads to rapid, significant changes in gene expression, the impact of oxidative stress on the underlying chromatin accessibility landscape remained unclear. Here, we asked whether ROS-responsive transcriptional reprogramming is accompanied by reprogramming of the chromatin environment in MCF7 human breast cancer cells. Using a time-course exposure to multiple inducers of oxidative stress, we determined that the widespread ROS-responsive changes in gene expression induced by ROS occur with minimal changes to the chromatin environment. While we did observe changes in chromatin accessibility, these changes were: (1) far less numerous than gene expression changes after oxidative stress, and (2) occur within pre-existing regions of accessible chromatin. Transcription factor (TF) footprinting analysis of our ATAC-seq experiments identified 5 TFs or TF families with evidence for ROS-responsive changes in DNA binding: NRF2, AP-1, p53, NFY, and SP/KLF. Importantly, several of these (AP-1, NF-Y, and SP/KLF factors) have not been previously implicated as widespread regulators in the response to ROS. In summary, we have characterized genome-wide changes in gene expression and chromatin accessibility in response to ROS treatment of MCF7 cells, and we have found that regulation of the large-scale transcriptional response to excess ROS is primarily constrained by the cell's pre-existing chromatin landscape.

New regulators of Drosophila eye development identified from temporal transcriptome changes.

In the last larval instar, uncommitted progenitor cells in the Drosophila eye primordium start to adopt individual retinal cell fates, arrest their growth and proliferation, and initiate terminal differentiation into photoreceptor neurons and other retinal cell types. To explore the regulation of these processes, we have performed mRNA-Seq studies of the larval eye and antennal primordial at multiple developmental stages. A total of 10,893 fly genes were expressed during these stages and could be adaptively clustered into gene groups, some of whose expression increases or decreases in parallel with the cessation of proliferation and onset of differentiation. Using in situ hybridization of a sample of 98 genes to verify spatial and temporal expression patterns, we estimate that 534 genes or more are transcriptionally upregulated during retinal differentiation, and 1367 or more downregulated as progenitor cells differentiate. Each group of co-expressed genes is enriched for regulatory motifs recognized by co-expressed transcription factors, suggesting that they represent coherent transcriptional regulatory programs. Using available mutant strains, we describe novel roles for the transcription factors SoxNeuro (SoxN), H6-like homeobox (Hmx), CG10253, without children (woc), Structure specific recognition protein (Ssrp), and multisex combs (mxc).

CrebA increases secretory capacity through direct transcriptional regulation of the secretory machinery, a subset of secretory cargo, and other key regulators.

Specialization of many cells, including the acinar cells of the salivary glands and pancreas, milk-producing cells of mammary glands, mucus-secreting goblet cells, antibody-producing plasma cells, and cells that generate the dense extracellular matrices of bone and cartilage, requires scaling up both secretory machinery and cell-type specific secretory cargo. Using tissue-specific genome-scale analyses, we determine how increases in secretory capacity are coordinated with increases in secretory load in the Drosophila salivary gland (SG), an ideal model for gaining mechanistic insight into the functional specialization of secretory organs. Our findings show that CrebA, a bZIP transcription factor, directly binds genes encoding the core secretory machinery, including protein components of the signal recognition particle and receptor, ER cargo translocators, Cop I and Cop II vesicles, as well as the structural proteins and enzymes of these organelles. CrebA directly binds a subset of SG cargo genes and CrebA binds and boosts expression of Sage, a SG-specific transcription factor essential for cargo expression. To further enhance secretory output, CrebA binds and activates Xbp1 and Tudor-SN. Thus, CrebA directly upregulates the machinery of secretion and additional factors to increase overall secretory capacity in professional secretory cells; concomitant increases in cargo are achieved both directly and indirectly.

The Drosophila MLR COMPASS complex is essential for programming cis-regulatory information and maintaining epigenetic memory during development.

The MLR COMPASS complex monomethylates H3K4 that serves to epigenetically mark transcriptional enhancers to drive proper gene expression during animal development. Chromatin enrichment analyses of the Drosophila MLR complex reveals dynamic association with promoters and enhancers in embryos with late stage enrichments biased toward both active and poised enhancers. RNAi depletion of the Cmi (also known as Lpt) subunit that contains the chromatin binding PHD finger domains attenuates enhancer functions, but unexpectedly results in inappropriate enhancer activation during stages when hormone responsive enhancers are poised, revealing critical epigenetic roles involved in both the activation and repression of enhancers depending on developmental context. Cmi is necessary for robust H3K4 monomethylation and H3K27 acetylation that mark active enhancers, but not for the chromatin binding of Trr, the MLR methyltransferase. Our data reveal two likely major regulatory modes of MLR function, contributions to enhancer commissioning in early embryogenesis and bookmarking enhancers to enable rapid transcriptional re-activation at subsequent developmental stages.

Pesticide Encapsulation at the Nanoscale Drives Changes to the Hydrophobic Partitioning and Toxicity of an Active Ingredient.

Given the costs associated with designing novel active ingredients, new formulations focus on the use of other ingredients to modify existing formulations. Nanosized encapsulated pesticides offer a variety of enhanced features including controlled release and improved efficacy. Despite the presence of nanosized capsules in current-use pesticide formulations, the analytical and toxicological implications of encapsulation are uncertain. To explore this issue quantitatively, we fractionated the capsules of a commercially available encapsulated insecticide formulation (γ-cyhalothrin active ingredient) into two size ranges: a large fraction (LF), with an average hydrodynamic diameter (HDD) of 758 nm, and a small fraction (SF), with an average HDD of 449 nm. We developed a novel extraction method demonstrating a time-dependent inhibition of γ-cyhalothrin from capsules for up to 48 h. An acute immobilization test with a freshwater macroinvertebrate (Ceriodaphnia dubia) revealed that the SF was significantly more toxic than both the LF and the free γ-cyhalothrin treatment (EC50 = 0.18 µg/L, 0.57 µg/L, and 0.65 µg/L, respectively). These findings highlight that encapsulation of γ-cyhalothrin mitigates hydrophobic partitioning in a time-dependent manner and influences toxicity in a size-dependent manner. Recognizing the analytical and toxicological nuances of various nanosized capsules can contribute to innovation in pesticide formulations and may lead to more comprehensive pesticide regulation.

Identification of a functional antioxidant response element at the HIF1A locus.

Reactive oxygen species (ROS), which are a byproduct of oxidative metabolism, serve as signaling molecules in a number of physiological settings. However, if their levels are not tightly maintained, excess ROS lead to potentially cytotoxic oxidative stress. Accordingly, several transcriptional regulatory networks have evolved to include components that are highly ROS-responsive. Depending on the context, these regulatory networks can leverage ROS to respond to nutrient conditions, metabolism, or other physiological signals, or to respond to oxidative stress. However, ROS signaling is complex, so regulatory interactions between various ROS-responsive transcription factors are still being mapped out. Here we show that the transcription factor NRF2, a key regulator of the adaptive response to oxidative stress, directly regulates expression of HIF1A, which encodes HIF1α, a key transcriptional regulator of the adaptive response to hypoxia. We used an integrative genomics approach to identify HIF1A as a ROS-responsive transcript and we found an NRF2-bound antioxidant response element (ARE) approximately 30 kilobases upstream of HIF1A. This ARE sequence is deeply conserved, and we verified that it is directly bound and activated by NRF2. In addition, we found that HIF1A is upregulated in breast and bladder tumors with high NRF2 activity. Taken together, our results demonstrate that NRF2 targets a functional ARE at the HIF1A locus, and reveal a direct regulatory connection between two important oxygen responsive transcription factors.

A distinct class of antioxidant response elements is consistently activated in tumors with NRF2 mutations.

NRF2 is a redox-responsive transcription factor that regulates expression of cytoprotective genes via its interaction with DNA sequences known as antioxidant response elements (AREs). NRF2 activity is induced by oxidative stress, but oxidative stress is not the only context in which NRF2 can be activated. Mutations that disrupt the interaction between NRF2 and KEAP1, an inhibitor of NRF2, lead to NRF2 hyperactivation and promote oncogenesis. The mechanisms underlying NRF2's oncogenic properties remain unclear, but likely involve aberrant expression of select NRF2 target genes. We tested this possibility using an integrative genomics approach to get a precise view of the direct NRF2 target genes dysregulated in tumors with NRF2 hyperactivating mutations. This approach revealed a core set of 32 direct NRF2 targets that are consistently upregulated in NRF2 hyperactivated tumors. This set of NRF2 "cancer target genes" includes canonical redox-related NRF2 targets, as well as target genes that have not been previously linked to NRF2 activation. Importantly, NRF2-driven upregulation of this gene set is largely independent of the organ system where the tumor developed. One key distinguishing feature of these NRF2 cancer target genes is that they are regulated by high affinity AREs that fall within genomic regions possessing a ubiquitously permissive chromatin signature. This implies that these NRF2 cancer target genes are responsive to oncogenic NRF2 in most tissues because they lack the regulatory constraints that restrict expression of most other NRF2 target genes. This NRF2 cancer target gene set also serves as a reliable proxy for NRF2 activity, and high NRF2 activity is associated with significant decreases in survival in multiple cancer types. Overall, the pervasive upregulation of these NRF2 cancer targets across multiple cancers, and their association with negative outcomes, suggests that these will be central to dissecting the functional implications of NRF2 hyperactivation in several cancer contexts.

cis-regulatory architecture of a short-range EGFR organizing center in the Drosophila melanogaster leg.

We characterized the establishment of an Epidermal Growth Factor Receptor (EGFR) organizing center (EOC) during leg development in Drosophila melanogaster. Initial EGFR activation occurs in the center of leg discs by expression of the EGFR ligand Vn and the EGFR ligand-processing protease Rho, each through single enhancers, vnE and rhoE, that integrate inputs from Wg, Dpp, Dll and Sp1. Deletion of vnE and rhoE eliminates vn and rho expression in the center of the leg imaginal discs, respectively. Animals with deletions of both vnE and rhoE (but not individually) show distal but not medial leg truncations, suggesting that the distal source of EGFR ligands acts at short-range to only specify distal-most fates, and that multiple additional 'ring' enhancers are responsible for medial fates. Further, based on the cis-regulatory logic of vnE and rhoE we identified many additional leg enhancers, suggesting that this logic is broadly used by many genes during Drosophila limb development.

A hypermorphic antioxidant response element is associated with increased MS4A6A expression and Alzheimer's disease.

Late onset Alzheimer's disease (AD) is a multifactorial disorder, with AD risk influenced by both environmental and genetic factors. Recent genome-wide association studies (GWAS) have identified genetic loci associated with increased risk of developing AD. The MS4A (membrane-spanning 4-domains subfamily A) gene cluster is one of the most significant loci associated with AD risk, and MS4A6A expression is correlated with AD pathology. We identified a single nucleotide polymorphism, rs667897, at the MS4A locus that creates an antioxidant response element and links MS4A6A expression to the stress responsive Cap-n-Collar (CNC) transcription factors NRF1 (encoded by NFE2L1) and NRF2 (encoded by NFE2L2). The risk allele of rs667897 generates a strong CNC binding sequence that is activated by proteostatic stress in an NRF1-dependent manner, and is associated with increased expression of the gene MS4A6A. Together, these findings suggest that the cytoprotective CNC regulatory network aberrantly activates MS4A6A expression and increases AD risk in a subset of the population.

A Role for Monomethylation of Histone H3-K27 in Gene Activity in Drosophila.

Polycomb repressive complex 2 (PRC2) is a conserved chromatin-modifying enzyme that methylates histone H3 on lysine-27 (K27). PRC2 can add one, two, or three methyl groups and the fully methylated product, H3-K27me3, is a hallmark of Polycomb-silenced chromatin. Less is known about functions of K27me1 and K27me2 and the dynamics of flux through these states. These modifications could serve mainly as intermediates to produce K27me3 or they could each convey distinct epigenetic information. To investigate this, we engineered a variant of Drosophila melanogaster PRC2 which is converted into a monomethyltransferase. A single substitution, F738Y, in the lysine-substrate binding pocket of the catalytic subunit, E(Z), creates an enzyme that retains robust K27 monomethylation but dramatically reduced di- and trimethylation. Overexpression of E(Z)-F738Y in fly cells triggers desilencing of Polycomb target genes significantly more than comparable overexpression of catalytically deficient E(Z), suggesting that H3-K27me1 contributes positively to gene activity. Consistent with this, normal genomic distribution of H3-K27me1 is enriched on actively transcribed Drosophila genes, with localization overlapping the active H3-K36me2/3 chromatin marks. Thus, distinct K27 methylation states link to either repression or activation depending upon the number of added methyl groups. If so, then H3-K27me1 deposition may involve alternative methyltransferases beyond PRC2, which is primarily repressive. Indeed, assays on fly embryos with PRC2 genetically inactivated, and on fly cells with PRC2 chemically inhibited, show that substantial H3-K27me1 accumulates independently of PRC2. These findings imply distinct roles for K27me1 vs. K27me3 in transcriptional control and an expanded machinery for methylating H3-K27.

Phylogenetic context determines the role of competition in adaptive radiation.

Understanding ecological mechanisms regulating the evolution of biodiversity is of much interest to ecologists and evolutionary biologists. Adaptive radiation constitutes an important evolutionary process that generates biodiversity. Competition has long been thought to influence adaptive radiation, but the directionality of its effect and associated mechanisms remain ambiguous. Here, we report a rigorous experimental test of the role of competition on adaptive radiation using the rapidly evolving bacterium Pseudomonas fluorescens SBW25 interacting with multiple bacterial species that differed in their phylogenetic distance to the diversifying bacterium. We showed that the inhibitive effect of competitors on the adaptive radiation of P. fluorescens decreased as their phylogenetic distance increased. To explain this phylogenetic dependency of adaptive radiation, we linked the phylogenetic distance between P. fluorescens and its competitors to their niche and competitive fitness differences. Competitive fitness differences, which showed weak phylogenetic signal, reduced P. fluorescens abundance and thus diversification, whereas phylogenetically conserved niche differences promoted diversification. These results demonstrate the context dependency of competitive effects on adaptive radiation, and highlight the importance of past evolutionary history for ongoing evolutionary processes.

A Polymorphic Antioxidant Response Element Links NRF2/sMAF Binding to Enhanced MAPT Expression and Reduced Risk of Parkinsonian Disorders.

The NRF2/sMAF protein complex regulates the oxidative stress response by occupying cis-acting enhancers containing an antioxidant response element (ARE). Integrating genome-wide maps of NRF2/sMAF occupancy with disease-susceptibility loci, we discovered eight polymorphic AREs linked to 14 highly ranked disease-risk SNPs in individuals of European ancestry. Among these SNPs was rs242561, located within a regulatory region of the MAPT gene (encoding microtubule-associated protein Tau). It was consistently occupied by NRF2/sMAF in multiple experiments and its strong-binding allele associated with higher mRNA levels in cell lines and human brain tissue. Induction of MAPT transcription by NRF2 was confirmed using a human neuroblastoma cell line and a Nrf2-deficient mouse model. Most importantly, rs242561 displayed complete linkage disequilibrium with a highly protective allele identified in multiple GWASs of progressive supranuclear palsy, Parkinson's disease, and corticobasal degeneration. These observations suggest a potential role for NRF2/sMAF in tauopathies and a possible role for NRF2 pathway activators in disease prevention.

DNA-binding sequence specificity of DUX4.

BACKGROUND: Misexpression of the double homeodomain transcription factor DUX4 results in facioscapulohumeral muscular dystrophy (FSHD). A DNA-binding consensus with two tandem TAAT motifs based on chromatin IP peaks has been discovered; however, the consensus has multiple variations (flavors) of unknown relative activity. In addition, not all peaks have this consensus, and the Pitx1 promoter, the first DUX4 target sequence mooted, has a different TAAT-rich sequence. Furthermore, it is not known whether and to what extent deviations from the consensus affect DNA-binding affinity and transcriptional activation potential. RESULTS: Here, we take both unbiased and consensus sequence-driven approaches to determine the DNA-binding specificity of DUX4 and its tolerance to mismatches at each site within its consensus sequence. We discover that the best binding and the greatest transcriptional activation are observed when the two TAAT motifs are separated by a C residue. The second TAAT motif in the consensus sequence is actually (T/C)AAT. We find that a T is preferred here. DUX4 has no transcriptional activity on "half-sites", i.e., those bearing only a single TAAT motif. We further find that DUX4 does not bind to the TAATTA motif in the Pitx1 promoter, that Pitx1 sequences have no competitive band shift activity, and that the Pitx1 sequence is transcriptionally inactive, calling into question PITX1 as a DUX4 target gene. Finally, by multimerizing binding sites, we find that DUX4 transcriptional activation demonstrates tremendous synergy and that at low DNA concentrations, at least two motifs are necessary to detect a transcriptional response. CONCLUSIONS: These studies illuminate the DNA-binding sequence preferences of DUX4.

Ribbon regulates morphogenesis of the Drosophila embryonic salivary gland through transcriptional activation and repression.

Transcription factors affect spatiotemporal patterns of gene expression often regulating multiple aspects of tissue morphogenesis, including cell-type specification, cell proliferation, cell death, cell polarity, cell shape, cell arrangement and cell migration. In this work, we describe a distinct role for Ribbon (Rib) in controlling cell shape/volume increases during elongation of the Drosophila salivary gland (SG). Notably, the morphogenetic changes in rib mutants occurred without effects on general SG cell attributes such as specification, proliferation and apoptosis. Moreover, the changes in cell shape/volume in rib mutants occurred without compromising epithelial-specific morphological attributes such as apicobasal polarity and junctional integrity. To identify the genes regulated by Rib, we performed ChIP-seq analysis in embryos driving expression of GFP-tagged Rib specifically in the SGs. To learn if the Rib binding sites identified in the ChIP-seq analysis were linked to changes in gene expression, we performed microarray analysis comparing RNA samples from age-matched wild-type and rib null embryos. From the superposed ChIP-seq and microarray gene expression data, we identified 60 genomic sites bound by Rib likely to regulate SG-specific gene expression. We confirmed several of the identified Rib targets by qRT-pCR and/or in situ hybridization. Our results indicate that Rib regulates cell growth and tissue shape in the Drosophila salivary gland via a diverse array of targets through both transcriptional activation and repression. Furthermore, our results suggest that autoregulation of rib expression may be a key component of the SG morphogenetic gene network.