Repeated morphological evolution through cis-regulatory changes in a pleiotropic gene

Microarray analysis of replicate populations selected against a wing-shape correlation in Drosophila melanogaster

We selected bidirectionally to change the phenotypic correlation between two wing dimensions in Drosophila melanogaster and measured gene expression differences in late third instar wing disks, using microarrays. We tested an array of 12 selected lines, including 10 from a Massachusetts population (5 divergently selected pairs) and 2 from a California population (1 divergently selected pair). In the Massachusetts replicates, 29 loci showed consistent, significant expression differences in all 5 line-pair comparisons. However, the significant loci in the California lines were almost completely different from these. The disparity between responding genes in different gene pools confirms recent evidence that surprisingly large numbers of loci can affect wing shape. Our results also show that with well-replicated selection lines, of large effective size, the numbers of candidate genes in microarray-based searches can be reduced to realistic levels.

Variation in neural V1aR predicts sexual fidelity and space use among male prairie voles in semi-natural settings

Although prairie voles (Microtus ochrogaster) are socially monogamous, males vary in both sexual and spatial fidelity. Most males form pairbonds, cohabit with one female, and defend territories. Wandering males, in contrast, have expansive home ranges that overlap many males and females. In the laboratory, pairing is regulated by arginine vasopressin and its predominant CNS receptor, vasopressin 1a receptor (V1aR). We investigated individual differences in forebrain V1aR expression of male prairie voles in mixed-sex seminatural enclosures. Individual differences in V1aR were compared with space use measured by radio telemetry and paternity determined with microsatellite markers. Animals engaging in extra-pair fertilizations (EPFs) as either wanderers or paired residents overlapped significantly more in same- and opposite-sex home ranges. Surprisingly, neither social fidelity measured by space use nor sexual fidelity measured by paternity was associated with V1aR expression in the ventral pallidum (VPall) or lateral septum, areas causally related to pairbond formation. In contrast, V1aR expression in the posterior cingulate/retrosplenial cortex (PCing) and laterodorsal thalamus (LDThal), areas implicated in spatial memory, strongly covaried with space use and paternity. Animals engaging in EPFs either as wanderers or paired residents exhibited low levels of LDThal and PCing V1aR expression. Individual differences in brain and behavior parallel differences between prairie voles and promiscuous congeners. The concordance among space use, paternity, and V1aR in spatial circuits suggests a common link between the mechanisms of spatial behaviors and success at EPF. The combined data demonstrate how organismal biology can inform our understanding of individual and species differences in behavioral mechanisms.

Evolutionary analysis of the cis-regulatory region of the spicule matrix gene SM50 in strongylocentrotid sea urchins

An evolutionary analysis of transcriptional regulation is essential to understanding the molecular basis of phenotypic diversity. The sea urchin is an ideal system in which to explore the functional consequence of variation in cis-regulatory sequences. We are particularly interested in the evolution of genes involved in the patterning and synthesis of its larval skeleton. This study focuses on the cis-regulatory region of SM50, which has already been characterized to a considerable extent in the purple sea urchin, Strongylocentrotus purpuratus. We have isolated the cis-regulatory region from 15 individuals of S. purpuratus as well as seven closely related species in the family Strongylocentrotidae. We have performed a variety of statistical tests and present evidence that the cis-regulatory elements upstream of the SM50 gene have been subject to positive selection along the lineage leading to S. purpuratus. In addition, we have performed electrophoretic mobility shift assays (EMSAs) and demonstrate that nucleotide substitutions within Element C affect the ability of nuclear proteins to bind to this cis-regulatory element among members of the family Strongylocentrotidae. We speculate that such changes in SM50 and other genes could accumulate to produce altered patterns of gene expression with functional consequences during skeleton formation.

cis-Regulatory changes in Kit ligand expression and parallel evolution of pigmentation in sticklebacks and humans

Dramatic pigmentation changes have evolved within most vertebrate groups, including fish and humans. Here we use genetic crosses in sticklebacks to investigate the parallel origin of pigmentation changes in natural populations. High-resolution mapping and expression experiments show that light gills and light ventrums map to a divergent regulatory allele of the Kit ligand (Kitlg) gene. The divergent allele reduces expression in gill and skin tissue and is shared by multiple derived freshwater populations with reduced pigmentation. In humans, Europeans and East Asians also share derived alleles at the KITLG locus. Strong signatures of selection map to regulatory regions surrounding the gene, and admixture mapping shows that the KITLG genomic region has a significant effect on human skin color. These experiments suggest that regulatory changes in Kitlg contribute to natural variation in vertebrate pigmentation, and that similar genetic mechanisms may underlie rapid evolutionary change in fish and humans.

Reliable prediction of regulator targets using 12 Drosophila genomes

Gene expression is regulated pre- and post-transcriptionally via cis-regulatory DNA and RNA motifs. Identification of individual functional instances of such motifs in genome sequences is a major goal for inferring regulatory networks yet has been hampered due to the motifs' short lengths that lead to many chance matches and poor signal-to-noise ratios. In this paper, we develop a general methodology for the comparative identification of functional motif instances across many related species, using a phylogenetic framework that accounts for the evolutionary relationships between species, allows for motif movements, and is robust against missing data due to artifacts in sequencing, assembly, or alignment. We also provide a robust statistical framework for evaluating motif confidence, which enables us to translate evolutionary conservation into a confidence measure for each motif instance, correcting for varying motif length, composition, and background conservation of the target regions. We predict targets of fly transcription factors and miRNAs in alignments of 12 recently sequenced Drosophila species. When compared to extensive genome-wide experimental data, predicted targets are of high quality, matching and surpassing ChIP-chip microarrays and recovering miRNA targets with high sensitivity. The resulting regulatory network suggests significant redundancy between pre- and post-transcriptional regulation of gene expression.

Simpler mode of inheritance of transcriptional variation in male Drosophila melanogaster

Sexual selection drives faster evolution in males. The X chromosome is potentially an important target for sexual selection, because hemizygosity in males permits accumulation of alleles, causing tradeoffs in fitness between sexes. Hemizygosity of the X could cause fundamentally different modes of inheritance between the sexes, with more additive variation in males and more nonadditive variation in females. Indeed, we find that genetic variation for the transcriptome is primarily additive in males but nonadditive in females. As expected, these differences are more pronounced on the X chromosome than the autosomes, but autosomal loci are also affected, possibly because of X-linked transcription factors. These differences may be of evolutionary significance because additive variation responds quickly to selection, whereas nonadditive genetic variation does not. Thus, hemizygosity of the X may underlie much of the faster male evolution of the transcriptome and potentially other phenotypes. Consistent with this prediction, genes that are additive in males and nonadditive in females are overrepresented among genes responding to selection for increased mating speed.

The limits of subfunctionalization

Background

The duplication-degeneration-complementation (DDC) model has been proposed as an explanation for the unexpectedly high retention of duplicate genes. The hypothesis proposes that, following gene duplication, the two gene copies degenerate to perform complementary functions that jointly match that of the single ancestral gene, a process also known as subfunctionalization. We distinguish between subfunctionalization at the regulatory level and at the product level (e.g within temporal or spatial expression domains).

Results

In contrast to what is expected under the DDC model, we use in silico modeling to show that regulatory subfunctionalization is expected to peak and then decrease significantly. At the same time, neofunctionalization (recruitment of novel interactions) increases monotonically, eventually affecting the regulatory elements of the majority of genes. Furthermore, since this process occurs under conditions of stabilizing selection, there is no need to invoke positive selection. At the product level, the frequency of subfunctionalization is no higher than would be expected by chance, a finding that was corroborated using yeast microarray time-course data. We also find that product subfunctionalization is not necessarily caused by regulatory subfunctionalization.

Conclusion

Our results suggest a more complex picture of post-duplication evolution in which subfunctionalization plays only a partial role in conjunction with redundancy and neofunctionalization. We argue that this behavior is a consequence of the high evolutionary plasticity in gene networks.

Identification and analysis of a conserved Tcfap2a intronic enhancer element required for expression in facial and limb bud mesenchyme

Tcfap2a, the gene encoding the mouse AP-2alpha transcription factor, is required for normal development of multiple structures during embryogenesis, including the face and limbs. Using comparative sequence analysis and transgenic-mouse experiments we have identified an intronic enhancer within this gene that directs expression to the face and limb mesenchyme. There are two conserved sequence blocks within this intron, and the larger of these directs tissue-specific activity and is found in all vertebrate Tcfap2a genes analyzed. To assess the role of the enhancer in regulating endogenous mouse Tcfap2a expression, we have deleted this cis-regulatory sequence from the genome. Loss of this element severely impairs Tcfap2a expression in the limb bud mesenchyme but generates only a modest reduction in the facial mesenchyme. The reduction in Tcfap2a transcription is accompanied by altered patterning of the forelimb, resulting in postaxial polydactyly. These results indicate that the major role for this enhancer resides within the limb bud, and it serves to maintain a level of Tcfap2a expression that limits the size of the hand plate and the associated number of digit primordia. The potential role of this cis-acting sequence in modeling the size and shape of the face and limbs during evolution is discussed.

Parallel evolution of facial stripe patterns in the Neolamprologus brichardi/pulcher species complex endemic to Lake Tanganyika

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past.

Proteomic analysis of N-glycosylation in mosquito dopachrome conversion enzyme

A novel dopachrome conversion enzyme (DCE) is present in insects and involved in their melanization pathway. DCE shares no sequence homology with any noninsect species from bacteria to humans. Several DCE sequences have been available, but enzyme structure and catalytic mechanism are unclear. This study concerns DCE PTMs, especially glycosylation. A mosquito DCE was purified and its monosaccharide composition, N-glycosylation site, and oligosaccharide structures were determined. Results showed that N-acetyl D-glucosamine and D-mannose are the major monosaccharides and L-fucose, D-xylose, and D-arabinose are the minor ones in mosquito DCE. Glycosylation site and oligosaccharide structures were elucidated from MS and MS/MS spectra of trypsin-digested DCE glycopeptides. A single N-glycosylation site (Asn285 -Glu-Thr) was identified in DCE and was proven to be fully glycosylated. Man3GlcNAc2, Man3(Fuc)1-2GlcNAc2, and their truncated structures were the dominant oligosaccharides. In addition, high mannose-type structures (Man4-7(Fuc)GlcNAc2) were also identified. Removal of DCE N-oligosaccharides with peptide N-glycosidase (PNGase F) decreased its activity and thermal stability. However, partial DCE deglycosylation with alpha-mannosidase or alpha-fucosidase somewhat stimulated its activity and improved its thermal stability. During mass spectrometric analysis of DCE glycopeptides, their CID patterns were highly intriguing, in that some glycopeptides underwent both C-terminal rearrangement and formation of dimeric structures during CID. Results of this study provide an interesting example in terms of potential complexity of the glycopeptide CID fragmentation pattern.

Constraints on utilization of the EDA-signaling pathway in threespine stickleback evolution

Many traits evolve in parallel in widely separated populations. The evolutionary radiation of threespine sticklebacks provides a powerful model for testing the molecular basis of parallel evolution in vertebrates. Although marine sticklebacks are completely covered with bony armor plates, most freshwater populations have dramatic reductions in plates. Recent genetic studies have shown that major changes in armor patterning are likely due to regulatory alterations in the gene encoding the secreted signaling molecule ectodysplasin (EDA). In mammals, mutations in many different components of the EDA-signaling pathway produce similar changes in hair, teeth, sweat glands, and dermal bones. To test whether other genes in the EDA pathway also control natural variation in armor plates, we identified and mapped stickleback EDA Receptor (EDAR), the EDAR-Associated Death Domain adaptor, Tumor Necrosis Factor Receptor (TNFR) SuperFamily member 19, its adaptor TNFR-Associated Factor 6, and the downstream regulator nuclear factor kappa B Essential Modulator (NEMO). In contrast to the diversity of genes underlying ectodermal dysplasia disease phenotypes in humans, none of these EDA pathway components map to chromosomes previously shown to modify armor plates in natural populations, though EDAR showed a small but significant effect on plate number. We further investigated whether these genes exhibit differences in copy number, target size, or genomic organization that might make them less suitable targets for evolutionary change. In comparison with EDA, all these genes have smaller surrounding noncoding (putative regulatory) regions, with fewer evolutionarily conserved regions. We suggest that the presence of highly modular cis-acting control sequences may be a key factor influencing the likelihood that particular genes will serve as the basis of major phenotypic changes in nature.

A multicassette Gateway vector set for high throughput and comparative analyses in ciona and vertebrate embryos

Background

The past few years have seen a vast increase in the amount of genomic data available for a growing number of taxa, including sets of full length cDNA clones and cis-regulatory sequences. Large scale cross-species comparisons of protein function and cis-regulatory sequences may help to understand the emergence of specific traits during evolution.

Principal Findings

To facilitate such comparisons, we developed a Gateway compatible vector set, which can be used to systematically dissect cis-regulatory sequences, and overexpress wild type or tagged proteins in a variety of chordate systems. It was developed and first characterised in the embryos of the ascidian Ciona intestinalis, in which large scale analyses are easier to perform than in vertebrates, owing to the very efficient embryo electroporation protocol available in this organism. Its use was then extended to fish embryos and cultured mammalian cells.

Conclusion

This versatile vector set opens the way to the mid- to large-scale comparative analyses of protein function and cis-regulatory sequences across chordate evolution. A complete user manual is provided as supplemental material.

The quest for natural selection in the age of comparative genomics

Continued genome sequencing has fueled progress in statistical methods for understanding the action of natural selection at the molecular level. This article reviews various statistical techniques (and their applicability) for detecting adaptation events and the functional divergence of proteins. As large-scale automated studies become more frequent, they provide a useful resource for generating biological null hypotheses for further experimental and statistical testing. Furthermore, they shed light on typical patterns of lineage-specific evolution of organisms, on the functional and structural evolution of protein families and on the interplay between the two. More complex models are being developed to better reflect the underlying biological and chemical processes and to complement simpler statistical models. Linking molecular processes to their statistical signatures in genomes can be demanding, and the proper application of statistical models is discussed.

Dating phylogenies with hybrid local molecular clocks

Background

Because rates of evolution and species divergence times cannot be estimated directly from molecular data, all current dating methods require that specific assumptions be made before inferring any divergence time. These assumptions typically bear either on rates of molecular evolution (molecular clock hypothesis, local clocks models) or on both rates and times (penalized likelihood, Bayesian methods). However, most of these assumptions can affect estimated dates, oftentimes because they underestimate large amounts of rate change.

Principal Findings

A significant modification to a recently proposed ad hoc rate-smoothing algorithm is described, in which local molecular clocks are automatically placed on a phylogeny. This modification makes use of hybrid approaches that borrow from recent theoretical developments in microarray data analysis. An ad hoc integration of phylogenetic uncertainty under these local clock models is also described. The performance and accuracy of the new methods are evaluated by reanalyzing three published data sets.

Conclusions

It is shown that the new maximum likelihood hybrid methods can perform better than penalized likelihood and almost as well as uncorrelated Bayesian models. However, the new methods still tend to underestimate the actual amount of rate change. This work demonstrates the difficulty of estimating divergence times using local molecular clocks.

Predicting evolutionary patterns of mammalian teeth from development

One motivation in the study of development is the discovery of mechanisms that may guide evolutionary change. Here we report how development governs relative size and number of cheek teeth, or molars, in the mouse. We constructed an inhibitory cascade model by experimentally uncovering the activator-inhibitor logic of sequential tooth development. The inhibitory cascade acts as a ratchet that determines molar size differences along the jaw, one effect being that the second molar always makes up one-third of total molar area. By using a macroevolutionary test, we demonstrate the success of the model in predicting dentition patterns found among murine rodent species with various diets, thereby providing an example of ecologically driven evolution along a developmentally favoured trajectory. In general, our work demonstrates how to construct and test developmental rules with evolutionary predictability in natural systems.

Mechanistic approaches to the study of evolution: the functional synthesis

An emerging synthesis of evolutionary biology and experimental molecular biology is providing much stronger and deeper inferences about the dynamics and mechanisms of evolution than were possible in the past. The new approach combines statistical analyses of gene sequences with manipulative molecular experiments to reveal how ancient mutations altered biochemical processes and produced novel phenotypes. This functional synthesis has set the stage for major advances in our understanding of fundamental questions in evolutionary biology. Here we describe this emerging approach, highlight important new insights that it has made possible, and suggest future directions for the field.

Functional conservation of Rel binding sites in drosophilid genomes

Evolutionary constraints on gene regulatory elements are poorly understood: Little is known about how the strength of transcription factor binding correlates with DNA sequence conservation, and whether transcription factor binding sites can evolve rapidly while retaining their function. Here we use the model of the NFKB/Rel-dependent gene regulation in divergent Drosophila species to examine the hypothesis that the functional properties of authentic transcription factor binding sites are under stronger evolutionary constraints than the genomic background. Using molecular modeling we compare tertiary structures of the Drosophila Rel family proteins Dorsal, Dif, and Relish and demonstrate that their DNA-binding and protein dimerization domains undergo distinct rates of evolution. The accumulated amino acid changes, however, are unlikely to affect DNA sequence recognition and affinity. We employ our recently developed microarray-based experimental platform and principal coordinates statistical analysis to quantitatively and systematically profile DNA binding affinities of three Drosophila Rel proteins to 10,368 variants of the NFKB recognition sequences. We then correlate the evolutionary divergence of gene regulatory regions with differences in DNA binding affinities. Genome-wide analyses reveal a significant increase in the number of conserved Rel binding sites in promoters of developmental and immune genes. Significantly, the affinity of Rel proteins to these sites was higher than to less conserved sites and was maintained by the conservation of the DNA binding site sequence (static conservation) or in some cases despite significantly diverged sequences (dynamic conservation). We discuss how two types of conservation may contribute to the stabilization and optimization of a functional gene regulatory code in evolution.

Genome clashes in hybrids: insights from gene expression

In interspecific hybrids, novel phenotypes often emerge from the interaction of two divergent genomes. Interactions between the two transcriptional networks are assumed to contribute to these unpredicted new phenotypes by inducing novel patterns of gene expression. Here we provide a review of the recent literature on the accumulation of regulatory incompatibilities. We review specific examples of regulatory incompatibilities reported at particular loci as well as genome-scale surveys of gene expression in interspecific hybrids. Finally, we consider and preview novel technologies that could help decipher how divergent transcriptional networks interact in hybrids between species.

Morphological evolution through multiple cis-regulatory mutations at a single gene

One central, and yet unsolved, question in evolutionary biology is the relationship between the genetic variants segregating within species and the causes of morphological differences between species. The classic neo-darwinian view postulates that species differences result from the accumulation of small-effect changes at multiple loci. However, many examples support the possible role of larger abrupt changes in the expression of developmental genes in morphological evolution. Although this evidence might be considered a challenge to a neo-darwinian micromutationist view of evolution, there are currently few examples of the actual genes causing morphological differences between species. Here we examine the genetic basis of a trichome pattern difference between Drosophila species, previously shown to result from the evolution of a single gene, shavenbaby (svb), probably through cis-regulatory changes. We first identified three distinct svb enhancers from D. melanogaster driving reporter gene expression in partly overlapping patterns that together recapitulate endogenous svb expression. All three homologous enhancers from D. sechellia drive expression in modified patterns, in a direction consistent with the evolved svb expression pattern. To test the influence of these enhancers on the actual phenotypic difference, we conducted interspecific genetic mapping at a resolution sufficient to recover multiple intragenic recombinants. This functional analysis revealed that independent genetic regions upstream of svb that overlap the three identified enhancers are collectively required to generate the D. sechellia trichome pattern. Our results demonstrate that the accumulation of multiple small-effect changes at a single locus underlies the evolution of a morphological difference between species. These data support the view that alleles of large effect that distinguish species may sometimes reflect the accumulation of multiple mutations of small effect at select genes.

Sex-specific expression of a HOX gene associated with rapid morphological evolution

Animal diversity is shaped by the origin and diversification of new morphological structures. Many examples of evolutionary innovations are provided by male-specific traits involved in mating and sexual selection. The origin of new sex-specific characters requires the evolution of new regulatory interactions between sex-determining genes and genes that control spatial patterning and cell differentiation. Here, we show that sex-specific regulation of the HOX gene Sex combs reduced (Scr) is associated with the origin and evolution of the Drosophila sex comb - a novel and rapidly diversifying male-specific organ. In species that primitively lack sex combs, Scr expression shows little spatial modulation, whereas in species that have sex combs, Scr is upregulated in the presumptive sex comb region and is frequently sexually dimorphic. Phylogenetic analysis shows that sex-specific regulation of Scr has been gained and lost multiple times in Drosophila evolution and correlates with convergent origin of similar sex comb morphologies in several independent lineages. Some of these transitions occurred on microevolutionary timescales, indicating that HOX gene expression can evolve with surprising ease. This is the first example of a sex-specific regulation of a HOX gene contributing to the development and evolution of a secondary sexual trait.

Cis-regulation of the amphioxus engrailed gene: Insights into evolution of a muscle-specific enhancer

To gain insights into the relation between evolution of cis-regulatory DNA and evolution of gene function, we identified tissue-specific enhancers of the engrailed gene of the basal chordate amphioxus (Branchiostoma floridae) and compared their ability to direct expression in both amphioxus and its nearest chordate relative, the tunicate Ciona intestinalis. In amphioxus embryos, the native engrailed gene is expressed in three domains - the eight most anterior somites, a few cells in the central nervous system (CNS) and a few ectodermal cells. In contrast, in C. intestinalis, in which muscle development is highly divergent, engrailed expression is limited to the CNS. To characterize the tissue-specific enhancers of amphioxus engrailed, we first showed that 7.8kb of upstream DNA of amphioxus engrailed directs expression to all three domains in amphioxus that express the native gene. We then identified the amphioxus engrailed muscle-specific enhancer as the 1.2kb region of upstream DNA with the highest sequence identity to the mouse en-2 jaw muscle enhancer. This amphioxus enhancer directed expression to both the somites in amphioxus and to the larval muscles in C. intestinalis. These results show that even though expression of the native engrailed has apparently been lost in developing C. intestinalis muscles, they express the transcription factors necessary to activate transcription from the amphioxus engrailed enhancer, suggesting that gene networks may not be completely disrupted if an individual component is lost.

The cis-regulatory logic of the mammalian photoreceptor transcriptional network

The photoreceptor cells of the retina are subject to a greater number of genetic diseases than any other cell type in the human body. The majority of more than 120 cloned human blindness genes are highly expressed in photoreceptors. In order to establish an integrative framework in which to understand these diseases, we have undertaken an experimental and computational analysis of the network controlled by the mammalian photoreceptor transcription factors, Crx, Nrl, and Nr2e3. Using microarray and in situ hybridization datasets we have produced a model of this network which contains over 600 genes, including numerous retinal disease loci as well as previously uncharacterized photoreceptor transcription factors. To elucidate the connectivity of this network, we devised a computational algorithm to identify the photoreceptor-specific cis-regulatory elements (CREs) mediating the interactions between these transcription factors and their target genes. In vivo validation of our computational predictions resulted in the discovery of 19 novel photoreceptor-specific CREs near retinal disease genes. Examination of these CREs permitted the definition of a simple cis-regulatory grammar rule associated with high-level expression. To test the generality of this rule, we used an expanded form of it as a selection filter to evolve photoreceptor CREs from random DNA sequences in silico. When fused to fluorescent reporters, these evolved CREs drove strong, photoreceptor-specific expression in vivo. This study represents the first systematic identification and in vivo validation of CREs in a mammalian neuronal cell type and lays the groundwork for a systems biology of photoreceptor transcriptional regulation.

The locus of evolution: evo devo and the genetics of adaptation

An important tenet of evolutionary developmental biology ("evo devo") is that adaptive mutations affecting morphology are more likely to occur in the cis-regulatory regions than in the protein-coding regions of genes. This argument rests on two claims: (1) the modular nature of cis-regulatory elements largely frees them from deleterious pleiotropic effects, and (2) a growing body of empirical evidence appears to support the predominant role of gene regulatory change in adaptation, especially morphological adaptation. Here we discuss and critique these assertions. We first show that there is no theoretical or empirical basis for the evo devo contention that adaptations involving morphology evolve by genetic mechanisms different from those involving physiology and other traits. In addition, some forms of protein evolution can avoid the negative consequences of pleiotropy, most notably via gene duplication. In light of evo devo claims, we then examine the substantial data on the genetic basis of adaptation from both genome-wide surveys and single-locus studies. Genomic studies lend little support to the cis-regulatory theory: many of these have detected adaptation in protein-coding regions, including transcription factors, whereas few have examined regulatory regions. Turning to single-locus studies, we note that the most widely cited examples of adaptive cis-regulatory mutations focus on trait loss rather than gain, and none have yet pinpointed an evolved regulatory site. In contrast, there are many studies that have both identified structural mutations and functionally verified their contribution to adaptation and speciation. Neither the theoretical arguments nor the data from nature, then, support the claim for a predominance of cis-regulatory mutations in evolution. Although this claim may be true, it is at best premature. Adaptation and speciation probably proceed through a combination of cis-regulatory and structural mutations, with a substantial contribution of the latter.

Patterns on the insect wing

The evolution of wings and the adaptive advantage they provide have allowed insects to become one of the most evolutionarily successful groups on earth. The incredible diversity of their shape, size, and color patterns is a direct reflection of the important role wings have played in the radiation of insects. In this review, we highlight recent studies on both butterflies and Drosophila that have begun to uncover the types of genetic variations and developmental mechanisms that control diversity in wing color patterns. In butterflies, these analyses are now possible because of the recent development of a suite of genomic and functional tools, such as detailed linkage maps and transgenesis. In one such study, extensive linkage mapping in Heliconius butterflies has shown that surprisingly few, and potentially homologous, loci are responsible for several major pattern variations on the wings of these butterflies. Parallel work on a clade of Drosophila has uncovered how cis-regulatory changes of the same gene correlate with the repeated gain and loss of pigmented wing spots. Collectively, our understanding of formation and evolution of color pattern in insect wings is rapidly advancing because of these recent breakthroughs in several different fields.

Canalization of segmentation and its evolution in Drosophila

Segmentation in Drosophila embryogenesis occurs through a hierarchical cascade of regulatory gene expression driven by the establishment of a diffusion-mediated morphogen gradient. Here, we investigate the response of this pattern formation process to genetic variation and evolution in egg size. Specifically, we ask whether spatial localization of gap genes Kruppel (Kr) and giant (gt) and the pair-rule gene even-skipped (eve) during cellularization is robust to genetic variation in embryo length in three Drosophila melanogaster isolines and two closely related species. We identified two wild-derived strains of D. melanogaster whose eggs differ by approximately 25% in length when reared under identical conditions. These two lines, a D. melanogaster laboratory stock (w1118), and offspring from crosses between the lines all exhibit precise scaling in the placement of gap and pair-rule gene expression along the anterior-posterior axis in relation to embryo length. Genetic analysis indicates that this scaling is maternally controlled. Maternal regulation of scaling must be required for consistent localization of segmentation gene expression because embryo size, a genetically variable and adaptive trait, is maternally inherited. We also investigated spatial scaling between these D. melanogaster lines and single lines of Drosophila sechellia and Drosophila simulans, the latter two differing by approximately 25% in egg length. In contrast to the robust scaling we observed within species, localization of gene expression relative to embryo length differs significantly between the three species. Thus, the developmental mechanism that assures robust scaling within a species does not prevent rapid evolution between species.

Molecular phylogeny of the Drosophila obscura species group, with emphasis on the Old World species

Background

Species of the Drosophila obscura species group (e.g., D. pseudoobscura, D. subobscura) have served as favorable models in evolutionary studies since the 1930's. Despite numbers of studies conducted with varied types of data, the basal phylogeny in this group is still controversial, presumably owing to not only the hypothetical 'rapid radiation' history of this group, but also limited taxon sampling from the Old World (esp. the Oriental and Afrotropical regions). Here we reconstruct the phylogeny of this group by using sequence data from 6 loci of 21 species (including 16 Old World ones) covering all the 6 subgroups of this group, estimate the divergence times among lineages, and statistically test the 'rapid radiation' hypothesis.

Results

Phylogenetic analyses indicate that each of the subobscura, sinobscura, affinis, and pseudoobscura subgroups is monophyletic. The subobscura and microlabis subgroups form the basal clade in the obscura group. Partial species of the obscura subgroup (the D. ambigua/D. obscura/D. tristis triad plus the D. subsilvestris/D. dianensis pair) forms a monophyletic group which appears to be most closely related to the sinobscura subgroup. The remaining basal relationships in the obscura group are not resolved by the present study. Divergence times on a ML tree based on mtDNA data are estimated with a calibration of 30–35 Mya for the divergence between the obscura and melanogaster groups. The result suggests that at least half of the current major lineages of the obscura group originated by the mid-Miocene time (~15 Mya), a time of the last developing and fragmentation of the temperate forest in North Hemisphere.

Conclusion

The obscura group began to diversify rapidly before invading into the New World. The subobscura and microlabis subgroups form the basal clade in this group. The obscura subgroup is paraphyletic. Partial members of this subgroup (D. ambigua, D. obscura, D. tristis, D. subsilvestris, and D. dianensis) form a monophyletic group which appears to be most closely related to the sinobscura subgroup.

Phylogenetic incongruence in the Drosophila melanogaster species group

Drosophila melanogaster and its close relatives are used extensively in comparative biology. Despite the importance of phylogenetic information for such studies, relationships between some melanogaster species group members are unclear due to conflicting phylogenetic signals at different loci. In this study, we use twelve nuclear loci (eleven coding and one non-coding) to assess the degree of phylogenetic incongruence in this model system. We focus on two nodes: (1) the node joining the Drosophila erecta-Drosophila orena, Drosophila melanogaster-Drosophila simulans, and Drosophila yakuba-Drosophila teissieri lineages, and (2) the node joining the lineages leading to the melanogaster, takahashii, and eugracilis subgroups. We find limited evidence for incongruence at the first node; our data, as well as those of several previous studies, strongly support monophyly of a clade consisting of D. erecta-D. orena and D. yakuba-D. teissieri. By contrast, using likelihood based tests of congruence, we find robust evidence for topological incongruence at the second node. Different loci support different relationships among the melanogaster, takahashii, and eugracilis subgroups, and the observed incongruence is not easily attributable to homoplasy, non-equilibrium base composition, or positive selection on a subset of loci. We argue that lineage sorting in the common ancestor of these three subgroups is the most plausible explanation for our observations. Such lineage sorting may lead to biased estimation of tree topology and evolutionary rates, and may confound inferences of positive selection.

A database of wing diversity in the Hawaiian Drosophila

BACKGROUND

Within genus Drosophila, the endemic Hawaiian species offer some of the most dramatic examples of morphological and behavioral evolution. The advent of the Drosophila grimshawi genome sequence permits genes of interest to be readily cloned from any of the hundreds of species of Hawaiian Drosophila, offering a powerful comparative approach to defining molecular mechanisms of species evolution. A key step in this process is to survey the Hawaiian flies for characters whose variation can be associated with specific candidate genes. The wings provide an attractive target for such studies: Wings are essentially two dimensional, and genes controlling wing shape, vein specification, pigment production, and pigment pattern evolution have all been identified in Drosophila.

METHODOLOGY/PRINCIPAL FINDINGS

We present a photographic database of over 180 mounted, adult wings from 73 species of Hawaiian Drosophila. The image collection, available at FlyBase.org, includes 53 of the 112 known species of "picture wing" Drosophila, and several species from each of the other major Hawaiian groups, including the modified mouthparts, modified tarsus, antopocerus, and haleakalae (fungus feeder) groups. Direct image comparisons show that major wing shape changes can occur even between closely related species, and that pigment pattern elements can vary independently of each other. Among the 30 species closest to grimshawi, diverse visual effects are achieved by altering a basic pattern of seven wing spots. Finally, we document major pattern variations within species, which appear to result from reduced diffusion of pigment precursors through the wing blade.

CONCLUSIONS/SIGNIFICANCE

The database highlights the striking variation in size, shape, venation, and pigmentation in Hawaiian Drosophila, despite their generally low levels of DNA sequence divergence. In several independent lineages, highly complex patterns are derived from simple ones. These lineages offer a promising model system to study the evolution of complexity.

The stars and stripes of animal bodies: evolution of regulatory elements mediating pigment and bristle patterns in Drosophila

Evolution has generated enormous morphological diversity in animals and one of the genetic processes that might have contributed to this is evolution of the cis-regulatory sequences responsible for the temporal and spatial expression of genes regulating embryonic development. This could be particularly relevant to pleiotropic genes with multiple independently acting regulatory modules. Loss or gain of modules enables altered expression without loss of other functions. Here I focus on recent studies correlating differences in morphological traits between related species of Drosophila to changes in cis-regulatory sequences. They show that ancestral regulatory modules have evolved to mediate different transcriptional outputs and suggest that evolution of cis-regulatory sequences might reflect a general mechanism driving evolutionary change.

Emerging principles of regulatory evolution

Understanding the genetic and molecular mechanisms governing the evolution of morphology is a major challenge in biology. Because most animals share a conserved repertoire of body-building and -patterning genes, morphological diversity appears to evolve primarily through changes in the deployment of these genes during development. The complex expression patterns of developmentally regulated genes are typically controlled by numerous independent cis-regulatory elements (CREs). It has been proposed that morphological evolution relies predominantly on changes in the architecture of gene regulatory networks and in particular on functional changes within CREs. Here, we discuss recent experimental studies that support this hypothesis and reveal some unanticipated features of how regulatory evolution occurs. From this growing body of evidence, we identify three key operating principles underlying regulatory evolution, that is, how regulatory evolution: (i) uses available genetic components in the form of preexisting and active transcription factors and CREs to generate novelty; (ii) minimizes the penalty to overall fitness by introducing discrete changes in gene expression; and (iii) allows interactions to arise among any transcription factor and downstream CRE. These principles endow regulatory evolution with a vast creative potential that accounts for both relatively modest morphological differences among closely related species and more profound anatomical divergences among groups at higher taxonomical levels.

Variable gene expression in eukaryotes: a network perspective

Changes in gene expression underlie phenotypic plasticity, variation within species, and phenotypic divergence between species. These expression differences arise from modulation of regulatory networks. To understand the source of expression differences, networks of interactions among genes and gene products that orchestrate gene expression must be considered. Here I review the basic structure of eukaryotic regulatory networks and discuss selected case studies that provide insight into how these networks are altered to create expression differences within and between species.

The evolutionary significance of cis-regulatory mutations

For decades, evolutionary biologists have argued that changes in cis-regulatory sequences constitute an important part of the genetic basis for adaptation. Although originally based on first principles, this claim is now empirically well supported: numerous studies have identified cis-regulatory mutations with functionally significant consequences for morphology, physiology and behaviour. The focus has now shifted to considering whether cis-regulatory and coding mutations make qualitatively different contributions to phenotypic evolution. Cases in which parallel mutations have produced parallel trait modifications in particular suggest that some phenotypic changes are more likely to result from cis-regulatory mutations than from coding mutations.

Comparative functional analysis of the lac operons in Streptococcus pyogenes

Having no known environmental reservoir, Streptococcus pyogenes, a bacterium responsible for a wider variety of human diseases than any other bacterial species, must rely on its host for metabolic substrates. Although a streptococcal aldolase, LacD.1, has been adapted to virulence gene regulation, both LacD.1 and a paralogous protein, LacD.2, are predicted to function in the tagatose 6-phosphate pathway for lactose and galactose utilization. In order to gain insight into the mechanism of the LacD.1 regulatory pathway and the role of genome context in the emergence of LacD.1's novel regulatory functions, we compared the function and regulation of the Lac.1 and Lac.2 loci. The Lac.1 operon is not inducible, and regulation by LacD.1 is independent of a functional tagatose 6-phosphate pathway and enhanced by the conserved truncation of upstream Lac.1 genes. In contrast, Lac.2 expression is sensitive to environmental carbohydrates, and LacD.2, not LacD.1, contributes to growth on galactose. Thus, we conclude that the Lac.1 locus has been specialized to participate in regulation, leaving efficient utilization of carbohydrate sources to the Lac.2 locus. The adaptation of LacD for transcription regulation may be an underappreciated strategy among prokaryotes, as homologues of this multifaceted enzyme are present in a broad range of species.

Two or four bristles: functional evolution of an enhancer of scute in Drosophilidae

Changes in cis-regulatory sequences are proposed to underlie much of morphological evolution. Yet, little is known about how such modifications translate into phenotypic differences. To address this problem, we focus on the dorsocentral bristles of Drosophilidae. In Drosophila melanogaster, development of these bristles depends on a cis-regulatory element, the dorsocentral enhancer, to activate scute in a cluster of cells from which two bristles on the posterior scutum arise. A few species however, such as D. quadrilineata, bear anterior dorsocentral bristles as well as posterior ones, a derived feature. This correlates with an anterior expansion of the scute expression domain. Here, we show that the D. quadrilineata enhancer has evolved, and is now active in more anterior regions. When used to rescue scute expression in transgenic D. melanogaster, the D. quadrilineata enhancer is able to induce anterior bristles. Importantly, these properties are not displayed by homologous enhancers from control species bearing only two posterior bristles. We also provide evidence that upstream regulation of the enhancer, by the GATA transcription factor Pannier, has been evolutionarily conserved. This work illustrates how, in the context of a conserved trans-regulatory landscape, evolutionary tinkering of pre-existing enhancers can modify gene expression patterns and contribute to morphological diversification.

Evolutionary change in function of the dorsocentral enhancer (DCE) of scute has resulted in altered bristle formation between two species of Drosophila.

Evolution of the ventral midline in insect embryos

The ventral midline is a source of signals that pattern the nerve cord of insect embryos. In dipterans such as the fruitfly Drosophila melanogaster (D. mel.) and the mosquito Anopheles gambiae (A. gam.), the midline is narrow and spans just 1-2 cells. However, in the honeybee, Apis mellifera (A. mel.), the ventral midline is broad and encompasses 5-6 cells. slit and other midline-patterning genes display a corresponding expansion in expression. Evidence is presented that this difference is due to divergent cis regulation of the single-minded (sim) gene, which encodes a bHLH-PAS transcription factor essential for midline differentiation. sim is regulated by a combination of Notch signaling and a Twist (Twi) activator gradient in D. mel., but it is activated solely by Twi in A. mel. We suggest that the Twi-only mode of regulation--and the broad ventral midline--represents the ancestral form of CNS patterning in Holometabolous insects.

The neural and genetic substrates of sexual behavior in Drosophila

fruitless (fru), originally identified with its mutant conferring male homosexuality, is a neural sex determination gene in Drosophila that produces sexually dimorphic sets of transcripts. In the nervous system, Fru is translated only in males. Fru proteins likely regulate the transcription of a set of downstream genes. The expression of Fru proteins is sufficient to induce male sexual behavior in females. A group of fru-expressing neurons called "mAL" neurons in the brain shows conspicuous sexual dimorphism. mAL is composed of 5 neurons in females and 30 neurons in males. It includes neurons with bilateral projections in males and contralateral projections in females. Terminal arborization patterns are also sexually dimorphic. These three characteristics are feminized in fru mutant males. The inactivation of cell death genes results in the production of additional mAL neurons that are of the male type in the female brain. This suggests that male-specific Fru inhibits mAL neuron death, leading to the formation of a male-specific neural circuit that underlies male sexual behavior. Fru orchestrates a spectrum of downstream genes as a master control gene to establish the maleness of the brain.

Large-scale turnover of functional transcription factor binding sites in Drosophila

The gain and loss of functional transcription factor binding sites has been proposed as a major source of evolutionary change in cis-regulatory DNA and gene expression. We have developed an evolutionary model to study binding-site turnover that uses multiple sequence alignments to assess the evolutionary constraint on individual binding sites, and to map gain and loss events along a phylogenetic tree. We apply this model to study the evolutionary dynamics of binding sites of the Drosophila melanogaster transcription factor Zeste, using genome-wide in vivo (ChIP-chip) binding data to identify functional Zeste binding sites, and the genome sequences of D. melanogaster, D. simulans, D. erecta, and D. yakuba to study their evolution. We estimate that more than 5% of functional Zeste binding sites in D. melanogaster were gained along the D. melanogaster lineage or lost along one of the other lineages. We find that Zeste-bound regions have a reduced rate of binding-site loss and an increased rate of binding-site gain relative to flanking sequences. Finally, we show that binding-site gains and losses are asymmetrically distributed with respect to D. melanogaster, consistent with lineage-specific acquisition and loss of Zeste-responsive regulatory elements.

Heliconius wing patterns: an evo-devo model for understanding phenotypic diversity

Evolutionary Developmental Biology aims for a mechanistic understanding of phenotypic diversity, and present knowledge is largely based on gene expression and interaction patterns from a small number of well-known model organisms. However, our understanding of biological diversification depends on our ability to pinpoint the causes of natural variation at a micro-evolutionary level, and therefore requires the isolation of genetic and developmental variation in a controlled genetic background. The colour patterns of Heliconius butterflies (Nymphalidae: Heliconiinae) provide a rich suite of naturally occurring variants with striking phenotypic diversity and multiple taxonomic levels of variation. Diversification in the genus is well known for its dramatic colour-pattern divergence between races or closely related species, and for Müllerian mimicry convergence between distantly related species, providing a unique system to study the development basis of colour-pattern evolution. A long history of genetic studies has showed that pattern variation is based on allelic combinations at a surprisingly small number of loci, and recent developmental evidence suggests that pattern development in Heliconius is different from the eyespot determination of other butterflies. Fine-scale genetic mapping studies have shown that a shared toolkit of genes is used to produce both convergent and divergent phenotypes. These exciting results and the development of new genomic resources make Heliconius a very promising evo-devo model for the study of adaptive change.

Parallel genetic origins of pelvic reduction in vertebrates

Despite longstanding interest in parallel evolution, little is known about the genes that control similar traits in different lineages of vertebrates. Pelvic reduction in stickleback fish (family Gasterosteidae) provides a striking example of parallel evolution in a genetically tractable system. Previous studies suggest that cis-acting regulatory changes at the Pitx1 locus control pelvic reduction in a population of threespine sticklebacks (Gasterosteus aculeatus). In this study, progeny from intergeneric crosses between pelvic-reduced threespine and ninespine (Pungitius pungitius) sticklebacks also showed severe pelvic reduction, implicating a similar genetic origin for this trait in both genera. Comparative sequencing studies in complete and pelvic-reduced Pungitius revealed no differences in the Pitx1 coding sequences, but Pitx1 expression was absent from the prospective pelvic region of larvae from pelvic-reduced parents. A much more phylogenetically distant example of pelvic reduction, loss of hindlimbs in manatees, shows a similar left-right size bias that is a morphological signature of Pitx1-mediated pelvic reduction in both sticklebacks and mice. These multiple lines of evidence suggest that changes in Pitx1 may represent a key mechanism of morphological evolution in multiple populations, species, and genera of sticklebacks, as well as in distantly related vertebrate lineages.

Regulation of body pigmentation by the Abdominal-B Hox protein and its gain and loss in Drosophila evolution

Hox genes have been implicated in the evolution of many animal body patterns, but the molecular events underlying trait modification have not been elucidated. Pigmentation of the posterior male abdomen is a recently acquired trait in the Drosophila melanogaster lineage. Here, we show that the Abdominal-B (ABD-B) Hox protein directly activates expression of the yellow pigmentation gene in posterior segments. ABD-B regulation of pigmentation evolved through the gain of ABD-B binding sites in a specific cis-regulatory element of the yellow gene of a common ancestor of sexually dimorphic species. Within the melanogaster species group, male-specific pigmentation has subsequently been lost by at least three different mechanisms, including the mutational inactivation of a key ABD-B binding site in one lineage. These results demonstrate how Hox regulation of traits and target genes is gained and lost at the species level and have general implications for the evolution of body form at higher taxonomic levels.

Nematode Development: New Tricks for Old Genes

The evolution of novel features requires evolutionary change in gene function. New work in nematodes shows how conserved genes can nonetheless contribute to innovative traits.