Low affinity binding site clusters confer hox specificity and regulatory robustness.

In animals, Hox transcription factors define regional identity in distinct anatomical domains. How Hox genes encode this specificity is a paradox, because different Hox proteins bind with high affinity in vitro to similar DNA sequences. Here, we demonstrate that the Hox protein Ultrabithorax (Ubx) in complex with its cofactor Extradenticle (Exd) bound specifically to clusters of very low affinity sites in enhancers of the shavenbaby gene of Drosophila. These low affinity sites conferred specificity for Ubx binding in vivo, but multiple clustered sites were required for robust expression when embryos developed in variable environments. Although most individual Ubx binding sites are not evolutionarily conserved, the overall enhancer architecture-clusters of low affinity binding sites-is maintained and required for enhancer function. Natural selection therefore works at the level of the enhancer, requiring a particular density of low affinity Ubx sites to confer both specific and robust expression.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila.

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

Phylogenomics provides insights into the evolution of cactophily and host plant shifts in Drosophila.

Cactophilic species of the Drosophila buzzatii cluster (repleta group) comprise an excellent model group to investigate genomic changes underlying adaptation to extreme climate conditions and host plants. In particular, these species form a tractable system to study the transition from chemically simpler breeding sites (like prickly pears of the genus Opuntia) to chemically more complex hosts (columnar cacti). Here, we report four highly contiguous genome assemblies of three species of the buzzatii cluster. Based on this genomic data and inferred phylogenetic relationships, we identified candidate taxonomically restricted genes (TRGs) likely involved in the evolution of cactophily and cactus host specialization. Functional enrichment analyses of TRGs within the buzzatii cluster identified genes involved in detoxification, water preservation, immune system response, anatomical structure development, and morphogenesis. In contrast, processes that regulate responses to stress, as well as the metabolism of nitrogen compounds, transport, and secretion were found in the set of species that are columnar cacti dwellers. These findings are in line with the hypothesis that those genomic changes brought about key mechanisms underlying the adaptation of the buzzatii cluster species to arid regions in South America.

Actors with Multiple Roles: Pleiotropic Enhancers and the Paradigm of Enhancer Modularity.

The current paradigm in the field of gene regulation postulates that regulatory information for generating gene expression is organized into modules (enhancers), each containing the information for driving gene expression in a single spatiotemporal context. This modular organization is thought to facilitate the evolution of gene expression by minimizing pleiotropic effects. Here we review recent studies that provide evidence of quite the opposite: (i) enhancers can function in multiple developmental contexts, implying that enhancers can be pleiotropic, (ii) transcription factor binding sites within pleiotropic enhancers are reused in different contexts, and (iii) pleiotropy impacts the structure and evolution of enhancers. Altogether, this evidence suggests that enhancer pleiotropy is pervasive in animal genomes, challenging the commonly held view of modularity.

Gene regulatory network architecture in different developmental contexts influences the genetic basis of morphological evolution.

Convergent phenotypic evolution is often caused by recurrent changes at particular nodes in the underlying gene regulatory networks (GRNs). The genes at such evolutionary 'hotspots' are thought to maximally affect the phenotype with minimal pleiotropic consequences. This has led to the suggestion that if a GRN is understood in sufficient detail, the path of evolution may be predictable. The repeated evolutionary loss of larval trichomes among Drosophila species is caused by the loss of shavenbaby (svb) expression. svb is also required for development of leg trichomes, but the evolutionary gain of trichomes in the 'naked valley' on T2 femurs in Drosophila melanogaster is caused by reduced microRNA-92a (miR-92a) expression rather than changes in svb. We compared the expression and function of components between the larval and leg trichome GRNs to investigate why the genetic basis of trichome pattern evolution differs in these developmental contexts. We found key differences between the two networks in both the genes employed, and in the regulation and function of common genes. These differences in the GRNs reveal why mutations in svb are unlikely to contribute to leg trichome evolution and how instead miR-92a represents the key evolutionary switch in this context. Our work shows that variability in GRNs across different developmental contexts, as well as whether a morphological feature is lost versus gained, influence the nodes at which a GRN evolves to cause morphological change. Therefore, our findings have important implications for understanding the pathways and predictability of evolution.

Morphological evolution caused by many subtle-effect substitutions in regulatory DNA.

Morphology evolves often through changes in developmental genes, but the causal mutations, and their effects, remain largely unknown. The evolution of naked cuticle on larvae of Drosophila sechellia resulted from changes in five transcriptional enhancers of shavenbaby (svb), a transcript of the ovo locus that encodes a transcription factor that governs morphogenesis of microtrichiae, hereafter called 'trichomes'. Here we show that the function of one of these enhancers evolved through multiple single-nucleotide substitutions that altered both the timing and level of svb expression. The consequences of these nucleotide substitutions on larval morphology were quantified with a novel functional assay. We found that each substitution had a relatively small phenotypic effect, and that many nucleotide changes account for this large morphological difference. In addition, we observed that the substitutions had non-additive effects. These data provide unprecedented resolution of the phenotypic effects of substitutions and show how individual nucleotide changes in a transcriptional enhancer have caused morphological evolution.

Phenotypic robustness conferred by apparently redundant transcriptional enhancers.

Genes include cis-regulatory regions that contain transcriptional enhancers. Recent reports have shown that developmental genes often possess multiple discrete enhancer modules that drive transcription in similar spatio-temporal patterns: primary enhancers located near the basal promoter and secondary, or 'shadow', enhancers located at more remote positions. It has been proposed that the seemingly redundant activity of primary and secondary enhancers contributes to phenotypic robustness. We tested this hypothesis by generating a deficiency that removes two newly discovered enhancers of shavenbaby (svb, a transcript of the ovo locus), a gene encoding a transcription factor that directs development of Drosophila larval trichomes. At optimal temperatures for embryonic development, this deficiency causes minor defects in trichome patterning. In embryos that develop at both low and high extreme temperatures, however, absence of these secondary enhancers leads to extensive loss of trichomes. These temperature-dependent defects can be rescued by a transgene carrying a secondary enhancer driving transcription of the svb cDNA. Finally, removal of one copy of wingless, a gene required for normal trichome patterning, causes a similar loss of trichomes only in flies lacking the secondary enhancers. These results support the hypothesis that secondary enhancers contribute to phenotypic robustness in the face of environmental and genetic variability.

Parental age influences developmental stability of the progeny in Drosophila.

The stochastic nature of biochemical processes is a source of variability that influences developmental stability. Developmental instability (DI) is often estimated through fluctuating asymmetry (FA), a parameter that deals with within-individual variation in bilateral structures. A relevant goal is to shed light on how environment, physiology and genotype relate to DI, thus providing a more comprehensive view of organismal development. Using Drosophila melanogaster isogenic lines, we investigated the effect of parental age, parental diet and offspring heterozygosity on DI. In this work, we have uncovered a clear relationship between parental age and offspring asymmetry. We show that asymmetry of the progeny increases concomitantly with parental age. Moreover, we demonstrate that enriching the diet of parents mitigates the effect of age on offspring symmetry. We show as well that increasing the heterozygosity of the progeny eliminates the effect of parental age on offspring symmetry. Taken together, our results suggest that diet, genotype and age of the parents interact to determine offspring DI in wild populations. These findings provide us with an avenue to understand the mechanisms underlying DI.

Conserved regulatory architecture underlies parallel genetic changes and convergent phenotypic evolution.

Similar morphological, physiological, and behavioral features have evolved independently in different species, a pattern known as convergence. It is known that morphological convergence can occur through changes in orthologous genes. In some cases of convergence, cis-regulatory changes generate parallel modifications in the expression patterns of orthologous genes. Our understanding of how changes in cis-regulatory regions contribute to convergence is hampered, usually, by a limited understanding of the global cis-regulatory structure of the evolving genes. Here we examine the genetic causes of a case of precise phenotypic convergence between Drosophila sechellia and Drosophila ezoana, species that diverged ~40 Mya. Previous studies revealed that changes in multiple transcriptional enhancers of shavenbaby (svb, a transcript of the ovo locus) caused phenotypic evolution in the D. sechellia lineage. It has also been shown that the convergent phenotype of D. ezoana was likely caused by cis-regulatory evolution of svb. Here we show that the large-scale cis-regulatory architecture of svb is conserved between these Drosophila species. Furthermore, we show that the D. ezoana orthologs of the evolved D. sechellia enhancers have also evolved expression patterns that correlate precisely with the changes in the phenotype. Our results suggest that phenotypic convergence resulted from multiple noncoding changes that occurred in parallel in the D. sechellia and D. ezoana lineages.

ASR1 mediates glucose-hormone cross talk by affecting sugar trafficking in tobacco plants.

Asr (for ABA, stress, ripening) genes are exclusively found in the genomes of higher plants, and the encoded proteins have been found localized both to the nucleus and cytoplasm. However, before the mechanisms underlying the activity of ASR proteins can be determined, the role of these proteins in planta should be deciphered. Results from this study suggest that ASR is positioned within the signaling cascade of interactions among glucose, abscisic acid, and gibberellins. Tobacco (Nicotiana tabacum) transgenic lines with reduced levels of ASR protein showed impaired glucose metabolism and altered abscisic acid and gibberellin levels. These changes were associated with dwarfism, reduced carbon dioxide assimilation, and accelerated leaf senescence as a consequence of a fine regulation exerted by ASR to the glucose metabolism. This regulation resulted in an impact on glucose signaling mediated by Hexokinase1 and Snf1-related kinase, which would subsequently have been responsible for photosynthesis, leaf senescence, and hormone level alterations. It thus can be postulated that ASR is not only involved in the control of hexose uptake in heterotrophic organs, as we have previously reported, but also in the control of carbon fixation by the leaves mediated by a similar mechanism.

Multiple layers of complexity in cis-regulatory regions of developmental genes.

Genomes contain the necessary information to ensure that genes are expressed in the right place, at the right time, and with the proper rate. Metazoan developmental genes often possess long stretches of DNA flanking their coding sequences and/or large introns which contain elements that influence gene expression. Most of these regulatory elements are relatively small and can be studied in isolation. For example, transcriptional enhancers, the elements that generate the expression pattern of a gene, have been traditionally studied with reporter constructs in transgenic animals. These studies have provided and will provide invaluable insights into enhancer evolution and function. However, this experimental approach has its limits; often, enhancer elements do not faithfully recapitulate native expression patterns. This fact suggests that additional information in cis-regulatory regions modulates the activity of enhancers and other regulatory elements. Indeed, recent studies have revealed novel functional aspects at the level of whole cis-regulatory regions. First, the discovery of "shadow enhancers." Second, the ubiquitous interactions between cis-regulatory elements. Third, the notion that some cis-regulatory regions may not function in a modular manner. Last, the effect of chromatin conformation on cis-regulatory activity. In this article, I describe these recent findings and discuss open questions in the field.

Pleiotropic Enhancers are Ubiquitous Regulatory Elements in the Human Genome.

Enhancers are regulatory elements of genomes that determine spatio-temporal patterns of gene expression. The human genome contains a vast number of enhancers, which largely outnumber protein-coding genes. Historically, enhancers have been regarded as highly tissue-specific. However, recent evidence has demonstrated that many enhancers are pleiotropic, with activity in multiple developmental contexts. Yet, the extent and impact of pleiotropy remain largely unexplored. In this study we analyzed active enhancers across human organs based on the analysis of both eRNA transcription (FANTOM5 consortium data sets) and chromatin architecture (ENCODE consortium data sets). We show that pleiotropic enhancers are pervasive in the human genome and that most enhancers active in a particular organ are also active in other organs. In addition, our analysis suggests that the proportion of context-specific enhancers of a given organ is explained, at least in part, by the proportion of context-specific genes in that same organ. The notion that such a high proportion of human enhancers can be pleiotropic suggests that small regions of regulatory DNA contain abundant regulatory information and that these regions evolve under important evolutionary constraints.

Can changes in 3D genome architecture create new regulatory landscapes that contribute to phenotypic evolution?

Animal genomes are compartmentalized into insulated regulatory units named topology-associated domains (TADs). TADs insulate gene promoters from enhancers that occupy neighboring TADs. Chromosomal rearrangements that disrupt TAD structure can generate new regulatory interactions between enhancers and promoters that were once separated into different TADs, which might lead to new gene expression patterns. On the one hand, TAD rearrangements are known to cause deleterious phenotypes, but, on the other hand, rearrangements can also create novel expression patterns that may be selected during evolution because they generate advantageous phenotypes. Here, we review recent studies that explore the effects of chromosomal rearrangements and genetic perturbations on TAD structure and gene regulation in the context of development and evolution. We discuss the possible contribution of evolutionary breakpoints (EBRs) that affect TAD structure to the evolution of gene regulation and the phenotype.

A complex gene regulatory architecture underlies the development and evolution of cuticle morphology in Drosophila.

The cuticle of insects is decorated with non-sensory hairs called trichomes. A few Drosophila species independently lost most of the dorso-lateral trichomes on first instar larvae. Genetic experiments revealed that this naked cuticle phenotype was caused by the evolution of enhancer function at the ovo/shavenbaby (ovo/svb) locus. Here we explore how this discovery catalyzed major new insights into morphological evolution in different developmental contexts, enhancer pleiotropy in gene regulation and the functionality and evolution of the Svb gene regulatory network (GRN). Taken together this highlights the importance of understanding the architecture and evolution of gene regulatory networks in detail and the great potential for further study of the Svb GRN.

Cis-regulatory variation in the shavenbaby gene underlies intraspecific phenotypic variation, mirroring interspecific divergence in the same trait.

Despite considerable progress in recent decades in dissecting the genetic causes of natural morphological variation, there is limited understanding of how variation within species ultimately contributes to species differences. We have studied patterning of the non-sensory hairs, commonly known as "trichomes," on the dorsal cuticle of first-instar larvae of Drosophila. Most Drosophila species produce a dense lawn of dorsal trichomes, but a subset of these trichomes were lost in D. sechellia and D. ezoana due entirely to regulatory evolution of the shavenbaby (svb) gene. Here, we describe intraspecific variation in dorsal trichome patterns of first-instar larvae of D. virilis that is similar to the trichome pattern variation identified previously between species. We found that a single large effect QTL, which includes svb, explains most of the trichome number difference between two D. virilis strains and that svb expression correlates with the trichome difference between strains. This QTL does not explain the entire difference between strains, implying that additional loci contribute to variation in trichome numbers. Thus, the genetic architecture of intraspecific variation exhibits similarities and differences with interspecific variation that may reflect differences in long-term and short-term evolutionary processes.

Unraveling the ecological processes modulating the population structure of Escherichia coli in a highly polluted urban stream network.

Escherichia coli dynamics in urban watersheds are affected by a complex balance among external inputs, niche modulation and genetic variability. To explore the ecological processes influencing E. coli spatial patterns, we analyzed its abundance and phylogenetic structure in water samples from a stream network with heterogeneous urban infrastructure and environmental conditions. Our results showed that environmental and infrastructure variables, such as macrophyte coverage, DIN and sewerage density, mostly explained E. coli abundance. Moreover, main generalist phylogroups A and B1 were found in high proportion, which, together with an observed negative relationship between E. coli abundance and phylogroup diversity, suggests that their dominance might be due to competitive exclusion. Lower frequency phylogroups were associated with sites of higher ecological disturbance, mainly involving simplified habitats, higher drainage infrastructure and septic tank density. In addition to the strong negative relationship between phylogroup diversity and dominance, the occurrence of these phylogroups would be associated with increased facilitated dispersal. Nutrients also contributed to explaining phylogroup distribution. Our study proposes the differential contribution of distinct ecological processes to the patterns of E. coli in an urban watershed, which is useful for the monitoring and management of fecal pollution.

Nucleotide polymorphism in the drought responsive gene Asr2 in wild populations of tomato.

The Asr gene family (named after abscicic acid [ABA], stress, ripening), exclusively present in plant genomes, is involved in transcriptional regulation. Its members are up-regulated in roots and leaves of water- or salt-stressed plants. In previous work, evidence of adaptive evolution (as inferred from synonymous and nonsynonymous divergence rates) has been reported for Asr2 in Solanum chilense and S. arcanum, two species dwelling in habitats with different precipitation regimes. In this paper we investigate patterns of intraspecific nucleotide variation in Asr2 and the unlinked locus CT114 in S. chilense and S. arcanum. The extent of nucleotide diversity in Asr2 differed between species in more than one order of magnitude. In both species we detected evidence of non-neutral evolution, which may be ascribed to different selective regimes, potentially associated to unique climatic features, or, alternatively, to demographic events. The results are discussed in the light of demographic and selective hypotheses.

Comprehensive Analysis of a cis-Regulatory Region Reveals Pleiotropy in Enhancer Function.

Developmental genes can have complex cis-regulatory regions with multiple enhancers. Early work revealed remarkable modularity of enhancers, whereby distinct DNA regions drive gene expression in defined spatiotemporal domains. Nevertheless, a few reports have shown that enhancers function in multiple developmental stages, implying that enhancers can be pleiotropic. Here, we have studied the activity of the enhancers of the shavenbaby gene throughout D. melanogaster development. We found that all seven shavenbaby enhancers drive expression in multiple tissues and developmental stages. We explored how enhancer pleiotropy is encoded in two of these enhancers. In one enhancer, the same transcription factor binding sites contribute to embryonic and pupal expression, revealing site pleiotropy, whereas for a second enhancer, these roles are encoded by distinct sites. Enhancer pleiotropy may be a common feature of cis-regulatory regions of developmental genes, and site pleiotropy may constrain enhancer evolution in some cases.

Genetic and Transgenic Reagents for Drosophila simulans, D. mauritiana, D. yakuba, D. santomea, and D. virilis.

Species of the Drosophila melanogaster species subgroup, including the species D. simulans, D. mauritiana, D. yakuba, and D. santomea, have long served as model systems for studying evolution. However, studies in these species have been limited by a paucity of genetic and transgenic reagents. Here, we describe a collection of transgenic and genetic strains generated to facilitate genetic studies within and between these species. We have generated many strains of each species containing mapped piggyBac transposons including an enhanced yellow fluorescent protein (EYFP) gene expressed in the eyes and a ϕC31 attP site-specific integration site. We have tested a subset of these lines for integration efficiency and reporter gene expression levels. We have also generated a smaller collection of other lines expressing other genetically encoded fluorescent molecules in the eyes and a number of other transgenic reagents that will be useful for functional studies in these species. In addition, we have mapped the insertion locations of 58 transposable elements in D. virilis that will be useful for genetic mapping studies.

Evolutionary history of the Asr gene family.

The Asr gene family is widespread in higher plants. Most Asr genes are up-regulated under different environmental stress conditions and during fruit ripening. ASR proteins are localized in the nucleus and their likely function is transcriptional regulation. In cultivated tomato, we identified a novel fourth family member, named Asr4, which maps close to its sibling genes Asr1-Asr2-Asr3 and displays an unshared region coding for a domain containing a 13-amino acid repeat. In this work we were able to expand our previous analysis for Asr2 and investigated the coding regions of the four known Asr paralogous genes in seven tomato species from different geographic locations. In addition, we performed a phylogenetic analysis on ASR proteins. The first conclusion drawn from this work is that tomato ASR proteins cluster together in the tree. This observation can be explained by a scenario of concerted evolution or birth and death of genes. Secondly, our study showed that Asr1 is highly conserved at both replacement and synonymous sites within the genus Lycopersicon. ASR1 protein sequence conservation might be associated with its multiple functions in different tissues while the low rate of synonymous substitutions suggests that silent variation in Asr1 is selectively constrained, which is probably related to its high expression levels. Finally, we found that Asr1 activation under water stress is not conserved between Lycopersicon species.