Developmental system drift and flexibility in evolutionary trajectories

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.

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.

Gene Duplication and the Evolution of Vertebrate Skeletal Mineralization

The mineralized skeleton is a critical innovation that evolved early in vertebrate history. The tissues found in dermal skeletons of ancient vertebrates are similar to the dental tissues of modern vertebrates; both consist of a highly mineralized surface hard tissue, enamel or enameloid, more resilient body dentin, and basal bone. Many proteins regulating mineralization of these tissues are evolutionarily related and form the secretory calcium-binding phosphoprotein (SCPP) family. We hypothesize here the duplication histories of SCPP genes and their common ancestors, SPARC and SPARCL1. At around the same time that Paleozoic jawless vertebrates first evolved mineralized skeleton, SPARCL1 arose from SPARC by whole genome duplication. Then both before and after the split of ray-finned fish and lobe-finned fish, tandem gene duplication created two types of SCPP genes, each residing on the opposite side of SPARCL1. One type was subsequently used in surface tissue and the other in body tissue. In tetrapods, these two types of SCPP genes were separated by intrachromosomal rearrangement. While new SCPP genes arose by duplication, some old genes were eliminated from the genome. As a consequence, phenogenetic drift occurred: while mineralized skeleton is maintained by natural selection, the underlying genetic basis has changed.

New candidate genes for sex-comb divergence between Drosophila mauritiana and Drosophila simulans

A large-effect QTL for divergence in sex-comb tooth number between Drosophila simulans and D. mauritiana was previously mapped to 73A-84AB. Here we identify genes that are likely contributors to this divergence. We first improved the mapping resolution in the 73A-84AB region using 12 introgression lines and 62 recombinant nearly isogenic lines. To further narrow the list of candidate genes, we assayed leg-specific expression and identified genes with transcript-level evolution consistent with a potential role in sex-comb divergence. Sex combs are formed on the prothoracic (front) legs, but not on the mesothoracic (middle) legs of Drosophila males. We extracted RNA from the prothoracic and mesothoracic pupal legs of two species to determine which of the genes expressed differently between leg types were also divergent for gene expression. Two good functional candidate genes, Scr and dsx, are located in one of our fine-scale QTL regions. In addition, three previously uncharacterized genes (CG15186, CG2016, and CG2791) emerged as new candidates. These genes are located in regions strongly associated with sex-comb tooth number differences and are expressed differently between leg tissues and between species. Further supporting the potential involvement of these genes in sex-comb divergence, we found a significant difference in sex-comb tooth number between co-isogenic D. melanogaster lines with and without P-element insertions at CG2791.

The frailty of adaptive hypotheses for the origins of organismal complexity

The vast majority of biologists engaged in evolutionary studies interpret virtually every aspect of biodiversity in adaptive terms. This narrow view of evolution has become untenable in light of recent observations from genomic sequencing and population-genetic theory. Numerous aspects of genomic architecture, gene structure, and developmental pathways are difficult to explain without invoking the nonadaptive forces of genetic drift and mutation. In addition, emergent biological features such as complexity, modularity, and evolvability, all of which are current targets of considerable speculation, may be nothing more than indirect by-products of processes operating at lower levels of organization. These issues are examined in the context of the view that the origins of many aspects of biological diversity, from gene-structural embellishments to novelties at the phenotypic level, have roots in nonadaptive processes, with the population-genetic environment imposing strong directionality on the paths that are open to evolutionary exploitation.

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.

Evolution and Development: Anchors away!

Developmental mechanisms can evolve even when the trait they produce does not, and the nematode vulva has become a model organ for detecting such "developmental system drift". A new study reveals what may be the very earliest stages of this process by experimentally modifying key vulval signaling pathways in different species of Caenorhabditis, and carefully quantifying the results.

Variation and constraint in plant evolution and development

The goal of this short review is to consider the interrelated phenomena of phenotypic variation and genetic constraint with respect to plant diversity. The unique aspects of plants, including sessile habit, modular growth and diverse developmental programs expressed at the phytomer level, merit a specific examination of the genetic basis of their phenotypic variation, and how they experience and escape genetic constraint. Numerous QTL studies with wild and domesticated plants reveal that most phenotypic traits are polygenic but vary in the number and effect of the loci contributing, from a few loci of large effects to many with small effects. Further, somatic mutations, developmental plasticity and epigenetic variation, especially gene methylation, can contribute to increases in phenotypic variation. The flip side of these processes, genetic constraint, can similarly be the result of many factors, including pleiotropy, canalization and genetic redundancy. Genetic constraint is not only a mechanism to prevent change, however, it can also serve to direct evolution along certain paths. Ultimately, genetic constraint often comes full circle and is released through events such as hybridization, genome duplication and epigenetic remodeling. We are just beginning to understand how these processes can operate simultaneously during the evolution of ecologically important traits in plants.

Lessons from leeches: a call for DNA barcoding in the lab

Many evolution of development labs study organisms that must be periodically collected from the wild. Whenever this is the case, there is the risk that different field collections will recover genetically different strains or cryptic species. Ignoring this potential for genetic variation may introduce an uncontrolled source of experimental variability, leading to confusion or misinterpretation of the results. Leeches in the genus Helobdella have been a workhorse of annelid developmental biology for 30 years. Nearly all early Helobdella research was based on a single isolate, but in recent years isolates from multiple field collections and multiple sites across the country have been used. To assess the genetic distinctness of different isolates, we obtained specimens from most Helobdella laboratory cultures currently or recently in use and from some of their source field sites. From these samples, we sequenced part of the mitochondrial gene cytochrome oxidase I (COI). Sequence divergences and phylogenetic analyses reveal that, collectively, the Helobdella development community has worked on five distinct species from two major clades. Morphologically similar isolates that were thought to represent the same species (H. robusta) actually represent three species, two of which coexist at the same locality. Another isolate represents part of a species complex (the "H. triserialis" complex), and yet another is an invasive species (H. europaea). We caution researchers similarly working on multiple wild-collected isolates to preserve voucher specimens and to obtain from these a molecular "barcode," such as a COI gene sequence, to reveal genetic variation in animals used for research.

Compensatory vs. pseudocompensatory evolution in molecular and developmental interactions

The evolution of molecules, developmental circuits, and new species are all characterized by the accumulation of incompatibilities between ancestors and descendants. When specific interactions between components are necessary at any of these levels, this requires compensatory coevolution. Theoretical treatments of compensatory evolution that only consider the endpoints predict that it should be rare because intermediate states are deleterious. However, empirical data suggest that compensatory evolution is common at all levels of molecular interaction. A general solution to this paradox is provided by plausible neutral or nearly neutral intermediates that possess informational redundancy. These intermediates provide an evolutionary path between coadapted allelic combinations. Although they allow incompatible end points to evolve, at no point was a deleterious mutation ever in need of compensation. As a result, what appears to be compensatory evolution may often actually be "pseudocompensatory." Both theoretical and empirical studies indicate that pseudocompensation can speed the evolution of intergenic incompatibility, especially when driven by adaptation. However, under strong stabilizing selection the rate of pseudocompensatory evolution is still significant. Important examples of this process at work discussed here include the evolution of rRNA secondary structures, intra- and inter-protein interactions, and developmental genetic pathways. Future empirical work in this area should focus on comparing the details of intra- and intergenic interactions in closely related organisms.

Comparative gene expression of steroidogenic factor 1 in Chrysemys picta and Apalone mutica turtles with temperature-dependent and genotypic sex determination

Characterizing the molecular network underlying temperature-dependent (TSD) and genotypic (GSD) sex determination, including patterns across closely related taxa, is crucial to elucidate the still enigmatic evolution of sex determining mechanisms in vertebrates. Here we examined the expression of an important gene for sexual differentiation common to both systems, Sf1, at male- and female-producing temperatures, in TSD (Chrysemys picta) and GSD turtles (Apalone mutica). We tested the hypotheses that Sf1 expression responds to temperature consistently across TSD turtles but is unaffected in GSD turtles, and that this differential expression starts no earlier than the onset of the thermosensitive period (TSP). As expected, Sf1 expression was thermally insensitive in A. mutica (GSD). Although Sf1 exhibited a differential expression by temperature in C. picta, the expression pattern differed from other TSD turtles (Trachemys scripta), perhaps reflecting divergence of the gene regulatory networks underlying sex determination over evolutionary time. Most notably, Sf1 was differentially expressed in C. picta (significantly higher at the male-producing temperature) before the onset of the TSP, implying that in TSD taxa significant thermal effects may occur early in development. This result may reconcile field observations where temperatures experienced prior to the TSP have an effect on sex ratios, thus challenging traditional TSP models. Importantly, the molecular factors that render TSD mechanisms thermosensitive remain unknown, and potential candidates are genes that express differentially before the onset of the TSP (genes shaping or opening the TSP-window rather those acting once the TSP window has opened). Therefore, our findings make Sf1 one such potential candidate.

Genome-wide patterns of expression in Drosophila pure species and hybrid males. II. Examination of multiple-species hybridizations, platforms, and life cycle stages

Species often produce sterile hybrids early in their evolutionary divergence, and some evidence suggests that hybrid sterility may be associated with deviations or disruptions in gene expression. In support of this idea, many studies have shown that a high proportion of male-biased genes are underexpressed, compared with non-sex-biased genes, in sterile F1 male hybrids of Drosophila species. In this study, we examined and compared patterns of misexpression in sterile F1 male hybrids of Drosophila simulans and 2 of its sibling species, Drosophila mauritiana and Drosophila sechellia, at both the larval and adult life stages. We analyzed hybrids using both commercial Drosophila melanogaster microarrays and arrays we developed from reverse transcriptase-polymerase chain reactions of spermatogenesis and reproduction-related transcripts from these species (sperm array). Although the majority of misexpressed transcripts were underexpressed, a disproportionate number of the overexpressed transcripts were located on the X chromosome. We detected a high overlap in the genes misexpressed between the 2 species pairs, and our sperm array was better at detecting such misexpression than the D. melanogaster array, suggesting possible weaknesses in the use of an array designed from another species. We found only minimal misexpression in the larval samples with the sperm array, suggesting that disruptions in spermatogenesis occur after this life stage. Further study of these misexpressed loci may allow us to identify precisely where disruptions in the spermatogenesis pathway occur.

Developmental constraints vs. variational properties: How pattern formation can help to understand evolution and development

This article suggests that apparent disagreements between the concept of developmental constraints and neo-Darwinian views on morphological evolution can disappear by using a different conceptualization of the interplay between development and selection. A theoretical framework based on current evolutionary and developmental biology and the concepts of variational properties, developmental patterns and developmental mechanisms is presented. In contrast with existing paradigms, the approach in this article is specifically developed to compare developmental mechanisms by the morphological variation they produce and the way in which their functioning can change due to genetic variation. A developmental mechanism is a gene network, which is able to produce patterns in space though the regulation of some cell behaviour (like signalling, mitosis, apoptosis, adhesion, etc.). The variational properties of a developmental mechanism are all the pattern transformations produced under different initial and environmental conditions or IS-mutations. IS-mutations are DNA changes that affect how two genes in a network interact, while T-mutations are mutations that affect the topology of the network itself. This article explains how this new framework allows predictions not only about how pattern formation affects variation, and thus phenotypic evolution, but also about how development evolves by replacement between pattern formation mechanisms. This article presents testable inferences about the evolution of the structure of development and the phenotype under different selective pressures. That is what kind of pattern formation mechanisms, in which relative temporal order, and which kind of phenotypic changes, are expected to be found in development.

Evolution of branched regulatory genetic pathways: directional selection on pleiotropic loci accelerates developmental system drift

Developmental systems are regulated by a web of interacting loci. One common and useful approach in studying the evolution of development is to focus on classes of interacting elements within these systems. Here, we use individual-based simulations to study the evolution of traits controlled by branched developmental pathways involving three loci, where one locus regulates two different traits. We examined the system under a variety of selective regimes. In the case where one branch was under stabilizing selection and the other under directional selection, we observed "developmental system drift": the trait under stabilizing selection showed little phenotypic change even though the loci underlying that trait showed considerable evolutionary divergence. This occurs because the pleiotropic locus responds to directional selection and compensatory mutants are then favored in the pathway under stabilizing selection. Though developmental system drift may be caused by other mechanisms, it seems likely that it is accelerated by the same underlying genetic mechanism as that producing the Dobzhansky-Muller incompatibilities that lead to speciation in both linear and branched pathways. We also discuss predictions of our model for developmental system drift and how different selective regimes affect probabilities of speciation in the branched pathway system.

Functional and evolutionary inference in gene networks: does topology matter?

The relationship between the topology of a biological network and its functional or evolutionary properties has attracted much recent interest. It has been suggested that most, if not all, biological networks are 'scale free.' That is, their connections follow power-law distributions, such that there are very few nodes with very many connections and vice versa. The number of target genes of known transcriptional regulators in the yeast, Saccharomyces cerevisiae, appears to follow such a distribution, as do other networks, such as the yeast network of protein-protein interactions. These findings have inspired attempts to draw biological inferences from general properties associated with scale-free network topology. One often cited general property is that, when compromised, highly connected nodes will tend to have a larger effect on network function than sparsely connected nodes. For example, more highly connected proteins are more likely to be lethal when knocked out. However, the correlation between lethality and connectivity is relatively weak, and some highly connected proteins can be removed without noticeable phenotypic effect. Similarly, network topology only weakly predicts the response of gene expression to environmental perturbations. Evolutionary simulations of gene-regulatory networks, presented here, suggest that such weak or non-existent correlations are to be expected, and are likely not due to inadequacy of experimental data. We argue that 'top-down' inferences of biological properties based on simple measures of network topology are of limited utility, and we present simulation results suggesting that much more detailed information about a gene's location in a regulatory network, as well as dynamic gene-expression data, are needed to make more meaningful functional and evolutionary predictions. Specifically, we find in our simulations that: (1) the relationship between a gene's connectivity and its fitness effect upon knockout depends on its equilibrium expression level; (2) correlation between connectivity and genetic variation is virtually non-existent, yet upon independent evolution of networks with identical topologies, some nodes exhibit consistently low or high polymorphism; and (3) certain genes show low polymorphism yet high divergence among independent evolutionary runs. This latter pattern is generally taken as a signature of positive selection, but in our simulations its cause is often neutral coevolution of regulatory inputs to the same gene.

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.

Gene regulation divergence is a major contributor to the evolution of Dobzhansky-Muller incompatibilities between species of Drosophila

The Dobzhansky-Muller model denotes incompatible gene interactions between diverging populations/species and is recognized as the basis of postzygotic reproductive isolation. Little is known about the molecular nature of such gene interactions. We have carried out comparative gene expression analyses in the testes of 3 closely related species of the Drosophila melanogaster subgroup and their hybrids (all of which are sterile). We show that in hybrids 1) a higher proportion of male-biased genes (i.e., genes with a higher level of expression in males) are underexpressed (or not expressed) compared with non-sex-biased genes, 2) the majority of the underexpressed genes appear to be under stabilizing selection by virtue of showing similar levels of expression in the parental species, and only a small proportion of genes show signs of directional selection, 3) very few of the misexpressed genes are shared between species pairs, suggesting that there may not be a "common" set of "speciation genes," and 4) expression of non-testes-specific genes is observed in the testes of interspecific hybrids, and the number of such genes is positively correlated with divergence time. These results suggest that gene regulation divergence of sex- and reproduction-related genes is a major contributor to the evolution of Dobzhansky-Muller incompatibilities between species of Drosophila.

Genetic flexibility in the convergent evolution of hermaphroditism in Caenorhabditis nematodes

The self-fertile hermaphrodites of C. elegans and C. briggsae evolved from female ancestors by acquiring limited spermatogenesis. Initiation of C. elegans hermaphrodite spermatogenesis requires germline translational repression of the female-promoting gene tra-2, which allows derepression of the three male-promoting fem genes. Cessation of hermaphrodite spermatogenesis requires fem-3 translational repression. We show that C. briggsae requires neither fem-2 nor fem-3 for hermaphrodite development, and that XO Cb-fem-2/3 animals are transformed into hermaphrodites, not females as in C. elegans. Exhaustive screens for Cb-tra-2 suppressors identified another 75 fem-like mutants, but all are self-fertile hermaphrodites rather than females. Control of hermaphrodite spermatogenesis therefore acts downstream of the fem genes in C. briggsae. The outwardly similar hermaphrodites of C. elegans and C. briggsae thus achieve self-fertility via intervention at different points in the core sex determination pathway. These findings are consistent with convergent evolution of hermaphroditism, which is marked by considerable developmental genetic flexibility.

Genetic basis for innovations in floral organ identity

Of the many innovations associated with the radiation of the angiosperms, the evolution of a petal identity program is among the best understood from a genetic standpoint. Although the existing data do indicate that similar genetic mechanisms control petal development across diverse taxa, there is also considerable evidence for variability in petal identity programs, likely due to a number of factors. These points are illustrated through a review of our current knowledge on the subject, integrating phylogenetic, morphological, and genetic studies. Comparative studies of petal identity highlight the complex nature of homology in plants and stand as a cautionary tale for the interpretation of gene expression data.

Phenogenetic drift in evolution: the changing genetic basis of vertebrate teeth

Vertebrate mineralized tissues are vital to the adaptive evolution of various traits. Among these traits is the tooth, which consists of two characteristic mineralized tissues, a highly mineralized surface layer (enamel in tetrapods and enameloid in fish) and a softer body (dentin), both supported by basal bone. However, enamel and enameloid are significantly different in development, and dentin shows many histological variations; hence their evolution has been intensively studied. Nevertheless, their genetic basis has been revealed only in tetrapods. We previously reported that many genes involved in tetrapod tissue mineralization arose from a common ancestor and constitute the secretory calcium-binding phosphoprotein (SCPP) gene family. Now we show that teleost fish also use many SCPPs for enameloid and dentin mineralization, but none of these directly corresponds to tetrapod SCPPs. This finding suggests that teleost and tetrapod SCPP genes have experienced independent parallel duplication histories. Thus, through phenogenetic drift, the tooth has remained a stable trait in jawed vertebrates, while evolving distinct genetic bases in teleosts and tetrapods. The characteristics of teleost SCPP genes and their expression domains in tooth development suggest the possibility that enameloid arose from dentin and enamel from enameloid more than once in vertebrate evolution. In fugu (puffer fish), expression of SCPP genes is also detected in an unusual beak-like structure that shelters numerous teeth. Their expression pattern suggests that the jaw consists of the dentin beak and supportive bone. These findings illustrate the complexity of the homology concept in understanding evolution, particularly the evolution of mineralized tissues.

Comparative Developmental Genetics and the Evolution of Arthropod Body Plans

The arthropods display a wide range of morphological diversity, varying tagmosis, as well as other aspects of the body plan, such as appendage and cuticular morphology. Here we review the roles of developmental regulatory genes in the evolution of arthropod morphology, with an emphasis on what is known from morphologically diverse species. Examination of tagmatic evolution reveals that these changes have been accompanied by changes in the expression patterns of Hox genes. In contrast, review of the modifications to wing morphology seen in insects shows that these body plan changes have generally favored alterations in downstream target genes. These and other examples are used to discuss the evolutionary implications of comparative developmental genetic data.

Hox gene function and interaction in the milkweed bug Oncopeltus fasciatus (Hemiptera)

Studies in genetic model organisms such as Drosophila have demonstrated that the homeotic complex (Hox) genes impart segmental identity during embryogenesis. Comparative studies in a wide range of other insect taxa have shown that the Hox genes are expressed in largely conserved domains along the anterior-posterior body axis, but whether they are performing the same functions in different insects is an open question. Most of the Hox genes have been studied functionally in only a few holometabolous insects that undergo metamorphosis. Thus, it is unclear how the Hox genes are functioning in the majority of direct-developing insects and other arthropods. To address this question, we used a combination of RNAi and in situ hybridization to reveal the expression, functions, and regulatory interactions of the Hox genes in the milkweed bug Oncopeltus fasciatus. Our results reveal many similarities and some interesting differences compared to Drosophila. We find that the gene Antennapedia is required for the identity of all three thoracic segments, while Ultrabithorax, abdominal-A and Abdominal-B cooperate to pattern the abdomen. The three abdominal genes exhibit posterior prevalence like in Drosophila, but apparently via some post-transcriptional mechanism. The functions of the head genes proboscipedia, Deformed, and Sex combs reduced were shown previously, and here we find that the complex temporal expression of pb in the labium is like that of other insects, but its regulatory relationship with Scr is unique. Overall, our data reveal that the evolution of insect Hox genes has included many small changes within general conservation of expression and function, and that the milkweed bug provides a useful model for understanding the roles of Hox genes in a direct-developing insect.

Compensatory cis-trans evolution and the dysregulation of gene expression in interspecific hybrids of Drosophila

Hybrids between species are often characterized by novel gene-expression patterns. A recent study on allele-specific gene expression in hybrids between species of Drosophila revealed cases in which cis- and trans-regulatory elements within species had coevolved in such a way that changes in cis-regulatory elements are compensated by changes in trans-regulatory elements. We hypothesized that such coevolution should often lead to gene misexpression in the hybrid. To test this hypothesis, we estimated allele-specific expression and overall expression levels for 31 genes in D. melanogaster, D. simulans, and their F1 hybrid. We found that 13 genes with cis-trans compensatory evolution are in fact misexpressed in the hybrid. These represent candidate genes whose dysregulation might be the consequence of coevolution of cis- and trans-regulatory elements within species. Using a mathematical model for the regulation of gene expression, we explored the conditions under which cis-trans compensatory evolution can lead to misexpression in interspecific hybrids.

The genetics of ray pattern variation in Caenorhabditis briggsae

BACKGROUND

How does intraspecific variation relate to macroevolutionary change in morphology? This question can be addressed in species in which derived characters are present but not fixed. In rhabditid nematodes, the arrangement of the nine bilateral pairs of peripheral sense organs (rays) in tails of males is often the most highly divergent character between species. The development of ray pattern involves inputs from hometic gene expression patterns, TGFbeta signalling, Wnt signalling, and other genetic pathways. In Caenorhabditis briggsae, strain-specific variation in ray pattern has provided an entree into the evolution of ray pattern. Some strains were fixed for a derived pattern. Other strains were more plastic and exhibited derived and ancestral patterns at equal frequencies.

RESULTS

Recombinant inbred lines (RILs) constructed from crosses between the variant C. briggsae AF16 and HK104 strains exhibited a wide range of phenotypes including some that were more extreme than either parental strain. Transgressive segregation was significantly associated with allelic variation in the C. briggsae homolog of abdominal B, Cb-egl-5. At least two genes that affected different elements of ray pattern, ray position and ray fusion, were linked to a second gene, mip-1. Consistent with this, the segregation of ray position and ray fusion phenotypes were only partially correlated in the RILs.

CONCLUSIONS

The evolution of ray pattern has involved allelic variation at multiple loci. Some of these loci impact the specification of ray identities and simultaneously affect multiple ray pattern elements. Others impact individual characters and are not constrained by covariance with other ray pattern elements. Among the genetic pathways that may be involved in ray pattern evolution is specification of anteroposterior positional information by homeotic genes.

Development of foregut and proboscis in the buccinid neogastropodNassarius mendicus: Evolutionary opportunity exploited by a developmental module

This article extends previous descriptions of foregut development and metamorphosis in neogastropods by providing data on the buccinid Nassarius mendicus, a species with a feeding larva. Histological sections showed that, like many other gastropods, the postmetamorphic buccal cavity and radular sac of N. mendicus differentiate during the larval stage from a ventral outpocketing of the distal larval esophagus. However, in N. mendicus the outpocketing also gives rise to the entire anterior esophagus and valve of Leiblein, suggesting that both these structures may be evolutionary derivatives of the gastropod buccal cavity. Scanning electron microscopy and three-dimensional reconstructions of section profiles revealed that the distal larval esophagus and larval mouth are completely destroyed at metamorphosis. The postmetamorphic mouth is formed as a new orifice. Furthermore, epithelia covering the proboscis and proboscis sac arise from preexisting epidermal epithelium of the larval head, an interpretation that contradicts an earlier suggestion on the origin of these epithelial elements in neogastropods with a feeding larval stage. These results, when compared to foregut development in other gastropods, lead me to propose that the gastropod buccal cavity and buccal mass is a developmental module. Canalized development of this module may have been important to the "evolvability" of the complex gastropod foregut, because it allowed a silent developmental novelty to arise (secondary formation of the postmetamorphic mouth) without disrupting development of the whole module. Nevertheless, this silent novelty might have subsequently facilitated dramatic evolutionary change by allowing the elaborate foregut structure of predatory, postmetamorphic neogastropods to arise in late stage larvae without compromising larval feeding.

Evolutionary diversification of specification mechanisms within the O/P equivalence group of the leech genus Helobdella

Developmental fates and cell lineage patterns are highly conserved in the teloblast lineages that give rise to the segmental ectoderm of clitellate annelids. But previous studies have shown that the pathways involved in specification of the ventrolateral O lineage and the dorsolateral P lineage differ to some degree in distantly related clitellate species such as the leeches Helobdella and Theromyzon, and the sludgeworm Tubifex. To examine this developmental variation at a lower taxonomic level, we have explored the specification pathways of the O and P lineages in the leech genus Helobdella. In leech, the O and P lineages arise from a developmental equivalence group of O/P teloblasts. In this study, we demonstrate that the cell-cell interactions involved in cell fate specification of the O/P equivalence group differ among three laboratory colonies of closely related species. In two populations, the Q lineage is necessary to specify the P fate in the dorsalmost O/P lineage, but in the third population the P fate can be specified by a redundant pathway involving the M lineage. We also observe interspecific variation in the role played by cell interactions within the O/P equivalence group, and in the apparent significance of extrinsic signals from the micromere cell lineages. Our data suggest that cell fate specification in the O/P equivalence group is a complex process that involves multiple cell-cell interactions, and that the developmental architecture of the O/P equivalence group has undergone evolutionary diversification in closely related species, despite maintaining a conserved morphology.

Conservation and evolution of cis-regulatory systems in ascomycete fungi

Relatively little is known about the mechanisms through which gene expression regulation evolves. To investigate this, we systematically explored the conservation of regulatory networks in fungi by examining the cis-regulatory elements that govern the expression of coregulated genes. We first identified groups of coregulated Saccharomyces cerevisiae genes enriched for genes with known upstream or downstream cis-regulatory sequences. Reasoning that many of these gene groups are coregulated in related species as well, we performed similar analyses on orthologs of coregulated S. cerevisiae genes in 13 other ascomycete species. We find that many species-specific gene groups are enriched for the same flanking regulatory sequences as those found in the orthologous gene groups from S. cerevisiae, indicating that those regulatory systems have been conserved in multiple ascomycete species. In addition to these clear cases of regulatory conservation, we find examples of cis-element evolution that suggest multiple modes of regulatory diversification, including alterations in transcription factor-binding specificity, incorporation of new gene targets into an existing regulatory system, and cooption of regulatory systems to control a different set of genes. We investigated one example in greater detail by measuring the in vitro activity of the S. cerevisiae transcription factor Rpn4p and its orthologs from Candida albicans and Neurospora crassa. Our results suggest that the DNA binding specificity of these proteins has coevolved with the sequences found upstream of the Rpn4p target genes and suggest that Rpn4p has a different function in N. crassa.

A systematic examination of the gene regulatory elements in ascomycete fungi reveals striking conservation along with some examples of the ways in which regulatory systems can evolve

Symmetry Breaking and the Evolution of Development

Because of its simplicity, the binary-switch nature of left-right asymmetry permits meaningful comparisons among many different organisms. Phylogenetic analyses of asymmetry variation, inheritance, and molecular mechanisms reveal unexpected insights into how development evolves. First, directional asymmetry, an evolutionary novelty, arose from nonheritable origins almost as often as from mutations, implying that genetic assimilation ("phenotype precedes genotype") is a common mode of evolution. Second, the molecular pathway directing hearts leftward-the nodal cascade-varies considerably among vertebrates (homology of form does not require homology of development) and was possibly co-opted from a preexisting asymmetrical chordate organ system. Finally, declining frequencies of spontaneous asymmetry reversal throughout vertebrate evolution suggest that heart development has become more canalized.

Sexual and Temporal Dynamics of Molecular Evolution in C. elegans Development

Dissection of the phenotypic and molecular details of development and differentiation is a centuries-old topic in evolutionary biology. However, an adequate understanding is missing for the molecular evolution of genes that are expressed differentially throughout development-across time, tissues, and the sexes. In this study, we investigate the dynamics of gene evolution across Caenorhabditis elegans ontogeny and among genes expressed differentially between each sex and gamete type. Using gene classes identified by genome-wide gene expression developmental time series and comparative sequence analysis with the congener C. briggsae, we demonstrate that genes expressed predominantly after reproductive maturity evolve more rapidly than genes expressed earlier in development and that genes expressed transiently during embryogenesis evolve faster than other embryonic transcripts. These results are indicative of relaxed selection on genes expressed after maturity, in accord with the mutation-accumulation model of aging. Furthermore, genes involved in spermatogenesis reveal more rapid evolution than other phenotypic classes of genes. Average rates of evolution among male soma-related genes indicates that selection acts to maintain males in these androdioecious species, despite their rarity, and the rapid evolution of sperm genes suggests that sexual selection acts on sperm development and function.

cis-Regulatory and protein evolution in orthologous and duplicate genes

The relationship between protein and regulatory sequence evolution is a central question in molecular evolution. It is currently not known to what extent changes in gene expression are coupled with the evolution of protein coding sequences, or whether these changes differ among orthologs (species homologs) and paralogs (duplicate genes). Here, we develop a method to measure the extent of functionally relevant cis-regulatory sequence change in homologous genes, and validate it using microarray data and experimentally verified regulatory elements in different eukaryotic species. By comparing the genomes of Caenorhabditis elegans and C. briggsae, we found that protein and regulatory evolution is weakly coupled in orthologs but not paralogs, suggesting that selective pressure on gene expression and protein evolution is quite similar and persists for a significant amount of time following speciation but not gene duplication. Additionally, duplicates of both species exhibit a dramatic acceleration of both regulatory and protein evolution compared to orthologs, suggesting increased directional selection and/or relaxed selection on both gene expression patterns and protein function in duplicate genes.

Common pattern of evolution of gene expression level and protein sequence in Drosophila

Sequence divergence scaled by variation within species has been used to infer the action of selection upon individual genes. Applying this approach to expression, we compared whole-genome whole-body RNA levels in 10 heterozygous Drosophila simulans genotypes and a pooled sample of 10 D. melanogaster lines using Affymetrix Genechip. For 972 genes expressed in D. melanogaster, the transcript level was below detection threshold in D. simulans, which may be explained either by sequence divergence between the primers on the chip and the mRNA transcripts or by down-regulation of these genes. Out of 6,707 genes that were expressed in both species, transcript level was significantly different between species for 534 genes (at P < 0.001). Genes whose expression is under stabilizing selection should exhibit reduced genetic variation within species and reduced divergence between species. Expression of genes under directional selection in D. simulans should be highly divergent from D. melanogaster, while showing low genetic variation in D. simulans. Finally, the genes with large variation within species but modest divergence between species are candidates for balancing selection. Rapidly diverging, low-polymorphism genes included those involved in reproduction (e.g., Mst 3Ba, 98Cb; Acps 26Aa, 63F; and sperm-specific dynein). Genes with high variation in transcript abundance within species included metallothionein and hairless, both hypothesized to be segregating in nature because of gene-by-environment interactions. Further, we compared expression divergence and DNA substitution rate in 195 genes. Synonymous substitution rate and expression divergences were uncorrelated, whereas there was a significant positive correlation between nonsynonymous substitution rate and expression divergence. We hypothesize that as a substantial fraction of nonsynonymous divergence has been shown to be adaptive, much of the observed expression divergence is likewise adaptive.

Evidence that Egfr contributes to cryptic genetic variation for photoreceptor determination in natural populations of Drosophila melanogaster

One objective of quantitative genetics is to identify the nucleotide variants within genes that contribute to phenotypic variation and susceptibility [1]. In an evolutionary context, this means characterizing the molecular polymorphisms that modify the penetrance and expressivity of perturbed traits. A survey of association between 267 SNPs in almost 11 kb of the D. melanogaster Egfr and the degree of eye roughening due to a gain-of-function Egfr(E1) allele crossed into 210 isogenic wild-type lines provides evidence that a handful of synonymous substitutions supply cryptic variation for photoreceptor determination. Ten sites exceed Bonferroni threshold for association in two sets of crosses to different Egfr(E1) backgrounds including a particularly significant cluster of sites in tight linkage disequilibrium toward the 3' end of the coding region. Epistatic interaction of this cluster with one other site enhances the expressivity of this haplotype. Replication of the strongest associations with an independent sample of 302 phenotypically extreme individuals derived from 1000 crosses of Egfr(E1) to freshly trapped males was achieved using modified case-control and transmission-disequilibrium tests. A tendency for the rarer alleles to have more disrupted eye development suggests that mutation-selection balance is a possible mechanism contributing to maintaining cryptic variation for Egfr.

A distinct patterning mechanism of O and P cell fates in the development of the rostral segments of the leech Helobdella robusta: implications for the evolutionary dissociation of developmental pathway and morphological outcome

Despite a high degree of homonomy in the segmental organization of the ectoderm, the body plan of the leech is divided into two zones based on the distinct cell lineage patterns that give rise to the O/P portion of the segmental ectoderm. In the midbody and caudal segments, each segmental repeat of ectoderm arises in part from one o' blast cell and one `p' blast cell. These two blast cells are positionally specified to distinct O and P fates,and give rise to differentiated descendant cells called O and P pattern elements, respectively. In the rostral segments, each segmental repeat of O and P pattern elements arises from a single `op' blast cell. Based on their developmental fates and their responses to the ablation of neighboring cells,the granddaughters of the primary op blast cell are categorized into two O-type cells and two P-type cells. The O-type cells do not require the presence of the rest of the op blast cell clone for their normal development. By contrast, normal development of the P-type cells depends upon interactions with the other OP sublineages. Additional experiments showed that the O-type cells are the source of a repressive signal involved in the normal fate specification of the P-type cells. Our data suggest that the cell interactions involved in fate specification differ substantially in the rostral and midbody segments, even though the set of differentiated descendants produced by the rostral OP pathway and the midbody O and P pathways are very similar.

Evolution of the APETALA3 and PISTILLATA lineages of MADS-box-containing genes in the basal angiosperms

The B class genes, including homologs of the Arabidopsis loci APETALA3 (AP3) and PISTILLATA (PI ), appear to play a conserved role in the determination of petal and stamen identity across core eudicot angiosperms. Understanding how and when these functions evolved is a critical component of elucidating the evolution of flowers, particularly the appearance of petaloid perianth organs. Before comparisons of gene expression patterns or functions can be made, however, it is necessary to establish the orthology of AP3 and PI homologs from basal angiosperms. Here, we report the identification and analysis of 29 new representatives of the B gene lineage from basal ANITA and magnoliid dicot angiosperms. These studies indicate that gene duplications have occurred at every phylogenetic level, both before and after the duplication that produced the separate AP3 and PI lineages. Comparison of genomic structure among PI homologs indicates that a 12-nucleotide deletion that had been considered synapomorphic for the whole PI lineage actually arose within the ANITA grade, after the split of the Nymphaeales but before the separation of the Austrobaileyales. Evidence for alternative splicing of the Nymphaea AP3 homolog is also presented. The implications of these findings for angiosperm systematics, the conservation of AP3 and PI gene function, and the evolution of the ABC program are discussed.

Functional tests of enhancer conservation between distantly related species

Expression patterns of orthologous genes are often conserved, even between distantly related organisms, suggesting that once established, developmental programs can be stably maintained over long periods of evolutionary time. Because many orthologous transcription factors are also functionally conserved, one possible model to account for homologous gene expression patterns, is conservation of specific binding sites within cis-regulatory elements of orthologous genes. If this model is correct, a cis-regulatory element from one organism would be expected to function in a distantly related organism. To test this hypothesis, we fused the green fluorescent protein gene to neuronal and muscular enhancer elements from a variety of Drosophila melanogaster genes, and tested whether these would activate expression in the homologous cell types in Caenorhabditis elegans. Regulatory elements from several genes directed appropriate expression in homologous tissue types, suggesting conservation of regulatory sites. However, enhancers of most Drosophila genes tested were not properly recognized in C. elegans, implying that over this evolutionary distance enough changes occurred in cis-regulatory sequences and/or transcription factors to prevent proper recognition of heterospecific enhancers. Comparisons of enhancer elements of orthologous genes between C. elegans and C. briggsae revealed extensive conservation, as well as specific instances of functional divergence. Our results indicate that functional changes in cis-regulatory sequences accumulate on timescales much shorter than the divergence of arthropods and nematodes, and that mechanisms other than conservation of individual binding sites within enhancer elements are responsible for the conservation of expression patterns of homologous genes between distantly related species.

Regulatory punctuated equilibrium and convergence in the evolution of developmental pathways in direct-developing sea urchins

We made hybrid crosses between closely and distantly related sea urchin species to test two hypotheses about the evolution of gene regulatory systems in the evolution of ontogenetic pathways and larval form. The first hypothesis is that gene regulatory systems governing development evolve in a punctuational manner during periods of rapid morphological evolution but are relatively stable over long periods of slow morphological evolution. We compared hybrids between direct and indirect developers from closely and distantly related families. Hybrids between eggs of the direct developer Heliocidaris erythrogramma and sperm of the 4-million year distant species H. tuberculata, an indirect developer, restored feeding larval structures and paternal gene expression that were lost in the evolution of the direct-developing maternal parent. Hybrids resulting from the cross between eggs of H. erythrogramma and sperm of the 40-million year distant indirect-developer Pseudoboletia maculata are strikingly similar to hybrids between the congeneric hybrids. The marked similarities in ontogenetic trajectory and morphological outcome in crosses of involving either closely or distantly related indirect developing species indicates that their regulatory mechanisms interact with those of H. erythrogramma in the same way, supporting remarkable conservation of molecular control pathways among indirect developers. Second, we tested the hypothesis that convergent developmental pathways in independently evolved direct developers reflect convergence of the underlying regulatory systems. Crosses between two independently evolved direct-developing species from two 70-million year distant families, H. erythrogramma and Holopneustes purpurescens, produced harmoniously developing hybrid larvae that maintained the direct mode of development and did not exhibit any obvious restoration of indirect-developing features. These results are consistent with parallel evolution of direct-developing features in these two lineages.

Descent with modification: the unity underlying homology and homoplasy as seen through an analysis of development and evolution

Homology is at the foundation of comparative studies in biology at all levels from genes to phenotypes. Homology is similarity because of common descent and ancestry, homoplasy is similarity arrived at via independent evolution. However, given that there is but one tree of life, all organisms, and therefore all features of organisms, share some degree of relationship and similarity one to another. That sharing may be similarity or even identity of structure and the sharing of a most recent common ancestor--as in the homology of the arms of humans and apes--or it may reflect some (often small) degree of similarity, such as that between the wings of insects and the wings of birds, groups whose shared ancestor lies deep within the evolutionary history of the Metazoa. It may reflect sharing of entire developmental pathways, partial sharing, or divergent pathways. This review compares features classified as homologous with the classes of features normally grouped as homoplastic, the latter being convergence, parallelism, reversals, rudiments, vestiges, and atavisms. On the one hand, developmental mechanisms may be conserved, even when a complete structure does not form (rudiments, vestiges), or when a structure appears only in some individuals (atavisms). On the other hand, different developmental mechanisms can produce similar (homologous) features. Joint examination of nearness of relationship and degree of shared development reveals a continuum within an expanded category of homology, extending from homology --> reversals --> rudiments --> vestiges --> atavisms --> parallelism, with convergence as the only class of homoplasy, an idea that turns out to be surprisingly old. This realignment provides a glimmer of a way to bridge phylogenetic and developmental approaches to homology and homoplasy, a bridge that should provide a key pillar for evolutionary developmental biology (evo-devo). It will not, and in a practical sense cannot, alter how homoplastic features are identified in phylogenetic analyses. But seeing rudiments, reversals, vestiges, atavisms and parallelism as closer to homology than to homoplasy should guide us toward searching for the common elements underlying the formation of the phenotype (what some have called the deep homology of genetic and/or cellular mechanisms), rather than discussing features in terms of shared or independent evolution.

Rapid coevolution of the nematode sex-determining genes fem-3 and tra-2

Unlike many features of metazoan development, sex determination is not widely conserved among phyla. However, the recent demonstration that one gene family controls sexual development in Drosophila, C. elegans, and vertebrates suggests that sex determination mechanisms may have evolved from a common pathway that has diverged radically since the Cambrian. Sex determination gene sequences often evolve quickly, but it is not known how this relates to higher-order pathways or what selective or neutral forces are driving it. In such a rapidly evolving developmental pathway, the fate of functionally linked genes is of particular interest. To investigate a pair of such genes, we cloned orthologs of the key C. elegans male-promoting gene fem-3 from two sister species, C. briggsae and C. remanei. We employed RNA interference to show that in all three species, the male-promoting function of fem-3 and its epistatic relationship with its female-promoting upstream repressor, tra-2, are conserved. Consistent with this, the FEM-3 protein interacts with TRA-2 in each species, but in a strictly species-specific manner. Because FEM-3 is the most divergent protein yet described in Caenorhabditis and the FEM-3 binding domain of TRA-2 is itself hypervariable, a key protein-protein interaction is rapidly evolving in concert. Extrapolation of this result to larger phylogenetic scales helps explain the dissimilarity of the sex determination systems across phyla.

Speciation despite gene flow when developmental pathways evolve

Evolutionary biologists assume that species formation requires a drastic reduction in gene exchange between populations, but the rate sufficient to prevent speciation is unknown. To study speciation, we use a new class of population genetic models that incorporate simple developmental genetic rules, likely present in all organisms, to construct the phenotype. When we allow replicate populations to evolve in parallel to a new, shared optimal phenotype, often their hybrids acquire poorly regulated phenotypes: Dobzhansky-Muller incompatibilities arise and postzygotic reproductive isolation evolves. Here we show that, although gene exchange does inhibit this process, it is the proportion of migrants exchanged (m) rather than the number of migrants (Nm) that is critical, and rates as high as 16 individuals exchanged per generation still permit the evolution of postzygotic isolation. Stronger directional selection counters the inhibitory effect of gene flow, increasing the speciation probability. We see similar results when populations in a standard two-locus, two-allele Dobzhansky-Muller model are subject to simultaneous directional selection and gene flow. However, in developmental pathway models with more than two loci, gene flow is more able to impede speciation. Genetic incompatibilities arise as frequent by-products of adaptive evolution of traits determined by regulatory pathways, something that does not occur when phenotypes are modeled using the standard, additive genetic framework. Development therefore not only constrains the microevolutionary process, it also facilitates the interactions among genes and gene products that make speciation more likely-even in the face of strong gene flow.

The best of all worlds or the best possible world? Developmental constraint in the evolution of beta-tubulin and the sperm tail axoneme

Through evolutionary history, some features of the phenotype show little variation. Stabilizing selection could produce this result, but the possibility also exists that a feature is conserved because it is developmentally constrained--only one or a few developmental mechanisms can produce that feature. We present experimental data documenting developmental constraint in the assembly of the motile sperm tail axoneme. The 9+2 microtubule architecture of the eukaryotic axoneme has been deeply conserved. We argue that the quality of motility supported by axonemes with this morphology explains their long conservation, rather than a developmental necessity for the 9+2 architecture. However, our functional tests in Drosophila spermatogenesis reveal considerable constraint in the coevolution of testis-specific beta-tubulin and the sperm tail axoneme. The evolution of testis beta-tubulins used in insect sperm tail axonemes is highly punctuated, indicating some pressure acting on their evolution. We provide a mechanistic explanation for their punctuated evolution by testing structure-function relationships between testis beta-tubulin and the motile axoneme in D. melanogaster. We discovered that a highly conserved sequence feature of beta-tubulins used in motile axonemes is needed to specify central pair formation. Second, our data suggest that cooperativity in the function of internal beta-tubulin amino acids is needed to support the long axonemes characteristic of Drosophila sperm tails. Thus, central pair formation constrains the evolution of the axoneme motif, and intramolecular cooperativity makes the evolution of the internal residues path dependent, which slows their evolution. Our results explain why a highly specialized beta-tubulin is needed to construct the Drosophila sperm tail axoneme. We conclude that these constraints have fixed testis-specific beta-tubulin identity in Drosophila.

Developmental evolution: getting robust about robustness

In the context of development, a process is robust if it can proceed normally despite the enormous capacity for perturbation inherent in all biological systems. A new mode of theoretical modeling of genetic networks holds great promise for increasing our understanding of both the quantitative mechanisms of robustness and its evolutionary impact.

Molecular evolution of the HoxA cluster in the three major gnathostome lineages

The duplication of Hox clusters and their maintenance in a lineage has a prominent but little understood role in chordate evolution. Here we examined how Hox cluster duplication may influence changes in cluster architecture and patterns of noncoding sequence evolution. We sequenced the entire duplicated HoxAa and HoxAb clusters of zebrafish (Danio rerio) and extended the 5' (posterior) part of the HoxM (HoxA-like) cluster of horn shark (Heterodontus francisci) containing the hoxa11 and hoxa13 orthologs as well as intergenic and flanking noncoding sequences. The duplicated HoxA clusters in zebrafish each house considerably fewer genes and are dramatically shorter than the single HoxA clusters of human and horn shark. We compared the intergenic sequences of the HoxA clusters of human, horn shark, zebrafish (Aa, Ab), and striped bass and found extensive conservation of noncoding sequence motifs, i.e., phylogenetic footprints, between the human and horn shark, representing two of the three gnathostome lineages. These are putative cis-regulatory elements that may play a role in the regulation of the ancestral HoxA cluster. In contrast, homologous regions of the duplicated HoxAa and HoxAb clusters of zebrafish and the HoxA cluster of striped bass revealed a striking loss of conservation of these putative cis-regulatory sequences in the 3' (anterior) segment of the cluster, where zebrafish only retains single representatives of group 1, 3, 4, and 5 (HoxAa) and group 2 (HoxAb) genes and in the 5' part of the clusters, where zebrafish retains two copies of the group 13, 11, and 9 genes, i.e., AbdB-like genes. In analyzing patterns of cis-sequence evolution in the 5' part of the clusters, we explicitly looked for evidence of complementary loss of conserved noncoding sequences, as predicted by the duplication-degeneration-complementation model in which genetic redundancy after gene duplication is resolved because of the fixation of complementary degenerative mutations. Our data did not yield evidence supporting this prediction. We conclude that changes in the pattern of cis-sequence conservation after Hox cluster duplication are more consistent with being the outcome of adaptive modification rather than passive mechanisms that erode redundancy created by the duplication event. These results support the view that genome duplications may provide a mechanism whereby master control genes undergo radical modifications conducive to major alterations in body plan. Such genomic revolutions may contribute significantly to the evolutionary process.

Microarrays in ecology and evolution: a preview

Microarray technology provides a new tool with which molecular ecologists and evolutionary biologists can survey genome-wide patterns of gene expression within and among species. New analytical approaches based on analysis of variance will allow quantification of the contributions of among individual variation, genotype, sex, microenvironment, population structure, and geography to variation in gene expression. Applications of this methodology are reviewed in relation to studies of mechanisms of adaptation and divergence; delineation of developmental and physiological pathways and networks; characterization of quantitative genetic parameters at the level of transcription ('quantitative genomics'); molecular dissection of parasitism and symbiosis; and studies of the diversification of gene content. Establishment of microarray resources is neither prohibitively expensive nor technologically demanding, and a commitment to development of gene expression profiling methods for nonmodel organisms could have a tremendous impact on molecular and genetic research at the interface of organismal and population biology.