Correlated evolution of the cis-acting regulatory elements and developmental expression of the Drosophila Gld gene in seven species from the subgroup melanogaster

Notch signaling generates the "cut here line" on the cuticle of the puparium in Drosophila melanogaster

During a molt or eclosion, insects shed their cuticle, an extracellular matrix made by underlying epidermal cells, by cleavage along a defined line. This means that the "cut here line" is pre-formed on the cuticle, and its formation is indispensable for insect life. Here, we show that the proper formation of the operculum ridge (OR), which is the "cut here line" on the puparium (pupal case) of Drosophila melanogaster, involves Notch signaling activation in the epidermal cells just beneath the future OR region (OR-forming cells). The inhibition of Notch signaling causes defects in eclosion due to failure in OR cleavage, the chitin organization and several cuticular proteins localization, glucose dehydrogenase (Gld) activity, and OR-forming cell shape. Our findings provide the first insight into the molecular basis of the structure and formation of the "cut here line" on the cuticle.

Selection and geography shape male reproductive tract transcriptomes in Drosophila melanogaster

Transcriptome analysis of several animal clades suggests that male reproductive tract gene expression evolves quickly. However, the factors influencing the abundance and distribution of within-species variation, the ultimate source of interspecific divergence, are poorly known. Drosophila melanogaster, an ancestrally African species that has recently spread throughout the world and colonized the Americas in the last roughly 100 years, exhibits phenotypic and genetic latitudinal clines on multiple continents, consistent with a role for spatially varying selection in shaping its biology. Nevertheless, geographic expression variation in the Americas is poorly described, as is its relationship to African expression variation. Here, we investigate these issues through the analysis of two male reproductive tissue transcriptomes [testis and accessory gland (AG)] in samples from Maine (USA), Panama, and Zambia. We find dramatic differences between these tissues in differential expression between Maine and Panama, with the accessory glands exhibiting abundant expression differentiation and the testis exhibiting very little. Latitudinal expression differentiation appears to be influenced by the selection of Panama expression phenotypes. While the testis shows little latitudinal expression differentiation, it exhibits much greater differentiation than the accessory gland in Zambia vs American population comparisons. Expression differentiation for both tissues is non-randomly distributed across the genome on a chromosome arm scale. Interspecific expression divergence between D. melanogaster and D. simulans is discordant with rates of differentiation between D. melanogaster populations. Strongly heterogeneous expression differentiation across tissues and timescales suggests a complex evolutionary process involving major temporal changes in the way selection influences expression evolution in these organs.

Conserved noncoding sequences provide insights into regulatory sequence and loss of gene expression in maize

Thousands of species will be sequenced in the next few years; however, understanding how their genomes work, without an unlimited budget, requires both molecular and novel evolutionary approaches. We developed a sensitive sequence alignment pipeline to identify conserved noncoding sequences (CNSs) in the Andropogoneae tribe (multiple crop species descended from a common ancestor ∼18 million years ago). The Andropogoneae share similar physiology while being tremendously genomically diverse, harboring a broad range of ploidy levels, structural variation, and transposons. These contribute to the potential of Andropogoneae as a powerful system for studying CNSs and are factors we leverage to understand the function of maize CNSs. We found that 86% of CNSs were comprised of annotated features, including introns, UTRs, putative cis-regulatory elements, chromatin loop anchors, noncoding RNA (ncRNA) genes, and several transposable element superfamilies. CNSs were enriched in active regions of DNA replication in the early S phase of the mitotic cell cycle and showed different DNA methylation ratios compared to the genome-wide background. More than half of putative cis-regulatory sequences (identified via other methods) overlapped with CNSs detected in this study. Variants in CNSs were associated with gene expression levels, and CNS absence contributed to loss of gene expression. Furthermore, the evolution of CNSs was associated with the functional diversification of duplicated genes in the context of maize subgenomes. Our results provide a quantitative understanding of the molecular processes governing the evolution of CNSs in maize.

Polymorphism and Divergence of Novel Gene Expression Patterns in Drosophila melanogaster

Transcriptomes may evolve by multiple mechanisms, including the evolution of novel genes, the evolution of transcript abundance, and the evolution of cell, tissue, or organ expression patterns. Here, we focus on the last of these mechanisms in an investigation of tissue and organ shifts in gene expression in Drosophila melanogaster. In contrast to most investigations of expression evolution, we seek to provide a framework for understanding the mechanisms of novel expression patterns on a short population genetic timescale. To do so, we generated population samples of D. melanogaster transcriptomes from five tissues: accessory gland, testis, larval salivary gland, female head, and first-instar larva. We combined these data with comparable data from two outgroups to characterize gains and losses of expression, both polymorphic and fixed, in D. melanogaster We observed a large number of gain- or loss-of-expression phenotypes, most of which were polymorphic within D. melanogaster Several polymorphic, novel expression phenotypes were strongly influenced by segregating cis-acting variants. In support of previous literature on the evolution of novelties functioning in male reproduction, we observed many more novel expression phenotypes in the testis and accessory gland than in other tissues. Additionally, genes showing novel expression phenotypes tend to exhibit greater tissue-specific expression. Finally, in addition to qualitatively novel expression phenotypes, we identified genes exhibiting major quantitative expression divergence in the D. melanogaster lineage.

What does it take to evolve an enhancer? A simulation-based study of factors influencing the emergence of combinatorial regulation

There is widespread interest today in understanding enhancers, which are regulatory elements typically harboring several transcription factor binding sites and mediating the combinatorial effect of transcription factors on gene expression. The evolution of enhancers poses interesting unanswered questions, for example, the evolutionary time taken for a typical enhancer to emerge or the factors shaping its evolution. Existing approaches to cis-regulatory evolution have often ignored the combinatorial nature and varied biochemical mechanisms of gene regulation encoded in enhancers. We report on our investigation of enhancer evolution through the use of PEBCRES, a framework for evolutionary simulation of enhancers that employs a mechanistic and well-supported sequence-to-expression model to assign fitness to the evolving enhancer genotype. We estimated the time necessary to evolve, from genomic background, enhancers capable of driving complex gene expression patterns similar to those involved in early development in Drosophila. We found the time-to-evolve to range between 0.5 and 10 Myr, and to vary greatly with the target expression pattern, complexity of the real enhancer known to encode that pattern, and the strength of input from specific transcription factors. To our knowledge, this is the first estimate of waiting times for realistic enhancers to evolve. The in silico evolved enhancers had, with a few interesting exceptions, site compositions similar to those seen in real enhancers for the same patterns. Our simulations also revealed that certain features of an enhancer might evolve not due to their biological function but as aids to the evolutionary process itself.

Evolution of cis-regulatory sequence and function in Diptera

Cis-regulatory sequences direct patterns of gene expression essential for development and physiology. Evolutionary changes in these sequences contribute to phenotypic divergence. Despite their importance, cis-regulatory regions remain one of the most enigmatic features of the genome. Patterns of sequence evolution can be used to identify cis-regulatory elements, but the power of this approach depends upon the relationship between sequence and function. Comparative studies of gene regulation among Diptera reveal that divergent sequences can underlie conserved expression, and that expression differences can evolve despite largely similar sequences. This complex structure-function relationship is the primary impediment for computational identification and interpretation of cis-regulatory sequences. Biochemical characterization and in vivo assays of cis-regulatory sequences on a genomic-scale will relieve this barrier.

Rate of promoter class turn-over in yeast evolution

BACKGROUND

Phylogenetic conservation at the DNA level is routinely used as evidence of molecular function, under the assumption that locations and sequences of functional DNA segments remain invariant in evolution. In particular, short DNA segments participating in initiation and regulation of transcription are often conserved between related species. However, transcription of a gene can evolve, and this evolution may involve changes of even such conservative DNA segments. Genes of yeast Saccharomyces have promoters of two classes, class 1 (TATA-containing) and class 2 (non-TATA-containing).

RESULTS

Comparison of upstream non-coding regions of orthologous genes from the five species of Saccharomyces sensu stricto group shows that among 212 genes which very likely have class 1 promoters in S. cerevisiae, 17 probably have class 2 promoters in one or more other species. Conversely, among 322 genes which very likely have class 2 promoters in S. cerevisiae, 44 probably have class 1 promoters in one or more other species. Also, for at least 2 genes from the set of 212 S. cerevisiae genes with class 1 promoters, the locations of the TATA consensus sequences are substantially different between the species.

CONCLUSION

Our results indicate that, in the course of yeast evolution, a promoter switches its class with the probability at least approximately 0.1 per time required for the accumulation of one nucleotide substitution at a non-coding site. Thus, key sequences involved in initiation of transcription evolve with substantial rates in yeast.

Is ectopic expression caused by deregulatory mutations or due to gene-regulation leaks with evolutionary potential?

It has long been thought that gene expression is tightly regulated in multicellular eukaryotes, so that expression profiles match functional profiles. This conception emerged from the assumption that gene activity is synonymous with gene function. This paradigm was first challenged by comparative protein electrophoresis studies showing extensive differences in expression patterns among related species. The paradigm is now being challenged by evolutionary transcriptomics using microarray technologies. Most gene expression profiles display features that lack any obvious functional significance. The so-called "ectopic" expression refers to the expression of genes at times and locations where the target gene is not known to have a function. However, ectopic expression might be associated with genuine function even if this function is not essential or has yet to be ascertained. Alternatively, ectopic expression might come about as a superfluous by-product of regulatory systems, which would call for a revision of prevailing ideas about the specificity of gene regulation. We herein review available evidence for ectopic expression and the hypotheses proposed for its origin and evolution. We propose that ectopic expression must be regarded as part of an integrated phenotypic whole. It seems likely that ectopic expression represents a leak in the evolution of regulatory systems, but one that is endowed with considerable evolutionary possibilities.

Molecular phylogeny of the Drosophila melanogaster species subgroup

Although molecular and phenotypic evolution have been studied extensively in Drosophila melanogaster and its close relatives, phylogenetic relationships within the D. melanogaster species subgroup remain unresolved. In particular, recent molecular studies have not converged on the branching orders of the D. yakuba-D. teissieri and D. erecta-D. orena species pairs relative to the D. melanogaster-D. simulans-D. mauritiana-D. sechellia species complex. Here, we reconstruct the phylogeny of the melanogaster species subgroup using DNA sequence data from four nuclear genes. We have employed "vectorette PCR" to obtain sequence data for orthologous regions of the Alcohol dehydrogenase (Adh), Alcohol dehydrogenase related (Adhr), Glucose dehydrogenase (Gld), and rosy (ry) genes (totaling 7164 bp) from six melanogaster subgroup species (D. melanogaster, D. simulans, D. teissieri, D. yakuba, D. erecta, and D. orena) and three species from subgroups outside the melanogaster species subgroup [D. eugracilis (eugracilis subgroup), D. mimetica (suzukii subgroup), and D. lutescens (takahashii subgroup)]. Relationships within the D. simulans complex are not addressed. Phylogenetic analyses employing maximum parsimony, neighbor-joining, and maximum likelihood methods strongly support a D. yakuba-D. teissieri and D. erecta-D. orena clade within the melanogaster species subgroup. D. eugracilis is grouped closer to the melanogaster subgroup than a D. mimetica-D. lutescens clade. This tree topology is supported by reconstructions employing simple (single parameter) and more complex (nonreversible) substitution models.

Gene co-option in physiological and morphological evolution

Co-option occurs when natural selection finds new uses for existing traits, including genes, organs, and other body structures. Genes can be co-opted to generate developmental and physiological novelties by changing their patterns of regulation, by changing the functions of the proteins they encode, or both. This often involves gene duplication followed by specialization of the resulting paralogous genes into particular functions. A major role for gene co-option in the evolution of development has long been assumed, and many recent comparative developmental and genomic studies have lent support to this idea. Although there is relatively less known about the molecular basis of co-option events involving developmental pathways, much can be drawn from well-studied examples of the co-option of structural proteins. Here, we summarize several case studies of both structural gene and developmental genetic circuit co-option and discuss how co-option may underlie major episodes of adaptive change in multicellular organisms. We also examine the phenomenon of intraspecific variability in gene expression patterns, which we propose to be one form of material for the co-option process. We integrate this information with recent models of gene family evolution to provide a framework for understanding the origin of co-optive evolution and the mechanisms by which natural selection promotes evolutionary novelty by inventing new uses for the genetic toolkit.

Evolution of functionally conserved enhancers can be accelerated in large populations: a population-genetic model

The evolution of cis-regulatory elements (or enhancers) appears to proceed at dramatically different rates in different taxa. Vertebrate enhancers are often very highly conserved in their sequences, and relative positions, across distantly related taxa. In contrast, functionally equivalent enhancers in closely related Drosophila species can differ greatly in their sequences and spatial organization. We present a population-genetic model to explain this difference. The model examines the dynamics of fixation of pairs of individually deleterious, but compensating, mutations. As expected, small populations are predicted to have a high rate of evolution, and the rate decreases with increasing population size. In contrast to previous models, however, this model predicts that the rate of evolution by pairs of compensatory mutations increases dramatically for population sizes above several thousand individuals, to the point of greatly exceeding the neutral rate. Application of this model predicts that species with moderate population sizes will have relatively conserved enhancers, whereas species with larger populations will be expected to evolve their enhancers at much higher rates. We propose that the different degree of conservation seen in vertebrate and Drosophila enhancers may be explained solely by differences in their population sizes and generation times.

Rapid evolution of cis-regulatory sequences via local point mutations

Although the evolution of protein-coding sequences within genomes is well understood, the same cannot be said of the cis-regulatory regions that control transcription. Yet, changes in gene expression are likely to constitute an important component of phenotypic evolution. We simulated the evolution of new transcription factor binding sites via local point mutations. The results indicate that new binding sites appear and become fixed within populations on microevolutionary timescales under an assumption of neutral evolution. Even combinations of two new binding sites evolve very quickly. We predict that local point mutations continually generate considerable genetic variation that is capable of altering gene expression.

Complex organization of promoter and enhancer elements regulate the tissue- and developmental stage-specific expression of the Drosophila melanogaster Gld gene

The Drosophila melanogaster Gld gene has multiple and diverse developmental and physiological functions. We report herein that interactions among proximal promoter elements and a cluster of intronically located enhancers and silencers specify the complex regulation of Gld that underlies its diverse functions. Gld expression in nonreproductive tissues is largely determined by proximal promoter elements with the exception of the embryonic labium where Gld is activated by an enhancer within the first intron. A nuclear protein, GPAL, has been identified that binds the Gpal elements in the proximal promoter region. Regulation of Gld in the reproductive organs is particularly complex, involving interactions among the Gpal proximal promoter elements, a unique TATA box, three distinct enhancer types, and one or more silencer elements. The three somatic reproductive organ enhancers each activate expression in male and female pairs of reproductive organs. One of these pairs, the male ejaculatory duct and female oviduct, are known to be developmentally homologous. We report evidence that the other two pairs of organs are developmentally homologous as well. A comprehensive model to explain the full developmental regulation of Gld and its evolution is presented.

Evolutionary developmental biology and the problem of variation

One of the oldest problems in evolutionary biology remains largely unsolved. Which mutations generate evolutionarily relevant phenotypic variation? What kinds of molecular changes do they entail? What are the phenotypic magnitudes, frequencies of origin, and pleiotropic effects of such mutations? How is the genome constructed to allow the observed abundance of phenotypic diversity? Historically, the neo-Darwinian synthesizers stressed the predominance of micromutations in evolution, whereas others noted the similarities between some dramatic mutations and evolutionary transitions to argue for macromutationism. Arguments on both sides have been biased by misconceptions of the developmental effects of mutations. For example, the traditional view that mutations of important developmental genes always have large pleiotropic effects can now be seen to be a conclusion drawn from observations of a small class of mutations with dramatic effects. It is possible that some mutations, for example, those in cis-regulatory DNA, have few or no pleiotropic effects and may be the predominant source of morphological evolution. In contrast, mutations causing dramatic phenotypic effects, although superficially similar to hypothesized evolutionary transitions, are unlikely to fairly represent the true path of evolution. Recent developmental studies of gene function provide a new way of conceptualizing and studying variation that contrasts with the traditional genetic view that was incorporated into neo-Darwinian theory and population genetics. This new approach in developmental biology is as important for microevolutionary studies as the actual results from recent evolutionary developmental studies. In particular, this approach will assist in the task of identifying the specific mutations generating phenotypic variation and elucidating how they alter gene function. These data will provide the current missing link between molecular and phenotypic variation in natural populations.

Generation of evolutionary novelty by functional shift

That biological features may change their function during evolution has long been recognized. Particularly, the acquisition of new functions by molecules involved in developmental pathways is suspected to cause important morphologic novelties. However, the current terminology describing functional changes during evolution (co-option or recruitment) fails to recognize important biologic distinctions between diverse evolutionary routes involving functional shifts. The main goal of our work is to stress the importance of an apparently trivial distinction: Whether or not the element that adopts a new function (anything from a morphologic structure to a protein domain) is a single or a duplicated element. We propose that natural selection must act in a radically different way, depending on the historic succession of co-option and duplication events; that is, co-option may provide the selective pressure for a subsequent gene duplication or could be a stabilizing factor that helps maintain redundancy after gene duplication. We review the evidence available on functional changes, focusing whenever possible on developmental molecules, and we propose a conceptual framework for the study of functional shifts during evolution with a level of resolution appropriate to the power of our current methodologies.