Cellular and clinical impact of protein phosphatase enzyme epigenetic silencing in multiple cancer tissues.

BACKGROUND: Protein Phosphatase Enzymes (PPE) and protein kinases simultaneously control phosphorylation mechanisms that tightly regulate intracellular signalling pathways and stimulate cellular responses. In human malignancies, PPE and protein kinases are frequently mutated resulting in uncontrolled kinase activity and PPE suppression, leading to cell proliferation, migration and resistance to anti-cancer therapies. Cancer associated DNA hypermethylation at PPE promoters gives rise to transcriptional silencing (epimutations) and is a hallmark of cancer. Despite recent advances in sequencing technologies, data availability and computational capabilities, only a fraction of PPE have been reported as transcriptionally inactive as a consequence of epimutations. METHODS: In this study, we examined promoter-associated DNA methylation profiles in Protein Phosphatase Enzymes and their Interacting Proteins (PPEIP) in a cohort of 705 cancer patients in five tissues (Large intestine, Oesophagus, Lung, Pancreas and Stomach) in three cell models (primary tumours, cancer cell lines and 3D embedded cancer cell cultures). As a subset of PPEIP are known tumour suppressor genes, we analysed the impact of PPEIP promoter hypermethylation marks on gene expression, cellular networks and in a clinical setting. RESULTS: Here, we report epimutations in PPEIP are a frequent occurrence in the cancer genome and manifest independent of transcriptional activity. We observed that different tumours have varying susceptibility to epimutations and identify specific cellular signalling networks that are primarily affected by epimutations. Additionally, RNA-seq analysis showed the negative impact of epimutations on most (not all) Protein Tyrosine Phosphatase transcription. Finally, we detected novel clinical biomarkers that inform on patient mortality and anti-cancer treatment sensitivity. CONCLUSIONS: We propose that DNA hypermethylation marks at PPEIP frequently contribute to the pathogenesis of malignancies and within the precision medicine space, hold promise as biomarkers to inform on clinical features such as patient survival and therapeutic response.

Epigenetic loss of the transfer RNA-modifying enzyme TYW2 induces ribosome frameshifts in colon cancer.

Transfer RNA (tRNA) activity is tightly regulated to provide a physiological protein translation, and tRNA chemical modifications control its function in a complex with ribosomes and messenger RNAs (mRNAs). In this regard, the correct hypermodification of position G37 of phenylalanine-tRNA, adjacent to the anticodon, is critical to prevent ribosome frameshifting events. Here we report that the tRNA-yW Synthesizing Protein 2 (TYW2) undergoes promoter hypermethylation-associated transcriptional silencing in human cancer, particularly in colorectal tumors. The epigenetic loss of TYW2 induces guanosine hypomodification in phenylalanine-tRNA, an increase in -1 ribosome frameshift events, and down-regulation of transcripts by mRNA decay, such as of the key cancer gene ROBO1. Importantly, TYW2 epigenetic inactivation is linked to poor overall survival in patients with early-stage colorectal cancer, a finding that could be related to the observed acquisition of enhanced migration properties and epithelial-to-mesenchymal features in the colon cancer cells that harbor TYW2 DNA methylation-associated loss. These findings provide an illustrative example of how epigenetic changes can modify the epitranscriptome and further support a role for tRNA modifications in cancer biology.

Characterization of Acute HCV Infection and Transmission Networks in People Who Currently Inject Drugs in Catalonia: Usefulness of Dried Blood Spots.

BACKGROUND AND AIMS: Accurate identification of recent HCV infections is critical for tracing the extent and mechanisms of ongoing transmission. We aimed to validate dried blood spot (DBS) samples for the assessment of Hepatitis C virus (HCV) genetic diversity and to determine epidemiological parameters including incidence, determinants of acute infection, and phylogenetic clustering in people who inject drugs (PWID). APPROACH AND RESULTS: HCV nonstructural protein 5B next-generation sequencing was performed from plasma and/or DBS in 220 viremic PWID from the HepCdetect II study. No significant differences were found in consensus sequences or Shannon entropy (SE) intrahost diversity estimate between paired plasma/DBS specimens. SE values were used to identify acute infections with 93.3% sensitivity (95% CI, 0.81-1.06) and 95.0% specificity (95% CI, 0.88-1.02) in a set of well-defined controls. An acute HCV infection (either primary infection or reinfection) was detected in 13.5% of viremic participants and was associated with age ≤30 years (OR, 8.09), injecting less than daily (OR, 4.35), ≤5 years of injected drug use (OR, 3.43), sharing cocaine snorting straws (OR, 2.89), and being unaware of their HCV status (OR, 3.62). Annualized HCV incidence was estimated between 31 and 59/100 person-years. On phylogenetic analysis, 46.8% of viremic cases were part of a transmission pair or cluster; age ≤30 years (OR, 6.16), acute infection (OR, 5.73), and infection with subtype 1a (OR, 4.78) were independently associated with this condition. CONCLUSIONS: The results obtained from plasma and DBS characterize PWID with acute infection and those involved in ongoing HCV transmission and allow estimating incidence from cross-sectional data. This information is critical for the design and assessment of targeted harm reduction programs and test-and-treat interventions and to facilitate monitoring of HCV elimination in this key population.

Epigenetic footprint enables molecular risk stratification of hepatoblastoma with clinical implications.

BACKGROUND & AIMS: Hepatoblastoma (HB) is a rare disease. Nevertheless, it is the predominant pediatric liver cancer, with limited therapeutic options for patients with aggressive tumors. Herein, we aimed to uncover the mechanisms of HB pathobiology and to identify new biomarkers and therapeutic targets in a move towards precision medicine for patients with advanced HB. METHODS: We performed a comprehensive genomic, transcriptomic and epigenomic characterization of 159 clinically annotated samples from 113 patients with HB, using high-throughput technologies. RESULTS: We discovered a widespread epigenetic footprint of HB that includes hyperediting of the tumor suppressor BLCAP concomitant with a genome-wide dysregulation of RNA editing and the overexpression of mainly non-coding genes of the oncogenic 14q32 DLK1-DIO3 locus. By unsupervised analysis, we identified 2 epigenomic clusters (Epi-CA, Epi-CB) with distinct degrees of DNA hypomethylation and CpG island hypermethylation that are associated with the C1/C2/C2B transcriptomic subtypes. Based on these findings, we defined the first molecular risk stratification of HB (MRS-HB), which encompasses 3 main prognostic categories and improves the current clinical risk stratification approach. The MRS-3 category (28%), defined by strong 14q32 locus expression and Epi-CB methylation features, was characterized by CTNNB1 and NFE2L2 mutations, a progenitor-like phenotype and clinical aggressiveness. Finally, we identified choline kinase alpha as a promising therapeutic target for intermediate and high-risk HBs, as its inhibition in HB cell lines and patient-derived xenografts strongly abrogated tumor growth. CONCLUSIONS: These findings provide a detailed insight into the molecular features of HB and could be used to improve current clinical stratification approaches and to develop treatments for patients with HB. LAY SUMMARY: Hepatoblastoma is a rare childhood liver cancer that has been understudied. We have used cutting-edge technologies to expand our molecular knowledge of this cancer. Our biological findings can be used to improve clinical management and pave the way for the development of novel therapies for this cancer.

Epigenetic loss of RNA-methyltransferase NSUN5 in glioma targets ribosomes to drive a stress adaptive translational program.

Tumors have aberrant proteomes that often do not match their corresponding transcriptome profiles. One possible cause of this discrepancy is the existence of aberrant RNA modification landscapes in the so-called epitranscriptome. Here, we report that human glioma cells undergo DNA methylation-associated epigenetic silencing of NSUN5, a candidate RNA methyltransferase for 5-methylcytosine. In this setting, NSUN5 exhibits tumor-suppressor characteristics in vivo glioma models. We also found that NSUN5 loss generates an unmethylated status at the C3782 position of 28S rRNA that drives an overall depletion of protein synthesis, and leads to the emergence of an adaptive translational program for survival under conditions of cellular stress. Interestingly, NSUN5 epigenetic inactivation also renders these gliomas sensitive to bioactivatable substrates of the stress-related enzyme NQO1. Most importantly, NSUN5 epigenetic inactivation is a hallmark of glioma patients with long-term survival for this otherwise devastating disease.

Multitrait genome association analysis identifies new susceptibility genes for human anthropometric variation in the GCAT cohort.

BACKGROUND: Heritability estimates have revealed an important contribution of SNP variants for most common traits; however, SNP analysis by single-trait genome-wide association studies (GWAS) has failed to uncover their impact. In this study, we applied a multitrait GWAS approach to discover additional factor of the missing heritability of human anthropometric variation. METHODS: We analysed 205 traits, including diseases identified at baseline in the GCAT cohort (Genomes For Life- Cohort study of the Genomes of Catalonia) (n=4988), a Mediterranean adult population-based cohort study from the south of Europe. We estimated SNP heritability contribution and single-trait GWAS for all traits from 15 million SNP variants. Then, we applied a multitrait-related approach to study genome-wide association to anthropometric measures in a two-stage meta-analysis with the UK Biobank cohort (n=336 107). RESULTS: Heritability estimates (eg, skin colour, alcohol consumption, smoking habit, body mass index, educational level or height) revealed an important contribution of SNP variants, ranging from 18% to 77%. Single-trait analysis identified 1785 SNPs with genome-wide significance threshold. From these, several previously reported single-trait hits were confirmed in our sample with LINC01432 (p=1.9×10-9) variants associated with male baldness, LDLR variants with hyperlipidaemia (ICD-9:272) (p=9.4×10-10) and variants in IRF4 (p=2.8×10-57), SLC45A2 (p=2.2×10-130), HERC2 (p=2.8×10-176), OCA2 (p=2.4×10-121) and MC1R (p=7.7×10-22) associated with hair, eye and skin colour, freckling, tanning capacity and sun burning sensitivity and the Fitzpatrick phototype score, all highly correlated cross-phenotypes. Multitrait meta-analysis of anthropometric variation validated 27 loci in a two-stage meta-analysis with a large British ancestry cohort, six of which are newly reported here (p value threshold <5×10-9) at ZRANB2-AS2, PIK3R1, EPHA7, MAD1L1, CACUL1 and MAP3K9. CONCLUSION: Considering multiple-related genetic phenotypes improve associated genome signal detection. These results indicate the potential value of data-driven multivariate phenotyping for genetic studies in large population-based cohorts to contribute to knowledge of complex traits.

Inosine modifications in human tRNAs are incorporated at the precursor tRNA level.

Transfer RNAs (tRNAs) are key adaptor molecules of the genetic code that are heavily modified post-transcriptionally. Inosine at the first residue of the anticodon (position 34; I34) is an essential widespread tRNA modification that has been poorly studied thus far. The modification in eukaryotes results from a deamination reaction of adenine that is catalyzed by the heterodimeric enzyme adenosine deaminase acting on tRNA (hetADAT), composed of two subunits: ADAT2 and ADAT3. Using high-throughput small RNA sequencing (RNAseq), we show that this modification is incorporated to human tRNAs at the precursor tRNA level and during maturation. We also functionally validated the human genes encoding for hetADAT and show that the subunits of this enzyme co-localize in nucleus in an ADAT2-dependent manner. Finally, by knocking down HsADAT2, we demonstrate that variations in the cellular levels of hetADAT will result in changes in the levels of I34 modification in all its potential substrates. Altogether, we present RNAseq as a powerful tool to study post-transcriptional tRNA modifications at the precursor tRNA level and give the first insights on the biology of I34 tRNA modification in metazoans.

Human mitochondrial disease-like symptoms caused by a reduced tRNA aminoacylation activity in flies.

The translation of genes encoded in the mitochondrial genome requires specific machinery that functions in the organelle. Among the many mutations linked to human disease that affect mitochondrial translation, several are localized to nuclear genes coding for mitochondrial aminoacyl-transfer RNA synthetases. The molecular significance of these mutations is poorly understood, but it is expected to be similar to that of the mutations affecting mitochondrial transfer RNAs. To better understand the molecular features of diseases caused by these mutations, and to improve their diagnosis and therapeutics, we have constructed a Drosophila melanogaster model disrupting the mitochondrial seryl-tRNA synthetase by RNA interference. At the molecular level, the knockdown generates a reduction in transfer RNA serylation, which correlates with the severity of the phenotype observed. The silencing compromises viability, longevity, motility and tissue development. At the cellular level, the knockdown alters mitochondrial morphology, biogenesis and function, and induces lactic acidosis and reactive oxygen species accumulation. We report that administration of antioxidant compounds has a palliative effect of some of these phenotypes. In conclusion, the fly model generated in this work reproduces typical characteristics of pathologies caused by mutations in the mitochondrial aminoacylation system, and can be useful to assess therapeutic approaches.

The chromosomal proteins JIL-1 and Z4/Putzig regulate the telomeric chromatin in Drosophila melanogaster.

Drosophila telomere maintenance depends on the transposition of the specialized retrotransposons HeT-A, TART, and TAHRE. Controlling the activation and silencing of these elements is crucial for a precise telomere function without compromising genomic integrity. Here we describe two chromosomal proteins, JIL-1 and Z4 (also known as Putzig), which are necessary for establishing a fine-tuned regulation of the transcription of the major component of Drosophila telomeres, the HeT-A retrotransposon, thus guaranteeing genome stability. We found that mutant alleles of JIL-1 have decreased HeT-A transcription, putting forward this kinase as the first positive regulator of telomere transcription in Drosophila described to date. We describe how the decrease in HeT-A transcription in JIL-1 alleles correlates with an increase in silencing chromatin marks such as H3K9me3 and HP1a at the HeT-A promoter. Moreover, we have detected that Z4 mutant alleles show moderate telomere instability, suggesting an important role of the JIL-1-Z4 complex in establishing and maintaining an appropriate chromatin environment at Drosophila telomeres. Interestingly, we have detected a biochemical interaction between Z4 and the HeT-A Gag protein, which could explain how the Z4-JIL-1 complex is targeted to the telomeres. Accordingly, we demonstrate that a phenotype of telomere instability similar to that observed for Z4 mutant alleles is found when the gene that encodes the HeT-A Gag protein is knocked down. We propose a model to explain the observed transcriptional and stability changes in relation to other heterochromatin components characteristic of Drosophila telomeres, such as HP1a.

GAGA factor repression of transcription is a rare event but the negative regulation of Trl is conserved in Drosophila species.

GAGA is a highly conserved Drosophila transcription factor encoded by the Trithorax-like (Trl) gene. While GAGA usually activates transcription, it represses its own promoter. Here we show that GAGA-mediated repression of Trl is conserved between two distant Drosophila species. A detailed promoter study showed that GAGA repressive activity can't be attributed to any discrete element in the Trl promoter. Genome-wide analysis of the transcriptome in S2 cells indicated that repression of Trl is very likely unique, being GAGA factor a transactivator for all the other promoters. Taken together, our results suggest a new mechanism to explain GAGA-mediated repression that involves a dose-dependent change in the architecture of the Trl promoter.

Novel Targeting of DNA Methyltransferase Activity Inhibits Ewing Sarcoma Cell Proliferation and Enhances Tumor Cell Sensitivity to DNA Damaging Drugs by Activating the DNA Damage Response.

DNA methylation is an important component of the epigenetic machinery that regulates the malignancy of Ewing sarcoma (EWS), the second most common primary bone tumor in children and adolescents. Coordination of DNA methylation and DNA replication is critical for maintaining epigenetic programming and the DNMT1 enzyme has been demonstrated to have an important role in both maintaining the epigenome and controlling cell cycle. Here, we showed that the novel nonnucleoside DNMT inhibitor (DNMTi) MC3343 induces a specific depletion of DNMT1 and affects EWS tumor proliferation through a mechanism that is independent on DNA methylation. Depletion of DNMT1 causes perturbation of the cell cycle, with an accumulation of cells in the G1 phase, and DNA damage, as revealed by the induction of γH2AX foci. These effects elicited activation of p53-dependent signaling and apoptosis in p53wt cells, while in p53 mutated cells, persistent micronuclei and increased DNA instability was observed. Treatment with MC3343 potentiates the efficacy of DNA damaging agents such as doxorubicin and PARP-inhibitors (PARPi). This effect correlates with increased DNA damage and synergistic tumor cytotoxicity, supporting the use of the DNMTi MC3343 as an adjuvant agent in treating EWS.

Look-alike humans identified by facial recognition algorithms show genetic similarities.

The human face is one of the most visible features of our unique identity as individuals. Interestingly, monozygotic twins share almost identical facial traits and the same DNA sequence but could exhibit differences in other biometrical parameters. The expansion of the world wide web and the possibility to exchange pictures of humans across the planet has increased the number of people identified online as virtual twins or doubles that are not family related. Herein, we have characterized in detail a set of "look-alike" humans, defined by facial recognition algorithms, for their multiomics landscape. We report that these individuals share similar genotypes and differ in their DNA methylation and microbiome landscape. These results not only provide insights about the genetics that determine our face but also might have implications for the establishment of other human anthropometric properties and even personality characteristics.

Genetically predicted telomere length and Alzheimer's disease endophenotypes: a Mendelian randomization study.

Telomere length (TL) is associated with biological aging, consequently influencing the risk of age-related diseases such as Alzheimer's disease (AD). We aimed to evaluate the potential causal role of TL in AD endophenotypes (i.e., cognitive performance, N = 2233; brain age and AD-related signatures, N = 1134; and cerebrospinal fluid biomarkers (CSF) of AD and neurodegeneration, N = 304) through a Mendelian randomization (MR) analysis. Our analysis was conducted in the context of the ALFA (ALzheimer and FAmilies) study, a population of cognitively healthy individuals at risk of AD. A total of 20 single nucleotide polymorphisms associated with TL were used to determine the effect of TL on AD endophenotypes. Analyses were adjusted by age, sex, and years of education. Stratified analyses by APOE-ɛ4 status and polygenic risk score of AD were conducted. MR analysis revealed significant associations between genetically predicted longer TL and lower levels of CSF Aß and higher levels of CSF NfL only in APOE-ɛ4 non-carriers. Moreover, inheriting longer TL was associated with greater cortical thickness in age and AD-related brain signatures and lower levels of CSF p-tau among individuals at a high genetic predisposition to AD. Further observational analyses are warranted to better understand these associations.

Transcriptional analysis of the HeT-A retrotransposon in mutant and wild type stocks reveals high sequence variability at Drosophila telomeres and other unusual features.

BACKGROUND: Telomere replication in Drosophila depends on the transposition of a domesticated retroelement, the HeT-A retrotransposon. The sequence of the HeT-A retrotransposon changes rapidly resulting in differentiated subfamilies. This pattern of sequence change contrasts with the essential function with which the HeT-A is entrusted and brings about questions concerning the extent of sequence variability, the telomere contribution of different subfamilies, and whether wild type and mutant Drosophila stocks show different HeT-A scenarios. RESULTS: A detailed study on the variability of HeT-A reveals that both the level of variability and the number of subfamilies are higher than previously reported. Comparisons between GIII, a strain with longer telomeres, and its parental strain Oregon-R indicate that both strains have the same set of HeT-A subfamilies. Finally, the presence of a highly conserved splicing pattern only in its antisense transcripts indicates a putative regulatory, functional or structural role for the HeT-A RNA. Interestingly, our results also suggest that most HeT-A copies are actively expressed regardless of which telomere and where in the telomere they are located. CONCLUSIONS: Our study demonstrates how the HeT-A sequence changes much faster than previously reported resulting in at least nine different subfamilies most of which could actively contribute to telomere extension in Drosophila. Interestingly, the only significant difference observed between Oregon-R and GIII resides in the nature and proportion of the antisense transcripts, suggesting a possible mechanism that would in part explain the longer telomeres of the GIII stock.

Genetic Predisposition to Alzheimer's Disease Is Associated with Enlargement of Perivascular Spaces in Centrum Semiovale Region.

This study investigated whether genetic factors involved in Alzheimer's disease (AD) are associated with enlargement of Perivascular Spaces (ePVS) in the brain. A total of 680 participants with T2-weighted MRI scans and genetic information were acquired from the ALFA study. ePVS in the basal ganglia (BG) and the centrum semiovale (CS) were assessed based on a validated visual rating scale. We used univariate and multivariate logistic regression models to investigate associations between ePVS in BG and CS with BIN1-rs744373, as well as APOE genotypes. We found a significant association of the BIN1-rs744373 polymorphism in the CS subscale (p value = 0.019; OR = 2.564), suggesting that G allele carriers have an increased risk of ePVS in comparison with A allele carriers. In stratified analysis by APOE-ε4 status (carriers vs. non-carriers), these results remained significant only for ε4 carriers (p value = 0.011; OR = 1.429). To our knowledge, the present study is the first suggesting that genetic predisposition for AD is associated with ePVS in CS. These findings provide evidence that underlying biological processes affecting AD may influence CS-ePVS.

Gene Amplification-Associated Overexpression of the Selenoprotein tRNA Enzyme TRIT1 Confers Sensitivity to Arsenic Trioxide in Small-Cell Lung Cancer.

The alteration of RNA modification patterns is emerging as a common feature of human malignancies. If these changes affect key RNA molecules for mRNA translation, such as transfer RNA, they can have important consequences for cell transformation. TRIT1 is the enzyme responsible for the hypermodification of adenosine 37 in the anticodon region of human tRNAs containing serine and selenocysteine. Herein, we show that TRIT1 undergoes gene amplification-associated overexpression in cancer cell lines and primary samples of small-cell lung cancer. From growth and functional standpoints, the induced depletion of TRIT1 expression in amplified cells reduces their tumorigenic potential and downregulates the selenoprotein transcripts. We observed that TRIT1-amplified cells are sensitive to arsenic trioxide, a compound that regulates selenoproteins, whereas reduction of TRIT1 levels confers loss of sensitivity to the drug. Overall, our results indicate a role for TRIT1 as a small-cell lung cancer-relevant gene that, when undergoing gene amplification-associated activation, can be targeted with the differentiation agent arsenic trioxide.

Molecular Features of Cancers Exhibiting Exceptional Responses to Treatment.

A small fraction of cancer patients with advanced disease survive significantly longer than patients with clinically comparable tumors. Molecular mechanisms for exceptional responses to therapy have been identified by genomic analysis of tumor biopsies from individual patients. Here, we analyzed tumor biopsies from an unbiased cohort of 111 exceptional responder patients using multiple platforms to profile genetic and epigenetic aberrations as well as the tumor microenvironment. Integrative analysis uncovered plausible mechanisms for the therapeutic response in nearly a quarter of the patients. The mechanisms were assigned to four broad categories-DNA damage response, intracellular signaling, immune engagement, and genetic alterations characteristic of favorable prognosis-with many tumors falling into multiple categories. These analyses revealed synthetic lethal relationships that may be exploited therapeutically and rare genetic lesions that favor therapeutic success, while also providing a wealth of testable hypotheses regarding oncogenic mechanisms that may influence the response to cancer therapy.

The DNA methylation landscape of human cancer organoids available at the American type culture collection.

One caveat in cancer research is the dependence of certain experimental systems that might not really reflect the properties of the primary tumours. The recent irruption of 3D cultured cells termed organoids could render a better representation of the original tumour sample. However, every laboratory has its own protocol and tissue-provider to establish these cancer models, preventing further dissemination and validation of the obtained data. To address this problem, the Human Cancer Models Initiative (HCMI) has selected the American Type Culture Collection (ATCC) to make available organoid models to the scientific community. In this regard, no epigenetic information is available for these samples and, overall, the DNA methylation profiles of human cancer organoids are largely unknown. Herein, we provide the DNA methylation landscape of 25 human cancer organoids available at the ATCC using a microarray that interrogates more than 850,000 CpG sites. We observed that the studied organoids retain the epigenetic setting of their original primary cancer type; that exhibit a DNA methylation landscape characteristic of transformed tissues excluding an overgrowth of normal-matched cells; and that are closer to the DNA methylation profiles of the corresponding primary tumours than to established 2D cell lines. Most importantly, the obtained DNA methylation results are freely available to everyone for further data mining. Thus, our findings support from the epigenetic standpoint that the ATCC human cancer organoids recapitulate many of the features of the disorder in the patient and are excellent tools to be shared among investigators for further tumour biology research.

General, negative feedback mechanism for regulation of Trithorax-like gene expression in vivo: new roles for GAGA factor in flies.

Expression of every gene is first regulated at the transcriptional level. While some genes show acute and discrete periods of expression others show a rather steady expression level throughout development. An example of the latter is Trithorax-like (Trl) a member of the Trithorax group that encodes GAGA factor in Drosophila. Among other functions, GAGA factor has been described to stimulate transcription of several genes, including some homeotic genes. Here we show that GAGA factor is continuously down-regulating the expression of its own promoter using a negative feedback mechanism in vivo. Like its expression, repression by GAGA factor is ubiquitous, prevents its accumulation, and takes place throughout development. Experimental alteration of GAGA factor dosage results in several unexpected phenotypes, not related to alteration of homeotic gene expression, but rather to functions that take place later during development and affect different morphogenetic processes. The results suggest that GAGA factor is essential during development, even after homeotic gene expression is established, and indicate the existence of an upper limit for GAGA factor dosage that should not be exceeded.