Designing realistic regulatory DNA with autoregressive language models.

Cis-regulatory elements (CREs), such as promoters and enhancers, are DNA sequences that regulate the expression of genes. The activity of a CRE is influenced by the order, composition, and spacing of sequence motifs that are bound by proteins called transcription factors (TFs). Synthetic CREs with specific properties are needed for biomanufacturing as well as for many therapeutic applications including cell and gene therapy. Here, we present regLM, a framework to design synthetic CREs with desired properties, such as high, low, or cell type-specific activity, using autoregressive language models in conjunction with supervised sequence-to-function models. We used our framework to design synthetic yeast promoters and cell type-specific human enhancers. We demonstrate that the synthetic CREs generated by our approach are not only predicted to have the desired functionality but also contain biological features similar to experimentally validated CREs. regLM thus facilitates the design of realistic regulatory DNA elements while providing insights into the cis-regulatory code.

IL-4-induced SOX9 confers lineage plasticity to aged adult lung stem cells.

Wound healing in response to acute injury is mediated by the coordinated and transient activation of parenchymal, stromal, and immune cells that resolves to homeostasis. Environmental, genetic, and epigenetic factors associated with inflammation and aging can lead to persistent activation of the microenvironment and fibrosis. Here, we identify opposing roles of interleukin-4 (IL-4) cytokine signaling in interstitial macrophages and type II alveolar epithelial cells (ATIIs). We show that IL4Ra signaling in macrophages promotes regeneration of the alveolar epithelium after bleomycin-induced lung injury. Using organoids and mouse models, we show that IL-4 directly acts on a subset of ATIIs to induce the expression of the transcription factor SOX9 and reprograms them toward a progenitor-like state with both airway and alveolar lineage potential. In the contexts of aging and bleomycin-induced lung injury, this leads to aberrant epithelial cell differentiation and bronchiolization, consistent with cellular and histological changes observed in interstitial lung disease.

Obligate role for Rock1 and Rock2 in adult stem cell viability and function.

The ability of stem cells to rapidly proliferate and differentiate is integral to the steady-state maintenance of tissues with high turnover such as the blood and intestine. Mutations that alter these processes can cause primary immunodeficiencies, malignancies and defects in barrier function. The Rho-kinases, Rock1 and Rock2, regulate cell shape and cytoskeletal rearrangement, activities essential to mitosis. Here, we use inducible gene targeting to ablate Rock1 and Rock2 in adult mice, and identify an obligate requirement for these enzymes in the preservation of the hematopoietic and gastrointestinal systems. Hematopoietic cell progenitors devoid of Rho-kinases display cell cycle arrest, blocking the differentiation to mature blood lineages. Similarly, these mice exhibit impaired epithelial cell renewal in the small intestine, which is ultimately fatal. Our data reveal a novel role for these kinases in the proliferation and viability of stem cells and their progenitors, which is vital to maintaining the steady-state integrity of these organ systems.

Senescent cells perturb intestinal stem cell differentiation through Ptk7 induced noncanonical Wnt and YAP signaling.

Cellular senescence and the senescence-associated secretory phenotype (SASP) are implicated in aging and age-related disease, and SASP-related inflammation is thought to contribute to tissue dysfunction in aging and diseased animals. However, whether and how SASP factors influence the regenerative capacity of tissues remains unclear. Here, using intestinal organoids as a model of tissue regeneration, we show that SASP factors released by senescent fibroblasts deregulate stem cell activity and differentiation and ultimately impair crypt formation. We identify the secreted N-terminal domain of Ptk7 as a key component of the SASP that activates non-canonical Wnt / Ca2+ signaling through FZD7 in intestinal stem cells (ISCs). Changes in cytosolic [Ca2+] elicited by Ptk7 promote nuclear translocation of YAP and induce expression of YAP/TEAD target genes, impairing symmetry breaking and stem cell differentiation. Our study discovers secreted Ptk7 as a factor released by senescent cells and provides insight into the mechanism by which cellular senescence contributes to tissue dysfunction in aging and disease.

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos.

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology.

An otopetrin family proton channel promotes cellular acid efflux critical for biomineralization in a marine calcifier.

Otopetrins comprise a family of proton-selective channels that are critically important for the mineralization of otoliths and statoconia in vertebrates but whose underlying cellular mechanisms remain largely unknown. Here, we demonstrate that otopetrins are critically involved in the calcification process by providing an exit route for protons liberated by the formation of CaCO3 Using the sea urchin larva, we examined the otopetrin ortholog otop2l, which is exclusively expressed in the calcifying primary mesenchymal cells (PMCs) that generate the calcitic larval skeleton. otop2l expression is stimulated during skeletogenesis, and knockdown of otop2l impairs spicule formation. Intracellular pH measurements demonstrated Zn2+-sensitive H+ fluxes in PMCs that regulate intracellular pH in a Na+/HCO3--independent manner, while Otop2l knockdown reduced membrane proton permeability. Furthermore, Otop2l displays unique features, including strong activation by high extracellular pH (>8.0) and check-valve-like outwardly rectifying H+ flux properties, making it into a cellular proton extrusion machine adapted to oceanic living conditions. Our results provide evidence that otopetrin family proton channels are a central component of the cellular pH regulatory machinery in biomineralizing cells. Their ubiquitous occurrence in calcifying systems across the animal kingdom suggest a conserved physiological function by mediating pH at the site of mineralization. This important role of otopetrin family proton channels has strong implications for our view on the cellular mechanisms of biomineralization and their response to changes in oceanic pH.

Cis-acting variation is common across regulatory layers but is often buffered during embryonic development.

Precise patterns of gene expression are driven by interactions between transcription factors, regulatory DNA sequences, and chromatin. How DNA mutations affecting any one of these regulatory "layers" are buffered or propagated to gene expression remains unclear. To address this, we quantified allele-specific changes in chromatin accessibility, histone modifications, and gene expression in F1 embryos generated from eight Drosophila crosses at three embryonic stages, yielding a comprehensive data set of 240 samples spanning multiple regulatory layers. Genetic variation (allelic imbalance) impacts gene expression more frequently than chromatin features, with metabolic and environmental response genes being most often affected. Allelic imbalance in cis-regulatory elements (enhancers) is common and highly heritable, yet its functional impact does not generally propagate to gene expression. When it does, genetic variation impacts RNA levels through two alternative mechanisms involving either H3K4me3 or chromatin accessibility and H3K27ac. Changes in RNA are more predictive of variation in H3K4me3 than vice versa, suggesting a role for H3K4me3 downstream from transcription. The impact of a substantial proportion of genetic variation is consistent across embryonic stages, with 50% of allelic imbalanced features at one stage being also imbalanced at subsequent developmental stages. Crucially, buffering, as well as the magnitude and evolutionary impact of genetic variants, is influenced by regulatory complexity (i.e., number of enhancers regulating a gene), with transcription factors being most robust to cis-acting, but most influenced by trans-acting, variation.

Lineage-Resolved Enhancer and Promoter Usage during a Time Course of Embryogenesis.

Enhancers are essential drivers of cell states, yet the relationship between accessibility, regulatory activity, and in vivo lineage commitment during embryogenesis remains poorly understood. Here, we measure chromatin accessibility in isolated neural and mesodermal lineages across a time course of Drosophila embryogenesis. Promoters, including tissue-specific genes, are often constitutively open, even in contexts where the gene is not expressed. In contrast, the majority of distal elements have dynamic, tissue-specific accessibility. Enhancer priming appears rarely within a lineage, perhaps reflecting the speed of Drosophila embryogenesis. However, many tissue-specific enhancers are accessible in other lineages early on and become progressively closed as embryogenesis proceeds. We demonstrate the usefulness of this tissue- and time-resolved resource to definitively identify single-cell clusters, to uncover predictive motifs, and to identify many regulators of tissue development. For one such predicted neural regulator, l(3)neo38, we generate a loss-of-function mutant and uncover an essential role for neuromuscular junction and brain development.

The Role of Chromatin Accessibility in cis-Regulatory Evolution.

Transcription factor (TF) binding is determined by sequence as well as chromatin accessibility. Although the role of accessibility in shaping TF-binding landscapes is well recorded, its role in evolutionary divergence of TF binding, which in turn can alter cis-regulatory activities, is not well understood. In this work, we studied the evolution of genome-wide binding landscapes of five major TFs in the core network of mesoderm specification, between Drosophila melanogaster and Drosophila virilis, and examined its relationship to accessibility and sequence-level changes. We generated chromatin accessibility data from three important stages of embryogenesis in both Drosophila melanogaster and Drosophila virilis and recorded conservation and divergence patterns. We then used multivariable models to correlate accessibility and sequence changes to TF-binding divergence. We found that accessibility changes can in some cases, for example, for the master regulator Twist and for earlier developmental stages, more accurately predict binding change than is possible using TF-binding motif changes between orthologous enhancers. Accessibility changes also explain a significant portion of the codivergence of TF pairs. We noted that accessibility and motif changes offer complementary views of the evolution of TF binding and developed a combined model that captures the evolutionary data much more accurately than either view alone. Finally, we trained machine learning models to predict enhancer activity from TF binding and used these functional models to argue that motif and accessibility-based predictors of TF-binding change can substitute for experimentally measured binding change, for the purpose of predicting evolutionary changes in enhancer activity.

Distribution and Clinical Utility of the 21-gene Recurrence Score in Pure Mucinous Breast Cancer Patients: a case-control study.

The 21-gene recurrence score (RS) is increasingly being used for patients with early stage, hormone receptor-positive, Her-2-negative breast cancer. However, these results are largely from populations with infiltrating ductal carcinoma (IDC). The clinical value of RS testing in mucinous carcinoma has not been well investigated. Pure mucinous breast cancer (PMBC) and paired pure IDC patients who underwent 21-gene RS were retrospectively reviewed and matched with tumor stage and molecular subtype. Clinic-pathological factors, treatment strategies, and RS distribution were compared between the PMBC and IDC patients. A total of 35 PMBC and 70 IDC patients were included. We found that RS was lower in the PMBC as compared with the IDC group: 21.26 vs. 24.40 (P=0.037). Regarding RS categories, PMBC patients had a relatively lower percentage of high RS patients than the IDC group: 8.57% vs. 22.86% (P = 0.048). Multivariate analysis showed that histologic type was an independent factor predicting RS distribution: IDC patients were associated with a higher RS as compared with PMBC patients (OR: 1.27, 95% CI: 1.03-2.13; P=0.014). Among genes in 21-gene RS testing, HER2, STMY3, STK15, and BAG1 were significantly different between the PMBC and IDC groups (P < 0.05). Two patients (5.71%) in the PMBC group, both with high RS, were recommended to receive adjuvant chemotherapy, much lower than patients with IDC (57.14%, P < 0.001). In multivariate analysis, histologic type of IDC was an independent factor for chemotherapy recommendation (OR = 22.00, 95% CI: 4.89-98.97, P<0.001). With a medium follow-up time of 24 months, one IDC patient had ipsilateral axillary lymph nodes recurrence and one PMBC patient had contralateral breast cancer. In conclusion, PMBC patients, mostly classified with low or intermediate RS category, were associated with lower RS as compared with IDC patients. PMBC and IDC had different genes expression patterns. Patients with high RS in the PMBC group might be recommended to receive adjuvant chemotherapy, which deserves further clinical evaluation.

Selective Filopodia Adhesion Ensures Robust Cell Matching in the Drosophila Heart.

The ability to form specific cell-cell connections within complex cellular environments is critical for multicellular organisms. However, the underlying mechanisms of cell matching that instruct these connections remain elusive. Here, we quantitatively explored the dynamics and regulation of cell matching processes utilizing Drosophila cardiogenesis. We found that cell matching is highly robust at boundaries between cardioblast (CB) subtypes, and filopodia of different CB subtypes have distinct binding affinities. Cdc42 is involved in regulating this selective filopodia binding adhesion and influences CB matching. Further, we identified adhesion molecules Fasciclin III (Fas3) and Ten-m, both of which also regulate synaptic targeting, as having complementary differential expression in CBs. Altering Fas3 expression changes differential filopodia adhesion and leads to CB mismatch. Furthermore, only when both Fas3 and Ten-m are lost is CB alignment severely impaired. Our results show that differential adhesion mediated by selective filopodia binding efficiently regulates precise and robust cell matching.

The cis-regulatory dynamics of embryonic development at single-cell resolution.

Understanding how gene regulatory networks control the progressive restriction of cell fates is a long-standing challenge. Recent advances in measuring gene expression in single cells are providing new insights into lineage commitment. However, the regulatory events underlying these changes remain unclear. Here we investigate the dynamics of chromatin regulatory landscapes during embryogenesis at single-cell resolution. Using single-cell combinatorial indexing assay for transposase accessible chromatin with sequencing (sci-ATAC-seq), we profiled chromatin accessibility in over 20,000 single nuclei from fixed Drosophila melanogaster embryos spanning three landmark embryonic stages: 2-4 h after egg laying (predominantly stage 5 blastoderm nuclei), when each embryo comprises around 6,000 multipotent cells; 6-8 h after egg laying (predominantly stage 10-11), to capture a midpoint in embryonic development when major lineages in the mesoderm and ectoderm are specified; and 10-12 h after egg laying (predominantly stage 13), when each of the embryo's more than 20,000 cells are undergoing terminal differentiation. Our results show that there is spatial heterogeneity in the accessibility of the regulatory genome before gastrulation, a feature that aligns with future cell fate, and that nuclei can be temporally ordered along developmental trajectories. During mid-embryogenesis, tissue granularity emerges such that individual cell types can be inferred by their chromatin accessibility while maintaining a signature of their germ layer of origin. Analysis of the data reveals overlapping usage of regulatory elements between cells of the endoderm and non-myogenic mesoderm, suggesting a common developmental program that is reminiscent of the mesendoderm lineage in other species. We identify 30,075 distal regulatory elements that exhibit tissue-specific accessibility. We validated the germ-layer specificity of a subset of these predicted enhancers in transgenic embryos, achieving an accuracy of 90%. Overall, our results demonstrate the power of shotgun single-cell profiling of embryos to resolve dynamic changes in the chromatin landscape during development, and to uncover the cis-regulatory programs of metazoan germ layers and cell types.

Upconverting nanoparticle micro-lightbulbs designed for deep tissue optical stimulation and imaging.

Optical methods for imaging and stimulation of biological events based on the use of visible light are limited to the superficial layers of tissue due to the significant absorption and scattering of light. Here, we demonstrate the design and implementation of passive micro-structured lightbulbs (MLBs) containing bright-emitting lanthanide-doped upconverting nanoparticles (UCNPs) for light delivery deep into the tissue. The MLBs are realized as cylindrical pillars made of Parylene C polymer that can be implanted deep into the tissue. The encapsulated UCNPs absorb near-infrared (NIR) light at λ = 980 nm, which undergoes much less absorption than the blue light in the brain tissue, and then locally emit blue light (1G4→3H6 and 1D2→3F4 transitions) that can be used for optogenetic excitation of neurons in the brain. The 3H4→3H6 transition will result in the emission of higher energy NIR photons at λ = 800 nm that can be used for imaging and tracking MLBs through thick tissue.

Promoter shape varies across populations and affects promoter evolution and expression noise.

Animal promoters initiate transcription either at precise positions (narrow promoters) or dispersed regions (broad promoters), a distinction referred to as promoter shape. Although highly conserved, the functional properties of promoters with different shapes and the genetic basis of their evolution remain unclear. Here we used natural genetic variation across a panel of 81 Drosophila lines to measure changes in transcriptional start site (TSS) usage, identifying thousands of genetic variants affecting transcript levels (strength) or the distribution of TSSs within a promoter (shape). Our results identify promoter shape as a molecular trait that can evolve independently of promoter strength. Broad promoters typically harbor shape-associated variants, with signatures of adaptive selection. Single-cell measurements demonstrate that variants modulating promoter shape often increase expression noise, whereas heteroallelic interactions with other promoter variants alleviate these effects. These results uncover new functional properties of natural promoters and suggest the minimization of expression noise as an important factor in promoter evolution.

Can Clinically Node-Negative Breast Cancer Patients with Suspicious Axillary Lymph Nodes at Ultrasound But Negative Fine-Needle Aspiration Be Approached as Having Node-Negative Disease?

BACKGROUND: Few data exist to elucidate whether patients with a suspicious axillary lymph node (ALN) at ultrasound but a negative fine-needle aspiration result (FNA group) can be managed as having ultrasound node-negative disease (AUN group). This study compared various ALN statuses between the AUN and FNA groups to guide further ALN management. METHODS: Patients who had clinical T1-2N0 breast cancer treated with sentinel lymph node (SLN) biopsy were retrospectively analyzed. The ALN metastasis status, SLN status, and non-SLN metastasis rates in the entire population and the patients with one or two positive SLNs were compared between the AUN and FNA groups. RESULTS: A total of 1007 patients (886 AUN and 121 FNA patients) were eligible for the final analysis: The incidence of ALN metastasis did not differ between the AUN group (16.5%) and the FNA group (21.5%) (P = 0.170). In addition, three or more metastases were found in only 2.4% of the AUN patients and 3.3% of the FNA patients (P = 0.405). The non-SLN metastasis rate was 22.6% (33/146) in the AUN group and 19.2% (5/26) in the FNA group (P = 0.699). For the patients with one or two positive SLNs, the rate of non-SLN metastasis was similar between the AUN group (19.6%, 27/138) and the FNA group (12.5%, 3/24) (P = 0.591). CONCLUSIONS: Patients with a suspicious ALN at ultrasound but a negative FNA result had ALN statuses similar to those of the ultrasound node-negative patients, indicating that these patients can be treated as having ultrasound node-negative disease.

Genetic variants regulating expression levels and isoform diversity during embryogenesis.

Embryonic development is driven by tightly regulated patterns of gene expression, despite extensive genetic variation among individuals. Studies of expression quantitative trait loci (eQTL) indicate that genetic variation frequently alters gene expression in cell-culture models and differentiated tissues. However, the extent and types of genetic variation impacting embryonic gene expression, and their interactions with developmental programs, remain largely unknown. Here we assessed the effect of genetic variation on transcriptional (expression levels) and post-transcriptional (3' RNA processing) regulation across multiple stages of metazoan development, using 80 inbred Drosophila wild isolates, identifying thousands of developmental-stage-specific and shared QTL. Given the small blocks of linkage disequilibrium in Drosophila, we obtain near base-pair resolution, resolving causal mutations in developmental enhancers, validated transcription-factor-binding sites and RNA motifs. This fine-grain mapping uncovered extensive allelic interactions within enhancers that have opposite effects, thereby buffering their impact on enhancer activity. QTL affecting 3' RNA processing identify new functional motifs leading to transcript isoform diversity and changes in the lengths of 3' untranslated regions. These results highlight how developmental stage influences the effects of genetic variation and uncover multiple mechanisms that regulate and buffer expression variation during embryogenesis.

Energy-Looping Nanoparticles: Harnessing Excited-State Absorption for Deep-Tissue Imaging.

Near infrared (NIR) microscopy enables noninvasive imaging in tissue, particularly in the NIR-II spectral range (1000-1400 nm) where attenuation due to tissue scattering and absorption is minimized. Lanthanide-doped upconverting nanocrystals are promising deep-tissue imaging probes due to their photostable emission in the visible and NIR, but these materials are not efficiently excited at NIR-II wavelengths due to the dearth of lanthanide ground-state absorption transitions in this window. Here, we develop a class of lanthanide-doped imaging probes that harness an energy-looping mechanism that facilitates excitation at NIR-II wavelengths, such as 1064 nm, that are resonant with excited-state absorption transitions but not ground-state absorption. Using computational methods and combinatorial screening, we have identified Tm(3+)-doped NaYF4 nanoparticles as efficient looping systems that emit at 800 nm under continuous-wave excitation at 1064 nm. Using this benign excitation with standard confocal microscopy, energy-looping nanoparticles (ELNPs) are imaged in cultured mammalian cells and through brain tissue without autofluorescence. The 1 mm imaging depths and 2 µm feature sizes are comparable to those demonstrated by state-of-the-art multiphoton techniques, illustrating that ELNPs are a promising class of NIR probes for high-fidelity visualization in cells and tissue.

Next-generation sequencing-based detection of germline L1-mediated transductions.

BACKGROUND: While active LINE-1 (L1) elements possess the ability to mobilize flanking sequences to different genomic loci through a process termed transduction influencing genomic content and structure, an approach for detecting polymorphic germline non-reference transductions in massively-parallel sequencing data has been lacking. RESULTS: Here we present the computational approach TIGER (Transduction Inference in GERmline genomes), enabling the discovery of non-reference L1-mediated transductions by combining L1 discovery with detection of unique insertion sequences and detailed characterization of insertion sites. We employed TIGER to characterize polymorphic transductions in fifteen genomes from non-human primate species (chimpanzee, orangutan and rhesus macaque), as well as in a human genome. We achieved high accuracy as confirmed by PCR and two single molecule DNA sequencing techniques, and uncovered differences in relative rates of transduction between primate species. CONCLUSIONS: By enabling detection of polymorphic transductions, TIGER makes this form of relevant structural variation amenable for population and personal genome analysis.

Genomic Characterization of the Evolutionary Potential of the Sea Urchin Strongylocentrotus droebachiensis Facing Ocean Acidification.

Ocean acidification (OA) is increasing due to anthropogenic CO2 emissions and poses a threat to marine species and communities worldwide. To better project the effects of acidification on organisms' health and persistence, an understanding is needed of the 1) mechanisms underlying developmental and physiological tolerance and 2) potential populations have for rapid evolutionary adaptation. This is especially challenging in nonmodel species where targeted assays of metabolism and stress physiology may not be available or economical for large-scale assessments of genetic constraints. We used mRNA sequencing and a quantitative genetics breeding design to study mechanisms underlying genetic variability and tolerance to decreased seawater pH (-0.4 pH units) in larvae of the sea urchin Strongylocentrotus droebachiensis. We used a gene ontology-based approach to integrate expression profiles into indirect measures of cellular and biochemical traits underlying variation in larval performance (i.e., growth rates). Molecular responses to OA were complex, involving changes to several functions such as growth rates, cell division, metabolism, and immune activities. Surprisingly, the magnitude of pH effects on molecular traits tended to be small relative to variation attributable to segregating functional genetic variation in this species. We discuss how the application of transcriptomics and quantitative genetics approaches across diverse species can enrich our understanding of the biological impacts of climate change.

Shadow Enhancers Are Pervasive Features of Developmental Regulatory Networks.

Embryogenesis is remarkably robust to segregating mutations and environmental variation; under a range of conditions, embryos of a given species develop into stereotypically patterned organisms. Such robustness is thought to be conferred, in part, through elements within regulatory networks that perform similar, redundant tasks. Redundant enhancers (or "shadow" enhancers), for example, can confer precision and robustness to gene expression, at least at individual, well-studied loci. However, the extent to which enhancer redundancy exists and can thereby have a major impact on developmental robustness remains unknown. Here, we systematically assessed this, identifying over 1,000 predicted shadow enhancers during Drosophila mesoderm development. The activity of 23 elements, associated with five genes, was examined in transgenic embryos, while natural structural variation among individuals was used to assess their ability to buffer against genetic variation. Our results reveal three clear properties of enhancer redundancy within developmental systems. First, it is much more pervasive than previously anticipated, with 64% of loci examined having shadow enhancers. Their spatial redundancy is often partial in nature, while the non-overlapping function may explain why these enhancers are maintained within a population. Second, over 70% of loci do not follow the simple situation of having only two shadow enhancers-often there are three (rols), four (CadN and ade5), or five (Traf1), at least one of which can be deleted with no obvious phenotypic effects. Third, although shadow enhancers can buffer variation, patterns of segregating variation suggest that they play a more complex role in development than generally considered.